Há uma cicatriz no céu do Hemisfério Norte que virou o objeto mais examinado da astronomia galáctica da última década. Ela se chama GD-1, tem quase cem graus de comprimento, é fina como um fio de cabelo esticado sobre a esfera celeste e carrega, em determinado ponto de sua extensão, uma falha na densidade de estrelas seguida de um pequeno apêndice que se projeta para fora da trilha principal. Desde 2018 essa falha e esse apêndice vêm sendo lidos como a impressão digital de algo que ninguém consegue ver: um grumo de matéria escura, com massa entre um milhão e cem milhões de vezes a do Sol, que teria passado raspando pela corrente estelar e arrancado um pedaço dela com a própria gravidade. Um trabalho publicado no The Astrophysical Journal em 27 de agosto de 2026 acaba de mostrar que a Via Láctea, sozinha, sabe fazer exatamente a mesma coisa.

O estudo é assinado por Arpit Arora, do Instituto DiRAC da Universidade de Washington, e por mais treze pesquisadores de instituições nos Estados Unidos, no Reino Unido, na Dinamarca e na Suíça. A estratégia foi de força bruta estatística: simular cerca de quinze mil correntes estelares originadas de aglomerados globulares, distribuídas em quatro galáxias do tamanho da Via Láctea extraídas de simulações cosmológicas, e deixá-las evoluir por cinco bilhões de anos dentro do campo gravitacional dessas galáxias. A parte decisiva do experimento está no que foi deliberadamente removido. Nenhum sub-halo de matéria escura foi incluído. Nenhuma nuvem molecular gigante. Nenhum perturbador de pequena escala de espécie alguma. Se as correntes ainda assim desenvolvessem lacunas, esporões e dobras, a culpa só poderia ser da própria galáxia hospedeira.

E elas desenvolveram. Cerca de três quartos das quinze mil correntes simuladas exibiram algum tipo de afastamento de uma trilha limpa e estreita. Mesmo entre as mais bem comportadas, aquelas que qualquer levantamento fotométrico classificaria como lisas e frias, a largura oscilou entre dez e vinte e cinco por cento ao longo do próprio corpo, e praticamente todas hospedaram sobredensidades e vazios em escalas angulares de cerca de dois graus — exatamente a faixa de um a cinco graus que a literatura atribui aos encontros com sub-halos escuros. De aproximadamente quinze mil correntes, apenas cerca de setenta chegaram ao fim da simulação sem nenhuma ondulação detectável em qualquer escala. O número é brutal na sua simplicidade: menos de meio por cento da população corresponde à corrente idealizada, fina, fria e uniforme, que serve de premissa para a maior parte das análises de subestrutura escura publicadas nos últimos vinte anos.

Para entender o peso desse resultado é preciso primeiro entender o que é uma corrente estelar e por que ela se tornou tão valiosa. Aglomerados globulares são bolas compactas de centenas de milhares de estrelas antigas, mantidas coesas pela própria gravidade, que orbitam o centro da Galáxia em trajetórias alongadas. A cada passagem pelo periélio galáctico, quando o aglomerado chega ao ponto mais próximo do centro, a maré gravitacional da Via Láctea puxa as estrelas das camadas externas com mais força do que puxa o núcleo do aglomerado. Algumas dessas estrelas escapam. Como escapam com velocidade orbital ligeiramente diferente da do progenitor, elas se espalham para frente e para trás ao longo da mesma órbita, formando dois braços — um que corre adiante do aglomerado e outro que fica para trás. Ao longo de bilhões de anos e dezenas de passagens, esses braços se esticam em um filamento tênue que pode dar meia volta no céu.

O valor científico dessa estrutura está na sua fragilidade. As estrelas de uma corrente têm dispersão de velocidades minúscula, da ordem de poucos quilômetros por segundo, o que significa que qualquer perturbação gravitacional externa deixa marca visível. É a razão pela qual as correntes viraram a ferramenta preferida para caçar aquilo que o modelo cosmológico padrão prevê mas nenhum telescópio consegue ver: os sub-halos de matéria escura de baixa massa. Abaixo de um certo limiar, esses grumos nunca acumularam gás suficiente para acender estrelas, então são completamente escuros. Não emitem, não absorvem, não refletem. A única maneira de detectá-los é pelo efeito gravitacional que exercem sobre outra coisa — e uma corrente estelar fria é, em tese, o detector perfeito.

O raciocínio padrão é elegante. Um sub-halo que passa perto de uma corrente aplica um empurrão gravitacional diferencial: as estrelas mais próximas do encontro recebem um impulso maior do que as mais distantes. O resultado, décadas depois, é um vazio na distribuição de estrelas ao longo da trilha, ladeado por acúmulos nas bordas, e frequentemente acompanhado de um apêndice de estrelas que foram jogadas para fora do plano da corrente. A largura angular desse vazio depende da massa do perturbador e da geometria do encontro, e os modelos preveem tipicamente algo entre um e cinco graus. Encontre lacunas dessa dimensão em uma corrente, conte quantas há, e você terá uma medida direta de quantos grumos escuros existem por unidade de volume na Galáxia — uma restrição observacional sobre a natureza da própria matéria escura.

Toda essa cadeia lógica repousa em uma única premissa, quase nunca testada de forma explícita: a de que uma corrente de aglomerado globular, deixada em paz, seria intrinsecamente lisa. A premissa faz sentido em um potencial gravitacional estático e simétrico, o tipo de modelo que domina a literatura teórica sobre correntes. Uma esfera de matéria escura suave, um disco axissimétrico, um bojo central — nesse cenário as estrelas escapadas se distribuem de maneira quase uniforme ao longo da trilha, com perfil transversal aproximadamente gaussiano, e qualquer irregularidade precisa vir de fora.

A Via Láctea real não é nada disso. Ela tem uma barra girando no centro, braços espirais que se formam e se desfazem, um disco massivo cuja assimetria evolui, um halo de matéria escura que é triaxial e inclinado em relação ao disco, e uma interação em curso com a Grande Nuvem de Magalhães que reorganiza o campo gravitacional de toda a periferia. Nada disso é estático em escalas de tempo comparáveis ao período orbital de uma corrente do halo, que gira em torno de algumas centenas de milhões de anos. Alguns desses efeitos já foram flagrados individualmente: a corrente associada a Pal 5 e a de Ophiuchus mostram variações de densidade e morfologia em leque atribuídas à barra galáctica; a corrente Orphan–Chenab exibe desalinhamentos de larga escala creditados à influência dependente do tempo da Grande Nuvem de Magalhães; a corrente ATLAS–Aliqa Uma tem uma dobra pronunciada na trilha que recebeu a mesma explicação.

O que faltava era juntar tudo. Os estudos anteriores costumavam acrescentar um componente por vez a um potencial idealizado — só a barra, só a Nuvem de Magalhães — ou recorrer a simulações de N-corpos com algum grau de simetria imposta. Ninguém havia perguntado, de forma sistemática e em escala populacional, quanta complexidade uma galáxia inteira, viva e evoluindo, imprime nas suas correntes sem nenhuma ajuda externa. Responder a isso exige duas coisas difíceis de conciliar: um modelo realista de galáxia, com toda a bagunça cosmológica de fusões e reorganizações, e um controle absoluto sobre o que está dentro e o que está fora do experimento.

A solução veio das simulações LATTE, parte do projeto FIRE-2, que reproduzem a formação de galáxias do tamanho da Via Láctea a partir de condições iniciais cosmológicas com física bariônica detalhada. Quatro halos foram escolhidos, cada um com uma história de montagem diferente. O halo m12i tem massa total de 1,18 trilhão de massas solares dentro de um raio de 336 quiloparsecs e é o mais isolado do conjunto, com apenas fusões menores nos últimos cinco bilhões de anos. O m12m, com 1,58 trilhão de massas solares e raio de 371 quiloparsecs, hospeda o disco mais antigo e mais massivo, com 112 bilhões de massas solares em estrelas. O m12f, o mais pesado com 1,71 trilhão de massas solares, sofreu duas fusões recentes com galáxias de massa comparável à de Sagitário. E o m12b, com 1,43 trilhão de massas solares, é o melhor análogo do sistema Via Láctea–Grande Nuvem de Magalhães, com um satélite de razão de massa de aproximadamente um para oito e periélio em torno de quarenta quiloparsecs.

A resolução dessas simulações define o que elas podem e não podem fazer. Cada partícula bariônica tem 7.100 massas solares e cada partícula de matéria escura, 35.000. Os comprimentos de suavização gravitacional são de aproximadamente quarenta parsecs para a matéria escura e quatro parsecs para as estrelas, com suavização adaptativa para o gás chegando ao mínimo de um parsec. Um aglomerado globular típico tem raio de escala entre quatro e dez parsecs e massa de algumas centenas de milhares de massas solares — ou seja, cabe em umas poucas partículas da simulação. Correntes de aglomerados globulares simplesmente não são resolvidas nesse tipo de cálculo. A saída foi injetá-las como populações de partículas-teste sem massa, que sentem a gravidade da galáxia mas não a produzem.

Para fazer isso com precisão, a equipe precisou converter as simulações em um potencial gravitacional analítico que evolui no tempo. O método usado, batizado de potencial multipolar de baixa ordem dependente do tempo, decompõe a distribuição de massa de cada instantâneo da simulação em uma soma de funções de base. O halo de matéria escura e o gás quente são representados por uma expansão em harmônicos esféricos sobre uma grade radial logarítmica; o bojo, o disco achatado e o gás frio, por uma expansão de Fourier sobre uma grade meridional bidimensional. Os coeficientes são ajustados a todas as partículas dentro de 600 quiloparsecs do centro e depois interpolados linearmente entre instantâneos consecutivos, salvos a cada 25 milhões de anos, produzindo uma força contínua no tempo.

O truque que torna o experimento possível está no corte da expansão. Ao restringir as ordens máximas a quatro tanto na parte esférica quanto na parte de Fourier, o modelo retém as estruturas de grande escala — a barra, os braços espirais, as deformações do disco, a resposta de baixa ordem do halo à passagem de satélites massivos — e descarta, por construção, tudo o que for de pequena escala. Sub-halos de matéria escura desaparecem. Nuvens moleculares gigantes desaparecem. O ruído de alta frequência da distribuição discreta de partículas é suavizado. É o equivalente gravitacional de tocar uma orquestra inteira usando apenas os instrumentos graves: a harmonia geral permanece, os detalhes agudos somem. O resultado é um campo de forças que contém toda a complexidade macroscópica de uma galáxia real e nenhum perturbador microscópico — a linha de base controlada que faltava.

As correntes foram então semeadas nesse campo. Os progenitores tiveram suas posições e velocidades atuais sorteadas de uma distribuição quase esférica e isotrópica baseada no perfil de traçadores de Dehnen, com raio de escala de seis quiloparsecs, que reproduz aproximadamente a distribuição do halo estelar. Cada órbita foi integrada cinco bilhões de anos para trás no potencial variável, e foram mantidos apenas os progenitores cujo periélio mínimo caiu entre dez e trinta quiloparsecs e cujo afélio máximo ficou entre vinte e cento e dez quiloparsecs. Sobraram cerca de 4.200 condições iniciais por halo. Nesse intervalo, uma corrente do halo interno completa aproximadamente dez órbitas, enquanto uma do halo externo completa apenas quatro, acumulando proporcionalmente menos passagens pelo periélio.

Cada progenitor foi modelado como um aglomerado com perfil de King, profundidade de potencial central fixa e raio de escala de quatro parsecs, com massas sorteadas de uma distribuição lognormal cobrindo de dez mil a dez milhões de massas solares, compatível com os sistemas de aglomerados globulares observados hoje. Um detalhe dessa escolha tem consequências que voltarão adiante: para as órbitas com periélio acima de vinte quiloparsecs, o raio de Jacobi — a fronteira além da qual a maré da Galáxia vence a gravidade do aglomerado — fica entre quarenta e quatrocentos e oitenta parsecs, com mediana de cerca de setenta e oito. Isso é dez a cem vezes maior que o raio do próprio aglomerado, o que significa que esses progenitores perdem pouquíssimo material a cada passagem e podem, na prática, nunca chegar a formar uma corrente observável.

A formação das correntes seguiu o método de aspersão de partículas, uma técnica que dispensa a simulação da dinâmica interna do aglomerado e libera as estrelas diretamente nos pontos de Lagrange do progenitor, onde a maré vence a coesão. Dez mil partículas-teste foram liberadas por corrente, duas a cada passo de integração, uma para o braço dianteiro e outra para o traseiro, com velocidades e deslocamentos de ejeção sorteados de distribuições empíricas calibradas contra simulações completas de N-corpos. A escolha de ejetar partículas uniformemente no tempo, e não concentradas nas passagens pelo periélio como acontece na realidade, é uma simplificação deliberada — e conservadora. A perda de massa real é episódica e concentrada, o que tenderia a produzir mais estrutura de densidade, não menos. Tudo o que o trabalho encontra é, portanto, um piso.

O primeiro contato com o resultado é visual e desconcertante. A Figura 1 mostra 250 correntes sorteadas ao acaso no halo m12i, projetadas em coordenadas galactocêntricas sobre um mapa fotométrico sintético da própria galáxia. Não há duas iguais. Algumas são traços finos e retos, como riscos de lápis; outras são manchas difusas, com bifurcações visíveis, envelopes largos em torno de um núcleo denso, interrupções abruptas. Nenhuma dessas estruturas foi produzida por um sub-halo, porque não existe sub-halo nenhum na simulação. Nos halos que passaram por fusões recentes, m12f e m12b, a tendência à bagunça é ainda mais acentuada.

Transformar essa impressão visual em número exigiu construir um conjunto de métricas do zero. Cada corrente foi primeiro colocada em um sistema de coordenadas próprio, ajustando um grande círculo à trilha e girando o sistema até minimizar a dispersão transversal — a projeção em que a corrente parece mais fina e mais coerente possível. Uma estimativa de densidade por núcleo ao longo da direção longitudinal define a região contígua que contém noventa por cento da massa, e é essa máscara que delimita a trilha principal. Dentro dela se medem o comprimento, tomado como a extensão em graus, e a largura, tomada como o intervalo que contém noventa por cento das estrelas na direção transversal. O exemplo mostrado no artigo tem 18,5 graus de comprimento e 0,30 grau de largura.

A escolha de trabalhar na projeção mais favorável é o que torna todo o resultado tão contundente. Um observador real não escolhe o ângulo de visada: vê a corrente do ponto onde o Sol está, e qualquer estrutura que exista na direção da linha de visada aparece projetada sobre a trilha, aumentando a aparente desordem. Ao medir tudo no referencial em que a corrente parece mais limpa, os autores garantem que nenhuma estrutura espúria foi introduzida pela geometria — e que a complexidade que encontram é um limite inferior sobre o que seria efetivamente observado.

As distribuições cumulativas de comprimento e largura já revelam diferenças entre os hospedeiros. Os comprimentos são parecidos nos quatro halos, refletindo a faixa comum de números de passagens pelo periélio, com a exceção do m12f, onde cerca de trinta por cento das correntes têm menos de dez graus, contra quinze a vinte por cento nos demais — assinatura das fusões recentes, que rompem correntes antes que elas se estiquem. As larguras variam mais. O halo isolado m12i produz correntes predominantemente finas, com setenta e cinco por cento abaixo de um grau e meio; no m12f a fração cai para a metade, e o m12m, dono do disco mais massivo, também engrossa suas correntes.

A métrica central do trabalho testa uma expectativa simples: se uma corrente é lisa, a distribuição de estrelas na direção transversal deve ser aproximadamente gaussiana em qualquer trecho curto da trilha. Um esporão, uma dobra, um envelope difuso — qualquer estrutura fora da trilha quebra essa gaussianidade, seja alargando o perfil, seja criando picos secundários, seja produzindo caudas assimétricas. Cada corrente foi dividida em oito segmentos de igual comprimento, e em cada segmento a distribuição transversal foi comparada com uma normal padrão usando a distância de Wasserstein, também chamada de distância do carregador de terra por medir literalmente quanto trabalho seria preciso para transportar uma distribuição de probabilidade até a forma da outra. O valor bruto foi normalizado por um limiar de significância de cinco por cento calibrado por Monte Carlo, de modo que uma pontuação abaixo de um significa compatibilidade com uma trilha gaussiana e valores maiores medem o quanto o segmento se afasta dela.

Dessa pontuação por segmento saem quatro diagnósticos. A perturbação global é a mediana ao longo dos oito segmentos e mede o nível típico de estrutura fora da trilha. A perturbação de pico é o valor máximo e captura o desvio localizado mais forte em qualquer ponto — o esporão isolado, exatamente o tipo de feição atribuída a um sub-halo. A fração gaussiana conta quantos segmentos passam no teste. E a variação de largura mede o desvio padrão relativo da espessura de segmento para segmento, sensível a qualquer perturbação que faça a corrente estreitar e alargar sem deslocar estrelas para fora da linha central.

Para a estrutura ao longo da trilha, e não transversal a ela, a receita foi outra. A densidade linear de estrelas foi estimada com um núcleo de largura de banda de 0,25 grau, comparável à resolução angular dos filtros usados nas buscas por correntes em levantamentos fotométricos de campo largo, e teve as tendências de grande escala removidas por um polinômio cúbico. O espectro de potência unidimensional do resíduo fracionário foi calculado pelo método de Welch e comparado com uma referência construída a partir de duas mil realizações de Monte Carlo que amostram uniformemente a mesma janela com o mesmo número de partículas. Onde a potência medida supera a flutuação puramente estatística por cinco desvios padrão, há estrutura real. A menor escala angular em que isso acontece é o diagnóstico que carrega o resultado mais perturbador do artigo.

Aplicadas à população inteira, restrita às correntes com pelo menos dez graus de comprimento — o que deixa 14.787 objetos, dos quais 3.909 no m12i, 4.093 no m12m, 3.493 no m12b e 3.292 no m12f —, essas métricas produzem uma classificação em quatro camadas. Correntes lisas, sem nenhuma estrutura significativa fora da trilha em nenhum segmento, somam vinte e seis por cento. Correntes lisas com uma feição isolada, trilha limpa mais um esporão ou dobra único, somam treze por cento. Correntes intermediárias, com trechos parcialmente desorganizados, somam vinte e sete por cento. E correntes bagunçadas, globalmente perturbadas, somam trinta e quatro por cento — a categoria mais numerosa de todas.

O que governa em qual dessas caixas uma corrente cai é, acima de tudo, a distância do periélio. Correntes cujos progenitores chegam mais perto do centro galáctico atravessam a região onde a barra, os braços espirais e o disco contribuem simultaneamente para o campo, e completam mais órbitas no mesmo intervalo, acumulando perturbação a cada passagem. A tendência é limpa e universal nos quatro halos: quanto maior o periélio, menor a perturbação global e maior a fração de correntes lisas. O limiar em torno de quinze quiloparsecs separa as duas populações de forma robusta, com variação entre hospedeiros — o isolado m12i chega ao regime liso já aos quinze quiloparsecs e tem cerca de oitenta por cento de correntes lisas além de vinte, enquanto o m12m só atinge o mesmo patamar perto de vinte quiloparsecs e os dois hospedeiros com fusões recentes decaem mais devagar e com dispersão maior.

A dependência da excentricidade é mais sutil e, curiosamente, troca de sinal com o raio. No halo interno, abaixo de quinze quiloparsecs, as órbitas quase circulares produzem correntes sistematicamente mais bagunçadas que as radiais — porque ficam confinadas à região complicada e nunca escapam dela, acumulando perturbação continuamente, enquanto uma órbita radial com o mesmo periélio passa boa parte do período longe, no afélio, diluindo o efeito integrado. Além de quinze quiloparsecs o padrão se inverte: a órbita circular permanece a uma distância aproximadamente constante em uma região tranquila, enquanto a excêntrica volta repetidamente para dentro e recolhe um novo empurrão a cada retorno. A escala de transição coincide com a extensão do disco massivo do m12m, o que sugere o disco como principal responsável, embora os autores se recusem a atribuir o limiar a um único componente sem os experimentos controlados que ainda estão por vir.

Um detalhe dentro da camada lisa merece atenção especial. Entre as correntes classificadas como globalmente lisas, aproximadamente metade contém pelo menos um segmento com estrutura localizada detectável fora da trilha. No m12i, cerca de setenta e cinco por cento das correntes lisas ficam abaixo do valor de referência a partir do qual essas feições se tornam visualmente identificáveis; no m12f a fração cai para cinquenta e cinco por cento. Traduzindo: uma corrente que parece impecável em um levantamento fotométrico, fina e coerente do começo ao fim, tem chance considerável de exibir um esporão isolado — a assinatura clássica do encontro com um sub-halo escuro — sem que nenhum sub-halo tenha existido.

A variação de largura conta uma história paralela. Entre as correntes lisas, a mediana cai de vinte e seis por cento nas órbitas de periélio entre dez e quinze quiloparsecs para vinte por cento naquelas entre vinte e cinco e trinta, e as distribuições se estreitam para fora. A interpretação é de choque de maré residual: cada passagem próxima aplica um forçamento impulsivo que deixa marca no perfil de espessura mesmo quando a distribuição transversal continua bem comportada. O ponto que interessa ao observador é outro, e está na distribuição cumulativa: nenhuma corrente lisa da amostra tem variação de largura abaixo de dez por cento, e a mediana fica em torno de vinte e dois por cento. Não existe corrente perfeitamente uniforme em espessura. Nenhuma. Modulações de largura desse tamanho já são detectáveis em levantamentos fotométricos profundos com contagem estelar suficiente, e vinham sendo interpretadas como evidência de impactos localizados.

O golpe mais duro vem da estrutura de densidade ao longo da trilha, medida apenas na camada mais limpa da população. Entre as correntes lisas com detecção significativa, a flutuação fracionária de densidade tem mediana de 0,17 — cerca de dezessete por cento de variação entre picos e vales ao longo do corpo da corrente — contra 0,07 no subconjunto sem detecção, uma diferença de fator dois e meio. A potência integrada acima do ruído estatístico separa ainda mais as duas populações, por um fator de aproximadamente sete. A dependência com os parâmetros orbitais existe, mas é secundária: as correntes de periélio menor têm amplitudes ligeiramente maiores, e além de vinte quiloparsecs os valores se estabilizam em torno de 0,16 sem queda adicional. Estrutura de densidade real aparece em todas as faixas orbitais.

E aparece na escala errada. A menor escala angular em que essa estrutura ultrapassa o ruído estatístico tem mediana de aproximadamente dois graus entre as correntes lisas, com a distribuição cobrindo de meio grau a cinco graus. É precisamente a janela onde os modelos preveem as lacunas produzidas por sobrevoos de sub-halos de matéria escura, que dependendo da massa e da geometria do encontro deixam marcas de um a cinco graus. Apenas cerca de dois por cento das correntes lisas, algo como setenta de três mil e oitocentas, não mostram nenhuma feição detectável em nenhuma escala. A conclusão que os autores extraem é explícita e não deixa margem: com base apenas na morfologia fotométrica, feições de densidade em escala de graus dentro de uma corrente lisa são indistinguíveis de perturbações induzidas por sub-halos.

Aqui entra a questão prática de quanto disso é efetivamente visível. Uma flutuação de dezessete por cento só se destaca se houver estrelas suficientes por elemento de resolução para que o ruído de Poisson fique abaixo dela. Com trinta estrelas por segmento, a flutuação puramente estatística já responde por cerca de dezoito por cento — o que coloca essas correntes exatamente no limiar de detecção dos levantamentos de campo largo atuais, cuja profundidade de referência gira em torno de magnitude 23,9. O Observatório Vera C. Rubin muda esse quadro: sua profundidade prevista, em torno de magnitude 26,75, é quase três magnitudes maior, o que multiplicará a contagem de estrelas por elemento de resolução para as correntes próximas e empurrará o piso de Poisson para baixo de dez por cento. Traduzido: a estrutura de densidade que este trabalho identifica na maioria das correntes simuladas vai se tornar fotometricamente recuperável na próxima geração de catálogos.

Sobram as setenta correntes. Cinco delas, escolhidas ao acaso e ordenadas por comprimento, aparecem no artigo como o análogo mais próximo do objeto idealizado que povoa a literatura de perturbação por sub-halos: comprimentos entre dez e trinta graus, larguras entre 0,3 e 0,5 grau, nenhuma estrutura de densidade detectável em escala alguma. Mesmo essas carregam variação de largura entre quinze e trinta e cinco por cento. São finas, são coerentes, são livres de lacunas — e ainda assim não são uniformes. O piso de complexidade não tem fundo perfeitamente liso em lugar nenhum da população.

A comparação com correntes reais fecha o argumento. A equipe vasculhou o catálogo simulado à procura de estruturas parecidas com as feições mais citadas do céu real e encontrou candidatas em ambos os casos. Uma corrente formada no halo isolado m12i, com periélio de cerca de dez quiloparsecs e afélio de vinte e quatro, exibe um esporão que, sobreposto às estrelas observadas de GD-1, coincide em escala angular e em forma. GD-1 tem periélio de quatorze quiloparsecs e afélio de vinte e dois — mesma vizinhança orbital. Outra corrente do mesmo halo, com periélio de quatorze quiloparsecs e afélio de vinte e cinco, reproduz qualitativamente a dobra da corrente ATLAS–Aliqa Uma, feição que vinha sendo creditada à influência da Grande Nuvem de Magalhães. Nenhuma das duas simuladas encostou em um sub-halo na vida.

Os casos não param aí. A estrutura bifurcada observada na corrente Jhelum tem paralelo nos perfis transversais bimodais da camada intermediária. Os desalinhamentos de larga escala da Orphan–Chenab encontram análogos naturais nos hospedeiros com fusões recentes, onde a passagem de um satélite massivo deixa marca em uma ampla faixa de raios orbitais. Envelopes difusos parecidos com os que cercam Jhelum e GD-1 aparecem por toda a camada bagunçada. O detalhe que mais incomoda é que nenhum dos quatro halos FIRE-2 desenvolveu uma barra forte e duradoura, e ainda assim as feições encontradas já se sobrepõem às estruturas ressonantes previstas para interações entre correntes e barra. Uma Via Láctea com a barra que ela de fato tem produziria um piso de complexidade ainda mais alto no halo interno.

Para verificar se toda essa desordem vem mesmo da galáxia e não de algum artefato do método de geração das correntes, os autores rodaram mil correntes adicionais em um potencial estático e axissimétrico da Via Láctea, com os mesmos cortes orbitais e o mesmo algoritmo de aspersão. O contraste é categórico em todos os eixos de diagnóstico. A perturbação global fica abaixo de um em todos os periélios. A variação de largura fica abaixo de vinte por cento para setenta e cinco por cento das correntes. A flutuação de densidade tem mediana de 0,07, menos da metade do valor da amostra cosmológica. E, o número decisivo, cerca de trinta e cinco por cento das correntes do potencial estático não mostram nenhuma estrutura detectável ao longo da trilha, contra apenas dois por cento na amostra FIRE-2. Somente quinze por cento delas apresentam estrutura em escalas de cerca de dois graus, contra sessenta por cento das correntes lisas cosmológicas.

A diferença entre esses dois experimentos é exatamente a diferença entre uma galáxia de livro-texto e uma galáxia de verdade. Evidências preliminares apontam para o disco como o principal culpado no halo interno: bifurcações observadas no m12i coincidem com o surgimento de estrutura não axissimétrica forte no disco, e desaparecem quando as mesmas correntes são reintegradas com o disco restrito à sua componente axissimétrica ou congelado na forma não axissimétrica que tinha — o que indica que tanto a assimetria quanto sua evolução temporal são necessárias. A decomposição sistemática dessa contribuição é objeto de um trabalho em preparação. Para o halo externo, onde a influência do disco enfraquece, os candidatos são a forma assimétrica e mutante do próprio halo de matéria escura e a natureza caótica das órbitas em um potencial triaxial. As interações com satélites massivos são o terceiro contribuinte, e a assinatura delas está clara na comparação entre os hospedeiros com fusões recentes e o isolado.

O rigor numérico do trabalho merece nota, porque um resultado desse tipo é exatamente o que se esperaria de um erro de simulação. Três apêndices se dedicam a descartar essa possibilidade. O ruído estatístico introduzido pelo ajuste das funções de base à distribuição finita de partículas foi calibrado por reamostragem, ficando entre cinco e sete décimos de milésimo da força total, plenamente coerente ao longo da seção transversal de uma corrente típica — ou seja, incapaz de produzir empurrões diferenciais. A integração das órbitas foi testada contra um integrador simplético de passo fixo, com resíduos de posição de mediana de dois milésimos de quiloparsec, vinte vezes menores que a espessura típica de uma corrente, e viés de velocidade acumulado ao longo de cinco bilhões de anos abaixo de um milésimo da velocidade orbital. E a fidelidade do potencial de base foi verificada contra o cálculo direto da aceleração a partir das partículas da simulação por um código de árvore, com concordância melhor que dois por cento em todos os raios acima de sete quiloparsecs.

A saída que os autores enxergam para o impasse não está na fotometria, e sim na cinemática. Movimentos próprios e velocidades radiais ao longo da trilha são sensíveis à coerência de velocidades e à estrutura no espaço de fases da perturbação, e há razão física para esperar que um encontro impulsivo com um sub-halo deixe uma assinatura cinemática diferente da produzida pelo campo de maré secular do hospedeiro. Testar isso é o próximo passo declarado: introduzir populações explícitas de sub-halos nos mesmos potenciais, gerar catálogos de correntes perturbadas e verificar se a informação completa de seis dimensões no espaço de fases consegue separar as duas origens. A linha de base estabelecida agora é a fundação sobre a qual essa separação vai ser tentada.

Enquanto isso, o que fica é um recado desconfortável para uma linha de pesquisa inteira. Restrições sobre a função de massa dos sub-halos de matéria escura vêm sendo derivadas de feições morfológicas em correntes individuais, com GD-1 no papel de estrela principal, sob a premissa de que uma corrente limpa não produz lacunas sozinha. Este trabalho mostra que produz, e que produz nas mesmas escalas angulares, com amplitudes comparáveis, em uma fração majoritária de uma população de quinze mil objetos. As feições não são raras nem exigem ajuste fino de parâmetros: emergem genericamente do campo de maré secular de qualquer galáxia com uma história de montagem realista. Os próprios autores registram que a complexidade combinada, vinda do potencial dependente do tempo e das passagens de sub-halos, provavelmente impede uma interpretação limpa das feições encontradas em qualquer corrente isolada. A separação, se for possível, terá caráter estatístico, obtida sobre grandes populações que cubram órbitas e ambientes diversos.

Isso não invalida as correntes estelares como sondas de matéria escura. Muda o que precisa ser feito antes de usá-las. Uma análise que compare a morfologia observada de GD-1 com a expectativa de um potencial estático está comparando com a coisa errada; a referência correta é uma população de correntes evoluída em um hospedeiro vivo, com barra, disco assimétrico, halo triaxial e satélites em queda. É o que este catálogo de quinze mil objetos oferece pela primeira vez. Convém notar que os números específicos — vinte e seis por cento de correntes lisas, trinta e quatro por cento de bagunçadas — estão condicionados à distribuição orbital amostrada e aos quatro potenciais escolhidos, e não devem ser lidos como previsão direta para a população real da Via Láctea. As tendências por parâmetro orbital, essas sim, são o produto transferível do trabalho.

O cronograma astronômico dá urgência ao assunto. O Legacy Survey of Space and Time, do Observatório Rubin, vai cobrir cerca de vinte mil graus quadrados com profundidade coadicionada em torno de magnitude 27,5. A missão Euclid, da Agência Espacial Europeia, imageará aproximadamente quinze mil graus quadrados no óptico e no infravermelho próximo. O Telescópio Espacial Nancy Grace Roman entregará fotometria profunda de alta resolução em cerca de três mil graus quadrados. Juntos, esses instrumentos vão resolver flutuações de densidade, variações de largura e feições fora da trilha em escalas de fração de grau para dezenas a centenas de correntes ao mesmo tempo. O censo de correntes conhecidas na Via Láctea, que saiu de meia dúzia para quase cem na última década, deve crescer de novo por um fator considerável.

Cada uma dessas correntes vai chegar cheia de detalhes morfológicos, e cada detalhe vai exigir uma decisão sobre sua origem. Sem uma linha de base como esta, a tentação de ler cada lacuna como um grumo escuro será grande e o erro sistemático resultante, silencioso. O trabalho de Arora e colaboradores não apaga a matéria escura do mapa das correntes estelares; apenas mostra que o mapa vinha sendo desenhado sem a escala. Ainda pode haver sub-halos escuros esculpindo GD-1. Só que agora será preciso provar que a Via Láctea, sozinha, com sua barra, seu disco torto e sua companheira massiva em queda livre, não teria feito o serviço primeiro.

Existe uma ironia elegante em tudo isso, e ela vale como fecho. Durante vinte anos, as correntes estelares foram tratadas como um instrumento — um sensor passivo, calibrado, cuja única função era registrar o que passasse perto delas. O que este catálogo de quinze mil filamentos mostra é que elas nunca foram instrumento nenhum. São parte da galáxia, feitas do mesmo material, sujeitas à mesma história, deformadas pelos mesmos processos que moldaram o disco e o halo que as abrigam. Ler uma corrente estelar é ler a biografia da Via Láctea escrita em uma tipografia que ninguém tinha decifrado até agora. As cicatrizes que os astrônomos vinham interpretando como marcas de agressores invisíveis podem ser, em boa parte, apenas rugas de uma galáxia que envelheceu.