计算化学公社

标题: 模拟一个MOF+CO2体系，V-rescale热浴出现问题 [打印本页]

作者
Author: Wenchuan 时间: 2023-10-25 22:32
标题: 模拟一个MOF+CO2体系，V-rescale热浴出现问题
模拟一个MOF+CO2体系

(, 下载次数 Times of downloads: 11) ，
模拟开始气体膨胀，将温度设为恒定300K (, 下载次数 Times of downloads: 10) ，为什么气体在膨胀过程中温度好像还会降低呢？温度看起来降低是因为气体膨胀后出现了气体聚集液化
(, 下载次数 Times of downloads: 10) 。
是因为V-rescale热浴的问题还是因为气体膨胀过程无法通过设定温度恒定300K来保持恒温，还是盒子太长出现bug了呢？

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

作者
Author: Wenchuan 时间: 2023-10-27 11:53

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

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

作者
Author: Wenchuan 时间: 2023-10-30 11:13

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

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

作者
Author: Wenchuan 时间: 2023-10-31 20:51

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)