Quasar Variability Research Articles

Probing the Origin of Changing-look Quasar Transitions with Chandra

Extremely variable quasars can also show strong changes in broad-line emission strength and are known as changing-look quasars (CLQs). To study the CLQ transition mechanism, we present a pilot sample of CLQs with X-ray observations in both the bright and faint states. From a sample of quasars with bright-state archival SDSS spectra and (Chandra or XMM-Newton) X-ray data, we identified five new CLQs via optical spectroscopic follow-up and then obtained new target-of-opportunity X-ray observations with Chandra. No strong absorption is detected in either the bright- or the faint-state X-ray spectra. The intrinsic X-ray flux generally changes along with the optical variability, and the X-ray power-law slope becomes harder in the faint state. Large-amplitude mid-infrared variability is detected in all five CLQs, and it echoes the variability in the optical with a time lag expected from the light-crossing time of the dusty torus for CLQs with robust lag measurements. The changing-obscuration model is not consistent with the observed X-ray spectra and spectral energy distribution changes seen in these CLQs. It is highly likely that the observed changes are due to the changing accretion rate of the supermassive black hole, so the multiwavelength emission varies accordingly, with promising analogies to the accretion states of X-ray binaries.

The variational slope of quasar light curves is not a distance indicator

ABSTRACT When the time difference quotients, or variational slopes, of quasar light curves are plotted against their absolute magnitudes, there is a tight positive correlation of ∼0.16 dex in the variational slope direction or ∼0.5 dex in the absolute magnitude direction. This finding resulted in suggestions that a variational slope–luminosity relation could be used as a distance indicator. However, I show that this relation can be explained almost entirely from self-correlation with luminosity. After properly accounting for the self-correlation component, the relation has a true scatter of ∼1.5 dex in luminosity, consistent with established correlations for quasar variability amplitudes. Given this large scatter, correlation with variational slope or variability amplitude and luminosity is not by itself a suitable distance indicator for quasars.

The need for obscured supermassive black hole growth to explain quasar proximity zones in the epoch of reionization

ABSTRACTProximity zones of high-redshift quasars are unique probes of supermassive black hole formation, but simultaneously explaining proximity zone sizes and black hole masses has proved to be challenging. We study the robustness of some of the assumptions that are usually made to infer quasar lifetimes from proximity zone sizes. We show that small proximity zones can be readily explained by quasars that vary in brightness with a short duty cycle of fduty ∼ 0.1 and short bright periods of ton ∼ 104 yr, even for long lifetimes. We further show that reconciling this with black hole mass estimates requires the black hole to continue to grow and accrete during its obscured phase. The consequent obscured fractions of ≳ 0.7 or higher are consistent with low-redshift measurements and models of black hole accretion. Such short duty cycles and long obscured phases are also consistent with observations of large proximity zones, thus providing a simple, unified model for proximity zones of all sizes. The large dynamic range of our simulation, and its calibration to the Lyman-α forest, allows us to investigate the influence of the large-scale topology of reionization and the quasar’s host halo mass on proximity zones. We find that incomplete reionization can impede the growth of proximity zones and make them smaller up to 30 per cent, but the quasar host halo mass only affects proximity zones weakly and indirectly. Our work suggests that high-redshift proximity zones can be an effective tool to study quasar variability and black hole growth.

Correction to: Optical variability of quasars with 20-year photometric light curves

Correction to: Optical variability of quasars with 20-year photometric light curves

Optical Variability of Quasars with 20-Year Photometric Light Curves

Optical Variability of Quasars with 20-Year Photometric Light Curves

UnTimely: a Full-sky, Time-domain unWISE Catalog

We present the unTimely Catalog, a deep time-domain catalog of detections based on Wide-field Infrared Survey Explorer (WISE) and NEOWISE observations spanning the 2010 through 2020 time period. Detections are extracted from “time-resolved unWISE coadds,” which stack together each biannual sky pass of WISE imaging to create a set of ∼16 all-sky maps (per band), each much deeper and cleaner than individual WISE exposures. unTimely incorporates the W1 (3.4 μm) and W2 (4.6 μm) channels, meaning that our data set effectively consists of ∼32 full-sky unWISE catalogs. We run the crowdsource crowded-field point-source photometry pipeline (Schlafly et al. 2018) on each epochal coadd independently, with low detection thresholds: S/N = 4.0 (2.5) in W1 (W2). In total, we tabulate and publicly release 23.5 billion (19.9 billion) detections at W1 (W2). unTimely is ∼1.3 mag deeper than the WISE/NEOWISE Single Exposure Source Tables near the ecliptic, with further enhanced depth toward higher ecliptic latitudes. The unTimely Catalog is primarily designed to enable novel searches for faint, fast-moving objects, such as Y dwarfs and/or late-type (T/Y) subdwarfs in the Milky Way’s thick disk or halo. unTimely will also facilitate other time-domain science applications, such as all-sky studies of quasar variability at mid-infrared wavelengths over a decade-long time baseline.

The optical variability analysis of quasars based on fractal Brownian motion in Sloan digital sky survey

<p indent="0mm">The optical variability analysis of quasars based on fractal Brownian motion in the Sloan Digital Sky Survey (SDSS). The violent activities of the high energy variability can be observed in multiple wavelength bands of quasars, and studying quasar variability can contribute to a deeper understanding of the physical mechanism of the emission from the central region. Many studies have shown that the great majority of quasar variability is considered to originate in the accretion disk’s thermal emission, but some people do not agree with this view. The Hurst exponent of SDSS quasar variability is calculated using the rescaled range analysis method to investigate the fractal properties of variability in quasars from the SDSS stripe 82. The Hurst exponent results of variability show that the quasar variability is self-similar and long-range correlated. The variability obeys the biased random walk model, and its variation is similar to fractal Brownian motion. The power spectral exponent obtained from the Hurst exponent conversion shows that most of the variabilities can be approximately fitted by the damped random walk model, but a small number of samples have large deviations.

Quasars with Periodic Variability: Capabilities and Limitations of Bayesian Searches for Supermassive Black Hole Binaries in Time-domain Surveys

Supermassive black hole binaries (SMBHBs) are an inevitable consequence of galaxy mergers. At sub-parsec separations, they are practically impossible to resolve, and the most promising technique is to search for quasars with periodic variability. However, searches for quasar periodicity in time-domain data are challenging due to the stochastic variability of quasars. In this paper, we used Bayesian methods to disentangle periodic SMBHB signals from intrinsic damped random walk (DRW) variability in active galactic nuclei light curves. We simulated a wide variety of realistic DRW and DRW+sine light curves. Their observed properties are modeled after the Catalina Real-time Transient Survey (CRTS) and expected properties of the upcoming Legacy Survey of Space and Time (LSST) from the Vera C. Rubin Observatory. Through a careful analysis of parameter estimation and Bayesian model selection, we investigated the range of parameter space for which binary systems can be detected. We also examined which DRW signals can mimic periodicity and be falsely classified as binary candidates. We found that periodic signals are more easily detectable if the period is short or the amplitude of the signal is large compared to the contribution of the DRW noise. We saw similar detection rates both in the CRTS and LSST-like simulations, while the false-detection rate depends on the quality of the data and is minimal in LSST. Our idealized simulations provide an excellent way to uncover the intrinsic limitations in quasar periodicity searches and set the stage for future searches for SMBHBs.

GPS Measurements of Precipitable Water Vapor Can Improve Survey Calibration: A Demonstration from KPNO and the Mayall z-band Legacy Survey

Dual-band Global Positioning Satellite (GPS) measurements of precipitable water vapor (PWV) at the Kitt Peak National Observatory predict the overall per-image sensitivity of the Mayall z-band Legacy Survey (MzLS). The per-image variation in the brightness of individual stars is strongly correlated with the measured PWV and the color of the star. Synthetic stellar spectra through TAPAS transmission models successfully predict the observed PWV-induced photometric variation. We find that PWV absorption can be well approximated by a linear relationship with (airmass × PWV)0.6 and present an update on the traditional treatment in the literature. The MzLS zero-point sensitivity in electrons s−1 varies with a normalized-mean absolute deviation of 61 mmag. PWV variation accounts 23 mmag of this zero-point variation. The MzLS per-image absolute sensitivity decreases by 40 mmag per effective mm of PWV. The overall gray offset portion of this variation is corrected by the calibration to a reference catalog. But the relative calibration error between blue (r − z < 0.5 mag) versus red (1.2 mag < r − z) stars increases by 0.3–2 mmag per effective mm of PWV. We argue that GPS systems provide more precise PWV measurements than using differential measurements of stars of different colors and recommend that observatories install dual-band GPS as a low-maintenance, low-cost, auxiliary calibration system. We extend our results of the need for well-calibrated PWV measurements by presenting the calculations of the PWV photometric impact on three science cases of interest: stellar photometry, supernova cosmology, and quasar identification and variability.

Optical variability of quasars with 20-yr photometric light curves

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.

Cosmological nanolensing by dense gas clouds

ABSTRACT We study the influence of a cosmological population of dense gas clouds on distant sources, with an emphasis on quasar optical variability. In addition to gravitational lensing, such clouds affect flux measurements via refraction in the neutral gas and via dust extinction, leading to a variety of possible light curves even in the low optical depth limit. We classify and illustrate the types of light curves that can arise. For sources as large as quasars, we show that gravitational lensing and extinction are the dominant effects, with gas refraction playing only a minor role. We find that clouds with mass $\sim 10^{-4.5\pm 0.5}\, \mathrm{M}_\odot$ can reproduce the observed distribution of quasar variation amplitudes, but only if such clouds make up a large fraction of the closure density. In that case, there may also be substantial extinction of distant optical sources, which can, in principle, be constrained by data on ‘standard candles’ such as type Ia supernovae. Unfortunately, that extinction is essentially grey, even when the material opacity is strongly wavelength dependent, making it difficult to distinguish from the influence of the background geometry. We propose a novel statistical test of the origin of quasar variability, based on the angular structure of the variation timescale for a large number of quasars distributed all over the sky. If quasar variability is primarily due to nanolensing, then angular structure is expected to include a quadrupole term of amplitude $\sim 5{{\ \rm per\ cent}}$, which ought to be measurable with future data from the Gaia mission.

Variability in quasar light curves: using quasars as standard candles

A relation between the variational slope, sF , and the mean absolute magnitude, 〈M〉, in the light curves of 58 spectroscopically confirmed quasars is measured with a dispersion of 0.15dex.Assuming it holds for quasars in general, not only does this relation add to our working knowledge of quasar variability but it also shows great promise at accurately measuring luminosity distance to a quasar in a model independent way.An accurate, model independent measure of the luminosity distance would allow quasars to be added to the cosmic distance ladder, easily extending the ladder out far beyond the redshifts accessible to type Ia supernovae where cosmological parameters can be better constrained.

The relation between X-ray and ultraviolet variability of quasars

ABSTRACT The relation between X-ray and UV/optical variability in AGNs has been explored in many individual sources, however a large sample study is yet absent. Through matching the XMM-Newton serendipitous X-ray and UV source catalogs with SDSS quasars, we build a sample of 802 epoch-pairs of 525 quasars showing clear variability in logFX − logFUV space. After correcting for the effect of photometric noise, we find 35.6 ± 2.1 per cent of the epoch-pairs show asynchronous variability between X-ray and UV (brightening in one band but dimming in the other). This indicates only in 28.8 ± 4.2 per cent of the epoch-pairs the X-ray and UV variability are intrinsically coordinated. The variability synchronicity exhibits no dependence on physical parameters of quasars or the time lag of the epoch-pairs, except for stronger variability tends to have stronger synchronicity. Switches between synchronous and asynchronous variability are also seen in individual sources. The poor coordination clearly contradicts both the X-ray reprocessing model and the accretion rate variation model for AGN variability. The ratios of the observed X-ray variability amplitude to that in UV span a broad range and peak at ∼2. The dominant fraction of the ratios appear too small to be attributed to X-ray reprocessing, and too large for accretion rate variation. The inhomogeneous disc model which incorporates both X-ray and UV/optical variability in AGNs is favoured by the observed stochastic relation between X-ray and UV variations.

Optical Variability of Quasars with 20-Year Photometric Light Curves

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 $\sim 1998-2020$ with SDSS, PanSTARRS-1, the Dark Energy Survey, and dedicated follow-up monitoring with Blanco 4m/DECam. With on average $\sim 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 $\lesssim 100$ epochs. We find that the average damping timescale $\tau_{\rm DRW}$ continues to rise with increased baseline, reaching a median value of $\sim 750$ days ($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 $\tau_{\rm DRW}$ (less than 20\% of the baseline), we confirm a weak wavelength dependence of $\tau_{\rm DRW}\propto \lambda^{0.51\pm0.20}$. We further quantify optical variability of these quasars over days to decades timescales using structure function (SF) and power spectrum density (PSD) analyses. The SF and PSD measurements qualitatively confirm the measured (hundreds of days) damping timescales from the DRW fits. However, the ensemble PSD is steeper than that of a DRW on timescales less than $\sim$ a month for these luminous quasars, and this second break point correlates with the longer DRW damping timescale.

Extreme Variability Quasars in Their Various States. I. The Sample Selection and Composite SDSS Spectra

Abstract Extremely variable quasars (EVQs) are a population of sources showing large optical photometric variability revealed by time-domain surveys. The physical origin of such extreme variability is yet unclear. In this first paper of a series, we construct the largest-ever sample of 14,012 EVQs using more than 15 yr of photometric data from Sloan Digital Sky Survey (SDSS) and Pan-STARRS1. We divide the EVQs into five subsamples according to the relative brightness of each EVQ during SDSS spectroscopic observation compared with the mean brightness from photometric observations. Corresponding control samples of normal quasars are built with matched redshift, bolometric luminosity, and supermassive black hole mass. We obtain the composite SDSS spectra of EVQs in various states and their corresponding control samples. We find EVQs exhibit clearly bluer SDSS spectra during bright states and clearly redder spectra during dim states, consistent with the “bluer-when-brighter” trend widely seen in normal quasars. We further find that the line equivalent widths (EWs) of broad Mg II, C IV and [O III] (but not broad Hβ, which is yet puzzling) gradually decreases from the dim state to the bright state, similar to the so-called intrinsic Baldwin effect commonly seen in normal active galactic nuclei. In addition, EVQs have consistently larger line EWs compared with the control samples. We also see that EVQs show slight excess in the very broad line component compared with control samples. Possible explanations for the discoveries are discussed. Our findings support the hypothesis that EVQs are in the tail of a broad distribution of quasar properties but are not a distinct population.

A study of the luminosity variability of distant quasars by near-infrared imaging with Nayuta/NIC

A study of the luminosity variability of distant quasars by near-infrared imaging with Nayuta/NIC

The Locus Algorithm: A novel technique for identifying optimised pointings for differential photometry

Studies of the photometric variability of astronomical sources from ground-based telescopes must overcome atmospheric extinction effects. Differential photometry by reference to an ensemble of reference stars which closely match the target in terms of magnitude and colour can mitigate these effects. This Paper describes the design, implementation, and operation of a novel algorithm – The Locus Algorithm – which enables optimised differential photometry. The Algorithm is intended to identify, for a given target and observational parameters, the Field of View (FoV) which includes the target and the maximum number of reference stars similar to the target. A collection of objects from a catalogue (e.g. SDSS) is filtered to identify candidate reference stars and determine a rating for each which quantifies its similarity to the target. The algorithm works by defining a locus of points around each candidate reference star, upon which the FoV can be centred and include the reference at the edge of the FoV. The Points of Intersection (PoI) between these loci are identified and a score for each PoI is calculated. The PoI with the highest score is output as the optimum pointing. The steps of the algorithm are precisely defined in this paper. The application of The Locus Algorithm to a sample target, SDSS1237680117417115655, from the Sloan Digital Sky Survey is described in detail. The algorithm has been defined here and implemented in software which is available online. The algorithm has also been used to generate catalogues of pointings to optimise Quasar variability studies and to generate catalogues of optimised pointings in the search for Exoplanets via the transit method.

Revealing the structure of the lensed quasar Q 0957+561

Aims. We aim to use signatures of microlensing induced by stars in the foreground lens galaxy to infer the size of the accretion disk in the gravitationally lensed quasar Q 0957+561. The long-term photometric monitoring of this system (which so far has provided the longest available light curves of a gravitational lens system) permits us to evaluate the impact of uncertainties on our recently developed method (controlled by the distance between the modeled and the experimental magnitude difference histograms between two lensed images), and thus to test the robustness of microlensing-based disk-size estimates. Methods. We analyzed the well-sampled 21-year GLENDAMA optical light curves of the double-lensed quasar and studied the intrinsic and extrinsic continuum variations. Using accurate measurements for the time delay between the images A and B, we modeled and removed the intrinsic quasar variability, and from the statistics of microlensing magnifications we used a Bayesian method to derive the size of the region emitting the continuum at λrest = 2558 Å. Results. Analysis of the Q 0957+561 R-band light curves show a slow but systematic increase in the brightness of the B relative to the A component during the past ten years. The relatively low strength of the magnitude differences between the images indicates that the quasar has an unusually big optical accretion disk of half-light radius: R1/2 = 17.6±6.1 √(M/0.3 M⊙) lt-days.

30 Years of Multifrequency Quasar Variability: A Personal Journey

I outline the history of and progress in observing and understanding quasar multi-frequency and multi-messenger variability from the point of view of someone who has been working in the field for over 30 years. I will present some important references for the evolution from optical monitoring to multi-frequency cooperative programs that revealed the true multi-frequency/multi-timescale nature of variability in these objects. Quasar observations began with separate radio and optical monitoring programs; then the optical and radio observations w ere combined. This was followed by expanding the analyses to include far IR, UV, X-rays, and finally adding gamma rays. This progression yielded simultaneous multi-frequency spectra of these objects and light curves over 15 decades in frequency. The future is adding particle (neutrino) and gravitational waves to the picture. I also present long-term (50 years) optical light curves, and discuss optical variability at all timescales from minutes to tens of years in some selected objects for which we have reliable long-term monitoring observations.

A Sample Bias in Quasar Variability Studies

Abstract When a flux-limited quasar sample is observed at later times, there will be more dimmed quasars than brightened ones, due to a selection bias induced at the time of sample selection. Quasars are continuously varying and there are more fainter quasars than brighter ones. At the time of selection, even symmetrical variability will result in more quasars with their instantaneous fluxes scattered above the flux limit than those scattered below, leading to an asymmetry in flux changes over time. The same bias would lead to an asymmetry in the ensemble structure function (SF) of the sample such that the SF based on pairs with increasing fluxes will be slightly smaller than that based on pairs with decreasing fluxes. We use simulated time-symmetric quasar light curves based on the damped random walk prescription to illustrate the effects of this bias. The level of this bias depends on the sample, the threshold of magnitude changes, and the coverage of light curves, but the general behaviors are consistent. In particular, the simulations matched to recent observational studies with decade-long light curves produce an asymmetry in the SF measurements at the few percent level, similar to the observed values. These results provide a cautionary note on the reported time asymmetry in some recent quasar variability studies.