High-Eccentricity Migration in Exoplanets

• HEM is an exoplanet orbital evolution mechanism where a planet on a highly eccentric orbit is gradually circularized by tidal dissipation, reducing its semimajor axis.
• It involves dynamical processes such as Lidov–Kozai oscillations, coplanar interactions, and chaotic tide effects that collectively sculpt the final orbital configuration.
• HEM leaves distinct population-level imprints in the period–mass, period–radius, and period–density distributions, providing empirical diagnostics for close-in giant and sub-Neptune planets.

Searching arXiv for recent HEM-related papers to ground the response. High-eccentricity migration (HEM) is an exoplanet-orbital evolution scenario in which a planet reaches the vicinity of its host star on a highly eccentric orbit and is then circularized, with a concomitant decrease in semimajor axis, by tidal dissipation during repeated close periastron passages (Giacalone et al., 2017). In the canonical picture, the eccentricity is excited by secular or scattering-driven dynamics, while the final short-period configuration is set by the interplay between angular-momentum conservation, tidal circularization, and tidal survival against Roche-limit or disruption constraints (Jackson et al., 2022, Owen et al., 2018). Although HEM was developed primarily to explain hot Jupiters, the same framework has been applied to warm-to-hot giant migration, Neptune-mass and sub-Saturn populations, and even rocky bodies around polluted white dwarfs (Dawson et al., 2021, Castro-González et al., 17 Apr 2026, O'Connor et al., 2020).

1. Core dynamical framework

The basic orbital geometry of HEM is set by approximately conserved orbital angular momentum during the dissipative phase. Two relations recur throughout the literature:

$a_{\rm final}=a(1-e^2),$

and, in the limit $e \simeq 1$,

$a_{\rm F}\simeq 2r_{\rm p}.$

These relations encode the fact that orbital energy is removed much more efficiently than angular momentum during close passages, so a highly eccentric orbit circularizes to a much smaller radius while retaining nearly the same angular-momentum content (Jackson et al., 2022, Owen et al., 2018).

This orbital mapping is physically meaningful only if the planet survives the periastron passages that make tides effective. A widely used survival scale is the tidal or Roche-like radius

$$r_{\rm tide}=\eta\,R_{\rm p}\left(\frac{M_\star}{M_{\rm p}}\right)^{1/3},$$

with survival commonly expressed as

$a_{\rm F}\ge 2\,r_{\rm tide}.$

In this formulation, HEM is not merely a migration mechanism; it is a selective filter that allows only a subset of highly eccentric trajectories to circularize without tidal stripping or disruption (Owen et al., 2018, Castro-González et al., 17 Apr 2026).

A central implication is that HEM predicts structured, not featureless, close-in exoplanet populations. It generically links the final period distribution to periastron survival, tidal efficiency, and the secular pathways that create extreme eccentricity in the first place. This suggests that HEM is best viewed as a family of excitation-plus-dissipation channels rather than as a single dynamical recipe.

2. Secular excitation channels and architectures

The most familiar HEM driver is Lidov–Kozai excitation by a distant inclined companion. In this channel, a stellar or planetary perturber drives coupled eccentricity–inclination oscillations, and tides become effective when the inner orbit attains sufficiently small periastron (Vick et al., 2018). Variants of this mechanism become especially rich in higher-order hierarchies. In stellar triples, 2+2 and 3+1 configurations increase hot-Jupiter formation efficiency relative to stellar binaries, but only by at most a few tens of per cent; the resulting hot Jupiters remain very similar to those formed in binaries, and no significant number of warm Jupiters is produced (Hamers, 2017). In the white-dwarf system WD1856+534b, a 2+2 secular inclination resonance broadens the extreme-LK window and permits $e \gtrsim 0.999$ over a wide range of initial inclinations, with an absolute pre-migration constraint $$a_1 \gtrsim 8\,{\rm AU}\,(a_3/1500\,{\rm AU})^{6/7}$$ and a likely successful range $30\,{\rm AU} \lesssim a_1 \lesssim 60\,{\rm AU}$ just before white-dwarf formation (O'Connor et al., 2020).

HEM does not, however, require large mutual inclinations. Coplanar High-eccentricity Migration (CHEM) is a secular two-planet channel operating for relatively low mutual inclinations, $i_{\rm tot}\lesssim 20^\circ$, in which an eccentric outer giant planet drives the inner orbit to extreme eccentricity. Two representative threshold families are emphasized: an initially circular inner orbit with $e_{\rm out}\gtrsim 0.67$ when $e \simeq 1$0, and a two-eccentric-planet configuration with both eccentricities $e \simeq 1$1 when the same control parameter is $e \simeq 1$2 (Petrovich, 2014). General relativity and the tidal quadrupole then cap the maximum eccentricity, setting a minimum pericenter of about $e \simeq 1$3–$e \simeq 1$4 AU and a minimum final hot-Jupiter semimajor axis of about $e \simeq 1$5 AU (Petrovich, 2014).

Warm-origin or short-distance HEM introduces an additional dynamical regime. When migration begins from a warm orbit rather than from beyond the ice line, general relativistic precession can reduce Kozai–Lidov amplitudes. In that regime, planets that do achieve HEM tend to end with near-polar spin-orbit alignments, $e \simeq 1$6–130 degrees, rather than with the more familiar bimodal concentrations at $e \simeq 1$7 and $e \simeq 1$8 degrees (Dawson et al., 2021). A related population-synthesis result is obtained when primordial stellar obliquity is generated during disk dispersal in binaries: ZLK-driven HEM then yields a predominantly retrograde obliquity distribution with a broad peak near $e \simeq 1$9, rather than the standard bimodal outcome (Vick et al., 2022).

3. Tidal dissipation, circularization, and disruption physics

The simplest HEM implementations use equilibrium tides. In this regime the circularization timescale is often parameterized by a modified or reduced tidal quality factor, with several population-level inferences clustering around $a_{\rm F}\simeq 2r_{\rm p}.$0 for giant planets in the circularization zone (Giacalone et al., 2017). An empirical calibration based on the eccentricity distribution of more than 500 Jovian-mass planets yields $a_{\rm F}\simeq 2r_{\rm p}.$1, where the reduced quantity $a_{\rm F}\simeq 2r_{\rm p}.$2 is denoted simply as $a_{\rm F}\simeq 2r_{\rm p}.$3 for convenience (Kawai et al., 19 Sep 2025). In this equilibrium-tide framework, tidal circularization must compete against disruption at small periastron and against insufficient damping at larger periastron.

A major development has been the inclusion of dynamical tides, especially tidally excited planetary $a_{\rm F}\simeq 2r_{\rm p}.$4-modes. In the chaotic-tide picture, successive periastron passages kick the mode, and the mode energy can grow diffusively if the phase change between passages becomes order unity. The usual criterion is

$a_{\rm F}\simeq 2r_{\rm p}.$5

with chaotic growth typically requiring $a_{\rm F}\simeq 2r_{\rm p}.$6 and $a_{\rm F}\simeq 2r_{\rm p}.$7–$a_{\rm F}\simeq 2r_{\rm p}.$8 for the standard giant-planet models considered (Vick et al., 2018). This mechanism creates transient very eccentric warm Jupiters on short timescales of a few to 100 Myrs, which then circularize efficiently into hot Jupiters because their periastra remain persistently small (Vick et al., 2018).

Nonlinear mode coupling strengthens this picture. Interactions between a parent $a_{\rm F}\simeq 2r_{\rm p}.$9-mode and daughter $$r_{\rm tide}=\eta\,R_{\rm p}\left(\frac{M_\star}{M_{\rm p}}\right)^{1/3},$$0- and $$r_{\rm tide}=\eta\,R_{\rm p}\left(\frac{M_\star}{M_{\rm p}}\right)^{1/3},$$1-modes induce an energy-dependent parent-frequency shift that randomizes the phase at pericenter and lowers the one-kick threshold for diffusive growth by about a factor of 5 relative to linear theory (Yu et al., 2021). This materially enlarges the portion of HEM parameter space in which rapid circularization is possible and was proposed as a way to mitigate the discrepancy between observed and predicted occurrence rates of close-in gas giants (Yu et al., 2021).

The same tidal physics can also destabilize the migrating planet’s internal structure. A coupled structure-and-orbital model for ridge Neptunes finds that circularizing a Neptune into the $$r_{\rm tide}=\eta\,R_{\rm p}\left(\frac{M_\star}{M_{\rm p}}\right)^{1/3},$$2–$$r_{\rm tide}=\eta\,R_{\rm p}\left(\frac{M_\star}{M_{\rm p}}\right)^{1/3},$$3 day ridge is difficult without runaway tidal inflation and likely atmospheric destruction if the tidal heat is deposited below the radiative-convective boundary (Hallatt et al., 26 Sep 2025). In those calculations, low-eccentricity ridge Neptunes can be emplaced by HEM only if they are metal-rich and have finely tuned tidal quality factors, or if dissipation operates in the upper reaches of the planet rather than in the deep interior (Hallatt et al., 26 Sep 2025).

4. Observational diagnostics and empirical tests

The cleanest direct test of standard steady-state HEM is the predicted existence of super-eccentric proto-hot Jupiters. If all hot Jupiters in the Kepler sample formed through standard HEM, the expected number of super-eccentric Jupiters is $$r_{\rm tide}=\eta\,R_{\rm p}\left(\frac{M_\star}{M_{\rm p}}\right)^{1/3},$$4, whereas the inferred observed number is $$r_{\rm tide}=\eta\,R_{\rm p}\left(\frac{M_\star}{M_{\rm p}}\right)^{1/3},$$5 (Jackson et al., 2022). The mismatch reaches 94.3% confidence, and the corresponding upper limit implies that HEM can account for at most about 62% of the hot Jupiters discovered by Kepler (Jackson et al., 2022). This does not refute HEM as a channel, but it does constrain its universality in its simplest equilibrium-tide, steady-state form.

A second diagnostic is the spatial eccentricity gradient across the circularization boundary. Explicit tidal-evolution models and backward integrations of observed systems show that circularization generally occurs outside the inner edge of the surviving distribution, while the observed period–mass and period–eccentricity planes display the predicted transition from mostly eccentric orbits at longer periods to mostly circular orbits at shorter periods (Giacalone et al., 2017). The same analysis finds tentative evidence that this gradient may extend to lower-mass planets, indicating that HEM may be relevant down to Neptune scales (Giacalone et al., 2017).

A third empirical signature is the distribution of orbital distance in units of the Roche limit. In a homogeneous radial-velocity analysis of 231 transiting giant planets, the distribution of $$r_{\rm tide}=\eta\,R_{\rm p}\left(\frac{M_\star}{M_{\rm p}}\right)^{1/3},$$6 for well-determined circular orbits peaks at 2.5, broadly consistent with the HEM expectation that circularized planets park near $$r_{\rm tide}=\eta\,R_{\rm p}\left(\frac{M_\star}{M_{\rm p}}\right)^{1/3},$$7 (Bonomo et al., 2017). The same study notes that the few circular planets with $$r_{\rm tide}=\eta\,R_{\rm p}\left(\frac{M_\star}{M_{\rm p}}\right)^{1/3},$$8 may instead have migrated through disc–planet interactions (Bonomo et al., 2017).

Spin-orbit geometry is informative but not decisive. Obliquity is often used as an HEM diagnostic because high obliquity frequently accompanies dynamical excitation, yet low obliquity does not rule out HEM because of coplanar channels and possible tidal realignment (Kawai et al., 19 Sep 2025). Kepler-1656 b illustrates the ambiguity: its sky-projected obliquity is consistent with alignment but also with moderate misalignment, with a most likely solution $$r_{\rm tide}=\eta\,R_{\rm p}\left(\frac{M_\star}{M_{\rm p}}\right)^{1/3},$$9 deg and $a_{\rm F}\ge 2\,r_{\rm tide}.$0 at 95% confidence (Rubenzahl et al., 2024). In one interpretation, these properties are compatible with CHEM if the mutual inclination is $a_{\rm F}\ge 2\,r_{\rm tide}.$1; in another, if the system is not relatively coplanar, the current state may be a rare snapshot of long-lived eccentricity oscillations that do not induce migration (Rubenzahl et al., 2024). An even stronger caution is provided by dynamical modeling of the same system, which argues that despite $a_{\rm F}\ge 2\,r_{\rm tide}.$2, Kepler-1656 b is not currently undergoing HEM and likely has not undergone HEM at all, but instead is a non-migrating, high-eccentricity survivor of gentle EKL forcing (Angelo et al., 2022). This directly addresses the common misconception that extreme eccentricity alone implies active migration.

5. Population-level imprints in the period–mass, period–radius, and period–density planes

HEM has been invoked not only for individual hot Jupiters but also for global close-in population structure. One prominent example is the sub-Jovian desert. In a two-process interpretation, photoevaporation shapes the low-mass or small-radius boundary, while HEM plus the tidal disruption barrier shapes the high-mass or large-radius boundary (Owen et al., 2018). The triangular geometry follows because photoevaporation becomes more effective inward, whereas more massive gas giants can survive closer in during HEM because their tidal radii are smaller (Owen et al., 2018).

A more recent extension applies the same survival logic to the Neptunian desert–ridge–savanna landscape. Mapping HEM survival constraints into the period–radius plane with empirically inferred mass–radius relations reproduces the slope of the Neptunian desert boundary across $a_{\rm F}\ge 2\,r_{\rm tide}.$3 with a single representative value $a_{\rm F}\ge 2\,r_{\rm tide}.$4, while the observed density dispersion broadens the disruption line into a finite survival band that traces the Neptunian ridge at $a_{\rm F}\ge 2\,r_{\rm tide}.$5–$a_{\rm F}\ge 2\,r_{\rm tide}.$6 d (Castro-González et al., 17 Apr 2026). In the period–density plane, the same framework yields a density-dependent survival threshold $a_{\rm F}\ge 2\,r_{\rm tide}.$7 and is reported to be consistent with a persistent ridge concentration near $a_{\rm F}\ge 2\,r_{\rm tide}.$8 (Castro-González et al., 17 Apr 2026).

Dynamical-tide models sharpen this demographic picture by distinguishing three regimes as a function of pericenter distance: stagnant, diffusive, and outflow (Zanazzi et al., 18 Jun 2026). In the diffusive regime, $a_{\rm F}\ge 2\,r_{\rm tide}.$9-modes shock but cool radiatively, allowing the envelope to survive while the orbit circularizes and bunches near the hot-Jupiter pile-up and Neptune ridge. In the outflow regime, the same shocks drive super-Eddington winds that unbind the envelope, leaving dense cores in the Neptune desert (Zanazzi et al., 18 Jun 2026). This suggests that the hot-Jupiter pile-up, Neptune ridge, and Neptune desert may be connected manifestations of a single HEM-plus-dynamical-tide framework rather than independent features.

6. Limits, extensions, and unresolved issues

HEM is neither a universal explanation for all close-in planets nor a mechanism restricted to ordinary hot-Jupiter systems. Around polluted white dwarfs, small bodies near the tidal disruption radius can be modeled as HEM products from initial distances of several astronomical units, with either tidal friction or ram-pressure drag in a compact disc providing the dissipative circularization (O'Connor et al., 2020). For tidal migration to work in that context, a viscosity comparable to molten rock is required, while disc migration implies a minimum disc column density and a total disc mass consistent with the metal budgets inferred for polluted white dwarfs (O'Connor et al., 2020).

HEM also has strong consequences for satellites. In binary-driven HEM of hot Jupiters, an exomoon acts as an additional short-range force that can quench ZLK excitation; if migration nevertheless proceeds, the moon typically spirals onto the planet, escapes, or is ejected (Trani et al., 2020). In the simulations reported there, efficient shielding occurs in only 0.6% of low-mass-moon cases and 9.1% of high-mass-moon cases, and there were no cases in which a hot Jupiter migrated inward while retaining its moon (Trani et al., 2020). This suggests that a confirmed long-lived exomoon around a hot Jupiter would disfavor at least the binary-driven HEM channel considered in that work.

Several controversies remain open. The Kepler super-eccentric deficit constrains standard steady-state equilibrium-tide HEM, but chaotic or diffusive tides, inflated eccentric migration, or HEM beginning only after some prior inward transport can shorten residence times in the super-eccentric phase and thereby weaken that tension (Jackson et al., 2022, Yu et al., 2021). Low obliquity does not uniquely imply disk migration, because CHEM and spin-orbit realignment can mimic aligned outcomes; conversely, high obliquity does not uniquely prove HEM, because primordial disk misalignment can generate it before migration (Kawai et al., 19 Sep 2025, Vick et al., 2022). The present literature therefore points toward a pluralistic interpretation in which HEM is important, sometimes dominant for specific subpopulations, but not exhaustive.

A plausible synthesis is that HEM is most secure where three ingredients coincide: a demonstrable route to extreme eccentricity, a tidal-dissipation mechanism fast enough to circularize before disruption, and a population-level imprint consistent with the associated survival boundary. Where one of these ingredients is absent—such as systems with extreme eccentricity but tidal timescales longer than the system age, or circular close-in giants with $e \gtrsim 0.999$0—other pathways, especially disk migration, remain necessary (Angelo et al., 2022, Kawai et al., 19 Sep 2025).