24小时热门版块排行榜    

Znn3bq.jpeg
查看: 563  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 307中医考研调剂 +9 于以采蘩 2026-04-14 9/450 2026-04-19 08:41 by 烟雨流涯
[考研] 327求调剂 +27 Xxjc1107. 2026-04-13 30/1500 2026-04-19 08:22 by cuisz
[考研] 291求调剂 +10 关忆北. 2026-04-14 10/500 2026-04-18 23:32 by 路病情
[考研] 一志愿沪9,326求生物学调剂 +12 刘墨墨 2026-04-13 12/600 2026-04-18 23:31 by 路病情
[考研] 320求调剂 +5 深郊akm 2026-04-17 5/250 2026-04-18 19:52 by 王珺璞
[考研] 化学070300 求调剂 +29 哈哈哈^_^ 2026-04-12 29/1450 2026-04-18 15:56 by Equinoxhua
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[考研] 一志愿中科大材料与化工,353分还有调剂学校吗 +10 否极泰来2026 2026-04-15 12/600 2026-04-17 17:54 by mapenggao
[考研] 294求调剂 +7 淡然654321 2026-04-17 8/400 2026-04-17 16:36 by wutongshun
[考博] 申博自荐 +3 Linxia林夏 2026-04-13 3/150 2026-04-16 12:55 by 墨荷之露
[考研] 药学求调剂 +14 喽哈加油 2026-04-14 16/800 2026-04-16 10:15 by beilsong20
[考研] 求调剂推荐 +8 小聂爱学习 2026-04-14 8/400 2026-04-16 07:22 by 学员JpLReM
[考研] 求调剂学校 +14 不会吃肉 2026-04-13 16/800 2026-04-15 21:59 by noqvsozv
[考研] 0854调剂 +13 长弓傲 2026-04-12 16/800 2026-04-15 13:45 by fenglj492
[考研] 085801电气专硕272求调剂 +19 电气李 2026-04-13 21/1050 2026-04-15 13:37 by 黑科技矿业
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[考研] 各位老师好,求调剂,本科211,一志愿天津大学生物与医药学硕,差两名录取。 +11 路六六jjj 2026-04-13 11/550 2026-04-14 16:01 by zs92450
[考研] 085408光电信息工程专硕355一志愿长春光机所调剂 +6 王ymaa 2026-04-13 13/650 2026-04-14 11:33 by 王ymaa
[考研] 302求调剂 +10 易!? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
信息提示
请填处理意见