Lags measured from correlated X-ray/UV/optical monitoring of AGN allow us to determine whether UV/optical variability is driven by reprocessing of X-rays or X-ray variability is driven by UV/optical ...seed photon variations. We present the results of the largest study to date of the relationship between the X-ray, UV and optical variability in an AGN with 554 observations, over a 750 d period, of the Seyfert 1 galaxy NGC 5548 with Swift. There is a good overall correlation between the X-ray and UV/optical bands, particularly on short time-scales (tens of days). The UV/optical bands lag the X-ray band with lags which are proportional to wavelength raised to the power 1.23 ± 0.31. This power is very close to the power (4/3) expected if short time-scale UV/optical variability is driven by reprocessing of X-rays by a surrounding accretion disc. The observed lags, however, are longer than expected from a standard Shakura–Sunyaev accretion disc with X-ray heating, given the currently accepted black hole mass and accretion rate values, but can be explained with a slightly larger mass and accretion rate, and a generally hotter disc. Some long-term UV/optical variations are not paralleled exactly in the X-rays, suggesting an additional component to the UV/optical variability arising perhaps from accretion rate perturbations propagating inwards through the disc.
The production of artificial light curves with known statistical and variability properties is of great importance in astrophysics. Consolidating the confidence levels during cross-correlation ...studies, understanding the artefacts induced by sampling irregularities, establishing detection limits for future observatories are just some of the applications of simulated data sets. Currently, the widely used methodology of amplitude and phase randomization is able to produce artificial light curves which have a given underlying power spectral density (PSD) but which are strictly Gaussian distributed. This restriction is a significant limitation, since the majority of the light curves, e.g. active galactic nuclei, X-ray binaries, gamma-ray bursts, show strong deviations from Gaussianity exhibiting 'burst-like' events in their light curves yielding long-tailed probability density functions (PDFs). In this study, we propose a simple method which is able to precisely reproduce light curves which match both the PSD and the PDF of either an observed light curve or a theoretical model. The PDF can be representative of either the parent distribution or the actual distribution of the observed data, depending on the study to be conducted for a given source. The final artificial light curves contain all of the statistical and variability properties of the observed source or theoretical model, i.e. the same PDF and PSD, respectively. Within the framework of Reproducible Research, the code and the illustrative example used in this paper are both made publicly available in the form of an interactive mathematica notebook.
Swift monitoring of NGC 4151 with an ∼6 hr sampling over a total of 69 days in early 2016 is used to construct light curves covering five bands in the X-rays (0.3-50 keV) and six in the ultraviolet ...(UV)/optical (1900-5500 Å). The three hardest X-ray bands (>2.5 keV) are all strongly correlated with no measurable interband lag, while the two softer bands show lower variability and weaker correlations. The UV/optical bands are significantly correlated with the X-rays, lagging ∼3-4 days behind the hard X-rays. The variability within the UV/optical bands is also strongly correlated, with the UV appearing to lead the optical by ∼0.5-1 days. This combination of 3 day lags between the X-rays and UV and 1 day lags within the UV/optical appears to rule out the "lamp-post" reprocessing model in which a hot, X-ray emitting corona directly illuminates the accretion disk, which then reprocesses the energy in the UV/optical. Instead, these results appear consistent with the Gardner & Done picture in which two separate reprocessings occur: first, emission from the corona illuminates an extreme-UV-emitting toroidal component that shields the disk from the corona; this then heats the extreme-UV component, which illuminates the disk and drives its variability.
The extensive use of the structure function (SF) in the field of blazar variability across the electromagnetic spectrum suggests that characteristics time-scales are embedded in the light curves of ...these objects. We argue that for blazar variability studies, the SF results are sometimes erroneously interpreted leading to misconceptions about the actual source properties. Based on extensive simulations, we caution that spurious breaks will appear in the SFs of almost all light curves, even though these light curves may contain no intrinsic characteristic time-scales, i.e. having a featureless underlying power spectral density (PSD). We show that the time-scales of the spurious SF breaks depend mainly on the length of the artificial data set and also on the character of the variability i.e. the shape of the PSD. The SF is often invoked in the framework of shot-noise models to determine the temporal properties of individual shots. We caution that although the SF may be fitted to infer the shot parameters, the resultant shot-noise model is usually inconsistent with the observed PSD. As any model should fit the data in both the time and the frequency domain, the shot-noise model, in these particular cases, cannot be valid. Moreover, we show that the lack of statistical independence between adjacent SF points, in the standard SF formulation, means that it is not possible to perform robust statistical model fitting following the commonly used least-squares fitting methodology. The latter yields uncertainties in the fitting parameters (i.e. slopes, breaks) that are far too small with respect to their true statistical scatter. Finally, it is also commonly thought that SFs are immune to the sampling problems, such as data gaps, which affects the estimators of the PSDs. However, we show that SFs are also troubled by gaps which can induce artefacts.
We present the first systematic physical modelling of the time-lag spectra between the soft (0.3-1 keV) and the hard (1.5-4 keV) X-ray energy bands, as a function of Fourier frequency, in a sample of ...12 active galactic nuclei which have been observed by XMM-Newton. We concentrate particularly on the negative X-ray time-lags (typically seen above 10−4 Hz), i.e. soft-band variations lag the hard-band variations, and we assume that they are produced by reprocessing and reflection by the accretion disc within a lamp-post X-ray source geometry. We also assume that the response of the accretion disc, in the soft X-ray bands, is adequately described by the response in the neutral Fe Kα line at 6.4 keV for which we use fully general relativistic ray-tracing simulations to determine its time evolution. These response functions, and thus the corresponding time-lag spectra, yield much more realistic results than the commonly used, but erroneous, top-hat models. Additionally, we parametrize the positive part of the time-lag spectra (typically seen below 10−4 Hz) by a power law. We find that the best-fitting black hole (BH) masses, M, agree quite well with those derived by other methods, thus providing us with a new tool for BH mass determination. We find no evidence for any correlation between M and the BH spin parameter, α, the viewing angle, θ, or the height of the X-ray source above the disc, h. Also on average, the X-ray source lies only around 3.7 gravitational radii above the accretion disc and θ is distributed uniformly between 20° and 60°. Finally, there is a tentative indication that the distribution of α may be bimodal above and below 0.62.
We present X-ray spectral variability of 24 local active galactic nuclei (AGN) from the Palomar sample of nearby galaxies, as observed mainly by Swift. From hardness ratio measurements, we find that ...18 AGN with low accretion rates show hardening with increasing count rate, converse to the softer-when-brighter behaviour normally observed in AGN with higher accretion rates. Two AGN show softening with increasing count rate, two show more complex behaviour, and two do not show any simple relationship. Sufficient data were available for the spectra of 13 AGN to be summed in flux-bins. In nine of these sources, correlated luminosity-dependent changes in the photon index (Γ) of a power-law component are found to be the main cause of hardness variability. For six objects, with a low accretion rate as a fraction of the Eddington rate (
$\dot{m}_\mathrm{{Edd}}$
), Γ is anticorrelated with
$\dot{m}_\mathrm{{Edd}}$
, i.e. ‘harder-when-brighter’ behaviour is observed. The three higher
$\dot{m}_\mathrm{{Edd}}$
-rate objects show a positive correlation between Γ and
$\dot{m}_\mathrm{{Edd}}$
. This transition from harder-when-brighter at low
$\dot{m}_\mathrm{{Edd}}$
to softer-when-brighter at high
$\dot{m}_\mathrm{{Edd}}$
can be explained by a change in the dominant source of seed-photons for X-ray emission from cyclo-synchrotron emission from the Comptonizing corona itself to thermal seed-photons from the accretion disc. This transition is also seen in the ‘hard state’ of black hole X-ray binaries (BHXRBs). The results support the idea that low-ionization nuclear emission-line regions are analogues of BHXRBs in the hard state and that Seyferts are analogues of BHXRBs in either the high-accretion rate end of the hard state or in the hard-intermediate state.
Abstract
We present an X-ray time lag analysis, as a function of Fourier frequency, for MCG–6-30-15 and Mrk 766 using long-term XMM–Newton light curves in the 0.5–1.5 and the 2–4 keV energy bands, ...together with some physical modelling of the corresponding time lag spectra. Both the time lag spectra of MCG–6-30-15 and Mrk 766 show negative values (i.e. soft band variations lag behind the corresponding hard band variations) at high frequencies, around 10−3 Hz, similar to those previously observed from 1H 0707−495. The remarkable morphological resemblance between the time lag spectra of MCG–6-30-15 and Mrk 766 indicate that the physical processes responsible for the observed soft time delays are very similar in the two sources, favouring a reflection scenario from material situated very nearby to the central black hole.
We present an X-ray time lag analysis, as a function of Fourier frequency, for MCG-6-30-15 and Mrk766 using long-termXMM-Newtonlight curves in the 0.5-1.5 and the 2-4keV energy bands, together with ...some physical modelling of the corresponding time lag spectra. Both the time lag spectra of MCG-6-30-15 and Mrk766 show negative values (i.e. soft band variations lag behind the corresponding hard band variations) at high frequencies, around 10-3Hz, similar to those previously observed from 1H0707-495. The remarkable morphological resemblance between the time lag spectra of MCG-6-30-15 and Mrk766 indicate that the physical processes responsible for the observed soft time delays are very similar in the two sources, favouring a reflection scenario from material situated very nearby to the central black hole.
Context. Theoretical modelling of time-lags between variations in the Fe Kα emission and the X-ray continuum might shed light on the physics and geometry of the X-ray emitting region in active ...galaxies (AGN) and X-ray binaries. We here present the results from a systematic analysis of time-lags between variations in two energy bands (5−7 vs. 2−4 keV) for seven X-ray bright and variable AGN. Aims. We estimate time-lags as accurately as possible and fit them with theoretical models in the context of the lamp-post geometry. We also constrain the geometry of the X-ray emitting region in AGN. Methods. We used all available archival XMM-Newton data for the sources in our sample and extracted light curves in the 5−7 and 2−4 keV energy bands. We used these light curves and applied a thoroughly tested (through extensive numerical simulations) recipe to estimate time-lags that have minimal bias, approximately follow a Gaussian distribution, and have known errors. Using traditional χ2 minimisation techniques, we then fitted the observed time-lags with two different models: a phenomenological model where the time-lags have a power-law dependence on frequency, and a physical model, using the reverberation time-lags expected in the lamp-post geometry. The latter were computed assuming a point-like primary X-ray source above a black hole surrounded by a neutral and prograde accretion disc with solar iron abundance. We took all relativistic effects into account for various X-ray source heights, inclination angles, and black hole spin values. Results. Given the available data, time-lags between the two energy bands can only be reliably measured at frequencies between ~5 × 10-5 Hz and ~10-3 Hz. The power-law and reverberation time-lag models can both fit the data well in terms of formal statistical characteristics. When fitting the observed time-lags to the lamp-post reverberation scenario, we can only constrain the height of the X-ray source. The data require, or are consistent with, a small (≲ 10 gravitational radii) X-ray source height. Conclusions. In principle, the 5−7 keV band, which contains most of the Fe Kα line emission, could be an ideal band for studying reverberation effects, as it is expected to be dominated by the X-ray reflection component. We here carried out the best possible analysis with XMM-Newton data. Time-lags can be reliably estimated over a relatively narrow frequency range, and their errors are rather large. Nevertheless, our results are consistent with the hypothesis of X-ray reflection from the inner accretion disc.