Mirage Sources and Large TeV Halo-Pulsar Offsets: Exploring the Parameter Space

Résumé Technique : Sources Mirages et Décalages de Halo Pulsar à l'échelle du TeV

Énoncé du Problème
Des observations récentes effectuées par le Large High Altitude Air Shower Observatory (LHAASO) et HAWC ont identifié des émissions de rayons gamma TeV étendues (« halos de pulsars ») autour de pulsars d'âge moyen. Un défi important dans l'interprétation de ces sources est l'observation fréquente de décalages spatiaux importants entre la position du pulsar et le centroïde de l'émission TeV. Les modèles traditionnels supposant une diffusion isotrope des paires électron-positron ($e^\pm$) dans le milieu interstellaire (MIS) ne peuvent pas expliquer naturellement ces décalages, car la diffusion isotrope n'est une approximation valide que lorsque l'échelle physique dépasse plusieurs longueurs de cohérence magnétique. Par conséquent, pour les électrons d'ultra-haute énergie ($>100$ TeV) présentant des taux de refroidissement rapides, l'échelle de propagation peut être plus petite que la longueur de cohérence du champ magnétique turbulent du MIS. En conséquence, le transport des particules peut être hautement asymétrique, suivant les lignes de champ magnétique locales plutôt que de diffuser uniformément. Cet article étudie si une telle propagation asymétrique peut générer des sources « mirages » — des pics d'émission apparents déplacés par rapport à l'accélérateur — et des décalages spatiaux significatifs.

Méthodologie
Les auteurs utilisent une simulation de particules de test sur GPU, basée sur les principes fondamentaux, pour suivre les trajectoires d'électrons de 100 TeV injectés dans un champ magnétique turbulent.

• Génération du Champ Magnétique : Le champ magnétique turbulent est modélisé comme une superposition de quatre champs indépendants à spectre de Kolmogorov ($n_{in} = -11/3$) générés via une méthode dans l'espace de Fourier. Cette approche multi-échelle (utilisant des espacements de grille allant de $4 \times 10^{-9}$ à $5 \times 10^{-4}$ kpc) garantit un intervalle inertiel réaliste et des champs sans divergence ($\nabla \cdot B = 0$). Le champ total inclut des composantes régulières ($B_r$) et turbulentes ($B_t$).
• Propagation des Particules : Les électrons sont propagés à l'aide d'un algorithme de pousseur de Boris pour résoudre l'équation du mouvement relativiste sous l'effet de la force de Lorentz. La simulation prend en compte les pertes d'énergie radiatives via le rayonnement synchrotron et la diffusion Compton inverse (IC) sur le fond diffus cosmologique (CMB) et le champ de rayonnement interstellaire (ISRF).
• Échantillonnage Statistique : Pour échantillonner efficacement le spectre d'énergie initial ($dN/dE \propto E^{-\alpha} \exp(-E/E_{cut})$), les auteurs utilisent une approche de macro-particules pondérées. 8 192 macro-particules sont injectées, chacune représentant un intervalle d'énergie logarithmique avec un poids proportionnel au nombre physique d'électrons dans cet intervalle.
• Synthèse Observationnelle : Les distributions d'électrons simulées sont converties en cartes de rayons gamma synthétiques (bande 25–40 TeV) via une matrice de diffusion Compton inverse (IC). Ces cartes sont convoluées avec une fonction de réponse impulsionnelle (PSF) gaussienne de $\sigma_{PSF} = 0,43^\circ$ (imitant LHAASO-KM2A) et ajoutées à un modèle de bruit de fond de Poisson.
• Identification des Sources : Un cadre de comparaison de modèles bayésien utilisant des chaînes de Markov Monte Carlo (MCMC) est appliqué. Des modèles avec un nombre croissant de sources gaussiennes sont ajustés aux données synthétiques. Une source est considérée comme « identifiée » si la statistique de test (TS) augmente de plus de 25 ($\Delta \log L > 12,5$) lors de l'ajout d'une composante.

Contributions Clés et Résultats
L'étude explore systématiquement l'espace des paramètres, faisant varier la distance du pulsar ($d$), l'indice spectral d'injection des électrons ($\alpha$), l'intensité du champ magnétique total ($B$), la longueur de cohérence de la turbulence ($L_c$) et le rapport entre le champ régulier et le champ turbulent ($B_r/B_t$).

1. Formation de Sources Mirages : Les simulations confirment qu'une propagation asymétrique peut créer plusieurs sources de rayons gamma distinctes (« sources mirages ») à partir d'un seul accélérateur. Cela se produit lorsque les électrons se propagent le long des lignes de champ magnétique qui, en raison d'effets de projection, semblent provenir d'un emplacement éloigné du point d'injection.
2. Multiplicité des Sources ($N$) :
   • Le nombre de sources identifiées est le plus fortement influencé par l'intensité du champ magnétique ($B$). Un $B$ élevé entraîne un refroidissement plus rapide, confinant les électrons à de plus petites échelles et réduisant la probabilité de former plusieurs sources distinctes ( $N$ plus faible).
   • Le nombre de sources est corrélé négativement avec la distance ($d$) et la longueur de cohérence ($L_c$), et positivement avec $B^2$.
3. Séparation des Sources et Rapport de Flux :
   • Séparation : La séparation maximale entre les sources identifiées est positivement corrélée avec $L_c$ et négativement corrélée avec $B^2$.
   • Rapport de Flux ($\chi$) : Le rapport entre le flux de la source principale et celui de la source mirage diminue avec l'augmentation de $B^2$. Dans les champs magnétiques forts, la source mirage peut parfois être plus brillante que la source principale.
   • Champs Réguliers vs Turbulents : La distance de séparation présente une dépendance non monotone vis-à-vis de $B_r/B_t$. La séparation maximale se produit à $B_r/B_t \approx 1$, où le champ régulier fournit une guidance directionnelle tandis que la turbulence permet une diffusion suffisante pour créer des décalages. La turbulence pure ou les champs réguliers dominants entraînent de plus petites séparations.
4. Décalages Spatiaux :
   • Les décalages entre le pulsar et la source la plus proche identifiée sont maximisés lorsque $B_r/B_t \approx 1$ et que $L_c$ est grand.
   • Dans les régions avec de grandes longueurs de cohérence (par exemple, les régions inter-bras ou le halo galactique), les décalages peuvent atteindre jusqu'à 50 pc.

Signification et Revendications
L'article prétend fournir une explication quantitative pour les sources LHAASO déroutantes où les halos TeV apparaissent spatialement déconnectés de leurs pulsars associés. Plus précisément, les auteurs suggèrent que des sources telles que 1LHAASO J0216+4237u, J0212+4254u et J0206+4302u pourraient être des halos mirages alimentés par un seul pulsar milliseconde (PSR J0218+4232), plutôt que des accélérateurs distincts.

L'étude souligne que l'approximation de la diffusion isotrope est insuffisante pour modéliser le transport des particules près des accélérateurs dans le MIS, particulièrement pour les leptons de haute énergie dont les temps de refroidissement sont plus courts que le temps de cohérence magnétique. En démontrant que les effets de projection dans les champs magnétiques turbulents peuvent naturellement produire de grands décalages et des morphologies de sources multiples, ce travail offre un mécanisme pour réconcilier le ciel TeV observé avec les modèles standards de nébuleuses de vent de pulsar. Les auteurs notent que bien que leur turbulence synthétique suppose l'isotropie (contrairement à la turbulence Alfvénique anisotrope de la MHD), elle capture qualitativement la nature filamentaire à grande échelle de la turbulence du MIS responsable de ces émissions asymétriques.