Dear cp2k users<div><br></div><div>Hi, I am wondering whether I can specify negative multiplicity for a particular system. I have put multiplicity as -2 for an anion molecule for a DFT calculation and I got the following error : </div><div><br></div><div><div> <b>*** WARNING in qs_environment.F:1295 :: More beta than alpha MOs          ***</b></div><div><b> *** requested. The number of beta MOs will be reduced to the number alpha ***</b></div><div><b> *** MOs.                                                                  ***</b></div><div><b><br></b></div><div><b><br></b></div><div><b> *******************************************************************************</b></div><div><b> *   ___                                                                       *</b></div><div><b> *  /   \                                                                      *</b></div><div><b> * [ABORT]                                                                     *</b></div><div><b> *  \___/                             CPASSERT failed                          *</b></div><div><b> *    |                                                                        *</b></div><div><b> *  O/|                                                                        *</b></div><div><b> * /| |                                                                        *</b></div><div><b> * / \                                                   qs_environment.F:1300 *</b></div><div><b> *******************************************************************************</b></div><div><b><br></b></div><div><b><br></b></div><div><b> ===== Routine Calling Stack =====</b></div><div><b><br></b></div><div><b>            2 qs_init_subsys</b></div><div><b>            1 CP2K</b></div></div><div><br></div><div>Based on the last discussion by Pavan in April last year, it seems like cp2k does not support a way to generate configuration where the unpaired electron is in the spin down state (beta). But I would like to ask once again if someone has come up with a way to specify beta spin state for the unpaired electron. I understand that alpha and beta states should be degenerate in the absence of magnetic field but my system needs to treat alpha and beta spin states separately.</div><div><br></div><div><b>Background of my system (just in case)</b></div><div><br></div><div>I am studying a system where an electron is transferred from an organic molecule to a oxygen molecule. Let A be the organic molecule (IDFBR) and B be the oxygen molecule and the electron transfer can be described by the following equation:</div><div><br></div><div>[A]- + [B]  --> [A] + [B]-</div><div><br></div><div>Since the ground state of oxygen molecule is triplet and anion organic molecule is doublet, I believe that the electron can only transferred from the anion organic molecule if the spin states of the electrons in A- and B are correctly aligned. (alpha electron at the SOMO of [A]- if electrons at the SOMOs of B are beta or  beta electrons at the SOMO of [A]- if electrons at the SOMOs of B are alpha).</div><div><br></div><div><br><img alt="electron transfer.png" width="128px" height="213px" data-atf="true" src="cid:b90bac8f-75af-4964-a0b3-0693ff8c72f3">                                                                <img alt="electron transfer2.png" width="128px" height="213px" data-atf="true" src="cid:fa97f6ef-e6ce-4c38-8ca0-a6f9b2d9b0f8"></div><div>Electron transfer not possible                                               Electron transfer possible</div><div><br></div><div><br></div><div>Therefore, I need to study the particular system where unpaired electron of A has opposite spin state as B.</div><div><br></div><div>Thank you so much for your time and I would really appreciate any comments / help on this.</div><div><br></div><div>kind regards</div><div><br></div><div>jae</div>