关于cp2k中SCAN泛函的输入

yoya · 发表于 Post on 2019-1-2 21:01:29

本帖最后由 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>

yoya · 发表于 Post on 2019-1-7 10:59:10

已解决

mizuchi · 发表于 Post on 2019-1-7 16:54:08

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

yoya · 发表于 Post on 2019-1-10 15:43:02

mizuchi 发表于 2019-1-7 16:54
Interface to LIBXC was changed since 6.1. You should split &LIBXC section:
&LIBXC
...

谢谢，是的，后面改成了这种格式就没有报错了

s345538842 · 发表于 Post on 2020-12-1 17:30:50

mizuchi 发表于 2019-1-7 16:54
Interface to LIBXC was changed since 6.1. You should split &LIBXC section:
&LIBXC
...

请问楼主，是要两个&LIBXC 的section 都写么？？

k64_cc · 发表于 Post on 2020-12-1 20:05:12

s345538842 发表于 2020-12-1 17:30
请问楼主，是要两个&LIBXC 的section 都写么？？

是的，分别是交换项和关联项……

s345538842 · 发表于 Post on 2020-12-2 08:40:00

本帖最后由 s345538842 于 2020-12-2 09:17 编辑

k64_cc 发表于 2020-12-1 20:05
是的，分别是交换项和关联项……

非常感谢解答，在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

s345538842 · 发表于 Post on 2020-12-2 09:26:06

本帖最后由 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

COORD_FILE_NAME last_frame_wrapped.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



elpa · 发表于 Post on 2022-7-6 16:00:47

s345538842 发表于 2020-12-2 09:26
在group里又找到了这个：（https://groups.google.com/g/cp2k/c/k0M3XuOdIHI/m/TEeDFRMwAAAJ）
1) Many pe ...

关于这个输入文件，我有个问题想问一下。盒子大小是24.45，但是VDW的截断半径却是40，远远大于盒子边长，这合理吗？

sobereva · 发表于 Post on 2022-7-26 06:00:03

elpa 发表于 2022-7-6 16:00
关于这个输入文件，我有个问题想问一下。盒子大小是24.45，但是VDW的截断半径却是40，远远大于盒子边长， ...

色散校正的截断没有这个限制

elpa · 发表于 Post on 2022-7-26 17:40:38

sobereva 发表于 2022-7-26 06:00
色散校正的截断没有这个限制

请教一下，对立方盒子边长20A，BLYP或PBE泛函，VDW截断半径设多少比较合适？

sobereva · 发表于 Post on 2022-7-26 18:08:40

elpa 发表于 2022-7-26 17:40
请教一下，对立方盒子边长20A，BLYP或PBE泛函，VDW截断半径设多少比较合适？

用程序默认的就行了

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 关于cp2k中SCAN泛函的输入

浏览过的版块