2025-11-19T18:49:13.470427

Depolarization and polarization transfer rates for the C$_2$ $(X ^1Î£^+_g, a ^3Î _u)$ + H$(^2S_{1/2})$ collisions in the solar photosphere

Qutub, Kalugina, Derouich

This paper is a continuation of a series of studies investigating collisional depolarization of solar molecular lines like those of MgH, CN and C$_2$. It is focused on the case of the solar molecule C$_2$ which exhibits striking scattering polarization profiles although its intensity profiles are inconspicuous and barely visible. In fact, interpretation of the C$_2$ polarization in terms of magnetic fields is incomplete due to the almost complete lack of collisional data. This work aims at accurately computing the collisional depolarization and polarization transfer rates for the C$_2$~$(X ^1Î£^+_g, a ^3Î _u)$ by isotropic collisions with hydrogen atoms H~$(^2S_{1/2})$. We also investigate the solar implications of our findings. We utilize the MOLPRO package to obtain potential energy surfaces (PESs) for the electronic states $X ^1Î£^+_g$ and $a^3Î _u$ of C$_2$, and the MOLSCAT code to study the quantum dynamics of the C$_2$~$(X ^1Î£^+_g, a ^3Î _u)$ + H$(^2S_{1/2})$ systems. We use the tensorial irreducible basis to express the resulting collisional cross-sections and rates. Furthermore, sophisticated genetic programming techniques are employed to determine analytical expressions for the temperature and total molecular angular momentum dependence of these collisional rates. We obtain quantum depolarization and polarization transfer rates for the C$_2$ $(X ^1Î£^+_g, a ^3Î _u)$ + H$(^2S_{1/2})$ collisions in the temperature range T=2,000--15,000~K. We also determine analytical expressions giving these rates as functions of the temperature and total molecular angular momentum. In addition, we show that isotropic collisions with neutral hydrogen can only partially depolarize the lower state of C$_2$ lines, rather than completely. This highlights the limitations of the approximation of neglecting lower-level polarization while modeling the polarization of C$_2$ lines.

academic

太陽光球におけるC $_2$ $(X ^1Σ^+_g, a ^3Π_u)$ + H $(^2S_{1/2})$ 衝突の脱偏光および偏光転移速度

基本情報

論文ID: 2501.00763
タイトル: Depolarization and polarization transfer rates for the C $_2$ $(X ^1Σ^+_g, a ^3Π_u)$ + H $(^2S_{1/2})$ collisions in the solar photosphere
著者: S. Qutub, Y.N. Kalugina, M. Derouich
分類: astro-ph.SR physics.atom-ph physics.chem-ph physics.plasm-ph
発表日: 2025年1月3日（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コードを用いて量子動力学を研究し、テンソル既約基底を用いて衝突断面積と速度を表現した。遺伝的プログラミング技術により、温度と分子全角運動量依存性の解析表現式が決定された。

研究背景と動機

問題背景

観測現象: 高精度分光偏光望遠鏡により、MgH、CN、C $_2$ 分子スペクトル線が前例のない空間および分光分解能で観測され、特にC $_2$ 分子は顕著な散乱偏光プロファイルを示す
磁場診断の必要性: 分子偏光線のハンレ効果は、空間未分解磁場の決定に良好な機会を提供し、特にC $_2$ Swan系（d $^3Π_u$ – a $^3Π_u$ ）の5141 Åの近傍のスペクトル線が該当する
データ欠落問題: 磁場強度の結果に差異が生じており、主な原因は衝突効果の考慮不足である

研究の重要性

理論の完善: 衝突と磁場のハンレ脱偏光効果は競合関係にあり、衝突を無視すると磁場の正確な測定に大きな影響を与える
太陽物理への応用: 水素原子は光球層で高い密度を有しており、その衝突効果はハンレ効果に極めて重要である
方法論への貢献: 分子スペクトル線偏光の磁場解釈に基礎的な衝突データを提供する

既存方法の限界

データ不足: C $_2$ 分子の衝突脱偏光および偏光転移速度はほぼ完全に欠落している
近似仮定: 既存モデルはしばしば下準位偏光を無視しており、この近似の限界を評価する必要がある
計算の複雑性: 同核分子の対称性制限により、衝突問題はより複雑になる

核心的貢献

初の計算: C $_2$ $(X ^1Σ^+_g, a ^3Π_u)$ + H $(^2S_{1/2})$ 衝突の量子脱偏光および偏光転移速度を初めて正確に計算した
ポテンシャルエネルギー面の構築: MOLPROソフトウェアパッケージを用いて3つの電子状態（1 $^2A'$ 、2 $^2A'$ 、2 $^2A''$ ）のポテンシャルエネルギー面を計算した
解析表現式: 遺伝的プログラミング技術により、温度（2000～15000 K）および角運動量依存の解析表現式を得た
太陽物理的意義: 等方性衝突がC $_2$ スペクトル線の下準位を完全には脱偏光できず、部分的にしか脱偏光できないことを実証した
データの公開: オンラインでアクセス可能なIOS断面積データを提供し、コミュニティの利用を促進した

方法の詳細

タスク定義

C $_2$ 分子が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）を用いて衝突問題を処理し、太陽温度条件に適用可能である。Corey & Alexander (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)$

技術的革新点

IOS方法の利点: 密結合法と比較して、IOSは包括的な一般化断面積表を生成でき、すべての遷移およびテンソル階数の衝突速度を効率的に計算できる
遺伝的プログラミング適合: 遺伝的プログラミング技術を革新的に使用して、二変数解析表現式を得た： $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}}$
同核分子対称性: C $_2$ 同核分子の対称性制限を考慮し、 $Δj$ が偶数の遷移のみを許可した

実験設定

計算パラメータ

エネルギー範囲: 50～40,000 cm $^{-1}$
温度範囲: 2,000～15,000 K
角運動量範囲: $j$ = 0-60
テンソル階数: $K$ = 0-158（偶数）
ポテンシャルエネルギー面グリッド: 84個の $R$ 値（ $a_0$ ≤ $R$ ≤ 50 $a_0$ ）、51個の $θ$ 値（0° ≤ $θ$ ≤ 90°）

ソフトウェアツール

MOLPRO 2010: ポテンシャルエネルギー面計算
MOLSCAT: 量子動力学計算
遺伝的プログラミング: 解析表現式適合

脱偏光速度は温度の増加とともに増加し、 $j$ の増加とともに減少する
$k=2$ （配列）の脱偏光速度は $k=1$ （方向）より大きい
遺伝的プログラミング適合誤差は5%未満である

2. a $^3Π_u$ 状態の結果

$j \gtrsim 5$ に対して、脱偏光速度は $j$ の増加とともに減少する
偏光転移速度は小さい $j$ 値で顕著に増加し、 $j \gtrsim 20$ で増加が緩やかになる
3つの $j$ 成分（ $j=N-1, N, N+1$ ）は異なる脱偏光挙動を示す

3. 数値例

表1は典型的な光球条件下での線形脱偏光速度と逆寿命の比較を示す：

スペクトル線	$j$	$B_{lu}I(λ_{lu})$ (10 $^5$ s $^{-1}$ )	$D^2$ ( $n_H$ =10 $^{15}$ cm $^{-3}$ )	$D^2$ ( $n_H$ =10 $^{16}$ cm $^{-3}$ )
R1(14)	14	0.846	0.279	2.79
R2(13)	13	0.918	0.254	2.54
R3(12)	12	1.014	0.350	3.50

衝突効果の分析

部分的脱偏光: 典型的な光球密度では、衝突はC $_2$ スペクトル線の下準位を部分的にしか脱偏光できない
密度依存性: $n_H = 10^{15}$ cm $^{-3}$ では脱偏光速度は逆寿命の約1/3； $n_H = 10^{16}$ cm $^{-3}$ では約3倍である
角運動量効果: R分枝線はP分枝線より衝突の影響を受けやすい。これは後者がより大きな $j$ 値を有するためである