Stochastic pole expansion method

Abstract

In this paper, we propose a new analytic continuation method to extract real frequency spectral functions from imaginary frequency Green's functions of quantum many-body systems. This method is based on the pole representation of Matsubara Green's function and a stochastic sampling procedure is utilized to optimize the amplitudes and locations of poles. In order to capture narrow peaks and sharp band edges in the spectral functions, a constrained sampling algorithm and a self-adaptive sampling algorithm are developed. To demonstrate the usefulness and performance of the new method, we at first apply it to study the spectral functions of representative fermionic and bosonic correlators. Then we employ this method to tackle the analytic continuation problems of matrix-valued Green's functions. The synthetic Green's functions, as well as realistic correlation functions from finite temperature quantum many-body calculations, are used as input. The benchmark results demonstrate that this method is capable of reproducing most of the key characteristics in the spectral functions. The sharp, smooth, and multi-peak features in both low-frequency and high-frequency regions of spectral functions could be accurately resolved, which overcomes one of the main limitations of the traditional maximum entropy method. More importantly, it exhibits excellent robustness with respect to noisy and incomplete input data. The causality of spectral function is always satisfied even in the presence of sizable noises. As a byproduct, this method could derive a fitting formula for the Matsubara data, which provides a compact approximation to the many-body Green's functions. Hence, we expect that this new method could become a pivotal workhorse for numerically analytic continuation and be broadly useful in many applications.

26 pages, 20 figures

References in corpus (19)

• Continuous-time Monte Carlo methods for quantum impurity models
• Electrodynamics of Correlated Electron Materials
• Nearly deconfined spinon excitations in the square-lattice spin-1/2 Heisenberg antiferromagnet
• Beyond extended dynamical mean-field theory: Dual boson approach to the two-dimensional extended Hubbard model
• Numerical analytic continuation: Answers to well-posed questions
• Using the average spectrum method to extract dynamics from quantum Monte Carlo simulations
• On the dangers of partial diagrammatic summations: Benchmarks for the two-dimensional Hubbard model in the weak-coupling regime
• $i$QIST: An open source continuous-time quantum Monte Carlo impurity solver toolkit
• Analytical Continuation of Matrix-Valued Functions: Carathéodory Formalism
• Analytical continuation of imaginary axis data using maximum entropy
• Pairing Glue in the Two Dimensional Hubbard Model
• $i$QIST v0.7: An open source continuous-time quantum Monte Carlo impurity solver toolkit
• Statistical and computational intelligence approach to analytic continuation in Quantum Monte Carlo
• The Average Spectrum Method for Analytic Continuation: Efficient Blocked Modes Sampling and Dependence on Discretization Grid
• Bosonic Nevanlinna Analytic Continuation
• A comparison between methods of analytical continuation for bosonic functions
• Fast and Efficient Stochastic Optimization for Analytic Continuation
• Maximum entropy analytic continuation for frequency-dependent transport coefficients with non-positive spectral weight
• Padé approximants and analytic continuation of Euclidean Phi-derivable approximations

Cited by in corpus (7)

• Minimal pole representation and analytic continuation of matrix-valued correlation functions
• $\text{Ab initio}$ dynamical mean-field theory with natural orbitals renormalization group impurity solver: Formalism and applications
• Reconstructing lattice QCD spectral functions with stochastic pole expansion and Nevanlinna analytic continuation
• Barycentric rational function approximation made simple: A fast analytic continuation method for Matsubara Green's functions
• Low-temperature transport in high-conductivity correlated metals: a density-functional plus dynamical mean-field study of cubic perovskites
• Autoencoder-based analytic continuation method for strongly correlated quantum systems
• The noiseless limit and improved-prior limit of the maximum entropy method and their implications for the analytic continuation problem