Issue |
A&A
Volume 521, October 2010
|
|
---|---|---|
Article Number | A2 | |
Number of page(s) | 6 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361/200912706 | |
Published online | 13 October 2010 |
The dynamical ages of two low-luminosity young radio sources
N. de Vries1 - I. A. G. Snellen1 - R. T. Schilizzi1,2 - K.-H. Mack3
1 - Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands
2 -
Square Kilometre Array Program Development Office, Jodrell Bank Centre for Astrophysics, Manchester M13 9PL, UK
3 -
Istituto di Radioastronomia - INAF, Via P. Gobetti 101, 40129 Bologna, Italy
Received 16 June 2009 / Accepted 1 June 2010
Abstract
Context. Young radio-loud active galactic nuclei form an
important tool for investigating the evolution of extragalactic radio
sources. To study the early phases of expanding radio sources, we
are using the CORALZ sample, which contains 25 compact (
)
radio sources associated with nearby galaxies.
Aims. In this paper we aim to determine the lobe-to-lobe
expansion speeds for two sources in the sample. Comparing the dynamical
ages of these relatively low-luminosity objects with ages of powerful
sources found in the literature, will serve to test radio source
evolution models.
Methods. We present our third epoch of VLBI observations, using
the EVN+VLBA array at 5.0 GHz, at four and eight years after
our previous measurements.
Results. We report the detection of a lobe-to-lobe expansion speed of 0.095
0.020 h-1 c for CORALZ J083139+460800, and a 1-
upper limit to the expansion speed of CORALZ J131739+411545 of
0.074 h-1 c. The detected expansion speed of CORALZ J083139+460800 corresponds to a dynamical age of 245
55 years.
Conclusions. The detected expansion speed and upper limit are
lower than any presented in the literature for powerful young AGN. This
confirms that the expansion speeds of young radio sources are
correlated with their radio luminosities.
Key words: galaxies: active - radio continuum: galaxies
1 Introduction
Young radio-loud active galactic nuclei (AGN) are key objects for constraining evolution models of classical double radio sources. Young radio sources can be recognised by their small sizes, by their inverted spectra at low frequencies, and by their often compact symmetric morphologies, with lobe-like structures on both sides of a central core. Depending on the observed characteristics, they are called compact steep spectrum (CSS) sources (Peacock & Wall 1982; Fanti et al. 1990; O'Dea 1998), gigahertz peaked spectrum (GPS) sources (O'Dea 1998), and/or compact symmetric objects (CSO, Wilkinson et al. 1994; Conway 2002).
Samples of young radio-loud AGN have been used in a statistical way to constrain the early evolution of radio sources. Relative number statistics indicate that at a young age, radio sources must have significantly higher radio luminosities than when they are older (O'Dea & Baum 1996; Readhead et al. 1996; Fanti & Fanti 2003). Some analytic models predict that this decrease in radio luminosity is preceded by a period of increasing luminosity at the earliest evolutionary stages (Snellen et al. 2000, 2003).
Very important results come from multi-epoch Very Long Baseline Interferometry (VLBI) observations of individual CSOs. For a handful of powerful objects, the proper motions of the lobes have been detected, from which their dynamical ages can be estimated (Owsianik & Conway 1998; Polatidis & Conway 2003, and references therein). Separation velocities between the opposite extremities of the sources of up to 0.5 h-1 c have been measured, corresponding to dynamical ages in the range of a few hundred to a few thousand years.
Snellen et al. (2004) selected a sample of 25 young radio sources, CORALZ (COmpact RAdio sources at Low Redshift), with the aim of obtaining an unbiased view of young radio-loud AGN in the nearby universe. The sources in the CORALZ sample are significantly closer than the archetypal CSO and GPS sources studied in the literature. Proper motions should be easier to detect in these nearby sources, since similar intrinsic expansion velocities would result in significantly higher angular motions than in previous studies. Moreover, in general these radio sources are less powerful than the archetypal CSOs and GPS sources, which will allow luminosity-dependent comparison studies.
The CORALZ sample was selected on flux density (
> 100 mJy) and angular size (
). Using the Very Large Array (VLA) Faint Images of the Radio Sky at Twenty-centimeters (FIRST) survey (White et al. 1997),
the optical Automated Plate Measuring machine (APM) catalogue of the
first Palomar Observatory Sky Survey (POSS-I) (McMahon & Irwin 1991), and follow-up observations, all radio sources identified with bright galaxies were selected (red magnitude of e < 16.5 mag or blue magnitude of o < 19.5 mag).
A detailed description of the selection process and additional radio observations can be found in Snellen et al. (2004).
Nearly all of the sources in the CORALZ sample are spectrally
classified as GPS and CSS sources, and morphologically classified
as CSOs, confirming that these are likely to be young radio
sources.
In our previous paper (de Vries et al. 2009) we presented radio maps, morphological classifications, and linear sizes of the CORALZ sample. Furthermore we derived a remarkably low upper limit to the expansion speed of one of the sources (CORALZ J083139+460800, which is re-examined in this paper). Since the CORALZ sample is of relatively low luminosity, we interpreted this as a hint that the expansion speeds of young radio sources are correlated with their radio luminosities. We now wish to investigate this further.
In this paper we present a new epoch of global VLBI
observations for two sources from the CORALZ sample. For these two
radio sources, for which in total three epochs of VLBI observations
have been obtained with the same telescope array at the same frequency,
we will estimate expansion velocities. The new observations of the
other sources from the CORALZ sample will be published at a later date,
when sufficient data (i.e. a third epoch of observations) has been
obtained for the kinematical study of those sources. In Sect. 2 the observations and data reduction are described; the data analysis is outlined in Sect. 3. The results are presented and discussed in Sect. 4, and we conclude with Sect. 5. Throughout the paper we adopt the cosmological parameters as found by WMAP5 (Komatsu et al. 2009;
,
,
).
2 Observations and data reduction
On March 13, 2008, we observed seven radio sources from the CORALZ
sample with global VLBI at 5.0 GHz with a bandwidth of
16 MHz, using the EVN (European VLBI Network) and the VLBA (Very
Long Baseline Array). A similar observing strategy was used as for
the previous global VLBI observations (de Vries et al. 2009),
in order to minimize differences in signal-to-noise and
resolution. Just like our observations in 2004, this global VLBI
experiment was conducted with an array of 16 stations, using the
EVN antennas Effelsberg, Jodrell Bank, Medicina, Noto, Torun,
Westerbork, and Onsala, and the VLBA antennas St. Croix, Hancock, North
Liberty, Fort Davis, Los Alamos, Pie Town, Kitt Peak, Owens Valley, and
Brewster. The largest projected baselines were typically 9000 km,
resulting in a resolution of about 2 mas. Each source was
typically observed for 4
10 min, spread in time to obtain optimal UV coverage.
4C 39.25 and J1657+5705 were observed as primary calibrators. The
data correlation was performed at the Joint Institute for VLBI in
Europe (JIVE). The Astronomical Image Processing System (AIPS) has been
used for editing, a-priori calibration, fringe-fitting,
self-calibration and imaging of the data. For fringe-fitting, the
AIPS task FRING was used with a point source model, a solution interval
of 6 min, and a standard signal-to-noise ratio cutoff
of 5. Several iterations of phase (or amplitude & phase)
self-calibration and imaging were performed using the AIPS tasks CALIB
and IMAGR, until the image quality converged, with negligible negative
structure and a low noise level of about 0.2 mJy.
3 Data analysis
The observed radio source structures were characterised using the
AIPS task JMFIT, which fits (by least-squares) up to four Gaussian
components to an image subsection. A first guess of the input
parameters (position, flux density, and shape of the components) was
made using the INPFIT procedure. To obtain error estimates of the
component positions we have used the formula from Fomalont (1999):
where


4 Results and discussion
We present the third epoch of VLBI observations for two sources from the CORALZ sample, CORALZ J083139+460800, and CORALZ J131739+411545. These have also been observed as part of the VLBA Imaging and Polarimetry Survey (VIPS, Helmboldt et al. 2007). In Figs. 4 and 5 we show all maps of CORALZ J083139+460800 and CORALZ J131739+411545, respectively, including our new epoch of observations. In the following sections the results from the kinematic analyses of these sources are presented and discussed.
4.1 CORALZ J083139+460800
CORALZ J083139+460800 is perfectly suitable for kinematic
analysis. It exhibits a compact double morphology with two strong,
unresolved components separated by about 4.4 mas. We
rereduced our May 2004 data for this source, since an additional
(possibly spurious) component was cleaned, not detected in the other
epochs, biasing the positions of the two radio lobes. The two hotspot
separation measurements differ with a 3 significance.
Note that the error estimate is based on the reduced map only, and does
not include possible biases introduced during the synthesis of the map,
as was noted by de Vries et al. (2009).
In Table 1 the
measured component separations for CORALZ J083139+460800
are presented, with in Col. 1 the observation date, in Col. 2
the post-fit rms error of the map, in Col. 3 the measured
hotspot separation in milliarcseconds, and in Col. 4 the 1- uncertainty of the hotspot separation. Figure 1
shows the measured component separations in milliarcseconds for
CORALZ J083139+460800 from March 2000 to
March 2008, with the vertical bars indicating the 1-
uncertainties. The data are in excellent agreement with a linear expansion.
The solid line in Fig. 1 is the result of a linear least-squares fit to the data, with a slope of 16
3 microarcsec per year, corresponding to a lobe-to-lobe, projected expansion rate of 0.095
0.020 h-1 c = 0.136
0.029 c at this redshift (z=0.127).
It is the lowest expansion speed measured in any GPS radio source
so far, supporting the dependence of the expansion speeds of young
radio sources on their radio luminosities. The projected expansion
speed of
0.029 c,
together with the measured projected hotspot separation of
4.4 milliarcseconds, corresponding to 33 light years, yields
a dynamical age for CORALZ J083139+460800 of 245
55 years. This confirms that CORALZ J083139+460800
is a genuinely young radio source, with an age similar to those found
for other compact symmetric objects (e.g. Polatidis & Conway 2003).
Table 1: Hotspot separation measurements and error estimates of CORALZ J083139+460800.
![]() |
Figure 1:
The measured hotspot separations in milliarcseconds for
CORALZ J083139+460800 from March 2000 to
March 2008, with the vertical bars indicating the 1- |
Open with DEXTER |
4.2 CORALZ J131739+411545
It is more of a challenge to analyse
CORALZ J131739+411545, because of its complicated
morphology. Its morphology was originally classified as complex
(CX), but considering all data currently available (Fig. 5),
we now tentatively classify this source as a compact double (CD) radio
source. The problem is that the strongest component is not resolved
well enough for detailed modelling of its structure, but too extended
to model it as a simple Gaussian component. Furthermore it is located
relatively close to the second component, making the construction of a
reliable radio map and the measurement of the component positions a
challenging effort. For this source JMFIT was run on each of the two
components separately, because otherwise the residuals of the strongest
component made the solution for the second component diverge. The
derived errors are crude estimates. On the one hand they are
overestimated because we estimate the centre of the components much
more accurately than the values of
indicate,
because these are dominated by symmetric residuals on either side of
the component. On the other hand the errors might be
underestimated, because the error estimates do not include possible
biases introduced during the synthesis of the map, which might
particularly arise in sources with a complex structure like this.
Table 2: Hotspot separation measurements and error estimates of CORALZ J131739+411545.
![]() |
Figure 2: The measured hotspot
separations in milliarcseconds for CORALZ J131739+411545
from March 2000 to March 2008, with the vertical bars
indicating the 1- |
Open with DEXTER |
![]() |
Figure 3:
Expansion velocities of compact symmetric objects in units of h-1 c, and radio luminosities at 5 GHz in
|
Open with DEXTER |
![]() |
Figure 4:
Contour maps of CORALZ J083139+460800, observed in
March 2000, in May 2004, in May 2006 with the VLBA as
part of VIPS (Helmboldt et al. 2007),
and in March 2008. Peak flux densities of each map are: 31.7,
38.9, 44.4, and 26.8 mJy. The rms noise levels are: 0.14,
0.14, 0.26, and 0.06 mJy. Contours are drawn at 3, 6, 9, 12,
and 24 mJy. All maps have been aligned and rotated, such that the
(lower) solid line, running through the main component's centre in each
map, is (horizontal). The dashed line runs parallel to this line,
while the (upper) solid line runs through the secondary component's
centre in each map. The lines slowly diverge, revealing a small but
significant increase in the separation between the two components. This
increase corresponds to a projected hotspot separation velocity of
0.095 |
Open with DEXTER |
![]() |
Figure 5: Contour maps of CORALZ J131739+411545, observed in March 2000 (top left), in May 2004 (top right), in May 2006 with the VLBA as part of VIPS (Helmboldt et al. 2007; bottom left), and in March 2008 (bottom right). Peak flux densities of each map are: 50.3, 63.8, 89.3, and 47.7 mJy. The rms noise levels are: 0.14, 0.14, 0.37, and 0.17 mJy. Contours are drawn at (-6, -3, 3, 6, 9, 12, 24, 48, 96, 200, 400, 800, 1600) times the rms noise level of the radio map. Negative contours are indicated with dashed lines. |
Open with DEXTER |
In Table 2 the
measured component separations for CORALZ J131739+411545
are presented, with in Col. 1 the observation date, in Col. 2
the post-fit rms error of the subimage containing the strongest
component, in Col. 3 the post-fit rms error of the subimage
of the second component, in Col. 4 the measured hotspot separation
in milliarcseconds, and in Col. 5 the 1- uncertainty of the hotspot separation. Figure 2
shows the measured hotspot separations in milliarcseconds for
CORALZ J131739+411545 from March 2000 to
March 2008, with the vertical bars indicating the 1-
uncertainties.
From this figure the error estimates appear reasonable. We perform a
linear least-squares fit to the data, resulting in a slope of 14
10 microarcsec per year, which corresponds to a 1-
upper limit to the projected expansion rate of
0.074 h-1 c = 0.106 c at a redshift of z=0.066
for CORALZ J131739+411545. This limit, together with the
measured projected hotspot separation of 3.3 milliarcsec,
corresponding to 14 light years, yields a lower limit to the
dynamical age of CORALZ J131739+411545 of
130 years. This result is in good agreement with the upper limit
of 4300 years, derived from the spectral ageing analysis of
CORALZ J131739+411545 (Mack et al. 2009).
4.3 Expansion velocity and radio power
The new results described above significantly increase the sparsely
available data on expansion velocities of young radio sources.
It is therefore important to review the correlation between
expansion speed and luminosity found by de Vries et al. (2009).
The expansion speed we determine here for
CORALZ J083139+460800, exceeds the previously derived
upper limit, but is still the lowest expansion speed measured in any
GPS radio source so far. Being a relatively weak radio source (
,
compared to a median radio luminosity of
for
sources with previously determined expansion speeds), this result
supports the correlation between the expansion speed and radio
luminosity. Further support comes from
CORALZ J131739+411545 (
), with an upper limit to the expansion speed lower than the value found for any radio source with a higher radio luminosity.
In Fig. 3 we show
the expansion velocities measured in young radio sources as a function
of radio luminosity. Both the radio luminosity and the lobe-to-lobe,
projected expansion velocity have been corrected for cosmological
effects, and represent rest-frame quantities. The arrows indicate upper
limits. The one indicated with the cross is the upper limit to the
expansion velocity of CORALZ J131739+411545 presented in
this paper. In de Vries et al. (2009)
we presented an analytic model of a young, expanding radio source with
a constant jet power in a power-law density environment. This model
predicts that
.
The solid line represents such a relation, with arbitrary scaling. The
dashed line represents this relation scaled down by a factor of
,
such that 90% of a sample of randomly oriented sources that
intrinsically would follow the solid line, should be located above the
dashed line. A Spearman rank correlation analysis shows that,
including the measured expansion speed presented in this paper and
ignoring all upper limits, the correlation is significant at
a 99% confidence level (n=10,
).
Including the upper limits in the statistical analysis tends to
increase the significance of the correlation, regardless of the precise
way these upper limits are treated, since the upper limits of the
low-luminosity radio sources are so stringent. Although it is still too
early to draw definite conclusions concerning the deviations of the
data from the model, it seems the data indicate a relation less
steep than the analytic model, suggesting the density
of the radio source environment is shallower than the assumed power-law slope of
(
), since the exponent of the radio luminosity
in this relation depends on the power-law slope of the environment as
(de Vries et al. 2009).
The random orientation of the radio sources only introduces a vertical
scatter, and does not affect the slope of the correlation.
5 Conclusions
In this paper we present a new epoch of global VLBI observations of
sources from the CORALZ sample. For
CORALZ J083139+460800 this resulted in the detection of
its projected hotspot separation velocity of 0.095
0.020 h-1 c,
representing the lowest expansion speed measured so far. This
result, together with the measured hotspot separation, yields a
dynamical age of 245
55 years. Despite the complicated structure of the double radio
source CORALZ J131739+411545, we were able to derive an
upper limit to the expansion speed of this source of
0.074 h-1 c,
yielding a lower limit for the dynamical age of 130 years. Both
expansion speeds are very low compared to CSOs with higher radio
luminosities, confirming the dependence of the expansion speeds of
young radio sources on their radio luminosities. These new measurements
raise the significance of the correlation above 99%, and indicate
the slope of this relation might be less steep than our model,
suggesting the radio source environments are shallower than
.
The European VLBI Network is a joint facility of European, Chinese, South African and other radio astronomy institutes funded by their national research councils (proposal codes: GS016, GS021, GV019). The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc (proposal codes: GS016, GS021, GV019).
References
- Conway, J. E. 2002, New Astron. Rev., 46, 263 [NASA ADS] [CrossRef] [Google Scholar]
- de Vries, N., Snellen, I. A. G., Schilizzi, R. T., Mack, K.-H., & Kaiser, C. R. 2009, A&A, 498, 641 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Fanti, C., & Fanti, R. 2003, Radio Astronomy at the Fringe, ed. J. A. Zensus, M. H. Cohen, & E. Ros, 300, 81 [Google Scholar]
- Fanti, R., Fanti, C., Schilizzi, R. T., et al. 1990, A&A, 231, 333 [NASA ADS] [Google Scholar]
- Fomalont, E. B. 1999, Synthesis Imaging in Radio Astronomy II, 180, 301 [NASA ADS] [Google Scholar]
- Helmboldt, J. F., Taylor, G. B., Tremblay, S., et al. 2007, ApJ, 658, 203 [NASA ADS] [CrossRef] [Google Scholar]
- Komatsu, E., Dunkley, J., Nolta, M. R., et al. 2009, ApJS, 180, 330 [NASA ADS] [CrossRef] [Google Scholar]
- Mack, K.-H., Snellen, I. A. G., Schilizzi, R. T., & de Vries, N. 2009, Astron. Nachr., 330, 217 [NASA ADS] [CrossRef] [Google Scholar]
- McMahon, R. G., & Irwin, M. J. 1991, in Proc. Conf. Digitised Optical Sky Surveys, ed. H. T. MacGillivray, & E. B. Thomson (Dordrecht: Kluwer Acad. Publ.), 417 [Google Scholar]
- O'Dea, C. P. 1998, PASP, 110, 493 [NASA ADS] [CrossRef] [Google Scholar]
- O'Dea, C. P., & Baum, S. A. 1996, in GPS and CSS radio sources, ed. I. A. G. Snellen, R. T. Schilizzi, H. J. A. Röttgering, & M. N. Bremer, 241 [Google Scholar]
- Owsianik, I., & Conway, J. E. 1998, A&A, 337, 69 [NASA ADS] [Google Scholar]
- Peacock, J. A., & Wall, J. V. 1982, MNRAS, 198, 843 [NASA ADS] [Google Scholar]
- Polatidis, A. G., & Conway, J. E. 2003, Publ. Astron. Soc. Aust., 20, 69 [Google Scholar]
- Readhead, A. C. S., Taylor, G. B., Xu, W., et al. 1996, ApJ, 460, 612 [NASA ADS] [CrossRef] [Google Scholar]
- Snellen, I. A. G., Schilizzi, R. T., Miley, G. K., et al. 2000, MNRAS, 319, 445 [NASA ADS] [CrossRef] [Google Scholar]
- Snellen, I. A. G., Mack, K.-H., Schilizzi, R. T., & Tschager, W. 2003, Publ. Astron. Soc. Aust., 20, 38 [Google Scholar]
- Snellen, I. A. G., Mack, K.-H., Schilizzi, R. T., & Tschager, W. 2004, MNRAS, 348, 227 [NASA ADS] [CrossRef] [Google Scholar]
- White, R. L., Becker, R. H., Helfand, D. J., & Gregg, M. D. 1997, ApJ, 475, 479 [NASA ADS] [CrossRef] [Google Scholar]
- Wilkinson, P. N., Polatidis, A. G., Readhead, A. C. S., Xu, W., & Pearson, T. J. 1994, ApJ, 432, L87 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Table 1: Hotspot separation measurements and error estimates of CORALZ J083139+460800.
Table 2: Hotspot separation measurements and error estimates of CORALZ J131739+411545.
All Figures
![]() |
Figure 1:
The measured hotspot separations in milliarcseconds for
CORALZ J083139+460800 from March 2000 to
March 2008, with the vertical bars indicating the 1- |
Open with DEXTER | |
In the text |
![]() |
Figure 2: The measured hotspot
separations in milliarcseconds for CORALZ J131739+411545
from March 2000 to March 2008, with the vertical bars
indicating the 1- |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
Expansion velocities of compact symmetric objects in units of h-1 c, and radio luminosities at 5 GHz in
|
Open with DEXTER | |
In the text |
![]() |
Figure 4:
Contour maps of CORALZ J083139+460800, observed in
March 2000, in May 2004, in May 2006 with the VLBA as
part of VIPS (Helmboldt et al. 2007),
and in March 2008. Peak flux densities of each map are: 31.7,
38.9, 44.4, and 26.8 mJy. The rms noise levels are: 0.14,
0.14, 0.26, and 0.06 mJy. Contours are drawn at 3, 6, 9, 12,
and 24 mJy. All maps have been aligned and rotated, such that the
(lower) solid line, running through the main component's centre in each
map, is (horizontal). The dashed line runs parallel to this line,
while the (upper) solid line runs through the secondary component's
centre in each map. The lines slowly diverge, revealing a small but
significant increase in the separation between the two components. This
increase corresponds to a projected hotspot separation velocity of
0.095 |
Open with DEXTER | |
In the text |
![]() |
Figure 5: Contour maps of CORALZ J131739+411545, observed in March 2000 (top left), in May 2004 (top right), in May 2006 with the VLBA as part of VIPS (Helmboldt et al. 2007; bottom left), and in March 2008 (bottom right). Peak flux densities of each map are: 50.3, 63.8, 89.3, and 47.7 mJy. The rms noise levels are: 0.14, 0.14, 0.37, and 0.17 mJy. Contours are drawn at (-6, -3, 3, 6, 9, 12, 24, 48, 96, 200, 400, 800, 1600) times the rms noise level of the radio map. Negative contours are indicated with dashed lines. |
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.