<div dir="ltr"><div>Thank you Matt, your answer was helpful and it worked!</div><div>The BS section allowed me to set up an excited state of each atom for the atomic guess, which is then rescaled to the number of electrons in each spin channel. However, I found that the resultant wavefunction can converge to some local minima. For example, while testing an O2 binding energy curve, there were some points that initially converged to ~1eV higher than the neighboring ones, however restarted form the wavefunctions of the distances adjacent to them they gave me energies consistent with the ones around them. <br></div><div><br></div><div>Do you have any idea how to make my calculations, and perhaps the initial guess more robust for molecular systems?</div><div><br></div><div>Thanks,</div><div>Tamas<br></div><br>On Sunday, 29 December 2019 16:58:15 UTC, Matt W  wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Yes, you need a broken symmetry guess. <div><br></div><div>You can use the BS section in the KIND section (many threads on how to do this it is non intuitive) or you can use the </div><div><br></div><div><a href="https://manual.cp2k.org/cp2k-6_1-branch/CP2K_INPUT/FORCE_EVAL/DFT/PRINT/WFN_MIX.html" target="_blank" rel="nofollow" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fmanual.cp2k.org%2Fcp2k-6_1-branch%2FCP2K_INPUT%2FFORCE_EVAL%2FDFT%2FPRINT%2FWFN_MIX.html\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNGcSNlzzP-9lzOgo18eVZmpNK_z8A';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fmanual.cp2k.org%2Fcp2k-6_1-branch%2FCP2K_INPUT%2FFORCE_EVAL%2FDFT%2FPRINT%2FWFN_MIX.html\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNGcSNlzzP-9lzOgo18eVZmpNK_z8A';return true;">https://manual.cp2k.org/cp2k-<wbr>6_1-branch/CP2K_INPUT/FORCE_<wbr>EVAL/DFT/PRINT/WFN_MIX.html</a></div><div><br></div><div>section to set up your calculation as you are attempting here. Again, it is not all that easy to use...</div><div><br></div><div>Matt<br><br>On Sunday, December 29, 2019 at 2:38:45 PM UTC, Tamas K. Stenczel wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><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 style="margin:4px 0px;padding:8px;font-size:12px;line-height:1.50001;white-space:pre-wrap;word-break:normal;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);color:rgb(29,28,29);font-style:normal;font-weight:400;letter-spacing:normal;text-align:left;text-indent:0px;text-transform:none;word-spacing:0px">******************************<wbr>******************************<wbr>*******************<br> *   ___                                                                       *<br> *  /   \                                                                      *<br> * [ABORT]                                                                     *<br> *  \___/    Number of occupied MOs on restart unit larger than allocated MOs. *<br> *    |                                                                        *<br> *  O/|                                                                        *<br> * /| |                                                                        *<br> * / \                                                          qs_mo_io.F:818 *<br> ******************************<wbr>******************************<wbr>*******************</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></blockquote></div></div></blockquote></div>