Reassessing the gallium anomaly using self-consistent electron wave functions

Este artigo reavalia a duradoura anomalia do gálio ao calcular as seções de choque de captura de neutrinos utilizando funções de onda eletrônicas autoconsistentes derivadas da equação de Dirac-Coulomb, atualizando, assim, a significância global da discrepância e sua interpretação em relação a neutrinos estéreis.

O essencial

Imagine que você está tentando contar quantas gotas de chuva caem em um balde específico. Você tem uma fórmula matemática muito precisa que prevê exatamente quantas gotas deveriam cair no balde com base no tamanho da chuva e no tamanho do balde. No entanto, toda vez que você realmente conta as gotas na vida real, descobre que há menos gotas do que sua fórmula previu. Essa chuva perdida é o que os físicos chamam de "Anomalia do Gálio".

Por mais de 30 anos, experimentos usando Gálio-71 (um tipo de metal) como um "balde" para capturar neutrinos (partículas minúsculas e fantasmagóricas vindas do sol ou de fontes radioativas) têm encontrado consistentemente menos neutrinos do que o esperado. O hiato entre a previsão e a realidade cresceu tanto que agora é considerado um grande mistério na física.

Este artigo de Cadeddu e colegas é como uma equipe de mecânicos mestres decidindo reconstruir o motor dessa fórmula de previsão do zero para ver se cometeram um erro de cálculo.

O Jeito Antigo vs. O Jeito Novo

O Motor Antigo (A Aproximação):
Anteriormente, os cientistas calculavam como os neutrinos interagem com os átomos de Gálio usando uma versão de "rascunho" da matemática. Eles tratavam os elétrons (as partículas minúsculas que orbitam o núcleo do átomo) como se fossem ondas simples e suaves que não mudavam muito dentro do átomo. Era como estimar a forma de uma estrada acidentada apenas olhando para um mapa plano. Eles assumiam que a onda do elétron era a mesma em todos os lugares dentro do núcleo minúsculo.

O Novo Motor (A Solução Exata):
Neste novo estudo, os autores decidiram parar de usar o mapa plano. Em vez disso, eles usaram um GPS de alta definição para resolver as equações exatas (chamadas de equação de Dirac-Coulomb) que descrevem como os elétrons realmente se comportam.

• A Analogia: Imagine que o núcleo é uma pista de dança lotada. O método antigo assumia que todos na pista de dança estavam parados em um círculo perfeito. O novo método realmente conta cada dançarino, levando em conta como eles esbarram uns nos outros e se movem no espaço lotado. Eles resolveram a matemática tanto para os elétrons presos ao átomo quanto para os que estão saindo, usando um programa de computador especializado para obter a forma exata da "onda" do elétron.

O Truque da "Média"

Outra mudança fundamental neste artigo é como eles lidam com o tamanho do núcleo.

• O Jeito Antigo: Eles escolhiam um único ponto no centro do núcleo (como medir a temperatura de uma sala colocando um termômetro exatamente no meio) e assumiam que isso representava todo o conjunto.
• O Novo Jeito: Eles perceberam que o núcleo tem um tamanho, então eles "tiraram a média" do comportamento do elétron por todo o volume do núcleo. É como tirar a temperatura em cada ponto da sala e encontrar a verdadeira média, em vez de apenas adivinhar com base no centro.

O Que Eles Descobriram?

Quando rodaram seus novos cálculos, mais precisos:

1. A Previsão Mudou: Sua nova fórmula, mais precisa, previu que menos neutrinos deveriam ser capturados do que a fórmula antiga previa.
2. O Hiato Aumentou: Como a nova previsão é menor, a diferença entre o que se espera ver e o que os experimentos realmente viram tornou-se ainda maior.
3. O Resultado: Os "neutrinos perdidos" agora são um problema de 5,5-sigma. No mundo da ciência, "sigma" é uma medida de confiança. Um resultado de 5-sigma é o padrão ouro para uma descoberta, o que significa que há menos de 1 chance em 3,5 milhões de que essa discrepância seja apenas um acaso aleatório.

A Hipótese do "Neutrino Estéril"

Os físicos têm uma teoria favorita para explicar essa chuva perdida: Neutrinos Estéreis.

• A Metáfora: Imagine que os neutrinos são como um bando de pássaros voando em direção ao balde. A teoria sugere que alguns desses pássaros são "invisíveis" (estéreis). Eles não interagem com o balde de forma alguma; eles apenas voam direto através dele. Se esses pássaros invisíveis existirem, eles explicariam por que o balde está mais vazio do que o esperado.

Os autores atualizaram a matemática para ver se essa teoria dos "pássaros invisíveis" ainda se encaixa.

• A Boa Notícia: A matemática ainda permite que esses pássaros invisíveis existam. Os dados do Gálio ainda apontam fortemente para a presença deles.
• A Má Notícia: Outros experimentos (observando neutrinos de reatores, neutrinos solares e medições de massa) estabeleceram "cercas" que dizem que esses pássaros invisíveis não deveriam ser capazes de voar da maneira que os dados do Gálio sugerem. Os dados do Gálio querem que os pássaros sejam muito ativos, mas as outras cercas dizem que eles devem ser muito tímidos.

A Conclusão

Os autores não resolveram o mistério; eles o tornaram mais misterioso. Ao usar uma matemática melhor e modelos de elétrons mais precisos, eles confirmaram que os "neutrinos perdidos" são um problema real e robusto, não um erro de cálculo.

Eles concluem que, embora a ideia do "Neutrino Estéril" ainda seja o principal suspeito, ela está atualmente em um impasse com outras evidências experimentais. O mistério permanece sem solução, e os autores sugerem que um novo experimento com um tipo diferente de detector pode ser necessário para finalmente capturar o culpado.

Resumo técnico

Resumo Técnico: Reavaliando a Anomalia do Gálio Usando Funções de Onda Eletrônicas Autoconsistentes

Definição do Problema
O artigo aborda a "anomalia do gálio", uma discrepância persistente observada em experimentos de captura de neutrinos utilizando $^{71}\text{Ga}$ como alvo. Experimentos empregando fontes radioativas de $^{51}\text{Cr}$ e $^{37}\text{Ar}$ (GALLEX, SAGE e, recentemente, BEST) mediram consistentemente taxas de decaimento beta inverso (IBD) induzido por neutrinos sistematicamente menores do que as previsões teóricas. Esse déficit, que agora excede uma significância de $4\sigma$, desafia a robustez dos cálculos de seção de choque teóricos e impulsionou hipóteses de física além do Modelo Padrão, como neutrinos estéreis. Os autores identificam que as avaliações teóricas anteriores basearam-se em aproximações de ordem líder (LO) para funções de onda eletrônicas e tratamentos inconsistentes de elementos de matriz nucleares e leptônicos.

Metodologia
Os autores realizam uma revisão abrangente da seção de choque de captura de neutrinos em $^{71}\text{Ga}$, abandonando as convenções de aproximações LO em favor de uma abordagem numérica autoconsistente. Os componentes metodológicos principais incluem:

1. Soluções Exatas de Função de Onda: Em vez de usar funções de onda aproximadas, os autores resolvem numericamente a equação de Dirac-Coulomb tanto para estados eletrônicos ligados (captura eletrônica, EC) quanto contínuos (IBD). Isso é alcançado através de um pacote de software customizado baseado na biblioteca RADIAL, empregando as equações de Dirac-Hartree-Fock-Slater (DHFS).
2. Potenciais Autoconsistentes: Os cálculos utilizam um potencial DHFS que inclui potenciais nucleares, eletrônicos e de troca. O potencial nuclear é modelado usando uma distribuição de Fermi de dois parâmetros, incorporando a primeira medição direta do raio de carga nuclear do $^{71}\text{Ge}$ ($R_{ch} = 4.032 \pm 0.002$ fm). O parâmetro de potencial de troca ($C_{ex}$) é calculado com uma média de $1.25$ para compensar variações de modelo entre as aproximações de Slater ($3/2$) e Kohn-Sham ($1$).
3. Procedimento de Média: Indo além da convenção de avaliar as funções de onda em um único ponto (por exemplo, $r_0 = 0$ ou $r_0 = R_{box}$), os autores propõem a média das funções de onda eletrônicas e da função de Fermi sobre a densidade de carga nuclear ($\rho_{ch}(r)$). Isso é aplicado tanto ao cálculo da seção de choque de IBD quanto à taxa de captura eletrônica (EC) usada para extrair o elemento de matriz nuclear via princípio do detalhamento do equilíbrio.
4. Entradas Atualizadas: A análise incorpora entradas experimentais recentes, incluindo:
   • Uma média ponderada de quatro medições independentes da meia-vida de EC do $^{71}\text{Ge}$ ($t_{1/2} = 11.465 \pm 0.003$ d).
   • Um valor Q atualizado para o processo de EC ($232.47 \pm 0.09$ keV).
   • Razões de probabilidade de captura eletrônica refinadas ($P_L/P_K$ e $P_M/P_K$) e energias de ligação.
   • Forças Gamow-Teller para estados excitados derivadas de dados de espalhamento $(p,n)$ e $(^3\text{He}, ^3\text{H})$.

Principais Contribuições e Resultados

• Seções de Choque Revisadas: O novo cálculo produz uma seção de choque de estado fundamental aproximadamente 2,8% menor do que as estimativas anteriores devido ao procedimento de média. Ao incluir contribuições de estados excitados, as seções de choque totais para as fontes de $^{51}\text{Cr}$ e $^{37}\text{Ar}$ são encontradas como aproximadamente 2% menores do que as relatadas na literatura mais recente (Ref. [27]).
• Significância da Anomalia: Ao comparar as previsões teóricas atualizadas com os dados experimentais de GALLEX, SAGE e BEST, os autores reavaliam a significância da anomalia:
  • Usando dados de espalhamento $(p,n)$ para contribuições de estados excitados, a discrepância atinge $5.5\sigma$ ($R = 0.835^{+0.030}_{-0.048}$).
  • Usando dados de espalhamento $(^3\text{He}, ^3\text{H})$, a discrepância é de $4.7\sigma$ ($R = 0.811^{+0.037}_{-0.035}$).
  • Uma estimativa conservadora usando apenas a seção de choque do estado fundamental resulta em uma anomalia de $3.8\sigma$.
• Autoconsistência: Os autores enfatizam que, ao utilizar o mesmo software e pressupostos subjacentes tanto para a função de Fermi (IBD) quanto para a densidade eletrônica no núcleo (EC), desvios sistemáticos se cancelam na razão, tornando o resultado mais robusto.

Significância e Interpretação
O artigo afirma que a anomalia do gálio é agora mais significativa do que anteriormente relatado, atingindo até $5.5\sigma$. Os autores reexaminam a interpretação desta anomalia em termos de oscilações de neutrinos ativos-estéreis de base curta (modelo 3+1). Seus resultados confirmam que os dados do gálio, incluindo as precisas medições do BEST, requerem um grande ângulo de mistura ativo-estéril ($0.16 \lesssim \sin^2 2\vartheta_{ee} \lesssim 0.54$ a $2\sigma$ para $\Delta m^2_{41} \gtrsim 0.98 \text{ eV}^2$).

No entanto, os autores observam que essa grande mistura necessária está em tensão com as restrições de dados de antineutrinos de reatores, limites de neutrinos solares e a medição de massa de neutrino do KATRIN. Consequentemente, o artigo conclui que, embora a significância estatística da anomalia tenha aumentado, a explicação de neutrino estéril permanece em tensão com outros dados experimentais. Os autores afirmam que a solução para a anomalia do gálio continua sendo um problema em aberto, potencialmente exigindo um novo experimento de fonte com um detector diferente para resolver a discrepância.