
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.

Composite dark matter candidates, which can arise from new stronglycoupled
sectors, are wellmotivated and phenomenologically interesting, particularly in
the context of asymmetric generation of the relic density. In this work, we
employ lattice calculations to study the electromagnetic form factors of
electroweakneutral darkmatter baryons for a threecolor, QCDlike theory with
Nf = 2 and 6 degenerate fermions in the fundamental representation. We
calculate the (connected) charge radius and anomalous magnetic moment, both of
which can play a significant role for direct detection of composite dark
matter. We find minimal Nf dependence in these quantities. We generate
massdependent crosssections for dark matternucleon interactions and use them
in conjunction with experimental results from XENON100, excluding dark matter
candidates of this type with masses below 10 TeV.

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.

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.

We present an adaptive multigrid Dirac solver developed for Wilson clover
fermions which offers orderofmagnitude reductions in solution time compared
to conventional Krylov solvers. The solver incorporates evenodd
preconditioning and mixed precision to solve the Dirac equation to double
precision accuracy and shows only a mild increase in time to solution for
decreasing quark mass. We show actual time to solution on production lattices
in comparison to conventional Krylov solvers and will also discuss the setup
process and its relative cost to the total solution time.

We present the first results for neutral kaon mixing using 2+1 flavors of
domain wall fermions. A new approach is used to extrapolate to the physical up
and down quark masses from our numerical studies with pion masses in the range
240  420 MeV; only $SU(2)_L \times SU(2)_R$ chiral symmetry is assumed and
the kaon is not assumed to be light. Our main result is $B_K^{\bar{\rm MS}}(2
\mathrm{GeV}) = 0.524(10)(28)$ where the first error is statistical and the
second incorporates estimates for all systematic errors.

We present results for the static interquark potential, light meson and
baryon masses, and light pseudoscalar meson decay constants obtained from
simulations of domain wall QCD with one dynamical flavour approximating the $s$
quark, and two degenerate dynamical flavours with input bare masses ranging
from $m_s$ to $m_s/4$ approximating the $u$ and $d$ quarks. We compare these
quantities obtained using the Iwasaki and DBW2 improved gauge actions, and
actions with larger rectangle coefficients, on $16^3\times32$ lattices. We seek
parameter values at which both the chiral symmetry breaking residual mass due
to the finite lattice extent in the fifth dimension and the Monte Carlo time
history for topological charge are acceptable for this set of quark masses at
lattice spacings above 0.1 fm. We find that the Iwasaki gauge action is best,
demonstrating the feasibility of using QCDOC to generate ensembles which are
good representations of the QCD path integral on lattices of up to 3 fm in
spatial extent with lattice spacings in the range 0.090.13 fm. Despite large
residual masses and a limited number of sea quark mass values with which to
perform chiral extrapolations, our results for light hadronic physics scale and
agree with experimental measurements within our statistical uncertainties.