[CP2K-user] [CP2K:20945] Energy from restart file not equal to energy at the end of geo opt
Bram Van der Linden
bramvdlinden98 at gmail.com
Wed Dec 4 23:43:48 UTC 2024
Dear all,
I was following this tutorial: Geometry Optimisation — CP2K documentation
<https://manual.cp2k.org/trunk/methods/optimization/geometry.html>
If I perform the geo opt, I get the snippet below at the end of the output
file.
If I do an energy calculation from the restart file, I obtain that
ENERGY| Total FORCE_EVAL ( QS ) energy [a.u.]:
-17.164620477328342a
So, the energy of the energy calculation is equal to the energy at the last
step (11), it is not equal to the energy that is recomputed after the last
geometry optimization step.
Could someone explain why this is the case?
Many thanks in advance!
Bram
-------- Informations at step = 11 ------------
Optimization Method = SD
Total Energy = -17.1646204773
Real energy change = -0.0000000529
Decrease in energy = YES
Used time = 2.886
Convergence check :
Max. step size = 0.0003393777
Conv. limit for step size = 0.0010000000
Convergence in step size = YES
RMS step size = 0.0001493511
Conv. limit for RMS step = 0.0010000000
Convergence in RMS step = YES
Max. gradient = 0.0001787811
Conv. limit for gradients = 0.0010000000
Conv. in gradients = YES
RMS gradient = 0.0000786768
Conv. limit for RMS grad. = 0.0010000000
Conv. in RMS gradients = YES
---------------------------------------------------
*******************************************************************************
*** GEOMETRY OPTIMIZATION COMPLETED
***
*******************************************************************************
Reevaluating energy at the minimum
Number of electrons:
8
Number of occupied orbitals:
4
Number of molecular orbitals:
4
Number of orbital functions:
23
Number of independent orbital functions:
23
Extrapolation method: ASPC
SCF WAVEFUNCTION OPTIMIZATION
Step Update method Time Convergence Total energy
Change
------------------------------------------------------------------------------
1 Pulay/Diag. 0.50E+00 0.0 0.00005610 -17.1646204769
-1.72E+01
2 Pulay/Diag. 0.50E+00 0.1 0.00000563 -17.1646347273
-1.43E-05
*** SCF run converged in 2 steps ***
Electronic density on regular grids: -8.0000016293
-0.0000016293
Core density on regular grids: 7.9999992519
-0.0000007481
Total charge density on r-space grids: -0.0000023775
Total charge density g-space grids: -0.0000023775
Overlap energy of the core charge distribution:
0.00000004555422
Self energy of the core charge distribution:
-43.83289054591484
Core Hamiltonian energy:
12.82175609863938
Hartree energy:
17.97395116248504
Exchange-correlation energy:
-4.12745148803461
Total energy:
-17.16463472727081
ENERGY| Total FORCE_EVAL ( QS ) energy [a.u.]:
-17.164634727270808
--
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/cp2k/7bd75aef-e648-41ca-b4c0-cab444f39165n%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20241204/043fb1ea/attachment.htm>
More information about the CP2K-user
mailing list