Thermal Disk Continua

A thermal disk continuum is not a single blackbody. It is the superposition of annuli with different effective temperatures, modified by scattering in the disk atmosphere and transported through the curved spacetime to an inclined observer. Its shape and normalization can constrain an apparent inner radius only when mass, distance, inclination, absorption, accretion state, detector response, and the inner-boundary model are carried together. The durable result is a forward model; a unique spin claim without those systematics lies outside this chapter.

Required background. Waves, Green Functions, and Geometric Optics supplies invariant intensity transport and the response-kernel viewpoint used to propagate emitted radiation to a detector.

Helpful background. Black-Hole Geometry and Mechanics supplies the Kerr orbital and inner-edge scales; Relativistic Matter and Compact Objects supplies radiative matter closure; Black-Hole Astrophysics supplies source states and inference context; and Computational Relativity supplies transfer-function and convergence checks.

The radial temperature field

For a steady, geometrically thin, optically thick Newtonian disk with a zero-torque inner edge $r_{\rm in}$, the radiative flux from one face is

$$F(r)=\frac{3GM\dot M}{8\pi r^3} \left[1-\left(\frac{r_{\rm in}}{r}\right)^{1/2}\right], \qquad \sigma T_{\rm eff}^4(r)=F(r).$$

Far outside the inner edge, $T_{\rm eff}\propto r^{-3/4}$. The bracket enforces the assumed boundary condition and makes the temperature vanish at $r_{\rm in}$; the maximum occurs outside it. In relativity, binding energy, angular momentum, frame dragging, photon capture, and redshift change the dimensionless radial flux. Page and Thorne derive the conserved-flux construction for a stationary, axisymmetric, equatorial thin disk in Page and Thorne 1974, §§ II–IV, pp. 499–506.

Those hypotheses are consequential. A thick or advective disk, appreciable stress across the nominal inner edge, a wind with radially varying $\dot M$, or irradiation of the photosphere changes the profile. Fitting a more flexible spectrum does not restore a violated thin-disk assumption.

Why a disk spectrum is broader than a blackbody

For a Newtonian flat disk seen at inclination $i$ and distance $D$, ignore atmosphere hardening and relativistic transfer for the moment. Adding the blackbody intensity from each visible annulus gives

$$F_\nu=\frac{2\pi\cos i}{D^2} \int_{r_{\rm in}}^{r_{\rm out}}B_\nu[T(r)]\,r\,dr.$$

There is one visible face in this observer flux, whereas the total luminosity includes both faces. In the region where $T\propto r^{-3/4}$, introduce $x=h\nu/[k_{\rm B}T(r)]$. Then $r\,dr\propto\nu^{-8/3}x^{5/3}\,dx$ while $B_\nu\propto\nu^3/(e^x-1)$. Consequently

$$F_\nu\propto\nu^{1/3} \int\frac{x^{5/3}}{e^x-1}\,dx.$$

When $k_{\rm B}T_{\rm out}\ll h\nu\ll k_{\rm B}T_{\max}$, the integral is approximately frequency independent and $F_\nu\propto\nu^{1/3}$. At still lower frequencies the entire finite disk approaches the Rayleigh–Jeans slope $\nu^2$; above the hottest annulus the spectrum turns down exponentially. The one-third slope is therefore an intermediate asymptotic result, not a fit to every thermal band.

At fixed $r/r_g$, spin, opacity, and $\dot m=\dot Mc^2/L_{\rm Edd}$, the flux scales as $\dot M/M^2\propto M^{-1}$, so $T_{\rm eff}\propto M^{-1/4}$. A $10^8M_\odot$ disk is about $56$ times cooler than a $10M_\odot$ disk at the same dimensionless rate and geometry. This explains why stellar black-hole disks can peak in X-rays while massive-black-hole disks peak in the ultraviolet.

Atmosphere hardening and relativistic transfer

Electron scattering shifts photons to a color temperature higher than the effective temperature. A common diluted-blackbody approximation is

$$T_{\rm col}=f_{\rm col}T_{\rm eff}, \qquad I_{\nu}^{\rm em}\simeq f_{\rm col}^{-4}B_\nu(T_{\rm col}),$$

where the dilution preserves bolometric flux. The hardening factor $f_{\rm col}$ is not universal: it depends on surface gravity, flux, composition, opacity, vertical dissipation, and illumination. Atmosphere grids replace the scalar approximation with angle- and frequency-dependent specific intensity.

Along a null geodesic, $I_\nu/\nu^3$ is invariant. If $g=\nu_{\rm obs}/\nu_{\rm em}$, the observed flux density is

$$F_{\nu_{\rm obs}}=\int_{\mathrm{image}}g^3 I^{\rm em}_{\nu_{\rm obs}/g}\,d\Omega_{\rm obs}.$$

The image integral contains gravitational redshift, Doppler boosting, lensing, and inclination, while the emitted intensity contains the atmosphere’s angular dependence. Returning radiation requires an additional disk-to-disk transport and reprocessing calculation that updates that emitted intensity; it is not supplied automatically by the disk-to-observer integral. McClintock, Narayan, and Steiner develop these inputs in McClintock et al. 2014, §§ 3.1–3.2, 4.1–4.2, and 5.1–5.4. The detector does not measure energy flux directly: convert to photon number flux before forward folding through its redistribution matrix and effective area.

What fixes an apparent inner radius

At fixed spectral shape, the normalization scales schematically as

$$N_{\rm disk}\propto \frac{R_{\rm in}^2\cos i}{D^2 f_{\rm col}^4},$$

before the full relativistic angular transfer is applied. Thus distance $D$, inclination $i$, hardening, and flux calibration are covariant with $R_{\rm in}$. Mass sets the conversion between a fitted length and $GM$. If $R_{\rm in}$ is then identified with the innermost stable circular orbit, a Kerr spin follows only under the disk-edge hypothesis. Reynolds summarizes why mass, distance, alignment, atmosphere, torque, and thickness must accompany such an inference in Reynolds 2021, §§ 3.3.1–3.3.3.

A useful internal check is constancy across independent disk-dominated epochs. If the source changes $\dot M$ while the inferred radius stays fixed, the data support a stable geometric scale. Failure of constancy can instead signal a changing edge, state contamination, an atmosphere error, or calibration drift; it cannot be assigned to spin evolution on observation timescales.

Observation and calibration boundary through 2026-08-18

The concrete evidence example is the paper-defined LMC X-3 RXTE/PCA sample of Steiner and collaborators. They grouped archival 1996–2009 pointings into roughly half-day spectra, used PCU-2, and selected 391 RXTE spectra within a larger multi-mission set Steiner et al. 2010, § II, § III.1, and Table 1. This page does not refit those events or claim a new posterior.

Inference layer	Published choice or data identity	Reproduction and limitation
Dataset and response	LMC X-3 archival RXTE/PCA PCU-2 spectra; historical mission response and Crab-based flux correction	The article defines the sample but does not provide a modern observation-level replay list; a rerun must record every ObsID, HEASoft version, PCA CALDB release, background model, response, and channel range
Selection	$0.05<L_{\rm disk}/L_{\rm Edd}<0.3$, fit-quality and Compton-tail cuts, with transition and hard-state spectra rejected	Selection precedes inference and depends on provisional mass and state classification
Forward model	Relativistic thin disk plus atmosphere hardening, returning radiation, limb law, absorption, and weak Compton tail	Inner stress, thickness, wind loss, alignment, $M$, $D$, $i$, and cross-calibration remain systematic uncertainties
Statistics	Joint physical parameters and nuisance priors evaluated through count-space residuals	A likelihood and posterior credible interval are conditional on response, calibration, prior ranges, and model family

The published sample found a nearly constant apparent radius within its selected regime Steiner et al. 2010, § IV and Figs. 1–2. That is a stability test, not a calibration-independent spin measurement. A defensible software reproduction would archive the observation list, response and background files, absorption prescription, atmosphere and relativistic-transfer versions, mass–distance–inclination priors, and scripts that regenerate posterior-predictive count residuals. It should vary torque, $f_{\rm col}$ or atmosphere grid, inclination, and luminosity cuts, and report the resulting parameter shifts.

Failure modes and cross-checks

Pile-up can harden a spectrum; background error matters when the source becomes faint; gain and effective-area errors distort curvature; absorption can imitate a soft-temperature change; and a Compton tail can steal thermal photons. Multi-epoch fitting should therefore share only physically stable quantities, leave epoch-dependent $\dot M$ and contamination free, and test cross-instrument normalization.

The strongest adversarial checks are a stable radius across selected epochs, consistent residuals in count space, recovery on simulated spectra passed through the same response, and agreement under alternative atmosphere, torque, and inclination assumptions. None alone proves that $R_{\rm in}=R_{\rm ISCO}$.

Common pitfalls

Fitting an unfolded spectrum. Unfolding depends on a reference model. The likelihood belongs in detector count space with the redistribution and effective area applied.

Converting radius to spin too early. A precise apparent radius can still carry systematic model dependence. The disk-edge identification and all external parameters must remain visible.

Exercise

At fixed observed color temperature, flux, distance, and inclination, a fit changes the hardening factor from $f_{\rm col}=1.70$ to $1.80$. By what factor does the inferred apparent radius change in the schematic normalization, and what does this teach about a spin claim?

Solution

Because $N_{\rm disk}\propto R_{\rm in}^2/f_{\rm col}^4$, holding the observed normalization fixed gives $R_{\rm in}\propto f_{\rm col}^2$. Hence

$$\frac{R_{\rm in,new}}{R_{\rm in,old}} =\left(\frac{1.80}{1.70}\right)^2=1.12.$$

The inferred radius grows by about $12\%$. Since the Kerr ISCO radius is mapped nonlinearly to spin, this atmosphere change need not produce a small or symmetric spin shift. The hardening prescription is therefore part of the inference, not a cosmetic correction.

References

• McClintock, J. E., Narayan, R., and Steiner, J. F. (2014). “Black Hole Spin via Continuum Fitting and the Role of Spin in Powering Transient Jets.” Space Science Reviews 183, 295–322. DOI.
• Page, D. N., and Thorne, K. S. (1974). “Disk-Accretion onto a Black Hole. I. Time-Averaged Structure of Accretion Disk.” The Astrophysical Journal 191, 499–506. DOI; Open PDF.
• Reynolds, C. S. (2021). “Observational Constraints on Black Hole Spin.” Annual Review of Astronomy and Astrophysics 59, 117–154. DOI; Open PDF.
• Steiner, J. F., McClintock, J. E., Remillard, R. A., Gou, L., Yamada, S., and Narayan, R. (2010). “The Constant Inner-Disk Radius of LMC X-3: A Basis for Measuring Black Hole Spin.” The Astrophysical Journal Letters 718, L117–L121. DOI; Open PDF.