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

Technische Samenvatting: Materie-geïnduceerde CPT-schending en Effecten van de Stratificatie van de Aardse Dichtheid in Long-Baseline Neutrino-oscillatie-experimenten

Probleemstelling
Long-baseline (LBL) neutrino-oscillatie-experimenten (T2K, NOνA, DUNE, Hyper-Kamiokake) beogen het nauwkeurig meten van de leptonische CP-schendende fase $\delta_{CP}$ en het bepalen van de neutrino-massahierarchie. Deze metingen worden echter bemoeilijkt door twee fysiek gekoppelde systematische fouten die voortvloeien uit de propagatie van neutrino's door materie van de Aarde:

1. Extrinsieke CPT-schending: De aanwezigheid van elektronen (maar niet positronen) in de materie van de Aarde creëert een materiepotentiaal ($V_f$) die van teken verandert tussen neutrino's en antineutrino's. Dit induceert een niet-nul asymmetrie ($A^{CPT}_{\mu e} \neq 0$) zelfs als de CPT-symmetrie in vacuüm exact standhoudt. Deze "valse" asymmetrie moet worden afgetrokken om intrinsieke CP-schending te isoleren.
2. Stratificatie van de Aardse Dichtheid: Voor baselines $L \gtrsim 5000$ km doorkruisen neutrino-trajecten de onderste mantel en de buitenste kern. De standaardbenadering van een constante padgemiddelde dichtheid faalt om deze variaties te vatten, wat potentieel systematische biases introduceert bij de reconstructie van $\delta_{CP}$.

Het kernprobleem dat wordt geadresseerd is dat beide effecten afhankelijk zijn van hetzelfde materie-potentiaal-Hamiltoniaan ($H_f \propto \rho(x)$). Hierdoor propageren fouten in de dichtheidsmodellering simultaan in zowel de schatting van $A^{CPT}_{\mu e}$ als de $\delta_{CP}$-reconstructie. Eerdere studies hebben deze systematische effecten geïsoleerd behandeld; dit artikel betoogt dat ze gezamenlijk gemodelleerd moeten worden.

Methodologie
De auteurs hanteren een verenigd theoretisch en statistisch kader:

• Theoretisch Kader: De analyse maakt gebruik van exacte drie-flavor matrix-exponentiële propagatoren om de evolutie-vergelijking $i d|\nu(x)\rangle/dx = H_f(x)|\nu(x)\rangle$ op te lossen. De Hamiltoniaan bevat de vacuümterm (met gebruik van NuFit 5.3 parameters) en de MSW-materiepotentiaal afgeleid van het Preliminary Reference Earth Model (PREM).
• Simulatie: Neutrino-oscillatiekansen worden berekend voor T2K, NOνA, DUNE en Hyper-Kamiokake. De studie evalueert de extrinsieke CPT-asymmetrie gedefinieerd als $A^{CPT}_{\mu e} = (P_{\mu e} - P_{\bar{e}\bar{\mu}}) / (P_{\mu e} + P_{\bar{e}\bar{\mu}})$.
• Dichtheidsonzekerheid: De gevoeligheid van $A^{CPT}_{\mu e}$ voor geofysische onzekerheden wordt gekwantificeerd door fractionele dichtheidsverschuivingen toe te passen waarbij $\rho \to \rho(1 + \epsilon)$ met $\epsilon = \pm 1\%, \pm 2\%, \pm 5\%$.
• Statistische Analyse: Een Poisson log-likelihood $\chi^2$-statistiek wordt geconstrueerd om de bias in $\delta_{CP}$-reconstructie te analyseren wanneer een constant-densiteitsmodel wordt gebruikt tegenover een "waarheid" gegenereerd met het volledige PREM-profiel. De analyse omvat nuisance-parameter pull-termen voor flux, dwarsdoorsnede en detector-efficiëntie. Cruciaal is dat de auteurs dit uitbreiden naar een "uitgebreide $\chi^2$" die een dichtheid-pull-term ($\xi_\rho$) bevat om de dichtheidsonzekerheid te behandelen als een gecorreleerde systematische fout die zowel $\delta_{CP}$ als $A^{CPT}_{\mu e}$ beïnvloedt.

Belangrijkste Bijdragen
Dit werk presenteert een van de eerste systematische studies die materie-geïnduceerde extrinsieke CPT-schending en de stratificatie van de aardse dichtheid als gecorreleerde systematische effecten binnen een enkel statistisch kader behandelt. Specifieke bijdragen zijn:

1. Gezamenlijke Behandeling: Een verenigd $\chi^2$-kader dat gelijktijdig rekening houdt met dichtheidsonzekerheden die zowel de CP-fase reconstructie als de CPT-asymmetrie schattingen beïnvloeden.
2. Hiërarchie-afhankelijkheid: Een kwantitatieve karakterisering van hoe de normale (NH) en geïnverteerde (IH) massaprioriteiten $A^{CPT}_{\mu e}$ beïnvloeden, waarbij een verschil van $\sim 9\%$ wordt vastgesteld bij de DUNE-baseline.
3. 3D-Oppervlakteanalyse: De berekening van het $A^{CPT}_{\mu e}(E, \delta_{CP})$ oppervlak, waarbij onderscheid wordt gemaakt tussen regimes waar intrinsieke CP en extrinsieke CPT-effecten verstrengeld zijn (lage energie) versus waar materie domineert (hoge energie).
4. Bias-kwantificering: Een rigoureuze beoordeling van de reconstructie-bias geïntroduceerd door de constant-densiteitsbenadering over een breed scala aan baselines, van huidige experimenten tot hypothetische zeer-lange-baseline (VLB) voorstellen.

Belangrijkste Resultaten

• Extrinsieke CPT-asymmetrie: Bij piekenergieën varieert $A^{CPT}_{\mu e}$ van 0,022 (T2K/Hyper-K) tot 0,180 (DUNE) onder een normale hiërarchie. De asymmetrie is gevoelig voor $\theta_{23}$, $\theta_{13}$ en $\Delta m^2_{31}$, met een $\sim 9\%$ verschil tussen NH en IH bij DUNE.
• Dichtheidsgevoeligheid: Een onzekerheid van $\pm 5\%$ in de materiedichtheid induceert een verschuiving $|\Delta A^{CPT}| \lesssim 0,008$ (ongeveer 3–4% relatief) bij DUNE. De ratio-definitie van $A^{CPT}_{\mu e}$ zorgt voor een gedeeltelijke zelf-annulatie van dichtheidsfouten, maar het effect blijft niet-verwaarloosbaar voor precisietesten.
• 3D-Oppervlakte-regimes: Bij DUNE ($L=1285$ km) is het $A^{CPT}_{\mu e}$-oppervlak vlak met betrekking tot $\delta_{CP}$ voor $E \gtrsim 1,5$ GeV, wat eenvoudige subtractie mogelijk maakt. Echter, voor $E \lesssim 0,5$ GeV helt het oppervlak significant met $\delta_{CP}$ door de zonnige oscillatieterm, wat een gezamenlijke CP–CPT likelihood-analyse vereist om bias te voorkomen.
• Reconstructie-bias:
  • Voor huidige experimenten ($L \le 1285$ km) is de bias in de $\delta_{CP}$-reconstructie door het gebruik van een constant-densiteitsmodel verwaarloosbaar ($< 0,3^\circ$), ruim onder de huidige experimentele precisiedoelen.
  • Voor voorgestelde VLB-faciliteiten ($L > 5000$ km) wordt de bias catastrofaal. Bij $L=7000$ km is de bias $17,8^\circ$ (3,6 keer de DUNE Fase II doelstelling); bij $L=12000$ km bereikt de bias $172,2^\circ$, wat effectief het teken van de gereconstrueerde fase omdraait.
• Statistisch Kader: De uitgebreide $\chi^2$ met een dichtheid-pull-term bevestigt dat voor huidige baselines de constant-densiteitsbenadering adequaat is. Echter, voor VLB-voorstellen is het behandelen van de dichtheidsonzekerheid als een gecorreleerde systematische fout verplicht om de totale systematische fouten niet te onderschatten.

Betekenis en Claims
Het artikel claimt dat voor de huidige generatie LBL-experimenten (T2K, NOνA, DUNE, Hyper-K), de constant-densiteitsbenadering voldoende is voor zowel CP- als CPT-analyses, en dat de materiedichtheidsonzekerheid slechts een beheersbare $\sim 1,3\%$ systematische fout bijdraagt aan $A^{CPT}_{\mu e}$.

De auteurs benadrukken echter dat voor toekomstige zeer-lange-baseline-voorstellen (bijv. transcontinentale baselines $L > 5000$ km), de aardse dichtheidsstratificatie en de materie-geïnduceerde CPT-schending geen onafhankelijke verfijningen zijn, maar "twee zijden van dezelfde dominante systematische fout". Het artikel concludeert dat elke toekomstige analyse van dergelijke faciliteiten het hier gepresenteerde rigoureuze, gecorreleerde statistische kader moet gebruiken, waarbij ruimtelijk opgeloste PREM-profielen en een gezamenlijke behandeling van dichtheidsonzekerheden noodzakelijk zijn om geldige CP- en CPT-metingen te garanderen.