Astronomy & Celestial Mechanics
the sky, computed
Concept

Lagrange coefficients — where it appears

Four functions, f, g and their time derivatives, that carry a body's position and velocity forward along a two-body orbit. Each new vector is a weighted sum of the initial ones, exact and free of orbital elements, so the method works for any conic.

Named by 3 essays across one field — each of them below, with the objects they name alongside it.

One admissible root, 0.01% from the truth. Gauss's reduction of three directions to a distance, drawn as the two relations whose intersection it is. Three observations of Ceres on days 0, 20, 40 of an arc, generated from its elements and used only as sight directions — no range, no radial velocity. The rising curve is geometry: the heliocentric distance a candidate at geocentric distance ρ₂ would have, r₂² = ρ₂² + 2ρ₂(R₂·L̂₂) + R₂², which contains no dynamics at all. The falling curve is dynamics: ρ₂ = A + µB/r₂³, with A and B built from the three sight vectors, the three observer positions and the three times, and containing no orbit. Eliminating ρ₂ between them gives r₂⁸ + a r₂⁶ + b r₂³ + c = 0 — an eighth-degree equation, from a problem with exactly as many equations as unknowns. Here they cross once at a positive ρ₂, at r₂ = 2.5893 AU against the true 2.5890. The other 2 real roots are rejected not by fitting but by sign: the ρ₂ each implies is negative, and an object behind the observer was not the thing observed.

Three observations and no orbit at all

Three directions in space give six numbers for the six elements of an orbit, which sounds like a solved problem. The algebra that solves it is of the eighth degree, and for a near-Earth asteroid three perfect observations can be consistent with three different orbits.

orbits · Orbit determination Four numbers that carry a position and a velocity round an orbit. The Lagrange coefficients f, g, ḟ and ġ — g and ḟ scaled by the mean motion so that all four are pure numbers — against time as a fraction of the period, for an ellipse of eccentricity 0.6 and periapsis 0.5 AU, each computed from one universal-variable solve starting at periapsis. The new position is f times the old position plus g times the old velocity, and the new velocity ḟ times the old position plus ġ times the old velocity: the state at any time is a linear combination of the state at the start, with no semi-major axis, eccentricity or anomaly formed. f and ġ return to one and g and ḟ to zero after a full period, which is the orbit closing. At every one of the 401 instants the combination f ġ − ḟ g equals one to 6·10⁻¹⁵, which is the conservation of angular momentum written as the area of a parallelogram.

Four numbers that carry a state without an orbit

Once the universal variable is solved, a body's position and velocity at any later time are a linear combination of its position and velocity now, with four coefficients and no orbital elements formed at all. The four obey one identity exactly — the conservation of angular momentum, written as an area — and its failure measures the error a calculation has made without knowing the answer.

orbits · Universal variables Errors the identity sees and errors it cannot. Propagations of an ellipse of eccentricity 0.6 through 0.37 of a period, spoiled four ways at sizes from 10⁻¹³ to a hundredth, plotted by their true relative position error across and by how far their coefficients miss f ġ − ḟ g = 1 up (both logarithmic; the bottom row is a miss below the arithmetic's resolution). A solve stopped early and coefficients jostled at random break the identity in proportion to the damage. A wrong time step and a wrong gravitational parameter do not break it at all, at any size — errors up to 0.009 of the position sit on the floor — because each produces an exact Kepler map, for a different question. One equation on four numbers leaves three directions unchecked, and two of them are the two inputs every propagation depends on.

A check that cannot see the clock

The four Lagrange coefficients that carry an orbit's state forward satisfy one identity exactly, and a propagation that breaks it is wrong by about as much as it breaks it. But one equation on four numbers leaves three directions unchecked, and two of them are the inputs every propagation depends on: the time and the strength of gravity. A state carried forward for the wrong time, or under the wrong gravitational parameter, satisfies the identity perfectly — and a solver can be written so that its own errors take exactly that form.

orbits · Universal variables

Named alongside it

The objects these essays reach for when they reach for this one.

Angular momentumEphemerisKepler equationOrbital elementsStumpff functionsUniversal variableAstrometryConditioningGauss's methodGeocentric distanceGravitational parameterInitial orbit determination

All concepts