<div dir="ltr">Dear all,<br>I am trying to set up calculation for my Cu+ system, however, even after reading this thread I cannot understand how the electrons are assigned to channels in CP2K. I wanna test two possible spin states (just for testing purposes) of Cu+:<br><br>Cu+ with electron configuration 3d9 4s1 (Cu1) and 3d10 4s0 (Cu2). I tried to build my input sections:<br><br>&KIND  Cu1<br>      ELEMENT Cu<br>      BASIS_SET  DZVP-MOLOPT-GTH<br>      POTENTIAL  GTH-BLYP-q11<br>      &BS<br>        &ALPHA<br>           &ALPHA<br>          NEL 0   0<br>          L    0   2<br>          N    4   3<br>        &END<br>        &BETA<br>          NEL -1  -2<br>          L    0   2<br>          N    4   3<br>      &END<br>    &END<br><br>and <br><br>&KIND  Cu2<br>      ELEMENT Cu<br>      BASIS_SET  DZVP-MOLOPT-GTH<br>      POTENTIAL  GTH-BLYP-q11<br>      &BS<br>        &ALPHA<br>           &ALPHA<br>          NEL -1   0<br>          L    0   2<br>          N    4   3<br>        &END<br>        &BETA<br>          NEL -1  0<br>          L    0   2<br>          N    4   3<br>      &END<br>    &END<br><br>Is this correct?<br><div><br></div><div>Also,regarding previous example on Cl1 - what is the meaning of half-electrons in alpha/beta channel? <br></div><div><br></div>Thank you very much for the help<br>Katarina<br><br><br>Dňa streda, 10. februára 2016 11:53:26 UTC+1 S Ling napísal(-a):<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">You cannot have half an electron on a p orbital?<div><br></div><div>In your case of Cl1, you have one extra electron in the alpha spin channel, and by doing that, you break the symmetry of alpha and beta spin channels, i.e. you MAY end up with different occupations in the alpha and beta spin channels after the calculation. Of course, this also depends on your specific system, that a symmetry-broken solution is energetically more stable.</div><div><br></div><div>SL</div><div><br></div></div><div><br><div class="gmail_quote">On 9 February 2016 at 07:49, tao liu <span dir="ltr"><<a href="javascript:" target="_blank" gdf-obfuscated-mailto="mK96HMPnGAAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">liu...@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear All,<br><br>I can understand the &BS setting for Cu1, but still confused on Cl and Cl1.<br><br>    &KIND  Cl<br>      BASIS_SET  DZVP-MOLOPT-GTH<br>      POTENTIAL  GTH-BLYP-q7<br>      &BS<br>        &ALPHA<br>          NEL 2<br>          L   1<br>          N   3<br>        &END<br>        &BETA<br>          NEL 2<br>          L   1<br>          N   3<br>        &END<br>      &END<br>    &END<br>    &KIND  Cl1<br>      ELEMENT Cl<br>      BASIS_SET  DZVP-MOLOPT-GTH<br>      POTENTIAL  GTH-BLYP-q7<br>      &BS<br>        &ALPHA<br>          NEL 2<br>          L   1<br>          N   3<br>        &END<br>        &BETA<br>          NEL 0<br>          L   1<br>          N   3<br>        &END<br>      &END<br>    &END<br><br>As suggested by Marcella, <br>for 'Cl' in the input file, Alpha is (5+2)/2=3.5, and same for Beta, in the end 'Cl' should have 3s-2 3p-7<br>for 'Cl1' in the input file, Alpha is (5+2)/2=3.5, and (5+0)/2 =2.5 for Beta, in the end 'Cl1' should have 3s-2 3p-6<br> <br>Could anyone tell me if I understand correctly or not ?<br><br>Thanks,<br>Tao<span><font color="#888888"><br><br></font></span></div><span><font color="#888888">

<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="javascript:" target="_blank" gdf-obfuscated-mailto="mK96HMPnGAAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp2k+uns...@googlegroups.<wbr>com</a>.<br>
To post to this group, send email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="mK96HMPnGAAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;">https://groups.google.com/<wbr>group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>.<br>
</font></span></blockquote></div><br></div>
</blockquote></div>