计算化学公社

标题: 求助:关于DFT-D3使用无法输出正确能量的问题 [打印本页]
作者
Author: fux    时间: 2023-10-25 10:59
标题: 求助:关于DFT-D3使用无法输出正确能量的问题
因为要考虑有机无机界面问题，所以加IVDW=12，但是输出能量值有问题
  1.   FREE ENERGIE OF THE ION-ELECTRON SYSTEM (eV)
  2.   ---------------------------------------------------
  3.   free  energy   TOTEN  = ****************** eV

  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)***********

  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
复制代码


作者
Author: fux    时间: 2023-10-27 16:39
本帖最后由 fux 于 2023-10-27 16:44 编辑

问题解决了，为了防止大家遇到类似的问题，还是打算回一下这个贴。因为这个问题很刁钻，也没有在前面提到过。
我算的体系含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元素的朋友们都能不踩这个坑！
作者
Author: 903795249    时间: 2024-6-3 18:54
请问一下，在vasp怎么样实现DFT-D3的方法啊？
作者
Author: 903795249    时间: 2024-6-3 19:20
903795249 发表于 2024-6-3 18:54
请问一下，在vasp怎么样实现DFT-D3的方法啊？

找到了



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