Скорости деполяризации и передачи поляризации для столкновений C$_2$ $(X ^1Σ^+_g, a ^3Π_u)$ + H$(^2S_{1/2})$ в фотосфере Солнца

Основная информация

• ID статьи: 2501.00763
• Название: Скорости деполяризации и передачи поляризации для столкновений C$_2$ $(X ^1Σ^+_g, a ^3Π_u)$ + H$(^2S_{1/2})$ в фотосфере Солнца
• Авторы: S. Qutub, Y.N. Kalugina, M. Derouich
• Классификация: astro-ph.SR physics.atom-ph physics.chem-ph physics.plasm-ph
• Дата публикации: 3 января 2025 г. (рукопись Astronomy & Astrophysics)
• Ссылка на статью: https://arxiv.org/abs/2501.00763

Аннотация

Настоящая работа является продолжением серии исследований деполяризации столкновениями молекулярных спектральных линий Солнца (таких как MgH, CN и C$_2$). Исследование сосредоточено на солнечной молекуле C$_2$, которая, хотя и имеет неярко выраженный профиль интенсивности и практически невидима, демонстрирует значительный профиль рассеянной поляризации. Из-за острого дефицита данных о столкновениях интерпретация поляризации C$_2$ на основе магнитного поля остаётся неполной. Целью данной работы является точный расчёт скоростей деполяризации и передачи поляризации при изотропных столкновениях C$_2$ $(X ^1Σ^+_g, a ^3Π_u)$ с атомами водорода H$(^2S_{1/2})$. Исследование использует пакет MOLPRO для получения поверхностей потенциальной энергии электронных состояний C$_2$, код MOLSCAT для изучения квантовой динамики и представление сечений столкновений и скоростей в неприводимом тензорном базисе. Аналитические выражения для температурной и молекулярной зависимости полного углового момента были определены с помощью техники генетического программирования.

Научный контекст и мотивация

Постановка проблемы

1. Наблюдаемые явления: Высокоточные спектрополяриметрические наблюдения телескопов выявили молекулярные спектральные линии MgH, CN и C$_2$ с беспрецедентным пространственным и спектральным разрешением, где молекула C$_2$ демонстрирует значительный профиль рассеянной поляризации
2. Потребность в диагностике магнитного поля: Эффект Ганле молекулярных поляризованных линий предоставляет благоприятную возможность для определения неразрешённых пространственных магнитных полей, особенно для линий системы Swan C$_2$ (d$^3Π_u$ – a$^3Π_u$) вблизи 5141 Å
3. Проблема отсутствия данных: Различия в результатах определения напряженности магнитного поля в основном обусловлены отсутствием учёта эффектов столкновений

Значимость исследования

1. Совершенствование теории: Столкновения конкурируют с эффектом Ганле-деполяризации магнитным полем, игнорирование столкновений существенно влияет на точное определение магнитного поля
2. Приложения в солнечной физике: Атомы водорода имеют высокую плотность в фотосфере, и эффекты их столкновений критичны для эффекта Ганле
3. Методологический вклад: Предоставление базовых данных о столкновениях для интерпретации магнитного поля молекулярных спектральных линий

Ограничения существующих методов

1. Дефицит данных: Скорости деполяризации и передачи поляризации столкновениями для молекулы C$_2$ практически полностью отсутствуют
2. Приблизительные предположения: Существующие модели часто игнорируют поляризацию нижних энергетических уровней, необходимо оценить ограничения такого приближения
3. Вычислительная сложность: Симметрия однородных молекул усложняет задачу столкновений

Основные вклады

1. Первый расчёт: Впервые произведён точный расчёт квантовых скоростей деполяризации и передачи поляризации для столкновений C$_2$ $(X ^1Σ^+_g, a ^3Π_u)$ + H$(^2S_{1/2})$
2. Построение поверхностей потенциальной энергии: Рассчитаны поверхности потенциальной энергии трёх электронных состояний (1 $^2A'$, 2 $^2A'$, 2 $^2A''$) с использованием пакета MOLPRO
3. Аналитические выражения: Получены аналитические выражения для температурной зависимости (2000–15000 К) и зависимости от углового момента с помощью техники генетического программирования
4. Значение для солнечной физики: Доказано, что изотропные столкновения могут только частично деполяризовать нижний уровень спектральных линий C$_2$, а не полностью их деполяризовать
5. Открытые данные: Предоставлены доступные в интернете данные сечений IOS для использования научным сообществом

Подробное описание методов

Определение задачи

Расчёт для молекулы C$_2$ в двух наиболее низких электронных состояниях (основное состояние X$^1Σ^+_g$ и первое возбуждённое состояние a$^3Π_u$) при столкновении с атомом водорода в основном состоянии:

• Скоростей деполяризации $D^k(j,T)$
• Скоростей передачи поляризации $D^k(j \to j',T)$

где $k$ — порядок тензора, $j$ — полный угловой момент молекулы, $T$ — температура.

Теоретическая база

1. Система координат и расчёт поверхностей потенциальной энергии

Использована система координат Якоби $(R, r_{C_2}, θ)$:

• $R$: вектор между центром масс C$_2$ и атомом водорода
• $θ$: угол вращения атома водорода вокруг молекулы C$_2$
• $r_{C_2} = 2.348 a_0$: длина связи C–C (приближение жёсткого ротатора)

Поверхности потенциальной энергии рассчитаны с использованием волновых функций многоконфигурационного взаимодействия конфигураций (MRCI+Q):

• C$_2$(X$^1Σ^+_g$) + H($^2S$) → состояние 1 $^2A'$
• C$_2$(a$^3Π_u$) + H($^2S$) → состояния 2 $^2A'$ и 2 $^2A''$

2. Обработка квантовой динамики

Применено приближение бесконечного порядка внезапности (IOS) для решения задачи столкновения, применимое при температурах Солнца. Использована схема развязки Кори и Александера (1985):

$N = R + L$$j = N + S_d$

где $R$ — угловой момент вращения, $L$ — орбитальный угловой момент электрона, $S_d$ — спин молекулы, $N$ — полный угловой момент вращения-орбиты.

3. Расчёт тензорных сечений

Выражение для сечения передачи поляризации в приближении IOS:

$σ^k_{IOS}(el, j → j', E) = \sum_K (-1)^{k+j+j'+K}(2R+1)(2R'+1) \times$$\frac{(2j'+1)(2j+1)(2N+1)(2N'+1)}{} \left\{ \begin{array}{ccc} j & j' & K \\ j' & j & k \end{array} \right\} \left\{ \begin{array}{ccc} R & R' & K \\ N' & N & L \end{array} \right\}^2 \times$$\left\{ \begin{array}{ccc} N & N' & K \\ j' & j & S_d \end{array} \right\}^2 \left( \begin{array}{ccc} R' & R & K \\ 0 & 0 & 0 \end{array} \right)^2 σ(el, 0 → K, E)$

Технологические инновации

1. Преимущества метода IOS: По сравнению с методом полной связи, IOS генерирует полный набор обобщённых таблиц сечений, эффективно рассчитывая скорости столкновений для всех переходов и порядков тензора
2. Аппроксимация генетическим программированием: Инновационное использование техники генетического программирования для получения двухпеременных аналитических выражений: $D^k_{GP} = n_H \times 10^{-10} \frac{\sum_i a^k_i j^{α^k_i} T^{β^k_i}}{\sum_i b^k_i j^{γ^k_i} T^{δ^k_i}}$
3. Симметрия однородной молекулы: Учтены ограничения симметрии однородной молекулы C$_2$, допускающие только переходы с чётным $Δj$

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

Параметры расчёта

• Диапазон энергий: 50 – 40 000 см$^{-1}$
• Диапазон температур: 2 000 – 15 000 К
• Диапазон углового момента: $j$ = 0–60
• Порядки тензора: $K$ = 0–158 (чётные)
• Сетка поверхности потенциальной энергии: 84 значения $R$ ($a_0$ ≤ $R$ ≤ 50 $a_0$), 51 значение $θ$ (0° ≤ $θ$ ≤ 90°)

Программные инструменты

• MOLPRO 2010: расчёт поверхностей потенциальной энергии
• MOLSCAT: расчёты квантовой динамики
• Генетическое программирование: аппроксимация аналитических выражений

Методы верификации

Точность метода IOS проверена сравнением с расчётами полной связи (CC), различие составляет менее 5% при $T$ = 1300 К.

Результаты исследования

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

1. Результаты для состояния X$^1Σ^+_g$

• Скорости деполяризации возрастают с увеличением температуры и уменьшаются с увеличением $j$
• Скорости деполяризации для $k=2$ (выравнивание) превышают скорости для $k=1$ (ориентация)
• Ошибка аппроксимации генетическим программированием менее 5%

2. Результаты для состояния a$^3Π_u$

• Для $j \gtrsim 5$ скорости деполяризации уменьшаются с увеличением $j$
• Скорости передачи поляризации значительно возрастают при малых значениях $j$, при $j \gtrsim 20$ рост замедляется
• Три компоненты $j$ ($j=N-1, N, N+1$) демонстрируют различное поведение деполяризации

3. Числовые примеры

Таблица 1 показывает сравнение линейных скоростей деполяризации с обратным временем жизни при типичных условиях фотосферы:

Анализ эффектов столкновений

1. Частичная деполяризация: При типичной плотности фотосферы столкновения могут только частично деполяризовать нижний уровень спектральных линий C$_2$
2. Зависимость от плотности: При $n_H = 10^{15}$ см$^{-3}$ скорость деполяризации составляет примерно 1/3 обратного времени жизни; при $n_H = 10^{16}$ см$^{-3}$ — примерно в 3 раза больше
3. Эффект углового момента: Линии R-ветви более подвержены влиянию столкновений, чем линии P-ветви, поскольку последние имеют большие значения $j$

Связанные работы

Исторический контекст

1. Развитие наблюдений: Высокоточные поляризационные наблюдения Stenflo (1994), Gandorfer (2000) и др.
2. Теоретические исследования: Теория магнитной диагностики Berdyugina & Fluri (2004), Milić & Faurobert (2012)
3. Предыдущие работы: Аналогичные исследования молекул MgH и CN, проведённые Qutub et al. (2020, 2021)

Развитие методологии

1. Метод IOS: Теоретические основы Pack (1972, 1974), Alexander и др.
2. Тензорная форма: Схема развязки Corey & Alexander (1985)
3. Молекулярные столкновения: Соответствующие исследования Green (1994), Dagdigian & Alexander (2009)

Заключение и обсуждение

Основные выводы

1. Значимость столкновений: Изотропные столкновения с атомами водорода являются фундаментальным элементом понимания поляризации C$_2$
2. Частичная деполяризация: Доказаны ограничения приближения, игнорирующего поляризацию нижних энергетических уровней
3. Температурная и угловая зависимость: Предоставлены полные аналитические выражения и база данных

Ограничения

1. Приближение IOS: Точность снижается при низких температурах, но достаточна при температурах Солнца
2. Жёсткий ротатор: Игнорирование возбуждения колебаний, но обоснованно при солнечных условиях
3. Развязка спина: Игнорирование эффектов спина атома водорода, необходимое приближение

Значение для солнечной физики

1. Магнитная диагностика: Предоставление необходимых данных о столкновениях для измерения магнитного поля системы Swan C$_2$
2. Улучшение моделей: Помощь в устранении различий в измерениях напряженности магнитного поля
3. Физика фотосферы: Углубление понимания молекулярных физических процессов в фотосфере

Глубокая оценка

Достоинства

1. Заполнение пробела: Впервые предоставлены систематические данные о столкновениях молекулы C$_2$
2. Передовые методы: Сочетание квантово-химических расчётов, квантовой динамики и техники генетического программирования
3. Практическая ценность: Предоставление онлайн-базы данных и аналитических выражений для применения
4. Достаточная верификация: Проверка надёжности метода IOS сравнением с методом полной связи

Недостатки

1. Диапазон температур: Ограниченная точность метода IOS при низких температурах
2. Приблизительные предположения: Несколько приближений (жёсткий ротатор, развязка спина) могут влиять на точность
3. Экспериментальная верификация: Отсутствие экспериментальных данных для прямой проверки расчётных результатов

Влияние

1. Научный вклад: Предоставление важных данных для солнечной физики и молекулярной астрофизики
2. Прикладная ценность: Прямое применение к наблюдениям магнитного поля Солнца и теоретическому моделированию
3. Методология: Служит примером для исследований аналогичных молекулярных систем

Области применения

1. Солнечная физика: Диагностика магнитного поля фотосферы и моделирование спектральных линий C$_2$
2. Звёздные атмосферы: Анализ молекулярных спектральных линий других звёзд
3. Теоретические исследования: Физика молекулярных столкновений и теория поляризации

Библиография

Статья содержит обширную библиографию соответствующих работ, включая:

• Gandorfer (2000): Наблюдения второго солнечного спектра
• Berdyugina & Fluri (2004): Теория магнитной диагностики C$_2$
• Corey & Alexander (1985): Теоретическая база молекулярных столкновений
• Qutub et al. (2020, 2021): Предыдущие исследования MgH и CN

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