<div dir="ltr">To provide more information I'm sending my output, sorry for not including it in my previous message.<br><div><br><div class="gmail_quote"><div dir="ltr">Em seg, 3 de ago de 2015 às 14:17, Henrique Junior <<a href="mailto:henri...@gmail.com">henri...@gmail.com</a>> escreveu:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">


        
        
        
        


        
        
        
        


<p style="margin-bottom:0cm;line-height:100%" lang="en-US">I'm
sorry if I'm bringing up a topic that was already discussed in
previous threads: Broken Symmetry, but things are still a little
confusing to me on how too do it in CP2K.</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">NWChem
approaches this problem by generating a state that has “an excess
of beta electrons”[1] and you move your unpaired electrons to beta
(spin down). At first I was thinking that CP2K was using a similar
approach with &ALPHA and &BETA in &BS. Since I'm working
with a system with no charge, S = 3, and two Cu(II) ( L-CuCl2-L) I
have to investigate the magnetic coupling between copper centers at
various possible distances in my unitary cell. Taking a look at the
examples I tried this:</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US"><br>

</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">   
&KIND  Cu1</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">     
ELEMENT Cu</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">     
BASIS_SET  DZVP-MOLOPT-SR-GTH</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">     
POTENTIAL  GTH-PBE-q11</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">     
&BS</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">       
&ALPHA</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">       
  NEL -1  -1    <b>(3d9 – One unpaired electron up)</b></p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">       
  L    0   2</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">       
  N    4   3</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">       
&END</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">       
&BETA</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">       
  NEL -1  -2   <b>(3d8 – No unpaired electrons)</b></p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">       
  L    0   2</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">       
  N    4   3</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">       
&END</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">     
&END</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">   
&END KIND</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">   
&KIND  Cu2</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">     
ELEMENT Cu</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">     
BASIS_SET  DZVP-MOLOPT-SR-GTH</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">     
POTENTIAL  GTH-PBE-q11</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">     
&BS</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">       
&ALPHA</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">       
  NEL -1  -2  <b>(3d8 – No unpaired electrons)</b></p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">       
  L    0   2</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">       
  N    4   3</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">       
&END</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">       
&BETA</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">       
  NEL -1  -1   <b>(3d9 – One unpaired electron down)</b></p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">       
  L    0   2</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">       
  N    4   3</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">       
&END</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">     
&END</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">   
&END KIND</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US"><br>

</p>
<p style="margin-bottom:0cm;line-height:100%" lang="en-US">Does
that make any sense?</p>

<p style="margin-bottom:0cm;line-height:100%"><br></p><p style="margin-bottom:0cm;line-height:100%">Thank you so much</p><p style="margin-bottom:0cm;line-height:100%"><br></p><p style="margin-bottom:0cm;line-height:100%">[1] --
<a href="http://www.nwchem-sw.org/index.php/Special:AWCforum/st/id1680/" target="_blank">http://www.nwchem-sw.org/index.php/Special:AWCforum/st/id1680/</a></p><p style="margin-bottom:0cm;line-height:100%"><br></p><p style="margin-bottom:0cm;line-height:100%">---------------</p><p style="margin-bottom:0cm;line-height:100%">Henrique C. S. Junior - UFRRJ<br></p><p style="margin-bottom:0cm;line-height:100%">CP2K - 2.6<br></p>

</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+uns...@googlegroups.com" target="_blank">cp2k+uns...@googlegroups.com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="http://groups.google.com/group/cp2k" target="_blank">http://groups.google.com/group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" target="_blank">https://groups.google.com/d/optout</a>.<br>
</blockquote></div></div></div><div dir="ltr">-- <br></div><p dir="ltr">Henrique C. S. Junior </p>