von Zeipel-Lidov-Kozai Migration

• von Zeipel-Lidov-Kozai migration is a secular mechanism where gravitational perturbations from a distant companion induce large oscillations in eccentricity and inclination, reducing the pericenter distance.
• The process couples conservative secular dynamics with dissipative mechanisms like tides or gravitational-wave emission to shrink and circularize orbits across varied astrophysical systems.
• Recent formulations extend classical quadrupole theory with octupole and Brown corrections to accurately capture non-linear effects and migration thresholds in both strongly and mildly hierarchical systems.

Searching arXiv for recent and foundational work on von Zeipel–Lidov–Kozai migration. von Zeipel–Lidov–Kozai migration is a secular high-eccentricity migration process in hierarchical or mildly hierarchical multiple systems in which gravitational perturbations from a distant companion drive large-amplitude, long-term exchanges between eccentricity and inclination in an inner orbit, reducing the pericenter distance until dissipative processes such as tides or gravitational-wave emission can shrink and circularize the orbit. In the conservative limit, the dynamics are migration only in eccentricity–inclination space, with semi-major axis fixed; true radial migration arises when dissipation is coupled to the secular dynamics (Beust et al., 2012, Bataille et al., 2018, Lei et al., 19 May 2025). Across planetary, stellar, compact-object, satellite, and small-body contexts, the mechanism is commonly described through libration or circulation of the argument of periapsis, conservation or approximate conservation of a Kozai-like angular-momentum component, and the existence of resonant phase-space structures that regulate access to extreme eccentricities (Beust et al., 2012, Lei et al., 2022, Lei et al., 2022).

1. Historical definition and dynamical scope

The modern term combines contributions associated with von Zeipel, Lidov, and Kozai and refers to the secular dynamics of a hierarchical triple in which a distant perturber excites coupled oscillations of eccentricity and inclination in an inner orbit (Marcos et al., 28 Jul 2025, Beust et al., 2012). In the simplest quadrupole, test-particle formulation, the secular Hamiltonian is integrable, the argument of periapsis $\omega$ may librate around $90^\circ$ or $270^\circ$, and the quantity

$\sqrt{1-e^2}\cos i$

is conserved, enforcing an anticorrelated exchange between $e$ and $i$ (Bataille et al., 2018, Beust et al., 2012, Muñoz et al., 2016). For initially circular inner orbits, the classical inclination window is

$$i \in \big[\arccos(\sqrt{3/5}),\,\arccos(-\sqrt{3/5})\big] \approx [39.2^\circ,140.8^\circ]$$

in the test-particle quadrupole limit (Bataille et al., 2018).

In migration problems, this conservative secular exchange is only the first stage. The second stage occurs when the maximum eccentricity becomes large enough that the pericenter

$q = a(1-e)$

or $q_1=a_1(1-e_1)$ becomes sufficiently small for tides, gravitational-wave emission, or related dissipative mechanisms to remove orbital energy during repeated pericenter passages (Beust et al., 2012, Bataille et al., 2018, Muñoz et al., 2016). The semi-major axis then shrinks while the orbit circularizes, producing hot Jupiters, close stellar binaries, or compact-object mergers depending on the physical setting (Beust et al., 2012, Bataille et al., 2018, Bhaskar et al., 18 Jun 2026).

The same secular architecture also appears outside the standard hierarchical binary-planet problem. Recent work treats low-hierarchy triples where the ZLK timescale approaches orbital periods (Lei et al., 19 May 2025), near-Earth asteroids affected by Jupiter (Marcos et al., 28 Jul 2025), resonant trans-Neptunian objects inside mean-motion resonances with Neptune (Lei et al., 2022), and white-dwarf planetary systems in which migration is delayed until after stellar evolution (Muñoz et al., 2020). This suggests that “von Zeipel–Lidov–Kozai migration” is best understood as a family of secular transport processes rather than a single specialized channel.

2. Conservative Hamiltonian structure

In the standard secular treatment of a hierarchical triple, the fast mean anomalies are averaged out and the dynamics reduce to slow evolution of orbital elements. For the inner orbit, a common canonical description uses Delaunay variables

$$L = \sqrt{\mu a},\quad G = L\sqrt{1-e^2},\quad H = G\cos i,$$

with conjugate angles $90^\circ$0, $90^\circ$1, and $90^\circ$2 (Lei et al., 19 May 2025). In the quadrupole test-particle problem, $90^\circ$3 is constant, so $90^\circ$4 is fixed, and axial symmetry makes $90^\circ$5 an integral, leaving a 1-degree-of-freedom Hamiltonian in $90^\circ$6 (Lei et al., 19 May 2025, Beust et al., 2012).

A representative quadrupole secular Hamiltonian is

$90^\circ$7

with the companion invariant

$90^\circ$8

or equivalently the “Kozai constant” $90^\circ$9 in stellar-triple notation (Beust et al., 2012, Bataille et al., 2018). These constants organize the familiar phase portraits in the $270^\circ$0 plane, separating circulation from libration islands (Beust et al., 2012, Bataille et al., 2018).

For low-hierarchy systems, the standard double-averaged quadrupole Hamiltonian can fail because nonlinear short-period effects contaminate the long-term evolution. A recent formulation extends the Brown Hamiltonian by carrying the von Zeipel procedure to second order in the disturbing function, producing a normalized long-term Hamiltonian

$270^\circ$1

where $270^\circ$2 is the standard double-averaged quadrupole term, $270^\circ$3 is the classical Brown Hamiltonian, and $270^\circ$4 is a new extended Brown term induced by inner-orbit fast oscillations (Lei et al., 19 May 2025). In element form,

$270^\circ$5

$270^\circ$6

with $270^\circ$7, and $270^\circ$8 adds higher harmonics such as $270^\circ$9 (Lei et al., 19 May 2025). The result remains integrable at quadrupole order in the restricted problem, but it modifies fixed points, resonant islands, precession rates, and the shape and period of ZLK cycles (Lei et al., 19 May 2025).

A complementary conceptual development identifies the eccentric ZLK effect at octupole order with apsidal secular resonance. In the restricted hierarchical planetary problem, after averaging over rotating ZLK cycles, one obtains an effective resonant Hamiltonian with slow angle

$\sqrt{1-e^2}\cos i$0

where $\sqrt{1-e^2}\cos i$1 is the apsidal angle difference. The resulting phase portraits contain multiple libration zones, separatrices, and flipping regions, and the paper concludes that the eccentric ZLK effect is dynamically equivalent to apsidal resonance in the test-particle limit (Lei et al., 2022). This rephrasing is especially relevant for migration because the extreme-eccentricity channels that enable dissipation are localized in specific resonant zones rather than being generic across phase space (Lei et al., 2022).

3. Secular timescales, hierarchy, and the onset of extreme eccentricity

A standard quadrupole ZLK timescale used in low-hierarchy analyses is

$\sqrt{1-e^2}\cos i$2

while in hierarchical stellar and planetary migration problems one often uses

$\sqrt{1-e^2}\cos i$3

or

$\sqrt{1-e^2}\cos i$4

depending on notation and context (Lei et al., 19 May 2025, Beust et al., 2012, Bhaskar et al., 18 Jun 2026). The common point is that ZLK is a secular process: much slower than orbital motion in strongly hierarchical triples, but not necessarily so in mildly hierarchical ones (Lei et al., 19 May 2025, Gao et al., 25 Jul 2025).

Recent low-hierarchy work introduces two correction parameters,

$\sqrt{1-e^2}\cos i$5

$\sqrt{1-e^2}\cos i$6

with $\sqrt{1-e^2}\cos i$7 essentially the single-averaging parameter and $\sqrt{1-e^2}\cos i$8 measuring how strongly inner-orbit fast variations contaminate the secular dynamics (Lei et al., 19 May 2025). For nearly circular outer orbits and $\sqrt{1-e^2}\cos i$9,

$e$0

so the hierarchy can be read directly from ratios of orbital and secular timescales (Lei et al., 19 May 2025). In the irregular-satellite examples treated there, $e$1 is comparable to $e$2 and only $e$3–$e$4 times $e$5, making the classical double average unreliable (Lei et al., 19 May 2025).

At higher order in the hierarchy expansion, the octupole parameter

$e$6

measures the strength of non-axisymmetric secular forcing (Lei et al., 19 May 2025). In hierarchical migration literature the corresponding quantity is often written

$e$7

or

$e$8

and controls orbit flips, very large eccentricity spikes, and chaotic modulation of the quadrupole ZLK cycles (Muñoz et al., 2016, Bhaskar et al., 18 Jun 2026). This suggests a practical separation of regimes: quadrupole-dominated ZLK for nearly circular outer companions and strong hierarchy, octupole-dominated or eccentric ZLK when $e$9 is significant, and Brown-corrected quasi-secular ZLK when $i$0 approaches orbital times (Lei et al., 19 May 2025, Gao et al., 25 Jul 2025).

4. Dissipative coupling and the mechanism of migration

The conservative Hamiltonian alone does not change the semi-major axis in the restricted secular problem; instead,

$i$1

in the absence of non-conservative forces (Lei et al., 19 May 2025). Migration therefore requires coupling the secular eccentricity pumping to a dissipative channel. In planetary systems this is usually equilibrium tides or, in some regimes, chaotic dynamical tides; in stellar triples it is tidal friction inside the inner stars; in compact-object problems it may be gravitational-wave emission (Beust et al., 2012, Bataille et al., 2018, Muñoz et al., 2020).

A standard high-eccentricity migration picture has two stages. First, ZLK cycles at large $i$2 drive the eccentricity toward $i$3, periodically reducing $i$4 while leaving $i$5 nearly fixed (Beust et al., 2012, Bataille et al., 2018). Second, dissipation concentrated near pericenter removes orbital energy. Since the orbital angular momentum is approximately conserved during this stage, the post-circularization orbit satisfies relations such as

$i$6

or, in the white-dwarf planetary case,

$i$7

depending on the adopted limit and notation (Bhaskar et al., 18 Jun 2026, Muñoz et al., 2020).

The migration threshold is commonly expressed as a required minimum pericenter. For hot-Jupiter production by ZLK in stellar binaries, a critical migration pericenter is

$i$8

with the corresponding migration eccentricity

$i$9

(Muñoz et al., 2016). The disruption threshold is similarly written

$$i \in \big[\arccos(\sqrt{3/5}),\,\arccos(-\sqrt{3/5})\big] \approx [39.2^\circ,140.8^\circ]$$0

so successful migration requires reaching $$i \in \big[\arccos(\sqrt{3/5}),\,\arccos(-\sqrt{3/5})\big] \approx [39.2^\circ,140.8^\circ]$$1 without exceeding $$i \in \big[\arccos(\sqrt{3/5}),\,\arccos(-\sqrt{3/5})\big] \approx [39.2^\circ,140.8^\circ]$$2 (Muñoz et al., 2016, Muñoz et al., 2020).

Short-range forces determine whether the required eccentricities are dynamically accessible. In quadrupole ZLK with short-range forces, the limiting eccentricity $$i \in \big[\arccos(\sqrt{3/5}),\,\arccos(-\sqrt{3/5})\big] \approx [39.2^\circ,140.8^\circ]$$3 is set by balancing the secular torque against apsidal precession from GR, tidal bulges, and rotation (Muñoz et al., 2016, Muñoz et al., 2020). In low-hierarchy or high-eccentricity planetary migration, a similar role is played by GR precession

$$i \in \big[\arccos(\sqrt{3/5}),\,\arccos(-\sqrt{3/5})\big] \approx [39.2^\circ,140.8^\circ]$$4

and tidal/rotational precession; once these dominate over the ZLK torque, the resonance is quenched and the orbit circularizes along a nearly constant-angular-momentum path (Beust et al., 2012, Lu et al., 2024).

The HAT-P-11 analysis shows that structural evolution of the migrating planet can be dynamically essential. There, the orbit-averaged tidal luminosity

$$i \in \big[\arccos(\sqrt{3/5}),\,\arccos(-\sqrt{3/5})\big] \approx [39.2^\circ,140.8^\circ]$$5

is used to update the planetary radius through a fitted luminosity–radius relation

$$i \in \big[\arccos(\sqrt{3/5}),\,\arccos(-\sqrt{3/5})\big] \approx [39.2^\circ,140.8^\circ]$$6

with the consequence that a larger $$i \in \big[\arccos(\sqrt{3/5}),\,\arccos(-\sqrt{3/5})\big] \approx [39.2^\circ,140.8^\circ]$$7 greatly enhances tidal precession because $$i \in \big[\arccos(\sqrt{3/5}),\,\arccos(-\sqrt{3/5})\big] \approx [39.2^\circ,140.8^\circ]$$8 (Lu et al., 2024). In that system, significant tidally driven radius inflation is required for ZLK cycles to be quenched at the semi-major axis that later circularizes to the observed orbit; fixed-radius integrations over-migrate the planet inward (Lu et al., 2024). This suggests that in some sub-Jovian migration problems, internal thermal response is part of the secular dynamics rather than a secondary correction.

5. Regime-dependent formulations and representative astrophysical settings

A major recent theme is that the secular backbone depends strongly on hierarchy and on whether the perturber is stellar, planetary, or embedded in a more complex multi-body architecture. The table summarizes regimes that are explicitly developed in recent work.

Regime	Representative formulation	Main consequence for migration
Strongly hierarchical triples	Double-averaged quadrupole or octupole secular Hamiltonian (Beust et al., 2012, Bataille et al., 2018, Muñoz et al., 2016)	Classical ZLK cycles; migration set by $$i \in \big[\arccos(\sqrt{3/5}),\,\arccos(-\sqrt{3/5})\big] \approx [39.2^\circ,140.8^\circ]$$9, GR, and tides
Mildly hierarchical triples	Brown-corrected secular Hamiltonians (Lei et al., 19 May 2025, Gao et al., 25 Jul 2025)	Periods, amplitudes, and flip boundaries can be misestimated by DA models
Planet–planet high-e migration	Octupole-order double-averaged Hamiltonian with tides and spin evolution (Bhaskar et al., 18 Jun 2026)	Short-period hot Jupiters preferentially arise from ZLK with inclined companions
Secular resonance plus ZLK in small-body dynamics	Direct $q = a(1-e)$0-body or semi-secular resonant models (Marcos et al., 28 Jul 2025, Lei et al., 2022)	Migration is mainly in perihelion distance or resonant action rather than circularization
White-dwarf systems	ZLK plus tides after stellar-evolution-driven orbital expansion (Muñoz et al., 2020)	Delayed migration becomes possible only in the WD phase

In stellar triples, one concrete application is the production of short-period inner binaries through ZLK cycles plus tidal friction. In a secular octupole-plus-GR-plus-tides framework, migrating systems preferentially have small $q = a(1-e)$1, where

$q = a(1-e)$2

and are not determined by inclination alone: the initial $q = a(1-e)$3 and $q = a(1-e)$4 also control whether the system reaches the large $q = a(1-e)$5 required for migration (Bataille et al., 2018). A key result is that systems starting near the center of a KL libration island can remain non-migrating because the eccentricity amplitude is too small, even though they are formally in the ZLK regime (Bataille et al., 2018).

For exoplanet migration, Gliese 436b provides a classic two-phase example. Starting near its current orbit, GR and tides suppress ZLK and the planet simply circularizes; starting $q = a(1-e)$6–$q = a(1-e)$7 times farther out allows active Kozai cycles, followed by a second phase in which the system leaves the resonance and rapidly shrinks at high eccentricity (Beust et al., 2012). The paper emphasizes that the total effective circularization time can exceed the standard constant-$q = a(1-e)$8 estimate by a factor larger than $q = a(1-e)$9, making it possible to observe the planet in an intermediate migration stage (Beust et al., 2012).

In planet–planet secular high-e migration, recent simulations identify two limiting channels: mutually inclined companions driving ZLK and nearly coplanar eccentric companions driving CHEM-like migration. The shortest-period hot Jupiters arise preferentially from ZLK with highly inclined companions and display a broad range of stellar obliquities, whereas longer-period hot Jupiters are produced over longer timescales by nearly coplanar high-e migration and remain aligned (Bhaskar et al., 18 Jun 2026). The paper derives an analytic envelope relating final mutual inclination and final semi-major axis,

$q_1=a_1(1-e_1)$0

which organizes the simulated obliquity–period anti-correlation (Bhaskar et al., 18 Jun 2026).

For near-Earth asteroids, the same label is used in a somewhat different sense. In the Miorita study, a von Zeipel–Lidov–Kozai secular resonance with Jupiter is identified numerically through long-period anticorrelated $q_1=a_1(1-e_1)$1–$q_1=a_1(1-e_1)$2 oscillations of period $q_1=a_1(1-e_1)$3 kyr, while a simultaneous apsidal resonance $q_1=a_1(1-e_1)$4 confines the semimajor axis and enables metastable low-perihelion states (Marcos et al., 28 Jul 2025). There, “migration” is chiefly migration of perihelion distance, not of semimajor axis, and some objects are ultimately driven into the Sun (Marcos et al., 28 Jul 2025).

6. Efficiency, limitations, and current controversies

The modern literature treats von Zeipel–Lidov–Kozai migration as robust but conditional. A central result from analytic work on giant planets in stellar binaries is that the total migration fraction

$q_1=a_1(1-e_1)$5

is much less sensitive to planet mass, host type, and tidal strength than the surviving hot-Jupiter fraction $q_1=a_1(1-e_1)$6, because the main geometric determinant is the fraction of initial inclinations that can reach the migration pericenter (Muñoz et al., 2016). For fiducial giant-planet populations at $q_1=a_1(1-e_1)$7 AU, the analytic calculations give $q_1=a_1(1-e_1)$8, while $q_1=a_1(1-e_1)$9 is only a few percent or less depending on mass and dissipation (Muñoz et al., 2016).

Population synthesis of stellar-binary KL migration yields a stronger negative conclusion. One large Monte Carlo study finds that KL migration produces hot Jupiters with semi-major axes generally smaller than observed and can explain the data only if both of the following hold: tidal dissipation at high eccentricities is at least $$L = \sqrt{\mu a},\quad G = L\sqrt{1-e^2},\quad H = G\cos i,$$0 times more efficient than the upper limit inferred from the Jupiter–Io interaction, and highly eccentric planets are tidally disrupted at distances $$L = \sqrt{\mu a},\quad G = L\sqrt{1-e^2},\quad H = G\cos i,$$1 AU (Petrovich, 2014). On the basis of occurrence rate and semi-major-axis distribution, that work argues that KL migration in stellar binaries can produce at most $$L = \sqrt{\mu a},\quad G = L\sqrt{1-e^2},\quad H = G\cos i,$$2 of observed hot Jupiters and almost no intermediate-period planets with $$L = \sqrt{\mu a},\quad G = L\sqrt{1-e^2},\quad H = G\cos i,$$3–$$L = \sqrt{\mu a},\quad G = L\sqrt{1-e^2},\quad H = G\cos i,$$4 AU (Petrovich, 2014). This remains one of the clearest statements of the channel’s limitations.

A related limitation is that low-hierarchy corrections can qualitatively change flip criteria and phase-space geometry. In mildly hierarchical triples, Brown corrections break the symmetry of flipping regions across $$L = \sqrt{\mu a},\quad G = L\sqrt{1-e^2},\quad H = G\cos i,$$5, and double-averaged quadrupole-plus-octupole models can predict flips where $$L = \sqrt{\mu a},\quad G = L\sqrt{1-e^2},\quad H = G\cos i,$$6-body integrations show none, or miss flips that do occur (Gao et al., 25 Jul 2025). The extended pendulum approximation developed there shows that Brown terms shift the resonant center away from $$L = \sqrt{\mu a},\quad G = L\sqrt{1-e^2},\quad H = G\cos i,$$7, so the flip region is no longer symmetric in prograde and retrograde initial conditions (Gao et al., 25 Jul 2025). For migration predictions that depend sensitively on the extreme-eccentricity tail, this is not a small correction.

The white-dwarf case adds a different limitation: the ZLK route is viable only after imposing survival through RGB and AGB evolution. For WD 1856b, the allowed planet masses and initial semimajor axes are constrained simultaneously by the requirement of survival during giant-star phases and the requirement of later WD-phase ZLK migration, leaving a narrow range $$L = \sqrt{\mu a},\quad G = L\sqrt{1-e^2},\quad H = G\cos i,$$8–$$L = \sqrt{\mu a},\quad G = L\sqrt{1-e^2},\quad H = G\cos i,$$9 and main-sequence semimajor axis $90^\circ$00–$90^\circ$01 AU (Muñoz et al., 2020). The inferred occurrence rate of Kozai-migrated planets around white dwarfs is only $90^\circ$02 (Muñoz et al., 2020).

The small-body literature raises a different interpretive issue. In the NEA and trans-Neptunian settings, ZLK migration often refers to inward or outward transport in perihelion distance or resonant action rather than tidal circularization (Marcos et al., 28 Jul 2025, Lei et al., 2022). This suggests that the phrase “ZLK migration” is used in at least two technically distinct ways: strict high-eccentricity migration producing compact circularized orbits, and secular transport in $90^\circ$03–$90^\circ$04 space that changes encounter geometry or perihelion without large semimajor-axis decay. A plausible implication is that context-specific clarification is necessary whenever the term is used outside the exoplanet and stellar-binary literature.

7. Prospects and open directions

Several current developments indicate where the subject is moving. First, the conservative backbone is being improved in the quasi-secular regime. The extended Brown Hamiltonian expresses both the Hamiltonian and the mean–osculating transformation in closed form with respect to the inner and outer eccentricities, avoiding small-$90^\circ$05 expansions and the Laplace-limit restriction that affects Hansen-series approaches (Lei et al., 19 May 2025). Because dissipation is applied to osculating orbits, this mean–osculating map is particularly relevant when the dissipation rate depends sensitively on pericenter distance (Lei et al., 19 May 2025).

Second, the octupole and Brown frameworks are beginning to converge. The mildly hierarchical eccentric-ZLK study extends Brown corrections to include octupole–octupole coupling, shows that quadrupole–octupole coupling vanishes in the chosen gauge, and demonstrates that Brown corrections are the key factor responsible for breaking the symmetry of flipping regions (Gao et al., 25 Jul 2025). This suggests a route toward secular migration calculations that remain accurate when $90^\circ$06, a regime relevant for planets in stellar binaries and for binaries near supermassive black holes (Gao et al., 25 Jul 2025).

Third, obliquity has become a diagnostic of migration channel. In planet–planet secular migration, the shortest-period hot Jupiters tend to arise from ZLK driven by highly inclined companions and thus preserve a broad range of stellar obliquities, while longer-period systems are preferentially aligned because they arise from nearly coplanar secular migration (Bhaskar et al., 18 Jun 2026). The paper predicts that the shortest-period hot Jupiters should have distant planetary companions with broadly distributed mutual inclinations, whereas the companions of longer-period systems should be nearly coplanar (Bhaskar et al., 18 Jun 2026). Gaia astrometry is explicitly identified there as a key observational test.

Finally, some systems now require coupled dynamical and thermodynamic modeling. The HAT-P-11 study shows that scattering, ZLK cycles, tidal dissipation, and radius inflation must be evolved together to reproduce the observed semimajor axis and eccentricity (Lu et al., 2024). This suggests that a purely orbital secular model may be insufficient for sub-Saturn and Neptune-mass planets undergoing high-e migration, especially when tidal luminosity can strongly alter the planetary radius and therefore the precession rate (Lu et al., 2024).

Taken together, the literature supports an overview in which von Zeipel–Lidov–Kozai migration is a secular transport mechanism whose conservative core is now understood at several levels of refinement—quadrupole, octupole, apsidal-resonant, Brown-corrected, and resonant-in-MMR—while the actual astrophysical outcome depends on how that core is coupled to GR, tides, spins, internal structure, stellar evolution, or additional resonances (Lei et al., 19 May 2025, Lei et al., 2022, Bataille et al., 2018, Bhaskar et al., 18 Jun 2026). This suggests that the central problem is no longer whether ZLK migration exists, but rather which specific formulation is appropriate for a given hierarchy, and which subset of systems can traverse the narrow corridor between insufficient eccentricity excitation and catastrophic disruption.