Question about XAS calculation code in VASP6.1

Question on input files/tags, interpreting output, etc.

Please check whether the answer to your question is given in the VASP online manual or has been discussed in this forum previously!

Moderators: Global Moderator, Moderator

Post Reply
Message
Author
zhaoru_sun1
Newbie
Newbie
Posts: 21
Joined: Wed Apr 08, 2020 2:12 pm

Question about XAS calculation code in VASP6.1

#1 Post by zhaoru_sun1 » Mon Aug 03, 2020 6:44 am

Dear everyone,

The XAS calculaion method of VASP6.1 described in VASPwiki is
Image
And the relevant code in src/core_con_mat.F is

Code: Select all

                                    EPSILON_HELP = PREFACTOR/FREQUENCY *                  &
                                                   WDES%WTKPT(KPOINT) * DELTAFUN *          &
                                                   MECC1 * CONJG(MECC2) * CH_AMPLIFICATION
"FREQUENCY" in the code should correspond to "ω" in the Wiki, but I am confused why not "FREQUENCY**2" like "ω^2" in wiki? What's wrong?

ferenc_karsai
Global Moderator
Global Moderator
Posts: 422
Joined: Mon Nov 04, 2019 12:44 pm

Re: Question about XAS calculation code in VASP6.1

#2 Post by ferenc_karsai » Mon Sep 14, 2020 8:20 am

Dear Zhaoru,

thank you very much for pointing out that problem. The equation is correct and unfortunately that was a bug in the code. It will be corrected in the next release.

Theoretically this could be more than a constant shift, but we tested the results for all of our calculations and fortunately only the overall intensities change visibly. Since one reports only arbitrary intensities and scales them all calculations until now should be correct.

zhaoru_sun1
Newbie
Newbie
Posts: 21
Joined: Wed Apr 08, 2020 2:12 pm

Re: Question about XAS calculation code in VASP6.1

#3 Post by zhaoru_sun1 » Sun Jul 11, 2021 9:14 am

ferenc_karsai wrote: Mon Sep 14, 2020 8:20 am thank you very much for pointing out that problem. The equation is correct and unfortunately that was a bug in the code. It will be corrected in the next release.
Thank you. I am reviewing the XAS theory. There is a problem with the relationship between the dielectric function and XAS that puzzles me.
In F. Karsai et al. Phys. Rev. B 98, 235205 (2018), the imaginary part of the dielectric function is:

Image

Image

While from EQ(3.10) in C. Brouder, J. Phys.: Condens. Matter 2, 701 (1990), the absorption cross section is:

Image

So I can roughly get:

Image

But this is not consistent with F. Karsai et al. Phys. Rev. B 98, 235205 (2018):
the transversal expression for the imaginary part of the dielectric function in the long-wavelength limit (q=0) is directly proportional to the absorption spectrum
This makes me confused.

ferenc_karsai
Global Moderator
Global Moderator
Posts: 422
Joined: Mon Nov 04, 2019 12:44 pm

Re: Question about XAS calculation code in VASP6.1

#4 Post by ferenc_karsai » Wed Jul 14, 2021 6:03 am

Your assessment look correct and it would have been better to write without the word "directly":

"the transversal expression for the imaginary part of the dielectric function in the long-wavelength limit (q=0) is proportional to the absorption spectrum"

zhaoru_sun1
Newbie
Newbie
Posts: 21
Joined: Wed Apr 08, 2020 2:12 pm

Re: Question about XAS calculation code in VASP6.1

#5 Post by zhaoru_sun1 » Wed Jul 14, 2021 3:51 pm

ferenc_karsai wrote: Wed Jul 14, 2021 6:03 am Your assessment look correct and it would have been better to write without the word "directly":

"the transversal expression for the imaginary part of the dielectric function in the long-wavelength limit (q=0) is proportional to the absorption spectrum"
Thank you. If so, for clarify I think the description in XANES_in_Diamond:
The XAS spectrum is the imaginary part of the frequency-dependent dielectric function
may need some modification and also the corresponding script plot_core_imdiel.sh.
Although the additional frequency factor has little effect on XAS within ~10 eV of the X-ray energy region.

ferenc_karsai
Global Moderator
Global Moderator
Posts: 422
Joined: Mon Nov 04, 2019 12:44 pm

Re: Question about XAS calculation code in VASP6.1

#6 Post by ferenc_karsai » Fri Jul 16, 2021 8:42 am

Thank you very much. Indeed the description of the tutorial is very sloppy. I will correct the tutorial and it's description.

zhaoru_sun1
Newbie
Newbie
Posts: 21
Joined: Wed Apr 08, 2020 2:12 pm

Re: Question about XAS calculation code in VASP6.1

#7 Post by zhaoru_sun1 » Tue Jul 20, 2021 6:01 am

ferenc_karsai wrote: Fri Jul 16, 2021 8:42 am
Thank you.
Is there any way for VASP to output the contributions of different unoccupied orbits in XAS calculation, that is, the discrete transition amplitudes before broadening. Just like what you did in Figure 2 and Figure 3 in your paper F. Karsai et al. Phys. Rev. B 98, 235205 (2018).

ferenc_karsai
Global Moderator
Global Moderator
Posts: 422
Joined: Mon Nov 04, 2019 12:44 pm

Re: Question about XAS calculation code in VASP6.1

#8 Post by ferenc_karsai » Mon Jul 26, 2021 4:31 am

Figure 2 is just an electron density plot, with the density calculated only in a given energy region (or for bands). Of course using the core-hole in the calculation. Please see the following for band decomposed charge densities:
wiki/index.php/Band_decomposed_charge_densities

Figure 3 is a BSE calculation and just shows the different directions of the diagonal elements of the dielectric function. Note that this calculation was done using WIEN2k.

There is no feature in VASP selecting only specific contributions from selected orbital quantum numbers. VASP always calculates all contributions in a given range.

Post Reply