<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">Dear Anup,</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">Getting DFT to come close is going to be a challenge. You may find DOI: 10.1103/PhysRevMaterials.2.025001 useful. Some general comments:</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">1) In my experience +U is too harsh for 4f electrons. In your case 6eV is not that far below E_F and may well be in the same energy range as states for other elements (e.g. O). +U will force the occupied 4f low, the unoccupied high and this may not be appropriate.</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">2) The width of the valence region will be too small, because you are not including core hole effects.</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">3) The tabulated cross-sections for 4f  (e.g. NIST) are horribly bad -- this is know experimentally but I do not know of any better calculations.</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">4) There is the problem of Spin-Orbit. One also has to worry about whether the correct parity of states is in the DFT, it is easy to get the wrong ones at least in the rare earth scandates.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Oct 20, 2019 at 1:49 PM 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 Prof. Laurence,<div><br></div><div>Thank you for your question. I have performed some magnetic measurements on the sample for which I am doing DFT calculations, which suggests Sm to be in Sm+3 state. At the same time I have performed ARPES measurements on the single crystalline material, and from the ARPES measurements (which I have performed) and the previous ARPES report on this sample it seems that there are no Sm 4f states near the Fermi energy. This also suggests that Sm exist in Sm+3 state. Now since Sm+3 state has an energy far below the Fermi level around -6 eV and only Sm +2 states lie near the Fermi level (JPS Conf. Proc. 3, 017038 (2014). So, I wanted to reproduce the band structure by treating Sm 4f as core states and compare the DFT results  and the ARPES results. The other reason is that when I am doing GGA+U calculations using various values of U, the DFT results are not matching with the observed ARPES spectra. Looking forward to hearing from you. </div><div><br></div><div>SIncerely, <div><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 </font></font></div></div></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=uVSeK0Qr8Kl3rdemU6WJuYxDEzFZOL0opAcdg9goxe0&s=l6hhlfacgvmaZaMhbLnEaFzV72c9l0ZidWpCm4SWCME&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=uVSeK0Qr8Kl3rdemU6WJuYxDEzFZOL0opAcdg9goxe0&s=l6hhlfacgvmaZaMhbLnEaFzV72c9l0ZidWpCm4SWCME&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=uVSeK0Qr8Kl3rdemU6WJuYxDEzFZOL0opAcdg9goxe0&s=DycLY_6qn0ps4Mt1fUJW1RYEVMItG9kUtP-RPtQ3Yz0&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=uVSeK0Qr8Kl3rdemU6WJuYxDEzFZOL0opAcdg9goxe0&s=DycLY_6qn0ps4Mt1fUJW1RYEVMItG9kUtP-RPtQ3Yz0&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>