HTML
--> --> -->In several countries, attempts are being made to introduce such projects in schools to stimulate interest in cosmic ray physics. In some countries, this endeavor has already been achieved on a smaller or larger scale, for example, the High School Project on Astrophysics Research with Cosmics (HiSPARC) [2] in the Netherlands, as well as WALTA [3], NALTA, ALTA [4], SALTA, CZELTA [5], SKALTA, CHICOS [6], CROP, CosMO [7], or Maze [8-10]. In several centers, work is in progress on the construction stage of prototypes or the research and development stage. From a technical persepective, such small arrays are not a serious challenge. A few typical scintillation detectors connected to a simple triggering, monitoring, and recording device would address this problem. However, because the idea of such projects is to popularize such arrays, and preferably to combine them into one large system, the cost of appropriate equipment for one school plays a significant role. In addition, the problem is complicated by the fact that teachers must be prepared to actively use the detectors and the entire network as part of their curricula and during after-hour activities. For this to work in practice, it is necessary to create a robust set of educational materials for both students and teachers. Students must be aware that they are carrying out observations of real physical phenomena. In addition, they should be aware that what they are observing are signals of cosmic rays reaching the Earth from distant cosmic spaces.
To achieve such awareness, they must have the tools required to interpret the signals recorded by their apparatus. The mere simultaneous appearance of some signals in several detectors is not fascinating to most, although it may be interesting to some. It starts to get really interesting when something tangible finally starts to emerge from such a phenomenon, or rather from a set of such phenomena, and from a registration that lasts continuously for days, weeks, or even months or years, which is actually the underlying knowledge of the surrounding reality. To go from direct measurement to its interpretation, appropriate tools are required, and such tools should be provided to students and teachers, together with clear and comprehensible documentation and user instructions.
Small instruments are inherently designed to record small showers. Certainly, it is possible to set the triggering thresholds high, and then once in a while after a long time period, we will observe a large EAS. Because, in our opinion, a very important component of the price of such a small shower array are scintillators, in the optimal version, they are possibly not going to be very large. There are also other considerations suggesting that the detectors should be "mobile" and able to work in different geometries, i.e., relatively small in size. In the case of small detectors, a trigger threshold of one minimal ionizing particle is necessary if we want to maintain the observed shower rate as high as possible. With such thresholds, the registrations will practically contain cases where the detectors register exactly one particle in coincidence. Clearly, we cannot omit the fact that a shower of primary energy 10
Such small showers obviously pose fundamental problems of interpretation. It is definitely impossible to localize the shower axes. In addition, with such a small number of detectors (3 or 4), localization, even at higher particle densities on detectors (which could be the case of larger showers), would be questionable regardless. For 3 – 4 detectors set up in small distances (approximately 10 m), it is also unrealistic to determine the directions of arrival of the showers. First, it would require precise measurements of time differences in the order of several nanoseconds, which would be possible in principle; however, even this approach would provide nothing tangible, considering the disk thickness of the small shower, which determines the fluctuations in the moment of appearance of the signal. In the absence of direction and axis position, event-by-event analysis does not provide any meaningful result. We can feasibly measure the observation frequency of various coincidences initially; then, if we equip the apparatus with amplitude measurement, we can determine the distribution of these amplitudes. The measurement times will allow us to study how the frequency of registration changes with time, which can be a basis for various studies conducted by groups of students, such as those on variability as a function of atmospheric parameters, day/night variability, and dependence of various "everyday life" quantities on the observed cosmic ray flux. To the best of our knowledge, the dependence of the Covid-19 infection rates on the cosmic ray flux is yet to be investigated.
From a fundamental physics perspective, it might be very interesting to combine several (the more, the better) instruments into one network and determine the existence of correlations in real time. The observation of large shower pairs at distant locations could indicate the existence of the Gerasimova-Zatsepin effect [13, 14] in the cosmic flux at the highest energies of heavy nuclei. The cosmic ray mass composition in this area remains a mystery waiting to be solved. The observation of other correlations could suggest the existence of more exotic objects, such as the Cosmic Ray Ensembles sought by the CREDO Collaboration.
In conclusion, building small, local, school-based EAS arrays, and networking them on a significantly larger, global scale has many positive aspects.
The analysis of data from small arrays requires specific methods. Methodological limitations do not allow to the adoption of conventional cosmic ray experiment methods for processing data. The interpretation of EAS registrations, in any case, is based on computer simulations describing the development of the shower and also on modeling the response of the array detectors to then compare the registrations with the predictions. The predictions incorporate parameters such as the mass of the original particle, its energy, and its arrival direction, along with the position of the shower axis (the intersection of the particle trajectory with the observation plane). The importance of these parameters is obvious. Apart from them, models and parameters of cosmic ray particle-air interaction also enter in an implicit way. These processes are not fully understood. Their descriptions are based on theoretically validated extrapolation of data of accelerator experiments.
The modeling of strong interactions has a long history. All models in the market today were developed a long time ago, and are constantly being refined and improved with the release of new data. These models need to be implemented into the geometrical structure of the transport through the Earth's atmosphere. Complementing them with a better known formalism of electromagnetic cascades and descriptions of other relatively important processes leads to the processes simulating the development of EAS. Today, one of the most widely adopted programs for this purpose is the CORSIKA program [15, 16], which was developed over 30 years ago in Karlsruhe for the KASCADE experiment [17, 18]. Since then, this program has been significantly extended and developed and is now also used for simulations at the highest observed energies (even up to 10
However, local arrays need to be simulated with energies at the other end of the cosmic-ray energy spectrum. This spectrum is very steep, and has a power-law character with an exponent (differential) of -2.7, which means that there are approximately 1000 times fewer particles with energies that are ten times larger. The low-energy end of the spectrum is truncated around energies of a few GeV/n, which is primarily owing to solar modulation. Certainly, atomic nuclei with kinetic energies of 1 GeV are not capable of causing what we might call an EAS or even a shower at all. The products of the interactions initiated by such particles would have no feasible chance of reaching the earth's surface.
The natural question, of course, is whether anything can reach the observation level (sea level) from a flux of cosmic ray particles with energies of 10 GeV/n or 100 GeV/n, at the top of the atmosphere. The only way to find an answer is to run appropriate simulations, which can be performed using CORSIKA. There is no limitation to determining the feasibility of observing a charged particle remaining from a proton initiated cascade, starting with the energy of 10 GeV. The possibility of observing anything is significantly negligible; however, we should recall that the flux of such particles is substantially large.
At energies of primary particles approximately 1 TeV and higher, we can expect that the particles from the showers initiated with them will arrive at the surface of the ground. The geometry of the shower development indicates that these particles will be scattered across large distances from the direction of the primary particle (from the axis of the shower). This is confirmed by simulations. These distances are on the order of hundreds of meters. In addition, a crucial complexity of the simulation challenges of showers at the lower end of the spectrum appears here. The densities of particles rarely exceed several per square meters. Such densities would require correspondingly larger primary energies, and as aforementioned, the energy spectrum is very steep.
2
A.Shower size
The first and undoubtedly most important parameter of a shower is its size, which is understood here as the number of particles at the observation level, as it is typical for surface shower arrays. In experiments analyzing the light (Cherenkov or fluorescent light) produced when a shower passes through the atmosphere, the shower size is defined as the number of charged particles in the shower maximum. CORSIKA allows us to precisely count the number of tracked electrons (and positrons), as well as positive and negative muons, that eventually reach the observation level. Certainly, all thinning options in the program have to be switched off. Repeating simulations continuously for the same particle with the same energy and the same angle of arrival, we will obtain a different result each time. Fluctuations in the size of the shower are an intrinsic property of the simulations, and they result from the probabilistic nature of multi-particle production processes. Figure 1 presents examples of electron and muon shower size spreads for cases of showers initiated by protons with energies of 10![](http://cpc.ihep.ac.cn/fileZGWLC/journal/article/zgwlc/2021/8/PIC/CPC-2021-0124-1.jpg)
For showers smaller than those illustrated in Fig. 1, when the sizes are approximately a few particles, the Log-normal distribution is obviously no longer the best. At small number of particles, we should expect an important correlation between them triggered by their common origin from subcascades, which developed accidentally just above the level of observation. In addition, in such cases, the Poisson-type statistical fluctuations are superimposed on fluctuations in their number caused by the probabilistic nature of the phenomenon itself.
The discrete nature of the size variable means that we expect a number of shower cases with zero charged particles. The analysis of the simulation outputs indicates that there are more of such cases than would originate from the Poisson nature of the process. The fraction of “empty” showers increases rapidly with a decrease in the energy of the primary particles. This relationship is important if you want to integrate the observed fluxes of particles appropriately. Obviously, this truncates the flux of primary cosmic ray particles on the low energy side. Fig. 2 shows this cut for vertical showers. It is represented as black dots for the results of simulations with the CORSIKA program for primary protons. As can be observed, the truncation for all charged particle size starts to work below the energy of 10
![](http://cpc.ihep.ac.cn/fileZGWLC/journal/article/zgwlc/2021/8/PIC/CPC-2021-0124-2.jpg)
A comparison of the truncation for the electron (Fig. 2a)) and muon (Fig. 2b)) sizes indicates that for very low energies, muons are the particles that manage to reach the ground, and they determine the counts in the individual detectors. In the following, we will analyze this effect quantitatively.
As already mentioned, the size of the shower initiated by a primary particle of a given energy is the most important (for integration of surface particle flux) parameter, which the CORSIKA program provides. We can define the size of simulated showers separately as electron and muon sizes, as the number of electrons or muons at the observation level. Both of these quantities are important for our purposes, and both are further analyzed in parallel.
Figure 3 presents the average values obtained from simulations for primary protons (filled symbols), where the proper energy scale is represented on the bottom axis and the values on the left axis, and for the iron nuclei empty symbols, on the top and right scales. First, it is important to note that, practically, in the entire energy range presented from 10
![](http://cpc.ihep.ac.cn/fileZGWLC/journal/article/zgwlc/2021/8/PIC/CPC-2021-0124-3.jpg)
2
B.Superposition model
First, the simplest assumption concerning the relationship between quantities describing showers initiated by protons and by the complex atomic nuclei is the assumption of simple superposition. According to this concept, a nucleus is considered a set of single nucleons, which behave like protons in their interaction with the nuclei of atoms of the Earth's atmosphere. Consequently, a shower initiated by, for example, an iron nucleus is the same as 56 proton showers. This assumption is relatively natural and correct to a large extent, as can be observed in Fig. 3, where the average sizes of proton showers are compared with the sizes of iron showers divided by 56 and shifted on the energy scale to the same energy per nucleon.The superposition model was also adopted to compare the number of showers with no particle at the observation level in the proton and iron events illustrated in Fig. 2. The possibility that none of the
Both of these observations (of the CORSIKA results) indicate that the superposition assumption is correct; however, this is not entirely and exactly true. Another important characteristic of EAS, which is indispensable for carrying out correct calculations of particle fluxes in small showers, is the size of fluctuations of shower sizes at small energies of particles that are initiating them. Examples of such fluctuations are presented in Fig. 1. As we have already demonstrated, in the first moments, the average values of these distributions agree with the superposition assumption; however, in the case of the second moments, this phenomenon is not as significant.
Figure 4 presents the dispersion of the logarithm of the electron (a) and muon (b) size distributions as a function of the energies of the protons and iron nuclei initiating the showers. The energy scales, abscissa, at the top and bottom are for the iron and proton showers, respectively, and the ordinates are depicted on the left for proton and on the right for iron showers (see Fig. 3). Energy scales for iron showers correspond to the same energy per nucleon as respective scales for proton induced showers. The ordinate for iron showers (right) is scaled down by
![](http://cpc.ihep.ac.cn/fileZGWLC/journal/article/zgwlc/2021/8/PIC/CPC-2021-0124-4.jpg)
As can be observed, the points from simulation calculations with the CORSIKA program for proton initiated showers do not overlap with the correspondingly shifted values for iron showers. We ignore the moment of discrepancies for very small showers with sizes below approximately 10 particles (for irons below
2
C.Radial distribution of small shower particles
It has been known for many years that the transverse distributions of particles in extensive air showers are well described by a simple formula proposed by Greisen [19]. Its validity was confirmed by the theoretical considerations and numerical calculations with respect to electromagnetic cascades by Kamata and Nishimura [20]; hence, its commonly accepted name: Nishimura-Kamata -Greisen (NKG) function $ \rho_{e/\mu}(r) = {N_{e/\mu} \over 2 \pi r_0^2}\: {\Gamma(4.5-s) \over \Gamma(s)\:\Gamma(4.5-2s)}\left({r\over r_0}\right)^{s-2} \left({1+{r\over r_0}}\right)^{s-4.5}\; \; \; , $ ![]() | (1) |
Figure 5 presents transverse distributions of electrons and muons in small vertical showers initiated by protons and iron nuclei with relatively high energies (10
![](http://cpc.ihep.ac.cn/fileZGWLC/journal/article/zgwlc/2021/8/PIC/CPC-2021-0124-5.jpg)
![](http://cpc.ihep.ac.cn/fileZGWLC/journal/article/zgwlc/2021/8/PIC/CPC-2021-0124-6.jpg)
For showers containing an average of approximately 1 particle, which corresponds to a (mean) central density of 10
For the purposes of this study, the lack of a statistically verified mean distribution of particles in very small showers is not particularly important, because the integrals of these distributions we eventually aim for will be normalized by the total number of particles in such showers, which will be approximately 0.01, as shown in Fig. 3.
for shower electrons and muons, separately. These quantities were obtained for different primary energies of cosmic ray particles, for a few (zenith) angles of the arrival of the particle, and for primary protons and iron nuclei
2
A.Superposition
First, we verified a simple superposition model that allows us to discard the mass of the primary particle as an independent variable from our descriptions. As we have already demonstrated, this model works well forThe simplest introduction of correlations among
It is well known that the main source of shower fluctuations is the height of the first interaction, which is related to the cross section of the interaction; hence, it is practical and theoretically justified to assume that the number of identical sub-shower components is proportional to the number of wounded nucleons in the interaction of the cosmic ray nucleus with the atomic nucleus of the atmosphere. Certainly, we are only referring to the wounded nucleon of the beam nucleus
In the case of iron and the atmospheric nucleus of nitrogen, the average number of wounded nucleons at the energies of interest in iron is approximately 12. Because we intend to use our superposition model also for other nuclei, we have approximately parameterized the number of wounded nucleons as a function of mass number, as presented in Fig. 7.
![](http://cpc.ihep.ac.cn/fileZGWLC/journal/article/zgwlc/2021/8/PIC/CPC-2021-0124-7.jpg)
The calculations performed indicate that the function presented in Fig. 7 does not provide a sufficient increase in the dispersion of the distribution (logarithmic) of the number of electrons and muons in iron showers. Because the assumption that the number of identical sub-showers depends on the number of wounded nucleons seems to be well-founded, in the next step, we assumed that it is not exactly equal but proportional to the number of wounded nucleons leaving the proportionality factor as a free parameter to be adjusted to the CORSIKA results. The actual integer number of identical sub-showers in each nucleus induced shower fluctuates according to a binomial distribution.
The results of the modified superposition model are shown in Fig. 4 for electrons and muons. As can be observed, in the region where the shower sizes are sufficiently large, we have obtained a sufficient agreement with CORSIKA simulations. As already mentioned, for very small showers, the integer number of particles in single showers and the combination of the physical spread with the effect from the Poisson distribution play a dominant role.
2
B.Small shower generator
The modified superposition method provides an opportunity to determine the average size (electrons and muons) of the shower initiated by any nucleus from H (protons) to Fe (iron), which dominate the cosmic ray flux of energies in the range of our interest from 10![](http://cpc.ihep.ac.cn/fileZGWLC/journal/article/zgwlc/2021/8/PIC/CPC-2021-0124-8.jpg)
There is a weak positive correlation between moun and electron sizes of CORSIKA showers, which is represented in Fig. 8 by solid lines and introduced in our small shower generator.
The dependence of the average number and dispersion of electrons and muons on the energy of the primary proton, shown in the figures above, was obtained for vertical showers. Simulations for inclined showers allowed to parametrize all these dependencies as functions of the arriving particle zenith angle.
For the particular size, the radial distribution of shower particles (for electrons and muons separately, with different scale paramerer
These findings are already sufficient to formulate a simple algorithm generating small EAS that will resemble the EAS generated by the CORSIKA program. Examples of showers generated in this manner are presented in Fig. 9.
![](http://cpc.ihep.ac.cn/fileZGWLC/journal/article/zgwlc/2021/8/PIC/CPC-2021-0124-9.jpg)
Eventually, the procedure is developed to generate the particular density of shower particles (electrons and muons) in the shower initiated by any cosmic ray primary particle of any energy, pointing to any point on the observation plane at any shower inclination. This allows to compare the results of our small shower generator, not only with the CORSIKA showers that it is based on but also with experimental results.
A.Shower particles density spectrum
The shower particles density spectrum was measured since at least the middle of the last century. The form of the spectrum determined agreed with a simple power law formula, for example the one measured by Cocconi, Loverdo, and Tongiorgi in 1946 for densities from approximately 10 to 1000 particles per meter squared [21, 22] $ N(x) = 700 \times x^{-1.47} , $ ![]() | (2) |
With our fast small shower generator, which not only reproduces the average shower characteristics determined for a fixed energy of the primary particle (including its mass and angle of arrival) but also considers respective dispersions, the multidimensional Monte Carlo integration can be carried out to obtain the shower particle density observed with a single small detector. The general formula is quite trivial:
$ \begin{aligned}[b] f(\varrho_{e/\mu}) =& \sum_A \: \int \limits_{10^9}^{10^{16}} {\rm d}E \:\: \Phi_A(E)\: \int\limits_0^{90}{\rm d}\phi\: \int {\rm d}{N_{e/\mu}}\:\\&\times p\left( N_{e/\mu},\langle N_{e/\mu}\rangle\right) (E,A,\phi) \\ &\times \:\int\limits_{\rm -2km}^{\rm 2km}{\rm d}x \:{\rm d}y\: 2\pi\:\sin(\phi)\:\cos(\phi)\:\varrho(r) , \end{aligned} $ ![]() | (3) |
The charged particle density is
The upper limit of the primary particle energy spectrum depends on the values of densities we attempt to study. This study is concerned with small densities and small showers; hence, if we limit ourselves to densities not exceeding a few hundred per m
Summing over the mass spectrum of cosmic ray particles, for practical reasons, as it is conventionally practiced, was changed to summing over a few groups of particles with similar masses: protons (
Certainly, if the total observed flux of primary cosmic radiation solely comprises iron nuclei, the flux of muons, as well as electrons, would be significantly larger than that of a purely proton composition; however, the differences in particle intensities at sea level for mixed mass spectrum compositions, PD and HD, are small, the calculated total muon flux at the sea level is 97 m
The integration results in Eq. (3) are depicted in Fig. 10, in comparison with measured results presented above [22-24]. As can be observed, the agreement exhibited is very good.
![](http://cpc.ihep.ac.cn/fileZGWLC/journal/article/zgwlc/2021/8/PIC/CPC-2021-0124-10.jpg)
Studies on the particle density spectra allow, to some extent, conclusions to be drawn on the energy spectrum of cosmic ray particles at the top of the atmosphere. In addition, averaging over the position of the shower axis (
2
B.1 m$ ^2 $![]()
![]()
detector
By integrating the spectra illustrated in Fig. 11 , we can obtain the fractions of corresponding observations and, for example, the rate of registering single muons.![](http://cpc.ihep.ac.cn/fileZGWLC/journal/article/zgwlc/2021/8/PIC/CPC-2021-0124-11.jpg)
If we measure the number of particles observed on a 1 m2 detector, each observed value corresponds to a different distribution of the primary particle energy. The results of the calculation are presented in Fig. 11, where the energy spectra of the primary particles leading to observations of 1, 2, 3, 4, and 5 electrons (solid lines) and muons (dashed lines) are illustrated. These results lead to some important conclusions:
By integrating the spectra shown in Fig. 11 , we can obtain the fractions of corresponding observations and, accordingly, the rate of registering single muons.
2
C.Small EAS array
As already mentioned in Sec. I, the main application of the small shower generator is to assist in the interpretation of data from small shower arrays, either for educational purposes or for its applications in integrated networks of such local stations.For example, let us assume that such stations would consist of four identical detectors located not far from each other (at a distance of 5 meters). We assume that the detector is ideal and generates no noise, which allows the interpretation of events when only one of the detectors has registered anything, usually one particle. The practical importance of cases where one detector registers several particles according to the results shown in Fig. 11 is negligible. Although these events are very rare, such cases can be analysed via the simulations of the small shower generator.
Let us further assume that every registration of a particle by any detector will trigger an event and that the binary state (hit/no hit) of each detector will be stored.
The small shower generator will help in answering the question on the energy of the primary particle required, or more precisely, the question on the energy distribution that should be associated with a given type of coincidence. Examples of such results are presented in Fig. 12.
![](http://cpc.ihep.ac.cn/fileZGWLC/journal/article/zgwlc/2021/8/PIC/CPC-2021-0124-12.jpg)
With the spectra illustrated in Fig. 12 , we can determine the rate of particular coincidence. The obtained results are presented in Table 1.
coincidence | charged | muons only | |||
1 m2 | 0.02 m2 | 1 m2 | 0.02 m2 | ||
single | 433 | 2.0 | 331 | 1.3 | |
2-fold | 5.6×10?1 | 2.5×10?3 | 2.7×10?2 | 1.5×10?5 | |
3-fold | 1.1×10?1 | 5.2×10?4 | 1.6×10?3 | 4.3×10?7 | |
4-fold | 6.0×10?2 | 2.9×10?4 | 4.1×10?4 | 4.4×10?8 |
Table1.Rates (s?1) of single detector registration as well as 2-fold, 3-fold, and 4-fold coincidences for stations of four detectors of 1 m2 and 0.02 m2 each. Results for muons alone are presented in the last two columns.
2
D.Individual detector size
As mentioned in Sec. I, the size of detectors in local, school arrays is a very important parameter in their design. It seems that the 1 mThe effect of the dominance of the electron component for higher rank coincidences, which we have shown in Fig. 12b), is substantially important for the planning of small, school shower arrays and their detailed location. Detectors placed under concrete roofs in physics labs will record shower events significantly less frequently. To obtain a four-fold coincidence once an hour, detectors with a size of 10 cm
We used CORSIKA to obtain the characteristics of the shower from the lowest primary energies allowed by this program. The average values of the quantities of interest were determined, as well as the shapes of their distributions, their second moments, and small prevalent correlations. The adjustments were done separately for the soft (electrons) and hard (muons) components. Accordingly, the transverse distributions of these particles, as a function of the primary energy of cosmic ray particle, its mass, and the zenith angle of its arrival, have been parameterized to build our small shower generator.
By analyzing electron and muon shower sizes for different atomic masses of the primary nuclei, we have examined the simple superposition hypothesis, and we infer that the extensive air shower produced in the interaction of a nucleus of mass A is a simple compound of A showers initiated by the nucleons (protons) of appropriately lower energy. The behavior of the average size values agrees with this assumption. In addition, we have also determined the consistency with the superposition rule in the range of the smallest showers that do not contain any particle (electrons nor muons).
However, by investigating dispersions of the shower size (electrons and muons), we observed that the assumption that the shower from nucleons, folding according to the superposition principle, fold independently of each other, does not correspond to reality, or at least, does not correspond to the reality of the CORSIKA program. In our shower generator, we introduced a correction to address the inconsistency of the simple superposition model by including a correlation between some of the constituent subshowers, the magnitude of which is proportional to the number of wounded nucleons in the interacting nucleus of the primary cosmic ray.
Using our small shower generator, it is possible to perform fast integrations of secondary particle fluxes at sea level and predict the registrations made by small school, local EAS arrays, as well as single detectors. These predictions, when confronted with the measured values, allow a deeper analysis of the local measurements and the properties of the detectors themselves.
Examples of the results presented in this paper, such as the spectra of primary particles illustrated in Figs. 11 and 12, demonstrate that signals from single detectors with low triggering thresholds primarily originate from single muons produced by primary particles of very low energies reaching tens or several tens of GeV. In contrast, the electron component of showers with energies of approximately 10
The small shower generator also allows the count rate in small EAS arrays to be estimated and the detector size and number of detectors to be optimized accordingly.
A small shower generator, packaged in a user-friendly manner and supplemented with appropriate documentation and instructions for use, can be an excellent research tool for allowing young people engaged in the study of modern physics during planned educational projects.