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

Résumé Technique : Transport des Particules Bêta et Thermalisation dans l'Éjecta de Kilonova avec une Microphysique Atomique Détaillée

Énoncé du Problème
L'émission des kilonovae (KN), alimentée par la désintégration radioactive de noyaux issus du processus $r$ dans l'éjecta de fusion d'étoiles à neutrons, est une sonde critique pour comprendre la nucléosynthèse, l'évolution des binaires et l'équation d'état de la matière dense. La modélisation précise des courbes de lumière des KN nécessite des calculs rigoureux du taux de chauffage fourni par la thermalisation des produits de désintégration ($\beta$-particules, $\alpha$-particules, $\gamma$-rayons et fragments de fission). Si le transport des $\gamma$-rayons est souvent traité par un transport de rayonnement complet, la thermalisation des particules chargées ($\beta$ et $\alpha$) a historiquement reposé sur des descriptions simplifiées. Ces simplifications introduisent des incertitudes significatives, particulièrement concernant l'efficacité du dépôt d'énergie, la distribution spatiale du chauffage et l'état d'ionisation de l'éjecta. De plus, le transport des particules chargées est sensible à la géométrie du champ magnétique, qui reste largement non contrainte. Cet article traite de l'absence de microphysique atomique détaillée et de la nécessité d'un cadre de transport relativiste pour mieux modéliser la thermalisation des $\beta$-particules dans les éjectas en expansion des kilonovae.

Méthodologie
Les auteurs emploient une approche multidimensionnelle combinant une physique atomique détaillée avec un code de transport Monte Carlo relativiste :

1. Données Microphysiques et Sections Efficaces :

   • Rendements Nucléaires et Spectres : Les rendements nucléaires pour les trajectoires du processus $r$ fortes, moyennes et faibles sont calculés à l'aide du réseau de réactions Prism. Les taux et spectres de désintégration $\beta$ sont issus de l'IAEA Livechart (ENSDF) et calculés via le modèle QRPA+HF, complétés par BetaShape pour les noyaux quasi-stables.
   • Interactions Électron-Atome : L'étude utilise la bibliothèque Evaluated Electron Data Library (EEDL) de la base de données EPICS2023 pour déterminer les sections efficaces pour l'ionisation par impact électronique (EII), l'excitation par impact électronique (EIE), le rayonnement de freinage (bremsstrahlung) et la diffusion de Mott.
   • Corrections d'Ionisation : Les sections efficaces des atomes neutres sont ajustées pour les cibles ionisées à l'aide de la mise à l'échelle empirique de Lotz et des énergies de liaison Hartree-Fock.
   • Électrons Secondaires : Les auteurs calculent la « cascade d'ionisation », en suivant la production d'électrons secondaires et tertiaires à partir des $\beta$-particules primaires pour déterminer le taux d'ionisation total.
   • Autres Interactions : La diffusion de Møller (électron-électron) est calculée de manière analytique. Les pertes par rayonnement synchrotron et Cherenkov sont évaluées et jugées négligeables par rapport aux pertes collisionnelles.

2. Modèle de Transport Global :

   • Géométrie de l'Éjecta : L'éjecta est modélisé comme un flux en expansion homologue à symétrie sphérique avec un profil de densité en loi de puissance ($\rho \propto r^{-\alpha}$).
   • Scénarios de Transport : Deux géométries limites de champ magnétique sont considérées :
     • Limite Piégée (Trapped Limit) : Correspond à des champs hautement emmêlés ou toroïdaux où les particules sont confinées à des coordonnées de vitesse spécifiques.
     • Limite de Champ Radial (Radial Field Limit) : Correspond à des champs radiaux cohérents où les particules circulent librement le long des lignes de champ. Les auteurs démontrent que pour les champs radiaux cohérents, les équations de transport 1D sont mathématiquement équivalentes au cas de libre parcours non magnétisé en raison de l'annulation de la force de miroir magnétique et des effets géométriques.
   • Simulation Monte Carlo : Un code Monte Carlo entièrement relativiste suit des paquets discrets de $\beta$-particules. Il prend en compte le refroidissement adiabatique dû à l'expansion, la perte d'énergie collisionnelle (pouvoir d'arrêt) et la diffusion (opacité de transport). Les particules sont terminées lors de la thermalisation ($E < 1$ eV), de l'échappement ($u > v_{ej}$), ou après avoir atteint un temps maximum (1 an).

3. Structure d'Ionisation :

   • Un modèle simplifié d'équation de taux est appliqué au Néodyme (Nd) pour illustrer l'impact du transport et de la cascade d'ionisation sur la température et l'état d'ionisation de l'éjecta, en négligeant la photoionisation et le couplage entre différents éléments.

Contributions Clés et Résultats

• Microphysique Atomique : Les auteurs fournissent des sections efficaces de transport et de pouvoir d'arrêt évaluées par espèce, ainsi que des sections efficaces EII, disponibles en ligne. Ils montrent qu'à des énergies de l'ordre du MeV, l'EII domine la perte d'énergie, tandis qu'à des énergies plus basses (phase de thermalisation), les interactions coulombiennes avec les électrons thermiques libres (diffusion de Møller) deviennent dominantes.
• Cascade d'Ionisation : L'étude quantifie la contribution des électrons secondaires au taux d'ionisation total. La cascade augmente significativement la section efficace d'ionisation, contribuant autant que les électrons primaires à 1 MeV.
• Efficacité de Thermalisation :
  • Effets Non-Locaux : Le dépôt d'énergie non local et l'échappement des particules réduisent l'efficacité de la thermalisation, particulièrement dans les parties les plus internes et les plus externes de l'éjecta. Cela rend l'efficacité plus uniforme à travers l'éjecta par rapport aux modèles de dépôt local.
  • Limites de Transport : Dans la limite du champ radial, l'échappement des particules est le principal moteur de l'inefficacité. Dans la limite du champ piégé, l'inefficacité est pilotée par les pertes par expansion adiabatique car les particules persistent plus longtemps dans l'éjecta.
  • Dépendance aux Paramètres : L'efficacité de la thermalisation est sensible à la masse de l'éjecta, à la vitesse, aux motifs d'abondance du processus $r$ et aux pentes du profil de densité. Elle est moins sensible à l'état d'ionisation, suggérant que le dépôt d'énergie peut être modélisé sans calculs d'ionisation auto-cohérents dans certains contextes.
• Prescription Analytique : Les auteurs dérivent une formule analytique pour l'efficacité de thermalisation spatialement dépendante, $f_{thm}(t) = \ln(1+\tau)/\tau$, paramétrée par un temps d'inefficacité ($t_{ineff}$) et une pente asymptotique ($n$). Ils fournissent des paramètres d'ajustement pour ces valeurs en fonction de la masse, de la vitesse et de la position de l'éjecta. Les résidus médians entre cette prescription et les résultats Monte Carlo sont de $\sim 0,04$ dex (radial) et $\sim 0,02$ dex (piégé).
• Structure d'Ionisation : Le calcul illustratif du Nd révèle une structure d'ionisation « inversée » dans la phase nébulaire (ionisation plus élevée dans les régions extérieures moins denses) en raison d'une recombinaison moins efficace. Les champs magnétiques radiaux conduisent à des états d'ionisation plus bas en raison de l'échappement des particules, tandis que l'inclusion de la cascade d'ionisation augmente l'état d'ionisation dans tout l'éjecta.

Signification et Revendications
L'article affirme que la prochaine génération de modèles de kilonova doit incorporer des données atomiques évaluées et un transport explicite des particules chargées pour réduire les incertitudes des prédictions des courbes de lumière. Les auteurs soutiennent que :

1. La Physique du Transport est Critique : Ignorer les effets de transport non locaux conduit à surestimer le chauffage dans les zones internes et externes et à obtenir des structures d'ionisation incorrectes.
2. Les Données Atomiques Comptent : L'utilisation de données atomiques évaluées (EEDL) et l'inclusion de la cascade d'ionisation sont nécessaires pour des taux d'ionisation précis, qui influencent à leur tour l'opacité et le refroidissement.
3. Utilité Analytique : Les prescriptions analytiques fournies offrent une approximation peu coûteuse en calcul de l'efficacité de thermalisation spatialement dépendante, adaptée aux modèles semi-analytiques et de transfert de rayonnement, bien qu'un transport explicite soit recommandé pour les simulations de haute fidélité.
4. Limites : Les auteurs notent avec modestie que leurs résultats sont basés sur une symétrie sphérique idéalisée et des géométries de champs magnétiques limites. Ils suggèrent que les morphologies de champs magnétiques réalistes, les éjectas non sphériques et les effets de plasma collectifs (par exemple, l'instabilité de flux) sont des domaines de recherche futurs. Ils soulignent également que leurs résultats se concentrent sur les $\beta$-particules et doivent être intégrés aux travaux futurs sur les $\alpha$-particules et les fragments de fission pour une image complète de la thermalisation.