计算化学公社

 找回密码 Forget password
 注册 Register

Gromacs跑有机物聚合过程,为什么有的跑散了

查看数: 14221 | 评论数: 10 | 收藏 Add to favorites 2
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2021-6-18 10:44

正文摘要:

本帖最后由 wudiazhu 于 2021-6-18 10:48 编辑 用acpype生成了我聚合物的itp文件,并用RESP电荷替代,用packmol构建50个我的聚合物和4w个水,能量极小化后用Berendsen压浴跑了5ns的动力学,想看我的聚合物在水环 ...

回复 Reply

zdworld 发表于 Post on 2023-8-11 09:54:06
你这是周期性盒子,把xyz方向的周期性镜像都开起来看看。
sobereva 发表于 Post on 2023-7-30 12:41:30
a-Student 发表于 2023-7-30 10:06
老师,打扰您一下,遇到这种情况请教您:
1. 如图1所示,GMX跑完查看最后一帧轨迹原本是一个有n=100组成 ...

1 本来就该这样。你画箭头的地方是被聚合物占的地方,自然没有水

2 先把轨迹弄成图2那样让聚合物保持完整,之后trjconv结合-center让选定的组居中,组选择聚合物部分
a-Student 发表于 Post on 2023-7-30 10:06:15
sobereva 发表于 2021-6-18 18:49
你这聚合物本来亲水性就不弱,没理由期望一定会堆到一起去

老师,打扰您一下,遇到这种情况请教您:
1. 如图1所示,GMX跑完查看最后一帧轨迹原本是一个有n=100组成的长链,跑完发现长链被断开了,然后我用命令gmx trjconv -f md.xtc -s md.tpr -pbc whole -o md.xtc优化了轨迹保持分子完整性,图2分子长链又变成聚集状态了,但发现水盒子在原来散开的聚合物片段位置出现了水分子缺失,水盒子看上去缺一块,老师这种情况我做了轨迹优化会影响结果分析吗?
2. 因为聚合物他跑出盒子了,所以我在gmx trjconv -f md.xtc -s md.tpr -pbc whole -o md.xtc生成的.xtc文件基础上又做了将聚合物放在盒子中心操作,gmx trjconv -f md.xtc -s md.tpr -pbc mol -o md.xtc,发现位置没有变化,请问老师怎么才能把它放在盒子里呢?

麻烦您啦,谢谢指教。

202307301005374131..png (236.25 KB, 下载次数 Times of downloads: 18)

图1

图1

202307301004183144..png (882.25 KB, 下载次数 Times of downloads: 16)

图2

图2
wudiazhu 发表于 Post on 2021-6-19 20:41:57
wbn 发表于 2021-6-18 17:51
那就不能保证他们全聚集在一起了。你的初始状态是聚集的吗?

初始状态不是聚集的,只是用packmol均匀分布在盒子中
wudiazhu 发表于 Post on 2021-6-19 20:41:17
sobereva 发表于 2021-6-18 18:49
你这聚合物本来亲水性就不弱,没理由期望一定会堆到一起去

好吧,谢谢sob老师,我再跑长时间看看吧
sobereva 发表于 Post on 2021-6-18 18:49:56
你这聚合物本来亲水性就不弱,没理由期望一定会堆到一起去
wbn 发表于 Post on 2021-6-18 17:55:59
你的第一个问题,想看聚合物在水环境中的堆叠过程,那怎么能去掉水。
第二个问题,你的体系都是小分子,你看镜像的距离干吗
wbn 发表于 Post on 2021-6-18 17:51:44
wudiazhu 发表于 2021-6-18 15:10
是50个聚合物,一个聚合物是4个单元的链

那就不能保证他们全聚集在一起了。你的初始状态是聚集的吗?
wudiazhu 发表于 Post on 2021-6-18 15:10:58
wbn 发表于 2021-6-18 11:48
盒子里是由50个聚合物还是1个50单体的聚合物链?

是50个聚合物,一个聚合物是4个单元的链
wbn 发表于 Post on 2021-6-18 11:48:29
盒子里是由50个聚合物还是1个50单体的聚合物链?

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

GMT+8, 2026-2-24 03:57 , Processed in 0.180803 second(s), 25 queries , Gzip On.

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