High-Energy Neutrino Tomography of the Earth's Interior with IceCube

Resumo Técnico: Tomografia de Neutrinos de Alta Energia do Interior da Terra com o IceCube

Problema e Motivação
A estrutura interna da Terra, especificamente sua distribuição de densidade radial, é fundamental para compreender a evolução geológica, o geodinamo e a formação planetária. Historicamente, essas propriedades foram inferidas principalmente através da sismologia (analisando a propagação de ondas sísmicas) e da gravimetria, que sondam a resposta macroscópica da matéria a forças elásticas e gravitacionais. Embora o Modelo de Referência Preliminar da Terra (PREM) forneça um padrão robusto baseado nesses métodos, ele depende de inferências indiretas. Uma abordagem física distinta envolve sondar a densidade da Terra usando a interação fraca. Neutrinos de alta energia, que interagem apenas via força fraca, atravessam a Terra, mas experimentam atenuação dependente da densidade de coluna de matéria que encontram. À medida que a energia do neutrino aumenta, sua seção de choque de interação cresce, tornando a Terra parcialmente opaca a neutrinos acima de $\sim$10 TeV. Essa atenuação é sensível ao perfil de densidade ao longo da trajetória do neutrino, oferecendo uma sonda novel e independente dos interiores planetários.

Metodologia
Este estudo utiliza 10,7 anos de dados (2011–2022) do Observatório de Neutrinos IceCube, no Polo Sul. A análise foca em 368.071 eventos de neutrinos de múon ascendentes na faixa de energia de 500 GeV a 100 TeV. Esses eventos incluem tanto neutrinos atmosféricos (produzidos por interações de raios cósmicos na atmosfera) quanto um fluxo astrofísico difuso subdominante.

A metodologia procede da seguinte forma:

1. Seleção de Eventos: O conjunto de dados compreende traços de múons de "início" (starting) e "atravessantes" (throughgoing). Traços de início (interações dentro do volume instrumentado) proporcionam uma reconstrução de energia superior, enquanto traços atravessantes (interações fora do detector) oferecem alta estatística.
2. Modelagem: O perfil de densidade da Terra é parametrizado como um conjunto de camadas concêntricas de densidade uniforme. O estudo utiliza primordialmente uma configuração de cinco camadas (núcleo interno, núcleo externo e camadas de manto/crosta) e valida os resultados com uma extensão de oito camadas.
3. Simulação e Verossimilhança: Simulações de Monte Carlo modelam a aceitação do detector IceCube, incertezas sistemáticas (incluindo propriedades do gelo glacial) e o fluxo de neutrinos esperado. A análise leva em conta a seção de choque neutrino-núcleon do Modelo Padrão, que é bem constrained por funções de distribuição de partons (PDFs) e dados de espalhamento profundo inelástico.
4. Inferência: Uma adequação conjunta (joint fit) é realizada usando uma verossimilhança de Poisson em bins para restringir as densidades das camadas, a normalização do fluxo, incertezas da seção de choque e parâmetros de incerteza (nuisance parameters) do detector. As distribuições posteriores são amostradas usando um algoritmo de Cadeia de Markov Monte Carlo (MCMC).
5. Análise Cega: Para evitar viés, a análise foi conduzida em duas etapas: um ajuste cego excluindo eventos próximos ao horizonte ($-1,0 < \cos\theta < -0,1$), onde as profundidades de coluna são menores e a modelagem errônea sistemática é mais provável, seguido por um ajuste utilizando todos os dados disponíveis ($-1,0 < \cos\theta < 0,0$).

Principais Resultados
A análise produz os seguintes resultados primários:

• Perfil de Densidade: O perfil de densidade radial inferido, derivado do modelo de cinco camadas, é consistente com o modelo PREM dentro dos intervalos de credibilidade de 68% de maior densidade posterior (HPD). Os resultados não mostram desvios significativos do perfil PREM, validando a sonda de interação fraca.
  • Camada 1 (Núcleo Interno): Moda entre 0,60–0,73 vezes a densidade PREM.
  • Camadas 2 & 3 (Núcleo Externo): As modas variam de 0,72 a 1,1 vezes a densidade PREM.
  • Camadas 4 & 5 (Manto/Crosta): As modas variam de 1,3 a 1,6 vezes a densidade PREM.
  • A inclusão de eventos próximos ao horizonte desloca ligeiramente a moda posterior da camada mais externa (Camada 5) para mais próximo da unidade, sugerindo que esses eventos contribuem construtivamente para as restrições.
• Grandezas Globais Derivadas:
  • Massa da Terra: A moda posterior para a massa da Terra é $M = 7,25 \times 10^{24}$ kg, com um intervalo HPD de 68% de $[6,31, 8,31] \times 10^{24}$ kg. Embora a moda seja superior ao valor de referência gravitacional ($5,97 \times 10^{24}$ kg), o valor de referência está contido no intervalo HPD de 95%. A hipótese de massa zero é excluída a $>5\sigma$.
  • Momento de Inércia Polar: A moda posterior é $I = 1,05 \times 10^{38}$ kg m$^2$ (68% HPD: $[9,04, 12,4] \times 10^{37}$ kg m$^2$). O valor de referência ($8,01 \times 10^{37}$ kg m$^2$) está contido no intervalo HPD de 95%. O momento de inércia normalizado $I/(MR_\oplus^2)$ é calculado como $\approx 0,357$, comparável ao PREM de referência de 0,331.

Significância e Alegações
O artigo afirma apresentar a determinação mais precisa do interior da Terra baseada na interação fraca até o momento. Ao utilizar um fator de dez vezes mais dados do que tentativas anteriores e implementar um tratamento de incerteza sistemática significativamente atualizado, o estudo demonstra que a atenuação de neutrinos pode servir como uma sonda viável e complementar à sismologia e gravimetria.

Os resultados confirmam que o perfil de densidade da Terra, conforme medido via atenuação de neutrinos, é consistente com modelos estabelecidos derivados de gravidade e ondas sísmicas. Os autores enfatizam que este trabalho estabelece os neutrinos como uma ferramenta física distinta para tomografia planetária. Embora as incertezas atuais sejam dominadas pela modelagem de fluxo, conhecimento da seção de choque e sistemáticas de detector/gelo, o artigo postula que melhorias futuras na calibração de detectores (ex: IceCube Upgrade) e telescópios de próxima geração maiores transformarão esta técnica de uma verificação complementar em uma sonda de precisão da estrutura planetária.