Cooling, conduction, compact objects: Gravothermal evolution of dissipative self-interacting dark matter halos

Résumé Technique : Évolution Gravothermique des Halos de Matière Noire Auto-Interagissante Dissipative

Énoncé du Problème
Les modèles de matière noire auto-interagissante (SIDM) prédisent souvent des processus radiatifs qui dissipent l'énergie, tels que le rayonnement de freinage (bremsstrahlung) ou la désintégration d'états liés collisionnellement excités. Si la SIDM élastique est connue pour induire une évolution gravothermique (formation de cœur suivie d'un effondrement), l'impact spécifique de la dissipation d'énergie sur ce processus reste insuffisamment exploré dans les simulations numériques. Les méthodes N-corps existantes pour les diffusions fréquentes de petits angles (fSIDM) supposent généralement des interactions élastiques, tandis que les modèles dissipatifs reposent souvent sur des approximations de fluide ou des hypothèses de sections efficaces isotropes qui peuvent ne pas capturer la physique des médiateurs légers (qui sont souvent orientés vers l'avant). De plus, des observations récentes d'un objet compact dans le système de lentille forte JVAS B1938+666 suggèrent un halo de SIDM en cours d'effondrement, mais les modèles élastiques nécessitent des sections efficaces irréalistement grandes pour expliquer l'effondrement à la redshift observée ($z=0,881$). Ce travail examine si les auto-interactions dissipatives peuvent modifier l'évolution gravothermique des halos isolés pour expliquer de telles observations avec des sections efficaces plus modérées.

Méthodologie
Les auteurs présentent la première extension du formalisme N-corps de la fSIDM (Fischer et al. 2021) pour inclure une dissipation effective.

• Implémentation Numérique : La méthode traite les diffusions fréquentes de petits angles comme une force de traînée et une diffusion de la quantité de mouvement transverse. Les auteurs dérivent une extension dissipative où la force de traînée est augmentée par un paramètre de dissipation $r_{\text{diss}}$, et où le coup de vitesse stochastique est modifié pour assurer que le système perde de l'énergie conformément à l'énergie rayonnée $k_0$ dans le processus de physique des particules sous-jacent. Le pas de temps est réduit par un facteur $r_{\text{diss}}$ pour maintenir la stabilité.
• Configuration de la Simulation : Ils simulent des halos NFW (Navarro-Frenk-White) isolés ($M \sim 10^{10} M_\odot$) en utilisant le code OpenGadget3. Deux séries de simulations sont effectuées :
  1. En faisant varier le taux de dissipation tout en maintenant constant le taux de conduction thermique.
  2. En faisant varier le taux de conduction thermique (via la section efficace de transfert $\sigma_T/m_\chi$) tout en maintenant constant le taux de refroidissement.
• Comparaisons : Les résultats sont comparés à :
  • Un modèle de fluide gravothermal (augmenté d'un terme de refroidissement) pour évaluer la validité.
  • Des simulations de "rares auto-interactions" (rSIDM) isotropes pour tester l'influence de la dépendance angulaire de la section efficace.
• Application Observationnelle : Les auteurs projettent les résultats de leurs simulations pour dériver des profils de densité de surface et des masses enfermées, en comparant cela à la distribution de masse inférée pour le perturbateur de JVAS B1938+666 par Vegetti et al. (2026).

Résultats Clés

1. Changement Qualitatif de l'Évolution : La dissipation ne se contente pas d'accélérer l'effondrement ; elle modifie qualitativement l'évolution. Dans le cas élastique, la conduction thermique finit par inverser sa direction (vers l'extérieur) à mesure que le cœur devient isotherme, provoquant l'effondrement. Dans le cas dissipatif, un refroidissement central suffisamment fort empêche la formation d'un cœur isotherme. Le gradient de dispersion de vitesse reste positif au centre, maintenant la conduction thermique dirigée vers l'intérieur tout au long de l'évolution.
2. Afflux de Masse vs Sortie de Masse : Dans les halos élastiques, le halo extérieur est chauffé par la conduction, ce qui entraîne une sortie de masse. Dans les halos dissipatifs, les régions extérieures se refroidissent efficacement (la dissipation dominant la conduction vers l'extérieur), ce qui conduit à un afflux de masse. Cela résulte en une région plus large d'afflux de masse et une indentation moins prononcée entre le cœur et le halo extérieur dans le profil de densité final.
3. Dépendance aux Paramètres :
   • Dissipation : Augmenter le taux de refroidissement (contrôlé par $r_{\text{diss}}$) accélère considérablement l'effondrement (par près de deux ordres de grandeur pour les paramètres testés) et produit des cœurs plus denses et moins étendus.
   • Conduction : L'effet de la section efficace de transfert est non monotone. À de faibles sections efficaces (long libre parcours moyen), l'augmentation de la conduction accélère l'effondrement. À des sections efficaces élevées (court libre parcours moyen), l'augmentation de la conduction supprime l'efficacité du transport de chaleur, retardant l'effondrement.
   • Dépendance Angulaire : Les auteurs constatent que l'évolution des halos isolés est largement insensible à la dépendance angulaire de la section efficace (orientée vers l'avant vs isotrope), à condition que les taux de refroidissement et de conduction soient appariés.
4. Application à JVAS B1938+666 : Les auteurs démontrent que les halos de SIDM faiblement dissipatifs peuvent reproduire le profil de masse inféré du perturbateur de JVAS B1938+666. Comparé au cas élastique, la SIDM dissipative peut expliquer l'observation avec :
   • Des temps d'évolution nettement plus courts (permettant l'effondrement avant $z=0,881$ avec des sections efficaces modérées).
   • Ou, équivalemment, des sections efficaces d'auto-interaction plus petites pour un même temps d'évolution.
   • Les sections efficaces requises ($\sigma_T/m_\chi \sim 10^2 \text{ cm}^2 \text{ g}^{-1}$) sont cohérentes avec les contraintes pour des halos de masse $M_{200c} \lesssim 10^7 M_\odot$.

Signification et Revendications
L'article affirme ouvrir une "nouvelle voie pour connecter la structure des halos et les objets compacts récemment rapportés à la microphysique du secteur sombre". En étendant le formalisme fSIDM pour inclure la dissipation, les auteurs fournissent un cadre pour tester des modèles présentant des médiateurs légers et des états liés. La principale signification réside dans le fait que la dissipation peut résoudre la tension entre l'objet compact observé dans JVAS B1938+666 et les grandes sections efficaces requises par les modèles élastiques. Les auteurs notent modestement que leurs conclusions reposent sur des conditions initiales spécifiques (halos NFW isolés) et que des simulations cosmologiques ainsi que des simulations s'étendant plus profondément dans la phase d'effondrement (au-delà de la rupture numérique) sont nécessaires pour contraindre pleinement les propriétés dissipatives. Ils reconnaissent également que la dépendance angulaire de la section efficace peut devenir plus pertinente pour les satellites et les fusions, un régime qui n'est pas couvert par cette étude de halo isolé.