Na constelação de Órion, a pouco mais de 1.300 anos-luz da Terra, existe uma nuvem escura que os astrônomos chamam pelo nome burocrático de Orion B. Dentro dela estão alguns dos cartões-postais mais reproduzidos da história da astronomia: a silhueta da Nebulosa Cabeça de Cavalo, o brilho da Nebulosa da Chama, a poeira iluminada de NGC 2023. O que interessa a um grupo de 29 pesquisadores liderados por Helena J. Mazurek, do Observatório de Paris, não são os cartões-postais. São os casulos de gás frio e denso escondidos atrás deles, mais de mil, cada um guardando a possibilidade de virar uma estrela. A pergunta que o grupo levou a um radiotelescópio de 30 metros no sul da Espanha é desconcertante de tão simples: esses casulos são todos iguais, como a astronomia supôs por décadas, ou cada um carrega a marca química do bairro onde nasceu?

Para a ciência, “casulo” tem nome técnico: núcleo pré-estelar. É a fase mais precoce da formação de uma estrela, quando uma bolha de gás molecular com algumas dezenas de milhares de unidades astronômicas de diâmetro, fria a ponto de beirar os 10 kelvins, começa a perder a batalha contra a própria gravidade. Durante muito tempo, o que se sabia sobre esses objetos veio de um punhado de exemplares famosos em nuvens próximas, dissecados com detalhe quase obsessivo: L1544, na nuvem de Touro; Barnard 1b, em Perseu; L1689N e IRAS 16293E, em Ofiúco; e o núcleo TMC-1, rico em moléculas de carbono. Esses objetos viraram gabaritos. Sempre que um núcleo novo era encontrado em algum canto da Galáxia, o reflexo era compará-lo a L1544 e supor que o resto seguiria o mesmo manual. O estudo publicado na revista Astronomy & Astrophysics, volume 713, artigo A137, recebido em 15 de abril de 2026 e aceito em 22 de junho, nasceu de uma desconfiança contra esse reflexo: e se o gabarito não representar a população?

A desconfiança tem base numérica. A maior parte do que se conhece sobre núcleos densos vem de mapas do contínuo de poeira, a emissão térmica dos grãos sólidos misturados ao gás, obtidos sobretudo pelo telescópio espacial Herschel, da Agência Espacial Europeia, no levantamento conhecido como Herschel Gould Belt Survey, e, mais recentemente, por programas do ALMA em nuvens distantes, como o ALMA-IMF e o ALMAGAL. Esses levantamentos produziram catálogos com milhares de entradas: posições, tamanhos, massas, estatísticas de distribuição. Um núcleo, nesse tipo de dado, é um máximo local no mapa de densidade de coluna, extraído por algoritmos especializados. Mas a poeira conta metade da história. Ela diz onde está a matéria e quanta matéria existe, e cala sobre o que o gás está fazendo e de que moléculas é feito. A química e a dinâmica, que decidem se um núcleo vai ou não colapsar, só aparecem nas linhas espectrais.

E foi atrás dessas linhas que o projeto ORION-B passou mais de uma década. A sigla vem de Outstanding Radio-Imaging of OrioN-B, e o programa, coordenado por Jérôme Pety, do Instituto de Radioastronomia Milimétrica (IRAM), e por Maryvonne Gerin, do Observatório de Paris, usou o radiotelescópio de 30 metros do IRAM, instalado no Pico Veleta, na Serra Nevada espanhola, para mapear cinco graus quadrados do céu. É uma região de 18 por 13 parsecs, algo como 59 por 42 anos-luz, na porção sul da nuvem molecular gigante, a uma distância adotada de 410 parsecs. A resolução angular típica das observações, 26 segundos de arco, corresponde a 0,05 parsec, pouco mais de dez mil vezes a distância entre a Terra e o Sol. A cobertura em frequência vai de 72 a 116 gigahertz, uma banda de 40 gigahertz em torno do comprimento de onda de três milímetros, com resolução espectral de 200 quilohertz, o equivalente a 0,5 quilômetro por segundo em velocidade.

Cada molécula do gás interestelar gira em frequências exatas, ditadas pela mecânica quântica, e ao saltar de um estado de rotação para outro emite um fóton de frequência precisa. É uma impressão digital. Na banda de três milímetros caem as transições fundamentais de dezenas de espécies, e o ORION-B selecionou 25 linhas detectadas em pelo menos uma dúzia de núcleos: as transições 1–0 de quatro isotopólogos do monóxido de carbono, versões da molécula em que um átomo é trocado por um isótopo mais pesado, ¹²CO, ¹³CO, C¹⁸O e C¹⁷O; os íons HCO⁺, H¹³CO⁺, DCO⁺ e N₂H⁺; espécies com carbono, nitrogênio e oxigênio como CCH, CN, HCN, HNC, H¹³CN, HN¹³C, DNC, HCO, H₂CO, o metanol CH₃OH, HNCO e o anel c-C₃H₂; quatro portadoras de enxofre, CS, SO, CCS e HCS⁺; e o monóxido de silício, SiO, que funciona como sensor de choques. Os cubos de DCO⁺, DNC, HNCO, HN¹³C, HCS⁺ e CCS nunca tinham sido publicados. Para a análise, todos foram reamostrados em pixels de 9 segundos de arco e suavizados para uma resolução comum de 36 segundos de arco, a do cubo de pior qualidade angular, centrado na linha do DCO⁺.

A lista de alvos veio de outro lugar. Em 2020, Vera Könyves e colaboradores publicaram o catálogo de núcleos de Orion B do Herschel Gould Belt Survey, obtido com o algoritmo getsources, um programa que decompõe as imagens em estruturas de escalas progressivamente maiores. Combinando a emissão a 160, 250, 350 e 500 micrômetros com um mapa de densidade de coluna de alta resolução, o catálogo identificou 1.768 núcleos sem estrela e 76 núcleos protoestelares, estes delatados pela emissão a 70 micrômetros de uma fonte interna que aquece a poeira ao redor. Dentro do campo coberto pelo ORION-B sobraram 1.001 objetos: 283 núcleos pré-estelares classificados como robustos, 198 candidatos a pré-estelares, 500 núcleos sem estrela e 20 objetos protoestelares. De cada um, os autores tomaram a posição, o raio observado e a classificação, sem corrigir nada, para que a amostra fosse exatamente aquela que o contínuo de poeira entrega.

A palavra “robusto” tem aqui um sentido preciso, e ele será posto à prova mais adiante. Um núcleo de gás pode ser imaginado como um balão em equilíbrio: a pressão térmica empurra para fora, a gravidade puxa para dentro. Na década de 1950, Ebert e Bonnor calcularam a massa máxima que uma esfera de gás isotérmica sustenta antes de desabar, e essa massa crítica de Bonnor-Ebert depende do raio, da temperatura e da pressão externa. O catálogo de 2020 chamou de robustos os núcleos cuja massa supera em pelo menos duas vezes a massa crítica de uma esfera do mesmo raio a 10 kelvins, ou 263 graus Celsius negativos. São, em tese, os que mais provavelmente vão formar estrelas. O detalhe que importa é que esse critério foi calibrado em nuvens frias e tranquilas como Touro e Perseu, e assume uma temperatura única, sem turbulência.

Antes de qualquer estatística, o artigo mostra quatro exemplos para educar o olhar. São mapas de momento zero, ou seja, mapas da intensidade de cada linha integrada ao longo da velocidade, para dez transições de alto sinal em torno de um núcleo sem estrela, o de número 284, de um pré-estelar, o 670, de um robusto, o 606, e de um protoestelar, o 991. A distribuição espacial da emissão muda de linha para linha e de objeto para objeto. O ¹³CO ocupa tudo, difuso, enquanto o N₂H⁺ praticamente não aparece no 284 e acende forte e concentrado no 606 e no 991. Essa diferença não é ruído: cada molécula se forma e se excita em condições específicas de densidade, temperatura e radiação, e por isso cada uma enxerga uma camada distinta da mesma estrutura. É esse mosaico de sensibilidades que o estudo pretende transformar em informação.

Extrair um número por linha e por núcleo exigiu cuidado. Para cada objeto, o espectro foi tomado como a média dentro do raio catalogado e integrado numa janela de velocidade definida pela emissão de HCO⁺, uma linha brilhante e onipresente, usada como régua comum para todas as outras. Quando o HCO⁺ não bastava, valia uma janela padrão de −4 a 16 quilômetros por segundo, estimada a partir do espectro médio de C¹⁸O de todo o campo, já que essa linha é a que melhor acompanha os filamentos onde os núcleos vivem. As linhas com estrutura hiperfina, cujas componentes se espalham em frequência, ganharam janelas próprias. Ao inspecionar os perfis de H¹³CO⁺, que não tem essa estrutura, a equipe percebeu que seis objetos da região B9 somavam duas camadas de gás em velocidades diferentes e os desdobrou em pares, o que elevou a lista a 1.007 medidas. Quatro núcleos sem dados em algumas linhas, os de número 1, 2, 3 e 982, foram descartados da etapa estatística, que ficou com 1.003 objetos.

A cada um deles o estudo associou três grandezas físicas tiradas de mapas independentes. A densidade de coluna de hidrogênio molecular e a temperatura da poeira vêm do trabalho de Marco Lombardi e colaboradores, de 2014, que ajustaram a distribuição espectral de energia do Herschel e do satélite Planck com um corpo negro modificado, converteram a opacidade a 850 micrômetros em extinção visual pela relação AV = 2,7 × 10⁴ τ₈₅₀ e esta em densidade de coluna, com 1,8 × 10²¹ átomos de hidrogênio por centímetro quadrado para cada magnitude de extinção, admitindo que todo o gás é molecular. Em torno de NGC 2024, onde os pixels do Herschel saturaram, foi aplicada uma correção adicional. A intensidade do campo ultravioleta distante, expressa pelo fator G₀ em unidades do campo médio interestelar na vizinhança solar, veio da luminosidade no infravermelho distante calculada por Miriam Santa-Maria e colegas em 2023. E a densidade volumétrica média do gás, pesada pela massa ao longo de cada linha de visada, foi tomada do mapa de Jan Orkisz e Jouni Kainulainen, de 2025, para que cada núcleo tivesse um único valor representativo.

O retrato geral das intensidades em função da densidade de coluna já diz muito. Os objetos cobrem um intervalo de 2,8 × 10²¹ a 2,4 × 10²³ moléculas de hidrogênio por centímetro quadrado, o que corresponde a extinções visuais de aproximadamente 3 a 240 magnitudes, de um véu fino a uma muralha que nenhum fóton visível atravessa. Quase todas as linhas ficam mais brilhantes onde há mais matéria, como esperado, mas a forma dessa relação varia de espécie para espécie. As linhas de ¹²CO, HCO, HCO⁺, HCN, CCH e CN exibem um gradiente nítido com a temperatura da poeira. Já N₂H⁺, DCO⁺, HN¹³C e DNC permanecem mudas até a coluna ultrapassar um limiar de cerca de 10²² por centímetro quadrado, dez sextilhões de moléculas em cada centímetro quadrado de céu, e só então disparam, sinal de uma mudança abrupta de abundância ou de excitação nesse patamar. As quatro linhas do CO saturam nas colunas mais altas, enquanto N₂H⁺, H₂CO e metanol crescem de forma quase linear, e HNC e SO ficam a meio caminho, com dispersão maior em torno da tendência.

Há personagens que fogem do roteiro. Três núcleos da região B9, os de número 301, 786 e 993, aparecem com razões de N₂H⁺ por densidade de coluna acima do normal, mas por um motivo mecânico: ali as camadas de velocidade se sobrepõem e a coluna de cada uma não pode ser separada. Já o protoestelar 996 e o núcleo 829, bem no centro do mapa, mostram o mesmo excesso sem camadas sobrepostas, e podem ter uma superabundância genuína da molécula. Espécies raras como CCS, HCS⁺ e HNCO só surgem em fontes submetidas a condições peculiares, o que impede qualquer estatística sobre elas. E o SiO, que se liberta dos grãos de poeira quando ondas de choque os despedaçam, revelou dois objetos com asas de alta velocidade típicas de jatos moleculares: o núcleo 606, perto de NGC 2023, e a protoestrela 995, a nordeste de B9.

Entre os isotopólogos do CO, dois se comportam como gêmeos. As intensidades de C¹⁷O seguem as de C¹⁸O com uma razão quase constante, de 0,291 ± 0,001 no corpo da análise e 0,292 ± 0,001 no apêndice que a detalha, com coeficiente de correlação de 0,93. A constância tem consequência física: se o C¹⁸O fosse opaco, a saturação distorceria a razão, e se as duas moléculas tivessem temperaturas de excitação distintas, a razão flutuaria. Como nada disso acontece, as duas linhas são opticamente finas e traçam o mesmo gás. Por isso o C¹⁸O, detectado na maioria dos núcleos com boa relação sinal-ruído e bem correlacionado com a coluna de poeira, foi escolhido como medidor de movimento. Seu perfil espectral reflete apenas o gás em deslocamento ao longo da linha de visada, sem necessidade de corrigir o alargamento por opacidade.

Para isso a equipe usou um segundo conjunto de dados de C¹⁸O, obtido com o correlacionador VESPA do mesmo telescópio, com resolução espectral de 0,1 quilômetro por segundo, cinco vezes melhor, e resolução angular de 23,6 segundos de arco. A cada espectro médio foi ajustada uma única gaussiana, inicializada pelo pico mais brilhante dentro da janela, com a rotina de ajuste do pacote spectral-cube e o detector de picos da biblioteca SciPy. Os ajustes ruins foram descartados por critérios explícitos: largura a meia altura acima de 5 quilômetros por segundo, que denuncia um ajuste à ondulação da linha de base; incerteza na dispersão de velocidade maior que 0,2 quilômetro por segundo, o dobro da largura do canal; e incerteza relativa acima de 30 por cento. Dos 1.007 objetos, cerca de 831 passaram no crivo. Alguns mostram duas ou mais componentes estreitas sobrepostas, que não indicam necessariamente objetos distintos, já que outras linhas, como o H¹³CO⁺, revelam uma feição só; esses casos ficaram na amostra, com a largura ampliada refletindo a estrutura média de velocidade.

O resultado vai contra a intuição herdada das nuvens tranquilas. As larguras a meia altura variam de 0,3 a 3 quilômetros por segundo, com média ponderada de 1,08 quilômetro por segundo, quase 3.900 quilômetros por hora. Todas superam a largura térmica do C¹⁸O para as temperaturas cinéticas conhecidas em Orion B, entre 10 e 60 kelvins. Mais revelador: a média é cerca do dobro da largura térmica do próprio hidrogênio molecular a 20 kelvins, de 0,5 quilômetro por segundo. Gás que se move além da sua agitação térmica está sendo empurrado por turbulência ou por movimentos sistemáticos, transônicos ou supersônicos. Só a ponta inferior da distribuição chega perto do térmico. Em Touro, a 10 kelvins, isso seria anomalia; em Orion B, uma região massiva que sofre radiação ultravioleta intensa e a retroalimentação dinâmica das estrelas que já formou, é o esperado. A nuvem tem ao menos três camadas de velocidade principais, cada uma com estruturas internas, segundo o mapeamento de Mathilde Gaudel e colaboradores em 2023. Um único perfil de linha fina não separa colapso de turbulência, porque ambos alargam a linha do mesmo jeito; separá-los exige comparar linhas finas e espessas e seguir o perfil através do núcleo.

Com as dispersões em mãos, os autores calcularam a massa virial de cada núcleo, a massa que o gás precisaria ter para que a gravidade equilibrasse a agitação interna, seguindo a prescrição de Bertoldi e McKee, de 1992. Para a maioria dos núcleos, a massa virial supera a massa de gás deduzida da poeira, como eles já tinham notado. Ela supera também as massas críticas de Bonnor-Ebert calculadas no catálogo de 2020, que assumiram ausência de turbulência e 10 kelvins. Os núcleos robustos têm massas de gás maiores que os sem estrela e os candidatos, e parecem mais próximos do equilíbrio virial, mas, no geral, as dispersões de velocidade e as massas viriais das três categorias não diferem de modo marcante. A conclusão incomoda quem usa a robustez como atalho: numa nuvem que abriga estrelas massivas, a temperatura cinética passa de 10 kelvins nas escalas observadas e o gás é transônico, e o critério de Bonnor-Ebert com temperatura única não basta para dizer quem vai colapsar. É preciso medir a temperatura ou a dispersão, e levar em conta a pressão externa e os movimentos internos, como argumentam Moon e Ostriker em trabalho de 2025.

A parte central do estudo, porém, não está na dinâmica, e sim na química, e exigiu uma ferramenta estatística para dar conta de 25 números por objeto. Com 25 intensidades, cada núcleo é um ponto num espaço de 25 dimensões, impossível de visualizar. A análise de componentes principais, ou PCA, da sigla em inglês principal component analysis, formalizada por Harold Hotelling em 1933, encontra as direções desse espaço ao longo das quais os pontos mais se espalham. Numa multidão, a primeira coisa que distingue as pessoas pode ser a altura, a segunda, a idade, a terceira, algo mais sutil. Cada eixo é uma combinação ponderada das 25 linhas, os pesos se chamam cargas, e a projeção de cada núcleo sobre o eixo é o seu escore. Os primeiros eixos concentram a maior parte da variação e, se tiverem leitura física, revelam os fatores que de fato governam a emissão. Pierre Gratier e colegas já tinham aplicado a técnica a Orion B em 2017, mas em todas as escalas da nuvem e numa área menor; agora ela foi apontada apenas para os núcleos.

Como as intensidades variam por ordens de grandeza e se distribuem de maneira muito diferente entre as linhas, com o ¹²CO espalhado até 200 kelvins vezes quilômetro por segundo e o N₂H⁺ ou o DNC amontoados perto de zero, os dados precisaram de preparo. A equipe aplicou a função seno hiperbólico inverso, que deixa os valores pequenos quase intocados e comprime os grandes como um logaritmo, e depois padronizou cada linha, subtraindo a média e dividindo pelo desvio padrão. A matriz final, de 1.003 objetos por 25 linhas, entrou na implementação da biblioteca scikit-learn, que opera sobre a matriz de correlação dos dados e por isso não deixa que uma linha brilhante domine as fracas só por ser brilhante.

Antes de interpretar, era preciso saber se o resultado dependia de uns poucos objetos esquisitos. Para isso o grupo fez uma análise leave-one-out, deixar um de fora: repetiu a PCA 1.003 vezes, cada vez retirando um núcleo, projetando-o no espaço definido pelos demais e medindo quão mal ele era reconstruído. Pelo critério clássico de John Tukey, que marca como discrepante tudo o que fica acima do terceiro quartil mais uma vez e meia o intervalo interquartil, 53 objetos se destacam. Dois deles passam de seis desvios padrão. O primeiro é o núcleo 606, aquele do jato de SiO, que também brilha em metanol de forma extensa e não tem par em todo o campo. O segundo é o 652, no coração de NGC 2024, a região de maior densidade de coluna, de maior densidade de volume e de formação de estrelas massivas de toda a nuvem. Outros discrepantes correspondem a fluxos negativos, provável sinal de dados de pior qualidade. O achado tranquilizador é que as três primeiras componentes quase não mudam com a retirada desses objetos; a partir da quarta, as cargas passam a oscilar, e por isso a interpretação ficou restrita às três primeiras.

Juntas, elas explicam em torno de 65 por cento da variância da amostra: 47,03 por cento para a primeira, 10,87 por cento para a segunda e 6,07 por cento para a terceira, com incertezas de centésimos. As outras 22 componentes dividem o restante em fatias cada vez menores, numa curva que só chega a 100 por cento quando todas são somadas. Três eixos, portanto, resumem dois terços de tudo o que distingue um núcleo de outro no espaço das 25 linhas, e cada um dos três tem uma leitura física clara, o que é raro em análises desse tipo, em que as componentes costumam misturar efeitos sem nome.

O primeiro eixo é o mais previsível. Todas as suas cargas são positivas, isto é, ele cresce quando todas as linhas crescem juntas, e seu escore acompanha a densidade de coluna de hidrogênio com coeficiente de correlação de 0,820. É a medida de quanta matéria existe ao longo da linha de visada. Trabalhos anteriores de Gratier e colaboradores, de 2017 e 2021, já tinham mostrado que essa componente fornece uma estimativa independente da coluna de gás, uma média ponderada das intensidades, e que a combinação de ¹²CO, ¹³CO, C¹⁸O, HCO⁺ e N₂H⁺, tratada por um método de floresta aleatória, recupera a coluna com precisão de regiões difusas a núcleos densos. Lucas Einig e colegas confirmaram, em 2024, com modelos teóricos e um critério de informação mútua, que a extinção total está de fato fortemente ligada à intensidade de linhas selecionadas, e que a seleção ideal depende do regime de extinção que se quer sondar.

O segundo eixo é o que dá título ao artigo. Suas cargas se dividem em dois campos de sinais opostos: positivas para ¹²CO, CN, HCO, CCH, HCO⁺, HCN, ¹³CO e c-C₃H₂; negativas para metanol, DCO⁺, N₂H⁺, DNC, H¹³CO⁺, HN¹³C, C¹⁷O e H¹³CN. Quem conhece astroquímica reconhece os dois lados. As moléculas com deutério, como DCO⁺ e DNC, só florescem em gás frio e blindado, de baixa ionização, onde o fracionamento isotópico empurra o deutério para dentro das moléculas: a 10 kelvins, as reações favorecem o isótopo pesado a tal ponto que a razão de deutério em certas espécies sobe várias ordens de grandeza acima da abundância cósmica, de cerca de uma parte em cem mil, como mostraram os modelos de Javier Goicoechea, em 2009, e de Evelyne Roueff, em 2015. Já os radicais CN e CCH nascem com eficiência nas camadas externas das regiões de fotodissociação, as superfícies das nuvens banhadas por ultravioleta, onde a ionização é maior. O segundo eixo, portanto, separa os núcleos densos, frios e protegidos dos núcleos expostos à radiação.

A grandeza física que o traduz é uma razão simples: G₀/n, a intensidade do campo ultravioleta dividida pela densidade do gás. Ela mede quantos fótons energéticos chegam para cada molécula capaz de absorvê-los. Valores baixos, abaixo de 0,01, indicam regiões densas e blindadas, de poeira fria; valores altos, acima de 0,1, indicam gás castigado pela radiação, seja porque o campo é mais forte, seja porque a densidade e a coluna são baixas o bastante para deixar os fótons penetrarem. O escore da segunda componente acompanha G₀/n com coeficiente de Pearson de 0,797, um valor alto para uma amostra de mil objetos espalhados por uma nuvem inteira. No mapa, os núcleos de escore positivo se concentram em torno de NGC 2024 e ao longo da crista ocidental da nuvem, a borda iluminada por σ Orionis, o sistema de estrelas quentes que recorta a Cabeça de Cavalo contra o fundo brilhante.

É aqui que a escolha de uma amostra selecionada pela poeira, e não pela química, mostra seu valor, porque não existe apenas um tipo de núcleo em Orion B. Existe a população clássica, mais numerosa, de núcleos frios e densos, brilhantes em H¹³CO⁺, N₂H⁺ e DCO⁺, com propriedades parecidas com as de L1544 e de seus congêneres de Touro e Perseu. E existe uma segunda população, de núcleos expostos ao ultravioleta, com poeira mais quente e emissão intensa de CN, HCN e HCO⁺, que um levantamento guiado pelos traçadores de gás frio jamais teria encontrado. As densidades de coluna e de volume das duas populações se sobrepõem em boa parte, embora os núcleos quentes não alcancem os valores mais extremos. Eles ficam nas fronteiras da região HII de NGC 2024, a bolha de hidrogênio ionizado inflada pelas estrelas jovens da Nebulosa da Chama, e na borda densa da nuvem voltada para σ Orionis. Têm química distinta, mas massas semelhantes e, ao que tudo indica, capacidade semelhante de formar estrelas.

Três fontes de NGC 2024, uma delas apontada como discrepante no teste de robustez, complicam a história. Elas têm poeira quente e, mesmo assim, escores negativos no segundo eixo, ao contrário de núcleos semelhantes em G₀/n ou em temperatura, que caem do lado positivo. São objetos mais densos e massivos que quase toda a amostra, com um padrão misto, com N₂H⁺, H¹³CO⁺ e CN intensos e déficit de DCO⁺. A explicação proposta passa pela história da região. A primeira geração de estrelas OB de NGC 2024 teria sido desencadeada pela colisão de duas nuvens, segundo Ohama e colaboradores em 2017, com o ponto de fusão coincidindo com a maior densidade de coluna do mapa, acima de 10²³ por centímetro quadrado. É ali que estão ao menos dois núcleos pré-estelares, o 646 e o 667, e uma protoestrela ativa, a 988, cuja formação pode ter sido disparada pela colisão. Os núcleos que cercam a região HII seriam, então, uma segunda geração, gatilhada pela retroalimentação das estrelas jovens, num processo sequencial parecido com o dos pilares da Nebulosa da Águia, onde Fukuda, Hanawa e Sugitani identificaram formação estelar em cadeia em 2002 e Wen e colaboradores contaram, em 2025, centenas de protoestrelas nas bordas.

O terceiro eixo tem um padrão que os autores classificam de interessante, e a palavra é merecida. Suas cargas positivas são justamente os isotopólogos mais fracos do CO, C¹⁸O, C¹⁷O e ¹³CO, acompanhados de SO e metanol; as negativas são N₂H⁺, CCH, HNC, DNC, DCO⁺, c-C₃H₂, HN¹³C, CCS e HCO. A grandeza associada é a densidade volumétrica média, mas a relação não é monotônica: o escore cresce com a densidade até perto de 3 mil moléculas por centímetro cúbico, 10 elevado a 3,5 na notação dos autores, e então decresce, com grande dispersão acima desse ponto. Abaixo do limiar, o coeficiente de correlação é de 0,490. A subida se explica por excitação colisional: as transições fundamentais de C¹⁸O e C¹⁷O têm densidade crítica moderada, em torno de 2 mil por centímetro cúbico, o valor a partir do qual as colisões com o hidrogênio molecular povoam os níveis com eficiência. Abaixo dela, mais densidade traz mais colisões e linhas mais brilhantes para a mesma quantidade de moléculas; acima, a linha para de responder.

A descida é outra física. Os escores mais negativos do terceiro eixo pertencem aos núcleos de maior densidade, alojados nas regiões mais frias da nuvem, B9 e a crista junto a NGC 2023. Esses objetos brilham em espécies deuteradas, DCO⁺ e DNC, e são pobres em C¹⁸O e C¹⁷O. O fenômeno se chama freeze-out, o congelamento das moléculas sobre os grãos de poeira: quando a temperatura cai abaixo de uns 20 kelvins, a temperatura de sublimação do CO, as moléculas que colidem com um grão ficam presas e formam mantos de gelo, como a geada que se deposita numa vidraça. O tempo típico para uma molécula encontrar um grão é de um bilhão de anos dividido pela densidade em moléculas por centímetro cúbico; a cem mil por centímetro cúbico, isso dá dez mil anos, menos que os cem mil anos de vida de um núcleo. O CO some do gás, e com ele some o principal destruidor do N₂H⁺ e das moléculas deuteradas, que então florescem. As temperaturas de poeira desses núcleos ficam, de fato, perto dos 20 kelvins, e o padrão repete o dos casos extremos de congelamento que Laurent Pagani e colaboradores estudaram em 2024 em nuvens próximas.

Duas razões de intensidade confirmam a leitura. Quando o escore da terceira componente cai abaixo de −3, a razão N₂H⁺/C¹⁸O sobe acima de 3, o valor esperado para o congelamento, com os casos mais extremos em B9: o núcleo 781 e as protoestrelas 993-1, 993-2 e 995. E um grupo peculiar de objetos, com escores entre −1 e −3 e razão G₀/n elevada, mostra excesso de ¹³CO em relação ao C¹⁸O. Esse excesso é a assinatura de dois processos ligados à radiação: o fracionamento do ¹³CO, favorecido pela química de baixa temperatura, e a fotodissociação seletiva do C¹⁸O, que, por ser raro, não consegue se autoblindar da luz ultravioleta como fazem os isotopólogos abundantes, conforme Antoine Roueff e colegas descreveram em 2021. A fonte dessa luz, ao que parece, é a região HII de NGC 2024. Até no eixo da densidade, a radiação deixa marca.

Para quantificar o contraste, os autores recortaram dois subconjuntos de 50 objetos cada, nos extremos do plano escore contra G₀/n. O primeiro reúne os núcleos mais expostos à radiação e mora exclusivamente em torno de NGC 2024: 22 sem estrela e 28 pré-estelares. O segundo reúne os mais blindados, quase todos em NGC 2023 e B9, com alguns em Cloak e na Chama, e é formado sobretudo por núcleos robustos, 39 deles, e por 9 protoestrelas. As distribuições de intensidade dos dois grupos, linha a linha, desenham dois repertórios químicos: o grupo irradiado domina em ¹²CO, HCO⁺, HCN e CN, o grupo protegido domina em N₂H⁺, DCO⁺, metanol e nas demais espécies de gás frio. São duas subcategorias submetidas às condições físicas mais extremas da nuvem, com marcadores químicos, histórias de formação e estágios evolutivos diferentes. Quase toda a literatura sobre núcleos pré-estelares trata do segundo grupo. O primeiro mal foi estudado.

Um apêndice do artigo trata com franqueza de uma limitação. As posições do catálogo de 2020 nem sempre coincidem com os picos da densidade de coluna nem com os picos de N₂H⁺, o traçador consagrado de núcleos frios. Comparando o catálogo com um mapa de N₂H⁺ limpo de ruído, a equipe identificou nove objetos com deslocamento evidente, os de número 453, 546, 547, 618, 652, 792, 829, 987 e 996. O fluxo de N₂H⁺ integrado dentro dos raios catalogados responde por aproximadamente 85 por cento do total do mapa, mas regiões de coluna muito alta, acima de 10²³, como NGC 2024 e NGC 2023, têm emissão estendida que o algoritmo de extração deixou de fora. A razão provável é que a emissão molecular acompanha a poeira mais fria, visível em comprimentos de onda acima de 300 micrômetros, enquanto o algoritmo foi otimizado entre 70 e 500 micrômetros com peso maior para as bandas curtas, de melhor resolução. O catálogo de Helen Kirk e colaboradores, de 2016, posiciona melhor esses casos. Depois de verificar os objetos com N₂H⁺ resolvido, os autores concluíram que as discrepâncias não afetam a metodologia nem as conclusões, e mantiveram as posições originais.

A análise completa, com as 25 componentes, também foi publicada, inclusive as cargas calculadas com e sem os objetos discrepantes e a média sobre todas as rodadas de validação. Em componentes de baixa ordem, como a quinta ou a décima sétima, a remoção dos discrepantes desloca as cargas de forma visível, o que era esperado: sem os pontos extremos, a distribuição de variância se aproxima do eixo verdadeiro. Os diagramas de escores mostram a nuvem de 1.003 objetos encolhendo componente após componente até virar um ponto. Os núcleos robustos e as protoestrelas são os que mais se afastam do centro, sinal de propriedades mais diversas, talvez por um estágio mais evoluído, talvez por influências ambientais distintas. A estatística aponta esses casos, mas, como os próprios autores admitem, não decide o estágio evolutivo de ninguém. Para isso é preciso desemaranhar a estrutura do gás ao longo da linha de visada, como Léontine Ségal e colegas fizeram para a Cabeça de Cavalo em 2024, e seguir com observações de resolução angular suficiente para alcançar as escalas físicas dos núcleos.

O que fica para a prática da astronomia é um parâmetro e um alerta. O parâmetro é G₀/n, barato de calcular sempre que há dados infravermelhos como os do Herschel: G₀ sai da luminosidade no infravermelho distante e a densidade se recupera da densidade de coluna. Ele resume a influência do ambiente sobre a química interna e pode ser acoplado ao critério de Bonnor-Ebert para identificar candidatos a colapso sem o viés de temperatura. O alerta é sobre os traçadores. Levantamentos que procuram núcleos pela emissão de N₂H⁺ ou de H₂D⁺ encontram, por construção, apenas os frios e densos; para completar a amostra com os objetos de G₀/n elevado, é preciso somar traçadores complementares como CN, HCN, CCH e CS. Para quem vai catalogar núcleos em nuvens distantes com o ALMA, onde a maioria dos ambientes se parece mais com Orion B do que com Touro, a lição é direta: só uma diversidade de linhas dá conta da diversidade de condições.

Os dados que sustentam tudo isso estão abertos. A primeira parte do ORION-B está no sistema de gerenciamento de dados do IRAM, e o arquivo deve ganhar, até o fim de 2026, os cubos que cobrem os cinco graus quadrados completos. A tabela da amostra e as intensidades medidas em cada núcleo, com incertezas, ficam disponíveis no Centro de Dados Astronômicos de Estrasburgo, o CDS. O trabalho reúne 29 autores de instituições da França, da Espanha, dos Estados Unidos, do Reino Unido e do Japão, entre eles Victor de Souza Magalhaes, do National Radio Astronomy Observatory, e usou observações de seis projetos do telescópio de 30 metros realizados entre 2013 e 2016, culminando no programa grande 124-16. Recebeu apoio da agência francesa de pesquisa, do programa nacional de física e química do meio interestelar do CNRS, do ministério espanhol da ciência e da NASA, por meio do Laboratório de Propulsão a Jato. O artigo é de acesso aberto sob a licença Creative Commons Atribuição 4.0.

Há uma ironia no fato de que a nuvem mais fotografada do céu tenha guardado por tanto tempo um segredo tão simples. Ao olhar a Cabeça de Cavalo, vemos a borda de Orion B queimando sob a luz de σ Orionis; não vemos os núcleos que, justo nessa borda, tentam formar estrelas com uma química que nenhum manual previa. O estudo de Mazurek e seus colegas não derruba L1544; apenas o devolve ao seu lugar, o de um exemplar entre muitos. Mil e tantos casulos de gás, lidos em 25 línguas moleculares, dizem que a Galáxia não fabrica estrelas com uma receita só. Ela usa a receita que o ambiente permite, e a diferença entre uma e outra está escrita, para quem souber ler, na luz que cada molécula emite no frio.