Probing AGN Accretion Physics through AGN Variability

Active Galactic Nuclei (AGN) exhibit large luminosity variations over the entire electromagnetic spectrum on timescales ranging from hours to years. The variations in luminosity are devoid of any periodic character and appear stochastic. While complex correlations exist between the variability observed in different parts of the electromagnetic spectrum, no frequency band appears to be completely dominant, suggesting that the physical processes producing the variability are exceedingly rich and complex. In the absence of a clear theoretical explanation of the variability, phenomenological models are used to study AGN variability. The stochastic behavior of AGN variability makes formulating such models difficult and connecting them to the underlying physics exceedingly hard. We study AGN light curves serendipitously observed by the NASA Kepler planet-finding mission. Compared to previous ground-based observations, Kepler offers higher precision and a smaller sampling interval resulting in potentially higher quality light curves. Using structure functions, we demonstrate that (1) the simplest statistical model of AGN variability, the damped random walk (DRW), is insufficient to characterize the observed behavior of AGN light curves; and (2) variability begins to occur in AGN on time-scales as short as hours. Of the 20 light curves studied by us, only 3-8 may be consistent with the DRW. The structure functions of the AGN in our sample exhibit complex behavior with pronounced dips on time-scales of 10-100 d suggesting that AGN variability can be very complex and merits further analysis. We examine the accuracy of the Kepler pipeline-generated light curves and find that the publicly available light curves may require re-processing to reduce contamination from field sources. We show that while the re-processing changes the exact PSD power law slopes inferred by us, it is unlikely to change the conclusion of our structure function study-Kepler AGN light curves indicate that the DRW is insufficient to characterize AGN variability. We provide a new approach to probing accretion physics with variability by decomposing observed light curves into a set of impulses that drive diffusive processes using C-ARMA models. Applying our approach to Kepler data, we demonstrate how the time-scales reported in the literature can be interpreted in the context of the growth and decay time-scales for flux perturbations and tentatively identify the flux perturbation driving process with accretion disk turbulence on length-scales much longer than the characteristic eddy size. Our analysis technique is applicable to (1) studying the connection between AGN sub-type and variability properties; (2) probing the origins of variability by studying the multi-wavelength behavior of AGN; (3) testing numerical simulations of accretion flows with the goal of creating a library of the variability properties of different accretion mechanisms; (4) hunting for changes in the behavior of the accretion flow by block-analyzing observed light curves; and (5) constraining the sampling requirements of future surveys of AGN variability.