Abstract 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 $$\hbox {C}_3\hbox {F}_8$$ 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. (JCAP 1509: 052, 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$$ σ SD ≲ 3 × 10 - 39 cm 2 ( $$6 \times 10^{-38} \mathrm {cm}^2$$ 6 × 10 - 38 cm 2 ) at $$\gtrsim 90\%$$ ≳ 90 % C.L. for a DM particle of mass 1 TeV annihilating into $$\tau ^+ \tau ^-$$ τ + τ - ( $$b\bar{b}$$ b b ¯ ) with a local density of $$\rho _{\mathrm{DM}} = 0.3~\mathrm {GeV/cm}^3$$ ρ DM = 0.3 GeV / cm 3 . The constraints scale inversely with $$\rho _{\mathrm{DM}}$$ ρ DM and are independent of the DM velocity distribution.

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_3F_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.

The identification of non-signal events is a major hurdle to overcome for bubble chamber dark matter experiments such as PICO-60. The current practice of manually developing a discriminator function to eliminate background events is difficult when available calibration data is frequently impure and present only in small quantities. In this study, several different discriminator input/preprocessing formats and neural network architectures are applied to the task. First, they are optimized in a supervised learning context. Next, two novel semi-supervised learning algorithms are trained, and found to replicate the Acoustic Parameter (AP) discriminator previously used in PICO-60 with a mean of 97% accuracy.

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 one 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 weakly interacting massive particle (WIMP)-proton spin-dependent cross section at 3.4×10^{-41}  cm^{2} for a 30-GeV c^{-2} WIMP, more than 1 order of magnitude improvement from previous PICO results.