Heavy-ion collision simulation with high performance computer

Der Vortrag stellt eine vorläufige Arbeit vor, die Schwerionenkollisionssimulationen mit den Modellen DJBUU und SQMD unter Verwendung von Hochleistungsrechnern (HPC) durchführt, um die komplexen Eigenschaften dichter Kernmaterie besser zu verstehen.

Das Wesentliche

Schwere Ionen-Kollisionen am Supercomputer: Ein Blick ins Innere des Atomkerns

Stellen Sie sich vor, Sie möchten herausfinden, wie ein Auto funktioniert, aber Sie dürfen es nicht auseinanderbauen. Stattdessen lassen Sie zwei Autos mit voller Wucht frontal zusammenprallen und analysieren die Trümmer. Genau das machen Physiker mit Atomkernen, nur dass sie statt Autos winzige, schwere Atomkerne (wie Gold oder Blei) aufeinander schießen.

Dieser Bericht beschreibt, wie ein Team von Wissenschaftlern aus Korea und Kanada diese „Atom-Unfälle" am Computer simuliert, um Geheimnisse des Universums zu lüften.

1. Warum machen wir das? (Die Motivation)

Das Universum ist voller Rätsel. Neutronensterne sind wie gigantische, superdichte Kugeln aus Materie, die so schwer sind wie die Sonne, aber so klein wie eine Stadt. Wie verhält sich Materie unter solch extremem Druck?

• Der kosmische Blick: Astronomen haben kürzlich Wellen im Raum-Zeit-Gewebe (Gravitationswellen) von kollidierenden Neutronensternen gemessen. Das gibt uns Hinweise auf die Dichte dieser Sterne.
• Der irdische Blick: Da wir keine Neutronensterne in unser Labor bringen können, bauen wir auf der Erde kleine „Neutronensterne" nach. Wir schießen schwere Atomkerne gegeneinander, um für einen winzigen Moment extrem dichte Materie zu erzeugen.

In Korea wird dafür eine riesige Maschine namens RAON gebaut. Bevor man dort echte Experimente durchführt, muss man am Computer wissen, was passieren wird.

2. Die Werkzeuge: Zwei verschiedene „Spielregeln"

Um diese Kollisionen zu simulieren, nutzen die Forscher zwei verschiedene Computer-Modelle. Man kann sich das wie zwei verschiedene Regelsätze für ein Videospiel vorstellen, die versuchen, die Realität nachzubilden:

• DJBUU (Der „Fließende Fluss"):
  Dieses Modell betrachtet die Teilchen wie eine große Menge Wasser oder einen dichten Nebel. Es berechnet, wie sich die Wahrscheinlichkeit bewegt, dass ein Teilchen hier oder dort ist. Es ist sehr gut darin, das große Ganze zu sehen, wie sich der „Nebel" ausbreitet.

  • Besonderheit: Sie nutzen hier ein neues mathematisches Werkzeug (das QMC-Modell), das berücksichtigt, dass die Teilchen im Inneren aus noch kleineren Bausteinen (Quarks) bestehen, ähnlich wie ein Haus, dessen Ziegel sich unter Druck verformen.

• SQMD (Der „Ballon-Schwarm"):
  Dieses Modell sieht die Teilchen eher wie einzelne Bälle oder Ballons, die durch den Raum fliegen. Wenn sie sich zu nahe kommen, stoßen sie ab oder prallen voneinander ab. Es ist sehr gut darin, zu sehen, wie sich die einzelnen Bälle zu neuen Gruppen (Fragmenten) zusammenfinden.

3. Der Supercomputer-Held

Diese Berechnungen sind so komplex, dass ein normaler Laptop sie in einer Ewigkeit nicht schaffen würde. Man braucht einen Supercomputer (den „NURION" in Korea).
Stellen Sie sich vor, Sie müssten Millionen von Kugeln gleichzeitig durch ein Labyrinth steuern und dabei jede Kollision, jeden Stoß und jede Abstoßung in Echtzeit berechnen. Das erfordert enorme Rechenleistung. Der Supercomputer ist wie ein riesiges Team von Millionen von Gehirnen, die gleichzeitig arbeiten.

4. Was haben sie herausgefunden? (Die Ergebnisse)

Die Forscher haben ihre beiden Modelle gegeneinander getestet, indem sie Kollisionen simulierten, die bald am RAON stattfinden werden.

• Bei langsamen Kollisionen (50 AMeV): Beide Modelle liefern sehr ähnliche Ergebnisse. Es ist, als würden zwei verschiedene Wettervorhersagen für einen ruhigen Tag fast identische Temperaturen vorhersagen.
• Bei schnellen Kollisionen (100 AMeV): Hier beginnen die Modelle, sich zu unterscheiden. Je schneller die Kollision, desto mehr zeigen sich die Unterschiede in den „Spielregeln".
• Das Problem mit instabilen Kernen: Ein spannendes Ergebnis betraf instabile Atomkerne (wie Natrium-20). In einem Modell blieb der Kern stabil, im anderen zerfiel er fast schon, bevor er kollidierte. Das zeigt: Um instabile Teilchen korrekt zu simulieren, müssen die Modelle noch verfeinert werden.
• Die Dichte: Mit dem neuen QMC-Modell in DJBUU erreichten die simulierten Kollisionen eine noch höhere Dichte im Zentrum als mit dem alten Modell. Das ist wichtig, weil die Dichte bestimmt, wie sich Materie verhält – ähnlich wie bei einem Schwamm, der unter Druck noch mehr Wasser aufnehmen kann.

5. Fazit: Warum ist das wichtig?

Diese Arbeit ist wie das Testen von Flugzeugen im Windkanal, bevor sie gebaut werden. Durch diese Simulationen können die Wissenschaftler:

1. Bessere Vorhersagen für die Experimente am RAON treffen.
2. Die Eigenschaften von Materie verstehen, die es nur in Neutronensternen gibt.
3. Die Verbindung zwischen dem winzigen Atomkern und den größten Objekten im Universum (Neutronensternen) herstellen.

Kurz gesagt: Sie nutzen Supercomputer, um die härtesten Kollisionen im Universum im Labor nachzubauen, damit wir verstehen können, woraus das Universum wirklich besteht.

Technische Zusammenfassung

Technische Zusammenfassung: Schwerionen-Kollisionssimulationen mit Hochleistungsrechnern

1. Problemstellung und Motivation
Das Verständnis der Eigenschaften von dichter Kernmaterie ist ein zentrales Ziel sowohl der Astrophysik als auch der Kernphysik. Neuere Entdeckungen, wie die Gravitationswellen GW170817 (Verschmelzung von Neutronensternen) und X-Ray-Beobachtungen durch NICER, haben gezeigt, dass astrophysikalische Beobachtungen wertvolle Einblicke in den Zustand dichter Materie liefern. Parallel dazu bieten terrestrische Schwerionen-Kollisionsexperimente die Möglichkeit, diese Materiezustände unter kontrollierten Bedingungen zu untersuchen, insbesondere bei Dichten jenseits der normalen Kernmateriedichte.

Ein zentrales zukünftiges Experiment ist der Aufbau des RAON-Beschleunigerkomplexes (Rare isotope Accelerator complex for ON-line experiment) in Südkorea. Um die Ergebnisse dieser Experimente, insbesondere mit instabilen Isotopenstrahlen, vorherzusagen und zu interpretieren, sind hochkomplexe numerische Simulationen erforderlich. Diese Simulationen sind aufgrund der großen Anzahl an Teilchen, der kurzen Zeitschritte und der Notwendigkeit statistisch signifikanter Ereignismengen extrem rechenintensiv und erfordern Hochleistungsrechnersysteme (HPC).

2. Methodik
Die Autoren verwenden zwei verschiedene Transportmodelle, die auf dem HPC-System NURION des Korea Institute of Science and Technology Information (KISTI) implementiert wurden:

• DJBUU (DaeJeon Boltzmann-Uehling-Uhlenbeck):

  • Dies ist ein Modell vom BUU-Typ, das die zeitliche Entwicklung der Phasenraumdichte $f(\vec{x}, \vec{p})$ mittels der relativistischen BUU-Gleichung beschreibt.
  • Es nutzt die Testpartikel-Methode, bei der Nukleonen als Ensemble von Testpartikeln behandelt werden.
  • Ein charakteristisches Merkmal ist die Verwendung von Polynomen (statt Gauß- oder Dreiecksprofilen) für die Profile der Testpartikel im Orts- und Impulsraum, was eine glattere Form und endliche Integrationsgrenzen ermöglicht.
  • Die Dynamik wird durch Relativistische Mean-Field (RMF) Theorien gesteuert. In dieser Arbeit wurde der ursprüngliche QHD-Ansatz (Quantum Hadrodynamics) durch das Quark-Meson-Coupling (QMC) Modell erweitert, das zusätzliche Effekte der "Bag-Konstante" berücksichtigt und eine realistischere Kompressibilität ($K_0$) ohne starke Selbstwechselwirkungspotenziale ermöglicht.

• SQMD (Sindong Quantum Molecular Dynamics):

  • Dies ist ein QMD-Modell, das die Dynamik von $n$-Körper-Systemen behandelt, im Gegensatz zur 1-Körper-Dynamik bei BUU.
  • Nukleonen werden als Gaußsche Wellenpakete mit fester Breite modelliert. Die Phasenraumdichte wird über die Wigner-Transformation abgeleitet.
  • Die Nukleon-Nukleon-Wechselwirkung wird durch ein dichteabhängiges Skyrme-Potenzial beschrieben.
  • Für die Fragmentierung und Clusterbildung wird der Minimal Spanning Tree (MST) Algorithmus verwendet.

Beide Modelle simulieren die Propagation von Teilchen in Mean-Fields sowie Kollisionen (elastisch und inelastisch), wobei Pauli-Blocking-Faktoren berücksichtigt werden.

3. Wichtige Beiträge und Ergebnisse

• Vergleichende Studie (DJBUU vs. SQMD):

  • Es wurden Kollisionen von $^{208}\text{Pb} + ^{40,48}\text{Ca}$ bei Strahlenergien von 50 und 100 AMeV simuliert.
  • Bei niedrigen Energien (50 AMeV) zeigen beide Modelle konsistente Ergebnisse bezüglich der Größe der größten Fragmente (Biggest Fragments, BFs).
  • Bei höheren Energien (100 AMeV) und kleinen Stoßparametern treten signifikante Unterschiede auf, bedingt durch die unterschiedlichen Zustandsgleichungen und Modelldefinitionen.
  • Instabile Projektile: Bei Simulationen mit dem instabilen Isotop $^{20}\text{Na}$ ($^{20}\text{Na} + ^{208}\text{Pb}$) zeigten sich große Diskrepanzen: Das DJBUU-Modell erzeugte etwa 30 % größere Fragmente als SQMD.
  • Ursacheanalyse: Stabilitätstests (Propagation ohne Kollision) offenbarten, dass sich die Dichteprofile der instabilen Kerne in den beiden Modellen unterschiedlich entwickeln: DJBUU zeigt eine Verdichtung der Zentraldichte, während SQMD eine Ausdehnung zeigt. Dies unterstreicht die Sensitivität der Ergebnisse gegenüber der Initialisierung instabiler Kerne.

• Implementierung des QMC-Modells in DJBUU:

  • Der Vergleich zwischen dem ursprünglichen DJBUU+QHD und dem neuen DJBUU+QMC bei $^{40}\text{Ca} + ^{40}\text{Ca}$ Kollisionen zeigte, dass das QMC-Modell korrekt implementiert wurde (ähnliche allgemeine Ergebnisse).
  • Ein entscheidender Unterschied ist jedoch die höhere maximale Nukleondichte im Zentrum der Kollision beim QMC-Modell im Vergleich zum QHD-Modell.
  • Dies hat potenzielle Auswirkungen auf Observablen, die empfindlich auf die Zustandsgleichung reagieren, wie z. B. die Pion-Multiplizität oder den direkten Fluss.

4. Bedeutung und Ausblick
Die Arbeit demonstriert erfolgreich die Anwendung von HPC-Ressourcen für die Simulation von Schwerionen-Kollisionen mit den Modellen DJBUU und SQMD.

• Für RAON: Die Ergebnisse liefern wichtige Vorhersagen für Experimente mit seltenen Isotopenstrahlen am RAON-Facility. Die identifizierten Diskrepanzen zwischen den Modellen bei instabilen Kernen zeigen, dass zukünftige Verbesserungen (z. B. Implementierung von Oberflächen-termen in den Feldgleichungen und zeitabhängige Wellenpaketbreiten) notwendig sind, um die Stabilität simulierter Kerne zu erhöhen und Kollisionen mit instabilen Projektilen zuverlässiger zu beschreiben.
• Für die Physik dichter Materie: Die erfolgreiche Integration des QMC-Modells in Transportcodes eröffnet neue Wege, um die Zustandsgleichung dichter Kernmaterie zu untersuchen. Dies verbindet direkt die Ergebnisse von Schwerionenexperimenten mit astrophysikalischen Beobachtungen von Neutronensternen und Gravitationswellen.

Zusammenfassend liefert diese Studie einen wichtigen Beitrag zur Validierung und Weiterentwicklung von Transportmodellen für die nächste Generation von Schwerionenexperimenten und trägt zum besseren Verständnis der fundamentalen Eigenschaften der Kernmaterie bei.