Lattice Calculation of Short-Range Contributions to Neutrinoless Double-Beta Decay $\pi^-\to\pi^+ ee$ at Physical Pion Mass

Dit artikel presenteert een berekening van de $\pi^- \to \pi^+ ee$-matrixelementen met behulp van rooster-QCD bij de fysische pionmassa met domeinwandfermionen, waarbij een nieuwe methode wordt toegepast om rond-de-wereld-effecten af te trekken en niet-perturbatieve renormalisatie wordt gebruikt om onzekerheden te verminderen en discrepanties in eerdere studies van kortafstandscontributies bij neutrinoloze dubbel-bètaverval op te lossen.

De kern

Het Grote Geheel: Een Kosmisch Mysterie Oplossen

Stel je het universum voor als een gigantische puzzel, en een van de ontbrekende stukjes is de aard van neutrino's. Dit zijn kleine, spookachtige deeltjes die door allesheen razen. Wetenschappers willen weten: zijn het "Dirac"-fermionen (zoals gewone elektronen, waarbij een deeltje en zijn antideeltje verschillend zijn) of "Majorana"-fermionen (waarbij een deeltje zijn eigen antideeltje is)?

De enige manier om dit mysterie op te lossen, is wachten op een zeer zeldzaam fenomeen dat neutrinoloze dubbel-bèta-verval wordt genoemd. Het is alsof je twee atomen ziet die spontaan veranderen in andere atomen en twee elektronen uitstoten, maar zonder dat er neutrino's vrijkomen. Als we dit zien gebeuren, bewijst het dat neutrino's hun eigen antideeltjes zijn.

Het Probleem: Een Ruig Signaal

Om te voorspellen of dit fenomeen zal plaatsvinden, moeten natuurkundigen zware wiskunde uitvoeren. Ze splitsen de berekening op in twee delen:

1. Het Lange-afstandsdeel: Alsof een fluistering door een kamer reist.
2. Het Korte-afstandsdeel: Alsof er direct naast je oor wordt geschreeuwd.

Dit artikel richt zich op het Korte-afstandsdeel. Specifiek berekenen ze hoe twee pionen (deeltjes gemaakt van quarks) met elkaar interageren om twee elektronen te produceren. Denk hierbij aan het meten van de "luidheid" van die schreeuw.

Het Conflict: Twee verschillende teams van wetenschappers hadden eerder geprobeerd deze "luidheid" te meten met supercomputers (zogenaamde Lattice QCD). Hun resultaten waren echter met een factor twee van elkaar verschillend. Het was alsof het ene team zei dat de schreeuw 60 decibel was, en het andere team dat het 120 decibel was. Deze enorme onenigheid maakte het moeilijk om de voorspellingen voor het neutrino-mysterie te vertrouwen.

De Oplossing: Een Nieuwe Manier om de Data Schoon te Maken

De auteurs van dit artikel besloten hun eigen experiment te runnen om de score te beslechten. Ze gebruikten een enorme supercomputer om de subatomaire wereld te simuleren. Maar ze stonden voor een specifiek technisch probleem: "Wereld-om-effecten".

De Analogie: Stel je voor dat je een gesprek opneemt in een kleine, echo-rijke kamer met een ronde muur. Als je in je handen klapt, reist het geluid vooruit, botst het tegen de muur, wikkelt het zich om de kamer en komt het van achteren bij je terug. In de computersimulatie is de "kamer" het rooster van ruimte-tijd. Omdat het rooster eindig is, kunnen de deeltjes helemaal om de lus heen reizen en de meting verstoren, waardoor een verwarrende "echo" ontstaat die de data bederft.

De Innovatie: Eerdere methoden probeerden te raden hoe ze deze echo's konden opheffen. Dit team bedacht een nieuwe aftrekmethode.

• In plaats van te raden, isoleerden ze het "echo"-signaal direct uit de data.
• Ze berekenden precies hoe sterk de echo was en trokken deze af van het hoofdsignaal.
• Het Resultaat: Het "ruis" verdween, waardoor een schoon, stabiel signaal (een "plateau") overbleef dat ze konden vertrouwen.

De Verificatie: De Liniaal Controleren

Om zeker te weten dat hun nieuwe methode niet kapot was, controleerden ze hun werk tegen een bekende standaard. Ze berekenden een specifieke waarde (een "zakparameter") die eerder door andere teams was gemeten.

• Hun resultaat kwam perfect overeen met de vertrouwde standaard.
• Toen ze hun resultaat vergeleken met het team dat het "factor van twee"-verschil had, ontdekten ze dat hun getallen precies dubbel zo groot waren als die van het andere team.
• De Conclusie: Het blijkt dat het andere team waarschijnlijk een iets andere "liniaal" (normalisatieconventie) had gebruikt voor hun metingen. Zodra je dat verschil in aanmerking neemt, vallen de datapunten eigenlijk perfect op elkaar. De methode van de auteurs bevestigt dat hun berekening correct is en lost de verwarring op.

Het Eindresultaat

Het team slaagde erin om de "kortebereik"-bijdrage aan het neutrinoloze dubbel-bèta-vervalproces te berekenen met veel hogere precisie dan voorheen.

• Ze verwijderden de "echo's" (wereld-om-effecten) die eerdere data verstoorden.
• Ze gebruikten twee verschillende wiskundige "lenzen" (renormalisatieschema's) om ervoor te zorgen dat hun wiskunde stevig was.
• Ze leverden een definitief, hoogprecisie getal dat natuurkundigen helpt te voorspellen of we dit zeldzame verval uiteindelijk in echte experimenten zullen zien.

Kortom: Ze bouwden een betere microscoop, maakten het statische ruis schoon en bevestigden dat het eerdere meningsverschil slechts een kwestie was van het gebruik van verschillende meetlinten. Nu heeft de wetenschappelijke gemeenschap een betrouwbaar getal om het mysterie van de neutrino-massa op te lossen.

Technische samenvatting

1. Probleemstelling

Neutrinoloze dubbel-bèta ($0\nu\beta\beta$) verval is een kritiek proces om te bepalen of neutrino's Majorana-fermionen zijn en om Leptongetal-schending (LNV) te onderzoeken. De vervalsnelheid hangt af van nucleaire matrixelementen (NME's). In de Standard Model Effective Field Theory (SMEFT) ontstaan kortafstandscontributies tot $0\nu\beta\beta$ uit dimensie-9 operatoren. Deze bijdragen worden gemedieerd door het proces $\pi^- \to \pi^+ ee$ via pionuitwisseling tussen nucleonen.

De Kernkwestie: Eerdere rooster-QCD-berekeningen van de kortafstands-matrixelementen voor $\pi^- \to \pi^+ ee$ (specifiek Refs. [30] en [31]) rapporteerden aanzienlijk discrepante resultaten. Voor de zakparameter $B_\pi^{MS}$ bij $\mu=3$ GeV rapporteerde Ref. [30] $\approx 0.42$, terwijl Ref. [31] $\approx 0.20$ rapporteerde—een discrepantie van ongeveer een factor twee. Deze inconsistentie belemmert betrouwbare theoretische voorspellingen voor $0\nu\beta\beta$ vervalsnelheden. Bovendien leden eerdere studies aan:

• Systematische Onzekerheden: Significante "around-the-world" (ATW) effecten door pion-achterwaartse propagatie in periodieke temporele roosters, wat de extractie van grondtoestands-matrixelementen verstoort.
• Chirale Extrapolatie: Sommige studies maakten gebruik van onfysische pionmassa's, wat extrapolatie naar het fysische punt vereiste, waardoor extra modelafhankelijkheid werd geïntroduceerd.

2. Methodologie

De auteurs voerden een hoogprecisie rooster-QCD-berekening uit met behulp van Domain Wall Fermion (DWF) ensemble's gegenereerd door de RBC/UKQCD-samenwerking.

A. Roosteropstelling

• Ensemble's: Gebruik gemaakt van $N_f = 2+1$ ensemble's, waaronder twee met fysische pionmassa's (48I en 64I) met verschillende roosterafstanden ($a^{-1} \approx 1.73$ en $2.36$ GeV) om een directe continuüm-extrapolatie mogelijk te maken. Onfysische massa-ensemble's werden ook gebruikt voor kruiscontroles.
• Operatoren: Berekening van matrixelementen voor vijf dimensie-9 vier-kwark operatoren ($O_1, O_2, O_3, O'_1, O'_2$) relevant voor de $\pi^- \to \pi^+ ee$ overgang.

B. Mitigatie van "Around-the-World" (ATW) Effecten

Een belangrijke innovatie van dit werk is een nieuwe methode om ATW-effecten te behandelen, die ontstaan wanneer pions zich voortplanten over de periodieke temporele grens, waardoor correlatiefuncties worden verontreinigd.

• Analyse: De auteurs analyseerden de temporele structuur van drie-punts correlatiefuncties, waarbij bijdragen werden geïdentificeerd van het gewenste diagram (A) en ATW-diagrammen (B, C, D).
• Nieuwe Subtraktiemethode: In plaats van te vertrouwen op verhoudingsmethoden (R1 of R2) die ATW-effecten slechts gedeeltelijk onderdrukken, stelden de auteurs een directe subtraktietechniek voor.
  • Zij identificeerden tijdsregio's waar de ATW-bijdrage van een enkele pion die de grens kruist (Diagram B) domineert.
  • Zij reconstrueerden de ATW-bijdrage direct uit de roosterdata met behulp van specifieke tijdscheidingen ($t_{\pi\pi}, t_\pi$).
  • Formule: Het gecorrigeerde matrixelement wordt verkregen door de gereconstrueerde ATW-term af te trekken:
    $$O^{(sub)}_i(t) = N_\pi^{-2} e^{2m_\pi t} \left[ C^i_3(t,t) - 2 C^i_3(t_{\pi\pi}, T - t_{\pi\pi} - t_\pi) e^{-m_\pi(T-t_\pi-2t_{\pi\pi})} \right]$$
  • Deze methode slaagde erin stabiele plateaus in de effectieve massa-plots te herstellen, wat de systematische fouten aanzienlijk verminderde in vergelijking met eerdere verhoudingsmethoden.

C. Renormalisatie

• Systeem: Niet-perturbatieve renormalisatie werd uitgevoerd met behulp van de RI/SMOM (Rome-Southampton) methode met niet-exceptionele kinematica.
• Schema's: Twee schema's werden ingezet om perturbatieve truncatiefouten te schatten: $(\gamma_\mu, \gamma_\mu)$ en $(\not{q}, \not{q})$.
• Matching: De resultaten werden omgezet naar het $\overline{MS}$-schema bij $\mu = 3$ GeV met behulp van één-lus matching-coëfficiënten. Het verschil tussen de twee RI/SMOM-schema's werd gebruikt om de systematische onzekerheid door perturbatieve truncatie te schatten.

3. Belangrijkste Resultaten

A. Blote Matrixelementen

De studie extraheerde blote matrixelementen $\langle \pi^+ | O_i | \pi^- \rangle$ bij de fysische pionmassa.

• Kruiscontrole met Zakparameters: De auteurs verifieerden hun normalisatie door het matrixelement van $O_3$ (gerelateerd aan de neutrale meson zakparameter $B_K$) te vergelijken met eerdere door FLAG-beoordeelde resultaten. De resultaten kwamen perfect overeen met Ref. [38] (die dezelfde normalisatieconventie gebruikte), maar waren ongeveer twee keer zo groot als de resultaten in Ref. [31].
• Normalisatie-discrepantie: De verhouding van de resultaten van de auteurs tot Ref. [31] is consistent $\approx 2.0$ voor alle vijf operatoren. Dit suggereert een globale normalisatiefactorfout in Ref. [31], en niet een fysische discrepantie.

B. Gerenormaliseerde Matrixelementen en Zakparameter

De uiteindelijke resultaten in het $\overline{MS}$-schema ($\mu = 3$ GeV) zijn:

• Zakparameter: $B_\pi^{MS}(3 \text{ GeV}) = 0.3769(24)_{stat}(76)_{PT}(1)_L$.
  • Deze waarde is consistent met Ref. [30] ($0.421$) binnen de onzekerheden, maar significant verschillend van Ref. [31] ($0.197$).
  • De verhouding tot Ref. [31] is $1.91(18)$, wat de normalisatieproblematiek van een factor 2 bevestigt.
• Continuüm-extrapolatie: Met behulp van de ensemble's met fysische massa (48I en 64I) voerden de auteurs een directe continuüm-extrapolatie uit, waardoor onzekerheden geassocieerd met chirale extrapolatie werden geëlimineerd.
• Onzekerheidsreductie: De totale onzekerheden in dit werk zijn aanzienlijk verminderd in vergelijking met eerdere studies door:
  1. Directe berekening bij fysische pionmassa.
  2. Verbeterde statistische precisie.
  3. Effectieve verwijdering van ATW-systematische fouten.
  4. Robuuste schatting van perturbatieve truncatiefouten.

4. Betekenis en Conclusie

• Oplossing van Discrepantie: Dit werk biedt een onafhankelijke kruiscontrole die sterk de normalisatie ondersteunt die wordt gebruikt in Ref. [30] en de FLAG-review, terwijl een waarschijnlijke normalisatiefout in Ref. [31] wordt geïdentificeerd. De discrepantie van een factor 2 wordt opgelost als een kwestie van normalisatieconventie in plaats van een fundamenteel fysisch meningsverschil.
• Methodologische Vooruitgang: De voorgestelde subtraktiemethode voor ATW-effecten is een robuuste techniek die kan worden toegepast op andere roosterberekeningen die lichte mesonen in periodieke volumes betreffen, waardoor de betrouwbaarheid van toekomstige $0\nu\beta\beta$-matrixelementberekeningen wordt verbeterd.
• Impact op de Fysica: Door precieze, continuüm-geëxtrapoleerde matrixelementen bij fysische pionmassa's te leveren, vermindert dit werk de niet-perturbatieve QCD-onzekerheden in de theoretische voorspelling van $0\nu\beta\beta$ vervalsnelheden. Dit is cruciaal voor het interpreteren van experimentele resultaten van experimenten van de volgende generatie (bijv. LEGEND, nEXO, CUPID) om de aard van neutrino-massa en de schaal van nieuwe fysica te bepalen.

Samenvattend stelt dit artikel een nieuwe benchmark vast voor kortafstandscontributies tot neutrinoloze dubbel-bèta verval, waarbij eerdere inconsistenties worden opgelost door middel van rigoureuze methodologie, simulaties met fysische massa's en een nieuwe behandeling van eindige-volume-artefacten.