Ruling out some predictions of deeply-bound light-heavy tetraquarks using lattice QCD

We discuss our lattice QCD calculations of a number of tetraquark channels with at least one heavy quark where some phenomenological models, already fully constrained by fits to the ordinary meson and baryon spectrum, predict deep binding. We find no evidence of deeply-bound tetraquarks, except in previously established strong-interaction stable $I=0$, $J^P=1^+$, $ud\bar{b}\bar{b}$ and $I=1/2$, $J^P=1^+$ $\ell s\bar{b}\bar{b}$ (where $\ell=u/d$) channels, allowing us to rule out models predicting deep binding. Preliminary results from an updated analysis of doubly-bottom tetraquarks are also presented.

Chiral dynamics in the low-temperature phase of QCD

We investigate the low-temperature phase of QCD and the crossover region with two light flavors of quarks. The chiral expansion around the point $(T,m=0)$ in the temperature vs. quark-mass plane indicates that a sharp real-time excitation exists with the quantum numbers of the pion. An exact sum rule is derived for the thermal modification of the spectral function associated with the axial charge density; the (dominant) pion pole contribution obeys the sum rule. We determine the two parameters of the pion dispersion relation using lattice QCD simulations and test the applicability of the chiral expansion. The time-dependent correlators are also analyzed using the Maximum Entropy Method, yie…

A lattice investigation of exotic tetraquark channels

We perform an $n_f=2+1$ lattice study of a number of channels where past claims exist in the literature for the existence of strong-interaction-stable light-heavy tetraquarks. We find no evidence for any such deeply-bound states, beyond the $J^P=1^+$, $I=0$ $ud\bar{b}\bar{b}$ and $I=1/2$ $ls\bar{b}\bar{b}$ states already identified in earlier lattice studies. We also describe a number of systematic improvements to our previous lattice studies, including working with larger $m_\pi L$ to better suppress possible finite volume effects, employing extended sinks to better control excited-state contamination, and expanding the number of operators used in the GEVP analyses. Our results also allow …

Study of the anomalous magnetic moment of the muon computed from the Adler function

We compute the Adler function on the lattice from vacuum polarization data with twisted boundary conditions using numerical derivatives. The study is based on CLS ensembles with two flavours of $O(a)$ improved Wilson fermions. We extrapolate the lattice data for the Adler function to the continuum limit and to the physical pion mass and analyze its dependence on the momentum transfer. We discuss the application of this method to the extraction of the $u,d$ contribution to $a_\mu^{\mathrm{HLO}}$.

Substructure of Multiquark Hadrons (Snowmass 2021 White Paper)

In recent years there has been a rapidly growing body of experimental evidence for existence of exotic, multiquark hadrons, i.e. mesons which contain additional quarks, beyond the usual quark-antiquark pair and baryons which consist of more than three quarks. In all cases with robust evidence they contain at least one heavy quark Q=c or b, the majority including two heavy quarks. Two key theoretical questions have been triggered by these discoveries: (a) how are quarks organized inside these multiquark states -- as compact objects with all quarks within one confinement volume, interacting via color forces, perhaps with an important role played by diquarks, or as deuteron-like hadronic molec…

Prospects for dilepton rates from lattice QCD

We discuss the prospects of computing thermal dilepton rates from first principles lattice QCD. The focus lies in the determination of the meson vector-vector current spectral function to estimate the electrical conductivity, heavy quark diffusion and quarkonium dissociation. We review and compare recent results from continuum-extrapolated, quenched calculations, as well as dynamical two-flavor setups.

Testing the strength of the UA(1) anomaly at the chiral phase transition in two-flavour QCD

We study the thermal transition of QCD with two degenerate light flavours by lattice simulations using $\mathcal{O}(a)$-improved Wilson quarks. Particular emphasis lies on the pattern of chiral symmetry restoration, which we probe via the static screening correlators. On $32^3$ volumes we observe that the screening masses in transverse iso-vector vector and axial-vector channels become degenerate at the transition temperature. The splitting between the screening masses in iso-vector scalar and pseudoscalar channels is strongly reduced compared to the splitting at zero temperature and is actually consistent with zero within uncertainties. In this proceedings article we extend our studies to …

Vector correlator and scale determination in lattice QCD

We implement a proposal made in [arXiv:1107.4388] to determine the lattice spacing by matching the lattice vector correlator at a reference distance scale with the same correlator obtained by a dispersion relation based on the $R$-ratio determined experimentally. We work with the isovector current, requiring a separation of the isovector hadronic final states on the phenomenological side. We also discuss the finite-size effect on the correlator, which must be controlled in order for the method to be applicable.

Charge transport and vector meson dissociation across the thermal phase transition in lattice QCD with two light quark flavors

We compute and analyze correlation functions in the isovector vector channel at vanishing spatial momentum across the deconfinement phase transition in lattice QCD. The simulations are carried out at temperatures $T/T_c=0.156, 0.8, 1.0, 1.25$ and $1.67$ with $T_c\simeq203$MeV for two flavors of Wilson-Clover fermions with a zero-temperature pion mass of $\simeq270$MeV. Exploiting exact sum rules and applying a phenomenologically motivated ansatz allows us to determine the spectral function $\rho(\omega,T)$ via a fit to the lattice correlation function data. From these results we estimate the electrical conductivity across the deconfinement phase transition via a Kubo formula and find eviden…

Thermal dilepton rates from quenched lattice QCD

We present new lattice results on the continuum extrapolation of the vector current correlation function. Lattice calculations have been carried out in the deconfined phase at a temperature of 1.1 Tc, extending our previous results at 1.45 Tc, utilizing quenched non-perturbatively clover-improved Wilson fermions and light quark masses. A systematic analysis on multiple lattice spacings allows to perform the continuum limit of the correlation function and to extract spectral properties in the continuum limit. Our current analysis suggests the results for the electrical conductivity are proportional to the temperature and the thermal dilepton rates in the quark gluon plasma are comparable for…

Vector screening masses in the quark–gluon plasma and their physical significance

Static and non-static thermal screening states that couple to the conserved vector current are investigated in the high-temperature phase of QCD. Their masses and couplings to the current are determined at weak coupling, as well as using two-flavor lattice QCD simulations. A consistent picture emerges from the comparison, providing evidence that non-static Matsubara modes can indeed be treated perturbatively. We elaborate on the physical significance of the screening masses.

Charmonium properties in hot quenched lattice QCD

We study the properties of charmonium states at finite temperature in quenched QCD on large and fine isotropic lattices. We perform a detailed analysis of charmonium correlation and spectral functions both below and above $T_c$. Our analysis suggests that both S wave states ($J/\psi$ and $\eta_c$) and P wave states ($\chi_{c0}$ and $\chi_{c1}$) disappear already at about $1.5 T_c$. The charm diffusion coefficient is estimated through the Kubo formula and found to be compatible with zero below $T_c$ and approximately $1/\pi T$ at $1.5 T_c\lesssim T\lesssim 3 T_c$.

The pion quasiparticle in the low-temperature phase of QCD

We extend our previous studies [PhysRevD.90.054509, PhysRevD.92.094510] of the pion quasiparticle in the low-temperature phase of two-flavor QCD with support from chiral effective theory. This includes the analysis performed on a finite temperature ensemble of size $20\times 64^3$ at $T\approx 151$MeV and a lighter zero-temperature pion mass $m_{\pi} \approx 185$ MeV. Furthermore, we investigate the Gell-Mann--Oakes-Renner relation at finite temperature and the Dey-Eletsky-Ioffe mixing theorem at finite quark mass.

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 {\pi}T/2 to 2{\pi}T. The relevant Euclidean vector current correlation functions are computed with $N_{\mathrm f}$ = 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 model-independent estimate from the Backus-Gilbert method to give an estimate for the p…

Investigating exotic heavy-light tetraquarks with 2+1 flavour lattice QCD

Static and non-static vector screening masses

Thermal screening masses of the conserved vector current are calculated both in a weak-coupling approach and in lattice QCD. The inverse of a screening mass can be understood as the length scale over which an external electric field is screened in a QCD medium. The comparison of screening masses both in the zero and non-zero Matsubara frequency sectors shows good agreement of the perturbative and the lattice results. Moreover, at $T\approx 508\mathrm{MeV}$ the lightest screening mass lies above the free result ($2\pi T$), in agreement with the $\mathcal{O}(g^2)$ weak-coupling prediction.

The leading disconnected contribution to the anomalous magnetic moment of the muon

The hadronic vacuum polarization can be determined from the vector correlator in a mixed time-momentum representation. We explicitly calculate the disconnected contribution to the vector correlator, both in the $N_f = 2$ theory and with an additional quenched strange quark, using non-perturbatively $O(a)$-improved Wilson fermions. All-to-all propagators are computed using stochastic sources and a generalized hopping parameter expansion. Combining the result with the dominant connected contribution, we are able to estimate an upper bound for the systematic error that arises from neglecting the disconnected contribution in the determination of $(g-2)_\mu$.

Adler function and hadronic vacuum polarization from lattice vector correlators in the time-momentum representation

We study a representation of the hadronic vacuum polarization based on the time-momentum representation of the vector correlator. This representation suggests a way to compute the hadronic vacuum polarization and the associated Adler function for any value of virtuality, irrespective of the flavor structure of the current. We present results on both of these phenomenologically important functions, derived from local-conserved two-point lattice vector correlation functions, computed on a subset of light two-flavor ensembles made available to us through the CLS effort.

On the strength of the U A (1) anomaly at the chiral phase transition in N f = 2 QCD

We study the thermal transition of QCD with two degenerate light flavours by lattice simulations using $O(a)$-improved Wilson quarks. Temperature scans are performed at a fixed value of $N_t = (aT)^{-1}=16$, where $a$ is the lattice spacing and $T$ the temperature, at three fixed zero-temperature pion masses between 200 MeV and 540 MeV. In this range we find that the transition is consistent with a broad crossover. As a probe of the restoration of chiral symmetry, we study the static screening spectrum. We observe a degeneracy between the transverse isovector vector and axial-vector channels starting from the transition temperature. Particularly striking is the strong reduction of the split…

A new representation of the Adler function for lattice QCD

We address several aspects of lattice QCD calculations of the hadronic vacuum polarization and the associated Adler function. We implement a representation derived previously which allows one to access these phenomenologically important functions for a continuous set of virtualities, irrespective of the flavor structure of the current. Secondly we present a theoretical analysis of the finite-size effects on our particular representation of the Adler function, based on the operator product expansion at large momenta and on the spectral representation of the Euclidean correlator at small momenta. Finally, an analysis of the flavor structure of the electromagnetic current correlator is perform…

Adler function and hadronic vacuum polarization from lattice vector correlators

We study a representation of the hadronic vacuum polarization based on the time-momentum representation of the vector correlator. This representation suggests a way to compute the hadronic vacuum polarization and the associated Adler function for any value of virtuality, irrespective of the flavor structure of the current. We present results on both of these phenomenologically important functions, derived from local-conserved two-point lattice vector correlation functions, computed on a subset of light two-flavor ensembles made available to us through the CLS effort.

QCD thermodynamics with two flavours of Wilson fermions on large lattices

We explore the phase diagram of two flavour QCD at vanishing chemical potential using dynamical O(a)-improved Wilson quarks. In the approach to the chiral limit we use lattices with a temporal extent of N_t=16 and spatial extent L=32,48 and 64 to enable the extrapolation to the thermodynamic limit with small discretisation effects. In addition to an update on the scans at constant \kappa, reported earlier, we present first results from scans along lines of constant physics at a pion mass of 290 MeV. We probe the transition using the Polyakov loop and the chiral condensate, as well as spectroscopic observables such as screening masses.

A lattice calculation of the hadronic vacuum polarization contribution to (g - 2)μ

We present results of calculations of the hadronic vacuum polarisation contribution to the muon anomalous magnetic moment. Specifically, we focus on controlling the infrared regime of the vacuum polarisation function. Our results are corrected for finite-size effects by combining the Gounaris-Sakurai parameterisation of the timelike pion form factor with the Lüscher formalism. The impact of quark-disconnected diagrams and the precision of the scale determination is discussed and included in our final result in two-flavour QCD, which carries an overall uncertainty of 6%. We present preliminary results computed on ensembles with Nf = 2 + 1 dynamical flavours and discuss how the long-distance …

Charmonium dissociation and heavy quark transport in hot quenched lattice QCD

We study the properties of charmonium states at finite temperature in quenched lattice QCD on large and fine isotropic lattices. We perform a detailed analysis of charmonium correlation and spectral functions both below and above Tc. Our analysis suggests that the S wave states disappear at about 1.5 Tc. The charm diffusion coefficient is estimated and found to be approximately 1/{\pi}T at 1.5Tc {\leq} T {\leq} 3Tc.

Lattice QCD study of the $H$ dibaryon using hexaquark and two-baryon interpolators

Physical review / D 99(7), 074505 (2019). doi:10.1103/PhysRevD.99.074505

The pion quasiparticle in the low-temperature phase of QCD

We investigate the properties of the pion quasiparticle in the low-temperature phase of two-flavor QCD on the lattice with support from chiral effective theory. We find that the pion quasiparticle mass is significantly reduced compared to its value in the vacuum, by contrast with the static screening mass, which increases with temperature. By a simple argument, near the chiral limit the two masses are expected to determine the quasiparticle dispersion relation. Analyzing two-point functions of the axial charge density at non-vanishing spatial momentum, we find that the predicted dispersion relation and the residue of the pion pole are simultaneously consistent with the lattice data at low m…

Lattice QCD estimate of the quark-gluon plasma photon emission rate

We present a computation of the photon emission rate of the quark-gluon plasma from two-flavor lattice QCD at a temperature of 254 MeV, which follows up on the work presented in [1]. We perform a continuum extrapolation of the vector-current correlator, and consider a linear combination of the Lorentz indices corresponding to a UV-finite spectral function. To extract the spectral function from the lattice correlators, an ill-posed inverse problem, we model the spectral function with a Pad\'e ansatz. We further constrain our analysis by simultaneously fitting data with different momenta. We present results for a multi-momentum fit including the three smallest momenta available from our latti…

A relation between screening masses and real-time rates

Thermal screening masses related to the conserved vector current are determined for the case that the current carries a non-zero Matsubara frequency, both in a weak-coupling approach and through lattice QCD. We point out that such screening masses are sensitive to the same infrared physics as light-cone real-time rates. In particular, on the perturbative side, the inhomogeneous Schrodinger equation determining screening correlators is shown to have the same general form as the equation implementing LPM resummation for the soft-dilepton and photon production rates from a hot QCD plasma. The static potential appearing in the equation is identical to that whose soft part has been determined up…

Antiscreening of the Ampère force in QED and QCD plasmas

The static forces between electric charges and currents are modified at the loop level by the presence of a plasma. While electric charges are screened, currents are not. The effective coupling constant at long distances is enhanced in both cases as compared to the vacuum, and by different amounts, a clear sign that Lorentz symmetry is broken. We investigate these effects quantitatively, first in a QED plasma and secondly using non-perturbative simulations of QCD with two light degenerate flavors of quarks.

Two-flavour lattice QCD correlation functions in the deconfinement transition region

We report on a lattice QCD calculation with two dynamical flavors of the isovector vector correlator in the high-temperature phase. We analyze the correlator in terms of the associated spectral function by performing a fit for the difference of the thermal and vacuum spectral functions, using also an exact sum rule that constrains this difference. Additonally we carry out a direct fit for the thermal spectral function, and obtain good agreement between the two analyses for frequencies below the two-pion threshold. Under the assumption that the spectral function is smooth in that region, we give an estimate of the electrical conductivity.

On the temperature dependence of the electrical conductivity in hot quenched lattice QCD

Extending our recent work, we report on a calculation of the vector current correlation function for light valence quarks in the deconfined phase of quenched QCD in the temperature range 1.16Tc<T<2.98Tc. After performing a systematic analysis of the in fluence of cut-off effects on light quark meson correlators using clover improved Wilson fermions, we discuss resulting constraints on the electrical conductivity in a quark gluon plasma.

The leading hadronic contribution to (g-2) of the muon: The chiral behavior using the mixed representation method

We extend our analysis of the leading hadronic contribution to the anomalous magnetic moment of the muon using the mixed representation method to study its chiral behavior. We present results derived from local-conserved two-point lattice vector correlation functions, computed on a subset of light two-flavor ensembles made available to us through the CLS effort with pion masses as low as 190 MeV. The data is analyzed also using the more standard four-momentum method. Both methods are systematically compared as the calculations approach the physical point.