Imagine uma estrela com duzentas vezes a massa do Sol, nascida em um Universo que ainda não tinha completado 700 milhões de anos de existência e que não continha um único átomo de ferro. Ela brilha por cerca de dois milhões de anos, um sopro em escala cósmica, e então morre de um jeito que nenhuma estrela de hoje consegue morrer: em vez de deixar para trás um buraco negro ou uma estrela de nêutrons, desfaz-se por inteiro, lançando no espaço cinco massas solares de ferro recém-fabricado. Esse ferro se dissolve no gás de um pequeno halo de matéria escura, entra na composição das estrelas seguintes e, treze bilhões de anos depois, ainda pode ser lido no espectro de estrelas anãs que orbitam a Via Láctea. Essa é a história que um grupo internacional de astrofísicos acaba de contar com um dos conjuntos de simulações mais detalhados já produzidos sobre a infância da nossa galáxia.

O ponto de partida é um enigma que incomoda os astrônomos há mais de uma década. As galáxias anãs ultrafracas, conhecidas pela sigla em inglês UFD, de ultra-faint dwarfs, são os menores sistemas estelares que merecem o nome de galáxia: por definição, brilham menos que 20 mil sóis somados, e algumas delas, como Segue 1, emitem menos luz que algumas centenas de estrelas como a nossa. Em galáxias maiores vale uma regra simples: quanto mais massa em estrelas, mais ferro no gás e nas estrelas, porque sistemas massivos retêm melhor os produtos de suas supernovas. Abaixo de 100 mil massas solares em estrelas, porém, essa regra para de funcionar. Todas as anãs ultrafracas conhecidas se aglomeram em torno do mesmo valor de ferro, tenham elas 100 mil ou apenas algumas centenas de massas solares em estrelas.

Os astrônomos expressam esse conteúdo com a notação [Fe/H], o logaritmo da razão entre ferro e hidrogênio de uma estrela comparado ao mesmo valor medido no Sol. Zero significa composição solar; cada unidade negativa representa um fator de dez a menos. O platô das anãs ultrafracas fica em torno de −2,5, o que equivale a cerca de 300 vezes menos ferro que o Sol. Um valor de −3 corresponde a mil vezes menos, e −4 a dez mil vezes menos. O ferro é a moeda corrente da arqueologia galáctica porque produz linhas espectrais fáceis de medir em estrelas antigas, enquanto o oxigênio, mais abundante, domina as medições de gás em galáxias distantes feitas pelo Telescópio Espacial James Webb. Essa diferença de elementos é uma das razões pelas quais os dois mundos, o das estrelas fósseis e o do Universo jovem, conversam tão pouco entre si.

O estudo que agora propõe uma solução para o platô é liderado por Martin Rey, da Universidade de Bath, no Reino Unido, ao lado de Harley Katz, da Universidade de Chicago, e Corentin Cadiou, do Instituto de Astrofísica de Paris, à frente de uma equipe de 26 pesquisadores de 14 instituições em sete países, entre elas Oxford, Cambridge, Lund, Lyon, Stanford, Yonsei, Gante e a Escola Politécnica Federal de Lausanne. O artigo, cujo título em tradução livre é “MEGATRON: como as primeiras estrelas podem criar um platô de metalicidade de ferro nas menores galáxias anãs”, circula no repositório arXiv sob o código 2510.05232 e teve sua versão revisada, depois do parecer de árbitros, divulgada em setembro de 2026 no formato do Open Journal of Astrophysics. Ele integra uma série de trabalhos coordenados que apresentaram o projeto MEGATRON em outubro de 2025.

MEGATRON é o nome de um conjunto de simulações cosmológicas de radiação-hidrodinâmica, isto é, simulações que acompanham ao mesmo tempo a gravidade, o movimento do gás, a propagação da luz e a química do meio interestelar em uma região do Universo primitivo. A região escolhida não é qualquer uma: trata-se de um volume que, levado até o presente, colapsa em um halo de matéria escura de um trilhão de massas solares, abrigando uma galáxia espiral com um disco de cerca de 50 bilhões de massas solares em estrelas. Em outras palavras, uma proto-Via Láctea. As simulações seguem essa região desde o Big Bang até um redshift, ou desvio para o vermelho, de aproximadamente 8, a medida que os astrônomos usam para distância e tempo cósmico, e que nesse caso corresponde a um Universo com perto de 650 milhões de anos.

A resolução é o que separa o MEGATRON da maioria dos trabalhos anteriores. Cada partícula de matéria escura representa 25 mil massas solares, e a malha que descreve o gás chega a células de 2,5 parsecs no início da simulação e de cerca de 5 parsecs em seu fim, algo como 16 anos-luz, pouco mais que a distância do Sol às estrelas vizinhas mais próximas. Isso basta para resolver individualmente os berçários das primeiras estrelas. O código empregado é o RAMSES-RTZ, uma versão do consagrado RAMSES, de malha adaptativa, estendida para transportar radiação e resolver a química fora do equilíbrio de mais de 80 espécies, entre átomos primordiais, moléculas e íons de elementos pesados, todos acoplados na mesma passagem de cálculo.

Acoplar radiação e química faz diferença porque o gás primordial, feito só de hidrogênio e hélio, tem um único refrigerante eficiente para perder calor e colapsar em estrelas: a molécula de hidrogênio, H₂. Essa molécula é frágil. Fótons ultravioleta com energias entre 11,2 e 13,6 elétron-volts, a chamada banda de Lyman-Werner, dissociam o H₂ sem ionizar o hidrogênio atômico. Quando uma região fica banhada por esse tipo de radiação, os menores halos de matéria escura perdem a capacidade de resfriar seu gás e deixam de formar estrelas. Apenas halos mais massivos, com poços gravitacionais mais fundos e gás mais denso, seguem adiante. O MEGATRON calcula esse fundo de radiação a partir das próprias estrelas que nascem na simulação, em vez de impor um valor vindo de fora.

As primeiras estrelas, a chamada População III, são tratadas uma a uma. Quando uma nuvem de gás primordial atinge as condições de colapso, o código sorteia estrelas individuais de uma função de massa inicial em forma de sino logarítmico centrada em 100 massas solares, seguindo a prescrição proposta por John Wise e colaboradores em 2012. O destino de cada estrela depende de sua massa. Entre 10 e 20 massas solares, ela termina como supernova de colapso de núcleo, com energia de 10⁵¹ erg, o valor típico de uma supernova comum. Entre 20 e 40, vira hipernova, cerca de vinte vezes mais energética. Entre 140 e 300 massas solares, explode como supernova de instabilidade de pares, liberando perto de 5 × 10⁵² erg, cinquenta vezes a energia de uma supernova ordinária. Fora dessas faixas, a estrela colapsa diretamente em um buraco negro, sem devolver nada ao meio.

A supernova de instabilidade de pares, ou pair-instability supernova, abreviada PISN, é um fenômeno exclusivo de estrelas muito massivas e pobres em metais. No núcleo dessas estrelas a temperatura sobe a ponto de os fótons de raios gama se converterem em pares de elétrons e pósitrons. Cada conversão rouba energia da radiação que sustenta a estrela contra a própria gravidade; a pressão despenca, o núcleo desaba e a queima explosiva de oxigênio e silício que se segue despedaça a estrela por completo. Não sobra remanescente algum. Entre as cinzas há níquel-56 em abundância, que decai em cobalto e depois em ferro em questão de meses. E a produção de ferro cresce com a massa: segundo as tabelas de Ken’ichi Nomoto, Chiaki Kobayashi e Nozomu Tominaga, de 2013, adotadas na simulação, uma progenitora de 140 massas solares rende cerca de 0,5 massa solar de ferro, enquanto uma de 180 massas solares produz perto de 5 massas solares, dez vezes mais.

Assim que o gás de uma região ultrapassa um milionésimo da metalicidade solar, limiar em que a poeira passa a favorecer a fragmentação das nuvens, a simulação muda para o regime de População II, as estrelas de segunda geração, representadas por partículas de 500 massas solares que descrevem populações inteiras. É nesse regime que os autores introduziram variações deliberadas. São quatro simulações da mesma região com a mesma física de População III, mas com receitas diferentes para a formação estelar e o feedback, o retorno de energia das estrelas ao gás que as cercou, no caso da População II. A primeira, batizada de Efficient SF, usa supernovas de 10⁵¹ erg; a segunda, Bursty SF, multiplica essa energia por cinco e produz histórias de formação estelar mais explosivas e intermitentes.

A terceira, Varying IMF, adota uma função de massa inicial que muda com a densidade do ambiente, tornando-se mais rica em estrelas massivas em regiões densas, e permite que estrelas de baixa metalicidade explodam como hipernovas em metade dos casos. A quarta, chamada HN com eficiência máxima, fixa em 100% a eficiência de conversão de gás em estrelas por tempo de queda livre e agrupa as estrelas em partículas de 2 mil massas solares, formando aglomerados mais compactos. O propósito dessa coleção não é encontrar a receita certa, mas testar a robustez do resultado: se o platô de ferro aparecer nas quatro, apesar de diferenças tão grandes na física das estrelas de segunda geração, a sua origem terá de estar em outro lugar. Ao todo, cada simulação forma centenas de galáxias anãs, nunca menos de 200, e as quatro somam mais de mil anãs fracas.

Faltava construir a ponte entre o Universo de 650 milhões de anos e o céu de hoje. Para isso, a equipe rodou uma versão só com matéria escura das mesmas condições iniciais até o presente, com a mesma resolução de partículas, e identificou os halos com o localizador ROCKSTAR. Cada galáxia anã da simulação completa foi casada com seu halo correspondente na versão sem gás. Em seguida, os pesquisadores marcaram os 5% de partículas de matéria escura mais fortemente ligadas de cada halo e as seguiram por 13 bilhões de anos, técnica que Rey havia desenvolvido com Tjitske Starkenburg em 2022. Se metade dessas partículas marcadas terminar dentro de uma estrutura ligada em z = 0, a anã é considerada uma sobrevivente, um satélite da Via Láctea simulada; caso contrário, foi despedaçada e dissolvida no halo estelar difuso da galáxia hospedeira.

O primeiro resultado é a reprodução da relação entre massa estelar e metalicidade observada nas anãs reais. Entre 10 mil e um milhão de massas solares em estrelas, as galáxias simuladas se sobrepõem aos dados do Universo local compilados por Andrew Pace em 2025, com a mesma inclinação, o mesmo nível e a mesma dispersão. Três das quatro simulações concordam com as observações até o topo da faixa de massa estudada. A exceção é a Efficient SF, que acima de um milhão de massas solares produz galáxias com ferro demais. A causa é conhecida: seu retorno de energia fraco gera ventos galácticos tímidos, incapazes de expulsar o ferro para fora do meio interestelar. É a demonstração, pela primeira vez sobre centenas de objetos e não sobre um punhado deles, de que a relação massa-metalicidade das anãs é um sensor direto da força dos fluxos que saem das galáxias.

Abaixo de 100 mil massas solares, acontece o que os dados exigiam. A moda da distribuição, isto é, o valor em torno do qual a maioria das galáxias se concentra, deixa de cair com a massa e se estabiliza perto de −2,5, exatamente onde fica o platô observado. Nas quatro simulações, 78% das anãs fracas têm metalicidade média entre −3 e −2, intervalo que os autores definem como o platô de ferro. As demais 22% formam algo que as observações ainda não mostram: uma cauda longa de galáxias com menos de um milésimo do ferro solar, com alguns objetos descendo até −5. O fato de o platô surgir igual nas quatro rodadas, apesar das receitas tão distintas para a População II, apontava a seta para o único ingrediente que elas compartilham, a física das primeiras estrelas.

Para dissecar o mecanismo, os pesquisadores dividiram a amostra de anãs fracas em três grupos conforme a posição em relação ao platô: 8 galáxias acima dele, 876 sobre ele e 303 abaixo, 1.187 ao todo. Depois contaram quantas estrelas de População III haviam explodido dentro de cada uma. O resultado é quase monótono de tão uniforme: a imensa maioria das galáxias, em qualquer dos três grupos, passou por exatamente uma explosão primordial. Uma só detonação basta para elevar o gás do minihalo acima do limiar de metalicidade e disparar a formação de estrelas de segunda geração. Apenas 6% das anãs registraram várias explosões, em geral quando um pequeno aglomerado de estrelas primordiais nasceu quase ao mesmo tempo, com idades espalhadas por menos de um milhão de anos, ou quando dois minihalos independentes se fundiram mais tarde na mesma galáxia.

Houve ainda um grupo minúsculo, 2% do total, de galáxias que começaram a formar estrelas de População II sem jamais ter abrigado uma estrela primordial. Elas foram poluídas de fora, por ejeções diluídas de vizinhas que atravessaram o espaço intergaláctico. Todas essas anãs enriquecidas externamente se encontram abaixo do platô, mas são tão poucas que não explicam a cauda inteira de galáxias pobres em ferro. Sua distribuição espacial acompanha a teia cósmica como a de todas as outras, e a química de suas primeiras estrelas é indistinguível do restante da população. Pelas contas dos autores, apenas 3 dessas 22 anãs sobreviveriam ao campo de maré da galáxia hospedeira até hoje, de modo que algumas das satélites reais da Via Láctea podem ter nascido assim, embora não exista um jeito óbvio de reconhecê-las.

A chave estava no tipo de estrela que explodiu. Quando a equipe separou as explosões por progenitora, verificou que quase nenhuma anã do MEGATRON passou por supernovas de colapso de núcleo ou hipernovas de População III, consequência direta de uma função de massa inicial inclinada para as estrelas gigantes. Praticamente todas as explosões primordiais foram PISNs. E aí a divisão ficou nítida: as galáxias sobre o platô foram enriquecidas quase exclusivamente por progenitoras entre 160 e 300 massas solares, enquanto as galáxias abaixo do platô receberam uma contribuição considerável de progenitoras na faixa estreita entre 140 e 160 massas solares. A diferença de dez vezes no ferro produzido pelas duas classes decide, sozinha, de que lado da linha cada galáxia vai cair. Em ambos os casos a explosão gera carbono, oxigênio e outros elementos em quantidade suficiente para levar o gás ao regime de População II; só o ferro muda.

Restava um problema que já havia derrubado ideias parecidas no passado. Uma PISN libera cinquenta vezes mais energia que uma supernova comum, e simulações de anãs isoladas mostravam que minihalos de 100 mil a um milhão de massas solares simplesmente não seguram os metais de uma explosão dessas: o gás enriquecido é soprado para o meio intergaláctico e se perde. No MEGATRON, as estrelas primordiais nascem em halos bem maiores, com no mínimo 5 milhões de massas solares e uma mediana entre 10 e 30 milhões. O responsável é o fundo de radiação de Lyman-Werner gerado pela protogaláxia central, intenso o bastante para suprimir o hidrogênio molecular nos halos pequenos e empurrar a formação estelar primordial para cima na escala de massas. Halos de 10 milhões de massas solares têm poços gravitacionais fundos e gás denso o suficiente para reter os metais mesmo de explosões de 10⁵³ erg.

Com esses dois ingredientes, a conta do platô cabe em um guardanapo. Um halo de 10 milhões de massas solares carrega, pela fração cósmica de bárions, pouco mais de um milhão de massas solares de hidrogênio. Dissolva nesse reservatório as 5 massas solares de ferro de uma PISN de alta massa, misturando tudo uniformemente dentro do raio virial, e o gás resultante terá [Fe/H] próximo de −2,5. É o platô. Repita a operação com a meia massa solar de ferro de uma PISN de baixa massa e o resultado cai para −3,5; com as 0,07 massas solares típicas de uma supernova de colapso de núcleo, para −4,5. A analogia é a de uma piscina: a quantidade de corante define a cor, mas o tamanho da piscina também, e no MEGATRON quase todas as piscinas têm o mesmo tamanho porque a radiação da vizinhança impede que as menores existam.

Os autores fazem questão de qualificar essa aritmética como aproximada. Ela supõe que todo o ferro fica retido e se mistura por igual, o que não ocorre nem na simulação completa nem na natureza, onde a mistura irregular ajuda a explicar a largura da distribuição em torno do platô. E ela revela uma degenerescência incômoda: uma única supernova comum em um minihalo de 100 mil massas solares também produziria −2,5, porque a piscina seria cem vezes menor. Os efeitos ainda se acoplam uns aos outros. Uma função de massa inicial mais carregada de estrelas gigantes aumenta o rendimento de ferro, mas também intensifica o fundo de Lyman-Werner e empurra a formação estelar para halos maiores, o que dilui mais o ferro. O platô, nas palavras do artigo, não é ajustado a dedo nem é óbvio; emerge das escolhas feitas para a População III combinadas ao ambiente específico de uma proto-Via Láctea.

O ferro não é a única assinatura. A equipe comparou também as abundâncias médias de carbono, nitrogênio e magnésio das anãs simuladas com as de estrelas individuais medidas em anãs ultrafracas reais. Carbono e magnésio mostram concordância razoável na forma das distribuições, embora as galáxias simuladas apareçam deslocadas para baixo em cerca de 0,5 dex, meio fator de dez, no carbono. Esse deslocamento supera as incertezas típicas das observações, mas cabe com folga nas incertezas dos modelos de nucleossíntese das PISNs, nos quais o tratamento da convecção e da mistura interna pode alterar a produção de carbono por um fator de alguns. O nitrogênio é outra história. Nas anãs simuladas, sua abundância média é um milhão de vezes menor que a solar, várias ordens de grandeza abaixo das poucas detecções disponíveis em estrelas de anãs ultrafracas.

A escassez de nitrogênio não é um acidente da simulação, e sim a impressão digital esperada das PISNs, cujas explosões produzem elementos de número atômico par em abundância e quase nada dos ímpares. Essa alternância par-ímpar é o sinal mais claro de enriquecimento por estrelas primordiais gigantes, e alegações de tê-la encontrado em estrelas do halo da Via Láctea foram derrubadas por análises posteriores. Como as estrelas das anãs ultrafracas respondem por menos de 1% do halo galáctico, a ausência do padrão ali não descarta o cenário. Para fechar a conta do nitrogênio, bastaria um canal adicional modesto: ventos das próprias estrelas primordiais antes de explodirem, ou a rotação, que nos modelos de Koh Takahashi e colaboradores eleva a produção de nitrogênio das PISNs em quatro ordens de grandeza. As galáxias ricas nesse elemento que o James Webb vem encontrando no Universo jovem, como GN-z11, mostram que ele guarda lacunas em qualquer modelo.

Comparar galáxias simuladas aos 650 milhões de anos do Universo com galáxias observadas hoje significa saltar 12 bilhões de anos de evolução cósmica, e a equipe dedica uma seção inteira a justificar o salto. O argumento começa pelas observações: todas as anãs ultrafracas conhecidas estão apagadas e formaram a maior parte de suas estrelas há mais de 10 bilhões de anos. Em z = 8, porém, a reionização ainda não havia terminado em nenhuma das simulações, e boa parte das anãs do MEGATRON seguia formando estrelas. Se essas galáxias continuassem a crescer e a se enriquecer por mais alguns bilhões de anos, deixariam de ser comparáveis aos fósseis locais. Era preciso mostrar que os traços encontrados, platô e cauda, não eram um retrato passageiro de uma fase intermediária.

O primeiro teste isolou as anãs já apagadas pela reionização em z = 8: galáxias centrais que não formavam uma estrela sequer havia 100 milhões de anos. A seleção derruba a amostra para 3, 28, 17 e 51 galáxias nas quatro simulações, com o menor número justamente na Efficient SF, a menos reionizada de todas. Mesmo assim, a sobreabundância de galáxias sobre o platô permanece, e um teste estatístico de Kolmogorov-Smirnov não encontra diferença entre a distribuição de ferro dessas anãs e a da população completa. Elas se concentram perto do centro do volume, onde a protogaláxia produz a maior parte dos fótons ionizantes; seu gás está a 11 mil kelvins, 20% acima da mediana, com 30% menos hidrogênio neutro. A massa máxima que seus halos atingem em toda a história cósmica, 800 milhões de massas solares na mediana, fica abaixo do limiar para voltar a acumular gás quando o fundo ultravioleta se torna uniforme. Elas estão fossilizadas para sempre.

Há um segundo ajuste necessário. Estrelas perdem massa ao envelhecer, e as mais massivas simplesmente deixam de existir como estrelas em 13 bilhões de anos. Para cada anã, os autores calcularam a massa estelar projetada para o presente, contando apenas as estrelas que ainda estariam na sequência principal pelo modelo de evolução adotado. A correção desloca as massas estelares para baixo em 0,2 a 0,3 dex, ou seja, as galáxias perdem entre cerca de 40% e 50% de sua massa luminosa, mas a forma do extremo de baixa massa da relação massa-metalicidade não muda. Nos gráficos do artigo, os pontos marcam a massa em z = 8 e os losangos a massa projetada para hoje, lado a lado, com o platô intacto nos dois casos.

O segundo grande processo entre lá e cá é o campo de maré da galáxia central, que desmonta satélites como quem desfia um novelo. Pelo rastreamento das partículas marcadas, apenas 25% das anãs fracas presentes em z = 8 podem ser identificadas como um remanescente ligado em z = 0; o restante foi triturado e hoje compõe o halo estelar difuso da Via Láctea simulada. Os números de sobreviventes são parecidos nas quatro simulações, 66, 76, 92 e 70, como seria de esperar, já que a história gravitacional do objeto central é a mesma em todas. E a distribuição de ferro das sobreviventes é indistinguível da população-mãe, com um valor p de 0,55 no teste de Kolmogorov-Smirnov. Platô e cauda sobrevivem à destruição seletiva e chegam ao presente como previsões robustas, não como transientes do Universo jovem.

A ponte tem fissuras que os autores não escondem. O rastreamento em uma simulação só de matéria escura ignora efeitos bariônicos como o arrancamento de estrelas pelo potencial do disco galáctico e as transições entre perfis de densidade que tornam as anãs mais ou menos frágeis às marés. O plano é refazer o mapeamento sobre uma simulação hidrodinâmica das mesmas condições iniciais levada até o presente com o modelo Vintergatan, nome que significa Via Láctea em sueco, desenvolvido pelo grupo de Oscar Agertz em Lund. E não há garantia de que a protogaláxia simulada reproduza a história exata da nossa, de modo que o fundo de radiação e a cronologia do enriquecimento podem ter sido diferentes aqui. Simulações em ambientes controlados, construídas com a técnica de modificação genética das condições iniciais ou com reconstruções do Grupo Local, deverão preencher essa lacuna.

É a cauda, mais que o platô, que transforma o trabalho em uma previsão testável. Se a distribuição real de metalicidades das anãs ultrafracas se parecer com a do MEGATRON, cerca de uma em cada cinco deveria ter em média menos de um milésimo do ferro solar. Nenhuma das 38 anãs ultrafracas com medições independentes de ferro exibe uma média tão baixa, e a probabilidade de isso ocorrer por acaso, se 20% delas de fato pertencessem à cauda, é de 0,8 elevado a 38, perto de 0,02%. Tomado ao pé da letra, o número desfavorece a existência da cauda. Os autores, contudo, apontam que o censo de ferro das anãs ultrafracas é esparso em duas direções: faltam galáxias medidas e faltam estrelas medidas em cada galáxia.

Pictor II, uma anã descoberta em 2016 nas imediações da Grande Nuvem de Magalhães, ilustra o problema. Sua primeira caracterização química, publicada por Pace e colaboradores em 2025, encontrou média de −2,99, encostada no limite da cauda. Em seguida, Anirudh Chiti e colegas identificaram na mesma galáxia uma estrela com [Fe/H] abaixo de −4, dez mil vezes menos ferro que o Sol, a primeira ultrapobre em ferro registrada em uma anã ultrafraca, e que nem entrou na média anterior. O caso mostra que teores tão baixos existem nessas galáxias e que médias calculadas sobre poucas estrelas podem estar enviesadas para cima. Nas anãs simuladas, mesmo as que ficam sobre o platô têm dispersão interna de cerca de 0,4 dex, idêntica à medida nas anãs reais, com estrelas individuais descendo a −3 e a −4. A cauda do MEGATRON pode ser apenas a extensão, à escala de galáxias inteiras, de uma cauda de estrelas que já está lá, escondida pelo número pequeno de espectros por galáxia.

A própria área é jovem. Até 2005, as galáxias satélites conhecidas da Via Láctea cabiam em uma lista de pouco mais de uma dezena, todas descobertas em placas fotográficas. Foi o Sloan Digital Sky Survey que revelou as primeiras anãs ultrafracas, começando por Willman 1 e Ursa Major I naquele ano, seguidas por Boötes I, Coma Berenices, Segue 1 e outras até 2007. O Dark Energy Survey acrescentou uma leva de candidatas em 2015, entre elas Reticulum II, e levantamentos como Pan-STARRS, DELVE e o Hyper Suprime-Cam do telescópio Subaru continuaram a engrossar o catálogo, que hoje soma cerca de seis dezenas de satélites confirmados ou candidatos. São objetos tão tênues que seus membros precisam ser separados das estrelas do primeiro plano um a um, com ajuda decisiva dos movimentos próprios medidos pelo satélite Gaia.

O que torna essas galáxias tão valiosas é justamente a razão de sua pobreza. Quando a reionização cósmica, completada por volta de um bilhão de anos após o Big Bang, aqueceu o gás intergaláctico a dezenas de milhares de kelvins, os halos mais leves perderam a capacidade de atrair gás novo, porque seus poços gravitacionais eram rasos demais para capturar um meio tão quente. A formação estelar cessou e nunca mais recomeçou. Desde então, essas galáxias são cápsulas do tempo, estratos de uma escavação arqueológica que preservam a química do Universo com menos de um bilhão de anos. Ler o ferro de suas estrelas é ler a história das primeiras explosões; por isso elas se tornaram o laboratório preferido de quem quer saber como eram as estrelas de População III, que ninguém jamais observou diretamente.

Modelos anteriores haviam tentado explicar o platô por caminhos diferentes. Myoungwon Jeon, Gurtina Besla e Volker Bromm já haviam apontado em 2017 para a física da População III como fonte do excesso de ferro. Outros grupos sugeriram uma função de massa inicial variável para as estrelas de segunda geração, um piso de metalicidade imposto ao meio intergaláctico por enriquecimento externo, ou detalhes do transporte de ferro pelos ventos galácticos. Simulações de anãs isoladas, como as de Mahsa Sanati e colaboradores em 2023, viam os metais das PISNs escaparem rapidamente para o meio intergaláctico. A diferença decisiva do MEGATRON é o ambiente: uma região superdensa que colapsa cedo e abriga mais estrelas do que a média cósmica para sua massa de halo, o que eleva permanentemente o fundo de Lyman-Werner e muda o lugar onde as primeiras estrelas podem nascer. Simular a anã sozinha, sem a protogaláxia ao redor, perdia o ingrediente principal.

Nada disso elimina as incertezas sobre as primeiras estrelas, e o artigo as enumera com franqueza. A função de massa inicial primordial adotada favorece estrelas de 100 massas solares ou mais, em linha com as simulações pioneiras de Tom Abel e Volker Bromm no início dos anos 2000, que mostravam nuvens primordiais colapsando em grumos de gás dessa massa. Trabalhos mais recentes encontram fragmentação adicional por turbulência e instabilidades durante o colapso, o que reduziria a massa característica, enquanto campos magnéticos amplificados e a radiação das próprias estrelas recém-nascidas podem inibir essa fragmentação e elevá-la de novo. A faixa de massas que produz PISNs também é disputada, com limiares inferiores entre 140 e 180 massas solares conforme o modelo, e tabelas alternativas de rendimento dão menos ferro para progenitoras de baixa massa, o que produziria uma cauda mais pronunciada e um platô mais tímido.

A fragilidade vira oportunidade quando se inverte o raciocínio. Se o platô é feito de PISNs de alta massa e a cauda, de PISNs de baixa massa, a razão entre o número de anãs em cada grupo mede diretamente a razão entre os dois tipos de progenitoras, isto é, a forma da função de massa inicial das primeiras estrelas. Um censo completo de ferro nas anãs ultrafracas se transformaria em uma régua para a distribuição de massas da População III, algo que nenhuma observação direta conseguiu oferecer até hoje. Se a cauda não existir, a função de massa teria de ser ainda mais carregada de gigantes, produzindo apenas PISNs de alta massa, ou o platô teria outra origem, como supernovas comuns em minihalos menores. Em qualquer dos casos, a resposta estaria escrita nas estrelas mais apagadas do céu.

As quatro simulações também divergem na velocidade com que reionizam seu entorno, e esse detalhe explica por que o número de anãs apagadas varia tanto entre elas. Todas estão mais ionizadas que a média cósmica inferida das observações, como se espera de uma região superdensa, mas a Efficient SF, que mais fótons ionizantes produz graças à sua formação estelar vigorosa, é a que menos reioniza, porque seu meio interestelar denso aprisiona a radiação. A HN com eficiência máxima, no extremo oposto, já está completamente reionizada em z = 8. A fuga dos fótons ionizantes é regulada pelo retorno mecânico das supernovas, que abre canais no gás, mais do que pela quantidade de luz produzida, em sintonia com outros estudos numéricos da última década.

O horizonte observacional está a favor dos autores. O Observatório Vera C. Rubin, já em operação no Chile, deverá revelar dezenas de novas candidatas a anãs ultrafracas em seu levantamento de dez anos, e a próxima geração de telescópios gigantes, com espelhos de 25 a 39 metros, terá a sensibilidade para obter espectros de estrelas individuais nesses objetos e medir seu ferro, carbono e outros elementos. Do lado teórico, a equipe planeja ressimular o mesmo volume com uma função de massa inicial primordial diferente, rastrear elementos relevantes para a espectroscopia estelar, como sódio, alumínio, manganês, cobre e zinco, e seguir todas as estrelas pobres em metais do volume, e não só as que terminam dentro de anãs. O objetivo declarado é usar ao mesmo tempo as linhas de emissão do James Webb e a arqueologia das estrelas locais para apertar o cerco à física das primeiras explosões.

Rodar o MEGATRON exigiu os maiores recursos computacionais acadêmicos do Reino Unido, com a infraestrutura DiRAC em Leicester e Durham e o serviço CSD3 de Cambridge, além dos aglomerados Infinity, em Paris, e Glamdring, em Oxford, para o processamento dos dados. Os autores declaram que os dados da simulação serão compartilhados sob pedido razoável, acompanhados da documentação e dos códigos de análise, e o trabalho foi construído sobre ferramentas de código aberto como pynbody, yt, NumPy e matplotlib. O projeto contou com financiamento da NASA, do Conselho Europeu de Pesquisa e de agências do Reino Unido, da Suécia, da Suíça, da Bélgica e da Coreia do Sul.

No fim, a imagem que fica é a de uma linha reta em um gráfico. Durante anos, ela foi um incômodo, um platô que nenhum modelo sabia desenhar. Agora pode ser lida como um epitáfio coletivo: cada ponto sobre ela seria uma galáxia que viveu uma única grande morte, a de uma estrela primordial com mais de 160 vezes a massa do Sol, dissolvida em um halo de dez milhões de massas solares exatamente grande o bastante para não deixar escapar o que ela produziu. O ferro que hoje tinge os espectros das estrelas mais apagadas da vizinhança da Via Láctea seria, em grande parte, o ferro dessa única explosão. Treze bilhões de anos depois, a primeira geração de estrelas, que ninguém jamais viu, ainda assina seu nome nas menores galáxias do Universo.