Heavy-ion collision simulation with high performance computer

Dit artikel presenteert voorlopige resultaten van zware-ionenbotsingssimulaties met behulp van de DJBUU- en SQMD-modellen op high-performance computers om de complexe eigenschappen van dichte kernmaterie te bestuderen.

De kern

Stel je voor dat je een enorme, onzichtbare Lego-constructie hebt: een atoomkern. Nu, in het heelal, gebeuren er soms dingen waarbij twee van deze enorme constructies met een snelheid van bijna het licht op elkaar botsen. Dit noemen we zware-ionenbotsingen.

Wetenschappers doen dit niet zomaar voor de lol. Ze willen weten hoe materie zich gedraagt als het extreem dicht en heet is. Denk aan de binnenkant van een neutronenster (een ster die zo zwaar is dat een theelepel ervan zo zwaar is als een berg) of de oerkracht die direct na de Big Bang heerste.

Hier is wat dit paper vertelt, vertaald naar alledaags taalgebruik:

1. Het Grote Doel: De "Super-Dichte" Materie

De wetenschappers (uit Korea en Canada) willen begrijpen hoe deze super-dichte materie werkt. Ze hebben twee manieren om dit te doen:

• Kijken naar de sterren: Door zwaartekrachtsgolven van botsende neutronensterren te meten.
• Het nabootsen op Aarde: Door in deeltjesversnellers (zoals de RAON in Zuid-Korea) atoomkernen met elkaar te laten botsen.

Het probleem? De wiskunde achter deze botsingen is zo ingewikkeld, dat je een supercomputer nodig hebt om het uit te rekenen. Het is alsof je probeert te voorspellen wat er gebeurt als je twee zakken vol met honderden knikkers met een hamer op elkaar slaat, maar dan in 3D en met quantum-wiskunde.

2. De Twee Simulatoren: DJBUU en SQMD

Om deze botsingen te simuleren, hebben de auteurs twee verschillende "virtuele laboratoria" gebouwd. Je kunt ze zien als twee verschillende soorten videospelletjes die proberen hetzelfde scenario te spelen, maar met een andere motor:

• DJBUU (De "Stroom" Simulator):
  Deze kijkt naar de botsing als een stromende vloeistof. Het behandelt de deeltjes (nucleonen) als een grote menigte die zich verplaatst volgens bepaalde regels. Het is alsof je kijkt naar een drukke menigte op een plein: je ziet niet elke individuele persoon, maar de stroming van de massa.

  • Het nieuwe trucje: Ze hebben een nieuwe regel toegevoegd (het QMC-model). In plaats van de deeltjes als simpele balletjes te zien, kijken ze nu alsof de deeltjes uit nog kleinere stukjes (quarks) bestaan die in een zakje zitten. Dit maakt de simulatie iets realistischer voor extreme situaties.

• SQMD (De "Balletjes" Simulator):
  Deze kijkt naar de botsing als een verzameling individuele balletjes die tegen elkaar aan stuiteren. Het is alsof je een potje biljart speelt, maar dan met duizenden ballen die ook nog eens een beetje "wazig" zijn (kwantummechanica).

  • Het verschil: Deze simulator is beter in het zien van hoe de balletjes in groepjes (fragmenten) bij elkaar blijven plakken na de botsing.

3. De Supercomputer: De "NURION"

Omdat deze berekeningen zo zwaar zijn, gebruiken ze een supercomputer genaamd NURION (van het Koreaanse instituut KISTI).

• Vergelijking: Als een gewone laptop een simpele taak doet (zoals een e-mail sturen), doet de NURION het alsof hij duizenden mensen tegelijkertijd een complexe puzzel laat oplossen.
• Ze draaien duizenden botsingen na, met miljarden kleine stappen in de tijd, om statistisch zeker te zijn van hun resultaten.

4. Wat Vonden Ze? (De Resultaten)

Toen ze de twee simulatoren lieten vechten tegen elkaar (met verschillende atoomkernen), vonden ze interessante dingen:

• Bij lage snelheid: Als de botsing niet te hard is, geven beide simulatoren ongeveer hetzelfde resultaat. Het is alsof twee verschillende voorspellers zeggen: "Het gaat regenen."

• Bij hoge snelheid: Als de botsing heel heftig is, beginnen de resultaten te verschillen. De ene simulator zegt: "De brokken worden groot," de andere: "Nee, ze worden kleiner."

  • Waarom? Omdat de onstabiele atoomkernen (zoals een instabiel isotoop van Natrium) in de ene simulator al uit elkaar beginnen te vallen voordat ze botsen, en in de andere nog heel blijven. Het is alsof je twee verschillende modellen van een instabiel torenblok bouwt: bij de ene valt hij al om als je er een beetje aan wrijft, bij de andere staat hij stevig.

• De Nieuwe Regels (QMC): Door de nieuwe "quark-regels" in de DJBUU-simulator te stoppen, zagen ze dat de materie in het midden van de botsing nog dichter werd dan voorheen. Dit is belangrijk, omdat de dichtheid bepaalt hoeveel nieuwe deeltjes (zoals pions) er ontstaan.

5. Waarom is dit belangrijk?

Dit onderzoek is als het bouwen van een betere "tijdmachine".

1. Voor de toekomst: Het helpt om de experimenten voor te bereiden die binnenkort in Korea (RAON) gaan plaatsvinden met nieuwe, zeldzame atoomkernen.
2. Voor het heelal: Het helpt ons te begrijpen wat er gebeurt in neutronensterren. Als we weten hoe materie zich gedraagt in een zware-ionenbotsing op Aarde, kunnen we beter begrijpen waarom neutronensterren niet instorten tot een zwart gat.

Kort samengevat:
De auteurs hebben twee verschillende computerspellen gemaakt om atoomkernen te laten botsen. Ze hebben deze op een supersterke computer laten draaien om te zien hoe de materie zich gedraagt. Ze ontdekten dat bij extreme snelheden de twee spellen verschillende uitkomsten geven, wat hen vertelt dat ze hun regels (de wiskunde) nog moeten verfijnen. Dit helpt ons uiteindelijk om de geheimen van de zwaarste objecten in het universum te ontrafelen.

Technische samenvatting

Technische Samenvatting: Zware-ionenbotsingssimulaties met HPC

Probleemstelling
Het begrijpen van de eigenschappen van dichte kernmaterie is cruciaal voor zowel de astrofysica (bijv. neutronensterren, geïmpliceerd door de GW170817-gravitationele golven en NICER-observaties) als de kernfysica. Terwijl terrestrische experimenten met zware ionen (zoals die gepland zijn bij de RAON-faciliteit in Zuid-Korea) deze eigenschappen kunnen testen, zijn de onderliggende systemen extreem complex. Het simuleren van de volledige tijdsontwikkeling van zware-ionenbotsingen vereist enorme rekenkracht vanwege het grote aantal deeltjes, de korte tijdstappen en de noodzaak om statistisch significante aantallen gebeurtenissen te accumuleren. Bestaande theoretische modellen moeten worden gekalibreerd en vergeleken om betrouwbare voorspellingen te doen voor experimenten met zeldzame isotoopstralen.

Methodologie
De auteurs hebben twee transportmodellen ontwikkeld en geïmplementeerd op het High Performance Computing (HPC)-systeem NURION van het Korea Institute of Science and Technology Information (KISTI):

1. DJBUU (DaeJeon Boltzmann-Uehling-Uhlenbeck):

   • Een semi-klassiek model dat de relativistische BUU-vergelijking oplost.
   • Gebruikt de "test particle"-methode waarbij nucleonen worden behandeld als een ensemble van testdeeltjes.
   • Een uniek kenmerk is het gebruik van polynoomprofielen (in plaats van Gaussians of driehoeken) voor de fase-ruimtedichtheid, wat zorgt voor een soepelere vorm met eindige eindpunten.
   • De interacties worden beschreven via een Relativistische Middenveld (RMF) theorie. In deze studie wordt een nieuwe variant geïntroduceerd: DJBUU+QMC, die het Quark-Meson Coupling (QMC) model gebruikt in plaats van het traditionele Quantum Hadrodynamics (QHD) model. Het QMC-model houdt rekening met de Bag-constante en berekent het effectieve nucleonmassa op quark-niveau.

2. SQMD (Sindong Quantum Molecular Dynamics):

   • Een QMD-type model dat de dynamica van $N$-lichamen simuleert (in tegenstelling tot de 1-lichaamsdistributie in BUU).
   • Nucleonen worden behandeld als Gaussische golfpakketten met een vaste breedte.
   • Gebruikt een dichtheidsafhankelijk Skyrme-potentieel voor nucleaire interacties.
   • Voor fragmentatie en clustering wordt het "Minimal Spanning Tree"-algoritme toegepast.

De simulaties omvatten botsingen zoals $^{208}\text{Pb} + ^{40,48}\text{Ca}$ en $^{20}\text{Na} + ^{208}\text{Pb}$ bij straalenergieën van 50 en 100 AMeV.

Belangrijkste Bijdragen

• HPC-Implementatie: Succesvolle uitvoering van zware-ionenbotsingssimulaties op het NURION-systeem, wat nodig is om de enorme rekenlast van duizenden deeltjes en miljoenen roostercellen te verwerken.
• Modelontwikkeling: Implementatie van het geavanceerde QMC-model binnen het DJBUU-framework, wat een brug slaat tussen quark-niveau fysica en macroscopische botsingsdynamica.
• Vergelijkende Analyse: Een grondige vergelijking tussen de twee fundamenteel verschillende transportbenaderingen (BUU vs. QMD) onder dezelfde experimentele condities.

Resultaten

• Vergelijking DJBUU vs. SQMD:

  • Bij lagere straalenergieën (50 AMeV) leveren beide modellen vergelijkbare resultaten op voor de grootte van de grootste fragmenten (BFs).
  • Bij hogere energieën (100 AMeV) en kleinere impactparameters ontstaan grotere verschillen, veroorzaakt door de verschillende toestandsvergelijkingen (EoS) en modeldefinities.
  • Bij het gebruik van een onstabiel projectiel ($^{20}\text{Na}$) vertonen de modellen aanzienlijke discrepanties: DJBUU produceert fragmenten die ongeveer 30% groter zijn dan die van SQMD.
  • Oorzaak van discrepantie: Analyse toont aan dat de onstabiele $^{20}\text{Na}$-kernen zich verschillend ontwikkelen in de modellen. DJBUU toont een toenemende centrale dichtheid (compacte kern), terwijl SQMD een afnemende en verspreidende dichtheid laat zien. Dit suggereert dat de initiële stabiliteit en vervorming van onstabiele kernen in de modellen nog verfijning behoeven.

• QMC vs. QHD in DJBUU:

  • De implementatie van het QMC-model (DJBUU+QMC) levert resultaten op die qua algemene evolutie vergelijkbaar zijn met het QHD-model (DJBUU+QHD).
  • Een kritiek verschil is dat DJBUU+QMC een hogere maximale nucleon-dichtheid bereikt in het centrum van de botsing.
  • Deze hogere dichtheid is significant omdat observabelen zoals pionmultipliciteit en directe stroming (direct flow) zeer gevoelig zijn voor de toestandsvergelijking van dichte kernmaterie.

Betekenis en Toekomstperspectief
Deze studie biedt een essentiële basis voor het interpreteren van toekomstige experimenten bij de RAON-faciliteit, met name die met stralen van zeldzame isotopen. De gevonden verschillen tussen de modellen benadrukken de noodzaak van verdere verfijning, zoals het implementeren van tijdsafhankelijke golfpakketbreedtes en oppervlaktermen in de veldvergelijkingen om de stabiliteit van onstabiele kernen beter te modelleren.

Door de koppeling tussen zware-ionenbotsingen, neutronensterren en gravitationele golven, dragen deze simulaties bij aan een beter begrip van de fundamentele eigenschappen van dichte kernmaterie. De resultaten wijzen erop dat het QMC-model een waardevol alternatief biedt dat mogelijk nieuwe inzichten biedt in de toestandsvergelijking bij extreme dichtheden.