Revisiting the Radio Lateral Distribution Function: An amplitude dependence on $X_{\rm max}$ and primary composition

Título: Revisitando a Função de Distribuição Lateral (LDF) de Rádio: Uma dependência de amplitude em $X_{max}$ e composição primária

Autores: Washington R. Carvalho Jr. e Lech Wiktor Piotrowski
Instituição: Faculdade de Física, Universidade de Varsóvia, Polônia.

1. O Problema

A detecção de raios cósmicos de ultra-alta energia (UHECR) via radiofrequência é uma técnica madura, mas a reconstrução precisa da composição de massa primária (se o raio cósmico é um próton leve ou um núcleo de ferro pesado) permanece um desafio.

• Contexto Atual: A maioria dos métodos de reconstrução baseia-se no parâmetro $X_{max}$ (a profundidade atmosférica onde a chuva de partículas atinge seu máximo). Tradicionalmente, a dependência da composição na detecção de rádio é atribuída principalmente à forma (inclinação) da Função de Distribuição Lateral (LDF) do sinal de rádio, causada por efeitos geométricos de projeção do cone de Cherenkov no solo.
• A Lacuna: O artigo identifica que existe uma forte dependência da amplitude do campo elétrico de rádio no solo em relação ao $X_{max}$, mesmo após corrigir as diferenças na energia eletromagnética (EM) das chuvas. Historicamente, essa dependência de amplitude foi negligenciada ou atribuída erroneamente apenas a variações na energia EM, especialmente com o advento de métodos de "caixa preta" (como o LOFAR) que focam na comparação de formas de LDF via $\chi^2$, ignorando a amplitude absoluta.

2. Metodologia

Os autores utilizaram uma abordagem combinando simulações físicas detalhadas e análise semi-quantitativa:

• Simulações: Utilizaram o pacote ZHAIRES (combinação de ZHS e AIRES) para simular chuvas de ar extensas (EAS).
  • Cenários: Simulações realizadas para dois locais (Observatório Pierre Auger e o futuro observatório GRAND) e várias energias primárias (1.25 EeV para GRAND, 5 EeV para Auger).
  • Composições: Comparação entre chuvas induzidas por prótons (leves, $X_{max}$ mais profundo) e ferro (pesado, $X_{max}$ mais superficial).
  • Normalização Crucial: Para isolar o efeito do $X_{max}$ da energia, todas as amplitudes de campo elétrico foram normalizadas pela energia eletromagnética (EM) de cada chuva simulada. Isso garante que qualquer diferença restante na amplitude seja puramente devido à geometria da chuva e à densidade atmosférica, não à energia total.
• Análise Teórica: Desenvolveram um modelo físico baseado em duas escalas concorrentes que governam a amplitude do campo elétrico medido no solo:
  1. Escala de distância ($1/R$).
  2. Escala de densidade ($1/\rho$) com correção de coerência ($J$).

3. Contribuições Principais e Mecanismo Físico

O trabalho propõe uma explicação física unificada para a dependência da amplitude de $X_{max}$, baseada na competição entre dois efeitos:

1. Escala de Distância ($1/R$):

   • $R$ é a distância do máximo da chuva ($X_{max}$) até o detector no solo.
   • Chuvas mais profundas (maior $X_{max}$, típicas de prótons) têm um $R$ menor.
   • Efeito: Quanto menor a distância, maior a amplitude do campo elétrico. Isso favorece prótons em ângulos zenitais baixos.

2. Escala de Densidade e Perda de Coerência ($1/\rho \cdot J(\theta)$):

   • $\rho$ é a densidade do ar em $X_{max}$.
   • A emissão geomagnética (principal mecanismo de rádio) é proporcional à corrente transversal, que escala com $1/\rho$. Chuvas mais profundas ocorrem em ar mais denso, o que reduz a velocidade de deriva das partículas e, teoricamente, reduziria a emissão.
   • Fator de Perda de Coerência ($J(\theta)$): Em densidades muito baixas (ângulos zenitais altos) e campos magnéticos fortes, as deflexões das partículas pela força de Lorentz aumentam, introduzindo atrasos de tempo entre as trajetórias das partículas. Isso reduz a coerência da emissão.
   • Efeito: O fator $J(\theta)$ (não linear) modera a escala $1/\rho$. Em campos magnéticos fortes (como no GRAND), a perda de coerência pode ser tão severa que inverte a tendência, fazendo com que a amplitude diminua mesmo com a diminuição da densidade.

Conclusão do Modelo:

• Ângulos Zenitais Baixos: A escala $1/R$ domina. Chuvas mais profundas (prótons) produzem amplitudes maiores.
• Ângulos Zenitais Altos: A escala de densidade (modificada por $J$) torna-se dominante. Dependendo da força do campo magnético local, isso pode favorecer chuvas mais superficiais (ferro) ou inverter a tendência de prótons.

4. Resultados Chave

• Dependência de Amplitude Confirmada: Mesmo com normalização por energia EM, existe uma correlação forte e bem comportada entre a amplitude máxima do LDF e o $X_{max}$. Prótons tendem a ter amplitudes maiores em ângulos baixos, enquanto em ângulos altos a diferença diminui ou inverte dependendo do local.
• Transição de Ângulo Zenital:
  • No Auger (campo magnético mais fraco), há uma transição clara em torno de $72^\circ$. Abaixo disso, prótons têm amplitudes maiores; acima, o ferro pode ter amplitudes maiores ou a dependência se torna nula.
  • No GRAND (campo magnético muito forte), a perda de coerência é tão intensa que a escala $1/R$ permanece dominante em todos os ângulos, mantendo prótons com amplitudes maiores, embora a diferença diminua com o aumento do ângulo.
• O "Ângulo Mágico": Em torno de $82^\circ$, a forma da LDF torna-se quase independente de $X_{max}$ (o tamanho da "pegada" no solo não muda), mas a amplitude ainda carrega forte informação sobre $X_{max}$. Isso é crítico, pois métodos que ignoram a amplitude falham nessa região.
• Validação com Aprendizado de Máquina: O estudo foi motivado por um classificador de Random Forest que, usando apenas distâncias e amplitudes, alcançou 82% de precisão na discriminação de massa. A análise de importância de características confirmou que o algoritmo estava explorando exatamente essa dependência de amplitude em $X_{max}$.

5. Significado e Impacto

• Reconstrução de Composição Direta: A dependência de amplitude permite inferir a composição primária evento a evento, possivelmente contornando a necessidade de reconstruir $X_{max}$ explicitamente antes da classificação.
• Melhoria em Métodos de $\chi^2$: Métodos atuais (como LOFAR/AERA) que usam fatores de escala de energia livres para ajustar dados a simulações podem estar "apagando" a informação de amplitude que contém dados sobre $X_{max}$. O artigo sugere que restringir ou corrigir esses fatores de escala poderia melhorar a resolução de $X_{max}$, especialmente em ângulos altos onde a forma da LDF é menos discriminativa.
• Viés na Estimativa de Energia: Se a fluência de rádio (usada para estimar a energia EM) também depende de $X_{max}$ e não for corrigida adequadamente, pode haver um viés sistemático na estimativa de energia dependendo da composição do raio cósmico.
• Seleção de Locais para Experimentos: Para experimentos focados em chuvas muito inclinadas (quase horizontais), a força do campo magnético local pode ser menos crítica do que se pensava, devido à alta perda de coerência que equaliza as amplitudes em campos fortes e fracos.

Em resumo, o trabalho demonstra que a amplitude do sinal de rádio é tão rica em informações sobre a composição e o desenvolvimento da chuva quanto a forma da distribuição lateral, e que ignorar essa dependência física pode levar a perdas de informação e viéses em reconstruções de raios cósmicos.