• We study spatial isovector meson correlators in $N_f=2$ QCD with dynamical domain-wall fermions on $32^3\times 8$ lattices at temperatures $T=220-380$ MeV. We measure the correlators of spin-one ($J=1$) operators including vector, axial-vector, tensor and axial-tensor. Restoration of chiral $U(1)_A$ and $SU(2)_L \times SU(2)_R$ symmetries of QCD implies degeneracies in vector--axial-vector ($SU(2)_L \times SU(2)_R$) and tensor--axial-tensor ($U(1)_A$) pairs, which are indeed observed at temperatures above $T_c$. Moreover, we observe an approximate degeneracy of all $J=1$ correlators with increasing temperature. This approximate degeneracy suggests emergent $SU(2)_{CS}$ and $SU(4)$ symmeries at high temperatures, that mix left- and right-handed quarks.
  • In situations where the low lying eigenmodes of the Dirac operator are suppressed one observed degeneracies of some meson masses. Based on these results a hidden symmetry was conjectured, which is not a symmetry of the Lagrangian but emerges in the quantization process. We show here how the difference between classes of meson propagators is governed by the low modes and shrinks when they disappear.
  • We present results from our recent lattice QCD study of $N\pi$ scattering in the positive-parity nucleon channel, where the puzzling Roper resonance $N^*(1440)$ resides in experiment. Using a variety of hadron operators, that include $qqq$-like, $N\pi$ in $p$-wave and $N\sigma$ in $s$-wave, we systematically extract the excited lattice spectrum in the nucleon channel up to 1.65 GeV. Our lattice results indicate that N$\pi$ scattering in the elastic approximation alone does not describe a low-lying Roper. Coupled channel effects between $N\pi$ and $N\pi\pi$ seem to be crucial to render a low-lying Roper in experiment, reinforcing the notion that this state could be a dynamically generated resonance. After giving a brief motivation for studying the Roper channel and the relevant technical details to this study, we will discuss the results and the conclusions based on our lattice investigation and in comparison with other lattice calculations.
  • We present a lattice QCD study of $N\pi$ scattering in the positive-parity nucleon channel, where the puzzling Roper resonance $N^*(1440)$ resides in experiment. The study is based on the PACS-CS ensemble of gauge configurations with $N_f=2+1$ Wilson-clover dynamical fermions, $m_\pi \simeq 156~$MeV and $L\simeq 2.9~$fm. In addition to a number of $qqq$ interpolating fields, we implement operators for $N\pi$ in $p$-wave and $N\sigma$ in $s$-wave. In the center-of-momentum frame we find three eigenstates below 1.65 GeV. They are dominated by $N(0)$, $N(0)\pi(0)\pi(0)$ (mixed with $N(0)\sigma(0)$) and $N(p)\pi(-p)$ with $p\simeq 2\pi/L$, where momenta are given in parentheses. This is the first simulation where the expected multi-hadron states are found in this channel. The experimental $N\pi$ phase-shift would -- in the approximation of purely elastic $N\pi$ scattering -- imply an additional eigenstate near the Roper mass $m_R\simeq 1.43~$GeV for our lattice size. We do not observe any such additional eigenstate, which indicates that $N\pi$ elastic scattering alone does not render a low-lying Roper. Coupling with other channels, most notably with $N\pi\pi$, seems to be important for generating the Roper resonance, reinforcing the notion that this state could be a dynamically generated resonance. Our results are in line with most of previous lattice studies based just on $qqq$ interpolators, that did not find a Roper eigenstate below $1.65~$GeV. The study of the coupled-channel scattering including a three-particle decay $N\pi\pi$ remains a challenge.
  • Working with a pion mass $m_\pi \approx 150$ MeV, we study $\pi\pi$ and $K\pi$ scattering using two flavours of non-perturbatively improved Wilson fermions at a lattice spacing $a\approx 0.071$ fm. Employing two lattice volumes with linear spatial extents of $N_s=48$ and $N_s=64$ points and moving frames, we extract the phase shifts for p-wave $\pi\pi$ and $K\pi$ scattering near the $\rho$ and $K^*$ resonances.Comparing our results to those of previous lattice studies, that used pion masses ranging from about 200 MeV up to 470 MeV, we find that the coupling $g_{\rho\pi\pi}$ appears to be remarkably constant as a function of $m_{\pi}$.
  • Lattice calculations for hadrons are now entering the domain of resonances and scattering, necessitating a better understanding of the observed discrete energy spectrum. This is a reviewing survey about recent lattice QCD results, with some emphasis on spectrum and scattering.
  • We discuss a recent lattice study of charmonium-like mesons with $J^{PC}=1^{++}$ and three quark contents $\bar cc\bar du$, $\bar cc(\bar uu + \bar dd)$ and $\bar cc\bar ss$, where the latter two can mix with $\bar cc$. In this quantum channel, the long known exotic candidate, X(3872), resides. This simulation employs $N_f=2$, $m_\pi=266~$MeV and a large basis of $\bar cc$, two-meson and diquark-antidiquark interpolating fields, with diquarks in both anti-triplet and sextet color representations. It aims at the possible signatures of four-quark exotic states. Along the way, we discuss the relations between the diquark-antidiquark operators and the two-meson operators via the Fierz transformations.
  • 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.
  • We perform a lattice study of charmonium-like mesons with $J^{PC}=1^{++}$ and three quark contents $\bar cc \bar du$, $\bar cc(\bar uu+\bar dd)$ and $\bar cc \bar ss$, where the later two can mix with $\bar cc$. This simulation with $N_f=2$ and $m_\pi=266$ MeV aims at the possible signatures of four-quark exotic states. We utilize a large basis of $\bar cc$, two-meson and diquark-antidiquark interpolating fields, with diquarks in both anti-triplet and sextet color representations. A lattice candidate for X(3872) with I=0 is observed very close to the experimental state only if both $\bar cc$ and $D\bar D^*$ interpolators are included; the candidate is not found if diquark-antidiquark and $D\bar D^*$ are used in the absence of $\bar cc$. No candidate for neutral or charged X(3872), or any other exotic candidates are found in the I=1 channel. We also do not find signatures of exotic $\bar cc\bar ss$ candidates below 4.3 GeV, such as Y(4140). Possible physics and methodology related reasons for that are discussed. Along the way, we present the diquark-antidiquark operators as linear combinations of the two-meson operators via the Fierz transformations.
  • Recently experimentalists have discovered several charged charmonium-like hadrons $Z_c^+$ with unconventional quark content $\bar cc\bar d u$. We perform a search for $Z_c^+$ with mass below $4.2~$GeV in the channel $I^G(J^{PC})=1^+(1^{+-})$ using lattice QCD. The major challenge is presented by the two-meson states $J/\psi\, \pi$, $\psi_{2S}\pi$, $\psi_{1D}\pi$, $D\bar D^*$, $D^*\bar D^*$, $\eta_c\rho$ that are inevitably present in this channel. The spectrum of eigenstates is extracted using a number of meson-meson and diquark-antidiquark interpolating fields. For our pion mass of 266~MeV we find all the expected two-meson states but no additional candidate for $Z_c^+$ below $4.2~$GeV. Possible reasons for not seeing an additional eigenstate related to $Z_c^+$ are discussed. We also illustrate how a simulation incorporating interpolators with a structure resembling low-lying two-mesons states seems to render a $Z_c^+$ candidate, which is however not robust after further two-meson states around $4.2~$GeV are implemented.
  • 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^3\times 64$, $m_\pi\approx 156$ MeV) a stochastic distillation variant is employed and for another ensemble ($16^3\times 32$, $m_\pi\approx 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.
  • We employ a variational basis with a number of $\bar{q}q$ and $\pi\pi$ lattice interpolating fields with quantum numbers of the $\rho$ resonance to extract the discrete energy spectrum in a finite volume. In the elastic region, this spectrum is related to the phase shift of the continuum scattering amplitude by L\"uscher's formula and the relation allows the extraction of resonance parameters from the spectrum calculation. The simulations are performed at three different total momenta of the coupled $\bar q q-\pi\pi$ system, which allows us to extract the p-wave scattering phase at five values of pion relative momenta near the resonance region. The effective range formula describes the phase-shift dependence nicely and we extract the resonance mass $m_\rho=792(7)(8)$ MeV and the coupling $g_{\rho\pi\pi}=5.13(20)$ at our $m_\pi\simeq 266 $MeV. The coupling $g_{\rho\pi\pi}$ is directly related to the width of the $\rho$ meson and our value is close to the value derived from the experimental width. The simulations are performed using dynamical gauge configurations with two mass-degenerate flavors of tree-level improved clover-Wilson fermions. Correlation functions are calculated using the recently proposed distillation method with Laplacian Heaviside (LapH) smearing of quarks, which enables flexible calculations, in many cases with unprecedented accuracy.
  • $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^3\times 32$ or $32^3\times 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.
  • 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_\pi\simeq 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.
  • A nonvanishing spectral density of the low-lying eigenmodes of the Dirac operator naturally is a signal for dynamical chiral symmetry breaking (D$\chi$SB) via the Banks-Casher relation. The low-lying eigenmodes alone saturate the pseudoscalar channel and the corresponding propagator successfully reproduces the pion mass. In this paper we investigate the effects on the mass generation of hadrons other than pions. The evolution of these masses upon inclusion of an increasing number of the low-lying eigenmodes is confronted with the hadron mass spectrum upon removal of such eigenmodes.
  • We extend our study of the $K\pi$ system to moving frames and present an exploratory extraction of the masses and widths for the $K^*$ resonances by simulating $K\pi$ scattering in p-wave with $I=1/2$ on the lattice. Using $K\pi$ systems with non-vanishing total momenta allows the extraction of phase shifts at several values of $K\pi$ 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.
  • Including the meson-baryon (5 quark) intermediate states in a lattice simulation is challenging. However, it is important in order to obtain the correct energy eigenstates and to relate them to scattering phase shifts. Recent results for the negative parity nucleon channel and the problem of baryonic resonances in lattice calculations are discussed.
  • In previous studies we have shown that hadrons, except for a pion, survive the removal of the lowest lying Dirac eigenmodes from the valence quark propagators. The low-modes are tied to the dynamical breaking of chiral symmetry and we found chiral symmetry to be restored by means of matching masses of chiral partners, like, e.g., the vector and axial vector currents. Here we investigate the influence of removing the lowest part of the Dirac spectrum on the locality of the Dirac operator. Moreover, we analyze the influence of low-mode truncation on the quark momenta and thereupon on the hadron spectrum and, finally, introduce a reweighting scheme to extend the truncation to the sea quark sector.
  • I discuss how masses and widths of hadron resonances are extracted from lattice QCD. Recent lattice results on the light, strange and charm meson resonances are reviewed. Their properties are revealed by simulating the corresponding scattering channels pi-pi, K-pi and D-pi on the lattice and extracting the scattering phase shifts. In particular we address the resonances rho, D0*(2400), D1(2430), K*, kappa and K0*(1430).
  • Results for the excited meson and baryon spectrum with two flavors of Chirally Improved sea quarks are presented. We simulate several ensembles with pion masses ranging from 250 to 600 MeV and extrapolate to the physical pion mass. Strange quarks are treated within the partially quenched approximation. Using the variational method, we investigate the content of the states. Among others, we discuss the flavor singlet/octet content of Lambda states. In general, our results compare well with experiment, in particular we get very good agreement with the Lambda(1405) and confirm its flavor singlet nature.
  • We study the coupled pion-nucleon system (negative parity, isospin 1/2) based on a lattice QCD simulation for nf=2 mass degenerate light quarks. Both, standard 3-quarks baryon operators as well as meson-baryon (4+1)-quark operators are included. This is an exploratory study for just one lattice size and lattice spacing and at a pion mass of 266 MeV. Using the distillation method and variational analysis we determine energy levels of the lowest eigenstates. Comparison with the results of simple 3-quark correlation studies exhibits drastic differences and a new level appears. A clearer picture of the negative parity nucleon spectrum emerges. For the parameters of the simulation we may assume elastic s-wave scattering and can derive values of the phase shift.
  • In a lattice QCD calculation with two light dynamical Chirally Improved (CI) quarks we determine ground state and some excited state masses in all four Lambda baryon channels 1/2^\pm and 3/2^\pm. We perform an infinite volume extrapolation and confirm the widely discussed Lambda(1405). We also analyze the amount of octet-singlet mixing, which is helpful in comparing states with the quark model.
  • We simulate K pi scattering in s-wave and p-wave for both isospins I=1/2, 3/2 using quark-antiquark and meson-meson interpolating fields. We extract the elastic phase shifts delta at several values of the K-pi relative momenta. The resulting phases exhibit qualitative agreement with the experimental phases in all four channels. We express the s-wave phase shifts near threshold in terms of the scattering length and the effective range. Our K pi system has zero total momentum and is simulated on a single ensemble with two dynamical quarks, so results apply for mpi=266 MeV and mK=552 MeV in our simulation. The backtracking contractions in both I=1/2 channels are handled by the use of Laplacian-Heavyside smeared quarks within the distillation method. Elastic phases are extracted from the energy levels using Luscher's relations. In all four channels we observe the expected K(n)pi(-n) scattering states, which are shifted due to the interaction. In both attractive I=1/2 channels we observe additional states that are related to resonances; we attribute them to K_0^*(1430) in s-wave and K*(892), K*(1410) and K*(1680) in p-wave.
  • Chiral symmetry breaking in Quantum Chromodynamics is associated with the low lying spectral modes of the Dirac operator according to the Banks-Casher relation. Here we study how removal of a variable number of low lying modes from the valence quark sector affects the masses of the ground states and first excited states of baryons and mesons in two flavor lattice QCD.
  • We study hadron correlators upon artificial restoration of the spontaneously broken chiral symmetry. In a dynamical lattice simulation we remove the lowest lying eigenmodes of the Dirac operator from the valence quark propagators and study evolution of the hadron masses obtained. All mesons and baryons in our study, except for a pion, survive unbreaking the chiral symmetry and their exponential decay signals become essentially better. From the analysis of the observed spectroscopic patterns we conclude that confinement still persists while the chiral symmetry is restored. All hadrons fall into different chiral multiplets. The broken U(1)_A symmetry does not get restored upon unbreaking the chiral symmetry. We also observe signals of some higher symmetry that includes chiral symmetry as a subgroup. Finally, from comparison of the \Delta - N splitting before and after unbreaking of the chiral symmetry we conclude that both the color-magnetic and the flavor-spin quark-quark interactions are of equal importance.