Relativistic jets from millisecond proto-magnetars

Résumé Technique : Jets Relativistes issus de Proto-Magnétars à Milliseconde

Énoncé du Problème
Les proto-magnétars à milliseconde — des étoiles à neutrons tournant rapidement ($P \sim 1$ ms) et fortement magnétisées ($B \gtrsim 10^{15}$ G) formées lors de supernovas à effondrement de cœur, de fusions d'étoiles à neutrons ou d'effondrement induit par accrétion (AIC) de naines blanches — sont depuis longtemps proposés comme moteurs centraux des sursauts gamma (GRB). Un défi théorique critique demeure : durant les premières secondes suivant leur naissance, ces objets émettent des flux de neutrinos intenses qui chauffent la surface et pilotent des vents riches en baryons. Ce « vent de neutrinos » charge l'écoulement en matière, ce qui peut supprimer la magnétisation ($\sigma$) et limiter le facteur de Lorentz de groupe ($\Gamma_\infty$) atteignable. Les modèles semi-analytiques précédents, calibrés sur des simulations unidimensionnelles et supposant des géométries de type monopole divisé à moyenne angulaire, prédisaient un $\sigma \lesssim 1$ aux premiers instants, suggérant que les proto-magnétars ne pourraient pas lancer les jets ultra-relativistes ($\Gamma \gtrsim 10-100$) nécessaires pour alimenter les sursauts gamma à durée courte. La question de savoir si des effets multidimensionnels pouvaient surmonter cette charge en baryons pour produire des jets relativistes dans les secondes suivant la formation n'avait pas été pleinement traitée.

Méthodologie
Les auteurs présentent les premières simulations en magnétohydrodynamique générale relativiste (GRMHD) en trois dimensions de vents de proto-magnétars chauffés par neutrinos. L'étude utilise une version modifiée du code open-source GRHydro (faisant partie de l'Einstein Toolkit) employant un schéma à volumes finis avec une reconstruction WENO-Z et un solveur de Riemann HLLE approximatif. Caractéristiques numériques clés :

• Transport de Neutrinos : Une approximation à un moment ($M_0$) de l'équation de Boltzmann en relativité générale est utilisée pour modéliser de manière auto-cohérente le chauffage et le refroidissement par neutrinos, incluant les processus $\beta$ à courant chargé, l'annihilation de paires et la diffusion.
• Évolution du Champ Magnétique : Le champ magnétique est évolué via le potentiel électromagnétique 4-vecteur dans la jauge de Lorenz généralisée, en utilisant un schéma de transport de contraintes (constraint transport) de type upwind pour maintenir la contrainte solénoïdale avec une précision machine.
• Conditions Initiales : Les modèles supposent une étoile à neutrons de $1,4 M_\odot$ avec une température centrale de $20$ MeV. Les simulations couvrent une gamme de périodes de rotation initiales ($P \approx 1 - 18$ ms) et de forces de champ magnétique de surface polaire ($B_P \approx 10^{13} - 3 \times 10^{15}$ G).
• Configuration : Les simulations utilisent une symétrie de réflexion par rapport au plan équatorial et utilisent un raffinement de maillage adaptatif (Berger-Oliger) avec la grille la plus interne résolvant la surface du proto-magnétar. Les modèles sont d'abord évolués vers un état stationnaire hydrodynamique piloté par les neutrinos avant d'initialiser un champ dipolaire poloidal.

Contributions Clés et Résultats
Les simulations révèlent que les effets multidimensionnels modifient fondamentalement la structure du vent par rapport aux modèles 1D ou à moyenne angulaire, particulièrement pour les modèles tournant rapidement ($P \approx 1$ ms) et fortement magnétisés ($B_P \approx 3 \times 10^{15}$ G) :

1. Stratification Angulaire de l'Écoulement : L'écoulement n'est pas uniforme mais se divise en trois régions distinctes :

   • Ceinture Équatoriale ($\theta \gtrsim 30^\circ$) : Les forces centrifuges renforcent significativement la perte de masse près de l'équateur de rotation. Cela crée un vent dense et sub-relativiste ($v \sim 0,1c$) avec une forte charge en baryons.
   • Jet Polaire ($\theta \lesssim 10^\circ$) : Les lignes de champ polaires subissent un renforcement minimal de la charge en baryons par force centrifuge. Par conséquent, cette région reste « propre » (faible charge en baryons) et hautement magnétisée.
   • Région de Transition : Latitudes intermédiaires où le vent accélère le plus efficacement.

2. Collimation Intrinsèque du Jet : Le vent équatorial dense et sub-relativiste agit comme un milieu de confinement, collimatant naturellement l'écoulement polaire moins chargé en baryons. Cette auto-collimation se produit à de petits rayons ($\lesssim 100$ km) sans nécessiter de confinement par des éjectas de supernova ou de fusion externes.

3. Atteinte de Conditions Ultra-Relativistes : Dans les modèles tournant le plus rapidement et étant les plus magnétisés (ex: modèle B3P1), le jet polaire atteint des valeurs de magnétisation de pointe $\sigma \sim 30 - 100$. Cela est suffisant pour atteindre des facteurs de Lorentz de groupe $\Gamma_\infty \sim 100$ sur de plus grandes échelles. L'angle d'ouverture polaire du jet est d'environ $8^\circ$, ce qui est cohérent avec les estimations observationnelles pour les GRB courts et longs.

4. Partition de l'Énergie : L'énergie totale du vent est répartie entre un jet polaire ultra-relativiste et un vent équatorial à large angle et sub-relativiste. Ce dernier porte la majeure partie de l'énergie cinétique et de la masse, alimentant potentiellement la supernova ou l'éjecta de kilonova, tandis que le premier alimente le GRB.

Signification et Revendications
L'article prétend démontrer que les proto-magnétars à milliseconde peuvent lancer des jets relativistes dans les secondes suivant leur formation, résolvant la tension de longue date entre la charge en baryons induite par les neutrinos et les exigences des moteurs de GRB. Les auteurs soutiennent que :

• La Multidimensionnalité est Essentielle : La charge en baryons n'empêche pas la formation de jets relativistes lorsque la dépendance latitudinale des forces centrifuges et du chauffage par neutrinos est prise en compte.
• Scénario de Moteur Unifié : Un seul moteur central peut simultanément alimenter un sursaut gamma à durée courte (via le jet polaire) et la supernova/kilonova associée (via le vent équatorial), expliquant la partition observée entre les énergies de GRB corrigées par le faisceau et les énergies cinétiques des éjectas.
• Cohérence Temporelle : Le mécanisme opère au cours des premières secondes de la vie du proto-magnétar, ce qui est cohérent avec la phase d'émission prompt des GRB courts.

Les auteurs notent que bien que leurs modèles supposent des champs dipolaires ordonnés et des luminosités de neutrinos fixes, les résultats suggèrent que l'émergence de jets relativistes est une conséquence générique de la rotation extrême et de la magnétisation dans les premières étapes de l'évolution des proto-étoiles à neutrons. Ils reconnaissent que les travaux futurs devront aborder les champs magnétiques générés par dynamo de manière auto-cohérente, le refroidissement dépendant du temps/le ralentissement (spin-down), ainsi que l'interaction avec des disques d'accrétion massifs entourant les rémanents de fusion.