Page 1 of 1
Phonon-Assisted Absorption Mismatch
Posted: Sat Aug 23, 2025 2:46 pm
by Agathyan
Dear all,
I’m studying Si optical absorption. Used scissor flag (0, 0.44 eV) in EPW instead of BGW eigenvalues; results don’t match tutorial. Attached: inputs, band structure, absorption plots. Kindly help me to match the results with the scissor flag.
Attached
1. Bandstructure
- Si_BS_comparison.jpg (138.7 KiB) Viewed 2391 times
2. Absorption
- Si_Abs_comparison.jpg (67 KiB) Viewed 2391 times
3. Inputs
Re: Phonon-Assisted Absorption Mismatch
Posted: Tue Aug 26, 2025 8:06 pm
by stiwari
Hello,
When you apply a scissor shift, you need to renormalize the oscillator strengths. This is not performed automatically in EPW and is an additional step you need to perform. Check page 11 of
https://journals.aps.org/prb/abstract/1 ... 109.195127 for details.
In a nutshell, you need to multiply your imaginary dielectric constant with ((omega/(omega-delta))^2; here omega is the photon frequency and delta is the scissor shift.
Best regards,
Sabya.
Re: Phonon-Assisted Absorption Mismatch
Posted: Wed Aug 27, 2025 6:29 pm
by xiaozha
Hi,
Just as to add to Sabya's point: this renormalization is performed in the case of reading an external eigenvalue. Therefore, it is mostly recommended to use that option as much as possible instead of using a scissor. If you would like use a scissor shift, one way to have EPW does the renormalization is to supply a .eig file where the scissor shift is already applied.
The other thing to be careful with is to make sure in either case (eig_read or a scissor), the Fermi level is correctly set. It looks possible that you Fermi level is set to something very close to the VBM or CBM in the scissor case of your plot. If that's the case, you may see some sub-gap contributions like what you are seeing in ~0.5 eV because of intraband absorption due to partially occupied bands.
Best!
Xiao
Re: Phonon-Assisted Absorption Mismatch
Posted: Fri Aug 29, 2025 12:21 pm
by Agathyan
Dear Sabya and Xiao,
Thank you so much for your valuable suggestions!
1. I learned about the necessity of renormalisation, but I couldn't understand how it was implemented from the paper. In Python, I tried renormalising and got the subsequent result (which I don't think is the proper way to do it).
If possible, could you tell or guide me with literature on how to find oscillator strengths?
Code: Select all
hcut = 6.582e-16
c = 3e10
n_r = 3.4
scissor = 0.4403
#Before Renormalisation
alpha = (data[:,0] / (hcut * c * n_r)) * data[:,7]
#After Renormalisation
alpha_renom = (data[:,0] / (hcut * c * n_r)) * data[:,7] * (data[:,0]/(data[:,0]-scissor))**2
- Scissor_forum_approach.jpg (116.66 KiB) Viewed 2294 times
2. As Xiao said, I had used the VBM value as the Fermi Energy. On shifting it properly (mid of VBM and CBM), I got the following result: close to the BGW eig file's absorption.
- Si_BS_scissor.jpg (125.46 KiB) Viewed 2294 times
- fermishift_forum_approach.jpg (132.54 KiB) Viewed 2294 times