计算化学公社
标题: 关于cp2k中SCAN泛函的输入 [打印本页]
作者Author: yoya 时间: 2019-1-2 21:01
标题: 关于cp2k中SCAN泛函的输入
本帖最后由 yoya 于 2019-1-4 19:56 编辑
大家好，我是cp2k的初学者，我想用SCAN（meta-GGA）泛函算下水分子，但是找不到该泛函在官网中输入文件的例子，请求各位帮助！
我自己写了一下SCAN部分的输入，但是报错了
输入部分为：
&XC
&XC FUNCTIONAL
&LIBXC
FUNCTIONAL XC_MGGA_X_SCAN XC_MGGA_C_SCAN
&END LIBXC
&END XC FUNCTIONAL
&XC
报错如下：found unexpected extra argument XC_MGGA_C_SCAN at , File: 'SCAN,inp', Line:55, Column:3, Chunk: <XC_MGGA_X_SCAN>
作者Author: yoya 时间: 2019-1-7 10:59
已解决
作者Author: mizuchi 时间: 2019-1-7 16:54
Interface to LIBXC was changed since 6.1. You should split &LIBXC section:
&LIBXC
FUNCTIONAL XC_MGGA_X_SCAN
&END LIBXC
&LIBXC
FUNCTIONAL XC_MGGA_C_SCAN
&END LIBXC
作者Author: yoya 时间: 2019-1-10 15:43
谢谢，是的，后面改成了这种格式就没有报错了
作者Author: s345538842 时间: 2020-12-1 17:30
请问楼主，是要两个&LIBXC 的section 都写么？？
作者Author: k64_cc 时间: 2020-12-1 20:05
是的，分别是交换项和关联项……
作者Author: s345538842 时间: 2020-12-2 08:40
本帖最后由 s345538842 于 2020-12-2 09:17 编辑
非常感谢解答，在cp2k的Google group找到一个CH4例子，用的还是PBE的[color=rgba(0, 0, 0, 0.87)]pseudopotentials，这样是可以的么？？？我的疑问是scan泛函对应的[color=rgba(0, 0, 0, 0.87)]pseudopotentials是啥？我对比了PBE和SCAN的计算结果（都用PBE的[color=rgba(0, 0, 0, 0.87)]pseudopotentials），相同体系ener文件中输出的potential差20多。
&XC
&XC_FUNCTIONAL
&LIBXC
FUNCTIONAL MGGA_X_SCAN
&END LIBXC
&LIBXC
FUNCTIONAL MGGA_C_SCAN
&END LIBXC
&END XC_FUNCTIONAL
&END XC
&PRINT
&MO_CUBES ON
NHOMO 1
NLUMO 1
&END MO_CUBES
&E_DENSITY_CUBE ON
&END E_DENSITY_CUBE
&END PRINT
&END DFT
&SUBSYS
&CELL
ABC 10.0 10.0 10.0
PERIODIC none
&END CELL
&COORD
UNIT bohr
C 0.00000000000000 0.00000000000000 0.00000000000000
H 1.18771160655551 -1.18771160655551 1.18771160655551
H -1.18771160655551 1.18771160655551 1.18771160655551
H -1.18771160655551 -1.18771160655551 -1.18771160655551
H 1.18771160655551 1.18771160655551 -1.18771160655551
&END COORD
&KIND C
BASIS_SET DZVP-MOLOPT-SR-GTH
ELEMENT C
POTENTIAL GTH-PBE-q4
&END KIND
&KIND H
ELEMENT H
BASIS_SET DZVP-MOLOPT-SR-GTH
POTENTIAL GTH-PBE-q1
&END KIND
作者Author: s345538842 时间: 2020-12-2 09:26
本帖最后由 s345538842 于 2020-12-2 09:30 编辑
在group里又找到了这个：（https://groups.google.com/g/cp2k/c/k0M3XuOdIHI/m/TEeDFRMwAAAJ）
1) Many people have used PBE pp previously.
SCAN optimized pp can be found at
https://github.com/juerghutter/GTH/blob/master/SCAN/POTENTIAL
2) Newly published SCAN parameters for D3 are available in the current
version of CP2K from Github. For older version you need to
add a line in the input with the parameters.
&XC
&XC_FUNCTIONAL
&LIBXC
FUNCTIONAL MGGA_X_SCAN
&END LIBXC
&LIBXC
FUNCTIONAL MGGA_C_SCAN
&END LIBXC
&END XC_FUNCTIONAL
&vdW_POTENTIAL
DISPERSION_FUNCTIONAL PAIR_POTENTIAL
&PAIR_POTENTIAL
R_CUTOFF 40.0
TYPE DFTD3
D3_SCALING 1.0 1.324 0.0
PARAMETER_FILE_NAME dftd3.dat
! REFERENCE_FUNCTIONAL SCAN
&END PAIR_POTENTIAL
&END vdW_POTENTIAL
&END XC
&END DFT
&SUBSYS
&CELL
ABC 24.450 24.450 24.450
PERIODIC XYZ
&END CELL
&TOPOLOGY
COORD_FILE_FORMAT XYZ
&END TOPOLOGY
&KIND Mg
BASIS_SET DZVP-MOLOPT-SR-GTH
POTENTIAL GTH-SCAN-q10
&END KIND
&KIND Cl
BASIS_SET DZVP-MOLOPT-SR-GTH
POTENTIAL GTH-SCAN-q7
&END KIND
&END SUBSYS
作者Author: elpa 时间: 2022-7-6 16:00
关于这个输入文件，我有个问题想问一下。盒子大小是24.45，但是VDW的截断半径却是40，远远大于盒子边长，这合理吗？
作者Author: sobereva 时间: 2022-7-26 06:00
色散校正的截断没有这个限制
作者Author: elpa 时间: 2022-7-26 17:40
请教一下，对立方盒子边长20A，BLYP或PBE泛函，VDW截断半径设多少比较合适？
作者Author: sobereva 时间: 2022-7-26 18:08
用程序默认的就行了
| 欢迎光临 计算化学公社 (http://bbs.keinsci.com/) |
Powered by Discuz! X3.3 |