How do Neptunes survive tidal migration without losing more energy than they can bind?

The energy budget that threatens migrating Neptunes

High eccentricity migration (HEM) is the idea that a planet starts on a distant, highly eccentric orbit and shrinks its orbit by dissipating tides raised at each close passage. The paper states the core difficulty plainly: shrinking an orbit by roughly one to two orders of magnitude requires a planet to dissipate more than about 10 times the gravitational binding energy held in its own structure. Because Neptunes have lower binding energies than gas giants, this budget problem is sharper for them, and the authors call it the Neptunian "final parsec" problem.

The framing matters because it separates two questions that are easy to conflate. One is whether tides can remove orbital energy fast enough; the other is whether the planet can absorb and dispose of that energy without being destroyed. The paper's target is the second half of the migration, from a semi-major axis of about 0.5 au down to about 0.05 au, after earlier work had shown that f-mode excitation at eccentricities near 0.99 could cover the first order of magnitude, from roughly 5 to 0.5 au.

Inertial waves carry the tidal response

The authors perform coupled orbit and structure calculations that track the planet's frequency-dependent tidal response under inertial waves. Inertial waves are oscillations restored by the Coriolis force in a rotating, convectively unstable layer; in this model they are hosted in the planet's convective envelope, and their dissipation is what drains orbital energy. The calculation solves for the wavelike part of the tidal response, computes a viscous dissipation rate, and converts it into the modified tidal quality factor Q-prime, where a smaller Q-prime means more energy lost per tidal period relative to the maximum energy stored in the tide.

Two structural variables control how strongly the planet responds. A larger planetary radius, or a smaller core relative to the envelope, lowers Q-prime because the tidal forcing overlaps more with the planet's global eigenmodes and because the dynamical frequency is shorter. The paper's frequency-averaged values make the size of this effect concrete: for a 30 Earth-mass planet with a 20 Earth-mass core, the averaged Q-prime falls from about 16,295 at entropy 6.0 to about 1,506 at entropy 9.0, while halving the core to 10 Earth-mass gives roughly 10 times lower averaged Q-prime at a given entropy.

Heating inflates the planet and feeds itself

The central mechanism is a feedback loop rather than a single dissipation event. Tidal heating puffs the planet up; a larger, more weakly bound planet is deformed more strongly by the tide; the stronger non-wave displacement seeds a larger inertial wave response; and that response dissipates more tidal energy. In the paper's fiducial case, a 30 Earth-mass Neptune with a 20 Earth-mass core starts at a semi-major axis of 0.3 au and eccentricity 0.896, with an initial entropy near 7 k_B/m_H and a cooling time of about 10 Gyr. Tidal heating then dwarfs radiative cooling for most of the evolution, producing runaway inflation on a timescale of about 10^7 yr, fast compared with typical system ages.

The authors explicitly test and reject a plausible hope: that the bursty, resonant character of inertial wave dissipation might let planets avoid runaway. Instead, the resonant frequencies in Q-prime are always excited strongly enough to drive vigorous dissipation. The spin period also jumps between local torque equilibria as the planet inflates, so the tidal luminosity dips at each readjustment and then grows again once the spin drifts away from the instantaneous minimum-torque rate.

Inflation brings f-modes back into resonance

The reason runaway inflation does not simply destroy the planet is that it re-engages a second dissipation channel. Fundamental modes, or f-modes, are global oscillations of the planet on a timescale comparable to its dynamical time; earlier work found they can drive orbital shrinkage at very high eccentricity but become inefficient once eccentricity drops below about 0.9. As the planet inflates at moderate eccentricity, its f-mode frequencies sweep back into resonance with the tidal forcing, and the orbit integral measuring the coupling between the tidal potential and the f-mode grows steeply with entropy.

The energy accounting is the key result. The paper estimates that reaching the entropy of about 8.75 k_B/m_H needed for strong coupling, starting from about 7 k_B/m_H, costs roughly 0.1 of the planet's binding energy, scaling with the entropy change, the integrated temperature profile, and the planet's radius and mass. That is a lower limit, since some energy leaves through radiative cooling, but it is still a small fraction of the binding energy — in contrast to the more than tenfold orbital energy the planet must shed. The authors therefore describe f-mode excitation as a route to safe circularization, conditional on the f-mode dissipating without harming the planet.

Gas giants absorb the same heat without drama

The paper's comparison case is an eccentric gas giant at eccentricity near 0.9. It circularizes over about 1 Gyr with comparatively little structural upheaval, because its radius changes far less with entropy than a Neptune's does; the flat radius-entropy relation lets radiative cooling luminosity grow enough to balance the tidal luminosity. The authors also check that periastron distances stay above the tidal disruption limit throughout, so the giant is not stripped.

This matters for how the result should be read. The gas giant still dissipates more orbital energy than its structure contains, and it survives by radiating the excess to space. The paper notes that at the eccentricity shown, the orbit integral grows from about 3 to 30 as the planet puffs up, which suggests the f-mode may dominate circularization rather than inertial waves; even so, the calculation shows f-mode dissipation is not the only route to survival. The authors also state that depositing tidal heat deep inside a gas giant, which maximizes radius inflation, still does not produce runaway inflation under inertial waves.

The mechanism only works inside about 0.1 au

The survival route is not universal, and the paper is explicit about the boundary. Inertial waves can only help circularize Neptunes that end up inside roughly 0.1 au; larger orbital angular momentum throttles the heating, because a longer periastron passage time spins the planet down, and slowly rotating planets sample fewer tidal harmonics and show a larger Q-prime. The authors describe the outcome as all-or-nothing: in their example, planets starting with spin periods slower than about 14 days, or with sufficiently large orbital angular momentum, avoid violent heating, while slightly faster spins or less angular momentum undergo runaway inflation.

Viscosity is a major uncertainty, since the damping mechanism for inertial waves is poorly understood. The authors therefore test a viscosity 10^5 times smaller than their fiducial choice and find the same resonant harmonics drive the evolution, with runaway taking about 10^7 yr longer because the resonant troughs become narrower. They conclude the main results are not sensitive to the adopted viscosity. The calculations also assume a large convective zone beneath a thin radiative layer; the authors argue the mechanism only needs some convective region to host inertial waves, and note that interior models consistent with gravity and magnetic field measurements do suggest sizable convection zones.

A testable dearth of eccentric warm Neptunes

The paper turns the timescale difference into an observational prediction. Because runaway inflation under inertial waves happens in about 10^7 yr for Neptunes but circularization takes about 1 Gyr for gas giants, warm Neptunes on a high eccentricity migration track — semi-major axis above about 0.1 au and eccentricity above about 0.8 — should be rare, roughly 100 times rarer than their gas giant counterparts. Moderately eccentric warm Neptunes, by contrast, should be common, because slow equilibrium spin throttles inertial wave dissipation for those planets.

The prediction comes with a stated exception that sharpens it: a highly eccentric warm Neptune could persist indefinitely if its rotation period exceeds about two weeks. Such an outlier should show negligible rotational oblateness, with an oblateness parameter below about 10^-4, and a low-entropy interior. The authors also suggest that if f-modes do not fully circularize the orbit after runaway inflation, hot Neptunes around K and M stars may keep small but non-zero eccentricities near 0.05 indefinitely, since inertial wave dissipation drops sharply as planets spin down toward synchronism. The most useful reading of this paper is therefore conditional rather than blanket: the survival mechanism is tied to low orbital angular momentum and final circularized distances inside about 0.1 au, and the cleanest observational test is the predicted scarcity of eccentric warm Neptunes on that track, paired with the slow-rotator exception.

Sources used in this answer