Tidal Resonances in Astrophysical Systems

• Tidal resonances are frequency‐matching phenomena where tidal forcing aligns with natural oscillation modes in various astrophysical systems.
• They enable precise tidal dissipation, resonance locking, and orbital evolution by channeling energy through discrete spectral interactions.
• Analytical models and simulations reveal that resonance signatures, such as phase shifts and dissipation spikes, serve as key diagnostics for gravitational wave and interior structure studies.

Tidal resonances are resonant responses produced when a tidal forcing frequency, or an integer combination of orbital frequencies associated with a tidal perturbation, matches a natural mode, a geometric commensurability, or a resonance-defined manifold of the system. In the literature this term covers several distinct but related phenomena: excitation of inertial modes in rotating convective envelopes, standing g-mode resonances in stellar radiative cores, gravito-inertial and crust-coupled mode excitations in neutron stars, resonance locking in eccentric stellar binaries, Laplace-type resonant states in compact multiplanet chains, transient tidal resonances in extreme-mass-ratio inspirals, resonant global ocean and core–basal magma ocean tides, and resonant tidal response functions in horizonless compact objects (Lin et al., 2021, Brasser et al., 2022, Gupta et al., 2021, Wei, 2020).

1. General framework and resonance conditions

Across these settings, the common structure is a forcing with discrete spectral content acting on a system that supports modes, characteristics, or commensurabilities. In rotating fluid bodies, the relevant band is often the inertial interval $|\omega|\le 2\Omega$, with inertial-wave dispersion

$$\omega=\pm 2\Omega \frac{k_z}{|k|}=\pm 2\Omega \cos\alpha,$$

and critical latitude

$\cos\theta_c=\frac{|\omega|}{2\Omega}.$

In compact planetary chains, the relevant conditions are two-body and three-body commensurabilities such as $(j+1)n'-jn\approx 0$ and

$k_1 n_1+k_2 n_2+k_3 n_3\approx 0.$

In EMRIs with an external tidal perturber, resonance occurs when

$m\omega_\phi+k\omega_\theta+n\omega_r=0.$

In binary neutron stars, the dominant quadrupolar resonances satisfy $m\Omega(t_{\rm res})=\omega_\alpha$, and for $l=2$, $m=\pm2$ one has $f_\alpha=f_{\rm GW}(t_{\rm res})$; in an early-Earth core–BMO system, resonance occurs when the semidiurnal forcing frequency is close to the natural frequency $$\omega=\pm 2\Omega \frac{k_z}{|k|}=\pm 2\Omega \cos\alpha,$$0 of the interfacial mode (Lin et al., 2021, Papaloizou, 2014, Gupta et al., 2021, Gittins et al., 4 Jun 2026, Kiernan et al., 15 Jun 2026).

Domain	Resonant object	Representative condition
Rotating convective envelope	Inertial mode	$$\omega=\pm 2\Omega \frac{k_z}{|k|}=\pm 2\Omega \cos\alpha,$$1
Three-planet chain	Laplace-type commensurability	$$\omega=\pm 2\Omega \frac{k_z}{|k|}=\pm 2\Omega \cos\alpha,$$2
EMRI with tidal perturber	Kerr harmonic	$$\omega=\pm 2\Omega \frac{k_z}{|k|}=\pm 2\Omega \cos\alpha,$$3
Binary neutron star	Stellar QNM	$$\omega=\pm 2\Omega \frac{k_z}{|k|}=\pm 2\Omega \cos\alpha,$$4
Core–BMO interface	Interfacial inertia–self-gravity mode	$$\omega=\pm 2\Omega \frac{k_z}{|k|}=\pm 2\Omega \cos\alpha,$$5

The resonant observable is likewise context-dependent. It may be a spike in viscous dissipation, a jump in orbital constants of motion, a departure from exact commensurability, a kink in the gravitational-wave phase, a crust-breaking strain, or a pole in a tidal response function. This suggests that “tidal resonance” is best understood as a family of frequency-matching phenomena rather than a single mechanism.

2. Wave and mode resonances in stars and giant planets

In rotating stars and planets with convective envelopes, excitation of inertial waves provides an important channel for tidal dissipation, but the dissipation depends erratically on the tidal frequency. The central result of Lin and Ogilvie is that resonant tidal peaks are produced by global inertial eigenmodes that possess a substantial large-scale, smooth flow component hidden beneath small-scale, localized wave beams. The tidal forcing couples strongly to this large-scale part, while viscous dissipation is concentrated in the beams; at representative peaks they measured $$\omega=\pm 2\Omega \frac{k_z}{|k|}=\pm 2\Omega \cos\alpha,$$6, $$\omega=\pm 2\Omega \frac{k_z}{|k|}=\pm 2\Omega \cos\alpha,$$7, and $$\omega=\pm 2\Omega \frac{k_z}{|k|}=\pm 2\Omega \cos\alpha,$$8. They also emphasized that such peaks are often mistakenly attributed to wave attractors, whereas attractor-dominated responses are broader and less eigenfrequency-locked (Lin et al., 2021).

In solar-type stellar cores, tidally excited internal gravity waves define a different resonant regime. For low-amplitude forcing, the response exhibits resonances with standing g-modes at particular frequencies. For high-amplitude forcing, the excited waves break promptly near the centre and spin up the core so that subsequent waves are absorbed in an expanding critical layer. For intermediate amplitude, linear damping gradually spins up the core, the system can evolve toward or away from resonance depending on the initial detuning, and eventually a critical layer forms. The torque then transitions toward the traveling-wave value, with $$\omega=\pm 2\Omega \frac{k_z}{|k|}=\pm 2\Omega \cos\alpha,$$9, and the classical picture of resonance locking in solar-type stars is revised because spinning up only $\cos\theta_c=\frac{|\omega|}{2\Omega}.$0 of $\cos\theta_c=\frac{|\omega|}{2\Omega}.$1 shifts the wave by one full resonance (Guo et al., 2022).

In hot Jupiters, strong stellar irradiation allows an interior radiative zone to form below the outer convective layer, and these radiative zones support g-modes that can be resonantly forced by the star’s gravitational tide. The mode response peaks when the Lorentzian denominator is minimized,

$\cos\theta_c=\frac{|\omega|}{2\Omega}.$2

and the paper argues that tidal dissipation in g-modes modifies the thermal structure in a way that tunes the g-mode eigenfrequencies toward resonance with the forcing. In that model, thermomechanical feedback, deep heat trapping, and the resulting entropy increase in the central convective region can drive radius inflation on Myr timescales (Jermyn et al., 2017).

In eccentric stellar binaries, resonance locking between axisymmetric g-modes and tidal harmonics provides sustained dissipation when

$\cos\theta_c=\frac{|\omega|}{2\Omega}.$3

For $\cos\theta_c=\frac{|\omega|}{2\Omega}.$4–$\cos\theta_c=\frac{|\omega|}{2\Omega}.$5 primaries, the work on resonance locking finds that orbital evolution via resonance locking occurs primarily during the pre-main-sequence phase, with the main-sequence phase contributing negligibly. Ignoring nonlinearity, resonance locking can circularize binaries out to circular periods of $\cos\theta_c=\frac{|\omega|}{2\Omega}.$6–$\cos\theta_c=\frac{|\omega|}{2\Omega}.$7 days; a saturated resonance lock reduces the circularization period by about a third, but resonance locking remains much more effective than the cumulative actions of equilibrium tides (Zanazzi et al., 2021).

3. Resonant chains and dissipative evolution in planetary systems

In compact multiplanet systems, tidal resonances are dynamical states in which resonant forcing from mean-motion resonances and Laplace-type three-body resonances continually excites small but nonzero eccentricities, while tidal dissipation within the planets removes orbital energy. The coupled result is system-wide, resonance-guided divergent migration at nearly conserved angular momentum: as eccentricities decrease, adjacent period ratios move slightly farther from exact commensurability, a process described as resonant repulsion. In TRAPPIST-1, this framework explains why the planets maintain small eccentricities and remain very near resonant period ratios while evolving along equilibrium curves in semimajor axis–eccentricity space; the inferred long-term tidal parameters for planets b–e are in the range $\cos\theta_c=\frac{|\omega|}{2\Omega}.$8 to $\cos\theta_c=\frac{|\omega|}{2\Omega}.$9 in $(j+1)n'-jn\approx 0$0 (Brasser et al., 2022).

For three-planet systems near first-order resonances, an analytic time-averaged treatment shows that once small-amplitude libration of the relevant resonant angles has been attained, tide-driven eccentricity damping removes orbital energy while leaving total angular momentum nearly constant, forcing the chain to spread. The departure from strict commensurability then follows the characteristic one-third law,

$(j+1)n'-jn\approx 0$1

so the offset from exact resonance grows as $(j+1)n'-jn\approx 0$2. Applied to Kepler-60, this provides bounds on the time average of $(j+1)n'-jn\approx 0$3 depending on whether a librating state has or has not been attained (Papaloizou, 2014).

Kepler-221 illustrates a related but distinct configuration: a pure three-body resonance in which

$(j+1)n'-jn\approx 0$4

librates while no two-body angles are observed to librate in the dynamical model. The measured frequency combination

$(j+1)n'-jn\approx 0$5

is exceptionally close to zero, and direct capture into such a pure three-body resonance is found to be practically impossible for super-Earth/sub-Neptune masses at those separations. The proposed origin is an earlier resonant chain followed by large-scale divergence under tidal dissipation, with obliquity tides identified as an excellent candidate for driving the divergence (Goldberg et al., 2021).

A broader survey of STIPs with three-body resonant chains reaches the same structural conclusion. For Kepler-80, Kepler-223, K2-138, TOI-178, and TRAPPIST-1, the trend of resonance offsets between adjacent pairs is reproduced by analytical estimates derived from the three-planet resonant dynamics, and N-body simulations show that tidal damping effects preserve the resonant patterns while moving the systems along the three-body resonances. The study finds that the trends in the offsets can be produced regardless of the considered value for the tidal factor, while the amplitudes for Kepler-80, K2-138, and TOI-178 can also be reproduced with appropriate $(j+1)n'-jn\approx 0$6 and the estimated system ages (Charalambous et al., 2023).

4. Relativistic inspirals and neutron-star seismology

In EMRIs perturbed by a nearby star or black hole, tidal resonances arise when a Fourier harmonic of the external tidal field becomes stationary with respect to the Kerr geodesic motion. The resonance condition

$(j+1)n'-jn\approx 0$7

produces a coherent jump in the orbit’s constants of motion over a resonance duration that is long compared to an orbital period but short compared to the radiation-reaction time. The instantaneous changes are small, but the lasting offsets in orbital frequencies accumulate into gravitational-wave phase shifts that can be significant over the many cycles in band; low-order resonances are common over the surveyed parameter space, and dephasings above $(j+1)n'-jn\approx 0$8 rad are detectable for a significant fraction of sources (Gupta et al., 2021).

For waveform modeling, this effect can be incorporated as a discrete jump at the resonance surface. The Resonance Model developed for transient EMRI resonances uses a step-function update to the actions and benchmarks it against osculating-orbit evolutions. In the reported comparisons, the mismatch between resonant waveforms built with this prescription and full osculating resonant waveforms remains at or below $(j+1)n'-jn\approx 0$9 across the inspiral, while neglecting resonances produces significantly larger mismatches and systematic biases that can exceed statistical errors at SNR $k_1 n_1+k_2 n_2+k_3 n_3\approx 0.$0 (Gupta et al., 2022).

In coalescing neutron-star binaries, the dominant quadrupolar tide can resonantly excite low-frequency stellar modes in the last seconds before merger. In a GR treatment with realistic EOSs consistent with GW170817 and modest stratification $k_1 n_1+k_2 n_2+k_3 n_3\approx 0.$1, $k_1 n_1+k_2 n_2+k_3 n_3\approx 0.$2 frequencies lie near $k_1 n_1+k_2 n_2+k_3 n_3\approx 0.$3–$k_1 n_1+k_2 n_2+k_3 n_3\approx 0.$4 Hz, so resonance occurs at orbital frequencies $k_1 n_1+k_2 n_2+k_3 n_3\approx 0.$5–$k_1 n_1+k_2 n_2+k_3 n_3\approx 0.$6 Hz. The amplitude scaling found numerically is

$k_1 n_1+k_2 n_2+k_3 n_3\approx 0.$7

and $k_1 n_1+k_2 n_2+k_3 n_3\approx 0.$8 modes with $k_1 n_1+k_2 n_2+k_3 n_3\approx 0.$9 can exceed a crust-breaking threshold $m\omega_\phi+k\omega_\theta+n\omega_r=0.$0 (Kuan et al., 2021).

The companion study on precursor flares shows that, for some combinations of stellar rotation, stratification, and magnetic field structure, resonantly excited $m\omega_\phi+k\omega_\theta+n\omega_r=0.$1- and $m\omega_\phi+k\omega_\theta+n\omega_r=0.$2-modes can trigger crustal failure and release $m\omega_\phi+k\omega_\theta+n\omega_r=0.$3–$m\omega_\phi+k\omega_\theta+n\omega_r=0.$4 erg, consistent with short-GRB precursor energetics. It also emphasizes that magnetic and rotational shifts can move $m\omega_\phi+k\omega_\theta+n\omega_r=0.$5 earlier than $m\omega_\phi+k\omega_\theta+n\omega_r=0.$6 s before merger or compress $m\omega_\phi+k\omega_\theta+n\omega_r=0.$7 into the last $m\omega_\phi+k\omega_\theta+n\omega_r=0.$8–$m\omega_\phi+k\omega_\theta+n\omega_r=0.$9 s (Kuan et al., 2021).

A more realistic mode calculation with a solid crust and compositional stratification substantially revises part of this picture. In that framework, the canonical interface mode associated with the crust-core boundary vanishes and is replaced by compositional gravity modes with mixed gravity-interfacial character; lower-frequency modes induce gravitational-wave phase shifts smaller than $m\Omega(t_{\rm res})=\omega_\alpha$0 for the EOS considered. At the same time, nonresonant fundamental and crustal shear modes can trigger crust breaking already near the first gravity-mode resonance, while gravity-mode resonance concentrates strain at the base of the crust and may marginally crack it (Gao et al., 29 Aug 2025).

Nonlinear mode evolution further complicates the late inspiral. Three-dimensional simulations of tidally excited gravito-inertial modes in a stably stratified Boussinesq spherical shell find axisymmetric differential rotation, parametric instability of the $m\Omega(t_{\rm res})=\omega_\alpha$1 mode, and stellar spin-up that extends the resonance window for any given mode; in a constant-density cavity, the induced differential rotation may theoretically be large enough to amplify an ambient magnetic field to $m\Omega(t_{\rm res})=\omega_\alpha$2 G (Reboul-Salze et al., 31 Mar 2025).

The strongest currently known resonance effect is the retrograde-spin f-mode resonance seen in long-term numerical-relativity simulations of binary neutron stars. There the resonance begins at $m\Omega(t_{\rm res})=\omega_\alpha$3 Hz, the linear resonance window lasts about $m\Omega(t_{\rm res})=\omega_\alpha$4 ms, the stellar spin changes by $m\Omega(t_{\rm res})=\omega_\alpha$5 in the linear regime and by $m\Omega(t_{\rm res})=\omega_\alpha$6 in the later nonlinear regime, and the cumulative gravitational-wave phase shift reaches $m\Omega(t_{\rm res})=\omega_\alpha$7 radians by merger (Kuan et al., 2024).

Third-generation detectors are projected to be sensitive to at least part of this phenomenology. A fully Bayesian study with the Einstein Telescope finds that resonant modes can be identified in about one-third of the $m\Omega(t_{\rm res})=\omega_\alpha$8 loudest simulated binary-neutron-star signals in one year, with sensitivity to phase shifts as small as $m\Omega(t_{\rm res})=\omega_\alpha$9 for favourable events. The same study shows that neglecting resonances can bias the inferred tidal deformabilities (Gittins et al., 4 Jun 2026).

5. Resonances in oceans, planetary interiors, and horizonless compact objects

On Earth, linear shallow-water calculations for a global uniform-depth ocean show that the semidiurnal forcing frequency

$l=2$0

sweeps through successive global eigenfrequencies as Earth’s rotation slows. For the adopted present-day depth and friction coefficient, the lunar semidiurnal tide is less than $l=2$1 metre at present, but the model predicts a past resonance with $l=2$2 m at $l=2$3 h $l=2$4 m and a future resonance with $l=2$5 m at $l=2$6 h $l=2$7 m; the solar semidiurnal tide has corresponding peaks of $l=2$8 m at $l=2$9 h $m=\pm2$0 m and $m=\pm2$1 m at $m=\pm2$2 h $m=\pm2$3 m (Wei, 2020).

In an early-Earth core overlain by a basal magma ocean, the Moon’s tidal potential acts directly on the density contrast at the core–BMO interface, exciting a degree-2 interfacial mode rather than only deforming a bounding envelope. The response transitions from an equilibrium tide to an inertia–self-gravity wave as forcing frequency increases, and the two regimes are separated by a resonance that occurs when the semidiurnal forcing frequency matches the natural frequency of the interfacial mode. In the inviscid limit, the core flow is formally identical to the canonical elliptical vortex; finite viscosity regularises the resonance, generates oscillatory boundary layers, and yields enhanced core-boundary ellipticity, core-flow speed, magnetic Reynolds number, and instability metrics, particularly near resonance. Parameterised histories of lunar recession and BMO crystallisation imply that Earth likely crossed this resonance between $m=\pm2$4 and $m=\pm2$5 Ga (Kiernan et al., 15 Jun 2026).

An altogether different realization appears in fuzzball microstate geometries. There the long throat and reflective cap create a cavity that traps perturbations and supports metastable bound states, so the tidal response function

$m=\pm2$6

develops a sequence of real-frequency poles associated with resonant quasi-bound modes. By contrast, black holes absorb perturbations at the horizon and exhibit a smooth, featureless dynamical response without real-frequency poles in the corresponding toy model (Russo et al., 2024).

6. Diagnostics, misconceptions, and unresolved structure

Several recurring diagnostics distinguish resonant tidal responses from nonresonant or geometry-dominated behavior. In rotating convective envelopes, true modal resonances are identified by discrete peaks aligned with least-damped eigenfrequencies, close spatial matching between the forced response and the eigenmode, and scale separation in which kinetic energy is concentrated at low spherical-harmonic degree while enstrophy and viscous dissipation reside at higher degree in the beams. This was explicitly introduced to correct the common attribution of sharp inertial-wave peaks to wave attractors (Lin et al., 2021).

In compact planetary systems, the comparable diagnostics are libration of resonant angles, motion along resonance-defined equilibrium curves, and correlated offsets between adjacent period ratios. The offset trend is not a pairwise property alone: it is fixed by the three-body architecture and remains robust under different assumptions about the tidal factor or about which inner planets experience direct damping (Brasser et al., 2022, Charalambous et al., 2023).

In gravitational-wave problems, the corresponding signatures are localized phase kinks, discrete resonance crossings, or cumulative dephasings that are not replicated by adiabatic templates. The EMRI studies therefore advocate resonance-aware waveform generation, and the binary-neutron-star studies show that missing resonance physics can bias tidal deformabilities or obscure interior seismology (Gupta et al., 2022, Gittins et al., 4 Jun 2026).

Several limitations recur across the literature. Many stellar and planetary fluid calculations assume uniform rotation, homogeneous or Boussinesq structure, stress-free boundaries, and linear dynamics; the solar-type core calculations show that wave-driven mean flows and critical layers can erase the sharp resonances assumed in classical resonance-locking pictures, while neutron-star studies show that crust elasticity, stratification, magnetic fields, and nonlinear saturation can qualitatively alter which modes matter most (Guo et al., 2022, Gao et al., 29 Aug 2025). A plausible implication is that the observational relevance of any given tidal resonance depends not only on frequency matching, but also on whether dissipation preserves a narrow eigenfrequency-locked response or instead drives the system toward critical layers, manifold drift, or nonlinear mode coupling.