24小时热门版块排行榜    

查看: 4438  |  回复: 16

sxjsn1

金虫 (小有名气)

veni vedi veci

[求助] vasp分子动力学报错 已有3人参与

最近开始学习做vasp分子动力学,用实验得到的坐标数据和别人提供的INCAR例子稍作修改后进行计算。

报错如下:

-----------------------------------------------------------------------------
|                                                                             |
|           W    W    AA    RRRRR   N    N  II  N    N   GGGG   !!!           |
|           W    W   A  A   R    R  NN   N  II  NN   N  G    G  !!!           |
|           W    W  A    A  R    R  N N  N  II  N N  N  G       !!!           |
|           W WW W  AAAAAA  RRRRR   N  N N  II  N  N N  G  GGG   !            |
|           WW  WW  A    A  R   R   N   NN  II  N   NN  G    G                |
|           W    W  A    A  R    R  N    N  II  N    N   GGGG   !!!           |
|                                                                             |
|      VASP found    297 degrees of freedom                                   |
|      the temperature will equal 2*E(kin)/ (degrees of freedom)              |
|      this differs from previous releases, where T was 2*E(kin)/(3 NIONS).   |
|      The new definition is more consistent                                  |
|                                                                             |
-----------------------------------------------------------------------------

POSCAR, INCAR and KPOINTS ok, starting setup
FFT: planning ...(           4 )
reading WAVECAR
prediction of wavefunctions initialized - no I/O
entering main loop
       N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1     0.764392063216E+04    0.76439E+04   -0.67877E+05   744   0.153E+03
RMM:   2     0.162604293910E+04   -0.60179E+04   -0.16132E+05   744   0.684E+02
RMM:   3     0.355026647917E+04    0.19242E+04   -0.99770E+03   744   0.477E+02
RMM:   4    -0.355257110459E+03   -0.39055E+04   -0.34173E+04  2074   0.368E+02
RMM:   5    -0.556354687631E+03   -0.20110E+03   -0.20557E+03  2181   0.685E+01    0.104E+02
RMM:   6    -0.807969295148E+03   -0.25161E+03   -0.40985E+03  1730   0.976E+01    0.105E+02
RMM:   7    -0.518920100251E+03    0.28905E+03   -0.64031E+02  1489   0.525E+01    0.353E+01
RMM:   8    -0.499452185825E+03    0.19468E+02   -0.83667E+01  1584   0.143E+01    0.153E+01
RMM:   9    -0.503489505756E+03   -0.40373E+01   -0.90584E+01  1512   0.131E+01    0.286E+01
RMM:  10    -0.496681907434E+03    0.68076E+01   -0.12014E+01  1574   0.532E+00    0.127E+01
RMM:  11    -0.494480696413E+03    0.22012E+01   -0.59178E+00  1503   0.397E+00    0.362E+00
RMM:  12    -0.494307140858E+03    0.17356E+00   -0.97898E-01  1544   0.131E+00    0.170E+00
RMM:  13    -0.494289008337E+03    0.18133E-01   -0.26473E-01  1628   0.610E-01    0.150E+00
RMM:  14    -0.494254078873E+03    0.34929E-01   -0.10911E-01  1551   0.494E-01    0.679E-01
RMM:  15    -0.494248439290E+03    0.56396E-02   -0.26568E-02  1536   0.246E-01    0.438E-01
RMM:  16    -0.494246774664E+03    0.16646E-02   -0.89946E-03  1457   0.138E-01    0.310E-01
RMM:  17    -0.494245830769E+03    0.94389E-03   -0.40946E-03  1375   0.905E-02    0.165E-01
RMM:  18    -0.494245656830E+03    0.17394E-03   -0.10657E-03  1378   0.426E-02    0.116E-01
RMM:  19    -0.494245540146E+03    0.11668E-03   -0.65529E-04  1240   0.436E-02    0.581E-02
RMM:  20    -0.494245488923E+03    0.51223E-04   -0.12036E-04   865   0.199E-02  
1 T=    NaN E=  NaN           F= -.49424549E+03 E0= -.49424406E+03  EK=  NaN     
   SP= 0.00E+00 SK=  NaN bond charge predicted
REAL_OPT: internal ERROR:         -48         -48         -48          32
REAL_OPT: internal ERROR:         -48         -48         -48          32
REAL_OPT: internal ERROR:         -48         -48         -48          32
REAL_OPT: internal ERROR:         -48         -48         -48          32
REAL_OPT: internal ERROR:         -48         -48         -48          32
REAL_OPT: internal ERROR:         -48         -48         -48          32
REAL_OPT: internal ERROR:         -48         -48         -48          32
REAL_OPT: internal ERROR:         -48         -48         -48          32
REAL_OPT: internal ERROR:         -48         -48         -48          32
REAL_OPT: internal ERROR:         -48         -48         -48          32
REAL_OPT: internal ERROR:         -48         -48         -48          32
REAL_OPT: internal ERROR:         -48         -48         -48          32
Job  /jobmgr/lsf8.0/8.0/linux2.6-glibc2.3-x86_64/bin/intelmpi_wrapper -np 12 /apps/soft/vasp/vasp.5.2/vasp.5.2/vasp
TID   HOST_NAME   COMMAND_LINE            STATUS            TERMINATION_TIME
===== ========== ================  =======================  ===================00000 c19b15     /apps/soft/vasp/  Done                     11/26/2013 16:48:30
00001 c19b15     /apps/soft/vasp/  Done                     11/26/2013 16:48:30
00002 c19b15     /apps/soft/vasp/  Done                     11/26/2013 16:48:30
00003 c19b15     /apps/soft/vasp/  Done                     11/26/2013 16:48:30
00004 c19b15     /apps/soft/vasp/  Done                     11/26/2013 16:48:30
00005 c19b15     /apps/soft/vasp/  Done                     11/26/2013 16:48:30
00006 c19b15     /apps/soft/vasp/  Done                     11/26/2013 16:48:30
00007 c19b15     /apps/soft/vasp/  Done                     11/26/2013 16:48:30
00008 c19b15     /apps/soft/vasp/  Done                     11/26/2013 16:48:30
00009 c19b15     /apps/soft/vasp/  Done                     11/26/2013 16:48:30
00010 c19b15     /apps/soft/vasp/  Done                     11/26/2013 16:48:30
00011 c19b15     /apps/soft/vasp/  Done                     11/26/2013 16:48:30
回复此楼

» 猜你喜欢

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

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

夜半酒醒人去后,更持红烛赏残花
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sxjsn1

金虫 (小有名气)

veni vedi veci

这是INCAR
Zr53 MD
# start parameters
   NWRITE = 2
   ISTART = 0
   ICHARG = 2
   ENCUT = 380


#Electroic  Relaxation
   PREC  = Medium
   NELM =  150
   NELMDL = 5
   EDIFF = 1.0e-4
   LREAL =  auto


#Ionic relaxation
  NSW  =  900
  IBRION  =  0
  ISIF      =   0
  ISYM    =   0
  ISMAR  =  1
  SIGMA  =  0.1
IALGO  =  48
  LDIAG  =  T
  VOSKOWN = 1
  NGX = 78
  NGY = 74
  NGZ = 78


# MD flags
  SMASS  =  2
  POTIM  =  2
  TEBEG  = 400
  TEEND  = 400

#Write flags
  LWAVE =    F
  LCHARG =  F
  LVTOT  =   F
  LELF   =    F
  LORBIT =  F


这是POSCAR部分数据
Cu64Zr36
   1.00000000000000
    12.0000000000000000    0.0000000000000000    0.0000000000000000
     0.0000000000000000   11.8000000000000007    0.0000000000000000
     0.0000000000000000    0.0000000000000000   12.5000000000000000
   Cu   Zr
  64  36
Direct
  0.0445316012624346  0.8767680807860145  0.6543323012345823
  0.0612530510210334  0.5975315534439858  0.2762958886825047
  0.0621291046130281  0.8904940274270459  0.4368290639755159
  0.0637359446944004  0.1744694999170520  0.8685319037139934
  0.0812315172936593  0.2980550432761078  0.5168033781711330
  0.1084071158034305  0.6662314568253791  0.5056049433440037
  0.1344183616036290  0.3962247685368325  0.3367521396043429
  0.1471929769311657  0.7068228518304380  0.9628619945375689
  0.1617505579197100  0.2573339277376760  0.1724685513823090
  0.1669119948037405  0.6910890948599757  0.1340729096298503
  0.1932835613752673  0.4156836164572724  0.6676752363990417
  0.2244594204021482  0.8978810500978298  0.1344598214099051
  0.2250764654001998  0.2997802501499507  0.9128368076884499
  0.2389981702790394  0.0355756381307657  0.5641253099168893
  0.2508580642317213  0.8239284177063908  0.7989103058139762
夜半酒醒人去后,更持红烛赏残花
2楼2013-11-26 17:29:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sxjsn1

金虫 (小有名气)

veni vedi veci

有哪位大能知道我哪儿错了吗?
我NSW=900,1000,3000;POTIM=3,5;SMASS=0,1,2,3都试过了

原子位置根本就不动
夜半酒醒人去后,更持红烛赏残花
3楼2013-11-26 17:31:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snail594

金虫 (小有名气)

建议先按手册后面的MD部分的参数设置,参数越少,错误也就越少。之后再按照自己需要一步步添加部分参数设置。
Inlearningwetrust
4楼2013-11-26 17:53:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sxjsn1

金虫 (小有名气)

veni vedi veci

好吧,简化INCAR后还是出现同样错误,简化后的INCAR如下
Zr53 MD
# start parameters
   ISTART = 0
   ICHARG = 2
   ENCUT = 380

#Electroic  Relaxation
   PREC  = Medium
   NELM =  150
   EDIFF = 1.0e-4
   LREAL =  auto


#Ionic relaxation
  NSW  =  900
  IBRION  =  0
  ISIF      =   0
  ISYM    =   0
  ISMAR  =  1
  SIGMA  =  0.1
  IALGO  =  48
  VOSKOWN = 1
  NGX = 78
  NGY = 74
  NGZ = 78


# MD flags
  SMASS  =  2
  POTIM  =  2
  TEBEG  = 400
  TEEND  = 400

#Write flags
  LWAVE =    F
  LCHARG =  F
夜半酒醒人去后,更持红烛赏残花
5楼2013-11-26 19:01:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sxjsn1

金虫 (小有名气)

veni vedi veci

现在报错  REAL_OPT: internal ERROR:          -8         -16          -8          52

INCAR如下,应该不能再简化了吧
   ISTART = 0
   ICHARG = 2
   ENCUT = 380
   EDIFF = 1.0e-4
   LREAL =  auto

#Ionic relaxation
  NSW  =  900
  IBRION  =  0
  ISIF      =   0
  ISYM    =   0
  ISMAR  =  1
  SIGMA  =  0.1
  IALGO  =  48

# MD flags
  SMASS  =  2
  POTIM  =  2
  TEBEG  = 400
  TEEND  = 400
  LWAVE =    F
夜半酒醒人去后,更持红烛赏残花
6楼2013-11-26 19:14:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sxjsn1

金虫 (小有名气)

veni vedi veci

没人知道吗?
夜半酒醒人去后,更持红烛赏残花
7楼2013-11-27 15:43:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wl405103906

木虫 (初入文坛)

菜鸟

楼主你的问题解决了吗?我也出现这种问题,尝试修改参数,都不行。我做的是分子动力学。
vasp
8楼2014-01-03 15:35:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

零下1℃

铁虫 (正式写手)

【答案】应助回帖


sunyang1988: 金币+1, 谢谢交流 2014-01-03 21:55:53
建议楼主将IBRION设置为3试一下。当IBRION=3时,适应于分子动力学计算,原子原始坐标位置不是很准确的时候。希望能够有帮助。
9楼2014-01-03 15:59:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sxjsn1

金虫 (小有名气)

引用回帖:
8楼: Originally posted by wl405103906 at 2014-01-03 15:35:17
楼主你的问题解决了吗?我也出现这种问题,尝试修改参数,都不行。我做的是分子动力学。

没。。。
10楼2014-01-03 20:47:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sxjsn1 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 297求调剂 +8 戏精丹丹丹 2026-03-17 8/400 2026-03-18 14:30 by laoshidan
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考博] 环境领域全国重点实验室招收博士1-2名 +3 QGZDSYS 2026-03-13 5/250 2026-03-18 11:13 by QGZDSYS
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[基金申请] 被我言中:新模板不强调格式了,假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 296求调剂 +5 大口吃饭 身体健 2026-03-13 5/250 2026-03-17 21:05 by 不惑可乐
[考研] 材料与化工求调剂 +6 为学666 2026-03-16 6/300 2026-03-17 20:15 by peike
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 26考研求调剂 +6 丶宏Sir 2026-03-13 6/300 2026-03-17 16:13 by 醉在风里
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
信息提示
请填处理意见