Exploring the diversity of kilonovae with 3D radiative transfer I. The polar direction

Resumen Técnico: Explorando la diversidad de las kilonovas mediante transferencia radiativa 3D

Título original: Exploring the diversity of kilonovae with 3D radiative transfer. I. The polar direction

1. El Problema (Contexto y Motivación)

Tras la detección de la señal de ondas gravitacionales GW170817 y la kilonova asociada AT2017gfo, la comunidad astrofísica ha intentado interpretar sus propiedades espectrales para comprender la nucleosíntesis de elementos pesados (proceso-$r$). Sin embargo, existen desafíos significativos:

• Incertidumbre en los datos atómicos: Los datos de los elementos del proceso-$r$ son incompletos o teóricos, lo que impide una identificación precisa de las líneas espectrales.
• Complejidad de la estructura del eyecta: Las kilonovas no son esféricas; la composición y la velocidad del material eyectado varían según la dirección de observación (polar vs. ecuatorial).
• Limitaciones de modelos previos: Muchos estudios se han limitado a modelos 1D o a un único modelo de fusión de estrellas de neutrones, lo que impide entender cómo la masa de las estrellas o la ecuación de estado (EoS) afectan la señal observada.

2. Metodología

El estudio emplea un enfoque de modelado avanzado y multidimensional:

• Transferencia Radiativa 3D: Utilizan el código de Monte Carlo artis para realizar simulaciones de transferencia radiativa dependientes del tiempo. El modelo asume equilibrio termodinámico local (LTE).
• Datos Atómicos Calibrados: Una contribución metodológica clave es la integración de nuevos datos atómicos calibrados para lantánidos (monocargados y bicargados) de Flörs et al. (2026), reemplazando datos teóricos menos precisos.
• Simulaciones Hidrodinámicas: Los modelos de transferencia radiativa se basan en simulaciones de hidrodinámica de partículas suavizadas (SPH) de fusiones de estrellas de neutrones. Se exploran diferentes configuraciones:
  • Ecuaciones de Estado (EoS): SFHO (relativamente "suave") y DD2 (relativamente "rígida").
  • Masas: Variaciones en la masa de las estrellas individuales y en la masa total del sistema, incluyendo sistemas con masas desiguales.
• Enfoque en la dirección polar: Se centran en la dirección $\theta < 37^\circ$, que es la orientación aproximada desde la cual se observó AT2017gfo.

3. Contribuciones Clave

• Ampliación del espacio de parámetros: Es uno de los primeros estudios que conecta directamente la diversidad de configuraciones de fusión (masas y EoS) con las propiedades espectrales y de luminosidad mediante transferencia radiativa 3D.
• Identificación espectral precisa: Gracias a los nuevos datos atómicos, el estudio logra identificar con mayor confianza las contribuciones de iones específicos en las características observadas de AT2017gfo.
• Análisis del eyecta dinámico: Demuestran que el eyecta dinámico por sí solo puede reproducir gran parte de la evolución espectral temprana, aunque con desfases temporales.

4. Resultados Principales

• Curvas de Luz: Las curvas de luz bolométricas simuladas son más tenues y decaen más rápido que las de AT2017gfo. Esto se debe a que solo se incluye el eyecta dinámico, omitiendo el eyecta secular (más masivo y lento). Se observa que una mayor masa de eyecta resulta en curvas de luz más brillantes y con un decaimiento más lento.
• Características Espectrales:
  • Se identifican con éxito las líneas de Sr II (triplete), que forman una característica tipo P-Cygni cerca del pico del SED.
  • Se detectan rasgos de La III y Gd III (mezclados) y de Ce III, los cuales coinciden con las observaciones de AT2017gfo.
  • Se observa que los iones de Cerio (Ce II y Ce III) juegan un papel crucial en la formación del continuo espectral.
  • En modelos específicos con baja fracción de lantánidos (como DD2 1.375–1.375), aparece una característica de Y II ($\lambda 7881$), posiblemente impulsada por fluorescencia.
• Desfase Temporal y de Velocidad: Las simulaciones coinciden con las observaciones de AT2017gfo en términos de forma espectral, pero ocurren en tiempos mucho más tempranos y presentan velocidades de línea (blueshift) generalmente más altas que las observadas.

5. Significancia y Conclusiones

El estudio demuestra que la señal espectral en la dirección polar es una predicción robusta de las fusiones de estrellas de neutrones, ya que la mayoría de los modelos muestran rasgos similares a pesar de las variaciones en la masa o la EoS.

La investigación subraya que, si bien el eyecta dinámico domina la evolución temprana, se requieren simulaciones futuras que incluyan el eyecta secular y efectos de no-equilibrio termodinámico (non-LTE) para resolver las discrepancias en la luminosidad y el exceso de absorción en longitudes de onda azules. Este trabajo sienta las bases para interpretar las futuras observaciones de kilonovas mediante la detección de estas firmas químicas específicas.