Search

The low-lying energy spectrum of the static-colour-source-anti-source system in a vacuum containing light and strange quarks is computed using lattice QCD for a range of different light quark masses. The resulting levels are described using a simple model Hamiltonian and the parameters in this model are extrapolated to the physical light-quark masses. In this framework, the QCD string tension is found to be $\sqrt{\sigma}=445(3)_{\rm stat}(6)_{\rm sys}$ MeV., Comment: 10 pages, 9 figures, small changes to match published version

Ongoing challenges in computing the spectrum of hadronic resonances and shallow bound-states from lattice QCD are reviewed. Since such states are identified as poles in the scattering matrix, nearby non-analyticities must be treated to analytically continue to complex center-of-mass energies. Significant lattice spacing effects have also been observed in some channels, necessitating a continuum limit. Recent achievements are also highlighted, including lattice investigations of states in the charm region, baryon-baryon scattering, and the first coupled channel meson-baryon amplitude in the $\Lambda(1405)$ channel., Comment: 6 pages, 4 figures, proceedings from talk at MENU 2023

Studies of the $\Delta$ baryon resonance and the $K_0^\ast(700)$ and $a_0(980)$ meson resonances using $N_f=2+1$ lattice QCD for pion masses near 200 MeV are presented. The $s$-wave scattering lengths for both the $I=1/2$ $N \pi$ and $I=3/2$ $N \pi$ channels and properties of the $\Delta$ resonance are identified from the finite-volume energy levels of the lattice simulation. The importance of a three-quark $\Delta$-operator in the $N\pi$ system and tetraquark operators in the mesonic systems is investigated., Comment: 10 pages, 5 figures, 3 tables

This work presents technical details of determining the finite-volume energy spectra for the scattering amplitude of the coupled-channel $\pi\Sigma - \bar{K}N$ from lattice QCD data. The importance of reliably extracting such spectra lies in the crucial dependence of the hadronic scattering amplitudes analysis on the energy spectrum when using L\"{u}scher's formalism. Results of the methods used are presented and the final finite-volume spectra are shown. The analysis of the scattering amplitude based on these results, exhibits a two-pole structure for the $\Lambda(1405)$, a virtual bound state below the $\pi\Sigma$ threshold and a resonance pole right below the $\bar{K}N$ threshold.

A novel application of lattice QCD spectral reconstruction is presented, in which euclidean correlation function data in a fixed time range are used to infer values outside the range, enabling a model-independent investigation of the asymptotic large-time behavior. Moments of the correlator are also determined, and reconstructed correlation matrices between different moments are included in a variational optimization similar to the standard Generalized Eigenvalue Problem (GEVP). These ideas are illustrated using a single-nucleon correlation function determined on an $N_{\rm}=2+1$ ensemble of gauge configurations at $m_{\pi} = 200{\rm MeV}$., Comment: 5 pages, 3 figures, proceedings from the HADRON2023 conference

Recent results studying the masses and widths of low-lying baryon resonances in lattice QCD are presented. The $S$-wave $N\pi$ scattering lengths for both total isospins $I = 1/2$ and $I = 3/2$ are inferred from the finite-volume spectrum below the inelastic threshold together with the $I = 3/2$ $P$-wave containing the $\Delta(1232)$ resonance. A lattice QCD computation employing a combined basis of three-quark and meson-baryon interpolating operators with definite momentum to determine the coupled channel $\Sigma\pi$-$N\overline{K}$ scattering amplitude in the $\Lambda(1405)$ region is also presented. Our results support the picture of a two-pole structure suggested by theoretical approaches based on $SU(3)$ chiral symmetry and unitarity., Comment: 5 pages, 3 figures, talk presented at the 20th International Conference on Hadron Spectroscopy and Structure (HADRON 2023) in Genova, Italy, June 5th to 9th 2023

A lattice QCD computation of the coupled channel $\pi\Sigma-\bar{K}N$ scattering amplitudes in the $\Lambda(1405)$ region is detailed. Results are obtained using a single ensemble of gauge field configurations with $N_{\rm f} = 2+1$ dynamical quark flavors and $m_{\pi} \approx 200$ MeV and $m_K\approx487$ MeV. Hermitian correlation matrices using both single baryon and meson-baryon interpolating operators for a variety of different total momenta and irreducible representations are used. Several parametrizations of the two-channel scattering $K$-matrix are utilized to obtain the scattering amplitudes from the finite-volume spectrum. The amplitudes, continued to the complex energy plane, exhibit a virtual bound state below the $\pi\Sigma$ threshold and a resonance pole just below the $\bar{K}N$ threshold., Comment: 18 pages, 11 figures, 14 tables. Final version accepted for publication in PRD with post-publication corrections

This letter presents the first lattice QCD computation of the coupled channel $\pi\Sigma-\bar{K}N$ scattering amplitudes at energies near $1405\,{\rm MeV}$. These amplitudes contain the resonance $\Lambda(1405)$ with strangeness $S=-1$ and isospin, spin, and parity quantum numbers $I(J^P)=0(1/2^-)$. However, whether there is a single resonance or two nearby resonance poles in this region is controversial theoretically and experimentally. Using single-baryon and meson-baryon operators to extract the finite-volume stationary-state energies to obtain the scattering amplitudes at slightly unphysical quark masses corresponding to $m_\pi\approx200$ MeV and $m_K\approx487$ MeV, this study finds the amplitudes exhibit a virtual bound state below the $\pi\Sigma$ threshold in addition to the established resonance pole just below the $\bar{K}N$ threshold. Several parametrizations of the two-channel $K$-matrix are employed to fit the lattice QCD results, all of which support the two-pole picture suggested by $SU(3)$ chiral symmetry and unitarity., Comment: 7 pages, 4 figures, and 1 table. Includes various corrections made during the last stages before publication

Substantial progress has been made recently in the generation of master-field ensembles. This has to be paired with efficient techniques to compute observables on gauge field configurations with a large volume. Here we present the results of the computation of hadronic observables, including hadron masses and meson decay constants, on large-volume and master-field ensembles with physical volumes of up to $(18\,\mathrm{fm})^4$ and $m_\pi L$ up to $25$, simulated using $N_{\mathrm{f}}=2+1$ stabilized Wilson fermions. We obtain sub-percent determinations from single gauge configurations with the combined use of position-space techniques, volume averages and master-field error estimation., Comment: 10 pages, 4 figures, 2 tables, talk presented at The 39th International Symposium on Lattice Field Theory, 8th-13th August, 2022, Rheinische Friedrich-Wilhelms-Universit\"at Bonn, Bonn, Germany

A recently re-discovered variant of the Backus-Gilbert algorithm for spectral reconstruction enables the controlled determination of smeared spectral densities from lattice field theory correlation functions. A particular advantage of this approach is the \emph{a priori} specification of the kernel with which the underlying spectral density is smeared, allowing for variation of its peak position, smearing width, and functional form. If the unsmeared spectral density is sufficiently smooth in the neighborhood of a particular energy, it can be obtained from an extrapolation to zero smearing-kernel width at fixed peak position. A natural application for this approach is scattering processes summed over all hadronic final states. As a proof-of-principle test, an inclusive rate is computed in the two-dimensional O(3) sigma model from a two-point correlation function of conserved currents. The results at finite and zero smearing radius are in good agreement with the known analytic form up to energies at which 40-particle states contribute, and are sensitive to the 4-particle contribution to the inclusive rate. The straight-forward adaptation to compute the $R$-ratio in lattice QCD from two-point functions of the electromagnetic current is briefly discussed., Comment: 16 pages, 10 figures, proceedings from the Lattice 2022 conference

Elastic nucleon-pion scattering amplitudes are computed using lattice QCD on a single ensemble of gauge field configurations with $N_{\rm f} = 2+1$ dynamical quark flavors and $m_{\pi} = 200~{\rm MeV}$. The $s$-wave scattering lengths with both total isospins $I=1/2$ and $I=3/2$ are inferred from the finite-volume spectrum below the inelastic threshold together with the $I=3/2$ $p$-wave containing the $\Delta(1232)$ resonance. The amplitudes are well-described by the effective range expansion with parameters constrained by fits to the finite-volume energy levels enabling a determination of the $I=3/2$ scattering length with statistical errors below $5\%$, while the $I = 1/2$ is somewhat less precise. Systematic errors due to excited states and the influence of higher partial waves are controlled, providing a pathway for future computations down to the physical light quark masses with multiple lattice spacings and physical volumes., Comment: 33 pages, 16 figures; version accepted for publication

The status and prospects for investigations of exotic and conventional hadrons with lattice QCD are discussed. The majority of hadrons decay strongly via one or multiple decay-channels, including most of the experimentally discovered exotic hadrons. Despite this difficult challenge, the properties of several hadronic resonances have been determined within lattice QCD. To further discern the spectroscopic properties of various hadrons and to help resolve their nature we present our suggestions for future analytic and lattice studies., Comment: 23 pages + references, Submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021); v2: few references added

Lattice QCD calculations of two-nucleon interactions have been underway for about a decade, but still haven't reached the pion mass regime necessary for matching onto effective field theories and extrapolating to the physical point. Furthermore, results from different methods, including the use of the Luscher formalism with different types of operators, as well as the HALQCD potential method, do not agree even qualitatively at very heavy pion mass. We investigate the role that different operators employed in the literature may play on the extraction of spectra for use within the Luscher method. We first explore expectations from Effective Field Theory solved within a finite volume, for which the exact spectrum may be computed given different physical scenarios. We then present preliminary lattice QCD results for two-nucleon spectra calculated using different operators on a common lattice ensemble., Comment: 11 pages, 7 figures, proceeding for The 38th International Symposium on Lattice Field Theory, LATTICE2021, 26th-30th July, 2021, Zoom/Gather@Massachusetts Institute of Technology

This work employs the spectral reconstruction approach of Ref. [1] to determine an inclusive rate in the $1+1$ dimensional O(3) non-linear $\sigma$-model, analogous to the QCD part of ${e}^+{e}^- \rightarrow \rm {hadrons}$. The Euclidean two-point correlation function of the conserved current $j$ is computed using Monte Carlo lattice field theory simulations for a variety of spacetime volumes and lattice spacings. The spectral density of this correlator is related to the inclusive rate for $j \rightarrow {\rm X}$ in which all final states produced by the external current are summed. The ill-posed inverse problem of determining the spectral density from the correlation function is made tractable through the determination of smeared spectral densities in which the desired density is convolved with a set of known smearing kernels of finite width $\epsilon$. The smooth energy dependence of the underlying spectral density enables a controlled $\epsilon \to 0$ extrapolation in the inelastic region, yielding the real-time inclusive rate without reference to individual finite-volume energies or matrix elements. Systematic uncertainties due cutoff effects and residual finite-volume effects are estimated and taken into account in the final error budget. After taking the continuum limit, the results are consistent with the known analytic rate to within the combined statistical and systematic errors. Above energies where 20-particle states contribute, the overall precision is sufficient to discern the four-particle contribution to the spectral density., Comment: 26 pages, 11 figures

We report on the first master-field simulations of QCD with 2+1 dynamical quark flavours using non-perturbatively improved stabilised Wilson fermions. Our simulations are performed at a lattice spacing of 0.094 fm with 96 and 192 points in each direction. On both lattices, the pion and kaon masses are equal to 270 and 450 MeV, respectively, and $m_\pi L$ thus reaches an unprecedented value of 25 on the larger lattice. This setup matches a single point on a chiral trajectory with fixed trace of the quark mass matrix and allows for comparisons to standard large-scale simulations. We present our algorithmic setup and performance measures, and report about our experience in thermalising large master-field lattices with fermions., Comment: 9 pages, 2 figures, 2 tables, talk presented at The 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021, Zoom/Gather@Massachusetts Institute of Technology

We present the current status of our efforts in search of $H$ dibaryon on $N_f$=2+1 CLS ensembles away from the $SU(3)$ flavor symmetric point. Utilizing the distillation framework (also known as LapH) in its exact and stochastic forms, we calculate two-point correlation matrices using large bases of bi-local two-baryon interpolators to reliably determine the low-energy spectra. We report the low lying spectrum on several moving frames for multiple ensembles with different lattice spacing and physical volumes. The status of finite-volume analysis to extract the scattering amplitudes is also discussed., Comment: 9 pages, 5 figures. Talk presented at the 38th International Symposium on Lattice Field Theory (Lattice 2021), 26-30 July, 2021, Zoom/Gather@Massachusetts Institute of Technology

Progress in computing various meson-baryon scattering amplitudes is presented on a single ensemble from the Coordinated Lattice Simulations (CLS) consortium with $m_\pi=200$ MeV and $N_f=2+1$ dynamical fermions. The finite-volume L\"uscher approach is employed to determine the lowest few partial waves from ground- and excited-state energies computed from correlation matrices rotated in a single pivot using a generalized eigenvector solution. This analysis requires evaluating matrices of correlation functions between single- and two-hadron interpolating operators which are projected onto definite spatial momenta and finite-volume irreducible representations. The stochastic LapH method is used to estimate all needed quark propagators. Preliminary results are presented for $I=\frac{1}{2}, \frac{3}{2}$ $N\pi$ amplitudes including the $\Delta(1232)$ resonance and the $I=0$ $S$-wave amplitude with unit strangeness relevant for the $\Lambda(1405)$., Comment: The 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021 Zoom/Gather@Massachusetts Institute of Technology 9 pages, 4 figures

The master-field approach to lattice QCD envisions performing calculations on a small number of large-volume gauge-field configurations. Substantial progress has been made recently in the generation of such fields, and this must be joined with measurement strategies that take advantage of the large volume. In these proceedings, we describe how to compute simple hadronic quantities efficiently and estimate their errors in the master-field approach, i.e. by studying cross-correlations of observables on a single configuration. We discuss the scaling of the uncertainty with the volume and compare extractions based on momentum-projected and position-space two-point functions. The latter show promising results, already at intermediate volumes, but come with additional technical complexities such as a more complicated manifestation of boundary effects, which we also address., Comment: 9 pages, 2 figures, talk presented at The 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021, Zoom/Gather@Massachusetts Institute of Technology

We report on the first application of the stochastic Laplacian Heaviside method for computing multi-particle interactions with lattice QCD to the two-nucleon system. Like the Laplacian Heaviside method, this method allows for the construction of interpolating operators which can be used to construct a positive definite set of two-nucleon correlation functions, unlike nearly all other applications of lattice QCD to two nucleons in the literature. It also allows for a variational analysis in which optimal linear combinations of the interpolating operators are formed that couple predominantly to the eigenstates of the system. Utilizing such methods has become of paramount importance in order to help resolve the discrepancy in the literature on whether two nucleons in either isospin channel form a bound state at pion masses heavier than physical, with the discrepancy persisting even in the $SU(3)$-flavor symmetric point with all quark masses near the physical strange quark mass. This is the first in a series of papers aimed at resolving this discrepancy. In the present work, we employ the stochastic Laplacian Heaviside method without a hexaquark operator in the basis at a lattice spacing of $a\sim0.086$~fm, lattice volume of $L=48a\simeq4.1$~fm and pion mass $m_\pi\simeq714$ MeV. With this setup, the observed spectrum of two-nucleon energy levels strongly disfavors the presence of a bound state in either the deuteron or dineutron channel., Comment: v2: version to be published in Phys. Rev. C.; v1: 13 pages plus figures and appendices

Calculations of the elastic $I=\frac{3}{2}$ nucleon-pion scattering phase shifts on two lattice QCD ensembles with $m_\pi=200\mathrm{MeV}$ and $280\mathrm{MeV}$ are presented. The ensembles both employ $N_\mathrm{f} = 2+1$ Wilson clover fermions. We determine the $\Delta(1232)$ resonance parameters from a finite volume scattering analysis. In one study the single partial wave simplification is employed to compute the $p$-wave amplitude while in the other we treat the partial wave mixing between $s$- and $p$-wave contributions. Fitting our data to a Breit-Wigner resonance model we find $m_\Delta/m_\pi = 7.13(9)$ and $4.75(5)$ on the two ensembles respectively, showing that for a lighter quark mass the resonance mass moves from near the $N\pi$ threshold to near the $N\pi\pi$ threshold, in agreement with experiment., Comment: 7 pages, 1 figure, presented at the 37th International Symposium on Lattice Field Theory

Lattice QCD calculations of resonant meson-meson scattering amplitudes have improved significantly due to algorithmic and computational advances. However, progress in meson-nucleon scattering has been slower due to difficulties in computing the necessary correlation functions, the exponential signal-to-noise problem, and the finite-volume treatment of scattering with fermions. Nonetheless, first benchmark calculations have now been performed. The status of lattice QCD calculations of meson-nucleon scattering amplitudes is reviewed together with comments on future prospects., Comment: 10 pages, 5 figures. Proceeding from a plenary talk at MENU 2019

Lattice QCD allows us to probe the low-lying hadron spectrum in finite-volume using a basis of single- and multi-hadron interpolating operators. Here we examine the effect of including tetraquark operators on the spectrum in the scalar meson sectors containing the $K_0^*(700)$ ($\kappa$) and the $a_0(980)$ in $N_f = 2 + 1$ QCD, with $m_\pi \approx 230$ MeV. Preliminary results of additional finite-volume states found using tetraquark operators are shown, and possible implications of these states are discussed., Comment: 6 pages, 4 figures, proceedings for The 15th International Conference on Meson-Nucleon Physics and the Structure of the Nucleon, MENU-2019

Lattice calculations allow us to probe the low-lying, non-perturbative spectrum of QCD using first principles numerical methods. Here we present the low-lying spectrum in the scalar sector with vacuum quantum numbers including, in fully dynamical QCD for the first time, the mixing between glueball, q-qbar, and meson-meson operators., Comment: 7 pages, 5 figures, Proceedings for The 15th International Conference on Meson-Nucleon Physics and the Structure of the Nucleon, MENU-2019

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., Comment: 22 pages, 2 figures, CERN-TH-2019-035, CP3-Origins-2019-006 DNRF90

The energy spectrum of a system containing a static quark anti-quark pair is computed for a wide range of source separations using lattice QCD with $N_\mathrm{f}=2+1$ dynamical flavours. By employing a variational method with a basis including operators resembling both the gluon string and systems of two separated static mesons, the first three energy levels are determined up to and beyond the distance where it is energetically favourable for the vacuum to screen the static sources through light- or strange-quark pair creation, enabling both these screening phenomena to be observed. The separation dependence of the energy spectrum is reliably parameterised over this saturation region with a simple model which can be used as input for subsequent investigations of quarkonia above threshold and heavy-light and heavy-strange coupled-channel meson scattering., Comment: 12 pages, 3 figures; version accepted for publication in Physics Letters B

The static potential $V(r)$ between a static quark and a static antiquark separated by a distance r is defined as the energy of the ground state of the system. As a consequence of confinement, the energy between the quark-antiquark pair is contained inside a color flux tube, which will break due to pair creation as soon as the energy is high enough. String breaking is manifested as a quantum-mechanical mixing phenomenon between different states containing two infinitely heavy quarks. We investigate this phenomenon with $N_\mathrm{f}=2+1$ flavors of dynamical fermions in the stochastic LapH framework, using an ensemble of gauge configurations generated through the CLS effort. We observe the effect of the third sea-quark flavor, which results in a second mixing-phenomenon., Comment: 7 pages, 2 figures, talk presented at the 36th Annual International Symposium on Lattice Field Theory, July 22-28, 2018, East Lansing, MI, USA

We present a study on ab-initio calculations of three-neutron correlators from Lattice QCD. We extend the method of baryon blocks to systems of three spacially displaced baryons. This allows the measurement of three-neutron $p$-wave correlators with total spin $S=1/2$ and $3/2$. In addition, we use automatic code generation that has high flexibility and allows for easy inclusion of additional channels in the future while optimizing the evaluation of contractions. Our measurements were performed on a newly generated $96\times 48^3$ Clover-Wilson gauge field ensemble with $m_{\pi} \approx 370\,\text{MeV}$. We present preliminary results of our calculations of one pion and nucleon as well as two nucleon ($2N$) and three neutron ($3n$) correlators., Comment: Proceedings from Lattice2018, 7 pages, 4 figures

A field theory involving two interacting scalar fields, previously studied by Rummukainen and Gottlieb, is revisited. Our study is not restricted to the limit of large quartic couplings, and a Symanzik-improved action is used so that continuum dispersion relations work well. The Metropolis method, combined with a local microcanonical updating algorithm, is employed in our Monte Carlo calculations. Isotropic lattices ranging from $16^3 \times 48$ to $53^3 \times 48$ are used, and scattering phase shifts are determined using a L\"uscher analysis with multiple partial waves., Comment: 7 pages, 7 figures, proceedings for The 36th Annual International Symposium on Lattice Field Theory

Elastic $I=1/2$, $s$- and $p$-wave $K\pi$ scattering amplitudes are simultaneously calculated using a L\"uscher style analysis on a single ensemble of dynamical Wilson-clover fermions at $m_\pi \sim 230$ MeV. Partial wave mixing due to the reduced rotational symmetries of the finite volume is included up to $\ell=2$. We also present finite-volume QCD spectra on two large anisotropic lattices ($32^3 \times 256$, $24^3 \times 128$) with $m_\pi \sim 230,\ 390$ MeV respectively. In each symmetry channel, a large basis of one- and two-hadron interpolating operators is employed with all-to-all quark propagation treated using the stochastic LapH method., Comment: 7 pages, proceedings for the 36th International Symposium on Lattice Field Theory (22-28 July 2018); references updated

The elastic $I=1$ $p$-wave $\pi\pi$ scattering amplitude is calculated together with the isovector timelike pion form factor using lattice QCD with $N_{\rm f}=2+1$ dynamical quark flavors. Wilson clover ensembles generated by the Coordinated Lattice Simulations (CLS) initiative are employed at four lattice spacings down to $a = 0.05\,\mathrm{fm}$, several pion masses down to $m_{\pi} = 200\,\mathrm{MeV}$, and spatial volumes of extent $L = 3.1-5.5\,\mathrm{fm}$. The set of measurements on these ensembles, which is publicly available, enables an investigation of systematic errors due to the finite lattice spacing and spatial volume. The $\pi\pi$ scattering amplitude is fit on each ensemble by a Breit-Wigner resonance lineshape, while the form factor is described better by a thrice-subtracted dispersion relation than the Gounaris-Sakurai parametrization., Comment: 32 pages, 12 figures. Typos fixed, small changes to the text, Refs. updated/added, agrees with published version