求助：溶液体系中给水分子赋初始速度后mdrun报错

查看数: 2277 | 评论数: 4 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

发布时间: 2022-12-27 22:20

老师好，有一个问题困扰了许久，多种尝试均失败，请老师指点：问题：体系是蛋白质分子周围分布球形溶剂分子盒子（如图），在gro文件中给溶剂分子（四点水）规定赋了指向球心方向的初始速度0.001nm/ps（速度格式正 ...

回复 Reply

mol 发表于 Post on 2022-12-29 22:31:18

请教lz，指向球心方向的速度矢量如何计算出来弄到gro里呢

sobereva 发表于 Post on 2022-12-29 03:07:33

韩式蛋炒范发表于 2022-12-28 21:22
感谢老师回复，我是想模拟水分子以指定速度方向对蛋白进行撞击，来寻找什么条件下可以发生蛋白分子结构的 ...

这种事情得反复去试，比如初速度设为0看能不能跑；如果能跑，再加大一丁点再试，反复尝试找原因
也包括尝试简化体系，比如只保留水层而不考虑蛋白去试
用了0.5 fs步长就没必要刻意换成刚性水（换了当然也没问题），此时不会因为水是刚性还是柔性而导致问题

sobereva 发表于 Post on 2022-12-28 03:35:06

用柔性水、不用constraints = hbonds还用2fs步长哪行
尝试1fs，还不行就0.5 fs

体系那么大真空区，根本不能用EnerPres

不知道你当前想模拟什么，如果模拟特定动能的水对蛋白质的撞击，根本不应当用热浴

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