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.

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}$.

In this paper we report on the characterization of the Hamamatsu VUV4 (S/N: S13370-6152) Vacuum Ultra-Violet (VUV) sensitive Silicon Photo-Multipliers (SiPMs) as part of the development of a solution for the detection of liquid xenon scintillation light for the nEXO experiment. Various SiPM features, such as: dark noise, gain, correlated avalanches, direct crosstalk and Photon Detection Efficiency (PDE) were measured in a dedicated setup at TRIUMF. SiPMs were characterized in the range $163 \text{ } \text{K} \leq \text{T}\leq 233 \text{ } \text{K}$. At an over voltage of $3.1\pm0.2$ V and at $\text{T}=163 \text{ }\text{K}$ we report a number of Correlated Avalanches (CAs) per pulse in the $1 \upmu\text{s}$ interval following the trigger pulse of $0.161\pm0.005$. At the same settings the Dark-Noise (DN) rate is $0.137\pm0.002 \text{ Hz/mm}^{2}$. Both the number of CAs and the DN rate are within nEXO specifications. The PDE of the Hamamatsu VUV4 was measured for two different devices at $\text{T}=233 \text{ }\text{K}$ for a mean wavelength of $189\pm7\text{ nm}$. At $3.6\pm0.2$ V and $3.5\pm0.2$ V of over voltage we report a PDE of $13.4\pm2.6\text{ }\%$ and $11\pm2\%$, corresponding to a saturation PDE of $14.8\pm2.8\text{ }\%$ and $12.2\pm2.3\%$, respectively. Both values are well below the $24\text{ }\%$ saturation PDE advertised by Hamamatsu. More generally, the second device tested at $3.5\pm0.2$ V of over voltage is below the nEXO PDE requirement. The first one instead yields a PDE that is marginally close to meeting the nEXO specifications. This suggests that with modest improvements the Hamamatsu VUV4 MPPCs could be considered as an alternative to the FBK-LF SiPMs for the final design of the nEXO detector.

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.

The PICASSO dark matter search experiment operated an array of 32 superheated droplet detectors containing 3.0 kg of C$_{4}$F$_{10}$ and collected an exposure of 231.4 kgd at SNOLAB between March 2012 and January 2014. We report on the final results of this experiment which includes for the first time the complete data set and improved analysis techniques including \mbox{acoustic} localization to allow fiducialization and removal of higher activity regions within the detectors. No signal consistent with dark matter was observed. We set limits for spin-dependent interactions on protons of $\sigma_p^{SD}$~=~1.32~$\times$~10$^{-2}$~pb (90\%~C.L.) at a WIMP mass of 20 GeV/c$^{2}$. In the spin-independent sector we exclude cross sections larger than $\sigma_p^{SI}$~=~4.86~$\times$~10$^{-5 }$~pb~(90\% C.L.) in the region around 7 GeV/c$^{2}$. The pioneering efforts of the PICASSO experiment have paved the way forward for a next generation detector incorporating much of this technology and experience into larger mass bubble chambers.

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.

The Sudbury Neutrino Observatory (SNO) has confirmed the standard solar model and neutrino oscillations through the observation of neutrinos from the solar core. In this paper we present a search for neutrinos associated with sources other than the solar core, such as gamma-ray bursters and solar flares. We present a new method for looking for temporal coincidences between neutrino events and astrophysical bursts of widely varying intensity. No correlations were found between neutrinos detected in SNO and such astrophysical sources.

We report results from a combined analysis of solar neutrino data from all phases of the Sudbury Neutrino Observatory. By exploiting particle identification information obtained from the proportional counters installed during the third phase, this analysis improved background rejection in that phase of the experiment. The combined analysis resulted in a total flux of active neutrino flavors from 8B decays in the Sun of (5.25 \pm 0.16(stat.)+0.11-0.13(syst.))\times10^6 cm^{-2}s^{-1}. A two-flavor neutrino oscillation analysis yielded \Deltam^2_{21} = (5.6^{+1.9}_{-1.4})\times10^{-5} eV^2 and tan^2{\theta}_{12}= 0.427^{+0.033}_{-0.029}. A three-flavor neutrino oscillation analysis combining this result with results of all other solar neutrino experiments and the KamLAND experiment yielded \Deltam^2_{21} = (7.41^{+0.21}_{-0.19})\times10^{-5} eV^2, tan^2{\theta}_{12} = 0.446^{+0.030}_{-0.029}, and sin^2{\theta}_{13} = (2.5^{+1.8}_{-1.5})\times10^{-2}. This implied an upper bound of sin^2{\theta}_{13} < 0.053 at the 95% confidence level (C.L.).

This paper details the solar neutrino analysis of the 385.17-day Phase-III data set acquired by the Sudbury Neutrino Observatory (SNO). An array of $^3$He proportional counters was installed in the heavy-water target to measure precisely the rate of neutrino-deuteron neutral-current interactions. This technique to determine the total active $^8$B solar neutrino flux was largely independent of the methods employed in previous phases. The total flux of active neutrinos was measured to be $5.54^{+0.33}_{-0.31}(stat.)^{+0.36}_{-0.34}(syst.)\times 10^{6}$ cm$^{-2}$ s$^{-1}$, consistent with previous measurements and standard solar models. A global analysis of solar and reactor neutrino mixing parameters yielded the best-fit values of $\Delta m^2 = 7.59^{+0.19}_{-0.21}\times 10^{-5}{eV}^2$ and $\theta = 34.4^{+1.3}_{-1.2}$ degrees.