| 查看: 2527 | 回复: 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 |
回复此楼» 猜你喜欢
求调剂0703
已经有5人回复
-
一志愿天津大学339材料与化工求调剂
已经有4人回复
-
一志愿华南师范361分，化学求调剂
已经有3人回复
-
282求调剂
已经有6人回复
-
一志愿 南京航空航天大学 ，080500材料科学与工程学硕
已经有10人回复
-
273求调剂
已经有13人回复
-
一志愿085600中科院宁波所276分求调剂
已经有21人回复
-
求调剂
已经有6人回复
-
339求调剂
已经有5人回复
-
265求调剂
已经有6人回复
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以内。（报错写的很清楚） 楼上说的方法可能可行。 若不行，还是仔细要看看你对接的结果 |
aizazadi5253楼2015-11-09 02:10:58
2楼2015-11-06 21:11:52
aizazadi525
超级版主
- 应助: 2 (幼儿园)
- 金币: 5372.7
- 散金: 5
- 红花: 7
- 帖子: 625
- 在线: 377.2小时
- 虫号: 1146537
- 注册: 2010-11-14
- 性别: GG
- 专业: 高分子合成化学
4楼2015-11-09 15:32:34
aizazadi525
主管区长
- 应助: 2 (幼儿园)
- 金币: 5372.7
- 散金: 5
- 红花: 7
- 帖子: 625
- 在线: 377.2小时
- 虫号: 1146537
- 注册: 2010-11-14
- 性别: GG
- 专业: 高分子合成化学
5楼2015-11-09 17:02:35
