<div dir="ltr">Hi, Florian,<div><br></div><div>Thanks great for your suggestion.</div><div><br></div><div>I have tried broken symmetry for Ce with 'scf_guess' set as 'atomic'. But the situation is same to what I got from calculation without Ce (total energy ~0.86 eV higher than the restart situation)</div><div>Here are the parameters I used. Could you check if they look good? First time to use broken symmetry. Also, I will try to add it for O.</div><div>For two Ce atoms next-nearby the oxygen vacancy, I tried the following setting to make it Ce3+. </div><div><div> &BS</div><div> &ALPHA</div><div> NEL -2 -1 1</div><div> L 0 2 3</div><div> N 6 5 4</div><div> &END</div><div> &BETA</div><div> NEL -2 -1 -1</div><div> L 0 2 3</div><div> N 6 5 4</div><div> &END</div><div> &END</div></div><div>But for other Ce atoms, I used the following setting to make it Ce4+.</div><div><div> &BS</div><div> &ALPHA</div><div> NEL -2 -1 -1</div><div> L 0 2 3</div><div> N 6 5 4</div><div> &END</div><div> &BETA</div><div> NEL -2 -1 -1</div><div> L 0 2 3</div><div> N 6 5 4</div><div> &END</div><div> &END</div></div><div><br></div><div>Thanks!</div><div><br></div><div>Best.</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Apr 7, 2015 at 12:37 AM, Florian Schiffmann <span dir="ltr"><<a href="mailto:flosch...@gmail.com" target="_blank">flosch...@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">Hi,<br><br>that makes it more clear. As I am not 100% sure what is going on I will have to do some guesswork.<br><br>From your input you run an unrestricted calculation. In the atomic guess, CP2K still uses the restricted atomic orbitals as initial guess for both alpha and beta. This might lead to a situation where you have 0.5 of an electron initially sitting in the orbital affected by the U correction. As far as I remember te U correction is symmetric at this point and it might be arbitrary whether the orbital gets further filled or emptied due to the penalty function (therefore there is some arbitraryness in your result). I think this is somewhat general that the U term can lock you in a given state depending on your initial guess.<br><br>Could you try to use a broken symmetry guess in your calculation where you initialize O as 2- and Ce as 4+ and see how the energy is affected?<br>For O this chould look like <br>&BS<br>&ALPHA<br> L 1<br> N 2<br> NEL 2<br>&END<br> &BETA<br> L 1<br> N 2<br> NEL 2<br>&END<br>&END<br><br>For Ce it is a bit more complicated. simply remeber, that cp2k puts a closed shell wfn in alpha and beta. Therefore NEL in this case is 2 as the closed shell representation of a 2p3 is a 2p6 wfn.<br><br>Hope that helps<br><br>Flo <br></div><div class="HOEnZb"><div class="h5">
<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>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">Best,<div>Lei</div><div><br></div><div><span style="color:rgb(73,73,73);font-family:'Lucida Grande',Verdana,Arial;font-size:12px;line-height:16px">University of Nebraska - Lincoln,<br>Department of Chemistry,<br>536 </span><span style="color:rgb(73,73,73);font-family:'Lucida Grande',Verdana,Arial;font-size:11px;line-height:16px">Hamilton Hall.</span></div><div><span style="color:rgb(73,73,73);font-family:'Lucida Grande',Verdana,Arial;font-size:11px;line-height:16px">Lincoln, NE 68588</span></div><div><span style="color:rgb(73,73,73);font-family:'Lucida Grande',Verdana,Arial;font-size:11px;line-height:16px">Phone: 402-853-8119</span></div><div><br></div></div>
</div>