24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 4588  |  回复: 19

转基因猴子

木虫 (正式写手)

[求助] gromacs蛋白计算结束后老是跑出盒子 已有3人参与

最近使用gromacs对我的目的蛋白进行分子模拟。box使用cubic,直径为1.0,在350K下50ns蛋白跑出了盒子。然后把直径调整为1.5nm,50ns跑完后蛋白依旧跑出了盒子。这是为何呀?是不是我参数哪里设错了?
下面附上我最后MD时的参数:
title               =  Yo
cpp                 =  /usr/bin/cpp
constraints         =  all-bonds
integrator          =  md
dt                  =  0.002        ; ps !
nsteps              =  25000000        ; total 10 ps.
nstxout             =  2000
nstvout             =  2000
nstfout             =  0
nstlog              =  2000
nstenergy           =  2000
ns_type             =  grid
rlist               =  1.0
coulombtype         =  PME
rcoulomb            =  1.0
vdwtype             =  cut-off
rvdw                =  1.4
fourierspacing      =  0.12
fourier_nx          =  0
fourier_ny          =  0
fourier_nz          =  0
pme_order           =  4
ewald_rtol          =  1e-5
optimize_fft        =  yes
; Berendsen temperature coupling is on in two groups
Tcoupl              =  v-rescale
tc-grps                    =  Protein        non-protein
tau_t               =  0.1        0.1
ref_t               =  450        450
; Energy monitoring
energygrps            =  Protein        non-protein
; Pressure coupling is on
Pcoupl              =  Parrinello-Rahman
Pcoupltype          =  isotropic
tau_p               =  2
compressibility     =  4.5e-5
ref_p               =  1.0
; Generate velocites is on at 450 K.
gen_vel             =  yes
gen_temp            =  450.0
gen_seed            =  173529
回复此楼

» 收录本帖的淘帖专辑推荐

MDs-Gromacs

» 猜你喜欢

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

» 本主题相关价值贴推荐,对您同样有帮助:

我还年轻,我要上路!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xulinan

木虫 (小有名气)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
月只蓝: 金币+1, 鼓励交流! 2014-04-02 16:49:55
转基因猴子: 金币+2, ★★★很有帮助 2014-04-04 07:02:25
跑出去也是正常的,反正有周期性边界条件在,你用trjconv 去周期后就好了,不会影响结果的
2楼2014-04-01 22:13:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

转基因猴子

木虫 (正式写手)

引用回帖:
2楼: Originally posted by xulinan at 2014-04-01 22:13:54
跑出去也是正常的,反正有周期性边界条件在,你用trjconv 去周期后就好了,不会影响结果的

能说详细点吗?如何去周期性条件?

[ 发自小木虫客户端 ]
我还年轻,我要上路!
3楼2014-04-02 06:12:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xulinan

木虫 (小有名气)

引用回帖:
3楼: Originally posted by 转基因猴子 at 2014-04-02 06:12:11
能说详细点吗?如何去周期性条件?
...

trjconv -f ** -s ** -o**-pbc nojump -center  也可以-pbc mol
4楼2014-04-02 06:48:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aileen2010

新虫 (初入文坛)

【答案】应助回帖


感谢参与,应助指数 +1
月只蓝: 金币+1, 鼓励交流! 2014-04-02 16:50:02
还可以使用editconf命令
5楼2014-04-02 09:17:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

转基因猴子

木虫 (正式写手)

引用回帖:
5楼: Originally posted by aileen2010 at 2014-04-02 09:17:54
还可以使用editconf命令

如何使用?请赐教。还有跑出去不影响计算准确性吗?
我还年轻,我要上路!
6楼2014-04-02 13:29:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

转基因猴子

木虫 (正式写手)

引用回帖:
4楼: Originally posted by xulinan at 2014-04-02 06:48:07
trjconv -f ** -s ** -o**-pbc nojump -center  也可以-pbc mol...

非常感谢。请问蛋白跑出盒子影响计算准确性吗?
我还年轻,我要上路!
7楼2014-04-02 13:29:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

转基因猴子

木虫 (正式写手)

引用回帖:
2楼: Originally posted by xulinan at 2014-04-01 22:13:54
跑出去也是正常的,反正有周期性边界条件在,你用trjconv 去周期后就好了,不会影响结果的

加入周期性边界条件是在最后MD中加入参数pbc=full这项吗?
我还年轻,我要上路!
8楼2014-04-02 13:34:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xueer小雯

金虫 (初入文坛)

【答案】应助回帖


感谢参与,应助指数 +1
月只蓝: 金币+1, 鼓励交流! 2014-04-07 10:05:23
引用回帖:
8楼: Originally posted by 转基因猴子 at 2014-04-02 13:34:18
加入周期性边界条件是在最后MD中加入参数pbc=full这项吗?...

跑出盒子没有关系,去除周期性边界就好了。
一般不影响计算精确度。
加周期性边界条件是 pbc=xyz 在三个方向都加
9楼2014-04-02 18:00:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

转基因猴子

木虫 (正式写手)

引用回帖:
9楼: Originally posted by xueer小雯 at 2014-04-02 18:00:09
跑出盒子没有关系,去除周期性边界就好了。
一般不影响计算精确度。
加周期性边界条件是 pbc=xyz 在三个方向都加...

也就是去掉pbc这项参数就可以了吗?

[ 发自小木虫客户端 ]
我还年轻,我要上路!
10楼2014-04-02 21:40:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 转基因猴子 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料调剂 +14 壹贰贰亿 2026-04-04 14/700 2026-04-05 23:31 by 来看流星雨10
[考研] 找调剂 +10 楚乔乔 2026-04-01 10/500 2026-04-05 22:19 by syh9288
[考研] 求调剂 +4 wos666 2026-04-03 4/200 2026-04-05 11:48 by arrow8852
[考研] 285求调剂 +11 哦呦呼o 2026-04-04 11/550 2026-04-05 08:15 by 544594351
[考研] 一志愿郑州大学材料与化工085600,求调剂 +24 吃的不少 2026-04-02 24/1200 2026-04-04 23:20 by 永字号
[考研] 求生物学调剂 +14 15172915737 2026-04-01 14/700 2026-04-04 20:13 by babysonlkd
[考研] 一志愿东北大学085901土木专硕345求调剂 +3 zxt11111 2026-04-04 3/150 2026-04-04 14:21 by 土木硕士招生
[考研] 278求调剂 +6 Yy7400 2026-04-03 6/300 2026-04-04 09:53 by zhangdingwa
[考研] 英一数一408,总分284,二战真诚求调剂 +13 12.27 2026-03-30 15/750 2026-04-03 14:41 by 氮气气气
[考研] 求调剂 +3 心想事成可 2026-04-03 3/150 2026-04-03 11:22 by wangjy2002
[考研] 一志愿中国科学院大学265求调剂 +9 恬淡ye 2026-03-31 10/500 2026-04-03 11:10 by txp1986
[考研] 一志愿厦门大学材料工程专硕354找调剂!!! +8 贝呗钡钡 2026-03-30 8/400 2026-04-03 09:41 by hypershenger
[考研] 318求调剂 +3 笃行致远. 2026-03-31 4/200 2026-04-02 15:56 by Jaylen.
[考研] 一志愿北交大材料工程,总分358 +4 cs0106 2026-04-01 4/200 2026-04-02 07:42 by 尚水阁主
[考研] 284求调剂 +12 小熊~~ 2026-03-31 12/600 2026-04-01 20:23 by 花??
[考研] 材料专业调剂 +5 啦啦啦哭 2026-03-31 6/300 2026-04-01 16:48 by JourneyLucky
[考研] 08工科275求调剂,可跨考。 +5 AaAa7420 2026-03-31 5/250 2026-04-01 15:21 by 159357hjz
[考研] 311求调剂 +10 李芷新1 2026-03-31 10/500 2026-04-01 14:38 by chenqifeng666
[考研] 08工科,295,接受跨专业调剂 +6 lmnlzy 2026-03-31 6/300 2026-04-01 11:02 by 逆水乘风
[考研] 一志愿中海洋320化学工程与技术学硕求调剂 +8 披星河 2026-03-30 8/400 2026-03-31 08:53 by lbsjt
信息提示
请填处理意见