Articles

TIME DEPENDENCE OF THE PROTON FLUX MEASURED BY PAMELA DURING THE 2006 JULY–2009 DECEMBER SOLAR MINIMUM

, , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , and

Published 2013 February 20 © 2013. The American Astronomical Society. All rights reserved.
, , Citation O. Adriani et al 2013 ApJ 765 91 DOI 10.1088/0004-637X/765/2/91

0004-637X/765/2/91

ABSTRACT

The energy spectra of galactic cosmic rays carry fundamental information regarding their origin and propagation. These spectra, when measured near Earth, are significantly affected by the solar magnetic field. A comprehensive description of the cosmic radiation must therefore include the transport and modulation of cosmic rays inside the heliosphere. During the end of the last decade, the Sun underwent a peculiarly long quiet phase well suited to study modulation processes. In this paper we present proton spectra measured from 2006 July to 2009 December by PAMELA. The large collected statistics of protons allowed the time variation to be followed on a nearly monthly basis down to 400 MV. Data are compared with a state-of-the-art three-dimensional model of solar modulation.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

Protons are the most abundant species in cosmic rays, representing about 90% of the total flux. They are believed to be mainly of primordial origin, accelerated during supernovae explosions and confined within our Galaxy, at least up to PeV energies. During propagation through the interstellar medium, cosmic rays undergo interactions with gas atoms and lose energy significantly modifying their injection spectra and composition.

Besides the effect of propagation in the Galaxy, cosmic-ray particles reaching the Earth are affected by the solar wind and the solar magnetic field. A continuous flow of plasma emanates from the solar corona and extends out far beyond Pluto, up to the heliopause where the interstellar medium is encountered. The solar wind travels at supersonic speeds of about 350–650 km s−1 up to the termination shock where it slows down to about 100–150 km s−1 to form the inner heliosheath.

The plasma carries the solar magnetic field out into the heliosphere to form the heliospheric magnetic field (HMF). When interstellar cosmic rays enter the heliosphere, they interact with the solar wind and the HMF that modify their energy spectra. These modifications are directly related to solar activity and the total effect is called solar modulation. Solar activity varies strongly with time, rising from a minimum level when the Sun is quiet (with cosmic rays then having a maximum intensity at Earth) to a maximum period (when cosmic rays have a minimum intensity) and then returning to a new minimum to repeat the cycle. This solar cycle spans approximately 11 years. Solar activity modulates cosmic rays with the same cycle as a function of position in the heliosphere. Protons with rigidities up to at least 30 GV are affected and the effect becomes progressively larger as the rigidity decreases (e.g., Potgieter 2011). At each solar activity maximum the polarity of the solar magnetic field reverses. When the solar magnetic field lines point outward in the northern hemisphere and inward in the southern hemisphere, the Sun is said to be in an A > 0 cycle, whereas during an A < 0 cycle the solar magnetic field lines point inward (outward) in the northern (southern) hemisphere. Direct cosmic-ray measurements above a few MeV have been performed from the inner to the outer heliosphere (see, e.g., Heber & Potgieter 2006; Webber et al. 2012) and particularly at Earth where solar modulation effects are at their largest.

Understanding the origin and propagation of cosmic rays requires knowledge of the cosmic-ray energy spectra in the interstellar medium, i.e., uninfluenced by the Sun's magnetic field. These spectra, especially for protons, can be inferred from gamma-ray data (e.g., see Neronov et al. 2012). However, uncertainties are large and direct cosmic-ray measurements are preferable. Voyager space missions (Decker et al. 2005; Richardson et al. 2008) have passed the termination shock and may soon start sampling the cosmic radiation in the local interstellar medium. They will provide fundamental data to understand the basic features of the heliosphere. However, these cosmic-ray measurements are limited to the very low energy part of the spectra (from tens to a few 100 MeV), with galactic protons, helium, and oxygen strongly contaminated by the anomalous component (Stone et al. 2008). Higher energies are needed to fully address fundamental questions in cosmic-ray physics. Moreover, the energy spectra of rare components like antiparticles may yield indirect information regarding new physics such as dark matter. For example, the antiproton flux between a few hundred MeV and a few GeV can be used (e.g., Cerdeno et al. 2012; Asano et al. 2012; Garny et al. 2012) to constrain models with light (∼10 GeV) dark matter candidates put forward to interpret recent measurements by direct-detection experiments (Bernabei et al. 2000; Aalseth et al. 2011; Angloher et al. 2012). However, solar modulation has a relevant effect on the cosmic-ray spectra at these energies that needs to be disentangled to allow a comprehensive picture to emerge.

Precise measurements of cosmic-ray spectra over a wide rigidity range, from a few hundred MV to tens of GV over an extended period of time (covering a significant fraction of the solar cycle) can be used to study the effect of solar modulation in greater detail, including the convective and adiabatic cooling effect of the expanding solar wind and the diffusive and particle drift effects of the turbulent HMF. The most recent solar minimum activity period (2006–2009) and the consequent minimum modulation conditions for cosmic rays were unusual. It was expected that the new activity cycle would begin early in 2008. Instead, minimum modulation conditions continued until the end of 2009 when the highest levels of cosmic rays were recorded at Earth since the beginning of the space age (see, e.g., Mewaldt et al. 2010). This period of prolonged solar minimum activity is well suited to study the various modulation processes that affect the propagation of galactic cosmic rays inside the heliosphere.

PAMELA is a satellite-borne instrument designed for cosmic-ray antimatter studies. The instrument is flying on board the Russian Resurs-DK1 satellite since 2006 June, in a semi-polar near-Earth orbit. PAMELA has provided important results on the antiproton (Adriani et al. 2009a) and positron (Adriani et al. 2009b) galactic abundances. The high-resolution magnetic spectrometer allowed hydrogen and helium spectral measurements up to 1.2 TV (Adriani et al. 2011), the highest limit ever achieved by this kind of experiment. We present here results on the long-term variations in the absolute flux of hydrogen nuclei, measured down to 400 MV.

The analysis, described in Section 3, is based on data collected between 2006 July and 2009 December. The analyzed period covers the unusually long most recent solar minimum. The large proton statistics collected by the instrument allowed the proton flux to be measured for each Carrington rotation23 (hereafter the flux measured on this time basis will be referred as Carrington flux).

2. THE INSTRUMENT

The PAMELA instrument, shown schematically in Figure 1, comprises a number of highly redundant detectors, capable of identifying particles by providing charge (Ze), mass, rigidity (R = pc/Ze, with p momentum of the particle and c speed of light), and velocity (β) measurements over a very wide energy range. Multiple detectors are placed around a permanent magnet with a silicon micro-strip tracking system, providing charge and track deflection (d = ±1/R) information. A scintillator system (TOF: S1, S2, and S3) provides the trigger for event acquisition, the time-of-flight measurement, and further charge information. A silicon–tungsten calorimeter is used to perform hadron/lepton separation; a scintillator (S4) and a neutron detector placed at the bottom of the apparatus increase the rejection power at high energy. An anti-coincidence system (AC: CAS, CAT, and CARD) registers hits produced by particles entering from outside the instrument acceptance. A detailed description of the apparatus can be found in Picozza et al. (2007).

Figure 1.

Figure 1. Schematic view of the PAMELA apparatus.

Standard image High-resolution image

3. ANALYSIS

3.1. Event Selection

In order to give a valid trigger for data acquisition, a particle must cross the S2 and S3 scintillators, located above and below the tracking system, respectively. The requirement of a hit on S3 sets a limit on the minimum detectable rigidity, which for vertical protons is on average 400 MV, i.e., ≃ 80 MeV of kinetic energy.

Clean events are selected by requiring a precision track in the spectrometer, a trajectory contained within the instrument acceptance, and no hits in the CARD and CAT scintillators. Protons are identified among all the selected particles by requiring a positive track curvature and ionization energy losses (dE/dX) in the tracking system consistent with a hydrogen nucleus. The dE/dX selection included deuterons, which in this work have been treated as being protons. Galactic particles have been selected by requiring that the measured rigidity was a factor 1.3 above the geomagnetic vertical cutoff, evaluated using the Störmer approximation, at a given orbital position. The satellite orbit has an inclination of 70°, which allows galactic particles to be studied down to the minimum detectable rigidity.

The analysis procedure was similar to previous work on high-energy proton and helium fluxes (Adriani et al. 2011). The major differences concern track-quality requirements and the fiducial volume definition. The procedure used in the proton and helium study aimed to maximize the spectrometer performance and minimize systematic uncertainties in the flux normalization, in order to perform a precise spectral measurement up to the highest achievable energy. The study presented in this paper aims to trace time variations in the low-energy particle flux, with the finest possible time binning. Hence, high-performing tracking was less important than high statistics. For this reason the tracking requirements were relaxed (at least three hits, instead of four, on both X and Y view, a track lever arm of at least four silicon planes in the tracker) and the fiducial volume was enlarged (bounded 1.5 mm from the magnet cavity walls, instead of 7 mm), allowing the collection of more particles with only a small degradation in the measurement precision in the energy range under study.

3.2. Flux Normalization

The absolute proton flux in the spectrometer was obtained by dividing the measured energy spectrum by the acquisition time, the geometrical acceptance, and the selection efficiency.

3.2.1. Acquisition Time

Flux evaluation was performed as a function of time, binning the data according to the Carrington rotations. The time period studied corresponds to Carrington rotation numbers 2045–2092. The acquisition time was evaluated for each energy bin as the total live time spent above the geomagnetic cutoff, outside the South Atlantic Anomaly and during periods when no solar events had occurred (hence 2006 December, when a large solar event took place, was excluded from the data). This work is based on data collected between 2006 July and 2009 December. The total live time was about 5 × 107 s (∼590 days) above ∼20 GV, reducing to about 13% of this value at 0.4 GV, due to the relatively short time spent by the satellite at low geomagnetic locations. The accuracy of the live-time determination was cross-checked by comparing different clocks available in flight, which showed a relative difference of less than 0.2%.

3.2.2. Geometrical Factor

The geometrical factor for this analysis (G = 19.9 cm2 sr) is determined by the requirement of containment within the fiducial volume: a reconstructed track had to be 1.5 mm away from the magnet walls and the TOF-scintillator edges. The geometrical factor is constant above 1 GV and slowly decreases by ∼2% to lower energies, where the curvature of the proton particle trajectory is no longer negligible.

The value of G has been evaluated with two independent numerical methods, which provided values consistent within 0.1%.

3.2.3. Selection Efficiency

The PAMELA instrument has been designed to provide redundant information on event topology and particle characteristics. This feature allowed the efficiencies of each detector to be measured directly from flight data. Due to the large proton statistics, detector behavior could be monitored over very short timescales. For this analysis, efficiencies were evaluated independently for each time interval, in order to follow any possible time variation of detector responses.

The major time-dependent effect was the sudden failure of some front-end chips in the tracking system. Chip failures started to occur few months after the launch of the satellite, initially with an almost constant rate. In the following years the rate of failures significantly decreased and the tracking system reached an almost stable condition near the end of 2009. The net effect was a progressive reduction of the tracking efficiency, since the number of hits available for track reconstruction decreased. No degradation in the signal-to-noise ratio and spatial resolution was observed.

The approach described in Adriani et al. (2011) for the evaluation of the tracking efficiency was used: a high energy average value was experimentally obtained using non-interacting minimum ionizing particles in the calorimeter from flight data and an energy-dependent correction factor was evaluated using Monte Carlo simulation. The average tracking efficiency was measured using a test sample of non-interacting relativistic protons in the imaging calorimeter. The resulting efficiency varied almost monotonically from ≃94% during the first months to ≃20% in the final months. The efficiency variation was correctly reproduced by the simulation, which included a map of dead channels as a function of time. Since the requirement of non-interacting minimum ionizing particles in the calorimeter selected essentially protons with rigidities ⩾2 GV, the determination of the tracking efficiency at lower rigidities was performed using the simulation. The simulation results were normalized to experimental values at high rigidities after performing an identical selection. This was required because, due to the relatively poor tracking resolution of the calorimeter, the average efficiency values obtained using the calorimeter selection were different from the average efficiency of an isotropically distributed particle sample (as was the case for triggered events). This normalization was applied to the simulation estimated efficiency that accounted for the full energy dependence of the tracking efficiency and was used for the flux estimation. When comparing the average experimental and simulated efficiencies, differences of less than 4% were found, except during the last few months of 2009 when the tracking efficiency fell below 40% and, consequently, the discrepancy became more significant, increasing up to 13%. These differences were included in the systematic uncertainties of the measurements. As a further check, two different simulation packages were used: GEANT3 (Brun et al. 1994) and GEANT4 (Agostinelli et al. 2003). The resulting efficiencies differed by ≃ 1.5% and this was included in the systematic uncertainty.

TOF and AC selection efficiencies were evaluated relative to the tracking selection. For well-contained events, the TOF and AC responses were not correlated and the selection efficiencies could be measured independently. At low rigidities TOF and AC efficiency measurements were affected by background contamination in the test samples of low-energy secondary particles produced by high-energy primaries interacting with the apparatus. This contamination resulted in an underestimation of the efficiency, since background interactions were often associated with multiple hits in both the TOF and AC scintillators. An upper limit on this contamination was determined by excluding the TOF and AC conditions from the test sample selection, thereby adding background events into the test sample. The difference between the results was found to be ≃6% in the lowest-energy bin. This value was conservatively treated as an asymmetric systematic uncertainty on the absolute flux.

3.3. Spectral Unfolding and Other Corrections

The proton spectrum measured in the spectrometer was corrected to the top of the payload by accounting for the effects of rigidity displacement due to finite spectrometer resolution and particle slow-down by ionization energy loss, the latter being more relevant at low energies. The correction was applied by means of an unfolding procedure, following the Bayesian approach proposed by D'Agostini (1995). The detector response matrix was obtained from the simulation. The resulting spectrum was corrected for the contamination of locally produced secondary particles and for attenuation due to particles interacting with the satellite and the instrument itself. In order to estimate the residual contamination of locally produced secondary particles, simulations of protons and helium nuclei isotropically hitting the PAMELA payload within a 2π solid angle were carried out. The estimated background was found to be less than 6% at 0.4 GV and inversely proportional to energy, becoming negligible above 1 GV. Flux attenuation, which was slightly energy dependent and amounted to 5%–6%, accounted for both particle loss by non-elastic interactions and loss and gain of particles from within the acceptance due to elastic scattering.

Finally, in order to correct for systematic time-dependent variations of the Carrington fluxes (see next section), these were normalized, at high energy (30–50 GeV), to the proton flux measured with the more stringent tracking conditions of Adriani et al. (2011) and averaged over the period 2006 July–2008 March.24

3.4. Systematic Uncertainties

By relaxing the tracking selection the proton statistics could be significantly increased, allowing the proton flux to be evaluated on a Carrington rotation basis up to the end of 2009. This approach also led to an increase of the systematic uncertainties on the measured absolute fluxes.

The reliability of the flux determination was verified by studying the time variation of the flux at high energy (above 30 GeV) since no significant modulation effects are expected in this energy region. As a reference, we considered the fluxes obtained with the stringent tracking conditions and averaged over the period 2006 July–2008 March. Figure 2 shows the time dependence of the ratio between the Carrington fluxes and the reference flux for two high-energy intervals at about 30 GV and 47 GV. Two different selection criteria were considered: stringent (open circles; as in Adriani et al. 2011) and relaxed tracking selections (filled circles; this work).

Figure 2.

Figure 2. Carrington proton fluxes shown relative to the flux obtained in Adriani et al. (2011; value averaged over the period 2006 July–2008 March) for stringent tracking criteria (open circles) and relaxed tracking criteria (filled circles). Two energy intervals are presented.

Standard image High-resolution image

Comparing the upper and lower panels of Figure 2 a correlation between the time profiles at different energies can be noticed, which indicates the presence of non-negligible systematic effects. The Carrington fluxes based on the analysis described in Adriani et al. (2011) fluctuate around a mean value consistent with the quoted ≃ 4% systematic uncertainty, while the fluxes for this analysis show a decreasing trend, with a relative variation over the whole time interval of ∼10%. This time-dependent systematic effect was corrected as discussed in the previous section. The assumption behind this correction implies that PAMELA results are not sensitive to solar modulation-induced flux variations smaller than about 5%.

Further contributions to the total uncertainty resulted from the unfolding procedure (1%–5%) and from the correction to the top of the payload (1%).

4. RESULTS

Figure 3 shows the time profile of the proton flux for three energy intervals; the error bars are statistical, while the shaded area represents the quadratic sum of all systematic uncertainties. The systematic uncertainty is dominated by the tracking efficiency and reaches a minimum above a few tens of GeV and for the first Carrington rotations, when the tracking efficiency is maximum. Additional time- and energy-dependent contributions to the overall uncertainty have been introduced in order to account for both the reduced accuracy of the simulation in reproducing the spectrometer behavior as the detection efficiency decreases in time, and the larger uncertainty on the simulated efficiency at low energy.

Figure 3.

Figure 3. Time profile of the proton flux for three energy intervals (from top to bottom): 0.597–0.625 GeV, 0.99–1.04 GeV, and 2.99–3.13 GeV. The error bars are statistical and the shaded area represents the quadratic sum of all systematic uncertainties.

Standard image High-resolution image

Figure 4 shows the evolution of the Carrington proton energy spectra from 2006 July (blue curve) to 2009 December (red curve). Around 100 MeV proton intensities increased by ≃ 20% over 12 months, from 2006 December to 2007 December, with a slightly larger increase of ≃ 25% from 2007 December to 2008 December. For 2009, when solar minimum conditions continued to prevail throughout the heliosphere, intensities increased by almost 40%. From these data it is clear that protons (evident across the greater part of the energy spectrum) reached maximum intensities at the end of 2009.

Figure 4.

Figure 4. Evolution of the proton energy spectrum as particle intensities approached the period of minimum solar activity, from 2006 July (violet) to 2009 December (red). The region between the blue and red curves indicates the spread in proton fluxes during this time.

Standard image High-resolution image

5. DATA INTERPRETATION AND DISCUSSION

Figure 5 and Table 1 present the galactic proton spectra measured by PAMELA over four time periods. These data illustrate how the proton spectra evolved from late 2006 to late 2009. The latter represents the highest proton flux observed by PAMELA. In Figure 5 the PAMELA proton spectra are overlaid with the corresponding computed spectra (solid lines). A full three-dimensional (3D) model was used to compute the differential intensity of cosmic-ray protons at Earth from 10 MeV to 30 GeV. The model is based on the numerical solution of the Parker transport equation (Parker 1965), including all four major modulation mechanisms: convection, particle drift caused by gradients and curvatures in the HMF, diffusion described by a full 3D diffusion tensor, and adiabatic energy changes. For these calculations, a local interstellar spectrum (LIS) was assumed (shown as a solid black line in Figure 5). The proton LIS was based on that by Langner & Potgieter (2004) modified at high energies to match PAMELA observations. It follows from the comparison between the model and observations that the proton spectra became progressively softer, as expected, from 2006 but requiring larger diffusion coefficients at lower energies than anticipated. At the same time full particle drifts occurred but with the modulation drift effects somewhat obscured by the effective diffusion. This caused the effect of the HMF wavy current sheet to be less well correlated to increases in cosmic-ray intensities at Earth. Details concerning the modulation model and theoretical assumptions and implications will be published in a subsequent paper (Potgieter et al. 2013).

Figure 5.

Figure 5. PAMELA proton spectra for 2006 November (blue), 2007 December (green), 2008 December (yellow), and 2009 December (red) overlaid with the corresponding computed spectra (solid lines; Potgieter et al. 2013). The LIS used for the computation is also shown (black solid line).

Standard image High-resolution image

Table 1. Hydrogen Flux Measured by PAMELA over Five Time Periods

Kinetic Energy Flux
(GeV) (m2 s sr MeV)−1
2006 Nov 13–2006 Dec 4 2007 Nov 30–2007 Dec 27 2008 Nov 19–2008 Dec 15 2009 Dec 6–2010 Jan 1
0.082–0.095 (0.721 ± 0.032+0.061−0.100) (1.045 ± 0.031+0.090−0.102 (1.191 ± 0.035+0.105−0.104) (1.730 ± 0.143+0.266−0.163)
0.095–0.10 (0.759 ± 0.036+0.063−0.096) (1.111 ± 0.030+0.094−0.098 (1.332 ± 0.034+0.116−0.100) (1.722 ± 0.156+0.263−0.160)
0.10–0.11 (0.818 ± 0.030+0.067−0.093) (1.151 ± 0.027+0.097−0.094) (1.401 ± 0.031+0.121−0.097) (2.016 ± 0.134+0.307−0.158)
0.11–0.12 (0.893 ± 0.034+0.072−0.089) (1.230 ± 0.027+0.102−0.091) (1.477 ± 0.031+0.126−0.094) (2.210 ± 0.151+0.336−0.157)
0.12–0.13 (0.928 ± 0.030+0.074−0.086) (1.334 ± 0.026+0.110−0.088) (1.543 ± 0.029+0.131−0.091) (2.162 ± 0.129+0.327−0.155)
0.13–0.15 (0.960 ± 0.028+0.076−0.084) (1.411 ± 0.025+0.115−0.086) (1.607 ± 0.027+0.135−0.089) (2.256 ± 0.123+0.341−0.154)
0.15–0.16 (1.030 ± 0.028+0.081−0.083) (1.453 ± 0.023+0.117−0.085) (1.675 ± 0.026+0.139−0.087) (2.447 ± 0.122+0.368−0.153)
0.16–0.17 (1.085 ± 0.027+0.084−0.081) (1.477 ± 0.022+0.118−0.083) (1.727 ± 0.025+0.142−0.085) (2.359 ± 0.112+0.354−0.152)
0.17–0.19 (1.123 ± 0.025+0.086−0.079) (1.545 ± 0.021+0.122−0.081) (1.805 ± 0.024+0.147−0.084) (2.277 ± 0.102+0.341−0.151)
0.19–0.21 (1.164 ± 0.024+0.088−0.078) (1.613 ± 0.020+0.126−0.080) (1.865 ± 0.023+0.150−0.082) (2.423 ± 0.097+0.362−0.150)
0.21–0.22 (1.194 ± 0.022+0.090−0.076) (1.634 ± 0.019+0.126−0.078) (1.899 ± 0.021+0.152−0.081) (2.592 ± 0.093+0.386−0.149)
0.22–0.24 (1.222 ± 0.022+0.091−0.075) (1.673 ± 0.018+0.128−0.077) (1.968 ± 0.021+0.156−0.080) (2.647 ± 0.091+0.393−0.149)
0.24–0.26 (1.278 ± 0.021+0.094−0.074) (1.732 ± 0.018+0.131−0.076) (2.031 ± 0.020+0.159−0.079) (2.682 ± 0.088+0.397−0.148)
0.26–0.29 (1.359 ± 0.020+0.098−0.073) (1.775 ± 0.017+0.133−0.075) (2.047 ± 0.019+0.159−0.078) (2.684 ± 0.082+0.396−0.148)
0.29–0.31 (1.371 ± 0.020+0.098−0.072) (1.788 ± 0.016+0.132−0.074) (2.025 ± 0.018+0.155−0.077) (2.603 ± 0.077+0.383−0.147)
0.31–0.34 (1.380 ± 0.019+0.098−0.071) (1.818 ± 0.016+0.133−0.073) (2.040 ± 0.017+0.155−0.076) (2.656 ± 0.075+0.390−0.147)
0.34–0.36 (1.381 ± 0.018+0.097−0.070) (1.822 ± 0.015+0.132−0.072) (2.036 ± 0.017+0.153−0.075) (2.633 ± 0.072+0.385−0.146)
0.36–0.39 (1.398 ± 0.017+0.097−0.069) (1.816 ± 0.014+0.130−0.072) (2.019 ± 0.016+0.150−0.074) (2.504 ± 0.067+0.366−0.146)
0.39–0.43 (1.410 ± 0.016+0.096−0.068) (1.778 ± 0.013+0.126−0.071) (1.983 ± 0.015+0.146−0.074) (2.415 ± 0.061+0.351−0.146)
0.43–0.46 (1.396 ± 0.016+0.094−0.068) (1.754 ± 0.013+0.123−0.070) (1.959 ± 0.014+0.143−0.073) (2.495 ± 0.060+0.362−0.145)
0.46–0.50 (1.391 ± 0.015+0.093−0.067) (1.715 ± 0.012+0.119−0.069) (1.941 ± 0.014+0.140−0.072) (2.383 ± 0.056+0.345−0.145)
0.50–0.54 (1.349 ± 0.014+0.089−0.066) (1.677 ± 0.012+0.115−0.068) (1.844 ± 0.013+0.132−0.071) (2.200 ± 0.052+0.318−0.144)
0.54–0.58 (1.328 ± 0.014+0.086−0.065) (1.643 ± 0.011+0.111−0.068) (1.789 ± 0.012+0.126−0.071) (2.213 ± 0.050+0.319−0.144)
0.58–0.62 (1.300 ± 0.009+0.084−0.064) (1.577 ± 0.010+0.105−0.067) (1.731 ± 0.013+0.121−0.070) (2.107 ± 0.032+0.303−0.144)
0.62–0.67 (1.254 ± 0.008+0.080−0.064) (1.527 ± 0.009+0.101−0.066) (1.641 ± 0.012+0.114−0.069) (1.966 ± 0.030+0.282−0.143)
0.67–0.72 (1.220 ± 0.008+0.077−0.063) (1.468 ± 0.009+0.096−0.065) (1.587 ± 0.011+0.109−0.068) (1.870 ± 0.028+0.268−0.143)
0.72–0.77 (1.173 ± 0.007+0.073−0.062) (1.403 ± 0.008+0.091−0.065) (1.510 ± 0.011+0.102−0.068) (1.769 ± 0.026+0.252−0.143)
0.77–0.83 (1.125 ± 0.007+0.069−0.061) (1.341 ± 0.007+0.086−0.064) (1.431 ± 0.010+0.096−0.067) (1.710 ± 0.025+0.243−0.142)
0.83–0.89 (1.064 ± 0.006+0.064−0.060) (1.272 ± 0.007+0.080−0.063) (1.357 ± 0.009+0.090−0.066) (1.550 ± 0.022+0.220−0.142)
0.89–0.96 (1.003 ± 0.006+0.060−0.060) (1.197 ± 0.007+0.075−0.062) (1.273 ± 0.009+0.084−0.066) (1.442 ± 0.021+0.204−0.142)
0.96–1.02 (0.964 ± 0.005+0.057−0.059) (1.135 ± 0.006+0.070−0.062) (1.192 ± 0.008+0.077−0.065) (1.382 ± 0.020+0.195−0.141)
1.02–1.09 (0.919 ± 0.005+0.053−0.058) (1.066 ± 0.006+0.065−0.061) (1.127 ± 0.007+0.072−0.064) (1.294 ± 0.018+0.183−0.141)
1.09–1.17 (0.860 ± 0.005+0.049−0.057) (0.993 ± 0.005+0.060−0.060) (1.041 ± 0.007+0.066−0.064) (1.178 ± 0.017+0.166−0.141)
1.17–1.25 (0.813 ± 0.004 ± 0.046) (0.934 ± 0.005 ± 0.056) (0.983 ± 0.006 ± 0.062) (1.105 ± 0.016 ± 0.155)
1.25–1.34 (0.753 ± 0.004 ± 0.042) (0.860 ± 0.005 ± 0.051) (0.903 ± 0.006 ± 0.056) (1.030 ± 0.014 ± 0.144)
1.34–1.42 (0.708 ± 0.004 ± 0.039) (0.813 ± 0.004 ± 0.047) (0.845 ± 0.006 ± 0.052) (0.931 ± 0.013 ± 0.130)
1.42–1.52 (0.650 ± 0.004 ± 0.035) (0.748 ± 0.004 ± 0.043) (0.782 ± 0.005 ± 0.048) (0.862 ± 0.012 ± 0.120)
1.52–1.62 (0.607 ± 0.003 ± 0.033) (0.693 ± 0.004 ± 0.039) (0.719 ± 0.005 ± 0.044) (0.789 ± 0.011 ± 0.110)
1.62–1.72 (0.562 ± 0.003 ± 0.030) (0.640 ± 0.003 ± 0.036) (0.670 ± 0.004 ± 0.040) (0.720 ± 0.010 ± 0.100)
1.72–1.83 (0.520 ± 0.003 ± 0.027) (0.590 ± 0.003 ± 0.033) (0.612 ± 0.004 ± 0.036) (0.675 ± 0.010 ± 0.094)
1.83–1.95 (0.481 ± 0.003 ± 0.025) (0.540 ± 0.003 ± 0.030) (0.562 ± 0.004 ± 0.033) (0.605 ± 0.009 ± 0.084)
1.95–2.07 (0.441 ± 0.002 ± 0.023) (0.501 ± 0.003 ± 0.027) (0.512 ± 0.004 ± 0.030) (0.547 ± 0.008 ± 0.076)
2.07–2.20 (0.412 ± 0.002 ± 0.021) (0.450 ± 0.002 ± 0.024) (0.461 ± 0.003 ± 0.027) (0.499 ± 0.008 ± 0.069)
2.20–2.33 (0.374 ± 0.002 ± 0.019) (0.409 ± 0.002 ± 0.022) (0.428 ± 0.003 ± 0.024) (0.450 ± 0.007 ± 0.062)
2.33–2.48 (0.347 ± 0.002 ± 0.017) (0.377 ± 0.002 ± 0.020) (0.388 ± 0.003 ± 0.022) (0.411 ± 0.006 ± 0.057)
2.48–2.62 (0.315 ± 0.002 ± 0.015) (0.348 ± 0.002 ± 0.018) (0.356 ± 0.003 ± 0.020) (0.390 ± 0.006 ± 0.054)
2.62–2.78 (0.286 ± 0.002 ± 0.014) (0.318 ± 0.002 ± 0.016) (0.322 ± 0.002 ± 0.018) (0.334 ± 0.005 ± 0.046)
2.78–2.94 (0.262 ± 0.002 ± 0.012) (0.288 ± 0.002 ± 0.015) (0.291 ± 0.002 ± 0.016) (0.303 ± 0.005 ± 0.042)
2.94–3.12 (0.239 ± 0.001 ± 0.011) (0.258 ± 0.002 ± 0.013) (0.262 ± 0.002 ± 0.014) (0.276 ± 0.005 ± 0.038)
3.12–3.30 (0.215 ± 0.001 ± 0.010) (0.234 ± 0.001 ± 0.012) (0.239 ± 0.002 ± 0.013) (0.251 ± 0.004 ± 0.034)
3.30–3.49 (0.195 ± 0.001 ± 0.009) (0.210 ± 0.001 ± 0.010) (0.215 ± 0.002 ± 0.012) (0.223 ± 0.004 ± 0.031)
3.49–3.69 (0.175 ± 0.001 ± 0.008) (0.189 ± 0.001 ± 0.009) (0.194 ± 0.002 ± 0.010) (0.204 ± 0.004 ± 0.028)
3.69–4.12 (151.4 ± 0.7 ± 6.9) × 10−3 (163.3 ± 0.7565 ± 8.050) × 10−3 (164.7 ± 0.9853 ± 8.790) × 10−3 (0.172 ± 0.002 ± 0.023)
4.12–4.59 (124.4 ± 0.6 ± 5.6) × 10−3 (132.0 ± 0.6391 ± 6.446) × 10−3 (134.3 ± 0.8364 ± 7.113) × 10−3 (0.141 ± 0.002 ± 0.019)
4.59–5.11 (99.15 ± 0.5 ± 4.4) × 10−3 (106.8 ± 0.5402 ± 5.168) × 10−3 (109.1 ± 0.7083 ± 5.731) × 10−3 (0.113 ± 0.002 ± 0.015)
5.11–5.68 (81.38 ± 0.4 ± 3.6) × 10−3 (85.34 ± 0.4537 ± 4.093) × 10−3 (86.97 ± 0.5944 ± 4.536) × 10−3 (0.091 ± 0.001 ± 0.012)
5.68–6.30 (65.27 ± 0.4 ± 2.9) × 10−3 (68.61 ± 0.3820 ± 3.262) × 10−3 (69.55 ± 0.4994 ± 3.601) × 10−3 (0.069 ± 0.001 ± 0.009)
6.30–6.99 (52.12 ± 0.3 ± 2.3) × 10−3 (54.86 ± 0.3211 ± 2.587) × 10−3 (55.60 ± 0.4197 ± 2.859) × 10−3 (56.69 ± 0.9533 ± 7.693) × 10−3
6.99–7.74 (41.66 ± 0.3 ± 1.8) × 10−3 (43.65 ± 0.2691 ± 2.042) × 10−3 (44.15 ± 0.3514 ± 2.255) × 10−3 (44.91 ± 0.7967 ± 6.088) × 10−3
7.74–8.57 (33.18 ± 0.2 ± 1.4) × 10−3 (34.31 ± 0.2241 ± 1.593) × 10−3 (34.65 ± 0.2923 ± 1.759) × 10−3 (34.70 ± 0.6572 ± 4.700) × 10−3
8.57–9.48 (26.26 ± 0.2 ± 1.1) × 10−3 (26.80 ± 0.1856 ± 1.235) × 10−3 (27.12 ± 0.2424 ± 1.368) × 10−3 (27.44 ± 0.5473 ± 3.713) × 10−3
9.48–10.48 (208.8 ± 1.5 ± 8.8) × 10−4 (214.3 ± 1.556 ± 9.872) × 10−4 (21.32 ± 0.2014 ± 1.075) × 10−3 (21.77 ± 0.4564 ± 2.945) × 10−3
10.48–11.57 (164.0 ± 1.2 ± 6.9) × 10−4 (167.1 ± 1.286 ± 7.697) × 10−4 (167.7 ± 1.672 ± 8.456) × 10−4 (16.86 ± 0.3755 ± 2.281) × 10−3
11.57–12.77 (128.9 ± 1.0 ± 5.5) × 10−4 (130.2 ± 1.061 ± 5.998) × 10−4 (134.6 ± 1.399 ± 6.785) × 10−4 (13.03 ± 0.3083 ± 1.763) × 10−3
12.77–14.09 (100.1 ± 0.8 ± 4.2) × 10−4 (103.4 ± 0.8845 ± 4.762) × 10−4 (102.8 ± 1.143 ± 5.182) × 10−4 (10.22 ± 0.2554 ± 1.384) × 10−3
14.09–15.54 (78.59 ± 0.7 ± 3.3) × 10−4 (79.72 ± 0.7250 ± 3.672) × 10−4 (78.55 ± 0.9319 ± 3.960) × 10−4 (8.339 ± 0.2150 ± 1.129) × 10−3
15.54–17.12 (62.71 ± 0.6 ± 2.7) × 10−4 (64.08 ± 0.6064 ± 2.952) × 10−4 (62.57 ± 0.7736 ± 3.154) × 10−4 (62.71 ± 1.735 ± 8.486) × 10−4
17.12–18.86 (47.92 ± 0.5 ± 2.0) × 10−4 (49.08 ± 0.4937 ± 2.261) × 10−4 (49.15 ± 0.6357 ± 2.478) × 10−4 (47.39 ± 1.400 ± 6.413) × 10−4
18.86–20.76 (37.16 ± 0.4 ± 1.6) × 10−4 (37.96 ± 0.4066 ± 1.748) × 10−4 (38.06 ± 0.5230 ± 1.919) × 10−4 (38.42 ± 1.181 ± 5.200) × 10−4
20.76–22.85 (28.92 ± 0.3 ± 1.2) × 10−4 (29.32 ± 0.3398 ± 1.351) × 10−4 (28.47 ± 0.4282 ± 1.435) × 10−4 (28.76 ± 9.688 ± 3.892) × 10−4
22.85–25.15 (225.8 ± 2.6 ± 9.6) × 10−5 (23.07 ± 0.2879 ± 1.063) × 10−4 (22.61 ± 0.3642 ± 1.140) × 10−4 (21.59 ± 0.8009 ± 2.922) × 10−4
25.15–27.66 (171.3 ± 2.2 ± 7.2) × 10−5 (173.8 ± 2.389 ± 8.006) × 10−5 (173.5 ± 3.049 ± 8.745) × 10−5 (18.35 ± 0.7055 ± 2.484) × 10−4
27.66–30.42 (138.0 ± 1.9 ± 5.8) × 10−5 (137.7 ± 2.031 ± 6.344) × 10−5 (138.8 ± 2.603 ± 6.997) × 10−5 (14.03 ± 0.5883 ± 1.898) × 10−4
30.42–33.44 (103.8 ± 1.5 ± 4.4) × 10−5 (105.6 ± 1.699 ± 4.862) × 10−5 (101.5 ± 2.127 ± 5.118) × 10−5 (10.95 ± 0.4965 ± 1.482) × 10−4
33.44–36.75 (81.42 ± 1.3 ± 3.4) × 10−5 (80.24 ± 1.416 ± 3.696) × 10−5 (81.51 ± 1.820 ± 4.109) × 10−5 (76.20 ± 0.3955 ± 1.031) × 10−4
36.75–40.39 (63.41 ± 1.1 ± 2.7) × 10−5 (62.97 ± 1.199 ± 2.900) × 10−5 (62.22 ± 1.520 ± 3.137) × 10−5 (61.51 ± 3.396 ± 8.323) × 10−5
40.39–44.37 (47.29 ± 0.9 ± 2.0) × 10−5 (49.30 ± 1.014 ± 2.271) × 10−5 (49.86 ± 1.300 ± 2.514) × 10−5 (47.38 ± 2.847 ± 6.412) × 10−5
44.37–48.74 (35.22 ± 0.8 ± 1.5) × 10−5 (33.67 ± 0.801 ± 1.551) × 10−5 (35.50 ± 1.049 ± 1.790) × 10−5 (36.99 ± 2.404 ± 5.006) × 10−5

Note. The first and second errors represent the statistical and systematic uncertainties, respectively.

Download table as:  ASCIITypeset images: 1 2

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 Foundation for Basic Research. Partial financial support was given by the South African National Research Council (NRF) for the numerical modeling.

Footnotes

  • 23 

    Mean synodic rotational period of the Sun surface, corresponding to about 27.28 days (Carrington 1863).

  • 24 

    An improved use of the simulation resulted in slightly higher selection efficiencies than previously estimated. Therefore, the proton spectrum of Adriani et al. (2011) should be lowered by 3.2%. Furthermore, the overall systematic uncertainty is slightly larger at low energies (4.5% at 1 GV).

Please wait… references are loading.
10.1088/0004-637X/765/2/91