We propose a method to reconstruct smeared spectral functions from two-point correlation functions measured on the Euclidean lattice. Arbitrary smearing function can be considered as far as it is smooth enough to allow an approximation using Chebyshev polynomials. We test the method with numerical lattice data of Charmonium correlators. The method provides a framework to compare lattice calculation with experimental data including excited state contributions without assuming quark-hadron duality.
A relation is presented between single-hadron long-range matrix elements defined in a finite Euclidean spacetime, and the corresponding infinite-volume Minkowski amplitudes. This relation is valid in the kinematic region where any number of two-hadron states can simultaneously go on shell, so that the effects of strongly-coupled intermediate channels are included. These channels can consist of non-identical particles with arbitrary intrinsic spins. The result accommodates general Lorentz structures as well as non-zero momentum transfer for the two external currents inserted between the single-hadron states. The formalism, therefore, generalizes the work by Christ et al.~[Phys.Rev. D91 114510 (2015)], and extends the reach of lattice quantum chromodynamics (QCD) to a wide class of new observables beyond meson mixing and rare decays. Applications include Compton scattering of the pion ($pi gamma^star to [pi pi, K overline K] to pi gamma^star$), kaon ($K gamma^star to [pi K, eta K] to K gamma^star$) and nucleon ($N gamma^star to N pi to N gamma^star$), as well as double-$beta$ decays, and radiative corrections to the single-$beta$ decay, of QCD-stable hadrons. The framework presented will further facilitate generalization of the result to studies of nuclear amplitudes involving two currents from lattice QCD.
Building upon our recent study arXiv:1709.04325, we investigate the feasibility of calculating the pion distribution amplitude from suitably chosen Euclidean correlation functions at large momentum. We demonstrate in this work the advantage of analyzing several correlation functions simultaneously and extracting the pion distribution amplitude from a global fit. This approach also allows us to study higher-twist corrections, which are a major source of systematic error. Our result for the higher-twist parameter $delta^pi_2$ is in good agreement with estimates from QCD sum rules. Another novel element is the use of all-to-all propagators, calculated using stochastic estimators, which enables an additional volume average of the correlation functions, thereby reducing statistical errors.
We propose the sparse modeling method to estimate the spectral function from the smeared correlation functions. We give a description of how to obtain the shear viscosity from the correlation function of the renormalized energy-momentum tensor (EMT) measured by the gradient flow method ($C(t,tau)$) for the quenched QCD at finite temperature. The measurement of the renormalized EMT in the gradient flow method reduces a statistical uncertainty thanks to its property of the smearing. However, the smearing breaks the sum rule of the spectral function and the over-smeared data in the correlation function may have to be eliminated from the analyzing process of physical observables. In this work, we demonstrate that the sparse modeling analysis in the intermediate-representation basis (IR basis), which connects between the Matsubara frequency data and real frequency data. It works well even using very limited data of $C(t,tau)$ only in the fiducial window of the gradient flow. We utilize the ADMM algorithm which is useful to solve the LASSO problem under some constraints. We show that the obtained spectral function reproduces the input smeared correlation function at finite flow-time. Several systematic and statistical errors and the flow-time dependence are also discussed.
We present a fast and simple algorithm that allows the extraction of multiple exponential signals from the temporal dependence of correlation functions evaluated on the lattice including the statistical fluctuations of each signal and treating properly backward signals. The basic steps of the method are the inversion of appropriate matrices and the determination of the roots of an appropriate polynomial, constructed using discretized derivatives of the correlation function. The method is tested strictly using fake data generated assuming a fixed number of exponential signals included in the correlation function with a controlled numerical precision and within given statistical fluctuations. All the exponential signals together with their statistical uncertainties are determined exactly by the algorithm. The only limiting factor is the numerical rounding off. In the case of correlation functions evaluated by large-scale QCD simulations on the lattice various sources of noise, other than the numerical rounding, can affect the correlation function and they represent the crucial factor limiting the number of exponential signals, related to the hadronic spectral decomposition of the correlation function, that can be properly extracted. The algorithm can be applied to a large variety of correlation functions typically encountered in QCD or QCD+QED simulations on the lattice, including the case of exponential signals corresponding to poles with arbitrary multiplicity and/or the case of oscillating signals. The method is able to to detect the specific structure of the multiple exponential signals without any a priori assumption and it determines accurately the ground-state signal without the need that the lattice temporal extension is large enough to allow the ground-state signal to be isolated.
We present a numerical pilot study of the meson correlation functions in the epsilon-regime of chiral perturbation theory. Based on simulations with overlap fermions we measured the axial and pseudo-scalar correlation functions, and we discuss the implications for the leading low energy constants in the chiral Lagrangian.