---
title: Binary-Perturbed Gravitational Atoms
url: https://www.emergentmind.com/topics/binary-perturbed-gravitational-atoms
type: topic
---

# Binary-Perturbed Gravitational Atoms

Binary-perturbed gravitational atoms are black-hole–boson-cloud systems in which a compact binary companion tidally drives the cloud away from its isolated superradiant evolution. In the nonrelativistic regime, the cloud is described by hydrogenic bound states around a Kerr black hole, but the companion’s time-dependent potential introduces three qualitatively distinct channels: resonant bound-to-bound transitions, bound-to-continuum ionization, and—when the companion is itself a black hole—accretion of cloud material onto the secondary. These processes backreact on the orbit and on the cloud occupation numbers, producing characteristic modifications of the emitted gravitational waves and, in some regimes, a direct gravitational-wave signal from the cloud transition itself [2112.14777], [2503.18121], [2412.12526].

## 1. Gravitational-atom framework

A gravitational atom consists of a Kerr black hole of mass $M$ and dimensionless spin $\tilde a$, surrounded by a bosonic cloud generated by superradiance. For a boson of mass $\mu$, the relevant coupling is the gravitational fine-structure constant
$$
\alpha \equiv \mu M \ll 1,
$$
and superradiance requires
$$
\omega < m\,\Omega_H, \qquad \Omega_H \equiv \frac{\tilde a}{2r_+}, \qquad r_+=M+\sqrt{M^2-a^2}.
$$
In this regime, the scalar field can be written in terms of a slowly varying envelope $\psi$ obeying
$$
i\partial_t\psi=\left[-\frac{\nabla^2}{2\mu}-\frac{\alpha}{r}\right]\psi,
$$
so the cloud admits hydrogenic bound states $\ket{n\ell m}$ with wavefunctions
$$
\psi_{n\ell m}(t,\mathbf r)=e^{-i(\omega_{n\ell m}-\mu)t}R_{n\ell}(r)Y_{\ell m}(\theta,\phi),
$$
and energies
$$
E_{n\ell m}\simeq \mu\left(1-\frac{\alpha^2}{2n^2}\right),
$$
up to higher-order corrections in $\alpha$ and spin-dependent splittings [2503.18121], [2412.12526].

The cloud’s characteristic size is the gravitational analog of a Bohr radius,
$$
r_c\equiv \frac{1}{\mu\alpha}\simeq \frac{M}{\alpha^2},
$$
and the cloud mass can grow to $\sim 0.1\,M$ in the fastest mode $\ket{211}$ [2305.15460]. Superradiant growth of an initial $\ket{211}$ cloud occurs on a timescale
$$
\tau_{\rm SR}\sim (\alpha^{9}\mu)^{-1},
$$
and clouds with $M_c\sim 10^{-4}$–$10^{-1}M$ can form around astrophysical black holes for $\mu\sim 10^{-13}$–$10^{-11}\ {\rm eV}$ [2112.14777]. The Kerr spin saturates at the critical value
$$
\tilde a_{\rm crit}=\frac{4m\alpha}{m^2+4\alpha^2},
$$
which fixes the endpoint of superradiant extraction in the hydrogenic regime [2503.18121].

## 2. Binary perturbation and tidal couplings

When the black hole–cloud system belongs to a binary, a companion of mass $M_*=qM$ on an orbit with separation $R_*(t)$, phase $\phi_*(t)$, and angular frequency $\Omega(t)$ generates a time-dependent tidal perturbation. In a Newtonian multipole expansion, the companion potential can be written as
$$
\Phi_{\rm ext}(r,\theta,\phi;t)= -q\alpha \sum_{\ell_*=0}^\infty\sum_{m_*=-\ell_*}^{\ell_*}\frac{4\pi}{2\ell_*+1}\,
Y_{\ell_*m_*}(\theta_*(t),\phi_*(t))\,Y_{\ell_*m_*}^*(\theta,\phi)\,F(r;\ell_*),
$$
which enters the Schrödinger Hamiltonian as $H'(t)=\mu\Phi_{\rm ext}$ [2412.12526]. For $R_*(t)\gg r_c$, the leading piece is quadrupolar,
$$
V_*(t,\mathbf x)=-m_c\,\alpha\sum_{m_*=-2}^{2}\mathcal E_{2m_*}(\phi_*)\,\frac{r^2Y_{2m_*}(\theta,\phi)}{R_*^3}+\cdots,
$$
and this leading-order perturbation mixes levels satisfying $\Delta m=\pm 2$ and $|\ell-\ell'|\le 2$ [2503.18121].

Matrix elements between cloud states are overlap integrals of radial hydrogenic functions and spherical harmonics. In the general expansion,
$$
\langle a|H'|b\rangle=\sum_g \eta^g e^{-ig\Omega t},
$$
where the Fourier coefficients $\eta^g$ depend slowly on $R_*(t)$ and encode the strength of the $g^{\rm th}$ harmonic [2412.12526]. The angular structure imposes the selection rules
$$
-m-m_*+m_b=0,\qquad \ell+\ell_*+\ell_b=\mathrm{even},\qquad |\ell_b-\ell|\le \ell_* \le \ell_b+\ell,
$$
while, in the quadrupolar approximation, the off-diagonal two-state coupling takes the dressed form
$$
V_{if}(t)=\eta\,e^{-i\Delta m\,\phi_*(t)},\qquad \Delta m=m_f-m_i
$$
[2503.18121].

This perturbative formulation provides a single framework for both discrete resonances and ionization. Bound-state mixing follows from matrix elements between discrete levels; ionization follows from the same tidal operator acting between a bound state and continuum states $\ket{k;\ell m}$ [2112.14777], [2412.12526].

## 3. Resonant bound-state transitions

A resonant transition occurs when the orbital driving matches a cloud level splitting. For two bound states $\ket{a}$ and $\ket{b}$, the resonance condition is
$$
(m_b-m_a)\Omega \approx E_b-E_a,
$$
or, in the notation used for generic harmonics,
$$
\Delta\epsilon=g\,\Omega_{\rm orb},
$$
with $g=m_b-m_a$ [2112.14777], [2403.03147]. In the two-level approximation, one passes to a dressed frame and obtains a Landau–Zener Hamiltonian. For quasi-circular inspiral with $\Omega(t)\simeq \Omega_0+\gamma t$, one form is
$$
H_{\rm LZ}(t)=
\begin{pmatrix}
(\Delta m/2)\gamma t & \eta \\
\eta & -(\Delta m/2)\gamma t-i\Gamma
\end{pmatrix},
$$
where $\Gamma=\mathrm{Im}\,\omega_f$ is the decay rate of the final state into the hole [2503.18121]. In the limit $\Gamma\to 0$, the standard transition probability is
$$
P_{\rm LZ}=1-\exp[-2\pi \zeta],\qquad \zeta=\frac{|\eta|^2}{|\Delta m|\,\gamma},
$$
equivalently $Z=\eta^2/(|g|\gamma)$ in the notation of the resonant-history analysis [2503.18121], [2403.03147].

The resonance radius has a characteristic scaling with $\alpha$,
$$
R_{\rm res}\sim M\,\alpha^{-p},
$$
with
$$
p=4\quad\text{(hyperfine)},\qquad p=10/3\quad\text{(fine)},\qquad p=2\quad\text{(Bohr)},
$$
depending on whether the transition changes $m$, $\ell$, or $n$ [2403.03147]. This hierarchy organizes the chronological sequence of resonances encountered during inspiral.

Generic eccentricity and inclination qualitatively modify whether a resonance operates. The resonant-history analysis introduces a backreaction parameter
$$
B=\frac{3M_cM}{q}\,\alpha^{4/3}\Omega_{\rm res}^{4/3},
$$
and shows that a floating resonance starts only if
$$
2\pi ZB \gtrsim 1.
$$
For eccentricity $e$, the effective chirp rate is modified by
$$
f(e)=\frac{1+\frac{73}{24}e^2+\frac{37}{96}e^4}{(1-e^2)^{7/2}},
$$
and the critical eccentricity satisfies $2\pi ZB=f(e_c)^{3/2}$; for inclination $\iota$, the coupling carries a Wigner-$d$ factor $\eta^g\propto d^{\ell_*}_{m_b-m_a,\;g}(\iota)$, which can suppress the resonance below threshold [2403.03147]. If the upper state has a nonzero horizon-decay rate $\Gamma$, the resonance can also break once
$$
|c_a|^2\lesssim \frac{\Gamma}{2ZB}.
$$

These results establish that binary-perturbed gravitational atoms are not described by a universal single-resonance picture. Whether a resonance is adiabatic, floating, broken by decay, or never activated depends on the sequence of crossings, the orbital geometry, and the cloud backreaction [2403.03147].

## 4. Ionization, dynamical friction, and accretion

Ionization is the bound-to-continuum channel in which the companion supplies enough orbital energy to unbind cloud quanta. In harmonic decomposition, the resonance condition for a continuum mode of wavenumber $k$ is
$$
\epsilon(k)-E_b \approx g\,\Omega(t),\qquad \epsilon(k)=\frac{k^2}{2\mu}.
$$
Integrating out continuum states in the Markov approximation yields the bound-state depletion law
$$
\frac{d\ln |c_b|^2}{dt}
=
-\sum_{\ell,g}\frac{\mu\,|\eta^{(g)}_{k_*b}|^2}{k_*}\,\Theta(k_*^2),
\qquad
k_*(t)=\sqrt{2\mu\,[g\Omega(t)+E_b]},
$$
and the ionized power
$$
P_{\rm ion}
=
\mu\sum_{\ell,g}\left[g\Omega(t)\frac{\mu\,|\eta^{(g)}_{k_*b}|^2}{k_*}\right]\Theta(k_*^2)\,|c_b|^2
$$
[2112.14777]. The same process can be derived from first-order perturbation theory and Fermi’s Golden Rule in the companion frame [2412.12526].

For typical parameters,
$$
M_c\lesssim 0.1M,\qquad \alpha\sim 10^{-1}\text{–}10^{-2},\qquad q\lesssim 10^{-2}\text{–}10^{-1},
$$
the orbital power lost to ionization satisfies
$$
\frac{P_{\rm ion}}{P_{\rm gw}}
\sim
\frac{\alpha^2(M_c/M)}{q}\left(\frac{r_c}{R_*}\right)^{\ell+2}\gg 1
$$
at separations $R_*\sim r_c$, so ionization can dominate over gravitational-wave radiation [2112.14777]. The ionization spectrum exhibits sharp jumps when $k_*\to 0$, producing threshold features or “kinks” that encode the cloud energy spectrum. In gravitational-wave frequency, the thresholds occur at
$$
\Omega^g=\frac{\alpha^3}{2gMn_b^2},
$$
or equivalently
$$
f_{\rm gw}^{(g)}=\frac{\alpha^2\mu}{\pi}\,g^{-1}n_b^{-2},
$$
depending on the notation adopted [2305.15460], [2412.12526].

A central result of the later literature is that ionization is equivalent to dynamical friction. The classical collisionless expression,
$$
P_{\rm df}
=
\frac{4\pi M_*^2\rho_c}{v}\left[\ln(v\mu b_{\max})+\gamma_e\right],
$$
matches the ionization power up to ${\cal O}(1)$ factors, so the cloud-induced backreaction on the orbit should be identified as dynamical friction [2305.15460]. In the small-$q$ limit,
$$
\frac{P_{\rm ion}}{P_{\rm gw}}
=
\frac{M_cM}{\alpha^5}\,\mathcal F(\alpha^2R_*/M),
$$
which explains why the effect is strongest when the companion traverses the cloud [2305.15460]. On open orbits, the additional energy dissipation can increase the dynamical-capture cross section by more than an order of magnitude, with typical $\sigma_{\rm tot}/\sigma_{\rm gw}\sim 10$–$10^2$ in the regime where cloud losses dominate [2305.15460].

If the companion is a black hole, cloud matter can also be accreted. In the nonrelativistic “fuzzy” regime $\mu M_*\ll 1$ and $kM_*\ll 1$, the mass-growth rate is
$$
\dot M_* = A_*\,\rho,
$$
with horizon area $A_*=8\pi M_*r_{+*}$ and local cloud density $\rho=M_c|\psi(R_*)|^2$ [2112.14777]. An equivalent particle-language expression is
$$
\dot M_* \simeq \pi b_{\max}^2\,v\,\rho,\qquad b_{\max}\sim \frac{4M_*}{v},
$$
and over the inspiral this accretion can increase the companion mass by $O(10\%)$ or more [2112.14777].

## 5. Backreaction on the binary and resonant histories

The orbit evolves under a combined energy budget rather than the vacuum quadrupole formula alone. In the Newtonian circular approximation,
$$
E_{\rm orb}=-\frac{qM^2}{2R},
\qquad
\dot E_{\rm orb}=-P_{\rm GW}-P_{\rm ion}-P_{\rm acc}-P_{\rm res},
$$
with
$$
P_{\rm GW}=\frac{32}{5}\frac{q^2M^5(1+q)}{R^5}
$$
[2503.23419]. Realistic waveform generation therefore requires integrating the coupled ODEs for the orbital separation and the cloud and black-hole masses, rather than adding a perturbative phase correction to an otherwise vacuum inspiral [2503.23419].

The different channels act on different timescales. Resonant transitions occur on
$$
\tau_{\rm res}\sim |\eta|/\sqrt{\dot\Omega}\sim (q\alpha)^{-1/2}\Omega^{-1},
$$
while ionization turns on when
$$
\Omega\sim \frac{|E_b|}{g}\sim \frac{\mu\alpha^2}{2n^2g},
$$
and then depletes the cloud on
$$
\tau_{\rm ion}\sim \frac{k_*}{\mu|\eta|^2}\lesssim 10^3 M,
$$
which can be much shorter than the gravitational-wave chirp time
$$
\tau_{\rm gw}\sim (q\alpha^2)^{-1}M
$$
[2112.14777]. This scale separation is one reason ionization often dominates the late cloud evolution once the companion enters the Bohr region.

Ionization and resonances affect orbital geometry differently. Ionization quickly circularizes the binary but barely affects the inclination angle: numerical solutions give $\tau_{\rm circ}\sim 10^3$–$10^5\,M$ when the semimajor axis is comparable to the cloud size, whereas $\tau_{\rm inc}\sim 10^6$–$10^8\,M$ and the total inclination change is $\lesssim {\cal O}(1^\circ)$ [2305.15460]. By contrast, resonant backreaction at generic eccentricity and inclination can produce dramatic changes in both $e$ and $\iota$, and the system’s long-term history bifurcates into two regimes [2403.03147].

In the first regime, “cloud destruction,” the first hyperfine or fine resonance is adiabatic and the cloud is fully transferred into a decaying state, after which the binary continues essentially as a vacuum inspiral but with finite kicks in eccentricity and inclination. In the second regime, “cloud survival,” near-counter-rotating systems or binaries entering too late for early resonances to become adiabatic retain the cloud until the Bohr region, where continuous ionization with $P_{\rm ion}\gg P_{\rm gw}$ produces a universal chirp together with nonadiabatic “sinking” jumps at each Bohr-resonant frequency. The associated phase jump satisfies
$$
\delta\Phi \gtrsim {\cal O}(10^3)\ {\rm rad},
$$
which is stated to be well above the LISA detectability threshold $\sim 2\pi$ [2403.03147].

## 6. Gravitational-wave phenomenology and detectability

The most immediate gravitational-wave imprint is dephasing of the binary inspiral. Because the orbit loses energy through ionization and, when applicable, through companion accretion, the gravitational-wave phase acquires an extra contribution
$$
\Delta\phi(f)=
-2\pi\int^f df'\,
\frac{P_{\rm ion}(f')+v^2\dot M_*(f')}{dE_{\rm orb}/df'},
$$
which is dominated by ionization when $P_{\rm ion}\gg P_{\rm gw}$ [2112.14777]. Each ionization kink changes the slope of $\phi_{\rm GW}(f)$ at the corresponding threshold frequency, yielding sharp spectroscopic features absent in vacuum templates [2112.14777]. In frequency-domain waveform models built from a reference IMR template, the environment-induced dephasing is
$$
\Delta\Phi(f)\equiv \Phi^{\rm cloud}(f)-\Phi^{\rm vacuum}(f),
$$
and amplitude corrections are typically subdominant and often neglected [2503.23419].

A second, distinct signal is the gravitational radiation emitted directly by a resonantly driven cloud transition. In the two-level approximation, the time-dependent cloud density contains an interference term $d_i^*d_f\psi_i^*\psi_f$ that sources a mass quadrupole. The resulting strains at distance $r_{\rm obs}$ are
$$
h_+(t)=h_0\frac{1+\cos^2\iota}{2}\,\mathrm{Re}\!\left[e^{-i2\Delta m\phi_*}Q(t)\right],
\qquad
h_\times(t)=h_0\cos\iota\,\mathrm{Im}\!\left[e^{-i2\Delta m\phi_*}Q(t)\right],
$$
with
$$
h_0=24M_c r_c^2\Omega_0^2/r_{\rm obs}
$$
[2503.18121]. In the stationary-phase approximation, the Fourier spectrum is peaked at
$$
f_p=2(\Omega_0-\zeta\Gamma)/\pi,
$$
with width $\Delta f\sim 4\Gamma/\pi$ in the small-$\zeta$ regime, and these observables can be combined with inspiral information to infer $\alpha=\mu M$, $\Gamma$, the cloud mass fraction, and thus the boson mass $\mu$ and black-hole spin $\tilde a$ [2503.18121].

Forecasts for detector sensitivity depend strongly on mass scale. For LISA-oriented systems, parameters
$$
\alpha\sim 0.1\text{–}0.3,\qquad
\mu\sim 10^{-14}\text{–}10^{-12}\ {\rm eV},\qquad
q\sim 10^{-2}\text{–}10^{-1}
$$
yield ionization kinks and an accretion-induced phase drift that are potentially observable with ${\rm SNR}\gtrsim 10$, and their combined detection would allow $\mu$ to be measured to percent-level accuracy [2112.14777]. For next-generation ground-based detectors, Fisher forecasts for the Einstein Telescope indicate that binaries with
$$
M\sim 40\,M_\odot,\qquad q\sim 0.1,\qquad M_c/M\gtrsim 10^{-3}\text{–}10^{-2},\qquad \alpha\sim 0.1\text{–}0.3
$$
can accumulate $\Delta\Phi\gg 2\pi$, with mismatches $\gtrsim {\rm few}\times 10^{-2}$ and relative uncertainties on $\alpha$ and $M_c$ at the $10^{-4}$–$10^{-2}$ level in Fisher analysis [2503.23419].

At much higher frequencies, primordial-black-hole gravitational atoms shift the discussion to the MHz–GHz band. A benchmark binary with
$$
M=6\times 10^{-10}M_\odot,\qquad \alpha=0.21,\qquad q=10^{-3},\qquad r=10\ {\rm kpc}
$$
has resonant orbital frequency $f_c\simeq 1.07\ {\rm GHz}$, within the ADMX Run-1 band $[0.645,1.4]\ {\rm GHz}$, but its characteristic strain is only $h_c(f_p)\lesssim 10^{-28}$, compared to $h_{\min}\sim 10^{-22}$; reaching the latter would require $r\lesssim 1\ {\rm AU}$ [2604.01407]. That analysis concludes that binary-driven high-frequency transients are generically too weak and too rare at astrophysical distances for current experiments, although isolated annihilation signals may be more promising targets [2604.01407].

A recurrent misconception is that binary perturbations are exhausted by Landau–Zener bound-state resonances. The literature instead supports a broader picture: ionization can dominate the energy budget, accretion can significantly change the companion mass, and the cloud state entering an observational band depends on the full resonant history, including orbital geometry and nonlinear backreaction [2112.14777], [2403.03147]. Under that more complete description, binary-perturbed gravitational atoms function as spectroscopic systems whose orbital and direct gravitational-wave signatures jointly probe ultralight bosons, black-hole spin extraction, and environmental modifications of compact-binary inspiral.

Source: https://www.emergentmind.com/topics/binary-perturbed-gravitational-atoms