The Average Spectrum Method for Analytic Continuation: Efficient Blocked Modes Sampling and Dependence on Discretization Grid

Abstract

The average spectrum method is a promising approach for the analytic continuation of imaginary time or frequency data to the real axis. It determines the analytic continuation of noisy data from a functional average over all admissible spectral functions, weighted by how well they fit the data. Its main advantage is the apparent lack of adjustable parameters and smoothness constraints, using instead the information on the statistical noise in the data. Its main disadvantage is the enormous computational cost of performing the functional integral. Here we introduce an efficient implementation, based on the singular value decomposition of the integral kernel, eliminating this problem. It allows us to analyze the behavior of the average spectrum method in detail. We find that the discretization of the real-frequency grid, on which the spectral function is represented, biases the results. The distribution of the grid points plays the role of a default model while the number of grid points acts as a regularization parameter. We give a quantitative explanation for this behavior, point out the crucial role of the default model and provide a practical method for choosing it, making the average spectrum method a reliable and efficient technique for analytic continuation.

12 pages, 10 figures

References in corpus (10)

• Simulation of truncated normal variables
• Orthogonal Polynomial Representation of Imaginary-Time Green's Functions
• A novel Bayesian approach to spectral function reconstruction
• Constrained sampling method for analytic continuation
• 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
• Analytical continuation of imaginary axis data for optical conductivity
• Maximum Entropy Analytic Continuation for Spectral Functions with Non-Positive Spectral Weight
• Analytic Continuation of Quantum Monte Carlo Data: Optimal Stochastic Regularization Approach

Cited by in corpus (18)

• Progress on stochastic analytic continuation of quantum Monte Carlo data
• Extending the average spectrum method: Grid points sampling and density averaging
• ACFlow: An open source toolkit for analytical continuation of quantum Monte Carlo data
• Learned Optimizers for Analytic Continuation
• Coulomb Engineering of two-dimensional Mott materials
• 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
• Training biases in machine learning for the analytic continuation of quantum many-body Green's functions
• Magnon, doublon and quarton excitations in 2D S=1/2 trimerized Heisenberg models
• Barycentric rational function approximation made simple: A fast analytic continuation method for Matsubara Green's functions
• Supersolid Phase in the Diluted Holstein Model
• Single-Particle Dispersion and Density of States of the Half-Filled 2D Hubbard Model
• Cross Validation in Stochastic Analytic Continuation
• Parameter-free analytic continuation for quantum many-body calculations
• A Parameter-Free Differential Evolution Algorithm for the Analytic Continuation of Imaginary Time Correlation Functions