Efficient calculation of two-neutrino double-beta-decay nuclear matrix elements

Technische Samenvatting: Efficiënte Berekening van Twee-Neutrino Dubbel-Bèta-Zerval-Nucleaire Matrixelementen

Probleemstelling
Betrouwbare nucleaire matrixelementen (NME's) zijn cruciaal voor het interpreteren van dubbel-bètaverval ($\beta\beta$) experimenten en het koppelen van gemeten halfwaardetijden aan fundamentele zwakke-interactiefysica. Terwijl de twee-neutrino modus ($2\nu\beta\beta$) een proces is binnen het Standaardmodel dat in verschillende kernen is waargenomen, blijft de neutrinovrije modus ($0\nu\beta\beta$) het primaire signaal voor schending van het leptongetal. Recente statistische shell-model analyses wijzen op een sterke correlatie tussen $2\nu\beta\beta$- en $0\nu\beta\beta$-NME's, wat suggereert dat nauwkeurige berekeningen van de eerste de laatstgenoemde kunnen beperken.

Het berekenen van $2\nu\beta\beta$-NME's ($M_{2\nu}$) vormt echter een aanzienlijke computationele uitdaging. De standaardaanpak vereist een expliciete sommatie over alle $1^+$-toestanden in de oneven-even kern. In grote modelruimtes (bijv. zware kernen of uitgebreide valentieruimtes zoals $sd-pf$) wordt het aantal bijdragende toestanden prohibitief groot, en de lopende som kan traag of niet-monotoon convergeren. Volledige diagonalisatie van de veel-deeltjes-Hamiltoniaan om deze toestanden te verkrijgen is vaak computationeel onhaalbaar, met name voor experimenten van de volgende generatie die streven naar sensitiviteiten tot $10^{28}$ jaar.

Methodologie
Het artikel presenteert een verbeterde implementatie van de sterktefunctie-methode gebaseerd op Lanczos-iteraties om $M_{2\nu}$ te berekenen zonder volledige diagonalisatie te vereisen. De aanpak volgt de formalisering van Ref. [24] en maakt gebruik van een numeriek voorschrift zoals gedetailleerd in Ref. [25].

1. Lanczos-iteraties: De methode start Lanczos-iteraties met behulp van Gamow-Teller (GT) "doorweg"-toestanden (doorweg states), geconstrueerd door de GT-operator toe te passen op de initiële ($0^+_i$) of finale ($0^+_f$) grondtoestanden.
2. Tridiagonaal Matrix: In plaats van de volledige eigenvectoren van de veel-deeltjes-Hamiltoniaan op te slaan, construeert het algoritme een veel kleinere tridiagonale Lanczos-matrix. De NME wordt benaderd door te sommeren over de eigenwaarden en eigenvectoren van deze kleine matrix.
3. Coëfficiënten: De expansiecoëfficiënten ($g^\pm_m$) worden afgeleid van de Lanczos-vectoren en de energiedenominatoren geassocieerd met de tussenliggende toestanden.
4. Efficiëntie: Deze aanpak vermijdt de opslag- en verwerkingskosten van het diagonaliseren van de volledige Hamiltoniaan. Het monitort convergentie iteratief met behulp van de Lancsjos-vectoren zelf, wat een computationele efficiëntiewinst van ongeveer een factor twee oplevert ten opzichte van eerdere sterktefunctie-implementaties.
5. Extensie naar Hogere Ordes: Het raamwerk wordt uitgebreid om hogere-orde NME's ($M_{2\nu-3}$ en $M_{2\nu-5}$) te berekenen die vereist zijn voor Taylor-expansie behandelingen van fase-ruimtefactoren en gerelateerde twee-neutrino elektron-vangst observabelen. Deze termen bevatten hogere machten van de energiedenominator, waardoor ze gevoeliger zijn voor laag-energetische toestanden.

Belangrijkste Bijdragen en Resultaten
De auteur past deze methode toe op verschillende experimenteel belangrijke isotopen ($^{48}\text{Ca}$, $^{76}\text{Ge}$, $^{82}\text{Se}$, $^{128,130}\text{Te}$, en $^{136}\text{Xe}$) met behulp van standaard effectieve Hamiltoniaanse (bijv. GXPF1A, KB3G, JUN45, GCN2850, GCN5082, SVD) in de overeenkomstige valentieruimtes.

• Benchmarking: De methode wordt rigoureus getest door de "lopende sommen" verkregen uit de Lanczos-expansie te vergelijken met expliciete sommaties over honderden of duizenden tussenliggende $1^+$-toestanden.
  • Voor $^{48}\text{Ca}$ (pf-schil) reproduceert de 25-iteratie Lanczos-resultaat de directe som over 1000 toestanden, waarbij de snelle stijging, veranderingen in sterktegebieden en het uiteindelijke stabiele plateau worden gevangen.
  • Voor zwaardere kernen zoals $^{82}\text{Se}$ en $^{136}\text{Xe}$, waar directe diagonalisatie van duizenden toestanden moeilijk is, reproduceert de Lanczos-methode het convergentiegedrag van directe sommen (bijv. over 2000 toestanden voor $^{136}\text{Xe}$).
  • De methode slaagt erin de bijdragen van de GT-reusachtige resonantie (giant resonance) regio te vangen en gaat om met de niet-triviale annuleringen die typerend zijn voor $2\nu\beta\beta$ matrixelementen.
• Hogere-orde Termen: Benchmarks voor $M_{2\nu-3}$ en $M_{2\nu-5}$ in $^{48}\text{Ca}$ tonen uitstekende overeenstemming met expliciete sommen. Dit bevestigt dat dezelfde Lanczos-coëfficiënten gebruikt voor de leidende NME gegeneraliseerd kunnen worden naar hogere inverse machten van de energiedenominator.
• Quenching-factoren: De studie bepaalt optimale GT quenching-factoren ($q_{opt}$) voor diverse Hamiltoniaanse om empirische matrixelementen te reproduceren. De resultaten benadrukken dat deze factoren fenomenologische kalibraties zijn die afhankelijk zijn van de specifieke interactie en valentieruimte, in plaats van universele constanten.

Betekenis en Claims
Het artikel beweert dat de voorgestelde Lanczos sterktefunctie-methode een praktische en efficiënte oplossing biedt voor het evalueren van $M_{2\nu}$ in modelruimtes waar directe sommatie computationeel te duur of onmogelijk is.

• Nauwkeurigheid: De methode behoudt de nauwkeurigheid van expliciete sommatie waar benchmarking mogelijk is, zelfs in gevallen met grote annuleringen en bijdragen van hoog-energetische toestanden.
• Nut voor $0\nu\beta\beta$: Gezien de sterke statistische correlatie tussen $2\nu\beta\beta$- en $0\nu\beta\beta$-NME's, kunnen betrouwbare en efficiënte $2\nu\beta\beta$-berekeningen helpen bij het identificeren van effectieve Hamiltoniaanse en operator-renormalisaties die gunstig zijn voor $0\nu\beta\beta$-voorspellingen.
• Schaalbaarheid: De efficiëntie van de methode maakt de exploratie van vergrote modelruimtes mogelijk (bijv. inclusief $2\hbar\omega$ excitaties) die voorheen onbereikbaar waren via directe diagonalisatie.
• Toekomstperspectief: De auteur suggereert dat dit raamwerk kan dienen als een complementaire grootschalige benchmark voor opkomende ab initio berekeningen, door statistische monsters en fenomenologische beperkingen te leveren die afgeleid zijn van gemeten halfwaardetijden.

Het werk stelt geen nieuwe experimentele signatures voor of claimt het $0\nu\beta\beta$-probleem direct op te lossen, maar positioneert de efficiënte berekening van $2\nu\beta\beta$ NME's als een noodzakelijk instrument om theoretische onzekerheden in de bredere zoektocht naar neutrinovrije verval te verminderen.