Euclidean Correlation Functions Research Articles

Towards the computation of inclusive decay rates using lattice QCD

We present a non-perturbative computation of inclusive rates of semileptonic decays of heavy mesons from lattice QCD simulations. The calculation is based on the extraction of smeared spectral functions obtained from four-point Euclidean correlation functions computed on configuration ensembles of the JLQCD and ETM collaborations. We compare our results for the inclusive decay rates with analytical predictions from the operator-product expansion, finding a good agreement for the calculation of the inclusive decay rate. This opens the path to the theoretical determination of the magnitude of the Cabibbo-Kobayashi-Maskawa matrix element V_{cb}Vcb to a level of precision competitive with the present experimental uncertainty.

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

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. Published by the American Physical Society 2024

Adiabatic Evolution of Low-Temperature Many-Body Systems

We consider finite-range, many-body fermionic lattice models and we study the evolution of their thermal equilibrium state after introducing a weak and slowly varying time-dependent perturbation. Under suitable assumptions on the external driving, we derive a representation for the average of the evolution of local observables via a convergent expansion in the perturbation, for small enough temperatures. Convergence holds for a range of parameters that is uniform in the size of the system. Under a spectral gap assumption on the unperturbed Hamiltonian, convergence is also uniform in temperature. As an application, our expansion allows us to prove closeness of the time-evolved state to the instantaneous Gibbs state of the perturbed system, in the sense of expectation of local observables, at zero and at small temperatures. As a corollary, we also establish the validity of linear response. Our strategy is based on a rigorous version of the Wick rotation, which allows us to represent the Duhamel expansion for the real-time dynamics in terms of Euclidean correlation functions, for which precise decay estimates are proved using fermionic cluster expansion.

Prospects for γ⋆γ⋆→ππ via lattice QCD

The $\gamma^\star \gamma^\star \to \pi \pi$ scattering amplitude plays a key role in a wide range of phenomena, including understanding the inner structure of scalar resonances as well as constraining the hadronic contributions to the anomalous magnetic moment of the muon. In this work, we explain how the infinite-volume Minkowski amplitude can be constrained from finite-volume Euclidean correlation functions. The relationship between the finite-volume Euclidean correlation functions and the desired amplitude holds up to energies where $3\pi$ states can go on shell, and is exact up to exponentially small corrections that scale like $\mathcal{O}(e^{-m_\pi L})$, where $L$ is the spatial extent of the cubic volume and $m_\pi$ is the pion mass. In order to implement this formalism and remove all power-law finite volume errors, it is necessary to first obtain $\pi \pi \to \pi\pi$, $\pi \gamma^\star \to \pi$, $\gamma^\star \to\pi\pi$, and $\pi\pi\gamma^\star \to\pi\pi$ amplitudes; all of which can be determined via lattice quantum chromodynamic calculations.

Pion spectral properties above the chiral crossover of QCD

Spectral functions encode a wealth of information about the dynamics of any given system, and the determination of their non-perturbative characteristics is a long-standing problem in quantum field theory. Whilst numerical simulations of lattice QCD provide ample data for various Euclidean correlation functions, the inversion required to extract spectral functions is an ill-posed problem. In this work, we pursue previously established constraints imposed by field locality at finite temperature T, namely that spectral functions possess a non-perturbative representation which generalises the well-known Källén-Lehmann spectral form to T > 0. Using this representation, we analyse lattice QCD data of the spatial pseudo-scalar correlator in the temperature range 220–960 MeV, and obtain an analytic expression for the corresponding spectral function, with parameters fixed by the data. From the structure of this spectral function we find evidence for the existence of a distinct pion state above the chiral pseudo-critical temperature Tpc, and contributions from its first excitation, which gradually melt as the temperature increases. As a non-trivial test, we find that the extracted spectral function reproduces the corresponding temporal lattice correlator data for T = 220 MeV.

Rethinking the ill-posedness of the spectral function reconstruction — Why is it fundamentally hard and how Artificial Neural Networks can help

Reconstructing hadron spectral functions through Euclidean correlation functions are of the important missions in lattice QCD calculations. However, in a Källen–Lehmann (KL) spectral representation, the reconstruction is observed to be ill-posed in practice. It is usually ascribed to the fewer observation points compared to the number of points in the spectral function. In this paper, by solving the eigenvalue problem of continuous KL convolution, we show analytically that the ill-posedness of the inversion is fundamental and it exists even for continuous correlation functions. We discussed how to introduce regulators to alleviate the predicament, in which include the Artificial Neural Networks (ANNs) representations recently proposed by the Authors in [28]. The uniqueness of solutions using ANNs representations is manifested analytically and validated numerically. Reconstructed spectral functions using different regularization schemes are also demonstrated, together with their eigen-mode decomposition. We observe that components with large eigenvalues can be reliably reconstructed by all methods, whereas those with low eigenvalues need to be constrained by regulators.

Real-time observables from Euclidean thermal correlation functions

In this work we apply a local quantum field theory approach in order to analyse the connection between real-time observables and Euclidean thermal correlation functions. In particular, using data generated from the functional renormalisation group in the quark-meson model, we demonstrate that in-medium effects can be directly extracted from the spatial momentum dependence of the Euclidean propagators, in contrast to conventional approaches, which rely on the reconstruction from different Matsubara frequencies. As an application, we determine the analytic features that arise from the discrete spectral contribution to the pion correlation function, and calculate the non-perturbative shear viscosity arising from these states.

Euclidean thermal correlation functions in local QFT

In this work we outline the general analytic characteristics satisfied by scalar correlation functions at finite temperature in local quantum field theory. We demonstrate that the locality of the fields in particular imposes significant constraints on the spectral structure of the theory, and that this enables the non-perturbative effects experienced by thermal particle states to be directly calculated from Euclidean correlation functions, avoiding the inverse problem.

Remarks about weighted energy integrals over Minkowski spectral functions from Euclidean lattice data

I make some simple observations about the calculation of weighted averages over energy of Minkowski space spectral densities from weighted averages over time of Euclidean space correlation functions, measured in latice simulations. The correlator of two vector currents is used as an example, where it appears that a determination of a weighted average of the spectral function near the rho pole at the five per cent level is possible from lattice simulations.

Computing light-front wave functions without light-front quantization: A large-momentum effective theory approach

Light-front wave functions play a fundamental role in the light-front quantization approach to QCD and hadron structure. However, a naive implementation of the light-front quantization suffers from various subtleties including the well-known zero-mode problem, the associated rapidity divergences which mixes ultraviolet divergences with infrared physics, as well as breaking of spatial rotational symmetry. We advocate that the light-front quantization should be viewed as an effective theory in which small ${k}^{+}$ modes have been effectively ``integrated out,'' with an infinite number of renormalization constants. Instead of solving light-front quantized field theories directly, we make the large momentum expansion of the equal-time Euclidean correlation functions in instant quantization as an effective way to systematically calculate light-front correlations, including the light-front wave function amplitudes. This large-momentum effective theory accomplishes an effective light-front quantization through lattice QCD calculations. We demonstrate our approach using an example of a pseudoscalar meson wave function.

A PDE Construction of the Euclidean Phi ^4_3 Quantum Field Theory

We present a new construction of the Euclidean Phi ^4 quantum field theory on {mathbb {R}}^3 based on PDE arguments. More precisely, we consider an approximation of the stochastic quantization equation on {mathbb {R}}^3 defined on a periodic lattice of mesh size varepsilon and side length M. We introduce a new renormalized energy method in weighted spaces and prove tightness of the corresponding Gibbs measures as varepsilon rightarrow 0, M rightarrow infty . Every limit point is non-Gaussian and satisfies reflection positivity, translation invariance and stretched exponential integrability. These properties allow to verify the Osterwalder–Schrader axioms for a Euclidean QFT apart from rotation invariance and clustering. Our argument applies to arbitrary positive coupling constant, to multicomponent models with O(N) symmetry and to some long-range variants. Moreover, we establish an integration by parts formula leading to the hierarchy of Dyson–Schwinger equations for the Euclidean correlation functions. To this end, we identify the renormalized cubic term as a distribution on the space of Euclidean fields.

Reconstruction of smeared spectral functions from Euclidean correlation functions

Abstract We propose a method to reconstruct smeared spectral functions from two-point correlation functions measured on the Euclidean lattice. An arbitrary smearing function can be considered as long 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.

Stochastic quantization of a self-interacting nonminimal scalar field in semiclassical gravity

We employ stochastic quantization for a self-interacting nonminimal massive scalar field in curved spacetime. The covariant background field method and local momentum space representation are used to obtain the Euclidean correlation function and evaluate multi-loop quantum corrections through simultaneous expansions in the curvature tensor and its covariant derivatives and in the noise fields. The stochastic correlation function for a quartic self-interaction reproduces the well-known one-loop result by Bunch and Parker and is used to construct the effective potential in curved spacetime in an arbitrary dimension D up to the first order in curvature. Furthermore, we present a sample of numerical simulations for D=3 in the first order in curvature. We consider the model with spontaneous symmetry breaking and obtain fully nonperturbative solutions for the vacuum expectation value of the scalar field and compare them with one- and two-loop solutions.

Scattering amplitudes from finite-volume spectral functions

A novel proposal is outlined to determine scattering amplitudes from finite-volume spectral functions. The method requires extracting smeared spectral functions from finite-volume Euclidean correlation functions, with a particular complex smearing kernel of width $\epsilon$ which implements the standard $i\epsilon$-prescription. In the $L \to \infty$ limit these smeared spectral functions are therefore equivalent to Minkowskian correlators with a specific time ordering to which a modified LSZ reduction formalism can be applied. The approach is presented for general $m \to n$ scattering amplitudes (above arbitrary inelastic thresholds) for a single-species real scalar field, although generalization to arbitrary spins and multiple coupled channels is likely straightforward. Processes mediated by the single insertion of an external current are also considered. Numerical determination of the finite-volume smeared spectral function is discussed briefly and the interplay between the finite volume, Euclidean signature, and time-ordered $i\epsilon$-prescription is illustrated perturbatively in a toy example.

To the sphere and back again: de Sitter infrared correlators at NTLO in 1/N

We analyze the infrared behavior of the two and four-point functions for the massless O(N) model in Lorentzian de Sitter spacetime, using the 1/N expansion. Our approach is based in the study of the Schwinger-Dyson equations on the sphere (Euclidean de Sitter space), using the fact that the infrared behavior in Lorentzian spacetime is determined by the pole structure of the Euclidean correlation functions. We compute the two-point function up to the NTLO in 1/N , and show that in the infrared it behaves as the superposition of two massive free propagators with effective masses of the same order, but not equal to, the dynamical mass mdyn. We compare our results with those obtained using other approaches, and find that they are equivalent but retrieved in a considerably simpler way. We also discuss the infrared behavior of the equal-times four-point functions.

Pion distribution amplitude from Euclidean correlation functions: Exploring universality and higher-twist effects

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.

Quantum Monte Carlo formalism for dynamical pions and nucleons

In most simulations of nonrelativistic nuclear systems, the wave functions found solving the many-body Schr\"odinger equations describe the quantum-mechanical amplitudes of the nucleonic degrees of freedom. In those simulations the pionic contributions are encoded in nuclear potentials and electroweak currents, and they determine the low-momentum behavior. In this work we present an alternative quantum Monte Carlo formalism in which both relativistic pions and nonrelativistic nucleons are explicitly included in the quantum-mechanical states of the system. We report the renormalization of the nucleon mass as a function of the momentum cutoff, an Euclidean time density correlation function that deals with the short-time nucleon diffusion, and the pion cloud density and momentum distributions. In the two-nucleon sector we show that the interaction of two static nucleons at large distances reduces to the one-pion exchange potential, and we fit the low-energy constants of the contact interactions to reproduce the binding energy of the deuteron and two neutrons in finite volumes. We show that the method can be readily applied to light-nuclei.

Stochastic reconstructions of spectral functions: Application to lattice QCD

We present a detailed study of the applications of two stochastic approaches, stochastic optimization method (SOM) and stochastic analytical inference (SAI), to extract spectral functions from Euclidean correlation functions. SOM has the advantage that it does not require prior information. On the other hand, SAI is a more generalized method based on Bayesian inference. Under mean field approximation SAI reduces to the often-used maximum entropy method (MEM), and for a specific choice of the prior SAI becomes equivalent to SOM. To test the applicability of these two stochastic methods to lattice QCD, firstly, we apply these methods to various reasonably chosen model correlation functions, and present detailed comparisons of the reconstructed spectral functions obtained from SOM, SAI and MEM. Next, we present similar studies for charmonia correlation functions obtained from lattice QCD computations using clover-improved Wilson fermions on large, fine, isotropic lattices at $0.75$ and $1.5T_c$, $T_c$ being the deconfinement transition temperature of a pure gluon plasma. We find that SAI and SOM give consistent results to MEM at these two temperatures.

Pion distribution amplitude from Euclidean correlation functions

Following the proposal in (Braun and Müller. Eur Phys J C55:349, 2008), we study the feasibility to calculate the pion distribution amplitude (DA) from suitably chosen Euclidean correlation functions at large momentum. In our lattice study we employ the novel momentum smearing technique (Bali et al. Phys Rev D93:094515, 2016; Bali et al. Phys Lett B774:91, 2017). This approach is complementary to the calculations of the lowest moments of the DA using the Wilson operator product expansion and avoids mixing with lower dimensional local operators on the lattice. The theoretical status of this method is similar to that of quasi-distributions (Ji. Phys Rev Lett 110:262002, 2013) that have recently been used in (Zhang et al. Phys Rev D95:094514, 2017) to estimate the twist two pion DA. The similarities and differences between these two techniques are highlighted.

An estimate for the thermal photon rate from lattice QCD

We estimate the production rate of photons by the quark-gluon plasma in lattice QCD. We propose a new correlation function which provides better control over the systematic uncertainty in estimating the photon production rate at photon momenta in the range πT/2 to 2πT. The relevant Euclidean vector current correlation functions are computed with Nf = 2 Wilson clover fermions in the chirally-symmetric phase. In order to estimate the photon rate, an ill-posed problem for the vector-channel spectral function must be regularized. We use both a direct model for the spectral function and a modelindependent estimate from the Backus-Gilbert method to give an estimate for the photon rate.