O estudo, assinado por Edward D. Young, Sarah P. Marcum, Aaron Werlen e Paula N. Wulff, do Departamento de Ciências da Terra, Planetárias e Espaciais da UCLA, saiu na edição de 1º de outubro de 2026 do The Astrophysical Journal, em acesso aberto, com o título “Ice Giants Revisited: Uranus and Neptune as Magma Ocean Worlds”. A tese é direta: raios, densidades médias, harmônicos gravitacionais, momentos de inércia, luminosidades internas e as principais feições da química atmosférica de Urano e Netuno são compatíveis com interiores formados por um oceano de magma supercrítico, rico em hidrogênio, coberto por um envelope de hidrogênio e hélio. E o modelo faz isso com apenas três parâmetros ajustáveis por planeta, uma economia rara num campo acostumado a estruturas de muitas camadas e muitas incógnitas.

O rótulo de gigantes de gelo nasceu de um raciocínio cosmoquímico, não de uma observação. Nos anos 1960, Reynolds e Summers e, na década seguinte, Morris Podolak com Alastair Cameron calcularam o que deveria condensar num disco protoplanetário de composição solar além da snow line, a linha de gelo, o raio a partir do qual a água congela. A conta indicava que água, metano e amônia, agrupados sob o nome genérico de “gelos”, deveriam superar a massa do material rochoso por um fator de aproximadamente dois. Daí o modelo canônico de três camadas, formalizado por William Hubbard e Joseph MacFarlane em 1980: um núcleo de rocha, um manto de gelos e, por cima, um envelope de hidrogênio e hélio. Urano e Netuno seriam, nessa visão, bolas de gelo com um pouco de pedra no meio e uma atmosfera fina por fora.

O problema é que os dados nunca exigiram essa arquitetura. O que se mede de fora de um planeta gigante é a massa, o raio, a forma do campo gravitacional e o calor que escapa. Modelos empíricos que ajustam o campo gravitacional sem pressupor composição alguma mostram que interiores com gradientes contínuos de densidade reproduzem as observações tão bem quanto camadas discretas, e que interiores dominados por rocha não podem ser descartados, como Ravit Helled e colaboradores vêm argumentando desde 2011 e Nicholas Teanby e colegas reforçaram em 2020. Com muitos parâmetros livres e poucos observáveis, surgem degenerescências: famílias inteiras de interiores diferentes produzem os mesmos números do lado de fora. Nos últimos dois anos, novas variantes com frações maiores de rocha, e até fases superiônicas exóticas, entraram na disputa.

O primeiro golpe dos autores no consenso vem de uma planilha de abundâncias solares. Para que gelo supere rocha na proporção de dois para um, não basta a água congelar: é preciso que o monóxido de carbono, o CO, também condense e seja incorporado aos planetesimais, o que só acontece além da linha de gelo do CO, muito mais fria e distante que a da água. Atribuindo todo o carbono de um gás de composição solar ao CO, o oxigênio restante primeiro aos silicatos e depois à água, e o ferro ao metal, a equipe chega a um inventário de aproximadamente 54% de CO, 23% de silicatos, 12% de água e 11% de ferro, em massa. Se todo esse CO virasse gelo, o material condensado além da linha do CO seria de fato 66% gelo. Entre as duas linhas de gelo, porém, o CO permanece gasoso e não contribui para os sólidos. A matéria-prima dos planetas, nessa região, fica em torno de 50% de silicatos, 27% de água e 23% de ferro.

A água sozinha, portanto, responde por cerca de um quarto do material planetário disponível fora da linha de gelo da água e dentro da linha do CO, bem longe da razão de dois para um que sustenta boa parte dos modelos de gigantes de gelo. A mesma conta tem um bônus elegante: ela explica a proporção de aproximadamente dois para um entre silicatos e ferro metálico que se observa na Terra. O momento da acreção, a condensação incompleta e a eficiência com que o gelo de CO seria capturado alteram os números exatos, mas o recado central resiste. Uma composição dominada por gelo exige que os progenitores de Urano e Netuno tenham acumulado grandes quantidades de material rico em CO vindo de além da linha de gelo desse composto.

E onde ficava essa linha? Ninguém sabe ao certo para o disco que formou o Sistema Solar, mas os discos ao redor de estrelas jovens parecidas com o Sol oferecem pistas. Em TW Hydrae, frequentemente tratada como análogo da nebulosa solar, a emissão do íon N₂H⁺ coloca a linha de gelo do CO perto de 30 unidades astronômicas, aproximadamente a posição atual de Netuno e além da órbita de Urano, enquanto observações de isotopólogos raros do CO a empurram para cerca de 20 unidades astronômicas, semelhante ao raio orbital de Urano. Em AS 209, IM Lup e GM Aur, modelos baseados em isotopólogos apontam valores de 12 a 30 unidades astronômicas, com incertezas de 5. Análises com N₂H⁺ em LkCa 15, GM Aur e DM Tau indicam raios de 58, 48 e 75 unidades astronômicas para a linha do CO e de 88, 78 e 145 para a do nitrogênio molecular. Durante a fase T Tauri do Sol, portanto, essas fronteiras podem ter ficado além das órbitas atuais dos dois planetas.

Os sobreviventes da formação planetária no Sistema Solar exterior contam a mesma história. Os objetos do Cinturão de Kuiper, talvez o melhor testemunho direto da composição daquela região, são cerca de 70% silicatos, segundo a análise de Carver Bierson e Francis Nimmo de 2019. Os gelos de Plutão são dominados por água; o CO, o metano e o nitrogênio congelados que a New Horizons fotografou em 2015 formam apenas um verniz fino na superfície. Os cometas da família de Júpiter, que amostram planetesimais ricos em voláteis nascidos no mesmo cenário, exibem razões entre poeira e gelo em torno de três, o que corresponde a frações de gelo de 25%, conforme as medidas da sonda Rosetta no cometa 67P/Churyumov-Gerasimenko analisadas por Martin Pätzold e por Marco Fulle. Nenhuma dessas testemunhas sustenta progenitores feitos majoritariamente de gelo.

O segundo argumento desloca a discussão do lugar de nascimento para a química do interior. Trabalhos recentes sobre os sub-Netunos, planetas com dois a quatro raios terrestres que não existem no Sistema Solar mas são os mais comuns da galáxia, mostraram que quantidades copiosas de água podem ser fabricadas dentro do planeta, pela reação entre uma atmosfera primordial de hidrogênio e o núcleo incandescente sob ela, como demonstraram Edwin Kite e Laura Schaefer em 2021, Hilke Schlichting e o próprio Young em 2022 e uma série de estudos posteriores. A quantidade de água de um planeta, portanto, não é um registro fiel de onde ele acretou. Corpos formados dentro da linha de gelo da água podem terminar com concentrações de água semelhantes às de corpos formados fora dela.

Mais do que isso, a memória química do gelo acretado tende a ser apagada. Os autores propõem um experimento mental: um planeta montado a partir de planetesimais com 25% de gelo de água, 50% de rocha na forma de MgSiO₃ e 25% de ferro, uma mistura parecida com a da Terra. Simulações de dinâmica molecular baseadas na teoria do funcional da densidade, publicadas por Tanja Kovačević e colegas em 2022, indicam que MgSiO₃ fundido e água se tornam inteiramente miscíveis a pressões e temperaturas elevadas. Os 25% de água contribuem com 22,2% de oxigênio e 2,8% de hidrogênio para a massa do planeta. O oxigênio se incorpora ao componente denso, semelhante a um óxido, da fase supercrítica; o hidrogênio engrossa o componente leve. O efeito de primeira ordem de acretar água, concluem, é aumentar o estoque efetivo de hidrogênio, não construir um manto de gelo.

Há ainda uma terceira frente. Quem defende que Urano e Netuno são mais rochosos do que se pensava precisa encarar o que “rocha” significa quando misturada a um reservatório imenso de hidrogênio a milhares de graus e a pressões de dezenas a centenas de gigapascais. Simulações ab initio, feitas a partir das equações fundamentais da mecânica quântica, do grupo de Young em 2024 e 2025 e de Gilmore e Lars Stixrude em 2026 mostram que silicato, hidrogênio e ferro são miscíveis nessas condições, e experimentos de laboratório de Francesca Miozzi e colaboradores, publicados na Nature em 2025, apoiam o resultado. Rocha, nesse regime, não é um sólido nem sequer um líquido de silicato convencional. É um fluido supercrítico em que ferro, silicato e hidrogênio se dissolvem uns nos outros sem fronteiras internas.

É desse cenário que emerge o conceito de gas dwarf, o anão gasoso, termo que a equipe usa para planetas formados por um oceano de magma coberto por um envelope de hidrogênio, mas pequenos demais para disparar o crescimento descontrolado que produz um Júpiter. A arquitetura resultante tem quatro camadas. No fundo, o oceano de magma supercrítico de silicato, ferro e hidrogênio. Logo acima, uma camada-limite onde a convecção é inibida, reforçada por gradientes de peso molecular. Depois, um envelope convectivo rico em hidrogênio. No topo, uma camada radiativa, onde a energia escapa por difusão de radiação em vez de ser carregada por correntes de gás. A fronteira entre o oceano e o envelope não é uma interface gradual: é uma transição de fase de primeira ordem, definida por uma curva chamada binodal, que depende da pressão e separa, no espaço de temperatura e composição, a mistura homogênea das duas fases que coexistem.

Uma imagem doméstica ajuda. Azeite e vinagre, batidos com vigor e aquecidos, formam um molho uniforme; deixados esfriar e descansar, separam-se em duas camadas com composições distintas. No interior dos dois planetas, a temperatura e a pressão são altas o bastante para que silicato e hidrogênio formem uma única solução. À medida que se sobe em direção à superfície e a temperatura cai, o fluido cruza a binodal e se divide em duas fases: uma rica em hidrogênio, que se torna gás a pressões mais baixas, e um condensado líquido de silicato. Os autores tratam esse condensado como uma chuva de rocha que precipita de volta à base do envelope, num processo que modelam como uma destilação do tipo Rayleigh. Essa chuva carrega magnésio, silício e oxigênio para baixo e deixa o gás acima progressivamente mais puro em hidrogênio, criando o gradiente de peso molecular que estabiliza a base do envelope contra a convecção.

O oceano de magma desse modelo pouco lembra a lava de um vulcão. Com hidrogênio dissolvido, a densidade cai muito abaixo da do MgSiO₃ puro, e o fluido fica mais compressível a alta pressão: a derivada do módulo de compressibilidade em relação à pressão, uma medida de quão depressa o material endurece quando comprimido, cai de 5,71 no silicato seco para 4,32 numa mistura com 4% de hidrogênio em massa, segundo as simulações de Marcum e colaboradores. A temperatura de fusão também despenca, cerca de 1.000 °C abaixo da do silicato puro, pelo efeito clássico da depressão do ponto de fusão por mistura.

O ferro, que num planeta terrestre afundaria para formar um núcleo metálico, permanece dissolvido no fundido de silicato e hidrogênio. O resultado é um interior inteiramente fundido, composto por uma única fase convectiva, sem núcleo e sem manto. Os próprios autores admitem, numa nota de rodapé, que a mistura é tão rica em hidrogênio que a palavra “magma” talvez já não seja de todo adequada; eles a mantêm para distinguir essa fase do envelope, ainda mais rico em hidrogênio.

Para testar a ideia, é preciso lembrar quão pouco se sabe de concreto sobre os dois planetas. Os raios no nível de pressão de 1 bar, 24.766 quilômetros para Netuno e 25.559 para Urano, e as primeiras determinações dos harmônicos gravitacionais J₂ e J₄ vieram das passagens da Voyager 2. Esses harmônicos descrevem como a distribuição de massa deforma o campo gravitacional de um corpo em rotação; o J₂ mede o achatamento, e o J₄ é mais sensível às camadas externas. Desde então, os valores foram refinados pela astrometria de longo prazo das órbitas dos satélites, com trabalhos de Robert Jacobson, Marina Brozović e, em 2023, de B. Wang e colegas para Netuno, e de Richard French e coautores, em 2024, para Urano. O momento de inércia normalizado, que registra quão concentrada a massa está no centro, assim como um patinador gira mais depressa ao encolher os braços, nunca foi medido diretamente em nenhum dos dois; existem apenas valores derivados de modelos, 0,2410 para Netuno e 0,2250 para Urano.

Há armadilhas até nos números básicos. O período de rotação de Netuno adotado nos modelos, 58.000 segundos, vem da periodicidade das emissões de rádio moduladas pelo campo magnético, medida pela Voyager 2, já que as nuvens visíveis giram em ritmos diferentes conforme a latitude por causa dos ventos, e sabe-se que ele é inconsistente com a forma achatada observada do planeta, uma fonte de erro sistemático que nenhuma incerteza astrométrica captura.

As luminosidades intrínsecas, o calor que cada planeta emite além do que recebe do Sol, dependem do albedo de Bond, a fração da luz solar refletida, que é a principal fonte de incerteza. Os dados do instrumento IRIS da Voyager indicavam razões entre potência emitida e absorvida de 1,06 ± 0,08 para Urano e 2,61 ± 0,28 para Netuno. Netuno irradia cerca de 3,3 × 10¹⁵ watts de calor interno; Urano, que durante décadas pareceu quase morto por dentro, teve sua luminosidade revisada para cima, para aproximadamente 6,3 × 10¹⁴ watts, em reanálises do albedo publicadas em 2025 por Patrick Irwin e por W. Wang e colegas. Mesmo assim, Netuno continua emitindo cinco vezes mais que o vizinho, um enigma que atravessa gerações de planetólogos.

A química das atmosferas acrescenta restrições de outra natureza. Ambas são dominadas por hidrogênio e hélio, com o metano em terceiro lugar. Em Netuno, a fração molar de metano na troposfera profunda, abaixo do nível de condensação de cerca de 1,5 bar, varia de 6% a 7% nas latitudes equatoriais a cerca de 3% perto dos polos; Urano exibe aproximadamente metade disso, de 1% a 4% nas baixas latitudes e cerca de 1% nos polos. Esses teores implicam enriquecimentos de carbono de cerca de 80 vezes o valor solar em Netuno e de 40 a 50 vezes em Urano.

O sulfeto de hidrogênio, H₂S, foi detectado entre 1 e 4 bar nos dois planetas, com frações de 0,4 a 0,8 parte por milhão em Urano e de 1 a 3 partes por milhão em Netuno. A amônia nunca foi detectada em nível algum de nenhum dos dois, o que implica uma razão entre enxofre e nitrogênio bem acima da solar nas regiões acessíveis. Qualquer modelo de interior que se preze precisa explicar por que o enxofre aparece e o nitrogênio some.

A ferramenta da equipe chama-se Planet_LAB, um código aberto, disponível no GitHub e depositado no Zenodo, que constrói um planeta esférico unidimensional a partir de três parâmetros: a pressão na binodal, que fixa onde o oceano de magma termina e o envelope começa; a fração de massa total de hidrogênio do planeta; e a razão entre o número de Rayleigh da atmosfera e seu valor crítico, uma medida do vigor da convecção que determina a espessura da camada-limite na base do envelope.

O código integra o equilíbrio hidrostático em 120 mil passos radiais, do centro até o nível de 1 bar, atravessando a fronteira definida pela binodal. A forma do planeta em rotação e os harmônicos J₂ e J₄ são calculados a posteriori pelo método dos esferoides de Maclaurin concêntricos, de William Hubbard, com um parâmetro de rotação de aproximadamente 0,026, igual para os dois planetas. Buscas por cadeias de Markov e Monte Carlo varrem o espaço tridimensional de parâmetros atrás da região de máxima verossimilhança para seis observáveis: J₂, J₄, o momento de inércia, o raio equatorial a 1 bar, a temperatura a 1 bar e a luminosidade intrínseca.

A filosofia é diferente da dos modelos clássicos. Na abordagem tradicional, impõem-se massa, raio, temperatura e rotação como condições de contorno e ajustam-se parâmetros internos, como as metalicidades de duas camadas de envelope e a massa do núcleo, até reproduzir o campo gravitacional, um problema inverso que fica subdeterminado por dois graus de liberdade quando se tenta encaixar quatro incógnitas em dois harmônicos. No modelo da UCLA, nada é imposto: para cada trio de parâmetros, o código prevê o conjunto completo de observáveis, e a concordância com as medidas é consequência de encontrar regiões de alta verossimilhança, não uma condição de partida. É a lógica bayesiana de avaliar a probabilidade de uma hipótese física diante dos dados, com precedentes nos trabalhos de Elizabeth Bailey e David Stevenson, de 2021, e de Ludwig Scheibe e colegas. E com apenas três botões para girar, o modelo é difícil de forçar a dizer o que se quer ouvir.

Um detalhe técnico que costuma ser varrido para debaixo do tapete recebeu atenção especial: os ventos. Os jatos zonais de Urano e Netuno, estes últimos os mais velozes do Sistema Solar, com rajadas de cerca de 2.000 quilômetros por hora, redistribuem massa e deixam uma marca dinâmica nos harmônicos gravitacionais, que precisa ser subtraída antes de comparar o campo medido com o de um planeta estático.

A maioria dos trabalhos anteriores assumia que os ventos só alteram o J₄; a equipe, seguindo os métodos de Wieland Dietrich e Johannes Wicht adaptados por Deniz Soyuer, permitiu que eles se estendessem até a base da atmosfera convectiva, com um perfil de decaimento em profundidade. Para Urano, com ventos que descem até 4.584 quilômetros, a correção em J₂ chega a −13,47 × 10⁻⁶, trinta vezes a incerteza de medida; em Netuno, com ventos até 3.838 quilômetros, a correção é de apenas 0,87 × 10⁻⁶. As correções em J₄ valem −3,46 × 10⁻⁶ e −3,88 × 10⁻⁶. Em Urano, portanto, os ventos reescrevem o número que se pretende ajustar.

Os resultados para Netuno saem de uma busca com 1.100 modelos. A solução de máxima verossimilhança tem a binodal a 5,42 gigapascais, o equivalente a cerca de 54 mil atmosferas, perto de cinquenta vezes a pressão no fundo da Fossa das Marianas, e uma fração de hidrogênio de 12,98% da massa do planeta, o que corresponde a mais de duas massas terrestres de hidrogênio distribuídas por um corpo de 17,15 massas terrestres. O oceano de magma se estende até 3,15 raios terrestres, e o envelope, até o raio de 1 bar, de 3,863 raios terrestres: a rocha derretida ocupa 82% do raio e abriga 16,45 das 17,15 massas terrestres, ou 96% de toda a massa. Na superfície do oceano, onde a mistura cruza a binodal, a temperatura é de 3.256 kelvins; no centro, 10.698 kelvins, quase o dobro da temperatura da fotosfera do Sol. O envelope gasoso de hidrogênio e hélio soma cerca de 4% da massa, perto de 0,7 massa terrestre.

Com essa configuração, o modelo reproduz o raio equatorial em 24.760,9 quilômetros, contra 24.766 medidos, o J₂ em 3.529,50 × 10⁻⁶, contra 3.528,05 × 10⁻⁶ do alvo corrigido pelos ventos, o J₄ em −31,50 × 10⁻⁶, contra −31,95 × 10⁻⁶, o momento de inércia em 0,2381, contra 0,2410, e a luminosidade intrínseca em 2,93 × 10¹⁵ watts, contra 3,30 × 10¹⁵, tudo dentro das incertezas observacionais. A única pedra no sapato é a temperatura a 1 bar: o modelo prevê 84,4 kelvins, 12 graus acima dos 72 ± 2 kelvins medidos por Gunnar Lindal a partir da ocultação de rádio da Voyager.

O qui-quadrado reduzido de 2,6, com três graus de liberdade, equivale a um valor-p de cerca de 0,05 e é causado apenas por esse excesso. Os autores o atribuem às incertezas nas opacidades cinzentas usadas no modelo atmosférico, e mostram que uma diferença de 2,5 vezes na opacidade média de Rosseland, dentro do intervalo plausível para a absorção induzida por colisões entre moléculas de hidrogênio, as bandas do metano e os aerossóis, deslocaria a temperatura em cerca de 15 kelvins, o suficiente para fechar a conta.

A luminosidade de Netuno nasce, nesse modelo, de um gargalo térmico. O calor do oceano de magma precisa atravessar duas barreiras: uma camada-limite finíssima, de 105 metros, logo acima da superfície do oceano, e uma camada estável pelo critério de Ledoux, com cerca de 120 quilômetros de espessura, na qual o gradiente de peso molecular criado pela chuva de silicato impede a convecção. A primeira reduz o fluxo de calor por um fator da ordem de 5; a segunda, por um fator da ordem de 50. É como uma panela fechada por duas tampas: a quantidade de vapor que escapa depende menos do fogo embaixo do que da vedação em cima.

O número de Rayleigh do envelope, que mede quão vigorosamente o fluido ferve, fica em torno de 10³², mas o que o modelo ajusta é a razão entre esse número e seu valor crítico, da ordem de 10¹² para Netuno, e os próprios autores insistem que esse valor é um parâmetro efetivo, não uma previsão de primeiros princípios da teoria da convecção em rotação.

Para Urano, a busca com 1.008 modelos converge para uma binodal a 7,04 gigapascais e uma fração de hidrogênio de 14,36% da massa, valores que incluem a parcela nominal de hélio, já que o código adota um peso molecular de 2,3 gramas por mol para o gás. O ajuste é notavelmente limpo: qui-quadrado reduzido de 0,19. O raio sai em 25.561,7 quilômetros, contra 25.559 observados; o J₂, em 3.522,80 × 10⁻⁶, praticamente idêntico ao alvo; a temperatura a 1 bar, em 75,0 kelvins, contra 76 ± 2; a luminosidade, em 6,65 × 10¹⁴ watts, contra 6,30 × 10¹⁴ das estimativas recentes; o momento de inércia, em 0,2275, contra 0,2250.

Cinco dos seis observáveis caem dentro das incertezas puramente observacionais; a exceção é o J₄, que se afasta do valor corrigido pelos ventos por 2,5 vezes a incerteza de medida, embora permaneça bem dentro da margem de erro do próprio modelo de ventos. Para um planeta que resistiu por décadas a modelos capazes de reproduzir simultaneamente densidade, luminosidade, estrutura térmica e harmônicos gravitacionais, é um resultado digno de nota.

Postos lado a lado, os dois planetas revelam diferenças sutis e coerentes. O oceano de magma de Urano se estende até 3,114 raios terrestres sob um raio total de 3,986, e abriga 13,89 das 14,54 massas terrestres do planeta. Sua densidade média, de 2.537 quilogramas por metro cúbico, é cerca de 12% menor que a de Netuno, de 2.899, porque há mais hidrogênio dissolvido no fundido: 10,41% contra 9,27% em massa. Nas palavras dos autores, Urano é mais “fofo” que Netuno.

A temperatura na superfície do oceano é de 3.071 kelvins, e a do centro, de 8.517, mais baixas que as do vizinho. E a camada estável de Urano é mais espessa: 186 quilômetros no total, contra 129 em Netuno, com a camada fina na base medindo 741 metros, contra 105. Essa tampa mais grossa é a explicação do modelo para a luminosidade menor de Urano, e a equipe faz questão de frisar que se trata de uma explicação, não de uma previsão: a espessura é prescrita pela razão de Rayleigh ajustada, não deduzida da física.

A pergunta sobre por que Urano brilha tão pouco, portanto, segue aberta. O modelo reproduz a luminosidade baixa porque ajusta a largura da camada-limite para isso, não porque a deduz. O que os cálculos fazem é dar peso a uma hipótese que vem ganhando terreno desde o trabalho de Allona Vazan e Ravit Helled em 2020: gradientes de composição que suprimem a convecção no envelope seriam a origem do fraco fluxo de calor. Uma diferença de um fator de aproximadamente 100 na criticalidade da convecção entre os dois planetas é fisicamente plausível, argumentam, se os envelopes profundos diferirem em viscosidade, composição detalhada ou excesso de temperatura superadiabática. E esse fator de 100 é, segundo eles, o resultado robusto do exercício, independente da calibração incerta da teoria de Chandrasekhar, extrapolada aqui por sete ordens de grandeza além do regime em que é conhecida.

As distribuições de probabilidade das buscas, apresentadas nos chamados corner plots, gráficos que mostram as correlações entre cada par de parâmetros, revelam a natureza integrada do modelo. A pressão na binodal e a fração de hidrogênio são fortemente correlacionadas: aumentar uma exige aumentar a outra para preservar o ajuste. O J₄, o harmônico mais sensível à estrutura externa, exibe uma anticorrelação forte com a concentração de hidrogênio no fundido, uma ligação que passa pela temperatura na base do envelope e pelas condições que a binodal impõe ali. Ou seja, um número medido pela deformação do campo gravitacional a milhares de quilômetros de distância responde à química de um oceano de rocha derretida que ninguém jamais verá. Mudanças razoáveis na equação de estado da mistura supercrítica deslocam os valores exatos de pressão e fração de hidrogênio sem degradar a qualidade do ajuste, o que os autores traduzem numa incerteza sistemática absoluta da ordem de 2% na fração de hidrogênio: 9,5% ± 2% para o fundido de Netuno, por exemplo.

A prova mais independente, contudo, vem da química. Para verificar se um oceano de magma em contato com um envelope de hidrogênio produz a atmosfera que se observa, a equipe aplicou o código GlobalChemicalEquilibrium, desenvolvido originalmente para super-Terras e sub-Netunos, às condições da binodal de seus modelos. O cálculo equilibra uma atmosfera rica em hidrogênio com um manto de silicato dominado por MgSiO₃, com fração molar de 0,9211, e traços de MgO, SiO₂, FeO, FeSiO₃, Na₂O e Na₂SiO₃, sem fase metálica separada, já que o ferro se dissolve no silicato. O gás começa com carbono, enxofre e nitrogênio em proporções solares. Quando o sistema atinge o equilíbrio, a mistura acima do oceano é dominada por hidrogênio, com fração molar de 0,9782, e água, com 0,0613; o metano, com 0,0149, é o principal portador de carbono, o H₂S, com 5,08 × 10⁻⁴, o principal portador de enxofre, e as espécies nitrogenadas praticamente desaparecem, com amônia e nitrogênio molecular em meros 4,16 × 10⁻⁹.

Traduzidos para a linguagem das observações, esses números batem com o que a Voyager e os telescópios terrestres mediram. O metano sai enriquecido em cerca de 28 vezes o valor solar em relação ao hidrogênio, na mesma faixa das 40 a 50 vezes de Urano e das 80 de Netuno. O enxofre aparece como H₂S com enriquecimento de 20 vezes, dentro das 20 a 30 vezes inferidas para os dois planetas por Imke de Pater e colegas em 1991.

O nitrogênio some porque, sob as condições redutoras impostas por um envelope espesso de hidrogênio em contato com o magma, ele é altamente solúvel no silicato fundido, um mecanismo estudado por Oliver Shorttle e colaboradores em 2024 e por Werlen e colegas em 2026. A água também se forma acima do oceano, mas fica presa numa armadilha fria: ao subir, a atmosfera satura e condensa a centenas de bar, onde o perfil de temperatura cruza a curva de condensação, e a pressão parcial de vapor fica limitada ao valor de saturação por mais vigorosa que seja a mistura vertical. O metano, com pressão de vapor alta mesmo a temperaturas baixas, sobrevive como gás detectável na região de 1 a 10 bar.

O padrão, metano e enxofre em abundância, nitrogênio ausente, é exatamente o que o Telescópio Espacial James Webb vem encontrando em sub-Netunos distantes. Em K2-18 b, Nikku Madhusudhan e colegas detectaram metano e dióxido de carbono em 2023, sem traço de amônia; em TOI-270 d, Björn Benneke e sua equipe encontraram a mesma assinatura em 2024. A ausência de amônia nesses mundos tem sido lida como a impressão digital de um oceano de magma sob a atmosfera. Se Urano e Netuno exibem a mesma marca, argumentam os autores, os gigantes de gelo do Sistema Solar pertencem à mesma classe de objetos que os anões gasosos caracterizados pelo Webb. Eles reconhecem, com cautela, que a não detecção de amônia nas camadas superiores não prova a existência de um oceano de magma hoje, nem a depleção de nitrogênio no planeta inteiro, como advertiram Cheng Li e colaboradores em 2017 e Tristan Guillot em 2022. Mas é consistente com o que se espera onde há magma.

A conexão com os sub-Netunos aparece também na quantidade de hidrogênio. As frações de hidrogênio inferidas para esses exoplanetas variam numa ordem de grandeza, de cerca de 1% a 10% da massa, segundo a compilação de Jacob Bean e colegas de 2021. Para entender onde Urano e Netuno se encaixam, a equipe montou um experimento numérico simples: planetas de 2, 5, 10, 15 e 20 massas terrestres são construídos pela adição sequencial de embriões sorteados de uma distribuição de massas entre 0,1 massa terrestre e metade da massa final, em discos cuja temperatura vai de 300 kelvins, para os menores, a 50 kelvins, para os maiores. A cada impacto, parte do envelope de hidrogênio preexistente é arrancada, e o embrião traz seu próprio hidrogênio, seguindo a escala de Sivan Ginzburg, Hilke Schlichting e Re’em Sari de 2016, em que a fração acretada cresce com a massa do embrião elevada a 0,8 e diminui com a temperatura do disco.

O resultado é uma distribuição aproximadamente log-normal cuja mediana sobe sistematicamente com a massa: de cerca de 1% de hidrogênio para planetas de 2 massas terrestres a cerca de 9% para os de 20. Dois efeitos empurram a curva. Planetas mais massivos engolem embriões maiores, e cada impacto contribui desproporcionalmente mais hidrogênio por unidade de massa sólida; e, no cenário adotado, eles se formam em regiões mais frias, o que favorece a retenção de gás. Os cerca de 13% de hidrogênio inferidos para Urano e Netuno ficam na extremidade superior da distribuição prevista, exatamente onde planetas de 14 a 17 massas terrestres deveriam estar. Nada no quadro exige um mecanismo de formação distinto para os gigantes de gelo; eles seriam a cauda massiva de uma população única que começa nos sub-Netunos.

O modelo também convive bem com um dos traços mais estranhos dos dois planetas: seus campos magnéticos. A Voyager 2 revelou que, ao contrário dos de Júpiter, Saturno e Terra, os campos de Urano e Netuno não são dominados por um dipolo; são multipolares, inclinados e deslocados do centro, como mostraram Norman Ness e colegas em 1986 e 1989. Estudos de dínamo indicam que a geração de campos multipolares não exige uma geometria específica da região condutora nem uma rotação mais lenta. Simulações de dinâmica molecular de François Soubiran e Burkhard Militzer, de 2018, mostraram que silicatos fundidos nas condições relevantes têm números de Reynolds magnéticos suficientes para sustentar um dínamo, e o trabalho de Dragulet e Stixrude, de 2025, indica que o ferro dissolvido aumenta ainda mais essa condutividade. Os autores apostam que o hidrogênio amplifica o efeito. Não há razão, concluem, para duvidar de que um oceano de magma convectivo gere os campos observados.

Um planeta assim teria mudado muito ao longo da história do Sistema Solar? A equipe calculou a evolução quase estática de Netuno nos últimos 4,5 bilhões de anos, impondo massa, fração de hidrogênio e temperatura de equilíbrio e deixando a pressão na binodal, que está unicamente correlacionada com a energia total e o estado térmico do planeta, evoluir sob o balanço de energia.

As mudanças são modestas: o raio encolhe pouco, a luminosidade intrínseca cai de mais de 4 × 10¹⁵ para os atuais 3 × 10¹⁵ watts, a temperatura central desce de aproximadamente 12.000 para 10.700 kelvins, e a pressão na superfície do oceano sobe de cerca de 3,6 para 5,4 gigapascais à medida que a binodal migra para dentro. A comparação entre a contração do raio e a fração de hidrogênio no fundido mostra que a exsolução do hidrogênio, sua saída do magma para o envelope durante o resfriamento, atua como freio no encolhimento. Eventos irreversíveis e violentos dos primeiros milhões de anos, como impactos gigantes, ficam fora do cálculo por necessidade.

Os autores dedicam uma seção inteira às limitações, e elas merecem atenção, porque são o que separa uma hipótese forte de uma afirmação exagerada. As soluções não são únicas: mudanças razoáveis nos parâmetros da equação de estado do fundido supercrítico produzem ajustes igualmente bons com pressões de binodal e frações de hidrogênio distintas, embora com curvas de pressão versus densidade praticamente iguais. A extrapolação linear da compressibilidade para concentrações de hidrogênio acima dos 4% simulados é admitida como aproximação, e um teste com 10,5% de hidrogênio mostrou que as simulações reproduzem as densidades extrapoladas, só que com 12% de hidrogênio, o que ilustra a incerteza absoluta de cerca de 2% na fração total. O modelo é unidimensional e esférico, com a forma rotacional calculada depois, sem permitir que a densidade relaxe em resposta à deformação, um efeito que os autores estimam ser menor que as incertezas observacionais. E a temperatura a 1 bar de Netuno continua 12 kelvins acima do medido.

Tampouco se trata do único modelo capaz de reproduzir Urano e Netuno. Modelos empíricos de densidade contínua, estruturas com camadas de água e de compostos de carbono, nitrogênio e hidrogênio separados por exsolução, como propôs Burkhard Militzer em 2024, e cenários com separação de fases entre hidrogênio e água, como os de Bailey e Stevenson e os de Cano Amoros e colegas, também se ajustam aos dados. O que a equipe da UCLA reivindica é parcimônia: a física das fases silicato-hidrogênio limita a dois os parâmetros necessários para reproduzir densidade, gravidade externa e momento de inércia, com um terceiro reservado à luminosidade. Modelos mais simples são mais fáceis de testar e menos sujeitos a se ajustar ao ruído, e os autores invocam o critério de informação bayesiano como régua para pesar complexidade contra parcimônia, sem fingir que menos parâmetros equivalem a verdade.

O que está em jogo vai além da nomenclatura. Se Urano e Netuno são anões gasosos crescidos, eles se tornam os únicos exemplares acessíveis da classe de planeta mais comum da galáxia, os sub-Netunos, que o telescópio Kepler revelou aos milhares e que o Webb começa a dissecar por espectroscopia. Dois laboratórios a 20 e 30 unidades astronômicas, com massas, raios, campos gravitacionais e luminosidades medidos, valem mais que qualquer estatística de exoplanetas para calibrar equações de estado e propriedades de materiais a pressões de centenas de gigapascais. A recíproca também vale: cada espectro de sub-Netuno sem amônia e com metano é um teste a mais para o modelo dos gigantes de gelo. A fronteira entre a ciência planetária do Sistema Solar e a dos exoplanetas, que durante décadas tratou os dois gigantes azuis como um caso à parte, se dissolve na mesma solução supercrítica que o modelo descreve.

As respostas definitivas dependerão de voltar lá. O relatório decenal de ciência planetária dos Estados Unidos, publicado em 2022, elegeu uma sonda orbital de Urano, acompanhada de uma cápsula atmosférica, como a missão de grande porte prioritária do período de 2023 a 2032, mas o projeto segue como conceito, ainda sem aprovação formal, com a janela de lançamento mais favorável no início da década de 2030 e chegada apenas na década de 2040.

Uma órbita prolongada permitiria medir harmônicos gravitacionais de ordem mais alta, como o J₆, necessários para separar de vez a contribuição dos ventos, determinar o momento de inércia pela precessão, mapear o campo magnético em detalhe e sondar a atmosfera profunda abaixo dos níveis que o sensoriamento remoto alcança. O modelo de oceano de magma faz previsões testáveis por esses instrumentos: um interior inteiramente fluido, sem núcleo sólido distinto, uma camada estável contra a convecção de 100 a 200 quilômetros sobre o magma, e uma atmosfera profunda enriquecida em metano e sulfeto de hidrogênio mas pobre em nitrogênio. Até lá, o laboratório mais próximo fica a quase 3 bilhões de quilômetros, no caso de Urano, e o mais distante, a mais de 4 bilhões, no de Netuno.

Há uma ironia histórica em tudo isso. Urano foi o primeiro planeta descoberto com um telescópio, por William Herschel, em 1781; Netuno foi previsto no papel antes de ser visto, em 1846, a partir das perturbações na órbita de Urano. Os dois nasceram para a ciência como triunfos da inferência sobre a observação direta, e o rótulo de gigantes de gelo, dado dois séculos depois, foi mais uma inferência, dessa vez cosmoquímica, que nunca teve uma sonda orbital para confirmá-la ou refutá-la. O que Young e seus colegas propõem é que, sob os azuis serenos que a Voyager fotografou, exista um mar de rocha derretida e hidrogênio dissolvido, aquecido a milhares de graus e revolto desde o nascimento do Sistema Solar, chovendo silicato sobre si mesmo há 4,5 bilhões de anos. Se estiverem certos, os dois planetas mais negligenciados da nossa vizinhança guardam, bem abaixo das nuvens de metano, a pedra de Roseta para ler os mundos mais comuns da galáxia. E o gelo que os batizou nunca esteve lá.