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

Resumen Técnico: Evolución de la Fusión de Estrellas de Neutrones Binarias y Enriquecimiento del Proceso-r en el Disco de la Vía Láctea

Planteamiento del Problema
El origen de aproximadamente la mitad de los elementos de captura rápida de neutrones (proceso-r) en el Universo sigue sin resolverse. Aunque las fusiones de estrellas de neutrones binarias (BNS) son un candidato principal debido a los efluentes ricos en neutrones observados en GW170817, estudios previos indican que las fusiones de BNS por sí solas tienen dificultades para reproducir la historia de enriquecimiento del proceso-r observada en las estrellas del disco de la Vía Láctea. Específicamente, los modelos estándar donde la tasa de fusión de BNS sigue la Historia de Formación Estelar (SFH) convolucionada con una Distribución de Tiempos de Retraso (DTD) fija, no logran coincidir con la tendencia decreciente de [Eu/Fe] frente a [Fe/H] observada en las estrellas del disco. Esta tensión se atribuye típicamente a los largos tiempos de retraso de las fusiones de BNS en relación con la formación estelar, lo que les impide enriquecer suficientemente el medio interestelar (ISM) temprano.

Metodología
Este trabajo extiende el marco bayesiano de mensajería múltiple desarrollado por H.-Y. Chen et al. (2025) y M. Saleem et al. (2026) para investigar si las fusiones de BNS pueden servir como el único canal de producción del proceso-r si su eficiencia de enriquecimiento evoluciona con el corrimiento al rojo más allá de la SFH.

• Modelo de Evolución Química: Los autores emplean un modelo de evolución química galáctica de una sola zona que rastrea el Hierro (Fe) y el Europio (Eu). La producción de Hierro está gobernada por Supernovas de Colapso del Núcleo (CCSNe) y Supernovas de Tipo Ia (SNe Ia) con tasas y rendimientos fijos. La producción de Europio se modela exclusivamente a través de fusiones de BNS.
• Poblaciones de BNS: La población de BNS se divide en dos subpoblaciones: un canal "Retrasado" (que se fusiona después de un tiempo de retraso en relación con la formación estelar, gobernado por una DTD de ley de potencia) y un canal de "Fusión Rápida" (FM) (que se fusiona de manera inmediata con la formación estelar).
• Eficiencia Evolutiva: La innovación central es la parametrización de una eficiencia de enriquecimiento del proceso-r evolutiva, $\epsilon(z)$, modelada como una ley de potencia rota en $(1+z)$. Esto permite que la contribución efectiva de las fusiones de BNS varíe de manera diferente en distintos corrimientos al rojo, impulsada potencialmente por una eficiencia de formación dependiente de la metalicidad o rendimientos promedio evolutivos.
• Inferencia Bayesiana: El estudio analiza conjuntamente mediciones de abundancia estelar de estrellas del disco de la Vía Láctea (Battistini & Bensby 2016), restricciones de ondas gravitacionales (GW) sobre las tasas de fusión locales (GWTC-4.0), DTDs de estallidos cortos de rayos gamma (sGRB) y observaciones de púlsares. Se prueban nueve escenarios distintos (M1–M9), variando la inclusión de canales Retrasados/FM y la presencia de evolución en sus eficiencias.

Resultados Clave

1. Preferencia por la Evolución: Los escenarios que incorporan evolución dependiente del corrimiento al rojo en la eficiencia de enriquecimiento de BNS son fuertemente preferidos sobre los escenarios no evolutivos. Los factores de Bayes ($B_{evolving}^{non-evolving}$) superan $10^{20}$, lo que indica evidencia decisiva de que el modelo estándar no evolutivo no puede reproducir las tendencias observadas de [Eu/Fe]–[Fe/H].
2. Reproducción de Tendencias de Abundancia: Los modelos con evolución (por ejemplo, M2, M4, M6–M9) reproducen con éxito el declive observado en [Eu/Fe] a altas metalicidades. Esto requiere una contribución temprana mejorada de las fusiones de BNS antes del pico de la formación estelar cósmica ($z \sim 2$).
3. Tensión con Observaciones de sGRB: Si bien los escenarios evolutivos ajustan las abundancias estelares, las densidades de tasas de fusión inferidas a alto corrimiento al rojo ($z \gtrsim 1$) aumentan más rápidamente que las tasas inferidas a partir de galaxias anfitrionas de sGRB (Zevin et al. 2022; Pracchia & Salafia 2026). Esto sugiere una tensión potencial entre las tasas de fusión requeridas para explicar la evolución química y las observadas mediante sGRB.
4. Consistencia con el Fondo Estocástico de Ondas Gravitacionales (SGWB): A pesar de las tasas de fusión mejoradas a alto corrimiento al rojo, los espectros de densidad de energía del SGWB resultante permanecen por debajo de los límites superiores actuales de LIGO–Virgo–KAGRA (LVK) (O1–O4a). El fondo predicho no descarta actualmente estos escenarios.
5. Tensión con la Síntesis de Poblaciones: La evolución requerida implica una eficiencia de formación de BNS que aumenta en un factor de unos pocos a un orden de magnitud a bajas metalicidades ([Fe/H] $\sim -1.5$). Esta tendencia contradice los hallazgos de varios estudios de síntesis de poblaciones (por ejemplo, van Son et al. 2024) que no encuentran una correlación fuerte entre la eficiencia de formación de BNS y la metalicidad, aunque otros estudios sugieren que tales correlaciones podrían existir.
6. Evolución del Rendimiento: Si la evolución observada se atribuye enteramente a cambios en el rendimiento promedio del proceso-r en lugar de la tasa de fusión, el rendimiento promedio debe aumentar en un factor de 2 a 7 para $z \sim 2$ en relación con los valores locales.

Significado y Afirmaciones
El artículo proporciona un marco cuantitativo para evaluar si las fusiones de BNS pueden explicar la historia de enriquecimiento del proceso-r de la Vía Láctea si su eficiencia evoluciona. Los autores afirman que, aunque los escenarios evolutivos son estadísticamente necesarios para ajustar los datos de abundancia estelar, introducen una tensión significativa con restricciones astrofísicas independientes, específicamente las tasas de sGRB y las predicciones de síntesis de poblaciones sobre la dependencia de la metalicidad.

El estudio concluye que la "evolución" podría manifestarse ya sea como una tasa de fusión dependiente del corrimiento al rojo (por ejemplo, formación dependiente de la metalicidad) o como un rendimiento del proceso-r evolutivo (por ejemplo, cambios en las distribuciones de masa/giro de las estrellas de neutrones). Los autores señalan que, aunque la evolución requerida se infiere robustamente a partir de datos químicos, su realización física sigue siendo incierta. Sugieren que futuras mejoras en la sensibilidad del SGWB (por ejemplo, la corrida O5) podrían proporcionar restricciones directas sobre estos escenarios, y que resolver la tensión con las observaciones de sGRB podría requerir reevaluar los canales de progenitores de sGRB o los modelos subyacentes de síntesis de poblaciones. El trabajo no afirma resolver definitivamente el problema del origen del proceso-r, sino más bien delinear los requisitos evolutivos específicos y las tensiones observacionales que cualquier modelo exitoso basado únicamente en BNS debe abordar.