Double Neutron Star Delay Times Across Cosmic Metallicities: The Role of Helium Star Progenitors

Resumen Técnico: Tiempos de Retraso de Fusiones de Estrellas de Neutrones Dobles a Través de Metallicitades Cósmicas

Enunciado del Problema
La distribución de tiempos de retraso (DTD) de las fusiones de estrellas de neutrones dobles (DNS) —el tiempo transcurrido entre la formación del sistema binario y la fusión— es un parámetro crítico para interpretar la contribución de las fusiones de DNS a la nucleosíntesis del proceso-r, a los estallidos cortos de rayos gamma (sGRB) y a las tasas de eventos de ondas gravitacionales. Aunque estudios anteriores han restringido la DTD utilizando datos observacionales de estrellas enriquecidas con elementos del proceso-r y galaxias anfitrionas de sGRB, la forma teórica subyacente de la DTD sigue siendo incierta. Específicamente, no está claro cómo la evolución dependiente de la metalicidad del sistema binario progenitor, particularmente la fase de estrella de helio (He), influye en la separación orbital al nacer y en las escalas de tiempo de fusión subsiguientes. Los modelos existentes a menudo asumen distribuciones de ley de potencia simplificadas para las DTD sin tener plenamente en cuenta la física compleja y dependiente de la metalicidad de la evolución de sistemas binarios aislados.

Metodología
Este trabajo emplea un enfoque guiado por la teoría que combina modelado detallado de evolución estelar con síntesis de poblaciones a gran escala para caracterizar la DTD de DNS a través de un rango cosmológico de metalicitades ($Z = [10^{-4}, 2] Z_{\odot}$).

1. Modelado de Evolución Estelar: Los autores utilizaron el código de evolución estelar 1-D MESA para modelar la evolución de estrellas de He individuales y binarias de estrella de He–NS. Los procesos físicos clave examinados incluyen:

   • Opacidad Radiativa: Analizar cómo la metalicidad afecta la opacidad de la envoltura y la expansión radial subsiguiente durante la secuencia principal de helio (HeMS).
   • Transferencia de Masa (MT): Modelar la fase de desbordamiento del lóbulo de Roche (RLO) (MT Caso BB) y el impacto de la pérdida de masa por viento en la evolución orbital. Las tasas de pérdida de masa por viento se modelaron utilizando la prescripción de Wolf–Rayet, con un ajuste específico para abordar las restricciones observacionales sobre estrellas despojadas (Y. Götberg et al. 2023).
   • Dinámica de Supernova (SN): Calcular las distribuciones orbitales post-SN asumiendo patadas natales isotrópicas de magnitudes variables para determinar cómo las patadas modifican la excentricidad y la separación.

2. Síntesis de Poblaciones: Los autores utilizaron el conjunto de evolución binaria POSYDON para simular la evolución completa de binarias masivas desde la Secuencia Principal de Edad Cero (ZAMS) hasta la fusión de DNS.

   • Mallas: Las simulaciones utilizaron mallas detalladas para binarias ricas en H, binarias de objeto compacto–ricas en H y binarias de objeto compacto–ricas en He a través de ocho metalicitades.
   • Variaciones Físicas: Se generaron cuatro modelos de población distintos para probar la robustez frente a incertidumbres en la física binaria, variando la eficiencia de eyección de la envoltura común (CE) ($\alpha_{CE}$) y las prescripciones de patadas natales (Eyecta Asimétrica vs. distribución Lognormal).
   • Corrección de Viento: Para abordar las sobreestimaciones en la pérdida de masa por viento para estrellas de He de baja masa en prescripciones estándar, los autores restringieron artificialmente la evolución de la fase estrella de He–NS en todas las metalicitades para seguir la trayectoria de la malla de $10^{-4} Z_{\odot}$, limitando efectivamente la expansión orbital no física.

Contribuciones y Resultados Clave

• Radio de la Estrella de He como Determinante: El estudio establece que el radio de la estrella de He durante la fase HeMS establece un límite inferior en la separación orbital de nacimiento de las DNS. La metalicidad modula este radio mediante la inflación de la envoltura impulsada por la opacidad; mayores metalicitades conducen a radios más grandes y, en consecuencia, a órbitas de nacimiento más amplias.
• Escalado del Tiempo de Fusión: Dado que el tiempo de espiral escala como $t_{insp} \propto a^4$, incluso diferencias modestas impulsadas por la metalicidad en los radios de las estrellas de He (por ejemplo, una diferencia de 5% entre $Z_{\odot}$ y $10^{-4} Z_{\odot}$) se traducen en variaciones significativas (20%) en los tiempos de fusión.
• Características Robustas de la DTD: A pesar de las variaciones en los supuestos de física binaria (eficiencia de CE, modelos de patadas), la síntesis de poblaciones revela dos características robustas en la DTD de DNS:
  1. Tiempo de Retraso Mínimo: Existe un límite inferior duro donde la mayor parte de las fusiones de DNS ocurren no antes de $\simeq 40$ Myr después de la formación estelar. El tiempo de formación absoluto mínimo (de ZAMS al nacimiento de DNS) es de $\simeq 20$ Myr, dictado por la escala de tiempo nuclear de la estrella primaria rica en H.
  2. Tiempo de Retraso Pico: La DTD alcanza un pico fuerte entre 80 y 250 Myr a través de todas las metalicitades.
• Estructura Dependiente de la Metalicidad: La forma de la DTD es compleja. En metalicitades intermedias ($Z \sim 10^{-2} Z_{\odot}$), el canal de formación de CE dominante se bifurca en dos subcanales (núcleo de He y núcleo de C/O), imprimiendo una estructura de doble pico en la DTD. En metalicitades bajas ($Z \lesssim 10^{-2} Z_{\odot}$), solo contribuye el canal de núcleo de C/O, mientras que en metalicitades solares y supersolares, solo el canal de núcleo de He produce DNS fusionadas.
• Patadas Natales: Si bien las orientaciones ideales de patadas pueden teóricamente producir tiempos de fusión muy cortos, los autores encuentran que las magnitudes de patada requeridas para producir la población observada de retrasos cortos son inconsistentes con las restricciones observacionales actuales. Los modelos de patadas típicos (por ejemplo, AsymEj) rara vez producen la contracción orbital extrema necesaria para retrasos inferiores a 20 Myr.
• Fracciones de Fusión:
  • $\approx 15\%$ de las DNS se fusionan dentro de 80 Myr, explicando potencialmente el enriquecimiento del proceso-r en entornos con historias de formación estelar breves.
  • $\gtrsim 20\%$ se fusionan en tiempos de retraso $> 1$ Gyr, proporcionando un mecanismo para sGRB en galaxias viejas y pobres en metales.
  • En metalicitades bajas ($Z \lesssim 10^{-2} Z_{\odot}$), $>70\%$ de las DNS se fusionan dentro de un tiempo de Hubble.

Significado y Afirmaciones
Los autores afirman que su trabajo proporciona un vínculo detallado y motivado físicamente entre los tiempos de retraso de las DNS y la evolución binaria de sus sistemas progenitores. Al demostrar que la DTD no es una simple ley de potencia, sino que posee una estructura compleja y dependiente de la metalicidad con un pico robusto entre 80–250 Myr, el estudio ofrece un marco refinado para evaluar la contribución de las fusiones de DNS al enriquecimiento cósmico del proceso-r y a eventos transitorios.

El artículo establece explícitamente que estos modelos ofrecen una base para futuros estudios que convuelvan las DTD derivadas con historias realistas de formación estelar cósmica para predecir tasas de fusión NS-NS, tasas de sGRB y transitorios de kilonova a través del tiempo cósmico. Los autores advierten que, aunque sus modelos sugieren que las fusiones de DNS pueden contribuir plausible al enriquecimiento en entornos de baja metalicidad, esto no es una prueba definitiva sin convolucionar con historias de formación estelar específicas. También señalan limitaciones respecto a la precisión del modelado de vientos de estrellas de He de baja masa y la necesidad de pruebas adicionales de las incertidumbres en la física binaria, manteniendo una postura modesta sobre la magnitud exacta de las variaciones en la forma de la DTD debidas a estos factores.