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

Resumo Técnico: Recuperações de resolução nativa de VHS 1256-1257 b abrangendo a faixa de comprimento de onda do JWST/NIRSpec

Problema e Contexto
O companheiro de massa planetária (PMC) VHS 1256-1257 b é um dos objetos subestelares mais variáveis conhecidos, exibindo variações de fluxo de 10–30%. Sua história de formação permanece debatida; embora seu amplo afastamento e eixo de rotação desalinhado sugiram uma formação "top-down" via fragmentação de núcleo ou filamento, sua massa situa-se na fronteira do limite de queima de deutério (~13 $M_{Jup}$), tornando ambíguo se é um planeta gigante ou uma anã marrom. Estudos anteriores utilizando dados do JWST/NIRSpec (ex: Miles et al. 2023; Gandhi et al. 2023) forneceram insights, mas foram limitados pelo agrupamento de dados (reduzindo a resolução para $R \sim 300$) ou pela cobertura de comprimento de onda restrita (ex: apenas a grade G395H). Essas limitações dificultaram a caracterização precisa do desequilíbrio químico, das razões isotópicas e da origem da extrema variabilidade do objeto.

Metodologia
Os autores apresentam uma recuperação atmosférica abrangente de VHS 1256-1257 b utilizando toda a faixa de comprimento de onda de 0,97–5,27 $\mu$m coberta pelas grades G140H, G235H e G395H do JWST/NIRSpec em suas resolões espectrais nativas ($R \sim 2000–4100$).

1. Redução de Dados: A equipe reprocessou os dados do programa de Ciência de Liberação Precoce (ERS) #1386 usando o pipeline mais recente do JWST (v1.20.2) e o contexto CRDS (1464). Esta nova redução corrige oscilações espúrias do continuum e valores atípicos presentes na redução inicial do ERS (Miles et al. 2023) e aplica correções de abertura para contabilizar o fluxo perdido fora do raio de extração.
2. Estrutura de Recuperação (Retrieval): Utilizando o código de transferência radiativa petitRADTRANS acoplado ao algoritmo de amostragem aninhada PyMultiNest, os autores realizaram recuperações bayesianas.
   • Estrutura Atmosférica: Eles modelaram uma atmosfera parcialmente nublada usando uma abordagem de duas colunas (uma nublada, uma mais clara) para considerar a alta variabilidade do objeto. As colunas são ponderadas por uma fração de cobertura ($CF$).
   • Nuvens: As nuvens foram modeladas com opacidades não-cinzas simplificadas para evitar o viés da solução, permitindo uma nuvem de ferro profunda e uma cobertura de silicato irregular.
   • Química: Um modelo de desequilíbrio químico foi empregado, calculando pontos de quenching para espécies principais (CO, CH$_4$, H$_2$O, NH$_3$, HCN, CO$_2$) baseados em escalas de tempo de mistura vertical ($K_{zz}$). O modelo constrange diretamente abundâncias elementares ([X/H]) e razões isotópicas ($^{12}$C/$^{13}$C, $^{16}$O/$^{18}$O, $^{16}$O/$^{17}$O, H/D).
   • Parâmetros: A recuperação ajusta simultaneamente massa, raio, velocidade radial, gradientes de temperatura, propriedades de nuvens e as abundâncias de C, O, N, S, F, K+Na, e Cr+Fe, bem como razões isotópicas.

Principais Resultados

• Ajuste do Modelo: O modelo de duas colunas fornece um ajuste excelente às observações, reduzindo os resíduos para $\sim 1\%$. O modelo confirma que os comprimentos de onda mais curtos são dominados por uma coluna atmosférica mais clara (cobrindo $\sim 21\%$ da superfície), enquanto os comprimentos de onda mais longos são dominados por uma coluna nublada ($\sim 79\%$).
• Mecanismo de Variabilidade: A recuperação sugere que pequenas mudanças na fração de cobertura de nuvens de $\pm 1–3\%$ são suficientes para reproduzir a observada variabilidade fotométrica de 10–30%, apoiando a hipótese de que a variabilidade surge da rotação de estruturas de nuvens não uniformes.
• Composição Química:
  • Metalicidade e Razões: A atmosfera exibe uma metalicidade semelhante à solar. A razão C/O é restringida a $0,567^{+0,001}_{-0,001}$, e a razão C/S a $20,1^{+0,3}_{-0,3}$, ambas consistentes com valores solares.
  • Isótopos: A razão $^{12}$C/$^{13}$C é de $89^{+2}_{-2}$, consistente com o vento solar e o meio interestelar (ISM) local. No entanto, a razão $^{16}$O/$^{18}$O é de $\sim 800$, indicando uma depleção significativa de $^{18}$O em relação ao Sol e ao ISM. A razão $^{16}$O/$^{17}$O é tentativamente restringida em $\sim 3000$.
  • Deutério: Não foi detectada nenhuma espécie deuterada significativa (HDO, CH$_3$D). A razão D/H recuperada é anormalmente alta ($\sim 16 \times 10^{-5}$), mas os autores alertam que isso é provavelmente um viés devido a não-detecções e degenerações de parâmetros, em vez de evidência para um regime de massa específico.
• Propriedades Globais: A recuperação produz uma massa de $12,9^{+0,1}_{-0,1} M_{Jup}$ e um raio de $1,433^{+0,002}_{-0,002} R_{Jup}$, com uma temperatura efetiva de $1137^{+1}_{-1}$ K. Os autores observam que esses valores estão fortemente correlacionados com outros parâmetros e podem ser enviesados.
• Mistura Vertical: Os dados indicam uma mistura vertical vigorosa na atmosfera profunda ($K_{zz, lower} \approx 1,6 \times 10^9$ cm$^2$ s$^{-1}$) e uma transição para um regime radiativo com eficiência de mistura reduzida em partes mais altas da atmosfera ($K_{zz, upper} \approx 2,8 \times 10^4$ cm$^2$ s$^{-1}$).

Significância e Alegações
O artigo afirma que este estudo demonstra o valor crítico de recuperações de resolução nativa pancromáticas para caracterizar atmosferas extrasolares complexas. Ao utilizar todo o intervalo espectral e uma redução de dados atualizada, os autores resolvem degenerações que afetaram estudos anteriores de menor resolução ou de cobertura parcial de comprimento de onda.

Especificamente, os autores destacam:

1. Impacto da Redução de Dados: A extração espectral atualizada altera significativamente a composição química inferida. Estudos anteriores (ex: Gandhi et al. 2023) sugeriam abundâncias isotópicas super-solares ($^{13}$C, $^{18}$O, $^{17}$O), enquanto este estudo encontra $^{13}$C solar e $^{18}$O depletado, atribuindo as discrepâncias anteriores a problemas na redução de dados inicial.
2. Restrições de Formação: A $^{18}$O depletada inferida desafia modelos simples de formação "top-down", onde o objeto herdaria a composição do ISM. Os autores sugerem que a fotodissoção seletiva de isótopos na nuvem de formação estelar poderia explicar a depleção, mas reconhecem que, sem medições da abundância da binária anã M hospedeira, o contexto químico permanece incompleto.
3. Avanço Metodológico: Este trabalho serve como uma prova de conceito para recuperações de alta resolução e múltiplas grades, mostrando que tais análises podem restringir simultaneamente a cobertura de nuvens, a estrutura de temperatura e a química detalhada, incluindo razões isotópicas, que são traçadores fundamentais de ambientes de formação.

Os autores permanecem modestos quanto à determinação da massa e à interpretação da depleção de $^{18}$O, observando que as degenerações de parâmetros (particularmente entre massa, raio e cobertura de nuvens) e a falta de dados de abundância da estrela hospedeira limitam conclusões definitivas sobre a história de formação e a classificação exata do objeto.