计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 分子模拟 (Molecular Modeling) › 计算一个胶束中的表面活性剂在溶液中扩散系数时，gmx MS ...

计算一个胶束中的表面活性剂在溶液中扩散系数时，gmx MSD曲线上出现拐点怎么办

查看数: 401 | 评论数: 5 | 收藏 Add to favorites 0

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

帖子模式

含光君

发布时间: 2025-3-25 09:31

正文摘要:

大家好，我最近在计算一个胶束中的表面活性剂在溶液中扩散系数时，发现MSD曲线上出现拐点且两边斜率不等，可能是因为模拟中胶束由盒子中心运动至盒子边界并穿越了边界。所用命令为： gmx msd -f md.xtc -s m ...

回复 Reply

JillW 发表于 Post on 2025-4-7 10:06:18

含光君发表于 2025-4-1 10:12
谈不上赐教互相学习

我这里确实不是穿越盒子的问题，按照卢老师说的用cluster处理后得到的msd图 ...

嗯嗯，是的，我也几乎把论坛里MSD相关的都看了一遍，感觉计算整段MSD再对前半段做拟合的误差会比较小

当时也看到卢老师在别的帖子底下回复说如果小分子个数比较少的话，就应当尽可能延长统计时间，因此我总时间设置的也比较长，一定程度上减少误差。

感谢回复交流

含光君 发表于 Post on 2025-4-1 10:12:38

JillW 发表于 2025-3-31 09:20
请问您结果怎么样呀？我发现我的体系在很长时间尺度下如1000ns时前半段会有很明显的线性段，但是我将msd ...

谈不上赐教互相学习

我这里确实不是穿越盒子的问题，按照卢老师说的用cluster处理后得到的msd图是一样的。

我觉得你的问题里小分子数量可能对统计的误差比较大，另外，从MSD的计算来看（见http://bbs.keinsci.com/thread-42063-1-1.html帖子下卢老师回复的PPT）, 只统计前半段和统计整段在取平均的区段上不是完全等同的，所以可能会有影响。如果你想计算前半段的D, 可以计算整段MSD再对前半段做拟合。

即 gmx msd -beginfit XXX -endfit XXX

JillW 发表于 Post on 2025-3-31 09:20:26

含光君发表于 2025-3-25 23:14
原来是这样，我来分段计算试试，非常感谢sob老师！

ps, 谢谢社长帮忙修改帖子标题

请问您结果怎么样呀？我发现我的体系在很长时间尺度下如1000ns时前半段会有很明显的线性段，但是我将msd计算范围缩小到一半时，线性段就消失了

由于后半段统计误差会比较大，是不是时间尺度长的时候取前半段拟合结果更可信呢？
期待您的回答，感谢不吝赐教

ps：体系中的小分子数量比较少，仅有1个

含光君 发表于 Post on 2025-3-25 23:14:54

本帖最后由含光君于 2025-3-25 23:17 编辑

sobereva 发表于 2025-3-25 23:11
这和穿越盒子没关系，gmx msd会自动考虑PBC。如果你就是非常怀疑是穿越盒子的问题，也可以trjconv结合-clus ...

gmx msd会自动考虑PBC

原来是这样，我来分段计算试试，非常感谢sob老师！

ps, 谢谢社长帮忙修改帖子标题

sobereva 发表于 Post on 2025-3-25 23:11:17

这和穿越盒子没关系，gmx msd会自动考虑PBC。如果你就是非常怀疑是穿越盒子的问题，也可以trjconv结合-cluster处理轨迹让簇的运动保持完整，然后再算算
可以尝试对轨迹不同区段计算MSD，试图搞清楚是什么原因

评分 Rate

参与人数 Participants 1	eV +3	收起理由 Reason
含光君	+ 3	谢谢

查看全部评分 View all ratings

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