<div dir="ltr"><span style="font-size:13px">Thanks for your reply.</span> <br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Jan 12, 2020 at 1:31 AM Patrick Gono <<a href="mailto:patri...@gmail.com">patri...@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Dear Khodayar,<div><br></div><div>The reevaluation of the energy at the end of the optimization might yield a different total energy value if your calculation is not converged. If that is the case, you might choose more stringent values for EPS_SCF inside the SCF subsection, which governs the convergence of the SCF loop, and EPS_DEFAULT, which is tied to the numerical accuracy of your calculation. From the excerpt you provide, I see that the energy changed by 2.71E-05 Hartree, or nearly one meV, in the last iteration before convergence. That's rather a lot, and hence I'd advise you to decrease EPS_SCF to, say, 1.0E-6.</div><div><br></div><div>However, in the end, the two values of total energy differ by less than 0.08 meV. For most intents and purposes, that's identical, and you can choose whichever.</div><div><br></div><div>Yours sincerely,</div><div>Patrick Gono</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, 12 Jan 2020 at 08:09, e khodayar <<a href="mailto:e.khod...@gmail.com" target="_blank">e.khod...@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><span style="font-size:small">Hi,</span><br><div><span style="color:rgb(0,0,0);font-family:"Times New Roman";font-size:medium">In the main output file</span> geometry optimization, <span style="color:rgb(0,0,0);font-family:"Times New Roman";font-size:medium">at the end of geometry optimisation step, I have </span><span style="color:rgb(0,0,0);font-family:"Times New Roman";font-size:medium">the following information:</span><span style="font-size:small"><br></span></div><div><div><font color="#000000" face="Times New Roman" size="3"> -------- Informations at step = 53 ------------</font></div><div><font color="#000000" face="Times New Roman" size="3"> Optimization Method = CG</font></div><div><font color="#000000" face="Times New Roman" size="3" style="background-color:rgb(255,255,0)"> Total Energy = -94.0112708604</font></div><div><font color="#000000" face="Times New Roman" size="3"> Real energy change = -0.0000036472</font></div><div><font color="#000000" face="Times New Roman" size="3"> Decrease in energy = YES</font></div><div><font color="#000000" face="Times New Roman" size="3"> Used time = 343.953</font></div><div><font color="#000000" face="Times New Roman" size="3"><br></font></div><div><font color="#000000" face="Times New Roman" size="3"> Convergence check :</font></div><div><font color="#000000" face="Times New Roman" size="3"> Max. step size = 0.0006162087</font></div><div><font color="#000000" face="Times New Roman" size="3"> Conv. limit for step size = 0.0010000000</font></div><div><font color="#000000" face="Times New Roman" size="3"> Convergence in step size = YES</font></div><div><font color="#000000" face="Times New Roman" size="3"> RMS step size = 0.0002737159</font></div><div><font color="#000000" face="Times New Roman" size="3"> Conv. limit for RMS step = 0.0010000000</font></div><div><font color="#000000" face="Times New Roman" size="3"> Convergence in RMS step = YES</font></div><div><font color="#000000" face="Times New Roman" size="3"> Max. gradient = 0.0005305121</font></div><div><font color="#000000" face="Times New Roman" size="3"> Conv. limit for gradients = 0.0010000000</font></div><div><font color="#000000" face="Times New Roman" size="3"> Conv. in gradients = YES</font></div><div><font color="#000000" face="Times New Roman" size="3"> RMS gradient = 0.0002356500</font></div><div><font color="#000000" face="Times New Roman" size="3"> Conv. limit for RMS grad. = 0.0010000000</font></div><div><font color="#000000" face="Times New Roman" size="3"> Conv. in RMS gradients = YES</font></div><div><font color="#000000" face="Times New Roman" size="3"> ---------------------------------------------------</font></div><div><font color="#000000" face="Times New Roman" size="3"><br></font></div><div><font color="#000000" face="Times New Roman" size="3"> *******************************************************************************</font></div><div><font color="#000000" face="Times New Roman" size="3"> *** GEOMETRY OPTIMIZATION COMPLETED ***</font></div><div><font color="#000000" face="Times New Roman" size="3"> *******************************************************************************</font></div><div><font color="#000000" face="Times New Roman" size="3"><br></font></div><div><font color="#000000" face="Times New Roman" size="3"> Reevaluating energy at the minimum</font></div><div><font color="#000000" face="Times New Roman" size="3"><br></font></div><div><font color="#000000" face="Times New Roman" size="3"> Number of electrons: 68</font></div><div><font color="#000000" face="Times New Roman" size="3"> Number of occupied orbitals: 34</font></div><div><font color="#000000" face="Times New Roman" size="3"> Number of molecular orbitals: 34</font></div><div><font color="#000000" face="Times New Roman" size="3"><br></font></div><div><font color="#000000" face="Times New Roman" size="3"> Number of orbital functions: 232</font></div><div><font color="#000000" face="Times New Roman" size="3"> Number of independent orbital functions: 232</font></div><div><font color="#000000" face="Times New Roman" size="3"><br></font></div><div><font color="#000000" face="Times New Roman" size="3"> Extrapolation method: ASPC</font></div><div><font color="#000000" face="Times New Roman" size="3"><br></font></div><div><font color="#000000" face="Times New Roman" size="3"><br></font></div><div><font color="#000000" face="Times New Roman" size="3"> SCF WAVEFUNCTION OPTIMIZATION</font></div><div><font color="#000000" face="Times New Roman" size="3"><br></font></div><div><font color="#000000" face="Times New Roman" size="3"> Step Update method Time Convergence Total energy Change</font></div><div><font color="#000000" face="Times New Roman" size="3"> ------------------------------------------------------------------------------</font></div><div><font color="#000000" face="Times New Roman" size="3"> 1 Pulay/Diag. 0.50E+00 3.7 0.00001087 -94.0112708604 -9.40E+01</font></div><div><font color="#000000" face="Times New Roman" size="3"> 2 Pulay/Diag. 0.50E+00 4.7 0.00001294 -94.0112951749 -2.43E-05</font></div><div><font color="#000000" face="Times New Roman" size="3"> 3 Pulay/Diag. 0.50E+00 4.3 0.00000783 -94.0112680643 2.71E-05</font></div><div><font color="#000000" face="Times New Roman" size="3"><br></font></div><div><font color="#000000" face="Times New Roman" size="3"> *** SCF run converged in 3 steps ***</font></div><div><font color="#000000" face="Times New Roman" size="3"><br></font></div><div><font color="#000000" face="Times New Roman" size="3"><br></font></div><div><font color="#000000" face="Times New Roman" size="3"> Electronic density on regular grids: -67.9999999983 0.0000000017</font></div><div><font color="#000000" face="Times New Roman" size="3"> Core density on regular grids: 68.0000000001 0.0000000001</font></div><div><font color="#000000" face="Times New Roman" size="3"> Total charge density on r-space grids: 0.0000000017</font></div><div><font color="#000000" face="Times New Roman" size="3"> Total charge density g-space grids: 0.0000000017</font></div><div><font color="#000000" face="Times New Roman" size="3"><br></font></div><div><font color="#000000" face="Times New Roman" size="3"> Overlap energy of the core charge distribution: 0.00000565490189</font></div><div><font color="#000000" face="Times New Roman" size="3"> Self energy of the core charge distribution: -223.82607238936723</font></div><div><font color="#000000" face="Times New Roman" size="3"> Core Hamiltonian energy: 69.37181702794591</font></div><div><font color="#000000" face="Times New Roman" size="3"> Hartree energy: 89.86177654697752</font></div><div><font color="#000000" face="Times New Roman" size="3"> Exchange-correlation energy: -29.41879490480421</font></div><div><font color="#000000" face="Times New Roman" size="3"><br></font></div><div><font color="#000000" style="background-color:rgb(255,255,0)"><font face="Times New Roman" size="3"> </font><font face="Times New Roman" size="3">Total energy: -94.01126806434611</font></font></div><div><font color="#000000" face="Times New Roman" size="3"><br></font></div><div><font color="#000000" face="Times New Roman" size="3"> ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -94.011268064346112</font></div><div><font color="#000000" face="Times New Roman" size="3"><br></font></div><div><span style="color:rgb(0,0,0);font-family:Arial,Helvetica,sans-serif;font-size:13px">Now my question is:</span></div><div><span style="color:rgb(0,0,0);font-family:Arial,Helvetica,sans-serif;font-size:13px">Which total energy i have to use for my data?and why this two total energy that with yellow color is d</span><font color="#000000">ifferent?</font><br></div></div><div>Thanks in advance.<font color="#000000"><br></font></div><div><font color="#000000"><br></font></div></div>
<p></p>
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/0996cb3f-acf7-4e34-b948-7be5936dba1f%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/0996cb3f-acf7-4e34-b948-7be5936dba1f%40googlegroups.com</a>.<br>
</blockquote></div>
<p></p>
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAPL1_v8JenHMs7NHOX1%3DWbZQwHwF7-%2Bu2Pbt1mYBbfAQ3gT3Gw%40mail.gmail.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/CAPL1_v8JenHMs7NHOX1%3DWbZQwHwF7-%2Bu2Pbt1mYBbfAQ3gT3Gw%40mail.gmail.com</a>.<br>
</blockquote></div>