| 查看: 1921 | 回复: 5 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
mdrun时出错
|
|||
|
mdrun -s aa.tpr -o aa.trr -c aa.gro & 出现了如下的提示错误： step 11: Water molecule starting at atom 180665 can not be settled. Check for bad contacts and/or reduce the timestep if appropriate. Wrote pdb files with previous and current coordinates step 13: Water molecule starting at atom 180665 can not be settled. Check for bad contacts and/or reduce the timestep if appropriate. Wrote pdb files with previous and current coordinates Stepsize too small, or no change in energy. Converged to machine precision, but not to the requested precision Fmax < 1000 Double precision normally gives you higher accuracy. You might need to increase your constraint accuracy, or turn off constraints alltogether (set constraints = none in mdp file) writing lowest energy coordinates. Steepest Descents converged to machine precision in 31 steps, but did not reach the requested Fmax < 1000. Potential Energy = -3.1653622e+06 Maximum force = 9.5183238e+05 on atom 180665 Norm of force = 2.3416899e+03 下面是我的em.mdp文件 em.mdp er spoel (236) ; Wed Nov 3 17:12:44 1993 ; Input file ; cpp = /lib/cpp constraints = none define = -DFLEX_SPC ; Parameters describing what to do, when to stop and what to save integrator = steep dt = 0.002 emtol = 1000.0 emstep = 0.01 nsteps = 3000 ; Parameters describing how to find the neighbors of each atom and how to calculate the interactions nstlist = 10 ns_type = grid rlist = 1.0 coulombtype = pme rcoulomb = 1.0 rvdw = 1.0 fourierspacing =0.12 pme_order =4 optimize_fft =yes freezegrps =Protein UNK POP freezedim =Y Y Y Y Y Y Y Y Y 我也尝试了改变em.mdp文件中的一些参数(只对一个参数做了修改)： 1.em.mdp文件中的freezegrps改为Protein UNK能量能很好地收敛，下面是它的显示结果： Steepest Descents: Tolerance (Fmax) = 1.00000e+03 Number of steps = 3000 writing lowest energy coordinates. Steepest Descents converged to Fmax < 1000 in 1678 steps Potential Energy = -4.2245720e+06 Maximum force = 9.5952905e+02 on atom 9032 Norm of force = 1.1612597e+01 2.em.mdp文件中的emstep改为0.00001（参照别人的回答), 结果显示为： Steepest Descents: Tolerance (Fmax) = 1.00000e+03 Number of steps = 3000 writing lowest energy coordinates. Steepest Descents converged to Fmax < 1000 in 32 steps Potential Energy = -3.0919205e+06 Maximum force = 9.1417908e+02 on atom 122811 Norm of force = 1.1197423e+02 这样的结果应该是有问题的吧，只跑了32步。 3.integrator改为md, 结果为： Writing final coordinates. Average load imbalance: 0.2 % Part of the total run time spent waiting due to load imbalance: 0.1 % Parallel run - timing based on wallclock. NODE (s) Real (s) (%) Time: 7515.940 7515.940 100.0 2h05:15 (Mnbf/s) (GFlops) (ns/day) (hour/ns) Performance: 63.678 3.946 0.115 208.755 gcq#290: "Oh, There Goes Gravity" (Eminem) 请问一下，我应该怎么做才是合理的？谢谢！ |
回复此楼» 猜你喜欢
国自然评审意见
已经有8人回复
-
国自然面上复盘~欢迎讨论
已经有14人回复
-
系统查不到
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有12人回复
-
有没有仍没收到信息的
已经有3人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有8人回复
-
基金不中，共勉
已经有11人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有9人回复
-
2026年叶企孙基金
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Gromacs中，mdrun出错
已经有5人回复
- 在gromacs中运行mdrun命令时出错
已经有10人回复
- vsap 5.2 环境和编译问题
已经有14人回复
- Gromacs做steep之后做cg出错
已经有5人回复
- 编译vasp出错
已经有11人回复
- gromacs做全原子模拟时候，突然到MD这一步出错
已经有6人回复
- 【求助】安装vasp出错 make: *** [fftmpi_map.o] 错误 1【已解决】
已经有6人回复
- 【求助】vasp安装编译错误，请大家帮忙看下那里出错了。
已经有6人回复
- 【求助成功】Vasp 5.2.11编译出错
已经有10人回复
- 【求助】在编译时出错
已经有6人回复
- 【求助】PBE0 计算Si、Ge能带结构
已经有9人回复
» 抢金币啦！回帖就可以得到:
-
大连理工大学王铁胜副教授招收2027年入学博士生
+1/168
-
2026年工程管理与智慧城市国际会议 (EMSC 2026)
+1/72
-
南京理工大学尚露教授团队联合中科院理化所张铁锐研究员团队招聘催化领域博士后
+1/34
-
中国海洋大学诚聘“人工智能（AI for Science）”方向青年英才
+1/17
-
山东大学-山东字节信息智能技术实验室2026年招聘博士后
+1/16
-
中国海洋大学诚聘“化学 + 人工智能（AI for Science）”方向青年英才
+1/15
-
中国海洋大学诚聘“化学 + 人工智能（AI for Science）”方向青年英才
+1/15
-
三乙基碳酸氢铵（TEAB）缓冲液的去除
+1/12
-
上海交通大学陆军课题组招聘博士后 - 计算模拟方向
+1/6
-
寻求合作【大学教师】兼职出版顾问｜无需坐班｜一单5000元
+1/5
-
北京大学袁琳洁课题组诚招药理及化学生物学方向27年入学博士生
+1/4
-
西安石油大学 发明专利 一种改性磷酸钾镁水泥及其应用 转让
+1/4
-
2027材料申博自荐C9本国科大硕
+1/4
-
北京大学招聘有机化学/药物化学/化学生物学方向博士后
+1/4
-
寻求合作【出版代理】兼职推广｜无需坐班｜一单5000元
+1/2
-
寻求合作【大学教师】兼职出版顾问｜无需坐班｜一单5000元
+1/2
-
项目助理招聘启事
+1/2
-
专利审查意见答复，是AI辅助生成中很难的一个功能...
+1/2
-
浙江大学医学院黄浩杰教授团队招聘4-5名博士后
+1/1
-
TGF‑β家族分子生物学机制与样本前处理检测关键技术
+1/1
4楼2011-04-19 09:07:48
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
Jfreda(金币+5): 2011-04-18 19:08:54
zh1987hs(金币+3): 谢谢 2011-04-18 21:33:02
小木虫(金币+0.5):给个红包，谢谢回帖
Jfreda(金币+5): 2011-04-18 19:08:54
zh1987hs(金币+3): 谢谢 2011-04-18 21:33:02
|
Steepest Descents converged to Fmax < 1000 in 32 steps Potential Energy = -3.0919205e+06 Maximum force = 9.1417908e+02 on atom 122811 Norm of force = 1.1197423e+02 这样的结果应该是有问题的吧，只跑了32步。 虽然只跑了32步，不过你的力已经基本收敛了，应该没有什么问题了。应该不用freegroup,如果要限制原子的话，可以尝试使用define. |
2楼2011-04-18 18:42:24
3楼2011-04-18 19:10:46
5楼2011-04-19 10:53:16
