<div dir="ltr"><div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">Dear Viktor,</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)">Let me add something about large RMTs. I have noticed in the past that, in general,  the gradient at the RMT calculated from both sides (x RMTCheck) is very similar, e.g. </div><div class="gmail_default" style="font-family:verdana,sans-serif;color:rgb(0,0,0)"><br></div><div class="gmail_default"><font color="#000000">Atom   1 O    | RMT Charge   0.489 Grad   0.690 | Step Charge  0.00136, 0.00000 Gradient   0.2655, -0.2655 O</font><br></div><div class="gmail_default"><font color="#000000"><br></font></div><div class="gmail_default"><font color="#000000">The last two numbers should sum to almost zero. However, with very large RMTs they can deviate from this, particularly when the gradient is small. I think this leads to anomalies at the RMT. (Peter would say that the linearization is breaking down with an inadequate basis set, or something similar.) Reducing the RMTs generally avoids this.</font></div><div class="gmail_default"><font color="#000000"><br></font></div><div class="gmail_default"><font color="#000000">While it is true that larger RMTs are "better" in terms of speed for a single iteration, slightly smaller RMTs which may take longer for a single iteration but are more stable so require fewer iterations can be nett faster. Speed is not the same as stability and convergence, and often they oppose each other. (Think about the tradeoff with condenser aperture size: smaller is more coherent but more drift...)</font></div><div class="gmail_default"><font color="#000000"><br></font></div><div class="gmail_default"><font color="#000000">N.B., are you really sure that the Nd distances are that large?</font></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jan 17, 2019 at 3:31 PM Peter Blaha <<a href="mailto:pblaha@theochem.tuwien.ac.at">pblaha@theochem.tuwien.ac.at</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">RMTs:  Make sure that you choose RMTs, which can be used in all <br>
elemental solids, but the same should also be used in your compound.<br>
<br>
Yes, you can use very large RMTs if you follow the recommendations <br>
during init_lapw:<br>
Use HDLOs for f and (as L.Marks mentioned) probably d.<br>
Use lvns=6 (or 8) in case.in1<br>
<br>
For a difficult system (and Nd with partially occupied 5f electrons is <br>
definitely a difficult systems), scf convergence can be more difficult. <br>
Thus, as mentioned before, use both -ec and -cc (at least 0.001; <br>
eventually a better cc is difficult to reach. Often you may need MORE <br>
than the default max of 40 scf cycles. So when runsp_lapw stops, make <br>
sure it does not stop beucause of the 40 it limits.<br>
Clearly, your energies with :dis=0.2 are wrong.<br>
<br>
For the compound but also the elements, make sure to use "consistent <br>
RKMAX" values.<br>
If you use: Al Ti Nd (2.3,2.5,2.7) and for your compound you use <br>
RKMAX=8, than   Al should be done with RKMAX=8, Ti with RKmax=8*2.5/2.3, <br>
Nd with 8*2.7/2.3<br>
<br>
Am 17.01.2019 um 21:19 schrieb Victor Zenou:<br>
> Dear Lauri<br>
> Thanks for your answer<br>
> I used only energy convergence criteria (10^-4). Still the charge <br>
> convergence was between 0.006 and 0.2 e.<br>
> The nearest neighbor distance is 6.856 a.u, I used RMT=2.7 a.u., summed <br>
> to 5.4 which is much lower than 6.856. Usually large RMT is preferred <br>
> to make the calculations as fast as possible. I plan to check a phases <br>
> in Al-Ti-Nd and I chose RMT equal to 2.3, 2.5 and 2.7 respectively.<br>
> Victor<br>
> <br>
> <br>
> <br>
> ‫בתאריך יום ה׳, 17 בינו׳ 2019 ב-11:16 מאת ‪Laurence Marks‬‏ <br>
> <‪<a href="mailto:L-marks@northwestern.edu" target="_blank">L-marks@northwestern.edu</a> <mailto:<a href="mailto:L-marks@northwestern.edu" target="_blank">L-marks@northwestern.edu</a>>‬‏>:‬<br>
> <br>
>     With such large RMT you certainly need HDLO for Nd, perhaps both d &<br>
>     f. I would not want to use such large RMTs.<br>
> <br>
>     Have you checked that the charge convergence is good?<br>
> <br>
>     _____<br>
>     Professor Laurence Marks<br>
>     "Research is to see what everybody else has seen, and to think what<br>
>     nobody else has thought", Albert Szent-Gyorgi<br>
>     <a href="http://www.numis.northwestern.edu" rel="noreferrer" target="_blank">http://www.numis.northwestern.edu</a> <<a href="http://www.numis.northwestern.edu" rel="noreferrer" target="_blank">http://www.numis.northwestern.edu</a>><br>
> <br>
>     On Thu, Jan 17, 2019, 03:04 Victor Zenou <<a href="mailto:zanov@post.bgu.ac.il" target="_blank">zanov@post.bgu.ac.il</a><br>
>     <mailto:<a href="mailto:zanov@post.bgu.ac.il" target="_blank">zanov@post.bgu.ac.il</a>> wrote:<br>
> <br>
>         Hi<br>
>         Here is updated question:<br>
> <br>
>         Dear Wien2k users<br>
> <br>
>         I'm using wien2k version 17.1 installed on Ubunto 18.04.<br>
> <br>
>         I need to calculate the enthalpy of formation of few<br>
>         intermetallic phases that include Nd. For that element I used<br>
>         RMT=2.7 a.u. and RMT*Kmax=8.5 for spin-polarization calculations<br>
>         via GGA-PBE. I didn't try spin-orbit coupling. Yet.<br>
> <br>
>         The first stage is to choose proper k-points, which will give<br>
>         accurate results (I used 1x10^-4 Ryd for energy convergence) on<br>
>         one hand, but on the other hand won't be expensive in terms of<br>
>         computing time.<br>
> <br>
>         I noticed that there are large energy fluctuations (1x10^-2 Ryd;<br>
>         which are 2 order of magnitude higher than the energy accuracy<br>
>         (defined by convergence criterion), even when I went to 15,000<br>
>         k-points. Here is a list of energy as function of k-points:<br>
> <br>
>         K-points          IBZ      E [Ryd]<br>
> <br>
>           100       8        -77040.4692<br>
> <br>
>           500     38        -77040.4780<br>
> <br>
>         1000    81        -77040.5062<br>
> <br>
>         2000    120      -77040.5061<br>
> <br>
>         3000    208      -77040.4391<br>
> <br>
>         4000    244      -77040.4699<br>
> <br>
>         5000    280      -77040.4878<br>
> <br>
>         7500    455      -77040.4707<br>
> <br>
>         10000  540      -77040.4881<br>
> <br>
>         15000  840      -77040.4694<br>
> <br>
>         I would be happy to get an idea what could have gone wrong.<br>
> <br>
>         Best regards, Victor<br>
> <br>
> <br>
> <br>
>         ‫בתאריך יום ד׳, 16 בינו׳ 2019 ב-13:59 מאת ‪Victor Zenou‬‏<br>
>         <‪<a href="mailto:zanov@post.bgu.ac.il" target="_blank">zanov@post.bgu.ac.il</a> <mailto:<a href="mailto:zanov@post.bgu.ac.il" target="_blank">zanov@post.bgu.ac.il</a>>‬‏>:‬<br>
> <br>
>             Dear Wien2k users<br>
> <br>
>             I'm using wien2k version 17.1 installed on Ubunto 18.04.<br>
> <br>
>             As a part of investigation of intermetallic phases in<br>
>             Al-Ti-Nd system, I need to calculate the enthalpy of<br>
>             formation of few phases in that system. For that purpose, I<br>
>             start calculating enthalpy (equal to energy at 0K) for each<br>
>             element and phase at their optimal relaxed state. I used<br>
>             spin-polarization calculations. (GGA-PBE) As advised in<br>
>             Wien2k website I used the same RMT for each elements: 2.3,<br>
>             2.5 and 2.7 a.u. for Al, Ti and Nd, respectively. Also<br>
>             proper RMT*Kmax (7.5 to Ti, 8.5 to Nd and 7.0 for the rest).<br>
> <br>
>             The first stage is to choose proper k-points, which give you<br>
>             accurate results (I used 1x10^-4 Ryd for energy convergence)<br>
>             on one hand, but also some save computing time.<br>
> <br>
>             The calculation for all elements and phases went OK, but Nd.<br>
>             From certain k-points the energy fluctuated between 2 values<br>
>             (differed by 1x10^-2). Attach a list of energy as function<br>
>             of k-points:<br>
> <br>
>             K-pointsIBZ E [Ryd]<br>
> <br>
>             1008-77040.4692<br>
> <br>
>             50038-77040.4780<br>
> <br>
>             100081-77040.5062<br>
> <br>
>             2000120-77040.5061<br>
> <br>
>             3000208-77040.4391<br>
> <br>
>             4000244-77040.4699<br>
> <br>
>             5000280-77040.4878<br>
> <br>
>             7500455-77040.4707<br>
> <br>
>             10000540-77040.4881<br>
> <br>
>             15000840-77040.4694<br>
> <br>
>             I tried to get smart and did volume optimization using 4000<br>
>             k-points, later 5000 k-points. I got strange results.<br>
> <br>
>             I would be happy to get an idea what could have gone wrong.<br>
> <br>
>             I didn't try spin-orbit coupling<br>
> <br>
>             Best regards, Victor<br>
> <br>
>     _______________________________________________<br>
>     Wien mailing list<br>
>     <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a> <mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>><br>
>     <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=eK9aqsUgRkfDy2bZdzRW1MQXAt1xlSX8wLRVAOOVSLk&s=-x4eR4rPflykk4gZARg89aJRk1lrjzh4KhQaak_JByQ&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=eK9aqsUgRkfDy2bZdzRW1MQXAt1xlSX8wLRVAOOVSLk&s=-x4eR4rPflykk4gZARg89aJRk1lrjzh4KhQaak_JByQ&e=</a><br>
>     SEARCH the MAILING-LIST at:<br>
>     <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=eK9aqsUgRkfDy2bZdzRW1MQXAt1xlSX8wLRVAOOVSLk&s=aP3iF7UUFlgDwQ3l1Ge7GzYpkF2bqtTRb4xZCL-kRhM&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=eK9aqsUgRkfDy2bZdzRW1MQXAt1xlSX8wLRVAOOVSLk&s=aP3iF7UUFlgDwQ3l1Ge7GzYpkF2bqtTRb4xZCL-kRhM&e=</a><br>
> <br>
> <br>
> _______________________________________________<br>
> Wien mailing list<br>
> <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
> <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=eK9aqsUgRkfDy2bZdzRW1MQXAt1xlSX8wLRVAOOVSLk&s=-x4eR4rPflykk4gZARg89aJRk1lrjzh4KhQaak_JByQ&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=eK9aqsUgRkfDy2bZdzRW1MQXAt1xlSX8wLRVAOOVSLk&s=-x4eR4rPflykk4gZARg89aJRk1lrjzh4KhQaak_JByQ&e=</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=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=eK9aqsUgRkfDy2bZdzRW1MQXAt1xlSX8wLRVAOOVSLk&s=aP3iF7UUFlgDwQ3l1Ge7GzYpkF2bqtTRb4xZCL-kRhM&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=eK9aqsUgRkfDy2bZdzRW1MQXAt1xlSX8wLRVAOOVSLk&s=aP3iF7UUFlgDwQ3l1Ge7GzYpkF2bqtTRb4xZCL-kRhM&e=</a><br>
> <br>
<br>
-- <br>
--------------------------------------------------------------------------<br>
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
Phone: +43-1-58801-165300             FAX: +43-1-58801-165982<br>
Email: <a href="mailto:blaha@theochem.tuwien.ac.at" target="_blank">blaha@theochem.tuwien.ac.at</a>    WIEN2k: <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.wien2k.at&d=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=eK9aqsUgRkfDy2bZdzRW1MQXAt1xlSX8wLRVAOOVSLk&s=yscbD2URDxKhBSNrsUvCXygszqazvEKzTMM4mptFXqY&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.wien2k.at&d=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=eK9aqsUgRkfDy2bZdzRW1MQXAt1xlSX8wLRVAOOVSLk&s=yscbD2URDxKhBSNrsUvCXygszqazvEKzTMM4mptFXqY&e=</a><br>
WWW: <br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.imc.tuwien.ac.at_tc-5Fblaha-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D&d=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=eK9aqsUgRkfDy2bZdzRW1MQXAt1xlSX8wLRVAOOVSLk&s=fi0SsUcgDnugWH41GhhT5XyalNBmpXl8Mik8Y2ueAkI&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.imc.tuwien.ac.at_tc-5Fblaha-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D&d=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=eK9aqsUgRkfDy2bZdzRW1MQXAt1xlSX8wLRVAOOVSLk&s=fi0SsUcgDnugWH41GhhT5XyalNBmpXl8Mik8Y2ueAkI&e=</a> <br>
<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=eK9aqsUgRkfDy2bZdzRW1MQXAt1xlSX8wLRVAOOVSLk&s=-x4eR4rPflykk4gZARg89aJRk1lrjzh4KhQaak_JByQ&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=eK9aqsUgRkfDy2bZdzRW1MQXAt1xlSX8wLRVAOOVSLk&s=-x4eR4rPflykk4gZARg89aJRk1lrjzh4KhQaak_JByQ&e=</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=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=eK9aqsUgRkfDy2bZdzRW1MQXAt1xlSX8wLRVAOOVSLk&s=aP3iF7UUFlgDwQ3l1Ge7GzYpkF2bqtTRb4xZCL-kRhM&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=eK9aqsUgRkfDy2bZdzRW1MQXAt1xlSX8wLRVAOOVSLk&s=aP3iF7UUFlgDwQ3l1Ge7GzYpkF2bqtTRb4xZCL-kRhM&e=</a><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="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>