Efficient dielectric matrix calculations using the Lanczos algorithm for fast many-body $G_0W_0$ implementations

We present a $G_0W_0$ implementation that assesses the two major bottlenecks of traditional plane-waves implementations, the summations over conduction states and the inversion of the dielectric matrix, without introducing new approximations in the formalism. The first bottleneck is circumvented by converting the summations into Sternheimer equations. Then, the novel avenue of expressing the dielectric matrix in a Lanczos basis is developed, which reduces the matrix size by orders of magnitude while being computationally efficient. We also develop a model dielectric operator that allows us to further reduce the size of the dielectric matrix without accuracy loss. Furthermore, we develop a scheme that reduces the numerical cost of the contour deformation technique to the level of the lightest plasmon pole model. Finally, the use of the simplified quasi-minimal residual scheme in replacement of the conjugate gradients algorithm allows a direct evaluation of the $G_0W_0$ corrections at the desired real frequencies, without need for analytical continuation. The performance of the resulting $G_0W_0$ implementation is demonstrated by comparison with a traditional plane-waves implementation, which reveals a 500-fold speedup for the silane molecule. Finally, the accuracy of our $G_0W_0$ implementation is demonstrated by comparison with other $G_0W_0$ calculations and experimental results.

[1]  Matthias Krack,et al.  Pseudopotentials for H to Kr optimized for gradient-corrected exchange-correlation functionals , 2005 .

[2]  Fabien Bruneval,et al.  Benchmarking the Starting Points of the GW Approximation for Molecules. , 2013, Journal of chemical theory and computation.

[3]  M. L. Tiago,et al.  Optical excitations in organic molecules, clusters, and defects studied by first-principles Green's function methods , 2006, cond-mat/0605248.

[4]  Ab initio calculations of electronic excitations: Collapsing spectral sums , 2010 .

[5]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[6]  S. Louie,et al.  Simple approximate physical orbitals for GW quasiparticle calculations. , 2011, Physical review letters.

[7]  Huy V. Nguyen,et al.  Improving accuracy and efficiency of calculations of photoemission spectra within the many-body perturbation theory , 2012 .

[8]  Lai‐Sheng Wang,et al.  High-resolution photoelectron imaging of cold C₆₀⁻ anions and accurate determination of the electron affinity of C₆₀. , 2014, The Journal of chemical physics.

[9]  L. Curtiss,et al.  Assessment of Gaussian-2 and density functional theories for the computation of enthalpies of formation , 1997 .

[10]  Matthieu Verstraete,et al.  First-principles computation of material properties: the ABINIT software project , 2002 .

[11]  D. Sánchez-Portal,et al.  An O(N3) implementation of Hedin's GW approximation for molecules. , 2011, The Journal of chemical physics.

[12]  A. Alavi,et al.  Efficient calculation of the exact exchange energy in periodic systems using a truncated Coulomb potential , 2008 .

[13]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[14]  D. Turnbull,et al.  Solid State Physics : Advances in Research and Applications , 1978 .

[15]  L. Reining,et al.  Efficient calculation of the polarizability: a simplified effective-energy technique , 2012 .

[16]  Michiel J. van Setten,et al.  The GW-Method for Quantum Chemistry Applications: Theory and Implementation. , 2013, Journal of chemical theory and computation.

[17]  K. Capelle,et al.  A bird's-eye view of density-functional theory , 2002 .

[18]  Stefan Goedecker,et al.  Fast Radix 2, 3, 4, and 5 Kernels for Fast Fourier Transformations on Computers with Overlapping Multiply-Add Instructions , 1997, SIAM J. Sci. Comput..

[19]  Kosmas Prassides,et al.  Crystal structure and bonding of ordered C60 , 1991, Nature.

[20]  W. Kohn,et al.  Self-Consistent Equations Including Exchange and Correlation Effects , 1965 .

[21]  V. Anisimov.  Strong Coulomb Correlations in Electronic Structure Calculations , 2000 .

[22]  M. Schlüter,et al.  Self-energy operators and exchange-correlation potentials in semiconductors. , 1988, Physical review. B, Condensed matter.

[23]  ELIMINATION OF UNOCCUPIED-STATE SUMMATIONS IN AB INITIO SELF-ENERGY CALCULATIONS FOR LARGE SUPERCELLS , 1997, cond-mat/9803228.

[24]  L. Steinbeck,et al.  The GW space-time method for the self-energy of large systems , 1999 .

[25]  Richard M. Martin.  Electronic Structure: Frontmatter , 2004 .

[26]  Friedhelm Bechstedt,et al.  Molecular electronic excitations calculated from a solid-state approach: Methodology and numerics , 2005 .

[27]  J. Neaton,et al.  Electronic energy level alignment at metal-molecule interfaces with a G W approach , 2011, 1111.2569.

[28]  Stefano de Gironcoli,et al.  Phonons and related crystal properties from density-functional perturbation theory , 2000, cond-mat/0012092.

[29]  G. Rignanese,et al.  Band offsets at the Si/SiO2 interface from many-body perturbation theory. , 2008, Physical review letters.

[30]  Stefan Goedecker,et al.  ABINIT: First-principles approach to material and nanosystem properties , 2009, Comput. Phys. Commun..

[31]  Stephen L. Adler,et al.  Quantum theory of the dielectric constant in real solids. , 1962 .

[32]  Andreas Frommer,et al.  BiCGStab(ℓ) for Families of Shifted Linear Systems , 2003, Computing.

[33]  S. Louie,et al.  GW method with the self-consistent Sternheimer equation , 2009, 0912.3087.

[34]  Franccois Bottin,et al.  Large scale ab initio calculations based on three levels of parallelization , 2007, 0707.3405.

[35]  M. Saunders,et al.  Solution of Sparse Indefinite Systems of Linear Equations , 1975 .

[36]  Teter,et al.  Separable dual-space Gaussian pseudopotentials. , 1996, Physical review. B, Condensed matter.

[37]  R. Needs,et al.  Metal-insulator transition in Kohn-Sham theory and quasiparticle theory. , 1989, Physical review letters.

[38]  Claudio Attaccalite,et al.  First-principles GW calculations for fullerenes, porphyrins, phtalocyanine, and other molecules of interest for organic photovoltaic applications , 2010, 1011.3933.

[39]  Kenneth W. Nebesny,et al.  Valence and core photoelectron spectroscopy of C60, buckminsterfullerene , 1991 .

[40]  Nathan Wiser,et al.  Dielectric Constant with Local Field Effects Included , 1963 .

[41]  S. Louie,et al.  Coulomb-hole summations and energies for GW calculations with limited number of empty orbitals: a modified static remainder approach , 2012, 1208.0266.

[42]  R. Gebauer,et al.  Efficient approach to time-dependent density-functional perturbation theory for optical spectroscopy. , 2005, Physical review letters.

[43]  W. Marsden.  I and J , 2012 .

[44]  W. Linden,et al.  Precise quasiparticle energies and Hartree-Fock bands of semiconductors and insulators. , 1988, Physical review. B, Condensed matter.

[45]  F. Gygi,et al.  Iterative calculations of dielectric eigenvalue spectra , 2009 .

[46]  Yasutake Toyoshima,et al.  Vacuum ultraviolet absorption cross sections of SiH4, GeH4, Si2H6, and Si3H8 , 1986 .

[47]  G. Rignanese,et al.  Electronic properties of interfaces and defects from many‐body perturbation theory: Recent developments and applications , 2011 .

[48]  Andrew G. Glen,et al.  APPL , 2001 .

[49]  R. Freund,et al.  Software for simplified Lanczos and QMR algorithms , 1995 .

[50]  B. Lundqvist,et al.  Single-particle spectrum of the degenerate electron gas , 1968 .

[51]  S. Blugel,et al.  Band convergence and linearization error correction of all-electron GW calculations: The extreme case of zinc oxide , 2011, 1102.3255.

[52]  S. Louie,et al.  Quasiparticle band gap of ZnO: high accuracy from the conventional G⁰W⁰ approach. , 2010, Physical review letters.

[53]  X. Gonze,et al.  Adiabatic density-functional perturbation theory. , 1995, Physical review. A, Atomic, molecular, and optical physics.

[54]  Reinhold Schneider,et al.  Daubechies wavelets as a basis set for density functional pseudopotential calculations. , 2008, The Journal of chemical physics.

[55]  L. Reining,et al.  Electronic excitations: density-functional versus many-body Green's-function approaches , 2002 .

[56]  Xavier Gonze,et al.  A brief introduction to the ABINIT software package , 2005 .

[57]  L. Hedin.  NEW METHOD FOR CALCULATING THE ONE-PARTICLE GREEN'S FUNCTION WITH APPLICATION TO THE ELECTRON-GAS PROBLEM , 1965 .

[58]  Srivastava,et al.  Electronic structure , 2001, Physics Subject Headings (PhySH).

[59]  Louie,et al.  Electron correlation in semiconductors and insulators: Band gaps and quasiparticle energies. , 1986, Physical review. B, Condensed matter.

[60]  Engel,et al.  Generalized plasmon-pole model and plasmon band structures of crystals. , 1993, Physical review. B, Condensed matter.

[61]  M. L. Cohen,et al.  High accuracy many-body calculational approaches for excitations in molecules. , 2000, Physical review letters.

[62]  S. Baroni,et al.  GW quasiparticle spectra from occupied states only , 2009, 0910.0791.

[63]  Y. Saad,et al.  Turbo charging time-dependent density-functional theory with Lanczos chains. , 2006, The Journal of chemical physics.

[64]  Stefan Goedecker,et al.  Efficient solution of Poisson's equation with free boundary conditions. , 2006, The Journal of chemical physics.

[65]  F. Giustino,et al.  Ab initio Sternheimer-GW method for quasiparticle calculations using plane waves , 2013 .

[66]  X. Gonze,et al.  Accurate GW self-energies in a plane-wave basis using only a few empty states: Towards large systems , 2008 .

[67]  A. Tkatchenko,et al.  Resolution-of-identity approach to Hartree–Fock, hybrid density functionals, RPA, MP2 and GW with numeric atom-centered orbital basis functions , 2012, 1201.0655.

[68]  F. Gygi,et al.  Efficient iterative method for calculations of dielectric matrices , 2008 .

[69]  R. Roberge,et al.  The far ultraviolet and HeI photoelectron spectra of alkyl and fluorine substituted silane derivatives , 1978 .

[70]  K. Sankaran,et al.  G 0 W 0 band gap of ZnO: Effects of plasmon-pole models , 2011 .

[71]  J. Neaton,et al.  Quantitative molecular orbital energies within a G0W0 approximation , 2012, 1204.0509.

[72]  Katharina Wagner,et al.  Solid State Physics Advances In Research And Applications , 2016 .

[73]  Robert H. Halstead,et al.  Matrix Computations , 2011, Encyclopedia of Parallel Computing.

[74]  Steven G. Louie,et al.  Excitonic Effects and the Optical Absorption Spectrum of Hydrogenated Si Clusters , 1998 .

[75]  R. Sternheimer,et al.  ELECTRONIC POLARIZABILITIES OF IONS FROM THE HARTREE-FOCK WAVE FUNCTIONS , 1954 .

[76]  Stefano Baroni,et al.  Optimal representation of the polarization propagator for large-scale GW calculations , 2009 .

[77]  Huy V. Nguyen,et al.  GW calculations using the spectral decomposition of the dielectric matrix: Verification, validation, and comparison of methods , 2013 .