Em 11 horas de espectros do NIRSpec, WISE 0855 oscilou até 10% de brilho e revelou duas camadas de variabilidade: nuvens de água no topo e gás rico em monóxido de carbono e fosfina trazido de 70 bar de profundidade

A7,4 anos-luz do Sol flutua um mundo que nenhuma estrela ilumina. Ele tem o tamanho de Júpiter, a massa de alguns Júpiteres e uma temperatura média em torno de oito graus Celsius negativos, mais ou menos o que marca o congelador de uma geladeira simples, e ainda assim quente o bastante para que a água forme nuvens, como as que atravessam o céu da Terra. Chama-se WISE 0855. É a anã marrom mais fria que conhecemos, o quarto sistema mais próximo do Sol e, no sentido mais literal, uma estrela que não deu certo.

Da noite de 1º de dezembro até a manhã de 2 de dezembro de 2023, no horário de Brasília, o Telescópio Espacial James Webb apontou para esse ponto quase invisível no céu e ficou ali, imóvel, por onze horas seguidas, registrando um espectro a cada quinze minutos. O que saiu dessa vigília é a primeira série temporal espectroscópica de resolução média já obtida da anã marrom mais fria conhecida e, com ela, o boletim meteorológico mais detalhado de um mundo com nuvens de água fora do Sistema Solar.

Anãs marrons são o elo entre planetas e estrelas. Nascem como estrelas, do colapso de uma nuvem de gás, mas não acumulam massa suficiente, menos de 75 a 80 vezes a de Júpiter, para acender a fusão do hidrogênio em seu núcleo. Sem essa fornalha, passam a vida esfriando, e as mais antigas e menos massudas chegam a temperaturas que já não lembram nada de estelar. Os astrônomos as classificam em sequências espectrais que descem a escada da temperatura: L, T e, no degrau mais baixo, Y. As anãs Y, reveladas pelo satélite WISE (Wide-field Infrared Survey Explorer) a partir de 2011, têm temperaturas efetivas abaixo de cerca de 450 kelvins, ou 177 graus Celsius. São tão frias que brilham quase exclusivamente no infravermelho, onde só telescópios especializados conseguem enxergá-las, e suas atmosferas, dominadas por hidrogênio, exibem bandas profundas de metano e amônia.

WISE 0855, nome completo WISE J085510.83-071442.5, foi identificada em 2014 pelo astrônomo Kevin Luhman, da Universidade Estadual da Pensilvânia, que combinou imagens do WISE com observações do telescópio espacial Spitzer e notou um ponto que se deslocava depressa demais contra o fundo de estrelas, sinal inequívoco de proximidade. A distância, hoje medida em 2,28 parsecs, coloca o objeto logo atrás do sistema Alfa Centauri, da Estrela de Barnard e do par Luhman 16 na lista dos vizinhos mais próximos do Sol. Sua massa, estimada entre três e dez vezes a de Júpiter, fica abaixo do limite em que o deutério entra em fusão, o que a situa num território ambíguo: tem massa de planeta, mas formou-se e vive sozinha, como uma estrela. A temperatura efetiva adotada pelo novo estudo é de 265 kelvins.

É essa temperatura que faz de WISE 0855 um laboratório sem equivalente. Em Júpiter, cuja temperatura efetiva é de cerca de 124 kelvins, a amônia congela e forma a camada de nuvens mais alta, escondendo a água, que só condensa muito mais embaixo. Em WISE 0855 acontece o contrário: a atmosfera é fria o suficiente para que o vapor d’água condense nas camadas superiores, mas quente demais para a amônia virar gelo. A consequência é que as nuvens de água, as mesmas que cobrem a Terra, ficam expostas na face visível do objeto, sem nada acima que as esconda. Os modelos teóricos preveem que anãs Y abaixo de cerca de 375 kelvins formem essas nuvens entre 1,2 e 2,0 bar de pressão, o equivalente a uma ou duas atmosferas terrestres. Havia, porém, um problema antigo: ninguém conseguia provar que elas estavam lá.

A disputa sobre a existência dessas nuvens tem mais de uma década. Em 2014, a primeira detecção de WISE 0855 no infravermelho próximo, por Jacqueline Faherty e colegas, mostrou cores compatíveis com modelos que incluíam nuvens de água e de sulfetos. No mesmo ano, Luhman e Taran Esplin argumentaram que modelos sem nuvens, com uma química fora de equilíbrio, também davam conta dos dados. Fotometria mais ampla, do óptico ao infravermelho médio, publicada em 2016 por Adam Schneider e, separadamente, por Luhman e Esplin, concluiu que nenhuma grade de modelos reproduzia as cores do objeto nem as da sequência inteira das anãs Y. Os primeiros espectros, obtidos no mesmo ano pela equipe de Andrew Skemer na banda M, em torno de 4,8 micrômetros, eram mais bem ajustados quando se adicionava uma fonte cinza de opacidade exatamente onde as nuvens de água deveriam estar. Era um indício, não uma prova.

Também em 2016, Esplin e colaboradores usaram o Spitzer para acompanhar o brilho de WISE 0855 ao longo de horas e detectaram variações de alguns por cento, a primeira evidência de que a atmosfera não é uniforme, de que algo gira para dentro e para fora do campo de visão. Eles descartaram manchas quentes e sugeriram que uma mudança de poucos por cento na cobertura de nuvens bastaria para explicar o sinal, mas admitiram que faltavam modelos melhores para interpretar a física por trás dele. Em 2018, Caroline Morley e colegas apresentaram o espectro na banda L e mostraram que nuvens de água explicavam os desvios do espectro em comprimentos de onda acima de 1,6 micrômetro, enquanto uma fonte de opacidade cinza profunda, perto de 11 bar, dava conta das bandas ópticas e do infravermelho próximo.

A chegada do James Webb mudou a escala do debate sem resolvê-lo. Em 2024, Luhman e colaboradores publicaram o espectro de 1 a 5 micrômetros e concluíram que um modelo sem nuvens, com química fora de equilíbrio e um índice adiabático ajustado, reproduzia as grandes feições espectrais, mas deixava discrepâncias entre 2 e 3 micrômetros. Melanie Rowland e colegas, no mesmo ano, fizeram análises de recuperação atmosférica, ajustes livres que inferem composição e perfil de temperatura diretamente dos dados, e explicaram a região de 4 a 5 micrômetros apenas com abundâncias moleculares fora de equilíbrio, ao preço de um perfil de temperatura que se afastava muito do esperado para uma atmosfera em equilíbrio radiativo-convectivo. Esse desvio apontava para uma fonte de opacidade ausente, que poderia ser, de novo, nuvens. O mesmo trabalho mediu pela primeira vez, numa anã marrom, as abundâncias de fosfina e de metano deuterado.

O espectro mais amplo já obtido do objeto, de 5 a 22 micrômetros, veio em 2025 com o instrumento MIRI do Webb, em um estudo liderado por Kühnle. As recuperações atmosféricas indicaram que a água está esgotada acima de 3 bar, compatível com condensação, mas sem prová-la, e o melhor ajuste com nuvens colocou uma camada a 10 bar, bem abaixo dos 1,2 a 2,0 bar previstos para a condensação da água. Todos esses estudos tropeçaram no mesmo detalhe: os raios inferidos para WISE 0855 ficavam entre 0,79 e 0,89 raio de Júpiter, pequenos demais para um objeto dessa temperatura com idade entre um e dez bilhões de anos, cujo raio esperado é o de Júpiter. Nenhum modelo direto ajustava o espectro inteiro, e os modelos com nuvens e química fora de equilíbrio de Brianna Lacy e Adam Burrows eram os que chegavam mais perto.

O novo trabalho, liderado por Brittany Miles, do Observatório Steward da Universidade do Arizona, com 26 coautores de instituições dos Estados Unidos e da Irlanda, ataca o problema por um ângulo que nenhum dos anteriores pôde usar: o tempo. Em vez de um único espectro médio, a equipe obteve 44 espectros consecutivos, um a cada quinze minutos, ao longo de 11 horas, dentro do programa 2327 do Ciclo 1 do Webb, proposto por Skemer sob o título “Nuvens de gelo de água e tempo na anã marrom mais fria”. A observação começou às 01h03 UTC de 2 de dezembro de 2023 e terminou às 12h09 do mesmo dia. Como a anã gira, regiões diferentes da atmosfera entram e saem do hemisfério voltado para nós; se o espectro muda, e muda de maneira diferente em cada comprimento de onda, a mudança carrega a assinatura do que varia na superfície.

O instrumento escolhido foi o NIRSpec, o espectrógrafo de infravermelho próximo construído pela Agência Espacial Europeia, operado no modo Bright Object Time-Series, série temporal de objeto brilhante, que usa uma abertura quadrada de 1,6 segundo de arco de lado para capturar toda a luz do alvo e minimizar perdas causadas por pequenos deslocamentos do telescópio. Não houve reposicionamentos durante a sequência, justamente para preservar a estabilidade fotométrica. A rede de difração G395M com o filtro F290LP cobriu de 2,87 a 5,27 micrômetros com resolução espectral em torno de mil, o bastante para separar as bandas moleculares. Cada integração de quinze minutos acumulou 996 leituras sucessivas do detector, e a relação sinal-ruído por integração ficou entre 80 e 100 nos comprimentos de onda mais brilhantes, com média de 73 entre 3,83 e 5,17 micrômetros, um nível impensável antes do Webb para um objeto que emite tão pouco.

Figura 1 — Espectro médio de WISE 0855 ao longo das 11 horas, em escala logarítmica de fluxo, com as barras de erro de u
Figura 1 — Espectro médio de WISE 0855 ao longo das 11 horas, em escala logarítmica de fluxo, com as barras de erro de uma única integração de 15 minutos em cinza (acima), e relação sinal-ruído por integração (abaixo). Acima de 3,6 micrômetros o sinal-ruído se mantém acima de 5; a média entre 3,83 e 5,17 micrômetros é 73.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327

Transformar essas leituras em espectros confiáveis exigiu a cadeia padrão de redução do Webb, com correção de viés, de linearidade, remoção de raios cósmicos e calibração de fluxo, seguida de procedimentos próprios para extrair o traço espectral. A equipe ajustou uma gaussiana em cada coluna do detector para localizar o centro do espectro, extraiu a luz num raio de quatro pixels e subtraiu o fundo estimado a partir dos pixels vizinhos. A estabilidade foi notável: ao longo das onze horas, a posição do objeto no detector mudou cerca de um décimo de pixel, e a solução de comprimento de onda deslocou-se 900 vezes menos do que a largura de um único elemento espectral. Pixels defeituosos em cinco comprimentos de onda e pixels quentes em sete dos 44 espectros foram mascarados. Abaixo de 3,6 micrômetros o objeto é fraco demais e o ruído domina; a análise concentrou-se acima desse limite.

Dividindo cada espectro pela média dos 44, os pesquisadores construíram um mapa de variação percentual: comprimento de onda no eixo horizontal, tempo no vertical, cor indicando se a atmosfera estava mais brilhante ou mais apagada que a média naquele instante. O mapa revelou que o objeto como um todo oscilou entre 1% e 10% de brilho, mas a oscilação não foi uniforme. A feição mais variável de todas fica entre 4,5 e 4,8 micrômetros, exatamente onde o monóxido de carbono absorve luz: ali, a diferença entre o momento mais escuro e o mais claro chegou a 10% de pico a pico. Nessa faixa, a anã ficou mais apagada em torno da segunda hora de observação e mais brilhante perto da sétima hora e meia, em contraste com a curva de luz integrada em todo o espectro.

Figura 2 — Espectros de brilho mínimo (rosa) e máximo (roxo) na região de alto sinal-ruído, sobrepostos às curvas de opa
Figura 2 — Espectros de brilho mínimo (rosa) e máximo (roxo) na região de alto sinal-ruído, sobrepostos às curvas de opacidade de metano, amônia, CO₂, fosfina, CO e água (acima); razão entre o máximo e o mínimo (meio); e mapa da variação percentual do espectro ao longo do tempo (abaixo). A mancha escura em torno de 4,6 micrômetros nas primeiras horas marca o enfraquecimento das bandas do monóxido de carbono.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327

Em curvas de luz calculadas em faixas de 0,1 micrômetro, o padrão se repete com nuances. Entre 3,75 e 5,1 micrômetros há mudanças amplas e de menor amplitude, uma espécie de respiração lenta de todo o espectro. Perto de 4,25 micrômetros aparece uma variação com o mesmo formato temporal da feição do monóxido de carbono. Acima de 4,8 micrômetros as amplitudes diminuem. E o mais revelador: as curvas de luz mantêm uma forma semelhante em quase todas as faixas, como se um mesmo relógio marcasse o ritmo, com a amplitude modulada pelo que cada molécula faz naquele comprimento de onda. É aqui que a espectroscopia de resolução média mostra seu poder, porque permite perguntar, molécula por molécula, quem está variando e quem está parado.

Figura 3 — Curvas de luz normalizadas em faixas de 0,1 micrômetro, com a curva de luz integrada em todo o espectro sobre
Figura 3 — Curvas de luz normalizadas em faixas de 0,1 micrômetro, com a curva de luz integrada em todo o espectro sobreposta em cinza. As faixas entre 4,5 e 4,8 micrômetros, dominadas pelo CO, têm as maiores amplitudes; abaixo de 3,6 micrômetros o sinal-ruído é baixo.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327

Para responder a essa pergunta, a equipe calculou a opacidade de cada gás, isto é, sua capacidade de absorver luz por unidade de volume, usando as abundâncias de água, metano, monóxido de carbono, dióxido de carbono, amônia e fosfina medidas por Rowland em 2024, supondo uma pressão de 1 bar e temperatura de 260 kelvins. A conta é simples na essência: a lei dos gases ideais dá a densidade total de partículas, a fração molar de cada espécie diz quantas delas são daquela molécula, e as seções de choque das bases de dados espectroscópicas compiladas por Roxana Lupu e colegas dizem quanto cada uma absorve. Sobrepondo essas curvas de opacidade ao mapa de variação, surge uma correspondência clara. Onde metano e amônia dominam, de 3,7 a 4,25 micrômetros e de 4,05 a 4,5 micrômetros, respectivamente, a variabilidade é menor que a da curva integrada.

Figura 4 — Região de 3,69 a 4,26 micrômetros: espectros mínimo e máximo, razão entre eles, opacidades de metano, amônia
Figura 4 — Região de 3,69 a 4,26 micrômetros: espectros mínimo e máximo, razão entre eles, opacidades de metano, amônia e fosfina a 1 bar e 260 K, e mapa de variação relativa à curva de luz integrada. Onde metano e amônia absorvem menos, as variações são mais fortes.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327

Nas janelas em que metano e amônia absorvem pouco, a história muda. Entre 4,2 e 4,5 micrômetros, onde o dióxido de carbono e a fosfina têm picos de opacidade encaixados nos vales das outras moléculas, as variações são fortes, e a diferença entre o espectro mais escuro e o mais claro é máxima. A sugestão é direta: ali, o que muda com o tempo é a quantidade aparente dessas moléculas. O auge da variabilidade, entre 4,5 e 4,8 micrômetros, coincide com as bandas do monóxido de carbono, o que faz dele o principal responsável pelas oscilações de brilho de WISE 0855. Uma molécula que, numa atmosfera tão fria, nem deveria existir em quantidade apreciável é a que mais pisca.

Figura 5 — Mesma montagem para 4,19 a 4,51 micrômetros, com as opacidades de CO₂ e água adicionadas. Os picos de CO₂ e f
Figura 5 — Mesma montagem para 4,19 a 4,51 micrômetros, com as opacidades de CO₂ e água adicionadas. Os picos de CO₂ e fosfina coincidem com variações mais intensas, o que sugere mudanças de abundância aparente dessas moléculas.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327
Figura 6 — Região de 4,38 a 4,68 micrômetros, onde as bandas do monóxido de carbono concentram as maiores amplitudes de
Figura 6 — Região de 4,38 a 4,68 micrômetros, onde as bandas do monóxido de carbono concentram as maiores amplitudes de variação de toda a série temporal.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327

A água, por sua vez, atua como um amortecedor. Seus picos de opacidade entre 4,5 e 5,18 micrômetros se sobrepõem às bandas do monóxido de carbono, e justamente onde o vapor absorve mais, as variações relativas são menores. Esse é um dos sinais mais importantes do estudo, e voltaremos a ele. Acima de 4,8 micrômetros, nos intervalos entre as linhas de água, aparecem variações que não se alinham com nenhuma molécula esperada, um indício de que parte da variabilidade de WISE 0855 não vem de mudanças de abundância química, mas de outra coisa: temperatura, ou nuvens. Os autores reconhecem que a comparação é aproximada, porque assume uma única temperatura e pressão para todo o espectro, e deixam a resposta para o confronto com os modelos atmosféricos.

Figura 7 — Região de 4,48 a 4,88 micrômetros: onde o vapor d'água tem opacidade maior que a do CO, as variações relativa
Figura 7 — Região de 4,48 a 4,88 micrômetros: onde o vapor d’água tem opacidade maior que a do CO, as variações relativas diminuem.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327
Figura 8 — Região de 4,68 a 5,18 micrômetros: entre as linhas de água aparecem variações que não se alinham com nenhuma
Figura 8 — Região de 4,68 a 5,18 micrômetros: entre as linhas de água aparecem variações que não se alinham com nenhuma molécula esperada, sinal de que parte da variabilidade não vem de mudanças de abundância.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327

Para separar essas contribuições sem depender de nenhum modelo, os autores recorreram à análise de componentes principais, técnica estatística que decompõe um conjunto de medidas em padrões independentes, ordenados por quanto cada um explica da variação total. A analogia mais próxima é a de um coral: ouvindo a gravação completa, é difícil dizer quem canta o quê, mas uma escuta sistemática consegue isolar as vozes que mais contribuem para a mudança de volume e timbre ao longo da música. Aplicada aos 44 espectros normalizados, a técnica identificou “autoespectros”, formas espectrais fixas cujas amplitudes variam no tempo. Cada autoespectro pode ser lido como a assinatura de uma estrutura de superfície não resolvida que gira com o objeto, e a curva de seu coeficiente ao longo das onze horas conta como essa estrutura entra e sai de cena.

O primeiro autoespectro responde por 60,5% de toda a variabilidade e tem um formato parecido com o do espectro médio: é a assinatura de uma mudança ampla, de banda larga, do tipo que uma variação de temperatura ou de cobertura de nuvens produziria. O segundo, com 4,7%, tem sua resposta mais forte exatamente nos comprimentos de onda do monóxido de carbono e da fosfina. O terceiro é ruidoso, com coeficientes espalhados em torno de zero, sinal de que dali em diante só resta ruído de medição. Os outros 34,7% da variabilidade se dispersam em componentes individuais de 1,5% ou menos, indício de um comportamento complexo que nenhum par de padrões esgota. O dado central está na forma das curvas: os coeficientes do primeiro e do segundo autoespectro evoluem de maneira distinta ao longo das onze horas, indicando que nascem de estruturas diferentes na superfície.

Figura 9 — Os três primeiros autoespectros da análise de componentes principais (acima, deslocados verticalmente) e o co
Figura 9 — Os três primeiros autoespectros da análise de componentes principais (acima, deslocados verticalmente) e o coeficiente de cada um ao longo do tempo (abaixo). O primeiro domina a amplitude; o segundo tem forma e período distintos; o terceiro é ruído.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327
Figura 10 — Ampliação do segundo autoespectro entre 4,1 e 4,9 micrômetros, comparado às opacidades de CO e fosfina: sua
Figura 10 — Ampliação do segundo autoespectro entre 4,1 e 4,9 micrômetros, comparado às opacidades de CO e fosfina: sua resposta se concentra exatamente onde essas duas moléculas absorvem.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327

Quanto tempo leva WISE 0855 para dar uma volta em torno de si? Ninguém sabe, e o novo estudo não fecha a questão, embora a cerque. Um periodograma de Lomb-Scargle, ferramenta que procura periodicidades numa série temporal, aplicado às curvas de luz acima de 3,7 micrômetros mostra potência significativa em períodos entre 5,4 e 5,6 horas, mas períodos mais longos também são prováveis, e nenhum sinal significativo sobrevive quando se usa apenas 75% da observação. A limitação é estrutural: com onze horas de dados, os períodos detectáveis com segurança ficam na metade desse intervalo. Um trabalho paralelo de Mikayla Wilson e colegas, divulgado em 2025, ajustou trechos da mesma curva de luz com um modelo de processo gaussiano e obteve períodos superiores a 14 horas.

Para situar o resultado, a equipe convolveu os espectros do Webb com os filtros dos canais 1 e 2 da câmera IRAC do Spitzer e os comparou às duas épocas observadas por Esplin em 2015, em 10 de março e 3 de agosto. A fotometria sintética ficou 37% mais fraca que a medida no canal 1, um desvio já visto em outras anãs frias e atribuído por Beiler e colegas a um vazamento de luz no NIRSpec que afeta objetos muito vermelhos; no canal 2 a diferença foi de 7%, dentro da dispersão esperada. O essencial, porém, se preservou: as amplitudes de pico a pico, de 4% a 6%, são semelhantes às do Spitzer.

Figura 11 — Fotometria de WISE 0855 em três épocas nos canais IRAC 1 (azul) e IRAC 2 (vermelho) do Spitzer: março e agos
Figura 11 — Fotometria de WISE 0855 em três épocas nos canais IRAC 1 (azul) e IRAC 2 (vermelho) do Spitzer: março e agosto de 2015 (dados de Esplin et al. 2016, reamostrados em 15 minutos) e dezembro de 2023, esta última sintetizada a partir dos espectros do Webb.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327

Nas épocas antigas, o canal 2 mostrou um período significativo de 9,33 horas em março de 2015 e apenas períodos mais longos que a própria observação em agosto; o canal 1 não exibiu período algum acima do limiar de confiança. Três épocas, três respostas diferentes. WISE 0855 modula seu brilho em escalas de 5 a 10 horas, e só observações com mais de 20 horas contínuas dirão se há um período estável ou uma mistura de rotação, ventos zonais e ondas planetárias de grande escala, como as que Dániel Apai e colegas identificaram em outras anãs marrons. Para efeito de comparação, Júpiter gira em pouco menos de dez horas, e anãs marrons bem estudadas têm períodos de poucas horas, de modo que os valores encontrados estão dentro do esperado para um gigante gasoso solitário.

Figura 12 — Periodograma das curvas de luz do Webb por faixa de comprimento de onda (esquerda), com as linhas tracejadas
Figura 12 — Periodograma das curvas de luz do Webb por faixa de comprimento de onda (esquerda), com as linhas tracejadas marcando a duração total da observação e sua metade; à direita, o período do pico de cada faixa, entre 5,2 e 5,7 horas.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327
Figura 13 — Periodogramas das três épocas de fotometria. A época 1 do IRAC 2 mostra um pico significativo em 9,33 horas;
Figura 13 — Periodogramas das três épocas de fotometria. A época 1 do IRAC 2 mostra um pico significativo em 9,33 horas; a época do Webb, picos em 5,4 e 5,5 horas; em todos os casos há potência em períodos mais longos que a própria observação.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327

A segunda metade do trabalho confronta os dados com modelos atmosféricos autoconsistentes, programas que resolvem o transporte de radiação, a química e a estrutura de temperatura de uma atmosfera a partir de poucos parâmetros globais. Duas famílias foram usadas. A primeira, PICASO, foi acoplada ao código de nuvens Virga, que descreve a condensação da água pelo parâmetro de eficiência de sedimentação, o fsed, a razão entre a velocidade com que as gotas ou cristais caem e a velocidade do gás que sobe; valores altos produzem nuvens finas e comprimidas, valores baixos produzem nuvens espessas e estendidas na vertical. A segunda, coolTLUSTY, mantida por Lacy e Burrows, descreve a nuvem pelo parâmetro CDAMPU, que fixa a razão entre a escala de altura do gás e a das partículas e controla a rapidez com que a nuvem se afina para cima: quanto maior o valor, menos estendida ela é.

Todas as grades fixaram a gravidade superficial em log g = 4 e a metalicidade solar, propriedades que não mudam em onze horas, e varreram temperaturas efetivas de 250 a 275 kelvins em passos de 5 kelvins no caso do PICASO, e de 200 a 450 kelvins em passos de 25 no coolTLUSTY. Houve versões com e sem nuvens, em equilíbrio químico e fora dele. A química fora de equilíbrio é parametrizada pelo coeficiente de difusão turbulenta, o Kzz, que mede o vigor da mistura vertical: a equipe adotou 10⁷ centímetros quadrados por segundo nos modelos PICASO, valor derivado do monóxido de carbono em trabalho anterior da própria Miles, e 10⁶ nas camadas radiativas do coolTLUSTY, com mistura bem mais vigorosa nas camadas convectivas. Grades derivadas ajustaram, uma a uma, as abundâncias de monóxido de carbono, dióxido de carbono e fosfina, para testar se os modelos acertavam a química.

Tabela 1 — Resumo das grades de modelos atmosféricos PICASO e coolTLUSTY usadas no ajuste, com intervalos de temperatura
Tabela 1 — Resumo das grades de modelos atmosféricos PICASO e coolTLUSTY usadas no ajuste, com intervalos de temperatura efetiva, parâmetros de nuvem (fsed, CDAMPU), de mistura (Kzz) e as abundâncias ajustadas de CO, CO₂ e fosfina.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327

O veredicto do espectro médio foi unânime num ponto: em todos os casos, os modelos com nuvens de água ajustaram melhor os dados do que suas versões sem nuvens, e a opacidade das nuvens foi o ingrediente que mais melhorou o ajuste da distribuição espectral de energia. Entre os modelos com nuvens, os coolTLUSTY venceram os PICASO, diferença que os autores atribuem à extensão vertical e à distribuição de partículas assumidas para a nuvem em cada código. O melhor ajuste global foi um modelo coolTLUSTY com nuvens e química fora de equilíbrio, com temperatura efetiva de 250 kelvins, CDAMPU igual a 6 e log Kzz igual a 6. O melhor modelo PICASO exigiu 265 kelvins, fsed igual a 6 e abundâncias ajustadas à mão: monóxido de carbono em log de fração molar −7,0, fosfina em −8,0 e dióxido de carbono em −9,0.

Esse detalhe esconde uma lição. Mesmo com a química fora de equilíbrio ligada, os modelos subestimam as quantidades de dióxido de carbono e fosfina necessárias para reproduzir o espectro; quando se ajusta uma das duas moléculas sozinha, a abundância inferida dispara para compensar a ausência da outra. E nenhum ajuste é perfeito: os valores de qui-quadrado reduzido ficam na casa das dezenas de milhares, o que, num conjunto de dados com barras de erro tão pequenas, revela resíduos que ainda carregam estrutura não capturada pela física dos modelos. Os autores são explícitos ao chamar o exercício de demonstração de quanto os modelos atuais conseguem explicar, e não de descrição definitiva, e pedem trabalho futuro de afinação específica para WISE 0855.

Figura 14 — Melhores ajustes das grades PICASO (vermelho) e coolTLUSTY (azul) ao espectro médio de WISE 0855 (preto), co
Figura 14 — Melhores ajustes das grades PICASO (vermelho) e coolTLUSTY (azul) ao espectro médio de WISE 0855 (preto), com os resíduos em cinza abaixo.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327
Tabela 2 — Melhor modelo de cada grade, com seus parâmetros e os valores de qui-quadrado e qui-quadrado reduzido. Modelo
Tabela 2 — Melhor modelo de cada grade, com seus parâmetros e os valores de qui-quadrado e qui-quadrado reduzido. Modelos com nuvens de água vencem os equivalentes sem nuvens em todos os casos.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327

O passo decisivo foi ajustar não o espectro médio, mas cada um dos 44 espectros individualmente, para ver que parâmetros precisam mudar de quinze em quinze minutos. Com o PICASO, a equipe costurou em 4,45 micrômetros a grade de monóxido de carbono ajustável com a grade de fosfina e dióxido de carbono ajustáveis, o que permitiu variar as três abundâncias ao mesmo tempo, e testou quatro cenários: mudar temperatura e as três abundâncias com nuvem fixa; mudar só a temperatura; mudar só a nuvem; mudar nuvem e temperatura. Com o coolTLUSTY, três cenários: só a temperatura, só a nuvem, ou as duas. Em cada espectro, um algoritmo de minimização encontrou a combinação de parâmetros que melhor reproduzia os dados, com o ruído de medição simulado por sorteios aleatórios a partir das incertezas.

Figura 15 — Mapas de variação percentual dos dados (a) e dos melhores ajustes PICASO quando se deixa variar temperatura
Figura 15 — Mapas de variação percentual dos dados (a) e dos melhores ajustes PICASO quando se deixa variar temperatura e abundâncias com nuvem fixa (b), só a temperatura (c), só a nuvem (d) e nuvem e temperatura (e).Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327
Figura 16 — Mapas equivalentes para o coolTLUSTY: dados (a), só temperatura com CDAMPU fixo (b), só CDAMPU com temperatu
Figura 16 — Mapas equivalentes para o coolTLUSTY: dados (a), só temperatura com CDAMPU fixo (b), só CDAMPU com temperatura fixa (c) e ambos variando (d).Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327

Os resultados desenham o retrato de um mundo em que a temperatura efetiva, entendida aqui como o calor aparente de um hemisfério heterogêneo em rotação, não como uma propriedade que muda de verdade, oscila 0,6 kelvin nos modelos PICASO e 1,3 kelvin nos coolTLUSTY, acompanhando a curva de luz integrada e a amplitude do primeiro autoespectro. A abundância de monóxido de carbono varia 0,05 dex, entre log −6,88 e −6,83, e a de fosfina entre −8,38 e −8,18; a de dióxido de carbono ficou presa em −9,0, o piso da grade, sinal de que a grade não cobre toda a dinâmica dos dados. Um detalhe merece segunda leitura: quando a temperatura cai, a química deveria favorecer o metano em detrimento do monóxido de carbono, mas o comportamento observado é o oposto. A curva de CO tem a forma do segundo autoespectro, não a da temperatura. Ela nasce em outro lugar.

Nos modelos coolTLUSTY, reproduzir a maior parte das mudanças relativas exigiu variar temperatura e espessura da nuvem juntas: o parâmetro de nebulosidade mudou em 0,6, de 4,0 para 4,6, o que corresponde a uma camada de nuvem que encolhe pela metade em altitude e perde opacidade. Os dois códigos discordam na temperatura média, 264 kelvins no PICASO contra 256 no coolTLUSTY, mas concordam na forma da tendência ao longo das onze horas. A diferença de altitude das nuvens tem causa identificável: no coolTLUSTY, a base da nuvem fica em pressões mais baixas, mais alto na atmosfera, porque o ponto em que a curva de condensação cruza o perfil de temperatura é definido pela versão sem nuvens do modelo. Isso reforça uma conclusão de James Mang de 2022: a forma como se descreve a nuvem decide o que o modelo consegue reproduzir.

Figura 17 — Evolução temporal dos parâmetros de melhor ajuste: temperatura efetiva e abundâncias de CO, CO₂ e fosfina no
Figura 17 — Evolução temporal dos parâmetros de melhor ajuste: temperatura efetiva e abundâncias de CO, CO₂ e fosfina no PICASO (esquerda); temperatura efetiva e parâmetro de nebulosidade CDAMPU no coolTLUSTY (direita). As abundâncias não acompanham a temperatura.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327

Para entender por que monóxido de carbono e fosfina existem numa atmosfera tão fria, é preciso um conceito que a química de Júpiter ensinou aos astrônomos: o quenching, ou congelamento químico. Nas profundezas quentes da anã, as reações entre as moléculas são rápidas e o monóxido de carbono é abundante. À medida que a convecção empurra o gás para cima e ele esfria, a química deveria convertê-lo em metano e água, mas as reações desaceleram exponencialmente com a queda da temperatura. Chega um nível em que a mistura ergue o gás mais depressa do que as reações conseguem transformá-lo, e a composição fica congelada no valor daquela altitude. O ponto em que o tempo químico iguala o tempo de mistura é o nível de quench, e a abundância observada no topo é uma fotografia da química que acontece lá embaixo.

As três reações químicas líquidas assumidas para o congelamento (quench) do monóxido de carbono, do dióxido de carbono e
As três reações químicas líquidas assumidas para o congelamento (quench) do monóxido de carbono, do dióxido de carbono e da fosfina.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327

Invertendo o raciocínio, a equipe usou as abundâncias ajustadas para calcular onde, no perfil do melhor modelo PICASO em equilíbrio, cada molécula teria sido congelada, e que valor de Kzz isso implica. As reações de referência são a conversão de CO em metano e água, a de CO2 em CO e água, e a oxidação da fosfina em óxido de fósforo. A fosfina congela entre 23,5 e 24,8 bar, o dióxido de carbono abaixo de 51 bar, e o monóxido de carbono entre 68,7 e 70,7 bar, todos dentro de uma mesma escala de altura de pressão. Os coeficientes de mistura inferidos, contudo, não batem: cerca de um milhão de centímetros quadrados por segundo para o CO, pelo menos 3 × 10⁸ para o CO2, quase duas ordens de grandeza acima, e ínfimos 23 a 44 para a fosfina, o que implicaria quase nenhuma mistura.

Figura 18 — Perfil de pressão e temperatura do modelo PICASO de 260 K com as curvas de condensação de água, NH₄H₂PO₄, Zn
Figura 18 — Perfil de pressão e temperatura do modelo PICASO de 260 K com as curvas de condensação de água, NH₄H₂PO₄, ZnS, KCl e Na₂S e as pressões de congelamento de CO e fosfina (esquerda); escalas de tempo químicas contra a pressão, com linhas de Kzz constante (centro); e perfis de abundância em equilíbrio e fora dele (direita).Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327
Tabela 3 — Pressões de congelamento e coeficientes de difusão turbulenta inferidos para fosfina, CO₂ e CO: as três moléc
Tabela 3 — Pressões de congelamento e coeficientes de difusão turbulenta inferidos para fosfina, CO₂ e CO: as três moléculas congelam na mesma escala de altura, mas os valores de Kzz não são consistentes entre si.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327

Num modelo unidimensional com mistura constante, moléculas que congelam a profundidades parecidas deveriam partilhar o mesmo Kzz. Que isso não aconteça diz que algo falta no quadro. No caso da fosfina, os autores lembram que a reação limitante e o sumidouro principal dos compostos de fósforo ainda são mal conhecidos. Para o dióxido de carbono, cuja abundância depende diretamente da do monóxido de carbono, a tensão é mais séria, e a suspeita recai sobre a metalicidade: a razão entre as duas moléculas é muito sensível à abundância de elementos pesados, e uma composição diferente da solar pode ser o que falta para reconciliar os números. A conclusão prudente é que, com um modelo de abundância solar, não se obtém um mesmo Kzz para todas as moléculas, e a intensidade da mistura não aumenta com a profundidade do congelamento.

Chama a atenção, ainda, que a fosfina esteja bem abaixo do esperado pela teoria, repetindo uma tendência de esgotamento já vista em outras anãs marrons por Beiler em 2024 e por Adam Burgasser em 2025. Em Júpiter e Saturno, a fosfina é o traçador clássico da mistura vertical, porque só sobrevive nas camadas visíveis se for reposta continuamente de baixo; em WISE 0855 ela aparece, mas em quantidade modesta, e varia em sincronia com o monóxido de carbono. Essa sincronia, confirmada pela análise de componentes principais, é uma das pistas mais fortes do artigo, porque restringe o que pode estar modulando as duas moléculas ao mesmo tempo. Chega-se então à pergunta do título: de que são feitas as nuvens?

Modelos de Morley de 2012 preveem que, mais fundo na atmosfera de WISE 0855, existam nuvens de sulfeto de zinco, cloreto de potássio e sulfeto de sódio, e os dados permitem interrogá-las. As pressões de congelamento de CO, CO2 e fosfina cruzam as curvas de condensação dessas três espécies; se houvesse uma camada de nuvem entre os níveis de congelamento do CO e do CO2, as duas moléculas não poderiam partilhar o mesmo Kzz. A fosfina congela um pouco acima dessas curvas; o CO, abaixo de todas. A próxima nuvem possível, de sulfeto de manganês, só condensa abaixo de cerca de 1.550 kelvins; como uma nuvem se estende por no máximo meia ordem de grandeza em pressão, ela só poderia modular o CO, não a fosfina. Os buracos nas nuvens superiores, concluem os autores, ou atravessam as camadas de sulfetos e cloreto de potássio, ou essas camadas não existem.

Como a fosfina e o monóxido de carbono variam juntos, as únicas nuvens capazes de produzir manchas que modulem as duas moléculas ao mesmo tempo são a de água e a de fosfato diácido de amônio, composto em que a fosfina pode condensar. A escala de tempo química da fosfina é incerta, mas ela deve ser a principal molécula portadora de fósforo. A partir desse quadro, os autores sugerem que as variações das moléculas fora de equilíbrio nascem de mudanças de temperatura em alta pressão, e não de nuvens profundas esburacadas, e observam que a nuvem a 10 bar inferida por Kühnle, além de não coincidir com nenhuma das espécies teoricamente esperadas, foi obtida num arcabouço que desfavorece mudanças no perfil de temperatura.

A água entra em cena pela porta dos fundos, pelas linhas que se recusam a variar. Nos modelos PICASO, a abundância de vapor é praticamente constante acima de 5 bar, de modo que as linhas de água não deveriam responder nem à convecção nem a mudanças de energia ao longo do perfil, e a grade falha em reproduzir as pequenas variações observadas dentro dessas linhas. Os modelos coolTLUSTY, em que a nuvem de água reside mais baixo na atmosfera, reproduzem melhor exatamente essa região: a água só deixa de ser constante onde começa a condensar.

A inferência é que as nuvens de WISE 0855 ficam em altitude menor do que o PICASO assume, e que a atmosfera é heterogênea por causa de nuvens de água de espessura variável. Daí a conclusão central do artigo: as nuvens superiores são provavelmente de água pura, e as regiões onde elas são mais finas expõem níveis mais profundos, com abundâncias realçadas de monóxido de carbono e fosfina.

Figura 19 — Região dominada pelo vapor d'água: opacidades de água e CO (acima) e razão entre os espectros mínimo e máxim
Figura 19 — Região dominada pelo vapor d’água: opacidades de água e CO (acima) e razão entre os espectros mínimo e máximo comparada aos modelos coolTLUSTY que variam só a temperatura (azul) ou temperatura e nebulosidade (vermelho). As nuvens de água amortecem as variações dentro das bandas de água.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327

Uma explicação alternativa, a de gradientes de mistura ao longo da longitude, com regiões em que a convecção é mais vigorosa, é descartada pelos autores com um argumento físico: a mistura convectiva deve escalar com a velocidade de rotação, que só muda com a latitude. Observações da sonda Juno, que resolveram espacialmente a atmosfera de Júpiter, apoiam o mesmo raciocínio. Restam duas fontes principais de variabilidade: mudanças de temperatura ou de nuvens em alta altitude, responsáveis pelo padrão de banda larga que domina 60,5% do sinal; e variações de temperatura em pressões profundas, no nível de congelamento do monóxido de carbono ou abaixo dele, que modulam CO e fosfina em conjunto. Oscilações de cerca de 1 kelvin, como as medidas, são compatíveis com as perturbações de temperatura previstas por Adam Showman e Yohai Kaspi em 2013 para uma anã marrom com período de dez horas.

Como está orientada essa anã em relação a nós? A amplitude elevada das modulações das moléculas fora de equilíbrio torna improvável que estejamos olhando para um dos polos, e a equipe conclui que WISE 0855 provavelmente é vista de lado, com o equador voltado para a Terra. Simulações de Mark Hammond e colegas, de 2023, indicam que anãs Y devem ser dominadas por jatos zonais, correntes de vento que circundam o objeto em faixas de latitude, como as bandas de Júpiter, e que vórtices estáveis de grande escala só sobrevivem perto dos polos, onde a Juno, aliás, encontrou os ciclones agrupados de Júpiter. Para ventos de 10, 100 e 1.000 metros por segundo, as escalas características de jatos e vórtices mostram que um vórtice só supera um jato acima de 78 graus de latitude, e apenas com os ventos mais rápidos.

Figura 20 — Raio de deformação de Rossby (vórtices) e escala de Rhines (jatos) em função da latitude, para um período de
Figura 20 — Raio de deformação de Rossby (vórtices) e escala de Rhines (jatos) em função da latitude, para um período de 11 horas e ventos de 10, 100 e 1.000 m/s: os jatos zonais dominam em quase todas as latitudes.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327

As contribuições do monóxido de carbono e da fosfina ocorrem em escalas de tempo curtas e têm forma aproximadamente senoidal, o que, dadas as escalas atmosféricas relevantes, leva os autores a atribuí-las a variações de temperatura ao longo de jatos zonais em latitudes equatoriais e médias. Regiões polares poderiam contribuir com variabilidade de longo prazo, mas não há, nos dados do Spitzer e do Webb, evidência de mudança no período. Decidir se as variações profundas de temperatura resultam de jatos deslocados por ventos de leste e oeste ou de uma tempestade estagnada em baixa latitude exigirá observações mais longas, em várias épocas, com o Webb. É a mesma receita que resolveria o período de rotação: mais de vinte horas seguidas sobre o alvo.

O esquema que resume o artigo é quase um corte geológico de uma atmosfera. No topo, em torno de 1 bar, uma camada de nuvens de água, heterogênea e, nos modelos usados, estática. Abaixo dela, uma região de gás rica em metano, amônia e vapor d’água. Entre cerca de 1 e 4 bar fica a faixa de pressões que o NIRSpec enxerga nesse modo de observação. E de muito mais fundo, de perto de 70 bar, a convecção traz gás enriquecido em monóxido de carbono, dióxido de carbono e fosfina até a fotosfera. Na face projetada, regiões cobertas por nuvens modulam o brilho em escalas de cerca de 1 kelvin, enquanto regiões de nuvens mais finas deixam a luz vir de níveis mais profundos, com feições mais intensas das moléculas fora de equilíbrio. Dois relógios, duas altitudes, uma mesma esfera girando.

Figura 21 — Esquema da atmosfera de WISE 0855: estrutura vertical com nuvens de água perto de 1 bar e gás enriquecido em
Figura 21 — Esquema da atmosfera de WISE 0855: estrutura vertical com nuvens de água perto de 1 bar e gás enriquecido em CO, CO₂ e fosfina trazido de cerca de 70 bar (a); face visível com regiões nebulosas (roxo) e de nuvens finas (azul) (b); e as duas curvas temporais associadas, de temperatura/nuvens e de moléculas profundas.Crédito: Miles et al. (2026), arXiv:2609.20664 · observações NASA/ESA/CSA JWST, programa GO 2327

A importância do resultado ultrapassa WISE 0855. Anãs Y são os melhores substitutos observáveis para os planetas gigantes frios que orbitam outras estrelas, cuja luz é quase sempre abafada pela da estrela hospedeira. Tudo o que se aprende aqui, como a mistura vertical opera, onde a água condensa, como nuvens parciais deixam escapar a assinatura das profundezas, vale como ensaio para o dia em que os telescópios conseguirem obter séries temporais de exoplanetas gelados. O trabalho também mostra que a estratégia de encarar um objeto por horas, em vez de minutos, abre uma janela que a fotometria do Spitzer apenas entreabria: é possível distinguir se o que muda é a temperatura, a nuvem ou a química, e a que profundidade cada mudança acontece.

Há arestas por lixar, e os autores as listam. A metalicidade não solar precisa ser explorada para resolver a tensão entre os valores de Kzz. A grade de dióxido de carbono precisa se estender a abundâncias menores. Os modelos de nuvem dos dois códigos colocam a água em altitudes diferentes, e as evidências favorecem a versão mais baixa. O raio inferido do objeto continua pequeno demais em comparação com os modelos evolutivos. E o período de rotação, peça-chave para qualquer mapa de superfície, permanece em aberto. Nada disso enfraquece a conclusão de que a anã marrom mais fria conhecida tem modulações atmosféricas complexas, guiadas por moléculas fora de equilíbrio e por nuvens de água; apenas desenha o programa de trabalho da próxima década.

Nos agradecimentos do artigo, Brittany Miles dedica o trabalho ao avô, Alvin Turner Jr. É um gesto pequeno num texto cheio de autoespectros e coeficientes de difusão, mas serve de lembrete de para que serve tudo isso. A 7,4 anos-luz de distância, um mundo que nunca foi estrela gira a cada punhado de horas sob um céu sem sol, com nuvens de água rasgadas aqui e ali por clareiras pelas quais sobe o hálito quente de suas entranhas. Ninguém jamais o verá de perto. Mas, por onze horas numa noite de dezembro, um telescópio a um milhão e meio de quilômetros da Terra acompanhou suas nuvens se abrindo e fechando, e trouxe de volta, em 44 espectros, a primeira previsão do tempo de um mundo que não é nosso.

Fonte
Miles, B. E.; Mang, J.; Morley, C. V.; Skemer, A. J. I.; Rowland, M.; Lacy, B.; Suarez, G.; Visscher, C.; Apai, D.; Bjoraker, G.; Carter, A. L.; Faherty, J. K.; Fortney, J. J.; Fuda, N.; Geballe, T.; Kothari, H.; Limbach, M. A.; Marley, M.; McCarthy, A. M.; Schneider, A.; Vos, J.; Cushing, M.; Freedman, R.; Line, M.; Lupu, R.; Martin, E.; Wilson, M. J. (2026). “Water Cloud and Chemical Modulations in the Coldest Brown Dwarf”. Pré-publicação arXiv:2609.20664v1 [astro-ph.EP], 17 de setembro de 2026. DOI: 10.48550/arXiv.2609.20664 · arxiv.org/abs/2609.20664