A parameterised approach to disequilibrium retrievals in the JWST era: Application to NIRCam observations of HD 189733b

Resumen Técnico: Un Enfoque Parametrizado para las Recuperaciones de Desequilibrio en la Era del JWST

Planteamiento del Problema
Las recuperaciones atmosféricas (retrievals) son esenciales para inferir las propiedades físicas y químicas de los exoplanetas a partir de espectros observados. Sin embargo, un supuesto simplificador común en estos análisis es que la atmósfera existe en equilibrio termoquímico. Si bien este supuesto se sostiene para los Júpiter calientes altamente irradiados, a menudo falla para objetos más fríos o aquellos con una fuerte mezcla vertical, donde el transporte induce el enfriamiento (quenching) que lleva las abundancias químicas fuera del equilibrio. Asumir el equilibrio en presencia de desequilibrio puede conducir a sesgos significativos en los parámetros atmosféricos globales recuperados, específicamente en la relación carbono-oxígeno (C/O) y la metalicidad ([M/H]). Los marcos de recuperación existentes enfrentan un compromiso: los modelos de "química libre" son flexibles pero propensos al sobreajuste (overfitting); los modelos de equilibrio son eficientes pero físicamente inexactos para escenarios de mezcla; y los modelos completos de radiación-convección-desequilibrio son físicamente consistentes pero computacionalmente costosos y dependientes de rejillas (grids) precomputadas.

Metodología
Este estudio presenta un marco de recuperación parametrizado que cierra la brecha entre la eficiencia y el realismo físico mediante el modelado explícito del enfriamiento inducido por el transporte. Los autores utilizan la herramienta de transferencia radiativa y recuperación NemesisPy, acoplada con FastChem 3.0 para la química de equilibrio y VULCAN para generar perfiles de abundancia de desequilibrio sintéticos.

1. Generación de Datos Sintéticos: Los autores generaron espectros de transmisión sintéticos del Telescopio Espacial James Webb (JWST) para el Júpiter caliente HD 189733 b. Estos espectros fueron producidos utilizando VULCAN para simular perfiles de abundancia química bajo variaciones en la fuerza de la mezcla vertical, parametrizada por el coeficiente de difusión de eddy ($K_{zz} = 10^6, 10^8, 10^{10}, 10^{12} \text{ cm}^2\text{s}^{-1}$).
2. Parametrización del Desequilibrio: Para incorporar efectos de desequilibrio sin el costo computacional de las rejillas cinéticas completas, los autores introdujeron parámetros libres para las presiones de enfriamiento (quench pressures). Se probaron dos enfoques:
   • Enfriamiento Elemental: Agrupar especies por constituyentes elementales (por ejemplo, enfriar todas las especies que contienen carbono en una única presión $P_{q,C}$ y las especies que contienen nitrógeno en una presión $P_{q,N}$).
   • Enfriamiento Molecular: Tratar a las moléculas individuales (por ejemplo, $\text{CH}_4$, $\text{NH}_3$) como si tuvieran presiones de enfriamiento independientes.
     La abundancia de una especie $X$ se define como el valor de equilibrio por debajo de la presión de enfriamiento ($P > P_{\text{quench}}$) y constante en el valor de equilibrio calculado en $P_{\text{quench}}$ por encima de ella ($P \le P_{\text{quench}}$).
3. Simulaciones de Recuperación: El marco fue probado contra observaciones sintéticas de los instrumentos JWST NIRISS/SOSS y NIRSpec/G395H. Las recuperaciones se realizaron asumiendo: (i) equilibrio termoquímico, (ii) enfriamiento elemental, y (iii) enfriamiento molecular.
4. Aplicación a Datos Reales: El marco se aplicó a las observaciones de transmisión publicadas de JWST/NIRCam para HD 189733 b (Fu et al. 2024). Se desarrolló una novedosa parametrización de $\text{H}_2\text{S}$ para modelar la región fotoquímicamente activa, definiendo una abundancia profunda y una presión de ruptura ($P_{\text{break}}$) por encima de la cual la abundancia declina mediante una ley de potencia.

Resultados Clave

• Sesgo en los Supuestos de Equilibrio: Las pruebas de recuperación sintética demostraron que asumir el equilibrio termoquímico cuando hay mezcla vertical conduce a sesgos significativos. Específicamente, las recuperaciones de equilibrio a menudo subestimaron los valores de C/O y [M/H], fallando en recuperar los parámetros de entrada reales a medida que aumentaba la fuerza de la mezcla.
• Recuperación de Propiedades Globales: Los marcos de recuperación de desequilibrio (tanto elemental como molecular) recuperaron con éxito las propiedades atmosféricas globales correctas (C/O y [M/H]) en la mayoría de los casos sintéticos, incluso cuando había una fuerte mezcla vertical. La inclusión de solo dos parámetros libres adicionales (presiones de enfriamiento) fue suficiente para corregir los sesgos inherentes a los modelos de equilibrio.
• Restricciones de la Presión de Enfriamiento: El marco recuperó con éxito las presiones de enfriamiento aproximadas para las especies de carbono y nitrógeno. Los resultados indicaron que la química del nitrógeno se enfría en niveles atmosféricos más profundos que la química del carbono, lo cual es consistente con las expectativas teóricas respecto a las escalas de tiempo de las reacciones.
• Aplicación a HD 189733 b: Al aplicarse a los datos reales de JWST, el modelo de desequilibrio proporcionó un mejor ajuste (mayor evidencia bayesiana) que el modelo de equilibrio químico. El modelo de mejor ajuste sugirió una fuerte mezcla vertical con una presión de enfriamiento de carbono de $\log P_{q,C} \approx 1.7$ bar.
• $\text{H}_2\text{S}$ y Fotoquímica: El estudio proporcionó evidencia tentativa de la región fotoquímicamente activa de la atmósfera. Al parametrizar el perfil vertical de $\text{H}_2\text{S}$, los autores restringieron la presión de ruptura ($P_{\text{break}}$) a $\log P_{\text{break}} \approx -3.25$ bar. Esto representa la primera restricción tentativa del gradiente vertical de un gas atmosférico en la atmósfera de un Júpiter caliente, ofreciendo un método indirecto para inferir la ubicación donde el azufre atómico se vuelve dominante.

Significancia y Reivindicaciones
El artículo afirma ofrecer una alternativa computacionalmente eficiente y físicamente motivada tanto para los métodos de recuperación de química libre como para los de rejilla completa. Al parametrizar las presiones de enfriamiento, los autores demuestran que es posible:

1. Inferir correctamente las propiedades atmosféricas globales (C/O y metalicidad) incluso en presencia de una mezcla vertical significativa, evitando los sesgos de los supuestos de equilibrio.
2. Sondear la estructura térmica profunda y las fuerzas de mezcla en las atmósferas de exoplanetas a través de la recuperación de las presiones de enfriamiento.
3. Restringir las regiones fotoquímicamente activas de las atmósferas de exoplanetas utilizando datos de JWST, específicamente a través de la novedosa parametrización de $\text{H}_2\text{S}$.

Los autores concluyen que los futuros estudios de recuperación de atmósferas de exoplanetas deberían considerar la inclusión de parámetros de presión de enfriamiento para especies moleculares clave para lograr una caracterización más físicamente precisa de la composición atmosférica y para dar cuenta de procesos de mezcla vertical desconocidos. Observan que, si bien sus modelos de desequilibrio y de equilibrio con $\text{H}_2\text{S}$ producen ajustes similares con los datos actuales, la degeneración entre ellos podría romperse con observaciones futuras de NIRISS/SOSS.