
We present results for several light hadronic quantities ($f_\pi$, $f_K$,
$B_K$, $m_{ud}$, $m_s$, $t_0^{1/2}$, $w_0$) obtained from simulations of 2+1
flavor domain wall lattice QCD with large physical volumes and nearlyphysical
pion masses at two lattice spacings. We perform a short, O(3)%, extrapolation
in pion mass to the physical values by combining our new data in a simultaneous
chiral/continuum `global fit' with a number of other ensembles with heavier
pion masses. We use the physical values of $m_\pi$, $m_K$ and $m_\Omega$ to
determine the two quark masses and the scale  all other quantities are outputs
from our simulations. We obtain results with subpercent statistical errors and
negligible chiral and finitevolume systematics for these light hadronic
quantities, including: $f_\pi$ = 130.2(9) MeV; $f_K$ = 155.5(8) MeV; the
average up/down quark mass and strange quark mass in the $\bar {\rm MS}$ scheme
at 3 GeV, 2.997(49) and 81.64(1.17) MeV respectively; and the neutral kaon
mixing parameter, $B_K$, in the RGI scheme, 0.750(15) and the $\bar{\rm MS}$
scheme at 3 GeV, 0.530(11).

We determine the neutral kaon mixing matrix element $B_K$ in the continuum
limit with 2+1 flavors of domain wall fermions, using the Iwasaki gauge action
at two different lattice spacings. These lattice fermions have near exact
chiral symmetry and therefore avoid artificial lattice operator mixing.
We introduce a significant improvement to the conventional NPR method in
which the bare matrix elements are renormalized nonperturbatively in the
RIMOM scheme and are then converted into the MSbar scheme using continuum
perturbation theory. In addition to RIMOM, we introduce and implement four
nonexceptional intermediate momentum schemes that suppress infrared
nonperturbative uncertainties in the renormalization procedure. We compute the
conversion factors relating the matrix elements in this family of RISMOM
schemes and MSbar at oneloop order. Comparison of the results obtained using
these different intermediate schemes allows for a more reliable estimate of the
unknown higherorder contributions and hence for a correspondingly more robust
estimate of the systematic error. We also apply a recently proposed approach in
which twisted boundary conditions are used to control the Symanzik expansion
for offshell vertex functions leading to a better control of the
renormalization in the continuum limit.
We control chiral extrapolation errors by considering both the NLO SU(2)
chiral effective theory, and an analytic mass expansion. We obtain
$B_K^{\msbar}(3 GeV) = 0.529(5)_{stat}(15)_\chi(2)_{FV}(11)_{NPR}$. This
corresponds to $\hat{B}_K = 0.749(7)_{stat}(21)_\chi(3)_{FV}(15)_{NPR}$. Adding
all sources of error in quadrature we obtain $\hat{B}_K =
0.749(27)_{combined}$, with an overall combined error of 3.6%.

We present the first results for the Kl3 form factor from simulations with
2+1 flavours of dynamical domain wall quarks. Combining our result, namely
f_+(0)=0.964(5), with the latest experimental results for Kl3 decays leads to
V_{us}=0.2249(14), reducing the uncertaintity in this important parameter.
For the O(p^6) term in the chiral expansion we obtain \Delta f=0.013(5).

We have simulated QCD using 2+1 flavors of domain wall quarks on a $(2.74
{\rm fm})^3$ volume with an inverse lattice scale of $a^{1} = 1.729(28)$ GeV.
The up and down (light) quarks are degenerate in our calculations and we have
used four values for the ratio of light quark masses to the strange (heavy)
quark mass in our simulations: 0.217, 0.350, 0.617 and 0.884. We have measured
pseudoscalar meson masses and decay constants, the kaon bag parameter $B_K$ and
vector meson couplings. We have used SU(2) chiral perturbation theory, which
assumes only the up and down quark masses are small, and SU(3) chiral
perturbation theory to extrapolate to the physical values for the light quark
masses. While nexttoleading order formulae from both approaches fit our data
for light quarks, we find the higher order corrections for SU(3) very large,
making such fits unreliable. We also find that SU(3) does not fit our data when
the quark masses are near the physical strange quark mass. Thus, we rely on
SU(2) chiral perturbation theory for accurate results. We use the masses of the
$\Omega$ baryon, and the $\pi$ and $K$ mesons to set the lattice scale and
determine the quark masses. We then find $f_\pi = 124.1(3.6)_{\rm
stat}(6.9)_{\rm syst} {\rm MeV}$, $f_K = 149.6(3.6)_{\rm stat}(6.3)_{\rm syst}
{\rm MeV}$ and $f_K/f_\pi = 1.205(0.018)_{\rm stat}(0.062)_{\rm syst}$. Using
nonperturbative renormalization to relate lattice regularized quark masses to
RIMOM masses, and perturbation theory to relate these to $\bar{\rm MS}$ we
find $ m_{ud}^{\bar{\rm MS}}(2 {\rm GeV}) = 3.72(0.16)_{\rm stat}(0.33)_{\rm
ren}(0.18)_{\rm syst} {\rm MeV}$ and $m_{s}^{\bar{\rm MS}}(2 {\rm GeV}) =
107.3(4.4)_{\rm stat}(9.7)_{\rm ren}(4.9)_{\rm syst} {\rm MeV}$.

We present results for light meson masses and pseudoscalar decay constants
from the first of a series of lattice calculations with 2+1 dynamical flavors
of domain wall fermions and the Iwasaki gauge action. The work reported here
was done at a fixed lattice spacing of about 0.12 fm on a 16^3\times32 lattice,
which amounts to a spatial volume of (2 fm)^3 in physical units. The number of
sites in the fifth dimension is 16, which gives m_{res} = 0.00308(4) in these
simulations. Three values of input light sea quark masses, m_l^{sea} \approx
0.85 m_s, 0.59 m_s and 0.33 m_s were used to allow for extrapolations to the
physical light quark limit, whilst the heavier sea quark mass was fixed to
approximately the physical strange quark mass m_s. The exact rational hybrid
Monte Carlo algorithm was used to evaluate the fractional powers of the fermion
determinants in the ensemble generation. We have found that f_\pi = 127(4) MeV,
f_K = 157(5) MeV and f_K/f_\pi = 1.24(2), where the errors are statistical
only, which are in good agreement with the experimental values.

We present results for the light meson masses, the bare strange quark mass
and preliminary nonperturbative renormalisation of B_K in 2+1 flavour domain
wall QCD. The ensembles used were generated with the Iwasaki gauge action and
have a volume of 16^3 x 32 with a fifth dimension size of 16 and an inverse
lattice spacing of 1.6 GeV. These ensembles have u and d masses as low as one
quarter of the strange quark mass. All data were generated jointly by the UKQCD
and RBC collaborations on QCDOC machines.

We present results for light meson masses and psedoscalar meson decay
constants in 2+1 flavour domain wall QCD with the DBW2 and Iwasaki gauge
actions, using lattices with linear sizes in the range 1.6 to 2.2fm and $u$ and
$d$ quark masses as low as one quarter of the strange quark mass. All data were
generated on the QCDOC machines at the University of Edinburgh and Brookhaven
National Laboratory. Despite large residual masses and a limited number of sea
quark mass values with which to perform chiral extrapolations, our results
agree with experiment and scale within errors.

We present results for some of the light baryon masses and their excited
states in 2+1 flavour domain wall QCD. We considered several lattice spacings,
with the DBW2 and Iwasaki gauge actions and different sea quark masses on a
volume of $16^3\times32$ and a fifth dimension of size 8. All data were
generated on the QCDOC machines. Despite large residual massses and a limited
number of sea quark mass values with which to perform chiral extrapolations,
our results are in reasonable agreement with experiment and scale within
errors. Finite size effects on most ensembles appear to be small.

We discuss some of the implications of simulating QCD when the action used
for the sea quarks is different from that used for the valence quarks. We
present exploratory results for the hadron mass spectrum and pseudoscalar meson
decay constants using improved staggered sea quarks and HYPsmeared overlap
valence quarks. We propose a method for matching the valence quark mass to the
sea quark mass and demonstrate it on UKQCD clover data in the simpler case
where the sea and valence actions are the same.

We present exploratory results for the hadron mass spectrum and pseudoscalar
meson decay constants using mixed actions. We use improved staggered sea quarks
and HYPsmeared overlap valence quarks. We obtain good signals on 10
configurations at one lattice spacing and two different sets of sea quark
masses.

We report on progress in setting up the International Lattice Data Grid. We
describe the aims and objectives of the ILDG, what has been achieved during its
first year of activity and invite feedback from the community.

The spectrum of orbitally excited $D_s$ mesons is computed in the continuum
limit of quenched lattice QCD. The results are consistent with the
interpretation that the narrow resonance in the $D_s \pi^0$ channel discovered
by the BABAR Collaboration is a $J^P=0^+$ $c\bar{s}$ meson. Furthermore, within
statistical errors, the $1^+1^$ and the $0^+0^$ mass splittings are equal,
in agreement with the chiral multiplet structure predicted by heavy hadron
chiral effective theory. On our coarsest lattice we present results from the
first study of orbitally excited $D_s$ mesons with two flavors of dynamical
quarks, with mass slightly larger than the strange quark mass. These results
are consistent with the quenched data.

We present results for the static interquark potential, lightest glueballs,
light hadron spectrum and topological susceptibility using a nonperturbatively
improved action on a $16^3\times 32$ lattice at a set of values of the bare
gauge coupling and bare dynamical quark mass chosen to keep the lattice size
fixed in physical units ($\sim 1.7$ fm). By comparing these measurements with a
matched quenched ensemble, we study the effects due to two degenerate flavours
of dynamical quarks. With the greater control over residual lattice spacing
effects which these methods afford, we find some evidence of charge screening
and some minor effects on the light hadron spectrum over the range of quark
masses studied ($M_{PS}/M_{V}\ge0.58$). More substantial differences between
quenched and unquenched simulations are observed in measurements of topological
quantities.

We present a comprehensive study of the masses of pseudoscalar and vector
mesons, as well as octet and decuplet baryons computed in O(a) improved
quenched lattice QCD. Results have been obtained using the nonperturbative
definition of the improvement coefficient c_sw, and also its estimate in
tadpole improved perturbation theory. We investigate effects of improvement on
the incidence of exceptional configurations, mass splittings and the parameter
J. By combining the results obtained using nonperturbative and tadpole
improvement in a simultaneous continuum extrapolation we can compare our
spectral data to experiment. We confirm earlier findings by the CPPACS
Collaboration that the quenched light hadron spectrum agrees with experiment at
the 10% level.

Recent results for the spectrum of light hadrons provide clear evidence for
the failure of quenched QCD and encouraging signs that simulations with
dynamical sea quarks rectify some of the discrepancies, although string
breaking has not yet been observed. The use of perturbation theory to match
lattice quark masses to continuum schemes remains questionable, but
nonperturbative methods are poised to remove this uncertainty. The inclusion
of dynamical sea quarks substantially reduces estimates of the light quark
masses. New results for the lightest glueball and the lightest exotic hybrid
state provide useful input to phenomenology, but still have limited or no
treatment of mixing. The $O(a)$improved Wilson quark action is
wellestablished in quenched QCD for $\beta\geq 5.7$, with most parameters
obtainable nonperturbatively, in which range scaling violations are small.
Progress has also been made with highorder improvement schemes for both Wilson
and staggered quarks.

We present the first results for the static quark potential and the light
hadron spectrum using dynamical fermions at $\beta=5.2$ using an O(a) improved
Wilson fermion action together with the standard Wilson plaquette action for
the gauge part. Sea quark masses were chosen such that the pseudoscalarvector
mass ratio, m_PS/m_V$, varies from 0.86 to 0.67. Finitesize effects are
studied by using three different volumes, 8^3\cdot 24, 12^3\cdot 24 and
16^3\cdot 24. Comparing our results to previous ones obtained using the
quenched approximation, we find evidence for sea quark effects in quantities
like the static quark potential and the vectorpseudoscalar hyperfine
splitting.

We present the results of the first lattice study of semileptonic decays of
baryons containing a bquark. Predictions for the decay distributions are given
and the IsgurWise functions for heavy baryons are computed, for values of the
velocity transfer up to about omega=1.2. The computations are performed on a
24^3 x 48 lattice at beta=6.2 using the SheikholeslamiWohlert action in the
quenched approximation.

We compute the light hadron mass spectrum at beta=5.7 using the O(a) improved
SheikholeslamiWohlert (SW) fermion action with two choices of the clover
coefficient: the classical value, c=1, and a meanfield or tadpoleimproved
estimate c=1.57. We compare our results with those of the GF11 Collaboration
who use the Wilson fermion action (c=0). We find that changing c from zero to 1
and 1.57 leads to significant differences in the masses of the chirally
extrapolated and strange pseudoscalar and vector mesons, the nucleon, the
Delta, and also in the Edinburgh plot. A number of other quantities, for
example m_V^2m_PS^2, J, am_K/am_\rho and am_{K^*}/am_\rho do not appear to
change significantly. We also investigate the effect of changing the lattice
volume from approximately (2 fm)^3 to (2.6 fm)^3. We find that the meson masses
are consistent to within one standard deviation and baryon masses are
consistent to within two standard deviations.

After a brief introduction to the HeavyQuark Effective Theory (HQET), I
review the extraction of the IsgurWise function from lattice QCD calculations
of the matrix elements for semileptonic decays of heavylight pseudoscalar
mesons both into pseudoscalar and into vector mesons. This work is beginning to
test the heavyquark spinflavour symmetries around the charm mass and to
indicate the size of $O(1/m_c)$ corrections. An alternative approach to put the
HQET on the lattice offers the prospect of computing the IsgurWise function
directly.

We compute the hyperfine splitting $m_{J/\psi}m_{\eta_c}$ on the lattice,
using both the Wilson and $O(a)$improved (clover) actions for quenched quarks.
The computations are performed on a $24^3\times48$ lattice at $\beta = 6.2$,
using the same set of 18 gluon configurations for both fermion actions. We find
that the splitting is 1.83\err{13}{15} times larger with the clover action than
with the Wilson action, demonstrating the sensitivity of the spinsplitting to
the magnetic moment term which is present in the clover action. However, even
with the clover action the result is less than half of the physical
masssplitting. We also compute the decay constants $f_{\eta_c}$ and
$f^{1}_{J/\psi}$, both of which are considerably larger when computed using
the clover action than with the Wilson action. For example for the ratio
$f^{1}_{J/\psi}/f^{1}_{\rho}$ we find 0.32\err{1}{2} with the Wilson action
and $0.48\pm 3$ with the clover action (the physical value is 0.44(2)).

We present a technique which permits the calculation of twopoint functions
of operators containing one heavy quark and an arbitrary number of light quarks
as analytic functions of the heavyquark mass. It is based on the standard
Jacobi linear solver used for the calculation of quark propagators. Results for
the heavylight pseudoscalar and vector meson masses are obtained on 16^3x48
lattices at beta = 6.2 using the Wilson fermion action, and agree with
published data. The incorporation of smeared operators and $O(a)$improved
actions presents no problems.

We present the first study of the light hadron spectrum and decay constants
for quenched QCD using an O(a)improved nearestneighbour Wilson fermion action
at \beta=6.2. We compare the results with those obtained using the standard
Wilson fermion action, on the same set of 18 gauge field configurations of a
24^3 times 48 lattice. For pseudoscalar meson masses in the range 330800 MeV,
we find no significant difference between the results for the two actions. The
scales obtained from the string tension and mesonic sector are consistent, but
differ from that derived from baryon masses. The ratio of the pseudoscalar
decay constant to the vector meson mass is roughly independent of quark mass as
observed experimentally, and in approximate agreement with the measured value.