This article presents an analysis and the resulting limits on light dark matter inelastically scattering off of electrons, and on dark photon and axion-like particle absorption, using a second-generation SuperCDMS high-voltage eV-resolution detector. The 0.93 gram Si detector achieved a 3 eV phonon energy resolution; for a detector bias of 100 V, this corresponds to a charge resolution of 3% of a single electron-hole pair. The energy spectrum is reported from a blind analysis with 1.2 gram-days of exposure acquired in an above-ground laboratory. With charge carrier trapping and impact ionization effects incorporated into the dark matter signal models, the dark matter-electron cross section $\bar{\sigma}_{e}$ is constrained for dark matter masses from 0.5--$10^{4} $MeV$/c^{2}$; in the mass range from 1.2--50 eV$/c^{2}$ the dark photon kinetic mixing parameter $\varepsilon$ and the axioelectric coupling constant $g_{ae}$ are constrained. The minimum 90% confidence-level upper limits within the above mentioned mass ranges are $\bar{\sigma}_{e}\,=\,8.7\times10^{-34}$ cm$^{2}$, $\varepsilon\,=\,3.3\times10^{-14}$, and $g_{ae}\,=\,1.0\times10^{-9}$.

Calculating the abundance of thermally produced dark matter particles has become a standard procedure, with sophisticated methods guaranteeing a precision that matches the percent-level accuracy in the observational determination of the dark matter density. Here, we point out that one of the main assumptions in the commonly adopted formalism, namely local thermal equilibrium during the freeze-out of annihilating dark matter particles, does not have to be satisfied in general. We present two methods for how to deal with such situations, in which the kinetic decoupling of dark matter happens so early that it interferes with the chemical decoupling process: i) an approximate treatment in terms of a coupled system of differential equations for the leading momentum moments of the dark matter distribution, and ii) a full numerical solution of the Boltzmann equation in phase-space. For illustration, we apply these methods to the case of Scalar Singlet dark matter. We explicitly show that even in this simple model the prediction for the dark matter abundance can be affected by up to one order of magnitude compared to the traditional treatment.

The thermal decoupling description of dark matter (DM) and co-annihilating partners is reconsidered. If DM is realized at around the TeV-mass region or above, even the heaviest electroweak force carriers could act as long-range forces, leading to the existence of meta-stable DM bound states. The formation and subsequent decay of the latter further deplete the relic density during the freeze-out process on top of the Sommerfeld enhancement, allowing for larger DM masses. While so far the bound-state formation was described via the emission of an on-shell mediator ($W^{\pm}$, $Z$, $H$, $g$, photon or exotic), we point out that this particular process does not have to be the dominant scattering-bound state conversion channel in general. If the mediator is coupled in a direct way to any relativistic species present in the Early Universe, the bound-state formation can efficiently occur through particle scattering, where a mediator is exchanged virtually. To demonstrate that such a virtually stimulated conversion process can dominate the on-shell emission even for all temperatures, we analyze a simplified model where DM is coupled to only one relativistic species in the primordial plasma through an electroweak-scale mediator. We find that the bound-state formation cross section via particle scattering can exceed the on-shell emission by up to several orders of magnitude.

Future direct searches for low-mass dark matter particles with germanium detectors, such as SuperCDMS SNOLAB, are expected to be limited by backgrounds from radioactive isotopes activated by cosmogenic radiation inside the germanium. There are limited experimental data available to constrain production rates and a large spread of theoretical predictions. We examine the calculation of expected production rates, and analyze data from the second run of the CDMS low ionization threshold experiment (CDMSlite) to estimate the rates for several isotopes. We model the measured CDMSlite spectrum and fit for contributions from tritium and other isotopes. Using the knowledge of the detector history, these results are converted to cosmogenic production rates at sea level. The production rates in atoms/(kg$\cdot$day) are 74$\pm$9 for $^3$H, 1.5$\pm$0.7 for $^{55}$Fe, 17$\pm$5 for $^{65}$Zn, and 30$\pm$18 for $^{68}$Ge.

We present the first limits on inelastic electron-scattering dark matter and dark photon absorption using a prototype SuperCDMS detector having a charge resolution of 0.1 electron-hole pairs (CDMS HVeV, a 0.93 gram CDMS HV device). These electron-recoil limits significantly improve experimental constraints on dark matter particles with masses as low as 1 MeV/$\mathrm{c^2}$. We demonstrate a sensitivity to dark photons competitive with other leading approaches but using substantially less exposure (0.49 gram days). These results demonstrate the scientific potential of phonon-mediated semiconductor detectors that are sensitive to single electronic excitations.

The Cryogenic Dark Matter Search low ionization threshold experiment (CDMSlite) searches for interactions between dark matter particles and germanium nuclei in cryogenic detectors. The experiment has achieved a low energy threshold with improved sensitivity to low-mass (<10 GeV/c$^2$) dark matter particles. We present an analysis of the final CDMSlite data set, taken with a different detector than was used for the two previous CDMSlite data sets. This analysis includes a data "salting" method to protect against bias, improved noise discrimination, background modeling, and the use of profile likelihood methods to search for a dark matter signal in the presence of backgrounds. We achieve an energy threshold of 70 eV and significantly improve the sensitivity for dark matter particles with masses between 2.5 and 10 GeV/c$^2$ compared to previous analyses. We set an upper limit on the dark matter-nucleon scattering cross section in germanium of 5.4$\times$10$^{-42}$ cm$^2$ at 5 GeV/c$^2$, a factor of $\sim$2.5 improvement over the previous CDMSlite result.

The Super Cryogenic Dark Matter Search experiment (SuperCDMS) at the Soudan Underground Laboratory studied energy loss associated with Frenkel defect formation in germanium crystals at mK temperatures using in situ $^{210}$Pb sources. We examine the spectrum of $^{206}$Pb nuclear recoils near its expected 103 keV endpoint energy and determine an energy loss of $\left(6.08\pm0.18\right)$ %, which we attribute to defect formation. From this result and using TRIM simulations, we extract the first experimentally determined average displacement threshold energy of $\left(19.7^{+0.6}_{-0.5}\right)$ eV for germanium. This has implications for the analysis thresholds of future germanium-based dark matter searches.

Traditional computations of the dark matter (DM) relic abundance, for models where attractive self-interactions are mediated by light force-carriers and bound states exist, rely on the solution of a coupled system of classical on-shell Boltzmann equations. This idealized description misses important thermal effects caused by the tight coupling among force-carriers and other charged relativistic species. We develop for the first time a comprehensive ab-initio derivation for the description of DM long-range interactions in the presence of a hot and dense plasma background directly from non-equilibrium quantum field theory. Most importantly, the scattering and bound states get strongly mixed in the thermal plasma environment, representing a characteristic difference from a pure vacuum theory computation. The main result of this work is a novel differential equation for the DM number density, written down in a form which is manifestly independent under the choice of what one would interpret as a bound or a scattering state at finite temperature. The collision term, unifying the description of annihilation and bound state decay, turns out to have in general a non-quadratic dependence on the DM number density. This generalizes the form of the conventional Lee-Weinberg equation which is typically adopted to describe the freeze-out process. We prove that our general number density equation is consistent with previous literature results under certain limits.

The Cryogenic Dark Matter Search (CDMS II) experiment aims to detect dark matter particles that elastically scatter from nuclei in semiconductor detectors. The resulting nuclear-recoil energy depositions are detected by ionization and phonon sensors. Neutrons produce a similar spectrum of low-energy nuclear recoils in such detectors, while most other backgrounds produce electron recoils. The absolute energy scale for nuclear recoils is necessary to interpret results correctly. The energy scale can be determined in CDMS II silicon detectors using neutrons incident from a broad-spectrum $^{252}$Cf source, taking advantage of a prominent resonance in the neutron elastic scattering cross section of silicon at a recoil (neutron) energy near 20 (182) keV. Results indicate that the phonon collection efficiency for nuclear recoils is $4.8^{+0.7}_{-0.9}$% lower than for electron recoils of the same energy. Comparisons of the ionization signals for nuclear recoils to those measured previously by other groups at higher electric fields indicate that the ionization collection efficiency for CDMS II silicon detectors operated at $\sim$4 V/cm is consistent with 100% for nuclear recoils below 20 keV and gradually decreases for larger energies to $\sim$75% at 100 keV. The impact of these measurements on previously published CDMS II silicon results is small.

Kinetic decoupling of dark matter typically happens much later than chemical freeze-out. In fact, local thermal equilibrium is an important assumption for the usual relic density calculations based on solving the Boltzmann equation (for its 0-th moment) describing the dark matter number density. But is this assumption always justified? Here we address this question and discuss the consequences of more accurate treatments. The first treatment is relying on the inclusion of higher moments of the Boltzmann equation and the second on solving the evolution of the phase-space distribution function fully numerically. For illustration, these methods are applied to the Scalar Singlet model, often referred to as the simplest WIMP DM possibility from a model-building perspective. It is explicitly shown that even in this simple model the prediction for the dark matter abundance can be affected by as much as one order of magnitude.

We explore the phenomenology of having a second epoch of dark matter annihilation into dark radiation long after the standard thermal freeze-out. Such a hidden reannihilation process could affect visible sectors only gravitationally. As a concrete realization we consider self-interacting dark matter (SIDM) with a light force mediator coupled to dark radiation. We demonstrate how resonantly Sommerfeld enhanced cross sections emerge to induce the reannihilation epoch. The effect is a temporally local modification of the Hubble expansion rate and we show that the Cosmic Microwave Background (CMB) measurements -- as well as other observations -- have a high sensitivity to observe this phenomenon. Special attention is given to the model region where late kinetic decoupling and strong self-interactions can alleviate several small-scale problems in the cold dark matter paradigm at the same time. Interestingly, we find that reannihilation might here also simultaneously lower the tension between CMB and low-redshift astronomical observations of $H_{0}$ and $\sigma_{8}$. Moreover, we identify reannihilation as a clear signature to discriminate between the phenomenologically otherwise almost identical vector and scalar mediator realizations of SIDM.

The SuperCDMS experiment is designed to directly detect weakly interacting massive particles (WIMPs) that may constitute the dark matter in our Galaxy. During its operation at the Soudan Underground Laboratory, germanium detectors were run in the CDMSlite mode to gather data sets with sensitivity specifically for WIMPs with masses ${<}$10 GeV/$c^2$. In this mode, a higher detector-bias voltage is applied to amplify the phonon signals produced by drifting charges. This paper presents studies of the experimental noise and its effect on the achievable energy threshold, which is demonstrated to be as low as 56 eV$_{\text{ee}}$ (electron equivalent energy). The detector-biasing configuration is described in detail, with analysis corrections for voltage variations to the level of a few percent. Detailed studies of the electric-field geometry, and the resulting successful development of a fiducial parameter, eliminate poorly measured events, yielding an energy resolution ranging from ${\sim}$9 eV$_{\text{ee}}$ at 0 keV to 101 eV$_{\text{ee}}$ at ${\sim}$10 eV$_{\text{ee}}$. New results are derived for astrophysical uncertainties relevant to the WIMP-search limits, specifically examining how they are affected by variations in the most probable WIMP velocity and the Galactic escape velocity. These variations become more important for WIMP masses below 10 GeV/$c^2$. Finally, new limits on spin-dependent low-mass WIMP-nucleon interactions are derived, with new parameter space excluded for WIMP masses $\lesssim$3 GeV/$c^2$

Dark Matter (DM) models providing possible alternative solutions to the small- scale crisis of standard cosmology are nowadays of growing interest. We consider DM interacting with light hidden fermions via well motivated fundamental operators showing the resultant matter power spectrum is suppressed on subgalactic scales within a plausible parameter region. Our basic description of the evolution of cosmological perturbations relies on a fully consistent first principles derivation of a perturbed Fokker-Planck type equation, generalizing existing literature. The cosmological perturbation of the Fokker-Planck equation is presented for the first time in two different gauges, where the results transform into each other according to the rules of gauge transformation. Furthermore, our focus lies on a derivation of a broadly applicable and easily computable collision term showing important phenomenological differences to other existing approximations. As one of the main results and concerning the small-scale crisis, we show the equal importance of vector and scalar boson mediated interactions between DM and light fermions.