The unit of the absorption spectrum calculated by BSE

[Guo_BIT]

Dear developers:
I am confused about alpha, the unit of the absorption spectrum.
I obtained the absorption spectrum by solving the BSE with Coulomb cutoff technique(o-2D_WR_WC.alpha_q1_diago_bse).
By modifying the value of BDmRange, the absorption spectrum values have changed significantly. Therefore, I don't know how to convert this result to other units, such as those shown in Figure 3 of this article：
https://journals.aps.org/prb/abstract/1 ... 102.155119

Best,

Jingda Guo

[Daniele Varsano]

Dear Jingda,
BDmRange governs the width of the Lorentzian functions used to simulate the spectrum. Note that as results, BSE provides excitation energies and strengths, so the meaningful values are the oscillator strengths and not the values of the peaks that of course increase when reducing the damping values.

In the paper you refer, it is reported the imaginary part of the dielectric function (which is unitless), but for a non-periodic system (e.g. 2D) this quantity is not well-defined as it depends on the volume of the used supercell and goes to (1,0) for large vacuum (i.e. isolated 2D system). This is explained e.g. in Cudazzo et al. PRB 84, 085406 (2011) (see Eq. 16 and discussion).

Anyway, you can compare your results with the one published in that paper as the shape and relative weight of the peaks are not affected as alpha is proportional to epsilon, and it is well-defined.

Alternatively you can avoid using the coulomb cutoff, and look at the epsilon, using the same cell used by the authors of the paper (if reported), note that convergence with vacuum will be slower.

Best,
Daniele

[Guo_BIT]

Dear Daniele:
Thank you for patiently explaining