Next Article in Journal
The Non-Thermal Radio Emissions of the Solar Transition Region and the Proposal of an Observational Regime
Previous Article in Journal
Neutron Star in Quantized Space-Time
Previous Article in Special Issue
Scintillating Bubble Chambers for Rare Event Searches
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Simple Model of the Energy Threshold for Snowball Chambers

1
Department of Physics, University at Albany SUNY, Albany, NY 12222, USA
2
Brookhaven National Laboratory, Upton, NY 11973, USA
*
Author to whom correspondence should be addressed.
Current Address: Department of Physics and Astronomy, University of Alabama, Tuscaloosa, AL 35487, USA.
Current Address: Department of Physics and Astronomy, University of California Los Angeles (UCLA), Los Angeles, CA 90095, USA.
§
Current Address: Department of Physics and Astronomy, Stony Brook University, Stony Brook, NY 11794, USA.
Universe 2024, 10(2), 81; https://doi.org/10.3390/universe10020081
Submission received: 18 December 2023 / Revised: 26 January 2024 / Accepted: 2 February 2024 / Published: 8 February 2024
(This article belongs to the Special Issue Recent Development and Prospects in Dark Matter Research)

Abstract

:
Cloud and bubble chambers have historically been used for particle detection, capitalizing on supersaturation and superheating, respectively. Here, we present new results from a prototype snowball chamber, in which an incoming particle triggers the crystallization of a purified, supercooled liquid. We demonstrate, for the first time, simulation agreement with our first results from 5 years ago: the higher temperature of the freezing of water and significantly shorter time spent supercooled compared to the control in the presence of a Cf-252 fission neutron source. This is accomplished by combining Geant4 modeling of neutron interactions with the Seitz nucleation model used in superheated bubble chambers, including those seeking dark matter. We explore the possible implications of using this new technology for GeV-scale WIMP searches, especially in terms of spin-dependent proton coupling, and report the first supercooling of WbLS (water-based liquid scintillator).

1. Introduction

The nature of dark matter remains an enduring enigma in cosmology and astroparticle physics. A continued lack of unambiguous evidence from any direct detection experiment of the traditional and physically motivated Weakly Interacting Massive Particle (WIMP) has led to an impetus to consider particle candidates with masses higher and lower than before, driven by many hypothetical models [1]. The main goal of our work is the development of inexpensive, scalable, supercooled water detectors for the low-mass dark matter search. Auxiliary purposes for particle physics are neutron detection and neutrino studies [2].
A water target has the advantage of containing hydrogen, ideal for seeking dark matter candidates O(1) GeV/c2 in mass due to the recoil kinematics. An additional advantage is the possibility of a high degree of purification, even on large scales [3]. While having no energy reconstruction, threshold detectors with metastable fluid targets are advantageous for dark matter experiments, due to their high degree of insensitivity to electron recoil backgrounds (BGs) in the search for nuclear recoil signals, as shown by bubble chamber experiments such as COUPP [4] and PICO [5]. Such detectors rely on particles depositing enough energy E within a critical distance, with the differential energy deposition d E / d x exceeding a critical value. By controlling the temperature T and the pressure P, the recoil E threshold can remain low while maintaining a high d E / d x threshold.
Supercooled water is a physics-motivated medium that has been thoroughly studied for its own sake [6,7,8]. Interest in it for particle detection existed years ago [9,10,11] and never for dark matter, until 2016. By contrast, bubble chambers of superheated water are possible but have too high an E threshold for WIMP detection [12]. Another motivation for supercooling over superheating is the fact freezing is exothermic not endothermic like boiling, so the phase transition is entropically favorable. The keV-scale recoil E thresholds achieved by dark matter experiments like PICO should be higher than those possible with supercooling, good for low mass. Existing theories suggest even a sub-eV threshold is likely at an achievable degree of supercooling [13,14,15], with a high d E / d x threshold for BG es at the same T [16].

2. Materials and Methods

Water, first deionized to remove metallic impurities, then filtered (150 nm pore size), was boiled. The resulting steam passed through a series of μ m-scale filters before the final one, a 20-nm PVDF thin-film membrane manufactured by Novamem (Schlieren, Switzerland) similar to that used by [17]. This last filter remained in place above the water during operation (Figure 1 right). A cylindrical, fused-quartz vessel with hemispherical bottom and top flange for sealing from Technical Glass Products (TGP, in Painesville Township, OH, USA) was ultrasonicated in an Alconox soap (Alconox, in White Plains, NY, USA) solution. The ultrasonic cleaning lasted for 15 min at 50 °C and 25 kHz. The vessel was then rinsed with the deionized and filtered water and allowed to dry before being sealed up in a Class-1000 cleanroom (Figure 1 left). A low-power vacuum pump reached ∼1 psia prior to steam flowing.
The quartz ultimately contained 22 ± 1 g of purified water and partial vacuum on top, 8.5 ± 0.5 psia of water vapor at room T. The final mass was limited by the poor throughput of the 20-nm filter, most likely due to particulate build-up. After filling, the quartz vessel was submerged in a ministat circulator for thermal regulation, manufactured by Huber (Raleigh, NC, USA) but sold by Chemglass Life Sciences (Vineland, NJ, USA). It was instrumented with three thermocouples for recording vessel T during the cooling process (T ramp) and during the freezing (exothermic-increase Ts [18]). These thermometers were located near the top (below the flange), middle (water line), and bottom (hemispherical tip).
The liquid water was continuously cooled in an ethanol bath to −35 °C at a linear cool-down rate of ∼2 °C/min, the best rate for the Huber circulating chiller. While introducing a lag in the water T, this had the advantage of reaching a low T and thus higher degree of supercooling rapidly [19], in an effort to reach low E threshold. The chiller sat on vibration-dampening pads for prevention of shock-induced nucleation. All data were read in using National Instruments hardware and their LabView software. Data were taken continuously day and night, alternating control (without source) runs versus source runs, to minimize systematics. An effort was made to ensure equal amounts of control and source data were taken. Our cycle is labeled in the phase diagram in the center plot in Figure 1.
The time spent by the water in the supercooled state was logged for 137Cs 662 keV γ -ray source runs, 252Cf neutron and γ source runs, and interwoven control runs. The T minima achieved in all cases, prior to exothermic rise, were recorded for all data sets as well. All events were included in the final analysis without non-blind data-quality cuts, to further mitigate bias. The time spent “active’’ by the water, Δ t a c t i v e , was defined as the time between crossing −15.5 °C and freezing. This value of −15.5 was determined from the T of the severe drop-off in trigger rate in every data set, although analyses using 0 °C as the border instead yielded consistent results in terms of a statistically significant difference between control and Cf. Less time spent supercooled compared to control implies sensitivity to incoming radiation, with freezing at higher T expected for a fixed cool-down rate.
Figure 1. (Left) The water in the quartz, right after filling (all droplets merged with the volume at bottom after the first cycle). See Figure 1 of [20] for a greater level of detail. (Middle) A phase diagram of water with our cycle in yellow with arrows. (Right) Simple diagram illustrating the salient features of the “non-linear” 20-nm filter: irregular, non-straight holes.
Figure 1. (Left) The water in the quartz, right after filling (all droplets merged with the volume at bottom after the first cycle). See Figure 1 of [20] for a greater level of detail. (Middle) A phase diagram of water with our cycle in yellow with arrows. (Right) Simple diagram illustrating the salient features of the “non-linear” 20-nm filter: irregular, non-straight holes.
Universe 10 00081 g001

3. Results

During a data-taking run of multiple cool-downs for which the 252Cf source was placed in a reproducible location near the liquid, it did not remain in a metastable, i.e., supercooled state as long, freezing also at correspondingly higher T. This main result from the seminal snowball chamber was already covered in [20], which the reader is encouraged to review. In this work, Geant4 (G4) [21,22] Monte Carlo simulations of the nuclear recoil (NR) and electronic recoil (ER) event rates induced by different particle sources placed near the water volume were performed for the purpose of replicating and explaining the data from [20], which has more detail on the Cs results, which did not differ significantly from the control.

3.1. G4 Monte Carlo

In this section we summarize the results of G4 Monte Carlo simulations, whereas in the next section we will apply (critical) thresholds in energy, d E / d x , and track length on top of these, along with a sigmoid-shaped efficiency for the NR-induced snowball nucleation. Table 1 contains all sim(ulation) results by source; its last two columns stem from Figure 2.
As seen on the right, the G4 ER rate is 40–50× higher (before/after thresholds) for 137Cs runs than for 252Cf, accounting for shielding and geometry (Figure 2 inset and Figure 2 in [20]). It is unlikely γ s can explain an enhanced probability for nucleation in 252Cf’s presence.
Results from [25] (60Co) and ref. therein can be explained with T. In these experiments γ s and β s were observed, but not n, so ER not NR were most likely. All thresholds are likely rather sensitive to T; the works with data around −20 °C did not quote errors. A small drop in the d E / d x threshold (probably only 1–2° colder) could steeply increase the Cs ER rate, from negligible to >100 Hz. Future work will test higher-E sources, to look for γ -induced n NR.
This simulated ER vs. NR discrimination power agrees well with our previous measurements presented in [20]. This power, coupled to low E threshold, would make supercooled water detectors suitable as dark matter experiments. Even if possessing only a keV-scale (not sub-keV) threshold, they would be an improvement over current technologies. A light element, even if only oxygen, does not require as low a threshold as heavier elements [26]. While we have the disadvantage of no visible (yet) proton recoils, upon which was based the idea of sensitivity to sub-GeV WIMPs, one must recall that most ER backgrounds are minimally ionizing (≈2 MeV/cm) so we should be able to decrease the d E / d x threshold safely, lowering T. If the snowball chamber acts like a bubble chamber except backwards, colder should imply a lower E threshold, as calculated by Barahona [13,14] and many others.
The results from control and 137Cs runs were nearly identical [20], despite a higher ER rate, suggesting that the snowball chamber, our novel detector, is capable of possessing ER “blindness” similar to PICO’s [5]. The older observations were of needle-like tracks [11,25], not spheres as we observed, implying high-energy es continuously losing E in the water, and/or multiple scattering from γ -rays such as a series of Compton interactions and δ -rays. In contrast, all of our results are significantly distinct, individual snowballs which are analogous to the individual bubbles of dark matter bubble chambers and most likely correspond to neutron elastic scatter vertices (see Figure 3, and Figure 5 inside of [20]).

3.2. Custom Simulation of Thresholds

It is unclear a priori what the best models are for thresholds and the critical radius of nucleation for supercooled fluids. After exploring many similar options, we were able to fit our 252Cf data using the Seitz “hot spike” model that has been used for superheated liquid bubble chambers for decades [29]. The model agrees well with the data, as seen in Figure 4. Figures 2–3 specifically of [16] provided a value for the critical radius of 20 nm, following S w = 0.97 (a quantity related to the pressure and sample geometry). This radius, held constant during fitting, is also a reasonable assumption because of our filter pore size. A PICASSO-modification-inspired Seitz model [30] (for a superheated droplet detector or SDD) that best fits our Cf data is codified by these equations, where subscript c is critical:
E > E c = 0.2 keV n r
After conservatively applying Equations (2) and (3), E c = 1.2 keV n r effectively
d E d x > E c r c = 200 eV 20 nm = 100 MeV / cm ;
l > ( 2 r c ) = 40 nm ;
E f f i c i e n c y = 1 / [ 1 + ( T / ( 252.8 ± 1.1 K ) ) 540 ± 150 ] .
Equation (3), ordinarily implied by (2), is necessary [28]: it ensures a particle traverses at least one critical diameter (l is track length, for a 1D approximation). Otherwise, a proto-snowball, like a proto-bubble, may collapse instead of expanding. Equation (2) was derived by combining Equation (3) with Figure 4 in [16]. Sub-keV NR has ill-defined d E / d x ; (2) is assumed to apply to ER.
Applying Equations (2) and (3) after G4, conservatively to NR too, all ER and NR below ∼1 keV are sub-threshold: see Section 3.1, where we presented the resulting d E / d x and E spectra for all of the possible recoiling species: e, H, and O. G4 was cross-checked using NIST for the first two [31,32] and SRIM [33] for oxygen, and found to be consistent. A severe degeneracy in fitting forced us to not float Equations (1)–(3) above, each already justified by the existing literature, so only a sigmoidal efficiency for nucleation, unrelated to [16], was a free (two-parameter) formula. The sigmoid was an expected addition to the Seitz model given the light elements involved, but it can also be expressed as a softer, sigmoidal energy threshold [34]. Our extreme steepness may explain the reduced nucleation at T −16 °C. This could address the fact that (modern) water supercooling researchers are seemingly unaware of the phenomenon of radiation-induced nucleation, known since at least the 1960s–70s if not earlier, as most of the research into supercooled water has taken place at significantly warmer Ts, far easier to achieve in the lab: less purity is required in terms of critical-radius-sized impurities causing heterogeneous nucleations, in large volumes.
Protons and es fall below the critical (threshold) d E / d x , but not O. This provides a plausible explanation for the lack of decreased supercooled time for Cs, which, despite the thick shielding, like steel flanges, still generated significantly more ER than other sources.
Figure 4 documents a good agreement between Monte Carlo simulation applying a preliminary model and real data on MeV-scale neutrons (which produce keV-scale recoils in H2O) leading to the solidification of water, assuming that it is first supercooled. The data are identical to those in [20], just binned differently. While not the main DAQ trigger, the top thermometer alone defined the reference Ts consistently, as in our seminal paper. It was the only thermocouple for which the JB Weld epoxy did not wear off due to repeated thermal cycling in ethanol. The exact T m i n to use was determined by checking for nucleation in pre-trigger images. (The T increase was the trigger for image recording). Usage of the top thermocouple was another conservative choice: it was farthest from the water and reacted the weakest/latest (in terms of exothermic rise in T). The middle and bottom thermometers exhibited greater statistical significance for the control versus 252Cf difference.
The control data could not be modeled effectively, so only the sim comparison with 252Cf is presented. It is not knowable what fraction of control events was due to residual particulate impurities on the order of the critical radius in size and/or wall/surface events (there was no 3D position information, given only one camera), and what fraction due to normal background radiation in the laboratory on the surface, including cosmic muon-induced neutrons. Thus, the control points could only be crudely subtracted. They created another source of systematic uncertainty. As it currently stands, our best fit to the mean time spent active by the water in 252Cf’s presence, 106.1 ± 2.0 s (data), is 108.3 ± 0.9 s (sim). The T m i n s agree as well: −18.62 ± 0.08 °C (data), −18.60 ± 0.03 °C (sim). (The errors on the simulated quantities are of course smaller, since the statistics can be increased arbitrarily). Not only the centroids, but the widths are reproduced, and even some of the non-Gaussian skewness, most evident below 50 s in Figure 4 left, although systematics from the G4 version and library choices do exist, so we only claim the threshold E is O(keV).
Our modified Seitz model was applied continuously as the T ramped down in simulation to mimic real chiller operation, with all threshold conditions dropping with T, to follow Figures 2–4 from [16]. Equations (1)–(3) show reference values for −20 °C. The fact that Equation (1) is orders of magnitude above even the most pessimistic past calculations of the nucleation activation E suggests a “quenching factor” for neutrons may exist, similar to other detector technologies, where NR take more E than ER to produce the same amount of visible signal [24,35]. On the other hand, metastable-liquid-based detectors should not exhibit this issue, with NR going directly into the signal channel (the phase transition). Further study is warranted, scanning in T and in neutron kinetic E with mono-energetic neutrons, to better calibrate the threshold or activation E.

3.3. Possible Energy Reconstruction

Using a continuous-spectrum source of neutrons like 252Cf incurs a penalty in terms of systematic uncertainty in the threshold energy and other conditions. But, in addition to future calibrations taking place at a mono-energetic neutron-beam facility (such as TUNL), another possible improvement is the addition of scintillation for direct reconstruction of energy. A scintillating snowball chamber should be analogous to the scintillating bubble chamber, invented by Dahl [36]. Water by itself produces Cerenkov radiation for sufficiently high-E particles [3] but, unlike noble elements, for example, such as xenon and argon, it does not scintillate on its own. Water-based liquid scintillator (WbLS), developed by Brookhaven National Laboratory (BNL), addresses this, with scintillation approximately proportional to the amount of E deposited, as expected [37]. An example from the preliminary qualitative results with a snowball chamber based on WbLS instead of water is presented in Figure 5: in the summer of 2023, a 1 mL sample of WbLS was supercooled at BNL, for the very first time anywhere to the best of our knowledge, without any radioactive sources yet, just to establish the baseline behavior of the fluid. WbLS was successfully cooled to −18.5 ± 0.2 °C repeatedly prior to the phase transition occurring.
The micelles of oil-based liquid scintillator inside of the WbLS were too small, being O(1 nm), to act as potential nucleators, thus allowing for the deep level of supercooling. No direct visual evidence of the scintillator coming out of the solution was observed after nine cycles, nor degradation in T m i n . In place of vacuum, 1 mL of immiscible, lower-density (mineral) oil on top of the WbLS served as a buffer liquid to maintain water purity and prevent surface nucleation events, in the style of bubble chambers, as rediscovered in the bio-medical community recently, as a means to more deeply supercool water for organ preservation [38,39]. Without oil, the WbLS could only be supercooled (inconsistently) to ≈−10–15 °C, while an identical volume of water of comparable purity in an identical plastic test tube with the same oil buffer froze at −18.9 ± 0.2 °C, so the control data point agreed with WbLS on T m i n . Another innovation was the first usage of a FLIR (Forward-Looking Infrared) camera, the Teledyne FLIR One Pro for iPhone, to capture the increase in T of supercooled water during its phase transition. It was operated at 8.5 FPS. Use of a plastic tube in place of quartz allowed for a greater thermal (IR) transmission.

4. Discussion

This article documents preliminary fits, based on simulations, for data from the 20 g UAlbany snowball chamber that first demonstrated that MeV-scale neutrons (generating keV-scale recoils) can solidify water, after supercooling. Neutrons, X-rays, and other forms of radiation have all been used to study the microscopic properties of supercooled water before in chemistry [40,41,42], but not its freezing properties for particle physics.
The feasibility of a full-scale water-based dark matter experiment has not yet been established, not until we can further lower the thresholds, so that hydrogen recoils become accessible. However, that should be very easily feasible by going colder through greater purity. Even the lowest T m i n s achieved in the first chamber were still far removed from the coldest possible: <−40 °C for a micro-droplet [43] and −33 °C for a more macroscopic volume [41,44]. Moreover, new filter technology now exists down to 5 nm in pore size.
An increase in the ER background from a reduced threshold in d E / d x could be offset by the replacement of water with WbLS, which may not have sufficient light yield to observe low-E NR from WIMPs or neutrons, but does create enough scintillation light (0.1–1 photons/keV, 10–100× below Xe or Ar) to at least act as a veto for a high-E background [45], as in [36,46].
α backgrounds such as those from radon would be mitigated by means of scintillation or through acoustic discrimination, akin to that discovered by PICASSO then confirmed by COUPP [28,47] (its use continues on PICO), since, unlike most liquids, water expands after it freezes, even given supercooling [44], making it possible that snowballs produce a faint ultrasonic signal, differing by interaction type, in a fashion analogous to expanding bubbles in superheated liquids. β and γ backgrounds have already been addressed in this article, while μ s can be rejected by Cerenkov light and scintillation, and by means of deep-underground deployment, with active or passive shielding acting as a cosmic-ray muon veto. Neutrons would be guaranteed to multiply scatter in a sufficiently large volume.
Spontaneous nucleation in the bulk can be divided into the categories of heterogeneous and homogeneous. The former can be addressed by robust purification, while the latter only happens near −40 °C [27], if at all (doubt has been cast because past works did not account for particle-induced background nucleation in surface laboratories). The problem of surface nucleations at the top of the water can be mitigated by vacuum, oil, or both. Wall events can be prevented by using radiopure smooth quartz, or hydrophobic vessels or coatings, with fiducialization added. The currently low livetime from melting is improvable by purity, modularity, and deployment deep underground to reduce the cosmic-induced event rate.
In light of all of these reasonable possibilities for addressing all of the possible background contributions, a strong sensitivity projection is justifiable even if speculative still. While assuming that proton recoil will be measurable, in a background-free environment, we humbly acknowledge our current lack of precise knowledge regarding the energy threshold, only concluding that its order of magnitude appears to be 1 keV n r based on the only data-sim comparisons currently available. See Figure 6, still conservative overall.

5. Conclusions

We have documented here an initial reproduction of snowball-chamber neutron data using the Seitz model for a superheated liquid, with a well-established modification for softening a step function in efficiency or threshold. MC comparison to the earliest data from 252Cf is already indicative of an energy threshold of ∼0.2–1.2 keV n r , at least for oxygen recoils. To best of the authors’ knowledge, this was the first attempt in the world, in any sub-field of research related to the supercooling of water, to infer the activation energy for nucleation directly from laboratory data through the introduction of interactions of known energies, as opposed to data on the rate of background events as a function of temperature.
Assuming that supercooling acts like superheating—except in reverse—much lower thresholds should be achievable with only slightly lower temperatures than those observed in this work, still a far cry from the lowest known to be possible, even without the modern purification techniques for reducing heterogeneous (non-radiation) nucleations. If possible to achieve without an overwhelming electron-recoil background appearing, sensitivity to recoils on hydrogen nuclei is a valid extrapolation based on this assumption.
We have also shown that supercooling of water-based liquid scintillator is not only possible but possible to the same low temperatures as for pure water, another world-first conclusion, which may lead to the creation of the first scintillating snowball chamber.
The combination of properties discussed here are highly suggestive of the utility of the snowball chamber for direct (WIMP) detection. However, to exceed the present state of the art in terms of spin-independent WIMP limit setting, a lower threshold energy must be established first. Fortunately, thresholds as low 1 eV and even lower in energy have been claimed for many years, with multiple different models of existing data on supercooled water [15,16]. If we conservatively assume that ionization is a necessity and use 12 eV, a snowball chamber would be sensitive to O ( 10 43 ) cm2 cross-section interactions for the mass range of 0.5–10 GeV/c2. However, as a sub-keV threshold is unproven, it would be premature to show an SI projection. We focus instead on the SD proton.

Author Contributions

Conceptualization, C.L. and M.S.; methodology, M.S.; software, M.S.; validation, C.L.; formal analysis, M.S.; investigation, G.J.H.; resources, M.V.D., A.E.B., R.R. and M.Y.; data curation, M.S.; writing—original draft preparation, M.S.; writing—review and editing, A.C.K.; visualization, J.M.; supervision, M.S. and C.L.; project administration, M.S.; funding acquisition, C.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the U.S. DOE Department of Energy (grant# DE-SC0015535) and by an NSF National Science Foundation ADVANCE grant# HRD2121620. The work conducted at Brookhaven National Laboratory was supported by the DOE under contract DE-AC02-98CH10886. Prof. Szydagis thanks BNL Instrumentation Division director Dr. Gabriella Carini for hosting him.

Data Availability Statement

The data presented in this study are available upon request.

Acknowledgments

The authors thank Peter Wilson for his assistance in the publication of the first snowball chamber paper. We also acknowledge Ernst Rehm and Claudio Ugalde of Argonne National Laboratory for useful discussions regarding their water bubble chamber for nuclear astrophysics. We thank Kathy Dunn of CNSE (College of Nanoscale Science & Engineering) for performing electron microscopy on the filtering membrane.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Cooley, J.; Lin, T.; Lippincott, W.H.; Slatyer, T.R.; Yu, T.-T.; Akerib, D.S.; Aramaki, T.; Baxter, D.; Bringmann, T.; Bunker, R.; et al. Report of the Topical Group on Particle Dark Matter for Snowmass 2021. arXiv 2022, arXiv:2209.07426. [Google Scholar]
  2. Akimov, D.; Albert, J.B.; An, P.; Awe, C.; Barbeau, P.S.; Becker, B.; Belov, V.; Brown, A.; Bolozdynya, A.; Cabrera-Palmer, B.; et al. Observation of coherent elastic neutrino-nucleus scattering. Science 2017, 357, 1123–1126. [Google Scholar] [CrossRef]
  3. Fukuda, S.; Fukuda, Y.; Hayakawa, T.; Ichihara, E.; Ishitsuka, M.; Itow, Y.; Kajita, T.; Kameda, J.; Kaneyuki, K.; Kasuga, S.; et al. The Super-Kamiokande detector. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrometers Detect. Assoc. Equip. 2003, 501, 418–462. [Google Scholar] [CrossRef]
  4. Behnke, E.; Behnke, J.; Brice, S.J.; Broemmelsiek, D.; Collar, J.I.; Conner, A.; Cooper, P.S.; Crisler, M.; Dahl, C.E.; Fustin, D.; et al. First dark matter search results from a 4-kg CF3I bubble chamber operated in a deep underground site. Phys. Rev. D 2012, 86, 052001. [Google Scholar] [CrossRef]
  5. Amole, C.; Ardid, M.; Arnquist, I.J.; Asner, D.M.; Baxter, D.; Behnke, E.; Bressler, M.; Broerman, B.; Cao, G.; Chen, C.J.; et al. Dark matter search results from the complete exposure of the PICO-60 C3F8 bubble chamber. Phys. Rev. D 2019, 100, 022001. [Google Scholar] [CrossRef]
  6. Dorsey, N.E. The Freezing of Supercooled Water. Trans. Am. Philos. Soc. 1948, 38, 247–328. [Google Scholar] [CrossRef]
  7. Mossop, S.C. The Freezing of Supercooled Water. Proc. Phys. Soc. Sect. B 1955, 68, 193. [Google Scholar] [CrossRef]
  8. Holten, V.; Sengers, J.V.; Anisimov, M.A. Equation of State for Supercooled Water at Pressures up to 400 MPa. J. Phys. Chem. Ref. Data 2014, 43, 043101. [Google Scholar] [CrossRef]
  9. Agarwal, V.; Gilra, N. Crystal density in supercooled radiation detector. Nucl. Instrum. Methods 1966, 40, 238–240. [Google Scholar] [CrossRef]
  10. Pisarev, A.F. On Ion Detection in Supercooled Liquids; Report number: CERN-TRANS-67-7; CERN: Meyrin, Switzerland, November 1967; 13p. [Google Scholar]
  11. Varshneya, N.C. Detecting Radiation with a Supercooled Liquid. Nature 1969, 223, 826–827. [Google Scholar] [CrossRef]
  12. Deitrich, L.W.; Connolly, T.J. A Study of Fission-Fragment-Induced Nucleation of Bubbles in Superheated Water. Nucl. Sci. Eng. 1973, 50, 273–282. [Google Scholar] [CrossRef]
  13. Barahona, D. Thermodynamic derivation of the activation energy for ice nucleation. Atmos. Chem. Phys. 2015, 15, 13819–13831. [Google Scholar] [CrossRef]
  14. Barahona, D. On the thermodynamic and kinetic aspects of immersion ice nucleation. Atmos. Chem. Phys. 2018, 18, 17119–17141. [Google Scholar] [CrossRef]
  15. Marcolli, C. Technical note: Fundamental aspects of ice nucleation via pore condensation and freezing including Laplace pressure and growth into macroscopic ice. Atmos. Chem. Phys. 2020, 20, 3209–3230. [Google Scholar] [CrossRef]
  16. Khvorostyanov, V.; Curry, J. The Theory of Ice Nucleation by Heterogeneous Freezing of Deliquescent Mixed CCN. Part I: Critical Radius, Energy, and Nucleation Rate. J. Atmos. Sci. 2004, 61, 2676–2691. [Google Scholar] [CrossRef]
  17. Limmer, D.T.; Chandler, D. Phase diagram of supercooled water confined to hydrophilic nanopores. J. Chem. Phys. 2012, 137, 044509. [Google Scholar] [CrossRef]
  18. Nevzorov, A. Some properties of metastable states of water. Phys. Wave Phenom. 2006, 1460, 92. [Google Scholar]
  19. Bigg, E. The Supercooling of Water. Proc. Phys. Soc. Sect. B 1953, 66, 688. [Google Scholar] [CrossRef]
  20. Szydagis, M.; Levy, C.; Huang, Y.; Kamaha, A.C.; Knight, C.C.; Rischbieter, G.R.C.; Wilson, P.W. Demonstration of Neutron Radiation-Induced Nucleation of Supercooled Water. Phys. Chem. Chem. Phys. 2021, 23, 13440–13446. [Google Scholar] [CrossRef] [PubMed]
  21. Agostinelli, S.; Allison, J.; Amako, K.; Apostolakis, J.; Araujo, H.; Arce, P.; Asai, M.; Axen, D.; Banerjee, S.; Barrand, G.; et al. Geant4: A Simulation Toolkit. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrometers Detect. Assoc. Equip. 2003, 506, 250–303. [Google Scholar] [CrossRef]
  22. Allison, J.; Amako, K.; Apostolakis, J.; Araujo, H.; Dubois, P.A.; Asai, M.; Barrand, G.; Capra, R.; Chauvie, S.; Chytracek, R.; et al. Geant4 Developments and Applications. IEEE Trans. Nucl. Sci. 2006, 53, 270–278. [Google Scholar] [CrossRef]
  23. Dahl, C.E. The Physics of Background Discrimination in Liquid Xenon, and First Results from XENON10 in the Hunt for WIMP Dark Matter. Ph.D. Thesis, Princeton University, Princeton, NJ, USA, 2009. [Google Scholar]
  24. Aprile, E.; Angle, J.; Arneodo, F.; Baudis, L.; Bernstein, A.; Bolozdynya, A.; Brusov, P.; Coelho, L.; Dahl, C.; DeViveiros, L.; et al. Design and performance of the XENON10 dark matter experiment. Astropart. Phys. 2011, 34, 679–698. [Google Scholar] [CrossRef]
  25. Varshneya, N. Theory of radiation detection through supercooled liquid. Nucl. Instrum. Methods 1971, 92, 147–150. [Google Scholar] [CrossRef]
  26. McCabe, C. Astrophysical uncertainties of dark matter direct detection experiments. Phys. Rev. D 2010, 82, 023530. [Google Scholar] [CrossRef]
  27. Atkinson, J.D.; Murray, B.J.; O’Sullivan, D. Rate of Homogenous Nucleation of Ice in Supercooled Water. J. Phys. Chem. A 2016, 120, 6513–6520. [Google Scholar] [CrossRef] [PubMed]
  28. Szydagis, M. Dark Matter Limits from a 15 kg Windowless Bubble Chamber. Ph.D. Thesis, The University of Chicago, Chicago, IL, USA, 2011. [Google Scholar] [CrossRef]
  29. Seitz, F. On the Theory of the Bubble Chamber. Phys. Fluids 1958, 1, 2–13. [Google Scholar] [CrossRef]
  30. Behnke, E.; Besnier, M.; Bhattacharjee, P.; Dai, X.; Das, M.; Davour, A.; Debris, F.; Dhungana, N.; Farine, J.; Fines-Neuschild, M.; et al. Final results of the PICASSO dark matter search experiment. Astropart. Phys. 2017, 90, 85–92. [Google Scholar] [CrossRef]
  31. NIST. ESTAR. 2019. Available online: https://physics.nist.gov/PhysRefData/Star/Text/ESTAR.html (accessed on 15 July 2020).
  32. NIST. PSTAR. 2019. Available online: https://physics.nist.gov/PhysRefData/Star/Text/PSTAR.html (accessed on 15 July 2020).
  33. Ziegler, J.F.; Ziegler, M.D.; Biersack, J.P. SRIM—The stopping and range of ions in matter. Nucl. Instrum. Methods Phys. Res. B 2010, 268, 1818–1823. [Google Scholar] [CrossRef]
  34. Behnke, E.; Benjamin, T.; Brice, S.J.; Broemmelsiek, D.; Collar, J.I.; Cooper, P.S.; Crisler, M.; Dahl, C.E.; Fustin, D.; Hall, J.; et al. Direct measurement of the bubble nucleation energy threshold in a CF3I bubble chamber. Phys. Rev. D 2013, 88, 021101. [Google Scholar] [CrossRef]
  35. Aprile, E.; Aalbers, J.; Agostini, F.; Alfonsi, M.; Althueser, L.; Amaro, F.D.; Antochi, V.C.; Angelino, E.; Arneodo, F.; Barge, D.; et al. Light Dark Matter Search with Ionization Signals in XENON1T. Phys. Rev. Lett. 2019, 123, 251801. [Google Scholar] [CrossRef]
  36. Baxter, D.; Chen, C.J.; Crisler, M.; Cwiok, T.; Dahl, C.E.; Grimsted, A.; Gupta, J.; Jin, M.; Puig, R.; Temples, D.; et al. First Demonstration of a Scintillating Xenon Bubble Chamber for Detecting Dark Matter and Coherent Elastic Neutrino-Nucleus Scattering. Phys. Rev. Lett. 2017, 118, 231301. [Google Scholar] [CrossRef]
  37. Alonso, J.R.; Barros, N.; Bergevin, M.; Bernstein, A.; Bignell, L.; Blucher, E.; Calaprice, F.; Conrad, J.M.; Descamps, F.B.; Diwan, M.V.; et al. Advanced Scintillator Detector Concept (ASDC): A Concept Paper on the Physics Potential of Water-Based Liquid Scintillator. arXiv 2014, arXiv:1409.5864. [Google Scholar]
  38. Huang, H.; Yarmush, M.; Usta, B. Long-term deep supercooling of large-volume water via surface sealing with immiscible liquids. Nat. Commun. 2018, 9, 3201. [Google Scholar] [CrossRef] [PubMed]
  39. de Vries, R.J.; Tessier, S.N.; Banik, P.D.; Nagpal, S.; Cronin, S.E.J.; Ozer, S.; Hafiz, E.O.A.; van Gulik, T.M.; Yarmush, M.L.; Markmann, J.F.; et al. Supercooling extends preservation time of human livers. Nat. Biotechnol. 2019, 37, 1131–1136. [Google Scholar] [CrossRef] [PubMed]
  40. Schiller, R.; Kules, I. Radiation chemistry of supercooled water. J. Phys. Chem. 1971, 75, 2997–2999. [Google Scholar] [CrossRef]
  41. Hare, D.E.; Sorensen, C.M. Raman spectroscopic study of bulk water supercooled to −33 °C. J. Chem. Phys. 1990, 93, 25–33. [Google Scholar] [CrossRef]
  42. Belitzky, A.; Mishuk, E.; Ehre, D.; Lahav, M.; Lubomirsky, I. The Source of Electro-Freezing of Super-Cooled Water by Polar Crystals. J. Phys. Chem. Lett. 2016, 7, 43–46. [Google Scholar] [CrossRef]
  43. Goy, C.; Potenza, M.A.C.; Dedera, S.; Tomut, M.; Guillerm, E.; Kalinin, A.; Voss, K.O.; Schottelius, A.; Petridis, N.; Prosvetov, A.; et al. Shrinking of Rapidly Evaporating Water Microdroplets Reveals their Extreme Supercooling. Phys. Rev. Lett. 2018, 120, 015501. [Google Scholar] [CrossRef]
  44. Hare, D.E.; Sorensen, C.M. The density of supercooled water. II. Bulk samples cooled to the homogeneous nucleation limit. J. Chem. Phys. 1987, 87, 4840–4845. [Google Scholar] [CrossRef]
  45. Bat, A.; Tiras, E.; Fischer, V.; Kamislioglu, M. Low energy neutrino detection with a compact water-based liquid scintillator detector. Eur. Phys. J. C 2022, 82, 734. [Google Scholar] [CrossRef]
  46. Alfonso-Pita, E.; Flores, L.J.; Peinado, E.; Vázquez-Jáuregui, E. New physics searches in a low threshold scintillating argon bubble chamber measuring coherent elastic neutrino-nucleus scattering in reactors. Phys. Rev. D 2022, 105, 113005. [Google Scholar] [CrossRef]
  47. Aubin, F.; Auger, M.; Genest, M.H.; Giroux, G.; Gornea, R.; Faust, R.; Leroy, C.; Lessard, L.; Martin, J.P.; Morlat, T.; et al. Discrimination of nuclear recoils from alpha particles with superheated liquids. New J. Phys. 2008, 10, 103017. [Google Scholar] [CrossRef]
  48. Essig, R.; Giovanetti, G.K.; Kurinsky, N.; McKinsey, D.N.; Ramanathan, K.; Stifter, K.; Yu, T.-T.; Aboubrahim, A.; Adams, D.; Alves, D.S.M.; et al. Snowmass2021 Cosmic Frontier: The landscape of low-threshold dark matter direct detection in the next decade. arXiv 2022, arXiv:2203.08297. [Google Scholar]
Figure 2. Stopping power spectra for each possible type of H2O recoil, in the original snowball setup, for tested sources. Corresponding initial species Es for which this is the mean d E / d x are along upper x-axes. Plots are, from left to right: e, p+ (hydrogen), and oxygen recoils, in green or yellow, orange, and cyan, respectively. A ∼100 MeV/cm threshold, a natural assumption as explained later in the text, explains a lack of discernible response from a γ source (diamonds and line near 100 Hz at upper left in the first pane). Lines are not fits but guides; errors are only bin widths (in x) or statistical (y). (Inset) G4 geometry: cross-sectional view. Steel flanges blocked γ s, but ethanol moderated neutrons.
Figure 2. Stopping power spectra for each possible type of H2O recoil, in the original snowball setup, for tested sources. Corresponding initial species Es for which this is the mean d E / d x are along upper x-axes. Plots are, from left to right: e, p+ (hydrogen), and oxygen recoils, in green or yellow, orange, and cyan, respectively. A ∼100 MeV/cm threshold, a natural assumption as explained later in the text, explains a lack of discernible response from a γ source (diamonds and line near 100 Hz at upper left in the first pane). Lines are not fits but guides; errors are only bin widths (in x) or statistical (y). (Inset) G4 geometry: cross-sectional view. Steel flanges blocked γ s, but ethanol moderated neutrons.
Universe 10 00081 g002
Figure 3. Images of the crystallization process: 20 mL of supercooled water (primarily black) turning into light-scattering snow (light blue, lit by blue LED) inside a quartz tube. This corresponds to a double nucleation from neutron data, likely due to multiple scattering. To the best of our knowledge, this was the first time multiple nucleation sites were observed within one sample, in the same image or adjacent frames. Circles mark software detections. Unlike with superheating, nucleations are slow here [27]. There are 150 ms between the frames depicted. The setup was similar to that of the COUPP 15 kg bubble chamber, with lighting orthogonal to camera line of sight [28]. Image artifacts included glare, reflections including off the thermocouples on the quartz wall, the wall itself, or the water line, and/or bubbles/dust within the circulating ethyl alcohol surrounding the TGP quartz.
Figure 3. Images of the crystallization process: 20 mL of supercooled water (primarily black) turning into light-scattering snow (light blue, lit by blue LED) inside a quartz tube. This corresponds to a double nucleation from neutron data, likely due to multiple scattering. To the best of our knowledge, this was the first time multiple nucleation sites were observed within one sample, in the same image or adjacent frames. Circles mark software detections. Unlike with superheating, nucleations are slow here [27]. There are 150 ms between the frames depicted. The setup was similar to that of the COUPP 15 kg bubble chamber, with lighting orthogonal to camera line of sight [28]. Image artifacts included glare, reflections including off the thermocouples on the quartz wall, the wall itself, or the water line, and/or bubbles/dust within the circulating ethyl alcohol surrounding the TGP quartz.
Universe 10 00081 g003
Figure 4. (Left) Histogram of 252Cf active time (10 s bins) with Gaussian fit, in yellow points and line, compared to control in gray for immediately adjacent runs only. Comparison to them lowers the statistical significance of the disagreement between Cf and control. (That was the comparison used conservatively for our original 5 σ claim of a discovery of a new physical-chemistry process involving neutrons freezing supercooled water). However, we also display all control data combined in black, avoiding bias from neutron activation contamination. G4 was used to determine E and d E / d x distributions, with an additional sim adding thresholds and efficiency on top of them: pink. (Right) Same data but looking instead at the T m i n achieved before crystallization, showing consistent results. In total, 74 Cf and 85 local control events were recorded (840 global control).
Figure 4. (Left) Histogram of 252Cf active time (10 s bins) with Gaussian fit, in yellow points and line, compared to control in gray for immediately adjacent runs only. Comparison to them lowers the statistical significance of the disagreement between Cf and control. (That was the comparison used conservatively for our original 5 σ claim of a discovery of a new physical-chemistry process involving neutrons freezing supercooled water). However, we also display all control data combined in black, avoiding bias from neutron activation contamination. G4 was used to determine E and d E / d x distributions, with an additional sim adding thresholds and efficiency on top of them: pink. (Right) Same data but looking instead at the T m i n achieved before crystallization, showing consistent results. In total, 74 Cf and 85 local control events were recorded (840 global control).
Universe 10 00081 g004
Figure 5. Thermal imaging capture of supercooled WbLS freezing. These three example frames are ordered chronologically from left to right and are taken from a 25-second-long clip. The color black is approximately −19 °C, while the purple coloration is 0° on the FLIR unit’s internally calibrated, automatically adjusting scale. The purple-colored artifact present in all images on the left is ice build-up on a thermocouple cable, in the open-air chiller used at BNL. One can see the hemispherical bottom of the test tube come into view very slowly, whereas it was entirely invisible (black) at the start, in thermal equilibrium with the cold air surrounding the vessel. (Note that the top oil was already frozen).
Figure 5. Thermal imaging capture of supercooled WbLS freezing. These three example frames are ordered chronologically from left to right and are taken from a 25-second-long clip. The color black is approximately −19 °C, while the purple coloration is 0° on the FLIR unit’s internally calibrated, automatically adjusting scale. The purple-colored artifact present in all images on the left is ice build-up on a thermocouple cable, in the open-air chiller used at BNL. One can see the hemispherical bottom of the test tube come into view very slowly, whereas it was entirely invisible (black) at the start, in thermal equilibrium with the cold air surrounding the vessel. (Note that the top oil was already frozen).
Universe 10 00081 g005
Figure 6. Projected sensitivities for the spin-dependent interaction of WIMPs with protons as dashed and dotted blue lines for two different snowball-chamber scenarios, compared with the currently world-leading 90% C.L. upper limit on the scattering cross-section in this parameter space, the latest result from the PICO bubble chamber experiment [5], depicted as the solid red line. Even 1.0 kg of water underground for only 1 year would probe (SD p) parameter space untouched by any current experiment, projected by a current experiment, or projected for a known planned experiment [48]. The turn-up due to the E threshold assumed occurs off the scale to the left for this plot, even for 1 keV n r , due to H’s lightness and its nature (single p+). The same SD model as PICO’s was used.
Figure 6. Projected sensitivities for the spin-dependent interaction of WIMPs with protons as dashed and dotted blue lines for two different snowball-chamber scenarios, compared with the currently world-leading 90% C.L. upper limit on the scattering cross-section in this parameter space, the latest result from the PICO bubble chamber experiment [5], depicted as the solid red line. Even 1.0 kg of water underground for only 1 year would probe (SD p) parameter space untouched by any current experiment, projected by a current experiment, or projected for a known planned experiment [48]. The turn-up due to the E threshold assumed occurs off the scale to the left for this plot, even for 1 keV n r , due to H’s lightness and its nature (single p+). The same SD model as PICO’s was used.
Universe 10 00081 g006
Table 1. Simulated event rates in the liquid water by source used. For EZAG’s 252Cf, total n rate is based on having a source identical to that in [23,24]. Source activity at run time is shown along with integrated G4 rates for the active volume, for both sources and recoil types. Parenthetical rates are for E > 1.2 keV, d E / d x > 100 MeV/cm. Note that NR is then strictly oxygen. Upper limits were based on the sim statistics: 1233 and 27.03 s of real time simulated respectively by row. These rates allow us to reproduce the data, in particular those of 252Cf, after more nucleation conditions are applied.
Table 1. Simulated event rates in the liquid water by source used. For EZAG’s 252Cf, total n rate is based on having a source identical to that in [23,24]. Source activity at run time is shown along with integrated G4 rates for the active volume, for both sources and recoil types. Parenthetical rates are for E > 1.2 keV, d E / d x > 100 MeV/cm. Note that NR is then strictly oxygen. Upper limits were based on the sim statistics: 1233 and 27.03 s of real time simulated respectively by row. These rates allow us to reproduce the data, in particular those of 252Cf, after more nucleation conditions are applied.
Calib SourceActivity [ μ Ci]Total RateNR [Hz] n+( γ ,n)ER Rate [Hz]
252Cf fission1.0 (all radiation)≈3000 n/s14.2 (2.7)13.4 (<8 × 10−4)
137Cs gammas10 (100% γ rays)3.7 × 10 5 γ /s<4 × 10−2570 (<4 × 10−2)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Szydagis, M.; Levy, C.; Bolotnikov, A.E.; Diwan, M.V.; Homenides, G.J.; Kamaha, A.C.; Martin, J.; Rosero, R.; Yeh, M. A Simple Model of the Energy Threshold for Snowball Chambers. Universe 2024, 10, 81. https://doi.org/10.3390/universe10020081

AMA Style

Szydagis M, Levy C, Bolotnikov AE, Diwan MV, Homenides GJ, Kamaha AC, Martin J, Rosero R, Yeh M. A Simple Model of the Energy Threshold for Snowball Chambers. Universe. 2024; 10(2):81. https://doi.org/10.3390/universe10020081

Chicago/Turabian Style

Szydagis, Matthew, Cecilia Levy, Aleksey E. Bolotnikov, Milind V. Diwan, George J. Homenides, Alvine C. Kamaha, Joshua Martin, Richard Rosero, and Minfang Yeh. 2024. "A Simple Model of the Energy Threshold for Snowball Chambers" Universe 10, no. 2: 81. https://doi.org/10.3390/universe10020081

APA Style

Szydagis, M., Levy, C., Bolotnikov, A. E., Diwan, M. V., Homenides, G. J., Kamaha, A. C., Martin, J., Rosero, R., & Yeh, M. (2024). A Simple Model of the Energy Threshold for Snowball Chambers. Universe, 10(2), 81. https://doi.org/10.3390/universe10020081

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop