
We study spatial isovector meson correlators in $N_f=2$ QCD with dynamical
domainwall fermions on $32^3\times 8$ lattices at temperatures $T=220380$
MeV. We measure the correlators of spinone ($J=1$) operators including vector,
axialvector, tensor and axialtensor. Restoration of chiral $U(1)_A$ and
$SU(2)_L \times SU(2)_R$ symmetries of QCD implies degeneracies in
vectoraxialvector ($SU(2)_L \times SU(2)_R$) and tensoraxialtensor
($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
righthanded 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 positiveparity 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 lowlying Roper. Coupled channel
effects between $N\pi$ and $N\pi\pi$ seem to be crucial to render a lowlying
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 positiveparity
nucleon channel, where the puzzling Roper resonance $N^*(1440)$ resides in
experiment. The study is based on the PACSCS ensemble of gauge configurations
with $N_f=2+1$ Wilsonclover 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
centerofmomentum 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 multihadron states are found
in this channel. The experimental $N\pi$ phaseshift 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 lowlying 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 coupledchannel scattering including a
threeparticle 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 nonperturbatively 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 pwave $\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 charmoniumlike 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$,
twomeson and diquarkantidiquark interpolating fields, with diquarks in both
antitriplet and sextet color representations. It aims at the possible
signatures of fourquark exotic states. Along the way, we discuss the relations
between the diquarkantidiquark operators and the twomeson 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 quarkantiquark and $B^{(*)}$ mesonKaon 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
isospinviolating 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 modelbased estimates.

We perform a lattice study of charmoniumlike 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 fourquark
exotic states. We utilize a large basis of $\bar cc$, twomeson and
diquarkantidiquark interpolating fields, with diquarks in both antitriplet
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
diquarkantidiquark 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
diquarkantidiquark operators as linear combinations of the twomeson operators
via the Fierz transformations.

Recently experimentalists have discovered several charged charmoniumlike
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 twomeson 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 mesonmeson and
diquarkantidiquark interpolating fields. For our pion mass of 266~MeV we find
all the expected twomeson 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 lowlying twomesons states seems to
render a $Z_c^+$ candidate, which is however not robust after further twomeson
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 charmstrange 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$ mesonmeson
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 pwave scattering phase at five values
of pion relative momenta near the resonance region. The effective range formula
describes the phaseshift 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 massdegenerate flavors of treelevel improved
cloverWilson 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 mesonmeson 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 lowlying eigenmodes of the Dirac
operator naturally is a signal for dynamical chiral symmetry breaking
(D$\chi$SB) via the BanksCasher relation. The lowlying 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 lowlying 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 pwave with $I=1/2$ on the lattice. Using
$K\pi$ systems with nonvanishing total momenta allows the extraction of phase
shifts at several values of $K\pi$ relative momenta. A BreitWigner 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 nonresonant $I=3/2$ channel, where the
phase is found to be negative and small, in agreement with experiment.

Including the mesonbaryon (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 lowmodes 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 lowmode 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 pipi, Kpi and Dpi 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 pionnucleon system (negative parity, isospin 1/2) based
on a lattice QCD simulation for nf=2 mass degenerate light quarks. Both,
standard 3quarks baryon operators as well as mesonbaryon (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 3quark 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 swave 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 octetsinglet mixing, which is helpful in comparing states with
the quark model.

We simulate K pi scattering in swave and pwave for both isospins I=1/2, 3/2
using quarkantiquark and mesonmeson interpolating fields. We extract the
elastic phase shifts delta at several values of the Kpi relative momenta. The
resulting phases exhibit qualitative agreement with the experimental phases in
all four channels. We express the swave 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
LaplacianHeavyside 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 swave and K*(892), K*(1410) and K*(1680) in pwave.

Chiral symmetry breaking in Quantum Chromodynamics is associated with the low
lying spectral modes of the Dirac operator according to the BanksCasher
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 colormagnetic and the flavorspin quarkquark interactions are
of equal importance.