Há uma molécula de álcool, do mesmo tipo que se obtém da destilação da madeira, flutuando no casulo morno de uma estrela recém-nascida a cerca de 950 anos-luz da Terra. Nada de extraordinário até aqui: o metanol é a molécula orgânica mais comum dos berçários estelares. O detalhe que fez quinze astrônomos de oito países se debruçarem sobre ela é outro. Em cada uma dessas moléculas, os quatro átomos de hidrogênio foram substituídos por deutério, o hidrogênio pesado que carrega um nêutron a mais no núcleo. É o metanol inteiramente deuterado, CD₃OD, e nunca antes alguém havia conseguido provar sua presença fora dos laboratórios da Terra. Agora há dez linhas espectrais, captadas no deserto do Atacama, que não deixam dúvida.

A prova chegou com o Atacama Large Millimeter/submillimeter Array, o ALMA, no âmbito de um levantamento espectral chamado COMPASS, sigla em inglês para Moléculas Orgânicas Complexas em Protoestrelas com Levantamentos Espectrais do ALMA. O artigo que a descreve, assinado por Arnaud Belloche, do Instituto Max Planck de Radioastronomia, em Bonn, e por mais quatorze colegas, foi aceito para publicação na revista Astronomy & Astrophysics e disponibilizado no repositório arXiv em 5 de outubro de 2026. É o sétimo trabalho da série COMPASS, e o título não deixa margem para interpretações: primeira detecção interestelar do metanol totalmente deuterado. Trata-se também, segundo os autores, da primeira identificação de uma molécula quadruplamente deuterada no meio interestelar.

O que torna o achado mais do que uma curiosidade química é o que ele revela sobre as regras do jogo. O deutério funciona, para os astroquímicos, como um marcador de origem: a proporção entre moléculas deuteradas e suas versões comuns guarda a memória da temperatura, da densidade e da idade do ambiente em que elas se formaram. Ao medir a abundância de todos os isotopólogos deuterados do metanol para os quais existem dados de laboratório, incluindo o CD₃OD, a equipe encontrou um padrão que nenhum modelo químico publicado na última década consegue reproduzir. Existe algo acontecendo na química fria das nuvens escuras que ainda não compreendemos, e a molécula mais pesada da família é justamente a que mais escancara o problema.

Para dimensionar a anomalia, é preciso começar pelo próprio deutério. Esse isótopo do hidrogênio foi forjado nos primeiros minutos após o Big Bang, e desde então as estrelas só fazem destruí-lo em seus interiores. No disco da Via Láctea, a proporção é de aproximadamente 1,6 átomo de deutério para cada 100 mil átomos de hidrogênio, ou um em cerca de 60 mil. Se a química fosse um simples sorteio, uma molécula com um átomo de deutério no lugar de um hidrogênio deveria ser dezenas de milhares de vezes mais rara do que a versão comum. Nas nuvens moleculares densas, essa expectativa falha por várias ordens de grandeza, e as moléculas deuteradas aparecem em quantidades que chegam a ser milhares de vezes maiores do que o sorteio permitiria.

A razão é uma reação química que funciona como uma catraca. Nos núcleos frios das nuvens, a temperaturas próximas de 10 kelvins, o íon H₃⁺ reage com a molécula HD, que carrega um átomo de hidrogênio e um de deutério, e produz H₂D⁺. A reação libera energia, e no frio extremo não há calor suficiente para fazê-la andar no sentido contrário. O deutério, antes trancado no HD, passa assim para outras moléculas do gás. O processo ganha força quando o monóxido de carbono congela sobre os grãos de poeira, porque o CO é justamente o principal destruidor do H₃⁺ e de seus parentes deuterados H₂D⁺, D₂H⁺ e D₃⁺. Menos CO no gás significa mais tempo de vida para esses íons e mais deutério redistribuído.

Os grãos de poeira participam do enriquecimento de outra maneira. Quando os íons deuterados do H₃⁺ se recombinam com elétrons, eles se despedaçam e liberam átomos de deutério soltos, que grudam nas superfícies geladas dos grãos. É ali, sobre essas cascas de gelo, que nasce o metanol: átomos de hidrogênio pousam sobre moléculas de CO congeladas e, de hidrogenação em hidrogenação, constroem CH₃OH. Como os átomos de deutério chegam misturados aos de hidrogênio, parte das moléculas sai da linha de montagem com um ou mais nêutrons extras a bordo. A proporção orto-para do hidrogênio molecular também entra na conta: o orto-H₂, mais energético, consegue reagir com o H₂D⁺ mesmo no frio e desfaz parte do trabalho da catraca, de modo que a deuteração só se torna eficiente quando essa proporção cai.

O metanol ocupa uma posição central nessa história. Ele é a mais abundante das moléculas orgânicas complexas, categoria que os astroquímicos definem como moléculas com carbono e pelo menos seis átomos, e é a matéria-prima de reações de superfície que levam a espécies maiores, como o formato de metila e o éter dimetílico. Por ter quatro hidrogênios, admite sete variantes deuteradas distintas, do CH₂DOH e do CH₃OD, com um único deutério, até o CD₃OD, com quatro. Cada uma dessas variantes, chamadas de isotopólogos, conta uma parte diferente da história: as que carregam deutério no grupo metila, o CH₃, e as que o carregam no grupo hidroxila, o OH, respondem a processos químicos distintos, e a comparação entre elas funciona como um teste cruzado das reações que operam no gelo.

A busca pelo metanol pesado tem quase quatro décadas. As versões com um único deutério, CH₃OD e CH₂DOH, foram detectadas pela primeira vez na região de formação de estrelas massivas Orion IRc2, em 1988 e 1993, com radiotelescópios de antena única. As versões duplamente e triplamente deuteradas, CHD₂OH e CD₃OH, apareceram depois, em 2002 e 2006, no sistema protoestelar de baixa massa IRAS 16293−2422, na nuvem de Ofiúco. Desde então, os quatro isotopólogos foram reencontrados em escalas menores e em outros objetos, incluindo os dois componentes da binária IRAS 16293 e as protoestrelas HOPS 373SW e BHR71-IRS1. Faltava o último degrau da escada, aquele em que não sobra nenhum hidrogênio comum na molécula.

Esse degrau resistia por dois motivos. O primeiro é a raridade: a cada substituição de hidrogênio por deutério, a molécula fica mais escassa, e o CD₃OD é a mais rara de todas. O segundo é a falta de um mapa. Para identificar uma molécula no espaço, é preciso conhecer de antemão, com precisão de laboratório, as frequências exatas em que ela emite ao girar. Só em 2023 um grupo liderado por Vadim Ilyushin publicou essa espectroscopia para o CD₃OD na faixa submilimétrica, e os primeiros indícios de sua presença surgiram quase no mesmo instante em IRAS 16293 e em BHR71-IRS1. Eram pistas, não provas: linhas fracas demais, ou misturadas demais à emissão de outras moléculas, para sustentar uma identificação robusta. Dois outros isotopólogos, o CH₂DOD e o CHD₂OD, continuam sem previsões espectroscópicas em bancos de dados públicos e permanecem, por isso, invisíveis.

Há uma razão adicional para perseguir essas moléculas com tanto empenho. O grau de fracionamento do deutério, isto é, a razão entre a quantidade de uma molécula deuterada e a de sua versão com hidrogênio, cresce com o tempo durante a fase pré-estelar de um núcleo denso, à medida que a temperatura cai, o CO congela e a proporção orto-para do H₂ diminui. Uma molécula muito deuterada formou-se, em geral, mais tarde e em condições mais frias que uma pouco deuterada. Existe evidência, por exemplo, de que a água deuterada se forma numa etapa posterior à da água comum. O padrão de deuteração funciona, portanto, como um registro fóssil da história de uma região de formação estelar, e cada isotopólogo novo acrescenta uma camada a esse registro.

O COMPASS foi desenhado para ler esse registro com uma resolução inédita. Trata-se de um large program do ALMA, categoria reservada aos projetos de maior porte do observatório, voltado a onze sistemas protoestelares de baixa massa que abrigam hot corinos, termo que designa as regiões quentes e compactas, com algumas dezenas de unidades astronômicas, onde o gelo dos grãos sublima e o gás se torna brilhante na emissão rotacional de moléculas orgânicas complexas. Entre janeiro de 2023 e agosto de 2024, o conjunto de antenas no deserto do Atacama varreu, sem interrupções, a faixa de frequências de 279,0 a 311,7 gigahertz, usando nove configurações espectrais com quatro janelas cada e uma resolução de 488 quilohertz, o equivalente a cerca de 0,5 quilômetro por segundo em velocidade.

A resolução angular do levantamento varia entre 0,35 por 0,20 segundo de arco e 0,62 por 0,37 segundo de arco, com um feixe mediano de 0,39 por 0,23 segundo de arco. À distância do alvo, isso corresponde a enxergar detalhes de pouco mais de cem unidades astronômicas, algo como três vezes a distância entre o Sol e Netuno. O ruído por canal fica entre 1,6 e 3,3 milijanskys por feixe, com mediana de 2,5, uma sensibilidade que permite captar linhas espectrais com apenas alguns milésimos da intensidade das transições mais brilhantes do metanol comum. Sem esse salto de sensibilidade, a caça ao isotopólogo mais raro do metanol não teria como sair do terreno das pistas.

O alvo escolhido foi IRAS4A, um sistema binário de protoestrelas na nuvem molecular NGC 1333, parte do complexo de Perseu, a 290 parsecs do Sol, ou cerca de 950 anos-luz. É uma das protoestrelas de Classe 0 mais estudadas do céu, tão jovem que ainda está envolta na maior parte do material do qual se forma. A escolha não foi aleatória. Em 2002, foi exatamente em IRAS4A que se identificou pela primeira vez a amônia totalmente deuterada, ND₃, numa região protoestelar. Se havia um lugar onde o metanol de quatro deutérios teria chance de aparecer, era ali. A análise concentrou-se no componente IRAS4A2, o mais rico em linhas moleculares dos dois.

A equipe extraiu o espectro na posição de pico de um mapa de contagem de canais, que destaca as regiões mais ricas em emissão molecular, excluídas as espécies simples. Essa posição fica a apenas 0,07 segundo de arco do pico de emissão contínua de IRAS4A2. O espectro é uma verdadeira floresta: para interpretá-lo, os pesquisadores construíram um modelo preliminar da emissão molecular da fonte que já inclui a contribuição de 154 entradas de catálogos espectroscópicos, reunidas a partir dos bancos de dados de Colônia (CDMS), de Lille (LSD) e do Laboratório de Propulsão a Jato (JPL), somadas a entradas privadas. É como tentar ouvir um único instrumento numa orquestra de centenas: antes de afirmar que o oboé está tocando, é preciso saber o que cada um dos outros músicos faz.

Para isso, o grupo recorreu ao pacote Weeds, uma extensão do programa CLASS da suíte GILDAS, que calcula a transferência radiativa ao longo da linha de visada sob a hipótese de equilíbrio termodinâmico local, ou LTE. Na prática, a hipótese supõe que as colisões no gás são frequentes o bastante para que uma única temperatura descreva a população dos níveis de energia de cada molécula. O modelo de cada espécie depende de cinco parâmetros: o tamanho da região emissora, a densidade de coluna, que mede quantas moléculas há num cilindro de um centímetro quadrado de seção ao longo da linha de visada, a temperatura rotacional, o deslocamento de velocidade e a largura das linhas. A emissão de cada molécula é modelada separadamente, e as contribuições são somadas no fim.

Quatro desses parâmetros foram fixados antes de qualquer ajuste. As temperaturas vieram de diagramas de população, que relacionam a intensidade de cada linha à energia do nível que a produz, usando apenas transições pouco opacas e pouco contaminadas. O resultado foi notavelmente uniforme: 207,4 kelvins para o metanol comum, 198,4 para o ¹³CH₃OH, 186,8 para o CH₂DOH, 197,2 para o CH₃OD, 193 para o CHD₂OH, 197,5 para o CD₃OH e 202 para o CD₃OD, com incertezas estatísticas que vão de 5 a 41 kelvins. A única exceção, os 122 kelvins do CH₃¹⁷OH, foi atribuída à contaminação de duas linhas brilhantes por espécies não identificadas. A equipe adotou 200 kelvins para todos, uma temperatura de cerca de 73 graus Celsius negativos que, na escala das nuvens moleculares, é um verdadeiro forno.

O tamanho da região emissora saiu do ajuste de funções gaussianas bidimensionais aos mapas de intensidade integrada de transições limpas. Descontado o efeito do feixe, os isotopólogos apresentam tamanhos de 0,20 a 0,22 segundo de arco e orientação de 90 a 100 graus a leste do norte, em bom acordo entre si. Como as linhas mais opacas tendem a parecer ligeiramente maiores, os autores adotaram um valor um pouco menor, de 0,18 segundo de arco, que a 290 parsecs equivale a cerca de 52 unidades astronômicas, dimensão comparável à da órbita de Netuno. As linhas opticamente finas forneceram o restante: uma largura de 2,0 quilômetros por segundo e uma velocidade sistêmica de 6,8 quilômetros por segundo ajustam bem todos os isotopólogos. Restou livre apenas a densidade de coluna de cada molécula, otimizada visualmente até o limite em que o espectro sintético não superestima nenhuma linha observada.

Foi nesse cenário que o CD₃OD apareceu. Dez linhas do isotopólogo, entre 288,9 e 301,8 gigahertz, passaram pelos três critérios da equipe: intensidade de pico acima de três vezes o ruído, coincidência de frequência com um pico claro no espectro observado e ausência de mistura severa com outras moléculas. Os níveis de energia envolvidos vão de 65,6 a 372,7 kelvins acima do estado fundamental, uma faixa larga o bastante para confirmar que o gás está mesmo a cerca de 200 kelvins. Duas das dez linhas, em 291.867 e 291.983 megahertz, são parcialmente contaminadas, a primeira por formato de metila e éter dimetílico em estados vibracionais excitados, a segunda por etilenoglicol e pelo próprio CHD₂OH. Cerca de metade das linhas tem razão sinal-ruído em torno de 6, e as demais são mais fortes.

Dez linhas podem parecer pouco, mas é a combinação que importa. O espectro rotacional de uma molécula é uma espécie de código de barras: a posição e a intensidade relativa de cada linha são determinadas pela geometria e pela massa da molécula, e nenhuma outra espécie reproduz exatamente o mesmo conjunto. Encontrar dez barras do código na frequência certa, com as intensidades previstas pelo mesmo modelo de temperatura e tamanho, é o que separa uma detecção de uma coincidência. A densidade de coluna resultante é de 2,0 × 10¹⁶ moléculas por centímetro quadrado, com incerteza conservadora de até 15%. O estado torsional excitado do CD₃OD não foi detectado, com limite superior de 6,0 × 10¹⁶. Os modelos preliminares das outras dez fontes do COMPASS indicam que a molécula não aparece em nenhuma delas.

Todos os isotopólogos do metanol com dados espectroscópicos públicos foram claramente identificados em IRAS4A2: o CH₃OH comum, com 97 linhas em três estados torsionais, o ¹³CH₃OH, com 49, o CH₃¹⁸OH, com 11 mais uma tentativa, o CH₃¹⁷OH, com 4, o CH₂DOH, com 47, o CH₃OD, com 46, o CHD₂OH, com 76, e o CD₃OH, com 74. As densidades de coluna contam a história em números: 5,0 × 10¹⁸ para o metanol, 5,5 × 10¹⁷ para a versão com carbono-13, 1,0 × 10¹⁷ para a de oxigênio-18, 2,5 × 10¹⁶ para a de oxigênio-17, 1,2 × 10¹⁸ para o CH₂DOH, 3,0 × 10¹⁷ para o CH₃OD, 1,0 × 10¹⁸ para o CHD₂OH e 4,5 × 10¹⁷ para o CD₃OH, sempre em moléculas por centímetro quadrado. É o inventário mais completo da família do metanol já feito numa única protoestrela.

A.1 · CH₃¹⁸OH

A.2 · CH₃¹⁷OH

A.3 · CH₂DOH

A.4 · CH₃OD

A.5 · CHD₂OH

A.6 · CD₃OH

Foi ao comparar esses números que surgiu uma segunda anomalia, inesperada e, à primeira vista, embaraçosa. A razão entre o metanol comum e sua versão com carbono-13 ficou em 9, quando a proporção isotópica ¹²C/¹³C esperada à distância do Sol ao centro galáctico é de cerca de 68. A razão entre o metanol e a versão com oxigênio-18 deu 50, contra os aproximadamente 560 esperados para a razão ¹⁶O/¹⁸O do meio interestelar local. Em ambos os casos, o metanol comum parece sete a onze vezes menos abundante do que deveria em relação a seus isotopólogos raros. Ou a química do carbono e do oxigênio em IRAS4A2 é profundamente diferente da do resto da vizinhança solar, ou alguma coisa está escondendo moléculas de metanol da análise.

A pista decisiva veio dos próprios isotopólogos raros, que concordam entre si. A razão ¹³CH₃OH/CH₃¹⁸OH é de 5,5, apenas um pouco abaixo dos cerca de 8 esperados, e a razão CH₃¹⁸OH/CH₃¹⁷OH é de 4,0, em pleno acordo com os valores de 4,11 ± 0,14 medidos em Ofiúco, de 4,0 ± 0,5 em Touro e de cerca de 3,8 previstos pelo gradiente galáctico. Um estudo publicado em 2025 pela equipe de L. A. Busch havia obtido razões igualmente baixas, entre 4 e 30, com o interferômetro NOEMA, para o metanol e o cianeto de metila em protoestrelas da vizinhança solar, incluindo fontes da mesma nuvem, e as interpretou como um enriquecimento real em carbono-13. Belloche e colegas discordam: para que essa leitura valesse em IRAS4A2, o metanol teria de estar enriquecido também em oxigênio-18 e oxigênio-17, nas mesmas proporções, algo improvável do ponto de vista químico.

A conclusão mais provável é que a análise subestima a verdadeira densidade de coluna do metanol. O suspeito é a estrutura da emissão em escalas menores que o feixe. A equipe assumiu uma fonte monolítica de 0,18 segundo de arco, mas uma emissão grumosa, feita de vários nós compactos espalhados pela mesma área, teria o mesmo tamanho aparente após a deconvolução e um fator de preenchimento efetivo menor. Isso elevaria a opacidade real das linhas e contaminaria os diagramas de população com transições mais espessas do que o modelo supõe. As seis linhas menos opacas do metanol no diagrama têm opacidades entre 0,06 e 0,3; para que todas se tornassem significativamente espessas, a grumosidade teria de multiplicar as opacidades por mais de 30. Só observações com resolução angular maior poderão confirmar ou descartar a hipótese.

Enquanto isso, os autores adotaram uma solução pragmática: usar o ¹³CH₃OH e o CH₃¹⁸OH como substitutos do metanol comum, multiplicando suas densidades de coluna pelos fatores 68 e 560. Todas as razões de deuteração do artigo foram calculadas das duas maneiras, e as duas contam a mesma história. O custo é uma incerteza maior. Com a calibração do ALMA, estimada em cerca de 10% a partir das variações do contínuo entre os conjuntos de dados, mais as incertezas de temperatura, tamanho e largura de linha, a incerteza total em cada densidade de coluna chega a cerca de 40%. Nas razões entre isotopólogos, boa parte desses erros se cancela: a incerteza fica em torno de 25% quando o carbono-13 serve de referência, 40% com o oxigênio-18, 18% para o CD₃OD em relação a outro isotopólogo deuterado e 15% nos demais casos.

Os números de deuteração, enfim. No grupo hidroxila, a razão CH₃OD/CH₃OH vale 0,80 ± 0,20% com a referência de carbono-13, ou 0,54 ± 0,21% com a de oxigênio-18: cerca de uma molécula em 150 carrega deutério no lugar do hidrogênio ligado ao oxigênio. Quando o metanol já tem o grupo metila completamente deuterado, o quadro muda. A razão CD₃OD/CD₃OH, que mede a deuteração da hidroxila nas moléculas de CD₃OH, sobe para 4,4 ± 0,8%, cinco a oito vezes mais do que no metanol comum. Esse valor não depende de qual substituto do metanol se usa, porque envolve apenas dois isotopólogos deuterados medidos diretamente, e por isso é um dos resultados mais sólidos do trabalho.

No grupo metila, a escada sobe degrau a degrau. A razão CH₂DOH/CH₃OH é de 3,2 ± 0,8%, ou 2,1 ± 0,9% com a referência de oxigênio-18; a razão CHD₂OH/CH₃OH, de 2,7 ± 0,7%; a razão CD₃OH/CH₃OH, de 1,2 ± 0,3%. No topo, a razão CD₃OD/CH₃OH fica em 0,053 ± 0,013%, ou 0,036 ± 0,014% com o oxigênio-18: o metanol de quatro deutérios é cerca de 2 mil vezes mais raro que o metanol comum, uma abundância relativa de aproximadamente 4 × 10⁻⁴. Entre os próprios isotopólogos deuterados, as razões são ainda mais eloquentes: CHD₂OH/CH₂DOH vale 83 ± 13%, CD₃OH/CHD₂OH vale 45 ± 7% e CD₃OD/CH₃OD vale 6,7 ± 1,2%. A cada degrau que se sobe, a queda de abundância é menor do que a anterior.

Esses percentuais ganham sentido quando traduzidos em razões D/H. Se cada um dos três hidrogênios do grupo metila tivesse, de forma independente, a mesma probabilidade de ser substituído por deutério, o número de moléculas com uma única substituição seria três vezes essa probabilidade, o de moléculas com duas seria proporcional ao seu quadrado e o de moléculas com três, ao seu cubo. Invertendo a conta, cada isotopólogo revela a probabilidade que lhe daria origem. O CH₂DOH aponta para uma razão D/H de 1,1 ± 0,3%; o CHD₂OH, para 9,4 ± 1,2%; o CD₃OH, para 23 ± 2%; e o CD₃OD em relação ao CH₃OD, para 41 ± 2%. Com a referência de oxigênio-18, os valores são 0,71%, 7,7% e 20%, e a razão derivada do CD₃OD em relação ao metanol fica entre 13,7% e 15,2%.

Para sentir o tamanho da discrepância, basta um exercício com dados de jogar. Se a chance de cada hidrogênio virar deutério fosse de 1,1%, como sugere o CH₂DOH, a chance de os quatro serem substituídos ao mesmo tempo seria o produto de quatro fatores de 0,011, ou cerca de 1,5 molécula em cem milhões. O ALMA mediu uma em cada 2 mil. É como se, ao lançar quatro dados de cem faces, o resultado de quatro faces iguais saísse dezenas de milhares de vezes mais do que a probabilidade permite. A conta é simplificada, mas o sentido é claro: quanto mais deutério uma molécula já carrega, maior parece ser a chance de ela receber ainda mais. O processo não é um sorteio; é uma espiral.

É exatamente isso que mostra a figura central do artigo, que reúne as razões D/H de IRAS4A2 e de outras protoestrelas de baixa massa. Se a deuteração do metanol seguisse uma distribuição puramente estatística a partir de uma razão D/H atômica fixa, todos os pontos de uma mesma fonte cairiam sobre uma linha horizontal. Em vez disso, os valores sobem de cerca de 0,01 para os isotopólogos simplesmente deuterados até perto de 0,1 ou mais para os múltipla e totalmente deuterados. A conclusão dos autores é direta: nas protoestrelas de baixa massa, a deuteração do metanol até sua forma completa não obedece a uma distribuição com probabilidade fixa. Alguma propriedade cinética ou termodinâmica da rede química privilegia as moléculas que já são pesadas.

A assimetria entre os dois grupos da molécula oferece uma pista sobre quais reações estão em jogo. A razão CH₂DOH/CH₃OD em IRAS4A2 é de 4, acima do valor 3 que uma repartição ao acaso produziria, já que o metila tem três posições e a hidroxila, uma. Esse excesso se repete em outras protoestrelas de massa baixa e intermediária e tem explicação experimental: experimentos com análogos de gelo interestelar, conduzidos em 2005 por Akira Nagaoka, Naoki Watanabe e Akira Kouchi, na Universidade de Hokkaido, mostraram que reações de abstração e substituição de hidrogênio por deutério favorecem o grupo metila e que, uma vez incorporado, o deutério resiste mais à remoção do que o hidrogênio. Em regiões de formação de estrelas massivas, a razão cai para perto de 1, e uma hipótese de 2015 atribui a diferença a trocas de hidrogênio e deutério entre o metanol e a água dos gelos durante o aquecimento.

Comparado a outras fontes, IRAS4A2 se encaixa num padrão coletivo. Seu perfil de deuteração é quase idêntico ao de BHR71-IRS1 e concorda dentro de um fator de 3 com os de IRAS 16293A, IRAS 16293B e HOPS 373SW. A exceção é V883 Ori, um objeto eruptivo cujo disco foi aquecido por um surto de acreção e que exibe razões muito menores. A semelhança com BHR71-IRS1 tem uma consequência prática: o limite superior de CD₃OD reportado nessa fonte em 2026 deve estar perto da densidade de coluna real, o que faz dela a próxima candidata a uma detecção. Em IRAS 16293B, onde todos os isotopólogos deuterados, exceto o CH₂DOH, são de 1,5 a 3 vezes menos abundantes do que em IRAS4A2, o CD₃OD deve estar algumas vezes abaixo do limite atual.

O passo seguinte foi confrontar esses números com a teoria, e é aí que o artigo se torna desconfortável para os modelos. Os valores medidos no hot corino refletem o gás de uma região onde os gelos já sublimaram, e a premissa dos autores é que esse gás ainda guarda, ao menos em parte, o padrão de deuteração alcançado nos mantos de gelo ao fim da fase pré-estelar, antes que a protoestrela acendesse e evaporasse tudo. Por isso a comparação incluiu tanto modelos de envelopes protoestelares quanto modelos de núcleos pré-estelares publicados ao longo da última década. Nenhum deles sobreviveu intacto ao confronto com os dados de IRAS4A2.

O modelo de Vianney Taquet e colaboradores, de 2014, acompanha um núcleo denso em colapso da fase pré-estelar à protoestelar e descreve os mantos de gelo camada por camada, como as cascas de uma cebola. Ele inclui a queda da razão orto-para do H₂ com o tempo, que aumenta a fração de deutério entre os átomos que pousam nos grãos, e as reações de abstração e substituição medidas em laboratório. Na comparação com IRAS4A2, a 50 unidades astronômicas de escala e ao longo de toda a fase de Classe 0, o modelo reproduz bem as razões do CH₂DOH e do CH₃OD, mas subestima as do CHD₂OH e do CD₃OH por uma a duas ordens de grandeza. O CD₃OD não foi calculado, e os autores consideram provável que ficasse igualmente abaixo do observado.

O modelo de Alexandre Faure e colaboradores, de 2015, testa outra ideia: trocas de isótopos de hidrogênio entre a água e o grupo hidroxila do metanol nos gelos, ativadas termicamente durante o aquecimento causado pela protoestrela. Esse modelo parte de uma distribuição estatística fixada na fase pré-estelar, com razão D/H de 12%, valor derivado de uma medida antiga do CH₂DOH em IRAS 16293 que, hoje se sabe, estava superestimada por um fator de 2. Com uma razão HDO/H₂O inicial de 2% nos gelos, ele acerta o CHD₂OH e o CH₃OD, superestima o CH₂DOH por um fator de 10 e subestima o CD₃OH e o CD₃OD por outro fator de 10. Ajustar a razão inicial aos 1,2% medidos em IRAS4A2 derrubaria todas as previsões, deixando as formas tri e tetradeuteradas duas ordens de grandeza abaixo do observado.

Um modelo cinético mais simples, publicado em 2019 por A. M. Cooper e Johannes Kästner, examina o estado estacionário a 30 kelvins usando novas constantes de velocidade para as reações de abstração de hidrogênio e deutério do grupo metila. Suas previsões dependem de um único parâmetro combinado, o produto entre a razão dos coeficientes de difusão do deutério e do hidrogênio sobre os grãos, estimada em cerca de 0,1, e a razão entre as abundâncias dos dois átomos na superfície. Se o modelo estivesse certo, as três razões medidas em IRAS4A2 deveriam cruzar as curvas teóricas num mesmo valor desse parâmetro. Elas cruzam em valores que se espalham por mais de uma ordem de grandeza. A abstração a partir do metila, sozinha, não dá conta do padrão.

O modelo mais elaborado de núcleo pré-estelar, apresentado em 2022 por Beatrice Kulterer e colegas, segue a química de gás e grãos por 300 mil anos e inclui, como o de Taquet, as reações de abstração e substituição ao longo da rota de hidrogenação do CO. A 15 kelvins, suas previsões para o CH₂DOH e o CH₃OD ficam apenas um fator de alguns abaixo das medidas. Para o CHD₂OH e o CD₃OH, a diferença é de duas ordens de grandeza. Para o CD₃OD, chega a cinco ordens de grandeza: cem mil vezes menos do que o ALMA encontrou. Se o gás do hot corino preserva a química dos gelos pré-estelares, falta a esse tipo de modelo algum ingrediente capaz de empurrar a deuteração até o fim.

Nem o modelo mais recente escapa. Desenvolvido por Wiebke Riedel e colegas entre 2023 e 2025 para o núcleo pré-estelar L1544, com temperatura de poeira caindo de 15 kelvins nas bordas a 6 no centro, ele incorpora dessorção reativa e processos não difusivos. A versão com reações de abstração produz razões CH₂DOH/CH₃OH de 5 a 7%, comparáveis às de IRAS4A2, mas prevê CHD₂OH/CH₂DOH de cerca de 10%, quase uma ordem de grandeza abaixo dos 83% medidos. Rotas de deuteração em fase gasosa a temperaturas mornas, de 50 a 70 kelvins, impulsionadas pelo íon CH₂D⁺, tampouco resolvem: suas escalas de tempo, superiores a 100 mil anos, são longas demais para alterar o padrão do metanol depois da sublimação. O balanço é unânime: todos os modelos subestimam, por várias ordens de grandeza, a produção das formas múltipla e totalmente deuteradas.

O problema não é exclusivo do metanol. O formaldeído, H₂CO, é uma etapa intermediária da mesma hidrogenação do CO, e seu padrão de deuteração conta uma história parecida. Em IRAS 16293B, as razões HDCO/H₂CO e D₂CO/HDCO valem 6,5% e 13%; o modelo de Taquet prevê 4,1% e 4,4%, acertando a primeira e errando a segunda por um fator de 3. Em IRAS4A, mapas publicados em 2024 indicam uma razão de intensidades D₂CO/HDCO de cerca de 0,6 na direção de IRAS4A2, o que, para emissão opticamente fina a 150 ou 200 kelvins, se traduz numa razão de densidades de coluna de cerca de 70%. A análise preliminar dos dados do COMPASS sugere que a deuteração simples e dupla do formaldeído acompanha a do metanol. Os modelos, mais uma vez, ficam para trás.

Para que serve, afinal, saber que o metanol de uma protoestrela de Perseu carrega deutério demais? A resposta está na palavra herança. As moléculas orgânicas que hoje encontramos em cometas, meteoritos e, por extensão, nos oceanos e na biosfera da Terra podem ter sido forjadas antes do nascimento do Sol, na mesma química gelada que produz o CD₃OD. A água dos oceanos terrestres tem uma razão D/H de 1,56 × 10⁻⁴, dez vezes maior que a do gás galáctico primordial, um excesso que só se explica por processos de fracionamento semelhantes aos descritos aqui. Medições da sonda Rosetta no cometa 67P/Churyumov-Gerasimenko indicaram que o metanol cometário carrega assinaturas de deuteração compatíveis com uma origem nas superfícies de grãos da fase pré-estelar. O deutério é a tinta com que a história dessa herança foi escrita, e o metanol pesado de IRAS4A2 é uma página inteira dela.

Os próximos capítulos já estão desenhados. Os autores consideram altamente provável que os dois isotopólogos ainda invisíveis, CH₂DOD e CHD₂OD, sejam detectados sem dificuldade em IRAS4A2 assim que seus espectros forem medidos em laboratório, o que transfere para os espectroscopistas a próxima jogada. Observações do ALMA em configurações mais estendidas poderão testar se a emissão do metanol é grumosa em escalas menores que 50 unidades astronômicas e resolver o enigma da densidade de coluna. Os modeladores, por sua vez, precisam encontrar o ingrediente que falta, seja uma reação ainda não considerada, seja um mecanismo que concentre o deutério nas moléculas já pesadas. E a comparação com o formaldeído, cuja análise completa os autores deixam para um trabalho futuro, promete um segundo teste independente.

O trabalho reúne quinze pesquisadores de doze instituições em oito países, da Alemanha ao Japão, passando por França, Suíça, Dinamarca, Coreia do Sul, Estados Unidos e Taiwan, e usa dados do projeto ALMA 2022.1.00316.L. O material suplementar, com os espectros completos de cada isotopólogo, está disponível publicamente no repositório Zenodo, e a cadeia de redução de dados do COMPASS foi igualmente arquivada em acesso aberto. O retorno do levantamento vai muito além desta detecção: com onze fontes observadas de forma homogênea na mesma faixa de frequências, o COMPASS se tornou a referência para comparar, lado a lado, a química de berçários estelares que até então eram estudados um de cada vez, com instrumentos, resoluções e métodos diferentes.

Há algo de vertiginoso em pensar no caminho dessa molécula. Cada um de seus quatro átomos de deutério foi arrancado de uma molécula de HD por uma reação que só funciona no escuro e no frio, viajou pelo gás como parte de um íon, foi liberado por um elétron, pousou sobre um grão de poeira a dez graus acima do zero absoluto e ali, átomo por átomo, ajudou a construir um álcool. Depois, uma estrela nasceu, o gelo evaporou e a molécula inteira, com seus quatro nêutrons extras, começou a girar e a emitir ondas de rádio que cruzaram 950 anos-luz até as antenas do Atacama. O que elas trazem não é apenas uma detecção inédita. É a prova de que o Universo escreve química num idioma que ainda estamos aprendendo a ler, e de que as palavras mais pesadas são, às vezes, as que mais dizem.