Stochastic approach to modeling the γ-ray variability of Fermi/LAT blazars

We studied the gamma-ray variability of 13 blazars observed with the Fermi Large Area Telescope (LAT). These blazars were among the brightest ones monitored during the first 4 years of the Fermi sky survey. We modelled their gamma-ray light curves with the Ornstein-Uhlenbeck (OU) process or mixed OU process. The power spectral density (PSD) of the OU process is a zero-centered Lorentzian function, proportional to 1/f(alpha) with a changing at a characteristic time scale, tau(0), from 0 (tau >> tau(0)) to 2 (tau << tau(0)). The PSD of the mixed OU process has in addition an intermediate part with 0 < alpha < 2 between the long and short characteristic time scales. We show that the OU model provides a good description of the Fermi/LAT light curves of three blazars in our sample. For the first time we provide constraints on the characteristic gamma-ray time scale of variability in two BL Lac sources, 3C 66A and PKS 2155-304. We find that the mixed OU process describes the light curves of the remaining 10 blazars better than the OU process. We infer that their Fermi/LAT PSD resemble power-law functions and constrain their PSD slopes.