SCF energy not identical to result of final energy of geometry optimization
Frank von Horsten
xxx.f... at gmail.com
Wed May 28 11:24:24 UTC 2014
Dear all,
I am recalculating the SCF energy of a previously optimized geometry.
However, the results of the single point calculation differ quite
significantly from the final energy of the optimization.
More specifically, the final energy of the geo opt reads
E=-3578.782615929515032. This is the final in the line "ENERGY| Total
FORCE_EVAL" line. Already here the energy is slightly different from the
subsequent "Reevaluating energy at the minimum" part.
There the SCF cycle gives -3578.782584007462901. Admittedly, that is only
0.084 kJ/mol but still I find this a bit strange.dir
When I use the same input file just changing GEO_OPT to ENERGY_FORCE and
taking the final structure of the optimization I get an energy of
-3578.778182659758386, which is lower by 11.5 kJ/mol.
Also the calculated forced differ. For the optization the last forces are
SUM OF ATOMIC FORCES -0.00006050 -0.00001418 -0.00343851
0.00343908
Whereas for the single point calculatin I get
SUM OF ATOMIC FORCES -0.00000811 -0.00001851 -0.01197202
0.01197204
The system is a 4x4 Cu(100) slab (64 Cu atoms) with a c(2x2) chloride
adlayer on top (8 Cl atoms) and eight Na atoms as counter ions in the
vacuum region.
Does anyone have an idea where these discrapancies come from?
Cheers,
Frank
PS: I have attached the input and output files in case someone wants to
take a closer look.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20140528/ccaa25c6/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: energy.in
Type: application/octet-stream
Size: 18869 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20140528/ccaa25c6/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: energy.out
Type: application/octet-stream
Size: 254582 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20140528/ccaa25c6/attachment-0001.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: opt.in
Type: application/octet-stream
Size: 18880 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20140528/ccaa25c6/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: opt.out
Type: application/octet-stream
Size: 735281 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20140528/ccaa25c6/attachment-0003.obj>
More information about the CP2K-user
mailing list