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

Sintesi Tecnica: Estrazione del Fattore di Forma Assiale del Nucleone tramite Lattice QCD utilizzando la Teoria della Perturbazione Chirale al NNLO

Problematica
Il fattore di forma assiale del nucleone, $F_A(Q^2)$, è una quantità fondamentale che governa le interazioni deboli coinvolgenti i nucleoni, in particolare le sezioni d'urto di scattering neutrino-nucleone, essenziali per gli esperimenti di oscillazione. Estrarre $F_A(Q^2)$ dai dati sperimentali è complesso a causa delle dipendenze dai modelli (ad esempio, parametrizzazioni dipolo vs. espansione $z$) e dei bias sperimentali. Sebbene la QCD su reticolo (LQCD) offra un approccio basato sui primi principi, le simulazioni sono influenzate da artefatti sistematici, tra cui il volume finito, l'intervallo di discretizzazione e masse dei quark leggeri non fisiche (masse del pione $M_\pi$). Raggiungere il punto fisico richiede un'estrapolazione chirale robusta. Le analisi precedenti spesso si sono affidate ad approssimazioni non relativistiche o a correzioni ai loop di ordine dominante ($O(p^3)$), che potrebbero essere insufficienti per descrivere la dipendenza di $M_\pi$ per la carica assiale ($g_A$) e il raggio ($\langle r_A^2 \rangle$) nell'intervallo dei dati LQCD disponibili.

Metodologia
Gli autori eseguono un'analisi globale dei recenti dati LQCD per il fattore di forma assiale utilizzando la Teoria della Perturbazione Chirale Relativistica (B$\chi$PT) fino al Next-to-Next-to-Leading Order (NNLO), corrispondente a $O(p^4)$.

• Framework Teorico: Il calcolo impiega lo schema di rinormalizzazione Extended On-Mass-Shell (EOMS) per preservare il conteggio delle potenze e le proprietà analitiche in presenza di barioni. La risonanza $\Delta(1232)$ è trattata come un grado di libertà esplicito utilizzando la small-scale expansion.
• Lagrangiana e Diagrammi: La Lagrangiana efficace include termini fino a $O(p^4)$ che coinvolgono i campi del nucleone ($N$) e del $\Delta$. Il calcolo incorpora diagrammi di Feynman fino a $O(p^4)$, incluse le contribuzioni da vertici $O(p^2)$ e inserzioni di massa barionica.
• Set di Dati: Il fit utilizza un set combinato di risultati LQCD delle collaborazioni RQCD, NME, Mainz, PNDME, PACS ed ETMC. L'analisi si concentra su ensemble con $M_\pi L \geq 3.5$ per trascurare gli effetti di volume finito.
• Strategia di Fit:
  • L'intervallo di fit è limitato a $M_\pi \leq 402$ MeV e $Q^2 \leq 0.36$ GeV$^2$ (determinato tramite plateau $\chi^2$/dof).
  • Gli effetti del passo di reticolo ($a$) sono corretti fenomenologicamente utilizzando termini di dipendenza specifici per diverse azioni di reticolo.
  • Gli errori di troncamento dell'espansione chirale sono stimati e aggiunti in quadratura agli errori statistici LQCD. Ciò avviene iterativamente per tenere conto della correlazione tra l'errore di troncamento e le Costanti a Bassa Energia (LEC) adattate.
  • Sono applicati priori Gaussiani di naturalness per ridurre le correlazioni, eccetto per la carica assiale nel limite chirale ($\mathring{g}_A$).

Contributi Chiave

1. Calcolo Relativistico NNLO: Lo studio fornisce la prima estrazione di $F_A(Q^2)$ utilizzando la B$\chi$PT relativistica al NNLO ($O(p^4)$) con gradi di libertà $\Delta$ espliciti. Questo migliora rispetto alle precedenti analisi $O(p^3)$ o non relativistiche.
2. Quantificazione dell'Incertezza Sistematica: Gli autori modellano esplicitamente l'errore di troncamento dell'espansione chirale, assegnando incertezze maggiori ai punti dati a $M_\pi$ e $Q^2$ più elevati, dove la convergenza è peggiore.
3. Ruolo della Risonanza $\Delta$: L'analisi dimostra che l'inclusione esplicita della risonanza $\Delta$ è necessaria per riprodurre la dipendenza di $M_\pi$ dei dati del reticolo sia per $g_A$ che per $\langle r_A^2 \rangle$, nonché la dipendenza dal momento del fattore di forma.
4. Parametrizzazione Indipendente dal Modello: Il lavoro offre una parametrizzazione di $F_A$ radicata nelle simmetrie della QCD piuttosto che in forme funzionali ad hoc, facilitando l'estrapolazione dei risultati LQCD al punto fisico.

Risultati

• Carica Assiale ($g_A$): Al punto fisico, l'analisi fornisce $g_A = 1.257 \pm 0.011$. Questo valore è in buon accordo con la revisione FLAG ($1.263 \pm 0.010$), ma mostra una lieve tensione ($\sim 1.5\sigma$) con il valore sperimentale dal decadimento $\beta$ ($1.2754 \pm 0.0013$). Gli autori osservano che correzioni radiative non tenute conto potrebbero influenzare questo confronto.
• Raggio Assiale ($\langle r_A^2 \rangle$): Il quadrato del raggio assiale estratto è $\langle r_A^2 \rangle = 0.312 \pm 0.037$ fm$^2$. Questo risultato favorisce valori più bassi rispetto ad alcune parametrizzazioni dipolo dei dati di scattering dei neutrini, ma è coerente con recenti analisi di espansione $z$ e determinazioni LQCD individuali.
• Costanti a Bassa Energia (LECs):
  • $d_{16} = -0.82 \pm 0.72$ GeV$^{-2}$ (governa la dipendenza da $M_\pi$ di $g_A$).
  • $d_{22} = 1.1 \pm 1.4$ GeV$^{-2}$ (entra nel raggio assiale).
  • $g_1 = 0.40 \pm 0.46$ (correlata con $d_{22}$).
  • Le grandi incertezze in $d_{16}$ e $d_{22}$ riflettono la lenta convergenza della serie chirale per $g_A(M_\pi)$ e le correlazioni tra le LEC.
• Convergenza: Il calcolo $O(p^4)$ con $\Delta$ esplicito fornisce una descrizione realistica dei dati fino a $M_\pi \approx 400$ MeV e $Q \approx 600$ MeV. Tuttavia, gli autori notano che la contribuzione $O(p^4)$ a $g_A$ per $M_\pi \gtrsim 200$ MeV è significativa, suggerendo che la piena convergenza potrebbe richiedere calcoli $O(p^5)$ (a due loop).

Significatività
Il paper sostiene di fornire una "parametrizzazione indipendente dal modello e sistematicamente migliorabile" della dipendenza dalla massa del pione e dal momento del fattore di forma assiale. Descrivendo con successo i dati LQCD fino a NNLO con contributi $\Delta$ espliciti, il lavoro valida la necessità di includere la risonanza $\Delta$ per accurate estrapolazioni chirali delle proprietà nucleoniche. I valori risultanti per $g_A$ e $\langle r_A^2 \rangle$ offrono un benchmark per migliorare le previsioni delle interazioni deboli a bassa energia che coinvolgono i nucleoni, in particolare per gli esperimenti di oscillazione dei neutrini. Gli autori sottolineano che il loro approccio evita la dipendenza dal modello delle tradizionali fit dipolo, fornendo al contempo una stima rigorosa delle incertezze teoriche.