- The paper develops a full-harmonic Maxwell viscoelastic model coupled to core thermal evolution, showing that tidal heating and cooling reach a quasi-steady balance at temperatures of several thousand kelvin.
- The paper finds that 3-, 4-, and 5-day final-period tracks circularize on similarly few-billion-year timescales, contrasting with the roughly 15-fold spread predicted by fixed-time-lag tides.
- The paper predicts a weak damping dependence of approximately tₑ ∝ Pᶠ⁰·³–⁰·⁵, helping explain eccentric hot Neptunes and linking their eccentricities to observable heating, radius inflation, and atmospheric escape.
Motivation and observational context
Close-in Neptune-like planets present a dynamical puzzle that standard tidal theory handles poorly. A population of warm Neptunes retains eccentricities of order e∼0.1–$0.2$ despite strong tidal forcing, and several close-in Neptunes and sub-Saturns show significant spin–orbit misalignments, both suggestive of high-eccentricity migration (HEM). In a preliminary archival analysis (NASA Exoplanet Archive, queried 2026 June 3), the authors find that the eccentric fraction of isolated hot Neptunes (10<Mp/M⊕<30) depends only weakly on the inferred circularization radius af=a(1−e2): in the innermost bin (af<0.04 AU), 5 of 10 hot Neptunes have e>0.1, versus only 4 of 48 hot Jupiters. Two-sample KS tests reinforce the contrast: the af distributions of eccentric and non-eccentric hot Neptunes are statistically indistinguishable (D=0.18, p=0.92), while for hot Jupiters eccentric systems are strongly shifted to larger af ($0.2$0, $0.2$1).
This behavior is difficult to reconcile with fixed-efficiency equilibrium tides. Constant-time-lag (CTL) and constant-$0.2$2 models predict circularization times scaling steeply as $0.2$3 along constant-angular-momentum tracks, which would erase eccentricity too efficiently at small separations and fail to damp it farther out. The key physical distinction the authors exploit is that Neptune-mass planets carry most of their mass in dense rocky or icy cores, where dissipation is viscoelastic and strongly temperature-dependent — unlike gas giants, whose extended envelopes are usually assumed to dominate.
The model treats the core as a self-gravitating Maxwell body of radius $0.2$4 and mass $0.2$5, computing the quadrupolar tidal heating as a full harmonic sum over azimuthal orders $0.2$6 and all integer harmonics $0.2$7:
$0.2$8
with $0.2$9 built from Hansen coefficients and forcing frequencies 10<Mp/M⊕<300 in the planetary rotating frame. The dissipative Love number follows the Maxwell form 10<Mp/M⊕<301, peaking near 10<Mp/M⊕<302. The relaxation time combines a material term 10<Mp/M⊕<303 and a gravity-controlled viscous term 10<Mp/M⊕<304; the Love-number contrast is set by the ratio of shear modulus to gravitational rigidity 10<Mp/M⊕<305.
Two technical points deserve emphasis. First, at high eccentricity the pseudo-synchronous spin is rapid (at 10<Mp/M⊕<306, 10<Mp/M⊕<307), so corotation for the dominant 10<Mp/M⊕<308 terms lies near 10<Mp/M⊕<309: heating is produced by many high-order harmonics overlapping the Maxwell response, not by a single frequency comparable to the mean motion. The authors define a dissipation-weighted effective forcing frequency af=a(1−e2)0 and show it can exceed af=a(1−e2)1 by more than an order of magnitude; replacing the spectrum by af=a(1−e2)2 can place a planet on the wrong side of the Maxwell peak. Second, along a constant-angular-momentum track, af=a(1−e2)3 decreases during circularization even though af=a(1−e2)4 increases — the relevant forcing frequency is set by the evolving harmonic spectrum and spin state, not by af=a(1−e2)5 alone.
The spin is assumed to sit on the highest-frequency stable zero-torque branch, computed numerically from the full torque sum. In the solid-like regime this branch exhibits discrete spin–orbit resonances at af=a(1−e2)6; at high eccentricity these become densely spaced, and the authors provide an envelope fit accurate to within roughly 10% for af=a(1−e2)7 that recovers Hut's pseudo-synchronous result exactly in the weak-friction limit.
Thermally regulated migration
The central calculation couples orbital evolution along af=a(1−e2)8 to a one-zone core temperature equation af=a(1−e2)9, with Arrhenius viscosity (af<0.040), melt-fraction-weakened shear modulus, and a power-law cooling luminosity af<0.041. Fiducial parameters correspond to a af<0.042 planet with a af<0.043, af<0.044 dissipating core.
The coupled evolution is self-regulated. Starting from cold cores (af<0.045 K), tidal heating rapidly lowers the viscosity and shifts the response from the elastic side toward the fluid-like side of the Maxwell peak. After a short transient, each track settles into a quasi-steady state in which af<0.046, with equilibrium temperatures of a few thousand kelvin. Notably, only the shortest-period track (af<0.047 d) briefly exceeds the adopted solidus temperature: the quasi-steady state does not require partial melting, which acts mainly as high-temperature regularization.
The headline numerical result: tracks with final circular periods af<0.048, 4, and 5 d reach low eccentricity on nearly indistinguishable few-Gyr timescales, despite final periods differing by nearly a factor of two. A fixed-af<0.049 model would predict e>0.10, a factor of e>0.11 spread between 3 and 5 d. The compensation mechanism is explicit: stronger forcing at smaller e>0.12 heats the core more, pushes the response farther into the fluid-like regime, and raises e>0.13, slowing subsequent damping.
Quasi-steady scaling
In the fluid-like limit (e>0.14), the weak-friction response reduces to a lag proportional to e>0.15, and thermal balance yields an implicit equation for the equilibrium temperature. Differentiating the balance condition gives a local power-law index
e>0.16
where e>0.17 measures the Arrhenius sensitivity. Because rocky/icy viscosity is strongly temperature-dependent (e>0.18) while equilibrium temperatures satisfy e>0.19–af0, the exponent collapses from the CTL value of 8 to af1, i.e., af2, robust across cooling indices af3–af4. This analytic flattening explains the numerical results and constitutes the paper's strongest quantitative claim.
Testable predictions and observational consistency
The model makes three predictions: (i) eccentric hot Neptunes should persist to smaller circularization radii than fixed-efficiency models allow; (ii) eccentricity should correlate with tidal-heating signatures such as radius inflation and enhanced atmospheric escape — GJ 3470 b, a polar warm Neptune with af5 and an evaporating atmosphere, is cited as suggestive; and (iii) actively migrating systems should produce an extended population of Neptune-mass planets with af6–af7 and small pericenters, for which Kepler-1656 b (af8) may be an example. The archival eccentric-fraction comparison supports prediction (i) qualitatively, though the authors are careful to label it a consistency check rather than a definitive test given the small sample.
Limitations and open questions
Several assumptions bound the results. The cooling prescription is deliberately one-zone and does not resolve heat transport through the core-envelope boundary; the gaseous envelope may regulate the cooling luminosity, introduce additional thermal timescales, and couple back through radius inflation or mass loss. The melt treatment is a phenomenological proxy rather than a proper interior model, and the spin-equilibrium fit applies only over af9 during the elastic-to-fluid transition. The mechanism's applicability to envelope-free rocky planets and super-Earths remains uncertain, since faster cooling could prevent long-lived eccentric phases; determining where the transition between rapidly cooling rocky planets and thermally buffered Neptune-like planets lies requires coupled interior-orbital models beyond this work. Finally, whether the observed weak D=0.180 dependence survives as the hot-Neptune sample grows is an open empirical question.
Conclusion
This work demonstrates that when tidal dissipation in eccentric migrating planets is controlled by a temperature-dependent viscoelastic core, the migration self-regulates: heating drives the core toward peak dissipation, then further heating suppresses it, producing a quasi-steady heating–cooling balance and long-lived eccentric phases. The resulting eccentricity damping time depends only weakly on final orbital distance — approximately D=0.181 versus D=0.182 for fixed-efficiency tides — offering a natural explanation for eccentric hot Neptunes across a broad range of separations. More broadly, the framework implies that present-day eccentricities, periods, radii, and ages can constrain interior properties such as viscosity, melting state, and cooling efficiency that are otherwise inaccessible (2606.14983).