<div dir="ltr">Dear Ivan, Victor and Marcella, <div><br></div><div>Thank you so much for your help. </div><div><br></div><div>I have been reading the bibliography you suggested and using your inputs and I almost optimised the heme group. I got to a point where the optimisation does not seem to converge but it is really close. </div><div><br></div><div><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">--------<span class="gmail-Apple-converted-space"> </span>Informations at step = <span class="gmail-Apple-converted-space"> </span>114 ------------</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Optimization Method<span class="gmail-Apple-converted-space"> </span>= <span class="gmail-Apple-converted-space"> </span>SD</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Total Energy <span class="gmail-Apple-converted-space"> </span>=<span class="gmail-Apple-converted-space"> </span>-517.5364879764</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Real energy change <span class="gmail-Apple-converted-space"> </span>= <span class="gmail-Apple-converted-space"> </span>0.0000000050</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Decrease in energy <span class="gmail-Apple-converted-space"> </span>= <span class="gmail-Apple-converted-space"> </span>NO</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Used time<span class="gmail-Apple-converted-space"> </span>=<span class="gmail-Apple-converted-space"> </span>460.425</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);min-height:13px"><span style="font-variant-ligatures:no-common-ligatures"></span><br></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Convergence check :</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Max. step size <span class="gmail-Apple-converted-space"> </span>= <span class="gmail-Apple-converted-space"> </span>0.0000000000</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Conv. limit for step size<span class="gmail-Apple-converted-space"> </span>= <span class="gmail-Apple-converted-space"> </span>0.0010000000</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Convergence in step size <span class="gmail-Apple-converted-space"> </span>=<span class="gmail-Apple-converted-space"> </span></span><span style="font-variant-ligatures:no-common-ligatures;background-color:rgb(224,228,9)">YES</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>RMS step size<span class="gmail-Apple-converted-space"> </span>= <span class="gmail-Apple-converted-space"> </span>0.0000000000</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Conv. limit for RMS step <span class="gmail-Apple-converted-space"> </span>= <span class="gmail-Apple-converted-space"> </span>0.0010000000</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Convergence in RMS step<span class="gmail-Apple-converted-space"> </span>=<span class="gmail-Apple-converted-space"> </span></span><span style="font-variant-ligatures:no-common-ligatures;background-color:rgb(224,228,9)">YES</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Max. gradient<span class="gmail-Apple-converted-space"> </span>= <span class="gmail-Apple-converted-space"> </span>0.0012644252</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Conv. limit for gradients<span class="gmail-Apple-converted-space"> </span>= <span class="gmail-Apple-converted-space"> </span>0.0010000000</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Conv. for gradients<span class="gmail-Apple-converted-space"> </span>= <span class="gmail-Apple-converted-space"> </span>NO</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>RMS gradient <span class="gmail-Apple-converted-space"> </span>= <span class="gmail-Apple-converted-space"> </span>0.0002927175</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Conv. limit for RMS grad.<span class="gmail-Apple-converted-space"> </span>= <span class="gmail-Apple-converted-space"> </span>0.0010000000</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>Conv. in RMS gradients <span class="gmail-Apple-converted-space"> </span>=<span class="gmail-Apple-converted-space"> </span></span><span style="font-variant-ligatures:no-common-ligatures;background-color:rgb(224,228,9)">YES</span></p>
<p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>---------------------------------------------------</span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><br></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">I have been playing a bit with the parameters but I haven't managed to make it work. </span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><br></span></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)">Thanks again, </p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)">Kind regards, </p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><br></p><p style="margin:0px;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)">Sara</p></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El jue, 15 sept 2022 a las 17:57, Marcella Iannuzzi (<<a href="mailto:marci.akira@gmail.com">marci.akira@gmail.com</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">Dear Sara<div><br></div><div>Before starting a geometry optimisation it is good practice to get a converged single point calculation. </div><div>If I understand correctly from the output you posted, this was not the case. </div><div><br></div><div>I played a bit with the settings for your system and I got a converged single point calculation at the PBE level of theory </div><div>for the molecule charged --2 (see the attached file).</div><div><br></div><div>Kind regards</div><div>Marcella</div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, September 13, 2022 at 4:00:25 PM UTC+2 <a href="mailto:salla...@gmail.com" target="_blank">salla...@gmail.com</a> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">Dear all, <div><br></div><div>I am having problems on a geometry optimisation of a heme group with Fe2+ using PBE in CP2K. My calculations run and the total energy is decreasing. However, the job has been running for almost a month and it doesn't seem to converge. </div><div><br></div><div>I am not an expert on CP2K nor on transition metals. I would really appreciate it if you could help me sort this out. I have attached the restart file and the output for clarity. </div><div><br></div><div>Thanks a lot, </div><div><br></div><div><br></div><div>Sara</div></blockquote></div>
<p></p>
-- <br>
You received this message because you are subscribed to a topic in the Google Groups "cp2k" group.<br>
To unsubscribe from this topic, visit <a href="https://groups.google.com/d/topic/cp2k/C7w_GP_CIxw/unsubscribe" target="_blank">https://groups.google.com/d/topic/cp2k/C7w_GP_CIxw/unsubscribe</a>.<br>
To unsubscribe from this group and all its topics, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com" target="_blank">cp2k+unsubscribe@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/bf21e618-f8d1-43d4-9e58-463e1b671947n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/bf21e618-f8d1-43d4-9e58-463e1b671947n%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:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAPHo29Y3BNicWoKMdkxERY4ppY_5qP-avWUxA-w3ruitG5H2gg%40mail.gmail.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/CAPHo29Y3BNicWoKMdkxERY4ppY_5qP-avWUxA-w3ruitG5H2gg%40mail.gmail.com</a>.<br />