Energy of atomic oxgyen and smearing
Natalie Austin
natalie... at gmail.com
Thu Jul 14 02:05:19 UTC 2016
Hello Matthias,
The input file is below:
&FORCE_EVAL
METHOD QS
&DFT
BASIS_SET_FILE_NAME BASIS_MOLOPT
POTENTIAL_FILE_NAME GTH_POTENTIALS
&MGRID
CUTOFF 400
REL_CUTOFF 60
NGRIDS 4
&END MGRID
&QS
EPS_DEFAULT 1.0E-14
MAP_CONSISTENT
&END QS
&SCF
SCF_GUESS ATOMIC
EPS_SCF 1.0E-7
MAX_SCF 500
CHOLESKY INVERSE
&DIAGONALIZATION
ALGORITHM STANDARD
&END DIAGONALIZATION
&MIXING
METHOD DIRECT_P_MIXING
ALPHA 0.1 *# I also ran a case where I increased the value
to 0.4 as you suggested. *
BETA 1.5
NBROYDEN 8
&END MIXING
&END SCF
&XC
&XC_FUNCTIONAL PBE
&PBE
PARAMETRIZATION REVPBE
&END PBE
&END XC_FUNCTIONAL
&VDW_POTENTIAL
DISPERSION_FUNCTIONAL PAIR_POTENTIAL
&PAIR_POTENTIAL
TYPE DFTD3
REFERENCE_FUNCTIONAL PBE
CALCULATE_C9_TERM .TRUE.
PARAMETER_FILE_NAME dftd3.dat
R_CUTOFF 15.0
VERBOSE_OUTPUT .TRUE.
&END PAIR_POTENTIAL
&END vdw_POTENTIAL
&END XC
MULTIPLICITY 3
UKS .TRUE.
&END DFT
&SUBSYS
&CELL
ABC 20.00000 20.00000 20.00000
ALPHA_BETA_GAMMA 90.00 90.00 90.00
PERIODIC NONE
&END CELL
&TOPOLOGY
COORD_FILE_NAME o.xyz
COORDINATE XYZ
&END TOPOLOGY
&KIND O
ELEMENT O
BASIS_SET DZVP-MOLOPT-SR-GTH-q6
POTENTIAL GTH-PBE-q6
&END KIND
&END SUBSYS
&END FORCE_EVAL
&GLOBAL
PROJECT o
RUN_TYPE ENERGY
PRINT_LEVEL MEDIUM
&END GLOBAL
&MOTION
&GEO_OPT
MAX_FORCE 0.0004
MAX_ITER 2000
OPTIMIZER BFGS
TYPE MINIMIZATION
&END GEO_OPT
&END MOTION
and the coordinate file:
1
o
O 7.46880 7.39185 7.47705
Natalie
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20160713/8e2d41f7/attachment.htm>
More information about the CP2K-user
mailing list