PyLIT: Reformulation and implementation of the analytic continuation problem using kernel representation methods

Abstract

Path integral Monte Carlo (PIMC) simulations are a cornerstone for studying quantum many-body systems. The analytic continuation (AC) needed to estimate dynamic quantities from these simulations is an inverse Laplace transform, which is ill-conditioned. If this inversion were surmounted, then dynamical observables (e.g. dynamic structure factor (DSF) $S(q,\omega )$) could be extracted from the imaginary-time correlation functions estimates. Although of important, the AC problem remains challenging due to its ill-posedness. To address this challenge, we express the DSF as a linear combination of kernel functions with known Laplace transforms that have been tailored to satisfy its physical constraints. We use least-squares optimization regularized with a Bayesian prior to determine the coefficients of this linear combination. We explore various regularization term, such as the commonly used entropic regularizer, as well as the Wasserstein distance and $L^2$-distance as well as techniques for setting the regularization weight. A key outcome is the open-source package PyLIT (\textbf{Py}thon \textbf{L}aplace \textbf{I}nverse \textbf{T}ransform), which leverages Numba and unifies the presented formulations. PyLIT's core functionality is kernel construction and optimization. In our applications, we find PyLIT's DSF estimates share qualitative features with other more established methods. We identify three key findings. Firstly, independent of the regularization choice, utilizing non-uniform grid point distributions reduced the number of unknowns and thus reduced our space of possible solutions. Secondly, the Wasserstein distance, a previously unexplored regularizer, performs as good as the entropic regularizer while benefiting from its linear gradient. Thirdly, future work can meaningfully combine regularized and stochastic optimization. (text cut for char. limit)

29 pages, 11 figures, submitted to Computer Physics Communications

References in corpus (48)

• Rényi Divergence and Kullback-Leibler Divergence
• Computational complexity and fundamental limitations to fermionic quantum Monte Carlo simulations
• Worm Algorithm and Diagrammatic Monte Carlo: A New Approach to Continuous-Space Path Integral Monte Carlo Simulations
• Worm Algorithm for Continuous-space Path Integral Monte Carlo Simulations
• The Uniform Electron Gas at Warm Dense Matter Conditions
• A novel Bayesian approach to spectral function reconstruction
• \textit{Ab Initio} Path Integral Monte Carlo Results for the Dynamic Structure Factor of Correlated Electrons: From the Electron Liquid to Warm Dense Matter
• Excitation spectrum of a supersolid
• The Fermion Sign Problem in Path Integral Monte Carlo Simulations: Quantum Dots, Ultracold Atoms, and Warm Dense Matter
• Algorithms for optimized maximum entropy and diagnostic tools for analytic continuation
• Ab-initio low-energy dynamics of superfluid and solid 4He
• Constrained sampling method for analytic continuation
• Sparse modeling approach to analytical continuation of imaginary-time quantum Monte Carlo data
• On the extraction of spectral densities from lattice correlators
• The Static Local Field Correction of the Warm Dense Electron Gas: An ab Initio Path Integral Monte Carlo Study and Machine Learning Representation
• Electronic Density Response of Warm Dense Matter
• Artificial Neural Network Approach to the Analytic Continuation Problem
• Numerical analytic continuation: Answers to well-posed questions
• Effective Static Approximation: A Fast and Reliable Tool for Warm Dense Matter Theory
• 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
• Analytic Continuation of Quantum Monte Carlo Data by Stochastic Analytical Inference
• Accurate Temperature Diagnostics for Matter under Extreme Conditions
• Numerical analytic continuation of Euclidean data
• Uniform electron gas at finite temperatures
• Analytic continuation via domain-knowledge free machine learning
• Analytical representation of the Local Field Correction of the Uniform Electron Gas within the Effective Static Approximation
• The Strongly Coupled Electron Liquid: ab initio Path Integral Monte Carlo Simulations and Dielectric Theories
• Emergence of an excitonic collective mode in the dilute electron gas
• Density Response of the Warm Dense Electron Gas beyond Linear Response Theory: Excitation of Harmonics
• Nonlinear density response from imaginary-time correlation functions: Ab initio path integral Monte Carlo simulations of the warm dense electron gas
• Finite-size effects in the reconstruction of dynamic properties from ab initio path integral Monte-Carlo simulations
• Ab initio path integral Monte Carlo approach to the momentum distribution of the uniform electron gas at finite temperature without fixed nodes
• Dynamic structure factor of liquid 4He across the normal-superfluid transition
• Sparse Modeling in Quantum Many-Body Problems
• Collective and single-particle excitations in 2D dipolar Bose gases
• X-ray Thomson scattering spectra from DFT-MD simulations based on a modified Chihara formula
• Correlation effects and collective excitations in bosonic bilayers: role of quantum statistics, superfluidity and dimerization transition
• Dynamic structure factor of superfluid He-4 from Quantum Monte Carlo: Maximum Entropy revisited
• Extraction of the frequency moments of spectral densities from imaginary-time correlation function data
• Spectral properties from Matsubara Green's function approach - application to molecules
• Model-free Rayleigh weight from x-ray Thomson scattering measurements
• Analysis of dynamical effects in the uniform electron liquids with the self-consistent method of moments complemented by the Shannon information entropy and the path-integral Monte-Carlo simulations
• ACFlow: An open source toolkit for analytical continuation of quantum Monte Carlo data
• A comparison between methods of analytical continuation for bosonic functions
• Estimates of the dynamic structure factor for the finite temperature electron liquid via analytic continuation of path integral Monte Carlo data
• Applying the Liouville-Lanczos Method of Time-Dependent Density-Functional Theory to Warm Dense Matter
• Conservative dielectric functions and electrical conductivities from the multicomponent Bhatnagar-Gross-Krook equation