求助:关于DFT-D3使用无法输出正确能量的问题

fux

因为要考虑有机无机界面问题，所以加IVDW=12，但是输出能量值有问题

1.   FREE ENERGIE OF THE ION-ELECTRON SYSTEM (eV)
   
2.   ---------------------------------------------------
   
3.   free  energy   TOTEN  = ****************** eV
   
4. 
   
5.   energy  without entropy=******************  energy(sigma->0) =******************
   

复制代码

查到说这种输出一般是因为值过大，在一些字段无法显示，OSZICAR可以看到：

1.       1 T=   300. E= 0.24773043E+99 F= 0.24773043E+99 E0= 0.24773043E+99  EK= 0.11789E+02 SP= 0.00E+00 SK= 0.00E+00
   

复制代码

这个值显然是有问题的，检查OUTCAR发现：

1.          DFTD3 V3.0 Rev 1
   
2. IVDW         = 12
   
3. DF pbe
   
4. parameters
   
5. VDW_S6       =    1.0000
   
6. VDW_S8       =    0.7875
   
7. VDW_A1       =    0.4289
   
8. VDW_A2       =    4.4407
   
9. k1-k3        =   16.0000    1.3333   -4.0000
   
10. VDW_RADIUS   =   50.2022 A
    
11. VDW_CNRADIUS =   21.1671 A
    
12. Edisp (eV)***********
    
13. 
    
14. E6    (eV) :***********
    
15. E8    (eV) :    -5.9189
    
16. % E8        : -0.00
    
17.     FORVDW:  cpu time     11.9764: real time     12.0124
    

复制代码

就相当于从这里开始就出问题了，如果是纯有机相计算就不会有问题，无机相和界面相都是有问题的，下面是测试的INCAR，请大家帮忙检查：

1. ENCUT  = 520  
   
2. ENAUG  = 780
   
3. LREAL  = .FALSE.   
   
4. IALGO  = 38
   
5. ADDGRID = .TRUE.
   
6. NSIM   = 4
   
7. ISMEAR = 0   
   
8. SIGMA  = 0.1
   
9. IBRION = 2
   
10. ISIF   = 2
    
11. NSW    = 1
    
12. NELMIN = 5
    
13. EDIFF  = 1E-5
    
14. EDIFFG = -0.01
    
15. LWAVE = .FALSE.
    
16. LCHARG = .FALSE.
    
17. POTIM = 0.8
    
18. IVDW = 12
    

复制代码

903795249

903795249 发表于 2024-6-3 18:54
请问一下，在vasp怎么样实现DFT-D3的方法啊？

找到了

903795249

请问一下，在vasp怎么样实现DFT-D3的方法啊？

fux

问题解决了，为了防止大家遇到类似的问题，还是打算回一下这个贴。因为这个问题很刁钻，也没有在前面提到过。
我算的体系含Zr，出错的计算中OUTCAR有一个非常不显眼的WARNING

1. WARNING: type information on POSCAR and POTCAR are incompatible
   
2. POTCAR overwrites the type information in POSCAR
   
3. typ   5 type information:  Zr r
   
4. LDA part: xc-table for Pade appr. of Perdew
   

复制代码

这里提醒POTCAR我才去看，然后发现POTCAR

1.   PAW_PBE Zr_sv 04Jan2005
   
2.    12.0000000000000
   
3. parameters from PSCTR are:
   
4.    VRHFIN =r: 4s4p5s4d
   
5.    LEXCH  = PE
   
6.    EATOM  =  1284.2219 eV,   94.3876 Ry
   

复制代码

VRHFIN 后面少了一个Z，添上就正常了，在其他计算中这个似乎不影响，这个计算应该是要靠这里识别的元素去对应色散修正中的一些参数，所以这个对应出了问题，导致计算值错误。


希望各位算含Zr元素的朋友们都能不踩这个坑！