Great, it's correct, thank you so much<br /><br />Guess for atomic kind: Cu<br /><br /> Electronic structure<br />    Total number of core electrons                                         18.00<br />    Total number of valence electrons                                       9.00<br />    Total number of electrons                                              27.00<br />    Multiplicity                                                         doublet<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      4.00<br /> <br />Based on what criteria did you assign NEL 0 -1 to alpha and NEL -2 -1 to beta?<br /><br />Thank you again, you have been very helpful.<div><br /></div><div>Lorenzo</div><div class="gmail_quote"><div dir="auto" class="gmail_attr">Il giorno giovedì 6 marzo 2025 alle 10:54:00 UTC+1 Krack Matthias ha scritto:<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="DE-CH" style="font-size:11.0pt">Dear Lorenzo<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="DE-CH" style="font-size:11.0pt"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt">Sorry, my suggestion for the &BS section was wrong, try this instead<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;font-family:Consolas;color:#212121"> 
</span><span lang="IT-CH" style="font-size:11.0pt;font-family:Consolas;color:#212121">&BS</span><span style="color:#212121"><u></u><u></u></span></p>
<p class="MsoNormal" style="font-variant-caps:normal;text-align:start;word-spacing:0px">
<span lang="IT-CH" style="font-size:11.0pt;font-family:Consolas;color:#212121">    &ALPHA</span><span style="color:#212121"><u></u><u></u></span></p>
<p class="MsoNormal" style="font-variant-caps:normal;text-align:start;word-spacing:0px">
<span lang="IT-CH" style="font-size:11.0pt;font-family:Consolas;color:#212121">      N    3  4</span><span style="color:#212121"><u></u><u></u></span></p>
<p class="MsoNormal" style="font-variant-caps:normal;text-align:start;word-spacing:0px">
<span lang="IT-CH" style="font-size:11.0pt;font-family:Consolas;color:#212121">      L    2  0</span><span style="color:#212121"><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" style="font-variant-caps:normal;text-align:start;word-spacing:0px">
<span lang="IT-CH" style="font-size:11.0pt;font-family:Consolas;color:#212121">      NEL  0 -1</span><span style="color:#212121"><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" style="font-variant-caps:normal;text-align:start;word-spacing:0px">
<span lang="IT-CH" style="font-size:11.0pt;font-family:Consolas;color:#212121">   
</span><span lang="EN-US" style="font-size:11.0pt;font-family:Consolas;color:#212121">&END ALPHA</span><span style="color:#212121"><u></u><u></u></span></p>
<p class="MsoNormal" style="font-variant-caps:normal;text-align:start;word-spacing:0px">
<span lang="EN-US" style="font-size:11.0pt;font-family:Consolas;color:#212121">    &BETA</span><span style="color:#212121"><u></u><u></u></span></p>
<p class="MsoNormal" style="font-variant-caps:normal;text-align:start;word-spacing:0px">
<span lang="EN-US" style="font-size:11.0pt;font-family:Consolas;color:#212121">     N     3  4</span><span style="color:#212121"><u></u><u></u></span></p>
<p class="MsoNormal" style="font-variant-caps:normal;text-align:start;word-spacing:0px">
<span lang="EN-US" style="font-size:11.0pt;font-family:Consolas;color:#212121">     L     2  0</span><span style="color:#212121"><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" style="font-variant-caps:normal;text-align:start;word-spacing:0px">
<span lang="EN-US" style="font-size:11.0pt;font-family:Consolas;color:#212121">     NEL  -2 -1</span><span style="color:#212121"><u></u><u></u></span></p>
<p class="MsoNormal" style="font-variant-caps:normal;text-align:start;word-spacing:0px">
<span lang="EN-US" style="font-size:11.0pt;font-family:Consolas;color:#212121">    &END BETA</span><span style="color:#212121"><u></u><u></u></span></p>
<p class="MsoNormal" style="font-variant-caps:normal;text-align:start;word-spacing:0px">
<span lang="EN-US" style="font-size:11.0pt;font-family:Consolas;color:#212121">  &END BS</span><span style="color:#212121"><u></u><u></u></span></p>
<p class="MsoNormal" style="font-variant-caps:normal;text-align:start;word-spacing:0px">
<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">Best<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>
<div>
<div></div></div></div></div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div><div><div><div>
<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="color:black">From: </span></b><span style="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 lorenzo briccolani <<a href data-email-masked rel="nofollow">lore.br...@gmail.com</a>><br>
<b>Date: </b>Wednesday, 5 March 2025 at 19:54<br>
<b>To: </b><a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>><br>
<b>Subject: </b>Re: [CP2K:21253] Broken Simmetry NEL<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt">Dear Matthias,<br>
<br>
Thank you very much for your help, as always—it is greatly appreciated.<br>
<br>
Just to ensure I fully understand, the first column of values defines the target orbital, while the second column corresponds to the initial orbital. Therefore, in the case of the alpha electron, I add one electron with spin-up (NEL +1) in the 3d orbital by
 removing it from the 4s orbital. For the beta electron, I do not add any electron (NEL -1), transferring it from the 4s to the 3d orbital.<br>
<br>
I apologize if I seem insistent, but I would like to clarify this point, as the manual refers to a single real value.<br>
<br>
Following your recommendations, I assigned an initial multiplicity of two to the Cu atom. However, I am unsure why the calculation still assigns 10 valence electrons overall, while for the beta electrons, it assigns 4.5 instead of 4. Because the values enclosed
 in square brackets, I believe, refer to the number of alpha or beta electrons contained in each orbital. This is why, in the case of beta electrons, I expected D 4.00 instead of 4.50. Below, I have included the relevant lines from the output:<br>
Guess for atomic kind: Cu<br>
<br>
 Electronic structure<br>
    Total number of core electrons                                         18.00<br>
    Total number of valence electrons                                      10.00<br>
    Total number of electrons                                              28.00<br>
    Multiplicity                                                         doublet<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      4.50<br>
<br>
Best Regards <br>
<br>
Lorenzo Briccolani<u></u><u></u></p>
</div>
<p class="MsoNormal" style="margin-left:36.0pt"><u></u> <u></u></p>
<div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt">On Wed, 5 Mar 2025 at 12:02, Krack Matthias <<a href data-email-masked rel="nofollow">matthia...@psi.ch</a>> wrote:<u></u><u></u></p>
</div>
<blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm">
<div>
<div>
<div>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="DE-CH" style="font-size:11.0pt">Dear Lorenzo</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="DE-CH" style="font-size:11.0pt"> </span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt">The valence electron configuration of the Cu-q11 pseudopotential employed by you is 3d10 4s1. The following &BS section should create an initial 3d9 doublet configuration for a Cu(2+) atom:</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt;font-family:Consolas"> </span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt;font-family:Consolas"> UKS</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt;font-family:Consolas">  MULTIPLICITY 2</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt"> </span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt;font-family:Consolas">  &BS</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="IT-CH" style="font-size:11.0pt;font-family:Consolas">    &ALPHA</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="IT-CH" style="font-size:11.0pt;font-family:Consolas">      N    3  4</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="IT-CH" style="font-size:11.0pt;font-family:Consolas">      L    2  0</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt;font-family:Consolas">      NEL +1 -1</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt;font-family:Consolas">    &END ALPHA</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt;font-family:Consolas">    &BETA</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt;font-family:Consolas">     N    3  4</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt;font-family:Consolas">      L    2  0</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt;font-family:Consolas">      NEL -1 -1</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt;font-family:Consolas">    &END BETA</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt;font-family:Consolas">  &END BS</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt">Note that these BS setting affects only the electron configuration of the initial atomic guess and the converged result might differ.</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt"> </span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt">HTH</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt"> </span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt">Matthias</span><u></u><u></u></p>
<p class="MsoNormal" style="margin-left:36.0pt">
<span lang="EN-US" style="font-size:11.0pt"> </span><u></u><u></u></p>
<div>
<div>
<div>
<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="color:black">From: </span></b><span style="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 lorenzo
 briccolani <<a href data-email-masked rel="nofollow">lore.br...@gmail.com</a>><br>
<b>Date: </b>Tuesday, 4 March 2025 at 19:29<br>
<b>To: </b>cp2k <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>><br>
<b>Subject: </b>[CP2K:21245] Broken Simmetry NEL</span><u></u><u></u></p>
</div>
<p class="MsoNormal" style="margin-left:72.0pt">
Dear Users, <br>
<br>
I am reaching out as I am uncertain about the correct assignment of NEL values in the BS routine. The system under investigation consists of a molecule with a central Cu2+ atom coordinating two ligands in a square planar geometry. In a standard spin density
 calculation (UPBE/DZVP-MOLOPT) with a multiplicity of two, the spin densities are well reproduced, with the unpaired electron localized precisely on the copper atom. The Hirshfeld output reports the following values:<br>
<br>
Element  Kind  Ref Charge     Population       Spin moment  Net charge<br>
           Cu     4      11.000        6.182   5.514            0.668     -0.695<br>
<br>
My objective is to replicate this result using the BS approach. For a Cu2+ ion, I assume the configuration should be as follows:<br>
&KIND Cu<br>
     &BS ON<br>
      &ALPHA<br>
       L = 0 2<br>
       N = 4 3<br>
      NEL = -2 ?<br>
      &END<br>
      &BETA<br>
       L = 0 2<br>
       N = 4 3<br>
      NEL = -2 ? <br>
      &END<br>
     &END<br>
<br>
I have conducted several tests; however, I am unable to obtain the correct number of valence electrons and the expected multiplicity in the output. I have attached the input file for reference (opt.inp).<br>
<br>
I would be sincerely grateful for any suggestions regarding the appropriate assignment of the NEL values.<br>
<br>
Best regards,<br>
<br>
Lorenzo Briccolani <u></u><u></u></p>
<p class="MsoNormal" style="margin-left:72.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
<a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>.<br>
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/7b7ad053-f176-49f2-a9ea-0580b9d0d584n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=it&q=https://groups.google.com/d/msgid/cp2k/7b7ad053-f176-49f2-a9ea-0580b9d0d584n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1741342498115000&usg=AOvVaw2qo2CTAhQx3NwrAcqLAh4_">
https://groups.google.com/d/msgid/cp2k/7b7ad053-f176-49f2-a9ea-0580b9d0d584n%40googlegroups.com</a>.<u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
<p class="MsoNormal" style="margin-left:36.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
<a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>.<br>
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/ZRAP278MB08276547DC4F2701A2644E5EF4CB2%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=it&q=https://groups.google.com/d/msgid/cp2k/ZRAP278MB08276547DC4F2701A2644E5EF4CB2%2540ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1741342498115000&usg=AOvVaw2JgwAyaBF_fCJRgkbSTuQd">
https://groups.google.com/d/msgid/cp2k/ZRAP278MB08276547DC4F2701A2644E5EF4CB2%40ZRAP278MB0827.CHEP278.PROD.OUTLOOK.COM</a>.<u></u><u></u></p>
</div>
</blockquote>
</div>
</div></div></div></div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div><div><div><div><p class="MsoNormal" style="margin-left:36.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
<a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>.<br></p></div></div></div></div></div><div lang="en-CH" link="blue" vlink="purple" style="word-wrap:break-word"><div><div><div><div><p class="MsoNormal" style="margin-left:36.0pt">
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/CAJPdbE-4Apgw%3DoyYG-LT%2BmkkwVrMJmfEw29VQ9VVqohHoCF85A%40mail.gmail.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=it&q=https://groups.google.com/d/msgid/cp2k/CAJPdbE-4Apgw%253DoyYG-LT%252BmkkwVrMJmfEw29VQ9VVqohHoCF85A%2540mail.gmail.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1741342498115000&usg=AOvVaw0xT2MH1zY_KYAKl5YZSAYq">
https://groups.google.com/d/msgid/cp2k/CAJPdbE-4Apgw%3DoyYG-LT%2BmkkwVrMJmfEw29VQ9VVqohHoCF85A%40mail.gmail.com</a>.<u></u><u></u></p>
</div>
</div>
</div>
</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 visit <a href="https://groups.google.com/d/msgid/cp2k/ded38ae6-21e6-4508-8043-572f01e461ddn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/ded38ae6-21e6-4508-8043-572f01e461ddn%40googlegroups.com</a>.<br />