模拟一个MOF+CO2体系，V-rescale热浴出现问题

Wenchuan · 发表于 Post on 2023-10-25 22:32:02

模拟一个MOF+CO2体系

，
模拟开始气体膨胀，将温度设为恒定300K

md.mdp (707 Bytes, 下载次数 Times of downloads: 10) ，为什么气体在膨胀过程中温度好像还会降低呢？温度看起来降低是因为气体膨胀后出现了气体聚集液化

。
是因为V-rescale热浴的问题还是因为气体膨胀过程无法通过设定温度恒定300K来保持恒温，还是盒子太长出现bug了呢？

Wenchuan · 发表于 Post on 2023-10-31 20:51:24

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

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

Wenchuan · 发表于 Post on 2023-10-30 11:13:56

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

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

Wenchuan · 发表于 Post on 2023-10-27 11:53:36

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

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

sobereva · 发表于 Post on 2023-10-26 07:54:23

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[GROMACS] 模拟一个MOF+CO2体系，V-rescale热浴出现问题