| 查看: 1784 | 回复: 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) 请问一下，我应该怎么做才是合理的？谢谢！ |
回复此楼» 猜你喜欢
调剂310
已经有11人回复
-
300求调剂，材料科学英一数二
已经有9人回复
-
296求调剂
已经有3人回复
-
一志愿南昌大学324求调剂
已经有3人回复
-
299求调剂
已经有9人回复
-
22408 359分调剂
已经有6人回复
-
085600 材料与化工 329分求调剂
已经有12人回复
-
一志愿北京化工大学材料与化工（085600）296求调剂
已经有14人回复
-
085601材料工程找调剂
已经有11人回复
-
总分293求调剂
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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/286
-
数一英一硕士调剂（双一流高校）
+1/183
-
五邑大学药学与食品工程学院 | 刘宇茸课题组 调剂招生
+5/95
-
赣南师范大学肖兴教授课题组招2026级电子信息硕士 4-6 名，欢迎调剂
+1/82
-
太原科技大学材料科学与工程学院招收金属材料类博士一名
+1/82
-
青岛大学固态电池研究团队招收材料、化学、物理等方向博士研究生
+2/36
-
【林业专业硕士调剂1名】中国林业科学研究院
+1/34
-
福建师范大学环境与资源学院环境电化学课题组招收2026硕士（含调剂）
+1/29
-
三峡大学材化学院李东升校长团队接收化学类硕士调剂
+1/21
-
年薪约32-60万！成会明院士团队诚聘材料理论计算博士后（长期有效）
+1/15
-
中国科学院上海硅酸盐研究所施剑林院士团队胡萍课题组招聘博士后
+1/14
-
双一流大学湘潭大学“化工过程模拟与强化”国家地方联合工程研究中心招收博士生
+1/14
-
西安工程大学能源化工新材料研发项目组接收2026年计算机类调剂研究生
+1/13
-
烟台大学先进陶瓷方向 招收材料专业调剂研究生
+1/11
-
08工学英一数二290 三个奖学金两个省三 一个省创优秀结项 B区求调剂
+1/10
-
欢迎各位同学申请参加澳门大学健康科学学院的2026年夏令营，报名截止4月12号
+1/4
-
武汉纺织大学全国重点实验室吕永钢教授课题组招收材料、化学、生物相关专业硕士研究生
+1/4
-
天津医科大学基础医学院张恒课题组博士后招聘
+1/3
-
26博士申请
+1/2
-
上海师范大学有机化学专业研究生招生
+1/1
★ ★ ★ ★
小木虫(金币+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
4楼2011-04-19 09:07:48
5楼2011-04-19 10:53:16
6楼2011-04-19 13:30:13
