Tentative detection of circularly polarized bursty radio emissions from the HD 189733 exoplanetary system using NenuFAR beamformed observations

Résumé technique : Détection provisoire d'émissions radio par rafales à polarisation circulaire du système exoplanétaire HD 189733

Problématique et contexte
La détection des champs magnétiques exoplanétaires demeure un défi majeur de la planétologie, l'émission radio aurorale via l'instabilité du maser cyclotronique (CMI) étant considérée comme la méthode de détection la plus prometteuse. Alors que des décennies d'observations au sol n'ont produit que des non-détections sans ambiguïté, des détections provisoires récentes dans les systèmes $\tau$ Boötis et HD 189733 ont renouvelé l'intérêt. Plus précisément, Zhang et al. (2025, ci-après Z25) ont rapporté une détection à 6$\sigma$ d'un sursaut radio hautement polarisé circulairement à 50 MHz provenant du système HD 189733 en utilisant des observations d'imagerie NenuFAR. Cependant, l'origine de cette émission reste ambiguë, avec quatre scénarios potentiels proposés : l'interaction étoile-planète sub-Alfvénique (SPI), l'émission aurorale magnétosphère-vent, la CMI stellaire, ou l'émission de plasma stellaire issue de éruptions. L'étude Z25 a utilisé des données d'imagerie mais manquait d'une vérification indépendante utilisant des données de formation de faisceaux (beamforming), qui offrent une sensibilité plus élevée aux sursauts atténués.

Méthodologie
Cette étude analyse les observations par formation de faisceaux de NenuFAR du système HD 189733, réalisées simultanément avec les données d'imagerie ayant produit la détection Z25 (28 septembre 2023, 19:31–22:47 UT). Le jeu de données comprend un « ON-beam » ciblant le système et trois « OFF-beams » ciblant des régions dépourvues de sources pour caractériser les fluctuations ionosphériques terrestres, les interférences radio (RFI) et les biais systématiques instrumentaux.

L'analyse a utilisé la pipeline BOREALIS (BeamfOrmed Radio Emission AnaLysIS), précédemment validée sur les données de Jupiter de LOFAR et les données de pulsars de NenuFAR. Le flux de travail comprenait :

1. Prétraitement : Génération de données de Niveau 2 (L2) calculant les paramètres de Stokes (I, Q, U, V), application de l'atténuation des RFI et correction des variations de gain instrumental.
2. Filtrage : Filtrage passe-haut des spectres dynamiques de Stokes V (polarisation circulaire) pour éliminer les systématiques à grande échelle (ex. changements de ligne de base et artefacts de lobes de réseau).
3. Analyse statistique : Calcul des observables de sursaut ($Q2$, $Q3$, $Q4$) comparant l'ON-beam aux OFF-beams. L'observable $Q2$ fournit une comparaison de séries temporelles, tandis que $Q3$ et $Q4$ offrent des mesures statistiques de la signification du sursaut par rapport à une courbe de référence de 10 000 tirages de bruit gaussien.
4. Paramètres de recherche : La recherche a couvert plusieurs plages de fréquences (21–58 MHz, avec un accent spécifique sur 27–40 MHz et la plage 47–52 MHz de Z25) et des résolutions temporelles ($\Delta\tau$ de 0,1, 1 et 8 secondes).

Résultats clés
L'analyse a produit une détection provisoire d'une émission par rafales à polarisation circulaire distincte du résultat d'imagerie Z25 tant en temps qu'en fréquence :

• Signification de la détection : Un signal fort a été identifié dans la plage 27–40 MHz entre 19:40 et 19:56 UT (environ 1 heure avant le sursaut d'imagerie Z25). La détection a atteint une signification de 7,4$\sigma$ dans la statistique $Q4f$ et de 10$\sigma$ dans la statistique $Q3f$ par rapport au bruit gaussien.
• Caractéristiques du signal : L'émission était polarisée à gauche (analyse V) avec une durée d'environ 16 minutes, composée de sursauts individuels d'environ 1 seconde. La densité de flux estimée est de ~3,9 Jy, ce qui correspond à une température de brillance ($T_B$) de $4–6 \times 10^{16}$ K.
• Comparaison avec Z25 : Contrairement à la détection par imagerie de Z25 (polarisation droite, 47,6–52,1 MHz, durée de 96s, 1,5 Jy), cette détection par formation de faisceaux s'est produite à une fréquence plus basse, plus tôt dans le temps, et présente des caractéristiques de polarisation et de durée différentes.
• Analyse du bruit : Les auteurs ont identifié un bruit corrélé excessif dans les données, visible sous la forme d'un nuage de points diagonal dans les graphiques de diagnostic (Figure 2A–B) et un comportement non gaussien dans les résidus des OFF-beams. Bien que les statistiques $Q4$ soient conçues pour ignorer les données diagonales corrélées, la présence de ce bruit empêche une affirmation concluante sur l'origine astrophysique.

Interprétation et contraintes
En supposant une origine astrophysique, les auteurs évaluent les quatre scénarios proposés :

1. Émission de plasma stellaire : Exclue comme mécanisme primaire car la température de brillance observée ($10^{16}$ K) dépasse largement la limite de saturation théorique pour l'émission de plasma ($\sim 10^{12}$ K).
2. SPI sub-Alfvénique et CMI stellaire : La plage de fréquences observée implique un rapport fréquence plasma/cyclotron ($f_{pe}/f_{ce}$) de 2,1–2,3, ce qui est trop élevé pour supprimer le mécanisme CMI, à moins que le champ magnétique stellaire ne soit exceptionnellement fort (150–200 G).
3. Interaction vent-magnétosphère : Ce scénario est jugé le plus cohérent avec les observations. La durée du sursaut d'environ 1 seconde ressemble aux S-bursts de Jupiter, et le timing s'aligne avec les modèles de transfert radiatif suggérant que l'émission radio planétaire s'échappe principalement près du transit primaire. Ce scénario implique une contrainte de champ magnétique planétaire de $\sim$12–18 G.

Signification et revendications
L'article revendique une détection provisoire d'une émission radio par rafales à polarisation circulaire provenant du système HD 189733. Les auteurs soulignent que, bien que le signal soit statistiquement significatif par rapport au bruit gaussien, la présence d'un bruit corrélé excessif (potentiellement des RFI de faible niveau ou des systématiques propres à NenuFAR) empêche une confirmation définitive.

L'étude met en évidence la puissance complémentaire des observations simultanées par formation de faisceaux et par imagerie. Les données par formation de faisceaux ont détecté un sursaut à un temps et une fréquence différents des données d'imagerie, suggérant que les émissions radio exoplanétaires peuvent être hautement variables et qu'une observation à époque unique peut manquer une activité significative. Les auteurs concluent que des observations radio à basse fréquence supplémentaires sont nécessaires pour confirmer la nature astrophysique du signal, rechercher la périodicité et distinguer définitivement les origines planétaires et stellaires.