Quarkonium-nucleus systems are composed of two interacting hadronic states without common valence quarks, which interact primarily through multi-gluon 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 strangeonium-nucleus and charmonium-nucleus 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 light-quark masses suggests that the binding energy of charmonium to nuclear matter is B < 40 MeV.

It has been argued that the leading scalar-isoscalar WIMP-nucleus interactions receive parametrically enhanced contributions in the context of nuclear effective field theories. These contributions arise from meson-exchange currents (MECs) and potentially modify the impulse approximation estimates of these interactions by 10--60%. We point out that these MECs also contribute to the quark mass dependence of nuclear binding energies, that is, nuclear \sigma-terms. In this work, we use recent lattice QCD calculations of the binding energies of the deuteron, He-3 and He-4 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 \sigma-terms for these light nuclei. For each nucleus, we find that the deviation of the corresponding nuclear \sigma-term from the single-nucleon 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 scalar-isoscalar WIMP-nucleus interactions arising from fundamental WIMP interactions with quarks do not suffer from significant uncertainties due to enhanced meson-exchange currents.

Composite dark matter candidates, which can arise from new strongly-coupled sectors, are well-motivated 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 electroweak-neutral dark-matter baryons for a three-color, QCD-like 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 mass-dependent cross-sections for dark matter-nucleon 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 low-energy nucleon-nucleon scattering are calculated in the limit of SU(3)-flavor symmetry at the physical strange-quark 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 two-nucleon 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 low-energy 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 two-nucleon system exhibiting Wigner's supermultiplet symmetry, which becomes exact in the limit of large-N_c. In both spin channels, the phase shifts change sign at higher momentum, near the start of the t-channel 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, di-neutron, H-dibaryon, 3He, Lambda 3He, Lambda 4He, and Lambda Lambda 4He, are calculated in the limit of flavor-SU(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 quark-action 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 low-energy neutron-Sigma^- 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 hyperon-nucleon 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 order-of-magnitude reductions in solution time compared to conventional Krylov solvers. The solver incorporates even-odd 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.09-0.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.