Imagine uma estrela com a massa de cem mil sóis, brilhando com a intensidade de quase quatro bilhões de estrelas como a nossa e tão inchada que seu raio equivale a mais de três vezes a distância entre o Sol e Plutão. Agora imagine que, nos últimos séculos de sua curta existência, esse colosso começa a tossir. A cada acesso, ele lança ao espaço dezenas ou centenas de massas solares de gás, em espasmos que se repetem por décadas e depois cessam por centenas de milhares de anos. O último ataque é o mais violento: ao longo de 282 anos, 348 sóis inteiros de matéria são expelidos e formam um casulo opaco que se estende por pouco mais de um ano-luz. Então, em menos de três horas, a estrela inteira desaba sobre si mesma e vira um buraco negro. Segundo um novo estudo teórico, é essa cena que o Telescópio Espacial James Webb pode estar flagrando toda vez que registra um dos misteriosos pontos vermelhos do universo primordial.

Esses objetos ganharam um apelido quase infantil, Little Red Dots, os pequenos pontos vermelhos, e desde 2023 se tornaram um dos maiores quebra-cabeças abertos pelo Webb. Eles aparecem às centenas nos campos profundos do telescópio como fontes compactas, avermelhadas, situadas em uma época em que o universo tinha entre 600 milhões e 1,5 bilhão de anos. Seus espectros exibem um contínuo vermelho na região óptica, linhas largas de emissão do hidrogênio da série de Balmer e uma estrutura pronunciada em torno da chamada quebra de Balmer, um degrau abrupto no brilho que costuma sinalizar atmosferas estelares em uma faixa específica de temperatura. A combinação não se encaixa com facilidade em nenhuma categoria conhecida, e a abundância dessas fontes superava em muito o que os censos de quasares daquela era previam.

A interpretação inicial, de que se tratava de núcleos galácticos ativos empoeirados, tropeçou em uma série de anomalias, como a fraqueza extrema em raios X e a ausência da variabilidade típica dos quasares. A partir de 2024, um consenso começou a se formar ao redor de outro ingrediente: gás muito denso, situado extremamente próximo da fonte central, moldando tanto o contínuo quanto as linhas. Trabalhos de Kohei Inayoshi e Roberto Maiolino, de Rohan Naidu e de Anna de Graaff, entre outros, mostraram que uma camada de hidrogênio comprimido pode absorver a luz de modo a produzir a quebra de Balmer exagerada, e um artigo publicado na Nature pela equipe de Vadim Rusakov argumentou que as linhas largas podem nascer do espalhamento da luz por elétrons nesse mesmo gás, e não de velocidades orbitais gigantescas. Surgiu daí a imagem da “estrela de buraco negro”, um buraco negro embrulhado em um envelope tão espesso que o conjunto se comporta como uma atmosfera estelar.

Só que essa imagem carrega um buraco lógico. De onde vem o casulo? Quase todos os modelos simplesmente o assumem como condição inicial, uma nuvem convenientemente colocada ao redor do motor central, sem explicar como ela chegou ali, quanto tempo sobrevive e por que possui a composição química que os espectros começam a revelar. É exatamente essa lacuna que um grupo de seis pesquisadores, liderado por Devesh Nandal, do Centro de Astrofísica Harvard e Smithsonian, se propõe a fechar em uma carta científica publicada em 20 de julho no Astrophysical Journal Letters. A equipe reúne Igor Chilingarian e Abraham Loeb, do mesmo centro, Chris Nagele, da Universidade Johns Hopkins, John Chisholm, da Universidade do Texas em Austin, e Franz Bauer, da Universidade de Tarapacá, no Chile. A proposta deles é que o casulo não é uma nuvem externa, e sim o hálito final de uma estrela supermassiva.

Para entender por que astrônomos levam a sério a ideia de estrelas com cem mil massas solares, é preciso voltar ao problema mais antigo da cosmologia dos buracos negros. Quasares alimentados por buracos negros com mais de um bilhão de massas solares já existiam quando o universo tinha menos de 700 milhões de anos. Um buraco negro nasce, no cenário convencional, do cadáver de uma estrela massiva, com algumas dezenas de massas solares. Para engordar até um bilhão, ele precisaria devorar matéria no ritmo máximo permitido pela pressão da própria radiação que emite, o chamado limite de Eddington, durante praticamente toda a idade do cosmos naquela época, sem nenhuma interrupção. As contas fecham por um triz, e qualquer pausa na alimentação as inviabiliza. Daí a atração das chamadas sementes pesadas: se o buraco negro inicial já nascesse com cem mil massas solares, o problema desapareceria.

A rota mais direta para uma semente desse porte é a estrela supermassiva. A ideia é mais antiga do que parece. Em 1963, Fred Hoyle e William Fowler propuseram objetos assim para explicar a energia dos quasares recém-descobertos, e em 1964 Subrahmanyan Chandrasekhar demonstrou que estrelas sustentadas pela pressão de radiação se tornam instáveis por efeito da relatividade geral quando atingem uma compactação crítica, condenando-se ao colapso. Nos anos 2000 e 2010, trabalhos de Volker Bromm e Abraham Loeb, Mitchell Begelman, Takashi Hosokawa, Tyrone Woods e Lionel Haemmerlé mostraram que nuvens de gás primordial, em halos onde a formação de estrelas comuns é suprimida, podem despejar matéria sobre um embrião estelar com tanta rapidez que ele cresce a dezenas ou centenas de milhares de massas solares antes de morrer. Estudos mais recentes ampliaram os cenários possíveis, incluindo aglomerados estelares densos e gás levemente enriquecido em elementos pesados.

Nos últimos dois anos, a conexão entre essas estrelas e os pontos vermelhos ganhou tração. Artigos de Nandal em parceria com Loeb e de John Chisholm mostraram que o espectro de uma estrela supermassiva reproduz várias das assinaturas espectrais que definem os LRDs, e outros, como o de Lorenz Zwick, propõem estrelas supermassivas embebidas em estruturas de acreção autogravitantes. Faltava, porém, a peça que as interpretações baseadas em gás denso exigem: se o objeto central é uma estrela, quem fornece o casulo? A resposta proposta na nova carta passa por um fenômeno pouco intuitivo, a perda de massa por pulsação, e por uma classe de oscilações que os astrofísicos batizaram, com sinceridade rara, de modos estranhos.

Estrelas pulsam por diferentes motivos. As cefeidas, faróis clássicos da astronomia, oscilam graças a um mecanismo de válvula: uma camada de hélio parcialmente ionizado retém calor quando comprimida e o libera quando se expande, empurrando as camadas externas no momento certo, como o pistão de uma máquina a vapor. Outras estrelas pulsam porque a taxa de queima nuclear responde à compressão. Os strange modes, ou modos estranhos, pertencem a uma terceira família. Eles surgem em envelopes tão dominados pela pressão de radiação que o gás vira um espectador, e em que o calor atravessa as camadas externas tão depressa que a relação habitual entre compressão e pressão se embaralha. A oscilação passa a se alimentar do próprio fluxo de energia que escapa da estrela, um pouco como uma bandeira tremulando ao vento: não é o tecido que gera o movimento, e sim a corrente que o atravessa, leve demais para resistir. Desde os anos 1990, trabalhos de Kiriakidis, Glatzel e Saio mostraram que essas oscilações aparecem em estrelas muito luminosas, e uma estrela supermassiva é o caso mais extremo imaginável.

Para testar a hipótese, a equipe recorreu ao GENEC, o código de evolução estelar desenvolvido em Genebra, e construiu cinco modelos de estrelas que crescem por acreção até cerca de cem mil massas solares. Cada um começa como um embrião de dez massas solares, inteiramente convectivo, e recebe gás a uma taxa quase constante de uma massa solar por ano, ou seja, engole um Sol inteiro a cada volta da Terra em torno do nosso. Os modelos diferem na metalicidade, a fração de elementos mais pesados que o hélio, que vai de zero, correspondente às estrelas de população III formadas do gás puro do Big Bang, até um centésimo do valor solar, passando por três degraus intermediários. A rede nuclear acompanha a queima de elementos até o grupo do ferro, e cada estrutura armazenada é descrita por mais de 1.200 camadas radiais, da superfície ao centro.

A história de acreção escolhida é deliberadamente simples, e os autores reconhecem que a realidade é mais turbulenta. Em ambientes reais, o fluxo de gás sobre a estrela varia bastante, e se a taxa cai abaixo de cerca de 0,025 massa solar por ano durante a fase de pré-sequência principal, a estrela se contrai e esquenta, o que pode disparar episódios adicionais de pulsação ainda durante o crescimento. A equipe optou por isolar o que acontece depois que a acreção termina, por uma razão clara: o casulo relevante para a fase de ponto vermelho é o último ejetado antes do colapso. Cascas expelidas cedo têm centenas de milhares de anos para se expandir e se tornam transparentes muito antes de o objeto ser observável como LRD.

O tratamento das pulsações combina dois níveis de análise. Primeiro, para cada modelo armazenado ao longo da evolução, os pesquisadores resolvem o problema das oscilações radiais adiabáticas na gravidade newtoniana, obtendo os modos de menor ordem, e realizam um cálculo complementar em relatividade geral, seguindo o método publicado por Saio e colegas em 2024, para monitorar a aproximação da instabilidade relativística. Em seguida, avaliam, camada por camada, quanto trabalho a radiação realiza sobre a oscilação ao longo de um ciclo, para saber se o modo ganha energia e cresce ou se a perde por amortecimento. Quando um modo cresce, a energia disponível é convertida, por conservação de energia, em uma estimativa de perda de massa que depende da velocidade de escape efetiva das camadas externas, com um teto imposto pela quantidade de matéria situada acima da região onde a pulsação é gerada.

Os próprios autores insistem em que essa última etapa é uma estimativa estruturada, e não uma previsão única. O trabalho não resolve a hidrodinâmica não linear da radiação, e o mapeamento entre a força da pulsação e a massa ejetada depende de fatores de eficiência que eles variam entre um ramo de referência e um ramo de limite superior. É como estimar quanta água uma bomba consegue elevar a partir de sua potência e da altura do reservatório, sem simular cada gota. O procedimento pode errar no detalhe, mas captura a física essencial, e o resultado principal, como se verá, não depende dos valores exatos: depende da diferença brutal entre a última ejeção e todas as anteriores.

Os pesquisadores concentraram a análise no modelo com um centésimo da metalicidade solar, por ser o que mais se aproxima do regime de baixa, mas não nula, abundância de metais associado aos pontos vermelhos e a seus análogos próximos. Nenhuma prescrição de vento estelar convencional foi imposta, de modo que toda a perda de massa discutida é de origem pulsacional. Durante a fase de crescimento por acreção, o modelo se mostra pulsacionalmente quieto, e a ação começa depois que o fluxo de gás cessa e a estrela, com cem mil massas solares acumuladas, precisa se reorganizar.

O que acontece então lembra uma respiração profunda. Sem o gás que mantinha o envelope inflado, a estrela se contrai na escala de tempo de Kelvin-Helmholtz, o tempo necessário para irradiar a energia liberada pela própria compressão gravitacional. No diagrama de Hertzsprung-Russell, o mapa que relaciona luminosidade e temperatura, a trajetória parte de uma temperatura superficial de cerca de 8.300 kelvin e caminha para o azul até atingir quase 21.000 kelvin, enquanto a luminosidade permanece praticamente inalterada, pouco abaixo de quatro bilhões de vezes a solar. Essa contração dura cerca de 135 mil anos e transcorre sem nenhuma perda de massa por pulsação.

A contração termina quando a fração de hidrogênio no centro ainda é de 69%. A essa altura, a queima pelo ciclo CNO, em que carbono, nitrogênio e oxigênio servem de catalisadores para fundir hidrogênio em hélio, já domina o balanço energético, e o fluxo de energia que sobe do núcleo passa a realizar trabalho sobre o envelope, inflando-o novamente. A estrela retorna a uma configuração semelhante à de uma supergigante vermelha, só que agora no ritmo lento da escala nuclear, e não no ritmo rápido da escala térmica. Ao longo dos 900 mil anos seguintes, ela desliza de volta para o vermelho até atingir cerca de 7.900 kelvin, com 1,1 milhão de anos de idade e apenas 10% de hidrogênio restante no centro. Durante esse deslizamento, o envelope dominado pela radiação fica cada vez menos ligado gravitacionalmente, e suas camadas externas permanecem fortemente não adiabáticas, o terreno ideal para os modos estranhos.

Os três primeiros episódios de ejeção acontecem ao longo desse lento deslocamento para o vermelho. A etapa final se passa a temperatura quase constante, entre 1,10 e 1,27 milhão de anos, no interior do que os autores chamam de corredor da quebra de Balmer, a faixa de temperaturas em que uma atmosfera produz o degrau espectral mais pronunciado. É ali que ocorre o quarto e último episódio. Aos 1,27 milhão de anos, com apenas 1% de hidrogênio no centro, o modelo satisfaz o critério de instabilidade relativística e entra em colapso dinâmico. Um detalhe chama atenção nesse percurso: apesar de a temperatura superficial variar por um fator superior a 2,5, a luminosidade permanece confinada a uma faixa de apenas 5% ao longo de toda a evolução pós-acreção, entre 3,5 e 3,7 bilhões de luminosidades solares.

A perda de massa não se espalha suavemente por esse milhão de anos. Ela se concentra em quatro surtos bem separados, centrados em 322 mil, 520 mil, 959 mil e 1,218 milhão de anos, com intervalos de quietude de centenas de milhares de anos entre eles. E a intensidade cresce de um surto para o outro. No ramo de referência, a taxa média de perda de massa sobe de 0,25 e 0,32 massa solar por ano nos dois primeiros episódios para 0,67 e 1,23 no terceiro e no quarto; no ramo superior, os valores são 0,84, 1,08, 2,25 e 2,13. As durações também se alongam: 41,1 e 60,6 anos nos dois primeiros, 151,2 e 281,9 anos nos dois últimos. Para dimensionar esses números, a Grande Erupção de Eta Carinae, no século XIX, um dos eventos mais violentos já testemunhados em uma estrela da Via Láctea, lançou ao espaço algo entre dez e vinte massas solares em cerca de duas décadas, uma taxa comparável, mas sustentada por um tempo quatorze vezes menor.

As massas ejetadas contam a mesma história de forma ainda mais nítida. No ramo de referência, os quatro episódios expelem 10,4, 19,7, 101,9 e 348,0 massas solares; no ramo superior, 34,6, 65,6, 339,5 e 599,9. Somados, os totais chegam a 480 massas solares na estimativa de referência e a 1.040 no limite superior. Mesmo o maior desses valores corresponde a apenas 1% da massa da estrela. As pulsações, portanto, não a destroem; removem uma fração modesta, porém decisiva, do envelope externo. E o último episódio domina o orçamento, respondendo por 73% da perda total no ramo de referência e por 58% no ramo superior. Esse fortalecimento não decorre de um reservatório maior de gás disponível, já que o teto de massa acessível permanece em torno de 600 massas solares nos quatro casos. A explicação mais natural é que, à medida que a estrela volta a uma configuração fria e estendida e se aproxima da instabilidade relativística, o acoplamento entre a pulsação e o envelope se torna muito mais eficiente.

Os tempos envolvidos revelam a natureza peculiar desses eventos. Os períodos de pulsação dos modos responsáveis são de 0,411, 0,606, 1,512 e 2,819 anos, e os tempos de crescimento da instabilidade, de 0,412, 0,406, 0,688 e 2,242 anos. Cada episódio dura, assim, cerca de cem períodos de pulsação, ou entre cem e duzentos tempos de crescimento. É como se um coração gigantesco batesse uma centena de vezes, com cada batida durando de cinco meses a quase três anos, e só então se aquietasse. Ao mesmo tempo, os episódios são muito mais curtos do que o meio intervalo entre modelos evolutivos vizinhos, que vai de 60 mil a 99 mil anos. Não são impulsos de um único ciclo nem ventos quase estacionários; são episódios finitos de pulsação, com começo, meio e fim.

O que decide se uma pulsação consegue arrancar gás da estrela é a comparação entre as velocidades envolvidas. A velocidade de escape na base do envelope cai ao longo da sequência, de 2.587 e 2.283 quilômetros por segundo nos dois primeiros episódios para 1.583 e 1.176 nos dois últimos, valores que ainda superam por larga margem os 618 quilômetros por segundo necessários para escapar do Sol. A velocidade característica mais simples da pulsação, o raio dividido pelo período, é pequena demais, caindo de 307 para 216 quilômetros por segundo. Só que uma oscilação senoidal com amplitude da ordem do raio atinge, no pico, uma velocidade 2π vezes maior. Essa escala sobe para 74%, 73% e 88% da velocidade de escape nos três primeiros episódios e chega a 115% no último. No derradeiro surto, e apenas nele, o ritmo natural da pulsação ultrapassa a velocidade que a gravidade exige, como um trampolim cujo impulso, pela primeira vez, lança o saltador para além do alcance do chão.

Há ainda duas réguas radiativas que ajudam a interpretar o fenômeno. A primeira mede quanta velocidade a pressão da luz conseguiria imprimir ao gás se cada fóton empurrasse a matéria uma única vez: 298, 232, 111 e 61,5 quilômetros por segundo nos quatro episódios, sempre muito abaixo do escape. A segunda mede o limite energético total, a velocidade que o gás alcançaria se toda a luminosidade fosse convertida em movimento, e fica entre 13.100 e 5.960 quilômetros por segundo, confortavelmente acima do escape. A conclusão é reveladora. A estrela tem energia de sobra para ejetar o envelope, mas a radiação sozinha, empurrando por pressão direta, não conseguiria fazê-lo. É a pulsação em modo estranho que funciona como uma engrenagem, acoplando a energia disponível ao gás fracamente ligado das camadas externas. Uma bateria carregada de nada serve sem um fio capaz de conduzir a corrente.

O passo seguinte é perguntar como fica a vizinhança da estrela no instante do colapso, quando um observador a veria como ponto vermelho. As velocidades de lançamento variam pouco de um episódio para outro, com um intervalo inferior entre 220 e 310 quilômetros por segundo e um superior entre 1.360 e 1.930. O que muda drasticamente é o tempo de voo. As três primeiras erupções tiveram centenas de milhares de anos para se afastar, enquanto a quarta nasce pouco antes do fim. No momento do colapso, as três cascas antigas ocupam faixas gigantescas, com bordas internas a 281, 190 e 59 parsecs da estrela e bordas externas a 1.766, 1.196 e 371 parsecs, distâncias que chegam a quase seis mil anos-luz. A quarta é de outra natureza: sua borda interna coincide com a superfície da estrela, a 0,000623 parsec, ou cerca de 130 unidades astronômicas, e sua borda externa alcança apenas 0,392 parsec, pouco mais de um ano-luz. O resultado é uma estrutura em camadas: três relíquias difusas e enormes envolvendo uma única ejeção compacta, colada ao astro.

A opacidade dessas cascas separa o relevante do irrelevante. Para as três primeiras, a profundidade óptica, a medida de quanto a luz é bloqueada ao atravessá-las, é desprezível, na casa de bilionésimos a décimos de milionésimo, mesmo em suas bordas internas mais densas. A quarta casca é opaca: sua profundidade óptica geométrica é de 8,04, o valor na borda interna chega a 5.060 e mesmo na borda externa, já bem diluída, permanece em 0,0128. Os autores tratam esses números como diagnósticos de ordem de grandeza, calculados com uma opacidade constante de espalhamento por elétrons, e não como soluções completas de transferência radiativa. Ainda assim, o veredito é inequívoco: no instante do colapso, apenas o último evento produz uma casca circunstelar capaz de reprocessar a luz. Um comportamento parecido é conhecido nas variáveis do tipo R Coronae Borealis, supergigantes pobres em hidrogênio que, ao ejetar cascas de poeira, chegam a esmaecer por um fator de até cinco mil.

As massas das cascas fecham o quadro. As três primeiras carregam 10,4, 19,7 e 101,9 massas solares diluídas em volumes de milhares de anos-luz; a última concentra 348 massas solares em uma bolha de um ano-luz. A casca observacionalmente relevante, escrevem os autores, não é o resíduo acumulado de toda a perda de massa anterior, e sim o produto direto do último surto em modo estranho. Isso tem uma consequência prática para quem interpreta espectros: as propriedades do casulo de um ponto vermelho devem refletir as camadas mais externas da estrela no momento em que ela estava prestes a morrer, e não uma média de toda a sua vida.

E essas camadas têm uma assinatura química própria. Na estimativa de referência, a casca final contém 243 massas solares de hidrogênio e 105 de hélio, acompanhadas de 0,00484 massa solar de carbono, 0,013 de nitrogênio e 0,011 de oxigênio. No limite superior, os números são 419 de hidrogênio, 180 de hélio, 0,00835 de carbono, 0,0224 de nitrogênio e 0,0204 de oxigênio. Em ambos os casos, o material é dominado por hidrogênio e hélio, mas o padrão dos elementos-traço é revelador: há mais nitrogênio do que oxigênio, e mais oxigênio do que carbono. Em razões logarítmicas de número de átomos, a casca de referência apresenta log(N/O) de 0,13 e log(C/O) de −0,23, e a casca de limite superior, 0,10 e −0,26. A razão hélio/hidrogênio fica em 0,108 e 0,107. Para comparação, no Sol o nitrogênio é sete vezes menos abundante que o oxigênio; nessa casca, é 35% mais abundante, um enriquecimento relativo de cerca de dez vezes.

O excesso de nitrogênio é uma impressão digital do ciclo CNO. Nesse ciclo, a etapa mais lenta é a captura de um próton pelo nitrogênio-14, o que faz os núcleos catalisadores se acumularem nesse elemento à custa do carbono e do oxigênio. Em uma estrela supermassiva, o núcleo convectivo abrange uma fração enorme da massa total, e o material processado alcança camadas próximas da superfície, de onde a pulsação o arranca. O ponto crucial, segundo o artigo, é que a assinatura é robusta: ela sobrevive à incerteza sobre a profundidade exata que a última erupção atinge no envelope, pois as duas estimativas, referência e limite superior, entregam praticamente o mesmo padrão. Em cálculos anteriores do próprio grupo, o enriquecimento em nitrogênio já havia emergido como marca recorrente das ejeções de estrelas supermassivas.

Esse detalhe conecta o modelo a uma tendência que a espectroscopia do Webb vem revelando. Estrelas supermassivas são discutidas há anos como fontes dos padrões químicos anômalos dos aglomerados globulares, em que carbono e nitrogênio, sódio e oxigênio, magnésio e alumínio variam de estrela para estrela de modo anticorrelacionado. Em 2023, a galáxia GN-z11, observada quando o universo tinha cerca de 430 milhões de anos, mostrou-se fortemente enriquecida em nitrogênio, reacendendo o interesse por esse tipo de poluição química em ambientes de protoaglomerados. Assinaturas semelhantes começam a aparecer em pontos vermelhos: o LRD ultraluminoso A2744-45924 exibe linhas intensas de nitrogênio duas e três vezes ionizado, o objeto CANUCS-LRD-z8.6 mostra nitrogênio triplamente ionizado em um espectro de baixa metalicidade semelhante ao de um núcleo ativo, e um LRD situado em redshift 6,98 apresenta uma linha estreita de nitrogênio quatro vezes ionizado com razões elevadas em relação ao carbono e ao hélio. Essas razões de linhas no ultravioleta não se traduzem diretamente em abundâncias, pois dependem da estrutura de ionização e da transferência ressonante, mas indicam que gás denso e rico em nitrogênio é um observável importante nessa população.

O modelo oferece também uma solução elegante para uma das características mais debatidas dos pontos vermelhos: a temperatura. Vários trabalhos recentes sustentam que ao menos parte do contínuo desses objetos vem de gás denso e frio, com temperaturas de poucos milhares de kelvin, muito abaixo do que se espera da superfície de uma estrela supermassiva ou de um disco de acreção. Segundo a equipe de Nandal, a temperatura observada não precisa ser a da superfície hidrostática da estrela. Ela pode ser fixada pela última casca opaca ejetada, que funciona como uma fotosfera aparente. A lei de Stefan-Boltzmann diz que, para reemitir uma dada luminosidade a uma temperatura menor, é preciso uma superfície maior. Com cerca de 3,7 bilhões de luminosidades solares, reemitir a luz a 4.000 kelvin exige um raio fotosférico de apenas 590 unidades astronômicas, e a 3.000 kelvin, de 1.050 unidades astronômicas. Esses raios são só algumas vezes maiores que o da estrela e ficam muito aquém da borda externa da casca, situada a mais de 80 mil unidades astronômicas.

O tempo necessário para montar essa fotosfera falsa é curtíssimo. Com as velocidades de lançamento do último episódio, entre 220 e 1.360 quilômetros por segundo, a casca alcança o raio correspondente a 4.000 kelvin em dois a treze anos, e o raio de 3.000 kelvin em quatro a 23 anos. Uma fotosfera aparente fria pode, portanto, se estabelecer poucos anos ou algumas décadas após a ejeção. A massa da casca, entre 350 e 600 massas solares, é grande o bastante para que apenas uma opacidade modesta seja necessária para manter parte do material com profundidade óptica próxima de um nessas distâncias, e tanto a expansão adiabática quanto a posterior formação de poeira reforçam a tendência. Uma estrela intrinsecamente mais quente, com 7.000 a 8.000 kelvin na superfície, pode assim exibir um contínuo com temperatura aparente de apenas 3.000 a 4.000 kelvin.

Esse arranjo escapa de uma armadilha em que caem os modelos que atribuem à fonte central uma fotosfera hidrostática de fato fria. Uma superfície estelar a 3.000 ou 4.000 kelvin se pareceria com a atmosfera de uma supergigante vermelha e deveria exibir bandas moleculares intensas, sobretudo de óxido de titânio, além de bandas de carbono diatômico e cianogênio, visíveis mesmo em espectros de baixa resolução. Uma fonte mais quente vista através de uma casca densa em expansão, ao contrário, produz um contínuo muito mais frio sem exigir que o objeto subjacente seja uma estrela fria, e um pouco de avermelhamento torna uma atmosfera de 7.000 kelvin parecida com uma mais fria, porém sem as bandas moleculares. O análogo de baixo redshift J1025+14 dá suporte a essa cautela: seu contínuo parece muito frio, mas suas linhas de absorção indicam temperaturas superficiais substancialmente mais altas. Estimativas de temperatura baseadas apenas no contínuo, alertam os autores, podem enviesar as fontes para valores artificialmente baixos.

O cenário rende uma previsão estatística que os levantamentos futuros poderão confrontar. Se uma estrela supermassiva só é classificada como ponto vermelho enquanto o gás denso lhe empresta uma fotosfera fria ou reprocessa uma fração substancial do contínuo, então a razão entre o número de estrelas envoltas e o número de estrelas nuas é dada, grosso modo, pela fração do tempo de vida em que a casca permanece opticamente relevante, multiplicada pela fração do céu ao redor da estrela que ela cobre. Com uma vida de 1,27 milhão de anos, uma ejeção final de 280 anos e uma possível sobrevida do casulo de até dez mil anos, essa razão fica entre dois décimos de milésimo e oito milésimos, corrigida pelo fator de cobertura. Em outras palavras, estrelas supermassivas pouco ou nada obscurecidas podem ser de mais de cem a alguns milhares de vezes mais comuns do que os pontos vermelhos, mas não seriam reconhecidas como tais, pois lhes faltaria o contínuo vermelho reprocessado. Ventos, rotação, assimetrias da casca e o gás natal podem alterar tanto o tempo de sobrevivência quanto a cobertura, e uma previsão precisa exigirá simulações radiativo-hidrodinâmicas da casca e uma modelagem cuidadosa da seleção observacional.

A sobrevida do casulo pode ser prolongada por um mecanismo adicional. Se parte do material ejetado permanecer marginalmente ligado à estrela, ou for freado pelo reservatório denso de gás natal e estacionar na escala de 0,392 parsec, ele retornará à estrela em cerca de 13 mil anos. Ejeções maiores e raios de estagnação mais amplos alongam esse tempo de retorno. O efeito é duplo: manter o casulo opaco por mais tempo e oferecer um novo reservatório de acreção. Se o gás que retorna guardar momento angular suficiente, parte dele pode se assentar em um disco, o que abriria um elo entre o colapso da estrela supermassiva e uma fonte de acreção envolta em gás e mais longeva, na linha dos cenários de quase-estrela propostos por Begelman em 2008 e das estrelas de buraco negro discutidas por Naidu e colegas. Os autores também ressalvam que não incluíram todas as formas de perda de massa: ventos impulsionados por linhas espectrais e pulsações de origem termonuclear, estudadas por Chris Nagele, podem liberar entre cem e mil massas solares adicionais. Quanto mais material é ejetado, mais leve o buraco negro remanescente, mas mais tempo o casulo permanece opaco e parecido com um LRD.

Enquanto o casulo se expande, a estrela caminha para o fim. Aqui entra a instabilidade prevista por Chandrasekhar. Em uma estrela sustentada quase inteiramente pela pressão de radiação, o equilíbrio é precário: a resposta da pressão à compressão fica no limiar exato da estabilidade. A relatividade geral acrescenta a esse balanço um pequeno excesso de atração, que cresce com a compactação do objeto, e para uma estrela de cem mil massas solares esse excesso basta para romper o equilíbrio quando o núcleo se contrai perto do fim da queima de hidrogênio. A equipe calculou o início da instabilidade por três critérios independentes. A análise linear adiabática relativística de Saio e a estimativa pós-newtoniana de Haemmerlé situam o ponto crítico em cerca de 1,3 milhão de anos, enquanto o método relativístico de Nagele, que contabiliza perturbações não lineares, encontra um início um pouco anterior, em 900 mil anos. Todos concordam que a estrela vive aproximadamente um milhão de anos e forma um buraco negro após perder apenas uma pequena fração de seu envelope.

Para confirmar o desfecho, os pesquisadores transferiram a estrutura instável para um código de hidrodinâmica relativística unidimensional, equipado com uma rede nuclear de 52 isótopos e resfriamento por neutrinos. O resultado é um colapso, não uma explosão. O buraco negro se forma dentro de cerca de dez mil segundos após o início da instabilidade, menos de três horas, e a vasta maioria da matéria estelar cai nele. Nas primeiras fases, os perfis de velocidade são quase homólogos, com toda a estrela caindo em conjunto; depois, as regiões internas aceleram mais que as externas e o colapso se concentra progressivamente no centro, até a formação do horizonte. A semente resultante conserva cerca de 99% da massa da estrela no momento da instabilidade, ou seja, nasce com aproximadamente cem mil massas solares, exatamente o tipo de objeto que os cenários de colapso direto invocam para explicar os quasares gigantes do universo jovem.

Os modelos apresentados não giram, e os autores tratam a rotação como a principal incerteza pendente. Ainda não existe um cálculo GENEC de pós-acreção com rotação para a faixa de cem mil massas solares. Os modelos rotantes disponíveis, restritos à sequência principal e a massas entre mil e dez mil sóis, mostram que a rotação estende a vida da estrela em 7%, 14% e 20% para mil, cinco mil e nove mil massas solares, e que esses objetos esbarram no limite combinado de rotação e luminosidade mesmo girando devagar, com velocidades entre 5% e 10% da crítica. Para estrelas em crescimento, trabalhos de Haemmerlé indicam que atingir a escala supermassiva exige acretar apenas cerca de 1% do momento angular kepleriano, o que resulta em rotação superficial lenta e deformação centrífuga desprezível, embora frações pequenas de momento angular possam elevar a massa máxima estável em uma ordem de grandeza. A expectativa da equipe é que a rotação altere o momento do colapso e talvez a distribuição angular das ejeções, mas não elimine o canal de perda de massa, cuja origem está na razão elevada entre luminosidade e massa, na baixa fração de pressão do gás e nos tempos térmicos curtos do envelope.

Outra escolha deliberada foi não impor um vento estelar contínuo nos modelos de referência, e a justificativa é honesta: nenhuma prescrição de vento foi calibrada para envelopes de cem mil massas solares no limiar de Eddington. A fórmula de Vink, amplamente usada, foi desenvolvida para ventos impulsionados por linhas em estrelas quentes dos tipos O e B; a relação de de Jager é um ajuste empírico a estrelas muito menos massivas. Aplicar qualquer uma delas exige extrapolar por ordens de grandeza em massa, luminosidade e fator de Eddington. Um teste anterior, conduzido por Nandal e Chon em uma estrela com um décimo de milésimo da metalicidade solar, ilustra a incerteza: partindo de 73.973 massas solares no meio da queima de hidrogênio, a prescrição de de Jager reduziu a massa final a 23.421 massas solares, enquanto a de Vink a levou a 62.348, e o vento à la Vink só se mostrou relevante durante breves fases quentes, acima de 25 mil kelvin. Se ventos reais forem tão fortes quanto as extrapolações mais agressivas sugerem, eles podem precondicionar o envelope e reduzir o reservatório disponível para as ejeções pulsacionais. Os resultados da carta devem ser lidos, por isso, como a componente estritamente pulsacional da perda de massa.

Embora o modelo de um centésimo da metalicidade solar seja o protagonista, os outros quatro mostram que o fenômeno não é uma peculiaridade dele. Todas as cinco sequências terminam com massas finais quase idênticas, entre 100.185 e 100.888 massas solares, e todas produzem ejeções pulsacionais de cascas, da população III ao modelo mais enriquecido. O que muda com a metalicidade é o modo como isso acontece. O modelo de um milésimo da metalicidade solar produz o maior número de episódios, 28, enquanto o de um centésimo de milésimo ejeta a maior massa acumulada, 6.262 massas solares na estimativa de referência e 14.656 no limite superior. O modelo de um centésimo, em contraste, exibe apenas quatro episódios e a menor massa total ejetada. Frequência de erupções e eficiência na remoção de massa, portanto, não são a mesma coisa. A dependência não monotônica é esperada, pois a instabilidade em modo estranho tem duas famílias, uma ligada à ionização do hélio, que atua mesmo com metalicidade quase nula, e outra ligada à opacidade dos elementos pesados, que se fortalece à medida que eles se tornam mais abundantes.

O estado final das estrelas divide a grade em dois regimes. O modelo de população III é o ponto fora da curva: termina muito mais quente e compacto, com temperatura superficial próxima de 74 mil kelvin e velocidade de escape de 7.731 quilômetros por segundo. Os modelos enriquecidos se agrupam em torno de 10 mil kelvin, com velocidades de escape entre 1.200 e 1.600 quilômetros por segundo, e todos mantêm o fator de Eddington próximo de um: a pressão da radiação quase equilibra a gravidade na superfície. Em toda a grade, a casca mais relevante se mantém compacta e opaca, com profundidades ópticas de espalhamento por elétrons entre 740 e 13 mil. A química também varia: a casca de população III é extrema, com razões carbono/oxigênio e nitrogênio/oxigênio que superam as solares por fatores da ordem de mil e de dezenas de milhares, enquanto o modelo de um milésimo entrega cascas robustamente ricas em nitrogênio, com log(N/O) entre 0,334 e 0,359. A razão hélio/hidrogênio varia muito menos que as razões de carbono, nitrogênio e oxigênio, confirmando que as ejeções expulsam o envelope externo, e não um núcleo inteiramente processado. O modelo de um centésimo se destaca não por ejetar mais massa, e sim por combinar o estado terminal enriquecido mais frio, uma casca compacta e opaca e uma composição rica em nitrogênio, o análogo mais claro de um LRD.

O artigo encerra com um catálogo de previsões observáveis. Uma estrela supermassiva com temperatura superficial de cerca de 7.000 kelvin, envolta em cascas empoeiradas, deveria se parecer com uma hipergigante amarela, classe rara de estrelas da Via Láctea e de galáxias vizinhas, só que ampliada em luminosidade e com outra química. Objetos como IRC+10420, uma das estrelas mais luminosas da Galáxia, e IRAS 17163-3907, a chamada Nebulosa do Ovo Frito, exibem espectros moderadamente avermelhados de tipos A tardio e F inicial, com linhas largas de emissão e absorção formadas em envelopes circunstelares complexos. Uma estrela supermassiva deveria gerar assinaturas semelhantes, com linhas de oxigênio, carbono e magnésio mais fracas e as de nitrogênio mais fortes; eventuais linhas de ferro refletiriam a composição do gás de nascimento, não produção durante a queima de hidrogênio. As múltiplas cascas, lançadas a velocidades entre cem e mil quilômetros por segundo, deveriam imprimir linhas de absorção deslocadas para o azul, com vários componentes ou perfis assimétricos, em espécies como sódio e potássio neutros. Ventos impulsionados por linhas, por sua vez, produziriam perfis do tipo P Cygni, a combinação de emissão com uma absorção deslocada para o azul que denuncia gás em expansão, em linhas como o H-alfa e o hélio neutro metaestável. Traços largos de hidrogênio e hélio desse tipo já foram detectados em pontos vermelhos tanto em baixo quanto em alto redshift.

O que torna a proposta atraente é a forma como ela costura ideias que costumam ser discutidas separadamente. No quadro do colapso direto, sementes massivas se formam do gás protogaláctico em colapso. Nos modelos de quase-estrela, um buraco negro cresce dentro de um envelope hidrostático maciço. Nas interpretações recentes dos LRDs como estrelas de buraco negro, um buraco negro está embebido em gás denso que molda o espectro emergente. A evolução de uma estrela supermassiva, sugerem os autores, pode ser a ponte entre esses cenários: ela primeiro monta a casca compacta exigida pela fase observável de ponto vermelho e depois colapsa quase intacta na própria semente pesada. Estrutura circunstelar, aparência de LRD e nascimento do buraco negro deixam de ser três problemas para se tornar três capítulos de uma única história. O próximo passo, indicam eles, é converter essas cascas em modelos diretos de contínuo e de linhas espectrais, para que a ejeção pulsacional possa ser testada frente aos espectros reais que o Webb continua a acumular.

Resta, claro, a dúvida de sempre na astrofísica teórica: a natureza terá seguido o roteiro? As estimativas de massa ejetada são aproximadas, a opacidade foi tratada como constante, a rotação ficou de fora, e os ventos permanecem um coringa. Mas o que a carta entrega é algo que a hipótese das estrelas supermassivas ainda não tinha, um mecanismo físico concreto, com números, para o casulo que os espectros exigem, e uma lista de assinaturas que podem derrubá-lo. Se ele estiver certo, cada ponto vermelho nos campos profundos do Webb é uma estrela que nunca deveria existir pelos padrões do universo atual, vista em seu último suspiro, vestida com o gás que ela mesma exalou nos derradeiros três séculos de vida, e prestes a se converter, em uma tarde, no buraco negro que um dia iluminará um quasar. O telescópio estaria fotografando, em um único pixel avermelhado, o atestado de óbito de uma estrela e a certidão de nascimento de um gigante.