Open Access
Issue
A&A
Volume 691, November 2024
Article Number A60
Number of page(s) 16
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/202348603
Published online 29 October 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

Over the past ten years, enormous observational efforts have been invested in the multi-wavelength monitoring of many nearby Seyferts with the Neil Gehrels Swift Observatory (hereafter, Swift). A large number of observational campaigns have already been conducted, monitoring the continuum variations from X-ray to ultraviolet (UV) and optical, at high cadence and over long periods of time (e.g. Shappee et al. 2014; McHardy et al. 2014, 2018; Fausnaugh et al. 2016; Cackett et al. 2020; Hernández Santisteban et al. 2020; Kara et al. 2021). The main objective of these observations is to estimate the cross-correlation between the variations seen in the X-ray and in the UV/optical and to constrain models for the continuum emission in active galactic nuclei (AGNs).

The hypothesis of X-ray illumination of the accretion disc was proposed many years ago to explain the detection of the Fe Kα emission line around 6.4 keV and the spectral hardening above ∼10 keV in the Ginga spectra of AGNs (e.g. Pounds et al. 1990; George & Fabian 1991). It was immediately recognised that if the accretion disc is illuminated by X-rays and the X-rays are variable, then the X-ray and disc emission should be variable in a correlated way – but with a delay that should be representative of the light travel time from the X-ray corona to the disc. This was indeed observed in the first multi-wavelength monitoring campaign of NGC 5548, which showed that the UV and X-ray flux variations were positively correlated and simultaneous to within a few days (Clavel et al. 1992).

We carried out a detailed study of the X-ray reverberation of the accretion disc in AGNs over the last few years. Kammoun et al. (2019) and Kammoun et al. (2021a, hereafter, K21a) studied the disc response to illuminating X-rays, taking into account all relativistic as well as the disc ionisation effects and then fitting (successfully) the time-lags as a function of wavelength in many sources (Kammoun et al. 2021b, 2023). Dovčiak et al. (2022) presented KYNSED, a new spectral model that can be used to fit the broadband spectral energy distributions (SEDs) of AGNs, from optical to X-rays, in the case where X-rays are illuminating the disc. These authors used the model to successfully fit the average SED of NGC 5548, using data from the STORM campaign (Fausnaugh et al. 2016). Papadakis et al. (2022) also used KYNSED and showed that X-ray illumination of the accretion disc can explain the discrepancy between the theoretical and the observed accretion disc sizes in micro-lensed quasars. Moreover, Kammoun et al. (2024) employed KYNSED and showed that X-ray illumination can explain both the time-averaged and time-resolved SEDs of NGC 5548.

Although the major objective of the monitoring campaigns so far has been the study of the delays between the optical/UV variations and those in the X-rays, the resulting light curves hold a wealth of additional information that can be used to constrain theoretical models. The time lags depend mainly on the width of the disc response function to a burst of X-ray radiation, but not on its amplitude. On the other hand, the variability amplitude in the optical/UV bands should depend on the amplitude of the disc response function and is an additional observational quantity that models should be able to explain. In fact, a model should be able to explain both the time lags and the variability amplitude of the observed light curves at the same time. This is was the approach adopted by Panagiotou et al. (2020, 2022) with respect to the STORM data of NGC 5548.

Our work is similar to Panagiotou et al. (2022, P22 hereafter), with the main difference being that P22 studied power spectra, whereas we focus on the variance of the light curves. In particular, we chose four sources from the sample of Kammoun et al. (2021b, K21b hereafter) and computed the excess variance of the light curves in various bands. Since the variance is the integral of the power spectrum, we used the P22 model transfer functions to compute the model variance in each band, assuming the model parameters that K21b determined by fitting the time-lags of the same sources. By comparing the model with the observed variances, we investigate whether the X-ray reverberation model can fit both the optical/UV time lags and observed variance in AGNs.

This work is organised as follows. In Section 2, we present the sample used in our work. In Section 3, we describe the X-ray reverberation model. In Sections 4 and 5, we explain how we compute the observed and model excess variance respectively. Finally, we compare the model variance spectra with the observations in Sect. 6. We summarise our work and discuss our results in Sect. 7.

2. The sample

The sources we study are taken from K21b and are listed in Table 1. We omitted NGC 5548, which was already studied by P22, NGC 4593, and NGC 7469 because the light curves are not available online. The luminosity distance, DL, the black hole mass, MBH, the 2−10 keV X-ray luminosity, LXobs, Edd, and photon index, Γ, are taken from K21b. The bolometric luminosity is used to compute the bending frequency in the X-ray power spectrum (see Sect. 5.2). It is computed assuming that L Bol = m ˙ Edd L Edd $ L_{\mathrm{Bol}}=\dot{m}_{\mathrm{Edd}}L_{\mathrm{Edd}} $, where LEdd is the Eddington luminosity and Edd is the accretion rate in Eddington units. We assumed the accretion rate estimate of Edelson et al. (2019) for Mrk 509, the values listed in Table 1 of K21b for NGC 4151 and NGC 2617, and 0.37 for Mrk 142 (with their computation explained in Sect. 6.6). These estimates are based on the use of the luminosity at a particular wavelength and various bolometric correction factors to transform it to LBol.

Table 1.

Redshift, z, luminosity distance, DL, BH mass, MBH, 2–10 keV X-ray luminosity, LXobs, Edd, and photon index, Γ, of the sources in our sample.

3. The X-ray reverberation model

We assumed a point-like X-ray source at height, h, on the rotational axis of the central black hole (BH). The source emits isotropically (in its rest frame) a spectrum of the form FX(t) = N(t)E−Γexp(−E/Ecut), where Γ and Ecut are assumed to be constant. X-rays illuminate the accretion disc, which co-rotates with the BH. A fraction of the X-rays falling on the disc will be re-emitted as X-rays (this is the ‘disc reflection component’), while the rest will be absorbed by the disc. They will thermalise in the disc and act as an extra source of heating. Consequently, the local disc temperature and emission will increase. As long as the X-rays are variable (which they indeed are in AGNs), the disc flux in the UV/optical wavebands will also be variable.

Following K21a, the disc flux in a wavelength λ, in the case of the X-ray illuminating discs, will be given by:

F λ ( t ) = F λ , NT + 0 Ψ λ ( t ) L Xobs , Edd ( t t ) d t , $$ \begin{aligned} F_{\lambda }(t) = F_{\rm \lambda , NT} + \int _{0}^{\infty } \Psi _{\lambda }(t\prime )L_{\rm Xobs,Edd}(t-t\prime )dt\prime , \end{aligned} $$(1)

where Fλ is the total disc flux at wavelength λ, Fλ, NT is the flux of a standard accretion disc when there is no illumination (Novikov & Thorne 1973), and LXobs, Edd is the observed X-ray luminosity in the 2–10 keV band, in units of Eddington luminosity (LEdd). Ψλ(t) is the disc response to an X-ray flash, defined by Eq. (3) in K21a. It is equal to the flux that the disc emits due to X-ray heating at time t after the X-ray flash. Therefore, the convolution integral in the right part of Eq. (1) is equal to the variable, X-ray thermally reprocessed flux of the disc when the X-rays are constantly variable.

3.1. The UV/optical power spectrum of X-ray illuminated discs

Equation (1) is significantly simplified when expressed in the frequency domain. It can be shown that (e.g. Papadakis et al. 2016, P22),

PSD λ ( ν ) = | Γ λ ( ν ) | 2 · PSD X , Edd ( ν ) , $$ \begin{aligned} \mathrm{PSD}_{\lambda }(\nu ) = |\Gamma _{\lambda }(\nu )|^{2} \cdot \mathrm{PSD}_{\rm X,Edd}(\nu ), \end{aligned} $$(2)

where PSDλ(ν) is the power spectrum at wavelength λ, PSDX, Edd is the X-ray power spectrum when the light curve is in units of LEdd, and Γλ(ν) is the transfer function, which is equal to the Fourier transform of the response function, Ψλ(t). Just as Ψλ(t), we know that Γλ(ν) depends on the geometric and physical properties of the X-ray source and the accretion disc; for example, MBH, Edd, X-ray luminosity, X-ray source height, and more (see P22). Consequently, by fitting the UV/optical PSDs when the X-ray PSD is known, we can determine Γλ(ν) and, hence, the physical parameters of the central source in AGNs.

This was the approach taken in P22 in their study of the optical and UV PSDs of NGC 5548. However, the estimation of the power spectrum is a non-trivial task, mainly because the light curves are not evenly sampled. Fitting the observed PSDs with Eq. (2) is not easy either, as shown in P22. We plan to perform a full PSD analysis of many sources in the near future (Panagiotou et al., in prep.). In the present work, we follow a simpler approach, namely: we measure the variance of the sources in the sample in the optical/UV bands and we fit the resulting ‘variability spectrum’ (i.e. variance versus photon frequency), assuming a wide range of model parameters. Then we compare our results with the best-fit results of K21b to investigate whether X-ray reverberation can explain both the time lags as well as the variability amplitude in AGNs with the same physical parameters.

3.2. The variance of X-ray illuminated discs in the UV/optical bands

The power spectrum of a random process, PSD(ν), is defined in such a way that

σ 2 = 0 PSD ( ν ) d ν , $$ \begin{aligned} \sigma ^{2} = \int _{0}^{\infty } \mathrm{PSD}(\nu ) d\nu , \end{aligned} $$(3)

where σ2 is the variance of the process. The variance we measure using observed light curves is not equal to the integral of the power spectrum over all frequencies, as we expect there to be intrinsic variations on long (and short) timescales that are not fully sampled by the observations. Allevato et al. (2013) showed that the variance (at a spectral band λ) measured using a light curve of duration T and bin size of Δt is an estimator of the so-called band variance, defined as:

σ λ , band 2 = 1 / T max 1 / T min PSD λ ( ν ) d ν , $$ \begin{aligned} \sigma ^{2}_{\lambda ,\mathrm{band}} = \int _{1/T_{\rm max}}^{1/T_{\rm min}} \mathrm{PSD}_{\lambda }(\nu ) d\nu , \end{aligned} $$(4)

where Tmax = T and Tmin = 2 × Δt, the longest and shortest sampled variability timescales respectively. Using Eq. (2), the equation above can be rewritten as follows

σ λ , band 2 = 1 / T max 1 / T min | Γ λ ( ν ) | 2 · PSD X , Edd ( ν ) d ν . $$ \begin{aligned} \sigma ^{2}_{\lambda ,\mathrm{band}} = \int _{1/T_{\rm max}}^{1/T_{\rm min}} |\Gamma _{\lambda }(\nu )|^{2} \cdot \mathrm{PSD}_{\rm X,Edd}(\nu ) d\nu . \end{aligned} $$(5)

This equation can be used to compute the model variance of the UV/optical light curves. The comparison between the model and the observed variance can test the validity of the model and can provide the set of model parameters that fit the data (i.e. the light curve variances) best.

4. Observed variance

4.1. The excess variance

One of the most common estimators of the intrinsic variance is the excess variance, σ XVobs 2 $ \sigma^{2}_{\mathrm{XVobs}} $,

σ XVobs 2 = 1 N i = 1 N [ ( x i x ¯ ) 2 σ err , i 2 ] , $$ \begin{aligned} \sigma ^{2}_{\rm XVobs} = \frac{1}{N} \sum _{i = 1}^{N}\left[(x_{i}-\bar{x})^2 - \sigma _{\mathrm{err}, i}^{2}\right], \end{aligned} $$(6)

according to Allevato et al. (2013), where xi, x ¯ $ \bar{x} $, and σerr, i are the count rate, the sample mean, and the error bar of the light curve points, respectively. We could compute the excess variance using the full light curve, however, single σXVobs2 measurements are not useful. Their distribution is very asymmetric and very broad, and their variance is unknown. For that reason, we divided the light curves into many segments, computing σXVobs2 for each one of them and then calculated their mean, σ 2 ¯ XVobs $ \overline{\sigma^2}_{\mathrm{XVobs}} $, and its error using Eqs. (12) and (13) in Allevato et al. (2013). According to these authors, as long as the signal-to-noise ratio of the light curves is higher than 3 (which is the case for all light curves, in all bands, in our case), the mean variance should be approximately Gaussian distributed, and its intrinsic standard deviation (i.e. its error) should be equal to the error of the mean variance, as long as the number of segments is ≳20.

Figures A.1, A.2, A.3, and A.4 show the light curves in X-rays and in at least one UV/optical band for all sources. We indicate the light curve segments we used to compute the variance, and we also plot σNXVobs2 of each segment. The σ 2 ¯ XVobs $ \overline{\sigma^2}_{\mathrm{XVobs}} $ measurements and their error are listed in Table 2. The first number in the parenthesis next to σ 2 ¯ XVobs $ \overline{\sigma^2}_{\mathrm{XVobs}} $ indicates the number of light curve segments we considered, while the second and third numbers list the duration of the segments, Tmax, and the average Δt of the light curves (note that, despite the excellent cadence of the Swift observations, the resulting light curves are not evenly sampled). The number of segments in some bands is smaller than 20. The shorter the Tmax, the larger the number of segments would be. However, if Tmax is too short, then the source variance is smaller than the Poisson noise variance and the σXVobs2 error bar is large. Therefore, the Tmax values we chose were a compromise between the need for the largest possible number of segments and the detection of the intrinsic variance with the smallest possible uncertainty.

4.2. Bias in the excess variance

In the case of red-noise-like PSDs (as is the case with AGNs), variations with timescales longer than Tmax can contribute to the measured variance, even if they are not fully sampled (this is known as the red-noise leakage). Consequently, the sample variance, σXVobs2, may overestimate the intrinsic variance, σ λ , band 2 $ \sigma^{2}_{\mathrm{\lambda, band}} $. This introduces a bias that we need to take into account when comparing the model band variance (computed using Eq. (5)) with σ 2 ¯ XVobs $ \overline{\sigma^2}_{\mathrm{XVobs}} $ (i.e. the mean of all variances we compute for each segment using Eq. (6)). According to Allevato et al. (2013) the bias depends on the PSD slope (say β) in the frequency range between 1/(2.1 Tmax) and 1/Tmax, in such a way that,

σ λ , band 2 = σ 2 ¯ XVobs × 0 . 48 ( β 1 ) . $$ \begin{aligned} \sigma ^{2}_{\rm \lambda , band} = \overline{\sigma ^{2}}_{\rm XVobs} \times 0.48^{(\beta - 1)}. \end{aligned} $$(7)

This means that before comparing the variance models with observations we need to take into account the correction factor of 0.48(β − 1). In all cases where it is necessary to correct for the bias of the measured variance, we can compute β by measuring the slope of a straight line between the model PSD points {log(1/(2.1 Tmax)), log[PSD(1/(2.1 Tmax))]} and {log(1/Tmax), log[PSD(1/Tmax)]}. A straight line can fit well the model PSDs in this relatively small frequency range in all cases we consider here. We note that β ∼ 1 − 3 for the sources studied here, with larger β values corresponding to larger wavelengths.

5. Model variance

According to Eq. (5) we need to compute the disc transfer function, Γλ(ν), and to know the X-ray power spectrum in order to compute the model variance that we will then compare with the observations. We describe below how we compute the model variance.

5.1. The disc transfer function

P22 studied in detail the accretion disc transfer functions in the case of X-ray-illuminated discs and showed that they can be well approximated by the following analytical formula:

| Γ λ ( ν ) | 2 = N 1 1 + ( ν ν b ) s , $$ \begin{aligned} |\Gamma _{\lambda }(\nu )|^{2} = N \frac{1}{1+ (\frac{\nu }{\nu _{b}})^s}, \end{aligned} $$(8)

where the normalisation N is in units of (erg/s/cm2/Å)2. According to the equation above, |Γλ(ν)|2 is constant at low frequencies and then it bends to a power-law shape, with a slope of s at frequencies higher than a characteristic frequency, νb. Equations (5) and (8) show that Γλ(ν) works like a low-pass frequency filter: all X-ray variations on time scales longer than ∼1/νb will also appear in the optical/UV bands. Shorter timescale X-ray variations are smeared out due to the X-ray reprocessing process. Consequently, the optical/UV PSD is suppressed at frequencies higher than ∼νb. The parameters N, s, and νb, are functions of the wavelength and the system’s physical parameters: MBH, Edd; the height and luminosity the X-ray source, h and LXobs, Edd; and the spectral photon index, Γ; the disc inclination θ; and the BH spin α* (see Appendix A of P22 for the exact forms of the analytical functions).

Equation (8) is a good approximation of the transfer function at low frequencies. It works well for the range of physical parameters and wave bands considered here, at frequencies below ∼4.3 day−1 (see Sect. 2.4 in P22), which corresponds to a sampling rate of around 0.23 day. This is comparable/smaller than the sampling rate of the sources in our sample (see Table 2).

Table 2.

Variance measurements, i.e. σ 2 ¯ HX , NXVobs $ \overline{\sigma^{2}}_{\mathrm{HX,NXVobs}} $, for the hard X-ray band (HX), and σ 2 ¯ XVobs $ \overline{\sigma^2}_{\mathrm{XVobs}} $ in the UV/optical bands, in units of 10−3 (mJy)2.

5.2. The X-ray power spectrum

We assume that the X-ray PSD has a bending power law shape, with a slope of −1 at low frequencies:

PSD X ( ν ) = A ν · 1 1 + ( ν ν b , X ) s X 1 , $$ \begin{aligned} \mathrm{PSD}_{\rm X}(\nu ) = \frac{A}{\nu } \cdot \frac{1}{1+ (\frac{\nu }{\nu _{\rm b,X}})^{s_{\rm X}-1}}, \end{aligned} $$(9)

where A is the normalisation, νb, X the bending frequency, and sX the spectral index above the bend. Equation (9) with sX = 2 fits the X-ray power spectra of AGNs well (e.g. Uttley et al. 2002; Markowitz et al. 2003; McHardy et al. 2004; González-Martín & Vaughan 2012). We note that the use of a different slope, from 1.5 up to 3, does not significantly affect our main results.

To compute the bend frequency, νb, X, we assume the McHardy et al. (2006) relation of,

log T b , X = a × log ( M BH , 6 ) b × log ( L Bol , 44 ) + c , $$ \begin{aligned} {\log }T_{\rm b,X} = a \times {\log }(M_{\rm BH,6}) - b \times {\log }(L_{\rm Bol,44}) + c, \end{aligned} $$(10)

where a = 2 . 17 0.25 + 0.32 $ a = 2.17^{+0.32}_{-0.25} $, b = 0 . 9 0.2 + 0.3 $ b = 0.9^{+0.3}_{-0.2} $, c = 2 . 42 0.25 + 0.22 $ c = -2.42^{+0.22}_{-0.25} $, Tb, X = 1/νb, X (in days), MBH, 6 is the black hole mass in units of 106 M, and LBol, 44 is the bolometric luminosity in units of 1044 erg/s.

As for the determination of the PSD normalisation, we can use the X-ray excess variance to compute A in Eq. (9) since, according to Eq. (4),

σ HX , band 2 = 1 / T max 1 / T min A ν · 1 1 + ( ν ν b , X ) s X 1 d ν . $$ \begin{aligned} \sigma ^{2}_{\rm HX,band} = \int _{1/T_{\rm max}}^{1/T_{\rm min}} \frac{A}{\nu } \cdot \frac{1}{1 + (\frac{\nu }{\nu _{\rm b,X}})^{s_{\rm X}-1}} d\nu . \end{aligned} $$(11)

The row ‘HX’1 in Table 2 lists the normalised excess variance of the X-ray light curves. Just like with the UV/optical light curves, we divided the X-ray light curves into many segments (their number indicated in parenthesis, next to the sample variances), and we computed the normalised excess variance of each one. The X-ray variance we list in Table 2 is the mean of the individual variance measurements of each light curve segment, σ 2 ¯ HX , NXVobs $ \overline{\sigma^{2}}_{\mathrm{HX,NXVobs}} $.

We can use the observed X-ray variance and the equation above to compute the PSD amplitude, A. However, as we mentioned in Sect. 4.2, σ HX , NXVobs 2 $ \sigma^{2}_{\mathrm{HX,NXVobs}} $ may be a biased estimator of σ HX , band 2 $ \sigma^{2}_{\mathrm{HX,band}} $, due to red noise leakage effects. We therefore used Eq. (7) and we corrected the observed variances for bias. Using then Eq. (11) (with sX = 2) it is straightforward to show that

A = σ HX , NXVobs 2 × 0 . 48 ( β 1 ) × [ ln ( ν b , X T max + 1 ν b , X T min + 1 ) ] 1 , $$ \begin{aligned} A = \sigma ^{2}_{\rm HX,NXVobs} \times 0.48^{(\beta -1)} \times \left[{\ln }\left(\frac{\nu _{\rm b,X} T_{\rm max} + 1}{\nu _{\rm b,X} T_{\rm min} + 1}\right)\right]^{-1}, \end{aligned} $$(12)

where β is the slope of the X-ray PSD, which was computed as described in Sect. 4.2. Computed in this way, A is the PSD normalisation for a light curve which is normalised to its mean. However, the PSD in Eq. (5) should be the power spectrum of the X-ray light curve in units of LEdd. The PSD in this case can be easily calculated by multiplying A by LXobs, Edd2.

When combining all the above, the final equation that gives the X-ray power spectrum is as follows:

PSD X , Edd ( ν ) = L Xobs , Edd 2 × σ HX , NXVobs 2 × 0 . 48 ( β 1 ) ln ( ν b , X T max + 1 ν b , X T min + 1 ) · ν 1 1 + ( ν ν b , X ) · $$ \begin{aligned} \mathrm{PSD}_{\rm X,Edd}(\nu ) = \frac{L_{\rm Xobs,Edd}^2\times \sigma ^{2}_{\rm HX,NXVobs} \times 0.48^{(\beta -1)}}{\ln \left(\frac{\nu _{\rm b,X} T_{\rm max} + 1}{\nu _{\rm b,X} T_{\rm min} + 1}\right)} \cdot \frac{\nu ^{-1}}{1 + (\frac{\nu }{\nu _{\rm b,X}})}\cdot \end{aligned} $$(13)

5.3. Calculation of the model band variance

We calculate the model band variance for the UV/optical bands, σλ, band2, using Eq. (5). The transfer function and the X-ray power spectrum are calculated as described in Sects. 5.1 and 5.2. However, there are a few more issues we must consider, before using Eq. (5) to compute σλ, band2.

First, since the variance in the UV/optical bands is given in physical units, we have also to consider the effects of Galactic (and host galaxy) interstellar absorption, which affects mainly the UV bands. To account for the effects of the Galactic interstellar reddening, we must multiply σ λ , band 2 $ \sigma^{2}_{\lambda,\mathrm{band}} $ by the square of the absorption factor fGal, λ = 10Aλ/2.5, where Aλ is the total extinction at a wavelength λ. Extinction was calculated using the reddening curves of Cardelli et al. (1989), with RV = 3.1. The Galactic E(B − V) values for each source in our sample were taken from Schlafly & Finkbeiner (2011). It is necessary to correct for reddening caused by the host galaxy in the case of NGC 4151, as we explain in Sect. 6. In this case, σ λ , band 2 $ \sigma^{2}_{\lambda,\mathrm{band}} $ should be multiplied by a second absorption factor as well, say f host , λ 2 $ f^{2}_{\mathrm{host},\lambda} $, which for simplicity was calculated like f Gal , λ 2 $ f^{2}_{\mathrm{Gal},\lambda} $ (for a different E(B − V), of course).

Second, the response functions of K21a that P22 used to compute the disc transfer functions assumed a distance of 1 Mpc. Therefore, σ λ , band 2 $ \sigma^{2}_{\lambda,\mathrm{band}} $ should be divided by D L 4 $ D^{4}_{\mathrm{L}} $, where DL is the luminosity distance (in Mpc).

Finally, to compare the model variance with σ 2 ¯ XVobs $ \overline{\sigma^2}_{\mathrm{XVobs}} $ we need to correct for bias. We use Eqs. (2), (8), and (13) to calculate the model power spectrum in each band λ and then we compute the slope β as described in Sect. 4.2. We then divide the model variance with a factor of 0.48(β − 1) to account for the bias of the variance measurements (given the model).

Taking into account the discussion in previous sections and the points above, the X-ray reverberation model band variance in the optical/UV bands becomes:

σ λ , band 2 = f Gal , λ 2 × f host , λ 2 D L 4 × 0 . 48 ( β 1 ) 1 / T max 1 / T min | Γ λ ( ν ) | 2 · PSD X , Edd ( ν ) · d ν , $$ \begin{aligned} \sigma ^{2}_{\lambda ,\mathrm{band}} = \frac{f^{2}_{\rm Gal,\lambda }\times f^{2}_{\rm host,\lambda }}{D^{4}_{\rm L} \times 0.48^{(\beta -1)}} \int _{1/T_{\rm max}}^{1/T_{\rm min}} |\Gamma _{\lambda }(\nu )|^{2} \cdot \mathrm{PSD}_{\rm X,Edd}(\nu ) \cdot d\nu , \end{aligned} $$(14)

where Γλ(ν) and PSDX, Edd(ν) are given by Eqs. (8) and (13), respectively. The model variance defined by the equation above is in units of (ergs/s/cm2/Å)2. It is computed using the transfer functions of P22, which are given in these units. However, we chose to compute the variance of the light curves in units of (mJy)2, because in this case, the variance versus photon frequency relation is flatter. This makes it easier to compare with the broadband spectral energy distribution of the objects. Therefore, in order to model the observed variance spectra, we need to transform the model σ λ , band 2 $ \sigma^{2}_{\lambda,\mathrm{band}} $ as computed above, to model variance in units of (mJy)2, σ ν , band 2 $ \sigma^{2}_{\nu,\mathrm{band}} $. This is done in the same way one transforms flux measurements from one to the other unit.

6. Fitting the model to the data

6.1. Data: Observed variance spectra

Black points in Figs. 1, 2, 3, and 4 show the observed variance spectra (i.e. σ 2 ¯ XVobs $ \overline{\sigma^2}_{\mathrm{XVobs}} $, taken from Table 2, plotted as a function of photon frequency), for Mrk 509, NGC 4151, NGC 2617, and Mrk 142, respectively. Points in red in Figs. 2 and 3 indicate variance measurements computed using different duration segments (the second number in parentheses in Table 2). The x-axis shows the rest frame frequency, that is, νrf = νobs(1 + z), where νobs = c/λobs, with λobs values listed in Table 2. We set the x-axis range to be the same in all sources, to indicate the difference between the spectral coverage of the available data for each source.

thumbnail Fig. 1.

Mrk 509 variance spectrum (i.e. the variance in units of (mJy)2 plotted as a function of photon frequency). Black points show the observed variance measurements calculated from light curve segments of duration, Tmax = 12 days (listed in Table 2). Open circles indicate data points excluded during the model fits. The black dashed line shows the best linear fit (in log-log space) to the data. The best-fit slope is b = 0.5 ± 0.4. The blue lines show the X-ray reverberation best-fit models for α* = 0 (left panel) and α* = 1 (right panel). The model’s M2 variance is lower than the adjacent W1 and W2 variances because of the ∼2200 Å feature in the Milky Way’s extinction law as determined by Cardelli et al. (1989). The bottom panels show the best-fit residuals, i.e. the observed variance and best fit model variance over the observed variance error, in the case of the X-ray reverberation best-fit models (blue lines). Grey squares on the top show the measured excess variance calculated using the full light curves, while the solid black line shows the best-fit model prediction for these measurements (see Sect. 6.4 for more details).

thumbnail Fig. 2.

NGC 4151 variance spectrum, similar to Fig. 1. The black points are calculated from light curve segments with duration Tmax = 3 days, while the red point with Tmax = 12 days. The blue line shows the best-fit variance model calculated with Tmax = 3. The slope of the best-fit black dashed line is b = −2.6 ± 0.4. The best-fit variance model for this source includes host galaxy absorption with E(B − V) = 0.35 (see Sect. 6.5 for more details). Grey points and solid black line as in Fig. 1.

thumbnail Fig. 3.

Same as in Fig. 1 for NGC 2617. We show σ 2 ¯ XVobs $ \overline{\sigma^2}_{\mathrm{XVobs}} $ with Tmax = 6 and 13 days with red triangle and black circle markers, respectively. We show the best-fit model variance for Tmax = 6 and 13 days with red and blue lines, respectively. The M2 variance is lower than the adjacent W1 and W2 variances for the same reason explained in Fig. 1. The slope of the best-fit dashed black line is b = 0.7 ± 0.8. The vertical dashed line indicates λ = 6000 Å. Grey points and solid black line as in Fig. 1.

thumbnail Fig. 4.

Same as in Fig. 1 for Mrk 142. The variance measurements have been calculated for Tmax = 6 days. Solid orange lines show the model predictions when assuming MBH = 2 × 106 and the parameters in K21b (see Sect. 6.6), while the blue lines show the best-fit model for MBH = 4 × 106 and the parameters listed in Tables 1 and 3. The slope of the best-fit dashed black line is b = 1.1 ± 0.2. Grey points and solid black line as in Fig. 1.

The figures show that the variance increases with increasing frequency (except for NGC 4151). This must be an intrinsic trend. It cannot be due to the increasing contribution of the host galaxy to the observed flux at lower frequencies because we computed the excess variance in physical units, and not normalised by the mean squared. Therefore, the square root of the sample variances indicates the average scatter of the light curve points around the mean in flux units. As such, it is equivalent to the amplitude of the variable part of the observed flux, only. The opposite trend in NGC 4151 is probably due to the fact that the active nucleus in this source is absorbed by gas and dust in the host galaxy, as we discuss below.

Firstly we fitted the observed variance spectra in the log-log space with a straight line of the form: log σ ν 2 = a + b ( log ν 15 ) $ \log\sigma^{2}_{\nu} = a + b (\log\nu - 15) $. During this model fit (and in all model fits in the next sections) we did not consider the SU band points in Mrk 509 and NGC 2617 (as variations of the Balmer continuum emitted by the broad line region, BLR, may contribute to the observed variations in this band), the measurements at wavelengths longer than 6000 Å (for reasons we explain in Sect. 7.2), and the SV σ 2 ¯ XVobs $ \overline{\sigma^2}_{\mathrm{XVobs}} $ in NGC 4151 (because the number of segments is small). Dashed lines in Figs. 1, 2, 3, and 4 show the best-fit lines. The best-fit slopes are: b = 0.5 ± 0.4, −2.6 ± 0.4, 0.7 ± 0.8, and 1.1 ± 0.2 for Mrk 509, NGC 4151, NGC 2617, and Mrk 142, respectively.

The variance spectra are best defined in the case of Mrk 142 (i.e. the best-fit slope is determined with the smallest error in this case). This is because the Mrk 142 light curves are long and well sampled and the source is highly variable even on relatively short timescales. This is not surprising, given the fact that Mrk 142 hosts the smallest BH mass. The variance measurements have larger errors, and the error of the best-fit slope is also large in NGC 2617. This is probably due to the fact that the Swift light curves of NGC 2617 are not sampled as densely as those of Mrk 142. The BH mass of NGC 4151 and NGC 2617 are comparable and although the duration of the NGC 4151 light curve segments is smaller than the respective duration in the NGC 2617 light curves, the NGC 4151 light curves are much better sampled, hence the variance spectrum is better defined. The Swift light curves of Mrk 509 are also well sampled, but the source is not highly variable on time scales of a few days (as the BH mass is large), so the error of the observed variance measurements (and of the best-fit slope) is relatively large.

A simple linear model fits the observed variance spectra well, but our objective is to fit the data with the X-ray reverberation model we have developed in the previous sections in order to constrain the physical parameters of the AGNs in the sample (i.e. accretion rate and X-ray corona height).

6.2. The model

We fitted the data plotted in Figs. 1, 2, 3, and 4 with the model σ ν , band 2 $ \sigma^{2}_{\nu,\mathrm{band}} $ as defined in Sect. 5. We note again that we choose to work in units of (mJy)2 because, in this way, the variance spectra are less steep and it is easier to observe differences between the model and the data visually. The DL are taken from Table 1, and fGal, λ is computed as explained in Sect. 5.3. We need to consider absorption by the host galaxy only in NGC 4151 (see Sect. 6.5). The disc transfer function, Γλ(ν), depends on α*, MBH, Edd, h, LXobs, Edd, Γ, and inclination, θ (see P22 for details). The BH mass, X-ray luminosity (in Eddington units) and the X-ray spectral photon index are fixed to the values listed in Table 1. As for the inclination, we assumed θ = 30° except for NGC 4151. Bentz et al. (2022) studied the geometry and kinematics of the BLR in this object. They found that the BLR is well described by a thick disc with similar opening angle and inclination angle of θ ∼ 57 − 58°, suggesting that our line of sight towards the central engine is just above the edge of the BLR. We therefore assumed θ = 60° for this object. Regarding the model X-ray power spectrum (i.e. PSDX, Edd(ν) in Eq. (14)), we used Eq. (13) with LXobs, Edd taken from Table 1, and σ 2 ¯ HX , NXVobs $ \overline{\sigma^2}_{\mathrm{HX,NXVobs}} $ from Table 2. We also computed νb from Eq. (10), using the BH mass and LBol values listed in Table 1. Consequently, the free parameters of the model are two: the accretion rate, Edd, and the corona height, h.

6.3. Best-fit results

We followed a standard χ2 minimisation procedure to fit the data in the case of α* = 0 and α* = 1. The data displayed with open points in Figs. 1, 2, 3, and 4 were not considered when fitting the data, for reasons explained in Sect. 6.1. The best-fit h and Edd values and their respective 1σ confidence region (in brackets), together with their χ2/degrees of freedom values, are listed in Table 3 in the third, fourth, and fifth columns respectively. Solid blue and red lines in Figs. 1, 2, 3, and 4 show the best-fit models over the full frequency range of the data in the plot, irrespective of whether we used the data points during the model fitting. We also show the best model fit residuals in the bottom panels of the same figures.

Table 3.

Best-fit results of this work.

The model fits the variance spectra of all sources well. The best-fit accretion rate in many cases is equal to 0.005, which is the smallest value we consider. This is the smallest accretion rate value for which the disc transfer functions of P22 are valid. It seems possible that the model could accommodate even lower Edd values while still providing a good fit to the data, albeit with different best-fit h values. This is because Edd and h are degenerate, as we explain below. However, our main objective is not to determine the parameter regions that explain the variance spectra on their own but to compare our best-fit results with those presented by K21b to investigate whether the X-ray reverberation hypothesis is sufficient to explain both the time lags and the variance spectra with the same physical parameter values. We can achieve this by comparing the confidence regions of the best-fit parameters as determined by us and by K21b.

thumbnail Fig. 5.

1σ confidence regions of (h, Edd) for each source. Blue and red regions indicate the confidence regions for α* = 0 and 1, respectively. Similarly, lighter blue and red shaded areas show the 1σ confidence regions as determined by K21b when fitting the time-lags, for α* = 0 and 1, respectively. In the case of Mrk 142 (lower right panel), the lighter shaded areas indicate the 1σ (h, Edd) confidence regions obtained from our best fit to the time-lags assuming MBH = 4 × 106 M (see Sect. 6.6).

The combined, 1σ confidence regions for Edd, and h are shown in Fig. 5 with blue-shaded and red-shaded regions (for the α* = 0 and α* = 1 best-fit models, respectively). The 1σ confidence regions were determined for a Δχ2 statistic of 2.3 which is equivalent to the 1σ confidence region for two interesting parameters (see e.g. Lampton et al. 1976). The figure shows that the 1σ confidence regions are rather large. This could be due to the fact that the sample variances have large errors in some cases, but this is probably not the main reason. For example, the variance spectrum of Mrk 142 has the smallest errors and the 1σ confidence region is still the largest. The large area of the confidence regions can be explained as follows.

The shape of the confidence regions indicates that the source height and the accretion rate are positively correlated. A larger accretion rate can fit the data equally well if the height of the X-ray source increases. This is because the amplitude of |Γλ(ν)|2 depends on h and Edd but in the opposite way: it increases with increasing h but it decreases with increasing Edd, (see Figs. 8 and 10 in P22, and the discussion in their Sects. 3.2 and 3.3). In fact, |Γλ(ν)|2, depends equally strongly on both parameters and it is for this reason that a rather large range of accretion rate and height values can fit the data equally well.

K21b fitted the time-lags spectrum of the sources we study using the same X-ray reverberation model. They also fitted the data keeping (h, Edd) as free parameters. The lighter red and blue shaded areas in Fig. 5 indicate their 1σ confidence region for these parameters (except for Mrk 142, where the respective confidence region is computed as explained in Sect. 6.6). The shape of the combined confidence regions in this case implies an anti-correlation between h and Edd (as opposed to the positive correlation implied from the model fits to the variance spectra). This is because the amplitude of the time-lags increases with increasing height and increasing accretion rate as well (see Figs. 17 and 18 in K21a). The range of the parameter values that can fit the time lags well is also large, for the same reason as above: the time lags depend equally strongly on the source height and the accretion rate; therefore, a large range of parameter values can provide an equally good fit to the time-lags.

However, the most important result of our work is that both confidence regions overlap. This result shows that X-ray reverberation can fit both the time lags and the variance spectra of these sources, with the same parameters. Furthermore, we can obtain much stronger constraints on the model parameters when fitting both the observed time lags and the variance spectra. The width of the projection of the intersection between the 1σ confidence regions, as defined from the time-lags and the variance spectra model fits, are listed in the third and fourth column in Table 3. The values in brackets in the first line below these results list the range of the projection of the 1σ confidence region resulting from the variance spectra model fits. For Mrk 142, we also list the 1σ range of the best-fit to the time lags (see Sect. 6.6). The size of the common area between the 1σ confidence regions should be representative of the 1σ confidence region of the model parameters that can fit both the time lags and the variance spectra. This is smaller than the range of the individual 1σ confidence regions, which shows that it is necessary to fit both the time lags and the variance spectra to constrain the model parameters as tight as possible in AGNs.

6.4. Excess variance in the full light curves

We calculated the excess variance of each source using the full light curves as well. The grey open squares in Figs. 1, 2, 3, and 4 show the results. In the case of Mrk 509 and NGC 4151, we used the full-length Swift light curves in all bands, as they are of the same duration. For NGC 2617, we used the Swift light curves from MJD 56640 to 56755 to avoid gaps. In the case of the r, i, and z bands, we used the light curve from the beginning of the ground-based observations until 115 days after the start, ensuring the same duration for the light curves in all wavebands. Similarly, for Mrk 142, we used light curves from MJD 58485 to 58600 in all wavebands.

The solid black lines in the same figures show the model predictions for the excess variance measurements when we consider the full light curves. We used the best-fit results of Table 3 and we took into account the larger Tmax, and the new σNXVobs2 measurements, computed from the full HX light curve. We note that the model black lines are less steep than the blue ones (in the same figures) since the bias correction term equals unity when the full light curves are used.

The excess variance measurements are in agreement with the model predictions in the case of NGC 2617 and they are larger than the model predictions by a factor of ∼10, 5.5, and 5 in the case of Mrk 509, NGC 4151, and Mrk 142, respectively. This is not surprising, given the fact that the distribution of the variance measurements from the full light curves is not Gaussian and resembles a χ2 distribution with low effective degrees of freedom (see e.g. the middle and bottom panels in Fig. 3 in Vaughan et al. 2003). In addition, the error of the full light curve variance measurements is unknown, as it depends on the value of the true variance of the intrinsic red-noise process (which is of course unknown; see e.g. Section 5.3.3 of Priestley 1981). These facts imply that we cannot properly compare the observed variance and the model predictions. Nevertheless, the overall shape of the model seems to be in good agreement with the observations, at least up to ∼6000 Å.

For example, the excess of variability amplitude we see in the SU band of Mrk 509 and NGC 2617 is present in the variance measurements computed using both the full light curves and the shorter segments. In contrast, no excess of variability amplitude in the same band appears in the variance spectra of NGC 4151 and Mrk 142. Furthermore, the long-term variance measurements at wavelengths longer than ∼6000 Å scatter erratically around the model predictions in NGC 2617 and in Mrk 142 (although with a small amplitude in the latter case). In NGC 4151 the measurements fall below the model predictions. These variations do not indicate a clear trend and, perhaps, they are due to the fact that the variability amplitude decreases with increasing wavelength in AGNs.

6.5. Notes on individual objects

Mrk 509: Figure 1 shows that X-ray reverberation can fit well the observed variance spectrum of this source both for α* = 0 and α* = 1, but for a different range of parameters. For example, the Edd 1σ region is (0.005−0.022) and (0.005−0.04) in the case of α* = 0 and 1, respectively. E19 list an accretion rate estimate of 0.05 for this source. They do not explain how they estimated Edd but it is probably computed with the use of a bolometric correction factor. This value is closer to the upper 1σṁEdd limit in the case of α* = 1. Although Edd estimates based on bolometric correction factors probably have large errors, this result could imply that the BH in Mrk 509 is rapidly spinning.

NGC 4151: Initially, we could not find model variance spectra with acceptable χ2 values for this source. However, it is well known that the central source in NGC 4151 is absorbed by neutral gas and dust (e.g. Ward et al. 1987). To take this into account, we modelled the host galaxy absorption using the reddening law given by Eq. (3) in Czerny et al. (2004), and we considered E(B − V) values in the range of 0−1, with a step of 0.05. For each E(B − V) value we computed fhost, λ2 (as explained in Sect. 5.3) and we applied it to the estimation of the model variance spectra (as shown in Eq. (14)). The best agreement with the data was achieved for E(B − V) = 0.35. Figure 2 shows that when accounting for intrinsic absorption the model can describe the data well.

An E(B − V) = 0.35 value translates to a hydrogen column density of N(H)≈2.9 × 1021 cm−2 (Liszt 2021). This is a rather small column density for the neutral absorber in NGC 4151 (e.g. Beuchert et al. 2017; Lubiński et al. 2016). However, N(H) is variable in this source. K21b fitted the mean X-ray spectrum of the source during the Swift monitoring campaign assuming an ionised absorber only. This implies that the column density of the neutral absorber may be relatively small during the Swift observations of NGC 4151 in 2016.

NGC 2617: Figure 3 shows the results for NGC 2617. The X-ray reverberation model fits the variance spectrum well. The number of the light curve segments we used to compute the excess variance is the smallest among the sources in our sample (in all bands). Consequently, the errors may not be Gaussian and the results from the model fitting when using χ2 statistics may not be accurate in this case. Figure 5 shows that the 1σ confidence region of the model parameters from the fits to the time-lags and the variance spectra barely overlap for α* = 1. This could be due to the fact that the error of the variance measurements may not be Gaussian (the range of the 1σ confidence region depends directly on the assumption of the errors being Gaussian). In any case, this is not a reason to claim that a maximally spinning BH is excluded for this source, as we consider the common area of the 1σ confidence regions only, which is not highly significant. On the other hand, F18 report an accretion rate estimate of 0.01 for this source. This is within the range of the common area between the 1σ regions that we report for α* = 0 in Table 3. This result could indicate a non-rotating BH in this AGN.

6.6. The case of Mrk 142

Black circles in Fig. 4 show the measured variance spectrum of Mrk 142. The orange-coloured solid lines show the model variance spectra calculated using the model parameters that K21b assumed for this source, namely: MBH = 2 × 106 M, LBol = 2.1 × 1044 erg/s, LXobs, Edd = 0.0246, and their best-fit h and Edd values, which are listed in their Table 3 (i.e. h = 51.1 and Edd = 0.64 for α* = 0, and h = 50.9 and Edd = 0.88 for α* = 1). The model predictions are significantly smaller than the observed variance spectrum. This disagreement could be a problem for the X-ray reverberation model, so we investigated possible reasons for this.

To increase the model variability amplitude we need to increase the normalisation of |Γλ(ν)|2. Figures 7–12 in P22 show how the transfer function amplitude depends on the various model parameters. The amplitude can increase with decreasing inclination (although the difference between the amplitude of |Γλ(ν)|2 with θ = 30° and θ = 0° is small), increasing Γ (however, with Γ already being very steep for this source), or increasing height. The amplitude of |Γ(ν)|2 will change by a factor of 2−3 if we increase the height from ∼50 to 100 Rg, which is the largest height that K21a and P22 considered; however, even this increase is not enough to explain the discrepancy between the model and the observed variance spectra in Fig. 4. So, we also considered the possibility that the BH mass is slightly larger than the value assumed by K21b. It is commonly assumed that the BH mass estimates in AGNs are affected by systematic errors of the order of ∼0.3 dex (i.e. by a factor of 2 or so, on average). We therefore investigated the possibility that MBH = 4 × 106 M in Mrk 142 (i.e. two times larger than the value assumed by K21b). We note that this value is consistent almost within 1σ with the BH mass measurements of Li et al. (2018) and Du et al. (2014).

We assumed a bolometric luminosity of LBol = 1.85 × 1044 erg/s, as computed by Cackett et al. (2020) when using the relation of L Bol = 9 λ L λ ( 5100 Å ) $ L_{\mathrm{Bol}} = 9\lambda L_{\lambda}(5100\,\AA) $ of Kaspi et al. (2000). We also considered the 2−10 keV X-ray flux measurement of K21b, but since we assumed a larger MBH, LXobs, Edd is now smaller. The chosen values for MBH, Lbol, and LXobs, Edd are listed in Table 1. For these parameters, the X-ray reverberation model can fit the variance spectrum of the source very well (best-fit results are listed in Sect. 6.3). The best-fit models are shown by the blue solid lines in Figure 4 and the 1σ confidence regions for h and Edd are shown in the lower-right panel of Figure 5. We note that the large difference between the blue and orange models in Fig. 4 is not only due to the change in mass. According to the results of Panagiotou et al. (2022), the amplitude of the transfer function (squared) increases with increasing BH mass as ∼MBH3. Therefore, the increase in BH mass implies that the blue line model should be about 8 times higher than the orange one. However, the increase in the BH mass implies a smaller accretion rate (for a fixed bolometric luminosity). This resulting decrease in the best-fit Edd explains the remaining difference between the two model lines.

Since the MBH and LXobs, Edd values that K21b assumed when fitting the Mrk 142 time-lags are different to the values we assumed when fitting the variance spectrum, we needed to investigate whether the X-ray reverberation model was also able to fit the observed time-lags of Mrk 142 when using the values listed in Table 1. To this end, we used the time lags model defined by Eqs. (1)–(9) of K21b and we fit the time lags spectrum by keeping the BH mass and the X-ray luminosity fixed to the values listed in Table 1. Figure 6 shows the observed time-lags with black points (data taken from C20) and the best-fit models with blue solid and red dashed lines for spin-zero and one, respectively. We did not consider the time-lags around ∼3000 Å when fitting the data (i.e. points plotted with grey open circles in Fig. 6) for the reasons explained above. The 1σ confidence region for (h, Edd) from the model fits to the time-lags is shown with the lighter blue and red shaded areas in the lower right panel of Fig. 5. The respective best-fit results, the 1σ confidence regions, and the χ2/d.o.f. values are listed in Table 3. If we consider the MBH and LBol values listed in Table 1, then Edd should be ∼0.4 for this source. This is higher than the upper limit of the common 1σ confidence region that we find for α* = 0. This could be an indication that just like Mrk 509, the BH in this source may be rotating; however, it is not rotating at a maximum rate, as the lower limit of the common 1σ region is larger than 0.4 for α* = 1.

thumbnail Fig. 6.

Time-lag spectrum of Mrk 142 (data taken from Cackett et al. 2020). The blue solid and red dashed lines show the best-fit model time-lags for α* = 0, 1, while the parameter values are listed in Table 1.

7. Summary and discussion

We studied the optical/UV variance spectra (i.e. the excess variance as a function of photon frequency) in four nearby Seyferts, namely: Mrk 509, NGC 4151, NGC 2617, and Mrk 142, using data from recent, long, multi-wavelength campaigns with Swift and ground-based telescopes (in two sources). We measured the excess variance in various spectral bands and compared it with the model variance predicted by the X-ray reverberation model. To compute the model variance, we employed the accretion disc transfer functions of Panagiotou et al. (2022) and assumed a bending power law for the X-ray power spectrum.

We found that the X-ray thermal reverberation model can fit well the variance spectra of the sources in the sample. However, the most important result of our work is that X-ray reverberation can explain both the variance spectra and the time-lags in these AGNs for the same physical parameters. To the best of our knowledge, this is the first time that the X-ray reverberation model has been shown to simultaneously fit well, with respect to both the time lags and variability spectra of an AGN.

7.1. On the X-ray reverberation model fits to the variance spectra and the time-lags

It has been stated many times in the past that X-ray reverberation predicts time lags that are smaller than the observed ones by a factor of approximately 2 to 3. This was first noticed by Edelson et al. (2015) and Fausnaugh et al. (2016), who studied the X-ray optical/UV times lags in NGC 5548. We have repeatedly shown in the past that X-ray reverberation can fit the observed time lags in many sources well. In fact, it can do so, for a large range of X-ray height and accretion rate values. In this study, we show that X-ray reverberation can also accurately account for the observed UV/optical variability amplitude in the four AGNs we studied (see also P22 for the case of NGC 5548). Furthermore, akin to the findings of K21b and Kammoun et al. (2023), we find that the X-ray reverberation model not only fits the variance spectra well but also does so across a broad range of h and Edd values (see Fig. 5).

K21b have discussed in detail what are the differences between the previous, analytical, X-ray reverberation time-lag models, and the X-ray reverberation models that were developed by K21a (and P22). A major difference that can significantly affect the model predictions for the time lags and variance spectra is that previous models assumed that the local ratio of external heating (due to X-ray illumination) to internal heating (due to accretion) is constant with the radius measurement. The observations show that this is not the case. One reason for the ability of the K21b and P22 models to fit well the observations is due to the fact that the model computes the external to internal heating at each radius, for any combination of h and Edd. This is important, as this ratio depends strongly on both of these physical parameters. As the accretion rate increases, the temperature of the disc rises as well. Consequently, for a constant X-ray luminosity and h, a smaller fraction of the observed emission is due to X-ray reprocessing (in each energy band); hence, the variability amplitude drops. On the other hand, as the corona height increases, the disc’s solid angle, as seen by the X-ray source, increases as well. As a result, for a constant Edd (and X-ray luminosity), the amount of X-ray emission being reprocessed in the disc is larger and the variability amplitude increases (in all bands, but not in the same way). It is exactly for these reasons that X-ray reverberation can fit the data for a large range of parameters (as we also explain in Sect. 6.3).

The main result of our work is that the X-ray reverberation model can fit well both the time lags and the variability amplitude of the AGNs we study here, for a reasonable range of parameters. However, we cannot be entirely certain about the range of the allowed parameter values, as we have used a particular X-ray reverberation model. Our modelling is based on the transfer functions of P22 which, in turn, were based on the disc response functions of K21b. The latter were computed only for two spin values, using a colour correction factor of fcol = 2.42 and assuming that the X-ray source is powered by an unknown source that is not related to the power released by the accretion process. As a result, the best-fit parameters reported here are only valid for the above assumptions as well. Kammoun et al. (2023) recently released a new code which can compute the disc transfer functions for any spin value and colour correction factor for when the X-ray source is powered by the accretion process or from a different source of power. We plan to use the new code for the study of the variability amplitude of AGNs in the future. We believe that our main result that X-ray reverberation can account for the variability amplitude and the time lags in these sources will remain valid, although the range of allowed parameters may change. For example, Kammoun et al. (2021b) found that the time lags of NGC 4593, and NGC 5548 can be aptly fitted with α* = 0 and 1 (assuming fcol = 2.4), while Kammoun et al. (2023) found that the time-lags of the same sources can be well fitted by also assuming α* = 0.7 and fcol = 1.7. According to the Kammoun et al. (2021b) results, the 1σ corona heights range in NGC 5548 and NGC 4593 is 23−59/33−78 Rg and 2.5−25/2.5−23 Rg (for α* = 0/1, respectively), while Kammoun et al. (2023) found heights greater than 64 Rg (for NGC 5548) and a range of 5.6−43 Rg for NGC 4593 for their solution of α* = 0, 0.7 and fcol = 1.7 (when 50% of the accretion power is transferred to the corona). Kammoun et al. (2024) showed that the X-ray reverberation model can also explain well the time-resolved SEDs of NGC 5548 for spin = 0 and fcol = 1.7. Although the range of the physical parameters may not be exactly equal to the values listed in Table 3 when we consider a different fcol and/or spin, the main result will remain the same: X-ray reverberation can explain both time-lags and the variance spectra in these sources.

As we argued above, the X-ray reverberation model can fit the time lags for various parameter values. In order to constrain the range of parameters that can explain the observational properties of AGNs, we should also consider additional observables. In a future publication, we plan to study time-lags, variance spectra and the mean, broadband SED, determined by the same data which are also used to compute time-lags and variance spectra.

7.2. Considering whether X-rays drive all the observed fast, optical/UV variations in AGNs

An excess of variability amplitude appears around ∼3000 Å in Mrk 509 and, perhaps, in NGC 2617 as well, but it is not statistically significant. An increase in the time-lags spectra3 in the same energy band is attributed to diffuse emission from gas in the much larger BLR (in terms of size), which should contribute to the observed flux between 2000−4000 Å (e.g. Netzer 2022; Korista & Goad 2019). If such a component is present, it should increase the variance in these wavelengths by an amount that should be representative of the percentage of this component’s flux over the total flux in the band. Our results show that such a component is not significant in the variance spectra of the four sources we study.

Figures 3 and 4 show that the measured variability amplitude may be larger than the model predictions approximately after λ = 6000 Å (marked with the vertical dashed lines) in NGC 2617 and Mrk 142. This could be an indication that X-ray reverberation cannot account for all the variable flux at large wavelengths. However, as K21a argued, their response functions at long wavelengths are truncated because they were computed for a disc with Rout = 104Rg. For example, Fig. 9 in K21a shows that the response functions at ∼10 000 Å for a 107 M black hole are truncated, even when Edd = 0.01. Such truncation effects could be more significant in NGC 2617 at long wavelengths (and even more significant if the accretion rate is larger than 0.01 of the Eddington ratio) and they could explain the slight disagreement between the model and the data variances at long wavelengths. An increased variability amplitude at longer wavelengths could also indicate that the disc is not entirely flat. If the disc height to radius ratio is constant, then we expect the disc height to increase at large radii. In this case, X-ray absorption will increase; hence, the amplitude of the variable flux component will be larger than that predicted by the P22 models, which assumed a flat disc at all radii. This effect should probably be more pronounced at longer wavelengths as the disc height will be quite large at large radii. It is for these reasons that we did not include σ 2 ¯ XVobs $ \overline{\sigma^2}_{\mathrm{XVobs}} $ measurements at wavelengths longer than 6000 Å in the model fits. However, we stress that the disagreement between the best-fit model and the variance measurements at long wavelengths is not significant.

We therefore conclude that all of the observed optical/UV variability in these AGNs, on the sampled time scales (i.e. a few days up to a week) could be due to X-ray reverberation. It is difficult to explain such fast variability with intrinsic disc variations, as the accretion disc characteristic time scales in AGNs are much longer. For example, the viscous and thermal time-scale at a distance of 100 Schwartzchild radii (so that any variation at this distance will affect both the UV and the optical emitting areas of the disc) is of the order of 800 and 80 days, respectively, for a black hole of 107 M mass. Our result of an intrinsically non-variable disc is in agreement with results from variability studies of Galactic X-ray black-hole binaries. The accretion disc in their soft state appears to be intrinsically constant (e.g. Churazov et al. 2001), as is the case for the AGNs we consider here. However, the accretion disc may be intrinsically variable on longer time scales. A variability analysis of light curves on longer time scales will allow us to address this possibility (in a future publication).

7.3. Height of the X-ray source and the BH spin in AGNs

Our results can put some interesting constraints on the height of the X-ray source and the BH and spin in AGNs. In agreement with K21b, our results indicate that the BH mass may be rotating in Mrk 509 and Mrk 142. This is because the accretion rate estimates for these objects (e.g. 0.05 for Mrk 509, from E19, and 0.4 for Mrk 142, computed as explained in Sect. 6.6) are larger than the upper limit of the 1σ confidence region that is listed in Table 3 (for α* = 0). In fact, we computed the common 3σ confidence regions for Edd and height, and we found that these accretion rate estimates are not consistent with the hypothesis of α* = 0 at the 3σ level. On the other hand, the accretion rate estimate of 0.01 for NGC 2617 is within the 1σ confidence region if we assume α* = 0. In this case, we can exclude a maximally spinning BH at the 2.5σ level.

Our results suggest that the height of the X-ray source is ∼10 − 20 Rg in Mrk 509, NGC 2617, and NGC 4151. These are relatively low heights, close to the lowest height limit below which it would be difficult to detect the UV/optical variations due to X-ray reverberation in these objects. Indeed, as K21a showed (see their Fig. 26), the fraction of the variable flux (due to X-ray reverberation) over the total UV/optical flux decreases significantly with decreasing height below ∼20 Rg. The X-ray source height in Mrk 142 turns out to be larger than 70 Rg. This is the source with the smallest BH mass and the largest accretion rate among the sources in the sample. However, our sample is too small for us to investigate whether the source height depends either on BH mass or accretion rate. We need a larger sample to start investigating these trends.

Our work highlights the importance of modelling the optical/UV variability in the long, densely sampled light curves from the recent, multi-wavelength monitoring campaigns of many AGNs. Most of the analysis so far has focussed on the study of the correlation between the X-ray and the optical/UV light curves. However, the study of either the variance (as we do in this work) and/or the power spectrum as a function of wavelength can also help constrain theoretical models. We believe that the study of both the time lags and of the variability amplitude as a function of wavelength, simultaneously, is the strongest test of any model that tries to explain the optical/UV variability in AGNs. We show that the X-ray reverberation can explain both the time lags and the variance spectra for a reasonable set of model parameter values. We plan to study more objects in the future, with the aim of determining the physical parameters of the inner engine in these objects.


1

HX refers to the ‘hard’ X-ray band, which is not the same in all sources. E19 extracted the Mrk 509 light curve in the 1.5–10 keV band, while we use the 5–10 keV band light curve of E17 for NGC 4151 (the source is absorbed in lower energies). As for NGC 2617 and Mrk 142, we consider the 0.3–10 keV band light curves of F18 and C20. Note: the broadband variance in AGNs does not depend significantly on the energy limits of the band (see e.g. Fig. 2 in Ponti et al. 2012).

2

This means that each disc element is assumed to emit as a modified blackbody, that is a blackbody normalised by a factor 1/fcol4 and with temperature boosted by a factor fcol, thus keeping the same total emitted flux.

3

A small excess at around 3000 Å appears in the time-lags spectra of the sources in our sample, except for NGC 4151 (see Fig. 3 in K21b). K21b did not take into account the time-lags at this wavelength when fitting the time-lags spectra, although the discrepancy between the measured time-lags at this wavelength. The best-fit model was mildly significant only in the case of Mrk 509.

Acknowledgments

We thank the referee for their comments which helped us improve the manuscript. M.P. acknowledges support from the European Research Council (ERC) under the European Union Horizon 2020 research and innovation program under grant agreement No. 771282.

References

  1. Allevato, V., Paolillo, M., Papadakis, I., & Pinto, C. 2013, ApJ, 771, 9 [NASA ADS] [CrossRef] [Google Scholar]
  2. Bentz, M. C., Williams, P. R., & Treu, T. 2022, ApJ, 934, 168 [CrossRef] [Google Scholar]
  3. Beuchert, T., Markowitz, A. G., Dauser, T., et al. 2017, A&A, 603, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Cackett, E. M., Gelbord, J., Li, Y.-R., et al. 2020, ApJ, 896, 1 [NASA ADS] [CrossRef] [Google Scholar]
  5. Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [Google Scholar]
  6. Churazov, E., Gilfanov, M., & Revnivtsev, M. 2001, MNRAS, 321, 759 [NASA ADS] [CrossRef] [Google Scholar]
  7. Clavel, J., Nandra, K., Makino, F., et al. 1992, ApJ, 393, 113 [NASA ADS] [CrossRef] [Google Scholar]
  8. Czerny, B., Li, J., Loska, Z., & Szczerba, R. 2004, MNRAS, 348, L54 [NASA ADS] [CrossRef] [Google Scholar]
  9. Dovčiak, M., Papadakis, I. E., Kammoun, E. S., & Zhang, W. 2022, A&A, 661, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Du, P., Hu, C., Lu, K.-X., et al. 2014, ApJ, 782, 45 [NASA ADS] [CrossRef] [Google Scholar]
  11. Edelson, R., Gelbord, J. M., Horne, K., et al. 2015, ApJ, 806, 129 [NASA ADS] [CrossRef] [Google Scholar]
  12. Edelson, R., Gelbord, J., Cackett, E., et al. 2017, ApJ, 840, 41 [NASA ADS] [CrossRef] [Google Scholar]
  13. Edelson, R., Gelbord, J., Cackett, E., et al. 2019, ApJ, 870, 123 [NASA ADS] [CrossRef] [Google Scholar]
  14. Fausnaugh, M. M., Denney, K. D., Barth, A. J., et al. 2016, ApJ, 821, 56 [Google Scholar]
  15. Fausnaugh, M. M., Starkey, D. A., Horne, K., et al. 2018, ApJ, 854, 107 [NASA ADS] [CrossRef] [Google Scholar]
  16. George, I. M., & Fabian, A. C. 1991, MNRAS, 249, 352 [Google Scholar]
  17. González-Martín, O., & Vaughan, S. 2012, A&A, 544, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Hernández Santisteban, J. V., Edelson, R., Horne, K., et al. 2020, MNRAS, 498, 5399 [CrossRef] [Google Scholar]
  19. Kammoun, E. S., Papadakis, I. E., & Dovčiak, M. 2019, ApJ, 879, L24 [NASA ADS] [CrossRef] [Google Scholar]
  20. Kammoun, E. S., Dovčiak, M., Papadakis, I. E., Caballero-García, M. D., & Karas, V. 2021a, ApJ, 907, 20 [NASA ADS] [CrossRef] [Google Scholar]
  21. Kammoun, E. S., Papadakis, I. E., & Dovčiak, M. 2021b, MNRAS, 503, 4163 [CrossRef] [Google Scholar]
  22. Kammoun, E. S., Robin, L., Papadakis, I. E., Dovčiak, M., & Panagiotou, C. 2023, MNRAS, 526, 138 [NASA ADS] [CrossRef] [Google Scholar]
  23. Kammoun, E., Papadakis, I. E., Dovčiak, M., & Panagiotou, C. 2024, A&A, 686, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Kara, E., Mehdipour, M., Kriss, G. A., et al. 2021, ApJ, 922, 151 [NASA ADS] [CrossRef] [Google Scholar]
  25. Kaspi, S., Smith, P. S., Netzer, H., et al. 2000, ApJ, 533, 631 [Google Scholar]
  26. Korista, K. T., & Goad, M. R. 2019, MNRAS, 489, 5284 [CrossRef] [Google Scholar]
  27. Lampton, M., Margon, B., & Bowyer, S. 1976, ApJ, 208, 177 [Google Scholar]
  28. Li, Y.-R., Songsheng, Y.-Y., Qiu, J., et al. 2018, ApJ, 869, 137 [CrossRef] [Google Scholar]
  29. Liszt, H. 2021, ApJ, 908, 127 [NASA ADS] [CrossRef] [Google Scholar]
  30. Lubiński, P., Beckmann, V., Gibaud, L., et al. 2016, MNRAS, 458, 2454 [Google Scholar]
  31. Markowitz, A., Edelson, R., Vaughan, S., et al. 2003, ApJ, 593, 96 [NASA ADS] [CrossRef] [Google Scholar]
  32. McHardy, I. M., Papadakis, I. E., Uttley, P., Page, M. J., & Mason, K. O. 2004, MNRAS, 348, 783 [NASA ADS] [CrossRef] [Google Scholar]
  33. McHardy, I. M., Koerding, E., Knigge, C., Uttley, P., & Fender, R. P. 2006, Nature, 444, 730 [NASA ADS] [CrossRef] [Google Scholar]
  34. McHardy, I. M., Cameron, D. T., Dwelly, T., et al. 2014, MNRAS, 444, 1469 [NASA ADS] [CrossRef] [Google Scholar]
  35. McHardy, I. M., Connolly, S. D., Horne, K., et al. 2018, MNRAS, 480, 2881 [NASA ADS] [CrossRef] [Google Scholar]
  36. Netzer, H. 2022, MNRAS, 509, 2637 [NASA ADS] [Google Scholar]
  37. Novikov, I. D., & Thorne, K. S. 1973, Black Holes (Les Astres Occlus), 343 [Google Scholar]
  38. Panagiotou, C., Papadakis, I. E., Kammoun, E. S., & Dovčiak, M. 2020, MNRAS, 499, 1998 [CrossRef] [Google Scholar]
  39. Panagiotou, C., Papadakis, I., Kara, E., Kammoun, E., & Dovčiak, M. 2022, ApJ, 935, 93 [NASA ADS] [CrossRef] [Google Scholar]
  40. Papadakis, I., Pecháček, T., Dovčiak, M., et al. 2016, A&A, 588, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Papadakis, I. E., Dovčiak, M., & Kammoun, E. S. 2022, A&A, 666, A11 [Google Scholar]
  42. Ponti, G., Papadakis, I., Bianchi, S., et al. 2012, A&A, 542, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Pounds, K. A., Nandra, K., Stewart, G. C., George, I. M., & Fabian, A. C. 1990, Nature, 344, 132 [NASA ADS] [CrossRef] [Google Scholar]
  44. Priestley, M. 1981, Spectral Analysis and Time Series, Probability and Mathematical Statistics (Academic Press) [Google Scholar]
  45. Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 [Google Scholar]
  46. Shappee, B. J., Prieto, J. L., Grupe, D., et al. 2014, ApJ, 788, 48 [Google Scholar]
  47. Uttley, P., McHardy, I. M., & Papadakis, I. E. 2002, MNRAS, 332, 231 [Google Scholar]
  48. Vaughan, S., Edelson, R., Warwick, R. S., & Uttley, P. 2003, MNRAS, 345, 1271 [Google Scholar]
  49. Ward, M. J., Geballe, T., Smith, M., Wade, R., & Williams, P. 1987, ApJ, 316, 138 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: The light curves and variance measurements of the sources

thumbnail Fig. A.1.

X-ray and Swift W2 light curves of Mrk 509 (data taken from E19). The green dashed and red solid lines indicate the starting and ending points of the light curve segments that we used to compute σNXVobs2 and σXVobs2. The segments we used to compute the excess variance in the other optical bands are very similar to the ones indicated in the panel which shows the W2 light curve. The black points in the bottom panels of each plot show the variance measurements of each segment. The grey dashed lines and the grey area show the mean variance and its error, respectively. We draw black dotted lines at σNXVobs2 = 0 and σXVobs2 = 0.

thumbnail Fig. A.2.

Same as in Fig. A.1 for NGC 4151 (data are taken from E17). In this case, we plot the Swift V band light curve as well, as the width of the segments in this light curve is different than the width of the segments in the Swift W2 and the other band light curves.

thumbnail Fig. A.3.

Same as in Fig. A.1 for NGC 2617 (data taken from F18). To extend the spectral coverage of the variance spectra, we also considered light curves in the r, i, and z−bands for this source (data from the same source). For this reason, we also plot the r−band light curve in this figure. The segments we used to compute the excess variance in the i and z−band light curves are very similar to the ones indicated in the r−band light curve.

thumbnail Fig. A.4.

Same as in Fig. A.3 for Mrk 142 (data taken from C20). To increase the spectra coverage and the number of variance data within the spectral width of the Swift filters we also used light curves from ground-based telescopes for this source. These light curves are well sampled and significantly longer than the Swift light curves. We plot the r−band light curve in this figure; the segments we used to compute the excess variance in the u, g, i and z−band light curves (data taken from C20 as well) are very similar to the ones indicated in the r−band light curve.

All Tables

Table 1.

Redshift, z, luminosity distance, DL, BH mass, MBH, 2–10 keV X-ray luminosity, LXobs, Edd, and photon index, Γ, of the sources in our sample.

Table 2.

Variance measurements, i.e. σ 2 ¯ HX , NXVobs $ \overline{\sigma^{2}}_{\mathrm{HX,NXVobs}} $, for the hard X-ray band (HX), and σ 2 ¯ XVobs $ \overline{\sigma^2}_{\mathrm{XVobs}} $ in the UV/optical bands, in units of 10−3 (mJy)2.

Table 3.

Best-fit results of this work.

All Figures

thumbnail Fig. 1.

Mrk 509 variance spectrum (i.e. the variance in units of (mJy)2 plotted as a function of photon frequency). Black points show the observed variance measurements calculated from light curve segments of duration, Tmax = 12 days (listed in Table 2). Open circles indicate data points excluded during the model fits. The black dashed line shows the best linear fit (in log-log space) to the data. The best-fit slope is b = 0.5 ± 0.4. The blue lines show the X-ray reverberation best-fit models for α* = 0 (left panel) and α* = 1 (right panel). The model’s M2 variance is lower than the adjacent W1 and W2 variances because of the ∼2200 Å feature in the Milky Way’s extinction law as determined by Cardelli et al. (1989). The bottom panels show the best-fit residuals, i.e. the observed variance and best fit model variance over the observed variance error, in the case of the X-ray reverberation best-fit models (blue lines). Grey squares on the top show the measured excess variance calculated using the full light curves, while the solid black line shows the best-fit model prediction for these measurements (see Sect. 6.4 for more details).

In the text
thumbnail Fig. 2.

NGC 4151 variance spectrum, similar to Fig. 1. The black points are calculated from light curve segments with duration Tmax = 3 days, while the red point with Tmax = 12 days. The blue line shows the best-fit variance model calculated with Tmax = 3. The slope of the best-fit black dashed line is b = −2.6 ± 0.4. The best-fit variance model for this source includes host galaxy absorption with E(B − V) = 0.35 (see Sect. 6.5 for more details). Grey points and solid black line as in Fig. 1.

In the text
thumbnail Fig. 3.

Same as in Fig. 1 for NGC 2617. We show σ 2 ¯ XVobs $ \overline{\sigma^2}_{\mathrm{XVobs}} $ with Tmax = 6 and 13 days with red triangle and black circle markers, respectively. We show the best-fit model variance for Tmax = 6 and 13 days with red and blue lines, respectively. The M2 variance is lower than the adjacent W1 and W2 variances for the same reason explained in Fig. 1. The slope of the best-fit dashed black line is b = 0.7 ± 0.8. The vertical dashed line indicates λ = 6000 Å. Grey points and solid black line as in Fig. 1.

In the text
thumbnail Fig. 4.

Same as in Fig. 1 for Mrk 142. The variance measurements have been calculated for Tmax = 6 days. Solid orange lines show the model predictions when assuming MBH = 2 × 106 and the parameters in K21b (see Sect. 6.6), while the blue lines show the best-fit model for MBH = 4 × 106 and the parameters listed in Tables 1 and 3. The slope of the best-fit dashed black line is b = 1.1 ± 0.2. Grey points and solid black line as in Fig. 1.

In the text
thumbnail Fig. 5.

1σ confidence regions of (h, Edd) for each source. Blue and red regions indicate the confidence regions for α* = 0 and 1, respectively. Similarly, lighter blue and red shaded areas show the 1σ confidence regions as determined by K21b when fitting the time-lags, for α* = 0 and 1, respectively. In the case of Mrk 142 (lower right panel), the lighter shaded areas indicate the 1σ (h, Edd) confidence regions obtained from our best fit to the time-lags assuming MBH = 4 × 106 M (see Sect. 6.6).

In the text
thumbnail Fig. 6.

Time-lag spectrum of Mrk 142 (data taken from Cackett et al. 2020). The blue solid and red dashed lines show the best-fit model time-lags for α* = 0, 1, while the parameter values are listed in Table 1.

In the text
thumbnail Fig. A.1.

X-ray and Swift W2 light curves of Mrk 509 (data taken from E19). The green dashed and red solid lines indicate the starting and ending points of the light curve segments that we used to compute σNXVobs2 and σXVobs2. The segments we used to compute the excess variance in the other optical bands are very similar to the ones indicated in the panel which shows the W2 light curve. The black points in the bottom panels of each plot show the variance measurements of each segment. The grey dashed lines and the grey area show the mean variance and its error, respectively. We draw black dotted lines at σNXVobs2 = 0 and σXVobs2 = 0.

In the text
thumbnail Fig. A.2.

Same as in Fig. A.1 for NGC 4151 (data are taken from E17). In this case, we plot the Swift V band light curve as well, as the width of the segments in this light curve is different than the width of the segments in the Swift W2 and the other band light curves.

In the text
thumbnail Fig. A.3.

Same as in Fig. A.1 for NGC 2617 (data taken from F18). To extend the spectral coverage of the variance spectra, we also considered light curves in the r, i, and z−bands for this source (data from the same source). For this reason, we also plot the r−band light curve in this figure. The segments we used to compute the excess variance in the i and z−band light curves are very similar to the ones indicated in the r−band light curve.

In the text
thumbnail Fig. A.4.

Same as in Fig. A.3 for Mrk 142 (data taken from C20). To increase the spectra coverage and the number of variance data within the spectral width of the Swift filters we also used light curves from ground-based telescopes for this source. These light curves are well sampled and significantly longer than the Swift light curves. We plot the r−band light curve in this figure; the segments we used to compute the excess variance in the u, g, i and z−band light curves (data taken from C20 as well) are very similar to the ones indicated in the r−band light curve.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.