Implications of low neutron star merger rates for gamma-ray bursts, r-process production and Galactic double neutron stars

Este estudio analiza la tasa actualizada de fusiones de estrellas de neutrones binarias derivada de los datos de GWTC-4 y revela una tensión creciente entre esta tasa y las estimaciones independientes basadas en estallidos de rayos gamma cortos, la producción de elementos del proceso-r y la población galáctica de sistemas de doble estrella de neutrones.

Lo esencial

Imagina que el universo es una inmensa fábrica de estrellas y explosiones, y nosotros, los astrónomos, somos los inspectores de calidad tratando de entender cuántos productos salen de esa fábrica.

Hace unos años, hicimos un descubrimiento histórico: detectamos por primera vez el "golpe" de dos estrellas de neutrones chocando (llamado GW170817). Fue como ver la primera vez que una fábrica de coches producía un modelo específico. Ese evento nos dijo dos cosas increíbles:

1. Esas colisiones crean los gamma-ray bursts (SGRBs), que son como los "fuegos artificiales" más brillantes del cosmos.
2. Son las "fábricas" que crean los elementos pesados, como el oro y el platino (lo que llamamos el proceso r).

Al principio, con solo ese primer evento, calculamos que la fábrica producía muchísimos de estos choques. Todo encajaba: el número de choques que veíamos coincidía con el número de fuegos artificiales que veíamos, con la cantidad de oro que hay en nuestra galaxia y con el número de parejas de estrellas de neutrones que tenemos en el vecindario.

Pero ahora, las cosas han cambiado.

La paper que me has pasado es como un informe de actualización de 2026. Los científicos han seguido escuchando al universo con sus "orejas" (los detectores de ondas gravitacionales LIGO y Virgo) durante más tiempo. Y aquí viene el giro de la trama: han encontrado muchos menos choques de los que esperaban.

Es como si, después de instalar cámaras de seguridad más potentes en la fábrica, te dieras cuenta de que, en realidad, solo salen dos camiones al día, y no los doscientos que pensabas.

El problema: La cuenta no cuadra

Los autores de este estudio dicen: "Oye, si solo hay 2 choques por cada 100 años en cada galaxia (una tasa muy baja), entonces deberíamos ver menos fuegos artificiales, menos oro y menos parejas de estrellas".

Pero cuando miramos el cielo, no es así.

• Los fuegos artificiales (SGRBs): Hay muchos más de los que deberíamos ver si solo hay dos choques. Es como si la fábrica produjera 2 camiones, pero en el cielo hubiera 10 fuegos artificiales. ¿De dónde salen los otros 8?
• El oro (Proceso r): Hay tanto oro en nuestra galaxia que, con tan pocos choques, parece que no habría habido tiempo suficiente para fabricarlo todo.
• Las parejas vecinas (DNS): En nuestra propia galaxia, contamos las parejas de estrellas de neutrones que sabemos que van a chocar en el futuro. El número que calculamos a partir de ellas es mucho más alto que el número de choques que hemos escuchado con los detectores.

¿Qué significa esto? (Las explicaciones posibles)

Los científicos proponen varias teorías para explicar por qué nuestra "cuenta" de la fábrica no coincide con la realidad:

1. Los fuegos artificiales son más grandes de lo que pensábamos:
   Quizás los choques no son como focos estrechos que solo iluminan un punto, sino como faros gigantes que iluminan todo el cielo. Si los "focos" (los chorros de energía) son más anchos, veríamos más fuegos artificiales aunque haya menos choques. Pero, hasta ahora, la mayoría de los que hemos visto parecen tener focos muy estrechos.

2. La fábrica ha cambiado con el tiempo:
   Tal vez la fábrica producía muchísimos más choques hace miles de millones de años (cuando el universo era joven) y ahora ha frenado. Los detectores de ondas gravitacionales solo escuchan lo que pasa "ahora" (cerca de nosotros), pero los fuegos artificiales y el oro acumulado son el resultado de toda la historia de la fábrica. Si la fábrica trabajaba a todo vapor antes, eso explicaría por qué hay tanto oro y tantos fuegos artificiales, aunque hoy veamos pocos choques.

3. Nuestros vecinos son una excepción:
   Quizás nuestra galaxia (la Vía Láctea) es especial y tiene más parejas de estrellas de neutrones que la galaxia promedio. Si es así, contar las parejas en nuestro vecindario nos daría un número alto, pero no representaría a todo el universo.

4. Estamos perdiendo algunos choques:
   Tal vez hay estrellas de neutrones muy pequeñas o muy ligeras que chocan, pero nuestros detectores son "sordos" a ellas porque son demasiado débiles. Si existiera una "fábrica fantasma" de estrellas pequeñas que no podemos ver, eso aumentaría el número total de choques y solucionaría el problema.

En resumen

Este estudio nos dice que el universo es más misterioso de lo que pensábamos. Tenemos una tasa de colisiones de estrellas de neutrones que parece demasiado baja para explicar todo lo que vemos: el oro, las explosiones de rayos gamma y las estrellas vecinas.

Es como si un detective llegara a una escena del crimen, contara solo dos balas en el suelo, pero viera que hay 500 agujeros en la pared. Algo no encaja. Los autores sugieren que necesitamos revisar nuestras suposiciones: ¿Son los fuegos artificiales más grandes? ¿La fábrica trabajaba más duro en el pasado? ¿O hay una parte de la fábrica que no podemos ver?

Esta investigación es crucial porque nos obliga a pensar de nuevo en cómo funcionan las estrellas, cómo se crea el oro que llevamos en nuestros anillos y cómo evoluciona el universo a lo largo del tiempo.

Resumen técnico

Aquí presento un resumen técnico detallado del artículo "Implications of low neutron star merger rates for gamma-ray bursts, r-process production and Galactic double neutron stars" (Implicaciones de las bajas tasas de fusión de estrellas de neutrones para los estallidos de rayos gamma, la producción de proceso-r y las estrellas de neutrones dobles galácticas), basado en la versión borrador del 8 de abril de 2026.

1. El Problema

Desde el descubrimiento multimensajero de GW170817, se estableció que las fusiones de estrellas de neutrones binarias (BNS) son la fuente de al menos algunos estallidos de rayos gamma cortos (SGRB) y la producción de elementos del proceso-r. Inicialmente, la tasa de fusión de BNS inferida de los primeros datos de LIGO-Virgo (110–3840 Gpc⁻³ yr⁻¹) era consistente con las tasas observadas de SGRB, la masa de proceso-r en la Vía Láctea y la población de sistemas de estrellas de neutrones dobles (DNS) galácticos.

Sin embargo, tras la actualización del catálogo de ondas gravitacionales GWTC-4 (que cubre los primeros 8 meses de la cuarta campaña de observación, O4), solo se han detectado dos eventos significativos de BNS (GW170817 y GW190425) y ninguno nuevo en O4. Esto ha llevado a una revisión a la baja de la tasa de fusión de BNS. El problema central de este trabajo es investigar si esta tasa de fusión de BNS reducida (inferida de las ondas gravitacionales) es consistente con las tasas independientes derivadas de:

1. La tasa cosmológica de SGRB.
2. La masa total de elementos del proceso-r en la Vía Láctea.
3. La población observada de sistemas DNS en nuestra galaxia.

2. Metodología

Los autores emplean un enfoque comparativo riguroso utilizando los datos más recientes de GWTC-4:

• Inferencia de Tasas de GW: En lugar de promediar sobre distribuciones de masa inciertas, los autores calculan las tasas de fusión en binarios de masa específicos (histogramas por tramos).
  • Dividen la distribución de masa en tres bins centrados en ~1.3 M☉, ~1.65 M☉ y ~2.1 M☉.
  • Incluyen límites superiores conservadores para fusiones de estrellas de neutrones y agujeros negros (NSBH) que podrían ser brillantes en el espectro electromagnético (EM-bright), considerando la masa del agujero negro y el espín.
  • Utilizan priores conservadores (Jeffreys para bins sin detecciones) y asumen distribuciones uniformes de masa y espín dentro de los bins.
• Comparación con SGRB: Comparan la tasa de BNS con la tasa de SGRB cosmológica, analizando la relación entre la tasa volumétrica y el ángulo de apertura del chorro ($\theta_{jet}$).
• Cálculo de Proceso-r: Recalculan la masa total de elementos del proceso-r en la Vía Láctea, incorporando por primera vez de manera explícita el medio circumgaláctico (CGM) y refinando la fracción de masa basada en isótopos pesados ($Z \ge 56$) para evitar incertidumbres del proceso-s.
• Población DNS Galáctica: Estiman la tasa de fusión de la Vía Láctea basándose en los sistemas DNS observados, aplicando correcciones por efectos de selección de radio y beaming de púlsares.

3. Contribuciones Clave y Resultados

A. Tasa de Fusión de BNS Actualizada (GWTC-4)

• La tasa total de fusión de BNS inferida es de $110^{+192}_{-82}$ Gpc⁻³ yr⁻¹ (90% de credibilidad).
• Desglosada por masa:
  • Eventos tipo GW170817 (~1.3 M☉): $53^{+176}_{-49}$ Gpc⁻³ yr⁻¹.
  • Eventos tipo GW190425 (~1.65 M☉): $30^{+99}_{-28}$ Gpc⁻³ yr⁻¹.
  • Eventos de alta masa (>1.6 M☉): Límite superior < 30 Gpc⁻³ yr⁻¹.
• Si no se detectan nuevos eventos al final de O4, la tasa podría reducirse a ~16–170 Gpc⁻³ yr⁻¹.

B. Tensión con Estallidos de Rayos Gamma Cortos (SGRB)

• Existe una tensión significativa. La tasa cosmológica de SGRB es un factor de 3.6 a 18 veces mayor que la tasa central de BNS de 100 Gpc⁻³ yr⁻¹.
• Para reconciliar las tasas, se requerirían ángulos de apertura de chorro muy anchos ($\theta_{jet} \gtrsim 10^\circ$) para la mayoría de los SGRB, lo cual es inconsistente con las observaciones de la mayoría de los jets (que suelen ser más estrechos).
• Esto sugiere que o bien la tasa de SGRB local está sobreestimada, o bien una gran fracción de los SGRB provienen de un canal de progenitores diferente (no fusiones de BNS), o bien existe una evolución de la tasa con el corrimiento al rojo (redshift) que no se está capturando.

C. Tensión con la Producción de Proceso-r

• Los autores recalculan la masa de proceso-r en la Vía Láctea, incluyendo el CGM, obteniendo un valor de $3800 \pm 800 M_\odot$ (aproximadamente el doble de estimaciones anteriores).
• Para producir esta masa, se requiere una tasa de fusión de NSM (BNS + NSBH) de 89–410 Gpc⁻³ yr⁻¹.
• Aunque la tasa de GW (110 Gpc⁻³ yr⁻¹) cae dentro del rango de incertidumbre, está en el extremo superior. Si la tasa de GW disminuye aún más, se requeriría una producción de material por evento mayor (>0.01 M☉) o una historia de formación estelar de la Vía Láctea con fusiones más tempranas (canales "prompt").

D. Tensión con Estrellas de Neutrones Dobles (DNS) Galácticas

• La tasa inferida de la población DNS galáctica es de 230–510 Gpc⁻³ yr⁻¹ (o 23–51 Myr⁻¹ por galaxia).
• Esto es un factor de 2.3 a 5.1 mayor que la tasa de GW de 100 Gpc⁻³ yr⁻¹.
• La tensión es particularmente aguda si se asume que solo las fusiones de baja masa (tipo GW170817) representan a los sistemas DNS observados. Esto podría indicar que la Vía Láctea tiene una tasa de fusión inusualmente alta o que las correcciones por beaming de púlsares están sobreestimadas.

4. Significado e Implicaciones Físicas

El trabajo destaca una tensión emergente entre las tasas de fusión medidas directamente por ondas gravitacionales (bajas) y las tasas inferidas de fenómenos electromagnéticos y químicos (altas). Las implicaciones físicas para resolver esta tensión incluyen:

1. Propiedades de los Jets: Es probable que los jets de los SGRB sean más anchos de lo pensado o que existan mecanismos de "asfixia" (choked jets) donde la fusión ocurre pero no produce un estallido observable.
2. Evolución de la Tasa (Redshift): Podría existir una fuerte evolución de la tasa de fusión con el corrimiento al rojo. Las fusiones con tiempos de retardo cortos (prompt) serían más comunes en el pasado ($z > 0$), lo que elevaría la tasa cosmológica de SGRB y de proceso-r, mientras que la tasa local ($z \approx 0$) medida por GW sería menor.
3. Distribución de Masas y Espines:
   • Si las estrellas de neutrones tienen espines bajos (<0.05), la tasa real de GW podría ser un factor de ~2 menor, exacerbando la tensión.
   • Si existe una población oculta de estrellas de neutrones de muy baja masa (<1 M☉), la tasa podría ser un factor de 2-3 mayor, aliviando la tensión.
4. Progenitores Alternativos: Es posible que una fracción significativa de los SGRB provenga de otros canales (ej. colapsares) o que las fusiones NSBH brillantes en EM sean más raras de lo que se pensaba.

Conclusión:
El artículo concluye que la tasa de fusión de BNS inferida por GWTC-4 es consistente con las estimaciones anteriores pero más baja, creando una tensión estadística con las observaciones independientes. Esta discrepancia no es necesariamente un error, sino una oportunidad para restringir parámetros astrofísicos fundamentales: la estructura de los jets de GRB, la historia de formación estelar de la Vía Láctea, la distribución de espines de las estrellas de neutrones y la naturaleza de los progenitores de los estallidos de rayos gamma. Se requiere una mayor precisión en las observaciones de GW (post-O4) y en las contrapartes electromagnéticas para resolver estas tensiones.