GW: Ionization Potential and Electron Affinity¶
Ionization potential (IP) and electron affinity (EA) of molecules can be predicted with reasonable accuracy using the G0W0 method 1 2, especially if as starting point a (range separated) hybrid functional is used. Note that functionals that use XCfun (like LCY-BP86 and CAMY-B3LYP) can not be used as starting point.
One can also do eigenvalue-only self-consistent GW (evGW) with ADF, which is typically a factor 5 to 10 more expensive than G0W0, but which has the advantage that it is less starting-point dependent. Starting from AMS2022 one can also do quasiparticle self-consistent GW (qsGW) with ADF, which is more expensive than evGW, but the results of a qsGW are starting-point independent, and is usually the most accurate one.
Next in AMS2022 it also possible to go beyond the GW approximation for the self-energy and take into account the screened interaction to second order, with an even more expensive method called G3W2. See the ADF manual on GW for more details.
In this tutorial we will show how to compute the IP and EA of Azulene with ADF using G0W0@PBE0.
Recommendations for GW calculations
Video tip of the week
In AMS2023.102 (or bug fixed AMS2023 binary with r112249 or higher) a bug was fixed regarding GW. An incorrect Fermi energy was used, which also affects BSE. See also the ADF manual on GW.
The results presented here use the bug fixed binary, which affected the results presented in this Tutorial by less than 0.15 eV.
Set up and run the calculation¶
Let’s start by importing the structure of the Azulene molecule from our molecules database:
Let’s now set up the calculation details in the main panel:
Typically for EA one may need a larger basis set than TZ2P, see the section Recommendations for GW calculations in the ADF manual.
To enable the G0W0 calculation:
Use Self consistency → EVGW for evGW, and use Self consistency → QSGW for qsGW. Use Self energy → G3W2 for G3W2, typically used in combination with qsGW. In case of (any variant of) GW, in order to avoid numerical issues, the dependency key is set in ADF (and the RIHartreeFock%DependencyThreshold is increased), see the Recommendations for GW calculations in the ADF manual. One can change these dependency values using the panel bar Details → Accuracy.
We are now ready to run the calculation:
The results of the calculation are written both to the text output file and to the binary
To see the G0W0 results in the output file:
In AMS2023.102 (or bug fixed AMS2023 binary with r112249 or higher) a bug was fixed regarding GW. An incorrect Fermi energy was used, which also affects BSE. The results presented here use the bug fixed binary, which affected the results presented in this Tutorial by less than 0.15 eV. See also the ADF manual on GW.
The ionization potential (IP) equals the negative of the G0W0 quasi-particle HOMO energy:
The electron affinity (EA) equals the negative of the G0W0 quasi-particle LUMO energy:
These values are in good agreement with the CCSD(T) results from 1 (IP=7.55 [eV] and EA=0.54 [eV]) and with experimental results cited in 1 (IP=7.42-7.44 [eV] and EA=0.69-0.80 [eV]).
Joseph W. Knight, Xiaopeng Wang, Lukas Gallandi, Olga Dolgounitcheva, Xinguo Ren, J. Vincent Ortiz, Patrick Rinke, Thomas Körzdörfer, and Noa Marom, Accurate Ionization Potentials and Electron Affinities of Acceptor Molecules III: A Benchmark of GW Methods, J. Chem. Theory Comput. 2016, 12, 2, 615–626
Arno Förster, Lucas Visscher, Low-order scaling G0W0 by pair atomic density fitting, arXiv.org (2020)