Dear Krack<div><br /></div><div>Thank you so much for your comment. I performed a calculation with this set up and the output now reads:</div><div><br /></div><div>Electronic structure<br />    Total number of core electrons                                         10.00<br />    Total number of valence electrons                                      13.00<br />    Total number of electrons                                              23.00<br />    Multiplicity                                                          sextet<br />    Alpha Electrons<br />    S   [  1.00  1.00] 1.00<br />    P   [  3.00] 3.00<br />    D      5.00<br />    Beta Electrons<br />    S   [  1.00  1.00] 1.00<br />    P   [  3.00] 3.00<br />    D      0.00</div><div><br /></div><div>The occupations of the p and d orbitals now look perfect, but if I am not wrong the occupation of the s orbitals should be 0. Is this ok for you or there is some problems here?</div><div><br /></div><div>Best,</div><div>Diego<br /><br /></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">El martes, 30 de mayo de 2023 a las 15:13:48 UTC+2, Krack Matthias escribió:<br/></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-CH" link="blue" vlink="purple" style="word-wrap:break-word">
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">Hi Diego<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">For Fe(+3) d^5 spin-up try<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Consolas">&ALPHA<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Consolas">N    4  3<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Consolas">L    0  2<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Consolas">NEL -2 +4<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Consolas">&END ALPHA
<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Consolas">&BETA<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Consolas">N    4  3<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Consolas">L    0  2<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Consolas">NEL -2 -6<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">&END BETA<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">Switch ALPHA and BETA for spin down.<u></u><u></u></span></p></div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">HTH<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">Matthias<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt"><u></u> <u></u></span></p>
<div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal" style="margin-right:0cm;margin-bottom:12.0pt;margin-left:36.0pt">
<b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt;color:black"><a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of Diego López <<a href data-email-masked rel="nofollow">diego...@gmail.com</a>><br>
<b>Date: </b>Tuesday, 30 May 2023 at 14:22<br>
<b>To: </b>cp2k <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>><br>
<b>Subject: </b>Re: [CP2K:18867] How to correctly set magnetization in CP2K<u></u><u></u></span></p>
</div>
<p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">Dear Krack,<br>
<br>
thank you a lot for your comments. The &BS section seems to be really usefull but I am struggling to obtain the actual electronic configuration of my system (Y3Fe5O12). I am trying to model theese Fe3+ (d5) ions, where Fe1 will be spin up and Fe2 will be spin
 down using the attatched &BS section on the input. <u></u><u></u></span></p>
<div>
<p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">&KIND Fe1<br>
      ELEMENT Fe<br>
      BASIS_SET SZV-MOLOPT-SR-GTH<br>
      POTENTIAL GTH-PBE-q16<br>
      &BS<br>
       &ALPHA<br>
       NEL -2 -1<br>
       L    0  2<br>
       N    4  3<br>
      &END<br>
      &BETA<br>
       NEL -2 -4<br>
       L    0  2<br>
       N    4  3<br>
      &END<br>
      &END<br>
    &END KIND<br>
    &KIND Fe2<br>
      ELEMENT Fe<br>
      BASIS_SET SZV-MOLOPT-SR-GTH<br>
      POTENTIAL GTH-PBE-q16<br>
      &BS<br>
       &ALPHA<br>
       NEL -2 -4<br>
       L    0  2<br>
       N    4  3<br>
      &END<br>
      &BETA<br>
       NEL -2 -1<br>
       L    0  2<br>
       N    4  3<br>
      &END<br>
      &END<br>
    &END KIND<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">The electronic configuration of Fe is [Ar] 4s2 3d6 and  the one of Fe3+ ions is [Ar] 3d5, so I tried to extract both 4s electrons and 1 of 3d in each spin orientation, but in the
 output I see that there is something wrong:<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">Guess for atomic kind: Fe1<br>
<br>
 Electronic structure<br>
    Total number of core electrons                                         10.00<br>
    Total number of valence electrons                                      11.50<br>
    Total number of electrons                                              21.50<br>
    Multiplicity                                                         triplet<br>
    Alpha Electrons<br>
    S   [  1.00  1.00] 1.00<br>
    P   [  3.00] 3.00<br>
    D      2.50<br>
    Beta Electrons<br>
    S   [  1.00  1.00] 1.00<br>
    P   [  3.00] 3.00<br>
    D      1.00<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">Guess for atomic kind: Fe2<br>
<br>
 Electronic structure<br>
    Total number of core electrons                                         10.00<br>
    Total number of valence electrons                                      11.50<br>
    Total number of electrons                                              21.50<br>
    Multiplicity                                                         triplet<br>
    Alpha Electrons<br>
    S   [  1.00  1.00] 1.00<br>
    P   [  3.00] 3.00<br>
    D      1.00<br>
    Beta Electrons<br>
    S   [  1.00  1.00] 1.00<br>
    P   [  3.00] 3.00<br>
    D      2.50<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">How could I configure the &BS section to obtain an accurate description of my system?<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">Thanks in advanced.<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">Best,<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-right:0cm;margin-bottom:12.0pt;margin-left:36.0pt">
<span style="font-size:11.0pt">Diego<u></u><u></u></span></p>
</div>
<div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">El martes, 30 de mayo de 2023 a las 13:31:11 UTC+2, Krack Matthias escribió:<u></u><u></u></span></p>
</div>
<blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt">
<div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt">Dear Diego</span><span style="font-size:11.0pt"><u></u><u></u></span></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><u></u><u></u></span></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt">Both keywords can be used to define the initial atomic states and the total magnetization of the system in conjunction with the MULTIPLICITY keyword. The &BS section provides a more detailed (orbital resolved) control.</span><span style="font-size:11.0pt"><u></u><u></u></span></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><u></u><u></u></span></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt">HTH</span><span style="font-size:11.0pt"><u></u><u></u></span></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><u></u><u></u></span></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt">Matthias</span><span style="font-size:11.0pt"><u></u><u></u></span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt"> </span><span style="font-size:11.0pt"><u></u><u></u></span></p>
<div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal" style="margin-bottom:12.0pt;margin-left:72.0pt">
<b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt;color:black"><span><a href data-email-masked rel="nofollow">cp...@googlegroups.com</a></span> <<span><a href data-email-masked rel="nofollow">cp...@googlegroups.com</a></span>> on behalf of Diego López <<span><a href data-email-masked rel="nofollow">diego...@gmail.com</a></span>><br>
<b>Date: </b>Monday, 29 May 2023 at 17:25<br>
<b>To: </b>cp2k <<span><a href data-email-masked rel="nofollow">cp...@googlegroups.com</a></span>><br>
<b>Subject: </b>[CP2K:18863] How to correctly set magnetization in CP2K</span><span style="font-size:11.0pt"><u></u><u></u></span></p>
</div>
<p class="MsoNormal" style="margin-left:72.0pt">
<span style="font-size:11.0pt">Dear CP2k users, I am interested on the study of Y3Fe5O12 and now I am trying to set the adequate magnetization for Fe atoms. The unit cell has 40 Fe atoms and the magnetic ground state experimentally observed is ferrimagnetic
 (24 Fe atoms spin up and 16 Fe atoms spin down).<u></u><u></u></span></p>
<div>
<p class="MsoNormal" style="margin-left:72.0pt">
<span style="font-size:11.0pt"> <u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:72.0pt">
<span style="font-size:11.0pt">Which is the proper way to specify the magnetization of the system? &MAGNETIZATION or &BS? Are they compatible? How would be the &BS specification for this case?<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:72.0pt">
<span style="font-size:11.0pt"> <u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:72.0pt">
<span style="font-size:11.0pt">Thanks in advanced. <u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:72.0pt">
<span style="font-size:11.0pt"> <u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:72.0pt">
<span style="font-size:11.0pt">Regards,<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:72.0pt">
<span style="font-size:11.0pt">Diego<u></u><u></u></span></p>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal" style="margin-left:72.0pt">
<span style="font-size:11.0pt">-- <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
<span><a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a></span>.<br>
To view this discussion on the web visit </span><a href="https://groups.google.com/d/msgid/cp2k/c7ef7446-40ba-46f0-a728-f4a9ef2309ddn%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=es&q=https://groups.google.com/d/msgid/cp2k/c7ef7446-40ba-46f0-a728-f4a9ef2309ddn%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1685540048347000&usg=AOvVaw3Ytz6VnTc-236hNj0n-iEK"><span style="font-size:11.0pt">https://groups.google.com/d/msgid/cp2k/c7ef7446-40ba-46f0-a728-f4a9ef2309ddn%40googlegroups.com</span></a><span style="font-size:11.0pt">.<u></u><u></u></span></p>
</div>
</div>
</blockquote>
</div>
</div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div><p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">--
<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
</span><a href data-email-masked rel="nofollow"><span style="font-size:11.0pt">cp2k+uns...@googlegroups.com</span></a></p></div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div><p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:11.0pt">.<br>
To view this discussion on the web visit </span><a href="https://groups.google.com/d/msgid/cp2k/911e3f9e-b68b-45bd-8f56-b1c88dba101bn%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=es&q=https://groups.google.com/d/msgid/cp2k/911e3f9e-b68b-45bd-8f56-b1c88dba101bn%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1685540048347000&usg=AOvVaw2wQSrsStnQYxSFExCeCv15"><span style="font-size:11.0pt">https://groups.google.com/d/msgid/cp2k/911e3f9e-b68b-45bd-8f56-b1c88dba101bn%40googlegroups.com</span></a><span style="font-size:11.0pt">.<u></u><u></u></span></p>
</div>
</div>

</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/b91ad0e7-83dc-459c-be77-62a081fe8b85n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/b91ad0e7-83dc-459c-be77-62a081fe8b85n%40googlegroups.com</a>.<br />