<div dir="ltr">Dear Matthias and Marcella<div> </div><div>thank you very much for your reply. </div><div>Now is much more clear.</div><div> </div><div>However, I still have one doubt.</div><div> </div><div>While in LSD we know the spin state of the system (i.e., S**2 and integrated spin mentioned above) at any step,</div><div>it seems to me that in SMEAR we lost such information: as Matthias said, .S**2 and integrated spin cannot be calculated for fractional occupation in SMEAR.</div><div> </div><div>Using SMEAR, is it possible to define some average quantities defining the spin state of the system? </div><div>I do not know, for example, the average occupation of the alpha and beta channels?</div><div>Or something else?</div><div> </div><div>Thank you very very much for your replies and help</div><div> </div><div>Best</div><div>Ivan </div><div> </div><div> On Tuesday, August 18, 2020 at 5:54:48 PM UTC+3, Marcella Iannuzzi wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">Dear Ivan, <div> </div><div>in addition to what Matthias already wrote, by default the multiplicity is not fixed with SMEAR + LSD. </div><div> </div><div>Working with fractional orbital occupations requires the use of diagonalisation schemes, which are in general more demanding in terms of computational costs if compared with the orbital transformation (OT) method. </div><div> </div><div>Kind regards</div><div>Marcella</div><div> </div><div class="gmail_quote"><div dir="auto">On Tuesday, August 18, 2020 at 4:42:07 PM UTC+2 Matthias Krack wrote: </div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"> <div lang="EN-US" link="blue" vlink="purple"> <div> Dear Ivan a) SMEAR allows for fractional orbital occupations whereas LSD with RELAX_MULTIPLICITY is meant for systems with integer occupations, e.g. to allow for the spin flip needed to switch automatically from a singlet to a possibly lower lying triplet state. Thus, these schemes are not equivalent. b) The expectation value of S**2 is only defined for calculations with integer orbital occupations Therefore this properties is not calculated/printed when SMEAR is active. HTH Matthias</div></div><div lang="EN-US" link="blue" vlink="purple"><div> From: <a rel="nofollow">c...@googlegroups.com</a> <<a rel="nofollow">c...@googlegroups.com</a>> On Behalf Of Ivan Gladich Sent: Dienstag, 18. August 2020 15:56 To: cp2k <<a rel="nofollow">c...@googlegroups.com</a>> Subject: [CP2K:13797] SMEAR vs. LSD+RELAX_MULTIPLICITY/MULTIPLICITY <div> <a name="02c52b4a-4c82-4970-8fa5-63fad726406fn@googlegroups.com_m_-8592122066790477932__GoBack" rel="nofollow"></a>Dear CP2K users, 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 I have two (maybe simple) questions: a) From this post <a href="https://groups.google.com/forum/#!topic/cp2k/GEKLwRQat_M" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/forum/#!topic/cp2k/GEKLwRQat_M';return true;" onclick="this.href='https://groups.google.com/forum/#!topic/cp2k/GEKLwRQat_M';return true;">https://groups.google.com/forum/#!topic/cp2k/GEKLwRQat_M</a> my understanding is that the SMEAR option populates the alpha and beta channel according to a Fermi-Dirac In LSD + MULTIPLICITY (or LSD + RELAX MULTIPLICITY) <a href="https://groups.google.com/forum/#!searchin/cp2k/smearing$20spin%7Csort:date/cp2k/ugkJLz4PSSI/C-ProBR6BwAJ" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/forum/#!searchin/cp2k/smearing$20spin%7Csort:date/cp2k/ugkJLz4PSSI/C-ProBR6BwAJ';return true;" onclick="this.href='https://groups.google.com/forum/#!searchin/cp2k/smearing$20spin%7Csort:date/cp2k/ugkJLz4PSSI/C-ProBR6BwAJ';return true;">https://groups.google.com/forum/#!searchin/cp2k/smearing$20spin%7Csort:date/cp2k/ugkJLz4PSSI/C-ProBR6BwAJ</a> 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 Please correct me if my understanding is mistaken. 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? b) If I use LSD +MULTIPLICTY in the output I get the “Integrated absolute spin density” and “Ideal and single determinant S**2” That are equations (2.271) and (2.272) in “Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory,Attila Bangha Szabo, Neil S. Ostlund “ <a href="https://groups.google.com/forum/#!topic/cp2k/flEEaseAd7c" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/forum/#!topic/cp2k/flEEaseAd7c';return true;" onclick="this.href='https://groups.google.com/forum/#!topic/cp2k/flEEaseAd7c';return true;">https://groups.google.com/forum/#!topic/cp2k/flEEaseAd7c</a> 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. Question: Is this expected? How I can easily get information about the spin state of my system when using SMEAR? Thank you very much for any possible help Best regards Ivan </div> </div></div><div lang="EN-US" link="blue" vlink="purple"><div>-- You received this message because you are subscribed to the Google Groups "cp2k" group. To unsubscribe from this group and stop receiving emails from it, send an email to <a rel="nofollow">cp...@googlegroups.com</a>. To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/80e157d7-5559-4066-bbf5-bca66ce52addo%40googlegroups.com?utm_medium=email&utm_source=footer" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/d/msgid/cp2k/80e157d7-5559-4066-bbf5-bca66ce52addo%40googlegroups.com?utm_medium\x3demail\x26utm_source\x3dfooter';return true;" onclick="this.href='https://groups.google.com/d/msgid/cp2k/80e157d7-5559-4066-bbf5-bca66ce52addo%40googlegroups.com?utm_medium\x3demail\x26utm_source\x3dfooter';return true;"> https://groups.google.com/d/msgid/cp2k/80e157d7-5559-4066-bbf5-bca66ce52addo%40googlegroups.com</a>. </div> </div> </blockquote></div></blockquote></div></div>