Modelling variability power spectra of active galaxies from irregular time series

Abstract

A common feature of Active Galactic Nuclei (AGN) is their random variations in brightness across the whole emission spectrum, from radio to $\gamma$-rays. Studying the nature and origin of these fluctuations is critical to characterising the underlying variability process of the accretion flow that powers AGN. Random timing fluctuations are often studied with the power spectrum; this quantifies how the amplitude of variations is distributed over temporal frequencies. Red noise variability -- when the power spectrum increases smoothly towards low frequencies -- is ubiquitous in AGN. The commonly used Fourier analysis methods, have significant challenges when applied to arbitrarily sampled light curves of red noise variability. Several time-domain methods exist to infer the power spectral shape in the case of irregular sampling but they suffer from biases which can be difficult to mitigate, or are computationally expensive. In this paper, we demonstrate a method infer the shape of broad-band power spectra for irregular time series, using a Gaussian process regression method scalable to large datasets. The power spectrum is modelled as a power-law model with one or two bends with flexible slopes. The method is fully Bayesian and we demonstrate its utility using simulated light curves. Finally, Ark 564, a well-known variable Seyfert 1 galaxy, is used as a test case and we find consistent results with the literature using independent X-ray data from XMM-Newton and Swift. We provide publicly available, documented and tested implementations in Python and Julia.

21 pages, 16 figures, accepted for publication in MNRAS, code available https://github.com/mlefkir/Pioran.jl

References in corpus (49)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole
• LSST: from Science Drivers to Reference Design and Anticipated Data Products
• Observational Evidence of AGN Feedback
• The Swift X-ray Telescope
• First Sagittarius A* Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole in the Center of the Milky Way
• The No-U-Turn Sampler: Adaptively Setting Path Lengths in Hamiltonian Monte Carlo
• Methods and results of an automatic analysis of a complete sample of Swift-XRT observations of GRBs
• Rank-normalization, folding, and localization: An improved $\hat{R}$ for assessing convergence of MCMC
• On characterising the variability properties of X-ray light curves from active galaxies
• Fast and scalable Gaussian process modeling with applications to astronomical time series
• An online repository of Swift/XRT light curves of GRBs
• Are the Variations in Quasar Optical Flux Driven by Thermal Fluctuations?
• Statistics: Handle with Care, Detecting Multiple Model Components with the Likelihood Ratio Test
• Active Galactic Nuclei as scaled-up Galactic black holes
• Active Galactic Nuclei: what's in a name?
• Non-linear X-ray variability in X-ray binaries and active galaxies
• Measuring the broadband power spectra of active galactic nuclei with RXTE
• Investigating a Fluctuating-accretion Model for the Spectral-timing Properties of Accreting Black Hole Systems
• Combined Long and Short Timescale X-ray Variability of NGC4051 with RXTE and XMM-Newton
• Juliet: a versatile modelling tool for transiting and non-transiting exoplanetary systems
• Is Quasar Optical Variability a Damped Random Walk?
• The flux-dependent amplitude of broadband noise variability in X-ray binaries and active galaxies
• False periodicities in quasar time-domain surveys
• Stingray: A Modern Python Library For Spectral Timing
• Flexible and Scalable Methods for Quantifying Stochastic Variability in the Era of Massive Time-Domain Astronomical Data Sets
• X-ray Fluctuation Power Spectral Densities of Seyfert 1 Galaxies
• X-ray variability of 104 active galactic nuclei. XMM-Newton power-spectrum density profiles
• Collaborative Nested Sampling: Big Data vs. complex physical models
• Validating Bayesian Inference Algorithms with Simulation-Based Calibration
• A Bayesian test for periodic signals in red noise
• Discovery of multiple Lorentzian components in the X-ray timing properties of the Narrow Line Seyfert 1 Ark 564
• A Stochastic Model for the Luminosity Fluctuations of Accreting Black Holes
• Bayesian Workflow
• Long-Term X-ray Spectral Variability in Seyfert 1 Galaxies
• Spectral variability and reverberation time delays in the Suzaku X-ray spectrum of NGC 4051
• Discovery of a ~5 day characteristic timescale in the Kepler power spectrum of Zw 229-15
• Efficient modeling of correlated noise II. A flexible noise model with fast and scalable methods
• A Degeneracy in DRW Modelling of AGN Light Curves
• Calculating Time Lags From Unevenly-Sampled Light Curves
• Studies in Astronomical Time Series Analysis: VII. An Enquiry Concerning Non-Linearity, the RMS-Mean Flux Relation, and log-Normal Flux Distributions
• Low frequency X-ray timing with Gaussian processes and reverberation in the radio-loud AGN 3C 120
• A Quasi-periodic Oscillation in the gamma-ray Emission from the Non-blazar Active Galactic Nucleus PKS 0521-36
• Timing evidence in determining the accretion state of the Seyfert galaxy NGC 3783
• The X-ray Power Spectral Density Function of the Seyfert Active Galactic Nucleus NGC 7469
• Deep State-Space Gaussian Processes
• Parameter Estimation and Confidence Regions in the Method of Light Curve Simulations for the Analysis of Power Density Spectra
• Search for X-ray Quasi-periodicity of six AGNs using Gaussian Process method

Cited by in corpus (3)

• Continuum optical-UV and X-ray variability of AGN: current results and future challenges
• Scalable and Robust Multiband Modeling of AGN Light Curves in Rubin-LSST
• The ensemble broad-frequency power spectrum of Stripe-82 quasars from multiple surveys