Extraction of the nucleon axial form factor from Lattice QCD using NNLO chiral perturbation theory

Resumo Técnico: Extração do Fator de Forma Axial do Núcleon usando Teoria de Perturbação Quiral de NNLO em QCD de Rede

Enunciado do Problema
O fator de forma axial do núcleon, $F_A(Q^2)$, é uma grandeza fundamental que governa interações fracas envolvendo núcleons, particularmente as secções transversais de espalhamento neutrino-núcleon, essenciais para experimentos de oscilação. A extração de $F_A(Q^2)$ a partir de dados experimentais é desafiadora devido a dependências de modelos (por exemplo, parametrizações dipolo vs. expansão $z$) e vieses experimentais. Embora a QCD de Rede (LQCD) ofereça uma abordagem de primeiros princípios, as simulações são afetadas por artefatos sistemáticos, incluindo volume finito, espaçamento de discretização e massas de quarks leves não físicas (massas de pione $M_\pi$). Alcançar o ponto físico requer uma extrapolação quiral robusta. Análises anteriores frequentemente dependiam de aproximações não relativísticas ou correções de loop de ordem líder ($O(p^3)$), que podem ser insuficientes para descrever a dependência de $M_\pi$ da carga axial ($g_A$) e do raio ($\langle r_A^2 \rangle$) através da gama de dados de LQCD disponíveis.

Metodologia
Os autores realizam uma análise global de dados recentes de LQCD para o fator de forma axial usando Teoria de Perturbação Quiral Bariônica Relativística (B$\chi$PT) até a Ordem de Próxima-Próxima-Liderança (NNLO), correspondente a $O(p^4)$.

• Estrutura Teórica: O cálculo emprega o esquema de renormalização Extended On-Mass-Shell (EOMS) para preservar a contagem de potência e as propriedades analíticas na presença de bárions. A ressonância $\Delta(1232)$ é tratada como um grau de liberdade explícito usando a expansão de pequena escala.
• Lagrangiana e Diagramas: A Lagrangiana efetiva inclui termos até $O(p^4)$ envolvendo campos de núcleon ($N$) e $\Delta$. O cálculo incorpora diagramas de Feynman até $O(p^4)$, incluindo contribuições de vértices $O(p^2)$ e inserções de massa bariônica.
• Conjunto de Dados: O ajuste utiliza um conjunto combinado de resultados de LQCD das colaborações RQCD, NME, Mainz, PNDME, PACS e ETMC. A análise foca em ensembles com $M_\pi L \geq 3.5$ para negligenciar efeitos de volume finito.
• Estratégia de Ajuste:
  • O intervalo de ajuste é restrito a $M_\pi \leq 402$ MeV e $Q^2 \leq 0.36$ GeV$^2$ (determinado via platôs de $\chi^2$/dof).
  • Efeitos de espaçamento de rede ($a$) são corrigidos fenomenologicamente usando termos de dependência específicos para diferentes ações de rede.
  • Erros de truncamento da expansão quiral são estimados e adicionados em quadratura aos erros estatísticos de LQCD. Isso é feito iterativamente para considerar a correlação entre o erro de truncamento e as Constantes de Baixa Energia (LECs) ajustadas.
  • Priores Gaussianos de naturalidade são aplicados às LECs livres para reduzir correlações, exceto para a carga axial no limite quiral ($\mathring{g}_A$).

Principais Contribuições

1. Cálculo Relativístico de NNLO: O estudo fornece a primeira extração de $F_A(Q^2)$ usando B$\chi$PT relativística de NNLO ($O(p^4)$) com graus de liberdade de $\Delta$ explícitos. Isso melhora análises anteriores de $O(p^3)$ ou não relativísticas.
2. Quantificação de Incerteza Sistemática: Os autores modelam explicitamente o erro de truncamento da expansão quiral, atribuindo incertezas maiores a pontos de dados em $M_\pi$ e $Q^2$ mais altos, onde a convergência é pior.
3. Papel da Ressonância $\Delta$: A análise demonstra que a inclusão explícita da ressonância $\Delta$ é necessária para reproduzir a dependência de $M_\pi$ dos dados de rede tanto para $g_A$ quanto para $\langle r_A^2 \rangle$, bem como a dependência de momento do fator de forma.
4. Parametrização Independente de Modelo: O trabalho oferece uma parametrização de $F_A$ fundamentada nas simetrias da QCD em vez de formas funcionais ad-hoc, facilitando a extrapolação dos resultados de LQCD para o ponto físico.

Resultados

• Carga Axial ($g_A$): No ponto físico, a análise resulta em $g_A = 1.257 \pm 0.011$. Este valor está em bom acordo com a revisão FLAG ($1.263 \pm 0.010$), mas mostra uma leve tensão ($\sim 1.5\sigma$) com o valor experimental de $\beta$-decaimento ($1.2754 \pm 0.0013$). Os autores observam que correções radiativas não contabilizadas poderiam afetar esta comparação.
• Raio Axial ($\langle r_A^2 \rangle$): O raio axial ao quadrado extraído é $\langle r_A^2 \rangle = 0.312 \pm 0.037$ fm$^2$. Este resultado favorece valores mais baixos comparados a algumas parametrizações dipolo de dados de espalhamento de neutrinos, mas é consistente com análises recentes de expansão $z$ e determinações individuais de LQCD.
• Constantes de Baixa Energia (LECs):
  • $d_{16} = -0.82 \pm 0.72$ GeV$^{-2}$ (governa a dependência de $M_\pi$ de $g_A$).
  • $d_{22} = 1.1 \pm 1.4$ GeV$^{-2}$ (entra no raio axial).
  • $g_1 = 0.40 \pm 0.46$ (correlacionado com $d_{22}$).
  • As grandes incertezas em $d_{16}$ e $d_{22}$ refletem a convergência lenta da série quiral e as correlações entre as LECs.
• Convergência: O cálculo $O(p^4)$ com $\Delta$ explícito fornece uma descrição realista dos dados até $M_\pi \approx 400$ MeV e $Q \approx 600$ MeV. No entanto, os autores observam que a contribuição $O(p^4)$ para $g_A$ em $M_\pi \gtrsim 200$ MeV é considerável, sugerindo que a convergência total pode exigir cálculos de $O(p^5)$ (dois loops).

Significância
O artigo afirma fornecer uma "parametrização independente de modelo e sistematicamente aprimorável" da dependência de massa de pione e de momento do fator de forma axial. Ao descrever com sucesso os dados de LQCD até NNLO com contribuições de $\Delta$ explícitas, o trabalho valida a necessidade de incluir a ressonância $\Delta$ para extrapolações quirais precisas de propriedades de núcleons. Os valores resultantes para $g_A$ e $\langle r_A^2 \rangle$ oferecem um referencial para melhorar as previsões de interações fracas de baixa energia envolvendo núcleons, particularmente para experimentos de oscilação de neutrinos. Os autores enfatizam que sua abordagem evita a dependência de modelo dos ajustes dipolo tradicionais, ao mesmo tempo em que fornece uma estimativa rigorosa das incertezas teóricas.