Skip to main content
Sandbox Physics
Kozai–Lidov cycles from torque to planetary orbits

AC09 · Paper guide

Read octupole dynamics, relativity and a proposed planet history

Sandbox Physics ·

The next step beyond a clean eccentricity cycle is to ask which approximation allowed it, which observations constrain a real system, and which added physics can change the outcome. Three fixed papers connect those questions through octupole dynamics, a highly eccentric planet and fast changes beyond double averaging.

Make a small model ledger before opening the figures

Write down four choices for every calculation: massless or finite inner body, circular or eccentric companion, averaged or resolved orbital phases, conservative or dissipative extra forces. This prevents a visually similar eccentricity curve from carrying an inapplicable conservation law.

Begin with the Lab’s circular-outer quadrupole model. Its axial angular momentum and averaged Hamiltonian are conserved. The full Newtonian reference keeps higher spatial orders and fast phases. The instantaneous-quadrupole control retains the fast phases while matching the spatial truncation. Keep those three rows available when reading the papers.

For each figure, identify what the axes measure, what varies across curves and what was held fixed. If the vertical axis measures the small distance from unit eccentricity on a logarithmic scale, a modest vertical change can represent a large change in periapsis. Convert one plotted point into a nearest distance before comparing physical outcomes.

q=a(1−e),ϵoct=m0−m1m0+m1a1a2e21−e22q=a(1-e),\qquad \epsilon_{\mathrm{oct}}=\frac{m_0-m_1}{m_0+m_1}\frac{a_1}{a_2}\frac{e_2}{1-e_2^2}
The octupole coefficient vanishes for a circular outer orbit or equal inner masses. For our test particle, the vanishing outer eccentricity sets this coefficient to zero.

Reading: Liu, Muñoz & Lai, Short-range forces in triple systems, arXiv:1409.6717v2 · Klein & Hamilton, Short-range forces beyond double averaging, arXiv:2606.19456v1

Liu, Muñoz and Lai 2015, fixed version two

Read equations 7 and 8 together: they separate quadrupole and octupole potentials. Then inspect the mass and outer-eccentricity factors in the octupole strength. Mark which symmetry from the circular-outer problem can now be lost.

Equations 32–33 define the averaged relativistic potential and its strength; equation 42 gives its periapsis-precession rate. Compare their normalization with L2 before inserting a coefficient. Figures 3–4 and 7 then show how additional short-range effects change extreme eccentricity and flips in the paper’s stated parameter families.

A useful reading task is to choose one curve and list every extra precession included. Relativistic, tidal-bulge and rotational-bulge precession are conservative contributions in this discussion; a tidal-migration history also needs a dissipative prescription. Keep that distinction when carrying a limiting eccentricity into a planet-formation argument.

Reading: Liu, Muñoz & Lai, Short-range forces in triple systems, arXiv:1409.6717v2

Gupta and colleagues 2024, observed planet and inferred history

Figures 1 and 2 show the transit and radial-velocity evidence for TIC 241249530 b. The reported eccentricity is about 0.9412, and the projected stellar spin–orbit angle is about 163.5 degrees. That projected angle is not the mutual inclination between the planet and the distant stellar orbit.

Figure 3 is a population diagram. For the proposed migration history, turn instead to Extended Data Figures 5 and 6. The latter is an illustrative secular evolution with octupole dynamics and an equilibrium-tide prescription. Its tidal parameters were adjusted to reach the present system near the adopted age.

The companion’s projected separation is not a precisely measured orbital semimajor axis. Write separate columns for observed quantities, uncertain orbital geometry and chosen dissipation parameters. The evolutionary scenario depends on those chosen inputs, while the observed orbit constrains its endpoint.

Reading: Gupta et al., A hot-Jupiter progenitor on a super-eccentric retrograde orbit, Nature 632 (2024)

Klein and Hamilton 2026, fixed preprint version one

The June 17 version is a preprint, and a later revision exists. Figure 1 separates double versus single averaging and inclusion versus omission of short-range forces. Read that map before treating the paper’s title as a universal statement about relativity.

Figures 2–3 follow changes between cycles when fast outer motion interacts with short-range precession. Effective quantities that would stay fixed under double averaging can jump during extreme episodes. Appendix B compares approximations with direct integrations; Appendix A supplies the small eccentricity seed behind shorthand circular initial conditions.

This is a useful challenge to an overgeneralization from L2. The Lab demonstrates suppression within its conservative doubly averaged model. It does not calculate the single-averaged relativistic examples in this paper. Reproducing them would require a new matched fast-phase force model and its own tests.

Reading: Klein & Hamilton, Short-range forces beyond double averaging, arXiv:2606.19456v1

Turn one claim into a reproducible question

Choose a claim narrow enough to test: for example, whether a cycle maximum changes when only averaging is altered at fixed quadrupole order. State the initial state, force terms, output cadence and stop rule before choosing an integrator. Preserve positions and velocities so another reader can define a different eccentricity or period statistic.

For a real planet, begin with the likelihood or published uncertainty of the measured orbit. Sample the unmeasured companion orientation and semimajor axis as assumptions, then compare several plausible dissipation prescriptions. A history that reaches the present state is an existence example; estimating how common it is requires the distribution of those inputs.

The public frozen experiment is a compact place to practice the first calculation. It provides an explicit force control and full state records. The observed system and frontier calculations offer the next questions once those distinctions are secure.

Reading: REBOUND 5.2.2 official Kozai example and IAS15 integrator · Gupta et al., A hot-Jupiter progenitor on a super-eccentric retrograde orbit, Nature 632 (2024)

Try the orbit exchange

  1. Use the instantaneous-quadrupole tutorial control to practice changing averaging while keeping spatial order fixed.
  2. Open the relativistic preset and state precisely which curve includes precession. Then locate the extra averaging choice required for the 2026 calculation.

Exercises with explanations

  1. Why does a circular outer orbit remove the octupole term even when the inner planet is massless?

    Read the explanation

    The outer-eccentricity factor is zero. The test-particle mass asymmetry remains, but multiplying it by zero removes the octupole coefficient in this hierarchical expansion.

  2. Which figure should be used to read the proposed migration track of TIC 241249530 b?

    Read the explanation

    Extended Data Figure 6. Main Figure 3 is the population diagram, and Figures 1–2 contain the transit and velocity evidence. Separate those roles before interpreting a historical track.

  3. Does an averaged relativistic suppression experiment contradict the 2026 preprint?

    Read the explanation

    The model scopes differ. The preprint investigates behavior after relaxing double averaging. A matched comparison must include the relevant fast-phase terms and the same short-range force model.

Sources and calculation scope

  1. Liu, Muñoz & Lai, Short-range forces in triple systems, arXiv:1409.6717v2

    Fixed 7 January 2015 version, DOI 10.1093/mnras/stu2396. Sections 2–4, quadrupole potential 7, octupole 8, GR 32–33/42 and Figures 3–4/7. The browser retains the circular-outer, test-particle quadrupole; the guide reads the wider model.

  2. REBOUND 5.2.2 official Kozai example and IAS15 integrator

    The official example motivates an independent direct reference. Frozen trajectories use REBOUND 5.2.2, source commit c34134716cd484ed9cdb9f7d54ef8d859bdf547d, IAS15 accuracy parameter one hundred-billionth, with separately stated initial conditions and stopping rules. REBOUND software is GPL-3.0-only; no publisher code is copied into the Lab.

  3. Gupta et al., A hot-Jupiter progenitor on a super-eccentric retrograde orbit, Nature 632 (2024)

    Version of record 17 July 2024, DOI 10.1038/s41586-024-07688-3, CC BY 4.0. Observed transits/RVs are Figures 1–2; Extended Data 5–6 give model-dependent initial conditions and an illustrative tidal-migration history. Full publisher-deposited XML, PMC11291287, was read for these distinctions.

  4. Klein & Hamilton, Short-range forces beyond double averaging, arXiv:2606.19456v1

    Fixed 17 June 2026 preprint version; an August v2 also exists. Figure 1 separates averaging/force choices, Figures 2–3 show cycle changes, and Appendix B compares approximations with direct N-body examples. This is a reading guide; the browser does not reproduce its single-averaged GR/tidal calculations.

FeedbackSomething broken, confusing, or missing on this page?Tell us
What kind of feedback?

Sent with this page’s address, its Lab ID, the site language and your device type. Please leave out personal details. Feedback is kept for up to 12 months.Privacy policy