计算化学公社

 找回密码 Forget password
 注册 Register

为什么MD中用gromacs的gmx rdf命令和vmd算的两个组分之间的配位数不一样

查看数: 15604 | 评论数: 12 | 收藏 Add to favorites 4
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2020-4-24 02:59

正文摘要:

本帖最后由 wgg1181543722 于 2020-4-24 13:30 编辑 老师,我算的体系是浸在水中的CNT,选取CNT中的 钙离子与氧,求Ca-O之间的配位数,我用gmx算出来的是5。而vmd算出来的是6。而通过vmd动图观察,钙离子全程与6 ...

回复 Reply

yl233 发表于 Post on 2023-11-3 20:22:12
wgg1181543722 发表于 2021-6-22 22:35
我是弄了半天,我自己整错了,你还是这个问题吗

你好,请问是什么原因呀,最终选择了哪个呀?
wgg1181543722 发表于 Post on 2021-6-22 22:35:25
amani73986 发表于 2021-6-4 12:49
楼主好,我也遇到了一样的问题,用gmx rdf计算出来的配位数总是偏小,请问找到原因了么。谢谢

我是弄了半天,我自己整错了,你还是这个问题吗
amani73986 发表于 Post on 2021-6-4 12:49:08
楼主好,我也遇到了一样的问题,用gmx rdf计算出来的配位数总是偏小,请问找到原因了么。谢谢
wgg1181543722 发表于 Post on 2020-4-24 16:52:29
liuyuje714 发表于 2020-4-24 15:44
看起来是那个峰高不一致导致的结果不同。我看你这rdf很粗糙啊,体系平衡好了么?帧数足够么?

是  很奇怪,峰的高度不一样,但最高峰的位置是一致的。因为太大了,我选取的是最后30ns,帧数有点少,所以 不是很光滑
liuyuje714 发表于 Post on 2020-4-24 15:44:30
wgg1181543722 发表于 2020-4-24 15:25
因为CNT是收尾相连的,管中的水分子是出不来的,index的编号和gro里的是一直对应的 也没变,配位数不是算 ...

看起来是那个峰高不一致导致的结果不同。我看你这rdf很粗糙啊,体系平衡好了么?帧数足够么?
liuyuje714 发表于 Post on 2020-4-24 14:37:08
本帖最后由 liuyuje714 于 2020-4-24 15:09 编辑

你确定你的gmx定义的index弄对了吗?CNT中的水分子你最好用gmx select -on进行选择,并且查看输出的index.ndx中不同时刻的水分子OW编号是否没有变过,即一直在CNT中。还有你说的5或者6是怎么计算得到的也要说一下,因为gmx和vmd好像并没有直接给出配位数是多少。我测试了一下CNT中的水,用gmx和vmd计算得到的一致。
wgg1181543722 发表于 Post on 2020-4-24 13:38:51
本帖最后由 wgg1181543722 于 2020-4-24 13:44 编辑
sobereva 发表于 2020-4-24 04:17
没具体信息不好说。VMD里你用了动态选区,而gmx rdf你用的是静态选区,考虑这个是否真的能对应上
并且尝试 ...

老师 我选取的这部分区域 是密封的,CNT内的水分子与钙是不会从里面出来的。这两个选取的区域 不都是动态的吗 怎么看静态 动态
wgg1181543722 发表于 Post on 2020-4-24 13:00:36
sobereva 发表于 2020-4-24 04:17
没具体信息不好说。VMD里你用了动态选区,而gmx rdf你用的是静态选区,考虑这个是否真的能对应上
并且尝试 ...

我算的是浸在水中的CNT中的水分子与钙离子,选取的组分都是一样的。我算了CA-OW,CL-OW,CA-CL,OW-OW四种相互作用,刚才对比了一下算出来的结果,rdf是完全重合的,唯一的不同,就是配位数,vmd算出来的要比gmx算出来的大一些。
wgg1181543722 发表于 Post on 2020-4-24 12:40:13
liuyuje714 发表于 2020-4-24 09:23
这明显算的不相同。

额 请问一下。两种算法,我都选取的是相同的物质,怎么不一样了?
liuyuje714 发表于 Post on 2020-4-24 09:23:57
这明显算的不相同。
sobereva 发表于 Post on 2020-4-24 04:17:28
没具体信息不好说。VMD里你用了动态选区,而gmx rdf你用的是静态选区,考虑这个是否真的能对应上
并且尝试把rdf曲线图作出来进行对比来分析原因

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

GMT+8, 2026-2-21 15:42 , Processed in 0.178714 second(s), 25 queries , Gzip On.

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