VLA Observations of a New Population of Blazars

Abstract: We present the first deep VLA radio images of flat-spectrum radio quasars (FSRQ) with multiwavelength emission properties similar to those of BL Lacs with synchrotron X-rays. Our observations of twenty-five of these sources show that their radio morphologies are similar to those of other radio quasars. However, their range of extended powers is more similar to that of BL Lacertae objects (BL Lacs) and extends down to the low values typical of FR I radio galaxies. Five out of our nine lobe-dominated sources have extended radio powers in the range typical of both FR I and FR II radio galaxies, but their extended radio structure is clearly FR II-like. Therefore, we have not yet found a large population of radio quasars hosted by FR Is. Two thirds of our sources have a core-dominated radio morpholgy and thus X-rays likely dominated by the jet. We find that their ratios of radio core to total X-ray luminosity are low and in the regime indicative of synchrotron X-rays. This result shows that also blazars with strong emission lines can produce jets of high-energy synchrotron emission and undermines at least in part the ``blazar sequence'' scenario which advocates that particle Compton cooling by an external radiation field governs the frequency of the synchrotron emission peak.

\section{Introduction}

Blazars are one of the most extreme classes of active galactic nuclei
(AGN). Their broad-band emission extends from radio up to gama-ray
frequencies, varies rapidly and irregularly, is strongly polarized,
and shows a core-dominated morphology. Furthermore, apparent
superluminal component speeds are often observed in these sources. The
properties of blazars can be best explained if we assume that, in
fact, they are classical double-lobed radio galaxies with their jets
seen at small angles relative to our line of sight and so subject to
strong relativistic beaming [see \citet{Urry95} for a review].

Based on their optical spectra, blazars are currently separated into
BL Lacertae objects (BL Lacs), defined to have no or only very weak
emission lines, and flat-spectrum radio quasars (FSRQ), defined to
have strong narrow and broad emission lines [but see \citet{L04} for a
suggested revision of the blazar classification scheme]. Within the
framework of unified schemes, BL Lacs are identified primarily with
beamed low-power Fanaroff-Riley type I \citep[FR I;][]{Fan74} radio
galaxies, whereas the high-luminosity Fanaroff-Riley type II (FR II)
radio galaxies are assumed to constitute the parent population of
radio quasars.

Traditionally, BL Lacs were discovered in radio and X-ray surveys, and
in recent years more efficiently through joint radio-X-ray selections,
whereas radio quasars had only been looked for in radio surveys. Thus,
it is not surprising that the synchrotron emission of BL Lac jets was
found to have its maximum within a wide range of frequencies, from
IR/optical to UV/soft-X-ray energies, whereas only radio quasars with
synchrotron emission maxima at relatively low energies, and thus
X-rays dominated by inverse Compton emission, were known. This, in
fact, led to the so-called ``blazar sequence'' which proposed that the
stronger the particle Compton cooling by an external radiation field
(such as the one produced by, e.g., the accretion disk or broad
emission line region, both stronger in quasars), the lower the
frequency of the synchrotron emission peak \citep{Sam96, Fos98,
Ghi98}. In other words, radio quasars with X-rays dominated by
synchrotron emission were not expected to exist.

Two recent surveys have drastically changed our picture of blazar jet
physics. About 10\% of the FSRQ discovered in the Deep X-ray Radio
Blazar Survey \citep[DXRBS;][]{Per98, L01} and $\sim 30\%$ of the ones
identified in the {\it ROSAT} All-Sky Survey (RASS) - Green Bank
\citep[RGB;][]{LM98, LM99} have been shown to have spectral energy
distributions (SEDs) and broad-band emission properties similar to
those of high-energy peaked BL Lacs \citep{P03}. Subsequent
investigations of the {\it EINSTEIN} Medium Sensitivity Survey (EMSS)
and the Slew survey proved that these ``X-ray-strong'' FSRQ had indeed
gone undetected in previous surveys \citep{Wol01b, Per01}.

In this paper we present the first deep radio interferometric
observations of this new class of blazars (Sections 2 and 3) and
discuss in the light of our observational results their nature and
their relation to ``standard'' FSRQ and BL Lacs (Section 4). Our
conclusions are presented in Section 5. For consistency with previous
work we have assumed throughout this paper cosmological parameters
$H_0 = 50$ km s$^{-1}$ Mpc$^{-1}$ and $q_0 = 0$. Spectral indices have
been defined as $S_\nu \propto \nu^{-\alpha}$.

\section{Observations and Data Reduction} \label{observe}

We have selected for observations with the NRAO\footnote{The National
Radio Astronomy Observatory is a facility of the National Science
Foundation operated under cooperative agreement by Associated
Universities, Inc.} Very Large Array (VLA) flat-spectrum ($\alpha_{\rm
r} \mathrel{\hbox{\rlap{\hbox{\lower4pt\hbox{$\sim$}}}\hbox{$<$}}} 0.5$) radio quasars from the DXRBS and RGB surveys with
radio-to-X-ray spectral indices $\alpha_{\rm rx} \mathrel{\hbox{\rlap{\hbox{\lower4pt\hbox{$\sim$}}}\hbox{$<$}}} 0.78$. The latter
constraint makes it likely that our sources have their X-rays
dominated by synchrotron rather than inverse Compton emission
\citep{P96}. We have preferably chosen sources at redshifts $z \mathrel{\hbox{\rlap{\hbox{\lower4pt\hbox{$\sim$}}}\hbox{$<$}}} 1$
in order to minimize the effects of surface brightness dimming
\citep{Per93} and to increase the probability of detecting large-scale
extended radio emission.

We were granted observing time for 25 sources. Their general
properties and log of observations are listed in Table
\ref{general}. The columns are: (1) object name; (2) redshift; (3)
total radio flux at 1.4 GHz from the NRAO/VLA Sky Survey (NVSS); (4)
radio spectral index between 1.4 and 5 GHz, calculated from the sum of
the fluxes of all NVSS sources within a $3'$ radius (corresponding
roughly to the beam size of the GB6 survey) and the total flux from
the GB6 and PMN surveys for northern and southern sources,
respectively; (5) unabsorbed {\it ROSAT} X-ray flux at 1 keV,
calculated using an X-ray spectral index derived from hardness ratios;
(6) $k$-corrected radio-to-X-ray spectral index between 5 GHz and 1
keV; (7) observation date; (8) VLA configuration; and (9) program
number. The observed sample contains 19/85 (or 22\%) low-redshift ($z
\mathrel{\hbox{\rlap{\hbox{\lower4pt\hbox{$\sim$}}}\hbox{$<$}}} 1$) FSRQ with $\alpha_{\rm rx} \mathrel{\hbox{\rlap{\hbox{\lower4pt\hbox{$\sim$}}}\hbox{$<$}}} 0.78$ from the combined RGB and
DXRBS samples. We expect the properties of the observed sources to be
representative of this entire subsample, since their total radio and
X-ray luminosities span similar ranges.

We observed with the VLA at 1.425 GHz in continuum mode. We imaged all
sources with the VLA in A configuration, and were able to gain C array
images for 14 of the 25 sources. Three or four scans of several
minutes length yielding a total exposure time of $\sim 45$ minutes
were interleaved with one minute scans on a suitable secondary VLA
flux calibrator for each source. Scans were spaced to optimize
coverage in the $(u,v)$ plane. Multiple observations of 3C 48, 3C 147,
or 3C 286 were used to flux-calibrate the maps.

The data were processed with the Astronomical Image Processing System
(AIPS; version 31DEC03) package. A model based on Clean components was
used to start the self-calibration process. Phase-only
self-calibration was used for the first three iterations and amplitude
and phase self-calibration for the last one or two iterations. The
task IMAGR with robust weighting (\mbox{ROBUST}=0.5) was used to
generate the maps and Clean components. We found it useful to combine
the self-calibrated data sets from the A and C arrays in order to
increase the sensitivity and to improve the $(u,v)$ plane coverage in
only 4/14 cases. Our results are shown in Figure 1 and the
corresponding map parameters are listed in Table \ref{map}, where the
columns are as follows: (1) object name; (2) VLA configuration; (3)
beam size; (4) position angle; (5) image rms; (6) peak flux; and (7)
corresponding panel in Figure 1. The dynamic ranges (peak/noise) of
the maps lie between $\sim 100-5000$.

The high spatial resolution provided by the A array data allowed a
good measure of the core emission, and thus in combination with the C
array data an accurate determination of the extended emission flux
density, morphology and size. In Table \ref{radioprop} we list the
radio properties of the sources as derived from our observations. The
columns are: (1) object name; (2) and (3) position of the unresolved
core, measured on A array map before start of self-calibration
process; (4) core flux density, measured on final A array map; (5)
extended flux density, derived from the total flux density measured on
C array map (where available and if the source was resolved, else
measured on A array map) substracting the core flux density given in
column (4); (6) core luminosity $k$-corrected with spectral index
$\alpha_{\rm r}=0$; (7) extended luminosity $k$-corrected with
spectral index $\alpha_{\rm r}=0.8$; (8) radio core dominance
parameter $R$, defined as $R=L_{\rm core}/L_{\rm ext}$, where $L_{\rm
core}$ and $L_{\rm ext}$ are the core and extended luminosities,
respectively; (9) ratio of NVSS flux to total flux density measured on
C array map (where available, else measured on A array map); and (10)
largest angular extent (LAS), measured from the core to the peak in
the extended radio structure [following \citet{Mur93}] on C array map
(where available, else on A array map). We give $1\sigma$ errors and
$2\sigma$ upper limits. The errors in the core flux densities are on
average $\sim 0.07$ mJy, close to the expected thermal noise. The
cumulative errors in the extended flux densities depend upon the solid
angular extent of the measured flux. For a considerable number of our
sources (11/25 objects) we have only A array data available, which
might miss extended flux. We have assessed this effect by using the
distribution of values in column (6), i.e., the ratio of NVSS flux to
total flux density measured. This has a median of 1.02 and
$\sigma=0.21$. Based on this we have assumed the NVSS flux as the
correct total flux for sources with ratios $\mathrel{\hbox{\rlap{\hbox{\lower4pt\hbox{$\sim$}}}\hbox{$>$}}} 1.5$ (the median value
plus $2 \sigma$). This applies to only one source, namely
WGAJ0435$-$0811, for which we list in column (5) the 'corrected'
extended flux.

\section{Discussion of Individual Sources}

We give now for the individual sources a brief discussion of their
general properties and a description of the radio images.

\vspace*{0.2cm}

{\sl RGB J0112$+$3818. $-$} An optical spectrum of this source was
published in \citet{LM98} based on which it was classified as a
quasar. This source is also part of the CLASS blazar survey
\citep{Marcha01, Cac02}. Faint extended radio emission is visible in
both images, and most pronounced in the A array image which shows two
faint lobes close to the core (Fig. 1(a)).

\vspace*{0.1cm}

{\sl RGB J0141$+$3923. $-$} Spectroscopic observations were published
in \citet{Marz03}. This source is also part of the CLASS blazar survey
\citep{Marcha01, Cac02}. The source remains unresolved in the radio on
the scales imaged (Fig. 1(b)).

\vspace*{0.1cm}

{\sl RGB J0254$+$3931. $-$} Optical spectra can be found in
\citet{Marcha96} and \citet{LM98} who classified this source as a
broad-line radio galaxy and quasar respectively. This source shows a
one-sided jet on the combined A$+$C image which is most prominent on
the A array image and appears to end in a lobe or hot spot (Figs. 1
(c) and (d)). A similar morphology is visible also on the map
published by \citet{Tay96} based on VLA B array snapshot observations
at 4.535 GHz.

\vspace*{0.1cm}

{\sl RGB J1629$+$4008. $-$} An optical spectrum of this source was
published in \citet{LM98} who classified it as a broad-line radio
galaxy based on an absolute optical magnitude of $M>-22$ mag. X-ray
observations of this source were obtained with {\sl Beppo}SAX by
\citet{P02}. The source remains unresolved in the radio on the scales
imaged (Fig. 1(e)). Similarly, the VLA B array map from the FIRST
survey at 20 cm shows only an unresolved core.

\vspace*{0.1cm}

{\sl RGB J2229$+$3057. $-$} \citet{LM98} presented an optical spectrum
of this source and classified it as a quasar. The source is extended
on both A and C array images (Figs. 1(f) and (g)). In particular, on
the A array image two well-defined lobes are visible at each side of
the nucleus, the closest lobe being fairly extended with a hot spot.

\vspace*{0.1cm}

{\sl RGB J2256$+$2618. $-$} An optical spectrum of this source was
published by \citet{LM98} who classified it as a broad-line radio
galaxy based on an absolute optical magnitude of $M>-22$ mag.
Extended emission on scales of several tens of arcseconds is detected
in the C array image (Fig. 1(h)).

\vspace*{0.1cm}

{\sl RGB J2308$+$2008. $-$} An optical spectrum can be found in
\citet{LM98}. Based on this and an absolute optical magnitude of
$M>-22$ mag these authors classified this source as a broad-line radio
galaxy. \citet{Falco98} give a slightly different redshift for this
source of $z=0.2342$ based on four emission lines. The source is
unresolved on C array scales, however, a faint lobe is visible close
to the core on the A array image (Fig. 1(i)).

\vspace*{0.1cm}

{\sl RGB J2318$+$3048. $-$} \citet{LM98} published an optical spectrum
of this source and classified it as a broad-line radio galaxy based on
an absolute optical magnitude of $M>-22$ mag. This source is clearly
extended on the C array image (Fig. 1(j)), however, only the core is
detected on A array scales.

\vspace*{0.1cm}

{\sl WGA J0106$-$1034. $-$} Optical spectra of this source can be
found in \citet{Wol98} and \citet{L01} who classify it as a
broad-emission line AGN and quasar, respectively. This source is also
part of the FIRST Bright Quasar Survey \citep{Beck01} and the Sloan
Digital Sky Survey \citep{Schn03}. The A array image of this source
(Fig. 1(k)) shows extended emission with hot spots on each side of the
nucleus, making it clearly an FR II source. However, one of the two
hot spots is considerably stronger than the other, and this source
appears unusual in that the radio jet we see points towards the weaker
of the two lobes. The source is unresolved on C array scales. Our ATCA
snapshot observations of this source give a steeper spectral index of
$\alpha_{\rm r}=0.67$ (between 4.8 and 8.6 GHz) than the one listed in
Table \ref{general} ($\alpha_{\rm r}=0.35$).

\vspace*{0.1cm}

{\sl WGA J0110$-$1647. $-$} Optical spectra of this source can be
found in \citet{Per98} and \citet{Cac00} who classify it as a quasar
and broad-line object, respectively. An extended jet-like feature is
clearly visible on the A array image (Fig. 1(l)) a few arcseconds away
from the core.

\vspace*{0.1cm}

{\sl WGA J0126$-$0500. $-$} An optical spectrum of this source was
published by \citet{L01} who classify it as a quasar. The A array
image (Fig. 1(m)) shows it to be a double-lobed source with some
extended emission also around the core.

\vspace*{0.1cm}

{\sl WGA J0227$-$0847. $-$} An optical spectrum of this source can be
found in \citet{L01} who classify it as a quasar. This source is also
part of the Sloan Digital Sky Survey \citep{Schn03}. Our ATCA snapshot
observations of this source give a spectral index of $\alpha_{\rm
r}=0.37$ between 4.8 and 8.6 GHz, which is steeper than the one listed
in Table \ref{general} ($\alpha_{\rm r}=-0.34$). This is either due to
variability or to the fact that this source is a Gigahertz
peaked-spectrum (GPS) candidate \citep[e.g.,][]{ODea98}. The source
remains unresolved in the radio on the scales imaged
(Fig. 1(n)). Similarly, the VLA B array map from the FIRST survey at
20 cm shows only an unresolved core.

\vspace*{0.1cm}

{\sl WGA J0259$+$1926. $-$} An optical spectrum of this source can be
found in \citet{Cac00} who classify it as a broad-line object. Our own
high-quality spectrum of this source will be published shortly. We
classify this source as a quasar. The A array image of this source
(Fig. 1(o)) shows a jet-like feature extending a few arcseconds away
from the core. A slightly extended morphology is visible even on the
scales imaged with the VLBA at 8.4 GHz \citep{Ojha04}.

\vspace*{0.1cm}

{\sl WGA J0304$+$0002. $-$} Optical spectra of this source were
published by \citet{Per98} and \citet{Cac00} who classified it as a
quasar and a broad-line object, respectively. This source is also part
of the Sloan Digital Sky Survey \citep{Schn03}. The source appears
extended on both the A and C array images (Fig. 1(p) and (q)) with two
well-defined lobes on each side of the nucleus and some fuzzy extended
emission on larger scales off one of the lobes. Our ATCA snapshot
observations of this source give a steeper spectral index of
$\alpha_{\rm r}=0.71$ (between 4.8 and 8.6 GHz) than the one listed in
Table \ref{general} ($\alpha_{\rm r}=0.40$).

\vspace*{0.1cm}

{\sl WGA J0435$-$0811. $-$} \citet{Per98} published an optical
spectrum and classified this source as a quasar. It is unresolved on A
array scales (Fig. 1(r)), however, most likely its extended emission
is resolved out. The NVSS image is slightly extended and the ratio of
NVSS to total measured flux on our image for this source is the
highest in the sample (see Table \ref{radio}, column (6)).

\vspace*{0.1cm}

{\sl WGA J0447$-$0322. $-$} Optical spectra of this source can be
found in \citet{Per98} and \citet{Cac00} who classified it as a quasar
and broad-line object, respectively. This source remains unresolved on
the scales imaged (Fig. 1(s)). Our ATCA snapshot observations give a
steeper spectral index of $\alpha_{\rm r}=0.91$ (between 4.8 and 8.6
GHz) than the one listed in Table \ref{general} ($\alpha_{\rm
r}=0.47$).

\vspace*{0.1cm}

{\sl WGA J0544$-$2241. $-$} \citet{Per98} published an optical
spectrum and classified this source as a quasar. The source remains
unresolved in the radio on the scales imaged (Fig. 1(t)). Our ATCA
snapshot observations give a spectral index of $\alpha_{\rm r}=0.31$
(between 4.8 and 8.6 GHz), which is steeper than the one listed in
Table \ref{general} ($\alpha_{\rm r}=-0.44$). This is either due to
variability or to the fact that this source is a Gigahertz
peaked-spectrum (GPS) candidate \citep[e.g.,][]{ODea98}.

\vspace*{0.1cm}

{\sl WGA J1026$+$6746. $-$} \citet{L01} published an optical spectrum
and classified this source as a quasar. The A array image (Fig. 1(u))
shows jet-like extended structure over several tens of arcseconds
emerging from the core to both sides and an extended lobe to the far
left.

\vspace*{0.1cm}

{\sl WGA J1457$-$2818. $-$} An optical spectrum of this source can be
found in \citet{L01} who classify it as a quasar. Slightly extended
emission away from the core is visible on A array scales (Fig. 1(v)).
In better accordance with the observed morphology, our ATCA snapshot
observations of this source give a steeper spectral index of
$\alpha_{\rm r}=0.72$ (between 4.8 and 8.6 GHz) than the one listed in
Table \ref{general} ($\alpha_{\rm r}=0.40$).

\vspace*{0.1cm}

{\sl WGA J2239$-$0631. $-$} \citet{L01} published an optical spectrum
and classified this source as a quasar. The source is extended on the
C array image (Fig. 1(w)) with two lobes to each side of the nucleus,
one of them having a stronger hot spot than the other. Only the
nucleus is detected on A array scales.

\vspace*{0.1cm}

{\sl WGA J2320$+$0032. $-$} An optical spectrum of this source can be
found in \citet{L01} who classify it as a quasar. This source is also
part of the Sloan Digital Sky Survey \citep{Schn03}. The source
remains unresolved in the radio on the scales imaged
(Fig. 1(x)). Similarly, the VLA B array map from the FIRST survey at
20 cm shows only an unresolved core.

\vspace*{0.1cm}

{\sl WGA J2322$+$2114. $-$} Optical spectra of this source can be
found in \citet{Wol98} and \citet{Per98} who classify it as a
broad-emission line AGN and quasar, respectively. The A array image
(Fig. 1(y)) shows clearly a double-lobed structure and a pronounced
core.

\vspace*{0.1cm}

{\sl WGA J2347$+$0852. $-$} \citet{Per98} published an optical
spectrum of this source and classified it as a quasar. This source is
extended in both our A and C array images (Figs. 1(z) and (aa)). The
X-shape of the galaxy is apparent only in the C array image, whereas
on A array scales we observe a central core with two almost
equidistant lobes. Both lobes have hot spots, however, one lobe is
brighter than the other.

\vspace*{0.1cm}

{\sl PKS 0256$-$005. $-$} The optical spectra of this source presented
by \citet{Rich80b} and \citet{Osm94} show several strong broad
emission lines. This source is also part of the FIRST Bright Quasar
Survey \citep{Beck01} and the Sloan Digital Sky Survey
\citep{Schn03}. The source remains unresolved in the radio on the
scales imaged (Fig. 1(ab)). Similarly, the VLA B array map from the
FIRST survey at 20 cm shows only an unresolved core.

\vspace*{0.1cm}

{\sl 0959$+$68W1. $-$} Optical spectra of this source are available
from the {\sl Hubble Space Telescope} archive. These show several
broad emission lines. This source is also part of the CLASS blazar
survey \citep{Marcha01, Cac02}. The A array map (Fig. 1(ac)) shows
jet-like features and jet blobs extending over several tens of
arcseconds on both sides of the pronounced core.

\section{Results and Discussion}

In this section we want to address three main questions: 1. How do the
general radio properties of these newly discovered strong-lined
blazars compare to those of ``classical'' radio quasars and BL Lacs?
2. Which type of radio galaxies form their parent population? and
3. Is the X-ray emission of these ``X-ray strong'' radio quasars
dominated by the jet and synchrotron in origin?

\subsection{The Comparison Samples} \label{samples}

For our studies we have selected for comparison other types of
radio-loud AGN, namely, steep-spectrum radio quasars (SSRQ),
``standard'' FSRQ and BL Lacs.

The comparison sample of SSRQ has been drawn from the Molonglo Quasar
Sample (MQS) excluding compact steep-spectrum (CSS) quasars and
GHz-peaked sources (74 objects). \citet{Kap98} published for this
sample total and core radio fluxes, largest angular sizes, and radio
spectral indices (between 0.4 and 5 GHz) from observations with the
VLA in BnA and CnB configurations at 1.4, 5, and 8.4 GHz, and with the
Parkes-Tidbinbilla Interferometer (PTI) at 2.3 GHz. Largest angular
sizes were measured between the brightest peaks in the extended radio
structure. Improved {\sl ROSAT} X-ray fluxes were published by
\citet{Bak95b}.

The comparison sample of FSRQ has been selected from the list of
\citet{Mur93}. These authors published core and extended radio fluxes
and largest angular sizes for 56 sources from the 1 Jy catalogue
\citep{Kuehr81} as obtained from observations with the VLA in A and B
configurations at 1.4 GHz. The largest angular sizes were measured
from the core to the peak in the extended radio structure. We have
compiled X-ray fluxes at 1 keV and radio spectral indices between 1.4
and 5 GHz for these sources from the multiwavelength catalogue of
\citet{P97}.

Additionally, we have chosen for comparison the well-known 1 Jy and
EMSS BL Lac samples, since these have complete dedicated VLA
observations \citep[e.g.,][]{Per93,Cas99,Rec00,Rec01} as well as in
depth {\sl ROSAT} observations \citep[e.g.,][]{Per96a,Urry96,Rec99}
available. \citet{Per93} give for the EMSS sample largest linear sizes
defined as the sum of two straight lines from the core to the
outermost $3\sigma$ contours, with each line intersecting the
brightest hotspot, if any, in the extended structure. Simultaneous
radio spectral indices between 1.4 and 5 GHz for this sample have been
taken from \citet{Cava04}. \citet{Rec01} give for the 1 Jy sample
largest linear sizes defined as in \citet{Per93}. Radio spectral
indices between 2.7 and 5 GHz for this sample have been taken from
\citet{Sti93c}.

We have transformed core and extended radio fluxes to 1.4 GHz using a
radio spectral index of $\alpha_{\rm r} = 0$ and 0.8,
respectively. These same indices have been used to $k$-correct the
luminosities. We have transformed band fluxes to an energy of 1 keV
and $k$-corrected the corresponding luminosities using an X-ray
spectral index of $\alpha_{\rm x} = 1.2$.

\subsection{The Radio Properties} \label{radio}

In this section we investigate the radio properties of the newly
discovered ``X-ray strong'' radio quasars and compare them to those of
``classical'' radio quasars and BL Lacs.

\subsubsection{Extended Radio Powers and the Radio Core Dominance Parameter}

In Fig. \ref{logRlext} we have plotted for our sources and the
comparison samples introduced in Section \ref{samples} the extended
radio power at 1.4 GHz versus the radio core dominance parameter at
this frequency. The extended radio power is believed to be radiated
isotropically and, therefore, unaffected by relativistic beaming. Thus
this quantity is a good measure of the source's intrinsic jet
power. Also marked (horizontal solid lines) are the ranges of extended
power typical of FR I and FR II radio galaxies \citep{Owen94}. For a
constant extended radio power an increase in radio core dominance
parameter is indicative of stronger relativistic beaming, since
beaming affects only the core emission. For this reason the radio core
dominance parameter $R$ is assumed to be a suitable orientation
indicator (the higher $R$, the smaller the viewing angle, i.e., the
angle between jet and observer's line of sight).

Fig. \ref{logRlext} shows that, as expected, ``classical'' radio
quasars (SSRQ and FSRQ, filled triangles and squares, respectively)
have FR II-like extended powers, with only one FSRQ from the sample of
\citet{Mur93} having an extended power more typical of an FR I. On the
other hand, the BL Lacs from the samples considered here span a range
in extended power larger than that of quasars, populating both the FR
I- and FR II-regimes. Our sources differ from the radio quasars
plotted here. Their range of intrinsic radio power is similar to that
of BL Lacs and extends down to the values typical of FR Is, reaching
about three orders of magnitude lower values than ``classical'' radio
quasars.

The low radio powers of our sources are partly due to the fact that
they are selected from surveys with factors of $\sim20-80$ lower radio
flux limits than the comparison radio quasars. This selection,
however, gives us a most important result. We have found either the
lowest luminosity examples of FR II radio galaxies or quasars hosted
by FR Is. The first possibility would indicate that morphology and
intrinsic radio power of radio galaxies are not as tightly connected
as first found by \citet{Fan74}. The second possibility would mean
that FR I radio galaxies can host quasar nuclei, i.e., they can have
broad emission line regions at their centers. To our knowledge, so far
only three such objects are known, J2114$+$820, a broad-line radio
galaxy \citep{Lara99}, E1821$+$643, a quasar associated with a
subarcsecond jet in an FR I radio structure \citep{Blu01}, and the
well-known FR I broad-line radio galaxy 3C 120
\citep[e.g.,][]{Walker87}. In any case, our finding of a large number
of radio quasars with strong emission lines but relatively low radio
powers can be regarded as another observational fact against a close
relationship between emission line luminosity and jet power in AGN
\citep{Raw91}, in addition to the existence of radio-quiet quasars
(i.e., sources with strong emission lines but extremely weak radio
jets) and BL Lacs (i.e., sources with very weak emission lines but
strong radio jets).

Fig. \ref{logRlext} also shows that our sources with extended powers
in the range populated by the comparison radio quasars have $R$ values
intermediate between those of SSRQ and FSRQ. This result, however,
does not mean that ``X-ray strong'' quasars {\it are} the less beamed
versions of ``classical'' FSRQ. For a correct interpretation we need
to take the above mentioned selection effects into account. FSRQ
selected in surveys with lower radio flux limits (such as our sources)
will have on average lower {\it total} radio luminosities than FSRQ
selected in surveys with higher radio flux limits. Then, considering
only a narrow range, i.e., an almost constant value, of relatively
high {\it extended} radio powers, the former will have on average
lower {\it core} luminosities and so will appear less beamed than the
latter. In fact, the sources in our sample reach higher $R$ values
than the FSRQ in the comparison sample, but at lower extended radio
powers.

\subsubsection{Largest Linear Sizes}

In Fig. \ref{logRLLS} we have plotted for our sources and the
comparison samples introduced in Section \ref{samples} the largest
linear size (LLS) versus the radio core dominance parameter $R$ at 1.4
GHz. We note that the largest linear sizes were measured differently
in each sample. In the FSRQ and our sources they were measured from
the core to the peak in the extended structure, whereas in the SSRQ
and BL Lac samples they were measured between the peaks of the
extended structure on either side of the nucleus. However, we did not
attempt to correct the values for the SSRQ and BL Lacs since, on the
one hand, such a correction is not straightforward (e.g., assumptions
about the intrinsic radio structure are required) and, on the other
hand, we are interested mainly in comparisons between our sources and
the FSRQ.

As expected from unified schemes, Fig. \ref{logRLLS} shows that the
largest linear sizes of radio-loud AGN decrease with increasing radio
core dominance parameter, i.e., with decreasing angle between jet and
observer's line of sight. Our sources follow this trend very
well. Significant LLS $- R$ correlations have been reported previously
by other authors for radio quasars, but not for BL Lacs
\citep[e.g.,][]{Hine80, Kap82, Hough89, Hough99}. Following their
example we plot in Fig. \ref{logRLLS} also the relation between LLS
and $R$ expected from relativistic beaming models (dashed, dotted and
solid lines). In view of the selection effects and incompleteness in
the data we restrict ourselves to a prediction of the expected upper
and lower envelopes and do not attempt to extract an average or median
value of LLS as a function of $R$ to test beaming models.

The fits were obtained as follows. If $\theta$ is the angle between
the jet axis and the line of sight, the observed values of LLS and $R$
are given by, ${\rm LLS} = L_0 \sin \theta$ and $R = (R_{90}
B(\theta))/2$, where $L_0$ is the intrinsic largest linear size,
$R_{90}$ the radio core dominance parameter that would be observed if
the source was oriented at $\theta=90^{\circ}$, and $B(\theta) = (1 -
\beta \cos\theta)^{-(2+\alpha_{\rm r})} + (1 + \beta
\cos\theta)^{-(2+\alpha_{\rm r})}$ for the case of a continuous jet
with $\beta=v/c$, the ratio of jet velocity to the speed of light. The
relation between LLS and $R$ results from their common dependence on
viewing angle and its simulation requires initial values for $L_0$,
$R_{90}$, the Lorentz factor $\Gamma = 1/\sqrt{1-\beta^2}$, and the
radio spectral index $\alpha_{\rm r}$.

For the radio spectral index we have assumed in all cases a value of
$\alpha_{\rm r} = 0.8$, which is typical of optically thin synchrotron
emission. Radio galaxies can have different intrinsic radio core
dominance values and we have assumed values of $R_{90} = 0.022$ and
0.003, the median values obtained by \citet{Morg97} for their samples
of FR I and narrow line FR II radio galaxies, respectively. In the
case of $R_{90} = 0.003$, the upper and lower envelopes in
Fig. \ref{logRLLS} can be fitted if we assume a bulk Lorentz factor of
$\Gamma = 8$ and values of $L_0 = 1900$ and 30 kpc, respectively
(dashed lines). In the case of $R_{90} = 0.022$, a slightly lower bulk
Lorentz factor of $\Gamma = 6$ and values of $L_0 = 1400$ and 20 kpc,
respectively, are required (dotted lines). We note that a fit with a
lower radio spectral index would require a much higher value of the
bulk Lorentz factor, e.g., in the case of $\alpha_{\rm r} = 0.2$ and
$R_{90} = 0.003$ a value of $\Gamma = 15$.

Previously presented LLS $-$ $R$ relations included only
lobe-dominated radio quasars, i.e., sources viewed at relatively large
angles, but did not account for this by restricting the range of
orientation accordingly. Thus it is not suprising that fits to their
envelopes required smaller Lorentz factors of $\Gamma = 2-5$ (for
values of $R_{90}$ similar to our FR II radio galaxy case). We note
also that the sources in the samples considered here reach higher LLS
values than those included in the previous studies. We have included
in Fig. \ref{logRLLS} both quasars and BL Lacs and it appears that
they follow the same relation, but their starting points are most
probably different, which we have accounted for here by adopting two
cases of $R_{90}$. The resulting maximum intrinsic LLS values of their
parent radio galaxies are not very much different and in accordance
with the result of \citet{Led02} who found maximum sizes for both FR I
and FR II radio galaxies around 2000 kpc (with only one source
exceding this value).

\subsubsection{The Selection of Blazars}

Past and on-going blazar surveys, such as, e.g., the 1 Jy BL Lac
survey \citep{Sti91}, CLASS blazar survey \citep{Marcha01, Cac02} and
DXRBS \citep{Per98, L01}, have selected their candidates based on
radio spectral index. A limit of $\alpha_{\rm r} = 0.5$ is generally
chosen, since radio sources are expected to become increasingly
core-dominated below this value \citep{Orr82}. Ideally, blazar
selection should be based on orientation, but methods to derive actual
viewing angles require long-term radio observations and/or
sophisticated multiwavelength data sets. A blazar selection based on
the value of the radio core dominance parameter of $R>1$ would then be
another possibility, but this would require dedicated radio
observations for a large number of sources, which are usually not
readily available. Therefore, this approach would be extremely time
consuming.

In order to assess the effect of blazar selection based on radio
spectral index relative to one possibly based on the radio core
dominance parameter we have plotted in Fig. \ref{logRalr} for our
sources and the samples introduced in Section \ref{samples} these two
quantities versus each other. Fig. \ref{logRalr} shows that indeed
radio-loud AGN appear to populate preferably the regions of either
steep radio spectral index ($\alpha_{\rm r}>0.5$) and lobe-dominance
(i.e., $R<1$) or flat radio spectral index ($\alpha_{\rm r}<0.5$) and
core-dominance (i.e., $R>1$). Only a few sources in blazar surveys
imposing a radio spectral index cut of $\alpha_{\rm r} = 0.5$ are
expected to be lobe-dominated, and similarly, only a few
core-dominated sources are expected to be missed by these surveys.

A considerable fraction of our sources (7/25 or 28\%), however, have a
flat radio spectral index $\alpha_{\rm r}\le0.5$ and {\it are}
lobe-dominated. Are we then justified to refer to these sources as
blazars? Variability undoubtedly plays a role in the explanation of
the existence of at least some of the lobe-dominated flat-spectrum
radio sources, however, it is unlikely that it can account for the
existence of all of them. For example, we have additional ATCA
snapshot observations available for 3 out of the 7 sources, and for
two of them, namely, WGAJ0106$-$1034 and WGAJ0304$+$0002, these yield
steep radio spectral indices of $\alpha_{\rm r} = 0.67$ and 0.71,
respectively. For WGAJ0126$-$0500 our ATCA observations give
$\alpha_{\rm r} = 0.39$, only slightly steeper than the value of 0.21
listed in Table \ref{general}. A more plausible explanation for the
existence of lobe-dominated flat-spectrum radio sources can instead be
gained from simulations of the relation between $\alpha_{\rm r}$ and
$R$.

From the definitions of the radio spectral index between two
frequencies $\nu_1$ and $\nu_2$ of $(S_{\nu_1}/S_{\nu_2}) \propto
(\nu_1/\nu_2)^{-\alpha_{\rm r}}$ and the relation between the radio
core dominance parameters at these two frequencies of $R_{\nu_2} =
R_{\nu_1} (\nu_1/\nu_2)^{\alpha_{\rm r_c}-\alpha_{\rm r_e}}$, where
$\alpha_{\rm r_c}$ and $\alpha_{\rm r_e}$ are the radio spectral
indices of the core and extended radio emission, respectively, we have
derived $\alpha_{\rm r} - \log R_{\nu_1}$ relations for different
values of $\alpha_{\rm r_c}$. We have assumed for $\nu_1$ and $\nu_2$
values of 1.4 and 5 GHz, respectively, and for the extended emission a
radio spectral index of $\alpha_{\rm r_e}=0.8$. Our results are shown
in Fig. \ref{logRalr} as dashed lines from top to bottom for a core
radio spectral index of $\alpha_{\rm r_c}=0.5$, 0, $-0.5$, and $-1.3$.

Firstly, we note that the relation between $\alpha_{\rm r}$ and $\log
R$ is not a linear one \citep[see also][]{Kem86}, contrary to what the
fit of \citet{Cava04} to their Fig. 4 suggests. Secondly, our results
show that lobe-dominated flat-spectrum radio sources can indeed exist,
if the core emission has a radio spectral index of $\alpha_{\rm
r_c}\la0$. In general, the lower the core radio spectral index, the
less core flux relative to the extended emission is required to render
the total radio spectral index of the source below a value of
$\alpha_{\rm r_c}=0.5$. Therefore, a blazar selection and definition
based on a radio core dominance parameter value of $R\ga1$ appears to
be more robust than one based on a radio spectral index of
$\alpha_{\rm r}\la0.5$. But, as Fig. \ref{logRalr} shows, a selection
based on a flat radio spectral index is expected to include
lobe-dominated sources rather than to miss blazars.

In this regard, we note that the result that FR I radio galaxies have
higher $R$ values than FR II radio galaxies \citep[e.g.,][]{Morg97}
could be due to selection effects. If the cores of FR II radio
galaxies had on average lower radio spectral indices than the ones of
FR I radio galaxies, in samples selecting at relatively high
frequencies and thus sensitive to flat-spectrum sources, the first
will appear on average less core-dominated than the latter.

\subsection{The Parent Population} \label{parent}

According to current unified schemes for radio-loud AGN, weak-lined
blazars, the BL Lacs, are harboured mainly by low-power radio
galaxies, the FR Is, whereas radio quasars are hosted by
high-luminosity FR II radio galaxies. In Section \ref{radio} we have
shown that our sources extend the range of intrinsic radio powers of
quasars down to the low values typical of FR Is and BL Lacs. What then
is the parent population of our newly discovered ``X-ray strong''
radio quasars?

In our sample, 9/25 sources are lobe-dominated and in principle a
morphological classification of their extended radio structure is
possible. Out of these, 4 sources have FR II-like extended radio
powers with two of them, WGAJ0106$-$1034 and WGAJ2322$+$2114, having
also a clear FR II-like morphology. They show extended lobes with
prominent hot spots in their outer regions. The morphology of
WGAJ0304$+$0002 is not clear. One of the lobes appears to have an
outer hot spot and so is indicative of an FR II. The other lobe,
however, might be more FR I-like. A higher quality map of this object
would be helpful in this case. The source WGAJ1026$+$6746 shows an
extended jet structure with some extended lobe emission. We
tentatively interpret its morphology as FR II-like. Another 5 sources
have extended radio powers typical of both FR Is and FR IIs. Out of
these, RGBJ2229$+$3057, WGAJ0126$-$0500, and WGAJ2347$+$0852 show a
clear FR II-like morphology, especially on the A array map. The
extended structure of WGAJ2239$-$0631 appears to be more FR II-like,
however, only a B array image can assess this with certainty. The
source WGAJ0435$-$0811 needs to be imaged in B or C array in order to
classify its extended emission, since this is most likely resolved out
on the scales imaged (see Section \ref{observe}).

Unfortunately, all sources in our sample with extended powers typical
of FR Is (8/25 objects) are core-dominated and so a morphological
classification of their extended emission is not possible. Sources
with extended powers in this low regime and lobe-dominated morphology
are the most interesting ones, since they can clarify if a large
population of low-luminosity FR II radio galaxies and/or quasars
hosted by FR Is exist. In this respect, it is vital to observe with
the VLA more ``X-ray strong'' radio quasars and in general radio
quasars with low total radio powers in order to search for this kind
of objects.

\subsection{The Origin of the X-ray Emission} \label{origin}

The X-ray emission of radio-loud AGN, contrary to their emission at
radio frequencies, is believed to be made up of at least three
components with different origins: an isotropic one, produced by the
ambient hot gas (presumably from the group or cluster associated with
the source or from the galaxy), and two anisotropic components, one
generally interpreted as inverse Compton emission from the accretion
disk corona, which is present also in radio-quiet quasars and is
expected to be seen only in sources with their nuclei unobscured by
the central dusty torus, and the other associated with the powerful
relativistic jet and thus subject to beaming.

In order to constrain the origin of the X-ray emission in our sources
we have plotted in Fig. \ref{lrclx} their ratio of radio core to total
X-ray luminosity versus the radio core dominance parameter $R$. A plot
such as the one displayed in Fig. \ref{lrclx} is most suitable for our
purpose. Relativistic beaming is believed to increase only the core
jet emission, both radio and X-ray, and leave any isotropically
emitted extended emission unaffected. If the total X-ray emission
$L_{\rm x}$ is dominated by an unbeamed component, the ratio $L_{\rm
core}/L_{\rm x}$ will then increase with $R$ until the beamed jet
emission starts to dominate the X-rays. Once this is the case, the
ratio $L_{\rm core}/L_{\rm x}$ will stay constant and reflect simply
the ratio intrinsic to the jet (assuming that the radio and X-ray
beaming parameters are the same).

To investigate if this behaviour is indeed observed we have included
for comparison in Fig. \ref{lrclx} also other types of radio-loud AGN,
selected as described in Section \ref{samples}. From Fig. \ref{lrclx}
we see that, as expected, the X-ray emission in SSRQ (filled
triangles), which are lobe-dominated sources and so are presumably
seen at larger angles, appears to be dominated by an unbeamed
component \citep{Kem93, Bak95b, Sie98}. Their ratios of radio core
emission to total X-ray luminosity $L_{\rm core}/L_{\rm x}$ increase
with $R$. On the other hand, in FSRQ (filled squares), which are
core-dominated sources and so are presumably seen at small angles, the
X-ray emission appears to be indeed dominated by the beamed jet. Their
ratios of $L_{\rm core}/L_{\rm x}$ are independent of $R$. In these
sources these ratios should then reflect the intrinsic jet SED, where
values of $\log L_{\rm core}/L_{\rm x} \ge 6$ and $<6$ are indicative
of X-rays dominated by inverse Compton and synchrotron emission,
respectively \citep{P96}.

However, so far, BL Lacs have been the only blazars known to have jets
with synchrotron X-rays. This is apparent in Fig. \ref{lrclx} where
these sources (open squares), almost all core-dominated, span a large
range in $L_{\rm core}/L_{\rm x}$ and reach relatively low values. On
the other hand, all but two of the ``classical'' FSRQ plotted here
have ratios of $L_{\rm core}/L_{\rm x} \ge 6$.

What is the origin of the X-ray emission in our newly discovered
``X-ray strong'' radio quasars? Their distribution in Fig. \ref{lrclx}
shows that a third (9/25 objects) fall in the range $R<1$ and follow
the behaviour exhibited by SSRQ. Therefore, the X-ray emission in
these sources is most likely {\it not} dominated by the jet, but comes
rather from the ambient hot gas and/or the the accretion disk. We have
performed correlation analysis for the following three cases, all SSRQ
(63 sources), SSRQ with detected X-ray emission (24 sources), and our
lobe-dominated sources. In the first case censored data was present
and we used the ASURV analysis package \citep{Iso86}. In all cases the
observed correlations are highly significant ($P>99.9\%$) with slopes
$0.79\pm0.11$, $0.88\pm0.09$, and $0.93\pm0.28$, respectively, not
significantly different from one (as expected for a roughly constant
ratio of $L_{\rm ext}/L_{\rm x}$). The best fits are shown as dashed
lines in Fig. \ref{lrclx}. At $R=1$, these correlations give the
average $\log L_{\rm core}/L_{\rm x}$ ratios intrinsic to the jets of
the corresponding population. Fig. \ref{lrclx} shows that this value
is for our sources $<6$, i.e., suggestive of synchrotron X-rays.

The majority of our sources (16/25 objects) have a core-dominated
radio morphology and, contrary to the majority of ``classical'' FSRQ,
their $L_{\rm core}/L_{\rm x}$ ratios extend well into the region
populated by high-energy peaked BL Lacs. In fact, almost all of our
core-dominated sources (14/16 objects) have ratios $\log L_{\rm
core}/L_{\rm x} \le 6$, indicative of synchrotron X-rays. The definite
proof of the synchrotron nature of their X-ray emission, however, will
require X-ray imaging and spectroscopy.

Dedicated X-ray observations exist or have been scheduled for 11/25
sources presented here. RGBJ1629$+$4008, the core-dominated radio
quasar with the lowest $L_{\rm core}/L_{\rm x}$ ratio in our sample,
has been observed with {\sl Beppo}SAX by \citet{P02} and found to have
X-rays dominated by synchrotron radiation with a peak frequency $\sim
2\times10^{16}$ Hz. A further four sources have been observed with
XMM-Newton in AO-3 (PI: Padovani) and six sources have been approved
for observation with {\sl Chandra} in Cycle 6 (PI: Landt). Out of
these 10 objects, 6 sources are core-dominated and have the lowest
$L_{\rm core}/L_{\rm x}$ ratios in our sample after RGBJ1629$+$4008.

\section{Summary and Conclusions}

A considerable fraction of flat-spectrum radio quasars discovered in
two recent blazar surveys (DXRBS and RGB) have multiwavelength
emission properties similar to those of BL Lacs with synchrotron
X-rays \citep{P03}. If it can be shown that these ``X-ray loud'' radio
quasars are indeed the strong-lined version of high-energy peaked BL
Lacs, it could mean that jet energetics are less dependent on emission
line luminosity than advocated so far \citep{Sam96, Fos98, Ghi98}.

In this paper we have presented deep VLA radio images obtained at 1.4
GHz in A and C configurations of 25 of these newly discovered blazars
(8 RGB and 17 DXRBS sources). We have compared the radio properties of
our sources with those of other classes of radio-loud AGN, namely,
``classical'' steep- and flat-spectrum radio quasars and BL Lacs. Our
main results can be summarized as follows.

(i) The radio morphologies of our sources are similar to those
generally observed for radio quasars and include prominent cores,
lobes and jet features. Their range of extended radio powers, however,
is similar to that of BL Lacs and extends down to the values typical
of FR I radio galaxies, reaching about three orders of magnitude lower
values than ``classical'' radio quasars. This result is partly due to
the fact that our sources are selected from surveys with relatively
low radio flux limits.

(ii) In our sample, 9/25 sources are lobe-dominated and in principle a
morphological classification of their extended radio structure is
possible. Out of these, 4 sources have FR II-like extended powers with
two having a clear and another two a tentative FR II-like
morphology. Another 5 sources have extended radio powers typical of
both FR Is and FR IIs. Out of these, 4 show a clear FR II-like
morphology, whereas one source requires further observations since its
extended emission is resolved out on the scales imaged.

(iii) As expected from unified schemes, we find for our sources and
the comparison samples of radio-loud AGN that the largest linear sizes
(LLS) decrease with increasing radio core dominance parameter
$R$. Fits based on the relativistic beaming model to the upper and
lower envelopes of the observed LLS$- R$ relation yield bulk Lorentz
factors of $\Gamma = 6-8$ and maximum and minimum intrinsic LLS values
of $\sim 2000$ and $\sim 20$ kpc, respectively.

(iv) The majority of our sources (16/25 objects) have a core-dominated
radio morphology and thus their total X-ray emission is expected to be
dominated by the jet. Out of these, 14/16 have ratios of radio core to
total X-ray luminosity $\log L_{\rm core}/L_{\rm x} \le 6$, indicative
of synchrotron X-rays \citep{P96}. The definite proof of the
synchrotron nature of their X-ray emission, however, will require
dedicated observations with {\sl Chandra} and XMM-Newton.

We plan to continue our VLA observing program of ``X-ray loud'' radio
quasars in order to find in particular lobe-dominated radio quasars
with low extended powers. These sources could then help us clarify if
a large population of quasars hosted by FR Is exists.

\acknowledgments

H.L. acknowledges financial support from the Deutsche Akademie der
Naturforscher Leopoldina grant number BMBF-LPD 9901/8-99. This
research has made use of the NASA/IPAC Extragalactic Database (NED)
which is operated by the Jet Propulsion Laboratory, California
Institute of Technology, under contract with the National Aeronautics
and Space Administration.

\input{ms2.bbl}

\include{tab1}

\include{tab2}

\include{tab3}

\begin{figure*}
\centerline{
\includegraphics[scale=0.43]{f1a.eps}
\includegraphics[scale=0.43]{f1b.eps}
}
\centerline{
\includegraphics[scale=0.43]{f1c.eps}
\includegraphics[scale=0.43]{f1d.eps}
}
\caption{\small (a) RGB J0112$+$3818, VLA A. Image rms is 0.04
mJy/beam. Image peak is 32.5 mJy/beam. (b) RGB J0141$+$3923, VLA
A. Image rms is 0.12 mJy/beam. Image peak is 104.5 mJy/beam. (c) RGB
J0254$+$3931, VLA A. Image rms is 0.08 mJy/beam. Image peak is 160.2
mJy/beam. (d) RGB J0254$+$3931, VLA A$+$C. Image rms is 0.07
mJy/beam. Image peak is 167.5 mJy/beam. In all images contours start
at 3 times the rms and positive values are spaced by factors of
$\sqrt{2}$ up to the image peak.}
\end{figure*}

\setcounter{figure}{0}

\begin{figure*}
\centerline{
\includegraphics[scale=0.43]{f1e.eps}
\includegraphics[scale=0.43]{f1f.eps}
}
\centerline{
\includegraphics[scale=0.43]{f1g.eps}
\includegraphics[scale=0.43]{f1h.eps}
}
\caption{\small (e) RGB J1629$+$4008, VLA A. Image rms is 0.04
mJy/beam. Image peak is 8.6 mJy/beam. (f) RGB J2229$+$3057, VLA
A. Image rms is 0.04 mJy/beam. Image peak is 44.2 mJy/beam. (g) RGB
J2229$+$3057, VLA A$+$C. Image rms is 0.04 mJy/beam. Image peak is
44.2 mJy/beam. (h) RGB J2256$+$2618, VLA C. Image rms is 0.06
mJy/beam. Image peak is 24.7 mJy/beam. In all images contours start at
3 times the rms and positive values are spaced by factors of
$\sqrt{2}$ up to the image peak.}
\end{figure*}

\setcounter{figure}{0}

\begin{figure*}
\centerline{
\includegraphics[scale=0.43]{f1i.eps}
\includegraphics[scale=0.43]{f1j.eps}
}
\centerline{
\includegraphics[scale=0.43]{f1k.eps}
\includegraphics[scale=0.43]{f1l.eps}
}
\caption{\small (i) RGB J2308$+$2008, VLA A. Image rms is 0.06
mJy/beam. Image peak is 218.0 mJy/beam. (j) RGB J2318$+$3048, VLA
C. Image rms is 0.06 mJy/beam. Image peak is 22.2 mJy/beam. (k) WGA
J0106$-$1034, VLA A. Image rms is 0.03 mJy/beam. Image peak is 105.5
mJy/beam. (l) WGA J0110$-$1647, VLA A. Image rms is 0.10
mJy/beam. Image peak is 72.0 mJy/beam. In all images contours start at
3 times the rms and positive values are spaced by factors of
$\sqrt{2}$ up to the image peak.}
\end{figure*}

\setcounter{figure}{0}

\begin{figure*}
\centerline{
\includegraphics[scale=0.43]{f1m.eps}
\includegraphics[scale=0.44]{f1n.eps}
}
\centerline{
\includegraphics[scale=0.43]{f1o.eps}
\includegraphics[scale=0.43]{f1p.eps}
}
\caption{\small (m) WGA J0126$-$0500, VLA A. Image rms is 0.13
mJy/beam. Image peak is 14.6 mJy/beam. (n) WGA J0227$-$0847, VLA
A. Image rms is 0.10 mJy/beam. Image peak is 76.7 mJy/beam. (o) WGA
J0259$+$1926, VLA A. Image rms is 0.05 mJy/beam. Image peak is 133.1
mJy/beam. (p) WGA J0304$+$0002, VLA A. Image rms is 0.03
mJy/beam. Image peak is 8.4 mJy/beam. In all images contours start at
3 times the rms and positive values are spaced by factors of
$\sqrt{2}$ up to the image peak.}
\end{figure*}

\setcounter{figure}{0}

\begin{figure*}
\centerline{
\includegraphics[scale=0.43]{f1q.eps}
\includegraphics[scale=0.43]{f1r.eps}
}
\centerline{
\includegraphics[scale=0.43]{f1s.eps}
\includegraphics[scale=0.43]{f1t.eps}
}
\caption{\small (q) WGA J0304$+$0002, VLA A$+$C. Image rms is 0.04
mJy/beam. Image peak is 29.1 mJy/beam. (r) WGA J0435$-$0811, VLA
A. Image rms is 0.04 mJy/beam. Image peak is 30.3 mJy/beam. (s) WGA
J0447$-$0322, VLA A. Image rms is 0.09 mJy/beam. Image peak is 68.0
mJy/beam. (t) WGA J0544$-$2241, VLA A. Image rms is 0.06
mJy/beam. Image peak is 149.6 mJy/beam. In all images contours start
at 3 times the rms and positive values are spaced by factors of
$\sqrt{2}$ up to the image peak.}
\end{figure*}

\setcounter{figure}{0}

\begin{figure*}
\centerline{
\includegraphics[scale=0.43]{f1u.eps}
\includegraphics[scale=0.43]{f1v.eps}
}
\centerline{
\includegraphics[scale=0.43]{f1w.eps}
\includegraphics[scale=0.43]{f1x.eps}
}
\caption{\small (u) WGA J1026$+$6746, VLA A. Image rms is 0.12
mJy/beam. Image peak is 84.7 mJy/beam. (v) WGA J1457$-$2818, VLA
A. Image rms is 0.16 mJy/beam. Image peak is 190.4 mJy/beam. (w) WGA
J2239$-$0631, VLA C. Image rms is 0.07 mJy/beam. Image peak is 52.6
mJy/beam. (x) WGA J2320$+$0032, VLA A. Image rms is 0.15
mJy/beam. Image peak is 65.6 mJy/beam. In all images contours start at
3 times the rms and positive values are spaced by factors of
$\sqrt{2}$ up to the image peak.}
\end{figure*}

\setcounter{figure}{0}

\begin{figure*}
\centerline{
\includegraphics[scale=0.43]{f1y.eps}
\includegraphics[scale=0.43]{f1z.eps}
}
\centerline{
\includegraphics[scale=0.43]{f1aa.eps}
\includegraphics[scale=0.43]{f1ab.eps}
}
\caption{\small (y) WGA J2322$+$2114, VLA A. Image rms is 0.04
mJy/beam. Image peak is 40.4 mJy/beam. (z) WGA J2347$+$0852, VLA
A. Image rms is 0.06 mJy/beam. Image peak is 19.4 mJy/beam. (aa) WGA
J2347$+$0852, VLA A$+$C. Image rms is 0.08 mJy/beam. Image peak is
52.3 mJy/beam. (ab) PKS 0256$-$005, VLA A. Image rms is 0.05
mJy/beam. Image peak is 252.1 mJy/beam. In all images contours start
at 3 times the rms and positive values are spaced by factors of
$\sqrt{2}$ up to the image peak.}
\end{figure*}

\setcounter{figure}{0}

\begin{figure*}
\centerline{
\includegraphics[scale=0.43]{f1ac.eps}
}
\caption{\small (ac) 0959$+$68W1, VLA A. Image rms is 0.06
mJy/beam. Image peak is 80.1 mJy/beam. Contours start at 3 times the
rms and positive values are spaced by factors of $\sqrt{2}$ up to the
image peak.}
\end{figure*}

\clearpage

\begin{figure}
\plotone{f2.eps}
\caption{
\label{logRlext}
The extended radio power versus the radio core dominance parameter $R$
at 1.4 GHz, where $R=L_{\rm core}/L_{\rm ext}$, with $L_{\rm core}$
and $L_{\rm ext}$ the radio core and extended luminosities,
respectively. Filled triangles and squares indicate steep-spectrum
(SSRQ) and flat-spectrum radio quasars (FSRQ), respectively, open
squares indicate BL Lacs. See text for details on the selected
samples. Large filled circles represent the sources in our
sample. Arrows indicate limits. The horizontal solid lines mark the
regions of extended radio powers typical of FR I and FR II radio
galaxies.}
\end{figure}

\begin{figure}
\plotone{f3.eps}
\caption{
\label{logRLLS}
The largest linear size (LLS) versus the radio core dominance
parameter $R$ at 1.4 GHz. Symbols are as in Fig. \ref{logRlext}. The
dotted and dashed lines represent LLS $-$ $R$ relations predicted by
the relativistic beaming model for FR I and FR II radio galaxies,
respectively, assuming different intrinsic largest linear sizes. See
text for details.}
\end{figure}

\begin{figure}
\plotone{f4.eps}
\caption{
\label{logRalr}
The radio spectral index versus the radio core dominance parameter $R$
at 1.4 GHz. Symbols are as in Fig. \ref{logRlext}. The vertical solid
line separates lobe- (left of the line) from core-dominated sources
(right of the line). The horizontal solid line indicates the locus of
constant $\alpha_{\rm r}=0.5$. Radio-loud AGN below this line are
currently defined as blazars. The dashed lines represent $\alpha_{\rm
r} - R$ relations for different values of the core radio spectral
index of 0.5, 0, $-0.5$, and $-1.3$ (from top to bottom).}
\end{figure}

\begin{figure}
\plotone{f5.eps}
\caption{\label{lrclx}
The ratio of radio core emission at 1.4 GHz to total X-ray luminosity
at 1 keV versus the radio core dominance parameter $R$ at 1.4
GHz. Symbols are as in Fig. \ref{logRlext}. The vertical solid line
separates lobe- (left of the line) from core-dominated sources (right
of the line). The horizontal solid line separates sources with jet
X-ray emission dominated by the synchrotron (above the line) and
inverse Compton process (below the line). The dashed lines represent
significant correlations for our lobe-dominated sources, SSRQ with
detected X-rays, and all SSRQ (from top to bottom).}
\end{figure}

\end{document}


Keck Infrared Observations of GRO J0422+32 in Quiescence

Abstract: We present Keck K-band photometry and low-resolution H & K-band spectroscopy of the X-ray nova GRO J0422+32 obtained while the system was in the quiescent state. No clear ellipsoidal modulation is present in the light curve, which is instead dominated by a strong flickering component. In the K-band we observe strong Br_gamma emission, with an equivalent width of 38 +/- 5 Angstroms. From this we conclude that the accretion disc is the most likely source of the observed photometric contamination, and that previous infrared-based attempts to constrain the mass of the putative black hole in this system are prone to considerable uncertainty. We finally proceed to show how it is possible to place meaningful constraints on some of the binary parameters of this system, even in the presence of a relatively high level of contamination from the disc.

\section{Introduction}

Low-mass X-ray binaries (LMXBs) are systems in which a degenerate
primary, typically a neutron star or black hole, accretes from a
low-mass (M$_2$ $ < $ 1 M$_{\sun}$), late-type companion (secondary) star.
X-ray novae (XRNe) form a subset of the LMXBs: these are transient
systems which undergo periods of greatly enhanced emission, by
factors of as much as $\sim 10^6$ or greater, before returning to
quiescence over a timescale of months. However, for many of these
systems, as they return to quiescence, the secondary (typically of K
or M type) begins to make a significant contribution to the optical
flux of the system. Observations of the secondary in this state are
the most reliable means for determining the mass of the compact
object. Thus far, 20 systems have been dynamically confirmed as having
a probable black hole primary; of these, 3 are persistent sources
while the other 17 are XRNe \citep{b37}; see \citet{b5} for a thorough
review of XRNe and \citet{b2} for a more recent review of black hole binaries.

GRO J0422+32 was first detected by the \textit{Compton Gamma
Ray Observatory} (GRO) in August 1992 by \citet{b3}. The optical
counterpart was identified soon thereafter and found to have
V $\approx$ 13 mag \citep{b4}. The retreat to quiescence was
punctuated by 3 further mini-outbursts of approximately 4 magnitudes each
(\citealp{b6,b7,b8}). By September 1994, the system was deemed to have
returned to quiescence \citep{b10}, with V $\approx$ 22.4 mag, over 2 years
after the initial outburst. In total the luminosity increased by over
9 magnitudes in V, making GRO J0422+32 the highest-amplitude XRN observed
to date. \citet{b4} searched archival plates and found no
indication of any previous outburst from this system since at least 1923.

Previous attempts to determine the mass of the compact object in this
system have been made by \citet{b18}, \citet{b19}, and \citet{b12};
through optical observations, they determined the mass of the primary
to be M$_x$ = 3.57 $\pm$ 0.34~M$_{\sun}$, 2.5$M_{\sun} < $M$_x < 5 M_{\sun}$,
and M$_x > 2.2 M_{\sun}$, respectively. Recent work in the infrared (IR)
by \citet{b13} has led to a value of M$_x$ = 3.97 $\pm$ 0.95~ M$_{\sun}$;
however, negligible contamination arising from the accretion disc was
assumed. The observations we describe here were undertaken in
an effort to constrain the IR contamination due to the accretion disc,
and hence to determine the reliability of the current mass estimate
for this black hole candidate.

\section{Data}

IR photometry and spectroscopy were undertaken using the Near Infrared
Camera (NIRC; \citealt{b1}) at the f/25 forward Cassegrain port of
the Keck-I 10-m telescope. NIRC is a 256 $\times$ 256 pixel InSb array
with an angular scale of 0.15${\arcsec}$ pixel$^{-1}$. Low-resolution
spectra were obtained using a grism (R $\approx$ 120) and the HK
filter, which covers the standard H and K-bands.

\subsection{Photometry}

Photometry of GRO J0422+32 was undertaken on the night of 1997
November 9 UT. The system was observed in the K-band for a
full orbital cycle of $\sim$5.1 hr. Conditions were excellent
throughout the night, with seeing no worse than 0.45${\arcsec}$ at
any time. Individual exposure times were 4~s, with 30 coadds per
image (i.e., 2~min exposure per image). The telescope pointing was
dithered (using a 3 $\times$ 3 grid) to allow for accurate
subtraction of the IR background; see Table \ref{obs-table} for
a log of observations. A single 9-image grid of observations of
similar duration was also taken nearly a year later (1998 September
28 UT), contemporaneously with the spectroscopy (see below).

\begin{figure}
\begin{center}
\includegraphics[height=72mm,width=84mm]{fig1}
\caption{The NIRC GRO J0422+32 field. The black arrow indicates
the position of GRO J0422+32. The exposure time was 2 min.}
\label{chart}
\end{center}
\end{figure}

\subsubsection{Data Reduction}

The data were dark-current subtracted, flat fielded, background
subtracted, and combined using standard IRAF routines\footnote{IRAF
is distributed by the National Optical Astronomy Observatories,
which are operated by the Association of Universities for Research
in Astronomy, Inc., under cooperative agreement with the National
Science Foundation}. Photometry of GRO J0422+32 and a number of
nearby stars was carried out using the profile-fitting task DAOPHOT
\citep{b11}. The standard stars FS12 and
FS16\footnote{www2.keck.hawaii.edu/inst/nirc/UKIRTstds.html}
\citep{b39} were also
observed to allow accurate calibration of the target frames. The
magnitudes were cross-checked using the bright star to the north
(2MASS J042142.57+325449.5), which is visible in the
2MASS\footnote{The Two Micron All Sky Survey is a joint project
of the University of Massachusetts and the Infrared Processing
and Analysis Center, California Institute of Technology, funded
by NASA and the National Science Foundation.} field of GRO J0422+32.
We display one of our frames in Fig. \ref{chart}. The resulting
light curve was then phased to the ephemeris of \citet{b12}. As our
data was taken $\sim$ 125 days after the ephemeris of \citet{b12}, any
errors will be minimal.

\begin{figure}
\begin{center}
\includegraphics[height=84mm,width=84mm,angle=-90]{fig2}
\caption{The NIRC K-band light curve of GRO J0422+32. The top
light curve is the 2-min resolution K-band light curve. We
plot this binned to a resolution of 18~min in the middle. As a
comparison, the best-fit model of \citet{b13} is plotted below.
The unbinned light curve is accurately calibrated, while the
middle and bottom light curves are shifted downward by 1 and 2 mag
(respectively) for clarity.}
\label{alllc}
\end{center}
\end{figure}

We find that our magnitudes are in agreement with those of
\citet{b13}, confirming that our observations were indeed taken in
quiescence. It is apparent from our images that there is a
line-of-sight star to the northwest, lying approximately 2.6$\arcsec$ from
GRO J0422+32, which may have been unresolved in previous observations of
the system (see Fig. \ref{chart}). Profile fitting was carried out on
this star simultaneously with our target star, GRO J0422+32. We find
this star to have a mean magnitude of 19.62 $\pm$ 0.06, which was
constant over the course of our observations.

The photometric accuracy of our GRO J0422+32 light curve (2min
resolution) was determined
by applying a linear fit to the light curve of the star to the northeast
(Fig. \ref{chart}), of comparable magnitude to GRO J0422+32, which was
assumed to be constant during the course of our observations. The
uncertainty was then taken to be the root-mean-square (rms) deviation
from this linear fit and was measured to be 0.024 magnitudes.

The 1998 observation was reduced in a manner identical
to that described above. When phased and calibrated similarly to the
previous epochs data, the magnitude was found to be consistent. This
confirms that the system was in a quiescent state similar to that
detected the previous year.

\begin{table*}
\caption{NIRC Observation Log.}
\label{obs-table}
\begin{center}
\begin{tabular}{ccccc}
\hline
JD$^a$ & Seeing [$\arcsec$] & Exp. time [s] & $\phi^b$ & Apparent
K-band mag$^c$ \\
\hline\hline
2450762.91995413 & 0.40$\arcsec$ & 1080 &  0.528 & 17.344 $\pm$ 0.014\\
2450762.93619041 & 0.35$\arcsec$ & 1080 &  0.605 & 17.389 $\pm$ 0.010\\
2450762.95500179 & 0.35$\arcsec$ & 1080 &  0.693 & 17.364 $\pm$ 0.015\\
2450762.97392475 & 0.30$\arcsec$ & 1080 &  0.783 & 17.397 $\pm$ 0.014\\
2450762.98945141 & 0.30$\arcsec$ & 1080 &  0.856 & 17.452 $\pm$ 0.015\\
2450762.00493642 & 0.30$\arcsec$ & 1080 &  0.929 & 17.406 $\pm$ 0.015\\
2450763.02045906 & 0.30$\arcsec$ & 1080 &  0.002 & 17.368 $\pm$ 0.015\\
2450763.03688440 & 0.30$\arcsec$ & 1080 &  0.079 & 17.311 $\pm$ 0.015\\
2450763.05382848 & 0.30$\arcsec$ & 1080 &  0.159 & 17.282 $\pm$ 0.015\\
2450763.06980607 & 0.35$\arcsec$ & 1080 &  0.234 & 17.345 $\pm$ 0.014\\
2450763.08537309 & 0.40$\arcsec$ & 1080 &  0.309 & 17.345 $\pm$ 0.015\\
2450763.10399760 & 0.40$\arcsec$ & 1080 &  0.392 & 17.389 $\pm$ 0.015\\
2450763.11857811 & 0.45$\arcsec$ & 1080 &  0.464 & 17.298 $\pm$ 0.023\\
2451085.99342259 & 0.35$\arcsec$ & 1080 &  0.310 & 17.264 $\pm$ 0.093\\
\hline
\end{tabular}
\end{center}
\medskip
{$^{a}$Exposure time consists of nine 2-min exposures; JD is at the
  midpoint. \\
$^{b} T_0$ = 2450274.4156, $P$ = 0.2121600 days \\
$^{c}$Magnitudes correspond to those of the binned light curve; see
  Fig. \ref{alllc}.}
\end{table*}

\subsection{Spectroscopy}
We obtained a number of low-resolution NIRC spectra of
GRO J0422+32 on the night of 1998 September 28 UT. The HK
filter was used (1.40--2.53~$\mu$m) in conjunction with
the gr120 grism and a 0.67$\arcsec$ slit, giving a resolution
R $\approx$ 120. Observing conditions were once again excellent,
with seeing as good as 0.30$\arcsec$, for the majority of the
observations. The system was observed from orbital phase 0.36 to
0.78. Individual exposure times were 200~s, dithered to 5
positions along the slit. In total 30 spectra of the system
were acquired.

\subsubsection{Data reduction}
The one-dimensional spectra were extracted in the standard manner
using IRAF. Of the 30 exposures, 28 contained useful spectra.
Spectra of the G0V stars BS8455 and BS1789 were also taken to
aid with the removal of telluric features. No arc spectra were
taken due to the availability of a known calibration, which was
obtained from the NIRC online
manual\footnote{http://www2.keck.hawaii.edu/inst/nirc/}.

The GRO J0422+32 and G0V star spectra were normalised using a
low-order spline. The region surrounding the Br$\gamma$ feature
in the G0V spectrum was masked and the GRO J0422+32 spectra were
then divided by the G0V spectra to remove the telluric features.
The resulting spectra were median combined so as to negate the effect
of cosmic-ray hits and other spurious effects, while at the same time
maximising the available signal-to-noise ratio (S/N).

No significant features were detected in the resulting H-band
spectrum. In the K-band we observe strong Br${\gamma}$ emission,
which we display in Fig. \ref{kspec}. The EW of the Br${\gamma}$
feature was measured to be 38 $\pm$ 5 {\AA}.

\begin{figure}
\includegraphics[height=84mm,width=64mm,angle=-90]{fig3}
\caption{The NIRC $K$-band spectrum of GRO J0422+32.}
\label{kspec}
\end{figure}

\section{Analysis}
\subsection{Colour}
We calculated the H-K colour index of the system from our
spectra. This was achieved by fitting a polynomial to the continuum of
the individual H\&K-band spectra and then integrating the flux across
the respective passband. The median colour was found to be H-K  -0.02
$\pm$ 0.01 mag, where the error is the rms deviation of the
data. Correcting for interstellar extinction using a value of E(B-V) =
0.3 $\pm$ 0.1 mag (\citealt{b36}) for the colour excess leads to an
extinction-corrected colour of (H-K)$_0$ = 0.07 $\pm$ 0.08 mag,
which equates to spectral types in the range $\sim$ B5 to $\sim$ K7
\citep{b20}. This is somewhat earlier than the expected spectral type
of $\sim$ M2 estimated by (\citealt{b18}), although the uncertainties
are large.

\subsection{Power Density Spectrum}
It is immediately apparent in the 2-minute resolution light curve
(Fig. \ref{alllc}) that  any ellipsoidal modulation in the data is
masked by a strong flickering component. In an attempt to decouple the
ellipsoidal variation of the secondary star from the flickering
component, we followed the method of \citet{b23}. As the data sampling
was not uniform a Lomb-Scargle periodogram \citep{b24} of the 2 minute
resolution data was calculated. A significant signal ($\sim$ 5$\sigma$) is
detected at a period $\sim$ 2.5 hrs, $\sim$ half the orbital period.

The amplitude of the modulation at this period is only 0.03
magnitudes. However, as can be clearly observed from the binned (18
minute resolution data: Fig.\ref{alllc}), this variation cannot be
ellipsoidal in origin: the minima for this modulation occur at phases
$\sim$ 0.3 \& 0.8 and not at phases 0.0 \& 0.5 as would be expected
for an ellipsoidal modulation. This phase shift is much greater than
the error in calculating the orbital phase from the known ephemeris
of \citet{b12}. Hence, we can use our data to place a conservative
upper limit of 0.03 magnitudes on the amplitude of any ellipsoidal
variation present in our data.

In an attempt to characterise the shorter timescale variability
present in our light curve, we computed the Power Density Spectrum
(PDS). A power law fit to our data of the form $f(x) \propto x^\beta$
yielded a power law index $\beta \approx$ -0.9: see Fig. \ref{logpsd}.
This is consistent with the variability timescales observed in the
optical in both A0620-00 and GRO J0422+32 by \citet{b23}, who measured
significant variability on timescales of up to 90 minutes.

\begin{figure}
\includegraphics[height=84mm,width=64mm,angle=-90]{fig4}
\caption{The IR power spectrum of GRO J0422+32. The solid line
  indicates the best power law fit to the data.}
\label{logpsd}
\end{figure}

\section{Discussion}
In the JHK photometry of \citet{b13}, an ellipsoidal modulation is
present in both sets of K-band observations (from 2000 and
2003). By contrast, our unbinned light curve (resolution $\sim$2 min),
while showing rapid variability, appears to display little evidence
of an ellipsoidal modulation. Binning the data to
$\sim$18 min resolution (this is equal to the duration of our 3
$\times$ 3 dither pattern and comparable to that of \citealt{b13}),
gives rise to a modulation that, although superficially ellipsoidal in
nature, clearly is not (see section 3.2). \citet{b14} observed this
system in 2002 with the 3.5m WIYN telescope, they also report a
modulation on a $\sim$ 5 hour period that is not ellipsoidal in nature.
It is unclear why our results differ from those of \citet{b13},
especially as the mean K-band magnitude we measure is consistent with
theirs (and hence there cannot have been a significant change in the
overall contribution of the flux from the disc to the light curve) but
it may be that their K-band modulation is similar to the one we
observe, but coincidentally at a phase that made it appear ellipsoidal
in origin.

In any case we now discuss the likely sources of the lightcurve we observe:
$(1)$ spots on the M-type secondary star, $(2)$ contamination originating in
a radio jet and/or  $(3)$ contamination from the quiescent accretion disc.

\subsection{Spots on the M-type Secondary}
Deviations from a pure ellipsoidal modulation have sometimes been
attributed to the presence of spots on the rapidly rotating
secondary star \citep{b61,b58,b57}.
We searched the literature in an attempt to find examples of M-type
dwarfs containing spots, but the number of candidate systems is few,
given that we require a rapidly rotating companion: $v_{rot} \sin(i)$
for GRO J0422+32 is close to 90 km s$^{-1}$ \citep{b50}.
The two best candidates are the isolated M1 stars HK Aqr and RE
1816+541, with rotational velocities of $\sim$ 70 km s$^{-1}$ and
$\sim$ 60 km s$^{-1}$, respectively.

Detailed V-band light curves of HK Aqr \citep{b27} show a $\sim$ 0.09
mag modulation due to spots. However, this modulation appears to be quite
uniform and smooth, in contrast to the rapid variability which we
observe in GRO J0422+32. Indeed, recent Doppler images \citep{b28}
reveal the presence of a large number of star spots on the surface
of these rapidly rotating M-type stars even though the observed variability
is quite low. Hence, we do not believe that star spots
offer a satisfactory explanation for the observed rapid variability of
GRO J0422+32 in the IR.

\subsection{Contamination from a Jet}
In recent years evidence has emerged that the synchrotron flux from a
radio jet could be contributing to the flux emitted at other
wavelengths, in particular the IR and optical \citep{b41}. IR
contamination, which has been attributed to a jet, has been observed
in the halo black hole candidate XTE J1118+480 during outburst
\citep{b53}. \citet{b43} claim that all black hole XRBs will become
jet dominated in the quiescent state. To date emission from a radio
jet in quiescence has been detected for two XRBs, V404 Cyg (GS
2023+338, \citealt{b54}) and A0620-00 \citep{b55}. The observation of
A0620 is the more pertinent here due to the similarity of its orbital
period (7.8 hrs) and quiescent X-ray luminosity ($3 \times 10^{30}$
erg s$^{-1}$) compared to those of GRO J0422+32 (5.1hrs and $8 \times
10^{30}$ erg s$^{-1}$ respectively). However, recent \textit{Spitzer}
observations \citep{b56} do not appear to be consistent with emission
from the jet in A0620-00 longwards of 8 microns: shortwards of 8
microns \citet{b56} argue that the emission appears to be consistent
with that originating from an optically thick blackbody and not the
photosphere of the secondary star. Hence, by analogy with A0620-00, we
think it unlikely that the flux we observe in the case of GRO J0422+32
originates in a jet.

\subsection{Contamination from the Accretion Disc}

\subsubsection{Photometric Constraints}
Previous estimates of the fractional contribution of the accretion
disc to the total flux range from $\sim$60\% \citep{b18} to as little
as 20\% \citep{b33} in the R-band. In a subsequent re-analysis of the
data from \citet{b18}, \citet{b50} measured the flux from the
accretion disc to contribute $\sim$ 40\% of the total optical flux in
quiescence. Similarly \citet{b12} measured the contamination from the
accretion disc to be $\sim$ 60\% in the I-band. These measurements
hint at an increasing level of contamination at longer wavelengths,
although the different template stars used in the respective analyses
(M2 vs M4) limit the conclusions that can be made.

In their analysis of GRO J0422+32, \citet{b13} assumed a negligible
contribution from the accretion disc in the K-band.
To check this assumption, we have used the ELC light curve modelling
code \citep{b21} to model the flux from an accretion
disc, assuming a temperature profile of the form $T(r) \propto
T_{inner}(r/r_{inner})^{\xi}$, with $\xi$ chosen so as to ensure a
steady state disc, and an outer disc radius of $\sim$ 0.55 R$_L$
(these values being consistent with observations of the quiescent disc
in the XRNe A0620-00 and XTE J1118+480 (\citealt{b25,b59,b60}).
We calculate that for any value of the accretion disc contribution to
the total flux in the R-band, the fractional contribution of the
accretion disc flux to the K-band is similar, if not greater.
This is in contrast to the canonical picture where the flux from the
disc falls off as we venture towards the IR. The
contributions from the accretion disc measured in the R \& I-band in
the case of GRO J0422+32 agree with this simple model and hence we
would expect a similar level of contamination in the K-band.  This
result is dependent on the inclination of the system, as the flux from
the accretion disc relative to that from the secondary increases with
lower inclination.

In Fig. \ref{rel-flux}, we plot the VRJHK spectral energy distribution
(SED) for various M-type secondary stars along with that of GRO
J0422+32. A distance of 2.6 $\pm$ 0.2 kpc and a reddening of E(B-V) =
0.3 $\pm$ 0.1 have been assumed (from \citet{b51}). The absolute
V-band magnitudes and the various colour relations of the M-type stars
are taken from \citet{b52}. The dashed line indicates the observed SED
of GRO J0422+32. The error bars account for both the uncertainty in
the reddening and the distance values quoted above. The observed
colours are broadly consistent with a spectral type of M1, and suggest
a disk contribution in the K-band of no more that $\sim$ 0.3 mag
(albeit at the 1$\sigma$ level). This value would be consistent with
that observed in the R-band, as expected from the modelling discussed
above. Alternatively, if there existed a
cool optically thick component to the accretion disc (as might be
expected theoretically - see \citealt{b63}), a contribution from a radio
jet (see Section 4.2 above), or some form of circumbinary disc
(\citealt{b56}), any of these contributions could lead to a flatter SED
than might otherwise be expected.

A spectral type of M1 is consistent with previous spectral type
determinations (\citealt{b19,b18,b13}), with the exception of
\citet{b12} who measured a spectral type between M4 \& M5.  To
reconcile the latter spectral type with the observed SED would require
that the distance estimate used above is significantly in error, or a
K-band contamination from the accretion disk of $\sim$ 3 magnitudes
(which we also regard as unlikely).

\subsubsection{Spectroscopic Measurements}

The Br${\gamma}$ feature at 2.16$\mu$m, with a equivalent width of
$\sim$ 38 \AA, dominates the K-band spectrum of GRO J0422+32.  No
other spectral lines are visible in either the H or K-bands. We also
note the absence of the CO bandheads from the companion star in this
spectrum, although it is possible that they are simply unresolved.

Previous observations of LMXBs support the conclusion that the
presence of strong Br${\gamma}$ emission is due to contamination by
the accretion disc.  IR spectroscopy of a sample of persistently
bright neutron star LMXB \citep{b46,b15,b49} revealed the presence of
Br${\gamma}$ emission. In the sources with the strongest emission (Sco
X-1, Sco X-2, GX 5-1), the presence of a large Br${\gamma}$ feature
(EW $\sim$ 23, 29, 45 \AA) was accompanied by the absence of the CO
bandheads or any other spectral features due to the secondary
star. In contrast, in the sources with lower levels of Br${\gamma}$
emission (GX 1+4, GX 13+1, EW $\sim$ 10\AA, 5\AA), the CO bandheads
and other lines due to the secondary star were detected.

Previous quiescent IR observation of the XRNe V404
Cyg and A0620-00 \citep{b47,b48} also reveal the presence of
Br${\gamma}$ emission. In V404 Cyg the equivalent width of the
Br${\gamma}$ line was measured to be only 2.9 \AA. In this case the CO
bandheads from the secondary star were also detected. The
Br${\gamma}$ line was detected with an EW $\sim$ 15 \AA~ in
A0620-00. In this case the CO bandheads were only marginally
detected (the S/N of the spectrum was very low). Using these
spectra the authors measured the IR contamination from the accretion
disc to be $\leq$ 14\% and 27\% respectively. Recent IR spectroscopy of
A0620-00 by \citet{b40} support the measurement of excess K-band
flux in this system.

In summary, our data, both spectroscopic and photometric, strongly
suggest a significant contribution to the K-band flux of GRO J0422+32
from the accretion disc. As previously discussed, these observations
are not the first to observe ``flickering'' in the light curve of a
quiescent XRN: \citet{b31} obtained fast photometry of the XRNe
A0620--00, Nova Mus 1991, and MM Vel in the optical and detected rapid
variability similar to that reported here. \citet{b23} obtained fast
photometry of 5 XRNe including GRO J0422+32; significant variability
was found, superposed on the ellipsoidal variations. We note that of
the 5 XRNe observed optically (V404 Cyg, A0620-00, GRO J0422+32, GS
2000+25 \& Cen X-4), the largest amplitude variability was observed in
the light curve of GRO J0422+32. \citet{b23} contend that the observed
optical flaring originates in the accretion disc and not as
chromospheric activity from the secondary. Our observations are the
first to show this kind of flickering variability in the K-band, where
such contamination had previously been thought to be minimal for
quiescent XRNe.

\begin{figure}
\includegraphics[height=84mm,width=64mm,angle=-90]{fig5}
\caption{VRJHK SED of various M-type stars and GRO J0422+32 (dashed line).}
\label{rel-flux}
\end{figure}

\subsection{Mass of the Compact Object}
As demonstrated above, there is no clear ellipsoidal modulation
present in our K-band data. However, for a given accretion disc
contamination level in the K-band, we can still hope to constrain
the mass ratio $q$, and orbital inclination, $i$, as follows. Firstly,
we assume that the true ellipsoidal modulation has an amplitude less
than the 0.03 magnitude variation discussed in Section 3.2.

Assuming that the excess flux from the accretion disc contributes
approximately 30\% of the flux in the K-band (see section 4.3.1), an
intrinsic ellipsoidal modulation of 0.04 magnitudes would be reduced to
the 0.03 magnitude variation limit above. In Fig. \ref{mfunc}, we
display the $q$ vs $i$  plot for GRO J0422+32: the dotted line
indicates the best fit to a 0.04 magnitude modulation. We can
represent the ellipsoidal modulation by a single amplitude because, at
these relatively low inclinations, the depth of the minima at phases
0.0 and 0.5 are very similar.

As this line represents an upper limit to any ellipsoidal modulation
present, allowed values of $q$ and $i$ are constrained to lie to the
right of and below this line. By combining this with the limit from
the radial velocity measurements (assuming a secondary mass of 0.45
M$_{\sun}$), as delineated by the solid line, we can constrain $q$ and
$i$ as follows: $i < 30 \degr$ and $q > 23$. Note that the inclination
constraint is an upper limit because we have only an upper limit on
the amplitude of the ellipsoidal modulation. The latter
yields a lower limit to the compact object mass of $\sim$ 10.4
M$_{\sun}$. This is in contrast to the inclination of 45$\degr \pm$
2$\degr$ (M$_x$ = 3.97 $\pm$ 0.95 M$_{\sun}$) obtained by \citet{b13},
bearing in mind the uncertainties discussed earlier.

\begin{figure}
\includegraphics[height=84mm,width=64mm,angle=-90]{fig6}
\caption{Mass ratio inclination plot for GRO J0422+32. The dashed
  lines indicate lines of constant primary mass. The solid lines
  indicates lines of constant mass appropriate for an M0, M1 \& M2-type
  secondaries. The dotted line indicate the best fit q,i pairs for a
  0.04 amplitude sine wave. }
\label{mfunc}
\end{figure}

\section{Conclusions}

The K-band properties of GRO J0422+32 are consistent with a
flux contaminated by the emission from the accretion disc, and not
solely from the secondary as previously supposed.

Previous authors have assumed negligible contamination of the IR flux
of quiescent XRNe by the accretion disc, and hence that the
ellipsoidal modulation and black hole mass could be most accurately
measured in the IR (in comparison to the optical). This conclusion was
generally based on an extrapolation of the contamination from the
optical to longer wavelengths (see \citealt{b22}). Our data show that
this is apparently an invalid assumption in the case of GRO J0422+32.

Our observations are the first to detect flickering from a quiescent
XRN accretion disc in the K-band, most likely because of the superior
S/N obtainable with the 10-m Keck telescope compared with the smaller
telescopes from which similar measurements have been previously
obtained. However. it is also possible that the K-band flickering is
present because GRO J0422+32 is a relatively low inclination system:
in this case, if the source of the flickering is located in the inner
disc, then the flickering would be more easily observable in a low
inclination system such as GRO J0422+32 in comparison to other (higher
inclination) quiescent XRNe.

In light of the observations discussed in this paper, it is clear that
higher time resolution photometry or higher S/N spectroscopic
observations of other XRNe in the IR are required to better constrain
the contribution of the accretion disc at these wavelengths, and to
determine the degree to which IR ellipsoidal variability measurements can be
used to reliably constrain the mass of black holes in quiescent
XRNe. Flickering of the amplitude reported here (0.2 magnitudes) would
have a significant effect on even the highest amplitude ellipsoidal
variation ($\sim$ 0.3 magnitudes: see for example the case of GRO
J1655-40 discussed by \citealt{b62}). Even when such contamination is
significant, however, it may still be possible to extract usefull
limits on system parameters, as long as the degree of such
contamination can be quantified (e.g. spectroscopically), as we have
discussed here.

\bigskip

This research made use of the SIMBAD database, operated at CDS,
Strasbourg, France, and NASA's Astrophysics Data System.
The data presented herein were obtained at the W. M. Keck
Observatory, which is operated as a scientific partnership
among the California Institute of Technology, the University
of California, and the National Aeronautics and Space Administration.
The Observatory was made possible by the generous financial support
of the W. M. Keck Foundation. A.V.F. is grateful for support
from NSF grant AST-0307894. M.T.R. \& P.J.C. acknowledge financial support
from Science Foundation Ireland. We thank an anonymous referee for
useful comments which improved the quality of this paper.