24小时热门版块排行榜    

查看: 545  |  回复: 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的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
把该问题发到分子模拟版块去吧。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2014-09-03 21:15:05
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 ddwy 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿东华大学化学070300,求调剂 +7 2117205181 2026-03-21 8/400 2026-03-22 22:55 by chixmc
[考研] 石河子大学(211、双一流)硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 308求调剂 +3 墨墨漠 2026-03-21 3/150 2026-03-22 16:54 by i_cooler
[考研] 一志愿中南化学(0703)总分337求调剂 +9 niko- 2026-03-19 10/500 2026-03-22 16:08 by ColorlessPI
[考研] 286求调剂 +10 Faune 2026-03-21 10/500 2026-03-21 23:34 by 314126402
[考研] 材料工程专硕 348分求调剂 +3 冬辞. 2026-03-17 5/250 2026-03-21 18:47 by 学员8dgXkO
[考研] 工科0856求调剂 +3 沐析汀汀 2026-03-21 3/150 2026-03-21 18:30 by 学员8dgXkO
[考研] 求助 +5 梦里的无言 2026-03-21 6/300 2026-03-21 17:51 by 学员8dgXkO
[考研] 297求调剂 +11 戏精丹丹丹 2026-03-17 12/600 2026-03-21 17:47 by ColorlessPI
[考研] 311求调剂 +3 勇敢的小吴 2026-03-20 3/150 2026-03-21 17:40 by ColorlessPI
[考研] 材料 271求调剂 +5 展信悦_ 2026-03-21 5/250 2026-03-21 17:29 by 学员8dgXkO
[考研] 303求调剂 +5 睿08 2026-03-17 7/350 2026-03-21 03:11 by JourneyLucky
[考研] 一志愿中国石油大学(华东) 本科齐鲁工业大学 +3 石能伟 2026-03-17 3/150 2026-03-21 02:22 by JourneyLucky
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 271材料工程求调剂 +8 .6lL 2026-03-18 8/400 2026-03-21 00:58 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
信息提示
请填处理意见