Beta-Particle Transport and Thermalization in Kilonova Ejecta with Detailed Atomic Microphysics

Resumen Técnico: Transporte de Partículas Beta y Termalización en Eyecta de Kilonova con Microfísica Atómica Detallada

Planteamiento del Problema
La emisión de las kilonovas (KN), impulsada por el decaimiento radiactivo de núcleos del proceso-$r$ en la eyecta de la fusión de estrellas de neutrones, es una sonda crítica para comprender la nucleosíntesis, la evolución binaria y la ecuación de estado de la materia densa. El modelado preciso de las curvas de luz de las KN requiere cálculos precisos de la tasa de calentamiento proporcionada por la termalización de los productos de decaimiento (partículas $\beta$, partículas $\alpha$, rayos $\gamma$ y fragmentos de fisión). Mientras que el transporte de rayos $\gamma$ suele tratarse con un transporte de radiación completo, la termalización de las partículas cargadas ($\beta$ y $\alpha$) ha dependido históricamente de descripciones simplificadas. Estas simplificaciones introducen incertidumbres significativas, particularmente respecto a la eficiencia de la deposición de energía, la distribución espacial del calentamiento y el estado de ionización resultante de la eyecta. Además, el transporte de partículas cargadas es sensible a la geometría del campo magnético, la cual permanece mayormente sin restringir. Este artículo aborda la falta de microfísica atómica detallada y la necesidad de un marco de transporte relativista para modelar mejor la termalización de partículas $\beta$ en la eyecta en expansión de las kilonovas.

Metodología
Los autores emplean un enfoque multifacético que combina física atómica detallada con un código de transporte de Monte Carlo relativista:

1. Datos Microfísicos y Secciones Eficaces:

   • Rendimientos Nucleares y Espectros: Los rendimientos nucleares para trayectorias del proceso-$r$ fuertes, medias y débiles se calculan utilizando la red de reacciones Prism. Las tasas y espectros de decaimiento $\beta$ se obtienen de la base de datos IAEA Livechart (ENSuff) y se computan mediante el modelo QRPA+HF, complementado con BetaShape para núcleos casi estables.
   • Interacciones Electrón-Átomo: El estudio utiliza la Biblioteca de Datos Electrónicos Evaluados (EEDL) de la base de datos EPICS2023 para determinar las secciones eficaces de ionización por impacto electrónico (EII), excitación por impacto electrónico (EIE), bremsstrahlung y dispersión de Mott.
   • Correcciones de Ionización: Las secciones eficaces de átomos neutros se ajustan para objetivos ionizados utilizando el escalamiento empírico de Lotz y energías de enlace de Hartree-Fock.
   • Electrones Secundarios: Los autores calculan la "cascada de ionización", rastreando la producción de electrones secundarios y terciarios a partir de las partículas $\beta$ primarias para determinar la tasa de ionización total.
   • Otras Interacciones: La dispersión de Møller (electrón-electrón) se calcula analíticamente. Las pérdidas por sincrotrón y Cherenkov se evalúan y se encuentran insignificantes en comparación con las pérdidas colisionales.

2. Modelo de Transporte Global:

   • Geometría de la Eyecta: La eyecta se modela como un flujo de expansión homóloga con simetría esférica y un perfil de densidad de ley de potencia ($\rho \propto r^{-\alpha}$).
   • Escenarios de Transporte: Se consideran dos geometrías límite de campo magnético:
     • Límite de Atrapamiento (Trapped): Corresponde a campos altamente enredados o toroidales donde las partículas están confinadas a coordenadas de velocidad específicas.
     • Límite de Campo Radial: Corresponde a campos radiales coherentes donde las partículas fluyen libremente a lo largo de las líneas de campo. Los autores demuestran que, para campos radiales coherentes, las ecuaciones de transporte 1D son matemáticamente equivalentes al caso no magnetizado de flujo libre debido a la cancelación de la fuerza de espejo magnético y los efectos geométricos.
   • Simulación de Monte Carlo: Un código de Monte Carlo totalmente relativista rastrea paquetes discretos de partículas $\beta$. Este tiene en cuenta el enfriamiento adiabático debido a la expansión, la pérdida de energía por colisión (poder de frenado) y la dispersión (opacidad de transporte). Las partículas terminan su trayectoria tras la termalización ($E < 1$ eV), la escape ($u > v_{ej}$) o al alcanzar un tiempo máximo (1 año).

3. Estructura de Ionización:

   • Se aplica un modelo simplificado de ecuaciones de tasa al Neodimio (Nd) para ilustrar el impacto del transporte y la cascada de ionización en la temperatura y el estado de ionización de la eyecta, omitiendo la fotoionización y el acoplamiento entre diferentes elementos.

Contribuciones Clave y Resultados

• Microfísica Atómica: Los autores proporcionan secciones eficaces de transporte y de frenado evaluadas por especie, incluyendo las secciones eficaces de EII, disponibles en línea. Demuestran que a energías de MeV, la EII domina la pérdida de energía, mientras que a energías más bajas (fase de termalización), las interacciones de Coulomb con electrones térmicos libres (dispersión de Møller) se vuelven dominantes.
• Cascada de Ionización: El estudio cuantifica la contribución de los electrones secundarios a la tasa de ionización total. La cascada aumenta significativamente la sección eficaz de ionización, contribuyendo tanto como los electrones primarios a 1 MeV.
• Eficiencia de Termalización:
  • Efectos No Locales: La deposición de energía no local y el escape de partículas reducen la eficiencia de la termalización, particularmente en las regiones más internas y externas de la eyecta. Esto hace que la eficiencia sea más uniforme en toda la eyecta en comparación con los modelos de deposición local.
  • Límites de Transporte: En el límite de campo radial, el escape de partículas es el principal motor de la ineficiencia. En el límite de atrapamiento, la ineficiencia es impulsada por las pérdidas de expansión adiabática a medida que las partículas persisten más tiempo en la eyecta.
  • Dependencia de Parámetros: La eficiencia de la termalización es sensible a la masa de la eyecta, la velocidad, los patrones de abundancia del proceso-$r$ y las pendientes del perfil de densidad. Es menos sensible al estado de ionización, lo que sugiere que la deposición de energía puede modelarse sin cálculos de ionización auto-consistentes en ciertos contextos.
• Prescripción Analítica: Los autores derivan una fórmula analítica para la eficiencia de termalización espacialmente dependiente, $f_{thm}(t) = \ln(1+\tau)/\tau$, parametrizada por un tiempo de ineficiencia ($t_{ineff}$) y una pendiente asintótica ($n$). Proporcionan parámetros de ajuste para estos valores en función de la masa de la eyecta, la velocidad y la posición. Los residuos medianos entre esta prescripción y los resultados de Monte Carlo son de $\sim0.04$ dex (radial) y $\sim0.02$ dex (atrapado).
• Estructura de Ionización: El cálculo ilustrativo de Nd revela una estructura de ionización "invertida" en la fase nebular (mayor ionización en las regiones exteriores menos densas) debido a una recombinación menos eficiente. Los campos magnéticos radiales conducen a estados de ionización más bajos debido al escape de partículas, mientras que la inclusión de la cascada de ionización eleva el estado de ionización en toda la eyecta.

Significancia y Reivindicaciones
El artículo afirma que la próxima generación de modelos de kilonovas debe incorporar datos atómicos evaluados y un transporte explícito de partículas cargadas para reducir las incertididades en las predicciones de las curvas de luz. Los autores sostienen que:

1. La Física de Transporte es Crítica: Ignorar los efectos de transporte no locales conduce a sobreestimar el calentamiento en las regiones interna y externa de la eyecta e incorrectas estructuras de ionización.
2. Los Datos Atómicos Importan: El uso de datos atómicos evaluados (EEDL) y la inclusión de la cascada de ionización son necesarios para tasas de ionización precisas, que a su vez afectan la opacidad y el enfriamiento.
3. Utilidad Analítica: Las prescripciones analíticas proporcionadas ofrecen una aproximación computacionalmente económica de la eficiencia de termalización espacialmente dependiente, apta para modelos semi-analíticos y de transferencia de radiación, aunque se recomienda el transporte explícito para simulaciones de alta fidelidad.
4. Limitaciones: Los autores señalan modestamente que sus resultados se basan en simetría esférica idealizada y geometrías de campo magnético límite. Sugieren que las morfologías de campo magnético realistas, la eyecta no esférica y los efectos de plasma colectivos (por ejemplo, la inestabilidad de streaming) son áreas para investigación futura. También enfatizan que sus resultados se centran en partículas $\beta$ y deben integrarse con trabajos futuros sobre partículas $\alpha$ y fragmentos de fisión para obtener una imagen completa de la termalización.