Planetary science

Planetary science studies planets, moons, small bodies, planetary systems, and the processes that form and evolve them.

A magnetosphere is the region around a body in which the body's magnetic interaction controls charged-particle motion and its response to the surrounding plasma. A planet with an intrinsic magnetic field can support a compressed dayside and an extended magnetotail in the solar wind. Magnetospheric currents must be distinguished from the internal field when inferring a planetary dynamo.

An induced magnetosphere is maintained by the interaction of an external magnetized plasma flow with a conducting ionosphere or atmosphere. Its currents and draped magnetic field do not establish an internal planetary dynamo. Venus is a central example; Mars combines this interaction with localized crustal fields.

A planetary magnetic field may include an intrinsic field maintained by an internal planetary dynamo, local crustal remanent magnetization, and externally driven fields associated with the solar wind. Separating these components requires spatial and temporal observations. A global internally sourced field, a localized crustal anomaly, and an induced magnetosphere are different evidence about the planet's present and past interior.

A planetary dynamo maintains an intrinsic magnetic field by motion of an electrically conducting fluid. The resistive induction equation balances stretching and transport against magnetic diffusion. A sufficiently large magnetic Reynolds number is necessary but does not ensure that the flow regenerates every required field component. Thermal or compositional convection supplies mechanical power; rotation organizes the flow through the Coriolis force, and the Lorentz force density eventually modifies the flow and limits growth. A fluid core can fail to operate a dynamo if it is stably stratified or inadequately driven.

The geodynamo is the maintenance of Earth's magnetic field by electrically conducting fluid motion in its liquid outer core. Rotation and buoyant convection organize helical flow and permit mean-field poloidal–toroidal regeneration. A scalar homogeneous alpha model explains a mechanism and a threshold but omits spherical boundaries, anisotropy, spatial variations and Lorentz force feedback. Reversals and saturation require dynamical or time-dependent transport information beyond a growing constant-coefficient kinematic mode.

A geomagnetic reversal changes the large-scale polarity of Earth's intrinsic magnetic field. Crustal remanent magnetization records past reversals. The reversal times are not a regular clock; their occurrence reflects time-dependent dynamo action rather than rigid rotation of a permanent magnet.

An asteroid is a small rocky or metallic body orbiting a star. Asteroid dynamics can involve mean-motion resonances, planetary scattering and secular perturbations. Solar-system asteroids illustrate how phase protection can prevent close encounters even when orbital paths approach a planet.

The astronomical unit is a length scale used for Solar System distances, approximately $1.5\times 10^{13}\mathrm{cm}$. Coarse estimates sometimes round it to $10^{13}\mathrm{cm}$. Numerical coefficients in planet formation estimates depend on which approximation is used.

Core accretion forms a planet by building a solid core through growth and accretion of solid material. A sufficiently massive core can then bind a gaseous envelope. It is distinct from direct fragmentation through gravitational instability of an astrophysical disk.

The minimum-mass solar nebula estimates a primordial gas-and-solid protoplanetary disk by spreading planetary material into annuli and restoring a solar mixture. A common idealized surface density is proportional to $r^{-3/2}$. This reconstruction is a model of minimum material requirements, not a unique history of the Solar System. Extrapolating it far beyond the planetary region need not give a physically consistent disk.

Accretion releases gravitational energy. Efficient radiation at an accretion shock leaves a lower-entropy cold start, while greater retained energy gives a higher-entropy hot start; intermediate outcomes are possible. Initial entropy changes a young planet's radius, luminosity and spectrum. Consequently an exoplanet direct imaging luminosity and age need not uniquely fix mass without a formation/evolution model.

A comet tail is an extended structure formed from material released by a comet. Its dust tail is shaped largely by radiation pressure, while an ion tail is shaped by the solar wind and magnetic field. The cometary dust tail is the solid-grain component.

A cometary dust tail is the extended distribution of solid grains released by a comet. Radiation pressure sorts its grains by size: high-$\beta$ grains are displaced rapidly in the anti-stellar direction, while larger low-$\beta$ grains remain closer to the projected orbit and form a more curved trail.

With uniformly distributed apsidal orientations, the phase-mixed radial probability of a Kepler orbit gives the surface probability per cometFor $Q\gg r>q$, this is ${\pi }^{-2}{Q}^{-3/2}{r}^{-1/2}(1-q/r{)}^{-1/2}$. Multiply by the population size for a number surface density of a disk.

The innermost planet has a large metallic core and a weak but global intrinsic magnetic field. Its observed field supports an active planetary dynamo in a fluid conducting part of that core. Slow rotation does not by itself exclude dynamo action.

A planet dominated by rocky and metallic material rather than a massive H/He envelope. Earth and Mars are examples. A super-Earth is a mass category and need not be a terrestrial planet.

A giant planet with a much larger heavy-element fraction than a typical gas giant, exemplified by Uranus and Neptune. Water, ammonia and other volatile-rich materials are called “ices” in formation terminology; they need not remain ordinary frozen solids in the hot deep interior.

A large planet whose bulk mass includes a substantial hydrogen–helium envelope, as in Jupiter and Saturn. The deep envelope is fluid rather than a dilute atmosphere with a sharp solid surface. Its cooling, compression and possible electron degeneracy pressure influence the mass-radius relation.

Earth's outer core is a liquid, electrically conducting region whose buoyant motion and rotation sustain the geodynamo. Its magnetic field both responds to the flow through induction and modifies that flow through the Lorentz force. A prescribed-velocity kinematic model describes onset rather than the entire nonlinear core dynamics.

The terrestrial atmosphere has regions with opposite signs of the atmospheric lapse rate: the troposphere usually cools upward, the stratosphere warms upward through ultraviolet heating, the mesosphere cools, and the thermosphere warms. A monotone atmospheric pressure-temperature profile can approximate one such interval, but cannot describe all their transitions with one fixed branch.

The thermosphere is a high atmospheric region where absorption of energetic stellar radiation can produce an upward-rising kinetic temperature. Its very low density distinguishes kinetic temperature from a large thermal energy content.

The mesosphere lies above the stratosphere and below the thermosphere. Its typical terrestrial temperature decreases upward.

The terrestrial stratosphere lies above the troposphere and has an upward-warming region associated with ultraviolet heating by the ozone layer.

The troposphere is the lowest main layer of the terrestrial atmosphere, with weather and generally decreasing temperature with altitude. Local inversions and moist convection complicate a single lapse-rate description.

A rogue planet has no bound host star. Without stellar irradiation, its luminosity can still be supplied by internal cooling and Kelvin-Helmholtz contraction; its emitting temperature need not equal the negligible irradiation-based planetary equilibrium temperature.

A natural satellite is a naturally occurring body orbiting a planet or another nonstellar body. The Moon is Earth's natural satellite; its synchronous rotation provides a familiar example of tidal locking.

The Moon is Earth's natural satellite. Its synchronous rotation keeps approximately the same near side of the Moon facing Earth, while both the near side of the Moon and the far side of the Moon experience sunlight and darkness. See illumination of a synchronously rotating moon for the geometry. ^[ref].

The near side of the Moon is the hemisphere directed approximately toward Earth. Libration changes the exact visible boundary, so synchronous rotation does not mean that an observer always sees precisely the same hemisphere.

The far side of the Moon is the hemisphere directed approximately away from Earth. “Far” describes the direction to Earth; it does not describe illumination by the Sun. The far side of the Moon experiences both sunlight and darkness during the lunar phase cycle. ^[ref].

A lunar phase is the appearance of the sunlit portion of the Moon as seen from Earth. Outside a lunar eclipse, the distant Sun illuminates approximately one hemisphere of the Moon, but an observer sees a changing overlap of that hemisphere with the hemisphere facing Earth. The lunar phase therefore changes even though synchronous rotation keeps approximately the same near side of the Moon facing Earth. ^[ref].

The cycle runs through new Moon, waxing crescent, first quarter, waxing gibbous, full Moon, waning gibbous, last quarter and waning crescent. “Waxing” means that the lunar illuminated fraction is increasing; “waning” means it is decreasing. A crescent has less than half its apparent disk sunlit, and a gibbous lunar phase has more than half. First and last quarter each show half a disk: “quarter” refers to progress around the cycle, rather than a quarter of the disk being lit. ^[ref].

An instructive circular orbit model puts Earth at the origin, the Sun far away along $+x$, and the Moon at $R(\cos \theta ,\sin \theta )$. The sunward position $x>0$, $\theta =0$, gives new Moon; the opposite position $x<0$, $\theta =\pi$, gives full Moon. Counterclockwise motion gives first quarter at $\theta =\pi /2$ and last quarter at $\theta =3\pi /2$. The Sun's direction is fixed during this snapshot of the cycle. Every small globe in the orbital view is illuminated on its $+x$ half; the adjacent apparent disks are separate views from Earth, rather than globes viewed from above.

The disk views choose north up, with the upward direction perpendicular to the orbital plane. In this convention waxing lunar phases are bright on the right and waning lunar phases on the left. Local horizon orientation and viewing from different locations on Earth can rotate the apparent disk; its lit edge always points toward the projected direction of the Sun. The lunar illuminated fraction is unaffected by this rotation.

The useful angle is the phase angle $\alpha =\angle SME$, with $S$ the Sun, $M$ the Moon and $E$ the observer at Earth's center: its vertex is at the Moon. Full Moon has $\alpha =0$, quarter has $\alpha =\pi /2$, and new Moon has $\alpha =\pi$ in the ideal model. By contrast, the elongation $ϵ=\angle SEM$ has its vertex at Earth. For a distant Sun, its rays are parallel, soHere $k$ is the lunar illuminated fraction: a fraction of projected disk area, not a fraction of the Moon's total surface area or its received radiant flux relative to full Moon. The phase angle and unsigned elongation each lie between $0$ and $\pi$, so either alone gives the illuminated amount without distinguishing waxing from waning. The orbital direction resolves that ambiguity.

The lunar illuminated fraction has an elegant geometric derivation. Model the Moon as a sphere of unit radius and view it along the $z$ direction. Take the Sun's unit direction to be $(\sin \alpha ,0,\cos \alpha )$. A visible surface point above apparent disk coordinates $(x,y)$ has $z=\sqrt{1-x^2-y^2}$. Direct illumination requires the outward normal's dot product with the sunward direction to be positive. The day-night terminator therefore obeys $x\sin \alpha +z\cos \alpha =0$. At a fixed $y$, set $q=\sqrt{1-y^2}$; the visible day-night terminator is $x_t=-q\cos \alpha$, and the illuminated strip extends from $x_t$ to the right edge $q$. Integrating its width givesThe integral is the area of half a unit circle. Thus a three-dimensional visibility problem reduces to a circle's area and a cosine. The derivation assumes parallel sunlight, a distant observer and a sphere, and ignores a lunar eclipse. Real observed brightness additionally depends on surface scattering of light and observing geometry, so $k=1/2$ does not promise half the full-Moon radiant flux.

For a worked example, elongation $ϵ=60^{\circ }$ gives phase angle $\alpha =120^{\circ }$ and $k=(1+\cos 120^{\circ })/2=1/4$: a crescent with $25\%$ of its apparent disk sunlit. At elongation $120^{\circ }$, the phase angle is $60^{\circ }$ and $k=3/4$: a gibbous lunar phase. At $90^{\circ }$, $k=1/2$. In each case the same $k$ occurs once while waxing and once while waning in the ideal cycle.

Ordinary lunar phases are not caused by Earth's shadow. A lunar eclipse occurs only when the Moon actually enters that shadow near full Moon. The orbital inclination of the Moon is about $5^{\circ }$ relative to Earth's orbital plane around the Sun; most full Moons pass above or below the shadow. This is why full Moon does not imply a lunar eclipse every month. ^[ref].

The phase cycle is a synodic month, about $29.5$ days, whereas the Moon's orbital period relative to distant stars is a sidereal month, about $27.3$ days. During one sidereal month, Earth advances around the Sun, so the Moon must move farther to repeat the sunward alignment. For coplanar circular orbits traversed in the same direction, subtract their mean motions: $n_{\mathrm{syn}}=n_{\mathrm{sid}}-n_{\mathrm{year}}$, with $n=2\pi /P$. HenceUsing $P_{\mathrm{sid}}=27.321661\mathrm{days}$ and $P_{\mathrm{year}}=365.256\mathrm{days}$ gives $P_{\mathrm{syn}}=29.530591\mathrm{days}$, consistent with the mean synodic month of $29.53059$ days. After one sidereal month, the sunward direction has advanced about $26.93^{\circ }$; catching up takes about $2.21$ more days. These are mean-period estimates: orbital eccentricity and changing geometry make individual synodic months vary. ^[ref]. ^[ref].

The lunar phase is also a practical clock for observing. Full Moon is opposite the Sun and rises roughly at sunset; first quarter is about a right angle from the Sun and is convenient in the evening, while last quarter is convenient before dawn. These are rough guides, not exact rise-time predictions for every place and season. ^[ref]. Historically, repeated lunar phases provided the month in lunar calendars, while half-lit disks gave Aristarchus of Samos a way to infer astronomical distances. ^[ref].

The astronomical phase angle is the angle at an illuminated body between the directions to its light source and to its observer. For a lunar phase it is the Sun–Moon–observer angle at the Moon, rather than the elongation measured at Earth. A small phase angle gives a nearly full apparent disk in the spherical model.

Elongation is the apparent angular separation of a celestial body from the Sun, measured at the observer. For a lunar phase, the elongation is the Sun–Earth–Moon angle; in the distant-Sun approximation it is supplementary to the phase angle. A given unsigned elongation can occur on either side of the Sun in the sky.

The lunar illuminated fraction is the directly sunlit projected area of the apparent lunar disk divided by its total projected area. For a sphere with parallel sunlight and a distant observer, outside a lunar eclipse, it is $k=(1+\cos \alpha )/2$, with phase angle $\alpha$ at the Moon. It measures geometry rather than the fraction of full-Moon radiant flux. ^[ref].

A sidereal month is the Moon's orbital period relative to distant stars, about $27.32166$ days. It is shorter than the synodic month because the sunward direction changes as Earth travels around the Sun. The lunar spin period associated with synchronous rotation matches the sidereal month in the mean. ^[ref].

A synodic month is the synodic period of the Moon relative to the Sun, about $29.53059$ days on average. A lunation is commonly measured from one new Moon to the next. Actual lunations vary, so the constant-speed relation derived in lunar phase is useful for understanding the mean cycle rather than predicting exact dates. ^[ref].

Aristarchus of Samos used the half-lit Moon to estimate the relative distances of the Sun and Moon. In the spherical illumination model a half-lit disk means phase angle $\alpha =\pi /2$, so the Sun–Moon–Earth triangle is right-angled at the Moon. Keeping the Sun at a finite distance, let $d_{EM}$ and $d_{ES}$ be the distances from Earth to the Moon and Sun. Then elongation $ϵ$ at Earth givesThe $87^{\circ }$ elongation attributed to Aristarchus of Samos gives a ratio about $19.1$, consistent with his conclusion that the Sun was between $18$ and $20$ times farther away. The example shows how lunar phase geometry becomes a distance measurement, but the ratio is extremely sensitive to an angle close to $90^{\circ }$. Unlike the distant-Sun approximation, this finite triangle has elongation slightly less than $90^{\circ }$ at half illumination. ^[ref].

A lunar eclipse occurs when Earth blocks some or all of the direct sunlight reaching the Moon. It requires alignment near full Moon; the orbital inclination usually carries the Moon clear of Earth's shadow. Ordinary lunar phases instead arise from viewing different amounts of the sunlit hemisphere. ^[ref].

Libration is the apparent rocking of a body as seen by an observer. For the Moon, longitudinal libration arises because orbital eccentricity makes its orbital angular speed vary while its spin is much more uniform: the direction to Earth advances ahead of, then falls behind, the mean facing direction. The Moon's tilted geometry also causes latitudinal libration; small actual rotational oscillations provide further corrections. Thus lunar synchronous rotation holds in the mean, rather than requiring exact instantaneous alignment. ^[ref]. ^[ref].

Planetary system dynamics studies the orbital evolution of planets, planetesimals, and dust under gravity, radiation, collisions, and dissipative forces.

Orbit crossing occurs when two geometrical orbital paths intersect. It permits collisions or close encounters when bodies also reach the intersection at compatible times; crossing alone does not imply a collision. For nearly circular coplanar nested orbits, a useful local test follows by differentiating their radial positions at fixed longitude.

For a cold coplanar disk with complex eccentricity $z(a,t)$, the first-order radial shape is $r=a[1-\mathrm{Re}(z{e}^{-i\theta })]$. Infinitesimally neighboring orbits first cross when $|z+a{\partial }_{a}z|=1$. Under secular forcing by a distant outer planet, $z=Z_f(1-e^{iAt})$, $Z_f\propto a$ and $A\propto a^{3/2}$, so the secular phase gradient gives $t_{\mathrm{cross}}\simeq 2/(3Z_{f}A)$. The finite-neighbor version requires the accumulated phase difference to remain small when replacing a difference by a derivative.

Isotropic stellar mass loss changes the central gravitational parameter $\mu (t)$ while preserving the central direction of the force, provided the escaping material causes no recoil or drag on the orbiting body. The relative equation is $\ddot{\mathbf{r}}=-\mu (t)\mathbf{r}/r^3$. Specific angular momentum remains exactly constant; specific orbital energy obeys $\dot{C}=-\dot{\mu }/r$. Orbital expansion is adiabatic only when the mass-loss law varies slowly on the orbital period.

Time averaging over a Kepler orbit gives $⟨\cos f⟩=-e$, as follows from $dt=(1-e\cos E)dE/n$ and $\cos f=(\cos E-e)/(1-e\cos E)$. Consequently slow isotropic stellar mass loss preserves the secular orbital eccentricity. Conserved specific angular momentum then gives $\mu a=\mathrm{constant}$, so pericentre, apocentre and semimajor axis expand inversely with the remaining gravitational parameter.

The signed instantaneous contribution of apsidal precession to azimuthal motion during isotropic stellar mass loss is $\dot{\varpi }/\dot{\theta }=-(\dot{\mu }/\mu n)(1-e^2{)}^{3/2}\sin f/[e(1+e\cos f{)}^2]$. Its maximum absolute magnitude is small when $n|\mu /\dot{\mu }|\gg F(e)$, where $F(e)\sim 1/e$ near a circular orbit and $F(e)\to 3\sqrt{3}/4$ as $e\to 1$. The low-orbital eccentricity coordinate singularity can be avoided by evolving the eccentricity vector.

For isotropic stellar mass loss the eccentricity vector obeys $\dot{\mathbf{e}}=-(\dot{\mu }/\mu )(\mathbf{e}+\hat{\mathbf{r}})$. With true anomaly $f$ and $e>0$, this gives $\dot{e}=-(\dot{\mu }/\mu )(e+\cos f)$ and $\dot{\varpi }=-(\dot{\mu }/\mu )\sin f/e$. The osculating pericentre satisfies $\dot{q}/q=-(\dot{\mu }/\mu )(1-\cos f)/(1+e)$, hence cannot decrease during mass loss and is instantaneously stationary at pericentre. The vector form is nonsingular at a circular orbit.

Anisotropic thermal radiation produces recoil that changes the semi-major axis of a rotating solid body. Finite thermal inertia provides the temperature lag. Large bodies have small acceleration per unit mass, while very small bodies become nearly isothermal. The Yarkovsky effect can supply bodies to dynamical escape regions.

Drift driven by the Yarkovsky effect can carry a body out of a belt or into an unstable orbital region. A useful two-limit removal law is $t_y=BD+C/D$, with slow escape at both large and small diameters.

A population supplied at differential number rate $n_b(D)/t_y(D)$ and surviving externally for a size-independent time $T$ has $n_{\mathrm{out}}(D)=T{n}_b(D)/t_y(D)$. Its differential number slope changes by the negative of the removal-time exponent, unless the belt is already removal-dominated, in which case it reflects the fragment injection spectrum.

For $t_c=(A/M)D^{1/2}$, the ratio $t_y/t_c$ has its minimum at $D_{\ast }=\sqrt{3C/B}$. Removal beats collisions somewhere exactly when $M<M_y=(3^{3/4}/4)A{B}^{-3/4}{C}^{-1/4}$, provided the relevant diameters are present. This criterion uses the collision-only cascade as a reference.

The law $t_y=BD+C/D$ for positive $B,C$ has a minimum at $D_0=\sqrt{C/B}$ with value $2\sqrt{BC}$. This size differs from the minimum of the ratio of removal and collision times.

A passing star perturbs a planetary or cometary system through its differential Newtonian gravitational field. A distant encounter is predominantly tidal; a sufficiently close one can substantially change specific orbital energy and angular momentum, eject a comet, or shift its periapsis into or out of the planetary region.

Planetary migration changes a planet's orbital semimajor axis through torques from a disc or exchanges of energy and angular momentum with other bodies, often followed by tidal dissipation. Short-period giant planets, eccentric progenitors, resonant architectures, and stellar obliquity provide complementary constraints rather than one unique migration diagnostic.

For a fixed mass on a Keplerian circular orbit, $L_s=M_s\sqrt{G{M}_{\ast }{r}_0}$ has derivative $M_s{r}_0\Omega /2$. The torque on the satellite is minus the torque $\Gamma$ on the disk. This gives the displayed drift law. In the density-slope satellite torque model, ${\dot{r}}_0=-(8\chi /9)\beta q({\Sigma }_0{r}_0^2/M_{\ast })r_0\Omega (H/r_0{)}^{-2}$.

The stellar obliquity is the angle between a star's spin and a planet's orbital angular momentum. The Rossiter-McLaughlin effect constrains the angle $\lambda$ between the sky projections of those two axes, rather than a component of $\psi$. If the orbital and stellar-spin inclinations to the line of sight are $i_{\mathrm{orb}}$ and $i_{\ast }$, then $\cos \psi =\cos i_{\ast }\cos i_{\mathrm{orb}}+\sin i_{\ast }\sin i_{\mathrm{orb}}\cos \lambda$; $\lambda$ alone generally does not determine $\psi$. Large misalignments constrain formation and planetary migration, but can also reflect a tilted birth disc.

A transiting body hides different velocity regions of a rotating stellar disc and causes an anomalous measured Doppler shift. Its time dependence constrains the sky-projected spin-orbit angle, which is not generally the full stellar obliquity.

Tidal locking settles a body's spin into a sustained relation with its orbital motion through tidal forces and tidal dissipation. The simplest case is synchronous rotation, a $1:1$ spin-orbit resonance. On an ideal circular orbit, with spin axis perpendicular to the orbital plane, the body keeps the same face toward its companion, as the Moon approximately does toward Earth. ^[ref].

This article retains the broader sense of tidal locking as capture into a spin-orbit resonance; the term is also commonly used specifically for synchronous rotation. Tidally evolved rotation need not be $1:1$: Mercury has a $3:2$ spin-orbit resonance with the Sun, turning three times in two orbital periods. The distinction matters when inferring illumination or climate from a spin period. ^[ref].

Tidal deformation is the change in a body's shape caused by differential gravity across it, that is, by a tidal force. The familiar two-bulge picture represents an elongated response along the companion direction, rather than a separate force pulling only the nearest surface. If spin and orbital motion differ, the forcing sweeps through the body's material. With tidal dissipation, the response is out of alignment with the instantaneous forcing, allowing a net torque. ^[ref].

A secular tidal torque is the torque remaining after averaging the changing tidal interaction over fast orbital and rotational motion. The misaligned tidal deformation lets the companion's gravity change the spin angular momentum; the orbit receives the opposite torque. Tidal dissipation converts mechanical energy into heat, while total angular momentum is conserved in an isolated system. Thus the spin-orbit coupling exchanges angular momentum rather than destroying it. ^[ref].

For intuition near synchronous rotation on a circular orbit, hold mean motion $n$ and moment of inertia $I$ fixed and adopt the illustrative linear dissipative lawHere $\omega$ is spin angular velocity, $\theta$ is a marked surface direction in an inertial frame, and $K$ is a response coefficient, not a universal material constant. A faster spin is slowed and a slower prograde spin is accelerated. The solution and the associated mechanical power loss areThe orbital contribution to this power is $-n⟨\tau ⟩$, and the spin contribution is $\omega ⟨\tau ⟩$; their sum becomes heat. Holding $n$ fixed is a local approximation to the coupled evolution. Actual secular tidal torque depends on the deformation response and orbital eccentricity, so this simple law does not establish which spin-orbit resonance an eccentric body captures into.

A prograde spin-orbit resonance relates mean spin angular velocity to mean motion by a ratio of small positive integers, $⟨\omega ⟩/n=p/q$: the body makes $p$ turns in $q$ orbital periods. Synchronous rotation is $1:1$; Mercury's $3:2$ spin-orbit resonance is nonsynchronous. The relevant combination of spin and orbital angles must stay bounded around a preferred orientation, not merely have a near-integer frequency ratio. For an elongated body on an eccentric orbit, the companion's gravitational torque on the permanent shape can sustain the resonance while tidal dissipation permits capture. ^[ref].

Synchronous rotation means one prograde turn about the body's own axis per orbital period, measured against fixed directions in an inertial frame. This $1:1$ spin-orbit resonance approximately describes the Moon's synchronous rotation toward Earth. ^[ref].

The striking geometric idea is that one turn in space can keep one face toward the companion: a fixed orientation in the frame turning with the orbit is a changing orientation in an inertial frame. This distinction helps us predict illumination: a synchronously rotating exoplanet orbiting its illuminating star can have persistent day and night hemispheres, whereas the Moon, synchronized to Earth, receives changing illumination from a different body, the Sun.

View a circular orbit along its perpendicular spin axis. Let $f$ be the orbital angle measured from a fixed direction, choosing an arbitrary zero for the circular limit of true anomaly, and let $\theta$ be the direction of an outward arrow painted on the body, measured from the same fixed direction. The direction from the body to its companion is $f+\pi$. Thus an inward-facing arrow satisfiesMore generally synchronous rotation keeps $\theta -f$ constant; its particular value selects which surface direction faces the companion. Over one orbital period $T$, both $f$ and the unwrapped $\theta$ increase by $2\pi$. This is the short geometric explanation for one turn per orbit. For a nonspinning body, $\dot{\theta }=0$ instead, and the companion direction sweeps through all surface longitudes. On an eccentric orbit, mean $1:1$ spin does not force $\theta -f$ to stay exactly constant: true anomaly advances nonuniformly, producing the longitudinal component of libration.

Facing a companion and facing the source of light are different relations. The Moon keeps its near side of the Moon approximately toward Earth, while the direction to the distant Sun is nearly fixed during the idealized orbit shown below. A fixed surface patch therefore turns toward, across and away from sunlight. The far side of the Moon has day and night too; its name does not mean permanent darkness. ^[ref].

If the Sun direction has fixed angle $\alpha$ and the marked surface patch has outward normal vector at angle $\theta$, its direct illumination is proportional to $\max (0,\cos (\theta -\alpha ))$. With synchronous rotation, $\theta =f+\pi$, so the patch goes through day and night even though it always faces Earth. In reality the Earth-Moon system advances around the Sun: the lunar spin and orbital cycle relative to distant stars is about $27.3$ Earth days, while the lunar phase cycle has a synodic period of about $29.5$ Earth days. These numbers describe the actual system, not the fixed-sunlight cartoon. ^[ref].

For a synchronously rotating exoplanet whose orbital companion is itself the illuminating star, the direction to the star turns with the orbit. On an ideal circular orbit, with perpendicular spin axis and no additional oscillations, this direction is fixed relative to the surface, giving persistent day and night hemispheres. Day-night heat redistribution then helps determine the nightside temperature; permanent lack of direct starlight does not imply zero temperature. Orbital eccentricity and departures from this ideal spin geometry modify that simple picture.

Tidal heating converts mechanical energy of tidal deformation into internal heat through tidal dissipation. Sustained orbital eccentricity or a spin axis tilted relative to the orbital plane's perpendicular can maintain changing tidal deformation and heat an exoplanet interior. Thus synchronous rotation alone does not eliminate tidal heating.

Tidal dissipation converts mechanical energy associated with a time-varying tidal deformation into heat. It can reduce orbital eccentricity, drive synchronous rotation, and contribute deep planetary heating. A circular orbit with exact synchronous rotation, perpendicular spin axis and no additional forcing has a static tidal deformation in the body's frame and need not sustain substantial tidal heating.

An astronomical conjunction occurs when two orbiting bodies have the same apparent longitude as viewed from their primary or from an observer.

For two uniform angular motions with $P_b<P_c$, an inner exoplanet transit at zero and next outer exoplanet transit at $\tau$, the next conjunction solves $(n_b-n_c)t+n_c\tau =2\pi$. Hence its longitude modulo $2\pi$ isFor an eccentric outer body this is only an approximation. Let $f_L,M_L$ be its line-of-sight anomalies and ${\Delta }_e(t)=f(t)-M(t)-(f_L-M_L)$. Exact true alignment instead requires $(n_b-n_c)t+n_c\tau -{\Delta }_e(t)=2\pi$. The correction follows directly from uniform mean anomaly and nonuniform true anomaly.

Keplerian shear is differential orbital motion caused by $n\propto a^{-3/2}$: neighboring material at smaller semi-major axis advances in longitude relative to material farther out.

Phase mixing disperses an initially localized ensemble because its members have slightly different orbital frequencies, causing their phases to spread with time.

A planetesimal is a solid body in a young or evolved planetary system whose collisions and gravitational interactions help build or erode larger bodies.

For an impactor number density per size $n(D_{\mathrm{im}})$, relative speed $v_{\mathrm{rel}}$, and collision cross-section $\sigma (D,D_{\mathrm{im}})$, the differential collision rate of one target is $n(D_{\mathrm{im}})v_{\mathrm{rel}}\sigma (D,D_{\mathrm{im}})d{D}_{\mathrm{im}}$.

The specific impact energy is an impactor's kinetic energy divided by the target mass. For equal-density spherical bodies,

A binary planetesimal is a gravitationally bound pair of planetesimals. Its internal two-body problem and its centre of mass orbit respond differently to external forces, including gas drag.

For linear gas drag with aerodynamic stopping times $t_1,t_2$, masses $m_1,m_2$, total mass $m=m_1+m_2$, and relative velocity $\mathbf{v}$, defineIf $\mathbf{w}$ is the gas velocity relative to the centre of mass, the relative drag acceleration is $-\kappa \mathbf{v}+\Delta \mathbf{w}$. The first term dissipates internal specific orbital energy; the second is a differential headwind that vanishes for equal aerodynamic stopping times.

For a weakly perturbed circular binary planetesimal, the constant differential headwind does zero work averaged over the unperturbed orbit. Orbit-averaged drag work givesWith constant aerodynamic stopping times, $a(t)=a_0{e}^{-t/\tau }$. This averaging assumes drag is slow compared with the orbital period and does not appreciably distort the circular orbit within one period.

For a nearly circular stellar orbit of radius $a_b$ in gas moving at $(1-\eta )$ times the local circular speed, the averaged tangential gas drag on the centre of mass is $-\eta {\kappa }_b{v}_K$. Orbit-averaged drag work givesIf $m_2/m_1=ϵ\ll 1$ and aerodynamic stopping time scales as $m^k$, $0<k<1$, then $\tau /{\tau }_b\simeq \eta {ϵ}^k\ll 1$. Internal contraction can outpace stellar migration. Tightly gas-coupled drift requires solving radial and azimuthal motion together instead of this weak-drag approximation.

For a small body with orbital elements $a,e,I$ relative to a planet on a circular orbit of semi-major axis $a_p$, the Tisserand parameter is approximatelyIt is derived from the Jacobi constant and is approximately conserved across separated weak encounters in the circular restricted three-body problem.

For bound planet-crossing Kepler orbits, $T=2/(q+Q)+2\sqrt{2qQ/(q+Q)}\cos I$, with the planet radius taken as one. At fixed $q<1$, $T$ is maximized by $I=0$ and $Q=1$. For $2<T<3$ this gives $q\ge (1-\sqrt{3-T}{)}^2/[1+2\sqrt{3-T}-(3-T)]$. This form remains regular at $T=\sqrt{8}$. It is an invariant-based accessibility bound, not a guarantee that a scattering history attains it; for $0<T\le 2$ it gives no positive universal pericentre floor.

In units where a circular planet has orbital radius and speed one, the particle's relative speed $u$ outside a close encounter satisfies $u^2=3-T$, to leading order in the planet-to-star mass ratio. This follows from the vis-viva equation and $u^2=v^2+1-2h_z$ at the encounter radius. A short conservative encounter changes the direction of this relative velocity while preserving its magnitude.

The disturbing function is the perturbing part of the gravitational potential used to derive the evolution of orbital elements in celestial mechanics.

For a particle far inside a planet's stellar orbit, the planet-centred disturbing function is $G{M}_{\star }/|{\mathbf{r}}_{\star }-\mathbf{r}|-G{M}_{\star }\mathbf{r}\cdot {\mathbf{r}}_{\star }/r_{\star }^3$. Expanding the first term cancels its dipole against the indirect disturbing function. Dropping the spatially constant monopole leaves $G{M}_{\star }{r}^2(3{\cos }^2\psi -1)/(2r_{\star }^3)$, where $\psi$ is the angle between the two offset vectors. The omitted instantaneous terms are $O(G{M}_{\star }{r}^3/r_{\star }^4)$. The gradient of the retained expression is the leading differential stellar tidal acceleration.

For two circular Kepler orbits with mutual inclination $I_m$, independently average the stellar quadrupole disturbing function over both mean anomalies. A circular ring with unit normal $\mathbf{k}$ has $⟨{\hat{r}}_i{\hat{r}}_j⟩=({\delta }_{ij}-k_i{k}_j)/2$. Contracting the two ring tensors gives $⟨{\cos }^2\psi ⟩=[1+(\mathbf{k}\cdot {\mathbf{k}}_p{)}^2]/4$, proving the displayed exact angular factor. For small orbital inclinations relative to a common plane, $I_m^2=I^2+I_p^2-2I{I}_p\cos ({\Omega }_p-\Omega )+O(I^4)$. Thus the secular angular factor is $1/4-3(I^2+I_p^2)/8+3I{I}_p\cos ({\Omega }_p-\Omega )/4$. This potential generates the slow stellar contribution to satellite nodal precession.

The Lagrange planetary equations express changes of osculating orbital elements in terms of derivatives of a disturbing function. Their small-eccentricity secular limit becomes a linear equation for complex eccentricity.

If retained Lagrange planetary equations have $\dot{a}=-2(p+q)a{e}^q{C}_r\sin \varphi$ and $\dot{e}=-q{e}^{q-1}{C}_r\sin \varphi$, then $d(e^2)/d\log a=q/(p+q)$. Thus $e^2-[q/(p+q)]\log (a/a_{\mathrm{ref}})$ is constant. A resonantly oscillating orbital eccentricity–semi-major axis trajectory retraces one curve rather than enclosing area in that projection.

For disturbing function $\mathcal{R}=(G{M}_2/a_2)f(\alpha )e\cos \varphi$, $\alpha =a/a_2$, the leading planar Lagrange planetary equations give $\dot{\lambda }=n+n(M_2/M_{\star })eg(\alpha )\cos \varphi$. The derivative with respect to $a$ holds the other osculating orbital elements fixed. If one writes $\lambda =n(t)t+ϵ(t)$, its literal derivative includes $t\dot{n}$; the physical mean-longitude equation cannot be obtained by silently dropping that term. A consistent accumulated-phase parametrization uses $\lambda ={\int }_0^{t}n(s)ds+{ϵ}_{\mathrm{int}}$.

For perturbing acceleration $R\hat{\mathbf{r}}+T\hat{\mathbf{\theta }}$ on an osculating orbital element system with fixed central parameter $\mu$,To prove it, differentiate the specific orbital energy $-\mu /(2a)$, obtaining $\dot{a}=(2a^2/\mu )(R\dot{r}+Th/r)$. Substitute $\dot{r}=\mu e\sin f/h$, $h/r=\mu (1+e\cos f)/h$ and $h^2=\mu a(1-e^2)$. Here $f$ is the true anomaly; a normal perturbing acceleration does no instantaneous work.