已解决：lammps reaxff力场模拟界面模型高温热分解时遇到hbondchk failed的问题

李天

各位老师，我在使用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这个错误，治标不治本。
所以具体问题出在哪里，想不通，麻烦老师们帮助解答，谢谢大家！

Stephen001

你好，我也是遇到了同样的问题，尝试调大收敛容差和迭代次数，也尝试过调整处理器的核数，但都不行，请问您解决了这个问题，能够传授下经验？

YQpolymer

老哥，我也遇到了这个问题，可以交流下么

NaNaNai

帖子可以开放查看么 我也遇到了这个问题

ddddgc

您好请问是如何解决的，能否共享一下，目前此贴仅作者可见，感谢

xuhuikai

我也遇到了同样的问题，想在这里学习一下解决办法都不行，还仅作者可见。。。。差评

蛾眉皓齿

您好，我也遇到了hbondchk failed问题，麻烦问一下您该怎么解决呢？

1651561651

大佬你好，请问这个问题怎么解决的呀？我也是这种错误ERROR on proc 8: step 4: hbondchk failed: H=2 end(H)=80 str(H+1)=75
(src/REAXFF/reaxff_forces.cpp:119)

蛾眉皓齿

求问，作者是如何解决hbondchk failed的呀

蛾眉皓齿

看到帖子回复是仅作者可见，所以求问作者，您是如何解决hbondchk failed这个问题的呀？

jwy20

我也模拟中也遇到这个问题了，因为帖子仅作者可见看不见，想请教下是怎么解决的呢？

Makarov

您好，我也在做reaxff热解相关，但ms模型只有烷烃分子→forcite优化→AC混合的过程，在300k到4000k升温模拟中遇到了与您帖子里相同的hbondchk failed问题，请问您最后是如何解决的呢

LaplaceL

您好，请问最后是怎样解决的呢

15210268712

楼主你好，我也遇见了这个问题，请问解决了吗？

lltll

可以分享一下怎么解决的吗？