arxiv:0711.4583v1 [astro-ph] 28 Nov 2007 Temperature and Humidity Dependence of Air Fluorescence Yield measured by AIRFLY AIRFLY Collaboration: M. Ave a, M. Bohacova b, B. Buonomo c, N. Busca a, L. Cazon a, S.D. Chemerisov d, M.E. Conde d, R.A. Crowell d, P. Di Carlo e, C. Di Giulio f, M. Doubrava g, A. Esposito c, P. Facal h, F.J. Franchini d, J. Hörandel i,1, M. Hrabovsky b, M. Iarlori e, T.E. Kasprzyk d, B. Keilhauer i, H. Klages j, M. Kleifges k, S. Kuhlmann d, G. Mazzitelli c, L. Nozka b, A. Obermeier i, M. Palatka b, S. Petrera e, P. Privitera f,, J. Ridky b, V. Rizi e, G. Rodriguez f, F. Salamida e, P. Schovanek b, H. Spinka d, E. Strazzeri f, A. Ulrich l, Z.M. Yusof d, V. Vacek g, P. Valente m, V. Verzi f, T. Waldenmaier j a University of Chicago, Enrico Fermi Institute, 5640 S. Ellis Ave., Chicago, IL 60637, United States b Institute of Physics of the Academy of Sciences of the Czech Republic, Na Slovance 2, CZ-182 21 Praha 8, Czech Republic c Laboratori Nazionali di Frascati dell INFN, INFN, Sezione di Frascati, Via Enrico Fermi 40, Frascati, Rome 00044, Italy d Argonne National Laboratory, Argonne, IL 60439 United States e Dipartimento di Fisica dell Università de l Aquila and INFN, Via Vetoio, I-67010 Coppito, Aquila, Italy f Dipartimento di Fisica dell Università di Roma Tor Vergata and Sezione INFN, Via della Ricerca Scientifica, I-00133 Roma, Italy g Czech Technical University, Technicka 4, 16607 Praha 6, Czech Republik h Departamento de Física de Partículas, Campus Sur, Universidad, E-15782 Santiago de Compostela, Spain i Universität Karlsruhe (TH), Institut für Experimentelle Kernphysik (IEKP), Postfach 6980, D - 76128 Karlsruhe, Germany j Forschungszentrum Karlsruhe, Institut für Kernphysik, Postfach 3640, D - 76021 Karlsruhe, Germany k Forschungszentrum Karlsruhe, Institut für Prozessdatenverarbeitung und Elektronik, Postfach 3640, D - 76021 Karlsruhe, Germany l Physik Department E12, Technische Universität Muenchen, James Franck Str. 1, D-85748 Garching, Germany m Sezione INFN di Roma 1, Ple. A. Moro 2, I-00185 Roma, Italy Abstract The fluorescence detection of ultra high energy cosmic rays requires a detailed knowledge of the fluorescence light emission from nitrogen molecules over a wide range of atmospheric parameters, corresponding to altitudes typical of the cosmic ray shower development in the atmosphere. We have studied the temperature and humidity dependence of the fluorescence light spectrum excited by MeV electrons in air. Results for the 313.6 nm, 337.1 nm, 353.7 nm and 391.4 nm bands are reported in this paper. We found that the temperature and humidity dependence of the quenching process changes the fluorescence yield by a sizeable amount (up to 20%) and its effect must be included for a precise estimation of the energy of ultra high energy cosmic rays. Key words: Air Fluorescence Detection, Ultra High Energy Cosmic Rays, Nitrogen Collisional Quenching PACS:, 96.50.S-, 96.50.sb, 96.50.sd, 32.50.+d, 33.50.-j, 34.50.Fa, 34.50.Gb Preprint submitted to 5th Fluorescence Workshop, Madrid, 2007 2 February 2008
1. Introduction The detection of ultra high energy ( 10 18 ev) cosmic rays using nitrogen fluorescence light emission from extensive air showers (EAS) is a well established technique, used by the Fly s Eye [1], HiRes [2], and Pierre Auger Observatory [3] experiments, and for the Telescope Array [4] under construction. It has also been proposed for the satellite-based EUSO and OWL projects. Excitation of atmospheric nitrogen by EAS charged particles induces fluorescence emission, mostly in the wavelength range between 300 to 430 nm. Information on the longitudinal EAS development can be obtained by fluorescence telescopes by recording the light intensity as a function of time and incoming direction. However, the fluorescence light yield from EAS charged particles must be well known at each point within the shower, and corrections applied for atmospheric effects between the shower and the telescope for an accurate primary energy determination. Thus, the intensities of the fluorescence bands should be measured over a range of air pressure and temperature corresponding to altitudes up to about 16 km, the typical elevation of EAS development in the atmosphere. The presence of humidity will also affect the fluorescence yield, the effect being more important for satellite experiments which will detect showers over the oceans. The AIRFLY (AIR FLuorescence Yield) collaboration is pursuing an extensive measurement program of the fluorescence light yield with significantly improved precision with respect to previous experiments [5]. A precise measurement of the pressure dependence of the air fluorescence band intensities has been published [6], and is summarized in [7]. Measurements of the fluorescence yield dependence on the electron kinetic energy from kev to GeV are presented in [8]. The status of AIRFLY measurement of the absolute fluorscence yield is presented in [9]. The data reported here address the measurement of the fluorescence yield spectrum dependence on air temperature and humidity. The fluorescence yield of a band of wavelength λ at a given pressure p and temperature T can be written as [6]: Y air (λ, p, T) = Y air (337, p 0, T 0 ) I λ (p 0, T 0 ) corresponding author Email address: priviter@roma2.infn.it ( P. Privitera). 1 now at Department of Astrophysics, Radboud University Nijmegen, Nijmegen, The Netherlands 1 + p 1 + 0 p air (λ,t0) p p air (λ,t0) T H λ (T 0 ) T 0 H λ (T), (1) where Y air (337, p 0, T 0 ) is the absolute yield of the 337 nm band at pressure p 0 and temperature T 0 (in photons emitted per MeV of energy deposited), I λ (p 0, T 0 ) is the λ band intensity relative to the 337 nm band and p air (λ, T 0) is the band quenching reference pressure. H λ (T) has been introduced to take into account a possible temperature dependence of the collisional cross sections. The effect of humidity in the fluorescence yield can be introduced by substituting in Eq. (1): 1 p air 1 p air ( 1 p ) + 1 p H 2O p, (2) where is the water vapour partial pressure and p H 2O is the water vapour collisional quenching pressure. The experimental method used for the measurement of the temperature and humidity dependence is described in Section 2. The results obtained for the humidity and temperature dependence of four major fluorescence bands (313.6 nm, 337.1 nm, 353.7 nm and 391.4 nm) are presented in Sections 3 and 4, respectively. Implications for the energy determination of ultra high energy cosmic rays of the AIRFLY measurements reported in this paper are discussed in Section 5. Conclusions are given in Section 6. 2. Experimental method Measurements of both the humidity and temperature dependence of the fluorescence light yield were performed at the Argonne National Laboratory Chemistry Division s Van de Graaff (VdG) electron accelerator [10] [11]. It was operated in the DC current mode with typical beam currents of approximately 10 µa and nominal beam kinetic energy of 3.0 MeV. After a 30 bend in an electromagnet, the beam was focused near the exit from the accelerator vacuum, which consisted of a 35 mm diameter and 0.152 mm thick dura-aluminum window. For beam tuning, the beam spot of about 6 mm diameter was viewed with a quartz plate and mirror close to the exit window, and a camera. The quartz plate and mirror were removed during data collection. For each humidity or temperature measurement, spectra from 284-429 nm were recorded with a spectrograph. The beam current was also measured by a Faraday cup located in the accelerator vacuum. 2
2.1. Humidity Measurements The humidity measurements occurred in the same setup shown in Fig. 1 of Ref. [6]. The electron beam entered the pressure chamber through a 0.50 mm thick beryllium window, traversed 378 mm of gas, and left through a 0.1 mm thick aluminum window. Fluorescence light produced in the gas left the pressure chamber through a quartz window; it was collected by a 10 m long, 1.5 mm diameter pure silica core optical fiber, which brought the light to a spectrograph. The aluminum spherical mirror in [6] was not used for these measurements. The spectrograph was located behind a concrete block wall with additional lead shielding to protect it from radiation produced by the VdG. A remotely-controlled vacuum and gas handling system was used with the pressure chamber. A dry scroll vacuum pump [12] was used, and the chamber pressure [13] and relative humidity (rh) (Sensorika probe HTP-9511533 [14]) were measured at the pump-out port. A high purity dry gas (78.0% nitrogen, 21.0% oxygen, 1.0% argon) passed through a bubbler containing high purity water before entering the pressure chamber. Air fluorescence spectra were recorded by an Oriel MS257 TM spectrograph [15], described in detail in [6]. Calibrations with a mercury pencil lamp (Oriel no. 60635) and a quartz tungsten halogen lamp (Oriel no. 63350) were performed. In a fluorescence run the spectrograph collected data first in the wavelength range 284-369 nm and then in the range 344-429 nm. For each range, 50 spectra of 1 second integration time were taken in an automated sequence. The beam current from the Faraday cup was recorded before and after each sequence of spectra collected. 2.2. Temperature Measurements A special temperature chamber was constructed for the temperature measurements; see Fig. 1. It consisted of a stainless steel pipe with two horizontal arms and one vertical arm. The electron beam entered and exited the chamber through 100 µm thick aluminum windows. One horizontal arm had a stainless steel plate blankoff, and the other had the same quartz window as the pressure chamber. The gas length traversed by the electron beam in the temperature chamber was approximately 220 mm. The vertical arm of the temperature chamber con- Fig. 1. The temperature chamber at the Van de Graaff: a) the chamber on the beam line as it appear before mounting the polysterene box, b) the chamber inside the polystyrene box, with the protective drum and pipe also in place. tained feed-thrus for five sensors: two relative humidity (moisture) sensors (Sensorika HS 2Ta [14]), a pressure sensor (Sensortechnics BT/PTU7000 [16]), and two temperature sensors (Sensorika Pt 10000 [14]). The humidity and pressure sensors were located near the vacuum flange with the feed-thrus, about 120 mm above the beam line. One temperature sensor was near beam height, but offset perpendicular to the beam and away from the quartz window about 3 cm, and the other was approximately 70 mm above beam height. Signals from all sensors were collected in a portable data acquisition system including a calibrated signal converter, and then sent to a PC to monitor the chamber conditions. The temperature chamber body was enclosed in a polystyrene box with aluminum supporting structure. The box was filled with dry ice (CO 2 ) for cooling of the chamber down to -40 C. A 2 m long strip heater (48 V, 325 W) wrapped around the chamber body was used to regulate the temperature inside. In order to avoid buildup of frost on the quartz window, a protective drum and pipe were added to the temperature chamber. These covered the window and end of the optical fiber. Boil-off gas from a liquid nitrogen dewar was continuously fed through the protective drum and pipe system. The same remotely-controlled vacuum and gas handling system, silica core optical fiber, Faraday cup, and spectrograph were used as for the humidity measurements. 3. Temperature dependence Fluorescence yield measurements were performed in the temperature range between 240 K and 310 K in dry air. After allowing for temperature to stabilize in the chamber, a fluorescence spectrum was taken. 3
rh (%) rh (%) (313.6) (arb. units) (353.7) (arb. units) 0.095 0.09 0.085 0.08 0.075 0.07 0.065 0.145 0.14 0.135 0.13 0.125 0.12 0.115 (337.1) (arb. units) (391.4) (arb. units) 0.7 0.68 0.66 0.64 0.62 0.6 0.58 0.17 0.165 0.16 0.155 0.15 0.145 0.14 (313.6) (arb. units) (353.7) (arb. units) 0.046 0.044 0.042 0.04 0.038 0.036 0.034 0.032 0.03 0.028 rh (%) 0.075 0.07 0.065 0.06 0.055 0.05 (337.1) (arb. units) (391.4) (arb. units) 0.36 0.34 0.32 0.3 0.28 0.26 0.24 0.22 rh (%) 0.095 0.09 0.085 0.08 0.075 0.07 0.065 0.06 0.055 Fig. 2. Temperature dependence of the air fluorescence signal for the 313.6 nm, 337.1 nm, 353.7 nm and 391.4 nm bands The fluorescence band intensities were estimated by integrating the CCD counts in a wavelength interval around each band. A detailed account of the procedure can be found in [6]. For each fluorescence band λ, the corresponding fluorescence signal (λ) was obtained from the ratio of the number of counts in the band interval to the beam current. The temperature dependence of the fluorescence signal of the 313.6 nm, 337.1 nm, 353.7 nm and 391.4 nm bands is shown in Fig. 2, together with the result of a fit to the data with the following ansantz in Eq. (1): ( ) αλ H λ (T) T H λ (T 0 ) =. (3) T 0 The fitted values of α λ are reported in Table 1. Notice that α λ = 0 indicates that collisional cross sections have no temperature dependence, which has been so far assumed in the application of fluorescence yield measurements to fluorescence detection of ultra high energy cosmic rays. The AIRFLY measurements show a sizable deviation of the temperature dependence of the fluorescence yield from this usual assumption. Also, the value of α λ appears to depend on the wavelength band, with the 391.4 nm band being significantly different from the others. 4. Humidity dependence Fluorescence yield measurements were performed for relative humidity ranging from 0 to almost 100% corresponding to water vapour partial pressures Fig. 3. Air fluorescence signal as a function of water vapour partial pressure (relative humidity on the top axis) for the 313.6 nm, 337.1 nm, 353.7 nm and 391.4 nm bands. up to about 25 hpa. After allowing for humidity to stabilize in the chamber, a fluorescence spectrum was taken. For each fluorescence band λ, the corresponding fluorescence signal (λ) was obtained from the ratio of the number of counts in the band interval to the beam current. The beam current was measured several times during the run with a Faraday cup placed inside the beam pipe. The measured fluorescence signal as a function of the water vapour partial pressure for the 313.6 nm, 337.1 nm, 353.7 nm and 391.4 nm bands is shown in Fig. 3, together with the result of a fit to the data using Eq. (2). The fitted values of the water vapour collisional quenching pressure p H 2O are reported in Table 1. Notice that the quenching effect of water vapour is not negligible, resulting in about 20% decrease of the fluorescence yield for rh=100%. λ (nm) α λ p H 2 O 313.6 0.09 ± 0.10 1.21 ± 0.13 337.1 0.36 ± 0.08 1.28 ± 0.08 353.7 0.21 ± 0.09 1.27 ± 0.12 391.4 0.80 ± 0.09 0.33 ± 0.03 Table 1 Summary of measured temperature and humidity dependence parameters for a selected group of air fluorescence bands.. 4
Y T / Y 1.1 Y H / Y 1.1 1 1 0.9 0.9 0.8 0.8 0.7 30 altitude (km) 0.7 0 2 4 6 8 10 altitude (km) Fig. 4. Ratio of fluorescence yield with measured α λ to the one with α λ = 0: dashed line 313.6 nm; full line 337.1 nm; dotted line 353.7 nm; dashed-dotted 391.4 nm 5. Implications for the fluorescence detection of ultra high energy cosmic rays Experiments employing the fluorescence technique for the detection of ultra high energy cosmic rays have so far assumed α λ = = 0 in the treatment of fluorescence yield. The AIRFLY results presented in Sections 3 and 4 do not support this assumption. In fact, neglecting temperature and humidity effects on the fluorescence yield introduces a sizeable bias in the energy estimation of EAS. In Fig. 4, the ratio of the fluorescence yield with α λ from Table 1 to the one with α λ = 0 is shown as a function of altitude. All yields were normalized to have the same value at ground level, and the U.S. 1976 Standard Atmosphere [17] was used to calculate the pressure and temperature at a given altitude. It appears from Fig. 4 that neglecting the temperature dependence results in an overestimation of the fluorescence yield by an amount going up to 20% for the 391.4 nm band. We also investigated the bias introduced by neglecting the humidity effect on the fluorescence yield. We expect the effect to be more important for satellite experiments which will detect showers over the oceans. In Fig. 5, the ratio of the fluorescence yield with p H 2O from Table 1 to the one with = 0 is shown as a function of altitude, assuming a humidity profile over oceans as given in [18]. The effect of humidity is sizeable up to about 5 km altitude. A full estimate of the effect of neglecting temperature and humidity in the fluorescence yield on Fig. 5. Ratio of fluorescence yield with measured p H 2 O to the one with = 0: dashed line 313.6 nm; full line 337.1 nm; dotted line 353.7 nm; dashed-dotted 391.4 nm the EAS energy measurement would need a detailed study, including the EAS energy deposit as a function of altitude, the wavelength dependent attenuation of the emitted fluorescence light in the atmosphere and the sensitivity of the fluorescence detectors. This goes beyond the scope of this paper. On the other hand, the experimental accuracy on the energy estimation of EAS with the fluorescence detection technique asks for a systematic uncertainty on the knowledge of the fluorescence yield better than 10%. Thus, biases of the order of 10%, as seen in Figs. 4 and 5, should not be neglected. 6. Conclusions The measurement of the EAS energy with the fluorescence detection technique requires a detailed knowledge of the air fluorescence emission as a function of pressure, temperature and humidity. The AIRFLY experiment has performed precise measurements of the fluorescence light spectrum excited by MeV electrons in dry air. The relative intensities of 34 fluorescence bands as well as their pressure dependence have been fully reported [6], and are summarized in a contribution to this workshop [7]. In this paper, we reported measurements of the temperature and humidity dependence of the 313.6 nm, 337.1 nm, 353.7 nm and 391.4 nm fluorescence bands in air. The temperature dependent measurements of the fluorescence yield, performed in the range between 240 K and 310 K, are the first reported in the literature. Our data show that 5
collisional cross sections are temperature dependent, contrary to the assumption so far made in the UHECR community. We also measured the quenching effect of water vapour, resulting in about 20% decrease of the fluorescence yield for rh=100%. The effect of the temperature and humidity dependence of the fluorescence yield on the energy estimation of EAS was studied, and we found that it should not be neglected. A full analysis of the AIRFLY data, with the measurement of the temperature and humidity dependence of 34 fluorescence bands, is in preparation. [12] Varian, Inc., 121 Hartwell Avenue, Lexington, MA 02421, model SH01001UNIV. [13] Oerlikon Leybold Vacuum GmbH, Bonner Strasse 498, Cologne D-50968, Germany, model Ceravac CTR 90. [14] Sensorika s.r.o., V. Záti sí 74/4, 147 00 Praha - Hodkovi cky, Czech Republic, [15] Spectra-Physics, 150 Long Beach Boulevard, Stratford, Ct 06615, USA. [16] Sensortechnics GmbH, Boschstr. 10, 82178 Puchheim, Germany [17] U.S. Standard Atmosphere 1976, U.S. Government Printing Office (Washington, D.C., 1976). [18] Liu et al, J. Climate 4 (1991) 1023. 7. Acknowledgments We thank the staff of Argonne National Laboratory for their support. This work was also supported by the grant of MSMT CR LC 527 and 1M06002 and ASCR grants AV0Z10100502 and AV0Z10100522. A. Obermeier and J. R. Hörandel acknowledge the support of VIHKOS, which made the participation at the measurement campaigns possible. References [1] R.M. Baltrusaitus et al., Nucl. Instrum. Meth. Phys. Res. A 240 (1985) 410. [2] T. Abu-Zayyad et al., Nucl Instrum. Meth. Phys. Res. A 450 (2000) 253. [3] J. Abraham et al., Nucl Instrum. Meth. Phys. Res. A 523 (2004) 50. [4] M. Fukushima, Prog. Theor. Phys. Suppl. 151 (2003) 206. [5] A. N. Bunner, PhD thesis, Cornell University, 1967; G. Davidson and R. O Neil, J. Chem. Phys. 41 (1964) 3946; F. Kakimoto et al., Nucl. Instrum. Meth. Phys. Res. A 372 (1996) 527; M. Nagano et al., Astropart. Phys. 20 (2003), 293; 22 (2004), 235; J. W. Belz et al., Astropart. Phys. 25 (2006) 129; P. Colin et al., Astropart. Phys. 27 (2007) 317. [6] AIRFLY Collaboration, M. Ave et al, Astropart. Phys. 28 (2007) 41. [7] AIRFLY Collaboration, M. Ave et al, Proceedings of this Workshop. [8] AIRFLY Collaboration, M. Ave et al, Proceedings of this Workshop. [9] AIRFLY Collaboration, M. Ave et al, Proceedings of this Workshop. [10] High Voltage Engineering Corporation model KS Van de Graaff. [11] An Eimac Y646 cathode electron source was used from Communications and Power Industries, Eimac Division, 607 Hansen Way, Palo Alto CA 94304 USA. 6