24小时热门版块排行榜    

Znn3bq.jpeg
查看: 4149  |  回复: 20
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zhangfan2192

金虫 (著名写手)

[交流] 【其他】vasp neb 计算讨论 已有11人参与

希望各位用vasp的虫友能够讨论下neb计算中的注意事项:

不知道为什么,我在用neb计算过度态的时候总是感觉计算结果不是很理想,很难找到自己需要的结果,比如说计算的8个image的能量是单调上升或下降的,而不是抛物线形状,不知道大家在计算的时候遇到这样的情况没,如果遇到的话大家都是怎么解决的?相信现在用vasp软件的不在少数,欢迎大家积极讨论建议!
回复此楼
相逢是缘!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangfan2192

金虫 (著名写手)

引用回帖:
Originally posted by fanchen021 at 2010-07-07 10:30:40:
最好能给出INCAR。

INCAR文件如下:

   SYSTEM = Sul-Fe
  Startparameter for this Run:
   NWRITE =     2     verbosity write-flag
   LPETIM=      F     write-flag & timer #
   ISTART =     0     0-new 1-cont 2-samecut job
   ISPIN  =     2     1-no, 2-yes
   INIWAV =     1     0-jellium WF; 1-random numbers (recomm.)

  Electronic Relaxation
!   ENCUT  =   400.00   #Ecut (eV)
   PREC   =    Medium  *  #precission: #Low/Medium/High
   NELM   =    200      max. # of electronic SC steps
   NELMDL  =    5      # of non-consistent steps at the beggining:>0 or <0
   EDIFF  = 1.0e-4    stopping-criterion for electronic upd.
   VOSKOWN = 1
   ISMEAR = 1
   SIGMA = 0.1
!   AMIX = 0.2
!   BMIX = 0.0001

  Ionic relaxation
   EDIFFG =  1.0e-3    stopping-conditions for ionic relax.loop (EDIFF*10)
   NSW    =    200      *  # of steps for ionic update (def:0)
   IBRION =    2   *  ionic relax: -1:no move; 0-MD;1-Newton;2-CG;3-damped;4-?
   ISIF   =    2   *  calc.stress:0-no;1-tr;2-7 yes; 3-change vol.4-change shape
   ISYM   =    0      symmetry:0-OFF, 1-ON,2-ON+memory conserving for CD

   POTIM  =   0.3     if IBRION=1,2,3:scaling constant forces;IBRION=0,time st.


   LPLANE = .TRUE.
   NPAR = 8
   IALGO  =   48    * algorithm: use only 8 (CG) or 48 (RMM-DIIS)
   LDIAG  = .TRUE.    sub-space diagonalisation
   GGA    =   91    * xc-type: PB, PW (86), LM, 91 (Perdew-Wang 91)

   LWAVE  =  .F.     * write WAVECAR
   LCHARG =  .T.    * write CHGCAR and CHG
   LVTOT  =  .F.      * write the local potential LOCPOT
   LELF   = .F.   create ELFCAR file
   LORBIT =   .F.      create PROOUT
   LREAL = Auto
! NEB
  IMAGES = 8              ! Number of images in the band, excluding endpoints
  SPRING = -5.0 eV/A2     ! Spring force between images; negative value turns on nudging
  LCLIMB = .False.         ! Flag to turn on the climbing image algorithm
相逢是缘!
6楼2010-07-07 10:40:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 21 个回答

stractor

金虫 (著名写手)

★ ★
lvjian8596(金币+1):感谢交流! 2010-07-05 22:38:22
xiaohunhun(金币+1):xiexie 2010-07-06 07:30:38
zhangfan2192(金币+5):谢谢 2010-07-06 09:20:03
能量表现为单调变化有可能出现。如始态和末态的某一个构型不是能量局域极小态的情况。
2楼2010-07-05 22:16:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主

★ ★ ★
gavinliu7390(金币+2):谢谢交流! 2010-07-06 11:25:33
xiaohunhun(金币+1):谢谢 2010-07-06 12:34:51
zhangfan2192(金币+5): 2010-07-06 14:08:52
做NEB要保证初态和末态处于能量极小值点(局域的或是全局的),也就是初态和末态都是经过优化之后的,在这个基础上再做中间的imagines可能就可以避免单调上升和下降了。个人观点,欢迎指正。

[ Last edited by uuv2010 on 2010-7-6 at 10:23 ]
3楼2010-07-06 10:20:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangfan2192

金虫 (著名写手)

引用回帖:
Originally posted by uuv2010 at 2010-07-06 10:20:14:
做NEB要保证初态和末态处于能量极小值点(局域的或是全局的),也就是初态和末态都是经过优化之后的,在这个基础上再做中间的imagines可能就可以避免单调上升和下降了。个人观点,欢迎指正。

[ Last edited ...

谢谢回复,我的构型都是经过优化之后才算NEB的,不知道除了能量极小点这个因素外还有没有其他因素。
相逢是缘!
4楼2010-07-06 14:10:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上和省基金B类撒花 +17 花田半亩~白 2026-04-21 17/850 2026-04-23 07:30 by john198197
[考研] 有没有学校收留 +3 蒋昌鹏qtj 2026-04-20 3/150 2026-04-22 20:25 by 学员JpLReM
[考研] 312求调剂 +3 山河似你温柔 2026-04-22 3/150 2026-04-22 20:17 by 学员JpLReM
[考博] 华师大读博 +3 xq83 2026-04-22 5/250 2026-04-22 10:42 by xq83
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +5 yexuqing 2026-04-19 5/250 2026-04-22 10:01 by easeheart
[论文投稿] 急需审稿人!!! +3 陆小果画大饼 2026-04-21 3/150 2026-04-21 23:54 by jzy_123456
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 8/400 2026-04-21 16:26 by 啃面包的小书虫
[考研] 295分求调剂 +6 ?要上岸? 2026-04-17 6/300 2026-04-21 08:18 by Equinoxhua
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[论文投稿] 有没有接收比较快的sci期刊呀,最好在一个月之内的,研三孩子求毕业 20+4 之护着 2026-04-16 7/350 2026-04-20 15:45 by 豆豆7758
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
信息提示
请填处理意见