<div dir="ltr">Many years ago I did some tests of the method in that paper. I remember  that the value for Ueff one gets does depend a bit upon the parameters used, e.g. RKMAX, RMTs, k-points. This is not really a surprise, and the value of RKMAX 5 in the original paper was probably a compromise for speed; larger RKMAX will localize the electrons more inside the RMTs.<div><br></div><div>I would suggest that you use parameters representative of what you will use in your +U calculations for the estimation of the +U value. </div><div><br></div><div>And, always remember that what you are doing is estimating the value for an approximate correction of the exchange-correlation potential for use in your computational experiments. It is not truly ab-initio, it is a heuristic approach.</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Mar 27, 2017 at 8:20 AM, karima Physique <span dir="ltr"><<a href="mailto:physique.karima@gmail.com" target="_blank">physique.karima@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>
<div dir="ltr">
<div>Thank you for your answer,<br>
<br>
</div>
<span id="m_6451582177433692495gmail-result_box" class="m_6451582177433692495gmail-" lang="en"><span>The two calculations converged but with more than 100 iteration (bad convergence)</span><br>
 <span>But in the paper presented by </span></span><span id="m_6451582177433692495gmail-result_box" class="m_6451582177433692495gmail-" lang="en"><span>G. K. H. Madsen and P. Novak, they used a single point and a value of RKmax = 5,</span><br>
<span class="m_6451582177433692495gmail-">Can I force the calculation to 1000 Point for example and RKmax = 8.</span></span></div>
<div class="gmail_extra"><br>
<div class="gmail_quote">2017-03-27 8:16 GMT+02:00 Fecher, Gerhard <span dir="ltr">
<<a href="mailto:fecher@uni-mainz.de" target="_blank">fecher@uni-mainz.de</a>></span>:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
your calciulataion is not converged with respect to the number of k-points and/or plane waves, is it ?<br>
<br>
Ciao<br>
Gerhard<br>
<br>
DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy:<br>
"I think the problem, to be quite honest with you,<br>
is that you have never actually known what the question is."<br>
<br>
==============================<wbr>======<br>
Dr. Gerhard H. Fecher<br>
Institut of Inorganic and Analytical Chemistry<br>
Johannes Gutenberg - University<br>
55099 Mainz<br>
and<br>
Max Planck Institute for Chemical Physics of Solids<br>
01187 Dresden<br>
______________________________<wbr>__________<br>
Von: Wien [<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at" target="_blank">wien-bounces@zeus.theochem.tu<wbr>wien.ac.at</a>] im Auftrag von karima Physique [<a href="mailto:physique.karima@gmail.com" target="_blank">physique.karima@gmail.com</a>]<br>
Gesendet: Montag, 27. März 2017 02:09<br>
An: <a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">wien@zeus.theochem.tuwien.ac.a<wbr>t</a><br>
Betreff: [Wien] Calculation of Ueff<br>
<div>
<div class="m_6451582177433692495h5"><br>
Dear Wien2k users,<br>
<br>
<br>
I used the method described by G. K. H. Madsen and P. Novak but I found two different values of Ueff when I changed the number of K-point and RKmax<br>
<br>
For K=1 Point and RKmax 5 i found Ueff=3.6 eV<br>
but For K=100 Points and RKmax 8.5 i found Ueff=5.4 eV<br>
<br>
Do I use one point only ? otherwise I want to know why I had this difference<br>
</div>
</div>
______________________________<wbr>_________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.a<wbr>t</a><br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwMFaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=O7-03TkSPbowoFoVcd14aXi40g0l5mDGsFeVmYbVYV4&s=KijjZx91su9ZvIQFZTQ0LMuOo95YPqMYfws06uC2TXQ&e=" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac<wbr>.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at:  <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=DwMFaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=O7-03TkSPbowoFoVcd14aXi40g0l5mDGsFeVmYbVYV4&s=VVnzkIzm0smIF0UlCuqlN5OvFXUVo_w4_lLFF_q9Rb0&e=" rel="noreferrer" target="_blank">
http://www.mail-archive.com/wi<wbr>en@zeus.theochem.tuwien.ac.at/<wbr>index.html</a><br>
</blockquote>
</div>
<br>
</div>
</div>

</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><span style="font-size:12.8px">Professor Laurence Marks</span><br></div><div dir="ltr"><span style="font-size:12.8px">"Research is to see what everybody else has seen, and to think what nobody else has thought", </span><span style="font-size:12.8px">Albert Szent-Gyorgi</span><br><a href="http://www.numis.northwestern.edu" target="_blank">www.numis.northwestern.edu</a> ; <span style="font-size:12.8px">Corrosion in 4D: </span><a href="http://MURI4D.numis.northwestern.edu" style="font-size:12.8px" target="_blank">MURI4D.numis.northwestern.edu</a><div><span style="font-size:12.8px">Partner of the CFW 100% program for gender equity, </span><a href="http://www.cfw.org/100-percent" style="font-size:12.8px" target="_blank">www.cfw.org/100-percent</a></div><div>Co-Editor, Acta Cryst A</div></div></div></div></div></div></div></div></div></div></div>
</div>