Proton-Neutron Pairing in N=Z Nuclei within the Quark-Meson-Coupling Energy Density Functional

Deze studie toont aan dat proton-neutron pairingcorrelaties, berekend met het Quark-Meson-Coupling-energie-dichtheidsfunctionaal en het quartetcondensatiemodel, een aanzienlijke bijdrage leveren aan de bindingsenergieën van even-even N=Z-kernen en zo de overeenstemming met experimentele data verbeteren.

De kern

Deel 1: De Kern als een drukke dansvloer

Stel je een atoomkern voor als een enorme, drukke dansvloer. Op deze vloer dansen twee soorten gasten: protonen (positief geladen) en neutronen (neutraal). Normaal gesproken dansen de protonen met elkaar en de neutronen met elkaar, maar in speciale atomen – waar het aantal protonen precies gelijk is aan het aantal neutronen – kan er iets bijzonders gebeuren: een proton en een neutron kunnen samen een koppel vormen en als een paar dansen.

Deze wetenschappers (T. Popa, N. Sandulescu en D. Gambacurta) hebben gekeken naar wat er gebeurt als deze paren zich gedragen als een echte, hechte groep. Ze noemen dit "koppelpaardjes" (pairing correlations).

Deel 2: Twee manieren om te dansen

In de wereld van de atoomkernen zijn er twee manieren waarop deze paren kunnen dansen:

1. De "T = 1" dans (Isovector): Dit is de standaarddans. Het is een beetje zoals een schoolfeest waar iedereen in zijn eigen groepje danst, maar protonen en neutronen kunnen hier ook samen dansen.
2. De "T = 0" dans (Isoscalar): Dit is de "deuteron-dans". Hierbij vormen een proton en een neutron een heel speciale, sterke eenheid, alsof ze een mini-atoomkern (een deuteron) binnen de grote kern hebben gevormd.

De grote vraag in de wetenschap was: Kunnen deze twee dansvormen tegelijkertijd bestaan, of moet je kiezen? En: Hoe sterk is deze "deuteron-dans" eigenlijk?

Deel 3: De nieuwe dansvloer-regels (QMC)

Vroeger gebruikten wetenschappers modellen die de deeltjes als kleine, harde balletjes zagen. Deze nieuwe studie gebruikt een veel geavanceerder model genaamd QMC (Quark-Meson Coupling).

• De analogie: Stel je voor dat de oude modellen de dansers zagen als ondoorgrondelijke balletjes. Het QMC-model kijkt echter onder de huid van de balletjes. Het ziet dat elke danser (proton of neutron) eigenlijk bestaat uit drie nog kleinere deeltjes (quarks) die aan elkaar hangen met elastiekjes (de mesonen).
• Waarom is dit belangrijk? Omdat deze kleine deeltjes zich anders gedragen als ze in een drukke menigte zitten (de kern) dan als ze alleen zijn. Het QMC-model houdt rekening met deze veranderingen. Het is alsof je niet alleen kijkt naar hoe de mensen dansen, maar ook hoe hun kleding en houding veranderen door de drukte in de zaal.

Deel 4: De "Kwartet" dans (QCM)

Om te berekenen hoe deze dansers zich gedragen, gebruiken de auteurs een methode genaamd QCM (Quartet Condensation Model).

• De analogie: In plaats van alleen naar paren te kijken, kijken ze naar kwartetten. Een kwartet bestaat uit twee protonen en twee neutronen die samen een groepje vormen.
• Het probleem met oude methoden: De oude methoden (zoals HFB) waren als een dansmeester die de regels een beetje losjes hanteerde. Ze wisten niet precies hoeveel dansers er waren (het aantal deeltjes was niet exact) en ze hielden geen rekening met de draaiing van de hele groep.
• De oplossing: De QCM-methode is als een strenge, perfecte dansmeester. Hij telt elke danser exact, zorgt dat de groep perfect draait en houdt rekening met de isospin (een soort "dans-identiteit"). Hierdoor zijn de berekeningen veel nauwkeuriger.

Deel 5: Wat hebben ze ontdekt?

De onderzoekers hebben gekeken naar atoomkernen met massagetallen tussen 16 en 120 (van zuurstof tot tin). Hier zijn de belangrijkste bevindingen, vertaald naar alledaags taal:

1. De dans maakt de kern sterker: Als je rekening houdt met het samendansen van protonen en neutronen, wordt de atoomkern veel steviger (meer gebonden energie). Dit lost een oud probleem op: oude berekeningen voorspelden dat deze kernen iets te "los" waren. Met de nieuwe methode komen ze perfect overeen met de echte wereld.
2. Twee dansen gaan samen: In tegenstelling tot wat sommigen dachten, kunnen de "T=1" en "T=0" dansvormen naast elkaar bestaan. Ze vechten niet om de dansvloer, maar vullen elkaar aan.
3. Geen magische "super-dans" bij zware kernen: Er was een theorie dat bij zware kernen (boven tin) de "deuteron-dans" (T=0) de overhand zou nemen en de hele kern zou domineren. De nieuwe berekeningen zeggen: Nee, dat is niet zo. De T=0-dans is wel aanwezig en sterk, maar hij neemt niet de volledige controle over. De oude theorieën die dit voorspelden, hadden waarschijnlijk te weinig rekening gehouden met de precieze bewegingen van de deeltjes.
4. De "dubbel-magische" kernen: Er zijn speciale kernen (zoals 16O, 40Ca) die zo stabiel zijn dat ze als "dubbel magisch" worden beschouwd. In oude modellen was er geen dans (geen koppelvorming) mogelijk in deze kernen. Maar met de nieuwe, nauwkeurige methode zien ze dat er toch een beetje dans plaatsvindt, zelfs in deze stabiele kernen. Dit betekent dat we onze oude berekeningen voor deze kernen moeten herzien.

Conclusie

Deze studie is als het gebruik van een superkrachtige bril om naar de binnenkant van een atoomkern te kijken. Door te kijken naar de kleine quarks waaruit de deeltjes bestaan en door een zeer nauwkeurige methode te gebruiken om de "dans" van protonen en neutronen te beschrijven, hebben de onderzoekers laten zien dat:

• Het samendansen van protonen en neutronen essentieel is om de massa van atomen correct te begrijpen.
• De natuur complexer is dan gedacht: verschillende soorten "dansvormen" kunnen samen bestaan.
• We onze modellen moeten aanpassen, vooral voor de zwaarste en lichtste atomen.

Kortom: de atoomkern is geen statisch blokje, maar een dynamische, dansende menigte waar protonen en neutronen in perfecte harmonie samenwerken.

Technische samenvatting

Probleemstelling

De rol van proton-neutron (pn) pairingcorrelaties in kernen nabij de $N=Z$-lijn (waar het aantal neutronen gelijk is aan het aantal protonen) is een langdurig onderwerp van debat in de nucleaire fysica. Specifiek wordt onderzocht of isovector ($T=1$) en isoscalar ($T=0$) pairing bijdragen aan de extra bindingsenergie van $N=Z$-kernen, en of er een condensaat van isoscalair, deuteron-achtige paren kan bestaan.

Bestaande theoretische benaderingen, zoals de Hartree-Fock-Bogoliubov (HFB) methode, hebben beperkingen:

1. Ze behouden doorgaans niet exact het deeltjesaantal, impulsmoment en isospin.
2. HFB-berekeningen voorspellen vaak dat isovector pairing domineert en dat isoscalar pairing slechts in zware kernen ($A > 100$) dominant kan worden, maar deze voorspellingen zijn sterk afhankelijk van de kernvervorming en zijn nog niet getest in volledig zelfconsistente berekeningen.
3. Traditionele energiedichtheidsfunctionalen (EDF's) zoals Skyrme of Gogny behandelen nucleonen als puntdeeltjes en missen de expliciete quark-substructuur.

Het doel van dit werk is om de impact van zowel isovector als isoscalar pn-pairing te onderzoeken op de grondtoestandseigenschappen van even-even $N=Z$-kernen (massa $16 \le A \le 120$) binnen het kader van het Quark-Meson Coupling (QMC) model, gekoppeld aan het Quartet Condensation Model (QCM).

Methodologie

1. Theoretisch Kader:

• QMC Model (Mean Field): In plaats van nucleonen als puntdeeltjes te behandelen, beschrijft het QMC-model nucleonen als clusters van drie opgesloten quarks die interageren via scalaire ($\sigma$), vector ($\omega$) en isovector ($\rho$) mesonvelden. Dit leidt tot een microscopisch afgeleid energiedichtheidsfunctionaal met een natuurlijke dichtheidsafhankelijkheid en automatisch gegenereerde spin-baan- en tensorinteracties. De berekeningen gebruiken de nieuwste versie van de QMC-EDF (QMC$\pi$-III).
• QCM (Pairing): In plaats van de BCS- of HFB-approximatie, wordt het Quartet Condensation Model gebruikt. De grondtoestand van even-even $N=Z$-kernen wordt beschreven als een condensaat van quartetten (twee neutronen en twee protonen) die gekoppeld zijn tot totale isospin $T=0$.
  • Voordeel: De QCM behoudt exact het deeltjesaantal, impulsmoment en isospin, wat cruciaal is voor een realistische beschrijving van $N=Z$-systemen.
  • De grondtoestand wordt benaderd als $|QCM\rangle = [Q^\dagger_{T=1} + (\Delta^\dagger_0)^2]^{n_q} |0\rangle$, waarbij $n_q$ het aantal quartetten is.

2. Pairing Interacties:

• Isovector ($T=1$): Een nul-bereik interactie met een dichtheidsafhankelijke term die consistent is afgeleid binnen het QMC-kader. De sterkte wordt aangepast (via schalingsfactoren $s$) om de experimentele oneven-even massaverschillen (OEMD) te reproduceren.
• Isoscalar ($T=0$): Omdat een directe afleiding ontbreekt, wordt aangenomen dat de isoscalar interactie evenredig is met de isovector interactie: $V^{T=0}_P = w V^{T=1}_P$. De factor $w=1.6$ wordt gekozen op basis van eerdere HFB+QRPA studies van Gamow-Teller toestanden.

3. Rekenprocedure:

• Zelfconsistente berekeningen worden uitgevoerd voor axiaal vervormde kernen.
• Het proces is cyclisch: het QMC-middenveld wordt berekend $\rightarrow$ matrixelementen van de pairing-interactie worden bepaald $\rightarrow$ de QCM-vergelijkingen worden opgelost om bezettingswahrscheinlijkheden te vinden $\rightarrow$ deze worden gebruikt om het dichtheidsprofiel en het middenveld te updaten.
• De berekeningen omvatten kernen in de $sd$-schil ($16 < A < 40$), $pfg$-schil ($40 < A < 100$) en zware kernen boven $^{100}\text{Sn}$.

Belangrijkste Bijdragen

1. Eerste zelfconsistente QMC+QCM studie: Dit is de eerste toepassing van het QCM-framework binnen het QMC-EDF-kader, waarbij proton-neutron pairing expliciet wordt meegenomen voor een breed scala aan $N=Z$-kernen.
2. Behoud van symmetrieën: Door het gebruik van QCM worden de beperkingen van HFB (niet-behoud van deeltjesaantal en isospin) overwonnen, wat leidt tot een robuustere beschrijving van de correlatie-energieën.
3. Analyse van Isoscalar vs. Isovector Competitie: Het werk biedt inzicht in hoe isoscalar en isovector pairing met elkaar concurreren en coëxisteren, iets wat in standaard HFB-berekeningen vaak moeilijk te onderscheiden is.
4. Onderzoek naar Isoscalar Condensatie: Het test de voorspelling dat zware $N=Z$-kernen een overgang naar een dominant isoscalar pairing-fase ondergaan, een hypothese die eerder in HFB-studies werd geopperd.

Resultaten

• Bindingsenergieën:
  • Zonder pn-pairing (QMC+BCS) worden de bindingsenergieën van $N=Z$-kernen systematisch onderschat (residuen van enkele MeV).
  • Het inklappen van isovector pairing (QCM1) verbetert de overeenkomst met experimentele data aanzienlijk (verbetering van ca. 2 MeV voor $sd$-schil en 1-2 MeV voor $pfg$-schil).
  • Het toevoegen van isoscalar pairing (volledige QCM) leidt tot slechts een marginale extra verbetering (300-600 keV) voor de meeste kernen, hoewel er lokale effecten zijn (bijv. in $^{20}\text{Ne}$ en $^{64}\text{Ge}$).
• Coëxistentie van Pairing:
  • In tegenstelling tot sommige HFB-resultaten, vinden de auteurs dat isovector en isoscalar pairingcorrelaties coëxisteren in alle onderzochte kernen.
  • In sommige kernen (zoals $^{20}\text{Ne}$, $^{24}\text{Mg}$, $^{64}\text{Ge}$) is de isoscalar pn-pairing-energie groter dan de isovector pn-pairing-energie, maar nooit groter dan de totale isovector pairing-energie. Dit wordt toegeschreven aan de oscillatie van de niet-diagonale matrixelementen in de isoscalar kanaal.
• Vervorming en Shape Coexistence:
  • De berekeningen voorspellen mogelijke vorm-coëxistentie (het bestaan van meerdere lokale energie-minima) in kernen zoals $^{20}\text{Ne}$, $^{32}\text{S}$, $^{64}\text{Ge}$ en $^{68}\text{Se}$.
  • De isoscalar pairing heeft een beperkt effect op de vervormingsafhankelijkheid van de bindingsenergie, behalve in specifieke gevallen zoals $^{20}\text{Ne}$ waar het een tweede minimaal creëert.
• Zware Kernen ($A > 100$):
  • Voor kernen boven $^{100}\text{Sn}$ (zoals $^{108}\text{Xe}$) vinden de auteurs geen bewijs voor een duidelijk dominant isoscalar pairing-fase, in tegenstelling tot eerdere HFB-voorspellingen. De overgang is niet abrupt, wat suggereert dat eerdere resultaten een artefact waren van het gebrek aan exacte deeltjesaantalbehoud in HFB.
• Dubbel Magische Kernen:
  • In dubbel magische kernen ($^{16}\text{O}$, $^{40}\text{Ca}$, $^{56}\text{Ni}$, $^{100}\text{Sn}$) is de pairing-energie in BCS/HFB nul. De QCM-berekeningen tonen echter aan dat pairingcorrelaties (door het behoud van deeltjesaantal) een significante bijdrage leveren (enkele MeV) aan de totale bindingsenergie. Dit heeft implicaties voor het opnieuw fitten van EDF-parameters.

Betekenis en Conclusie

Dit onderzoek bevestigt dat proton-neutron pairing een essentiële rol speelt in het correct beschrijven van de bindingsenergieën van $N=Z$-kernen. De combinatie van het QMC-model (voor het middenveld met quark-substructuur) en het QCM (voor exact behoud van symmetrieën bij pairing) levert een nauwkeuriger en fysiek consistenter beeld op dan traditionele HFB-benaderingen.

De belangrijkste conclusies zijn:

1. Isovector pairing is de dominante bijdrage aan de extra binding van $N=Z$-kernen.
2. Isoscalar pairing is aanwezig en coëxisteert met isovector pairing, maar veroorzaakt geen globale fase-overgang naar een puur isoscalar condensaat in de onderzochte massa-gebieden.
3. De voorspelling van een dominante isoscalar fase in zware kernen door eerdere HFB-studies wordt niet bevestigd wanneer exacte deeltjesaantalbehoud en zelfconsistentie worden toegepast.
4. De bevindingen suggereren dat bij het gebruik van EDF's in combinatie met geavanceerde pairing-methoden (beyond-BCS), een herfiting van de parameters noodzakelijk is om dubbel telling van correlatie-energieën te voorkomen.

De auteurs plannen om dit kader uit te breiden naar kernen met een klein neutron-overschot om de overgang van $N=Z$ naar neutron-rijke systemen verder te onderzoeken.