Efficient temperature-dependent Green's function methods for realistic systems: using cubic spline interpolation to approximate Matsubara Green's functions

Abstract

The popular, stable, robust and computationally inexpensive cubic spline interpolation algorithm is adopted and used for finite temperature Green's function calculations of realistic systems. We demonstrate that with appropriate modifications the temperature dependence can be preserved while the Green's function grid size can be reduced by about two orders of magnitude by replacing the standard Matsubara frequency grid with a sparser grid and a set of interpolation coefficients. We benchmarked the accuracy of our algorithm as a function of a single parameter sensitive to the shape of the Green's function. Through numerous examples, we confirmed that our algorithm can be utilized in a systematically improvable, controlled, and black-box manner and highly accurate one- and two-body energies and one-particle density matrices can be obtained using only around 5% of the original grid points. Additionally, we established that to improve accuracy by an order of magnitude, the number of grid points needs to be doubled, whereas for the Matsubara frequency grid an order of magnitude more grid points must be used. This suggests that realistic calculations with large basis sets that were previously out of reach because they required enormous grid sizes may now become feasible.

References in corpus (5)

• Continuous-time Monte Carlo methods for quantum impurity models
• Hybridization expansion impurity solver: General formulation and application to Kondo lattice and two-orbital models
• Truncated Configuration Interaction expansions as solvers for correlated quantum impurity models and dynamical mean field theory
• Systematically improvable multi-scale solver for correlated electron systems
• The correlation potential in density functional theory at the GW-level: spherical atoms

Cited by in corpus (33)

• Towards the solution of the many-electron problem in real materials: equation of state of the hydrogen chain with state-of-the-art many-body methods
• Sparse sampling approach to efficient ab initio calculations at finite temperature
• Finite temperature quantum embedding theories for correlated systems
• Self-energy embedding theory (SEET) for periodic systems
• Generalized self-energy embedding theory
• Rigorous ab initio quantum embedding for quantum chemistry using Green's function theory: screened interaction, non-local self-energy relaxation, orbital basis, and chemical accuracy
• Exploring connections between statistical mechanics and Green's functions for realistic systems. Temperature dependent electronic entropy and internal energy from a self-consistent second-order Green's function
• Chebyshev polynomial representation of imaginary time response functions
• Comparing self-consistent GW and vertex corrected G0W0 (G0W0Γ) accuracy for molecular ionization potentials
• Finite Temperature Auxiliary Field Quantum Monte Carlo in the Canonical Ensemble
• Legendre-spectral Dyson equation solver with super-exponential convergence
• Testing self-energy embedding theory in combination with GW
• Effect of propagator renormalization on the band gap of insulating solids
• A wave function perspective and efficient truncation of renormalised second-order perturbation theory
• Spectral properties from Matsubara Green's function approach - application to molecules
• Finite-temperature coupled cluster: Efficient implementation and application to prototypical systems
• Spin-unrestricted self-energy embedding theory
• A regularized second-order correlation method from Green's function theory
• Unveiling the Finite Temperature Physics of Hydrogen Chains via Auxiliary Field Quantum Monte Carlo
• Minimal Pole Representation and Controlled Analytic Continuation of Matsubara Response Functions
• Minimal pole representation and analytic continuation of matrix-valued correlation functions
• Challenges with relativistic GW calculations in solids and molecules
• Nonconvergence of the Feynman-Dyson diagrammatic perturbation expansion of propagators
• Tensor hypercontraction for fully self-consistent imaginary-time GF2 and GWSOX methods: theory, implementation, and role of the Green's function second-order exchange for intermolecular interactions
• Combining Density Functional Theory and Green's Function Theory: Range-Separated, Non-local, Dynamic, and Orbital-Dependent Hybrid Functional
• Minimal pole representation for spectral functions
• Normal ordered exponential approach to thermal properties and time-correlation functions: General theory and simple examples
• Causal optimization method for imaginary-time Green's functions in interacting electron systems
• Heating and cooling in self-consistent many-body simulations
• Compact representation and long-time extrapolation of real-time data for quantum systems using the ESPRIT algorithm
• Electronic specific heat capacities and entropies from density matrix quantum Monte Carlo using Gaussian process regression to find gradients of noisy data
• Piecewise Interaction Picture Density Matrix Quantum Monte Carlo
• Modified discrete Laguerre polynomials for efficient computation of exponentially bounded Matsubara sums