---
title: Multi-Band Continuum Reverberation Mapping
url: https://www.emergentmind.com/topics/multi-band-continuum-reverberation-mapping
type: topic
---

# Multi-Band Continuum Reverberation Mapping

Multi-band continuum reverberation mapping is a reverberation-based methodology for active galactic nuclei in which correlated variability between multiple continuum bands is used to infer light-travel delays, and hence geometric scales, within the optical/UV continuum-emitting region of the central engine. In its modern form, it generalizes photometric reverberation mapping from line-dominated versus continuum-dominated bands to a broader program of measuring inter-band continuum lags, reconstructing the lag–wavelength relation, and testing accretion-disk structure while accounting for broad-line region contamination, diffuse continuum, and host dilution [1206.4486][1510.05648][2207.06432].

## 1. Physical basis and formalism

Continuum reverberation mapping is usually framed within the lamp-post reprocessing picture: a compact central source illuminates a geometrically thin, optically thick disk, and progressively larger radii re-emit at progressively longer wavelengths after a light-travel delay. In a standard Shakura–Sunyaev disk, the effective temperature scales as
\[
T(R) \propto R^{-3/4},
\]
which implies
\[
R(\lambda) \propto \lambda^{4/3},
\qquad
\tau(\lambda) \propto \lambda^{4/3}.
\]
A commonly used empirical parameterization is
\[
\tau(\lambda) = \alpha\left[\left(\frac{\lambda}{\lambda_0}\right)^{\beta} - 1\right],
\]
with the thin-disk expectation \(\beta = 4/3\) [1510.05648][2307.16568].

Several campaigns interpret \(\alpha\) or an equivalent \(R_{\lambda_0}\) as the normalization of the continuum-emitting size scale. In irradiated-disk treatments, the normalization depends on \(M\), \(\dot M\), \(\eta\), and any external heating term, and is often compared against a thin-disk prediction to quantify the now standard disk-size discrepancy [1510.05648][1912.10319].

A more general time-domain formulation writes the variable continuum at wavelength \(\lambda\) as a convolution of a driver \(X_{\rm b}(t)\) with a wavelength-dependent transfer function:
\[
X(\lambda,t) = \int_0^{\infty} \psi(\lambda,\tau) X_{\rm b}(t-\tau)\,d\tau.
\]
In this language, the measured lag is not a property of a single annulus, but of a responsivity-weighted transfer function whose centroid, peak, and shape depend on wavelength and geometry [2507.21671]. Simulations of thin-disk lamp-post transfer functions show that multi-band light curves are not only shifted in time but also distorted, because the transfer functions are skewed and broaden with wavelength; this “twisted” behavior implies that methods that ignore the transfer-function shape can underestimate accretion-disk sizes by \(\sim 30\%\), whereas CREAM recovers unbiased sizes in the LSST ugrizy simulations [1909.08638].

## 2. Observational architectures and filter design

The observational design of multi-band continuum reverberation mapping is driven by the expected lag amplitude, the need for high S/N light curves, and the requirement to minimize line contamination. A recurrent strategy is to define at least one reference band that is as line-free as possible, then to sample progressively redder continuum windows with either broad-band or custom narrow-band filters [1206.4486][2111.07385].

The broad-band photometric reverberation mapping campaign on NGC 4395 is formally a line-lag experiment, but it provides a methodological prototype. Using the Wise Observatory 1m telescope with SDSS \(g'\), \(r'\), and \(i'\), monitored for 9 consecutive nights with about one data point every \(\approx 10\)–15 min in each band, the \(i'\) band was taken as pure continuum, \(r'\) as continuum + H\(\alpha\), and \(g'\) as essentially continuum-dominated with weak H\(\beta\) contamination [1206.4486]. This same logic underlies continuum RM: establish a clean continuum reference, identify contaminated bands, and use multi-band information to separate components.

Other campaigns use specially designed narrow-band filters to suppress BLR contamination. PG 2130+099 employed four custom narrow-band filters placed in line-free continuum windows in the rest frame, with daily cadence over \(\sim 6\) months, and detected significant inter-band continuum lags of up to 3 days [2111.07385]. Mrk 509 used narrow bands centered at 4311, 5688, 6208, and 7018 Å with sub-day time sampling and typical flux uncertainties of \(1\%\), allowing lags of up to \(\sim 2\) days across the optical range to be measured [1912.10319]. MCG 08-11-011 extended this narrow-band strategy from \(\sim 4000\) to \(8000\) Å with daily sampling [2302.08402].

The AGN STORM campaign on NGC 5548 exemplifies the high-leverage UV–optical architecture: HST continuum windows at 1158, 1367, 1479, and 1746 Å, Swift UVOT bands, and ground-based \(BVRI\) plus \(ugriz\), giving contiguous coverage from 1158 Å to \(\simeq 9160\) Å with cadences from \(\sim 0.3\) to 1 day [1510.05648]. Survey-scale implementations then trade wavelength leverage for sample size: ZTF uses \(gri\) at \(\sim 3\)-day cadence over multiple years [2207.06432], while the Dan Zowada Memorial Observatory sample provides \(ugriz\) over four years with typical 1–3 day cadence [2511.05669].

Target selection can be decisive. In the IMBH candidate J0249−0815, high-cadence \(g,r,i,z\) monitoring on 2 m telescopes yielded no significant variability (\(<1.4\%\)) over 6–10 hr, primarily because the host fraction at 5100 Å was \(94.5\%\)–\(98.3\%\); the study argues that future IMBH campaigns should favor targets with higher AGN-to-host flux ratios and use tailored simulations before intensive monitoring [2405.11750].

## 3. Lag extraction, interpolation, and forward modeling

The most widely used lag estimators in multi-band continuum RM are the interpolated cross-correlation function, DRW-based forward models, and transfer-function reconstructions. ICCF remains standard: lags are usually defined as centroids above a threshold such as \(0.8\,r_{\max}\), with uncertainties estimated through FR/RSS Monte Carlo realizations [1510.05648][2111.07385]. Variants include zDCF, the von Neumann estimator, MICA, and multivariate correlation approaches [2111.07385][2302.08402].

JAVELIN models the driving light curve as a damped random walk and each responding band as a shifted, scaled, and top-hat-smoothed version of that driver. In NGC 5548, ICCF and JAVELIN lags agree well within uncertainties, although JAVELIN errors are systematically smaller, and ICCF values were adopted as conservative [1510.05648]. PyROA instead constructs a shared driving light curve and infers per-band shifts, scale factors, and smoothing through MCMC; in Mrk 876, PyCCF, JAVELIN, and PyROA gave consistent lag sequences once long-term trends were removed [2307.02630].

Detrending can be methodologically essential. In Mrk 876, slow variations broadened the cross-correlation function significantly; subtracting a moving boxcar average with a fiducial width of 100 days narrowed the CCFs and stabilized the inter-band lags, which then became consistent across methods [2307.02630]. By contrast, the NGC 4395 photometric RM experiment found that FR alone gave reliable lags with errors \(\approx \pm 1\) hour, whereas FR+RSS tended to overestimate lags by up to \(50\%\) in a densely sampled light curve [1206.4486]. This suggests that the optimal uncertainty model depends strongly on cadence and window function.

A separate methodological issue is whether the inference engine represents the true transfer function. In simulated LSST-like ugrizy light curves, PyCS and JAVELIN underestimate disk sizes by \(\sim 20\%\), JAVELIN thin-disk fits give sizes \(\sim 30\%\) too small, and only CREAM yields unbiased measurements because it explicitly models the skewed, wavelength-dependent transfer functions of a lamp-post irradiated thin disk [1909.08638]. This result directly affects the interpretation of high-S/N survey light curves, where transfer-function systematics can dominate over statistical errors.

## 4. BLR contamination, diffuse continuum, and host dilution

A central complication in multi-band continuum reverberation mapping is that broad-band “continuum” filters are often not pure continuum. In NGC 5548, spectral decomposition with LBT/MODS showed that the U and \(u\) bands contain \(\sim 19\%\) Balmer continuum, while \(r\) and \(R\) contain \(15\)–\(20\%\) H\(\alpha\) emission [1510.05648]. Simulations based on those fractions found upward biases of \(\sim 0.6\)–1.2 d in the U/\(u\) lags and \(\lesssim 0.6\)–1.0 d in \(r/R\), leading to the practical recommendation that U/\(u\) be excluded from disk-size fits in that dataset [1510.05648].

The ZTF \(gri\) survey extends this contamination problem to ensemble statistics. For 94 type 1 AGN with significant and consistent inter-band lags, the inferred continuum-emitting sizes are larger than the standard thin-disk prediction, the size scales as \(R_{\mathrm{CE}} \propto L^{0.48\pm0.04}\) with a scatter of 0.2 dex, and the observed relation can be explained by a photoionization model in which \(\sim 23\%\) of the total flux comes from diffuse BLR emission [2207.06432]. In that interpretation, optical multi-band continuum RM often measures a composite disk + diffuse BLR radius rather than a pure accretion-disk scale.

Direct evidence for this outer continuum component was presented for six AGNs using the ICCF-Cut method. In the Swift U band, an outer component was extracted with lags consistent with the predicted Balmer continuum lag, about half of the H\(\beta\) lag value, and similar results were obtained with the JAVELIN photometric reverberation mapping model for four of the six objects [2403.14091]. This establishes that Balmer diffuse continuum can contribute significantly to CRM lags in the rest-frame \(u/U\) band.

Broad-line contamination can also bias individual red bands. In Mrk 876, the \(i\)-band lag shows an excess attributed to variable H\(\alpha\) broad-line emission, and the flux–flux analysis finds both a variable spectrum following \(f_\nu \propto \lambda^{-1/3}\) and an \(i\)-band excess that points to strong variable H\(\alpha\) emission in that band [2307.02630]. In MCG 08-11-011, the very steep \(\tau \propto \lambda^{4.74}\) relation is explicitly interpreted as a possible consequence of diffuse continuum emission [2302.08402].

Host dilution is an additional systematic that can suppress the variability amplitude without directly altering the lag kernel. J0249−0815 is an extreme case: the host contributes \(98.3\%\) and \(94.5\%\) of the 5100 Å flux in the SDSS and P200 spectra, respectively, and the optical AGN fraction is \(\lesssim 5\%\), rendering hour-scale continuum RM infeasible at the achieved photometric precision [2405.11750]. This suggests that target vetting by AGN-to-host contrast is as important as cadence.

## 5. Empirical results from representative campaigns

The observational literature now spans both intensive single-object campaigns and large survey samples. The table lists representative results that have shaped the field.

| Campaign | Bandpass and sampling | Main result |
|---|---|---|
| NGC 5548 STORM [1510.05648] | 1158 Å to \(\simeq 9160\) Å; UV–optical; nearly daily cadence | Lag–wavelength trend broadly consistent with \(\tau \propto \lambda^{4/3}\); best-fit normalization implies a disk radius 3 times larger than the standard thin-disk prediction for \(L=0.1L_{\rm Edd}\) |
| Mrk 509 [1912.10319] | Custom optical narrow bands; sub-day sampling over two years | Inter-band delays up to \(\sim 2\) days; lags consistent with \(\tau \propto \lambda^{4/3}\); disk size factor 1.8 larger than standard thin-disk theory |
| PG 2130+099 [2111.07385] | Four line-free narrow bands; daily cadence over six months | Significant lags up to 3 days; derived disk sizes are a factor of 2–6 larger than theoretical disk sizes predicted from the AGN luminosity estimate |
| MCG 08-11-011 [2302.08402] | Five narrow bands from \(\sim 4000\) to \(8000\) Å; daily cadence | Lags up to \(\sim 7\) days; observed lags larger than standard thin-disk predictions by a factor of \(\sim 3\)–7; \(\tau \propto \lambda^{4.74}\) |
| Mrk 876 [2307.02630] | \(ugriz\) over 3.3 years with robotic observatories | Lag of around 13 days from \(u\) to \(z\); lags roughly follow \(\lambda^{4/3}\), but normalization is approximately a factor of 3 longer than expected; \(i\)-band excess attributed to H\(\alpha\) |
| ZTF parent/core samples [2207.06432] | \(gri\), \(\sim 3\)-day cadence, 94 AGN parent sample | Continuum-emitting sizes larger than standard thin-disk sizes; \(R_{\mathrm{CE}} \propto L^{0.48\pm0.04}\) with a scatter of 0.2 dex in the core+local sample |

These campaigns are complemented by source-specific results at the low- and high-accretion-rate ends. For the SEAMBH source IRAS 04416+1215, a thin-disk fit with \(\beta=1.33\) gives \(R_{\lambda_0}^{\rm (JAVELIN)} = 2.96^{+3.28}_{-2.03}\) light days at the \(g\)-band rest wavelength, approximately 4 times larger than the Shakura–Sunyaev prediction [2307.16568]. Over a larger baseline, a 4-year optical CRM study of 18 AGN found \(\tau_{5100} \propto L_{5100}^{0.4}\), broadly consistent with the \(\tau \propto L^{1/2}\) expectation for continuum reverberation from either the accretion disk or the BLR [2511.05669].

## 6. Interpretive tensions, scalability, and emerging applications

A persistent interpretive tension in multi-band continuum RM is that the shape of the lag–wavelength relation is often close to the thin-disk expectation, while the normalization is too large. NGC 5548, Mrk 509, PG 2130+099, Mrk 876, IRAS 04416+1215, and the ZTF sample all exhibit this pattern to varying degrees [1510.05648][1912.10319][2111.07385][2307.02630][2307.16568][2207.06432]. One explanatory route emphasizes genuine departures from standard thin-disk assumptions, including vertical structure, irradiation geometry, or super-Eddington effects; another emphasizes BLR diffuse continuum and line contamination; the ensemble evidence suggests that both classes of effects can operate simultaneously.

This ambiguity matters for downstream inferences. In NGC 5548, the optical continuum lags at wavelengths longer than the \(V\) band are equal to or greater than the lags of high-ionization-state emission lines such as He II \(\lambda 1640\) and \(\lambda 4686\), implying that the optical continuum-emitting source has a physical size comparable to the inner BLR [1510.05648]. The same paper argues that if H\(\beta\) lags are measured relative to the optical continuum, the true ionizing radius is larger by \(\approx 0.2\) d, so BLR radii can be underestimated when optical rather than UV drivers are used [1510.05648]. By contrast, the Mrk 509 analysis shows that a larger black hole mass due to the unknown geometry scaling factor \(f\) can reconcile the observed continuum lag normalization with thin-disk theory for that source [1912.10319]. This suggests that the disk-size discrepancy is not necessarily identical in origin from object to object.

On scalability, photometric continuum RM is intrinsically survey-compatible. The NGC 4395 broad-band study explicitly noted that LSST will monitor \(\ge 10^4\) quasars over a decade and that broad-band photometric RM can “increase by several orders of magnitude” the number of reverberation-mapped objects [1206.4486]. ZTF has already demonstrated that \(\sim 3\)-day cadence \(gri\) light curves can yield significant inter-band lags for a parent sample of 94 AGN [2207.06432], while the four-year 18 AGN program shows that medium-sized homogeneous optical samples can constrain the global \(\tau_{5100}\)–\(L_{5100}\) relation [2511.05669]. A plausible implication is that future progress will depend less on merely increasing sample size than on improving decomposition of disk, diffuse BLR continuum, and line contributions.

The method is also being generalized beyond single-SMBH disks. In a model for low-mass-ratio supermassive black hole binaries, the presence of a low-density cavity between an active mini-disk and a circumbinary disk produces a lag–wavelength relation that is flat at short wavelengths and transitions to \(\lambda^{4/3}\) at long wavelength, rather than following a single \(\lambda^{4/3}\) law throughout [2507.21671]. Applied to PG 1302−102, this SMBHB model reproduces the inter-band time lags and yields an inferred total mass and orbital period consistent with values derived from other methods [2507.21671]. This extends multi-band continuum RM from a disk-sizing tool to a structural discriminator for more complex accretion flows.

Multi-band continuum reverberation mapping therefore occupies a dual role. At one level it is a precision timing method for measuring \(\tau(\lambda)\) and testing \(\tau \propto \lambda^{4/3}\). At another, it has become a diagnostic of contaminating continua, BLR stratification, host dilution, transfer-function asymmetry, and non-standard accretion geometries. The contemporary literature does not reduce it to a single measurement of a single radius; it treats the continuum lag spectrum as a composite observable whose physical interpretation depends on filter design, cadence, decomposition strategy, and the chosen forward model.

Source: https://www.emergentmind.com/topics/multi-band-continuum-reverberation-mapping