Inferring the role of binary neutron star mergers in r-process nucleosynthesis with multi-messenger observations using Cosmic Explorer and Einstein Telescope

Este trabajo propone un método que utiliza detectores de ondas gravitacionales de tercera generación (Cosmic Explorer y Einstein Telescope) para medir la contribución fraccional de las fusiones de estrellas de neutrones binarias a la nucleosíntesis de procesos-r cósmicos mediante el análisis de la correlación dependiente del corrimiento al rojo entre las tasas de eventos de fusión y las abundancias de procesos-r, logrando una precisión de aproximadamente 5–6% tanto para los escenarios de "sirenas brillantes" como de "sirenas oscuras" de mensajería múltiple.

Lo esencial

El Gran Misterio: ¿De dónde provienen los elementos pesados?

Imagina el universo como una cocina gigante. Los "ingredientes" para los elementos pesados en nuestros cuerpos (como el oro, el platino y el uranio) se crean en un proceso llamado proceso-r. Los científicos han sabido durante mucho tiempo que esto ocurre durante eventos cósmicos violentos, pero aún discuten sobre qué evento es el chef principal.

Hay dos principales sospechosos:

1. Fusiones de Estrellas de Neutrones Binarias (BNS): Dos estrellas muertas y súper densas chocando entre sí.
2. Supernovas Raras: Estrellas que explotan girando increíblemente rápido o colapsando en agujeros negros.

Sabemos que el primer sospechoso (las fusiones de estrellas de neutrones) puede cocinar estos elementos porque lo vimos ocurrir una vez (el famoso evento GW170817). Sin embargo, no estamos seguros si ellos hacen todo el trabajo de cocina, o solo una pequeña parte. Algunas evidencias sugieren que podrían ser demasiado lentos para explicar por qué los elementos pesados existían en el universo muy temprano.

La Nueva Idea: Una verificación de "recibo" cósmico

Los autores de este artículo proponen una nueva forma de resolver este misterio. Sugieren que dejemos de mirar solo un evento y empecemos a observar la historia completa del universo como un gran libro de contabilidad.

Piénsalo así:

• El Detector de Ondas Gravitacionales (GW) es un contador. Cuenta cuántas veces chocan dos estrellas de neutrones en todo el universo.
• El Telescopio es un espectroscopista. Mide cuánto "sopa de elementos pesados" (abundancia del proceso-r) existe en las galaxias en diferentes momentos del pasado.

El método de los autores es comparar estas dos listas. Si las fusiones de estrellas de neutrones son la única fuente de elementos pesados, el número de choques y la cantidad de elementos pesados deberían subir y bajar en perfecta sincronía. Si no son la única fuente, las dos listas se desviarán entre sí.

La Máquina del Tiempo: Mirando hacia atrás en el tiempo

Para hacer esto, los científicos necesitan mirar muy atrás en el tiempo. Proponen usar futuras "máquinas del tiempo" superpoderosas (telescopios y detectores) llamadas Cosmic Explorer y el Telescopio Einstein.

Estas máquinas podrán:

1. Contar los choques: Detectar miles de colisiones de estrellas de neutrones de hace miles de millones de años (alto desplazamiento al rojo).
2. Probar la sopa: Medir los elementos pesados en galaxias que existían hace miles de millones de años.

El artículo simula qué pasaría si usáramos estas máquinas durante un año de observación. Crearon una base de datos "ficticia" (un universo falso) para probar sus matemáticas.

Los Dos Escenarios: Sirenas Brillantes vs. Oscuras

El artículo prueba dos formas diferentes de recopilar datos, usando una analogía de "sirena":

1. La "Sirena Brillante" (El Caso Ideal):

   • Imagina que ocurre un choque y también emite un destello de luz (como una kilonova o un estallido de rayos gamma) que nuestros telescopios pueden ver.
   • Esta luz nos dice exactamente dónde ocurrió el choque y qué tan lejos está. Es como ver un choque de coches y la matrícula claramente.
   • Resultado: Esto proporciona datos muy precisos.

2. La "Sirena Oscura" (El Caso Más Difícil):

   • Imagina que ocurre un choque, pero no hay destello de luz. Solo escuchamos el "sonido" (ondas gravitacionales) pero no podemos ver la fuente.
   • Tenemos que adivinar la distancia basándonos solo en el sonido, lo cual es más difuso. Es como escuchar un choque en la oscuridad y adivinar dónde ocurrió.
   • Resultado: Esto es menos preciso, pero el artículo muestra que aún funciona lo suficientemente bien.

¿Qué descubrieron?

Usando su verificación de "recibo" matemática, los autores encontraron que:

• Precisión: Incluso en el escenario de "Sirena Oscura" (sin luz), pudieron determinar cuánto de los elementos pesados proviene de las fusiones de estrellas de neutrones con una precisión de aproximadamente 94–95% (lo que significa un margen de error del 5–6%).
• El Factor "Retraso": También pudieron calcular cuánto tiempo tardan las estrellas de neutrones en fusionarse después de nacer. ¿Chocan inmediatamente o esperan miles de millones de años? Su método puede medir este "tiempo de espera" bastante bien, aunque es ligeramente más difícil que medir la cantidad total de elementos.
• El Veredicto: Si las fusiones de estrellas de neutrones son responsables de una parte significativa (más del 10%) de los elementos pesados del universo, este método puede demostrarlo.

La Conclusión

Este artículo no afirma haber resuelto el misterio aún porque todavía no tenemos estos super-telescopios construidos. En cambio, es un plan.

Dice: "Si construimos estos detectores de próxima generación y comenzamos a medir elementos pesados en galaxias distantes, finalmente podremos contar exactamente cuánto del oro y el uranio del universo proviene de estrellas de neutrones chocando versus estrellas explotando."

Convierte la pregunta de "¿Quién es el chef?" en un problema matemático que realmente podemos resolver comparando el número de choques con la cantidad de comida en la mesa.

Resumen técnico

Resumen Técnico: Inferencia del Papel de las Fusiones de Estrellas de Neutrones Binarias en la Nucleosíntesis del Proceso-r con Observaciones Multimensajero

Planteamiento del Problema
El origen cósmico de los elementos de captura rápida de neutrones (proceso-r) sigue siendo una pregunta fundamental abierta en astrofísica. Aunque las fusiones de estrellas de neutrones binarias (BNS) están confirmadas como sitios del proceso-r (por ejemplo, GW170817), su dominio exclusivo es cuestionado por estudios de evolución química. Estos estudios destacan discrepancias respecto a los tiempos de retraso de las fusiones de BNS en relación con la formación estelar, las tasas observadas de BNS por los detectores actuales de ondas gravitacionales (GW) y los problemas de retención en halos de baja masa. Por el contrario, clases raras de supernovas de colapso del núcleo (CCSNe), como las supernovas magnetorotacionales y los colapsares, ofrecen un enriquecimiento "rápido" pero carecen de firmas observacionales definitivas de producción del proceso-r. Los métodos actuales que dependen de la arqueología estelar local o de modelos de evolución química enfrentan degeneraciones entre la historia de formación estelar, las tasas de fusión y las suposiciones sobre el rendimiento.

Metodología
Los autores proponen una técnica de correlación novedosa para medir la contribución acumulada fraccional de las fusiones de BNS ($F_{\text{BNS},z_0}$) a la abundancia cósmica total del proceso-r. El método explota la correlación dependiente del corrimiento al rojo entre el número total de eventos GW de BNS observados ($N_{\text{GW}}$) y las abundancias promedio del proceso-r ($Ab_{\text{El}}$) de galaxias en corrimientos al rojo $z \lesssim 1$.

1. Marco Teórico: El método asume que tanto $N_{\text{GW}}(z)$ como $Ab_{\text{El}}(z)$ son funciones de la historia de formación estelar cósmica ($\psi_{\text{SF}}$). Al definir una función de correlación $K(\theta) = N_{\text{GW}}(z) / Ab_{\text{El}}(z)$, se elimina la dependencia de la historia de formación estelar, haciendo que la correlación sea sensible principalmente a los parámetros de la distribución de tiempos de retraso (DTD) de las BNS (retraso mínimo $t_{\text{min}}$ e índice de ley de potencias $b$) y a la contribución fraccional $F_{\text{BNS},z_0}$.
2. Modelado:
   • Eventos GW: La tasa de fusión de BNS se modela como una convolución de la tasa de formación estelar con una DTD de ley de potencias ($D(t) \propto t^b$ para $t > t_{\text{min}}$).
   • Abundancias: Se emplea un modelo de evolución química cósmica (CCE), asumiendo que el enriquecimiento del proceso-r surge de dos fuentes: fusiones de BNS (retrasadas) y CCSNe del proceso-r "rápidas" (rCCSNe).
   • Escenarios: El análisis considera dos escenarios observacionales límite para detectores GW de tercera generación (3G) (Explorador Cósmico y Telescopio Einstein):
     • Sirenas Brillantes: Eventos GW con contrapartes electromagnéticas (EM), que permiten una determinación precisa del corrimiento al rojo de la galaxia anfitriona.
     • Sirenas Oscuras: Eventos GW sin contrapartes EM, que dependen de la conversión de distancia de luminosidad a corrimiento al rojo con mayores incertidumbres.
3. Análisis Estadístico: Los autores emplean proyecciones de la matriz de información de Fisher para estimar la precisión de la recuperación de parámetros. Generan datos simulados para $N_{\text{GW}}$ y $Ab_{\text{Eu}}$ (utilizando Europio como trazador) incorporando incertidumbres observacionales esperadas, incluidas estadísticas de conteo, errores de corrimiento al rojo y dispersión intrínseca en las mediciones de abundancia.

Contribuciones y Resultados Clave
Utilizando proyecciones de Fisher con un escenario fiducial de dos detectores Explorador Cósmico (CE) y un Telescopio Einstein (ET) observando durante un año, el artículo demuestra lo siguiente:

• Precisión en $F_{\text{BNS},z_0}$: El método puede estimar la contribución acumulada fraccional de las fusiones de BNS al proceso-r cósmico con una precisión de $\lesssim 5-6\%$ a $1\sigma$ para escenarios fiduciales donde $F_{\text{BNS},z_0} \gtrsim 0.1$. Esta precisión se mantiene tanto para escenarios de sirenas brillantes como de sirenas oscuras, ya que la estimación está dominada por la incertidumbre en el proxy de abundancia y no por la incertidumbre del corrimiento al rojo de las GW.
• Parámetros DTD: El método produce estimaciones para el índice de ley de potencias de la DTD, $b$, comparables a las restricciones existentes de estudios de estallidos cortos de rayos gamma (GRB) (precisión $\lesssim 10-20\%$). Sin embargo, el tiempo de retraso mínimo $t_{\text{min}}$ está menos restringido (incertidumbre $\approx 100-150\%$) en el valor fiducial de $t_{\text{min}} = 0.2$ Gyr, debido a la baja sensibilidad del tiempo de retraso promedio a $t_{\text{min}}$ cuando $t_{\text{min}} \ll t_{\text{max}}$.
• Requisitos de Corrimiento al Rojo: La relación señal-ruido (SNR) para restringir $F_{\text{BNS},z_0}$ dentro de $\lesssim 20\%$ requiere observaciones de abundancias del proceso-r y de GW hasta $z \lesssim 0.7$. Las restricciones sobre los parámetros DTD $b$ y $t_{\text{min}}$ mejoran con observaciones que se extienden hasta $z \lesssim 1.0$.
• Independencia de la Historia de Formación Estelar: Se demuestra que la técnica de correlación es en gran medida independiente de la historia de formación estelar cósmica subyacente, aislando la física de las fuentes del proceso-r.

Significado y Afirmaciones
Los autores afirman que este método proporciona una manera robusta y libre de modelos de cuantificar el papel de las fusiones de BNS en la nucleosíntesis cósmica del proceso-r utilizando futuros datos GW de 3G y electromagnéticos.

• Resolución de Degeneraciones: A diferencia de los estudios previos de evolución química que luchan con degeneraciones entre tasas de fusión, rendimientos e historias de formación estelar, este método de correlación vincula directamente la tasa de eventos con el rendimiento químico.
• Caso Científico para Abundancias de Alto Corrimiento al Rojo: El artículo argumenta que, aunque la arqueología estelar local puede reconstruir tendencias de abundancia, este método proporciona un caso científico convincente para identificar firmas de elementos de captura de neutrones en galaxias más allá del Universo local ($z \gtrsim 0.1$).
• Implicaciones Más Amplias: La determinación precisa de los parámetros DTD de las BNS mediante este método mejoraría la identificación de galaxias anfitrionas para el seguimiento multimensajero, refinaría las restricciones sobre parámetros cosmológicos mediante sirenas estándar e informaría modelos de canales de formación de BNS (por ejemplo, fases de envoltura común).

Los autores permanecen modestos respecto a la viabilidad de obtener los datos de abundancia necesarios a altos corrimientos al rojo, señalando que, aunque el progreso en espectroscopía de continuo y la inferencia de abundancias en fase gaseosa (por ejemplo, mediante JWST) es prometedor, la realización de estas mediciones sigue siendo un desafío abierto. El método se presenta como una herramienta para aplicar una vez que tales datos estén disponibles, o como una motivación para desarrollar las capacidades observacionales necesarias.