Multi-Epoch X-Ray Detection of SLSN-I 2018bsz: Constraints on the Powering Mechanism and Ejecta Structure

Résumé technique : Détection de rayons X multi-époques de la SLSN-I 2018bsz

Problématique et motivation
Les supernovas superlumineuses (SLSN) sont de rares explosions stellaires dont les énergies sont environ 100 fois supérieures à celles des supernovas de l'effondrement de cœur ordinaires. Les mécanismes alimentant ces événements font l'objet de débats, les candidats principaux incluant une forte interaction entre l'éjecta et le milieu circumstellaire (CSM), ou l'injection d'énergie par une étoile à neutrons rapidement rotative et hautement magnétisée (magnétar). Bien que la modélisation optique ait réussi à reproduire les courbes de lumière de nombreuses SLSN-I à l'aide de modèles de magnétar, distinguer les sources de puissance est difficile car différents mécanismes peuvent produire des signatures optiques similaires. Les observations de rayons X offrent une sonde critique pour lever cette dégénérescence ; cependant, les détections sont rares, particulièrement aux époques tardives lorsque l'absorption par l'éjecta diminue et que l'émission du magnétar pourrait devenir visible. Avant cette étude, seule une SLSN-I (SCP 06F6) avait été définitivement détectée en rayons X.

Méthodologie
Les auteurs présentent une étude spectroscopique de rayons X multi-époques de la SN 2018bsz, la plus proche SLSN-I dépouillée connue (distance ~111 Mpc). Le jeu de données comprend quatre observations de Chandra (ACIS-S/S3) réalisées entre 87 et 304 jours après l'explosion, et une observation de XMM-Newton (EPIC-pn/MOS) réalisée à 1253 jours après l'explosion.

La réduction des données a suivi les procédures standards utilisant CIAO (v4.13) pour Chandra et SAS (v22.1.0) pour XMM. L'analyse spectrale a été effectuée avec XSPEC (v12.15.0c) en utilisant la statistique C en raison des faibles taux de comptage.

• Identification de la source : Les images de Chandra ont révélé la source SN et trois sources de rayons X plus faibles et proches (S1–S3) à l'intérieur de la fonction d'étalement du point (PSF) de XMM.
• Modélisation spectrale : Les spectres ont été ajustés avec une absorption galactique ($N_H = 2,02 \times 10^{21}$ cm$^{-2}$). Les modèles incluaient des lois de puissance (représentant le spin-down du magnétar ou une interaction non thermique) et des modèles de plasma thermique (mekal, représentant l'interaction du choc avec le CSM).
• Gestion de la contamination : Pour l'observation tardive de XMM, les auteurs ont modélisé la contribution des sources proches (S1–S3) en utilisant un « anneau de contamination » dérivé des données de Chandra. Ils ont ajusté le spectre de XMM avec un modèle composite : tbabs(powCont.A + pow) ou tbabs(powCont.A + mekal), où la composante de contamination était fixée sur la base des résultats de Chandra.
• Contraintes sur la source de puissance : Les auteurs ont testé deux scénarios :
  1. Magnétar : Un ajustement par la méthode de Monte Carlo par chaînes de Markov (MCMC) a été effectué pour contraindre la période de rotation initiale ($P_0$) et le champ magnétique ($B$) en supposant un modèle de spin-down dipolaire. Ils ont comparé ces résultats aux ajustements de la courbe de lumière optique de Gomez et al. (2024) et considéré les effets de l'absorption par l'éjecta et du breakout d'ionisation.
  2. Interaction CSM : Les taux de perte de masse ($\dot{M}$) ont été estimés en utilisant la relation luminosité-perte de masse de Chevalier & Fransson (2017).

Résultats clés

1. Détection : La SN 2018bsz a été détectée avec une signification >3$\sigma$ lors des quatre époques de Chandra. Une détection ténue a été faite dans l'observation XMM à 1253 jours, bien que celle-ci soit sensible au niveau de contamination des sources proches. Cela fait de SN 2018bsz la deuxième SLSN-I détectée en rayons X et la troisième SLSN détectée en rayons X au total.
2. Évolution spectrale :
   • Époques précoces (Chandra) : Les spectres sont bien décrits par des lois de puissance dures avec des indices de photons $\Gamma \approx 0,9–1,9$ (le meilleur ajustement étant $\Gamma = 1,32$ pour un ajustement simultané).
   • Époque tardive (XMM) : Le spectre s'adoucit de manière significative. À 1253 jours, les données sont cohérentes avec une loi de puissance douce ($\Gamma < 4,0$) ou un modèle de plasma thermique avec $kT \approx 0,20$ keV.
3. Évolution de la luminosité : La luminosité en rayons X reste relativement plate au cours du temps ($L_{0,5-8 \text{ keV}} \approx 1,1 \times 10^{40}$ erg s$^{-1}$), montrant peu de déclin entre 87 et 304 jours.
4. Contraintes du magnétar :
   • Les ajustements MCMC aux données de rayons X (en supposant une absorption par l'éjecta négligeable) produisent un champ magnétique faible ($\log B \approx 13,1$) et une période de rotation initiale longue ($P_0 \approx 16$ ms).
   • Ces paramètres sont incohérents avec les paramètres de « meilleur ajustement » du magnétar dérivés des courbes de lumière optiques par Gomez et al. (2024) ($P_0 \approx 0,9$ ms, $\log B \approx 14,4$).
   • En utilisant les paramètres dérivés de l'optique ($B_{opt}, P_{opt}$), le modèle de magnétar sous-estime les luminosités de rayons X observées à toutes les époques, sauf pour les plus précoces, et ne parvient pas à reproduire la courbe de lumière plate.
   • La prise en compte de l'absorption par l'éjecta accentue cette divergence. Les modèles de breakout d'ionisation suggèrent que même avec des paramètres favorables, l'émission du magnétar ne s'échapperait de l'éjecta qu'entre ~2 et 72 ans après l'explosion, bien au-delà de la fenêtre d'observation.
5. Contraintes de l'interaction CSM :
   • L'adoucissement spectral et l'évolution plate de la luminosité sont cohérents avec une interaction initiale forte éjecta-CSM qui s'affaiblit avec le temps.
   • Les taux de perte de masse dérivés sont de l'ordre de $5,4 \times 10^{-3}$ à $1,4 \times 10^{-2}$ $M_{\odot}$ an$^{-1}$. Ces taux dépassent ceux des supernovas dépouillées typiques et de la plupart des SLSN-I, se rapprochant davantage des SLSN-II ou des supernovas dépouillées qui évoluent vers des systèmes en interaction via une perte de masse éruptive.
   • L'époque estimée de l'épisode de perte de masse du progéniteur est de 5 à 18 ans avant l'explosion, ce qui est cohérent avec une structure de CSM de type coquille ou disque proposée dans des études optiques antérieures.

Signification et conclusions
L'article conclut que l'émission de rayons X de la SN 2018bsz est peu probable d'être alimentée par un moteur central de magnétar à milliseconde standard, car les paramètres requis pour ajuster les données de rayons X contredisent ceux dérivés de la modélisation optique, et que les échelles de temps de breakout d'ionisation sont trop longues.

Au lieu de cela, les observations favorisent un scénario où l'émission de rayons X est dominée par l'interaction entre l'éjecta et un CSM dense et asymétrique (de type disque). Cette interprétation est soutenue par :

• L'évolution temporelle plate de la luminosité en rayons X.
• L'évolution spectrale de dur à mou.
• Les taux de perte de masse élevés cohérents avec une activité éruptive du progéniteur.
• Les résultats optiques précédents indiquant une émission d'hydrogène et une interaction avec un CSM de type disque.

Les auteurs avancent que la SN 2018bsz représente une sous-classe distincte de SLSN-I où l'interaction avec le CSM est cruciale pour la forte émission, et où un magnétar, s'il est présent, ne domine pas la production de rayonnement à haute énergie. Cette étude souligne la nécessité des observations multi-longueurs d'onde pour distinguer les mécanismes de puissance dans les supernovas extrêmes.