Keine geschlossene Formel bedeutet nicht absolutes Chaos
In Liu Cixins Science-Fiction-Klassiker Die drei Sonnen und der breiten öffentlichen Wahrnehmung gilt das Dreikörperproblem als Inbegriff absoluter Unvorhersehbarkeit. Doch dahinter verbirgt sich ein über hundert Jahre altes Missverständnis. In den 1880er Jahren wies Henri Poincaré nach, dass für die Bewegung dreier Himmelskörper keine allgemeine geschlossene analytische Lösung (keine algebraischen ersten Integrale) existiert – eine Erkenntnis, die den Grundstein für die moderne Chaostheorie legte. Doch das Fehlen einer universellen Formel bedeutet keineswegs, dass in chaotischen dynamischen Systemen keine exakten periodischen Lösungen existieren.
Die Webplattform Three Body Orbits hat 3.915 exakte periodische Lösungen aus 13 wissenschaftlichen Publikationen in einem interaktiven Atlas zusammengetragen. Drei Himmelskörper identischer Masse umkreisen einander unter präzise kalibrierten Anfangsbedingungen und kehren nach Ablauf einer vollen Periode exakt an ihre Ausgangspositionen und zu ihren Ausgangsgeschwindigkeiten zurück. Dass die theoretische Mathematik über keine geschlossene Formel verfügt, ist etwas grundlegend anderes als Unvorhersehbarkeit in der numerischen Ingenieurpraxis. Steht ausreichend Rechenleistung zur Verfügung, lassen sich selbst im Ozean des Chaos Inseln mathematischer Ordnung kartieren.
Abbildung: Gesamtansicht der Startseite von Three Body Orbits. Jede periodische Bahn nimmt einen Punkt auf der Himmelskarte ein; topologisch verwandte Trajektorien bilden Inselstrukturen. Quelle: Screenshot von Three Body Orbits
An den Grenzen von 32-Stellen-Arithmetik: Warum auf Echtzeitsimulation im Browser verzichtet wurde
Hinter den 3.915 im Atlas erfassten Bahnen steht eine immense ingenieurtechnische Rekonstruktionsleistung. Der Entwickler übernahm keineswegs blind die in wissenschaftlichen Fachaufsätzen abgedruckten Anfangsbedingungen. Aus typografischen Gründen veröffentlichen Paper ihre Parameter typischerweise nur auf 10 bis 15 signifikante Stellen gerundet. In einem chaotischen System wie dem Dreikörperproblem genügt jedoch bereits ein minimaler Rundungsfehler in der Größenordnung von 10⁻¹⁵, der sich innerhalb weniger Umläufe exponentiell aufschaukelt und das System zur sofortigen Divergenz und Auflösung bringt.
Um jede Bahn nahtlos und geschlossen von Anfang bis Ende zu verbinden, rekonvergierte der Autor die Poincaré-Schnittabbildungen (Poincaré return maps) der zugrunde liegenden Differentialgleichungen mithilfe von Arbitrary-Precision-Bibliotheken und dem Newton-Raphson-Verfahren. Durch iterative Kalibrierung der Anfangspositionen und Geschwindigkeiten auf 32 Dezimalstellen konnte der Median des Schließungsfehlers auf beispiellose 1,4 × 10⁻²⁴ gedrückt werden. In dieser Größenordnung ist die Positionsabweichung der Himmelskörper nach einem vollen Umlauf weitaus kleiner als der Durchmesser eines einzelnen Atoms. 3.889 der Bahnen weisen eine Schließungsgenauigkeit von besser als 10⁻⁹ auf.
Diese unter extremem rechnerischen Aufwand generierten Trajektorien stellten das Frontend-Engineering vor ein Dilemma. Hätte man sich für eine Echtzeitsimulation im Browser entschieden – bei der JavaScript lokal mit nativen 64-Bit-Gleitkommazahlen mit doppelter Genauigkeit integriert –, wären die geschlossenen Bahnen durch akkumulierte Rundungsfehler innerhalb von Sekunden auseinandergedriftet.
Daher wählte die Website den Architekturansatz “Vorberechnete Splines + Wiedergabe im Client”. Das Backend komprimiert die hochpräzisen Integrationskoordinaten in über 9,7 Millionen Spline-Knoten höherer Ordnung und überträgt diese an den Client. Das Frontend muss diese Knotenpunkte lediglich anhand von Zeitstempeln interpolieren und flüssig rendern. Indem Vorberechnungs-Speicherplatz gegen Laufzeitgenauigkeit getauscht und die rechenintensive Last auf das Backend verlagert wurde, kann der Browser seine gesamte Leistung für eine ruckelfreie Darstellung mit 60 Bildern pro Sekunde nutzen.
Abbildung: Detailansicht der berühmten Achterbahn (Figure-8), bei der drei Körper auf verschlungenen Pfaden umeinander kreisen. Quelle: Screenshot von Three Body Orbits
Drei Forschergenerationen rechnen 3.915 exakte Lösungen heraus
In diesen 3.915 präzisen Bahnen spiegelt sich ein jahrzehntelanger wissenschaftlicher Kraftakt der nichtlinearen Dynamik wider. Es handelt sich nicht um das Werk eines einzelnen Instituts, sondern um eine systematische Klassifikationskarte, die von internationalen Forschungsteams über Generationen hinweg auf Supercomputern erarbeitet wurde.
Das Team um Shijun Liao an der Jiaotong-Universität Shanghai (Li–Liao-Gruppe) steuerte mit 1.964 Bahnen den Löwenanteil bei. Mithilfe von Supercomputer-Clustern und ihrer Methode der sauberen numerischen Simulation (Clean Numerical Simulation, CNS) spürten sie systematisch eine Fülle zuvor unbekannter topologischer Strukturen auf. Ein Forschungsteam der Universität Sofia in Bulgarien steuerte 709 komplexe Bahnen bei und bereicherte die Familien höherer Harmonischer erheblich. Auf die in der Mathematik berühmten Choreographien (Choreographies) entfallen 583 Bahnen – Konfigurationen, bei denen die drei Körper wie Artisten auf einem Seil tanzen und dieselbe geschlossene Raumkurve teilen.
Hinzu kommen 311 Speziallösungen mit Freifall-Charakteristik (Free-fall) sowie 230 geometrisch eindrucksvolle Lösungen der Belgrader Forschungsgruppe in Serbien (Šuvakov und Dmitrašinović). Sogar 52 klassische Bahnen, die Roger Broucke bereits 1975 herleitete, sind vertreten. Zusammen mit 66 durch Crowdsourcing auf der Website neu entdeckten Lösungen formen sie ein lückenloses dynamisches Kompendium. Hinter jedem Koordinatensatz stehen Hunderttausende auf Rechenclustern verbrannte CPU-Stunden.
Ein Tritt per Schieberegler: Die Floquet-Stabilitätsgrenze im Stresstest
1.593 Bahnen im Atlas sind explizit als “linear stabil” klassifiziert – ein Anteil von 40,7 %. In der Theorie dynamischer Systeme ist der Floquet-Multiplikator die maßgebliche Kennzahl, um zu beurteilen, wie robust eine periodische Bahn gegenüber Störungen ist. Er misst, wie sich winzige Abweichungen über eine volle Periode hinweg fortpflanzen.
Liegt der Betrag eines Floquet-Multiplikators über 1, wächst der Fehler exponentiell an und das System bricht rasch zusammen. Beträgt der maximale Floquet-Multiplikator hingegen exakt 1,000, befindet sich das System an der Stabilitätsgrenze. Dieser Wert von 1,000 bedeutet, dass Störungen lediglich begrenzte periodische Schwingungen hervorrufen, während die Gesamttopologie der Bahn intakt bleibt. Ein Multiplikator von 1,000 ist der mathematische Anker, der die Existenz dieser Bahnen im Ozean des Chaos absichert.
Um den mechanischen Unterschied zwischen stabilen und instabilen Bahnen unmittelbar erfahrbar zu machen, bietet die Website auf der rechten Seite einen “Kick”-Schieberegler. Verändert man diesen im Browser, um einem der Körper einen winzigen Geschwindigkeitsimpuls zu versetzen, übernimmt sofort der lokale Echtzeitintegrator die Steuerung von den vorberechneten Splines.
Bei linear stabilen Bahnen wie Broucke R12 führt eine kleine Störung lediglich dazu, dass die Körper kurz um ihre Sollbahn oszillieren; dank intrinsischer Rückstellkräfte findet das Dreikörpersystem rasch ins Gleichgewicht zurück. Bei einer instabilen Bahn hingegen genügt der kleinste Schubs des Reglers: Innerhalb weniger Zyklen zerfällt das System und die Himmelskörper werden unkontrolliert in die Weiten des Alls geschleudert. Dieses minimalistische Interaktionsdesign übersetzt die abstrakte Eigenwertanalyse von Differentialgleichungen in eine greifbare physikalische Intuition. In den Diskussionsforen von Hacker News entbrannten unter Physikern und Softwareentwicklern lebhafte Debatten über die Grenzwerte zulässiger Perturbationen.
Vom Rechenaufwand zur dezentralen Crowdsourcing-Suche
Kürzlich führte die Plattform eine browserbasierte Rechenbeteiligung ein: Über den “Hunt”-Modus können Nutzer ungenutzte Rechenkapazität bereitstellen, um innerhalb definierter Parameterbereiche nach bislang unentdeckten Bahnen zu suchen. Konvergiert ein Rechner auf eine neue, linear stabile periodische Lösung, darf der Finder dem algorithmischen Planetensystem einen Namen geben. So wurde aus einem trockenen mathematischen Suchlauf ein dezentrales Geek-Crowdsourcing.
In der Populärkultur gilt die Unlösbarkeit des Dreikörperproblems oft als unabwendbares Dogma des Chaos. Dieser Atlas beweist jedoch: Wenn moderne Algorithmen geschlossene Lösungen auf 24 Dezimalstellen genau verifizieren können, treten aus dem Chaos fast viertausend Inseln exakter mathematischer Ordnung hervor. Das Dreikörperproblem hat die Schranken philosophischer Ausweglosigkeit längst hinter sich gelassen und ist zu einem lebendigen Testfeld für datengestützte Spitzenrechenleistung geworden.
Referenzen:
- Three Body Orbits Website
- Physical Review Letters 110, 114301
- Science China Physics, Mechanics & Astronomy 60, 129511
- Sofia University Research Group Publications
- Celestial Mechanics 12, 439-462
- Hacker News Diskussion (item?id=49670852)
