Reconstructing lattice QCD spectral functions with stochastic pole expansion and Nevanlinna analytic continuation

Li Huang,Shuang Liang

DOI: https://doi.org/10.1103/physrevd.109.054508

IF: 5.407

2024-03-19

Physical Review D

Abstract:The reconstruction of spectral functions from Euclidean correlation functions is a well-known, yet ill-conditioned inverse problem in the fields of many-body and high-energy physics. In this paper, we present a comprehensive investigation of two recently developed analytic continuation methods, namely stochastic pole expansion and Nevanlinna analytic continuation, for extracting spectral functions from mock lattice QCD data. We examine a range of Euclidean correlation functions generated by representative models, including the Breit-Wigner model, the Gaussian mixture model, the resonance-continuum model, and the bottomonium model. Then, we apply the two methods to reconstruct spectral functions from charmonium correlation functions computed using lattice QCD simulations at a finite temperature. Our findings demonstrate that the stochastic pole expansion method, when combined with the constrained sampling algorithm and the self-adaptive sampling algorithm, successfully recovers the essential features of the spectral functions and exhibits excellent resilience to noise of input data. In contrast, the Nevanlinna analytic continuation method suffers from numerical instability, often resulting in the emergence of spurious peaks and significant oscillations in the high-energy regions of the spectral functions, even with the application of the Hardy basis function optimization algorithm. https://doi.org/10.1103/PhysRevD.109.054508 Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI. Funded by SCOAP 3 . Published by the American Physical Society

astronomy & astrophysics,physics, particles & fields

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the classical but difficult inverse problem of reconstructing spectral functions from Euclidean correlation functions in lattice calculations of quantum chromodynamics (QCD). Specifically, since QCD lattice calculations are usually carried out on a discrete Euclidean space - time grid, it becomes very challenging to directly extract spectral functions from the Euclidean correlation functions obtained by numerical calculations. This is because:

1. **Discretization**: The Euclidean correlation function is evaluated only at a finite number of points, which is caused by the space - time discretization in QCD lattice calculations.
2. **Noise**: QCD lattice simulations rely on random Monte Carlo sampling, so the obtained Euclidean correlation functions are essentially noisy.

These problems make the extraction of spectral functions from Euclidean correlation functions through numerical analytic continuation a ill - conditioned inverse problem, that is, small input deviations or fluctuations will lead to significant uncertainties in the spectral functions. To meet these challenges, researchers have developed a variety of analytic continuation methods, such as the maximum entropy method (MaxEnt), the stochastic analytic continuation method (SAC), etc.

In this paper, the author focuses on two newly proposed analytic continuation methods: stochastic pole expansion (SPX) and Nevanlinna analytic continuation (NAC), and applies them to simulated lattice QCD data to evaluate their performance in extracting spectral functions. By comparing with the traditional maximum entropy method, the author hopes to verify the effectiveness and robustness of these two new methods, especially their performance in dealing with noisy data.