<div><br></div><div><br></div>Dear Senja, <div><br></div><div>With smearing activated, by default the multiplicity is not conserved.</div><div>Have you already tried the keyword</div><div><ul><li><a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/SCF/SMEAR.html#FIXED_MAGNETIC_MOMENT">FIXED_MAGNETIC_MOMENT</a></li></ul>in the section </div><div><a href="https://manual.cp2k.org/trunk/CP2K_INPUT.html">CP2K_INPUT</a> / <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL.html">FORCE_EVAL</a> / <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT.html">DFT</a> / <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/SCF.html">SCF</a> / <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/SCF/SMEAR.html">SMEAR</a>?<br></div><div><br></div><div>Regards</div><div>Marcella</div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Thursday, February 24, 2022 at 9:26:37 AM UTC+1 senja.r...@gmail.com wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Dear CP2K developers and users,<br><br>I have a question regarding spin polarized calculations and the use of the multiplicity and relax_multilplicity tags. <br><br>I would like to calculate the H-abstraction energy on my Si-terminated 4H-SiC (0001) surface. <br><br>Initial structure: each Si surface atom has one dangling bond which is passivated by hydrogen atom (non magnetic groundstate, magnet moment = 0).<br><br>Final structure: one Si-H bond is broken, the  H atom is moved away from the surface. The dangling bond and H atom now have a spin polarized ground state in which both unbound electrons occupy one spin channel (spin polarized groundstate, magnetic moment = 2).<br><br>I have plotted my results for VASP and CP2K in the attached figure. The nonmagnetic calculations show good overlap (blue and purple curve).<br><br>For vasp, one can fix the total magnetic moment with NUPDOWN = 0/2. If one sets this to 2 and starts with a ferro magnetic initial guess for all atoms (MAGMOM = N*1), the orange curve is retrieved.  One can see that  the groundstate is nonmagnetic until 3.4A and then becomes ferromagnetic. Additionally, one can not set NUPDOWN and let it relax (red and green curve for two different initial guesses). In this case, VASP is close to the correct transition but not perfect. <br><br>In CP2K, once I turn on the magnetization with UKS and MAGNETIZATION = 1 in the &KIND section, the curves are all over the place and hop in and out of different magnetic states. I tried to tune this with the MULTIPLICITY and RELAX_MULTIPLICITY tag without success.  <br><br>My questions are:<br><br>1. How can I fix the total magnetic moment so that I  could reproduce the orange VASP curve with CP2K?<br>2. What's the influence of the multiplicity on this calculation (I tried the default 1 and the triplet/3 settings). I suppose my system is similar to the O2 example where the Si-dangling bond and H electron are similar to the triplet state.<br>3.What does relax_multiplicity do?  It seems to improve some data points (grey vs. brown curve) but not all.<br>4. Which print settings should I use to get an overview of the magnetic moment total (like printed in OSZICAR VASP) and on each atom (comparable to LORBIT = 11 in VASP) <br><br>I would like to start using CP2K more but need to understand the spin polarization. Can I improve my results to come closer to VASP (or maybe even better)?<br><br>Thank you in advance!<br><br>Best wishes,<br>Senja Ramakers (PhD student at Ruhr-Universität Bochum)<br><br> </blockquote></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/d4f93141-cf1f-423e-8538-5d4e3d9168a1n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/d4f93141-cf1f-423e-8538-5d4e3d9168a1n%40googlegroups.com</a>.<br />