Космические лучи предельно высоких энергий тема автореферата и диссертации по физике, 01.04.23 ВАК РФ

Шустова, Ольга Павловна АВТОР
кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Москва МЕСТО ЗАЩИТЫ
2012 ГОД ЗАЩИТЫ
   
01.04.23 КОД ВАК РФ
Диссертация по физике на тему «Космические лучи предельно высоких энергий»
 
Автореферат диссертации на тему "Космические лучи предельно высоких энергий"

На правах рукописи

Шустова Ольга Павловна

Космические лучи предельно высоких энергий: интерпретация экспериментальных данных и перспективы дальнейших исследований

01.04.23 - физика высоких энергии

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических паук

■ 3 Mi.il 2012

Москва 2012

005019040

005019040

Работа выполнена на кафедре физики космоса физического факультета Московского государственного университета имени М. В. Ломоносова.

Научный руководитель:

доктор физ.-мат. наук, профессор Калмыков Николай Николаевич

Официальные оппоненты:

Роганова Татьяна Михайловна, доктор физ.-мат. наук, Научно-исследовательский институт ядерной физики имени Д. В. Скобельцына МГУ имени М. В. Ломоносова (г. Москва), зав. лабораторией теории электронно-фотонных ливней отдела излучений и вычислительных методов

Рубцов Григорий Игоревич, кандидат физ.-мат. наук, Учреждение Российской академии наук Институт ядерных исследований РАН (г. Москва), научный сотрудник отдела теоретической физики

Ведущая организация:

Государственное учреждение Институт земного магнетизма, ионосферы и распространения радиоволн имени Н. В. Пушкова Российской академии наук (г. Троицк)

Защита состоится ШСьМ, 2012 года в часов ¿¿'минут на заседании диссертационного совета Д501.001.77, Московский государственный университет имени М.В. Ломоносова, 119234, г. Москва, ГСП-1, Ленинские горы, д. 1, стр. 5 (19-й корпус НИИЯФ МГУ), ауд. 2-15.

С диссертацией можно ознакомиться в библиотеке Научно-исследовательского института ядерной физики имени Д. В. Скобельцына МГУ имени М. В. Ломоносова.

Автореферат разослан «оЩ» года.

Ученый секретарь диссертационного сове доктор физ.-мат. наук, профессор

Общая характеристика работы

Актуальность темы

Космические лучи (КЛ) предельно высоких энергий (ПВЭ) привлекают особое внимание начиная с 1966 г., когда Грейзен, Зацепин и Кузьмин предсказали обрезание спектра протонов от далеких источников при энергиях выше ~ 6 • 1019 эВ (ГЗК-порог) в результате взаимодействия с реликтовыми фотонами.

Поток КЛ с энергиями > 1019 эВ составляет примерно одну частицу на квадратный километр в год. Регистрация столь редких событий прямыми методами невозможна. Однако земная атмосфера является прекрасным инструментом для косвенного измерения характеристик КЛ ПВЭ. При взаимодействии первичной частицы КЛ с ядрами элементов, входящих в состав атмосферы, рождается каскад вторичных частиц, так называемый широкий атмосферный ливень (ШАЛ), и первоначальная энергия Ео распределяется между ними. Именно вторичные частицы ШАЛ, а также оптические излучения, создаваемые заряженными частицами ливня, и используются для регистрации КЛ ПВЭ, при этом создаются гигантские установки с площадью в сотни-тысячи квадратных километров. Процедура восстановления типа и энергии первичной частицы по экспериментальным данным требует знания процессов взаимодействия адронов в области предельно высоких энергий. На сегодняшний день таких сведений нет: спектр КЛПВЭ протягивается до ~ 3-102°эВ, что примерно на три порядка выше эквивалентной энергии, которая получена на самом крупном ускорителе - Большом Адронном Коллайдере. Поэтому современные модели адронных взаимодействий опираются на данные, известные при более низких энергиях, и результаты, предсказанные в рамках той или иной модели, могут существенно различаться.

До сих пор остается под вопросом происхождение КЛ ПВЭ. Ограничения, накладываемые на мощность источников в рамках известных моделей ускорения, позволяют с уверенностью говорить о том, что это

внегалактические объекты. Но прямая идентификация источников по направлениям прихода частиц возможна только в случае протонов и легких ядер (при условии, что напряженность магнитных полей в межгалактической среде < Ю-9 Гс). Действительно, ларморовский радиус Ri частицы с зарядом Z для протонов с энергией Ю20эВ, движущихся в галактическом магнитном поле с типичной напряженностью ~ 2 цГс, составляет ~ 50 кпк, что примерно в два раза больше размера Галактики. Однако для ядер железа Ri ~ 2 кпк, поэтому они сильно отклоняются от первоначального направления. Необходимо также учитывать, что при распространении в межгалактической среде частицы КЛПВЭ теряют свою энергию во взаимодействиях с фоновыми излучениями (главным образом, это процессы фоторасщепления ядер и фоторождения пионов протонами). Это приводит к искажению как первоначального состава, так и энергетического спектра частиц, ускоренных в источниках. Таким образом, для того чтобы получить представление о возможных источниках КЛПВЭ, необходимо проводить комплексные расчеты, с учетом процессов взаимодействия частиц KJ1 в межгалактической среде и их отклонений в магнитных полях Галактики.

К настоящему времени в физике КЛПВЭ достигнуты определенные успехи: построены гигантские установки High Resolution Fly's Eye (HiRes), Pierre Auger Observatory (Auger) и Telescope Array (ТА), разработаны комплексные методы для оценки энергии и массового состава первичных частиц, считается подтвержденным обрезание спектра при энергиях выше ГЗК-порога. Тем не менее между экспериментальными данными существуют расхождения.

Несмотря на огромные площади установок для регистрации КЛ ПВЭ, понадобятся еще многие годы для того, чтобы набрать значимую статистику событий с энергиями выше ГЗК-порога. Поэтому в настоящее время готовятся новые эксперименты ТУС, КЛПВЭ и JEM-EUSO для изучения КЛ ПВЭ путем регистрации флуоресцентного излучения ШАЛ.

Преимущество таких детекторов заключается в возможности исследования огромных участков земной атмосферы по всей небесной сфере. Предполагается, что установки космического базирования существенно превзойдут по эффективности современные установки Auger и ТА.

Цель работы

Целью диссертационной работы является анализ и интерпретация современных экспериментальных данных о космических лучах предельно высоких энергий, полученных гигантскими установками, а также исследование возможности регистрации черепковского света ШАЛ детектором космического базирования.

Основные результаты работы

1. Впервые реализован эффективный численный алгоритм решения задачи о распространении в межгалактической среде ядер KJIПВЭ вплоть до железа в диапазоне энергий от 5-101ЭэВ до 2-1О20эВ, позволяющий учитывать все существенные каналы реакций фоторасщепления и обеспечивать в ходе расчета сохранение полного числа нуклонов с точностью не хуже Ю-4 вплоть до максимально исследуемых расстояний (400 Мпк).

2. В рамках модели ускорения частиц, учитывающей их энергетические потери в источнике, получены спектры и массовый состав KJI ПВЭ у Земли для различных показателей степенных спектров инжекции в источниках. Показано, что противоречивые экспериментальные данные установок High Resolution Fly's Eye и Pierre Auger Observatory можно интерпретировать как обусловленные различием в показателях инжекционных спектров. При этом можно удовлетворительно воспроизвести как спектр, так и массовый состав, наблюдаемые на обеих установках.

3. В рамках «жесткостной» модели ускорения проведен анализ вкла-

дов от источников, находящихся на различных расстояниях, в поток КЛПВЭ у Земли. В результате сравнения с экспериментальными данными показано, что интервал расстояний до наиболее вероятных источников сужается до нескольких десятков мегапарсек.

4. В рамках «жесткостной» модели ускорения, на примере близкого источника - аналога радиогалактики Сеп А, подтверждена возможность существования одного или нескольких источников тяжелых ядер, наблюдаемых в Южном полушарии.

5. Впервые получены оценки интенсивности различных составляющих черенковского света, возникающего при развитии ШАЛ над поверхностью океана. С учетом волновой поверхности океана наложены ограничения на зенитные углы ливней, при которых возможна регистрация вспышки отраженных черепковских фотонов детектором космического базирования.

Научная новизна работы

В настоящее время интерес к проблемам энергетического спектра и массового состава КЛПВЭ чрезвычайно возрос. В данной работе был предложен оригинальный подход, состоящий в применении эффективного алгоритма решения задачи о распространении частиц КЛ ПВЭ в межгалактической среде и ориентирующийся на область энергий выше ГЗК-порога. Сравнительный анализ полученных результатов и данных гигантских установок HiR.es и Auger позволил сделать интересный вывод о том, что противоречивые экспериментальные данные можно объяснить различием режимов ускорения в источниках, которые вносят основной вклад в потоки КЛПВЭ, регистрируемые в разных полушариях. Исследование возможности регистрации вспышки черенковских фотонов, отраженных от поверхности океана, установкой космического базирования проведено впервые.

Практическая значимость результатов

Исследования КЛ ПВЭ за ГЗК-порогом проводятся в условиях острого дефицита информации относительно их возможных источников и характера ускорения. Поэтому последовательное рассмотрение достаточно реалистичной модели представляется полезным для решения проблемы происхождения КЛПВЭ. Развитие вычислительных методов, позволяющих эффективно и быстро верифицировать справедливость тех или иных гипотез, также является необходимым этапом продвижения к окончательному решению проблемы. Наконец, результаты, относящиеся к исследованию условий возникновения вспышки черепковских фотонов при отражении от поверхности океана, представляют интерес для развития космических проектов по регистрации оптических излучений ШАЛ.

Вклад автора

Автором созданы компьютерные программы расчета распространения протонов и тяжелых ядер в межгалактической среде с учетом их взаимодействий с фоновыми излучениями, генерации и распространения черепковских фотонов в воздухе и воде с учетом отражения на границе раздела в рамках двух моделей волновой поверхности океана и проведены необходимые вычисления. Все публикации подготовлены автором.

Достоверность результатов

Достоверность результатов работы обусловлена использованием современных данных о плотности инфракрасного излучения и сечений, полученных с помощью хорошо апробированных программ ТАЬУЭ и Сеап14, а также тщательным контролем устойчивости созданных численных схем и обеспечением высокой точности расчетов

Апробация работы

Основные результаты работы докладывались на научных семинарах НИИЯФ МГУ и Ломоносовских чтениях 2010 и 2011 гг., а также на конференциях:

• 12th International Conference on Topics in Astroparticle and Underground Physics, 5-Q September 2011, Munich, Germany;

• 15th International Lomonosov Conference on Elementary Particle Physics, 18-24 August 2011, Moscow;

• International Conference on Advanced Technology and Particle Physics, 7-8 October 2010, Como, Italy;

• SI1"1 Всероссийская конференция но космическим лучам, 5-9 июля 2010 г., Москва;

• XVII Международная конференция студентов, аспирантов и молодых ученых по фундаментальным наукам «Ломоносов-2010», 9-13 апреля 2010 г., Москва.

Публикации

По материалам диссертации опубликовано 5 печатных работ, в том числе 4 статьи в рецензируемых российских и зарубежных изданиях:

1. О. П. Шустова, Н. Н. Калмыков, Б. А. Хренов «Моделирование распространения черенковского излучения от ШАЛ над поверхностью океана» // Вестник Московского Университета. Серия 3. Физика. Астрономия, 2011, № 5, С. 63;

2. О.P. Shustova, N.N. Kalmykov, В.A, Khrenov «On the possibility of registering UHE EAS Cherenkov light by the TUS detector» // Astroparticle, Particle, Space Physics, Radiation Interaction, Detectors and Medical Physics Applications, 2011, V. 6, P. 290;

3. О. П. Шустова, H. H. Калмыков, Б. А. Хренов «Возможность регистрации детектором космического базирования черенковского света от ШАЛ ультравысокой энергии, проникающего в воду океана» // Известия РАН. Серия Физическая, 2011, Т. 75, С. 410;

4. О. П. Шустова, Н. Н. Калмыков, А. В. Урысон «Влияние фонового излучения на состав космических лучей предельно высоких энергий при их распространении от источника до установки» // Известия РАН. Серия Физическая, 2011, Т. 75, С. 342;

N.N. Kalmykov, О. P. Shustova, A.V. Uryson «Spectra and mass composition of ultrahigh-energy cosmic rays from point sources» // [arXiv:1112.5523];

5. О. П. Шустова «Возможность регистрации черенковского света от широких атмосферных ливней от частиц предельно высоких энергий детектором космического базирования» // Сборник тезисов XVII Международной конференции «Ломоносов-2010» Секция «Физика», 2010, Т. 1, С. 38.

Структура и объем диссертации

Диссертация изложена на 98 страницах и включает в себя введение, три главы, заключение, два приложения и список литературы из 124 наименований. Работа содержит 36 рисунков и 5 таблиц.

Содержание работы

Во введении дана общая характеристика диссертационной работы: обоснована актуальность темы, сформулирована цель работы, раскрыты научная новизна работы и практическая значимость результатов, перечислены этапы апробации, указана структура диссертации.

В первой главе представлен обзор существующих сведений о KJI ПВЭ. Экспериментальные данные, полученные гигантскими наземными установками HiRes и Auger, носят противоречивый характер. Во-первых, хотя резкое ослабление потока KJI в области энергий выше ГЗК-порога подтверждается измерениями обеих установок, налицо существенное различие в форме энергетических спектров KJ1 ПВЭ. Данные HiRes хорошо аппроксимируются степенными функциями с показателями ~2.8 и ~5.4

до и после излома при энергии ~ 5.6 • 1019эВ [1]. В то же время результаты Auger показывают, что ослабление потока начинается уже при энергии ~2.9'1019эВ, а спектр хорошо описывается сглаженной функцией [2]. Еще одно противоречие касается данных о массовом составе К Л вблизи ГЗК-порога. Согласно измерениям HiR.es [3], в составе К Л с энергиями выше 1.6 ■ 1018 эВ преобладают протоны, тогда как данные установки Auger [4] указывают на переход от протонов к тяжелым ядрам. Существуют расхождения и в распределениях направлений прихода событий. В данных Auger при энергиях > 5.5 • 1019 эВ [5] найдена значимая корреляция с активными ядрами галактик (АЯГ) из каталога [6] внутри сферы с радиусом ~ 75 Мпк на угловых масштабах 3.1°. Кроме того, обнаружен значимый избыток событий в направлении на радиогалактику СепА. Напротив, данные HiRes хорошо согласуются с гипотезой об изотропном распределении [7].

Одним из возможных объяснений указанных расхождений является существование одного или нескольких ближних источников тяжелых ядер, которые вносят значительный вклад в поток КЛ ПВЭ, регистрируемый установкой Auger в южном полушарии [8]. В пользу этой гипотезы говорят недавно опубликованные и согласующиеся с результатами HiRes данные установки ТА [9], которая, как и HiRes, расположена на территории штата Юта, США. Тем не менее современная статистика событий КЛПВЭ не позволяет сделать окончательных выводов. Однако создание наземных установок с площадью большей, чем у Auger, уже невыгодно. На сегодняшний день решением являются космические эксперименты по регистрации флуоресцентного излучения ШАЛ (КЛПВЭ [10], JEM-EUSO [11]), которые существенно увеличат экспозицию за разумный интервал времени наблюдения.

Вторая глава посвящена исследованию, целью которого явилась попытка объяснить различия между данными HiRes и Auger относительно формы энергетического спектра и массового состава, а также определить

возможные режимы ускорения в источниках KJIПВЭ. Саму задачу можно сформулировать следующим образом: при заданных параметрах режима ускорения в источниках, их расстояниях и интеисивностях рассчитать энергетический спектр и массовый состав KJ1 у Земли в интервале энергий от 5-101ЭэВ до 2-102ОэВ и путем сравнения с экспериментальными данными HiR.es и Auger определить оптимальные параметры режима ускорения, которые позволяют описать результаты обеих установок.

Основные модельные предположения заключаются в следующем. В качестве источников рассматриваются «ближние» астрофизические объекты - ядра сейфертовских галактик в радиусе ~ 40 Мпк из каталога [6], причем их интенсивность в KJI, состав и условия ускорения считаются одинаковыми. Такой выбор мотивирован работами [12,13], где данные объекты были отождествлены как возможные источники КЛПВЭ. Кроме того, для них была предложена модель ускорения частиц с учетом энергетических потерь на синхротронное излучение [14,15]. В расчетах максимальные энергии Етах, до которых ускоряются частицы в источнике, пропорциональны полученным в рамках этой модели значениям, а коэффициент пропорциональности выбирается из условия наилучшего согласия расчетных спектров с экспериментальными данными. Состав частиц в источнике соответствует содержанию элементов в галактической среде, а инжекционные спектры описываются обратным степенным законом с показателями 7о и 7о + 4 при энергиях ниже и выше Етах. Поскольку в случае степенного спектра количество частиц с массовым числом А и энергиями из диапазона (Е, E+dE) пропорционально Л70-1 [16], доля тяжелых ядер с заданной энергией в источнике существенно возрастает при достаточно больших значениях 7о-

Энергетический спектр и массовый состав частиц, ускоренных в источниках, искажаются в результате взаимодействия с фоновыми излучениями и космическими магнитными полями. Поэтому при интерпретации экспериментальных данных необходимо детальное моделирование

прохождения частиц в межгалактической среде. В интересующем диапазоне энергий важны процессы фоторасщепления ядер и фоторождения пионов протонами при взаимодействиях с реликтовым и инфракрасным фоновым излучениями. Однако для оценки вклада ядер в поток вторичных протонов учитывается также фоторождение пионов ядрами. Следует отметить, что предположение о ближних источниках КЛПВЭ позволяет пренебречь процессом образования электрон-позитронных пар протонами и ядрами. Кроме того, в расчетах не учитывается ослабление потока КЛ ПВЭ за счет отклонений частиц в магнитных полях.

Приблизительное сохранение Лоренц-фактора в учитываемых процессах позволяет свести задачу о распространении ядер, ускоренных в источнике, к решению цепочки дифференциальных уравнений для вторичных ядер. В целом рассматривается свыше 100 изотопов. Уравнение, описывающее эволюцию ядра с зарядом 2 и массовым числом А, имеет вид:

-НЦх) \\%2) + £ £ А^+0+

к,1 т

"и*) +

п=0,1

где Аг - обратная средняя длина свободного пробега для г-го канала реакции (фоторасщепление с испусканием одного или нескольких нуклонов, а-чаегиц и т.д., /3-распад и фоторождение пиона), а А'04 - обратная полная длина свободного пробега. Программа расчета позволяет учесть все значимые каналы реакций фоторасщепления.

При вычислении спектров ядра объединяются в группы в зависимости от их массового числа. Это группа гелия Не, а также группы легких Ь, средних М, тяжелых Н и сверхтяжелых УН ядер.

Результаты расчета показывают (см. рис. 1), что в рамках принятых модельных предположений можно получить хорошее согласие с аппроксимациями спектров обеих установок, причем единственным парамет-

т

йх

ром является показатель инжекционного спектра: 7о = 2.2 (HiRes) и 4.3 (Auger).

Вследствие недостаточной статистики событий КЛ ПВЭ массовый состав измерен только до ГЗК-порога. Поэтому при сравнении предполагается, что резкого изменения состава в области более высоких энергий не происходит. На рис. 3 представлен средний логарифм массового числа частиц, соответствующего спектру на рис.1. Здесь же указаны значения (In Л), извлеченные из данных о средней глубине максимума ливней (Хтах) при самой высокой энергии Ес (~ 5 • 1019эВ для HiRes [3], ~3.5 • 1019эВ для Auger [4|) путем интерполяции между предсказаниями моделей адронных взаимодействий EPOSvl.99 и QGSJET01 для «протонного» и «железного» составов. Данные модели дают экстремальные предсказания. Вообще, интерпретация массового состава сильно зависит от выбранной модели взаимодействия. Тем не менее очевидно, что вблизи ГЗК-порога массовый состав, измеренный HiRes, является более легким, поэтому результаты расчета качественно согласуются с данными обеих установок.

Энергетические спектры и массовый состав КЛПВЭ от источников, равномерно распределенных в слое от 1 до 40 Мпк, также согласуются с данными обеих установок (см. рис. 2, 3). Это свидетельствует о том, что полученные результаты не слишком чувствительны к распределению ближних источников.

В данной главе также исследуется вопрос о расстояниях до возможных источников КЛПВЭ. При этом рассматривается равномерное распределение источников, а максимальные энергии ядер определяются в рамках «жесткостной» модели ускорения. Расчеты средних расстояний, с которых группы ядер могут достигнуть нашей Галактики, показывают, что источники КЛ ПВЭ должны находиться в пределах нескольких десятков мегапарсек. Сравнение спектров от источников, расположенных в пределах 40, 100 и 400 Мпк, с экспериментальными данными также

19.8 20.0 lg(E/eV)

Рис. 1 — Дифференциальные спектры групп ядер от ближних источников при 7о = 2.2 (слева) и 4.3 (справа), р (Не) - поток вторичных протонов, образованных ядрами гелия. Полный спектр, представленный серыми кружками, нормирован на единицу при 5-1019эВ. Экспериментальные точки установок HiRes и Auger, а также их аппроксимации (штриховые линии) показаны на графиках слева и справа, соответственно.

10

19.6 19.8 20.0 20.2 IgCE/eV)

о Auger

у -4.3

10"

10

19.6 19.8 20.0 20.2 lg(E/eV)

■Р РОВД

* Не

□ L

+ М

i Н

Л VH

Ф tot

Рис. 2 — Дифференциальные спектры групп ядер от источников, равномерно распределенных в слое от 1 Мпк до 40 Мпк, при 70 = 2.2 (слева) и 4.3 (справа). См. пояснение к рис. 1.

19.7 19.9 20.1 lg(E/eV)

19.7 19.9 20.1 lg(E/eV)

Рис. 3 — Энергетическая зависимость среднего логарифма массового числа частиц, соответствующего спектрам на рис. 1 (светло-серые треугольники) и 2 (темно-серые кружки). Для сравнения представлены результаты, извлеченные из данных HiRes (слева) и Auger (справа) о величине (Хтах) при энергии Ес; точки, определенные из предсказаний моделей EPOSvl.99 (звездочки) и QGSJET01 (квадратики), смещены в противоположные стороны относительно Ес на Д lg£' = 0.01.

демонстрирует, что ослабление потока от далеких источников наступает слишком быстро. Следует подчеркнуть, что использование двух разных моделей ускорения дает качественно одинаковые результаты, а именно, для объяснения данных установок HiRes и Auger требуются различные значения инжекционного индекса.

Сама по себе интерпретация противоречивого характера результатов измерений как обусловленного сменой режимов ускорения кажется несколько странной, поскольку приводит к фундаментальному различию между источниками, находящимися в поле зрения каждой установки. Более естественным является предположение о том, что один или несколько внегалактических объектов в поле зрения установки Auger являются источниками тяжелых ядер и дают существенный вклад в регистрируемый поток KJI [8). В данной главе эта гипотеза рассматривается по отношению к мощной радиогалактике Сеп А, которая широко обсуждается как возможный источник KJIПВЭ (см., например, [17,18]). Спектр, полученный при условии, что источник на расстоянии ~4 Мпк (аналог Сеп А) с инжекционным индексом 70 = 4.6, с одной стороны, и остальные АЯГ

из каталога [б) в радиусе 40 Мпк с 70 = 2.2, с другой, дают приблизительно одинаковые вклады в поток KJIПВЭ, хорошо согласуется с гладкой аппроксимацией Auger.

В третьей главе представлены оценки различных компонент черепковского света, возникающих при развитии ливня над поверхностью океана и направленных в верхнюю полусферу, и результаты исследования условий, при которых окажется возможной регистрация детектором ТУС вспышки черенковских фотонов ШАЛ, отраженных от волновой поверхности океана.

В ближайшее время планируется осуществление космического эксперимента ТУС [10,19], нацеленного на исследование КЛПВЭ путем регистрации флуоресцентного излучения, возникающего при прохождении ШАЛ. Это излучение изотропно, что позволяет регистрировать фотоны, идущие под произвольными углами к оси ливня. Одно из преимуществ такого эксперимента состоит в возможности проводить поиск событий КЛ ПВЭ по всей небесной сфере одним детектором.

Помимо флуоресценции ШАЛ сопровождается мощным потоком че-ренковского излучения, которое концентрируется в основном вдоль оси ливня. Тем не менее при прохождении ШАЛ над поверхностью с высоким альбедо отраженные черепковские фотоны можно регистрировать на больших высотах. Наилучшие условия для такой регистрации возникают при отражении от снега (нормальное альбедо чистого снега достигает ~90 %). Еще в 1971 г. Чудаков указал на возможность регистрировать ШАЛ по черенковскому свету, отраженному от заснеженной поверхности [20). Данный метод положен в основу успешно развивающегося эксперимента СФЕРА-2 [21].

Однако для установок космического базирования особый интерес представляет возможность регистрации черенковского света, отраженного от водной поверхности, поскольку ~70% земного шара покрыто водой. Хотя нормальное альбедо воды составляет всего ~2%, поток черенковских

фотонов ШАЛ примерно в пять раз больше флуоресцентного потока, поэтому данная проблема является весьма актуальной в преддверии запуска детектора ТУС. Регистрация вспышки черепковских фотонов в конце флуоресцентного трека существенно повысила бы надежность экспериментальных данных, позволив уточнить траекторию ливня.

В потоке черепковского света, возникающего при развитии ШАЛ над водной поверхностью, можно условно выделить несколько компонент, направленных в верхнюю полусферу. Часть черенковских фотонов идет вверх в результате рассеяния в атмосфере, хотя значительно большая доля излучения достигает поверхности океана. При этом часть фотонов отражается и также идет в верхнюю полусферу, а остальные фотоны входят в воду. Однако вследствие рассеяния в воде некоторая доля фотонов выходит на поверхность. Кроме того, остаток ливня, имеющий энергию порядка десятков процентов от первичной, при попадании в воду генерирует большое число черенковских фотонов, часть из которых также выходит на поверхность воды. Для будущих космических проектов расчет указанных компонент может оказаться полезным. Например, предполагается, что рассеянные черепковские фотоны будут создавать существенный «шум», величину которого необходимо знать для выделения флуоресцентного сигнала.

Генерация фотонов и их распространение в среде моделировались методом Монте-Карло. В табл. 1 при различных значениях длины волны фотона А и зенитного угла ливня в¡¡^ представлены: И^ея - доля фотонов, отраженных от поверхности воды; И'ьаск _ доля рассеянных в воздухе фотонов, которые достигли высоты 400 км; У/ц,-™ ~ доля фотонов, попавших в воду и вышедших в результате рассеяния на поверхность;

~ доля фотонов, образованных в воде остатком ливня и вышедших на поверхность раздела в результате рассеяния. Расчеты для Wr<,f^, И^-«,, И^тй выполнены в модели абсолютно гладкой границы раздела сред.

При наблюдениях на больших высотах важную роль играет угловое

А, нм Л^/Ю13 нм-1 % И'Ьаск % % Л7"7ю12 нм-1 %

300 5.44 1.6 1.4 0.02 26.1 0.2

0° 340 4.24 1.6 1.3 0.03 20.4 0.3

390 3.22 1.6 1.0 0.05 15.5 0.4

300 6.43 1.7 1.8 0.02 12.6 0.3

30° 340 5.00 1.7 1.6 0.04 9.83 0.4

390 3.81 1.7 1.3 0.06 7.47 0.4

300 7.11 1.9 2.5 0.02 2.81 0.5

45° 340 5.54 2.1 2.4 0.04 2.20 0.6

390 4.21 2.2 1.9 0.06 1.67 0.7

300 7.20 2.5 4.0 0.02 0.02 0.9

60° 340 5.60 3.3 4.1 0.03 0.01 1.0

390 4.26 3.9 3.4 0.06 0.01 1.2

Таблица 1 - Вклады компонент черенковского излучения, возникающих при развитии ливня над поверхностью океана и направленных в верхнюю полусферу (см. пояснение в тексте). Также представлено число черенковских фотонов

образованных в воздухе (воде), по отношению к которому получены доли IV««, 1УЬасЪ И/а_„ (1У„а1).

распределение фотонов, поэтому необходим учет волн на поверхности океана. В настоящей работе рассматриваются две модели волновой поверхности океана: 1) периодической структуры и 2) хаотической структуры. В первой модели поверхность представляется в виде неподвижной волны с длиной Аш и амплитудой вдоль оси Ох, форма которой описывается параметрическими уравнениями укороченной циклоиды. Во второй модели поверхность разбивается на множество мелких плоских ячеек высотой г). Величины 1]и = дт]/дх, г\с = дг)1ду (индексы «и» и «с» указывают на направление по ветру и поперек) следуют двумерному нормальному распределению с дисперсиями а\ = 3.16-10-3и„,, а% = 1.92-10~3Иг„, где иш - скорость ветра в м/с [22]. В обеих моделях предполагается, что на границе раздела фотоны отражаются и преломляются согласно законам геометрической оптики. Принятые значения параметров волновой поверхности океана воспроизводят усредненную картину, наблюдаемую в спокойную погоду, что должно благоприятствовать регистрации ШАЛ

О 4 8 12 16 20 24 28 Л, км

0 4 8 12 16 20 24 28 Я, км

Рис. 4 - Число отраженных черенковских фотонов, которые попадают на зеркало детектора в зависимости от величины Л, для ливней с зенитными углами в,и = 0° (а) и 20° (б). Каждое кольцо разбивается на 4 равных сегмента, в зависимости от азимутального угла фотонов: =45° —135°, <¿>2 = 135° — 225°, <¿>3=225° —315°, <Р4 = 315°—45°.

из космоса.

Основным параметром, определяющим возможность регистрации вспышки отраженных черенковских фотонов, является зенитный угол ливня. Кроме того, важно учитывать расположение детектора относительно оси ливня. На рис. 4 представлено число фотонов, которые попадают на зеркало детектора ТУС в зависимости от величины Я в случае ливней с зенитными углами 0^ = 0° и 20° в модели периодической структуры (Д соответствует радиусу круга на высоте 400 км с центром в проекции точки пересечения осью ливня гладкой поверхности океана). Как видно, наблюдается также зависимость от азимутального угла фотона <р, поскольку в модели периодической структуры угловое распределение отраженных фотонов имеет «выбросы» поперек фронта волны. Учитывая временное разрешение детектора ТУС (~1 мкс), можно считать, что все отраженные черепковские фотоны будут регистрироваться одновременно.

В рамках модели периодической структуры волновой поверхности океана с параметрами волны hw = 0.7 м, \w = 40 м показано, что регистрация вспышки отраженных фотонов в конце флуоресцентного трека возможна в ливнях с зенитными углами < 20°. Анализ, выполненный в модели хаотической структуры при скорости ветра uw = 6 м/с, дает ограничение 6sh<25°.

В заключении перечислены основные результаты диссертационной работы.

Цитируемая литература

[1] P. Sokolsky (The High Resolution Fly's Eye Collaboration) «Observation of the GZK cutoff by the HiRes Experiment» // Nuclear Physics В Proceedings Supplements, 2009, V. 196, P. 67.

[2j J. Abraham et al. (The Pierre Auger Collaboration) «Measurement of the energy spectrum of cosmic rays above 1018eV using the Pierre Auger Observatory» // Physics Letters B, 2010, V. 685, P. 239 [arXiv: 1002.1975].

[3] R. U. Abbasi et al. (The HiRes Collaboration) «Indications of protondominated cosmic-ray composition above 1.6 EeV» // Physical Review Letters, 2010, V. 104, P. 161101 [arXiv:0910.4184],

[4] J. Abraham et al. (The Pierre Auger Collaboration) «Measurement of the depth of maximum of extensive air showers above 1018eV» // Physical Review Letters, 2010, V. 104, P. 091101 [arXiv:1002.0699].

[5] P. Abreu et al. (The Pierre Auger Collaboration) «Update on the correlation of the highest energy cosmic rays with nearby extragalactic matter» // Astroparticle Physics, 2010, V. 34, P. 314 [arXiv: 1009.1855].

[6] M.-P. Veron-Cetty, P. Veron «А catalogue of quasars and active nuclei: 12th edition» // Astronomy and Astrophysics, 2006, V. 455, P. 773.

18

[7] R. U. Abbasi et al. (The High Resolution Fly's Eye Collaboration) «Search for correlations between HiRes stereo events and active galactic nuclei» // Astroparticle Physics, 2008, V. 30, P. 175 [arXiv:0804.0382],

[8] B. Schwarzschild «The highest-energy cosmic rays may be iron nuclei» // Physics Today, 2010, May, P. 15.

[9] Y. Tsunesada (The Telescope Array Collaboration) «Highlights from Telescope Array» // Proceedings of the 32nd International Cosmic Ray Conference, August 2011, Beijing, China [arXiv:l 111.2507].

[10] B. A. Khrenov et al. «KLYPVE/TUS space experiments for study of ultrahigh-energy cosmic rays» // Physics of Atomic Nuclei, 2004, V. 67, P. 2058.

[11] Y. Takahashi (The JEM-EUSO Collaboration) «The JEM-EUSO mission» // New Journal of Physics, 2009, V. 11, P. 065009 [arXiv:0910.4187].

[12] A. V. Uryson «Possible sources of ultrahigh-energy cosmic protons» // Journal of Experimental and Theoretical Physics Letters, 1996, V. 64, P. 77.

[13] A. V. Uryson «Identification of sources of ultrahigh energy cosmic rays» // Astronomical Reports, 2001, V. 45, P. 591.

[14] A.V. Uryson «Ultrahigh energy cosmic ray acceleration in Seyfert galactic nuclei» // Astronomy Letters, 2001, V. 27, P. 775.

[15] A.V. Uryson «Seyfert nuclei as sources of ultrahigh-energy cosmic rays» // Astronomical Reports, 2004, V. 48, P. 81.

[16] V. L. Ginzburg, S. I. Syrovatskii «The origin of cosmic rays» // Pergamon Press, Oxford, 1964.

[17] D.S. Gorbunov, P.G. Tinyakov, I.I. Tkachev, S.V. Troitsky «On the interpretation of the cosmic-ray anisotropy at ultra-high energies» // [arXiv:0804.1088].

[18] T. Piran «A new limit on the distances of nuclei UHECRs sources» // [arXiv:1005.3311].

[19] V. Abrashkin et al. «Space detector TUS for extreme energy cosmic ray study» // Nuclear Physics В Proceedings Supplements, 2007, V. 166, P. 68.

[20] A. E. Чудаков «Возможный метод регистрации ШАЛ по черенков-скому излучению, отраженному от заснеженной поверхности Земли» // Материалы Всесоюзного симпозиума «Экспериментальные методы исследования космических лучей сверхвысоких энергий», 1972, Якутск, С. 69.

[21] R. А. Antonov et al. «Antarctic balloon-borne detector of high-energy cosmic rays (SPHERE project)» // Radiation Physics and Chemistry, 2006, V. 75, P. 887.

[22] C. D. Mobley «Light and water: radiative transfer in natural waters» // London, Academic Press, 1994.

Подписано в печать: 20.04.12

Объем: 1,5усл.п.л. Тираж: 100 экз. Заказ № 761 Отпечатано в типографии «Реглет» 119526, г. Москва, ул. Бауманская д.ЗЗ (495) 979-96-99; www.reglet.ru

 
Текст научной работы диссертации и автореферата по физике, кандидата физико-математических наук, Шустова, Ольга Павловна, Москва

61 12-1/773

МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ имени М. В. ЛОМОНОСОВА

ФИЗИЧЕСКИЙ ФАКУЛЬТЕТ

На правах рукописи

Шустова Ольга Павловна

Космические лучи предельно высоких энергий: интерпретация экспериментальных данных и перспективы

дальнейших исследований

01.04.23 — физика высоких энергий

ДИССЕРТАЦИЯ на соискание ученой степени кандидата физико-математических наук

НАУЧНЫЙ РУКОВОДИТЕЛЬ доктор физ.-мат. наук, профессор Н. Н. Калмыков

Москва - 2012

Содержание

Введение 4

1 Космические лучи предельно высоких энергий 9

1.1 Экспериментальные установки............................................11

1.1.1 High Resolution Fly's Eye..........................................12

1.1.2 Pierre Auger Observatory..........................................14

1.1.3 Telescope Array......................................................16

1.1.4 Космические проекты ТУС и JEM-EUSO......................17

1.2 Энергетический спектр....................................................18

1.3 Массовый состав . ......................................22

1.4 Анизотропия................................................................24

2 Энергетический спектр и массовый состав KJI ПВЭ от ближних источников 29

2.1 Распространение KJI ПВЭ в межгалактической среде..................30

2.1.1 Взаимодействия протонов ........................................30

2.1.2 Взаимодействия ядер..............................................32

2.2 Постановка задачи и модельные предположения ......................35

2.3 Решение задачи ............................................................38

2.3.1 Ядра..................................................................39

2.3.2 Протоны ............................................................41

2.3.3 Спектры протонов и групп ядер..................................43

2.4 Спектры и состав КЛ ПВЭ от одиночного источника..................44

2.5 Сравнение результатов расчета с данными HiRes и Auger............48

2.6 Расстояния до возможных источников KJI ПВЭ........................53

2.7 Сеп А как возможный источник тяжелых ядер KJ1 ПВЭ..............56

3 Регистрация черенковского излучения ШАЛ от KJI ПВЭ детектором космического базирования 58

3.1 Постановка задачи и модельные предположения ......................59

3.2 Моделирование генерации и распространения черенковских фотонов 61 3.2.1 Число черенковских фотонов в ливне............................62

3.2.2 Схема генерации фотона..........................................63

3.2.3 Схема распространения фотона в среде..........................66

3.2.4 Доли компонент черенковского излучения ШАЛ..............70

3.3 Волновая поверхность океана ............................................71

3.3.1 Модель периодической структуры................................72

3.3.2 Модель хаотической структуры..................................73

3.3.3 Параметры моделей................................................73

3.4 Ограничения на зенитные углы ливней..................................74

Заключение 82

Приложение 1 84

Приложение 2 85

Список литературы 86

Введение

Актуальность темы

Космические лучи (КЛ) предельно высоких энергий (ПВЭ) привлекают особое внимание начиная с 1966 г., когда Грейзен [1], Зацепин и Кузьмин [2] предсказали обрезание спектра протонов от далеких источников при энергиях выше ~ 6-1019 эВ (ГЗК-порог) в результате взаимодействия с реликтовыми фотонами.

Поток КЛ с энергиями > 1019 эВ составляет примерно одну частицу на квадратный километр в год. Регистрация столь редких событий прямыми методами невозможна. Однако земная атмосфера является прекрасным инструментом для косвенного измерения характеристик КЛ ПВЭ. При взаимодействии первичной частицы КЛ с ядрами элементов, входящих в состав атмосферы, рождается каскад вторичных частиц, так называемый широкий атмосферный ливень (ШАЛ), и первоначальная энергия Е0 распределяется между ними. При Еа > 1019 эВ значительное число ливневых частиц достигает земной поверхности, и для их регистрации используют гигантские установки с площадью в сотни-тысячи квадратных километров, что само по себе является довольно сложной, с технической точки зрения, задачей. Кроме того, процедура восстановления типа и энергии первичной частицы по экспериментальным данным требует знания процессов взаимодействия ад-ронов в области предельно высоких энергий. На сегодняшний день таких сведений нет: спектр КЛПВЭ протягивается до ~ 3 • Ю20 эВ, что примерно на три порядка выше эквивалентной энергии, которая получена на самом крупном ускорителе — Большом Адронном Коллайдере (ЬНС). Поэтому современные модели адрон-ных взаимодействий опираются на данные, известные при более низких энергиях, и результаты, предсказанные в рамках той или иной модели, могут существенно различаться.

До сих пор остается под вопросом происхождение КЛПВЭ. Ограничения, накладываемые на мощность источников в рамках известных моделей ускорения, позволяют с уверенностью говорить о том, что это внегалактические объекты. Однако прямая идентификация источников по направлениям прихода частиц возможна только в случае протонов и легких ядер (при условии, что напряженность магнитных полей в межгалактической среде < 10~9 Гс). Действительно, ларморов-

ский радиус частицы с зарядом Z и энергией Е в магнитном поле напряженностью В равен

1 Е 2 уТс

и для протонов с энергией 102ОэВ, движущихся в галактическом магнитном поле с типичной напряженностью ~ 2 /.¿Гс, составляет ~ 50 кпк, что примерно в два раза больше размера Галактики. Однако для ядер железа с той же энергией Rl ~ 2 кпк, поэтому они сильно отклоняются от первоначального направления. Необходимо также учитывать, что при распространении в межгалактической среде частицы KJIПВЭ теряют свою энергию во взаимодействиях с фоновыми излучениями (главным образом, это процессы фоторасщепления ядер и фоторождения пионов протонами). Это приводит к искажению как первоначального состава, так и энергетического спектра частиц, ускоренных в источниках. Таким образом, для того чтобы получить представление о возможных источниках КЛПВЭ, необходимо проводить комплексные расчеты, с учетом процессов взаимодействия частиц КЛ в межгалактической среде и их отклонений в магнитных полях Галактики.

К настоящему времени в физике КЛПВЭ достигнуты определенные успехи. Так, например, построены гигантские установки High Resolution Fly's Eye (HiRes), Pierre Auger Observatory (Auger) и Telescope Array (ТА), разработаны комплексные методы для оценки энергии и массового состава первичных частиц, считается подтвержденным обрезание спектра при энергиях выше ГЗК—порога. Тем не менее между экспериментальными данными существует расхождения.

Несмотря на огромные площади наземных установок для регистрации КЛ ПВЭ, понадобятся еще многие годы для того, чтобы набрать значимую статистику событий с энергиями выше ГЗК-порога. Поэтому в настоящий момент готовятся космические эксперименты ТУС, КЛПВЭ и JEM-EUSO для изучения КЛПВЭ путем регистрации флуоресцентного излучения ШАЛ. Преимущество таких детекторов заключается в возможности исследования огромных участков земной атмосферы по всей небесной сфере. Предполагается, что установки космического базирования существенно превзойдут по эффективности современные установки Auger и ТА.

Цель работы

Целью диссертационной работы является анализ и интерпретация современных экспериментальных данных о космических лучах предельно высоких энергий, полученных гигантскими установками, а также исследование возможности регистрации черепковского света ШАЛ детектором космического базирования.

Научная новизна работы

В настоящее время интерес к проблемам энергетического спектра и массового состава КЛ ПВЭ чрезвычайно возрос. В данной работе был предложен оригинальный подход, состоящий в применении эффективного алгоритма решения задачи о распространении частиц КЛ ПВЭ в межгалактической среде и ориентирующийся на область энергий выше ГЗК-порога. Сравнительный анализ полученных результатов и данных гигантских установок HiRes и Auger позволил сделать интересный вывод о том, что противоречивые экспериментальные данные можно объяснить различием режимов ускорения в источниках, которые вносят основной вклад в потоки К Л ПВЭ, регистрируемые в разных полушариях. Исследование возможности регистрации вспышки черенковских фотонов, отраженных от поверхности океана, установкой космического базирования проведено впервые.

Практическая значимость результатов

Исследования КЛ ПВЭ за ГЗК-порогом проводятся в условиях острого дефицита информации относительно их возможных источников и характера ускорения. Поэтому последовательное рассмотрение достаточно реалистичной модели представляется полезным для решения проблемы происхождения К Л ПВЭ. Развитие вычислительных методов, позволяющих эффективно и быстро верифицировать справедливость тех или иных гипотез, также является необходимым этапом продвижения к окончательному решению проблемы. Наконец, результаты, относящиеся к исследованию условий возникновения вспышки черенковских фотонов при отражении от поверхности океана, представляют интерес для развития космических проектов по регистрации оптических излучений ШАЛ.

Апробация работы

Основные результаты работы докладывались на научных семинарах НИИЯФ МГУ и Ломоносовских чтениях 2010 и 2011гг., а также на конференциях:

• 12th International Conference on Topics in Astroparticle and Underground Physics, 5-9 September 2011, Munich, Germany;

• 15th International Lomonosov Conference on Elementary Particle Physics, 18-24 August 2011, Moscow;

• International Conference on Advanced Technology and Particle Physics, 7-8 October 2010, Como, Italy;

• 31ая Всероссийская конференция по космическим лучам, 5-9 июля 2010 г., Москва;

• XVII Международная конференция студентов, аспирантов и молодых ученых по фундаментальным наукам «Ломоносов-2010», 9-13 апреля 2010 г., Москва.

По материалам диссертации опубликовано 5 печатных работ:

1) О. П. Шустова, Н. Н. Калмыков, Б. А. Хренов «Моделирование распространения черенковского излучения от ШАЛ над поверхностью океана» // Вестник Московского Университета. Серия 3. Физика. Астрономия, 2011, № 5, С. 63 [arXiv:1110.4870];

2) О. P. Shustova, N. N. Kalmykov, В. A. Khrenov «On the possibility of registering UHE EAS Cherenkov light by the TUS detector» // Astroparticle, Particle, Space Physics, Radiation Interaction, Detectors and Medical Physics Applications, 2011, V. 6, P. 290;

3) О. П. Шустова, H. H. Калмыков, Б. А. Хренов «Возможность регистрации детектором космического базирования черенковского света от ШАЛ ультравысокой энергии, проникающего в воду океана» // Известия РАН. Серия Физическая, 2011, Т. 75, С. 410 [arXiv:1110.2974];

4) О.П. Шустова, H.H. Калмыков, A.B. Урысон «Влияние фонового излучения на состав космических лучей предельно высоких энергий при их распространении от источника до установки» // Известия РАН. Серия Физическая, 2011, Т. 75, С. 342;

N. N. Kalmykov, О. P. Shustova, А. V. Uryson «Spectra and mass composition of ultrahigh-energy cosmic rays from point sources» // [arXiv:1112.5523];

5) О. П. Шустова «Возможность регистрации черенковского света от широких атмосферных ливней от частиц предельно высоких энергий детектором космического базирования» // Сборник тезисов XVII Международной конференции «Ломоносов-2010» Секция «Физика», 2010, Т. 1, С. 38.

Структура и объем диссертации

Диссертация изложена на 98 страницах и включает в себя введение, три главы, заключение, два приложения и список литературы из 124 наименований. Работа содержит 36 рисунков и 5 таблиц.

1 Космические лучи предельно высоких энергий

На Землю из космоса приходит различные типы излучений. В данной работе под термином космическщ луча будут пониматься потоки заряженных частиц от протонов до ядер железа.

ч

Fluxes of Cosmic Ravs

__(I particle per л""'-second)

«пег

(1 particle per r^-year)

\ /

10

\

10

-25 1 10

Ankle

(1 particle per iarr-veor) S,

V

Energy {eV)

Рис. 1. Компиляция измерений дифференциального Э1к>ргеги-теского спектра космических лучей. 1 Ь'нктирпом линией представлен степенной закон Е 3 it указаны энергий, соответствующие приблизительным интегральным потокам 1 частица/(м^сек), 1 частица/(м2год) л 1 частица/(к м^грд). (Из работы [3| )

Энергетический спектр КЛ простирается от lfsB до ~ Ю20 эВ и описывается степенной функцией вида Е"в (рисЛ). В области колена (knee), при энергии Е — 1015-5эВ. показатель в возрастает от до 3.0. Еще один излом наблюдается при Е щ. 10 эВ: спектр выполажйвается до Е~2Л, формируя лодыжку (ankle). Начиная с энергий > 5. Ю13 зВ поток КЛ становится настолько малым,

10

-"С

что современные экспериментальные данные не позволяют с уверенностью судить о поведении спектра в данной области. Можно, однако, утверждать, что предсказанное в 1966 г. Грейзеном |1|, Зацепиным и Кузьминым [2] обрезание спектра КЛ ПВЭ, действительно, существует. В рамтах проекта ШНеэ было убедительно показано резкое ослабление потока КЛ выше ГЗК-порога ¡4]. Ныне действующие установки Ан<>,ег и ТА также подтверждают этот результат [5], |6].

Изучению основных характеристик К Л ПВЭ (энергетический спектр, массовый состав, распределение направлении прихода) сопутствует ряд ограничений и трудностей. Во-первых, малость потока КЛ с анергиями > 1019 эВ 1 событие на квадратный километр в год) означает неприменимость прямых методов регистрации. Возможность получать ипформапшо о КЛ косвенным методом обусловлена существованием земной атмосферы, в которой первичная частица, взаимодействуй с ядрами элементов, входящих в состав воздуха, развивает каскад из большого числа вторичных частиц — широкий атмосферный ливень. Размеры ливней от частиц с энергиям» > Ю1оэВ настолько велики, что часть вторичных частиц доходит до поверхности и может быть зарегистрирована наземными детекторами. Для получения статистически значимых результатов в случае ливней от частиц предельно высоких энергий необходимо строить огромные по площади установки и набирать статистику в течение многих лет. Кроме того, определение энергетического спектра и массового состава КЛ путем регистрации вторичных частиц ШАЛ включает в себя задачу о восстановлении параметров первичной частицы по откликам детекторов, входящих и состав установки. Подобная процедура подразумевает использование моделей адронных взаимодействий частиц при энергиях, которые недоступны в лабораторных условиях. Так. при взаимодействии протона с энергией 1019 эВ и покоящегося нуклона в атмосфере их энергия в системе центра, масс достигает сотен ТэВ, тогда как максимальная энергия современных ускорителей составляет около 7 ТэВ. Поскольку все используемые модели адронных взаимодействий являются чисто феноменологическими, то, строго говоря, надежность их предсказаний вне пределов области энергий, внутри которой были определены параметры модели, нельзя гарантировать. Более того, пи одна из этих моделей не позволяет описать все характеристики, измеренные на ЬНС в р—р столкновениях при эквивалентной энергии 3-ЮшэВ [7].

На сегодняшний день яспользуются различные методы Исследования ШАЛ. Исторически первым был предложен метод регистрации заряженных частиц. С помощью детекторов электронов и мюонов можно оценить энергию, перешедшую, соответственно. в электронно-фотонную и мюониую компоненты ливня. Сравнение полученных данных с усредненными характеристиками ливня на определенной глубине и с заданным зенитным углом позволяет восстановить параметры первичной частицы. Кроме того, широкое распространение получили методы изучения ШАЛ. основанные на регистрации сопутствующих ливню оптических излучений. Важное преимущество данных методов по сравнению с регистрацией заряженных частиц ШАЛ состоит в том, что оценки величины потоков черепковского и флуоресцентного излучений в меньшей степени подвержены модельной зависимости. Однако их реализация возможна лишь в ясные безлунные ночи, что сокращает рабочее время до 5 — 10%.

В настоящей главе дается краткий обзор существующих сведений о КЛ ПВЭ. В разделе 1.1 приводится описание гигантских наземных установок для регистрации КЛ ПВЭ, а также будущих космических проектов. В разделах 1.2-1,4 обсуждаются основные характеристики КЛ —■ энергетический спектр, массовый состав и анизотропия — в области энергий вблизи и выше ГЗК—порога и имеющиеся противоречия между экспериментальными данными гигантских установок.

1.1 Экспериментальные установки

По величине эффективной площади (или апертуры) установки, регистрирующие ШАЛ с первичными энергиями > 1О10эВ, условно можно разделить на несколько поколений. К первому поколению относятся наземные решетки Volcano Rancli [8]. Haverah Park [9]. Sidney University Giant Air-shower Recorder (SUGAR.) jlOj, Akeno Giant Air Shower Array (AGASA) [11|> Якутская комплексная установка ШАЛ |12j и флуоресцентный телескоп Fly's Eye |13], со сравнительно небольшой эффективной площадью ß дееятки-сотни квадратных километров. Следующим этапом и развитии установок для регистрации ШАЛ стал проект High Resolution Fly's Eye [14| с эффективной апертурой ~ 10 ООО км2стер при энергии частиц ~ 10 эВ. Третье поколение включает ныне действующие гибридные установки Pierre Anger Observatory [15] и Telescope Array [IG), объединяющие наземные pe-

шеткп и флуоресцентные телескопы. В настоящее время также идет подготовка детекторов нового поколения: в рамках космических экспериментов ТУС (17(. КЛГ1ВЭ |18) и .ТЕМ-ЕШО 119| планируется изучение КЛПВЭ путем регистрации флуоресцентного излучения ШАЛ с борта орбитального спутника Земли. Эффективная площадь детекторов КЛПВЭ и ЛЕМ-ЕИЗО будет достигать десятков тысяч квадратных километров. Еще одним преимуществом окажется возможность обзора всей небесной сферы одним детектором. На рис. 2 представлена эволюция экспозиции ряда установок для регистраци�