Energy estimation of cosmic rays by air shower radio signals

Resumo Técnico: Estimativa de Energia de Raios Cósmicos por Sinais de Radio de Chuva Atmosférica

Definição do Problema
Determinar a energia primária dos raios cósmicos é essencial para compreender sua origem, aceleração e mecanismos de propagação. Embora as Extensas Chuvas Atmosféricas (EAS) iniciadas por raios cósmicos destruam a informação direta da partícula primária, as partículas secundárias e as emissões de rádio oferecem camros de reconstrução. Experimentos anteriores (ex: AERA, Tunka-Rex, LOPES) utilizaram a Função de Distribuição Lateral (LDF) de sinais de rádio, correlacionando a intensidade do sinal a distâncias específicas do núcleo da chuva com a energia primária. No entanto, esses métodos frequentemente exigem a identificação de uma distância ideal específica para a geometria de cada estação e podem ser sensíveis à localização do núcleo da chuva. O experimento SURA, um arranjo de rádio de quatro antenas em Semnan, Irã, reconstruiu com sucesso as direções de chegada e as localizações dos núcleos, mas requer um método robusto para estimar a energia primária que minimize a sensibilidade às variações da posição do núcleo.

Metodologia
Os autores propõem um método de reconstrução baseado na comparação de sinais de rádio experimentais com simulações de alta fidelidade utilizando a extensão CoREAS do CORSIKA. O cerne da metodologia envolve as seguintes etapas:

1. Configuração da Simulação:

   • Configuração SURA: As simulações foram realizadas para corresponder aos parâmetros específicos do experimento SURA, incluindo seus quatro Log-Periodic Dipole Arrays (LPDA), componentes do campo magnético local ($B_x = 39,53 \, \mu T$, $B_z = 28,1 \, \mu T$) e altitude (1130 m).
   • Arranjo de Referência Denso: Um arranjo virtual denso de 12.342 antenas com espaçamento de 1 metro foi simulado para gerar uma referência de alta resolução para as intensidades de sinal de rádio. Este arranjo serve como uma linha de base para comparação.
   • Geração de Eventos: Chuvas induzidas por prótons foram simuladas com energias primárias variando de $10^{17}$ eV a $7 \times 10^{17}$ eV. A direção de chegada primária foi definida para um ângulo zenital de $35^\circ$ e azimute de $40^\circ$ (baseado na taxa de detecção de pico do SURA), com simulações adicionais para ângulos zenitais de $60^\circ$ e $65^\circ$ para testar efeitos de inclinação.

2. Calibração de Dados:

   • O estudo incorpora um processo de calibração de dois estágios: filtragem FPGA em tempo real (filtros passa-alta, passa-banda, notch) e pós-processamento (reconstrução da direção de chegada, aplicação de resposta de frequência, redução de ruído por wavelet e up-sampling de sinal para 0,2 ns de resolução).

3. A Abordagem do Fator de Escala ($C_{ij}$):

   • Um fator de escala, $C_{ij}$, é definido como a razão entre a intensidade de pico média dos sinais simulados do arranjo denso ($\bar{A}_{sim}$) e a intensidade de pico dos sinais da configuração SURA ($\bar{A}_{exp}$):
     $$C_{ij} = \frac{\bar{A}_{sim}}{\bar{A}_{exp}}$$
   • Este fator é derivado pela minimização de uma função $\chi^2$ que compara os sinais experimentais (SURA simulado) e os sinais do arranjo denso através do arranjo de antenas.
   • O método postula que, embora a intensidade absoluta do sinal dependa da localização do núcleo, a razão (fator de escala) entre o arranjo esparso do SURA e o arranjo de referência denso exibe uma correlação específica com a energia primária.

Principais Contribuições

• Independência da Localização do Núcleo: O artigo demonstra que o fator de escala $C_{ij}$ é amplamente independente da localização do núcleo da chuva de raios cósmicos. Mesmo para posições de núcleo distantes do centro do arranjo (até $\pm 30$ m e $\pm 125$ m em diferentes testes), o fator de escala permanece efetivamente constante para uma determinada energia.
• Correlação de Energia: Os autores estabelecem uma clara correlação inversa entre o fator de escala e a energia primária. À medida que a energia primária aumenta, o fator de escala diminui. Esta relação é modelada por uma função inversa de primeira ordem:
  $$C_{ij} = y_0 + \frac{a}{E}$$
  onde $E$ é a energia primária normalizada por $10^{17}$ eV, e $y_0$ e $a$ são constantes de ajuste dependentes da direção de chegada.
• Robustez à Direção de Chegada: O método foi testado em chuvas inclinadas (ângulos zenitais de $60^\circ$ e $65^\circ$). A correlação inversa entre $C_{ij}$ e a energia mantém-se nestes ângulos, embora os parâmetros de ajuste ($a$ e $y_0$) variem ligeiramente com o ângulo zenital.

Resultos

• Precisão de Reconstrução: Utilizando o fator de escala derivado e a função inversa ajustada, a energia primária dos raios cósmicos simulados foi reconstruída. O método alcançou um erro de reconstrução máximo de aproximadamente 11% (observado em $5 \times 10^{17}$ eV para o ângulo zenital de $35^\circ$).
• Análise de Erro: O erro relativo ($\Delta E / E$) foi calculado como sendo tão baixo quanto 2% para certos níveis de energia (ex: $2 \times 10^{17}$ eV) e manteve-se geralmente na faixa de 1%–11% em todo o espectro de energia testado.
• Chuvas Inclinadas: Para chuvas mais inclinadas ($60^\circ$ e $65^\circ$), o método permaneceu eficaz, embora a precisão tenha variado. A menor precisão foi observada em $1 \times 10^{17}$ eV e $5 \times 10^{17}$ eV para o ângulo de $65^\circ$, enquanto o erro máximo para o ângulo de $60^\circ$ ocorreu em $7 \times 10^{17}$ eV.

Significância e Alegações
O artigo afirma que esta técnica oferece uma vantagem distinta sobre abordagens anteriores ao reduzir a sensibilidade à localização do núcleo da chuva, uma fonte comum de incerteza na reconstrução de energia baseada em rádio. Ao utilizar um arranjo de referência denso para normalizar os dados esparsos experimentais, o método fornece um estimador confiável para a energia primária sem exigir a determinação precisa da distância ótima do núcleo para cada evento específico.

Os autores concluem que o método é válido tanto para chuvas verticais quanto inclinadas e pode ser aplicado aos dados do experimento SURA assim que a coleta de dados reais for concluída. Eles observam que trabalhos futuros envolverão a aplicação desta técnica aos dados experimentais reais do SURA e a expansão das simulações para incluir chuvas com localizações de núcleo ainda mais distantes e direções de chegada variadas para validar ainda mais a robustez do método. O estudo não pretende substituir os métodos existentes, mas oferece uma abordagem complementar especificamente adaptada à geometria e às restrições do arranjo SURA.