<div dir="ltr">Thank you for your answer but i I noticed that Ueff is very sensitive with the change of RKmax and Kpoint it is for this I evoked this subject abd for example<div>For K = 1 and RKmax 5 i found Ueff = 3.6 eV</div><div>But For K = 100 Points and RKmax 8.5 i found Ueff = 5.4 eV</div><div><br></div><div>5.4 and 3.6 are relatively far </div><div><br></div><div>What I really ask is that the choice of one point is it mandatory? (as is indicated in the paper of <span style="color:rgb(80,0,80);font-size:12.8px">G. K. H. Madsen and P. Novak</span>) or else I have to force these parameters (a large number of Kpoint and rkmax> 7) which of these two choices is good?</div></div><div class="gmail_extra"><br><div class="gmail_quote">2017-03-27 17:08 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">you should increase the number of k points (and/or rmtkmax) until the U value does not change any longer within the limit that you set yourself (say 0.1 eV)<br>
that is what we call convergence.<br>
<br>
.. and as Laurence told, you should remember that you are doing an approximation to estimate some value<br>
<span class=""><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">wien-bounces@zeus.theochem.<wbr>tuwien.ac.at</a>] im Auftrag von karima Physique [<a href="mailto:physique.karima@gmail.com">physique.karima@gmail.com</a>]<br>
</span>Gesendet: Montag, 27. März 2017 16:53<br>
<span class="">An: A Mailing list for WIEN2k users<br>
Betreff: Re: [Wien] Calculation of Ueff<br>
<br>
</span><span class="">Indeed, the calculation converges quickly when I use 100 points with rkmax 8 with an iteration number that does not exceed 25 but with one point the number of iteration exceeds 100. In your opinion I take 100 points or one point as indicated in the paper of G. K. H. Madsen and P. Novak?<br>
<br>
thank you in advance<br>
<br>
</span>2017-03-27 16:25 GMT+02:00 Fecher, Gerhard <<a href="mailto:fecher@uni-mainz.de">fecher@uni-mainz.de</a><mailto:<a href="mailto:fecher@uni-mainz.de">fe<wbr>cher@uni-mainz.de</a>>>:<br>
<span class="">"The two calculations converged but with more than 100 iteration (bad convergence)"<br>
==> did you start them from scratch ? Try to run a calculation with the unchanged inputs (regular setting of electron occupation)<br>
then use the converged potential to run the calculations with changed input files, (I checked some of my old calculations, it never took more than 20-25 cycles)<br>
<br>
you also may increase the number of k-points (plane waves) and rerun the calculation with the preconverged charge densities.<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>
</span><span class="">Von: Wien [<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.<wbr>tuwien.ac.at</a><mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-<wbr>bounces@zeus.theochem.tuwien.<wbr>ac.at</a>>] im Auftrag von karima Physique [<a href="mailto:physique.karima@gmail.com">physique.karima@gmail.com</a><<wbr>mailto:<a href="mailto:physique.karima@gmail.com">physique.karima@gmail.<wbr>com</a>>]<br>
</span><span class="">Gesendet: Montag, 27. März 2017 15:41<br>
An: A Mailing list for WIEN2k users<br>
Betreff: Re: [Wien] Calculation of Ueff<br>
<br>
thank you very much Prof.Laurence Marks for your answer.<br>
<br>
</span>2017-03-27 15:33 GMT+02:00 Laurence Marks <<a href="mailto:L-marks@northwestern.edu">L-marks@northwestern.edu</a><<wbr>mailto:<a href="mailto:L-marks@northwestern.edu">L-marks@northwestern.<wbr>edu</a>><mailto:<a href="mailto:L-marks@northwestern.edu">L-marks@<wbr>northwestern.edu</a><mailto:<a href="mailto:L-marks@northwestern.edu">L-<wbr>marks@northwestern.edu</a>>>>:<br>
<span class="">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.<br>
<br>
I would suggest that you use parameters representative of what you will use in your +U calculations for the estimation of the +U value.<br>
<br>
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.<br>
<br>
</span><span class="">On Mon, Mar 27, 2017 at 8:20 AM, karima Physique <<a href="mailto:physique.karima@gmail.com">physique.karima@gmail.com</a><<wbr>mailto:<a href="mailto:physique.karima@gmail.com">physique.karima@gmail.<wbr>com</a>><mailto:<a href="mailto:physique.karima@gmail.com">physique.karima@<wbr>gmail.com</a><mailto:<a href="mailto:physique.karima@gmail.com">physique.<wbr>karima@gmail.com</a>>>> wrote:<br>
Thank you for your answer,<br>
<br>
The two calculations converged but with more than 100 iteration (bad convergence)<br>
But in the paper presented by G. K. H. Madsen and P. Novak, they used a single point and a value of RKmax = 5,<br>
Can I force the calculation to 1000 Point for example and RKmax = 8.<br>
<br>
</span>2017-03-27 8:16 GMT+02:00 Fecher, Gerhard <<a href="mailto:fecher@uni-mainz.de">fecher@uni-mainz.de</a><mailto:<a href="mailto:fecher@uni-mainz.de">fe<wbr>cher@uni-mainz.de</a>><mailto:<a href="mailto:fecher@uni-mainz.de">fech<wbr>er@uni-mainz.de</a><mailto:<a href="mailto:fecher@uni-mainz.de">fecher@<wbr>uni-mainz.de</a>>>>:<br>
<span class="">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>
</span>Von: Wien [<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.<wbr>tuwien.ac.at</a><mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-<wbr>bounces@zeus.theochem.tuwien.<wbr>ac.at</a>><mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@<wbr>zeus.theochem.tuwien.ac.at</a><<wbr>mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.<wbr>theochem.tuwien.ac.at</a>>>] im Auftrag von karima Physique [<a href="mailto:physique.karima@gmail.com">physique.karima@gmail.com</a><<wbr>mailto:<a href="mailto:physique.karima@gmail.com">physique.karima@gmail.<wbr>com</a>><mailto:<a href="mailto:physique.karima@gmail.com">physique.karima@<wbr>gmail.com</a><mailto:<a href="mailto:physique.karima@gmail.com">physique.<wbr>karima@gmail.com</a>>>]<br>
<span class="">Gesendet: Montag, 27. März 2017 02:09<br>
</span>An: <a href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.<wbr>at</a><mailto:<a href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.<wbr>tuwien.ac.at</a>><mailto:<a href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@<wbr>zeus.theochem.tuwien.ac.at</a><<wbr>mailto:<a href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.<wbr>tuwien.ac.at</a>>><br>
<span class="">Betreff: [Wien] Calculation of Ueff<br>
<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>
______________________________<wbr>_________________<br>
Wien mailing list<br>
</span><a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.<wbr>at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.<wbr>tuwien.ac.at</a>><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@<wbr>zeus.theochem.tuwien.ac.at</a><<wbr>mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.<wbr>tuwien.ac.at</a>>><br>
<span class=""><a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.<wbr>ac.at/mailman/listinfo/wien</a><<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">ht<wbr>tps://urldefense.proofpoint.<wbr>com/v2/url?u=http-3A__zeus.<wbr>theochem.tuwien.ac.at_mailman_<wbr>listinfo_wien&d=DwMFaQ&c=<wbr>yHlS04HhBraes5BQ9ueu5zKhE7rtNX<wbr>t_d012z2PA6ws&r=U_<wbr>T4PL6jwANfAy4rnxTj8IUxm818jnvq<wbr>KFdqWLwmqg0&m=O7-<wbr>03TkSPbowoFoVcd14aXi40g0l5mDGs<wbr>FeVmYbVYV4&s=<wbr>KijjZx91su9ZvIQFZTQ0LMuOo95YPq<wbr>MYfws06uC2TXQ&e=</a>><br>
SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/<wbr>wien@zeus.theochem.tuwien.ac.<wbr>at/index.html</a><<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">https://<wbr>urldefense.proofpoint.com/v2/<wbr>url?u=http-3A__www.mail-<wbr>2Darchive.com_wien-40zeus.<wbr>theochem.tuwien.ac.at_index.<wbr>html&d=DwMFaQ&c=<wbr>yHlS04HhBraes5BQ9ueu5zKhE7rtNX<wbr>t_d012z2PA6ws&r=U_<wbr>T4PL6jwANfAy4rnxTj8IUxm818jnvq<wbr>KFdqWLwmqg0&m=O7-<wbr>03TkSPbowoFoVcd14aXi40g0l5mDGs<wbr>FeVmYbVYV4&s=<wbr>VVnzkIzm0smIF0UlCuqlN5OvFXUVo_<wbr>w4_lLFF_q9Rb0&e=</a>><br>
<br>
<br>
<br>
<br>
--<br>
Professor Laurence Marks<br>
"Research is to see what everybody else has seen, and to think what nobody else has thought", Albert Szent-Gyorgi<br>
</span><a href="http://www.numis.northwestern.edu" rel="noreferrer" target="_blank">www.numis.northwestern.edu</a><<a href="http://www.numis.northwestern.edu" rel="noreferrer" target="_blank">htt<wbr>p://www.numis.northwestern.edu</a><wbr>><<a href="http://www.numis.northwestern.edu" rel="noreferrer" target="_blank">http://www.numis.<wbr>northwestern.edu</a>> ; Corrosion in 4D: <a href="http://MURI4D.numis.northwestern.edu" rel="noreferrer" target="_blank">MURI4D.numis.northwestern.edu</a><<a href="http://MURI4D.numis.northwestern.edu" rel="noreferrer" target="_blank"><wbr>http://MURI4D.numis.<wbr>northwestern.edu</a>><<a href="http://MURI4D.numis.northwestern.edu" rel="noreferrer" target="_blank">http://<wbr>MURI4D.numis.northwestern.edu</a>><br>
Partner of the CFW 100% program for gender equity, <a href="http://www.cfw.org/100-percent" rel="noreferrer" target="_blank">www.cfw.org/100-percent</a><<a href="http://www.cfw.org/100-percent" rel="noreferrer" target="_blank">http:/<wbr>/www.cfw.org/100-percent</a>><<a href="http://www.cfw.org/100-percent" rel="noreferrer" target="_blank">http<wbr>://www.cfw.org/100-percent</a>><br>
<span class="">Co-Editor, Acta Cryst A<br>
<br>
______________________________<wbr>_________________<br>
Wien mailing list<br>
</span><a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.<wbr>at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.<wbr>tuwien.ac.at</a>><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@<wbr>zeus.theochem.tuwien.ac.at</a><<wbr>mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.<wbr>tuwien.ac.at</a>>><br>
<div class="HOEnZb"><div class="h5"><a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.<wbr>ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/<wbr>wien@zeus.theochem.tuwien.ac.<wbr>at/index.html</a><br>
<br>
<br>
______________________________<wbr>_________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.<wbr>at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.<wbr>tuwien.ac.at</a>><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.<wbr>ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/<wbr>wien@zeus.theochem.tuwien.ac.<wbr>at/index.html</a><br>
<br>
______________________________<wbr>_________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.<wbr>at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.<wbr>ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/<wbr>wien@zeus.theochem.tuwien.ac.<wbr>at/index.html</a><br>
</div></div></blockquote></div><br></div>