24小时热门版块排行榜    

查看: 4545  |  回复: 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的回帖

转基因猴子

木虫 (正式写手)

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

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

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

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 268求调剂 +7 一定有学上- 2026-03-14 8/400 2026-03-17 13:10 by zz820
[考研] 332求调剂 +5 Zz版 2026-03-13 5/250 2026-03-17 11:01 by JineShine
[考研] 290求调剂 +6 @将就将就看 2026-03-10 11/550 2026-03-17 10:22 by @将就将就看
[考研] 材料专硕274一志愿陕西师范大学求调剂 +5 薛云鹏 2026-03-13 5/250 2026-03-17 10:15 by Sammy2
[考研] 304求调剂 +7 小熊joy 2026-03-14 7/350 2026-03-17 08:53 by 雾散后相遇lc
[考研] 341求调剂 +5 捣蛋猪猪 2026-03-11 5/250 2026-03-17 06:56 by Muziwl
[基金申请] 国自科面上基金字体 +6 iwuli 2026-03-12 7/350 2026-03-16 21:18 by sculhf
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 301求调剂 +6 Liyouyumairs 2026-03-11 6/300 2026-03-13 20:11 by JourneyLucky
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
[考研] 298求调剂 +3 Vv呀! 2026-03-10 3/150 2026-03-10 22:40 by 剑诗杜康
信息提示
请填处理意见