Anna Pancoast

Modelling reverberation mapping data - II. Dynamical modelling of the Lick AGN Monitoring Project 2008 data set

Author(s): Pancoast, A; Brewer, BJ; Treu, T; Park, D; Barth, AJ; Bentz, MC; Woo, JH | Abstract:

SPACE TELESCOPE AND OPTICAL REVERBERATION MAPPING PROJECT. II. SWIFT AND HST REVERBERATION MAPPING OF THE ACCRETION DISK OF NGC 5548

Recent intensive Swift monitoring of the Seyfert 1 galaxy NGC 5548 yielded 282 usable epochs over 125 days across six UV/optical bands and the X-rays. This is the densest extended active galactic nucleus (AGN) UV/optical continuum sampling ever obtained, with a mean sampling rate <0.5 day. Approximately daily Hubble Space Telescope UV sampling was also obtained. The UV/optical light curves show strong correlations (r max =0.57-0.90) and the clearest measurement to date of interband lags. These lags are well-fit by a τ ∝ λ4/3 wavelength dependence, with a normalization that indicates an unexpectedly large disk radius of ∼0.35 ± 0.05 lt-day at 1367 A, assuming a simple face-on model. The U band shows a marginally larger lag than expected from the fit and surrounding bands, which could be due to Balmer continuum emission from the broad-line region as suggested by Korista and Goad. The UV/X-ray correlation is weaker (rm < 0.45) and less consistent over time. This indicates that while Swift is beginning to measure UV/optical lags in general agreement with accretion disk theory (although the derived size is larger than predicted), the relationship with X-ray variability is less well understood. Combining this accretion disk size estimate with those from quasar microlensing studies suggests that AGN disk sizes scale approximately linearly with central black hole mass over a wide range of masses.

The Galaxy Zoo survey for giant AGN‐ionized clouds: past and present black hole accretion events

Some active galactic nuclei (AGN) are surrounded by extended emission-line regions (EELRs), which trace both the illumination pattern of escaping radiation and its history over the light-travel time from the AGN to the gas. From a new set of such EELRs, we present evidence that the AGN in many Seyfert galaxies undergo luminous episodes 0.2–2×10 5 years in duration. Motivated by the discovery of the spectacular nebula known as Hanny’s Voorwerp, ionized by a powerful AGN which has apparently faded dramatically within � 10 5 years, Galaxy Zoo volunteers have carried out both targeted and serendipitous searches for similar emission-line clouds around lowredshift galaxies. We present the resulting list of candidates and describe spectroscopy identifying 19 galaxies with AGN-ionized regions at projected radii rproj > 10 kpc. This search recovered known EELRs (such as Mkn 78, Mkn 266, and NGC 5252) and identified additional previously unknown cases, one with detected emission to r = 37 kpc. One new Sy 2 was identified. At least 14/19 are in interacting or merging systems, suggesting that tidal tails are a prime source of distant gas out of the galaxy plane to be ionized by an AGN. We see a mix of one- and two-sided structures, with observed cone angles from 23–112 ◦ . We consider the energy balance in the ionized clouds, with lower and upper bounds on ionizing luminosity from recombination and ionizationparameter arguments, and estimate the luminosity of the core from the far-infrared data. The implied ratio of ionizing radiation seen by the clouds to that emitted by the nucleus, on the assumption of a nonvariable nuclear source, ranges from 0.02 to > 12; 7/19 exceed unity. Small values fit well with a heavily obscured AGN in which only a small fraction of the ionizing output escapes to be traced by surrounding gas. However, large values may require that the AGN has faded over tens of thousands of years, giving us several examples of systems in which such dramatic long-period variation has occurred; this is the only current technique for addressing these timescales in AGN history. The relative numbers of faded and non-faded objects we infer, and the projected extents of the ionized regions, give our estimate (0.2–2×10 5 years ) for

GEOMETRIC AND DYNAMICAL MODELS OF REVERBERATION MAPPING DATA

We present a general method to analyze reverberation (or echo) mapping data that simultaneously provides estimates for the black hole mass and for the geometry and dynamics of the broad-line region (BLR) in active galactic nuclei (AGNs). While previous methods yield a typical scale size of the BLR or a reconstruction of the transfer function, our method directly infers the spatial and velocity distribution of the BLR from the data, from which a transfer function can be easily derived. Previous echo mapping analysis requires an independent estimate of a scaling factor known as the virial coefficient to infer the mass of the black hole, but this is not needed in our more direct approach. We use the formalism of Bayesian probability theory and implement a Markov Chain Monte Carlo algorithm to obtain estimates and uncertainties for the parameters of our BLR models. Fitting of models to the data requires knowledge of the continuum flux at all times, not just the measured times. We use Gaussian Processes to interpolate and extrapolate the continuum light curve data in a fully consistent probabilistic manner, taking the associated errors into account. We illustrate our method using simple models of BLR geometry and dynamics and show that we can recover the parameter values of our test systems with realistic uncertainties that depend upon the variability of the AGN and the quality of the reverberation mapping observing campaign. With a geometry model we can recover the mean radius of the BLR to within ~0.1 dex random uncertainty for simulated data with an integrated line flux uncertainty of 1.5%, while with a dynamical model we can recover the black hole mass and the mean radius to within ~0.05 dex random uncertainty, for simulated data with a line profile average signal-to-noise ratio of 4 per spectral pixel. These uncertainties do not include modeling errors, which are likely to be present in the analysis of real data, and should therefore be considered as lower limits to the accuracy of the method.

The lick agn monitoring project 2011: Fe II reverberation from the outer broad-line region

The prominent broad Fe II emission blends in the spectra of active galactic nuclei have been shown to vary in response to continuum variations, but past attempts to measure the reverberation lag time of the optical Fe II lines have met with only limited success. Here we report the detection of Fe II reverberation in two Seyfert 1 galaxies, NGC 4593 and Mrk 1511, based on data from a program carried out at Lick Observatory in Spring 2011. Light curves for emission lines including Hβ and Fe II were measured by applying a fitting routine to decompose the spectra into several continuum and emission-line components, and we use cross-correlation techniques to determine the reverberation lags of the emission lines relative to V-band light curves. In both cases, the measured lag (τcen) of Fe II is longer than that of Hβ, although the inferred lags are somewhat sensitive to the choice of Fe II template used in the fit. For spectral decompositions done using the Fe II template of Veron-Cetty et al., we find τcen (Fe II)/τcen (Hβ) = 1.9 ± 0.6 in NGC 4593 and 1.5 ± 0.3 in Mrk 1511. The detection of highly correlated variations between Fe II and continuum emission demonstrates that the Fe II emission in these galaxies originates in photoionized gas, located predominantly in the outer portion of the broad-line region.

The lick AGN monitoring project 2011: Dynamical modeling of the broad-line region in Mrk 50

We present dynamical modeling of the broad-line region (BLR) in the Seyfert 1 galaxy Mrk 50 using reverberation mapping data taken as part of the Lick AGN Monitoring Project (LAMP) 2011. We model the reverberation mapping data directly, constraining the geometry and kinematics of the BLR, as well as deriving a black hole mass estimate that does not depend on a normalizing factor or virial coefficient. We find that the geometry of the BLR in Mrk 50 is a nearly face-on thick disk, with a mean radius of 9.6^(+1.2)_(–0.9) light days, a width of the BLR of 6.9^(+1.2)_(–1.1) light days, and a disk opening angle of 25 ± 10 deg above the plane. We also constrain the inclination angle to be 9^(+7)_(–5) deg, close to face-on. Finally, the black hole mass of Mrk 50 is inferred to be log_(10)(M_(BH)/M_☉) = 7.57^(+0.44)_(–0.27). By comparison to the virial black hole mass estimate from traditional reverberation mapping analysis, we find the normalizing constant (virial coefficient) to be log_(10) f = 0.78^(+0.44)_(–0.27), consistent with the commonly adopted mean value of 0.74 based on aligning the M_(BH)-σ* relation for active galactic nuclei and quiescent galaxies. While our dynamical model includes the possibility of a net inflow or outflow in the BLR, we cannot distinguish between these two scenarios.

BROAD-LINE REVERBERATION IN THE KEPLER-FIELD SEYFERT GALAXY Zw 229-015

The Seyfert 1 galaxy Zw 229-015 is among the brightest active galaxies being monitored by the Kepler mission. In order to determine the black hole mass in Zw 229-015 from Hβ reverberation mapping, we have carried out nightly observations with the Kast Spectrograph at the Lick 3 m telescope during the dark runs from 2010 June through December, obtaining 54 spectroscopic observations in total. We have also obtained nightly V-band imaging with the Katzman Automatic Imaging Telescope at Lick Observatory and with the 0.9 m telescope at the Brigham Young University West Mountain Observatory over the same period. We detect strong variability in the source, which exhibited more than a factor of two change in broad Hβ flux. From cross-correlation measurements, we find that the Hβ light curve has a rest-frame lag of 3.86+0.69 –0.90 days with respect to the V-band continuum variations. We also measure reverberation lags for Hα and Hγ and find an upper limit to the Hδ lag. Combining the Hβ lag measurement with a broad Hβ width of σline = 1590 ± 47 km s–1 measured from the rms variability spectrum, we obtain a virial estimate of M BH = 1.00+0.19 –0.24 × 107 M ☉ for the black hole in Zw 229-015. As a Kepler target, Zw 229-015 will eventually have one of the highest-quality optical light curves ever measured for any active galaxy, and the black hole mass determined from reverberation mapping will serve as a benchmark for testing relationships between black hole mass and continuum variability characteristics in active galactic nuclei.

The Lick AGN Monitoring Project 2011: Reverberation Mapping of Markarian 50

The Lick AGN Monitoring Project 2011 observing campaign was carried out over the course of 11 weeks in spring 2011. Here we present the first results from this program, a measurement of the broad-line reverberation lag in the Seyfert 1 galaxy Mrk 50. Combining our data with supplemental observations obtained prior to the start of the main observing campaign, our data set covers a total duration of 4.5 months. During this time, Mrk 50 was highly variable, exhibiting a maximum variability amplitude of a factor of ~4 in the U-band continuum and a factor of ~2 in the Hβ line. Using standard cross-correlation techniques, we find that Hβ and Hγ lag the V-band continuum by τ_(cen) = 10.64^(+0.82)_(–0.93) and 8.43^(+1.30)_(–1.28) days, respectively, while the lag of He II λ4686 is unresolved. The Hβ line exhibits a symmetric velocity-resolved reverberation signature with shorter lags in the high-velocity wings than in the line core, consistent with an origin in a broad-line region (BLR) dominated by orbital motion rather than infall or outflow. Assuming a virial normalization factor of f = 5.25, the virial estimate of the black hole mass is (3.2 ± 0.5) × 10^7 M_☉. These observations demonstrate that Mrk 50 is among the most promising nearby active galaxies for detailed investigations of BLR structure and dynamics.

The Mass of the Black Hole in Arp 151 from Bayesian Modeling of Reverberation Mapping Data

Supermassive black holes are believed to be ubiquitous at the centers of galaxies. Measuring their masses is extremely challenging yet essential for understanding their role in the formation and evolution of cosmic structure. We present a direct measurement of the mass of a black hole in an active galactic nucleus (Arp 151) based on the motion of the gas responsible for the broad emission lines. By analyzing and modeling spectroscopic and photometric time series, we find that the gas is well described by a disk or torus with an average radius of 3.99 ± 1.25 light days and an opening angle of 68.9+21.4 – 17.2 deg, viewed at an inclination angle of 67.8 ± 7.8 deg (that is, closer to face-on than edge-on). The black hole mass is inferred to be 106.51 ± 0.28 M ☉. The method is fully general and can be used to determine the masses of black holes at arbitrary distances, enabling studies of their evolution over cosmic time.

SPACE TELESCOPE AND OPTICAL REVERBERATION MAPPING PROJECT. I. ULTRAVIOLET OBSERVATIONS OF THE SEYFERT 1 GALAXY NGC 5548 WITH THE COSMIC ORIGINS SPECTROGRAPH ON HUBBLE SPACE TELESCOPE

We describe the first results from a six-month long reverberation-mapping experiment in the ultraviolet based on 171 observations of the Seyfert 1 galaxy NGC 5548 with the Cosmic Origins Spectrograph on the Hubble Space Telescope. Significant correlated variability is found in the continuum and broad emission lines, with amplitudes ranging from ∼30% to a factor of two in the emission lines and a factor of three in the continuum. The variations of all the strong emission lines lag behind those of the continuum, with He II λ1640 lagging behind the continuum by ∼2.5 days and Lyα λ1215 ,C IV λ1550, and Si IV λ1400 lagging by ∼5–6 days. The relationship between the continuum and emission lines is complex. In particular, during the second half of the campaign, all emission-line lags increased by a factor of 1.3–2 and differences appear in the detailed structure of the continuum and emissionline light curves. Velocity-resolved cross-correlation analysis shows coherent structure in lag versus line of sight velocity for the emission lines; the high-velocity wings of C IV respond to continuum variations more rapidly than the line core, probably indicating higher velocity broad-line region clouds at smaller distances from the central