Low rank approximation in G0W0 calculations
Abstract
The single particle energies obtained in a Kohn-Sham density functional theory (DFT) calculation are generally known to be poor approximations to electron excitation energies that are measured in tr ansport, tunneling and spectroscopic experiments such as photo-emission spectroscopy. The correction to these energies can be obtained from the poles of a single particle Green’s function derived from a many-body perturbation theory. From a computational perspective, the accuracy and efficiency of such an approach depends on how a self energy term that properly accounts for dynamic screening of electrons is approximated. The G0W0 approximation is a widely used technique in which the self energy is expressed as the convolution of a noninteracting Green’s function (G0) and a screened Coulomb interaction (W0) in the frequency domain. The computational cost associated with such a convolution is high due to the high complexity of evaluating W 0 at multiple frequencies. In this paper, we discuss how the cost of G0W0 calculation can be reduced by constructing a low rank approximation to the frequency dependent part of W 0 . In particular, we examine the effect of such a low rank approximation on the accuracy of the G0W0 approximation. We also discuss how the numericalmore »
- Authors:
-
- Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). Computational Research Division
- Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). Computational Research Division; Univ. of California, Berkeley, CA (United States). Dept. of Mathematics
- Central Univ. of Finance and Economics, Beijing (China). School of Statistics and Mathematics
- Univ. of California, Berkeley, CA (United States). Dept. of Physics; Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). Materials Science Division
- Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center (NERSC)
- Publication Date:
- Research Org.:
- Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR)
- OSTI Identifier:
- 1379538
- Grant/Contract Number:
- AC02-05CH11231; 11171232
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Science China Mathematics
- Additional Journal Information:
- Journal Volume: 59; Journal Issue: 8; Journal ID: ISSN 1674-7283
- Publisher:
- Science China Press - Springer
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 97 MATHEMATICS AND COMPUTING; density functional theory; G0W0 approximation; Sternheimer equation; contour deformation; low rank approximation
Citation Formats
Shao, MeiYue, Lin, Lin, Yang, Chao, Liu, Fang, Da Jornada, Felipe H., Deslippe, Jack, and Louie, Steven G. Low rank approximation in G0W0 calculations. United States: N. p., 2016.
Web. doi:10.1007/s11425-016-0296-x.
Shao, MeiYue, Lin, Lin, Yang, Chao, Liu, Fang, Da Jornada, Felipe H., Deslippe, Jack, & Louie, Steven G. Low rank approximation in G0W0 calculations. United States. https://doi.org/10.1007/s11425-016-0296-x
Shao, MeiYue, Lin, Lin, Yang, Chao, Liu, Fang, Da Jornada, Felipe H., Deslippe, Jack, and Louie, Steven G. Sat .
"Low rank approximation in G0W0 calculations". United States. https://doi.org/10.1007/s11425-016-0296-x. https://www.osti.gov/servlets/purl/1379538.
@article{osti_1379538,
title = {Low rank approximation in G0W0 calculations},
author = {Shao, MeiYue and Lin, Lin and Yang, Chao and Liu, Fang and Da Jornada, Felipe H. and Deslippe, Jack and Louie, Steven G.},
abstractNote = {The single particle energies obtained in a Kohn-Sham density functional theory (DFT) calculation are generally known to be poor approximations to electron excitation energies that are measured in tr ansport, tunneling and spectroscopic experiments such as photo-emission spectroscopy. The correction to these energies can be obtained from the poles of a single particle Green’s function derived from a many-body perturbation theory. From a computational perspective, the accuracy and efficiency of such an approach depends on how a self energy term that properly accounts for dynamic screening of electrons is approximated. The G0W0 approximation is a widely used technique in which the self energy is expressed as the convolution of a noninteracting Green’s function (G0) and a screened Coulomb interaction (W0) in the frequency domain. The computational cost associated with such a convolution is high due to the high complexity of evaluating W 0 at multiple frequencies. In this paper, we discuss how the cost of G0W0 calculation can be reduced by constructing a low rank approximation to the frequency dependent part of W 0 . In particular, we examine the effect of such a low rank approximation on the accuracy of the G0W0 approximation. We also discuss how the numerical convolution of G0 and W0 can be evaluated efficiently and accurately by using a contour deformation technique with an appropriate choice of the contour.},
doi = {10.1007/s11425-016-0296-x},
journal = {Science China Mathematics},
number = 8,
volume = 59,
place = {United States},
year = {Sat Jun 04 00:00:00 EDT 2016},
month = {Sat Jun 04 00:00:00 EDT 2016}
}
Web of Science
Works referenced in this record:
Accurate self-energies in a plane-wave basis using only a few empty states: Towards large systems
journal, August 2008
- Bruneval, Fabien; Gonze, Xavier
- Physical Review B, Vol. 78, Issue 8
Solution of Sparse Indefinite Systems of Linear Equations
journal, September 1975
- Paige, C. C.; Saunders, M. A.
- SIAM Journal on Numerical Analysis, Vol. 12, Issue 4
Applied Numerical Linear Algebra
book, January 1997
- Demmel, James W.
- Society for Industrial and Applied Mathematics
Numerical integration for ab initio many-electron self energy calculations within the GW approximation
journal, April 2015
- Liu, Fang; Lin, Lin; Vigil-Fowler, Derek
- Journal of Computational Physics, Vol. 286
The GW -Method for Quantum Chemistry Applications: Theory and Implementation
journal, December 2012
- van Setten, M. J.; Weigend, F.; Evers, F.
- Journal of Chemical Theory and Computation, Vol. 9, Issue 1
Phonons and related crystal properties from density-functional perturbation theory
journal, July 2001
- Baroni, Stefano; de Gironcoli, Stefano; Dal Corso, Andrea
- Reviews of Modern Physics, Vol. 73, Issue 2
Recycling Krylov Subspaces for Sequences of Linear Systems
journal, January 2006
- Parks, Michael L.; de Sturler, Eric; Mackey, Greg
- SIAM Journal on Scientific Computing, Vol. 28, Issue 5
Improving accuracy and efficiency of calculations of photoemission spectra within the many-body perturbation theory
journal, February 2012
- Nguyen, Huy-Viet; Pham, T. Anh; Rocca, Dario
- Physical Review B, Vol. 85, Issue 8
Toward the Optimal Preconditioned Eigensolver: Locally Optimal Block Preconditioned Conjugate Gradient Method
journal, January 2001
- Knyazev, Andrew V.
- SIAM Journal on Scientific Computing, Vol. 23, Issue 2
BerkeleyGW: A massively parallel computer package for the calculation of the quasiparticle and optical properties of materials and nanostructures
journal, June 2012
- Deslippe, Jack; Samsonidze, Georgy; Strubbe, David A.
- Computer Physics Communications, Vol. 183, Issue 6
New Method for Calculating the One-Particle Green's Function with Application to the Electron-Gas Problem
journal, August 1965
- Hedin, Lars
- Physical Review, Vol. 139, Issue 3A
Optical excitations in organic molecules, clusters, and defects studied by first-principles Green’s function methods
journal, May 2006
- Tiago, Murilo L.; Chelikowsky, James R.
- Physical Review B, Vol. 73, Issue 20
Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems
journal, March 1992
- van der Vorst, H. A.
- SIAM Journal on Scientific and Statistical Computing, Vol. 13, Issue 2
QMR: a quasi-minimal residual method for non-Hermitian linear systems
journal, December 1991
- Freund, Roland W.; Nachtigal, No�l M.
- Numerische Mathematik, Vol. 60, Issue 1
Resolution-of-identity approach to Hartree–Fock, hybrid density functionals, RPA, MP2 and GW with numeric atom-centered orbital basis functions
journal, May 2012
- Ren, Xinguo; Rinke, Patrick; Blum, Volker
- New Journal of Physics, Vol. 14, Issue 5
Krylov subspace recycling for sequences of shifted linear systems
journal, July 2014
- Soodhalter, Kirk M.; Szyld, Daniel B.; Xue, Fei
- Applied Numerical Mathematics, Vol. 81
The GW method
journal, March 1998
- Aryasetiawan, F.; Gunnarsson, O.
- Reports on Progress in Physics, Vol. 61, Issue 3
Electron correlation in semiconductors and insulators: Band gaps and quasiparticle energies
journal, October 1986
- Hybertsen, Mark S.; Louie, Steven G.
- Physical Review B, Vol. 34, Issue 8, p. 5390-5413
Self-energy operators and exchange-correlation potentials in semiconductors
journal, June 1988
- Godby, R. W.; Schlüter, M.; Sham, L. J.
- Physical Review B, Vol. 37, Issue 17
Implementation of an all-electron GW approximation based on the projector augmented wave method without plasmon pole approximation: Application to Si, SiC, AlAs, InAs, NaH, and KH
journal, April 2003
- Lebègue, S.; Arnaud, B.; Alouani, M.
- Physical Review B, Vol. 67, Issue 15
Solution of generalized shifted linear systems with complex symmetric matrices
journal, July 2012
- Sogabe, Tomohiro; Hoshi, Takeo; Zhang, Shao-Liang
- Journal of Computational Physics, Vol. 231, Issue 17
GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems
journal, July 1986
- Saad, Youcef; Schultz, Martin H.
- SIAM Journal on Scientific and Statistical Computing, Vol. 7, Issue 3
Time-Dependent Density Functional Response Theory for Molecules
book, November 1995
- Casida, Mark E.
- Recent Advances in Density Functional Methods
Ground and Low-Lying Excited States of Interacting Electron Systems: A Survey and Some Critical Analyses
book, December 1999
- Farid, Benham
- Electron Correlation in the Solid State
Green’s-function approach to linear response in solids
journal, May 1987
- Baroni, Stefano; Giannozzi, Paolo; Testa, Andrea
- Physical Review Letters, Vol. 58, Issue 18
KSSOLV—a MATLAB toolbox for solving the Kohn-Sham equations
journal, March 2009
- Yang, Chao; Meza, Juan C.; Lee, Byounghak
- ACM Transactions on Mathematical Software, Vol. 36, Issue 2
Electronic excitations in light absorbers for photoelectrochemical energy conversion: first principles calculations based on many body perturbation theory
journal, January 2013
- Ping, Yuan; Rocca, Dario; Galli, Giulia
- Chemical Society Reviews, Vol. 42, Issue 6
GW method with the self-consistent Sternheimer equation
journal, March 2010
- Giustino, Feliciano; Cohen, Marvin L.; Louie, Steven G.
- Physical Review B, Vol. 81, Issue 11
On the Early History of the Singular Value Decomposition
journal, December 1993
- Stewart, G. W.
- SIAM Review, Vol. 35, Issue 4
Large Scale GW Calculations
journal, May 2015
- Govoni, Marco; Galli, Giulia
- Journal of Chemical Theory and Computation, Vol. 11, Issue 6
GW quasiparticle spectra from occupied states only
journal, March 2010
- Umari, P.; Stenuit, Geoffrey; Baroni, Stefano
- Physical Review B, Vol. 81, Issue 11
Implementation and performance of the frequency-dependent method within the PAW framework
journal, July 2006
- Shishkin, M.; Kresse, G.
- Physical Review B, Vol. 74, Issue 3
An O ( N 3 ) implementation of Hedin's GW approximation for molecules
journal, August 2011
- Foerster, D.; Koval, P.; Sánchez-Portal, D.
- The Journal of Chemical Physics, Vol. 135, Issue 7
Works referencing / citing this record:
Identifying substitutional oxygen as a prolific point defect in monolayer transition metal dichalcogenides
journal, July 2019
- Barja, Sara; Refaely-Abramson, Sivan; Schuler, Bruno
- Nature Communications, Vol. 10, Issue 1