VASP 结构优化收敛后进行单点计算不收敛

Jiongci-21

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

            
这个是吸附的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. ###################################
   
7. 
   
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. ###################################
    
12. 
    
13. # lda + u method #
    
14. #LDAU =
    
15. #LDAUTYPE =
    
16. #LDAUL =
    
17. #LDAUU
    
18. #LDAUJ
    
19. ##################################
    
20. 
    
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. #################################
    
32. 
    
33. # ionic relaxation #
    
34. EDIFFG = -0.02
    
35. NSW = 1000
    
36. IBRION = 2
    
37. ISIF = 2
    
38. ISYM = 0
    
39. POTIM = 0.15
    
40. #################################
    
41. 
    
42. # smearing method #
    
43. ISMEAR = 0
    
44. SIGMA = 0.05
    
45. ################################
    
46. 
    
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. #################################
    
59. 
    
60. # out put #
    
61. LWAVE = .FALSE.
    
62. LCHARG = .FALSE.
    
63. #LVTOT
    
64. #LELF
    
65. #################################
    
66. 
    
67. IVDW = 12
    
68. 
    
69. #IALGO
    

复制代码

1. fre-calculation
   
2. 
   
3. # start setting #
   
4. ISTART = 1
   
5. #ICHARG
   
6. #INIWAV
   
7. ###################################
   
8. 
   
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. ###################################
    
13. 
    
14. # lda + u method #
    
15. #LDAU =
    
16. #LDAUTYPE =
    
17. #LDAUL =
    
18. #LDAUU
    
19. #LDAUJ
    
20. ##################################
    
21. 
    
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. #################################
    
33. 
    
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. #################################
    
43. 
    
44. # smearing method #
    
45. ISMEAR = 0
    
46. SIGMA = 0.05
    
47. ################################
    
48. 
    
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. #################################
    
61. 
    
62. # out put #
    
63. LWAVE = .FALSE.
    
64. LCHARG = .FALSE.
    
65. #LVTOT
    
66. #LELF
    
67. #################################
    
68. 
    
69. IVDW = 12
    
70. 
    
71. #IALGO
    

复制代码

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

Jiongci-21

jiangning198511 发表于 2021-12-20 11:25
这么大的晶胞 111应该可以；111和221吸附能差0.1 eV有点大，是不是存在收敛的基态不同？再者类似的计算应 ...

好的，谢谢解答。

jiangning198511

Jiongci-21 发表于 2021-12-18 13:33
K点测试过，目前我用的111（晶胞17*17*18），看体系绝对能量的话，331K点会收敛到小数点2位。但是用吸附 ...

这么大的晶胞 111应该可以；111和221吸附能差0.1 eV有点大，是不是存在收敛的基态不同？再者类似的计算应该有很多，你可以参考一下他们的计算设置

Jiongci-21

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

K点测试过，目前我用的111（晶胞17*17*18），看体系绝对能量的话，331K点会收敛到小数点2位。但是用吸附氮气的相对能量做参考的话，111和221的吸附能差距是不到0.1 eV。但是时间差了快有4倍。我就选了前者。

jiangning198511

K点测试过吗？ 估计是包含过渡金属 导致每次收敛到不同的基态了 只能每次都读取WAVECAR 保证初始猜测一样了，