Earth-Density Effects in Long Baseline Neutrino Experiments

Resumen Técnico: Efectos de la Densidad Terrestre en Experimentos de Neutrinos de Larga Distancia

Enunciado del Problema
La determinación precisa de la fase de violación de CP ($\delta_{CP}$) en experimentos de oscilación de neutrinos de larga distancia (LBL) depende en gran medida de la modelización precisa de los efectos de la materia terrestre. A medida que los neutrinos se propagan a través de la Tierra, experimentan el efecto Mikheyev–Smirnov–Wolfenstein (MSW), donde la dispersión coherente hacia adelante sobre electrones crea un potencial efectivo que modifica las probabilidades de oscilación. Si bien los análisis actuales suelen emplear una aproximación de "densidad constante" o "densidad promediada a lo largo del trayecto" para simplificar los cálculos, este artículo investiga si tales simplificaciones siguen siendo válidas a medida que los experimentos entran en una era de precisión. Los autores postulan que una modelización incorrecta de la estructura de densidad interna de la Tierra, específicamente la transición de la corteza y el manto al núcleo, puede introducir sesgos sistemáticos que imiten u oculten la violación de CP intrínseca, particularmente para distancias que superan las escalas experimentales actuales.

Metodología
Los autores utilizan un marco completo de oscilación de neutrinos de tres sabores para simular las probabilidades de aparición $\nu_\mu \to \nu_e$. El estudio compara dos enfoques de modelado distintos a lo largo de distancias que van desde $L = 1000$ km hasta $L = 12000$ km:

1. Verdad Fundamental: Un perfil de densidad resuelto espacialmente basado en el Modelo de Referencia Preliminar de la Tierra (PREM), que tiene en cuenta las capas estratificadas de la Tierra (corteza, manto superior/inferior, núcleo externo/interno).
2. Aproximación: Un modelo de densidad constante que utiliza la densidad promedio ponderada por la longitud del trayecto ($\rho_{avg}$) para cada trayectoria específica.

La simulación emplea una técnica de exponenciación de matrices para resolver la ecuación de evolución tipo Schrödinger en la base de sabores, incorporando el Hamiltoniano completo de tres sabores (términos de vacío derivados de la matriz PMNS y diferencias de masa al cuadrado, más el potencial de materia). El análisis asume un ordenamiento de masas normal y una fase de CP real de $\delta_{CP} = -90^\circ$, un valor elegido para maximizar la asimetría de CP leptónica y alinearse con las indicaciones experimentales actuales de T2K y NO$\nu$A.

Para cuantificar el impacto, los autores calculan el sesgo absoluto en el $\delta_{CP}$ reconstruido ($|\Delta \delta_{CP}|$) minimizando una estadística $\chi^2$ de Poisson entre las tasas de eventos "reales" generadas por el PREM y las tasas predichas por el modelo de densidad constante. También analizan los perfiles de $\Delta \chi^2$ y los espectros de tasas de eventos para examinar cómo la modelización incorrecta de la densidad distorsiona el paisaje estadístico y la dependencia energética de la señal de oscilación.

Resultados Clave
El estudio revela un umbral crítico en la longitud de la distancia donde la aproximación de densidad constante pasa de ser una simplificación negligible a una fuente de error sistemático significativo:

• Distancias Cortas a Medias ($L \leq 5000$ km): La aproximación de densidad constante introduce un sesgo negligible de menos de $0.3^\circ$ en el $\delta_{CP}$ reconstruido. A estas distancias, los neutrinos atraviesan principalmente el manto superior y la corteza, donde las variaciones de densidad son lo suficientemente lentas como para que un valor promedio sea suficiente.
• Distancias Largas ($L > 5000$ km): Más allá de los 5000 km, el sesgo aumenta bruscamente. A $L = 7000$ km, el sesgo alcanza $17.8^\circ$. A $L = 12000$ km, donde la trayectoria penetra el núcleo externo de alta densidad, el sesgo se dispara a $172.2^\circ$.
• Comportamiento No Monótono: El sesgo no es estrictamente monótono; disminuye de $17.8^\circ$ a 7000 km a $8.8^\circ$ a 9000 km antes de volver a aumentar. Los autores atribuyen esto a degeneraciones de fase de oscilación donde el mínimo de $\chi^2$ del modelo de densidad constante coincide accidentalmente con el valor real.
• Divergencia Estructural: A $L = 7000$ km, el modelo de densidad constante predice un pico de oscilación más agudo y pronunciado cerca de 3.0 GeV en comparación con el pico más amplio y suprimido predicho por el perfil PREM. Por el contrario, por encima de 4 GeV, el perfil PREM produce tasas de eventos significativamente mayores debido al potencial acumulado del manto inferior más denso. Esta divergencia supera las incertidumbres estadísticas de $1\sigma$, lo que indica que los modelos son estadísticamente distinguibles.
• Error Fundamental: El análisis demuestra que el sesgo surge de distorsiones dependientes de la energía causadas por la estratificación interna de la Tierra. Estas distorsiones crean degeneraciones entre los efectos inducidos por la materia y la violación de CP intrínseca que no pueden resolverse marginalizando sobre un único parámetro de densidad efectiva.

Significado y Afirmaciones
El artículo concluye que la aproximación de densidad constante no es una "simplificación conservadora", sino una fuente de error sistemático fundamental para los experimentos de larga distancia. Los autores argumentan que, a medida que los experimentos avancen hacia una mayor precisión, el perfil de densidad de la Tierra debe tratarse como una incertidumbre sistemática resuelta espacialmente en lugar de una simple constante numérica.

Si bien los autores señalan que ningún experimento LBL basado en aceleradores propuesto actualmente supera una distancia de aproximadamente 3000 km (donde la aproximación sigue siendo válida), enfatizan la relevancia de sus hallazgos para:

1. Propuestas Futuras de Muy Larga Distancia: Experimentos que planean distancias superiores a 5000 km.
2. Análisis de Neutrinos Atmosféricos: Detectores como IceCube-Upgrade y KM3NeT/ORCA, donde los neutrinos atraviesan una distribución continua de distancias y ángulos cenitales, haciendo que la modelización precisa de la densidad sea esencial para la determinación simultánea de $\delta_{CP}$ y el ordenamiento de masas.
3. Marcos Sistemáticos: La necesidad de incorporar perfiles PREM continuos y propagar las incertidumbres de los modelos terrestres (incluyendo desviaciones regionales y zonas de transición) en los marcos de análisis de experimentos de próxima generación para asegurar que las asimetrías observadas entre neutrinos y antineutrinos reflejen física fundamental y no suposiciones ambientales.