MOJAVE: Monitoring of Jets in AGN with VLBA Experiments - II. First-Epoch 15 GHz Circular Polarization Results

Abstract: We report first-epoch circular polarization results for 133 active galactic nuclei in the MOJAVE program to monitor the structure and polarization of a flux limited sample of extra-galactic radio jets with the VLBA at 15 GHz. We found strong circular polarization ($\geq 0.3$%) in approximately 15% of our sample. The circular polarization was usually associated with jet cores; however, we did find a few strong jet components to be circularly polarized. The levels of circular polarization were typically in the range of $0.3-0.5$% of the local Stokes-$I$. We found no strong correlations between fractional circular polarization of jet cores and source type, redshift, EGRET detections, linear polarization, or other observed parsec-scale jet properties. There were differences between the circular-to-linear polarization ratios of two nearby galaxies versus more distant quasars and BL Lac objects. We suggest this is because the more distant sources either have (1) less depolarization of their linear polarization, and/or (2) poorer effective linear resolution and therefore their VLBA cores apparently contain a larger amount of linearly polarized jet emission. The jet of 3C 84 shows a complex circular polarization structure, similar to observations by Homan & Wardle five years earlier; however, much of the circular polarization seems to have moved, consistent with a proper motion of 0.06$c$. The jet of 3C 273 also has several circularly polarized components, and we find their fractional circular polarization decreases with distance from the core.

\section{Introduction}
\label{s:intro}

Only a small fraction of the radio emission from
extra-galactic radio jets associated with Active Galactic Nuclei (AGN)
is circularly polarized; however, this weak component is a potentially
important probe of the particle population
and magnetic field structure of the jets.  Circular polarization (CP)
can be produced either as an intrinsic component of the emitted
synchrotron radiation, or via Faraday conversion
of linear polarization into circular \citep{JOD77}.  If the intrinsic
mechanism dominates, jets must be a predominantly electron-proton
plasma and must contain a significant fraction of unidirectional
ordered magnetic field.  If the conversion mechanism dominates,
circular polarization probes the low energy end of the relativistic
particle distribution \citep[e.g.][]{WH03}, a key parameter in
studying the bulk kinetic luminosity of AGN jets and
their particle content \citep{CF93, WHOR98}.
Faraday conversion is expected to be the dominant mechanism on
theoretical grounds \citep[e.g.][]{J88}; however, there is still
only limited observational evidence supporting this view
\citep{WHOR98, HW04}.

Early integrated measurements of CP from AGN jets \citep{WdP83, K84}
revealed circular polarization to be only a very small fraction,
typically $< 0.1$\%, of their total flux.  However, the Very Long Baseline
Array (VLBA)\footnote{The VLBA is operated by the National Radio Astronomy
Observatory, a facility of the National Science Foundation operated under
cooperative agreement by Associated Universities, Inc.}, has made it possible to
study circular polarization from extragalactic jets at sub-milliarcsecond
resolution.  \citet{WHOR98} and \citet[hereafter HW99]{HW99} found circular
polarization in cores of four powerful AGN jets at {\em local} fractional levels of
$0.3$ to $1$\% of Stokes-$I$ with the VLBA.  Circular polarization has since
been detected with various instruments in other AGN jets \citep{RNS00, HAW01}
and a wide variety of other synchrotron emitting sources:  low-luminosity AGN such
as Sagittarius A* and M81* \citep*{BFB99,SM99,BBFM01,BFM02}, intra-day
variable sources \citep{MKRJ00}, and galactic micro-quasars
\citep{Fend00,Fend02}.

Despite the detection of the circular polarization in a variety of sources,
observations of circular polarization from AGN jets are fragmented and incomplete.
A large temporal gap exists between
the integrated measurements of the 1970's to early 1980's and the high
resolution VLBA measurements from the mid 1990s onward.  Only a few
AGN jets have high resolution, multi-epoch circular polarization
observations that are matched in observing frequency and sensitivity.
An even smaller number of sources have multi-frequency circular polarization
observations suitable for measuring a spectrum.  Additionally, the
small, {\em ad hoc} samples studied to date have not allowed for an
analysis of correlations between the appearance of strong circular
polarization and other source properties.  While these previous studies
raised a number of interesting possibilities of how circular polarization
could be used to constrain the magnetic field properties and
particle populations of AGN jets, a systematic survey of a large
complete sample is necessary to make further progress.

We have exploited the VLBA's unique imaging capabilities by studying the
milliarcsecond-scale circular polarization of a complete AGN jet
sample at 15 GHz as part of the MOJAVE (Monitoring Of Jets in AGN with
VLBA Experiments) program.
The MOJAVE sample
of 133 AGN is complete and flux limited on the
basis of VLBA-scale emission.  The selection criteria are described
in detail in \citet{LH05} which also reports first-epoch linear
polarization results; however, in brief, for inclusion in MOJAVE,
a source must have had a total 15 GHz VLBA flux density of at least
1.5 Jy ($> 2$ Jy for sources south of the celestial equator) at any
epoch from 1994--2003.  When completed, MOJAVE will have obtained
at least four epochs of full Stokes VLBA images at 15 GHz of all
sources in the sample.  In addition, multi-frequency follow-up
programs will study the circular polarization spectrum of the strongly
circularly polarized sources identified by MOJAVE.

In this paper, we report only the first-epoch circular polarization
results from the MOJAVE program.  Future papers will explore other
aspects of this rich data set including multi-epoch circular polarization
variability.  The details of our observations and calibration appear
in section 2.  We present and discuss our results in section 3, and
our conclusions appear in \S{4}.

\section{Observations and Calibration}
\label{s:obs}
Each MOJAVE observing session targets
18 sources over a 24 hour period, and this paper presents results drawn from
the first eleven sessions, from 2002 May 31 through 2004 February 11.
During each 24 hour session, observations of the 18 sources are broken up
into highly interleaved scans, each approximately seven minutes in duration.
The schedule is designed to maximize coverage in the $(u,v)$-plane; however, as
discussed below, the highly interleaved nature of the observations also improves our
ability to calibrate the antenna gains to detect circular polarization.  We
observed each source for approximately 65 minutes of total integration time. No additional
calibrators were observed as all of our target sources were strong and
compact.

The data were recorded at each VLBA antenna at a central observing frequency
of 15.366 GHz using 1-bit sampling and were correlated at the VLBA correlator
in Socorro, NM.  The correlator output was 2 s integrations for all four
cross-correlations (RR, RL, LR, LL).  The correlated data contained four
intermediate frequencies (IFs) and 16 channels per IF, yielding a total bandpass
of 32 MHz in both RR and LL.  Tapes of the correlated data were shipped to NRAO
in Charlottesville,
VA and Purdue University where they were loaded into NRAO's Astronomical Image
Processing System (AIPS) for further processing.

The basic calibration procedure followed the standard methods outlined
in the AIPS cookbook and the procedure by HW99 for measuring
circular polarization with the VLBA using the {\em gain transfer} technique.
Below we briefly describe our basic calibration procedure, in
\S{\ref{s:gains}} we discuss some improvements to the {\em gain transfer}
procedure developed for these observation, and in \S{\ref{s:phase}}
we discuss some checks on our calibration procedure.

The data were loaded into AIPS using the VLBALOAD procedure.  The procedure
VLBACALA was then used to do a-priori amplitude calibration, running both the
ACCOR and APCAL tasks, including an atmospheric opacity correction.  The parallactic angles were
removed using the procedure VLBAPANG and the pulse calibration was applied
using VLBAPCOR.  The procedure VLBAFRNG was then used to perform fringe fitting.
All of our sources were strong enough to serve as their own fringe calibrator, and
no special interpolation from strong to weak sources was necessary.
We used the procedure CRSFRING to remove any residual multi-band delay
difference between the right and left-hand systems of the array.  Finally we ran
BPASS to apply a bandpass correction before averaging across the channels within
each IF and 'splitting' the data off into individual source files.

Following a 2 s time-scale point-source phase calibration to increase
coherence, the data were written out to the CALTECH VLBI program, DIFMAP
\citep{S97}.  In DIFMAP, the data were averaged to 30 seconds and edited
in a station-based manner.  An edited version of the data was then saved without
any self-calibration applied.  We will refer to this version of the data as the
edited/unself-calibrated data.  Normal self-calibration and imaging techniques
were then used within DIFMAP to obtain the best Stokes-$I$ model for
each source.

For each source, the best Stokes-$I$ model and the edited/unself-calibrated data
were then read back into AIPS, and the data were self-calibrated against the model
using the task CALIB in amplitude and phase (SOLMODE = 'A\&P') on a 30 second
timescale.  The antenna feed-leakage term corrections from LPCAL (see below) were applied
using
SPLIT and the parameter DOPOL = 3 to guarantee that the corrections were applied
to both the parallel and cross-hand correlations.  A final pass of CALIB self-calibrated
the data again against the Stokes-$I$ model to remove any errors
by the initial self-calibration prior to D-term removal.  For circular polarization
calibration via {\em gain transfer} it is important that both passes of
CALIB described above were performed only making the rigorous assumption\footnote{For
the total intensity and linear polarization imaging described in \citet{LH05} we
assumed there was no Stokes-$V$ in these CALIB steps by assuming
$RR = I_{model}$ and $LL = I_{model}$.}
that $(RR+LL)/2 = I_{model}$.  This left as un-calibrated
any gain offsets between the right (R) and left (L) hand feeds at each antenna.
The {\em gain transfer} calibration described in \S{\ref{s:gains}} was
used to remove these offsets.  After correction of the $R/L$ antenna gains,
final images were made in all four Stokes parameters.

It is important to note that the antenna feed leakage corrections were initially
determined on each source individually using LPCAL, and the final corrections (which
were applied to the data) were determined by taking the median of these
after removing obvious outliers.  Our typical feed leakage terms (D-terms) had amplitude of
about 0.02 with the largest being 0.05, and we estimate that we have made these
corrections accurate to about $\pm0.005$ .  Following HW99, we estimate any uncertainty
in our fractional circular polarization measurements, $m_c$, due to leakage
correction errors to be of the order $\Delta m_c < 0.001 m_L + 1\times 10^{-5}$, where $m_L$ is the
fractional linear polarization of the map peak (typically a few percent or
less).  Thus any uncertainty due to leakage feed corrections is very much smaller
than the uncertainty due to the gain calibration and can be neglected.

\subsection{R/L Gain Calibration}
\label{s:gains}
As described in HW99 and \citet{HAW01}, the {\em gain transfer} technique
requires that we initially make no assumption about circular polarization during
self-calibration, i.e., we only assume $(RR+LL)/2 = I_{model}$ as discussed above.
This will leave as un-calibrated a small relative R/L complex gain ratio at each VLBA
antenna. The essential idea is to calibrate this gain ratio with a smoothed set of
R/L complex antenna gains determined from all 18 sources in each 24 hour observing
session.

To determine the smoothed set of antenna gains that will ultimately be applied
to the data, we first solve for the R/L complex gain ratios by
self-calibrating each source assuming no circular polarization in that
source ($RR = I_{model}$ and $LL = I_{model}$).   These solutions are not
applied to the data, but rather they are collected into a single calibration table
which we smooth by using a six-hour running median with all sources equally
weighted.  Solutions from sources with
low gain SNR or very strong circular polarization (greater than 0.5\%) are
prevented from contributing to the running median\footnote{Such strongly
polarized sources were very rare; however, when we did find one, we simply
repeated the calibration with the source in question prevented from contributing
to the median gain solution}.  To account
for the possibility of smaller levels of real circular polarization showing up in the
raw R/L antenna gains of individual sources, we subtract from the raw R/L gain ratios for
each source a single mean offset (taken across all times, antennas, and IFs
for that source) from the running median.  The sum of the applied offsets for all
sources in a given 24 hour observing session is constrained to be zero.  This constraint
guarantees quick
convergence of an iterative algorithm which generates a trial set of smoothed gains,
computes mean offsets for each source, takes those offsets out of the data, and then
repeats the process.  The constraint also underscores the fundamental assumption of
the gain transfer technique: that any single observing session will have an average
circular polarization of zero across all the sources in the sample\footnote{AGN as a
whole have no preferred sign of circular polarization, and thus we expect the net
circular polarization from a large sample of AGN to be close to zero.}.  This assumption
may introduce a small overall gain bias in any individual epoch; however, our Monte
Carlo method of estimating the uncertainty due to the gain calibration (described below)
naturally accounts for this possibility.

We note that the procedure described above for accounting for sources with small levels
of real circular polarization differs from the procedure used in
HW99 and \citet{HAW01}.  In those papers, we prevented such sources
from contributing to the smoothed gains through an iterative process where the strongest
polarized sources were removed first, the smoothed gains recomputed, and the process
repeated until no source contributing to the smooth gains appeared to have real circular
polarization.  In those papers, we found the final results were insensitive to the
detailed steps of this process; however, the process was time consuming and difficult to
automate.  For the MOJAVE program, we needed a procedure that was automated, reproducible,
and could be easily tested by a Monte Carlo simulation we developed.  As a check that the
new procedure gave the same results as our previous methods, we reprocessed the antenna gain
solutions from the 40 source observation analyzed in \citet{HAW01}.  For the circularly
polarized sources in that experiment, we found an RMS R/L gain difference of $0.017$\%
between our technique and the original.  This difference is approximately one-third of the
estimated gain uncertainty of that experiment, verifying that the two techniques indeed
give effectively the same results.

\subsubsection{Uncertainty in Measured CP due to Gain Transfer}
In the appendix of \citet{HAW01}, we laid out a prescription for estimating
the uncertainty in circular polarization measurements made by gain transfer.
Here we improve on that work by using a Monte Carlo simulation to test our
methods and obtain more accurate uncertainty estimates.  As in \citet{HAW01},
we focus on amplitude gain errors, as those
dominate the uncertainty in the technique.  Phase errors will not contribute
to first order for point-like sources, and for sources with significant extended
structure, phase errors will produce anti-symmetric structure around the phase
center.  In section \S{\ref{s:phase}} we discuss possible phase
effects and tests to constrain/resolve those issues.

The Monte Carlo simulation we developed generates 18 point sources (Most sources in
our sample are dominated by a point source, e.g., \citet{LH05}.), and for each
source generates a circular polarization
with magnitude drawn at random from the distribution of observed
mean gain offsets in the real gain table.  The sign ($+$ or $-$) of the CP signal
for each source is determined at random.  Finally, a Gaussian deviate representative of our
measurement uncertainty is added to each value.  Each of these generated CP values is
then stored for comparison to the ``measured'' value at the end of a Monte Carlo trial.

To generate the gain table for the simulation, a random slowly-varying gain curve
is created for each antenna/IF with amplitudes and timescales similar to those
in the real gain curves.  To this slowly-varying gain curve is added randomly-generated
scan-to-scan variations of the same range of amplitudes as seen for the corresponding
source in the real gain table.  Finally, the generated CP values are added to the
gain curves as offsets for each source.  This step is important to mimic the real gain
table, which was generated from the assumption that each source has no circular
polarization, so the real raw gain table contains 'corrections' to the circular
polarization of each source in addition to the actual gain variations.

Finally, this randomly-generated 'raw' gain table is processed by exactly the same
procedure we used on the real gain curve.  A smoothed version of the raw table is
produced by our standard procedures and compared with the raw gain table.  If the procedures
work correctly, the only differences between the raw and smoothed tables should be
the circular polarizations set at the beginning of the simulation trial.  Therefore,
by comparing the two tables, we can estimate the uncertainty
due to the gain calibration procedure.  We repeat this simulation for
1000 independent trials to generate an estimate of the uncertainty for each source.

To safe-guard against the possibility that an unusual combination of circularly
polarized sources in our data has made our results more uncertain than the above
simulation indicates, we also ran a second version of the Monte Carlo described above.
The only difference is that the initial circular polarization of each source in the
simulation is forced to match the mean gain offset of the corresponding source in the real
gain table (plus a Gaussian deviate to represent measurement uncertainty).
To these values, we add a randomly-determined overall bias (the same for all 18 sources)
to account for the possibility that the circular polarizations in the sample do not
average to zero.  The size of this gain bias is determined from the RMS of the mean
offset values in the real gain table.  The larger of the uncertainties from the first and
second version of the Monte Carlo simulations are used to represent the total gain uncertainty for a
given source; however, in practice the two Monte Carlo simulations give very nearly the
same estimate of uncertainty.

Figure \ref{f:gains} shows a plot of uncertainty in fractional circular polarization
measurements due to gain calibration as a function of the peak Stokes-$I$ flux for
the 133 sources in the MOJAVE sample. The trend
of larger uncertainties for weaker sources is expected and is driven by the larger
self-calibration gain uncertainties in the weaker sources; however, it will be
important to consider
this trend when interpreting the results in \S{\ref{s:res}}.  Note that most sources
in the sample have uncertainty contributions from gain calibration very close to
0.1\%, and given that the
gain uncertainties dominate the total uncertainty (the much smaller RMS map noise is
added in quadrature), our procedures are able to robustly detect CP signals as small
as 0.3\% of the Stokes-$I$ peak at the 3-$\sigma$ level.

\subsection{Phase Calibration for Circular Polarization}
\label{s:phase}
Because there is no phase connection between sources in our sample (they are well separated
in the sky and fringe-fit separately), the gain calibration technique described above is
only effective for calibrating the antenna amplitude gains.  We do allow the possibility of
a phase correction in the gain transfer technique (the median phase difference between
the right and left hands is treated in a similar fashion to the median amplitude ratio.); however
the applied corrections are generally $\lesssim 0.1$ degrees for each antenna/IF combination.

For effective relative R/L phase correction at each antenna we rely first on the pulse
calibration corrections, then
on fringe fitting and the initial phase calibration in CALIB to increase coherence before
averaging the data in time in DIFMAP.  Each of these steps provides separate phase corrections
for the left and right-hand feeds at each antenna.  The pulse calibration makes no assumptions
about circular polarization in the source; however both fringe fitting and the initial phase
calibration to increase coherence make their corrections assuming the source is point-like.
Any real circular polarization on a point source will not appear as a gain error in the antenna
phases (HW99), so for most sources in our sample phase corrections found in this fashion
will remove only genuine R/L phase gain errors from the data.

If a source has significant extended structure in Stokes-$I$, the right and left hand systems of the array
may receive slightly different solutions in the algorithm's attempt to find point source solutions
to the right and left hand data separately.  Most likely this will generate a small difference in the
phase centers of the right and left hand data.  While this small difference will not affect
the Stokes-$I$ image, it may be enough to generate an anti-symmetric stripe of circular
polarization across the strongest parts of the source. The effect is only apparent in the
circular polarization image because Stokes-$V$ is formed by taking the difference of the
right and left hand data ($V = (RR-LL)/2$), whereas Stokes-$I$ is formed from their average.
See panels (a) and (c) of Figure \ref{f:pz_nocp} for examples of this effect.

To account for this effect and other possible phase errors in our data, we have generated two
circular polarization images of each source.  The first image is generated from
the gain transfer calibration alone and may contain unaccounted-for phase errors.  The second
image is generated after an additional round of self-calibration of the phases alone.  During this
self-calibration of the phases, we assume there is no circular polarization in the data so
both the $RR$ and $LL$ data are compared to $I_{model}$ to find the antenna phase corrections
for the right and left-hand feeds separately.  This self-calibration will very effectively
remove any genuine R/L phase gain errors remaining in the data, as shown in panels (b) and (d) of
Figure \ref{f:pz_nocp}.  The large majority
of the sources in our sample are core-dominated and thus show no significant differences
between the images produced before and after this phase self-calibration; however, an
important minority of sources with significant extended structure do show differences
(usually both in their noise levels and in their apparent circular polarization) between
the two sets of images.

Figure \ref{f:pz_cp} illustrates two cases where the antenna phase errors
have been removed to reveal real circular polarization.  It is important to know to what
extent (if any) the resulting circular polarization in these images has been modified by
the phase calibration assuming zero circular polarization.  In general, phase differences
between RR and LL due to real circular polarization do not show up in a consistent antenna-based
fashion (as would phase gain errors); however it is possible that some combination of
antenna-based phase gains could be found by the self-calibration algorithm which minimizes
these real phase differences and thus modifies the real circular polarization of the source.

We have tested for this effect by (1) removing all circular polarization
from a real source (by setting $RR=LL=I$), (2) adding known amounts of circular
polarization, and then (3) self-calibrating the phases of that data assuming no
circular polarization.  We have done this for several sources from our observations,
and in general we find that the known circular polarization is not modified very
much by this procedure.
We usually see some reduction in amplitude from a few percent up to about 10\% of the
known CP value, and the location of the CP peak may be shifted a fraction of a beam-width to
better align with the map peak (assuming it is near the map
peak to begin with).  In some sources with a broad core-region (core $+$ barely resolved
component), a single component of CP is spread out by this self-calibration to encompass
the whole region.
Large shifts in CP location (more
than a half a beam-width) or large reductions in amplitude (more than 20\%) do not seem
possible with this procedure.  For example, our tests on the images in Figure \ref{f:pz_cp}
show that real circular polarization that began on the first jet component would remain there
after phase self-calibration and would not be shifted to the core; thus, we can be confident
that the real circular polarization in this epoch of 3C\,273 is actually associated with
the core as revealed in panel (b).

\section{Results and Discussion}
\label{s:res}

Our first-epoch circular polarization measurements are summarized in Tables \ref{t:cp_tab1},
\ref{t:cp_tab2}, and \ref{t:cp_tab3}.  Table 1 lists the linear and circular polarization
measurements at the peak of each map.  In most cases the Stokes-I map peak corresponds
to the unresolved VLBI core; however, in a handful of cases the core and map peak do not
coincide.  In most of these cases, the peak appears at the location of a strong component
within the jet, and we take the core to be the compact component at the base of the jet.
The core of 0923$+$392 is taken at the location identified by \citet{AMZ93}.  The core of 0316$+$413
(3C\,84) is a special case and is discussed further in \S{\ref{s:sources}}.
The off-peak jet cores are listed in Table 2, along with measurements of their linear
and circular polarization at the core location.  Table 3 lists the circular
polarization measured on all jet components in our sample with peak fluxes stronger than
300 mJy beam$^{-1}$.  We chose 300 mJy beam$^{-1}$ as a reasonable cutoff for the jet component
analysis as our typical uncertainty levels would still allow very strong circular
polarization of $m_c = 0.5$\% to be detected on such a component.

The linear polarization was always measured at the location of the Stokes-I peak (or at
the location of the VLBI core or jet component in the cases of Tables 2 and 3 respectively).
Limits are given for the linear polarization when the measured value does not exceed five
times the map RMS noise.  Circular polarization was also measured at the location of the
Stokes-I peak (or at the location of the VLBI core or jet component for Tables 2 or 3)
unless there was a CP peak within two pixels of the Stokes-I measurement location.  In that
case, the value of the CP peak was used instead.  This procedure was used because with our
phase calibration (see section 2), there may be a small positional shift in the
location of the circular polarization relative to the Stokes-I.  The raw circular polarization
measurements and their uncertainties are given in Tables 1, 2, and 3.  The uncertainties
includes both the gain transfer uncertainty from our Monte Carlo simulation (see \S{2})
and the RMS noise\footnote{The RMS noise was calculated over a large number of pixels away
for source structure.  For each pixel, we used largest value of CP within 2 pixels of
that location.  This procedure increases our estimate of the RMS noise, but it is consistent
with the procedure we use for measuring CP from a Stokes-I peak.} added in quadrature.  It is important
to note that the fractional circular polarization values reported in the Tables include
both measurements and limits.  While we do actually measure a CP value in every case, most
of these measurements are not significant, and thus we have chosen a $2\sigma$ cutoff, below
which we report a limit on the value of $|m_c|$.  The quoted limits on $|m_c|$ are the larger of
$\left(|\mathrm{V}| + 1\sigma_\mathrm{V}\right)$ or $2\sigma_\mathrm{V}$.  The raw values of $m_c$
in the cases that are less significant than $2 \sigma_\mathrm{V}$ can be computed from
the $V$ and $I$ intensities reported in the Tables.

We note that we have applied exactly the same procedures to the cores and jet components in
Tables 2 and 3 as we applied to the map peak values; and, as argued in \S{\ref{s:gains}}, the
interpretation of the gain uncertainties is most robust for the map peak.  Therefore, we are
somewhat less confident of the uncertainties reported in Tables 2 and 3; however, we believe
that the gain uncertainties will still accurately account for the possibility of any overall
gain bias in the CP map. The gain uncertainties do not account for the possibility of
random (but not net) gain errors creating small amounts of spurious circular polarization off
the map peak; however, the RMS noise level in the map naturally takes this into account.  We notice
no particular trend for map noise to be scattered preferentially onto weaker source structure,
and we therefore believe that the combination of gain uncertainties and map RMS noise provides
a reasonable uncertainty estimate for the off-peak circular polarization reported in Tables 2
and 3.

\subsection{Circular Polarization on VLBI Cores}
\label{s:cores}

Figure \ref{f:core_dist} shows the distribution of core circular polarization for all the sources
in our sample.  For the purposes of this figure we use measured values rather than limits
in every case.  We detect circular polarization in 17 out of 133 jet cores at the $3\sigma$ level
or higher significance.  Of the 115 sources with overall uncertainties less than $0.15$\%, we
measure circular polarization stronger than $0.3$\% of the local Stokes-I in 17 of those
sources (all but two of those at $\geq 3\sigma$ significance).  The existence of strong
circular polarization ($\geq 0.3$\%) in about 15\% of our sample is somewhat higher than
the rate of strong circular polarization observed by \citet{HAW01} with the VLBA at 5 GHz.
They found CP stronger than 0.3\% in only 2 out of 36 sources (6\%); however, these rates of
strong circular polarization detection are compatible according to Poisson statistics.

We examined the circular polarization distribution between source types \citep[using the
optical classifications given by][]{LH05}, and a Gehan's
Generalized Wilcoxon  test from the ASURV package \citep{LIF92} shows no significant difference
between the distributions for quasars and BL Lac objects.
Unfortunately there are simply too few galaxies in the MOJAVE sample for a robust test
of their distribution.  A Gehan's test also shows no significance difference between
the EGRET and non-EGRET detected sources \citep[as defined in][]{LH05} in terms of
their CP levels.

Due to the large number of limits in our sample, we used a Kendall-Tau test for censored data
\citep{LIF92} to test for correlations between circular polarization and various other source properties.
Table \ref{t:tau} collects the results of these tests.  The Kendall-Tau test gives the probability that
{\em no correlation} exists between fractional circular polarization, $|m_c|$, and the property
in question.  Small values in the Table indicate that a correlation may exist, and we take $p < 0.05$
as suggestive of a correlation.  Of the twenty properties tested in Table \ref{t:tau}, we find
only three with $p \lesssim 0.05$: 15 GHz jet flux ($S_{jet}$), 15 GHz core/jet flux ratio ($R_{jet}$),
and 15 GHz core/total flux ratio ($R_{total}$) as defined by \citet{LH05}.  These three properties
are closely related through the prominence of the VLBI scale jet, so we really only have one property
with a possible correlation at the $p < 0.05$ level, consistent with pure chance.  Removing 3C\,273 from
consideration raises $p$ to $0.12$ for $S_{jet}$ and $0.09$ for $R_{total}$, so we don't find
strong evidence that circular polarization is correlated with any of these properties.

While Galactic latitude, $|b|$, is not correlated with circular polarization, it is interesting
that of the 12 sources near the galactic plane ($|b| < 10$ degrees), none of them have
circular polarization at $\geq 2\sigma$.  Treating these sources and those with $|b| > 10$ degrees as
separate samples, a Kolomogorov-Smirnov gives only a 5.3\% chance they are drawn
from the same distribution.  Again, this
low probability may be spurious as there seems little reason the galactic plane should
interfere with circular polarization signals.  In particular, any galactic Faraday
rotation should have no noticeable effect on circular polarization \citep[e.g.,][]{HAW01}.

One property that might naturally be expected to correlate with core circular polarization is
core linear polarization.  In the most common models of circular polarization production \citep{WH03},
both linear and circular polarization increase with increasing field order; however, the two may not
be correlated in a few scenarios: (1) if circular polarization is produced through a high internal
Faraday depth model \citep[e.g.,][]{RB02,  BF02}, (2) if the linear polarization is externally
depolarized \citep[e.g.,][]{HAW01}, or (3) if the circular and linear polarization we observe are
produced in two distinct regions of the jet but are unresolved in what we see as the VLBI core.

Figure \ref{f:core_mcvml} is a plot of absolute fractional core circular polarization, $|m_c|$, versus
fractional core linear polarization, $m_l$.  There is no clear relation between the two, and the
obvious lack of
correlation between linear and circular polarization is similar to the 5 GHz VLBA results of \citet{HAW01},
and the integrated 5 GHz measurements from the ATCA by \citet{RNS00}.  The persistence of no trend
at higher frequency suggests that the complicating factors (lack of
resolution/Faraday effects) are to some degree self-similar and scale-up with observing frequency.

Although there is no correlation between linear and circular polarization for jet cores, Figure
\ref{f:core_mcvml} shows that for quasars and BL Lac objects, linear polarization almost always
exceeds circular polarization.  This does not seem to be the case for radio galaxies, which
generally have weak linear polarization.  \citet{BBFM01} have noted
that the nearby very low luminosity radio sources Sgr A* and M81* are circularly polarized but
show no linear polarization. These authors have suggested that the $m_c/m_l$ ratio is distinctly
different for low luminosity AGN ($m_c/m_l >> 1$) and high luminosity AGN ($m_c/m_l << 1$).
This relation was
tested further by \citet{BFM02}, who failed to find additional LLAGN with circular polarization;
however, if we consider the lowest redshift (and hence lowest luminosity) sources in our sample
(which are admittedly quite different in character from the Sgr A* and M81*), we find that
$m_c/m_l >> 1$ for M87 (1228$+$126) and 3C\,84 (0316$+$413), the only two radio galaxies
for which we can do the comparison.
Both M87 and 3C\,84 are known to have strong external Faraday screens which can depolarize the
linear polarization \citep{W71, ZT02}, and in fact, all the radio galaxies in the MOJAVE sample
have very low levels of core
linear polarization, likely due to external depolarization.  Higher redshift sources may
have smaller $m_c/m_l$ ratios not because a different production mechanism for CP, but rather because
either (1) they have less depolarization, and/or (2) we study them with less effective linear
resolution and therefore polarization from more of the jet contributes to what we see as the "core".
If (2) is the case, it also naturally explains the lack of correlation between $m_c$ and $m_l$ in
Figure \ref{f:core_mcvml}.

\subsection{Circular Polarization of Jet Components}

Table 3 lists circular polarization measurements for all strong jet components in our sample with
peak Stokes-I $\geq 300$ mJy beam$^{-1}$.  Of the 33 jet components that meet this criteria, we detect
circular polarization in six (18\%) of them at the $3\sigma$ level or higher significance.
Five of these six
detections are on the relatively nearby sources 3C\,84 (0316$+$413) and 3C\,273 (1226$+$023) which
are discussed in more detail
in section \S{\ref{s:sources}}, and we note that 3C\,84 is a special case as the
distinction between the VLBI core and jet is not well defined and the circular polarization measured is
unusually strong ($\gtrsim 1$\%).  For these reasons, we must be cautious about drawing conclusions
about the circular polarization characteristics of jet components in general; however, it is clear
that strong jet components can be circularly polarized, and the levels of circular polarization found in
strong jet components are not clearly different from jet cores.

\subsection{Discussion of Selected Sources}
\label{s:sources}

Here we discuss a few sources in detail, and we also compare our first
epoch results with previous VLBA observations of circular polarization
(HW99; \citet{HAW01}; \citet{HW03};
\citet{HW04}) and from historical integrated measurements (1970s through
early 1980s: \citet{WdP83}, hereafter WdP83; \citet{K84}, hereafter K84;
\citet{AAP03}).  Because the latter are at lower frequency
(typically 5 GHz and below) and are integrated measurements, we only
compare to the historical circular polarization results in those cases
where the source was repeatedly (three times or more) detected in
circular polarization.

Because there is a large overlap of our sample with the sources observed
at 5 GHz with the VLBA on 1996 December 15 by \citet{HAW01}, we list results
from the latter study alongside our 15 GHz results in Table \ref{t:ba18}. We
discuss some of the sources common to both of these studies in more detail
below.  The large number of limits from both programs makes the results in
Table \ref{t:ba18} difficult to interpret; however, it appears that CP of the
core can vary with time and/or frequency.

\paragraph{0316$+$413 (3C\,84):}

The left-hand panels of Figure \ref{f:3c84_comp} show the central 5 milli-arcseconds
from our March 2003 MOJAVE observations of 3C\,84.  The image reveals four distinct
spots of circular polarization within this central region of the source.
Without a clear component to identify as the core, we have taken spot "D" as
the location of the VLBI core. \citet{HW04} found that even the locations
corresponding to spots "C" and "B" were optically thick in December of 1997, so
the precise definition of the jet core in this source is not clear.  With only
one frequency, we don't know where the $\tau\sim 1$ surface is in 2003, so it
is possible that all three spots of negative circular polarization: B, C, and
D, may be part of an extended optically thick region; however, for the purposes
of this paper, we have chosen to treat spots B and C as locations in the
jet, along with spot A.

Both spots D and C are similarly circularly polarized at the $-0.4$\% level; however,
the strength of the circular polarization increases sharply in spot B to
$-1.4$\%.  A little further out in the jet, but at quite a different position
angle, spot A is more than $2$\% circularly polarized with {\em positive}
circular polarization.  \citet{HW04} attribute this change in sign to
the stark change in opacity they observe between spots B and A.

The circular polarization we observe is very similar to that seen
by \citet{HW04} in observations more than 5 years earlier; however,
spot D is entirely new and was not present in December of 1997.  A version of
the 15 GHz image from \citet{HW04} (convolved with a matching beam to
our MOJAVE observations) appears in the right-hand panels of Figure
\ref{f:3c84_comp}.  The 1997 image is registered with the 2003 image by
matching contours on the northern-most end of the central region shown here.
The red dots in the right-hand panels of Figure \ref{f:3c84_comp} have
the same relative positions as the blue stars which mark the locations of
spots A, B, and C in 2003, but they have been shifted by 0.2 milli-arcseconds
to the North and by 0.2 milli-arcseconds to the West.  It is clear that
the {\em relative} positions of these circularly polarized spots has not
changed significantly in 5$+$ years; however, they do appear to move outward
together at a sub-luminal rate of roughly $0.06c$, almost directly
South-East.  This speed is, of course, dependent on our alignment of
the images; however, it is similar to the speed found by \citet{DKR98} at 43 GHz
for two of the three component they studied within this region, although
those two components were moving almost due South.

\citet{HW04} found spot A to have $m_c = +3.2\pm0.1$\% which is somewhat
larger than the $m_c = +2.3\pm0.3$\% we measure five years later.  The intensity
at the location of spot A has faded by a factor of three, from 1.2 Jy beam$^{-1}$ to
0.4 Jy beam$^{-1}$, so presumably the region has become more optically thin.  Multi-frequency
observations are necessary to test if this is indeed the case, but if so, it is
remarkable that this location has remained as strongly circularly polarized as
it has.  \citet{HW04} argued that Faraday conversion is responsible for the
high levels of circular polarization detected here.  Faraday conversion
is strongly dependent on optical depth, and we would eventually expect a
significant decline in circular polarization from this mechanism if this
region of the jet is becoming increasingly optically thin \citep[e.g.,][]{WH03}.
It will be particularly interesting to follow the changes in spot A in future
MOJAVE epochs.

\citet{HW04} found spots B and C to have $m_c = -0.7\pm0.1$\% and
$m_c = -0.6\pm0.1$\% respectively at 15 GHz.  Spot B has essentially doubled
in circular polarization percentage in 5 years, and Spot C has not
changed significantly within our uncertainties.  Neither location has
changed much in Stokes-$I$ intensity.

If the alignment between the epochs is correct, the lack of large changes in
the circular polarization distribution between spots A, B, and C is surprising
given their motion.  The magnetic field order and particles responsible for
producing this circular polarization must be in motion, however, this motion
seems to have preserved the particular field characteristics responsible for
the strong levels of circular polarization we observe.

\paragraph{0430$+$052 (3C\,120):}
We find a limit of $|m_c| < 0.29$\% in the core at 15 GHz in February of 2003.
HW99 also found limits between 0.2 and 0.3\%  in the core at 15 GHz throughout
1996.  Both signs of circular polarization (with a clear preference
for positive CP) were detected intermittently in the historical integrated
measurements (WdP83 and K84) at 8.9 GHz and below.

\paragraph{0528$+$134:}
We find a limit of $|m_c| < 0.18$\% in the core at 15 GHz in February of 2003.
HW99 detected multiple epochs of significant circular polarization at about the
$m_c = -0.5$\% level in the core at 15 GHz during 1996.  During 1996,
0528$+$134 was undergoing a large outburst, peaking with an integrated
core flux 4.6 times larger than observed in February 2003.

\paragraph{0607$-$157:}
We find a limit of $|m_c| < 0.21$\% in the core at 15 GHz in March of
2003. \citet{HW03} measured only $m_c = -0.23\pm0.12$\% in the core at 15 GHz in
January of 1998, but they simultaneously detected $m_c = -0.65\pm0.07$\% in
the core at 8 GHz.  \citet{HAW01} detected circular polarization in core of
$m_c = -0.18\pm0.05$ \% at 5 GHz thirteen months earlier in
December of 1996.  \citet{AAP03} found strong circular polarization,
$m_c = -0.23\pm0.04$\% at 4.8 GHz in integrated measurements with the
UMRAO during the period 2001--2002.

\paragraph{0735$+$178:}
We measure $m_c = -0.30\pm0.11$ \% in the core at 15 GHz in November of
2002.  HW99 did not detect core circular polarization at 15 GHz at any epoch
during 1996 with a typical limit of $|m_c| < 0.4$\%.

\paragraph{0851$+$202 (OJ287):}
We measure $m_c = -0.20\pm0.08$ \% in the core at 15 GHz in October of 2002.
HW99 did not detect core circular polarization at 15 GHz at any epoch during
1996, with a typical limit of $|m_c| < 0.2$\%.  \citet{HAW01} also found a limit of
$|m_c| < 0.10$\% for the core circular polarization at 5 GHz in December of
1996.  In older integrated measurements during the period 1978--1983,
\citet{AAP03} found strong CP at 4.8 and 8.0 GHz with $m_c = -0.37\pm0.12$\%
and $m_c = -0.22\pm0.03$\% respectively.

\paragraph{0923$+$392 (4C 39.25):}
We find a limit of $|m_c| < 0.18$\% in the core at 15 GHz in October of 2002.
In older integrated measurements during the period 1978--1983,
\citet{AAP03} found strong CP at 4.8 GHz with $m_c = -0.21\pm0.05$\%.

\paragraph{1127$-$145:}
We measure $m_c = -0.24\pm0.10$ \% in the core at 15 GHz in October of
2002.  Negative circular polarization was repeatedly detected at 5 GHz and
below in the historical integrated measurements cataloged by WdP83 and
K84.

\paragraph{1222$+$216:}
We find a limit of $|m_c| < 0.35$\% in the core at 15 GHz in February of 2003.
HW99 also found limits between 0.1 and 0.2\%  in the core at 15 GHz throughout
1996.

\paragraph{1226$+$023 (3C\,273):}
Figure \ref{f:3c273cp} shows the first 10 milli-arcseconds of the jet of 3C\,273
at 15 GHz in October of 2002.  The top panel is Stokes I, the middle panel is
linear polarization, and the bottom panel is circular polarization.
Positions where circular polarization measurements or limits were taken
are marked and numbered on the Stokes I map.  Figure \ref{f:3c273_cpvr} shows
the fractional circular polarization at these locations as a function of
distance from the jet core.  The absolute fractional level of the circular
polarization  appears to decrease as distance from the core increases.  This
decrease could be tied to decreasing opacity along the jet or a decrease in poloidal
magnetic field strength.

The core of 3C\,273 was found by HW99 to have $m_c \simeq -0.5$\% at three epochs during
1996.  The circular polarization appeared with the emergence of a new jet
component.  The core was very opaque at this time, and HW99 showed that the
appearance of the CP was tied to the flattening of the core's spectral
index due to the emergence of the new component.  As cataloged in WdP83
and K84, 3C\,273 was repeatedly detected in integrated measurements
during the 1970s and early 80s at 8 GHz and below where it had a consistently
negative sign of circular polarization.

\paragraph{1228$+$126 (M87):}
Figure \ref{f:m87ipv.ps} shows our 15 GHz VLBA image of the jet of M87
in February of 2003.  No linear polarization is detected, but strong
circular polarization, $m_c = -0.49\pm0.10$\% is seen at the peak of
the core.

M87 was studied at the VLA in April of 2000 by \citet{BFM02}, who
found no significant circular polarization at 8 GHz,
$m_c = 0.01\pm0.10$; however, they did find linear polarization
at the $1.7$\% level.  These results are the opposite of our finding
of strong circular polarization and no linear polarization.  It is
important to consider that their VLA results include contributions
from the jet well beyond the region imaged by the VLBA, and the
linear polarization they observe must come from beyond the
depolarizing screen which covers most of the VLBA scale jet
\citep{ZT02}.  The fractional circular polarization they
measure may be diluted by emission beyond the VLBA scale, or the
circular polarization may simply vary with frequency or time.
Regardless, the differing linear polarizations observed between the
VLBA and VLA scale images underscores our point from \S{\ref{s:cores}}
that at least some of the lack of correlation between linear
and circular polarization for jet cores could be due to the two
quantities being produced in different regions of the unresolved
jet core.

M87 has been detected intermittently in integrated circular
polarization (with both signs) at 5 GHz and below in the
historical observations cataloged by WdP83 and K84.

\paragraph{1253$-$055 (3C\,279):}
Figure \ref{f:3c279cp} shows the first 5 milli-arcseconds of the jet of 3C\,279
at 15 GHz in November of 2002.  The top panel is Stokes I, the middle panel is
linear polarization, and the bottom panel is circular polarization.
Positions where circular polarization measurements or limits were taken
are marked and numbered on the Stokes I map.  We detect strong circular
polarization of $m_c=+0.30\pm0.08$\% only on the jet core.  The sign and magnitude
of the circular polarization is consistent with five epochs of CP detections
by HW99 at 15 GHz during 1996.  At 5 GHz, \citet{HAW01} did not detect circular
polarization with a limit of $|m_c| < 0.11$\%.  Historically, 3C\,279 has
repeatedly and consistently shown a positive sign for its circular
polarization in integrated measurements at 8 GHz and below during the 1970s
and early 80s (WdP83 and K84); however, \citet{AAP03} found $m_c = -0.17\pm0.01$\%
at 4.8 GHz with the UMRAO during 2001-2002.  This change in sign at the lower
frequency is particularly interesting as no change has yet occurred at 15 GHz.

\paragraph{1308$+$326:}
We find a limit of $|m_c| < 0.21$\% in the core at 15 GHz in November of 2003.
HW99 also found limits between 0.1 and 0.2\%  in the core at 15 GHz throughout
1996.

\paragraph{1510$-$089:}
We measure $m_c = +0.20\pm0.09$ \% in the core at 15 GHz in November of
2002.  HW99 did not detect circular polarization in the core at 15 GHz
during 1996: three epochs had limits of 0.1 to 0.2\%, and
for two epochs they reported peak CP values of $+0.2\pm0.2$ \%.  \citet{HAW01} find a
limit of $|m_c| < 0.24$\% in the core at 5 GHz in December of 1996.

\paragraph{1641$+$399 (3C\,345):}
We find a limit of $|m_c| < 0.19$\% in the core at 15 GHz in November of
2002.  \citet{HW03} detected circular polarization in the core
at 15 GHz, $m_c = -0.38\pm012$\%, in January of 1998.  The same location
was optically thick at 8 GHz where they did not detect any circular
polarization in simultaneous measurements.  3C\,345 was detected
intermittently (and with both signs) at 8 GHz and below in the historical
integrated measurement compiled by WdP83.

\paragraph{1730$-$130:}
We find a limit of $|m_c| < 0.17$\% in the core at 15 GHz in October of
2002.  The historical integrated measurements (WdP83 and K84) show
several intermittent detections of positive circular polarization
at 5 GHz and below.

\paragraph{1749$+$096:}
We find a limit of $|m_c| < 0.14$\% in the core at 15 GHz in May of 2002.
HW99 also found typical limits of 0.2\%  in the core at 15 GHz throughout
1996.

\paragraph{1928$+$738:}
We measure $m_c = +0.19\pm0.09$ \% in the core at 15 GHz in June of
2002.  HW99 also measured circular polarization in the core 15 GHz at
the limit of their sensitivity.  They found $m_c =+0.3\pm0.1$ \% in
two epochs, $m_c =+0.2\pm0.1$ \% in one epoch, and limits of 0.1
and 0.2\% in their other two epochs during 1996.

\paragraph{2134$+$004:}
Figure \ref{f:2134cp} shows our observations of 2134$+$004 at 15 GHz in May of 2002.
The top panel is Stokes I, the middle panel is linear polarization, and the
bottom panel is circular polarization.  Positions where circular polarization
measurements or limits were taken are marked and numbered on the Stokes I map.
We detect strong circular polarization of $m_c=+0.35\pm0.12$\% on the jet core,
we measure a similar value at location number 1 of $m_c=+0.29\pm0.12$\%, and
we find a limit on the CP at location 2 to be $|m_c| < 0.34$\%.

In recent integrated measurements from 2001--2002, \citet{AAP03} found strong
circular polarization at 4.8 GHz with $m_c = +0.62\pm0.06$\%.  2134$+$004 was
also detected repeatedly in the historical integrated measurements at 5 GHz
and below, always with positive circular polarization (WdP83 and K84).

\paragraph{2145$+$067:}
We find a limit of $|m_c| < 0.26$\% in the core at 15 GHz in March of
2003.  In recent integrated measurement from 2001--2002, \citet{AAP03}
measured $m_c = -0.50\pm0.07$\% at 4.8 GHz with the UMRAO.
The historical integrated measurements (WdP83 and K84) show
consistent, repeated detections of negative circular polarization
at 5 GHz and below.  Given this historical record of detection, it is
interesting that the 'limit' reported here came from measurements that
had a significance of $1.98\sigma$, just barely less than the
$2\sigma$ cutoff we had for defining measurements versus limits.
If we allow $1.98\sigma$ as a measurement, then our fractional level
for the core is $m_c = -0.17\pm0.09$ \%.

\paragraph{2200$+$420:}
We find a limit of $|m_c| < 0.21$\% in the core at 15 GHz in June of 2002.
WdP83 reported intermittent detections of both signs of circular polarization
at frequencies below 5 GHz in the historical integrated measurements.

\paragraph{2230$+$114 (CTA102):}
We find a limit of $|m_c| < 0.19$\% in the core at 15 GHz in October of 2002.
\citet{HAW01} measured $m_c = -0.09\pm0.04$\% in the core at 5 GHz in December
of 1996.  CTA102 was intermittently detected (with both signs of CP) in the
historical integrated measurements (WdP83; K84) at 5 GHz and below.

\paragraph{2251$+$158 (3C\,454.3):}
Figure \ref{f:2251cp} shows the first 4 milli-arcseconds of the jet of
3C\,454.3 at 15 GHz in March of 2003.  Positions where circular polarization
measurements or limits were taken are marked and numbered on the
Stokes I map.  The first jet component (position 1) appears slightly
more strongly circularly polarized, $m_c = +0.38\pm0.10$\%, than
the core, $m_c = +0.23\pm0.10$\%; however, an alternative calibration
assuming no circular polarization for calibrating the phases (see
\S{\ref{s:phase}})  spreads the circular polarization between the core
and component 1 more evenly with both having $m_c \simeq 0.3$\%.

\citet{HAW01} measured $m_c = +0.08\pm0.04$\% in the core of 3C\,454.3 at 5 GHz
in December of 1996, and the source was repeatedly detected with positive
circular polarization at 8 GHz and below in the integrated historical
measurements compiled by WdP83 and K84.

\paragraph{2345$-$167:}
We detect $m_c = +0.40\pm0.11$\% in the core at 15 GHz in October of 2002.
\citet{HAW01} found a limit of $|m_c| < 0.11$\% in the core at 5 GHz in December of
1996.  WdP83 and K84 reported intermittent detections (but always with positive
sign) of circular polarization in integrated measurements at frequencies of 5 GHz
and below.

\section{Conclusions}
\label{s:conclude}

In this paper, we report first-epoch circular polarization results for
each of the 133 AGN in the MOJAVE program to monitor the VLBA structure and
polarization of a flux limited sample of AGN jets at 15 GHz.  We found strong
circular polarization ($\geq 0.3$\%) in approximately 15\% of our sample.  The
circular polarization we detected was usually associated with jet cores, typically
in the range of $0.3-0.5$\% of the local Stokes-$I$.
Our ability to detect circular polarization on jet components was limited
by their relatively weak flux density; however, we did find a few strong jet
components to be circularly polarized.

We found no significant correlations between fractional circular polarization
of jet cores and optical source type, redshift, EGRET detections, linear polarization,
or a number of other pc-scale jet properties listed in Table 4.  While there was no
clear relation between linear and circular polarization of jet cores, quasars
and BL Lac objects consistently had more linear than circular polarization, and
the two galaxies with detected circular polarization had much more circular
than linear polarization.  We suggested that higher redshift sources may have
smaller $m_c/m_l$ ratios, not because they produce circular polarization
differently from low redshift sources, but rather because either (1) they
have less depolarization, and/or (2) we study them with less effective
linear resolution and therefore polarization from more of the jet
contributes to what we see as the ``core''.

The lack of clear correlations between circular polarization and other
source properties is intriguing, particularly given that some sources,
3C\,273 and 3C\,279 in particular, show consistently high levels of
circular polarization over years and decades of observation, even though
they are both known to be highly variable with regular ($\sim$ annual)
ejections of new jet components.  It seems highly unlikely that these
sources would be circularly polarized in a very similar fashion from
epoch to epoch due only to stochastic factors (such as a chance favorable
magnetic field configuration).  It may be that we simply detect too few sources in
circular polarization to find subtle trends between circular polarization
and other source properties, or perhaps, even with the VLBA we are
not resolving the core region well enough to see where the circular
polarization is being generated.

In addition to looking for overall trends in the circular polarization of
our sample, we also examined the structure of a few sources in detail.
3C\,84 was particularly interesting, showing the strongest circular
polarization in our sample, up to $2$\% in one location.  Our image
revealed essentially the same structure as \citet{HW04} found in data
taken five years earlier; however, a new circularly polarized spot
appeared in our data, and the older spots of circular polarization,
while maintaining the same relative positions, appear to have moved outward
in the jet at a rate of approximately 0.06$c$.  3C\,273 also showed multiple
locations of circular polarization in its jet, and we were able to plot circular
polarization measurements for several strong jet components.  This plot revealed
that, in 3C\,273, fractional circular polarization gets weaker at greater
distances from the jet core.  Multi-frequency observations will be necessary
for both sources so we can use opacity measurements to help understand the
trends of circular polarization with jet position.

\acknowledgments

The authors wish to acknowledge the contributions of the other members of the
MOJAVE team: Hugh and Margo Aller, Tigran Arshakian, Marshall Cohen, Matthias Kadler,
Ken Kellermann, Yuri Kovalev, Andrei Lobanov, Eduardo Ros, Rene Vermeulen, and Tony Zensus.
This work has been supported by Denison University, the National Radio Astronomy
Observatory, and NSF grant 0406923-AST.  This research was also supported by
an award from Research Corporation.