Plasma, Matter, and Environmental Trapping

An environment can alter wave propagation, but it does not automatically turn black-hole superradiance into an instability. Plasma dispersion, a disk or torus, a companion’s tidal field, and a material cavity can each create turning points or reflection; the same environment can also absorb, dephase, or advect the field. The decisive question is global: does a mode satisfying the future-horizon condition return coherently with net gain after every loss channel is included? The claim boundaries below were reviewed on 2026-08-12 against the sources cited on this page; this is not an exhaustive census of every post-2020 environmental calculation. No universal verdict applies to “an astrophysical environment” without a specified field, background profile, boundary condition, and dissipation model.

Helpful background. Black-Hole Geometry and Mechanics supplies Kerr horizon data; PDEs and Green Functions for Gravitation supplies turning-point and non-self-adjoint spectral methods; Formulations of General Relativity supplies field equations in media; Relativistic Matter and Compact Objects supplies plasma and material scales; Black-Hole Perturbation Theory supplies flux-normalized scattering.

Local dispersion and global modes

For a cold, collisionless plasma in its local rest frame, a transverse electromagnetic disturbance is often introduced through the schematic dispersion relation

$$\omega_{\mathrm{loc}}^2=k_{\mathrm{loc}}^2+\omega_p^2, \qquad \omega_p^2=\frac{4\pi e^2n_e}{m_e}$$

in Gaussian units. The resemblance to a massive field is useful but limited. The plasma frequency depends on position and time, the medium selects a four-velocity, longitudinal and transverse responses differ, and collisions make the dielectric response complex. A constant Proca mass in vacuum is therefore not generally equivalent to an inhomogeneous magnetized plasma. The possible use of a plasma frequency as an effective photon mass in superradiance estimates, together with the qualifications on realistic astrophysical plasmas, is summarized in Brito, Cardoso, and Pani 2020, §5.8.3, p. 134.

The wavenumber and the distance in a tunneling integral must be conjugate. For a static radial patch with metric $ds^2=-f\,dt^2+f^{-1}dr^2+\cdots$, the proper distance and tortoise coordinate satisfy

$$d\ell=\frac{dr}{\sqrt f}=\sqrt f\,dr_*, \qquad \Omega=\frac{\omega}{\sqrt f},\qquad k_*=\sqrt f\,k_{\mathrm{loc}}.$$

Thus, neglecting geometric potentials only for this local illustration,

$$k_{\mathrm{loc}}^2=\Omega^2-\omega_p^2,\qquad k_*^2=\omega^2-f\omega_p^2.$$

Across a smooth forbidden interval, leading WKB leakage has the equivalent expressions

$$\mathcal T_{\mathrm{leak}}\sim \exp\left[-2\int |k_{\mathrm{loc}}|\,d\ell\right] =\exp\left[-2\int |k_*|\,dr_*\right].$$

Putting a proper wavenumber directly beside $dr_*$ misses the redshift. One turning point still does not make a bound state. If the density falls to zero outside a finite barrier, a positive frequency can escape into an outer propagating region. Reflection is then imperfect unless another boundary is imposed.

A dissipative slab with both losses computed

To isolate this issue, consider a one-dimensional scalarized medium, not Maxwell–fluid theory or magnetohydrodynamics. In a constant-redshift patch $f=f_0$, use proper time $\tau=\sqrt{f_0}\,t$ and proper position $x$. A real field obeys

$$u_{\tau\tau}+2\nu u_\tau-u_{xx}+\Omega_p^2u=0 \qquad(0\lt x\lt w),$$

with $\nu=\Omega_p=0$ outside the slab. Here $\nu\geq0$ is an amplitude damping rate. At both sharp interfaces impose continuity of $u$ and $u_x$. These are exact interface conditions for this specified equation; no WKB connection formula is applied to a discontinuous edge. The energy identity is

$$\partial_\tau\left[\frac{u_\tau^2+u_x^2+\Omega_p^2u^2}{2}\right] +\partial_x(-u_\tau u_x)=-2\nu u_\tau^2.$$

It establishes passivity and fixes the loss normalization. For $u=\operatorname{Re}[\psi(x)e^{-i\Omega\tau}]$, $\Omega\gt 0$, write the incident/reflected field as $e^{ikx}+r_s e^{-ikx}$ and the transmitted field as $t_s e^{ik(x-w)}$, with $k=\Omega$. Inside the slab,

$$q^2=\Omega^2-\Omega_p^2+2i\nu\Omega.$$

Choose the passive branch by continuation from $\nu\gt0$: $\operatorname{Im}q\geq0$, with positive real $q$ in the lossless propagating limit. Expressions at the cutoff $q=0$ mean their continuous limits. Matching the two fields and their derivatives gives

$$\begin{aligned} D_s&=(k^2+q^2)\sin(qw)+2ikq\cos(qw),\\ r_s&=\frac{(k^2-q^2)\sin(qw)}{D_s}, \qquad t_s=\frac{2ikq}{D_s},\\ \psi(x)&=(1+r_s)\cos(qx) +\frac{ik(1-r_s)}q\sin(qx). \end{aligned}$$

Both exterior channels have the same wavenumber, so the fractions of one incident flux are

$$\mathcal R=|r_s|^2,\quad \mathcal T=|t_s|^2,\quad \mathcal A_s=1-\mathcal R-\mathcal T =\frac{2\nu\Omega}{k}\int_0^w|\psi(x)|^2\,dx.$$

The last equality follows by integrating the dissipative current identity and is an independent way to compute absorption. The evanescent intensity exponent is $2w\operatorname{Im}q$; it describes attenuation inside the slab. It is not the complete transmission coefficient, which also includes reflection at both interfaces.

Width, phase and net gain

Choose a proper reference length $L_0$ and prescribe

$$f_0=0.8,\qquad \Omega_p L_0=2,\qquad \nu L_0=0.004,\qquad Z_H=0.004.$$

The inner boundary is an idealized, frequency-independent active reflector with amplitude $\sqrt{1+Z_H}$ and zero reflection phase. It represents a given small gain, not a newly calculated Kerr horizon. Put the slab a proper distance $L_c=32.46270879537195L_0$ from it. The real-axis cavity phase is

$$2\Omega L_c+\arg r_s(\Omega)=2\pi n,\qquad n=10.$$

Choose the continuous branch of $\arg r_s$ near $-2.09$. This length makes the width-$2L_0$ case phase matched at $\Omega L_0=1$. Keep the length fixed when doubling the width and solve the same phase equation. The computed values are rounded below.

Quantity	Width $2L_0$	Width $4L_0$
Phase-matched $\Omega L_0$	1	1.0000125863
Interior exponent $2w\operatorname{Im}q$	6.9282094	13.8563606
Leakage $\mathcal T$	0.0029295675	0.000002873320
Absorption $\mathcal A_s$	0.0023199303	0.0023068248
Return gain $\mathcal G_E$	0.9987295042	1.0016810631

Here the correct combination is

$$\mathcal G_E=(1+Z_H)\mathcal R =(1+Z_H)(1-\mathcal T-\mathcal A_s).$$

Leakage and absorption in this table refer to the same incident energy. They must be added, not multiplied as $(1-\mathcal T)(1-\mathcal A_s)$. An equivalent sequential convention uses conditional absorption $\mathcal A=\mathcal A_s/(1-\mathcal T)$.

This example has two useful failure tests. Recalculating the narrow slab with $\nu=0$ gives $\mathcal G_E\simeq1.0010518551$: dropping the imaginary response changes a damping prediction into growth. That recalculation also changes the dispersive scattering coefficients; it is not simply setting a measured absorption entry to zero. Conversely, the wider slab has gain above one but at the old $\Omega L_0=1$ its phase mismatch is $-0.0008317053$ radians. Gain alone has not selected a cavity mode.

A growth-time estimate and its boundary

For the analytic loop factor $\rho=\sqrt{1+Z_H}\,r_s(\Omega)e^{2i\Omega L_c}$, the complete complex-mode condition is $\rho=1$. We have computed its real-axis phase and losses, not that complex root. Applying the local linearization derived on Black-Hole Bombs and Reflecting Boundaries gives, at each phase-matched frequency,

$$\begin{aligned} a_1&=\operatorname{Re}\frac{r_s'}{r_s},\qquad \tau_g=2L_c+\operatorname{Im}\frac{r_s'}{r_s}\gt 0,\\ \Gamma_E^{(t)} &\simeq \sqrt{f_0}\,\frac{\log\mathcal G_E\,\tau_g}{a_1^2+\tau_g^2}, \qquad \tau_g^{(t)}=\frac{\tau_g}{\sqrt{f_0}} . \end{aligned}$$

The prime here means $d/d\Omega$. Local analyticity and a valid first-order frequency expansion are required. The amplitude-to-phase slope ratios are $1.2193\times10^{-4}$ and $4.0904\times10^{-5}$, both below a declared $10^{-3}$ diagnostic threshold. Thus neglecting $a_1^2$ hardly changes this linearized result; that diagnostic is not an error bound on a true complex pole or on plasma physics. The Killing-time delays are respectively $73.88065L_0$ and $73.87982L_0$, and the linearized energy rates are

$$\Gamma_E^{(t)} L_0\simeq -1.7207531\times10^{-5},\qquad +2.2734918\times10^{-5}.$$

For the growing case the simpler group-delay estimate $\tau_E=\tau_g^{(t)}/\log\mathcal G_E$ is $43985.2L_0$. Ten energy e-folds would take about $439852L_0$, below a prescribed environmental lifetime $10^6L_0$. This is a conservative opportunity test if the assumed environment persists and its later change terminates this mode. It does not determine an astrophysical time, formation probability, nonlinear saturation, or survival under an evolving medium. The narrow slab fails even this first test because its net rate is negative.

A common gain-and-loss language

Let $Z_H$ denote the energy gain on one superradiant horizon encounter. Let $\mathcal L$, $\mathcal A$, and $\mathcal M$ denote, respectively, leakage to infinity, irreversible environmental absorption, and loss from the selected mode through mixing. Define each as a conditional fraction of the energy left by the preceding process. For a fixed-delay one-channel model,

$$\mathcal G_E=(1+Z_H)(1-\mathcal L)(1-\mathcal A)(1-\mathcal M).$$

Coherent fixed-delay modes grow when $\mathcal G_E\gt 1$. For frequency-dependent returns use the complex mode equation or the justified weak-envelope approximation above. With weak losses and a common effective delay, the first-order balance becomes

$$\Gamma_{\mathrm{net}}\simeq \Gamma_{\mathrm{SR}}-\Gamma_{\mathrm{leak}} -\Gamma_{\mathrm{abs}}-\Gamma_{\mathrm{mix}}.$$

These relations are meaningful only after every term has been computed for the same mode and energy normalization. An effective-potential well supplies neither $Z_H$ nor the loss terms. Likewise, a positive local scattering gain measured for a wave incident from infinity does not determine the overlap of a quasi-bound environmental mode with the horizon.

A complex material response contributes to both sides of the problem. Its real part changes propagation and may strengthen reflection; its imaginary part absorbs energy and broadens the resonance. Increasing density can raise a barrier while simultaneously increasing collisional damping. Treating only the real dispersion relation systematically biases the result toward growth. The exact mirror model in Black-Hole Bombs and Reflecting Boundaries is valuable because its reflectivity can be varied independently and the leakage limit can be checked.

Environmental trapping results through 2026-08-12

The entries below classify mechanisms, not all members of a broad astrophysical category. “Formal trapping” means that a stated linear model admits a potential well or reflecting boundary. “Established growth” would add a globally solved mode with positive net rate and controlled losses.

Mechanism	Trapping input	Principal loss or ambiguity	What is established in the stated model	Status through 2026-08-12
Constant vacuum boson mass	Evanescence for $\omega_R\lt \mu_b$	Horizon absorption outside the superradiant band	Quasi-bound spectrum; scalar and vector growing sectors on Kerr	Controlled linear mechanism
Prescribed cold-plasma profile	Position-dependent cutoff or barrier	Escape where density falls, collisions, mode conversion	Formal turning points and model-dependent resonances	No environment-independent growth verdict
Material mirror or cavity	Imposed reflection coefficient and phase	Imperfect reflectivity and absorption	Growing fixed-delay modes or a solved complex return equation	Controlled toy boundary problem
Disk or torus	Density, pressure, magnetic, or geometric structure	Absorption, turbulence, advection, profile evolution	Possible reflection or mode modification for specified models	Strongly model-dependent
Companion tide	Time-periodic mode coupling	Sideband leakage, decoherence, orbital evolution	Coupling and resonant transitions can be computed	Not confinement by itself
Asymptotically AdS boundary	Global boundary condition	Horizon absorption and boundary-condition dependence	Growing sectors for particular rotating AdS black holes	Different global spacetime problem

The vacuum massive-field baseline and the distinction between confinement and horizon amplification are reviewed in Brito, Cardoso, and Pani 2020, §§5.4–5.8, pp. 99–135. For comparison, a mathematically specified reflecting AdS problem with a matched unstable scalar spectrum is given in Cardoso and Dias 2004, §III, pp. 3–5, PDF, for small, slowly rotating Kerr–AdS holes, long wavelength and fast boundary falloff. Neither source licenses replacing a dynamical disk or plasma by a perfect mirror without an error estimate.

Competing scope statements

One position starts from formal mode theory: if the environment produces a long-lived trapped state whose real frequency lies in the superradiant band, the horizon can make its imaginary part positive. This is correct as a conditional statement. It motivates searches for plasma-induced photon clouds, disk cavities, or companion-assisted resonances.

A more conservative position asks whether the same environment that creates the trap also destroys coherence or absorbs the wave faster than the horizon amplifies it. This is also correct as a conditional statement. Realistic plasmas can be inhomogeneous, collisional, magnetized, and time-dependent; disks advect and dissipate; companions introduce sidebands and a finite orbital evolution time. Neither position wins without a calculation that closes the gain-and-loss budget.

The apparent disagreement often comes from different questions. Existence of a formal quasi-bound solution in a prescribed, stationary, nondissipative profile is a spectral result. Survival of that solution in a self-consistent plasma or magnetohydrodynamic environment is a robustness result. An astrophysical population constraint additionally needs formation and observation models. These three levels should not be reported as one claim.

Environmental counterchecks

Three counterchecks should accompany a positive result.

Remove or open the trap. Continue the outer profile to its measured asymptotic behavior rather than ending it at a numerical wall. Increase the outer domain and impose outgoing conditions. A genuine leaky-mode frequency should converge; a box mode will track the artificial endpoint.

Restore passive losses. Add the imaginary part of the response, absorption by matter, and mode conversion permitted by the same model. Report the individual fluxes and verify that their sum equals the change of exterior energy. If an instability disappears for any plausible damping, the result is a conditional upper-growth estimate, not a robust prediction.

Evolve the background on its own timescale. Compare the growth time with density evolution, inflow, orbital motion, and magnetic or thermal reconfiguration. A stationary rate is not a reliable long-time prediction when its profile changes before several e-folds. Evolution can preserve or enhance growth as well as destroy it, so a frozen-profile lifetime test is not a universal necessary condition for amplification. Superradiance in astrophysical settings and the hierarchy of competing processes are surveyed in Brito, Cardoso, and Pani 2020, §§6.1–6.5, pp. 160–190.

What would change this assessment

The frontier conclusion is deliberately limited: controlled mirrors and vacuum massive fields demonstrate the mechanism, while generic plasma, matter, and companion environments have no single settled net-growth answer. Whether any particular realistic environmental trap remains unstable after self-consistent damping and evolution is unresolved; broad claims beyond a specified model are conjectural, and there is no universal verdict. Stronger conclusions would require a self-consistent plasma or magnetohydrodynamic calculation that simultaneously resolves the mode, horizon flux, outer leakage, collisional or kinetic damping, and background evolution; independent reproduction of a claimed environmental instability across codes and boundary treatments; or observations constraining the relevant density and dissipation profiles well enough to close the rate budget. The present references do not constitute that self-consistent calculation or a comprehensive post-2020 literature search.

Nonlinear saturation and the endpoint of a growing cloud remain outside this linear page. The next page compares linear growth with accretion, binary, environmental, and nonlinear timescales in Instability Timescales and Astrophysical Relevance.

Exercises

A coherent fixed-delay environmental cavity has horizon gain $Z_H=0.0040$ per round trip, leakage $\mathcal L=0.0015$, and material absorption $\mathcal A=0.0030$ conditional on surviving leakage. Neglect mode mixing. Does it grow? Explain why the answer would be misleading if absorption were omitted.

Solution

The complete multiplier is

$$\mathcal G_E=(1.0040)(1-0.0015)(1-0.0030) =(1.0040)(0.9985)(0.9970)\simeq0.9994865.$$

It is below unity, so the mode is damped even though the horizon encounter is superradiant and an outer trap exists. If absorption were omitted, one would obtain

$$\mathcal G_E^{(\mathcal A=0)}=(1.0040)(0.9985)=1.002494\gt 1$$

and predict growth. The sign reversal is exactly why the dispersive real part of an environmental response cannot be retained while its dissipative imaginary part is discarded. The calculation still assumes a coherent phase; mode conversion or temporal profile changes could add further losses.

References

• Brito, R., Cardoso, V., and Pani, P. (2020). Superradiance: New Frontiers in Black Hole Physics. Lecture Notes in Physics 971. Springer. DOI.
• Cardoso, V., and Dias, Ó. J. C. (2004). “Small Kerr–anti-de Sitter black holes are unstable.” Physical Review D 70, 084011. DOI. Open PDF.