<div dir="ltr">Thanks great, Matt!<div><br></div><div>I will continue work on the calibration and see which method is most suitable in my case.</div><div><br></div><div><br></div><div>Best,</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Apr 7, 2015 at 2:50 PM, Matt W <span dir="ltr"><<a href="mailto:MattWa...@gmail.com" target="_blank">MattWa...@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"><div>Hi again,</div><span class=""><div> </div><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div></div><div>Thanks for your great suggestion. </div><div><br></div><div>Core correction indeed can help localize electron density correctly. But next step I will run molecular dynamics simulation with adsorbed molecules. In that case, I will not be able to use the core correction.</div><div><br></div></div></blockquote><div> </div></span><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div></div><div>So question is: if the electronic properties can be correctly predicted by restarting from the wavefunction obtained from the calculation with core correction, are the parameters (basis set, potential, U value, etc.) good enough to be used in MD simulation?</div><div><br></div></div></blockquote><div><br></div><div>sorry can't answer that. This would allow you to systematically get different solutions and compare them. To be clear, I only suggest this "hack" with the core correction to get an initial guess - then you remove it. At that point, when the correction is removed, you get what your density functional gives you.</div><div><br></div><div>For MD you need to decide a starting point - then what wants to happen -> happens, if you see what I mean, e.g. <a href="http://journals.aps.org/prl/abstract/10.1103/PhysRevLett.105.146405" target="_blank">http://journals.aps.org/prl/abstract/10.1103/PhysRevLett.105.146405</a></div><div><br></div><div>Whether the basis sets/pseudopotentials/density functional you are using are good enough to describe what you want ...you need to decide.</div><div><br></div><div>Matt</div><div><br></div><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div></div><div>Sorry for so many questions since I don't want to start a big project with some unknown questions. </div><div><br></div><div>Thanks great!</div></div><div><div class="h5"><div><br><div class="gmail_quote">On Tue, Apr 7, 2015 at 12:10 PM, Matt W <span dir="ltr"><<a rel="nofollow">Matt...@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"><div>Hi,</div><div><br></div>you should check in the atomic output at the beginning of the run whether the changes to the atomic states have worked as you want - it shows the orbital occupations of both spin channels for each atomic kind. Also, because of .. well it is a feature .. the changes in electron number will be divided by two (there is a post from Marcella about this somewhere).<br><br>I'd still go with what I suggested a while back - add a CORE_CORRECTION to your desired Ce3+ to make the electrons localize, then restart with that wavefunction without the extra core charge. <div><br></div><div>Run in a triplet (or maximum multiplicity for all aligned Ce3+) state - there should be no noticable exchange interaction between pairs of Ce3+, if you get splittings between singlet/triplet spin states etc then something is wrong. Running in max spin state automatically breaks some symmetry and makes it easier to get a proper solution.<div><br></div><div>Matt<div><div><br><br>On Tuesday, April 7, 2015 at 6:00:17 PM UTC+2, stclair wrote:</div></div><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><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></div><div><br><div class="gmail_quote"><span>On Tue, Apr 7, 2015 at 12:37 AM, Florian Schiffmann <span dir="ltr"><<a rel="nofollow">flos...@gmail.com</a>></span> wrote:<br></span><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span><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></span><div><div><span>

<p></p>

-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br></span>
To unsubscribe from this group and stop receiving emails from it, send an email to <a rel="nofollow">cp2k+uns...@googlegroups.<u></u>com</a>.<br>
To post to this group, send email to <a rel="nofollow">cp...@googlegroups.com</a>.<span><br>
Visit this group at <a href="http://groups.google.com/group/cp2k" rel="nofollow" target="_blank">http://groups.google.com/<u></u>group<u></u>/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" rel="nofollow" target="_blank">https://groups.google.com/d/<u></u>op<u></u>tout</a>.<br>
</span></div></div></blockquote></div><span><br><br clear="all"><div><br></div>-- <br><div>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: <a value="+14028538119">402-853-8119</a></span></div><div><br></div></div>
</span></div>
</blockquote></div></div></div><div><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 rel="nofollow">cp2k+uns...@googlegroups.<u></u>com</a>.<br>
To post to this group, send email to <a rel="nofollow">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="http://groups.google.com/group/cp2k" rel="nofollow" target="_blank">http://groups.google.com/<u></u>group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" rel="nofollow" target="_blank">https://groups.google.com/d/<u></u>optout</a>.<br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div>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: <a href="tel:402-853-8119" value="+14028538119" target="_blank">402-853-8119</a></span></div><div><br></div></div>
</div>
</div></div></blockquote></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>