We present the results of optical photometric observations of three extreme TeV blazars, 1ES 0229$+$200, 1ES 0414$+$009, and 1ES 2344$+$514, taken with two telescopes (1.3 m Devasthal Fast Optical Telescope, and 1.04 m Sampuranand Telescope) in India and two (1.4 m Milankovi\'{c} telescope and 60 cm Nedeljkovi\'{c} telescope) in Serbia during 2013--2019. We investigated their flux and spectral variability on diverse timescales. We examined a total of 36 intraday $R-$band light curves of these blazars for flux variations using the power-enhanced {\it F}-test and the nested ANOVA test. No significant intraday variation was detected on 35 nights, and during the one positive detection the amplitude of variability was only 2.26 per cent. On yearly timescales, all three blazars showed clear flux variations in all optical wavebands. The weighted mean optical spectral index ($\alpha_{BR}$), calculated using $B - R$ color indices, for 1ES 0229$+$200 was 2.09 $\pm$ 0.01. We also estimated the weighted mean optical spectral indices of 0.67 $\pm$ 0.01 and 1.37 $\pm$ 0.01 for 1ES 0414$+$009, and 1ES 2344$+$514, respectively, by fitting a single power-law ($F_{\nu} \propto \nu^{-\alpha}$) in their optical ({\it VRI}) spectral energy distributions. A bluer-when-brighter trend was only detected in the blazar 1ES 0414$+$009. We briefly discuss different possible physical mechanisms responsible for the observed flux and spectral changes in these blazars on diverse timescales.

We report the detection of a simultaneous quasi--periodic oscillation (QPO) in $\gamma$-ray and optical light curves of the blazar 3C 454.3. Periodic flux modulations were detected in both of these wavebands with a dominant period of $\sim 47$ days. This is the first multi-band simultaneous QPO detection in any blazar, and more generally in any AGN, on a timescale of a few tens of days. The QPO lasted for over 450 days (from MJD 56800 to 57250) resulting in over nine observed cycles which is among the highest number of periods ever detected in a blazar light curve. An auto regressive integrated moving average (ARIMA) modelling of the light curve revealed a significant, exponentially decaying, trend in the light curve during the QPO; along with the $47$ days periodicity. We explore several physical models to explain the origin of this transient quasi-periodic modulation and the overall trend in the observed flux with a month-like period. These scenarios include a binary black hole system, a hotspot orbiting close to the innermost stable circular orbit of the supermassive black hole, and precessing jets. We conclude that the most likely scenario involves a region of enhanced emission moving helically inside a curved jet. The helical motion gives rise to the QPO and the curvature ($\sim 0.05^{\circ}$ pc$^{-1}$) of the jet is responsible for the observed trend in the light curve.

We report the results of our optical (VRI) photometric observations of the TeV blazar 1ES 0806$+$524 on 153 nights during 2011-2019 using seven optical telescopes in Europe and Asia. We investigated the variability of the blazar on intraday as well as on long-term timescales. We examined eighteen intraday light curves for flux and color variations using the most reliable power-enhanced F-test and the nested ANOVA test. Only on one night was a small, but significant, variation found, in both $V$ band and $R$ band light curves. The $V-R$ color index was constant on every one of those nights. Flux density changes of around 80 % were seen over the course of these eight years in multiple bands. We found a weighted mean optical spectral index of 0.639$\pm$0.002 during our monitoring period by fitting a power law ($F_{\nu} \propto \nu^{-\alpha}$) in 23 optical ($VRI$) spectral energy distributions of 1ES 0806$+$524. We discuss different possible mechanisms responsible for blazar variability on diverse timescales.

We present the results of flux density, spectral index, and polarization intra-night monitoring studies of a sample of eight optically bright blazars, carried out by employing several small to moderate aperture (0.4\,m to 1.5\,m diameter) telescopes fitted with CCDs and polarimeters located in Europe, India, and Japan. The duty cycle of flux variability for the targets is found to be $\sim 45$ percent, similar to that reported in earlier studies. The computed two-point spectral indices are found to be between 0.65 to 1.87 for our sample, comprised of low- and intermediate frequency peaked blazars, with one exception; they are also found to be statistically variable for about half the instances where `confirmed' variability is detected in flux density. In the analysis of the spectral evolution of the targets on hourly timescale, a counter-clockwise loop (soft-lagging) is noted in the flux-spectral index plane on two occasions, and in one case a clear spectral flattening with the decreasing flux is observed. In our data set, we also observe a variety of flux-polarization degree variability patterns, including instances with a relatively straightforward anti-correlation, correlation, or counter-clockwise looping. These changes are typically reflected in the flux-polarization angle plane: the anti-correlation between the flux and polarization degree is accompanied by an anti-correlation between the polarization angle and flux, while the counter-clockwise flux-PD looping behaviour is accompanied by a clockwise looping in the flux-polarization angle representation. We discuss our findings in the framework of the internal shock scenario for blazar sources.

In order to find clues to the origin of the "winged" or "X-shaped" radio galaxies (XRGs) we investigate here the parent galaxies of a large sample of 106 XRGs for optical-radio axes alignment, interstellar medium, black hole mass, and large-scale environment. For 41 of the XRGs it was possible to determine the optical major axis and the primary radio axis and the strong tendency for the two axes to be fairly close is confirmed. However, several counter-examples were also found and these could challenge the widely discussed backflow diversion model for the origin of the radio wings. Comparison with a well-defined large sample of normal FR II radio galaxies has revealed that: (i) XRGs possess slightly less massive central black holes than the normal radio galaxies (average masses being log$M_{\rm BH} \sim$ 8.81 $M_{\odot}$ and 9.07 $M_{\odot}$, respectively); (ii) a much higher fraction of XRGs ($\sim$ 80%) exhibits red mid-IR colors ($W2 - W3 > 1.5$), indicating a population of young stars and/or an enhanced dust mass, probably due to relatively recent galaxy merger(s). A comparison of the large-scale environment (i.e., within $\sim$ 1 Mpc) shows that both XRGs and FRII radio galaxies inhabit similarly poor galaxy clustering environments (medium richness being 8.94 and 11.87, respectively). Overall, the origin of XRGs seems difficult to reconcile with a single dominant physical mechanism and competing mechanisms seem prevalent.

We performed a long-term optical (B, V, R bands), infra-red (J and K bands) and radio band (15, 22, 37 GHz band) study on the flat spectrum radio quasar, 3C 454.3, using the data collected over a period of more than 8 years (MJD 54500--57500). The temporal variability, spectral properties and inter-waveband correlations were studied by dividing the available data into smaller segments with more regular sampling. This helped us constrain the size and the relative locations of the emission regions for different wavebands. Spectral analysis of the source revealed the interplay between the accretion disk and jet emission. The source predominantly showed a redder-when-brighter trend, though we observed a bluer-when-brighter trend at high flux levels which could be signatures of particle acceleration and radiative cooling. Significant correlations with near-zero lag were seen between various optical/infra-red bands, indicating that these emission regions are co-spatial. Correlations with a time lag of about 10--100 days are seen between optical/infra-red and radio bands indicating these emissions arise from different regions. We also observe the DCF peak lag change from year to year. We try to explain these differences using a curved jet model where the different emission regions have different viewing angles resulting in a frequency dependent Doppler factor. This variable Doppler factor model explains the variability timescales and the variation in DCF peak lag between the radio and optical emissions in different segments. Lags of 6-180 days are seen between emissions in various radio bands, indicating a core-shift effect

We present X-ray flux and spectral analyses of the three pointed Suzaku observations of the TeV high synchrotron peak blazar Mrk 421 taken throughout its complete operational duration. The observation taken on 5 May 2008 is, at 364.6 kiloseconds (i.e., 101.3 hours), the longest and most evenly sampled continuous observation of this source, or any blazar, in the X-ray energy 0.8 - 60 keV until now. We found large amplitude intra-day variability in all soft and hard bands in all the light curves. The discrete correction function analysis of the light curves in soft and hard bands peaks on zero lag, showing that the emission in hard and soft bands are cospatial and emitted from the same population of leptons. The hardness ratio plots imply that the source is more variable in the harder bands compared to the softer bands. The source is harder-when-brighter, following the general behavior of high synchrotron peak blazars. Power spectral densities of all three light curves are red noise dominated, with a range of power spectra slopes. If one assumes that the emission originates very close to the central super massive black hole, a crude estimate for its mass, of ~ 4 * 10^{8} M_{\odot}, can be made; but if the variability is due to perturbations arising there that are advected into the jet and are thus Doppler boosted, substantially higher masses are consistent with the quickest seen variations. We briefly discuss the possible physical mechanisms most likely responsible for the observed flux and spectral variability.

We present a detailed spectral and temporal study of the intermediate-type blazar ON 231 during the TeV outburst phase in 2008 June with observations performed by Swift and XMM-Newton. The X-ray flux of the source, which was significantly dominated by the soft photons (below $3-4$ keV), varies between 27$\%$ and 38$\%$ on day timescales, while mild variations were observed in the optical/UV emissions. We found a maximum soft lag of $\sim 1$ hr between the UV and soft X-ray bands, which can be understood if the magnetic field of the emitting region is $\sim 5.6~ \delta^{-1/3}$ G. The $0.6-10$ keV spectra can be well represented by a broken power-law model, which indicates the presence of both synchrotron and inverse Compton components in the studied X-ray regime. The synchrotron part of the SEDs constructed with simultaneous optical/UV and X-ray data follows a log-parabolic shape. A time-resolved spectral analysis shows that the break energy varies significantly between 2.4 and 7.3 keV with the changing flux state of the source, and the similar variations of the spectral slopes of the two components support the SSC scenario. The synchrotron tail, following a log-parabolic function, shows that the peak frequency ($\nu_{p}$) varies by two orders of magnitude ($\sim 10^{14}-10^{16}$ Hz) during the event. A significantly positive $E_{p}-\beta$ relation is observed from both SED and time-resolved spectral analyses. The most feasible scenario for the observed trend during the flaring event could be associated with a magnetic-field-driven stochastic process evolving toward an equilibrium energy level.

This article focuses on certain variability and emission characteristics of Active Galactic Nuclei (AGN), especially their radio-loud subset consisting of quasars, BL Lacs and $\gamma$-ray detected narrow-line Seyfert 1 galaxies, all of which exhibit relativistically beamed jets of nonthermal radiation. Several striking trends and correlations, including some that have received scant attention, drawn from the available comparatively recent literature are highlighted. These can provide very useful inputs to models of AGN and be probed at a deeper level using the optical telescopes recently set up at Devasthal (Nainital).

We present a catalog of 290 `winged' or X-shaped radio galaxies (XRGs) extracted from the latest (17 Dec. 2014) data release of the `Very Large Array Faint Images of the Radio Sky at Twenty centimeter' (VLA FIRST survey). We have combined these radio images with their counterparts in the TIFR GMRT sky survey (TGSS) at 150 MHz (TGSS alternative data release, TGSS_ADR1), in an attempt to identify any low surface-brightness radio emission present in these sources. This has enabled us to assemble a sample of 106 `strong' XRG candidates and 184 `probable' XRG candidates whose XRG designation must be verified by further observations. The present sample of 290 XRG candidates is almost twice as large as the number of XRGs presently known. Twenty-five of our 290 XRG candidates (9 `strong' and 16 `probable') are identified as quasars. Double-peaked narrow emission lines are seen in the optical spectra of three of the XRG candidates (2 `strong' and 1 `probable'). Nearly 90% of the sample is located in the FR II domain of the Owen-Ledlow diagram. A few of the strong XRG candidates have a rather flat radio spectrum (spectral index $\alpha$ flatter than -0.3) between 150 MHz and 1.4 GHz, or between 1.4 GHz and 5 GHz. Since this is not expected for lobe-dominated extragalactic radio sources (like nearly all known XRGs), these sources are particularly suited for follow-up radio imaging and near-simultaneous measurement of the radio spectrum.

We present and analyze the variability of three Seyfert galaxies on intra-day timescales. We have analyzed in a uniform manner the 38 longest ( $>30$ ks) observations made for NGC 4051, MCG$-06-30-15$ and NGC 4151 by XMM-Newton between 2000 and 2015. The nuclei were quite active during essentially all of these observations and the overall X-ray fluxes (0.3-10 keV) varied by an order of magnitude. Most of the observations do appear to show characteristic timescales, estimated through their auto-correlation functions, ranging between $\sim 2.9$ ks and $\sim 45.3$ ks. The hard (2-10 keV) and soft (0.3-02 keV) bands are very well correlated but consideration of their hardness ratios show that the sources typically soften during flares. We also provide new estimates of the central black hole masses for these three AGN that support the hypothesis that Narrow Line Seyfert 1 galaxies have relatively small central black holes.

The latest flare of the regular $\sim$ 12 years quasi-periodic optical outbursts in the binary SMBH candidate system OJ 287 occurred in December 2015. Following this, the source has exhibited enhanced multi-wavelength (MW) variability in spectral, temporal and polarization domains with new features never seen before. Our MW investigation show that the overall MW variability can be divided into two-phase, (i) November 2015 -- May 2016 with variability from near-infrared (NIR) to Fermi-LAT $\rm \gamma$-ray energies (0.1 -- 300 GeV), and (ii) September 2016 -- July 2017 with intense NIR to X-ray variability but without any activity in the Fermi-LAT band, and the very first detection at very high energies (VHEs, E $>$ 100 GeV) by VERITAS. The broadband SEDs during the first phase show a thermal bump in the NIR-optical region and a hardening in the $\rm \gamma$-ray spectra with a shift in its peak. The thermal bump like feature is consistent with the description of the standard accretion-disk associated with the primary SMBH of mass $\sim 1.8\times10^{10} M_\odot$ while the $\rm \gamma$-ray emission can be naturally reproduced by inverse Compton scattering of photons from the broad line region which has been seen during the close encounter duration of the binary SMBHs, thereby suggesting a sub-parsec scale origin. The SEDs during the second phase (VHE detection) is a mixture of typical OJ 287 SED with hardened $\rm \gamma$-ray spectra and an HBL SED and can be explained in a two-zone model, one located at sub-parsec scales and other at parsec scales. During both the phases, the MW variability is simultaneous and almost always accompanied by changes in the polarization properties, exhibiting random and systematic variations, suggesting a strong role of magnetic field and turbulence.

We present the results of photometric observations of three TeV blazars, 3C 66A, S5 0954+658 and BL Lacertae, during the period 2013--2017. Our extensive observations were performed in a total of 360 nights which produced $\sim$6820 image frames in BVRI bands. We study flux and spectral variability of these blazars on these lengthy timescales. We also examine the optical Spectral Energy Distributions of these blazars, which are crucial in understanding the emission mechanism of long-term variability in blazars. All three TeV blazars exhibited strong flux variability during our observations. The colour variations are mildly chromatic on long timescales for two of them. The nature of the long-term variability of 3C 66A and S5 0954+658 is consistent with a model of a non-thermal variable component that has a continuous injection of relativistic electrons with power law distributions around 4.3 and 4.6, respectively. However, the long-term flux and colour variability of BL Lac suggests that these can arise from modest changes in velocities or viewing angle toward the emission region, leading to variations in the Doppler boosting of the radiation by a factor ~1.2 over the period of these observations.

We present the results of our optical (VRI) observations of the TeV blazar PG 1553+113 over eight nights in 2016 April. We monitored the blazar quasi-simultaneously in V and R bands each night and examined the light curves (LCs) for intraday flux and color variations using two of the most powerful tests; the power-enhanced F-test and the nested ANOVA test. The source was found to be significantly (> 99%) variable in both V and R bands only on April 13, while clear variations only in R band LCs were seen on April 8 and 12. No temporal variation was seen in the color during the observation period. We did not find any significant correlation between V-R color index and R magnitude on any observing night. We found a mean optical spectral index of $\sim 0.83 \pm 0.02 $ with a maximum variation of 0.21 by fitting a power law ($F_{\nu} \propto \nu^{-\alpha}$) in the optical (VRI) spectral energy distribution of PG 1553+113. We briefly discuss the possible physical processes responsible for the observed flux and spectral variability.

We report on a recent multi-band optical photometric and polarimetric observational campaign of the blazar OJ 287 which was carried out during September 2016 -- December 2017. We employed nine telescopes in Bulgaria, China, Georgia, Japan, Serbia, Spain and the United States. We collected over 1800 photometric image frames in BVRI bands and over 100 polarimetric measurements over ~175 nights. In 11 nights with many quasi-simultaneous multi-band (V, R, I) observations, we did not detect any genuine intraday variability in flux or color. On longer timescales, multiple flaring events were seen. Large changes in color with respect to time and in a color--magnitude diagram were seen, and while only a weak systematic variability trend was noticed in color with respect to time, the color--magnitude diagram shows a bluer-when-brighter trend. Large changes in the degree of polarization, and substantial swings in the polarization angle were detected. The fractional Stokes parameters of the polarization showed a systematic trend with time in the beginning of these observations, followed by chaotic changes and then an apparently systematic variation at the end. These polarization changes coincide with the detection and duration of the source at very high energies as seen by VERITAS. The spectral index shows a systematic variation with time and V-band magnitude. We briefly discuss possible physical mechanisms that could explain the observed flux, color, polarization, and spectral variability.

Since 2016 October, the active galaxy PKS 2247-131 has undergone a gamma-ray outburst, which we studied using data obtained with the Fermi Gamma-ray Space Telescope. The emission arises from a relativistic jet in PKS 2247-131, as an optical spectrum only shows a few weak absorption lines, typical of the BL Lacertae sub-class of the blazar class of active galactic nuclei. Here we report a ~34.5 day quasi-periodic oscillation (QPO) in the emission after the initial flux peak of the outburst. Compared to one-year time-scale QPOs, previously identified in blazars in Fermi energies, PKS 2247-131 exhibits the first clear case of a relatively short, month-like oscillation. We show that this QPO can be explained in terms of a helical structure in the jet, where the viewing angle to the dominant emission region in the jet undergoes periodic changes. The time scale of the QPO suggests the presence of binary supermassive black holes in PKS 2247-131.

We detected a possible quasi-periodic oscillation (QPO) of ~ 71 days in the 0.1 -- 300 GeV gamma-ray Fermi-LAT light curve of the high redshift flat spectrum radio quasar B2 1520+31. We identify and confirm that quasi-period by Lomb Scargle periodogram (LSP), and weighted wavelet z-transform (WWZ) analyses. Using this QPO period, and assuming it originates from accretion-disc fluctuations at the innermost stable circular orbit, we estimate the central supermassive black hole mass to range between ~ 5.4 * 10$^{9}$ M$_{\odot}$ for a non-rotating black hole and ~ 6.0 * 10$^{10}$ M$_{\odot}$ for a maximally rotating black hole. We briefly discuss other possible radio-loud active galactic nuclei emission models capable of producing a gamma-ray QPO of such a period in a blazar.

We present a multi-wavelength (MW) spectral and temporal study of the recent activity of the claimed super-massive binary black hole system OJ 287 since December 2015. The overall MW activity can be divided into two durations: December 2015 - April 2016 (MJD: 57360 - 57500), showing strong activity from near-infrared (NIR) to $\gamma$-rays and September 2016 - June 2017 (MJD: 57650 - 57930), showing intense NIR to X-ray variability concurrent with detection at very high energies (VHE) by VERITAS, but without any signatures of variability in the Fermi-LAT band. In the first duration, the variations are almost simultaneous and the SEDs show new components in NIR-optical and optical-UV region. The NIR-optical bump is consistent with standard accretion-disk (AD) description while the optical-UV appears consistent with contributions from the broad-line region. The extracted broadband SEDs also show a clear shift in gamma-ray SED peak and can be explained with inverse Compton scattering of photons from broad line region. In the second period, the variations are also simultaneous except for one duration during which X-ray leads the optical/UV by $\sim$ 5-6 day. The broadband SEDs, on the other hand, show mixture of a typical OJ 287 SED and an HBL SED, consistent with an origin from two different zones, one located at sub-parsec scales and other at parsec scales.

The enormous sizes and variability of emission of radio-loud AGNs arise from the relativistic flows of plasma along two oppositely directed jets. We use the Athena hydrodynamics code to simulate an extensive suite of 54 propagating three-dimensional relativistic jets with wide ranges of input jet velocities and jet-to-ambient matter density ratios. We determine which parameter sets yield unstable jets that produce jet dominated FR I type radio galaxy morphologies and which tend to produce stable jets with hot-spots and FR II morphologies. Nearly all our simulations involve jets with internal pressures matched to those of the ambient medium but we also consider over-pressured jets and discuss differences from the standard ones. We also show that the results are not strongly dependent on the adiabatic index of the fluid. We focus on simulations that remain stable for extended distances (60-240) times the initial jet radius. Scaled to the much smaller sizes probed by VLBI observations, the fluctuations in such simulated flows yield variability in observed emissivity on timescales from months. Adopting results for the densities, pressures and velocities from these simulations we estimate normalized rest frame synchrotron emissivities from individual cells in the jets. The observed emission from each cell is strongly dependent upon its variable Doppler boosting factor. We sum the fluxes from thousands of zones around the primary reconfinement shock. The light curves an power-spectra, with red-noise slopes between -2.1 and -2.5, so produced are similar to those observed from blazars.

We found a possible ~ 1 hour quasi-periodic oscillation (QPO) in a ~ 55 ks X-ray observation of the narrow-line Seyfert 1 galaxy MCG--06--30--15 made with the XMM-Newton EPIC/pn detector in the energy range 0.3 -- 10 keV. We identify a total modulation of ~ 16% in the light curve and find a $\simeq$ 3670~s quasi-period using Lomb-Scargle periodogram (LSP) and weighted wavelet Z-transform (WWZ) techniques. Our analyses of eight light curves of MCG--06--30--15, indicated the possible presence of an oscillation during one of them. The LSP indicates a statistically significant ($\simeq$ 3$\sigma$) QPO detection. A WWZ analysis shows that the signal at this possible roughly 3670s period is present, and rather persistent, throughout the observation; however, a signal around 8735s is more persistent. We briefly discuss models that can produce X-ray QPOs with such periods in narrow line Seyfert 1 galaxies, as both other claimed QPO detections in this class of AGN had very similar periods.

We present an extensive study of 72 archival Chandra light curves of the high-frequency-peaked type blazar Mrk 421, the first strong extragalactic object to be detected at TeV energies. Between 2000 and 2015 Mrk 421 often displayed intraday variability in the 0.3-10.0 keV energy range, as quantified through fractional variability amplitudes that range up to 21.3 per cent. A variability duty cycle of ~84 per cent is present in these data. Variability timescales, with values ranging from 5.5 to 30.5 ks, appear to be present in seven of these observations. Discrete correlation function analyses show positive correlations between the soft (0.3-2.0 keV) and hard (2.0-10.0 keV) X-ray energy bands with zero time lags, indicating that very similar electron populations are responsible for the emission of all the X-rays observed by Chandra. The hardness ratios of this X-ray emission indicate a general "harder-when-brighter" trend in the spectral behaviour of Mrk 421. Spectral index-flux plots provide model independent indications of the spectral evolution of the source and information on the X-ray emission mechanisms. Brief discussions of theoretical models that are consistent with these observations are given.

We present the results of our power spectral density analysis for the BL Lac object OJ\,287, utilizing the {\it Fermi}-LAT survey at high-energy $\gamma$-rays, {\it Swift}-XRT in X-rays, several ground-based telescopes and the {\it Kepler} satellite in the optical, and radio telescopes at GHz frequencies. The light curves are modeled in terms of continuous-time auto-regressive moving average (CARMA) processes. Owing to the inclusion of the {\it Kepler} data, we were able to construct \emph{for the first time} the optical variability power spectrum of a blazar without any gaps across $\sim6$ dex in temporal frequencies. Our analysis reveals that the radio power spectra are of a colored-noise type on timescales ranging from tens of years down to months, with no evidence for breaks or other spectral features. The overall optical power spectrum is also consistent with a colored noise on the variability timescales ranging from 117 years down to hours, with no hints of any quasi-periodic oscillations. The X-ray power spectrum resembles the radio and optical power spectra on the analogous timescales ranging from tens of years down to months. Finally, the $\gamma$-ray power spectrum is noticeably different from the radio, optical, and X-ray power spectra of the source: we have detected a characteristic relaxation timescale in the {\it Fermi}-LAT data, corresponding to $\sim 150$\,days, such that on timescales longer than this, the power spectrum is consistent with uncorrelated (white) noise, while on shorter variability timescales there is correlated (colored) noise.

We present a multi-wavelength spectral and temporal investigation of OJ 287 emission during its strong optical-to-X-ray activity between July 2016 - July 2017. The daily $\gamma$-ray fluxes from \emph{Fermi}-LAT are consistent with no variability. The strong optical-to-X-ray variability is accompanied by a change in power-law spectral index of the X-ray spectrum from $< 2$ to $>2$, with variations often associated with changes in optical polarization properties. Cross-correlations between optical-to-X-ray emission during four continuous segments show simultaneous optical-ultraviolet (UV) variations while the X-ray and UV/optical are simultaneous only during the middle two segments. In the first segment, the results suggest X-rays lag the optical/UV, while in the last segment X-rays lead by $\sim$ 5-6 days. The last segment also shows a systematic trend with variations appearing first at higher energies followed by lower energy ones. The LAT spectrum before the VHE activity is similar to preceding quiescent state spectrum while it hardens during VHE activity period and is consistent with the extrapolated VHE spectrum during the latter. Overall, the broadband spectral energy distributions (SEDs) during high activity periods are a combination of a typical OJ 287 SED and an HBL SED, and can be explained in a two-zone leptonic model, with the second zone located at parsec scales, beyond the broad line region, being responsible for the HBL-like spectrum. The change of polarization properties from systematic to chaotic and back to systematic, before, during and after the VHE activity, suggest dynamic roles for magnetic fields and turbulence.

We report the first results of timing and spectral studies of \textit{Nuclear Spectroscopic Telescope Array (NuSTAR)} observations of six TeV emitting high-frequency peaked blazars: 1ES 0347$-$121, 1ES 0414$+$009, RGB J0710$+$591, 1ES 1101$-$232, 1ES 1218$+$304 and H 2356-309. Two out of these six TeV blazars, 1ES 1101$-$232 and 1ES 1218$+$304, showed strong evidence of intraday variations in the 3--79 keV energy range during those observations. We also found a hint of an intraday variability timescale of 23.5 ks in the light curve of 1ES 1218$+$304 using an autocorrelation function analysis. We obtained magnetic field $B \sim 0.03$ G, electron Lorentz factor $\gamma \sim 2.16 \times 10^6 $ and emission region size $R \sim 1.19 \times 10^{16} $ cm for 1ES 1218$+$304 using that variability timescale. The other blazars' light curves do not show any variability timescales shorter than their observation lengths; however, we note that the data were both noisier and sparser for them. We also investigated the spectral shape of these TeV blazars and found that the spectrum of 1ES 0414$+$009 is well described by a single power-law with photon index $\Gamma \sim 2.77$. The spectra of the other five HBLs are somewhat better represented by log-parabola models with local photon indices (at 10 keV) $\alpha \sim 2.23-2.67$ and curvature parameters, $\beta \sim 0.27-0.43$.

We analyzed Kepler light curves of the blazar W2R 1926+42 that provided nearly continuous coverage from quarter 11 through quarter 17 (589 days between 2011 and 2013) and examined some of their flux variability properties. We investigate the possibility that the light curve is dominated by a large number of individual flares and adopt exponential rise and decay models to investigate the symmetry properties of flares. We found that those variations of W2R 1926+42 are predominantly asymmetric with weak tendencies toward positive asymmetry (rapid rise and slow decay). The durations (D) and the amplitudes (F0) of flares can be fit with log-normal distributions. The energy (E) of each flare is also estimated for the first time. There are positive correlations between logD and logE with a slope of 1.36, and between logF0 and logE with a slope of 1.12. Lomb-Scargle periodograms are used to estimate the power spectral density (PSD) shape. It is well described by a power law with an index ranging between -1.1 and -1.5. The sizes of the emission regions, R, are estimated to be in the range of 1.1*10^15 cm - 6.6*10^16 cm. The flare asymmetry is difficult to explain by a light travel time effect but may be caused by differences between the timescales for acceleration and dissipation of high-energy particles in the relativistic jet. A jet-in-jet model also could produce the observed log-normal distributions.

We monitored BL Lacertae frequently during 2014 - 2016 when it was generally in a high state. We searched for intra-day variability for 43 nights using quasi-simultaneous measurements in the B, V, R, and I bands (totaling 143 light curves); the typical sampling interval was about eight minutes. On hour-like timescales, BL Lac exhibited significant variations during 13 nights in various optical bands. Significant spectral variations are seen during most of these nights such that the optical spectrum becomes bluer when brighter. The amplitude of variability is usually greater for longer observations but is lower when BL Lac is brighter. No evidence for periodicities or characteristic variability time-scales in the light curves was found. The color variations are mildly chromatic on long timescales.

This overview provides a historical perspective highlighting the pioneering role which the fairly modest observational facilities of ARIES have played since the 1990s in systematically characterizing the optical variability on hour-like time scale (intra-night optical variability, or INOV) of several major types of high-luminosity Active Galactic Nuclei (AGN). Such information was previously available only for blazars. Similar studies have since been initiated in at least a dozen countries, giving a boost to AGN variability research. Our work has, in particular, provided strong indication that mild INOV occurs in radio-quiet QSOs (amplitude up to $\sim 3-5\%$ and duty cycle $\sim 10\%$) and, moreover, has demonstrated that similarly mild INOV is exhibited even by the vast majority of radio-loud quasars which possess powerful relativistic jets (even including many that are beamed towards us). The solitary outliers are blazars, the tiny strongly polarized subset of powerful AGN, which frequently exhibit a pronounced INOV. Among the blazars, BL Lac objects often show a bluer-when-brighter chromatic behavior, while the flat spectrum radio quasars seem not to. Quantifying any differences of INOV among the major subclasses of non-blazar type AGNs will require dedicated monitoring programs using $2-3$ metre class telescopes.

We present a multi-wavelength spectral and temporal analysis of the blazar OJ 287 during its recent activity between December 2015 -- May 2016, showing strong variability in the near-infrared (NIR) to X-ray energies with detection at $\gamma$-ray energies as well. Most of the optical flux variations exhibit strong changes in polarization angle and degree. All the inter-band time lags are consistent with simultaneous emissions. Interestingly, on days with excellent data coverage in the NIR--UV bands, the spectral energy distributions (SEDs) show signatures of bumps in the visible--UV bands, never seen before in this source. The optical bump can be explained as accretion-disk emission associated with the primary black hole of mass $\sim \rm 1.8 \times10^{10} M_{\odot}$ while the little bump feature in the optical-UV appears consistent with line emission. Further, the broadband SEDs extracted during the first flare and during a quiescent period during this span show very different $\gamma$-ray spectra compared to previously observed flare or quiescent spectra. The probable thermal bump in the visible seems to have been clearly present since May 2013, as found by examining all available NIR-optical observations, and favors the binary super-massive black hole model. The simultaneous multi-wavelength variability and relatively weak $\gamma$-ray emission that shows a shift in the SED peak is consistent with $\gamma$-ray emission originating from inverse Compton scattering of photons from the line emission that apparently contributes to the little blue bump.

The blazar 3C454.3 exhibited a strong flare seen in gamma-rays, X-rays, and optical/NIR bands during 3--12 December 2009. Emission in the V and J bands rose more gradually than did the gamma-rays and soft X-rays, though all peaked at nearly the same time. Optical polarization measurements showed dramatic changes during the flare, with a strong anti-correlation between optical flux and degree of polarization (which rose from ~ 3% to ~ 20%) during the declining phase of the flare. The flare was accompanied by large rapid swings in polarization angle of ~ 170 degree. This combination of behaviors appear to be unique. We have cm-band radio data during the same period but they show no correlation with variations at higher frequencies. Such peculiar behavior may be explained using jet models incorporating fully relativistic effects with a dominant source region moving along a helical path or by a shock-in-jet model incorporating three-dimensional radiation transfer if there is a dominant helical magnetic field. We find that spectral energy distributions at different times during the flare can be fit using modified one-zone models where only the magnetic field strength and particle break frequencies and normalizations need change. An optical spectrum taken at nearly the same time provides an estimate for the central black hole mass of ~ 2.3 * 10^9 M_sun. We also consider two weaker flares seen during the $\sim 200$ d span over which multi-band data are available. In one of them, the V and J bands appear to lead the $\gamma$-ray and X-ray bands by a few days; in the other, all variations are simultaneous.

We have examined 40 NuSTAR light curves (LCs) of five TeV emitting high synchrotron peaked blazars: 1ES 0229+200, Mrk 421, Mrk 501, 1ES 1959+650 and PKS 2155-304. Four of the blazars showed intraday variability in the NuSTAR energy range of 3-79 keV. Using an auto correlation function analysis we searched for intraday variability timescales in these LCs and found indications of several between 2.5 and 32.8 ks in eight LCs of Mrk 421, a timescale around 8.0 ks for one LC of Mrk 501, and timescales of 29.6 ks and 57.4 ks in two LCs of PKS 2155-304. The other two blazars' LCs do not show any evidence for intraday variability timescales shorter than the lengths of those observations, however, the data was both sparser and noisier, for them. We found positive correlations with zero lag between soft (3-10 keV) and hard (10-79 keV) bands for most of the LCs, indicating that their emissions originate from the same electron population. We examined spectral variability using a hardness ratio analysis and noticed a general "harder-when-brighter" behavior. The 22 LCs of Mrk 421 observed between July 2012 and April 2013 show that this source was in a quiescent state for an extended period of time and then underwent an unprecedented double peaked outburst while monitored on a daily basis during 10 - 16 April 2013. We briefly discuss models capable of explaining these blazar emissions.

We analyze the 2.5--10 keV X-ray spectra of the luminous quasar 3C 273 and simultaneous observations in UV wavelengths from XMM--Newton between 2000 and 2015. The lowest flux level ever was observed in 2015. The continuum emission from 3C 273 is generally best described by an absorbed power-law but during extremely low states the addition of fluorescence from the K-shell iron line improves the fit. We study the spectral evolution of the source during its extended quiescent state and also examine connections between the X-ray and ultraviolet emissions, which have been seen in some, but not all, previous work. We detect a possible anti-correlation between these two bands during the low state that characterized 3C 273 for most of this period; however, this was not present during a flaring state. A harder-when-brighter trend for the X-ray spectrum was observed in these long-term observations of 3C 273 for the first time. We suggest that the X-ray emission in 3C 273 is the result of inverse Compton scattering of soft UV seed photons (emitted from the local environment of the AGN), most likely in a thermal corona. We can explain the significant temporal variation of the spectral continuum as an outcome of changing optical depth of the comptonizing medium, along the lines of the wind-shock model proposed by Courvoisier and Camenzind (1989).

We present the results of our power spectral analysis for the BL Lac object PKS 0735+178 utilizing the Fermi-LAT survey at high-energy $\gamma$-rays, several ground-based optical telescopes, and single-dish radio telescopes operating at GHz frequencies. The novelty of our approach is that, by combining long-term and densely sampled intra-night light curves in the optical regime, we were able to construct for the first time the optical power spectrum of the blazar for a time domain extending from 23 years down to minutes. Our analysis reveals that: (i) the optical variability is consistent with a pure red noise, for which the power spectral density can well be approximated by a single power-law throughout the entire time domain probed; (ii) the slope of power spectral density at high-energy $\gamma$-rays ($\sim 1$), is significantly flatter than that found at radio and optical frequencies ($\sim 2$) within the corresponding time variability range; (iii) for the derived power spectra we did not detect any low-frequency flattening, nor do we see any evidence for cut-offs at the highest frequencies down to the noise floor levels due to measurement uncertainties. We interpret our findings in terms of a model where the blazar variability is generated by the underlying single stochastic process (at radio and optical frequencies), or a linear superposition of such processes (in the $\gamma$-ray regime). Along with the detailed PSD analysis, we also present the results of our extended (1998-2015) intra-night optical monitoring program and newly acquired optical photo-polarimetric data for the source.

We present recent optical photometric observations of the blazar OJ 287 taken during September 2015 -- May 2016. Our intense observations of the blazar started in November 2015 and continued until May 2016 and included detection of the large optical outburst in December 2016 that was predicted using the binary black hole model for OJ 287. For our observing campaign, we used a total of 9 ground based optical telescopes of which one is in Japan, one is in India, three are in Bulgaria, one is in Serbia, one is in Georgia, and two are in the USA. These observations were carried out in 102 nights with a total of ~ 1000 image frames in BVRI bands, though the majority were in the R band. We detected a second comparably strong flare in March 2016. In addition, we investigated multi-band flux variations, colour variations, and spectral changes in the blazar on diverse timescales as they are useful in understanding the emission mechanisms. We briefly discuss the possible physical mechanisms most likely responsible for the observed flux, colour and spectral variability.

We use the Athena hydrodynamics code to model propagating two-dimensional relativistic jets as approximations to the growth of radio-loud active galactic nuclei for various input jet velocities and jet-to-ambient matter density ratios. Using results from these simulations we estimate the changing synchrotron emission by summing the fluxes from a vertical strip of zones behind the reconfinement shock, which is nearly stationary, and from which a substantial portion of the flux variability should arise. We explore a wide range of time scales by considering two light curves from each simulation; one uses a relativistic turbulence code with bulk velocities taken from our simulations as input, while the other uses the bulk velocity data to compute fluctuations caused by variations in the Doppler boosting due to changes in the direction and the speed of the flow through all zones in the strip. We then calculate power spectral densities (PSDs) from the light curves for both turbulent and bulk velocity origins for variability. The range of the power-law slopes of the PSDs for the turbulence induced variations is -1.8 to -2.3, while for the bulk velocity produced variations this range is -2.1 to -2.9; these are in agreement with most observations. When superimposed, these power spectra span a very large range in frequency (about five decades), with the turbulent fluctuations yielding most of the shorter timescale variations and the bulk flow changes dominating the longer periods.

We present our optical photometric observations of three TeV blazars, PKS 1510-089, PG 1553+113 and Mrk 501 taken using two telescopes in India, one in Bulgaria, one in Greece and one in Serbia during 2012 - 2014. These observations covered a total of 95 nights with a total of 202 B filter frames, 247 images in V band, 817 in R band while 229 images were taken in the I filter. This work is focused on multi-band flux and colour variability studies of these blazars on diverse timescales which are useful in understanding the emission mechanisms. We studied the variability characteristics of above three blazars and found all to be active over our entire observational campaigns. We also searched for any correlation between the brightness of the sources and their colour indices. During the times of variability, no significant evidence for the sources to display spectral changes correlated with magnitude was found on timescales of a few months. We briefly discuss the possible physical mechanisms most likely responsible for the observed flux variability.

We analyzed the multi-band optical behaviour of the BL Lacertae object, S5 0716+714, during its outburst state from 2014 November - 2015 March. We took data on 23 nights at three observatories, one in India and two in Bulgaria, making quasi-simultaneous observations in B, V, R, and I bands. We measured multi-band optical fluxes, colour and spectral variations for this blazar on intraday and short timescales. The source was in a flaring state during the period analyzed and displayed intense variability in all wavelengths. R band magnitude of 11.6 was attained by the target on 18 Jan 2015, which is the brightest value ever recorded for S5 0716+714. The discrete correlation function method yielded good correlation between the bands with no measurable time lags, implying that radiation in these bands originate from the same region and by the same mechanism. We also used the structure function technique to look for characteristic timescales in the light curves. During the times of rapid variability, no evidence for the source to display spectral changes with magnitude was found on either of the timescales. The amplitude of variations tends to increase with increasing frequency with a maximum of $\sim$ 22% seen during flaring states in B band. A mild trend of larger variability amplitude as the source brightens was also found. We found the duty cycle of our source during the analyzed period to be $\sim$ 90%. We also investigated the optical spectral energy distribution of S5 0716+714 using B, V, R, and I data points for 21 nights. We briefly discuss physical mechanisms most likely responsible for its flux and spectral variations.

We have analyzed XMM-Newton observations of the high energy peaked blazar, PKS 2155-304, made on 24 May 2002 in the 0.3 - 10 keV X-ray band. These observations display a mini-flare, a nearly constant flux period and a strong flux increase. We performed a time-resolved spectral study of the data, by dividing the data into eight segments. We fitted the data with a power-law and a broken power-law model, and in some of the segments we found a noticeable spectral flattening of the source's spectrum below 10 keV. We also performed time-resolved cross-correlation analyses and detected significant hard and soft lags (for the first time in a single observation of this source) during the first and last parts of the observation, respectively. Our analysis of the spectra, the variations of photon-index with flux as well as the correlation and lags between the harder and softer X-ray bands indicate that both the particle acceleration and synchrotron cooling processes make an important contribution to the emission from this blazar. The hard lags indicate a variable acceleration process. We also estimated the magnetic field value using the soft lags. The value of the magnetic field is consistent with the values derived from the broad-band SED modeling of this source.

We observed the prototype blazar, BL Lacertae, extensively in optical and radio bands during an active phase in the period 2010--2013 when the source showed several prominent outbursts. We searched for possible correlations and time lags between the optical and radio band flux variations using multifrequency data to learn about the mechanisms producing variability. During an active phase of BL Lacertae, we searched for possible correlations and time lags between multifrequency light curves of several optical and radio bands. We tried to estimate any possible variability timescales and inter-band lags in these bands. We performed optical observations in B, V, R and I bands from seven telescopes in Bulgaria, Georgia, Greece and India and obtained radio data at 36.8, 22.2, 14.5, 8 and 4.8 GHz frequencies from three telescopes in Ukraine, Finland and USA. Significant cross-correlations between optical and radio bands are found in our observations with a delay of cm-fluxes with respect to optical ones of ~250 days. The optical and radio light curves do not show any significant timescales of variability. BL Lacertae showed many optical 'mini-flares' on short time-scales. Variations on longer term timescales are mildly chromatic with superposition of many strong optical outbursts. In radio bands, the amplitude of variability is frequency dependent. Flux variations at higher radio frequencies lead the lower frequencies by days or weeks. The optical variations are consistent with being dominated by a geometric scenario where a region of emitting plasma moves along a helical path in a relativistic jet. The frequency dependence of the variability amplitude supports an origin of the observed variations intrinsic to the source.

We present the results of extensive multi-band intra-night optical monitoring of BL Lacertae during 2010--2012. BL Lacertae was very active in this period and showed intense variability in almost all wavelengths. We extensively observed it for a total for 38 nights; on 26 of them observations were done quasi-simultaneously in B, V, R and I bands (totaling 113 light curves), with an average sampling interval of around 8 minutes. BL Lacertae showed significant variations on hour-like timescales in a total of 19 nights in different optical bands. We did not find any evidence for periodicities or characteristic variability time-scales in the light curves. The intranight variability amplitude is generally greater at higher frequencies and decreases as the source flux increases. We found spectral variations in BL Lacertae in the sense that the optical spectrum becomes flatter as the flux increases but in several flaring states deviates from the linear trend suggesting different jet components contributing to the emission at different times.

We study the core shift effect in the parsec scale jet of the blazar 3C 454.3 using the 4.8 GHz - 36.8 GHz radio light curves obtained from three decades of continuous monitoring. From a piecewise Gaussian fit to each flare, time lags $\Delta t$ between the observation frequencies $\nu$ and spectral indices $\alpha$ based on peak amplitudes $A$ are determined. From the fit $\Delta t \propto \nu^{1/k_r}$, $k_r = 1.10 \pm 0.18$ indicating equipartition between the magnetic field energy density and the particle energy density. From the fit $A \propto \nu^\alpha$, $\alpha$ is in the range $-0.24$ to $1.52$. A mean magnetic field strength at 1 pc, $B_1 = 0.5 \pm 0.2$ G, and at the core, $B_{\rm core} = 46 \pm 16$ mG, are inferred, consistent with previous estimates. The measure of core position offset is $\Omega_{r\nu} = 6.4 \pm 2.8$ pc GHz$^{1/k_r}$ when averaged over all frequency pairs. Based on the statistical trend shown by the measured core radius $r_{\rm core}$ as a function of $\nu$, we infer that the synchrotron opacity model may not be valid for all cases. A Fourier periodogram analysis yields power law slopes in the range $-1.6$ to $-3.5$ describing the power spectral density shape and gives bend timescales in the range $0.52 - 0.66~$yr. This result, and both positive and negative $\alpha$, indicate that the flares originate from multiple shocks in a small region. Important objectives met in our study include: the demonstration of the computational efficiency and statistical basis of the piecewise Gaussian fit; consistency with previously reported results; evidence for the core shift dependence on observation frequency and its utility in jet diagnostics in the region close to the resolving limit of very long baseline interferometry observations.

To search for optical variability on a wide range of timescales, we have carried out photometric monitoring of two flat spectrum radio quasars, 3C 454.3 and 3C 279, plus one BL Lac, S5 0716+714, all of which have been exhibiting remarkably high activity and pronounced variability at all wavelengths. CCD magnitudes in B, V, R and I pass-bands were determined for $\sim$ 7000 new optical observations from 114 nights made during 2011 - 2014, with an average length of $\sim$ 4 h each, at seven optical telescopes: four in Bulgaria, one in Greece, and two in India. We measured multiband optical flux and colour variations on diverse timescales. Discrete correlation functions were computed among B, V, R, and I observations, to search for any time delays. We found weak correlations in some cases with no significant time lags. The structure function method was used to estimate any characteristic time-scales of variability. We also investigated the spectral energy distribution of the three blazars using B, V, R, I, J and K pass-band data. We found that the sources almost always follows a bluer-when-brighter trend. We discuss possible physical causes of the observed spectral variability.

We have undertaken a nearly simultaneous optical/UV and X-ray variability study of the flat spectrum radio quasar, 3C 273 using data available from the XMM$-$Newton satellite mission from June 2000 to July 2012. Here we focus on the multi-wavelength flux variability on both intra-day and long time scales of this very well known radio-loud source. We found high flux variability over long time scales in all bands for which observations were made. The optical/UV variability amplitude was more than twice than that in the X-ray bands. There is some frequency dependence of the variability in optical/UV bands in the sense that the variability amplitude increases with increasing frequency; however, the X-ray emissions disagree with this trend as the variability amplitude decreases from soft to hard X-ray bands. On intra-day time scales 3C 273 showed small amplitude variability in X-ray bands. A hardness ratio analysis in the X-ray regime indicates that the particle acceleration mechanism dominates the cooling mechanism during most of the ~12 year span of these observations.

We present a numerical model developed to calculate observed fluxes of relativistic jets in active galactic nuclei. The observed flux of each turbulent eddy is dependent upon its variable Doppler boosting factor, computed as a function of the relativistic sum of the individual eddy and bulk jet velocities and our viewing angle to the jet. The total observed flux is found by integrating the radiation from the eddies over the turbulent spectrum. We consider jets that contain turbulent eddies that have either standard Kolmogorov or recently derived relativistic turbulence spectra. We also account for the time delays in receiving the emission of the eddies due to their different simulated positions in the jet, as well as due to the varying beaming directions as they turn over. We examine these theoretical light curves and compute power spectral densities (PSDs) for a range of viewing angles, bulk velocities of the jet, and turbulent velocities. These PSD slopes depend significantly on the turbulent velocity and are essentially independent of viewing angle and bulk velocity. The flux variations produced in the simulations for realistic values of the parameters tested are consistent with the types of variations observed in radio-loud AGN as, for example, recently measured with the Kepler satellite, as long as the turbulent velocities are not too high.

We introduce a novel formalism to investigate the role of the spin angular momentum of astrophysical black holes in influencing the behaviour of low angular momentum general relativistic accretion. We propose a metric independent analysis of axisymmetric general relativistic flow, and consequently formulate the space and time dependent equations describing the general relativistic hydrodynamic accretion flow in the Kerr metric. The associated stationary critical solutions for such flow equations are provided and the stability of the stationary transonic configuration is examined using an elegant linear perturbation technique. We examine the properties of infalling material for both prograde and retrograde accretion as a function of the Kerr parameter at extremely close proximity to the event horizon. Our formalism can be used to identify a new spectral signature of black hole spin, and has the potential of performing the black hole shadow imaging corresponding to the low angular momentum accretion flow.

We have used all 20 archival XMM-Newton observations of PKS 2155-304 with simultaneous X-ray and UV/optical data to study its long term flux and spectral variability. We find significant variations, in all bands, on time-scales of years with an rms amplitude of ~ 35-45 per cent, through the optical/UV variations are not correlated with those in the X-ray. We constructed SEDs that span more than three orders of magnitude in frequency and we first fitted them with a log-parabolic model; such models have been applied many times in the past for this, and other, blazars. These fits were poor, so we then examined combined power-law and log-parabolic fits that are improvements. These models indicate that the optical/UV and X-ray flux variations are mainly driven by model normalization variations, but the X-ray band flux is also affected by spectral variations, as parametrized with the model 'curvature' parameter, b. Overall, the energy at which the emitted power is maximum correlates positively with the total flux. As the spectrum shifts to higher frequencies, the spectral 'curvature' increase, in contrast to what is expected if a single log-parabolic model were an acceptable representation of the broad band SEDs. Our results suggest that the optical/UV and X-ray emissions in this source may arise from different lepton populations.

We present the first results of a systematic search for the rare extragalactic radio sources showing an inverted (integrated) spectrum, with spectral index $\alpha \ge +2.0$, a previously unexplored spectral domain. The search is expected to yield strong candidates for $\alpha \ge +2.5$, for which the standard synchrotron self-absorption (characterized by a single power-law energy distribution of relativistic electron population) would not be a plausible explanation, even in an ideal case of a perfectly homogeneous source of incoherent synchrotron radiation. Such sharply inverted spectra, if found, would require alternative explanations, e.g., free-free absorption, or non-standard energy distribution of relativistic electrons which differs from a power-law (e.g., Maxwellian). The search was carried out by comparing two sensitive low-frequency radio surveys made with sub-arcminute resolution, namely, the WISH survey at 352 MHz and TGSS/DR5 at 150 MHz. The overlap region between these two surveys contains 7056 WISH sources classified as `single' and brighter than 100 mJy at 352 MHz. We focus here on the seven of these sources for which we find $\alpha > +2.0$. Two of these are undetected at 150 MHz and are particularly good candidates for $\alpha > +2.5$. Five of the seven sources exhibit a `Gigahertz-Peaked-Spectrum' (GPS).

We have used photometry from the Kepler satellite to characterize the variability of four radio-loud active galactic nuclei (AGN) on timescales from years to minutes. The Kepler satellite produced nearly continuous high precision data sets which provided better temporal coverage than possible with ground based observations. We have now accumulated eleven quarters of data, eight of which were reported in our previous paper. In addition to constructing power spectral densities (PSDs) and characterizing the variability of the last three quarters, we have linked together the individual quarters using a multiplicative scaling process, providing data sets spanning ~2.8 years with >98% coverage at a 30 minute sampling rate. We compute PSDs on these connected data sets that yield power law slopes at low frequencies in the approximate range of -1.5 to -2.0, with white noise seen at higher frequencies. These PSDs are similar to those of both the individual quarters and to those of ground based optical observations of other AGN. We also have explored a PSD binning method intended to reduce a bias toward shallow slope fits by evenly distributing the points within the PSDs. This tends to steepen the computed PSD slopes, especially when the low frequencies are relatively poorly fit. We detected flares lasting several days in which the brightness increased by ~15-20% in one object, as well a smaller flare in another. Two AGN showed only small, ~1-2%, fluctuations in brightness.

We discuss the origin of the low-energy cutoff, or LEC, seen in the radio spectra of many extragalactic jets and relate this to the spin of the supermassive black holes that presumably power them. Pion decay via proton-proton collisions is a possible mechanism to supply a secondary positron population with a low energy limit. We expect that pion production would occur in advection dominated accretion flows or ADAFs. In radiatively inefficient ADAFs the heat energy of the accreting gas is unable to radiate in less than the accretion time and the particle temperature could be high enough so that thermal protons can yield such pion production. Strong starbursts are another option for the injection of a truncated particle population into the jet. The role of both mechanisms is discussed with respect to the black hole spin estimate. The energy demanded to produce the pion decay process involves a minimum threshold for kinetic energy of the interacting protons. Therefore the mean proton speed in the flow can determine whether a LEC is generated. In ADAFs the random velocity of the protons can exceed the minimum speed limit of pion production around the jet launching region in the innermost part of the flow. Finally we summarize the additional work needed to put the model assumptions on a more rigorous basis.

We present the results of photometric (V band) and polarimetric observations of the blazar BL Lac during 2008--2010 using TRISPEC attached to the KANATA 1.5-m telescope in Japan. The data reveal a great deal of variability ranging from days to months with detection of strong variations in fractional polarization. The V band flux strongly anti-correlates with the degree of polarization during the first of two observing seasons but not during the second. The direction of the electric vector, however, remained roughly constant during all our observations. These results are consistent with a model with at least two emission regions being present, with the more variable component having a polarization direction nearly perpendicular to that of the relatively quiescent region so that a rising flux can produce a decline in degree of polarization. We also computed models involving helical jet structures and single transverse shocks in jets and show that they might also be able to agree with the anti-correlations between flux and fractional polarization.

We present a new type of observation relating the gravitational wave emission of supermassive black hole mergers to their electromagnetic counterparts. This dual emission involves variability of a relativistic jet arising from the spin-orbit precession of the supermassive black hole binary at its base.

The incidence of intra-night optical variability (INOV) is known to to differ significantly among different classes of powerful active galactic nuclei (AGN). A number of statistical methods have been employed in the literature for testing the presence of INOV in the light curves, sometimes leading to discordant results. In this paper we compare the INOV characteristics of six prominent classes of AGN, as evaluated using three commonly used statistical tests, namely the $\chi^2-$test, the modified $C-$test and the $F-$test, which has recently begun to gain popularity. The AGN classes considered are: radio-quiet quasars (RQQs), radio-intermediate quasars (RIQs), lobe-dominated quasars (LDQs), low optical polarization core-dominated quasars (LPCDQs), high optical polarization core-dominated quasars (HPCDQs), and TeV blazars. Our analysis is based on a large body of AGN monitoring data, involving 262 sessions of intra-night monitoring of a total 77 AGN, using 1-2 metre class optical telescopes located in India. In order to compare the usefulness of the statistical tests, we have also subjected them to a `sanity check' by comparing the number of false positives yielded by each test with the corresponding statistical prediction. The present analysis is intended to serve as a benchmark for future INOV studies of AGN of different classes.

It is important to quantify the underestimation of rms photometric errors returned by the commonly used $\emph APPHOT$ algorithm in the $\emph IRAF$ software, in the context of differential photometry of point-like AGN, because of the crucial role it plays in evaluating their variability properties. Published values of the underestimation factor, $\eta$, using several different telescopes, lie in the range 1.3 - 1.75. The present study aims to revisit this question by employing an exceptionally large data set of 262 differential light curves (DLCs) derived from 262 pairs of non-varying stars monitored under our ARIES AGN monitoring program for characterizing the intra-night optical variability (INOV) of prominent AGN classes. The bulk of these data were taken with the 1-m Sampurnanad Telescope (ST). We find $\eta$ = 1.54$\pm$0.05 which is close to our recently reported value of $\eta$ = 1.5. Moreover, this consistency holds at least up to a brightness mismatch of 1.5 mag between the paired stars. From this we infer that a magnitude difference of at least up to 1.5 mag between a point-like AGN and comparison star(s) monitored simultaneously is within the same CCD chip acceptable, as it should not lead to spurious claims of INOV.

We have used Kepler photometry to characterize variability in four radio-loud active galactic nuclei (three quasars and one object tentatively identified as a Seyfert 1.5 galaxy) on timescales from minutes to months, comparable to the light crossing time of the accretion disk around the central supermassive black hole or the base of the relativistic jet. Kepler's almost continuous observations provide much better temporal coverage than is possible from ground-based observations. We report the first such data analyzed for quasars. We have constructed power spectral densities using 8 Kepler quarters of long-cadence (30-minute) data for three AGN, 6 quarters for one AGN and 2 quarters of short-cadence (1-minute) data for all four AGN. On timescales longer than about 0.2-0.6 day, we find red noise with mean power-law slopes ranging from -1.8 to -1.2, consistent with the variability originating in turbulence either behind a shock or within an accretion disk. Each AGN has a range of red noise slopes which vary slightly by month and quarter of observation. No quasi-periodic oscillations of astrophysical origin were detected. We detected several days-long flares when brightness increased by 3% - 7% in two objects. No flares on timescales of minutes to hours were detected. Our observations imply that the duty cycle for enhanced activity in these radio-loud AGN is small. These well-sampled AGN light curves provide an impetus to develop more detailed models of turbulence in jets and instabilities in accretion disks.

We report the results of optical monitoring for a sample of 11 blazars including 10 BL Lacs and 1 Flat Spectrum Radio Quasar (FSRQ). We have measured the multiband optical flux and colour variations in these blazars on intra-day and short-term timescales of months and have limited data for 2 more blazars. These photometric observations were made during 2009 to 2011, using six optical telescopes, four in Bulgaria, one in Greece and one in India. On short-term timescales we found significant flux variations in 9 of the sources and colour variations in 3 of them. Intra-day variability was detected on 6 nights for 2 sources out of the 18 nights and 4 sources for which we collected such data. These new optical observations of these blazars plus data from our previous published papers (for 3 more blazars) were used to analyze their spectral flux distributions in the optical frequency range. Our full sample for this purpose includes 6 high-synchrotron-frequency-peaked BL Lacs (HSPs), 3 intermediate-synchrotron-frequency-peaked BL Lacs (ISPs) and 6 low-synchrotron-frequency-peaked BL Lacs (LSPs; including both BL Lacs and FSRQs). We also investigated the spectral slope variability and found that the average spectral slopes of LSPs show a good accordance with the Synchrotron Self-Compton (SSC) loss dominated model. Our analysis supports previous studies that found that the spectra of the HSPs and FSRQs have significant additional emission components. The spectra of all these HSPs and LSPs get flatter when they become brighter, while for FSRQs the opposite appears to hold. This supports the hypothesis that there is a significant thermal contribution to the optical spectrum for FSRQs.

We report results from a 1 week multi-wavelength campaign to monitor the BL Lac object S5 0716+714 (on December 9-16, 2009). In the radio bands the source shows rapid (~ (0.5-1.5) day) intra-day variability with peak amplitudes of up to ~ 10 %. The variability at 2.8 cm leads by about 1 day the variability at 6 cm and 11 cm. This time lag and more rapid variations suggests an intrinsic contribution to the source's intraday variability at 2.8 cm, while at 6 cm and 11 cm interstellar scintillation (ISS) seems to predominate. Large and quasi-sinusoidal variations of ~ 0.8 mag were detected in the V, R and I-bands. The X-ray data (0.2-10 keV) do not reveal significant variability on a 4 day time scale, favoring reprocessed inverse-Compton over synchrotron radiation in this band. The characteristic variability time scales in radio and optical bands are similar. A quasi-periodic variation (QPO) of 0.9 - 1.1 days in the optical data may be present, but if so it is marginal and limited to 2.2 cycles. Cross-correlations between radio and optical are discussed. The lack of a strong radio-optical correlation indicates different physical causes of variability (ISS at long radio wavelengths, source intrinsic origin in the optical), and is consistent with a high jet opacity and a compact synchrotron component peaking at ~= 100 GHz in an ongoing very prominent flux density outburst. For the campaign period, we construct a quasi-simultaneous spectral energy distribution (SED), including gamma-ray data from the FERMI satellite. We obtain lower limits for the relativistic Doppler-boosting of delta >= 12-26, which for a BL\,Lac type object, is remarkably high.

{Abridged} Rapid variations in optical flux are seen in many quasars and all blazars. The amount of variability in different classes of Active Galactic Nuclei has been studied extensively but many questions remain unanswered. We present the results of a long-term programme to investigate the intra-night optical variability (INOV) of powerful flat spectrum radio core-dominated quasars (CDQs), with a focus on probing the relationship of INOV to the degree of optical polarization. We observed a sample of 16 bright CDQs showing strong broad optical emission lines and consisting of both high and low optical polarization quasars (HPCDQs and LPCDQs). We employed ARIES, IIA, IGO telescopes, to carry out {\it R}-band monitoring on a total of 47 nights. Combining these INOV data with those taken from the literature, we were able to increase the sample size to 21 CDQs(12 LPCDQs and 9 HPCDQs) monitored on a total of 73 nights. As the existence of a prominent flat-spectrum radio core signifies that strong relativistic beaming is present in all these CDQs, the definitions of the two sets differ primarily in fractional optical polarization, the LPCDQs showing a very low median$ P_{op} \simeq$ 0.4 per cent. Our study yields an INOV duty cycle (DC) of $\sim$28 per cent for the LPCDQs and $\sim 68$ percent for HPCDQs. If only strong INOV with fractional amplitude above 3 per cent is considered, the corresponding DCs are $\sim$ 7 per cent and $\sim$ 40 per cent, respectively.From this strong contrast between the two classes of luminous, relativistically beamed quasars, it is apparent that relativistic beaming is normally not a sufficient condition for strong INOV and a high optical polarization is the other necessary condition.

We discuss the nature of the multi-component radio continuum and HI emission associated with the nearby galaxy group comprised of two dominant ellipticals, NGC 5898 and NGC 5903, and a dwarf lenticular ESO514-G003. Striking new details of radio emission are unveiled from the 2nd Data Release of the ongoing TIFR.GMRT.SKY.SURVEY (TGSS) which provides images with a resolution of ~ 24 arcsec x 18 arcsec and a typical rms noise of 5 mJy at 150 MHz. Previous radio observations of this compact triplet of galaxies include images at higher frequencies of the radio continuum as well as HI emission, the latter showing huge HI trails originating from the vicinity of NGC 5903 where HI is in a kinematically disturbed state. The TGSS 150 MHz image has revealed a large asymmetric radio halo around NGC 5903 and also established that the dwarf SO galaxy ESO514-G003 is the host to a previously known bright double radio source. The radio emission from NGC 5903 is found to have a very steep radio spectrum ({\alpha} ~ -1.5) and to envelope a network of radio continuum filaments bearing a spatial relationship to the HI trails. Another noteworthy aspect of this triplet of early-type galaxies highlighted by the present study is that both its radio loud members, namely NGC 5903 and ESO514-G003, are also the only galaxies that are seen to be connected to an HI filament. This correlation is consistent with the premise that cold gas accretion is of prime importance for triggering powerful jet activity in the nuclei of early-type galaxies.

We report the results of photometric observations of the blazars Mrk 421 and 3C 454.3 designed to search for intraday variability (IDV) and short-term variability (STV). Optical photometric observations were spread over eighteen nights for Mrk 421 and seven nights for 3C 454.3 during our observing run in 2009-2010 at the 1.04 m telescope at ARIES, India. Genuine IDV is found for the source 3C 454.3 but not for Mrk 421. Genuine STV is found for both sources. Mrk 421 was revealed by the MAXI X-ray detector on the International Space Station to be in an exceptionally high flux state in 2010 January - February. We performed a correlation between the X-ray and optical bands to search for time delays and found a weak correlation with higher frequencies leading the lower frequencies by about ten days. The blazar 3C 454.3 was found to be in high flux state in November-December 2009. We performed correlations in optical observations made at three telescopes, along with X-ray data from the MAXI satellite and public release gamma-ray data from the Fermi space telescope. We found strong correlations between the gamma-ray and optical bands at a time lag of about four days but the X-ray flux is not correlated with either. We briefly discuss the possible reasons for the time delays between these bands within the framework of existing models for X-ray and gamma-ray emission mechanisms.

We report the results of quasi-simultaneous two filter optical monitoring of two high-energy peaked blazars, 1ES 1959+650 and 1ES 2344+514, to search for microvariability and short-term variability (STV). We carried out optical photometric monitoring of these sources in an alternating sequence of B and R pass-bands, and have 24 and 19 nights of new data for these two sources, respectively. No genuine microvariability (intra-night variability) was detected in either of these sources. This non-detection of intra-night variations is in agreement with the conclusions of previous studies that high-energy peaked BL Lacs are intrinsically less variable than low-energy peaked BL Lacs in the optical bands. We also report the results of STV studies for these two sources between July 2009 and August 2010. Genuine STV is found for the source 1ES 1959+650 but not for 1ES 2344+514. We briefly discuss possible reasons for the difference between the intra-night variability behaviour of high- and low-energy peaked blazars.

After a brief, critical review of the leading explanations proposed for the small but important subset of radio galaxies showing an X-shaped morphology (XRGs) we propose a generalized model, based on the jet-shell interaction and spin-flip hypotheses. The most popular scenarios for this intriguing phenomenon invoke either hydrodynamical backflows and over-pressured cocoons or rapid jet reorientations, presumably from the spin-flips of central engines following the mergers of pairs of galaxies, each of which contains a supermassive black hole (SMBH). We confront these models with a number of key observations and thus argue that none of the models is capable of explaining the entire range of salient observational properties of XRGs, although some of the arguments raised in the literature against the spin-flip scenario are probably not tenable. We then propose here a new scenario which also involves galactic mergers but would allow the spin of the central engine to maintain its direction. Motivated by the detailed multi-band observations of the nearest radio galaxy, Centaurus A, this new model emphasizes the role of interactions between the jets and the shells of stars and gas that form and rotate around the merged galaxy and can cause temporary deflections of the jets, occasionally giving rise to an X-shaped radio structure. Although each of the models is likely to be relevant to a subset of XRGs, the bulk of the evidence indicates that most of them are best explained by the jet-shell interaction or spin-flip hypotheses.

In the context of AGN unification scheme rapid variability properties play an important role in understanding any intrinsic differences between sources in different classes. In this respect any clue based on spectral properties will be very useful toward understanding the mechanisms responsible for the origin of rapid small scale optical variations, or microvariability. Here we have added spectra of 46 radio-quiet quasars (RQQSOs) and Seyfert 1 galaxies to those of our previous sample of 37 such objects, all of which had been previously searched for microvariability. We took new optical spectra of 33 objects and obtained 13 others from the literature. Their \hbeta and \mgii emission lines were carefully fit to determine line widths (FWHM) as well as equivalent widths (EW) due to the broad emission line components. The line widths were used to estimate black hole masses and Eddington ratios, $\ell$. Both EW and FWHM are anticorrelated with $\ell$. Nearly all trends were in agreement with our previous work, although the tendency for sources exhibiting microvariability to be of lower luminosity was not confirmed. Most importantly, this whole sample of EW distributions provides no evidence for the hypothesis that a weak jet component in the radio quiet AGNs is responsible for their microvariability.

We present the results of our study of spectral energy distributions (SEDs) of a sample of ten low- to intermediate-synchrotron-peaked blazars. We investigate some of the physical parameters most likely responsible for the observed short-term variations in blazars. To do so, we focus on the study of changes in the SEDs of blazars corresponding to changes in their respective optical fluxes. We model the observed spectra of blazars from radio to optical frequencies using a synchrotron model that entails a log-parabolic distribution of electron energies. A significant correlation among the two fitted spectral parameters ($a$, $b$) of log-parabolic curves and a negative trend among the peak frequency and spectral curvature parameter, $b$, emphasize that the SEDs of blazars are fitted well by log-parabolic curves. On considering each model parameter that could be responsible for changes in the observed SEDs of these blazars, we find that changes in the jet Doppler factors are most important.

In this first systematic attempt to characterise the intranight optical variability (INOV) of TeV detected blazars, we have monitored a well defined set of 9 TeV blazars on total 26 nights during 2004-2010. In this R (or V)-band monitoring programme only one blazar was monitored per night for a minimum duration of 4 hours. Using the CCD, an INOV detection threshold of ~ 1-2 % was achieved in the densely sampled DLCs. We have further expanded the sample by including another 13 TeV blazars from literature. This enlarged sample of 22 TeV blazars, monitored on a total of 116 nights (including 55 nights newly reported here), has enabled us to arrive at the first estimate of the INOV duty cycle of TeV detected blazars. Applying the C-test, the INOV DC is found to be 59 %, which decreases to 47 % if only INOV fractional amplitudes above 3 % are considered. These observations also permit, for the first time, a comparison of the INOV characteristics of the two major subclasses of TeV detected BL Lacs, namely LBLs and HBLs, for which we find the INOV DCs to be ~ 63 % and ~ 38 %, respectively. This demonstrates that the INOV differential between LBLs and HBLs persists even when only their TeV detected subsets are considered. Despite dense sampling, the intranight light curves of the 22 TeV blazars have not revealed even a single feature on time scale substantially shorter than 1 hour, even though the inner jets of TeV blazars are believed to have exceptionally large bulk Lorentz factors (and correspondingly stronger time compression). An intriguing feature, clearly detected in the light curve of the HBL J1555+1111, is a 4 per cent `dip' on a 1 hour timescale. This unique feature could have arisen from absorption in a dusty gas cloud, occulting a superluminally moving optical knot in the parsec scale jet of this relatively luminous BL Lacs object.

(Abridged) With the largest homogeneous quasar sample with optical spectra and robust radio morphology classifications assembled to date, we investigate quasar radio and optical properties with unprecedented statistical power. The data consist of 4714 radio quasars from FIRST with S_{20}>2mJy and SDSS spectra. Radio morphology classes include core-only (core), core-lobe (lobe), core-jet (jet), lobe-core-lobe (triple), and double-lobe. We examine the optical colors of radio-morphology subsamples and find that radio quasars with core emission unresolved by FIRST (on 5" scale) have a redder color distribution than radio-quiet quasars (S_{20}<1mJy); other classes of radio quasars have optical color distributions similar to the radio-quiet quasars. This analysis also suggests that optical colors of z<2.7 SDSS quasars are not strongly (<0.1mag) biased blue. We show that the radio core-to-lobe flux density ratio (R) and the radio-to-optical (i-band) ratio of the quasar core (RI) are correlated, suggesting that both parameters are indicative of line-of-sight orientation. We investigate spectral line equivalent widths as a function of R and RI, including the [OIII] narrow line doublet and the CIV {\lambda}1549 and MgII {\lambda}2799 broad lines. We find that the rest equivalent widths (EWs) of the broad lines correlate positively with RI at the 4-8{\sigma} level. But we find no strong dependence of EW on R, in contrast to previous results. One interpretation is that EWs increase as the line-of-sight angle to the radio-jet axis decreases. These results are in stark contrast to commonly accepted orientation-based theories, which suggest that continuum emission should increase as the angle to the radio-jet axis decreases, resulting in smaller EW of emission lines (assumed isotropic). Finally, we find that the Baldwin effect in our sample does not depend on quasar radio morphology.

We have carried out optical (R band) intraday variability (IDV) monitoring of a sample of ten bright low energy peaked blazars (LBLs). Forty photometric observations, of an average of ~ 4 hours each, were made between 2008 September and 2009 June using two telescopes in India. Measurements with good signal to noise ratios were typically obtained within 1-3 minutes, allowing the detection of weak, fast variations using N-star differential photometry. We employed both structure function and discrete correlation function analysis methods to estimate any dominant timescales of variability and found that in most of the cases any such timescales were longer than the duration of the observation. The calculated duty cycle of IDV in LBLs during our observing run is ~ 52%, which is low compared to many earlier studies; however, the relatively short periods for which each source was observed can probably explain this difference. We briefly discuss possible emission mechanisms for the observed variability.

We present a study of the extended radio emission in a sample of 8434 low redshift (z < 0.35) broad line active galactic nuclei (AGN) from the Sloan Digital Sky Survey (SDSS). To calculate the jet and lobe contributions to the total radio luminosity, we have taken the 846 radio core sources detected in our previous study of this sample and performed a systematic search in the Faint Images of the Radio Sky at Twenty-centimeters (FIRST) database for extended radio emission that is likely associated with the optical counterparts. We found 51 out of 846 radio core sources have extended emission (> 4" from the optical AGN) that is positively associated with the AGN, and we have identified an additional 12 AGN with extended radio emission but no detectable radio core emission. Among these 63 AGN, we found 6 giant radio galaxies (GRGs), with projected emission exceeding 750 kpc in length, and several other AGN with unusual radio morphologies also seen in higher redshift surveys. The optical spectra of many of the extended sources are similar to that of typical broad line radio galaxy spectra, having broad H$\alpha$ emission lines with boxy profiles and large M_BH. With extended emission taken into account, we find strong evidence for a bimodal distribution in the radio-loudness parameter R, where the lower radio luminosity core-only sources appear as a population separate from the extended sources, with a dividing line at log(R) $\approx 1.75$. This dividing line ensures that these are indeed the most radio-loud AGN, which may have different or extreme physical conditions in their central engines when compared to the more numerous radio quiet AGN.

The emission from blazars is known to be variable at all wavelengths. The flux variability is often accompanied by spectral changes. Spectral energy distribution (SED) changes must be associated with changes in the spectra of emitting electrons and/or the physical parameters of the jet. Meaningful modeling of blazar broadband spectra is required to understand the extreme conditions within the emission region. Not only is the broadband SED crucial, but also information about its variability is needed to understand how the highest states of emission occur and how they differ from the low states. This may help in discriminating between models. Here we present the results of our SED modeling of the blazar S5 0716+714 during various phases of its activity. The SEDs are classified into different bins depending on the optical brightness state of the source.

We present optical light curves of 19 radio quiet (RQ) broad absorption line (BAL) QSOs and study their rapid variability characteristics. Systematic CCD observations, aided by a careful data analysis procedure, have allowed us to clearly detect any such microvariability exceeding 0.01--0.02 mag. Our observations cover a total of 13 nights (~72 hours) with each quasar monitored for about 4 hours on a given night. Our sample size is a factor of three larger than the number of radio-quiet BALQSOs previously searched for microvariability. We introduce a scaled F-test statistic for evaluating the presence of optical microvariability and demonstrate why it is generally preferable to the statistics usually employed for this purpose. Considering only unambiguous detections of microvariability we find that ~11 per cent of radio-quiet BALQSOs (two out of 19 sources) show microvariability for an individual observation length of about 4 hr. This new duty cycle of 11 per cent is similar to the usual low microvariability fraction of normal RQQSOs with observation lengths similar to those of ours. This result provides support for models where radio-quiet BALQSO do not appear to be a special case of the RQQSOs in terms of their microvariability properties.

The ideal gas equation of state with a constant adiabatic index, although commonly used in relativistic hydrodynamics, is a poor approximation for most relativistic astrophysical flows. Here we propose a new general equation of state for a multi-component relativistic gas which is consistent with the Synge equation of state for a relativistic perfect gas and is suitable for numerical (special) relativistic hydrodynamics. We also present a multidimensional relativistic hydrodynamics code incorporating the proposed general equation of state, based on the HLL scheme, which does not make use of a full characteristic decomposition of the relativistic hydrodynamic equations. The accuracy and robustness of this code is demonstrated in multidimensional calculations through several highly relativistic test problems taking into account nonvanishing tangential velocities. Results from three-dimensional simulations of relativistic jets show that the morphology and dynamics of the relativistic jets are significantly influenced by the different equation of state and by different compositions of relativistic perfect gases. Our new numerical code, combined with our proposed equation of state is very efficient and robust, and unlike previous codes, it gives very accurate results for thermodynamic variables in relativistic astrophysical flows.

Over the course of three hours on 27 December 2008 we obtained optical (R-band) observations of the blazar S5 0716+714 at a very fast cadence of 10 s. Using several different techniques we find fluctuations with an approximately 15-minute quasi-period to be present in the first portion of that data at a > 3 sigma confidence level. This is the fastest QPO that has been claimed to be observed in any blazar at any wavelength. While this data is insufficient to strongly constrain models for such fluctuations, the presence of such a short timescale when the source is not in a very low state seems to favor the action of turbulence behind a shock in the blazar's relativistic jet.

Ultra high energy cosmic rays (UHECRs), with energies above ~6 x 10^19 eV, seem to show a weak correlation with the distribution of matter relatively near to us in the universe. It has earlier been proposed that UHECRs could be accelerated in either the nucleus or the outer lobes of the nearby radio galaxy Cen A. We show that UHECR production at a spatially intermediate location about 15 kpc northeast from the nucleus, where the jet emerging from the nucleus is observed to strike a large star-forming shell of gas, is a plausible alternative. A relativistic jet is capable of accelerating lower-energy heavy seed cosmic rays (CRs) to UHECRs on timescales comparable to the time it takes the jet to pierce the large gaseous cloud. In this model many CRs arising from a starburst, with a composition enhanced in heavy elements near the knee region around PeV, are boosted to ultra-high energies by the relativistic shock of a newly oriented jet. This model matches the overall spectrum shown by the Auger data and also makes a prediction for the chemical composition as a function of particle energy. We thus predict an observable anisotropy in the composition at high energy in the sense that lighter nuclei should preferentially be seen toward the general direction of Cen A. Taking into consideration the magnetic field models for the Galactic disk and a Galactic magnetic wind, this scenario may resolve the discrepancy between HiRes and Auger results concerning the chemical composition of UHECRs.

We selected a sample of 24 XMM-Newton light curves (LCs) of four high energy peaked blazars, PKS 0548-322, ON 231, 1ES 1426+428 and PKS 2155-304. These data comprise continuous light curves of 7.67h to 18.97h in length. We searched for possible quasi-periodic oscillations (QPO) and intra-day variability (IDV) timescales in the LCs of these blazars. We found a likely QPO in one LC of PKS 2155-304 which was reported elsewhere (Lachowicz et al. 2009). In the remaining 23 LCs we found hints of possible weak QPOs in one LC of each of ON 231 and PKS 2155-304, but neither is statistically significant. We found IDV timescales that ranged from 15.7 ks to 46.8 ks in 8 LCs. In 13 LCs any variability timescales were longer than the length of the data. Assuming the possible weak QPO periods in the blazars PKS 2155-304 and ON 231 are real and are associated with the innermost portions of their accretion disk, we can estimate that their central black hole masses exceed 1.2 $\times$ 10$^{7}$ M$_{\odot}$. Emission models for radio-loud active galactic nuclei (AGN) that could explain our results are briefly discussed.

We obtained SDSS spectra for a set of 37 radio-quiet quasars (RQQSOs) that had been previously examined for rapid small scale optical variations, or microvariability. Their H-beta and Mg II emission lines were carefully fit to determine line widths (FWHM) as well as equivalent widths (EW) due to the broad emission line components. The line widths were used to estimate black hole masses and Eddington ratios, ell. Both EW and FWHM are anticorrelated with ell. The EW distributions provide no evidence for the hypothesis that a weak jet component in the RQQSOs is responsible for their microvariability.

We have measured multi-band optical flux and colour variations for a sample of 12 low energy peaked blazars on short, day-to-month, timescales. Our sample contains six BL Lacertae objects and six flat spectrum radio quasars. These photometric observations, made during September 2008 to June 2009, used five optical telescopes, one in India and four in Bulgaria. We detected short term flux variations in eleven of these blazars and colour variability in eight of them. Our data indicate that six blazars (3C 66A, AO 0235+164, S5 0716+714, PKS 0735+178, OJ 287 and 3C 454.3) were observed in pre- or post-outburst states, that five (PKS 0420-014, 4C 29.45, 3C 279, PKS 1510-089 and BL Lac) were in a low state, while one (3C 273) was in an essentially steady state. The duty cycles for flux and colour variations on short timescales in these low energy peaked blazars are ~ 92 percent and ~ 33 percent, respectively. The colour vs magnitude correlations seen here support the hypothesis that BL Lac objects tend to become bluer with increase in brightness; however, flat spectrum radio quasars may show the opposite trend, and there are exceptions to these trends in both categories of blazar. We briefly discuss emission models for active galactic nuclei that might explain our results.

We report the results of our intensive intranight optical monitoring of 8 `radio-intermediate quasars' (RIQs) having flat or inverted radio spectra. The monitoring was carried out in {\it R-} band on 25 nights during 2005-09. An intranight optical variability (INOV) detection threshold of $\sim$ 1--2% was achieved for the densely sampled differential light curves (DLCs). These observations amount to a large increase over those reported hitherto for this rare and sparsely studied class of quasars which can, however, play an important role in understanding the link between the dominant varieties of powerful AGN, namely the radio-quiet quasars (RQQs), radio-loud quasars (RLQs) and blazars. Despite the probable presence of relativistically boosted nuclear jets, clear evidence for INOV in our extensive observations was detected only on one night. These results demonstrate that as a class, RIQs are much less extreme in nuclear activity compared to blazars. The availability in the literature of INOV data for another 2 RIQs conforming to our selection criteria allowed us to enlarge the sample to 10 RIQs (monitored on a total of 42 nights for a minimum duration of $\sim 4$ hours per night). The absence of large amplitude INOV $(\psi > 3%)$ persists in this enlarged sample. This extensive database has enabled us to arrive at the first estimate for the INOV Duty Cycle (DC) of RIQs. The DC is found to be small ($\sim$ 9%). The corresponding value is known to be $\sim 60%$ for BL Lacs and $\approx 15%$ for RLQs and RQQs. On longer-term, the RIQs are found to be fairly variable with typical amplitudes of $\approx$ 0.1-mag. The light curves of these RIQs are briefly discussed in the context of a theoretical framework proposed earlier for linking this rare kind of quasars to the much better studied dominant classes of quasars.

We report a possible detection of an ~4.6-hour quasi-periodic oscillation (QPO) in the 0.3-10 keV emission of the high-energy peaked blazar PKS 2155-304 from a 64 ks observation by the XMM-Newton EPIC/pn detector. We identify a total modulation of ~5% in the light curve and confirm that nominal period by periodogram, structure function and wavelet analyses. The limited light curve duration allows the capture of only 3.8 cycles of this oscillation and thus precludes a very strong claim for this QPO, despite a nominally high (>3 sigma) statistical significance. We briefly discuss models capable of producing an X-ray QPO of such a period in a blazar.

We withdraw our claim that a component in an XMM-Newton satellite light curve of the BL Lacertae object S5 0716 + 714 shows quasi-periodic oscillations (QPOs) of $\sim$30 minutes. Although both our original periodogram and wavelet analyses gave consistent results, the data do not lead to a statistically significant result once red-noise at low frequencies is properly taken into account for periodogram analyses.

The recently discovered apparent dramatic expansion in the effective radii of massive elliptical galaxies from $z \simeq 2$ to $z \simeq 0.1$ has been interpreted in terms of either galaxy mergers or the rapid loss of cold gas due to AGN feedback. In examining the latter case we have quantified the extent of the expansion, which is uncertain observationally, in terms of the star formation parameters and time of the expulsion of the cold gas. In either case, the large global decrease in stellar density should translate into a major drop in the ISM density and pressure, and a much steeper radial decline in those quantities with cosmic epoch. These cosmological changes are expected to have a major influence on the gas accretion mode, which will shift from `cold' thin disk accretion at high redshifts toward `hot' Bondi fed ADAF accretion at low redshifts. The decline of angular momentum inflow would then lead to a spin down of the black hole, for which we have calculated more precise time scales; a value of about 0.2 Gyr is typical for a $10^9 M_\odot$ central black hole. These results have implications for the different cosmological evolutionary patterns found for the luminosity functions of powerful and weak radio galaxies.

We present low-frequency observations with the GMRT of three giant radio sources (J0139+3957, J0200+4049 and J0807+7400) with relaxed diffuse lobes which show no hotspots and no evidence of jets. The largest of these three, J0200+4049, exhibits a depression in the centre of the western lobe, while J0139+3957 and J0807+7400 have been suggested earlier by Klein et al. and Lara et al. respectively to be relic radio sources. We estimate the spectral ages of the lobes. All three sources have compact radio cores. Although the radio cores suggest that the sources are currently active, we suggest that the lobes in these sources could be due to an earlier cycle of activity.

We have performed a structure function analysis of the Rossi X-ray Timing Explorer All Sky Monitor data to search for variability in 24 blazars using data trains that each exceed 12 years. Although 20 of them show nominal periods though this technique, the great majority of these `periods' are clearly related to yearly variations arising from the instrument.Nonetheless, an apparently real periodic component of about 17 days was detected for the blazar AO 0235+164 and it was confirmed by discrete correlation function and periodogram analyses. For 1ES 2321+419 a component of variability with a near periodicity of about 420 days was detected by all of these methods. We discuss several possible explanations for these nearly periodic components and conclude that they most likely arise from the intersections of a shock propagating down a relativistic jet that possesses a helical structure.

Quasi-Periodic Oscillations (QPOs) are very puzzling since they remain totally unexplained by popular earlier models of accretion disks. The significant rms value in power density spectrum implies that the oscillation involves in the dynamical and non-linear variation of certain region of the accretion disk itself. The nature of the energy dependence implies that the region which produces Comptonized hard tail is also responsible for QPOs. Similarly, the occurrences of the QPOs are strongly related to the jet formation and the spectral states. These features are the natural consequences of the advective disk paradigm that we are advocating. In the mid 90s, some of the present authors first pointed out that the QPOs in all possible types of black holes may be simply due to the oscillations of the CENBOL, the CENtrifugal pressure supported BOundary Layer which is formed in the sub-Keplerian flows around a black hole. This CENBOL could be axi-symmetric as well as non-axisymteric in nature since its boundary, namely, the centrifugally driven shocks could be axi-symmetric or non-axisymmetric. In addition, we pointed out that the transition radius where the flow becomes Keplerian to sub-Keplerian, as well as the location of the inner sonic point can also oscillate and produce the QPOs. Since the shock locations are functions of the specific angular momentum ($\lambda$) and specific energy (${\cal E}$) of the flow, our model naturally predicts that the QPO frequency should vary with mass, spin, $\lambda$ and ${\cal E}$. The QPO frequencies with specific ratios, such as, 2:3 must be due to non-axisymmetric effects when the shock switches between the two-armed and the three-armed spirals. We also discuss the possible effects that the disk inclination might have with the line of sight.

This work introduces a novel formalism to investigate the role of the spin of astrophysical black holes in determining the behaviour of matter falling onto such accretors. Equations describing the general relativistic hydrodynamic accretion flow in the Kerr metric are formulated, and stationary solutions for such flow equations are provided. The accreting matter may become multi-transonic, allowing a stationary shock to form for certain initial boundary conditions. Such a shock determines the disc geometry and can drive strong outflows. The properties of matter extremely close to the event horizon are studied as a function of the Kerr parameter, leading to the possibility of detecting a new spectral signature of black hole spin.

We present low-frequency observations with the Giant Metrewave Radio Telescope (GMRT) of three giant radio sources (GRSs, J0139+3957, J0200+4049 and J0807+7400) with relaxed diffuse lobes which show no hotspots and no evidence of jets. The largest of these three, J0200+4049, exhibits a depression in the centre of the western lobe, while J0139+3957 and J0807+7400 have been suggested earlier by Klein et al. and Lara et al. respectively to be relic radio sources. We estimate the ages of the lobes. We also present Very Large Array (VLA) observations of the core of J0807+7400, and determine the core radio spectra for all three sources. Although the radio cores suggest that the sources are currently active, we explore the possibility that the lobes in these sources are due to an earlier cycle of activity.

We describe a method for incorporating ambipolar diffusion in the strong coupling approximation into a multidimensional magnetohydrodynamics code based on the total variation diminishing scheme. Contributions from ambipolar diffusion terms are included by explicit finite difference operators in a fully unsplit way, maintaining second order accuracy. The divergence-free condition of magnetic fields is exactly ensured at all times by a flux-interpolated constrained transport scheme. The super time stepping method is used to accelerate the timestep in high resolution calculations and/or in strong ambipolar diffusion. We perform two test problems, the steady-state oblique C-type shocks and the decay of Alfv\'en waves, confirming the accuracy and robustness of our numerical approach. Results from the simulations of the compressible MHD turbulence with ambipolar diffusion show the flexibility of our method as well as its ability to follow complex MHD flows in the presence of ambipolar diffusion. These simulations show that the dissipation rate of MHD turbulence is strongly affected by the strength of ambipolar diffusion.

The question as to whether the distribution of radio-loudness in active galactic nuclei (AGN) is actually bimodal has been discussed extensively in the literature. Futhermore, there have been claims that radio-loudness depends on black hole mass and Eddington ratio. We investigate these claims using the low redshift broad line AGN sample of Greene & Ho (2007), which consists of 8434 objects at z < 0.35 from the Sloan Digital Sky Survey Fourth Data Release (SDSS DR4). We obtained radio fluxes from the Very Large Array Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) survey for the SDSS AGN. Out of the 8434 SDSS AGN, 821 have corresponding observed radio fluxes in the FIRST survey. We calculated the radio-loudness parameter (R) for all objects above the FIRST detection limit (1 mJy), and an upper limit to R for the undetected objects. Using these data, the question of radio bimodality is investigated for different subsets of the total sample. We find no clear demarcation between the radio-loud (RL, R > 10) and radio-quiet (RQ, R < 10) objects, but instead fill in a more radio-intermediate population in a continuous fashion for all subsamples. We find that 4.7% of the AGN in the flux-limited subsample are RL based on core radio emission alone. We calculate the radio-loud fraction (RLF) as both a function of black hole mass and Eddington ratio. The RLF decreases (from 13% to 2%) as Eddington ratio increases over 2.5 order of magnitude. The RLF is nearly constant (~5%) over 4 decades in black hole mass, except for an increase at masses greater than 10^8 solar masses. We find for the FIRST detected subsample that 367 of the RL AGN have black hole masses less than 10^8 solar masses, a large enough number to indicate that RL AGN are not a product of only the most massive black holes in the local universe.

We present results of a periodicity search of 20 intra-day variable optical light curves of the blazar S5 0716+714, selected from a database of 102 light curves spanning over three years. We use a wavelet analysis technique along with a randomization test and find strong candidates for nearly periodic variations in eight light curves, with probabilities ranging from 95% to >99%. This is the first good evidence for periodic, or more-precisely, quasi-periodic, components in the optical intra-day variable light curves of any blazar. Such periodic flux changes support the idea that some active galactic nuclei variability, even in blazars, is based on accretion disk fluctuations or oscillations. These intra-day variability time scales are used to estimate that the central black hole of the blazar S5 0716+714 has a mass > 2.5 \times 10^6$ M$_{\odot}$. As we did not find any correlations between the flux levels and intra-day variability time scales, it appears that more than one emission mechanism is at work in this blazar.

We present a sample of 16 radio galaxies, each of which is characterized by a wide, elongated emission gap with fairly sharp and straight edges between the two radio lobes. This particular subset of the 'superdisk' radio galaxies is chosen because of a highly asymmetric location of the host elliptical galaxy relative to the gap's central axis. In addition to posing a considerable challenge to the existing models, such a morphology also means that the two jets traverse highly unequal distances through the superdisk material. One thus has a possibility to directly investigate if the marked asymmetry between the two jets' interaction with the (much denser) ambient medium, during their propagation, has a significant import for the brightness of the hotspot forming near each jet's extremity. We also propose a new explanation for the formation of superdisks through the merger of a smaller elliptical galaxy with the massive host, in which the gas attached to the infalling galaxy deposits its angular momentum into the host's circumgalactic gas, thereby causing it to flatten into a fat pancake, or superdisk. The asymmetric location of the host galaxy can be assisted by the kick imparted to it during the merger. We also suggest a physical link between these radio galaxies and those with X-shaped and Z-symmetric radio lobes, commonly believed to arise from mergers of two galactic nuclei, each harboring a supermassive black hole.

The origin of the dichotomy of radio loudness among quasars can be explained using recent findings that the mass of the central supermassive black hole (SMBH) in extended radio-loud quasars is systematically a few times that of their counterparts in radio-quiet quasars. This sensitive dependence of radio jet ejection upon SMBH mass probably arises from the blockage of jets by the presence of substantial quantities of gas tidally stripped from stars by the central BH. This disruptive gas, however, will only be available around BHs with masses less than $M_c ~\gtrsim ~10^8\Msun$, for which the tidal disruption radius lies outside the SMBH's event horizon. Consequently, we find that AGN with $M_{BH} > M_c$ can successfully launch jets with a wide range of powers, thus producing radio-loud quasars. The great majority of jets launched by less massive BHs, however, will be truncated in the vicinity of the SMBH due to mass loading from this stellar debris. This scenario also can naturally explain the remarkable dearth of extended radio structures in quasars showing broad absorption line spectra.

Taking a clue from their sharp-edged (strip-like) morphology observed in several cases, a new mechanism is proposed for the formation of the emission gaps seen between the radio lobes of many powerful extragalactic double radio sources. Canonical understanding of the radio gaps invokes either blocking of the back-flowing lobe plasma by the denser interstellar medium (ISM) of the host galaxy, or "squeezing" of the radio bridge in the middle through buoyancy force exerted by either the ISM or the surrounding intra-cluster medium (ICM). These pictures encounter difficulties in explaining situations where the sharp-edged radio gaps associated with non-cluster radio galaxies have widths running into several tens (even hundreds) of kiloparsecs. More particularly, the required dense high-pressure ISM/ICM is likely to be lacking at least in the case of high-redshift radio galaxies. We propose here that radio emission gaps in at least such cases could arise from a dynamical interaction between the powerful thermal wind outflowing from the active galactic nucleus and the back-flowing synchrotron plasma in the two radio lobes, which occurs once the rapidly advancing jets have crossed out of the wind zone into the intergalactic medium. A simple analytical scheme is presented to explore the plausibility of the side-ways confinement of the thermal wind by the radio lobe pair, which would "freeze" pancake shaped conduits in the space, along which the hot, metal enriched wind from the AGN can escape (roughly orthogonal to the radio axis). Some other possible consequences of this scenario are pointed out.

We present analytical modelling of conical relativistic jets, in order to evaluate the role of the jet opening angle on certain key parameters that are inferred from VLBI radio observations of blazar nuclear jets. The key parameters evaluated are the orientation angle (i.e., the viewing angle) of the jet and the apparent speed and Doppler factor of the radio knots on parsec scales. Quantitative comparisons are made of the influence of the jet opening angle on the above parameters of the radio knots, as would be estimated for two widely discussed variants of relativistic nuclear jets, namely, those having uniform bulk speed and those in which the bulk Lorentz factor of the flow decreases with distance from the jet axis (a `spine--sheath' flow). Our analysis shows that for both types of jet velocity distributions the expectation value of the jet orientation angle at first falls dramatically with increases in the (central) jet Lorentz factor, but it levels off at a fraction of the opening angle for extremely relativistic jets. We also find that the effective values of the apparent speeds and Doppler factors of the knots always decline substantially with increasing jet opening angle, but that this effect is strongest for ultra-relativistic jets with uniform bulk speed. We suggest that the paucity of highly superluminal parsec-scale radio components in TeV blazars can be understood if their jets are highly relativistic and, being intrinsically weaker, somewhat less well collimated, in comparison to the jets in other blazars.

We had earlier shown that for a constant velocity jet the discrepancy between the low speeds indicated by VLBI knot motions and the high Doppler factors inferred from emission of TeV photons could be reconciled if ultrarelativistic jets possessed modest opening angles. Here we evaluate the (flux-weighted) viewing angles of the jet and the apparent speeds and Doppler factors of the radio knots on parsec scales. The influence of the jet opening angle on these radio knot parameters are found for the usually considered types of relativistic nuclear jets: those with uniform bulk speeds and those where the bulk Lorentz factor of the flow decreases with distance from the jet axis, known as `spine--sheath' flows. For both types of jet velocity structures the expectation value of the jet orientation angle at first falls dramatically with increases in the (central) jet Lorentz factor, but for extremely relativistic jets it levels off at a fraction of the opening angle. The effective values of the apparent speeds and Doppler factors of the knots always decline substantially with increasing jet opening angle. The rarity of highly superluminal parsec-scale radio components in TeV blazars can be understood if their jets are both highly relativistic and intrinsically weaker, so probably less well collimated, than the jets in ordinary blazars.

We investigate aspects of the cosmological evolution of FR II radio galaxies, focusing first on the abilities of models to match data for linear-sizes, radio powers, redshifts and spectral indices. Here we consider modifications to the theoretical models we had treated earlier, primarily by accounting for the growth of the radius of hotspots with source size. Better fits to the distributions of most of the data in three low-frequency surveys can be found with sensible choices of model parameters but no model yet considered gives a good match to all of the survey data simultaneously, nor does any do a good job of producing the spectral index distributions. The observational datasets are too small to completely discriminate among the models. We calculate the volume fraction of the ``relevant universe'' cumulatively occupied by the expanding radio galaxy lobes over the quasar era, when these powerful radio galaxies were much more common, and when they have been argued to play an important role in triggering galaxy formation and spreading magnetic fields and metals. We found the cumulative relevant volume filling factor of radio galaxies to be $\sim 5 %$, so we conclude that these impacts are smaller than previously estimated but that they are still significant.

We have studied three-dimensional hydrodynamic interactions of relativistic extragalactic jets with two-phase ambient media. These jets propagate through a denser homogeneous gas and then impact clouds with densities 100 to 1000 times higher than the initial beam density. The deflection angle of the jet is influenced more by the density contrast of the cloud than by the beam Mach number of the jet. A relativistic jet with low relativistic beam Mach number can eventually be slightly bent after it crosses the dense cloud; however, we have not seen permanently bent structures in the interaction of a high relativistic beam Mach number jet with a cloud. The relativistic jet impacts on dense clouds do not necessarily destroy the clouds completely, and much of the cloud body can survive as a coherent blob. This enhancement of cloud durability is partly due to the geometric influence of the off-axis collisions we consider and also arises from the lower rate of cloud fragmentation through the Kelvin-Helmholtz instability for relativistic jets. To compare our simulations with observed extragalactic radio jets, we have computed the approximate surface distributions of synchrotron emission at different viewing angles. These surface intensity maps show relativistic jets interacting with clouds can produce synchrotron emission knots similar to structures observed in many VLBI-scale radio sources. We find that the synchrotron emission increases steeply at the moment of impact and the emission peaks right before the jet passes through the cloud.

We seek to develop an essentially analytical model for the evolution of Fanaroff-Riley Class II radio galaxies as they age individually and as their numbers vary with cosmological epoch. Such modeling is required in order to probe in more detail the impact of radio galaxies on the growth of structures in the universe, which appears likely to have been quite significant at z > 1. In this first paper of a series we compare three rather sophisticated analytical models for the evolution of linear size and lobe power of FR II radio galaxies, those of Kaiser, Dennett-Thorpe & Alexander (1997), Blundell, Rawlings & Willott (1999), and Manolakou & Kirk (2002). We perform multi-dimensional Monte Carlo simulations in order to compare the predictions of each model for radio powers, sizes, redshifts and spectral indices with data. The observational samples used here are the low frequency radio surveys, 3CRR, 6CE and 7CRS, which are flux-limited and complete. We search for and describe the best parameters for each model, after doing statistical tests on them. We find that no existing model can give acceptable fits to all the properties of the surveys considered, although the Kaiser, Dennett-Thorpe & Alexander (1997) model gives better overall results than do the Manolakou & Kirk (2002) or Blundell, Rawlings & Willott (1999) models for most of the tests we performed. We suggest ways in which these models may be improved.

We first review some basic properties of the most important classes of active galactic nuclei (AGN), including quasars, blazars, Seyfert galaxies and radio galaxies. The most commonly accepted type of scheme designed to individually unify the radio-loud and radio-quiet categories of AGN is based upon three parameters: black hole mass, accretion rate, and our orientation to the accretion flow. Some recent evidence from optical microvariability of several classes of AGN points in favor of a strong unification scheme that unites both radio-loud and radio-quiet categories. An important question concerning the nature of blazars and other AGN whose jet emission appears to dominate their spectral energy distributions involves the velocities of those flows. A variety of apparently contradictory observations can be reconciled if such flows are ultrarelativistic but have an opening angle of a few degrees. Radio galaxies (RGs) were much more numerous at redshifts ~2 than they are today. Combining this fact with the realization that older RGs at such redshifts are very difficult to detect, and with cosmological simulations of the growth of structure in the universe has led us to propose that RG lobes have impacted a significant fraction of the cosmic web of baryons. These impacts may have triggered extensive star formation and perhaps even engendered new galaxies; they also probably played important roles in the spreading of magnetic fields and heavier elements into the intergalactic medium.

Allowing for the conical shape of ultrarelativistic blazar jets with opening angles of a few degrees on parsec-scales we show that their bulk Lorentz factors and viewing angles can be much larger than the values usually inferred by combining their flux variability and proper motion measurements. This is in accord with our earlier finding that such ultrarelativistic (Lorentz factor, Gamma > 30) conical jets can reconcile the relatively slow apparent motions of VLBI knots in TeV blazars with the extremely fast flows implied by their rapid gamma-ray variability. This jet geometry also implies that de-projected jet opening angles will typically be significantly underestimated from VLBI measurements. In addition, de-projected jet lengths will be considerably overestimated if high Lorentz factors and significant opening angles are not taken into account.

We report results of multiband optical monitoring of two well known blazars, S5 0716+714 and BL Lacertae, carried out in 1996 and 2000-01 with an aim to study optical variations on time scales from minutes to hours and longer.The light curves were derived relative to comparison stars present on the CCD frames. Night to night flux variations of >0.1 mag were observed in S5 0716+714 during a campaign of ~2 weeks in 1996.A good correlation between the lightcurves in different optical bands was found for both inter-night and intra-night observations. Two prominent events of intra-night optical variability were detected in S5 0716+714.Each of these rapidly varying segments of the lightcurves trace an exponential flux profile whose rate of variation is the same in both cases. Our long-term monitoring data of S5 0716+714 showed a distinct flare around JD 2451875 which can be identified in the BVRI bands.This flare coincides with the brightest phase recorded during 1994-2001 in the long-term lightcurves reported by Raiteri et al.(2003). No evidence for the bluer when brighter trend was noticed on inter-night and intra-night time scales. On the other hand, our nearly simultaneous multiband observations of BL Lacertae in October 2001 showed flux variations that were not achromatic. This blazar was found to become bluer when brighter on intra-night time scales and there is a hint of the same trend on inter-night time scales. Based on five nights of observations during a week, BL Lacertae showed a peak night-to-night variability of ~0.6 mag in B. Thus, we found that the present observations of the two blazars, reveal a contrasting behaviour in terms of the dependence of spectral hardening with increasing brightness, at least on intra-night, and possibly also on inter-night, time scales.

Although blazar variability is probably dominated by emission from relativistic jets, accretion disks should be present in all blazars. These disks produce emission over most of the electromagnetic spectrum; various unstable processes operate in those disks which lead to variable emission. Here I summarize some of the most relevant disk mechanisms for AGN variability. I also discuss some aspects of jet variability, focusing on the possibility that ultrarelativisitic jets of modest opening angle can reconcile TeV blazar emission with the many subluminal VLBI knots seen in those sources. Finally, I present recently illuminated characteristics of optical microvariability of different classes of AGN which have important implications for the dominant processes involved.

As Compact Steep Spectrum radio sources have been shown to be more asymmetrical than larger sources of similar powers, there is a high probability that they interact with an asymmetric medium in the central regions of the host elliptical galaxy. We consider a simple analytical model of the propagation of radio jets through a reasonable asymmetric environment and show that they can yield the range of arm-length and luminosity asymmetries that have been observed. We then generalize this to allow for the effects of orientation, and quantify the substantial enhancements of the asymmetries that can be produced in this fashion. We present two-dimensional and three-dimensional simulations of jets propagating through multi-phase media and note that the results from the simulations are also broadly consistent with the observations.

We compare earlier estimates of the volumes filled by lobes of radio galaxies during the quasar era based upon non-evolving ambient media with new ones assuming a strong cosmological evolution of the ambient medium. If the sources remain active for over 10^8 years the volumes filled by them are found to be comparable for the two scenarios. This strengthens our earlier inference that much of the cosmic web of gaseous filaments, the site of galaxy formation, was probably permeated by radio lobes during the quasar era and this could have triggered extensive star formation and made large contributions to the spread of magnetic fields and metals through the universe by z ~ 2.

This review provides the background and current perspective of some of the striking asymmetries exhibited by extragalactic double-lobed radio sources, paying special attention to the reported interrelationships among some of them. We particularly aim to highlight some less discussed aspects pertaining to the radio lobes and hot-spots for which any relativistic beaming effects are milder than those imposed on the jets. Various attributes of these features, such as their geometry, intensity, spectrum, and polarization provide important underpinnings to models of radio galaxies and quasars, as well as to theories that attempt to unify them. These contraints are sharpened by exploiting the asymmetry properties of the extended optical emission-line regions associated with powerful radio sources, including the giant halos of Lyman-alpha emission associated with them. We suggest that highly extended disks of dusty gaseous material, which could also be parts of cosmic filaments, around the host galaxies of radio sources may hold the key to understanding the correlated radio--optical asymmetry and perhaps also the depolarization asymmetry of the twin radio lobes. We also discuss the nature and possible implications of the extreme type of asymmetry witnessed in the enigmatic quasar 3C273, whose entire radio emission is observed on just one side of the nucleus. In this context, theoretical scenarios for the production of one-sided jets are also briefly recounted. Taking a clue from theobserved one-sidedness of the extended radio jets even in `weak-headed-quasars', we also propose a scenario for the late evolutionary stage of double radio sources. We end with a list of key questions concerning radio source asymmetries which still require final answers. (Abridged)

We compute the effective values of apparent transverse velocity and flux boosting factors for the VLBI radio knots of blazar jets, by integrating over the angular distributions of these quantities across the widths of jets with finite opening angles but constant velocities. For high bulk Lorentz factors (Gamma > 10) variations across the jet can be quite large if the opening angle, omega, is even a few degrees on sub-parsec scales. The resulting apparent speeds are often much lower than those obtained from the usual analyses that ignore the finite jet opening angles. We can thus reconcile the usually observed subluminal or mildly superluminal speeds with the very high (>~ 20) Gamma factors, required by the inverse Compton origin and rapid variability of TeV fluxes, as well as by intraday radio variability. Thus it is possible to associate the VLBI radio knots directly with shocks in the ultra-relativistic main jet flow, without invoking very rapid jet deceleration on parsec scales, or extremely unlikely viewing angles.

We report optical monitoring observations of 20 high-luminosity AGN, 12 of which are radio-quiet quasars (RQQs). Intra-night optical variability (INOV) was detected for 13 of the 20 objects, including 5 RQQs. The variations are distinctly stronger and more frequent for blazars than for the other AGN classes. By combining these data with results obtained earlier in our program, we have formed an enlarged sample consisting of 9 BL Lacs, 19 RQQs and 11 lobe-dominated radio-loud quasars. The moderate level of rapid optical variability found for both RQQs and radio lobe-dominated quasars argues against a direct link between INOV and radio-loudness. We supplemented the present observations of 3 BL Lacs with additional data from the literature. In this extended sample of 12 well observed BL Lacs, stronger INOV is found for the EGRET detected BL Lacs.

We propose a novel approach to the investigation of how black hole spin influences the behavior of astrophysical accretion. We examine the terminal behavior of general relativistic, multi-transonic, hydrodynamic advective accretion in the Kerr metric, by analyzing various dynamical as well as thermodynamic properties of accreting matter extremely close to the event horizon, as a function of Kerr parameter $a$. The examination of the transonic properties of such flows can be useful in self-consistently discriminating between prograde and retrograde relativistic accretion, and in studying the spectral signatures of black hole rotation.

We positionally match the 6 cm GB6, 20 cm FIRST and NVSS, and 92 cm WENSS radio catalogs and find 16,500 matches in ~3,000 deg2 of sky. Using this unified radio database, we construct radio "color-magnitude-morphology" diagrams and find that they display a clear structure, rather than a random scatter. We propose a simple, yet powerful, method for morphological classification of radio sources based on FIRST and NVSS measurements. For a subset of matched sources, we find optical identifications using the SDSS Data Release 1 catalogs, and separate them into quasars and galaxies. Compact radio sources with flat radio spectra are dominated by quasars, while compact sources with steep spectra, and resolved radio sources, contain substantial numbers of both quasars and galaxies.

Although the possibility of a low energy cut-off (LEC) of the relativistic electron population (E_{min} ~ 0.1 GeV) in the different components of radio galaxies has been discussed in the literature, from both theoretical and observational viewpoints, its possible role in causing a distinct spectral flattening, and thus reducing the apparent brightness temperature of the radio nucleus, has thus far not been explicitly discussed. Here we point out that such an effect may in fact be quite significant for the parsec-scale, superluminal emission regions associated with the inner radio jets. This is because the observed frequency of the spectral turnover arising from the LEC increases linearly with the bulk Doppler factor (\delta) of the jet flow, whereas the frequency of synchrotron self-absorption (the well known cause of spectral flattening and turnover at low frequencies) shows a mild decrease as \delta rises. For an observationally relevant range of source parameters, we investigate here the role of the LEC. We argue that some statistical trends revealed by the recent VLBI surveys can in fact be understood in terms of effects arising from LEC, including the faster superluminal motion found in the VLBI surveys at higher frequencies, and the apparent inability of most radio cores to even attain the brightness temperatures consistent with the equipartition condition. We also point out some possible implications of the LEC for X-ray observations of the kiloparsec scale relativistic jets of quasars.

We present results of a multi-epoch intra-night optical monitoring of eleven blazars consisting of six BL Lac objects and five radio core-dominated quasars (CDQs). These densely sampled and sensitive R-band CCD observations, carried out from November 1998 through May 2002 during a total of 47 nights with an average of 6.5 hours/night, have enabled us to detect variability amplitudes as low as $\sim 1%$ on intra-night time scales. A distinction is found for the first time between the intra-night optical variability (INOV) properties of the these two classes of relativistically beamed radio-loud AGNs. BL Lacs are found to show a duty cycle (DC) of INOV of $\sim$60%, in contrast to CDQs, which show a much smaller INOV DC of $\sim$20%, the difference being attributable mainly to the weakly polarized CDQs. On longer time scales (i.e., between a week to a few years) variability is seen from all the CDQs and BL Lacs in our sample. The results reported here form part of our long-term programme to understand the intra-night optical variability characteristics of the four main classes of luminous AGNs, i.e., radio-quiet quasars (RQQs) and radio lobe-dominated quasars (LDQs), as well as CDQs and BL Lac objects.

We present an investigation of new optical and ultraviolet spectra of the mass donor star in the massive X-ray binary LS 5039. The optical band spectral line strengths indicate that the atmosphere is N-rich and C-poor, and we classify the stellar spectrum as type ON6.5 V((f)). The N-strong and C-weak pattern is also found in the stellar wind P Cygni lines of N V 1240 and C IV 1550. We suggest that the N-enrichment may result from internal mixing if the O-star was born as a rapid rotator, or the O-star may have accreted N-rich gas prior to a common-envelope interaction with the progenitor of the supernova. We re-evaluated the orbital elements to find an orbital period of P=4.4267 +/- 0.0010 d. We compared the spectral line profiles with new non-LTE, line-blanketed model spectra, from which we derive an effective temperature T_eff = 37.5 +/- 1.7 kK, gravity log g = 4.0 +/- 0.1, and projected rotational velocity V sin i = 140 +/- 8 km/s. We fit the UV, optical, and IR flux distribution using a model spectrum and extinction law with parameters E(B-V)= 1.28 +/- 0.02 and R= 3.18 +/- 0.07. We confirm the co-variability of the observed X-ray flux and stellar wind mass loss rate derived from the H-alpha profile, which supports the wind accretion scenario for the X-ray production in LS 5039. Wind accretion models indicate that the compact companion has a mass M_X/M_sun = 1.4 +/- 0.4, consistent with its identification as a neutron star. The observed eccentricity and runaway velocity of the binary can only be reconciled if the neutron star received a modest kick velocity due to a slight asymmetry in the supernova explosion (during which >5 solar masses was ejected).

Using the surface photometric analysis data for a well defined sample of 79 nearby radio galaxies (z < 0.12), it is shown that the level of radio emission associated with an elliptical galaxy during the radio-loud phase is related to the presence of large amounts of isophotal twist in its optical image. In particular, radio galaxies with P_408 > 3 x 10^{25} W/Hz appear to show a preference to be associated with elliptical hosts exhibiting an isophotal twist in excess of ~20 degrees.

Existing radio images of a few X-shaped radio galaxies reveal Z-symmetric morphologies in their weaker secondary lobes which cannot be naturally explained by either the galactic merger or radio-lobe backflow scenarios, the two dominant models for these X-shaped radio sources. We show that the merger picture can explain these morphologies provided one takes into account that, prior to the coalescence of their supermassive black holes, the smaller galaxy releases significant amounts of gas into the ISM of the dominant active galaxy. This rotating gas, whose angular momentum axis will typically not be aligned with the original jets, is likely to provide sufficient ram pressure at a distance ~10 kpc from the nucleus to bend the extant jets emerging from the central engine, thus producing a Z-symmetry in the pair of radio lobes. Once the two black holes have coalesced some 10^7 yr later, a rapid reorientation of the jets along a direction close to that of the orbital angular momentum of the swallowed galaxy relative to the primary galaxy would create the younger primary lobes of the X-shaped radio galaxy. This picture naturally explains why such sources typically have powers close to the FR I/II break. We suggest that purely Z-symmetric radio sources are often en route to coalescence and the concomitant emission of substantial gravitational radiation, while X-shaped ones have already merged and radiated.

We present results of a programme of multi-epoch, intra-night optical monitoring of a sample of non-blazar type AGN, which includes seven radio-quiet QSOs (RQQs) and an equal number of radio-loud, lobe-dominated quasars (LDQs), covering a redshift range from about 0.2 to 2.0. These two sets of optically bright and intrinsically luminous QSOs are well matched in the redshift-optical luminosity ($z - M_B$) plane. Our CCD monitoring covered a total of 61 nights with an average of 6.1 hours of densely sampled monitoring of just a single QSO per night, thereby achieving a typical detection threshold of $\sim 1$% variation over the night. Unambiguous detection of intra-night variability(INOV) amplitudes in the range 1$-$3% on day-like or shorter time scales were thus made for both RQQs and LDQs. Based on these clear detections of INOV, we estimate duty cycles of 17% and 9% for RQQs and LDQs, respectively; inclusion of the two cases of probable variations of LDQs would raise the duty cycle to 15% for LDQs. The similarity in the duty cycle and amplitude of INOV for the RQQs and LDQs suggests, firstly, that the radio loudness alone does not guarantee an enhanced INOV in QSOs and,secondly, that as in LDQs, relativistic jets may also be present in RQQs. We argue that, as compared to BL Lacs, the conspicuously milder, rarer and possibly slower INOV of RQQs and LDQs, can in fact be readily understood in terms of their having optical synchrotron jets which are modestly misaligned from us,but are otherwise intrinsically as relativistic and active as the jets in BL Lacs.

We present the results of a comparative study of the intra-night optical variability (INOV) characteristics of radio-loud and radio-quiet quasars, which involves a systematic intra-night optical monitoring of seven sets of high luminosity AGNs covering the redshift range {\it z} $\simeq 0.2$ to {\it z} $\simeq 2.2$. The sample, matched in the optical luminosity -- redshift (M$_B$ -- z) plane, consists of seven radio-quiet quasars (RQQs), eight radio lobe-dominated quasars (LDQs), six radio core-dominated quasars (CDQs) and five BL Lac objects (BLs). Systematic CCD observations, aided by a careful data analysis procedure, have allowed us to detect INOV with amplitudes as low as 1%. Present observations cover a total of 113 nights (720 hours) with only a single quasar monitored as continuously as possible on a night. Considering cases of only unambiguous detections of INOV we have estimated duty cycles (DCs) of 17%, 12%, 20% and 72% respectively for RQQs, LDQs, CDQs, and BLs. The low amplitude and low DC of INOV shown by RQQs compared to BLs can be understood in terms of their having optical synchrotron jets which are modestly misdirected from us. From our fairly extensive dataset, no unambiguous general trend of a correlation between the INOV amplitude and the apparent optical brightness of the quasar is noticed.

The compact steep spectrum (CSS) and gigahertz peaked spectrum (GPS) sources are widely believed to be young radiosources, with ages <10^6 yr. If the activity in the nucleus is fuelled by the supply of gas, one might find evidence of this gas by studying the structural and polarisation characteristics of CSS sources and their evolution through this gas. In this paper we discuss some of the possible `smoking-gun' evidence of this gas which may have triggered and fuelled the radio source.

We present new clues to the problem of the radio loudness dichotomy arising from an extensive search for intranight optical variability in seven sets of optically luminous radio-quiet quasars and (radio-loud) BL Lacertae objects, which are matched in optical luminosity and redshift. Our monitoring of radio-quiet quasars has for the first time clearly detected such intranight variability, with peak-to-peak amplitudes ~1%, occurring with a duty cycle of ~ 1/6. The matched BL Lacs have both higher variability amplitudes and duty cycles when observed in the same fashion. We show that the much less pronounced intranight variability of the radio-quiet quasars relative to BL Lacs can be understood in terms of a modest misalignment of the jets in radio-quiet quasars from the line-of-sight. We thus infer that relativistic particle jets may well also emerge from radio-quiet quasars, but while traversing the short optical-emitting distances, they could be snuffed out, possibly through inverse Compton losses in the nuclear region.

Observed radio galaxies had a much higher comoving density during the `quasar era', at z ~ 2-3, but these sources are only detectable for small fractions of their active lifetimes at such high z due to expansion losses and increased inverse Compton losses against the cosmic microwave background. Using recent models for the evolution of the size and luminosity of powerful double radio sources, as well as LCDM simulations of the cosmic web of baryonic material, we argue that during the quasar era a high volume fraction of this web was occupied by the lobes of double radio sources. They could have seeded the IGM with an average magnetic field approaching 10^{-8} G. Further, these advancing overpressured lobes could compress the denser interstellar gas clouds of the galaxies engulfed by them and thus trigger starbursts. This can probably account for much of the intense star-formation activity witnessed beyond z ~ 1.5. Also, the sweeping up of the ISM of the gas-rich galaxies by the rapidly advancing radio lobes may well be responsible for the widespread metal pollution of the IGM and proto-galaxies at high redshifts.

Multi-wavelength observations made in the last decade suggest that the universe underwent an intense phase of star formation in the past (z > 1). This intensive activity is commonly attributed to a higher galaxy merger rate when the universe was a fraction of its present age. We examine the role of the powerful radio sources whose comoving density is known to be a few orders of magnitude higher at z ~ 2, the `quasar era'. Taking into account recent models for the temporal evolution of the size and luminosity of a powerful double radio source, as well as Lambda-CDM simulations of the cosmic web of baryonic material at different redshifts, we argue that during the quasar era a high fraction of the volume of the web was occupied by the lobes of double radio sources. Widespread compression of protostellar clouds, triggered by the high pressure of the synchrotron plasma of the radio lobes, can thus be expected to have played a significant role in the global star formation history of the universe. These lobes can also yield a rather high level of magnetization of the intergalactic medium at these early cosmic epochs.

We have investigated the phenomenon of explosive chromospheric evaporation from an accretion disk as a mechanism for fast variability in accreting sources such as low mass X-ray binaries and active galactic nuclei. This has been done in the context of advection dominated accretion flows, allowing both high and low states to be considered. This mechanism can in principle produce sub-millisecond timescales in binaries and sub-minute timescales in active galaxies. However, even considering the possibility that large numbers of these microflares may be present simultaneously, the power emitted from these microflares probably amounts to only a small fraction of the total X-ray luminosity.

We present the results of a spectroscopic monitoring program (from 1998 to 2002) of the H-alpha emission strength in HDE 226868, the optical counterpart of the black hole binary, Cyg X-1. The H-alpha profiles consist of (1) a P Cygni component associated with the wind of the supergiant, (2) emission components that attain high velocity at the conjunctions and that probably form in enhanced outflows both towards and away from the black hole, and (3) an emission component that moves in anti-phase with the supergiant's motion. We argue that the third component forms in accreted gas near the black hole, and the radial velocity curve of the emission is consistent with a mass ratio of M_X / M_opt = 0.36 +/- 0.05. We find that there is a general anti-correlation between the H-alpha emission strength and X-ray flux in the sense that when the H-alpha emission is strong (W_\lambda < -0.5 Angstroms) the X-ray flux is weaker and the spectrum harder. On the other hand, there is no correlation between H-alpha emission strength and X-ray flux when H-alpha is weak. During the low/hard X-ray state, the strong wind is fast and the accretion rate is relatively low, while in the high/soft state the weaker, highly ionized wind attains only a moderate velocity and the accretion rate increases. We argue that the X-ray transitions from the normal low/hard to the rare high/soft state are triggered by episodes of decreased mass loss rate in the supergiant donor star.

We discuss plausible mechanisms to produce bullet-like ejecta from the precessing disk in the SS 433 system. We show that non-steady shocks in the sub-Keplerian accretion flow can provide the basic timescale of the ejection interval while the magnetic rubber-band effect of the toroidal flux tubes in this disk can yield flaring events.

We present results from new optical and UV spectroscopy of the unusual binary system SS 433, and we discuss the relationship of the particular spectral components we observe to the properties of the binary. (1) The continuum spectrum which we associate with flux from the super-Eddington accretion disk and the dense part of its wind. (2) H-alpha moving components which are formed far from the binary orbital plane in the relativistic jets. (3) H-alpha and He I "stationary" emission lines which we suggest are formed in the disk wind in a volume larger than the dimensions of the binary. (4) A weak "stationary" emission feature we identify as a C II 7231,7236 blend that attains maximum radial velocity at the orbital quadrature of disk recession. (5) Absorption and emission features from outflowing clumps in the disk wind (seen most clearly in an episode of blue-shifted Na I emission). (6) We found no clear evidence of the absorption line spectrum of the optical star, although we point out the presence of He I absorption features (blended with the stationary emission) with the expected radial velocity trend at the orbital and precessional phases when the star might best be seen. (7) A rich interstellar absorption spectrum of diffuse interstellar bands. The results suggest that the binary is embedded in an expanding thick disk (detected in recent radio observations) which is fed by the wind from the super-Eddington accretion disk.

Using a large set of optically detected hotspots in powerful extragalactic double radio sources, we examine the basic question of whether the detection of optical synchrotron radiation requires in-situ acceleration of relativistic electrons within the hotspots/lobes. For this, we take into account the jet's bulk relativistic motion, as well as its likely misalignment from the plane of the sky. Together, both these factors can drastically reduce the apparent range of the ultra-relativistic electrons ejected from the nucleus in the form of a jet. The conventionally adopted parameter space for the fundamental variables, namely, the hotspot magnetic field, radio source orientation angle relative to the line-of-sight and the bulk speed of the jet plasma, is considered. We find that the observed optical/near-IR synchrotron emission of the hotspots can be explained even if the radiating relativistic electrons were accelerated exclusively within the nuclear region, provided the energy losses incurred by the electrons during their transport down the jet are dominated by inverse Compton upscatterings of the cosmic microwave background photons. Under this circumstance, in-situ acceleration of relativistic electrons inside the hotspots or lobes is not found to be mandated by their reported optical/near-infrared detections.

Evidence for extended active lifetimes (> 10^8 yr) for radio galaxies implies that many large radio lobes were produced during the `quasar era', 1.5 < z < 3, when the comoving density of radio sources was 2 -- 3 dex higher than the present level. However, inverse Compton losses against the intense microwave background substantially reduce the ages and numbers of sources that are detected in flux-limited surveys. The realization that the galaxy forming material in those epochs was concentrated in filaments occupying a small fraction of the total volume then leads to the conclusion that radio lobes permeated much of the volume occupied by the protogalactic material during that era. The sustained overpressure in these extended lobes is likely to have played an important role in triggering the high inferred rate of galaxy formation at z > 1.5 and in the magnetization of the cosmic network of filaments.

We show that a model for radio source dynamics we had earlier proposed can readily reproduce the relationship between the radio power division separating the two Fanaroff-Riley classes of extragalactic radio sources and the optical luminosity of the host galaxy, as found by Owen and Ledlow (1994). In our scenario, when less powerful jets eventually slow down to the point that the advance of the working surface (i.e., hotspot) becomes subsonic with respect to the external gas, the jet's collimation is severely weakened. This criterion distinguishes the powerful and well collimated FR II sources from the weaker sources producing the less collimated FR I type morphologies.

Extragalactic radio sources, including quasars, are now typically understood as being produced by a pair of nearly symmetric, oppositely directed relativistic jets. While some these sources span megaparsecs, and are thus the largest physically connected structures in the universe, emitting regions identified as jets have now been found on all scales down to fractions of a parsec, and jets appear to be a common element of most (maybe all) types of active galactic nuclei (AGN). We first summarize key observations of different classes of cosmic radio jets, and describe how they may be connected. Theoretical models for the launching and propagation of extragalactic jets are briefly described. All of these models assume a magnetized plasma, which typically amounts to only a small fraction of the accreted gas, is ejected from the vicinity of a supermassive black hole. The extreme complexity of the relevant physics has demanded numerical simulations to examine non-linear effects on the stability of propagating jets, and some recent results from these efforts are summarized.

A small fraction of double radio sources show a peculiar and striking hybrid morphology; they have a distinctly FR I structure on one side of the nucleus, and a FR II structure on the other. We argue that the mere existence of these HYMORS is quite incompatible with the theoretical explanations for the Fanaroff-Riley dichotomy that are based upon the nature of the jet plasma, or those invoking an intrinsic property of the central engine. Rather, these HYMORS strongly support models that explain the difference between FR I and FR II sources in terms of asymmetry of interaction of the jets with the external environments. We further show that a model for radio source dynamics we had earlier proposed can neatly reproduce the observed dependence of the radio power dividing the two FR classes on the optical luminosity of the host galaxy, as found by Owen & White and Ledlow & Owen

We provide observational and theoretical perspectives on the currently much debated issue of the Fanaroff-Riley (FR) morphological dichotomy of extragalactic radio sources. In this context we introduce a new, albeit rare, class of double radio sources in which the two lobes exhibit clearly different FR morphologies. It is argued that such `HYbrid MOrphology Radio Sources', or HYMORS, could be used to effectively constrain the theoretical mechanisms proposed for the FR dichotomy. Basically, the existence of HYMORS supports explanations for the FR dichotomy based upon jet interaction with the medium external to the central engine, and appears quite difficult to reconcile with the class of explanations that posit fundamental differences in the central engine, such as black hole spin or jet composition, to be responsible for the two FR classes of double radio sources.

The possibility that some of the observed X-ray and optical variability in active galactic nuclei and galactic black hole candidates is produced in accretion disks through the development of a self-organized critical state is reconsidered. New simulations, including more complete calculations of relativistic effects, do show that this model can produce light-curves and power-spectra for the variability which agree with the range observed in optical and X-ray studies of AGN and X-ray binaries. However, the universality of complete self-organized criticality is not quite achieved. This is mainly because the character of the variations depend quite substantially on the extent of the unstable disk region. If it extends close to the innermost stable orbit then a physical scale is introduced and the scale-free character of self-organized criticality is vitiated. Significant dependence of the power spectrum density slope on the type of diffusion within the disk and a weaker dependence on the amount of differential rotation are noted. When general relativistic effects are incorporated in the models, additional substantial differences are produced if the disk is viewed from directions far from the accretion disk axis.

We show that a certain class of flare models for variability from accretion disk coronae are subject to beam-plasma instabilities. These instabilities can prevent significant direct acceleration and greatly reduce the variable X-ray emission argued to arise via inverse Compton scattering involving relativistic electrons in beams and soft photons from the disk.

A key argument in favor of orientation based unification schemes is the finding that among the most powerful 3CRR radio sources the (apparent) median linear size of quasars is smaller than that of radio galaxies, which supports the idea that quasars are a subset of radio galaxies, distinguished by being viewed at smaller angles to the line of sight. Recent measurements of radio sizes for a few other low frequency samples are, however, not in accord with this trend, leading to the claim that orientation may not be the main difference between radio galaxies and quasars. We point out that this ``inconsistency'' can be removed by making allowance for the temporal evolution of sources in both size and luminosity, as inferred from independent observations. This approach can also readily explain the other claimed ``major discrepancy'' with the unified scheme, namely, the difference between the radio luminosity--size correlations for quasars and radio galaxies.