<div dir="ltr">Dear Du,<div><br></div><div>it's a matter of trial and error and finding the right combination of total charge and spin. To me it's not clear what OH- is. OH* is a radical, hydronium, on the other hand, has a closed shell. SO4- is also a strange beast. If you set total spin to 2 and charge to -2, you'll likely get OH* and SO4(2-). There is an option to play with the guess:</div><div><br></div><div>https://manual.cp2k.org/cp2k-6_1-branch/CP2K_INPUT/FORCE_EVAL/SUBSYS/KIND/BS.html</div><div><br></div><div>Yours,</div><div><br></div><div>Vladimir<br><br>среда, 15 апреля 2020 г., 3:39:58 UTC+2 пользователь Du написал:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear Vladimir,<div><br></div><div>Thank you! </div><div><br></div><div>Here is the thing, I want to simulate the SO4- radical and the OH- radical, in a bulk of 40 water molecules. I don't know whether cp2k recognize them as radicals or just ions<br>? As you mentioned, I've already set the UKS T in the input profile. </div><div><br></div><div>Yours,</div><div><br></div><div>Du</div><div>在 2020年4月14日星期二 UTC+8下午6:54:02，Vladimir Rybkin写道：<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear Du,<div><br></div><div>if both radicals are doublet and the product is singlet you should set total spin to 0 while keeping UKS true. This should work automatically.</div><div><br></div><div>Yours,</div><div><br></div><div>Vladimir<br><br>понедельник, 13 апреля 2020 г., 14:25:56 UTC+2 пользователь Du написал:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello All,<div><br></div><div>I have some doubts about the MD simulation, when both the reactant are radical, is there any special parameter should be set in the input profile? or Can cp2k recognize the radical species automatically?</div><div><br></div><div>Thanks!</div><div><br></div><div>Du</div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div>