[CP2K-user] How to set up an input file that uses the PBE0 functional: my version runs extremely slowly.
Nat
sf1... at gmail.com
Wed Feb 26 17:36:59 UTC 2020
Dear all:
I would like to use CP2K to run NPT water simulations that use different
functionals. I have so far used CP2K tutorials to set up simulations based
on the PBE and PBE0 functionals. The simulations that used the PBE
functional ran quite efficiently and I'm happy with the results. The
problem is that the simulations based on PBE0 are running so inefficiently
that I must have made a mistake somewhere in my input file (i.e., it will
take decades for the simulations to finish). If anyone out there is
experienced in using PBE0, I would be very grateful for any advice.
I used this tutorial to set up the simulation.
https://www.cp2k.org/exercises:2017_uzh_cp2k-tutorial:hybrid
Here is my input file. I have colored the exchange–correlation functional
part blue. If anything looks a little suspect, it probably is because I
don't really know what I'm doing as far as optimizing the solver settings
or choosing the best PBE0 specifications.
&GLOBAL
PROJECT pbe0
RUN_TYPE MD
&END GLOBAL
&FORCE_EVAL
STRESS_TENSOR ANALYTICAL
METHOD Quickstep
&DFT
BASIS_SET_FILE_NAME BASIS_MOLOPT
POTENTIAL_FILE_NAME POTENTIAL_SET
CHARGE 0
MULTIPLICITY 1
&MGRID
CUTOFF [Ry] 800
REL_CUTOFF [Ry] 60
&END
&QS
METHOD GPW
EPS_DEFAULT 1.0E-10
EXTRAPOLATION ASPC
EXTRAPOLATION_ORDER 4
&END QS
&POISSON
PERIODIC XYZ
&END
&PRINT
&E_DENSITY_CUBE OFF
&END E_DENSITY_CUBE
&MO_CUBES
NLUMO 4
NHOMO 4
WRITE_CUBE .FALSE.
&EACH
MD 10
&END
&END
&END
&SCF
SCF_GUESS ATOMIC
MAX_SCF 30
EPS_SCF 1.0E-6
&OT
PRECONDITIONER FULL_SINGLE_INVERSE
MINIMIZER DIIS
&END OT
&OUTER_SCF
MAX_SCF 10
EPS_SCF 1.0E-6
&END
&PRINT
&RESTART OFF
&END
&END
&END SCF
* ! Exchange-correlation approach*
&XC
! Use PBE0
&XC_FUNCTIONAL
&PBE
! 75% GGA exchange
SCALE_X 0.75
! 100% GGA correlation
SCALE_C 1.0
&END PBE
&END XC_FUNCTIONAL
&XC_GRID
XC_DERIV SPLINE2_SMOOTH
&END XC_GRID
&HF
! 25 % HFX exchange
FRACTION 0.25
&SCREENING
EPS_SCHWARZ 1.0E-6
SCREEN_ON_INITIAL_P TRUE
&END
&INTERACTION_POTENTIAL
POTENTIAL_TYPE TRUNCATED
CUTOFF_RADIUS 6.0
T_C_G_DATA ../t_c_g.dat
&END
&MEMORY
MAX_MEMORY 4000
EPS_STORAGE_SCALING 0.1
&END MEMORY
&END HF
&VDW_POTENTIAL
POTENTIAL_TYPE PAIR_POTENTIAL
&PAIR_POTENTIAL
PARAMETER_FILE_NAME dftd3.dat
TYPE DFTD3
REFERENCE_FUNCTIONAL PBE0
R_CUTOFF [angstrom] 16
&END
&END VDW_POTENTIAL
&END XC
&END DFT
&SUBSYS
&CELL
ABC 13.66 13.66 13.66
&END CELL
&COORD
O -7.6415350291705231E+00 6.0494971538864633E+00
8.1919697993062783E+00 H2O
H -6.6418362232982036E+00 5.9231018862335736E+00
8.4056822626120979E+00 H2O
H -8.0518670684031459E+00 6.3158684817228137E+00
9.0061271153555484E+00 H2O
O 2.6825296782034535E+01 2.0466107496669252E+01
1.7802574406580192E+01 H2O
H 2.6447741120699803E+01 2.0407102905755600E+01
1.6912157191231316E+01 H2O
H 2.6439964847387159E+01 2.1295072206802320E+01
1.8100927322697789E+01 H2O
! I have omitted most of the atoms for clarity. I have 300.
&END COORD
&TOPOLOGY
CONNECTIVITY GENERATE
&GENERATE
BONDLENGTH_MAX 7
&END
&END
&KIND H
BASIS_SET DZVP-MOLOPT-SR-GTH
POTENTIAL GTH-NLCC-PBE-q1
&END KIND
&KIND O
BASIS_SET DZVP-MOLOPT-SR-GTH
POTENTIAL GTH-NLCC-PBE-q6
&END KIND
&END SUBSYS
&END FORCE_EVAL
&MOTION
&MD
ENSEMBLE NPT_I
TEMPERATURE [K] 400.
TIMESTEP [fs] 0.5
STEPS 8000000
&THERMOSTAT
&NOSE
LENGTH 3
YOSHIDA 3
TIMECON 1000
MTS 2
&END NOSE
&END THERMOSTAT
&BAROSTAT
PRESSURE 1000.
&END BAROSTAT
&END MD
&PRINT
&CELL
&END
&RESTART_HISTORY OFF
&END
&END PRINT
&END MOTION
Thank you.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20200226/1795d173/attachment.htm>
More information about the CP2K-user
mailing list