Adopting the Standard Halo Model (SHM) of an isotropic Maxwellian velocity distribution for dark matter (DM) particles in the Galaxy, the most stringent current constraints on their spin-dependent scattering cross-section with nucleons come from the IceCube neutrino observatory and the PICO-60 C$_3$F$_8$ superheated bubble chamber experiments. The former is sensitive to high energy neutrinos from the self-annihilation of DM particles captured in the Sun, while the latter looks for nuclear recoil events from DM scattering off nucleons. Although slower DM particles are more likely to be captured by the Sun, the faster ones are more likely to be detected by PICO. Recent N-body simulations suggest significant deviations from the SHM for the smooth halo component of the DM, while observations hint at a dominant fraction of the local DM being in substructures. We use the method of Ferrer et al. (2015) to exploit the complementarity between the two approaches and derive conservative constraints on DM-nucleon scattering. Our results constrain $\sigma_{\mathrm{SD}} \lesssim 3 \times 10^{-39} \mathrm{cm}^2$ (6 $ \times 10^{-38} \mathrm{cm}^2$) at $\gtrsim 90\%$ C.L. for a DM particle of mass 1~TeV annihilating into $\tau^+ \tau^-$ ($b\bar{b}$) with a local density of $\rho_{\mathrm{DM}} = 0.3~\mathrm{ GeV/cm}^3$. The constraints scale inversely with $\rho_{\mathrm{DM}}$ and are independent of the DM velocity distribution.

We report direct-detection constraints on light dark matter particles interacting with electrons. The results are based on a method that exploits the extremely low levels of leakage current of the DAMIC detector at SNOLAB of 2-6$\times$10$^{-22}$ A cm$^{-2}$. We evaluate the charge distribution of pixels that collect $<10~\rm{e^-}$ for contributions beyond the leakage current that may be attributed to dark matter interactions. Constraints are placed on so-far unexplored parameter space for dark matter masses between 0.6 and 100 MeV$c^{-2}$. We also present new constraints on hidden-photon dark matter with masses in the range $1.2$-$30$ eV$c^{-2}$.

We present a novel unifying interpretation of excess event rates observed in several dark matter direct-detection experiments that utilize single-electron threshold semiconductor detectors. Despite their different locations, exposures, readout techniques, detector composition, and operating depths, these experiments all observe statistically significant excess event rates of $\sim$ 10 Hz/kg. However, none of these persistent excesses has yet been reported as a dark matter signal because individually, each can be attributed to different well-motivated but unmodeled backgrounds, and taken together, they cannot be explained by dark matter particles scattering elastically off detector nuclei or electrons. We show that these results can be reconciled if the semiconductor detectors are seeing a collective inelastic process, consistent with exciting a plasmon. We further show that plasmon excitation could arise in two compelling dark matter scenarios, both of which can explain rates of existing signal excesses in germanium and, at least at the order of magnitude level, across several single-electron threshold detectors. At least one of these scenarios also yields the correct relic density from thermal freeze-out. Both dark matter scenarios motivate a radical rethinking of the standard interpretations of dark matter-electron scattering from recent experiments.

We respond to Robinson and Michaud's (RM) comment (arXiv:2002.08893) on our recent preprint arXiv:2002.06937, in which we discuss recent excesses in low-threshold dark matter searches, and offer a potential unifying dark matter interpretation. We thank RM for their feedback, which highlights the critical need for future measurements to directly calibrate plasmon charge yields for low $\sim$ 10 eV energy depositions. RM objected to our assertion that plasmons generated at energy scales below 100~eV may have a large branching fraction into phonons. As we argue below, the points raised by RM do not invalidate our primary conclusions, as they pertain to a much different energy scale than we discuss in our paper.

The primary advantage of moderately superheated bubble chamber detectors is their simultaneous sensitivity to nuclear recoils from WIMP dark matter and insensitivity to electron recoil backgrounds. A comprehensive analysis of PICO gamma calibration data demonstrates for the first time that electron recoils in C$_3$F$_8$ scale in accordance with a new nucleation mechanism, rather than one driven by a hot-spike as previously supposed. Using this semi-empirical model, bubble chamber nucleation thresholds may be tuned to be sensitive to lower energy nuclear recoils while maintaining excellent electron recoil rejection. The PICO-40L detector will exploit this model to achieve thermodynamic thresholds as low as 2.8 keV while being dominated by single-scatter events from coherent elastic neutrino-nucleus scattering of solar neutrinos. In one year of operation, PICO-40L can improve existing leading limits from PICO on spin-dependent WIMP-proton coupling by nearly an order of magnitude for WIMP masses greater than 3 GeV c$^{-2}$ and will have the ability to surpass all existing non-xenon bounds on spin-independent WIMP-nucleon coupling for WIMP masses from 3 to 40 GeV c$^{-2}$.

There are currently several existing and proposed experiments designed to probe sub-GeV dark matter (DM) using electron ionization in various materials. The projected signal rates for these experiments assume that this ionization yield arises only from DM scattering directly off electron targets, ignoring secondary ionization contributions from DM scattering off nuclear targets. In this work we investigate the validity of this assumption and show that if sub-GeV DM couples with comparable strength to both protons and electrons, as would be the case for a dark photon mediator, the ionization signal from atomic scattering via the Migdal effect scales with $Z^2 (m_e/m_N)^2 \mathbf{q}^2$ where $m_N$ is the mass of the target nucleus and $\mathbf{q}$ is the 3-momentum transferred to the atom. The result is that the Migdal effect is always subdominant to electron scattering when the mediator is light, but that Migdal-induced ionization can dominate over electron scattering for heavy mediators and DM masses in the hundreds of MeV range. We put these two ionization processes on identical theoretical footing, address some theoretical uncertainties in the choice of atomic wavefunctions used to compute rates, and discuss the implications for DM scenarios where the Migdal process dominates, including for XENON10, XENON100, and the recent XENON1T results on light DM scattering.

Final results are reported from operation of the PICO-60 C$_3$F$_8$ dark matter detector, a bubble chamber filled with 52 kg of C$_3$F$_8$ located in the SNOLAB underground laboratory. The chamber was operated at thermodynamic thresholds as low as 1.2 keV without loss of stability. A new blind 1404-kg-day exposure at 2.45 keV threshold was acquired with approximately the same expected total background rate as the previous 1167-kg-day exposure at 3.3 keV. This increased exposure is enabled in part by a new optical tracking analysis to better identify events near detector walls, permitting a larger fiducial volume. These results set the most stringent direct-detection constraint to date on the WIMP-proton spin-dependent cross section at 2.5 $\times$ 10$^{-41}$ cm$^2$ for a 25 GeV WIMP, and improve on previous PICO results for 3-5 GeV WIMPs by an order of magnitude.

New results are reported from the operation of the PICO-60 dark matter detector, a bubble chamber filled with 52 kg of C$_3$F$_8$ located in the SNOLAB underground laboratory. As in previous PICO bubble chambers, PICO-60 C$_3$F$_8$ exhibits excellent electron recoil and alpha decay rejection, and the observed multiple-scattering neutron rate indicates a single-scatter neutron background of less than 1 event per month. A blind analysis of an efficiency-corrected 1167-kg-day exposure at a 3.3-keV thermodynamic threshold reveals no single-scattering nuclear recoil candidates, consistent with the predicted background. These results set the most stringent direct-detection constraint to date on the WIMP-proton spin-dependent cross section at 3.4 $\times$ 10$^{-41}$ cm$^2$ for a 30-GeV$\thinspace$c$^{-2}$ WIMP, more than one order of magnitude improvement from previous PICO results.

A 30-g xenon bubble chamber, operated at Northwestern University in June and November 2016, has for the first time observed simultaneous bubble nucleation and scintillation by nuclear recoils in a superheated liquid. This chamber is instrumented with a CCD camera for near-IR bubble imaging, a solar-blind photomultiplier tube to detect 175-nm xenon scintillation light, and a piezoelectric acoustic transducer to detect the ultrasonic emission from a growing bubble. The time of nucleation determined from the acoustic signal is used to correlate specific scintillation pulses with bubble-nucleating events. We report on data from this chamber for thermodynamic "Seitz" thresholds from 4.2 to 15.0 keV. The observed single- and multiple-bubble rates when exposed to a $^{252}$Cf neutron source indicate that, for an 8.3-keV thermodynamic threshold, the minimum nuclear recoil energy required to nucleate a bubble is $19\pm6$ keV (1$\sigma$ uncertainty). This is consistent with the observed scintillation spectrum for bubble-nucleating events. We see no evidence for bubble nucleation by gamma rays at any of the thresholds studied, setting a 90% C.L. upper limit of $6.3\times10^{-7}$ bubbles per gamma interaction at a 4.2-keV thermodynamic threshold. This indicates stronger gamma discrimination than in CF$_3$I bubble chambers, supporting the hypothesis that scintillation production suppresses bubble nucleation by electron recoils while nuclear recoils nucleate bubbles as usual. These measurements establish the noble-liquid bubble chamber as a promising new technology for the detection of weakly interacting massive particle dark matter and coherent elastic neutrino-nucleus scattering.

New data are reported from a second run of the 2-liter PICO-2L C$_3$F$_8$ bubble chamber with a total exposure of 129$\,$kg-days at a thermodynamic threshold energy of 3.3$\,$keV. These data show that measures taken to control particulate contamination in the superheated fluid resulted in the absence of the anomalous background events observed in the first run of this bubble chamber. One single nuclear-recoil event was observed in the data, consistent both with the predicted background rate from neutrons and with the observed rate of unambiguous multiple-bubble neutron scattering events. The chamber exhibits the same excellent electron-recoil and alpha decay rejection as was previously reported. These data provide the most stringent direct detection constraints on weakly interacting massive particle (WIMP)-proton spin-dependent scattering to date for WIMP masses $<$ 50$\,$GeV/c$^2$.

New data are reported from the operation of the PICO-60 dark matter detector, a bubble chamber filled with 36.8 kg of CF$_3$I and located in the SNOLAB underground laboratory. PICO-60 is the largest bubble chamber to search for dark matter to date. With an analyzed exposure of 92.8 livedays, PICO-60 exhibits the same excellent background rejection observed in smaller bubble chambers. Alpha decays in PICO-60 exhibit frequency-dependent acoustic calorimetry, similar but not identical to that reported recently in a C$_3$F$_8$ bubble chamber. PICO-60 also observes a large population of unknown background events, exhibiting acoustic, spatial, and timing behaviors inconsistent with those expected from a dark matter signal. These behaviors allow for analysis cuts to remove all background events while retaining $48.2\%$ of the exposure. Stringent limits on weakly interacting massive particles interacting via spin-dependent proton and spin-independent processes are set, and most interpretations of the DAMA/LIBRA modulation signal as dark matter interacting with iodine nuclei are ruled out.

New data are reported from the operation of a 2-liter C$_3$F$_8$ bubble chamber in the 2100 meter deep SNOLAB underground laboratory, with a total exposure of 211.5 kg-days at four different recoil energy thresholds ranging from 3.2 keV to 8.1 keV. These data show that C3F8 provides excellent electron recoil and alpha rejection capabilities at very low thresholds, including the first observation of a dependence of acoustic signal on alpha energy. Twelve single nuclear recoil event candidates were observed during the run. The candidate events exhibit timing characteristics that are not consistent with the hypothesis of a uniform time distribution, and no evidence for a dark matter signal is claimed. These data provide the most sensitive direct detection constraints on WIMP-proton spin-dependent scattering to date, with significant sensitivity at low WIMP masses for spin-independent WIMP-nucleon scattering.