<p style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);">Hello,</p><p style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);">I want to perform single-point energy and force calculations on a condensed dodecane system (684 atoms, 18 chains) for geometries including a few strained/broken C–C bonds. While the SCF converges for geometries with single broken bond, I encounter convergence issues when multiple bonds (8) are broken. Input and output files for the multiple-bond case are attached.</p><p style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);">I am using max_SCF 20 for inner SCF and 70 for outer SCF. After ~3 outer SCF iterations, the energy mostly stabilizes but then drops by 20 Hartrees at iteration 33. Are the any parameters I should change (so far I've tried changing the energy_gap and linesearch values)? Apologies if this is a basic issue—I am new to CP2K and DFT. For context, the energies and forces will be used to train a machine learning potential for polyethylene to capture C–C bond breaking, which is why I use a relatively strict convergence criteria and high grid cutoff.</p><p style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);">Thank you!</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 visit <a href="https://groups.google.com/d/msgid/cp2k/f5696b77-8cc9-4a66-8cb6-fb14f27923cbn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/f5696b77-8cc9-4a66-8cb6-fb14f27923cbn%40googlegroups.com</a>.<br />