<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><div dir='ltr'><font style="" face="Courier New,sans-serif" color="#000000">Hi, S Ling, thank you for your answer, things are starting to make much more sense now.<br>After trying the values I still have the NaN issue:<br><br>Spin 1<br><br> Guess for atomic kind: Cu1<br><br> Electronic structure<br>    Total number of core electrons                                         18.00<br>    Total number of valence electrons                                      15.00<br>    Total number of electrons                                              33.00<br>    Multiplicity                                                   not specified<br>    S   [  2.00  2.00  2.00]<br>    P   [  6.00  6.00]<br>    D     15.00<br><br><br> *******************************************************************************<br>                  Iteration          Convergence                     Energy [au]<br> *******************************************************************************<br>                          1         0.00000                                  NaN<br><br> Energy components [Hartree]           Total Energy ::                       NaN<br>                                        Band Energy ::                       NaN<br>                                     Kinetic Energy ::                       NaN<br>                                   Potential Energy ::                       NaN<br>                                      Virial (-V/T) ::                       NaN<br>                                        Core Energy ::                       NaN<br>                                     Coulomb Energy ::                       NaN<br>                       Total Pseudopotential Energy ::                       NaN<br>                       Local Pseudopotential Energy ::                       NaN<br>                    Nonlocal Pseudopotential Energy ::                       NaN<br>                                        Confinement ::                       NaN<br><br> Orbital energies  State     L     Occupation   Energy[a.u.]          Energy[eV]<br><br>                       1     2         15.000            NaN                 NaN<br><br>---------------------------------------------------------------------------------------<br>  ----------------------------------- OT ---------------------------------------<br><br>  Step     Update method      Time    Convergence         Total energy    Change<br>  ------------------------------------------------------------------------------<br><br>  Trace(PS):                                             NaN<br>  Electronic density on regular grids:                   NaN                 NaN<br>  Core density on regular grids:              449.9998255599       -0.0001744401<br>  Total charge density on r-space grids:                 NaN<br>  Total charge density g-space grids:                    NaN<br><br>I was thinking that CP2k was starting from d10 in &ALPHA instead of d5 and from d10 in &BETA too. In this case, since I'm still getting the NaN, is it possible that something like this:<br>    &KIND  Cu1
<br>      ELEMENT Cu
<br>      BASIS_SET  DZVP-MOLOPT-SR-GTH
<br>      POTENTIAL  GTH-PBE-q11
<br>      &BS
<br>        &ALPHA
<br>          NEL -1   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>    &END KIND
<br>    &KIND  Cu2
<br>      ELEMENT Cu
<br>      BASIS_SET  DZVP-MOLOPT-SR-GTH
<br>      POTENTIAL  GTH-PBE-q11
<br>      &BS
<br>        &ALPHA
<br>          NEL -1  -2
<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>    &END KIND<br><br>Can give me the desired 5 electrons spin up and 4 spin down (and vice versa)?<br></font><br><div><font style="" face="Courier New">----------</font></div><font style="" face="Courier New" color="#666666"><b>Henrique C. S. Junior</b><br></font><div><font style="" face="Courier New" color="#666666">Químico Industrial - UFRRJ<br>Mestrando em Química Inorgânica - UFRRJ<br></font></div><div><font style="" face="Courier New" color="#666666">Centro de Processamento de Dados - PMP</font></div><br><br><div><hr id="stopSpelling">Date: Wed, 5 Aug 2015 00:53:39 +0100<br>Subject: Re: [CP2K:6868] NaN Results with &BS<br>From: lingsa...@gmail.com<br>To: cp...@googlegroups.com<br><br><div dir="ltr"><div>Hi</div><div><br></div><div>For the d9 electron configuration in a spin-polarised calculation, you should have 5 d electrons in the spin ALPHA channel, and 4 d electrons in the spin BETA channel (or 4 spin ALPHA and 5 spin BETA). Therefore, for Cu1, you should use</div><div><br></div><div>&BS<br> &ALPHA<br>   NEL -1 5  #((5+5)/2=5)<br>   L 0 2<br>   N 4 3<br> &END<br> &BETA<br>   NEL -1 3  #((5+3)/2=4)<br>   L 0 2<br>   N 4 3<br> &END<br>&END BS</div><div><br></div><div>and for Cu2, you should use</div><div><br></div><div>&BS<br> &ALPHA<br>   NEL -1 3  #((5+3)/2=4)<br>   L 0 2<br>   N 4 3<br> &END<br> &BETA<br>   NEL -1 5  #((5+5)/2=5)<br>   L 0 2<br>   N 4 3<br> &END<br>&END BS</div><div><br></div><div>SL</div></div><div class="ecxgmail_extra"><br><div class="ecxgmail_quote">On 4 August 2015 at 21:29, Henrique Junior <span dir="ltr"><<a href="mailto:henri...@gmail.com" target="_blank">henri...@gmail.com</a>></span> wrote:<br><blockquote class="ecxgmail_quote" style="border-left:1px #ccc solid;padding-left:1ex;"><div dir="ltr">


        
        
        
        <p style="line-height:100%;" lang="en-US">Dear
CP2K group,</p>
<p style="line-height:100%;" lang="en-US">Trying
to perform a broken symmetry in a structure with two Cu(II) centers,
as I have understood from here [1] (S Ling) and here [2] (Marcella)
there is a scaling factor = 2 to take into consideration when setting
my &ALPHA and &BETA channels. Copper is [Ar] 3d10 4s1 and
since I want to get Cu(II) my goal is [Ar] 3d9 (with one unpaired
electron).</p>
<p style="line-height:100%;" lang="en-US">Taking
into account the scaling factor: 
</p>
<p style="line-height:100%;" lang="en-US"><br>

</p>
<p style="line-height:100%;" lang="en-US">   
&KIND  Cu1</p>
<p style="line-height:100%;" lang="en-US">     
ELEMENT Cu</p>
<p style="line-height:100%;" lang="en-US">     
BASIS_SET  DZVP-MOLOPT-SR-GTH</p>
<p style="line-height:100%;" lang="en-US">     
POTENTIAL  GTH-PBE-q11</p>
<p style="line-height:100%;" lang="en-US">     
&BS</p>
<p style="line-height:100%;" lang="en-US">       
&ALPHA</p>
<p style="line-height:100%;" lang="en-US">       
  NEL -1   8  <b>(because (10 + 8)/2 = 9 and I want 3d9, unpaired e-
spin up) </b>
</p>
<p style="line-height:100%;" lang="en-US">       
  L    0   2</p>
<p style="line-height:100%;" lang="en-US">       
  N    4   3</p>
<p style="line-height:100%;" lang="en-US">       
&END</p>
<p style="line-height:100%;" lang="en-US">       
&BETA</p>
<p style="line-height:100%;" lang="en-US">       
  NEL -1  -10   <b>(because (10 - 10)/2 = 0 and I want 3d0, no e-) </b>
</p>
<p style="line-height:100%;" lang="en-US">       
  L    0   2</p>
<p style="line-height:100%;" lang="en-US">       
  N    4   3</p>
<p style="line-height:100%;" lang="en-US">       
&END</p>
<p style="line-height:100%;" lang="en-US">     
&END</p>
<p style="line-height:100%;" lang="en-US">   
&END KIND</p>
<p style="line-height:100%;" lang="en-US">   
&KIND  Cu2</p>
<p style="line-height:100%;" lang="en-US">     
ELEMENT Cu</p>
<p style="line-height:100%;" lang="en-US">     
BASIS_SET  DZVP-MOLOPT-SR-GTH</p>
<p style="line-height:100%;" lang="en-US">     
POTENTIAL  GTH-PBE-q11</p>
<p style="line-height:100%;" lang="en-US">     
&BS</p>
<p style="line-height:100%;" lang="en-US">       
&ALPHA</p>
<p style="line-height:100%;" lang="en-US">       
  NEL -1  -10 <b>(because (10 - 10)/2 = 0 and I want 3d0, no e-) </b>
</p>
<p style="line-height:100%;" lang="en-US">       
  L    0   2</p>
<p style="line-height:100%;" lang="en-US">       
  N    4   3</p>
<p style="line-height:100%;" lang="en-US">       
&END</p>
<p style="line-height:100%;" lang="en-US">       
&BETA</p>
<p style="line-height:100%;" lang="en-US">       
  NEL -1   8   <b>(because (10 + 8)/2 = 9 and I want 3d9, unpaired e-
spin down)</b></p>
<p style="line-height:100%;" lang="en-US">       
  L    0   2</p>
<p style="line-height:100%;" lang="en-US">       
  N    4   3</p>
<p style="line-height:100%;" lang="en-US">       
&END</p>
<p style="line-height:100%;" lang="en-US">     
&END</p>
<p style="line-height:100%;" lang="en-US">   
&END KIND</p>
<p style="line-height:100%;" lang="en-US"><br>

</p>
<p style="line-height:100%;" lang="en-US">The
results (attached) is a myriad of NaNs and 0.000000000</p>
<p style="line-height:100%;" lang="en-US">Any
ideas will be very much appreciated.</p>
<p style="line-height:100%;" lang="en-US"><br>

</p>
<p style="line-height:100%;" lang="en-US">[1] -
<a href="https://groups.google.com/forum/#%21searchin/cp2k/Broken$20symmetry/cp2k/ghhbtIxUzKk/uhwgb5PMTsgJ" target="_blank">https://groups.google.com/forum/#!searchin/cp2k/Broken$20symmetry/cp2k/ghhbtIxUzKk/uhwgb5PMTsgJ</a></p>
<p style="line-height:100%;" lang="en-US">[2] -
<a href="https://groups.google.com/forum/#%21searchin/cp2k/$26BS$20Marci/cp2k/8fTVlCEjSME/iUem9cK1-McJ" target="_blank">https://groups.google.com/forum/#!searchin/cp2k/$26BS$20Marci/cp2k/8fTVlCEjSME/iUem9cK1-McJ</a></p><span class="ecxHOEnZb"><font color="#888888">

</font></span></div><span class="ecxHOEnZb"><font color="#888888">

<BR>

-- <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>
</font></span></blockquote></div><br></div>

<BR>

-- <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">cp2k+uns...@googlegroups.com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com">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></div>                                       </div></body>
</html>