计算化学公社

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

[GROMACS] 求助:如何解决多条纤维素链体系RMSD图的波动问题?

[复制链接 Copy URL]

25

帖子

0

威望

475

eV
积分
500

Level 4 (黑子)

各位老师好,我的动力学体系中包含5条纤维素链,我想要做一下纤维素的RMSD,结果图像出现巨大的突变。我参考sob老师的博文,输入了以下命令,结果仍然没有改变。
第一步使分子保持完整
gmx trjconv -s md.tpr -f md.xtc -o md_mol.xtc -pbc mol
第二步让指定的组在轨迹中一直保持完整
gmx trjconv -s md.tpr -f md_mol.xtc -o md_mol_cluster.xtc -pbc cluster
最后再用命令将纤维素链进行叠合
gmx trjconv -s md.tpr -f md_mol_cluster.xtc -o new.xtc -fit rot+trans
最后将得到new.xtc文件来做RMSD,结果如图。
RMSD图任然是这种上上下下的状态,我想知道该如何解决这个问题?

rmsd.png (15.83 KB, 下载次数 Times of downloads: 6)

rmsd.png

310

帖子

0

威望

1492

eV
积分
1802

Level 5 (御坂)

2#
发表于 Post on 2023-10-20 18:17:49 | 只看该作者 Only view this author
多条链的情况有时候不适合用-pbc mol,因为mol要求成键的原子距离不能太远,这样只能保持一条链不会断开,因为不同链之间没有成键。
这种情况可以使用-pdc nojump,-s选择一个链没有断开的结构文件。

25

帖子

0

威望

475

eV
积分
500

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2023-10-20 18:41:12 | 只看该作者 Only view this author
Frozen-Penguin 发表于 2023-10-20 18:17
多条链的情况有时候不适合用-pbc mol,因为mol要求成键的原子距离不能太远,这样只能保持一条链不会断开, ...

问题已解决,感谢。

25

帖子

0

威望

173

eV
积分
198

Level 3 能力者

4#
发表于 Post on 7 day ago | 只看该作者 Only view this author
huashang 发表于 2023-10-20 18:41
问题已解决,感谢。

老师,请问最后是怎么解决的哇,遇到了相同的困惑

本版积分规则 Credits rule

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

GMT+8, 2025-8-12 23:08 , Processed in 0.166721 second(s), 23 queries , Gzip On.

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