Results on the primary CR spectrum and composition reconstructed with the SPHERE-2 detector

Results on the primary CR spectrum and composition reconstructed with the SPHERE-2 detector

R A Antonov, S P Beschapov, E A Bonvech, D V Chernov, T A Dzhatdoev, Mir Finger, Mix Finger, V I Galkin, N V Kabanova, A S Petkun, D A Podgrudkov, T M Roganova, S B Shaulov and T I Sysoeva

First preliminary results of the balloon-borne experiment SPHERE-2 on the all-nuclei primary cosmic rays (PCR) spectrum and primary composition are presented. The primary spectrum in the energy range eV was reconstructed using characteristics of Vavilov-Cherenkov radiation of extensive air showers (EAS), reflected from a snow surface. Several sources of systematic uncertainties of the spectrum were analysed. A method for separation of the primary nuclei’ groups based on the lateral distribution function’ (LDF) steepness parameter is presented. Preliminary estimate of the mean light nuclei’ fraction at energies eV was performed and yielded = (2111) %.

  • Lomonosov Moscow State University Skobeltsyn Institute of Nuclear Physics, Moscow, Russia

  • P.N. Lebedev Physical Institute of the Russian Academy of Sciences, Moscow, Russia

  • Charles University, Prague, Czech Republic

  • Faculty of Physics, Lomonosov Moscow State University, Moscow, Russia

  • E-mail:

1 Introduction

Despite more than 50 years of the cosmic rays’ (CR) energy spectrum measurements, a considerable uncertainty of the spectral shape in the ultrahigh energy domain ( eV= 1 PeV) still exists (for recent results compilation see [abb12], fig. 15). The difference between spectra reconstructed by the various experimental groups is sometimes quite pronounced, and cannot be attributed to the statistical issues only. For example, the CR all-nuclei spectrum measured by the GAMMA group [gar08] exhibits statistically significant () sharp peak at 810 eV. A similar feature exists in the spectrum based on the 2009–2010 run of the Tunka-133 experiment [ant11a][ant12a], but for the 2010–2011 run of the same experiment it was not confirmed [ant12a]. Also, the peak is not seen by other CR experiments such as Akeno [nag92], Tibet-III [ame08] and KASCADE-Grande [ape11][ape12].

The situation with the primary CR composition at superhigh energies is even more enigmatic, as the consensus between the different measurements, such as [ant05][tok08], is virtually absent (see fig. 1 of [tsu08]). More recent review [kam12] claims that the results on the mean logarithmic mass number, , are still consistent with a trend of growing at energies 3-4 PeV, where the primary CR spectrum exhibits a pronounced break, the so-called knee [kul58]. But anyway, the scatter of values measured by the different experiments is high, and more detailed information than just is needed for the detailed comparison with the CR acceleration and propagation models.

These examples highlight the importance of the all-nuclei CR spectrum measurement and composition studies in various experiments and by the diverse methods. In present work we present the first results of the SPHERE-2 experiment. The SPHERE-2 experiment is a Cherenkov telescope situated above the snow surface of Lake Baikal and registers reflected Vavilov-Cherenkov radiation (”Cherenkov light”). The basic idea of such an experiment was proposed in [chu72].

2 The data

The SPHERE-2 experiment uses a mosaic of 109 PMTs for recording the light signal reflected from the spherical mirror. Observation level is typically 400–900 m above the snow surface. In the present work we use the data from the 2011 and 2012 experimental runs. Time sampling step is 25 ns for the 2011 run and 12.5 ns for the 2012 run. For more detailed description of the SPHERE-2 experiment see [ant11b][ant12b].

In the four flights of the 2011 run 230 events consistent with EAS Cherenkov light reflected from the snow surface were found. For the 2012 run we use the data taken in the first four flights, and the number of events classified as showers is 301. Relatively small number of the showers compared to the ground-based arrays is due to a high threshold of the EAS observation, as, in average, only 1 out of each 10 Cherenkov photons reflected from the snow surface reaches the detector’s mirror.

3 A method for the primary spectrum and composition reconstruction

In order to reconstruct the primary CR spectrum and composition we applied an approach based on detector response simulation by means of the full Monte Carlo method. The procedure of the all-nuclei spectrum measurement consists of the several stages.

1. As the first step, the full Monte Carlo simulation of the lateral distribution function (LDF) of the EAS Cherenkov light, as well as the time structure of the shower was performed using the standard code CORSIKA6.500 [hec98] (for more details see [ano09][dzh11]). All results presented here are based on the QGSJET-I high energy hadronic model [kal97] and the GHEISHA low energy hadronic model [fes85].

2. Then, the geometrical and optical effects smearing the observed signal were accounted for. For the 2011 run reconstruction we use the simplified detector response model developed in [ant09][dzh11] and written in C++. In the case of the 2012 run the more advanced model is used, which was implemented using the standard code Geant4 [ago03].

3. The simulation of the instrumental acceptance was performed. The result is most sensitive to the observation altitude’s change, as well as to the variations of the thresholds of the detector’s channels.

4. For the sample of the EASs registered in the experiment the primary energy was estimated. We use the method of [ded04], where the energy estimation was performed together with the axis position reconstruction by normalising of the experimental LDFs to the model LDFs with known energy.

5. The all-nuclei CR spectrum was reconstructed using the sample of the showers with estimated energy and results of the acceptance calculation.

The primary composition study in the present work is based on the LDF steepness parameter, defined as the ratio of the number of Cherenkov photons in the circle with the radius of 67 m to the same number in the ring with the radii of 67 m and 134 m [ant09]. The parameter allows to select some fraction of the light nuclei [ano09][ant09]. In present work protons stand for light nuclei, and Iron for heavy ones.

4 Preliminary results

The all-nuclei CR spectrum reconstructed using 2011–2012 runs data of the SPHERE-2 experiment is shown in fig. 1 (red stars). Statistical (red bars) and systematic (red dashed lines) uncertainties of the spectrum are shown. The main sources of systematic uncertainty of the spectrum are as follows: a) bin-to-bin migration due to the energy estimation error (dominant factor at E20 PeV); b) acceptance estimation inaccuracy (not more than 5 %); c) the spectal shape vs. the primary composition dependence (dominates at E20 PeV). With respect to the direct predecessor of the present work [ant97], we achieved a considerable progress of the reconstruction methodology. Also the results of some other experiments are shown in fig. 1a for comparison: Akeno [nag92] (green circles), KASCADE-Grande [ape11] (black triangles), and Tunka-133 [ber11] (blue squares). Statistical errors of the Akeno spectrum are comparable with the circle’s diameter; systematic uncertainties of the KASCADE-Grande result are shown by black dashed lines.

Estimated fraction of light nuclei vs. is shown in fig. 2 (red circles) with statistical uncertainties of the values (red bars). Several independent reconstructions of the experimental showers’ LDFs, and then the vs. dependence were performed, and induced systematic uncertainty of the value was estimated (black lines). In the first two bins the value of is strongly modified by the threshold effects; an additional uncertainty assotiated with them is shown in the picture by blue arrows. The data presented in fig. 2 makes it possible to estimate the mean (not intensity-weighted!) fraction of light nuclei in the energy region 30–150 PeV . The corresponding value of (comparable with values allowed by [kam12]; see fig. 17).

All presented results are subject to additional recheck and develompent of the methodology, and the detailed study of the systematic uncertainties that might change the results presented at figs. 1–2 is underway.

\hb@xt@ . c.eps \hb@xt@ . missing
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description