24小时热门版块排行榜    

查看: 659  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ddwy

木虫 (小有名气)

[求助] 分子动力学模拟出错 已有1人参与

大家好:
我用DL_POLY_4.03做分子动力学模拟,CONTROL输入文件是:

temperature 298.15
# pressure 0.001
ensemble nvt hoover  0.5

integrat velocity
restart
steps 3000000
equilibration 1000000
mult  1
scale  5
stats 20
rdf            10
print rdf

traj 1000001 100 0

densvar 20

timestep 0.002
cutoff        12.0

ewald precision 1.0E-6

minimise energy 50 5.0E-03

shake 1.0E-6
mxshak  10000
quaternion tolerance  1.0E-05

job time       99999999999999.0
close time     200000.0

finish
但是计算到第步的时候出现错误,OUTPUT文件为:
----------------------------------------------------------------------------------------------------------------------------------
   174600 -2.8822E+04  2.9815E+02 -3.1606E+04  3.9295E+03 -3.5876E+04  1.7169E+02  1.0398E+02  6.5269E+01  0.0000E+00
349.20000 -2.9405E+04  0.0000E+00  9.0503E+03 -9.0711E+04  3.5875E+04  8.3469E+01  0.0000E+00  6.3803E+04  0.0000E+00
119751.3  4.9579E+04  0.0000E+00  0.0000E+00  0.0000E+00  9.0000E+01  9.0000E+01  9.0000E+01  0.0000E+00 -1.2860E+00

  rolling -2.6814E+04  2.9487E+02 -2.9563E+04  4.1468E+03 -3.4341E+04  3.6908E+02  1.8861E+02  7.4187E+01  0.0000E+00
averages -2.6525E+04  0.0000E+00  6.3658E+03 -9.3434E+04  3.4340E+04  1.4041E+03  0.0000E+00  6.4056E+04  0.0000E+00
           4.9579E+04  0.0000E+00  0.0000E+00  0.0000E+00  9.0000E+01  9.0000E+01  9.0000E+01  0.0000E+00 -8.0166E-02
----------------------------------------------------------------------------------------------------------------------------------

*** warning - global constraint number         6 , with particle numbers:      2885 &      2886 ,
converges to a length of     7.94 Angstroms with a factor  3.1039E+01 ,contributes towards next error !!! ***

*** warning - global constraint number         7 , with particle numbers:      2885 &      2887 ,
converges to a length of    12.88 Angstroms with a factor  8.2412E+01 ,contributes towards next error !!! ***

*** warning - global constraint number         8 , with particle numbers:      2886 &      2887 ,
converges to a length of    18.15 Angstroms with a factor  1.0006E+02 ,contributes towards next error !!! ***

DL_POLY_4 terminated due to error   105

error - shake algorithm (constraints_shake) failed to converge
"OUTPUT" 64939L, 6075639C                     


请各位指点一下,谢谢
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

我喜欢收集微笑没有原因……
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddwy

木虫 (小有名气)

引用回帖:
2楼: Originally posted by 月只蓝 at 2014-09-03 21:15:05
把该问题发到分子模拟版块去吧。

谢谢
我喜欢收集微笑没有原因……
3楼2014-09-04 08:50:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
把该问题发到分子模拟版块去吧。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2014-09-03 21:15:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-09-02 00:19 by xiangy672
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
信息提示
请填处理意见