The effect of inverse Compton losses on particle acceleration in three-dimensional relativistic reconnection

Résumé technique : L'effet des pertes par Compton inverse sur l'accélération des particules dans la reconnexion relativiste tridimensionnelle

Énoncé du problème
La reconnexion magnétique relativiste est un mécanisme primaire de dissipation de l'énergie magnétique et d'accélération des particules dans les environnements astrophysiques tels que les magnétosphères d'objets compacts et les jets relativistes. Des simulations récentes de type « particle-in-cell » (PIC) fondées sur les premiers principes ont établi un modèle d'accélération des particules dans la reconnexion 3D (Zhang et al. 2023), distinguant une « phase libre » où les particules accélèrent activement dans la région amont, et une « phase piégée » où les particules sont confinées dans des cordes de flux en aval sans énergisation significative. Ce modèle prédit un spectre en loi de puissance universel ($dN/d\gamma \propto \gamma^{-1}$) pour les particules libres et un spectre plus raide ($dN/d\gamma \propto \gamma^{-2}$) pour les particules piégées en l'absence de refroidissement radiatif.

Cependant, le comportement de ce mécanisme d'accélération en présence de refroidissement radiatif — spécifiquement les pertes par Compton inverse (IC) — reste sous-exploré. Il n'est pas clair si les pertes par refroidissement modifient le taux de reconnexion, la structure de la couche ou la physique fondamentale des phases libre et piégée, particulièrement dans le régime de « faible refroidissement » où le facteur de Lorentz limité par la réaction de rayonnement ($\gamma_{\text{rad}}$) excède la magnétisation ($\sigma$).

Méthodologie
Les auteurs ont réalisé une série de simulations PIC 3D en utilisant le code TRISTAN-MP pour étudier la reconnexion relativiste avec des pertes par refroidissement IC. Les simulations ont utilisé une configuration de feuillet de courant de Harris avec un champ guide ($B_g = 0,1 B_0$) et des conditions aux limites périodiques le long de la direction du courant.

Les paramètres clés de la simulation comprenaient :

• Composition du plasma : Paires électron-positron, avec un refroidissement radiatif appliqué uniquement aux électrons. Les positrons n'ont pas été refroidis afin de simuler le comportement des ions lourds dans les plasmas électron-ion et pour assurer une altération minimale de l'équilibre des forces de la couche, permettant une comparaison directe avec le modèle non refroidi de Zhang et al. (2023).
• Magnétisation : Une magnétisation fixe de $\sigma = 10$.
• Régime de refroidissement : L'étude s'est concentrée sur le régime de faible refroidissement ($\gamma_{\text{rad}} > \sigma$). Trois valeurs spécifiques pour le facteur de Lorentz limité par la réaction de rayonnement ont été testées : $\gamma_{\text{rad}} = 30, 43$ et $60$. Celles-ci correspondent à des ruptures de refroidissement ($\gamma_{\text{cool}}$) d'environ $3,5, 7,2$ et $14,2$, respectivement.
• Analyse : Les auteurs ont suivi de grands échantillons de particules (0,5 million d'électrons par exécution) pour distinguer les populations « libres » (amont) et « piégées » (aval) en utilisant un critère de facteur de mélange. Ils ont analysé les taux d'accélération, les temps d'échappement et les spectres d'énergie, en les comparant aux prédictions analytiques dérivées du modèle de Zhang et al. (2023) étendu pour inclure le refroidissement.

Contributions clés et résultats

1. Robustesse de la dynamique de reconnexion : Les simulations démontrent que les pertes par refroidissement électronique n'altèrent pas de manière significative le taux de reconnexion ou la structure macroscopique de la couche de courant. Le taux de reconnexion en régime permanent se stabilise à $\eta_{\text{rec}} \approx 0,05$, ce qui est cohérent avec les cas non refroidis. Bien qu'un refroidissement plus intense conduise à des couches légèrement plus fines et à des plasmoïdes plus petits, la physique fondamentale de la couche reste inchangée à travers les intensités de refroidissement testées.

2. Validation de la physique de la phase libre : L'étude confirme que la « phase libre » d'accélération est insensible aux pertes par refroidissement. Les électrons libres accélèrent systématiquement au taux optimal ($\dot{\gamma}_{\text{acc}} \sim \eta_{\text{rec}}\omega_c$) et présentent un spectre en loi de puissance de $dN_{\text{free}}/d\gamma \propto \gamma^{-1}$, identique au cas non refroidi. La coupure de ce spectre est déterminée par $\gamma_{\text{rad}}$ plutôt que par la taille du système.

3. Modification du spectre de la phase piégée : Une découverte critique est la modification du spectre des particules piégées en présence de refroidissement.

   • Dans le cas non refroidi, le spectre des particules piégées suit une loi de type $dN/d\gamma \propto \gamma^{-2}$.
   • Dans le régime de faible refroidissement, pour les énergies supérieures à la rupture de refroidissement ($\gamma_{\text{cool}} < \gamma < \gamma_{\text{rad}}$), le spectre des électrons piégés s'accentue vers $dN/d\gamma \propto \gamma^{-3}$.
   • Cet accentuation confirme l'hypothèse fondamentale du modèle de Zhang et al. (2023) : aucune énergisation significative ne se produit pendant la phase piégée. L'accentuation provient du fait que les particules piégées se refroidissent sur une échelle de temps comparable ou plus courte que le temps d'advection, réduisant ainsi leur population aux énergies plus élevées par rapport au cas non refroidi.

4. Cinématique des particules : Les simulations vérifient que les particules libres se déplacent presque le long du champ électrique de reconnexion ($\beta_z \approx -1$), facilitant une accélération rapide, tandis que les particules piégées présentent une distribution symétrique de $\beta_z$ autour de zéro, cohérente avec l'absence d'une énergisation systématique.

Signification et implications astrophysiques
L'article valide le modèle analytique proposé par Zhang et al. (2023) pour la reconnexion relativiste 3D en présence d'un faible refroidissement radiatif. La principale implication est que les modèles d'émission radiative pour les sources alimentées par la reconnexion (tels que les sursauts gamma et les jets de blazars) doivent employer une structure à deux zones :

• Zone 1 (Libre) : Particules accélérant rapidement avec un spectre $dN/d\gamma \propto \gamma^{-1}$.
• Zone 2 (Piégée) : Particules subissant un refroidissement passif avec un spectre qui s'accentue de $\gamma^{-2}$ à $\gamma^{-3}$ au-dessus de la rupture de refroidissement.

Les auteurs notent que bien que la majorité des particules émettrices soient piégées, les particules libres deviennent plus abondantes près du facteur de Lorentz maximal ($\gamma_{\text{rad}}$). Comme les particules libres sont préférentiellement dirigées le long du champ électrique de reconnexion, leur émission pourrait produire des signatures observationnelles uniques, telles qu'une variabilité temporelle rapide.

Limites et mises en garde
Les auteurs précisent explicitement que leurs résultats s'appliquent au régime de faible refroidissement où $\gamma_{\text{rad}} > \sigma$. Le modèle n'est pas applicable aux scénarios de refroidissement fort où $\gamma_{\text{rad}} < \sigma$. De plus, la plage dynamique entre $\gamma_{\text{cool}}$ et $\gamma_{\text{rad}}$ dans les simulations est plus petite que dans les systèmes astrophysiques réels. Des travaux futurs sont suggérés pour étendre ces conclusions aux plasmas électron-proton avec des rapports de masse réalistes et à la reconnexion avec des champs guides significatifs.