Matter-Induced CPT Violation and Earth-Density Stratification Effects in Long-Baseline Neutrino Oscillation Experiments

Résumé Technique : Violation de la CPT induite par la matière et effets de stratification de la densité terrestre dans les expériences d'oscillation de neutrinos à longue ligne de base

Énoncé du Problème
Les expériences d'oscillation de neutrinos à longue ligne de base (LBL) (T2K, NO$\nu$A, DUNE, Hyper-Kamiokande) visent à mesurer avec précision la phase de violation de CP leptonique $\delta_{CP}$ et à déterminer la hiérarchie de masse des neutrinos. Cependant, ces mesures sont compliquées par deux systématiques physiquement couplées découlant de la propagation des neutrinos à travers la matière terrestre :

1. Violation de la CPT Extrinsèque : La présence d'électrons (mais pas de positrons) dans la matière terrestre crée un potentiel de matière ($V_f$) qui change de signe entre les neutrinos et les antineutrinos. Cela induit une asymétrie non nulle ($A^{CPT}_{\mu e} \neq 0$) même si la symétrie CPT est respectée exactement dans le vide. Cette asymétrie « factice » doit être soustraite pour isoler la violation de CP intrinsèque.
2. Stratification de la Densité Terrestre : Pour des lignes de base $L \gtrsim 5000$ km, les trajectoires des neutrinos traversent le manteau inférieur et le noyau externe. L'approximation standard d'une densité moyenne constante échoue à capturer ces variations, ce qui peut introduire des biais systématiques dans la reconstruction de $\delta_{CP}$.

Le problème central abordé est que les deux effets dépendent du même hamiltonien de potentiel de matière ($H_f \propto \rho(x)$). Par conséquent, les erreurs de modélisation de la densité se propagent simultanément dans l'estimation de $A^{CPT}_{\mu e}$ et dans la reconstruction de $\delta_{CP}$. Les études antérieures ont traité ces systématiques de manière isolée ; cet article soutient qu'elles doivent être modélisées conjointement.

Méthodologie
Les auteurs emploient un cadre théorique et statistique unifié :

• Cadre Théorique : L'analyse utilise des propagateurs exacts par exponentiation de matrice à trois saveurs pour résoudre l'équation d'évolution $i d|\nu(x)\rangle/dx = H_f(x)|\nu(x)\rangle$. L'hamiltonien inclut le terme de vide (utilisant les paramètres NuFit 5.3) et le potentiel de matière MSW dérivé du modèle de référence de la Terre (PREM).
• Simulation : Les probabilités d'oscillation des neutrinos sont calculées pour T2K, NO$\nu$A, DUNE et Hyper-Kamiokande. L'étude évalue l'asymétrie CPT extrinsèque définie par $A^{CPT}_{\mu e} = (P_{\mu e} - P_{\bar{e}\bar{\mu}}) / (P_{\mu e} + P_{\bar{e}\bar{\mu}})$.
• Incertitude de Densité : La sensibilité de $A^{CPT}_{\mu e}$ aux incertitudes géophysiques est quantifiée en appliquant des décalages fractionnaires de densité $\rho \to \rho(1 + \epsilon)$ où $\epsilon = \pm 1\%, \pm 2\%, \pm 5\%$.
• Analyse Statistique : Une statistique $\chi^2$ de log-vraisemblance de Poisson est construite pour analyser le biais dans la reconstruction de $\delta_{CP}$ lors de l'utilisation d'un modèle à densité constante par rapport à une « vérité » générée avec le profil complet PREM. L'analyse inclut des termes de tirage (pull terms) de paramètres de nuisance pour le flux, la section efficace et l'efficacité du détecteur. Crucialement, les auteurs étendent cela à un « $\chi^2$ étendu » qui inclut un terme de tirage de densité ($\xi_\rho$) pour traiter l'incertitude de densité comme une systématique corrélée affectant à la fois $\delta_{CP}$ et $A^{CPT}_{\mu e}$.

Contributions Clés
Ce travail présente l'une des premières études systématiques traitant la violation de la CPT extrinsèque induite par la matière et la stratification de la densité terrestre comme des systématiques corrélées au sein d'un cadre statistique unique. Les contributions spécifiques sont :

1. Traitement Conjoint : Un cadre $\chi^2$ unifié qui prend en compte simultanément les incertitudes de densité affectant à la fois la reconstruction de la phase CP et les estimations d'asymétrie CPT.
2. Dépendance de la Hiérarchie : Une caractérisation quantitative de la manière dont les ordonnancements de masse normal (NH) et inversé (IH) affectent $A^{CPT}_{\mu e}$, révélant une différence d'environ $9\%$ à la ligne de base de DUNE.
3. Analyse de Surface 3D : Le calcul de la surface $A^{CPT}_{\mu e}(E, \delta_{CP})$, identifiant des régimes distincts où les effets de CP intrinsèque et de CPT extrinsèque sont entremêlés (basse énergie) par rapport aux régimes où la matière domine (haute énergie).
4. Quantification du Biais : Une évaluation rigoureuse du biais de reconstruction introduit par l'approximation de densité constante à travers une large gamme de lignes de base, allant des expériences actuelles aux propositions hypothétiques de très longues lignes de base (VLB).

Résultats Clés

• Asymétrie CPT Extrinsèque : Aux énergies de pic, $A^{CPT}_{\mu e}$ varie de 0,022 (T2K/Hyper-K) à 0,180 (DUNE) sous la hiérarchie normale. L'asymétrie est sensible à $\theta_{23}$, $\theta_{13}$ et $\Delta m^2_{31}$, avec une différence d'environ $9\%$ entre NH et IH à DUNE.
• Sensibilité de la Densité : Une incertitude de $\pm 5\%$ sur la densité de la matière induit un décalage $|\Delta A^{CPT}| \lesssim 0,008$ (environ 3–4 % relatif) à DUNE. La définition par rapport au ratio de $A^{CPT}_{\mu e}$ permet une auto-annulation partielle des erreurs de densité, mais l'effet reste non négligeable pour les tests de précision.
• Régimes de Surface 3D : À DUNE ($L=1285$ km), la surface $A^{CPT}_{\mu e}$ est plate par rapport à $\delta_{CP}$ pour $E \gtrsim 1,5$ GeV, permettant une soustraction directe. Cependant, pour $E \lesssim 0,5$ GeV, la surface s'incline significativement avec $\delta_{CP}$ en raison du terme d'oscillation solaire, nécessant une analyse de vraisemblance conjointe CP–CPT pour éviter tout biais.
• Biais de Reconstruction :
  • Pour les expériences actuelles ($L \le 1285$ km), le biais dans la reconstruction de $\delta_{CP}$ provenant de l'utilisation d'un modèle à densité constante est négligeable ($< 0,3^\circ$), bien en dessous des cibles de précision expérimentale actuelle.
  • Pour les installations VLB proposées ($L > 5000$ km), le biais devient catastrophique. À $L=7000$ km, le biais est de $17,8^\circ$ (3,6 fois la cible de la Phase II de DUNE) ; à $L=12000$ km, le biais atteint $172,2^\circ$, inversant efficacement le signe de la phase reconstruite.
• Cadre Statistique : Le $\chi^2$ étendu avec un terme de tirage de densité confirme que pour les lignes de base actuelles, l'approximation de densité constante est adéquate. Cependant, pour les propositions VLB, traiter l'incertitude de densité comme une systématique corrélée est impératif pour éviter de sous-estimer le total des systématiques.

Signification et Revendications
L'article affirme que pour la génération actuelle d'expériences LBL (T2K, NO$\nu$A, DUNE, Hyper-K), l'approximation de densité constante est suffisante pour les analyses CP et CPT, et que l'incertitude de la densité de matière ne contribue qu'à un système de $\sim 1,3\%$ pour $A^{CPT}_{\mu e}$.

Cependant, les auteurs soulignent que pour les futures propositions de très longues lignes de base (ex: baselines transcontinentales $L > 5000$ km), la stratification de la densité terrestre et la violation de la CPT induite par la matière ne sont pas des raffinements indépendants mais « deux faces d'une même systématique dominante ». L'article conclut que toute analyse future de ces installations doit utiliser le cadre statistique rigoureux et corrélé présenté ici, en employant des profils PREM spatialement résolus et un traitement conjoint des incertitudes de densité pour garantir des mesures de CP et de CPT valides.