Page 1 of 1

How to read the result of GW+BSE calculation?

Posted: Wed Apr 09, 2025 7:57 am
by Jike.Wang

Hello everyone,
I am trying to calculate the optical properties of Al2O3 by GW+BSE. I already have the calculation result file(OUTCAR and vasprun.xml), but I am not sure where the relevant data is stored, so I would like to ask you for some information:
1. In which file is the data of transition energy stored? How to read the oscillator strength corresponding to the transition energy?
2. How to read the exciton binding energy?
3. How to read the dielectric function?


Re: How to read the result of GW+BSE calculation?

Posted: Wed Apr 09, 2025 11:10 am
by michael_wolloch

Hi Jike Wang,

There are a bunch of resources available that I can recommend for you to understand more about GW and BSE:

Now to your specific questions:

  1. In which file is the data of transition energy stored? How to read the oscillator strength corresponding to the transition energy?

You want to search your vasprun.xml file for "opticaltransitions". There you will find two columns. The first one corresponds to BSE eigenvalues, the second to the oscillator strength

  1. How to read the exciton binding energy?

You will have to subtract the BSE eigenvalues from the independent particle (IP) transition energies, both of which you can find in the BSEFATBAND file. A detailed tutorial on excitons analysis is available in exercise 10 of the BSE tutorials.

  1. How to read the dielectric function?

The frequency-dependent dielectric function is again stored in the vasprun.xml file. Search for "dielectricfunction" to find the tensor data. They are right below the optical transition data. The layout is: "energy, xx, yy, zz, xy, yz,zx"
Exercise 9.2 of the BSE tutorials explains how to compute the dielectric function. There also are some scripts that could help you with plotting.

Please let me know if this answers your questions,
Michael