计算化学公社

标题: 已解决：lammps reaxff力场模拟界面模型高温热分解时遇到hbondchk failed的问题 [打印本页]

作者
Author: 李天 时间: 2023-8-13 21:41
标题: 已解决：lammps reaxff力场模拟界面模型高温热分解时遇到hbondchk failed的问题
本帖最后由李天于 2023-8-16 19:19 编辑

各位老师，我在使用lammps reaxff力场模拟界面模型高温热分解时遇到hbondchk failed的问题。
我的模型是由炸药分子和聚合物粘结剂组成，具体建模流程为：
①（ms中）炸药分子单胞compass优化-->切面-->扩胞--->forcite compass优化
②（ms中）画出聚合物粘结剂单链-->forcite compass优化-->AC模块compass力场混合多条粘结剂链
③利用ms中的build layer模块组合上下两层界面模型

此时ms中的模型密度达不到混合物的计算理论密度，且通过调整上下两层的真空层距离和晶胞大小也达不到理论密度。于是我ms建立完模型后，输出data，利用lammps npt系综，在界面上下两层z轴方向施加压力，压缩盒子到模型理论密度。压缩后弛豫先跑了nvt 300k 10ps和npt 300k 0bar 10ps，模型已经达到平衡，且模型结构和密度未发现问题，力场适合。
接着跑nvt高温热分解，2500k 2750k 3000k 等温度。发现2500k并行20核可以成功跑完300ps，但是只是换了一个温度，例如3000k，就会报hbondchk failed的错误。

查找网络给出回答是可以在in文件中添加safezone mincap命令，但是回答说这会忽略hbondchk failed这个错误，治标不治本。
所以具体问题出在哪里，想不通，麻烦老师们帮助解答，谢谢大家！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)