24小时热门版块排行榜    

查看: 4539  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 转基因猴子 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 326求调剂 +4 上岸的小葡 2026-03-15 5/250 2026-03-16 08:39 by Linda Hu
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 材料专硕326求调剂 +4 墨煜姒莘 2026-03-15 4/200 2026-03-15 11:02 by dyw
[考研] 309求调剂 +4 花与叶@ 2026-03-10 4/200 2026-03-14 21:26 by a不易
[考研] 一志愿天大化工(085600)调剂总分338 +6 蔡大美女 2026-03-09 6/300 2026-03-14 02:46 by JourneyLucky
[考研] 085600求调剂 +3 a邵星池 2026-03-09 3/150 2026-03-14 01:32 by JourneyLucky
[考研] 一志愿华中农业大学071010,总分三百二,求调剂 +3 困困困困坤坤 2026-03-10 3/150 2026-03-14 00:35 by JourneyLucky
[考研] 285 求调剂 资源与环境 一志愿北京化工大学 +3 未名考生 2026-03-10 3/150 2026-03-13 23:04 by JourneyLucky
[考研] 332求调剂 +3 zjy101327 2026-03-11 6/300 2026-03-13 22:48 by JourneyLucky
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 求调剂(材料与化工327) +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 一志愿西南交大,材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-12 5/250 2026-03-13 10:56 by houyaoxu
[考研] 274求调剂0856材料化工 +12 z2839474511 2026-03-11 13/650 2026-03-13 10:39 by peike
[考研] 296求调剂 +3 大口吃饭 身体健 2026-03-13 3/150 2026-03-13 10:31 by 学员8dgXkO
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
[考研] 298求调剂 +3 Vv呀! 2026-03-10 3/150 2026-03-10 22:40 by 剑诗杜康
[考研] 294 英二数二物化 求调剂 +6 米饭团不好吃 2026-03-09 6/300 2026-03-09 23:55 by barlinike
信息提示
请填处理意见