Formation and tidal evolution of hot super-Earths in multiple planetary systems

Hot super-Earths are exoplanets with masses < 10 Earth masses and orbital periods < 20 days. Around 8 hot super-Earths have been discovered in the neighborhood of solar system. In this lecture, we review the mechanisms for the formation of hot super-Earths, dynamical effects that play important roles in sculpting the architecture of the multiple planetary systems. Two example systems (HD 40307 and GJ 436) are presented to show the formation and evolution of hot super-Earths or Neptunes.

💡 Research Summary

The paper provides a comprehensive review of the formation pathways and subsequent tidal evolution of hot super‑Earths—planets with masses below 10 M⊕ and orbital periods shorter than 20 days—focusing on how these processes operate within multi‑planet systems. It begins by summarizing the observational landscape: eight nearby hot super‑Earths have been detected, displaying a range of masses (1–9 M⊕), radii, eccentricities, and host‑star metallicities. These empirical constraints already rule out simple, single‑planet formation scenarios and demand models that can simultaneously explain compact orbital architectures, resonant chains, and sometimes surprisingly high eccentricities.

Two broad formation paradigms are examined. The first, in‑situ accretion, posits that solid material (planetesimals, pebbles, and dust) is sufficiently concentrated inside ~0.1 AU to allow rapid core growth to super‑Earth size before the gas disc dissipates. The authors discuss the required surface density (Σ ≳ 10³–10⁴ g cm⁻²), the role of pebble accretion efficiencies, and the temperature‑dependent condensation of silicates and metals, which together set an upper mass limit near 10 M⊕. The second paradigm involves formation at larger radii followed by Type I migration. Analytic expressions for the migration timescale (τ_I ∝ (M⋆/Mₚ)(M⋆/Σa²)(h/a)²Ω⁻¹) are presented, and the paper emphasizes that migration is typically faster than disc dispersal, causing multiple embryos to converge toward the inner disc edge. Convergent migration naturally produces mean‑motion resonances (2:1, 3:2, etc.), creating resonant chains that are observed in many compact systems. When the gas density drops, migration stalls, and the remaining solid material continues to collide, merging embryos into a handful of super‑Earths with final masses of 5–10 M⊕.

After formation, long‑term tidal interactions dominate orbital evolution. The authors adopt the constant‑time‑lag model, describing tidal dissipation through the planetary quality factor Q and Love number k₂. For rocky planets, typical values are Q≈10–100 and k₂≈0.3, leading to tidal decay timescales τ_tide ≈ (a⁸/G M⋆³ Rₚ⁵)(Q/k₂). This formulation predicts that a 5 M⊕ planet at 0.05 AU will circularize and spin‑synchronize within a few hundred Myr to a few Gyr, depending on Q. The paper also discusses how tidal heating can affect interior structure, potentially inflating radii or altering Q over time.

Two case studies illustrate the theory. HD 40307 hosts three super‑Earths with periods of 4.3, 9.6, and 20.4 days and minimum masses of 4–8 M⊕. The authors run N‑body simulations coupled to a viscously evolving disc model. Starting with a swarm of 0.1–0.5 M⊕ embryos inside 0.2 AU, the embryos undergo convergent Type I migration, lock into a resonant chain, and later experience tidal damping that breaks the resonances, reproducing the observed non‑resonant configuration. The simulations also show that modest variations in disc viscosity (α≈10⁻³–10⁻⁴) can shift the final spacing, matching the observed period ratios.

GJ 436b, a Neptune‑mass (≈22 M⊕) planet on a 2.64‑day orbit with a surprisingly high eccentricity (e≈0.15), challenges a simple tidal circularization picture. Using the same tidal framework, the authors find that achieving the present eccentricity after several Gyr would require an unusually high planetary Q (≈10⁵) or a continuous source of eccentricity pumping. They explore the latter by inserting a low‑mass, non‑transiting companion in an outer orbit. Secular interactions, especially near a Laplace‑Lagrange resonance, can maintain e≈0.15 for GJ 436b over Gyr timescales, provided the companion’s mass exceeds ~2 M⊕ and its semi‑major axis lies between 0.05 and 0.1 AU. This scenario also predicts a small but detectable radial‑velocity signal, offering a testable prediction.

In the concluding section, the authors synthesize the findings: hot super‑Earths likely arise from a combination of high‑density inner discs, rapid Type I migration, resonant capture, and subsequent tidal damping. The relative importance of each process depends on disc properties (viscosity, temperature gradient, surface density), planetary interior parameters (Q, k₂), and the presence of additional planets that can sustain eccentricities. The paper highlights that future high‑precision transit timing variations (TTVs), radial‑velocity monitoring, and asteroseismic constraints on host‑star ages will be crucial for refining Q estimates and testing the proposed multi‑planet tidal evolution models. Ultimately, the integration of disc‑driven dynamics with tidal physics offers a coherent framework for interpreting the diverse architectures of observed hot super‑Earth systems.