计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 分子模拟 (Molecular Modeling) › 求助为何gromacs计算的MSD线出现转折以及扩散常数偏差太 ...

求助为何gromacs计算的MSD线出现转折以及扩散常数偏差太大？

查看数: 11292 | 评论数: 3 | 收藏 Add to favorites 1

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

帖子模式

背影里的微笑

发布时间: 2020-10-3 14:16

正文摘要:

本帖最后由背影里的微笑于 2020-10-3 14:19 编辑小弟想计算一下重油分子在超临界水中的扩散常数（体系中只有一个油分子，其余都是水分子），但发现在拟合好的MSD并非是一条斜线，而是发生了弯曲，并且扩散常数 ...

回复 Reply

sobereva 发表于 Post on 2021-1-16 08:19:16

背影里的微笑发表于 2021-1-9 18:52
多谢解答，我说的扩散常数偏差是指D[ SA1] = 83.6452 (+/- 22.9934) (1e-5 cm^2/s)括号中的正 ...

误差太大。
你让gmx msd只针对msd的2000~7000 ps比较平直的那段拟合就不会有那么大误差

背影里的微笑 发表于 Post on 2021-1-9 18:52:39

bobosiji 发表于 2020-10-7 09:48
就你展示的msd plot而言，拟合0-7000ps的数据求D应该是可以的。
不过你的系统只有一个重油分子，得多跑（ ...

多谢解答

，我说的扩散常数偏差是指D[ SA1] = 83.6452 (+/- 22.9934) (1e-5 cm^2/s)括号中的正负偏差

bobosiji 发表于 Post on 2020-10-7 09:48:33

本帖最后由 bobosiji 于 2021-1-9 21:42 编辑

就你展示的msd plot而言，拟合3000-7000ps的数据求D应该是可以的。
不过你的系统只有一个重油分子，得多跑（比如10+次；或者跑很长时间比如几百ns）模拟求
D的平均值才比较可信。
另外“扩散常数的偏差也特别大”是啥意思？

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