<div dir="ltr">Dear Professor Giannozzi,<div><br></div><div>Thank you very much for your help, indeed the spinsq routine does what I need (and testing it within CP gives me good values for <S^2>). If it is not too much to ask, what would be your best quick guess as to where the routine could be interfaced with PWSCF without too much tinkering? I am unfortunately not very familiar with the code structure of QE.</div><div><br></div><div>I see from spinsq that it requires the psi, <psi|beta> and rhor. I can't quite figure out the second variable; in cpr.f90 spinsq is called as </div><div><br></div><div>IF ( tstdout) CALL spinsq ( c0_bgrp, bec_bgrp, rhor )<br></div><div><br></div><div>I believe the first variable is usually called evc in pwscf but I am not sure if bec_bgrp here is the one from becmod in /Modules?</div><div><br></div><div>Any help is much appreciated!</div><div><br></div><div>With best wishes,</div><div>Chris </div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, 4 Aug 2020 at 21:02, Paolo Giannozzi <<a href="mailto:p.giannozzi@gmail.com">p.giannozzi@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>There is a simple routine in CP that estimates the amount of "spin contamination": CPV/src/spinsq.f90. Hope this helps (but I am not sure it will!)<br></div><div><br></div><div>Paolo<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Jul 26, 2020 at 8:36 AM Christoph Wolf <wolf.christoph@qns.science> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Dear all,<div><br></div><div>I am navigating in circles here, I hope someone can "kick" me in the right direction ;)</div><div><br></div><div>When estimating the "spin-contamination" literature (e.g. JP Malrieu, Chem Rev 2014, 114, 429 ff) mentions that the degree of contamination is estimated by the deviation of <BS|S^2|BS> and <HS|S^2|HS> from their ideal (Ising wave-function type) values. I have a system in which the two spin-centers seem spatially well-separated and indeed I find that the LDOS indicates negligible overlap of the magnetic orbitals on the two centers. However I also wanted to check if the Noodleman condition is satisfied, e.g. if <S^2>_{HS}-<S^2>_{BS}~(S_max)^2.</div><div><br></div><div>However, I can't figure out if this is calculated in PWSCF at some point? Or does literature suggest to calculate this "by hand" possibly from the output of the Hubbard U routine where I have account of the individual electrons?</div><div><br></div><div>Any help is appreciated, have a great weekend!</div><div><br></div><div>Best,</div><div>Chris <br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr">Postdoctoral Researcher<br>Center for Quantum Nanoscience, Institute for Basic Science<br>Ewha Womans University, Seoul, South Korea<blockquote type="cite" style="font-size:12.8px"><div dir="ltr"><div><div dir="ltr"></div></div></div></blockquote></div></div></div></div>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div><br clear="all"><br>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr">Postdoctoral Researcher<br>Center for Quantum Nanoscience, Institute for Basic Science<br>Ewha Womans University, Seoul, South Korea<blockquote type="cite" style="font-size:12.8px"><div dir="ltr"><div><div dir="ltr"></div></div></div></blockquote></div></div>