---
title: Dynamical Tidal Love Numbers
url: https://www.emergentmind.com/topics/dynamical-tidal-love-numbers
type: topic
---

# Dynamical Tidal Love Numbers

Dynamical tidal Love numbers are frequency-dependent coefficients that quantify the linear tidal response of a gravitating body to a time-dependent external field. In the simplest formulation, they generalize the static or hydrostatic Love numbers by retaining the forcing frequency and therefore the coupling to internal degrees of freedom such as normal modes, inertial motions, or higher-derivative worldline operators. Across the contemporary literature, the same underlying idea appears in several technically distinct settings: Newtonian mode sums for rotating stars and giant planets, membrane and viscoelastic descriptions of satellites and solid bodies, Teukolsky- or Regge–Wheeler–based perturbation theory for compact objects and black holes, and point-particle effective field theory descriptions in which dynamical response is encoded by frequency-dependent Wilson coefficients [2103.06186], [2205.07577], [2110.12129], [2511.02372], [2507.22994], [2603.12331].

## 1. Definition and conceptual scope

A dynamical tidal Love number is a response coefficient at finite forcing frequency. In a rotating-fluid planetary formulation, for a given spherical-harmonic component \((l,m)\), the Love number is defined by
\[
k_{lm}=\frac{\delta\Phi}{U}\Big|_{r=R},
\]
where \(U\) is the imposed tidal potential and \(\delta\Phi\) is the self-gravitational response at the surface [2103.06186]. In rapidly rotating planets and stars, the same quantity is written componentwise as
\[
k_{\ell m}=\frac{\delta\Phi_{\ell m}}{U_{\ell m}},
\]
with forcing frequency
\[
\omega_m=m(\Omega_o-\Omega),
\]
so that the response is explicitly evaluated in the rotating frame [2110.12129].

In relativistic response theory, the induced multipole moments are commonly written as
\[
I_{\ell m}\left(\omega\right) =-\frac{\left(\ell-2\right)!}{(2\ell-1)!!}k_{\ell m}(\omega) r_+^{2\ell+1} \mathcal{E}_{\ell m}\left(\omega\right),
\]
and the frequency-dependent coefficient is expanded as
\[
k_{\ell m} \simeq \kappa_{\ell m} + {\rm i}\nu_{\ell m}\left(\omega - m \Omega\right) + \dots,
\]
with \(\kappa_{\ell m}\) the conservative part and \(\nu_{\ell m}\) the dissipative part [2407.07156]. Closely related EFT treatments instead use retarded response kernels \(K_\ell^{(E,B)}(\omega)\), whose low-frequency expansions define static, dissipative, and dynamical coefficients order by order in \(\omega\) [2511.02372].

The static or hydrostatic Love number is the zero-frequency limit only when that limit is regular. Several recent works stress that the \(\omega\to0\) limit of a dynamical response need not coincide with the result of solving the exactly static problem. For Kerr-like compact objects with reflectivity, the zero-frequency limit of the frequency-dependent tidal Love numbers is discontinuous, so the strictly static TLNs differ from the static limit of the frequency-dependent TLNs [2310.06023]. This suggests that “dynamical” and “static” are not merely different computational approximations, but can correspond to inequivalent response coefficients.

## 2. Response formalisms and mode decompositions

In Newtonian fluid theory, the dynamical Love number is naturally expressed as a sum over normal modes. For Jupiter, the mode-amplitude equation in the corotating frame leads to the central formula
\[
k_{lm}={2\pi \over 2l+1}\sum_\alpha {\bar Q_{\alpha,lm}^2\over \bar\varepsilon_\alpha (\bar\omega_\alpha-\bar\omega)},
\]
where \(Q_{\alpha,lm}\) is the tidal overlap, \(\varepsilon_\alpha\) a rotationally modified normalization, and \(\omega_\alpha-\omega\) the detuning from forcing [2103.06186]. The same logic appears in the Newtonian theory of spinning stars, where the effective Love number is written as a mode sum over the stellar oscillation spectrum, but only after a phase-space expansion is used to respect rotating-star symplectic orthogonality [2205.07577].

For a nonrotating star, the effective Love number takes the form
\[
k_{lm} = \frac{2\pi G}{(2l+1)R^{2l+1}} \sum_{\alpha'} \frac{I_{\alpha'}^2}{\mathcal A_{\alpha'}^2\left[\omega_{\alpha'}^2-(m\Omega_{\rm orb})^2\right]},
\]
where \(I_{\alpha'}\) is the mode mass multipole moment and \(\mathcal A_{\alpha'}\) the mode normalization [2205.07577]. In the rotating case, the full response contains contributions from both prograde and retrograde mode branches, and the Love-number expression acquires explicit \(l\)-mixing because rotation breaks spherical symmetry [2205.07577].

For rapidly rotating, centrifugally distorted planets and stars, the non-dissipative response is similarly written as
\[
k_{\ell m}^\alpha =\frac{2\pi}{(2\ell+1)} \sum_{\ell'} \frac{Q_{\ell m}^\alpha Q_{\ell' m}^\alpha} {\epsilon_{\alpha}(\omega_\alpha-\omega_m)} \left(\frac{ U_{\ell' m}}{ U_{\ell m}}\right),
\]
and the total Love number is \(k_{\ell m}=\sum_\alpha k_{\ell m}^\alpha\) [2110.12129]. The explicit \(\ell'\)-sum is the centrifugal-distortion correction: in a spherical body only \(\ell'=\ell\) survives, but in an oblate rotator a sectoral forcing component can drive a tesseral response [2110.12129].

In relativistic EFT language, the same frequency dependence is encoded in retarded kernels rather than explicit fluid modes. For Schwarzschild black holes, the odd- and even-parity kernels are expanded as
\[
K_{\ell}^{(B)}(\omega)= \frac{1}{\ell!}\left[ \frac{\lambda_{0,\ell}^B}{2}+i\omega r_s \frac{\lambda_{1,\ell}^B}{2}+ (\omega r_s)^2\frac{\lambda_{2,\ell}^B}{2}+\cdots\right],
\]
\[
K_{\ell}^{(E)}(\omega)= \frac{1}{\ell!}\left[ \lambda_{0,\ell}^E+i\omega r_s\lambda_{1,\ell}^E+ (\omega r_s)^2\lambda_{2,\ell}^E+\cdots\right],
\]
so that \(\lambda_{0,\ell}\) are static TLNs, \(\lambda_{1,\ell}\) dissipative coefficients, and \(\lambda_{2,\ell}\) dynamical Love numbers in the paper’s convention [2511.02372]. A plausible implication is that the mode-sum and EFT viewpoints are complementary low- and high-level descriptions of the same linear response structure.

## 3. Fluid planets and stars: rotation, resonances, and mode content

The most direct astrophysical application is to giant planets and rotating stars, where the Love numbers are dominated by a small subset of global modes. For Jupiter, the measured quadrupolar Love number
\[
k_2=0.565\pm 0.006
\]
is below the hydrostatic value
\[
k_2^{\rm (hs)}=0.590,
\]
and the finite-frequency response of rotation-modified \(f\)-modes was proposed as a simple explanation of the \(\sim4\%\) anomaly [2103.06186]. In an \(n=1\) polytropic model applied to Jupiter–Io, the calculation gives
\[
k_2^{\rm (hs)}=0.550,\quad k_2=0.528=0.960\,k_2^{\rm (hs)},
\]
so the dynamical correction is a 4% reduction relative to the hydrostatic value [2103.06186].

The same paper argues that the Love number is usually dominated by the response of the rotation-modified \(f\)-modes, while inertial-mode contributions are negligible because their overlaps are suppressed as \(Q_\pm\propto\Omega_s^2\) [2103.06186]. It also shows that sufficiently strong stratification in a large region of the interior can produce \(g\)-mode resonances that significantly alter \(k_{lm}\), turning dynamical Love numbers into probes of internal stratification [2103.06186].

For spinning Newtonian stars, the conceptual result is that rotation affects the dynamical and static problems differently. In a slow-rotation expansion,
\[
\text{dynamical } k_{lm} = k_{lm}^{(0)} + O(\Omega), \qquad \text{static } k_l = k_l^{(0)} + O(\Omega^2),
\]
so the dynamical tide is corrected already at first order in spin, whereas the static tide changes only at second order [2205.07577]. This first-order sensitivity arises because finite-frequency forcing distinguishes prograde and retrograde propagation; in the static limit those terms cancel [2205.07577].

Rapid rotation amplifies these effects. In centrifugally distorted \(n=1\) polytropes rotating up to \(\simeq90\%\) of breakup, \(f\)-modes remain dominant at nearly all rotation rates, but retrograde inertial modes can produce strong tesseral resonances, and at \(\gtrsim70\%\) of breakup prograde \(f\)- and \(i\)-modes mix strongly [2110.12129]. For \(\Omega/\Omega_d\simeq0.30\), the hydrostatic values listed for an \(n=1\) polytrope are
\[
k_{31}=0.18,\quad k_{22}=0.56,\quad k_{42}=0.32,\quad k_{33}=0.22,\quad k_{53}=0.16,\quad k_{44}=0.12,\quad k_{64}=0.10,
\]
whereas at the orbital frequencies of the Galilean satellites the tesseral Love numbers become much larger, for example
\[
k_{42}: 1.54\ {\rm (Io)},\ 3.81\ {\rm (Europa)},\ 9.65\ {\rm (Ganymede)},\ 29.89\ {\rm (Callisto)}
\]
[2110.12129]. This suggests that, in rapid rotators, tesseral Love numbers are especially sensitive to rotational coupling and inertial resonances.

A distinct but related result concerns Jupiter’s high-degree tesseral response. An analytic perturbative treatment of \(k_{42}\) shows that the hydrostatic \(k_{42}\) is dominated by the tidal response at \(\ell=m=2\) coupled into the spherical harmonic \(\ell,m=4,2\) by the planet’s oblate figure, and that high-degree tesseral Love numbers are dominated by lower-degree Love numbers when they are primarily hydrostatic [2112.05901]. The same abstract states that, after including the coupling from the well-understood \(\ell=2\) dynamical tides \((\Delta k_2 \approx -4\%)\), Jupiter’s hydrostatic \(k_{42}\) requires an unknown dynamical effect to produce a fractional correction \(\Delta k_{42}\approx-11\%\) in order to fit Juno’s observation within \(3\sigma\) [2112.05901].

## 4. Ocean worlds and viscoelastic bodies

In icy satellites with subsurface oceans, the dynamical correction is carried primarily by the ocean rather than the solid shell. In membrane theory for Europa- and Titan-like bodies, the Love numbers are written as
\[
\left( h_n , l_n , k_n \right) = \left( g y_1(R) , g y_3(R) , y_5(R)-1 \right),
\]
and the dynamical modification enters through
\[
\Lambda_T=\Lambda+\Lambda_\chi+\Lambda_\rho+\Lambda_\omega,
\]
with the key inertial term
\[
\Lambda_\omega = - q_\omega \, \frac{y_3^F(R_\varepsilon)}{y_1(R_\varepsilon)}
\]
[1504.04574]. For a homogeneous incompressible ocean over a rigid mantle,
\[
\Lambda_\omega = - \frac{q_\omega}{n(n+1)} \frac{n+1+n z^{2n+1}}{1-z^{2n+1}},
\]
and in the shallow-ocean limit
\[
\Lambda_\omega \sim - \frac{q_\omega}{n(n+1)} \frac{R_\varepsilon}{D},
\]
so the correction scales like \(1/D\) and becomes large when the ocean is shallow [1504.04574]. The same work identifies the resonance as the ocean’s surface gravity mode, and for Europa with a 10 km crust gives a representative resonant thickness
\[
D\sim 160\ \mathrm{m}
\]
[1504.04574].

In solid viscoelastic bodies, the relevant dynamical quantity is usually the dissipative combination \(k_2(\omega)\sin\epsilon_2(\omega)\), identified with the imaginary part of the complex Love number:
\[
{k}_l(\chi)\;\sin \epsilon_l(\chi)\;=\;-\;{\rm Im} [\bar{k}_l(\chi)].
\]
A damped mass-spring simulation of a Kelvin–Voigt self-gravitating sphere reproduces the expected “kink” shape of this function as a function of frequency: it is proportional to frequency at small \(\tilde\chi\), peaks near \(\tilde\chi\approx1\), and declines at higher \(\tilde\chi\) [1601.08222]. The inversion formula used in the simulations is
\[
[k_2 \sin \epsilon_2](\tilde \chi)\;= \;-\;\frac{\dot a}{3na} \left( \frac{M}{M^*} \right)\left( \frac{a}{R} \right)^5,
\]
so the dynamical Love response is inferred directly from tidal drift [1601.08222]. The paper states that the measured quality function is about 30% larger than the analytic prediction, with \(qf_{ratio}\sim 1.3\), showing that the frequency dependence is robust but the absolute normalization is only approximate [1601.08222].

A common misconception is that any paper relating Love numbers to an instability criterion is therefore computing dynamical Love numbers. The AdS-bubble study explicitly does not compute genuinely frequency-dependent Love numbers \(k_l(\omega)\); it works with static, time-independent perturbations and uses the sign of the resulting static TLNs to compare with a separate dynamical instability criterion [2606.13014]. This distinction is terminologically important because “dynamical” may refer either to the response coefficient itself or to an instability of the background system.

## 5. Compact objects and black holes

For compact objects, the dynamical problem is most naturally formulated in terms of asymptotic coefficients of wave-equation solutions. In Kerr-like compact objects with a reflective surface, the response is packaged into a frequency-domain function
\[
F_{\ell m}(\omega)=2k_{\ell m}(\omega)+i\omega \tau_{0}\nu_{\ell m}(\omega),
\qquad
k_{\ell m}=\frac12 \operatorname{Re} F_{\ell m},
\]
and is extracted from the decaying term in the asymptotic expansion of \(\Psi_4\) [2310.06023]. For non-rotating objects, the quadrupolar response can be written as
\[
k_{2}=\textrm{Re}\Bigg[\left(\frac{iM\omega}{30}\right)\left(1+16M^{2}\omega^{2}\right)\left(1+4M^{2}\omega^{2}\right) \left(\frac{1+\mathcal{R}(\omega)G(\omega)}{1-\mathcal{R}(\omega)G(\omega)}\right)\Bigg],
\]
where \(\mathcal R(\omega)\) is the reflectivity and \(G(\omega)\to1\) as \(\omega\to0\) [2310.06023]. For generic \(\mathcal R_0\neq1\), the zero-frequency dynamical TLN vanishes; only when
\[
\mathcal{R}(\omega)=1+iM\omega \mathcal{R}_{1}
\]
does the limit remain finite [2310.06023]. The same paper shows explicitly that the strictly static TLNs,
\[
k_{2}^{\rm polar}=\frac{8}{5\left(7+3\ln \epsilon \right)}, \qquad
k_{2}^{\rm axial}=\frac{32}{5\left(25+12 \ln \epsilon\right)},
\]
differ from the \(\omega\to0\) limit of the dynamical response [2310.06023].

For Kerr black holes, the static Love numbers vanish identically, but several finite-frequency analyses find nontrivial behavior. A low-frequency hidden-conformal treatment gives a response coefficient \(k_{\ell m}(\omega)\) proportional to a gamma-function ratio times
\[
\log\!\left(\frac{r_+-r_-}{r}\right),
\]
and argues that the corresponding dynamical tidal coefficients are generically non-zero and exhibit logarithmic behavior [2310.03660]. The Schwarzschild limit of that formula yields an \(\mathcal O(\omega^2)\) response structure proportional to \(\log(2M/r)\) [2310.03660].

A more directly perturbative Schwarzschild calculation in advanced null coordinates finds
\[
{}_{-2}\mathcal{F}_2(\omega)=\frac{i m\omega}{15} + \frac{2m^2\omega^2}{15} \left[ \ln\left(\frac{r}{2m}\right)+\frac{137}{21} \right],
\]
so that the linear term is purely dissipative and the quadratic term contains both the logarithmic running and a finite conservative piece [2507.22994]. The corresponding Love number inferred from the real part is
\[
{}_{-2}k_2 = \frac{m^2\omega^2}{15} \left[ \ln\left(\frac{r}{2m}\right)+\frac{137}{21} \right],
\]
although the same work cautions that proper EFT matching is still required before this can be treated as a definitive observable prediction [2507.22994].

The EFT-matched Schwarzschild analysis goes further in organizing the response. It defines the retarded kernels \(K^{(E,B)}_\ell(\omega)\), shows that the static black-hole coefficients vanish, and computes the leading nontrivial small-frequency terms. For example,
\[
\frac{45}{\pi}K^{(B)}_{\ell=2}(\omega)=i\omega r_s^6+\omega^2r_s^7\left(\frac{787}{2520}-\log (\mu r_s)\right),
\]
\[
\frac{45}{2\pi}K^{(E)}_{\ell=2}(\omega)= i\omega r_s^6+\omega^2r_s^7\left(\frac{797}{1260}-\log(\mu r_s)\right),
\]
so the \(i\omega\) term is dissipative while the \(\omega^2\log(\mu r_s)\) term is the universal logarithmic running of the conservative dynamical response [2511.02372]. This paper emphasizes that the finite constants are scheme dependent, whereas the logarithmic running is universal [2511.02372].

Extremal Kerr behaves differently again. In a Teukolsky–MST treatment, the low-frequency Love number is
\[
k_{\ell m} = - \frac{(-1)^s\,i\,m^{2\ell+1}\,(\ell-s)!(\ell+s)!} {2(2\ell)!(2\ell+1)!} \left(1+Q_\ell\,\omega+\mathcal O(\omega^2)\right),
\]
so the leading-order response is purely imaginary and therefore purely dissipative, while a conservative real part appears at order \(\omega\) [2412.19699]. A plausible implication is that finite-frequency black-hole response is best viewed as a hierarchy of conservative and dissipative coefficients rather than as a single yes-or-no statement about whether black holes “have Love numbers.”

## 6. Effective field theory, renormalization, and definitional issues

The EFT viewpoint turns dynamical Love numbers into Wilson coefficients. For compact stars, the point-particle action can be organized as
\[
S= S_{p.p}-M(GM)^4 \sum_{n=0}^{\infty}\sum_{l=2}^{\infty}\int d\tau \Bigg[(-1)^n(GM)^n\Lambda_{\omega^n}^{E,l}\frac{d^n\nabla^{i_3}..\nabla^{i_l}E^{ij}}{d\tau^n}\nabla_{i_3}..\nabla_{i_l}E_{ij} + E \leftrightarrow B\Bigg],
\]
so that \(\Lambda_{\omega^0}^{E,l}\) is the static TLN, odd-\(n\) coefficients are dissipative, and even-\(n\) coefficients are conservative dynamical TLNs [2603.12331]. In the electric quadrupole sector, the response function is written as
\[
\mathcal{F}^E(\omega) = \Lambda_{\omega^0}^E + i (GM\omega) H_{\omega^1}^E+ (GM\omega)^2 \Lambda_{\omega^2}^E,
\]
which makes the distinction between static, dissipative, and dynamical pieces explicit [2603.12331]. The same paper states that for non-viscous neutron stars the dissipative term vanishes, so the main new result is the NNLO conservative dynamical TLN and its renormalization-group equation [2603.12331].

A closely related PN/EFT treatment of gravitoelectric quadrupolar dynamical tides introduces explicit oscillator degrees of freedom \(Q_{\mu\nu}\) and finds that the first post-adiabatic tidal Wilson coefficient \(\kappa\) must be renormalized at 3PN. The renormalization-group equation is
\[
\beta(\kappa_{(a)}) = R \frac{d\kappa_{(a)}}{dR} = -\frac{214}{105},
\]
with solution
\[
\kappa_{(a)}(R)=\kappa_{(a)}(R_0)-\frac{214}{105} \log\left(\frac{R}{R_0}\right)
\]
[2308.01865]. This paper is careful that \(\kappa\) is the post-adiabatic Love number, not the ordinary adiabatic Love number \(\lambda\), so the renormalized coefficient is a higher-order dynamical response parameter rather than the static deformability itself [2308.01865].

A recurring conceptual issue is the relation between Love numbers and Green’s functions. One response-theory analysis stresses that the retarded Green’s function contains not only the instantaneous tidal response but also radiation-reaction effects, especially non-analytic tail terms such as \(\log\omega\), so dynamical Love numbers should not be naively identified with the full retarded Green’s function in asymptotically flat black-hole spacetimes [2407.07156]. By contrast, in BTZ black holes, where the absence of radiative modes eliminates tails, the Green’s function can be linked directly to the Love coefficient [2407.07156]. This suggests that, in four-dimensional asymptotically flat gravity, part of the current controversy is definitional rather than purely computational.

Another terminological point concerns low-frequency fluid stars in general relativity. The gravitomagnetic response of an irrotational fluid star is not a computation of modern frequency-dependent \(k_\ell(\omega)\), but it shows that even in an adiabatic slow-tide limit, internal motions driven by relativistic circulation conservation can qualitatively change the Love numbers. In particular, the gravitomagnetic Love numbers are positive in strict hydrostatic equilibrium and negative in the irrotational state [1504.06606]. This is not itself a dynamical Love-number calculation in the modern sense, but it identifies a low-frequency internal-motion effect that any truly dynamical theory should reproduce in the appropriate limit [1504.06606].

Taken together, these developments indicate that “dynamical tidal Love number” now denotes a family of related but not identical objects: mode-sum susceptibilities in fluid bodies, frequency-expanded Wilson coefficients in EFT, asymptotic response coefficients in relativistic perturbation theory, and sometimes low-frequency conservative pieces after radiation-reaction subtraction. The common content is finite-frequency tidal response; the precise observable depends on the formulation and on how source, response, dissipation, and running are separated [2407.07156], [2511.02372], [2603.12331].

Source: https://www.emergentmind.com/topics/dynamical-tidal-love-numbers