High-Energy Neutrino Tomography of the Earth's Interior with IceCube

Resumen Técnico: Tomografía de Neutrinos de Alta Energía del Interior de la Tierra con IceCube

Problema y Motivación
La estructura interna de la Tierra, específicamente su distribución radial de densidad, es fundamental para comprender la evolución geológica, el geodinamo y la formación planetaria. Históricamente, estas propiedades se han inferido principalmente mediante la sismología (analizando la propagación de ondas sísmicas) y la gravimetría, que sondean la respuesta macroscópica de la materia a las fuerzas elásticas y gravitatorias. Si bien el Modelo de Referencia Preliminar de la Tierra (PREM) proporciona un estándar robusto basado en estos métodos, este depende de inferencias indirectas. Un enfoque físico distinto consiste en sondear la densidad de la Tierra utilizando la interacción débil. Los neutrinos de alta energía, que interactúan únicamente a través de la fuerza débil, atraviesan la Tierra pero experimentan una atenuación dependiente de la densidad de columna de la materia que encuentran. A medida que aumenta la energía del neutrino, la sección eficaz de interacción aumenta, lo que hace que la Tierra sea parcialmente opaca para los neutrinos por encima de $\sim$10 TeV. Esta atenuación es sensible al perfil de densidad a lo largo de la trayectoria del neutrino, ofreciendo una sonda novedosa e independiente de los interiores planetarios.

Metodología
Este estudio utiliza 10.7 años de datos (2011–2022) del Observatorio de Neutrinos IceCube en el Polo Sur. El análisis se centra en 368,071 eventos de neutrinos de muones ascendentes en el rango de energía de 500 GeV a 100 TeV. Estos eventos incluyen tanto neutrinos atmosféicos (producidos por interacciones de rayos cósmicos en la atmósfera) como un flujo difuso subdominante astrofísico.

La metodología procede de la siguiente manera:

1. Selección de Eventos: El conjunto de datos comprende trazas de muones "de inicio" (starting) y "transversales" (throughgoing). Las trazas de inicio (interacciones dentro del volumen instrumentado) proporcionan una reconstrucción de energía superior, mientras que las trazas transversales (interacciones fuera del detector) ofrecen una alta estadística.
2. Modelado: El perfil de densidad de la Tierra se parametriza como un conjunto de capas concéntricas de densidad uniforme. El estudio emplea principalmente una configuración de cinco capas (núcleo interno, núcleo externo y capas de manto/corteza) y valida los resultados con una extensión de ocho capas.
3. Simulación y Verosimilitud: Las simulaciones de Monte Carlo modelan la aceptancia del detector IceCube, las incertidumbres sistemáticas (incluyendo las propiedades del hielo glacial) y el flujo de neutrinos esperado. El análisis tiene en cuenta la sección eficaz de interacción neutrino-nucleón del Modelo Estándar, la cual está bien constreñida por las funciones de distribución de partones (PDF) y los datos de dispersión inelástica profunda.
4. Inferencia: Se realiza un ajuste conjunto utilizando una verosimilitud de Poisson por intervalos para constreñir las densidades de las capas, la normalización del flujo, las incertididades de la sección eficaz y los parámetros de nulidad del detector. Las distribuciones posteriores se muestrean utilizando un algoritmo de Monte Carlo de Cadenas de Markov (MCMC).
5. Análisis Ciego: Para evitar sesgos, el análisis se llevó a cabo en dos etapas: un ajuste ciego que excluye eventos casi horizontales ($-1.0 < \cos\theta < -0.1$), donde las profundidades de columna son más pequeñas y el modelado de errores sistemáticos es más probable, seguido de un ajuste utilizando todos los datos disponibles ($-1.0 < \cos\theta < 0.0$).

Resultados Clave
El análisis arroja los siguientes resultados primarios:

• Perfil de Densidad: El perfil de densidad radial inferido, derivado del modelo de cinco capas, es consistente con el modelo PREM dentro de los intervalos de credibilidad del 68% de la densidad más alta (HPD). Los resultados no muestran desviaciones significativas del perfil PREM, validando la sonda de la interacción débil.
  • Capa 1 (Núcleo Interno): Modo de 0.60–0.73 veces la densidad PREM.
  • Capas 2 y 3 (Núcleo Externo): Los modos oscilan entre 0.72 y 1.1 veces la densidad PREM.
  • Capas 4 y 5 (Manto/Corteza): Los modos oscilan entre 1.3 y 1.6 veces la densidad PREM.
  • La inclusión de eventos casi horizontales desplaza ligeramente el modo posterior de la capa más externa (Capa 5) más cerca de la unidad, lo que sugiere que estos eventos contribuyen constructivamente a las restricciones.
• Cantidades Globales Derivadas:
  • Masa de la Tierra: El modo posterior para la masa de la Tierra es $M = 7.25 \times 10^{24}$ kg, con un intervalo HPD del 68% de $[6.31, 8.31] \times 10^{24}$ kg. Aunque el modo es superior al valor de referencia gravitacional ($5.97 \times 10^{24}$ kg), el valor de referencia se encuentra dentro del intervalo HPD del 95%. La hipótesis de masa cero se excluye a $>5\sigma$.
  • Momento de Inercia Polar: El modo posterior es $I = 1.05 \times 10^{38}$ kg m$^2$ (68% HPD: $[9.04, 12.4] \times 10^{37}$ kg m$^2$). El valor de referencia ($8.01 \times 10^{37}$ kg m$^2$) está contenido dentro del intervalo HPD del 95%. El momento de inercia normalizado $I/(MR_\oplus^2)$ se calcula en $\approx 0.357$, comparable al valor de referencia de PREM de 0.331.

Significancia y Reivindicaciones
El artículo afirma presentar la determinación más precisa de la estructura interna de la Tierra basada en la interacción débil hasta la fecha. Al utilizar un orden de magnitud más de datos que intentos previos e implementar un tratamiento de incertidumbre sistemática significativamente actualizado, el estudio demuestra que la atenuación de neutrinos puede servir como una herramienta viable y complementaria a la sismología y la gravimetría.

Los resultados confirman que el perfil de densidad de la Tierra, medido a través de la atenuación de neutrinos, es consistente con los modelos establecidos derivados de la gravedad y las ondas sísmicas. Los autores enfatizan que este trabajo establece a los neutrinos como una herramienta física distinta para la tomografía planetaria. Si bien las incertidumbres actuales están dominadas por el modelado del flujo, el conocimiento de la sección eficaz y las sistemáticas del detector/hielo, el artículo postula que las mejoras futuras en la calibración de detectores (por ejemplo, IceCube Upgrade) y los telescopios de próxima generación más grandes transformarán esta técnica de una verificación complementaria a una sonda de precisión de la estructura planetaria.