We consider gravitationally bound states of asymmetric dark matter (ADM stars), and the impact of ADM capture on the stability of neutron stars. We derive and interpret the equation of state for ADM with both attractive and repulsive interactions, and solve the Tolman-Oppenheimer-Volkoff equations to find equilibrium sequences and maximum masses of ADM stars. Gravitational wave searches can utilize our solutions to model exotic compact objects (ECOs). Our results for attractive interactions differ substantially from those in the literature, where fermionic ADM with attractive self-interactions was employed to destabilize neutron stars more effectively than non-interacting fermionic ADM. By contrast, we argue that fermionic ADM with an attractive force is no more effective in destabilizing neutron stars than fermionic ADM with no self-interactions.

Nuggets---very large stable bound objects arising in the presence of a sufficiently attractive and long-range force and in the absence of a dark Coulomb force---are a smoking gun signature for Asymmetric Dark Matter (ADM). The cosmology of ADM nuggets is both generic and unique: nuggets feature highly exothermic fusion processes, which can impact the shape of the core in galaxies, as well as give rise to rare dark star formation. We find, considering the properties of nuggets in a generic extended nuclear model with both attractive and repulsive forces, that self-interaction constraints place an upper bound on nugget masses at the freeze-out of synthesis in the ballpark of $M_{\rm fo} \lesssim 10^{16}$ GeV. We also show that indirect detection strongly constrains models where the scalar mediator binding the nuggets mixes with the Higgs.

We compute the mass function of bound states of Asymmetric Dark Matter--nuggets--synthesized in the early Universe. We apply our results for the nugget density and binding energy computed from a nuclear model to obtain analytic estimates of the typical nugget size exiting synthesis. We numerically solve the Boltzmann equation for synthesis including two-to-two fusion reactions, estimating the impact of bottlenecks on the mass function exiting synthesis. These results provide the basis for studying the late Universe cosmology of nuggets in a future companion paper.

Identifying the true theory of dark matter depends crucially on accurately characterizing interactions of dark matter (DM) with other species. In the context of DM direct detection, we present a study of the prospects for correctly identifying the low-energy effective DM-nucleus scattering operators connected to UV-complete models of DM-quark interactions. We take a census of plausible UV-complete interaction models with different low-energy leading-order DM-nuclear responses. For each model (corresponding to different spin-, momentum-, and velocity-dependent responses), we create a large number of realizations of recoil-energy spectra, and use Bayesian methods to investigate the probability that experiments will be able to select the correct scattering model within a broad set of competing scattering hypotheses. We conclude that agnostic analysis of a strong signal (such as Generation-2 would see if cross sections are just below the current limits) seen on xenon and germanium experiments is likely to correctly identify momentum dependence of the dominant response, ruling out models with either "heavy" or "light" mediators, and enabling downselection of allowed models. However, a unique determination of the correct UV completion will critically depend on the availability of measurements from a wider variety of nuclear targets, including iodine or fluorine. We investigate how model-selection prospects depend on the energy window available for the analysis. In addition, we discuss accuracy of the DM particle mass determination under a wide variety of scattering models, and investigate impact of the specific types of particle-physics uncertainties on prospects for model selection.

We examine the effect of nuclear response functions, as laid out in [Fitzpatrick et al, arXiv:1203.3542], on dark matter (DM) direct detection in the context of well-motivated UV completions, including electric and magnetic dipoles, anapole, spin-orbit, and pseudoscalar-mediated DM. Together, these encompass five of the six nuclear responses extracted from the non-relativistic effective theory of [Fitzpatrick et al, arXiv:1203.3542] (with the sixth difficult to UV complete), with two of the six combinations corresponding to standard spin-independent and -dependent responses. For constraints from existing direct detection experiments, we find that only the COUPP constraint, due to its heavy iodine target with large angular momentum and an unpaired spin, and its large energy range sensitivity, is substantially modified by the new responses compared to what would be inferred using the standard form factors to model the energy dependence of the response. For heavy targets such as xenon and germanium, the behavior of the new nuclear responses as recoil energy increases can be substantially different than that of the standard responses, but this has almost no impact on the constraints derived from experiments such as LUX, XENON100 and CDMS since the maximum nuclear recoil energy detected in these experiments is relatively low. We simulate mock data for 80 and 250 GeV DM candidates utilizing the new nuclear responses to highlight how they might affect a putative signal, and find the new responses are most important for momentum-suppressed interactions such as the magnetic dipole or pseudoscalar-mediated interaction when the target is relatively heavy (such as xenon and iodine).

We examine the consistency of light dark matter (DM) elastic scattering in CoGeNT, DAMA, and CDMS-Silicon in light of constraints from XENON, CDMS, LUX, PICASSO and COUPP. We consider a variety of operators that have been employed to reconcile anomalies with constraints, including anapole, magnetic dipole, momentum-dependent, and isospin-violating DM. We find that elastic scattering through these alternative operators does not substantially reduce the tension between the signals and the null constraints for operators where at least two of the three purported signals map onto a common space in the DM mass--scattering cross-section plane. Taking a choice of the scintillation efficiency that lies at the $-1 \sigma$ region of the Manzur et al measurement relieves tension between signals and the LUX constraint---in particular for a magnetic dipole interaction and a xenophobic interaction (though for the latter the signal regions do not substantially overlap). We also find that modest changes in the halo model does not alter this result. We conclude that, even relaxing the assumption about the type of elastic scattering interaction and taking a conservative choice for the scintillation efficiency, LUX and the results from other null experiments remain in tension with a light DM elastic scattering explanation of direct detection anomalies.

We examine cosmological perturbations in a dynamical theory of inflation in which an Abelian gauge field couples directly to the inflaton, breaking conformal invariance. When the coupling between the gauge field and the inflaton takes a specific form, inflation becomes anisotropic and anisotropy can persist throughout inflation, avoiding Wald's no-hair theorem. After discussing scenarios in which anisotropy can persist during inflation, we calculate the dominant effects of a small persistent anisotropy on the primordial gravitational wave and curvature perturbation power spectra using the "in-in" formalism of perturbation theory. We find that the primordial power spectra of cosmological perturbations gain significant direction dependence and that the fractional direction dependence of the tensor power spectrum is suppressed in comparison to that of the scalar power spectrum.

We study the classical stability of an anisotropic space-time seeded by a spacelike, fixed norm, dynamical vector field in a vacuum-energy-dominated inflationary era. It serves as a model for breaking isotropy during the inflationary era. We find that, for a range of parameters, the linear differential equations for small perturbations about the background do not have a growing mode. We also examine the energy of fluctuations about this background in flat-space. If the kinetic terms for the vector field do not take the form of a field strength tensor squared then there is a negative energy mode and the background is unstable. For the case where the kinetic term is of the form of a field strength tensor squared we show that perturbations about the background have positive energy at lowest order.