Where Do Hot Jupiters Come From? Revisiting Tidal Disruption and Ejection in High-Eccentricity Migration

TL;DR

This study uses 3D hydrodynamic simulations to show that dense cores in giant planets prevent total disruption during high-eccentricity migration, altering previous models of hot Jupiter formation.

Contribution

It demonstrates that core-bearing giant planets are more resilient to tidal disruption, revises the tidal exclusion zone, and clarifies energy changes in highly eccentric orbits.

Findings

No total disruptions occur within 2.7 tidal radii for core-bearing planets.

Planets can survive deep encounters by downsizing or ejection, contributing to free-floating planets.

Mass loss is minimal at wider encounters, allowing gradual circularization into hot Jupiters.

Abstract

The origin of hot Jupiters remains a key open question. In the high-eccentricity migration scenario, traditional coreless models predict a strict tidal exclusion zone within $\sim 2.7$ tidal radii $r_\textrm{t}$, in which giant planets are either fully disrupted or ejected. We revisit this limit using three-dimensional hydrodynamic simulations of giant planets with realistic dense cores (10 - 20 $M_\oplus$). We find that even a few-percent-mass core fundamentally changes the outcome: \textbf{no total disruptions} occur within the previously suggested destruction zone ($\lesssim 2.7 \, r_\textrm{t}$). For deep encounters ($\lesssim 1.7 \, r_\textrm{t}$) planets suffer severe envelope stripping and are either progressively downsized to dense remnants or ejected after a few close encounters, possibly contributing to the free-floating planet population. In the intermediate regime ($ \sim 1.7 $--$2.0, r_\mathrm{t}$), planets experience significant partial mass loss over repeated encounters. For wider encounters ($ \gtrsim 2.0\, r_\mathrm{t} $), mass loss is minimal, allowing the planets gradually circularize into hot Jupiters. Furthermore, we show that for highly eccentric orbits ($e\gtrsim 0.9$), the change in specific orbital energy $ \Delta E_{\mathrm{orb}} $ depends primarily on periastron distance $ r_\mathrm{p} $ rather than semi-major axis $ a $. This enables us to extrapolate our fixed-$ a $ results across a broad ($a$, $e$) parameter space and identify a well-defined tidal ejection zone whose sharp boundaries converge asymptotically. Our results highlight the crucial role of planetary internal structure in high-eccentricity migration and suggest that the survival and transformation of core-bearing giant planets are far more common than previously thought.