A modern Fortran library for SU(3) coupling and recoupling coefficients

De auteurs hebben een moderne Fortran-bibliotheek ontwikkeld voor het berekenen van SU(3)-koppelings- en herschakelingscoëfficiënten die, in vergelijking met eerdere implementaties, een groter bereik aan kwantumgetallen ondersteunt en nauwkeurigere resultaten levert.

De kern

Stel je voor dat het universum een gigantische, ingewikkelde legpuzzel is. In de kernfysica proberen wetenschappers te begrijpen hoe atoomkernen (de centrale stukjes van die puzzel) in elkaar zitten. Om dit te doen, gebruiken ze een wiskundig gereedschapskistje genaamd SU(3).

Dit klinkt misschien als een saaie wiskundige code, maar je kunt het zien als een speciale taal om te beschrijven hoe de bouwstenen van de natuur (zoals quarks en protonen) met elkaar dansen en samenwerken.

Hier is wat de auteurs van dit artikel hebben gedaan, vertaald naar alledaags Nederlands:

1. Het Probleem: De Verouderde Rekenmachine

In de wereld van de kernfysica moeten wetenschappers vaak berekenen hoe twee groepen deeltjes samenkomen tot een nieuwe groep, of hoe ze van de ene dansvorm naar de andere veranderen. Hiervoor hebben ze specifieke getallen nodig: koppelingscoëfficiënten.

Voor jarenlang gebruikten wetenschappers een oude, beroemde rekenmachine (een computerprogramma genaamd de AD-bibliotheek) om deze getallen te berekenen.

• Het probleem: Deze oude machine was gebouwd met verouderde technologie (oud Fortran).
• De beperking: Als je probeerde om te rekenen met heel grote, complexe puzzelstukken (grote kwantumgetallen), begon de oude machine te haperen. De resultaten werden onnauwkeurig, of de machine gaf gewoon de verkeerde antwoorden. Het was alsof je probeerde een gigantisch gebouw te tekenen met een potlood dat steeds uitdroogde.

2. De Oplossing: De Nieuwe "Super-Rekenmachine"

De auteurs (Jakub Herko en zijn team) hebben een moderne, krachtige nieuwe bibliotheek gebouwd genaamd ndsu3lib.

Je kunt dit zien als het vervangen van die oude potloodrekenmachine door een moderne, digitale supercomputer met de volgende verbeteringen:

• Modernere taal: Ze hebben het geschreven in een moderne programmeertaal (Fortran 2003), die beter werkt met de snelle computers van vandaag.
• Preciezer: De nieuwe machine kan veel grotere en complexere puzzels oplossen zonder dat de resultaten "wazig" worden.
• Flexibeler: Hij kan kiezen tussen verschillende soorten "rekenkracht". Voor simpele taken gebruikt hij snelheid, maar voor de zwaarste taken schakelt hij automatisch over op "meerdere decimalen" (meerdere precisie), zodat er geen enkele foutje in de berekening komt.

3. Hoe werkt het? (De Analogie van de Dans)

Stel je voor dat je twee dansgroepen hebt die samenkomen om een nieuwe dans te maken.

• Koppelen: De oude methode kon alleen zeggen "ze dansen samen". Maar als er veel groepen waren, wist de oude methode niet precies welke danspasjes (de "buitenste veelvoudigheid") de juiste waren.
• De nieuwe aanpak: De nieuwe bibliotheek gebruikt een slimme truc (het "bouwproces" van Draayer en Akiyama). In plaats van te raden, bouwt hij de oplossing stap voor stap op, zoals een architect die eerst de fundering legt en dan laag voor laag de muren optrekt. Hierdoor weet hij precies welke danspasjes bij elkaar horen, zelfs bij heel grote groepen.

4. Waarom is dit belangrijk?

Dit klinkt misschien als puur wiskunde, maar het heeft grote gevolgen:

• Betere modellen: Wetenschappers kunnen nu atoomkernen modelleren die veel zwaarder en complexer zijn dan voorheen.
• Nieuwe inzichten: Het helpt bij het begrijpen van hoe sterren werken en hoe zware elementen in het universum ontstaan.
• Toekomstproof: Omdat de code modern is, kunnen andere wetenschappers hem makkelijk koppelen aan hun eigen software (zelfs als die in C++ of Python is geschreven).

Samenvatting

Kortom: De auteurs hebben een nieuwe, snellere en nauwkeurigere rekenmachine gebouwd voor de taal van de atoomkernen. Ze hebben de oude, haperende software vervangen door een modern systeem dat het universum beter kan "ontcijferen", zelfs als de puzzelstukken enorm groot en ingewikkeld zijn. Het is alsof ze de blauwdrukken voor de bouw van de materie hebben opgefrist van handgetekend papier naar een 3D-ontwerp op een krachtige computer.

Technische samenvatting

Titel: Een moderne Fortran-bibliotheek voor SU(3)-koppelings- en herkoppelingscoëfficiënten

Auteurs: Jakub Herko, Mark A. Caprio, Anna E. McCoy, Patrick J. Fasano.
Publicatie: The European Physical Journal A (voorlopig).

---

1. Het Probleem

De symmetriegroep SU(3) is van fundamenteel belang in diverse takken van de fysica, waaronder kernfysica, deeltjesfysica en kwantumoptica. Berekeningen binnen deze domeinen vereisen vaak de transformatie tussen gekoppelde en ongekoppelde producten van irreducibele representaties (irreps) van SU(3). Dit vereist de berekening van:

• Koppelingscoëfficiënten (Clebsch-Gordan coëfficiënten) voor de groepsketens $SU(3) \supset U(1) \times SU(2)$ (kanoniek) en $SU(3) \supset SO(3)$ (impulsmoment).
• Herkoppelingscoëfficiënten (analoog aan 6j- en 9j-symbolen) voor producten van drie of vier irreps die in verschillende volgorde zijn gekoppeld.

Bestaande software, met name de Akiyama-Draayer (AD) bibliotheek (oorspronkelijk geschreven in ouder Fortran), heeft ernstige beperkingen:

1. Numerieke onnauwkeurigheid: Bij grotere kwantumgetallen verliest de bibliotheek precisie en levert deze soms volledig onjuiste resultaten op (fouten van orde 1).
2. Beperkt bereik: Dit beperkt de toepasbaarheid in moderne kernstructuurberekeningen waar grote modelruimtes nodig zijn.
3. Oude architectuur: De code is niet geoptimaliseerd voor moderne computerarchitecturen en ondersteunt geen moderne functies zoals efficiënte multithreading of multiprecisie-aritmetiek op een robuuste manier.

2. Methodologie

De auteurs hebben ndsu3lib ontwikkeld, een moderne Fortran 2003-bibliotheek die de algoritmen van Draayer, Akiyama en Millener opnieuw implementeert. De kern van de methodologie omvat:

• Algoritme-implementatie: De bibliotheek implementeert het Draayer-Akiyama (DA) algoritme voor de berekening van gereduceerde koppelingscoëfficiënten (RCC's) en Millener's algoritmen voor herkoppelingscoëfficiënten.
• Biedenharn-Louck-Hecht (BLH) voorschrift: Om het probleem van de "externe multipliciteit" (waarbij meerdere lineair onafhankelijke irreps dezelfde symmetrie hebben) op te lossen, wordt strikt het BLH-voorschrift toegepast. Dit garandeert een unieke en reproduceerbare basis.
• Opbouw-proces (Building-up process): Het algoritme gebruikt een recursieve aanpak waarbij coëfficiënten voor een complexe koppeling worden afgeleid van eenvoudigere koppelingen met lagere multipliciteit, gebruikmakend van een "seed" (zaad) en recursierelaties gebaseerd op infinitesimale generatoren.
• Modernisering en Optimalisatie:
  • Gebruik van moderne Fortran 2003 voor betere leesbaarheid, onderhoudbaarheid en optimalisatie voor huidige CPU's.
  • Ondersteuning voor OpenMP voor multithreaded berekeningen.
  • Integratie van multiprecisie-aritmetiek (via de MPFUN2020 bibliotheek) om numerieke stabiliteit te behouden bij zeer grote kwantumgetallen.
  • Flexibele precisie: De bibliotheek schakelt automatisch over tussen dubbele (double), viervoudige (quadruple) en multiprecisie, afhankelijk van de grootte van de kwantumgetallen ($\lambda, \mu, L$), om zowel snelheid als nauwkeurigheid te optimaliseren.
  • Ondersteuning voor externe bibliotheken voor $SU(2)$-koppeling (GSL of WIGXJPF), waarbij WIGXJPF wordt aanbevolen voor hoge precisie.

3. Belangrijkste Bijdragen

• ndsu3lib: Een volledig nieuwe, open-source (MIT-licentie) Fortran-bibliotheek die zowel de kanonieke ($U(1) \times SU(2)$) als de impulsmoment ($SO(3)$) groepsketens ondersteunt.
• Verbeterde Nauwkeurigheid: De implementatie levert aanzienlijk nauwkeurigere resultaten dan de AD-bibliotheek, vooral voor $SU(3) \supset SO(3)$ koppelingscoëfficiënten met grote kwantumgetallen.
• Uitgebreid Bereik: De bibliotheek werkt correct voor een veel groter bereik van kwantumgetallen dan de AD-bibliotheek, waarbij de AD-bibliotheek faalt bij $\Sigma_w \geq 66$ (waarbij $\Sigma_w$ de som van de kwantumgetallen is).
• Documentatie en Afleiding: Het artikel biedt een zelfstandige afleiding van het DA-algoritme en documenteert de gebruikte formules, inclusief correcties op eerdere literatuur.
• Gebruiksgemak: Biedt wrappers voor C/C++ en is compatibel met moderne build-systemen (CMake).

4. Resultaten en Validatie

De auteurs hebben de bibliotheek uitgebreid getest en vergeleken met de AD-bibliotheek en een recente C++-implementatie (SU3lib):

• Validatie: De geldigheid van de berekende coëfficiënten is getest via de orthonormaliteitsrelaties en de vergelijkingen van de methode van infinitesimale generatoren.
• Precisie:
  • Voor $SU(3) \supset U(1) \times SU(2)$: ndsu3lib behoudt een foutenmarge van $\approx 10^{-15}$ tot $\Sigma_w = 81$, terwijl de AD-bibliotheek faalt bij $\Sigma_w \geq 66$.
  • Voor $SU(3) \supset SO(3)$: Met viervoudige of multiprecisie-aritmetiek behoudt ndsu3lib een hoge nauwkeurigheid ($< 10^{-12}$) over een breed bereik. De AD-bibliotheek toont exponentieel groeiende fouten.
  • Voor herkoppelingscoëfficiënten (U, Z, 9-(λ,µ)): De precisie van de drie bibliotheken is vergelijkbaar, maar ndsu3lib blijft stabiel waar de AD-bibliotheek onjuiste resultaten produceert.
• Snelheid:
  • De snelheid van ndsu3lib is vergelijkbaar met SU3lib en de AD-bibliotheek.
  • Bij grotere kwantumgetallen ($\Sigma_w > 20$) is ndsu3lib vaak sneller dan de AD-bibliotheek (tot een factor 3 sneller bij $\Sigma_w = 46$ voor $SO(3)$), voornamelijk omdat de AD-bibliotheek meer iteraties nodig heeft of faalt.
  • De gebruikte multiprecisie-aritmetiek heeft een acceptabel snelheidsverlies, maar is noodzakelijk voor de nauwkeurigheid.

5. Betekenis en Toekomstperspectief

De ndsu3lib bibliotheek is een cruciale tool voor de moderne theoretische kernfysica, met name voor symmetrie-gestuurde ab initio kernstructuurberekeningen (zoals het symplectische no-core configuration interaction framework).

• Het stelt onderzoekers in staat om grotere modelruimtes te behandelen met grotere nucleonenmassa's, wat essentieel is voor het begrijpen van zware kernen.
• Het biedt een robuustere en nauwkeurigere basis dan bestaande legacy-code, waardoor fouten door numerieke instabiliteit worden geëlimineerd.
• De open-source aard en de compatibiliteit met C/C++ en moderne compilers maken het breed toepasbaar in diverse wetenschappelijke projecten.

Kortom, dit artikel presenteert een noodzakelijke modernisering van een fundamenteel wiskundig hulpmiddel, waardoor de grenzen van wat berekenbaar is in de SU(3)-gebaseerde kernfysica worden verlegd.