<div dir="ltr">
<style type="text/css">p { margin-bottom: 0.25cm; line-height: 120%; }a:link { }</style><p style="margin-bottom: 0cm; line-height: 100%" lang="en-US">Dear
CP2K group,</p>
<p style="margin-bottom: 0cm; 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="margin-bottom: 0cm; line-height: 100%" lang="en-US">Taking
into account the scaling factor:
</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 8 <b>(because (10 + 8)/2 = 9 and I want 3d9, unpaired e-
spin 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 -10 <b>(because (10 - 10)/2 = 0 and I want 3d0, no e-) </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 -10 <b>(because (10 - 10)/2 = 0 and I want 3d0, no e-) </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 8 <b>(because (10 + 8)/2 = 9 and I want 3d9, unpaired e-
spin 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">The
results (attached) is a myriad of NaNs and 0.000000000</p>
<p style="margin-bottom: 0cm; line-height: 100%" lang="en-US">Any
ideas will be very much appreciated.</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">[1] -
https://groups.google.com/forum/#!searchin/cp2k/Broken$20symmetry/cp2k/ghhbtIxUzKk/uhwgb5PMTsgJ</p>
<p style="margin-bottom: 0cm; line-height: 100%" lang="en-US">[2] -
https://groups.google.com/forum/#!searchin/cp2k/$26BS$20Marci/cp2k/8fTVlCEjSME/iUem9cK1-McJ</p>
</div>