On the Nonlinear Dependence of Underground Muon Rate on Atmospheric Temperature Observed at Daya Bay

Resumo Técnico: Sobre a Dependência Não Linear da Taxa de Múons Subterrâneos em Relação à Temperatura Atmosférica Observada em Daya Bay

Definição do Problema
A modulação do fluxo de múons de raios cósmicos subterrâneos pela temperatura atmosférica é um fenômeno bem estabelecido, tradicionalmente explicado pelas teorias de Barrett, Gaisser e outros. Essas teorias postulam uma correlação positiva entre a taxa de múons e a temperatura atmosférica, impulsionada pela competição entre o decaimento (píons e kaons) e a interação de mésons. No entanto, os dados do Experimento de Neutrinos de Daya Bay, particularmente nas salas experimentais mais rasas EH1 e EH2, exibem um efeito de modulação não linear: a correção parece mais forte durante os meses quentes do que nos meses frios. Os arcabouços teóricos existentes, que frequentemente dependem de suposições de atmosfera isotérmica ou de projeções assintóticas de perturbações de temperatura em pesos de profundidade unidimensionais, falham em explicar essa não linearidade quando aplicados a perfis de temperatura realistas e não isotérmicos.

Metodologia
Os autores propõem um arcabouço teórico revisado e uma estratégia de verificação numérica:

1. Solução Analítica Aperfeiçoada: O artigo deriva uma solução aprimorada para as equações de evolução de cascata atmosférica para os espectros de energia de mésons e múons. Diferente de aproximações anteriores que vinculam os limites de alta e baixa energia ou assumem condições isotérmicas, esta solução retém a dependência não local completa do perfil de temperatura atmosférica. Especificamente, ela trata o comprimento de decaimento como um funcional explícito da relação densidade-temperatura, preservando o integral do histórico de temperatura ao longo da trajetória do méson, do ponto de produção ao decaimento.
2. Novas Definições: Com base nesta solução aperfeiçoada, os autores redefinem a temperatura efetiva ($T_{eff}$), sua variação ($\Delta T_{eff}$) e o coeficiente de temperatura ($\alpha_T$). Em vez de usar quantidades projetadas isotérmicas, estas são definidas via a derivada funcional da taxa de múons em relação ao perfil de temperatura real e variante no tempo $T(Y)$. Esta abordagem decompõe a resposta em um termo local (probabilidade de decaimento na camada de produção) e um termo de propagação (histórico de sobrevivência do méson).
3. Simulação Numérica (MCEq): Para validar a teoria, os autores utilizam a ferramenta Matrix Cascade Equations (MCEq). Eles inserem dados atmosféricos reais do banco de dados ERA5 (abrangendo 2011–2021) e dados de terreno específicos para as salas experimentais de Daya Bay (EH1, EH2, EH3).
   • A simulação considera o fluxo de raios cósmicos primários (modelo Hillas Gaisser 2012), interações hadrônicas (SIBYLL2.3C) e a sobreposição de rocha (rock overburden).
   • O estudo reproduz as variações da taxa de múons ao longo de 11 anos.
   • Dois caminhos de cálculo são comparados: um utilizando o método de ponderação convencional de Grashorn (que projeta perturbações em um peso de profundidade 1D) e outro utilizando o novo método de derivada funcional (calculado numericamente via diferenças finitas no MCEq).

Principais Resultados

• Reprodução da Não Linearidade: Utilizando soluções de cascata convencionais e o método de ponderação padrão de Grashorn, os autores reproduzem com sucesso a correlação não linear entre a mudança relativa na taxa de múons subterrâneos e a mudança relativa na temperatura efetiva observada em Daya Bay. Demonstra-se que o coeficiente de temperatura $\alpha_T$ é maior em temperaturas mais altas.
• Restauração da Linearidade: Quando a variação da taxa de múons é calculada usando os novos pesos baseados na derivada funcional (que consideram o histórico completo de temperatura não local), a relação entre a taxa de múons e a temperatura efetiva torna-se linear.
• Validação: O cálculo numérico da derivada funcional usando MCEq confirma as previsões teóricas derivadas na Seção 3. A dependência linear é recuperada através de diferentes modelos de interação (SIBYLL2.3C, 2.1, 2.3E) e limiares de energia.

Significância e Alegações
O artigo afirma que a não linearidade observada em Daya Bay não é uma quebra fundamental da física da correlação temperatura-múon, mas sim um artefato das aproximações usadas nas definições convencionais de temperatura efetiva. Ao tratar o perfil de temperatura como um input funcional em vez de projetá-lo sobre um peso simplificado, os autores demonstram que a relação física subjacente é linear.

O trabalho fornece:

1. Uma solução mais geral para as equações de cascata que descreve plenamente como todo o perfil de temperatura influencia a taxa final de múons.
2. Definições refinadas para temperatura efetiva e o coeficiente de temperatura que são válidas para atmosferas não isotérmicas.
3. Um arcabouço de cálculo que pode explicar os efeitos não lineares encontrados em Daya Bay e oferece um método mais preciso para calcular coeficientes de temperatura para outros experimentos subterrâneos.

Os autores concluem que seu método recupera a dependência linear esperada, sugerindo que as interpretações anteriores de modulação não linear foram limitadas pelas suposições de projeção isotérmica inerentes aos tratamentos analíticos padrão.