<div dir="ltr"><a name="_GoBack"></a>Dear CP2K users,<o:p></o:p><o:p> </o:p>I would like some clarifications about SMEAR vs. LSD +RELAX_MULTIPLICTY in cp2k. I am trying to simulate a metallic nanoparticle (few tens of atoms) in water<o:p></o:p><o:p> </o:p>I have two (maybe simple) questions:<o:p></o:p><o:p> </o:p>a) From this post <o:p></o:p><a href="https://groups.google.com/forum/#!topic/cp2k/GEKLwRQat_M" style="color: rgb(149, 79, 114); text-decoration: underline;">https://groups.google.com/forum/#!topic/cp2k/GEKLwRQat_M</a><o:p></o:p><o:p> </o:p>my understanding is that the SMEAR option populates the alpha and beta channel according to a Fermi-Dirac <o:p></o:p><o:p> </o:p>In LSD + MULTIPLICITY (or LSD + RELAX MULTIPLICITY)<o:p></o:p><a href="https://groups.google.com/forum/#!searchin/cp2k/smearing$20spin%7Csort:date/cp2k/ugkJLz4PSSI/C-ProBR6BwAJ" style="color: rgb(149, 79, 114); text-decoration: underline;">https://groups.google.com/forum/#!searchin/cp2k/smearing$20spin%7Csort:date/cp2k/ugkJLz4PSSI/C-ProBR6BwAJ</a><o:p></o:p><o:p> </o:p>my understating is that after energy levels are calculated, electrons are placed according to Aufbau principle, from the lowest to the highest in an unrestricted scheme. If MULTIPLICITY is set, the multiplicity is fixed, with RELAX_MULTIPLICITY no<o:p></o:p><o:p> </o:p>Please correct me if my understanding is mistaken.<o:p></o:p><o:p> </o:p>Question: are the two approaches (SMEAR vs. LSD +RELAX_MULTIPLICTY) equivalent? One approach is more appropriate than the other for particular systems? Is the computational cost the same?<o:p></o:p><o:p> </o:p>b) If I use LSD +MULTIPLICTY in the output I get the<o:p></o:p><o:p> </o:p>“Integrated absolute spin density” and “Ideal and single determinant S**2”<o:p></o:p><o:p> </o:p>That are equations (2.271) and (2.272) in “Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory,Attila Bangha Szabo, Neil S. Ostlund “<o:p></o:p><o:p> </o:p><a href="https://groups.google.com/forum/#!topic/cp2k/flEEaseAd7c" style="color: rgb(149, 79, 114); text-decoration: underline;">https://groups.google.com/forum/#!topic/cp2k/flEEaseAd7c</a><o:p></o:p><o:p> </o:p><o:p> </o:p>If I see the output of my calculation with SMEAR I do not see such or any other information about the spin state of my system. <o:p></o:p><o:p> </o:p>Question: Is this expected? How I can easily get information about the spin state of my system when using SMEAR?<o:p></o:p><o:p> </o:p>Thank you very much for any possible help<o:p></o:p><o:p> </o:p>Best regards<o:p></o:p>Ivan<o:p></o:p></div>