<p style="line-height: 107%; margin-top: 10pt; margin-bottom: 8pt; margin-left: 0in; text-indent: 0in; text-align: justify; direction: ltr; unicode-bidi: embed; word-break: normal;"><span style="font-size: 13pt; font-family: Aptos; color: black;">Hi
everyone</span><span style="font-size: 13pt; font-family: Aptos; color: black;">,</span></p>

<p style="line-height: 107%; margin-top: 10pt; margin-bottom: 8pt; margin-left: 0in; text-indent: 0in; text-align: justify; direction: ltr; unicode-bidi: embed; word-break: normal;"><span style="font-size: 13pt; font-family: Aptos; 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 </span><span style="font-size: 13pt; font-family: Aptos; color: black;">chemistry</span><span style="font-size: 13pt; font-family: Aptos; color: black;">.</span><span style="font-size: 13pt; font-family: Aptos; color: black;"> </span><span style="font-size: 13pt; font-family: Aptos; color: black;">We
run optimisations for minima and transition states (TS), whose character we
confirm by running phonon calculations.</span></p>

<p style="line-height: 107%; margin-top: 10pt; margin-bottom: 8pt; margin-left: 0in; text-indent: 0in; text-align: justify; direction: ltr; unicode-bidi: embed; word-break: normal;"><span style="font-size: 13pt; font-family: Aptos; color: black;">Now,
to the issue. I am working with a crystal structure of an iridium
organometallic salt, [</span><span style="font-size: 13pt; font-family: Aptos; color: black;">Ir</span><span style="font-size: 13pt; font-family: Aptos; color: black;">(PONOP)(H)Me][BAr</span><span style="font-size: 13pt; font-family: Aptos; color: black; vertical-align: super;">F</span><span style="font-size: 13pt; font-family: Aptos; color: black; vertical-align: sub;">4</span><span style="font-size: 13pt; font-family: Aptos; color: black;">],
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></p>

<p style="line-height: 107%; margin-top: 10pt; margin-bottom: 8pt; margin-left: 0in; text-indent: 0in; text-align: justify; direction: ltr; unicode-bidi: embed; word-break: normal;"><span style="font-size: 13pt; font-family: Aptos; color: black;">The
convergence criteria that I am currently using is MAX_FORCE 1.0 x 10</span><span style="font-size: 13pt; font-family: Aptos; color: black; vertical-align: super;">-4</span><span style="font-size: 13pt; font-family: Aptos; color: black;">. I have attempted to improve this by
tightening the criteria (MAX_FORCE, MAX_DR, RMS_DR and RMS_FORCE) to 1.0 x 10</span><span style="font-size: 13pt; font-family: Aptos; color: black; vertical-align: super;">-6</span><span style="font-size: 13pt; font-family: Aptos; color: black;">, 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></p>

<p style="line-height: 107%; margin-top: 10pt; margin-bottom: 8pt; margin-left: 0in; text-indent: 0in; text-align: justify; direction: ltr; unicode-bidi: embed; word-break: normal;"><span style="font-size: 13pt; font-family: Aptos; 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></p>

<p style="line-height: 107%; margin-top: 10pt; margin-bottom: 8pt; margin-left: 0in; text-indent: 0in; text-align: justify; direction: ltr; unicode-bidi: embed; word-break: normal;"><span style="font-size: 13pt; font-family: Aptos; 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></p>

<p style="line-height: 107%; margin-top: 10pt; margin-bottom: 8pt; margin-left: 0in; text-indent: 0in; text-align: justify; direction: ltr; unicode-bidi: embed; word-break: normal;"><span style="font-size: 13pt; font-family: Aptos; color: black;">Has
anyone seen this before? Any advice on how to proceed?</span></p>

<p style="line-height: 107%; margin-top: 10pt; margin-bottom: 8pt; margin-left: 0in; text-indent: 0in; text-align: justify; direction: ltr; unicode-bidi: embed; word-break: normal;"><span style="font-size: 13pt; font-family: Aptos; color: black;">Thanks
in advance,</span></p>

<p style="line-height: 107%; margin-top: 10pt; margin-bottom: 8pt; margin-left: 0in; text-indent: 0in; text-align: justify; direction: ltr; unicode-bidi: embed; word-break: normal;"><span style="font-size: 13pt; font-family: Aptos; color: black;">Daniel.</span></p>

<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/d50cb0ef-42e9-465d-81c7-2fdacdab6d41n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/d50cb0ef-42e9-465d-81c7-2fdacdab6d41n%40googlegroups.com</a>.<br />