计算化学公社

标题: 求助：循环每一帧轨迹文件进行小模拟，结果press为-nan [打印本页]

作者
Author: pumpkin111 时间: 2024-9-23 09:50
标题: 求助：循环每一帧轨迹文件进行小模拟，结果press为-nan
本帖最后由 pumpkin111 于 2024-9-25 10:15 编辑

我用lammps做热解模拟，用的势文件是laspnn，体系只有CHO三种元素。因为不能用fix reaxff/species命令但是想得到species文件，所以想循环每一帧轨迹文件用reaxff做一个非常小小的模拟过程，来输出产物文件。然后我用轨迹文件中的某一帧做模拟，运行报错如下。但是用MS构建的初始模型没有报错并且输出了想要的Species文件。

WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 0 step (../fix_qeq_reax.cpp:735)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 0 step (../fix_qeq_reax.cpp:735)
Per MPI rank memory allocation (min/avg/max) = 73.27 | 77.51 | 82.9 Mbytes
Step Temp E_pair E_mol TotEng Press
0 1500 110428.05 0 115252.49 -nan

电荷平衡不收敛且压力显示不出来，应该是模型的几何结构不合理。请问各位老师，是应该修改in文件吗？还是说输出的轨迹不能作为模型文件呢？
in文件如下

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