Espectros elementales de rayos cósmicos: dos fuentes galácticas y transición extragaláctica más allá de la rodilla ✦
Los espectros energéticos de elementos individuales en rayos cósmicos en el rango de 100×Z GeV a 1000×Z PeV muestran dos quiebres de «rodilla» con dependencia de rigidez: a energías de alrededor de 15×Z TeV y 3×Z PeV, así como tres «tobillos» — aplanamientos. El primer tobillo (≈500×Z GeV) también depende de la rigidez, pero el segundo (≈150 TeV para protones, ~1 PeV para helio) y probablemente el tercero (≈100 PeV para protones) no. Gracias a la alta precisión de las mediciones en los experimentos DAMPE, LHAASO y IceTop, se logró por primera vez establecer que todas las características espectrales se describen bien mediante la suma de tres componentes: una galáctica de baja energía (con una forma convexa del espectro, que decae después de la rodilla de TeV), una galáctica de alta energía (incluyendo la rodilla de PeV) y una extragaláctica. No se requiere ninguna componente galáctica adicional en el intervalo de 10 PeV a 1 EeV.
Context
La naturaleza de la rodilla en el espectro de rayos cósmicos galácticos a ~3 PeV sigue siendo uno de los misterios centrales de la astrofísica de partículas. Históricamente se supuso que podría estar relacionada con la energía máxima de los aceleradores, efectos de propagación o procesos fotonucleares. Ya en la década de 1990 aparecieron indicios de que el espectro está formado por varias poblaciones de fuentes. La idea clásica de Fritz Zwicky (1934) de que los remanentes de supernovas son los principales aceleradores sigue siendo dominante, aunque la configuración exacta de las fuentes y su contribución en diferentes rangos de energía no estaban claras hasta la llegada de mediciones de precisión de nueva generación.
Methods
En el trabajo se utilizaron datos de espectros de protones, helio y hierro de los experimentos DAMPE (mediciones directas en satélite), LHAASO e IceTop (registros en tierra de cascadas atmosféricas extensas). Gracias a la espectroscopía de alta resolución se pudieron identificar elementos individuales. La función de ajuste consistió en la suma de tres componentes (dos galácticas y una extragaláctica). Para describir el espectro de una fuente individual se propuso una nueva forma fenomenológica «doble log-parábola con corte exponencial» (EDLP), que reproduce el aplanamiento observado antes del quiebre. La población de fuentes se modeló mediante la convolución del espectro de una fuente individual con la función de luminosidad, tras lo cual los parámetros se ajustaron a los datos mediante los métodos MINUIT (MIGRAD) y MCMC (emcee). Adicionalmente, se tuvieron en cuenta las pérdidas por fragmentación para los núcleos pesados utilizando una dependencia modelada del recorrido libre con la rigidez.
Results
Se muestra que los espectros de hidrógeno, helio y hierro en el rango de 40 GeV a 500 PeV se describen satisfactoriamente mediante un modelo de tres componentes. La primera componente galáctica tiene una forma convexa con un índice espectral de baja energía de 2.76 (el mismo para todos los núcleos) y termina después de la rodilla de TeV a una rigidez de ~15×Z TeV. La segunda componente domina hasta la rodilla de PeV (3×Z PeV) y es responsable del quiebre principal. La tercera componente, extragaláctica, aparece a ~100 PeV para los protones y se vuelve dominante por encima de 1 EeV. Se encontró que los llamados «tobillos» — quiebres hacia un aplanamiento a 500×Z GeV y ~150 TeV (protones) o ~1 PeV (helio) — son consecuencia de la intersección de diferentes componentes y no siempre obedecen una dependencia de rigidez debido a las diferencias en la composición química de las poblaciones. El ajuste formal del espectro de helio requiere un corte de alta energía de la primera componente a 700 TeV, aunque este detalle no es obligatorio si se tienen en cuenta las incertidumbres sistemáticas. Es importante destacar que para describir los datos en la región de 10 PeV a 1 EeV no es necesaria la introducción de una tercera población galáctica adicional.
Implications
La interpretación propuesta implica que la transición de rayos cósmicos galácticos a extragalácticos ocurre «sin costuras» justo después de la rodilla de PeV, y no en el tradicionalmente considerado «tobillo» del espectro completo a 4 EeV. Esto apoya el modelo de composición «decepcionante» (disappointing model) para energías ultra altas y subraya que la Galaxia deja de ser un acelerador eficaz de protones ya a energías de ~100 PeV. La forma convexa de la primera componente podría indicar efectos no lineales en los procesos de aceleración (aceleración por choque difusivo) o la contribución de fuentes locales en etapas tardías de la evolución de supernovas.
Future development
Futuras mediciones de precisión de los espectros de núcleos del grupo CNO (carbono, oxígeno) y hierro en la región de la rodilla con las instalaciones LHAASO y la planeada SWGO permitirán verificar de forma independiente las predicciones del modelo de tres componentes y la dependencia de rigidez de los quiebres. El análisis conjunto de anisotropía y espectros en el rango de 0.1–10 PeV, especialmente teniendo en cuenta la fase del dipolo, ayudará a identificar la contribución de fuentes cercanas. La simulación numérica de la aceleración en remanentes de supernova con retroalimentaciones no lineales realistas y composición variable del entorno será un paso importante hacia la fundamentación teórica de la imagen observada.
Impact
Los resultados tendrán un impacto en la física de rayos cósmicos, la astronomía gamma de altas energías (búsqueda de PeVatrones galácticos) y la interpretación de datos de telescopios de neutrinos, ya que la transición abrupta a protones extragalácticos cambia las expectativas del fondo difuso de neutrinos.
Next steps
Verificar la dependencia predicha de la transición «corrediza» con la carga nuclear en nuevos datos de CNO y hierro de LHAASO; refinar las secciones eficaces de fragmentación de núcleos pesados en la región de TeV; modelar los espectros de núcleos secundarios (boro, berilio) en el marco del modelo de tres componentes.
Key open problems
El modelo propuesto vincula directamente el origen de la rodilla con el agotamiento de las capacidades de los aceleradores galácticos y la entrada inmediata en escena de las fuentes extragalácticas, lo que requiere una revisión de las ideas tradicionales sobre el rango «galáctico» hasta 1 EeV. Además, la forma convexa de la primera componente y su posible relación con la aceleración no lineal en los frentes de choque de los remanentes de supernovas sigue siendo objeto de intensos debates, conectando el problema de la rodilla con los detalles de la física de plasmas.
🎯 Un protón con una energía de 1 EeV en un campo magnético galáctico típico de 3 microgauss tiene un radio de Larmor de alrededor de 300 pársecs, casi la mitad del espesor del disco galáctico, por lo que la Galaxia no puede retener tales partículas, y casi con certeza provienen de más allá de sus fronteras.
Key numbers
- Energía de la rodilla de TeV para protones: 15 TeV
- Energía de la rodilla de PeV para protones: 3 PeV
- Energía del primer quiebre (tobillo) para protones: 500 GeV
- Energía del segundo quiebre para protones: 150 TeV
- Índice espectral de baja energía de la primera componente: 2.76