谈谈怎么判断分子动力学模拟是否达到了平衡

查看数: 47497 | 评论数: 30 | 收藏 Add to favorites 96

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

发布时间: 2021-12-25 23:53

谈谈怎么判断分子动力学模拟是否达到了平衡On the judgement whether molecular dynamics simulation has reached equilibrium 文/Sobereva@北京科音 2021-Dec-24 0 前言时常有人在思想家公社Q ...

回复 Reply

sobereva 发表于 Post on 2025-11-29 05:31:58

cheng1 发表于 2025-11-28 22:20
请教一下老师，比如1000个分子团簇的单体放在盒子里做动力学平衡，具体如何判定是否平衡呢？通过统计簇的数 ...

没法用RMSD判断
结合体系具体特征、团簇的动力学行为考虑什么能拿来说事。如果簇的数目不固定、不是始终就只有一个大的簇出现，可以统计簇的数目。认真看当前帖子，认真领会判断平衡的思路、思想

cheng1 发表于 Post on 2025-11-28 22:20:52

请教一下老师，比如1000个分子团簇的单体放在盒子里做动力学平衡，具体如何判定是否平衡呢？通过统计簇的数目可以吗？针对该体系可不可以理解为RMSD平稳代表体系平衡，而RMSD不平稳不代表体系不平衡？

raku 发表于 Post on 2024-11-12 11:10:14

牛啊来学习了，刚接触MD时间不长，发现论坛真是好工具

sobereva 发表于 Post on 2024-6-16 05:32:44

qikongwanli 发表于 2024-6-16 01:32
社长老师好，我是新手，能帮忙看下我的结果吗，这种还要继续延长模拟时间吗，谢谢。

我模拟的是蛋白质和 ...

既然肉眼都看出来逐渐远离了，当然没平衡

qikongwanli 发表于 Post on 2024-6-16 01:32:52

本帖最后由 qikongwanli 于 2024-6-16 01:36 编辑

社长老师好，我是新手，能帮忙看下我的结果吗，这种还要继续延长模拟时间吗，谢谢。

我模拟的是蛋白质和小分子配体化合物的结合，我看轨迹图，模拟的后期，配体的位置逐渐原理初始的结合位点。

屏幕截图 2024-06-16 010124.png (16.61 KB, 下载次数 Times of downloads: 104)

ldx022 发表于 Post on 2023-8-27 16:34:27

请问社长如果用了增强采样的话，怎么来判断是否平衡了呢？

sobereva 发表于 Post on 2022-10-19 23:57:10

luzujia 发表于 2022-10-19 09:45
好的，谢谢社长。借楼再问一个问题，就是xtb在做动力学的时候，设定了初始体系的单电子数，后来体系发生 ...

xtb程序里的GFN-xTB默认做的是限制性的有限温度计算，你说的情况未必不能描述。比如两个离得很远的甲基自由基，GFN-xTB的这种计算直接能得到两个占据数几乎为1的轨道，都不用设对称破缺初猜。

luzujia 发表于 Post on 2022-10-19 09:45:22

sobereva 发表于 2022-10-19 01:09
xtb本来就不适合做周期性体系的动力学模拟，不是周期性体系也就没有密度一说。GFN1-xTB级别做周期性体系 ...

好的，谢谢社长。借楼再问一个问题，就是xtb在做动力学的时候，设定了初始体系的单电子数，后来体系发生了分解形成自由基，此时单电子数会增加，与开始设置的不符，这种情况是不是没法太较真啊。

sobereva 发表于 Post on 2022-10-19 01:09:11

luzujia 发表于 2022-10-19 00:04
社长，如果用的是xtb做动力学，各种属性怎么输出啊。我现在只知道轨迹文件里输出有能量，利用vmd可以统计氢 ...

xtb本来就不适合做周期性体系的动力学模拟，不是周期性体系也就没有密度一说。GFN1-xTB级别做周期性体系动力学模拟推荐用CP2K

luzujia 发表于 Post on 2022-10-19 00:04:51

社长，如果用的是xtb做动力学，各种属性怎么输出啊。我现在只知道轨迹文件里输出有能量，利用vmd可以统计氢键和RMSD等，其它的属性如密度之类的怎么看

yangzaia 发表于 Post on 2022-9-24 20:16:02

sobereva 发表于 2022-9-24 00:57
得结合具体问题
如果这是个生物分子，再跑个20ns，红线还是平稳的话就可以断定稳定了

是生物大分子，纤维素酶，谢谢sob老师

sobereva 发表于 Post on 2022-9-24 00:57:53

yangzaia 发表于 2022-9-23 15:01
老师，这是我的体系RMSD，我理解的是图中37ns后基本稳定了，但是我不知道如何判断50ns的模拟时间是否足够 ...

得结合具体问题
如果这是个生物分子，再跑个20ns，红线还是平稳的话就可以断定稳定了

yangzaia 发表于 Post on 2022-9-23 15:01:23

sobereva 发表于 2022-9-14 16:10
没图没法说
通常怎么判断在当前帖子里已经说得没法更具体了

老师，这是我的体系RMSD，我理解的是图中37ns后基本稳定了，但是我不知道如何判断50ns的模拟时间是否足够

yangzaia 发表于 Post on 2022-9-23 14:54:53

大头攒毛发表于 2022-9-14 17:01
是不是可以看看37ns前后构象之间的区别

我现在是结合轨迹动画看的，谢谢你的回复

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