| 查看: 3468 | 回复: 7 | ||
xueer小雯金虫 (初入文坛)
|
[求助]
求助：Gromacs平衡时蛋白跑出盒子了，怎么办？已有2人参与
|
|
图片3是进行md之前的蛋白与盒子状态，即npt结束后。 图片2是md 20ns后的状态，蛋白几乎全都出去了，盒子中的是小分子。 md.mdp 如下，请高人帮我看看怎么回事，我之前从来没遇到过这样的问题。。。 title = Protein-ligand complex NVT equilibration ; Run parameters integrator = md ; leap-frog integrator nsteps = 10000000 ; 2 * 10000000 = 20000 ps (20 ns) dt = 0.002 ; 2 fs ; Output control nstxout = 1000 ; suppress .trr output nstvout = 1000 ; suppress .trr output nstenergy = 1000 ; save energies every 2 ps nstlog = 1000 ; update log file every 2 ps nstxtcout = 1000 ; write .xtc trajectory every 2 ps energygrps = Protein NAD ; Bond parameters continuation = yes ; first dynamics run constraint_algorithm = lincs ; holonomic constraints constraints = all-bonds ; all bonds (even heavy atom-H bonds) constrained lincs_iter = 1 ; accuracy of LINCS lincs_order = 4 ; also related to accuracy ; Neighborsearching ns_type = grid ; search neighboring grid cells nstlist = 5 ; 10 fs rlist = 0.9 ; short-range neighborlist cutoff (in nm) rcoulomb = 0.9 ; short-range electrostatic cutoff (in nm) rvdw = 1.4 ; short-range van der Waals cutoff (in nm) ; Electrostatics coulombtype = PME ; Particle Mesh Ewald for long-range electrostatics pme_order = 4 ; cubic interpolation fourierspacing = 0.16 ; grid spacing for FFT ; Temperature coupling tcoupl = V-rescale ; modified Berendsen thermostat tc-grps = Protein_NAD Water_and_ions ; two coupling groups - more accurate tau_t = 0.1 0.1 ; time constant, in ps ref_t = 298 298 ; reference temperature, one for each group, in K ; Pressure coupling pcoupl = Parrinello-Rahman ; pressure coupling is on for NPT pcoupltype = isotropic ; uniform scaling of box vectors tau_p = 2.0 ; time constant, in ps ref_p = 1.0 ; reference pressure, in bar compressibility = 4.5e-5 ; isothermal compressibility of water, bar^-1 ; Periodic boundary conditions pbc = xyz ; 3-D PBC ; Dispersion correction DispCorr = EnerPres ; account for cut-off vdW scheme ; Velocity generation gen_vel = no ; assign velocities from Maxwell distribution 多谢啊！  2.PNG  3.PNG |
回复此楼» 猜你喜欢
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有3人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
跑完Gromacs的动力学模拟之后的轨迹文件.trr里面没有内容，为什么？
已经有4人回复- gromacs安装出错求助
已经有4人回复
- 【求助】这是我安装gromacs过程，但是出错阿，这个怎么办呢？谢谢
已经有11人回复
wangyan10
铁杆木虫 (正式写手)
- 应助: 36 (小学生)
- 金币: 8272.5
- 散金: 58
- 红花: 4
- 帖子: 368
- 在线: 254.3小时
- 虫号: 1089510
- 注册: 2010-09-04
- 性别: GG
- 专业: 理论和计算化学
2楼2014-02-26 16:53:55
xueer小雯
金虫 (初入文坛)
- 应助: 3 (幼儿园)
- 金币: 2240.9
- 红花: 1
- 帖子: 36
- 在线: 359.1小时
- 虫号: 2318904
- 注册: 2013-03-05
- 专业: 药物设计与药物信息
3楼2014-02-26 19:24:49
xueer小雯
金虫 (初入文坛)
- 应助: 3 (幼儿园)
- 金币: 2240.9
- 红花: 1
- 帖子: 36
- 在线: 359.1小时
- 虫号: 2318904
- 注册: 2013-03-05
- 专业: 药物设计与药物信息
4楼2014-03-03 15:58:53
5楼2014-03-03 17:18:55
xueer小雯
金虫 (初入文坛)
- 应助: 3 (幼儿园)
- 金币: 2240.9
- 红花: 1
- 帖子: 36
- 在线: 359.1小时
- 虫号: 2318904
- 注册: 2013-03-05
- 专业: 药物设计与药物信息
6楼2014-03-03 17:21:53
7楼2014-08-19 09:07:14
|
本帖内容被屏蔽 |
8楼2021-07-08 18:13:01