Ich habe mich früher auch eine Zeit lang mit der Langzeitstabilität des Sonnensystems beschäftigt, unter anderem auch mit der Arbeit von Laskar. Eine Kollegin von mir aus Wien  hat sogar ihre Diplomarbeit über das chaotische Verhalten von Merkur geschrieben (sie liest hier mit soweit ich weis – vielleicht schreibt sie ja einen Kommentar dazu 😉 ). Diese Arbeit von Laskar war zwar sehr wichtig und äußerst interessant – aber es gab einiges dabei zu kritisieren.

Auch hier wurde natürlich eine Näherung benutzt. Normalerweise benutzt man für diese Art von Simulation die newtonschen Gravitationsgleichungen die dann am Computer numerisch gelöst werden. Um aber einen Zeitraum von einigen Milliarden Jahren zu untersuchen brauchen selbst die besten Computer enorm lang. Deshalb entschied sich Laskar, vereinfachte Bewegungsgleichungen zu benutzen. Diese haben aber den Nachteil, sehr ungenau zu werden, wenn die Exzentrizitäten der Planeten sehr groß werden. Es war also nicht klar, wie exakt Laskars Rechnungen wirklich waren. Zusätzlich handelte es sich nicht um eine durchgängige Simulation. Die Rechungen wurden “zusammengestückelt” und stellten nur eine Art “Worst-case”-Szenario dar. Laskar untersuchte nicht nur die Bahn der “echten” Planeten sondern auch die von 4 benachbarten “Klonen” (die Unterschiede in der Anfangsposition waren aber sehr klein – nur 150 Meter). Die Bahn dieser 5 Objekte wurden nun für einige Zeit berechnet – und dann diejenige ausgewählt, die den stärksten Anstieg der Exzentrizität verzeichnete. Diese Bahn zu diesem Zeitpunkt diente dann als neuer Ausgangspunkt; es wurden wieder 4 Klone hergestellt und das ganze Spiel wiederholt. Am Ende ergab sich dann insgesamt der oben erwähnte enorme Anstieg der Exzentrizität des Merkur. Für diese Methodik wurde Laskar des öfteren kritisiert; was er aber zweifelsfrei zeigte, war das potentiell chaotische Verhalten des inneren Sonnensystems bei der Betrachtung von sehr langen Zeiträumen.

Die neuen Ergebnisse

In seiner neuen Arbeit hat Laskar nun die alte Arbeit weitergeführt. Die besseren Computer machten es möglich sehr viel mehr Simulationen durchzuführen als 1994. Deswegen konnten diesmal auch genaue Wahrscheinlichkeiten berechnet werden, die angeben, wie groß die Chancen sind, dass ein bestimmter Planet einen bestimmten Anstieg in der Exzentrizität zeigt. Laut diesen Rechnungen besteht eine Wahrscheinlichkeit von 1 bis 2% das Merkurs Exzentrizität innerhalb der nächsten 5 Milliarden Jahre auf Werte von über 0.6 steigt. Das ist eine relativ große Wahrscheinlichkeit – und ein so hoher Anstieg der Exzentrizität birgt die Möglichkeit einer Kollision mit Venus! Allerdings hat Laskar auch hier wieder die vereinfachten Bewegungsgleichungen verwendet. Kollegen aus Japan hatten ähnliche Simulationen durchgeführt; diesmal aber mit den ungekürzten Gleichungen (Ito & Tanikawa: “Long-term integrations and stability of planetary orbits in our Solar System“, MNRAS 336, 483, 2002). Sie fanden, dass die Exzentrizität von Merkur zwar größer werden kann, aber immer kleiner als 0.35 bleibt – ein Wert, bei dem keine Gefahr einer Kollision besteht. Ito & Tanikawa hatten allerdings in ihrer Simulation keine Effekte der allgemeinen Relativitätstheorie berücksichtigt und auch den gravitativen Einfluss des Mondes vernachlässigt (Das ist nicht unüblich, denn beide Effekte sind sehr, sehr klein). Laskar argumentiert nun aber, dass gerade bei solch langen Zeiträumen diese kleinen Effekte durchaus berücksichtigt werden müssen. Deswegen hat er selbst die Simulationen von Ito & Tanikawa wiederholt – und fand deutlich höhere Werte für die Exzentrizität des Merkur! (Es wurden Exzentrizitäten von bis zu 0.8 erreicht). Er zeigte auch (diesmal wieder mit den vereinfachten Gleichungen) das die Effekte der allgemeinen Relativitätstheorie einen großen Einfluss haben: Bei Berücksichtigung dieser Effekte ist das System wesentlich weniger chaotisch; die Werte von Merkurs Exzentrizität bleiben kleiner.

Alles in allem wieder eine etwas verwirrende Angelegenheit. Wenn die relativistischen Effekte das Chaos abschwächen, dann müsste die Simulation von Ito & Tanikawa größere Exzentrizitäten liefern als die von Laskar – es ist aber genau umgekehrt. Eine Wiederholung der Simulation der Japaner konnte ihre Ergebnisse auch nicht bestätigen. Und wie genau sind die von Laskar benutzten vereinfachten Gleichungen nun wirklich? Laskar selbst führt den Unterschied zwischen den Ergebnissen auf unterschiedliche Anfangswerte zurück. Das ist durchaus möglich, denn in chaotischen Systemen spielen die eine wichtige Rolle. Um die Gültigkeit seiner vereinfachten Gleichungen einzuschätzen müsste man eine Simulation mit den kompletten Gleichung unter Berücksichtigung der relativistischen Effekte und des Mondeinflusses zur Verfügung haben. Die gibt es aber – für solch lange Zeiträume – leider noch nicht. Laskar hat abschließend angekündigt, diese in nächster Zeit selbst durchführen zu wollen.

1 / 2 / 3

Kommentare (12)

  1. #1 L. Carone
    28. April 2008
  2. #2 florian
    29. April 2008

    Danke! SpOn hab ich irgendwie übersehen…

    Die haben ja sogar das gleiche unpassende Bild wie ich verwendet 😉 (Eigentlich stammt das von ner NASA Presseaussendung. Dort wird beschrieben wie in ner Staubscheibe um nen anderen Stern zwei Pluto-große Körper kollidieren und dabei Staub erzeugen, denn man beobachten konnte)

  3. #3 Birgit
    29. April 2008

    Natürlich les ich hier mit und wenn man schon mal “zitiert” wird:

    Ich find es interessant, dass Jaques Laskar nach so langer Zeit doch wieder sein altes Thema aufgegriffen hat und mit den neuen Informationen auch vom japanischen Team hätt ich große Lust meinen alten – leider sehr langsamen – Linux Rechner wieder anzuwerfen und auch meine Ergebnisse nochmal anzusehen und eventuell zu wiederholen bzw. auszuweiten. Mal schaun wann mir meine Hauptbeschäftigung dazu Zeit lässt 😉

    Das Problem der vereinfachten Bewegungsgleichungen habe wir damals schon unserer Arbeitsgruppe diskutiert. Ich bin auf jeden Fall schon gespannt, was die weiteren Berechnungen von Laskar bringen!

  4. #4 Andylee
    1. Mai 2008

    Also ich habe das noch nicht so ganz verstanden.

    Sind diese Sachen so schwer zu berechnen, weil die Daten über die Planeten noch nicht genau genug sind?
    Sind die “Regeln” dieser Simulationen noch nicht perfekt?
    Oder sind es einfach zu viele Daten, als dass man das genau berechnen könnte, zuviele Variablen und Einflüsse?

    Andylee

  5. #5 florian
    2. Mai 2008

    Naja – Daten über die Planeten hat man alle die man braucht. Es liegt vor allem erstmal an den enorm langen Zeitäumen. Bei so einer numerischen Integration muss man sich ja immer Schritt für Schritt von der Gegenwart in die Zukunft vortasten. Die Schrittweite muss dabei klein genug sein, um vernünfige Ergebnisse zu kriegen. Merkur z.b. ist der innerste Planet und bewegt sich deswegen am schnellsten um die Sonne (eine Umrundung dauert ~115 Tage). Um diese Bewegung einigermassen vernünftig auflösen zu können, sollte die Schrittweite in etwa ein paar Tage betragen. Und wenn ich nun eben z.B. 10 Milliarden Jahre bei einer Schrittweite von 2 Tagen berechnen will, dann sind das knapp 2 Billionen Schritte – und in jedem Schritt müssen die Bewegungsgleichungen für alle Planeten numerisch gelöst werden. Das dauert dann halt einfach ein bisschen – selbst bei den schnellen Computern 😉 Und wenns fertig ist, dann hat man erst eine Simulation – um vernünftige Aussagen machen zu können muss man das ganze öfters mit verschiedenen Parametern wiederholen (Laskar hat 1001 Simulationen durchgeführt).

  6. #6 T. T.
    5. Mai 2008

    Also Florian,
    ich bitte Dich, wenn Du nicht ein mal die Infrarotfotos von Uranus (von Miranda ganz zu schweigen) einigermassen logisch wahrnehmen kannst, was willst Du eigentlich noch in der Astronomie? Willst Du eigentlich nur die alten Kamele verbreiten?
    Hast Du keine andere Möglichkeit deine Brötchen zu verdienen?

    Gruß
    TT

  7. #7 florian
    5. Mai 2008

    Ah – der Herr Tumalski! Schön, dass auch sie mein Blog gefunden haben!
    Die Beleidigungen ignorier ich jetzt mal – und auch darüber, wer von uns zur logischen Wahrnehmung unfähig ist sag ich lieber mal nix.

    Übrigens: einer meiner nächsten Blogbeiträge beschäftigt sich mit der Entstehung des Mondes. Ich hoffe, sie ertragen es, wenn ich ihre “Theorie” dann nur am Rande erwähnen werde 😉

  8. #8 L. Carone
    5. Mai 2008

    Ich glaub, womit Andylee hier Probleme hat, ist, dass wir mit chaotischen Systemen rechnen. Also, dass selbst winzigste Unsicherheiten in der Messung am Anfang auf langen Zeitskalen zu immensen Unterschieden führen können. Da aber immer und jede Messung mit Fehlern und Unsicherheiten behaftet ist und immer behaftet sein wird, kann man immer “nur” sagen, welche Szenarien möglich und welche wahrscheinlicher sind als andere. Die Planetenkollision ist eine Möglichkeit, aber eine sehr unwahrscheinliche.

    Das gleiche gilt übrigens auch für Asteroidenbahnen und die Kollisionswahrscheinlichkeit mit der Erde. Für kurze Zeiten kannst Du eine relativ genaue Vorhersage treffen, aber wenn es über längere Zeiträume geht, dann kannst Du mit Hilfe von Computern “ausprobieren”, was alles möglich ist (Simulation), indem Du mit leicht unterschiedlichen Anfangswerten rumrechnest und das machst Du dann ganz oft und daraus bestimmst Du dann, wie wahrscheinlich bestimmte Szenarien sind. Welches von den durchgespielten wirklich eintritt, das kann Dir genausowenig jemand sagen, wie jemand vorhersagen kann, was Du als nächstes würfeln wirst. Ich kann Dir aber sagen, dass Du mit einem sechsseitigen Würfel unmöglich eine 7 würfeln kannst und dass Du bei 1000 Würfen höchstwahrscheinlich die Zahl sechs ungefähr 167 mal werfen wirst.

    Willkommen in der modernen Wissenschaft, wo wir uns schon vor Jahrzehnten vom reinem Determinismus verabschiedet haben! Aber immerhin sind wir erwachsen genug, damit konstruktiv weiterzuarbeiten, anstatt wie kleine verzogene Gören rumzumaulen, dass das ja alles Sch*** ist mit der Wissenschaft 😉

  9. #9 florian
    5. Mai 2008

    @Ludmila: also besser hätte ichs auch nicht erklären können 😉

    Und die “Theorien” des Herrn Tumalski sind ja eigentlich eher dein Revier: “The photos below prove, that the idea about the spherical symmetry of the inner
    structures of the celestial bodies was one of the biggest mistakes made
    by science in the XX century. “
    – für mich als Dynamiker sind die Himmelskörper sowieso nur Massenpunkte; mit den inneren Aufbau müssen sich die Planetologen rumschlagen… 😉

  10. #10 L. Carone
    5. Mai 2008

    @Florian: Wah! Bist Du des Wahnsinns?

    Naja, dann will ich mal nicht so sein:
    Fotos-Schmotos! Willst Du den inneren Aufbau der Körper bestimmen, gibt es verschiedene Ansätze:
    a) die Rotation des Körpers über einen längeren Zeitraum beobachten, Abplattung, Keiselbewegung etc.
    b) bei einem Stern machst Du Helioseismologie
    c) auf der Erde stellst Du Erdbebendetektoren auf und analysierst die Laufzeit und Wellenform.
    d) Du misst das Schwerefeld des Körpers und das in sehr, sehr hoher Auflösung.
    e) Du simulierst die Entstehung von Planeten und Sternen.

    Hmm, und wer macht d) ? Jepp, wir machen das u.a.

    Mind. 5 verschiedene Ansätze, die alle das gleiche sagen: Spherical symmetry, Baby! And no shit! 😉

    Hmm, wär eigentlich mal wieder ein eigener Blogeintrag wert…

  11. #11 Kurt Kowalsky
    12. Dezember 2010

    So ganz kann ich die Weltuntergangberuhigungsparole nicht nachvollziehen. In etwa 0,9 Milliarden Jahren überschreitet die mittlere Temperatur auf der Erdoberfläche den für höhere Lebewesen kritischen Wert von 30 °C. (Nach einigen Klimaaktivisten ist das aber bereits in zehn Jahren der Fall.) Und jetzt mache ich mir doch ernsthaft Sorgen, dass vielleicht schon 500 Jahre früher, dieser Merkur zu spinnen anfängt. Den stabilisieren wir mit einem Sonnensegel. Wie kann man da so ruhig bleiben? Alle Eigenheimbesitzer werden mir sicherlich zustimmen. 🙂

  12. #12 Thomas J
    12. Dezember 2010

    @Kowalsky

    Du hast keine Ahnung, vor was für Problemen eigenheimbesitzer sonst noch stehen 🙂
    https://www.br-online.de/br-alpha/alpha-centauri/alpha-centauri-mallorca–2002-ID1208358755584.xml