arxiv:astro-ph/ v2 10 Jan 2007

Transcription

1 Measurement of the specific activity of 39 Ar in natural argon arxiv:astro-ph/ v2 10 Jan 2007 P. Benetti a, F. Calaprice b, E. Calligarich a, M. Cambiaghi a, F. Carbonara c, F. Cavanna d, A. G. Cocco c, F. Di Pompeo d, N. Ferrari e,1, G. Fiorillo c, C. Galbiati b, L. Grandi a, G. Mangano c, C. Montanari a, L. Pandola e,, A. Rappoldi a, G. L. Raselli a, M. Roncadelli a, M. Rossella a, C. Rubbia a, R. Santorelli c, A. M. Szelc f, C. Vignoli a and Y. Zhao b a INFN and Department of Physics, University of Pavia, Pavia, Italy b Department of Physics, Princeton University, Princeton NJ, USA c INFN and Department of Physics, University of Napoli Federico II, Napoli, Italy d INFN and Department of Physics, University of L Aquila, L Aquila, Italy Abstract e INFN, Gran Sasso National Laboratory, Assergi, Italy f Institute of Nuclear Physics PAN, Kraków, Poland We report on the measurement of the specific activity of 39 Ar in natural argon. The measurement was performed with a 2.3-liter two-phase (liquid and gas) argon drift chamber. The detector was developed by the WARP Collaboration as a prototype detector for WIMP Dark Matter searches with argon as a target. The detector was operated for more than two years at Laboratori Nazionali del Gran Sasso, Italy, at a depth of 3,400 m w.e. The specific activity measured for 39 Ar is 1.01±0.02(stat)±0.08(syst) Bq per kg of nat Ar. Key words: 39 Ar specific activity, low-background experiments, cosmogenic activation PACS: s, z, Mc, d INFN, Laboratori Nazionali del Gran Sasso, S.S. 17 bis km , I Assergi (AQ), Italy. Telephone number: , Fax number: address: pandola@lngs.infn.it (L. Pandola). 1 Deceased Preprint submitted to Nucl. Instrum. Methods A 5 February 2008

2 1 Introduction A 2.3-liter two-phase (liquid and gas) argon drift chamber [1] was developed and built by the WARP collaboration as a prototype detector for WIMP Dark Matter searches with argon as a target [2]. The detector was operated for more than two years by the WARP collaboration at Laboratori Nazionali del Gran Sasso, Italy, at a depth of 3,400 m w.e. One important by-product of the operation of the prototype WARP detector was the precise determination of the 39 Ar specific activity in natural argon. 39 Ar and 85 Kr are two radioactive nuclides whose activity in the atmosphere is of the order of 10mBq/m 3 and 1Bq/m 3, respectively [3,4]. As a result of the liquid argon production process, they are both present in abundant quantities and are the two most significant radioactive contaminations in liquid argon. The two isotopes decay primarily by β emission, and their presence can limit the sensitivity of experiments looking for low energy rare events (WIMP Dark Matter interactions, neutrinoless double beta decay) using liquid argon either as a target or as a shielding material. 85 Kr is not a pure β emitter, owing to the presence of a 0.43% branching ratio for decay with β emission on a metastable state of 85 Rb, which then decays by emitting a γ-ray of energy 514 kev, with a half-life of 1.01µs [5]. The coincidence between the β and γ emitted in a fraction of the 85 Kr decays may ease in some cases the task of determining experimentally the activity of 85 Kr in low-background detectors. The determination of the specific activity of 39 Ar is intrinsically more challenging and not as widely discussed in the literature [3,6]. A theoretical estimate is presented in Ref. [7]. In the last twenty years liquid argon technology has acquired great relevance for astroparticle physics applications. Several experimental techniques, employing liquid argon as sensitive medium, have been proposed especially for rare events detection [2,8,9,10,11,12,13,14]. For WIMP Dark Matter direct detection, the discrimination of nuclear recoils from the β-γ induced background plays a crucial role. The discrimination provided by the experimental technique must be sufficient to reduce the radioactive background below the very low interaction rates foreseeable for WIMP Dark Matter. A precise determination of the intrinsic specific activity of 39 Ar is therefore of significant interest for the design of WIMP Dark Matter detectors employing argon as a target. 2 The 2.3-liter WARP detector The detector consists of a two-phase argon drift chamber with argon as a target. Two-phase argon drift chamber was first introduced within the ICARUS program [9] in the framework of a wide-range study of the properties of noble 2

3 Fig. 1. Schematic view of the WARP 2 liter detector. gases. The drift chamber (see Figure 1) is operated at the argon boiling point (86.7K) at the atmospheric pressure of the Gran Sasso Laboratory (about 950 mbar) [15]. The cooling is provided by a thermal bath of liquid argon, contained in an open stainless steel dewar, in which the chamber is fully immersed. The pressure of the gas phase on the top of the chamber is naturally equalized to the surrounding atmospheric pressure. Ionizing events inside the liquid argon volume produce scintillation VUV light mainly at 128nm. The scintillation light is shifted to wavelengths in the blue by an organic wavelength shifter (TetraPhenilButadiene) covering the walls, and collected by photomultipler tubes (PMTs) located in the gas phase and facing the liquid volume below. The 2-inch PMTs are manufactured by Electron Tubes Ltd (model D757UFLA) and have a special photocathode that ensures functionality down to liquid argon temperatures. The quantum efficiency at the emission wavelength of TPB is about 18%. The material of the PMTs has been selected for high radiopurity: according to the supplier s specifications, the total γ activity above 100 kev is 0.2 Bq/PMT, dominated by the 232 Th and 238 U chains. A series of field-shaping rings surrounding the liquid phase superimpose an electric field of 1kV/cm. The electrons are drifted toward the anode (located atop the chamber) and then extracted from the liquid to the gaseous phase by a local extraction field provided by a couple of grids. The electrons are linearly multiplied in the gas phase by a second, stronger local field. The PMTs detect the primary scintillation light (directly produced by the ionizing event) 3

4 and also the secondary scintillation light (produced by the electron multiplication process in the gas phase). The PMT signals are summed and sent to a multi-channel analyzer recording the pulse height spectrum. The liquid argon contained in the chamber is Argon 6.0 supplied by Rivoira S.p.A. The liquid argon is successively purified from electronegative impurities down to an equivalent contamination of less than 0.1ppb of O 2 by using the chemical filter Hopkalit TM from Air Liquid. The purity from electronegative elements is actively maintained by means of continuous argon recirculation through the chemical filter. The experimental set-up is located in the Laboratori Nazionali del Gran Sasso underground laboratory (3,400 m w.e. of rock coverage). The flux of cosmic ray muons is suppressed by a factor of 10 6 with respect to the surface (residual flux 1.1µ/(m 2 h), average muon energy 320GeV [16]). The detector is shielded by 10cm of lead, to reduce the external γ background. The sensitive volume of the detector has the shape of frustum of cone and it is delimited by a stainless-steel cathode. The sensitive volume for the configuration under analysis is 1.86±0.07 liter. The density of liquid argon in the operating conditions (950mbar and 86.7K) is1.399 g/cm 3 [15]. The sensitive volume is viewed by seven PMTs, whose responses have been equalized in gain. Daily calibrations ensure the long-term stability and the linearity of the response. The sensitive volume and the argon thermal bath are contained in a stainless steel dewar, 50 cm internal diameter and 200 cm internal height. 3 Data analysis For the measurements described in this work, the electric fields were switched off and the chamber was operated as a pure scintillation detector. The gain of the PMTs has been set up to optimize the data acquisition in the typical energy range of the environmental γ-ray background, namely up to 3 MeV. The energy threshold for data acquisition is about 40 kev. The threshold used for analysis is 100 kev, in order to exclude events from electronic noise. The response of the detector to γ radiation was studied using different γ-ray sources ( 57 Co, 60 Co, 137 Cs) placed outside the chamber. The spectra obtained with the 57 Co and 137 Cs sources are shown in Figure 2. Typical values of the resolution observed with the calibration sources are σ(e)/e = 13% at 122 kev ( 57 Co) and σ(e)/e = 6% at 662 kev ( 137 Cs). The correlation between energy and primary scintillation light detected was linear within the range tested with our sources. The energy resolution of the detector can be described empirically by the following parametrization: σ(e) = a a 1 E + (a 2 E) 2, (1) 4

5 (a) (b) Fig. 2. Energy spectra taken with external γ-ray sources, superimposed with the corresponding Monte Carlo simulations. (a) 57 Co source (E = 122 kev, B.R. 85.6%, and 136 kev, B.R. 10.7%), (b) 137 Cs source (E = 662 kev). where a 0 = 9.5 kev, a 1 = 1.2 kev and a 2 = The three terms take into account the effects from non-uniform light collection (a 2 term), statistical fluctuations in the light production (a 1 term) and electronic noise (a 0 term). The β-γ spectrum in the detector (35 hours of live time) is displayed in Figure 3. The total counting rate is about 6 Hz (4.2 Hz above the analysis threshold). A simulation based on the Geant4 toolkit [17] has been developed to reproduce and understand the observed features. The simulation takes into account the following two components of background: External radiation from radioactive contaminants in the materials surrounding the liquid argon volume (stainless steel, thermal bath, PMTs). Only γ- emitters are taken into account, that originate an effective γ-ray flux through the surface of the sensitive volume. This includes: 238 U and daughters (especially 222 Rn dissolved in the external liquid argon), 232 Th and daughters, 60 Co and 40 K. The presence of these sources has been confirmed with a portable NaI γ-spectrometer inserted into the empty dewar. Bulk contaminations in the liquid argon of the chamber. In this case, both β and γ emitters are relevant. The most important contributions are from 222 Rn, 39 Ar and 85 Kr. 2 The signal from internal 222 Rn and its daughters can be monitored by counting the α decays of the isotopes 222 Rn, 218 Po and 214 Po in the energy region Other radioactive isotopes of Ar, as 37 Ar and 41 Ar, are short-lived (T 1/2 are 35 days and 109 min, respectively) and their cosmogenic production is negligible in the underground laboratory. The long-lived 42 Ar (T 1/2 = 32.9 y) is expected to be present in natural argon because of thermonuclear tests in the atmosphere; however, its concentration in nat Ar is negligible, < g/g (at 90% CL) [18], corresponding to less than 85 µbq/liter in liquid argon. 5

6 Table 1 Characteristics of the β decays of 39 Ar and 85 Kr. Isotope Half-life β end-point β mean energy (kev) (kev) 39 Ar 269 y Kr 10.8 y MeV. After a new filling with freshly produced liquid argon the chamber shows a total 222 Rn decay rate of the order of 1 2 Hz. Due to the 222 Rn half-life of 3.8 days, the observed rate decreases to few tens of events/day four weeks after the filling. Therefore, β decays from the Rn daughters 214 Pb and 214 Bi can be neglected, provided the measurement is performed a few weeks after the filling of the chamber. The counting rate due to the decay of 14 C (dissolved in the liquid argon or located in the surrounding plastics) is estimated to be much less than 50 mhz. Most of the 14 C events occur close to the chamber walls. The main characteristics of the 39 Ar and 85 Kr β decays are summarized in Table 1. Since both decays are classified as forbidden unique β transitions ( I π = 2 ), the β spectrum is not described by the usual Fermi function. For the present work, we assumed the β spectra from Ref. [5]. The Geant4-based simulation is used to generate normalized spectra s i (E) from the different radioisotopes, taking into account the energy resolution of the detector. The isotopes from the natural radioactive chains are treated independently. A χ 2 fit of the experimental spectrum F(E) is then performed in the energy range from 100 kev to 3 MeV with a linear combination of the single components, i.e. F(E) = i w i s i (E). (2) The coefficients w i are treated as free parameters and represent the counting rates induced by the single sources. In Figure 3 we show the experimental spectrum, superimposed with the output of the fit (i.e. w i s i (E)). The fit is satisfactory in all the energy range considered. The signals from the most important external γ-rays radioactivity sources and from internal contaminations are shown in Figure 4, as derived from the analysis of the experimental spectrum. 4 Discussion Figure 4 shows that the energy region 2 3 MeV is dominated by interactions of γ-rays from 232 Th daughters, the region MeV by γ-rays from 238 U daughters, and the region MeV by γ-rays from 60 Co and 40 K. Below 6

7 Fig. 3. Black histogram: background energy spectrum observed running the 2 liters detector deep underground at LNGS inside a 10 cm thick Pb shielding. The superimposed red histogram is the result of a fit with a Monte Carlo simulated signal (see text). 0.5 MeV the main contribution comes from β decays from internal contaminations of 39 Ar and 85 Kr; the two isotopes account for 65% of the total counting rate between 100 and 500 kev. The cosmogenically-originated 39 Ar contamination of nat Ar in the troposphere was measured in Ref. [3] to be (7.9 ± 0.3) g/g; the quoted error was statistical only 3. Since the liquid argon used for the experiment is produced from the atmospheric gas, a similar 39 Ar/ nat Ar ratio is expected to be present in our sample. 85 Kr is mainly produced as fission product of uranium and plutonium. Its abundance in the atmosphere is of the order of 1 Bq/m 3, corresponding to about g( 85 Kr)/g( nat Ar) in air. Nevertheless, the distillation procedure for the production of liquid argon substantially reduces the 85 Kr fraction. The residual 85 Kr in liquid argon may vary in different batches of liquid. In order to better show the Ar and Kr signals, Figure 5 displays the spectrum obtained from the experimental data after subtracting the fitted contribution from the other sources. The single 39 Ar and 85 Kr contributions can be disentangled from the different end-point energies, 565 kev and 687 kev respectively. Since 85 Kr and 39 Ar decays populate the same energy region of the 3 The 39 Ar/ nat Ar ratio was measured for dating purposes. The knowledge of the absolute 39 Ar specific activity was hence not necessary. 7

8 Fig. 4. The energy spectrum observed with the 2 liter detector can be reproduced by: (a) an external component dominated by interactions of γ-rays coming from U, Th, 60 Co and 40 K radioactivity of the materials surrounding the liquid argon; (b) an internal component dominated by 39 Ar and 85 Kr β contaminations inside the liquid argon. spectrum their estimates are anti-correlated, as displayed in Figure 6. The specific activity of 39 Ar in liquid argon resulting from the analysis is 1.41±0.02(stat)±0.11(syst) Bq/liter (1σ) (see below for the discussion of systematic uncertainties). The corresponding 39 Ar/ nat Ar mass ratio is (8.0±0.6) g/g, with errors summed in quadrature. The result is in excellent agreement with the atmospheric determination of Ref. [3]. From the fit it is also obtained that the 85 Kr activity in the sample under investigation is (0.16±0.13) Bq/liter (1σ). In a previous measurement performed with a sample of Argon % from Air Liquid, it was found a 85 Kr activity three times larger. This indicates that a non-negligible 85 Kr contamination may be found in commercial liquid argon samples. In this case, an additional fractional distillation could be required to reduce the radioactive background for the WARP experiment. The systematic uncertainties are summarized in Table 2. The dominating item is related to the energy calibration of the detector response: since the discrimination between 39 Ar and 85 Kr is based upon the β end-points, the fit result is sensitive to the energy calibration and resolution. The uncertainty on the energy calibration in the range of interest was evaluated to be 2% (1σ) from the meaurements with the γ-ray sources; the corresponding 39 Ar systematic error is 6.5%. The second important contribution is related to the uncertainty in the active volume of the chamber. The filling level can be determined with accuracy of about 1 mm, and the diameter of the teflon container with the 8

9 Fig. 5. Residual energy spectrum observed in the WARP chamber after subtraction of the expected contribution from external γ-rays. The spectrum is well described by a combination (dashed line) of β decays of 39 Ar and 85 Kr. Table 2 Systematic uncertainties on the 39 Ar specific activity. Item Relative error Absolute error (Bq/liter) Energy calibration ±6.5% ±0.092 Energy resolution ±1.3% ±0.018 Sensitive mass ±3.8% ±0.054 Total ± 0.11 reflector fixed on it is known with a precision of about 2 mm. This corresponds to an uncertainty on the sensitive mass of 3.8%. 5 Conclusions The best estimate of the 39 Ar specific activity in the liquid argon is (1.41 ± 0.11) Bq/liter, or (1.01 ± 0.08) Bq/kg of natural Ar, or (8.0 ± 0.6) g( 39 Ar)/g( nat Ar). The value is consistent with the previous determination by H. Loosli [3]. The uncertainty in our measurement is mainly due to systematics. The liquid argon sample under investigation shows a contamination of 85 Kr, 9

10 Fig. 6. Confidence regions (1σ, 2σ and 3σ) for the 39 Ar and 85 Kr specific activities. Systematic uncertainties have been added in quadrature. 0.16±0.13 Bq/liter (1σ). 6 Acknowledgments We wish to dedicate this work to the memory of our friend and colleague Nicola Ferrari, co-author of the paper, who prematurely passed away on July, We also thank Prof. H. Loosli for helpful communications concerning his 39 Ar paper. L. P. acknowledges the support by the EU FP6 project Ilias. A. M. S. has been in part supported by a grant of the President of the Polish Academy of Sciences, by the Ilias contract Nr. RII3-CT and by the MNiSW grant 1P03B References [1] R. Brunetti et al., New Astron. Rev. 49, 265 (2005). [2] WARP Collaboration, R. Brunetti et al., WARP: WIMP Argon Programme, Proposal for WARP to INFN, March 2004, available at 10

11 [3] H. H. Loosli, Earth and Planetary Science Letters, 63, 51 (1983). [4] J. A. Formaggio and C. J. Martoff, Ann. Rev. Nucl. Part. Sci. 54, 361 (2004). [5] LBNL Isotopes Project Nuclear Data Dissemination Home Page (available at [6] W. Kutschera et al., Nucl. Instr. Meth. B 92, 241 (1994). [7] P. Cennini et al., Nucl. Instr. Meth. A 356, 526 (1995). [8] C. Rubbia, CERN-EP/77-08 (1977). [9] P. Benetti et al., Nucl. Instr. and Meth. A 327, 203 (1993). [10] P. Aprili et al. [ICARUS Collaboration], CERN/SPSC (2002). [11] D. B. Cline et al., Nucl. Instr. Meth. A 503, 136 (2003). [12] M. G. Boulay and A. Hime, Astropart. Phys. 25, 179 (2006) [13] A. Rubbia, [arxiv::hep-ph/ ], Proceedings of the XI International Conference on Calorimetry in High Energy Physics (CALOR2004), Perugia, Italy (2004) [14] A. Rubbia, Journal of Physics, Conf. Series 39, 129 (2005) [arxiv::hep-ph/ ] [15] NIST Chemistry WebBook, NIST Standard Reference Database Number 69, June 2005, Eds. P. J. Linstrom and W. G. Mallard. [16] MACRO Collaboration, M. Ambrosio et al., Astropart. Phys. 10, 11 (1999). [17] Geant4 Collaboration, S. Agostinelli et al., Nucl. Instr. Meth. A 506, 250 (2003); Geant4 Collaboration, J. Allison et al., IEEE Trans. Nucl. Scie. 53, 270 (2006). [18] V. D. Ashitkov et al., Nucl. Instr. Meth. A 416, 179 (1998). 11