X-ray Reverberation Research Articles

Multi-time-scale X-ray reverberation mapping of accreting black holes

Accreting black holes show characteristic reflection features in their X-ray spectrum, including an iron K$\alpha$ line, resulting from hard X-ray continuum photons illuminating the accretion disk. The reverberation lag resulting from the path length difference between direct and reflected emission provides a powerful tool to probe the innermost regions around both stellar-mass and supermassive black holes. Here, we present for the first time a reverberation mapping formalism that enables modeling of energy dependent time lags and variability amplitude for a wide range of variability timescales, taking the complete information of the cross-spectrum into account. We use a pivoting power-law model to account for the spectral variability of the continuum that dominates over the reverberation lags for longer time scale variability. We use an analytic approximation to self-consistently account for the non-linear effects caused by this continuum spectral variability, which have been ignored by all previous reverberation studies. We find that ignoring these non-linear effects can bias measurements of the reverberation lags, particularly at low frequencies. Since our model is analytic, we are able to fit simultaneously for a wide range of Fourier frequencies without prohibitive computational expense. We also introduce a formalism of fitting to real and imaginary parts of our cross-spectrum statistic, which naturally avoids some mistakes/inaccuracies previously common in the literature. We perform proof-of-principle fits to Rossi X-ray Timing Explorer data of Cygnus X-1.

X-ray lags in PDS 456 revealed by Suzaku observations

X-ray reverberation lags from the vicinity of supermassive black holes have been detected in almost 30 AGN. The soft lag, which is the time delay between the hard and soft X-ray light curves, is usually interpreted as the time difference between the direct and reflected emission, but is alternatively suggested to arise from the direct and scattering emission from distant clouds. By analysing the archival Suzaku observations totalling an exposure time of ~ 770 ks, we discover a soft lag of $10\pm3.4$ ks at $9.58\times10^{-6}$ Hz in the luminous quasar PDS 456, which is the longest soft lag and lowest Fourier frequency reported to date. In this study, we use the maximum likelihood method to deal with non-continuous nature of the Suzaku light curves. The result follows the mass-scaling relation for soft lags, which further supports that soft lags originate from the innermost areas of AGN and hence are best interpreted by the reflection scenario. Spectral analysis has been performed in this work and we find no evidence of clumpy partial-covering absorbers. The spectrum can be explained by a self-consistent relativistic reflection model with warm absorbers, and spectral variations over epochs can be accounted for by the change of the continuum, and of column density and ionization states of the warm absorbers.

Evolution of the reverberation lag in GX 339–4 at the end of an outburst

We studied X-ray reverberation lags in the BHXRB GX 339-4 at the end of the 2014-2015 outburst. We analysed data from a XMM-Newton campaign covering the end of the transition from the soft to the hard state, and the decrease of luminosity in the hard state. During all the observations we detected, at high frequencies, significant disc variability, responding to variations of the power law emission with an average time delay of $\sim0.009\pm 0.002$ s. These new detections of disc thermal reverberation add to those previously obtained and suggest the lag to be always present in hard and hard-intermediate states. Our study reveals a net decrease of lag amplitude as a function of luminosity. We ascribe this trend to variations of the inner flow geometry. A possible scenario implies a \emph{decrease} of the inner disc truncation radius as the luminosity increases at the beginning of the outburst, followed by an \emph{increase} of the inner disc truncation radius as the luminosity decreases at the end of the outburst. Finally, we found hints of FeK reverberation ($\sim3\sigma$ significance) during the best quality observation of the XMM monitoring. The lag at the FeK energy has similar amplitude as that of the thermally reprocessed component, as expected if the same irradiated region of the disc is responsible for producing both the thermalized and reflected component. This finding suggests FeK reverberation in BHXRBs to be at the reach of current detectors provided observations of sufficiently long exposure are available.

ARE THE kHz QPO LAGS IN NEUTRON STAR 4U 1608–52 DUE TO REVERBERATION?

X-ray reverberation lags have recently been discovered in both active galactic nuclei (AGN) and black hole X-ray binaries. A recent study of the neutron star low-mass X-ray binary 4U 1608-52 has also shown significant lags, whose properties hint at a reverberation origin. Here, we adapt general relativistic ray tracing impulse response functions used to model X-ray reverberation in AGN for neutron star low-mass X-ray binaries. Assuming relativistic reflection forms the broad iron line and associated reflection continuum, we use reflection fits to the energy spectrum along with the impulse response functions to calculate the expected lags as a function of energy over the range of observed kHz QPO frequencies in 4U 1608-52. We find that the lag energy spectrum is expected to increase with increasing energy above 8 keV, while the observed lags in 4U 1608-52 show the opposite behavior. This demonstrates that the lags in the lower kHz QPO of 4U 1608-52 are not solely due to reverberation. We do note, however, that the models appear to be more consistent with the much flatter lag energy spectrum observed in the upper kHz QPO of several neutron star low-mass X-ray binaries, suggesting that lower and upper kHz QPOs may have different origins.

THE REVERBERATION LAG IN THE LOW-MASS X-RAY BINARY H1743-322

ABSTRACT The evolution of the inner accretion flow of a black hole X-ray binary during an outburst is still a matter of active research. X-ray reverberation lags are powerful tools for constraining disk–corona geometry. We present a study of X-ray lags in the black hole transient H1743-322. We compared the results obtained from analysis of all the publicly available XMM-Newton observations. These observations were carried out during two different outbursts that occurred in 2008 and 2014. During all the observations the source was caught in the hard state and at similar luminosities ( ). We detected a soft X-ray lag of ∼60 ms, most likely due to thermal reverberation. We did not detect any significant change of the lag amplitude among the different observations, indicating a similar disk–corona geometry at the same luminosity in the hard state. On the other hand, we observe significant differences between the reverberation lag detected in H1743-322 and in GX 339-4 (at similar luminosities in the hard state), which might indicate variations of the geometry from source to source.

A global look at X-ray time lags in Seyfert galaxies

X-ray reverberation, where light-travel time delays map out the compact geometry around the inner accretion flow in supermassive black holes, has been discovered in several of the brightest, most variable and well-known Seyfert galaxies. In this work, we expand the study of X-ray reverberation to all Seyfert galaxies in the XMM-Newton archive above a nominal rms variability and exposure level (a total of 43 sources). 50 per cent of source exhibit iron K reverberation, in that the broad iron K emission line responds to rapid variability in the continuum. We also find that on long timescales, the hard band emission lags behind the soft band emission in 85 per cent of sources. This `low-frequency hard lag' is likely associated with the coronal emission, and so this result suggests that most sources with X-ray variability show intrinsic variability from the nuclear region. We update the known iron K lag amplitude vs. black hole mass relation, and find evidence that the height or extent of the coronal source (as inferred by the reverberation time delay) increases with mass accretion rate.

Relativistic X-ray reverberation modelling of the combined time-averaged and lag-energy spectra in AGN

General relativistic ray tracing simulations of the time-averaged spectrum and energy-dependent time delays in AGN are presented. We model the lamp-post geometry in which the accreting gas is illuminated by an X-ray source located on the rotation axis of the black hole. The spectroscopic features imprinted in the reflection component are modelled using REFLIONX. The associated time delays after the direct continuum, known as reverberation lags, are computed including the full effects of dilution and ionization gradients on the disc. We perform, for the first time, simultaneous fitting of the time-averaged and lag-energy spectra in three AGN: Mrk 335, IRAS 13224-3809 and Ark 564 observed with XMM-Newton. The best fitting source height and central mass of each AGN partly agree with those previously reported. We find that including the ionization gradient in the model naturally explains lag-energy observations in which the 3 keV and 7-10 keV bands precede other bands. To obtain the clear 3 keV and 7-10 keV dips in the lag-energy profile, the model requires either a source height > 5$r_g$, or a disc that is highly ionized at small radii and is colder further out. We also show that fitting the lag or the mean spectra alone can lead to different results and interpretations. This is therefore important to combine the spectral and timing data in order to find the plausible but self-consistent fits which is achievable with our model.

Extensive X-ray variability studies of NGC 7314 using longXMM–Newtonobservations

We present a detailed X-ray variability study of the low mass Active Galactic Nuclei (AGN) NGC 7314 using the two newly obtained XMM-Newton observations ($140$ and $130$ ks), together with two archival data sets of shorter duration ($45$ and $84$ ks). The relationship between the X-ray variability characteristics and other physical source properties (such as the black hole mass) are still relatively poorly defined, especially for low-mass AGN. We perform a new, fully analytical, power spectral density (PSD) model analysis method, which will be described in detail in a forthcoming paper, that takes into consideration the spectral distortions, caused by red-noise leak. We find that the PSD in the $0.5-10$ keV energy range, can be represented by a bending power-law with a bend around $6.7\times10^{-5}$ Hz, having a slope of $0.51$ and $1.99$ below and above the bend, respectively. Adding our bend time-scale estimate, to an already published ensemble of estimates from several AGN, supports the idea that the bend time-scale depends linearly only on the black hole mass and not on the bolometric luminosity. Moreover, we find that as the energy range increases, the PSD normalization increases and there is a hint that simultaneously the high frequency slope becomes steeper. Finally, the X-ray time-lag spectrum of NGC 7314 shows some very weak signatures of relativistic reflection, and the energy resolved time-lag spectrum, for frequencies around $3\times10^{-4}$ Hz, shows no signatures of X-ray reverberation. We show that the previous claim about ks time-delays in this source, is simply an artefact induced by the minuscule number of points entering during the time-lag estimation in the low frequency part of the time-lag spectrum (i.e. below $10^{-4}$ Hz).

Signatures of X-ray reverberation in the power spectra of AGN

We compute fully relativistic disc response functions in the case of the "lamp-post" geometry using the full observed reflection spectrum for various X-ray source heights, disc inclination, and spin values of the central black hole. Since the observed PSD is equal to the product of the intrinsic power spectrum with the "transfer function" (i.e. the Fourier transform of the disc response function), we are able to predict the observed PSDs in the case of X-ray illumination of the inner disc. The observed PSD should show a prominent dip at high frequencies and an oscillatory behaviour, with a decreasing amplitude, at higher frequencies. The reverberation "echo" features should be more prominent in energy bands where the reflection component is more pronounced. The frequency of the dip is independent of energy, and it is mainly determined by the black hole mass and the X-ray source height. The amplitude of the dip increases with increasing black hole spin and inclination angle, as long as the height of the "lamp" is smaller than ~10 gravitational radii. The detection of the X-ray reverberation signals in the PSDs can provide further evidence for X-ray illumination of the inner disc in AGN. Our results are largely independent of the assumed geometry of the disc-corona system, as long as it does not change with time, and the disc response function is characterized by a sharp rise, a "plateau", and a decline at longer times. Irrespective of the geometry, the frequency of the main dip should decrease with increasing "mean time" of the response function, and the amplitude of the dip should increase with increasing reflection fraction.

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.

Strong constraints on a super-Eddington accretion flow:XMM–Newtonobservations of an intermediate-mass black hole

RX J1140.1+0307 (hereafter RX1140) is a Narrow Line Seyfert 1 (NLS1) with one of the lowest black hole masses known in an AGN (M ${\le} 10^6$ M$_{\odot}$). We show results from two new {\it XMM-Newton} observations, showing soft 2-10~keV spectra, a strong excess at lower energies, and fast X-ray variability as is typical of this class. The soft excess can be equally well fit by either low temperature Comptonisation or highly smeared, ionised reflection models, but we use a covariance analysis of the fast X-ray variability as well as lag and coherence spectra to show that the low temperature Comptonisation model gives a better description of the break in variability properties between soft and hard X-rays. Both models also require an additional component at the softest energies, as expected from the accretion disc. However, this inner disc spectrum does not join smoothly onto the variable optical and far UV emission (which should be produced in the outer disc) unless the mass is underestimated by an order of magnitude. The variable optical and far UV emission instead suggests that $L/L_{Edd}\sim 10$ through the outer disc, in which case advection and/or wind losses are required to explain the observed broadband spectral energy distribution. However, the similarity of the X-ray properties of RX1140 to other simple NLS1 such as PG 1244+026, RE J1034+396 and RX~J0136 means it is likely that these are also super-Eddington sources. This means their spectral energy distribution cannot be used to determine black hole spin despite appearing to be disc dominated. It also means that the accretion geometry close to the black hole is unlikely to be a flat disc as assumed in the new X-ray reverberation mapping techniques.

Iron K and Compton hump reverberation in SWIFT J2127.4+5654 and NGC 1365 revealed by NuSTAR and XMM–Newton

In the past five years, a flurry of X-ray reverberation lag measurements of accreting supermassive black holes have been made using the XMM-Newton telescope in the 0.3-10 keV energy range. In this work, we use the NuSTAR telescope to extend the lag analysis up to higher energies for two Seyfert galaxies, SWIFT J2127.4+5654 and NGC 1365. X-ray reverberation lags are due to the light travel time delays between the direct continuum emission and the reprocessed emission from the inner radii of an ionised accretion disc. XMM-Newton has been particularly adept at measuring the lag associated with the broad Fe K emission line, where the gravitationally redshifted wing of the line is observed to respond before the line centroid at 6.4 keV, produced at larger radii. Now we use NuSTAR to probe the lag at higher energies, where the spectrum shows clear evidence for Compton reflection, known as the Compton 'hump'. The XMM-Newton data show Fe K lags in both SWIFT J2127.4+5654 and NGC 1365. The NuSTAR data provide independent confirmation of these Fe K lags, and also show evidence for the corresponding Compton hump lags, especially in SWIFT J2127.4+5654. These broadband lag measurements confirm that the Compton hump and Fe K lag are produced at small radii. At low-frequencies in NGC 1365, where the spectrum shows evidence for eclipsing clouds in the line of sight, we find a clear negative (not positive) lag from 2--10 keV, which can be understood as the decrease in column density from a neutral eclipsing cloud moving out of our line of sight during the observation.

Pinpointing the base of the AGN jets through general relativistic X-ray reverberation studies

AbstractMany theoretical models of Active Galactic Nuclei (AGN) predict that the X-ray corona, lying above the black hole, constitutes the base of the X-ray jet. Thus, by studying the exact geometry of the close black hole environment, we can pinpoint the launching site of the jet. Detection of negative X-ray reverberation time delays (i.e. soft band X-ray variations lagging behind the corresponding hard band X-ray variations) can yield significant information about the geometrical properties of the AGN, such as the location of the X-ray source, as well as the physical properties of the the black hole, such as its mass and spin. In the frame-work of the lamp-post geometry, I present the first systematic X-ray time-lag modelling results of an ensemble of 12 AGN, using a fully general relativistic (GR) ray tracing approach for the estimation of the systems' response functions. By combing these state-of-the art GR response models with statistically innovative fitting routines, I derive the geometrical layout of the close BH environment for each source, unveiling the position of the AGN jet-base.

X-ray reverberation around accreting black holes

Luminous accreting stellar mass and supermassive black holes produce power-law continuum X-ray emission from a compact central corona. Reverberation time lags occur due to light travel time-delays between changes in the direct coronal emission and corresponding variations in its reflection from the accretion flow. Reverberation is detectable using light curves made in different X-ray energy bands, since the direct and reflected components have different spectral shapes. Larger, lower frequency, lags are also seen and are identified with propagation of fluctuations through the accretion flow and associated corona. We review the evidence for X-ray reverberation in active galactic nuclei and black hole X-ray binaries, showing how it can be best measured and how it may be modelled. The timescales and energy-dependence of the high frequency reverberation lags show that much of the signal is originating from very close to the black hole in some objects, within a few gravitational radii of the event horizon. We consider how these signals can be studied in the future to carry out X-ray reverberation mapping of the regions closest to black holes.

OBSERVATIONS OF MCG-5-23-16 WITHSUZAKU,XMM-NEWTONANDNUSTAR: DISK TOMOGRAPHY AND COMPTON HUMP REVERBERATION

MCG-5-23-16 is one of the first AGN where relativistic reverberation in the iron K line originating in the vicinity of the supermassive black hole was found, based on a short XMM-Newton observation. In this work, we present the results from long X-ray observations using Suzaku, XMM-Newton and NuSTAR designed to map the emission region using X-ray reverberation. A relativistic iron line is detected in the lag spectra on three different time-scales, allowing the emission from different regions around the black hole to be separated. Using NuSTAR coverage of energies above 10 keV reveals a lag between these energies and the primary continuum, which is detected for the first time in an AGN. This lag is a result of the Compton reflection hump responding to changes in the primary source in a manner similar to the response of the relativistic iron K line.

Modelling the broad Fe Kα reverberation in the AGN NGC 4151

The recent detection of X-ray reverberation lags, especially in the Fe Kα line region, around active galactic nuclei (AGN) has opened up the possibility of studying the time-resolved response (reflection) of hard X-rays from the accretion disc around supermassive black holes. Here, we use general relativistic transfer functions for reflection of X-rays from a point source located at some height above the black hole to study the time lags expected as a function of frequency and energy in the Fe Kα line region. We explore the models and the dependence of the lags on key parameters such as the height of the X-ray source, accretion disc inclination, black hole spin and black hole mass. We then compare these models with the observed frequency- and energy-dependence of the Fe Kα line lag in NGC 4151. Assuming the optical reverberation mapping mass of 4.6 × 107 M⊙, we get a best fit to the lag profile across the Fe Kα line in the frequency range (1-2) × 10− 5 Hz for an X-ray source located at a height |$h = 7^{+2.9}_{-2.6} R_G$| with a maximally spinning black hole and an inclination i < 30°.

The spin of supermassive black holes

Black hole spin is a quantity of great interest to both physicists and astrophysicists. We review the current status of spin measurements in supermassive black holes (SMBH). To date, every robust SMBH spin measurement uses x-ray reflection spectroscopy, so we focus almost exclusively on this technique as applied to moderately-luminous active galactic nuclei (AGN). After describing the foundations and uncertainties of the method, we summarize the current status of the field. At the time of writing, observations by XMM-Newton, Suzaku and NuSTAR have given robust spin constraints on 22 SMBHs. We find a significant number of rapidly rotating SMBHs (with dimensionless spin parameters a > 0.9) although, especially at the higher masses (M > 4 × 107M⊙), there are also some SMBHs with intermediate spin parameters. This may be giving us our first hint at a mass-dependent spin distribution which would, in turn, provide interesting constraints on models for SMBH growth. We also discuss the recent discovery of relativistic x-ray reverberation which we can use to ‘echo map’ the innermost regions of the accretion disc. The ultimate development of these reverberation techniques, when applied to data from future high-throughput x-ray observatories such as LOFT, ATHENA+, and AXSIO, will permit the measurement of black hole spin by a characterization of strong-field Shapiro delays. We conclude with a brief discussion of other electromagnetic methods that have been attempted or are being developed to constrain SMBH spin.

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.

DIRECT MEASUREMENT OF THE X-RAY TIME-DELAY TRANSFER FUNCTION IN ACTIVE GALACTIC NUCLEI

The origin of the observed time lags, in nearby active galactic nuclei (AGN), between hard and soft X-ray photons is investigated using new XMM-Newton data for the narrow-line Seyfert I galaxy Ark 564 and existing data for 1H0707-495 and NGC 4051. These AGN have highly variable X-ray light curves that contain frequent, high peaks of emission. The averaged light curve of the peaks is directly measured from the time series, and it is shown that (i) peaks occur at the same time, within the measurement uncertainties, at all X-ray energies, and (ii) there exists a substantial tail of excess emission at hard X-ray energies, which is delayed with respect to the time of the main peak, and is particularly prominent in Ark 564. Observation (i) rules out that the observed lags are caused by Comptonization time delays and disfavors a simple model of propagating fluctuations on the accretion disk. Observation (ii) is consistent with time lags caused by Compton-scattering reverberation from material a few thousand light-seconds from the primary X-ray source. The power spectral density and the frequency-dependent phase lags of the peak light curves are consistent with those of the full time series. There is evidence for non-stationarity in the Ark 564 time series in both the Fourier and peaks analyses. A sharp `negative' lag (variations at hard photon energies lead soft photon energies) observed in Ark 564 appears to be generated by the shape of the hard-band transfer function and does not arise from soft-band reflection of X-rays. These results reinforce the evidence for the existence of X-ray reverberation in type I AGN, which requires that these AGN are significantly affected by scattering from circumnuclear material a few tens or hundreds of gravitational radii in extent.

Relativistic iron K X-ray reverberation in NGC 4151

Recent X-ray observations have enabled the study of reverberation delays in AGN for the first time. All the detections so far are in sources with a strong soft excess, and the measured delay is between the hard (1-3 keV) direct continuum and the soft excess (0.5-1 keV), interpreted as the reflection continuum smeared by relativistic effects. There is however an inherent ambiguity in identifying and studying the details of the lines in the soft excess. Here we report the first detection of reverberation in the iron K band in any AGN. Using XMM-Newton observations of NGC 4151, we find delays of order 2000 s on time-scales of 10e5 s between the 5-6 keV band and 2-3 and 7-8 keV bands, with a broad lag profile resembling a relativistically-broadened iron line. The peak of the lag spectra shifts to lower energies at higher frequencies, consistent with the red wing of the line being emitted at smaller radii, as expected from reflection off the inner accretion disk. This is a first detection of a broad iron line using timing studies.