计算化学公社

标题: 用gromacs模拟己酸甲脂150摄氏度界面时nvt出现真空区怎么解决？ [打印本页]

作者
Author: shanxitaiyuan 时间: 2018-6-3 18:15
标题: 用gromacs模拟己酸甲脂150摄氏度界面时nvt出现真空区怎么解决？
模拟过程1.首先建立一个5*5*5的盒子
2.能量最小化100000步
3.NPT步长0.002ps，跑20000000步，计算长程力采用cut-off方法，采用Berendsen压力耦合。温度采用V-rescale控温跑后盒子 4.87489*4.87489*4.87489
4.增大盒子x,y方向不变，z方向增大到10
5.NVT模拟，步长0.002ps，15000000步，计算长程采用cut-off方法，也采用Berendsen压力耦合
出现了真空区，不知道什么原因，是因为盒子增大的不够吗？

作者
Author: sobereva 时间: 2018-6-4 04:43
增大盒子跑NVT不就是为了能有一部分真空区

作者
Author: wbn 时间: 2018-6-4 05:34
NVT哪来的压力耦合

作者
Author: shanxitaiyuan 时间: 2018-6-4 11:53

sobereva 发表于 2018-6-4 04:43
增大盒子跑NVT不就是为了能有一部分真空区

是盒子中间有真空区

作者
Author: shanxitaiyuan 时间: 2018-6-4 11:54

wbn 发表于 2018-6-4 05:34
NVT哪来的压力耦合

是的没有，写错了

作者
Author: sobereva 时间: 2018-6-4 13:57

shanxitaiyuan 发表于 2018-6-4 11:53
是盒子中间有真空区

没图说不清楚
本身有PBC条件，中间有真空区和盒子上方有真空区本身是一码事

作者
Author: wbn 时间: 2018-6-5 00:40

shanxitaiyuan 发表于 2018-6-4 11:53
是盒子中间有真空区

一般盒子没平衡好的话有可能中间个把原子离得过近或者初始动能过大，导致放开体积后slab弹开。不清楚为什么你NPT了2ns还会出现这种情况。建议试一下在z方向伸长之前做两ns的NVT。记得mdp文件中不要重新生成velocity.T-couple建议改用Nose-Hoover

作者
Author: k64_cc 时间: 2018-6-7 18:01

wbn 发表于 2018-6-5 00:40
一般盒子没平衡好的话有可能中间个把原子离得过近或者初始动能过大，导致放开体积后slab弹开。不清楚为什 ...

CSVR怎么看都比Nose-Hoover强好吧……

作者
Author: 一条君 时间: 2020-12-10 11:58
请问最后怎么解决真空问题呢？另外可以求一些用Berendsen压力耦合的文献吗？可以用它来发文章吗，谢谢

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