• 回答数

    4

  • 浏览数

    1802

  • 收藏数

    0

作者:一抹烟红 发表于 2018-11-26 15:00:28
跳转到指定楼层
请教大家,最近用gromcs计算扩散系数,使用命令为g_msd -f md.trr s md.tpr -o msd.xvg,整个体系为一个纯体系,包括5000个原子,直接用上述命令求出来的扩散系数可信吗?看相关帖子说拟合的时候需要仔细选取拟合的起点和终点,那么整个起点和终点怎么确定?还有manual中提到-trestart ,The time between additional starting points for the MSD calculation is set with -trestart 意思是 -trestart 标定了起始点的时间间隔。 请问这里的起始点(starting points)什么意思?是不是计算msd的初始位置点? g_msd 是不是只计算两个起始点时间间隔内的msd,然后在绘制一条组合的msd曲线(上一间隔的终点为下一间隔的起点,问题较多,请各位帮忙解答!谢谢
分享:
回复

使用道具

该用户从未签到

新手上路

Rank: 1

积分
30
极客币
51
主题
8
帖子
23
注册时间
2018-10-31
在线时间
1 小时
性别
保密
发表于 2018-11-26 15:00:40 | 显示全部楼层
结果是可信的,就是要选取平衡后的时间段,trestart是每隔多长时间取点计算msd
回复

使用道具 举报

该用户从未签到

新手上路

Rank: 1

积分
31
极客币
55
主题
9
帖子
26
注册时间
2018-10-31
在线时间
1 小时
性别
保密
 楼主| 发表于 2018-11-26 15:00:47 | 显示全部楼层
谢谢,那意思是我可以先做出体系的rmsd,确定平衡的时间段,然后再这个时间段用g_msd -f xx.trr -s xx.tpr -beginfitt1 -endfit t2-o msd.xvg-trestart 100 去求扩散系数?终端及msd文件里的D即是要求的扩散系数吧。谢谢了
回复

使用道具 举报

该用户从未签到

新手上路

Rank: 1

积分
30
极客币
51
主题
8
帖子
23
注册时间
2018-10-31
在线时间
1 小时
性别
保密
发表于 2018-11-26 15:00:59 | 显示全部楼层
是的      
回复

使用道具 举报

该用户从未签到

新手上路

Rank: 1

积分
31
极客币
55
主题
9
帖子
26
注册时间
2018-10-31
在线时间
1 小时
性别
保密
 楼主| 发表于 2018-11-26 15:01:05 | 显示全部楼层
回复

使用道具 举报

高级模式 评论
您需要登录后才可以回帖 登录 | 立即注册 微信登录
关于作者
一抹烟红

用户组:新手上路

  • 主题

    9

  • 帖子

    26

  • 关注者

    1