24小时热门版块排行榜    

查看: 976  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 312求调剂 +3 陌宸希 2026-03-16 4/200 2026-03-16 15:06 by peike
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 考研材料与化工,求调剂 +8 戏精丹丹丹 2026-03-09 8/400 2026-03-14 01:14 by JourneyLucky
[考研] 调剂 +3 13853210211 2026-03-10 3/150 2026-03-14 00:47 by JourneyLucky
[考研] 341求调剂 +4 番茄头--- 2026-03-10 4/200 2026-03-13 23:12 by JourneyLucky
[考研] 285 求调剂 资源与环境 一志愿北京化工大学 +3 未名考生 2026-03-10 3/150 2026-03-13 23:04 by JourneyLucky
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] 一志愿中科院,化学方向,295求调剂 +4 一氧二氮 2026-03-11 4/200 2026-03-13 22:35 by JourneyLucky
[考研] 求调剂(材料与化工327) +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] (081700)化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 工科278分求调剂 +5 周慢热啊 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
[考研] 土木第一志愿276求调剂,科研和技能十分丰富,求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
信息提示
请填处理意见