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

Technische Zusammenfassung: Materieinduzierte CPT-Verletzung und Effekte der Erdkerndichte-Stratifizierung in Langbasis-Neutrino-Oszillationsexperimenten

Problemstellung
Langbasis-Neutrino-Oszillationsexperimente (T2K, NO$\nu$A, DUNE, Hyper-Kamiokande) zielen darauf ab, die leptonische CP-verletzende Phase $\delta_{CP}$ präzise zu messen und die Neutrinomassehierarchie zu bestimmen. Diese Messungen werden jedoch durch zwei physikalisch gekoppelte Systematiken erschwert, die aus der Ausbreitung von Neutrinos durch die Materie der Erde resultieren:

1. Extrinsische CPT-Verletzung: Das Vorhandensein von Elektronen (aber nicht Positronen) in der Erdmaterie erzeugt ein Materiepotenzial ($V_f$), das das Vorzeichen zwischen Neutrinos und Antineutrinos wechselt. Dies induziert eine nicht verschwindende Asymmetrie ($A^{CPT}_{\mu e} \neq 0$), selbst wenn die CPT-Symmetrie im Vakuum exakt gilt. Diese „falsche“ Asymmetrie muss subtrahiert werden, um die intrinsische CP-Verletzung zu isolieren.
2. Erdkerndichte-Stratifizierung: Für Basislängen $L \gtrsim 5000$ km durchqueren Neutrino-Trajektorien den unteren Erdmantel und den äußeren Erdkern. Die Standardapproximation einer konstanten, pfadgemittelten Dichte versagt bei der Erfassung dieser Variationen, was potenziell systematische Verzerrungen bei der Rekonstruktion von $\delta_{CP}$ einführt.

Das Kernproblem, das hier adressiert wird, ist, dass beide Effekte vom selben Materiepotenzial-Hamiltonian ($H_f \propto \rho(x)$) abhängen. Fehler in der Dichtemodellierung pflanzen sich infolgedessen gleichzeitig sowohl in die Schätzung von $A^{CPT}_{\mu e}$ als auch in die $\delta_{CP}$-Rekonstruktion fort. Zuvorige Studien haben diese Systematiken isoliert behandelt; dieses Paper argumentiert, dass sie gemeinsam modelliert werden müssen.

Methodik
Die Autoren verwenden einen vereinheitlichten theoretischen und statistischen Rahmen:

• Theoretischer Rahmen: Die Analyse verwendet exakte Drei-Flavour-Matrix-Exponentiation-Propagatoren, um die Evolutionsgleichung $i d|\nu(x)\rangle/dx = H_f(x)|\nu(x)\rangle$ zu lösen. Der Hamiltonian enthält den Vakuumterm (unter Verwendung von NuFit 5.3 Parametern) und das MSW-Materiepotenzial, das aus dem Preliminary Reference Earth Model (PREM) abgeleitet ist.
• Simulation: Neutrino-Oszillationswahrscheinlichkeiten werden für T2K, NO$\nu$A, DUNE und Hyper-Kamiokande berechnet. Die Studie evaluiert die extrinsische CPT-Asymmetrie, definiert als $A^{CPT}_{\mu e} = (P_{\mu e} - P_{\bar{e}\bar{\mu}}) / (P_{\mu e} + P_{\bar{e}\bar{\mu}})$.
• Dichtevarianz: Die Sensitivität von $A^{CPT}_{\mu e}$ gegenüber geophysikalischen Unsicherheiten wird quantifiziert, indem fraktionale Dichteverschiebungen $\rho \to \rho(1 + \epsilon)$ mit $\epsilon = \pm 1\%, \pm 2\%, \pm 5\%$ angewendet werden.
• Statistische Analyse: Eine Poisson-Log-Likelihood-$\chi^2$-Statistik wird konstruiert, um die Verzerrung (Bias) der $\delta_{CP}$-Rekonstruktion zu analysieren, wenn ein Modell mit konstanter Dichte gegenüber einem mit dem vollen PREM-Profil generierten „Wahrheitsmodell“ verwendet wird. Die Analyse umfasst Nuisance-Parameter-Pull-Terme für Fluss, Wirkungsquerschnitt und Detektoreffizienz. Entscheidend ist, dass die Autoren dies auf ein „erweitertes $\chi^2$“ ausweiten, das einen Dichte-Pull-Term ($\xi_\rho$) enthält, um die Dichteunsicherheit als korrelierte Systematik zu behandeln, die sowohl $\delta_{CP}$ als auch $A^{CPT}_{\mu e}$ beeinflusst.

Zentrale Beiträge
Diese Arbeit präsentiert eine der ersten systematischen Studien, die Materie-induzierte extrinsische CPT-Verletzung und Erdkerndichte-Stratifizierung innerhalb eines einzigen statistischen Rahmens als korrelierte Systematiken behandelt. Spezifische Beiträge sind:

1. Gemeinsame Behandlung: Ein vereinheitlichter $\chi^2$-Rahmen, der gleichzeitig Dichteunsicherheiten berücksichtigt, welche sowohl die CP-Phasen-Rekonstruktion als auch die CPT-Asymmetrie-Schätzungen beeinflussen.
2. Hierarchie-Abhängigkeit: Eine quantitative Charakterisierung dessen, wie die normale (NH) und die invertierte (IH) Massenordnung $A^{CPT}_{\mu e}$ beeinflussen, wobei eine Differenz von $\sim 9\%$ bei der DUNE-Basislänge festgestellt wurde.
3. 3D-Oberflächenanalyse: Die Berechnung der $A^{CPT}_{\mu e}(E, \delta_{CP})$-Oberfläche, die distinkte Regime identifiziert, in denen intrinsische CP- und extrinsische CPT-Effekte verschränkt sind (niedrige Energie) versus jene, in denen die Materie dominiert (hohe Energie).
4. Bias-Quantifizierung: Eine rigorose Bewertung der Rekonstruktionsverzerrung, die durch die Konstantdichte-Approximation über ein breites Spektrum von Basislängen hinweg eingeführt wird, von aktuellen Experimenten bis hin zu hypothetischen Very-Long-Baseline (VLB)-Vorschlägen.

Wichtigste Ergebnisse

• Extrinsische CPT-Asymmetrie: Bei Spitzenenergien liegt $A^{CPT}_{\mu e}$ zwischen 0,022 (T2K/Hyper-K) und 0,180 (DUNE) unter normaler Hierarchie. Die Asymmetrie ist sensitiv gegenüber $\theta_{23}$, $\theta_{13}$ und $\Delta m^2_{31}$, mit einer Differenz von $\sim 9\%$ zwischen NH und IH bei DUNE.
• Dichte-Sensitivität: Eine $\pm 5\%$ Unsicherheit in der Materiedichte induziert eine Verschiebung $|\Delta A^{CPT}| \lesssim 0,008$ (ca. 3–4 % relativ) bei DUNE. Die Verhältnis-Definition von $A^{CPT}_{\mu e}$ ermöglicht eine teilweise Selbst-Kompensation von Dichtefehlern, dennoch bleibt der Effekt nicht vernachlässigbar für Präzisionstests.
• 3D-Oberflächen-Regime: Bei DUNE ($L=1285$ km) ist die $A^{CPT}_{\mu e}$-Oberfläche bezüglich $\delta_{CP}$ für $E \gtrsim 1,5$ GeV flach, was eine einfache Subtraktion ermöglicht. Jedoch neigt sich die Oberfläche für $E \lesssim 0,5$ GeV aufgrund des solaren Oszillationsterms signifikant mit $\delta_{CP}$, was eine gemeinsame CP–CPT-Likelihood-Analyse erfordert, um Bias zu vermeiden.
• Rekonstruktionsverzerrung (Bias):
  • Für aktuelle Experimente ($L \le 1285$ km) ist die Verzerrung in der $\delta_{CP}$-Rekonstruktion durch die Verwendung eines Konstantdichte-Modells vernachlässigbar ($< 0,3^\circ$), weit unter den aktuellen experimentellen Präzisionszielen.
  • Für vorgeschlagene VLB-Anlagen ($L > 5000$ km) wird die Verzerrung katastrophal. Bei $L=7000$ km beträgt der Bias $17,8^\circ$ (das 3,6-fache des DUNE Phase II Ziels); bei $L=12000$ km erreicht der Bias $172,2^\circ$, was effektiv das Vorzeichen der rekonstruierten Phase umkehrt.
• Statistischer Rahmen: Das erweiterte $\chi^2$ mit einem Dichte-Pull-Term bestätigt, dass für aktuelle Basislängen die Konstantdichte-Approximation adäquat ist. Für VLB-Vorschläge ist jedoch die Behandlung der Dichteunsicherheit als korrelierte Systematik zwingend erforderlich, um die Gesamtsystematik nicht zu unterschätzen.

Bedeutung und Behauptungen
Das Paper behauptet, dass für die aktuelle Generation von LBL-Experimenten (T2K, NO$\nu$A, DUNE, Hyper-K) die Konstantdichte-Approximation sowohl für CP- als auch für CPT-Analysen ausreichend ist und die Materiedichte-Unsicherheit nur zu einer handhabbaren $\sim 1,3\%$ Systematik zu $A^{CPT}_{\mu e}$ beiträgt.

Die Autoren betonen jedoch, dass für zukünftige Very-Long-Baseline-Vorschläge (z. B. transkontinentale Basislängen $L > 5000$ km) die Erdkerndichte-Stratifizierung und die materieinduzierte CPT-Verletzung keine unabhängigen Verfeinerungen sind, sondern „zwei Seiten derselben dominanten Systematik“. Das Paper kommt zu dem Schluss, dass jede zukünftige Analyse solcher Anlagen den hier vorgestellten rigorosen, korrelierten statistischen Rahmen nutzen muss, wobei räumlich aufgelöste PREM-Profile und eine gemeinsame Behandlung der Dichteunsicherheiten einzusetzen sind, um valide CP- und CPT-Messungen zu gewährleisten.