求助：reaxff/lg力场模拟金属配合物

Kmetsch · 发表于 Post on 2023-2-3 19:35:04

各位老师好，我想要将手上的一批金属有机配合物进行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文件中的指令有错误？

Kmetsch · 发表于 Post on 2023-2-5 12:44:01

我又尝试了1-5000范围的Clg值，依然无法运行

Kmetsch · 发表于 Post on 2023-2-7 20:21:21

求助老师们解答

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Lammps] 求助：reaxff/lg力场模拟金属配合物

浏览过的版块