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

Basic Information

• Paper ID: 2501.00763
• Title: Depolarization and polarization transfer rates for the C$_2$ $(X ^1Σ^+_g, a ^3Π_u)$ + H$(^2S_{1/2})$ collisions in the solar photosphere
• Authors: S. Qutub, Y.N. Kalugina, M. Derouich
• Classification: astro-ph.SR physics.atom-ph physics.chem-ph physics.plasm-ph
• Publication Date: January 3, 2025 (Astronomy & Astrophysics manuscript)
• Paper Link: https://arxiv.org/abs/2501.00763

Abstract

This paper continues a series of investigations on collisional depolarization of solar molecular spectral lines (such as MgH, CN, and C$_2$). The research focuses on the solar molecule C$_2$, which, although exhibiting inconspicuous intensity profiles and being nearly invisible, demonstrates significant scattering polarization profiles. Due to severe scarcity of collision data, the interpretation of C$_2$ polarization based on magnetic fields remains incomplete. This work aims to precisely calculate depolarization and polarization transfer rates for isotropic collisions between C$_2$ $(X ^1Σ^+_g, a ^3Π_u)$ and hydrogen atoms H$(^2S_{1/2})$. The study employs the MOLPRO software package to obtain potential energy surfaces of C$_2$ electronic states, uses the MOLSCAT code to investigate quantum dynamics, and employs tensor irreducible basis representation for collision cross sections and rates. Analytical expressions for temperature and molecular total angular momentum dependence are determined through genetic programming techniques.

Research Background and Motivation

Problem Background

1. Observational Phenomena: High-precision spectropolarimetric telescope observations of MgH, CN, and C$_2$ molecular spectral lines with unprecedented spatial and spectral resolution, where C$_2$ molecules exhibit significant scattering polarization profiles
2. Magnetic Field Diagnostic Requirements: The Hanle effect of molecular polarization lines provides an excellent opportunity to determine spatially unresolved magnetic fields, particularly for C$_2$ Swan system lines (d$^3Π_u$ – a$^3Π_u$) near 5141 Å
3. Data Deficiency Problem: Discrepancies in magnetic field strength results are primarily due to lack of consideration of collision effects

Research Significance

1. Theoretical Refinement: Collisional and magnetic field Hanle depolarization effects compete; neglecting collisions significantly affects precise magnetic field determination
2. Solar Physics Applications: Hydrogen atoms have high density in the photosphere, and their collision effects are crucial to the Hanle effect
3. Methodological Contribution: Provides fundamental collision data for magnetic field interpretation of molecular spectral line polarization

Limitations of Existing Methods

1. Data Scarcity: Collision depolarization and polarization transfer rates for C$_2$ molecules are almost entirely lacking
2. Approximate Assumptions: Existing models often neglect lower level polarization, and the limitations of such approximations require assessment
3. Computational Complexity: Symmetry constraints of homonuclear molecules make collision problems more complex

Core Contributions

1. First Calculation: First precise calculation of quantum depolarization and polarization transfer rates for C$_2$ $(X ^1Σ^+_g, a ^3Π_u)$ + H$(^2S_{1/2})$ collisions
2. Potential Energy Surface Construction: Calculated potential energy surfaces for three electronic states (1 $^2A'$, 2 $^2A'$, 2 $^2A''$) using MOLPRO software package
3. Analytical Expressions: Obtained analytical expressions for temperature (2000-15000 K) and angular momentum dependence through genetic programming techniques
4. Solar Physics Significance: Demonstrated that isotropic collisions can only partially depolarize the lower state of C$_2$ spectral lines, not completely depolarize them
5. Open Data Access: Provided online-accessible IOS cross-section data for community use

Detailed Methodology

Task Definition

Calculate for C$_2$ molecules in two lowest electronic states (X$^1Σ^+_g$ ground state and a$^3Π_u$ first excited state) colliding with hydrogen atoms in ground state:

• Depolarization rates $D^k(j,T)$
• Polarization transfer rates $D^k(j \to j',T)$

where $k$ is the tensor rank, $j$ is the molecular total angular momentum, and $T$ is temperature.

Theoretical Framework

1. Coordinate System and Potential Energy Surface Calculation

Employs Jacobi coordinates $(R, r_{C_2}, θ)$:

• $R$: intermolecular vector connecting C$_2$ center of mass with hydrogen atom
• $θ$: rotation angle of hydrogen atom around C$_2$ molecule
• $r_{C_2} = 2.348 a_0$: C-C bond length (rigid rotor approximation)

Uses multi-reference configuration interaction wavefunctions (MRCI+Q) to calculate potential energy surfaces:

• C$_2$(X$^1Σ^+_g$) + H($^2S$) → 1 $^2A'$ state
• C$_2$(a$^3Π_u$) + H($^2S$) → 2 $^2A'$ and 2 $^2A''$ states

2. Quantum Dynamics Treatment

Employs infinite order sudden approximation (IOS) for collision problems, applicable to solar temperature conditions. Uses Corey & Alexander (1985) decoupling scheme:

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

where $R$ is rotational angular momentum, $L$ is electronic orbital angular momentum, $S_d$ is molecular spin, and $N$ is total rotational-orbital angular momentum.

3. Tensor Cross Section Calculation

IOS polarization transfer cross section expression:

$σ^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)$

Technical Innovations

1. IOS Method Advantages: Compared to close-coupling methods, IOS generates comprehensive generalized cross section tables and efficiently calculates collision rates for all transitions and tensor ranks
2. Genetic Programming Fitting: Innovatively uses genetic programming techniques to obtain bivariate analytical expressions: $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. Homonuclear Molecule Symmetry: Accounts for symmetry constraints of C$_2$ homonuclear molecules, allowing only even $Δj$ transitions

Experimental Setup

Computational Parameters

• Energy Range: 50 - 40,000 cm$^{-1}$
• Temperature Range: 2,000 - 15,000 K
• Angular Momentum Range: $j$ = 0-60
• Tensor Ranks: $K$ = 0-158 (even numbers)
• Potential Energy Surface Grid: 84 $R$ values ($a_0$ ≤ $R$ ≤ 50 $a_0$), 51 $θ$ values (0° ≤ $θ$ ≤ 90°)

Software Tools

• MOLPRO 2010: Potential energy surface calculations
• MOLSCAT: Quantum dynamics calculations
• Genetic Programming: Analytical expression fitting

Verification Methods

Validates IOS method accuracy through comparison with close-coupling (CC) calculations, with differences less than 5% at $T$ = 1300 K.

Experimental Results

Main Results

1. X$^1Σ^+_g$ State Results

• Depolarization rates increase with temperature and decrease with $j$
• Depolarization rates for $k=2$ (alignment) exceed those for $k=1$ (orientation)
• Genetic programming fitting errors less than 5%

2. a$^3Π_u$ State Results

• For $j \gtrsim 5$, depolarization rates decrease with increasing $j$
• Polarization transfer rates increase significantly at small $j$ values, growing slowly for $j \gtrsim 20$
• Three $j$ components ($j=N-1, N, N+1$) exhibit different depolarization behavior

3. Numerical Examples

Table 1 shows linear depolarization rates compared with inverse lifetimes under typical photospheric conditions:

Collision Effects Analysis

1. Partial Depolarization: Under typical photospheric densities, collisions can only partially depolarize the lower state of C$_2$ spectral lines
2. Density Dependence: At $n_H = 10^{15}$ cm$^{-3}$, depolarization rates are approximately 1/3 of inverse lifetimes; at $n_H = 10^{16}$ cm$^{-3}$ approximately 3 times
3. Angular Momentum Effects: R-branch lines are more susceptible to collisions than P-branch lines because the latter have larger $j$ values

Related Work

Historical Background

1. Observational Development: High-precision polarimetric observations by Stenflo (1994), Gandorfer (2000), and others
2. Theoretical Research: Magnetic field diagnostic theory by Berdyugina & Fluri (2004), Milić & Faurobert (2012)
3. Prior Work: Similar studies on MgH and CN molecules by Qutub et al. (2020, 2021)

Methodological Development

1. IOS Method: Theoretical foundations by Pack (1972, 1974), Alexander et al.
2. Tensor Formalism: Decoupling scheme by Corey & Alexander (1985)
3. Molecular Collisions: Related research by Green (1994), Dagdigian & Alexander (2009)

Conclusions and Discussion

Main Conclusions

1. Collision Importance: Isotropic hydrogen atom collisions are fundamental elements for understanding C$_2$ polarization
2. Partial Depolarization: Demonstrates limitations of approximations neglecting lower level polarization
3. Temperature and Angular Momentum Dependence: Provides complete analytical expressions and database

Limitations

1. IOS Approximation: Reduced accuracy at low temperatures, but sufficiently accurate at solar temperatures
2. Rigid Rotor: Neglects vibrational excitation, but reasonable under solar conditions
3. Spin Decoupling: Neglects hydrogen atom spin effects, a necessary approximation

Solar Physics Significance

1. Magnetic Field Diagnostics: Provides necessary collision data for magnetic field measurements of C$_2$ Swan system
2. Model Improvement: Helps eliminate discrepancies in magnetic field strength measurements
3. Photospheric Physics: Enhances understanding of molecular physics processes in the photospheric layer

In-Depth Evaluation

Strengths

1. Filling Gaps: First systematic collision data for C$_2$ molecules
2. Advanced Methods: Combines quantum chemistry calculations, quantum dynamics, and genetic programming techniques
3. Practical Value: Provides online database and analytical expressions for convenient application
4. Thorough Verification: Validates IOS method reliability through comparison with close-coupling methods

Weaknesses

1. Temperature Range: IOS method has limited accuracy at low temperatures
2. Approximate Assumptions: Multiple approximations (rigid rotor, spin decoupling) may affect accuracy
3. Experimental Verification: Lacks direct experimental data to verify computational results

Impact

1. Academic Contribution: Provides important data for solar physics and molecular astrophysics
2. Application Value: Directly serves solar magnetic field observations and theoretical modeling
3. Methodology: Provides exemplary approach for similar molecular systems

Applicable Scenarios

1. Solar Physics: Photospheric magnetic field diagnostics and C$_2$ spectral line modeling
2. Stellar Atmospheres: Molecular spectral line analysis in other stars
3. Theoretical Research: Molecular collision physics and polarization theory

References

The paper cites abundant relevant literature, including:

• Gandorfer (2000): Solar second spectrum observations
• Berdyugina & Fluri (2004): C$_2$ magnetic field diagnostic theory
• Corey & Alexander (1985): Molecular collision theory framework
• Qutub et al. (2020, 2021): Prior research on MgH and CN

This research provides an important theoretical foundation for solar molecular spectroscopy, particularly valuable for magnetic field diagnostic applications. Through precise quantum calculations and innovative data processing methods, it makes significant contributions to understanding physical processes in the solar photospheric layer.