Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “detector simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,045 records · Page 58

Superconducting Qubits for Dark Matter Detection

Developments over the last decade have pushed the search for particle dark matter (DM) to new frontiers, including the keV-scale lower mass limit for thermally-produced DM. Galactic DM at this mass is kinematically matched with the energy needed to break a Cooper pair in common superconductors (~meV). Quantum sensors such as superconducting qubits are sensitive to these broken Cooper pairs, and can potentially be exploited as low-threshold detectors for particle-like DM scattering. The Quantum Science Center group at Fermilab exploring the use of qubit-based detectors for particle detection in the LOUD surface dilution fridge facility. This talk will discuss recent R&D efforts to understand qubit response to energy depositions through a combination of measurements in LOUD and low-energy G4CMP simulations of our superconducting qubit chips.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Search for Long-Lived Particles with Muon Detector Shower Signature in the CMS Run-3 data

Many beyond standard model theories predict the existence of long-lived particles (LLPs). These LLPs can have sizable lifetimes and decay several meters from their production vertex. In this poster/talk, we present the analysis strategy for searching for LLPs using the Compact Muon Solenoid (CMS) Experiment. The proton-proton collision data used in this analysis were collected from 2022 to 2024 at a center-of-mass energy of 13.6 TeV, corresponding to an integrated luminosity of 170 fb^-1. The LLP decays are reconstructed as a high-multiplicity cluster of detector hits in the cathode strip chambers (CSC) of the muon system endcap. This signature is referred to as the Muon Detector Showers (MDS). This search requires events to contain at least one MDS cluster. The analysis focuses on LLP hadronic decays and LLP masses up to a few tens of GeV. We present the signal properties in MonteCarlo simulation, event selection, background modeling, and evaluation of the expected sensitivity. The results are interpreted under the Twin Higgs model benchmark.

Agyemang-Duah, Andrews [Grambling State U.]↗

Studies on the scintillation light detection in the ProtoDUNE Dual Phase liquid-argon TPC and its capability for the supernova trigger in DUNE

The Deep Underground Neutrino Experiment (DUNE) is a long-baseline neutrino oscillation experiment that aims at addressing key questions in neutrino physics in the next decades. Its scientific program includes the detection of the neutrino flux from a core-collapse supernova. The DUNE far detector will have four 17-kt mass liquid-argon (LAr) time-projection chamber (TPC) modules. ProtoDUNE Dual Phase (DP), a dual-phase LAr TPC with 300 t of active mass and 6 m of drift distance, was operated with cosmic muons in 2019-2020 as part of an R&D program at the CERN Neutrino Platform to demonstrate the feasibility of the technology at such a large scale. In a LAr TPC, the photon detection system (PDS) provides fun- damental timing information and trigger capabilities. The PDS of ProtoDUNE-DP, which consisted of 36 photomultiplier tubes (PMTs), counted on a dedicated light calibration system (LCS) to monitor the PMT response. In this dissertation, the characterization and validation of the ProtoDUNE-DP PDS and LCS components before their installation will be reviewed, highlighting the results of general interest for experiments that use liquid noble gasses as target medium. The results from the stable performance of both systems in the detector during 15 months will be presented next as well as the studies on the scintillation light detection in ProtoDUNE-DP, where the collection of light produced in LAr at 7 m from the photosensors has been achieved for the first time. It is worth pointing out that the excellent LAr purity and the large size of the detector have enabled to develop a unique data-driven investigation on aspects that are critical for LAr- based experiments but that are not completely understood. The analyses cover the characterization of the low-energy background detected by the PDS, the quantification of the electric field impact on the light yield, the evaluation of the Rayleigh scattering affecting the light propagation, and the analysis of the PMT detection efficiency. The effect of the VUV reflectivity of the detector materials will be also discussed. In addition, the estimation of the cosmic muon flux crossing the TPC and the study of the observed light yield by the PDS will be reported. Finally, the results from the simulation-based study of the supernova burst trigger capability with the PDS of a 12.1-kt active mass dual-phase LAr TPC as the one proposed for DUNE will be summarized. Several configurations of reflective foils installed in the TPC to enhance the light collection will be compared

Gallego-Ros, Ana↗

Modeling the Effect of the Interplanetary Magnetic Field on Cosmic Ray Muon Shadows

Cosmic ray (CR) muon shadows were located using muon data from the MINOS Far Detector in Northern Minnesota. The shadows were located independently across 3 time periods; near solar minimum, solar maximum, and over the entire 13 year span of the data. A muon distribution for each shadow was then sampled to simulate CR motions near the Sun using the Parker Spiral model of the Sun’s IMF and a dipole model of the Earth’s GMF. The resulting particle distributions were then compared to their position with respect to the Sun. An analytical relationship between the primary CR cutoff energy E0 and the associated shadow deflections was observed. The minimum cutoff energy of the primary CR energy spectrum consistent with the Parker Spiral model was determined to be $E_0 = 22.8 \pm 1.2$~TeV. This value is 2-3 times higher than the experimentally determined value. Complexities that can be added to the Parker Spiral Model and affect the overall shift of the the Sun’s cosmic ray shadow are discuss ed to explain discrepancies. Results show that the Parker Spiral model is most consistent with the the solar minimum shadow and least consistent with the solar maximum shadow, as expected.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Studies on the scintillation light detection in the ProtoDUNE Dual Phase liquid-argon TPC and its capability for the supernova trigger in DUNE

The Deep Underground Neutrino Experiment (DUNE) is a long-baseline neutrino oscillation experiment that aims at addressing key questions in neutrino physics in the next decades. Its scientific program includes the detection of the neutrino flux from a core-collapse supernova. The DUNE far detector will have four 17-kt mass liquid-argon (LAr) time-projection chamber (TPC) modules. ProtoDUNE Dual Phase (DP), a dual-phase LAr TPC with 300 t of active mass and 6 m of drift distance, was operated with cosmic muons in 2019-2020 as part of an R&D program at the CERN Neutrino Platform to demonstrate the feasibility of the technology at such a large scale. In a LAr TPC, the photon detection system (PDS) provides fundamental timing information and trigger capabilities. The PDS of ProtoDUNE-DP, which consisted of 36 photomultiplier tubes (PMTs), counted on a dedicated light calibration system (LCS) to monitor the PMT response. In this dissertation, the characterization and validation of the ProtoDUNE-DP PDS and LCS components before their installation will be reviewed, highlighting the results of general interest for experiments that use liquid noble gasses as target medium. The results from the stable performance of both systems in the detector during 15 months will be presented next as well as the studies on the scintillation light detection in ProtoDUNE-DP, where the collection of light produced in LAr at 7 m from the photosensors has been achieved for the first time. It is worth pointing out that the excellent LAr purity and the large size of the detector have enabled to develop a unique data-driven investigation on aspects that are critical for LAr-based experiments but that are not completely understood. The analyses cover the characterization of the low-energy background detected by the PDS, the quantification of the electric field impact on the light yield, the evaluation of the Rayleigh scattering affecting the light propagation, and the analysis of the PMT detection efficiency. The effect of the VUV reflectivity of the detector materials will be also discussed. In addition, the estimation of the cosmic muon flux crossing the TPC and the study of the observed light yield by the PDS will be reported. Finally, the results from the simulation-based study of the supernova burst trigger capability with the PDS of a 12.1-kt active mass dual-phase LAr TPC as the one proposed for DUNE will be summarized. Several configurations of reflective foils installed in the TPC to enhance the light collection will be compared.

Gallego Ros, Ana↗

Charge carrier motion and effect of fixed oxide charge in a microstructured silicon radiation detector

Signal formation in a microstructured semiconductor neutron detector is more complex than in planar diode geometry. Three-dimensional microstructures are laterally smaller than the ionization cloud length, and the electric fields may be weak enough to exhibit plasma time effects. This work is the first detailed treatment of charge carrier motion in these complex semiconductor devices to replicate the time profile and signal magnitude. Simulations were performed using COMSOL Multiphysics to investigate various parameters that affect the propagation of the charge cloud. It was observed that the size of the simulated three-dimensional structure had an impact on the induced current pulse, indicating the importance of simulation geometry optimization to accurately simulate charge cloud expansion. COMSOL Multiphysics was used to replicate accurate charge creation profiles using energy deposition information imported from radiation transport codes. A detailed simulation methodology is presented to benchmark preamplifier event pulses along with complexities in modeling the charge carrier motion along the etched microstructured trenches with Si–SiO 2 boundary conditions, including fixed oxide charge and interface trapping.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Determining the bubble nucleation efficiency of low-energy nuclear recoils in superheated C 3 F 8 dark matter detectors

The bubble nucleation efficiency of low-energy nuclear recoils in superheated liquids plays a crucial role in interpreting results from direct searches for weakly interacting massive particle (WIMP) dark matter. Here, the PICO collaboration presents the results of the efficiencies for bubble nucleation from carbon and fluorine recoils in superheated C 3 F 8 from calibration data taken with five distinct neutron spectra at various thermodynamic thresholds ranging from 2.1 to 3.9 keV. Instead of assuming any particular functional forms for the nuclear recoil efficiency, a generalized piecewise linear model is proposed with systematic errors included as nuisance parameters to minimize model-introduced uncertainties. A Markov chain Monte Carlo routine is applied to sample the nuclear recoil efficiency for fluorine and carbon at 2.45 and 3.29 keV thermodynamic thresholds simultaneously. The nucleation efficiency for fluorine was found to be ≥ 50% for nuclear recoils of 3.3 keV (3.7 keV) at a thermodynamic Seitz threshold of 2.45 keV (3.29 keV), and for carbon the efficiency was found to be ≥ 50% for recoils of 10.6 keV (11.1 keV) at a threshold of 2.45 keV (3.29 keV). Simulated datasets are used to calculate a p value for the fit, confirming that the model used is compatible with the data. The fit paradigm is also assessed for potential systematic biases, which although small, are corrected for. Additional steps are performed to calculate the expected interaction rates of WIMPs in the PICO-60 detector, a requirement for calculating WIMP exclusion limits.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantitative Modeling of High-Energy Electron Scattering in Thick Samples Using Monte Carlo Techniques

Cryo-electron microscopy (cryo-EM) is a powerful tool for imaging biological samples but is typically limited by sample thickness, which is restricted to a few hundred nanometers depending on the electron energy. However, there is a growing need for imaging techniques capable of studying biological samples up to 10 µm in thickness while maintaining nanoscale resolution. This need motivates the use of mega-electron-volt scanning transmission electron microscopy (MeV-STEM), which leverages the high penetration power of MeV electrons to generate high-resolution images of thicker samples. In this study, we employ Monte Carlo simulations to model electron–sample interactions and explore the signal decay of imaging electrons through thick specimens. By incorporating material properties, interaction cross-sections for energy loss, and experimental parameters, we investigate the relationship between the incident and transmitted beam intensities. Key factors such as detector collection angle, convergence semi-angle, and the material properties of samples were analyzed. Our results demonstrate that the relationship between incident and transmitted beam intensities follows the Beer–Lambert law over thicknesses ranging from a few microns to several tens of microns, depending on material composition, electron energy, and collection angles. The linear depth of silicon dioxide reaches 3.9 µm at 3 MeV, about 6 times higher than that at 300 keV. Meanwhile, the linear depth of amorphous ice reaches 17.9 µm at 3 MeV, approximately 11.5 times higher than that at 300 keV. These findings are crucial for advancing the study of thick biological and semiconductor samples using MeV-STEM.

36 MATERIALS SCIENCE↗

Enhancing strontium clock atom interferometry using quantum optimal control

Strontium clock atom interferometry is a promising new technology, with multiple experiments under development around the world to explore its potential for dark matter and gravitational wave detection. In these detectors, large momentum transfer using sequences of many laser pulses is necessary, and thus high-fidelity pulses are important since small errors become magnified. Quantum optimal control (QOC) is a framework for developing control pulse waveforms that achieve high fidelity and are robust against experimental imperfections. Resonant single-photon transitions using the narrow clock transition of strontium involve significantly different quantum dynamics than more established atom interferometry methods based on far-detuned two-photon Raman or Bragg transitions, which leads to new opportunities and challenges when applying QOC. Here, we study in simulation QOC pulses for strontium clock interferometry and demonstrate their advantage over basic square pulses (primitive pulses) and composite pulses in terms of robustness against multiple noise channels. In conclusion, this could improve the scale of large momentum transfer in Sr clock interferometers, paving the way to achieving these scientific goals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fuel-rich oxidation of gasoline surrogate components in an atmospheric flow reactor

Fuel-rich oxidation of three typical gasoline surrogate components, toluene, isooctane, and n-heptane, was investigated in an atmospheric-pressure flow reactor at mean gas temperatures from 1050 to 1350 K, equivalence ratio of 9.0, and residence times of 0.45 and 1.2 s. Not only polycyclic aromatic hydrocarbons (PAHs) up to 3 ring structure but also small intermediate products from C 1 to C 5 were quantified by a gas chromatograph mass spectrometry coupled with photon ionization and gas chromatograph with flame ionization detector, respectively. The kinetic model recently developed by Lawrence Livermore National Laboratory was revised to reflect the results of many recent investigations. Basically, the updated model could satisfactorily reproduce the experimental mole fractions of many species. The experimental and simulated results showed that PAH mole fractions produced were in the order of toluene, isooctane, and n-heptane. The kinetic analysis using the model was carried out to explore PAH formation pathways, especially focusing on naphthalene, acenaphthylene, and phenanthrene. Through rate of production analysis, it was found that the main formation pathways of many PAHs were affected by the fuels. Although resonantly stabilized radicals, such as benzyl and fulvenallenyl radicals, played a crucial role in the formation pathways of many PAHs in every fuel, they were produced through hydrogen elimination of toluene in toluene fuel, while they were formed from small intermediate products in isooctane and n-heptane fuels. Sensitivity analysis revealed the difference and similarity of the reactions with large positive coefficients according to the fuels studied here. In conclusion, the molecular growth reactions of aromatic species were influential in PAH production in every fuel, whereas the ring formation reactions from small species and the reactions involving toluene had large positive sensitivity coefficients in isooctane/n-heptane and toluene fuels, respectively.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Photon detector response function methodology using MCNP and shift hybrid radiation transport code for wide-area contamination assay applications

Here, radiation transport modeling using the Monte Carlo N-Particle (MCNP) radiation transport code and Monte Carlo code, Shift, were employed to model detector responses for a variety of wide-area photon contamination scenarios. In this study, 2" × 2" and 3" × 3" cylindrical NaI(Tl) scintillation detector configurations at source detector-distances of 0.5 cm, 1 cm, 2.54 cm, 10 cm, and 30 cm were modeled. Media of soil, concrete, and steel were evaluated for contamination depths ranging from surface to a depth of an infinite thickness in each medium for photon energies ranging from 20 keV to 3 MeV, which correspond to the energies that current detectors can discern. Monoenergetic photon surface contamination detector responses for each of the media, source–detector distances, and detectors were estimated using MCNP v6.2. Shift was harnessed for improved variance reduction of particle transport in highly attenuating media to obtain average cell fluxes in the two MCNP NaI(Tl) scintillation detector configurations. Average cell flux values in Shift were coupled with detector responses from MCNP to convert average cell flux in a void to energy distribution of pulses in the NaI(Tl) scintillation detector crystal of interest. An optimized detector response function methodology was developed by coupling the MCNP radiation transport method with the Consistent Adjoint Driven Importance Sampling (CADIS) hybrid radiation transport method built into Shift to significantly decrease the runtime of thousands of MCNP pulse height simulations. The methodology may be utilized to quickly and accurately facilitate the assessment of a broad range of wide-area environmental contamination assay and decommissioning cleanup applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Experimental and computational evaluation of alpha particle production from laser-driven proton–boron nuclear reaction in hole-boring scheme

The majority of studies on laser-driven proton–boron nuclear reaction is based on the measurement of α-particles with solid-state nuclear tracks detector (Cr39). However, Cr39's interpretation is difficult due to the presence of several other accelerated particles which can bias the analysis. Furthermore, in some laser irradiation geometries, cross-checking measurements are almost impossible. In this case, numerical simulations can play a very important role in supporting the experimental analysis. In our work, we exploited different laser irradiation schemes (pitcher–catcher and direct irradiation) during the same experimental campaign, and we performed numerical analysis, allowing to obtain conclusive results on laser-driven proton–boron reactions. A direct comparison of the two laser irradiation schemes, using the same laser parameters is presented.

Huault, M. (ORCID:0000000186999166)↗

Evaluation of radiation hardness of High-Voltage silicon vertical JFETs

In the future ATLAS Inner Tracker, each silicon strip module will be equipped with a switch able to separate the high voltage supply from the sensor in case the latter becomes faulty. The switch, placed in between the HV supply and the sensor, needs to sustain a high voltage in its OFF state, to offer a low resistance path for the sensor leakage current in the ON state, and be radiation hard up to 1.2 * 10 15 neq/cm 2 along with other requirements. While GaN JFETs have been selected as suitable rad-hard switch, a silicon vertical HV-JFET was developed by Brookhaven National Laboratory as an alternative option. Pre-irradiation results showed the functionality of the device and proved that the silicon HV-JFET satisfied the pre-irradiation requirements for the switch. To assess its suitability after irradiation, a few p-type HV-JFETs have been neutron irradiated at Jozef Stefan Institute (JSI, Ljubljana, Slovenia). Overall, this paper reports the static characterization of these irradiated devices and the TCAD numerical simulations used to get an insight of the physics governing the post-irradiation behaviour.

47 OTHER INSTRUMENTATION↗

Nuclear recoil detection with color centers in bulk lithium fluoride

We present initial results on the detection of nuclear recoils in lithium fluoride (LiF) through the fluorescence of color centers created by particle interactions in the crystal lattice. Using light-sheet fluorescence microscopy, we image nuclear recoil tracks from both fast and thermal neutron interactions deep within a cubic-centimeter-scale sample. Automated three-dimensional feature extraction based on machine-learning tools enables the identification and classification of individual events. We observe that the fluorescence response of LiF to gamma irradiation is strongly suppressed, by a factor of 30–50 compared to neutron exposure, demonstrating intrinsic insensitivity to electromagnetic backgrounds. The observed and simulated event characteristics are consistent, including their number, size, and topology. These results establish the feasibility of LiF as a scalable detection medium for rare nuclear-recoil events and constitute a first step toward 10–1000 g scale detectors with single-event sensitivity for applications in reactor-neutrino detection, neutron spectroscopy, and dark matter searches.

Aroujo, G R [University of Zurich]↗

ADRIANO2 Calorimeter performance from 2022 Prototypes

A novel high-granularity dual-readout calorimetric technique was developed as part of the T1604 collaboration research program. The building block of ADRIANO2 Calorimeter consists of a pair of optically isolated, small sized tiles made of scintillating plastic and lead glass. The two components of the dual-readout energy compensation technique. Furthermore,Cˇerenkovlightfromtheleadglassareexploitedtoforhighresolutiontimingmeasurements,whilehigh granularity from scintillating plastic can be used to probe the spatial component of the particle shower. This setup works for excellent energy resolution and particle identification for REDTOP as it is crucial for a calorimeter to detect and identify the decay products of eta/eta-prime mesons. Measurements were collected on several ADRIANO2 prototypes between February to December 2022 to evaluate the detector performance. The key metrics extracted from the analyses are presented in this paper are the efficiencies and light-yi eld for various tile configurations. The measured values will be used as input parameters for improved REDTOP Montecarlo simulation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Measurements of Beam Spin Asymmetries in p+p0 and p´p0 Dihadron Production at CLAS12

Semi-Inclusive Deep Inelastic Scattering (SIDIS) is a powerful experimental tool for studying the internal structure and dynamics of the proton, revealing how quarks and gluons are distributed and interact within it. SIDIS describes a process where an elec tron scatters off one of the constituent quarks within the proton, causing it to undergo hadronization, creating multiple hadrons in the final state. Through factorization, the full process can be split into probabilistic components: one which describes the internal structure of the proton using Parton Distribution Functions (PDFs), and another which describes the hadronization process using Fragmentation Functions (FFs). These functions are non-perturbative quantities of Quantum Chromodynamics (QCD), meaning they cannot be calculated directly from first principles and must instead be extracted from experimental measurements. Acommon approach for accessing PDFs and FFs using SIDIS is to measure asymmetries. In this context, asymmetries correspond to subtle differences in the angular distribution of outgoing particles that arise when the spin orientation of the incoming beam or target is reversed. Because many of these effects only appear when spin is involved, they isolate specific, nuanced properties of the proton’s spin-structure that are otherwise hidden in spin averaged measurements. In practice, they show up as specific azimuthal modulations (e.g., sin ¿R, sin(¿h ´ ¿R)), whose amplitudes isolate convolutions of PDFs and FFs at leading and subleading twist. Non-zero asymmetries of these angular distributions can be traced back to unique combinations of PDFs and FFs, offering a way to probe them directly. In this work, we measure SIDIS by analyzing high energy electron-proton scattering events using the CLAS12 detector at Jefferson Lab. This study focuses on subset of SIDIS referred to as dihadron SIDIS, where pairs of hadrons — here p+p0 and p´p0 — are observed. We analyzed these dihadrons using detector data collected during Fall 2018 and Spring 2019, where longitudinally polarized electrons from the CEBAF accelerator were incident on a liquid hydrogen target. A photon classifier using a Gradient Boosted Trees (GBTs) architecture was trained using Monte Carlo simulations to reduce the amount of iv false combinatorial background p0’s. When deployed on experimental data, the model in creases our dihadron statistics by up to five-fold compared to previous CLAS12 p0 analyses. This work reports the first measurements of beam spin asymmetries for p+p0 and p´p0 dihadron production in SIDIS. The measured asymmetries offer new insights to the spin-dependent structure and dynamics within the proton, as well as the spin-dependent properties of quark fragmentation. Non-zero twist-3 sin¿R amplitudes are observed, pro viding sensitivity to the subleading twist PDF e(x). The PDF e(x) encodes quark-gluon correlations within the proton — a property that is otherwise inaccessible at leading twist. Additionally, this work measured significant twist-2 modulations carried by sin(¿h ´ ¿R) and sin(2¿h ´2¿R), providing experimental access to the helicity dihadron fragmentation function (DiFF) GK 1 . Because there is no equivalent quark helicity-dependent FF in single pion SIDIS, the DiFF GK 1 offers a unique lens into novel spin-dependent fragmentation. For instance, the twist-2 modulations observed in this study are enhanced by vector mesons created during fragmentation — a behavior predicted by phenomenological models. This study broadens our understanding of dihadron fragmentation, revealing new details about the flavor and charge dependence of hadronization.

Matousek, Gregory [Duke Univ., Durham, NC (United ↗

Electromagnetic Shower Reconstruction and Identification in FASER's Emulsion Detector for LHC Forward Neutrino Measurements

We present methods for electromagnetic shower reconstruction and identification in the FASERnu emulsion detector using 100 GeV and 200 GeV electron test-beam data from the CERN SPS H4 beamline. The reconstruction employs a clustering-based algorithm without energy-dependent tuning to determine shower axes. A multi-level identification chain comprising track pre-selection, a cut-based selection, and a BDT classifier achieves combined background rejection rates of 99.99% (100 GeV) and 99.94% (200 GeV). The method reaches total reconstruction and identification efficiencies of 58.9% (100 GeV) and 70.8% (200 GeV) evaluated from simulated samples. Energy reconstruction using the total number of reconstructed segments as the calorimetric estimator yields relative biases of +0.6% (100 GeV) and -0.8% (200 GeV), with resolutions of 25.4% and 22.6%, respectively. Systematic uncertainties on the energy reconstruction are dominated by variations in emulsion film detection efficiency, contributing (+10.9%/-8.2%) at 100 GeV and (+10.3%/-6.9%) at 200 GeV. The methodology provides a validated framework for electron neutrino identification with the FASERnu detector at the LHC.

Mammen Abraham, Roshan [UC, Irvine]↗

Determining the bubble nucleation efficiency of low-energy nuclear recoils in superheated C$_3$F$_8$ dark matter detectors

The bubble nucleation efficiency of low-energy nuclear recoils in superheated liquids plays a crucial role in interpreting results from direct searches for weakly interacting massive particle (WIMP) dark matter. The PICO Collaboration presents the results of the efficiencies for bubble nucleation from carbon and fluorine recoils in superheated C$_3$F$_8$ from calibration data taken with 5 distinct neutron spectra at various thermodynamic thresholds ranging from 2.1 keV to 3.9 keV. Instead of assuming any particular functional forms for the nuclear recoil efficiency, a generalized piecewise linear model is proposed with systematic errors included as nuisance parameters to minimize model-introduced uncertainties. A Markov-Chain Monte-Carlo (MCMC) routine is applied to sample the nuclear recoil efficiency for fluorine and carbon at 2.45 keV and 3.29 keV thermodynamic thresholds simultaneously. The nucleation efficiency for fluorine was found to be $\geq 50\, \%$ for nuclear recoils of 3.3 keV (3.7 keV) at a thermodynamic Seitz threshold of 2.45 keV (3.29 keV), and for carbon the efficiency was found to be $\geq 50\, \%$ for recoils of 10.6 keV (11.1 keV) at a threshold of 2.45 keV (3.29 keV). Simulated data sets are used to calculate a p-value for the fit, confirming that the model used is compatible with the data. The fit paradigm is also assessed for potential systematic biases, which although small, are corrected for. Additional steps are performed to calculate the expected interaction rates of WIMPs in the PICO-60 detector, a requirement for calculating WIMP exclusion limits.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗