Roche-Grenze

Dieser Artikel wurde in die Qualitätssicherung der Redaktion Physik eingetragen. Wenn du dich mit dem Thema auskennst, bist du herzlich eingeladen, dich an der Prüfung und möglichen Verbesserung des Artikels zu beteiligen. Der Meinungsaustausch darüber findet derzeit nicht auf der Artikeldiskussionsseite, sondern auf der Qualitätssicherungs-Seite der Physik statt.

Die Roche-Grenze [ʀɔʃ-] beschreibt diejenige Umlaufbahn zweier idealisierter, theoretischer Körper, bei der sich ein perfektes Gleichgewicht von nach innen gerichteten Gravitationskräften und den nach außen wirkenden Gravitationskräften („Gezeitenkräfte“) einstellt. Innerhalb dieser Zone wird angenommen, dass die nach außen wirkenden Kräfte die Oberhand haben und es somit für Materie schwer ist eine kohäsive Form einzunehmen bzw. zu halten. Die Roche-Grenze ist nach Édouard Albert Roche benannt, der sie 1850 im Zusammenhang mit seinen Überlegungen zur Erklärung der Bildung planetarer Ringsysteme formulierte.

Die mathematisch konstruierbare Umlaufbahn stellt den frühestmöglichen Abstand dar, ab dem sich größere Materie-Ansammlungen bilden können, die zu eigenständigen Körpern, den Trabanten oder Monden, werden können. Dies ist nicht gleichbedeutend mit dem umgekehrten Fall, dass sich jedwede, annähernde Materie beim Übertreten der Roche-Grenze auflöst. Nur wenn ein realer Körper innerhalb dieser Zone eine stabile Umlaufbahn einnehmen und halten kann, ist es theoretisch denkbar, dass die Zeit in kosmischen Maßstäben ausreichend ist, damit größere Materie-Klumpen zu kleineren zerfallen können, bis schließlich nur noch kleinste Partikel vorhanden sind. Partikel innerhalb der Roche-Grenze würden sich in ringförmigen Bahnen um den Zentralkörper ausbilden, etwa so wie es bei den Ringen des Saturn zu sehen ist.

In der Realität spielen jedoch sehr viele weitere Faktoren eine Rolle, die bestimmen ob oder ab wann ein Körper beginnen könnte zu zerfallen. Diese werden unter dem Begriff „innere Stabilität“ zusammengefasst. Sie beinhaltet alle Kräfte die den Körper über die bloßen Gravitationskräfte hinaus zusammenhalten, bei der Betrachtung bzw. Berechnung der Roche-Grenze aber nicht berücksichtigt werden. Es ist daher wahrscheinlicher, dass sich ein Himmelskörper (der in jedem Fall ein Maß an innerer Stabilität aufweist das von Null verschieden ist), auf Grund der Anziehungskräfte des Zentralkörpers, auf einer Abwärtsspiralbahn und Kollisionskurs befindet. Die Roche-Grenze beschreibt dann den niedrigsten Orbit einer vollständigen Umrundung, bevor der Einschlag (nicht die vollständige Auflösung des Körpers) auf den Zentralkörper stattfindet. Ein solcher Fall konnte 1994 beobachtet werden, als der Komet Shoemaker-Levy 9, nach etwa 40 Jahren der Umrundung des Jupiters in mehrere große Stücke zerbrach, die mit Jupiter kollidierten. Das Auseinanderbrechen geschah (weit) innerhalb der theoretisch ermittelbaren Roche-Grenze, jedoch nicht exakt an der Roche-Grenze selbst.

Das im Allgemeinen betrachtete System aus Zentralkörper und ihn umkreisenden Körper, lässt sich auf viele Szenarien übertragen. Der naheliegendste Fall ist ein Planet-Mond-System. Es lassen sich aber auch Sonne-Planet, Mond-Asteroid, Partikel-Partikel, Planet-Planet, Sonne-Materiewolke, Sonne-Sonne etc. betrachten.

Der Zentralkörper ist dabei der größere bzw. massereichere Körper und übt daher starke Gezeitenkräfte aus. Der kleinere Körper wird als der umkreisende und durch diese Kräfte verformende Körper angenommen.

Ursache der Gezeitenkräfte ist der Umstand, dass die Anziehungskraft durch den Zentralkörper auf der ihm zugewandten Seite des Himmelskörpers größer ist als auf der abgewandten. Daher kommt es zu inneren Spannungen oder Verformungen, die bis zur Auflösung des Himmelskörpers führen können.

Die Roche-Grenze eines Himmelskörpers, der einen Hauptkörper umkreist, ist die Entfernung, in der der Körper aufgrund der Gezeitenkräfte, die auf ihn wirken, zerrissen wird.

Es wird dabei angenommen, dass der Körper nur von den eigenen Gravitationskräften zusammengehalten wird und seine mechanische Festigkeit vernachlässigbar ist. Für reale Festkörper ist diese Annahme umso besser erfüllt, je größer der Körper ist. Daher können künstliche Satelliten problemlos auch innerhalb der Roche-Grenze kreisen, während große Objekte wie Monde und Planeten dort nicht existieren können.

Befindet sich Material, das sich noch nicht zu einem Einzelkörper zusammengeballt hat, auf einer Umlaufbahn um den Hauptkörper, so wird sich dieses Material innerhalb der Roche-Grenze ringförmig auf dem Orbit verteilen, während es außerhalb der Grenze einen Klumpen bildet.

Tatsächlich befinden sich fast alle bekannten Planetenringe innerhalb der Roche-Grenze ihrer Planeten. Sie könnten sich daher entweder aus der protoplanetaren Akkretionsscheibe direkt geformt haben, da die Gezeitenkräfte verhindert haben, dass sich Monde aus diesem Material formen, oder Bruchstücke zerstörter Monde sein, die sich von außen über die Roche-Grenze bewegt haben. Alle größeren Monde des Sonnensystems befinden sich dagegen weit außerhalb der Roche-Grenze, kleinere Monde sind jedoch in der Lage, sich auch innerhalb der Roche-Grenze aufzuhalten. So befinden sich die Bahnen von Jupiters Mond Metis und Saturns Mond Pan innerhalb der Roche-Grenze für so genannte flüssige Körper. Die mechanische Festigkeit dieser Körper wirkt dabei zum einen direkt gegen die Gezeitenkräfte, die auf den Körper einwirken, und zum anderen bewirkt die Festigkeit auch, dass diese Körper starr bleiben, d. h. ihre Form nicht verändern – ein Effekt, der weiter unten beschrieben wird und die Gezeitenkräfte zusätzlich verstärkt. Dieser Effekt wird besonders anschaulich durch die Tatsache beschrieben, dass ein Objekt, das auf die Oberfläche eines solchen Mondes „gelegt“ würde, sich nicht auf dem Mond halten würde, sondern durch die Gezeitenkräfte von der Oberfläche fortgezogen würde. Ein Körper mit geringerer mechanischer Festigkeit, beispielsweise ein Komet, würde in diesen Regionen zerstört werden, wie am Beispiel des Kometen Shoemaker-Levy 9 zu sehen war, dessen Orbit im Juli 1992 Jupiters Roche-Grenze durchstieß, woraufhin der Kern des Kometen in zahlreiche Fragmente zerfiel. Bei der nächsten Annäherung an den Planeten im Jahre 1994 kollidierten diese Fragmente dann mit dem Planeten.

Die Roche-Grenze hängt von der Verformbarkeit des Satelliten ab, der sich dem Hauptkörper nähert. Zur Berechnung dieser Grenze werden daher zwei extreme Fälle betrachtet. Im ersten Fall nimmt man an, dass der Körper absolut starr bleibt, bis der Körper von den Gezeitenkräften zerrissen wird. Der gegenteilige Fall ist ein so genannter „flüssiger Körper“, d. h. ein Satellit, der sich der Verformung überhaupt nicht widersetzt und daher bei der Annäherung an die Roche-Grenze zunächst länglich verformt wird und dann zerreißt. Der zweite Fall liefert erwartungsgemäß den größeren Abstand zum Planeten als Roche-Grenze.

Beim starren Satelliten wird angenommen, dass die inneren Kräfte die Form des Körpers stabil halten, der Körper aber trotzdem nur durch seine Eigengravitation zusammengehalten wird. Weitere Idealisierungen sind die Vernachlässigung eventueller Verformungen des Hauptkörpers durch Gezeitenkräfte oder dessen Eigenrotation sowie die Eigenrotation des Satelliten. Die Roche-Grenze ${\displaystyle d}$ ist in diesem Fall

${\displaystyle d=R\cdot {\sqrt[{3}]{\frac {2\,\rho _{M}}{\rho _{m}}}}\approx 1{,}25992\cdot R\cdot {\sqrt[{3}]{\frac {\rho _{M}}{\rho _{m}}}}}$

wobei ${\displaystyle R}$ der Radius und ${\displaystyle \rho _{M}}$ die Dichte des Hauptkörpers ist, und ${\displaystyle \rho _{m}}$ die Dichte des Satelliten beschreibt.

Man bemerkt an der obigen Formel, dass die Roche-Grenze innerhalb des Hauptkörpers liegt, wenn der Satellit ein starrer Körper ist und seine Dichte mindestens doppelt so hoch ist wie die Dichte des Hauptkörpers. Dieser Fall tritt z. B. bei vielen felsigen Monden der Gasriesen unseres Sonnensystems auf. Solche Satelliten werden also auch bei nächster Annäherung an den Hauptkörper nicht durch dessen Gezeitenkräfte zerrissen.

Herleitung der Formel

Um die obige Formel herzuleiten, nehmen wir an, dass eine kleine Masse ${\displaystyle u}$ auf der Oberfläche des Satelliten an dem Punkt liege, der dem Hauptkörper am nächsten ist. Der Satellit selbst wird bei dieser Herangehensweise als sphärisch betrachtet und hat einen Radius r sowie eine Masse m. Auf die kleine Masse m' , die auf der Oberfläche liegt, wirken nun zwei Kräfte:

• die Gravitationskraft, mit der der Satellit die auf seiner Oberfläche liegende Masse m' anzieht:

${\displaystyle F_{G}={\frac {G\,m\,m'}{r^{2}}}}$

• die Gezeitenkraft, die auf die Masse m' wirkt, da sie vom Hauptkörper angezogen wird, sich aber nicht im Schwerpunkt des Satelliten befindet, der sich im freien Fall (Orbit) um den Hauptkörper bewegt. Im mit dem Satelliten mitbewegten Bezugssystem kann man diese Gezeitenkraft auch als Differenz der Gravitationskraft, die vom Hauptkörper auf die Masse m' ausgeübt wird, und der Trägheitskraft betrachten. Für sie ergibt sich in 1. Näherung

${\displaystyle F_{T}={\frac {2\,G\,M\,m'\,r}{d^{3}}}}$

Die Roche-Grenze wird erreicht, wenn der kleine Testkörper auf der Oberfläche des Satelliten anfängt zu schweben, d. h., wenn die Gravitationskraft ${\displaystyle F_{G}}$ und die Gezeitenkraft ${\displaystyle F_{T}}$ den gleichen Betrag annehmen. Für diesen Fall erhält man aus den obigen Gleichungen die Beziehung

${\displaystyle d=r\,{\sqrt[{3}]{\frac {2\,M}{m}}}}$

die die Testmasse m' nicht mehr enthält. Drückt man die Massen der beiden Himmelskörper durch ihre durchschnittlichen Dichten ${\displaystyle \rho _{M}}$ und ${\displaystyle \rho _{m}}$ und ihre Radien ${\displaystyle R}$ und ${\displaystyle r}$ aus, so erhält man die obige von Masse und Radius des Satelliten unabhängige Beziehung.

Das Modell eines flüssigen Satelliten, der den Hauptkörper umkreist, bildet den im Vergleich zu einem starren Satelliten entgegengesetzten Grenzfall. Flüssig bedeutet dabei, dass der Satellit sich der Verformung durch Gezeitenkräfte überhaupt nicht entgegensetzt. (Außer durch die von ihm erzeugte Gravitation. Oberflächenspannung und anderes wird vernachlässigt.) Die Gezeitenkräfte führen dann zu einer länglichen Verformung des Satelliten in Richtung der Verbindungslinie von Satellit und Hauptkörper. Tatsächlich ist dies genau der Effekt, den wir auf der Erde als Gezeiten kennen, bei dem sich die flüssigen Ozeane auf der Oberfläche der Erde in Richtung der Verbindungslinie zum Mond verformen und zwei Flutberge bilden. Da die Stärke der Gezeitenkraft mit der Ausdehnung des Körpers in die Richtung der Verbindungslinie wächst, sorgt eine starke Verformung des Satelliten allerdings für eine noch größere Gezeitenkraft. Daher ist die Roche-Grenze für den Bahnradius eines flüssigen Satelliten wesentlich größer, als wir im starren Modell berechnet haben, und zwar:

${\displaystyle d\approx 2{,}423\cdot R\cdot {\sqrt[{3}]{\frac {\rho _{M}}{\rho _{m}}}}}$

d. h. ungefähr doppelt so groß wie die des starren Modells. Roche hat diese Grenzdistanz bereits um das Jahr 1850 (siehe Literatur) berechnet, den numerischen Faktor in der Formel mit 2,44 aber etwas zu hoch angesetzt. Die Roche-Grenze realer Satelliten liegt zwischen den beiden Grenzmodellen und hängt von Größe und Starrheit des Satelliten ab.

Herleitung der Formel

Um die oben angegebene Formel herzuleiten, ist wesentlich mehr Aufwand vonnöten als im Falle des starren Körpers. Zunächst müssen wir den Begriff des flüssigen Körpers spezifizieren. Gemeint ist damit ein Körper, der aus einer inkompressiblen Flüssigkeit besteht, die also unabhängig von den äußeren und inneren Kräften eine vorgegebene Dichte ρ_m und ein vorgegebenes Volumen V besitzt. Weiterhin gehen wir davon aus, dass sich der Satellit in gebundener Rotation auf einer Kreisbahn bewegt, d. h., sein Schwerpunkt ruht in einem mit fester Winkelgeschwindigkeit ω rotierenden Bezugssystem mit Ursprung im Schwerpunkt des Gesamtsystems. Die Winkelgeschwindigkeit ist dabei durch das dritte Kepler’sche Gesetz gegeben:

${\displaystyle \omega ^{2}={\frac {G(M+m)}{d^{3}}}}$

In diesem Bezugssystem heißt gebundene Rotation des Satelliten, dass sich die Flüssigkeit, aus der der Satellit besteht, nicht bewegt, das Problem kann also als statisch angesehen werden. Daher spielen auch die Viskosität und Reibung der Flüssigkeit in diesem Modell keine Rolle, da diese Größen nur bei Bewegung des Fluids in die Rechnung eingehen würden.

Auf die Flüssigkeit des Satelliten wirken nun im rotierenden Bezugssystem folgende Kräfte:

• die Gravitationskraft des Hauptkörpers
• die Zentrifugalkraft als Scheinkraft im rotierenden Bezugssystem
• die Gravitationskraft des Satelliten selbst

Da alle vorkommenden Kräfte konservativ sind, können sie alle durch ein Potential dargestellt werden. Die Oberfläche des Satelliten nimmt dabei die Form einer Äquipotentialfläche des Gesamtpotentials an, da es andernfalls eine Horizontalkomponente der Kraft an der Oberfläche geben würde, der Teile der Flüssigkeit folgen würden. Welche Form der Satellit bei gegebener Distanz zum Hauptkörper annehmen muss, damit diese Forderung erfüllt ist, soll nun diskutiert werden.

Wir wissen bereits, dass sich die Gravitationskraft des Hauptkörpers und die Zentrifugalkraft im Schwerpunkt des Satelliten aufheben, da dieser sich auf einer (frei fallenden) Kreisbahn bewegt. Die äußere Kraft, die auf die Flüssigkeitspartikel wirkt, ist daher vom Abstand zum Schwerpunkt abhängig und ist die schon im starren Modell benutzte Gezeitenkraft. Für kleine Körper ist der Abstand der Flüssigkeitspartikel vom Schwerpunkt klein gegenüber dem Abstand d zum Hauptkörper, und die Gezeitenkraft kann linearisiert werden, wodurch sie die oben angegebene Formel für F_T ergibt. Als Abstand vom Schwerpunkt wurde im starren Modell nur der Radius des Satelliten r betrachtet, betrachtet man aber nun einen beliebigen Punkt auf der Oberfläche des Satelliten, so hängt die dort wirksame Gezeitenkraft vom Abstand Δd des Punktes zum Schwerpunkt in radialer Richtung (d. h. parallel zur Verbindungslinie vom Satelliten zum Hauptkörper) ab. Da die Gezeitenkraft linear im radialen Abstand Δd ist, ist ihr Potential in dieser Variable quadratisch, und zwar ergibt sich (für ${\displaystyle m\ll M}$):

${\displaystyle V_{T}=-{\frac {3\,GM}{2\,d^{3}}}\Delta d^{2}\,}$

Wir suchen nun also eine Form für den Satelliten, so dass sich sein Eigengravitationspotential gerade diesem Gezeitenpotential so überlagert, dass das Gesamtpotential auf der Oberfläche konstant wird. Ein solches Problem ist im Allgemeinen sehr schwierig zu lösen, durch die einfache quadratische Abhängigkeit des Gezeitenpotentials vom Abstand vom Schwerpunkt kann die Lösung dieses Problems aber glücklicherweise durch geschicktes Raten gefunden werden.

Da sich das Gezeitenpotential nur in eine Richtung, nämlich in die Richtung zum Hauptkörper, verändert, ist es naheliegend, dass der Satellit bei seiner Verformung axialsymmetrisch um diese Verbindungslinie bleibt, also einen Rotationskörper bildet. Das Eigenpotential eines solchen Rotationskörpers auf der Oberfläche kann dann nur vom radialen Abstand zum Schwerpunkt abhängen, da die Schnittfläche eines solchen Körpers bei festem radialen Abstand ja gerade eine Kreisscheibe ist, deren Rand sicherlich konstantes Potential besitzt. Sollen nun die Summe des Eigenpotentials und des Gezeitenpotentials an jedem Punkt der Oberfläche gleich sein, muss das Eigenpotential genau wie das Gezeitenpotential eine quadratische Abhängigkeit vom radialen Abstand besitzen. Es zeigt sich, dass man dann als Form ein prolates (zigarrenförmiges) Rotationsellipsoid wählen muss. Bei vorgegebener Dichte und Volumen ist das Eigenpotential eines solchen Ellipsoids abhängig von der numerischen Exzentrizität ε des Ellipsoids:

${\displaystyle V_{S}=V_{S0}+G\pi \rho _{m}f(\epsilon )\Delta d^{2}}$

wobei ${\displaystyle V_{S0}}$ das konstante Eigenpotential auf dem kreisförmigen Rand der zentralen Symmetrieebene bei Δd = 0 ist. Die dimensionslose Funktion f ist aus der exakten Lösung des Potentials eines Rotationsellipsoids zu bestimmen und ergibt sich zu

${\displaystyle f(\epsilon )={\frac {\epsilon ^{2}-1}{\epsilon ^{3}}}\left(3\epsilon +\left(\epsilon ^{2}-3\right)\cdot \operatorname {arsinh} {\frac {\epsilon }{\!\!{\sqrt {1-\epsilon ^{2}}}}}\!\right)}$

und hängt erstaunlicherweise nicht vom Volumen des Satelliten ab.

So kompliziert die Abhängigkeit der Funktion f von der Exzentrizität ist, brauchen wir nun dennoch nur noch den geeigneten Wert für die Exzentrizität zu bestimmen, damit ${\displaystyle V_{T}+V_{S}}$ konstant in der einzigen Ortsvariable Δd ist. Dies ist genau dann der Fall, wenn

${\displaystyle {\frac {2\,G\pi \rho _{M}R^{3}}{d^{3}}}=G\pi \rho _{m}f(\epsilon )}$

ist, eine Gleichung, die jeder Computer numerisch leicht lösen kann. Wie man am Verlauf der Funktion f in nebenstehender Grafik entnehmen kann, hat diese Gleichung im Allgemeinen zwei Lösungen, wobei die kleinere Lösung, also die geringere Exzentrizität, die stabile Gleichgewichtslage darstellt. Diese Lösung der Gleichung gibt daher die Exzentrizität des Gezeitenellipsoids an, das sich bei einer festgelegten Entfernung zum Hauptkörper einstellt.

Die Roche-Grenze entsteht nun durch die Tatsache, dass die Funktion f, die man als die Stärke der Kraft, die das Ellipsoid in die kugelförmige Gestalt zurückformen will, ansehen kann, nicht beliebig groß werden kann. Es gibt eine gewisse Exzentrizität, bei der diese Kraft maximal wird. Da die Gezeitenkraft aber bei Annäherung an den Hauptkörper über alle Grenzen steigen kann, ist klar, dass es einen Grenzabstand gibt, bei dem das Ellipsoid zerrissen wird.

Die maximale Exzentrizität des Gezeitenellipsoids errechnet man numerisch aus der Nullstelle der Ableitung der Funktion f, die in der Grafik dargestellt ist. Man erhält

${\displaystyle \epsilon _{\text{max}}\approx 0{,}85883075}$

was einem Achsenverhältnis von ${\displaystyle 1\!:\!{\tfrac {1}{\sqrt {1-{\epsilon _{\text{max}}}^{2}}}}\approx 1\!:\!1{,}9521}$ entspricht. Setzt man diesen Wert in die Funktion f ein, kann man den minimalen Abstand berechnen, in dem ein solches Gezeitenellipsoid existiert – die Roche-Grenze:

${\displaystyle d\approx 2{,}423\cdot R\,{\sqrt[{3}]{\frac {\rho _{M}}{\rho _{m}}}}}$

Die oben angegebenen Werte werden nun benutzt, um die Roche-Grenzen für das starre Modell und das flüssige Modell zu berechnen. Als mittlere Dichte eines Kometen wird 500 kg/m³ angenommen. Die wahre Roche-Grenze hängt von der Flexibilität des jeweiligen Satelliten, aber auch von zahlreichen anderen Parametern, wie der Verformung des Hauptkörpers und der genauen Dichteverteilung innerhalb des Satelliten, ab und liegt normalerweise zwischen den beiden angegebenen Werten.

Man sieht an der obigen Tabelle, dass für besonders dichte Satelliten, die einen weit weniger dichten Hauptkörper umkreisen, die Roche-Grenze innerhalb des Hauptkörpers liegen kann (z. B. beim System Sonne-Erde). In der nächsten Tabelle werden noch einige weitere Beispiele vorgestellt, wobei der tatsächliche Abstand des Satelliten in Prozent der Roche-Grenze angegeben wird. Man sieht z. B., dass der Neptunmond Naiad besonders nahe der Roche-Grenze des starren Modells liegt und daher wohl seiner tatsächlichen physikalischen Roche-Grenze bereits recht nahe ist.

In diesem Fall bezeichnet die Roche-Grenze den kleinstmöglichen Abstand, bei welchem ein kleiner Himmelskörper, der einen größeren, massereicheren Hauptkörper umkreist, noch existieren kann, ohne von den Gezeitenkräften, also der Anziehungskraft des Hauptkörpers, zerrissen zu werden.

Umkreist ein Stern einen Partner, so wird er durch die Gezeitenkräfte deformiert. Ist der Stern groß und nah genug, so nimmt er eine Tropfenform an mit einer Spitze, die dem Partner zugewandt ist. Befindet er sich in einer Expansionsphase, wie beispielsweise in der Übergangsphase zu einem roten Riesen, so kann er nicht weiter wachsen, sondern es fließt Material über diese Spitze zum Partner. Diese Tropfenform wird ebenfalls als Roche-Grenze bezeichnet. Da dieser Masseverlust die Roche-Grenze (für die Form des Umkreisenden) verkleinert, kann das ganze System instabil werden und der Stern komplett zu seinem Partner hinüberfließen.

Handelt es sich bei dem Partner um ein kompaktes Objekt, wie einen weißen Zwerg, einen Neutronenstern oder ein schwarzes Loch, so spielen sich beim Materialtransfer dramatische Prozesse ab. Siehe dazu Novae und Röntgendoppelsterne.

Die Roche-Grenze des Gesamtsystems setzt sich aus den zwei tropfenförmigen Äquipotentialflächen zusammen, die sich an den Spitzen berühren und so die Form einer Acht bilden. Diese Spitze ist der so genannte Lagrange-Punkt L1 des Systems. Diese Potentialfläche muss für ein mitrotierendes Koordinatensystem berechnet werden. Es handelt sich um ein effektives Potential, das neben den Gravitationskräften auch die Zentrifugalkräfte berücksichtigt. Sobald sich Material in diesem System bewegt, erfährt es zusätzlich eine Corioliskraft, die jedoch nur durch ein geschwindigkeitsabhängiges Potential beschrieben werden kann.^[1]

• Hill-Sphäre

• Édouard Roche: La figure d’une masse fluide soumise à l’attraction d’un point éloigné. Acad. des sciences de Montpellier, Vol. 1 (1847–50), S. 243
• Roche, É. (1873). Essai sur la constitution et l'origine du système solaire. Frankreich: Gauthier-Villars.

• Detaillierte Ableitung der Roche-Grenze Scienceworld Wolfram (allerdings sehr unübersichtlich)

1. ↑ Roche-Grenze auf Spektrum.de. In: Lexikon der Physik, Spektrum.de. Abgerufen am 5. Dezember 2021.