Infrared Variation of Radio Selected BL Lacertae Objects

Abstract: In this paper, the historical infrared (JHK) data compiled from the published literature are presented in electronic form for 40 radio selected BL Lacertae objects (RBLs) for the first time. Largest variations are found and compared with the largest optical variation. Relations between color index and magnitude and between color-color indices are discussed individually. For the color-magnitude relation, some objects (0048-097, 0735+178, 0851+202, 1215+303, 1219+285, 1749+096, and perhaps 0219+428, 0537-441, 1514-241) show that color index increases with magnitude indicating that the spectrum flattens when the source brightens while some other objects (0754+100, 1147+245, 1418+546, and 1727+502) perhaps show an opposite behaviour, while remaining objects do not show any clear tendency; For color-color relation, it is found common that (J-K) is closely correlated with (J-H) while (J-H) is not correlated with (H-K).

\section{Introduction}
 While the nature of active galactic nuclei (AGNs) is still an open problem,
 the study of AGNs variability  can yield valuable information about
 their nature.  Photometric observations of AGNs are important to
 construct light curves and to study  variability behavior over
 different time scales. In AGNs, some long-term optical variations have been
 observed and in some cases claimed to be periodic (see Fan et al. 1998a).

 Before $\gamma$-ray observations were available, Impey \& Neugebauer (1988)
 found that the infrared emission  (1-100 $\mu m$) dominates the bolometric
 luminosities of blazars. The infrared emissions are also an  important
 component for the luminosity even when the $\gamma$-ray emissions are
 included ( von Montigny 1995). Study of the infrared will provide much
 information of the emission mechanism. The long-term optical variations
 have been discussed in the papers: Webb et al. (1988); Pica et al. (1988);
 Bozyan et al. (1990); Stickel et al. (1993); Takalo (1994); Heidt \& Wagner
 (1996) and references therein. The long-term infrared variations have been
 presented for some selected blazars in a paper by Litchfield et al. (1994).

 BL Lacertae objects are an extreme subclass of AGNs showing rapid and large
 amplitude variability, high and variable polarization, no or weak emission
  line features (EW $ < 5\AA)$. From surveys, BL Lacertae objects can be
 divided into radio selected BL Lac objects (RBLs) and X-ray selected BL
 Lac objects (XBLs); From the differences in their overall spectral energy
 distributions as represented by their location in the  $\alpha_{ro}$ versus
 $\alpha_{ox}$, they can be divided into quasar-like BL Lacertae objects
 (Q-BLs) and X-ray strong BL Lacertae objects (X-BLs) (Giommi et al. 1990);
 From the difference in the peak frequency of the spectral energy
 distribution (SED), they can be divided into high frequency peaked BL Lac
 objects (HBLs) and low frequency peaked BL Lac objects (LBLs) (Padovani
 \& Giommi 1995; Urry \& Padovani 1995).  In the last catalogue of active
 galaxies, Veron \& Veron (1998) list 357 BL Lac objects and BL Lac object
 candidates.

 There are some obvious differences between  RBLs and  XBLs (see Fan 1998
 in preparation, for a brief review). Some authors claim that they are the
 same type of objects with XBLs
 having wider viewing angle than RBLs (Ghisellini \& Maraschi, 1989;
 Georganoppoulos \& Marscher 1996; Fan et al. 1997); Other claim
 that there is a continuous spectral sequence from XBLs to RBLs (Sambruna
 et al. 1996) while some claim that they are different classes
 (Lin et al. 1997).

 Infrared observations have been done for BL Lac objects for more than 20
 years. Observations in the infrared were obtained originally on BL Lacertae
 by Oke et al. (1969). Strittmatter et al. (1972) observed the continuity
 of spectral-flux distribution from the optical to infrared wavelengths.
 Afterwards, some other authors also did similar observations (see Stein
 1978 for summary).  Using the 1.26-m infrared telescope at Xing Long
 station, Beijing observatory, we obtained some infrared observations of a
 few BL Lac objects ( Xie et al. 1994; Fan et al. 1998b). But long-term
 infrared variations are not available in the literature.
 In this paper, we mainly present the long-term infrared (J,H, and K bands )
 light curves for some RBLs and discuss their variability properties.
  The paper has been arranged as follows: In section 2, we review the
 literature data and  light curves; in section 3, we discuss the
 variations and give some remarks for each object,  and a short  summary.

\section{Near-infrared Light Curves}

\subsection{Data from Literature }

 Infrared observations are available since the end of the 1960s. We
 compiled the data  from 58 publications. They are listed in table 1 \& 2,
  which gives the observers in Col. 1.  and the telescope(s) used in Col.2.
 For the J, H, and K magnitudes, they are listed in an table 4 in
 electronic form. In table 4, Col. 1 gives the name; Col. 2. the JD time;
 Col. 3, the J magnitude; Col. 4 the uncertainty for J; Col. 5, the
 H magnitude; Col. 6, the uncertainty for H; Col. 7, K magnitude;
 Col. 8, the uncertainty for K.

\subsection{Data Analysis}

 The flux densities from the literature had been converted back to magnitudes
 using the original formulae.  In the literature, different telescopes
 are used, different telescopes use photometers with slightly
 different filter profiles, resulting in slightly  different calibration
 standards and zero-points, but the  uncertainty aroused by the different
 systems is no more than a few percent.

 The magnitude is dereddened using
 $A(\lambda)=A_{V}(0.11\lambda^{-1}+0.65\lambda^{-3}-0.35\lambda^{-4})$
 for $\lambda > 1 \mu m$ and  $A_V$ = 0.165(1.192- $\tan{b})\csc{b}$ for
 $ |b| \leq 50^{\circ}$ and $A_V$ = 0.0 for  $ |b| > 50^{\circ}$
 (Cruz-Gonzalez \& Huchra 1984, hereafter CH (84); Sandage 1972;
 Fan et al. 1998c).  It is clear that some RBLs have many observations
 while others have only a few data points, in Figures 2-22, we  show only
 those with  more observations.

 For data shown in Figs. d - i of each object we have
 performed linear fit with uncertainties  in both coordinates
 considered (see Press et al. 1992 for detail),
 \begin{equation}
 y(x) = a + bx
 \end{equation}
 In principle, $a$ and $b$ can be determined  by
 minimizing the $\chi^{2}$ merit function (i.e. equation 15.3.2 in
 Press et al. book)
 \begin{equation}
\chi^{2}(a,b)= {\sum\limits_{i=1}^{n} {(y_{i}-a-bx_{i})^{2}w_{i} }}
 \end{equation}
 where ${\frac{1}{w_{i}}}=\sigma^{2}_{y_{i}}+b^{2}\sigma^{2}_{x_{i}}$,  $\sigma_{y_{i}}$ and $\sigma_{x_{i}}$ are the $x$ and $y$ standard
 deviations for the $i$th point
 Unfortunately, the occurrence of $b$ in the denominator of the above
 $\chi^{2}$ merit function, resulting in equation
 ${\frac{\partial \chi^{2}}{\partial b}}=0$ being
 nonlinear makes the  fit very complex, although we can get a
 formula for $a$ from
  ${\frac{\partial \chi^{2}}{\partial a}}=0$:
 \begin{equation}
  a={\frac{\sum\limits_{i}w_{i}(y_{i}-bx_{i})}{\sum\limits_{i}w_{i}}}
 \end{equation}
  Minimizing the $\chi^{2}$ merit function, equation (2), with respect to
 $b$ and using the equation (3) for $a$ at each stage to ensure that the
 minimum with respect to $b$ is also minimized with respect to $a$, we
 can obtain $a$ and $b$. As Press et al. stated, the linear correlation
 coefficient $r$ then can be  obtained by  weighting the relation (14.5.1)
  terms of Press et al. (1992).
 Finding the  uncertainties, $\sigma_{a}$ and $\sigma_{b}$, in $a$ and $b$ is
 more complicated (see Press et al. 1992). Here, we have  not performed this.
 For the data whose uncertainties were not given in the original literature
 are not included in our linear fit or in Fig. d-i, but they
 are included in the light curves (Fig. a - c).

 As discussed in the paper of Massaro \& Trevese (1996), there is a
 statistical bias in the  spectral index-flux density correlation. Following
 their suggestion, we considered the relation between magnitude in one
 band and the color index obtained fro two other  bands to avoid this
 bias. The relations are shown in Fig. d-f.

\section{Discussion}

\subsection{Variations}

 BL Lac objects are variable at all wavelengths. For the RBLs with more
 infrared data, we presented their variability properties in Table 3, which
 gives the name in Col. 1, redshift in Col. 2, Reddening correction, $A_{V}$,
 in Col. 3, largest optical polarization ($P_{opt}$) in Col. 4, largest
 amplitude optical variation ($m_{opt}$) in Col. 5, largest amplitude
 variations in J, H, and K in Col. 6, 7, and 8, the averaged color index
 (J-H) and (H-K) in Col. 9, and 10, the uncertainty in Col. 9  and 10
 is the one sigma
 deviation. The infrared variations are very large for some objects.
 The color indices are also variable and associated with the magnitude
 for some objects,  indicating that the spectrum flattens when the source
 brightens while some objects show different behaviours. For
 variation, it is reasonable for larger variation to correspond to shorter
 wavelength, but this paper shows that some objects do show a different
 behaviour.
 The objects showing large optical variation also show large infrared
 variations and high polarizations.  Large variation and high polarization
 are associated. For the color-indices, it is common that there is a
 correlation for (J-K) and (J-H) while there is almost no correlation for
 (J-H) and (H-K) although some objects show different properties.
 But, when the averaged color indices are taken into account, there are
 mutual correlations between any two color indices of (J-H), (J-K) and
 (H-K) with the correlation between (J-H) and (H-K) being not so
 close as compared with other two correlations. (J-K) = 2.26(J-H) - 0.28
 with correlation coefficient  $r$ = 0.90 and the probability of chance
 correlation  $p = 6.3\times 10^{-11}$, (H-K) = 1.44(J-H) -
 0.41 with $r$ = 0.80 and $p = 4.4\times 10^{-7}$, and (J-K) = 1.62(H-K) +
 0.33  with $r$ = 0.93 and $p = 6.1\times 10^{-13}$. The (J-H) and (H-K)
 relation is shown in Fig. 1.  In the next section, individual objects
 are discussed.

\subsection{Remarks}
\subsubsection{PKS 0048-097, $OB-081, MC3, PHL856$}

 There is no redshift available for this polarized ( $P_{3.7cm}=9.3\%$,
 Wardle 1978;  $P_{opt}=$1$\%$-27 $\%$, CH(84);  Brindle et al. 1986 )
 object. Very faint infrared  fluxes have been observed by Gear (1993),
 which leads
 the object to have the largest infrared variation among the objects
 discussed in the present paper. The infrared variations  are $\Delta J$ =
 6.55,  $\Delta H$ = 6.23, and $\Delta K$ = 6.07 while that  in the
 optical band is 2.7 mag (Angel \& Stockman 1980). We can expect that the
 optical  monitoring should get a similar amplitude.  For
 color-color indices,  there are correlations for (J-K) and (J-H) and for
 (J-K) and (H-K):
 (J-K) = 1.20(J-H) + 0.61 with $r$ =0.69 and  $p = 5.\times 10^{-4} $,
 (J-K) = 4.7(H-K) - 2.13  with $r = 0.64$ and $p = 2.\times 10^{-3}$,
 but no  correlation was found for (J-H) and (H-K). For the color indices
 and the  magnitudes,  there are correlation indications of (H-K) increasing
 with J ($p \sim 1\%$) and  (J-K) with H ($p \sim 0.3\%$) (see Fig 2d, 2e).

\subsubsection{PKS 0118-272}

 An absorption feature at $z$ = 0.559 was found by Falomo (1991). It was
 observed  in the infrared during the period JD 244000+ 4773-7753. For
 this object,  there are no correlations  for color-color indices or for
 color index and  magnitude.  It is a highly polarized object
 ( $P_{opt} \sim 17 \%$) and  shows $ dP/d\nu > $ 0 (Mead et al. 1990).

\subsubsection{PKS 0138-097}

 It was observed by Falomo et al (1993) at the visual and 2.2$\mu m$
 wavelengths,  the spectral flux distribution is well described by a single
  power law of
 $\alpha$ = 1.0 ($f_{\nu} \sim \nu^{-\alpha}$). It shows variable and high
 optical polarization ($P = 6\%-29.25\%$)(Mead et al. 1990). The variation in
 the infrared is almost the same as that in the  optical band ($\Delta m$ =
 1.52 on a time scale of one week, Bozyan et al. 1990).
  There is a correlation for (J-K) and (J-H) but not for other colors.

\subsubsection{0219+428, $3C 66A$}

 A redshift of 0.444 has been derived from a possible  MgII  $\lambda$2800
 feature (Miller et al. 1978). A rapid decline of 1.2 mag in 6 days was
 observed by Takalo et al. (1992). Different properties for the color and
 brightness are found in this source: Takalo et al. (1992) found that
 although the spectral index is considerably flat when the source was at
 its observed maximum state, the (J-H) and (H-K) against K color-magnitude
 diagrams are essentially scattered; Massaro et al. (1995) found that there
 is a positive correlation between the infrared spectral index and the
  brightness at more than the 98 $\%$ level while De Diego et al. (1997)
 found a different behaviour.  From the long-term data compiled in the
 present paper no clear correlation has been found for color index
 and magnitude except that there is a tendency for (H-K) to increase with
 J (Fig. 3e). There are close correlations for (J-K) and (J-H) (Fig. 3h)
 and for (J-K) and (H-K) (Fig. 3i) but not for (J-H) and (H-K):
 (J-K) = 1.38(J-H) + 0.33 with $r = 0.76$ and $p = 10^{-7}$,
 (J-K) = 1.60(H-K) + 0.18 with $r = 0.45$ and $p = 5.0 \times 10^{-3}$.
 It shows two lower states in the light curves (Fig. 3 a,b,c). Its
 largest variation of 2.0 mag (Miller \& McGimsey 1978) obtained in the
 optical band is similar to those found in the infrared band (1.67mag in
 the K  for instance). Polarizations of $P_{opt}$ = 6$\%$-18$\%$ and
 $P_{Rad.}$ = 5$\%$ are reported in the papers of (CH(84)) and Takalo (1991).

\subsubsection{AO 0235+164}

 Absorption redshifts of $z_{abs}$ = 0.524 and 0.852 were found by Burbidge
 et al. (1976) and an emission line redshift of $z_{em} = 0.94$ was found by
 Nilsson et al.  (1996). It is one of the most extensively studied objects
 in the infrared  (Impey et al. 1982; Gear et al. 1986a; Takalo et al. 1992),
 variations of  $\Delta J$ = 1.9, $\Delta H$ = 1.8, and $\Delta K$ = 1.9 are
  found by Takalo et al. (1992).  The spectral index is found to be
 correlated with J flux at a significance level of
 more than 90$\%$ by Brown et al. (1989). It is the most highly  polarized
 ($P(V) = 44 \%$,$P_{IR} = 36.4 \%$) and the reddest ($\alpha$ = 4.61) object
 (Impey et al. 1982; Stickel et al. 1993).  The infrared light curves show
 three clear outbursts, and variations  as large as the optical ones
 ( $\Delta m$ = 5.3 mag, Stickel et al. 1993). There are positive correlations
 between color indices against magnitude when it is brighter than J = 14.,
 H = 12.5 and K = 12.0 magnitudes, but the correlation is complex in fainter
 states. There are correlations for (J-K) and (J-H) and
 (H-K): (J-K) = 1.30(J-H) + 0.64 with $r$ = 0.56 and $p = 3 \times 10^{-5}$,
 (J-K) = 1.42(H-K) + 0.56  with $r$ = 0.53 and $p = 9 \times 10^{-5}$.
 (see Fig. 4)

\subsubsection{0300+471, $4C47.08$}

 This object has been observed in the infrared  only once  (Puschell
 \& Stein 1980) with H = 13.44 and K = 12.46, which gives a color index of
 (H-K) = 0.98$\pm$ 0.07.  Its flux distribution from 0.36$\mu m$ to 2.3$\mu m$
 can be expressed by a power law with  $\alpha=1.0$
 ( $f_{\nu} \propto \nu^{-\alpha}$). It is  11.5$\%$ polarized at 3.7cm
 (Wardle 1978) and 24$\pm$5$\%$ polarized at 0.44$\mu m$ (Puschell \& Stein
 1980).

\subsubsection{0306+103}

 This variable ($\Delta m=2.2$ mag, Leacock et al 1976) object has been
 observed in K for only two times (K =11.82 and 11.92 mag).

\subsubsection{PKS 0537-441}

 Large variation of $\Delta m$ = 5.4 mag (Stickel et al. 1993) in the
 optical and a variation of $\Delta H$ = 2.0 mag (CH(84)) are found in
 this highly polarized ($P_{opt} = 18.7\%$, Impey \& Tapia 1988) object.
 From our compiled data, a large variation of 3 mag in J and weak
 correlations  between  (J-K) and  (J-H) and (H-K) are found:
 (J-K) = 1.92(J-H) + 0.04 with $r$ = 0.54 and $p = 2.7\%$,
 (J-K) = 1.42(H-K) + 0.47 with $r$ = 0.547 and $p = 2.7\%$.
 Figs. 5d,e show that there is a tendency for the color index to increase
 with the magnitude.

\subsubsection{PKS 0735+178, $VRO 17.07.02, OI158$}

 It has been observed in the optical band for about 90 years (Fan et al.
 1997) and in the infrared for about 20 years (Lin \& Fan 1998). The largest
 optical variation of $\Delta m = 4.6$ (Fan et al. 1997) is greater than
 that in the infrared (Lin \& Fan 1998). Rapid variability of 0.4 mag in
 about 15 minutes and a  positive correlation between infrared spectral
 index and brightness are found (Massaro et al. 1995; Brown et al. 1989).
 It shows correlations between H and (J-K) and between K and (J-H) for
 long-term data: (J-K) = 0.12H + 0.17 with $r$ = 0.447 and
 $p \sim 5.3 \times 10^{-5}$);
 (J-H) = 0.04K + 0.32 $r$ = 0.30 and $p \sim 7.8 \times 10^{-3}$ )
 (Fig. 6d and 6e, see also Lin \& Fan 1998). (J-K) is found to be correlated
 with both (J-H) and (H-K) (see Lin \& Fan 1998).
 It is highly polarized in the radio ($P_{3.7cm} = 4.7\%$, Wardle 1978),
 in the infrared ($P_{IR} = 32.6\%$, Impey et al. 1984) and in the optical
 ($P_{opt}=35\%$, see Lin \& Fan 1998) bands.

\subsubsection{PKS 0754+100, OI 090.4}

 PKS 0754+100 is  observed to show rapid variability in the infrared
 ($P_{IR} = 4\%-19\%$)  and in the visual  ($P_{Opt.} = 4-26 \% $)
 polarizations, a variation of  $6\%$ over 4 hours in optical polarization
 was observed by Impey et al. (1982, 1984).  The polarization is independent
 of wavelength between the visual and the  2.2$\mu m$ regions
 (Craine et al. 1978). An optical variation  of 2.0 magnitude
 was found from the examination of the plate archives by Baumert (1980)
 and a variation of 1 mag in the infrared was reported (CH84). It is the only
 object in Massaro et al (1995) sample which, despite a recorded
 variability of  about 1 mag, does not show statistically significant
 changes of the spectral index.
 0754+100 shows a decreasing brightness since the first infrared observation
 although there are some brightness fluctuations in the light curves.  It is
 interesting that the variation in H is greater than those in J and K with
  variability ranges of J = 12.31 to 14.08, H = 11.21 to 13.27 and
 K = 10.55 to 12.43. Besides, although an anti-correlation between H and (J-K)
 appears in Fig. 7d, it is far from any firm conclusion. (J-K) is
 correlated  with both (J-H) and (H-K):
 (J-K) = 1.3(J-H) + 0.54 with $r$ = 0.78 and $p = 2.4 \times 10^{-8}$,
 (J-K) = 2.16(H-K) - 0.15 with $r$ = 0.61 and $p = 8.9 \times 10^{-5}$.

\subsubsection{S4 0814+42, $OJ425$}

 This object is observed to show variations of 1.62, 1.63, and 1.52 mags
 in J, H, and K respectively (Gear 1993). Mean color indices of (J-H)
  = 0.96$\pm$0.05, (J-K) = 1.90$\pm$0.04, and (H-K) = 0.95$\pm$0.07 can
 be obtained from those observations. Polarizations of $P_{opt} = 12 \%$
 (Impey et al. 1991) and $P_{3.7cm}=9.8$ (Wardle 1978) are reported.

\subsubsection{PKS 0818-128, OJ-131}

 It shows optical variation from 14.6 to 17.5 magnitude. The
 linear polarization is high and variable with $P_{opt} = 8\%$-36$\%$
 (Craine et al. 1978) and $P_{IR} = 9\%-14\%$ (Impey et al. 1984).
 The infrared variation of 2.0 mag is a little smaller than  in the
 optical band ($\Delta m = 2.9$, Craine et al. 1978). There are indications
 of correlations between (J-K) and (J-H) and (H-K), but no definite
 conclusion can be drawn from the few points with large error bars. The
 correlation between the color index and magnitude is also complex as
 shown in Fig 8.

\subsubsection{PKS 0820+22, $4C22.21, MG0823+2222$}

 This RBL was observed twice and does not show any clear variation in the
 infrared: J = 17.63$\pm$0.16-17.76$\pm$0.07, H = 16.54$\pm$0.07, and
 K = 15.59$\pm$0.04-15.75$\pm$0.14 (Gear 1993). Averaged color indices of
 (J-H) = 1.22$\pm$0.10, (J-K) = 2.02$\pm$0.17, and (H-K) = 0.95$\pm$0.08
 can be obtained. It is at a redshift of 0.951 (Stickel et al. 1991) and
 $P_{opt} = 5.2\%$ polarized in optical band (Kuhr \& Schmidt 1990).

\subsubsection{PKS 0823+03, $OJ038, MG0825+0309$}

 This highly polarized ($P_{opt}=22.9\%$, Impey \& Tapia 1988) object shows
 a variation of about 0.8 mag in the infrared: J = 13.73-14.49,
 H = 12.98-13.79, and K = 12.39-13.24 (Gear 1993).  Mean color indices are
 (J-H) = 0.72$\pm$0.04, (J-K) = 1.29$\pm$0.07 and (H-K) = 0.57$\pm$0.03.
 Its redshift is 0.506 (Stickel et al. 1991).

\subsubsection{S4 0828+493,$OJ448, BP077$}

 An infrared variation of about 0.5 mag (J = 15.81-16.39, H = 14.80-15.44,
 K = 14.00-14.53) is observed (Gear 1993) from this polarized ($P_{opt}
 = 7.9\%$,  Kuhr \& Schmidt 1990) object.  Mean color indices are
  (J-H) = 0.98$\pm$0.04,  (J-K) = 1.79$\pm$0.04, and (H-K) = 0.81$\pm$0.14.
 The recorded  optical variation is $\Delta R$ = 2.0 $mag$
 (Bozyan et al. 1990).

\subsubsection{PKS 0829+046, $OJ 049$}

 It is also an extremely variable object ( $\Delta m$ = 5.0 Angel \&
 Stockman 1980).  Polarizations of $P_{3.7cm}=4.8\%$ ( Wardle 1978) and
 $P_{opt} = 12\%$ (CH(84))  are reported. It is hard to draw any conclusion
 for the relation between color  index and  magnitude.  It shows a
  correlation between (J-K) and (J-H):
 (J-K) = 1.32(J-H) + 0.48 with $r$ = 0.76 and $p = 4.3 \times 10^{-5}$.
 Between  (J-K) and (H-K), a close correlation with $r$ = 0.68 and
 $p = 5.4\times 10^{-4}$ can be obtained once the outlying point is excluded
 (see Fig. 9).

\subsubsection{PKS 0851+202, $VRO 20.08.01, OJ287$}

 OJ287 is one of the very few AGNs for which continuous light curve over more
 than one hundred years has been observed.  A redshift of $z$ = 0.306 is
 derived from  [OIII]$\lambda$5007 by Miller et al. (1978)
 and confirmed by Sitko \& Junkkarinen (1985). Its observed properties from
 the radio to the X-ray have been reviewed by Takalo (1994). OJ287 shows
 a 12-year  period in the optical (Sillanpaa et al. 1988;
 Kidger et al. 1992)  and in the infrared (Fan et al. 1998b). It is highly
 polarized in the radio  ($P_{3.7cm} = 16.6\%$, Wardle 1978), infrared
 ($P_{IR} = 14.3\%$, Impey et al.  1982) and optical ($P_{opt} = 37.2\%$,
 Smith et al. 1987) bands.  Its optical variation of 6.0 $mag$ (Fan et al.
 1998b) is larger than the  infrared variation. There is a tendency for the
 color index to increase  with magnitude. (J-K) is correlated with
 both (J-H) and (H-K):
 (J-K) = 1.62(J-H) + 0.35  with $r$ = 0.76,
 (J-K) = 1.54(H-K) + 0.35  with $r_{s}$ = 0.65 (see Fig. 10 and Fan et al.
 1998b).

\subsubsection{B2 1101+384, $Mkn 421,  OM303$ }

 This object is found to have nebulosity with a redshift of 0.0308 by
 Ulrich et al. (1975). A large optical variation of 4.6 mag was observed
 (Stein  et al. 1976). The observed optical, infrared
 and radio polarizations are, respectively,  0$\%$-7.9$\%$ (Rieke et al.
 1977), 2$\%$-4$\%$ (CH(84)) and $3.7\%$ (Wardle 1978).
 (J-K) is correlated with (J-H):
 (J-K) = 0.72(J-H) + 0.69 with $r$ = 0.68 and $p = 6.1 \times 10^{-4}$.
 (H-K) shows an anti-correlation with (J-H), (H-K) = -0.68(J-H) + 0.94
 ($p = 6.0 \times 10^{-3}$).  For (J-K) and (H-K), it is hard to draw a
 definite conclusion as  shown in Fig.11i.  No definite color index-magnitude
 relation can be obtained. The infrared variation  found here is much
 smaller than the optical band one.

\subsubsection{PKS 1144-379}

 It is a high-redshift object with $z$ = 1.048 ( Stickel et al. 1989). An
 optical  variation of 1.92 mag  (Bozyan et al. 1990) and an optical
 polarization of $P_{opt} = 8.5\%$  (Impey \& Tapia 1988) are reported.
 The infrared variation is larger at lower frequency
 and is greater than in the optical.

\subsubsection{B2 1147+245, $OM280, MG1150+2417$}

 No redshift is published  for this source. Similar infrared ($P_{IR} =
 14.8\%$)  and optical ($P_{opt} = 1.5\%-13\%$) polarizations (Impey et al.
 1982), and  an  optical variation of 0.77mag (Pica et al. 1988) are reported.
 The variation in H is found greater than those in J and K. Fig. 12f shows
 that the  spectrum flattens when the source dims, a similar result was
 also found by Lorenzetti et al. (1989). (J-K) is associated
 with (J-H):
 (J-K) = 1.2 (J-H) + 0.67 with $r$ = 0.90 and $p = 1 \times 10^{-6}$
 (see Fig. 12).

\subsubsection{B2 1215+303, $ ON325$}

 For this polarized ($P_{3.7cm} = 6.3\%$ and $P_{opt} = 14\%$, Wardle 1978)
 object, we found (J-H)  to be correlated with  K:
 (J-H) = 0.24K - 2.13  with $r$ = 0.59 and $p = 7.2 \times 10^{-3}$. The
 infrared variation found here is smaller than that in the optical band
 ($\Delta m$ = 3.3 mag, Zekl et al. 1981). (J-K) is correlated with (J-H) and
 (H-K): (J-K) = 1.06(J-H) + 0.60 with $r$ = 0.69 and $p = 1.2\times 10^{-3}$,
 (J-K) = 1.68(H-K) + 0.27 with $r$ = 0.42  and $p = 5.3\%$. (see Fig. 13)

\subsubsection{B2 1219+285, $W Com, ON231, VR 28.12.02, GC$}

 This extremely variable object ($\Delta V$ = 5 mag) was observed to be
 polarized at $P_{3.7cm} = 5.8\%$ (Wardle 1978) and  $P_{opt} = 2\%-10\%$
 ( CH(84)). A correlation between H and (J-K) was found: (J-K) = 0.28H - 2.01
 with $r$ = 0.62 and $p \sim 1\%$. A correlation is also found between (J-K)
  and (J-H): (J-K) = 1.04(J-H) + 0.72 with $r$ = 0.61 and $p = 1.1\%$
 (see Fig. 14h).
 For (J-K) and (H-K), there is a close relation if the two outlying points
 are excluded.  In table 3, we can see that the variation in J is smaller
 than those in H and K. The reason is that there are only  H and K data
 for the faint state (Glass 1981).  The variations in J and H  are the
 same, and greater than in K if we do not consider the low state points.

\subsubsection{B2 1308+326, $ OP 313, GC$}

 An emission line redshift of 0.996 comes from the identification of a fairly
 strong, broad emission feature at 5586$\AA$ as MgII$\lambda$2800, while the
 absorption redshift of 0.879 is from a weak doublet at $\lambda$5252.4$\AA$
 and
 5264.4 $\AA$ (Miller et al. 1978). B2 1308+326 shows high and variable
 polarization  indicating variability of $\Delta P = 10\%$ on a time scale
 of 24 hours (Impey et al. 1984) and a decreasing of 12$\%$ in two days
 (Moore et al. 1980).  The recorded maximum polarizations of
 $P_{IR} = 19.6\%$ and $P_{opt} = (28\pm 5)\%$  at 0.44$\mu m$ was also
 observed by Puschell et al. (1979). An optical variation
 of $\Delta m$ = 5.6 mag has been observed (Angel \& Stockman, 1980). When the
 infrared emission increases the spectral index does not change (Sitko et al.
 1983),  but an indication of a positive correlation between the infrared
 spectral index  and brightness is found by Brown et al. (1989). The
 1977/1978 optical outbursts  show a double peak structure separated by
 1 year. Following the 1978 peak,  Puschell et al (1979) found that the
 infrared fluxes did not decrease with time  as rapidly as the visual fluxes
 and proposed that the more rapid decline  of visual compared with infrared
 light is a general characteristic  of these events.  From the light curve,
 we can see that the brightness decreases with time with some flux
 fluctuations showing up. No definite relation can be obtained for the
 spectral index and magnitude  for the long-term data, but there is a
 correlation between (J-K) and (J-H):
 (J-K) = 1.50(J-H) + 0.48 with $r$ = 0.596 and $p = 1.8 \times 10^{-3}$
 (see Fig. 15).

\subsubsection{MC 1400+162, $4C 16.39, OO100$}

 This variable object ($\Delta m$ = 2.8, Zekl et al. 1981 ) is found in
 a cluster  of galaxies (Miller 1978b) with a redshift of
 $z_{em}$ $\sim$ 0.245  from the  [OII], $H_{\beta}$, [OIII],
 [OII] (Miller et al. 1978) and show a  variable polarization
 ($P_{opt}=4.4\%-14\%$, Angel \& Stockman 1980; Impey et al. 1984). The
 infrared observations show a 0.54 mag variation  in K and
  (H-K) = 0.89$\pm$0.11.

\subsubsection{B2 1418+546, $OQ530$}

 1418+546 shows a large optical variation of 4.8 magnitude ranging from
 11.3 to  16.1 magnitude in Miller's (1978) study; its  optical polarization
 from 5$\%$ to 24$\%$ (Marchenko 1980) is greater than  the infrared
 polarization of $P_{IR}=19\%$ (Impey et al. 1984). In the infrared, some
 properties have been shown in the papers of O'Dell et al. (1978) and
 Puschell \& Stein (1980). No correlation was found between  spectral
 index and brightness, but the flux variations are found larger
 at lower frequency by Massaro et al. (1995) and Takalo et al (1992).
 For the spectral  shape and the brightness of the source,  the long-term
 data show  (see Fig. 16):  (J-K) = -0.08H + 2.65 with
  $r$ = -0.40 and $p = 2.6\%$, and (J-H) = -0.08K + 1.75 with $r$ = -0.41
 and $p = 2.2\%$   indicating that the spectrum steepens when the source
 brightens. No similar correlation has been found in the visual band
 (Worrall et al 1984a).  The infrared light curves show two bursts separated
 by about 8 years and the  variation at lower frequency to be greater.
 The infrared variations are smaller than  in the optical band.
 The (J-K) is found related with (J-H) and (H-K) (see Fig. 16):
 (J-K) = 1.56(J-H) + 0.43 with $r$ = 0.856 and $p = 1.6\times 10^{-8}$,
 (J-K) = 1.98(H-K) + 0.03 with $r$ = 0.61  and $p = 7.0\times 10^{-4}$.

\subsubsection{PKS 1514-241, $AP ~Librae$}

 This variable ( $\Delta m = 3.0$mag, Branly et al. 1996) object is located
 at the center of an elliptical galaxy with an emission line redshift of
 $z = 0.049$  (Disney et al. 1974), polarizations of $P_{3.7cm} = 4.4\%$
 (Wardle 1978) and  $P_{IR} = 2\%-7.4\%$ and $P_{opt} = 7\%$ (Impey et al.
 1982) are reported.  The infrared data show  a correlation between (J-K)
 and (H-K):  (J-K) = 0.18(H-K) + 0.84 with $r$ = 0.928 and
 $p = 4.5\times 10^{-5}$.  For color  index and magnitude, there is an
 indication of (J-K) increasing with H but this is far from any certain
 conclusion (see Fig. 17).

\subsubsection{PKS 1519-273}

 This polarized (  $P_{opt} = 11.4\%$, Impey \& Tapia 1988) object has
 only been
 observed  twice by Gear (1993) and Glass (1981).  The observations
 indicate color indices of  (J-H) = 0.45$\pm$0.21, (J-K) = 1.35$\pm$0.17,
 and (H-K) = 0.90$\pm$0.24.  The optical variation is $\Delta m$ = 2.43
 (Fan 1997).

\subsubsection{PKS 1538+149,$4C14.60, OR165$}

 Polarizations of $P_{3.7cm} = 9.3\%$ (Wardle 1978) and $P_{opt} = 32.8\%$
 (Brindle et al 1986) are reported.  The infrared variations found here from
 few data points  is much smaller than in the optical bands
 ( $\Delta m$ = 3.7, Stickel et al. 1993).

\subsubsection{B2 1652+398, $Mkn 501, 4C39.49, OS387$}

 B2 1652+398 was shown to have a redshift of 0.0337 by Ulrich et al. (1975).
 It has effective spectral indices similar to that of  X-ray selected
 BL Lac objects but has been classified as an RBLs. The largest optical
 variation of  1.3 mag (Stickel et al. 1993) is smaller than in J.
 Its optical and  infrared emissions are polarized at $P_{opt} = 2\%-7\%$
 and $P_{IR} = 3 \%$ ( CH(84)).  During the flare the colors are found to be
 very
 much bluer than in quiescent state,  reaching (H-K) = 0.13 and (J-H)=0.48
 (Takalo et al. 1992), but there are almost  no correlation between
  color index  and magnitude  or for color-color  indices (see Fig. 18).

\subsubsection{1727+502, $I ZW 187$}

 IZW 187 was observed spectroscopically by Oke (1978) to show a redshift of
 0.0554. It is interesting that the optical polarization ($P_{opt} = 4\%-6\%$)
 is smaller than the radio polarization ($P_{3.7cm} = 11.3\%$) (Kinman, 1976).
 IZW 187 shows the highest ratio of X-ray luminosity to optical luminosity,
 $L_{X-ray}/L_{opt.}$,  among BL Lac objects observed by $Einstein$
 observatory (Ku et al. 1982). An optical  variation of 2.1 mag was reported
 by Scott (1976) (see also Bozyan et al. 1990).  Recently, it was
 found to have a massive black hole ($5.4\times 10^{8}M_{\odot}$) at its
 center  in our previous paper (Fan 1995).  Variations of $\Delta H $ =
 0.57mag,  and $\Delta K $ = 0.82 mag are found.
 There is an indication of color index (J-K) decreasing  with H. No
 correlation for the color-color indices is really found (see Fig. 19).

\subsubsection{PKS 1749+096, $OT 081, 4C 09.57$}

 This highly polarized ($P_{3.7cm} = 7.2\%$, Wardle 1978; $P_{opt} = 32 \%$,
 Brindle et al. 1986 ) BL Lac object  was found to show an optical variation
 of 2.7 mag  (Stickel et al. 1993). The infrared data show an indication of
 a correlation  between J and (H-K) with  $r_{S}$ = 0.83 and
 $p = 1.1\times 10^{-3}$. (J-K) is found correlated with (J-H) and (H-K):
 (J-K) = 1.44(J-H) + 0.48 with
 $r$ = 0.79 and $p = 2.4\times 10^{-3}$, (J-K) = 3.31 (H-K) -1.13 with
 $r$ = 0.75 and $p = 5\times 10^{-3}$ (see Fig. 20).

\subsubsection{S5 1803+78}

 There are only observations of  $f_{J } = 5.1 mJy$, $f_{H} = 5.5  mJy$,
 and $f_{K} = 6.7 mJy$,  which gives J = 13.68, H = 13.13, and K = 12.42 mag
 in the paper of Heckman et al. (1983) for this highly polarized
 ($P_{opt} = 35.2\%$, Kuhr \& Schmidt, 1990) object.

\subsubsection{S4 1807+698, $3C371.0, 4C 69.24$}

 The N galaxy 3C371.0 ($z_{em}$ = 0.0506 Miller, 1975) is related to BL
 Lac object  because of its optical variability ($\Delta m \sim $ 2.0 mag)
  and polarization  ($P \sim 12 \%$) (Shaffer 1978; Kinman 1976; CH(84)).
 Miller (1975) suggests  that the polarization increases significantly in
 the UV band if the visual polarization were diluted by
 starlight. Capps \& Knacke (1978)  reported an expected lower-than-optical
 infrared polarization ($P_{IR} = 5.3\pm1.5\%$)  when the starlight diluted
 heavily the infrared.  There are only 7-night of infrared data in the
 literature, giving  $f_{J} = 13.8  mJy$,
 $f_{H} = 16.0 mJy$, and $f_{K} = 19.6 mJy$ (Heckman et al. 1983), included,
 which show a variation of 1 mag in K (from 10.86 to 11.88 mag).

\subsubsection{S4 1823+568, $4C56.27$}

 There are only 2-night of infrared data for this polarized
 ($P_{opt} = 16.8\%$,
 Impey et al. 1991) object: J = 14.87-14.92, H = 13.92, and K = 13.02-13.08
 (Gear 1993), averaged color indices of (J-H) = 0.97$\pm$0.03,
 (J-K) = 1.85$\pm$0.07, and (H-K) = 0.87$\pm$0.04 can be obtained.

\subsubsection{PKS 2005-489}

 This is a moderately variable object ($\Delta m$ = 0.53, Bozyan et al. 1990).
 The limited infrared data only show a variation of about 0.25 mag. No
 polarization is published.

\subsubsection{PKS 2131-021, $4C 02.81$}

 There are only two night of infrared data showing a variation of 0.8 mag
 (J = 16.06-16.88, H = 15.05-15.97, K = 14.23-14.88),  which is comparable
 with the optical variation of $\Delta m$ = 1.2 mag (Stickel et al. 1993).
 Mean infrared color indices of (J-H) = 0.96$\pm$0.07,
 (J-K) = 1.91$\pm$0.12,  and (H-K) = 0.96$\pm$0.16 can be obtained.
 $P_{opt} = 16.9\%$ is reported by  Kuhr \& Schmidt (1990).

\subsubsection{PKS 2200+420, $VRO 42.22.01, BL Lacertae$}

 The prototypical object of its class lies in a giant elliptical galaxy at a
 redshift of $\sim$ 0.07. Maximum radio polarization of $P_{3.7cm} = 6.4\%$
 (Wardle 1978), infrared polarization of $P_{IR} = 15.1 \%$
 (Impey et al. 1984)  and optical polarization of $P_{opt} = 23\%$
 (Kuhr \& Schmidt 1990) are obtained.  It has been observed for about 27
 years in the infrared (see Fan et al. 1998c) and about 100 years in the
 optical (Fan et al. 1998a). A 14-year period  and a maximum optical variation
 of $\Delta B = 5.31$ are found from the B light curve (Fan et al. 1998a).
 The optical variation  is greater  than infrared variation.  Correlations
 have been found between (J-K) and (J-H)  and (H-K) but not between
 color index and magnitude (Fig 21, also see Fan et al. 1998c).

\subsubsection{PKS 2240-260}

 There are only five-night of infrared data, showing a variation of
 $\Delta K$ = 1.7 mag.  There is an indication of a correlation between
 J-H and K,
 but it should be  investigated with more observations. $P_{opt} = 15.1\%$
 is reported for  it (Impey \& Tapia 1988).

\subsubsection{PKS 2254+074, $OY 091$}

 The highly polarized object ($P_{3.7cm} = 8.1\%$, Wardle 1978,
 $P_{opt} = 14\%-21\%,$  $P_{IR}=17.2\%$, Impey et al. 1982) shows an
 optical variation of $\Delta m_{B}$ = 2.37 ( Bozyan et al. 1990).
 The infrared data show correlations between color indices (see Fig.22g-i).

\subsubsection{2335+031,$4C03.59, OZ061$}

 There are only one night of infrared data for this object: J = 16.68,
 H = 15.77, and  K = 14.89, which show color indices of (J-H) = 0.91$\pm$0.06,
 (J-K) = 1.79$\pm$0.06, and (H-K) = 0.88$\pm$0.06. An optical variation of
 0.8 mag and a radio  polarization of $P_{R} = 6\%$ are reported (CH(84)).

\subsection{Summary}

 In this paper, the infrared variations are presented for 40 RBLs, the
 light curves  are shown in Fig. 2 to 22 for 22 objects with more
 observations. The amplitude variation in the optical and infrared
 bands are compared. Correlations between color-color indices and between
 color and magnitude are investigated: (J-K) is always found correlated
 with (J-H). For the color-magnitude relation, some objects ( 0048-097,
 0735+178, 0851+202, 1215+303, 1219+285,  1749+096, and perhaps 0219+428,
 0537-441, 1514-241) show a color index increasing with magnitude, indicating
 that the spectrum flattens when the source brightens while some other objects
 ( 0754+100, 1147+245, 1418+546, and 1727+502) perhaps show an opposite
 behaviour, the rest objects do not show any clear tendency between color
 index and magnitude.

\begin{acknowledgements}{The authors are grateful to the anonymous referee
 for his/her comments and to Dr. G. Adam for his critical reading. J.H. F
 thanks Y. Copin for his help with the figures. This work
 is supported by the National Pan Deng  Project of China}
\end{acknowledgements}
\cite{}


High Energy Emission from AGN and Unified Schemes

Abstract: Active Galactic Nuclei (AGN) are now known to be strong gamma-ray emitters. After briefly describing AGN classification and the main ideas behind unified schemes, I summarize the main properties of blazars (that is BL Lacs and flat-spectrum radio quasars) and their connection with relativistic beaming. Finally, I address the question of why blazars, despite being extreme and very rare objects, are the only AGN detected at high (E > 100 MeV) energies, and touch upon the relevance of TeV astronomy for AGN research.

\section{Active Galactic Nuclei}\label{sec:agn}

Active Galactic Nuclei (AGN) are extragalactic sources, in many cases clearly
associated with nuclei of galaxies (although in the most distant objects the
host galaxy is too faint to be seen), whose emission is dominated by
non-stellar processes in some waveband(s) (typically but not exclusively the
optical). One important feature of AGN is the fact that their emission covers
the whole electromagnetic spectrum, from the radio to the $\gamma$-ray band,
sometimes over almost 20 orders of magnitude in frequency.

It is now well established that AGN are strong $\gamma$-ray ($E > 100$ MeV)
emitters. To be more specific: 1. at least 40\% of all EGRET sources are AGN
(Thompson et al. 1995, 1996; some more AGN are certainly present amongst the
still unidentified sources) and these make up almost 100\% of all
extragalactic sources (the only exceptions being the Large Magellanic Cloud
and possibly Centaurus A); 2. {\it all} detected AGN are blazars, that is BL
Lacertae objects (BL Lacs) or flat-spectrum radio quasars
(FSRQ).\footnote{The term ``blazar'' is here given a wider meaning than the
one sometimes implied, which is restricted to highly polarized quasars (HPQ)
and/or optically violently variable (OVV) quasars. The reason is that there is
increasing evidence that these categories and the flat-spectrum radio quasars,
which reflect different empirical definitions, refer to more or less the same
class of sources. That is, the majority of flat-spectrum radio quasars tend to
show rapid variability and high polarization.} To appreciate the relevance of
the latter point, we will first have to tackle the subject of AGN
classification.

\subsection{The AGN Zoo}\label{sec:class}

The large number of classes and subclasses appearing in AGN literature could
disorientate physicists or even astronomers working in other fields. A
simplified classification, however, can be made based on only two parameters,
that is radio-loudness and the width of the emission lines, as summarized in
Table 1 (see, e.g., Urry and Padovani 1995).

\begin{table}[t]
\vspace{0.4cm}
\begin{center}
\begin{tabular}{|l|l|l|l|}
\hline
{\bf Radio Loudness} & \multicolumn{3}{c|}{\bf Optical Emission Line
Properties} \\ \hline

 & {\bf Type 2} (Narrow Line) & {\bf Type 1} (Broad Line) & {\bf Type 0}
(Unusual) \\
 & & & \\
Radio-quiet: & Seyfert 2 & Seyfert 1 &  \\
 & & & \\
             &           & QSO       &  \\
 & & & \\
\hline
 & & & \\
Radio-loud: & NLRG $\cases {{\rm FR~I} \cr ~ \cr {\rm FR~II} \cr}$ &
BLRG & Blazars $\cases {{\rm BL~Lacs} \cr ~ \cr {\rm (FSRQ)} \cr}$ \\

            &    & SSRQ & \\
            &    & FSRQ & \\
 & & & \\
\hline
 & \multicolumn{3}{c|}{decreasing angle to the line of sight $\longrightarrow$} \\
\hline
\end{tabular}
  \caption{\em {AGN Taxonomy: A Simplified Scheme}}
\end{center}
\end{table}

Although it was the strong radio emission of some quasars that led to their
discovery about 35 years ago, it soon became evident that the majority of
quasars were actually radio-quiet, that is most of them were not detected by
the radio telescopes of the time. It then turned out that radio-quiet did not
mean radio-silent, that is even radio-quiet AGN can be detected in the radio
band. Why then the distinction? If one plots radio luminosity versus optical
luminosity for complete samples of optically selected sources, it looks like
there are two populations, the radio quiet one having, for the same optical
power, a radio power which is about $3 - 4$ orders of magnitudes smaller.
The distribution of the luminosity ratio $L_{\rm r}/L_{\rm opt}$ for complete
samples, including the upper limits on the radio luminosity, appears to be
bimodal, with a dividing line at a value $L_{\rm r}/L_{\rm opt} \sim 10$
(e.g., Stocke et al. 1992; both luminosities are in units of power/Hz).
It would therefore be better to call the two classes ``radio-strong'' and
``radio-weak'' but the original names are still used. Note that only about $10
- 15\%$ of AGN are radio-loud.

The other main feature used in AGN classification is the width, in case they
are present, or the absence, of emission lines. These are produced by the
recombination of ions of various elements (most notably H, He, C, N, O, Ne,
Mg, Fe). Their width is due to the Doppler effect thought to result from
more or less ordered motion around the central object. AGN are then divided in
Type 1 (broad-lined) and Type 2 (narrow-lined) objects according to their
line-widths, with 1000 km/s (full width half maximum) being the dividing
value. Some objects also exist with unusual emission line properties, such as
BL Lacs, which have very weak emission lines with typical equivalent widths (a
measure of the ratio between line and continuum luminosity) $< 5$ \AA.

As illustrated in Table 1, we then have radio-quiet Type 2 and Type 1 AGN,
that is Seyfert 2 galaxies and Seyfert 1 galaxies/radio-quiet quasars (QSO)
respectively. Radio-loud Type 2 AGN are radio galaxies (sometimes also called
narrow-line radio galaxies [NLRG] to distinguish them from the broad-lined
ones), classified as Fanaroff-Riley (Fanaroff and Riley 1974) I and II (FR I
and II) according to their radio morphology (connected with their radio
power), while radio-loud Type 1 AGN are broad-line radio galaxies (BLRG) and
radio quasars.  Finally, radio-loud sources with very weak emission lines are
known as BL Lacertae objects, from the name of the class prototype, which was
originally presumed to be a variable star in the Lacerta constellation.

Concentrating on the radio-loud sources, to which most of this paper is
devoted, the BLRG are, at least in my view, simply local versions of radio
quasars where we can detect the host galaxy, as Seyfert 1 galaxies are local
versions of QSO (the possible reasons why we do not see the high-redshift
counterparts of Seyfert 2 galaxies are discussed by Padovani 1998). Radio
quasars are generally divided into steep-spectrum radio quasars (SSRQ) and
flat-spectrum radio quasars (FSRQ), according to the value of their radio
spectral index at a few GHz ($\alpha_{\rm r} = 0.5$ being usually taken as the
dividing line, with $f_{\nu} \propto \nu^{-\alpha}$). This distinction
reflects the size of the radio emitting region. In fact, radio emission in
these sources is explained in terms of synchrotron radiation (that is
radiation from relativistic particles moving in a magnetic field), which for
extended regions has a relatively steep spectrum ($\alpha_{\rm r} \sim 0.7$).
On the other hand, nuclear, compact emission has a flatter spectrum, thought
to be the result of the superposition of various self-absorbed components. The
flat radio spectrum then indicates that nuclear emission dominates over the
more extended emission, generally associated with the so-called
``radio-lobes.''  In fact, flat-spectrum quasars are generally core-dominated
in the radio band, that is emission from the core is much stronger than
emission from the extended regions, unlike for example SSRQ or narrow-line
radio galaxies which are both lobe-dominated. Note that even though FSRQ have
strong broad lines they are also included in the ``Type 0'' column in Table 1
because their multifrequency spectra are dominated by non-thermal emission as
in BL Lac objects.

\subsection{Unified Schemes}\label{sec:unif}

\begin{figure}[p]
\hskip 1.5truecm
\psfig{figure=fig1_mor.ps,height=18truecm}
\caption{A schematic (and highly idealized) diagram of the current paradigm
for radio-loud AGN (not to scale). Surrounding the central black hole is a
luminous accretion disk.  Broad emission lines are produced in clouds (dark
spots) orbiting above the disk and perhaps by the disk itself. A thick dusty
torus (or warped disk) obscures the broad-line region from transverse lines of
sight; some continuum and broad-line emission can be scattered into those
lines of sight by hot electrons (black dots) that pervade the region. A hot
corona above the accretion disk may also play a role in producing the hard
X-ray continuum.  Narrow lines are produced in clouds (grey spots) much
farther from the central source. Radio jets, shown here as the diffuse jets
characteristic of low-luminosity, or FR~I-type, radio sources, emanate from
the region near the black hole, initially at relativistic speeds (Urry and
Padovani 1995; copyright Astronomical Society of the Pacific, reproduced with
permission).}
\label{fig:unif}
\end{figure}

All this might seem complicated, but in recent years we have developed a
consistent scenario which at least explains the Type 0/1/2 distinction. We
have in fact come to understand that some classes of apparently different (and
therefore classified under different names) AGN might actually be
intrinsically the same class of objects seen at different angles with the line
of sight (see for example Antonucci 1993 and Urry and Padovani 1995 and
references therein).

The main idea, based on various observations and summarized in Fig. 1, is that
emission in the inner parts of AGN is highly anisotropic. The current paradigm
for AGN includes a central engine, possibly a massive black hole, surrounded
by an accretion disk and by fast-moving clouds, probably under the influence
of the strong gravitational field, emitting Doppler-broadened lines. More
distant clouds emit narrower lines.  Absorbing material in some flattened
configuration (usually idealized as a toroidal shape) obscures the central
parts, so that for transverse lines of sight only the narrow-line emitting
clouds are seen (Type 2 AGN), whereas the near-infrared to soft-X-ray nuclear
continuum and broad-lines are visible only when viewed face-on (Type 1
AGN). In radio-loud objects we have the additional presence of a relativistic
jet, roughly perpendicular to the disk, which produces strong anisotropy and
amplification of the continuum emission (``relativistic beaming''), which I
discuss in more detail in Sect. \ref{subsec:prop}. For reasons still unclear,
BL Lac objects have extremely weak emission lines, and their continuum is
very strong and non-thermal (i.e., due to synchrotron and, at high
energies, inverse Compton emission or perhaps hadronic processes).

This axisymmetric model of AGN implies widely different observational
properties (and therefore classifications) at different aspect angles. Hence
the need for ``Unified Schemes'' which look at intrinsic, isotropic
properties, to unify fundamentally identical (but apparently different)
classes of AGN. Seyfert 2 galaxies have therefore been ``unified'' with
Seyfert 1 galaxies, whilst low-luminosity (FR I) and high-luminosity (FR II)
radio galaxies have been unified with BL Lacs and radio quasars respectively
(see Antonucci 1993 and Urry and Padovani 1995 and references
therein). In other words, BL Lacs are thought to be FR I radio galaxies with
their jets at relatively small ($\simlt 20 - 30^{\circ}$) angles w.r.t. the
line of sight. Similarly, we believe FSRQ to be FR II radio galaxies oriented
at small ($\simlt 15^{\circ}$) angles, while SSRQ should be at angles
in between those of FSRQ and FR II's ($15 \simlt \theta \simlt 40^{\circ}$).
Blazars are then a special class of AGN which we think have their jets
practically oriented towards the observer.

In general, different AGN components are important at different wavelengths.
Na\-me\-ly: 1. the jet emits non-thermal radiation, via electromagnetic
(synchrotron and inverse Compton) and perhaps hadronic processes, all the way
from the radio to the $\gamma$-ray band (Mastichiadis 1998; Dar 1998);
2. the accretion disk probably emits thermal radiation, peaked in
optical/ultraviolet/soft-X-ray band; 3. the obscuring material (torus) will
emit predominantly in the infrared. These different components are apparent,
for example, in the multifrequency spectrum of 3C 273 (Lichti et al. 1995)
the first quasar to be discovered and one of the best studied.

At this point one might ask: what has all this to do with $\gamma$-ray
emission? The answer in the next section.

\section{Blazars as $\gamma$-ray Sources}\label{sec:role}

According to Unified Schemes, blazars are that special class of radio-loud AGN
with their jets pointing more or less towards the observer, and therefore
constitute a relatively rare class of objects. Radio-loud AGN make up only
$\sim 10 - 15\%$ of all AGN (e.g., Kellermann et al. 1989), while a generous
upper limit to the fraction of blazars amongst radio sources is 50\% (as
inferred, for example, from the fraction of FSRQ and BL Lacs in the 1 Jy
catalogue [Stickel et al. 1994] which, being a high-frequency radio catalogue,
is biased towards flat-spectrum sources). It then follows that blazars make up
at most 5\% of all AGN, but more likely even less than that.

Mukherjee et al. (1997) have identified 51 high-confidence EGRET sources
(mainly from the Second EGRET catalogue; Thompson et al. 1995, 1996) with AGN,
all of them blazars. If the probability of detecting an AGN with EGRET were
independent of the class, then in this list we would expect at maximum 3
blazars, with most sources being associated with radio-quiet AGN. Instead, we
have 100\% blazars and 0\% other sources. In particular, no radio-quiet AGN
has been detected so far by EGRET. Note that blazar $\gamma$-ray luminosities
are in the range $10^{45} - 10^{49}$ erg/s (under the assumption of
isotropy) and in many cases the output in
$\gamma$-rays dominates the total (bolometric) luminosity.

To find out what is so special about blazars we need to have a closer look at
their properties.

\subsection{Blazar Properties and Relativistic Beaming}\label{subsec:prop}

The main properties of blazars can be summarized as follows:

\begin{itemize}
\item radio loudness;
\item rapid variability (high $\Delta L/\Delta t$);
\item high and variable optical polarization ($P_{\rm opt} > 3\%$);
\item smooth, broad, non-thermal continuum;
\item compact, flat-spectrum radio emission ($f_{\rm core} \gg f_{\rm
extended}$);
\item superluminal motion in sources with multiple-epoch Very Large Baseline
Interferometry (VLBI) maps.
\end {itemize}

The last property might require some explanation. The term ``superluminal
motion'' describes proper motion of source structure (traditionally mapped at
radio wavelengths) that, when converted to an apparent speed $v_{\rm app}$,
gives $v_{\rm app} > c$. This phenomenon occurs for emitting regions moving at
very high (but still subluminal) speeds at small angles to the line of sight
(Rees 1966). Relativistically moving sources ``run after'' the photons they
emit, strongly reducing the time interval separating any two events in the
observer's frame and giving the impression of faster than light motion.

Analytically, the observed transverse velocity of an emitting blob, $v_{\rm a}
= \beta_{\rm a} c$, is related to its true velocity, $v = \beta c$, and the
angle to the line of sight $\theta$ by
\begin{equation}
\beta_{\rm a} = {\beta \sin\theta \over 1 - \beta \cos\theta} \quad.
\label{eq:beta_app}
\end{equation}
It can be shown that if $\beta > 1/\sqrt{2} \simeq 0.7$, then for some
orientations superluminal motion (that is, $\beta_{\rm a} > 1$) is
observed. The maximum value of the apparent velocity, $\beta_{\rm a,max} =
\sqrt{\gamma^2 -1}$, where $\gamma = (1-\beta^2)^{-1/2}$ is the Lorentz
factor, occurs when $\cos \theta = \beta$ or $\sin \theta = \gamma^{-1}$. This
implies a minimum value for the Lorentz factor $\gamma_{\rm min} =
\sqrt{\beta_{\rm a}^2 + 1}$. For example, if one detects superluminal motion
in a source with $\beta_{\rm a} = 5$, the Lorentz factor responsible for it
has to be at least 5.1.

All these properties are consistent with {\it relativistic beaming}, that is
with bulk relativistic motion of the emitting plasma towards the observer.
There are by now various arguments in favor of relativistic beaming in
blazars, summarized for example by Urry and Padovani (1995). Beaming has
enormous effects on the observed luminosities. Adopting the usual definition
of the relativistic Doppler factor $\delta = [\gamma (1 - \beta \cos
\theta)]^{-1}$
and applying
simple relativistic transformations, it turns out that the {\it observed}
luminosity at a given frequency is related to the {\it emitted} luminosity in
the rest frame of the source via

\begin{equation}
L_{\rm obs} = \delta^p L_{\rm em} \quad ,
\end{equation}

with $p = 2+\alpha$ or $3+\alpha$ respectively in the case of a continuous jet
or a moving sphere (Urry and Padovani 1995; $\alpha$ being the spectral
index), although other values are also possible (Lind and Blandford 1985).
For $\theta \sim 0^{\circ}$, $\delta \sim 2 \gamma$ (Fig. 2)
and the observed luminosity can be amplified by factors of thousands (for
$\gamma \sim 5$ and $p \sim 3$, which are typical values). That is, for jets
pointing almost towards us we can overestimate the emitted luminosity
typically by three orders of magnitude.  Apart from this amplification,
beaming also gives rise to a strong collimation of the radiation, which is
larger for higher $\gamma$ (Fig. 2): $\delta$ decreases by a
factor $\sim 2$ from its maximum value at $\theta \sim 1/\gamma$ and
consequently the inferred luminosity goes down by $2^p$. For example, if
$\gamma \sim 5$ the luminosity of a jet pointing $\sim 11^{\circ}$ away from
our line of sight is already about an order of magnitude smaller (for $p = 3$)
than that of a jet aiming straight at us.

\begin{figure}
\hskip 2.5truecm
\psfig{figure=fig2_mor.ps,width=9.9truecm}
\caption{The dependence of the Doppler factor $\delta$ on the angle to the
line of sight. Different curves correspond to different Lorentz factors: from
the top down, $\gamma = 15$ (solid line), $\gamma = 10$ (dotted line), $\gamma
= 5$ (short-dashed line), $\gamma = 2$ (long-dashed line).}
\label{fig:doppler}
\end{figure}

All this is very relevant to the issue of $\gamma$-ray emission from blazars.
In fact, if blazars were not beamed, we would not see any $\gamma$-ray photons
from them! The qualitative explanation is relatively simple: in sources as
compact as blazars all $\gamma$-ray photons would be absorbed through
photon-photon collisions with target photons in the X-ray band. The end
product would be electron-positron pairs. But if the radiation is beamed then
the luminosity/radius ratio, which is the relevant parameter, is smaller by a
factor $\delta^{p+1}$ and the $\gamma$-ray photons manage to escape from the
source. More formally, it can be shown (Maraschi et al. 1992) that the
condition that the optical depth to photon-photon absorption $\tau_{\gamma
\gamma}(x)$ is less than 1 implies (under the assumption that the X-ray and
$\gamma$-ray photons are produced in the same region)

\begin{equation}
\delta > C \left({L_{48}\over \Delta t_{\rm d}}\right)^{1/
(p+1)} \left({x \over 10^4}\right)^{\alpha_{\rm x}/(p+1)} \quad ,
\end{equation}

where $L_{48} \equiv L_{\gamma}/(10^{48}$ erg/s), $\Delta t_{\rm d}$ is the
$\gamma$-ray variability time scale in days (which is used to estimate the
source size), $x \equiv h\nu/m_{\rm e} c^2$, $\alpha_{\rm x}$ is the X-ray
spectral index, and $C$ is a numerical constant $\approx 10$. In other words,
transparency for the $\gamma$-ray photons {\it requires} a relatively large
Doppler factor for most blazars (Dondi and Ghisellini 1995) and therefore
relativistic beaming.

\section{The Importance of Being a Blazar}\label{sec:impo}

We can now address the main question of this paper: why have blazars been
detected by EGRET? There are at least three reasons, which have to
do with the fact that blazars are characterized by:

\begin{enumerate}
\item high-energy particles, which can produce GeV photons;
\item relativistic beaming, to avoid photon-photon collision and amplify the
flux;
\item strong non-thermal (jet) component.
\end{enumerate}

Point number one is obvious. We know that in blazars synchrotron emission
reaches at least the infrared/optical range, which reveals the presence of
high-energy electrons which can produce $\gamma$-rays via inverse Compton
emission (although hadronic processes can also be important or even dominant;
Mannheim 1993). Very recent BeppoSAX observations have shown that synchrotron
emission can actually reach the X-ray band, namely around 10 keV for 1ES
2344+514 (Giommi, Padovani, and Perlman 1998) and 100 keV for MKN 501 (Pian et
al. 1998). Pian et al. fit a synchrotron-self-Compton (SSC) model to the
multifrequency spectrum of MKN 501 and infer a maximum value for the Lorentz
factor of the electrons $\gamma_{\rm max} \sim 3 \times 10^6$. With such high
values of $\gamma_{\rm max}$, Pian et al. have been able to reproduce the
observed TeV flux of this source (see Sect. \ref{sec:tev}) with an SSC model.
Point number two is vital, as described in the previous section, not only to
enable the $\gamma$-ray photons to escape from the source, but also to amplify
the flux and therefore make the source more easily detectable. Point number
three is also very important. $\gamma$-ray emission is clearly non-thermal
(although we still do not know for sure which processes are responsible for
it) and therefore related to the jet component. The stronger the jet
component, the stronger the $\gamma$-ray flux.

Having understood why blazars have been detected by EGRET, one could also ask:
why have not {\it all} blazars been detected? Many blazars with radio
properties similar to those of the detected sources, in fact, still have only
upper limits in the EGRET band. This problem has been addressed, for example,
by von Montigny et al. (1995) who suggest as possible solutions variability
(only objects flaring in the $\gamma$-ray band can be detected) or a
$\gamma$-ray beaming cone which either points in a different direction or is
more narrowly beamed than the radio one (see also Salamon and Stecker 1994 and
Dermer 1995). These can certainly be viable explanations, but one should also
note that even a moderate dispersion in the values of the parameters required
for $\gamma$-ray emission described above (particle energy, Doppler factor,
and non-thermal component strength) can easily imply the non-detection of some
sources and the detection of others with similar radio properties.

Do the points discussed above also explain why EGRET has not detected any of
the more numerous radio-quiet AGN? Yes, although not all of them might be
essential in this case. As radio-emission (at least in radio-loud AGN) is
certainly non-thermal, while the optical/ultraviolet emission might be thermal
emission associated with the accretion disk (at least in radio-quiet AGN),
then the ratio $L_{\rm r}/L_{\rm opt}$ could actually be related to the ratio
$L_{\rm non-thermal}/L_{\rm thermal}$. Furthermore, $L_{\gamma}$ seems to
scale with $L_{\rm r}$, although the details of this dependence are still
under debate (see e.g., Mattox et al. 1997 and references therein). It
could then be that even radio-quiet AGN are $\gamma$-ray emitters, although
scaled down by their ratio between radio and optical powers, that is at a
level $3 - 4$ orders of magnitude below that of blazars, with typical fluxes
$F_{\gamma} \approx 10^{-10}$ photons cm$^{-2}$ s$^{-1}$. In other words,
radio-quiet AGN would fulfill requirements number one and two\footnote{One
could argue that the relativistic beaming requirement would probably not be
very important in radio-quiet AGN as the luminosity/radius ratio in the
$\gamma$-ray band in these sources would be much lower anyway and $\gamma$-ray
photons would escape even without beaming. However, GeV photons collide
preferentially with X-ray photons, which are plentiful in these sources: some
beaming might then be required for the (putative) $\gamma$-ray emission in
radio-quiet AGN as well.} but not number three.

Alternatively, it could be that for some reason the emission mechanisms at
work in radio-loud sources are simply not present in the radio-quiet ones,
either because there is no jet at all in radio-quiet AGN or because, for
example, there is no accelerating mechanism. In this case, either condition
number one or number three (or both) would be missing (number two would now be
unimportant), and no $\gamma$-ray emission would be expected.

Unfortunately, it might not be possible to test these two alternatives on the
basis of $\gamma$-ray data for some time: even in the first case, in fact, the
expected $\gamma$-ray fluxes are below the sensitivity of currently planned
future $\gamma$-ray missions, like GLAST (Morselli 1998). If radio-quiet AGN
emit $\gamma$-rays, however, GLAST might be able to detect some of the
nearest sources.

\section{AGN as TeV Sources}\label{sec:tev}

So far, only emission up to a few GeV has been considered. However, TeV
astronomy is now in full swing, as we have heard at this meeting, and it is
therefore interesting to consider the situation at these energies.

Three extragalactic sources have been detected at $E > 0.3$ TeV (Lorenz 1998)
and these are all BL Lacs. That is, even at energies above those of EGRET (and
exactly for the same reasons) the only $\gamma$-ray emitting AGN are still
blazars! The difference here is that, unlike the situation in the EGRET band
where the majority of detected blazars are flat-spectrum radio quasars, only
BL Lacs have been detected so far. Furthermore, these BL Lacs are the three
nearest confirmed BL Lacs in the recent catalogue by Padovani and Giommi
(1995) namely MKN 421 (redshift $z = 0.031$), MKN 501 ($z = 0.055$) and 1ES
2344+514 ($z = 0.044$). The fact that only relatively nearby BL Lacs have been
detected is probably related to absorption of TeV photons by the infrared
background (see, e.g., Biller et al. 1995 and references therein).
Additionally, the position of the synchrotron peak might anticorrelate with
bolometric luminosity in blazars (Fossati et al. 1997). One would then expect
the less luminous, and therefore nearer, objects to have the synchrotron peak
at UV/X-ray energies and the peak of the inverse Compton emission (within the
SSC model) in the TeV band. These objects, therefore, would be, on average,
stronger TeV sources.

Why have no flat-spectrum radio quasars been detected at TeV energies? These
sources are typically at higher redshifts and so the effect on them of the
cosmological absorption by infrared photons is more severe. However, there are
at least four strong radio sources classified as flat-spectrum quasars at $z <
0.1$, including 3C 120 and 3C 111, the latter having been looked at by the
Whipple experiment, with negative results (Kerrick et al. 1995; 3C 111,
however, although superluminal [Vermeulen and Cohen 1994] is
lobe-dominated [$f_{\rm
core}/f_{\rm extended} \simeq 0.2$; Hes et al. 1995], which suggests it is an
unlikely blazar.  Also, it has not been detected by EGRET.)

This is certainly small number statistics and definite conclusions should only
be drawn after a larger number of relatively local flat-spectrum radio quasars
have been observed at TeV energies. However, the non-detection of
flat-spectrum radio quasars could simply mean that {\it internal} absorption
is significant in these sources. In fact, the cross-section for photon-photon
interaction for $\sim 1$ TeV photons is maximum at $\sim 10^{14}$ Hz or $\sim
2.5 \mu$ and quasars have a larger photon density than BL Lacs at these
frequencies, be it emission from the obscuring torus\footnote{It is not clear
if the presence of an obscuring torus is required in BL Lacs as well as in
radio quasars: see discussion in Urry and Padovani (1995) and Padovani (1997)
and references therein.} or even the accretion disk (see e.g., Protheroe and
Biermann 1997).

The new, more sensitive projects which are underway in the field of TeV
astronomy (e.g., Krennrich 1998) will certainly shed light on these issues and
on the processes which are responsible for the GeV/TeV emission in blazars.
The origin of such emission is in fact still debated as being due to inverse
Compton radiation, either SSC or Comptonization of external radiation (e.g.,
Sikora et al. 1994), or to hadronic processes, that is pion production from
accelerated protons, with subsequent pion decay and $\gamma$-ray production
(e.g., Mannheim 1993). In the latter case, these so-called ``proton blazars''
would also be sources of cosmic rays and of neutrinos and could be detectable
by planned km$^3$ neutrino detectors (Halzen and Zas 1997).

\section{Summary}\label{sec:sum}
The main conclusions are as follows:

\begin{enumerate}
\item The only AGN detected at GeV and even TeV energies are blazars,
that is a special class of sources which includes BL Lacertae objects and
radio quasars with a relatively flat radio spectrum. These sources are thought
to have their jets moving at relativistic speeds almost directly towards the
observer, a phenomenon which goes under the name of ``relativistic beaming''
and causes strong amplification and collimation of the radiation in our rest
frame.
\item As blazars make up {\it at most} 5\% of all AGN, they must have some
peculiar characteristics which favor their $\gamma$-ray detection. I have
shown that the presence of high-energy particles, relativistic beaming, and a
strong non-thermal (jet) component play, in fact, a fundamental role in making
these sources detectable at $\gamma$-ray energies.
\item $\gamma$-ray missions $\sim 1,000$ times more sensitive than EGRET
might
also detect the bulk of the more common radio-quiet AGN, under the assumption
that they also possess, on much smaller scales, a non-thermal engine.
\item TeV astronomy, a very young branch of astronomy which has already
produced some very exciting results, will likely play an important role in the
near future in constraining blazar models. More sensitive TeV telescopes
are clearly needed.
\end{enumerate}

In summary, there exists a tight connection between unified schemes and
$\gamma$-ray emission in AGN, as they both depend on relativistic beaming, the
former as a mechanism to produce a strong angle dependence of the observed
properties, the latter as a way to let $\gamma$-ray photons escape from the
source.

\section{References}

\hang\noindent Antonucci, R.: 1993, {Annual Review of Astronomy and Astrophysics} {\bf 31}, p. 473.

\hang\noindent Biller, S., et al.: 1995, {Astrophysical Journal} {\bf 445}, p. 227.

\hang\noindent Dar, A.: 1998, these proceedings.

\hang\noindent Dermer, C.D.: 1995, {Astrophysical Journal (Letters)} {\bf 446}, p. 63.

\hang\noindent Dondi, L., Ghisellini, G.: 1995, {Monthly Notices of the Royal Astronomical Society} {\bf 273}, p. 583.

\hang\noindent Fanaroff, B.L, Riley, J.M.: 1974, {Monthly Notices of the Royal Astronomical Society} {\bf 167}, 31p.

\hang\noindent Fossati, G., Celotti, A., Ghisellini, G., Maraschi, L.: 1997, {Monthly Notices of the Royal Astronomical Society}
{\bf 289}, p. 136.

\hang\noindent Giommi, P., Padovani, P., Perlman, E.: 1998, {Monthly Notices of the Royal Astronomical Society}, submitted.

\hang\noindent Halzen, F., Zas, E.: 1997, {Astrophysical Journal} {\bf 488}, p. 669.

\hang\noindent Hes, R., Barthel, P.D., Hoekstra, H.: 1995, {Astronomy and Astrophysics} {\bf 303}, p. 8.

\hang\noindent Kellermann, K.I., Sramek, R., Schmidt, M., Shaffer, D.B., Green, R.:
1989,  {Astronomical Journal} {\bf 98}, p. 1195.

\hang\noindent Kerrick, A.D., et al.: 1995, {Astrophysical Journal} {\bf 452}, p. 588.

\hang\noindent Krennrich, F.: 1998, these proceedings.

\hang\noindent Lichti, G.G., et al.: 1995, {Astronomy and Astrophysics} {\bf 298}, p. 711.

\hang\noindent Lind, K.R., Blandford, R.D.: 1985, {Astrophysical Journal} {\bf 295}, p. 398.

\hang\noindent Lorenz, E.: 1998, these proceedings.

\hang\noindent Mannheim, K.: 1993, {Astronomy and Astrophysics}, {\bf 269}, p. 67.

\hang\noindent Maraschi, L., Ghisellini, G., Celotti, A.: 1992, {Astrophysical Journal (Letters)} {\bf 397}, p. 5.

\hang\noindent Mastichiadis, A.: 1998, these proceedings.

\hang\noindent Mukherjee, R., et al.: 1997, {Astrophysical Journal} {\bf 490}, p. 116.

\hang\noindent Mattox, J.R., Schachter, J., Molnar, L., Hartman, R.C., Patnaik, A.R.:
1997, {Astrophysical Journal} {\bf 481}, p. 95.

\hang\noindent von Montigny, C., et al.: 1995, {Astronomy and Astrophysics} {\bf 299}, p. 680.

\hang\noindent Morselli, A.: 1998, these proceedings.

\hang\noindent Padovani, P.: 1997, From the Micro- to the Mega-Parsec,
eds. A. Comastri, T. Venturi, M. Bellazzini, {\SAIT} {\bf 68}, p. 47.

\hang\noindent Padovani, P.: 1998, New Horizons from Multi-Wavelength Sky
Surveys, eds. B. McLean et al., p. 257.

\hang\noindent Padovani, P., Giommi, P.: 1995, {Monthly Notices of the Royal Astronomical Society} {\bf 277}, p. 1477.

\hang\noindent Pian, E., et al.: 1998, {Astrophysical Journal (Letters)} {\bf 492}, p. 17.

\hang\noindent Protheroe, R.J., Biermann, P.L.: 1997, {\APP} {\bf 6}, p. 293.

\hang\noindent Rees, M.J.: 1966, {Nature} {\bf 211}, p. 468.

\hang\noindent Salamon, M.H., Stecker, F.W.: 1994, {Astrophysical Journal (Letters)} {\bf 430}, p. 21.

\hang\noindent Sikora, M., Begelman, M.C., Rees, M.J.: 1994, {Astrophysical Journal}, {\bf 421}, p. 153.

\hang\noindent Stickel, M., Meisenheimer, K., K\"uhr, H.: 1994, {Astronomy and Astrophysics Supplement Series} {\bf 105}, p. 211.

\hang\noindent Stocke, J.T., Morris, S.L., Weymann, R.J., Foltz, C.B.: 1992,
{Astrophysical Journal} {\bf 396}, p. 487.

\hang\noindent Thompson, D.J., et al.: 1995, {Astrophysical Journal Supplement Series} {\bf 101}, p. 259.

\hang\noindent Thompson, D.J., et al.: 1996, {Astrophysical Journal Supplement Series} {\bf 107}, p. 227.

\hang\noindent Urry, C.M., Padovani, P.: 1995, {Publications of the Astronomical Society of the Pacific} {\bf 107}, p. 803.

\hang\noindent Vermeulen, R.C., Cohen, M.H.: 1994, {Astrophysical Journal} {\bf 430} p. 467.

\end{document}.


Extragalactic Synchrotron Sources at Low Frequencies

Abstract: The LWA will be well suited to address many important questions about the physics and astrophysics of extragalactic synchrotron sources. Good low-frequency data will enable major steps forward in our understanding of radio galaxy physics, of the plasma in clusters of galaxies, and of active objects in the high-redshift universe. Such data will also be important in answering some basic questions about the physics of synchrotron-emitting plasmas.

\section{General Introduction}

Synchrotron emission is one of our most important tools for understanding
the active extragalactic universe.  Energetic processes in many
astrophysical settings create a relativistic, magnetized plasma which
we can observe {\it via} its synchrotron emission.  We can study  the
high-energy activity directly, for instance in the creation of a radio
jet.  We can use radio sources as probes of larger questions, for instance
learning about the plasma atmosphere surrounding a radio tail by studying
the dynamics of the radio tail, or  understanding the dynamical state of the
dark matter in a cluster of galaxies by studying the synchrotron halo in
the cluster.  We can study the population of massive black holes in
the early universe through detections of high-redshift synchrotron
sources.

High quality images at low frequencies will be critical to all of
these applications.  In this paper we discuss some important
contributions the LWA will make.  One contribution is simple
exploration.  Many interesting sources are steep spectrum, and
therefore can best  be detected and studied at low
frequencies. Multi-frequency imaging will be another important
contribution, combining high-resolution images from the LWA with those
from  higher-frequency  instruments such as the EVLA.   We  need good
images  across a broad frequency range, including the very low
frequencies which the LWA will make available, in order to understand
the nature and dynamics of radio jets on large and small scales.
Finally, some fundamental questions regarding particle acceleration
and the physics of synchrotron emission are still unresolved;   good
low-frequency data will be an important part of the answers.

\section{Basic Synchrotron Physics}

If we want to use synchrotron emission to study  active
extragalactic sources, we need to understand the physics of the emitting
plasma.  Two important issues (at least) remain obscure, the low-frequency
``injection'' spectrum and the high-frequency turndown known as
 ``spectral aging''.
We remind the reader that the electron energy, $\gamma m_e c^2$, maps
to the emission frequency in synchrotron radiation
as $\nu_{sy} \propto \gamma^2 B$.  Thus, for a given magnetic field,
lower emission frequencies correspond to lower electron energies.

\subsection{Low Frequency Spectra}

We still do not know  how relativistic electrons are accelerated in
radio sources.  Options include shock acceleration, turbulent acceleration,
and electrodynamic acceleration (such as in a reconnection layer).  These
different mechanisms  have different signatures in the energy range
 and spectrum of the accelerated
electrons.  In particular, whether or not the low-energy electron spectrum
is a power law is often assumed but has not been established from the data.
Some authors have disputed this assumption based on the data
(Katz-Stone {\it et al.\,}1993), and simple models of acceleration can lead to
alternatives ({\it e.g.}, Borovsky \& Eilek 1986, Brunetti {\it et
  al.\,}2001).  The particle distribution at low
energies is also important to the energy budget of the system;
small differences in the usual assumptions (power laws and cutoffs)
can have large consequences on the interpretation ({\it e.g.}, Harris 2005).

Because radiative losses affect the highest energy electrons
first, the low frequency synchrotron spectrum, coming from the
low energy electrons, should be the most ``pristine''.  If we can
measure the low-frequency spectrum  {\it in a homogeneous region}, we
can  determine both the spectrum and low-energy cutoff of the
electrons.  This will allow critical tests of the different acceleration
models.   We note that the only unambiguous sign of a low-energy
cutoff in the electron distribution is the characteristic  $\nu^{1/3}$
low-frequency spectrum;  steeper low-frequency spectra are due to
absorption effects.

\subsection{Spectral Aging}

Another important question is the interpretation of
 high-frequency spectral steepening.   This common effect
is illustrated by 3C465 in Figure 1, which also shows how different a radio
source can look at low and high frequencies.   Other striking
examples include the lobes of Cyg A (Carilli {\it et al.\,}1991) and the long
tails recently detected in Hyd A (Lane {\it et al.\,}2004).

\begin{figure}[ht]

\begin{center}
\scalebox{0.224}{\includegraphics{Eilek1fig1a.ps}}
\scalebox{0.224}{\includegraphics{Eilek1fig1b.ps}}
\scalebox{0.224}{\includegraphics{Eilek1fig1c.ps}}
\end{center}
\vspace{-2ex}
\caption{Three VLA images of 3C465:  left, 327 MHz;  center, 1.4 GHz;
  right, 4.8 GHz;  from Eilek \& Owen 2002 and references
  therein. This shows the  rapid decay of the radio tails at high
  frequency, that is the steepening of the radio spectrum going out
  along the tail.  This effect is usually ascribed to synchrotron
  losses of the electrons, but this interpretation may be too
  simple.  These images
  also illustrate how much of a source can be hidden at high
  frequencies; without good low-frequency   images we cannot
  understand the full story  of the radio source.  }
\end{figure}

Spectral steepening is often interpreted very simply,
namely, that the high-energy electrons have been depleted by
synchrotron losses (thus the name ``spectral aging'').   The spectrum
steepens at the frequency radiated by electrons whose
radiative lifetime equals the age of the source.
 If the magnetic field is known (or assumed, say
from equipartition analyses), the age of the radiating plasma can be
derived.  However, the results here are often inconsistent with other
estimates which find a longer age for the radio source ({\it e.g.}
O'Donoghue {\it et al.\,}1993;  Blundell \& Rawlings 2000;  Eilek
2004).

Because source ages are important to many aspects of source modelling,
the simple interpretation of spectral  steepening must be tested,
not just assumed.  We suggest the best way to do this is to develop
careful dynamic models of some nearby, bright radio sources. Doing this
requires good images over a broad frequency range, {\it including the
lowest possible frequencies}:  in what parts of the source does the
spectrum steepen, and what are the likely plasma dynamics there? Is simple
aging plausible, or must we consider {\it in situ} acceleration or
magnetic field variability, or both?

\subsection{Cautionary Note}

Both of these analyses will be  complicated by inhomogeneities in the
source, especially the magnetic field.  On small scales, we expect the
field to be filamented; on larger scales the dynamics of the source
will determine the magnetic field structure.  We know that magnetic
inhomogeneities can complicate  interpretation of synchrotron spectra
({\it e.g.}, a single electron energy in a power law magnetic field
distribution  can also give a power law synchrotron spectrum;  Eilek
\& Arendt 1996).   In some applications this complication can be dealt
with in terms of mean fields or a two-field approximation (Eilek {\it
  et al.\,}1997);  but
observations should also be designed with the highest possible
resolution and using nearby sources.

\section{Radio Galaxy Physics}

We have already pointed out that radio galaxies look quite different
at low and high frequencies.  Combining LWA images with higher
frequency images from other telescopes will enable us to probe
the life cycles and internal physics of these sources.  In addition,
the LWA by itself will add to our understanding of these galaxies.
Two examples of the latter are low-frequency absorption studies and
searches for old radio galaxies.

\subsection{Probing the Gas in the Parent Galaxy}

The LWA opens up an exciting new possibility:  studying the interstellar
medium in the gas in the parent galaxy {\it via} detection of
free-free absorption from that gas.  The $\nu^{-2.1}$
behavior of the opacity means the LWA will be able to detect absorbing
clouds that would be invisible at higher frequencies.
Free-free absorption will be easy to distinguish from a
low-energy cutoff in the electron distribution, by its different radiation
spectrum, as well as by its
  patchy nature (in a well resolved source).
Most of our knowledge of the interstellar medium (ISM) in an
elliptical galaxy comes from X-rays, which detect the hot phase
(typically $T \sim 10^7$K, $n \sim 0.1$cm$^{-3}$;  Brighenti \& Mathews
1997).  Detection of free-free absorption would allow us to study
quite different, cooler components of the the ISM.  For example,
consider cool gas at $10^4$K in pressure balance with the hot ISM.
Cool clouds in the galactic core would be  just opaque
 at 50 MHz if they are 1 pc in size;  in the lower-pressure
outer regions of the galaxy, kpc-sized clouds will be just opaque.
Such clouds would easily be detectable in the  spectrum
of the background source, and possibly in  high-resolution images.

\subsection{Where are the old sources?}

With good low-frequency data we should be able to answer a
longstanding problem: where are the old radio galaxies?  Current
models of radio source physics predict  that the radio spectrum will
steepen as the source ages.  This holds whether or not the jet remains
active,  long as the magnetic field in the source
remains somewhat steady  (Eilek \& Shore 1989).   The models also
predict that the sources should remain bright at low frequencies,
where synchrotron losses have not yet had an impact.

These models therefore predict the existence of
a large population of radio sources which
have steep spectra around 1 GHz but are bright at 100 MHz.  However,
such sources are uncommon;  we have not found nearly as many as the
models predict.  Where can they be?  Are the models seriously wrong?
One possibility is that magnetic field dissipates rapidly when the
central engine shuts down.  This could happen, for instance, if the
magnetic field is supported by jet-driven turbulence in the extended
source. If this is the case, then old, steep-spectrum radio
sources would be much fainter than the current models predict.  The
LWA will easily be able to detect these sources. If such sources are
searched for and not found, then the ``last chance'' will be gone, and
will definitely know that our general model of radio galaxies needs
serious revision.

\section{Radio Haloes in Clusters of Galaxies}

An important use of the LWA will be the  direct study of steep-spectrum
synchrotron sources.  Radio haloes in clusters of galaxies are
one interesting example. The plasma atmosphere in a
cluster of galaxies contains a significant  nonthermal component, just
as the ISM in our galaxy does.
The nonthermal component is seen {\it via} its synchrotron
emission.  As illustrated in Figure 2, the synchrotron halo is usually
a diffuse, cluster-wide source,  often obscured by small, bright radio
galaxies within the cluster.

\begin{figure}[!ht]
\begin{center}
\scalebox{0.23}{\includegraphics{Eilek1fig2a.ps}}
\scalebox{0.242}{\includegraphics{Eilek1fig2b.ps}}
\scalebox{0.21}{\includegraphics{Eilek1fig2c.ps}}
\end{center}
\vspace{-4ex}
\caption{Radio haloes reveal nonthermal plasma in the intraclutser
  medium, when they can be detected.  Abell 773, from Markovi\'c (2004), is an
  example. Left, a high-resolution 1.4 GHz VLA image showing the
  embedded galactic sources in the cluster.  Center, a low-resolution
  1.4 GHz image of the same cluster,  with the embedded sources removed
  to reveal the diffuse halo.  Right, an archival CHANDRA image of the
  cluster, illustrating that the extent of the radio halo is
  similar to that of the thermal, X-ray loud intracluster medium.}
\end{figure}

The existence of this nonthermal component
tells us the cluster is not fully relaxed.  Both
the relativistic particles and the (turbulent) magnetic field must be
re-energized over the life of the cluster.
Until recently it was thought that haloes or relics were rare, and
thus connected to some unusual event, such as a major merger.  We now
have evidence that they are more common (Markovi\'c, 2004) and that
they occur in quiescent
clusters as well as in actively merging ones.   This suggests
that the dynamic history of clusters is not as simple as has been thought;
there are ongoing energy sources  even in apparently quiescent clusters.
Fully understanding this synchrotron emission  will tell
us much about the dynamical state and history of the cluster.

In order to pursue this issue we need good images of radio haloes in
a large number of clusters, both merging and quiescent.
However, with current instruments and much hard work we are just at
the limits of being able to address this question.  The haloes are
faint, and contaminated by embedded radio galaxies in the cluster.
We need sensitive observations on a wide range of spatial scales to
sort them out.  Because haloes are usually steep spectrum, the LWA
will be both necessary and ideal for this work.

\section{The High Redshift Universe}

There is still much to learn about high-redshift synchrotron sources.
Because they are likely to have steep
radio spectra, high-$z$  sources are best studied --- or simply
detected --- at low frequencies.   We expect them to have steep
spectra because inverse Compton
losses on the microwave background increase with redshift $\propto (1
+ z)^4$, and because high-frequency spectral breaks are redshifted to
lower frequencies. This is
consistent with radio galaxies we know about now.  They are generally
the brighter, FRII-type objects, and do tend to have steep
radio spectra ({\it e.g.} de Breuck {\it et al.\,}2000).  If the nearby
universe can be used as an analogy, FRI-type radio
galaxies and strong starbursts should also exist at high $z$. We
might expect these objects to have weaker magnetic fields than the
FRII's (as they do in the nearby universe), which would shift their
spectral break to even lower frequencies.  We thus expect a large
population of faint, high-$z$ synchrotron sources is waiting to be
studied.  This work calls for the high resolution and high
sensitivity of the LWA.

One important question in this area is  the demographics of synchrotron
sources in the early universe.  How abundant are
star-forming galaxies and massive, accreting black holes at early
epochs?  To answer this question,  we must first detect the objects;
follow-up multiband spectra should enable us to separate the black
holes from the starburst galaxies.  But if the initial radio searches
are carried out only at GHz frequencies, as they must be today, we may
well be missing a large number of the objects.

Another question is  the physics of radio galaxies at high
redshift.  Are high-$z$ sources just distant  analogues of nearby
sources, subject to stronger Compton losses, or is  the situation more
complex?   FRII sources are not typical of the nearby radio galaxy
population;  tailed, FRI types are much more common.  Will more
sensitive searches, at low frequency, find large numbers of faint,
steep-spectrum FRI sources?  We pointed out above that spectral
steepening is not yet understood in nearby radio sources, particularly
FRI's.  How  will current models survive when tested in the harsher
environment at high redshift?

\section{Closing comments}

In this paper we have  discussed a few of the important
contributions the LWA can make to extragalactic synchrotron physics;
many more good ideas will no doubt appear.  Much of the
science, however,
depends on getting the highest possible resolution, down to
a few arcseconds.  We thus look forward to creative solutions to the
imaging problems which the ionosphere will impose on us.