The HST/WFC3 Transmission Spectrum of AU Mic b Part I: An Atmosphere Obscured by Contamination and Systematics

Resumen Técnico: El Espectro de Transmisión de HST/WFC3 de AU Mic b Parte I

Planteamiento del Problema
Los progenitores de sub-Neptunos jóvenes que orbitan enanas M ofrecen una ventana crítica hacia la formación y evolución temprana de planetas, revelando potencialmente las condiciones prístinas antes de que ocurra una pérdida atmosférica significativa. AU Mic b, un sub-Neptuno de $\sim$20 Ma que orbita una estrella de tipo M0Ve de pre-secuencia principal cercana (9.7 pc), es un candidato ideal para este estudio. Sin embargo, la caracterización de las atmósferas de planetas alrededor de enanas M jóvenes y activas se complica debido a la intensa variabilidad estelar. Específicamente, el efecto de la "Fuente de Luz de Tránsito" (TLS, por sus siglas en inglés)—donde las manchas estelares y las fáculas no ocultas alteran la profundidad aparente del tránsito—puede imitar u ocultar las características atmosféricas. Además, las observaciones del Telescopio Espacial Hubble (HST) de AU Mic b se vieron afectadas por un escaneo inestable del telescopio, lo que introdujo una Función de Propagación de Punto (PSF) variable que generó sistemáticas significativas de órbita a órbita. El desafío principal abordado en este trabajo es desenredar el espectro de transmisión planetaria de estos contaminantes instrumentales y estelares combinados para determinar si las características atmosféricas pueden ser restringidas.

Metodología
Los autores analizaron dos tránsitos de HST/WFC3 de AU Mic b observados en agosto de 2021 (grisma G141, 1.14–1.64 $\mu$m) y abril de 2022 (grisma G102, 0.79–1.13 $\mu$m). La reducción y el análisis de los datos siguieron un enfoque de múltiples etapas:

1. Reducción de Datos y Agrupación (Binning): Utilizando el pipeline PACMAN, los autores redujeron los datos brutos. Para mitigar la PSF variable causada por el escaneo inestable, emplearon un esquema de agrupación de longitud de onda personalizado. En lugar de grupos uniformes, colocaron los bordes de los grupos basándose en el RMS residual de los ajustes lineales de las curvas de luz fuera de tránsito (OOT) y el desfase de flujo normalizado entre órbitas. Esto redujo la resolución espectral pero minimizó el impacto de la PSF variable en la línea base.
2. Caracterización Estelar (SED OOT): Los autores modelaron la distribución espectral de energía (SED) de la estrella fuera de tránsito utilizando un modelo de descomposición espectral multitemperatura (biblioteca BT-Settl AGSS2009). Probaron modelos de 1, 2 y 3 temperaturas para restringir la temperatura de la fotósfera ($T_{phot}$), la temperatura de las manchas ($T_{spot}$) y el factor de llenado de manchas ($f_{spot}$).
3. Modelado de Curvas de Luz: Las curvas de luz de luz blanca se modelaron utilizando dos enfoques: un modelo "con manchas" (usando el código fleck para contabilizar la modulación rotacional y las ocultaciones de manchas) y un modelo de línea base lineal. El modelo con manchas identificó cruces de manchas específicos en la cuerda del tránsito, los cuales fueron recortados de las curvas de luz espectroscópicas. Las sistemáticas instrumentales (efectos de rampa y de respiración) se modelaron mediante una combinación de rampas exponenciales y polinomios de 4º orden en la fase orbital del HST.
4. Recuperación Atmosférica (Retrieval): Se realizaron recuperaciones bayesianas en el espectro de transmisión resultante utilizando el código PLATON. Se compararon tres modelos:
   • Modelo Completo (Full Model): Incluye parámetros atmosféricos (metalicidad, C/O, temperatura, nubes, brumas) y el efecto TLS (parametrizado por $T_{phot}$, $T_{spot}$ y $f_{spot}$).
   • Solo TLS (TLS-Only): Asume un espectro atmosférico plano donde todas las características se atribuyen al efecto TLS.
   • Función de Escalón (Step Function): Un espectro sin características con solo un desfase de conjunto de datos.

Resultados Clave

• Parámetros Estelares: El modelo de 3 temperaturas fue fuertmente preferido para la SED estelar. Los autores restringieron la temperatura de la fotósfera a $3891 \pm 37$ K y la temperatura de las manchas a $3020 \pm 69$ K, con un factor de llenado de manchas de $0.33 \pm 0.05$. Se detectó un tercer componente más caliente ($\sim$6000 K), pero probablemente compensa sistemáticas no resueltas en lugar de representar una mancha caliente física.
• Sistemáticas Instrumentales: El escaneo inestable creó una PSF variable que se manifestó como desplazamientos de la línea base de órbita a órbita. Aunque la agrupación mitigó esto, los efectos residuales siguieron siendo la principal limitación en la precisión.
• Espectro de Transmisión: El espectro recuperado está dominado por el efecto TLS. Los datos prefieren fuertemente los casos de "Modelo Completo" y "Solo TLS" sobre la "Función de Escalón" (evidencia decisiva), pero los datos no pueden distinguir entre los casos de "Modelo Completo" y "Solo TLS". Los parámetros atmosféricos recuperados (masa, metalicidad, temperatura) están pobremente restringidos.
• Escala de Altura Atmosférica: A pesar de las débiles restricciones, la recuperación establece un límite superior estricto en la escala de altura atmosférica de $<185$ km (a $3\sigma$). Esto es aproximadamente 1/5 de la escala de altura de los sub-Neptunos jóvenes que orbitan estrellas tipo G (por ejemplo, V1298 Tau b, HIP 67522 b).
• Consistencia: El factor de llenado de manchas derivado de la recuperación atmosférica ($0.36 \pm 0.06$ para el Modelo Completo) es consistente con el valor derivado independientemente de la SED OOT ($0.33 \pm 0.05$), lo que sugiere que el modelado del TLS es consistente consigo mismo.

Significancia y Reivindicaciones
El artículo concluye que el espectro de transmisión de AU Mic b está actualmente dominado por la contaminación estelar (el efecto TLS) y las sistemáticas del HST, lo que hace que las características atmosféricas estén débilmente restringidas. Los autores afirman modestamente que los datos son insuficientes para detectar características moleculares como H$_2$O o CH$_4$ con alta confianza.

La significancia de este trabajo radica en:

1. Demostrar la Dominancia del TLS: Proporciona un caso de estudio claro donde el efecto TLS sobrepasa las señales atmosféricas en un sistema de enana M joven, advirtiendo contra la sobreinterpretación de pendientes o características espectrales en sistemas similares sin un modelado estelar robusto.
2. Caracterizar la Naturaleza de AU Mic b: El límite superior derivado de la escala de altura sugiere que AU Mic b es distinto de los planetas jóvenes de atmósfera clara y baja densidad que orbitan estrellas G. Probablemente posee una masa más alta, una metalicidad más alta o una capa gruesa de nubes/brumas, lo que implica vías de formación o evolución diferentes para los sub-Neptunos alrededor de enanas M.
3. Marco Metodológico: El artículo establece un marco riguroso para manejar escaneos inestables del HST y actividad estelar, preparando el terreno para el "Papel II", que presentará un modelo empírico para decontaminar el espectro.
4. Predicciones para JWST: Los autores predicen que las observaciones de JWST (NIRCam F322W2 y F444W) probablemente también estarán dominadas por el efecto TLS, con solo contribuciones atmosféricas menores, y que las profundidades de tránsito ópticas podrían ser $\sim$25% más altas de lo medido en el infrarrojo cercano debido al efecto TLS.

El trabajo subraya los desafíos en la caracterización de sistemas planetarios jóvenes y activos, y motiva el desarrollo de mejores modelos empíricos para las fotósferas estelares para desbloquear las propiedades atmosféricas de estos mundos primordiales.