J-OCTA在材料研发领域的应用

J-OCTA软件作为分子动力学计算软件，通过对材料从原子级别到微米级别的模拟计算实现对几乎所有材料的分析。本文主要讲述该软件在材料研发领域的表现，着重介绍J-OCTA通过建模、求解、分析等功能在材料多尺度仿真方面的重要作用。

（图1：离子液体的MD模拟）

（图2：全原子及粗粒化模型示意图）

J-OCTA作为计算软件优于其他软件的地方在于J-OCTA中所有输入数据和计算结果都对用户开放，便于定制运算。其他商业软件（如Materials Studio,简称MS）的数据通常以二进制格式存取，不会公开，只可作为黑箱使用。

J-OCTA通过开源多尺度仿真平台用UDF文本传递数据、利用并行求解器+功能扩展及外部扩展器解析，将所有所有功能集成到GUI界面中，为用户使用提供了较大便利。

（图3：软件架构）

除此外，J-OCTA在材料研发领域还有很多区别于其他同类软件的优势。

1. 速度快：VSOP和LAMMPS在同类软件中的计算速度较快，用户可以使用VSOP和LAMMPS进行大规模的计算(原子数超过100K+)，计算效率是MS软件的N倍（实测N＞5）。

（图4：由1200万个粒子组成的大规模复合材料模型）

2. 精度高：J-OCTA的内置工具可以通过量子化学的计算结果对力场参数进行优化。此外，支持EAM（金属）、PCFF（聚合物）、混合（EAM+LJ）等多种力场。

（图5：分子二面角势能优化）

3. 反应快：根据量子化学计算结果，可以在MD(VSOP)中进行反应模拟。J-OCTA内置GUI可以简单便捷地设置反应位点。更重要的是，计算效率是MS软件的几十倍（1天vs2个月），物超所值！

（图6：环氧树脂的固化反应）

4. 界面：J-OCTA利用SIESTA(基于原子轨道基底的DFT计算软件)作为内置的界面能计算工具，可以预测分子和固体表面之间的界面能，以此来确定LJ势。

（图7：界面能和势能优化的DFT计算）

5. 建模：J-OCTA可以通过将珠簧模型(Kremer-Grest)反向映射为全原子模型来建立无定形结构。该功能可以显著提高聚合物的构建效率。

（图8：利用反向映射建立无定形结构）

6. 流变：具有独特的“流管模型”，可用来模拟流变特性。只能在J-OCTA中通过PASTA、NAPLES和VSOP-DPD实现。利用VSOP-DPD，可以用相分离和填料分散来评价熔融聚合物的动力学特性。

（图9：利用原语链表网络模型预测熔融聚合物粘弹性）

7. 介观：对于介观尺度上的模拟，尤其是SCF或DPD，J-OCTA显著优于其他软件，其求解器功能强大且效率很高。

8. 粗粒化：J-OCTA支持粗粒化势能的自动推算功能。同时支持反向映射函数功能，可以从SCF缩放成MD或FEM模型。其中从SCF到MD的反向映射很少见，该功能可以简单便捷地构建界面结构。

（图10：利用SCF得到的相分离结构导入为FEM模型）

9. 有限元：对于有限元建模，J-OCTA可以使用相分离结构数据创建STL数据、LS-DYNA和Nastran的体素网格。LS-DYNA可以直接在J-OCTA中调用。

（图11：相分离结构、体素网格以及使用LS-DYNA模拟单轴拉伸）

10. RVE : J-OCTA可以用于构建复合材料的复杂结构，该结构称为代表性体积单元(RVE)，这一功能支持建立高填充填料的分散结构模型。

（图12：代表性体积单元(RVE)示例）

11. 费用低：除了VSOP，J-OCTA中大部分求解器都是开源的，建立模型后用户可以在服务器上免费进行运算。

12. 机器学习：对于数据科学，J-OCTA具有基于机器学习(ML)的定量构效-属性关系(QSPR)预测功能。该功能可以通过学习化学结构和材料属性数据之间的关系来预测物理特性（如图13）。

（图13：使用ML-QSPR法预测聚合物密度 ）

13. 批处理：为了分析庞大的数据，利用MD建模接口和Python脚本，支持高通量的分子动力学模拟，可以在命令行中构建分子模型并连续执行MD模拟命令。同时，利用Python脚本可以进行多种后处理分析。

14. 后处理：可以将其他软件（LAMMPS等）的计算结果转换为J-OCTA格式，以便灵活地输出想要的结果。

15. 便于学习的帮助文档：包含案例手册及用户手册。

（图14：帮助文档界面）

多达10多种的软件优势，在实际应用中是如何工作的，以下将做部分应用举例。

1. 高分子的交联结构模拟

（1）利用J-OCTA/ VSOP，模拟DGEBA、44DDS交联为环氧树脂的过程（如图15）

（图15，Ref: T. Okabe,“Atomistic simulation of curing and mechanical properties of Epoxy resin”, J-OCTA Users Conference 2014）

（2）弹性模量的模拟（如图16）

（图16，Ref: V. Sundararaghavan, A. Kumar, International Journal of Plasticity, 47, 111 (2013)）

• NPT系综、伸长率为10.22m/s

• 施加5%应变，拟合应力-应变曲线得到弹性模量

• 模拟得到的弹性模量:4.22GPa(实验值:3.76GPa)

2. 机器学习（如图17）

（图17）

3. 代表性体积单元（RVE）µ-FEM（如图18）

• 从分散相到STL界面或体素网格

• 从分散填料到STL或位置信息

（图18）

4. 分散相结构

• 非线性结构分析：采用J-OCTA获得橡胶、树脂的相分离（如图19、20）

（图19）

（图20）

5. 粗粒化（CG）模型的典型案例（如图21）

• 具有交联结构的橡胶

• 交联、填料-聚合物作用的影响

（图21：左图为无交联网络、不同变形速度下的应力应变曲线；右图为存在交联网络（变形速度为0.001）的应力应变曲线）

6. 循环拉伸变形下的填充橡胶模拟

（1）从形态学（Morphology）的角度讨论动态响应差异的来源（如图22）

（图22，数据来源：）

（2）滞后机理（如图23）

（图23，数据来源：）

7. 缠结熔体聚合物的动力学行为（如图24）

（图24）

如果您对分子动力学或J-OCTA软件感兴趣，可以关注我们的官方公众号（庭田科技）了解更多内容。同时我们诚挚邀请您参与我们于10月24日及10月26日开展的针对J-OCTA软件专题研讨会，详情请关注官网及官方公众号（庭田科技），期待您的到来！

相关产品链接：http://www.anscos.com/jocta.html

J-OCTA专题研讨会：邀请函- 庭田科技 (anscos.com)

如需更多技术咨询，请随时与我们联系：

全国热线：400 633 6258

官方邮箱：info@anscos.com