<div dir="auto">My three cents. I think an agreement of 0.1eV should be considered as fortuitous. There are many issues which are glossed over even with the miraculous exact functional:<div dir="auto"><br></div><div dir="auto">1) The Slater method is a very clever use of the mean value theorem for an integral. However, it is only 1 value. You can check the literature, I remember seeing papers where people use a range of holes to more accurately do the integral.</div><div dir="auto"><br></div><div dir="auto">2) The simple dft-based calculations assume that the final states are plane waves. Rigorously the exiting photoelectron in XPS is an evanescent Bloch wave (for a crystal). There is literature on this, but I doubt that it has been combined with DFT.</div><div dir="auto"><br></div><div dir="auto">3) In experiments you have to worry about photoelectron diffraction, and there will be some shifts to higher apparent binding energy due to phonon inelastic scattering. And you have to worry about charging and band bending for insulators, chemisorption induced work function changes (how clean is your XPS?)....</div><div dir="auto"><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">www.numis.northwestern.edu</a></div></div></div><br><div class="gmail_quote"><div dir="ltr">On Thu, Jan 17, 2019, 02:29 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Sorry for the confusion. The quoted values in the pdf are probably from <br>
a lousy calculation and are ment only to demonstrate the effect (if I <br>
remember correctly, only a 4x bigger P supercell, but for sure not <br>
converged, I also don't remember which functional, ....) and only <br>
accidentally match experiment.<br>
<br>
My general experience is that core-eigenvalues (taken with respect to EF <br>
!!) are 10-15% off, while slater TS gives 1-2 %, i.e. an order of <br>
magnitude better.<br>
<br>
<br>
On 1/17/19 9:15 AM, Pavel Ondračka wrote:<br>
> Dear Wien2k mailing list,<br>
> <br>
> I'm looking for some advice regarding the calculation of core level<br>
> binding energies (to compare with XPS experiments). First of all there<br>
> is this nice lecture where prof. Blaha actually shows some calculations<br>
> <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__susi.theochem.tuwien.ac.at_reg-5Fuser_textbooks_WIEN2k-5Flecture-2D&d=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=tkgGA6DHWOQZpcOauDl3Lx8BZ_vGdV0W35tV5es6rQw&s=RWPHx3CucXJZx8K7anGE8WYmvRBxkjijDa1up4sQ7Rg&e=" rel="noreferrer noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__susi.theochem.tuwien.ac.at_reg-5Fuser_textbooks_WIEN2k-5Flecture-2D&d=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=tkgGA6DHWOQZpcOauDl3Lx8BZ_vGdV0W35tV5es6rQw&s=RWPHx3CucXJZx8K7anGE8WYmvRBxkjijDa1up4sQ7Rg&e=</a><br>
> notes_2011/Blaha_xas_eels.pdf of core levels with perfect results. For<br>
> example with TiN the deltaSCF method gets 397.1eV for the N1s level as<br>
> compared to 397.0eV experiment. The trouble is that I'm not able to<br>
> reproduce this.<br>
> <br>
> I've done some calculations before and I was never really happy with<br>
> the absolute values which were always few eV off but I've always<br>
> thought this is just the limitation of xc functional or methodology.<br>
> Hence seeing the nice results in the lecture surprised me. However, I'm<br>
> not able to reproduce the values even for metals from the example. For<br>
> the TiN I'm getting values of 404.8eV with the slaters transition state<br>
> approach and 404.6eV with delta-scf (here I'm using the formula for<br>
> metals E_b = E^tot_initialstate(N) - E^tot_finalstate(N), i.e. placing<br>
> the core-electron in the valence band and with PBE). I have thought<br>
> that this is maybe functional difference, since while taking LDA<br>
> instead of PBE shifts the results differ almost by 4eV (to 400.9eV).<br>
> However with the PBE I get the core energy ε_i as 377.4eV (consistent<br>
> with the mentioned pdf where it is 377.5eV) so maybe this is not just<br>
> about functional?. I've already checked convergence with supercell size<br>
> as well as numerical parameters and I'm actually out of ideas.<br>
> <br>
> To be honest, I'm not much concerned personally about the discrepancy<br>
> since the chemical shifts seem to be reasonable even if the absolute<br>
> values are not. I just think that if it is possible to get the absolute<br>
> values right (or at least closer to experiment) as in the lecture pdf,<br>
> the results would of course look way better, therefore I'll be grateful<br>
> for any comments and help.<br>
> <br>
> Best regards<br>
> Pavel<br>
> <br>
> <br>
> <br>
> <br>
> _______________________________________________<br>
> Wien mailing list<br>
> <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank" rel="noreferrer">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=tkgGA6DHWOQZpcOauDl3Lx8BZ_vGdV0W35tV5es6rQw&s=CdYQl6ILK-N4jpQvI_tObJBKSKSG7cLccPjdQ2OlXTk&e=" rel="noreferrer 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=tkgGA6DHWOQZpcOauDl3Lx8BZ_vGdV0W35tV5es6rQw&s=CdYQl6ILK-N4jpQvI_tObJBKSKSG7cLccPjdQ2OlXTk&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=tkgGA6DHWOQZpcOauDl3Lx8BZ_vGdV0W35tV5es6rQw&s=Jb9VWdpdCBzntaGLtD1ZQIA49GTWcrfpXPyEcA51d64&e=" rel="noreferrer 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=tkgGA6DHWOQZpcOauDl3Lx8BZ_vGdV0W35tV5es6rQw&s=Jb9VWdpdCBzntaGLtD1ZQIA49GTWcrfpXPyEcA51d64&e=</a><br>
> <br>
<br>
-- <br>
<br>
                                       P.Blaha<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" rel="noreferrer">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=tkgGA6DHWOQZpcOauDl3Lx8BZ_vGdV0W35tV5es6rQw&s=OuguOKewmJ4Pb2vzPpRwBjAB04HvTZ4MMCl22prGhgA&e=" rel="noreferrer noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.wien2k.at&d=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=tkgGA6DHWOQZpcOauDl3Lx8BZ_vGdV0W35tV5es6rQw&s=OuguOKewmJ4Pb2vzPpRwBjAB04HvTZ4MMCl22prGhgA&e=</a><br>
WWW:   <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.imc.tuwien.ac.at_TC-5FBlaha&d=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=tkgGA6DHWOQZpcOauDl3Lx8BZ_vGdV0W35tV5es6rQw&s=7TX18mPbP30USAQH5hslaat3rhIFwtyw2Jw6u2BeU0w&e=" rel="noreferrer noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.imc.tuwien.ac.at_TC-5FBlaha&d=DwIGaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=tkgGA6DHWOQZpcOauDl3Lx8BZ_vGdV0W35tV5es6rQw&s=7TX18mPbP30USAQH5hslaat3rhIFwtyw2Jw6u2BeU0w&e=</a><br>
--------------------------------------------------------------------------<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank" rel="noreferrer">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=tkgGA6DHWOQZpcOauDl3Lx8BZ_vGdV0W35tV5es6rQw&s=CdYQl6ILK-N4jpQvI_tObJBKSKSG7cLccPjdQ2OlXTk&e=" rel="noreferrer 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=tkgGA6DHWOQZpcOauDl3Lx8BZ_vGdV0W35tV5es6rQw&s=CdYQl6ILK-N4jpQvI_tObJBKSKSG7cLccPjdQ2OlXTk&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=tkgGA6DHWOQZpcOauDl3Lx8BZ_vGdV0W35tV5es6rQw&s=Jb9VWdpdCBzntaGLtD1ZQIA49GTWcrfpXPyEcA51d64&e=" rel="noreferrer 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=tkgGA6DHWOQZpcOauDl3Lx8BZ_vGdV0W35tV5es6rQw&s=Jb9VWdpdCBzntaGLtD1ZQIA49GTWcrfpXPyEcA51d64&e=</a><br>
</blockquote></div>