Two-neutrino $ββ$ decay to excited states at next-to-leading order

De auteurs berekenen de kernmatrixelementen voor het dubbel-bètaverval naar aangeslagen toestanden in verschillende kernen met behulp van het schillenmodel en chiral effectieve veldtheorie, waarbij ze vaststellen dat de bijdragen van de volgende-orde-correxties over het algemeen klein zijn, maar significant kunnen worden door cancelaties, en dat de voorspelde halveringstijden consistent zijn met recente experimentele indicaties.

De kern

Titel: Het langzaamste dansje van het universum: Een verhaal over atoomkernen en dubbel-decay

Stel je voor dat atoomkernen niet statische blokken zijn, maar levende dansers op een podium. Soms willen deze dansers van outfit veranderen. In de meeste gevallen doen ze dit door een neutron om te zetten in een proton, waarbij ze een elektron en een antineutrinootje afschieten. Dit noemen we gewone radioactieve verval.

Maar er is een heel speciale, extreem zeldzame dans: de dubbel-beta-verval (of $\beta\beta$-verval). Hierbij probeert een atoomkern twee neutronen tegelijkertijd om te zetten in twee protonen. Het is alsof je twee dansers tegelijkertijd van het podium wilt laten springen.

Dit gebeurt op twee manieren:

1. De "spook" versie (Neutrinoloos): Ze gooien geen antineutrino's weg. Dit is nog nooit gezien en zou betekenen dat deeltjes hun eigen tegenpool zijn. Dit is de heilige graal van de fysica.
2. De "echte" versie (Twee-neutrino): Ze gooien wél twee antineutrino's weg. Dit gebeurt wel, maar het is zo traag dat het de langzaamste beweging is die we ooit hebben gemeten. Het kan meer dan $10^{21}$ jaar duren! Dat is miljarden malen ouder dan het heelal zelf.

Wat hebben deze onderzoekers gedaan?

De auteurs van dit paper (een team van de Universiteit van Barcelona) hebben zich gericht op de twee-neutrino versie, maar dan met een twist. Ze kijken niet alleen naar de dansers die op de grond blijven staan (de "grondtoestand"), maar ook naar diegenen die even een sprong maken naar een opgewekte toestand (een hogere energieladder) voordat ze landen.

Ze hebben dit voor verschillende atoomkernen gedaan (zoals Germanium, Selenium, Xenon), die gebruikt worden in grote experimenten om de "spook" versie te vinden.

Hoe hebben ze het berekend? (De Analogieën)

Om te voorspellen hoe snel deze dansjes plaatsvinden, moeten ze een ingewikkelde wiskundige formule oplossen. De onderzoekers hebben dit op een nieuwe, nog nauwkeurigere manier gedaan:

1. De "Bare" vs. "Kostuum" aanpak:
   Stel je voor dat je de kracht van de dansers meet. Soms meet je ze "naakt" (zonder extra effecten), en soms trek je ze een speciaal kostuum aan (een "renormaliseerde operator") dat rekening houdt met alle andere deeltjes in de kern die meedansen. De onderzoekers hebben beide methoden geprobeerd. Ze ontdekten dat de resultaten vaak lijken, maar dat het "kostuum" soms grote verschillen maakt als de dansers in de weg van elkaar lopen (wiskundige "annuleringen").

2. De "Nieuwe Stappen" (NLO):
   In de wiskunde van de deeltjesfysica werken ze met benaderingen. Ze beginnen met de belangrijkste stap (Leading Order) en voegen daarna kleinere, subtielere stappen toe (Next-to-Leading Order of NLO).

   • Vergelijking: Het is alsof je een routebeschrijving geeft. Eerst zeg je: "Ga rechtdoor." Dat is de basis. De NLO-stap is: "En oh ja, er staat een klein kuilje op de weg, pas je pas op."
   • De ontdekking: In de meeste gevallen is dat kuilje zo klein dat het de reistijd nauwelijks beïnvloedt (minder dan 5% verschil). Maar! Als de basisroute al heel erg onzeker is (door wiskundige annuleringen), kan dat kleine kuilje ineens heel groot worden en de hele voorspelling veranderen.

3. De Vorm van de Kern (Deformatie):
   Atomen zijn niet altijd perfecte ballen. Sommige zijn meer als een rugbybal (langwerpig) of zelfs als een ei met een knikje (triaxiaal).

   • De ontdekking: De onderzoekers ontdekten dat hoe meer de vorm van de startkern verschilt van de vorm van de eindkern, hoe moeilijker het is om die dubbele dans te maken. Het is alsof je probeert een dansstap uit te voeren terwijl je van een balletje op een rugbybal springt: het kost meer moeite en duurt langer. Ze keken zelfs naar hoe "zacht" of "stijf" deze vormen zijn.

Wat is het resultaat?

• De voorspellingen: Ze hebben berekend hoe lang het duurt voordat deze atoomkernen veranderen. Voor de meeste kernen is hun voorspelling nog steeds veel langer dan wat we nu in het lab kunnen meten. Het is dus nog niet waarschijnlijk dat we dit snel gaan zien.
• De uitzonderingen: Voor twee atomen (Germanium-76 en Selenium-82) zitten hun voorspellingen dichter bij de huidige meetgrenzen.
  • Bij Selenium-82 lijkt hun berekening zelfs te kloppen met een heel recente, nog niet bevestigde meting.
  • Bij Germanium-76 is het een beetje een raadsel. Sommige berekeningen zeggen "het is heel snel", andere zeggen "het is heel traag". Dit komt door de onzekerheid over de exacte "dansstijl" (de kernstructuur) van deze atomen.

Waarom is dit belangrijk?

Je vraagt je misschien af: "Waarom doen we dit als we de 'spook' versie (neutrinoloos) nog niet hebben gezien?"

Het antwoord is simpel: Om de 'spook' versie te vinden, moeten we eerst de 'echte' versie perfect begrijpen.
De twee versies zijn als tweelingbroers; ze delen dezelfde familiegeschiedenis (de kernstructuur). Als we de "echte" dans (twee-neutrino) niet goed kunnen voorspellen, kunnen we ook niet zeggen of een signaal dat we zien bij de "spook" versie echt nieuw is, of gewoon een foutje in onze berekening.

Conclusie

De onderzoekers hebben een zeer gedetailleerde kaart getekend van hoe atoomkernen zich gedragen tijdens dit zeldzame proces. Ze hebben laten zien dat de vorm van de atoomkern en de subtiele wiskundige correcties (de NLO-stappen) cruciaal zijn voor de nauwkeurigheid.

Hoewel ze nog niet de "spook" versie hebben gevonden, helpen ze de zoektocht door de "echte" versie beter te begrijpen. Het is alsof ze de achtergrondruis van een radio hebben verduidelijkt, zodat we straks misschien eindelijk het flauwe piepje van het nieuwe signaal kunnen horen.

Technische samenvatting

Probleemstelling

Het artikel richt zich op de theoretische berekening van de halfwaardetijden voor tweeneutrino-dubbel bèta-verval ($2\nu\beta\beta$) van de grondtoestand naar de eerste aangeslagen $0^+_2$-toestand in verschillende kernen die relevant zijn voor experimenten: $^{48}$Ca, $^{76}$Ge, $^{82}$Se, $^{124}$Sn, $^{130}$Te en $^{136}$Xe.

Hoewel $2\nu\beta\beta$-verval naar de grondtoestand ($0^+_{gs} \to 0^+_{gs}$) uitgebreid is bestudeerd, zijn berekeningen voor overgangen naar aangeslagen toestanden uitdagender. De kernproblemen zijn:

1. Overestimatie van NME's: Bestaande veeldeeltjesmethoden, zoals het nucleaire schillenmodel (NSM), de quasipartikel-random-fasebenadering (QRPA) en het interactieve bosonmodel (IBM), neigen om de nucleaire matrixelementen (NME's) voor $2\nu\beta\beta$ te overestimeren. Dit wordt vaak toegeschreven aan ontbrekende nucleaire correlaties en het ontbreken van twee-nucleon stromen.
2. Afhankelijkheid van Hamiltonian: De resultaten zijn sterk afhankelijk van de gekozen nucleaire Hamiltonian en de manier waarop de operator wordt behandeld (naakt vs. gereduceerd/renormaliseerd).
3. Nieuwe theoretische ontwikkelingen: Er is behoefte aan berekeningen die rekening houden met hogere-orde correcties uit de chirale effectieve veldtheorie ($\chi$EFT), specifiek de "next-to-leading order" (NLO) termen, die recent zijn geïntroduceerd.

Methodologie

De auteurs gebruiken het Nucleaire Schillenmodel (NSM) als basis voor hun berekeningen. De aanpak omvat de volgende stappen:

1. Hamiltonianen en Schillen: Voor elke kern worden specifieke interacties gebruikt die de spectroscopie van de initiële en finale kernen goed beschrijven:

   • $^{48}$Ca: KB3G en GXPF1A in de $pf$-schil.
   • $^{76}$Ge en $^{82}$Se: JJ4BB, GCN2850, JUN45 en RG in de $p_{3/2}, f_{5/2}, p_{1/2}, g_{9/2}$-orbitals.
   • $^{124}$Sn, $^{130}$Te, $^{136}$Xe: GCN5082 en QX in de $sdg$-schil.
   • Berekeningen zijn uitgevoerd met de codes ANTOINE en NATHAN.

2. Behandeling van Operatoren:

   • Er worden twee soorten Gamow-Teller (GT) operatoren gebruikt: de naakte GT-operator ($GT_{bare}$) en een effectieve GT-operator ($GT_{eff}$) die is gereduceerd binnen de valentie-schil.
   • Om de systematische overestimatie van het NSM te corrigeren, wordt een quenching-factor ($q_{\beta\beta}$) toegepast op de naakte operator. Deze factor wordt bepaald door de theoretische halfwaardetijden voor de grond-to-grond overgang af te stemmen op experimentele waarden. Voor de effectieve operator wordt soms geen quenching toegepast, maar wordt de impact van ontbrekende twee-nucleon stromen indirect meegenomen via de quenching die nodig is om grondtoestanden te reproduceren.

3. Lepton-uitbreiding en NLO-bijdragen:

   • De auteurs gebruiken een geavanceerd formalisme waarbij de energie van de leptonen wordt uitgebreid (Taylor-reeks), wat leidt tot hogere-orde termen in de NME's ($M^{(-3)}_{GT}$ en $M^{(-5)}_{GT}$).
   • Voor het eerst worden NLO-bijdragen uit $\chi$EFT meegenomen voor $2\nu\beta\beta$-verval. Dit omvat zwakke magnetisme en diagrammen met één-pion-uitwisseling. De berekening houdt rekening met drie termen in de expansie van de energiedeler.
   • De fase-ruimtefactoren (PSF) worden numeriek opgelost via de Dirac-vergelijking met de code RADIAL.

4. Structuuranalyse:

   • Er wordt een diepgaande analyse uitgevoerd van nucleaire vervorming (deformatie), inclusief triaxialiteit, met behulp van vorminvarianten ($Q_n$).
   • De relatie tussen het verschil in vervorming tussen initiële en finale toestanden en de grootte van de NME's wordt onderzocht.
   • De seniority-structuur van de toestanden wordt geanalyseerd om te begrijpen waarom bepaalde Hamiltonianen sterk afwijkende resultaten geven.

Belangrijkste Resultaten

1. Halfwaardetijden en Onzekerheden:

   • De berekende halfwaardetijden voor $0^+_{gs} \to 0^+_2$ overgangen variëren sterk afhankelijk van de gebruikte Hamiltonian. De onzekerheid wordt gedomineerd door de keuze van de nucleaire interactie, niet door de operatorkeuze (naakt vs. effectief).
   • In de meeste gevallen is de bijdrage van de NLO-termen aan de halfwaardetijd kleiner dan 5%.
   • Uitzonderingen: Bij kernen waar de leidende orde (LO) NME's sterke destructieve interferentie (cancellaties) vertonen (zoals $^{82}$Se en $^{136}$Xe), kunnen de NLO-bijdragen aanzienlijk toenemen (tot 10-30% of zelfs dominant worden).

2. Vergelijking met Experiment:

   • $^{76}$Ge: De voorspelde halfwaardetijden liggen dicht bij de huidige experimentele limiet, maar de theoretische spreiding is groot (factor 100, van $2 \cdot 10^{24}$ tot $180 \cdot 10^{24}$ jaar), afhankelijk van de Hamiltonian.
   • $^{82}$Se: De voorspellingen (bereik $0.4 - 190 \cdot 10^{23}$ jaar) zijn consistent met de zeer recente experimentele indicatie voor dit verval.
   • Andere kernen: Voor $^{48}$Ca, $^{124}$Sn, $^{130}$Te en $^{136}$Xe liggen de voorspelde halfwaardetijden minstens twee ordes van grootte boven de huidige experimentele ondergrenzen, wat detectie op korte termijn onwaarschijnlijk maakt.

3. Invloed van Nucleaire Structuur:

   • Er is een lineair verband gevonden tussen het verschil in vervorming ($\delta_{def}$) tussen initiële en finale kernen en de grootte van de NME's: kleinere verschillen in vorm leiden tot grotere NME's.
   • Triaxialiteit speelt een cruciale rol; het negeren van de $\gamma$-parameter leidt tot verkeerde correlaties.
   • De seniority-structuur is even belangrijk als de vervorming. Bijvoorbeeld, voor $^{82}$Se verklaart de afwezigheid van typische cancellaties tussen seniority-componenten in de JJ4BB-interactie de grote NME's, terwijl JUN45 wel cancellaties vertoont.

4. Quenching en Operatoren:

   • De resultaten met de naakte operator (gequenchd) en de effectieve operator (vaak ongequenchd) zijn over het algemeen vergelijkbaar, wat aangeeft dat de effectieve operator veel van de ontbrekende correlaties inherent bevat.

Significantie en Conclusie

Dit werk biedt de meest complete NSM-berekening tot nu toe voor $2\nu\beta\beta$-verval naar aangeslagen toestanden, inclusief de eerste toepassing van NLO-correcties uit $\chi$EFT.

• Theoretische Vooruitgang: Het artikel demonstreert dat hoewel NLO-termen in de meeste gevallen klein zijn, ze essentieel kunnen zijn wanneer LO-termen elkaar opheffen. Het benadrukt de noodzaak om nucleaire structuur (deformatie, triaxialiteit, seniority) nauwkeurig te modelleren om betrouwbare NME's te verkrijgen.
• Experimentele Implicatie: De resultaten zijn direct relevant voor experimenten zoals GERDA/LEGEND ($^{76}$Ge) en SNO+ ($^{130}$Te), en voor de interpretatie van recente signalen in $^{82}$Se. De grote onzekerheid in de voorspellingen voor $^{76}$Ge onderstreept de noodzaak van verdere theoretische verfijning en experimentele metingen.
• Koppeling aan Neutrinoloos Verval: Een beter begrip van $2\nu\beta\beta$-verval is cruciaal voor het beperken van onzekerheden in de voorspellingen voor neutrinoloos dubbel bèta-verval ($0\nu\beta\beta$), aangezien de NME's voor beide processen sterk gecorreleerd zijn.

De auteurs concluderen dat toekomstig werk gericht moet zijn op het opnemen van kort-bereik NLO-termen (waarvoor de koppelingsconstanten nog onbekend zijn) en twee-nucleon stromen, evenals op het verder verkennen van de relatie tussen nucleaire structuur en NME's.