ABSTRACT
The propagation of cosmic rays inside our galaxy plays a fundamental role in shaping their injection spectra into those observed at Earth. One of the best tools to investigate this issue is the ratio of fluxes for secondary and primary species. The boron-to-carbon (B/C) ratio, in particular, is a sensitive probe to investigate propagation mechanisms. This paper presents new measurements of the absolute fluxes of boron and carbon nuclei as well as the B/C ratio from the PAMELA space experiment. The results span the range 0.44–129 GeV/n in kinetic energy for data taken in the period 2006 July to 2008 March.
Export citation and abstract BibTeX RIS
1. INTRODUCTION
Propagation in the interstellar medium (ISM) significantly affects the spectrum of galactic cosmic rays. After being accelerated by high-energy astrophysical processes such as supernovae explosions, cosmic rays are injected into the interstellar space, propagate through it and eventually reach the Earth where they are detected. The multitude of physical processes that cosmic rays undergo during propagation (e.g., diffusion, spallation, emission of synchrotron radiation, etc.) shape the injection spectra and chemical composition into the observed values. A detailed knowledge of these processes is therefore needed in order to interpret the experimental data in terms of source parameters, or in estimating the expected background when searching for contributions from new sources.
There is still a relatively high degree of uncertainty regarding the physical processes relevant to propagation of cosmic rays and the impact of experimental uncertainties on the determination of propagation parameters (see Maurin et al. 2010 and references therein). The propagation is usually modeled in terms of a diffusive transport equation (Ptuskin 2012). The equation contains terms that account for diffusion in the irregular galactic magnetic field, convection due to the galactic wind, energy losses, re-acceleration (modeled as diffusion in momentum space), spallation and radioactive decay, and source terms. Some parameters of the equation are simply related to directly measurable quantities unrelated to cosmic rays, and thus they can be obtained from independent measurements (e.g., the density of atomic hydrogen in the ISM, which is needed in order to estimate the spallation rate, can be measured by means of 21 cm radio surveys). Other parameters are obtained by fitting distributions derived from numerical propagation models like GALPROP (Strong & Moskalenko 1998; Vladimirov 2012) or DRAGON (Gaggero et al. 2013) to direct cosmic ray measurements.
In order to test and tune the propagation models, a particularly useful measurable quantity is the secondary to primary flux ratio. Primary nuclei are those accelerated by cosmic ray sources such as supernova remnants, whereas secondaries are those produced in interactions of primaries with the ISM during propagation. The boron to carbon flux ratio (B/C) has been widely studied. Since boron is produced in negligible quantities by stellar nucleosynthesis processes (Bethe 1939), almost all of the observed boron is believed to be from spallation reactions of CNO primaries on atomic and molecular H and He present in the ISM. The B/C flux ratio is therefore a clean and direct probe of propagation mechanisms, and it is considered as the "standard tool" for studying propagation models (Strong et al. 2007; Obermeier et al. 2012).
The B/C flux ratio, as well as the absolute boron and carbon fluxes, have been measured by balloon-borne (Freier et al. 1959; Panov et al. 2008; Ahn et al. 2008; Obermeier et al. 2011) and by space-based experiments (Engelmann et al. 1990; Swordy et al. 1990; Webber et al. 2002; Aguilar et al. 2010; Lave et al. 2013; Oliva et al. 2013), with different techniques and spanning various energy ranges from about 80 MeV/n up to a few TeV/n. Even if the spread in the measurements and their associated errors makes it difficult to clearly discriminate between the various models or to tightly constrain model parameters, there is a general consensus about several points. The relative abundance of the light elements Li, Be, and B in cosmic rays is significantly higher than in the solar system (de Nolfo et al. 2006). This supports the idea of creation by spallation reactions in ISM. The B/C flux ratio has a peak value at ∼1 GeV/n, which can favor a model with distributed stochastic re-acceleration (Letaw et al. 1993). The B/C flux ratio decreases at high energies and its shape, in diffusive models, is mainly determined by the energy dependence of the diffusion coefficient (Castellina & Donato 2005).
In this paper, a new set of measurements of boron and carbon fluxes as well as the B/C flux ratio obtained with the PAMELA instrument in the kinetic energy range 0.44–129 GeV/n during the solar minimum period spanning from 2006 July to 2008 March are presented. The study of solar modulation effects on the low-energy component of the spectra over a longer time interval will be the subject of a future publication. After a brief description of the PAMELA detector system, the analysis techniques and an evaluation of systematic uncertainties are presented, followed by a discussion of the results.
2. THE PAMELA DETECTOR
A schematic view of the PAMELA detector system (Picozza et al. 2007) is shown in Figure 1. The design was chosen to meet the main scientific goal of precisely measuring the light components of the cosmic ray spectrum in the energy range starting from tens of MeV up to 1 TeV (depending on particle species), with a particular focus on antimatter. To this end, the design is optimized for |Z| = 1 particles and to provide a high lepton–hadron discrimination power. The core of the instrument is a magnetic spectrometer (Adriani et al. 2007) made by six double-sided silicon microstrip tracking layers placed in the bore of a permanent magnet. The read-out pitch of the silicon sensors is 51 μm in the X (bending) view and 66.5 μm in the Y view. The spectrometer provides information about the magnetic rigidity ρ = pc/(Ze) of the particle (where p and Z are the particle momentum and the electric charge, respectively). Six layers of plastic scintillator paddles arranged in three X–Y planes (S1, S2, and S3 in Figure 1) placed above and below the magnetic cavity constitute the time-of-flight (TOF) system (Barbarino et al. 2008; Osteria & Russo 2008). The flight time of particles is measured with a time resolution of 250 ps for |Z| = 1 particles and about 70 ps for boron and carbon nuclei (Campana et al. 2009). This allows albedo particles to be rejected and, in combination with the track length information obtained from the tracking system, precise measurement of the particle velocity, β = v/c. The TOF scintillators can identify the absolute particle charge up to oxygen by means of six independent ionization measurements. The tracking system and the upper TOF system are shielded by an anticoincidence system (AC; Pearce et al. 2003) made of plastic scintillators and arranged in three sections (CARD, CAT, and CAS in Figure 1), which allows spurious triggers generated by secondary particles to be rejected during offline data analysis. A sampling electromagnetic calorimeter (Boezio et al. 2002; Bonvicini et al. 2009) is placed below S3. It consists of 22 modules, each comprising a tungsten converter layer placed between two layers equipped with single-sided silicon strip detectors with orthogonal read-out strips. The total depth of the calorimeter is 16.3X0, while the readout pitch of the strips is 2.44 mm. The calorimeter measures the energy of electrons and positrons, and gives a lepton/hadron rejection power of ∼105 by means of topological shower analysis, thanks to its fine lateral and longitudinal segmentation. A tail-catcher scintillating detector (S4) and a neutron detector placed below the calorimeter help to further improve the rejection power.
Figure 1. Schematic view of the PAMELA apparatus.
Download figure:
Standard image High-resolution imageThe geometric factor of the apparatus, defined by the magnetic cavity, is energy dependent because of the track curvature induced by the magnetic field, and increases as the energy of the particle increases. However, for rigidities above 1 GV it varies only by a few per mil, reaching the value of 21.6 cm2 sr at the highest rigidity.
The PAMELA apparatus was launched on 2006 June 15, and has been continuously taking data since then. It is hosted as a piggyback payload on the Russian satellite Resurs-DK1, which executes a 70° semi-polar orbit. The orbit was elliptical with variable height between 350 and 620 km up to 2010, after which it was converted to the current circular orbit with height about 600 km.
3. DATA ANALYSIS
3.1. Data Processing
The event reconstruction routines require silicon strips to be gathered into clusters. A "seed" strip is defined as a strip with a signal-to-noise ratio (S/N) greater than 7; it is grouped with its neighboring "signal" strips with S/N > 4 to form a cluster. For each cluster an estimate of the particle impact point is obtained by means of an analog position finding algorithm (Adriani et al. 2007). The original reconstruction routines were conceived and tuned to deal with |Z| ∼ 1 particles. However the higher ionization energy losses of boron and carbon in the silicon layers of the tracking system saturate the front-end electronics, leading to a degradation of the performance of position finding with respect to |Z| ∼ 1 particles. A different position finding algorithm has thus been implemented for this saturation regime. For each cluster of silicon strips, the saturated strips have been treated as if read-out system was digital, and the impact point has been evaluated as the geometric center. The associated spatial resolution can be approximated as the readout pitch over , which translates to ∼14 μm for the X (bending) view and ∼19 μm for the Y view. The associated maximum detectable rigidity (MDR24) is ∼250 GV.
Prior to event reconstruction, the clusters with an associated energy release less than 5 MIP25 have been removed. This helps to eliminate clusters associated with delta rays and light secondary particles, e.g., backscattered particles from the calorimeter. There is a twofold effect: the tracking efficiency is increased since the tracking algorithm has less clusters to deal with, and the energy dependence of the tracking efficiency is reduced at high energies by removing backscattering clusters, which are mainly produced by high-energy primaries interacting in the calorimeter.
3.2. Event Selection
In order to be able to reliably measure the magnetic rigidity, events with a single track in the spectrometer containing at least four hits in the X view and three hits in the Y view have been selected. A good χ2 value for the fitted track was required. The χ2 distribution is energy dependent and thus the selection criterion has been calibrated in order to obtain a constant efficiency of about 90% over the whole energy range, in particular at low energies where multiple scattering leads to generally higher χ2 values. Reconstructed tracks were required to lie entirely inside a fiducial volume with bounding surfaces 0.15 cm from the magnet walls. Galactic events were selected by imposing that the lower edge of the rigidity bin to which the event belongs exceeds the critical rigidity, ρc, defined as 1.3 times the cutoff rigidity ρSVC computed in the Störmer vertical approximation (Shea et al. 1987) as ρSVC = 14.9/L2, where L is the McIlwain L-shell parameter (McIlwain 1961) obtained by using the Resurs-DK1 orbital information and the IGRF magnetic field model (MacMillan & Maus 2005). The South Atlantic Anomaly region has been included in the analysis. Reconstructed particle trajectories were required to be down-going according to the TOF. No selections on the hit pattern in the TOF paddles or AC were made, since this can lead to very low efficiencies due to the production of delta rays in the aluminum dome of the pressurized vessel in which PAMELA is hosted. This introduces a contamination from secondaries produced in hadronic interactions of primaries in the dome. This effect has been accounted for using Monte Carlo simulations.
Boron and carbon events have been selected by means of ionization energy losses in the TOF system. Charge consistency has been required between S1226 and 〈S2〉 and 〈S3〉 (the arithmetic mean of the ionizations for the two layers constituting S2 and S3, respectively). Requiring charge consistency above and below the tracking system rejected events interacting in the silicon layers. The selection bands as functions of the rigidity measured by the spectrometer are shown in Figure 2.
Figure 2. Charge selection bands for S12, 〈S2〉, and 〈S3〉 as a function of rigidity. The red vertical dotted lines denote the upper and lower rigidity limits of this analysis. The absence of relativistic protons in this sample is due to the 5 MIP cluster selection described in Section 3.1.
Download figure:
Standard image High-resolution imageIn order to assess the presence of possible contamination in the selected samples, the above selection cuts have been applied to boron and carbon samples independently selected by means of S11 (the upper layer of S1) and the first silicon layer of the calorimeter. The probabilities of misidentifying a carbon nucleus as boron and vice versa are about 3 × 10−4 and 10−3, respectively, over the whole energy range considered in this analysis. Stricter analysis criteria were imposed by narrowing the selection bands. When properly corrected by the selection efficiency (see Section 3.3), the event counts showed no statistically significant deviation from that obtained using the standard selection. The contamination is therefore assumed to be negligible. Selected events have been binned according to the rigidity measured by the magnetic spectrometer.
3.3. Efficiencies
The tracking efficiency has been evaluated with flight data and Monte Carlo simulations using a methodology similar to that described in (Adriani et al. 2013). Two samples of boron and carbon were selected by means of a β dependent requirement on ionization energy losses in the TOF system. Fiducial containment was verified using calorimeter information. Firstly, non-interacting events penetrating deeply into the calorimeter were identified, and a straight track fitted. Then, the rigidity of the nucleus was derived from the β measured by the TOF and used to back-propagate the track through the spectrometer magnetic field up to the top of the apparatus. The containment criteria were applied to this back-extrapolated track. The tracking efficiency was determined for this sample of non-interacting nuclei as a function of the rigidity derived from β. The 70 ps resolution of the TOF system for carbon leads to Δβ/β ∼ 2% at β = 0.9 (Campana et al. 2009). Bin folding effects on the efficiency have therefore been neglected.
Due to the calorimeter selection criteria described above, the efficiency is measured for a non-isotropically distributed sample, while the fluxes impinging on PAMELA are isotropic. Moreover, a possible energy dependence of the efficiency at relativistic energies cannot be accounted for by an efficiency measured as a function of β. To account for these effects, a simulation of the PAMELA apparatus based on GEANT4 (Agostinelli et al. 2003; Allison et al. 2006) has been used to estimate the isotropic, rigidity-dependent tracking efficiency which is subsequently divided by a Monte Carlo efficiency obtained using the same procedure as the experimental efficiency. The resulting ratio, which has an almost constant value of about 0.97, has been used as the correction factor for the experimental efficiency. The constancy of the ratio results from an isotropic efficiency that is also almost constant above 10 GV because of the data processing procedures described in Section 3.1.
The efficiencies for the selection of down-going particles and for charge selection have been estimated using flight data exclusively. The down-going requirement is 100% efficient due to the 70 ps resolution of the TOF system. To evaluate the charge selection efficiency, the redundancy of the PAMELA subdetectors has been exploited. Two samples of boron and carbon have been tagged requiring charge consistency on S11 and on the first silicon layer of the calorimeter. These two detectors are placed at the two extrema of the apparatus, so this selection rejects interactions which change the reconstructed charge of the incident particle. The resulting efficiencies have a peak value of ∼75% at 3 GV and then decrease at high energies toward an almost constant value of about 50% for boron and 60% for carbon above some tens of GV.
The tracking and the charge selection efficiencies are shown in Figure 3 together with the total selection efficiency.
Figure 3. Selection efficiencies as functions of rigidity. The dashed line is a fit of the charge selection efficiency above 3 GV with a power law at low rigidities and a constant value at high rigidities. The slope, the break point and the normalization are free parameters of the fit. The fitted charge selection efficiency is used to compute the total efficiency above 3 GeV/n (about 7.6 GV for C and 10B and 8.4 GV for 11B) in order to smooth the statistical fluctuations.
Download figure:
Standard image High-resolution imageThe measurement of the charge selection efficiency sets the lower rigidity limit for fluxes to 2 GV, corresponding to about 0.44 GeV/n for 10B and 12C. Below this threshold charge confusion in the calorimeter selection becomes too large to be able to reliably tag pure boron and carbon samples for an efficiency measurement.
The effects of a possible contamination in the efficiency samples tagged with S11 and the calorimeter (S11+CALO tag) have been investigated by considering a single TOF layer and measuring the charge selection efficiency both on the event set tagged with S11+CALO and on a purer sample obtained by adding the other TOF planes to the S11+CALO tag. The two efficiencies were found to be consistent within statistical errors for each layer. No effect due to contamination in the S11+CALO tagged set was observed.
3.4. Corrections
The selected boron and carbon samples are contaminated by secondaries produced during fragmentation processes occurring in the aluminum dome on top of the pressurized vessel hosting PAMELA. This effect has been studied with a Monte Carlo calculation based on the FLUKA code (Battistoni et al. 2007) by simulating the cosmic spectra for C and O, which are the main contributors to the contamination. The resulting contamination is of the order of 10−3 for carbon, whereas for boron it ranges from about 5% at some GV up to about 20% at ∼200 GV, coming mainly from spallation of carbon.
After subtracting the contamination, the rigidity distributions of boron and carbon events have been corrected for folding effects using a Bayesian procedure (D'Agostini 1995), in order to obtain the distributions at the top of payload. These effects include possible rigidity displacements at high energies due to the finite position resolution of the silicon tracking layers and the energy loss of low-energy nuclei traversing the apparatus. The smearing matrix was derived using the GEANT4 simulations.
Interactions with the aluminum dome also remove primaries from the selected samples. Elastic scattering processes can remove primaries from the instrument acceptance or slow them down so that they are swept out by the magnetic field; inelastic scattering can destroy the primary. A correction factor for these effects has been evaluated using the FLUKA simulations, and applied to the unfolded event count. The correction is almost flat above 10 GV and amounts to 15% for carbon and 14% for boron, increasing at lower energies because of energy loss. These numbers have been treated as corrections to the geometrical factor for the two nuclear species. The resulting geometrical factors are shown in Figure 4.
Figure 4. Effective geometrical factors including the fiducial containment criterion and the correction for interactions of primary particles above the tracker. The dashed lines are fits used to obtain asymptotic values at high energy.
Download figure:
Standard image High-resolution imageEnergy loss in the apparatus may lower the measured rigidity below the critical rigidity, leading to rejection of galactic nuclei with initial rigidity above the critical one. A "cutoff correction factor" for each nuclear species was computed by assigning a random cutoff value (distributed as observed for in-flight values) to events simulated with GEANT4 and deriving the fraction of rejected events. This correction factor rises from about 0.97 at 2 GV to unity (i.e., no correction) at 3 GV and above.
3.5. Live Time
The live time of the apparatus is measured by on-board clocks and has been evaluated as a function of the vertical cutoff as the time spent in regions where the critical rigidity is below the lower limit of the rigidity bin. The total live time is constant at a value of ∼3.14 × 107 s for rigidities above 20 GV and decreases at lower rigidities because of the shorter time spent by the satellite in high latitude (i.e., low cutoff) regions down to ∼1.00 × 107 s at 2 GV. The overall error on live time determination is less than 0.2%, and has therefore been neglected.
3.6. Geometrical Factor
Due to the requirement of track containment inside a fiducial volume (see Section 3.2), the effective geometrical factor turns out to be lower than the nominal one, and assumes a constant value of 19.9 cm2 sr above 1 GV. This value has been cross-checked using two different numerical methods. The first one is a numerical computation of the integral defining the geometrical factor (Sullivan 1971), taking into account the curvature of the track due to the magnetic field, while the second method relies on a Monte Carlo simulation (Sullivan 1971). The two methods yield results differing by less than 0.1%. This error has also been neglected.
3.7. Flux Computation
The fluxes have been computed both as functions of rigidity and as functions of kinetic energy per nucleon. For each bin i, the event count corrected for the effects described in Section 3.4 was divided by the live time ΔTi, the effective geometrical factor
, the total selection efficiency
i and the bin width Δρi or ΔEi. The flux expressed as a function of rigidity is computed as
![](https://content.cld.iop.org/journals/0004-637X/791/2/93/revision1/apj498735ufd1.gif)
while as a function of kinetic energy per nucleon:
![](https://content.cld.iop.org/journals/0004-637X/791/2/93/revision1/apj498735ufd2.gif)
For boron, the latter formula needs to properly account for isotopic composition, as explained in the next section.
3.8. Isotopic Composition for Boron
In cosmic rays, both the isotopes 10B and 11B are present in comparable quantities. Since the event selection did not distinguish between them and since the events are binned according to their rigidity, a given value for the isotopic composition of boron must be assumed in order to perform the measurements as functions of kinetic energy per nucleon. Large uncertainties plague the available estimates of the isotopic composition of boron. Direct measurements are available only at relatively low energies (Ahlen et al. 2000; Hams et al. 2004; Aguilar et al. 2011). Galactic propagation models predict a high-energy value for the 10B fraction (i.e., 10B/(10B+11B)) which is weakly dependent on kinetic energy per nucleon and whose consensus value is . This value has been used in this analysis for the whole energy range.
The boron flux has been evaluated considering two different hypotheses: pure 10B and pure 11B. Assuming a binning in kinetic energy per nucleon, the corresponding binning in rigidity for each of the two hypotheses has been derived. Event selection, efficiency measurements, flux computation and corrections have then been performed in the same way for the two binnings. The two boron fluxes are combined to obtain the final flux, considering that each bin of each flux distribution contains 10B and 11B events with the same rigidity but different energy due to the different masses. Consequently, in each bin the isotopic fraction does not resemble the usual fraction expressed as a function of kinetic energy per nucleon, and a simple bin-by-bin linear combination of the two fluxes using as the weight would lead to an incorrect result. A fraction FB(ρ) has been derived by means of Monte Carlo simulations and used as a weight in order to linearly combine the two boron fluxes bin by bin and obtain a final flux. A detailed description of the calculation is presented in Appendix B.
4. RESULTS
The observed number of selected boron and carbon events, the absolute fluxes and the B/C flux ratio are reported in Tables 1 and 2. The quoted systematic uncertainties are discussed in detail in Appendix A. The fluxes and the B/C ratio are also shown in Figures 5 and 6 along with measurements from other experiments and a theoretical calculation based on GALPROP. The details of the calculation are described in Section 5. The mean kinetic energy 〈E〉 and the mean rigidity 〈ρ〉 for each bin have been computed according to (Lafferty & Wyatt 1995) using an iterative procedure starting from the middle point of each bin. The resulting mean energies and rigidities for boron and carbon differ by less than 1% and have been considered to be equal.
Figure 5. Absolute boron and carbon fluxes multiplied by E2.7 (upper panel) and B/C flux ratio (lower panel) as measured by PAMELA, together with results from other experiments (AMS02 (Oliva et al. 2013), CREAM (Ahn et al. 2008), TRACER (Obermeier et al. 2011), ATIC-2 (Panov et al. 2008), HEAO (Engelmann et al. 1990), AMS01 (Aguilar et al. 2010), CRN (Swordy et al. 1990)) and a theoretical calculation based on GALPROP (see Section 5), as functions of kinetic energy per nucleon. For PAMELA data the error bars represent the statistical error and the shaded area is the overall systematic uncertainty summarized in Appendix A.
Download figure:
Standard image High-resolution imageFigure 6. Absolute boron and carbon fluxes multiplied by ρ2.7 (upper panel) and B/C flux ratio (lower panel) as measured by PAMELA, together with a theoretical calculation based on GALPROP (see Section 5), as functions of rigidity. The error bars represent the statistical error and the shaded area is the overall systematic uncertainty summarized in Appendix A, except for the boron mixing error which does not affect the rigidity-dependent boron flux.
Download figure:
Standard image High-resolution imageTable 1. Observed Number of Events, Absolute Fluxes, and the B/C Flux Ratio as Function of Kinetic Energy per Nucleon
Kinetic Energy | 〈E〉 | C Events | C Flux | 10B Events | 11B Events | B Flux | B/C |
---|---|---|---|---|---|---|---|
at Top of Payload | Value ± Stat. ± Syst. | Value ± Stat. ± Syst. | Value ± Stat. ± Syst. | ||||
(GeV/n) | (GeV/n) | (GeV/n m2 s sr)−1 | (GeV/n m2 s sr)−1 | ||||
0.44–0.58 | 0.49 | 5146 | (5.26 ± 0.08 ± 0.26) | 1566 | 1795 | ![]() |
![]() |
0.58–0.76 | 0.65 | 6651 | (4.27 ± 0.05 ± 0.21) | 1955 | 2092 | ![]() |
![]() |
0.76–1.00 | 0.85 | 7359 | (3.30 ± 0.04 ± 0.16) | 2300 | 2320 | ![]() |
![]() |
1.00–1.30 | 1.13 | 7578 | (2.45 ± 0.03 ± 0.12) | 2351 | 2248 | ![]() |
![]() |
1.30–1.71 | 1.50 | 7033 | (1.612 ± 0.019 ± 0.078) | 2281 | 2166 | ![]() |
![]() |
1.71–2.24 | 1.94 | 6369 | (1.057 ± 0.013 ± 0.051) | 1960 | 1737 | ![]() |
![]() |
2.24–2.93 | 2.53 | 5673 | (6.70 ± 0.09 ± 0.32) × 10−1 | 1691 | 1553 | ![]() |
![]() |
2.93–3.84 | 3.34 | 4795 | (3.99 ± 0.06 ± 0.20) × 10−1 | 1350 | 1202 | (1.03 ± 0.03 ± 0.07) × 10−1 | ![]() |
3.84–5.03 | 4.36 | 3990 | (2.32 ± 0.04 ± 0.12) × 10−1 | 1078 | 945 | (5.8 ± 0.2 ± 0.4) × 10−2 | ![]() |
5.03–6.60 | 5.73 | 3270 | (1.31 ± 0.02 ± 0.07) × 10−1 | 811 | 704 | ![]() |
![]() |
6.60–8.65 | 7.49 | 2717 | (7.32 ± 0.14 ± 0.38) × 10−2 | 612 | 540 | (1.56 ± 0.07 ± 0.12) × 10−2 | ![]() |
8.65–11.3 | 9.81 | 2048 | (3.65 ± 0.08 ± 0.19) × 10−2 | 454 | 369 | (7.8 ± 0.4 ± 0.6) × 10−3 | (2.128 ± 0.12 ± 0.20) × 10−1 |
11.3–14.9 | 12.9 | 1337 | (1.81 ± 0.05 ± 0.10) × 10−2 | 253 | 217 | (3.6 ± 0.2 ± 0.3) × 10−3 | ![]() |
14.9–19.5 | 16.9 | 851 | (9.0 ± 0.3 ± 0.5) × 10−3 | 149 | 121 | ![]() |
(1.73 ± 0.15 ± 0.17) × 10−1 |
19.5–25.5 | 22.1 | 571 | (4.6 ± 0.2 ± 0.3) × 10−3 | 85 | 69 | (6.7 ± 0.7 ± 0.6) × 10−4 | (1.45 ± 0.16 ± 0.14) × 10−1 |
25.5–43.8 | 32.6 | 590 | (1.67 ± 0.07 ± 0.07) × 10−3 | 79 | 65 | (2.1 ± 0.2 ± 0.1) × 10−4 | (1.22 ± 0.13 ± 0.09) × 10−1 |
43.8–75.3 | 55.7 | 225 | (3.8 ± 0.3 ± 0.2) × 10−4 | 31 | 24 | (4.2 ± 0.6 ± 0.3) × 10−5 | (1.11 ± 0.18 ± 0.08) × 10−1 |
75.3–129 | 95.6 | 86 | (8.5 ± 0.8 ± 0.4) × 10−5 | 9 | 7 | (8.4 ± 1.5 ± 0.5) × 10−6 | (10 ± 2 ± 0.7) × 10−2 |
Note. Both the event counts for pure 10B and pure 11B hypotheses are reported.
Download table as: ASCIITypeset image
Table 2. Observed Number of Events, Absolute Fluxes, and the B/C Flux Ratio as Function of Rigidity
Rigidity | 〈ρ〉 | C Events | C Flux | B Events | B Flux | B/C |
---|---|---|---|---|---|---|
at Top of Payload | Value ± Stat. ± Syst. | Value ± Stat. ± Syst. | Value ± Stat. ± Syst. | |||
(GV) | (GV) | (GV m2 s sr)−1 | (GV m2 s sr)−1 | |||
2.02–2.38 | 2.19 | 5146 | (2.01 ± 0.03 ± 0.10) | 1566 | (6.26 ± 0.16 ± 0.29) × 10−1 | (3.12 ± 0.09 ± 0.21) × 10−1 |
2.38–2.82 | 2.57 | 6651 | (1.73 ± 0.02 ± 0.08) | 1955 | (5.49 ± 0.13 ± 0.25) × 10−1 | (3.17 ± 0.08 ± 0.21) × 10−1 |
2.82–3.37 | 3.06 | 7359 | (1.413 ± 0.017 ± 0.068) | 2300 | (4.72 ± 0.10 ± 0.21) × 10−1 | (3.34 ± 0.08 ± 0.22) × 10−1 |
3.37–4.06 | 3.67 | 7578 | (1.093 ± 0.013 ± 0.053) | 2351 | (3.67 ± 0.08 ± 0.17) × 10−1 | (3.35 ± 0.08 ± 0.22) × 10−1 |
4.06–4.93 | 4.45 | 7033 | (7.44 ± 0.09 ± 0.36) × 10−1 | 2281 | (2.52 ± 0.05 ± 0.11) × 10−1 | (3.39 ± 0.08 ± 0.22) × 10−1 |
4.93–6.06 | 5.44 | 6369 | (5.00 ± 0.06 ± 0.24) × 10−1 | 1960 | (1.56 ± 0.04 ± 0.07) × 10−1 | (3.12 ± 0.08 ± 0.21) × 10−1 |
6.06–7.50 | 6.70 | 5673 | (3.23 ± 0.04 ± 0.16) × 10−1 | 1691 | (9.9 ± 0.2 ± 0.5) × 10−2 | (3.06 ± 0.08 ± 0.21) × 10−1 |
7.50–9.36 | 8.34 | 4795 | (1.95 ± 0.03 ± 0.10) × 10−1 | 1350 | (5.52 ± 0.15 ± 0.34) × 10−2 | (2.83 ± 0.09 ± 0.23) × 10−1 |
9.36–11.8 | 10.4 | 3990 | (1.143 ± 0.018 ± 0.058) × 10−1 | 1078 | (3.19 ± 0.10 ± 0.20) × 10−2 | (2.79 ± 0.10 ± 0.23) × 10−1 |
11.8–15.0 | 13.2 | 3270 | (6.49 ± 0.11 ± 0.33) × 10−2 | 811 | (1.70 ± 0.06 ± 0.11) × 10−2 | (2.61 ± 0.10 ± 0.22) × 10−1 |
15.0–19.1 | 16.8 | 2717 | (3.64 ± 0.07 ± 0.19) × 10−2 | 612 | (8.8 ± 0.4 ± 0.6) × 10−3 | (2.43 ± 0.11 ± 0.21) × 10−1 |
19.1–24.5 | 21.4 | 2048 | (1.82 ± 0.04 ± 0.10) × 10−2 | 454 | (4.4 ± 0.2 ± 0.3) × 10−3 | (2.42 ± 0.13 ± 0.21) × 10−1 |
24.5–31.5 | 27.6 | 1337 | (9.1 ± 0.3 ± 0.5) × 10−3 | 253 | (2.02 ± 0.13 ± 0.15) × 10−3 | (2.24 ± 0.15 ± 0.20) × 10−1 |
31.5–40.8 | 35.6 | 851 | (4.51 ± 0.16 ± 0.24) × 10−3 | 149 | (8.9 ± 0.7 ± 0.7) × 10−4 | (1.96 ± 0.18 ± 0.18) × 10−1 |
40.8–52.9 | 46.1 | 571 | (2.32 ± 0.10 ± 0.13) × 10−3 | 85 | (3.9 ± 0.4 ± 0.3) × 10−4 | (1.7 ± 0.2 ± 0.16) × 10−1 |
52.9–89.5 | 67.1 | 590 | (8.4 ± 0.4 ± 0.3) × 10−4 | 79 | (1.18 ± 0.14 ± 0.07) × 10−4 | (1.41 ± 0.18 ± 0.10) × 10−1 |
89.5–152 | 113 | 225 | (1.92 ± 0.14 ± 0.08) × 10−4 | 31 | (2.5 ± 0.5 ± 0.1) × 10−5 | (1.3 ± 0.3 ± 0.09) × 10−1 |
152–260 | 193 | 86 | (4.3 ± 0.4 ± 0.2) × 10−5 | 9 | (4.7 ± 1.1 ± 0.3) × 10−6 | (1.1 ± 0.3 ± 0.08) × 10−1 |
Download table as: ASCIITypeset image
The discrepancies with other experiments at low energies can be reasonably ascribed to solar modulation effects. The data used for this analysis were taken by PAMELA during an unusually quiet solar minimum period, resulting in an enhanced flux of galactic cosmic rays at low energies in the heliosphere, which has already been observed for protons (Adriani et al. 2013) and nuclei (Mewaldt et al. 2010). Above 6 GeV/n the fluxes are in overall agreement with the other available measurements, especially with those from HEAO and CREAM. A power-law fit above 20 GeV/n results in a spectral index γB = 3.01 ± 0.13 for boron and γC = 2.72 ± 0.06 for carbon.
5. DISCUSSION
A comprehensive and detailed study of the results presented above is beyond the scope of this paper. The following discussion is intentionally limited to a single propagation model in order to compute an estimate of the most significant propagation parameters from the PAMELA boron and carbon data. Results may vary when considering different models or propagation software packages.
The data presented in the previous section as a function of kinetic energy per nucleon has been fitted with a diffusive cosmic ray propagation model using the GALPROP code interfaced with the MIGRAD minimizer in the MINUIT2 minimization package distributed within the ROOT framework (Brun & Rademakers 1997). Only a few parameters have been left free because of the high computation time required for multiple GALPROP runs. The values for the other parameters have been taken from (Vladimirov 2012). The diffusion coefficient is found to have a fitted slope value of δ = 0.397 ± 0.007 and a normalization factor D0 = (4.12 ± 0.04) × 1028 cm2 s−1. Other fitted parameters are the solar modulation parameter in the force-field approximation Φ = (0.40 ± 0.01) GV and the overall normalization of the fluxes N = 1.04 ± 0.03. The result of the fit is shown in Figures 5 and 6. A contour plot of the confidence intervals for δ and D0 is shown in Figure 7.
Figure 7. Contour plot of the 1σ, 2σ and 3σ confidence levels for δ and D0.
Download figure:
Standard image High-resolution imageThe fitted value for δ falls between the predicted values for Kolmogorov (δ = 1/3) and Kraichnan (δ = 1/2) diffusion types, thus the PAMELA data cannot distinguish between these two types.
We acknowledge support from The Italian Space Agency (ASI), Deutsches Zentrum für Luft-und Raumfahrt (DLR), The Swedish National Space Board, The Swedish Research Council, The Russian Space Agency (Roscosmos), and The Russian Science Foundation.
APPENDIX A: SYSTEMATIC UNCERTAINTIES
The following contributions to the systematic uncertainty have been considered.
- 1.Selection efficiencies. The measurement of the tracking and charge selection efficiencies from flight data is performed using samples of finite size. The associated statistical error has been propagated to the flux as a systematic uncertainty.
- 2.Fiducial containment. The finite tracking resolution of the calorimeter can lead to a contamination of the tracking efficiency sample by events coming from outside the fiducial acceptance, and possibly also crossing the magnet walls. These can in principle be eliminated by further restricting the fiducial volume for both event selection and efficiency measurement, but this would significantly reduce the sample sizes. The chosen approach is to use protons from both flight and simulated data to measure the tracking efficiency for both the fiducial volume defined in Section 3.2 and a more restrictive one. Their relative difference is taken as an estimate of the systematic uncertainty, which is about 2%. Monte Carlo simulations give results for boron and carbon which are consistent with the one obtained with protons. The uncertainty is propagated to the final flux.
- 3.Monte Carlo correction factor for the tracking efficiency. This correction factor should introduce only relatively small errors, since it is computed as the ratio of two Monte Carlo efficiencies. Systematic effects should largely cancel out. The correction factor is constant at 0.97 for both boron and carbon. That this factor remains constant at high rigidity is due to the isotropic efficiency being constant at relativistic rigidities. A conservative factor of 3% has been taken as an estimate of the systematic uncertainty on the flux because of this correction factor.
- 4.Residual coherent misalignment of the spectrometer. The spectrometer alignment procedure results in a residual coherent misalignment producing a systematic shift in the measured rigidity. The error estimation procedure is described in the Supporting Online Material of Adriani et al. (2011). This error has been propagated to the measured flux. It is negligible at low energy and increases up to about 2% at 250 GV.
- 5.Cutoff, contamination and geometrical factor corrections. All these factors have been evaluated on finite-size samples, so they are affected by a statistical error which has been propagated to the flux as a systematic uncertainty.
- 6.Unfolding. The unfolding error has been assessed by means of the procedure described in the Supporting Online Material of (Adriani et al. 2011), comparing a given initial spectrum and an unfolded Monte Carlo simulation. The two were found to be in agreement within 3%, so this value has been taken as the unfolding contribution to the flux error.
- 7.Isotopic composition of boron. The uncertainty associated with this poorly known parameter has been propagated to the flux by assuming the extreme values of 0.2 and 0.5 for the 10B fraction and taking the difference between these fluxes and the one obtained with
as the estimated upper and lower errors on the flux. This error affects only the measurement expressed as a function of kinetic energy per nucleon.
The overall uncertainty has been estimated as the quadratic sum of the above terms in the hypothesis of uncorrelated errors. A summary plot is shown in Figure 8.
Figure 8. Systematic uncertainties for absolute fluxes. The total contribution is computed as the quadratic sum of the individual terms. The track selections term is the quadratic sum of the contributions from statistics and from fiducial containment. The contributions of the track efficiency correction and of the unfolding have been slightly shifted apart from their 3% value for the sake of readability.
Download figure:
Standard image High-resolution imageAPPENDIX B: ISOTOPIC COMPOSITION OF BORON
In this analysis the events have been binned according to their rigidity as measured by the magnetic spectrometer. Given that the event selection does not distinguish between the two isotopes 10B and 11B, each bin is populated by 10B and 11B events with approximately the same rigidity (within the bin limits) but different kinetic energy per nucleon because of the different mass numbers. Consequently, the isotopic composition in a given bin is not described by the 10B fraction FB expressed as a function of kinetic energy per nucleon E:
![Equation (B1)](https://content.cld.iop.org/journals/0004-637X/791/2/93/revision1/apj498735fd1.gif)
where and
are the fluxes of 10B and 11B respectively. A fraction expressed as a function of rigidity must then be derived in order to correctly account for the isotopic composition in each bin:
![Equation (B2)](https://content.cld.iop.org/journals/0004-637X/791/2/93/revision1/apj498735fd2.gif)
Using rigidity bins of finite size leads to
![Equation (B3)](https://content.cld.iop.org/journals/0004-637X/791/2/93/revision1/apj498735fd3.gif)
where FB(ρi) is the 10B fraction for the ith rigidity bin centered at ρi, while and
are the 10B and 11B event count for the same bin, respectively.
can be rewritten using the 10B fraction in kinetic energy:
![Equation (B4)](https://content.cld.iop.org/journals/0004-637X/791/2/93/revision1/apj498735fd4.gif)
Here denotes the 11B event count in a bin in kinetic energy per nucleon whose limits are obtained by converting the limits in rigidity of the ith bin to kinetic energy assuming the mass and the charge of 11B.
is the kinetic energy per nucleon of a 11B nucleus of rigidity ρi. Then, by construction, the first equality in the above equation follows. The second equality follows from the definition of FB(Ei) which is the equivalent of Equation (B3) for kinetic energy bins. Note that
![Equation (B5)](https://content.cld.iop.org/journals/0004-637X/791/2/93/revision1/apj498735fd5.gif)
since the limits of the energy and rigidity bins do not correspond for 10B. Converting the bin limits in energy back to rigidity but assuming now the mass and the charge of 10B yields
![Equation (B6)](https://content.cld.iop.org/journals/0004-637X/791/2/93/revision1/apj498735fd6.gif)
is then the rigidity of a 10B nucleus having the same kinetic energy per nucleon
of a 11B nucleus of rigidity ρi (the same relation holds between the limits of the bins centered in
and ρi). To obtain the explicit relationship between
and ρ, write
as
![Equation (B7)](https://content.cld.iop.org/journals/0004-637X/791/2/93/revision1/apj498735fd7.gif)
where Z is the atomic number of boron, A10 is the mass number of 10B and mp is the proton mass, and then as a function of ρi:
![Equation (B8)](https://content.cld.iop.org/journals/0004-637X/791/2/93/revision1/apj498735fd8.gif)
with A11 the mass number of 11B. It follows that
![Equation (B9)](https://content.cld.iop.org/journals/0004-637X/791/2/93/revision1/apj498735fd9.gif)
The final form of the rigidity-dependent 10B fraction is then
![Equation (B10)](https://content.cld.iop.org/journals/0004-637X/791/2/93/revision1/apj498735fd10.gif)
where the approximated energy-independent value has been used.
Generally speaking, the fraction expressed as a function of rigidity is not constant and depends on the spectral shape. To account for this a toy Monte Carlo simulation of realistic 10B and 11B spectra taken from a galactic propagation model has been set up, the resulting event counts have been trimmed to reproduce and finally the events have been binned according to their rigidity for both the pure 10B and pure 11B hypotheses. Knowing the fraction in each rigidity bin of the two binnings one can express the final boron flux in the ith energy bin as
![Equation (B11)](https://content.cld.iop.org/journals/0004-637X/791/2/93/revision1/apj498735fd11.gif)
where and
are the 10B fraction obtained from the toy Monte Carlo in the ith rigidity bin for pure 10B and pure 11B hypotheses respectively, and ϕ10(Ei) and ϕ11(Ei) are the experimental fluxes for the pure 10B and pure 11B hypotheses respectively (see Section 3.8).
To assess the difference between the 10B fraction as a function of kinetic energy per nucleon and as a function of rigidity, Equation (B10) can be computed at high energies. Above few GeV/n, where the spectrum can be well described with a power-law function with index γ, Equation (B10) gives a 10B fraction
![Equation (B12)](https://content.cld.iop.org/journals/0004-637X/791/2/93/revision1/apj498735fd12.gif)
which is in agreement with the value obtained from the toy Monte Carlo and differs from by about 18%.
Footnotes
- 24
The MDR is defined as the rigidity with an associated 100% error due to the finite spatial resolution of the spectrometer.
- 25
1 MIP is defined here as the most probable energy release of a |Z| = 1 minimum ionizing particle.
- 26
S12 is the lowest of the two layers constituting S1; the upper layer S11 was used for efficiency measurement as explained in Section 3.3.