One intriguing BSM particle is the QCD axion, which could simultaneously provide a solution to the Strong CP problem and account for some, if not all, of the dark matter density in the universe. This particle is a pNGB of the conjectured Peccei-Quinn (PQ) symmetry of the Standard Model. Its mass and interactions are suppressed by a heavy symmetry breaking scale, $f_a$, whose value is roughly greater than $10^{9}$ GeV (or, conversely, the axion mass, $m_a$, is roughly less than $10^4\ \mu \text{eV}$). The density of axions in the universe, which cannot exceed the relic dark matter density and is a quantity of great interest in axion experiments like ADMX, is a result of the early-universe interplay between cosmological evolution and the axion mass as a function of temperature. The latter quantity is proportional to the second derivative of the QCD free energy with respect to the CP-violating phase, $\theta$. However, this quantity is generically non-perturbative and previous calculations have only employed instanton models at the high temperatures of interest (roughly 1 GeV). In this and future works, we aim to calculate the temperature-dependent axion mass at small $\theta$ from first-principle lattice calculations, with controlled statistical and systematic errors. Once calculated, this temperature-dependent axion mass is input for the classical evolution equations of the axion density of the universe. Due to a variety of lattice systematic effects at the very high temperatures required, we perform a calculation of the leading small-$\theta$ cumulant of the theta vacua on large volume lattices for SU(3) Yang-Mills with high statistics as a first proof of concept, before attempting a full QCD calculation in the future. From these pure glue results, the misalignment mechanism yields the axion mass bound $m_a \geq (14.6\pm0.1) \ \mu \text{eV}$ when PQ-breaking occurs after inflation.

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.

We consider deuterium compressed to higher than atomic, but lower than nuclear densities. At such densities deuterium is a superconducting quantum liquid. Generically, two superconducting phases compete, a "ferromagnetic" and a "nematic" one. We provide a power counting argument suggesting that the dominant interactions in the deuteron liquid are perturbative (but screened) Coulomb interactions. At very high densities the ground state is determined by very small nuclear interaction effects that probably favor the ferromagnetic phase. At lower densities the symmetry of the theory is effectively enhanced to SU(3), and the quantum liquid enters a novel phase, neither ferromagnetic nor nematic. Our results can serve as a starting point for investigations of the phase dynamics of deuteron liquids, as well as exploration of the stability and dynamics of the rich variety of topological objects that may occur in phases of the deuteron quantum liquid, which range from Alice strings to spin skyrmions to Z_2 vortices.

We point out that the exchange of a Goldstone pseudo-scalar can provide an enhancement in the dark matter annihilation rate capable of explaining the excess flux seen in high energy cosmic ray data. The mechanism of enhancement involves the coupling of s and d waves through the tensor force that is very strong and, in fact, singular at short distances. The results indicate that large enhancements require some amount of fine tuning. We also discuss the enhancement due to other singular attractive potentials, such as WIMP models with a permanent electric dipole.