QUEST-DMC: Background Modelling and Resulting Heat Deposit for a Superfluid Helium-3 Bolometer (2024)

[2,3]\fnmE. \surLeason

1]\orgdivDepartment of Physics, \orgnameLancaster University, \orgaddress \cityLancaster \postcodeLA1 4YB, \countryUK2] \orgdivDepartment of Physics, \orgnameRoyal Holloway University of London, \orgaddress\cityEgham, \postcodeTW20 0EX, \countryUK3] \orgdivDepartment of Physics, \orgnameUniversity of Oxford, \orgaddress\cityOxford, \postcodeOX1 3PJ, \countryUK4] \orgdivDepartment of Physics, \orgnameUniversity of Sussex, \orgaddress\cityBrighton, \postcodeBN1 9QH, \countryUK

\fnmS. \surAutti  \fnmA. \surCasey  \fnmN. \surEng  \fnmN. \surDarvishi  \fnmP. \surFranchini  \fnmR.P. \surHaley  \fnmP.J. \surHeikkinen  \fnmA. \surKemp  elizabeth.leason@rhul.ac.uk  \fnmL.V. \surLevitin  \fnmJ. \surMonroe  \fnmJ. \surMarch-Russel  \fnmM.T. \surNoble  \fnmJ.R. \surPrance  \fnmX. \surRojas  \fnmT. \surSalmon  \fnmJ. \surSaunders  \fnmR. \surSmith  \fnmM.D. \surThompson  \fnmV. \surTsepelin  \fnmS.M. \surWest  \fnmL. \surWhitehead  \fnmK. \surZhang  \fnmD.E. \surZmeev[[[[

Abstract

We report the results of radioactivity assays and heat leak calculations for a range of common cryogenic materials, considered for use in the QUEST-DMC superfluid 3He dark matter detector. The bolometer, instrumented with nanomechanical resonators, will be sensitive to energy deposits from dark matter interactions. Events from radioactive decays and cosmic rays constitute a significant background and must be precisely modelled, using a combination of material screening and Monte Carlo simulations. However, the results presented here are of wider interest for experiments and quantum devices sensitive to minute heat leaks and spurious events, thus we present heat leak per unit mass or surface area for every material studied. This can inform material choices for other experiments, especially if underground operation is considered – where the radiogenic backgrounds will dominate even at shallow depths.

1 Introduction

The nature of dark matter remains an open question in fundamental physics, with extensive direct, indirect and collider searches all returning null results. These searches have typically focused on GeV/c2 - TeV/c2 mass particle dark matter. An increasing number of experiments are also investigating ultra-light boson dark matter with masses much below eV/c2, using techniques to search for wavelike phenomena. However, low mass particle dark matter in the intermediate mass range is not well constrained.

Superfluid helium is an attractive target for low mass dark matter searches due to good kinematic matching, intrinsic radiopurity and small superfluid energy gap.The HeRALD [1] and DELight [2] collaborations are investigating the use of 4He target for a dark matter search. The QUEST-DMC collaboration is exploring the complementary use of superfluid 3He to search for spin-dependent dark matter interactions in the sub-GeV mass range [3]. With a projected energy threshold of 0.51 eV for nuclear recoil interactions, we expect to be able to probe dark matter masses down to 25similar-toabsent25\sim 25∼ 25 MeV/c2MeVsuperscriptc2\rm MeV/c^{2}roman_MeV / roman_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT and spin-dependent cross sections down to 1036similar-toabsentsuperscript1036\sim 10^{36}∼ 10 start_POSTSUPERSCRIPT 36 end_POSTSUPERSCRIPT cm2 with a 4.9 g day exposure, from a 6 month run.

1.1 QUEST-DMC Experiment

The idea of using 3He as a bolometer for particle detection dates back to 1988 [4] and was explored by the MACHe3 [5] and ULTIMA projects [6]. In the QUEST-DMC experiment the superfluid 3He target will be enclosed in a similar-to\sim 1 cm3 transparent box instrumented with a nanomechanical resonators (NEMS) [7]. This is surrounded by a secondary superfluid volume with connection via a similar-to\sim 1 mm2 hole in the bolometer wall. Energy deposition following a dark matter scattering interaction with 3He leads to the production of quasiparticles (broken Cooper pairs) and scintillation photons (following excitation and ionization processes). The quasiparticles are detected as a damping force on the NEMS driven on resonance. Superfluid 3He in the bolometer is cooled to around 100µKtimes100µK100\text{\,}\mathrm{\SIUnitSymbolMicro}\mathrm{K}start_ARG 100 end_ARG start_ARG times end_ARG start_ARG roman_µ roman_K end_ARG to ensure a small thermal population of quasiparticles and maximise sensitivity of the detector to generated quasiparticles. Scintillation photons can be detected using photon sensors surrounding the bolometer. For a complete description of the detector and operation see Ref. [3].

2 Background Modelling

Energy deposition from particles such as cosmic rays or radioactive decay products interacting with target atoms can mimic a dark matter interaction — a significant background in a rare event search. Modelling those events using Monte Carlo simulation and material screening is important for experiment design, assessment of the projected dark matter sensitivity and eventual limit setting. In the context of QUEST-DMC, detailed modelling of energy deposits in the system will also be useful for future studies of superfluid helium physics.

Background sources can be external to the experiment — cosmic rays, neutrons and γ𝛾\gammaitalic_γ rays coming from the surroundings, or internal — radioisotope decays in the detector materials, surfaces or the target itself. The incoming particles can transfer energy to the target through interactions with either electrons or nucleons, which result in the production of quasiparticles and photons described above. External backgrounds depend on the experiment location and can be mitigated using shielding or external veto tagging detectors. Radiogenic backgrounds from detector materials can be minimised using careful design choices in material selection and detector geometry. Internal radiogenic backgrounds from intrinsic contaminants can be minimised by improving material purity and surface contaminants can be minimised using strict cleaning protocols. Once mitigations are implemented it is important to accurately understand the expected background levels, in order to establish or rule out the presence of any candidate dark matter signal.

Superfluid 3He is intrinsically radiopure — at this operating temperature impurities will have frozen out before entering the bolometer. The only other possible isotope is 4He, but the low solubility and preferential adsorption on the cell walls mean that no 4He atoms are expected in the bulk liquid at sub-millikelvin temperatures. Therefore, our assessment of radiogenic backgrounds will focus on naturally occurring radioisotopes embedded in the detector materials. The most common radioisotopes are: 238U, 235U, 232Th, 40K, 60Co and 137Cs. Uranium and thorium isotopes and their progeny form chains which decay through multiple α𝛼\alphaitalic_α and β𝛽\betaitalic_β (and subsequent γ𝛾\gammaitalic_γ) emissions to eventually form stable lead isotopes. The isotopes 40K, 60Co and 137Cs undergo single β𝛽\betaitalic_β decays, with subsequent γ𝛾\gammaitalic_γ ray emission.

The 238U and 232Th decay chains are typically assumed to be in secular equilibrium, due to the long lifetimes of parent nuclei relative to their daughters, allowing a measurement of activity in one part of the chain to determine the activity of the rest. However, secular equilibrium can be broken in both chains by enrichment or removal of radium. In the 238U chain there is a simple equilibrium break at 226Ra, which has a half life of 1600 years so any change will take thousands of years to be restored. The chain can be divided into “early” — for isotopes above 226Ra and “late” — for 226Ra and below, both of which are in secular equilibrium. Similarly, the 232Th chain is split into “early”, above 224Ra, and “late”, including and below 224Ra.

2.1 Material Screening

A range of spectroscopic assay techniques are used to measure radioisotope activity from materials, with each method sensitive to different radiation types and energy ranges. Commonly used measurements include high purity germanium (HPGe) spectrometry, inductively coupled plasma mass spectrometry (ICP-MS), radon emanation and alpha detection. For the uranium and thorium chains at keV-MeV energies HPGe is the most relevant technique and has the advantage of being non-destructive. HPGe assays use a Ge crystal for gamma spectroscopy, to determine levels of naturally occurring radioisotopes in a sample through detection of γ𝛾\gammaitalic_γ rays associated with their decay. This technique cannot distinguish between decays happening on the material surface or in the bulk.

The Boulby UnderGround Screening (BUGS) facility, located 1.1 km underground in Boulby Mine, was used to perform HPGe measurements of materials for the QUEST-DMC experiment. The BUGS facility, originally dedicated to HPGe, contains seven HPGe detectors in a class 1000 cleanroom, with ICP-MS, alpha detection and radon emanation facilities added later [8, 9]. For these measurements the ultra-low background detectors Chaloner and Lunehead and speciality ultra-low background detector Roseberry were used. The detectors are housed inside multi-layer castles, consisting of 10 cm high-purity copper and 10 cm lead, to shield them from environmental γ𝛾\gammaitalic_γ rays. Detector materials used inside the castle are specially selected based on low radioactivity and the manufacturing process is carefully controlled to minimise contaminants. The castle is purged using N2 gas to remove airborne radon, with residual radon in the N2 removed using charcoal traps.

Detectors in the BUGS facility have a range of different types and configurations to allow for a range of different sample geometries and cover a large range of γ𝛾\gammaitalic_γ ray energies. Roseberry is a Mirion BE6530 Broad Energy Germanium (BEGe) planar detector with a 65 cm2 face and 30 mm thickness, giving high sensivity to low energy γ𝛾\gammaitalic_γ rays. Chaloner is a Mirion BE5030 BEGe planar detector with a 50 cm2 face and 30 mm thickness, again giving high efficiency for low energy gammas, but a small volume (150 cm2) best suited for small samples. Lunehead is a p-type coaxial Ortec GEM-XX240-S detector, with a larger (370 cm2) volume but reduced sensitivity to the lowest energy γ𝛾\gammaitalic_γ rays. The 46.5 keV γ𝛾\gammaitalic_γ ray emission from 210Pb cannot be detected by Lunehead, preventing measurement of the 210Pb activity.

Measured activity [mBq/kg]
SampleMass [g]DetectorUearly238superscriptsubscriptUearly238\rm{}^{238}U_{early}start_FLOATSUPERSCRIPT 238 end_FLOATSUPERSCRIPT roman_U start_POSTSUBSCRIPT roman_early end_POSTSUBSCRIPTUlate238superscriptsubscriptUlate238\rm{}^{238}U_{late}start_FLOATSUPERSCRIPT 238 end_FLOATSUPERSCRIPT roman_U start_POSTSUBSCRIPT roman_late end_POSTSUBSCRIPTPb210superscriptPb210\rm{}^{210}Pbstart_FLOATSUPERSCRIPT 210 end_FLOATSUPERSCRIPT roman_PbThearly232superscriptsubscriptThearly232\rm{}^{232}Th_{early}start_FLOATSUPERSCRIPT 232 end_FLOATSUPERSCRIPT roman_Th start_POSTSUBSCRIPT roman_early end_POSTSUBSCRIPTThlate232superscriptsubscriptThlate232\rm{}^{232}Th_{late}start_FLOATSUPERSCRIPT 232 end_FLOATSUPERSCRIPT roman_Th start_POSTSUBSCRIPT roman_late end_POSTSUBSCRIPTK40superscriptK40\rm{}^{40}Kstart_FLOATSUPERSCRIPT 40 end_FLOATSUPERSCRIPT roman_K
Stainless steel544.2Roseberry16(8)16816(8)16 ( 8 )2.5(0.9)2.50.92.5(0.9)2.5 ( 0.9 )82(27)822782(27)82 ( 27 )3.1(1.2)3.11.23.1(1.2)3.1 ( 1.2 )3.9(0.9)3.90.93.9(0.9)3.9 ( 0.9 )<6.2absent6.2<6.2< 6.2
Al 6061-O642.6Lunehead8330(270)83302708330(270)8330 ( 270 )15.3(3.9)15.33.915.3(3.9)15.3 ( 3.9 )-356(12)35612356(12)356 ( 12 )334.4(8.2)334.48.2334.4(8.2)334.4 ( 8.2 )56(8)56856(8)56 ( 8 )
Painted Al923.0Chaloner25680(230)2568023025680(230)25680 ( 230 )16.2(3.2)16.23.216.2(3.2)16.2 ( 3.2 )60480(540)6048054060480(540)60480 ( 540 )259.2(8.3)259.28.3259.2(8.3)259.2 ( 8.3 )342.2(6.2)342.26.2342.2(6.2)342.2 ( 6.2 )21.8(9.6)21.89.621.8(9.6)21.8 ( 9.6 )
Brass107.0Roseberry<7.6absent7.6<7.6< 7.64(1)414(1)4 ( 1 )14990(350)1499035014990(350)14990 ( 350 )<1absent1<1< 1<1.1absent1.1<1.1< 1.1<7.3absent7.3<7.3< 7.3
Silver sinters37.1Roseberry<90absent90<90< 90<36absent36<36< 36430(320)430320430(320)430 ( 320 )<27absent27<27< 27<28absent28<28< 28<385absent385<385< 385
Vespel38.3Chaloner87±66plus-or-minus876687\pm 6687 ± 6690(14)901490(14)90 ( 14 )418±85plus-or-minus41885418\pm 85418 ± 85111(25)11125111(25)111 ( 25 )64(14)641464(14)64 ( 14 )430(240)430240430(240)430 ( 240 )
Fiberglass6.02Chaloner32580(640)3258064032580(640)32580 ( 640 )15154(62)151546215154(62)15154 ( 62 )68600(1400)68600140068600(1400)68600 ( 1400 )11400(100)1140010011400(100)11400 ( 100 )12005(62)120056212005(62)12005 ( 62 )23520(440)2352044023520(440)23520 ( 440 )
Araldite161.9Roseberry<3.6absent3.6<3.6< 3.6<4.8absent4.8<4.8< 4.814.5(9.7)14.59.714.5(9.7)14.5 ( 9.7 )<3.4absent3.4<3.4< 3.4<2.2absent2.2<2.2< 2.2<25.5absent25.5<25.5< 25.5
Stycast131.5Chaloner<10.5absent10.5<10.5< 10.5<9.5absent9.5<9.5< 9.5<14.9absent14.9<14.9< 14.9<12.9absent12.9<12.9< 12.9<6.2absent6.2<6.2< 6.2<122.2absent122.2<122.2< 122.2
GRP106.9Lunehead5700(1000)570010005700(1000)5700 ( 1000 )7460(120)74601207460(120)7460 ( 120 )-7840(160)78401607840(160)7840 ( 160 )7350(100)73501007350(100)7350 ( 100 )4900(570)49005704900(570)4900 ( 570 )
PEN35.1Roseberry<4.2absent4.2<4.2< 4.26.4(2.7)6.42.76.4(2.7)6.4 ( 2.7 )26(13)261326(13)26 ( 13 )<3.4absent3.4<3.4< 3.4<2.4absent2.4<2.4< 2.4<22.8absent22.8<22.8< 22.8
\hdashline
OFHC Cu8800.0GeMPI I<1.9absent1.9<1.9< 1.9<0.13absent0.13<0.13< 0.13-0.090(4)0.09040.090(4)0.090 ( 4 )0.090(4)0.09040.090(4)0.090 ( 4 )0.5(2)0.520.5(2)0.5 ( 2 )
Kapton Cu PCB-Ge Paquito<184.9absent184.9<184.9< 184.9<5.7absent5.7<5.7< 5.7-<4.84absent4.84<4.84< 4.84<2.91absent2.91<2.91< 2.91<17.6absent17.6<17.6< 17.6

In the process of designing the QUEST-DMC experiment and evaluating the sensitivity to dark matter interactions eleven materials, commonly used in ultra-low temperature physics experiments, have been screened for 1-2 weeks at the BUGS facility. Surfaces of the samples were cleaned with lint-free wipes using isopropyl alcohol, to minimise surface contamination. Samples, detector details and measured radioisotope activity levels in these materials are shown in Table 1. The samples screened consisted of cryostat metal parts, other cryostat materials and candidate materials for the experimental cell, described below:

  • stainless steel grade 304, used for the vacuum can – 15 cm square sheet sample, <0.5absent0.5<0.5< 0.5 cm thickness

  • aluminium 6061-O, with and without paint, used to make the helium dewar surrounding the experiment – machined disks, 14-19 cm diameter

  • brass grade CZ121, used for radiation shield cap – machined disk, 6 cm diameter

  • silver sinters, immersed in the 3He for thermal coupling – four blocks, 5 cm length

  • vespel pillars, used as thermally-insulating mechanical supports between different stages of the dilution refrigerator – three hollow cylindrical pillars, 6 cm length

  • copper coated fiberglass PCB – (20×4)204(20\times 4)( 20 × 4 ) cm sheet, <0.5absent0.5<0.5< 0.5 cm thickness

  • Stycast 1266 epoxy manufactured by Henkel (parts A and B, with 100:28 mixing ratio by mass), used for experimental cell – single cuboid (8×6×2.5)similar-toabsent862.5\sim(8\times 6\times 2.5)∼ ( 8 × 6 × 2.5 ) cm

  • Araldite epoxy – two cylindrical pieces, 4 cm diameter, combined height 10similar-toabsent10\sim 10∼ 10 cm

  • glass reinforced plastic (GRP) black nylon 66 with 30% glass reinforcement, possible experimental cell material – three cylindrical pieces, 4 cm diameter, 1similar-toabsent1\sim 1∼ 1 cm height

  • polyethylene naphthalate (PEN), possible cell material or wavelength shifter – 1 m folded sheet

In addition, we can make use of assay results previously reported by other groups, which are extensively catalogued in the SNOLAB radiopurity database [10].Two examples of candidate materials screened by other dark matter experiments are; oxygen-free high conductivity (OFHC) copper C10100 by the XENON1T collaboration [11] and kapton copper monolayer printed PCBs by the TREX collaboration [12]. The screening results for these materials, chosen specifically for their low activities, are shown in Table 1 for comparison. It is also important to note that the level of contaminants in composite materials such as alloys and epoxies can vary significantly between grades or batches. For example, different brass grades can show significant differences in 210Pb activity, as demonstrated in screening results reported in Ref. [13]. Stycast and Araldite epoxy variation between batches depends on the ratios and mixing process e.g. Araldite screening results reported on [10] have Uearly238superscriptsubscriptUearly238\rm{}^{238}U_{early}start_FLOATSUPERSCRIPT 238 end_FLOATSUPERSCRIPT roman_U start_POSTSUBSCRIPT roman_early end_POSTSUBSCRIPT activities varying from 22.2±2.5plus-or-minus22.22.522.2\pm 2.522.2 ± 2.5 mBq/kg to 119.8±6.2plus-or-minus119.86.2119.8\pm 6.2119.8 ± 6.2 mBq/kg.

2.2 Heat Leak

Emitted Power [pW/kg]
SampleAlphaBetaGamma
Unpainted Al20.1(3)20.1320.1(3)20.1 ( 3 )1.15(3)1.1531.15(3)1.15 ( 3 )0.252(4)0.25240.252(4)0.252 ( 4 )
Painted Al110.0(5)110.05110.0(5)110.0 ( 5 )7.21(4)7.2147.21(4)7.21 ( 4 )0.374(3)0.37430.374(3)0.374 ( 3 )
Stainless steel0.14(2)0.1420.14(2)0.14 ( 2 )0.009(2)0.00920.009(2)0.009 ( 2 )0.0036(4)0.003640.0036(4)0.0036 ( 4 )
Brass12.8(3)12.8312.8(3)12.8 ( 3 )1.43(2)1.4321.43(2)1.43 ( 2 )0.022(4)0.02240.022(4)0.022 ( 4 )
Silver sinters0.9(3)0.930.9(3)0.9 ( 3 )0.06(2)0.0620.06(2)0.06 ( 2 )0.04(1)0.0410.04(1)0.04 ( 1 )
Vespel1.4(1)1.411.4(1)1.4 ( 1 )0.07(1)0.0710.07(1)0.07 ( 1 )0.082(9)0.08290.082(9)0.082 ( 9 )
Fiberglass262(1)2621262(1)262 ( 1 )1.5(1)1.511.5(1)1.5 ( 1 )12.51(3)12.51312.51(3)12.51 ( 3 )
Araldite0.06(1)0.0610.06(1)0.06 ( 1 )0.0027(7)0.002770.0027(7)0.0027 ( 7 )0.0037(6)0.003760.0037(6)0.0037 ( 6 )
Stycast0.13(2)0.1320.13(2)0.13 ( 2 )0.004(1)0.00410.004(1)0.004 ( 1 )0.000(2)0.00020.000(2)0.000 ( 2 )
GRP152(2)1522152(2)152 ( 2 )10.0(2)10.0210.0(2)10.0 ( 2 )8.13(5)8.1358.13(5)8.13 ( 5 )
PEN0.07(1)0.0710.07(1)0.07 ( 1 )0.0035(9)0.003590.0035(9)0.0035 ( 9 )0.006(6)0.00660.006(6)0.006 ( 6 )
\hdashline
OFHC Cu0.005(2)0.00520.005(2)0.005 ( 2 )0.0003(1)0.000310.0003(1)0.0003 ( 1 )0.00015(3)0.0001530.00015(3)0.00015 ( 3 )
Kapton Cu PCB4.48(2)4.4824.48(2)4.48 ( 2 )1.12(2)1.1221.12(2)1.12 ( 2 )0.0049(9)0.004990.0049(9)0.0049 ( 9 )
\botrule

The measured activities from material screening can be converted into a heat leak for each material in pW per kg. This is done by considering all α𝛼\alphaitalic_α, β𝛽\betaitalic_β and γ𝛾\gammaitalic_γ emitting processes for the U, Th decay chains and individual radioisotopes. Sensitivity to 235U activity is very low in measurements with the sample masses and exposure times above. When no 235U activity could be measured the theoretical ratio of natural abundance 235U/238U = 0.007257 is used in the heat leak calculations [14]. The power P𝑃Pitalic_P, is calculated by taking the product of energy Eisubscript𝐸𝑖E_{i}italic_E start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, activity aisubscript𝑎𝑖a_{i}italic_a start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, branching ratio bri𝑏subscript𝑟𝑖br_{i}italic_b italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT for every decay in each chain then the summing over all decays of a given type:

P=iEi×ai×bri.𝑃subscript𝑖subscript𝐸𝑖subscript𝑎𝑖𝑏subscript𝑟𝑖P=\sum_{i}E_{i}\times a_{i}\times br_{i}.italic_P = ∑ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_E start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT × italic_a start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT × italic_b italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT .(1)

For the α𝛼\alphaitalic_α and γ𝛾\gammaitalic_γ emissions Eisubscript𝐸𝑖E_{i}italic_E start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT is the discrete energy of the decay, whilst for β𝛽\betaitalic_β emission the mean energy of the emitted spectrum is used. Only decays with branching ratios greater than 1% are considered, branching ratio and decay energies are taken from the ENSDF database [15]. Resulting powers per unit mass of sample are shown in Table 2 for the different decay types.

SampleAlpha power [pW/m2]
Al0.85(1)0.8510.85(1)0.85 ( 1 )
Painted Al5.76(3)5.7635.76(3)5.76 ( 3 )
Stainless0.015(2)0.01520.015(2)0.015 ( 2 )
Araldite0.0021(5)0.002150.0021(5)0.0021 ( 5 )
Stycast0.004(1)0.00410.004(1)0.004 ( 1 )
GRP6.37(50)6.37506.37(50)6.37 ( 50 )
Brass1.15(3)1.1531.15(3)1.15 ( 3 )
Silver sinters0.10(3)0.1030.10(3)0.10 ( 3 )
Vespel1.05(1)1.0511.05(1)1.05 ( 1 )
Fiberglass19.0(5)19.0519.0(5)19.0 ( 5 )
PEN0.0024(5)0.002450.0024(5)0.0024 ( 5 )
\hdashline
OFHC Cu0.0004(1)0.000410.0004(1)0.0004 ( 1 )
Kapton Cu PCB0.016(4)0.01640.016(4)0.016 ( 4 )
\botrule

It is important to note that the stopping powers of the three radiation types will differ, so their expected ranges will vary. For example an α𝛼\alphaitalic_α particle with kinetic energy of 1 MeV will have a range of 3µmsimilar-toabsenttimes3micrometer\sim$3\text{\,}\mathrm{\SIUnitSymbolMicro m}$∼ start_ARG 3 end_ARG start_ARG times end_ARG start_ARG roman_µ roman_m end_ARG in aluminium, whilst a 1 MeV β𝛽\betaitalic_β will have a similar-to\sim 2 mm range and 1 MeV γ𝛾\gammaitalic_γ a similar-to\sim 6 cm range [16]. Therefore many of the emitted α𝛼\alphaitalic_α particles will not escape source material which is more than µmtimesabsentmicrometer\text{\,}\mathrm{\SIUnitSymbolMicro m}start_ARG end_ARG start_ARG times end_ARG start_ARG roman_µ roman_m end_ARG thickness, so we can also report the α𝛼\alphaitalic_α heat leak per unit surface area. Tabulated mass range values (g/cm2) as a function of α𝛼\alphaitalic_α energy are taken from the ASTAR database [16] for different materials and interpolated to find the range corresponding to the decay product energy. The power corresponding to each decay is again calculated from the product of energy, activity, branching ratio and range. The sum is then taken over all decays to give the total power emitted per unit surface area, reported in Table 3.

These results allow for a comparison of relative radioactive emission by different materials, however the decay products can be attenuated or stopped by surrounding materials and the heat leak in a given volume will depend strongly on the geometry and combination of materials used. Power emitted by certain popular construction materials makes them unusable for ultra-low temperature experiments. It is careful choice of materials that ultimately defines the lowest achievable experimental temperatures.

2.3 Background Simulations

QUEST-DMC: Background Modelling and Resulting Heat Deposit for a Superfluid Helium-3 Bolometer (1)

In order to build up a picture of the full energy deposited in an experimental cell a detailed model of the surrounding materials is required. This can be made using the GEANT4 software [17, 18] which simulates interactions of particles with matter across a wide range of energies, tracking their interactions and energy deposits. As described in Ref. [3] a detailed model of the QUEST-DMC detector and cryostat volumes has been constructed and for each volume 105superscript10510^{5}10 start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT-1010superscript101010^{10}10 start_POSTSUPERSCRIPT 10 end_POSTSUPERSCRIPT primary decays are simulated per isotope, depending on distance from the cell. The resulting energy in the cell is recorded and normalised using screening or previous results from the SNOLAB radiopurity database [10]. Energy spectra for all isotopes are summed to find the total radiogenic background expected.

Since the experiment will be located above ground we also expect a significant background from cosmic ray interactions in the cell. This is simulated using GEANT4, plus the CRY library [19] as a particle generator for incident cosmic rays. The cosmic ray flux at the Earth’s surface is normalised to 0.017/cm2/s0.017superscriptcm2s\rm 0.017/cm^{2}/s0.017 / roman_cm start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / roman_s [20], where the uncertainty arising from the measured flux is much smaller than the statistical uncertainty in simulations. Many dark matter experiments operate underground in order to minimise the cosmic ray induced background, for example a muon flux of (4.09±0.15)×108/cm2/splus-or-minus4.090.15superscript108superscriptcm2s(4.09\pm 0.15)\times 10^{-8}\rm/cm^{2}/s( 4.09 ± 0.15 ) × 10 start_POSTSUPERSCRIPT - 8 end_POSTSUPERSCRIPT / roman_cm start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / roman_s has been measured at a depth of 1.1 km in Boulby mine [21].

Figure 1 shows the energy spectra expected from radiogenic decays and cosmic rays in a 0.315 cm3 experimental cell in the Lancaster cryostat described in Ref. [3], operated at saturated vapour pressure and 0.12 Tcsubscript𝑇𝑐T_{c}italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT, where Tcsubscript𝑇𝑐T_{c}italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT is the superfluid transition temperature. The inset plot shows the spectra in the energy region below 100 keV, which is most interesting for dark matter searches. The dominant contribution to radiogenic backgrounds changes across different energy ranges. In the low energy region, below 100 keV, low energy γ𝛾\gammaitalic_γ rays arising from higher activity materials further from the target dominate. In these spectra individual γ𝛾\gammaitalic_γ peaks cannot be distinguished due to coarse binning and statistical fluctuations in the simulations. Rare α𝛼\alphaitalic_α emissions from materials adjacent to the cell will dominate at high energies, above 1 MeV, as they are the only particles that can deposit this amount of energy in the cell. At intermediate energies β𝛽\betaitalic_β emissions from materials close to the cell become important. For cosmic ray backgrounds most of the energy deposits arise from secondary electrons generated in the detector materials. Cosmic muons also deposit energy which depends on path length through the cell, giving the peak seen at similar-to\sim 70 keV. Table 4 shows the average power in the bolometer cell resulting from these energy deposits, where the cosmic ray induced power is shown for both a surface experiment and one located underground e.g. at Boulby.

Events/cell/dayAverage power/cell [pW]
Radiogenic490(20)49020490(20)490 ( 20 )1.15(9)×1041.159superscript1041.15(9)\times 10^{-4}1.15 ( 9 ) × 10 start_POSTSUPERSCRIPT - 4 end_POSTSUPERSCRIPT
Cosmic ray surface5220(70)5220705220(70)5220 ( 70 )6.3(1)×1046.31superscript1046.3(1)\times 10^{-4}6.3 ( 1 ) × 10 start_POSTSUPERSCRIPT - 4 end_POSTSUPERSCRIPT
Cosmic ray underground1.3(1)×1021.31superscript1021.3(1)\times 10^{-2}1.3 ( 1 ) × 10 start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT1.5(1)×1091.51superscript1091.5(1)\times 10^{-9}1.5 ( 1 ) × 10 start_POSTSUPERSCRIPT - 9 end_POSTSUPERSCRIPT
\botrule

3 Conclusion

The QUEST-DMC programme aims to utilise superfluid 3He instrumented with nanomechanical resonators as a bolometer for dark matter detection. Design and realisation of a search for rare interactions requires detailed knowledge of potential background events. These can be modelled using Monte Carlo simulations, normalised using extensive radioassay measurement campaigns. Here, the background modelling efforts for the QUEST-DMC experiment are reported, including germanium screening results for materials commonly used in ultra-low temperature cryostats and comprehensive GEANT4 simulations of both radiogenic and cosmic ray backgrounds. These simulation results have been used to select materials for the design of the QUEST-DMC experiment and evaluate the dark matter sensitivity, as reported in Ref. [3].

Since heating in the experimental cell is of interest beyond the dark matter community, the screening results have also been converted into heat leaks per unit mass or surface area for the different materials. The expected energy spectrum, event rate and heat leak in a single QUEST-DMC cell is shown as an example for a typical ultra-low temperature cryostat. For an experiment on the Earth’s surface, with no dedicated shielding, cosmic ray backgrounds are expected to dominate, particularly at low energies. However, if such an experiment is operated underground e.g. at a depth of 1.1 km in Boulby mine the cosmic ray background is reduced by 6similar-toabsent6\sim 6∼ 6 orders of magnitude. Modelling of muon energy loss at small depths based on Ref. [22], validated using Ref. [23], shows that cosmic muon flux is reduced by more than an order of magnitude for depths greater than 18 m in standard rock (density 2.65 g cm-3). At greater depths the radiogenic heat leak will dominate over cosmic ray energy deposits, limiting the experimental sensitivity – so choice of radiopure materials and construction techniques becomes critical for any underground operation.

An increasing number of cryogenic experiments rely on isolation from interactions with the environment, specifically energy deposits. Successful operation of nuclear demagnetisation cryostats at temperatures of the order 0.1similar-toabsent0.1\sim 0.1∼ 0.1 mK depends on minimising heat leaks to below the pW level [24]. Once thermal and vibrational isolation has been optimised it may be important to consider heat generated by ionising radiation from radioactivity and cosmic rays, as done in Ref. [25]. In recent years superconducting circuits and qubit technology have improved sufficiently that this is also reaching the point of being limited by energy deposits from cosmic rays or radioactive decays [26], which must be well understood to enable robust error correction [27].

\bmhead

AcknowledgementsThis work was funded by UKRI EPSRC and STFC (Grants ST/T006773/1, EP/P024203/1, EP/W015730/1 and EP/W028417/1), as well the European Union’s Horizon 2020 Research and Innovation Programme under Grant Agreement no 824109 (European Microkelvin Platform). M.D.T acknowledges financial support from the Royal Academy of Engineering (RF\201819\18\2), We thank Boulby Underground Laboratory for radioassaying the materials.

References

  • \bibcommenthead
  • [1]S.A. Hertel, A.Biekert, J.Lin, V.Velan, D.N. McKinsey, Direct detection ofsub-GeV dark matter using a superfluid He4superscriptHe4{}^{4}\mathrm{He}start_FLOATSUPERSCRIPT 4 end_FLOATSUPERSCRIPT roman_He target.Phys. Rev. D 100, 092007 (2019).10.1103/PhysRevD.100.092007.URL https://link.aps.org/doi/10.1103/PhysRevD.100.092007
  • [2]B.von Krosigk, etal., Delight: A direct search experiment for light darkmatter with superfluid helium.SciPost Phys. Proc. p. 016 (2023).10.21468/SciPostPhysProc.12.016.URL https://scipost.org/10.21468/SciPostPhysProc.12.016
  • [3]S.Autti, A.Casey, N.Eng, N.Darvishi, P.Franchini, R.P. Haley, P.J.Heikkinen, A.Jennings, A.Kemp, E.Leason, L.V. Levitin, J.Monroe,J.March-Russel, M.T. Noble, J.R. Prance, X.Rojas, T.Salmon, J.Saunders,R.Smith, M.D. Thompson, V.Tsepelin, S.M. West, L.Whitehead, V.V. Zavjalov,D.E. Zmeev, Quest-dmc superfluid 3he detector for sub-gev dark matter.Eur. Phys. J. C 84(3) (2024).10.1140/epjc/s10052-024-12410-8.URL http://dx.doi.org/10.1140/epjc/s10052-024-12410-8
  • [4]G.R. Pickett, Proceedings of the Second European Workshop on neutrinosand Dark Matter Detectors (1988), p. 377
  • [5]F.Mayet, D.Santos, G.Perrin, Y.M. Bunkov, H.Godfrin, Design optimizationof MACHe3, a project of superfluid He-3 detector for direct dark mattersearch.Nucl. Instrum. Meth. A 455, 554–563 (2000).10.1016/S0168-9002(00)00540-4
  • [6]C.B. Winkelmann, J.Elbs, Y.M. Bunkov, E.Collin, H.Godfrin, M.Krusius,Bolometric calibration of a superfluid 3He detector for Dark Mattersearch: Direct measurement of the scintillated energy fraction for neutron,electron and muon events.Nucl. Instrum. Meth. A 574(2), 264–271 (2007).10.1016/j.nima.2007.01.180
  • [7]S.Autti, A.Casey, M.Connelly, N.Darvishi, P.Franchini, J.Gorman, R.P.Haley, P.J. Heikkinen, A.Kemp, E.Leason, J.March-Russell, J.Monroe,T.Noble, G.R. Pickett, J.R. Prance, X.Rojas, T.Salmon, J.Saunders,J.Slater, R.Smith, M.D. Thompson, S.M. West, L.Whitehead, V.V. Zavjalov,K.Zhang, D.E. Zmeev.Long nanomechanical resonators with circular cross-section (2023).URL https://arxiv.org/abs/2311.02452
  • [8]P.Scovell, E.Meehan, H.Araújo, J.Dobson, C.Ghag, H.Kraus,V.Kudryavtsev, X.R. Liu, P.Majewski, S.Paling, R.Preece, R.Saakyan,A.Tomás, C.Toth, L.Yeoman, Low-background gamma spectroscopy at theboulby underground laboratory.Astropart. Phys. 97, 160–173 (2018).10.1016/j.astropartphys.2017.11.006.URL https://doi.org/10.1016%2Fj.astropartphys.2017.11.006
  • [9]P.Scovell, E.Meehan, S.Paling, M.Thiesse, X.Liu, C.Ghag, M.Ginsz,P.Quirin, D.Ralet, Ultra-low background germanium assay at the boulbyunderground laboratory.Journal of Instrumentation 19(01), P01017 (2024).10.1088/1748-0221/19/01/P01017.URL https://dx.doi.org/10.1088/1748-0221/19/01/P01017
  • [10]J.C. Loach, J.Cooley, G.A. Cox, Z.Li, K.D. Nguyen, A.W.P. Poon, A databasefor storing the results of material radiopurity measurements.Nucl. Instrum. Meth. A 839, 6–11 (2016).https://doi.org/10.1016/j.nima.2016.09.036.URL https://www.sciencedirect.com/science/article/pii/S0168900216309639
  • [11]E.Aprile, etal., Material radioassay and selection for the XENON1t darkmatter experiment.Eur. Phys. J. C 77(12) (2017).10.1140/epjc/s10052-017-5329-0.URL https://doi.org/10.1140%2Fepjc%2Fs10052-017-5329-0
  • [12]J.Castel, S.Cebrián, I.Coarasa, T.Dafni, J.Galán, F.J. Iguaz,I.G. Irastorza, G.Luzón, H.Mirallas, A.O. deSolórzano,E.Ruiz-Chóliz, Background assessment for the TREX dark matterexperiment.Eur. Phys. J. C 79(9) (2019).10.1140/epjc/s10052-019-7282-6.URL https://doi.org/10.1140%2Fepjc%2Fs10052-019-7282-6
  • [13]E.Armengaud, etal., Performance of the edelweiss-iii experiment for directdark matter searches.J. Instrum. 12(08), P08010–P08010 (2017).10.1088/1748-0221/12/08/p08010.URL http://dx.doi.org/10.1088/1748-0221/12/08/P08010
  • [14]C.Liu, B.Hu, J.Shi, J.Li, X.Zhang, H.Chen, Determination of uraniumisotopic ratio (235u/238u) using extractive electrospray ionization tandemmass spectrometry.J. Anal. At. Spectrom. 26, 2045–2051 (2011).10.1039/C1JA10054H.URL http://dx.doi.org/10.1039/C1JA10054H
  • [15]NNDC, Evaluated nuclear structure data file (ensdf).Nucl. Instrum. Meth. A 10.18139/nndc.ensdf/1845010.URL https://www.osti.gov/biblio/1845010
  • [16]M.Berger, J.Coursey, M.Zucker.Estar, pstar, and astar: Computer programs for calculatingstopping-power and range tables for electrons, protons, and helium ions(version 1.21) (1999).URL https://physics.nist.gov/PhysRefData/Star/Text/ASTAR.html
  • [17]S.Agostinelli, etal., Geant4 – a simulation toolkit.Nucl. Instrum. Meth. A 506(3), 250–303 (2003).https://doi.org/10.1016/S0168-9002(03)01368-8.URL https://www.sciencedirect.com/science/article/pii/S0168900203013688
  • [18]J.Allison, etal., Recent developments in Geant4.Nucl. Instrum. Meth. A 835, 186–225 (2016).https://doi.org/10.1016/j.nima.2016.06.125
  • [19]Cosmic-ray shower library (cry) version 1.7.URL https://nuclear.llnl.gov/simulation/main.html
  • [20]P.K.F. Grieder, Cosmic rays at Earth : researcher’s reference manual anddata book / Peter K.F. Grieder., first edition. edn. (Elsevier Science Ltd.,Amsterdam ;, 2001)
  • [21]M.Robinson, V.Kudryavtsev, R.Lüscher, J.McMillan, P.Lightfoot,N.Spooner, N.Smith, I.Liubarsky, Measurements of muon flux at verticaldepth in the boulby underground laboratory.Nucl. Instrum. Meth. A 511(3), 347–353 (2003).10.1016/s0168-9002(03)01973-9.URL https://doi.org/10.1016%2Fs0168-9002%2803%2901973-9
  • [22]R.L. Workman, etal., Review of Particle Physics.PTEP 2022, 083C01 (2022).10.1093/ptep/ptac097
  • [23]A.I. Barbouti, B.C. Rastin, A study of the absolute intensity of muons at sealevel and under various thicknesses of absorber.J. Phys. G-Nucl. 9(12), 1577 (1983).10.1088/0305-4616/9/12/018.URL https://dx.doi.org/10.1088/0305-4616/9/12/018
  • [24]D.I. Bradley, A.M. Guénault, V.Keith, C.J. Kennedy, I.E. Miller,S.G. Mussett, G.R. Pickett, W.P. Pratt, New methods for nuclearcooling into the microkelvin regime.Journal of Low Temperature Physics 57(3-4), 359–390 (1984).10.1007/BF00681199
  • [25]E.Nazaretski, V.Kostroun, S.Dimov, R.Pohl, J.Parpia, Heat inputs to sub-mktemperature cryostats and experiments from γ𝛾\gammaitalic_γ-radiation and cosmic raymuons.J. Low Temp. Phys. 137(5), 609–623 (2004).10.1007/s10909-004-0896-5.URL https://doi.org/10.1007/s10909-004-0896-5
  • [26]L.Cardani, etal., Reducing the impact of radioactivity on quantum circuits ina deep-underground facility.Nat. Commun. 12(1) (2021).10.1038/s41467-021-23032-z.URL https://doi.org/10.1038%2Fs41467-021-23032-z
  • [27]C.Wilen, S.Abdullah, N.Kurinsky, C.Stanford, L.Cardani, G.D’Imperio,C.Tomei, L.Faoro, L.Ioffe, C.Liu, A.Opremcak, B.Christensen, J.DuBois,R.McDermott, Correlated charge noise and relaxation errors insuperconducting qubits.Nature 594, 369–373 (2021).10.1038/s41586-021-03557-5
QUEST-DMC: Background Modelling and Resulting Heat Deposit for a Superfluid Helium-3 Bolometer (2024)
Top Articles
Latest Posts
Article information

Author: Kareem Mueller DO

Last Updated:

Views: 6074

Rating: 4.6 / 5 (46 voted)

Reviews: 93% of readers found this page helpful

Author information

Name: Kareem Mueller DO

Birthday: 1997-01-04

Address: Apt. 156 12935 Runolfsdottir Mission, Greenfort, MN 74384-6749

Phone: +16704982844747

Job: Corporate Administration Planner

Hobby: Mountain biking, Jewelry making, Stone skipping, Lacemaking, Knife making, Scrapbooking, Letterboxing

Introduction: My name is Kareem Mueller DO, I am a vivacious, super, thoughtful, excited, handsome, beautiful, combative person who loves writing and wants to share my knowledge and understanding with you.