Addressing uncertainties of model predictions for extensive air showers initiated by high energy cosmic rays

Problema de Pesquisa
Estudos experimentais de raios cósmicos de ultra-alta energia (UHECRs) dependem da reconstrução das propriedades das partículas primárias a partir das características de chuvas atmosféricas extensas (EAS), especificamente a profundidade média do máximo da chuva ($X_{max}$) e o número de múons ao nível do solo ($N_\mu$). Essas reconstruções dependem fortemente de geradores de Monte Carlo (MC) que modelam interações hadrônicas na atmosfera. Um desafio significativo reside em quantificar as incertezas desses geradores. Embora a comparação entre diferentes modelos existentes ofereça algum insight, os modelos atuais podem compartilhar deficiências comuns ou omitir física crítica, o que pode levar a incertezas subestimadas ou desvios sistemáticos da realidade. O Observatório Pierre Auger, por exemplo, indica que os modelos atuais podem exigir valores de $X_{max}$ substancialmente maiores para serem consistentes com a composição de massa observada, enquanto um "problema do múon" sugere que os modelos subestimam o número de múons em chuvas induzidas por UHECRs.

Metodologia
Os autores utilizam um novo gerador de MC de colisões hadrônicas, o QGSb, que é baseado em um formalismo teórico transparente (modelo de String de Quark-Gluon), mas possui liberdade de parâmetros suficiente para explorar as incertezas do modelo. O estudo envolve a modificação sistemática de parâmetros específicos do modelo para alterar as previsões de $X_{max}$ e $N_\mu$, aderindo estritamente às restrições impostas por dados de aceleradores (ex: LHC, TOTEM, ALICE, NA61/SHINE, EHS-NA22).

A investigação foca em duas áreas principais:

1. Modificando $X_{max}$: Os autores exploram mecanismos para aumentar ou diminuir a profundidade prevista do máximo da chuva. Isso envolve o ajuste da dependência de energia da seção de choque inelástica próton-ar, a inelasticidade das interações ($K_{inel}$) e as distribuições de momento dos partões constituintes. As modificações específicas incluem:
   • O ajuste da sobrecriticidade do "Pomeron semihardo" ($\Delta_{sh}$).
   • A alteração do expoente da distribuição de momento dos quarks (anti)mar de (anti)quarks ($\alpha^{(sh)}_{sea}$).
   • A implementação de um mecanismo de "quebra de diquark", onde a coerência dos quarks de valência em um núcleon é destruída durante a espalhamento múltiplo, potencialmente aumentando a inelasticidade.
2. Modificando $N_\mu$: Os autores investigam formas de aumentar o conteúdo de múons alterando o balanço de energia entre hádrons secundários estáveis (núcleons, kaons) e píons neutros ($\pi^0$) em interações píon-ar. Como o $N_\mu$ é dominado pela produção frontal de hádrons estáveis, o estudo foca em:
   • Aumentar a troca de píon no canal-$t$ para favorecer píons carregados sobre os neutros.
   • Aumentar a produção frontal de (anti)bárens e kaons na região de fragmentação do projétil, modificando os parâmetros de fragmentação de string (ex: probabilidades para criar pares de quarks estranhos ou pares diquark-antidiquark a partir do vácuo).

Principais Contribuições e Resultados

• Aumentando $X_{max}$: O estudo descobre que a forma mais eficaz de aumentar o $X_{max}$ é assumir distribuições de momento mais suaves para os partões constituintes acoplados às strings de campo de cor. Ao aumentar o expoente $\alpha^{(sh)}_{sea}$ de 0,5 (padrão) para 0,8 e 0,95 (Tunes 3 e 4), a inelasticidade das colisões próton-nitrogênio é reduzida em até 11%, resultando em um aumento de $X_{max}$ de até $\approx 40$ g/cm$^2$ a $10^{20}$ eV.

  • Restrição: Essas modificações extremas são desmentidas por dados recentes do ALICE a $\sqrt{s}=13$ TeV. O experimento ALICE observa uma forte anticorrelação entre a multiplicidade de hádrons carregados centrais e o rendimento de nêutrons frontais. As Tunes 3 e 4 preveem uma anticorrelação muito mais fraca, contradizendo os dados. Além disso, essas modificações levam a uma profundidade máxima de produção de múons ($X^\mu_{max}$) que exigiria núcleos primários mais pesados que o ferro para coincidir com os dados de Pierre Auger, o que é fisicamente improvável.
  • Alternativa: Uma redução moderada na sobrecriticidade do Pomeron ($\Delta_{sh} = 0,21$ vs. 0,22, Tune 2) produz um aumento menor em $X_{max}$ ($\approx 7$ g/cm$^2$) e permanece consistente com as correlações do ALICE, embora subestime ligeiramente a perda de energia em nêutrons líderes.

• Diminuindo $X_{max}$: Os autores testaram o mecanismo de "quebra de diquark" (Tune 5), onde o espalhamento múltiplo destrói a coerência do diquark, teoricamente aumentando a inelasticidade. Apesar de uma previsão de aumento de 3% na inelasticidade em altas energias, este mecanismo não produziu nenhuma mudança perceptível no $X_{max}$ médio (dentro das incertezas estatísticas de 1–2 g/cm$^2$). Os autores concluem que as previsões do modelo padrão para $X_{max}$ provavelmente já estão próximas do limite inferior, pois as interações periféricas (que dominam a média) permanecem inalteradas por este mecanismo.

• Aumentando $N_\mu$: O estudo confirma que o $N_\mu$ é altamente sensível à produção frontal de hádrons estáveis.

  • O aumento da produção frontal de kaons (Tune 7) aumentou o $N_\mu$ em até $\approx 2\%$. No entanto, esta tune superestima os rendimentos de kaons frontais em outros sistemas de colisão (ex: $K^+$ em colisões $\pi p$) e falha em corresponder a certos espectros experimentais.
  • O aumento da produção frontal de (anti)núcleons (Tune 8), modificando a hadronização de quarks "remanescentes" para incluir mais conteúdo bariônico, aumentou o $N_\mu$ em até $\approx 5\%$. Embora esta tune melhore o acordo com alguns dados do NA61/SHINE, ela leva a um excesso de fator de dois na produção de prótons/antiprótons em comparação com os dados do LEBC-EHS a 360 GeV/c.
  • O estudo destaca que a incerteza em $N_\mu$ é atualmente limitada por sérias contradições entre diferentes medições de aceleradores de rendimentos de (anti)prótons na região de fragmentação.

Significância e Alegações
O artigo afirma fornecer um mapeamento detalhado do "espaço de parâmetros" disponível para o ajuste de modelos de interação hadrônica dentro dos limites dos dados atuais de aceleradores. Sua principal significância reside em:

1. Descartar cenários extremos: Demonstra que as modificações mais drásticas ao $X_{max}$ (via distribuições de partões suaves) são incompatíveis com os dados de física frontal do LHC (ALICE) e com os dados de profundidade de múon de Auger.
2. Identificar limites na redução de $X_{max}$: Desafia alegações da literatura de que a divisão de diquark reduz significativamente o $X_{max}$, mostrando que tais mecanismos têm impacto negligível no máximo da chuva médio quando restringidos por espectros de prótons de energia fixa.
3. Quantificar a incerteza do "Problema do Múon": Estabelece que aumentar o $N_\mu$ em 5% é teoricamente possível dentro do framework QGSb, mas requer suposições sobre a produção frontal de (anti)núcleons que estão atualmente em conflito com os dados experimentais existentes.

Os autores concluem que, embora o gerador QGSb permita explorar essas incertezas, as restrições atuais dos dados de aceleradores (particularmente em relação à produção de partículas frontais e correlações) limitam significativamente a capacidade de resolver as discrepâncias entre as observações de EAS e as previsões dos modelos sem introduzir novos mecanismos de física ainda não considerados.