Native-resolution retrievals of VHS 1256-1257 b spanning the JWST/NIRSpec wavelength range: Chemical composition of a partially cloudy atmosphere

Resumen Técnico: Recuperaciones de resolución nativa de VHS 1256-1257 b que abarcan el rango de longitud de onda de JWST/NIRSpec

Problema y Contexto
El compañero de masa planetaria (PMC) VHS 1256-1257 b es uno de los objetos subestelares más variables conocidos, exhibiendo variaciones de flujo del 10–30%. Su historia de formación sigue siendo objeto de debate; si bien su amplia separación y el eje de rotación desalineado sugieren una formación "top-down" (de arriba hacia abajo) mediante fragmentación de núcleo o filamento, su masa se sitúa en el límite de la quema de deuterio (~13 $M_{Jup}$), lo que genera ambigüedad sobre si es un planeta gigante o una enana marrón. Estudios previos utilizando datos de JWST/NIRSpec (p. ej., Miles et al. 2023; Gandhi et al. 2023) proporcionaron información, pero estuvieron limitados por el agrupamiento de datos (reduciendo la resolución a $R \sim 300$) o por una cobertura de longitud de onda restringida (por ejemplo, solo la rejilla G395H). Estas limitaciones dificultaron la caracterización precisa del desequilibrio químico, las proporciones isotópicas y el origen de la extrema variabilidad del objeto.

Metodología
Los autores presentan una recuperación atmosférica integral de VHS 1256-1257 b utilizando todo el rango de longitud de onda de 0.97–5.27 $\mu$m cubierto por las rejillas JWST/NIRSpec G140H, G235H y G395H en sus resoluciones espectrales nativas ($R \sim 2000–4100$).

1. Reducción de Datos: El equipo re-redujo los datos del programa ERS (Early Release Science) #1386 utilizando el último pipeline de JWST (v1.20.2) y el contexto CRDS (1464). Esta nueva reducción corrige las oscilaciones espurias del continuo y los valores atípicos presentes en la reducción inicial de ERS (Miles et al. 2023) y aplica correcciones de apertura para contabilizar el flujo perdido fuera del radio de extracción.
2. Marco de Recuperación (Retrieval): Utilizando el código de transferencia radiativa petitRADTRANS junto con el algoritmo de muestreo anidado PyMultiNest, los autores realizaron recuperaciones bayesianas.
   • Estructura Atmosférica: Modelaron una atmósfera parcialmente nubosa utilizando un enfoque de dos columnas (una nubosa, una más clara) para dar cuenta de la alta variabilidad del objeto. Las columnas están ponderadas por una fracción de cobertura ($CF$).
   • Nubes: Las nubes se modelaron con opacidades simplificadas no grises para evitar sesgar la solución, permitiendo una nube de hierro profunda y una cubierta de silicatos irregular (patchy).
   • Química: Se empleó un modelo de desequilibrio químico, calculando los puntos de extinción (quench points) para las especies principales (CO, CH$_4$, H$_2$O, NH$_3$, HCN, CO$_2$) basadas en las escalas de tiempo de mezcla vertical ($K_{zz}$). El modelo restringe directamente las abundancias elementales ([X/H]) y las proporciones isotópicas ($^{12}$C/$^{13}$C, $^{16}$O/$^{18}$O, $^{16}$O/$^{17}$O, H/D).
   • Parámetros: La recuperación ajusta simultáneamente la masa, el radio, la velocidad radial, los gradientes de temperatura, las propiedades de las nubes y las abundancias de C, O, N, S, F, K+Na, y Cr+Fe, así como las proporciones isotópicas.

Resultados Clave

• Ajuste del Modelo: El modelo de dos columnas proporciona un ajuste excelente a las observaciones, reduciendo los residuales a $\sim 1\%$. El modelo confirma que las longitudes de onda más cortas están dominadas por una columna atmosféica más clara (que cubre $\sim 21\%$ de la superficie), mientras que las longitudes de onda más largas están dominadas por una columna nubosa ($\sim 79\%$).
• Mecanismo de Variabilidad: La recuperación sugiere que cambios pequeños en la fracción de cobertura de nubes de $\pm 1–3\%$ son suficientes para reproducir la variabilidad fotométrica observada de 10–30%, apoyando la hipótesis de que la variabilidad surge de la rotación de estructuras de nubes no uniformes.
• Composición Química:
  • Metalicidad y Proporciones: La atmósfera exhibe una metalicidad similar a la solar. La relación C/O está restringida a $0.567^{+0.001}_{-0.001}$, y la relación C/S a $20.1^{+0.3}_{-0.3}$, ambas consistentes con los valores solares.
  • Isótopos: La relación $^{12}$C/$^{13}$C es de $89^{+2}_{-2}$, consistente con el viento solar y el medio interestelar (ISM) local. Sin embargo, la relación $^{16}$O/$^{18}$O es de $\sim 800$, indicando una deplesión significativa de $^{18}$O en comparación con el Sol y el ISM. La relación $^{16}$O/$^{17}$O está restringida tentativamente en $\sim 3000$.
  • Deuterio: No se encontró una detección significativa de especies deuteradas (HDO, CH$_3$D). La relación D/H recuperada es anomalamente alta ($\sim 16 \times 10^{-5}$), pero los autores advierten que esto es probablemente un sesgo debido a las no detecciones y las degeneraciones de parámetros, más que evidencia de un régimen de masa específico.
• Propiedades Globales: La recuperación arroja una masa de $12.9^{+0.1}_{-0.1} M_{Jup}$ y un radio de $1.433^{+0.002}_{-0.002} R_{Jup}$, con una temperatura efectiva de $1137^{+1}_{-1}$ K. Los autores señalan que estos valores están fuertemente correlacionados con otros parámetros y pueden estar sesgados.
• Mezcla Vertical: Los datos indican una mezcla vertical vigorosa en la atmósfera profunda ($K_{zz, lower} \approx 1.6 \times 10^9$ cm$^2$ s$^{-1}$) y una transición a un régimen radiativo con una eficiencia de mezcla reducida en la parte alta de la atmósfera ($K_{zz, upper} \approx 2.8 \times 10^4$ cm$^2$ s$^{-1}$).

Significado y Reivindicaciones
El artículo afirma que este estudio demuestra el valor crítico de las recuperaciones pancráticas de resolución nativa para caracterizar atmósferas extrasolares complejas. Al utilizar el rango espectral completo y una reducción de datos actualizada, los autores resuelven las degeneraciones que afectaron a estudios previos de menor resolución o de cobertura de longitud de onda parcial.

Específicamente, los autores destacan:

1. Impacto de la Reducción de Datos: La extracción espectral actualizada altera significativamente la composición química inferida. Estudios previos (p. ej., Gandhi et al. 2023) sugerían abundancias isotópicas super-solares ($^{13}$C, $^{18}$O, $^{17}$O), mientras que este estudio encuentra $^{13}$C de tipo solar y $^{18}$O deplesionado, atribuyendo las discrepancias previas a problemas en la reducción de datos inicial.
2. Restricciones de Formación: La deplesión de $^{18}$O inferida desafía los modelos simples de formación "top-down" donde el objeto heredaría la composición del ISM. Los autores sugieren que la fotodisociación selectiva de isótopos en la nube de formación estelar podría explicar la deplesión, pero reconocen que sin las mediciones de abundancia de la binaria de la estrella enana M anfitriona, el contexto químico permanece incompleto.
3. Avance Metodológico: Este trabajo sirve como prueba de concepto para las recuperaciones de alta resolución con múltiples rejillas, mostrando que tales análisis pueden restringir simultáneamente la cobertura de nubes, la estructura de temperatura y la química detallada, incluyendo las proporciones isotópicas, que son trazadores clave de los entornos de formación.

Los autores mantienen la modestia respecto a la determinación de la masa y la interpretación de la deplesión de $^{18}$O, señalando que las degeneraciones de parámetros (particularmente entre masa, radio y cobertura de nubes) y la falta de datos de abundancia de la estrella anfitriona limitan las conclusiones definitivas sobre la historia de formación y la clasificación exacta del objeto.