Co-evolution of cosmic ray energy spectra, composition, and anisotropies

Resumen Técnico: Co-evolución de los Espectros de Energía, Composición y Anisotropías de los Rayos Cósmicos

Planteamiento del Problema
El origen de los rayos cósmicos (RC) sigue siendo un problema fundamental sin resolver en la astrofísica, principalmente debido a la desviación de las partículas cargadas por los campos magnéticos, que oscurece las ubicaciones de las fuentes, y a la dificultad de determinar la composición de masa mediante mediciones indirectas desde tierra. Aunque los avances recientes en experimentos espaciales y terrestres han proporcionado datos de alta precisión sobre los espectros de energía, la composición de masa (masa logarítmica media, $\langle \ln A \rangle$) y las anisotropías a gran escala, estos observables a menudo se han analizado de forma independiente. Existe una brecha crítica en la comprensión de la evolución correlacionada de estos observables multi-mensajero a través de un vasto rango de energías (desde $\mathcal{O}(10)$ GeV hasta $10^{11}$ GeV). Específicamente, los modelos anteriores han tenido dificultades para explicar simultáneamente las estructuras espectrales finas (características de endurecimiento y ablandamiento), la variación en la composición y las complejas inversiones de fase en la anisotropía sin invocar múltiples poblaciones de fuentes distintas.

Metodología
Los autores proponen un modelo unificado de cuatro componentes para interpretar conjuntamente los espectros de energía, la composición y las anisotropías de los RC. El modelo descompone el flujo total de RC en cuatro poblaciones de fuentes distintas:

1. Fuentes de Fondo Galácticas (Componente A): Una población distribuida (por ejemplo, restos de supernova) modelada con un escenario de propagación dependiente espacialmente (SDP) para tener en cuenta las variaciones de difusión entre el disco y el halo galácticos. El espectro de inyección es una ley de potencia con corte exponencial en rigidez.
2. Fuente Galáctica Cercana (Componente B): Una fuente local única, ilustrada por el púlsar Geminga, inyectada como una fuente puntual con una edad y distancia específicas. Su contribución se calcula utilizando un método de función de Green.
3. Fuentes Extragalácticas de Baja Energía (Componente C): Una población responsable de la transición del dominio galáctico al extragaláctico.
4. Fuentes Extragalácticas de Alta Energía (Componente D): Una segunda población extragaláctica responsable del corte de energía más alta.

El modelo emplea el código DRAGON para resolver la ecuación de transporte para los componentes galácticos, incorporando difusión, convección y re-aceleración. Para los componentes extragalácticos, se utiliza un coeficiente de difusión fenomenológico para ajustar los datos de anisotropía por encima de energías de EeV. El espectro observado total es la suma algebraica de los cuatro componentes, mientras que la anisotropía es la suma vectorial de sus direcciones individuales de flujo. Los parámetros del modelo (índices espectrales, rigideces de corte, normalizaciones) se ajustan para reproducir los grupos de "endurecimiento-ablandamiento" observados en los datos, la evolución de $\langle \ln A \rangle$, y la amplitud y fase del dipolo de las anisotropías.

Contribuciones y Resultados Clave
El estudio identifica cuatro grupos correlacionados de características espectrales de "endurecimiento-ablandamiento" a través del espectro de energía, vinculando cada uno a transiciones de fuentes específicas y límites físicos:

• Grupo 1 ($\sim 10$ GeV – $100$ TeV): Se caracteriza por un endurecimiento en $\sim 200$ GV y un ablandamiento en $\sim 15$ TV. El modelo atribuye esto a la superposición del fondo galáctico (A) y la fuente cercana (B). La fuente cercana domina el flujo y la dirección de la anisotropía entre $\sim 100$ GeV y $100$ TeV, explicando la inversión de fase observada en la anisotropía dipolar (desplazándose desde la dirección del Centro Galáctico hacia la dirección de la fuente cercana).
• Grupo 2 ($\sim 100$ TeV – $10$ PeV): Incluye la "rodilla" en $\sim 3\text{--}4$ PeV. Esta característica surge a medida que disminuye la contribución de la fuente cercana (B), permitiendo que el fondo galáctico (A) domine nuevamente, hasta alcanzar la energía máxima de aceleración de protones por fuentes galácticas.
• Grupo 3 ($\sim 10$ PeV – $1$ EeV): Presenta un endurecimiento en $\sim 30$ PeV y la "segunda rodilla" en $\sim 200$ PeV. El endurecimiento es causado por la aparición del primer componente extragaláctico (C), mientras que la segunda rodilla corresponde a la energía máxima de aceleración de núcleos de hierro por fuentes galácticas (A). La transición en la fase de anisotropía ocurre alrededor de $1$ EeV a medida que el flujo del componente C supera al componente A.
• Grupo 4 ($\sim 1$ EeV – $100$ EeV): Incluye el "tobillo" en $\sim 5$ EeV y la supresión espectral en $\sim 50$ EeV. El tobillo se interpreta como la transición del componente extragaláctico C al componente D. La supresión de energía más alta se atribuye al límite de aceleración del componente D, y no al proceso GZK.

Correlaciones de Composición y Anisotropía
El modelo reproduce con éxito las estructuras observadas de "bulto-hundimiento" en la masa logarítmica media ($\langle \ln A \rangle$) que se correlacionan con el endurecimiento y ablandamiento espectrales. Los autores argumentan que a medida que aparecen nuevos componentes (típicamente dominados por protones), la masa media se aligera, y a medida que se desvanecen o alcanzan cortes, la composición se vuelve más pesada. De manera similar, la amplitud y fase de la anisotropía a gran escala muestran transiciones distintas en las energías características de estos cambios de fuente. Cabe destacar que el modelo sugiere que la transición del origen galáctico al extragaláctico ocurre alrededor de $\mathcal{O}(10^8)$ GeV, lo cual es más bajo que la energía tradicional del tobillo.

Significado y Afirmaciones
El artículo afirma que la co-evolución de los espectros, la composición y las anisotropías proporciona una fuerte evidencia de un origen multi-poblacional de los rayos cósmicos. Al unificar estos observables, el modelo de cuatro componentes ofrece una explicación coherente para las estructuras complejas en el espectro de RC sin depender de una única población de fuentes homogéneas. Los autores enfatizan que, aunque las mediciones actuales de composición y anisotropía a altas energías (por encima de decenas de PeV) y bajas energías (por debajo de TeV) conllevan incertidumbres significativas, las tendencias correlacionadas observadas respaldan el escenario propuesto. El trabajo destaca la necesidad de futuras mediciones de alta precisión tanto de experimentos de detección directa (por ejemplo, HERD, VLAST) como de detección indirecta para restringir aún más las poblaciones de fuentes y validar las energías de transición específicas y los comportamientos de fase de anisotropía predichos por el modelo.