MD用GULP热解时，设置nvt系综但是跑完温度却变化

zaozaojue

各位老师好，最近在用MS的GULP模块进行分子动力学模拟热解仿真，但是在进行热解的时候，明明选取的是nvt系综，但是跑出来的温度趋势却在变化，请问各位老师知道是什么原因吗
初始模型：
在600k下的gulp运算：
整个运算温度变化：
热解后模型：
在文献中600k不应该全都热解完全，但是我的模型全都热解成小分子了