Proton isovector helicity PDF at NNLO and the twist-3 moment $\tilde{d}_2$ from lattice QCD at physical quark masses

In dit artikel presenteren de auteurs de eerste berekening op het rooster van QCD bij fysische quarkmassa's van de isovector-heliciteit PDF van de proton en de twist-3 moment $\tilde{d}_2$ in het $\overline{\rm MS}$-schema, gebruikmakend van het LaMET-framework met NNLO-nauwkeurigheid.

De kern

Stel je voor dat een proton (het deeltje in de kern van een atoom) een enorm drukke stad is. In deze stad wonen kleine deeltjes: quarks en gluonen. Jarenlang dachten wetenschappers dat ze precies wisten hoe deze stad eruitzag en hoe de inwoners zich bewogen. Maar toen ze de "spin" (een soort interne draaiing) van het proton bestudeerden, ontdekten ze een raadsel: de draaiing van de quarks verklaarde maar een klein deel van de totale draaiing van het proton. Waar zat de rest? Dit heet het "proton-spinraadsel".

Dit wetenschappelijke artikel is als een nieuwe, superkrachtige kaart die we hebben getekend om deze stad beter te begrijpen. Hier is wat de auteurs hebben gedaan, vertaald naar alledaags Nederlands:

1. De Microscoop: Lattice QCD

Om deze deeltjes te zien, kunnen we ze niet gewoon onder een gewone microscoop leggen. Ze zijn te klein en bewegen te snel. In plaats daarvan bouwen de onderzoekers een digitale traliewerk (een "lattice") in een computer. Het is alsof ze de hele stad op een gigantisch schaakbord zetten, waarbij elke vierkantje een stukje ruimte en tijd is. Ze simuleren hoe de quarks en gluonen zich gedragen op dit bord, precies zoals ze dat in de echte natuur doen.

2. De Snelle Auto: LaMET

Een groot probleem is dat we in de computer niet direct kunnen kijken naar hoe de quarks zich bewegen in een "snelle" proton (zoals in een deeltjesversneller). In de computer zitten ze vaak stil.
De auteurs gebruiken een slimme truc genaamd LaMET (Large Momentum Effective Theory).

• De analogie: Stel je voor dat je een foto maakt van een rennende atleet. Als je de camera te langzaam instelt, wordt de atleet wazig. Maar als je de camera razendsnel kunt fotograferen (een "boost" geven aan het proton in de computer), kun je de details van de rennende benen (de quarks) scherp zien.
• Ze hebben het proton in hun simulatie een flinke duw gegeven, zodat het met hoge snelheid (tot 1,53 GeV) door de computer "reed". Hierdoor konden ze de verdeling van de quarks met hun spin (helicity) in kaart brengen.

3. De Kracht van de Kleefband: Renormalisatie

In hun digitale traliewerk zitten echter ook "ruis" en fouten die door de computer zelf worden gegenereerd (zoals wazigheid door de lens).

• De analogie: Stel je voor dat je een foto maakt van een schilderij, maar je hebt een vieze, plakkerige tape (de "Wilson line") over het schilderij geplakt om het vast te houden. Die tape maakt de foto wazig.
• De auteurs hebben een nieuwe manier bedacht om die plakkerige tape te "schoonmaken" (renormalisatie). Ze gebruiken een techniek die lijkt op het wegpoetsen van de plakkerigheid zodat je weer het echte schilderij ziet, zonder de kunstzinnige details te verliezen. Ze hebben zelfs een speciale "reinigingsvloeistof" (renormalon resummation) ontwikkeld om de lastigste vlekken weg te halen.

4. De Verborgen Kracht: De Twist-3 Moment (˜d2)

Naast het tekenen van de kaart van de rennende quarks, hebben ze iets heel speciaals ontdekt: een verborgen kracht.

• De analogie: Stel je voor dat je in de proton-stad een autootje (een quark) hebt dat door de straten rijdt. Er is een onzichtbare wind (de gluonen) die tegen de auto duwt. De onderzoekers hebben gemeten hoe hard die wind de auto duwt.
• Ze hebben voor het eerst precies berekend hoeveel kracht deze "kleur-Lorentz-kracht" is. Het resultaat? De kracht is bijna nul.
• Dit is een groot nieuws! Het betekent dat de quarks in het proton niet zo hard worden "geknepen" of geduwd door de gluonen als sommigen dachten. Het is alsof je ontdekt dat de wind in de stad veel stiller is dan je dacht.

5. Het Eindresultaat: Een Scherpere Kaart

De auteurs hebben twee dingen geproduceerd:

1. Een gedetailleerde kaart: Ze hebben laten zien hoe de quarks met hun spin zich verdelen over de verschillende snelheden in het proton. Hun kaart toont dat er in het midden van de snelheidsverdeling (waar de meeste quarks zitten) meer activiteit is dan eerdere schattingen suggereerden.
2. Een nieuwe wetenschappelijke wet: Ze hebben de "krachtmeting" (˜d2) gedaan die laat zien dat de interne krachten in het proton heel specifiek zijn.

Waarom is dit belangrijk?

Voorheen waren wetenschappers afhankelijk van experimenten met deeltjesversnellers (zoals de LHC of de toekomstige Electron-Ion Collider) om deze kaarten te tekenen. Dat is als proberen de stad te begrijpen door alleen naar de uitlaatgassen van auto's te kijken.
Met deze berekening hebben ze de stad van binnenuit bestudeerd, rechtstreeks vanuit de wetten van de natuurkunde, zonder afhankelijk te zijn van experimenten.

Kort samengevat:
De onderzoekers hebben een supercomputer gebruikt om een virtueel proton te bouwen, het een flinke duw gegeven om de details te zien, en vervolgens de "ruis" weggepoetst. Ze hebben niet alleen een betere kaart gemaakt van hoe de proton draait, maar hebben ook ontdekt dat een bepaalde interne kracht die ze zochten, bijna niet bestaat. Dit helpt ons eindelijk het mysterie van de "proton-spin" op te lossen.

Technische samenvatting

Probleemstelling

De interne structuur van de protonspin is een fundamenteel vraagstuk in de QCD. Hoewel experimenten aangeven dat quark-spin ongeveer 30% bijdraagt aan de totale protonspin, blijft de exacte verdeling van de resterende spin tussen gluonen en orbitale hoekmomenten onzeker. Een cruciale component hiervoor is de heliciteit-partonverdelingsfunctie (PDF), $g_1(x)$, die het verschil beschrijft tussen quarks met spin parallel en antiparallel aan de protonspin.

Daarnaast biedt de twist-3 PDF, $g_T(x) = g_1(x) + g_2(x), inzicht in quark-gluon correlaties. Een specifieke grootheid, de gereduceerde twist-3 moment$\tilde{d}_2$, is evenredig met de gemiddelde transversale kleur-Lorentzkracht die op een quark in de proton werkt. Experimentele bepalingen van $\tilde{d}_2$ zijn echter moeilijk en onnauwkeurig.

De uitdaging voor Gitter-QCD (Lattice QCD) is dat PDFs gedefinieerd zijn op de lichtkegel (niet-lokale operatoren) en dus niet direct in Euclidische ruimtetijd kunnen worden berekend. Hoewel de Large Momentum Effective Theory (LaMET) een oplossing biedt door "quasi-PDFs" te gebruiken, blijven er uitdagingen bestaan:

1. Berekeningen moeten vaak worden uitgevoerd bij fysische quarkmassa's en in de continue limiet, wat rekenkundig zeer intensief is.
2. De matching tussen gitter-data en de MS-scheme (het standaard perturbatieve schema) vereist hoge orde correcties (NNLO) en resummatie van grote logaritmen voor nauwkeurigheid.
3. Het extraheren van twist-3 grootheden zoals $\tilde{d}_2$ is complex door vermenging met lagere-dimensionale operatoren en de noodzaak van niet-perturbatieve aftrekkingen.

Methodologie

De auteurs voeren een berekening uit op een gitter met fysische lichte en vreemde quarkmassa's, gebruikmakend van de volgende technieken:

• Gitter-Setup: Gebruik van een ensemble gegenereerd door de HotQCD samenwerking met $N_f=2+1$ HISQ-acties. De gitterruimte is $64^3 \times 64$ met een rasterafstand van $a = 0.076$ fm.
• Quasi-PDF Matrix-elementen: Berekening van niet-lokale quark-bilineaire operatoren verbonden door een Wilson-lijn in protontoestanden die zijn geboost naar impulsen $P_z = \{0, 0.25, 1.02, 1.53\}$ GeV.
• Renormalisatie: Toepassing van een hybride renormalisatieschema met leading renormalon resummation (LRR). Dit schema elimineert lineaire divergenties in de Wilson-lijn en reguleert de renormalon-ambiguïteit, waardoor een consistente matching naar de MS-scheme mogelijk wordt.
• Operator Product Expansion (OPE): Gebruik van OPE op renormalisatiegroep-invariante verhoudingen van matrix-elementen om de Mellin-momenten van de PDFs direct te extraheren zonder Fourier-transformatie.
• Twist-3 Analyse: Constructie van een specifieke verhouding voor de transversale correlator om de twist-2 bijdrage (Wandzura-Wilczek benadering) te subtraheren en de zuivere twist-3 bijdrage ($\tilde{d}_2$) te isoleren via de $\lambda$-afhankelijkheid (Ioffe-tijd) op korte afstanden.
• Matching en Resummatie: De quasi-PDFs worden perturbatief gematched naar de lichtkegel-PDFs in de MS-scheme met NNLO nauwkeurigheid. Om grote logaritmen te behandelen bij $x \to 0$ en $x \to 1$, worden Renormalization Group (RG) en threshold resummaties toegepast.
• Extrapolatie: Asymptotische extrapolatie van de matrix-elementen naar grote afstanden ($z$) om de Fourier-transformatie naar impulsruimte mogelijk te maken, gevolgd door een modellering van de eindpunt-regio's ($x \to 0, 1$) via Short-Distance Factorization (SDF).

Belangrijkste Bijdragen

1. Eerste berekening bij fysische massa's met NNLO: Dit is een van de eerste berekeningen van de protonheliciteit-PDF die zowel op een fijne gitter bij fysische quarkmassa's wordt uitgevoerd als gebruikmaakt van volledige NNLO-matching en resummatie.
2. Eerste extractie van $\tilde{d}_2$ in de MS-scheme: De auteurs extraheren voor het eerst de twist-3 moment $\tilde{d}_2$ direct in de $\overline{MS}$-scheme vanuit gitter-data, in plaats van in een tussenliggend momentum-subtractieschema. Dit omzeilt problemen met vermenging met lagere-dimensionele operatoren.
3. Geavanceerde theoretische behandeling: Integratie van leading renormalon resummation in de hybride renormalisatie en toepassing van resummatie op de matching-kernen om de theoretische onzekerheden te minimaliseren.

Resultaten

• Mellin-momenten: De auteurs extraheren de eerste paar momenten van de isovector helicity PDF ($\langle x \rangle$ en $\langle x^2 \rangle$) bij een schaal van $\mu = 2$ GeV. De resultaten zijn consistent met, maar iets hoger dan, waarden uit globale fits (zoals JAM22).
  • $\langle x \rangle^{u-d} = 0.291(7)$
  • $\langle x^2 \rangle^{u-d} = 0.101(10)$
• Twist-3 Moment $\tilde{d}_2$: De berekende waarde voor de gereduceerde twist-3 moment is:
  $$\tilde{d}_2^{u-d} (2 \text{ GeV}) = 0.0024(46)$$
  Dit resultaat is statistisch consistent met nul. Dit impliceert dat de zuivere quark-gluon correlatie (de twist-3 bijdrage die niet door de Wandzura-Wilczek relatie wordt beschreven) sterk onderdrukt is in de isovector kanaal op deze schaal.
• PDF Verdeling ($x$-afhankelijkheid): De berekende PDF in het gebied $x \in [0.25, 0.75]$ toont een verhoogde verdeling in het midden-$x$ gebied ($x \in [0.4, 0.75]$) vergeleken met bestaande globale fits. De eindpunt-regio's ($x < 0.25$ en $x > 0.75$) worden gemodelleerd en gekoppeld aan de gitter-data, wat leidt tot een volledige $x$-afhankelijke PDF.
• Vergelijking: De resultaten tonen goede overeenkomst met recente berekeningen in de Coulomb-gauge voor $x < 0.6$, maar vertonen afwijkingen voor $x > 0.6$, wat mogelijk te wijten is aan macht-correcties bij de lagere impulsen die in deze studie worden gebruikt.

Significantie

Deze studie markeert een belangrijke stap in de precisie van gitter-QCD berekeningen van nucleaire structuur:

• Het levert een eerste-principes test van de spin-decompositie van het proton en bevestigt de onderdrukking van specifieke twist-3 quark-gluon correlaties.
• De methode om $\tilde{d}_2$ direct in de MS-scheme te extraheren biedt een robuust alternatief voor eerdere methoden die gevoelig waren voor niet-perturbatieve subtraheringen.
• De toepassing van NNLO-matching en resummatie verlaagt de theoretische onzekerheid aanzienlijk, wat essentieel is voor het interpreteren van toekomstige experimentele data van faciliteiten zoals het Electron-Ion Collider (EIC).
• De bevinding dat de twist-3 bijdrage klein is, ondersteunt het idee dat de Wandzura-Wilczek benadering een goede eerste orde beschrijving is van $g_2(x)$, hoewel de kleine afwijkingen (zoals gemeten in de $\tilde{d}_2$) cruciaal zijn voor een volledig begrip van de nucleaire dynamica.

Kortom, dit werk combineert geavanceerde gittertechnieken met state-of-the-art perturbatieve QCD om een nauwkeuriger en completer beeld te schetsen van de spin-structuur van het proton dan ooit tevoren mogelijk was vanuit eerste principes.