The noiseless limit and improved-prior limit of the maximum entropy method and their implications for the analytic continuation problem

Abstract

Quantum Monte Carlo (QMC) methods are uniquely capable of providing exact simulations of quantum many-body systems. Unfortunately, the applications of a QMC simulation are limited because extracting dynamic properties requires solving the analytic continuation (AC) problem. Across the many fields that use QMC methods, there is no universally accepted analytic continuation algorithm for extracting dynamic properties, but many publications compare to the maximum entropy method. We investigate when entropy maximization is an acceptable approach. We show that stochastic sampling algorithms reduce to entropy maximization when the Bayesian prior is near to the true solution. We investigate when is Bryan's controversial optimization algorithm [Bryan, Eur. Biophys. J. 18, 165-174 (1990)] for entropy maximization (sometimes known as the maximum entropy method) appropriate to use. We show that Bryan's algorithm is appropriate when the noise is near zero or when the Bayesian prior is near to the true solution. We also investigate the mean squared error, finding a better scaling when the Bayesian prior is near the true solution than when the noise is near zero. We point to examples of improved data-driven Bayesian priors that have already leveraged this advantage. We support these results by solving the double Gaussian problem using both Bryan's algorithm and the newly formulated dual approach to entropy maximization [Chuna et al., J. Phys. A: Math. Theor. 58, 335203 (2025)].

6 pages, 3 figures, 26 equations

References in corpus (23)

• Worm Algorithm and Diagrammatic Monte Carlo: A New Approach to Continuous-Space Path Integral Monte Carlo Simulations
• Conformal field theories at non-zero temperature: operator product expansions, Monte Carlo, and holography
• Electronic Density Response of Warm Dense Matter
• Numerical analytic continuation: Answers to well-posed questions
• Progress on stochastic analytic continuation of quantum Monte Carlo data
• Ab Initio Path Integral Monte Carlo Approach to the Static and Dynamic Density Response of the Uniform Electron Gas
• Analytical representation of the Local Field Correction of the Uniform Electron Gas within the Effective Static Approximation
• Nonlinear density response from imaginary-time correlation functions: Ab initio path integral Monte Carlo simulations of the warm dense electron gas
• Collective and single-particle excitations in 2D dipolar Bose gases
• Sparse Modeling in Quantum Many-Body Problems
• Analytical continuation of imaginary axis data using maximum entropy
• X-ray Thomson scattering spectra from DFT-MD simulations based on a modified Chihara formula
• Physics-Driven Learning for Inverse Problems in Quantum Chromodynamics
• Rethinking the ill-posedness of the spectral function reconstruction -- why is it fundamentally hard and how Artificial Neural Networks can help
• ACFlow: An open source toolkit for analytical continuation of quantum Monte Carlo data
• Fast and Efficient Stochastic Optimization for Analytic Continuation
• Stochastic pole expansion method
• Applying the Liouville-Lanczos Method of Time-Dependent Density-Functional Theory to Warm Dense Matter
• Estimates of the dynamic structure factor for the finite temperature electron liquid via analytic continuation of path integral Monte Carlo data
• Dual formulation of the maximum entropy method applied to analytic continuation of quantum Monte Carlo data
• Enhancing the Efficiency of Time-Dependent Density Functional Theory Calculations of Dynamic Response Properties
• High-precision Quantum Monte-Carlo study of charge transport in a lattice model of molecular organic semiconductors
• A Momentum-Resolved X-ray Thomson Scattering Benchmark of Electronic-Response Models in Warm Dense Aluminium