Denoising of Imaginary Time Response Functions with Hankel projections

Abstract

Imaginary-time response functions of finite-temperature quantum systems are often obtained with methods that exhibit stochastic or systematic errors. Reducing these errors comes at a large computational cost -- in quantum Monte Carlo simulations, the reduction of noise by a factor of two incurs a simulation cost of a factor of four. In this paper, we relate certain imaginary-time response functions to an inner product on the space of linear operators on Fock space. We then show that data with noise typically does not respect the positive definiteness of its associated Gramian. The Gramian has the structure of a Hankel matrix. As a method for denoising noisy data, we introduce an alternating projection algorithm that finds the closest positive definite Hankel matrix consistent with noisy data. We test our methodology at the example of fermion Green's functions for continuous-time quantum Monte Carlo data and show remarkable improvements of the error, reducing noise by a factor of up to 20 in practical examples. We argue that Hankel projections should be used whenever finite-temperature imaginary-time data of response functions with errors is analyzed, be it in the context of quantum Monte Carlo, quantum computing, or in approximate semianalytic methodologies.

5 figures

References in corpus (32)

• Electronic Structure Calculations with Dynamical Mean-Field Theory: A Spectral Density Functional Approach
• Continuous-time Monte Carlo methods for quantum impurity models
• Quantum Cluster Theories
• A continuous-time solver for quantum impurity models
• The ALPS project release 2.0: Open source software for strongly correlated systems
• Continuous Time Quantum Monte Carlo method for fermions
• Nonlocal Dynamical Correlations of Strongly Interacting Electron Systems
• Hybridization expansion impurity solver: General formulation and application to Kondo lattice and two-orbital models
• Full nonperturbative QCD simulations with 2+1 flavors of improved staggered quarks
• Real-time path integral approach to nonequilibrium many-body quantum system
• Polaron Problem by Diagrammatic Quantum Monte Carlo
• Diagrammatic Monte Carlo simulation of non-equilibrium systems
• TRIQS/CTHYB: A Continuous-Time Quantum Monte Carlo Hybridization Expansion Solver for Quantum Impurity Problems
• Continuous-time auxiliary field Monte Carlo for quantum impurity models
• Taming the dynamical sign problem in real-time evolution of quantum many-body problems
• Orthogonal Polynomial Representation of Imaginary-Time Green's Functions
• Deviations from Fermi-liquid behavior above $T_c$ in 2D short coherence length superconductors
• A Quantum Monte Carlo algorithm for non-local corrections to the Dynamical Mean-Field Approximation
• w2dynamics: Local one- and two-particle quantities from dynamical mean field theory
• Bold Diagrammatic Monte Carlo: When Sign Problem is Welcome
• Nevanlinna Analytical Continuation
• Recent developments in Quantum Monte-Carlo simulations with applications for cold gases
• Determinant Diagrammatic Monte Carlo in the Thermodynamic Limit
• Learning Feynman Diagrams with Tensor Trains
• Analytical Continuation of Matrix-Valued Functions: Carathéodory Formalism
• Diagrammatic expansion for positive spectral functions beyond GW: Application to vertex corrections in the electron gas
• Continuous-time hybridization expansion quantum impurity solver for multi-orbital systems with complex hybridizations
• Robust analytic continuation of Green's functions via projection, pole estimation, and semidefinite relaxation
• Quantum Monte Carlo in the steady-state
• A Tensor Train Continuous Time Solver for Quantum Impurity Models
• Magnetic and charge susceptibilities in the half-filled triangular lattice Hubbard model
• Minimal Pole Representation and Controlled Analytic Continuation of Matsubara Response Functions

Cited by in corpus (3)

• Adaptive variational quantum computing approaches for Green's functions and nonlinear susceptibilities
• Compact representation and long-time extrapolation of real-time data for quantum systems using the ESPRIT algorithm
• Autoencoder-based analytic continuation method for strongly correlated quantum systems