Forschungsbericht 2017 - Max-Planck-Institut für Astrophysik

Die Kugelsternhaufen-Simulation DRAGON: eine Million Sterne, Schwarze Löcher und Gravitationswellen

Autoren
Naab, Thorsten; Spurzem, Reiner; Wang, Long für die DRAGON collaboration
Abteilungen
Forschungsgruppe „Computergestützte Strukturbildung“
Zusammenfassung
Ein Team von Experten aus Europa und China hat die ersten Kugelsternhaufen-Simulationen mit einer Million Sternen auf dem High-Performance GPU-Cluster der Max Planck Computing and Data Facility durchgeführt. Diese – bis jetzt – größten und realistischsten Simulationen reproduzieren die beobachteten Eigenschaften von Kugelsternhaufen in bisher unerreichter Genauigkeit und erlauben einen Blick in die dunkle Welt der Schwarzen Löcher. Sie prognostizieren zentrale Ansammlungen von einzelnen und binären Schwarzen Löchern.

Kugelsternhaufen sind in der Tat rätselhafte Objekte. Sie bestehen aus Hunderttausenden leuchtenden Sternen sowie Sternüberresten in einer Region von nur wenigen Dutzend Parsec (bis zu 100 Lichtjahre) – sie sind die dichtesten und ältesten gravitativ gebundenen Sternsysteme im Universum. Ihre Sterndichte im Zentrum kann eine Million Mal größer sein als die Sterndichte in der Nähe unserer Sonne. Etwa 150 Kugelsternhaufen umkreisen die Milchstraße; massereiche Galaxien können sogar über 10.000 gravitativ gebundene Kugelsternhaufen besitzen. Da sich ihre Sterne meist zur gleichen Zeit gebildet haben, allerdings mit unterschiedlichen Massen, sind Kugelsternhaufen ideale Laboratorien um Stellardynamik und Sternentwicklung zu untersuchen.

Dynamik und Sternentwicklung in Kugelsternhaufen

Die dynamische Entwicklung der Kugelsternhaufen ist sehr komplex. Anders als in Galaxien ist die Sterndichte so hoch, dass Sterne im nahen Vorbeiflug gravitativ wechselwirken oder sogar physisch miteinander kollidieren. Aufgrund dieser Wechselwirkungen gibt es mehr fest gebundene Doppelsterne als bei normalen galaktischen Feldsternen. Darüber hinaus sinken massereiche Sterne durch den Prozess der Massensegregation ins Zentrum des Systems.

Die Evolution eines Kugelsternhaufens als Ganzes wird durch den Lebenszyklus von Einzel- und Doppelsternen verkompliziert. Bei jungen massereichen Sternen (mit mehr als 8 Sonnenmassen) kommt es zu einem erheblichen Massenverlust in ihrer Sternwindphase; am Ende ihres Lebens explodieren diese Sterne als Kernkollaps-Supernova. Die Überreste dieser längst erloschenen Sterne sind Neutronensterne oder Schwarze Löcher. Letztere können Massen im Bereich von 10 bis 50 Sonnenmassen besitzen. Sie sind für normale elektromagnetische Beobachtungen unsichtbar und konnten bis vor Kurzem nur indirekt nachgewiesen werden.

Das Licht der Kugelsternhaufen wird von wenigen hundert sehr stark leuchtenden roten Riesensternen dominiert. Die meisten anderen Sterne im System haben eine wesentlich geringere Masse als unsere Sonne und eine viel geringe Leuchtkraft. Aus diesem Grund ist das Hubble Weltraumteleskop ein bevorzugtes Instrument um die Sternpopulation von Kugelsternhaufen zu beobachten. Farben-Helligkeits-Diagramme (CMD) mit Daten von Hubble weisen im Vergleich zu erdgebundenen Instrumenten aufgrund der sehr niedrigen photometrischen Fehler (dadurch ergeben sich scharfe Strukturen wie die Hauptsequenz, der Riesen- oder der Weißer-Zwerg-Zweig) und der sehr hohen Empfindlichkeit eine bessere Qualität auf. So konnten mit Hubble zum ersten Mal Eigenschaften von weißen Zwergen mit geringer Leuchtkraft und die Hauptsequenz von Sternen geringer Masse in hoher Qualität beobachtet werden.

Hochpräzise Computersimulationen

original
Original 1508158463
Abb. 1: Oben: Der Hydra-Supercomputer (1,7 Petaflop/s) betrieben von der Max Planck Computing and Data Facility ist ausgestattet mit 676 Kepler K20 GPGPU-Beschleunigern (1 Petaflop/s, unten links). Dieser Supercomputer wurde verwendet, um die DRAGON-Simulationen durchzuführen. Unten rechts: Der Laoh-Supercomputer der Nationalen Astronomischen Observatorien, Chinesische Akademie der Wissenschaften in Peking (96 Teraflop/s) betrieben am Zentrum für Information und Computing ist mit 64 Kepler K20 GPGPU-Beschleunigern ausgestattet.
Abb. 1: Oben: Der Hydra-Supercomputer (1,7 Petaflop/s) betrieben von der Max Planck Computing and Data Facility ist ausgestattet mit 676 Kepler K20 GPGPU-Beschleunigern (1 Petaflop/s, unten links). Dieser Supercomputer wurde verwendet, um die DRAGON-Simulationen durchzuführen. Unten rechts: Der Laoh-Supercomputer der Nationalen Astronomischen Observatorien, Chinesische Akademie der Wissenschaften in Peking (96 Teraflop/s) betrieben am Zentrum für Information und Computing ist mit 64 Kepler K20 GPGPU-Beschleunigern ausgestattet.

Seit vielen Jahren stellt es eine große Herausforderung dar, die Entwicklung eines massereichen Kugelsternhaufens mit selbstkonsistenten numerischen Simulationen nachzuverfolgen. Zum ersten Mal gelang es nun einem Team internationaler Experten am Max-Planck-Institut für Astrophysik (MPA), der Chinesischen Akademie der Wissenschaften und der Peking-Universität die – bisher – realistischsten Simulationen der Evolution eines Kugelsternhaufen mit einer Million Sterne im Gezeitenfeld unserer Milchstraße über einen Zeitraum von rund 12 Milliarden Jahren zu simulieren. Die Simulationen, die auf dem Hydra-Supercomputer der Max Planck Computing and Data Facility (MPCDF) (Abb. 1) im Rahmen des internationalen DRAGON-Projektes durchgeführt wurden, setzen einen neuen Standard bei der Modellierung von Kugelsternhaufen.

Emöglicht wurden diese Simulationen durch deutliche Verbesserungen der Software auf dem Laohu-Supercomputer des Center of Information and Computing at National Astronomical Observatories der Chinesischen Akademie der Wissenschaften. Der Code kann gleichzeitig Mehrfachknoten-Parallelisierung, OpenMP auf den Knoten und die Kepler K20-Grafikkarten-Parallelisierung (GPGPUs) nutzen, um die Gravitationskräfte zwischen den Sternen zu berechnen. Eine typische DRAGON-Sternhaufensimulation verwendet 8 Hydra-Knoten mit 160 CPU-Kernen und mehr als 32k GPU-Threads für eine Gesamtrechenzeit von rund einem Jahr (8.000 Stunden).

original
Original 1508158464
Abb. 2: Künstliches Farbbild (BVI) aller Sterne eines simulierten Kugelsternhaufens (zentrales Bild, etwa 60 pc im Durchmesser), wie sie sich nach 12 Milliarden Jahren entwickelt haben.
Abb. 2: Künstliches Farbbild (BVI) aller Sterne eines simulierten Kugelsternhaufens (zentrales Bild, etwa 60 pc im Durchmesser), wie sie sich nach 12 Milliarden Jahren entwickelt haben.

Die Entwicklung der Sternpopulation eines Kugelsternhaufens kann nun im Detail nachverfolgt werden, über alle dynamischen und Sternentwicklungsphasen hinweg, einschließlich dem Verlust von Sternen im Gezeitenfeld der Milchstraße. Die Hauptentwicklungsphasen von Einzel- und Doppelsternen mit ganz unterschiedlichen Massen (0,08–100 Sonnenmassen) werden ebenfalls simuliert (Abb. 2). Die DRAGON-Simulationen können auch dazu verwendet werden, synthetische Farben-Helligkeits-Diagramme (CMD) wie mit Hubble (Abb. 3) zu erstellen.

Verschmelzende Schwarze Löcher und Graviationswellen

In der DRAGON-Simulation bilden die schwarzen Löcher – 10 bis 50 Sonnenmassen schwere Überreste massereicher Sterne – im Zentrum des Systems einen dichten Kern (Abb. 2, auf weißem Hintergrund). In der klassischen Astronomie können dieses Schwarzen Löcher nur indirekt durch den Einfluss ihrer Gravitation auf leuchtende Sterne beobachtet werden. Einige Dutzend Schwarze Löcher bilden Doppelsysteme und verlieren durch Gravitationsstrahlung Energie; dieser Prozess ist in unseren Simulationen enthalten.

original
Original 1508158464
Abb. 3: Vergleich des Farben-Helligkeits-Diagramms des mit dem Hubble Space Telescope beobachteten Sternhaufens NGC4372 mit denen von zwei simulierten Haufen. Zur Simulation der Beobachtungen wurde eine typische Entfernung zu einem galaktischen Kugelsternhaufen angenommen sowie die Spezifikation der Kameras an Bord des Hubble-Weltraumteleskops mit CACAO.
Abb. 3: Vergleich des Farben-Helligkeits-Diagramms des mit dem Hubble Space Telescope beobachteten Sternhaufens NGC4372 mit denen von zwei simulierten Haufen. Zur Simulation der Beobachtungen wurde eine typische Entfernung zu einem galaktischen Kugelsternhaufen angenommen sowie die Spezifikation der Kameras an Bord des Hubble-Weltraumteleskops mit CACAO.

Vor Kurzem konnte die LIGO-Kollaboration Gravitationswellenemission aus dem Verschmelzen zweier Schwarzer Löcher (mit Massen von 36 und 29 Sonnenmassen) in einem Abstand von 410 Mpc nachweisen. Unsere DRAGON-Kugelsternhaufen erzeugen solche Verschmelzungen mit ähnlichen Parametern; etwa zehn Ereignisse in jedem Haufen. Daher erwarten wir, dass in den kommenden Monaten oder Jahren weitere derartige Ereignisse beobachtet werden. An einer genaueren Vorhersage für Gravitationswellen-Ereignisse in unseren Modellen wird gearbeitet. Dies hängt nicht nur von der Entwicklung der Haufen ab, sondern auch von der Anzahl und der Verteilung der Kugelsternhaufen im Universum. Schon jetzt können wir aber vorhersagen, dass Kugelsternhaufen – ähnlich unseren DRAGON-Haufen – eine mögliche Quelle der kürzlich beobachteten spektakulären Gravitationswellen-Ereignisses sind.

Die jetzt nachgewiesene Verschmelzung zweier Schwarzer Löcher ist wahrscheinlich nur die Spitze des Eisbergs. Die dynamische Entwicklung der zentralen Bereiche der simulierten Haufen wird durch Hunderte (wenn nicht Tausende) von einzelnen Schwarzen Löchern und solchen in Doppelsternsystemen beherrscht. Zukünftige Studien werden zeigen, ob solche Ansammlungen von stellaren Schwarzen Löchern in den Zentren der meisten Kugelsternhaufen existieren.

original
Original 1508158466
Abb. 4: Kumulative Massenverteilung der einzelnen Komponenten aus Abb. 2. Schwarze Löcher populieren das Zentrum des Haufens (schwarze Linie). Hauptreihensterne geringer Masse (cyanfarbene Linie) dominieren die Gesamtmasse des Systems und ihre Verteilung ist ausgedehnter. Die Punkte zeigen den Halbmassenradius der jeweiligen Population.
Abb. 4: Kumulative Massenverteilung der einzelnen Komponenten aus Abb. 2. Schwarze Löcher populieren das Zentrum des Haufens (schwarze Linie). Hauptreihensterne geringer Masse (cyanfarbene Linie) dominieren die Gesamtmasse des Systems und ihre Verteilung ist ausgedehnter. Die Punkte zeigen den Halbmassenradius der jeweiligen Population.

Das DRAGON-Projekt

Das DRAGON-Projekt ist eine Supercomputing-Initiative von NAOC/CAS und Kiaa/PKU (Peking), MPA (Garching), CAMK (Warschau) und IoA ( Cambridge), um die Entwicklung von Kugelsternhaufen mit GPU-gestützten High-Performance-Simulationen zu untersuchen.

Das Team besteht aus Rainer Spurzem (unterstützt vom Thousand-Talents-Programm der Volksrepublik China an den Nationalen Astronomischen Observatorien von China (NAOC), Peking und der Universität Heidelberg), Long Wang und Thijs (M.B.N.) Kouwenhoven (Kavli-Institut für Astronomie und Astrophysik, Universität Peking), Peter Berczik (NAOC und Hauptsternwarte der nationalen Akademie der Wissenschaften der Ukraine in Kiew), Sverre Aarseth (Institut für Astronomie, Cambridge), Mirek Giersz und Abbas Askar (Nicolaus Copernicus Astronomical Center, Warschau), Thorsten Naab und Riko Schadow (Max-Planck-Institut für Astrophysik, Garching) sowie weiteren Studenten und Mitarbeitern. Die Berechnungen werden an der MPCDF durchgeführt, Software-Entwicklung und vorläufige Simulationen am KIAA und am Laohu-Cluster von NAOC.

Literaturhinweise

1.
Wang, L.; Spurzem, R.; Aarseth, S.; Nitadori, K.; Berczik, P.; Kouwenhoven, M. B.N.; Naab, T.
NBODY6++GPU: ready for the gravitational million-body problem
Monthly Notices Royal Astronomy Society 450, 4070-4080 (2015)
2.
Wang, L.; Spurzem, R.; Aarseth, S.; Giersz, M.; Askar, A.; Berczik, P.; Naab, T.; Kouwenhoven, M. B. N.; Schadow, R.
The DRAGON simulations: globular cluster evolution with a million stars
Monthly Notices Royal Astronomy Society 458, 1450-1465 (2016)
Zur Redakteursansicht