计算水合氢离子内聚能的问题

jiangua

各位老师好，用material studio的Forcite模块计算H3O+内聚能时出现问题，具体如下
使用力场：COMPASS，自动分配力场，力场类型与电荷分别为，三个H原子：h1+，0.403，氧原子：O3, -0.209。分子总电荷为+1
计算方法：100个H3O+放入周期性盒子中，结构优化，然后500 ps 的NPT 和500 ps的NVT平衡，然后计算最后10帧的内聚能
结果：H3O+的能聚能为5.93e+010，相比较H2O的 1.94e+090（相同的方法计算），能聚能高非常多。而且由于内聚能高，H3O+的密度达到了1.90 g/cm3，相比水分子的高了一倍。
问题：对于结果难以理解，水合氢离子由于吸收质子而呈现正电性，H3O+与H3O+的相互作用应该低于H2O，也就是相互作用能应该时低于水分子的，但是结果表现为H3O+比H2O更容易团聚。
查了很多论文和资料都没有头绪，猜测是否是力场分配不合理

请有经验的老师给予指导，十分感谢


左图时平衡后水盒子的构象，密度0.95 g/cm3
右图是H3O+的构象，密度为1.90 g/cm3，明显可以看出H3O+相互团聚

sobereva

“100个H3O+放入周期性盒子”这种模拟有无意义我持怀疑态度，放在现实中，这样的盒子里带高密度正电荷的体系立马就炸了。另外，这样的氢键属于带电荷氢键，非常特殊（H与两边的O都形成一定共价作用），靠COMPASS这样的力场的模拟结果我也不认为有什么可信性。