Simulating observed point-symmetric core-collapse supernova morphologies with the jittering jets explosion mechanism

Résumé technique : Simulation des morphologies de supernovas à effondrement de cœur à symétrie ponctuelle observées avec le mécanisme d'explosion par jets oscillants (Jittering Jets)

Énoncé du problème
L'article traite du débat en cours concernant le mécanisme d'explosion primaire des supernovas à effondrement de cœur (CCSN). Deux cadres théoriques concurrents existent : le mécanisme piloté par les neutrinos (explosion retardée par les neutrinos) et le mécanisme d'explosion par jets oscillants (JJEM - jittering jets explosion mechanism). Bien que le mécanisme piloté par les neutrinos soit au centre de nombreuses simulations récentes, les auteurs soutiennent qu'il ne parvient pas à rendre compte des morphologies à symétrie ponctuelle observées dans environ 20 rémanents de supernova à effondrement de cœur (CCSNR). Ces morphologies, caractérisées par des paires de structures opposées (anneaux, lobes, buses) qui ne partagent pas un axe unique, sont proposées comme un observatif décisif pour distinguer les deux mécanismes. L'étude vise à démontrer que le JJEM, impliquant spécifiquement des paires multiples de jets, peut reproduire ces structures complexes à symétrie ponctuelle, même lorsque les jets sont étouffés profondément au sein du cœur stellaire.

Méthodologie
Les auteurs ont réalisé deux simulations hydrodynamiques tridimensionnelles (3D) à l'aide du code FLASH (version 4.8) avec un solveur hydrodynamique non fractionné (unsplit).

• Conditions initiales : Les simulations ont modélisé l'effondrement d'un cœur stellaire massif (progéniteur de 15 $M_{\odot}$) qui avait été dépouillé de ses enveloppes d'hydrogène et d'hélium (modèle Wolf–Rayet). La grille s'étendait à $\pm 6 \times 10^{10}$ cm, avec un raffinement de maillage adaptatif (AMR) fournissant une haute résolution (niveau 10) dans la région centrale ($r < 5 \times 10^8$ cm) pour résoudre l'injection des jets et l'interaction initiale avec le cœur.
• Injection de jets : Trois paires de jets inclinés et opposés ont été lancées dans le cœur en l'espace d'environ une seconde. Les jets ont été injectés à un rayon de $5 \times 10^7$ cm avec une vitesse initiale de $5 \times 10^4$ km s$^{-1}$ (Mach $\approx 6,5$).
• Cas de simulation :
  • Simulation D3 : Présentait trois paires de jets inégaux. Au sein de chaque paire, les jets opposés différaient par leur taux d'injection de masse, leur angle d'ouverture ou leur durée, entraînant des énergies cinétiques asymétriques.
  • Simulation E3 : Présentait trois paires de jets égaux, où les jets opposés d'une même paire possédaient des propriétés identiques, ne différant que par leur direction.
• Durée : Les simulations ont suivi les éjectas pendant plus de 10 secondes, dépassant largement le temps de rupture de l'onde de choc ($t \approx 4$ s) afin d'analyser l'évolution à long terme de la morphologie.
• Analyse : Les auteurs ont analysé les cartes de densité, les rapports d'énergie cinétique à thermique, la vorticité, les taux de croissance de l'instabilité de Rayleigh-Taylor (RTI) et les intégrales d'émission. Ils ont également généré des cartes de vitesse Doppler pour comparer les résultats avec les données observationnelles de CCSNR spécifiques.

Contributions clés et résultats
L'étude a réussi à reproduire plusieurs caractéristiques morphologiques de CCSNR observés que des paires de jets uniques ou des instabilités seules ne peuvent expliquer. Les principales conclusions sont les suivantes :

1. Morphologie multipolaire à symétrie ponctuelle : Malgré le fait que les jets aient été étouffés profondément à l'intérieur de l'étoile et aient cessé toute activité bien avant la rupture de l'onde de choc, les deux simulations ont produit une morphologie multipolaire (à symétrie ponctuelle) prononcée dans les éjectas externes. Dans la Simulation E3, cela s'est manifesté par trois paires distinctes de lobes alignés avec les trois axes des jets.
2. Filaments et amas : Les interactions jet-étoile ont généré des vortex et des instabilités (spécifiquement l'RTI), formant de petits amas denses et des filaments étroits. Bien que la formation d'amas individuels soit stochastique, les zones instables ont collectivement formé une morphologie à symétrie ponctuelle, ressemblant aux caractéristiques observées dans des rémanents tels que SNR 0540-69.3, la Vela et Cassiopeia A.
3. Structure de renflement (Simulation D3) : Le jet le plus énergétique de la Simulation D3 a créé un large renflement (blowout) de faible densité devant les éjectas, avec des filaments traînant derrière lui. Cette structure ressemble qualitativement au renflement sud de la CCSNR du Cygnus Loop. Les auteurs attribuent le bord brillant observé dans le Cygnus Loop réel à l'interaction de ce renflement avec le milieu ambiant, une caractéristique qui n'est pas pleinement capturée car la simulation ne s'est pas étendue dans le milieu circumstellaire.
4. Structure des éjectas internes (Simulation E3) : Les éjectas internes présentaient une paire d'anneaux et une paire de buses le long de deux axes de jets perpendiculaires. Cette configuration spécifique correspond qualitativement aux éjectats internes de SNR J0450.4-7050.
5. Écoulement bipolaire Doppler incliné : Une découverte critique fut la divergence entre l'axe morphologique et l'axe de vitesse décalé par l'effet Doppler. Le matériau le plus dense (métaux lourds) a été compressé en deux blocs entre les axes des jets, plutôt que le long de ceux-ci. Par conséquent, les cartes Doppler ont révélé un axe d'écoulement bipolaire Doppler fortement incliné par rapport aux axes de symétrie morphologique. Cela a permis de résoudre avec succès l'« énigme » des CCSNR comme W49B et SNR G292.0+1.8, où l'axe bipolaire Doppler observé est significativement incliné par rapport à l'axe de symétrie morphologique.

Signification et affirmations
L'article affirme que ces résultats apportent un soutien robuste au mécanisme d'explosion par jets oscillants (JJEM) en tant que mécanisme primaire des CCSN. Les auteurs soutiennent que l'identification de morphologies à symétrie ponctuelle dans une grande fraction de CCSNR bien résolus constitue une avancée majeure qui remet en question le mécanisme piloté par les neutrinos, qu'ils affirment incapable de produire naturellement ces structures à axes multiples spécifiques.

Plus précisément, l'étude soutient l'affirmation que :

• Le JJEM est capable de façonner des morphologies complexes à symétrie ponctuelle, même lorsque les jets sont étouffés profondément dans le progéniteur.
• La présence de trois ou plusieurs paires de jets énergétiques est nécessaire pour expliquer les caractéristiques structurelles spécifiques (anneaux, buses et axes Doppler inclinés) observées dans des rémanents comme J0450.4-7050, W49B et G292.0+1.8.
• L'accord qualitatif entre les simulations et les observations de ces quatre rémanents spécifiques renforce l'idée que le JJEM, plutôt que le mécanisme piloté par les neutrinos, déclenche l'explosion de la plupart, sinon de toutes, les CCSN.

Les auteurs maintiennent un ton modeste concernant les limites, notant que les simulations ne reproduisent pas toutes les propriétés des rémanents observés (par exemple, l'étroitesse exacte des filaments ou les bords brillants des renflements) en raison des limites de résolution numérique et de l'exclusion du refroidissement radiatif et des interactions avec le milieu ambiant dans la configuration actuelle. Cependant, ils affirment que la reproduction des principales caractéristiques qualitatives est suffisante pour valider le cadre du JJEM.