Acerca de 212 mil anos-luz da Terra, na direção da constelação do Boieiro, existe um borrão de luz tão tênue que ninguém o enxergou até 2006, quando um levantamento digital do céu o arrancou do ruído estatístico. Boötes I não tem braços espirais, não tem gás, não tem berçários estelares. O que tem é um punhado de estrelas velhíssimas, nascidas quando o universo ainda não completava o primeiro bilhão de anos, e uma quantidade de matéria escura que supera em cerca de cem vezes a massa de todas as suas estrelas somadas. É um fóssil. E foi nesse fóssil que o Telescópio Espacial James Webb acaba de contar, uma por uma, mais de dez mil estrelas, descendo até astros com apenas 15% da massa do Sol, para responder a uma pergunta que atravessa toda a astrofísica: o universo sempre fabricou estrelas do mesmo jeito?

A pergunta tem nome técnico. Astrônomos chamam de função de massa inicial, ou IMF, sigla em inglês para initial mass function, a distribuição que descreve quantas estrelas de cada massa nascem numa mesma geração. Pense numa fornada de pães saindo do forno: a receita diz quantos pãezinhos pequenos e quantos pães grandes saem a cada leva. Em 1955, Edwin Salpeter mediu essa receita na vizinhança solar e encontrou uma lei simples, uma potência com expoente de −2,35, que significa que estrelas pequenas são muitíssimo mais numerosas que as grandes. Para cada estrela com dez vezes a massa do Sol, nascem centenas de astros menores que ele. Desde então, a forma dessa curva nas massas mais baixas virou uma das medidas mais disputadas da astronomia.

Essa curva não é uma curiosidade de catálogo. Na ponta das massas altas, ela determina quantas supernovas explodirão numa galáxia, quantas explosões de raios gama ela produzirá e quantos buracos negros acabarão se fundindo. Na ponta das massas baixas, ela governa algo ainda mais básico: a estimativa da massa estelar de qualquer galáxia que não conseguimos resolver em estrelas individuais, o que inclui praticamente todas as galáxias distantes fotografadas pelo Webb nos últimos anos. Quando um astrônomo anuncia que uma galáxia do universo primitivo tem tantos bilhões de massas solares em estrelas, ele está assumindo uma IMF, em geral a da Via Láctea. Se a receita for diferente lá longe e lá atrás no tempo, essas massas podem estar erradas por fatores consideráveis.

Há bons motivos teóricos para suspeitar que a receita mude. Estrelas nascem do colapso de nuvens de gás frio, e o tamanho dos fragmentos que colapsam depende de quão bem a nuvem consegue se livrar do calor. Em gás rico em elementos pesados, como o das nuvens atuais da Via Láctea, grãos de poeira irradiam a energia com eficiência e permitem que o gás se quebre em pedaços pequenos, favorecendo estrelas de baixa massa. Em gás quase puro de hidrogênio e hélio, como o das primeiras galáxias, o resfriamento depende de moléculas de hidrogênio, muito menos eficazes. A física da chamada massa de Jeans, o limiar acima do qual a gravidade vence a pressão térmica, prevê fragmentos maiores e, portanto, menos estrelas pequenas nos ambientes pobres em metais.

Nem todos os modelos concordam. Simulações de alta resolução conduzidas por grupos como o de Matthew Bate, em várias rodadas entre 2014 e 2019, e trabalhos mais recentes de Tanvir e Krumholz, em 2024, sugerem que a dependência da IMF em relação à metalicidade é fraca, difícil de detectar e passível de ser apagada por efeitos dinâmicos aleatórios dentro dos próprios berçários, como ejeções de protoestrelas e competição por gás. Entre uma teoria que prevê variação e outra que prevê invariância, a única arbitragem possível é a observacional: encontrar um lugar onde estrelas tenham nascido de gás muito pobre em metais e contar quantas de cada massa sobreviveram até hoje.

A contagem é viável por um motivo elegante. Estrelas com menos de oito décimos da massa do Sol queimam seu hidrogênio com tamanha parcimônia que vivem mais do que a idade atual do universo. Nenhuma delas morreu ainda. Uma população que nasceu há 13 bilhões de anos preserva intacta, na faixa de baixa massa, a proporção com que foi fabricada: basta contar as estrelas que estão lá e converter o brilho de cada uma em massa, com a ajuda de modelos de evolução estelar. Antes do Telescópio Espacial Hubble, isso só era possível dentro da Via Láctea, porque nenhum outro instrumento alcançava a profundidade necessária em sistemas externos. O Hubble abriu a porta extragaláctica, e o Webb a escancarou.

Os primeiros resultados fora da nossa galáxia foram ambíguos. Em 2002, Rosemary Wyse e colaboradores mediram com o Hubble a IMF da galáxia anã de Ursa Menor e a encontraram indistinguível da dos aglomerados globulares da Via Láctea, apesar da quantidade colossal de matéria escura do sistema. Em 2013, Jason Kalirai e colegas obtiveram na Pequena Nuvem de Magalhães uma lei de potência com expoente −1,9, mais rasa que a de Salpeter, sem conseguir detectar nenhum pico na distribuição. O próprio Webb reacendeu a disputa na mesma galáxia: em 2025, uma análise da periferia da Pequena Nuvem concluiu que nenhuma quebra na curva era necessária, enquanto outra, publicada em 2026, encontrou uma quebra em 0,49 massa solar. A resposta, ao que parece, dependia muito de onde e de quão fundo se olhava.

É aqui que entram as galáxias anãs ultratênues, ou ultrafaint dwarf galaxies, sistemas tão apagados que alguns deles brilham menos do que uma única estrela supergigante. Elas são os objetos mais pobres em metais, mais antigos e mais dominados por matéria escura que se conhece. Modelos cosmológicos elaborados por Mia Bovill e Massimo Ricotti, um dos coautores do novo estudo, indicam que essas galáxias se formaram antes ou durante a época da reionização, quando a radiação das primeiras estrelas e galáxias ionizou o hidrogênio do cosmos e aqueceu o gás a ponto de sufocar a formação estelar nos sistemas menores. Elas pararam de fabricar estrelas ali, há mais de 12 bilhões de anos, e desde então carregam, congeladas, as propriedades das galáxias de alto desvio para o vermelho.

Há uma segunda vantagem, menos óbvia. Num aglomerado globular, as estrelas estão tão apertadas que trocam energia em encontros gravitacionais ao longo de bilhões de anos; as mais pesadas afundam para o centro e as mais leves são empurradas para fora ou expulsas, num processo chamado segregação de massa que distorce a distribuição original. Nas anãs ultratênues, a densidade estelar é tão baixa e a matéria escura tão dominante que o tempo necessário para esses encontros alterarem alguma coisa ultrapassa a idade do universo. O sistema é, na prática, sem colisões. A distribuição de massas que se mede hoje é a distribuição com que as estrelas nasceram, sem nenhum embaralhamento dinâmico pelo caminho.

Por isso essas galáxias viraram alvos preferenciais do Hubble na última década. Em 2013, Marla Geha e colaboradores mediram a IMF de Hércules e Leo IV até cerca de meia massa solar e encontraram inclinações mais rasas que a de Salpeter. Em 2018, Mario Gennaro, hoje no Space Telescope Science Institute e segundo autor do novo trabalho, reanalisou esses dados e acrescentou mais quatro anãs, chegando a 0,45 massa solar e notando uma tendência discreta de inclinações mais rasas em sistemas mais pobres em metais. Em 2022, Safarzadeh, Simon e Loeb estudaram Reticulum II e concluíram que uma lei de potência rasa era necessária para explicar a função de massa. Nenhum desses estudos, contudo, conseguiu rejeitar estatisticamente a IMF da Via Láctea.

O problema de todos eles era o mesmo: a profundidade. A IMF da Via Láctea, na descrição log-normal proposta por Gilles Chabrier em 2003, tem um pico característico em torno de 0,2 massa solar, e nenhuma dessas observações chegava perto disso. Tentativas posteriores de Gennaro, em Coma Berenices, e de Carrie Filion, em Boötes I em 2022 e em outras quatro anãs em 2024, alcançaram entre 0,3 e 0,2 massa solar, mas com amostras pequenas e incertezas grandes. Kareem El-Badry, também coautor do novo estudo, havia mostrado em 2017 que ajustar uma lei de potência única a uma distribuição que na verdade tem um pico produz inclinações artificialmente mais rasas quanto mais fundo se observa, o que explicaria a dispersão de resultados na literatura sem exigir nenhuma variação física.

Boötes I reunia as condições para resolver o impasse. Descoberta em 2006 por Vasily Belokurov e colaboradores nos dados do Sloan Digital Sky Survey, ela está a 65 quiloparsecs de distância e, com magnitude absoluta visual de −6,3, é uma das mais luminosas da categoria, o que significa mais estrelas para contar. Medidas espectroscópicas recentes com o espectrógrafo DEIMOS, no telescópio Keck, fixaram sua razão massa-luminosidade em 165 ± 45, uma assinatura de domínio da matéria escura, e sua metalicidade média em [Fe/H] = −2,4 ± 0,05. Traduzindo a notação: suas estrelas contêm cerca de 250 vezes menos ferro, em proporção ao hidrogênio, do que o Sol. A distribuição de metalicidades é uma curva em forma de sino, com pico em −2,4 e dispersão de 0,46.

A história química simples, documentada por Anna Frebel em 2016 e confirmada em 2026, indica que a galáxia teve poucos episódios de enriquecimento antes de parar. Thomas Brown, do Space Telescope Science Institute, mostrou em 2014, com imagens profundas do Hubble, que sua formação estelar foi extinta pela reionização. Trabalhos independentes estimam que cerca de 90% de suas estrelas nasceram antes da primeira queda da galáxia no campo gravitacional da Via Láctea, o que a protege de contaminação por gerações posteriores. Massimo Ricotti e Nickolay Gnedin a classificaram, já em 2005, como um fóssil verdadeiro: um pedaço do universo com menos de um bilhão de anos que sobreviveu, quase intocado, até a vizinhança da nossa galáxia.

O programa de observação, identificado como PID 3849, apontou a câmera NIRCam do Webb para os mesmos campos que Brown havia fotografado com o Hubble. A NIRCam registra simultaneamente dois comprimentos de onda, e a equipe escolheu o filtro F150W, centrado em 1,5 micrômetro, no canal de ondas curtas, e o filtro F322W2, uma banda larga em torno de 3,2 micrômetros, no canal de ondas longas. A combinação não foi casual: a diferença de brilho entre as duas bandas, que os astrônomos chamam de cor, separa com nitidez as estrelas da galáxia das galáxias de fundo, muito mais avermelhadas nessa faixa. O mosaico cobriu três fileiras e duas colunas de campos, com vinte posições de apontamento levemente deslocadas e um tempo total de exposição de 18.896,715 segundos, pouco mais de cinco horas e um quarto.

As imagens passaram pela versão 1.17.1 do processamento oficial do Webb, com ajustes específicos para remover o ruído eletrônico de baixa frequência dos detectores e os reflexos difusos conhecidos como wisps, manchas de luz espalhada que afetam quatro dos detectores da câmera. A fotometria, isto é, a medição do brilho de cada fonte, foi feita com o pacote DOLPHOT, que ajusta a cada estrela o perfil de luz característico do instrumento, a chamada point spread function, e por isso consegue separar fontes próximas umas das outras. A equipe dividiu o trabalho em 48 tarefas independentes, uma para cada detector do canal de ondas curtas em cada campo, depois de perceber que uma imagem combinada introduzia pequenas variações artificiais entre os detectores.

O resultado bruto foi um catálogo de 1.376.174 fontes detectadas, a esmagadora maioria ruído, galáxias de fundo ou medições defeituosas. Para separar o joio do trigo, a equipe aplicou uma bateria de critérios de qualidade: a classificação do objeto pelo DOLPHOT, que distingue estrelas de fontes alongadas ou extensas; a nitidez, que compara o tamanho angular da fonte com o perfil de luz do instrumento e denuncia galáxias; o parâmetro de aglomeração, que mede quanto o brilho de uma estrela seria contaminado pelas vizinhas; e limiares de relação sinal-ruído de 12 no filtro F150W e de 6 no F322W2. Um último crivo empírico rejeitou fontes cuja incerteza fotométrica destoava da curva típica para seu brilho. Sobraram 12.889 fontes pontuais confiáveis, menos de 1% do catálogo original.

Essas fontes foram dispostas num diagrama cor-magnitude, ou color-magnitude diagram, o retrato de família da astronomia estelar, em que cada estrela é um ponto posicionado pelo brilho e pela cor. Nele, as estrelas de Boötes I desenham uma faixa estreita e inclinada, a sequência principal, onde astros queimam hidrogênio no núcleo. As galáxias de fundo não resolvidas se agrupam numa nuvem avermelhada, com cor F150W − F322W2 acima de 0,8, exatamente como previam modelos construídos a partir do Campo Ultraprofundo do Hubble. As estrelas da Via Láctea que estão no caminho, anãs M e anãs marrons do halo segundo a simulação TRILEGAL, formam uma faixa vertical em torno de 0,6. Cortes em cor e brilho isolaram 11.046 candidatas a membros, distribuídas de maneira quase uniforme pelo campo, com variação de densidade de apenas 3% entre detectores.

Contar estrelas, porém, exige saber quantas se deixou de contar. Nenhum levantamento é completo: estrelas muito fracas se perdem no ruído, outras caem sobre pixels defeituosos ou nas bordas das imagens. Para medir essas perdas, a equipe recorreu aos testes de estrelas artificiais, ou artificial star tests: injetou nas próprias imagens 1,2 milhão de estrelas falsas com brilhos e cores conhecidos, um milhão delas distribuídas ao longo da sequência principal de Boötes I e outras 200 mil espalhadas uniformemente pelo diagrama, e depois rodou sobre elas exatamente o mesmo processamento aplicado aos dados reais. A fração recuperada em cada brilho fornece a completude; a diferença entre o brilho injetado e o medido fornece os erros fotométricos reais, incluindo os efeitos de aglomeração.

A completude cai a 50% na magnitude 27,9 no filtro F150W e em 27,3 no F322W2. Nos modelos estelares adotados, isso corresponde a estrelas de cerca de 0,15 massa solar, ou 150 vezes a massa de Júpiter, não muito acima do limiar de 0,08 massa solar abaixo do qual um objeto nem sequer consegue fundir hidrogênio. É um recorde para uma anã ultratênue. A variação da completude entre campos e detectores ficou em 4%, e pequenos desvios de calibração entre os detectores, de 0,005 a 0,022 magnitude, foram corrigidos com um procedimento estatístico de reamostragem repetido cinco mil vezes. Como a galáxia é esparsa, com dez mil estrelas espalhadas por detectores de mais de quatro milhões de pixels cada, a profundidade é limitada pelo ruído, e não pelo empilhamento de fontes.

Havia ainda um obstáculo nos modelos. Para converter brilhos em massas, os astrônomos usam isócronas, curvas teóricas que descrevem como aparece, no diagrama cor-magnitude, uma população de estrelas de mesma idade e composição. Quando a equipe sobrepôs às observações as isócronas da biblioteca BaSTI com a metalicidade de Boötes I, o ponto em que as estrelas deixam a sequência principal, o turn-off, apareceu mais fraco do que o previsto. Para reproduzi-lo, seria preciso uma idade maior que a do universo. A explicação está no oxigênio. Estrelas muito pobres em ferro costumam ter oxigênio em excesso, e esse elemento acelera o ciclo CNO de fusão, encurtando a vida na sequência principal: um aumento de 0,5 dex em [O/Fe] faz uma isócrona parecer cerca de um bilhão de anos mais jovem.

A solução foi adotar os modelos de Victoria–Regina enriquecidos em oxigênio que Brown já havia usado em 2014, válidos para estrelas acima de 0,4 massa solar, e costurá-los aos modelos BaSTI abaixo desse limite, com uma transição suave entre 0,4 e 0,6 massa solar controlada por uma função matemática que garante continuidade de valor e de inclinação. Os brilhos sintéticos nos filtros do Webb vieram de espectros das atmosferas MARCS. A população foi tratada como uma única geração de 13,3 bilhões de anos, nascida num surto breve, a uma distância de 65 quiloparsecs, com extinção interestelar quase nula, e com metalicidades sorteadas da distribuição espectroscópica do Keck, fixa em torno de −2,4.

O coração do trabalho é um código chamado Starwave, desenvolvido por Gennaro, pela primeira autora Keyi Ding, da Universidade de Maryland, e por Vedant Chandra, de Harvard. Ele funciona por modelagem direta, ou forward modeling: em vez de tentar extrair a IMF das observações por algum atalho analítico, o programa sorteia um conjunto de parâmetros, fabrica uma galáxia sintética com eles, aplica a essa galáxia os erros e as perdas medidos nos testes de estrelas artificiais e compara o diagrama cor-magnitude resultante com o real. Repetindo o processo milhares de vezes, aprende quais parâmetros produzem diagramas parecidos com o observado. A novidade está em como essa comparação é feita: por inferência baseada em simulação, ou simulation-based inference, em que uma rede neural aprende diretamente a distribuição de probabilidade dos parâmetros a partir das simulações.

Comparar dois diagramas com milhares de pontos cada um não é trivial, porque o número de estrelas varia e a ordem delas não importa. A equipe resolveu isso com uma impressão digital matemática: cada estrela é projetada num espaço abstrato por uma função de similaridade, e a soma das projeções de todas as estrelas vira um vetor de cinquenta números que resume o diagrama inteiro, insensível à ordem e ao tamanho da amostra. Dois diagramas parecidos têm impressões digitais próximas. Antes das comparações, duas calibrações empíricas aproximaram modelo e dados: uma correção de cor dependente do brilho, de no máximo 0,01 magnitude, que alinhou as cristas das sequências principais sintética e observada a menos de 0,005 magnitude, e um ruído adicional de 0,015 a 0,035 magnitude para reproduzir a largura real da sequência.

A inferência rodou em cinco rodadas de três mil simulações cada, com a rede neural refinando a cada etapa a região do espaço de parâmetros onde concentra os sorteios seguintes. Para garantir que o método não inventa resultados, a equipe o submeteu a testes de convergência: gerou galáxias sintéticas com parâmetros conhecidos, fez o Starwave tentar recuperá-los e varreu uma grade de configurações da impressão digital até encontrar aquela que devolvia os valores verdadeiros dentro das incertezas. As incertezas que sobram ao final desses testes são intrínsecas, isto é, não podem ser reduzidas com o método e os dados disponíveis, e por isso servem como piso para os erros do resultado real.

Três formas matemáticas para a IMF foram colocadas à prova. A primeira é a lei de potência única de Salpeter, uma rampa reta que sobe sem parar na direção das massas menores. A segunda é uma lei de potência quebrada, no espírito do modelo de Pavel Kroupa de 2001, uma rampa que muda de inclinação numa massa de quebra. A terceira é a função log-normal de Chabrier, um morro com um pico numa massa característica e uma largura. Em todos os casos a equipe mediu a IMF de sistemas, tratando cada binária não resolvida como um único objeto com a massa somada, e deixou livres tanto o número total de estrelas quanto a fração de binárias, porque ambos se misturam com a forma da curva.

A rampa reta foi a primeira a tropeçar. O melhor ajuste de lei de potência única tem expoente −1,60 ± 0,04, incompatível com os −2,35 de Salpeter com 99% de confiança. Os autores são cautelosos sobre o significado disso: a lei de Salpeter foi derivada de uma amostra que nem chegava perto das massas aqui observadas, e uma potência única diverge nas massas baixas, prevendo um número infinito de estrelas minúsculas. O que o resultado mostra é precisamente o efeito antecipado por El-Badry: quando a distribuição verdadeira tem um pico e a observação passa abaixo dele, a reta forçada sobre a curva fica rasa. A inclinação de −1,60 não descreve uma física diferente; descreve um modelo inadequado.

A lei de potência quebrada se saiu bem melhor. Abaixo da massa de quebra, a inclinação ficou em −1,31 ± 0,07; acima, em −1,99, com incertezas assimétricas de +0,16 e −0,22; e a quebra se assentou em 0,47 massa solar, com erros de +0,09 e −0,07. Os valores são estatisticamente consistentes com os da Via Láctea. A equipe restringiu a priori a massa de quebra ao intervalo de 0,35 a 0,65 massa solar e explica o motivo com franqueza: quando a quebra é deixada livre para percorrer toda a faixa de massas observada, o ajuste a empurra para as bordas e o modelo degenera numa potência única com a segunda inclinação solta. Um teste com dados sintéticos confirmou a degenerescência, produzindo uma distribuição com dois picos nas extremidades do intervalo permitido.

O resultado central veio da função log-normal. A massa característica, o pico do morro, ficou em 0,18 ± 0,02 massa solar, e a largura em 0,56 ± 0,06. Os valores de referência da Via Láctea, medidos por Chabrier na vizinhança solar, são 0,20 massa solar e 0,55. A concordância cabe dentro do intervalo de 68% de confiança, e a incerteza na massa característica, de 11%, é a menor já obtida para uma galáxia anã ultratênue. Uma galáxia com 250 vezes menos ferro que o Sol, que parou de fabricar estrelas antes de o universo completar seu primeiro bilhão de anos, usou a mesma receita que as nuvens moleculares da nossa vizinhança usam hoje. A fração de binárias ficou entre 51% e 60%, dependendo do modelo.

Para checar a qualidade dos ajustes, os autores sortearam mil conjuntos de parâmetros da distribuição final de cada modelo, geraram mil galáxias sintéticas e compararam a função de luminosidade acumulada de cada uma, isto é, a contagem de estrelas até cada brilho, com a observada. As curvas das leis quebrada e log-normal acompanham os dados ao longo de toda a faixa de magnitudes, com desvios de poucos por cento. A curva da potência única se afasta bem mais no trecho das estrelas fracas, entre as magnitudes 26 e 28, porque uma rampa que sobe sem limite na direção das massas menores não consegue acompanhar a inflexão que os dados mostram em torno da massa característica. É nessa faixa, exatamente a que só o Webb conseguiu alcançar, que a diferença entre os modelos se decide.

A comparação visual foi além da contagem por brilho. Os autores construíram histogramas bidimensionais, dividindo o diagrama cor-magnitude em celas pequenas, e calcularam em cada uma a diferença entre o número de estrelas observadas e o número previsto, normalizada pelo ruído estatístico esperado. O resultado é um mapa colorido do acerto do modelo: amarelo onde ele reproduz os dados, azul onde o céu tem mais estrelas do que o modelo prevê, vermelho onde o modelo coloca estrelas que não existem. Como as celas são menores que a incerteza fotométrica típica, nenhum modelo acerta cela por cela, e os três mapas parecem mosaicos ruidosos. Mas, no mapa da potência única, as discrepâncias se concentram ao longo da base da sequência principal, justamente onde a forma da curva mais importa.

O teste quantitativo usou a própria impressão digital dos diagramas. Para cada modelo, a equipe mediu a distância média entre a impressão digital do diagrama observado e as das galáxias sintéticas sorteadas da distribuição final, e comparou esse número com a distribuição de distâncias entre as próprias galáxias sintéticas, num critério herdado de trabalhos anteriores de Gennaro. Se o modelo descreve bem os dados, o diagrama real se comporta como mais uma simulação, nem mais distante nem mais próxima das demais. Foi o que aconteceu com a lei quebrada e com a log-normal, cuja consistência com as observações ficou mais apertada que o nível de 68%. A potência única, ao contrário, caiu fora dessa faixa: o diagrama real ficou mais longe das simulações do que a maioria delas fica umas das outras.

Faltava mostrar que o método conseguiria, de fato, distinguir uma rampa reta de um morro com os dados disponíveis. A equipe gerou cem galáxias sintéticas a partir de cada um dos dois modelos, ajustou a cada uma tanto a potência única quanto a log-normal e verificou qual das duas o procedimento preferia. Quando a galáxia havia sido fabricada com uma potência única, o método reconheceu a forma correta em 97% dos casos; quando havia sido fabricada com a log-normal, acertou em 100%. A matriz de confusão resultante é quase diagonal. Somado ao mau ajuste da potência única aos dados reais, o teste indica que a preferência por uma IMF com pico característico não é artefato do método, mas uma propriedade da população estelar de Boötes I.

Os números conversam com os estudos anteriores da mesma galáxia. A análise de Gennaro de 2018, limitada a 0,45 massa solar, havia encontrado uma potência única de −1,84 e parâmetros log-normais de 0,33 massa solar e largura 0,59. Filion, em 2022, chegando a 0,3 massa solar, obteve −1,95 para a potência única, com incertezas próximas de uma unidade, e massa característica de 0,17 massa solar, ajustando porém a IMF de estrelas individuais, e não de sistemas, o que impede comparação numérica direta. A conclusão qualitativa, essa sim, é a mesma: sempre que se adota um modelo com pico, Boötes I se parece com a Via Láctea. A diferença é que, agora, a observação atravessa o pico em vez de extrapolar até ele, e as barras de erro encolheram a uma fração das anteriores.

A comparação com as outras anãs ultratênues estudadas em profundidade pelo Hubble, Hércules, Leo IV, Coma Berenices, Canes Venatici II, Ursa Maior I, Reticulum II e Ursa Maior II, revela um padrão claro. Todas apresentam inclinações de potência única mais rasas que a de Salpeter, entre −0,8 e −2,2, e essas inclinações parecem acompanhar a profundidade de cada observação, o que reforça a leitura de que o achatamento é um efeito da faixa de massas observada, e não uma variação física entre galáxias. As massas características log-normais, por sua vez, vão de 0,15 a 0,62 massa solar, mas com incertezas tão grandes nos sistemas menos profundos que a comparação direta entre elas tem valor limitado.

Quando a massa característica de cada galáxia é colocada contra sua metalicidade, surge uma tendência tímida: sistemas mais pobres em metais parecem ter picos em massas mais altas, como prevê a teoria do resfriamento. Os autores são explícitos em dizer que a tendência não é estatisticamente significativa e que se confunde com a profundidade das observações. Há ainda uma possibilidade intrigante. Piyush Sharda e Mark Krumholz propuseram em 2022 que a transição entre o regime primordial, dominado pelo resfriamento por hidrogênio molecular, e o regime moderno, dominado pela poeira, ocorre em torno de um centésimo da metalicidade solar. Boötes I, com [Fe/H] = −2,4, está na ponta mais rica das anãs ultratênues conhecidas e pode simplesmente estar acima desse limiar. Galáxias ainda mais pobres poderiam ter uma IMF diferente, mas são também menores, com menos estrelas para contar.

Os autores também testaram se o resultado depende das hipóteses sobre a galáxia. Numa segunda rodada de inferência, deixaram livres a distância, entre 55 e 72,5 quiloparsecs, a extinção, a idade, entre 13 e 14 bilhões de anos, e a fração de binárias, ajustando tudo junto com a IMF. A massa característica se moveu para 0,19 ± 0,04 massa solar e a largura para 0,49 ± 0,06, deslocamentos de meio desvio padrão e pouco mais de um, respectivamente. A fração de binárias se mostrou fortemente entrelaçada com a distância, já que mais binárias ou uma galáxia mais próxima deslocam as estrelas para brilhos maiores de maneira parecida, mas essa ambiguidade quase não contamina a forma da IMF. Como a análise trabalha abaixo do turn-off, a idade também pouco interfere.

O significado do resultado ultrapassa a galáxia. Se Boötes I é o que os modelos dizem que é, um sistema com menos de cem mil massas solares em estrelas que parou de evoluir em desvios para o vermelho acima de 6, então ela é o análogo local das galáxias mais comuns do universo primitivo, aquelas que o Webb só enxerga como pontos avermelhados no limite da detecção. Medir sua IMF estrela por estrela equivale a medir a IMF daquela época. E a medida diz que o processo de formação de estrelas de baixa massa já estava calibrado no primeiro bilhão de anos, do mesmo jeito que opera hoje. Os autores sugerem que esse processo se autorregula, por mecanismos como o retorno de energia das próprias estrelas recém-nascidas, e por isso pouco se importa com a época cósmica ou com diferenças modestas de composição química.

Para quem interpreta as galáxias distantes do Webb, a notícia é tranquilizadora. Nos últimos anos, parte do debate sobre galáxias aparentemente massivas demais para sua época girou em torno da possibilidade de uma IMF diferente no universo jovem, mais pesada no topo, que faria cada galáxia brilhar mais por unidade de massa e levaria os astrônomos a superestimar suas massas estelares. O estudo de Boötes I não entra nessa discussão e não fecha essa porta, porque trata da ponta de baixa massa e de um único sistema, mas retira força da ideia de que a receita mudou de forma drástica. Ao menos em ambientes com 250 vezes menos ferro que o Sol, até onde as estrelas pequenas conseguem testemunhar, o universo primitivo e o atual contam a mesma história.

Para Boötes I, os autores consideram que as observações atuais estão perto do ótimo: nas massas muito baixas, um pequeno intervalo de massa se espalha por uma faixa enorme de brilho, de modo que mais tempo de exposição mal melhoraria o limite inferior. O passo seguinte é repetir a medida em outras anãs ultratênues, cada uma com sua própria compensação entre distância, área no céu e número de estrelas: Hércules, mais distante, exigirá exposições mais longas; Coma Berenices, mais próxima e mais espalhada, exigirá mosaicos mais largos. Com meia dúzia de galáxias medidas com a mesma precisão, será possível decidir estatisticamente entre uma IMF universal e uma IMF que varia com a metalicidade. As imagens e os catálogos do programa estão públicos no arquivo MAST, e o Starwave está disponível em código aberto.

Há algo de vertiginoso em contar dez mil estrelas que nasceram antes de existir a Via Láctea tal como a conhecemos, e descobrir que elas obedecem à mesma aritmética das estrelas que nascem hoje em Órion. A luz que o Webb coletou durante cinco horas saiu dessas anãs há 212 mil anos, mas as próprias estrelas são testemunhas de um tempo sessenta mil vezes mais remoto, e sobreviveram justamente por serem pequenas, frugais, lentas. Elas não explodem, não enriquecem o meio interestelar, não chamam atenção. Apenas ficam, queimando seu hidrogênio devagar, carregando no número exato em que existem a memória de como o universo aprendeu a fazer estrelas. Pelo que elas dizem, aprendeu cedo e nunca mais mudou de ideia.