多晶气体水合物模拟能量异常

luogaoyang123

各位晚上好！

我在做气体水合物分子动力学模拟时遇到了一些问题，想请教大家。

模拟体系中，水分子采用TIP4P/ice力场，甲烷用OPLS-AA力场，二氧化碳用TraPPE力场，积分步长1fs。模拟刚开始能量最小化阶段就出现异常，具体情况如下：

• 能量异常
  在能量最小化第30步，系统出现严重能量发散，输出了极端负值（例如 -1.6e+16），并伴随警告信息：
  WARNING: Inconsistent image flags
  我不太理解这里的极端负能量，以前遇到的大多是因原子距离太近而出现正的NaN，为什么会是负的NaN？这种极端负值该如何理解？
  
  

附上log文件，方便大家分析。谢谢！