Joint Estimation of Properties of the Lunar Subsurface and Galactic Foregrounds with LuSEE-Night

Titre

Estimation conjointe des propriétés du sous-sol lunaire et des avant-plans galactiques avec LuSEE-Night

1. Problématique

Le projet LuSEE-Night (Lunar Surface Electromagnetics Experiment) est un télescope radio basse fréquence qui sera déployé sur la face cachée de la Lune en 2027. Son objectif principal est de mesurer le signal de l'hydrogène neutre (21 cm) provenant de l'Âge Sombre et de l'Aube Cosmique, un signal extrêmement faible (de l'ordre de quelques dizaines de mK).

Le défi majeur pour la calibration de cet instrument réside dans les propriétés diélectriques inconnues du sous-sol lunaire au site d'atterrissage. Le sous-sol lunaire agit comme un milieu diélectrique absorbant et réfléchissant qui modifie le diagramme de rayonnement (le « faisceau ») de l'antenne.

• Impact : Les réflexions et l'absorption par le sous-sol peuvent altérer le faisceau primaire de 10 à 20 %, affectant l'amplitude, la position et la largeur de la résonance de l'antenne (autour de 30 MHz).
• Complexité : Il est difficile de distinguer les variations spectrales causées par le sous-sol lunaire de celles provenant des avant-plans galactiques (le bruit de fond radio de la Voie lactée), qui sont intrinsèquement plus brillants que le signal cosmologique recherché.

2. Méthodologie

Les auteurs ont développé une chaîne de traitement basée sur l'inférence bayésienne pour estimer simultanément les paramètres du sous-sol et ceux de l'avant-plan galactique.

A. Simulation du faisceau et du sous-sol

• Outil : Utilisation du logiciel Ansys HFSS (Méthode des Éléments Finis) pour simuler le champ lointain des antennes de LuSEE-Night (quatre monopôles en dipôle croisé).
• Modèle du sous-sol : Le sous-sol est modélisé comme un milieu stratifié à deux couches avec des constantes diélectriques ($\epsilon_{r,top}$, $\epsilon_{r,bottom}$) et une profondeur d'interface ($L$).
• Grille de paramètres : 216 configurations différentes du sous-sol ont été simulées (basées sur les données inversées du site d'atterrissage d'Apollo 17) pour créer une grille de paramètres ($6 \times 6 \times 6$).

B. Modélisation de l'avant-plan galactique

• Le spectre de brillance de la Voie lactée est modélisé par une loi de puissance avec une courbure logarithmique (équation 2), utilisant la carte de référence ULSA (Ultralong-wavelength Sky Model).
• Les paramètres libres sont l'indice spectral ($\beta$), la courbure ($\gamma$), et des facteurs d'amplitude et d'offset ($k_1, k_2$).

C. Génération des données simulées

• Les mesures de densité spectrale de puissance (PSD) de tension sont simulées en convoluant le diagramme de rayonnement de l'antenne (dépendant du sous-sol) avec la carte du ciel (dépendant de l'avant-plan).
• Le bruit de fond lunaire est supposé calibré (température de 0 K pour simplifier), et un bruit gaussien de 5 % est ajouté aux données simulées.

D. Pipeline d'inférence

• Émulateur : Pour éviter de lancer des simulations HFSS coûteuses à chaque itération, un émulateur basé sur des fonctions de base radiales (RBF) est construit pour interpoler les résultats entre les points de la grille de simulation.
• Algorithme : Un algorithme de Monte Carlo Préconditionné (PMC) est utilisé pour maximiser la fonction de vraisemblance et récupérer les paramètres a posteriori.
• Objectif : Estimer conjointement les 7 paramètres : 3 pour le sous-sol ($\epsilon_{r,top}, \epsilon_{r,bottom}, L$) et 4 pour l'avant-plan ($\beta, \gamma, k_1, k_2$).

3. Contributions Clés

1. Preuve de concept de l'estimation conjointe : L'article démontre qu'il est possible de séparer les effets du sous-sol lunaire de ceux de l'avant-plan galactique, car leurs signatures spectrales sont orthogonales.
   • Les effets du sous-sol sont localisés autour de la résonance de l'antenne (~30 MHz) et affectent la structure fine du spectre.
   • Les effets de l'avant-plan galactique produisent des variations spectrales lisses et à large bande (pentes et courbures).
2. Analyse des modes propres : L'analyse par décomposition en valeurs singulières (SVD) confirme que les modes de variation dominants pour le sous-sol et le ciel sont distincts, justifiant l'approche d'inférence conjointe.
3. Modélisation réaliste de l'impédance : Contrairement à des approches simplistes utilisant uniquement la température d'antenne, cette étude intègre l'effet du sous-sol sur l'impédance de l'antenne, qui est directement liée à la tension mesurée.

4. Résultats

• Convergence : L'algorithme PMC converge vers une distribution a posteriori qui englobe correctement les paramètres d'entrée « vrais » utilisés pour générer les données simulées.
• Précision : Les paramètres du sous-sol et de l'avant-plan sont récupérés avec une précision d'environ 1 % à 5 % (écart-type faible).
  • Exemple : $\epsilon_{r,top}$ récupéré à $3.8012 \pm 0.0016$ (vrai : 3.8).
  • Exemple : $L$ récupéré à $2.0035 \pm 0.0048$ m (vrai : 2.0 m).
• Corrélations : Bien que des corrélations existent entre certains paramètres (ex: entre les constantes diélectriques et la profondeur), l'information spectrale indépendante permet de contraindre l'ensemble du modèle.
• Limites de l'émulateur : Une analyse de sensibilité (Annexe A) montre que si les paramètres réels ne se trouvent pas exactement sur la grille de simulation ($6 \times 6 \times 6$), les erreurs d'interpolation de l'émulateur augmentent, en particulier près de la résonance à 30 MHz. Cela indique que la grille actuelle est suffisante pour une démonstration de principe mais devra être densifiée pour une application réelle.

5. Signification et Perspectives

• Calibration instrumentale : Ce travail est crucial pour la mission LuSEE-Night. Il montre que l'instrument peut non seulement mesurer le signal cosmologique, mais aussi caractériser in situ les propriétés diélectriques du régolithe lunaire à son site d'atterrissage, ce qui est essentiel pour une calibration précise.
• Robustesse : Malgré un modèle de sous-sol simplifié (2 couches) et un avant-plan idéalisé, les résultats suggèrent que la méthode est robuste car les signatures spectrales sont fondamentalement différentes.
• Travaux futurs :
  • Développer des modèles de sous-sol plus complexes (multi-couches) et des stratégies d'échantillonnage plus denses (ex: Latin Hypercube Sampling) pour améliorer la précision de l'émulateur.
  • Intégrer des modèles analytiques pour le diagramme de rayonnement.
  • Prendre en compte les bruits instrumentaux réels et les effets systématiques non idéalisés.

En conclusion, cette étude valide la faisabilité d'une approche d'inférence conjointe pour démêler les effets géophysiques locaux (Lune) des effets astrophysiques globaux (Galaxie), ouvrant la voie à une exploitation optimale des données de LuSEE-Night pour sonder l'Univers primordial.