Modeling the Effect of the Heliospheric Magnetic Field on Cosmic Ray Muon Shadows

Resumo Técnico: Modelagem do Efeito do Campo Magnético Heliosférico nas Sombras de Múons de Raios Cósmicos

Enunciado do Problema
Os raios cósmicos (RCs), principalmente prótons e partículas alfa, são defletidos por campos eletromagnéticos ao atravessarem o sistema solar. Embora o campo magnético geomagnético (GMF) da Terra seja bem mapeado, o Campo Magnético Heliosférico (HMF) é complexo e difícil de modelar. O Sol e a Lua absorvem os RCs incidentes, criando "sombras" no fluxo de raios cósmicos observáveis da Terra. No entanto, essas sombras não são estáticas; elas são deslocadas e distorcidas pelo HMF e pelo GMF. A Colaboração MINOS+ utilizou dados de múons do detector distante MINOS para localizar essas sombras solares durante diferentes fases do ciclo solar (mínimo solar, máximo solar e o intervalo completo de 13 anos). O problema central abordado é se os deslocamentos e distorções observados na sombra solar de múons são consistentes com o modelo padrão da espiral de Parker do HMF combinado com um modelo dipolar do GMF e, caso contrário, quais modificações nos modelos ou parâmetros são necessárias para conciliar a teoria com a observação.

Metodologia
O estudo utilizou 13 anos de dados de múons (julho de 2003 a julho de 2016) do detector distante MINOS, localizado a 2070 metros de equivalente em água subterrâneos na Mina Soudan. A análise focou em eventos de múons únicos com ângulos zenitais incidentes menores que 80° e comprimento de trajetória suficiente para reconstruir trajetórias.

1. Localização da Sombra: A equipe empregou um teste de probabilidade $\Lambda$ (uma razão de verossimilhança logarítmica) para localizar a sombra. Isso envolveu comparar as contagens de partículas observadas ($N_{obs}$) próximas ao Sol contra um fundo plano teórico ($N_{back}$) e um modelo de sombra teórico ($N_{th}$). O modelo levou em conta o tamanho angular do Sol (0,52°) e os efeitos de espalhamento na sobrecarga de rocha.
2. Validação Estatística: Simulações de Monte Carlo foram usadas para gerar histogramas de fundo aleatórios para determinar a significância das localizações de sombra observadas. A probabilidade de flutuações aleatórias mimetizarem a sombra foi calculada, resultando em níveis de confiança elevados (até 6$\sigma$ para o conjunto de dados completo).
3. Modelagem Computacional: Para interpretar os deslocamentos da sombra, os autores simularam as trajetórias reversas no tempo dos RCs do detector de volta ao Sol.
   • Modelos de Campo: A simulação superpôs o modelo da espiral de Parker para o HMF e um modelo dipolar para o GMF. Os parâmetros da espiral de Parker incluíram uma velocidade constante do vento solar ($u_{sowi} = 400$ km/s) e frequência de rotação solar. O GMF foi modelado usando um momento dipolar e convertido em um sistema de coordenadas centrado no Sol.
   • Dinâmica de Partículas: A equação da força de Lorentz foi integrada usando um esquema Runge-Kutta de 4ª ordem. As partículas foram tratadas como relativísticas, com energias seguindo uma distribuição de lei de potência ($E^{-2.75}$).
   • Cenários: As simulações testaram várias razões de polaridade do campo magnético (50/50, 33/67, 67/33) e dois tipos de modelo de sombra: uma sombra circular simples e uma sombra "distorcida" derivada dos histogramas reais de probabilidade $\Lambda$.
   • Variação de Parâmetros: O corte de energia dos raios cósmicos primários ($E_0$) foi tratado como um parâmetro livre para determinar se um espectro de energia diferente poderia conciliar o modelo com as posições de sombra observadas.

Principais Contribuições

• Caracterização da Sombra: O artigo localizou com sucesso sombras solares de múons com alta significância estatística em três períodos distintos: o conjunto de dados completo de 13 anos, o mínimo solar (2007–2009) e o máximo solar (2012–2014).
• Desenvolvimento de Modelos: Os autores introduziram um método para criar modelos de amostragem de sombra "distorcidos" baseados nas distribuições de probabilidade $\Lambda$ observadas, indo além da suposição de uma sombra perfeitamente circular.
• Teste de Modelos: O estudo testou rigorosamente o modelo da espiral de Parker contra dados observacionais, quantificando o acordo entre os deslocamentos de sombra simulados e observados em Ascensão Reta ($\Delta RA$) e Declinação ($\Delta DEC$).
• Análise do Espectro de Energia: O artigo investigou a sensibilidade das posições das sombras ao espectro de energia dos RCs primários, testando especificamente o impacto do corte de baixa energia ($E_0$).

Resultados

• Deslocamentos da Sombra: As sombras observadas mostraram deslocamentos distintos dependendo da fase do ciclo solar. A sombra do mínimo solar foi deslocada de forma diferente da sombra do máximo solar, e o conjunto de dados completo representou uma combinação desses deslocamentos.
• Consistência do Modelo:
  • O modelo da espiral de Parker mostrou o melhor acordo com os dados do mínimo solar.
  • O modelo mostrou o menor acordo com os dados do máximo solar, particularmente no que diz respeito ao deslocamento $\Delta DEC$.
  • Modelos de sombra distorcida performaram melhor para os períodos de máximo solar e dados completos, enquanto modelos redondos simples foram mais eficazes para dados de mínimo solar.
• Efeitos de Polaridade: As simulações confirmaram que a razão de polaridade do campo magnético afeta diretamente a posição $\Delta DEC$ da sombra. Os dados pareciam favorecer um HMF dominado por polaridade positiva, especialmente durante o máximo solar.
• Discrepância no Corte de Energia: Para alinhar as médias das sombras simuladas com o centro teórico (0,0) usando o modelo da espiral de Parker, o corte de energia dos RCs primários ($E_0$) teve que ser ajustado para 22,8 $\pm$ 1,2 TeV. Isso é significativamente maior (2–3 vezes) do que o valor experimentalmente determinado de $\approx$10 TeV.
• Correção Zenital: Uma correção para a dependência do ângulo zenital do corte de energia foi aplicada. Embora isso tenha ajustado o espectro para corresponder melhor à distribuição de energia dos múons observada, foi insuficiente para explicar totalmente o deslocamento $\Delta RA$ observado.

Significância e Conclusões
O artigo conclui que as sombras solares de múons de raios cósmicos são observáveis e podem ser usadas para sondar a estrutura do HMF. Os resultados indicam que a sombra do mínimo solar é a mais consistente com o modelo padrão da espiral de Parker combinado com um GMF dipolar. No entanto, o modelo falha em reproduzir totalmente a sombra do máximo solar e o conjunto de dados completo sem invocar um corte de energia unrealisticamente alto.

Os autores afirmam que a discrepância sugere que o modelo da espiral de Parker é uma simplificação excessiva para as condições presentes durante o máximo solar. Eles propõem que fatores não modelados, como a inclinação da folha de corrente neutra (o "Modelo Bailarina"), frentes de choque e variações complexas de polaridade, provavelmente contribuem para os deslocamentos observados. O estudo estabelece que, embora a direção geral da deflexão seja consistente com a espiral de Parker, a estrutura detalhada do HMF durante o máximo solar requer modelagem mais complexa do que a fornecida pelo modelo espiral simples. O trabalho demonstra a utilidade das observações de sombras de múons para restringir modelos de HMF, mas destaca a necessidade de simulações mais sofisticadas para contabilizar as variações do ciclo solar.