Rendern von Spiralgalaxien

Bild 1: Spiralgalaxie M51. Die Galaxie besitzt eine deutlich sichtbare Spiralstruktur mit einer Vielzahl an H-II-Regionen die durch neu entstanden Sterne angeregt, in strahlendem Rot leuchten. Sie befindet sich in etwa 20-30 Millionen Lichtjahre Entfernung von der Erde.

Spiralgalaxien zählen zu den schönsten astronomischen Objekten am Nachthimmel. Das Bild rechts zeigt die Whirlpool-Galaxie, fotografiert vom NASA/ESA Weltraumteleskop Hubble. Es zeigt mit großem Detailreichtum wie eine Spiralgalaxie aufgebaut ist. Erkennbar wird eine deutliche, helle Spiralarmstruktur mit einer Vielzahl an rot leuchtenden H-II-Regionen. Das ganze wird von vielen feinen dunklen Streifen durchzogen, die von interstellaren Staubwolken gebildet werden. Diese absorbieren das Licht der dahinter liegenden Sterne und erscheinen dadurch dunkler. Sie liefern aber auch Gas und schwere Elemente, die benötigt werden um neue Sterne und Planetensysteme zu erzeugen, denn die Spiralarme sind auch der Ort, an dem neue Sterne entstehen.

So schön, wie es ist, es ist nur ein Bild. Wie würde es wohl aussehen, wenn die Galaxie rotiert? Das naheliegendste ist sicherlich, eine N-Körpersimulation durchzuführen, allerdings ist so etwas wesentlich komplizierter und rechenintensiver als hier wünschenswert wäre. Wir brauchen also eine Vereinfachung, etwas das uns die Berechnung der komplizierten Interaktionen zwischen tausenden von Sternen erspart und das trotzdem nahe an der Realität ist. Diese Vereinfachung wird uns die Dichtewellentheorie liefern.

Ungeduldige Leser können dem Bildlink unten folgen um direkt zu einer Implementierung des hier beschriebenen Algorithmus zu gelangen, für alle anderen erkläre ich erst einmal die Grundlagen.

Die Dichtewellentheorie

Bereits 1925 erkannte der schwedische Astronom Bertil Lindblad, dass unsere Milchstraße rotiert. Er erkannte weiterhin, dass die Winkelgeschwindigkeit dieser Rotation nach außen hin abnimmt. Damit war klar, dass die Spiralarme keine bloße statische Anhäufung von Sternen und Staub sein konnten, denn infolge der differentiellen Rotation der Milchstraße würden sich diese in kurzer Zeit auflösen (siehe "Winding Problem"; Animation 2). Als Lösung dieses Problems entwickelte er eine erste Version der Dichtewellentheorie. Grundlage ist die Annahme, dass im Gravitationsfeld einer Galaxie eine Dichtewelle mit konstanter Geschwindigkeit umläuft. Diese Winkelgeschwindigkeit stimmt im allgemeinen nicht mit der Rotationsgeschwindigkeit des restlichen Systems überein. Diese Ansätze wurden 1960 Jahren von den chinesischen Astronomen C.C. Lin und Frank Shu weiterentwickelt [1].

Die Dichtewellentheorie soll hier Ausgangspunkt für die Modellierung einer Galaxie werden. Sie stellt eine vereinfachtes Modell dar, das es ermöglicht die Dynamik einer Galaxie zu beschreiben ohne auf N-Körpersimulationen zurückgreifen zu müssen.

Klassifikation von Spiralgalaxien

Um die Vielzahl neu entdeckter Galaxien klassifizieren zu können, entwarf Edwin Hubble 1936 ein Ordnungsschema, das noch heute verwendet wird. Hubble unterschied dabei zwischen elliptischen- und Spiralgalaxien. Elliptische Galaxien werden mit den Bezeichnern E0 bis E7 versehen, wobei die Abplattung der Galaxie bis E7 immer mehr zunimmt. Der Typ S0 ist eine Übergangsform hin zu den Spiralgalaxien. Die Spiralgalaxien selbst werden nochmals in Spiralgalaxien (Sa bis Sc) und Spiralbalkengalaxien (SBa bis SBc) untergliedert.

Galaxientypen nach Edwin Hubble; Bild: Ville Koistinen (CC BY-SA 3.0)

Für die weitere Betrachtungen werden nur noch die Spiralgalaxietypen von Bedeutung sein. Diese unterscheiden sich wie folgt voneinander:

• Die Spiralarme werden von a nach c hin offener
• Der Anteil von Staub und jungen Sternen nimmt von a nach c hin zu.
• Der runde Zentralbereich der Galaxie wird im Verhältnis zur Größe der Galaxie von a nach c hin kleiner.

Bahnkurven

Bild 3: Gegeneinander um einen kleinen Winkel verschobene Elliptische Bahnkurven. Bild 4: Kernbereich und Scheibengröße für die in Bild 3 dargestellten Bahnkurven. Der Kernbereich gelb, die Grenze der Scheibe grün markiert.

Die erste Teilaufgabe der Modellierung besteht darin, die Form der Galaxie möglichst realistisch wiederzugeben. Die Dichtewellentheorie sagt im Prinzip voraus, dass die Spiralarme keine statischen Strukturen sind, sondern Gebiete größerer Sternendichte, die aufgrund der gegeneinander geneigten Bahnkurven der Sterne entstehen. Wir müssen also zunächst ein Modell für die Bahnkurven entwerfen.

Es liegt nahe, uns am Sonnensystem zu orientieren und einfach anzunehmen, das die Sterne sich auf Bahnen bewegen, die ähnlich dem ersten Keplerschen Gesetz sind. Die folgenden Eigenschaften scheinen für die Modellierung geeignet:

• Die Sterne bewegen sich auf elliptischen Bahnen um das galaktische Zentrum in deren Brennpunkt sich das galaktische Zentrum befindet (1. Keplersche Gesetz)
• Benachbarte Sterne bewegen sich in etwa in die gleiche Richtung, die Bahnkurven der Sterne kreuzen sich nicht.
• Mit zunehmendem Radius werden die Bahnkurven gegeneinander leicht gedreht, ähnlich einer Periheldrehung.

In der Realität sind die Bahnkurven keine Ellipsen. Auch die Annahme der Gültigkeit des ersten Keplerschen Gesetzes ist nicht ganz korrekt, denn die Masse einer Galaxie ist nicht punktförmig in deren Zentrum konzentriert, sondern flächig über die gesamte Scheibe verteilt. Es geht hier aber nicht um die möglichst exakte Berechnung von Bahnkurven, sondern nur um den optischen Eindruck und dafür reicht diese Näherung.

Heute weiß man, dass ein bedeutender Teil der Masse im Zentrum konzentriert ist. So hat vermutlich jede Spiralgalaxie ein superschweres Schwarzes Loch im Zentrum und dieses kann durchaus die Millionen- bis Milliardenfache Masse unserer Sonne haben. Zusätzlich ist im Kern die Sternendichte am größten, d.h. auf einem relativ kleinen Gebiet befindet sich ein Großteil der Gesamtmasse. Die Annahme einer Punktmasse im Zentrum ist also nicht ganz richtig, aber als grobe Näherung für unsere Zwecke akzeptabel. Bild 3 zeigt elliptische Bahnkurven, die mit zunehmendem Radius leicht gegeneinander verschoben sind. Es wird deutlich, dass bereits dieses einfache Modell der gegeneinander verschobenen Bahnkurven eine deutliche Spiralstruktur erzeugt. Die Spiralarme starten allerdings direkt im Zentrum und der Zentralbereich der Galaxie tritt nicht besonders deutlich hervor.

Beobachtungen zeigen, dass der Zentralbereich (engl.: Bulge) einer Galaxie häufig rund ist und teilweise ausgeprägte Balkenstrukturen besitzt. Daher sollte das Modell auch solche Formen erlauben. Bei den Balkengalaxien fangen die Spiralarme erst ab einem bestimmten Radius an. Es scheint also logisch einen Kernbereich zu definieren innerhalb dessen andere Bahnexzentrizitäten verwendet werden können. Für die Modellierung wird angenommen, dass der innerste Kern immer kreisförmig ist und die Exzentrizitäten bis zum Rand des Kernbereichs hin zunehmen. Diese Korrektur erlaubt ansatzweise Balkenstrukturen, ist allerdings nicht astronomisch begründet sind sondern rein optischer Natur (siehe Bild 4). Die Simulation von Balkenstrukturen ist eine Wissenschaft für sich und soll hier nicht Thema sein. Darüber hinaus wird aus optischen Gründen in der Modellierung ein Außenbereich definiert, der wiederum kreisrund ist und die komplette Spiralgalaxie einschließt. Je weiter man von einer Galaxie entfernt ist, um so besser kann deren Masse als Punktmasse genähert werden und es scheint zumindest nicht gänzlich unplausibel für Objekte, die eine Galaxie umkreisen nahezu kreisförmige Bahnen anzunehmen. Solche Objekte wären z.B. Kugelsternhaufen. (Die Bahnkurven des Außenbereichs wurden in Bild 4 weggelassen)

Flächenhelligkeitsprofile von Galaxien

Mit den Bahnkurven kann man nun daran gehen die Galaxie mit Sternen zu füllen. Dies sollte nach Möglichkeit so geschehen, daß die Sternverteilung die Oberflächenhelligkeit real beobachteter Galaxien wiedergibt. Im wesentlichen muss hier beachtet werden, dass die Flächenhelligkeit des Zentralbereichs anderen Gesetzmäßigkeiten unterliegt als die der galaktischen Scheibe. Die Flächenhelligkeit im Zentralbereich einer Spiralgalaxie folgt den selben Gesetzmäßigkeiten, wie die Flächenhelligkeit von elliptischen Galaxien. Sie kann mit dem de Vaucouleurs-Gesetz beschrieben werden:

wobei \(I_0\) die zentrale Intensität ist. Die Verteilung der Oberflächenhelligkeit einer galaktischen Scheibe folgt dabei nach Freeman [3] in guter Näherung einem Exponentialgesetz:

wobei \(I_0\) die Intensität an der Grenze zum Zentralbereich und \(R_D\) die Scheibenskalenlänge ist. (Länge auf der die Intensität um die Hälfte abfällt) Den Gesamthelligkeitsverlauf kann man durch Kombination der beiden Gleichungen ermitteln. Für die Modellierung wird angenommen, dass dieser Helligkeitsverlauf in direktem Zusammenhang mit der Sternhäufigkeit steht. Er dient daher als Grundlage für die Bildung einer Verteilungsfunktion für die Berechnung der Sternhäufigkeiten in Abhängigkeit des Galaxienradius. Die Sterne in der simulierten Galaxie werden so verteilt, dass deren Verteilungsfunktion der hier abgebildeten Kurve entspricht.

Bild 5.1: Simulierte Galaxie mit einer Sternenverteilung, die der theoretisch zu erwartenden Oberflächenhelligkeitsfunktion entspricht.

Bild 5.2: Schematische Darstellung des Oberflächenhelligkeitsverlaufes einer Spiralgalaxie. Das Diagramm kombiniert die von Freeman und de Vaucouleurs ermittelten Flächenhelligkeitsprofile für den Zentral- und Scheibenbereiche einer Galaxie.

Farben von Sternen

Bis jetzt wurde die Form der Galaxie und deren Sternverteilung betrachtet. Die Farbgebung der Galaxie wurde außer acht gelassen. Für eine möglichst realistische Wiedergabe der Galaxie ist die Farbgebung jedoch wesentlich, denn die Sterne können nur bestimmte, genau definierte Farben besitzen. Die Farbe eines Sternes hängt dabei primär von dessen Oberflächentemperatur ab. Das zugrunde liegende Modell ist die Schwarzkörperstrahlung.

Bild 6: Farbtemperatur nach dem planckschen Strahlungsgesetz.

Für eine möglichst realistische Farbgebung sollten also nur solche Farben verwendet werden, die physikalisch für Sterne möglich sind. Für diese Modellierung werden die Temperaturen der Sterne zufällig im Bereich zwischen 3000 und 9000 Kelvin verteilt. In der Realität durchläuft ein Stern während seiner Lebensspanne verschiedene Phasen, die sich auch in unterschiedlichen Oberflächentemperaturen äußern. Die Farbe eines Sterns ist also veränderlich und gibt unter anderem Hinweise auf dessen Alter.

Die Berechnung dieser Farben ist ein verhältnismäßig komplizierter Prozess, der hier nicht im Detail erläutert werden kann. Hier wird für die Berechnung der Farben auf das Program specrend.c von John Walker [5] zurückgegriffen, der sich auf seiner Webseite ausführlich mit diesem Thema auseinander gesetzt hat. Sein Programm ermöglicht die Umrechnung von Temperaturwerten in die Rot, Grün und Blauwerte die ein Computer für die Farbdarstellung benötigt. Eine C#-Version seines Algorithmus stelle ich auf GitHub zum download bereit.

Wendet man diese Farbgebung auf Bild 5 an, so erhält man die in Bild 7 abgebildete Galaxie. Es zeigt sich nun besser, dass die Sternhäufigkeit im Zentrum exponentiell zunimmt. Die Galaxie wirkt allerdings immer noch wie eine Ansammlung von Einzelpunkten. Die Darstellung besteht aus ca. 60000 einzelnen Punkten, Galaxien bestehen jedoch aus mehreren hundert Milliarden Sternen.

Einzelne Sterne in anderen Galaxien als der Milchstraße oder ihrer unmittelbaren Begleiter können in Fotografien von Galaxien aufgrund der Entfernung in der Regel gar nicht aufgelöst werden. Was man sieht ist eine Art Nebel, der sich erst beim starken hereinzoomen in einzelne Sterne auflösen würde. Aus diesem Grund benötigen wir in der Grafik eine weitere Komponente, die diesem Umstand Rechnung trägt. Zusätzlich zu den Sternen werden dafür tausende diffuser Nebelwölkchen verteilt. In OpenGL kann man dafür beispielsweise die Point sprite Extensions verwenden. Der Einfachheit halber werden diese Nebelelemente gleichmäßig über die galaktische Scheibe verteilt (siehe Bild 8).

Bild 7: Darstellung einer Galaxie mit Sternenfarben auf Basis des Planckschen Strahlungsgesetzes.

Bild 8: Galaxie mit Sternen und Nebelkomponente auf Basis von OpenGL Pointsprites. Die Spiralstrukturen treten nun deutlicher hervor,

Bild 9: Das Endergebnis zeigt eine Galaxie mit Sternen, Nebeln und H-II Regionen, die sich in den Spiralarmen konzentieren.

Mit Bild 8 wird die Gesamtstruktur schon gut erfaßt, aber es fehlt noch ein entscheidendes Merkmal. In Bild 1 sieht man deutlich eine Vielzahl roter Nebel, die sich in den Spiralarmen der Galaxie konzentrieren. Dabei handelt es sich um H-II Regionen. Das sind Sternentstehungsgebiete in ausgedehnten Wasserstoffwolken, die durch die darin entstehenden Sterne in einem typischen roten Licht erscheinen. Die hier modellierte Galaxie soll auch H-II Regionen enthalten, diese sollen sich dabei, wie in realen Galaxien, im wesentlichen in den Spiralarmen befinden. Für die Modellierung werden die H-II Regionen zunächst zufällig auf der galaktischen Scheibe verteilt. Um deren Sichtbarkeit auf die Spiralarme einzuschränken wird ein einfacher Trick angewendet. Die Größe mit der die H-II Region gezeichnet wird ist proportional zur Entfernung zweier benachbarter Dichtewellen an der Stelle, an der sie sich befindet.

Rotationskurve einer Spiralgalaxie

Bis jetzt wurde das Aussehen der Galaxie im statischen Zustand betrachtet. Im letzten Abschnitt geht es darum, die modellierte Galaxie wie ihre realen Vorbilder rotieren zu lassen.

Bild 10: Rotationskurve einer Galaxie. Bild: Phil Hibbs (CC BY-SA 3.0) Das Endergebnis der Modellierung zeigt eine Galaxie mit H-II Regionen, mit einem Rotationsverhalten, das real existierenden Galaxien nachempfunden ist. (Ingo Berg CC BY 4.0)

Dafür muss man die Rotationskurve von Galaxien näher betrachten. Rotationskurven zeigen, wie schnell die Sterne und das Gas einer Galaxie in Abhängigkeit von der Entfernung zum Zentrum rotiert. Beobachtungen zeigen, dass im inneren Bereich die Rotationskurve zunächst steil ansteigt. Dieses Verhalten entspricht in etwa der Rotation eines starren Körpers. Direkt im Anschluss daran flacht sich die Kurve ab und bleibt nahezu konstant (Bild 10; Kurve B).

Dies ist jedoch ein Widerspruch zu Vorhersagen auf Basis der Masseverteilung der sichtbaren Materie. Diese Vorhersagen gehen von einer nach außen hin abfallenden Rotationskurve aus (Bild 10; Kurve A). Der Unterschied zwischen theoretisch berechneten und real beobachteten Rotationskurven führte zur Hypothese der dunklen Materie, mit deren Masse die flach verlaufenden Kurve B erklärbar wäre. Für die Modellierung wird eine Flache Rotationskurve wie in Kurve B dargestellt verwendet. Das Endergebnis der Modellierung ist in in der folgenden Animation dargestellt. Es zeigt eine Galaxie mit zwei Spiralarmen und H-II Regionen, die in etwa einer Galaxie von Hubble-Typus Sb entspricht. Ihre Rotation entspricht der Rotation einer Galaxie unter dem Einfluß von dunkler Materie.

Download

Das Beispielprojekt ist in C++ geschrieben. Zum öffnen des Projektes wird die NetBeans IDE mit C++ Plugin oder Visual Studio benötigt.

Galaxienrenderer herunterladen (C++ Version)

Galaxienrenderer herunterladen (Typescript Version)

Quellenverzeichnis

1. C.C. Lin, Yuan, C., and F.H. Shu, "On the Spiral Structure of Disk Galaxies III. Comparison with Observations", Ap.J. 155, 721 (1969). (SCI)
2. Carroll, Bradley W. and Dale A. Ostlie, "An Introduction to Modern Astrophysics." (2007). Addison Wesley. pp. 966.
3. Freeman, K. C., "On the Disks of Spiral and S0 Galaxies" (1970). Astrophysical Journal, vol. 160, p.811
4. Density wave theory. (2010, Augst 03). In Wikipedia, the free encyclopedia. Retrieved August 03, 2010, from https://secure.wikimedia.org/wikipedia/en/wiki/Density_wave_theory
5. John Walker "Colour Rendering of Spectra." (1996, April 25). Webpage, retrieved from http://www.fourmilab.ch/documents/specrend/