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  , 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.
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:
- 1. Task → Single Point2. XC functional → Hybrid → PBE03. Basis set → TZ2P4. Frozen core → None
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:
- 1. In the panel bar select Properties → GW2. In the GW panel, check the Calculate GW quasi-particle energies checkbox3. Self consistency → G0W0
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:
- 1. File → Save and give the job an appropriate name, e.g. “Azulene”2. File → Run3. Wait for the calculation to finish (you can monitor the progress via the logfile, SCM → Logfile)
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:
- 1. SCM → Output2. In AMSoutput, search for “GW Quasi-Particle Spectrum”
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:
|||(1, 2, 3) 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)|