The universal power spectrum of quasars in optical wavelengths

Abstract

Context. The optical variability of quasars is one of the few windows through which we can explore the behaviour of accretion discs around supermassive black holes. Aims. We aim to establish the dependence of variability properties, such as characteristic timescales and the variability amplitude, on basic quasar parameters such as black hole mass and the accretion rate, controlling for the rest-frame wavelength of emission. Methods. Using large catalogues of quasars, we selected the g-band light curves for 4770 objects from the Zwicky Transient Facility archive. All the selected objects fall into a narrow redshift bin, 0.6 < z < 0.7, but cover a wide range of accretion rates in Eddington units (REdd) and black hole masses (M). We grouped these objects into 26 independent bins according to these parameters, calculated low-resolution g-band variability power spectra for each of these bins, and approximated the power spectra with a simple analytic model that features a break at a timescale, tb. Results. We find a clear dependence of the break timescale, tb, on REdd, on top of the known dependence of tb on the black hole mass, M. In our fits, tb ∝ M0.65 − 0.55REdd0.35−0.3, where the ranges in the exponents correspond to the best-fitting parameters of different power spectrum models. This mass dependence is slightly steeper than that found in other studies. Scaling tb to the orbital timescale of the innermost stable circular orbit (ISCO), tISCO, results approximately in tb/tISCO ∝ (REdd/M)0.35. In the standard thin disc model, (REdd/M) ∝ Tmax4, where Tmax is the maximum disc temperature, so that tb/tISCO appears to scale approximately with the maximum temperature of the disc to a small power. The observed values of tb are ∼10 longer than the orbital timescale at the light-weighted average radius of the disc region emitting in the (observer frame) g-band. The different scaling of the break frequency with M and REdd shows that the shape of the variability power spectrum cannot be solely a function of the quasar luminosity, even for a single rest-frame wavelength. Finally, the best-fitting models have slopes above the break in the range between −2.5 and −3. A slope of −2, as in the damped random walk models, fits the data significantly worse.