计算化学公社

 找回密码 Forget password
 注册 Register

计算多条相同链的温敏特性,周期性边界问题处理

查看数: 556 | 评论数: 3 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2024-6-22 10:57

正文摘要:

各位老师好,我计算的体系是多条相同链的温敏特性,请问我使用gmx trjconv 对轨迹进行处理代码如下 gmx trjconv -s eq.tpr -f eq.xtc  -o mol-eq.xtc -pbc mol   输出选择0系统 观察轨迹出现 ...

回复 Reply

byj5213344 发表于 Post on 2024-6-23 09:03:56
本帖最后由 byj5213344 于 2024-6-23 09:19 编辑

感谢sob老师的建议,我下次发帖注意,
有个问题,我用-pbc cluster 选择的是非水部分就是我不想被截断的索引组,因为我模拟的体系是水+多条纤维素链,不想其中的纤维素链穿越盒子引起计算纤维素链rmsd时较大的波动,在第二步选择cluster出现好多纤维素链在盒子外面的情况,我第三步进行居中出现了169帧和170帧有一条链穿越盒子的问题,既然我已经-pbc cluster了,为什么居中还会出现穿越盒子的行为?
sobereva 发表于 Post on 2024-6-22 23:26:52
那叫命令,不叫代码,又不是写源代码

搞清楚-pbc cluster是什么意思,哪能整个体系当做一个cluster,此时和-nojump效果几乎一样

把你不想被盒子截断的部分一起定义为一个索引组,用-pbc cluster时选上这个组
sobereva 发表于 Post on 2024-6-22 23:09:57
不要自己在标体系里手写[GROMACS],给你去了,以后注意!

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

GMT+8, 2024-11-24 00:58 , Processed in 0.182559 second(s), 26 queries , Gzip On.

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