lammps从NVT或者NPT切换到NVE时温度持续上升

最近，尝试计算twin-graphene的热导率，为了简化，直接先使用石墨烯的tersoff势能进行EMD下热导率的计算，测试发现NVT和NPT下势能表现正常，但是从系统NVT或NPT切换到NVE开始热导率计算时，观测到温度发散，最终材料解体。

针对这个问题，修改了很多计算参数。最终发现原因如下：

猜测可能来源于势能本身的特征，应用于twin graphene这种内部有局域双层的结构时，可能存在部分原子距离急剧变小的情况，势能没有更好的描述，导致原子在分动模拟中，如果步长太大，会有反常的持续加速的结果（类似于原子靠的很近的势能表现），表现就是温度持续上升，除非有温度控制的系综，解决方法有以下：

1.就是使用很小的步长来准确描述势能急剧变化情况下原子的运动，经测试0.1fs步长，温度能够稳定。但是EMD用这个步长计算太耗时间。

2.当然可以选择专门针对相应体系的其他更好的势能模型，不过没有合适势能的情况下，经测试发现可以加入LJ势能，来避免原子太近或势能不合理的地方。但是tersoff和lj无法直接hybrid，一对原子要么采用tersoff，要么采用LJ描述，只能使用一种，因此，需要对晶胞内的原子类型做区分，并修改原来的tersoff势能，将局域双层结构中的上下层原子区分，关掉相应的tersoff部分，用LJ替代，经过测试发现修正后的tersoff势能能够很好的描述twin-graphene这种结构。

3.使用fix NVE/limit也能够回复正常，这个是系综直接控制原子每步偏离不能超过某个临界值，只是对热输运计算的影响需要进一步思考。