24小时热门版块排行榜    

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

转基因猴子

木虫 (正式写手)

引用回帖:
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的回帖

转基因猴子

木虫 (正式写手)

引用回帖:
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的回帖

转基因猴子

木虫 (正式写手)

引用回帖:
11楼: Originally posted by xueer小雯 at 2014-04-03 19:17:29
nonono
相反是得加这个参数,才能保证你的体系是在一个无限大的空间内重复,而你的模拟只是取了其中的一个而已。
蛋白没有跑出盒子的意思是它只是有一部分跑到了相邻的盒子里,这是不影响模拟的,因为每一个盒子 ...

既然不影响模拟结果为何还要用trjconv来纠正呢?

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

转基因猴子

木虫 (正式写手)

送红花一朵
引用回帖:
13楼: Originally posted by xueer小雯 at 2014-04-04 10:10:16
因为蛋白可能会跳在两个盒子之间,修正轨迹后才能在一个盒子里计算,修正是为了计算方便。否则例如跨盒子的两个原子之间的距离它可能告诉你的就变成了该盒子中这两个原子的距离,这时候就有了一个盒子的距离差,你 ...

我明白了,非常感谢你的解答

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

转基因猴子

木虫 (正式写手)

引用回帖:
15楼: Originally posted by QQ_Q_Q_QQ at 2014-04-04 14:48:02
comm_mode              = Linear
comm_grps              = Protein Non-protein

这两个参数指征什么意思?
我还年轻,我要上路!
16楼2014-04-04 15:38:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

转基因猴子

木虫 (正式写手)

引用回帖:
18楼: Originally posted by huyingdamon at 2014-06-27 13:18:05
您好 看到你的帖子 我想问一下我模拟环糊精跟另一个小分子,NVT后发现这两个分子都有某些时候跑出了盒子,看到你对别人的回复说用trjconv 去周期就好了,我对gromacs还不深,所以没看懂,是NVT模拟之后再去除周期, ...

在mdp中加上pdb项
我还年轻,我要上路!
19楼2014-06-27 13:30:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 转基因猴子 的主题更新
信息提示
请填处理意见