Energy estimation of cosmic rays by air shower radio signals

Resumen Técnico: Estimación de la Energía de los Rayos Cósmicos mediante Señales de Radio de Lluvia Atmosférica

Planteamiento del Problema
Determinar la energía primaria de los rayos cósmicos es esencial para comprender su origen, aceleración y mecanismos de propagación. Dado que las Lluvias Atmosféricas Extensas (EAS, por sus siglas en inglés) iniciadas por rayos cósmicos destruyen la información directa de la partícula primaria, las partículas secundarias y las emisiones de radio ofrecen vías de reconstrucción. Experimentos previos (por ejemplo, AERA, Tunka-Rex, LOPES) han utilizado la Función de Distribución Lateral (LDF) de las señales de radio, correlacionando la intensidad de la señal a distancias específicas del núcleo de la lluvia con la energía primaria. Sin embargo, estos métodos a menudo requieren identificar una distancia óptima específica para la geometría de cada estación y pueden ser sensibles a la ubicación del núcleo de la lluvia. El experimento SURA, un arreglo de radio de cuatro antenas en Semnan, Irán, ha reconstruido con éxito las direcciones de llegada y las ubicaciones de los núcleos, pero requiere un método robusto para estimar la energía primaria que minimice la sensibilidad a las variaciones en la posición del núcleo.

Metodología
Los autores proponen un método de reconstrucción basado en la comparación de las señales de radio experimentales con simulaciones de alta fidelidad utilizando la extensión CoREAS de CORSIKA. El núcleo de la metodología consiste en los siguientes pasos:

1. Configuración de la Simulación:

   • Configuración de SURA: Las simulaciones se realizaron para coincidir con los parámetros específicos del experimento SURA, incluyendo sus cuatro Antenas de Dipolo Log-Periódico (LPDA), componentes del campo magnético local ($B_x = 39.53 \, \mu T$, $B_z = 28.1 \, \mu T$) y altitud (1130 m).
   • Arreglo de Referencia Denso: Se simuló un arreglo denso virtual de 12,342 antenas con un espaciamiento de 1 metro para generar una referencia de alta resolución de las intensidades de la señal de radio. Este arreglo sirve como base para la comparación.
   • Generación de Eventos: Se simularon lluvias inducidas por protones con energías primarias que oscilan entre $10^{17}$ eV y $7 \times 10^{17}$ eV. La dirección de llegada primaria se fijó en un ángulo cenital de $35^\circ$ y un azimut de $40^\circ$ (basado en la tasa máxima de detección de SURA), con simulaciones adicionales para ángulos cenitales de $60^\circ$ y $65^\circ$ para probar los efectos de la inclinación.

2. Calibración de Datos:

   • El estudio incorpora un proceso de calibración de dos etapas: filtrado de FPGA en tiempo real (filtros de paso alto, paso de banda y muesca) y post-procesamiento (reconstrucción de la dirección de llegada, aplicación de la respuesta de frecuencia, reducción de ruido por wavelet y sobremuestreo de la señal a una resolución de 0.2 ns).

3. El Enfoque del Factor de Escala ($C_{ij}$):

   • Se define un factor de escala, $C_{ij}$, como la razón entre la intensidad media de pico de las señales simuladas del arreglo denso ($\bar{A}_{sim}$) y la intensidad de pico de las señales de la configuración de SURA ($\bar{A}_{exp}$):
     $$C_{ij} = \frac{\bar{A}_{sim}}{\bar{A}_{exp}}$$
   • Este factor se deriva minimizando una función $\chi^2$ que compara las señales experimentales (SURA simulado) y las señales del arreglo denso a través del arreglo de antenas.
   • El método postula que, si bien la intensidad absoluta de la señal depende de la ubicación del núcleo, la razón (factor de escala) entre el arreglo disperso de SURA y el arreglo de referencia denso exhibe una correlación específica con la energía primaria.

Contribuciones Clave

• Independencia de la Ubicación del Núcleo: El artículo demuestra que el factor de escala $C_{ij}$ es ampliamente independiente de la ubicación del núcleo de la lluvia de rayos cósmicos. Incluso para posiciones de núcleo alejadas del centro del arreglo (hasta $\pm 30$ m y $\pm 125$ m en diferentes pruebas), el factor de escala permanece efectivamente constante para una energía dada.
• Correlación de Energía: Los autores establecen una clara correlación inversa entre el factor de escala y la energía primaria. A medida que la energía primaria aumenta, el factor de escala disminuye. Esta relación se modela mediante una función inversa de primer orden:
  $$C_{ij} = y_0 + \frac{a}{E}$$
  donde $E$ es la energía primaria normalizada por $10^{17}$ eV, y $y_0$ y $a$ son constantes de ajuste que dependen de la dirección de llegada.
• Robustez ante la Dirección de Llegada: El método fue probado en lluvias inclinadas (ángulos cenitales de $60^\circ$ y $65^\circ$). La correlación inversa entre $C_{ij}$ y la energía se mantiene a través de estos ángulos, aunque los parámetros de ajuste ($a$ y $y_0$) varían ligeramente con el ángulo cenital.

Resultados

• Precisión de Reconstrucción: Utilizando el factor de escala derivado y la función inversa ajustada, se reconstruyó la energía primaria de los rayos cósmicos simulados. El método logró un error de reconstrucción máximo de aproximadamente 11% (observado en $5 \times 10^{17}$ eV para el ángulo cenital de $35^\circ$).
• Análisis de Error: El error relativo ($\Delta E / E$) se calculó siendo tan bajo como 2% para ciertos niveles de energía (por ejemplo, $2 \times 10^{17}$ eV) y generalmente se mantuvo dentro del rango del 1% al 11% a través del espectro de energía probado.
• Lluvias Inclinadas: Para lluvias más inclinadas ($60^\circ$ y $65^\circ$), el método siguió siendo efectivo, aunque la precisión varió. La menor precisión se observó en $1 \times 10^{17}$ eV y $5 \times 10^{17}$ eV para el ángulo de $65^\circ$, mientras que el error máximo para el ángulo de $60^\circ$ ocurrió en $7 \times 10^{17}$ eV.

Significancia y Reivindicaciones
El artículo afirma que esta técnica ofrece una ventaja distintiva sobre enfoques anteriores al reducir la sensibilidad a la ubicación del núcleo de la lluvia, una fuente común de incertidumbre en la reconstrucción de energía basada en radio. Al utilizar un arreglo de referencia denso para normalizar los datos dispersos experimentales, el método proporciona un estimador confiable para la energía primaria sin requerir la determinación precisa de la distancia óptima desde el núcleo para cada evento específico.

Los autores concluyen que el método es válido tanto para lluvias verticales como inclinadas y puede aplicarse a los datos del experimento SURA una vez que la recolección de datos reales esté completa. Señalan que el trabajo futuro consistirá en aplicar esta técnica a los datos reales del experimento SURA y expandir las simulaciones para incluir lluvias con ubicaciones de núcleo aún más distantes y diversas direcciones de llegada para validar aún más la robustez del método. El estudio no pretende reemplazar los métodos existentes, sino ofrecer un enfoque complementario diseñado específicamente para la geometría y las restricciones del arreglo SURA.