Infrared spectroscopy of gas-phase hydrogenated and methylated pyrenes: from laboratory spectra to the simulated 3.4 $\mu$m emission band

Problema
Las Bandas Infrarrojas Aromáticas (AIBs) son características de emisión prominentes en el medio interestelar, atribuidas principalmente a Hidrocarburos Aromáticos Policíclicos (PAHs). Si bien la banda de 3,3 µm está bien establecida como la firma del estiramiento C-H aromático, el rango de 3,4–3,6 µm contiene características satélite (notablemente en 3,40, 3,46, 3,51 y 3,57 µm) y una meseta ancha. Estas se asignan generalmente a vibraciones de estiramiento C-H alifático en PAHs hidrogenados (H-PAHs) y metilados (Me-PAHs). Sin embargo, los portadores moleculares específicos de estas características, particularmente el componente prominente de 3,403 µm, siguen siendo objeto de debate. Además, los datos teóricos y experimentales existentes a menudo carecen de las correcciones de anharmonicidad necesarias para simular con precisión los espectros de emisión de PAHs calientes mientras se enfrían tras la excitación por UV, lo que limita la comparación directa con las observaciones de alta resolución del Telescopio Espacial James Webb (JWST).

Metodología
Los autores llevaron a cabo un estudio exhaustivo combinando espectroscopia de laboratorio de alta temperatura, cálculos teóricos y modelado astrofísico:

1. Espectroscopia Experimental: Se registraron espectros de absorción en el infrarrojo medio (1,4–25 µm) para pireno, pireno metilado (1-metilpireno) y pirenos hidrogenados (dihidro-, tetrahidro- y hexahidropireno) en fase gaseosa a temperaturas que oscilan entre 373 K y 673 K. Las mediciones se realizaron utilizando una celda de alta temperatura dedicada en el montaje ESPOIRS en IRAP.
2. Cálculos Teóricos: Se realizaron cálculos de la Teoría del Funcional de la Densidad (DFT) utilizando el funcional B3LYP-D3 con el conjunto de bases N-07D para obtener campos de fuerza semicuárticos. El código AnharmoniCaOs se utilizó para calcular espectros de absorción anharmónicos a 0 K para el estado vibracional fundamental, apoyando la asignación de las bandas experimentales.
3. Extracción de Anarmonicidad: Se utilizó una herramienta interactiva (cosmicPAHmfit) para realizar ajustes de múltiples componentes de las bandas experimentales. Esto permitió la derivación de leyes de anharmonicidad empíricas que describen la evolución de las posiciones y anchuras de las bandas en función de la temperatura. El ensanchamiento rotacional se estimó y se restó para aislar los efectos vibracionales.
4. Modelado de Emisión: Se empleó un código cinético de Monte Carlo para simular las cascadas de enfriamiento de los PAHs tras la absorción de fotones UV (hasta el límite de disociación, ~6 eV). El modelo incorporó los parámetros de anharmonicidad empírica derivados para generar espectros de emisión sintéticos que tengan en cuenta los desplazamientos y el ensanchamiento de las bandas dependientes de la temperatura.
5. Comparación Astrofísica: Los espectros sintéticos resultantes se compararon con observaciones del JWST de la Región de Barrera de Orión (programa PDRs4All) y datos terrestres de NGC 1333.

Contribuciones Clave

• Nuevos Datos Experimentales: Proporcionó espectros de infrarrojo en fase gaseosa de última generación para pirenos hidrogenados y metilados a temperaturas astrofísicas relevantes, cubriendo la región crítica de 3,3–3,6 µm.
• Leyes de Anharmonicidad Empírica: Derivó parámetros dependientes de la temperatura para las posiciones y anchuras de las bandas para estas especies que contienen grupos alifáticos, abordando un vacío en las bases de datos existentes que carecían de datos para PAHs con enlaces C-H alifáticos.
• Integración de Bases de Datos: Todos los espectros experimentales, parámetros de anharmonicidad y cálculos teóricos se han puesto a disposición pública a través de CosmicPAH-IRDB y la base de datos espectral teórica Cagliari-Toulouse.
• Identificación Refinada de Portadores: Propuso un portador molecular específico para la característica de emisión de 3,403 µm basándose en una comparación rigurosa de espectros simulados y observados.

Resultados

• Características Espectrales: El estudio confirmó que la hidrogenación altera significativamente la región de estiramiento C-H. Mientras que el pireno metilado contribuye principalmente a la banda de 3,3 µm y a una meseta débil, las especies fuertemente hidrogenadas (específicamente el 6H-pireno) dominan la región de 3,4 µm.
• Efectos de Temperatura: A medida que la temperatura aumenta, las bandas alifáticas generalmente sufren un desplazamiento hacia el rojo (redshift) y se ensanchan. El estudio observó que el 2H- y el 4H-pireno son térmicamente inestables a temperaturas más altas, regresando a pireno, mientras que el 6H-pireno exhibe una mayor resistencia térmica.
• Identificación del Portador: El espectro de emisión sintética del 1,2,3,6,7,8-hexahidropireno (6H-pireno), excitado por fotones UV de ~4 eV, proporciona una coincidencia detallada con el componente de 3,403 µm observado en la Barrera de Orión. El modelo reproduce con éxito el perfil y la posición de la banda sin requerir parámetros no físicos.
• Pireno Metilado: Aunque el 1-metilpireno contribuye a la meseta de emisión subyacente y a la banda de 3,3 µm, su falta de una banda de infrarrojo fuerte y distinta en el rango de 3,4–3,6 µm complica su detección como un portador específico para las características satélite en comparación con el 6H-pireno.
• Precisión del Modelo: La inclusión de funciones de anharmonicidad empírica mejoró significativamente la concordancia entre los espectros simulados y observados en comparación con los modelos que utilizan únicamente espectros de laboratorio medidos a una sola temperatura (por ejemplo, 523 K).

Significancia
El artículo afirma que la identificación del 6H-pireno como el portador de la banda de 3,403 µm redefine la interpretación de la característica de 3,4 µm en observaciones astronómicas. Sugiere que el componente de 3,403 µm es un trazador robusto del contenido alifático en regiones blindadas donde los enlaces alifáticos frágiles pueden sobrevivir a una excitación UV suave (~4 eV). En consecuencia, los autores argumentan que utilizar la relación de bandas 3,4-a-3,3 µm como un trazador simple del grado general de superhidrogenación en las poblaciones de PAHs es una aproximación, ya que el componente de 3,395 µm (y la banda de 3,3 µm) puede originarse de diferentes portadores o condiciones de excitación. Este trabajo proporciona un marco más preciso para interpretar los datos del JWST en diversos entornos, desde regiones de fotodisociación hasta discos protoplanetarios, al vincular estructuras moleculares específicas con características espectrales observadas mediante restricciones teóricas y de laboratorio rigurosas.