Wenchuan 发表于 2023-10-30 11:13
后面的模拟又发现系统温度降低是盒子z方向太长的缘故（800 nm），看来盒子某一方向尺寸太大会出现bug。

这个bug是在NVT情况下出现，NPT下情况可能好点。

sobereva 发表于 2023-10-26 07:54
控温不是立刻把温度拉到期望的温度，有一定速度，取决于tau_t。跑足够长时间，最终早晚能维持到期望的温度

后面的模拟又发现系统温度降低是盒子z方向太长的缘故（800 nm），看来盒子某一方向尺寸太大会出现bug。

sobereva 发表于 2023-10-26 07:54
控温不是立刻把温度拉到期望的温度，有一定速度，取决于tau_t。跑足够长时间，最终早晚能维持到期望的温度

感谢解答，温度降低可能是因为气体膨胀时粒子在z方向的速度很大，计算时控温算法将此速度加入到温度计算中，所以认为温度很高并降温。但是模拟很长时间粒子温度依旧很低，气体发生团聚，温度似乎无法回升。可能解决方案是要么不要模拟气体膨胀过程，要么不对膨胀方向控温。

控温不是立刻把温度拉到期望的温度，有一定速度，取决于tau_t。跑足够长时间，最终早晚能维持到期望的温度