Issue |
A&A
Volume 512, March-April 2010
|
|
---|---|---|
Article Number | A74 | |
Number of page(s) | 10 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361/200913343 | |
Published online | 07 April 2010 |
Multifrequency observations of a sample of very low frequency peaked BL Lacertae objects
A. Maselli1 - E. Massaro2 - R. Nesci2 - S. Sclavi3 - C. Rossi2 - P. Giommi4
1 - Istituto di Astrofisica Spaziale e Fisica Cosmica, INAF, via U. La
Malfa 153, 90146 Palermo, Italy
2 - Dipartimento di Fisica, Università La Sapienza, Piazzale A. Moro 2,
00185 Roma, Italy
3 - Stazione Astronomica di Vallinfreda, via del Tramonto, Vallinfreda
(RM), Italy
4 - ASI Science Data Center, ASI, ESRIN, via G. Galilei, 00040
Frascati, Italy
Received 24 September 2009 / Accepted 16 December 2009
Abstract
Context. BL Lacertae objects with an
extremely low value of the peak frequency (
Hz)
of the synchrotron component (VLBL objects) are poorly known.
Aims. We selected a small sample of these candidate
objects on the basis of their radio/optical flux ratio to evaluate
their variability, define the low energy part of their SED and obtain a
reliable estimate of their
through a fit with a log-parabolic law.
Methods. They were the object of a ground-based
observational campaign in the infrared and optical bands; data from
several pointings by the Swift satellite
in the UV and X-ray band were also obtained for half the sample. We
complemented our data with a large number of literature data.
Results. These sources showed a marked flux
variability, as expected for BL Lacertae objects of
the low energy peaked BL Lac objects (LBL) and intermediate BL Lac
objects type; all of them have a
value around 1013 Hz,
but this is low enough only for a few sources to consider them
genuine VLBL objects. We report the overall properties of
individual sources and discuss some aspects of this class of rare
objects, which could correspond to the transition between flat spectrum
radio quasars and LBLs.
Key words: galaxies: active - BL Lacertae objects: general - radiation mechanisms: non-thermal
1 Introduction
Blazars are a well known class of active galactic nuclei (AGNs) which
include BL Lac objects and flat spectrum radio quasars
(FSRQs). Their spectral energy distributions (SED), based on data from
the radio band to the highest energy -rays, are characterised by
two main broad bumps.
Following the scenario originally proposed by Blandford & Rees (1978),
the bump peaking at lower frequency is generally explained in
terms of synchrotron emission from relativistic electrons which move
down a jet pointing close to the observers' line of sight. The second
feature, peaking at higher frequencies, is likely due to the inverse
Compton upscattering of low energy photons by the same electron
population. The so-called low energy peaked
BL Lac objects (LBL) (Padovani & Giommi 1995) have the
maximum of the first component in the IR-optical range, where they
generally show a broad and fast variability, while those peaking in the
UV-X rays are named high energy peaked BL Lac objects
(HBL). LBL sources observed in the X-ray band show flat
spectra with an energy index (
)
smaller than unity and mostly exhibit a rather moderate variability in
this band. The SEDs of HBL sources peak as a rule above
a few keV, and their X-ray energy indices are close
to or smaller than unity. Their variability in this band is high in
contrast with the optical and radio data, where their flux is more
stable (Fiorucci et al. 2004).
In addition to these two classes of BL Lacs a finer taxonomy,
always based on the frequency of the synchrotron peak
,
is used in the literature. Following Padovani & Giommi (1995), sources
with log
< 14.5
(measured in Hz) are LBL objects, those with
log
> 16.5 are
HBLs, and those with 14.5 < log
< 16.5 are
classified as intermediate BL Lac objects (IBL). Two more
classes have also been considered in the literature:
the extreme HBL with log
> 18
(Costamante et al. 2001)
and the very low energy peaked BL Lacs (VLBL) with
log
< 13
(Antón & Browne 2005).
We stress that the actual distribution of
can be continuous and that the frequency borders of these classes,
historically influenced by an observational approach, are rather
arbitrary. Our aim therefore is to investigate the properties of some
sources at the low-frequency end of the synchrotron peak range.
The SEDs of FSRQs are mainly similar to those of LBL objects, and only very few peculiar sources like ROXA J081009.9+384757.0 (Giommi et al. 2007) are thought to have the synchrotron peak in the UV-X range. The main difference between LBL objects and FSRQs is therefore the occurrence of broad prominent emission lines in the optical spectra. Some LBL sources in faint states occasionally showed emission lines, which modified the definition of BL Lac, allowing for lines with an equivalent width (EW) smaller than 5 Å (Stickel et al. 1991). However, this limit was exceeded in some occasions by well-known sources like BL Lac itself (Vermeulen et al. 1995) and a few others that have been classified as ``transition objects'' between the two classes of blazars.
For a better understanding of the relationship between BL Lacs and FSRQs we decided to study the spectral properties and the variability of some VLBL candidates. The present observational knowledge of VLBL objects is poor: by definition the radio/optical flux ratio must be high for a VLBL source so that in a radio flux limited sample VLBLs are optically fainter than the other flavours of BL Lacs and therefore less present in optical monitoring samples. Moreover, this picture can be further complicated by considerable brightness changes often associated with an increase of the synchrotron peak frequency in the bright states. Sources potentially classified as VLBLs can therefore occasionally appear as LBLs when brighter and vice-versa.
About the nature of VLBL sources, one possibility is that they
are normal LBLs with a low
because of a loose alignment of the jet to the line of sight, and
therefore with a relatively low Doppler boosting. In this
case, we expect a rather weak optical variability associated
with the absence of high superluminal motions in the jet. The detection
of components in the jet with a high proper motion will be useful for a
reliable estimate of the beaming factor and can be
used to test this interpretation.
There is no well-established and efficient way to select
VLBL candidates. The best approach is of course to estimate
the synchrotron peak from well-sampled spectra, possibly derived from
nearly simultaneous observations. This selection criterion is however
limited to a very small number of
bright sources. The obvious alternative possibility is to use archive
data. Nieppola et al. (2006) fitted
the synchrotron peak of a very substantial sample of radio sources,
reported as BL Lacs in the literature, by means of a parabola
in a double-log plot (i.e. a log-normal distribution,
hereafter named log-parabola) using new radio measurements combined
with archive data to estimate .
According to their results only nine of the 308 sources for
which best fits were computed have log
in the interval from 12.6 to 13.0. There are some
problems of misclassification though, and only a few of them can be
considered as genuine BL Lacs. In fact, three of
these nine sources have spectra with broad lines and are FSRQs, two
have a very low flux density at 1.4 GHz, with radio/optical
and optical/X flux ratios similar to those of HBL sources,
others are classified as blazars of an uncertain type in the Roma-BZCAT
(Massaro et al. 2009).
There is furthermore the possibility that peak frequencies were not
always well-estimated because of the scarcity of data and of badly
resulting best fits, so that some sources with
slightly higher than 1013 Hz are actually
VLBL. In any case, apart from the classification of individual
objects one can reasonably expect that VLBL sources are indeed rare
objects on the basis of the existing data.
To investigate the properties of some VLBL candidates, we
selected a small sample of BL Lacs and performed new optical
and near IR photometry to improve the knowledge of their SED
at frequencies higher than the expected .
We present the results of our observations and estimate the peak
frequency of the synchrotron component.
Table 1: The list of VLBL candidates and the corresponding values of some interesting parameters taken from the literature.
2 Selection of VLBL candidates
The major problem in the search for VLBL objects is the lack of data in
the range of 1012-1014 Hz,
which makes it difficult to estimate the
.
A first and quite simple criterion to select
VLBL candidates is the assumption that the overall SED shape
between radio and optical frequencies is described by a log-parabola
(see Sect. 5). In fact, considering that for
a VLBL object
Hz
and assuming for the ``curvature'' parameter b
(see Eq. (2)) a typical value in the range of
0.15-0.25, it is easy to derive that the expected highest
value of the ratio between the optical flux in the R bandpass
and the radio flux at 5 GHz should be of the order of
5
10-4.
Thus the lowest R magnitude should be
17 for a
radio flux density of 1 Jy.
Another useful quantity to classify BL Lac objects is the ratio between the X-ray and the radio flux. This ratio was used in the Sedentary Survey (Giommi et al. 2005) to select HBL objects and was found to work well also for wider BL Lac samples (Maselli 2009). In the Roma-BZCAT (Massaro et al. 2009) the X-ray flux in the 0.1-2.4 keV band is given in units of 10-12 erg cm-2 s-1, while the radio flux density at 1.4 GHz is in mJy.
To express this ratio in a convenient adimensional number, we
multiplied the radio flux density by a bandwidth of
,
assumed equal to 1 GHz, and the value of the ratio was
opportunely
multiplied by a scale factor equal to 10-3
In this way the value of FX/S1.4 = 10-11 erg cm-2 s-1 Jy-1 corresponds to

We selected therefore from the Roma-BZCAT a small sample of VLBL candidates in the northern sky according to the following criteria:
- BL Lac classification from a featureless optical spectrum;
- S (5 GHz) > 0.25 Jy;
- V > 17.0 mag;
- X-ray to radio (1.4 GHz) flux ratio:
.


The list of the eight selected sources is reported in
Table 1
together with some interesting parameters taken from the literature.
Throughout the paper we named the sources using the
J2000 coordinates as in the Roma-BZCAT;
their usual names are also given to simplify their identification.
Radio fluxes of BZB J0929+8612 and
BZB J1223+8040 are from the S5 survey (Kühr
et al. 1981b)
and for the other sources from GB6 (Gregory et al. 1996); R magnitudes
are from USNO-B1 (Monet et al. 2003); the
reddening values are from NED; the values of
(Eq. (1))
are calculated using RASS data (Voges et al. 1999); redshift
values of BZB J0818+4222 and
BZB J0823+2223 are from Sowards-Emmerd
et al. (2005)
and Stickel et al. (1993)
respectively; estimates of the peak frequency and of the Doppler factor
are from Nieppola et al. (2006) and Wu
et al. (2007),
respectively.
The estimates of
given by Nieppola et al. (2006) are all
higher than
1013 Hz, although many of them are
close to this threshold. We stress that our scientific goal is limited
to the study of the main properties of some VLBL candidates
based on a multifrequency analysis and to verify whether the resulting
estimates of
support this classification.
We also intend to check whether the adopted criteria work well for a
VLBL selection to enrich the number of these relatively rare
AGNs. The sample is obviously too small for statistical studies on the
BL Lacs' population which require analyses of much richer and
possibly complete samples.
Table 2: Magnitudes of VLBL candidates from our observational campaign; the error associated to the B and R filters is 0.05 mag.
For many of these VLBL candidates very little is known about the optical variability. They are not expected to be bright in this range: consequently, they are not included among the main targets of monitoring programmes. Nevertheless, one can expect that they may undergo large optical flares because the emission in this band corresponds to the high frequency tail of the synchrotron jet emission. Therefore we performed an accurate photometry of archive data to derive information on their possible brightness changes over time scales longer than a decade.
Table 3: Magnitudes of VLBL candidates observed by Swift in the six different UVOT filters.
3 Observations and data reduction
We aimed to obtain a better characterisation of the SED in the frequency range where the decay of the synchrotron component is expected for VLBLs. For this reason we carried out a ground-based photometrical campaign in the infrared and optical band for our selected BL Lacs. A posteriori we found for half of them observations by Swift (Gehrels et al. 2004) and performed the reduction of optical-UV and X-ray data acquired by the UVOT (Roming et al. 2004) and XRT (Burrows et al. 2005) telescopes aboard this mission.
3.1 Optical-NIR observations and data reduction
Photometric optical observations were performed in several runs from January 2005 to November 2006 at the Loiano (aperture = 152 cm) and Asiago (aperture = 182 cm) observatories in Italy, both equipped with liquid nitrogen-cooled CCD cameras. Standard B (Johnson) and R (Cousins) filters were used. Exposure times, depending upon the brightness of the source and weather conditions, were 5-10 min long for the R filter and 15-20 min long for the B filter to have a typical noise level on the comparison stars of 0.01 mag or less.
Bias and flat field corrections of the optical frames were
performed with IRAF
tasks. Differential photometry with respect to some comparison stars in
the same field of view of the sources was performed with the DAOPHOT
task in IRAF. The same circular aperture, with a radius of 5
,
was used for the photometry of targets and of the comparison stars.
Magnitude uncertainties in the R and B filters were estimated from the dispersion of measurements of the reference stars in the field of each source. In each field we typically selected five stars with B and R magnitudes in the GSC2 catalogue. Some of them were chosen with a brightness nearly comparable to the target and some others, substantially brighter, were used to put the instrumental magnitudes from different runs on a common scale. The typical dispersion of all the reference stars was always found to be smaller than or equal to 0.05 mag, and therefore we decided to assume this value as the magnitude uncertainty for both the R and B filter.
On three nights simultaneous observations were carried out at Campo Imperatore with the AZT24 telescope (aperture = 110 cm), which is equipped with an infrared CCD liquid nitrogen cooled camera. Standard J, H and K filters were used; typical exposure times were about 2-3 min long. Infrared images were obtained from an appropriate stack of several short exposures performed with a locally ad hoc developed software. The J, H and K magnitudes for the reference stars were taken from the 2MASS catalogue. Each source was observed three times during the night, and the rms deviation from the mean was adopted as photometric error. The results of our photometry in the infrared and optical band are reported in Table 2.
3.2 Optical-UV observations and data reduction
The UVOT instrument (Roming et al. 2004) is able to
obtain images in the V,
B,
U, W1, M2, W2 filters;
still, each sequence of observation may contain only some of them. The
data analysis was performed with dedicated tasks included in the
HEASOFT 6.6.3 package. The instrument typically takes a series
of frames in each filter. We summed up each series with the UVOTIMSUM
task to obtain a single frame in each filter, then performed photometry
using the
UVOTSOURCE task. We followed
UVOT Team prescriptions (Immler et al. 2008) and
adopted an aperture radius of 3
for the photometry, independently of the image filter. We generally
took the background region in the form of an annulus around the
extraction region, with
inner and outer radii chosen in a way to avoid the inclusion of
spurious sources, which are particularly numerous in UV filter
frames; typical values are 18
and 23
,
respectively. The obtained magnitudes were then de-reddened, using the
appropriate value of E(B-V)
derived by Schlegel et al. (1998) and
adopting the extinction curve given in Seaton (1979) with RV
= 3.2, and finally converted to specific fluxes: the results are
reported in Table 3.
Table 4:
The values of the radio spectral index
,
the peak frequency
,
the corresponding
in 10-12 erg cm-2 s-1
and the curvature parameter b
characterising VLBL candidates.
3.3 X-ray observations and data reduction
The XRT observations were carried out in the most sensitive photon counting readout mode (see Hill et al. 2004, for a description of readout modes). Data were reduced with the XRTDAS software package (version 2.4.3) developed at the ASDC (ASI Science Data Center) and distributed within the HEASOFT 6.6.3 package by the NASA High Energy Astrophysics Archive Research Center (HEASARC).
We used the XIMAGE task to detect
the source count rate and the centroid, and to choose a nearby
source-free region to extract the background spectrum. In all the
sources considered here the count rate was well below the photon
pile-up threshold ( cts/s).
Spectral data were then extracted in a circular region within
a 20 pixel radius
(1 pixel =
)
from the centroid; the background spectrum was estimated in a circular
region of 50 pixels. We then obtained all the files necessary
for the spectral analysis with the XRTPRODUCTS
task; the exposure map was taken into account to correct for hot and
damaged CCD pixels and for the vignetting.
In the spectral file, instrumental channels were combined to
include at least 20 counts in each new energy bin.
The spectral analysis was carried out in the 0.3-10 keV energy
range using XSPEC 11.3.2. Four sources of our sample have been
observed by Swift and more than one
observation is generally available for each of them, but only a few are
provided with a suitable exposure time to produce good quality spectra.
These were fitted using an absorbed power law model
after fixing the hydrogen column density at the Galactic value derived
by Dickey & Lockman (1990);
the corresponding value of the photon index
is reported in Table 4.
4 Optical variability
To extend the variability analysis to the longest time scales we retrieved all the scanned images available on-line from the POSS-I and POSS-II surveys and the UK-Schmidt plates in the red (E or F emulsions) and blue (J emulsion) bands for each object. Differential aperture photometry on targets and reference stars for these scanned images was carried out following the same criteria described for our CCD images. Conversion from photographic magnitudes to the CCD magnitude scale was made with a calibration curve for each image by means of the reference stars; the CCD magnitude of the target was determined by interpolation with an uncertainty of 0.2 mag. The results of this photometry are reported in Table 5.
From the combined analysis of historic plates and of our two years' campaign it follows that the candidates show a marked optical variability. In particular four sources (BZB J0818+4222, BZB J0832+4913, BZB J0929+8612 and BZB J2009+7229) varied by more than 1.4 mag, and in three other cases (BZB J0753+5352 and BZB J1927+6117) the variation was confined between 0.8 and 1.4 mag; only BZB J0823+2223 had a moderate variation of 0.4 mag.
Table 5: Magnitudes of VLBL candidates from our photometry on POSS-I, POSS-II and UK-Schmidt plates in the optical band.
5 Spectral energy distributions
In the radio band the SED is generally well-described by a
power law with spectral index
up to frequencies in the range of 1010-1011 Hz.
At higher frequencies, up to the IR-optical band,
it shows a curvature well-described by a log-parabolic law
that can be written as
![]() |
(2) |
where



![]() |
Figure 1: Spectral energy distributions of selected VLBL sources. Power law best fits in the radio and X-ray band are indicated by solid lines, while dashed lines mark the log-parabolic models of the SED. |
Open with DEXTER |
The spectral energy distributions of selected VLBL objects are
plotted in Fig. 1.
Filled coloured symbols are our data and black open symbols are archive
data; red circles mark our ground-based campaign, while blue triangles
refer to POSS data; magenta squares indicate data from UVOT
and XRT aboard Swift. Power law best fits
in the radio and X-ray band are indicated by solid lines, while dashed
lines mark the log-parabolic models of the SED. The low frequency
portion of the log-parabolic fit was limited to the point where it
assumed the same slope of
,
and we used a power law below this frequency. For
sources with non-simultaneous data showing a significant variability at
radio and optical frequencies, we computed the log-parabolic fit to
data points selecting low and high states in both ranges.
A different choice, for instance that of combining a low state
in the radio band with a high optical state, would introduce a bias
which
would increase the value of
and vice-versa. Best-fit parameter values for all the sources are
reported in Table 4.
6 Comments on individual sources
For an accurate description of the synchrotron component in the SED we collected a large number of non-simultaneous multifrequency data. We took measurements from well-known catalogues, easily accessible through NED and SIMBAD databases, and integrated them with those extracted from literature articles where available. In this section we briefly describe the main properties of the sources in our sample.
6.1 BZB J0753+5352 (OI 582)
This source has been classified as a BL Lac object by
Kühr & Schmidt (1990)
because of a featureless optical spectrum and high optical
polarisation. The absence of any relevant emission feature in the
optical spectrum was later confirmed by Henstock et al. (1997) and also
by Stickel & Kühr (1993),
who fixed the lower limit for the redshift at
due to the undetection of the host galaxy. Radio data, which are
available up to 37 GHz (Teräsranta et al. 2005;
Nieppola et al. 2007),
evidence a very fast rising of the synchrotron emission component (
= -0.60).
The source was detected by ROSAT and is included both in the Faint
Source Catalog (FSC, Voges et al. 2000) and in the
RGB sample (Laurent-Muehleisen et al. 1999). Reich
et al. (2000)
reported an integrated flux density of
10-13 erg cm-2 s-1
in the 0.1-2.4 keV band, with a power law
photon index
.
Three recent observations, separated by a very short time interval,
have been carried out by Swift -XRT: the
corresponding X-ray spectra show a very similar behaviour, and we
reported in the SED only data relative to the observation with the
higher statistics. Our optical data indicate a moderate variability,
and the two best fits of the faint and high state (Fig. 1) are rather
similar with
higher
than 1013 Hz: therefore it
cannot be properly considered a VLBL object.
6.2 BZB J0818+4222 (OJ 425)
This source, observed in the radio and submillimetric band up to
230 GHz (Steppe et al. 1992) shows a
flat radio spectrum with a flux density at 5 GHz of about
3.3 Jy (Kühr et al. 1981a). The
combined analysis of data from the POSS plates and our campaign
indicates a variability range of more than two magnitudes. Available
optical spectra (Stickel et al. 1993; Lawrence
et al. 1996;
Rector & Stocke 2001;
SDSS-DR6) are featureless and redshift estimates are very uncertain
because the host galaxy is undetected, suggesting the lower limit to be
(Falomo et al. 1997).
Its VLBL classification is controversial: we computed two SED
best fits for two brightness states but excluded the high flaring
values observed by us (Fig. 1), and we found
the
value in the low
state to exceed 1013 Hz, while it was
below this threshold in the high state. Such a behaviour
(redder when brighter) is in contrast with other well studied
LBL objects, but we cannot exclude that this could be due to
data selection. BZB J0818+4222 was detected in the
X-ray band by ROSAT (Urry et al. 1996) with
1 keV flux of 0.05
Jy, but the spectral slope is very uncertain with
a photon index of
1.16+0.67-0.90
in the energy range of 0.1-2.0 keV. The recent
observation by Swift XRT reported here
caught this source in a somewhat higher state, with a flux of
0.07
Jy
and a steeper spectral index in the 0.3-10 keV energy range.
6.3 BZB J0823+2223 (OJ 233)
It is a bright radio source with a rather unusual radio structure
(Gabuzda et al. 2001).
Its uncertain high redshift (z=0.951) was
estimated from a possible
MgII emission line (Stickel et al. 1993). The
results of our photometry confirm a low brightness level. Considering
also archive values, they point out a mild variability of about
0.4 mag. The radio spectrum is the steepest one, and one
cannot exclude that this value is due to the presence of some extended
emission regions, as shown by the above quoted radio images.
Both values are very
similar and lower than 1013 Hz:
for this reason this source can be considered an interesting
VLBL candidate.
6.4 BZB J0832+4913 (OJ 448)
The redshift z=0.548 was estimated from a possible
[OIII] emission line (Stickel et al. 1993).
A large variability in the optical band of about two
magnitudes was reported
by Walsh et al. (1984),
who found it at a magnitude of R
16.9 mag. Our observations confirmed such a bright state with
a variation range of
1.7 mag,
whereas the SDSS photometry corresponds to a much fainter
state. Radio data allow us to establish two rather well-defined
brightness states with similar spectral indices, and therefore we
computed best
fits for both conditions (Fig. 1). Spectral
curvatures were similar, but the peak frequency in the high state is
above our VLBL threshold: consequently, this source does
not properly belong to this class of BL Lacs, particularly if
the high z value is confirmed by
new spectroscopic data.
6.5 BZB J0929+8612 (S5 0916+86)
Literature data in the radio band are few and quite scattered, and our
photometry is limited only to the optical band. Consequently, the
spectrum cannot be well-constrained at low frequencies and the
log-parabolic best fits of the SED can be more uncertain than for the
other sources. The resulting
are always lower than 1013 Hz, and the
high value of
given by Nieppola et al. (2006) is not
confirmed. This source can be considered as a good
VLBL candidate.
6.6 BZB J1223+8040 (S5 1221+80)
A featureless optical spectrum has been reported by Henstock
et al. (1997).
Below 1 GHz
the spectrum is steeper than at higher frequencies: this can be
considered as an indication of more components with a spectrum
different from the radio core, likely the dominant
structure at high frequencies. For these reasons, radio data at
frequencies lower than 4
108 Hz were not included in
the fit. We observed this source at a few epochs and found an
indication for a variability of about one magnitude, but the
photometric accuracy in the faint
state is rather poor. Values of
higher
than 1014 Hz, as derived by Nieppola
et al. (2006),
have not been found by us: the peak frequency of the
unique best fit (Fig. 1)
is slightly lower than 1013 Hz, hence
this source could be a VLBL object, but this needs to be
confirmed by other data, possibly in different brightness states.
6.7 BZB J1927+6117 (S4 1926+61)
This source was classified as a BL Lac object by Kühr
& Schmidt (1990).
The absence of detectable features both in emission and in absorption
was later confirmed by all subsequent spectra. No host galaxy
has ever been evidenced, suggesting z>0.2.
The radio spectrum is very flat across a large frequency range, which
makes it detectable in the microwave range by WMAP (Bennett
et al. 2003).
In the X rays, Britzen et al. (2007) reported
a flux of
10-12 erg cm-2 s-1
in the 0.1-2.4 keV range. We observed it on several occasions
and found a rather stable brightness with a variation range within
0.3 mag. However, it was characterised by a faint B-R colour,
which implies a flat optical spectrum, which in turn does not follow
the log-parabolic law connecting radio points. A similar flat
spectral distribution follows from SDSS photometry. Assuming
that the log-parabola must fit the K flux
(Fig. 1)
we obtained a value of the curvature parameter higher than the other
sources and a peak frequency below the VLBL limit.
6.8 BZB J2009+7229 (S5 2010+72)
This source was generally very faint, and only on a few occasions its R magnitude was brighter than 19 mag. It was undetected in all the bandpasses of Swift-UVOT, and archive data indicate a comparable luminosity level, considering their large uncertainties. Two possible SEDs are shown in Fig. 1 with peak frequencies around the limit value of 1013 Hz. It is therefore a possible VLBL candidate, but more data in different luminosity states are necessary to confirm its classification.
7 Discussion
We have presented new observational results covering the frequency interval from the near infrared to the X-ray band of a small sample of BL Lacs, which were selected to be VLBL candidates. We have also considered a collection of archive data useful to complete the description of their spectral energy distributions and, more specifically, to obtain a good estimate of the peak frequencies of the synchrotron component.
The analysis of recent Swift-XRT observations of four sources of the sample (one of them, BZB J2009+7229, was not included in the RASS) confirms that the optical/X-ray ratios and the X-ray spectral indices are similar to those expected from LBL sources. This result supports that the X-ray emission can be explained by inverse Compton scattering in the framework of synchrotron-self Compton models.
Table 6:
Absolute magnitudes MR
of sample sources calculated adopting R1
values from USNO-B1 (Table 1) and rest-frame
peak frequencies
adopting values relative to the faint state, as reported
in Table 4.
Optical monitoring and comparison with historic data show that all our candidate VLBL sources have a significant variability comparable with LBL objects: this implies that a jet disalignment with respect to the line of sight should not be the most likely explanation for their high radio/optical flux ratio.
Out of our eight candidates, four showed a
below 1013 Hz, supporting their
identification as VLBL sources. For three other sources the
peak frequency was found in one
observation lower and in the other one slightly higher than this limit;
only one source was found (slightly) above it. The selection
criteria adopted to build our sample have therefore a
satisfactory efficiency. We cannot identify any systematic trend
of
with the brightness level of the sources.
An important point to be taken into account is that the
redshift for the majority of our sources is unknown and that the few
estimated values (see Table 1) are relatively
high for BL Lac objects. The distribution of
known z for the BL Lacs in the Roma-BZCAT
(Massaro et al. 2009)
peaks around 0.3 and decays fast for higher values.
If the typical redshifts of our sources were actually higher
than 0.3, their rest-frame frequencies
would result close to or higher than the VLBL limit.
We calculated the absolute magnitudes of the sources included
in our VLBL sample adopting the R magnitudes
from USNO-B1 and assuming z=0.4 for all sources
with unknown redshift. To compare our results with the study
of the well-known sample of host galaxies of
BL Lac objects with the Hubble Space Telescope (Urry
et al. 2000),
we also assumed H0 =
50 km s-1 Mpc-1
and q0=0. With this
comparison we expect that the absolute magnitudes of our sources are
substantially brighter than the typical host galaxy. Two points seem
remarkable from the results which are reported in Table 6:
i) there is no clear correlation between
and MR;
ii) four sources have an MR
substantially fainter than -23.7, the typical absolute
magnitude of the host galaxies of BL Lacs derived by Urry
et al. (2000).
This second point may imply either that the actual redshifts of our
sources are even higher than 0.4, in which case the
rest-frame
would be higher than 1013 Hz,
or that we are dealing with sources with a below-the-average host
galaxy. This could be related to the low peak frequency of these AGNs,
but a more substantial sample of sources with known
redshifts would be necessary to shed light on this topic.
A comparison between the absolute luminosity of our
VLBL candidates and FSRQs is difficult from a statistical
point of view due to the lack of actual redshift information for most
of our sources and to the observational biases present in the
AGN samples. Assuming the radio/optical flux ratio
as a proxy of
,
we selected from the Roma-BZCAT a sample of
150 FSRQs with z<1 and
103; radio fluxes were taken from NVSS (Condon
et al. 1998)
and R1 optical magnitudes from USNO-B1. The median
value of FSRQs
was found to be brighter than that of our small sample
;
still, one has to take into account that the rms scatter of
the former sample (
mag) is
quite high.
Wu et al. (2007)
and Nieppola et al. (2008)
argued that the Doppler factor
in AGNs is anticorrelated with
,
albeit a large scatter: on average, QSOs are more boosted than
LBLs, which in turn are more boosted than HBLs. Therefore one would
expect that VLBL sources have the highest Doppler factors
among BL Lacs and could be the transition case between BL Lacs
and FSRQs. According to the measurements reported in Hovatta
et al. (2009,
their Table 2)
the median Doppler boosting factor for BL Lac Objects
is
.
We found an estimate of this factor for seven of our sources in Wu
et al. (2007).
These were derived from the correlation between the core and the total
low-frequency radio power for a sample of FRI galaxies
(Giovannini et al. 2001)
which holds also for a sample of BL Lacs as shown by Giroletti
et al. (2004).
We reported these estimates in Table 1: the median
value
agrees reasonably well with expectations, given the possible biases
intrinsic to the method and the small size of the sample. The
evolutionary model of blazars proposed by Cavaliere & D'Elia (2001)
considers BL Lacs as FSRQs in which the broad line region is
depleted. In the hypothesis that VLBLs are FSRQs approaching the last
phase of depletion we would expect a median redshift moderately higher
than BL Lac's and a Doppler factor lower than FSRQ's,
consistent with our findings. Clearly a larger and better defined
statistical sample will be necessary to confirm these indications.
Very useful for a reliable identification of a population of
VLBL sources will be the data obtained with the recent PLANCK
and HERSCHEL space observatories. These instruments will in fact cover
the frequency range where the SED peak frequency of VLBLs is expected,
and therefore their data will for the first time provide a direct
information on it. Finally we stress that one cannot exclude
that a weak -ray
emission from some VLBL candidates could also be detected in
the coming years by the LAT instrument on board Fermi-GST.
A detection will help to establish the nature of
these AGNs.
8 Conclusions
The main results of this work can be summarised as follows:
- VLBL sources are very rare in blazar catalogues: this can be due to observational selection effects, classification issues or a real scarcity of sources;
- we selected a small sample of VLBL candidates from radio/optical fluxes and made an optical-NIR and NUV-X-ray monitoring;
- the variability behaviour of the sources was similar to that of normal LBL objects;
- the peak frequency of the synchrotron emission
of our sample objects is below or very close to 1013 Hz, as expected from our selection criteria, confirming their VLBL nature;
- the redshifts of most of our sources are unknown but they must be higher than the typical values of classified BL Lacs, if their host galaxies are similar to those of the known BL Lacs;
- if the anticorrelation between
and Doppler beaming holds for VLBLs, these may be transition objects between FSRQ and LBL; in this scenario a relatively high redshift would naturally be expected for these sources; a fast transition could explain their rarity;
- given the large cosmic scatter of the above correlation, a number of low-power VLBL sources may still be hidden in poorly studied radio galaxies.
We thank the referee Elina Nieppola for constructive criticism which helped to improve the paper. 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, and of the SIMBAD database which is operated at CDS, Strasbourg, France. The authors acknowledge the financial support from ``Sapienza'' - Università di Roma.
References
- Antón, S., & Browne, I. W. A. 2005, MNRAS, 356, 225 [NASA ADS] [CrossRef] [Google Scholar]
- Bennett, C. L., Hill, R. S., Hinshaw, G., et al. 2003, ApJS, 148, 97 [NASA ADS] [CrossRef] [Google Scholar]
- Blandford, R. D., & Rees, M. J. 1978, Proceedings of Pittsburgh Conference on BL Lac Objects [Google Scholar]
- Britzen, S., Brinkmann, W., Campbell, R.M. et al., 2007, A&A, 476, 759 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Burrows, D. N., Hill, J. E., Nousek, J. A., et al. 2005, Space Sci. Rev., 120, 165 [NASA ADS] [CrossRef] [Google Scholar]
- Cavaliere, A., & D'Elia, V. 2002, ApJ, 571, 226 [NASA ADS] [CrossRef] [Google Scholar]
- Condon, J. J., Cotton, W. D., Greisen, E. W., et al. 1998, AJ, 115, 1693 [NASA ADS] [CrossRef] [Google Scholar]
- Costamante, L., Ghisellini, G., Giommi, P., et al. 2001, A&A, 371, 512 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Dickey, J., & Lockman, F. 1990, ARA&A, 28, 215 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- Falomo, R., Urry, M. C., Pesce, J. E., et al. 1997, ApJ, 476, 113 [NASA ADS] [CrossRef] [Google Scholar]
- Fiorucci, M., Ciprini, S., & Tosti, G. 2004, A&A, 419, 25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gabuzda, D. C., Pushkarev, A. B., & Garnich, N. N. 2001, MNRAS, 327, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Gehrels, N., Chincarini, G., Giommi, P. et al., 2004, ApJ, 611, 1005 [NASA ADS] [CrossRef] [Google Scholar]
- Giommi, P., Piranomonte, S., Perri, M. & Padovani, P., 2005, A&A, 434, 385 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Giommi, P., Massaro, E., Padovani, P. et al., 2007, A&A, 468, 97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Giovannini, G., Cotton, W. D., Feretti, L., Lara, L., & Venturi, T. 2001, ApJ, 552, 508 [NASA ADS] [CrossRef] [Google Scholar]
- Giroletti, M., Giovannini, G., Taylor, G. B., & Falomo, R. 2004, ApJ, 613, 752 [NASA ADS] [CrossRef] [Google Scholar]
- Gregory, P. C., Scott, W. K., Douglas, K., & Condon, J. J. 1996, ApJS, 103, 427 [NASA ADS] [CrossRef] [Google Scholar]
- Henstock, D. R., Browne, I. W. A., Wilkinson, P. N., & McMahon, R. G. 1997, MNRAS, 290, 380 [NASA ADS] [Google Scholar]
- Hill, J. E., Burrows, D. N., Nousek, J. A., et al. 2004, Proc. SPIE, 5165, 217 [NASA ADS] [CrossRef] [Google Scholar]
- Hovatta, T., Valtaoja, E., Tornikoski, M., & Lähteenmäki, A. 2009, A&A, 498, 723 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Immler, S., et al. 2006, The Swift -UVOT Software guide, version 2.2 [Google Scholar]
- Kühr, H., & Schmidt, G. D. 1990, AJ, 99, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Kühr, H., Witzel, A., Pauliny-Toth, I. I. K., & Nauber, U. 1981a, A&AS, 45, 367 [NASA ADS] [Google Scholar]
- Kühr, H., Pauliny-Toth, I. I. K., Witzel, A., & Schmidt, J. 1981b, AJ, 86, 6 [Google Scholar]
- Lamer, G., Brunner, H., & Staubert, R. 1996, A&A, 311, 384 [NASA ADS] [Google Scholar]
- Landau, R., Golisch, B., Jones, T. J. et al., 1986, ApJ, 308, 78 [NASA ADS] [CrossRef] [Google Scholar]
- Laurent-Muehleisen, S. A., Kollgaard, R. I., Feigelson, E. D., et al. 1999, ApJ, 525, 127 [NASA ADS] [CrossRef] [Google Scholar]
- Lawrence, C. R., Zucker, J. R., Readhead, A. C. S., et al. 1996, ApJS, 107, 541 [NASA ADS] [CrossRef] [Google Scholar]
- Marchã, M. J., Caccianiga, A., Browne, I. W. A., & Jackson, N. 2001, MNRAS, 326, 1455 [NASA ADS] [CrossRef] [Google Scholar]
- Maselli, A. 2009, Ph.D. Thesis, Sapienza, Università di Roma [Google Scholar]
- Massaro, E., Perri, M., Giommi, P., & Nesci, R. 2004a, A&A, 413, 489 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Massaro, E., Perri, M., Giommi, P., et al. 2004b, A&A, 422, 103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Massaro, E., Giommi, P., Leto, C. et al. 2009, A&A, 495, 691 (http://www.asdc.asi.it/bzcat) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Monet, D. G., Levine, S. E., Canzian, B., et al. 2003, AJ, 125, 984 [NASA ADS] [CrossRef] [Google Scholar]
- Nieppola, E., Tornikoski, M., & Valtaoja, E. 2006, A&A, 445, 441 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nieppola, E., Tornikoski, M., Lähteenmäki, A., et al. 2007, AJ, 133, 1947 [NASA ADS] [CrossRef] [Google Scholar]
- Nieppola, E., Valtaoja, E., Tornikoski, M., et al. 2008, A&A, 488, 867 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Padovani, P., & Giommi, P. 1995, ApJ, 444, 567 [NASA ADS] [CrossRef] [Google Scholar]
- Rector, T. A., & Stocke, J. T. 2001, AJ, 122, 565 [Google Scholar]
- Reich, W., Fürst, E., Reich, P., et al. 2000, A&A, 363, 141 [NASA ADS] [Google Scholar]
- Roming, P. W. A., Hunsberger, S. D., Mason, K. O., et al. 2004, in X-Ray and Gamma-Ray Instrumentation for Astronomy XIII, Proc. SPIE, 5165, 262 [Google Scholar]
- Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 [NASA ADS] [CrossRef] [Google Scholar]
- Seaton, M. J. 1979, MNRAS, 187, 73 [Google Scholar]
- Sowards-Emmerd, D., Romani, R. W., Michelson, P. F., et al. 2005, ApJ, 626, 95 [NASA ADS] [CrossRef] [Google Scholar]
- Steppe, H., Liechti, S., Mauersberger, R., et al. 1992, A&AS, 96, 441 [NASA ADS] [Google Scholar]
- Stickel, M., & Kühr, H. 1993, A&AS, 101, 521 [NASA ADS] [Google Scholar]
- Stickel, M., Fried, J. W., & Kühr, H. 1989, A&AS, 80, 103 [Google Scholar]
- Stickel, M., Padovani, P., Urry, C. M., et al. 1991, ApJ, 374, 431 [NASA ADS] [CrossRef] [Google Scholar]
- Stickel, M., Fried, J. W., & Kühr, H. 1993, A&AS, 98, 393 [Google Scholar]
- Teräsranta, H., Wiren, S., Koivisto, P., et al. 2005, A&A, 440, 409 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Urry, M. C., Sambruna, R. M., Worrall, D. M., et al. 1996, ApJ, 463, 424 [NASA ADS] [CrossRef] [Google Scholar]
- Urry, M. C., Scarpa, R., O'Dowd, M., et al. 2000, ApJ, 532, 816 [NASA ADS] [CrossRef] [Google Scholar]
- Vermeulen, R. C., Ogle, P. M., Tran, H. D., et al. 1995, ApJ, 452, L5 [NASA ADS] [CrossRef] [Google Scholar]
- Voges, W., Aschenbach, B., Boller, T., et al. 1999, A&A, 349, 389 [NASA ADS] [Google Scholar]
- Voges, W., Aschenbach, B., Boller, T., et al. 2000, IAUC, 7432, 1 [Google Scholar]
- Walsh, D., Beckers, J. M., Carswell, R. F., & Weymann, R. J. 1984, MNRAS, 211, 105 [NASA ADS] [Google Scholar]
- Wu, Z., Jiang, D. R., Gu, M., & Liu, Y. 2007, A&A, 466, 63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Footnotes
- ... IRAF
- Image Reduction and Analysis Facility, distributed by NOAO, operated by AURA, Inc. under agreement with the US NSF.
All Tables
Table 1: The list of VLBL candidates and the corresponding values of some interesting parameters taken from the literature.
Table 2: Magnitudes of VLBL candidates from our observational campaign; the error associated to the B and R filters is 0.05 mag.
Table 3: Magnitudes of VLBL candidates observed by Swift in the six different UVOT filters.
Table 4:
The values of the radio spectral index
,
the peak frequency
,
the corresponding
in 10-12 erg cm-2 s-1
and the curvature parameter b
characterising VLBL candidates.
Table 5: Magnitudes of VLBL candidates from our photometry on POSS-I, POSS-II and UK-Schmidt plates in the optical band.
Table 6:
Absolute magnitudes MR
of sample sources calculated adopting R1
values from USNO-B1 (Table 1) and rest-frame
peak frequencies
adopting values relative to the faint state, as reported
in Table 4.
All Figures
![]() |
Figure 1: Spectral energy distributions of selected VLBL sources. Power law best fits in the radio and X-ray band are indicated by solid lines, while dashed lines mark the log-parabolic models of the SED. |
Open with DEXTER | |
In the text |
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.