Hi Matthias,<div><br /></div><div>I did have a look at the cutoff and EPS_DEFAULT at the start of my project but that was for a different issue. So I haven't tried running consecutive optimisations with these parameters but I will give that a go and get back to you.</div><div><br /></div><div>Thanks for your help,</div><div><br /></div><div>Daniel</div><div><br /></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, October 8, 2024 at 2:35:01 PM UTC+1 Krack Matthias wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">





<div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word">
<div>
<p class="MsoNormal"><span lang="DE-CH" style="font-size:11.0pt">Hi Daniel<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="DE-CH" style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">Did you try to decrease EPS_DEFAULT and to increase the PW cutoff? These input parameters mainly determine the accuracy of the atomic forces rather than the convergence
 thresholds for the force calculation like MAX_FORCE or MAX_DR which can also show convergence for poor forces.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">HTH<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">Matthias<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><u></u> <u></u></span></p>
<div>
<div>
<div></div></div></div></div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div><div><div><div>
<div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal" style="margin-right:0cm;margin-bottom:12.0pt;margin-left:36.0pt">
<b><span style="color:black">From: </span></b><span style="color:black">'daniel Storm' via cp2k <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>><br>
<b>Date: </b>Tuesday, 8 October 2024 at 14:50<br>
<b>To: </b>cp2k <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>><br>
<b>Subject: </b>[CP2K:20756] Inconsistent geometry optimisations<u></u><u></u></span></p>
</div>
<p style="margin-right:0cm;margin-bottom:8.0pt;margin-left:36.0pt;text-align:justify;line-height:106%">
<span style="font-size:13.0pt;line-height:106%;color:black">Hi everyone,</span><u></u><u></u></p>
<p style="margin-right:0cm;margin-bottom:8.0pt;margin-left:36.0pt;text-align:justify;line-height:106%">
<span style="font-size:13.0pt;line-height:106%;color:black">Just to provide some context, in our group we use CP2k to model molecular reactivity in the solid state, so slightly different than “standard” solid state chemistry. We run optimisations for minima
 and transition states (TS), whose character we confirm by running phonon calculations.</span><u></u><u></u></p>
<p style="margin-right:0cm;margin-bottom:8.0pt;margin-left:36.0pt;text-align:justify;line-height:106%">
<span style="font-size:13.0pt;line-height:106%;color:black">Now, to the issue. I am working with a crystal structure of an iridium organometallic salt, [Ir(PONOP)(H)Me][BAr<sup>F</sup><sub>4</sub>], however, I am encountering issues with achieving consistent
 optimisations. Despite starting from the same experimental structure, I am getting very different SCF energies (up to 3 kcal/mol) and I have observed minor differences in the resulting geometries. I think that those geometries are close enough to not show
 that big of an SCF difference, but my main worry is the fact that the starting point is always the same.</span><u></u><u></u></p>
<p style="margin-right:0cm;margin-bottom:8.0pt;margin-left:36.0pt;text-align:justify;line-height:106%">
<span style="font-size:13.0pt;line-height:106%;color:black">The convergence criteria that I am currently using is MAX_FORCE 1.0 x 10<sup>-4</sup>. I have attempted to improve this by tightening the criteria (MAX_FORCE, MAX_DR, RMS_DR and RMS_FORCE) to 1.0 x
 10<sup>-6</sup>, but this didn’t fix the issue. All the optimisations are done with the PBE-D3 functional and DZVP-MOLOPT-SR-GTH basis set.</span><u></u><u></u></p>
<p style="margin-right:0cm;margin-bottom:8.0pt;margin-left:36.0pt;text-align:justify;line-height:106%">
<span style="font-size:13.0pt;line-height:106%;color:black">So, to sum up, I run optimisations with the same starting geometry that converge to different structures, and I am not sure why or which criterion to use to select one structure over the other, etc.
</span><u></u><u></u></p>
<p style="margin-right:0cm;margin-bottom:8.0pt;margin-left:36.0pt;text-align:justify;line-height:106%">
<span style="font-size:13.0pt;line-height:106%;color:black">I have worked with other organometallic salts in the past, and this is the first time I find this issue. I am using CP2k version 2023.2 in the HPC cluster Archer-2.</span><u></u><u></u></p>
<p style="margin-right:0cm;margin-bottom:8.0pt;margin-left:36.0pt;text-align:justify;line-height:106%">
<span style="font-size:13.0pt;line-height:106%;color:black">Has anyone seen this before? Any advice on how to proceed?</span><u></u><u></u></p>
<p style="margin-right:0cm;margin-bottom:8.0pt;margin-left:36.0pt;text-align:justify;line-height:106%">
<span style="font-size:13.0pt;line-height:106%;color:black">Thanks in advance,</span><u></u><u></u></p>
<p style="margin-right:0cm;margin-bottom:8.0pt;margin-left:36.0pt;text-align:justify;line-height:106%">
<span style="font-size:13.0pt;line-height:106%;color:black">Daniel.</span><u></u><u></u></p>
</div></div></div></div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div><div><div><div><p class="MsoNormal" style="margin-left:36.0pt">-- <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 data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/d50cb0ef-42e9-465d-81c7-2fdacdab6d41n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/d50cb0ef-42e9-465d-81c7-2fdacdab6d41n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1728551269960000&usg=AOvVaw3449S2CFrjjQc7XxOciAHa">
https://groups.google.com/d/msgid/cp2k/d50cb0ef-42e9-465d-81c7-2fdacdab6d41n%40googlegroups.com</a>.<u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>

</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/07149c81-16d2-414a-9816-aeafe65595afn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/07149c81-16d2-414a-9816-aeafe65595afn%40googlegroups.com</a>.<br />