Efficient dielectric matrix calculations using the Lanczos algorithm for fast many-body implementations
arXiv:1408.3036 · doi:10.1103/PhysRevB.91.125120
Abstract
We present a 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 corrections at the desired real frequencies, without need for analytical continuation. The performance of the resulting 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 implementation is demonstrated by comparison with other calculations and experimental results.
19 pages, 2 figures
References in corpus (21)
- Phonons and related properties of extended systems from density-functional perturbation theory
- Separable Dual Space Gaussian Pseudo-potentials
- Resolution-of-identity approach to Hartree-Fock, hybrid density functionals, RPA, MP2, and \textit{GW} with numeric atom-centered orbital basis functions
- First-principles GW calculations for fullerenes, porphyrins, phtalocyanine, and other molecules of interest for organic photovoltaic applications
- Daubechies wavelets as a basis set for density functional pseudopotential calculations
- Turbo charging time-dependent density-functional theory with Lanczos chains
- A new and efficient approach to time-dependent density-functional perturbation theory for optical spectroscopy
- Efficient solution of Poisson's equation with free boundary conditions
- Optical excitations in organic molecules, clusters and defects studied by first-principles Green's function methods
- GW quasi-particle spectra from occupied states only
- High Accuracy Many-Body Calculational Approaches for Excitations in Molecules
- Coulomb-hole summations and energies for GW calculations with limited number of empty orbitals: a modified static remainder approach
- The GW space-time method for the self-energy of large systems
- Large scale ab initio calculations based on three levels of parallelization
- Band convergence and linearization error correction of all-electron GW calculations: The extreme case of zinc oxide
- Band Offsets at the Si/SiO Interface from Many-Body Perturbation Theory
- GW method with the self-consistent Sternheimer equation
- An O(N^3) implementation of Hedin's GW approximation for molecules
- Electronic energy level alignment at metal-molecule interfaces with a GW approach
- Quantitative Molecular Orbital Energies within a Approximation
- Elimination of unoccupied state summations in ab-initio self-energy calculations for large supercells
Cited by in corpus (9)
- Quasi-Particle Self-Consistent for Molecules
- Towards Fully Converged GW Calculations for Large Systems
- All-electron Gaussian-based for Valence and Core Excitation Energies of Periodic Systems
- Accurate description of charged excitations in molecular solids from embedded many-body perturbation theory
- GPU Acceleration of Large-Scale Full-Frequency GW Calculations
- Critical assessment of calculations for 2D materials: the example of monolayer MoS
- Plasmon-pole approximation for many-body-effects in extrinsic graphene
- Efficient Full-frequency GW Calculations using a Lanczos Method
- Many-body perturbation theory with hybrid density functional theory starting points accelerated by adaptively compressed exchange