Energy Computation
August Melcher
august... at gmail.com
Sat Mar 8 00:40:44 UTC 2014
Hello, I am trying to compute the energy of a single water molecule, but
the I just keep getting a sum of isolated oxygen and two hydrogen atom
energies.
Here is my input:
&FORCE_EVAL
METHOD Quickstep
&DFT
UKS
BASIS_SET_FILE_NAME /usr/common/usg/cp2k/2.1.534/tests/QS/GTH_BASIS_SETS
POTENTIAL_FILE_NAME
/usr/common/usg/cp2k/2.1.534/tests/QS/GTH_POTENTIALS
&MGRID
NGRIDS 6
CUTOFF 300
&END MGRID
&QS
WF_INTERPOLATION ASPC
EXTRAPOLATION_ORDER 3
&END QS
&SCF
EPS_SCF 1.E-6
SCF_GUESS ATOMIC
MAX_SCF 400
&DIAGONALIZATION ON
ALGORITHM STANDARD
&END DIAGONALIZATION
&MIXING T
METHOD BROYDEN_MIXING
ALPHA 0.4
NBROYDEN 8
&END MIXING
&SMEAR ON
METHOD FERMI_DIRAC
ELECTRONIC_TEMPERATURE [K] 300
&END SMEAR
ADDED_MOS 10
&END SCF
&XC
&XC_FUNCTIONAL PBE
&END XC_FUNCTIONAL
&END XC
&PRINT
&END PRINT
&END DFT
&SUBSYS
&CELL
A 5.000000000 0.000000000 0.000000000
B 0.000000000 5.000000000 0.000000000
C 0.000000000 0.000000000 5.000000000
&END CELL
&COORD
O 2.790689766 2.0 2.0
H 2.0 2.0 2.612217330
H 3.581379532 2.0 2.612217330
&END COORD
&KIND O
BASIS_SET DZVP-GTH
POTENTIAL GTH-PBE-q6
&END KIND
&KIND H
BASIS_SET DZVP-GTH
POTENTIAL GTH-PBE-q1
&END KIND
&END SUBSYS
&END FORCE_EVAL
&GLOBAL
PROJECT water
RUN_TYPE ENERGY
PRINT_LEVEL LOW
&END GLOBAL
Both the box dimensions and coordinates are meant to be in angstroms. The
energy should be around -76 Eh, but this input is giving energies around
-16.8Eh. Any hints as to what I am doing wrong would be greatly appreciated.
Thanks,
August
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20140307/78593245/attachment.htm>
More information about the CP2K-user
mailing list