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

Resumen Técnico: Violación de CPT Inducida por Materia y Efectos de Estratificación de la Densidad Terrestre en Experimentos de Oscilación de Neutrinos de Línea de Base Larga

Planteamiento del Problema
Los experimentos de oscilación de neutrinos de línea de base larga (LBL) (T2K, NOνA, DUNE, Hyper-Kamiokande) pretenden medir con precisión la fase de violación de CP leptónica $\delta_{CP}$ y determinar la jerarquía de masas de los neutrinos. Sin embargo, estas mediciones se ven complicadas por dos sistemáticas físicamente acopladas que surgen de la propagación de los neutrinos a través de la materia terrestre:

1. Violación de CPT Extrínseca: La presencia de electrones (pero no de positrones) en la materia terrestre crea un potencial de materia ($V_f$) que cambia de signo entre neutrinos y antineutrinos. Esto induce una asimetría no nula ($A^{CPT}_{\mu e} \neq 0$) incluso si la simetría CPT se cumple exactamente en el vacío. Esta asimetría "falsa" debe ser restada para aislar la violación de CP intrínseca.
2. Estratificación de la Densidad Terrestre: Para líneas de base $L \gtrsim 5000$ km, las trayectorias de los neutrinos atraviesan el manto inferior y el núcleo externo. La aproximación estándar de una densidad promedio constante falla al capturar estas variaciones, lo que puede introducir sesgos sistemáticos en la reconstrucción de $\delta_{CP}$.

El problema central abordado es que ambos efectos dependen del mismo hamiltoniano de potencial de materia ($H_f \propto \rho(x)$). En consecuencia, los errores en el modelado de la densidad se propagan simultáneamente tanto en la estimación de $A^{CPT}_{\mu e}$ como en la reconstrucción de $\delta_{CP}$. Estudios previos han tratado estas sistemáticas de forma aislada; este artículo argumenta que deben modelarse de forma conjunta.

Metodología
Los autores emplean un marco teórico y estadístico unificado:

• Marco Teórico: El análisis utiliza propagadores de exponencial de matriz de tres sabores exactos para resolver la ecuación de evolución $i d|\nu(x)\rangle/dx = H_f(x)|\nu(x)\rangle$. El hamiltoniano incluye el término de vacío (usando parámetros NuFit 5.3) y el potencial de materia MSW derivado del Modelo de Referencia Preliminar de la Tierra (PREM).
• Simulación: Las probabilidades de oscilación de neutrinos se computan para T2K, NOνA, DUNE y Hyper-Kamiokande. El estudio evalúa la asimetría CPT extrínseca definida como $A^{CPT}_{\mu e} = (P_{\mu e} - P_{\bar{e}\bar{\mu}}) / (P_{\mu e} + P_{\bar{e}\bar{\mu}})$.
• Incertidumbre de Densidad: La sensibilidad de $A^{CPT}_{\mu e}$ a las incertidumbres geofísicas se cuantifica aplicando desplazamientos fraccionales de densidad $\rho \to \rho(1 + \epsilon)$ donde $\epsilon = \pm 1\%, \pm 2\%, \pm 5\%$.
• Análisis Estadístico: Se construye un estadístico $\chi^2$ de log-verosimilitud de Poisson para analizar el sesgo en la reconstrucción de $\delta_{CP}$ al utilizar un modelo de densidad constante frente a un "valor real" generado con el perfil completo de PREM. El análisis incluye términos de pull de parámetros de molestia para flujo, sección eficaz y eficiencia del detector. Crucialmente, los autores extienden esto a un "$\chi^2$ extendido" que incluye un término de pull de densidad ($\xi_\rho$) para tratar la incertidumbre de la densidad como una sistemática correlacionada que afecta tanto a $\delta_{CP}$ como a $A^{CPT}_{\mu e}$.

Contribuciones Clave
Este trabajo presenta uno de los primeros estudios sistemáticos que trata la violación de CPT extrínseca inducida por la materia y la estratificación de la densidad terrestre como sistemáticas correlacionadas dentro de un único marco estadístico. Las contribuciones específicas son:

1. Tratamiento Conjunto: Un marco $\chi^2$ unificado que contabiliza simultáneamente las incertidumbres de densidad que afectan tanto a la reconstrucción de la fase de CP como a las estimaciones de la asimetría CPT.
2. Dependencia de la Jerarquía: Una caracterización cuantitativa de cómo los ordenamientos de masa normal (NH) e invertido (IH) afectan a $A^{CPT}_{\mu e}$, revelando una diferencia de $\sim 9\%$ en la línea de base de DUNE.
3. Análisis de Superficie 3D: El cálculo de la superficie $A^{CPT}_{\mu e}(E, \delta_{CP})$, identificando regímenes distintos donde los efectos de CP intrínseco y CPT extrínseco están entrelazados (baja energía) frente a donde la materia domina (alta energía).
4. Cuantificación del Sesgo: Una evaluación rigurosa del sesgo de reconstrucción introducido por la aproximación de densidad constante a través de una amplia gama de líneas de base, desde experimentos actuales hasta propuestas hipotéticas de línea de base muy larga (VLB).

Resultos Clave

• Asimetría CPT Extrínseca: En las energías de pico, $A^{CPT}_{\mu e}$ oscila entre 0.022 (T2K/Hyper-K) y 0.180 (DUNE) bajo jerarquía normal. La asimetría es sensible a $\theta_{23}$, $\theta_{13}$ y $\Delta m^2_{31}$, con una diferencia de $\sim 9\%$ entre NH e IH en DUNE.
• Sensibilidad de la Densidad: Una incertidumbre del $\pm 5\%$ en la densidad de la materia induce un desplazamiento $|\Delta A^{CPT}| \lesssim 0.008$ (aproximadamente 3–4% relativo) en DUNE. La definición de la razón de $A^{CPT}_{\mu e}$ proporciona una autocancelación parcial de los errores de densidad, pero el efecto sigue siendo no despreciable para pruebas de precisión.
• Regímenes de Superficie 3D: En DUNE ($L=1285$ km), la superficie de $A^{CPT}_{\mu e}$ es plana respecto a $\delta_{CP}$ para $E \gtrsim 1.5$ GeV, lo que permite una resta directa. Sin embargo, para $E \lesssim 0.5$ GeV, la superficie se inclina significativamente con $\delta_{CP}$ debido al término de oscilación solar, requiriendo un análisis de verosimilitud conjunto CP–CPT para evitar el sesgo.
• Sesgo de Reconstrucción:
  • Para experimentos actuales ($L \le 1285$ km), el sesgo en la reconstrucción de $\delta_{CP}$ al usar un modelo de densidad constante es insignificante ($< 0.3^\circ$), muy por debajo de los objetivos de precisión experimental actuales.
  • Para las instalaciones propuestas de VLB ($L > 5000$ km), el sesgo se vuelve catastrófico. En $L=7000$ km, el sesgo es de $17.8^\circ$ (3.6 veces el objetivo de la Fase II de DUNE); en $L=12000$ km, el sesgo alcanza los $172.2^\circ$, revirtiendo efectivamente el signo de la fase reconstruida.
• Marco Estadístico: El $\chi^2$ extendido con un término de pull de densidad confirma que, para las líneas de base actuales, la aproximación de densidad constante es adecuada. Sin embargo, para las propuestas de VLB, tratar la incertidumbre de la densidad como una sistemática correlacionada es obligatorio para evitar subestimar las sistemáticas totales.

Significancia y Reivindicaciones
El artículo sostiene que para la generación actual de experimentos LBL (T2K, NOνA, DUNE, Hyper-K), la aproximación de densidad constante es suficiente tanto para los análisis de CP como de CPT, y que la incertidumbre de la densidad de la materia contribuye con solo un manejable $\sim 1.3\%$ sistemático a $A^{CPT}_{\mu e}$.

No obstante, los autores enfatizan que para las futuras propuestas de líneas de base muy largas (por ejemplo, líneas de base transcontinentales $L > 5000$ km), la estratificación de la densidad terrestre y la violación de CPT inducida por la materia no son refinamientos independientes, sino "dos caras de la misma sistemática dominante". El artículo concluye que cualquier análisis futuro de tales instalaciones debe utilizar el riguroso marco estadístico correlacionado presentado aquí, empleando perfiles de PREM espacialmente resueltos y un tratamiento conjunto de las incertididades de densidad para asegurar mediciones válidas de CP y CPT.