Cosmogenic activation of xenon and copper

Cosmogenic activation of xenon and copper

Abstract

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.

1\thankstext

e1e-mail: lbaudis@physik.uzh.ch \thankstexte2e-mail: alexander.kish@physik.uzh.ch \thankstexte3e-mail: fpiastra@physik.uzh.ch \thankstexte4e-mail: marc.schumann@lhep.unibe.ch

1 Introduction

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
Abundance [%] 0.09(1) 0.09(1) 1.92(3) 26.44(24) 4.08(2) 21.18(3) 26.89(6) 10.44(10) 8.87(16)
T [y] 1.610 [5-12]10 5.810 2.16510
Cu mass number 63 65
Abundance [%] 69.17(3) 30.83(3)
Table 1: Isotopic composition of natural xenon and copper (16). The half-lives of the double-beta decay isotopes Xe, Xe and Xe are from the NuDat 2.6 database (17), the range for the double electron capture isotope Xe is a theoretical prediction from (18).
Location Altitude [m] Atm. depth [g/cm] Vert. flux [mssr] Scaling factor
Sea level 0 1030 2.6 1.0
Lauterbrunnen 795 954 4.7 1.8
LNGS 985 936 5.4 2.1
Jungfraujoch 3470 728 29.0 11.2
Table 2: Altitude, atmospheric depth, and vertical flux of protons and neutrons, as well as scaling factors to relate the flux to sea level values for the various locations in our study.

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.

Figure 1: (a) Relation between atmospheric depth and altitude from the U.S. Standard Atmosphere Model (21). We adopt a sea level value of 1030 g/cm (22); (23). (b) Vertical nucleon flux as a function of the atmospheric depth (24). (c) The cosmic ray spectrum from Activia and Cosmo, based on parametrization from Refs. (30); (31). Only the total flux is used in the calculations as protons and neutrons interact similarly at these energies.

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.

Figure 2: Pre- and post-activation spectra of the 2.04 kg xenon sample. One can identify the Xe lines at 202.9 keV and 375.0 keV, and the I line at 388.6 keV. Other prominent lines are from radioactive contaminations in the stainless steel bottle containing the xenon (primordial U and Th chains, K, cosmogenic Mn, Co).
Figure 3: Pre- and post-activation spectra of the 10.35 kg OFHC copper sample. The entire post-activation spectrum is dominated by the cosmogenic activation products.

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

(1)

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

(2)

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 44) 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).

Figure 4: (a) Posterior PDF for the detected Xe activity, clearly visible in Figure 2. The green region represents the shortest 68.3% C.I., which is taken as 1 uncertainty. (b) Sb, (c) Rh and (d) Be are detected as well, but with larger uncertainties. (e) For Te a signal is present, but an upper limit is reported since the intensity is too low to be determined. (f) An upper limit is reported for the activity of Rh. The right edge of the yellow region represents the 95.5% upper limit (or the lower 95.5% C.I.).

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

(3)

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

(4)

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

(5)

The specific production (activation) rate is equal to the specific saturation activity

(6)

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:

(7)

5 Results

5.1 Copper

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)
Sc 44
V 15.97
Mn 312.12 376 321
Fe 44.50      21   75   
Co 77.24 81    163  153  231
Co 271.74 376 350  376   1022 858 650
Co 70.86 656 632 655  441 1840 1430
Co 1925.28 112 1130  641 537
Table 3: Results for the specific saturation activity of natural copper at sea level, derived from our measurements of the cosmogenic activation. These are compared to our predictions from Activia and Cosmo, using semi-empirical formulae for the cross sections. We also compare to a measurement performed at LNGS (39) (scaled to sea level by a factor 2.1 and corrected for an evaluation error in case of Mn and Fe (40)), to predictions with Activia (11) using () the same semi-empirical formulae and () the MENDL-2P database for the cross sections, to semi-analytical calculations (41) using cosmic ray spectra from () Ziegler and () Gordon et al., and to predictions using TALYS (42). Deviations from our measured values beyond 1 and 1 are indicated by bold or italic font styles, respectively.

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).

5.2 Xenon

Isotope T Xenon: specific saturation activity at sea level [Bq/kg]
[days] This work Literature values
Measurement Calculations Measurement Calculation
Activia Cosmo LUX (47) TALYS (42)
Be 53.3 6.4 6.4
Sr 64.8
Zr 83.4
Nb 62.0
Rh 15.0
Rh 1205.3 16.6 15.3 0.5
Ag 252.0
Sn 115.0
Sb 986.0 0.2 13.5 0.5
Te 154.0
Te 119.7
I 13.0
I 8.04
Xe 36.4 415 555
Xe 8.89
Xe 11.77
Xe 5.25
Cs 6.47 166 164
Table 4: Results for the cosmogenic activation of natural xenon. The specific saturation activities at sea level are compared to our predictions based on Activia and Cosmo, to a measurement by LUX (47) and to a calculation using the TALYS code (42). The half-lifes refer to the numbers used in the software packages. Deviations from our measured values beyond 1 and 1 are indicated by bold or italic font styles, respectively.

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.

Acknowledgements

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.

Footnotes

  1. journal: Eur. Phys. J. C

References

  1. E. Aprile et al. (XENON Collaboration), Astropart. Phys. 35, 573 (2012).
  2. K. Abe et al. (XMASS Collaboration), Nucl. Instr. Meth. A 716, 78 (2013).
  3. D.S. Akerib et al. (LUX Collaboration), Nucl. Instr. Meth. A 704, 111 (2013).
  4. E. Aprile et al. (XENON Collaboration), arXiv:1206.6288.
  5. M. Schumann et al., arXiv:1506.08309.
  6. M. Auger et al. (EXO Collaboration), Phys. Rev. Lett. 109, 032505 (2012).
  7. 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).
  8. J. Martin-Albo et al. (NEXT Collaboration), J. of Phys.: Conf. Ser. 460, 012010 (2013).
  9. G. Plante et al., Phys. Rev. C 84, 045805 (2011).
  10. A. Manalaysay et al., Rev. Sci. Instr. 81, 073303 (2010).
  11. J.J. Back and Y.A. Ramachers, Nucl. Instr. Meth. A 586, 286 (2008).
  12. C.J. Martoff and P.D. Lewin, Comput. Phys. Comm. 72, 96 (1992).
  13. E.Aprile et al. (XENON Collaboration), Phys. Rev. D 83, 082001 (2011).
  14. A. Kish, PhD dissertation, University of Zurich (2011), http://opac.nebis.ch/ediss/20121322.pdf
  15. P. Gondolo, Nucl. Data Sheets 120, 175 (2014).
  16. K.J.R. Rosman and P.D.P. Taylor (IUPAC), Pure and Appl. Chem. 70, 217 (1998).
  17. Nuclear Structure and Decay Database (NuDat 2.6), http://www.nndc.bnl.gov/nudat2/
  18. A. Shukla, P.K. Raina, P.K. Rath, J. Phys. G: Nucl. Phys. 34, 549 (2007),
  19. E. Aprile et al. (XENON Collaboration), Astropart. Phys. 35, 43 (2011).
  20. http://www.ifjungo.ch/jungfraujoch/
  21. U.S. Standard Atmosphere 1976,
    http://www.pdas.com/atmos.html
  22. J.F. Ziegler, IBM J. Res. Develop. 40, 19 (1996).
  23. T. Stanev, High Energy Cosmic Rays, 2nd ed., Chichester UK, Praxis Publishing Ltd. (2010).
  24. J. Beringer et al. (Particle Data Group), Phys. Rev. D86, 010001 (2012).
  25. J.J. Beatty et al., Phys. Rev. D70, 092005 (2004).
  26. B. Heisinger et al., Earth and Planetary Sci. Lett. 200, 345 (2002).
  27. T. Hagner et al., Astropart. Phys. 14, 33 (2000).
  28. 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).
  29. C.J. Martoff, private communication (2015).
  30. T.M. Armstrong, K.C. Chandler, and J. Barish, J. Geophys. Res. 78, 2715 (1973).
  31. N. Gehrels, Nucl. Instr. Meth. A 239, 324 (1985).
  32. L. Baudis et al., JINST 6, P08010 (2011).
  33. S.P. Ahlen et al. (MACRO Collaboration), Phys. Lett. B 249, 149 (1990).
  34. K. Weise and M. Matzke, Nucl. Instr. Meth. A 280, 103 (1989).
  35. K. Weise and W. Woger, Meas. Sci. Techn. 3, 1 (1992).
  36. I. Lira and D. Grientschnig, Metrologia 47, R1 (2010).
  37. S. Agostinelli et al. (GEANT4), Nucl. Instr. Meth. A 506, 250 (2003).
  38. A. Caldwell et al., Comp. Phys. Comm. 180, 2197 (2009).
  39. M. Laubenstein and G. Heusser, Appl. Rad. Isot. 67, 750 (2009).
  40. M. Laubenstein, private communication (2015).
  41. S. Cebrian et al., Astropart. Phys. 33, 316 (2010).
  42. D.-M. Mei and Z.-B. Yin, Astropart. Phys. 31, 417 (2009).
  43. http://www.talys.eu
  44. Yu.N. Shubin et al., MENDL-2P: Proton Reaction Data Library for Nuclear Activation, IAEA-NDS-204 (1998).
  45. J.F. Ziegler, IBM J. Res. Develop. 42, 117 (1998).
  46. M.S. Gordon et al., IEEE Trans. Nucl. Sci. 51, 3427 (2004).
  47. D.S Akerib et al. (LUX Collaboration), Astropart. Phys. 62, 33 (2015).
  48. D.C. Malling et al. (LZ Collaboration), arXiv:1110.0103.
  49. E. Aprile et al. (XENON1T Collaboration), JINST 9, P11006 (2014).
  50. L. Baudis et al., JCAP 01, 044 (2014).
190769
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
Edit
-  
Unpublish
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel
Comments 0
Request comment
""
The feedback must be of minumum 40 characters
Add comment
Cancel
Loading ...

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description