Grey-Body Factors in Black Hole Physics

• Grey-body factors are frequency- and multipole-dependent transmission probabilities that determine how waves are filtered by a black hole’s effective potential.
• They are computed using a Schrödinger-like equation and high-order WKB methods, with numerical techniques validating their accuracy across various field types.
• They bridge Hawking radiation and quasinormal modes, offering insights into the stability and dynamics of black hole spacetimes.

Grey-body factors are the frequency- and multipole-dependent transmission probabilities that quantify how waves and particles produced near a black-hole horizon are filtered by the spacetime before reaching infinity. In the standard radial reduction, perturbations satisfy a one-dimensional Schrödinger-like equation in the tortoise coordinate, and the grey-body factor is the squared transmission amplitude through the associated effective potential barrier. They therefore encode the deviation of Hawking radiation from an ideal black-body spectrum, and, in recent WKB formulations, they are also closely related to quasinormal modes (QNMs) of the same background (Konoplya et al., 2024, Bolokhov, 19 Mar 2026).

1. Scattering definition and basic formalism

For static, spherically symmetric black holes, scalar, electromagnetic, Dirac, and gravitational perturbations reduce to

$$\frac{d^2\Psi}{dr_*^2} + \bigl(\omega^2 - V(r)\bigr)\Psi = 0,$$

with tortoise coordinate defined by

$\frac{dr_*}{dr}=\frac{1}{f(r)}$

or, in more general parametrized metrics, by $dr_*=(f(r)g(r))^{-1/2}dr$ (Bolokhov, 19 Mar 2026, Dubinsky et al., 2024). In the scattering problem one uses real frequency $\Omega$ and imposes asymptotic boundary conditions of the form

$$\Psi \sim T\,e^{-i\Omega r_*}, \qquad r_*\to-\infty,$$

at the horizon, and

$$\Psi \sim e^{-i\Omega r_*}+R\,e^{+i\Omega r_*}, \qquad r_*\to+\infty,$$

at spatial infinity, with transmission amplitude $T$ and reflection amplitude $R$. Conservation of the Wronskian gives

$|R|^2+|T|^2=1,$

so the grey-body factor is

$\Gamma_\ell(\Omega)\equiv |T|^2=1-|R|^2,$

with $\frac{dr_*}{dr}=\frac{1}{f(r)}$0 (Bolokhov, 19 Mar 2026, Konoplya et al., 2024).

Its physical role is immediate in Hawking radiation. The near-horizon emission is approximately thermal, but the flux at infinity is multiplied by $\frac{dr_*}{dr}=\frac{1}{f(r)}$1, so the spectrum is “grey” rather than black. Representative formulas used across the literature have the form

$\frac{dr_*}{dr}=\frac{1}{f(r)}$2

or, mode by mode, $\frac{dr_*}{dr}=\frac{1}{f(r)}$3 (Bolokhov, 19 Mar 2026, Malik, 19 Sep 2025, Konoplya et al., 2020).

For massive fields the asymptotic structure changes. In the massive scalar Hayward analysis, the effective potential tends to $\frac{dr_*}{dr}=\frac{1}{f(r)}$4 at infinity, $\frac{dr_*}{dr}=\frac{1}{f(r)}$5, and modes with $\frac{dr_*}{dr}=\frac{1}{f(r)}$6 are evanescent, so

$\frac{dr_*}{dr}=\frac{1}{f(r)}$7

(Dubinsky, 2 Nov 2025). In asymptotically de Sitter black holes the same logic applies, but the scattering region is the static patch between the event horizon and cosmological horizon; the transmission coefficient then measures propagation from the cosmological horizon toward the black hole or vice versa (Lütfüoğlu et al., 24 May 2026).

2. Effective potentials and field content

The effective potential $\frac{dr_*}{dr}=\frac{1}{f(r)}$8 is the central geometrical object. For a massless scalar field in a static, spherically symmetric background,

$\frac{dr_*}{dr}=\frac{1}{f(r)}$9

while for the electromagnetic field,

$dr_*=(f(r)g(r))^{-1/2}dr$0

For a massless Dirac field one obtains supersymmetric partner potentials

$dr_*=(f(r)g(r))^{-1/2}dr$1

with $dr_*=(f(r)g(r))^{-1/2}dr$2 in the four-dimensional conventions used in several applications (Bolokhov, 19 Mar 2026, Lütfüoğlu et al., 19 Sep 2025).

Massive scalar fields introduce an additional $dr_*=(f(r)g(r))^{-1/2}dr$3 term. In both the Hayward and generalized Proca analyses the scalar potential takes the form

$dr_*=(f(r)g(r))^{-1/2}dr$4

so the field mass raises the asymptotic value and typically also raises and broadens the barrier (Dubinsky, 2 Nov 2025, Lütfüoğlu et al., 24 May 2026).

For axial gravitational perturbations of static regular black holes supported by effective anisotropic fluids, several papers use

$dr_*=(f(r)g(r))^{-1/2}dr$5

which differs from the Schwarzschild Regge–Wheeler potential because the background geometry is not vacuum Einstein gravity in the usual sense (Malik, 26 Aug 2025, Dubinsky, 14 Sep 2025, Bolokhov et al., 9 Jul 2025). In five-dimensional Schwarzschild–Tangherlini spacetime, scalar-, vector-, and tensor-type gravitational perturbations have distinct Kodama–Ishibashi potentials, and the tensor sector exists only for $dr_*=(f(r)g(r))^{-1/2}dr$6 (Han et al., 18 Aug 2025).

The dilatonic GMGHS spacetime is exceptional because the gravitational, electromagnetic, and scalar perturbations form a coupled system whose diagonalization produces two axial wave equations and three polar wave equations,

$dr_*=(f(r)g(r))^{-1/2}dr$7

with the polar potentials sufficiently complicated that Ferrari, Pauri, and Piazza did not write them explicitly (Dubinsky, 2024). This already signals that grey-body factors are not universal functions of mass and spin alone; they are sensitive to the detailed effective barrier generated by the background fields.

3. Semiclassical and numerical computation

The dominant practical method is high-order WKB barrier penetration. For a single smooth peak with maximum $dr_*=(f(r)g(r))^{-1/2}dr$8 and curvature $dr_*=(f(r)g(r))^{-1/2}dr$9, the transmission probability is written as

$\Omega$0

or equivalently in forms differing by sign conventions for $\Omega$1, with

$\Omega$2

where the $\Omega$3 are higher-order WKB corrections involving higher derivatives of $\Omega$4 at the peak (Bolokhov, 19 Mar 2026, Dubinsky, 2 Nov 2025, Lütfüoğlu et al., 19 Sep 2025). Sixth-order WKB is repeatedly used for grey-body factors, while QNMs are often computed with higher-order WKB plus Padé resummation (Bolokhov, 19 Mar 2026, Lütfüoğlu et al., 19 Sep 2025, Malik, 26 Aug 2025).

Several studies validate WKB against independent calculations. Time-domain integration followed by Prony analysis is used to check QNM spectra in the Einasto-supported regular black hole and in the massive Hayward background (Bolokhov, 19 Mar 2026, Dubinsky, 2 Nov 2025). The proper-time asymptotic-safety analysis reports that WKB grey-body factors are accurate within a few percent over the parameter ranges studied (Lütfüoğlu et al., 19 Sep 2025). In five-dimensional Schwarzschild–Tangherlini spacetime, QNMs are computed numerically by Leaver’s continued-fraction method and grey-body factors are computed independently by direct numerical scattering, providing a stringent benchmark for the QNM-based reconstruction (Han et al., 18 Aug 2025). For non-commutative black holes under Einstein-coupled scalar fields, grey-body factors are obtained with the partial-wave method rather than WKB, and the QNM correspondence is then tested against that direct result (Fan et al., 20 Nov 2025).

The WKB approximation is nevertheless conditional. It is tailored to a single-peaked barrier and becomes less reliable for very low multipoles, for shallow or distorted barriers, and for multi-peak potentials. Several papers explicitly note that the QNM–grey-body correspondence can fail when the effective potential develops more than one peak, because secondary scattering and interference effects are then not captured by the standard one-barrier WKB expansion (Malik, 26 Aug 2025, Dubinsky, 2 Nov 2025).

4. Correspondence with quasinormal modes

Grey-body factors and quasinormal modes arise from the same master equation but with different boundary conditions. QNMs impose purely ingoing behavior at the horizon and purely outgoing behavior at infinity, producing a discrete complex spectrum $\Omega$5; grey-body factors impose ingoing behavior at the horizon and mixed incoming-plus-outgoing behavior at infinity, producing a continuous transmission function $\Omega$6 (Konoplya et al., 2024, Bolokhov, 19 Mar 2026).

A recent WKB development rewrites the barrier-transmission problem directly in terms of QNM data. In the eikonal regime,

$\Omega$7

where $\Omega$8 is the fundamental QNM of the same effective potential (Konoplya et al., 2024, Lütfüoğlu et al., 19 Sep 2025, Bolokhov, 19 Mar 2026). Beyond the leading eikonal approximation, the first overtone $\Omega$9 enters, and the resulting expression for the WKB scattering parameter depends only on $$\Psi \sim T\,e^{-i\Omega r_*}, \qquad r_*\to-\infty,$$0 and $$\Psi \sim T\,e^{-i\Omega r_*}, \qquad r_*\to-\infty,$$1, not on higher overtones (Konoplya et al., 2024, Malik, 26 Aug 2025, Skvortsova, 2024).

This correspondence has now been tested across a wide class of geometries. It is verified for asymptotically safe, proper-time quantum-corrected black holes with scalar, electromagnetic, and Dirac fields (Lütfüoğlu et al., 19 Sep 2025); for axial gravitational perturbations of the Hayward spacetime, where percent-level agreement is found for low multipoles and essentially exact agreement for higher ones (Malik, 26 Aug 2025); for Bonanno–Reuter regular black holes (Bolokhov et al., 9 Jul 2025); for GMGHS dilaton black holes (Dubinsky, 2024); for regular black holes with sub-Planckian curvature (Tang et al., 27 Mar 2025); for five-dimensional Schwarzschild–Tangherlini black holes, including the tensor channel (Han et al., 18 Aug 2025); and for rotating black holes in the KRZ parametrization, provided the perturbation equation is separable and the waves are non-superradiant (Konoplya et al., 2024).

The restriction to non-superradiant modes is fundamental in the rotating case. The WKB transmission formula used in the correspondence yields only positive grey-body factors, whereas superradiant scattering requires negative $$\Psi \sim T\,e^{-i\Omega r_*}, \qquad r_*\to-\infty,$$2 in the convention $$\Psi \sim T\,e^{-i\Omega r_*}, \qquad r_*\to-\infty,$$3. The rotating extension therefore does not describe superradiant amplification (Konoplya et al., 2024). A further limitation is that the derivation presumes a single smooth barrier. Where the potential has several peaks, the correspondence ceases to be quantitatively reliable (Malik, 26 Aug 2025).

One consequence repeatedly emphasized is that grey-body factors are often more stable than higher QNM overtones under small geometric deformations. Since the correspondence uses only the fundamental mode and first overtone, it naturally explains why transmission probabilities are much less sensitive than high overtones to localized near-horizon changes (Konoplya et al., 2024, Skvortsova, 2024).

5. Dependence on geometry, matter content, and couplings

Applications show that grey-body factors respond in a structured but not universal way to the field content and to the background geometry.

In the Einasto-supported regular black hole, grey-body factors are much less sensitive to the halo environment than QNMs. The main effect is a moderate increase of the effective potential near the horizon as the halo parameter grows, leading to a mild suppression of transmission at low frequencies, while at higher frequencies the transmission probabilities remain close to the Schwarzschild case (Bolokhov, 19 Mar 2026). A similar conclusion is reached for the Dymnikova regular black hole: variations of the quantum parameter affect the effective potential only near the horizon, leading to minor deviations of grey-body factors and absorption cross-sections from Schwarzschild (Dubinsky, 14 Sep 2025).

Field mass is a much stronger effect. In the quantum-corrected Hayward black hole, increasing the scalar mass $$\Psi \sim T\,e^{-i\Omega r_*}, \qquad r_*\to-\infty,$$4 lifts the asymptotic value of the potential to $$\Psi \sim T\,e^{-i\Omega r_*}, \qquad r_*\to-\infty,$$5, imposes the threshold $$\Psi \sim T\,e^{-i\Omega r_*}, \qquad r_*\to-\infty,$$6 for $$\Psi \sim T\,e^{-i\Omega r_*}, \qquad r_*\to-\infty,$$7, shifts the transmission peak toward higher frequencies, and strongly suppresses the low-frequency part of the spectrum (Dubinsky, 2 Nov 2025). In the asymptotically de Sitter generalized Proca black hole, increasing the scalar mass raises and broadens the scalar barrier, suppresses transmission at fixed frequency, and shifts efficient transmission to higher frequencies (Lütfüoğlu et al., 24 May 2026).

Quantum and nonlinear couplings can either suppress or enhance transmission, depending on how they reshape the barrier. In proper-time asymptotic safety, near-horizon quantum corrections make the effective potentials taller and broader near the peak, so grey-body factors are suppressed compared to the classical case, with the strongest suppression in the near-extremal regime and for low multipoles (Lütfüoğlu et al., 19 Sep 2025). In the Hayward axial gravitational sector, the quantum parameter $$\Psi \sim T\,e^{-i\Omega r_*}, \qquad r_*\to-\infty,$$8 raises the barrier and suppresses both grey-body factors and the absorption cross section (Malik, 26 Aug 2025). By contrast, in Einstein–Euler–Heisenberg black holes, the nonlinear QED coupling $$\Psi \sim T\,e^{-i\Omega r_*}, \qquad r_*\to-\infty,$$9 lowers the effective potential barrier for neutral scalar and Dirac fields and systematically enhances the transmission probability, while increasing the electric charge $$\Psi \sim e^{-i\Omega r_*}+R\,e^{+i\Omega r_*}, \qquad r_*\to+\infty,$$0 raises the barrier and reduces $$\Psi \sim e^{-i\Omega r_*}+R\,e^{+i\Omega r_*}, \qquad r_*\to+\infty,$$1 (Malik, 19 Sep 2025).

The dilaton black hole provides an especially strong suppression effect. For GMGHS perturbations, larger charge $$\Psi \sim e^{-i\Omega r_*}+R\,e^{+i\Omega r_*}, \qquad r_*\to+\infty,$$2 leads to smaller grey-body factors, and this suppression becomes substantial as the charge approaches its extreme value (Dubinsky, 2024). The paper attributes this to the increase of the real part of the fundamental QNM with $$\Psi \sim e^{-i\Omega r_*}+R\,e^{+i\Omega r_*}, \qquad r_*\to+\infty,$$3 while the imaginary part remains almost unchanged, together with the rise of the effective potential barrier (Dubinsky, 2024).

Environmental deformations appear to be comparatively weak unless they generate a secondary barrier comparable to the primary peak. In “Dirty Black Holes, Clean Signals,” far-zone perturbations modeling environmental matter have only a minor impact on grey-body factors and Hawking radiation unless the additional barrier becomes comparable in height to the main black-hole peak (Konoplya et al., 2 Jul 2025). The exact galactic-center solution with a Hernquist-like matter distribution leads to the same qualitative conclusion: in the regime fitting galactic behavior, the influence of the environment on classical and quantum radiation must be relatively small (Konoplya, 2021).

A broader parametrized analysis reaches an important structural conclusion. For generic parametrized spherically symmetric and asymptotically flat black holes, the primary parameter determining grey-body factors is the deviation of the event horizon radius from its Schwarzschild value, while the higher-order coefficients of the parametrization, which govern the near-horizon geometry, are much less significant (Dubinsky et al., 2024). This supports the empirical observation that transmission probabilities are considerably more stable against small deformations of the near-horizon geometry than quasinormal modes (Dubinsky et al., 2024, Konoplya et al., 2024).

6. Hawking radiation, absorption, and broader significance

Grey-body factors enter every partial-wave computation of black-hole absorption and Hawking emission. They determine effective absorption cross-sections through standard partial-wave sums and multiply thermal occupation factors in the energy and particle fluxes (Malik, 26 Aug 2025, Dubinsky, 14 Sep 2025, Lütfüoğlu et al., 24 May 2026). For that reason, even when the Hawking temperature is known, the observable spectrum cannot be specified without $$\Psi \sim e^{-i\Omega r_*}+R\,e^{+i\Omega r_*}, \qquad r_*\to+\infty,$$4.

The relative importance of temperature and grey-body filtering depends on the model. In the Dymnikova regular black hole, the Hawking radiation spectrum is governed mainly by the modified Hawking temperature, with grey-body factors providing only subleading corrections, because the transmission curves stay close to Schwarzschild while the temperature changes substantially (Dubinsky, 14 Sep 2025). In four-dimensional Einstein–Gauss–Bonnet gravity, by contrast, the coupling $$\Psi \sim e^{-i\Omega r_*}+R\,e^{+i\Omega r_*}, \qquad r_*\to+\infty,$$5 changes both the temperature and the transmission barriers; the study concludes that a positive coupling constant leads to much smaller evaporation rate and longer lifetime, while a negative one enhances Hawking radiation (Konoplya et al., 2020). In Einstein–Euler–Heisenberg theory, the enhancement of grey-body factors by nonlinear QED corrections implies a less strongly filtered Hawking spectrum for neutral scalar and Dirac fields than in Reissner–Nordström with the same charge (Malik, 19 Sep 2025). In the GMGHS case, the strong suppression of grey-body factors near extremality indicates a correspondingly suppressed gravitational and electromagnetic Hawking output (Dubinsky, 2024).

Grey-body factors have also become relevant to black-hole spectroscopy. Recent work relates the Fourier-domain ringdown amplitude to $$\Psi \sim e^{-i\Omega r_*}+R\,e^{+i\Omega r_*}, \qquad r_*\to+\infty,$$6 at frequencies above the fundamental mode, thereby linking real-frequency scattering data to gravitational-wave observables (Konoplya et al., 2024, Dubinsky et al., 2024). This does not collapse the distinction between scattering and resonance problems, but it reinforces the interpretation of grey-body factors as geometric observables that are closely related to, though not identical with, the QNM spectrum (Konoplya et al., 2 Jul 2025).

A recurring methodological and conceptual theme is robustness. Across regular black holes, parametrized deviations from Schwarzschild, galactic environments, and several quantum-corrected spacetimes, grey-body factors typically remain closer to the classical result than higher QNM overtones do (Konoplya et al., 2024, Dubinsky, 14 Sep 2025, Dubinsky et al., 2024). This suggests that transmission probabilities are controlled mainly by the integrated shape of the primary barrier, whereas higher overtones probe finer details of the near-horizon geometry. A plausible implication is that ringdown spectroscopy and grey-body-filtered Hawking spectra offer complementary diagnostics: the former is often more sensitive to localized structural changes, while the latter is often more stable and more directly connected to the real-frequency scattering problem.