No evidence of "gray" dust from composite quasar spectra

Abstract: Two recent studies based on composite reddened quasar spectra have indicated the presence of `gray' dust in quasar environments. This gray dust has a relatively flat extinction law in the UV, consistent with the theoretical expectation of a lack of small dust grains close to a quasar. In contrast, individual reddened quasars in the Sloan Digital Sky Survey tend to have steep extinction laws in the UV, similar to that in the SMC. We analyze the method used in determining extinction laws from composite quasar spectra in order to resolve this discrepancy. We show that quasars reddened by SMC-type dust that are present in quasar samples have a negative correlation between E_{B-V} and redshift, due to selection effects. The fact that the highest redshift quasars (which contribute to the UV part of a composite spectrum) are less extincted leads to shallower extinction in the UV. We construct a composite quasar spectrum from a simulated sample of quasars reddened by SMC-type dust and show that the extinction curve derived from the composite does not recover the intrinsic extinction law. We conclude there is no evidence of gray dust in quasar environments.

\section{Introduction}

Only $<0.01$\% of matter in galaxies is composed of dust --
coagulations of atoms and molecules with sizes ranging from a few
angstroms to a few microns (Seaquist et al. 2004). However, the
combined effects of extinction, scattering and thermal emission from
dust have a profound impact on our view of the universe. For example,
about half of the energy emitted as visible/UV radiation by stars in
distant galaxies has been reprocessed by dust into infrared emission.

Understanding the properties of this dust is vital. While dust in our
Galaxy and the Magellanic Clouds has been heavily studied (see,
e.g., the review by Draine 2003), there has been only limited work in
external galaxies, especially at high redshift. Of particular interest
is the possibility that `gray' dust is commonplace at high
redshifts. Gray dust is characterized by a flat extinction curve in
the UV and as such could be hard to identify from rest-frame UV data
alone (i.e. observed-frame optical for sources at high
redshift). Large surveys underway to measure cosmological parameters
from the Hubble diagram of Type Ia supernovae could be seriously
affected by the presence of unknown quantities of gray dust at high
redshift (Aguirre 1999).

One of the best ways to study dust extinction in distant galaxies is
via the effect that dust has on background quasars. Quasars provide a
bright UV source that can still be detected at high significance even
after significant absorption. However, a caveat is that if the dust
extinction occurs close to the active nucleus, then the unusual gas
density and radiation field could lead to atypical dust properties.

In the past, radio-loud quasar surveys have proved most fruitful for
searching for reddened quasars due to the fact that radio emission is
unaffected by dust (Webster et al. 1995; Baker \& Hunstead 1995; Gregg
et al. 2002; White et al. 2003). However the vast area and generous
color selection criteria employed by the Sloan Digital Sky Survey
(SDSS) mean that it has discovered hundreds of reddened quasars
(Richards et al. 2003). In an analysis of the continuum shapes of 9566
SDSS quasars, Hopkins et al. (2004) showed that the dust extinction is
best fit by dust similar to that in the SMC, i.e., extinction
increasing sharply in the UV. This work included 1886 SDSS quasars
with 2MASS matches which give the longer wavelength baseline necessary
to better distinguish between different extinction laws.

In stark contrast to the results from Hopkins et al. (2004), there
have been two recent studies that claim that the UV extinction curves
in quasars are `gray' or flatten in the UV (Czerny et al. 2004;
Gaskell et al. 2004). Both these works use the ratios of reddened and
unreddened composite quasars to derive quasar extinction laws. In this
Letter we analyze their methodology to understand the discrepancy
between these works and that of Hopkins et al.

\section{Using composite quasar spectra to derive extinction laws}

The existence of large digital quasar surveys has led to the now
commonplace technique of combining spectra of many quasars into a
composite spectrum (Francis et al. 1991; Baker \& Hunstead 1995; Zheng
et al. 1997; Brotherton et al. 2001; Vanden Berk et al. 2001). These
composites improve the signal-to-noise ratio and wavelength baseline and
also allow ``typical'' quasar properties to be determined. However,
the range in {\em rest-frame} wavelengths as a function of redshift
probed by optical spectroscopy means there may be very little overlap
between the quasars that formed the rest-frame optical and UV parts
of the spectrum. Given that luminosity and redshift are always
correlated in flux-limited samples, this means that there is also a
luminosity-dependence of the composite as a function of
wavelength.

Czerny et al. (2004) used the five SDSS composites generated by
Richards et al. (2003) to determine the quasar extinction law.  These
composites were generated from quasars separated into five bins based
on the {\em relative} $g'-i'$ colors (the colors with respect to the
median as a function of redshift, see Richards et al. for more
details); composites 1 and 2 lie on the blue side of the $g'-i'$ distribution
and 3 and 4 lie on the red side. The quasars for composite 5 (the
``reddened composite'') were defined by a relative $g'-i'$ that was a
function of redshift to incorporate the redshift-dependent effect of
dust reddening on observed frame-colors.  Richards et al. considered
the asymmetric color distribution and the emission line properties of
composites 1-5 and concluded that the differences between 1-4 were due
to a combination of dust reddening and {\em intrinsic} spectral slopes
and that the intrinsic differences dominate. Czerny et al., however,
assume that the bluest quartile SDSS composite 1 represents an
unreddened quasar and composites 3-5 differ from 1 solely due to dust
reddening.  An extinction curve is generated by dividing composites
3-5 by composite 1 and averaging the result. For this reason alone, we
do not expect the Czerny et al. curve to truly represent dust
reddening.  As we will show, there is a further effect that casts
doubt on their results.

Gaskell et al. (2004) used composite quasars from Baker \& Hunstead
(1995) to determine a quasar extinction law. Baker \& Hunstead
generated composite quasar spectra for four groups of radio-loud
quasars from the Molonglo Quasar Survey (MQS). The groups were defined
by the ratio of the radio fluxes of the core and lobe, $\Re$. $\Re$ is
an orientation indicator and Baker \& Hunstead found that the
composite spectrum of low $\Re$ quasars is redder than the high $\Re$
composite, presumably due to dust reddening by a torus perpendicular
to the radio jet axis. Baker \& Hunstead also produced a composite
compact-steep-spectrum quasar and this also appeared highly reddened,
although in this case the explanation is more likely to be an
evolutionary effect. Gaskell et al. used a similar method to Czerny et
al. and compared the $\Re<0.1$ composite to the $\Re \geq 1$ composite
and the compact-steep-spectrum composite to the $0.1 \leq \Re < 1$
composite. The extinction curves from these two pairs agree very well.
Gaskell et al. also compared the $\Re \geq 1$ composite to the
UV-selected composite spectrum of the Large Bright Quasar Survey
(LBQS; Francis et al. 1991), which showed a similar extinction
curve. However, since the $\Re \geq 1$ composite is the least reddened
of all the radio-loud composites it is not clear that differences
between this and the LBQS composite can be attributed to reddening
rather than intrinsic spectral differences associated with the
different selection methods (radio-loud vs radio-quiet).

Both Czerny et al. and Gaskell et al.  found the result that the
extinction laws are flatter in the UV than the extinction laws of the
SMC, LMC or Milky Way.  Both extinction curves are plotted in
Fig.\,\ref{fig:gas} along with the curves for the SMC and Milky
Way. The curve of Gaskell et al. is flat at all wavelengths below
$0.4\,\mu$m. The Czerny et al. curve is similar to that of the SMC
down to $0.25\,\mu$m and below this it flattens to become much flatter
than that of either the SMC or the Milky Way at $<0.13\,\mu$m. These
flat UV extinction curves are particularly interesting since they
agree with theoretical predictions based on the grain size
distribution close to an AGN with the destruction of small grains
due to grain growth and sputtering in a high-density,
high-radiation environment (e.g. Maiolino et al. 2001a,b).

As discussed above, there are serious concerns about combining quasar
spectra that span a wide range of rest-frame wavelengths into a
composite. This is especially true when dealing with reddened quasars.
Richards et al. (2003) and White et al. (2003) showed that the
majority of the quasars with reddening $E_{B-V}$\,$\gtsimeq 0.2$ in their
surveys (which are optical and optical+radio selected, respectively)
are found at redshifts $z<1$. Although similarly reddened quasars
do exist at higher redshifts they are less likely to be found in
quasar samples with an optical magnitude limit. We outline here
several reasons why, in practice, quasar surveys contain more highly
reddened quasars at lower redshifts:

\begin{figure}
\resizebox{0.48\textwidth}{!}{\includegraphics{f1.eps}}
\caption{Wavelength dependence of extinction for various extinction
  laws: SMC (Pei 1992; {\it solid line}), Milky Way (MW; Pei 1992;
  {\it dotted line}), AGN extinction laws derived from composite
  quasar spectra (Czerny et al. 2004, {\it dash-triple-dotted line};
  Gaskell et al. 2004, {\it dash-dotted line}). The two dashed lines show the
  resulting extinction curves determined from a simulated sample of
  quasars reddened by the SMC law (Sec.\,3). It rapidly diverges from
  the SMC law in the UV due to the redshift-dependent extinction
  present in typical quasar samples.
\label{fig:gas}
}
\end{figure}

\begin{enumerate}

\item{At high redshift the selection is based on the rest-frame UV
  spectrum and $A_{UV}\,\gg$\,$A_V$\, for SMC or Milky Way extinction
  laws. Therefore the $A_V$\, necessary to exclude quasars from samples
  due to optical magnitude limits is a function of redshift.}

\item{There is mounting evidence of an increase in the fraction of
  obscured AGN at lower luminosities (e.g. Ueda et al. 2003; Simpson
  2005). Due to the optical magnitude limit, quasars at lower
  redshifts are typically less luminous. Therefore the ratio of the
  space densities of lightly reddened to unobscured quasars is
  greater at low redshifts than at high redshifts.}

\item{Low redshift, low luminosity quasars can have significant
  contamination of the red end of their spectra by light from the host
  galaxy. If this is not identified as host galaxy contamination it
  could lead to the quasar being classified as reddened. At high
  redshift and luminosities there is little contamination by the
  hosts. For example, the SDSS composite quasar of Vanden Berk et
  al. (2001) shows a spectral steepening longward of 0.5\,$\mu$m due to
  host galaxy contamination in low redshift quasars.}

\item{Red quasar selection is sometimes based on the observed frame
  color or spectral index. For an extinction law that induces
  spectral curvature (such as the SMC), the relationship between
  observed color or spectral index and $E_{B-V}$\, is also a function of
  redshift, such that high-redshift quasars have an $E_{B-V}$\, lower than
  that of low-redshift quasars with the same color. Richards et
  al. (2003) attempted to take account of this bias in their composite
  5 spectrum, however they did not for composites 3 and 4 which Czerny
  et al. also used to represent reddened quasars.}

\end{enumerate}
The consequence of all these points is that the quasars going into a
composite spectrum will have a negative correlation between $E_{B-V}$\, and
redshift.

\begin{figure}
\resizebox{0.48\textwidth}{!}{\includegraphics{f2.eps}}
\caption{Optical spectral index, $\alpha_{\rm opt}$, against redshift for the MQS
sample used to create the composites (Baker et al. 1999; {\it filled
circles}). Note that the reddest quasars in the observed-frame are
located preferentially at low redshifts. Also plotted are curves
showing how a power law in the observed-frame optical with
$\alpha_{\rm opt}$\,$=0.4$ is affected by SMC ({\it solid line}) or Milky Way
({\it dotted line}) dust reddening as a function of redshift. This
shows that for a given $\alpha_{\rm opt}$, the $A_V$\, is greater for low redshift
quasars. Combining these two effects means that the $A_V$\, (or $E_{B-V}$\,) of
the quasars is a strong function of redshift in this sample. For the
Gaskell et al. extinction law ({\it dash-dotted line}), unfeasibly high values
of extinction are necessary to redden the high-redshift quasars.
\label{fig:mqs}
}
\end{figure}

To illustrate the significance of these effects in the light of the
current discussion, in Fig.\,\ref{fig:mqs} we plot optical spectral
index, $\alpha_{\rm opt}$, against redshift for the MQS sample of Baker et
al. (1999) which was used to generate the composites in Baker \&
Hunstead (1995). It can be seen that in this sample, most of the
reddest quasars ($\alpha_{\rm opt}$$>1.5$) are at lower redshifts. This is a
consequence of points 1--3 outlined above. However, there is not a
significant correlation between redshift and $\alpha_{\rm opt}$\ in the full
sample because the bluest quasars are also predominantly at low
redshift.

The fact that there is no strong positive or negative correlation of
$\alpha_{\rm opt}$\ with redshift does not mean however that there will be no
bias when making a composite from these quasars, because of point 4.
On Fig.\,\ref{fig:mqs} we also plot curves showing how the observed
spectral index changes with reddening as a function of redshift. For
simplicity we represent a quasar by a single power law with
$\alpha_{\rm opt}$\,$=0.4$ and determine the observed $\alpha_{\rm opt}$\ by the ratio of
fluxes in observed-frame wavelength bins $0.8-0.85\,\mu$m and
$0.35-0.4\,\mu$m (roughly corresponding to the extreme useful range of
optical spectroscopy). The power law is reddened by $A_V$\,$=0.5,1$ and $1.5$
adopting the SMC and Milky Way (excluding the $0.22\,\mu$m feature not
usually seen in quasars) extinction laws (Pei 1992) and $A_V$\,$=2,4$ and $6$
adopting the law of Gaskell et al. (2004). For the SMC and Milky Way
laws there is a strong correlation of $A_V$\, with redshift in this
sample.  For the SMC dust, almost all the quasars with
$A_V>0.5$ are at $z<1$. Even for the Milky Way law which does not
induce such strong curvature, a correlation is seen at $z>2$.

This plot also illustrates a fundamental problem with the extinction
law of Gaskell et al. (2004). The extreme flatness in the UV means
that a huge $A_V$\, is necessary to get any reddening at all of the UV
continuum. Many surveys, including the MQS and the SDSS, contain
reddened quasars at $z>2$ where only the UV is visible in optical
spectroscopy. According to the extinction law of Gaskell et al. the
quasars with $\alpha_{\rm opt}$$\approx 1.5$ in the MQS must be reddened by
$A_V$\,$\sim 20$ and therefore have intrinsic absolute magnitudes of
$M_B\sim -45$! What is actually found for high-redshift quasars in the
SDSS that have rest-frame optical data available is that the reddening
in the optical is much lower than in the UV, similar to SMC dust
(Hopkins et al. 2004).

\section{Simulating the biases of composite reddened quasars}

It is clear that combining the quasars in Fig.\,\ref{fig:mqs} into
composites could lead to serious problems for deriving a dust
extinction law. In particular, there is an inconsistency in the
extinction law of Gaskell et al. (derived from the MQS) and the
existence of red high-redshift quasars in the MQS. To understand this
quantitatively we now consider what would happen if one were to
generate a composite from a quasar sample selected on the basis of
their red observed optical colors (or equivalently, steep
$\alpha_{\rm opt}$). This analysis only includes the bias due to point 4 of the
previous section.

We consider a hypothetical sample of 211 quasars distributed smoothly
over the redshift range $0.4 \leq z \leq 2.5$. This redshift range is
typical of quasar samples such as the MQS and the SDSS study of
Richards et al. (2003). It is also the redshift range where the
observed-frame wavelengths of optical spectroscopy (assumed
$0.38-0.92\,\mu$m as in the SDSS) are covered by our composite
spectrum ($0.1-0.8\,\mu$m). We assume that the intrinsic spectrum of
each quasar is represented by a composite spectrum specifically
constructed from the SDSS to represent a typical unreddened quasar
(Willott et al. in prep.).

Each quasar is reddened by some amount of SMC type dust until its
observed $g'-i'$ color equals 0.7. We also repeated the analysis
using relative $g'-i'$ colors equal to 0.5. Note that the use of
relative colors takes out the effect of emission lines and is a
closer approximation to selecting quasars with a constant $\alpha_{\rm opt}$. As
we will show, there is little difference in the results when using
$g'-i'$ or relative $g'-i'$. To get an idea of the amounts of
reddening required to produce these colors as a function of redshift,
$A_V$\,$=0.8$ at $z=0.4$ and $A_V$\,$=0.19$ at $z=2.5$ for relative
$g'-i'=0.5$. For reference to the MQS data in Fig.\,\ref{fig:mqs},
this corresponds to observed $\alpha_{\rm opt}$\,$\approx\,1.5$. We repeated our
analysis with different values for the colors and found only minor
differences in the results.

The individual reddened quasars are then combined using the same
procedure as in Vanden Berk et al. (2001). Because we use the
geometric mean composite spectrum, the spectral shape, extinction law
and mean extinction should all be preserved in the absence of any
biases (Reichard et al. 2003).  Therefore we now divide the composite
reddened quasar spectrum by the intrinsic quasar spectrum to determine
the extinction curve in the same manner as Czerny et al. and Gaskell
et al.

The results are shown in Fig.\,\ref{fig:gas} where we plot the curves
generated from our two composites (using constant $g'-i'$ and constant
relative $g'-i'$ reddening). The first thing to notice is that our
derived composite extinction curves are nothing like the original SMC
curve which was used to redden the individual quasars.  The typically
lower reddening in the high redshift quasars that dominates the UV
part of the spectrum lead to a much flatter extinction curve than the
SMC. In particular, at $\lambda^{-1}> 6\,\mu{\rm m}^{-1}$, the curves
are much flatter than either the SMC or the Milky Way, which is a
characteristic of extinction curves generated from composite quasars.

Our simulated extinction curves are even flatter than that of Czerny
et al., but have a shape remarkably similar to theirs. The fact that
one of three reddened composites used by Czerny et al. was
selected taking account of this redshift dependent bias may account
for the fact that their curve is somewhat closer to the SMC curve than
our derived curves. Therefore we find that the difference between the
extinction laws of Czerny et al. and the SMC can be fully accounted
for by this bias.

The Gaskell et al. law is extremely flat in the UV, so if their curve
really does come from quasars extincted with SMC type dust then a
stronger bias than represented by an approximately constant $\alpha_{\rm opt}$\,
must exist. We can only speculate on exactly why this is without full
rest-frame UV and optical spectroscopy of the MQS quasars that went into
the composites. We note that Baker \& Hunstead (1995) scaled their
quasar spectra before combining to a common flux-density at
0.3\,$\mu$m. If we repeat our analysis using this scaling, rather than
the method of ordering each quasar by its redshift and scaling to the
common wavelength region of the previous mean (Vanden Berk et
al. 2001), then we find somewhat flatter curves in the UV than shown
in Fig.\,\ref{fig:gas}. Also, we note that the number of quasars in
each MQS composite was small (in the range 13-18) and at the extreme
wavelength ranges of the composites even fewer quasars were
contributing. Therefore object-to-object differences could have had a
large impact. Gaskell et al. claimed that the fact that the broad
emission lines revealed the same extinction curve as the continuum was
independent evidence of the flat extinction curve. The broad
emission lines are subject to exactly the same bias as the continuum,
so it is not in fact independent.

\section{Discussion}

Our analysis shows that the gray UV extinction laws derived using
composite quasars are an erroneous artifact of the technique. This is
due to the fact that the quasars contributing to the composite in the
UV have typically lower reddening than those that contribute in the
optical.

There have been other studies that have claimed that gray dust may be
present in AGN environments. Maiolino et al. (2001a,b) presented a
number of arguments for a different grain size distribution, devoid of
small grains. However their arguments are based largely on
mid-infrared extinction (which is dominated by very large grains) and
the different observed gas-to-dust absorption ratios in Seyfert
nuclei. However, Willott et al. (2004) showed that this difference in
gas-to-dust absorption ratios is also present in X-ray selected
quasars whose rest-frame UV and optical reddening can be fit with SMC
or Milky Way extinction laws, not gray dust.

We argue that evidence based on rest-frame UV and optical spectral
energy distributions for individual reddened quasars (as in Hopkins et
al. 2004) is much more reliable than extinction laws based on
composite spectra. We conclude that the SMC curve
remains a useful description of the reddening in quasar environments
and that there is little evidence supporting extinction curves that
flatten in the UV.

\acknowledgments Thanks to Jo Baker, Daniel Fryer, Matt Jarvis, Ross
McLure and the anonymous referee for assistance and useful discussions.


The Case of the 300 kpc Long X-ray Jet in PKS 1127-145 at z=1.18

Abstract: The complex X-ray morphology of the 300 kpc long X-ray jet in PKS1127-145 (z=1.18 quasar) is clearly discerned in a ~100 ksec Chandra observation. The jet X-ray surface brightness gradually decreases by an order of magnitude going out from the core. The X-ray spectrum of the inner jet is relatively flat with alpha_X=0.66+/-0.15 and steep in the outer jet with alpha_X=1.0+/-0.2. The X-ray and radio jet intensity profiles are strikingly different, with the radio emission peaking strongly at the two outer knots while the X-ray emission is strongest in the inner jet region. We discuss the constraints implied by these data on the X-ray emission models and conclude that ``one-zone'' models fail and that at least a two-component model is needed to explain the jet's broadband emission. We propose that the X-ray emission originates in the jet proper while the bulk of the radio emission comes from a surrounding jet sheath. We also consider intermittent jet activity as a possible cause of the observed jet morphology.

\section{Introduction}

The origin of the X-ray emission in quasar jets is puzzling
\citep[see][for review]{harris2006} because a straight
extrapolation of the synchrotron radio-to-optical continuum into the
X-rays underpredicts the luminosity of powerful {\it Chandra} large
scale jets. Thus, the same single power-law population of electrons
cannot produce the radio, optical, and X-ray emission in the framework
of a homogeneous one-emission zone approximation \citep[e.g.][]{sambruna2004}
The synchrotron self-Compton (SSC) process cannot easily
explain the data because it does not produce enough X-rays at the
equipartition fields and therefore requires large departures from the
minimum-power condition \citep{chartas2000,harris2002,kataoka2005}.
The alternative scenario involves
inverse-Compton scattering of the Cosmic Microwave Background (IC/CMB)
photons implying large jet bulk Lorentz factors ($\Gamma$) at hundreds
of kpc from the active nuclei \citep{tavecchio2000,celotti2001,schwartz2002}.
In this paper we discuss the constraints on
the X-ray emission models given by the {\it Chandra} and VLA
observations of the PKS~1127-145 jet.

\section{Jet Properties}

\begin{figure}[t]
\begin{center}
\includegraphics[width=7.8cm]{img_grey_0.3-8_g2555.ps}
\caption{{\it Chandra} ACIS-S exposure corrected image of PKS1127--145
(E=0.3-7~keV) and smoothed with the Gaussian kernel
(1$\sigma$=0.615~arcsec).   Contours are factor of two in log scale starting at
10$^{-10}$ photons~cm$^{-2}$~s$^{-1}$~pixel$^2$. The two final contours
within the core are at 10$^-6$ and 10$^-7$
photons~cm$^{-2}$~s$^{-1}$~pixel$^2$ (1 pixel=0.123'').}
\end{center}
\label{xrays}
\end{figure}

One of the longest X-ray jets known (see
\texttt{http://hea-www.harvard.edu/XJET/}) is associated with the
redshift z=1.18 radio-loud quasar PKS~1127-145 (Fig.~1). This jet
poses several questions for X-ray emission models
\citep[see][]{siemi2002,siemi2007}. The jet morphology is
complex with several knots seen in the X-ray image (Fig.~1) with the
furthest knot C being the weakest. The prominent knots A and B are
connected by continuous X-ray emission that stops beyond knot B at
$\sim 22$~arcsec from the core. The VLA radio maps show low brightness
emission along the jet (Fig.~2). The one-sided jet shows an X-ray
surface brightness which declines with the distance from the core,
while the radio brightness increases (Fig.~3[left]).  The X-ray to
radio intensity ratio decreases along the jet. Fig.3[right] shows the
changes to spectral index along the jet. The radio spectrum is steep
and X-ray spectrum is flat in the inner jet regions.  There are only
upper limits to the knots optical brightness, but they are too high
for constraining the emission models.  Thus only radio and X-ray
observations can be used to study emission processes in this jet.  We
discuss the constraints on the models below.

\begin{figure}
\begin{center}
\includegraphics[width=7.6cm]{f6.eps}
\caption{VLA 1.4 and 4.9 GHz images of PKS~1127--145. The lowest contour levels
plotted are  0.75 mJy beam$^{-1}$ (1.4 GHz) and
0.39 mJy beam$^{-1}$ (4.9 GHz) and increase by
factors of $\sqrt{2}$ up to peaks of 5.46 and 4.14 Jy/beam,
respectively. The (uniform weighted) beam sizes are:
2.71\arcsec$\times$1.62\arcsec\ at PA=84.4$^{\circ}$,
and 2.77\arcsec$\times$1.66\arcsec\ at PA=84.0$^{\circ}$, respectively.}
\end{center}
\label{radio}
\end{figure}

\begin{figure}[t]
\includegraphics[width=7cm]{f9.eps}
\includegraphics[width=7cm]{f10.eps}
\caption{{\bf Left:} Flux along the jet for each observed frequency:
long dashed line -- hard X-rays (2-10~keV); short dashed line - soft
X-rays (0.5-2~keV); dotted line - R1 = 1.4 GHz; short-long-dashed - R2
= 5~GHz; and dash-dot - R3 = 8.5~GHz.  Lower panel shows the flux
ratio (1keV/1.4GHz) along the jet. {\bf Right:} Spectral index
($S_{\nu} \sim \nu^{-\alpha}$) along the jet: $\alpha_X$- X-ray
squares, and
long-short dash; $\alpha_{R3-X}$ - 8.5~GHz to 1~keV
circles,
dot-dashed line; $\alpha_{R2-X}$ - 5~GHz to 1~keV,
circles,
dashed line; $\alpha_{R1-X}$ - 1.4~GHz to 1keV, circles,
solid line; $\alpha_{R1-R2}$ - 1.4~GHz to 5~GHz, triangles,
dotted line, $\alpha_{R2-R3}$ - 5~GHz to 8.5~GHz, triangles
dashed line.  }
\label{indices}
\end{figure}

\section{Jet Emission Models}

A detailed discussion of the model constraints can be found in
\citet{siemi2007}. Here we give a summary of the main results:

\subsection{One Zone Model Failures}

\subsubsection{Synchrotron} X-ray emission requires a continuous
acceleration of high energy particles along the jet.  The expected
distribution of the freshly accelerated electrons gives $\alpha_X=
0.5$ as observed in the inner jet, but the radio spectrum is steep
with $\alpha_R > 1.0$ (Fig.3). A different cooling mechanism or
difference in electron populations emitting radio and X-rays is
needed.

\subsubsection{IC/CMB} models in the simplest version involve
an adiabatically expanding jet. In this case the observed
inverse-Compton (X-rays) and synchrotron (radio) intensity ratio is
different than the predicted one.  In the case of an efficient jet
confinement with some deceleration and B-field amplification it is
possible to get an agreement between the predicted and the observed
luminosity ratios, but the required B-field amplification is high
[$B(r_C)/B(r_o)$]$^2 \sim 50$. In this case for the equipartition at
$r_o$ there will be no equipartition at $r_C$.  The expected radio
spectrum should steepen more significantly than the X-ray
spectrum. The observed behavior is different: $\alpha_R$ decreases,
while $\alpha_X$ increases (Fig.3[right]).

\subsection{Two-Component Model: The Jet and a Sheath}

Radio and X-ray emissions are produced in two separate regions. The
X-ray emission comes from the proper relativistically moving jet and
the radio emission from the sheath - a slow moving radial extension of
the jet boundary layer.  The X-ray emission is due to synchrotron or
IC/CMB process.

\subsubsection{Proper Jet:}

The observed smooth decrease of the X-ray flux, when unconnected with
the increase of the radio one, is expected as a result of the decrease
in the inverse-Compton or synchrotron luminosities along decelerating
and/or expanding jet.

\subsubsection{Sheath:}

The sheath dominates the jet's radio emission and has two cooling
regions: (1) radiative cooling at the inner sheath: $t_{rad}<10$~Myr$<
t_{dyn}= 30 R_{10}$~Myr.  Due to the frequency dependent losses
a steep radio spectrum is expected in this regime.  (2) Adiabatic
cooling in the outer sheath: $t_{rad} > t_{dyn} < 10R_{10}$Myr. Because
of the frequency independent losses a flat spectrum is expected.

\subsection{Modulated Jet Activity?}

Can the morphology of the $>$300~kpc long jet be a result of the
intermittent activity of the quasar?  The knots are potentially too
large ($>$10~kpc) to be considered a result of extended shock waves
formed within a continuous jet outflow. They might form during
episodes of separate continuous activity. Then the size gives a
duration for each epoch of $\sim 10^5$ years.  The radio core can be
related to blazar phenomena \citep{blazej2004}, but the core also
resembles a compact GPS object at this time.

\section{Summary}

The key results of this project:

\begin{itemize}

\item The jet is long: radio and X-ray jet emission is detected up to $\sim$300~kpc
(projected distance) from the core.

\item The jet X-ray brightness is decreasing while the radio
brightness is increasing with distance from the quasar.

\item The one component X-ray/radio emission models failed to explain the
observations.

\item X-ray and radio jet properties suggests two separate
components for the jet emission.

\item Intermittent quasar activity might be
reflected in the observed jet morphology: each knot represents a
continuous jet activity with a duration of 10$^5$years.

\end{itemize}

\acknowledgements

This research is funded in part by NASA contract NAS8-39073 and
Chandra Awards GO5-6113X, GO4-5131X, GO5-6118X, GO5-6111X. The
National Radio Astronomy Observatory is operated by Associated
Universities, Inc. under an agreement with the National Science
Foundation.


Angular Diameters of Compact Planetary Nebulae

Abstract: We have obtained H$\alpha$ CCD images of 21 compact planetary nebulae in the galactic bulge. All objects are resolved and, after deconvolving the seeing disk, we find diameters of 1--2$''$ with typical uncertainties of 20\%. The values are generally in good agreement with radio measurements but are smaller than older optical measurements. Based on the new values, we find that the Shklovsky method for distance determination greatly overestimates the actual distance for all of our objects.

\section*{References}\list
 {}{\setlength\labelwidth{1.4em}\leftmargin\labelwidth
 \setlength\parsep{0pt}\setlength\itemsep{0pt}
 \setlength{\itemindent}{-\leftmargin}
 \usecounter{enumi}}
 \def\hskip .11em plus .33em minus -.07em{\hskip .11em plus .33em minus -.07em}
 \sloppy
 \sfcode`\.=1000\relax}
\let\endthebibliography=\endlist

\typeout{Document Style `l-aa' version 3.0}

\frenchspacing
\def~{\penalty\@M\kern3pt}

\def\@ptsize{0} \@namedef{ds@11pt}{\def\@ptsize{1}}
\@namedef{ds@12pt}{\def\@ptsize{2}}
\def\ds@twoside{\@twosidetrue \@mparswitchtrue}
\def\ds@draft{\overfullrule 5\p@}
\def\ds@noframe{\let\picplace\vspace}
\newif\if@referee \@refereefalse
\def\ds@referee{\@refereetrue}

\def\picplace#1{\vbox{\hrule\@height 0.4pt\@width\hsize
\hbox to\hsize{\vrule\@width 0.4pt\@height#1\hfil
\vrule\@width 0.4pt\@height#1}\hrule\@height 0.4pt\@width\hsize}}

\@options
\ds@twoside

\widowpenalty=10000
\clubpenalty=10000
\def\thisbottomragged{\def\@textbottom{\vskip\phantom{0}@ plus.0001fil
\global\let\@textbottom\relax}}

\lineskip 1pt \normallineskip 1pt
\def1{1}

\def\@normalsize{\@setsize\normalsize{12pt}\xpt\@xpt
\abovedisplayskip=3 mm plus6pt minus 4pt
\belowdisplayskip \abovedisplayskip
\abovedisplayshortskip=0mm plus6pt
\belowdisplayshortskip=2 mm plus4pt minus 4pt
\def\@listi{\topsep 4pt plus 2pt minus 2pt
\leftmargin \leftmargini
\parsep 0pt
\itemsep \parsep}}
\def\small{\@setsize\small{11pt}\ixpt\@ixpt
\abovedisplayskip 8.5pt plus 3pt minus 4pt
\belowdisplayskip \abovedisplayskip
\abovedisplayshortskip \phantom{0}@ plus 2pt
\belowdisplayshortskip 4pt plus2pt minus 2pt
\def\@listi{\topsep 3pt plus 2pt minus 2pt
\leftmargin \leftmargini
\parsep 0pt
\itemsep \parsep}}
\let\footnotesize=\small
\def\@setsize\scriptsize{8pt}\viipt\@viipt{\@setsize\@setsize\scriptsize{8pt}\viipt\@viipt{8pt}\viipt\@viipt}
\def\@setsize\tiny{6pt}\vpt\@vpt{\@setsize\@setsize\tiny{6pt}\vpt\@vpt{6pt}\vpt\@vpt}
\def\@setsize\large{17pt}\xivpt\@xivpt{\@setsize\@setsize\large{17pt}\xivpt\@xivpt{17pt}\xivpt\@xivpt}
\def\@setsize\Large{17pt}\xivpt\@xivpt{\@setsize\@setsize\Large{17pt}\xivpt\@xivpt{17pt}\xivpt\@xivpt}
\def\@setsize\LARGE{20pt}\xviipt\@xviipt{\@setsize\@setsize\LARGE{20pt}\xviipt\@xviipt{20pt}\xviipt\@xviipt}
\def\@setsize\huge{25pt}\xviipt\@xviipt{\@setsize\@setsize\huge{25pt}\xviipt\@xviipt{25pt}\xviipt\@xviipt}
\def\@setsize\Huge{30pt}\xviipt\@xviipt{\@setsize\@setsize\Huge{30pt}\xviipt\@xviipt{30pt}\xviipt\@xviipt}
\@normalsize\rm

\parskip 0pt plus 1pt
\topsep 10pt plus2pt minus5pt
\partopsep 0pt
\itemsep 0pt plus 1pt
\parindent 1.5em
\leftmargin    \parindent
\leftmargini   \leftmargin
\leftmarginii  \leftmargin
\leftmarginiii \leftmargin
\leftmarginiv  \leftmargin
\leftmarginv   \leftmargin
\leftmarginvi  \leftmargin
\labelwidth\leftmargini
\labelsep 4pt
\advance\labelwidth-\labelsep
\def\@listii{\leftmargin\leftmarginii
 \labelwidth\leftmarginii\advance\labelwidth-\labelsep
 \topsep \phantom{0}@}
\def\@listiii{\leftmargin\leftmarginiii
 \labelwidth\leftmarginiii\advance\labelwidth-\labelsep
 \topsep \phantom{0}@}
\def\@listiv{\leftmargin\leftmarginiv
 \labelwidth\leftmarginiv\advance\labelwidth-\labelsep}
\def\@listv{\leftmargin\leftmarginv
 \labelwidth\leftmarginv\advance\labelwidth-\labelsep}
\def\@listvi{\leftmargin\leftmarginvi
 \labelwidth\leftmarginvi\advance\labelwidth-\labelsep}

\oddsidemargin 63pt \evensidemargin 63pt \marginparwidth 90pt
\marginparsep 11pt
\voffset -3mm
\hoffset -25mm

\topmargin 27pt \headheight 12pt \headsep 14pt \footheight 12pt
\footskip 30pt
\textheight=56pc
\textwidth=180mm
\columnsep=4mm
\columnseprule=0pt
\hfuzz=2pt
\vfuzz=2pt
\tolerance=1000

\footnotesep 8pt
\skip\footins 9pt plus 4pt minus 2pt
\floatsep 12pt plus 2pt minus 2pt
\textfloatsep 20pt plus 2pt minus 4pt
\intextsep 20pt plus 2pt minus 2pt
\@maxsep 20pt
\dblfloatsep 12pt plus 2pt minus 2pt
\dbltextfloatsep 20pt plus 2pt minus 4pt
\@dblmaxsep 20pt
\@fptop 0pt plus 1fil
\@fpsep 8pt plus 2fil
\@fpbot 0pt plus 1fil
\@dblfptop 0pt plus 1fil
\@dblfpsep 8pt plus 2fil
\@dblfpbot 0pt plus 1fil
\marginparpush 5pt

\@lowpenalty 51 \@medpenalty 151 \@highpenalty 301 \@beginparpenalty
-\@lowpenalty \@endparpenalty -\@lowpenalty \@itempenalty -\@lowpenalty

\def\vec#1{\ifmmode\mathchoice{\mbox{\boldmath$\displaystyle#1$}}
{\mbox{\boldmath$\textstyle#1$}}
{\mbox{\boldmath$\scriptstyle#1$}}
{\mbox{\boldmath$\scriptscriptstyle#1$}}\else
\hbox{\boldmath$\textstyle#1$}\fi}

\def\tens#1{\ifmmode\mathchoice{\mbox{$\sf\displaystyle#1$}}
{\mbox{$\sf\textstyle#1$}}
{\mbox{$\sf\scriptstyle#1$}}
{\mbox{$\sf\scriptscriptstyle#1$}}\else
\hbox{$\sf\textstyle#1$}\fi}

\def\par \addvspace{4ex} \@afterindentfalse \secdef\@part\@spart{\par \addvspace{4ex} \@afterindentfalse \secdef\@part\@spart}
\def\@part[#1]#2{\ifnum \c@secnumdepth >\m@ne \refstepcounter{part}
\addcontentsline{toc}{part}{\Roman{part}} \def\thesection {\arabic{section} \hspace{1em}#1}\else
\addcontentsline{toc}{part}{#1}\fi { \parindent 0pt \raggedright \ifnum
\c@secnumdepth >\m@ne \@setsize\Large{17pt}\xivpt\@xivpt \bf Part \Roman{part}} \def\thesection {\arabic{section} \par \nobreak \fi \@setsize\huge{25pt}\xviipt\@xviipt
\bf #2\par } \nobreak \vskip 3ex \@afterheading }
\def\@spart#1{{\parindent 0pt \raggedright \@setsize\huge{25pt}\xviipt\@xviipt \bf #1\par} \nobreak
\vskip 3ex \@afterheading }

\def\section{\@startsection {section}{1}{\phantom{0}@}{-3.5ex plus -1ex minus
 -.2ex}{1.5ex plus .2ex}{\normalsize\bf\boldmath}}
\def\subsection{\@startsection{subsection}{2}{\phantom{0}@}{-3.25ex plus -1ex
minus -.2ex}{1.5ex plus .2ex}{\normalsize\it}}
\def\subsubsection{\@startsection{subsubsection}{3}{\phantom{0}@}{-3.25ex plus
-1ex minus -.2ex}{1.5ex plus .2ex}{\normalsize}}
\def\paragraph{\@startsection {paragraph}{4}{\phantom{0}@}{3.25ex plus 1ex minus
.2ex}{-.6em}{\normalsize\it}}
\def\subparagraph#1{\typeout{AandA Warning: You should not use
\noexpand\protect\noexpand\subparagraph \space in this style.}\vskip0.5cm
You should not use $\backslash${\tt subparagraph} in this
style.\vskip0.5cm}

\def\@sect#1#2#3#4#5#6[#7]#8{\ifnum #2>\c@secnumdepth
     \def\@svsec{}\else
     \refstepcounter{#1}\edef\@svsec{\csname the#1\endcsname.\ }\fi
     \@tempskipa #5\relax
      \ifdim \@tempskipa>\phantom{0}@
 \begingroup #6\relax
   \@hangfrom{\hskip #3\relax\@svsec}{\interlinepenalty \@M #8\par}
 \endgroup
       \csname #1mark\endcsname{#7}\addcontentsline
  {toc}{#1}{\ifnum #2>\c@secnumdepth \else
        \noexpand\protect\noexpand\numberline{\csname the#1\endcsname.}\fi
      #7}\else
 \def\@svsechd{#6\hskip #3\@svsec #8\csname #1mark\endcsname
        {#7}\addcontentsline
      {toc}{#1}{\ifnum #2>\c@secnumdepth \else
        \noexpand\protect\noexpand\numberline{\csname the#1\endcsname.}\fi
         #7}}\fi
     \@xsect{#5}}

\def\@xsect#1{\@tempskipa #1\relax
      \ifdim \@tempskipa>\phantom{0}@
       \par \nobreak
       \addvspace{\@tempskipa}
       \@afterheading
    \else \global\@nobreakfalse \global\@noskipsectrue
       \everypar{\if@noskipsec \global\@noskipsecfalse
     \clubpenalty\@M \hskip -\parindent
     \begingroup \@svsechd \endgroup \unskip
     \hskip -#1
    \else       \clubpenalty \@clubpenalty
      \everypar{}\fi}\fi\ignorespaces}

\setcounter{secnumdepth}{3}

\def