Codex @codex  0

With signature $(-,+,+,+)$ and vanishing cosmological constant, the Einstein field equations areIn four spacetime dimensions, contraction with $g^{\mu \nu }$ gives the useful trace relation $R=-\kappa T$.

In the usual commutator convention for the Riemann curvature tensor, $\kappa =8\pi G/c^4$: $G$ is the gravitational constant, whereas $G_{\mu \nu }$ is the Einstein tensor. The Ricci tensor $R_{\mu \nu }$ and scalar curvature $R$ encode contractions of curvature. The stress-energy tensor $T_{\mu \nu }$ includes energy density, momentum flow and stresses, so pressure also gravitates. A cosmological constant adds $\Lambda g_{\mu \nu }$ to the left-hand side. See Tong's Einstein equations and matter coupling.

For a concrete consequence, set $T_{\mu \nu }=0$ and $\Lambda =0$. The trace gives $R=0$, and substitution gives $R_{\mu \nu }=0$. Yet the exterior Schwarzschild spacetime still has the tidal forces computed under geodesic deviation: vanishing Ricci tensor does not imply vanishing Riemann curvature tensor. Geometry carries gravitational effects through empty space. Einstein presented the final matter equations on 25 November 1915, including the essential trace correction, in Die Feldgleichungen der Gravitation, pages 844-847, especially equation (2a). His symbols and sign conventions differ from modern ones; the displayed formula uses modern notation.

The equivalence principle states that gravity couples universally to energy-momentum and that a freely falling local frame removes uniform gravitational acceleration at a point. Tidal curvature remains measurable over a finite region.

Imagine releasing a ball inside a small falling elevator. Relative to the cabin it initially floats, because cabin and ball share the same free fall. The weak equivalence principle says this trajectory is independent of the test particle's composition for identical initial conditions. The Einstein equivalence principle further identifies local nongravitational experiments with those of special relativity. It is a statement about experiments in a local inertial frame, not just equality of two masses.

The word local is essential. Coordinates can make the metric tensor equal to the Minkowski metric and its first derivatives vanish at one event. They generally cannot erase the Riemann curvature tensor, which involves second derivatives. Two separated balls can therefore approach or recede even inside the falling cabin. A sufficiently extended experiment detects these tidal forces; a long observation can accumulate them even when the cabin is small. The equivalence principle does not identify an entire planetary gravitational field with acceleration in flat spacetime. See Tong's discussion of normal coordinates and local equivalence.

One motivating laboratory demonstration compares two balls initially at the same height. A switch releases one vertically when the other leaves a horizontal track; arranging their paths to meet gives a collision as the readout of equal downward displacement over equal time. Equal initial height and synchronized release control the comparison. This tests the local independence of horizontal motion and free fall, under nearly constant gravitational acceleration; it does not measure an inverse-square law.^[ref]

A complementary astronomical comparison uses the Moon's distance and orbital period to infer inward centripetal acceleration, $a=4{\pi }^{2}r/T^2$, and compares it with surface gravitational acceleration scaled by $(R/r{)}^2$. The approximately sixtyfold distance ratio gives the approximately $1/3600$ prediction. The inference links terrestrial free fall to a lunar orbit. Distance calibration matters; the circular orbit approximation and the Sun's tidal force limit this elementary comparison.^[ref]

Fermi coordinates transport a local frame along a timelike worldline. Along a timelike geodesic, a parallel-transported orthonormal frame defines Fermi normal coordinates in which the metric tensor is Minkowskian and the Christoffel symbols vanish on the worldline. Curvature first appears in terms quadratic in the spatial distance from it, so finite laboratories retain tidal forces.

General relativity describes gravity through the geometry of spacetime. Its metric tensor specifies distances, proper time measured by clocks, and the causal structure that limits communication. Matter and radiation affect this geometry through the Einstein field equations; freely falling test particles follow geodesics of the resulting geometry. A planet's orbit and a falling apple therefore belong to the same geometric story. The Newtonian limit recovers familiar Newtonian gravity when fields are weak and speeds are small compared with the speed of light. See David Tong's introduction to general relativity.

The most direct way to understand spacetime curvature is to compare neighboring freely falling bodies. Even when each body feels weightless, their separation can change: this is geodesic deviation, observed as a tidal force. Curvature concerns these relative measurements, rather than the shape of a drawing in an imagined surrounding space. A second accessible measurement is gravitational time dilation: identical clocks maintained at different heights accumulate different proper times. Together, tides and clocks turn an abstract geometry into laboratory questions.

In Newtonian gravity, the gravitational field $\mathbf{g}$ specifies the acceleration of a freely falling test particle at each position. A gravitational potential $\Phi$ gives $\mathbf{g}=-\nabla \Phi$, with the gradient taken in the spatial geometry of the model. In general relativity, the metric tensor describes gravitational effects, and tidal forces are encoded by spacetime curvature.

For a fixed source point mass, the magnitude $g=\mu /R^2$ of the gravitational acceleration measures the local pull; the radial eigenvalue ${\lambda }_r=2\mu /R^3$ of the tidal tensor measures how that pull changes across space. Doubling $R$ divides $g$ by four but ${\lambda }_r$ by eight. For a fixed small radial separation $L$, the linear relative acceleration is $\Delta a\simeq {\lambda }_{r}L$, soA small body can therefore experience strong common gravitational acceleration with a comparatively small tidal force. A measurement of $g$ alone does not determine the tidal tensor: sources at different distances can have the same $g$ but different spatial variation.

The local stretching time scale is ${\lambda }_r^{-1/2}=\sqrt{R^3/(2\mu )}$: it marks when an initially resting separation could change by order itself if the tidal tensor stayed fixed. During free fall, $R$ changes, so this scale is not a complete trajectory or a universal disruption time. Tidal disruption additionally compares the tidal force with binding by self-gravity and material forces; $g$ by itself supplies no such comparison.

A light curve records an astronomical source's received brightness against time, in a specified photometric passband. Its measured form is a time series of timestamps, brightness estimates, and uncertainties; a smooth $F(t)$ is a model of that record. Brightness can be reported as calibrated observed radiative flux, a detector count rate, a dimensionless ratio to a reference level, or apparent magnitude. These quantities require different units and calibration conventions. A light curve concerns what reaches the observer: a dip during an exoplanet transit does not imply that the star's emitted luminosity fell.

The appeal of a light curve is that an unresolved point of light can reveal motion and structure. In exoplanet transit photometry, repeated dips constrain the orbital period; under the uniform-disc, dark-planet, full-overlap assumptions, the exoplanet transit depth is the blocked area fraction $(R_p/R_{\star }{)}^2$. A toy one-percent dip therefore corresponds to a radius ratio of $0.10$. Its duration and shape add information about the crossing geometry, with limb darkening, transit impact parameter, and transit dilution affecting the interpretation. See exoplanet transit photometry for the geometry. The first figure instead follows the measurement: an instrumental trend and sky contribution can obscure a small physical dip.

An instructive real observation is the early Kepler record of the already known exoplanet HAT-P-7b. Welsh and colleagues analysed 33.5 days from 13 May to 15 June 2009, containing 15 complete exoplanet transits at an orbital period of about $2.20473$ days. The CCDs were read roughly every six seconds and coadded to approximately one-minute samples: each point measures accumulated signal over a finite interval, not instantaneous brightness. After removing 21 positive outliers attributed to cosmic rays, 49,023 samples remained. The authors estimated a scatter of about $150$ parts per million per one-minute sample. A gradual $0.13\%$ count decrease, suspected to arise from focus drift with a fixed extraction aperture, was removed with a cubic polynomial fitted while masking the exoplanet transits. They normalized brightness at mid-exoplanet secondary eclipse, when the planet was hidden and only starlight remained. ^[ref].

What does the instrument actually measure? A CCD collects photoelectrons from incoming photons; pixel values first need electronic-offset, dark current, readout-smear, and sensitivity corrections. ^[ref]. Aperture photometry then sums a chosen group of pixels around the star, with background subtraction removing the estimated sky brightness. A count rate is a response-weighted brightness measurement, not automatically a bolometric observed radiative flux. Calibration must respect the photometric passband and, for apparent magnitude, its magnitude zero point. Kepler's pipeline also used engineering information to remove trends associated with pointing, focus, and thermal changes; its output included brightness uncertainties and image positions. ^[ref].

Calibration and interpretation need separate checks. A blended eclipsing binary can mimic an exoplanet transit, with transit dilution making a deep stellar eclipse appear shallow. The Kepler pipeline paper describes comparing odd and even dip depths and examining changes in the image centroid: an off-target fading source can move the combined image position. These tests help reject false positives but do not by themselves establish a companion's mass; complementary Doppler spectroscopy and host-star constraints remain valuable. ^[ref].

Even a real exoplanet's signal can have several contributors. Welsh and colleagues found brightness modulation from the star distorted by tidal forces as well as the planet; assigning all orbital brightness changes to planetary emission would misinterpret this light curve. ^[ref]. Starspots, limb darkening, and transit dilution likewise belong in the relevant physical model. A flexible trend fitted through a dip could absorb that dip; the masked fit in the real example illustrates why the fitting procedure is part of the measurement. Exposure-averaged light curves and phase folding of a light curve make two further limitations explicit.

Neighboring freely falling comoving particles in a homogeneous Friedmann-Lemaître-Robertson-Walker metric have separation $\mathbf{r}=a(t)\mathbf{\chi }$ and relative acceleration $\ddot{\mathbf{r}}=(\ddot{a}/a)\mathbf{r}$, where $t$ is cosmic time and $\mathbf{\chi }$ is constant. This is a local tidal force effect, of order $H^{2}r$ when the dimensionless background parameters are of order one. In a weak-field central-mass model retaining the homogeneous background tide, one writes $\ddot{r}=-GM/r^2+{\ell }^2/r^3+(\ddot{a}/a)r$. The detailed residual tide depends on the cosmological embedding; a bound orbit is not required to satisfy $\dot{r}=Hr$. Comparing $H^2$ with $GM/r^3$ estimates whether that background tide can affect its dynamics. A locally vacuum region with a cosmological constant has the corresponding outward term $\Lambda c^{2}r/3$.

The Newtonian potential of a point mass $M$ at the origin isIt is the free-space Green function of the Laplacian with the normalization of the Poisson equation for Newtonian gravity. Differentiating gives inward gravitational acceleration $-GM/r^2$. With $z=r/r_0$, the normalized Newtonian gravitational potential is $\Phi /(GM/r_0)=-1/z$; at $z=1,2,4$ it equals $-1,-1/2,-1/4$.

For a radially outward launch from $r_0$, conservation of mechanical energy in the fixed-source test particle approximation gives the minimum escape velocity:At this threshold the speed tends to zero only at infinite separation. The model neglects air resistance, external tidal forces, and general relativity. ^[ref]

A circular orbit at $r_0$ instead requires tangential speed $v_c=\sqrt{GM/r_0}$. Choosing units with $G=M=r_0=1$ gives $g_0=1$, $v_c=1$, and $v_{\mathrm{esc}}=\sqrt{2}\simeq 1.41421$. A radial launch at speed $1$ reaches $r=2$ and returns; the same speed directed tangentially gives a circular orbit. ^[ref]

An encounter between galaxies or clusters mainly couples their collective gravitational fields; it need not involve physical impacts between individual stars. Direct stellar collision rates depend on number density, relative speed and gravitational focusing. For two stars the collision cross-section is approximatelyIn an ordinary galaxy, the small geometrical stellar filling factor makes such impacts rare even while the galaxy undergoes a major encounter. Dense stellar cores are an important exception, not a reason to treat all galaxies as collisional fluids.

Tidal forces depend on the gradient of the perturber's field. For a distant fast encounter of mass $M_p$, impact parameter $b$ and speed $V$, the impulse scale for a star at displacement $r$ isA more specific result illustrates the geometry. Put the straight-line perturber trajectory in the $(x,y)$ plane along $y$, with closest approach at $x=b$. In the tidal and impulse limits,for an initially spherical, nonstreaming target with randomly sampled internal phases. The cross term averages to zero, and $⟨x^2+z^2⟩=2⟨r^2⟩/3$. Heating, expansion and tidal stripping preferentially affect weakly bound outer stars. A slow encounter can be adiabatically shielded rather than impulsive; a penetrating or strongly resonant encounter needs more than this distant straight-line approximation.

Disc rotation and orbital orientation matter. Prograde passages can drive stronger coherent distortions than corresponding retrograde passages. Tidal bridges, long tails, asymmetric arms and shells can record the encounter. Restricted gravitational calculations already reproduced prominent bridge and tail structures without requiring individual stellar collisions; full self-consistent models additionally follow the gravity of the responding discs and haloes.^[ref]

A first high-speed passage need not bind the pair. Capture or merging requires orbital energy and angular momentum to be redistributed into internal motions, escaping material, gas radiation or the surrounding halo. Dynamical friction transfers orbital energy to a gravitational wake and can bring a bound satellite inward. In a major merger, the potential changes on orbital times and violent relaxation redistributes stellar energies; phase mixing then smooths fine orbital structure. These are collective collisionless processes, distinct from slow two-body relaxation.

Gas behaves differently from stars: gas streams can shock, dissipate energy, lose angular momentum and fuel a central starburst or active nucleus. Gas-poor mergers are correspondingly more nearly dissipationless. Outcomes depend on mass ratio, orbital angular momentum, relative velocity, halo structure and gas fraction: they range from a modest flyby and tidal debris to disruption of a satellite or a merged remnant. Slow encounters in a group are generally more favorable to merging than otherwise comparable very fast unbound passages, but no single speed or mass-ratio threshold follows without an explicit structural model.

The strong equivalence principle extends universality of free fall to bodies with appreciable gravitational self-energy, and states that the outcome of any local experiment, including a gravitational experiment, in a freely falling laboratory is independent of its location and velocity. In a sufficiently small freely falling laboratory the external gravitational influence can be removed to first order by choosing a local inertial frame.

This does not remove tidal forces or the Riemann curvature tensor. The laboratory must be small compared with the length scale over which the external gravitational field varies. These qualifications distinguish a local equivalence statement from the incorrect assertion that gravity can be transformed away throughout an extended curved region.

Outside the core, $f^{\prime \prime }=0$ and the metric is locally flat: set $\tilde{\varphi }=(1-\lambda )\varphi$. Static test particles have no Newtonian-like attraction because $g_{00}$ is constant and ${\Gamma }^r{}_{00}=0$. Local trajectories are straight lines in locally unwrapped coordinates, so a very thin dark core is hard to identify through ordinary orbital accelerations or local tidal forces.

The global angular range is nevertheless $2\pi (1-\lambda )$ rather than $2\pi$. The missing angle $\delta =2\pi \lambda$ can give gravitational lensing with two nearly identical, essentially undistorted images of a background galaxy, and multiple such pairs aligned along the projected object. For a small deficit, a straight string and favourable geometry, their angular separation is of order $\delta D_{ls}/D_s$, with an additional inclination factor for a tilted string. The matching spectra and common line of image pairs distinguish this from an ordinary compact lens.

A moving cosmic string can also induce a step in the cosmic microwave background temperature across its path. That signature needs motion in addition to the static geometry supplied here. These are consequences of the exterior conical geometry; interpreting the inferred deficit as a positive mass per length requires the regular-axis/core assumptions distinguished in part b.

Approximate Earth by an isolated, nonrotating, spherically symmetric body. By Birkhoff's theorem, the exterior vacuum metric tensor is the Schwarzschild metric. With $x^0=ct$ it isHere $r$ is the areal radius, $\theta$ and $\varphi$ are spherical angles, $t$ is Schwarzschild time normalized to the proper time of a static observer at infinity, $M$ is Earth's mass, $G$ the gravitational constant, and $c$ the speed of light. The nonzero diagonal metric tensor components are $g_{00}=-f$, $g_{rr}=f^{-1}$, $g_{\theta \theta }=r^2$, $g_{\varphi \varphi }=r^2{\sin }^2\theta$.

The apple is a test particle; neglect its effect on the metric tensor, air resistance, Earth's rotation and nonspherical structure, and external tidal forces. Near Earth's surface, $GM/(R_E{c}^2)\simeq 7\times 10^{-10}$, $h/R_E\ll 1$, and the falling speed remains much smaller than $c$. Thus the Newtonian limit of this metric tensor is appropriate, while the exact radial geodesic is available as a check. A static observer under the tree has proper time $T=\sqrt{f(R_E)}t$, up to an additive constant. Such an observer is supported by the ground and has nonzero proper acceleration; the released apple instead has zero proper acceleration.

General relativity models spacetime by a smooth four-dimensional Lorentzian manifold whose metric tensor governs clock readings, causal cones and free-particle motion. The Einstein equivalence principle motivates local special relativity and universal geodesic free fall for test bodies; it does not by itself select a unique gravitational field equation. We assume a Levi-Civita connection, a metric tensor as the gravitational field with no additional dynamical gravitational variables, and matter with a covariantly conserved stress-energy tensor. Units have $c=1$.

The existence of local inertial frames is a local statement. At any regular event choose an orthonormal frame and Riemann normal coordinates, so that $g_{ab}={\eta }_{ab}$ and ${\Gamma }^a{}_{bc}=0$ at the event, with metric signature $(+---)$. Metric compatibility then gives ${\partial }_c{g}_{ab}=0$ there. The geodesic equation at that point is the special-relativistic inertial equation. Along a timelike geodesic, a parallel-transported orthonormal frame gives Fermi normal coordinates, eliminating the connection along the central worldline. These constructions do not in general eliminate second derivatives of the metric tensor in a neighborhood. The Riemann curvature tensor is built from the derivatives of the connection coefficients as well as their quadratic terms; nonzero Riemann curvature tensor components can remain when all connection coefficients vanish at a chosen point. Conversely, Polar coordinates in Minkowski spacetime can have nonzero connection coefficients but zero Riemann curvature tensor. Connection coefficients alone are therefore not an invariant measure of gravity.

For the tensorial content, consider a two-parameter family of geodesics parametrized by an affine parameter, with tangent $U$ and connecting field $J$. Their coordinate construction gives $[U,J]=0$, and torsion-freeness gives ${\nabla }_{U}J={\nabla }_{J}U$. Since ${\nabla }_{U}U=0$, commuting derivatives in the convention of Question 1 yields geodesic deviation:The left side is covariant relative acceleration, rather than the second coordinate derivative of an arbitrary separation. In an orthonormal frame in free fall with $U=(1,0,0,0)$, its spatial components are $D^2{J}^i/d{\tau }^2=-R^i{}_{00j}{J}^j$. This is a measurable tidal force, and a coordinate transformation cannot make a nonzero Riemann curvature tensor disappear. It is the distinction between gravity at one event and the behavior of an extended freely falling laboratory.

In a static Newtonian limit, $g_{00}=1+2\Phi$ and ${\Gamma }^i{}_{00}={\partial }_i\Phi$, so the test-body acceleration is $-\nabla \Phi$. The curvature-sign convention in Killing derivative identities gives $R^i{}_{00j}\simeq {\partial }_i{\partial }_j\Phi$. Hence relative acceleration is the negative Hessian of the gravitational potential. For the Newtonian potential of a point mass $\Phi =-GM/r$, the Hessian has radial eigenvalue $-2GM/r^3$ and two tangential eigenvalues $GM/r^3$: radial separations stretch and tangential ones compress. Outside the source the Ricci tensor can vanish while the Riemann curvature tensor remains nonzero. Its vacuum tidal part is encoded by the Weyl tensor; vacuum does not mean flatness. Gravitational waves likewise carry Riemann curvature tensor disturbances without a local material source.

To obtain the field equations, make the additional dynamical assumptions explicit: locality and diffeomorphism invariance of general relativity; a symmetric metric tensor field equation using at most second derivatives of the metric tensor; a geometrical left-hand side linear in the Riemann curvature tensor with constant coefficients; and compatibility with local matter conservation. These restrictions exclude higher-curvature or extra-field alternatives. Under them the available rank-two curvature contractions give the ansatzThe contracted Bianchi identity and metric compatibility giveFor this to vanish identically for arbitrary metric tensors, $b=-a/2$. Assuming $a\ne 0$ and absorbing its normalization into the coupling givesThe permitted constant term is the cosmological constant. Its value is not fixed by the equivalence principle. Conservation of $T$ follows for generally covariant matter on its equations of motion: varying the matter action under a compactly supported infinitesimal diffeomorphism, then integrating its term $T^{ab}{\nabla }_a{\xi }_b$ by parts, gives ${\nabla }_a{T}^{ab}=0$. The divergence-free Einstein tensor is consistent with this identity.

The coupling sign must be fixed in the actual curvature-sign convention in Killing derivative identities, with $R_{ab}=R^c{}_{acb}$. For the static weak-field approximation, $R_{00}\simeq -\Delta \Phi$. Set $\Lambda =0$ for the local asymptotically flat Newtonian comparison. In four dimensions the trace equation gives $R=-\kappa T$, so $R_{00}=\kappa (T_{00}-\frac{1}{2}{g}_{00}T)\simeq \kappa \rho /2$ for slowly moving matter. Matching $\Delta \Phi =4\pi G\rho$ forces $\kappa =-8\pi G$. Thus, with this convention,Under the opposite curvature-sign convention in Killing derivative identities the coupling sign reverses, as in Einstein-equation coupling under reversed curvature. This sign statement is compatible with the supplied linearized handout. Local inertial frames remove the connection at an event; tidal Riemann curvature tensor effects remain tensorial, and the field equations require dynamical assumptions beyond the equivalence principle.

Take $z$ positive outward and define the separation as upper minus lower position. The Newtonian potential of a point mass gives the acceleration $a(z)=-{\partial }_z\varphi =-GM/z^2$. Expanding its difference across a separation $\zeta \ll z$ gives the tidal force equationThe upper particle falls slightly less rapidly, so the separation increases. If the upper particle carries the electric charge and the electric field points outward, the extra relative acceleration is $qE/m$. Reversing the field or placing the charge on the lower particle changes its sign. Write this signed acceleration as $a_e$ and put $k=2GM/z^3$. With initial relative velocity zero and nearly constant $k$ and $a_e$, the ordinary differential equation has solutionFor $k{t}^2\ll 1$, the leading changes due to the tidal force and electric field are respectivelywith the stated sign convention. Thus, for a nonzero fixed electric acceleration,A laboratory of size $L$ has $|{\zeta }_0|\le L$, so this ratio tends to zero as $L\to 0$. This is local tidal acceleration scaling: gravitational relative acceleration is proportional to separation, whereas the electric difference persists even at coincident positions.

The duration also matters. At fixed time, $\Delta {\zeta }_g\to 0$ in absolute size, but $\Delta {\zeta }_g/{\zeta }_0\simeq (GM/z^3)t^2$ need not tend to zero. To keep the electrically displaced particles in the shrinking laboratory, take $|a_e|t^2\lesssim L$; then $|\Delta {\zeta }_g|/L\lesssim GML/(z^3|a_e|)\to 0$. Alternatively a stationary terrestrial laboratory has a free-fall time of order $\sqrt{L/g}$, giving $|\Delta {\zeta }_g|/|{\zeta }_0|=O(L/z)$. Both make precise the small spacetime-region limit. An extended, long-running experiment can still measure arbitrarily small tidal forces.

The weak equivalence principle says that freely falling test particles with the same initial data have composition-independent motion. In particular, a freely falling frame removes the common gravitational acceleration of both particles. The electric accelerations depend on charge and cannot be removed for charged and uncharged particles simultaneously by the same change of frame.

The Einstein equivalence principle extends the freely falling laboratory statement to all local nongravitational experiments: their outcomes obey special relativity, with local Lorentz invariance and local position invariance. The strong equivalence principle also includes local gravitational experiments and bodies with significant gravitational binding energy. It asserts universality of their free fall and independence of local experimental outcomes from the laboratory's location and velocity. The two-particle experiment illustrates the local distinction between gravity and a nongravitational force; by itself it does not establish this stronger claim about self-gravitating bodies.

In a metric theory, the common freely falling trajectories are timelike geodesics of a Lorentzian metric. Normal coordinates make $g_{ab}={\eta }_{ab}$ and the Christoffel symbols zero at one point. Fermi normal coordinates do the same along a reference timelike geodesic. The residual effects are controlled by the Riemann curvature tensor through geodesic deviation, and enter the local metric at quadratic order in distance. Curvature is a physical tidal field, so no coordinate choice removes it over a finite region.

A relativistic theory should consequently couple matter universally to the metric tensor, recover nongravitational laws in local inertial frames and recover the Newtonian limit. General relativity implements this with the Einstein field equations; the contracted Bianchi identity is consistent with stress-energy conservation. The equivalence principle motivates this geometric description but does not alone uniquely determine the field equations. Extra gravitational fields may satisfy the Einstein equivalence principle for ordinary matter while failing the strong equivalence principle because local gravitational measurements or self-gravitating trajectories depend on those fields.

Synchronous radius and its limits. Kepler's third law for a circumplanetary orbit givesHere $G$ is the gravitational constant. A planet-synchronous orbit repeats after one spin period. A planet-stationary orbit additionally requires a circular orbit, zero orbital inclination to the equator and prograde motion. Then an antenna fixed on the ground points continuously at the same satellite; merely matching the orbital period does not give this property.

Write the planetary radius as $R_p=(3M/(4\pi \rho ){)}^{1/3}$. Requiring the semi-major axis to exceed $R_p$ gives $P^2>3\pi /(G\rho )$. The stellar tidal force restricts the circumplanetary orbit to the Hill sphere, of radius $r_H=a(M/(3M_{\star }){)}^{1/3}$. ThusThese are the surface and Hill-radius constraints in the idealized spherical, small-$M/M_{\star }$ model. Long-term prograde stability generally requires a radius appreciably inside the Hill sphere; for nonzero orbital eccentricity, the surface constraint applies to the periapsis, not just to the semi-major axis.

Collision time and the launch population. A phase-mixed isotropic swarm occupies a shell of radial scale $e{a}_s$. Its number density scales as $N_s/(e{a}_s^3)$, its relative speed as $a_s/P$, and its geometric collision cross-section as $r^2$. Consequently the total collision rate scales as $N_s^2{r}^2/(Pe{a}_s^2)$. To give the numerical normalization used here, adopt an effective shell volume $V=4\pi e{a}_s^3$. At a fixed position the velocity directions are uniformly distributed in the tangent plane, so their mean relative speed is $4v_K/\pi =8a_s/P$. For an unordered pair, the collision cross-section is $4\pi r^2$, givingThus the mean interval between collisions is $A/N_s^2$ in this large-population kinetic model. The effective shell width fixes an order-one coefficient: small orbital eccentricity and random planes alone do not specify a unique radial probability density. Phase mixing, negligible gravitational focusing, $r\ll e{a}_s$, and uncorrelated encounters are implicit in this estimate. Exactly identical orbital periods with perfectly fixed phases do not themselves produce a memoryless collision process.

For nearly planet-stationary orbits with $e\ll I\ll 1$, the swarm volume is smaller by a factor of order $I$, while the relative speed is smaller by the same factor, since vertical motion of scale $I{v}_K$ dominates the eccentric motion. The two changes cancel in the rate $n\sigma v$: there is no parametric factor $I$ in the collision time within the same phase-mixed kinetic approximation. Numerical factors and phase correlations can differ. This is not an argument that bringing all satellites into one nearly circular plane makes their phases random.

With $N_s=Rt$, an Inhomogeneous Poisson process has cumulative hazard functionThe survival function of the first collision is $\exp [-H(t)]$. Setting the expected number of collisions to one givesThis is a characteristic first-event population, with probability $1-e^{-1}$ of an earlier event. It is not a median: the median has an additional factor $(\log 2{)}^{1/3}$.

Which population collides next? Immediately after the first disruption let $S=N_{s0}-2$ and $F=2x^{-3}$. The latter follows from mass conservation for equal-density spherical pieces. Using the same unordered-pair counting and geometric collision cross-sections as above, define $b=(1+x{)}^2/2$; thenHere $S^2,F^2$ stand for the usual large-population approximations to $S(S-1),F(F-1)$. The factor two distinguishing identical and different species is essential. The probability that the next collision is fragment–fragment isIt is the most likely type if $S<4/[x(1+x{)}^2]$ and $S<2/x^2$. It has probability greater than one half ifFor $x\ll 1$, a strongly fragment-dominated next event therefore requires $Sx\ll 4$; the fragment–satellite comparison is more restrictive than the satellite–satellite comparison. If $S=0$, the next event must be fragment–fragment, provided fragments remain.

Population equations and the normalization discrepancy. Every satellite–satellite event produces $2x^{-3}$ fragments and destroys two satellites. Every satellite–fragment event produces a net $x^{-3}-1$ fragments and destroys one satellite; every fragment–fragment event destroys two fragments. Consistent collision counting therefore givesFor $x\ll 1$, $b\simeq 1/2$ and $x^{-3}-1\simeq x^{-3}$. The last two terms in $A{\dot{N}}_f$ are then twice those in the printed equation. This is a genuine factor-of-two inconsistency: equal-size fragment pairs have a collision cross-section smaller by $x^2$, so their event rate must be $x^2{N}_f^2/A$ if the satellite event rate is $N_s^2/A$; destroying both fragments necessarily gives the sink $-2x^2{N}_f^2/A$.

If the printed equation is taken as a prescribed approximate rate model instead, its implicit event rates are $N_s^2/A$, $N_f{N}_s/(4A)$ and $x^2{N}_f^2/(2A)$. Under precisely that mixed normalization, its corresponding satellite equation isThe printed source term also neglects the consumed fragment in a satellite–fragment event, a legitimate relative $O(x^3)$ approximation. That approximation does not repair the pair-counting discrepancy.

The ensuing cascade. This is a two-size fragmentation cascade. First use the printed approximate model, dropping satellite–satellite events as requested. Put $S=N_s$, $F=N_f$, $k=4x^2$ and take $x<1/2$. ThenIntegrating this linear differential equation, with $F_0=2x^{-3}$ and $S_0=N_{s0}-2>0$, givesInitially the collisional cascade grows if $S_0>8x^2$, with approximate early exponential growth time $4A{x}^3/S_0$ when the satellite population is nearly fixed. The fragment population reaches its maximum atAfterwards satellites are depleted and fragment–fragment losses dominate. For $S_0>0$ the continuum solution has $S\to 0$, $F\to 0$, withWhen no satellites remain initially, $F(t)=F_0/[1+x^2{F}_{0}t/A]$ directly. Small integer populations eventually invalidate these deterministic differential equations.

The consistently counted model has, to leading order in $x$, exactly the same curve $F(S)$ and peak, but both retained time derivatives are twice as large. Its growth time is $2A{x}^3/S_0$, and $F(t)\sim A/(2x^{2}t)$. Keeping $(1+x{)}^2$ and the consumed fragment replaces $k$ by $4x^2/(1+x{)}^2$ and $x^{-3}$ in the curve by $x^{-3}-1$. This explicitly separates the physical collision bookkeeping from the printed normalization while giving the evolution under both conventions.

At $x_n=0$ and $y_n=hn$, every radial pair displacement vanishes, so $f_x^n=0$. The particles with indices $n+l$ and $n-l$ give opposite azimuthal accelerations,which cancel pairwise. The acceleration sum has absolute convergence because ${\sum }_{l\ge 1}{l}^{-2}<\infty$. With zero velocities and $x_n=0$, both the Coriolis force and radial tidal force vanish. The equally spaced row is therefore an equilibrium of the infinite gravitating particle lattice. Its absolute gravitational potential would require an additional reference or regularization; the acceleration used here is well defined.

Radiative drag is only one of several loss mechanisms. Around a star, radiation-pressure blowout can eject small fragments; its $\beta =1/2$ threshold applies specifically to zero-kick release from a circular parent orbit. Stellar-wind drag and gas drag can drive planetary migration, while sublimation destroys grains approaching high-temperature regions. Collisional cascades destroy or fragment grains and can feed the unbound size range. Planetary scattering can cause ejection, collision with a planet, or a stellar impact; resonant trapping of dust can instead delay planetary migration.

For circumplanetary orbits, collisions with the planet or its satellites, disruption in collisions, and escape under stellar tidal forces are additional losses. Orbits near or outside the Hill sphere need not remain planet-bound. Radiation pressure on circumplanetary dust can excite planetocentric orbital eccentricity or unbind very small grains; it need not act only through slow Poynting–Robertson drag. For charged grains, the Lorentz force in stellar or planetary magnetic fields can alter or destabilize an orbit. Shadowing of circumplanetary dust changes the radiation-force average and can reduce the quoted decay rate. Which mechanism dominates depends on grain size and composition, environment, orbit orientation and the available collision or gas density.

Work in the planet's rest frame, neglecting its recoil and tidal forces during the close encounter. A hyperbolic Kepler orbit has positive specific orbital energy $\epsilon =v_{\infty }^2/2$ and specific angular momentum $h=b{v}_{\infty }$, where $b$ is the impact parameter. Its orbital eccentricity isAt the pericentre distance, $r_{\min }=h^2/[\mu (1+e_h)]$. Substituting $h^2={\mu }^2(e_h^2-1)/v_{\infty }^2$ yieldsThe two asymptotes have true anomalies $\pm f_{\infty }$, where $\cos f_{\infty }=-1/e_h$. Their directions of motion give the gravitational scattering angle $\theta =2f_{\infty }-\pi$, so $\sin (\theta /2)=-\cos f_{\infty }=1/e_h$. ThereforeAs a useful independent form, $b^2=r_{\min }^2+2\mu r_{\min }/v_{\infty }^2$, displaying gravitational focusing. Here $\mu =G{M}_p$ to the stated small-comet accuracy; the exact relative-coordinate formula replaces it by $G(M_p+M)$.