| 查看: 816 | 回复: 0 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 伊里斯K 的 5 个金币 | |||
[交流]
LAMMPS求助：reaxff/lg力场如何模拟金属配合物
|
|||
|
各位老师好，我想要将手上的一批金属有机配合物进行reaxff/lg模拟，但没有找到使用包含铁氮的lg力场相关的文献，于是选择将包含铁氮的非lg力场，按照文献(10.1021/jp201599t)的方法手动加入铁的两个lg参数。Fe的Re依照文献选用92年的UFF中的数据，Clg通过DFT计算出长程能量校正后，依照文献中给出的公式2.3算出。 然而我在实际进行模拟时发现，无论是作为试验用的Fe24N10晶体，还是我想要模拟的铁氮配合物，使用lg力场时都会报错：non-numeric pressure, simulation unstable。我起初认为这是我的初始结构不合理，但即使经过cg/minimize后也依然会出现相同的错误。我尝试着手动更改了这两个数值，貌似无论取多大的值都是一样的报错。但同样的模型，使用非lg力场时却可以正常地模拟，这是否证明我的模型本身没有问题呢。 请问问题是出现在我手动添加的两个数值上吗？还是我的in文件中的指令有错误？ 我的测试文件如下链接：https://pan.baidu.com/s/13l-CTCXRvIjU-u0KmweqMQ?pwd=17zy 提取码：17zy 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有6人回复
