Forming Double Neutron Stars using Detailed Binary Evolution Models with POSYDON: Comparison to the Galactic Systems

Utilizando el código de síntesis de población binaria POSYDON, este estudio revela que las estrellas de neutrones dobles galácticas se forman a través de dos subcanales distintos de envoltura común (Caso B y Caso C) que requieren definiciones de núcleo específicas o eficiencias de eyección altas y patadas de supernova bajas, explicando así las divisiones observadas en los periodos orbitales y restringiendo la física evolutiva detallada de estos sistemas.

Lo esencial

La visión general: Una historia de detectives cósmicos

Imagina el universo como una pista de baile gigante y caótica. La mayor parte del tiempo, las estrellas bailan solas. Pero a veces, dos estrellas masivas nacen juntas, tomadas de la mano. A lo largo de millones de años, pasan por un baile salvaje y dramático que puede terminar de dos maneras: o se separan para siempre, o chocan entre sí y se fusionan.

Este artículo es una historia de detectives sobre las Estrellas de Neutrones Dobles (DNS). Estas son los "pesos pesados" de la pista de baile cósmica: dos estrellas increíblemente densas y muertas (estrellas de neutrones) orbitando entre sí. Conocemos unas 25 de estas parejas en nuestra propia galaxia, la Vía Láctea. Los autores utilizaron una simulación de supercomputadora llamada POSYDON para averiguar exactamente cómo se forman estos pares y por qué lucen de la forma en que lo hacen.

El principal descubrimiento: Dos caminos diferentes hacia el mismo destino

La mayor sorpresa de este artículo es que no hay una sola forma de crear una Estrella de Neutrones Doble. Los autores descubrieron que estos pares se forman a través de dos "subcanales" distintos, como dos autopistas diferentes que conducen a la misma ciudad.

La diferencia depende de la "edad" y la "madurez" de la segunda estrella cuando comienza el drama:

1. El camino "Joven" (El canal del núcleo de helio):

   • El escenario: Imagina una estrella joven que no ha tenido mucho tiempo para crecer. Todavía tiene un núcleo hecho de helio (como un adolescente que aún no ha terminado la secundaria).
   • El drama: Cuando esta estrella joven intenta expandirse, se acerca demasiado a su compañera (la primera estrella de neutrones). Se enredan en una nube gigante y desordenada de gas llamada Envoltura Común (Common Envelope).
   • El resultado: Para sobrevivir a este enredo, el sistema necesita mucha energía para expulsar la nube de gas. Si lo logran, terminan en una órbita muy apretada y rápida. Estos son los pares que eventualmente chocarán entre sí y crearán ondas gravitacionales (los que se "fusionan").

2. El camino "Viejo" (El canal del núcleo de carbono-oxígeno):

   • El escenario: Ahora imagina una estrella más vieja, más madura. Ha tenido tiempo de quemar su combustible y tiene un núcleo hecho de carbono y oxígeno (como un jubilado que ha terminado su carrera).
   • El drama: Esta estrella más vieja también se enreda en una nube de gas con su compañera.
   • El resultado: Debido a que la nube de gas de la estrella más vieja es más laxa y fácil de empujar, el sistema sobrevive al enredo mucho más fácilmente. Sin embargo, terminan en una órbita mucho más amplia y lenta. Estos pares están demasiado lejos para chocar entre sí dentro de la vida útil del universo.

La analogía: Piensa en ello como dos parejas tratando de detener una pelea en una habitación llena de gente.

• El Pareja Joven (núcleo de helio) está peleando en un espacio estrecho y apretado. Se requiere mucho esfuerzo (energía) para apartar a la multitud, pero si lo logran, terminan abrazándose muy fuerte.
• La Pareja Vieja (núcleo de carbono-oxígeno) está peleando en una habitación espaciosa. Es más fácil apartar a la multitud, pero terminan parados lejos el uno del otro.

El problema del "Patada" (Kick): Por qué la pista de baile es resbaladiza

Cuando una estrella muere, explota como una supernova. Esta explosión a menudo le da a la nueva estrella de neutrones una "patada", como una bola de billar siendo golpeada.

• La teoría antigua: Los científicos solían pensar que estas patadas eran enormes y aleatorias (como una bola de billar siendo golpeada con un mazo).
• El nuevo hallazgo: Los autores descubrieron que si las patadas fueran tan grandes, la mayoría de los pares saldrían volando y nunca se convertirían en Estrellas de Neutrones Dobles.
• La solución: Para coincidir con lo que realmente vemos en el cielo, las "patadas" deben ser mucho más suaves (como un toque ligero). El artículo sugiere que la segunda estrella en el par recibe una patada muy suave porque ha sido despojada de la mayor parte de su peso antes de explotar. Esto permite que el par permanezca unido.

El "Reciclaje" de la primera estrella

Antes de que la segunda estrella explote, la primera estrella de neutrones (la más vieja) recibe un "cambio de imagen". A medida que la segunda estrella desprende su gas, la primera estrella se lo come. Esto es como un plan de dieta cósmica que, de hecho, hace que la estrella gire más rápido y brille más, convirtiéndola en un "púlsar reciclado".

El artículo muestra que el "Camino Joven" (núcleo de helio) proporciona una mejor dieta para la primera estrella, acelerando su rotación. Esto explica por qué los pares que se fusionan que vemos en el cielo (Subpoblación i) giran muy rápido, mientras que los pares que no se fusionan (Subpoblación ii) giran más lento.

¿Qué pasa con los extraños?

El artículo señala que hay un pequeño grupo de tres Estrellas de Neutrones Dobles que no encajan en ninguno de estos dos caminos. Tienen órbitas extrañas que son difíciles de explicar. Los autores admiten que su modelo no puede explicar fácilmente estos tres, sugiriendo que podrían haberse formado de una manera totalmente distinta (quizás chocando entre sí en un cúmulo estelar concurrido, en lugar de nacer como un par).

Resumen de resultados

• La división: Las Estrellas de Neutrones Dobles se forman a través de dos caminos principales: uno para donantes "jóvenes" (que conduce a pares que se fusionan) y otro para donantes "viejos" (que conduce a pares que no se fusionan).
• La energía: Para que el "camino joven" funcione, el universo necesita ser muy eficiente en la expulsión de las nubes de gas.
• Las patadas: Las explosiones que crean estas estrellas deben ser suaves, no violentas, para mantener los pares unidos.
• Las fusiones: Los pares que eventualmente chocarán y crearán ondas gravitacionales (como el famoso evento GW170817) provienen casi con seguridad del "Camino Joven".

Lo que el artículo NO dice

• No predice cuándo ocurrirá la próxima fusión.
• No afirma explicar cada estrella extraña en el universo (específicamente las tres "extrañas" mencionadas anteriormente).
• No incluye un modelo de cómo las estrellas se "frenan" o se desvanecen con el tiempo, lo que limita qué tan perfectamente pueden coincidir con las observaciones del mundo real en este momento.

En resumen, este artículo utiliza una simulación computacional detallada para mostrar que el universo tiene dos recetas distintas para fabricar Estrellas de Neutrones Dobles, y solo una de esas recetas produce los pares que eventualmente colisionarán y sacudirán el tejido del espacio-tiempo.

Resumen técnico

Resumen Técnico: Formación de Estrellas de Neutrones Dobles mediante Modelos Detallados de Evolución Binaria con POSYDON

Planteamiento del Problema
Las estrellas de neutrones dobles (DNS, por sus siglas en inglés) sirven como laboratorios críticos para comprender la evolución de binarias masivas, los mecanismos de supernovas y las fuentes de ondas gravitacionales. Aunque se han identificado más de dos docenas de DNS en la Vía Láctea, sus canales de formación siguen siendo complejos y objeto de debate. Los estudios previos de síntesis de población se han basado mayoritariamente en códigos de síntesis de población binaria (BPS) rápidos que utilizan fórmulas de ajuste o librerías de trayectorias de estrellas individuales. Estos enfoques suelen carecer de un tratamiento consistente de las interacciones binarias, particularmente durante la fase de envolvente común (CE) y las supernovas de colapso de núcleo (SN). Además, la población observada de DNS en la Galaxia exhibe una bifurcación distintiva en los períodos orbitales y las excentricidades que los modelos estándar luchan por reproducir sin un ajuste fino. Este artículo aborda la necesidad de restringir las etapas evolutivas detalladas que conducen a la formación de DNS utilizando un marco más riguroso y autoconsistente.

Metodología
Los autores utilizan POSYDON (POpulation SYnthesis with Detailed binary-evolution simulatiONs), un código de síntesis de población binaria de vanguardia que integra modelos detallados de evolución estelar y binaria basados en el código 1D MESA. A diferencia de los códigos BPS rápidos, POSYDON emplea mallas de evolución estelar precomputadas y cálculos en tiempo real para manejar la física compleja, incluyendo el frenado magnético, la radiación de ondas gravitacionales, la evolución de la envolvente común (CE) y los retrocesos (kicks) de supernova.

El estudio simula $2 \times 10^6$ poblaciones binarias a través de 25 variaciones distintas de modelos (resumidas en la Tabla 2). Las variables clave manipuladas incluyen:

1. Tratamientos de la Envolvente Común (CE): Variaciones en la eficiencia de eyección de la CE ($\alpha_{CE}$), la definición del límite núcleo-envolvente (fracciones de masa de hidrógeno de 1%, 10% y 30%) y el tratamiento de la pérdida de masa post-CE (transferencia de masa estable frente a pérdida por viento o fase única).
2. Modelos de Supernova (SN): Diferentes prescripciones para los resultados de colapso de núcleo (Sukhbold et al. 2016; Patton & Sukhbold 2020; Fryer et al. 2012), distinguiendo entre colapso de núcleo de hierro (CCSN) y supernovas de captura electrónica (ECSN).
3. Prescripciones de Retroceso Natal (Natal Kick): Comparación de la distribución estándar de Hobbs-Maxwellian ($\sigma \approx 265$ km/s) frente a modelos basados en la física:
   • AsymEj: Retrocesos derivados de mecanismos de eyección asimétrica (Janka 2017), que escalan con la masa de la eyecta.
   • Linear: Retrocesos que escalan linealmente con la masa de la eyecta (Richards et al. 2023).
   • Low Maxwellian: Un modelo de dispersión reducida ($\sigma \approx 20$ km/s).

Los modelos asumen metalicidad solar e una historia de formación estelar constante. Las poblaciones sintéticas resultantes se comparan cualitativamente con la muestra observada de 25 DNS galácticas (Tabla 1), centrándose en las distribuciones de período orbital ($P_{orb}$) y excentricidad ($e$), así como en las tasas de fusión.

Contribuciones Clave y Resultados

• Bifurcación de los Canales de Formación: El estudio identifica una bifurcación fundamental en el canal de formación por CE, determinada por el estado evolutivo de la estrella donante al inicio del desbordamiento del Lóbulo de Roche (RLO) con la primera estrella de neutrones nacida:

  • Canal de Núcleo de He (RLO de Caso B): La donante inicia la transferencia de masa mientras posee un núcleo de helio. Estos sistemas suelen tener períodos orbitales más cortos y energías de enlace de la envolvente más altas.
  • Canal de Núcleo de C/O (RLO de Caso C): La donante evoluciona más allá para desarrollar un núcleo de carbono-oxígeno antes del RLO. Estos sistemas poseen envolventes débilmente ligadas y resultan en períodos orbitales más amplios.
  • Significancia: Esta bifurcación explica naturalmente la división observada en la población de DNS galácticas entre sistemas de período corto que se fusionan (Subpoblación i) y sistemas de período amplio que no se fusionan (Subpoblación ii).

• Restricciones sobre la Física de la Envolvente Común (CE): La formación del canal de núcleo de He (que produce las DNS que se fusionan) es altamente sensible a la eficiencia de la CE. Los modelos revelan que este canal solo sobrevive si:

  • Se utiliza una definición de núcleo generosa (fracción de H $\simeq 30\%$), o
  • Se adopta una alta eficiencia de eyección de la CE ($\alpha_{CE} \gtrsim 1.2$).
    Sin estas condiciones, los sistemas de núcleo de He se fusionan dentro de la CE. El canal de núcleo de C/O es más robusto y sobrevive incluso con valores de $\alpha_{CE}$ más bajos.

• Restricciones de Retroceso Natal (Natal Kick): La distribución de retroceso estándar de Hobbs-Maxwellian ($\sigma \approx 265$ km/s) queda descartada, ya que desintegra la mayoría de las binarias durante la primera supernova, resultando en tasas de fusión órdenes de magnitud menores a las observadas. La población de DNS galácticas observada se reproduce mejor mediante prescripciones de retroceso de baja velocidad ($\lesssim 50$ km/s). Los modelos "basados en la física" AsymEj y Linear, que producen naturalmente retrocesos más bajos para progenitores despojados (stripped), reproducen con éxito tanto la distribución orbital como la tasa de densidad de fusión local ($O(10) - O(100)$ Gpc$^{-3}$ año$^{-1}$).

• Reciclaje de Púlsares y Acreción de Masa: El estudio encuentra que el reciclaje de púlsares (el giro rápido de la primera NS nacida) ocurre predominantemente durante la fase de transferencia de masa estable de Caso BB con la compañera de estrella de helio, en lugar de durante la fase de CE. Los sistemas formados a través del canal de núcleo de He experimentan un reciclaje más eficiente, lo cual es consistente con los períodos de rotación más cortos observados en la Subpoblación (i).

• Tipos de Supernova y Progenitores:

  • La primera supernova es predominantemente un colapso de núcleo de hierro (CCSN).
  • La segunda supernova a menudo involucra un progenitor ultra-despojado (ultra-stripped). Aproximadamente el 60% de los progenitores de DNS que se fusionan están desbordando sus lóbulos de Roche en el momento de la segunda SN.
  • Las DNS que se fusionan nacen de progenitores con masas de envolvente pre-SN medianas de $\sim 0.2 M_\odot$. Esto sugiere que una fracción significativa de las DNS que se fusionan podría ser observada como supernovas de Tipo Ib o Tipo Ic.

• Tasas de Fusión y Fuentes de GW: Los modelos predicen que solo las DNS formadas a través del canal de núcleo de He se fusionan dentro de un tiempo de Hubble. Por consiguiente, se infiere que eventos de ondas gravitacionales como GW170817 y GW190425 se originan de este canal. Las tasas de fusión locales predichas para los modelos con retrocesos basados en la física se alinean con las restricciones de LIGO-Virgo-Kagra y la muestra de DNS galácticas.

Significancia y Reivindicaciones
El artículo sostiene que los modelos evolutivos detallados son esenciales para resolver la historia de formación de las DNS. Al ir más allá de las fórmulas de ajuste hacia mallas autoconsistentes, los autores demuestran que la diversidad observada de las DNS galácticas es una consecuencia natural del estado evolutivo de la estrella donante al inicio de la fase de CE.

Los autores señalan modestamente que sus conclusiones se basan en comparaciones cualitativas debido a la ausencia actual de un modelo de evolución de púlsares dentro de POSYDON (que daría cuenta de los efectos de selección y el frenado de rotación). También reconocen que sus resultados son específicos para la metalicidad solar. Sin embargo, el estudio establece que:

1. Una bifurcación en los canales de formación (núcleo de He vs. núcleo de C/O) es una característica robusta de la evolución binaria que explica la división del período orbital observada.
2. Una eyección de CE eficiente (mediante definiciones de núcleo generosas o un $\alpha_{CE}$ alto) es necesaria para formar la población que se fusiona.
3. Los retrocesos natales bajos son una condición necesaria para la supervivencia de las binarias DNS, apoyando los cambios teóricos recientes que se alejan de las distribuciones de retroceso de alta velocidad.
4. La transferencia de masa inestable de Caso BB no produce DNS que se fusionen en estos modelos detallados, desafiando hipótesis previas sobre la formación de sistemas de fusión ultra-rápidos.

El trabajo proporciona una base para futuros análisis bayesianos cuantitativos una vez que se integren los modelos de evolución de púlsares en el marco de POSYDON.