АСТРОНОМИЯ, АСТРОФИЗИКА И КОСМОЛОГИЯ
Диагностика области ядерных взаимодействий нейтронов в процессе солнечной вспышки 20 января 2005 г. посредством гамма-излучения
в линиях
Е. В. Троицкая1,0, И. В. Архангельская2,6, А. И. Архангельский2
1 Московский государственный университет имени М. В. Ломоносова, Научно-исследовательский институт ядерной физики имени Д. В. Скобельцына. Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2.
2 НИЯУ МИФИ, кафедра экспериментальной физики и космофизики.
Россия, 115409, Москва, Каширское ш., д. 31.
E-mail: а [email protected], ь [email protected] Статья поступила 07.10.2011, подписана в печать 24.11.2011.
На основе анализа временного профиля скорости счета гамма-излучения в линии 2.223 МэВ, зарегистрированного во время экстремального солнечного события 20 января 2005 г., исследованы плотность окружающей среды и содержание изотопа 3Не в области генерации гамма-линии 2.223 МэВ. К анализу были привлечены данные, полученные на борту спутника КОРОНАС-Ф аппаратурой АВС-Ф с детектора СОНГ-Д о скорости счета в гамма-линии с энергией 2.223 МэВ, а также в линиях в диапазоне 4-7 МэВ. Показано, что для построения согласованной с имеющимися данными модели временного профиля потоков гамма-излучения в линии 2.223 МэВ в рассматриваемом событии необходимо допустить повышенное содержание изотопа 3Не в области вспышки относительно водорода 1Н. В области генерации линии 2.223 МэВ отношение концентраций n(3He)/rt(1H) должно составлять (1.4±0.15)- 10^4, тогда как в среднем по данным от различных вспышек обычно оно составляет около 2-10^5. Повышенное содержание 3Не служит также косвенным подтверждением возможного наличия продуктов реакции радиационного поглощения нейтронов изотопом 3Не в составе обнаруженного ранее по данным КОРОНАС-Ф/АВС-Ф небольшого возрастания над фоном гамма-излучения в диапазоне 15-21 МэВ.
Ключевые слова: солнечная вспышка, гамма-излучение, ядерные реакции, плотность солнечной атмосферы, содержание изотопа 3Не.
УДК: 523.985-866. PACS: 96.60.qe, 95.85.Pw.
Введение
Экстремальное солнечное событие 20 января 2005 г. было зарегистрировано в период аномально высокой вспышечной активности Солнца в январе 2005 г. в фазе глубокого спада солнечного цикла. По рентгеновским данным космического аппарата GOES в диапазоне 0.1-8 нм, вспышка 20 января 2005 г. (класс Х7.1, балл ЗВ, гелиокоординаты N14, W61) началась в 06:36 UT, достигла максимума в 07:01 UT и закончилась в 07:26 UT. Источником вспышки была активная область NOAA 10720. Эта вспышка сравнима по мощности с наиболее сильными из наблюдавшихся ранее. Она сопровождалась корональным выбросом (СМЕ), а также мощным протонным событием и выдающимся GLE (ground level enhancement) [1].
Гамма-излучение во время изучаемой вспышки наблюдалось приборами спутников RHESSI [2], INTEGRAL [3], а также на борту ИСЗ КОРОНАС-Ф [4, 5] детектором СОНГ-Д [6] в двух параллельных экспериментах: собственно прибором СОНГ [7, 8] и аппаратурой АВС-Ф (амплитудно-временная спектрометрия Солнца) [9-13], в которых осуществлена регистрация гамма-излучения в широком энергетическом диапазоне (вплоть до 250 МэВ).
На борту спутника КОРОНАС-Ф аппаратурой
АВС-Ф с детектора СОНГ-Д были получены данные о временном профиле скорости счета фотонов в гамма-линии с энергией 2.223 МэВ, возникающей при радиационном поглощении вспышечных нейтронов водородом солнечной атмосферы в период солнечной вспышки 20 января 2005 г. Наблюдаемый временной профиль совместно с данными о потоках излучения в ядерных линиях в диапазоне 4-7 МэВ [9, 10] были предварительно исследованы в работе [14] с помощью метода, предложенного ранее в НИИЯФ МГУ. При этом были рассчитаны модельные временные профили линии 2.223 МэВ в предположении бесселевского энергетического спектра ускоренных частиц для трех значений спектрального индекса аТ = 0.005, 0.03, 0.1 (рис. 1). Расчеты были проведены в различных предположениях о вертикальном профиле плотности солнечной плазмы, изотропного испускания нейтронов в нижнюю полусферу и обычно принимаемого в расчетах предположения о содержании изотопа 3Не в области взаимодействия нейтронов с окружающей средой: к = и(3Не)/и(1Н) = 2 • 10^5.
В отличие от четырех исследованных ранее солнечных событий, в случае события 20 января 2005 г. моделирование с достаточной точностью оказалось невозможным в рамках обычных предположений о свойствах солнечной атмосферы и параметрах вспышки. Начиная
Одним из видов энерговыделения в процессе солнечной вспышки является ускорение заряженных частиц — электронов, протонов и ионов в короне над структурой магнитных петель. С точки зрения дальнейших ядерных процессов в атмосфере Солнца наибольший интерес представляют энергичные протоны и ионы, которые первоначально направлены в нижнее (по отношению к Солнцу) полупространство. Наиболее значительное 7-излучение в линиях производят в конечном итоге реакции генерации нейтронов, а также реакции возбуждения ядерных уровней ряда наиболее распространенных средних по массе элементов, главным образом, 12 С и 16 О. Наиболее вероятная область первоначальных ядерных реакций соответствует концентрациям 10'1 — 1012 частиц-см~3. Возбужденные ядра |2С* и |60* практически мгновенно, через время ~ 10~9 с, испускают 7-кванты с энергией соответственно 4.44 и 6.13 МэВ. Нейтроны же, обладающие значительным средним временем жизни до распада (886 с) и относительно большими сечениями упругого взаимодействия с водородом (зависящими от энергии), замедляются при упругих столкновениях с ядрами окружающей среды в среднем за времена ^0.05-0.8 с до тепловых скоростей, проникая в глубокие слои хромосферы и фотосферу. Далее нейтроны начинают более эффективно захватываться ядрами водорода окружающей среды с образованием дейтрона и испусканием 7-кванта с энергией 2.223 МэВ. Сечения последнего процесса малы, что приводит к временам захвата от нескольких десятков до сотни секунд.
В НИИЯФ МГУ ранее была предложена и развита методика моделирования временного профиля потока гамма-излучения в линии 2.223 МэВ, возникающего во время солнечных вспышек при захвате нейтронов ядрами водорода окружающей среды [15, 16]. Методом статистического моделирования мы прослеживаем траекторию каждого нейтрона от его генерации до одного из возможных процессов ее окончания: поглощения
Отсч./с
0 200 400 600 800
Время, с
Рис. 1. Сравнение наилучших аппроксимаций с данными наблюдений [8, 9] для временного профиля потока гамма-квантов с энергией 2.223 МэВ от вспышки 20 января 2005 г. при четырех моделях плотности (т) (см. табл. 1), трех значениях спектрального индекса ускоренных частиц аТ и величине отношения к = 2-10~5 [14]. Ноль времени на этом и последующих графиках соответствует 06:43:00.499 иТ. Расчетная кривая для модели плотности т=3 на рисунке не показана, так как в рассматриваемом случае она практически сливается с кривой для т = 2
аТ= 0.005
т = 1
т =2
т = 4 модели т = 5
20.01.2005 данные АВС-Ф/СОНГ-Д
аТ= 0.03
аТ= 0.1
с фазы максимума наблюдаемые скорости счета убывают быстрее со временем, чем моделирующие скорости счета. В настоящей работе мы обсуждаем две наиболее вероятные причины отклонения стандартно рассчитанного временного профиля от наблюдательных данных: повышенное содержание 3 Не в области взаимодействий нейтронов и характер геометрии испускания генерированных во вспышке нейтронов. Для этого мы включаем эти параметры в новые модельные расчеты временных профилей. Новые расчеты основываются на тех же наблюдательных данных. Целью работы является нахождение наиболее вероятного значения таких параметров вспышки, как содержание 3Не в области генерации нейтронов и характер вертикального профиля плотности солнечной атмосферы во время вспышки. Мы также оцениваем спектральный индекс первоначально ускоренных протонов и исследуем влияние характера геометрии испускания генерированных во вспышке нейтронов на временной профиль 7-излучения в линии 2.223 МэВ.
1. Метод
протоном и излучением 7-кванта, безрадиационным поглощением ядром изотопа 3Не, распадом нейтрона или его выходом от Солнца.
Для расчетов траекторий нейтронов необходима информация как о свойствах среды, так и о начальных характеристиках рождающихся в процессе вспышки вторичных нейтронах. Некоторые из необходимых величин при этом могут являться параметрами расчета. Начальный спектр нейтронов мы задаем, основываясь на работах [17, 18], в которых спектр генерированных во вспышке нейтронов рассчитан исходя из спектров ускоренных протонов в предположении бесселевского спектра последних, что характерно для стохастического ускорения. Таким образом, вид временного профиля 7-линии 2.223 МэВ опосредованно связан со спектром первоначальных протонов, что позволяет получить оценку их спектра.
Для задания начального временного профиля генерированных нейтронов используется общепринятый метод, впервые примененный в работе [19]. При этом профиль генерированных нейтронов принимается соответствующим временному профилю суммы гамма-излучения в узких «мгновенных» гамма-линиях. Такая возможность основана на следующих соображениях. Для генерации нейтронной линии 7-излучения основную роль играет участок спектра нейтронов от 10 до 100 МэВ. Характерные черты энергетической зависимости сечений реакций генерации, с одной стороны, нейтронов, а с другой — ядерных 7-линий, возникающих при высвечивании возбужденных уровней ядер, в условиях падающего вида спектра первичных ускоренных протонов приводят к подобию, в основных чертах, спектральных областей их генерации. Наиболее сильными ядерными 7-линиями являются линии от кислорода 16 О и углерода 12 С с энергиями 6.13 и 4.44 МэВ соответственно. Их суммарный временной профиль вместе с другими, более слабыми линиями
в диапазоне от 4 до 7 МэВ обычно и принимается пропорциональным временному профилю генерируемых нейтронов. При расчетах было учтено также поглощение генерированных 7-квантов с энергией 2.223 МэВ в солнечной атмосфере.
Кроме того, рассмотрение отдельных временных интервалов дает возможность по временному профилю 7-излучения находить изменения некоторых параметров вспышки со временем. В работах [20, 21] было дано подробное описание методики и ее применения к изучению свойств солнечной атмосферы и физики вспышек.
Для обеспечения возможности нахождения вида зависимости плотности солнечной плазмы от глубины в атмосфере Солнца в 1996 г. были предложены дополнительно к известной модели плотности спокойной солнечной атмосферы различные дополнительные модели [22]. Модели подобраны таким образом, чтобы были представлены случаи как повышенной, так и пониженной плотности атмосферы [20] (табл. 1). В работах [16, 22] впервые была исследована зависимость временного профиля 7-излучения 2.223 МэВ от характера высотного профиля плотности солнечной атмосферы. Приложение методики к анализу крупных солнечных событий, сопровождавшихся измеренным 7-излучением, 22 марта 1991, 6 ноября 1997, 16 декабря 1988 и 28 октября 2003 г. [14, 20 (и цит. лит.), 21] для нахождения ряда параметров вспышечного процесса привело к следующим важным, с точки зрения настоящей работы, выводам.
1. В период вспышки плотность солнечной атмосферы не описывается моделью спокойной атмосферы. Плотность изменяется более резко, чем в рамках модели спокойного Солнца и достигает более высоких значений либо в фотосфере, либо в подфотосферных слоях. При этом наиболее вероятной моделью оказывается модель с уплотнением в нижней хромосфере
Таблица 1
Характеристики принятых моделей плотности солнечной атмосферы; т — порядковый номер модели. За верхний уровень фотосферы принят уровень, на котором оптическая глубина для длины волны 5000 Ä
составляет т = 0.005
т Характер изменения плотности с глубиной в атмосфере Солнца
Кратко Подробно
1 Комбинированная модель спокойной атмосферы: модель HSRA для нижней хромосферы и фотосферы [23] в комбинации с моделью Spruit [24] для конвективной зоны Плавно возрастает от 1.5- 1016 см^3 на вершине фотосферы до 2 • 1017 см^3 на уровне, лежащем на 330 км ниже, где т= 1, и резко возрастает до т= 10 на протяжении 60 км вглубь
2 Повышена относительно 1 внутри и под фотосферой Увеличена до 8 • 1017 см^3 на глубинах ~500 км под верхним уровнем фотосферы, т. е. в глубоких подфотосферных слоях
3 То же Увеличена до 6 • 1017 см^3 на тех же глубинах
4 Понижена относительно 1 начиная с нижней хромосферы и глубже Понижение плотности начинается еще над фотосферой, на ее вершине плотность достигает 3 • 1015 см^3, а на уровне, лежащем ниже на 330 км, составляет 2 • 1016 см^3
5 Скомбинированная из трех участков модель: в хромосфере — модель HSRA; излом и резкий рост в нижней хромосфере, повышена по всей толще фотосферы; глубже — модель Spruit Повышенные плотности в толще фотосферы — по всей ее толще и несколько выше ее, приблизительно на 300 км и ниже плотность составляет 2 • 1017 см^3
и фотосфере (т = 5, табл. 1). Эта модель представляет собой композицию из начального (верхнего) участка в хромосфере являющегося областью действия известной модели НБИА [23], далее, после излома, участка в области нижней хромосферы и фотосферы с резким ростом плотности с глубиной, выходом на плато и наиболее глубоким участком, описываемым моделью Спруита [24]. Модели в верхнем и нижнем участках основаны на данных по оптическим излучениям, в среднем — на нашем анализе данных по 7-излучению в процессе вспышек, исследованных ранее.
Надо отметить, что впоследствии работу, связанную также с выбором модели атмосферы, позволяющей объяснить гамма-излучение солнечной вспышки, выполнили авторы [25] в 2007 г., которые также исследовали зависимость временного профиля 7-излучения в линии 2.223 МэВ от характера модели плотности солнечной атмосферы. Они пришли к выводу о том, что моделью, наилучшим образом описывающей пока только одну исследованную ими в отношении модели вспышку (4 июня 1991 г.), является, как и в нашем случае, составная модель, имеющая некоторые общие черты с нашей моделью т = 5 по своей конфигурации. Модель авторов [25] скомбинирована из модели АэЬ-\уапс1еп (1996, основана на данных по рентгеновскому излучению), с переходом глубже (в точке излома) в модель АугеМ (1981, основана на оптических данных по области, лежащей над пятном) (см. [25] и цит. лит.).
2. Исследование событий 16 декабря 1988 и 28 октября 2003 г. показало, что наша модель т = 5 с уплотнением в фотосфере начинает реализовываться, лишь спустя некоторое время от начала генерации 7-излучения в линии 2.223 МэВ: в первом случае и 140 с, а во втором и 200 с. На самом первом этапе нарастания интенсивности более вероятной оказывается стандартная модель спокойной атмосферы.
3. Исследование этих же событий (16 декабря 1988 и 28 октября 2003 г.) показало, что спектр ускоренных частиц становится более жестким за время наблюдения 7-излучения в линии 2.223 МэВ.
2. Данные наблюдений
В работах [9, 10] представлен временной профиль скорости счета в энергетическом диапазоне 1.7-2.3 МэВ, преимущественный вклад в который дает 7-линия 2.223 МэВ. Этот профиль был «очищен» от локального фона космического аппарата и фона тормозного 7-излучения релятивистских электронов во время вспышки.
В качестве исходных данных в области энергий, соответствующих ядерным линиям, мы имели скорости счета в двух энергетических диапазонах: в аппаратурный энергетический диапазон 3.2-5.0 МэВ основной вклад вносят 7-кванты с энергией 4.44 МэВ от высвечивания возбужденных протонами ядер 12 С, а в диапазон 5.3-6.9 МэВ — 7-кванты с энергией 6.13 МэВ от высвечивания возбужденных ядер 16О. Для дальнейшей работы с этими ядерными линиями, помимо обеих составляющих фона, было вычтено также излучение той части континуума, которая возникает из-за доплеровского уширения при «обратных» ядерных реакциях. Последний эффект возникает из-за возбуж-
дения ускоренных ионов (главным образом 12 С и 16 О) при их столкновениях с ядрами водорода окружающей среды и последующим излучением с их возбужденных уровней. Процедура «очистки» данных по линиям 4.44 и 6.13 МэВ описана в работах [9, 10].
3. Приложение метода
Для применения рассматриваемой методики к анализу данного события модельные расчеты временного профиля нейтронной линии выполнены со следующими параметрами.
1. Спектральный индекс ускоренных частиц аТ, который, как показано в разд. 2, связан со спектром генерированных нейтронов, которые, в свою очередь, участвуют в формировании временного профиля нейтронной линии. Как и ранее (рис. 1), сделано предположение о том, что спектр ускоренных протонов можно аппроксимировать функцией Бесселя. Рассмотрены случаи с характерными значениями параметра стохастического ускорения аТ = 0.005, 0.03, 0.1.
2. Модель высотной зависимости плотности солнечной атмосферы, пг — номер модели в табл. 1. В новых расчетах использованы два различных вертикальных профиля плотности солнечной плазмы: модель пг = 5 с уплотнением в нижней хромосфере и фотосфере до 2 • 1017 см^3 и модель (пг = 1) спокойной солнечной атмосферы.
3. Содержание 3Не относительно водорода к. С целью исследования возможной роли повышенного содержания изотопа 3Не в процессе данной вспышки выполнены расчеты временных профилей потоков гамма-излучения в линии 2.223 МэВ при условии повышенного содержания 3Не для значений параметра к, с шагом 3■10~5, начиная с 2•10^5 и до 2•10^4.
В качестве временного профиля инжекции нейтронов используется суммарный временной профиль 7-линий, выделенных после процедуры очистки от всех трех составляющих фона излучения в диапазонах 3.2-5.0 МэВ и 5.3-6.9 МэВ [9, 10] (разд. 2). Основные расчеты выполнены для случая изотропного начального испускания нейтронов в нижнее полупространство.
4. Результаты моделирования
В табл. 2 в строках с 1-й по 7-ю представлены суммы квадратов отклонений рассчитанных аппроксимаций от наблюдательных данных. Результаты представлены для двух моделей: модели спокойной атмосферы (пг = 1) и модели в период вспышки (пг = 5) по трем значениям параметра аТ. На рис. 2, 3 приведены расчетные временные профили вместе с наблюдательными данными при двух значениях к.
По всем рассмотренным относительным значениям содержания 3Не суммы квадратов отклонений минимальны в случае использования модели пг = 5. Анализ показывает, что с увеличением содержания 3Не приближение к расчетным данным в случае модели плотности пг = 5 и спектрального индекса аТ = 0.1 улучшается, сумма квадратов проходит через минимум при к = (1.4 ± 0.15) • 10""4, далее возрастает.
Таблица 2
Суммы квадратов отклонений наилучших аппроксимирующих кривых от данных наблюдений временного профиля линии 2.223 МэВ во время вспышки 20 января 2005 г. для двух случаев испускания нейтронов: равномерной инжекции в нижнее полупространство (первые семь строк данных) и веерообразной инжекции в последних трех строках, к = л(3Не)/л(1Н)
аТ = 0.005 аТ = 0.03 а7 = 0.1
тп = 1 тп = 5 тп = 1 тп = 5 тп = 1 тп = 5
к = 2•10^5
1.258- 105 2.880 • 104 8.076 - 104 1.991 - 104 7.322 - 104 1.862- 104
к = 5•10^5
9.419- 104 1.366- 104 5.200 - 104 8.470 - 103 4.702 - 104 7.897 - 103
к = 8•10^5
6.882 • 104 1.159- 104 3.754 - 104 5.910- 103 3.283 - 104 5.264 - 103
к= 1.1 • 10^4
5.803 • 104 6.482 - 103 2.713- 104 4.353 - 103 2.292 - 104 4.110- 103
к= 1.4- 10^4
3.946 • 104 5.035 - 103 1.974- 104 3.857 - 103 1.695- 104 3.779 - 103
к= 1.7- 10^4
3.636 • 104 6.687 - 103 1.536- 104 5.103- 103 1.286- 104 5.007 - 103
к = 2.0 • 10^4
2.827- 104 7.329- 103 1.172- 104 6.201 - 103 9.747 - 103 6.155- 103
к = 2 • 10^5, ф = 89°
1.881 ■ 105 4.376 - 104 1.421 - 105 3.249 - 104 1.331 - 105 3.056 - 104
к= 1.1 • 10^4, ф = 89°
9.833 • 104 1.504- 104 6.074 - 104 9.054 - 103 5.483 - 104 8.286- 103
к= 1.4 • 10^4, ф = 89°
7.477 • 104 1.366- 104 4.532 - 104 7.482 - 103 4.070 - 104 6.860 - 103
5. Обсуждение результатов и выводы
Отметим, что на начальной стадии гамма-излучения вспышки 20 января 2005 г. (фаза роста) наилучшая аппроксимация может достигаться при меньших значениях содержания 3Не (рис. 1 ,а). Динамика роста содержания 3Не со временем может быть объяснена постепенным его накоплением в фотосфере. При этом можно предположить, что увеличенное содержание 3Не в области взаимодействия нейтронов со средой может быть обусловлено преимущественным ускорением ионов 3Не в короне. Например, в работе [26] предложен механизм ионно-звуковой турбулентности в качестве ответственного за подобный процесс. Затем ионы распространяются с потерей энергии в солнечной атмосфере вниз, в нижнюю хромосферу и фотосферу, где и могут накапливаться.
К тому же анализ тонкой структуры временных профилей исследуемой вспышки в диапазонах, которые соответствуют выделенным в работах [9, 10] особенностям в энергетическом спектре, позволяет предполагать возможность многократного ускорения частиц в течение вспышки. Действительно, в суммарных энергетических спектрах этой вспышки выделены четыре
особенности в диапазоне 0.1-8 МэВ, соответствующие позитронной и ядерным линиям и линии от захвата нейтронов, которые наблюдались на всем протяжении регистрации гамма-излучения. При этом на временных профилях 7-излучения, построенных в диапазонах, соответствующих выделенным линиям, присутствует тонкая структура с характерными масштабами 33-92 с на уровне значимости 99% (см. также [3]).
С другой стороны, если исходить из наилучшей аппроксимации, найденной по полному времени излучения нейтронной линии (рис. 3, в, см. также рис. 2, в), то предварительный вывод можно сделать о характере плотности солнечной атмосферы на начальной стадии вспышки. Наиболее вероятной моделью плотности в это время, в отличие от фазы спада, оказывается стандартная модель спокойной атмосферы Солнца (т= 1).
Для продвижения в вопросе о динамике как содержания 3Не, так и плотности плазмы, необходимы более детальные исследования, связанные с приложением метода к отдельным временным интервалам вспышечного 7-излучения в линии 2.223 МэВ.
Наряду с повышенным содержанием изотопа 3Не существуют и другие возможности, которые могли бы
Отсч./с О. Отсч./с Я
100
100
100
данные
• 20.01.2005 данные АВС-Ф/СОНГ-Д
100
10
т = 5
данные
• 20.01.2005 данные АВС-Ф/СОНГ-Д
0 200 400 600 800
Время, с
Рис. 2. Наблюдаемые временные профили скорости счета от 7-линии 2.223 МэВ во время вспышки 20 января 2005 г. [8, 9] и результаты моделирования со спектральными индексами аТ = 0.005, 0.03 и 0.1 для двух моделей плотности при к = я(3Не)/я('Н) = 8 • 10^5
0 200 400 600 800
Время, с
Рис. 3. Наблюдаемые временные профили скорости счета от 7-линии 2.223 МэВ во время вспышки 20 января 2005 г. [8, 9] и результаты моделирования со спектральными индексами аТ = 0.005, 0.03 и 0.1 для двух моделей плотности при к = rt(3He)/«(1H) = 1.4- 10-4
рассматриваться в качестве причины, воздействующей на временной профиль 7-линии 2.223 МэВ от вспышки, в частности, наиболее реально можно ожидать некоторого изменения временного профиля потока 7-из-лучения в рассматриваемой линии при отклонении распределения первоначальных направлений нейтронов в нижнюю полусферу от изотропного.
В связи с этим обстоятельством мы рассмотрели случай веерообразного углового распределения первоначальных нейтронов, которое может реализоваться для определенной доли нейтронов и при распространении ускоренных первичных частиц в магнитной вспышечной петле [27]. Действительно, при распространении ускоренных протонов и других ионов в магнитном поле вспышечной петли частицы в своем ларморовском движении отражаются от оснований арки со сходящимся магнитным полем, где в силу большей плотности вещества солнечной атмосферы становятся более вероятными ядерные реакции генерации нейтронов. Поскольку образовавшиеся нейтроны в основном сохраняют направления породивших их частиц, то направления нейтронов распределены равномерно в плоскости сечения петли. Отметим, что эта доля зависит как от спектра протонов, так и от условий в магнитной петле.
Данные табл. 2 (строки 8-10) демонстрируют, что и в этом случае минимум суммы квадратов отклонений достигается при повышенном содержании 3Не, хотя соответствие с наблюдательными данными становится хуже как в предположении обычного, ранее наблюдавшегося в других вспышках содержания изотопа 3Не, так и в случае повышенного.
Поскольку прямое варьирование большого количества параметров затруднено в рамках настоящего метода, связанного с решением обратной задачи, то окончательное решение вопроса об увеличенном содержании изотопа 3 Не могло бы базироваться на дополнительных данных, полученных в изучаемой вспышке.
В частности, в энергетическом спектре этой вспышки обнаружено слабое возрастание скорости счета 7-излучения в диапазоне 15-21 МэВ [10, 11], которое представляет собой превышение над фоном на ^ 2.5а в период 06:44:52-06:51:16 ОТ и на Зст в период максимума (06:47:00-06:49:08 11Т), где а — стандартная ошибка. Хотя среднеквадратичное отклонение Зст находится на минимальном уровне, условно допустимом для утверждения о статистически значимом отклонении, тем не менее, более подробный анализ показывает, что эффект наблюдается в шести независимых спектрах, идущих последовательно по времени, начиная с 06:42:43 ОТ 25 января 2005 г. Тренд увеличивает статистическую значимость наблюдаемого возрастания.
В работе [28] показано, что эта спектральная особенность может быть вызвана комбинацией двух 7-ли-ний с энергиями 15.11 МэВ и 20.58 МэВ. Одна из них испускается возбужденным ядром 12 С, а вторая обусловлена радиационным поглощением нейтронов ядрами 3Не (точнее, энергия 20.58 МэВ представляет собой нижнюю границу энергии испускаемых в этом процессе 7-квантов). Предварительные расчеты привели к оценке отношения концентраций снизу: 3Не/1Н > 1 • 10^4 в данной вспышке. Этот расчет был
основан на ядерной реакции радиационного поглощения нейтронов гелием-3: п + 3Не -¥ 4Не + 7. Расчет носит оценочный характер и основан только на сведениях о сечениях ядерных реакций и спектрах нейтронов. Найдена лишь нижняя граница возможных значений отношения 3Не/1Н.
Расчет же концентрации, проведенный в настоящей работе, связан с другим взаимодействием нейтрона с гелием-3 — безрадиационным поглощением нейтронов этим же нуклидом: п + 3Не -¥ 3Н + р. Эта реакция является конкурирующей при расчетах радиационного захвата нейтронов водородом с генерацией 7-линии с энергией кванта 2.223 МэВ. В расчете, связанном с двумя последними процессами, в отличие от первого способа, учтены все основные обстоятельства, влияющие на результирующее 7-излучение в нейтронной линии и при этом также и на сопутствующий процесс безрадиационного поглощения нейтронов ядрами 3Не. В качестве одного из результатов мы находим конкретное значение содержания 3Не вместе с погрешностью. Несмотря на различие методов, важно, что результаты обоих методов определения содержания 3Не в исследуемой вспышке не противоречат друг другу.
В работе [28] приведены и другие свидетельства, указывающие на повышенное содержание ядер 3 Не в процессе этой вспышки. В частности, сравнимая по интенсивности с линией железа 0.84 МэВ особенность в спектре по данным RHESSI для той же вспышки в области 0.9-1.0 МэВ может быть интерпретирована как линия 0.937 МэВ из комплекса (0.937, 1.04 и 1.08 МэВ), образующегося в прямой реакции на 3 Не [29]: 3Не +160 -¥ р + 18F*, и ее присутствие также указывает на повышенное содержание 3Не в области вспышки.
Имеются, кроме того, данные прямых измерений содержания 3Не в солнечном ветре и в составе энергичных солнечных частиц. Так, во время протонного события 20 января 2005 г. по измерениям на КА АСЕ, SAM-РЕХ и GOES-11 у Земли было получено отношение 4Не/1Н~0.2 при энергии частиц £>0.5 МэВ/нук-лон [1]. С другой стороны, по измерениям на КА АСЕ (точка Лагранжа LI) был оценен верхний предел для отношения 3Не/4Не < 7.7 • 10^4 для интервала 0.5-2.0 МэВ/нуклон [1, 30]. Эти оценки дают отношение 3Не/1Н < 1.54 • 10^4 вблизи орбиты Земли. Этот верхний предел согласуется с найденным в настоящей работе отношением, полученным из анализа временного профиля линии 2.223 МэВ вспышки 20 января 2005 г. Найденное соотношение 3Не/1Н = (1.4 ±0.15) • 10^4 можно рассматривать как независимое от метода примененного в [28] подтверждение повышенного содержания 3Не в процессе вспышки 20 января 2005 г. на основании анализа 7-излучения в линиях 2.223, 4.44 и 6.12 МэВ. Таким образом, заключение о повышенном содержании 3Не не только объясняет наблюдаемый временной профиль 7-линии 2.223 МэВ, но и согласуется с аналогичными выводами, полученными на основании анализа других данных.
Заключение
Изучение временного профиля 7-излучения солнечной вспышки 20 января 2005 г. в линии 2.223 МэВ
показало, что для его интерпретации необходимо допустить повышенное содержание изотопа 3Не в области ядерных реакций до отношения 3Не/1Н = (1.4 ±0.15) х х 10^4. При этом показано, что исследуемая вспышка 20 отличается помимо своей мощности определенными особенностями, изучение которых открывает возможности для более глубокого анализа электромагнитных и ядерных процессов в атмосфере Солнца.
Авторы благодарны В. Т. Ворончеву за консультацию по реакциям радиационного захвата между нуклонами и лёгкими ядрами и JI. И. Мирошниченко за прочтение первоначального варианта рукописи и советы.
Работа выполнена при финансовой поддержке РФФИ (грант 08-02-92208).
Список литературы
1. Mewaldt, R.A., Looper, M. D., Cohen С.M.S. et al. // Proc. 29th ICRC. Pune, 2005. 1. SH12. P. 101
2. Share G., Murphy R.J. Smith D. et al. // В AAS. 2006. 38. P. 255.
3. Струминский A.B. // Письма в Астрой, жури. 2006. 32, № 10. С. 767.
4. Кузнецов В.Д. // Астрой, вестник. 2005. 39, № 6. С. 485.
5. Солнечно-земная физика / Под ред. В. Д. Кузнецова. М., 2009.
6. Kuznetsov S.N., Kudela К., Ryumin S. P. et al. 11 Adv. Space Res. 2002. 30, N 7. P. 1857.
7. Кузнецов С.H., Курт В.Г., Юшков Б.Ю. и др. // Солнечно-земная физика / Под ред. В. Д. Кузнецова. М., 2009. С. 308.
8. Kuznetsov S.N., Kurt V.G., Yushkov B.Yu. et al. 11 Contributions of the Astronomical Observatory Skalnaté Pleso. 2006. 36, N 2. P. 85.
9. Arkhangelskaja I.V., Kotou Yu.D., Arkhangelsky A.I., Glyanenko A.S. 11 Adv. Space Res. 2009. 43, N 4. P. 542.
10. Котов Ю.Д., Архангельская И.В., Архангельский А.И. и др. // Солнечно-земная физика. / Под ред. В. Д. Кузнецова. М., 2009. С. 178.
11. Arkhangelskaja I.V., Arkhangelsky А.I., Kotov Yu.D. et al. // Adv. Space Res. 2009. 43, N 4. P. 589.
12. Архангельская И.В., Архангельский А.И., Котов Ю.Д. и др. // Астрон. вестник. 2008. 42, № 4. С. 373.
13. Гляненко A.C., Котов Ю. Д., Павлов В.Н. и др. // Приб. и техн. эксперимента. 1999. № 5. С. 16.
14. Troitskaia Е. V, Arkhangelskaja I. V, Miroshnichenko L.I, Arkhangelsky A.I 11 Adv. Space Res. 2009. 43, N 4. P. 547.
15. Кужевский Б.М., Троицкая E.B. Распространение нейтронов от солнечной вспышки в атмосфере Солнца. Препринт НИПЯФ МГУ 89-28/105. М., 1989.
16. Kuzhevskij В.М., Kuznetsov S.N., Troitskaia E.V. 11 31st COSPAR Scientific Assembly. Birmingham, 1996. 11 Adv. Space Res. 1998. 22. P. 1141.
17. Hua X.-M., Lingenfelter R.E. 11 Solar Phys. 1987. 107. P. 351.
18. Hua X.-M., Kozlovsky В., Lingenfelter R.E. et al. 11 Astrophys. J. Suppl. 2002. 140. P. 563.
19. Prince T.A., Forrest DJ., Chupp E.L. et al. 11 Proc. 18th Int. Cosmic Ray Conf. 1983. 4. P. 9.
20. Кужевский Б.М., Мирошниченко Л. И., Троицкая E.B. 11 Астрон. журнал. 2005. 82, № 7. С. 637.
21. Troitskaia E.V., Gan W.Q., Kuzhevskij B.M., Miroshnichenko L.I. // Solar Phys. 2007. 242. P. 87.
22. Kuzhevskij B.M., Kuznetsov S.N., Troitskaia E.V. // J. Moscow Phys. Soc. 1996. 6, N 4. P. 381.
23. Gingerich O., Noyes R.W., Kalkofen W., Сипу Y. // Solar Phys. 1971. 18. P. 347.
24. Spruit H.C. // Solar Phys. 1974. 34. P. 277.
25. Murphy R.J., Kozlovsky В., Share G.H. et al. // Astrophys. J. Suppl. 2007. 168, N 1. P. 167.
26. Kocharov L.G., Kocharov G.E. 11 Space Sei. Rev. 1984. 38. P. 89.
27. Hua X.-M., Kozlovsky В., Lingenfelter R.E. et al. // Astrophys. J. Suppl. 2002. 140. P. 563.
28. Arkhangelskaja I.V., Arkhangelsky A.I, Troitskaya E.V, Miroshnichenko L.I // Adv. Space Res. 2009. 43, N 4. P. 594.
29. Gan W.Q. // Chinese J. Astron. Astrophys. 2002. 26. P. 255.
30. Desai M.I, Mason G.M., Gold R.E. et al. // Astrophys. J. 2006. 649, N 1. P. 470.
Diagnostics of neutron interaction region in the process of 20 January 2005 solar flare by means of flare gamma-lines
E.V. Troitskaia'", I.V. Arkhangelskaja2,b, A.I. Arkhangelsky2
1D. V. Skobeltsyn Institute of Nuclear Physics, M. V. Lomonosov Moscow State University, Moscow 119991, Russia.
2 Department of Experimental Physics and Cosmophysics, National Research Nuclear University «MEPhh,
31 Kashirskoe road, Moscow 115409, Russia.
E-mail: a [email protected], b [email protected].
Basing on the analysis of the capture line (2.223 MeV) counting rate from the extreme solar event of 20 January 2005, we investigate the density of surrounding medium and 3He content in the area of 2.223 MeV gamma-line production. For analysis we have used AVS-F apparatus data of SONG-D detector onboard CORONAS-F satellite on counting rates of gamma-quanta in the line 2.223 MeV as well as in the lines in the range of 7.4 MeV. It is shown that for creating the appropriate time profile model of this solar flare it is necessary to suppose an enlarged content of 3He. In the area of 2.223 MeV gamma-line production the ratio of n(3He)/n('H) is found to be of (1.4±0.15) • 10^4, although the averaged value on the usual flares is about 2- 10^5. Enlarged content of 3He may be also considered as an additional indirect argument for possible existence of gamma-emission from reaction of radiative absorption of neutrons by 3He isotope.
Keywords: solar flare, gamma rays, nuclear reactions, density of solar atmosphere, 3 He isotope content. PACS: 96.60.qe, 95.85.Pw. Received 7 October 2011.
English version: Moscow University Physics Bulletin 3(2012).
Сведения об авторах
1. Троицкая Евгения Викторовна — канд. физ.-мат. наук, науч. сотрудник; тел.: (495) 939-58-68, e-mail: [email protected].
2. Архангельская Ирина Владимировна — ассистент; тел.: (495) 788-56-99, доб. 83-39, e-mail: [email protected].
3. Архангельский Андрей Игоревич — вед. инженер; тел.: (495) 788-56-99, доб. 96-79, e-mail: [email protected].