计算化学公社

 找回密码 Forget password
 注册 Register
Views: 127|回复 Reply: 4
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] 蛋白质配体体系RMSD图像异常

[复制链接 Copy URL]

15

帖子

0

威望

45

eV
积分
60

Level 2 能力者

如题,RMSD显示的图像呈现进行性下降,md.mdp文件已附录(模拟过程中先完成了没有energygrps的mdrun,完成后又加上了energygrps重跑mdrun),感谢各位老师


md.mdp

2.19 KB, 下载次数 Times of downloads: 1

251

帖子

0

威望

426

eV
积分
677

Level 4 (黑子)

A Student

2#
发表于 Post on 2024-9-17 10:58:14 | 只看该作者 Only view this author
本帖最后由 student0618 于 2024-9-17 11:01 编辑

gmx rms -s 用了哪个结构? 看上去是用了最后一个Frame作reference算RMSD。
Yet to be strong in theory, yet to have enough practical skills.
Still I am having fun with MD simulation.

15

帖子

0

威望

45

eV
积分
60

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2024-9-17 11:19:55 | 只看该作者 Only view this author
本帖最后由 一条空空的河 于 2024-9-17 11:25 编辑
student0618 发表于 2024-9-17 10:58
gmx rms -s 用了哪个结构? 看上去是用了最后一个Frame作reference算RMSD。

您好!-s用的是加跑能量组后的md.tpr,“看上去是用了最后一个Frame作reference算RMSD”这句没太听懂,能详细说一下吗?

251

帖子

0

威望

426

eV
积分
677

Level 4 (黑子)

A Student

4#
发表于 Post on 2024-9-17 12:03:08 | 只看该作者 Only view this author
本帖最后由 student0618 于 2024-9-17 12:12 编辑

RMSD是对比一个参考帧(reference)的值。你的图看上去用作reference的结构是轨迹的最后一帧,因此RMSD会越来越小,到最后因为用了同一结构RMSD数值便会是零了。
RMSD更详尽的解说可参考社长的帖子http://bbs.keinsci.com/thread-27122-1-1.html

你的md.tpr 是否用了跑完md的gro?gmx rms -s 试一下用模拟开始的tpr吧。
Yet to be strong in theory, yet to have enough practical skills.
Still I am having fun with MD simulation.

15

帖子

0

威望

45

eV
积分
60

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2024-9-17 12:48:36 | 只看该作者 Only view this author
本帖最后由 一条空空的河 于 2024-9-17 12:53 编辑
student0618 发表于 2024-9-17 12:03
RMSD是对比一个参考帧(reference)的值。你的图看上去用作reference的结构是轨迹的最后一帧,因此RMSD会越来 ...


您好!确实是用的md01.gro(不带能量组跑出的gro)做的md02.tpr,再用md02.tpr作的rmsd。顺带问一个问题,重跑能量组md02时要用npt.gro(即md01用到的gro)?附带:使用md01.tpr(无能量组mdrun得出的)作的rmsd

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-25 14:31 , Processed in 0.312651 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list