24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
查看: 663  |  回复: 0
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 伊里斯K 的 5 个金币

伊里斯K

新虫 (初入文坛)

[交流] 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客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 伊里斯K 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见