Astronomy & Celestial Mechanics
the sky, computed
Series

Kozai–Lidov — the series

6 essays on one idea, from the one that introduces it to the one that assumes the rest.
  1. An eccentricity and an inclination trading, at 65° of mutual tilt. The secular equations integrated from a nearly circular orbit (e = 0.02) inclined at 65° 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.8380 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.8380, 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.83 Kozai times. What the figure cannot show is what happens at the top of the cycle in a real system: at e = 0.838 the pericentre is 0.1620 of the semi-major axis, where tides, general relativity or a stellar surface all intervene, and the quadrupole picture ends.

    An inclination that turns into an eccentricity

    A distant companion cannot change an orbit's size or its energy. It can take a circular orbit tilted past 39.23 degrees and drive it to an eccentricity near one, and back, over and over — and the companion's mass and distance set only the clock.

    part 1 · orbits
  2. 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.

    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.

    part 2 · orbits
  3. A cycle and a radiation, racing for the same pericentre. Two integrations of the same compact triple, in units of the cycle's own timescale at the starting separation. The dashed curves are the semi-major axis and the solid ones the eccentricity. Radiation alone shrinks the orbit slowly and smoothly, reaching a tenth of its starting size after 1136 cycle times — which the closed form for a circular inspiral predicts to within a few per cent, and which is the check that the integrator is right. With the cycle running, the distant companion drives the eccentricity to 0.9104, the pericentre falls by two orders of magnitude, and the radiation — whose rate goes as the inverse seventh power of one minus the square of the eccentricity — empties the orbit in 18. That is a factor of 62, obtained without changing a single parameter of either mechanism. The two do not merely add: the cycle manufactures, at its own peak, exactly the configuration in which the radiation is most efficient, and the radiation then shrinks the orbit until the relativistic precession switches the cycle off.

    A cycle that builds what the wave needs

    Gravitational radiation empties an orbit at a rate that goes as the inverse seventh power of one minus the square of the eccentricity. A distant companion drives exactly that eccentricity, over and over. The two are not independent mechanisms that happen to act on the same binary — the first manufactures, at its own peak, the configuration the second is most efficient in.

    part 3 · orbits
  4. A ceiling the outer orbit imposes. The greatest eccentricity a cycle can drive an inner orbit to, against the mutual inclination it started at, for four ratios of the inner orbit's angular momentum to the outer one's. The bottom curve is the test-particle case and it is the familiar closed form, √(1 − (5/3)cos²i₀), reaching one at ninety degrees — recovered here by a solver that knows nothing about it, which is the check. Every other curve is lower. The reason is conservation: the cycle works by moving angular momentum out of the inner orbit, and it has nowhere to put it but the outer one. A test particle can give up all of its angular momentum because the outer orbit is infinite by comparison; an inner binary carrying 1 times as much as its companion cannot, because tilting the outer orbit by the required amount changes the mutual inclination and shuts the cycle down. At ninety degrees the ceiling falls from 1.000 to 0.866. The pericentre a mechanism can deliver is bounded by what the perturber can absorb, and the test-particle formula is an upper bound rather than a prediction.

    A ceiling the outer orbit imposes

    The closed form for how far a cycle can drive an eccentricity assumes the inner orbit is a test particle, so the angular momentum it gives up goes nowhere. A real inner binary of comparable mass has to be paid for: the outer orbit absorbs what the inner one loses, tilts, and changes the very inclination the cycle depends on.

    part 4 · orbits
  5. An inner orbit's inclination, with and without the octupole. The inclination of an inner orbit to the plane of a distant perturber against time, in units of the quadrupole timescale, starting at 65 degrees with a nearly circular orbit, for octupole strengths ε = 0 and 0.03. At ε = 0 — a circular outer orbit — the inclination oscillates in the regular Kozai–Lidov cycle and never crosses ninety degrees, because the angular momentum along the outer orbit's axis is conserved and cannot change sign. With the outer orbit eccentric, the inclination drifts from cycle to cycle and crosses ninety degrees, the orbit turning from prograde to retrograde, first at t ≈ 120, 2 times in all within the time shown. The integration conserves the secular energy to 7·10⁻⁶.

    An orbit that turns over

    A distant third body makes an inner orbit trade inclination for eccentricity in slow, regular cycles, and one quantity — the inner orbit's angular momentum along the outer orbit's axis — keeps the trade bounded and forbids the orbit from ever tipping past ninety degrees. That conservation holds only if the outer orbit is a circle. Make it eccentric and the quantity drifts, the cycles stop repeating, and at a moment of extreme eccentricity the inner orbit turns over, from prograde to retrograde, with its pericentre pressed a thousand times closer than the regular cycle allowed.

    part 5 · orbits
  6. How much of a population is driven how close. The fraction of an isotropic population of inner orbits whose Kozai–Lidov cycling brings 1 − e below a given value within 150 secular timescales, against that value, both logarithmic, for octupole strengths ε = 0, 0.01, 0.03; each curve is 128 integrated orbits with random pericentres and nodes. The thin line is the quadrupole closed form, √((3/5)(1 − (1 − q)²)), which the ε = 0 population follows. At 1 − e = 0.01 the fractions are 9.4, 20.3, 40.6 per cent; at 0.001, 3.1, 18.0, 38.3. The octupole matters most where the quadrupole is weakest — at the extreme closeness a planet needs to be captured by tides, or a pair of compact objects needs to merge — because it adds orbits whose quadrupole cycle stopped far short of it.

    Eleven in a hundred by geometry, and the rest by waiting

    A distant companion drives an inner orbit through cycles of eccentricity, and whether a given system reaches a dangerously close approach depends on its inclination. Counted over a population with random orientations, the answer at quadrupole order is pure geometry — about eleven per cent reach one per cent of their semi-major axis, all at the same inclination. An eccentric outer orbit multiplies that share, spreads the inclinations, and makes the count depend on how long the systems have been waiting.

    part 6 · orbits

All series