<div dir="auto">N.N.B., instead of a smaller ecut one can use .lcore. Sometimes this also improves convergence.<br><br><div data-smartmail="gmail_signature">_____<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><a href="http://www.numis.northwestern.edu">www.numis.northwestern.edu</a></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Nov 13, 2017 2:12 PM, "Laurence Marks" <<a href="mailto:L-marks@northwestern.edu">L-marks@northwestern.edu</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="auto"><span style="font-family:sans-serif">"besides the statement: be careful with the energy parameter (very fine E-search</span><span style="font-family:sans-serif"> in case.in1)"</span><div dir="auto"><font face="sans-serif">N.B.</font>, one does not need a fine search in the more recent Wien2k, only in the older versions.</div><div dir="auto"><br></div><div dir="auto">I seem to remember that GMax 20 mainly cancels out for LDA/GGA -- I would always use it (or more) for a meta GGA.<br><br><div data-smartmail="gmail_signature" dir="auto">_____<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><a href="http://www.numis.northwestern.edu" target="_blank">www.numis.northwestern.edu</a></div><div class="gmail_extra" dir="auto"><br><div class="gmail_quote">On Nov 13, 2017 2:00 PM, "Peter Blaha" <<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@theochem.tuwien.ac.at</a>> wrote:<br type="attribution"><blockquote class="m_-6420275147523541327quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="m_-6420275147523541327quoted-text">> 1- Computing the atomic energies of atoms like N and P in an FCC cell is<br>
> ok, however for O atom the high symmetry of the FCC cell results in 1/3<br>
> occupancies (for the 4th p electron of O) in the spin down case. Only<br>
> using a lower symmetry cell (orthorhombic) for O atom eliminates this issue.<br>
<br>
</div>This is true, however, for O this is a very small (negligible) energy<br>
difference. It is a bit more important (also to help convergence) for<br>
transition metals or even lanthanides.<br>
<div class="m_-6420275147523541327quoted-text"><br>
> 2- For very small RMTs (e.g. 1.1 for O) one needs a large value of Gmax<br>
> (~20).<br>
<br>
</div>Again true. GMAX for small spheres has some effect on E-tot. I have not<br>
cross checked if this truncation cancels out in the energy differences.<br>
<div class="m_-6420275147523541327quoted-text"><br>
> Using RMT = 1.1 for O in a compound will require one to use too small<br>
> RMTs for large atoms. For example if we use RMT = 1.1 for O in compounds<br>
> like SnO2, this will require one to use an RMT ~ 1.6 for Sn (RMTs<br>
> difference should not be > 30 %). This means that the ecut should be set<br>
> < -9.8 Ry which includes the 4s electrons to valence states. What are<br>
> the consequences?<br>
<br>
</div>Without test: I would "break" my own rules and set the Sn sphere still<br>
to 2.0 or 2.2 (such that no core leakage occurs at ecut=-6). I would not<br>
expect ghostbands coming up. This rule usually applies to cases where<br>
one atom could have very large RMT (2.5 or even larger) and probably<br>
converges rapidly (eg. alkali(earth) atoms in TM perovskites, ...)<br>
<br>
But in principle there is nothing to say against Sn 4s as valence,<br>
besides the statement: be careful with the energy parameter (very fine<br>
E-search in case.in1) and make sure that RKmax is sufficient for such<br>
localized states.<br>
<br>
<br>
><br>
> Best Regards,<br>
><br>
> Alay<br>
<div class="m_-6420275147523541327quoted-text">><br>
><br>
><br>
> ______________________________<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>
</div>> <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=AQ_8qeipj47Nvnx3rH6ZWBnmVNH79mLnft7Udb_swzA&s=zF60VSvgkhXPgAFuKbK1JGP3NGj05D16OuNwQfn_Biw&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.<wbr>com/v2/url?u=http-3A__zeus.the<wbr>ochem.tuwien.ac.at_mailman_lis<wbr>tinfo_wien&d=DwICAg&c=yHlS04Hh<wbr>Braes5BQ9ueu5zKhE7rtNXt_<wbr>d012z2PA6ws&r=U_T4PL6jwANfAy4r<wbr>nxTj8IUxm818jnvqKFdqWLwmqg0&m=<wbr>AQ_8qeipj47Nvnx3rH6ZWBnmVNH79m<wbr>Lnft7Udb_swzA&s=zF60VSvgkhXPgA<wbr>FuKbK1JGP3NGj05D16OuNwQfn_Biw&<wbr>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=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=AQ_8qeipj47Nvnx3rH6ZWBnmVNH79mLnft7Udb_swzA&s=HPtZ_3ky-3VP12W-LspS4XUBxVKHyqea_oR9uKgVGyE&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.<wbr>com/v2/url?u=http-3A__www.mail<wbr>-2Darchive.com_wien-40zeus.<wbr>theochem.tuwien.ac.at_index.<wbr>html&d=DwICAg&c=yHlS04HhBraes5<wbr>BQ9ueu5zKhE7rtNXt_d012z2PA6ws&<wbr>r=U_T4PL6jwANfAy4rnxTj8IUxm818<wbr>jnvqKFdqWLwmqg0&m=AQ_8qeipj47N<wbr>vnx3rH6ZWBnmVNH79mLnft7Udb_<wbr>swzA&s=HPtZ_3ky-3VP12W-LspS4XU<wbr>BxVKHyqea_oR9uKgVGyE&e=</a><br>
<div class="m_-6420275147523541327quoted-text">><br>
<br>
--<br>
------------------------------<wbr>------------------------------<wbr>--------------<br>
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
Phone: +43-1-58801-165300             FAX: +43-1-58801-165982<br>
</div>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=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=AQ_8qeipj47Nvnx3rH6ZWBnmVNH79mLnft7Udb_swzA&s=0xx31L8x_pX0tSE6QCwnu0ieVyBv4AN1-QfKUvUvjnk&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.<wbr>com/v2/url?u=http-3A__www.wien<wbr>2k.at&d=DwICAg&c=yHlS04HhBraes<wbr>5BQ9ueu5zKhE7rtNXt_<wbr>d012z2PA6ws&r=U_T4PL6jwANfAy4r<wbr>nxTj8IUxm818jnvqKFdqWLwmqg0&m=<wbr>AQ_8qeipj47Nvnx3rH6ZWBnmVNH79m<wbr>Lnft7Udb_swzA&s=0xx31L8x_pX0tS<wbr>E6QCwnu0ieVyBv4AN1-QfKUvUvjnk&<wbr>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=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=AQ_8qeipj47Nvnx3rH6ZWBnmVNH79mLnft7Udb_swzA&s=GjL51v86wzTeoP28vT7izbE-GVvMlSxptO2qD8t58To&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.<wbr>com/v2/url?u=http-3A__www.imc.<wbr>tuwien.ac.at_tc-5Fblaha-2D-2D-<wbr>2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-<wbr>2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-<wbr>2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-<wbr>2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-<wbr>2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-<wbr>2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-<wbr>2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-<wbr>2D&d=DwICAg&c=yHlS04HhBraes5BQ<wbr>9ueu5zKhE7rtNXt_d012z2PA6ws&r=<wbr>U_T4PL6jwANfAy4rnxTj8IUxm818jn<wbr>vqKFdqWLwmqg0&m=AQ_8qeipj47Nvn<wbr>x3rH6ZWBnmVNH79mLnft7Udb_swzA&<wbr>s=GjL51v86wzTeoP28vT7izbE-GVvM<wbr>lSxptO2qD8t58To&e=</a><br>
<div class="m_-6420275147523541327quoted-text"><br>
______________________________<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>
</div><a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=AQ_8qeipj47Nvnx3rH6ZWBnmVNH79mLnft7Udb_swzA&s=zF60VSvgkhXPgAFuKbK1JGP3NGj05D16OuNwQfn_Biw&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.<wbr>com/v2/url?u=http-3A__zeus.the<wbr>ochem.tuwien.ac.at_mailman_lis<wbr>tinfo_wien&d=DwICAg&c=yHlS04Hh<wbr>Braes5BQ9ueu5zKhE7rtNXt_<wbr>d012z2PA6ws&r=U_T4PL6jwANfAy4r<wbr>nxTj8IUxm818jnvqKFdqWLwmqg0&m=<wbr>AQ_8qeipj47Nvnx3rH6ZWBnmVNH79m<wbr>Lnft7Udb_swzA&s=zF60VSvgkhXPgA<wbr>FuKbK1JGP3NGj05D16OuNwQfn_Biw&<wbr>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=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=AQ_8qeipj47Nvnx3rH6ZWBnmVNH79mLnft7Udb_swzA&s=HPtZ_3ky-3VP12W-LspS4XUBxVKHyqea_oR9uKgVGyE&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.<wbr>com/v2/url?u=http-3A__www.mail<wbr>-2Darchive.com_wien-40zeus.<wbr>theochem.tuwien.ac.at_index.<wbr>html&d=DwICAg&c=yHlS04HhBraes5<wbr>BQ9ueu5zKhE7rtNXt_d012z2PA6ws&<wbr>r=U_T4PL6jwANfAy4rnxTj8IUxm818<wbr>jnvqKFdqWLwmqg0&m=AQ_8qeipj47N<wbr>vnx3rH6ZWBnmVNH79mLnft7Udb_<wbr>swzA&s=HPtZ_3ky-3VP12W-LspS4XU<wbr>BxVKHyqea_oR9uKgVGyE&e=</a><br>
</blockquote></div><br></div></div></div>
</blockquote></div></div>