Predicting positive parity $B_{s}$ mesons from lattice QCD

We determine the spectrum of $B_s$ 1P states using lattice QCD. For the $B_{s1}(5830)$ and $B_{s2}^*(5840)$ mesons, the results are in good agreement with the experimental values. Two further mesons are expected in the quantum channels $J^P=0^+$ and $1^+$ near the $BK$ and $B^{*}K$ thresholds. A combination of quark-antiquark and $B^{(*)}$ meson-Kaon interpolating fields are used to determine the mass of two QCD bound states below the $B^{(*)}K$ threshold, with the assumption that mixing with $B_s^{(*)}eta$ and isospin-violating decays to $B_s^{(*)}pi$ are negligible. We predict a $J^P=0^+$ bound state $B_{s0}$ with mass $m_{B_{s0}}=5.711(13)(19)$ GeV. With further assumptions motivated theoretically by the heavy quark limit, a bound state with $m_{B_{s1}}= 5.750(17)(19)$ GeV is predicted in the $J^P=1^+$ channel. The results from our first principles calculation are compared to previous model-based estimates.

Low lying charmonium states at the physical point

We present results for the mass splittings of low-lying charmonium states from a calculation with Wilson clover valence quarks with the Fermilab interpretation on an asqtad sea. We use five lattice spacings and two values of the light sea quark mass to extrapolate our results to the physical point. Sources of systematic uncertainty in our calculation are discussed and we compare our results for the 1S hyperfine splitting, the 1P-1S splitting and the P-wave spin orbit and tensor splittings to experiment.

Searching for the $X(3872)$ and $Z_c^+(3900)$ on HISQ lattices

We present preliminary simulation results for the I = 0 charmonium state $X(3872)(1^{++})$ and the I = 1 charmonium state $Z_c^+(3900)(1^{+-})$. The study is performed on gauge field configurations with 2+1+1 flavors of highly improved staggered sea quarks (HISQ) with clover (Fermilab interpretation) charm quarks and HISQ light valence quarks. Since the $X(3872)$ lies very close to the open charm $D bar D^*$ threshold, we use a combination of $bar c c$ and $D bar D^* + bar D D^*$ interpolating operators. For the $Z_c^+(3900)$ we use a combination of $J/psi pi$ and $D bar D^* + bar D D^*$ channels. This is the first such study with HISQ sea quarks and light valence quarks. To this end, we describe a variational method for treating staggered quarks that incorporates both oscillating and non-oscillating components.

$DK$ and $D^* K$ scattering near threshold

We study the three $D_s$ quantum channels $J^P = 0^+$, $1^+$ and $2^+$ where experiments have identified the charm-strange states $D^*_{s0} (2317)$, $D_{s1}(2460)$, $D_{s1}(2536)$ near the $DK$ and $D^*K$ thresholds, and $D^*_{s2}(2573)$. We consider correlation functions for sets of $overline q q$ operators and, for $J^P = 0^+$, $1^+$, also the $DK$ and $D^*K$ meson-meson interpolators and determine for these cases values of the elastic scattering amplitude. Constructing the full set of correlators requires propagators which connect any pair of lattice sites. For one ensemble of gauge configurations ($32^3times 64$, $m_piapprox 156$ MeV) a stochastic distillation variant is employed and for another ensemble ($16^3times 32$, $m_piapprox 266$ MeV) we use the full distillation method. Both, $D^*_{s0} (2317)$ and $D_{s1}(2460)$, are found as bound states below threshold, whereas $D_{s1}(2536)$, and $D^*_{s2}(2573)$ are identified as narrow resonances close to the experimental masses.

$D_{s}$ mesons with $DK$ and $D^{*}K$ scattering near threshold

$D_s$ mesons are studied in three quantum channels ($J^P=0^+$, $1^+$ and $2^+$), where experiments have identified the very narrow $D_{s0}^*(2317)$, $D_{s1}(2460)$ and narrow $D_{s1}(2536)$, $D_{s2}^*(2573)$. We explore the effect of nearby $DK$ and $D^*K$ thresholds on the subthreshold states using lattice QCD. Our simulation is done on two very different ensembles of gauge configurations (2 or 2+1 dynamical quarks, Pion mass of 266 or 156 MeV, lattice size $16^3times 32$ or $32^3times 64$). In addition to $bar{q}q$ operators we also include meson-meson interpolators in the correlation functions. This clarifies the identification of the states above and below the scattering thresholds. The ensemble with $m_pi simeq 156~$MeV renders the $D_{s1}(2460)$ as a strong interaction bound state 44(10)MeV below $D^*K$ threshold, which is in agreement with the experiment. The $D_{s0}^*(2317)$ is found 37(17)MeV below $DK$ threshold, close to experiment value of 45MeV. The narrow resonances $D_{s1}(2536)$ and $D_{s2}^*(2573)$ are also found close to the experimental masses.

K pi scattering in moving frames

We extend our study of the $Kpi$ system to moving frames and present an exploratory extraction of the masses and widths for the $K^*$ resonances by simulating $Kpi$ scattering in p-wave with $I=1/2$ on the lattice. Using $Kpi$ systems with non-vanishing total momenta allows the extraction of phase shifts at several values of $Kpi$ relative momenta. A Breit-Wigner fit of the phase renders a $K^*(892)$ resonance mass and $K^*to K pi $ coupling compatible with the experimental numbers. We also determine the $K^*(1410)$ mass assuming the experimental $K^*(1410)$ width. We contrast the resonant $I=1/2$ channel with the repulsive non-resonant $I=3/2$ channel, where the phase is found to be negative and small, in agreement with experiment.

$D_{s0}^*(2317)$ Meson and $D$-Meson-Kaon Scattering from Lattice QCD

The scalar meson $D_{s0}^*(2317)$ is found 37(17)MeV below DK threshold in a lattice simulation of the $J^P=0^+$ channel using, for the first time, both DK as well as $bar sc$ interpolating fields. The simulation is done on $N_f=2+1$ gauge configurations with $m_pisimeq 156 $MeV, and the resulting $M_{D_{s0}^*}-tfrac{1}{4}(M_{D_s}+3M_{D_s^*})=266(16)$ MeV is close to the experimental value 241.5(0.8)MeV. The energy level related to the scalar meson is accompanied by additional discrete levels due to DK scattering states. The levels near threshold lead to the negative DK scattering length $a_0=-1.33(20)$ fm that indicates the presence of a state below threshold.