We present cross-correlation studies of gamma-ray (0.1-300 GeV), X-ray (0.2-10 keV) and optical (R-band) variability of a sample of 26 blazars during 2008-2016. The light curves are from Fermi-LAT, Swift-XRT, and the Yale-SMARTS blazar monitoring program. We stack the discrete cross-correlation functions of the blazars such that the features that are consistently present in a large fraction of the sample become more prominent in the final result. We repeat the same analysis for two subgroups, namely, low synchrotron peaked (LSP) and high synchrotron peaked (HSP) blazars. We find that, on average, the variability at multiple bands is correlated, with a time lag consistent with zero in both subgroups. We describe this correlation with a leptonic model of non-thermal emission from blazar jets. By comparing the model results with those from the actual data we find that the inter-band cross-correlations are consistent with an emission region of size nearly 0.1 pc within the broad line region for LSP blazars. We rule out large changes of magnetic field (> 0.5 Gauss) across the emission region or small values of magnetic field (e.g. 0.2 Gauss) for this population. We also find that the observed variability of the HSP blazars can be explained if the emission region is much larger than the distance to the broad line region from the central black hole.

The power spectral density (PSD) of the X-ray emission variability from the accretion disc-corona region of black hole X-ray binaries and active galactic nuclei has a broken power law shape with a characteristic break time-scale. If the disc and the jet are connected, the jet variability may also contain a characteristic time-scale related to that of the disc-corona. Recent observations of the blazar Mrk 421 have confirmed the broken power law shape of the PSD of its jet X-ray variability. We model the time variability of a blazar, in which emitting particles are assumed to be accelerated by successive shock waves flowing down the jet with a varying inter-shock time-scale. We investigate the possible relation between the characteristic time-scales in the disc and jet variability based on the above model, along with mathematically and physically simulated disc variability. We find that both the PSD of the jet and disc variability may have a broken power law shape but the break time-scales are not related in general except only in systems with a small range of BH mass. The break in the jet and the disc PSD are connected to the interval between large amplitude outbursts in the jet (inter-shock time-scale) and to the viscous time-scale in the disc, respectively. In frequency bands where multiple emission processes are involved or emission is from lower energy particles, the break in the PSD may not be prominent enough for detection.

We model the projected angular two-point correlation function (2PCF) of obscured and unobscured quasars selected using the Wide-field Infrared Survey Explorer (WISE), at a median redshift of $z \sim 1$ using a five-parameter Halo Occupation Distribution (HOD) parameterization, derived from a cosmological hydrodynamic simulation by Chatterjee et al. The HOD parameterization was previously used to model the 2PCF of optically selected quasars and X-ray bright active galactic nuclei (AGN) at $z \sim 1$. The current work shows that a single HOD parameterization can be used to model the population of different kinds of AGN in dark matter halos suggesting the universality of the relationship between AGN and their host dark matter halos. Our results show that the median halo mass of central quasar hosts increases from optically selected ($4.1^{+0.3}_{-0.4} \times 10^{12} \; h^{-1} \; {M_{sun}}$) and infra-red (IR) bright unobscured populations ($6.3^{+6.2}_{-2.3} \times 10^{12} \; h^{-1} \; {M_{sun}}$) to obscured quasars ($10.0^{+2.6}_{-3.7} \times 10^{12} \; h^{-1} \; {M_{sun}}$), signifying an increase in the degree of clustering. The projected satellite fractions also increase from optically bright to obscured quasars and tend to disfavor a simple `orientation only' theory of active galactic nuclei unification. Our results also show that future measurements of the small-scale clustering of obscured quasars can constrain current theories of galaxy evolution where quasars evolve from an IR- bright obscured phase to the optically bright unobscured phase.

We perform a Halo Occupation Distribution (HOD) modeling of the projected two-point correlation function (2PCF) of quasars that are observed in the Wide-field Infrared Survey Explorer (WISE) telescope with counter-parts in the Sloan Digital Sky Survey (SDSS) Data Release (DR)-8 quasar catalog at a median redshift of $z\sim 1.04 (\pm 0.58)$. Using a four parameter HOD model we derive the host mass scales of WISE selected quasars. Our results show that the median halo masses of central and satellite quasars lie in the range $M_{\mathrm{cen}} = (5 \pm 1.0) \times 10^{12} M_{\odot}$ and $M_{\mathrm{sat}} = 8 (^{+7.8} _{-4.8}) \times 10^{13} M_{\odot}$, respectively. The derived satellite fraction is $f_{\mathrm{sat}}= 5.5 (^{+35} _{-5.0})\times 10^{-3}$. Previously Richardson et al.\ used the SDSS DR7 quasar clustering data to obtain the halo mass distributions of $z\sim 1.4$ quasars. Our results on the HOD of central quasars are in excellent agreement with Richardson et al.\ but the host mass scale of satellite quasars for the WISE sample, is lower than that of Richardson et al.\ resulting in an order of magnitude higher satellite fraction for the WISE sample. We note that our sample of quasars are systematically brighter in the WISE frequency bands compared to the full quasar sample of SDSS. We discuss the implication of this result in the context of current theories of galaxy evolution.