求助：mdp文件中耦合时间常数、约束键等参数的设置问题

查看数: 2055 | 评论数: 6 | 收藏 Add to favorites 4

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

发布时间: 2023-5-9 15:54

本帖最后由 wmy1217 于 2023-5-9 15:57 编辑最近在做关于金属离子与蛋白质的分子动力学模拟，金属离子用的是虚拟原子模型，参考了这篇文章的虚拟原子模型参数（https://pubs.acs.org/doi/10.1021/acs.jctc.6b002 ...

回复 Reply

zdworld 发表于 Post on 2023-5-23 12:09:10

wmy1217 发表于 2023-5-21 14:30
非常感谢，我想请问如果不对这个键进行约束，对于整个模拟会有什么影响吗

这个不一定，我自己的体系不约束就崩了，你可以跑个短时间的试试

wmy1217 发表于 Post on 2023-5-21 14:30:46

zdworld 发表于 2023-5-18 00:45
这个提示就是这个键的振动周期太短，步长一般要求低于10%最小振动周期，步长改小就行，说实话步长改小计算 ...

非常感谢，我想请问如果不对这个键进行约束，对于整个模拟会有什么影响吗

zdworld 发表于 Post on 2023-5-18 00:45:40

这个提示就是这个键的振动周期太短，步长一般要求低于10%最小振动周期，步长改小就行，说实话步长改小计算量要翻倍，不知道能不能单独约束这个键的长度。

sobereva 发表于 Post on 2023-5-10 22:06:22

wmy1217 发表于 2023-5-10 12:57
非常感谢老师的解答！我还想请教一下，使用shake算法与lincs相比有何不同呢？

wmy1217 发表于 Post on 2023-5-10 12:57:12

sobereva 发表于 2023-5-9 22:18
用线性升温
甭用朗之万动力学，非主流，毫无用的必要
始终用constrains = hbonds，有问题该怎么解决怎么 ...

非常感谢老师的解答！我还想请教一下，使用shake算法与lincs相比有何不同呢？

sobereva 发表于 Post on 2023-5-9 22:18:56

用线性升温
甭用朗之万动力学，非主流，毫无用的必要
始终用constrains = hbonds，有问题该怎么解决怎么解决，别试图用all-bonds逃避
若无特殊情况，v-rescale热浴用tau_t = 0.2，PR压浴用tau_p = 2.0

别盲目效仿文献

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