计算化学公社

 找回密码 Forget password
 注册 Register

gromacsr软件模拟后,周期性设置的影响?

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

正文摘要:

本帖最后由 努力向上的boy 于 2021-10-10 14:54 编辑 现在对于周期性带来的影响不是很明白现在呢,先讲一下我的体系:是这个样子的,我的体系是溶质(聚合物晶体,聚合物由一条条相同的链组成的)+溶剂(小分子) ...

回复 Reply

努力向上的boy 发表于 Post on 2021-10-28 19:15:39
sobereva 发表于 2021-10-22 02:58
VMD里把相邻镜像显示出来,一看便知

老师您的意思是说,用-pbc nojump调整的轨迹,再用gmx trjconv提取之后的gro不是按照-pbc nojump合适-pbc mol调整的?
sobereva 发表于 Post on 2021-10-22 02:58:22
努力向上的boy 发表于 2021-10-21 16:54
请问老师,这个该怎么判断呢?判断是否合理是否与镜像接触了?

VMD里把相邻镜像显示出来,一看便知
努力向上的boy 发表于 Post on 2021-10-21 16:54:51
sobereva 发表于 2021-10-20 02:30
gro文件默认是按照-pbc mol方式记录的轨迹,只有当分子质心越过盒子时才会被整体挪到另一头

请问老师,这个该怎么判断呢?判断是否合理是否与镜像接触了?
sobereva 发表于 Post on 2021-10-20 02:30:18
努力向上的boy 发表于 2021-10-20 01:35
谢谢老师。
老师,我想请教一下个问题,在模拟的时候,盒子的物质,不是从一边跑出盒子就会从另外一边跑 ...

gro文件默认是按照-pbc mol方式记录的轨迹,只有当分子质心越过盒子时才会被整体挪到另一头
努力向上的boy 发表于 Post on 2021-10-20 01:35:11
sobereva 发表于 2021-10-15 07:29
可能只是相对于整体结构而言的很微小分散,导致回转半径数值上体现不出来明显差异
可以从X/Y/Z方向最大 ...

谢谢老师。
老师,我想请教一下个问题,在模拟的时候,盒子的物质,不是从一边跑出盒子就会从另外一边跑进来吗,从原始轨迹提取的结构中,几条链相邻,一点残基在盒子外面,但是为什么利用-pbc nojump处理轨迹后,会有几条链在盒子内且相邻,另外几条链分别在盒子相邻的不同的镜像盒子中,这是为什么呢
sobereva 发表于 Post on 2021-10-15 07:29:46
努力向上的boy 发表于 2021-10-14 19:29
1、前边的聚合物轨迹总的来说还是保持相对规整,但是边端已经有明显的分散,层间或者链间还是相对紧密, ...

可能只是相对于整体结构而言的很微小分散,导致回转半径数值上体现不出来明显差异
可以从X/Y/Z方向最大减最小坐标衡量一下跨距的变化
努力向上的boy 发表于 Post on 2021-10-14 19:29:55
sobereva 发表于 2021-10-14 13:45
1 结合动力学图像判断是否正常

2、3 a中任意一个原子都必须和镜像的a的任意一个原子距离间距超过截断 ...

1、前边的聚合物轨迹总的来说还是保持相对规整,但是边端已经有明显的分散,层间或者链间还是相对紧密,请问老师还有别的方法辅助判断一下吗?
sobereva 发表于 Post on 2021-10-14 13:45:57
努力向上的boy 发表于 2021-10-14 11:56
1、看起来是一条水平的线,数值的话在不断的变化,开始大概是1.65,然后在左右徘徊,然后是1.64,然后在 ...

1 结合动力学图像判断是否正常

2、3 a中任意一个原子都必须和镜像的a的任意一个原子距离间距超过截断半径
原则上a也不应当遇到其它的链的镜像,毕竟和实际情况不符。但如果a运动到接触其它的链的镜像之前,被考察的问题已经(如溶胀过程)已经足够充分体现了,就不用管之后的轨迹
努力向上的boy 发表于 Post on 2021-10-14 11:56:20
sobereva 发表于 2021-10-14 00:52
不可能是严格水平的
看看具体数值怎么变的

1、看起来是一条水平的线,数值的话在不断的变化,开始大概是1.65,然后在左右徘徊,然后是1.64,然后在左右浮动,然后1.63,然后1.65,这样的,请问老师,这个正常吗?
2、老师,针对个溶质,它不是有好几条链吗,比如有三条,每条链分别叫做a、b、c,在做模拟的时候溶质要和镜像的距离大于截断半径,是不是说a要a的镜像大于截断半径,并且是a中的相应原子之间的距离大于截断半径?
3、如果a中的原子1与镜像中a上的2、3、4等原子的距离小于截断半径,或者a链与c链有些原子小于截断半径,这都是没问题的是吧?
sobereva 发表于 Post on 2021-10-14 00:52:36
努力向上的boy 发表于 2021-10-13 19:02
老师我想请教一下,利用gmx gyrate计算回旋半径,开始的时候,是规整排列的有序的,但是结构已将发生了变 ...

不可能是严格水平的
看看具体数值怎么变的
sobereva 发表于 Post on 2021-10-13 04:35:10
努力向上的boy 发表于 2021-10-13 02:07
2、老师,如果不对轨迹进行处理,溶质有一部分在跑出盒子,然后跑出的那部分从另外一端进入到盒子,这样 ...

当前情况不适合把所有聚合物整体当做一个组用mindist说事,对每一条聚合物链分别用mindist测量它与镜像最短距离才有意义
努力向上的boy 发表于 Post on 2021-10-13 02:07:38
sobereva 发表于 2021-10-11 20:03
1 最好用-nojump,保证聚合物分子运动的连续性,这种情况下跑出去完全不用管,是什么样就是什么样

2  ...

2、老师,如果不对轨迹进行处理,溶质有一部分在跑出盒子,然后跑出的那部分从另外一端进入到盒子,这样在利用gmx mindist -f *.xtc -s *.tpr -n *.ndx -pi计算溶质和镜像的最短距离的时候,距离就会显示的很小,这怎么办呢?
3、我主要是观察溶质由紧密的结构到了松散分开的结构的程度
sobereva 发表于 Post on 2021-10-11 20:03:57
努力向上的boy 发表于 2021-10-11 11:04
溶剂里面有水有无机盐有有机小分子,一共三种组分
1、利用-pbc -center -fit等处理轨迹后,如果还有部分原 ...

1 最好用-nojump,保证聚合物分子运动的连续性,这种情况下跑出去完全不用管,是什么样就是什么样

2 不用

3 合理
可加可不加。如果你对绕着不同主轴的回转半径分量感兴趣,加了便于考察

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

GMT+8, 2026-2-23 16:56 , Processed in 0.187323 second(s), 25 queries , Gzip On.

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