кандидат физико-математических наук;
candidate of physical and mathematical sciences;
В настоящее время существует значительное разнообразие методик для исследования характеристик потоков как первичных, так и вторичных космических лучей (КЛ). Основными являются экспериментальные методы, с использованием разных детекторов, однако, благодаря современным компьютерным технологиям, экспериментальные методы все чаще совмещаются с численными. Данный шаг позволяет получить наибольшее количество информации о интересующем процессе (явлении) и сделать, соответственно, наиболее точные выводы. В работе представлена концепция нового программного комплекса RUSCOSMIC, основанного на пакете GEANT4 и представляющего собой набор различных численных моделей для исследования прохождения КЛ через вещество различных систем (детекторы излучения, атмосфера Земли). Получены и приведены функции откликов основных детекторов излучения, а также некоторые типовые характеристики потоков вторичных КЛ. Также представлены результаты, демонстрирующие работу модуля верификации вычислений экспериментальными данными.
At present, cosmic ray (CR) physics uses a considerable variety of methods for studying CR characteristics, and both primary and secondary fluxes, respectively. Experimental methods make the main contribution, using various types of detectors, but numerical methods increasingly complement it due to the active development in computer technology. This approach provides researchers with the most extensive information about details of the process or phenomenon and allows us to make the most competent conclusions. This paper presents a concept of the RUSCOSMIC software package based on the GEANT4 toolkit and representing a range of different numerical models for studying the CR passage through a matter of different systems. The obtained results represent response functions of the main radiation detectors as well as some typical characteristics of secondary CR fluxes. Comparative results also show the operation of the module verifica-tion of calculations with experimental data.
ВВЕДЕНИЕВ настоящее время исследования космических лучей (КЛ) проводятся разными способами. Первоначально это были в основном экспериментальные методы, широко использующиеся до сих пор. Традиционно запускаются шары-зонды с установленными на них детекторами [Bazilevskaya et al., 2008], а международная сеть нейтронных мониторов является одной из самых современных и надежных наземных систем регистрации космических лучей. По данным этой сети определяются характеристики релятивистских солнечных протонов, зарегистрированных во время событий солнечных космических лучей (СКЛ), посредством моделирования возрастаний приземного фона вторичных космических лучей (события GLE [Vashenyuk et al., 2009; Shea, Smart, 1982]).С начала компьютеризации активно развивается методика изучения распространения КЛ в веществе с помощью численного метода Монте-Карло. Нами был разработан программный комплекс RUSCOSMIC, включающий в себя модели детекторов разных типов, а также модель прохождения первичных КЛ через атмосферу Земли для получения информации о каскадах вторичных КЛ. Основой является пакет GEANT4 [Agostinelli et al., 2003], из которого наследуются классы, отвечающие за отображение процесса взаимодействия частиц с веществом, построение геометрии, задание начальных параметров, сбор информации о ходе моделирования, состоянии частиц и др.Первым этапом данной работы являлось детальное изучение свойств детектирующего оборудования на станции космических лучей в г. Апатиты. С этой целью были созданы программные модули, представляющие собой модели нейтронного монитора (НМ), счетчиков Гейгера, сцинтилляционных детекторов с геометрией и свойствами материалов, максимально приближенным к реальным.При моделировании прохождения потока нейтронов через нейтронный монитор была получена функция отклика, а также детально изучены свойства замедлителя и свинцового генератора. Сравнение наших результатов с результатами работ, проводившихся ранее и ведущихся в настоящее время, показало хорошее согласие. Современные сечения взаимодействий, используемые при моделировании, позволили более детально изучить функцию отклика нейтронного монитора в диапазоне энергий от 10 МэВ до эпитермальных. Проведенное параллельно моделирование множественности в нейтронном мониторе подтвердило предположение, что феномен вызывается не единственной частицей, попадающей в свинцовый генератор НМ, а целым облаком частиц.Моделирование сцинтилляционных детекторов гамма-квантов стало необходимым, поскольку на станции КЛ уже в течение нескольких лет функционирует система мониторинга рентгеновского излучения, благодаря которой были выявлены его возрастания во время осадков. С помощью GEANT4 вычислены функции отклика NaI-детекторов различной геометрии. Подобные работы, разумеется, проводились и ранее, однако наше исследование отличают современный подход и совершенно новая концепция, позволяющая учесть больше факторов при расчете прохождения частиц через вещество.Вторым этапом было создание собственного программного модуля для расчета прохождения частиц различных энергий через атмосферу Земли. Это позволило не только использовать для моделирования постоянно обновляющиеся данные, но и расширить границы применения, поскольку без труда можно менять состояние самой модели (изменять физические свойства вещества, добавлять различные поля и т. д.). Далее проводились расчеты прохождения галактических космических лучей (ГКЛ) через атмосферу Земли и развития каскадов частиц во время некоторых событий GLE. В результате создана база данных энергетических спектров вторичного космического излучения для различных высот. Данные расчеты могут применяться как для оценки скорости ионизации при вычислении эквивалентной дозы, так и для детального исследования развития каскадов во время GLE и поиска их новых особенностей. Результаты моделирования сверялись с результатами подобных исследований, проводившихся с использованием пакета PLANETOCOSMICS [Maurchev et al., 2011], а также с полученными с помощью шаров-зондов экспериментальными данными, с которыми наблюдается хорошее согласие.