Daqui a cerca de sete bilhões de anos, o Sol não será mais a esfera amarela e estável que conhecemos. Inchado a ponto de engolir Mercúrio e, muito provavelmente, Vênus, avermelhado, com a superfície fervendo em células de convecção do tamanho de estrelas inteiras, ele pulsará como um coração cansado. De tempos em tempos, uma pluma de gás se soltará dessa superfície e partirá para o espaço, carregando algo como um décimo de milésimo da massa solar. A cada despedida dessas, a estrela sentirá um empurrão minúsculo na direção contrária, de alguns metros por segundo. Repetido milhares de vezes, em direções escolhidas pelo acaso, esse empurrão vai se acumulando. E lá fora, Júpiter, Saturno, Urano e Netuno, que atravessaram incólumes bilhões de anos de história, começarão a sentir suas órbitas vagar.

A pergunta sobre se os planetas são estáveis é mais antiga que a própria mecânica celeste moderna. Isaac Newton, ao perceber que Júpiter e Saturno se perturbam mutuamente, com um lento desvio das órbitas que Johannes Kepler já havia notado, desconfiou de que a ordem planetária acabaria se desfazendo sem uma intervenção divina ocasional. Pierre-Simon de Laplace e Joseph-Louis Lagrange responderam a ele no século XVIII, demonstrando que as perturbações seculares deixam os semieixos maiores das órbitas invariantes e apenas fazem as excentricidades oscilar. A explicação da chamada Grande Desigualdade entre Júpiter e Saturno, descrita por Laplace como uma modulação periódica quase ressonante, virou o emblema do resgate da estabilidade do Sistema Solar. Henri Poincaré mostrou depois que o resgate jamais poderia ser completo, e a era dos computadores deu o veredito final em forma estatística.

Esse veredito, até agora, era tranquilizador. A órbita de Mercúrio carrega uma probabilidade de cerca de 1% de se desestabilizar antes que o Sol deixe a sequência principal, com a relatividade geral desempenhando um curioso papel estabilizador. Para os planetas externos, a resposta parecia ainda mais segura: os quatro gigantes são fracamente caóticos, mas estáveis em escalas de tempo que ultrapassam em muito a idade do Universo. Em 1999, Norman Murray e Matthew Holman estimaram que a difusão caótica levaria cerca de 10¹⁸ anos para desestabilizar Urano, algo como setenta milhões de vezes os 13,8 bilhões de anos transcorridos desde o Big Bang. Em horizontes assim, a conclusão lógica era que fatores externos agiriam primeiro, e que o fim do Sistema Solar exterior seria remoto e infligido de fora.

Um artigo publicado em 1º de outubro de 2026 no periódico The Astrophysical Journal Letters derruba essa conclusão pela base. Assinado por Konstantin Batygin e Jim Fuller, ambos do Instituto de Tecnologia da Califórnia, e por Fred C. Adams, da Universidade de Michigan, o trabalho tem um título que não deixa margem para dúvida: “Instabilidade terminal do Sistema Solar desencadeada pela perda de massa estocástica do Sol”. A tese central é que as estimativas anteriores repousavam sobre uma premissa nunca examinada, a de que o Sol perderá massa de maneira suave e contínua. Quando essa premissa é trocada por uma perda de massa em parcelas discretas, compatível com as observações mais recentes de anãs brancas, a vida dinâmica dos planetas gigantes despenca de 10¹⁸ anos para aproximadamente um bilhão de anos depois da formação da anã branca. O algoz não é nenhuma estrela intrusa. É o próprio Sol.

Os números das simulações são contundentes. Nas realizações que adotam a granularidade de perda de massa mais compatível com os dados, cruzamentos de órbitas aparecem ainda na fase de gigante vermelha, quando o Sol conserva 89% de sua massa atual. Cerca de 40% dos sistemas simulados sofrem ruptura ou espalhamento violento antes mesmo que a anã branca se forme, e perto de 90% se autodestroem dentro de três bilhões de anos. Saturno pode ser ejetado em poucos milhões de anos, Urano e Netuno podem trocar de lugar, periélios podem ser empurrados para dentro da órbita de Júpiter. Segundo os autores, as consequências não são uma correção perturbativa à resposta antiga, e sim a sua substituição completa: a duração do Sistema Solar exterior passa a ser ditada pela granularidade da perda de massa de sua estrela.

Para entender o que muda, é preciso lembrar o que o Sol fará nos próximos bilhões de anos. Com cerca de 4,6 bilhões de anos de idade, a estrela tem hidrogênio no núcleo para mais uns cinco ou seis bilhões de anos. Depois disso, o núcleo se contrai, as camadas externas se expandem e o Sol se transforma numa gigante vermelha; mais tarde, passa pelo chamado ramo assintótico das gigantes, uma fase de pulsações e ventos intensos, até expulsar por completo o envelope e deixar para trás uma anã branca, o caroço degenerado e quente que esfriará lentamente por toda a eternidade. No modelo estelar usado no estudo, calculado com o código MESA, o Sol termina com 0,538 massa solar, ou seja, perde 46% do que tem hoje. Toda a arquitetura planetária precisa sobreviver a essa dieta brutal.

Quando a perda de massa é lenta e suave, a mecânica celeste sabe muito bem o que acontece. A atração do Sol enfraquece aos poucos, e cada planeta migra para fora de forma proporcional, como um peso preso a um elástico que vai afrouxando devagar. Os semieixos maiores crescem pelo mesmo fator para todos, cerca de 1,85 vezes, e as razões entre os períodos orbitais ficam congeladas; Júpiter e Saturno, por exemplo, mantêm a razão de 2,481 entre seus períodos, medida no estudo por um indicador de movimento médio livre de oscilações espúrias. As excentricidades e inclinações preservam amplitudes e fases, e os ângulos ressonantes de todas as comensurabilidades próximas continuam a circular. O Sistema Solar pós-sequência principal, nesse cenário, é uma cópia uniformemente ampliada do atual, e igualmente estável.

Com as órbitas quase duas vezes maiores, o sistema passa a oferecer uma seção de choque maior a estrelas que passam nas vizinhanças, e foi com base nisso que Jon Zink, Batygin e Adams argumentaram, em 2020, que encontros estelares desmontariam o conjunto de planetas gigantes em algo entre 30 e 100 bilhões de anos. Era uma escada de escalas de tempo cada vez mais longas: 10¹⁸ anos para a instabilidade intrínseca, cem bilhões de anos para a erosão por sobrevoos. A morte do Sistema Solar exterior, nessa visão, viria de fora e tarde, depois que a Via Láctea e Andrômeda provavelmente já tivessem se encontrado e que a maior parte das estrelas que hoje vemos tivesse se apagado.

O que derrubou a escada foi um resultado observacional que, à primeira vista, nada tem a ver com planetas. A astrometria do satélite Gaia, aplicada a pares de estrelas em órbitas muito largas, revelou que anãs brancas recém-nascidas recebem um kick, um chute de velocidade, de cerca de 0,75 quilômetro por segundo no momento em que se formam. O efeito foi medido por Kareem El-Badry e Hans-Walter Rix em 2018 e confirmado por Hsiang-Chih Hwang e Nadia Zakamska em 2025, e é corroborado pela escassez e pela distribuição espacial estendida de anãs brancas em aglomerados estelares, onde o chute arremessa parte delas para longe da região central. Uma anã branca não deveria sair do lugar ao nascer. Se sai, algo a empurrou.

A explicação que vem ganhando corpo é que o envelope de uma gigante vermelha não se despede num vento contínuo e simétrico, mas em uma sequência de ejeções discretas e assimétricas. Fisicamente, essas ejeções seriam plumas lançadas pela convecção assistida por pulsações, estruturas que simulações tridimensionais preveem e que imagens interferométricas de estrelas evoluídas resolvem como nuvens de poeira individuais em evolução. Cada pluma que parte numa direção dá à estrela um recuo na direção oposta, proporcional à fração de massa ejetada multiplicada pela velocidade de escape. Como as direções são aleatórias, os recuos não se somam de forma ordenada: eles executam um random walk, um passeio aleatório, cujo deslocamento líquido cresce com a raiz quadrada do número de passos até atingir a amplitude observada pelo Gaia.

A imagem clássica do passeio aleatório é a de alguém parado sobre um lago congelado que arremessa pedras em direções sorteadas. Cada arremesso o empurra um pouco para trás; depois de mil arremessos, ele não está onde começou, embora nenhum deles tenha sido grande. O ponto crucial, para os planetas, é que cada chute impulsivo no Sol desloca instantaneamente todas as órbitas heliocêntricas. Ventos suaves, mesmo que anisotrópicos, se diluem na média ao longo das órbitas e não deixam marca; impulsos descorrelacionados entre si, ao contrário, escapam dessa diluição secular. É o mesmo mecanismo pelo qual forçamentos estocásticos reestruturam ressonâncias em discos turbulentos ou na migração de Netuno empurrada por planetesimais, com a diferença de que agora a fonte de ruído é a própria estrela.

Alexander Stephan e Keivan Stassun já haviam mostrado, em 2026, que o chute observado, aplicado como um único impulso terminal, é suficiente para desestabilizar sistemas de múltiplos planetas em órbitas largas. O novo trabalho dá um passo além: decompõe esse impulso em seus constituintes físicos, os milhares de recuos individuais, e com isso transfere o início da instabilidade para dentro da própria época de perda de massa. Nenhuma das análises clássicas de estabilidade, dos teoremas de invariabilidade de Laplace e Lagrange às estimativas de erosão por sobrevoos, havia contemplado um Sol que chacoalha. A ligação entre a física estelar e a dinâmica planetária, aqui, deixa de ser um detalhe e passa a ser o eixo do problema.

Para testar a hipótese, os autores recorreram a ferramentas padrão da comunidade: o pacote de N corpos REBOUND, com o integrador adaptativo de alta ordem IAS15 durante a fase de perda de massa e o integrador híbrido TRACE, com passo de dois anos, para as integrações de bilhões de anos. A evolução do Sol veio de uma trajetória calculada com o código MESA para uma estrela de uma massa solar e metalicidade de 2%, com as mesmas prescrições de vento usadas no trabalho de 2020. O Sol e os quatro gigantes foram inicializados a partir dos vetores de estado do sistema JPL Horizons, e não de elementos orbitais médios tabelados, um cuidado que faz diferença: a configuração reproduz a razão de movimentos médios observada entre Saturno e Júpiter, 2,4833, e um período de circulação da Grande Desigualdade de 976 anos, coerente com a modulação de 900 a 1.000 anos do par real.

O primeiro passo foi estabelecer a linha de base. Os quatro gigantes foram integrados através da perda de massa suave em três níveis de realismo: um decaimento linear ao longo de mil a um milhão de anos, perfis em dois estágios terminando num supervento rápido, e a trajetória completa do MESA da sequência principal até a anã branca de 0,538 massa solar. Em mais de sessenta integrações cobrindo todos esses perfis e todas as fases seculares, o resultado foi qualitativamente idêntico: expansão uniforme, razões de período congeladas, estrutura secular preservada, estabilidade por vários bilhões de anos. Nesse processo, os autores revisitaram a captura de Júpiter e Saturno na ressonância 5:2 relatada no artigo de 2020, do qual Batygin e Adams eram coautores, e concluíram que ela não pode ser reproduzida: trata-se, ao que tudo indica, de um artefato no tratamento de coordenadas na passagem entre duas fases daquela simulação. Corrigir o próprio trabalho é parte do método científico, e aqui isso foi feito às claras.

Tudo muda quando o envelope parte em pedaços. Seguindo o arcabouço de ejeção episódica proposto por Fuller em 2026, a perda de massa foi modelada como uma sequência de eventos discretos distribuídos ao longo da trajetória do MESA em proporção à taxa local de perda de massa, cada um retirando do Sol uma parcela fixa e aplicando um recuo impulsivo em direção isotrópica, com módulo igual à fração de massa ejetada vezes a velocidade de escape da estrela naquele instante. Cada impulso dura um ou dois anos, muito menos que qualquer período orbital. A maior parte dos eventos se concentra perto das pontas do ramo das gigantes vermelhas e do ramo assintótico, onde a perda de massa é mais intensa e os chutes individuais ficam entre 5 e 9 metros por segundo; os maiores, de cerca de 30 metros por segundo, vêm da minoria de eventos ocorridos nas fases compactas, quando a velocidade de escape é maior.

A discretude é indispensável, e os autores fizeram questão de provar isso com experimentos de controle. Em simulações com ejeções cronometradas ao acaso, mas esfericamente simétricas, sem recuo algum, a evolução se reduz ao caso adiabático e o sistema permanece estável por três bilhões de anos. O mesmo ocorre quando os decrementos de massa são aplicados instantaneamente ou atrasados pelo tempo de travessia da casca ejetada. Tanto a discretude temporal quanto a assimetria de cada evento são, portanto, estritamente necessárias para a excitação em larga escala. Quem desestabiliza os planetas não é o fato de o Sol emagrecer, e sim o fato de emagrecer aos solavancos, em direções que ninguém escolhe.

A grade de simulações cobre sete valores de massa por ejeção, espaçados logaritmicamente entre 10⁻⁷ e 10⁻⁴ massa solar, com 48 sementes aleatórias para cada valor e duas variantes: numa, as ejeções ocorrem ao longo de todas as fases de gigante; na outra, apenas no ramo assintótico, quando a massa do Sol já caiu abaixo de 0,75 massa solar, que é a fase em que o comportamento episódico está mais bem calibrado. São 672 realizações no total. Ao longo dessa faixa, a velocidade acumulada de chute da estrela vai de 0,018 a 0,58 quilômetro por segundo; no caso fiducial de 10⁻⁴ massa solar, fica entre 0,33 e 0,58 quilômetro por segundo, abaixo do valor nominal de 0,75, mas dentro da faixa permitida pela calibração do Gaia, que vai de cerca de 0,25 a 4 quilômetros por segundo, com progenitoras de massa solar esperadas na ponta inferior, em torno de 0,5.

No caso fiducial, o Sol recebe cerca de 4.600 recuos de alguns metros por segundo, com mediana próxima de 7 metros por segundo, a velocidade de um ciclista apressado. Individualmente são minúsculos diante dos 13 quilômetros por segundo com que Júpiter percorre sua órbita; coletivamente são transformadores. Como recuos sucessivos atingem os planetas em fases orbitais descorrelacionadas, cada elemento orbital executa seu próprio passeio aleatório, e a dispersão acumulada da razão de períodos cresce com a raiz quadrada do número de eventos multiplicada pela massa de cada um, o que equivale à raiz quadrada da massa por ejeção. A grade de 672 realizações confirma essa lei de escala ao longo de três ordens de grandeza, uma concordância que dá confiança de que o modelo captura a física essencial.

Os efeitos sobre a razão de períodos de Saturno e Júpiter dão a medida do estrago. No valor fiducial da massa por ejeção, a dispersão final dessa razão, cujo valor não perturbado é 2,481, chega a 0,23 quando as ejeções ficam restritas ao ramo assintótico e a cerca de 0,6 quando operam por todas as fases de gigante. É uma ordem de grandeza além de qualquer coisa que o cenário adiabático permita. Enquanto isso, a excentricidade de Netuno é bombeada até valores próximos de 0,5, o que significa uma órbita cujo afélio fica três vezes mais longe que o periélio, e as inclinações mútuas se difundem por vários graus. As trajetórias das razões de período, acompanhadas evento a evento, deixam de ser linhas retas e viram feixes que se abrem como um leque.

A arquitetura concêntrica e bem-comportada que todas as análises de estabilidade anteriores tomavam como condição inicial é substituída por uma distribuição ampla de configurações excêntricas, reordenadas e comprimidas. Em alguns estados finais, os quatro planetas ainda desenham anéis quase circulares, mas carregam instabilidades que só vão eclodir centenas de milhões ou bilhões de anos depois; em outros, Urano e Netuno já aparecem com as posições trocadas no exato momento em que a anã branca se forma, ou com órbitas que se cruzam. Um sistema visualmente intacto, no instante da morte da estrela, não é garantia de nada. Essa é uma das lições mais inquietantes do estudo: a aparência de ordem pode esconder um relógio de contagem regressiva.

A consequência central do forçamento estocástico é temporal. O início da instabilidade, antes relegado a 10¹⁸ anos pela dinâmica intrínseca ou a centenas de bilhões de anos pelos sobrevoos, irrompe durante a própria evolução do Sol como gigante. Na variante fiducial que estende as ejeções a todas as fases, configurações de cruzamento de órbitas surgem antes que a perda de massa esteja completa em 37 das 48 realizações, ou 77%, com os primeiros cruzamentos ocorrendo quando o Sol ainda é uma gigante vermelha de 0,89 massa solar. Quando a anã branca se forma, os cruzamentos já progrediram para ruptura completa ou espalhamento violento em 40% dos casos, 19 em 48. A instabilidade, nessas realizações, não é apenas acelerada pela morte da estrela: é contemporânea a ela.

As histórias individuais que as simulações contam variam de dramáticas a lentas. Numa delas, as órbitas se cruzam enquanto o Sol ainda está perdendo massa, com Urano e Netuno se entrelaçando antes mesmo do fim da expansão orbital. Em outra, Saturno é espalhado e ejetado em apenas 2,5 milhões de anos após a formação da anã branca, lançado para uma órbita de dez mil unidades astronômicas com excentricidade tendendo a um, enquanto Urano e Netuno ficam para trás travados numa ressonância 2:1 protegida por fase. Com massa por ejeção dez vezes menor, de 10⁻⁵ massa solar, os cruzamentos entre Urano e Netuno começam por volta de 100 milhões de anos e Netuno é ejetado aos 200 milhões; na variante restrita ao ramo assintótico, os dois planetas trocam de ordem aos 1,24 bilhão de anos e Netuno é expulso aos 1,32 bilhão.

A variante restrita ao ramo assintótico, livre de chutes durante o ramo das gigantes vermelhas por construção, não produz nenhuma ruptura nessa fase, mas converge para o mesmo destino. Quando os 672 estados finais são integrados por mais três bilhões de anos sem qualquer perturbação externa, a fração de sistemas intrinsecamente instáveis sobe de menos de 2%, para massas por ejeção de até 3 × 10⁻⁷ massa solar, até 96% na variante completa e 83% na variante restrita, no valor fiducial. Os instantes de início da instabilidade se espalham de poucos milhões de anos a mais de 2,5 bilhões, com distribuição ampla em escala logarítmica, e a mediana do tempo até a instabilidade na variante restrita fica em torno de 40 milhões de anos.

A vulnerabilidade de cada sistema está escrita na sua arquitetura, e a chave de leitura é a separação mínima entre planetas vizinhos medida em raios de Hill mútuos. O raio de Hill é o tamanho da região em torno de um planeta onde a gravidade dele vence a do Sol, e a distância entre órbitas vizinhas contada nessa unidade é o termômetro clássico da estabilidade. Abaixo de três raios de Hill, a sobrevivência cai para cerca de 3%; entre três e cinco há uma transição ampla; acima de cinco, mais de 85% dos sistemas sobrevivem. O Sistema Solar atual, sem chutes, está em 5,5, confortavelmente do lado seguro. As estrelas que passam, pilar da explicação antiga, simplesmente perdem a corrida: o sistema chutado se destrói por causas internas em escalas de tempo cem vezes mais curtas que a espera prevista pelo primeiro sobrevoo disruptivo.

O balanço completo da grade ajuda a dimensionar o fenômeno. Das 672 realizações, 448 permanecem estáveis e 224 se desestabilizam; entre essas, o primeiro sinal registrado é cruzamento de órbitas em 151 casos, espalhamento em larga escala em 49 e ejeção direta em 24. Para resolver o que acontece depois do primeiro cruzamento, os autores seguiram cada um desses sistemas por mais 100 milhões de anos: 113 dos 151 ejetam pelo menos um planeta nessa janela, 30 evoluem para espalhamento de grande amplitude e apenas 8 continuam formalmente em cruzamento sem evolução macroscópica. Nenhum deles se acomoda numa configuração protegida por fase, sem cruzamentos. Uma vez que as órbitas se tocam, o caminho de volta não existe.

Nem todo rearranjo é destrutivo. O passeio aleatório tem um efeito colateral sutil: a captura ocasional em ressonâncias de movimento médio, algo inatingível no limite adiabático. Uma ressonância desse tipo ocorre quando os períodos de dois planetas guardam uma razão de números inteiros simples, como 2:1 ou 5:2, de modo que os encontros se repetem sempre nos mesmos pontos das órbitas, como quem empurra um balanço no instante certo. A captura vira um problema de primeira passagem: o passeio precisa entrar no domínio ressonante devagar o bastante para ficar preso, em vez de atravessá-lo. A probabilidade, por isso, não cresce monotonicamente com a massa por ejeção; ela atinge um pico onde a dispersão acumulada é comparável à soma da distância da ressonância à razão inicial com a largura do domínio.

Dois canais operam na mesma faixa de massas por ejeção: a ressonância 5:2 entre Júpiter e Saturno, com probabilidade máxima de captura de 16,7%, e a 2:1 entre Urano e Netuno, com 20,8% na variante restrita ao ramo assintótico, esta favorecida por uma geometria que se retroalimenta, pois os chutes bombeiam a excentricidade de Netuno e assim alargam o próprio domínio do qual o par se aproxima. Os pares capturados libram com amplitudes próximas da separatriz, fixados na comensurabilidade com uma razão de períodos de 2,000 ± 0,002. Em oito das 672 realizações, os dois pares ficam presos simultaneamente numa configuração duplamente ressonante, 5:2 mais 2:1, que o Sistema Solar jamais possuiu e que, no regime de ejeções suaves, persiste pelos três bilhões de anos integrados. Cadeias ressonantes, conhecidas como fósseis de captura dissipativa em sistemas jovens como TRAPPIST-1, podem evidentemente também ser montadas pelo ruído em torno de uma estrela moribunda.

A ressonância, contudo, não é santuário. No regime granular, os sistemas capturados sobrevivem um pouco menos que os não ressonantes, 13 em 40 contra 162 em 344, numa tendência de 1,8 desvio padrão, porque a captura estocástica sem dissipação deposita os sistemas perto da fronteira caótica do domínio ressonante, onde são instáveis mesmo sem perturbação externa alguma. O quadro geral se organiza em três regimes. Massas por ejeção de 10⁻⁷ massa solar ou menores preservariam em grande parte a réplica adiabática congelada; valores intermediários, entre 10⁻⁷ e 10⁻⁶, rendem uma loteria de 15% a 20% de captura ressonante num sistema apenas moderadamente excitado; valores de 10⁻⁵ massa solar para cima produzem arquiteturas quentes que se desmontam dentro de bilhões de anos ou, na calibração fiducial, enquanto o Sol ainda brilha como gigante.

A evidência disponível favorece o ramo violento. Sob a interpretação episódica, mesmo o limite inferior conservador da faixa de chutes permitida, 0,25 quilômetro por segundo, corresponde a uma massa por ejeção de cerca de 2 × 10⁻⁵ massa solar na variante completa, ou 6 × 10⁻⁵ na restrita, já dentro do regime violento. O valor nominal de 0,75 quilômetro por segundo corresponde a algo entre 1,7 e 5 × 10⁻⁴ massa solar, no topo da grade ou acima dele. Tomada ao pé da letra, a distribuição observada de chutes de anãs brancas fixa, portanto, a escala de tempo da dissolução do Sistema Solar. As contagens das simulações são eloquentes: no valor fiducial, 46 de 48 realizações da variante completa e 40 de 48 da restrita se desestabilizam, e uma população substancial sequer sobrevive à perda de massa.

Resultados tão fortes exigem testes de robustez, e o artigo os faz em série. A primeira pergunta é se tratar os recuos como instantâneos é legítimo, já que o lançamento de uma pluma leva um ou dois anos e a parcela ejetada continua exercendo gravidade sobre cada planeta até cruzar sua órbita, o que a 30 quilômetros por segundo leva de 1,6 a 9 anos. Os autores repetiram baterias de simulações espalhando cada recuo por 3, 10 e 30 anos. Com três anos, a duração fisicamente motivada, nenhuma dispersão muda além do ruído estatístico, em acordo com correções previstas de 1% a 4%. No extremo de 30 anos, o único canal ao qual a teoria de transferência por planeta se aplica diretamente, o passeio da razão de períodos de Júpiter e Saturno, mostra exatamente a supressão prevista, 0,39 medido contra 0,48 previsto, enquanto as excentricidades seguem intactas, porque o déficit de momento angular é absorvido sobretudo por Urano e Netuno, cujos períodos de 84 a 570 anos tornam qualquer evento plausível profundamente impulsivo.

Outros controles delimitam a geometria da ejeção. Confinar todos os chutes ao plano equatorial da estrela, a configuração que maximiza o forçamento no plano das órbitas, aumenta as dispersões de excentricidade pelo fator esperado de cerca de 1,22, com 1,27 medido para Netuno, sem alterar a dispersão das razões de período nem a fração de instabilidade em três bilhões de anos, 23 em 24 contra 22 em 24. As fases compactas da estrela, com raio menor que 30 raios solares, contribuem com 42% da difusão total, mas suprimir os recuos nessas fases não elimina o canal imediato de instabilidade: os cruzamentos durante a perda de massa caem de 19 para 15 em 24 realizações, diferença sem significância estatística. Achatamento solar, precessão relativística, marés, planetas terrestres, maré galáctica e sobrevoos foram omitidos, por serem irrelevantes entre 5 e 60 unidades astronômicas nas escalas de tempo consideradas.

Uma comparação controlada esclarece o que a decomposição do chute acrescenta. Quando o vetor líquido acumulado de cada realização, com módulo de 0,5 ± 0,25 quilômetro por segundo, é aplicado como um único impulso ao final da perda de massa suave, os sistemas alcançam 83% de instabilidade em três bilhões de anos, contra 92% para os passeios aleatórios. A instabilidade bruta está, assim, quase saturada nos dois cenários a essas amplitudes, mas metade das instabilidades do passeio aleatório ocorre antes que a anã branca sequer se forme, algo que um chute único no fim jamais produziria. A fração de instabilidades imediatas, em menos de 10 milhões de anos, no cenário de chute único é de 46%, bem próxima dos cerca de 47% encontrados por Stephan e Stassun, o que sugere coerência entre as duas abordagens.

As implicações vão muito além do Sistema Solar, porque o raciocínio se estende a sistemas planetários em torno de toda estrela destinada a virar anã branca, o que inclui quase todas elas. Análogos verdadeiros de Júpiter e Saturno são intrinsecamente raros, com ocorrência de planetas gigantes entre 10% e 15% das estrelas, mas planetas da classe dos gigantes de gelo parecem estar entre os desfechos mais comuns da formação planetária. E a suscetibilidade à ruptura por chutes cresce com a raiz quadrada da distância orbital e diminui com a raiz cúbica da massa do planeta, para um mesmo espaçamento em raios de Hill. Isso decorre de dois fatos: o recuo da estrela perturba a órbita de todos os planetas igualmente, independentemente da massa deles, e órbitas mais fracamente ligadas são mais fáceis de excitar, enquanto a excentricidade necessária para alcançar o vizinho escala com o raio de Hill, proporcional à raiz cúbica da massa.

Os sistemas de menor massa, que dominam o censo de exoplanetas, são por isso mais frágeis que os gigantes gasosos do Sistema Solar a qualquer distância orbital. Agrava o quadro o fato de os sistemas planetários parecerem se auto-organizar rumo à beira da estabilidade, um princípio insinuado pelo modelo de Nice para o Sistema Solar exterior, pelo cenário de quebra de cadeias ressonantes nos sistemas compactos de super-Terras e pelo espaçamento marginal em raios de Hill medido nos sistemas de múltiplos planetas conhecidos. Para arquiteturas assim, a perda de massa estocástica fornece um empurrão final universal: uma estabilidade marginal mantida por bilhões de anos na sequência principal não garante nada diante de uma morte granular. Os estados dinâmicos dos sistemas planetários em torno de anãs brancas tornam-se, assim, um registro fóssil da granularidade da perda de massa terminal de suas estrelas.

Um caso concreto ilustra o ponto. Em 2020, Andrew Vanderburg e colaboradores descobriram um planeta gigante numa órbita compacta em torno da anã branca WD 1856+534, a cerca de 80 anos-luz da Terra, completando uma volta a cada 1,4 dia. Para os autores do novo estudo, a interpretação mais natural é um espalhamento interno entre planetas, que pode entregar um gigante ao raio de circularização por marés. De forma mais ampla, os diagnósticos de conjunto são monotônicos na massa por ejeção, dispersão de razões de período, excitação de excentricidade e fração de destruição crescem todas com ela, enquanto a ocupação de ressonâncias atinge o pico em valores intermediários. As estatísticas dos sistemas planetários de anãs brancas poderão um dia, portanto, restringir a granularidade da perda de massa por dois lados.

Uma segunda consequência é que sistemas chutados perdem planetas em profusão, começando poucos milhões de anos após a formação da anã branca, e com isso alimentam uma população inicialmente morna de planetas errantes que se movem no espaço junto com anãs brancas jovens. Esses mundos sem estrela, que levantamentos de microlentes gravitacionais vêm detectando desde a década de 2010, teriam ganhado uma nova fábrica. Nas simulações do Sistema Solar, Netuno e Saturno aparecem entre os expulsos. É uma ideia estranha e poderosa: um dia, um gigante gelado nascido do Sol pode estar vagando sozinho pela Via Láctea, sem lembrança da estrela que o criou.

A terceira consequência conecta o resultado a um dos enigmas mais estudados das anãs brancas: a poluição de suas atmosferas por elementos pesados, que deveriam afundar em pouco tempo sob a gravidade intensa dessas estrelas e que, mesmo assim, aparecem numa fração considerável delas, exigindo um suprimento contínuo de material. As rupturas retardadas das simulações fornecem exatamente a agitação dinâmica tardia que as estatísticas de poluição pedem, com inícios distribuídos amplamente em escala logarítmica ao longo dos três bilhões de anos de idade de resfriamento simulados. E fazem isso sem recorrer às arquiteturas primordiais marginalmente empacotadas que os modelos de instabilidade por empacotamento precisam supor, condição que o Sistema Solar exterior real, com seu espaçamento atual, nem sequer satisfaz. Os chutes fabricam instabilidade onde antes não havia como obtê-la.

O desfecho dos cruzamentos de órbita é governado por um parâmetro chamado número de Safronov, o quadrado da razão entre a velocidade de escape da superfície do planeta e sua velocidade orbital. Quando esse número é grande, um encontro próximo tende a arremessar o corpo para fora do sistema antes que uma colisão aconteça; quando é pequeno, as colisões dominam. Para os gigantes do Sistema Solar expandido, o número de Safronov fica entre 34 e 72, de modo que ejeções prevalecem sobre colisões, como de fato ocorre nas integrações, alimentando a população de planetas errantes. Sistemas compactos, com valores da ordem de um ou menores, canalizariam suas instabilidades para colisões e trajetórias que roçam a estrela, justamente a rota que entrega material à anã branca. Dentro de um mesmo sistema, a lógica favorece colisões para os planetas de menor massa e ejeções para os gigantes, o que prevê uma poluição alimentada preferencialmente por corpos rochosos, em concordância com as composições predominantemente rochosas inferidas nas anãs brancas poluídas.

Nada disso dispensa cautela, e os próprios autores listam os pontos frágeis. A massa de cada ejeção é incerta por uma ordem de grandeza, e as observações restringem o chute líquido acumulado, não a granularidade que o produz; a interpretação episódica é a mais natural, mas é uma interpretação. A extensão do comportamento episódico ao ramo das gigantes vermelhas é uma hipótese, e as duas variantes da grade servem exatamente para delimitar suas consequências. Idealizações do modelo, como parcelas de massa idêntica, acoplamento total do momento à velocidade de escape e a omissão da interação gravitacional direta entre pluma e planeta, limitada a menos de 25% da contribuição do recuo, são absorvidas na incerteza de ordem de grandeza coberta pela grade. E trajetórias caóticas só podem ser interpretadas de forma distribucional: reduzir o passo de integração à metade reproduz 11 de 12 veredictos de estabilidade, com o único caso discrepante refletindo a não convergência esperada de órbitas individuais.

Quem lê tudo isso pensando na Terra deve ter em mente que os planetas terrestres ficaram fora das simulações, por serem irrelevantes para a dinâmica entre 5 e 60 unidades astronômicas. O destino deles, de todo modo, se decide antes: Mercúrio e muito provavelmente Vênus serão engolidos pela gigante vermelha, e a sobrevivência da Terra ao sopro final do Sol é uma questão disputada há décadas. Muito antes disso, daqui a cerca de um bilhão de anos, o aumento gradual do brilho solar já terá evaporado os oceanos. O que o novo estudo muda não é o futuro da vida na Terra, e sim o retrato do que restará do Sistema Solar como estrutura quando a estrela que o organizou tiver se apagado.

Colocadas lado a lado, as escalas dizem tudo. A estabilidade do Sistema Solar exterior era uma escada de degraus cada vez mais altos: 10¹⁸ anos para a difusão caótica intrínseca, cerca de cem bilhões de anos para a erosão por estrelas passantes. O novo trabalho acrescenta um degrau na base, dramaticamente mais curto, cuja altura é fixada por uma única grandeza da física estelar, a granularidade da perda de massa. A vida dinâmica dos gigantes passa a ser contada em um bilhão de anos depois da formação da anã branca, um intervalo menor que o tempo que o Sol já levou para chegar até aqui. O fim deixa de ser um evento remoto e externo e volta ao território astrofísico familiar.

O artigo termina onde o assunto começou. Newton, contemplando as órbitas errantes de Júpiter e Saturno, desconfiou de que a ordem planetária era mortal, e três séculos de mecânica celeste, dos teoremas de invariabilidade às integrações numéricas de bilhões de anos, adiaram progressivamente esse veredito, mais recentemente para muito além da idade do Universo. Os resultados de Batygin, Fuller e Adams trazem a dissolução do Sistema Solar de volta para perto e trocam o culpado: não o lento gotejar do caos, nem o encontro casual com uma estrela de passagem, mas o próprio Sol, que ao morrer não apenas amplia o sistema planetário que construiu; ele o sacode e, na maioria das vezes, o derrama. A instabilidade que Newton imaginou é real. Ele só se enganou quanto ao autor.

Fica a imagem de um futuro distante em que uma anã branca, pequena como a Terra e cada vez mais fria, gira cercada por um séquito reordenado: talvez Júpiter e Saturno ainda ali, presos numa ressonância que nunca tiveram, talvez Urano ocupando o lugar de Netuno, talvez um anel de destroços onde antes havia planetas. E, em algum ponto escuro da galáxia, um gigante gelado seguindo sozinho, filho de uma estrela que já não o segura. Daqui a bilhões de anos, astrônomos de outros mundos poderão ler nos sistemas das anãs brancas o modo como cada estrela morreu. Entre esses fósseis dinâmicos estará, com suas órbitas embaralhadas pelos últimos suspiros do Sol, o que um dia chamamos de Sistema Solar.