MFGSB (ver. 1.0): Computer code for self-consistent mean-field calculations of atomic nuclei using Gaussian expansion method

Dit artikel kondigt de beschikbaarheid aan van MFGSB (versie 1.0), een computercode voor zelfconsistente middelveldberekeningen van atoomkernen met behulp van de Gauss-expansiemethode, die nu toegankelijk is via het repository van de Universiteit van Chiba.

De kern

Stel je voor dat je een gigantische, ingewikkelde puzzel probeert op te lossen. De stukjes van deze puzzel zijn de atoomkernen, de kleine hartjes van alle materie in het universum. De wetenschappers willen precies begrijpen hoe deze stukjes samenwerken, hoe ze bewegen en waarom ze stabiel zijn.

Het document dat je hierboven ziet, is eigenlijk een gebruikershandleiding en een uitnodiging voor een nieuwe, krachtige computercode genaamd MFGSB. Hier is wat dit allemaal betekent, vertaald naar gewoon Nederlands:

1. Wat is MFGSB eigenlijk?

MFGSB is een speciaal computerprogramma geschreven in een oude maar betrouwbare programmeertaal (FORTRAN77). Het is ontworpen om de "gemiddelde kracht" te berekenen die de deeltjes in een atoomkern op elkaar uitoefenen.

In plaats van te proberen elke deeltjes-botsing afzonderlijk te volgen (wat als het proberen is om elke druppel regen in een storm te tellen), kijkt het programma naar het gemiddelde gedrag. Het is alsof je niet elke auto op de snelweg individueel bestudeert, maar de stroom van het verkeer als geheel bekijkt om te zien hoe snel de file beweegt.

2. De "Magische" Methode: De Gaussische Expansie

Het programma gebruikt een slimme truc genaamd de "Gaussian expansion method" (GEM).

• De Analogie: Stel je voor dat je een schilderij wilt maken van een berg. Je kunt proberen elke steen en elke boom te tekenen, maar dat duurt eeuwen. Of je kunt de berg benaderen met een paar grote, zachte penseelstreken die de vorm perfect nabootsen.
• Het Voordeel: De code gebruikt deze "zachte penseelstreken" (wiskundige functies) om de vorm van de atoomkern te beschrijven. Het mooie is: deze penseelstreken werken voor elke berg (of atoomkern), of het nu een kleine steen is of een enorme rots. Je hoeft de instellingen niet telkens aan te passen als je van de ene kern naar de andere gaat. Dat bespaart enorm veel tijd en moeite!

3. Waarom is dit programma speciaal?

De auteur, een professor van de Universiteit van Chiba in Japan, noemt een paar dingen die dit programma uniek maken:

• Het kan alles: Het kan verschillende soorten krachten tussen de deeltjes simuleren, zelfs de heel lastige en complexe krachten die vaak worden genegeerd in andere programma's.
• Geen rommelige benaderingen: Veel andere programma's maken "ruwe schattingen" om het rekenwerk sneller te maken. MFGSB doet dit niet; het rekent de zware wiskundige details (zoals de Coulomb-kracht) precies uit, zonder te knoeien.
• Een universele sleutel: Omdat de basis-instellingen voor alle atoomkernen hetzelfde zijn, heb je maar één bestand nodig met de rekenregels om bijna alle atoomkernen in het heelal te bestuderen.

4. Wat kun je er precies mee doen?

Het programma is als een multifunctionele gereedschapskist. Je kunt er drie hoofdsoorten berekeningen mee doen:

1. HF (Hartree-Fock): De basisberekening van hoe de deeltjes zich gedragen.
2. HF+BCS: Voegt toe hoe de deeltjes paren vormen (een beetje zoals danspartners die samen bewegen).
3. HFB: De meest geavanceerde versie, die zowel de beweging als de paren perfect combineert.

Je kunt het programma ook instellen om te kijken naar kernen die perfect bol zijn, of juist langgerekt en vervormd (zoals een rugbybal).

5. Wat heb je nodig om het te gebruiken?

• Een sterke computer: Omdat de berekeningen zwaar zijn, heb je veel geheugen nodig. Voor de standaardinstelling heb je ongeveer 8 GB RAM nodig. Wil je nog preciezer zijn? Dan moet je het geheugen verdubbelen voor elke stap die je verder gaat.
• De software: Je kunt het programma gratis downloaden van de website van de Universiteit van Chiba.
  • Er is een groot pakket (ongeveer 16,6 GB) met alles erin.
  • Er is een licht pakket als je de zware data-bestanden apart wilt downloaden (handig als je internettraject beperkt is).

Samenvattend

Dit document kondigt een nieuw, krachtig gereedschap aan voor natuurkundigen. Het is als het krijgen van een nieuwe, super-snelle GPS voor het universum. In plaats van raden hoe atoomkernen werken, geeft MFGSB een nauwkeurige kaart, gebaseerd op een slimme methode die werkt voor bijna elk type kern zonder dat je de instellingen hoeft te veranderen.

Als je resultaten publiceert die met dit programma zijn berekend, moet je de auteurs natuurlijk even bedanken door te verwijzen naar hun eerdere werk. Maar voor nu: de code is klaar, staat online en wacht op de volgende ontdekking!

Technische samenvatting

Probleemstelling

De kernfysica staat voor de uitdaging om de structuur van atoomkernen nauwkeurig te beschrijven via zelfconsistente middelveldberekeningen (Self-Consistent Mean-Field, SCMF). Traditionele methoden hebben vaak te kampen met beperkingen, zoals:

• De noodzaak om basisfuncties specifiek af te stemmen (tunen) op elke specifieke nuclide of vervorming.
• Moeilijkheden bij het hanteren van complexe interacties, zoals de Yukawa-interactie inclusief tensor-kanaal, zonder extra benaderingen.
• Inefficiënte beschrijvingen van de asymptotisch gedrag van golf functies die energie-afhankelijk zijn.
• De complexiteit van het correct behandelen van de beweging van het zwaartepunt (center-of-mass, c.m.) en de Coulomb-uitwisselterm.

Methodologie

Het paper introduceert MFGSB, een computercode geschreven in FORTRAN77 die SCMF-berekeningen uitvoert voor atoomkernen met behulp van de Gaussian Expansion Method (GEM).

De kern van de methodologie omvat:

• Gaussian Expansion Method (GEM): In plaats van vaste roosters of standaard harmonische oscillator-basisfuncties, gebruikt de code een expansie in Gaussische functies. De parameters van deze basisfuncties zijn robuust en niet gevoelig voor specifieke nucliden of kernenvervormingen.
• Zelfconsistente Berekeningen: De code lost de Hartree-Fock (HF), Hartree-Fock-BCS (HF+BCS) en Hartree-Fock-Bogolyubov (HFB) vergelijkingen op.
• Symmetrieën: Berekeningen kunnen worden uitgevoerd onder verschillende symmetrievoorwaarden:
  • Sferische symmetrie (J) met pariteit (P) en tijdomkering (T).
  • Axiale symmetrie (rotatie rond de z-as) met P, T en reflectiesymmetrie rond de y-as.
  • Alleen $J_z$ symmetrie met P en $RT$ symmetrie.
• Numerieke Oplossing: Er zijn twee methoden beschikbaar om zelfconsistentie te bereiken:
  1. Iteratieve methode: De Hamiltoniaan wordt bij elke iteratie diagonaliseerd.
  2. Conjugate Gradient methode: De Hamiltoniaan wordt pas aan het einde diagonaliseerd (efficiënter voor bepaalde gevallen).
• Interacties: De code ondersteunt diverse effectieve twee-nucleon interacties, met name de Yukawa-interactie tot aan het tensor-kanaal. Ook de Coulomb-uitwisselterm en de twee-lichaams-term voor de zwaartepuntsbeweging worden exact behandeld zonder extra benaderingen.

Belangrijkste Bijdragen

De belangrijkste bijdragen van MFGSB (versie 1.0) zijn:

1. Unieke Capabiliteit: Het is de enige SCMF-code die de Yukawa-interactie volledig kan hanteren, inclusief het tensor-kanaal, zonder extra benaderingen.
2. Universele Basisparameters: Een unieke eigenschap is dat de parameters van de basisfuncties onafhankelijk zijn van het specifieke nuclide of de vervorming. Hierdoor is één set parameters toepasbaar op een groot deel van het nucleaire kaart (nuclear chart), wat de noodzaak tot parameter-tuning elimineert.
3. Efficiënte Asymptotiek: De code beschrijft het energie-afhankelijke asymptotische gedrag van de enkel-deeltje (s.p.) of quasipartikel (q.p.) golf functies efficiënt binnen een redelijke precisie.
4. Flexibiliteit: De code biedt opties voor niet-zelfconsistente berekeningen (Woods-Saxon of harmonische oscillator potentiaal) naast de volledige SCMF-modi (HF, HF+BCS, HFB).
5. Beschikbaarheid: De volledige code (inclusief interactiegegevens) is openbaar beschikbaar via de Chiba University Repository, hoewel de interactiebestanden (die het grootste deel van de 16,6 GB beslaan) optioneel apart kunnen worden gedownload.

Resultaten en Specificaties

Hoewel het paper voornamelijk een technische handleiding en beschrijving van de code is, worden de volgende technische specificaties en resultaten van de implementatie gepresenteerd:

• Systeemvereisten: Voor de standaardwaarde $\ell_{cut} = 7$ wordt ongeveer 8 GB RAM geheugen vereist. Bij elke verhoging van $\ell_{cut}$ met 1 verdubbelt het geheugengebruik ongeveer.
• Afhankelijkheden: De code vereist de installatie van de BLAS en LAPACK bibliotheken.
• Bestandsomvang: De volledige distributie (mfgsb.zip) is 16,6 GB groot, voornamelijk door de interactiebestanden. Een "light" versie is beschikbaar waarbij de interactiebestanden (int_data_l7) apart worden gedownload.
• Berekeningsmogelijkheden: De code kan kwadrupoolmomenten constraineerbaar maken (voor axiale en $J_z$ symmetrie) en biedt keuzevrijheid in de behandeling van de zwaartepuntsbeweging.

Significantie

MFGSB is een significante toevoeging aan het instrumentarium van de theoretische kernfysica om de volgende redenen:

• Verhoogde Nauwkeurigheid: Door de exacte behandeling van de Coulomb-uitwisseling en de zwaartepuntsbeweging, evenals de mogelijkheid om tensor-krachten te gebruiken, kunnen fysisch realistischere modellen worden opgesteld.
• Efficiëntie en Schaalbaarheid: Het vermogen om één set basisparameters te gebruiken voor een breed scala aan kernen maakt het mogelijk om systematische studies uit te voeren over het hele periodiek systeem zonder tijdrovende aanpassingen voor elke nieuwe kern.
• Open Science: Door de code en de data openbaar beschikbaar te stellen (via DOI), faciliteert het onderzoekers wereldwijd om zelfconsistente middelveldberekeningen uit te voeren met geavanceerde interacties, wat de reproduceerbaarheid en voortgang in het veld bevordert.

Het paper dient primair als een technische documentatie en gebruikershandleiding voor wetenschappers die deze geavanceerde berekeningen willen uitvoeren, met verwijzingen naar eerdere werken [1-4] waarin de numerieke methoden zijn ontwikkeld.