Optical Variability of Active Galactic Nuclei
- Astronomical Observatory, University of Warsaw, Warsaw, Poland
Variability studies of active galactic nuclei (AGNs) typically use either power spectral density (PSD) and structure function (SF) analyses or direct modeling of light curves with the damped random walk (DRW) and the continuous autoregressive moving average (CARMA) models. A fair fraction of research publications on the subject are flawed, and simply report incorrect results, because they lack a deep understanding of where these methods originate from and what their limitations are. For example, SF analyses typically lack or use a wrong noise subtraction procedure, leading to flat SFs. DRW, on the other hand, can only be used if the experiment length is sufficient, at least ten times the signal decorrelation time scale τ, and if the data show the power-law SF slope of γ ≡ 0.5.
2. Structure Functions
The structure function (SF) analysis is a model-independent technique of converting an active galactic nucleus (AGN) light curve into a different space, the variability amplitude–timescale space. The basic approach behind the SF analysis is as follows. Data points yi in an AGN light curve are, in the simplest case, a sum of the variable signal si (with the variance ) and the observational noise ni (with the variance ), so yi = si + ni. SF originates from simple mathematical properties of the covariance of the light curve (index i) with a shifted copy of itself (index j) by the timelag Δt = ti − tj, via (MacLeod et al., 2010; Kozłowski, 2016b)
where SF(Δt) is typically measured from data as
In order to measure the true AGN variability, so in fact cov(si, sj) in Equation (1), one needs to subtract the full noise term () from the SF in Equation (2). This is either rarely done in recent works or done incorrectly, as commonly only a fraction of the noise term () is subtracted. This leads to flat power-law SF slopes of γ = 0.1–0.4 at short timescales Δt (SF(Δt) ∝ Δtγ) (e.g., Vanden Berk et al., 2004; de Vries et al., 2005), but when correctly measured, the SF slope in optical is significantly steeper γ = 0.55 ± 0.08, based on ~9,200 SDSS AGN from Stripe 82 (Kozłowski, 2016b) and γ ≈ 0.45 in mid-IR (Kozłowski et al., 2010a, 2016). An equally important variability observable to the SF slope is the decorrelation timescale τ, a timescale at which the SF changes slope from the red noise (γ = 0.5) to the white noise (γ = 0.0). It seems to be about one year rest-frame, again based on ~9,200 SDSS AGN from Stripe 82 (Kozłowski, 2016b). I recently proposed a new method of the measurement of the unbiased decorrelation timescale τ from SFs (Kozłowski, 2017a). Another SF observable is the AGN variability amplitude measured at 1 year (rest-frame) with the value of 0.20 ± 0.06 mag in optical bands, while the asymptotic variability amplitude at long timescales (Δt ≫ τ, so Δt ≫ 1 year rest-frame) is 0.25±0.06 mag (Kozłowski, 2016b). The SF amplitude at 1 year may be affected, while the asymptotic variability amplitude is not, by the bias due to the unknown underlying stochastic process for short datasets (Kozłowski, 2017a).
3. The Damped Random Walk
AGN light curves can be modeled and interpolated using the damped random walk (DRW) stochastic process. DRW modeling (Kelly et al., 2009; Kozłowski et al., 2010b; MacLeod et al., 2010) by definition assumes an exponential covariance matrix of the signal of the form
that again by definition produces a fixed SF power-law slope of γ ≡ 0.5 at timelags Δt = ti − tj shorter than the signal decorrelation timescale τ (Kozłowski, 2016b). If the variability signal is due to a different stochastic process, where the SF slope is shallower/steeper than γ = 0.5, DRW will obtain a reasonable fit, however, it will report biased measurements (Kozłowski, 2016a). As of now, there is no statistical correction available to this problem, however, using the information on the SF/PSD slope one can modify the DRW model covariance matrix (Equation 3) and model the light curves with the modified DRW model to obtain correct parameters. DRW should be used to model an AGN light curve if one is convinced that the SF slope for a light curve is γ = 0.5 (or equivalently the PSD slope is −2). There is another issue with DRW: if the light curve is shorter than 10τ (~10 years rest-frame), it will simply report meaningless variability parameters (Kozłowski, 2017b). DRW is the simplest of the CARMA models [i.e., DRW ≡ CARMA(1, 0)], therefore the whole CARMA model family is plausibly affected by biases or problems reported above.
Constraining the SF, PSD, and DRW (or more generally the Gaussian processes) parameters typically require long and well-sampled AGN light curves. Such tight constraints may soon be available from the OGLE Sky Survey (Udalski et al., 2015), that has been monitoring the sky for 25 years, and in particular from its 20-year-long monitoring of nearly 1,000 AGNs (each with ~1,000 epochs), discovered mostly by the Magellanic Quasars Survey (Kozłowski et al., 2013). Similarly to the results from the SDSS Stripe 82, preliminary results from OGLE point to the mean SF slope γ ≳ 0.5. For sparsely sampled or short light curves some corrections to improve biases in PSD/SF/excess variance measurements are available (Vaughan et al., 2003; Allevato et al., 2013), although in a statistical (ensemble) sense, rather than for individual objects.
The author confirms being the sole contributor of this work and approved it for publication.
SK acknowledges the financial support of the Polish National Science Center through the OPUS grant 2014/15/B/ST9/00093 and MAESTRO grant 2014/14/A/ST9/00121.
Conflict of Interest Statement
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Kozłowski, S., Kochanek, C. S., Ashby, M. L. N., Assef, R. J., Brodwin, M., Eisenhardt, P. R., et al. (2016). Quasar variability in the mid-infrared. Astrophys. J. 817:119. doi: 10.3847/0004-637X/817/2/119
Kozłowski, S., Kochanek, C. S., Stern, D., Ashby, M. L. N., Assef, R. J., Bock, J. J., et al. (2010a). Mid-infrared variability from the spitzer deep wide-field survey. Astrophys. J. 716, 530–543. doi: 10.1088/0004-637X/716/1/530
Kozłowski, S., Kochanek, C. S., Udalski, A., Wyrzykowski, Ł., Soszyński, I., Szymański, M. K., et al. (2010b). Quantifying quasar variability as part of a general approach to classifying continuously varying sources. Astrophys. J. 708, 927–945. doi: 10.1088/0004-637X/708/2/927
Kozłowski, S., Onken, C. A., Kochanek, C. S., Udalski, A., Szymański, M. K., Kubiak, M., et al. (2013). The magellanic quasars survey. III. Spectroscopic confirmation of 758 active galactic nuclei behind the magellanic clouds. Astrophys. J. 775:92. doi: 10.1088/0004-637X/775/2/92
MacLeod, C. L., Ivezić, Ž., Kochanek, C. S., Kozłowski, S., Kelly, B., Bullock, E., et al. (2010). Modeling the time variability of SDSS stripe 82 quasars as a damped random walk. Astrophys. J. 721, 1014–1033. doi: 10.1088/0004-637X/721/2/1014
Vanden Berk, D. E., Wilhite, B. C., Kron, R. G., Anderson, S. F., Brunner, R. J., Hall, P. B., et al. (2004). The ensemble photometric variability of ~25,000 quasars in the sloan digital sky survey. Astrophys. J. 601, 692–714. doi: 10.1086/380563
Keywords: accretion, accretion disks, galaxies: active, methods: data analysis, quasars: general
Citation: Kozłowski S (2017) Optical Variability of Active Galactic Nuclei. Front. Astron. Space Sci. 4:14. doi: 10.3389/fspas.2017.00014
Received: 04 August 2017; Accepted: 31 August 2017;
Published: 21 September 2017.
Edited by:Mauro D'Onofrio, University of Padua, Italy
Reviewed by:Maurizio Paolillo, Dipartimento di Fisica Ettore Pancini, Università degli Studi di Napoli Federico II, Italy
Copyright © 2017 Kozłowski. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Szymon Kozłowski, email@example.com