Orbits

The precession that switches the cycle off

A distant companion can trade an orbit's inclination for its eccentricity, over and over, and drive a pericentre almost onto the central body. General relativity's own precession competes with the mechanism for the same pericentre, and when it wins the cycle stops — sharply, at a ratio of one, which puts a floor on how close anything can be delivered.

Assumes Kozai–Lidov, Precession and Relativistic orbits.

An inclination can turn into an eccentricity, and the previous rung of this anchor established the mechanism: a distant companion, averaged over both orbits, conserves one combination of the inner orbit’s eccentricity and inclination and lets the two exchange. Above a critical inclination the exchange becomes an oscillation that can drive a nearly circular orbit to an eccentricity of 0.99 and back.

That is the mechanism in isolation. Nothing in astronomy is in isolation, and the thing that most reliably interferes with it is the inner orbit’s own relativity.

Relativity switches the cycle off, halving its reach at a ratio of 0.80. The greatest eccentricity a Kozai–Lidov cycle reaches, against the strength of the orbit's own relativistic pericentre precession, measured in units of the cycle's own precession rate at zero eccentricity. The horizontal axis is logarithmic and spans three decades. At the left the relativistic term is negligible and the cycle reaches 0.838, which is the closed-form value for a start at 65 degrees and is what the integration is checked against. At the right it is gone. The mechanism depends on the pericentre staying put while the outer body pulls on the same side of the orbit for a whole cycle, and the relativistic precession is a competing rotation of that same pericentre; when it is faster, the pull averages away. The threshold sits near one by construction and the transition is sharp rather than gradual, with the reach halved at 0.80. What makes it matter is where the relativistic term is largest: it grows as the pericentre falls, so it strengthens exactly as the cycle drives the orbit inward, and it therefore sets a floor on the pericentre distance that this mechanism can deliver a body to.
Fig. 1 The greatest eccentricity a cycle reaches, against the strength of the orbit’s own relativistic pericentre precession, measured in units of the cycle’s own precession rate at zero eccentricity. At the left the relativistic term is negligible and the cycle reaches its full closed-form value; at the right it is gone. The threshold sits near one and the transition is sharp rather than gradual, with the reach halved at a ratio of about 0.8.

Why the mechanism needs the pericentre to stay put

The cycle is not a resonance between periods. It is a secular effect: both orbits are averaged over, so the individual positions of the bodies do not matter, and what is left is a slow exchange between the shapes of two rings of mass.

What makes the exchange coherent is the argument of pericentre. The outer body’s averaged quadrupole field is not spherically symmetric about the inner orbit’s plane, so the torque it exerts on the inner ring depends on where that ring’s long axis points. If the long axis stays put, the torque acts in the same sense for a long time and the eccentricity climbs. If the long axis rotates, the torque reverses as it goes round and averages away.

So the mechanism is a competition between two rotations of the same pericentre: the one the outer body drives, and every other one.

The phase plane at Θ = 0.4226, and the curve that divides it. Level curves of the averaged Hamiltonian in the pericentre argument and the eccentricity, at fixed √(1−e²)cos i = 0.4226 — the whole of the dynamics, because the averaged problem has one degree of freedom. Two kinds of curve, separated by one: outside the shaded boundary ω circulates, running through every value while the eccentricity wobbles a little; inside it ω librates about 90° or 270° and the eccentricity swings between 0.67-scale extremes. The dividing curve is the level through e = 0, and it is exactly that because F contains ω only multiplied by e² — a circular orbit is a level curve of its own, which is why an orbit that starts circular starts on the separatrix and always reaches the same maximum. The fixed point sits at ω = 90°, e = 0.6741, where j⁴ = (5/3)Θ²; an orbit placed exactly there never changes at all.
Fig. 2 The phase portrait the competition takes place in. The horizontal axis is the argument of pericentre and the vertical is the eccentricity, and the closed curves are trajectories of the averaged problem. Above the critical inclination there are fixed points at ninety degrees, and trajectories circulate around them: a body starting nearly circular runs up the left side of an island, reaches high eccentricity, and comes back down. Anything that adds a steady drift in the horizontal direction shears those islands, and enough drift destroys them.

The competitor that cannot be avoided

An orbit precesses for many reasons. A planet’s oblateness precesses the orbits of its moons; the other planets precess each other’s, through the secular exchange that gives no planet an eccentricity of its own; and a tidal bulge on either body precesses the pair. Those depend on the system, and in principle a system can be found without them.

Relativity is not optional. Any two bodies in orbit have a pericentre that advances, at a rate

ω˙GR=3G(M1+M2)nc2a(1e2),\dot\omega_{\rm GR} = \frac{3 G (M_1+M_2)\, n}{c^2 a (1 - e^2)},

and the rate has two properties that matter here. It is always in the same direction, prograde. And it grows as the pericentre falls — the factor of one over one minus the square of the eccentricity is exactly the term that blows up as the cycle drives the orbit inward.

The relativistic advance, computed for five orbits. The general-relativistic apsidal advance, in arcseconds per century, computed from each body's own semi-major axis, eccentricity and period. It falls steeply with distance and rises with eccentricity, which is why Mercury at 42.98 is the only planet whose value was within reach of nineteenth-century astrometry — and why Icarus, a near-Earth asteroid on a far more eccentric orbit, was measured for the same effect in the 1960s.
Fig. 3 The precession budget of a real orbit, with the contributions separated. For Mercury the relativistic term is the smallest of them and it is the forty-three arcseconds a century that broke Newtonian gravity. For a compact inner binary in a triple the ordering is reversed: the relativistic term dominates everything, because it scales as the inverse five halves power of the semi-major axis and nothing else in the budget rises that fast.

The threshold, and why it is sharp

Putting the relativistic precession into the averaged equations adds one term to the rate of change of the argument of pericentre, and the hero figure is the result of integrating the modified system for a range of values of that term.

Two features of the result are worth separating. The first is that the transition is at a ratio of order one — which is not a surprise, since the only dimensionless quantity in the problem is the ratio of the two precession rates. The second is that it is sharp: the greatest eccentricity falls from its full value to a third of it over less than a decade in the ratio.

The sharpness comes from the shape of the phase portrait. Adding a drift shifts the fixed points and shrinks the islands; when the drift exceeds a threshold the islands disappear entirely, and there is no intermediate configuration in which a small island supports a small oscillation. What is on the other side is not a weaker cycle but no cycle.

An eccentricity and an inclination trading, at 70° of mutual tilt. The secular equations integrated from a nearly circular orbit (e = 0.02) inclined at 70° to a distant perturber's plane, over three oscillations. Above, the eccentricity; below, the inclination, with the constant √(1−e²)cos i drawn as the flat line it is. The eccentricity climbs to 0.8972 and the inclination falls to 39.25° at the same instant, and neither is a coincidence: the product is fixed, so one can only rise as the other falls. That floor is the same for every starting tilt — at maximum eccentricity j = √(5/3)Θ, so cos i = √(3/5) and the inclination arrives at 39.23° whether the orbit began at 50° or at 89°. The closed form for a circular start is e_max = √(1 − (5/3)cos²i₀) = 0.8972, which contains nothing about the perturber — not its mass, not its distance. Those set the clock and not the amplitude, and the period here is 4.54 Kozai times. What the figure cannot show is what happens at the top of the cycle in a real system: at e = 0.897 the pericentre is 0.1028 of the semi-major axis, where tides, general relativity or a stellar surface all intervene, and the quadrupole picture ends.
Fig. 4 An unquenched cycle in the time domain, for a start at seventy degrees. The eccentricity runs from nearly zero to near one and back, and the inclination runs the other way, with the conserved combination holding their product fixed. The period is set by the outer body’s mass and distance and is typically thousands to millions of the inner orbit’s own periods — long enough that the relativistic precession, which is small per orbit, accumulates to something comparable over one cycle.

The floor it puts on delivery

The reason anyone cares is that the mechanism is the leading candidate for making close-in giant planets and close binaries, by a route that does not require anything to migrate through a disc.

The route is: a distant companion drives the inner orbit to high eccentricity; at the pericentre of that eccentric orbit, tides raised on the inner body dissipate energy; the orbit shrinks and circularises at roughly twice its pericentre distance. It is called high-eccentricity migration, and unlike disc migration it can happen long after the gas has gone and it naturally produces the misaligned orbits that are observed.

Everything close in is circular, and nothing else has to be. Orbital eccentricity against period for fifteen real planets, with the tidal circularisation boundary computed from τ_e = (2/21)(Q′/n)(M_p/M⋆)(a/R_p)⁵ for a Jupiter with Q′ = 1e+6 and an age of 5 billion years. It falls at 6.2 days, and the reason it is a wall rather than a slope is the fifth power: at half the period the timescale is 91 times shorter. Nothing inside it has a measurable eccentricity, and outside it eccentricities run to 0.95 — which is the number to hold on to, because a planet on a 0.93 orbit at 111 days comes within 0.030 AU of its star at periastron, closer than Mercury, and is being circularised as it is observed.
Fig. 5 The tidal circularisation boundary the route delivers planets to: inside the crossing, dissipation removes the eccentricity within the system’s age; outside it, an eccentric orbit survives. High-eccentricity migration works by driving the pericentre inside this boundary while the apocentre is still far outside it, so the orbit is circularised from a highly eccentric state — which leaves the planet at about twice the pericentre it reached, and nowhere else.

The quench sets a floor on that pericentre. As the cycle drives the eccentricity up, the relativistic precession rate rises as one over one minus the square of the eccentricity, so a mechanism that was subdominant at the start of the cycle can become dominant partway through and stop the cycle before it reaches the eccentricity that would have delivered the planet.

The result is a maximum eccentricity that depends on the inner orbit’s semi-major axis, and hence a minimum pericentre. Planets cannot be delivered arbitrarily close by this route: below a certain final orbital period, relativity would have stopped the cycle that was supposed to produce them.

Relativity switches the cycle off, halving its reach at a ratio of 0.57. The greatest eccentricity a Kozai–Lidov cycle reaches, against the strength of the orbit's own relativistic pericentre precession, measured in units of the cycle's own precession rate at zero eccentricity. The horizontal axis is logarithmic and spans three decades. At the left the relativistic term is negligible and the cycle reaches 0.975, which is the closed-form value for a start at 80 degrees and is what the integration is checked against. At the right it is gone. The mechanism depends on the pericentre staying put while the outer body pulls on the same side of the orbit for a whole cycle, and the relativistic precession is a competing rotation of that same pericentre; when it is faster, the pull averages away. The threshold sits near one by construction and the transition is sharp rather than gradual, with the reach halved at 0.57. What makes it matter is where the relativistic term is largest: it grows as the pericentre falls, so it strengthens exactly as the cycle drives the orbit inward, and it therefore sets a floor on the pericentre distance that this mechanism can deliver a body to.
Fig. 6 The same quench curve for a steeper initial inclination, where the unquenched cycle reaches much higher eccentricity. The threshold ratio hardly moves — it is set by the competition between two rates and not by how far the cycle would have gone — but the unquenched value it falls from is larger, so the same relativistic term costs more in absolute eccentricity. Systems whose geometry would have delivered the deepest pericentres are the ones the quench most affects.

What sets the ratio in practice

Written out, the ratio of the two rates is

ϵGR    3GM12a23(1e22)3/2a14c2M2,\epsilon_{\rm GR} \;\sim\; \frac{3 G M_1^2\, a_2^3 (1-e_2^2)^{3/2}}{a_1^4\, c^2\, M_2},

with the subscripts distinguishing the inner pair from the outer companion. The dependences are steep in every variable. A closer inner orbit quenches harder, as the inverse fourth power. A more distant or lighter companion quenches harder, because it drives the cycle more weakly.

The practical consequence is a region of parameter space rather than a line. Triples with a wide, light outer companion and a compact inner binary are quenched; triples with a close, heavy companion and a wide inner binary are not. Both exist.

What was actually measured

None of this has been watched. A Kozai cycle in a stellar triple takes ten thousand to a million years, and the quench is a statement about what does not happen.

Three lines of evidence bear on it, and all are statistical.

The first is the period distribution of close binaries. A substantial fraction of short-period binaries have a distant third companion — a much higher fraction than for wide binaries — which is what the mechanism predicts, since the tight pairs are supposed to have been made by it. Establishing that a pair has a third member at all often requires resolving two components a Fourier transform can barely tell apart.

The second is the obliquity distribution of hot Jupiters. Planets delivered by high-eccentricity migration should arrive on orbits misaligned with their star’s spin, because the mechanism works by tilting the orbit; planets delivered through a disc should arrive aligned. Both populations are observed, in a proportion that depends on stellar type in a way tidal realignment explains. The third is the pile-up in orbital period. Hot Jupiters cluster near three days, and twice the tidal circularisation radius for a Jupiter around a Sun-like star is about that — a pile-up visible in the occurrence-rate census rather than in any individual system. Whether the inner edge of that pile-up shows the floor this essay predicts is a question about a few dozen objects, and the answer is not yet decisive.

Three rates, and which one is largest where

It helps to put numbers on the competition for two systems at opposite ends of the range.

Take a Jupiter at five astronomical units from a solar-mass star, with a stellar companion of half a solar mass at five hundred astronomical units. The cycle’s period is of order a million years. The relativistic precession of the inner orbit is about a hundredth of an arcsecond a century at zero eccentricity — utterly negligible — but at an eccentricity of 0.99 the factor of one over one minus the square of the eccentricity is fifty, and the pericentre has fallen by a factor of a hundred, so the rate has risen by five thousand. The cycle is unquenched at the start and quenched near the top, which is exactly the case the floor argument is about.

Now take a pair of white dwarfs at a hundredth of an astronomical unit with the same outer companion. The relativistic precession is a degree a year at zero eccentricity, and the cycle’s own rate, which falls as the cube of the ratio of the two semi-major axes, is nothing at all. This system never cycles.

The lesson is that the same triple geometry gives completely different answers depending on the inner orbit’s size, and that the dependence is a fourth power rather than a linear one. There is no gradual regime.

Relativity switches the cycle off, halving its reach at a ratio of 1.05. The greatest eccentricity a Kozai–Lidov cycle reaches, against the strength of the orbit's own relativistic pericentre precession, measured in units of the cycle's own precession rate at zero eccentricity. The horizontal axis is logarithmic and spans three decades. At the left the relativistic term is negligible and the cycle reaches 0.764, which is the closed-form value for a start at 60 degrees and is what the integration is checked against. At the right it is gone. The mechanism depends on the pericentre staying put while the outer body pulls on the same side of the orbit for a whole cycle, and the relativistic precession is a competing rotation of that same pericentre; when it is faster, the pull averages away. The threshold sits near one by construction and the transition is sharp rather than gradual, with the reach halved at 1.05. What makes it matter is where the relativistic term is largest: it grows as the pericentre falls, so it strengthens exactly as the cycle drives the orbit inward, and it therefore sets a floor on the pericentre distance that this mechanism can deliver a body to.
Fig. 7 The quench curve at an inclination only just above the critical value, where the unquenched cycle is weak to begin with. The threshold is in the same place — it is a ratio of rates — but the curve falls from a much lower starting value, so a system near the inclination threshold is quenched by a relativistic term that a strongly inclined system would shrug off. The two thresholds, in inclination and in precession ratio, therefore interact: a system can fail to cycle for either reason and the boundary between them is a curve rather than a pair of lines.

The other things that precess

Relativity is the quencher this essay is named for and it is not the only one. Anything that turns the argument of pericentre faster than the cycle does will do the same job, and in most real systems something else gets there first.

Tidal bulges. Each body raises a bulge on the other, and a bulge is a mass distribution that is not a point, so the orbit precesses. The rate falls as the eighth power of the separation — far steeper than the relativistic term’s fourth — so at very small pericentre distances the tidal term dominates everything. For a Jupiter approaching a solar-type star inside about ten stellar radii, tidal precession quenches the cycle before relativity does.

Rotational flattening. A rapidly rotating star is oblate, and an oblate primary precesses an orbit around it by the same mechanism that regresses a satellite’s node about the Earth. Young stars rotate fast, so the term is largest exactly when the systems in question are being assembled, and it also falls steeply with separation.

A third orbit. In a system with more than one planet, each planet precesses the others, and those secular rates are often faster than a distant companion’s cycle. That is the standard reason a compact multi-planet system is immune to a stellar companion that would otherwise wreck it: the planets protect each other by precessing each other, and removing one of them can destabilise the rest.

The competition is therefore between four rates with four different dependences on separation — the cycle itself scaling with the cube of the ratio of the two semi-major axes, relativity as the inverse fourth power of pericentre, tides as the inverse eighth, and mutual planetary precession as whatever the system’s own architecture gives. Which one wins is a different answer at every separation, and the boundaries between them are sharp for the reason the whole essay is about: a ratio of rates crossing one is not a gradual transition.

How high the eccentricity would have gone

The quench is a statement about a cycle that would otherwise have run, and quantifying it needs the unquenched amplitude — which at quadrupole order has a closed form worth having.

For an inner orbit starting nearly circular at inclination i0i_0 to the outer one, the cycle exchanges inclination for eccentricity while conserving 1e2cosi\sqrt{1-e^2}\cos i, and the maximum reached is

emax=153cos2i0.e_{\max} = \sqrt{1 - \tfrac{5}{3}\cos^2 i_0}.

Two things follow immediately. The cycle does nothing at all unless cos2i0<3/5\cos^2 i_0 < 3/5, which is the critical inclination band from 39.2° to 140.8°; and at exactly 90° the formula gives emax=1e_{\max} = 1, meaning the inner orbit is driven to a collision.

The consequence is the whole reason the mechanism matters. A pair starting at a separation of five astronomical units and 85° of mutual inclination reaches an eccentricity of 0.998, and its pericentre falls to a hundredth of an astronomical unit — from a wide orbit to a grazing one, with no dissipation and no change in the semi-major axis, purely by moving angular momentum between two orbits.

That is a factor of five hundred in pericentre distance from a geometric condition on an angle, and it is why the quench is decisive rather than a correction. The relativistic precession at the top of that excursion is thousands of times its value at the start, because the rate goes as the inverse fourth power of pericentre — so the mechanism generates, at its own maximum, exactly the term that switches it off.

It also explains why the mechanism is so often invoked and so rarely demonstrated in a specific system. The excursion happens once, over a million years, and leaves behind an orbit that carries no label saying how it got there — so the evidence is always a population’s shape rather than an individual’s history.

The same reasoning applies to the systems where nothing happened. A wide, coplanar pair with a distant companion is not evidence that the mechanism is weak; it is evidence that the inclination was outside the band, which is a statement about how the system formed rather than about how it evolved.

The wider lesson

The quench is a specific instance of a general rule about secular mechanisms, and the rule is worth stating on its own.

A secular effect works by accumulating a small perturbation coherently over many orbits. Anything that rotates the geometry of the accumulation faster than the effect itself destroys the coherence, and the competition is decided by a ratio of rates rather than by a ratio of forces. So a mechanism can be switched off by something far weaker than itself, provided the weaker thing acts on the right angle.

That is why precession budgets are worth constructing at all. A term that is negligible in the force is not negligible in the argument of pericentre, and the second is what secular theory is about.

There is a second reading of the same rule that is worth keeping, because it runs the other way. If a weak competitor can switch off a strong mechanism, then measuring whether a mechanism operated is a measurement of the competitor. A triple observed to have a tight, misaligned inner pair is a triple in which the relativistic term was small enough to lose, which bounds the inner orbit’s size at the time the cycle acted — a bound on a configuration that no longer exists, obtained from one that does.

That kind of inference is characteristic of this whole anchor. Secular mechanisms leave no trace of their operation except the final state, so every argument about them is a reconstruction: this arrangement is reachable and that one is not, and the observed population is the intersection of what is reachable with what survives. The quench is valuable precisely because it is a sharp boundary rather than a gradual falling-off, and a sharp boundary in a reconstruction is worth several soft ones.

One practical consequence deserves stating because it is the reason the quench appears in the literature far more often than its size would suggest. Numerical work on triples is usually done with the secular equations rather than by integrating every orbit, because the cycle takes millions of inner orbital periods and a direct integration is unaffordable. The relativistic term has to be added to those secular equations by hand — it is not there unless somebody puts it there — and for two decades a good deal of published work omitted it. The results were not wrong so much as answers to a different question: they described what a triple would do in a universe governed by Newtonian gravity, and the systems most affected were exactly the compact ones the calculations were aimed at. The correction is one line of code and it changes the predicted merger rate of compact binaries formed this way by an order of magnitude, which is a large return on a term nobody can measure directly.

One more starting argument of periapsis shows that the quenching does not depend on the phase.

Relativity switches the cycle off, halving its reach at a ratio of 0.69. The greatest eccentricity a Kozai–Lidov cycle reaches, against the strength of the orbit's own relativistic pericentre precession, measured in units of the cycle's own precession rate at zero eccentricity. The horizontal axis is logarithmic and spans three decades. At the left the relativistic term is negligible and the cycle reaches 0.897, which is the closed-form value for a start at 70 degrees and is what the integration is checked against. At the right it is gone. The mechanism depends on the pericentre staying put while the outer body pulls on the same side of the orbit for a whole cycle, and the relativistic precession is a competing rotation of that same pericentre; when it is faster, the pull averages away. The threshold sits near one by construction and the transition is sharp rather than gradual, with the reach halved at 0.69. What makes it matter is where the relativistic term is largest: it grows as the pericentre falls, so it strengthens exactly as the cycle drives the orbit inward, and it therefore sets a floor on the pericentre distance that this mechanism can deliver a body to.
Fig. 8 The same quenching started from a different argument of periapsis. The relativistic precession still switches the cycle off at the same separation, because the competition is between two rates and neither of them depends on where in the cycle the system happens to be.

Where the ladder goes

The next rung is the octupole term, which is what happens when the outer orbit is itself eccentric. At quadrupole order the inner orbit’s angular momentum along the outer orbit’s axis is conserved and the inclination cannot pass through ninety degrees; at octupole order it is not, and the inner orbit can flip from prograde to retrograde. That extends the reach of the mechanism enormously and is much harder to quench.

The other direction is the competition with gravitational radiation in compact triples, where the cycle drives the pericentre down and the radiation shrinks the orbit — a race whose outcome is the rate at which black hole binaries in dense environments merge.

About the same objects

Not linked from either essay — found by the objects both name.

What links here

Essays that link to this one from their own argument.

The objects this essay names

Each one links to every other essay that touches it.

Adiabatic invariantApsidal precessionEccentricityGeneral relativityHierarchical tripleHigh-eccentricity migrationHot jupiterInclinationKozai lidov mechanismSecular perturbationTidal circularisation