Silicon Detector Dark Matter Results from the Final Exposure of CDMS II

Extensive evidence indicates that a large fraction of the matter in the universe is nonluminous, nonbaryonic, and “cold” — nonrelativistic at the time matter began to dominate the energy density of the universe. Weakly Interacting Massive Particles (WIMPs) are an excellent candidate for nonbaryonic, cold dark matter. Minimal supersymmetry provides a natural WIMP candidate in the form of the lightest superpartner, with a typical mass Mδ ~ 100 GeV c-2 . WIMPs are expected to have collapsed into a roughly isothermal, spherical halo within which the visible portion of our galaxy resides. They would scatter off nuclei via the weak interaction, potentially allowingtheir direct detection. The Cryogenic Dark Matter Search (CDMS) employs Ge and Si detectors to search for WIMPs via their elastic-scatteringinteractions with nuclei while discriminatingagainst interactions of background particles. The former yield nuclear recoils while the latter produce electron recoils. The ionization yield (the ratio of ionization production to recoil energy in a semiconductor) of a particle interaction differs greatly for nuclear and electron recoils. CDMS detectors measure phonon and electron-hole-pair production to determine recoil energy and ionization yield for each event and thereby discriminate nuclear recoils from electron recoils. This dissertation reports new limits on the spin-independent WIMP-nucleon elastic-scattering cross section that exclude unexplored parameter space above 10 GeV c-2 WIMP mass and, at > 75% CL, the entire 3σ allowed region for the WIMP signal reported by the DAMA experiment. The experimental apparatus, detector performance, and data analysis are fully described.

We discuss a novel dark matter signature relevant for directional detection of Weakly Interacting Massive Particles (WIMPs). For heavy enough WIMPs and low enough recoil energies, the maximum of the recoil rate is not in the direction of the average WIMP arrival direction but in a ring around it at an angular radius that increases with the WIMP mass and can approach 90° at very low energies. The ring is easier to detect for smaller WIMP velocity dispersion and larger average WIMP velocities relative to the detector. In principle the ring could be used as an additional indication of the WIMP mass range.

While evidence of non-baryonic dark matter has been accumulating for decades, its exact nature continues to remain a mystery. Weakly Interacting Massive Particles (WIMPs) are a well motivated candidate which appear in certain extensions of the Standard Model, independently of dark matter theory. If such particles exist, they should occasionally interact with particles of normal matter, producing a signal which may be detected. The DarkSide-50 direct dark matter experiment aims to detect the energy of recoiling argon atoms due to the elastic scattering of postulated WIMPs. In order to make such a discovery, a clear understanding of both the background and signal region is essential. This understanding requires a careful study of the detector's response to radioactive sources, which in turn requires such sources may be safely introduced into or near the detector volume and reliably removed. The CALibration Insertaion System (CALIS) was designed and built for this purpose in a j oint effort between Fermi National Laboratory and the University of Hawaii. This work describes the design and testing of CALIS, its installation and commissioning at the Laboratori Nazionali del Gran Sasso (LNGS) and the multiple calibration campaigns which have successfully employed it. As nuclear recoils produced by WIMPs are indistinguishable from those produced by neutrons, radiogenic neutrons are both the most dangerous class of background and a vital calibration source for the study of the potential WIMP signal. Prior to the calibration of DarkSide-50 with radioactive neutron sources, the acceptance region was determined by the extrapolation of nuclear recoil data from a separate, dedicated experiment, ScENE, which measured the distribution of the pulse shape discrimination parameter, $f_{90}$, for nuclear recoils of known energies. This work demonstrates the validity of the extrapolation of ScENE values to DarkSide-50, by direct comparison of the $f_{90}$ distributio n of nuclear recoils from ScENE and an AmBe calibration sour! ce. The combined acceptance as defined by ScENE and the \textit{in-situ} AmBe calibration were used to establish the best WIMP exclusion limit on an argon target. Unfortunately, radioactive sources used for the calibration of DarkSide-50 are universally accompanied by gamma decays, which obscure the low energy region where most WIMP interactions are expected to occur and seem to make continuing dependence on an external measurement such as ScENE inevitable. However, this work presents a novel method of nuclear recoil calibration employing event selection, unique to the design of DarkSide-50, which produces a nearly pure sample of nuclear recoils. Further, it describes the execution of a neutron calibration campaign, from planning to analysis, which yielded a valuable data set for defining the acceptance region. Together with the event selection techniques, this allows for the definition of the acceptance region independent of ScENE values. Two analytical models of the $f_{90 }$ distribution are described and their results for nuclear recoils are compared. Finally, a detailed study of integrated noise in nuclear and electron recoil

A wide variety of astrophysical observations indicate that approximately 85% of the in the universe is nonbaryonic and nonluminous. Understanding the nature of this matter is one of the most important outstanding questions in cosmology. Weakly Interacting Massive Particles (WIMPs) are a leading candidate for dark since they would be thermally produced in the early universe in the correct abundance to account for the observed relic density of dark matter. If WIMPs account for the dark matter, then rare interactions from relic WIMPs should be observable in terrestrial detectors. Recently, unexplained excess events in the DAMA/LIBRA, CoGeNT, and CRESST-II experiments have been interpreted as evidence of scattering from WIMPs with masses ~10 GeV and spin-independent scattering cross sections of 10-41-10-40 cm2. The Cryogenic Dark Matter Search (CDMS II) attempts to identify WIMP interactions using an array of cryogenic germanium and silicon particle detectors located at the Soudan Underground Laboratory in northern Minnesota. In this dissertation, data taken by CDMS II are reanalyzed using a 2 keV recoil energy threshold to increase the sensitivity to WIMPs with masses ~10 GeV. These data disfavor an explanation for the DAMA/LIBRA, CoGeNT, and CRESST-II results in terms of spin-independent elastic scattering of WIMPs with masses ≲12 GeV, under standard assumptions. At the time of publication, they provided the strongest constraints on spin-independent elastic scattering from 5-9 GeV, ruling out previously unexplored parameter space. To detect WIMPs or exclude the remaining parameter space favored by the most popular models will ultimately require detectors with target masses ≳1 ton, requiring an increase in mass by more than two orders of magnitude over CDMS II. For cryogenic detectors such as CDMS, scaling to such large target masses will require individual detector elements to be fabricated more quickly and cheaply, while maintaining the nearly background-free operation of the existing experiment. We describe the development of athermal phonon mediated particle detectors using Microwave Kinetic Inductance Detectors (MKIDs), which could provide a simpler path to extending the CDMS detector technology to the ton scale. Results from prototype devices have demonstrated energy resolutions as good as σ = 0.55 keV at 30 keV, comparable to existing CDMS II detectors. Such designs can be scaled to kg-scale detector elements, while reducing the complexity of the detector fabrication and cryogenic readout electronics relative to existing designs. Since MKIDs are naturally multiplexed in the frequency domain, MKID-based designs also allow much finer pixelization of the phonon sensor, which is expected to enhance background rejection for large detectors while simultaneously reducing the number of wires needed to read out the detectors.

A wide variety of astrophysical observations indicate that approximately 85% of the matter in the universe is nonbaryonic and nonluminous. Understanding the nature of this "dark matter" is one of the most important outstanding questions in cosmology. Weakly Interacting Massive Particles (WIMPs) are a leading candidate for dark matter since they would be thermally produced in the early universe in the correct abundance to account for the observed relic density of dark matter. If WIMPs account for the dark matter, then rare interactions from relic WIMPs should be observable in terrestrial detectors. Recently, unexplained excess events in the DAMA/LIBRA, CoGeNT, and CRESST-II experiments have been interpreted as evidence of scattering from WIMPs with masses ~10 GeV and spin-independent scattering cross sections of 10-41-10-40 cm2. The Cryogenic Dark Matter Search (CDMS II) attempts to identify WIMP interactions using an array of cryogenic germanium and silicon particle detectors located at the Soudan Underground Laboratory in northern Minnesota. In this dissertation, data taken by CDMS II are reanalyzed using a 2 keV recoil energy threshold to increase the sensitivity to WIMPs with masses ~10 GeV. These data disfavor an explanation for the DAMA/LIBRA, CoGeNT, and CRESST-II results in terms of spin-independent elastic scattering of WIMPs with masses ≲12 GeV, under standard assumptions. At the time of publication, they provided the strongest constraints on spin-independent elastic scattering from 5-9 GeV, ruling out previously unexplored parameter space. To detect WIMPs or exclude the remaining parameter space favored by the most popular models will ultimately require detectors with target masses ≳1 ton, requiring an increase in mass by more than two orders of magnitude over CDMS II. For cryogenic detectors such as CDMS, scaling to such large target masses will require individual detector elements to be fabricated more quickly and cheaply, while maintaining the nearly background-free operation of the existing experiment. We describe the development of athermal phonon mediated particle detectors using Microwave Kinetic Inductance Detectors (MKIDs), which could provide a simpler path to extending the CDMS detector technology to the ton scale. Results from prototype devices have demonstrated energy resolutions as good as σ = 0.55 keV at 30 keV, comparable to existing CDMS II detectors. Such designs can be scaled to kg-scale detector elements, while reducing the complexity of the detector fabrication and cryogenic readout electronics relative to existing designs. Since MKIDs are naturally multiplexed in the frequency domain, MKID-based designs also allow much finer pixelization of the phonon sensor, which is expected to enhance background rejection for large detectors while simultaneously reducing the number of wires needed to read out the detectors.

Although dark matter appears to constitute over 80% of the matter in the Universe, its composition is a mystery. Astrophysical observations suggest that the luminous portions of the Galaxy are embedded in a halo of darkmatter particles. Weakly Interacting Massive Particles (WIMPs) are the most studied class of dark-matter candidates and arise naturally within the context of many weak-scale supersymmetric theories. Direct-detection experiments like the Cryogenic Dark Matter Search (CDMS) strive to discern the kinetic energy of recoiling nuclei resulting from WIMP interactions with terrestrial matter. This is a considerable challenge in which the low (expected) rate of WIMP interactions must be distinguished from an overwhelming rate due to known types of radiation. An incontrovertible positive detection has remained elusive. However, a few experiments have recorded data that appear consistent with a low-mass WIMP. This thesis describes an attempt to probe the favored parameter space. To increase sensitivity to low-mass WIMPs, a low-threshold technique with improved sensitivity to small energy depositions is applied to CDMS shallowsite data. Four germanium and two silicon detectors were operated between December 2001 and June 2002, yielding 118 days of exposure. By sacrificing some of the CDMS detectors’ ability to discriminate signal from background, energy thresholds of ~1 and ~2 keV were achieved for three of the germanium and both silicon detectors, respectively. A large number of WIMP candidate events are observed, most of which can be accounted for by misidentification of background sources. No conclusive evidence for a low-mass WIMP signal is found. The observed event rates are used to set upper limits on the WIMPnucleon scattering cross section as a function of WIMP mass. Interesting parameter space is excluded for WIMPs with masses below ~9GeV/c2. Under standard assumptions, the parameter space favored by interpretations of other experiments’ data as low-mass WIMP signals is partially excluded, and new parameter space is excluded for WIMP masses between 3 and 4GeV/c2.

We provide a first systematic and quantitative discussion of the phenomenology of the nonrelativistic effective Hamiltonian describing the nuclear scattering process for a weakly interacting massive particle (WIMP) of arbitrary spin ${j}_{\ensuremath{\chi}}$. To this aim we obtain constraints from a representative sample of present direct detection experiments assuming the WIMP-nucleus scattering process to be driven by each one of the 44 effective couplings that arise for ${j}_{\ensuremath{\chi}}\ensuremath{\le}2$. We find that a high value of the multipolarity $s\ensuremath{\le}2{j}_{\ensuremath{\chi}}$ of the coupling, related to the power of the momentum transfer $q$ appearing in the scattering amplitude, leads to a suppression of the expected rates and pushes the expected differential spectra to large recoil energies ${E}_{R}$. For $s\ensuremath{\le}4$ the effective scales probed by direct detection experiments can be suppressed by up to five orders of magnitude compared to the case of a standard spin-independent interaction. For operators with large $s$ the expected differential spectra can be pushed to recoil energies in the MeV range, with the largest part of the signal concentrated at ${E}_{R}\ensuremath{\gtrsim}100\text{ }\text{ }\mathrm{keV}$ and a peculiar structure of peaks and minima arising when both the nuclear target and the WIMPs are heavy. As a consequence the present bounds on the effective operators can be significantly improved by extending the recoil energy intervals to higher recoil energies. Our analysis assumes effective interaction operators that are irreducible under the rotation group. Such operators drive the interactions of high-multipole dark matter candidates, i.e., states that possess only the highest multipole allowed by their spin. As a consequence our analysis represents also the first phenomenological study of the direct detection of quadrupolar, octupolar, and hexadecapolar dark matter.

Direct detection of Weakly Interacting Massive Particle (WIMP) candidates of Dark Matter (DM) is studied within the context of a self-consistent truncated isothermal model of the finite-size dark halo of the Galaxy. The halo model, based on the ``King model'' of the phase space distribution function of collisionless DM particles, takes into account the modifications of the phase-space structure of the halo due to the gravitational influence of the observed visible matter in a self-consistent manner. The parameters of the halo model are determined by a fit to a recently determined circular rotation curve of the Galaxy that extends up to ∼ 60 kpc. Unlike in the Standard Halo Model (SHM) customarily used in the analysis of the results of WIMP direct detection experiments, the velocity distribution of the WIMPs in our model is non-Maxwellian with a cut-off at a maximum velocity that is self-consistently determined by the model itself. For our halo model that provides the best fit to the rotation curve data, the 90% C.L. upper limit on the WIMP-nucleon spin-independent cross section from the recent results of the CDMS-II experiment, for example, is ∼ 5.3 × 10−8 pb at a WIMP mass of ∼ 71 GeV. We also find, using the original 2-bin annual modulation amplitude data on the nuclear recoil event rate seen in the DAMA experiment, that there exists a range of small WIMP masses, typically ∼ 2–16 GeV, within which DAMA collaboration's claimed annual modulation signal purportedly due to WIMPs is compatible with the null results of other experiments. These results, based as they are on a self-consistent model of the dark matter halo of the Galaxy, strengthen the possibility of low-mass (≲10 GeV) WIMPs as a candidate for dark matter as indicated by several earlier studies performed within the context of the SHM. A more rigorous analysis using DAMA bins over smaller intervals should be able to better constrain the ``DAMA regions'' in the WIMP parameter space within the context of our model.

Direction sensitive direct detection of Weakly Interacting Massive Particles (WIMPs) as dark matter would provide an unambiguous non-gravitational signature of dark matter (DM). The diurnal variation of DM signal due to earth's rotation around its own axis can be a significant signature for galactic WIMPs. Because of particular orientation of earth's axis of rotation with respect to WIMP wind direction, the apparent direction of WIMP wind as observed at a detector can alter widely over a day. In this work we calculate the directional detection rates with their daily and yearly modulations in earth-bound dark matter experiments considering detailed features of the geometry and dynamics of the earth-sun system along with the solar motion in galactic frame. A separate halo model namely the dark disc model other than the usual standard halo model for dark matter halo is also considered and the results for two models are compared. We demonstrate the results for two types of gas detectors namely DRIFT (target material CS2) and NEWAGE (target material CF4) that use Time Projection Chamber techniques for measuring directionality of the recoil nucleus. The WIMP mass and recoil energy dependence of the daily variation of event rates are computed for specific detector and the sensitive ranges of mass and recoil energies for the considered detector are probed.

The CDEX-1 experiment conducted a search of low-mass (< 10 GeV/c2) Weakly Interacting Massive Particles (WIMPs) dark matter at the China Jinping Underground Laboratory using a p-type point-contact germanium detector with a fiducial mass of 915 g at a physics analysis threshold of 475 eVee. We report the hardware set-up, detector characterization, data acquisition and analysis procedures of this experiment. No excess of unidentified events are observed after subtraction of known background. Using 335.6 kg-days of data, exclusion constraints on the WIMP-nucleon spin-independent and spin-dependent couplings are derived.

The XENONnT experiment searches for weakly interacting massive particle (WIMP) dark matter scattering off a xenon nucleus. In particular, XENONnT uses a dual-phase time projection chamber with a 5.9-ton liquid xenon target, detecting both scintillation and ionization signals to reconstruct the energy, position, and type of recoil. A blind search for nuclear recoil WIMPs with an exposure of 1.1 ton-years (4.18 t fiducial mass) yielded no signal excess over background expectations, from which competitive exclusion limits were derived on WIMP-nucleon elastic scatter cross sections, for WIMP masses ranging from 6 GeV / c 2 up to the TeV / c 2 scale. This work details the modeling and statistical methods employed in this search. By means of calibration data, we model the detector response, which is then used to derive background and signal models. The construction and validation of these models is discussed, alongside additional purely data-driven backgrounds. We also describe the statistical inference framework, including the definition of the likelihood function and the construction of confidence intervals.

The next generation of dark matter (DM) direct detection experiments and neutrino telescopes will probe large swaths of dark matter parameter space. In order to interpret the signals in these experiments, it is necessary to have good models of both the halo DM streaming through the Solar System and the population of DM bound to the Solar System. In this paper, the first in a series of three on DM in the Solar System, we present simulations of orbits of DM bound to the Solar System by solar capture in a toy solar system consisting of only the Sun and Jupiter, assuming that DM consists of a single species of weakly interacting massive particle (WIMP). We describe how the size of the bound WIMP population depends on the WIMP mass m_chi, spin-independent cross section sigma_p^(SI), and spin-dependent cross section sigma_p^(SD). Using a standard description of the Galactic DM halo, we find that the maximum enhancement to the direct detection event rate, consistent with current experimental constraints on the WIMP-nucleon cross section, is <1% relative to the event rate from halo WIMPs, while the event rate from neutrinos from WIMP annihilation in the center of the Earth is unlikely to meet the threshold of next-generation, km^3-sized (IceCube, KM3NeT) neutrino telescopes.

The count rate of weakly interacting massive particle (WIMP) dark matter candidates in direct detection experiments experiences an annual modulation due to the Earth's motion around the Sun. In the standard isothermal halo model, the signal peaks near June 2nd at high recoil energies; however, the signal experiences a phase reversal and peaks in December at low energy recoils. We show that this phase reversal may be used to determine the WIMP mass. If an annual modulation were observed with the usual phase (i.e., peaking on June 2nd) in the lowest accessible energy recoil bins of the DAMA, CDMS-II, CRESST-II, EDELWEISS-II, GENIUS-TF, ZEPLIN-II, XENON, or ZEPLIN-IV detectors, one could immediately place upper bounds on the WIMP mass of 103, 48, 6, 97, 10, 52, 29, and 29 GeV, respectively. In addition, detectors with adequate energy resolution and sufficiently low recoil energy thresholds may determine the crossover recoil energy at which the phase reverses, thereby obtaining an independent measurement of the WIMP mass. We study the capabilities of various detectors, and find that CRESST-II, ZEPLIN-II, and GENIUS-TF should be able to observe the phase reversal in a few years of runtime, and can thus determine the mass of the WIMP if it is $O(100\mathrm{GeV}).$ Xenon based detectors with 1000 kg (XENON and ZEPLIN-IV) and with energy recoil thresholds of a few keV require 25 kg yr exposure, which will be readily attained in upcoming experiments.

Competitive limits on the weakly interacting massive particle (WIMP) spin-independent scattering cross section are currently being produced by 76Ge detectors originally designed to search for neutrinoless double beta decay, such as the Heidelberg-Moscow and IGEX experiments. In the absence of background subtraction, limits on the WIMP interaction cross section are set by calculating the upper confidence limit on the theoretical event rate, given the observed event rate. The standard analysis technique involves calculating the 90% upper confidence limit on the number of events in each bin, and excluding any set of parameters (WIMP mass and cross-section) which produces a theoretical event rate for any bin which exceeds the 90% upper confidence limit on the event rate for that bin. We show that, if there is more than one energy bin, this produces exclusion limits that are actually at a lower degree of confidence than 90%, and are hence erroneously tight. We formulate criteria which produce true 90% confidence exclusion limits in these circumstances, including calculating the individual bin confidence limit for which the overall probability that no bins exceeds this confidence limit is 90% and calculating the 90% minimum confidence limit on the number of bins which exceed their individual bin 90% confidence limits. We then compare the limits on the WIMP cross-section produced by these criteria with those found using the standard technique, using data from the Heidelberg-Moscow and IGEX experiments.

Weakly Interacting Massive Particles (WIMPs) are one of the leading candidates for Dark Matter. So far the usual procedure for constraining the WIMP-nucleon cross sections in direct Dark Matter detection experiments have been to fit the predicted event rate based on some model(s) of the Galactic halo and of WIMPs to experimental data. One has to assume whether the spin-independent (SI) or the spin-dependent (SD) WIMP-nucleus interaction dominates, and results of such data analyses are also expressed as functions of the as yet unknown WIMP mass. In this article, I introduce methods for extracting information on the WIMP-nucleon cross sections by considering a general combination of the SI and SD interactions. Neither prior knowledge about the local density and the velocity distribution of halo WIMPs nor about their mass is needed. Assuming that an exponential-like shape of the recoil spectrum is confirmed from experimental data, the required information are only the measured recoil energies (in low energy ranges) and the number of events in the first energy bin from two or more experiments.