Recommended conventions for reporting results from direct dark matter searches

The field of dark matter detection is a highly visible and highly competitive one. In this paper, we propose recommendations for presenting dark matter direct detection results particularly suited for weak-scale dark matter searches, although we believe the spirit of the recommendations can apply more broadly to searches for other dark matter candidates, such as very light dark matter or axions. To translate experimental data into a final published result, direct detection collaborations must make a series of choices in their analysis, ranging from how to model astrophysical parameters to how to make statistical inferences based on observed data. While many collaborations follow a standard set of recommendations in some areas, for example the expected flux of dark matter particles (to a large degree based on a paper from Lewin and Smith in 1995), in other areas, particularly in statistical inference, they have taken different approaches, often from result to result by the same collaboration. We set out a number of recommendations on how to apply the now commonly used Profile Likelihood Ratio method to direct detection data. In addition, updated recommendations for the Standard Halo Model astrophysical parameters and relevant neutrino fluxes are provided. The authors of this note include members of the DAMIC, DarkSide, DARWIN, DEAP, LZ, NEWS-G, PandaX, PICO, SBC, SENSEI, SuperCDMS, and XENON collaborations, and these collaborations provided input to the recommendations laid out here. Wide-spread adoption of these recommendations will make it easier to compare and combine future dark matter results

Simulations of events for the LUX-ZEPLIN (LZ) dark matter experiment

The LUX-ZEPLIN dark matter search aims to achieve a sensitivity to the WIMP-nucleon spin-independent cross-section down to (1–2)×10−12 pb at a WIMP mass of 40 GeV/c2. This paper describes the simulations framework that, along with radioactivity measurements, was used to support this projection, and also to provide mock data for validating reconstruction and analysis software. Of particular note are the event generators, which allow us to model the background radiation, and the detector response physics used in the production of raw signals, which can be converted into digitized waveforms similar to data from the operational detector. Inclusion of the detector response allows us to process simulated data using the same analysis routines as developed to process the experimental data

Projected WIMP sensitivity of the LUX-ZEPLIN dark matter experiment

LUX-ZEPLIN (LZ) is a next-generation dark matter direct detection experiment that will operate 4850 feet underground at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, USA. Using a two-phase xenon detector with an active mass of 7 tonnes, LZ will search primarily for low-energy interactions with weakly interacting massive particles (WIMPs), which are hypothesized to make up the dark matter in our galactic halo. In this paper, the projected WIMP sensitivity of LZ is presented based on the latest background estimates and simulations of the detector. For a 1000 live day run using a 5.6-tonne fiducial mass, LZ is projected to exclude at 90% confidence level spin-independent WIMP-nucleon cross sections above 1.4×10-48 cm2 for a 40 GeV/c2 mass WIMP. Additionally, a 5σ discovery potential is projected, reaching cross sections below the exclusion limits of recent experiments. For spin-dependent WIMP-neutron(-proton) scattering, a sensitivity of 2.3×10-43 cm2 (7.1×10-42 cm2) for a 40 GeV/c2 mass WIMP is expected. With underground installation well underway, LZ is on track for commissioning at SURF in 2020

Projected sensitivity of the LUX-ZEPLIN experiment to the 0νββ decay of Xe 136

The LUX-ZEPLIN (LZ) experiment will enable a neutrinoless double β decay search in parallel to the main science goal of discovering dark matter particle interactions. We report the expected LZ sensitivity to Xe136 neutrinoless double β decay, taking advantage of the significant (>600 kg) Xe136 mass contained within the active volume of LZ without isotopic enrichment. After 1000 live-days, the median exclusion sensitivity to the half-life of Xe136 is projected to be 1.06×1026 years (90% confidence level), similar to existing constraints. We also report the expected sensitivity of a possible subsequent dedicated exposure using 90% enrichment with Xe136 at 1.06×1027 years

Improved Limits on Scattering of Weakly Interacting Massive Particles from Reanalysis of 2013 LUX Data

We present constraints on weakly interacting massive particles (WIMP)-nucleus scattering from the 2013 data of the Large Underground Xenon dark matter experiment, including 1.4×104 kg day of search exposure. This new analysis incorporates several advances: single-photon calibration at the scintillation wavelength, improved event-reconstruction algorithms, a revised background model including events originating on the detector walls in an enlarged fiducial volume, and new calibrations from decays of an injected tritium β source and from kinematically constrained nuclear recoils down to 1.1 keV. Sensitivity, especially to low-mass WIMPs, is enhanced compared to our previous results which modeled the signal only above a 3 keV minimum energy. Under standard dark matter halo assumptions and in the mass range above 4 GeV c-2, these new results give the most stringent direct limits on the spin-independent WIMP-nucleon cross section. The 90% C.L. upper limit has a minimum of 0.6 zb at 33 GeV c-2 WIMP mass

Results on the Spin-Dependent Scattering of Weakly Interacting Massive Particles on Nucleons from the Run 3 Data of the LUX Experiment

We present experimental constraints on the spin-dependent WIMP (weakly interacting massive particle)-nucleon elastic cross sections from LUX data acquired in 2013. LUX is a dual-phase xenon time projection chamber operating at the Sanford Underground Research Facility (Lead, South Dakota), which is designed to observe the recoil signature of galactic WIMPs scattering from xenon nuclei. A profile likelihood ratio analysis of 1.4×104 kg day of fiducial exposure allows 90% C.L. upper limits to be set on the WIMP-neutron (WIMP-proton) cross section of σn=9.4×10-41 cm2 (σp=2.9×10-39 cm2) at 33 GeV/c2. The spin-dependent WIMP-neutron limit is the most sensitive constraint to date

Chromatographic separation of radioactive noble gases from xenon

The Large Underground Xenon (LUX) experiment operates at the Sanford Underground Research Facility to detect nuclear recoils from the hypothetical Weakly Interacting Massive Particles (WIMPs) on a liquid xenon target. Liquid xenon typically contains trace amounts of the noble radioactive isotopes 85 Kr and 39 Ar that are not removed by the in situ gas purification system. The decays of these isotopes at concentrations typical of research-grade xenon would be a dominant background for a WIMP search experiment. To remove these impurities from the liquid xenon, a chromatographic separation system based on adsorption on activated charcoal was built. 400 kg of xenon was processed, reducing the average concentration of krypton from 130 ppb to 3.5 ppt as measured by a cold-trap assisted mass spectroscopy system. A 50 kg batch spiked to 0.001 g/g of krypton was processed twice and reduced to an upper limit of 0.2 ppt

Measurement of the gamma ray background in the Davis cavern at the Sanford Underground Research Facility

Deep underground environments are ideal for low background searches due to the attenuation of cosmic rays by passage through the earth. However, they are affected by backgrounds from γ-rays emitted by 40K and the 238U and 232Th decay chains in the surrounding rock. The LUX-ZEPLIN (LZ) experiment will search for dark matter particle interactions with a liquid xenon TPC located within the Davis campus at the Sanford Underground Research Facility, Lead, South Dakota, at the 4850-foot level. In order to characterise the cavern background, in-situ γ-ray measurements were taken with a sodium iodide detector in various locations and with lead shielding. The integral count rates (0–3300 keV) varied from 596 Hz to 1355 Hz for unshielded measurements, corresponding to a total flux from the cavern walls of 1.9 ± 0.4 γ cm−2s−1. The resulting activity in the walls of the cavern can be characterised as 220 ± 60 Bq/kg of 40K, 29 ± 15 Bq/kg of 238U, and 13 ± 3 Bq/kg of 232Th

Cosmogenic production of {37}^Ar in the context of the LUX-ZEPLIN experiment

We estimate the amount of {37}^Ar produced in natural xenon via cosmic-ray-induced spallation, an inevitable consequence of the transportation and storage of xenon on the Earth’s surface. We then calculate the resulting {37}^Ar concentration in a 10-tonne payload (similar to that of the LUX-ZEPLIN experiment) assuming a representative schedule of xenon purification, storage, and delivery to the underground facility. Using the spallation model by Silberberg and Tsao, the sea-level production rate of {37}^Ar in natural xenon is estimated to be 0.024 atoms/kg/day. Assuming the xenon is successively purified to remove radioactive contaminants in 1-tonne batches at a rate of 1 tonne/month, the average {37}^Ar activity after 10 tons are purified and transported underground is 0.058 - 0.090 μ Bq/kg, depending on the degree of argon removal during above-ground purification. Such cosmogenic {37}^Ar will appear as a noticeable background in the early science data, while decaying with a 35-day half-life. This newly noticed production mechanism of {37}^Ar should be considered when planning for future liquid-xenon-based experiments

Simulations of events for the LUX-ZEPLIN (LZ) dark matter experiment

The LUX-ZEPLIN dark matter search aims to achieve a sensitivity to the WIMP-nucleon spin-independent cross-section down to (1–2)x10-12 pb at a WIMP mass of 40  GeV/c2. This paper describes the simulations framework that, along with radioactivity measurements, was used to support this projection, and also to provide mock data for validating reconstruction and analysis software. Of particular note are the event generators, which allow us to model the background radiation, and the detector response physics used in the production of raw signals, which can be converted into digitized waveforms similar to data from the operational detector. Inclusion of the detector response allows us to process simulated data using the same analysis routines as developed to process the experimental data