Cosmogenic activation of xenon and copper
Rare event search experiments using liquid xenon as target and detection medium require ultra-low background levels to fully exploit their physics potential. Cosmic ray induced activation of the detector components and, even more importantly, of the xenon itself during production, transportation and storage at the Earth’s surface, might result in the production of radioactive isotopes with long half-lives, with a possible impact on the expected background. We present the first dedicated study on the cosmogenic activation of xenon after 345 days of exposure to cosmic rays at the Jungfraujoch research station at 3470 m above sea level, complemented by a study of copper which has been activated simultaneously. We have directly observed the production of Be, Rh, Sb, I and Xe in xenon, out of which only Sb could potentially lead to background for a multi-ton scale dark matter search. The production rates for five out of eight studied radioactive isotopes in copper are in agreement with the only existing dedicated activation measurement, while we observe lower rates for the remaining ones. The specific saturation activities for both samples are also compared to predictions obtained with commonly used software packages, where we observe some underpredictions, especially for xenon activation.
e1e-mail: firstname.lastname@example.org \thankstexte2e-mail: email@example.com \thankstexte3e-mail: firstname.lastname@example.org \thankstexte4e-mail: email@example.com
Liquid xenon (LXe) is used as detection medium in current and future rare event search experiments, such as direct dark matter (1); (2); (3); (4); (5) and neutrinoless double-beta decay searches (6); (7); (8). It features a high scintillation and ionization yield (9); (10), as well as a high radio-purity. Apart from long-lived double-beta emitters, such as Xe, Xe, Xe and Xe, where only the decay of Xe has been observed so far, xenon has no unstable isotopes. However, the exposure to cosmic rays during production, transportation and storage aboveground can produce instable radio-isotopes via nuclear activation processes.
For a given exposure time and altitude, the activation yield of materials can be predicted using software packages such as Activia (11) and Cosmo (12). Calculations with both codes were performed for the XENON100 experiment (1), however the predicted production rates were too high to be compatible with the measured background rates (13); (14); (15).
In this work we present the first dedicated experimental measurement of the cosmogenic activation of a natural xenon sample, which we exposed to cosmic rays for 345 days at an altitude of 3470 m. In order to provide a benchmark for our activation measurements and predictions, a sample of oxygen-free high thermal conductivity (OFHC) copper has been simultaneously exposed at the same location. Due to its very high purity, and thus low radioactivity levels, OFHC copper is one of the most-frequently used materials to construct low-background detectors, and a well validated material in the software packages.
|Xe mass number||124||126||128||129||130||131||132||134||136|
|Cu mass number||63||65|
|Location||Altitude [m]||Atm. depth [g/cm]||Vert. flux [mssr]||Scaling factor|
The intrinsic radioactivity of the samples, initially stored underground at the Laboratori Nazionali del Gran Sasso (LNGS, Italy) for more than 1.5 y, was measured by means of a high-purity germanium -spectrometer at LNGS, before and after the exposure to cosmic rays. Section 2 describes the samples and their handling, Section 3 the cosmic activation procedure and its modelling, and Section 4 provides details on the -spectrometer and the data analysis. The results of the measurement and the comparison with predictions are presented in Section 5, while conclusions are drawn in Section 6.
2 Samples and preparation
The xenon sample consisted of 2.04 kg research-grade xenon (impurities 10 ppm) from Praxair with natural isotopic composition (see Table 1). It was contained in a 1 liter stainless steel bottle at a pressure of 100 bar (at C). To avoid the measurement results to be dominated by activation products in the bottle ( kg), two identical bottles were used for the -screening and for the activation procedure. The first one always remained underground, where the gas transfer took place as well. Both bottles were evacuated and baked at 100C for several days before being filled with xenon.
The 10.35 kg copper sample consisted of 5 blocks of OFHC copper (purity 99.99%, isotopic composition in Table 1) from Norddeutsche Affinerie (now Aurubis). It came from the batch used to construct inner parts of the XENON100 detector (sample 6 in (19)). Before each measurement, surface contaminations were removed in an ultrasonic bath filled with the acid detergent Elma Clean 60 diluted with deionized water. The sample was then rinsed and wiped with pure ethanol (%), and stored under boil-off N atmosphere for several days in order to let the Rn diffuse out and decay.
3 Activation by cosmic rays
The activation took place in a controlled manner at the High Altitude Research Station Jungfraujoch (20), at an altitude of 3470 m above sea level, from October 31, 2012 to October 15, 2013. The cumulated activation time was 345.0 days. The initial transport from the underground laboratory of LNGS to the Jungfraujoch took less than 5 days and is neglected in the analysis. Following the activation, the samples were brought to Lauterbrunnen (795 m), where they were stored for about 4 days before being transported by car to LNGS for the underground -measurement. Transport to and storage at the LNGS aboveground laboratory (985 m) lasted about 1 day. The cool-down times, from the time when the samples were brought underground until the start of the measurement, were 2.5 days and 14.8 days for the xenon and copper samples, respectively.
For comparison with theoretical predictions and with other measurements, we have to convert our high-altitude activation results into specific saturation activities at sea level. We thus calculate scaling factors that relate the nucleon cosmic ray flux at the various locations, as shown in Table 2. The atmospheric depth at different altitudes was obtained using the U.S. Standard Atmosphere 1976 model (21), shown in Figure 1, with a sea level value of 1030 g/cm (22); (23). The vertical flux of protons and neutrons was then calculated for a given atmospheric depth based on Ref. (24), see Figure 1. Due to their lower flux (25) and lower cross section (26); (27) muons only play a sub-dominant role in activation processes and are neglected here.
The cosmogenic activation was predicted using Activia (11) and an updated and revised version of Cosmo (12); (29). The cosmic ray spectrum encoded in the programs is based on the parametrization from Refs. (30); (31) and is shown in Figure 1. It was sampled from 10 MeV to 10 GeV in 10 MeV intervals: the lower boundary is chosen because the energy thresholds of the nuclear excitation functions are above this energy. The cosmic ray flux above 10 GeV does not affect our results due to the exponential behavior of the energy spectrum. Both packages calculate the cross-sections of the relevant nuclear processes, such as spallation, fission, and evaporation, using semi-empirical formulae developed by Silberberg and Tsao (28), with identical parameters for protons and neutrons. To distinguish between spallation, fission, and nuclear breakdown reactions, including contributions from fast fragmentation processes and transition areas, these formulae are defined in separate regions, depending on the masses of the target and product nuclides (light, intermediate and heavy isotopes). Above a nucleon energy of 3 GeV, the cross-sections are assumed to be energy-independent (11); (28).
Our calculations of the production rates were performed for an exposure to cosmic rays at sea level. To convert into saturation activities, we assumed an activation time of 100 y and no cool-down time.
4 Measurement and data analysis
Before and after activation, the intrinsic radioactivity of the xenon and copper samples was measured with the high-purity germanium -spectrometer Gator (32), operated underground at LNGS. The detector has a background rate of 230 counts/day in the 100-2700 keV interval. Due to the 3500 m water equivalent of rock shielding from cosmic radiation, the hadronic component of the cosmic radiation is completely absent and the atmospheric muon flux is suppressed by 6 orders of magnitude with respect to the aboveground laboratory (33). Hence it is safe to assume that no significant cosmogenic activation took place once the samples were stored underground.
To establish the detection sensitivity to the -lines of the expected activation products as a function of measuring time, the xenon and the copper samples were measured pre-activation for 26.5 days and 34.3 days, respectively. The post-activation spectra were acquired for 11.5 days for xenon and for 4.0 days for copper. The specific activities of the produced radio-isotopes are inferred from the intensities of their most prominent full absorption peaks in the post-activation spectra. The information in the pre-activation measurements is not used in the analysis as no cosmogenic isotopes were present. Both, pre-activation and post-activation spectra are shown in Figures 2 and 3.
A Bayesian method, described in Refs. (34); (35); (36), was employed to infer the activity of an isotope using all its -lines with a sufficiently large branching ratio. The spectrum around each line is divided into three regions, their width is related to the energy-dependent resolution of the detector which is given in (32). The signal region is defined as 3 around the mean position of the full absorption peak. The background is inferred from the count rates in two control regions, +3 above and 3 below the signal region. The total count rate in the signal region is given by
where is the sample mass, is the specific activity of the sample, are the background rates in the left and right control regions, respectively, and are the widths of the three regions. Hence the background rate in the signal region is the interpolation of the background rates in the control regions. The product of branching ratio BR and detection efficiency of the -line is calculated in Monte Carlo simulations, based on a detailed implementation of the detector and sample geometry in GEANT4 (37). The decays of the radio-isotopes of interest are simulated using the G4RadioactiveDecay class, where branching ratios and directional correlations are taken into account.
Considering -lines for the same isotope, the likelihood of the model is
where the are the counts in the three regions for each line , and is the measurement time. is the Poisson distribution.
We use flat priors for the parameters and and generate the posterior probability density function (PDF) of the specific activity by Markov Chain Monte Carlo methods implemented in the Bayesian Analysis Toolkit (BAT) (38). The marginalised posterior PDF is used to decide whether we can claim the detection of a line: if the global mode of the posterior PDF and the left edge of the shortest 68.3% credibility interval (C.I., green region in Figures 4–4) are positive, we calculate an activity, using the mode of the posterior PDF as its estimator and the shortest 68.3% C.I. as 1 uncertainty. If only the global mode is positive, but the left edge of the shortest 68.3% C.I. is zero, we report an upper limit as the signal is too weak to be determined. The upper limit is given as the 95.5% quantile of the posterior PDF (the right edge of the yellow region in Figures 4 and 4).
Because the half-lives of several examined isotopes are comparable to the integral time of the measurement, the measured mean specific activity can significantly differ from the one at the beginning of the measurement
where is the isotope’s mean lifetime, and is the measurement time.
The production rates at Jungfraujoch are calculated from the specific activities at the start of the measurement , taking into account the activation time and the cool-down time . For short-lived nuclides, we must also take into account the rather short time interval , where the samples were exposed to a reduced cosmic ray flux during storage and transportation at lower altitudes, leading to a lower activation rate . The corrected specific activity is
where is the ratio of vertical nucleon fluxes at Lauterbrunnen and the Jungfraujoch. We combine the storage and transport time to days and assume an altitude of 795 m for the whole period. The specific activity of a sample of mass after an activation time is
The specific production (activation) rate is equal to the specific saturation activity
for long activation times. In order to compare with other measurements and calculations, we convert to the specific saturation activity at sea level by scaling it with the factor given in Table 2:
|Isotope||T||Copper: specific saturation activity at sea level [Bq/kg]|
|[days]||This work||Literature values|
|Measurement||Calculations||Measurement||Activia (11)||Calc. (41)||Calculation|
|Activia||Cosmo||LNGS (39)||TALYS (42)|
The radioisotopes of interest were selected based on their half-life (T5 days), the expected production rate, and on their -spectrum: we require at least one line at keV with BR10%, which is separated from a prominent background line by more than 6 .
Table 3 presents the results on the specific saturation activity at sea level for copper, derived from our activation sample. The numbers are compared with our own predictions using the Activia and Cosmo codes, with another measurement performed at LNGS (39), and with additional predictions from the literature: these are based on the Activia package (11), semi-analytical calculations (41), and the TALYS 1.0 code (42); (43). Deviations from our measurement beyond the 1 level are indicated by bold (too high) or italic (too low) font styles.
All cosmogenic radionuclides identified in the activated copper sample are well-known: they are produced in spallation reactions from the stable isotopes Cu and Cu (11); (39); (41). The overall agreement between our measurement and the Activia/Cosmo calculations is remarkable, with most of the predicted activities within the 68% C.I. The highest specific saturation activity, about 0.8 mBq/kg, is measured for Co which has a half-life of 71 d. This value is 20% higher than the prediction. The only other isotope where we measure a higher saturation activity (30%) than predicted is Co. For most of the isotopes the calculations with Cosmo yield systematically lower activities (10%), with the exception of V ( higher than Activia) and Mn ( higher). The general good agreement between measurement and predictions indicates the validity of our implementation of the cosmic ray flux at different altitudes, and can be considered as a benchmark for the comparison of the measurement and predictions for the xenon sample.
For five out of eight isotopes, our results for copper agree with the only other measurement available in the literature, which has been performed at LNGS (39). For Mn, Co and Co, we observe production rates which are , , and times lower, respectively. We obtain identical results with the Activia calculations in (11), when we use the same semi-empirical formulae to calculate the excitation functions (case ). The predictions using the MENDL-2P database (44) (case ) tend to underpredict the production rates. The semi-analytical study (41) predicts much higher production rates than observed for both tested cosmic ray spectra (cases (45), (46)). The TALYS-based work (42) yields reasonable values, which are 20-60% higher than measured, and show better agreement with our measurement than with the one of Ref. (39).
|Isotope||T||Xenon: specific saturation activity at sea level [Bq/kg]|
|[days]||This work||Literature values|
|Activia||Cosmo||LUX (47)||TALYS (42)|
For the xenon sample, we present the measured saturation activities at sea level in Table 4, together with our Activia/Cosmo-based predictions. We also compare the results to measurement performed by the LUX Collaboration (47), and to predictions using the TALYS code (42). To re-scale the LUX-numbers to saturation activities at sea level, we use the procedure described in Section 3 and the information provided in Ref. (47): we assume that all xenon was activated at sea level for 150 days, followed by 49 days (7 days) activation of 50% (50%) of the inventory at the SURF aboveground laboratory. The subsequent cool-down time underground was 90 days (132 days). The atmospheric depth at the SURF altitude of 1600 m is 881 g/cm, corresponding to a vertical nucleon flux of 8.5 mssr. This yields a conversion factor of 3.3, in agreement with the number given by LUX (47).
After activation, we have detected -lines from the following isotopes: Be, Rh, Sb, I and Xe. While our measurement and predictions for I agree within the statistical and systematic errors, this is not the case for the light isotope Be, where the measurement is 50 times higher than the prediction, and for Sb, where we observe a () times higher activity than predicted by Activia and a () times higher activity than the Cosmo prediction. The observed production rate of Xe agrees with the measurement in the LUX detector (47), however, the predictions are about a factor 4 too low. Our sensitivity did not allow us to detect the short-lived xenon isotopes Xe, Xe, and Xe. This also holds for the various other isotopes predicted by Activia and Cosmo, or by the study using TALYS (42). We note that the cosmogenic production rates for xenon isotopes predicted by Cosmo are systematically higher than the ones from Activia.
6 Discussion and Conclusions
We have carried out the first experiment to directly study the cosmogenic activation of xenon, which is employed as target and detection medium for rare event searches, and the second dedicated measurement on the activation of OFHC copper, often used as ultra-pure detector construction material. Both samples were activated by cosmic rays in a controlled manner at the Jungfraujoch research station (3470 m) to maximize the incident nucleon flux. The activation products were measured with a low-background germanium detector in the Gran Sasso underground laboratory. We have compared the measurement results to our own predictions obtained with the Activia (11) and Cosmo (12) software packages, as well as to other measurements and predictions found in the literature.
From the eight detected isotopes in the OFHC copper sample, the production rates of five agree with a measurement performed at LNGS in 2009 (39), which is the only dedicated measurement available in the literature. We measured up to a factor of 3 times lower values for the other three isotopes, see Table 3. While the agreement with our Activia/Cosmo predictions is satisfactory, supporting the validity of our comic ray model for the xenon study, discrepancies with other published calculations are present.
After cosmic activation of the xenon sample, we detected five isotopes, four of which were not directly measured before: Be, Rh, Sb, I. The measured saturation activity of Xe agrees with the measurement by the LUX collaboration (47). In general, the agreement with predictions from Activia and Cosmos is unsatisfactory, the production rate for only one isotope, I, is calculated correctly. The predicted rates for most isotopes are too low by a factor of a few to 10, which is potentially worrisome, as some of these have rather long half-lives (3 years). This conclusion also holds for the short-lived xenon isotopes Xe, Xe and Xe which we could not detect in our study: all predicted numbers are significantly lower than the values measured by LUX (47).
The main sources of background in multi-ton scale LXe-based dark matter detectors are the ones which are uniformly distributed in the target, such as Rn and Kr, or interactions from solar neutrinos (50). Among the observed cosmogenic activation products of xenon, only Rh and Sb have half-lives which are long enough to affect an experiment with a foreseen operation time-scale of 5 years. Both elements have a relatively high electronegativity and could possibly be removed by the xenon purification systems that use hot zirconium getters. However, this has not yet been demonstrated experimentally. The decay of Rh does neither produce low-energetic electrons nor low-energetic -rays not accompanied by a prompt (1 ns) second signal, which shifts the initial low energy signal to high energies by pile-up. We therefore conclude that Rh will not lead to dangerous background. The situation is different for Sb, where 13.6% of the beta-decays end up in a long-lived excited state of the Te daughter ( d) and are therefore unaccompanied by a subsequent -ray. The low-energy tail of these electrons will lead to single-scatter electronic recoil background, while all other decay paths do not contribute, in agreement with the Geant4 prediction. Assuming that Sb is not removed by the purification system, and using the activation times quoted by LUX (47) together with our measured activation rate, a background contribution of events keV kg d would be expected. This rate is too large compared to the published total background level of events keV kg d (47). We therefore conclude that the true activation rate is either close to the lower end of our quoted credibility interval, or Sb is removed by the getter or plates out at surfaces.
In addition, the chemically inert noble gas isotopes could also affect next-generation dark matter searches. We have only observed Xe, which has a relatively short half-life ( d), similarly to other cosmogenically produced Xe-isotopes found in (47). These isotopes will not pose a problem for next-generation LXe detectors, as their contamination will be reduced by an order of magnitude just after a few months of storage below ground. However, we note that for multi-ton scale LXe detectors such as LZ (48), XENONnT (49) and eventually DARWIN (5); (50), the cosmogenic production of radioactive isotopes by the muon flux present at the underground location must be studied as well, as in this case also the decays of the short-lived isotopes might provide a source of backgrounds.
We gratefully acknowledge support from the staff and the International Foundation High Altitude Research Stations Jungfraujoch and Gornergrat (HFSJG) operating the Jungfraujoch laboratory (20), where the activation with cosmic rays was performed, the SNF grants 200020-149256 and 20AS21-136660, the ITN Invisibles (Marie Curie Actions, PITN-GA-2011-289442), Dr. D. Coderre, L. Bütikofer, and Dr. A.D. Ferella for help with transportation and handling the samples at LNGS, and Prof. Dr. C.J. Martoff for discussions on cosmogenic activation.
- journal: Eur. Phys. J. C
- E. Aprile et al. (XENON Collaboration), Astropart. Phys. 35, 573 (2012).
- K. Abe et al. (XMASS Collaboration), Nucl. Instr. Meth. A 716, 78 (2013).
- D.S. Akerib et al. (LUX Collaboration), Nucl. Instr. Meth. A 704, 111 (2013).
- E. Aprile et al. (XENON Collaboration), arXiv:1206.6288.
- M. Schumann et al., arXiv:1506.08309.
- M. Auger et al. (EXO Collaboration), Phys. Rev. Lett. 109, 032505 (2012).
A. Gando et al. (KamLAND-Zen Collaboration), Phys. Rev. C 85, 045504 (2012);
A. Gando et al. (KamLAND-Zen Collaboration), Phys. Rev. Lett. 110, 062502 (2013).
- J. Martin-Albo et al. (NEXT Collaboration), J. of Phys.: Conf. Ser. 460, 012010 (2013).
- G. Plante et al., Phys. Rev. C 84, 045805 (2011).
- A. Manalaysay et al., Rev. Sci. Instr. 81, 073303 (2010).
- J.J. Back and Y.A. Ramachers, Nucl. Instr. Meth. A 586, 286 (2008).
- C.J. Martoff and P.D. Lewin, Comput. Phys. Comm. 72, 96 (1992).
- E.Aprile et al. (XENON Collaboration), Phys. Rev. D 83, 082001 (2011).
- A. Kish, PhD dissertation, University of Zurich (2011), http://opac.nebis.ch/ediss/20121322.pdf
- P. Gondolo, Nucl. Data Sheets 120, 175 (2014).
- K.J.R. Rosman and P.D.P. Taylor (IUPAC), Pure and Appl. Chem. 70, 217 (1998).
- Nuclear Structure and Decay Database (NuDat 2.6), http://www.nndc.bnl.gov/nudat2/
- A. Shukla, P.K. Raina, P.K. Rath, J. Phys. G: Nucl. Phys. 34, 549 (2007),
- E. Aprile et al. (XENON Collaboration), Astropart. Phys. 35, 43 (2011).
U.S. Standard Atmosphere 1976,
- J.F. Ziegler, IBM J. Res. Develop. 40, 19 (1996).
- T. Stanev, High Energy Cosmic Rays, 2nd ed., Chichester UK, Praxis Publishing Ltd. (2010).
- J. Beringer et al. (Particle Data Group), Phys. Rev. D86, 010001 (2012).
- J.J. Beatty et al., Phys. Rev. D70, 092005 (2004).
- B. Heisinger et al., Earth and Planetary Sci. Lett. 200, 345 (2002).
- T. Hagner et al., Astropart. Phys. 14, 33 (2000).
R. Silberberg and C.H. Tsao, Astrophys. J. Suppl. 220 (I) 25, 315 (1973),
R. Silberberg and C.H. Tsao, Astrophys. J. Suppl. 220 (II) 25, 335 (1973),
R. Silberberg, C.H. Tsao, and A.F. Barghouty, Astrophys. J. 501, 911 (1998).
- C.J. Martoff, private communication (2015).
- T.M. Armstrong, K.C. Chandler, and J. Barish, J. Geophys. Res. 78, 2715 (1973).
- N. Gehrels, Nucl. Instr. Meth. A 239, 324 (1985).
- L. Baudis et al., JINST 6, P08010 (2011).
- S.P. Ahlen et al. (MACRO Collaboration), Phys. Lett. B 249, 149 (1990).
- K. Weise and M. Matzke, Nucl. Instr. Meth. A 280, 103 (1989).
- K. Weise and W. Woger, Meas. Sci. Techn. 3, 1 (1992).
- I. Lira and D. Grientschnig, Metrologia 47, R1 (2010).
- S. Agostinelli et al. (GEANT4), Nucl. Instr. Meth. A 506, 250 (2003).
- A. Caldwell et al., Comp. Phys. Comm. 180, 2197 (2009).
- M. Laubenstein and G. Heusser, Appl. Rad. Isot. 67, 750 (2009).
- M. Laubenstein, private communication (2015).
- S. Cebrian et al., Astropart. Phys. 33, 316 (2010).
- D.-M. Mei and Z.-B. Yin, Astropart. Phys. 31, 417 (2009).
- Yu.N. Shubin et al., MENDL-2P: Proton Reaction Data Library for Nuclear Activation, IAEA-NDS-204 (1998).
- J.F. Ziegler, IBM J. Res. Develop. 42, 117 (1998).
- M.S. Gordon et al., IEEE Trans. Nucl. Sci. 51, 3427 (2004).
- D.S Akerib et al. (LUX Collaboration), Astropart. Phys. 62, 33 (2015).
- D.C. Malling et al. (LZ Collaboration), arXiv:1110.0103.
- E. Aprile et al. (XENON1T Collaboration), JINST 9, P11006 (2014).
- L. Baudis et al., JCAP 01, 044 (2014).