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

Technische Samenvatting: Extractie van de Nucleon Axiale Formfactor uit Lattice QCD met behulp van NNLO Chiral Pertische Perturbatietheorie

Probleemstelling
De nucleale axiale formfactor, $F_A(Q^2)$, is een fundamentele grootheid die de zwakke interacties betrokken bij nucleonen bepaalt, met name de cross-secties van neutrino-nucleonverstrooiing die essentieel zijn voor oscillatie-experimenten. Het extraheren van $F_A(Q^2)$ uit experimentele data is uitdagend vanwege modelafhankelijkheden (bijv. dipool versus $z$-expansie parametrisaties) en experimentele biases. Hoewel Lattice QCD (LQCD) een benadering vanuit de eerste principes biedt, worden simulaties beïnvloed door systematische artefacten, waaronder eindige volume-effecten, discretisatie-afstand en niet-fysische lichte quark-massa's (pion-massa's $M_\pi$). Het bereiken van het fysieke punt vereist een robuuste chirale extrapolatie. Eerdere analyses vertrouwden vaak op niet-relativistische benaderingen of leading-order loop-correcties ($O(p^3)$), wat mogelijk onvoldoende is om de $M_\pi$-afhankelijkheid van de axiale lading ($g_A$) en de radius ($\langle r_A^2 \rangle$) over het bereik van beschikbare LQCD-data te beschrijven.

Methodologie
De auteurs voeren een globale analyse uit van recente LQCD-data voor de axiale formfactor met behulp van Relativistische Baryon Chiral Pertische Perturbatietheorie (B$\chi$PT) tot aan de Next-to-Next-to-Leading Order (NNLO), wat overeenkomt met $O(p^4)$.

• Theoretisch Kader: De berekening maakt gebruik van het Extended On-Mass-Shell (EOMS) renormalisatieschema om de power counting en analytische eigenschappen in aanwezigheid van baryonen te behouden. De $\Delta(1232)$-resonantie wordt behandeld als een expliciete vrijheidsgraad met behulp van de small-scale expansie.
• Lagrangiaan en Diagrammen: De effectieve Lagrangiaan bevat termen tot $O(p^4)$ die betrokken zijn bij nucleon ($N$) en $\Delta$-velden. De berekening incorporeert Feynman-diagrammen tot $O(p^4)$, inclusief bijdragen van $O(p^2)$-vertices en baryon massa-inserties.
• Dataset: De fit maakt gebruik van een gecombineerde set LQCD-resultaten van de RQCD, NME, Mainz, PNDME, PACS en ETMC collaboraties. De analyse richt zich op ensembles met $M_\pi L \geq 3.5$ om eindelijke volume-effecten te verwaarlozen.
• Fit Strategie:
  • Het fit-bereik wordt beperkt tot $M_\pi \leq 402$ MeV en $Q^2 \leq 0.36$ GeV$^2$ (bepaald via $\chi^2$/dof plateaus).
  • Lattice-afstand ($a$) effecten worden fenomenologisch gecorrigeerd met behulp van specifieke afhankelijkheidstermen voor verschillende lattice acties.
  • Truncatiefouten van de chirale expansie worden geschat en in kwadratuur toegevoegd aan de LQCD statistische fouten. Dit gebeurt iteratief om rekening te houden met de correlatie tussen de truncatiefout en de gefitte Low-Energy Constants (LECs).
  • Naturalness Gaussian priors worden toegepast op vrije LEC's om correlaties te verminderen, behalve voor de axiale lading in het chirale limiet ($\mathring{g}_A$).

Kernbijdragen

1. NNLO Relativistische Berekening: De studie biedt de eerste extractie van $F_A(Q^2)$ met behulp van NNLO ($O(p^4)$) relatieve ChPT met expliciete $\Delta$-vrijheidsgraden. Dit is een verbetering ten opzichte van eerdere $O(p^3)$ of niet-relativistische analyses.
2. Systematische Onzekerheidskwantificering: De auteurs modelleren expliciet de truncatiefout van de chirale expansie, waarbij grotere onzekerheden worden toegewezen aan datapunten bij hogere $M_\pi$ en $Q^2$ waar de convergentie slechter is.
3. Rol van de $\Delta$-resonantie: De analyse toont aan dat de expliciete inclusie van de $\Delta$-resonantie noodzakelijk is om de $M_\pi$-afhankelijkheid van de lattice data voor zowel $g_A$ als $\langle r_A^2 \rangle$, evenals de momentum-afhankelijkheid van de formfactor, te reproduceren.
4. Model-onafhankelijke Parametrisatie: Het werk biedt een parametrisatie van $F_A$ die geworteld is in QCD-symmetrieën in plaats van ad-hoc functionele vormen, wat de extrapolatie van LQCD-resultaten naar het fysieke punt vergemakkelijkt.

Resultaten

• Axiale Lading ($g_A$): Op het fysieke punt levert de analyse $g_A = 1.257 \pm 0.011$. Deze waarde komt in goede overeenstemming met de FLAG review ($1.263 \pm 0.010$), maar vertoont een milde spanning ($\sim 1.5\sigma$) met de experimentele waarde van $\beta$-verval ($1.2754 \pm 0.0013$). De auteurs merken op dat niet-verantwoorde radiatieve correcties deze vergelijking kunnen beïnvloeden.
• Axiale Radius ($\langle r_A^2 \rangle$): De geëxtraheerde axiale radius gekwadrateerd is $\langle r_A^2 \rangle = 0.312 \pm 0.037$ fm$^2$. Dit resultaat geeft de voorkeur aan lagere waarden vergeleken met sommige dipool-parametrisaties van neutrino-verstrooiingsdata, maar is consistent met recente $z$-expansie analyses en individuele LQCD-bepalingen.
• Low-Energy Constants (LECs):
  • $d_{16} = -0.82 \pm 0.72$ GeV$^{-2}$ (bepaalt de $M_\pi$-afhankelijkheid van $g_A$).
  • $d_{22} = 1.1 \pm 1.4$ GeV$^{-2}$ (komt binnen de axiale radius).
  • $g_1 = 0.40 \pm 0.46$ (gecorreleerd met $d_{22}$).
  • De grote onzekerheden in $d_{16}$ en $d_{22}$ weerspiegelen de trage convergentie van de chirale reeks voor $g_A(M_\pi)$ en de correlaties tussen de LEC's.
• Convergentie: De $O(p^4)$ berekening met expliciete $\Delta$ biedt een realistische beschrijving van de data tot $M_\pi \approx 400$ MeV en $Q \approx 600$ MeV. De auteurs merken echter op dat de $O(p^4)$ bijdrage aan $g_A$ bij $M_\pi \gtrsim 200$ MeV aanzienlijk is, wat suggereert dat volledige convergentie mogelijk een $O(p^5)$ (two-loop) berekening vereist.

Betekenis
Het artikel claimt een "model-onafhankelijke en systematisch verbeterbare parametrisatie" te bieden van de pion-massa en momentum afhankelijkheid van de axiale formfactor. Door succesvol de LQCD-data te beschrijven tot NNLO met expliciete $\Delta$-bijdragen, valideert dit werk de noodzaak om de $\Delta$-resonantie op te nemen voor nauwkeurige chirale extrapolaties van nucleale eigenschappen. De resulterende waarden voor $g_A$ en $\langle r_A^2 \rangle$ bieden een benchmark voor het verbeteren van voorspellingen van lage-energie zwakke interacties betrokken bij nucleonen, in het bijzonder voor neutrino-oscillatie-experimenten. De auteurs benadrukken dat hun aanpak de modelafhankelijkheid van traditionele dipool-fits vermijdt terwijl het een rigoureuze schatting van theoretische onzekerheden biedt.