计算化学公社
标题: CP2K中LIBXC泛函测试(与高斯比较) [打印本页]
作者Author: jiangning198511 时间: 2016-11-4 10:03
标题: CP2K中LIBXC泛函测试(与高斯比较)
本帖最后由 jiangning198511 于 2016-11-4 10:17 编辑
测试体系：H2O 基组：6-31g**
总能比较：(单位：au)
可以看出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
作者Author: genies 时间: 2016-11-7 19:58
pbe和pbe0的差别主要来自于基组。G09用的是6D，而CP2K用的是5D，所以G09算出来的能量要低一些。如果都用5D，能量的差别可以减小到10-5次方以下，按我的标准，应该是10-8次分以下才对。实际上Cp2k的积分格点要比G09粗不少，但似乎CP2K提高积分精度之后，反而SCF收敛更难了，甚至水分子都不收敛。
TPSS等泛函等能量差就完全不能用积分精度来解释了。总之，肯定是某个地方不对。
作者Author: jiangning198511 时间: 2016-11-8 16:17
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:论坛不好上，国内的墙太厉害了
作者Author: genies 时间: 2016-11-14 21:13
thanks
作者Author: yoya 时间: 2019-1-5 12:25
您好，想问您一下cp2k中关于SCAN泛函的输入
作者Author: mizuchi 时间: 2019-1-5 19:09
&LIBXC
FUNCTIONAL XC_MGGA_X_SCAN
&END LIBXC
&LIBXC
FUNCTIONAL XC_MGGA_C_SCAN
&END LIBXC
作者Author: Libesi 时间: 2019-1-15 11:00
感谢
| 欢迎光临 计算化学公社 (http://ccc.keinsci.com/) |
Powered by Discuz! X3.3 |