Search

Using lattice QCD simulations with $O(a)$-improved Wilson quarks and physical up, down and strange quark masses, we investigate the properties of thermal QCD matter at the temperatures $T=\{128,154,192\}\mathrm{MeV}$ with a fixed lattice spacing $a=0.064\mathrm{fm}$ and volume $V=(6.1\text{fm})^3$. We find that the pion quasiparticle, defined as the low-energy pole in the two-point function of the axial charge, becomes lighter as the temperature increases and give an argument based on hydrodynamics as to why the pole becomes purely diffusive above the chiral crossover. We study the thermal modification of the isovector vector spectral function using the Backus-Gilbert method, finding an enhancement at low energies and a depletion at energies around $1\mathrm{GeV}$.The analogous study of the axial-vector channel reveals a larger enhancement at energies below $1\mathrm{GeV}$, and we show that these findings are consistent with rigorous spectral sum rules. The difference between vector and axial-vector correlators, an order parameter for chiral symmetry, turns out to be overall suppressed by more than an order of magnitude at the crossover., Comment: 47 pages, 20 figures. arXiv admin note: text overlap with arXiv:2211.15558

The baryon masses on CLS ensembles are used to determine the LEC that characterize QCD in the flavor-SU(3) limit with vanishing up, down, and strange quark masses. Here we reevaluate some of the baryon masses on flavor-symmetric ensembles with much-improved statistical precision, in particular for the decuplet states. These additional results then lead to a more significant chiral extrapolation of the Lattice data set to its chiral SU(3) limit. Our results are based on the chiral Lagrangian with baryon octet and decuplet fields considered at the one-loop level. Finite-box and discretization effects of the Lattice data are considered systematically. While in our global fit of the data we insist on large-Nc sum rules for the LEC that enter at N3LO, all other LEC are unconstrained. In particular, we obtain values for the chiral limit of the pion decay constant and the isospin-limit of the quark-mass ratio compatible with the FLAG report., Comment: 30 pages, 6 figures and 11 tables, further two figures and tables. arXiv admin note: text overlap with arXiv:2301.06837

In this paper we determine the masses of $I(J^P)=0\left(3/2^+\right)$ and $0\left(3/2^-\right)$ $\Omega$-baryon ground states using lattice QCD. We utilise Wilson-clover ensembles with $2+1$ dynamical quark flavours generated by the CLS consortium along a trajectory with a constant trace of the quark-mass matrix. We show that N$^3$LO $\text{SU}(3)_f$ chiral perturbation theory expressions describe the ground-state masses with positive-parity well, and we use them to set the lattice scale. Methodologically, our combination of gauge-fixed wall sources and the generalized Pencil of Functions allows for high-precision determinations of the lattice spacing at a relative error of around $0.3\%$ with controlled excited-state contamination. The fit we perform allows for the continuum value of $t_0$ to vary, thereby determining this quantity with a comparable level of precision to that of the lattice scale. Using the resulting scales our measurement of the negative-parity $\Omega^{3/2^{-}}$ state is found to be consistent with the recently-discovered $\Omega(2012)^-$, which can therefore be assigned the quantum numbers $I(J^P)=0\left(3/2^-\right)$., Comment: 29 pages, 7 figures

The thermal photon emissivity at the QCD chiral crossover is investigated using imaginary momentum correlators. These have been measured on a newly generated $20 \times 96^3$ lattice-QCD ensemble with $\mathcal{O}(a)$-improved Wilson quarks and physical up, down and strange quark masses at a temperature $T=154$\,MeV near the pseudo-critical temperature. In order to realize the photon on-shell condition, the spatially transverse Euclidean correlators have to be evaluated at imaginary spatial momenta. Employing a bounding method, we present a preliminary result on the quantity $H_E(\omega_1)$, which corresponds to an energy-moment of the photon spectral function $\sigma(\omega)/\omega$ defined by the weight function $1/(\omega^2+\omega_n^2)$, the $\omega_n$ being integer multiples of $2\pi T$., Comment: 12 pages, 4 Figures, 1 Table

Following the publication of the new measurement of the anomalous magnetic moment of the muon, the discrepancy between experiment and the theory prediction from the $g-2$ theory initiative has increased to $4.2\,\sigma$. Recent lattice QCD calculations predict values for the hadronic vacuum polarization contribution that are larger than the data-driven estimates, bringing the Standard Model prediction closer to the experimental measurement. Euclidean time windows in the time-momentum representation of the hadronic vacuum polarization contribution to the muon $g-2$ can help clarify the discrepancy between the phenomenological and lattice predictions. We present our calculation of the intermediate distance window contribution using $N_\mathrm{f}=2+1$ flavors of O$(a)$ improved Wilson quarks. We employ ensembles at six lattice spacings below $0.1\,$fm and pion masses down to the physical value. We present a detailed study of the continuum limit, using two discretizations of the vector current and two independent sets of improvement coefficients. Our result at the physical point displays a tension of $3.9\,\sigma$ with a recent evaluation of the intermediate window based on the data-driven method., Comment: 10 pages, 5 figures, Proceedings of the 39th International Symposium on Lattice Field Theory, 8th-13th August 2022, Bonn, Germany

Hadronic light-by-light scattering is one of the virtual processes that causes the gyromagnetic factor $g$ of the muon to deviate from the value of two predicted by Dirac's theory. This process makes one of the largest contributions to the uncertainty of the Standard Model prediction for the muon $(g-2)$. Lattice QCD allows for a first-principles approach to computing this non-perturbative effect. In order to avoid power-law finite-size artifacts generated by virtual photons in lattice simulations, we follow a coordinate-space approach involving a weighted integral over the vertices of the QCD four-point function of the electromagnetic current carried by the quarks. Here we present in detail the semi-analytical calculation of the QED part of the amplitude, employing position-space perturbation theory in continuous, infinite four-dimensional Euclidean space. We also provide some useful information about a computer code for the numerical implementation of our approach that has been made public at https://github.com/RJHudspith/KQED., Comment: 84 pages, 13 figures; text matches the published version

Euclidean time windows in the integral representation of the hadronic vacuum polarization contribution to the muon $g-2$ serve to test the consistency of lattice calculations and may help in tracing the origins of a potential tension between lattice and data-driven evaluations. In this paper, we present results for the intermediate time window observable computed using O($a$) improved Wilson fermions at six values of the lattice spacings below 0.1\,fm and pion masses down to the physical value. Using two different sets of improvement coefficients in the definitions of the local and conserved vector currents, we perform a detailed scaling study which results in a fully controlled extrapolation to the continuum limit without any additional treatment of the data, except for the inclusion of finite-volume corrections. To determine the latter, we use a combination of the method of Hansen and Patella and the Meyer-Lellouch-L\"uscher procedure employing the Gounaris-Sakurai parameterization for the pion form factor. We correct our results for isospin-breaking effects via the perturbative expansion of QCD+QED around the isosymmetric theory. Our result at the physical point is $a_\mu^{\mathrm{win}}=(237.30\pm0.79_{\rm stat}\pm1.22_{\rm syst})\times10^{-10}$, where the systematic error includes an estimate of the uncertainty due to the quenched charm quark in our calculation. Our result displays a tension of 3.9$\sigma$ with a recent evaluation of $a_\mu^{\mathrm{win}}$ based on the data-driven method., Comment: 43 pages, 9 figures, 10 tables; version accepted for publication: extended discussion of finite-volume corrections. Results and conclusions unchanged

We compute the hadronic light-by-light scattering contribution to the muon $g-2$ from the charm quark using lattice QCD. The calculation is performed on ensembles generated with dynamical $(u,d,s)$ quarks at the SU(3)$_{\rm f}$ symmetric point with degenerate pion and kaon masses of around 415 MeV. It includes the connected charm contribution, as well as the leading disconnected Wick contraction, involving the correlation between a charm and a light-quark loop. Cutoff effects turn out to be sizeable, which leads us to use lighter-than-physical charm masses, to employ a broad range of lattice spacings reaching down to 0.039 fm and to perform a combined charm-mass and continuum extrapolation. We use the $\eta_c$ meson to define the physical charm-mass point and obtain a final value of $a_\mu^{\rm HLbL,c} = (2.8\pm 0.5) \times 10^{-11}$, whose uncertainty is dominated by the systematics of the extrapolation. Our result is consistent with the estimate based on a simple charm-quark loop, whilst being free of any perturbative scheme dependence on the charm mass. The mixed charm-light disconnected contraction contributes a small negative amount to the final value., Comment: 21 pages, 8 figures, 9 tables

We present work designed to compute baryon masses on $N_f = 2 + 1$ CLS ensembles including isospin-breaking effects due to non-degenerate light quark masses and electromagnetic interactions. These effects are determined at leading order via a perturbative expansion around the iso-symmetric theory. We furthermore apply a group-theoretical operator construction for the various interpolators describing the different members of the baryon octet and decuplet based on a classification by spin, parity, and flavor content.

We compute the hadronic light-by-light scattering contribution to the muon $g-2$ from the up, down, and strange-quark sector directly using lattice QCD. Our calculation features evaluations of all possible Wick-contractions of the relevant hadronic four-point function and incorporates several different pion masses, volumes, and lattice-spacings. We obtain a value of $a_\mu^{\text{Hlbl}} = 106.8(14.7) \times 10^{-11}$ (adding statistical and systematic errors in quadrature), which is consistent with current phenomenological estimates and a previous lattice determination. It now appears conclusive that the hadronic light-by-light contribution cannot explain the current tension between theory and experiment for the muon $g-2$., Comment: 42 pages, 21 figures

We perform a lattice QCD calculation of the hadronic light-by-light contribution to $(g-2)_\mu$ at the SU(3) flavor-symmetric point $m_\pi=m_K\simeq 420\,$MeV. The representation used is based on coordinate-space perturbation theory, with all QED elements of the relevant Feynman diagrams implemented in continuum, infinite Euclidean space. As a consequence, the effect of using finite lattices to evaluate the QCD four-point function of the electromagnetic current is exponentially suppressed. Thanks to the SU(3)-flavor symmetry, only two topologies of diagrams contribute, the fully connected and the leading disconnected. We show the equivalence in the continuum limit of two methods of computing the connected contribution, and introduce a sparse-grid technique for computing the disconnected contribution. Thanks to our previous calculation of the pion transition form factor, we are able to correct for the residual finite-size effects and extend the tail of the integrand. We test our understanding of finite-size effects by using gauge ensembles differing only by their volume. After a continuum extrapolation based on four lattice spacings, we obtain $a_\mu^{\rm hlbl} = (65.4\pm 4.9 \pm 6.6)\times 10^{-11}$, where the first error results from the uncertainties on the individual gauge ensembles and the second is the systematic error of the continuum extrapolation. Finally, we estimate how this value will change as the light-quark masses are lowered to their physical values., Comment: 19 figures, 39 pages; improved references, in particular concerning the eta exchange; no figures or results changed

The measurement of the anomalous magnetic moment of the muon and its prediction allow for a high-precision test of the Standard Model (SM). In this proceedings article we present ongoing work combining lattice QCD and continuum QED in order to determine an important SM contribution to the magnetic moment, the hadronic light-by-light contribution. We compute the quark-connected contribution in the Mainz position-space approach and investigate the long-distance part of our data using calculations of the $\pi^0$-pole and charged pion loop contributions., Comment: 7 pages, 10 figures

We continue our study of heavy-light four-quark states and find evidence from lattice QCD for the existence of a strong-interaction-stable $I(J^P)=0(1^+)$ $ud\bar{c}\bar{b}$ tetraquark with mass in the range of 15 to 61 MeV below $\bar{D}B^*$ threshold. Since this range includes the electromagnetic $\bar{D}B\gamma$ decay threshold, current uncertainties do not allow us to determine whether such a state would decay electromagnetically, or only weakly. We also perform a study at fixed pion mass, with NRQCD for the heavy quarks, simulating $qq^\prime \bar{b}^\prime \bar{b}$ and $q q^\prime \bar{b}^\prime\bar{b}^\prime$ tetraquarks with $q,\, q^\prime =ud$ or $\ell s$ and variable, unphysical $m_{b^\prime}$ in order to investigate the heavy mass-dependence of such tetraquark states. We find that the dependence of the binding energy follows a phenomenologically-expected form and that, though NRQCD breaks down before $m_{b^\prime}=m_c$ is reached, the results at higher $m_{b^\prime}$ clearly identify the $ud\bar{b}^\prime \bar{b}$ channel as the most likely to support a strong-interaction-stable tetraquark state at $m_{b^\prime}=m_c$. This observation serves to motivate the direct $ud\bar{c}\bar{b}$ simulation. Throughout we use dynamical $n_f=2+1$ ensembles with pion masses $m_\pi=$415, 299, and 164 MeV reaching down almost to the physical point, a relativistic heavy quark prescription for the charm quark, and NRQCD for the bottom quark(s)., Comment: 24 pages, 4 figures

As a simple model for dark matter, we propose a QCD-like theory based on $\rm{SU}(2)$ gauge theory with one flavor of dark quark. The model is confining at low energy and we use lattice simulations to investigate the properties of the lowest-lying hadrons. Compared to QCD, the theory has several peculiar differences: there are no Goldstone bosons or chiral symmetry restoration when the dark quark becomes massless; the usual global baryon number symmetry is enlarged to $\rm{SU}(2)_B$, resembling isospin; and baryons and mesons are unified together in $\rm{SU}(2)_B$ iso-multiplets. We argue that the lightest baryon, a vector boson, is a stable dark matter candidate and is a composite realization of the hidden vector dark matter scenario. The model naturally includes a lighter state, the analog of the $\eta^\prime$ in QCD, for dark matter to annihilate into to set the relic density via thermal freeze-out. Dark matter baryons may also be asymmetric, strongly self-interacting, or have their relic density set via $3 \to 2$ cannibalizing transitions. We discuss some experimental implications of coupling dark baryons to the Higgs portal., Comment: 26 pages, 16 figures

We present a determination of the QCD coupling constant, $\alpha_s$, obtained by fitting lattice results for the flavour $ud$ hadronic vacuum polarisation function to continuum perturbation theory. We use $n_f=2+1$ flavours of Domain Wall fermions generated by the RBC/UKQCD collaboration and three lattice spacings $a^{-1}=1.79,\, 2.38$ and $3.15\ \text{GeV}$. Several sources of potential systematic error are identified and dealt with. After fitting and removing expected leading cut-off effects, we find for the five-flavour $\overline{\text{MS}}$ coupling the value $\alpha_s(M_Z)=0.1181(27)^{+8}_{-22}$, where the first error is statistical and the second systematic.

We report on our progress in studying exotic, heavy tetraquark states, $qq\prime \bar Q\bar Q\prime$. Using publicly available dynamical $n_f =2+1$ Wilson-Clover gauge configurations, generated by the PACS-CS collaboration, with pion masses $\simeq$164, 299 and 415 MeV, we extend our previous analysis to heavy quark components containing heavier than physical bottom quarks $\bar Q\bar Q\prime=\bar b\prime\bar b\prime$ or $\bar Q\bar Q\prime=\bar b\bar b\prime$, charm and bottom quarks $\bar c\bar b$ and also only charm quarks $\bar c\bar c$. Throughout we employ NRQCD and relativistic heavy quarks for the heavier than bottom, bottom and charm quarks. Using our previously established diquark-antidiquark and meson-meson operator basis we comment in particular on the dependence of the binding energy on the mass of the heavy quark component $\bar Q\bar Q$, with heavy quarks ranging from $m_Q=0.85\ldots 6.3\cdot m_b$. In the heavy flavor non-degenerate case, $\bar Q\bar Q\prime$, and especially for the tetraquark channel $ud\bar c\bar b$, we extend our work to utilize a $3\times 3$ GEVP to study the ground and threshold states thereby enabling a clear identification of possible binding. Finally, we present initial work on the $\bar Q\bar Q\prime=\bar c\bar c$ system where a much larger operator basis is available in comparison to flavor combinations with NRQCD quarks., Comment: 8 pages, 5 figures, proceedings contribution to "Lattice 2017. 35th International Symposium on Lattice Field Theory", 18th-24th June 2017, Granada, Spain

We revisit the puzzle of $|V_{us}|$ values obtained from the conventional implementation of hadronic-$\tau$-decay-based flavor-breaking finite-energy sum rules lying $>3 \sigma$ below the expectations of three-family unitarity. Significant unphysical dependences of $|V_{us}|$ on the choice of weight, w, and upper limit, $s_0$, of the experimental spectral integrals entering the analysis are confirmed, and a breakdown of assumptions made in estimating higher dimension, $D>4$, OPE contributions is identified as the main source of these problems. A combination of continuum and lattice results is shown to suggest a new implementation of the flavor-breaking sum rule approach in which not only $|V_{us}|$, but also $D>4$ effective condensates, are fit to data. Lattice results are also used to clarify how to reliably treat the slowly converging D=2 OPE series. The new sum rule implementation is shown to cure the problems of the unphysical w- and $s_0$-dependence of $|V_{us}|$ and to produce results $\sim 0.0020$ higher than those of the conventional implementation. With B-factory input, and using, in addition, dispersively constrained results for the $K\pi$ branching fractions, we find $\vert V_{us}\vert =0.2231(27)_{exp}(4)_{th}$, in excellent agreement with the result from $K_{\ell 3}$, and compatible within errors with the expectations of three-family unitarity, thus resolving the long-standing inclusive $\tau$ $|V_{us}|$ puzzle., Comment: 15 pages, 4 figures. Updated treatment of strange hadronic tau decay data; expanded discussion of the breakdown of the conventional implementation approach. Final version to appear in Physics Letters B

We compute the hadronic matrix elements of the four-quark operators relevant for $K^0-{\bar K^0}$ mixing beyond the Standard Model. Our results are from lattice QCD simulations with $n_f=2+1$ flavours of domain-wall fermion, which exhibit continuum-like chiral-flavour symmetry. The simulations are performed at two different values of the lattice spacing ($a\sim0.08$ and $a\sim 0.11 \, \fm $) and with lightest unitary pion mass $\sim 300\, \MeV$. For the first time, the full set of relevant four-quark operators is renormalised non-perturbatively through RI-SMOM schemes; a detailed description of the renormalisation procedure is presented in a companion paper. We argue that the intermediate renormalisation scheme is responsible for the discrepancies found by different collaborations. We also study different normalisations and determine the matrix elements of the relevant four-quark operators with a precision of $\sim 5\%$ or better., Comment: 38 pages

We investigate the possibility of $qq^\prime \bar b \bar b$ tetraquark bound states using $n_f=2+1$ lattice QCD ensembles with pion masses $\simeq 164$, $299$, and $415$ MeV. Motivated by observations from heavy baryon phenomenology, we consider two lattice interpolating operators both of which are expected to couple efficiently to tetraquark states: one with diquark-antidiquark and one with a meson-meson structure. Using nonrelativistic QCD to simulate the bottom quarks, we study the $ud\bar b \bar b$, $\ell s\bar b \bar b$ channels with $\ell=u,d$, and find unambiguous signals for strong-interaction-stable $J^P=1^+$ tetraquarks. These states are found to lie $189(10)$ and $98(7)$ MeV below the corresponding free two-meson thresholds., Comment: 5 pages, 3 figures; published version

The QCD coupling appears in the perturbative expansion of the current-current two-point (vacuum polarization) function. Any lattice calculation of vacuum polarization is plagued by several competing non-perturbative effects at small momenta and by discretization errors at large momenta. We work in an intermediate region, computing the vacuum polarization for many off-axis momentum directions on the lattice. Having many momentum directions provides a way to monitor and account for lattice artifacts. Our results are competitive with, and have certain systematic advantages over, the alternate phenomenological determination of the strong coupling from the same light quark vacuum polarization produced by sum rule analyses of hadronic tau decay data., Comment: 7 pages, 4 figures, Proceedings of the 33rd International Symposium on Lattice Field Theory (Lattice 2015), 14-18 July 2015, Kobe International Conference Center, Kobe, Japan

We report on our computation of the hadronic matrix elements of the four-quark operators needed for the study of $K^0-{\bar K^0}$ mixing beyond the Standard Model (SM). We consider nf=2+1 Domain-Wall fermions on Iwasaki gauge action with lightest unitary pion of 290 MeV and a single lattice spacing a=0.086 fm. The renormalization is performed non-perturbatively through the RI-MOM scheme and our results are converted perturbatively to MSbar. We have estimated the various systematic errors. Our results confirm a previous quenched study, where large ratios of non-SM to SM matrix elements were obtained., Comment: 7 pages, 3 figures. Contribution to the 30th International Symposium on Lattice Field Theory, June 24-29, 2012, Cairns, Australia