Extending the average spectrum method: Grid points sampling and density averaging

Abstract

Analytic continuation of imaginary time or frequency data to the real axis is a crucial step in extracting dynamical properties from quantum Monte Carlo simulations. The average spectrum method provides an elegant solution by integrating over all non-negative spectra weighted by how well they fit the data. In a recent paper, we found that discretizing the functional integral as in Feynman's path-integrals, does not have a well-defined continuum limit. Instead, the limit depends on the discretization grid whose choice may strongly bias the results. In this paper, we demonstrate that sampling the grid points, instead of keeping them fixed, also changes the functional integral limit and rather helps to overcome the bias considerably. We provide an efficient algorithm for doing the sampling and show how the density of the grid points acts now as a default model with a significantly reduced biasing effect. The remaining bias depends mainly on the width of the grid density, so we go one step further and average also over densities of different widths. For a certain class of densities, including Gaussian and exponential ones, this width averaging can be done analytically, eliminating the need to specify this parameter without introducing any computational overhead.

arXiv comment10 pages, 10 figures

References in corpus (16)

• Machine learning for many-body physics: The case of the Anderson impurity model
• A reliable Pade analytical continuation method based on a high accuracy symbolic computation algorithm
• Algorithms for optimized maximum entropy and diagnostic tools for analytic continuation
• Constrained sampling method for analytic continuation
• Implementation of the Maximum Entropy Method for Analytic Continuation
• Artificial Neural Network Approach to the Analytic Continuation Problem
• Analytic Continuation of Quantum Monte Carlo Data by Stochastic Analytical Inference
• Using the average spectrum method to extract dynamics from quantum Monte Carlo simulations
• Identifying the maximum entropy method as a special limit of stochastic analytic continuation
• Analytic continuation via domain-knowledge free machine learning
• Analytic continuation by averaging Padé approximants
• Analytical continuation of imaginary axis data for optical conductivity
• Analytical continuation of imaginary axis data using maximum entropy
• Maximum Entropy Analytic Continuation for Spectral Functions with Non-Positive Spectral Weight
• The Average Spectrum Method for Analytic Continuation: Efficient Blocked Modes Sampling and Dependence on Discretization Grid
• One particle spectral function and analytic continuation for many-body implementation in the EMTO method

Cited by in corpus (13)

• Progress on stochastic analytic continuation of quantum Monte Carlo data
• ACFlow: An open source toolkit for analytical continuation of quantum Monte Carlo data
• Analytical solution for time-integrals in diagrammatic expansions: application to real-frequency diagrammatic Monte Carlo
• Learned Optimizers for Analytic Continuation
• Noise Enhanced Neural Networks for Analytic Continuation
• Stochastic pole expansion method
• Robust Extraction of Thermal Observables from State Sampling and Real-Time Dynamics on Quantum Computers
• Reconstructing lattice QCD spectral functions with stochastic pole expansion and Nevanlinna analytic continuation
• Connecting Tikhonov regularization to the maximum entropy method for the analytic continuation of quantum Monte Carlo data
• Magnon, doublon and quarton excitations in 2D S=1/2 trimerized Heisenberg models
• Single-Particle Dispersion and Density of States of the Half-Filled 2D Hubbard Model
• Barycentric rational function approximation made simple: A fast analytic continuation method for Matsubara Green's functions
• Cross Validation in Stochastic Analytic Continuation