Binary Neutron Star Merger Evolution and r-Process Enrichment in the Milky Way Disk

Resumo Técnico: Evolução de Fusões de Estrelas de Nêutrons Binárias e Enriquecimento do Processo-r no Disco da Via Láctea

Declaração do Problema
A origem de aproximadamente metade dos elementos de captura rápida de nêutrons (processo-r) no Universo permanece sem resolução. Embora fusões de Estrelas de Nêutrons Binárias (BNS) sejam um candidato principal devido aos ejectos ricos em nêutrons observados no GW170817, estudos anteriores indicam que fusões de BNS sozinhas lutam para reproduzir a história de enriquecimento por processo-r observada em estrelas do disco da Via Láctea. Especificamente, modelos padrão onde a taxa de fusão de BNS segue a História de Formação Estelar (SFH) convoluída com uma Distribuição de Tempo de Atraso (DTD) fixa não conseguem corresponder à tendência decrescente de [Eu/Fe] versus [Fe/H] observada em estrelas do disco. Esta tensão é tipicamente atribuída aos longos tempos de atraso das fusões de BNS em relação à formação estelar, o que impede que elas enriqueçam suficientemente o meio interestelar (ISM) primordial.

Metodologia
Este trabalho estende o quadro bayesiano de multi-mensageiros desenvolvido por H.-Y. Chen et al. (2025) e M. Saleem et al. (2026) para investigar se fusões de BNS podem servir como o único canal de produção de processo-r se a sua eficiência de enriquecimento evoluir com o desvio para o vermelho (redshift) além da SFH.

• Modelo de Evolução Química: Os autores empregam um modelo de evolução química galáctica de uma única zona que rastreia Ferro (Fe) e Európio (Eu). A produção de Ferro é governada por Supernovas de Colapso do Núcleo (CCSNe) e Supernovas Tipo Ia (SNe Ia) com taxas e rendimentos fixos. A produção de Európio é modelada exclusivamente através de fusões de BNS.
• Populações de BNS: A população de BNS é dividida em duas subpopulações: um canal "Atrasado" (que funde após um atraso de tempo em relação à formação estelar, governado por uma DTD de lei de potência) e um canal "Fusão Rápida" (FM) (que funde prontamente com a formação estelar).
• Eficiência Evolutiva: A inovação central é a parametrização de uma eficiência de enriquecimento por processo-r evolutiva, $\epsilon(z)$, modelada como uma lei de potência quebrada em $(1+z)$. Isso permite que a contribuição efetiva das fusões de BNS varie de maneira diferente em diferentes desvios para o vermelho, potencialmente impulsionada por eficiência de formação dependente da metalicidade ou rendimentos médios evolutivos.
• Inferência Bayesiana: O estudo analisa conjuntamente medições de abundância estelar de estrelas do disco da Via Láctea (Battistini & Bensby 2016), restrições de ondas gravitacionais (GW) sobre taxas de fusão locais (GWTC-4.0), DTDs de explosões de raios gama curtas (sGRB) e observações de pulsares. Nove cenários distintos (M1–M9) são testados, variando a inclusão dos canais Atrasado/FM e a presença de evolução em suas eficiências.

Principais Resultados

1. Preferência por Evolução: Cenários que incorporam evolução dependente do desvio para o vermelho na eficiência de enriquecimento de BNS são fortemente preferidos em relação a cenários não evolutivos. Os fatores de Bayes ($B_{evolving}^{non-evolving}$) excedem $10^{20}$, indicando evidência decisiva de que o modelo padrão não evolutivo não pode reproduzir as tendências observadas de [Eu/Fe]–[Fe/H].
2. Reprodução de Tendências de Abundância: Modelos com evolução (por exemplo, M2, M4, M6–M9) reproduzem com sucesso o declínio observado em [Eu/Fe] em altas metalicidades. Isso requer uma contribuição aprimorada inicial de fusões de BNS antes do pico da formação estelar cósmica ($z \sim 2$).
3. Tensão com Observações de sGRB: Embora os cenários evolutivos se ajustem às abundâncias estelares, as densidades de taxas de fusão inferidas em altos desvios para o vermelho ($z \gtrsim 1$) aumentam mais rapidamente do que as taxas inferidas a partir de galáxias hospedeiras de sGRB (Zevin et al. 2022; Pracchia & Salafia 2026). Isso sugere uma tensão potencial entre as taxas de fusão necessárias para explicar a evolução química e aquelas observadas via sGRBs.
4. Consistência com o Fundo Estocástico de Ondas Gravitacionais (SGWB): Apesar das taxas de fusão aprimoradas em altos desvios para o vermelho, os espectros de densidade de energia do SGWB resultante permanecem abaixo dos limites superiores atuais do LIGO–Virgo–KAGRA (LVK) (O1–O4a). O fundo previsto não exclui atualmente esses cenários.
5. Tensão com Síntese Populacional: A evolução requerida implica uma eficiência de formação de BNS que aumenta por um fator de alguns até uma ordem de magnitude em baixas metalicidades ([Fe/H] $\sim -1.5$). Esta tendência contradiz as descobertas de vários estudos de síntese populacional (por exemplo, van Son et al. 2024) que não encontram correlação forte entre eficiência de formação de BNS e metalicidade, embora outros estudos sugiram que tais correlações podem existir.
6. Evolução de Rendimentos: Se a evolução observada for atribuída inteiramente a mudanças no rendimento médio de processo-r em vez da taxa de fusão, o rendimento médio deve aumentar por um fator de 2 a 7 em $z \sim 2$ em relação aos valores locais.

Significado e Alegações
O artigo fornece uma estrutura quantitativa para avaliar se fusões de BNS podem explicar a história de enriquecimento por processo-r da Via Láctea se a sua eficiência evoluir. Os autores afirmam que, embora cenários evolutivos sejam estatisticamente necessários para ajustar os dados de abundância estelar, eles introduzem tensão significativa com restrições astrofísicas independentes, especificamente taxas de sGRB e previsões de síntese populacional relacionadas à dependência da metalicidade.

O estudo conclui que a "evolução" pode manifestar-se como uma taxa de fusão dependente do desvio para o vermelho (por exemplo, formação dependente da metalicidade) ou um rendimento de processo-r evolutivo (por exemplo, mudanças nas distribuições de massa/rotação de estrelas de nêutrons). Os autores notam que, embora a evolução requerida seja robustamente inferida a partir de dados químicos, sua realização física permanece incerta. Eles sugerem que melhorias futuras na sensibilidade do SGWB (por exemplo, a corrida O5) poderiam fornecer restrições diretas sobre esses cenários, e que resolver a tensão com observações de sGRB pode exigir reavaliar os canais de progenitores de sGRBs ou os modelos subjacentes de síntese populacional. O trabalho não afirma resolver definitivamente o problema da origem do processo-r, mas sim delimita os requisitos evolutivos específicos e as tensões observacionais que qualquer modelo bem-sucedido apenas de BNS deve abordar.