Modeling the Effect of the Heliospheric Magnetic Field on Cosmic Ray Muon Shadows

Resumen Técnico: Modelado del Efecto del Campo Magnético Heliosférico en las Sombras de Muones de Rayos Cósmicos

Enunciado del Problema
Los rayos cósmicos (RC), principalmente protones y partículas alfa, son desviados por campos electromagnéticos al atravesar el sistema solar. Si bien el campo magnético geomagnético (GMF) de la Tierra está bien cartografiado, el Campo Magnético Heliosférico (HMF) es complejo y difícil de modelar. El Sol y la Luna absorben los RC incidentes, creando "sombras" en el flujo de rayos cósmicos observables desde la Tierra. Sin embargo, estas sombras no son estáticas; son desplazadas y distorsionadas por el HMF y el GMF. La Colaboración MINOS+ utilizó datos de muones del detector lejano de MINOS para localizar estas sombras solares durante diferentes fases del ciclo solar (mínimo solar, máximo solar y el lapso completo de 13 años). El problema central abordado es si los desplazamientos y distorsiones observados de la sombra solar de muones son consistentes con el modelo estándar de espiral de Parker del HMF combinado con un modelo dipolar del GMF, y de no ser así, qué modificaciones a los modelos o parámetros son necesarias para conciliar la teoría con la observación.

Metodología
El estudio utilizó 13 años de datos de muones (julio de 2003–julio de 2016) del detector lejano de MINOS, ubicado a 2070 metros de equivalente en agua bajo tierra en la Mina Soudan. El análisis se centró en eventos de muones individuales con ángulos cenitales incidentes menores a 80° y longitud de trayectoria suficiente para reconstruir las trayectorias.

1. Localización de la Sombra: El equipo empleó una prueba de probabilidad $\Lambda$ (una relación de verosimilitud logarítmica) para localizar la sombra. Esto implicó comparar los conteos de partículas observados ($N_{obs}$) cerca del Sol contra un fondo plano teórico ($N_{back}$) y una plantilla de sombra teórica ($N_{th}$). La plantilla tuvo en cuenta el tamaño angular del Sol (0.52°) y los efectos de dispersión en la sobrecarga de roca.
2. Validación Estadística: Se utilizaron simulaciones de Monte Carlo para generar histogramas de fondo aleatorios y determinar la significancia de las ubicaciones de sombra observadas. Se calculó la probabilidad de que fluctuaciones aleatorias imiten la sombra, obteniéndose niveles de confianza altos (hasta 6$\sigma$ para el conjunto de datos completo).
3. Modelado Computacional: Para interpretar los desplazamientos de la sombra, los autores simularon las trayectorias de los RC en reverso desde el detector hasta el Sol.
   • Modelos de Campo: La simulación superpuso el modelo de espiral de Parker para el HMF y un modelo dipolar para el GMF. Los parámetros de la espiral de Parker incluyeron una velocidad constante del viento solar ($u_{sowi} = 400$ km/s) y la frecuencia de rotación solar. El GMF se modeló utilizando un momento dipolar y se convirtió a un sistema de coordenadas centrado en el Sol.
   • Dinámica de Partículas: La ecuación de la fuerza de Lorentz se integró utilizando un esquema Runge-Kutta de cuarto orden. Las partículas se trataron como relativistas, con energías siguiendo una distribución de ley de potencia ($E^{-2.75}$).
   • Escenarios: Las simulaciones probaron diversas relaciones de polaridad del campo magnético (50/50, 33/67, 67/33) y dos tipos de plantillas de sombra: una sombra circular simple y una sombra "distorsionada" derivada de los histogramas reales de probabilidad $\Lambda$.
   • Variación de Parámetros: El corte de energía de los rayos cósmicos primarios ($E_0$) se trató como un parámetro libre para determinar si un espectro de energía diferente podía conciliar el modelo con las posiciones de sombra observadas.

Contribuciones Clave

• Caracterización de la Sombra: El artículo localizó con éxito sombras solares de muones con alta significancia estadística en tres marcos temporales distintos: el conjunto de datos completo de 13 años, el mínimo solar (2007–2009) y el máximo solar (2012–2014).
• Desarrollo de Plantillas: Los autores introdujeron un método para crear plantillas de muestreo de sombra "distorsionadas" basadas en las distribuciones de probabilidad $\Lambda$ observadas, superando la suposición de una sombra perfectamente circular.
• Prueba de Modelos: El estudio probó rigurosamente el modelo de espiral de Parker contra datos observacionales, cuantificando el acuerdo entre los desplazamientos de sombra simulados y observados en Ascensión Recta ($\Delta RA$) y Declinación ($\Delta DEC$).
• Análisis del Espectro de Energía: El artículo investigó la sensibilidad de las posiciones de la sombra al espectro de energía de los RC primarios, probando específicamente el impacto del corte de baja energía ($E_0$).

Resultados

• Desplazamientos de la Sombra: Las sombras observadas mostraron desplazamientos distintos dependiendo de la fase del ciclo solar. La sombra del mínimo solar se desplazó de manera diferente a la sombra del máximo solar, y el conjunto de datos completo representó una combinación de estos desplazamientos.
• Consistencia del Modelo:
  • El modelo de espiral de Parker mostró el mejor acuerdo con los datos del mínimo solar.
  • El modelo mostró el menor acuerdo con los datos del máximo solar, particularmente con respecto al desplazamiento $\Delta DEC$.
  • Las plantillas de sombra distorsionadas funcionaron mejor para los períodos de máximo solar y datos completos, mientras que las plantillas redondas simples fueron más efectivas para los datos del mínimo solar.
• Efectos de Polaridad: Las simulaciones confirmaron que la relación de polaridad del campo magnético afecta directamente la posición $\Delta DEC$ de la sombra. Los datos parecían favorecer un HMF dominado por polaridad positiva, especialmente durante el máximo solar.
• Discrepancia del Corte de Energía: Para alinear las medias de sombra simuladas con el centro teórico (0,0) utilizando el modelo de espiral de Parker, el corte de energía de los RC primarios ($E_0$) tuvo que ajustarse a 22.8 $\pm$ 1.2 TeV. Esto es significativamente más alto (2–3 veces) que el valor experimentalmente determinado de $\approx$10 TeV.
• Corrección Cenital: Se aplicó una corrección para la dependencia del ángulo cenital del corte de energía. Si bien esto ajustó el espectro para coincidir mejor con la distribución de energía de muones observada, fue insuficiente para explicar completamente el desplazamiento $\Delta RA$ observado.

Significado y Conclusiones
El artículo concluye que las sombras solares de muones de rayos cósmicos son observables y pueden utilizarse para sondear la estructura del HMF. Los resultados indican que la sombra del mínimo solar es la más consistente con el modelo estándar de espiral de Parker combinado con un GMF dipolar. Sin embargo, el modelo falla en reproducir completamente la sombra del máximo solar y el conjunto de datos completo sin invocar un corte de energía unrealmente alto.

Los autores afirman que la discrepancia sugiere que el modelo de espiral de Parker es una simplificación excesiva para las condiciones presentes durante el máximo solar. Proponen que factores no modelados, como la inclinación de la lámina de corriente neutra (el "Modelo de la Bailarina"), frentes de choque y variaciones complejas de polaridad, probablemente contribuyen a los desplazamientos observados. El estudio establece que, si bien la dirección general de la desviación es consistente con la espiral de Parker, la estructura detallada del HMF durante el máximo solar requiere un modelado más complejo que el que proporciona el modelo de espiral simple. El trabajo demuestra la utilidad de las observaciones de sombras de muones para restringir modelos de HMF, pero destaca la necesidad de simulaciones más sofisticadas para dar cuenta de las variaciones del ciclo solar.