A&A 466, 63-73 (2007)
DOI: 10.1051/0004-6361:20066754
Zhongzu Wu1,2,3 - D. R. Jiang1,2 - Minfeng Gu1,2 - Yi Liu1,2
1 - Shanghai Astronomical Observatory, Chinese Academy of Sciences, Shanghai 200030, PR China
2 - Joint Institute for Galaxy and Cosmology (JOINGC) of SHAO
and USTC
3 - Graduate School of the Chinese Academy of Sciences, Beijing 100039, PR China
Received 16 November 2006 / Accepted 27 December 2006
Abstract
We present EVN observations of seven BL Lac objects selected from
the RGB sample. To investigate the intrinsic radiation property of
BL Lac objects, we estimated the Doppler factor with the VLA or MERLIN core and the total 408 MHz luminosity for a sample of 170 BL Lac objects. The intrinsic (comoving) synchrotron peak frequency was then calculated by using the estimated Doppler factor. Assuming a
Lorentz factor of 5, the viewing angle of jets was constrained. The
high-resolution VLBI images of seven sources all show a core-jet
structure. We estimated the proper motions of three sources with the
VLBI archive data, and find that the apparent speed increases with
the distance of components to the core for all of them. In our BL
Lacs sample, the Doppler factor of LBLs is systematically larger
than that of IBLs and HBLs. We find a significant anti-correlation
between the total 408 MHz luminosity and the intrinsic synchrotron
peak frequency. However, the scatter is much larger than for the
blazar sequence. Moreover, we find a significant positive correlation between the viewing angle and the intrinsic synchrotron peak frequency. The BL Lac objects show a continuous distribution on the viewing angle. While LBLs have a smaller viewing angle than that of IBLs and HBLs, IBLs are comparable to HBLs. We conclude that the intrinsic synchrotron peak frequency is not only related to the intrinsic radio power (though with a large scatter), but also to the viewing angle for the present sample.
Key words: BL Lacertae objects: general - galaxies: active - galaxies: jets - galaxies: nuclei - radio continuum: galaxies
The so-called "blazar sequence'' was proposed to link the shape of
the SED and the synchrotron peak frequency to the source luminosity
for blazars, which consists of more luminous flat-spectrum radio
quasars (FSRQs) and BL Lac objects (Fossati et al. 1998; Ghisellini
et al. 1998). The most powerful sources have relatively small
synchrotron peak frequencies, and the least powerful ones have the
highest
values. This anti-correlation can be
theoretically explained by the cooling processes. In more powerful
sources, the energy density is higher and the emitting particles
have a higher probability of losing energy so are subjected to more
cooling, resulting in a lower value for
.
However,
this sequence is based on the absence of high-luminosity HBLs and
low-luminosity LBLs, therefore at least part of this systematic
trend can result from selection effects (Antón & Browne 2005).
Indeed, the evidence of low-power LBLs has been recently discovered
(Padonavi et al. 2003; Caccianiga & Marchã 2004; Antón &
Browne 2005), and the possible discovery of high luminosity HBLs is
also reported in the Sedentary survey (Giommi et al. 2005). In addition,
the high-power-high-
FSRQs were found in the Deep
X-ray Radio Blazar Survey (DXRBS), which is both X-ray and
radio-selected, though they do not reach the extreme
values of HBLs. From all these discoveries, it seems
that the blazar sequence in its simplest form cannot be valid
(Padovani 2006).
Instead of focusing on one or two limited surveys,
Nieppola et al. (2006) have examined the properties of BL Lacs through
constructing the SEDs for a large, heterogeneous sample of BL Lacs
taken from the Veron-Cetty & Veron BL Lac catalogue and visible
from the Metshähovi radio observatory. To our knowledge, this is
the largest BL Lacs sample to revisit the blazar sequence up to now,
though it is heterogeneous. The authors have found an
anti-correlation between the radio power and
,
with
a huge scatter, reaching 5 orders of magnitude in power, and many
outliers were also found, especially in the low-power-low-
region. However, it is well known that the radio power
of BL Lacs is affected by the beaming effect of radio jets, since BL
Lacs are generally believed to have a relativistic jet aligned close
to the line of sight (Urry & Padovani 1995). As the blazar sequence
was originally proposed as the anti-correlation between the
intrinsic luminosity and peak frequency, it nevertheless might be
important to revisit the correlation using the intrinsic radio
power. Moreover, the peak frequency
,
where B is the magnetic field,
the Doppler factor, and
a characteristic
electron energy that is determined by a competition between
accelerating and cooling processes. While blazars as a whole
population suffer the Doppler boosting, the difference in Doppler
factor can exist between various subsets, which implies the
importance of investigating the intrinsic peak frequency after
excluding the Doppler factor. As Kollgaard et al. (1996b) find, the average
angles of radio jets to the line of sight is approximately
for XBLs, but
for RBLs. When the BL Lac
objects show a continuous distribution of synchrotron emission-peak
energies and most recent works focus on the correlation between the
power and peak frequency, it is not clear whether they also show a
continuous distribution for the angle to the line of sight, and HBLs
and LBLs occupy the opposite ends of this distribution.
Apart from investigating the SED of BL Lac objects, the high-resolution radio observation is an important way to explore the radio structure of BL Lacs, from which we can obtain the physical information, such as the viewing angle, Lorentz factor and magnetic field. From the VLBI observations, the different subsets of BL Lacs can be compared based on the compact radio structure. In particular, it enables us to explore the jet orientation for BL Lac objects, because the properties of the parsec scale structure are strongly dependent on jet orientation. Kollgaard et al. (1996a) show that the jets in XBLs fade more quickly than in RBLs. The VLBI observations show that most of LBLs display the rapid superluminal apparent motions (Jorstad et al. 2001), while the TEV blazar sources (most are HBLs) display subluminal or mildly relativistic (Giroletti et al. 2004b,2006; Piner & Edwards 2004). Rector et al. (2003) have observed 15 HBLs and 3 BL Lacs from the RGB sample and find that the HBLs, like most LBLs, show parsec-scale core-jet morphologies with complex kilo-parsec scale morphologies. Moreover, the jets in HBLs are more well-aligned, suggesting that the jets of HBLs are either intrinsically straighter or are seen further off-axis than LBLs. Giroletti et al. (2004a) selected 30 low-redshift BL Lac objects and confirmed that parsec and kilo-parsec scale jets are oriented at the same P.A. in a large fraction of HBLs. The HBLs show less distortion and therefore are expected to be oriented at larger angles than the LBL sources.
Laurent-Muehleisen et al. (1999) present a sample of 127 BL Lac objects from the ROSAT ALL-Sky Survey-Green Bank catalogue (RGB), which exhibits properties intermediate between high- and low-energy-peaked BL Lac objects. In this paper, we present EVN observations of seven BL Lac objects selected from the RGB sample. Our goal is to investigate their compact radio structure. In combination with the previous observations and the archive data, we aim to explore the possible structure variation and measure the proper motion for our sources. Moreover, to explore BL Lac objects as a whole population, we investigated the distribution of the intrinsic radio power, the Doppler factor, and the viewing angle along the intrinsic synchrotron peak frequency for a sample of BL Lac objects selected from Nieppola et al. (2006).
The layout of this paper is as follows. In Sect. 2, the observations
and data reduction for our seven sources are given. The compact
radio structure are investigated in Sect. 3, where the proper motion
are estimated for three sources as well. In Sect. 4, the intrinsic
radio low-frequency luminosity, the Doppler factor, and the viewing
angle are estimated for a sample of BL Lacs selected from
Nieppola et al. (2006), and their distribution with the intrinsic
synchrotron peak frequency is also explored. The discussions are
shown in Sect. 5, while the conclusions are drawn in Sect. 6.
Throughout this work, we adopt the following values for the
cosmological parameters:
,
,
and
,
except an
otherwise stated. The spectral index
is defined as
.
All seven BL Lac objects were observed with EVN at 5 GHz on 22 November 1999. The raw data were obtained with the MKIII VLBI recording system, with an effective bandwidth of 28 MHz, and correlated in Bonn at the Max-Planck-Institute. To study the structure and proper motion of these sources, we collected EVN and VLBA archive data. All seven sources are listed in Table 1: Col. (1) Source name, Col. (2) redshift, Cols. (3) the observational date and (4) the half-power beamwidth (HPBW) of the weighted beam, Col. (5) the noise of the image, Col. (6) the peak flux of the image, Col. (7) VLBI array. All the data were reduced using the NRAO Astronomical Image Processing System (AIPS) package. The system temperatures and gains were obtained from the VLBA website. After the initial reduction (including editing, amplitude calibration, instrumental phase corrections, and antenna-based fringe-fitting), we imported the calibrated data into DIFMAP package (Shepherd et al. 1994, 1995) to make images and fit the sources with a number of discrete circular Gaussian components. The fitting was done directly on the final, self-calibrated visibility data. The modelfit parameters are listed in Table 2.
The VLBI images show core-jet structure for all sources we selected,
and some of them are in good agreement with previous observations.
The MERLIN archive images and available VLA images also show that
these sources posse very straight jets. Model fitting of the radio
structures in each image was performed, using the MODELFIT of
Caltech DIFMAP package (uvweight = 0, -1). The modelfit parameters are
presented in Table 2. To determine the uncertainty of
the component position, we used the Difwrap package (Lovell 2000).
The errors were estimated by perturbing the fit parameters until the
resulting residual maps were unacceptable and taking the most
extreme accepted points. A visual inspection was done to determine
the goodness of the fit, and we used the best-fit total
as
the upper limit for every fit of the components.
In our seven BL Lac objects, previous VLBI observations at 6 cm are available for five of them. Combining our observations with the archive data, we got at least four-epoch observations for three sources: 1727+502, 1133+704, and 1741+196. The weighted linear fits to the component distances from the core as a function of the observing time were used to estimate the proper motion of the components. The results are given in Table 3, in which Col. (1) is the source name, (2) the components label, (3) the proper motion, and (4) the apparent speed. From Table 3, we find that all three sources have superluminal motion, and the apparent speed increases with distance from the core. For 1011+496 and 1424+240, only three-epoch data are available, therefore we do not measure their proper motion and only present their images in this work. For 4C+37.46 and 1542+614, there are no previous VLBI images and archive data available, so we present their first image.
1011+496 - Snapshot VLA observations have found hints of
extended structure around a bright core (Machalski & Condon 1983). This was
confirmed by Kollgaard et al. (1996b), who found diffuse emission extending
around
to the north from their VLA B-configuration
data. Augusto et al. (1998) have modeled this source with three components,
which clearly showed that the jet of this source extends to about a
few hundreds mas of the core with position angle
.
The
MERLIN archive image (see Fig. 1) shows that the jet
extends in nearly the same direction with our EVN and VLBA images
(see Fig. 2), which show a jet in the opposite direction
on the parsec scale. As the jet direction is similar in MERLIN, EVN,
and VLBA images, this source is well-aligned.
![]() |
Figure 1:
The MERLIN archive image of 1011+496 at 1.6 GHz on May 6, 1999. The contour levels are -1, 1, 2, 4, 8, 16, 32, 64... times 0.001742 Jy beam-1 (3![]() |
Open with DEXTER |
![]() |
Figure 2: Multi-epoch EVN and VLBA images at 5 GHz. The minimum contour level is 3 times the rms noise given in Table 1. The contour levels are multiples (-1, 1, 2, 4, 8, 16, 32, 64...) of the minimum contour level and the axes are labeled in milliarcseconds. The top three panels are for 1011+496 at 1997.055, 1998.493, and 1999.893 ( from left to right); the middle for 1424+240 at 1999.696, 1999.893, and 2003.447 ( from left to right); while the bottom for 1542+614 and 4C+37.46 at 1999.893. |
Open with DEXTER |
1133+704 - The images from three epochs (see Fig. 3) all show that the jet extends at a structural
position angle of
,
and the outer components extend
to the northeast around
.
These features are similar to
the parsec scale jet found by Kollgaard et al. (1996a) and are well-aligned
with the kilo-parsec scale halo and short inner jet presented by
Giroletti et al. (2004a) and Giroletti et al. (2006). It can be seen from Fig. 3 that five components are detected in 1133+704, which
is labeled as M0, M1, M2, M3, and M4, respectively. The weighted
linear fits of components M1, M2, and M3 are presented in Fig. 4.
As expected from the fitted line, M3 can be close to M2 at epoch
1991.433; however, it is not detected in the archive data.
As M4 is only detected at two epochs, we do not perform the fit to
this component. It can be seen from Table 3 that the
inner components move more slowly, and the outer ones move faster.
The variability of this source has been investigated on other
wavebands. Xie et al. (2004) show that the source did not show any
noticeable variability, and it was quiet in the optical band during
their observations. Albert et al. (2006) have detected very high-energy
-Rays from this source, and an optical outburst in 2006
March was also observed but no evidence of flaring was detected by
the University of Michigan Radio Observatory (UMRAO).
![]() |
Figure 3: EVN and VLBA images of 1133+704 at 5 GHz. The epoch of observation from left to right are: epoch 1997.055, epoch 1998.493, epoch 1999.893. The axes are labeled in milliarcseconds. Contours are drawn at -1, 1, 2, 4, 8, 16,... times the noise level. Numerical parameters of the images are given in Table 1. |
Open with DEXTER |
1424+240 - From the images of three epoch VLBI observations
(Fig. 2), the component can not be well-resolved;
however, the jet clearly extends to the southeast with position
angle of about
.
This result is similar with
Fey & Charlot (2000), who modeled this source with two components at 2.3 GHz
and 8.5 GHz, all with jets extending in the nearly same direction.
However, this direction is a little misaligned with the VLA map
found by Rector et al. (2003), which shows a compact structure consisting of
a core and either a halo or roughly collinear jets extending north
and south
.
The MERLIN
archive data (see Fig. 5) shows a very compact structure, and the jet direction is not evident.
1542+614 - There have no previous observations and archive
data for this source. The redshift is unknown. We present our
observations in Fig. 2. The source shows a compact core
and straight jet extending to
on a parsec scale.
4C+37.46 - We present our image in Fig. 2, in
which the source clearly shows one jet towards the northeast
.
The MERLIN archive image (see Fig. 6) also
shows a jet in the same direction. The jet direction is consistent
with Marecki et al. (2006), in which this source could be core-jet-lobe
source, either with a steep spectrum or some diffuse structure. This
source may also be a well aligned source.
1727+502 - The large amplitude variability of this source has
been detected on several occasions. A violent variation of 2.1 mag
in the optical band was reported by Scott et al. (1976). Fan et al. (1999)
detected the variations of
and
mag. This source was observed by Kollgaard et al. (1996a) in
1991 at 6 cm and showed a core-jet morphology. Giroletti et al. (2004a) and
Giroletti et al. (2006) show that this source displays a remarkable alignment
from the parsec to the kilo-parsec scale. The multi-epoch images in
Fig. 7 show a core and a jet extending to northeast with
a position angle of about
.
The modelfit parameters of
this source were presented in Table 2. It can be seen
from Fig. 7 that six components were detected, which are
labeled with C0, C1, C2, C3, C4, and C5. The weighted linear fit of
components C1, C2, C3, C4, and C5 are shown in Fig. 8,
and their measured proper motion and apparent speed are listed in
Table 3. Interestingly, the component C5 was not detected
at epoch 1991.433 (data from Kollgaard et al. 1996a) and 2002.403.
From the fitted line for C5, the position of C5 at epoch 1991.433 is
expected to be close to that of component C4; therefore, it is
possible that these components are mixed together, causing
difficulty in detection. At 2002.403, C5 very likely moved far away
from the core and became very weak, so it was not detected. Due to
the poor uv-coverage of our data and the fact that the beam size at
this epoch is larger than other observations, we did not detect
component C1 at 1999.893 (see Fig. 8). It is apparent
that C3 and C4 at epoch 1998.493 was closer to the core than the
expected distance from the fit. Actually, it seems that these two
components move slowly (or were even stationary) before 1998.493,
but move quickly afterwards, with the proper motion of C4 comparable
to that of C5. However, the present data do not allow us to
investigate this possibility further. From Table 3, we
find that the apparent speed of components increases with the
distance from the core. When looking in detail, components C4 and C5
show superluminal motion, while C1, C2 and C3 show subluminal
motions. C1 is the closest component to the core, and it also shows
the smallest proper motion. Actually, the error in the proper motion
of C1 is larger than its proper motion. Therefore, it is possible
that this component is stationary within the available observations.
![]() |
Figure 4: Positions of components with respect to the core at different epochs from model fitting for 1133+704. The data point at 1991.433 is taken from Kollgaard et al. (1996a). All the other epochs are those presented in this paper, and the lines represent the linear fitting of the motion for each jet component. |
Open with DEXTER |
![]() |
Figure 5:
The MERLIN archive image of 1424+240 at 1.6 GHz observed on 1993 Apr. 1. The contour levels are -1, 1, 2, 4, 8, 16, 32, 64... times 0.001745 Jy beam-1(3![]() |
Open with DEXTER |
![]() |
Figure 6:
The MERLIN archive image of 4C+37.46 at 5 GHz observed on 1999 May 6. The contour levels are -1, 1, 2, 4, 8, 16, 32, 64... times 0.001745 Jy beam-1 (3![]() |
Open with DEXTER |
1741+196 - The six-epoch images of this source in Fig. 9 all show a very straight jet extending to the east
with
.
This is similar to the result in
Rector et al. (2003), in which one epoch VLBA map was presented, and the jet
direction is nearly aligned with the VLA snapshot in Perlman et al. (1996)
with
.
Five components are detected in Fig. 9, which are
labeled as E0, E1, E2, E3, and E4. The weighted linear fit of E1,
E2, E3, and E4 are shown in Fig. 10, and their measured
proper motion and apparent speed are listed in Table 3.
When we perform fits to the components, we exclude the data at epoch
2002.400, because the beam size of this epoch is much larger than
others.
While E4 shows a superluminal motion, E1, E2, E3 all show a rather
small apparent speed. For components E2 and E3, the error in proper
motion is actually larger than their proper motion. In view of this
fact and the rather low value of proper motion, E2 and E3 could be
stationary (within the error) in the time duration covered by the
presenting available observations.
Nieppola et al. (2006) collected a large amount of multi-frequency data
for the objects in the Mets
hovi radio observatory BL Lac
sample. The Mets
hovi BL Lac sample includes 381 objects
selected from the Veron-Cetty & Veron BL Lac Catalogue (Veron-Cetty
& Veron 2000, hereafter VCV2000), and 17 objects from the
literature, of which many sources are from the well-known BL Lac
samples like 1Jy, S4, S5, Einstein Medium Sensitivity Survey (EMSS),
Einstein Slew Survey, and DXRBS. The authors argue that this sample
is supposed to have no selection criteria (other than declination)
in addition to the ones in the original surveys. Based on the
multi-frequency data, the SED of each source were constructed in the
representation. The synchrotron
component of the SED was fitted with a parabolic function
y=Ax2+Bx+C | (1) |
We collected the available data at 330 MHz, 360 MHz, 408 MHz, and
1.4 GHz from the Astrophysical Catalogues Support System (CATs) maintained by the Special
Astrophysical Observatory, Russia, and also the available VLA or
MERLIN core and extended flux for their sample, resulting in a
sample of 170 BL Lac objects. The redshift is known for most of the
sources from NED and the literature. For the sources without
redshift, we adopt the average redshift of 0.473, 0.302, 0.271 for
LBLs, IBLs, and HBLs, respectively. The sample of BL Lac objects is
listed in Table 4: Col. 1 the source IAU name (J2000); Col. 2 the
source alias name; Col. 3 the redshift; Col. 4 the intrinsic
synchrotron peak frequency (see Sect. 4.4); Col. 5 the total 408 MHz radio power; Col. 6 the 5 GHz core luminosity; Col. 7 the 5 GHz or 1.4 GHz extended flux density; Col. 8 the references; Col. 9 the Doppler factor (see Sect. 4.4); Cols. 10-12 the viewing angle
estimated assuming
,
5, and 10, respectively (see Sect. 4.5). In the following sections, we investigate the various properties, e.g. the Doppler factor, viewing angle for this sample.
![]() |
Figure 7: EVN and VLBA images of 1727+502 at 5 GHz. The top from left to right are: epoch 1995.534, epoch 1997.055, epoch 1998.493. The bottom from left to right are: epoch 1999.893, epoch 2002.403. The axes are labeled in milliarcseconds. Contours are drawn at -1, 1, 2, 4, 8, 16,... times the noise level. Numerical parameters of the images are given in Table 1. |
Open with DEXTER |
![]() |
Figure 8: Positions of components with respect to the core at different epochs from model fitting for 1727+502. The data point at 1991.433 is taken from Kollgaard et al. (1996a). All the other epochs are those presented in this paper, and the lines represent the linear fitting of the motion for each jet component. |
Open with DEXTER |
In this work, we use the method suggested by Giovannini et al. (2001) and
Giroletti et al. (2004a). Generally, the low-frequency radio power is a
reliable indicator of the intrinsic radio power as it is only
marginally affected by Doppler-boosted compact components
(Bondi et al. 2001), while the relativistic boosting at the base of the
jet affects the observed radio power of the core. Giovannini et al. (1988)
and Giovannini et al. (2001) find a general correlation between the intrinsic
core radio power and total radio power for a sample of radio
galaxies
![]() |
(2) |
![]() |
Figure 9: EVN and VLBA images of 1741+196 at 5 GHz. The top from left to right are: epoch 1997.055, epoch 1997.378, epoch 1998.493. The bottom from left to right are: epoch 1999.893, epoch 2002.400, epoch 2003.447. The axes are labeled in milliarcseconds. Contours are drawn at -1, 1, 2, 4, 8, 16,... times the noise level. Numerical parameters of the images are given in Table 1. |
Open with DEXTER |
![]() |
Figure 10: Positions of components with respect to the core at different epochs from model fitting for 1741+196. All the data points are presented in this paper, the lines represent the linear fitting of the motion for each jet component except the data point at 2002.400, as the baseline of this epoch is much shorter than the others. |
Open with DEXTER |
The anti-correlation between the synchrotron peak frequency and
the intrinsic power was first suggested by Fossati et al. (1998) and
Ghisellini et al. (1998) for samples of blazars.
Nieppola et al. (2006) find an anti-correlation between the luminosity
at 5 GHz and 37 GHz with
for a large sample of BL
Lac objects. However, the 5 GHz and 37 GHz luminosity can be
severely affected by the beaming effect from the radio jets and
can't be used to indicate of the intrinsic radio power.
Bondi et al. (2001) and Giroletti et al. (2004a) suggest that the low-frequency
radio luminosity of BL Lac objects can indicate the intrinsic power
from the radio jet. However, most of the objects in their works are
HBLs and IBLs. For LBLs, especially those in the 1 Jy sample, as
there is a selection criterion of the radio spectral index being
less than 0.5 (
;
),
the low-frequency radio emission may still be Doppler boosted
(Cao et al. 2003). Indeed, Liu et al. (2006) find that the low-frequency
radio emission of sources with a flat spectrum are likely to be
Doppler-boosted, although only radio-loud quasars are considered in
their work.
To test whether the low-frequency radio emission is Doppler-boosted,
we selected 61 BL Lac objects from Nieppola et al. (2006), which have
VLA observations and extended flux. These sources are listed in
Table 4, including all powerful sources of the 1 Jy
sample, and some other type of BL Lac objects. After extrapolating
from the 1.4 GHz extended flux using
for 15 sources
without 5 GHz extended flux, we plot the 5 GHz extended flux with
the total 408 MHz flux in Fig. 11, which is directly
from the available measurements, or extrapolated from 330 MHz, 360
MHz, or the total 1.4 GHz radio flux (see below). We find a
generally good correlation. However, the scatter is
obviously seen.
This can be caused by the difference in the spectral index of two
bands between sources, while both emissions are not influenced by
the beaming effect. However, the possibility that the low-frequency
radio emission in some sources is indeed Doppler boosted can't be
excluded completely with present data. Nevertheless, we believe that
the low-frequency radio emission is intrinsic, as most of the
sources follow the correlation.
![]() |
Figure 11: The relationship between the radio 5 GHz extended flux and the radio low-frequency 408 MHz flux. |
Open with DEXTER |
To calculate the Doppler factor using the method in Sect. 4.2, we
had to obtain the total 408 MHz luminosity. Unfortunately, not all
sources in our sample had an observed 408 MHz flux. We adopted three
methods of calculating 408 MHz luminosity according to the available
data. For BL lac objects with the extended flux available, we
estimated the 408 MHz luminosity by extrapolating from the 5 GHz or
1.4 GHz extended flux assuming a spectral index of
.
For
the sources with observed low-frequency radio flux, we estimated the
total 408 MHz luminosity after a K-correction using spectral index
.
For the remaining sources, neither the extended flux
nor the low-frequency flux was available, so we extrapolated from
the total 1.4 GHz flux density with a spectral index of
,
which is the average spectral index between 408 MHz
and 1.4 GHz of the BL Lac objects in our sample.
The estimated total 408 MHz radio power is listed in Col. 5 of Table 4. To compare with the FR I radio galaxies sample of
Zirbel & Baum (1995), we converted the total 408 MHz radio luminosity into
the same cosmological frame (
and q0=0), as used in their calculation. The distribution of
this total 408 MHz radio power is shown in Fig. 12.
Although our sample contains about two thirds of the BL Lac objects
from the 1Jy sample, we find that the mean and median values of the
408 MHz radio power (mean
and median
)
are consistent with those of the low-luminosity
FR I radio galaxies in Zirbel & Baum (1995), which is mean
and median
.
This, on
one hand, supports the scenario that FR I radio galaxies are the
parent population of BL Lac objects; and on the other hand, it
indicates that the estimated total 408 MHz radio power of our BL Lac
objects is intrinsic, as it is comparable to FR I radio galaxies.
![]() |
Figure 12: The distribution of 408 MHz luminosity for BL Lac objects selected from Nieppola et al. (2006). |
Open with DEXTER |
The synchrotron peak frequency
,
therefore it may not be intrinsic since
the Doppler factor can vary between sources in our sample. Once the
Doppler factor can be estimated, it would be fundamental in
calculating the intrinsic peak frequency. In this work, we calculate
the intrinsic synchrotron peak frequency
,
where
is from
Nieppola et al. (2006). The factor (1+z) is included since the
of Nieppola et al. (2006) is determined from the log
panel in the observer's frame.
The distribution of the Doppler factor and the intrinsic peak frequency is shown in Fig. 13.
While all IBLs and HBLs have a Doppler factor below 10, LBLs have a
wide spread for the Doppler factor, with most sources below 20, and
it can be over 30 in a few cases. It seems that the Doppler factor
systematically decreases with increasing peak frequency; however, it
is similar for IBLs and HBLs. Indeed, the Spearman rank correlation
analysis shows an anti-correlation between the Doppler factor and
at
99.99 per cent significance level. This
may indicate that the Doppler boosting is systematically more
pronounced in LBLs than in IBLs and HBLs, although a fraction of
LBLs are comparable to IBLs and HBLs.
![]() |
Figure 13:
The Doppler factor versus the intrinsic synchrotron peak frequency
![]() |
Open with DEXTER |
The estimation of Doppler factor enables us to investigate the
relationship between the intrinsic radio power and the intrinsic
synchrotron peak frequency. The relationship between the intrinsic
peak frequency and 408 MHz radio power is shown in Fig. 14. It can be seen that LBLs are statistically more
powerful than other BL Lac objects. There is a significant
anti-correlation between luminosity and
at
99.99% confidence level from the Spearman Rank Correlation
analysis.
Although we used the low-frequency radio emission to indicate the
intrinsic radio power and correct the peak frequency for Doppler
factor, our results in Fig. 14 are very similar to
those of radio power at 5 GHz in Nieppola et al. (2006). Several
low-luminosity LBLs are apparently seen, and they even reach lower
luminosities than any of the HBLs.
We note that the 408 MHz radio power in our sample is calculated in inhomogenous ways, and the redshift of some sources is actually unknown and tentatively taken as the average value of subclasses. To verify how the correlation depends on the data, we check with a sample that has the measured 408 MHz flux and known redshift. Still, these twenty-one sources (Fig. 14) show a strong anti-correlation at about a 98 per cent confidence level. However, due to the small number of these sources, we add the sources with known redshift and 408 MHz radio power extrapolated from either 330 MHz or 360 MHz radio flux (Fig. 14), resulting in a sample of 56 sources, of which the uncertainty from extrapolation is minimized. It is clearly seen that these sources follow the general trend of whole sample. The Spearman rank correlation analysis shows a significant anti-correlation at about a 99.9 per cent confidence level. We are thus confident that the inhomogenous data derivation does not influence our results.
Despite much more scatter in our figure than in the corresponding
one of Fossati et al. (1998), the anti-correlation between luminosity and
is still significant and possibly caused by the
intrinsic physical differences (such as magnetic fields, electron
energies), as suggested by Sambruna et al. (1996). As jet formation is
closely connected with the central engine (Liu et al. 2006), this
correlation can also be caused by the common relationship between
both the radio luminosity
and peak frequency and the central engine, as suggested by Wang et al. (2002)
that the peak frequency is significantly correlated with the
accretion rate.
Although our sample is relatively large, still only a small fraction of low-luminosity LBLs are included in the sample. As a matter of fact, there is no evidence of very high-luminosity HBLs in the sample of Nieppola et al. (2006), which is our parent sample.
![]() |
Figure 14:
The low-frequency 408 MHz luminosity versus the intrinsic
synchrotron-peak frequency
![]() |
Open with DEXTER |
Usually, the viewing angles can be estimated by using the Doppler
factor and the measurements of proper motion, however, the latter
are not available for most of sources. Using the measured proper
motion and the Doppler factor, we calculated the Lorentz factor for
the three sources we observed. We find that the estimated Lorentz
factor is generally consistent with the assumption of ,
and the range of
between 3 and 10 (see Table 3).
This indicates that the angle got by assuming a Lorentz factor of 5 can be reasonable for our sample. Certainly, the adoption of a
Lorentz factor of 5 only makes sense for our statistical sample
investigation and can't be used for an individual source. We believe
that the exact value of Lorentz factor will not change our results
significantly. From Fig. 15, it seems that the BL Lac
objects as a single population possess a continuous range of viewing
angle and a positive correlation between the viewing angle and the
.
Similar to the Doppler factor distribution,
the viewing angle increases with the increasing peak frequency from
LBLs to IBLs, but, it smooths out upward.
In Sect. 3.1, our observations have shown that all the jets of seven BL Lac objects are very straight, and they do not exhibit large bending between the small and large scales. This result is similar with the observations of HBLs in Rector et al. (2003), in which two of our sources are included. This indicates that these sources are either intrinsically straighter or are seen further off-axis than LBL jets. In three sources with measured proper motion (Sect. 3.1), it shows that different components of the same source have different apparent speeds, and the apparent speed increases with the distance from the core. This is similar to the results in Piner et al. (2006), in which BL Lac object 0235+164 and two FSRQs were studied. This phenomena will need a change in either bulk Lorentz factor or viewing angle. However, these changes would also yield predictable changes in the brightness of components due to changing Doppler boosting, and these changes are not seen.
While recent discoveries of low-luminosity LBLs and possible
high-luminosity HBLs (Padovani et al. 2003; Caccianiga & Marchã
2004; Giommi et al. 2005) seems to undermine the blazar sequence,
an anti-correlation between radio power and peak frequency
was found for a large sample of BL Lac objects (Nieppola et al. 2006),
though with a large scatter.
We improved their work by using the intrinsic low-frequency radio
power and the intrinsic synchrotron peak frequency, instead.
It appears that LBLs are intrinsically more powerful than other
types of BL Lac objects, and the distribution of the low-frequency
power of BL Lac objects as a whole is comparable to those of FR I
radio galaxies. Although there is a correlation between the
intrinsic low-frequency power and intrinsic peak frequency, the
scatter of this relation seems large, which contrast with the much
tighter relation of "blazar sequence''. More important, some LBLs
have low radio power, which is difficult to explain by the radiation
cooling process, which was used to explain the blazar sequence. In
fact, these sources tend to weaken the blazar sequence; however,
they are not enough to eliminate the correlation in this work.
While the blazar sequence was suggested for a sample of blazars
including FSRQs, of which the peak frequency can be lower than those
of LBLs, our sample only includes BL Lac objects, although our
sample is much larger. Moreover, the original peak frequency used in
this work is directly adopted from Nieppola et al. (2006), which was
derived by simply fitting the synchrotron component of the SED with
a parabolic function. The determination of peak frequency is thus
influenced by the data quality and sampling. As Nieppola et al. (2006)
discuss, the peak frequencies of the most extreme objects can be
exaggerated by using a simple parabolic function in the fitting
procedure. Furthermore, the large scatter in Fig. 14
implies that other parameters may be at work, and the scatter can be
much lower once this parameter is included. Besides the significant
anti-correlation between the intrinsic low-frequency luminosity and
intrinsic peak frequency, there is a positive correlation between
the viewing angle and intrinsic peak frequency. LBLs apparently have
smaller viewing angle, while the viewing angle of IBLs is comparable
to that of HBLs. We conclude that the intrinsic synchrotron peak
frequency is related not only to the intrinsic luminosity, but maybe
also to the viewing angle. In fact, some HBLs like Mrk 421 and Mrk
501 are supposed to have
,
which was estimated by
Celotti et al. (1998) based on their TeV flux variability. However, the
studies of jet morphology, core dominance, the fit to the trend of
the jet brightness, and Full-Width Half-Maximum (FWHM) all show that
the radio jet emission is oriented to larger angles with respect to
the line of sight (Giroletti et al. 2004b; Giroletti et al. 2006).
![]() |
Figure 15:
The viewing angle versus the intrinsic synchrotron peak frequency
![]() |
Open with DEXTER |
There are only a small fraction of low-luminosity LBLs in our
sample, and actually no evidence of very high-luminosity HBLs
(Nieppola et al. 2006). This precludes us from draws any firm
conclusion about the validity of the blazar sequence. Actually, the
FSRQs with a high
found in DXRBS do not reach the
extreme
values of HBLs. It is still unclear whether
this fact indicates an intrinsic, physical limit to this parameter.
An alternative scenario is one where really
high-power-high-
blazars have their thermal emission
swamped by the non-thermal, featureless jet emission, which makes
their redshift
determination impossible.
While our results seem to undermine the blazar sequence, the
anti-correlation between the intrinsic emission and intrinsic peak
frequency needs further investigation. It is important to compile a
larger sample including various subclasses of blazars, especially
blue quasars, low-luminosity LBLs, and high-luminosity HBLs and to
calculate the accurate peak frequency and intrinsic luminosity. Only
then, the physical mechanism of the different peak frequencies can
be explored. While the peak frequency is believed to be determined
by the magnetic strength and the energy distribution of relativistic
electrons,
the shift of the peak frequency in blazars can give us clues to this
information.
It might also be important to investigate the relationship between
the SEDs of blazars and the central engine, such as accretion rate
(Wang et al. 2002) and/or black hole mass.
Acknowledgements
We thank the anonymous referee for insightful comments and constructive suggestions. This work is supported by the NSFC under grants 10373019, 10333020, and 10543002. The European VLBI Network is a joint facility of European, Chinese, and other radio astronomy institutes funded by their national research councils. The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under cooperative agreement with the National Science Foundation. MERLIN is a National Facility operated by the University of Manchester at Jodrell Bank Observatory on behalf of PPARC. This research 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. This work also made use of Astrophysical Catalogues Support System (CATS) maintained by the Special Astrophysical Observatory, Russia.
Table 1: The VLBI observational log.
Table 2: Model parameters.
Table 3: Proper motions.
Table 4: The sample of 170 BL Lac objects selected from Nieppola et al. (2006).