补充：我之前以为是步长的原因。我的体系都是轻原子，仍然将步长设为了0.5 fs。上图就是步长为0.5 fs时的温度图。

timecon设为10，3000K的NVT上下波动1千多K

温度波动跟体系大小关系比较大，AIMD的时候体系一般很小，温度波动大很正常

我跑2000K的温度也这样，上下浮动200左右，应该正常吧。timecon设的小点温度耦合越强，可能波动就小一点。但是不知道设置多大合适。希望有大佬解答。