| 查看: 1739 | 回复: 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) 请问一下，我应该怎么做才是合理的？谢谢！ |
回复此楼» 猜你喜欢
酰胺脱乙酰基
已经有9人回复
-
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有3人回复
-
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优
已经有5人回复
-
有70后还继续奋斗在职场上的吗？
已经有6人回复
-
博士延得我，科研能力直往上蹿
已经有7人回复
-
退学或坚持读
已经有27人回复
-
面上基金申报没有其他的参与者成吗
已经有5人回复
-
遇见不省心的家人很难过
已经有22人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
» 抢金币啦！回帖就可以得到:
-
南方医科大学中药学院 申请考核博士一名 （天然药化方向，天然产物分离经验优先）
+1/270
-
湖南大学材料院陶益杰老师招收2026年秋季入学博士生一名及联合培养硕士学生一名
+1/195
-
西南科技大学曹克课题组招收2026级申请考核制有机化学博士研究生
+1/188
-
香港岭南大学郭瑛课题组博士招生 (2026年秋季入学）
+1/180
-
多功能 电子微生物生长分析仪 及 微生物快检技术开发服务
+2/140
-
诚聘 有机光电材料计算方向 博士后、科研助理、访问学生
+1/75
-
结构动力学与结构健康监测方向欧盟玛丽居里全奖博士招聘
+1/64
-
山东征女友，坐标济南
+1/64
-
深圳大学信息功能电子材料方向“申请-考核制”博士生招生
+2/44
-
新加坡 南洋理工大学- 智能光子/ 传感 PHD 全奖一名 2026 - 8 月入学
+1/24
-
南京医科大学国家级高层次青年人才团队招收博士研究生
+1/11
-
哈工大 张乃庆课题组招收博士快响计划（名额充足），通过后随时入学
+1/10
-
江汉大学轩亮教授课题组招博士研究生/博士后
+1/10
-
中国科学院深圳先进技术研究院——招聘客座研究生
+3/7
-
中国科学院深圳先进技术研究院——招聘博士后
+3/7
-
江汉大学轩亮教授课题组招博士研究生/博士后
+1/7
-
【博士招生】哈工大（深圳）智能学部机器人与先进制造学院 陆文杰老师课题组
+1/4
-
山东第一医科大学第一附属医院招聘事业编制科研岗
+1/2
-
武汉工程大学董志兵教授课题组招收博士/硕士研究生（长期有效）
+1/1
-
武汉理工大学数学与统计学院张秀军教授课题组招收2026级秋季博士研究生
+1/1
3楼2011-04-18 19:10:46
★ ★ ★ ★
小木虫(金币+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
4楼2011-04-19 09:07:48
5楼2011-04-19 10:53:16