24小时热门版块排行榜    

查看: 4566  |  回复: 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
回复此楼

» 猜你喜欢

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

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

QQ_Q_Q_QQ

新虫 (正式写手)

引用回帖:
16楼: Originally posted by 转基因猴子 at 2014-04-04 15:38:26
这两个参数指征什么意思?...

看手册
17楼2014-04-04 16:35:00
已阅   回复此楼   关注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的回帖

转基因猴子

木虫 (正式写手)

引用回帖:
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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +6 墨墨漠 2026-03-20 6/300 2026-03-21 09:23 by 天青等烟雨
[考研] 346求调剂[0856] +4 WayneLim327 2026-03-16 7/350 2026-03-21 04:02 by JourneyLucky
[考研] 机械专硕299求调剂至材料 +3 kkcoco25 2026-03-16 4/200 2026-03-21 03:52 by JourneyLucky
[考研] 301求调剂 +10 yy要上岸呀 2026-03-17 10/500 2026-03-21 03:14 by JourneyLucky
[考研] 085600材料与化工 +8 安全上岸! 2026-03-16 8/400 2026-03-20 22:13 by luoyongfeng
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 5/250 2026-03-20 22:11 by 云游重阳
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 生物学调剂招人!!! +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 6/300 2026-03-19 17:18 by fei626-918
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 304求调剂 +12 小熊joy 2026-03-14 13/650 2026-03-18 12:34 by Linda Hu
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
信息提示
请填处理意见