计算化学公社

标题: VASP 结构优化收敛后进行单点计算不收敛 [打印本页]
作者
Author: Jiongci-21    时间: 2021-12-18 10:36
标题: VASP 结构优化收敛后进行单点计算不收敛
      各位老师和同学好，本人最近拿金属原子嵌入的氮掺杂石墨烯做氮还原反应，一些吸附中间体做了结构优化之后，直接拿CONTCAR ——> POSCAR 设置频率计算打算进行热力学矫正，某些结构的第一步的单点计算会不收敛，或者出现收敛了，但是stdout里的mag和E0会出现变化。如下图：

(, 下载次数 Times of downloads: 21)        (, 下载次数 Times of downloads: 27)      
这个是吸附的NH2中间体的结构优化（前）和频率计算（后）。


拿优化后的结构打开LWAVE算单点也会遇到不收敛或者收敛但E0和mag改变的情况（结构优化和单点计算的收敛标准一致不变(EDIFF = 1E-06; EDIFFG = -0.02) , 频率计算我会把收敛标准收紧一些(EDIFF = 1E-07; EDIFFG = -0.01)）
我重新打开LWAVE并进行结构优化可以再次得到相同结构，用输出的wavecar做频率计算就一切正常了，得到的频率里是没有虚频的。



请问如果出现这样的情况，那我得到优化后的结构它可靠吗？单点计算不收敛现象出现的原因是什么呢？（甚至有些结构，我拿优化后的结构再次优化，也可能会出现结构、能量、磁矩的改变（所有参数均不变只将CONTCAR——>POSCAR））


以下附上对应的INCAR参数：

  1. <div>OPT
  2. </div># start setting #
  3. ISTART = 1
  4. #ICHARG
  5. #INIWAV
  6. ###################################

  8. # spin polarized calculation #
  9. ISPIN = 2
  10. MAGMOM = 100*0.0 7*0.0 2*2.0 1*4.0 2*0.0
  11. ###################################

  13. # lda + u method #
  14. #LDAU =
  15. #LDAUTYPE =
  16. #LDAUL =
  17. #LDAUU
  18. #LDAUJ
  19. ##################################

  21. # electronic relaxation #
  22. ENCUT = 450
  23. NELM = 250
  24. NELMIN = 4
  25. EDIFF = 1E-06
  26. #GGA
  27. PREC = Normal
  28. LREAL = Auto
  29. #NBANDS
  30. #NELECT
  31. #################################

  33. # ionic relaxation #
  34. EDIFFG = -0.02
  35. NSW = 1000
  36. IBRION = 2
  37. ISIF = 2
  38. ISYM = 0
  39. POTIM = 0.15
  40. #################################

  42. # smearing method #
  43. ISMEAR = 0
  44. SIGMA = 0.05
  45. ################################

  47. # mixing method #
  48. #ADDGIRD
  49. ALGO = Very Fast
  50. #AMIX
  51. #BMIX
  52. #AMIX_MAG
  53. #BMIX_MAG
  54. LPLANE = .TRUE.
  55. LSCALU = .FALSE.
  56. NCORE = 4
  57. #NSIM
  58. #################################

  60. # out put #
  61. LWAVE = .FALSE.
  62. LCHARG = .FALSE.
  63. #LVTOT
  64. #LELF
  65. #################################

  67. IVDW = 12

  69. #IALGO
复制代码

  1. fre-calculation

  3. # start setting #
  4. ISTART = 1
  5. #ICHARG
  6. #INIWAV
  7. ###################################

  9. # spin polarized calculation #
  10. ISPIN = 2
  11. MAGMOM = 100*0.0 7*0.0 2*2.0 1*4.0 2*0.0
  12. ###################################

  14. # lda + u method #
  15. #LDAU =
  16. #LDAUTYPE =
  17. #LDAUL =
  18. #LDAUU
  19. #LDAUJ
  20. ##################################

  22. # electronic relaxation #
  23. ENCUT = 450
  24. NELM = 250
  25. NELMIN = 4
  26. EDIFF = 1E-07
  27. #GGA
  28. PREC = Normal
  29. LREAL = Auto
  30. #NBANDS
  31. #NELECT
  32. #################################

  34. # ionic relaxation #
  35. EDIFFG = -0.01
  36. NSW = 1000
  37. IBRION = 5
  38. ISIF = 2
  39. ISYM = 0
  40. POTIM = 0.015
  41. NFREE = 2
  42. #################################

  44. # smearing method #
  45. ISMEAR = 0
  46. SIGMA = 0.05
  47. ################################

  49. # mixing method #
  50. #ADDGIRD
  51. ALGO = Very Fast
  52. #AMIX
  53. #BMIX
  54. #AMIX_MAG
  55. #BMIX_MAG
  56. LPLANE = .TRUE.
  57. LSCALU = .FALSE.
  58. NCORE = 4
  59. #NSIM
  60. #################################

  62. # out put #
  63. LWAVE = .FALSE.
  64. LCHARG = .FALSE.
  65. #LVTOT
  66. #LELF
  67. #################################

  69. IVDW = 12

  71. #IALGO
复制代码

第一次发帖格式和表述可能会有问题，请各位批评指正。


作者
Author: jiangning198511    时间: 2021-12-18 10:57
本帖最后由 jiangning198511 于 2021-12-18 11:12 编辑

K点测试过吗？ 估计是包含过渡金属 导致每次收敛到不同的基态了 只能每次都读取WAVECAR 保证初始猜测一样了，
作者
Author: Jiongci-21    时间: 2021-12-18 13:33
jiangning198511 发表于 2021-12-18 10:57
K点测试过吗？ 估计是包含过渡金属 导致每次收敛到不同的基态了 只能每次都读取WAVECAR 保证初始猜测一样了 ...

K点测试过，目前我用的111（晶胞17*17*18），看体系绝对能量的话，331K点会收敛到小数点2位。但是用吸附氮气的相对能量做参考的话，111和221的吸附能差距是不到0.1 eV。但是时间差了快有4倍。我就选了前者。
作者
Author: jiangning198511    时间: 2021-12-20 11:25
Jiongci-21 发表于 2021-12-18 13:33
K点测试过，目前我用的111（晶胞17*17*18），看体系绝对能量的话，331K点会收敛到小数点2位。但是用吸附 ...

这么大的晶胞 111应该可以；111和221吸附能差0.1 eV有点大，是不是存在收敛的基态不同？再者类似的计算应该有很多，你可以参考一下他们的计算设置
作者
Author: Jiongci-21    时间: 2021-12-20 15:03
jiangning198511 发表于 2021-12-20 11:25
这么大的晶胞 111应该可以；111和221吸附能差0.1 eV有点大，是不是存在收敛的基态不同？再者类似的计算应 ...

好的，谢谢解答。



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3