<div dir="ltr"><div>Dear cp2k developers and users,</div><div><br></div><div>I am working modelling the combustion of small hydrocarbons and have issue with spin. I am using CP2K 5.1 right now.<br></div><div><br></div><div>I would like to calculate my system in UKS singlet state, with broken symmetry of the alpha-beta channels.  I have understood that this is the expected behavior of a large system with many oxygen molecules and the result should be better than the singlet with RKS or with spin-symmetric UKS (these two give me the same results actually).</div><div><br></div><div>Is this possible with the OT methods, how to do it? I have been using OT%DIIS for its performance with SCF convergence. Diag and mixing methods were not as reliable as this so far, but can perform relaxation of multiplicity, which has given me asymetric spins but I cannot use the wavefunctions of that to restart OT from, due to the error:</div><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div><pre class="c-mrkdwn__pre" data-stringify-type="pre" style="box-sizing: inherit; margin: 4px 0px; padding: 8px; --saf-0:rgba(var(--sk_foreground_low,29,28,29),0.13); font-size: 12px; line-height: 1.50001; font-variant-ligatures: none; white-space: pre-wrap; overflow-wrap: break-word; word-break: normal; tab-size: 4; font-family: Monaco, Menlo, Consolas, "Courier New", monospace !important; border: 1px solid var(--saf-0); border-radius: 4px; background: rgba(var(--sk_foreground_min,29,28,29),0.04); counter-reset: list-0 0 list-1 0 list-2 0 list-3 0 list-4 0 list-5 0 list-6 0 list-7 0 list-8 0 list-9 0; color: rgb(29, 28, 29); font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: 2; text-align: left; text-indent: 0px; text-transform: none; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration-style: initial; text-decoration-color: initial;">*******************************************************************************<br style="box-sizing: inherit;"> *   ___                                                                       *<br style="box-sizing: inherit;"> *  /   \                                                                      *<br style="box-sizing: inherit;"> * [ABORT]                                                                     *<br style="box-sizing: inherit;"> *  \___/    Number of occupied MOs on restart unit larger than allocated MOs. *<br style="box-sizing: inherit;"> *    |                                                                        *<br style="box-sizing: inherit;"> *  O/|                                                                        *<br style="box-sizing: inherit;"> * /| |                                                                        *<br style="box-sizing: inherit;"> * / \                                                          qs_mo_io.F:818 *<br style="box-sizing: inherit;"> *******************************************************************************</pre></div></blockquote><div><br></div><div><br></div><div>some of the input. I have attached an input file of an O2 molecule for example.</div><div>&DFT<br>    POTENTIAL_FILE_NAME GTH_POTENTIALS<br>    UKS TRUE<br>    MULTIPLICITY 1<br>    BASIS_SET_FILE_NAME COMBINED_BASIS_TKS32<br>    &SCF<br>      MAX_SCF 25<br>      EPS_SCF 0.0001<br>      SCF_GUESS RESTART<br>      &OT<br>        MINIMIZER DIIS<br>      &END OT<br>      &OUTER_SCF<br>        EPS_SCF 0.0001<br>        MAX_SCF 4<br>      &END OUTER_SCF<br>    &END SCF</div><div><br></div><div><br></div><div>Thank you very much,</div><div>Tamas K Stenczel<br></div></div>