# Cosmic ray primary mass composition above the knee: deduction from lateral distribution of electrons

###### Abstract

Influence of shower fluctuations on the shape of lateral distribution of electrons in EAS of fixed size measured by scintillation counters is analyzed in framework of scaling formalism. Correction factors for the mean square radius of electrons are calculated for the experimental conditions of KASCADE array. Possible improvement of the primary mass discrimination by analysis of lateral distribution of EAS electrons is discussed in detail.

\runtitleCosmic ray primary mass composition above the knee \runauthorR. I. Raikin et al

## 1 Introduction

The determination of the primary cosmic ray chemical composition from extensive air showers (EAS) observations is an open problem. Various techniques based on different EAS characteristics measured by different experiments including multi-component methods do not exhibit consistent results in estimation neither primary mass in the case of individual showers nor the mean mass composition at a certain energy.

One of the key EAS quantities necessary for basic shower parameters reconstruction is the lateral distribution of charged particles at fixed observation depth. The exact form of lateral distribution function (LDF) is still uncertain. The majority of analytical parameterizations of LDF of different EAS components is traditionally based on the well known Nishimura-Kamata-Greizen (NKG) function [1]:

(1) |

Here is the particle density at radial distance from the core position in shower with primary energy and the age parameter , – total number of particles at the observation depth, – shower scale radius, which does not depend on primary particle type and energy (originally – the Mòliere unit). Various modifications of NKG form, such as introducing an additional fixed or age-dependent scale coefficient or a local age parameter and also generalizations of the function by using third power-law term were suggested.

A different theoretically motivated approach (scaling formalism) was proposed in our papers [2, 3]:

(2) |

where is the observation depth. The scaling function can be described as follows:

(3) |

For electron densities we find , , , , – root mean square radius of electrons:

(4) |

According to our calculations, the scaling formalism allows to reproduce electron LDF with 10% uncertainty for eV, g/cm, . The last condition corresponds to the radial distance range from m to km depending on the shower age. This limitation makes scaling approach inadequate for shower size and core position estimation, but accurate enough for description of the shape of LDF measured by ground-based shower arrays far from the core.

The method for mean primary mass deduction based on scaling formalism was developed in [4, 5]. The advantage of this method is its applicability to the experimental data of both compact and giant air shower arrays in wide primary energy range and also relatively weak sensitivity of the conclusions to variations of basic parameters of hadronic interaction model implemented in calculations. Unfortunately, the shape of charged particle LDF measured experimentally is affected by the experimental method of shower classification. This effect can also be described in the framework of scaling formalism, but the variation of the root mean square radius compared with the data obtained theoretically for the fixed primary energy should be taken into account in case of relatively low primary energies, when shower selection is made by the total number of electrons, e.g. for example for KASCADE and Moscow State University air shower arrays.

In this paper we examine thoroughly the influence of shower fluctuations on the shape of lateral distribution of electrons in EAS of fixed size under the experimental conditions of KASCADE array in order to improve the reliability of mean primary mass deduction from the experimentally measured LDFs.

## 2 Calculation Methods

We made simulations of extensive air showers initiated by protons and iron nuclei of vertical incidence assuming power-law differential energy spectrum of primaries with exponent , and also with sharp knee from to at GeV. We used the semi-analytical code [5] with full Monte-Carlo treatment of hadronic part of cascade based on quark-gluon string model and analytical expressions of pure electromagnetic subshowers keeping all the basic sources of fluctuations. The fluctuations of different EAS components calculated by our code are in good agreement with CORSIKA/QGSjet results.

According to [6] we simulated shower classification procedure used at KASCADE array and evaluated lateral distributions and root mean square radiuses of electrons in eight bins of shower size. The number of showers in each bin amounts from for lower energies () to for higher energies ().

## 3 Results

The primary energy distributions in four from eight bins of shower size is shown in fig. 1. One can see that the energies largely overlap in different bins though the selected bins are not neighboring.

We evaluated the correction factors defined as the ratio of root mean square radius calculated for certain shower size bin to that for corresponding average primary energy: . These correction factors for vertical proton initiated showers at sea level calculated with above mentioned assumptions about primary energy spectrum are shown in fig. 2. It is clear, that values of approach to 1 with energy as shower fluctuations decrease. The same effect takes place for a heavier primary nuclear or for smaller shower size bins. At the same time the correction, which should be made for an adequate comparison of theoretical and experimentally estimated root mean square radiuses is essential for all considered shower size bins.

Function (3) | Modified NKG (s=1.65) | Polinomial (5) | |
---|---|---|---|

Assuming the validity of scaling approach we made one-parametric fitting of experimental lateral distribution of electrons obtained by KASCADE array [6] using formula (2) with three different scaling functions:

1) theoretically proved function (3) with mean square radius of electrons as radial scale parameter;

2) modified NKG-function [6] with fixed shower age parameter (s=1.65) and variable ;

3) polinomial function:

(5) |

with and parameters being fixed for all bins independently from . Some discrimination of experimental data at small radial distances was done in order to eliminate points with .

All the fitting functions give satisfactory overall fit of experimental data of KASCADE with residuals not exceeding 10% for polynomial and modified NKG functions and 15% for function (3). Though polynomial and NKG functions give considerably better accuracy in considered radial distance range, they both lead to incorrect predictions for extremely large core distances, while theoretically motivated function (3) remains realistic up to .

The values of radial scale parameters obtained by fitting the experimental data are summarized in Table 1. It is not surprising, that different fitting functions correspond to significantly different values of . An additional bias in can be related with the insufficiency of radial distance range well covered by the array or some other systematic errors in data processing. So it is worth to compare the rate of change of radial scale factors with energy which obviously reflects the rate of change of mean primary mass. In fig. 3 values of superposed with each other by the appropriate factors are shown in comparison with root mean square radius calculated for vertical proton and iron initiated showers at sea level taking into account the correction factor for primary spectrum with the knee. As it is seen from the figure the rate of change of obtained using different functions is consistent.

## 4 Conclusions

1. Correction by a factor should be made when comparison of lateral distributions of electrons measured by ground-based experimental arrays to LDFs calculated theoretically for fixed primary energy is carried out.

2. Absolute values of radial scale factor contain systematic errors and could be biased depending on the form of lateral distribution function chosen for experimental data processing and final fitting. However, if one concerns the rate of change of with primary energy, which is a good measure for primary mass composition variation, then different scaling functions used in the framework of scaling formalism do not contradict each other.

3. Basic analysis using different assumptions about the form of scaling LDF leads to consistent model insensitive conclusion that average primary particle mass above the knee increases with energy, that is in reasonable agreement with the large number of experimental observations (see e.g. [7]) and also with recent results of the anomalous diffusion model [8].

## References

- [1] K. Greizen. Ann. Rev. Nucl. Sci. 10 (1960), p.63.
- [2] A. A. Lagutin et al. Proc. 25 ICRC, Durban, 1997, Vol. 6, p. 285.
- [3] A. A. Lagutin, R. I. Raikin. Nucl. Phys. B (Proc. Suppl.), 97B (2001). P.274.
- [4] R. I. Raikin et al. Proc. 27th ICRC, Hamburg, Germany, 2001, V.1, P.290.
- [5] A. A. Lagutin et al. J. Phys. G: Nucl. Part. Phys., 28 N6 (2002), p.1259.
- [6] T. Antoni et al. Astroparticle Physics, 14 (2001), P.245.
- [7] J. R. Hörandel. Astropart. Phys, 19 (2003), P.193; arXiv:astro-ph/0407554 v1.
- [8] A. A. Lagutin, A. G. Tyumentsev. Proc. 29th RCCR, Moscow (2006); these proceedengs.