skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Characterization of optical light curves of extreme variability quasars over a ∼16-yr baseline
ABSTRACT We study the optical light curves – primarily probing the variable emission from the accretion disc – of ∼900 extreme variability quasars (EVQs, with maximum flux variations more than 1 mag) over an observed-frame baseline of ∼16 yr using public data from the SDSS Stripe 82, PanSTARRS-1 and the Dark Energy Survey. We classify the multiyear long-term light curves of EVQs into three categories roughly in the order of decreasing smoothness: monotonic decreasing or increasing (3.7 per cent), single broad peak and dip (56.8 per cent), and more complex patterns (39.5 per cent). The rareness of monotonic cases suggests that the major mechanisms driving the extreme optical variability do not operate over time-scales much longer than a few years. Simulated light curves with a damped random walk model generally under-predict the first two categories with smoother long-term trends. Despite the different long-term behaviours of these EVQs, there is little dependence of the long-term trend on the physical properties of quasars, such as their luminosity, BH mass, and Eddington ratio. The large dynamic range of optical flux variability over multiyear time-scales of these EVQs allows us to explore the ensemble correlation between the short-term (≲6 months) variability and the seasonal-average flux across the decade-long baseline (the rms-mean flux relation). We find that unlike the results for X-ray variability studies, the linear short-term flux variations do not scale with the seasonal-average flux, indicating different mechanisms that drive the short-term flickering and long-term extreme variability of accretion disc emission. Finally, we present a sample of 16 EVQs, where the approximately bell-shaped large amplitude variation in the light curve can be reasonably well fit by a simple microlensing model.  more » « less
Award ID(s):
1715579
PAR ID:
10161999
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
494
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
3686 to 3698
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT The variability of quasars across multiple wavelengths is a useful probe of physical conditions in active galactic nuclei. In particular, variable accretion rates, instabilities, and reverberation effects in the accretion disc of a supermassive black hole are expected to produce correlated flux variations in ultraviolet (UV) and optical bands. Recent work has further argued that binary quasars should exhibit strongly correlated UV and optical periodicities. Strong UV–optical correlations have indeed been established in small samples of (N ≲ 30) quasars with well-sampled light curves, and have extended the ‘bluer-when-brighter’ trend previously found within the optical bands. Here, we further test the nature of quasar variability by examining the observed-frame UV–optical correlations among bright quasars extracted from the Half Million Quasars (HMQ) catalogue. We identified a large sample of 1315 quasars in HMQ with overlapping UV and optical light curves from the Galaxy Evolution Explorer and the Catalina Real-time Transient Survey, respectively. We find that strong correlations exist in this much larger sample, but we rule out, at ∼95 per cent confidence, the simple hypothesis that the intrinsic UV and optical variations of all quasars are fully correlated. Our results therefore imply the existence of physical mechanism(s) that can generate uncorrelated optical and UV flux variations. 
    more » « less
  2. ABSTRACT We study the optical gri photometric variability of a sample of 190 quasars within the SDSS Stripe 82 region that have long-term photometric coverage during ∼1998−2020 with SDSS, PanSTARRS-1, the Dark Energy Survey, and dedicated follow-up monitoring with Blanco 4m/DECam. With on average ∼200 nightly epochs per quasar per filter band, we improve the parameter constraints from a Damped Random Walk (DRW) model fit to the light curves over previous studies with 10–15 yr baselines and ≲ 100 epochs. We find that the average damping time-scale τDRW continues to rise with increased baseline, reaching a median value of ∼750 d (g band) in the rest frame of these quasars using the 20-yr light curves. Some quasars may have gradual, long-term trends in their light curves, suggesting that either the DRW fit requires very long baselines to converge, or that the underlying variability is more complex than a single DRW process for these quasars. Using a subset of quasars with better-constrained τDRW (less than 20 per cent of the baseline), we confirm a weak wavelength dependence of τDRW∝λ0.51 ± 0.20. We further quantify optical variability of these quasars over days to decades time-scales using structure function (SF) and power spectrum density (PSD) analyses. The SF and PSD measurements qualitatively confirm the measured (hundreds of days) damping time-scales from the DRW fits. However, the ensemble PSD is steeper than that of a DRW on time-scales less than ∼ a month for these luminous quasars, and this second break point correlates with the longer DRW damping time-scale. 
    more » « less
  3. ABSTRACT NK22 proposed a new method to reconstruct the temperature perturbation map (as functions of time and disc radius) of active galactic nuclei (AGN) accretion discs using multiwavelength photometric light curves. We apply their technique to 100 quasars at z = 0.5–2 from the Sloan Digital Sky Survey Reverberation Mapping project, using multi-epoch spectroscopy that covers rest-frame UV-optical continuum emission from the quasar and probes days to months time-scales. Consistent with NK22 for low-redshift AGNs, we find that the dominant pattern of disc temperature perturbations is either slow inward/outward moving waves with typical amplitudes $$\delta T/T_0\sim 10~{{\ \rm per \, cent}}$$ traveling at ∼0.01–0.1c, with a typical radial frequency of ∼ 0.5 dex in log R, or incoherent perturbations. In nearly none of the cases do we find clear evidence for coherent, fast outgoing temperature perturbations at the speed of light, reminiscent of the lamppost model; but such lamppost signals may be present in some quasars for limited periods of the monitoring data. Using simulated data, we demonstrate that high-fidelity temperature perturbation maps can be recovered with high-quality monitoring spectroscopy, with limited impact from seasonal gaps in the data. On the other hand, reasonable temperature perturbation maps can be reconstructed with high-cadence photometric light curves from the Vera C Rubin Observatory Legacy Survey of Space and Time. Our findings, together with NK22, suggest that internal disc processes are the main driver for temperature fluctuations in AGN accretion discs over days to months time-scales. 
    more » « less
  4. ABSTRACT The scarce optical variability studies in spectrally classified Type 2 active galactic nuclei (AGNs) have led to the discovery of anomalous objects that are incompatible with the simplest unified models (UMs). This paper focuses on the exploration of different variability features that allow to distinguish between obscured, Type 2 AGNs and the variable, unobscured Type 1s. We analyse systematically the Zwicky Transient Facility, 2.5-yr-long light curves of ∼15 000 AGNs from the Sloan Digital Sky Survey Data Release 16, which are generally considered Type 2s due to the absence of strong broad emission lines (BELs). Consistent with the expectations from the UM, the variability features are distributed differently for distinct populations, with spectrally classified weak Type 1s showing one order of magnitude larger variances than the Type 2s. We find that the parameters given by the damped random walk model lead to broader H α equivalent width for objects with τg > 16 d and long-term structure function SF∞, g > 0.07 mag. By limiting the variability features, we find that ∼11 per cent of Type 2 sources show evidence for optical variations. A detailed spectral analysis of the most variable sources (∼1 per cent of the Type 2 sample) leads to the discovery of misclassified Type 1s with weak BELs and changing-state candidates. This work presents one of the largest systematic investigations of Type 2 AGN optical variability to date, in preparation for future large photometric surveys. 
    more » « less
  5. ABSTRACT Sgr A* exhibits regular variability in its multiwavelength emission, including daily X-ray flares and roughly continuous near-infrared (NIR) flickering. The origin of this variability is still ambiguous since both inverse Compton and synchrotron emission are possible radiative mechanisms. The underlying particle distributions are also not well constrained, particularly the non-thermal contribution. In this work, we employ the GPU-accelerated general relativistic magnetohydrodynamics code H-AMR to perform a study of flare flux distributions, including the effect of particle acceleration for the first time in high-resolution 3D simulations of Sgr A*. For the particle acceleration, we use the general relativistic ray-tracing code bhoss to perform the radiative transfer, assuming a hybrid thermal+non-thermal electron energy distribution. We extract ∼60 h light curves in the sub-millimetre, NIR and X-ray wavebands, and compare the power spectra and the cumulative flux distributions of the light curves to statistical descriptions for Sgr A* flares. Our results indicate that non-thermal populations of electrons arising from turbulence-driven reconnection in weakly magnetized accretion flows lead to moderate NIR and X-ray flares and reasonably describe the X-ray flux distribution while fulfilling multiwavelength flux constraints. These models exhibit high rms per cent amplitudes, $$\gtrsim 150{{\ \rm per\ cent}}$$ both in the NIR and the X-rays, with changes in the accretion rate driving the 230 GHz flux variability, in agreement with Sgr A* observations. 
    more » « less