Lattice determination of the neutrino background for $J/ψ\rightarrow γ+ \textrm{invisible}$

Este artigo apresenta o primeiro cálculo de QCD em rede do fundo irredutível do Modelo Padrão para o decaimento $J/\psi \to \gamma + \text{invisível}$, determinando a fração de ramificação para $J/\psi \to \gamma\nu\bar{\nu}$ como sendo $1,00(9)(7) \times 10^{-10}$ para fornecer um referencial crítico para buscas de matéria escura.

O essencial

Imagine que você é um detetive tentando encontrar um fantasma em uma sala muito cheia e barulhenta. O "fantasma" nesta história é a Matéria Escura, uma substância misteriosa que compõe a maior parte do universo, mas que se recusa a interagir com a luz ou com a matéria comum. Cientistas querem vislumbrar essa matéria observando partículas pesadas chamadas J/ψ (pronuncia-se "J-psi") decair. Especificamente, eles estão procurando por um J/ψ que se transforme em um único flash de luz (um fóton) e depois desapareça completamente. Se ele desaparecer, pode ter se transformado em uma partícula de matéria escura.

No entanto, há um problema: Neutrinos.

Os neutrinos são partículas minúsculas e fantasmagóricas que fazem parte do Modelo Padrão da física. Eles também fazem o J/ψ desaparecer no nada quando decai. Para o detector, um neutrino se parece exatamente com a matéria escura. É como tentar encontrar um pássaro raro em uma floresta, mas toda vez que você olha, vê um pombo comum que parece exatamente igual. Se você não souber exatamente quantos pombos existem lá, não poderá ter certeza se encontrou o pássaro raro.

A Missão do Artigo
Este artigo é a primeira vez que cientistas usaram uma simulação matemática superpoderosa (chamada QCD em Rede ou Lattice QCD) para contar exatamente quantos "pombos" (neutrinos) estão escondidos na floresta. Eles queriam calcular a taxa exata na qual um J/ψ se transforma em um fóton e um par de neutrinos ($J/\psi \to \gamma \nu \bar{\nu}$).

Como Eles Fizeram Isso: O "Universo Pixelado"
Para fazer isso, os pesquisadores não usaram um telescópio; eles usaram um computador para construir uma grade 3D (uma rede ou lattice) que representa o espaço e o tempo.

• A Grade: Imagine uma rede de pesca invisível e gigante estendida pelo universo. Eles colocaram a partícula J/ψ nesta rede.
• A Simulação: Eles observaram como o J/ψ interagia com a grade, emitindo um fóton e um par de neutrinos. Como a força forte que mantém o J/ψ unido é incrivelmente complexa (como um novelo de lã emaranhado), eles não podiam usar apenas matemática simples. Eles tiveram que simular o "emaranhado" da lã se dando nós e desatando na grade.
• Limpando o Sinal: Eles tiveram que ser muito cuidadosos para garantir que estivessem vendo apenas o J/ψ e não "ecos" de versões mais pesadas e excitadas da partícula. Eles usaram uma técnica chamada "ajuste de múltiplos estados" (multi-state fit), que é como sintonizar um rádio para filtrar a estática e ouvir apenas a estação clara.
• A Escala: Eles rodaram esta simulação em três tamanhos diferentes de grades (fina, média e grossa) para garantir que seus resultados não fossem apenas um artefato do tamanho da grade. Eles então combinaram matematicamente esses resultados para prever o que aconteceria no mundo real e contínuo.

O Resultado
A equipe calculou a "fração de ramificação" (branching fraction), que é essencialmente a probabilidade de este evento específico acontecer.

• O Número: Eles descobriram que, para cada 10 bilhões de partículas J/ψ que decaem, cerca de 1 delas se transformará em um fóton e neutrinos.
• A Precisão: O cálculo deles é extremamente preciso: $1,00 \times 10^{-10}$. Eles até forneceram uma "margem de erro" para mostrar o quão confiantes estão.

Por Que Isso Importa
O artigo explica que experimentos futuros, como a Instalação Super Tau Charm (STCF), estão sendo construídos para serem tão sensíveis que serão capazes de detectar sinais exatamente neste nível ($10^{-10}$).

Antes deste artigo, os cientistas não tinham um número preciso para o "fundo de neutrinos" (neutrino background). Era como tentar pesar uma pena em uma balança que já estava vibrando com uma quantidade desconhecida de vento. Agora, eles têm uma medição precisa do "vento" (os neutrinos).

A Conclusão
Ao fornecer este número exato, o artigo dá aos experimentalistas um referencial (baseline). Quando eles realizarem seus experimentos no futuro, poderão subtrair esse fundo de neutrinos conhecido de seus dados. Se houver qualquer sinal restante após a subtração dos neutrinos, esse sinal restante pode ser a elusiva Matéria Escura.

Em resumo: Este artigo não encontrou a Matéria Escura, mas construiu a régua perfeita para medir o ruído para que, no futuro, possamos finalmente ouvir o sussurro do escuro.

Resumo técnico

Resumo Técnico: Determinação em Rede do Fundo de Neutrinos para $J/\psi \to \gamma + \text{invisível}$

Definição do Problema
A busca por matéria escura leve e neutrinos estéreis através de decaimentos radiativos de quarks pesados (ex: $J/\psi \to \gamma + \text{invisível}$) é um objetivo primário na física de partículas moderna. Embora os limites superiores experimentais para essas frações de ramificação tenham melhorado significamente, instalações futuras como o Super Tau Charm Facility (STCF) projetam alcançar sensibilidades na ordem de $10^{-10}$. Nesta sensibilidade, o fundo irredutível do Modelo Padrão (SM) proveniente do decaimento $J/\psi \to \gamma \nu \bar{\nu}$ torna-se dominante. Para distinguir potenciais sinais de nova física do fundo do SM, é indispensável uma determinação precisa e independente de modelo do fundo de neutrinos do SM. Estimativas fenomenológicas anteriores basearam-se em aproximações (ex: modelos de singlete de cor não relativísticos) que introduziam incertezas sistemáticas não controladas.

Metodologia
Este trabalho apresenta o primeiro cálculo ab initio de QCD em rede (Lattice QCD) da largura de decaimento $J/\psi \to \gamma \nu \bar{\nu}$. A metodologia emprega as seguintes estratégias para garantir precisão e controle sobre os efeitos sistemáticos:

1. Fatoração e Formadores de Momento (Form Factors): A amplitude de decaimento é fatorada em uma parte leptônica perturbativa e uma parte hadrônica não perturbativa. A interação hadrônica é codificada em uma função $H_{\mu\nu\alpha}$, que é parametrizada por um único formador de momento $F_{\gamma\nu\bar{\nu}}(q^2)$.
2. Método da Função Escalar: Em vez de calcular os formadores de momento em momentos discretos de rede e realizar interpolação (o que introduz dependência de modelo), os autores utilizam o "método da função escalar". Isso envolve a construção de uma função escalar $I(E_\gamma, \Delta t)$ integrando a função de correlação de três pontos euclidiana sobre as coordenadas espaciais com um núcleo de função de Bessel ($j_0$). Isso permite a extração do formador de momento com cobertura total da distribuição de $q^2$ no espaço de fase.
3. Supressão de Estados Excitados: Para mitigar a contaminação de estados excitados, os autores empregam um ajuste de múltiplos estados usando um ansatz de dois estados. Eles analisam a dependência temporal das funções de correlação e extraem a contribuição do estado fundamental no limite de grande separação temporal ($\Delta t \to \infty$).
4. Controle de Volume Finito e Discretização:
   • Volume Finito: Uma integral de volume espacial com um intervalo de truncamento $R$ é utilizada para monitorar os efeitos de volume finito. Os autores demonstram que esses efeitos são exponencialmente suprimidos e negligenciáveis para o sistema de charmonium.
   • Extrapolação para o Contínuo: Os cálculos são realizados em três conjuntos de gauge com diferentes espaçamentos de rede ($a \approx 0.067, 0.085, 0.098$ fm) gerados pela Extended Twisted Mass Collaboration (ETMC). Os resultados são extrapolados para o limite contínuo ($a \to 0$) assumindo um comportamento linear em $a^2$, consistente com a melhoria automática $O(a)$ de fermiões de massa torcida (twisted-mass).
5. Renormalização: As correntes eletromagnética e fraca são renormalizadas usando constantes $Z_V$ e $Z_A$ determinadas em trabalhos anteriores e via esquema RI-MOM.

Principais Resultados
Os autores computam a razão adimensional $R_f = \Gamma_{\gamma\nu\bar{\nu}} / f_{J/\psi}$ e realizam a extrapolação para o contínuo para obter a fração de ramificação física. O resultado final é:
$$\text{Br}(J/\psi \to \gamma \nu \bar{\nu}) = 1.00(9)(7) \times 10^{-10}$$
onde a primeira incerteza é estatística (incluindo incertezas da extrapolação para o contínuo) e a segunda é a incerteza sistemática estimada.

• Comparação com Modelos: Este resultado é comparado com uma previsão fenomenológica anterior baseada no modelo de singlete de cor não relativístico, que estimou $0.7 \times 10^{-10}$. O resultado de rede fornece um benchmark mais preciso e livre de parâmetros, diferindo ligeiramente, mas permanecendo na mesma ordem de magnitude.
• Verificações Sistemáticas: Os autores verificaram a robustez de seu resultado excluindo o conjunto de rede mais grosseiro ($a_{98}$) da extrapolação, obtendo $1.07(13) \times 10^{-10}$. A consistência entre os conjuntos de dados total e reduzido apoia a ausência de erros de discretização $O(a)$ residuais significativos.
• Contribuições Negligenciadas: Diagramas desconectados são observados como negligenciáveis devido à supressão Okubo-Zweig-Iizuka (OZI) no sistema de charmonium.

Significância e Alegações
O artigo afirma entregar o primeiro benchmark de QCD em rede ab initio para o fundo de neutrinos em decaimentos $J/\psi \to \gamma + \text{invisível}$. Sua significância reside em:

1. Viabilização de Buscas Futuras: Ao fornecer uma previsão precisa (nível de $10^{-10}$), este trabalho permite que experimentos como o STCF subtraiam com precisão o fundo do SM, permitindo que investiguem a nova física no setor invisível sem serem limitados pela incerteza teórica.
2. Generalizabilidade: A metodologia é apresentada como generalizável para outros canais de quarkonium, tais como $\Upsilon \to \gamma + \text{invisível}$ e $\phi \to \gamma + \text{invisível}$, fornecendo um roteiro teórico para buscas de matéria escura em diferentes escalas de energia.
3. Rigor Teórico: O cálculo evita as aproximações inerentes aos modelos fenomenológicos, oferecendo uma determinação não perturbativa rigorosa da largura de decaimento, livre de incertezas sistemáticas não controladas associadas a pressupostos de modelos.