Variations in scaling behavior in the flux and emissions of distant astronomical sources with respect to their cosmic time are important l phenomena that can provide valuable information about the dynamics within the sources and their cosmological evolution with time. Different studies have been applying linear analysis to understand and model quasars' light curves. Here, we study the multifractal behavior of selected quasars' radio emissions in their observed frame (at 22 and 37 GHz bands) and and their rest frame. To this end, we apply the wavelet transform-based multifractal analysis formalism called wavelet transform modulus maxima. In addition, we verify whether the autoregressive integrated moving average (ARIMA) models fit our data or not. In our work, we observe strong multifractal behavior for all the sources. Additionally, we find that the degree of multifractality is strongly similar for each source and significantly different between sources at 22 and 37 GHz. This similarity implies that the two frequencies have the same radiation region and mechanism, whereas the difference indicates that the sources have intrinsically different dynamics. Furthermore, we show that the degree of multifractality is the same in the observed and rest frames of the quasars, i.e., multifractality is an intrinsic property of radio quasars. Finally, we show that the ARIMA models fit the 3C 345 quasar at 22 GHz and partially fit most of the time series with the exception of the 3C 273 and 3C 279 quasars at 37 GHz, for which the models are found to be inadequate.

We present time-resolved broad-band observations of the quasar 3C 279 obtained from multi-wavelength campaigns conducted during the first two years of the Fermi Gamma-ray Space Telescope mission. While investigating the previously reported gamma-ray/optical flare accompanied by a change in optical polarization, we found that the optical emission appears delayed with respect to the gamma-ray emission by about 10 days. X-ray observations reveal a pair of `isolated' flares separated by ~90 days, with only weak gamma-ray/optical counterparts. The spectral structure measured by Spitzer reveals a synchrotron component peaking in the mid-infrared band with a sharp break at the far-infrared band during the gamma-ray flare, while the peak appears in the mm/sub-mm band in the low state. Selected spectral energy distributions are fitted with leptonic models including Comptonization of external radiation produced in a dusty torus or the broad-line region. Adopting the interpretation of the polarization swing involving propagation of the emitting region along a curved trajectory, we can explain the evolution of the broad-band spectra during the gamma-ray flaring event by a shift of its location from ~ 1 pc to ~ 4 pc from the central black hole. On the other hand, if the gamma-ray flare is generated instead at sub-pc distance from the central black hole, the far-infrared break can be explained by synchrotron self-absorption. We also model the low spectral state, dominated by the mm/sub-mm peaking synchrotron component, and suggest that the corresponding inverse-Compton component explains the steady X-ray emission.