# Roche limit

## Information about Roche limit

 Consider an orbiting mass of fluid held together by gravity, here viewed from above the orbital plane. Far from the Roche limit the mass is practically spherical. Closer to the Roche limit the body is deformed by tidal forces. Within the Roche limit the mass's own gravity can no longer withstand the tidal forces, and the body disintegrates. Particles closer to the primary move more quickly than particles farther away, as represented by the red arrows. The varying orbital speed of the material eventually causes it to form a ring.
The Roche limit, sometimes referred to as the Roche radius, 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.[1] 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 calculated this theoretical limit in 1848.[2]

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 such satellites, which hold together because of their tensile strength. In extreme cases, objects resting on the surface of such a satellite could actually be lifted away by tidal forces. A weaker satellite, such as a comet, could be broken up when it passes within its Roche limit.

Since tidal forces overwhelm gravity within the Roche limit, no large satellite can coalesce out of smaller particles within that limit. Indeed, almost all known planetary rings are located within their Roche limit (Saturn's E-Ring being a notable exception). 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.

(Note that the Roche limit should not be confused with the concept of the Roche lobe or Roche sphere, which are also named after Édouard Roche. The Roche lobe describes the limits at which an object which is in orbit around two other objects will be captured by one or the other. The Roche sphere approximates the gravitational sphere of influence of one astronomical body in the face of perturbations from another heavier body around which it orbits.)

## Determining the Roche limit

The Roche limit depends on the rigidity of the satellite. At one extreme, a completely rigid satellite will maintain its shape until tidal forces break it apart. At the other extreme, a highly fluid satellite gradually deforms leading to increased tidal forces, causing the satellite to elongate, further compounding the tidal forces and causing it to break apart more readily. Most real satellites are somewhere between these two extremes, with internal friction, viscosity, and tensile strength rendering the satellite neither perfectly rigid nor perfectly fluid.

### Rigid satellites

To calculate the rigid body Roche limit for a spherical satellite, the cause of the rigidity is neglected but the body is assumed to maintain its spherical shape while being held together only by its own self-gravity. Other effects are also neglected, such as tidal deformation of the primary, rotation of the satellite, and its irregular shape. These somewhat unrealistic assumptions greatly simplify the Roche limit calculation.

The Roche limit, , for a rigid spherical satellite orbiting a spherical primary is

,

where is the radius of the primary, is the density of the primary, and is the density of the satellite.

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 on the surface of the satellite closest to the primary. There are two forces on this mass : 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.

Derivation of the Roche limit
The gravitational pull on the mass towards the satellite with mass and radius can be expressed according to Newton's law of gravitation.

The tidal force on the mass towards the primary with radius and a distance between the centers of the two bodies can be expressed approximately as

.

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

or

,

which quickly gives the Roche limit, , as

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 can be written as

where is the radius of the primary.

And likewise

where is the radius of the satellite.

Substituting for the masses in the equation for the Roche limit, and cancelling out gives

,

which can be simplified to the Roche limit:

.

### Fluid satellites

A more accurate 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 it (into a prolate spheroid).

The calculation is complex and its result cannot be represented as an algebraic formula. Historically, Roche himself derived the following numerical solution for the Roche Limit:

However, a better approximation that takes into account the primary's oblateness and the satellite's mass is:

where is the oblateness of the primary. The numerical factor is calculated with the aid of a computer.

The fluid solution is appropriate for bodies that are only loosely held together, such as a comet. For instance, 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. Shoemaker-Levy 9 was first observed in 1993, but its orbit indicated that it had been captured by Jupiter a few decades prior. [1]

#### Derivation of the formula

As the fluid satellite case is more delicate than the rigid one, the satellite is described with some simplifying assumptions. First, assume the object consists of incompressible fluid that has constant density and volume that do not depend on external or internal forces.

Second, assume the satellite moves in a circular orbit and it remains in synchronous rotation. This means that the angular speed at which it rotates around its center of mass is the same as the angular speed at which it moves around the overall system barycenter.

The angular speed is given by Kepler's third law:

The synchronous rotation implies that the liquid does not move and the problem can be regarded as a static one. Therefore, the viscosity and friction of the liquid in this model do not play a role, since these quantities would play a role only for a moving fluid.

Given these assumptions, the following forces should be taken into account:
• The force of gravitation due to the main body;
• the centrifugal force in the rotary reference system; and
• the self-gravitation field of the satellite.
Since all of these forces are conservative, they can be expressed by means of a potential. Moreover, the surface of the satellite is an equipotential one. Otherwise, the differences of potential would give rise to forces and movement of some parts of the liquid at the surface, which contradicts the static model assumption. Given the distance from the main body, our problem is to determine the form of the surface that satisfies the equipotential condition.

As the orbit has been assumed circular, we know that the total gravitational force and centrifugal force acting on the main body cancel. Therefore, the force that affects the particles of the liquid is the tidal force, which depends on the position with respect to the center of mass (already considered in the rigid model). For small bodies, the distance of the liquid particles from the center of the body is small in relation to the distance d to the main body. Thus the tidal force can be linearized, resulting in the same formula for FT as given above. While this force in the rigid model depends only on the radius r of the satellite, in the fluid case we need to consider all the points on the surface and the tidal force depends on the distance Δd from the center of mass to a given particle projected on the line joining the satellite and the main body. We call Δd the radial distance (see the picture). Since the tidal force is linear in Δd, the related potential is proportional to the square of the variable and for we have

We want to determine the shape of the satellite for which the sum of the self-gravitation potential and is constant on the surface of the body. In general, such a problem is very difficult to solve, but in this particular case, it can be solved by a skillful guess due to the square dependence of the tidal potential on the radial distance Δd

Since the potential VT changes only in one direction (i.e. the direction to the main body), the satellite can be expected to take an axially symmetric form. More precisely, we may assume that it takes a form of a solid of revolution. The self-potential on the surface of such a solid of revolution can only depend on the radial distance to the center of mass. Indeed, the intersection of the satellite and a plane perpendicular to the line joining the bodies is a disc whose boundary by our assumptions is a circle of constant potential. Should the difference between the self-gravitation potential and VT be constant, both potentials must depend in the same way on Δd. In other words, the self-potential has to be proportional to the square of Δd. Then it can be shown that the equipotential solution is an ellipsoid of revolution. Given a constant density and volume the self-potential of such body depends only on the eccentricity ε of the ellipsoid:

where is the constant self-potential on the intersection of the circular edge of the body and the central symmetry plane given by the equation Δd=0.

The dimensionless function f is to be determined from the accurate solution for the potential of the ellipsoid

and, surprisingly enough, does not depend on the volume of the satellite.

The graph of the dimensionless function f which indicates how the strength of the tidal potential depends on the eccentricity ε of the ellipsoid

Although the explicit form of the function f looks complicated, it is clear that we may and do choose the value of ε so that the potential VT is equal to VS plus a constant independent of the variable Δd. By inspection, this occurs when

This equation can easily be solved numerically. The graph indicates that there are two solutions and thus the smaller one represents the stable equilibrium form (the ellipsoid with the smaller eccentricity). This solution determines the (eccentricity of) the tidal ellipsoid as a function of the distance to the main body. The derivative of the function f has a zero where the maximal eccentricity is attained. This corresponds to the Roche limit.
The derivative of f determines the maximal eccentricity. This gives the Roche limit.

More precisely, the Roche limit is determined by the fact that the function f, which can be regarded as a (nonlinear) measure of the force squeezing the ellipsoid towards a spherical shape, is bounded so that there is an eccentricity at which this contracting force becomes maximal. Since the tidal force increases when the satellite approaches the main body, it is clear that there is a critical distance at which the ellipsoid is torn up.

The maximal eccentricity can be calculated numerically as the zero of the derivative of f' (see the diagram). One obtains

which corresponds to the ratio of the ellipsoid axes 1:1.95. Inserting this into the formula for the function f one can determine the minimal distance at which the ellipsoid exists. This is the Roche limit,

## 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)
Sun1408696,000,000
Jupiter132671,492,000
Earth55136,378,137
Moon33461,738,100
Saturn687.360,268,000
Uranus131825,559,000
Neptune163824,764,000

Using these data, the Roche Limits for rigid and fluid bodies can easily be calculated. The average density of comets is taken to be 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 density and rigidity.

Body Satellite Roche limit (rigid) Roche limit (fluid)
Distance (km) R Distance (km) R
EarthMoon9,4961.4918,2612.86
Earthaverage Comet17,8802.8034,3905.39
SunEarth554,4000.801,066,3001.53
SunJupiter890,7001.281,713,0002.46
SunMoon655,3000.941,260,3001.81
Sunaverage Comet1,234,0001.782,374,0003.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.

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. Both rigid and fluid body calculations are given. Note Pan and Naiad in particular, which may be quite close to their actual break-up points.

In practice, the densities of most of the inner satellites of giant planets are not known. In these cases (shown in italics), likely values have been assumed, but their actual Roche limit can vary from the value shown.

Primary Satellite Orbital Radius vs. Roche limit
(rigid) (fluid)
SunMercury104:154:1
EarthMoon41:121:1
MarsPhobos172%89%
Deimos451%234%
JupiterMetis~186%~94%
Amalthea175%88%
Thebe254%128%
SaturnPan142%70%
Atlas156%78%
Prometheus162%80%
Pandora167%83%
Epimetheus200%99%
Janus195%97%
UranusCordelia~154%~79%
Ophelia~166%~86%
Bianca~183%~94%
Cressida~191%~98%
Desdemona~194%~100%
Juliet~199%~102%
Thalassa~145%~75%
Despina~152%~78%
Galatea153%79%
Larissa~218%~113%
PlutoCharon12.5:16.5:1

## References

1. ^ Eric W. Weisstein (2007). Eric Weisstein's World of Physics - Roche Limit (English). scienceworld.wolfram.com. Retrieved on September 5, 2007.
2. ^ NASA. What is the Roche limit? (English). NASA - JPL. Retrieved on September 5, 2007.

## Sources

• Édouard Roche: La figure d'une masse fluide soumise à l'attraction d'un point éloigné, Acad. des sciences de Montpellier, Vol. 1 (1847–50) p. 243

Gravitation is a natural phenomenon by which all objects with mass attract each other. In everyday life, gravitation is most familiar as the agency that endows objects with weight.
tidal force is a secondary effect of the force of gravity and is responsible for the tides. It arises because the gravitational field is not constant across a body's diameter.
ORBit is a CORBA compliant Object Request Broker (ORB). The current version is called ORBit2 and is compliant with CORBA version 2.4. It is developed under the GPL license and is used as middleware for the GNOME project.
For other uses, see Coalescence (disambiguation).

Coalescence is the process by which two or more droplets or particles merge during contact to form a single daughter droplet (or bubble).
Édouard Albert Roche (1820-1883) was a French scientist, who is best known for his work in the field of celestial mechanics. He gave his name to the concepts of the Roche sphere, Roche limit and Roche lobe.
Motto
Liberté, Égalité, Fraternité
"Liberty, Equality, Fraternity"
Anthem
"La Marseillaise"

Astronomy is the scientific study of celestial objects (such as stars, planets, comets, and galaxies) and phenomena that originate outside the Earth's atmosphere (such as the cosmic background radiation).
18th century - 19th century - 20th century
1810s  1820s  1830s  - 1840s -  1850s  1860s  1870s
1845 1846 1847 - 1848 - 1849 1850 1851

:
Subjects:     Archaeology - Architecture -
satellite is an object which has been placed into orbit by human endeavor. Such objects are sometimes called artificial satellites to distinguish them from natural satellites such as the Moon.
ORBit is a CORBA compliant Object Request Broker (ORB). The current version is called ORBit2 and is compliant with CORBA version 2.4. It is developed under the GPL license and is used as middleware for the GNOME project.
A natural satellite is an object that orbits a planet or other body larger than itself and which is not man-made. Such objects are often called moons. Technically, the term could also refer to a planet orbiting a star, or even to a star orbiting a galactic center, but these
satellite is an object which has been placed into orbit by human endeavor. Such objects are sometimes called artificial satellites to distinguish them from natural satellites such as the Moon.
Jupiter

This processed color image of Jupiter was produced in 1990 by the U.S. Geological Survey from a Voyager image captured in 1979. The colors have been enhanced to bring out detail.
Metis

Image of Metis was taken by Galileo's solid state imaging system between November 1996 and June 1997.
Discovery
Discovered by: S.
Saturn

Saturn, as seen by Cassini
Orbital characteristics[1][2]
Epoch J2000
Aphelion distance: 1,513,325,783 km
10.11595804 AU
Perihelion distance: 1,353,572,956 km
9.
Pan

Pan amid the rings of Saturn. The 'side' view gives Pan the appearance of being embedded in the rings, although it actually travels in the empty Encke Division.
Discovery
Discovered by: M. R.
Tensile strength , or measures the force required to pull something such as rope, wire, or a structural beam to the point where it breaks.

## Explanation

comet is a small body in the solar system that orbits the Sun and (at least occasionally) exhibits a coma (or atmosphere) and/or a tail — both primarily from the effects of solar radiation upon the comet's nucleus, which itself is a minor body composed of rock, dust, and
A planetary ring is a ring of dust and other small particles orbiting around a planet in a flat disc-shaped region. The most spectacular and famous planetary rings are those around Saturn, but the other three gas giants of the solar system (Jupiter, Uranus, Neptune) possess ring
accretion disc (or accretion disk) is a structure formed by diffuse material in orbital motion around a central body. The central body is typically either a young star, a protostar, a white dwarf, a neutron star, or a black hole.
Roche lobe is the region of space around a star in a binary system within which orbiting material is gravitationally bound to that star. If the star expands past its Roche lobe, then the material outside of the lobe will fall into the other star.
A Hill sphere approximates the gravitational sphere of influence of one astronomical body in the face of perturbations from another heavier body around which it orbits. It was defined by the American astronomer George William Hill based upon the work of the French astronomer
Gravitation is a natural phenomenon by which all objects with mass attract each other. In everyday life, gravitation is most familiar as the agency that endows objects with weight.
Astronomical objects are significant physical entities, associations or structures which current science has confirmed to exist in space. This does not necessarily mean that more current science will not disprove their existence.
ORBit is a CORBA compliant Object Request Broker (ORB). The current version is called ORBit2 and is compliant with CORBA version 2.4. It is developed under the GPL license and is used as middleware for the GNOME project.