32ND INTERNATIONAL COSMIC RAY CONFERENCE,BEIJING 2011 Measurement of Deuterium and 3He component in cosmic rays with PAMELA experiment VALERIO FORMATO1,2 ON BEHALF OF THE PAMELA COLLABORATION 1Universita` degli studi di Trieste 2INFN - Sezione di Trieste [email protected] DOI: 10.7529/ICRC2011/V06/0842 Abstract: PAMELA is a satellite borne experiment designed to study with great accuracy cosmic rays of galactic, solar, and trapped nature, with particular focus on the antimatter component. The detector consists of a permanent magnet spectrometer core to provide rigidity and charge sign information, a Time-of-Flight system for velocity and charge in- formation, a Silicon-Tungsten calorimeter and a Neutron detector for lepton/hadron identification. The beta and rigidity information allow to identify isotopes for Z =1and Z =2particles in the energy range 100 MeV/n to 1 GeV/n. In this work we will present observation of cosmic rays D and He isotopes at the 23rd solar minimum from 2006 to 2008. Specifically we will discuss the ratios D/H and 3He/4He as a function of time and compare the results with existing measurements and models. Keywords: 1 Introduction in CR, see [4]. The core of the instrument is a perma- nent magnet with six planes of silicon microstrip detectors Hydrogen and helium isotopes in cosmic rays are gener- for tracking particles. A three-scintillators system provides ally believed to be of secondary origin, resulting mainly trigger, charge and time of flight information. A silicon- from the nuclear interactions of primary cosmic-ray pro- tungsten calorimeter is used for hadron/lepton separation tons and 4He with the interstellar medium. A precise mea- with a shower tail catcher and a neutron detector at the bot- surement of the isotopic composition of hydrogen and he- tom of the apparatus to help increase this separation. An lium could provide important information about the propa- anticounter system is used to reject spurious events in the gation of cosmic rays in the interstellar space, since they of- off-line phase. The readout electronics, the interfaces with fer better statistics than other secondaries. This can be cru- the CPU and all primary and secondary power supplies are cial when analyzing positron and antiproton spectra, which housed around the detectors. All systems (power supply, is the primary goal of several current cosmic-ray experi- readout boards etc.) are redundant with the exception of ments [1, 2]. the CPU which is more tolerant to failures. The apparatus The PAMELA experiment has been measuring since 2006 is enclosed in a pressurized container located on one side both hydrogen and helium isotopes over a wide range of of the Resurs-DK1 satellite. Total weight of PAMELA is time with reduced instrumental uncertainty and with no en- 470 kg; power consumption is 355 W. A more detailed de- vironmental systematics, such as - for instance - the resid- scription of the instruments and the data handling can be ual atmosphere present in balloon experiments. found in [3]. In this paper we present the isotopic composition of helium in the energy range from 0.1 to about 1 GeV nucleon-1 (cor- 3 Analysis responding to rigidities from 800 MV to 2.5 GV) ob- ∼ ∼ tained from the flight data recorded until December 2008. The high inclination (70◦) orbit of the Resurs-DK1 satellite and the long mission duration allow PAMELA to study low 2 Pamela detector energy particles (down to 50 MeV) investigating several aspects of the solar and terrestrial environment. PAMELA [3] is constituted by a number of redundant de- To select the primary (galactic) component we evaluated tectors (see Fig. 1) capable of identifying particles pro- the local geomagnetic cutoff G in the Stormer¨ approxi- 2 viding charge, mass, rigidity (ρ = p/Ze) and velocity mation [5]. The value of G =14.9/L - valid for verti- (β = v/c) over a very wide energy range. Its main sci- cally incident particles - is estimated using the IGRF mag- entific objective is the study of the antimatter component netic field model along the orbit; from this the McIlwain Vol. 6, 79 V. F ORMATO et al. DEUTERIUM AND 3 HEIN COSMIC RAYS WITH PAMELA EXPERIMENT ! 30 104 $% & 25 103 20 "# 102 15 4He 10 dE/dx (arbitrary units) 10 3He $% ' 2H 1 5 1H 0 -1 -1 10 10 1 ρ (GV) Figure 2: Energy loss in tracker (mean in all planes hit) vs tracker rigidity for positively charged particles. The Pro- $% ( ton and Helium bands are clearly visible. The black lines represent cuts used to select H and He nuclei. Hydrogen β 1.2 1 102 1 0.8 H Figure 1: A schematic overview of the PAMELA satellite 0.6 experiment. The experiment stands 1.3 m high and, from 10 ∼ 0.4 2H top to bottom, consists of a time-of-flight (ToF) system (S1, 0.2 S2, S3 scintillator planes), an anticoincidence shield sys- 0 1 0 0.5 1 1.5 2 2.5 tem, a permanent magnet spectrometer (the magnetic field ρ (GV) runs in the y-direction), a silicon-tungsten electromagnetic Helium β calorimeter, a shower tail scintillator (S4) and a neutron de- 1.2 tector. 1 102 0.8 L shell is calculated [6]. Particles where selected requiring 0.6 3He 10 ρ>1.3 G to remove any effect due to directionality in the 4 0.4 He detector and Earth’s penumbral regions. 0.2 For the measurement of ratios as functions of kinetic en- 0 1 0 0.5 1 1.5 2 2.5 ergy the geomagnetic cutoff cut is applied considering the ρ (GV) same geographical region for the two isotopes requiring ρ∗(E) > 1.3G, where ρ∗(E) is the lowest rigidity between Figure 3: Mass separation for Z =1(top) and Z =2 the two isotopes having kinetic energy E. (bottome) particles with the β vs. rigidity method. The We have selected events that do not produce secondary par- black lines represent the theoretical curves from Eq. 1 ticles by requiring a single track fitted within the spectrom- eter fiducial acceptance and a maximum of one paddle hit relativistic also at low rigidities and the background of pion in the two top scintillators of the ToF system. The hits in and secondary particles. The black lines show the energy the two scintillators must match the extrapolated trajectory dependent cuts used to select the helium sample. from the tracker. Albedo particles are rejected requiring dE/dx β>0. Cuts in the energy loss ( ) vs rigidity remove positrons, pions and particles with Z>2 as shown Fig 2. 3.1 Hydrogen and Helium selection The measurement of the average energy released in the 3.2 Isotope separation tracker planes for a given event at a given rigidity can be 2 1 3 4 used to discriminate between different particles since en- To separate H from H and He from He we have taken 1/β ergy loss of a charged particle through matter follows the the distribution (described by a Gauss function) at Bethe Block formula, dE/dx Z2/β2 (neglecting log- fixed rigidity. As shown in Fig. 3 the different masses ∝ arithmic terms). In Fig. 2 is shown the energy loss in the of the isotopes guarantee a separation (decreasing with en- 1/β tracker as a function of the particle rigidity. The top band is ergy) in the distribution according to due to helium nuclei which have energy loss in the tracker 2 1 A2 m2 Z =4times the protons, identified in the central band. = 1+ p (1) The band at the bottom left of Fig. 2 is due to positrons, β Z2 ρ2 Vol. 6, 80 32ND INTERNATIONAL COSMIC RAY CONFERENCE,BEIJING 2011 no longer valid when measuring kinetic energy isotopic ra- tios since, due to the different mass, at a given kinetic en- 1 Counts 4 1.25 GV < ρ < 1.32 GV Counts 4 H 1.89 GV < ρ < 2.00 GV ergy value correspond different rigidities for the isotopes; 10 1H 10 thus the selection efficiency has been corrected by a factor 103 103 equal to the ratio of the efficiencies estimated by the analy- 2H sis on proton and helium absolute fluxes. This correction is 2 2 -1 10 2H 10 less than 10% below 300 MeV nucleon and less than 1% above 300 MeV nucleon-1. 10 10 1 1 3.3 Contribution of Secondary Particles - Top of the Payload Correction 10-1 10-1 1 1.5 2 2.5 3 3.5 4 1 1.5 2 2.5 3 3.5 4 1/β 1/β Helium and Hydrogen nuclei may be lost due to hadronic interactions in the 2 mm Al thick pressurized container and the top scintillator detectors; the correction to the ratio due 4He 4He 2 3 Counts Counts to particles lost and the contribution of H and He from 3 1.51 GV < ρ < 1.62 GV 3 1.74 GV < ρ < 1.86 GV 4 10 10 3He He inelastic scattering has been estimated to be less than 3He 1% and has been included in the systematic uncertainty of the measurement. 102 102 10 10 4 Results The results on the 2H/1H and 3He/4He ratios will be pre- 1 1 sented during the conference. 1 1.2 1.4 1.6 1.8 2 2.2 2.4 1 1.2 1.4 1.6 1.8 2 2.2 2.4 1/β 1/β 5 Conclusions Figure 4: Examples of isotopic separation between 2H and 1H(top), and between 3He and 4He (bottom) at fixed rigid- In this work we have shown the capability of PAMELA to 2 1 3 4 ity range.