GW approximation of Hedin's equations: Difference between revisions

From VASP Wiki
No edit summary
Line 100: Line 100:


== Single Shot: G<sub>0</sub>W<sub>0</sub> ==  
== Single Shot: G<sub>0</sub>W<sub>0</sub> ==  
Performing only one GW iteration step is commonly referred to the G<sub>0</sub>W<sub>0</sub> method. Here the self-energy <math>\Sigma^{(0)}</math> is determined and the eigenvalue equation<ref name="HybertsenLouie"/>
<span id="EVG0W0">
<math>
(T+V_{ext}+V_h)\phi_{n{\bf k}}({\bf r})+\int d{\bf r}\Sigma^{(0)}\left({\bf r},{\bf r}',\omega=E^{(0)}_{n{\bf k}}\right)\phi_{n{\bf k}}({\bf r}') = E^{(0)}_{n{\bf k}}\phi_{n{\bf k}}({\bf r})
</math>
</span>
is solved for the quasi-particle energies <math>E^{(0)}_{n\mathbf{k}}</math>. The orbitals, in contrast, remain unchanged.
Although self-consistency is completely neglected, G<sub>0</sub>W<sub>0</sub> often yields satisfactory results with band-gaps close to experimental measurements.<ref name="shishkin-PRB74"/><ref name="shishkin-PRB75"/>


== Partially self-consistent: GW<sub>0</sub> or EVGW<sub>0</sub> ==
== Partially self-consistent: GW<sub>0</sub> or EVGW<sub>0</sub> ==


== Self-consistent Quasi-particle approximation: QPGW<sub>0</sub> ==
== Self-consistent Quasi-particle approximation: QPGW<sub>0</sub> ==
Line 114: Line 126:
<ref name="shishkin:prl:07">[http://link.aps.org/doi/10.1103/PhysRevLett.99.246403 M. Shishkin, M. Marsman, and G. Kresse, Phys. Rev. Lett. 99, 246403 (2007).]</ref>
<ref name="shishkin:prl:07">[http://link.aps.org/doi/10.1103/PhysRevLett.99.246403 M. Shishkin, M. Marsman, and G. Kresse, Phys. Rev. Lett. 99, 246403 (2007).]</ref>
<ref name="rojas">[https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.74.1827 H. N. Rojas, R. W. Godby and R. J. Needs, Phys. Rev. Lett. 74, 1827 (1995)]</ref>
<ref name="rojas">[https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.74.1827 H. N. Rojas, R. W. Godby and R. J. Needs, Phys. Rev. Lett. 74, 1827 (1995)]</ref>
<ref name="shishkin-PRB74">[http://link.aps.org/doi/10.1103/PhysRevB.74.035101 M. Shishkin and G. Kresse, Phys. Rev. B 74, 035101 (2006).]</ref>
<ref name="shishkin-PRB75">[http://link.aps.org/doi/10.1103/PhysRevB.75.235102 M. Shishkin and G. Kresse, Phys. Rev. B 75, 235102 (2007).]</ref>
</references>
</references>


----
----
[[Category:Many-Body Perturbation Theory]][[Category:Theory]]
[[Category:Many-Body Perturbation Theory]][[Category:Theory]]

Revision as of 10:13, 24 July 2019

Green's functions

The GW method can be understood in terms of the following eigenvalue equation[1]

Here is the kinetic energy, the external potential of the nuclei, the Hartree potential and the quasiparticle energies with orbitals . In contrast to DFT, the exchange-correlation potential is replaced by the many-body self-energy and should be obtained together with the Green's function , the irreducible polarizability , the screened Coulomb interaction and the irreducible vertex function in a self-consistent procedure. For completeness, these equations are[2]

Here the common notation was adopted and denotes the bare Coulomb interaction. Note, that these equations are exact and provide an alternative to the Schrödinger equation for the many-body problem. Nevertheless, approximations are necessary for realistic systems. The most popular one is the GW approximation and is obtained by neglecting the equation for the vertex function and using the bare vertex instead:

This means that the equations for the polarizability and self-energy reduce to

while the equations for the Green's function and the screened potential remain the same.

However, in practice, these equations are usually solved in reciprocal space in the frequency domain

In principle Hedin's equations have to be solved self-consistently, where in the first iteration is the non-interacting Green's function

with being a set of one-electron orbitals and the corresponding energies. Afterwards the polarizability is determined, followed by the screened potential and the self-energy . This means that GW calculations require a first guess for the one-electron eigensystem, which is usually taken from a preceding DFT step.

In principle, one has to repeat all steps by the updating the Green's function with the Dyson equation given above in each iteration cycle until self-consistency is reached. In practice, this is hardly ever done due to computational complexity on the one hand (in fact fully self-consistent GW calculations are available as of VASP 6 only).

On the other hand, one observes that by keeping the screened potential in the first iteration to the DFT level one benefits from error cancelling,[3] which is the reason why often the screening is kept on the DFT level and one aims at self-consistency in Green's function only.

Following possible approaches are applied in practice.

Single Shot: G0W0

Performing only one GW iteration step is commonly referred to the G0W0 method. Here the self-energy is determined and the eigenvalue equation[1]

is solved for the quasi-particle energies . The orbitals, in contrast, remain unchanged.

Although self-consistency is completely neglected, G0W0 often yields satisfactory results with band-gaps close to experimental measurements.[4][5]

Partially self-consistent: GW0 or EVGW0

Self-consistent Quasi-particle approximation: QPGW0

Low-scaling GW: The Space-time Formalism

The scaling with system size (number of electrons) of GW calculations can be reduced[6] by performing a so-called Wick-rotation to imaginary time .

References