Laura Ferrarese Rutgers University Lecture 4 Beyond the

Laura Ferrarese Rutgers University Lecture 4: Beyond the Resolution Limit SIGRAV Graduate School in Contemporary Relativity and Gravitational Physics

Lecture Outline 1. The Innermost Regions of Active Galactic Nuclei 2. Variability in Active Galactic Nuclei h h Continuum variability Emission line variability 3. Reverberation mapping: h making sense of the continuum-emission line connection 4. Potential problems…. 5. … and how to solve them: the Transfer Function 6. Where the Observations Stand

Schematic View of an AGN

Reverberation Mapping in a Nutshell d Measuring masses of SBHs using Reverberation Mapping is based on the assumption that the size r and velocity of the Broad Line Region (BLR) clouds are connected by a simple virial relationship: d According to the standard model, Broad Line Region (BLR) clouds are h many (107 -8, Arav et al. 1997, 1998, Dietrich et al. 1999) h small, dense (Ne ~ 109 -11 cm-3) h cold (Te ~ 2 104 K) h photoionized (Ferland et al. 1992) h localized within a volume of a few to several tens of light weeks in diameter around the central ionization source. d As such, the BLR is, and will remain, spatially unresolved at optical wavelengths even using space based instrumentation, and its size cannot be determined using conventional images.

Reverberation Mapping in a Nutshell d If the BLR is photoionized, the broad lines should respond to continuum variations. The line response contains a wealth of information regarding the spatial and kinematic structure of the BLR; therefore monitoring programs for AGNs started in the early ‘ 80 s in the attempt to quantify the nature of the continuum and emission line variations (if any!) d If everything works as planned, the time delay between continuum and line variation is simply (? ) related to the size of the BLR.

Reverberation Mapping in a Nutshell d Advantages: h reverberation mapping probes regions very close to the central engine (103 RSch), a factor at least 1000 smaller than allowed by “traditional’ methods which relay on resolved kinematics. h This leaves little doubt that the measured mass (if correct!) is in the form of a supermassive black hole d Disadvantages: h The observations are difficult. Close time monitoring at very closely spaced intervals and multiple frequencies is necessary. h For the virial approximation to be applicable, the kinematics must be dominated by gravity. The presence of outflows or not gravitational motions (to which gas might be prone) would undermine the method entirely. h The geometry of the BLR are not known.

Observational Requirements for Monitoring Programs d Temporal Sampling h observations must be closely spaced in time relative to the physical timescale of interest (generally a problem for the early monitoring programs) h difficulties in scheduling the observations d S/N of the data h S/N must be high relative to the magnitude of the flux variations (i. e. Fvar must be >> 0) h e. g. S/N=30 is necessary to detect 10% variations in continuum flux at the 3 confidence level. d Aperture effects h variations in seeing at the time of the observations, as well as pointing and guiding errors can cause variations in the amount of light entering the spectrograph, since both the NLR and the host galaxy are extended. This can cause spurious spectral variations. h Datasets must be as homogenous as possible, ideally using a single instrument in a stable configuration

Variability in AGNs d AGNs are variable at all wavelengths at which they have been studied, not only in the continuum, but also in the emission lines. h Typical quasars vary at the 0. 3 - 0. 5 mag level over timescales of a few months, with extreme cases varying on timescales as short as a few days. h The variability in Seyfert galaxies is less dramatic and was not discovered until the late ‘ 60 s. h Causality arguments imply that the emitting region is less than a few light days across d Periodicity in the light curve have been searched for but never found: variations are aperiodic and have variable amplitudes

Emission Line Variability d Broad Emission lines in AGN spectra can vary in both flux and profile. d Narrow lines fluxes do not vary! This is due to the fact that in the NLR both the light crossing time and the recombination time are large (>100 years), therefore short-term variability is smeared out.

Line-Continuum Variations H line flux against the continuum flux measured at the same time (left) and 15 days earlier (right), for the Seyfert 1 galaxy Mrk 335. The emission line fluxes are better correlated with the earlier rather than current continuum fluxes (Peterson et al. 1998, Ap. J 501, 82)

Observational Results d AGNs with lags for multiple lines show that highest ionization emission lines respond most rapidly ionization stratification

The Basis of Reverberation Mapping d The fact that emission lines vary in response to the optical/UV continuum variation immediately implies that: h The line emitting clouds are close to the continuum source h the line emitting clouds are optically thick h The ionizing continuum ( < 912Å) is closely related to the observable optical/UV continuum. d Therefore, our hopes are realized: by characterizing the emission line response to continuum variations, the kinematics and geometry of the BLR can be constrained: h the time delay between continuum and emission line variations are ascribed to light travel time effects within the BLR: the emission lines ‘echo’ or ‘reverberate’ to the continuum changes (Blandford & Mc. Kee 1982, Ap. J, 255, 419).

Reverberation Mapping Assumptions d The continuum originates in a single central source. Typical scalelengths are: Accretion disk (for 107 – 108 M SBH): 1013– 14 cm ¢ Broad Line Region: 1016 cm h To all effects, as seen from the BLR, the continuum source can be treated as pointlike h The continuum is not required to be emitted isotropically (although isotropy is usually assumed) ¢ d The most important timescale is the light-travel time. h the cloud response to a change in the continuum flux is instantaneous. ¢ Light travel time: ¢ Timescale to re-establish photoionization equilibrium: ¢ Timescale it takes a Lyman photon to diffuse outward through the BLR: ne = electron density; U = ionization parameter; B= recombination coefficient; Rion = depth to which the BLR is completely ionized

Reverberation Mapping Assumptions d The structure of the BLR does not change on the variability time scale (or the timescale over which the experiment is conducted). h Dynamical (cloud-crossing) time: d There is a simple, though not necessarily linear, relationship between the observed continuum and the ionizing continuum.

Reverberation Mapping Assumptions d Once the (responsivity weighted) size r of the BLR is known, the AGN central mass can be obtained through the virial relationship: where f is a dimensionless factor of order unity that depends on the geometry and kinematics of the BLR, and is the emission line velocity dispersion. d The velocity width of the lines is measured in the rms spectrum: h the rms spectrum only contains information on the variable part of the lines; constant components do not contribute.

Potential Problems d What is f? h circular, coplanar orbits: mean-square line-of-sight velocity is GMsin 2 i/(2 r), therefore f=2/sin 2 i. f could therefore take any value between 2 and . h random, isotropic circular orbits: mean-square line-of-sight velocity is GM/(3 r), therefore f=3 h random, isotropic parabolic orbits: mean-square line-of-sight velocity is 2 GM/(3 r), therefore f=3/2 d These potential problems add to the systematics arising from the (generally) inadequate temporal sampling of the observations, and the (generally) short duration of the experiments.

Potential Problems: the Virial Hypothesis d How can we test the virial hypothesis? h If the motion of the gas is gravitational, using BLR sizes and velocity derived from different emission lines in the same AGN must produce the same estimate of the central mass. h NGC 5548: highest ionization lines have smallest lags and largest Doppler widths, such that virial product r V 2 is constant. Virial relationship with M = 6 107 M. 1989 data from IUE and ground-based telescopes. 1993 data from HST and IUE.

Potential Problems: the Virial Hypothesis d There a total four AGNs for which lag measurements for multiple emission lines exist, all supporting the virial approximation (Onken & Peterson 2002) h NGC 7469: 8. 4 106 M h NGC 3783: 8. 7 106 M h NGC 5548: 5. 9 107 M h 3 C 390. 3: 3. 2 108 M

Potential Problems d To summarize, all of our problems would be solved if the geometry and kinematics of the broad line region were completely determined: h One of the major remaining mysteries of AGN astrophysics h We need to know this to understand systematic uncertainties in AGN masses. d Can we determine the BLR geometry and kinematics from the observations? YES! d BUT… this will require a leap in data quality. h Accurate mapping requires a number of characteristics (nominal values follow): ¢ ¢ High time resolution ( 0. 2 day) Long duration (several months) Moderate spectral resolution ( 600 km s-1) High homogeneity and signal-to-noise (~100) d Given these data, we could not just restrict ourselves to measuring time lags, but we could measure the complete transfer function.

The Transfer Function d The transfer function determines the relation between continuum and emission lines variations: Emission-line flux at line of sight velocity Vz Transfer Function Continuum Light Curve d The transfer function is simply the time-smeared emission-line response to a function outburst in the continuum. In other words, the transfer function can be interpreted as a ‘velocity delay map’. d Solving for the transfer function is a classical inversion problem. In practice, it requires extremely well sampled, high quality data.

The Transfer Function Integrate over time delay to get the line profile Integrate over velocity to get the delay map d In the best case, the data so far only allows us to solve for the velocity independent (or 1 -d) transfer equation where both ( ) and L(t) represent integrals over the emission line width:

Transfer Function: Keplerian Disk d Transfer function for a thin keplerian disk at a 45 degree inclination.

Transfer Function: BLR Outflows d Transfer functions in the case of a BLR in spherical (left) and biconical (right) outflow.

Complex Transfer Functions d Recovering complex transfer functions requires mapping at multiple emission lines.

Recovering Velocity-Delay Maps from Real Data Notice little response from material along our line of sight to the continuum source Transfer function recovered from the H emission in NGC 5548. Caution should be exercised since the data spans a period longer than the BLR dynamical timescale. Transfer function recovered from the CIV emission in NGC 5548. The data has been interpreted as 1) evidence of no outflows; 2) evidence of radial outflows; 3) evidence of radial inflow (!).

Observational Results d Although no experiment yet has recovered a reliable velocity-delay map, emission-line lags have been measured in 37 AGNs, in some cases for multiple emission lines. The H response in NGC 5548 has been measured for 14 individual observing seasons. Measured lags range from 6 to 26 days.

Reverberation Mapped AGNs From Kaspi et al. 2000, Ap. J, 533, 631

Mass-Luminosity Relationship d The measured masses correlate, although with very large scatter, with the continuum luminosity, in the sense that brighter AGNs have larger SBHs. M L 0. 3± 0. 1 QSOs (Kaspi et al. 2000) Seyfert 1 s (Wandel, Peterson, Malkan 1999) Narrow-line AGNs NGC 4051 (NLS 1)

“Secondary” Mass Estimators d Reverberation mapping opens the way to calibrate a “secondary” mass estimators since, to first order, we expect the broad line region size to correlate with the ionizing continuum luminosity: h. Photoionization equilibrium models are parameterized by the shape of the ionizing continuum, the elemental abundances, and the ionization parameter U: where Q(H) is the number of hydrogen ionizing photons ( =13. 6 e. V) emitted per second by the central source: h. U characterizes the ionization balance within the cloud, since Q(H)/r 2 is proportional to the number of ionizations occurring per unit area, while ne is proportional to the recombination rate. h. To first order, AGN spectra all look alike, i. e. they have the same ionization parameter and electron density (typical values are: Q(H) ~ 1054 h 0 -2 photons s-1; ne ~ 1011 cm-3; U ~ 0. 1). Therefore, we expect

BLR Scaling with Luminosity d This is close to what we observe! For the 37 AGNs which have been reverberation mapped, the BLR radius, measured from the H time lag, correlates (although with large scatter) with the continuum luminosity. r(H ) L 0. 6± 0. 1 QSOs (Kaspi et al. 2000) Seyfert 1 s (Wandel, Peterson, Malkan 1999) Narrow-line AGNs NGC 4051 (NLS 1)

Suggested Readings d Review: Peterson, B. M. 2001, “Variability of Active Galactic Nuclei”, in The AGN Connection, World Scientific (astro-ph/0109495). d Criticism: Krolik 2001, Ap. J, 551, 72 Starburst-