No Arabic abstract
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 discuss, how to study $I = 0$ quarkonium resonances decaying into pairs of heavy-light mesons using static potentials from lattice QCD. These static potentials can be obtained from a set of correlation functions containing both static and light quarks. As a proof of concept we focus on bottomonium with relative orbital angular momentum $L = 0$ of the $bar{b} b$ pair corresponding to $J^{P C} = 0^{- +}$ and $J^{P C} = 1^{- -}$. We use static potentials from an existing lattice QCD string breaking study and compute phase shifts and $mbox{T}$ matrix poles for the lightest heavy-light meson-meson decay channel. We discuss our results in the context of corresponding experimental results, in particular for $Upsilon (10860)$ and $Upsilon (11020)$.
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.
Single state saturation of the temporal correlation function is a key condition to extract physical observables such as energies and matrix elements of hadrons from lattice QCD simulations. A method commonly employed to check the saturation is to seek for a plateau of the observables for large Euclidean time. Identifying the plateau in the cases having nearby states, however, is non-trivial and one may even be misled by a fake plateau. Such a situation takes place typically for the system with two or more baryons. In this study, we demonstrate explicitly the danger from a possible fake plateau in the temporal correlation functions mainly for two baryons ($XiXi$ and $NN$), and three and four baryons ($^3{rm He}$ and $^4{rm He})$ as well, employing (2+1)-flavor lattice QCD at $m_{pi}=0.51$ GeV on four lattice volumes with $L=$ 2.9, 3.6, 4.3 and 5.8 fm. Caution is given for drawing conclusion on the bound $NN$, $3N$ and $4N$ systems only based on the temporal correlation functions.
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.
In this work we present the first non-perturbative determination of the hadronic susceptibilities that constrain the form factors entering the semileptonic $B to D^{(*)} ell u_ell $ transitions due to unitarity and analyticity. The susceptibilities are obtained by evaluating moments of suitable two-point correlation functions obtained on the lattice. Making use of the gauge ensembles produced by the Extended Twisted Mass Collaboration with $N_f = 2+1+1$ dynamical quarks at three values of the lattice spacing ($a simeq 0.062, 0.082, 0.089$ fm) and with pion masses in the range $simeq 210 - 450$ MeV, we evaluate the longitudinal and transverse susceptibilities of the vector and axial-vector polarization functions at the physical pion point and in the continuum and infinite volume limits. The ETMC ratio method is adopted to reach the physical $b$-quark mass $m_b^{phys}$. At zero momentum transfer for the $b to c$ transition we get $chi_{0^+}(m_b^{phys}) = 7.58,(59) cdot 10^{-3}$, $chi_{1^-}(m_b^{phys}) = 6.72,(41) cdot 10^{-4}$ GeV$^{-2}$, $chi_{0^-}(m_b^{phys}) = 2.58,(17) cdot 10^{-2}$ and $chi_{1^+}(m_b^{phys}) = 4.69,(30) cdot 10^{-4}$ GeV$^{-2}$ for the scalar, vector, pseudoscalar and axial susceptibilities, respectively. In the case of the vector and pseudoscalar channels the one-particle contributions due to $B_c^*$- and $B_c$-mesons are evaluated and subtracted to improve the bounds, obtaining: $chi_{1^-}^{sub}(m_b^{phys}) = 5.84,(44) cdot 10^{-4}$ GeV$^{-2}$ and $chi_{0^-}^{sub}(m_b^{phys}) = 2.19,(19) cdot 10^{-2}$.