<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">I will ask why you consider it critical to put 4f electrons in the core. While there are some cases when they are isolated, there are others when this is not the case. By forcing them into the core you are deciding what the results should be which is very, very bad science.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Oct 20, 2019 at 9:47 AM Anup Shakya <<a href="mailto:npshakya31@gmail.com">npshakya31@gmail.com</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"><div dir="ltr"><div>Dear Xavier,</div><div><br></div><div>Thank you for the suggestion. I have actually performed GGA+U calculations using various values of U.  But even with U = 6eV, Sm 4f states lies at an energy of -0.6 eV which shifts a bit to -0.7 eV when I increase the U to 10 eV.  There are some papers in the literature where people have performed the calculations by keeping the rare earth 4f states in the core region using VASP.  A recent example for Pr is shown here in this paper (Phys. Rev. B <b>99</b>, 245131). I had previously performed calculations by treating rare earth 4f states in the core region using VASP which is pretty easy as there is already a defined pseudo potential for this specific purpose only. But since I don't have access to VASP right now I want to know whether this can be done with Wien2k also or not? Any suggestions would be much appreciated.</div><div><br></div><div>Sincerely, <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div dir="ltr" style="font-size:12.8px"><font color="#333333"><font face="tahoma, sans-serif">Anup Pradhan Sakhya<br></font></font></div></div></div></div></div></div></div>
_______________________________________________<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=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=7tXWAtMh0obwafTByO6OYcHQsytLI1lVv5X-6FOwAo4&s=pIhpbkK3wFYsjCLhGSInL17oWR5iMoo58DBp3sxPK6A&e=" rel="noreferrer" target="_blank">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=7tXWAtMh0obwafTByO6OYcHQsytLI1lVv5X-6FOwAo4&s=pIhpbkK3wFYsjCLhGSInL17oWR5iMoo58DBp3sxPK6A&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=7tXWAtMh0obwafTByO6OYcHQsytLI1lVv5X-6FOwAo4&s=PbCvY1xk3r-VnJVt0d9BUEn-9fr078I_SpNnjLrIFUo&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=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=7tXWAtMh0obwafTByO6OYcHQsytLI1lVv5X-6FOwAo4&s=PbCvY1xk3r-VnJVt0d9BUEn-9fr078I_SpNnjLrIFUo&e=</a> <br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr">Professor Laurence Marks<br>Department of Materials Science and Engineering<br>Northwestern University<br><a href="http://www.numis.northwestern.edu/" target="_blank">www.numis.northwestern.edu</a><div>Corrosion in 4D: <a href="http://www.numis.northwestern.edu/MURI" target="_blank">www.numis.northwestern.edu/MURI</a><br>Co-Editor, Acta Cryst A<br>"Research is to see what everybody else has seen, and to think what nobody else has thought"<br>Albert Szent-Gyorgi</div></div></div>