
The spinindependent and transversity generalised form factors (GFFs) of the
$\phi$ meson are studied using lattice QCD calculations with light quark masses
corresponding to a pion mass $m_\pi\sim450(5)$ MeV. One transversity and three
spinindependent GFFs related to the lowest moments of leadingtwist
spinindependent and transversity gluon distributions are obtained at six
nonzero values of the momentum transfer up to 1.2 GeV$^2$. These quantities
are compared with the analogous spinindependent quark GFFs and the
electromagnetic form factors determined on the same lattice ensemble. The
results show quantitative distinction between the spatial distribution of
transversely polarised gluons, unpolarised gluons, and quarks, and point the
way towards further investigations of the gluon structure of nucleons and
nuclei.

We renew our intent to submit a proposal to perform a search for a nonzero
value of the unmeasured hadronic double helicity flip structure function
$\Delta(x,Q^2)$, predicted to be sensitive to gluons in the nucleus. This would
be performed with an unpolarized electron beam and transversely polarized,
spin1, nuclear target. This structure function was first identified by Jaffe
and Manohar in 1989 as "a clear signature for exotic gluonic components in the
target," and a recent lattice QCD result by our collaborators has prompted
renewed interest in the topic. An inclusive search with deep inelastic
scattering, below $x$ of 0.3, via single spin tensor asymmetries may be
feasible using the CEBAF 12 GeV electron beam and JLab/UVa solid polarized
target, and would represent the first experimental exploration of this
quantity.

We present an exploratory study of the gluonic structure of the $\phi$ meson
using lattice QCD (LQCD). This includes the first investigation of gluonic
transversity via the leading moment of the twisttwo doublehelicityflip
gluonic structure function $\Delta(x,Q^2)$. This structure function only exists
for targets of spin $J\ge1$ and does not mix with quark distributions at
leading twist, thereby providing a particularly clean probe of gluonic degrees
of freedom. We also explore the gluonic analogue of the Soffer bound which
relates the helicity flip and nonflip gluonic distributions, finding it to be
saturated at the level of 80%. This work sets the stage for more complex LQCD
studies of gluonic structure in the nucleon and in light nuclei where
$\Delta(x,Q^2)$ is an 'exotic glue' observable probing gluons in a nucleus not
associated with individual nucleons.

Quarkoniumnucleus systems are composed of two interacting hadronic states
without common valence quarks, which interact primarily through multigluon
exchanges, realizing a color van der Waals force. We present lattice QCD
calculations of the interactions of strange and charm quarkonia with light
nuclei. Both the strangeoniumnucleus and charmoniumnucleus systems are found
to be relatively deeply bound when the masses of the three light quarks are set
equal to that of the physical strange quark. Extrapolation of these results to
the physical lightquark masses suggests that the binding energy of charmonium
to nuclear matter is B < 40 MeV.

It has been argued that the leading scalarisoscalar WIMPnucleus
interactions receive parametrically enhanced contributions in the context of
nuclear effective field theories. These contributions arise from mesonexchange
currents (MECs) and potentially modify the impulse approximation estimates of
these interactions by 1060%. We point out that these MECs also contribute to
the quark mass dependence of nuclear binding energies, that is, nuclear
\sigmaterms. In this work, we use recent lattice QCD calculations of the
binding energies of the deuteron, He3 and He4 at pion masses near 500 MeV and
800 MeV, combined with the experimentally determined binding energies at the
physical point, to provide approximate determinations of the \sigmaterms for
these light nuclei. For each nucleus, we find that the deviation of the
corresponding nuclear \sigmaterm from the singlenucleon estimate is at the
few percent level, in conflict with the conjectured enhancement. As a
consequence, lattice QCD calculations currently indicate that the cross
sections for scalarisoscalar WIMPnucleus interactions arising from
fundamental WIMP interactions with quarks do not suffer from significant
uncertainties due to enhanced mesonexchange currents.

We present the results of lattice QCD calculations of the magnetic moments of
the lightest nuclei, the deuteron, the triton and ${}^3$He, along with those of
the neutron and proton. These calculations, performed at quark masses
corresponding to $m_\pi \sim 800$ MeV, reveal that the structure of these
nuclei at unphysically heavy quark masses closely resembles that at the
physical quark masses. In particular, we find that the magnetic moment of
${}^3$He differs only slightly from that of a free neutron, as is the case in
nature, indicating that the shellmodel configuration of two spinpaired
protons and a valence neutron captures its dominant structure. Similarly a
shellmodellike moment is found for the triton, $\mu_{{}^3{\rm H}} \sim
\mu_p$. The deuteron magnetic moment is found to be equal to the nucleon
isoscalar moment within the uncertainties of the calculations.

The scattering lengths and effective ranges that describe lowenergy
nucleonnucleon scattering are calculated in the limit of SU(3)flavor symmetry
at the physical strangequark mass with Lattice Quantum Chromodynamics. The
calculations are performed with an isotropic clover discretization of the quark
action in three volumes with spatial extents of L \sim 3.4 fm, 4.5fm and 6.7
fm, and with a lattice spacing of b \sim 0.145 fm. With determinations of the
energies of the twonucleon systems (both of which contain bound states at
these up and down quark masses) at rest and moving in the lattice volume,
Luscher's method is used to determine the lowenergy phase shifts in each
channel, from which the scattering length and effective range are obtained. The
scattering parameters, in the 1S0 channel are found to be m_pi a^(1S0) =
9.50^{+0.78}_{0.69}^{+1.10}_{0.80} and m_pi r^(1S0) =
{4.61^{+0.29}_{0.31}^{+0.24}_{0.26}, and in the 3S1 channel are m_pi a^(3S1)
= 7.45^{+0.57}_{0.53}^{+0.71}_{0.49} and m_pi r^(3S1) =
3.71^{+0.28}_{0.31}^{+0.28}_{0.35}. These values are consistent with the
twonucleon system exhibiting Wigner's supermultiplet symmetry, which becomes
exact in the limit of largeN_c. In both spin channels, the phase shifts change
sign at higher momentum, near the start of the tchannel cut, indicating that
the nuclear interactions have a repulsive core even at the SU(3)symmetric
point.

The binding energies of a range of nuclei and hypernuclei with atomic number
A <= 4 and strangeness s <= 2, including the deuteron, dineutron,
Hdibaryon, 3He, Lambda 3He, Lambda 4He, and Lambda Lambda 4He, are calculated
in the limit of flavorSU(3) symmetry at the physical strange quark mass with
quantum chromodynamics (without electromagnetic interactions). The nuclear
states are extracted from Lattice QCD calculations performed with n_f=3
dynamical light quarks using an isotropic clover discretization of the
quarkaction in three lattice volumes of spatial extent L ~ 3.4 fm, 4.5 fm and
6.7 fm, and with a single lattice spacing b ~ 0.145 fm.

An analysis of the pion mass and pion decay constant is performed using
mixedaction Lattice QCD calculations with domainwall valence quarks on
ensembles of rooted, staggered n_f = 2+1 MILC configurations. Calculations were
performed at two lattice spacings of b~0.125 fm and b~0.09 fm, at two strange
quark masses, multiple light quark masses, and a number of lattice volumes. The
ratios of light quark to strange quark masses are in the range 0.1 <= m_l / m_s
<= 0.6, while pion masses are in the range 235 < m_\pi < 680 MeV. A twoflavor
chiral perturbation theory analysis of the Lattice QCD calculations constrains
the GasserLeutwyler coefficients bar{l}_3 and bar{l}_4 to be bar{l}_3 =
4.04(40)(+7355) and bar{l}_4 = 4.30(51)(+8460). All systematic effects in the
calculations are explored, including those from the finite lattice spacetime
volume, the finite lattice spacing, and the finite fifth dimension in the
domainwall quark action. A consistency is demonstrated between a chiral
perturbation theory analysis at fixed lattice spacing combined with a leading
order continuum extrapolation, and the mixedaction chiral perturbation theory
analysis which explicitly includes the leading order discretization effects.
Chiral corrections to the pion decay constant are found to give f_\pi / f =
1.062(26)(+4240) where f is the decay constant in the chiral limit. The most
recent scale setting by the MILC Collaboration yields a postdiction of f_\pi =
128.2(3.6)(+4.46.0)(+1.23.3) MeV at the physical pion mass.

The lowenergy neutronSigma^ interactions determine, in part, the role of
the strange quark in dense matter, such as that found in astrophysical
environments. The scattering phase shifts for this system are obtained from a
numerical evaluation of the QCD path integral using the technique of Lattice
QCD. Our calculations, performed at a pion mass of m_pi ~ 389 MeV in two large
lattice volumes, and at one lattice spacing, are extrapolated to the physical
pion mass using effective field theory. The interactions determined from QCD
are consistent with those extracted from hyperonnucleon experimental data
within uncertainties, and strengthen theoretical arguments that the strange
quark is a crucial component of dense nuclear matter.

The pi+pi+ swave scattering phaseshift is determined below the inelastic
threshold using Lattice QCD. Calculations were performed at a pion mass of
m_pi~390 MeV with an anisotropic n_f=2+1 clover fermion discretization in four
lattice volumes, with spatial extent L~2.0, 2.5, 3.0 and 3.9 fm, and with a
lattice spacing of b_s~0.123 fm in the spatial direction and b_t b_s/3.5 in the
time direction. The phaseshift is determined from the energyeigenvalues of
pi+pi+ systems with both zero and nonzero total momentum in the lattice volume
using Luscher's method. Our calculations are precise enough to allow for a
determination of the threshold scattering parameters, the scattering length a,
the effective range r, and the shapeparameter P, in this channel and to
examine the prediction of twoflavor chiral perturbation theory: m_pi^2 a r =
3+O(m_pi^2/Lambda_chi^2). Chiral perturbation theory is used, with the Lattice
QCD results as input, to predict the scattering phaseshift (and threshold
parameters) at the physical pion mass. Our results are consistent with
determinations from the Roy equations and with the existing experimental phase
shift data.

Results of a highstatistics, multivolume Lattice QCD exploration of the
deuteron, the dineutron, the Hdibaryon, and the XiXi system at a pion mass
of m ~ 390 MeV are presented. Calculations were performed with an anisotropic
n_f = 2+1 Clover discretization in four lattice volumes of spatial extent L ~
2.0, 2.5, 3.0 and 4.0 fm, with a lattice spacing of b_s ~ 0.123 fm in the
spatialdirection, and b_t ~ b_s/3.5 in the timedirection. The XiXi is found
to be bound by B_{XiXi} = 14.0(1.4)(6.7) MeV, consistent with expectations
based upon phenomenological models and lowenergy effective field theories
constrained by nucleonnucleon and hyperonnucleon scattering data at the
physical lightquark masses. We find weak evidence that both the deuteron and
the dineutron are bound at this pion mass, with binding energies of B_d =
11(05)(12) MeV and B_{nn} = 7.1(5.2)(7.3) MeV, respectively. With an increased
number of measurements and a refined analysis, the binding energy of the
Hdibaryon is B_H = 13.2(1.8)(4.0) MeV at this pion mass, updating our previous
result.

The current constraints from lattice QCD on the existence of the Hdibaryon
are discussed. With only two significant lattice QCD calculations of the
Hdibaryon binding energy at approximately the same lattice spacing, the forms
of the chiral and continuum extrapolations to the physical point are not
determined. In this brief report, we consider the constraints on the Hdibaryon
imposed by two simple chiral extrapolations. In both instances, the
extrapolation to the physical pion mass allows for a bound Hdibaryon or a
nearthreshold scattering state. Further lattice QCD calculations are required
to clarify this situation.

Chiral perturbation theory makes definitive predictions for the extrinsic
behavior of hadrons in external electric and magnetic fields. Near the chiral
limit, the electric and magnetic polarizabilities of pions, kaons, and nucleons
are determined in terms of a few wellknown parameters. In this limit, hadrons
become quantum mechanically diffuse as polarizabilities scale with the inverse
squareroot of the quark mass. In some cases, however, such predictions from
chiral perturbation theory have not compared well with experimental data.
Ultimately we must turn to first principles numerical simulations of QCD to
determine properties of hadrons, and confront the predictions of chiral
perturbation theory. To address the electromagnetic polarizabilities, we
utilize the background field technique. Restricting our attention to
calculations in background electric fields, we demonstrate new techniques to
determine electric polarizabilities and baryon magnetic moments for both
charged and neutral states. As we can study the quark mass dependence of
observables with lattice QCD, the lattice will provide a crucial test of our
understanding of lowenergy QCD, which will be timely in light of ongoing
experiments, such as at COMPASS and HI\gamma S.

We present evidence for the existence of a bound Hdibaryon, an I=0, J=0,
s=2 state with valence quark structure uuddss, at a pion mass of m_pi ~ 389
MeV. Using the results of Lattice QCD calculations performed on four ensembles
of anisotropic clover gaugefield configurations, with spatial extents of L ~
2.0, 2.5, 3.0 and 3.9 fm at a spatial lattice spacing of b ~ 0.123 fm, we find
an Hdibaryon bound by B = 16.6 + 2.1 + 4.6 MeV at a pion mass of m_pi ~ 389
MeV.

The volume dependence of the octet baryon masses and relations among them are
explored with Lattice QCD. Calculations are performed with n_f=2+1 clover
fermion discretization in four lattice volumes, with spatial extent L ~ 2.0,
2.5, 3.0 and 3.9 fm, with an anisotropic lattice spacing of b_s ~ 0.123 fm in
the spatial direction, and b_t = b_s/3.5 in the time direction, and at a pion
mass of m_pi ~ 390 MeV. The typical precision of the groundstate baryon mass
determination is ~0.2%, enabling a precise exploration of the volume dependence
of the masses, the GellMannOkubo mass relation, and of other mass
combinations. A comparison of the volume dependence with the predictions of
heavy baryon chiral perturbation theory is performed in both the SU(2)_L X
SU(2)_R and SU(3)_L X SU(3)_R expansions. Predictions of the threeflavor
expansion for the hadron masses are found to describe the observed volume
dependences reasonably well. Further, the DeltaNpi axial coupling constant is
extracted from the volume dependence of the nucleon mass in the twoflavor
expansion, with only small modifications in the threeflavor expansion from the
inclusion of kaons and etas. At a given value of m_pi L, the finitevolume
contributions to the nucleon mass are predicted to be significantly smaller at
m_pi ~ 140 MeV than at m_pi ~ 390 MeV due to a coefficient that scales as ~
m_pi^3. This is relevant for the design of future ensembles of lattice
gaugefield configurations. Finally, the volume dependence of the pion and kaon
masses are analyzed with twoflavor and threeflavor chiral perturbation
theory.

We review recent progress toward establishing lattice Quantum Chromodynamics
as a predictive calculational framework for nuclear physics. A survey of the
current techniques that are used to extract lowenergy hadronic scattering
amplitudes and interactions is followed by a review of recent twobody and
fewbody calculations by the NPLQCD collaboration and others. An outline of the
nuclear physics that is expected to be accomplished with Lattice QCD in the
next decade, along with estimates of the required computational resources, is
presented.

We discuss the extraction of the electromagnetic and spin polarisabilities of
nucleons from lattice QCD. We show that the external field method can be used
to measure all the electromagnetic and spin polarisabilities including those of
charged particles. We then turn to the extrapolations required to connect such
calculations to experiment in the context of chiral perturbation theory,
finding a strong dependence on the lattice volume and quark masses.

We discuss the extraction of the electromagnetic and spin polarisabilities of
nucleons from lattice QCD. We show that the external field method can be used
to measure all the electromagnetic and spin polarisabilities including those of
charged particles. We then turn to the extrapolations required to connect such
calculations to experiment in the context of finite volume chiral perturbation
theory. We derive results relevant for lattice simulations of QCD,
partiallyquenched QCD and quenched QCD. Our results for the polarisabilities
show a strong dependence on the lattice volume and quark masses, typically
differing from the infinite volume limit by ~10% for current lattice volumes
and quark masses.

We study the influence of the strong nuclear vector potential, treated using
the meanfield approximation, in deep inelastic scattering. A consistent
treatment of the electromagnetic current operator, combined with the use of the
operator product expansion is presented and discussed.

We show that hadronic matrix elements can be extracted from lattice
simulations with background fields that arise from operator exponentiation.
Importantly, flavoursinglet matrix elements can be evaluated without requiring
the computation of disconnected diagrams, thus facilitating a calculation of
the quark contribution to the spin of the proton and the singlet axial
coupling, $g_A^0$. In the two nucleon sector, a background field approach will
allow calculation of the magnetic and quadrupole moments of the deuteron and an
investigation of the EMC effect directly from lattice QCD. Matrix elements
between states of differing momenta are also analysed in the presence of
background fields.

The results of different nonperturbative studies agree on a power law as the
infrared behavior of the Landau gauge gluon propagator. This propagator
violates positivity and thus indicates the absence of the transverse gluons
from the physical spectrum, i.e. gluon confinement. A simple analytic structure
for the gluon propagator is proposed capturing all of its features. We comment
also on related investigations for the Landau gauge quark propagator.

We derive double distributions for the proton in a simple model that contains
scalar as well as axialvector diquark correlations. The model parameters are
tuned so that the experimentally measured electromagnetic form factors are
reproduced for small momentum transfer. Resulting generalized parton
distributions satisfy known constraints, including the positivity bounds.

We explore the analytic structure of the gluon and quark propagators of
Landau gauge QCD from numerical solutions of the coupled system of renormalized
DysonSchwinger equations and from fits to lattice data. We find sizable
negative norm contributions in the transverse gluon propagator indicating the
absence of the transverse gluon from the physical spectrum. A simple analytic
structure for the gluon propagator is proposed. For the quark propagator we
find evidence for a masslike singularity on the real timelike momentum axis,
with a mass of 350 to 500 MeV. Within the employed Green's functions approach
we identify a crucial term in the quarkgluon vertex that leads to a positive
definite Schwinger function for the quark propagator.

In Landau gauge QCD the infrared behavior of the propagator of transverse
gluons can be analytically determined to be a power law from DysonSchwinger
equations. This propagator clearly shows positivity violation, indicating the
absence of the transverse gluons from the physical spectrum, i.e. gluon
confinement. A simple analytic structure for the gluon propagator is proposed
capturing all important features. We provide arguments that the Landau gauge
quark propagator possesses a singularity on the real timelike axis. For this
propagator we find a positive definite Schwinger function.