Constraints on dark matter-nucleon effective couplings in the presence of kinematically distinct halo substructures using the DEAP-3600 detector

DEAP-3600 is a single-phase liquid argon detector aiming to directly detect Weakly Interacting Massive Particles (WIMPs), located at SNOLAB (Sudbury, Canada). After analyzing data taken during the first year of operation, a null result was used to place an upper bound on the WIMP-nucleon spin-independent, isoscalar cross section. This study reinterprets this result within a Non-Relativistic Effective Field Theory framework, and further examines how various possible substructures in the local dark matter halo may affect these constraints. Such substructures are hinted at by kinematic structures in the local stellar distribution observed by the Gaia satellite and other recent astronomical surveys. These include the Gaia Sausage (or Enceladus), as well as a number of distinct streams identified in recent studies. Limits are presented for the coupling strength of the effective contact interaction operators $\mathcal{O}_1$, $\mathcal{O}_3$, $\mathcal{O}_5$, $\mathcal{O}_8$, and $\mathcal{O}_{11}$, considering isoscalar, isovector, and xenonphobic scenarios, as well as the specific operators corresponding to millicharge, magnetic dipole, electric dipole, and anapole interactions. The effects of halo substructures on each of these operators are explored as well, showing that the $\mathcal{O}_5$ and $\mathcal{O}_8$ operators are particularly sensitive to the velocity distribution, even at dark matter masses above 100 GeV/$c^2$

The liquid-argon scintillation pulseshape in DEAP-3600

DEAP-3600 is a liquid-argon scintillation detector looking for dark matter. Scintillation events in the liquid argon (LAr) are registered by 255 photomultiplier tubes (PMTs), and pulseshape discrimination (PSD) is used to suppress electromagnetic background events. The excellent PSD performance of LAr makes it a viable target for dark matter searches, and the LAr scintillation pulseshape discussed here is the basis of PSD. The observed pulseshape is a combination of LAr scintillation physics with detector effects. We present a model for the pulseshape of electromagnetic background events in the energy region of interest for dark matter searches. The model is composed of (a) LAr scintillation physics, including the so-called intermediate component, (b) the time response of the TPB wavelength shifter, including delayed TPB emission at O(ms) time-scales, and c) PMT response. TPB is the wavelength shifter of choice in most LAr detectors. We find that approximately 10% of the intensity of the wavelength-shifted light is in a long-lived state of TPB. This causes light from an event to spill into subsequent events to an extent not usually accounted for in the design and data analysis of LAr-based detectors

Constraints on dark matter-nucleon effective couplings in the presence of kinematically distinct halo substructures using the DEAP-3600 detector

DEAP-3600 is a single-phase liquid argon detector aiming to directly detect weakly interacting massive particles (WIMPs), located at SNOLAB (Sudbury, Canada). After analyzing data taken during the first year of operation, a null result was used to place an upper bound on the WIMP-nucleon, spin-independent, isoscalar cross section. This study reinterprets this result within a nonrelativistic effective field theory framework and further examines how various possible substructures in the local dark matter halo may affect these constraints. Such substructures are hinted at by kinematic structures in the local stellar distribution observed by the Gaia satellite and other recent astronomical surveys. These include the Gaia Sausage (or Enceladus), as well as a number of distinct streams identified in recent studies. Limits are presented for the coupling strength of the effective contact interaction operators O1, O3, O5, O8, and O11, considering isoscalar, isovector, and xenonphobic scenarios, as well as the specific operators corresponding to millicharge, magnetic dipole, electric dipole, and anapole interactions. The effects of halo substructures on each of these operators are explored as well, showing that the O5 and O8 operators are particularly sensitive to the velocity distribution, even at dark matter masses above 100 GeV=c

Pulse-shape discrimination against low-energy Ar-39 beta decays in liquid argon with 4.5 tonne-years of DEAP-3600 data

The DEAP-3600 detector searches for the scintillation signal from dark matter particles scattering on a 3.3 tonne liquid argon target. The largest background comes from 39Ar beta decays and is suppressed using pulse-shape discrimination (PSD). We use two types of PSD estimator: the prompt-fraction, which considers the fraction of the scintillation signal in a narrow and a wide time window around the event peak, and the log-likelihood-ratio, which compares the observed photon arrival times to a signal and a background model. We furthermore use two algorithms to determine the number of photons detected at a given time: (1) simply dividing the charge of each PMT pulse by the mean single-photoelectron charge, and (2) a likelihood analysis that considers the probability to detect a certain number of photons at a given time, based on a model for the scintillation pulse shape and for afterpulsing in the light detectors. The prompt-fraction performs approximately as well as the log-likelihood-ratio PSD algorithm if the photon detection times are not biased by detector effects. We explain this result using a model for the information carried by scintillation photons as a function of the time when they are detected

SiPM-matrix readout of two-phase argon detectors using electroluminescence in the visible and near infrared range

Proportional electroluminescence (EL) in noble gases is used in two-phase detectors for dark matter searches to record (in the gas phase) the ionization signal induced by particle scattering in the liquid phase. The “standard” EL mechanism is considered to be due to noble gas excimer emission in the vacuum ultraviolet (VUV). In addition, there are two alternative mechanisms, producing light in the visible and near infrared (NIR) ranges. The first is due to bremsstrahlung of electrons scattered on neutral atoms (“neutral bremsstrahlung”, NBrS). The second, responsible for electron avalanche scintillation in the NIR at higher electric fields, is due to transitions between excited atomic states. In this work, we have for the first time demonstrated two alternative techniques of the optical readout of two-phase argon detectors, in the visible and NIR range, using a silicon photomultiplier matrix and electroluminescence due to either neutral bremsstrahlung or avalanche scintillation. The amplitude yield and position resolution were measured for these readout techniques, which allowed to assess the detection threshold for electron and nuclear recoils in two-phase argon detectors for dark matter searches. To the best of our knowledge, this is the first practical application of the NBrS effect in detection science

Design and construction of a new detector to measure ultra-low radioactive-isotope contamination of argon

Large liquid argon detectors offer one of the best avenues for the detection of galactic weakly interacting massive particles (WIMPs) via their scattering on atomic nuclei. The liquid argon target allows exquisite discrimination between nuclear and electron recoil signals via pulse-shape discrimination of the scintillation signals. Atmospheric argon (AAr), however, has a naturally occurring radioactive isotope, 39Ar, a β emitter of cosmogenic origin. For large detectors, the atmospheric 39Ar activity poses pile-up concerns. The use of argon extracted from underground wells, deprived of 39Ar, is key to the physics potential of these experiments. The DarkSide-20k dark matter search experiment will operate a dual-phase time projection chamber with 50 tonnes of radio-pure underground argon (UAr), that was shown to be depleted of 39Ar with respect to AAr by a factor larger than 1400. Assessing the 39Ar content of the UAr during extraction is crucial for the success of DarkSide-20k, as well as for future experiments of the Global Argon Dark Matter Collaboration (GADMC). This will be carried out by the DArT in ArDM experiment, a small chamber made with extremely radio-pure materials that will be placed at the centre of the ArDM detector, in the Canfranc Underground Laboratory (LSC) in Spain. The ArDM LAr volume acts as an active veto for background radioactivity, mostly γ-rays from the ArDM detector materials and the surrounding rock. This article describes the DArT in ArDM project, including the chamber design and construction, and reviews the background required to achieve the expected performance of the detector

First Direct Detection Constraints on Planck-Scale Mass Dark Matter with Multiple-Scatter Signatures Using the DEAP-3600 Detector

Dark matter with Planck-scale mass (?1019 GeV/c2) arises in well-motivated theories and could be produced by several cosmological mechanisms. A search for multiscatter signals from supermassive dark matter was performed with a blind analysis of data collected over a 813 d live time with DEAP-3600, a 3.3 t single-phase liquid argon-based detector at SNOLAB. No candidate signals were observed, leading to the first direct detection constraints on Planck-scale mass dark matter. Leading limits constrain dark matter masses between 8.3×106 and 1.2×1019 GeV/c2, and Ar40-scattering cross sections between 1.0×10-23 and 2.4×10-18 cm2. These results are interpreted as constraints on composite dark matter models with two different nucleon-to-nuclear cross section scalings

Erratum: Constraints on dark matter-nucleon effective couplings in the presence of kinematically distinct halo substructures using the DEAP-3600 detector (Physical Review D (2020) 102 (082001) DOI: 10.1103/PhysRevD.102.082001)

In the article, the Non-Relativistic Effective Field Theory (NREFT) rate calculations were determined using the wimpy_nreft software [1], which was updated on September 29, 2021, to include a previously missing (q/mN)2 factor in the implementation. This update affects the results related to the O3 operator that now scales as (q/mN)4 instead of (q/mN)2. The corrections to Figs. 2, 6, 9, 10, and 11 are presented below. The couplings to O3 constrained by this analysis are higher than those reported in the article. Additionally: (i) In Sec. V A, operator O3 is suppressed at low recoil energies, exhibiting now a peak around 50 keV (Fig. 2). (ii) The third paragraph in Sec. V B should read as follows: “The operator O3 is proportional to (q/mN)4, while O11 goes as (q/mN)2. O3 is described by the F'' multipole operator [discussed in Eqs. (9) and (10)], while O11 is described by M. Since the former operator is related to spin-orbit coupling, it couples to the two unpaired neutrons and proton holes in 40 Ar , rather than to all 40 nucleons. This leads to a suppression of ~10 2 in addition to the extra q2 suppression.” (iii) In Sec. V F, the statement “Operators that introduce a factor of q2 to the DM response function, such as O3, O5, and O11 change the shape of the recoil energy spectrum, compared to O1” should read “Operators that introduce a factor of q2 or q4 to the DM response function, such as O3, O5, and O11 change the shape of the recoil energy spectrum, compared to O1.” (iv) The sentence in Sec. VI “Constraints on operators proportional to v? are weaker than those proportional to q, which are weaker than those proportional to neither” should read “Constraints on operators proportional to vn? are weaker than those proportional to q raised to the same power, which in turn are weaker than constant couplings.” (v) In Sec. VI, exclusion curves on O3 and data to reproduce its recoil energy spectra were uploaded to a new Zenodo version [2]. (Figure Presented)

Pulse-shape discrimination against low-energy Ar-39 beta decays in liquid argon with 4.5 tonne-years of DEAP-3600 data

The DEAP-3600 detector searches for the scintillation signal from dark matter particles scattering on a 3.3 tonne liquid argon target. The largest background comes from 39Ar beta decays and is suppressed using pulse-shape discrimination (PSD). We use two types of PSD estimator: the prompt-fraction, which considers the fraction of the scintillation signal in a narrow and a wide time window around the event peak, and the log-likelihood-ratio, which compares the observed photon arrival times to a signal and a background model. We furthermore use two algorithms to determine the number of photons detected at a given time: (1) simply dividing the charge of each PMT pulse by the mean single-photoelectron charge, and (2) a likelihood analysis that considers the probability to detect a certain number of photons at a given time, based on a model for the scintillation pulse shape and for afterpulsing in the light detectors. The prompt-fraction performs approximately as well as the log-likelihood-ratio PSD algorithm if the photon detection times are not biased by detector effects. We explain this result using a model for the information carried by scintillation photons as a function of the time when they are detected