大家好！
小弟目前在做沸腾的分子动力学模拟，遇到了一个问题，想了很久也没有想明白。
我的模型下方黄色为固体原子，底面尺寸为36.15A的正方形，上方为液体水分子。使用“fix heat heat 1 0.0082”命令使热源持续向系统内输入热流，模拟过程中单位选择为metal，故热流为0.0082eV/ps，换算后的热流密度为10^8W/m2。而后在统计液体能量的变化过程中，计算得到dE/dt的值为0.13eV/ps左右，是输入热流的约16倍，在整个模拟过程中，只有热源热流的输入，无其他能量输入，所以这个结果是不符合能量守恒定律的。在单位换算过程中没有错误，所以我想不明白是什么原因导致的这个现象，还请各位帮忙看一下是什么原因。
模拟过程单位步长为0.001ps，模拟先在NVT下进行10W步弛豫，后为250W时间步的NVE，只设置了下方的热源，未设置冷源。
欢迎大家前来讨论！
模型图
液体能量变化

不加热的nve有试过么

你好，不知道你这个问题解决了吗？最近也遇到了fix NVE的时候能量增加的情况

你好，我和你现在出现了差不多的问题，能帮我解决一下吗

楼主你好。我的模型是银原子向水传热，很奇怪，fix heat 或者fix ehex给银添加热流 银原子固定不动 ，但是显示温度升高 想了两周 不知道问题出在哪里

楼主你好，NVE能量不守恒问题解决了吗？我也遇到了相同问题，一直搞不定。可以指点一下吗？