24小时热门版块排行榜    

查看: 4433  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085601材料工程专硕求调剂 +5 慕寒mio 2026-03-16 5/250 2026-03-17 21:31 by hmn_wj
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 10/500 2026-03-17 15:07 by 一鸭鸭哟
[考研] 材料与化工专硕调剂 +5 heming3743 2026-03-16 5/250 2026-03-17 14:03 by 勇敢太监王公公
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 6/300 2026-03-16 15:58 by 漾漾123sun
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +7 Ncdx123456 2026-03-13 8/400 2026-03-16 12:15 by karry wen
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 工科,求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
信息提示
请填处理意见