about energy calculation error
chun
chun... at gmail.com
Fri Oct 5 20:51:08 UTC 2007
Hi,
I use cp2k to do all-electron energy calculation on a single
water
molecule with periodic boundry condition and unit cell size
9.4Angstrom X 9.4Angstrom X 9.4Angstrom (see the input file in the
end). When I move the molecule from the center (4.7,4.7,4.7) of the
cell to the corner (2.7,2.7,2.7) with the same structure, there is an
energy difference about 0.5 mili-Hartree between them. Could anyone
tell me what causes the energy difference?
Thanks.
Chun-Min
&FORCE_EVAL
METHOD Quickstep
&DFT
&MGRID
CUTOFF 158
&END MGRID
&QS
METHOD GAPW
FULL_GAPW
&END QS
&SCF
SCF_GUESS ATOMIC
&END SCF
&XC
&XC_FUNCTIONAL BLYP
&END XC_FUNCTIONAL
&END XC
&END DFT
&SUBSYS
&CELL
ABC 9.4 9.4 9.4
UNIT ANGSTROM
&END CELL
&COORD
O 4.700000 4.700000 4.634413
H 4.700000 3.942864 5.220545
H 4.700000 5.457136 5.220545
# O 2.700000 2.700000 2.634413
# H 2.700000 1.942864 3.220545
# H 2.700000 3.457136 3.220545
&END COORD
&KIND H
BASIS_SET 6-31G*
POTENTIAL ALL
&END KIND
&KIND O
BASIS_SET 6-31G*
POTENTIAL ALL
&END KIND
&END SUBSYS
&END FORCE_EVAL
&GLOBAL
PROJECT H2O
PRINT_LEVEL MEDIUM
&END GLOBAL
More information about the CP2K-user
mailing list