Элементные спектры космических лучей: два галактических источника и внегалактический переход за коленом ✦
Энергетические спектры отдельных элементов в космических лучах в диапазоне от 100×Z ГэВ до 1000×Z ПэВ демонстрируют два излома-«колена» с жесткостной зависимостью: при энергиях около 15×Z ТэВ и 3×Z ПэВ, а также три «лодыжки» — выполаживания. Первая лодыжка (≈500×Z ГэВ) также зависит от жесткости, а вот вторая (≈150 ТэВ для протонов, ~1 ПэВ для гелия) и, вероятно, третья (≈100 ПэВ для протонов) — нет. Благодаря высокой точности измерений в экспериментах DAMPE и LHAASO, а также IceTop, впервые удалось установить, что все спектральные особенности хорошо описываются суммой трёх компонент: низкоэнергетической галактической (с выпуклой формой спектра, затухающей после ТэВ-колена), высокоэнергетической галактической (включая ПэВ-колено) и внегалактической. Никакой дополнительной галактической компоненты в интервале 10 ПэВ – 1 ЭэВ не требуется.
Контекст
Природа колена в спектре галактических космических лучей при ~3 ПэВ остаётся одной из центральных загадок астрофизики частиц. Исторически предполагалось, что оно может быть связано с максимальной энергией ускорителей, эффектами распространения или фотоядерными процессами. Ещё в 1990-е годы появились указания на то, что спектр формируется несколькими популяциями источников. Классическая идея Фрица Цвикки (1934 г.) о том, что остатки сверхновых являются главными ускорителями, остаётся доминирующей, однако точная конфигурация источников и их вклад в разные энергетические диапазоны были неясны до появления прецизионных измерений нового поколения.
Методы
В работе использованы данные по спектрам протонов, гелия и железа с экспериментов DAMPE (прямые измерения на спутнике), LHAASO и IceTop (наземные регистрации широких атмосферных ливней). Благодаря спектроскопии высокого разрешения удалось выделить индивидуальные элементы. Функция фита представляла сумму трёх компонент (две галактические и одна внегалактическая). Для описания спектра одиночного источника предложена новая феноменологическая форма «двойная лог-парабола с экспоненциальным обрезанием» (EDLP), которая воспроизводит наблюдаемое уплощение перед обрывом. Популяция источников моделировалась свёрткой спектра одиночного источника с функцией светимости, после чего параметры подгонялись под данные методами MINUIT (MIGRAD) и MCMC (emcee). Дополнительно учитывались потери на фрагментацию для тяжёлых ядер с использованием модельной зависимости пробега от жёсткости.
Результаты
Показано, что спектры водорода, гелия и железа в диапазоне 40 ГэВ – 500 ПэВ удовлетворительно описываются трёхкомпонентной моделью. Первая галактическая компонента имеет выпуклую форму с низкоэнергичным показателем наклона 2.76 (одинаковым для всех ядер) и заканчивается после ТеВ-колена при жёсткости ~15×Z ТэВ. Вторая компонента доминирует вплоть до ПэВ-колена (3×Z ПэВ) и ответственна за главный излом. Третья, внегалактическая, компонента появляется при ~100 ПэВ для протонов и становится доминирующей выше 1 ЭэВ. Обнаружено, что так называемые «лодыжки» — изломы в сторону уплочения при 500×Z ГэВ и ~150 ТэВ (протоны) или ~1 ПэВ (гелий) — являются следствием пересечения разных компонент и не всегда подчиняются жёсткостной зависимости из-за различий в химическом составе популяций. Формальный фит спектра гелия требует высокоэнергичного обрезания первой компоненты на 700 ТэВ, однако эта деталь не является обязательной при учёте систематических погрешностей. Важно, что для описания данных в области 10 ПэВ – 1 ЭэВ не требуется введения дополнительного третьего галактического населения.
Значение
Предложенная интерпретация означает, что переход от галактических к внегалактическим космическим лучам происходит «бесшовно» сразу за ПэВ-коленом, а не при традиционно рассматриваемой «лодыжке» в полном спектре при 4 ЭэВ. Это поддерживает модель «разочаровывающего» состава (disappointing model) для ультравысоких энергий и подчёркивает, что Галактика перестаёт быть эффективным ускорителем протонов уже при энергиях ~100 ПэВ. Выпуклая форма первого компонента может указывать на нелинейные эффекты в ускорительных процессах (диффузионное ударное ускорение) либо на вклад локальных источников на поздних стадиях эволюции сверхновых.
Развитие
Будущие прецизионные измерения спектров ядер группы CNO (углерод, кислород) и железа в области колена на установках LHAASO и планируемой SWGO позволят независимо проверить предсказания трёхкомпонентной модели и жёсткостную зависимость изломов. Совместный анализ анизотропии и спектров в диапазоне 0.1–10 ПэВ, особенно с учётом фазы диполя, поможет выделить вклад близких источников. Численное моделирование ускорения в остатках сверхновых с реалистичными нелинейными обратными связями и переменным составом окружения станет важным шагом к теоретическому обоснованию наблюдаемой картины.
Влияние
Результаты окажут влияние на физику космических лучей, гамма-астрономию высоких энергий (поиск галактических PeVatron'ов) и интерпретацию данных нейтринных телескопов, поскольку резкий переход к внегалактическим протонам меняет ожидания от диффузного нейтринного фона.
Следующие шаги
Проверка предсказанной зависимости «бегущего» перехода от заряда ядра на новых данных по CNO и железу с LHAASO; уточнение сечений фрагментации тяжёлых ядер в области ТэВ; моделирование спектров вторичных ядер (бор, бериллий) в рамках трёхкомпонентной модели.
Ключевые проблемы
Предложенная модель напрямую связывает происхождение колена с исчерпанием возможностей галактических ускорителей и немедленным выходом на сцену внегалактических источников, что требует пересмотра традиционных представлений о «галактическом» диапазоне вплоть до 1 ЭэВ. Кроме того, выпуклая форма первого компонента и её возможная связь с нелинейным ускорением на фронтах ударных волн в остатках сверхновых остаётся предметом активных дебатов, соединяя проблему колена с тонкостями физики плазмы.
🎯 Протон с энергией 1 ЭэВ в типичном галактическом магнитном поле (3 мкГс) имеет ларморовский радиус около 300 парсек — почти половину толщины диска, поэтому Галактика не может удерживать такие частицы, и они почти наверняка приходят из-за её пределов.
Ключевые числа
- энергия ТеВ-колена для протонов: 15 ТэВ
- энергия ПэВ-колена для протонов: 3 ПэВ
- энергия первого излома (лодыжки) для протонов: 500 ГэВ
- энергия второго излома для протонов: 150 ТэВ
- низкоэнергичный спектральный индекс первой компоненты: 2.76