24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 4478  |  回复: 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的回帖

snail594

金虫 (小有名气)

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

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 286求调剂 +12 PolarBear11 2026-03-26 12/600 2026-03-28 12:14 by zllcz
[考研] 张芳铭-中国农业大学-环境工程专硕-298 +4 手机用户 2026-03-26 4/200 2026-03-28 07:17 by mmm just
[考研] 295求调剂 +5 1428151015 2026-03-27 6/300 2026-03-28 04:04 by fmesaito
[考研] 295材料工程专硕求调剂 +6 1428151015 2026-03-27 6/300 2026-03-28 03:51 by fmesaito
[考研] 一志愿 西北大学 总分282 英语一62 求调剂 +7 18419759900 2026-03-25 8/400 2026-03-27 16:38 by 18419759900
[考研] 材料与化工085600,总分304,本科有两篇sci参与,求调剂 +10 幸运的酱酱 2026-03-22 12/600 2026-03-27 16:08 by muchong357
[考研] 化学308分求调剂 +8 你好明天你好 2026-03-23 9/450 2026-03-27 14:01 by 杨光于青云
[考研] 303求调剂 +7 安忆灵 2026-03-22 8/400 2026-03-27 11:46 by sanrepian
[考研] 085600材料与化工306 +10 z1z2z3879 2026-03-21 11/550 2026-03-27 11:31 by wangjy2002
[考研] 349求调剂 +4 李木子啊哈哈 2026-03-25 4/200 2026-03-26 22:49 by fmesaito
[考研] 调剂 +4 柚柚yoyo 2026-03-26 4/200 2026-03-26 20:43 by fmesaito
[考研] 085602 289分求调剂 +8 WWW西西弗斯 2026-03-24 8/400 2026-03-26 16:33 by 不吃魚的貓
[考研] 一志愿哈工大,085400,320,求调剂 +4 gdlf9999 2026-03-24 4/200 2026-03-25 23:01 by boxking200
[考研] 考研调剂 +6 来好运来来来 2026-03-21 7/350 2026-03-25 22:43 by 418490947
[考研] 318求调剂 +5 plum李子 2026-03-21 8/400 2026-03-25 09:26 by aa331100
[考研] 求调剂 +6 研研,接电话 2026-03-24 7/350 2026-03-24 17:01 by barlinike
[考研] 300求调剂,材料科学英一数二 +5 leaflight 2026-03-24 5/250 2026-03-24 16:25 by laoshidan
[考研] 一志愿河北工业大学0817化工278分求调剂 +7 jhybd 2026-03-23 12/600 2026-03-24 09:03 by jhybd
[基金申请] 请教下大家 2026年国家基金申请是双盲审吗? +3 lishucheng1 2026-03-22 5/250 2026-03-24 08:22 by gltch
[考研] 328求调剂 +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
信息提示
请填处理意见