Search for a solar-bound axion halo using the Global Network of Optical Magnetometers for Exotic physics searches

Résumé Technique : Recherche d'un halo d'axions lié au Soleil à l'aide de GNOME

Problématique et Motivation
La nature de la matière noire demeure une question centrale et ouverte en astrophysique et en physique des particules. Les bosons ultra-légers, spécifiquement les axions ou les particules de type axion (ALP) avec des masses $m_a \ll 1$ eV, sont des candidats fortement motivés. Des travaux théoriques récents suggèrent que des corps astrophysiques denses, tels que le Soleil, peuvent capturer gravitationnellement des axions. Si les axions possèdent des auto-interactions quartiques, ces interactions peuvent fournir le mécanisme de dissipation d'énergie nécessaire pour permettre aux axions de se thermaliser dans le puits de potentiel gravitationnel fondamental du Soleil. Ce processus forme un « halo d'axions solaire », une configuration cohérente de l'état fondamental du champ d'axions centrée sur le Soleil.

Contrairement au halo de matière noire galactique standard, qui est virialisé et stochastique, un halo d'axions solaire dans l'état fondamental ($n=1, \ell=0$) présenterait un profil de densité distinct ($\rho \propto e^{-2r/R_\star}$) et une oscillation monochromatique cohérente à la fréquence de Compton de l'axion $\omega_a = m_a c^2 / \hbar$. La présence d'un tel halo impliquerait une surdensité significative d'axions à proximité de la Terre par rapport à la densité moyenne de la matière noire galactique, ce qui pourrait améliorer la sensibilité de détection. Le champ d'axion dans ce halo se couple aux spins atomiques via des interactions de gradient (linéaires ou quadratiques), générant un champ pseudo-magnétique ($B_p$) qui oscille selon un motif de modulation journalière spécifique dû à la rotation de la Terre et au mouvement orbital de la Terre par rapport au Soleil.

Méthodologie
Cette étude utilise les données du Réseau Mondial de Magnétomètres Optiques pour les recherches de physique exotique (GNOME), plus précisément les données de sa cinquième campagne scientifique (SR5), qui comprenait 69 jours de données provenant de 12 magnétomètres à blindage magnétique répartis géographiquement.

1. Modèle de Signal : L'analyse suppose que les axions sont dans l'état fondamental du potentiel gravitationnel solaire. Le champ d'axion résultant $a(r,t)$ génère un champ pseudo-magnétique $B_p$ par le biais de couplages de gradient linéaires ($H_l \propto S \cdot \nabla a$) et quadratiques ($H_q \propto S \cdot \nabla a^2$) aux spins des protons. Le signal se compose de deux composantes : une composante radiale (due au gradient spatial vers le Soleil) et une composante transversale (due au « vent d'axion » provenant du mouvement orbital de la Terre). La force relative de ces composantes dépend de la masse de l'axion (fréquence).
2. Prétraitement des Données : Les séries temporelles brutes de chaque station ont été filtrées pour éliminer le bruit de ligne électrique (50/60 Hz) et les dérives lentes (filtre passe-haut à 0,0003 Hz). Les données ont été regroupées par intervalles (binning) à un taux d'échantillonnage de 64 Hz. Des impulsions d'étalonnage horaires ont été utilisées pour corriger les dérives de gain.
3. Corrélation Croisée Décalée dans le Temps : Pour distinguer un signal global d'un bruit local, l'analyse emploie un pipeline de corrélation croisée. Comme la modulation journalière de l'amplitude du signal varie pour chaque station en fonction de sa localisation géographique et de son orientation, les données de différentes stations sont décalées temporellement ($\Delta t_j$) afin d'aligner les pics de signal prédits avant la corrélation croisée.
4. Analyse Spectrale : Les corrélations croisées par paires sont calculées dans des fenêtres de 20 minutes. Les signaux résultants sont convertis en densités spectrales de puissance (PSD), normalisés, puis combinés en un spectre au niveau du réseau pondéré.
5. Sélection et Ajustement des Candidats :
   • Étape 1 : Les fréquences où la PSD normalisée dépasse un seuil ($P_\star = 3,35$) chaque jour de la campagne sont identifiées comme candidates.
   • Étape 2 : Pour les fréquences candidates, le motif de modulation journalière de la PSD combinée est ajusté au modèle théorique dérivé de la simulation du halo d'axions. Une statistique de test $Z$ compare la qualité de l'ajustement d'un modèle d'axion par rapport à un modèle de bruit constant.
   • Signification : Un candidat n'est considéré comme significatif que si l'amplitude ajustée est non nulle et si la valeur $p$ de la statistique de test $Z$ tombe en dessous d'un seuil de signification globale ($p_{crit}$) dérivé de la distribution empirique des valeurs de $Z$.

Principales Contributions

• Première Recherche Dédiée : Ce travail présente la première recherche dédiée à un halo d'axions solaire lié par la gravitation en utilisant un réseau mondial de magnétomètres.
• Modélisation du Signal : Les auteurs ont développé un modèle de signal complet pour les champs pseudo-magnétiques générés par les couplages de gradient proton-axion linéaires et quadratiques dans un halo solaire à état fondamental, en tenant explicitement compte de la modulation journalière et de l'interaction entre les composantes de champ radial et transversale à travers le spectre de fréquences (0,05–20 Hz).
• Pipeline d'Analyse : L'étude met en œuvre un pipeline de corrélation croisée robuste qui utilise le décalage temporel pour aligner les signaux globaux à travers un réseau mondial, supprimant efficacement le bruit local non corrélé tout en préservant la cohérence d'un signal lié au soleil.

Résultats

• Absence de Détection : L'analyse de 69 jours de données provenant de 12 stations GNOME n'a révélé aucun signal candidat statistiquement significatif cohérent avec un halo d'axions solaires. Aucune valeur $p$ observée pour la statistique de test $Z$ n'est tombée en dessous du seuil critique de signification globale.
• Limites Supérieures : En l'absence de détection, l'étude établit des limites supérieures de confiance à 95 % sur l'amplitude du champ pseudo-magnétique induit par les axions sur la plage de fréquences $\approx 0,05 - 20$ Hz.
• Contraintes de Couplage : Ces limites de champ ont été traduites en contraintes sur les constantes de couplage proton-axion linéaires ($1/f_l$) et quadratiques ($1/f_q$). Les contraintes ont été dérivées pour deux scénarios de densité de halo de référence :
  1. La surdensité maximale autorisée par les éphémérides planétaires.
  2. La densité prédite par le modèle de capture gravitationnelle par auto-interaction quartique.
• Sensibilité du Couplage Quadratique : Pour le cas du couplage quadratique, les limites dérivées surpassent les limites astrophysiques existantes (par exemple, de la Supernova 1987A) de plus de deux ordres de grandeur sur une grande partie de l'espace des paramètres accessible. Pour les couplages linéaires, les limites sont compétitives ou complémentaires aux limites existantes, selon la densité de halo supposée.

Signification et Revendications
L'article affirme avoir établi les premières contraintes expérimentales sur l'existence d'un halo d'axions solaires formé via la capture gravitationnelle et la thermalisation. Les résultats démontrent que GNOME est capable de sonder des couplages de spin exotiques avec une approche de réseau mondial, en ciblant spécifiquement la signature de modulation journalière unique d'un champ lié au soleil.

Les auteurs notent avec modestie que leurs contraintes sont conditionnelles à l'hypothèse que les axions sont dans l'état fondamental et que le mécanisme de formation (auto-interaction quartique) est efficace. Ils déclarent que pour les couplages quadratiques, leurs résultats représentent une amélioration significative par rapport aux limites astrophysiques précédentes. Le travail souligne également que les données pourraient être réanalysées pour d'autres cibles de physique exotique, telles que les défauts topologiques ou les champs stochastiques, et que les futures itérations de GNOME utilisant des co-magnétomètres alcalins-nobles offriront une sensibilité accrue tant pour les spins des protons que pour ceux des neutrons.