roche limit

The Roche limit is the distance within which a celestial body held together only by its own gravity will disintegrate due to a second celestial body's tidal forces exceeding the first body's gravitational self-attraction. Inside the Roche limit, orbiting material will tend to disperse and form rings, while outside the limit, material will tend to coalesce. The term is named after douard Roche, the French astronomer who first discovered this theoretical limit in 1848. The Roche limit should not be confused with the concept of the Roche lobe which is also named after douard Roche and which describes the limits at which an object which is in orbit around two other objects will be captured by one or the other. Typically, the Roche limit applies to a satellite disintegrating due to tidal forces induced by its primary, the body about which it orbits. Some real satellites, both natural and artificial, can orbit within their Roche limits because they are held together by forces other than gravitation. Jupiter's moon Metis and Saturn's moon Pan are examples of natural satellites which are able to hold together despite being within their fluid Roche limits. They hold together partly because of their tensile strength, and partly because they are not actually fluid. In such cases, it is possible for an object resting on the surface of such a satellite to be pulled away by tidal forces, depending on where it is: tidal forces are most repulsive along the line of centers between the satellite and primary. A weaker satellite, such as a comet, could be broken up when it passes within its Roche limit. Comet Shoemaker-Levy 9's decaying orbit around Jupiter passed within its Roche limit in July 1992, causing it to fragment into a number of smaller pieces. On its next approach in 1994 the fragments crashed into the planet. Since tidal forces overwhelm gravity within the Roche limit, no large satellite can coalesce out of smaller particles within that limit. Indeed, all known planetary rings are located within their Roche limit. They could either be remnants from the planet's proto-planetary accretion disc that failed to coalesce into moonlets, or conversely have formed when a moon passed within its Roche limit and broke apart.

Determining the Roche limit

The Roche limit depends on the rigidity of the satellite. At one extreme, a rigid satellite will maintain its shape until tidal forces break it apart. At the other extreme, a highly fluid satellite gradually deforms with increasing tidal forces until it breaks apart. For a rigid spherical satellite, the cause of the rigidity is neglected, in that the material constituting the satellite is still treated as though held together only by its own self-gravity. Other effects are also neglected, such as tidal deformation of the primary, and rotation of the satellite. The Roche limit,

d

, is then the following:

d = R\left( 2\;\frac {\rho_M} {\rho_m} \right)^{\frac{1}{3}} \approx  1.260R\left( \frac {\rho_M} {\rho_m} \right)^{\frac{1}{3}}

where

R

is the primary's radius,

\rho_M

is the primary's density and

\rho_m

is the satellite's density. For a fluid satellite, tidal forces cause the satellite to elongate, further compounding the tidal forces and causing it to break apart more readily. The calculation is complex and cannot be solved exactly, but a close approximation is the following:

d \approx  2.423R\left( \frac {\rho_M} {\rho_m} \right)^{\frac{1}{3}}

which indicates that a fluid satellite will disintegrate at almost twice the distance of a rigid sphere of similar density. Most real satellites are somewhere between these two extremes, with internal friction, viscosity, and chemical bonds rendering the satellite neither perfectly rigid nor perfectly fluid.

Rigid satellites

As stated above, the formula for calculating the Roche limit,

d

, for a rigid spherical satellite orbiting a spherical primary is:

d = R\left( 2\;\frac {\rho_M} {\rho_m} \right)^{\frac{1}{3}}

where

R

is the radius of the primary,

\rho_M

is the density of the primary, and

\rho_m

is the density of the satellite. As described below, this rigid-body approximation does not take into account the deformation of the satellite's spherical shape due to tidal effects and is only an approximation of what a real satellite's Roche limit would be. Notice that if the satellite is more than twice as dense as the primary (as can easily be the case for a rocky moon orbiting a gas giant) then the Roche limit will be inside the primary and hence not relevant.

Derivation of the formula

In order to determine the Roche limit, we consider a small mass

u

on the surface of the satellite closest to the primary. There are two forces on this mass

u

: the gravitational pull towards the satellite and the gravitational pull towards the primary. Since the satellite is already in orbital free fall around the primary, the tidal force is the only relevant term of the gravitational attraction of the primary. The gravitational pull

F_G

on the mass

u

towards the satellite with mass

m

and radius

r

can be expressed according to Newton's law of gravitation.

F_G = \frac{Gmu}{r^2}

The tidal force

F_T

on the mass

u

towards the primary with radius

R

and a distance

d

between the center of the two bodies can be expressed as:

F_T = \frac{2GMur}{d^3}

The Roche limit is reached when the gravitational pull and the tidal force cancel each other out.

F_G = F_T

\frac{Gmu}{r^2} = \frac{2GMur}{d^3}

Which quickly gives the Roche limit, d, as:

d = r \left( 2 M / m \right)^{\frac{1}{3}}

However, we don't really want the radius of the satellite to appear in the expression for the limit, so we re-write this in terms of densities. For a sphere the mass

M

can be written as:

M = \frac{4\pi\rho_M R^3}{3}

where

R

is the radius of the primary.

And likewise:

m = \frac{4\pi\rho_m r^3}{3}

where

r

is the radius of the satellite.

Substiting for the masses in the equation for the Roche limit, and cancelling out

4\pi/3

gives:

d = r \left( \frac{ 2 \rho_M R^3 }{ \rho_m r^3 } \right)^{1/3}

which can be simplified to the Roche limit:

d = R\left( 2\;\frac {\rho_M} {\rho_m} \right)^{\frac{1}{3}}

Fluid satellites

A more correct approach for calculating the Roche Limit takes the deformation of the satellite into account. An extreme example would be a tidally locked liquid satellite orbiting a planet, where any force acting upon the satellite would deform the satellite. In this case, the satellite is deformed into a prolate spheroid. The calculation is complex and cannot be solved exactly. Historically, Roche himself derived the following numerical solution for the Roche Limit:

d \approx  2.44R\left( \frac {\rho_M} {\rho_m} \right)^{1/3}

However, with the aid of a computer a better numerical solution is:

d \approx 2.423 R\left( \frac {\rho_M} {\rho_m} \right)^{1/3} \left( \frac{(1+\frac{m}{3M})+\frac{c}{3R}(1+\frac{m}{M})}{1-c/R} \right)^{1/3}

where

c/R

is the oblateness of the primary.

Roche limits for selected examples

The table below shows the mean density and the equatorial radius for selected objects in our solar system.

Primary !! Density (kg/m³) !! Radius (m)
Sun	align="center"\| 1,400	align="right"\|695,000,000
Jupiter	align="center"\|1,330	align="right"\|71,500,000
Earth	align="center"\| 5,515	align="right"\| 6,376,500
Moon	align="center"\| 3,340	align="right"\|1,737,400

Using these data, the Roche Limits for rigid and fluid satellites can easily be calculated. The average density of comets is around 500 kg/m³. The table below gives the Roche limits expressed in metres and in primary radii. The true Roche Limit for a satellite depends on its flexibility, and will be somewhere between the rigid and fluid Roche Limits given below.

rowspan="2"\| Body !!rowspan="2"\| Satellite !!colspan="2"\| Roche limit (rigid) !!colspan="2"\| Roche limit (fluid)
Distance (m) !! Radii !! Distance (m) !! Radii
Earth	Moon	align="right"\| 9,495,665	align="center"\| 1.49	align="right"\| 18,261,459	align="center"\| 2.86
Earth	Comet	align="right"\| 17,883,432	align="center"\| 2.80	align="right"\| 34,392,279	align="center"\| 5.39
Sun	Earth	align="right"\| 554,441,389	align="center"\| 0.80	align="right"\| 1,066,266,402	align="center"\| 1.53
Sun	Jupiter	align="right"\| 890,745,427	align="center"\| 1.28	align="right"\| 1,713,024,931	align="center"\| 2.46
Sun	Moon	align="right"\| 655,322,872	align="center"\| 0.94	align="right"\| 1,260,275,253	align="center"\| 1.81
Sun	Comet	align="right"\| 1,234,186,562	align="center"\| 1.78	align="right"\| 2,373,509,071	align="center"\| 3.42

If the primary is less than half as dense as the satellite, the rigid-body Roche Limit is less than the primary's radius, and the two bodies may collide before the Roche limit is reached. For example, the Sun-Earth Roche Limit indicates that the Earth would collide with the Sun before disintegrating due to tidal forces. How close are the solar system's moons to their Roche limits? The table below gives each inner satellite's orbital radius divided by its own Roche radius, for both the rigid and fluid cases.

rowspan="2"\| Primary !!rowspan="2"\| Satellite !!colspan="2"\| Orbital Radius vs. Roche limit
(rigid) !! (fluid)
Sun	Mercury	align="center"\| 104:1	align="center"\| 54:1
Earth	Moon	align="center"\| 41:1	align="center"\| 21:1
owspan="2"\| Mars	Phobos	align="right"\| 172%	align="right"\| 89%
Deimos	align="right"\| 451%	align="right"\| 233%
owspan="4"\| Jupiter	Metis	align="right"\| 186%	align="right"\| 93%
Adrastea	align="right"\| 220%	align="right"\| 110%
Amalthea	align="right"\| 228%	align="right"\| 114%
Thebe	align="right"\| 260%	align="right"\| 129%
owspan="5"\| Saturn	Pan	align="right"\| 174%	align="right"\| 85%
Atlas	align="right"\| 182%	align="right"\| 89%
Prometheus	align="right"\| 185%	align="right"\| 90%
Pandora	align="right"\| 185%	align="right"\| 90%
Epimetheus	align="right"\| 198%	align="right"\| 97%
owspan="4"\| Uranus	Cordelia	align="right"\| 155%	align="right"\| 79%
Ophelia	align="right"\| 167%	align="right"\| 86%
Bianca	align="right"\| 184%	align="right"\| 94%
Cressida	align="right"\| 192%	align="right"\| 99%
owspan="5"\| Neptune	Naiad	align="right"\| 140%	align="right"\| 72%
Thalassa	align="right"\| 149%	align="right"\| 77%
Despina	align="right"\| 153%	align="right"\| 78%
Galatea	align="right"\| 184%	align="right"\| 95%
Larissa	align="right"\| 220%	align="right"\| 113%
Pluto	Charon	align="center"\| 14:1	align="center"\| 7.2:1