24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhangfan2192

金虫 (著名写手)

应助: 9 (幼儿园)
金币: 1007.3
散金: 62
红花: 3
帖子: 1170
在线: 171.3小时
虫号: 648132
注册: 2008-11-07
性别: GG
专业: 理论和计算化学

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

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

不知道为什么，我在用neb计算过度态的时候总是感觉计算结果不是很理想，很难找到自己需要的结果，比如说计算的8个image的能量是单调上升或下降的，而不是抛物线形状，不知道大家在计算的时候遇到这样的情况没，如果遇到的话大家都是怎么解决的？相信现在用vasp软件的不在少数，欢迎大家积极讨论建议！

回复此楼

» 收录本帖的淘帖专辑推荐

VASP	first principle	第一性原理和电化学	vasp计算总结
计算资料收集byFrank

» 猜你喜欢

AI论文写作工具：是科研加速器还是学术作弊器？已经有5人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有7人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有5人回复
2025冷门绝学什么时候出结果已经有7人回复

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

vasp计算过渡态问题已经有13人回复
新手请教用NEB方法计算energy barrier. 已经有15人回复
vasp用neb计算迁移势垒已经有11人回复
vasp NEB 如何选取原胞大小？已经有10人回复
vasp 过渡态计算线性插值方法。已经有8人回复
【求助】vasp计算过渡态结果出错已经有17人回复
【文献学习】JACS上一篇纯DFT计算的文献： CI-NEB计算，Leuwdin布局分析，态密度分析已经有30人回复
【求助】关于vasp NEB 已经有8人回复
【求助】NEB算过渡态出现致命错误已经有5人回复
【求助成功】怎么样提高NEB计算速度? 已经有19人回复
【讨论】NEB寻找过渡态专题讨论已经有49人回复

相逢是缘！

1楼2010-07-05 21:46:51

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zhangfan2192

金虫 (著名写手)

应助: 9 (幼儿园)
金币: 1007.3
散金: 62
红花: 3
帖子: 1170
在线: 171.3小时
虫号: 648132
注册: 2008-11-07
性别: GG
专业: 理论和计算化学

引用回帖:

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