Non-linear Evolution of Dark Plasma Subhalos

Résumé Technique : Évolution Non Linéaire des Sous-halos de Plasma Sombre

Énoncé du Problème
Cet article étudie l'évolution dynamique des sous-halos de matière noire (DM) en orbite dans le potentiel de la Voie Lactée (MW) sous l'hypothion que la DM constitue un « plasma sombre ». Dans les modèles où la DM interagit avec elle-même via des forces à longue portée médiées par un photon caché (un secteur $U(1)$ sombre), la dynamique ressemble à celle des plasmas du Modèle Standard. Plus précisément, lorsqu'un sous-halo de DM se déplace à travers le plasma sombre ambiant de la MW, la vitesse de dérive relative entre les deux populations peut excéder leurs dispersions de vitesse internes. Cette configuration est susceptible de générer des instabilités de plasma (par exemple, les instabilités de deux flux et de faisceau-plasma), qui produisent des champs électrostatiques et électromagnétiques sombres. Ces champs diffusent les particules, modifient la distribution de phase et poussent le système vers un équilibre thermique. Le problème central abordé est l'évolution non linéaire de ces instabilités et leur impact sur la perte de masse des sous-halos, en comparant spécifiquement le chauffage induit par le plasma au décapage tidale gravitationnel standard.

Méthodologie
Les auteurs emploient une approche hybride semi-analytique et numérique pour modéliser l'évolution des sous-halos :

1. Cadre Semi-Analytique :

   • Conditions Initiales : Les sous-halos sont modélisés comme des plasmas sphériques à neutralité nette, présentant des profils de densité NFW et des distributions de vitesse de Maxwell-Boltzmann. Deux masses de référence sont considérées : $10^7 M_\odot$ et $10^9 M_\odot$.
   • Dynamique Orbitale : Les orbites des sous-halos sont intégrées à l'aide du package galpy dans un potentiel de MW (MWPotential2014). Les orbites sont initialisées au rayon viriel avec des eccentricités variables ($\epsilon = 0,1$ à $0,9$).
   • Couplage de la Perte de Masse : Le profil de densité du sous-halo est discrétisé en couches radiales concentriques. L'évolution de chaque couche est suivie sur des pas de temps ($\delta t = 7$ Myr) jusqu'à $1,2 \times t_{\text{peri}}$ (temps du premier passage au périastre).
   • Mécanismes : La perte de masse est calculée via deux canaux distincts :
     • Décapage Tidal (Tidal Stripping) : Modélisé à l'aide d'un paramètre d'efficacité standard ($\alpha = 0,55$) appliqué aux couches situées à l'extérieur du rayon tidal instantané.
     • Chauffage par Plasma : La fraction de particules liées ($f_{\text{bound}}$) dans chaque couche est déterminée par interpolation des résultats d'une suite de simulations de type Particle-in-Cell (PIC). Les particules chauffées au-dessus de la vitesse d'échappement locale ($v_{\text{esc}} \approx 4\sigma_{\text{SH}}$) sont éliminées, et la masse restante liée se revirialise.

2. Simulations Particle-in-Cell (PIC) :

   • Code : L'étude utilise le code TRISTAN-MP, adapté aux plasmas sombres via l'invariance de mise à l'échelle des coordonnées des équations de Vlasov-Poisson.
   • Configuration : Des simulations électrostatiques unidimensionnelles (1D) sont effectuées avec des conditions aux limites périodiques. La configuration modélise des faisceaux en contre-flux représentant le sous-halo et le fond de la MW.
   • Paramètres : Les simulations font varier le rapport de densité ($\rho_{\text{SH}}/\rho_{\text{MW}} \in [10^{-3}, 10^2]$) et la vitesse de streaming relative ($v_{\text{SH}}/\sigma_{\text{MW}}$). Le rapport charge/masse de référence est $q_\chi/m_\chi = 10^{-14} e/m_p$.
   • Dimensionnalité : Bien que les principaux résultats soient en 1D, les auteurs réalisent une simulation 2D-2V (deux dimensions spatiales, deux dimensions de vitesse) limitée pour vérifier que l'efficacité du chauffage augmente dans des dimensions supérieures, suggérant que les résultats 1D sont conservateurs.

Contributions Clés

• Première Étude Non Linéaire : Ce travail présente la première étude de l'évolution non linéaire des instabilités électrostatiques dans les sous-halos de DM en utilisant des simulations cinétiques PIC complètes.
• Quantification du Chauffage par Plasma : Les auteurs quantifient l'efficacité de la perte de masse induite par le plasma, démontrant qu'elle peut excéder significativement le décapage tidal, particulièrement pour les sous-halos de faible masse sur des orbites excentriques.
• Modèle Semi-Analytique : Un nouveau cadre est développé, couplant la croissance locale de l'instabilité du plasma (dérivée des données PIC) avec l'évolution orbitale globale pour prédire l'évolution de la masse et du profil de densité du sous-halo.

Résultats

• Efficacité de la Perte de Masse : Pour les sous-halos sur des orbites hautement excentriques ($\epsilon = 0,9$), le chauffage par plasma est le mécanisme de perte de masse dominant. Au moment du premier périastre ($t = 1,2 t_{\text{peri}}$) :
  • Un sous-halo de $10^7 M_\odot$ perd environ $97 \%$ de sa masse initiale.
  • Un sous-halo de $10^9 M_\odot$ perd environ $84 \%$ de sa masse initiale.
  • En revanche, le décapage tidal seul ne retire que $\sim 30 \%$ de la masse pour les deux cas.
• Dépendance Radiale (« Outside-In ») : La perte de masse progresse de l'extérieur vers l'intérieur. Les instabilités croissent plus efficacement là où le contraste de densité entre le sous-halo et la MW est modeste ($\rho_{\text{SH}}/\rho_{\text{MW}} \sim 1$) et la vitesse relative est élevée. Par conséquent, les couches externes sont chauffées et décapées avant le cœur interne. Pour les orbites très excentriques, le cœur est finalement affecté près du périastre ; pour les orbites circulaires, la perte de masse est confinée aux régions externes.
• Implications sur la Fonction de Masse : Le décapage préférentiel des sous-halos de faible masse (en raison de leurs dispersions de vitesse internes et vitesses d'échappement plus faibles) suggère une suppression de l'extrémité de faible masse de la Fonction de Masse des Sous-halos (SHMF). Les auteurs estiment que pour les paramètres de référence, la SHMF pourrait présenter une pente plus faible ou une coupure nette en dessous de $M_{\text{SH}} \sim 10^7 M_\odot$.
• Excentricité Orbitale : L'effet est très sensible à l'excentricité. Pour les orbites circulaires ($\epsilon = 0,1$), le chauffage par plasma ne retire que $\sim 50 \%$ ($10^7 M_\odot$) et $\sim 30 \%$ ($10^9 M_\odot$) de la masse, ce qui est comparable ou légèrement supérieur aux effets tidaux.

Signification et Revendications
L'article affirme que la physique du plasma sombre offre un mécanisme distinct pour l'évolution des sous-halos capable de laisser des signatures observables dans la population de satellites de la MW. Plus précisément, les auteurs soutiennent que le chauffage et le décapage induits par le plasma peuvent :

1. Remodeler la SHMF : En détruisant préférentiellement les sous-halos de faible masse et « froids », les modèles de plasma sombre peuvent résoudre ou altérer le problème des « satellites manquants » ou le problème du « too big to fail », menant potentiellement à une pente supprimée de la faible masse dans la population de sous-halos observée.
2. Fournir des Contraintes Observationnelles : La suppression des sous-halos de faible masse et la modification résultante de la SHMF offrent une nouvelle voie pour contraindre les paramètres du secteur sombre (spécifiquement le rapport charge/masse $q_\chi/m_\chi$) en utilisant les données de relevés tels que ELVES, SPARC, SAGA et Merian.
3. Souligner la Dynamique Non Linéaire : L'étude souligne que les effets de plasma collectifs, pilotés par les instabilités de flux, peuvent opérer sur des échelles de temps comparables aux échelles de temps dynamiques galactiques, conduisant à une évaporation rapide que les simulations CDM sans collisions classiques ne capturent pas.

Les auteurs maintiennent un ton modeste concernant la portée de leurs résultats, notant que leurs conclusions reposent sur des simulations 1D (bien que des tests 2D suggèrent que la 1D est conservatrice) et qu'une prédiction complète au niveau de la population nécessite des travaux futurs intégrant la gravité et les effets de plasma de manière auto-cohérente dans des dimensions supérieures.