Num mapa de velocidades produzido por um espectrógrafo moderno, uma galáxia deixa de ser uma mancha difusa de luz e vira um relógio em movimento: metade do disco aparece pintada de azul, porque aquele lado vem na nossa direção, e a outra metade aparece em vermelho, porque se afasta. Entre as duas cores passa uma linha invisível, o eixo em torno do qual bilhões de estrelas e nuvens de gás executam uma volta que dura centenas de milhões de anos. A pergunta que atravessa a cosmologia há mais de meio século é de onde veio esse giro. Uma equipe internacional de astrônomos acaba de responder com um número que muda o patamar do debate: a direção em que cada galáxia gira hoje foi gravada antes de existirem estrelas, e essa gravação ainda pode ser lida no céu com uma confiança estatística de 7 sigma.

O resultado foi publicado na última quarta-feira, 5 de agosto, na revista Nature Astronomy, num artigo liderado por Ming-Jie Sheng e Hao-Ran Yu, da Universidade de Xiamen, na China, com colaboradores de instituições chinesas, de Taiwan e do Canadá, incluindo o Instituto Perimeter e a Universidade de Toronto. O grupo comparou o giro medido de centenas de galáxias próximas com o giro previsto a partir do campo de densidade primordial do nosso pedaço do Universo, reconstruído por computador. Entre todas as populações testadas, o gás das galáxias elípticas massivas centrais entregou o sinal mais límpido: uma correlação de direção que se afasta do acaso no nível de 7 sigma, a evidência observacional mais forte já obtida para a chamada teoria do torque de maré.

Momento angular não é um detalhe técnico na vida de uma galáxia. Ele governa a forma, o tamanho e a cinemática interna desses sistemas, e ajuda a explicar por que o céu abriga desde discos espirais finos como lâminas até esferoides avermelhados e quase parados. Para a cosmologia, a rotação carrega um atrativo extra: por ser uma grandeza vetorial, com direção e sentido, ela guarda informações que os mapas de densidade, feitos apenas de números escalares, não conseguem registrar. Nessa camada adicional de informação podem estar escondidas respostas sobre a massa dos neutrinos, sobre possíveis violações de paridade nas leis físicas, sobre a natureza da energia escura e sobre alinhamentos de larga escala que os levantamentos do céu começam agora a medir.

Existe, porém, um abismo entre a teoria e o telescópio. O roteiro padrão da formação de estruturas diz que as galáxias nascem e evoluem dentro de halos de matéria escura, concentrações invisíveis que respondem pela maior parte da massa de cada sistema. São os halos que adquirem rotação primeiro, e as galáxias herdam parte desse movimento ao se formarem em seu interior. Só que nenhum instrumento observa diretamente o giro de um halo de matéria escura. O que os astrônomos conseguem medir é o giro do material luminoso, estrelas e gás, e precisam usá-lo como procurador de algo que jamais verão. Toda a arquitetura do novo estudo foi montada para atravessar esse abismo com segurança.

A teoria posta à prova tem endereço e data de nascimento. Em 1969, o físico James Peebles, que décadas depois receberia o Nobel, publicou o cálculo fundador sobre a origem do momento angular das galáxias; no ano seguinte, o soviético Andrei Doroshkevich refinou a ideia, e Simon White a consolidou em 1984. O mecanismo proposto ficou conhecido como teoria do torque de maré e parte de uma observação simples: as sementes das galáxias, os chamados proto-halos, não nasceram esféricas. Eram borrões alongados e irregulares de matéria, imersos num campo gravitacional que também não era uniforme, porque a distribuição de matéria ao redor puxava cada região com forças de intensidades e direções diferentes, o mesmo tipo de efeito diferencial que a Lua exerce sobre os oceanos da Terra.

Imagine uma tábua comprida boiando atravessada num rio. Se a correnteza fosse perfeitamente paralela à tábua, ela apenas deslizaria. Como a água empurra com mais força uma ponta do que a outra, e num ângulo que não coincide com o eixo da madeira, a tábua começa a rodar. Foi exatamente isso que aconteceu com os proto-halos no Universo jovem: o desalinhamento entre a orientação do borrão de matéria, descrita pelos físicos por um objeto matemático chamado tensor de inércia, e a orientação do campo de marés gravitacional ao redor, descrita pelo tensor de maré, produziu um torque. Esse torque imprimiu uma rotação coerente, que cresceu enquanto o proto-halo ainda se expandia junto com o cosmos e se congelou quando ele colapsou e se virializou, tornando-se um halo propriamente dito.

A previsão central da teoria é que essa herança não se perde por completo. Simulações de computador mostram que, mesmo no Universo atual, a rotação dos halos de matéria escura continua correlacionada com o campo de marés primordial, tanto em direção quanto em magnitude. Mostram também que as galáxias centrais tendem a girar alinhadas com seus halos hospedeiros, e que ambos preservam memória do ambiente gravitacional inicial. Se tudo isso estiver certo, o eixo de rotação que medimos hoje numa galáxia qualquer funciona como uma agulha de bússola apontando para condições físicas de mais de 13 bilhões de anos atrás. Faltava demonstrar, com dados reais e rigor estatístico, que a agulha de fato aponta.

Tentativas anteriores chegaram perto. Os giros tridimensionais de galáxias espirais podem ser inferidos a partir de suas elipticidades, ângulos de projeção, do sentido de enrolamento dos braços, de assimetrias Doppler e da absorção por poeira, e a espectroscopia de campo integral, tradução do termo em inglês integral field unit, técnica que cobre a galáxia inteira com uma malha de fibras ópticas e obtém um espectro em cada ponto, permitiu mapear diretamente os campos de velocidade. Com esse arsenal, um estudo de 2021 liderado por Pavel Motloch, e que já contava com Hao-Ran Yu e Ue-Li Pen entre os autores, encontrou uma correlação entre giros de espirais e as condições iniciais reconstruídas, mas no nível de cerca de 3 sigma, um indício sugestivo que ainda não configurava detecção robusta.

Para dar o salto, o novo trabalho apoiou-se em dois pilares de dados. O primeiro é o levantamento MaNGA, sigla de Mapping Nearby Galaxies at Apache Point Observatory, parte da quarta geração do Sloan Digital Sky Survey. O MaNGA acoplou feixes hexagonais de fibras ópticas ao telescópio de Apache Point, no Novo México, e produziu mapas bidimensionais da velocidade ao longo da linha de visada para cerca de dez mil galáxias próximas, cobrindo tanto o componente estelar quanto o gás ionizado. A equipe trabalhou com a liberação final desses dados, a mesma que integra o décimo sétimo lançamento público do Sloan, e extraiu dali a matéria-prima observacional do estudo: o giro projetado de cada galáxia no plano do céu.

Medir esse giro exige cuidado artesanal. Os pesquisadores ajustaram, com o pacote computacional PAFIT, o ângulo de posição cinemático dos campos de velocidade, isto é, a orientação da linha que separa o lado que se aproxima do lado que se afasta, tanto para o gás, rastreado pela linha de emissão Hα do hidrogênio, quanto para as estrelas. O vetor de rotação foi definido perpendicular a essa linha, com o sentido determinado pela regra física de que a matéria circula do lado avermelhado para o lado azulado. Entraram na análise apenas medidas confiáveis: pixels espectrais de gás com razão sinal-ruído mínima de 3, campos estelares com sinal mediano equivalente, exclusão de galáxias sem linhas de emissão suficientes e corte de todos os ajustes com incerteza angular acima de 20 graus.

O segundo pilar é uma proeza computacional chamada ELUCID, projeto conduzido ao longo da última década por Huiyuan Wang, Xiaohu Yang e colaboradores. A ideia soa quase abusada: pegar o filme da formação de estruturas do nosso canto do Universo e rodá-lo de trás para a frente, até o primeiro fotograma. O ponto de partida foi um catálogo de grupos de galáxias construído sobre os dados do Sloan, restrito à calota galáctica norte, a distâncias correspondentes a desvios para o vermelho entre 0,01 e 0,12, e a grupos com massas acima de um trilhão de massas solares. Distorções observacionais clássicas, como o efeito conhecido pelo apelido em inglês fingers of God, o alongamento aparente dos aglomerados na direção da linha de visada causado pelos movimentos internos das galáxias, foram suprimidas ou corrigidas.

Sobre esse mapa do presente, o ELUCID aplicou um algoritmo de Monte Carlo hamiltoniano acoplado a uma dinâmica de partícula-malha com vinte passos de tempo, buscando o campo de densidade inicial que, evoluído pela gravidade, reproduzisse as estruturas observadas hoje. O volume reconstruído é um cubo de 500 megaparsecs por h de lado, o que equivale a quase 700 megaparsecs, ou cerca de 2,3 bilhões de anos-luz, discretizado em 500 células por dimensão com resolução espacial de 1 megaparsec por h. Testes com simulações de N-corpos mostraram que as condições iniciais recuperadas regeneram as estruturas reais até escalas de 2 megaparsecs por h, o que corresponde a halos com cerca de 3 × 10¹² massas solares.

Com as condições iniciais em mãos, a equipe rodou sua própria simulação restrita usando o código CUBE2, desenvolvido pelo grupo de Yu para maximizar eficiência de memória. A cosmologia adotada seguiu os parâmetros do satélite WMAP em sua análise de cinco anos, com densidade de matéria de 0,258, densidade de energia escura de 0,742, constante de Hubble adimensional de 0,72, amplitude de flutuações σ8 de 0,796 e índice espectral primordial de 0,963. Foram 125 milhões de partículas, posicionadas inicialmente numa grade regular, deslocadas conforme a aproximação de Zel’dovich a partir do desvio para o vermelho 100 e evoluídas até o presente com um solucionador gravitacional do tipo partícula-partícula-malha.

Restava um obstáculo conceitual espinhoso. Na sua formulação original, a teoria do torque de maré calcula o momento angular inicial a partir do produto entre o tensor de inércia do proto-halo e o tensor de maré, e o tensor de inércia de um objeto que existiu há 13 bilhões de anos não pode ser medido a partir de princípios primeiros em observação nenhuma. A saída veio de um artigo de 2020 na Physical Review Letters, assinado por Yu, Pen e colegas, que propôs um substituto engenhoso batizado de reconstrução de spin, palavra inglesa que a comunidade usa como sinônimo compacto para a direção do momento angular. Em vez do tensor de inércia inacessível, o método usa o próprio campo de marés, suavizado em duas escalas diferentes.

A intuição por trás do truque é elegante. Quando se olha o campo gravitacional primordial em torno de um ponto com lentes de zoom distintas, os eixos principais do tensor de maré giram conforme a escala aumenta, porque camadas cada vez maiores de matéria entram na conta. A torção entre o campo visto na escala menor, que rastreia a anisotropia do próprio proto-halo, e o campo visto numa escala ligeiramente maior, que caracteriza o ambiente ao redor, funciona como aproximação do giro inicial. Simulações mostraram que a receita funciona melhor quando a escala menor vale cerca de 0,8 vez o raio lagrangiano do halo, a dimensão da região primordial que continha toda a sua massa, e quando a segunda escala tende à primeira por cima.

Havia ainda que ligar cada galáxia de hoje ao seu ponto de origem no Universo bebê. Os autores cruzaram as galáxias do MaNGA com os grupos do catálogo usado pelo ELUCID, identificando o halo hospedeiro de cada uma, e aplicaram uma técnica de remapeamento do espaço lagrangiano: na simulação restrita, a trajetória de cada partícula define um campo de deslocamentos, e a inversão desse campo, interpolada na posição atual de um objeto, devolve a posição que a matéria correspondente ocupava no início dos tempos. Nesses locais primordiais, e só neles, o giro previsto foi calculado com a fórmula das duas escalas, depois projetado no plano do céu de cada galáxia para permitir comparação justa com o dado observado.

A amostra final tem números precisos. Ficaram no estudo apenas galáxias centrais, as dominantes de seus grupos, habitando halos com massa mínima de um trilhão de sóis e dentro do volume reconstruído pelo ELUCID, com morfologias classificadas por um catálogo de aprendizado profundo do próprio Sloan. Isso resultou em 781 elípticas centrais e 815 espirais centrais com giro do gás medido com confiança, além de 635 elípticas e 812 espirais com giro estelar confiável, conjuntos que não se sobrepõem por completo porque os cortes de qualidade foram impostos separadamente em cada componente. Para cada objeto, dois vetores no plano do céu ficaram frente a frente: o giro observado e o giro previsto pelo campo primordial.

Um exemplo concreto ajuda a visualizar a corrente inteira. A galáxia catalogada como 1-436574 no MaNGA é uma elíptica situada a um desvio para o vermelho de 0,058, o que a coloca a cerca de 800 milhões de anos-luz, com massa estelar próxima de 2 × 10¹¹ massas solares, o triplo da Via Láctea em estrelas, dentro de um halo de 4 × 10¹³ massas solares. No campo primordial reconstruído, o contorno do seu proto-halo aparece desenhado sobre o potencial gravitacional suavizado, com uma seta indicando o giro previsto pela torção das marés. No mapa de velocidades do gás da mesma galáxia, obtido bilhões de anos depois, a seta medida aponta para uma direção notavelmente parecida.

Casos isolados, claro, não provam teoria nenhuma; o jogo se decide na estatística. Para cada galáxia, os autores calcularam o cosseno do ângulo entre o vetor observado e o vetor previsto, um número que vale 1 para alinhamento perfeito, menos 1 para giros opostos e se distribui em torno de zero quando não existe relação alguma entre as duas direções. A correlação de uma amostra inteira foi resumida pela média ponderada desses cossenos, e a escolha do peso revela o refinamento do método: nem todos os cantos do Universo primordial torciam suas sementes com a mesma eficiência, e tratar todos igualmente diluiria o sinal.

Entra em cena o parâmetro de ambiente de maré, escrito como βTT, uma quantidade adimensional extraída do próprio campo primordial que mede quanta torção as marés conseguiam gerar na escala de cada proto-halo. Ele se anula onde nenhum torque era possível e cresce com a força do desalinhamento entre as duas escalas de suavização. Simulações indicam que proto-halos com βTT alto têm giros previstos com mais confiabilidade, e por isso cada galáxia entrou na média com peso igual à raiz quadrada do seu βTT, escolha conservadora cuja variação, testada entre os expoentes 0,5 e 2, muda a amplitude do resultado sem alterar sua significância. As distribuições desse parâmetro são quase idênticas para elípticas e espirais, garantia de que a ponderação não favoreceu artificialmente um time.

Com a máquina montada, os pesquisadores variaram a escala de suavização do campo de marés de zero a 10 megaparsecs por h e mediram a correlação para quatro combinações de população e traçador: elípticas com gás, elípticas com estrelas, espirais com gás e espirais com estrelas. O veredito visual é imediato. A curva do gás das elípticas centrais dispara acima das demais e atinge o pico de 0,13 na escala ótima de 3,4 megaparsecs por h, cerca de 15 milhões de anos-luz. As espirais exibem correlações bem mais tímidas, e o componente estelar fica atrás do gasoso nas duas morfologias. Diante disso, a análise principal se concentrou no gás das elípticas centrais.

O passo seguinte transformou uma previsão qualitativa da teoria em teste quantitativo. Se o torque nasce na escala do proto-halo, e proto-halos mais massivos ocupavam regiões maiores, a escala ótima de suavização deve crescer com a massa do halo, acompanhando a relação teórica de 0,8 vez o raio lagrangiano. Ao dividir a amostra em faixas de massa com números parecidos de objetos, foi isso que os dados mostraram para todos os halos acima de cerca de 3 × 10¹² massas solares, com a busca da escala ótima restrita ao intervalo físico entre 0,5 e 1,1 raio lagrangiano para evitar picos espúrios. Abaixo desse limiar de massa abre-se a chamada região de decorrelação, onde os erros do remapeamento lagrangiano e as incertezas da reconstrução crescem e sufocam o sinal.

Nenhuma correlação vale muito sem saber o quanto o acaso conseguiria imitá-la. A resposta veio de 40.000 realizações de Monte Carlo construídas sob a hipótese nula de que giros observados e previstos nada têm a ver entre si: em cada realização, os giros medidos foram embaralhados aleatoriamente dentro de cada faixa de massa, preservando a distribuição intrínseca das direções, o perfil de massas e a seleção da amostra, com a mesma ponderação por βTT aplicada ao embaralhamento. Para cada faixa definiu-se uma significância padronizada, comparando o valor medido à média e ao desvio padrão da distribuição nula, e um índice global somou os desvios de todas as faixas. Como cada halo pertence a uma única faixa, as amostras são disjuntas, o que reduz a covariância entre elas, embora um resíduo fraco, vindo do uso do mesmo campo reconstruído, faça os autores tratarem o número final como medida aproximada.

Os números finais impõem respeito. No caso de referência, com a correlação avaliada na escala ótima de cada faixa de massa, o sinal se afastou da expectativa nula em 7,0 sigma, um desvio cuja probabilidade de surgir por flutuação aleatória é inferior a um caso em algumas centenas de bilhões. Dois testes deliberadamente mais duros confirmaram o quadro: fixando a escala no valor teórico de 0,8 raio lagrangiano, sem qualquer otimização, a significância ficou em 6,3 sigma; permitindo que cada realização embaralhada otimizasse sua própria escala, competindo em igualdade com os dados reais, ela ficou em 5,9 sigma. A faixa de massa mais baixa, dentro da região de decorrelação, foi excluída dessas contas, e os autores registram que incluí-la reduziria o total sem tocar na conclusão.

Traduzindo o estatisticês: as galáxias de hoje giram, de forma mensurável, na direção ditada pelas marés gravitacionais que atuavam sobre suas sementes antes da formação das primeiras estrelas. A hipótese rabiscada por Peebles em 1969 deixou de ser um indício de 3 sigma para virar uma detecção no padrão que a física de partículas exige para anunciar descobertas, algo raro na astronomia extragaláctica. O céu local, com suas centenas de elípticas aparentemente banais, comportou-se como um arquivo que preservou, por 13,8 bilhões de anos, a assinatura vetorial das condições iniciais do cosmos.

Resta explicar o aparente paradoxo de o sinal mais limpo morar justamente nas elípticas, e não nas espirais que dominaram as tentativas anteriores. A primeira razão é demográfica. A amostra elíptica habita halos sistematicamente mais massivos que a espiral, e massa ajuda duas vezes: halos grandes retêm conexão intrínseca mais forte com o campo primordial, segundo as simulações, e ocupam regiões lagrangianas maiores, exatamente onde a reconstrução do ELUCID e o remapeamento funcionam com mais precisão. Boa parte dos halos que hospedam espirais cai abaixo do limiar de 3 × 10¹² massas solares, dentro da zona onde o método perde o fôlego.

Um teste com massa controlada confirmou que a demografia não conta a história toda. Comparando gás de elípticas e de espirais em duas faixas idênticas de massa do halo, entre 6 × 10¹² e 10¹³ massas solares e entre 10¹³ e 1,8 × 10¹³, as espirais passaram a exibir correlações detectáveis, sinal de que sua fraqueza na amostra completa vinha em parte dos halos leves. Ainda assim as elípticas seguiram na frente: 0,25 contra 0,10 na faixa inferior, com significâncias de 3,8 contra 1,7 sigma, e 0,18 contra 0,08 na superior, com 3,4 contra 1,2 sigma. Observações recentes indicam que o alinhamento entre galáxias centrais e seus halos depende da morfologia, mais forte nas elípticas, o que ajudaria a explicar a diferença residual.

As outras duas razões apontam para o gás. Trabalhos anteriores do mesmo grupo, baseados em simulações hidrodinâmicas, já indicavam que o componente gasoso guarda mais memória do campo de marés primordial do que o estelar, cujo momento angular é remexido por fusões, barras e trocas internas ao longo da vida da galáxia. Existe também um motivo puramente instrumental: nas elípticas da amostra, os ajustes do ângulo de posição cinemático do gás carregam incertezas sistematicamente menores que os das estrelas, com a maioria abaixo de 2 graus, prova de que o campo de velocidades do gás ionizado é mais coerente e rende medidas de giro mais firmes.

Há uma ironia deliciosa nesse arranjo. Elípticas massivas são tratadas no jargão como galáxias vermelhas e mortas, cemitérios estelares que pararam de formar sóis há bilhões de anos. Pois é justamente o tênue gás ionizado que ainda circula nesses sistemas, às vezes assentado em discos finos e bem-comportados, o melhor fóssil dinâmico do Universo primordial identificado até agora. O material mais discreto das galáxias mais quietas carrega a mensagem mais antiga.

A escolha de trabalhar com halos, e não diretamente com as galáxias, também foi testada e se mostrou decisiva. Em princípio, seria possível remapear a posição de cada galáxia de volta às condições iniciais usando o mesmo campo de deslocamentos, caminho adotado no estudo de 2021. Os halos, contudo, correspondem a escalas físicas e lagrangianas maiores, onde o erro relativo do remapeamento diminui, e a simulação restrita do ELUCID, ancorada num catálogo de grupos, reproduz halos com mais fidelidade do que objetos individuais. Quando os autores tentaram o remapeamento direto das galáxias, a correlação despencou; a ponte pelos halos hospedeiros foi parte essencial do sinal recuperado.

Convém dosar o entusiasmo com uma régua honesta sobre o que o resultado não é. Um cosseno médio de 0,13 está longe de 1, e a cadeia de degradações é conhecida das simulações: a correlação entre o giro do proto-halo e o do halo atual fica em torno de 0,7, cai para a faixa de 0,5 a 0,6 quando o giro inicial é estimado pela teoria ou pelo método das duas escalas, e desce para 0,3 a 0,5 quando o traçador é o material bariônico das galáxias, com o gás um pouco acima das estrelas. Sobre isso ainda pesam erros de reconstrução do ELUCID, incertezas do remapeamento, ruído nos ângulos medidos, efeitos de projeção e atribuições equivocadas de halo. O sinal detectado é uma correlação residual de conjunto, não uma máquina de reconstruir o giro de galáxias individuais.

Outra limitação assumida com franqueza envolve a dimensionalidade. Nas espirais, o giro tridimensional pode ser recuperado combinando o ângulo cinemático com a inclinação estimada pela razão entre os eixos aparentes, sob a hipótese de disco fino e circular. Nas elípticas o caminho é traiçoeiro: mesmo quando o gás se acomoda num disco rotante, o corpo estelar costuma ser triaxial, com cinemática interna complexa que frustra estimativas confiáveis de orientação plena. Em nome da consistência, todos os giros do estudo foram tratados em duas dimensões, projeção que descarta um grau de liberdade direcional sem apagar a correlação intrínseca codificada nas direções completas, tarefa que observações futuras de maior resolução e modelagem dinâmica aprimorada poderão resgatar.

Por que tanto esforço para medir um cosseno? Porque o momento angular abre uma janela que os mapas de densidade não alcançam. Como grandeza de paridade ímpar, que distingue o Universo de sua imagem no espelho, o giro das galáxias é sensível a fenômenos que trocam de sinal sob reflexão: uma eventual quiralidade primordial, violações de paridade nas interações fundamentais, assinaturas de energia escura dinâmica. Neutrinos massivos, ao viajarem quase à velocidade da luz pelo cosmos jovem, alteram sutilmente a aquisição de rotação dos halos, de modo que estatísticas de giro em grandes amostras podem apertar os limites sobre a soma das massas dessas partículas fantasmagóricas, um dos números mais perseguidos da física atual.

Existe ainda um gargalo metodológico que o novo observável promete atacar. As técnicas de reconstrução das condições iniciais, o ELUCID entre elas, apoiam-se em posições e densidades de galáxias e estacionaram há mais de uma década na barreira das pequenas escalas, abaixo de uns 2 megaparsecs por h, justamente onde se esconde informação cosmológica rica. Incorporar os giros como canal independente de dados, algo que só faz sentido com amostras gigantes, é uma das rotas mais promissoras para furar esse bloqueio, e o presente trabalho funciona como prova de conceito de que a informação está lá, gravada e legível.

O futuro imediato tem nome e cronograma. O instrumento espectroscópico DESI, no Arizona, está observando milhões de galáxias brilhantes centrais de aglomerados, e os autores calculam que apenas duas a quatro fibras bem posicionadas por objeto bastam para determinar seu giro de forma robusta, mesmo sem a cobertura densa dos feixes do MaNGA. A meta desta fase é consolidar a detecção direcional em volumes cada vez maiores; estender o teste à magnitude do momento angular, e não só à sua orientação, vai exigir medidas mais precisas e uma teoria mais completa da evolução rotacional de cada componente galáctico. Reconstruções mais acuradas e amostras maiores de espectroscopia de campo integral devem levar a detecção das elípticas massivas para halos mais leves e populações mais amplas, espirais incluídas.

Quem quiser conferir cada vírgula pode fazê-lo hoje mesmo. Os giros observados e reconstruídos foram depositados no repositório Zenodo junto com um caderno de Python que reproduz os gráficos de diagnóstico e as medidas de correlação, os dados do MaNGA e os catálogos de morfologia e de grupos são públicos, e tanto o código de simulação CUBE2 quanto as rotinas de reconstrução de giro e remapeamento estão abertos no GitHub. Numa era de resultados extraordinários que às vezes se escondem atrás de dados fechados, a transparência aqui é total, à altura da ambição da alegação.

Meio século separa o lápis de Peebles e Doroshkevich, rabiscando torques sobre borrões de matéria num Universo que a cosmologia mal começava a mapear, deste 7 sigma obtido com dez mil espectros, um cubo de 2,3 bilhões de anos-luz reconstruído por Monte Carlo e 125 milhões de partículas simuladas. Entre um ponto e outro, a teoria do torque de maré atravessou gerações como hipótese elegante à espera de veredito. O veredito chegou, e com ele um programa de pesquisa inteiro: integrar o momento angular aos arcabouços de observação e reconstrução das próximas décadas, transformando uma propriedade antes tratada como ruído astrofísico em régua de precisão para pesar neutrinos, testar a paridade do cosmos e sondar a energia escura.

Da próxima vez que você vir um daqueles mapas azuis e vermelhos de uma galáxia girando, experimente enxergá-los como o que este estudo demonstrou que são: páginas de um manuscrito escrito antes da primeira estrela. Cada eixo de rotação é uma frase preservada do capítulo inicial do Universo, redigida pelas marés gravitacionais que esculpiram nossa vizinhança cósmica quando ela não passava de névoa de matéria quase uniforme. Durante 13,8 bilhões de anos, galáxias nasceram, colidiram, apagaram suas estrelas e seguiram girando, fiéis à torção original. Nós apenas aprendemos, finalmente, a ler o que sempre esteve escrito no movimento.