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

Tasas de despolarización y transferencia de polarización para colisiones C2_2 (X1Σg+,a3Πu)(X ^1Σ^+_g, a ^3Π_u) + H(2S1/2)(^2S_{1/2}) en la fotosfera solar

Información Básica

  • ID del artículo: 2501.00763
  • Título: Tasas de despolarización y transferencia de polarización para colisiones C2_2 (X1Σg+,a3Πu)(X ^1Σ^+_g, a ^3Π_u) + H(2S1/2)(^2S_{1/2}) en la fotosfera solar
  • Autores: S. Qutub, Y.N. Kalugina, M. Derouich
  • Clasificación: astro-ph.SR physics.atom-ph physics.chem-ph physics.plasm-ph
  • Fecha de publicación: 3 de enero de 2025 (manuscrito de Astronomy & Astrophysics)
  • Enlace del artículo: https://arxiv.org/abs/2501.00763

Resumen

Este artículo continúa una serie de investigaciones sobre despolarización por colisiones de líneas moleculares solares (como MgH, CN y C2_2). La investigación se centra en la molécula solar C2_2, que aunque presenta perfiles de intensidad poco notables y es casi invisible, exhibe perfiles de polarización por dispersión significativos. Debido a la grave escasez de datos de colisiones, la interpretación de la polarización de C2_2 basada en campos magnéticos sigue siendo incompleta. Este trabajo tiene como objetivo calcular con precisión las tasas de despolarización y transferencia de polarización para colisiones isotrópicas de C2_2 (X1Σg+,a3Πu)(X ^1Σ^+_g, a ^3Π_u) con átomos de hidrógeno H(2S1/2)(^2S_{1/2}). El estudio utiliza el paquete MOLPRO para obtener superficies de energía potencial de estados electrónicos de C2_2, el código MOLSCAT para investigar dinámicas cuánticas, y representación en base de tensores irreducibles para secciones transversales y tasas de colisión. Se determinaron expresiones analíticas para dependencias de temperatura y momento angular molecular total mediante técnicas de programación genética.

Antecedentes e Importancia de la Investigación

Contexto del Problema

  1. Fenómenos observados: Observaciones de telescopios de polarización espectral de alta precisión de líneas moleculares MgH, CN y C2_2 con resolución espacial y espectral sin precedentes, donde la molécula C2_2 exhibe perfiles de polarización por dispersión significativos
  2. Necesidad de diagnóstico de campos magnéticos: El efecto Hanle de líneas polarizadas moleculares proporciona una buena oportunidad para determinar campos magnéticos no resueltos espacialmente, particularmente líneas del sistema Swan de C2_2 (d3Πu^3Π_u – a3Πu^3Π_u) cerca de 5141 Å
  3. Problema de datos faltantes: Existen discrepancias en los resultados de intensidad de campo magnético, principalmente debido a la falta de consideración de efectos de colisión

Importancia de la Investigación

  1. Perfeccionamiento teórico: Las colisiones compiten con el efecto Hanle de despolarización magnética; ignorar las colisiones afecta significativamente la determinación precisa del campo magnético
  2. Aplicación en física solar: Los átomos de hidrógeno tienen densidades muy altas en la capa fotosférica, y sus efectos de colisión son cruciales para el efecto Hanle
  3. Contribución metodológica: Proporciona datos de colisión fundamentales para la interpretación de campos magnéticos en líneas moleculares polarizadas

Limitaciones de Métodos Existentes

  1. Escasez de datos: Las tasas de despolarización y transferencia de polarización por colisiones de la molécula C2_2 están casi completamente ausentes
  2. Suposiciones aproximadas: Los modelos existentes a menudo ignoran la polarización del nivel inferior, siendo necesario evaluar las limitaciones de esta aproximación
  3. Complejidad computacional: Las restricciones de simetría de moléculas homonucleares hacen que el problema de colisión sea más complejo

Contribuciones Principales

  1. Cálculo por primera vez: Cálculo preciso por primera vez de tasas cuánticas de despolarización y transferencia de polarización para colisiones C2_2 (X1Σg+,a3Πu)(X ^1Σ^+_g, a ^3Π_u) + H(2S1/2)(^2S_{1/2})
  2. Construcción de superficies de energía potencial: Cálculo de superficies de energía potencial para tres estados electrónicos (1 2A′^2A', 2 2A′^2A', 2 2A′′^2A'') utilizando el paquete MOLPRO
  3. Expresiones analíticas: Obtención de expresiones analíticas dependientes de temperatura (2000-15000 K) y momento angular mediante técnicas de programación genética
  4. Significado en física solar: Demostración de que las colisiones isotrópicas solo pueden despolarizar parcialmente el estado inferior de las líneas de C2_2, no completamente
  5. Acceso abierto a datos: Provisión de datos de secciones transversales IOS accesibles en línea para uso de la comunidad

Descripción Detallada de Métodos

Definición de la Tarea

Cálculo de las siguientes cantidades para colisiones de la molécula C2_2 en dos estados electrónicos más bajos (estado base X1Σg+^1Σ^+_g y primer estado excitado a3Πu^3Π_u) con átomos de hidrógeno en estado base:

  • Tasas de despolarización Dk(j,T)D^k(j,T)
  • Tasas de transferencia de polarización Dk(j→j′,T)D^k(j \to j',T)

donde kk es el orden del tensor, jj es el momento angular rotacional molecular total, y TT es la temperatura.

Marco Teórico

1. Sistema de Coordenadas y Cálculo de Superficies de Energía Potencial

Se adopta el sistema de coordenadas de Jacobi (R,rC2,θ)(R, r_{C_2}, θ):

  • RR: vector intermolecular conectando el centro de masa de C2_2 con el átomo de hidrógeno
  • θθ: ángulo de rotación del átomo de hidrógeno alrededor de la molécula C2_2
  • rC2=2.348a0r_{C_2} = 2.348 a_0: longitud de enlace C-C (aproximación de rotor rígido)

Se utilizan funciones de onda de interacción de configuraciones multireferencia (MRCI+Q) para calcular superficies de energía potencial:

  • C2_2(X1Σg+^1Σ^+_g) + H(2S^2S) → estado 1 2A′^2A'
  • C2_2(a3Πu^3Π_u) + H(2S^2S) → estados 2 2A′^2A' y 2 2A′′^2A''

2. Tratamiento de Dinámicas Cuánticas

Se adopta la aproximación de orden infinito repentino (IOS) para tratar el problema de colisión, aplicable en condiciones de temperatura solar. Se utiliza el esquema de desacoplamiento de Corey & Alexander (1985):

N=R+LN = R + Lj=N+Sdj = N + S_d

donde RR es el momento angular rotacional, LL es el momento angular orbital electrónico, SdS_d es el espín molecular, y NN es el momento angular rotacional-orbital total.

3. Cálculo de Secciones Transversales de Tensores

Expresión de secciones transversales de transferencia de polarización IOS:

σIOSk(el,j→j′,E)=∑K(−1)k+j+j′+K(2R+1)(2R′+1)×σ^k_{IOS}(el, j → j', E) = \sum_K (-1)^{k+j+j'+K}(2R+1)(2R'+1) \times(2j′+1)(2j+1)(2N+1)(2N′+1){jj′Kj′jk}{RR′KN′NL}2×\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{NN′Kj′jSd}2(R′RK000)2σ(el,0→K,E)\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)

Puntos de Innovación Técnica

  1. Ventajas del método IOS: En comparación con métodos de acoplamiento cerrado, IOS genera tablas de secciones transversales generalizadas completas, calculando eficientemente tasas de colisión para todas las transiciones y órdenes de tensor
  2. Ajuste de programación genética: Uso innovador de técnicas de programación genética para obtener expresiones analíticas bivariadas: DGPk=nH×10−10∑iaikjαikTβik∑ibikjγikTδikD^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. Simetría de moléculas homonucleares: Consideración de restricciones de simetría de la molécula homonuclear C2_2, permitiendo solo transiciones con ΔjΔj par

Configuración Experimental

Parámetros de Cálculo

  • Rango de energía: 50 - 40,000 cm−1^{-1}
  • Rango de temperatura: 2,000 - 15,000 K
  • Rango de momento angular: jj = 0-60
  • Órdenes de tensor: KK = 0-158 (números pares)
  • Malla de superficie de energía potencial: 84 valores de RR (a0a_0 ≤ RR ≤ 50 a0a_0), 51 valores de θθ (0° ≤ θθ ≤ 90°)

Herramientas de Software

  • MOLPRO 2010: Cálculo de superficies de energía potencial
  • MOLSCAT: Cálculo de dinámicas cuánticas
  • Programación genética: Ajuste de expresiones analíticas

Métodos de Verificación

Verificación de la precisión del método IOS mediante comparación con cálculos de acoplamiento cerrado (CC), con diferencias menores al 5% a TT = 1300 K.

Resultados Experimentales

Resultados Principales

1. Resultados del Estado X1Σg+^1Σ^+_g

  • Las tasas de despolarización aumentan con la temperatura y disminuyen con jj
  • La despolarización de k=2k=2 (alineación) es mayor que la de k=1k=1 (orientación)
  • El error de ajuste de programación genética es menor al 5%

2. Resultados del Estado a3Πu^3Π_u

  • Para j≳5j \gtrsim 5, las tasas de despolarización disminuyen con el aumento de jj
  • Las tasas de transferencia de polarización aumentan significativamente en valores pequeños de jj, con crecimiento lento para j≳20j \gtrsim 20
  • Los tres componentes de jj (j=N−1,N,N+1j=N-1, N, N+1) exhiben comportamientos de despolarización diferentes

3. Ejemplos Numéricos

La Tabla 1 muestra comparaciones de tasas de despolarización lineal con vidas medias inversas bajo condiciones fotosféricas típicas:

LíneajjBluI(λlu)B_{lu}I(λ_{lu}) (105^5s−1^{-1})D2D^2 (nHn_H=1015^{15}cm−3^{-3})D2D^2 (nHn_H=1016^{16}cm−3^{-3})
R1(14)140.8460.2792.79
R2(13)130.9180.2542.54
R3(12)121.0140.3503.50

Análisis de Efectos de Colisión

  1. Despolarización parcial: En densidades fotosféricas típicas, las colisiones solo pueden despolarizar parcialmente el estado inferior de las líneas de C2_2
  2. Dependencia de densidad: A nH=1015n_H = 10^{15} cm−3^{-3}, la tasa de despolarización es aproximadamente 1/3 de la vida media inversa; a nH=1016n_H = 10^{16} cm−3^{-3} es aproximadamente 3 veces
  3. Efecto del momento angular: Las líneas de rama R son más susceptibles a efectos de colisión que las líneas de rama P, porque estas últimas tienen valores de jj más grandes

Trabajos Relacionados

Antecedentes Históricos

  1. Desarrollo observacional: Observaciones de polarización de alta precisión de Stenflo (1994), Gandorfer (2000) y otros
  2. Investigación teórica: Teoría de diagnóstico de campos magnéticos de Berdyugina & Fluri (2004), Milić & Faurobert (2012)
  3. Trabajos anteriores: Investigaciones similares de Qutub et al. (2020, 2021) sobre moléculas MgH y CN

Desarrollo Metodológico

  1. Método IOS: Fundamentos teóricos de Pack (1972, 1974) y Alexander et al.
  2. Forma de tensor: Esquema de desacoplamiento de Corey & Alexander (1985)
  3. Colisiones moleculares: Investigaciones relacionadas de Green (1994), Dagdigian & Alexander (2009)

Conclusiones y Discusión

Conclusiones Principales

  1. Importancia de colisiones: Las colisiones isotrópicas de átomos de hidrógeno son elementos fundamentales para entender la polarización de C2_2
  2. Despolarización parcial: Demostración de las limitaciones de la aproximación que ignora la polarización del nivel inferior
  3. Dependencia de temperatura y momento angular: Provisión de expresiones analíticas completas y base de datos

Limitaciones

  1. Aproximación IOS: Precisión reducida a temperaturas bajas, pero suficiente a temperaturas solares
  2. Rotor rígido: Ignora excitación vibracional, pero es razonable bajo condiciones solares
  3. Desacoplamiento de espín: Ignora efectos de espín del átomo de hidrógeno, siendo una aproximación necesaria

Significado en Física Solar

  1. Diagnóstico de campos magnéticos: Proporciona datos de colisión necesarios para mediciones de campos magnéticos del sistema Swan de C2_2
  2. Mejora de modelos: Ayuda a eliminar discrepancias en mediciones de intensidad de campo magnético
  3. Física de la fotosfera: Mejora la comprensión de procesos de física molecular en la capa fotosférica

Evaluación Profunda

Fortalezas

  1. Llenado de vacíos: Primera provisión de datos de colisión sistemáticos para la molécula C2_2
  2. Métodos avanzados: Combinación de cálculos de química cuántica, dinámicas cuánticas y técnicas de programación genética
  3. Valor práctico: Provisión de base de datos en línea y expresiones analíticas para aplicaciones
  4. Verificación suficiente: Validación de la confiabilidad del método IOS mediante comparación con métodos de acoplamiento cerrado

Insuficiencias

  1. Rango de temperatura: Precisión limitada del método IOS a temperaturas bajas
  2. Suposiciones aproximadas: Múltiples aproximaciones (rotor rígido, desacoplamiento de espín) pueden afectar la precisión
  3. Verificación experimental: Falta de datos experimentales para verificar directamente los resultados de cálculo

Impacto

  1. Contribución académica: Proporciona datos importantes para física solar y astrofísica molecular
  2. Valor de aplicación: Servicio directo a observaciones de campos magnéticos solares y modelado teórico
  3. Metodología: Proporciona un ejemplo para investigaciones de sistemas moleculares similares

Escenarios Aplicables

  1. Física solar: Diagnóstico de campos magnéticos de la fotosfera y modelado de líneas de C2_2
  2. Atmósferas estelares: Análisis de líneas moleculares de otras estrellas
  3. Investigación teórica: Física de colisiones moleculares y teoría de polarización

Referencias

El artículo cita referencias abundantes y relevantes, incluyendo:

  • Gandorfer (2000): Observaciones del segundo espectro solar
  • Berdyugina & Fluri (2004): Teoría de diagnóstico de campos magnéticos de C2_2
  • Corey & Alexander (1985): Marco teórico de colisiones moleculares
  • Qutub et al. (2020, 2021): Investigaciones anteriores sobre MgH y CN

Esta investigación proporciona una base teórica importante para la espectroscopia molecular solar, con valor particular en aplicaciones de diagnóstico de campos magnéticos. Mediante cálculos cuánticos precisos y métodos innovadores de procesamiento de datos, realiza contribuciones importantes a la comprensión de procesos físicos en la capa fotosférica solar.