Reverberation Signatures Research Articles

Tracing the evolving X-ray reverberation lags within an individual AGN light curve

Abstract We present the Granger causality (GC) test for the X-ray reverberation analysis of Active Galactic Nuclei (AGN). If the light curves in the continuum-dominated band help predict (Granger cause) those dominated by reflection, the Granger lags that associate to the intrinsic reverberation lags can be inferred. We focus on six AGN observed by XMM-Newton, including the sources well-known to exhibit clear X-ray reverberation lags (IRAS 13224–3809 and 1H 0707–495) and those in which reverberation signatures are not well confirmed (MCG–6-30-15, IZW1, Mrk 704 and Mrk 1040). We employ the sliding-window algorithm and estimate the Granger (intrinsic) Fe-L lags along the light curve as the window moves through. This reveals the evolving lags towards the end of some individual observations, suggesting that the corona varies progressively. Occasionally, we observe two clearly separate lags that suggest an extended corona consisting of two zones while producing competing reverberation of two lags. While the GC test is purely hypothetical and might not explain true causality, our conclusion is that the lags are present and could be understood as reverberation lags. Assuming the lags changing solely with the corona, we find that the IRAS 13224–3809 corona varies between ∼10–25 rg and sometimes move to ≳ 50 rg. The corona of 1H 0707–495 and MCG–6-30-15 may be analogous to that of IRAS 13224–3809, while in IZw1, Mrk 704 and Mrk 1040 a more compact corona is expected.

X-ray timing and spectral analysis of reverberating active galactic nuclei

ABSTRACT We use the publicly available XMM–Newton archive to select a sample of 20 active galactic nuclei (AGNs) known to exhibit reverberation signatures caused by the reflection of X-rays from the corona off the accretion disc that feeds the central black hole engine. Inverse Compton scattering by energetic electrons, coupled with accretion disc fluctuations give rise to the highly variable observed X-ray spectrum, the behaviour of which is still not fully understood. We use 121 observations in 3–4 distinct spectral states for each source and calculate the time lags as a function of frequency. We fit the relativistic reflection model RELXILL and explore parameter correlations. The known scaling relationship between the black hole mass and time lag is well recovered and the continuum flux is coupled strongly to the disc reflection flux. We also find that 1H 0707−495 and IRAS 13224−3809 are well described using reflection and absorption modelling in a variety of flux states. The reflection fraction is strongly coupled to the power-law photon index and may be linked to dynamics of the emitting region. The data reveal hints of the power law evolutionary turnover when the 2–10 keV Eddington fraction is ∼0.02, the origin of which is not fully understood. Finally, we report the covering fraction is inversely correlated with the flux and power-law photon index in IRAS 13224−3809. These findings support recent studies of 1H 0707−495 where the covering fraction may contribute to the observed variability via flux modulations from non-uniform orbiting clouds.

Velocity-resolved Reverberation Mapping of Changing-look AGN NGC 2617

Abstract NGC 2617 attracted a lot of attention after the detection of changes in its spectral type; the geometry and kinematics of the broad-line region (BLR) are still ambiguous. In this paper, we present the high cadence (∼2 days) reverberation mapping campaign of NGC 2617 from 2019 October to 2020 May undertaken at the Lijiang 2.4 m telescope. For the first time, the velocity-resolved reverberation signature of the object was successfully detected. Both Hα and Hβ show an asymmetrical profile with a peak in the velocity-resolved time lags. For both of the lines, the lag of the line core is longer than those of the relevant wings, and the peak of the velocity-resolved lags is slightly blueshifted. These characteristics are not consistent with the theoretical prediction of the inflow, outflow or Keplerian disk model. Our observations give the time lags of Hα, Hβ, Hγ, and He i, with a ratio of τ Hα : τ Hβ : τ Hγ : τ He I = 1.27:1.00:0.89:0.20, which indicates a stratified structure in the BLR of the object. It is the first time that the lags of Hα and He i are obtained. Assuming a virial factor of f = 5.5 for the dispersion width of the line, the masses of the black holes derived from Hα and Hβ are and , respectively. Our observed results indicate the complexity of the BLR of NGC 2617.

Characterizing continuum variability in the radio-loud narrow-line Seyfert 1 galaxy IRAS 17020+4544

ABSTRACT We present results of temporal and spectral analyses on four XMM–Newton EPIC pn observations of IRAS 17020+4544, a narrow-line Seyfert 1 galaxy with evidence of a radio jet. Analysis of the light curves reveals that this radio-loud source does not behave like the bulk population of its radio-quiet counterparts. A trend of spectral hardening with increased flux is found. Variability is found to increase with energy, though it decreases as the spectrum hardens. The first 40 ks of the most recent observation behave uniquely among the epochs, exhibiting a softer spectral state than at any other time. Possible non-stationarity at low energies is found, with no such effect present at higher energies, suggesting at least two distinct spectral components. A reverberation signature is confirmed, with the lag-frequency, lag-energy, and covariance spectra changing significantly during the soft-state epoch. The temporal analysis suggests a variable power law in the presence of a reflection component, thus motivating such a fit for the 0.3−10 keV EPIC pn spectra from all epochs. We find an acceptable spectral fit using the timing-motivated parameters and report the detection of a broad Fe K emission line, requiring an additional model component beyond the reflection spectrum. We discuss links between this source and other narrow-line Seyfert 1 sources that show evidence of jet activity, finding similarities among this currently very limited sample of interesting objects.

Modeling the X-Ray Power Spectra of AGN by Using Reprocessing Echoes from an Extended Corona

Abstract Characteristic signatures that X-ray reverberation from an extended corona can manifest in the observed PSD of AGN are investigated. The presence of two X-ray blobs illuminating an accretion disk can cause the interference between two reprocessing-echo components and produce distinct physical features in the PSD. The oscillatory structures (e.g., dips and humps) are seen, but, contrarily to the lamp-post case, the strongest dip is not always the one at the lowest frequency. Instead, we find that the frequency where the strongest dip is seen is associated with the lower-source height while the lowest frequency where the first dip appears is linked to the upper-source height. This is because the reverberation timescales increase with the source height. Accurate modeling of the PSD then helps put constraints on the lower and upper limits of the corona extent. Furthermore, the reverberation signatures are less pronounced with increasing number of sources that do not produce reflection (e.g., additional X-rays from fast, relativistic outflows). The amplitude of the oscillations also depends on the amount of dilution contributed by the X-ray sources, and thus encodes information about their relative brightness. Due to stronger dilutions, robust detection of these signatures with current observations will become even more difficult if the corona is extended. Future observations made by Athena will enable us to fit these characteristics with statistically significant details, and to reveal the nature of the disk–corona system.

X-Ray Reverberation from Black Hole Accretion Disks with Realistic Geometric Thickness

Abstract X-ray reverberation in active galactic nuclei, believed to be the result of the reprocessing of corona photons by the underlying accretion disk, has allowed us to probe the properties of the innermost regions of the accretion flow and the central black hole. This process is modeled via raytracing in the Kerr metric, with the disk thickness almost ubiquitously assumed to be negligible (“razor thin”) and the corona commonly approximated as a point source located along the polar axis (a lamppost). In this work, we use the new raytracing suite, Fenrir, to explore the effect that accretion disk geometry has on reverberation signatures, assuming a lamppost configuration but allowing for a finite disk scale height. We characterize the signatures of finite disk thickness in the reverberation transfer function and calculate how they might manifest in observed lag-frequency spectra. We also show that a disk-hugging corona (approximated by off-axis point-like flares) exhibits characteristics that are qualitatively different from observation, thus providing further evidence for a flaring corona that is separated from the underlying disk material.

Velocity-resolved Reverberation Mapping of Five Bright Seyfert 1 Galaxies

Abstract We present the first results from a reverberation-mapping campaign undertaken during the first half of 2012, with additional data on one active galactic nucleus (AGN) (NGC 3227) from a 2014 campaign. Our main goals are (1) to determine the black hole masses from continuum-Hβ reverberation signatures, and (2) to look for velocity-dependent time delays that might be indicators of the gross kinematics of the broad-line region. We successfully measure Hβ time delays and black hole masses for five AGNs, four of which have previous reverberation mass measurements. The values measured here are in agreement with earlier estimates, though there is some intrinsic scatter beyond the formal measurement errors. We observe velocity-dependent Hβ lags in each case, and find that the patterns have changed in the intervening five years for three AGNs that were also observed in 2007.

Quantifying the diffuse continuum contribution of BLR Clouds to AGN Continuum Inter-band Delays

Disk reverberation mapping of a handful of nearby AGN suggest accretion disk sizes which are a factor few too large for their luminosities, apparently at odds with the standard model. Here, we investigate the likely contribution to the measured delay signature of diffuse continuum emission arising from broad line region gas. We start by constructing spherically symmetric pressure-law BLR models (i.e., $P(r)\propto r^{-s}$) that approximately reproduce the observed emission line fluxes of the strong UV--optical emission-lines in the best-studied source, NGC~5548. We then determine the contribution of the diffuse continuum to the measured continuum flux and inter-band delays, accounting for the observed variability behaviour of the ionizing nuclear continuum. Those pressure-law models that approximately reproduce the observed emission-line luminosities unavoidably produce substantial diffuse continuum emission. This causes a significant contamination of the disk reverberation signature (i.e., wavelength-dependent continuum delays). Qualitatively, the diffuse continuum delay signatures produced by our models resemble that observed for NGC~5548, including the deviation of the lag spectrum above that of a simple power-law in wavelength, short-ward of the Balmer and Paschen jumps. Furthermore, for reasonable estimates of the BLR covering fraction, the delay induced by diffuse continuum emission causes elevated inter-band delays over the entire UV--optical regime; for these pressure-law models, there are no `disk-dominated' wavelength intervals. Thus, the diffuse continuum contribution must be taken into account in order to correctly infer AGN accretion disk sizes based on inter-band continuum delays.

The reverberation signatures of rotating disc winds in active galactic nuclei

The broad emission lines (BELs) in active galactic nuclei (AGN) respond to ionizing continuum variations. The time and velocity dependence of their response depends on the structure of the broad-line region: its geometry, kinematics and ionization state. Here, we predict the reverberation signatures of BELs formed in rotating accretion disc winds. We use a Monte Carlo radiative transfer and ionization code to predict velocity-delay maps for representative high- (C$~IV$) and low-ionization (H$\alpha$) emission lines in both high- and moderate-luminosity AGN. Self-shielding, multiple scattering and the ionization structure of the outflows are all self-consistently taken into account, while small-scale structure in the outflow is modelled in the micro-clumping approximation. Our main findings are: (1) The velocity-delay maps of smooth/micro-clumped outflows often contain significant negative responses. (2)~The reverberation signatures of disc wind models tend to be rotation dominated and can even resemble the classic "red-leads-blue" inflow signature. (3) Traditional "blue-leads-red" outflow signatures can usually only be observed in the long-delay limit. (4) Our models predict lag-luminosity relationships similar to those inferred from observations, but systematically underpredict the observed centroid delays. (5) The ratio between "virial product" and black hole mass predicted by our models depends on viewing angle. Our results imply that considerable care needs to be taken in interpreting data obtained by observational reverberation mapping campaigns. In particular, basic signatures such as "red-leads-blue", "blue-leads-red" and "blue and red vary jointly" are not always reliable indicators of inflow, outflow or rotation. This may help to explain the perplexing diversity of such signatures seen in observational campaigns to date.

Reverberation Mapping of Optical Emission Lines in Five Active Galaxies

Abstract We present the first results from an optical reverberation mapping campaign executed in 2014 targeting the active galactic nuclei (AGNs) MCG+08-11-011, NGC 2617, NGC 4051, 3C 382, and Mrk 374. Our targets have diverse and interesting observational properties, including a “changing look” AGN and a broad-line radio galaxy. Based on continuum-Hβ lags, we measure black hole masses for all five targets. We also obtain Hγ and He ii λ4686 lags for all objects except 3C 382. The He ii λ4686 lags indicate radial stratification of the BLR, and the masses derived from different emission lines are in general agreement. The relative responsivities of these lines are also in qualitative agreement with photoionization models. These spectra have extremely high signal-to-noise ratios (100–300 per pixel) and there are excellent prospects for obtaining velocity-resolved reverberation signatures.

Towards modelling X-ray reverberation in AGN: piecing together the extended corona

Models of X-ray reverberation from extended coronae are developed from general relativistic ray tracing simulations. Reverberation lags between correlated variability in the directly observed continuum emission and that reflected from the accretion disc arise due to the additional light travel time between the corona and reflecting disc. X-ray reverberation is detected from an increasing sample of Seyfert galaxies and a number of common properties are observed, including a transition from the characteristic reverberation signature at high frequencies to a hard lag within the continuum component at low frequencies, as well a pronounced dip in the reverberation lag at 3keV. These features are not trivially explained by the reverberation of X-rays originating from simple point sources. We therefore model reverberation from coronae extended both over the surface of the disc and vertically. Causal propagation through its extent for both the simple case of constant velocity propagation and propagation linked to the viscous timescale in the underlying accretion disc is included as well as stochastic variability arising due to turbulence locally on the disc. We find that the observed features of X-ray reverberation in Seyfert galaxies can be explained if the long timescale variability is dominated by the viscous propagation of fluctuations through the corona. The corona extends radially at low height over the surface of the disc but with a bright central region in which fluctuations propagate up the black hole rotation axis driven by more rapid variability arising from the innermost regions of the accretion flow.

The closest look at 1H0707−495: X-ray reverberation lags with 1.3 Ms of data

Reverberation lags in AGN were first discovered in the NLS1 galaxy, 1H0707-495. We present a follow-up analysis using 1.3 Ms of data, which allows for the closest ever look at the reverberation signature of this remarkable source. We confirm previous findings of a hard lag of ~100 seconds at frequencies v ~ [0.5 - 4] e-4 Hz, and a soft lag of ~30 seconds at higher frequencies, v ~ [0.6 - 3] e-3 Hz. These two frequency domains clearly show different energy dependences in their lag spectra. We also find evidence for a signature from the broad Fe K line in the high frequency lag spectrum. We use Monte Carlo simulations to show how the lag and coherence measurements respond to the addition of Poisson noise and to dilution by other components. With our better understanding of these effects on the lag, we show that the lag-energy spectra can be modelled with a scenario in which low frequency hard lags are produced by a compact corona responding to accretion rate fluctuations propagating through an optically thick accretion disc, and the high frequency soft lags are produced by short light-travel delay associated with reflection of coronal power-law photons off the disc.

Ionized X-ray reverberation mapping: testing a light-bending model for 1H0707−495

The narrow-line Seyfert 1 galaxy 1H0707−495 has shown the first clear evidence of X-ray reverberation in which the soft, disc-reflection and Fe-L-dominated 0.3–1 keV band lags behind the hard, continuum-dominated 1–4 keV band by ≈ 30 s. A simple reflection model, however, cannot explain the different amplitudes of variability between the continuum and reflection components in time-resolved spectra. Here we apply a light-bending model (LBM) to 1H0707−495 to determine whether or not it can fully account for both the time-averaged spectrum and observed X-ray variability and time lags. A point-like source is assumed to move vertically from 3.5rg to 8 rg approximately along the rotation axis of the black hole. Photon paths are integrated via a numerical method for each source position to calculate the corresponding direct and ionized reflection spectra from the accretion disc. The source is assumed to move up and down in such a way as to reproduce the 1–4 keV light curve from the XMM–Newton observations. We find that, even though the reproduced time-averaged spectrum generally agrees with the observations, this LBM cannot account for the observed time lags between the continuum and reflection spectrum. Furthermore, timing properties such as coherence are sensitive to the ionization state of the inner disc. We conclude that a more complex model of the source geometry and intrinsic variability are required, and probably one in which fluctuations propagating through the disc are directly associated with the X-ray sources. Our ionized reverberation mapping model is completely general, allowing other scenarios to be investigated in the future. Reverberation signatures are beginning to constrain the true nature of the X-ray source of not only 1H0707−495, but also other 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 dataset 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-beta line. Using standard cross-correlation techniques, we find that H-beta and H-gamma lag the V-band continuum by tau_cen = 10.64(-0.93,+0.82) and 8.43(-1.28,+1.30) days, respectively, while the lag of He II 4686 is unresolved. The H-beta 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 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 solar masses. These observations demonstrate that Mrk 50 is among the most promising nearby active galaxies for detailed investigations of broad-line region structure and dynamics.

Discovery of universal outflow structures above and below the accretion disc plane in radio-quiet quasars

59 quasars in the background of the Magellanic Clouds had brightness records monitored by the MACHO project during the years 1992–99. Because the circumpolar fields of these quasars had no seasonal sampling defects, their observation produced data sets well suited to further careful analysis. Following a preliminary report wherein we showed the existence of reverberation in the data for one of the radio-quiet quasars in this group, we now show that similar reverberations have been seen in all of the 55 radio-quiet quasars with adequate data, making possible the determination of the quasar inclination to the observer’s line of sight. The reverberation signatures indicate the presence of large-scale elliptical outflow structures similar to that predicted by the Elvis and ‘dusty torus’ models of quasars, whose characteristic sizes vary within a surprisingly narrow range of scales. More importantly, the observed opening angle relative to the polar axis of the universal elliptical outflow structure present was consistently found to be on the order of 78°.

Iron Line Reverberation Mapping withConstellation‐X

The broad X-ray iron line seen in the spectra of many AGN is thought to originate from the inner regions of the putative black hole accretion disk, and hence provides a rare probe of that central region. In principle, future high throughput X-ray spectrometers should be able to examine the temporal response of this fluorescent line to flares in the X-ray continuum source (which energizes this emission line) - i.e. iron line reverberation mapping will be possible. It has been previously found that there are robust reverberation signatures of the black hole spin, mass and the X-ray flare location above the accretion disk. Here, we simulate observations of a bright Seyfert nucleus with the proposed NASA mission Constellation-X and demonstrate the feasibility of detecting these reverberation signatures with this mission. Hence, starting with XMM in 1999, and maturing with Constellation-X in c2010, iron line reverberation will open a new window on the innermost regions of AGN.

On the Lack of X‐Ray Iron Line Reverberation in MCG −6‐30‐15: Implications for the Black Hole Mass and Accretion Disk Structure

We use the method of Press, Rybicki & Hewitt (1992) to search for time lags and time leads between different energy bands of the RXTE data for MCG-6-30-15. We tailor our search in order to probe any reverberation signatures of the fluorescent iron Kalpha line that is thought to arise from the inner regions of the black hole accretion disk. In essence, an optimal reconstruction algorithm is applied to the continuum band (2-4keV) light curve which smoothes out noise and interpolates across the data gaps. The reconstructed continuum band light curve can then be folded through trial transfer functions in an attempt to find lags or leads between the continuum band and the iron line band (5-7keV). We find reduced fractional variability in the line band. The spectral analysis of Lee et al. (1999) reveals this to be due to a combination of an apparently constant iron line flux (at least on timescales of few x 10^4s), and flux correlated changes in the photon index. We also find no evidence for iron line reverberation and exclude reverberation delays in the range 0.5-50ksec. This extends the conclusions of Lee et al. and suggests that the iron line flux remains constant on timescales as short as 0.5ksec. The large black hole mass (>10^8Msun) naively suggested by the constancy of the iron line flux is rejected on other grounds. We suggest that the black hole in MCG-6-30-15 has a mass of M_BH~10^6-10^7Msun and that changes in the ionization state of the disk may produce the puzzling spectral variability. Finally, it is found that the 8-15keV band lags the 2-4keV band by 50-100s. This result is used to place constraints on the size and geometry of the Comptonizing medium responsible for the hard X-ray power-law in this AGN.

X-ray iron line variability for the model of an orbiting flare above a black hole accretion disc

The broad X-ray iron line, detected in many active galactic nuclei, is likely to be produced by fluorescence from the X-ray-illuminated central parts of an accretion disc close to a supermassive black hole. The time-averaged shape of the line can be explained most naturally by a combination of special and general relativistic effects. Such line profiles contain information about the black hole spin and the accretion disc, as well as the geometry of the emitting region, and may help to test general relativity in the strong gravity regime. In this paper we embark on the computation of the temporal response of the line to the illuminating flux. Previous studies concentrated on the calculation of reverberation signatures from static sources illuminating the disc. In this paper we focus on the more physically justified case of flares located above the accretion disc and corotating with it. We compute the time-dependent iron line, taking into account all general relativistic effects, and show that its shape is of a very complex nature, and we also present light curves accompanying the iron line variability. We suggest that present and future X-ray satellites like XMM or Constellation-X may be capable of detecting features present in the computed reverberation maps.

X‐Ray Iron Line Reverberation from Black Hole Accretion Disks

The relativistically broad X-ray iron line seen in many AGN spectra is thought to originate from the central regions of the putative black hole accretion disk. Both the line profile and strength will vary in response to rapid variability of the primary X-ray continuum source. The temporal response of the line contains information on the accretion disk structure, the X-ray source geometry, and the spin of the black hole. Since the X-ray source will have a size comparable to the fluorescing region of the accretion disk, the general reverberation problem is not invertible. However, progress can be made since, empirically, AGN light curves are seen to undergo dramatic short timescale variability which presumably corresponds to the creation of a single new active region within the distributed X-ray source. The iron line response to these individual events can be described using linear transfer theory. We consider the line response to the activation/flaring of a new X-ray emitting region. Most of our detailed calculations are performed for the case of an X-ray source on the symmetry axis and at some height above the disk plane around a Kerr black hole. We also present preliminary calculations for off-axis flares. We suggest ways in which future, high-throughput X-ray observatories such as XMM and the Constellation X-ray Mission may use these reverberation signatures to probe both the mass and spin of AGN black holes, as well as the X-ray source geometry.