CP2K中LIBXC泛函测试(与高斯比较)

jiangning198511

测试体系：H2O 基组：6-31g**

总能比较：(单位：au)

	G09	CP2K
PBE	-76.3327487	-76.3310874
TPSS	-76.4204004	-76.4107116
PBE0	-76.3369363	-76.3357600

可以看出PBE和PBE0的结果差别比较小，TPSS的差别较大。估计差别来源于积分精度。

欢迎讨论

输入文件例子：

G09:

  1 #pbepbe/6-31G**

  2

  3 H2O

  4

  5 0 1

  6     O   0.000000    0.000000   -0.065587

  7     H   0.000000   -0.757136    0.520545

  8     H   0.000000    0.757136    0.520545

CP2K:

&FORCE_EVAL

  METHOD Quickstep

  &DFT

    BASIS_SET_FILE_NAME EMSL_BASIS_SETS

    POTENTIAL_FILE_NAME POTENTIAL

    &MGRID

      CUTOFF 250

      REL_CUTOFF 50

    &END MGRID

    &QS

      METHOD GAPW

    &END QS

    &POISSON

      PERIODIC NONE

      PSOLVER MT

    &END

    &SCF

      EPS_SCF 1.0E-6

      SCF_GUESS ATOMIC

      MAX_SCF 100

    &END SCF

#TPSS-function

      &XC

      &XC_FUNCTIONAL

        &LIBXC

          FUNCTIONAL XC_GGA_X_PBE XC_GGA_C_PBE

        &END LIBXC

      &END XC_FUNCTIONAL

    &END XC

  &END DFT

  &SUBSYS

    &CELL

      ABC 6.0 6.0 6.0

      PERIODIC NONE

    &END CELL

    &COORD

    O   0.000000    0.000000   -0.065587

    H   0.000000   -0.757136    0.520545

    H   0.000000    0.757136    0.520545

    &END COORD

    &KIND H

      BASIS_SET 6-31Gxx

      POTENTIAL ALL

    &END KIND

    &KIND O

      BASIS_SET 6-31Gxx

      POTENTIAL ALL

    &END KIND

  &END SUBSYS

&END FORCE_EVAL

&GLOBAL

  PROJECT H2O-TPSS

#  TRACE

  PRINT_LEVEL MEDIUM

&END GLOBAL

genies

pbe和pbe0的差别主要来自于基组。G09用的是6D，而CP2K用的是5D，所以G09算出来的能量要低一些。如果都用5D，能量的差别可以减小到10-5次方以下，按我的标准，应该是10-8次分以下才对。实际上Cp2k的积分格点要比G09粗不少，但似乎CP2K提高积分精度之后，反而SCF收敛更难了，甚至水分子都不收敛。
TPSS等泛函等能量差就完全不能用积分精度来解释了。总之，肯定是某个地方不对。

jiangning198511

genies 发表于 2016-11-7 19:58
pbe和pbe0的差别主要来自于基组。G09用的是6D，而CP2K用的是5D，所以G09算出来的能量要低一些。如果都用5D ...

CP2K官方论坛的解释
Hi

so this is the GAPW bug for meta-functionals. As I have mentioned
before, there is an inconsistency in the way the kinetic energy density
is handled in the GAPW one-center terms.

PS:论坛不好上，国内的墙太厉害了

genies

jiangning198511 发表于 2016-11-8 16:17
CP2K官方论坛的解释
Hi

thanks

yoya

您好，想问您一下cp2k中关于SCAN泛函的输入

mizuchi

&LIBXC
  FUNCTIONAL  XC_MGGA_X_SCAN
&END LIBXC

&LIBXC
  FUNCTIONAL  XC_MGGA_C_SCAN
&END LIBXC

Libesi

jiangning198511 发表于 2016-11-8 16:17
CP2K官方论坛的解释
Hi

感谢