Spectra and Ionization Efficiencies of Charged Decay Particles in Kilonova Ejecta

Resumen Técnico: Espectros y Eficiencias de Ionización de Partículas de Decaimiento Cargadas en Ejecta de Kilonova

Enunciado del Problema
El estado de ionización de las ejectas de kilonova durante la fase nebulosa es un determinante crítico de la evolución del color observado y de la intensidad de las líneas de absorción y emisión. Si bien estudios anteriores han modelado extensamente los presupuestos de calentamiento de las kilonovas impulsados por decaimiento radiactivo (específicamente decaimiento $\beta$, decaimiento $\alpha$ y fisión espontánea), los tratamientos de ionización han dependido en gran medida de tasas calibradas para electrones de decaimiento $\beta$. Sin embargo, las partículas $\alpha$ y los fragmentos de fisión poseen propiedades físicas distintas: son significativamente más masivos, portan energías cinéticas iniciales más altas (de MeV a cientos de MeV) y exhiben comportamientos de frenado diferentes en comparación con los electrones $\beta$. La falta de un marco autoconsistente que vincule la termalización de estas partículas cargadas pesadas con sus secciones eficaces de ionización específicas genera incertidumbre respecto al estado de ionización de las ejectas donde el decaimiento $\alpha$ y la fisión dominan el presupuesto de calentamiento (por ejemplo, en escalas de tiempo de $\sim$10–100 días).

Metodología
Los autores desarrollan un marco teórico autoconsistente para calcular las historias de termalización, los espectros de degradación de partículas y las eficiencias de ionización para partículas $\alpha$ y fragmentos de fisión que se propagan a través de ejectas de kilonova en expansión. La metodología integra varios componentes físicos clave:

1. Cálculos del Poder de Frenado: La pérdida de energía por unidad de longitud de trayectoria se calcula sumando las contribuciones de las colisiones con electrones ligados y electrones libres del plasma.

   • Electrones Ligados: La ionización se trata como el canal dominante de pérdida de energía utilizando el formalismo clásico de Gryziński para colisiones binarias. Las secciones eficaces diferenciales se derivan asumiendo que la masa del objetivo es mucho menor que la masa del proyectil.
   • Electrones Libres: El frenado debido a electrones térmicos se calcula utilizando la fórmula de Bohr, aplicable dada la alta velocidad de los fragmentos de fisión y $\alpha$ en relación con los electrones térmicos en la fase nebulosa.
   • Carga Efectiva: Para los fragmentos de fisión, que retienen electrones ligados, la carga efectiva ($Z_{\text{eff}}$) se modela como una evolución de equilibrio basada en tasas de intercambio de carga (Peter & Meyer-ter Vehn), mientras que las partículas $\alpha$ se tratan como completamente ionizadas ($Z=2$).
   • Estructura Atómica: Las energías de enlace y las funciones de onda orbitales se derivan utilizando cálculos del Hamiltoniano relativista X2C (PySCF) para garantizar la consistencia entre los modelos de poder de frenado e ionización.

2. Termalización y Espectros:

   • La evolución de la energía de las partículas tiene en cuenta tanto las pérdidas colisionales como las pérdidas de energía adiabáticas debidas a la expansión de las ejectas. Esto conduce a una solución analítica para la energía $E(t)$ en función del tiempo.
   • Espectros de Partículas: Los autores derivan los espectros de energía instantáneos de las partículas utilizando dos enfoques:
     • El formalismo Spencer–Fano (ecuación integral) para determinar el espectro de degradación en estado estacionario para partículas pesadas donde las contribuciones de electrones secundarios son despreciables.
     • Un enfoque de ecuación de continuidad para modelar el espectro físico dependiente del tiempo, incorporando el perfil de inyección de núcleos en descomposición y el enfriamiento adiabático.

3. Eficiencia de Ionización:

   • Se calcula el "trabajo por par iónico" ($w$), definido como la energía cinética inicial dividida por el número total de ionizaciones producidas.
   • El estudio compara estos resultados con los electrones $\beta$, para los cuales los espectros de degradación incluyen contribuciones significativas de electrones secundarios y procesos Auger, modelados utilizando secciones eficaces de HULLAC y distribuciones de energía secundarias específicas.

Contribuciones y Resultados Clave

• Marco Unificado de Termalización e Ionización: El artículo proporciona el primer cálculo autoconsistente de poderes de frenado y secciones eficaces de ionización específicamente adaptados a partículas $\alpha$ y fragmentos de fisión en condiciones de kilonova, superando la suposición de que las calibraciones de decaimiento $\beta$ son aplicables universalmente.
• Robustez del Trabajo por Par Iónico: Un hallazgo central es que el trabajo por par iónico ($w$) es notablemente similar en todos los canales de descomposición ($\beta$, $\alpha$ y fisión) y especies objetivo.
  • Para iones de elementos pesados (por ejemplo, Xe, Nd, Sr), la eficiencia de ionización normalizada es casi idéntica independientemente de si el proyectil es un electrón $\beta$, una partícula $\alpha$ o un fragmento de fisión.
  • Los valores de $w$ generalmente oscilan entre $2.5I$ y $18.1I$ (donde $I$ es el potencial de ionización), dependiendo del estado de carga del objetivo y de la fracción de abundancia de electrones ($\chi_e$).
• Relaciones de Escala: El estudio identifica que, para elementos pesados, el trabajo por par iónico escala principalmente con el número efectivo de electrones en el ion objetivo y la fracción de electrones ambiental, en lugar del tipo específico de partícula de descomposición.
• Evolución Espectral: El artículo detalla cómo evolucionan los espectros de partículas con el tiempo. En tiempos tempranos (alta densidad), las pérdidas colisionales dominan y los espectros siguen la forma instantánea de Spencer–Fano. En tiempos tardíos (baja densidad), el enfriamiento adiabático se convierte en el mecanismo dominante de pérdida de energía, aplanando la porción de alta energía del espectro.
• Comparación con Electrones $\beta$: Si bien los electrones $\beta$ exhiben formas espectrales diferentes debido a las cascadas de electrones secundarios, sus valores resultantes de trabajo por par iónico son cuantitativamente similares a los de las partículas pesadas para objetivos de elementos pesados. Solo se observan discrepancias para elementos ligeros (He I–II), donde las mayores energías de enlace y las capas electrónicas limitadas de los objetivos afectan significativamente más a las partículas pesadas que a los electrones $\beta$.

Significado
Los autores concluyen que el estado de ionización de las ejectas de kilonova en tiempos tardíos depende débilmente del canal de decaimiento radiactivo dominante. Incluso en ejectas enriquecidas con actínidos donde el decaimiento $\alpha$ y la fisión dominan el presupuesto de calentamiento, la eficiencia de ionización resultante sigue siendo comparable a la del decaimiento $\beta$. Esta robustez surge porque la física colisional subyacente que gobierna tanto el poder de frenado como las secciones eficaces de ionización hace que la relación entre ionización y pérdida de energía permanezca aproximadamente constante en los rangos de energía relevantes. En consecuencia, los modelos de espectros de kilonova y evolución de color pueden basarse en un tratamiento unificado de la eficiencia de ionización sin necesidad de ajustar drásticamente los parámetros según la mezcla específica de isótopos radiactivos, siempre que la composición de las ejectas esté dominada por elementos pesados del proceso-r.