Beta-Particle Transport and Thermalization in Kilonova Ejecta with Detailed Atomic Microphysics
Abstract
When two neutron stars collide, they eject material containing heavy nuclei formed by the rapid neutron capture process (r-process). As these nuclei decay, they power a bright optical/near-infrared transient known as a kilonova (KN). Modeling KN emission is a complex problem involving atomic opacities, radiation transport, and heating powered by the thermalization of radioactive decay products like γ-rays, α-particles, and β-particles. For heating by γ-rays, many KN modeling codes do full radiation transport calculations. However, heating by α- and β-particles relies on simplified descriptions of collisions and transport, and remains an important source of uncertainty in KN models. In this paper, we study the thermalization and transport of β-particles. To study thermalization, we use evaluated atomic physics data to estimate per-species contributions to energy deposition, scattering, and electron impact ionization, which we make available online. To include non-local effects, we develop a fully relativistic framework for charged particle transport in a spherically symmetric, homologously expanding ejecta, considering two limiting magnetic-field geometries. Non-local energy deposition and escape reduce thermalization efficiency, especially in the innermost and outermost ejecta, lowering the ejecta temperature and ionization state compared to local deposition models. Coulomb scattering partially offsets these effects by trapping particles at intermediate times ∼ 3−30 days. Ionization by secondary electrons significantly enhances the overall ionization rate. We provide analytic prescriptions for the spatially dependent thermalization efficiency for use in future light-curve calculations. Our results demonstrate that evaluated atomic data and charged-particle transport should be incorporated into the next generation of KN models.
Related
- Automated optimization of solar energetic particle acceleration and transport parameters against five-spacecraft observations of the 9 October 2021 event
- Automated optimization of solar energetic particle acceleration and transport parameters against five-spacecraft observations of the 9 October 2021 event
- Multispacecraft Constrained Modeling of Solar Energetic Particle Transport across Different Heliocentric Radii and Heliolongitudes
- STORM: RDMA-based Monte Carlo Transport Scheme for Distributed-Memory Particle Simulations
- A comprehensive simulation framework for multi-modal kilonova observations from all-sky surveys
BibTeX
@misc{andalman2026betaparticle,
title = {{Beta-Particle Transport and Thermalization in Kilonova Ejecta with Detailed Atomic Microphysics}},
author = {Zachary L. Andalman and Christopher L. Fryer and Christopher J. Fontes and Matthew R. Mumpower and Ryan T. Wollaeger},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2607.01228},
eprint = {2607.01228},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2607.01228},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Beta-Particle Transport and Thermalization in Kilonova Ejecta with Detailed Atomic Microphysics
AU - Zachary L. Andalman
AU - Christopher L. Fryer
AU - Christopher J. Fontes
AU - Matthew R. Mumpower
AU - Ryan T. Wollaeger
PY - 2026
DA - 2026/10/06
PB - arXiv
DO - 10.48550/arxiv.2607.01228
UR - https://arxiv.org/abs/2607.01228
AB - When two neutron stars collide, they eject material containing heavy nuclei formed by the rapid neutron capture process (r-process). As these nuclei decay, they power a bright optical/near-infrared transient known as a kilonova (KN). Modeling KN emission is a complex problem involving atomic opacities, radiation transport, and heating powered by the thermalization of radioactive decay products like γ-rays, α-particles, and β-particles. For heating by γ-rays, many KN modeling codes do full radiation transport calculations. However, heating by α- and β-particles relies on simplified descriptions of collisions and transport, and remains an important source of uncertainty in KN models. In this paper, we study the thermalization and transport of β-particles. To study thermalization, we use evaluated atomic physics data to estimate per-species contributions to energy deposition, scattering, and electron impact ionization, which we make available online. To include non-local effects, we develop a fully relativistic framework for charged particle transport in a spherically symmetric, homologously expanding ejecta, considering two limiting magnetic-field geometries. Non-local energy deposition and escape reduce thermalization efficiency, especially in the innermost and outermost ejecta, lowering the ejecta temperature and ionization state compared to local deposition models. Coulomb scattering partially offsets these effects by trapping particles at intermediate times ∼ 3−30 days. Ionization by secondary electrons significantly enhances the overall ionization rate. We provide analytic prescriptions for the spatially dependent thermalization efficiency for use in future light-curve calculations. Our results demonstrate that evaluated atomic data and charged-particle transport should be incorporated into the next generation of KN models.
ER -