<div dir="ltr">Hi Olivia & Marcella, <div><br></div><div>I allow myself to jump back on this problem/answer, and ask the CP2K team/users:</div><div><br></div><div>For a MD calculation (or maybe other, geometry for example here), is there a way to stop it if the SCF hasn't converged (after the inner/outer SCF) ? </div><div><br></div><div>I try to find in the manual but no infos (or i looked at the wrong place ?)</div><div><br></div><div>If one SCF didn't converge, forces are still calculated and MD goes on to next step. So I kill it manually, but making it automatically save CPU usage and tell you right away: "No convergence, please check".</div><div><br></div><div>Thanks</div><div><br></div><div>Best regards,</div><div><br></div><div>Rolf</div><div><br></div><div>On Thursday, July 16, 2015 at 1:59:54 PM UTC+2, Marcella Iannuzzi wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><div><br></div><div><br></div>Dear Olivia,<div><br></div><div>The SCF in your run is never converging. Forces  calculated from a not converged SCF are affected by large errors,  which explains the behaviour of the geometry optimisation. </div><div>One shouldn't  carry on a geometry optimisation run when the SCF does not converge.</div><div>Allowing a larger number of SCF iterations might help (see OUTER_SCF / MAX_SCF).</div><div><br></div><div>Kind regards</div><div>Marcella <br></div><div><br>On Thursday, July 16, 2015 at 12:41:44 PM UTC+2, Olivia Lynes wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi all,<div><br></div><div>I'm having a problem with a few of my geometry optimisations not converging after several thousand steps and having looked at the energy changes through the main output file it's not decreasing but is fluctuating a lot in a certain range. I'm at a loss to why as I've got similar calculations completed with similar input and just a different box size and they have converged without this fluctuation in a few hundred steps.</div><div><br></div><div>I've attached the input, coordination files and the first hundred or so steps of the output file.</div><div><br></div><div>These are geometry optimisations of an Mg 2+ ion in a box of 64 waters, with a cell size of 13.41 angstroms, so the overall system is charged 2. I've used the GTH DZVP basis sets and potentials. </div><div><br></div><div>Having done a cell optimisation on just the 64 waters, we introduced the Mg ion and charged the system. In lieu of being able to do cell optimisations on a charged system, then geometry optimisations are being done at varying box sizes by 10ths of an angstrom as a sort of manual cell optimisation to try to find both a minimum energy structure and a minimum energy box size. </div><div><br></div><div>I've done similar calculations using Ca and Sr and haven't had this issue with any of my jobs, they all converge in a few hundred steps. </div><div><br></div><div>Is this potentially a problem using charged systems? Is it necessary for me to put a counter ion in which would let me do a standard cell optimisation? </div><div><br></div><div>Or is this approach completely wrong and we would be better off doing MD? </div><div><br></div><div>The aim is to do ab initio MD on these systems to look at the coordination of the ion in a bulk of water, so reaching a minimum energy structure and box size isn't essential but it would be nice to have the regular DFT data to back it up. </div><div><br></div><div>Thanks in advance for any insights.</div><div><br></div><div>Olivia</div><div><br></div><div><br></div><div><br></div></div></blockquote></div></div></blockquote></div></div>