Dear CP2K developers and users,<div><br></div><div>Here I report a bug of CP2K (6.1, 7.1, and 8.2), or maybe it's not a bug. Any responses and answers to the conversation are appreciated.</div><div><br></div><div>Descriptions: </div><div>a) test models: ceria (stoichiometric CeO2)</div><div>b) occurring in Geometry optimation or molecular dynamic simulation</div><div>c) bug or error: in the first few steps, QS-SCF gives the correct wavefunction (no spin population value on Ce atoms that are zero for stoichiometric ceria ). But in the further GEO-OPT or MD, QS-SCF suddenly needs more SCF-steps to reach convergence and give a completely wrong spin on Ce and O atoms (see output file attached). Eventually, the job will die due to the SCF-convergence problem or converger to a wrong state.</div><div>d) such BUG only exist in the Orthorhombic system, for the Rhombohedron there is no such BUG</div><div>e) currently, my solution for this bug in the orthorhombic system is to restart the calculation without reading wavefunction again and again. The QS-OT can obtain the correct wavefunction when it starts SCF from atomic guess.</div><div><br></div><div>I have been trying to solve this problem by running hundreds of tests in which I tried to change so many parameters but I never get it solved.  By the way, such BUG only exists in OT methods, it's completely OK for the diagonalization method.</div><div><br></div><div>Best regards,</div><div>YONG</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/c9e8d155-2ddc-45cd-8291-5a19c040b07dn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/c9e8d155-2ddc-45cd-8291-5a19c040b07dn%40googlegroups.com</a>.<br />