Модельный комплекс для исследования космических лучей
Автор: Маурчев Е.А., Балабин Ю.В.
Журнал: Солнечно-земная физика @solnechno-zemnaya-fizika
Статья в выпуске: 4 т.2, 2016 года.
Бесплатный доступ
В настоящее время существует значительное разнообразие методик для исследования характеристик потоков как первичных, так и вторичных космических лучей (КЛ). Основными являются экспериментальные методы, с использованием разных детекторов, однако, благодаря современным компьютерным технологиям, экспериментальные методы все чаще совмещаются с численными. Данный шаг позволяет получить наибольшее количество информации о интересующем процессе (явлении) и сделать, соответственно, наиболее точные выводы. В работе представлена концепция нового программного комплекса RUSCOSMIC, основанного на пакете GEANT4 и представляющего собой набор различных численных моделей для исследования прохождения КЛ через вещество различных систем (детекторы излучения, атмосфера Земли). Получены и приведены функции откликов основных детекторов излучения, а также некоторые типовые характеристики потоков вторичных КЛ. Также представлены результаты, демонстрирующие работу модуля верификации вычислений экспериментальными данными.
Космические лучи, техника эксперимента, численное моделирование, метод монте-карло, детекторы излучения, взаимодействие частиц с веществом
Короткий адрес: https://sciup.org/142103618
IDR: 142103618 | DOI: 10.12737/21289
Текст научной статьи Модельный комплекс для исследования космических лучей
Cosmic rays (CRs) are currently being studied in different ways. Initially these were mainly experimental methods widely used so far. Sounding balloons equipped with detectors are launched [Bazilevskaya et al., 2008], and the worldwide network of neutron monitors is one of the most advanced and reliable ground-based systems for detecting cosmic rays. Data from this network are used to determine characteristics of relativistic solar protons, recorded during events of solar cosmic rays (SCR), by simulating increases in the ground-level background of secondary cosmic rays (GLE events [Vashenyuk et al., 2009; Shea, Smart, 1982]).
From the beginning of computerization, a technique for studying CR propagation through a medium by the Monte Carlo method has been rapidly developed. We have devised a software package RUSCOS-MIC©, which includes models of detectors of different types as well as a model of primary CR propagation through Earth’s atmosphere, to gain information on a cascade of secondary cosmic rays. The software package is based on a toolkit GEANT4 [Agostinelli et al., 2003], which provides classes describing particle interaction with matter, geometry, selection of initial parameters, collection of information on simulation, particle conditions, etc.
At the first stage of this work, we examined the properties of detection equipment at the cosmic ray station located in the town of Apatity. For this purpose, we worked out software modules representing models of neutron monitor (NM), Geiger counters, scintillation detectors with geometry and properties of materials as close to reality as possible.
When simulating propagation of a neutron flux through a neutron monitor, we found a response function and examined properties of a decelerator and lead generator. Comparison of our results with results of earlier and currently ongoing works showed good agreement between them. Up-to-date interaction cross-sections used in modeling enabled a more detailed study of the response function of the neutron monitor in a range from 10 MeV to epithermal energies. The simultaneously performed modeling of multiplicity in the neutron monitor confirmed the assumption that the phenomenon is caused by a cloud of particles rather than a single particle entering the NM lead generator.
The modeling of scintillation gamma-quantum detectors became necessary because there is an X-ray monitoring system at the CR station, which has been operating for several years already. This system allowed us to reveal X-ray increases during precipitation. GEANT4 was used to calculate response functions of NaI detectors of different geometry. Of course, such works have also been carried out before, but our study features a modern approach and a completely new concept enabling us to take into account more factors in calculating propagation of particles through a medium.
The second step was to create our own software module to compute propagation of particles of different energies through Earth’s atmosphere. This allowed us not only to use continuously updated data for simulation, but also to make it more broadly applicable since we can easily change input parameters of the model itself (to change physical properties of matter, to add various fields, etc.). Then, we calculated propagation of galactic cosmic rays (GCR) through Earth’s atmosphere and development of particle cascades during some GLE events. The result is a database containing data on energy spectra of secondary cosmic rays for different heights. These calculations can be used both to measure ionization rate in estimating an equivalent dose and to examine development of cascades during GLE events and to find their new features. The simulation results were compared with the results of similar studies carried out with the PLANETOCOSMICS software package [Maurchev et al., 2011] as well as with experimental balloon data, which show good agreement with our results.
MODEL OF SCINTILLATION NAI DETECTOR
In detectors of this type, gamma-quantum energy and intensity are determined from secondary charged particles (electrons and positrons), which result from interactions of gamma-quantums with matter. As mentioned in Introduction, all models of the complex are based on the Monte Carlo method and a selected value x e ( x 1 , x 2 ) of the normalized probability density function f ( x ):
n f (x ) = Z af (x) gi(x), (1) i=1
where α i > 0 is the probability sample of the normalized probability density function f i ( x ) and 0≤ g i ( x ) ≤1. Then, having selected a random integer with a probability proportional to α i and having chosen x 0 from the distribution f i ( x ), we can calculate x = x 0 with a probability of g i ( x 0 ). If the value deviates, the calculation scheme is repeated from the beginning. The detailed mathematical description of the model is too cumbersome and therefore not given in this paper. The above concept was used to design two scintillation detectors with geometrical dimensions Ø 6.3×2 cm and Ø 15×10 cm. The simulation of propagation of monoenergetic gamma-quantum beams through the scintillation detectors with the above geometry allowed us to estimate gamma-quantum detection efficiency (Figure 1).

Figure 1. Gamma-quantum detection efficiency obtained from the simulation of interactions of particles with the scintillation detectors of standard sizes Ø6.3×2 ( a ) and Ø15×10 cm ( b )
NEUTRON MONITOR MODEL МОНИТОРА
Geometrical parameters of the model are identical to the parameters of the real NM, the model is illustrated in Figure 2; the standard tracking is shown in Figure 3. A flux of monoenergetic primary neutrons is incident perpendicular to the upper plane of the detector; the particles are uniformly distributed; the number of particles in the flux is 500 000 per energy value. The response function obtained from the simulation is shown in Figure 4.
SIMULATION OF CR PROPAGATION THROUGH EARTH’S ATMOSPHERE
One of the most complex model in the software package RUSCOSMIC© is a model for calculating propagation of primary CR protons through Earth’s atmosphere and for studying formations of secondary CR cascades.
The model is based on the concept of the so-called plane geometry when a column of air is separated at a given latitude and longitude, its sizes are specified, and physical parameters of each layer (percentage of masses of constituent elements, density, pressure, and temperature) is calculated using NRLMSISE-00 [Picone et al., 2002]. At the boundary, we determine a model particle source with a preset intensity and spectral characteristic. Figure 5 depicts the spectrum of primary CR pro-5
tons corresponding to the period of maximum solar activity in the absence of SCR. Parameters of secondary CR cascades obtained from the simulation of propagation of primary particles with such energy characteristics through Earth’s atmosphere are shown in Figures 6 and 7.
Vertical profiles are usually determined from data acquired by balloons with detectors sensitive to the charged component. The total count can be expressed by the formula
С = С\ + С\ + С+0.01-С, (2)
сумм e+e- и+ц- p . у ’ where Сe, Cи|| , Cp, CY are counts of electron-positron, muon, and proton components, and gammaquantums respectively.
Accordingly, summing up the vertical profiles determined from the simulation, we can compare them with experimental data for certain periods of time. An example of such verification is shown in Figure 8.
The sum of average and internal errors in the method when simulating the vertical profiles can conveniently be represented as a line width (see Figure 9).

Figure 2. The model of the standard NM, which consists of a polyethylene decelerator, a lead generator, a polyethylene reflector, and proportional 10BF 3 counters, illustrated to present the main components of the detector

Figure 3. A part of the front projection of the model of one section in the standard NM. A 300 MeV neutron falls on the polyethylene from above. Inelastic collision with lead (we employ the model of intranuclear cascades [Bertini, 1969;. Heikkinen et al, 2003], point “1”) causes neutron multiplication. We can clearly see processes of neutron drift with multiple elastic collisions (point “2”), stop, and absorption in the polyethylene (point “3”). All these processes are calculated using the ENDF/B-VII.1 toolkit [Chadwick et al., 2011]. The process of neutron capture by the counter is not illustrated here, but it should be noted that it is described by the standard reaction n + 10 B → 7 Li + α and is simulated also using ENDF / B-VII.1

Figure 4. NM response function obtained from the simulation of propagation of monoenergetic neutron fluxes through the detector system. It is compared with a function [Clem, Dorman, 2000] derived previously with the FLUKA toolkit – they show close agreement. The difference in the so-called function “tail” ( E = 10–7 –10–2 MeV) is likely to arise from the use of different models of interactions of low energy neutrons as well as different interaction cross-sections

Figure 5. The spectrum of primary GCR protons from 400 to 105 MeV (~ 1 GW geomagnetic cutoff rigidity) in the absence of its SCR modulation, used as one of the model particle source characteristics when calculating secondary CR fluxes in Earth’s atmosphere

Figure 6. Energy spectra of protons (left) and gamma-quantums (right) of secondary CR obtained from the simulation of propagation of primary CR particles with the spectrum corresponding to the period of maximum solar activity in the absence of its SCR modulation through Earth’s atmosphere

Figure 7. Vertical profiles of electron-positron (left) and muon (right) secondary CR components obtained from the simulation of propagation of primary CR with the spectrum corresponding to the period of maximum solar activity in the absence of its SCR modulation through Earth’s atmosphere

Figure 8. Comparison between experimental data (September 04, 2009 (left) and December 21, 2009 (right)) and data obtained from the simulation of GCR propagation through Earth’s atmosphere and processed in accordance
with Formula (2)

Figure 9. Vertical profiles of particle fluxes presented with a calculation error range, which were obtained from the simulation of GCR propagation through Earth’s atmosphere. The line width indicates the total deviation from mean and includes both the standard and internal errors of the method
CONCLUSION
In this work, we have carried out a wide range of tasks – from simple simulation of standard detectors, widely used in experiments, to accurate calculations of secondary CR fluxes, generated by propagation of primary CR particles through Earth’s atmosphere. We obtained the response function for devices of different types, which then worked well in calibration. A detailed study of characteristics of secondary CR fluxes is a separate task, and its description is beyond the scope of this work. It should be noted, however, that a distinctive property of the calculations of propagation of primary CR protons through Earth’s atmosphere is the integration of spectra, obtained in [Vashenyuk et al., 2011], into the model. This allows us to quickly and accurately acquire energy characteristics and vertical profiles of both neutron and charged components (protons, muons, electrons, positrons) not only for ordinary GCR, but also for different instants of SCR events. These results are used to calculate the rate of ionization of Earth’s atmosphere under different initial conditions. In addition, we have turned from the simple concept of plane geometry to a more complex global model, which enables us to use in modeling full information on the primary source of particles and to draw maps for different heights. To implement this model fully, it is supplemented with methods to increase the rate of computation, which along with the simulation results will be described comprehensively in subsequent studies.
A great advantage is the fact that the results are to be compared with experimental data, for which purpose there is a huge database accumulated over years. In conclusion, it is worth noting that the very first comparison gave good agreement with a small, within ~ 10% from mean, deviation.
Список литературы Модельный комплекс для исследования космических лучей
- Agostinelli S., Allison J., Amako K., et al. Geant4 -a simulation toolkit//Nuclear Instruments and Methods in Physics Research. Section A: Accelerators, Spectrometers, Detectors and Associated Equipment. 2003. V. 506, N 3. P. 250 DOI: 10.1016/S0168-9002(03)01368-8
- Basilevskaya G.A., et al. Cosmic ray induced ion production in the atmosphere//Space Sci. Rev. 2008. V. 137. P. 149-173.
- Bertini H.W. Intranuclear-cascade calculations of the secondary nucleon spectra from nucleon-nucleus inter-actions in the energy range 340 to 2900 MeV and comparison with experiment//Phys. Rev. 1969. V. 188. P. 1711-1730.
- Chadwick M.B., Herman M., Obložinský P., et al. ENDF/B-VII.1 Nuclear Data for science and technology: Cross sections, covariances, fission product yields and decay data//Nuclear Data Sheets. 2011. V. 112, iss. 12. P. 2887-2996 DOI: 10.1016/j.nds.2011.11.002
- Clem J.M., Dorman L.I. Neutron monitor response function//Space Sci. Rev. 2000. V. 93. P. 335-359.
- Heikkinen A., Stepanov N., Wellisch J.P. Bertini intra-nuclear cascade implementation in Geant4//Computing in High Energy and Nuclear Physics. 24-28 March 2003, La Jolla, California. MOMT008.PDF.
- Maurchev E.A., Balabin Yu.V., Vashenyuk E.V., Makhmutov V.S. Simulation of the transport of solar protons through the atmosphere in the 13 December 2006 GLE//Physics of Auroral Phenomena: Proc. XXXIV Annual Seminar. Apatity, 2011. P. 110-113.
- Picone J.M., Hedin A.E., Drob D.P., Aikin A.C. NRLMSISE-00 empirical model of the atmosphere: Statistical comparisons and scientific issues//J. Geophys. Res. 2002. V. 107, N A12. P. 1468 DOI: 10.1029/2002JA009430
- Shea M.A., Smart D.F. Possible evidence for a rigidity dependent release of relativistic protons from the solar corona//Space Sci. Rev. 1982. V. 32. P. 251-271.
- Vashenyuk E.V., Balabin Yu.V., Gvozdevsky B.B. Relativistic solar cosmic ray dynamics in large ground level events//Proc. 21st ECRS, Kosice, Slovakia, 9-12 September 2008. Inst. of Exp. Phys Slovak Academy of Sci., 2009. P. 264-268.
- Vashenyuk E.V. Balabin Yu.V., Gvozdevsky B.B. Features of relativistic solar proton spectra derived from ground level enhancement events (GLE) modeling//Astrophys. Space Sci. Trans. 2011. V. 7. P. 459-463.