Assessment of stochastic and deterministic models of 6304 quasar lightcurves from SDSS Stripe 82

The optical light curves of many quasars show variations of tenths of a magnitude or more on time scales of months to years. This variation often cannot be described well by a simple deterministic model. We perform a Bayesian comparison of over 20 deterministic and stochastic models on 6304 QSO light curves in SDSS Stripe 82. We include the damped random walk (or Ornstein-Uhlenbeck [OU] process), a particular type of stochastic model which recent studies have focused on. Further models we consider are single and double sinusoids, multiple OU processes, higher order continuous autoregressive processes, and composite models. We find that only 29 out of 6304 QSO lightcurves are described significantly better by a deterministic model than a stochastic one. The OU process is an adequate description of the vast majority of cases (6023). Indeed, the OU process is the best single model for 3462 light curves, with the composite OU process/sinusoid model being the best in 1706 cases. The latter model is the dominant one for brighter/bluer QSOs. Furthermore, a non-negligible fraction of QSO lightcurves show evidence that not only the mean is stochastic but the variance is stochastic, too. Our results confirm earlier work that QSO light curves can be described with a stochastic model, but place this on a firmer footing, and further show that the OU process is preferred over several other stochastic and deterministic models. Of course, there may well exist yet better (deterministic or stochastic) models which have not been considered here.

[1]  Ronald Christensen,et al.  Testing Fisher, Neyman, Pearson, and Bayes , 2005 .

[2]  A. J. Drake,et al.  Variability-selected Quasars in MACHO Project Magellanic Cloud Fields , 2002, astro-ph/0209513.

[3]  B. Peterson,et al.  Characteristic Ultraviolet/Optical Timescales in Active Galactic Nuclei , 2001 .

[4]  M. Rees.  BLACK HOLE MODELS FOR ACTIVE GALACTIC NUCLEI , 1984 .

[5]  D. Gillespie.  The mathematics of Brownian motion and Johnson noise , 1996 .

[6]  Zeljko Ivezic,et al.  Sloan Digital Sky Survey Standard Star Catalog for Stripe 82: The Dawn of Industrial 1% Optical Photometry , 2007, astro-ph/0703157.

[7]  R. Engle.  Autoregressive conditional heteroscedasticity with estimates of the variance of United Kingdom inflation , 1982 .

[8]  J. Berger.  Could Fisher, Jeffreys and Neyman Have Agreed on Testing? , 2003 .

[9]  Brandon C. Kelly,et al.  Are the Variations in Quasar Optical Flux Driven by Thermal Fluctuations , 2009 .

[10]  Joseph C. Shields,et al.  THE COLOR VARIABILITY OF QUASARS , 2011, 1109.6653.

[11]  Gregory Dobler,et al.  SELECTING QUASARS BY THEIR INTRINSIC VARIABILITY , 2010, 1002.2642.

[12]  R. Kron,et al.  Continuum Variability of Active Galactic Nuclei in the Optical-Ultraviolet Range , 2000, astro-ph/0012408.

[13]  Usa,et al.  QUANTIFYING QUASAR VARIABILITY AS PART OF A GENERAL APPROACH TO CLASSIFYING CONTINUOUSLY VARYING SOURCES , 2009, 0909.1326.

[14]  Nathaniel R. Butler,et al.  OPTIMAL TIME-SERIES SELECTION OF QUASARS , 2010, 1008.3143.

[15]  Pavlos Protopapas,et al.  A Refined QSO Selection Method Using Diagnostics , 2011, Proceedings of the International Astronomical Union.

[16]  S. Baliunas,et al.  A Prescription for period analysis of unevenly sampled time series , 1986 .

[17]  E. Bullock,et al.  MODELING THE TIME VARIABILITY OF SDSS STRIPE 82 QUASARS AS A DAMPED RANDOM WALK , 2010, 1004.0276.

[18]  S. Mineshige,et al.  Optical Variability in Active Galactic Nuclei: Starbursts or Disk Instabilities? , 1997, astro-ph/9712006.

[19]  Paul S. Smith,et al.  Long-term optical variability properties of the Palomar–Green quasars , 1999 .

[20]  John A. Nelder,et al.  A Simplex Method for Function Minimization , 1965, Comput. J..

[21]  Richard A. Davis,et al.  Introduction to time series and forecasting , 1998 .

[22]  P. Protopapas,et al.  An improved quasar detection method in EROS-2 and MACHO LMC data sets , 2012, 1304.0401.

[23]  C. Bailer-Jones.  A Bayesian method for the analysis of deterministic and stochastic time series , 2012, 1209.3730.

[24]  C. Kochanek,et al.  IS QUASAR OPTICAL VARIABILITY A DAMPED RANDOM WALK? , 2012, 1202.3783.

[25]  T. Bollerslev,et al.  Generalized autoregressive conditional heteroskedasticity , 1986 .

[26]  Brandon C. Kelly,et al.  A STOCHASTIC MODEL FOR THE LUMINOSITY FLUCTUATIONS OF ACCRETING BLACK HOLES , 2010, 1009.6011.

[27]  G. Uhlenbeck,et al.  On the Theory of the Brownian Motion , 1930 .

[28]  Pavlos Protopapas,et al.  QSO Selection Algorithm Using Time Variability and Machine Learning: Selection of 1,620 QSO Candidates from MACHO LMC Database , 2011, 1101.3316.

[29]  Mamoru Doi,et al.  Exploring the Variable Sky with the Sloan Digital Sky Survey , 2007, 0704.0655.

[30]  Pavlos Protopapas,et al.  QUASI-STELLAR OBJECT SELECTION ALGORITHM USING TIME VARIABILITY AND MACHINE LEARNING: SELECTION OF 1620 QUASI-STELLAR OBJECT CANDIDATES FROM MACHO LARGE MAGELLANIC CLOUD DATABASE , 2011 .

[31]  Pavlos Protopapas,et al.  A REFINED QSO SELECTION METHOD USING DIAGNOSTICS TESTS: 663 QSO CANDIDATES IN THE LARGE MAGELLANIC CLOUD , 2011, 1110.5632.