2021MNRAS.505.6103P

Decade-long monitoring of blazars at optical and infrared (OIR) wavelengths with the Small and Moderate Aperture Research Telescope System (SMARTS) in Chile and in γ-rays with the Fermi -Large Area Telescope (LAT) has enabled the systematic study of their multiwavelength long-term variability. In this work, we investigate, from a theoretical perspective, the long-term variability properties of blazar emission by introducing an observationally motivated time-dependence to four main parameters of the one-zone leptonic model: injection luminosity of relativistic electrons, strength of magnetic field, Doppler factor, and external photon field luminosity. For the first time, we use both the probability density function and the power spectral density of the 10-yr-long Fermi-LAT light curves to create variation patterns for the model parameters. Using as test beds two bright blazars from the SMARTS sample (PKS 2155-304 and 3C 273), we compute 10-yr-long OIR, X-ray, and γ-ray model light curves for different varying parameters. We compare the findings of our theoretical investigation with multiwavelength observations using various measures of variability. While no single-varying parameter simulation can explain all multiwavelength variability properties, changes in the electron luminosity and external radiation field in PKS 2155-304 and 3C 273, respectively, can account for most of them. Our results motivate future time-dependent studies with coupling between two or more physical parameters to describe the multiwavelength long-term blazar variability.