24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1000  |  回复: 12
当前主题已经存档。

sijifengsd

银虫 (小有名气)

[交流] 【求助】Gromacs 模拟小分子出来怎么办?

请教:

  我在用Gromacs做动力学模拟的时候,发现随着模拟的进行,小分子一短时间跑到水盒子外面,一会又回来了重复了好几个ns;而且,最后,发现大分子也有一部分露在盒子的外面。不知道怎么办,恳请各位高手指点迷津。

[ Last edited by wuli8 on 2009-11-13 at 21:40 ]
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

似水年华
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuhongaaa

金虫 (正式写手)


sijifengsd(金币+1):谢谢参与
使用了周期性边界条件么?
用周期性边界条件是不会有这种问题的,在就是把盒子做大一点。。。。
分子模拟的主页http://varmilion.tk/
2楼2009-09-25 18:17:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wally8962

木虫 (著名写手)

★ ★
sijifengsd(金币+1):谢谢参与
yjcmwgk(金币+1,VIP+0): 9-26 13:06
周期边界就会这样。不奇怪。

任何一个分子,包括水,当运动到盒子边缘时,都会一部分在盒内,一部分在盒外。实际上就是到盒子的另一侧去了。
3楼2009-09-25 22:16:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sijifengsd

银虫 (小有名气)

谢谢

引用回帖:
Originally posted by zhuhongaaa at 2009-9-25 18:17:
使用了周期性边界条件么?
用周期性边界条件是不会有这种问题的,在就是把盒子做大一点。。。。

我使用了周期性边界条件,设定盒子边缘到分子的值为 -c -d 0.9 。我的大分子和小分子不是一起动的,有时候出现小分子脱离大分子的情况,这个怎么解决啊?
似水年华
4楼2009-09-26 10:06:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sijifengsd

银虫 (小有名气)

请教

引用回帖:
Originally posted by wally8962 at 2009-9-25 22:16:
周期边界就会这样。不奇怪。

任何一个分子,包括水,当运动到盒子边缘时,都会一部分在盒内,一部分在盒外。实际上就是到盒子的另一侧去了。

这样怎么解决呢?当大分子或者小分子跑到外边的时候,做它们的rmsd时候结果偏离很大。那结果还能用么?
似水年华
5楼2009-09-26 10:10:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wally8962

木虫 (著名写手)

★ ★
sijifengsd(金币+1,VIP+0):感谢回复 9-26 11:04
yjcmwgk(金币+1,VIP+0): 9-26 13:06
周期边界,你得先把坐标处理一下。
比如周期盒范围是0-10,你原始坐标为0.5,后来跑到周期盒另一边去了,坐标为9.5。考虑到周期盒的因素,它的坐标实际上是9.5-10=-0.5。
ok
6楼2009-09-26 10:27:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sijifengsd

银虫 (小有名气)

请教

引用回帖:
Originally posted by wally8962 at 2009-9-26 10:27:
周期边界,你得先把坐标处理一下。
比如周期盒范围是0-10,你原始坐标为0.5,后来跑到周期盒另一边去了,坐标为9.5。考虑到周期盒的因素,它的坐标实际上是9.5-10=-0.5。
ok

非常感谢,我还的请教您。您说的是在Gromacs 的这一步改动么:editconf -bt cubic -f *.gro -o *.gro -c -d 0.9改动后面的数字么?
似水年华
7楼2009-09-26 10:40:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sijifengsd

银虫 (小有名气)

请教

引用回帖:
Originally posted by wally8962 at 2009-9-26 10:27:
周期边界,你得先把坐标处理一下。
比如周期盒范围是0-10,你原始坐标为0.5,后来跑到周期盒另一边去了,坐标为9.5。考虑到周期盒的因素,它的坐标实际上是9.5-10=-0.5。
ok

我可能理解的有问题,我的问题在于小分子和大分子在md的时候分离的比较远,不形成复合物了,直接导致RMSD偏离较大。我一共跑了30ns的动力学,前20个ns的时候小分子和大分子还是一起的;但是到最后发现小分子和大分子分离了,并且小分子离开了盒子:您看还是周期性的问题么?
还能用trjconv解决么?
似水年华
8楼2009-09-26 11:03:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wally8962

木虫 (著名写手)

引用回帖:
Originally posted by sijifengsd at 2009-9-26 11:03:

我可能理解的有问题,我的问题在于小分子和大分子在md的时候分离的比较远,不形成复合物了,直接导致RMSD偏离较大。我一共跑了30ns的动力学,前20个ns的时候小分子和大分子还是一起的;但是到最后发现小分子和大 ...

gromacs我不熟。editconf应该是模拟前设置系统的,事后分析怎么用我就不知道了。

我只能说是那个原理。至于gromacs具体怎么处理就不清楚了。
9楼2009-09-26 13:37:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuhongaaa

金虫 (正式写手)

★ ★
sijifengsd(金币+1,VIP+0):谢谢 9-27 08:48
yjcmwgk(金币+1,VIP+0): 10-1 22:35
引用回帖:
Originally posted by sijifengsd at 2009-9-26 11:03:

我可能理解的有问题,我的问题在于小分子和大分子在md的时候分离的比较远,不形成复合物了,直接导致RMSD偏离较大。我一共跑了30ns的动力学,前20个ns的时候小分子和大分子还是一起的;但是到最后发现小分子和大 ...

感觉这样的话可能是体系的问题,相互作用能太小?如果结合的较紧的话是不会分开的。。
分子模拟的主页http://varmilion.tk/
10楼2009-09-26 18:39:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sijifengsd 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +3 yexuqing 2026-04-19 3/150 2026-04-20 11:29 by jurkat.1640
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 5/250 2026-04-20 10:47 by YuY66
[考研] 通信工程求调剂!!! +7 zlb770521 2026-04-14 7/350 2026-04-19 20:56 by Equinoxhua
[考博] 申博自荐 +4 Linxia林夏 2026-04-13 4/200 2026-04-19 19:55 by Equinoxhua
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 291求调剂 +11 关忆北. 2026-04-14 11/550 2026-04-19 17:16 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 307中医考研调剂 +9 于以采蘩 2026-04-14 9/450 2026-04-19 08:41 by 烟雨流涯
[考研] 327求调剂 +27 Xxjc1107. 2026-04-13 30/1500 2026-04-19 08:22 by cuisz
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 0854求调剂 +23 门路摸摸 2026-04-15 27/1350 2026-04-19 01:59 by 烟雨流涯
[考研] 300求调剂 +12 橙a777 2026-04-15 12/600 2026-04-18 23:51 by 路病情
[考研] 320求调剂 +5 深郊akm 2026-04-17 5/250 2026-04-18 19:52 by 王珺璞
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[考研] 一志愿华中农业071010,320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 一志愿沪9,生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
[考研] 药学求调剂 +14 喽哈加油 2026-04-14 16/800 2026-04-16 10:15 by beilsong20
[考研] 考研调剂 +13 长弓傲 2026-04-13 14/700 2026-04-14 14:44 by zs92450
信息提示
请填处理意见