| 查看: 2373 | 回复: 14 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
aizazadi525铁杆木虫 (正式写手)
|
[求助]
能量优化之后还是大于1000 已有2人参与
|
|
|
小虫刚接触gromacs，现在在对一个自己对接的蛋白复合物进行动力学模拟，在能量优化过程中，使用如下.mdp文件 ; Title of run ; Parameters describing what to do, when to stop and what to save integrator = steep ; Algorithm (steep = steepest descent minimization) emtol = 1000.0 ; Stop minimization when the maximum force < 10.0 kJ/mol emstep = 0.01 ; Energy step size nsteps = 50000 ; Maximum number of (minimization) steps to perform energygrps = system ; Which energy group(s) to write to disk ; Parameters describing how to find the neighbors of each atom and how to calculate the interactions nstlist = 1 ; Frequency to update the neighbor list and long range forces ns_type = grid ; Method to determine neighbor list (simple, grid) rlist = 1.0 ; Cut-off for making neighbor list (short range forces) coulombtype = PME ; Treatment of long range electrostatic interactions rcoulomb = 1.0 ; long range electrostatic cut-off rvdw = 1.0 ; long range Van der Waals cut-off pbc = xyz ; Periodic Boundary Conditions (yes/no) 但是运行总是很快就结束了，返回结果说： Stepsize too small, or no change in energy. Converged to machine precision, but not to the requested precision Fmax< 1000 试着调了几个参数，如把nsteps调整为100000，或把emstep调整为0.02或0.03，但还是出现同样的结果。想着不报错就继续运行，但是在NVT平衡的时候，就会出现LINCE WARNING,出现太多程序就自动停止了。 有没有哪位大神知道一下，这种情况该如何处理。谢谢了！！  2015-11-06_172833.png  2015-11-06_172803.png  2015-11-06_172600.png |
回复此楼» 猜你喜欢
全日制（定向）博士
已经有5人回复
-
假如你的研究生提出不合理要求
已经有10人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
aizazadi525
铁杆木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 5372.7
- 散金: 5
- 红花: 7
- 帖子: 625
- 在线: 377.2小时
- 虫号: 1146537
- 注册: 2010-11-14
- 性别: GG
- 专业: 高分子合成化学
送红花一朵 4楼2015-11-09 15:32:34
angjzy
木虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 1745.2
- 散金: 1491
- 红花: 5
- 帖子: 251
- 在线: 172.3小时
- 虫号: 647898
- 注册: 2008-11-06
- 专业: 药物材料
2楼2015-11-06 21:11:52
awaken2013
铜虫 (小有名气)
- 模拟EPI: 2
- 应助: 26 (小学生)
- 金币: 137.9
- 红花: 15
- 帖子: 243
- 在线: 149.6小时
- 虫号: 2854550
- 注册: 2013-12-07
- 性别: GG
- 专业: 生物大分子结构与功能
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
aizazadi525: 金币+10, ★★★很有帮助 2015-11-09 11:07:28
感谢参与，应助指数 +1
aizazadi525: 金币+10, ★★★很有帮助 2015-11-09 11:07:28
|
只能说你自己对接的这个构象，无法被能量最小化到1000以内。（报错写的很清楚） 楼上说的方法可能可行。 若不行，还是仔细要看看你对接的结果 |
» 本帖已获得的红花（最新10朵）
3楼2015-11-09 02:10:58
aizazadi525
铁杆木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 5372.7
- 散金: 5
- 红花: 7
- 帖子: 625
- 在线: 377.2小时
- 虫号: 1146537
- 注册: 2010-11-14
- 性别: GG
- 专业: 高分子合成化学
5楼2015-11-09 17:02:35