A new implementation of estimating the two-to-two $K$-matrix from finite-volume energies based on the Luescher formalism is described. The method includes higher partial waves and multiple decay channels, and the fitting procedure properly includes a
ll covariances and statistical uncertainties. The method is also simpler than previously used procedures. Formulas and software for handling total spins up to $S=2$ and orbital angular momenta up to $L=6$ are presented.
An implementation of estimating the two-to-two $K$-matrix from finite-volume energies based on the Luscher formalism and involving a Hermitian matrix known as the box matrix is described. The method includes higher partial waves and multiple decay ch
annels. Two fitting procedures for estimating the $K$-matrix parameters, which properly incorporate all statistical covariances, are discussed. Formulas and software for handling total spins up to $S=2$ and orbital angular momenta up to $L=6$ are obtained for total momenta in several directions. First tests involving $rho$-meson decay to two pions include the $L=3$ and $L=5$ partial waves, and the contributions from these higher waves are found to be negligible in the elastic energy range.
We non-perturbatively determine the renormalization factor of the axial vector current in lattice QCD with $N_f=3$ flavors of Wilson-clover fermions and the tree-level Symanzik-improved gauge action. The (by now standard) renormalization condition is
derived from the massive axial Ward identity and it is imposed among Schr{o}dinger functional states with large overlap on the lowest lying hadronic state in the pseudoscalar channel, in order to reduce kinematically enhanced cutoff effects. We explore a range of couplings relevant for simulations at lattice spacings of $approx 0.09$ fm and below. An interpolation formula for $Z_A(g_0^2)$, smoothly connecting the non-perturbative values to the 1-loop expression, is provided together with our final results.
The $I=1$ $p$-wave and $I=2$ $s$-wave elastic $pi$-$pi$ scattering amplitudes are calculated from a first-principles lattice QCD simulation using a single ensemble of gauge field configurations with $N_{mathrm{f}} = 2+1$ dynamical flavors of anisotro
pic clover-improved Wilson fermions. This ensemble has a large spatial volume $V=(3.7mathrm{fm})^3$, pion mass $m_{pi} = 230mathrm{MeV}$, and spatial lattice spacing $a_s = 0.11mathrm{fm}$. Calculation of the necessary temporal correlation matrices is efficiently performed using the stochastic LapH method, while the large volume enables an improved energy resolution compared to previous work. For this single ensemble we obtain $m_{rho}/m_{pi} = 3.350(24)$, $g_{rhopipi} = 5.99(26)$, and a clear signal for the $I=2$ $s$-wave. The success of the stochastic LapH method in this proof-of-principle large-volume calculation paves the way for quantitative study of the lattice spacing effects and quark mass dependence of scattering amplitudes using state-of-the-art ensembles.
We report on progress applying the stochastic LapH method to estimate all-to-all propagators required in correlation functions of multi-hadron operators relevant for pion-pion scattering. Large-volume results for $I=2$ and $I=1$ pion-pion scattering
phase shifts with good statistical precision are obtained from an $N_{rm f} = 2+1$ anisotropic Wilson clover ensemble with $m_{pi} = 240mathrm{MeV}$. We also present a preliminary determination of the $I=1$ pion-pion scattering phase shift and timelike pion form factor on an isotropic $N_{rm f}=2+1$ flavour ensemble generated by the Coordinated Lattice Simulation (CLS) community effort.
Recent results in computing excited-state energies and meson-meson scattering phase shifts in lattice QCD are presented. A stochastic method of treating the low-lying modes of quark propagation that exploits Laplacian Heaviside quark-field smearing m
akes such studies possible now on large 32^3 x 256 and 48^3 x 128 lattices at near physical pion masses. Levels are identified using a variety of probe interpolating operators, which include both single-hadron and a large number of two-hadron operators.
We report on a non-perturbative computation of the renormalization factor Z_A of the axial vector current in three-flavour O(a) improved lattice QCD with Wilson quarks and tree-level Symanzik improved gauge action and also recall our recent determina
tion of the improvement coefficient c_A. Our normalization and improvement conditions are formulated at constant physics in a Schrodinger functional setup. The normalization condition exploits the full, massive axial Ward identity to reduce finite quark mass effects in the evaluation of Z_A and correlators with boundary wave functions to suppress excited state contributions in the pseudoscalar channel.
The coefficient c_A required for O(a) improvement of the axial current in lattice QCD with N_f=3 flavors of Wilson fermions and the tree-level Symanzik-improved gauge action is determined non-perturbatively. The standard improvement condition using S
chroedinger functional boundary conditions is employed at constant physics for a range of couplings relevant for simulations at lattice spacings of ~ 0.09 fm and below. We define the improvement condition projected onto the zero topological charge sector of the theory, in order to avoid the problem of possibly insufficient tunneling between topological sectors in our simulations at the smallest bare coupling. An interpolation formula for c_A(g_0^2) is provided together with our final results.
Progress in calculating scattering phase shifts on $N_f=2+1$ anisotropic clover Wilson lattices is described. The stochastic LapH method facilitates computations in large volumes and for light pion masses. Results for pion masses down to 240 MeV, keeping $m_pi L > 4$, are presented.
Progress in computing the spectrum of excited baryons and mesons in lattice QCD is described. Results in the zero-momentum bosonic I=1/2, S=1, T1u symmetry sector of QCD using a correlation matrix of 58 operators are presented. All needed Wick contra
ctions are efficiently evaluated using a stochastic method of treating the low-lying modes of quark propagation that exploits Laplacian Heaviside quark-field smearing. Level identification using probe operators is discussed.
