A VLA-GMRT Look at 11 Powerful FR II Quasars
S. Vaddi, P. Kharb, R.A. Daly, C.P. O'Dea, S.A. Baum, D.K. Deo, T.C., Barbusca, C. Murali

TL;DR
This study uses VLA and GMRT radio observations to analyze 11 powerful FR II quasars, exploring their orientation, environment, and jet power, revealing complex morphologies and intrinsic effects influencing their radio properties.
Contribution
It provides new insights into the environmental and intrinsic factors affecting FR II quasars and estimates their jet kinetic power using low-frequency observations.
Findings
Radio-loud unification largely holds but environmental effects are significant.
Quasars are observed at smaller angles than radio galaxies, indicating orientation effects.
Some sources show complex lobe morphologies suggestive of restarted AGN activity.
Abstract
We present results from 1.4 and 5 GHz observations at matched resolution with the Karl G. Jansky Very Large Array (VLA) of 11 powerful 3C FR II quasars. We examine the 11 quasars along with a sample of 13 narrow-line FR II radio galaxies and find that radio-loud unification largely holds but environmental effects cannot be ignored. The radio core prominence, largest linear size, and axial ratio parameter values indicate that quasars are at relatively smaller angles compared to the radio galaxies and thus probe orientation. Lack of correlation between statistical orientation indicators such as misalignment angle and radio core prominence, and larger lobe distortions in quasars compared to radio galaxies suggest that intrinsic/environment effects are also at play. Some of 150 MHz observations with the TGSS-GMRT reveal peculiar lobe morphologies in these FR II sources, suggesting complex…
Click any figure to enlarge with its caption.
Figure 1
Figure 10
Figure 11
Figure 12
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7
Figure 8
Figure 9
Figure 13
Figure 14
Figure 15
Figure 16
Figure 17
Figure 18
Figure 19
Figure 20
Figure 21
Figure 22
Figure 23
Figure 24
Figure 25
Figure 26
Figure 27
Figure 28
Figure 29
Figure 30
Figure 31
Figure 32
Figure 33
Figure 34
Figure 35
Figure 36
Figure 37
Figure 38
Figure 39
Figure 40| Name | RA | DEC | z | log | S | S | S | LLS | Arm-length | Misalignment | ||
| J2000 | J2000 | M⊙ | (Jy) | (Jy) | (Jy) | (kpc) | Ratio | Angle (deg) | ( erg s-1) | |||
| (1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) | (11) | (12) | (13) |
| 3C14 | 00:36:06.5 | +18:37:59 | 1.469 | 9.4 | 1.970.03 | 0.680.01 | 13.512.03 | 201.9 | -2.650.05 | 1.59 | 9.3, 24.8 | 9.710.33 |
| 3C47 | 01:36:24.4 | +20:57:27 | 0.425 | 9.2 | 3.590.01 | 1.010.01 | 36.334.46 | 374.1 | -1.860.05 | 1.23 | 6.4 | 3.630.05 |
| 3C109 | 04:13:40.4 | +11:12:14 | 0.306 | 8.3 | 3.410.02 | 1.260.01 | 20.722.73 | 407.6 | -1.180.05 | 1.04 | 13.7 | 1.300.04 |
| 3C204 | 08:37:44.9 | +65:13:35 | 1.112 | 9.5 | 1.200.02 | 0.350.01 | 13.881.86 | 282.8 | -1.680.05 | 1.40 | 3.3, 3.4 | 5.290.24 |
| 3C205 | 08:39:06.4 | +57:54:17 | 1.534 | 9.6 | 2.290.08 | 0.550.02 | 21.533.20 | 130.5 | -2.330.05 | 1.17 | 0, 7.5 | 15.30.30 |
| 3C208 | 08:53:08.8 | +13:52:55 | 1.11 | 9.4 | 2.390.1 | 0.480.02 | 23.083.56 | 88.60 | -1.970.05 | 1.28 | 0 | 10.600.23 |
| 3C249.1 | 11:04:13.7 | +76:58:58 | 0.312 | 9.3 | 2.120.01 | 0.680.01 | 18.262.60 | 199.2† | -1.530.05 | 2.04 | 7.6 | 1.150.05 |
| 3C263 | 11:39:57.0 | +65:47:49 | 0.646 | 9.1 | 2.960.05 | 1.030.02 | 20.373.07 | 320.5 | -1.530.05 | 1.72 | 1.6 | 4.420.12 |
| 3C336 | 16:24:39.1 | +23:45:12 | 0.927 | 9.2 | 2.640.06 | 0.380.01 | 14.472.23 | 172.0 | -2.610.05 | 2.00 | 13.6 | 5.830.21 |
| 3C351 | 17:04:41.4 | +60:44:31 | 0.372 | 9.5 | 2.650.03 | 1.080.01 | 16.942.36 | 316.8 | -2.250.05 | 2.23 | 5.8, 17 | 1.500.06 |
| 3C432 | 21:22:46.2 | +17:04:38 | 1.785 | 10.1 | 1.620.07 | 0.340.01 | 14.642.30 | 109.8 | -2.870.05 | 1.33 | 11.9 | 11.700.38 |
|
Col.(1-3) give the source name and position. Col(4)- Redshift from NED. Col.(5)-Logarithm of black hole mass from McLure et al. (2006). Col(6-8)- Total flux density at 1.4 GHz, 5 GHz, and 150 MHz respectively in Jy Beam-1. Col(9)- Largest linear size which is the distance between the two hotspots, measured using IMDIST on 5 GHz images. Col(10)- Logarithm of the radio core prominence. Col(11)- Arm-length ratio which is the distance between the core to the hotspot, measured using TVDIST on 5 GHz images. Col(12)- Misalignment angle in degree. Where multiple hotspots are present, angle to each of the hotspots is mentioned separated by a comma. An angle of zero indicates that the hotspots and core all lie on a single line. 5 GHz images are used to get the angles. The uncertainty on the misalignment angles is about a tenth of a degree. Col(13)- Jet kinetic power.
† Since this source has only one hotspot, the LLS for this source is measured using TVDIST and choosing the end points by eye. |
||||||||||||
| Name | Side | S | ||
| (mJy beam-1) | (arcsec) | |||
| (1) | (2) | (3) | (4) | (5) |
| 3C14 | N1 | 81.90 | -0.980.16 | 9.440 |
| N2 | 59.00 | -1.150.18 | 6.960 | |
| S | 111.7 | -0.860.14 | 14.97 | |
| 3C47 | N | 20.6 | -0.770.15 | 31.19 |
| S | 163.0 | -0.380.16 | 38.2 | |
| 3C109 | N | 31.9 | -0.450.11 | 46.09 |
| S | 211.0 | -0.660.12 | 48.12 | |
| 3C204 | E | 63.90 | -0.720.12 | 18.13 |
| W1 | 85.60 | -0.620.10 | 13.01 | |
| W2 | 24.10 | -1.040.17 | 17.12 | |
| 3C205 | N | 104.8 | -1.220.19 | 7.230 |
| S1 | 211.3 | -0.940.15 | 8.480 | |
| S2 | 207.0 | -1.100.18 | 8.650 | |
| 3C208 | E | 1356. | -0.620.10 | 6.220 |
| W | 205.3 | -0.460.07 | 4.860 | |
| 3C249.1 | E | 83.20 | -0.750.12 | 7.620 |
| W | 85.80 | -0.810.13 | 15.51 | |
| 3C263 | E | 27.50 | -0.520.08 | 27.88 |
| W | 486.7 | -0.660.11 | 16.24 | |
| 3C336 | N | 110.7 | -1.360.21 | 15.01 |
| S | 65.90 | -1.370.22 | 7.530 | |
| 3C351 | N1 | 330.0 | -0.730.12 | 27.29 |
| N2 | 194.0 | -0.760.13 | 25.02 | |
| S | 3.366 | -0.910.16 | 60.79 | |
| 3C432 | N | 130.0 | -1.310.21 | 5.700 |
| S | 109.0 | -0.910.15 | 7.590 | |
| Col.(2): N, S, E, W denote north, south, east, and west hotspots respectively;N2, S2, W2 are the second brightest hotspots in the north, west, and south sides. Col.(3):Peak flux density of the hotspot at 5 GHz estimated using TVMAXFIT. The errors are of the order of a few per cent (5%). Col(4): Spectral index of the hotspot between 1.4 and 5 GHz corresponding to the peak flux density position at the hotspot. This is obtained using IMVAL on the spectral images. Col(5): Distance from core to the hotspot in arcsec obtained using TVDIST on 5 GHz intensity images. The uncertainties are 1%. | ||||
| Name | S | S | |
| (mJy) | (mJy) | ||
| (1) | (2) | (3) | (4) |
| 3C14 | 24.32.7 | 10.91.6 | 0.680.10 |
| 3C47 | 58.26.5 | 69.27.5 | +0.080.09 |
| 3C109 | 248.727.8 | 274.030.1 | 0.040.05 |
| 3C204 | 58.16.5 | 50.35.8 | 0.290.08 |
| 3C205 | 40.44.5 | 26.63.3 | 0.340.05 |
| 3C208 | 59.46.6 | 51.05.9 | 0.030.02 |
| 3C249.1 | 76.08.5 | 78.98.8 | +0.140.05 |
| 3C263 | 148.116.6 | 137.615.6 | +0.050.02 |
| 3C336 | 38.94.4 | 11.12.2 | 0.620.09 |
| 3C351 | 35.54.0 | 20.32.7 | 0.450.12 |
| 3C432 | 10.21.2 | 6.10.8 | 0.470.08 |
| Col(2-3): Core flux density at 1.4 and 5 GHz respectively. Col(4):Core spectral index between 1.4 and 5 GHz. This is obtained using TVWIN + IMSTAT on the spectral index image. | |||
| Name | Bmin | Pmin | Axial Ratio | ||
|---|---|---|---|---|---|
| N,S | N,S | G | 10-11dyne cm-2 | N,S | |
| (1) | (2) | (3) | (4) | (5) | (6) |
| 3C14 | -1.390.14, -1.390.15 | 2.40.3, 1.80.3 | 109.218.3, 95.218.3 | 59.34, 45.11 | 4.00.5, 8.51.4 |
| 3C47 | -1.280.18, -1.280.18 | 31.04.7, 35.24.9 | 9.51.7, 9.51.7 | 0.45, 0.44 | 1.00.2, 1.10.1 |
| 3C109 | -1.180.18, -1.570.22 | 7.70.8, 15.61.6 | 16.92.4, 16.32.4 | 1.43, 1.32 | 6.00.6, 3.10.3 |
| 3C204‡ | -1.420.17, -1.50.17 | 5.90.7, 4.90.6 | 35.67.0, 50.77.0 | 6.29, 12.76 | 3.10.4, 2.70.3 |
| 3C205 | -1.420.15, -1.130.11 | 2.10.3, 1.80.3 | 116.422.1, 110.522.1 | 67.44, 60.74 | 3.40.5, 4.70.8 |
| 3C208‡ | -0.750.09, -0.710.08 | 2.00.3, 1.30.3 | 79.77.9, 72.77.9 | 31.56, 26.28 | 3.10.4, 3.71.0 |
| 3C249.1‡ | …, -1.030.11 | …, 7.60.8 | …, 20.82.0 | …, 2.15 | …, 2.00.2 |
| 3C263‡ | -1.060.14, -0.880.09 | 12.21.5, 9.11.1 | 18.32.2, 15.02.2 | 1.67, 1.11 | 2.30.3, 1.80.2 |
| 3C336 | -2.010.2, -1.890.20 | 6.30.7, 4.90.6 | 93.725.6, 93.225.6 | 43.68, 43.22 | 2.40.3, 1.50.2 |
| 3C351 | -1.100.18, -0.920.18 | 14.31.7, 7.40.8 | 15.42.1, 7.22.1 | 1.19, 0.26 | 1.90.2, 8.20.9 |
| 3C432 | -1.200.13, -1.350.15 | 2.80.4, 2.60.4 | 103.516.2, 98.416.2 | 53.27, 48.16 | 2.10.3, 2.90.5 |
|
Col(1): Source name. Col(2): Spectral index of the lobe between 1.4 GHz and 5 GHz for the northern and southern lobes. Col(3): Lobe width in arcsec obtained using the 1.4 GHz images. Col(4): Magnetic field at minimum pressure condition obtained using Equation 4. Col(5): Minimum internal pressure in the lobe obtained using Equation 3. Col(6) Axial Ratios i.e., ratio of the lobe length (core-to-hotspot distance) to lobe width.
‡ The estimates are for the eastern and western lobe respectively. |
|||||
| Name | Jet Side | Bright HS | Ratio of HS | Shorter |
|---|---|---|---|---|
| Side | peak flux | Side | ||
| (1) | (2) | (3) | (4) | (5) |
| 3C14 | S | S | 1.36 | N |
| 3C47 | S | S | 7.91 | N |
| 3C109 | S† | S | 6.61 | S |
| 3C204 | W | W | 1.34 | W |
| 3C205 | S | S | 2.02 | N |
| 3C208 | E(?) | E | 6.62 | W |
| 3C249.1# | - | - | - | N |
| 3C263 | E | E | 17.70 | E |
| 3C336 | S† | N | 1.68 | S |
| 3C351 | N(?) | N | 98.03 | N |
| 3C432 | S(?) | N | 1.19 | N |
|
Not observed in the current data, but observed in other high-resolution VLA archival images (Giovannini et al., 1987).
# This is a hybrid structure and so jet cannot be defined. Jet features that are ambiguous are denoted with a ‘?’. |
||||
| Name | z | Arm-length | log Rc | LLS | S | S | Pmin(N,S) | Qjet |
| Ratio | (kpc) | Jy | Jy | 10-11dyne cm-2 | 1045 erg s-1 | |||
| (1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) |
| 3C6.1 | 0.804 | 1.15 | -2.90 | 241 | 3.8 | 17.86 | 15.70, 7.51 | 11.48 |
| 3C13 | 1.351 | 1.36 | -4.19 | 276 | 1.92 | 17.11 | … | 28.46 |
| 3C34 | 0.69 | 1.14 | -3.47 | 346 | 1.78 | 16.03 | 7.46, 2.81 | 6.74 |
| 3C41 | 0.794 | 1.11 | -3.40 | 204 | 3.51 | 12.16 | 15.10, 9.27 | 8.21 |
| 3C44 | 0.66 | 1.63 | -3.94 | 502 | 1.37 | 10.98 | 4.24, 5.27 | 4.82 |
| 3C54 | 0.8274 | 1.07 | -3.96 | 458 | 2.03 | 11.48 | 3.87, 3.16 | 8.21 |
| 3C114 | 0.815 | 1.09 | -2.53 | 451 | 1.01 | 11.56 | 5.36, 2.22 | 6.35 |
| 3C142.1 | 0.4061 | 1.87 | -3.57 | 299 | 3.10 | 28.48 | 4.88, 1.42 | 3.00 |
| 3C169.1 | 0.633 | 1.44 | -3.40 | 367 | 1.297 | 8.00 | 2.18, 1.21 | 3.12 |
| 3C172 | 0.5191 | 1.13 | -4.18 | 676 | 2.84 | 18.85 | 1.45, 0.65 | 4.54 |
| 3C441 | 0.707 | 2.45 | -4.32 | 299 | 2.53† | 14.56 | 1.49, 2.74 | 7.29 |
| 3C469.1 | 1.336 | 1.07 | -3.40 | 712 | 1.72 | 12.89 | 13.20, 4.15 | 25.79 |
| 3C470 | 1.653 | 1.36 | -3.97 | 255 | 1.72 | 11.92 | … | 37.56 |
|
Taken from Condon et al. (1998)
Col.(1)- Name of the source. Col(2)- Redshift from NED. Col(3)- Arm-length ratio. Col(4)- K-corrected radio core prominence parameter . The 1.4 GHz flux density is taken from White & Becker (1992) and the 5 GHz flux density is estimated from 8 GHz assuming . Col.(5)- Projected largest linear separation. Col.(6)- Total flux density of the source at 1.4 GHz taken from White & Becker (1992) in units of (Jy). Col(7)- 150 MHz flux density from TGSS-GMRT in units of (Jy). Col(8)- Minimum internal pressure in the lobe obtained using Equation 3. Col(9)- Jet kinetic power estimated as described in Section 4.1. |
||||||||
| Property 1 | Property 2 | Correlation | Probability |
|---|---|---|---|
| coefficient | p-value | ||
| Arm-length Ratio | 0.2 | 0.4 | |
| -0.4 | 0.003 | ||
| -0.5 | 0.02 | ||
| Arm-length Ratio | -0.4 | 0.05 | |
| LLS | -0.3 | 0.2 | |
| Axial Ratio | 0.3 | 0.01 | |
| Misalignment Angle | 0.03 | 0.8 |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
A VLAGMRT Look at 11 Powerful FR Quasars
S. Vaddi,1 P. Kharb, 1 R. A. Daly, 2 C. P. O’Dea, 3,4 S. A. Baum, 3,4 D. K. Deo, 5,7 T.C. Barbusca, 2 C. Murali 6,7
1National Centre for Radio Astrophysics - Tata Institute of Fundamental Research, Ganeshkhind, Pune 411007, India
2Penn State University, Berks Campus, P. O. Box 7009, Reading, PA 19610-6009
3Physics and Astronomy, University of Manitoba, Winnipeg, MB R3T 2N2 Canada
4Rochester Institute of Technology, Rochester, New York, 14623, USA
5University of Missouri-Kansas City, Kansas City, Missouri, 64110, USA
6University of Texas at Dallas, Texas, 75080, USA
7Indian Institute of Astrophysics, II Block, Koramangala, Bangalore 560034, India E-mail: [email protected] towards this project was performed at this institute.
(Accepted XXX. Received YYY; in original form ZZZ)
Abstract
We present results from 1.4 and 5 GHz observations at matched resolution with the Karl G. Jansky Very Large Array (VLA) of 11 powerful 3C FR quasars. We examine the 11 quasars along with a sample of 13 narrow-line FR radio galaxies and find that radio-loud unification largely holds but environmental effects cannot be ignored. The radio core prominence, largest linear size, and axial ratio parameter values indicate that quasars are at relatively smaller angles compared to the radio galaxies and thus probe orientation. Lack of correlation between statistical orientation indicators such as misalignment angle and radio core prominence, and larger lobe distortions in quasars compared to radio galaxies suggest that intrinsic/environment effects are also at play. Some of 150 MHz observations with the TGSS-GMRT reveal peculiar lobe morphologies in these FR sources, suggesting complex past lives and possibly restarted AGN activity. Using the total 150 MHz flux density we estimate the time-averaged jet kinetic power in these sources and this ranges from erg s*-1*, with 3C 470 having the highest jet kinetic power.
keywords:
quasars: general — radio continuum: galaxies
††pubyear: 2018††pagerange: A VLAGMRT Look at 11 Powerful FR Quasars–A
1 Introduction
Active galactic nuclei (AGNs) are the centres of a special class of galaxies that consist of actively accreting central supermassive black hole (SBH). High-speed collimated outflows are present in a small fraction of AGNs. AGNs show a dichotomy in their radio power - these are referred to as radio-loud (RL) and radio-quiet (RQ) AGNs (Strittmatter et al., 1980; Miller et al., 1993). Historically, AGNs with higher radio-to-optical flux where are termed as RL AGNs and RQ otherwise (Kellermann et al., 1989). This distinction however, is suggested to be valid only for broad-lined, unobscured AGNs since AGN luminosity estimates may be affected by dust, and therefore Padovani (2017) suggested a nomenclature based on a fundamental physical difference, namely the presence or absence of relativistic jets.
Based on the radio morphology and power, extended RL AGNs are classified into FR and FR (Fanaroff & Riley, 1974). FR sources have radio structures consisting of a core, collimated jets, lobes, and hotspots located at the edges of these lobes. FR s show extended plume-like structures and tails with no distinct collimated jets or terminal hotspots (more detailed differences in Bridle & Perley, 1984). Based on the optical spectra RL AGNs are classified into Type 1 (have broad and narrow emission lines) and Type 2 (have only narrow emission lines). Type 2 AGNs comprise narrow-line radio galaxies (NLRG); Type 1 AGNs comprise broad-line radio galaxies (BLRG) at low luminosity and radio-loud quasars at high luminosity. The latter are further divided into flat spectrum radio-loud quasars (FSRQ) () and steep spectrum radio-loud quasars (SSRQ) () (see Urry & Padovani (1995); Tadhunter (2008); Netzer (2015); Tadhunter (2016) for a review). SSRQs tend to have lobe-dominated radio structure while FSRQs often have core-dominated structure. Although we see different types of AGN, it is hypothesised that these are intrinsically similar objects but appear different depending on their orientation with respect to the line of sight. This has come to be known as orientation-based AGN Unification Scheme (Antonucci, 1993; Urry & Padovani, 1995).
There is compelling evidence in the literature that is consistent with the predictions of the unification scheme. Some of the evidences include detection of broad emission lines in the polarized intensity spectra of Type 2 AGNs indicating a hidden broad line region obscured from direct view (Antonucci, 1984; Antonucci & Miller, 1985; Tran et al., 1995; Young et al., 1996; Ogle et al., 1997; Cohen et al., 1999), detection of ionisation cones of the narrow line region in optical images of Type 2 AGNs (Pogge, 1988; Wilson et al., 1993; Jackson et al., 1998), detection of excess near-IR emission and polarized light perpendicular to the jet axis (Antonucci & Barvainis, 1990), apparently superluminal velocities in VLBI observations of compact radio sources (Cohen et al., 1977; Kellermann et al., 2007), detection of extended halos around compact sources representing radio lobes viewed face-on (Browne et al., 1982; Antonucci & Ulvestad, 1985), and observed anisotropy in the radio structure (Laing et al., 1983; Barthel, 1989; Smith & Heckman, 1990; Mullin et al., 2008). According to unification, FSRQ form the beamed counterparts of FR II radio galaxies, and SSRQ (and BLRG) are oriented at angles that are in between NLRG and FSRQ.
This paper is the first in a series of papers in our study of FR II radio-loud AGNs. Here we present the results from our VLA study of 11 SSRQ that show classical double lobes at GHz frequencies. We have analysed the properties of these quasars along with a sample of 13 FR radio galaxies from Kharb et al. (2008). Interestingly, several of these quasars and radio galaxies reveal extended, winged or asymmetric radio morphologies in the TIFR GMRT Sky Survey (TGSS) at 150 MHz. This is in contrast to the symmetric structures detected with the VLA at GHz frequencies.
Throughout this paper, we assume a cosmology with H0 = 73 km s*-1* Mpc*-1*, = 0.27, = 0.73. Spectral index, , is defined such that flux density at frequency is .
2 The Sample
The 11 quasars in this study were selected to augment the sample of 55 radio sources with spin determinations described by Daly (2011). The parent population is powerful FR II sources that are identified as quasars, where the original published optical classification was used to identify them as quasars, and this identification was confirmed by the follow-up observations of McLure et al. (2006); these 11 sources were identified based on the following criterion: (i) the quasar has a black hole mass determination (from McLure et al. (2006)); (ii) at least one side of the source does not exhibit a radio jet, so that a spectral ageing analysis can be carried out on the non-jetted side; (iii) the non-jetted lobe(s) have angular sizes large enough to allow the spectra of the lobe to be measured at several locations.
2.1 Observations and Data Reduction
These 11 quasars were observed with the VLA at 1.4 GHz (L-band) and 5 GHz (C-band) with A and CnB-array configurations, respectively (Project ID: 10C-178) for a total of 22 hours. 5 GHz observations were not scheduled for 3C 336111One of the observations (for 3C432) had technical issues and was therefore re-observed. However, due to confusion in the sequencing of observations, we missed obtaining data for 3C336. This was realised late and could not be rectified. Hence, we used archival VLA data222Project ID: AB0454; C-band B-array configuration for this source. Also, due to the noisy 1.4 GHz data on 3C 351, we chose archival data333Project ID: AD0429; L-band A-array configuration for this source as well. For 3C 208 at 5 GHz, since there were issues with amplitude calibration using 3C 48, the source flux was rescaled using the aips task RESCALE and the core flux density from Yuan & Wang (2012).
The data in the SDMset format were pre-calibrated using the casa calibration pipeline444casa release 4.5.3 with pipeline release 4.6.0. The final imaging and self-calibration were done in the Astronomical Image Processing System (aips), by iteratively running the tasks IMAGR and CALIB. The 1.45 GHz spectral index images were made using task COMB after first creating L and C-band images with identical circular beams (at the poorer L-band resolution), which were then positionally aligned using the AIPS task OGEOM. We used 3 r.m.s. noise of the respective L and C-band images as the clipping level input to COMB. Final flux density and spectral index values were estimated using the AIPS verbs TVWIN and IMSTAT, as well as the Gaussian-fitting task JMFIT. All extents were measured in AIPS using IMDIST. The basic properties and derived parameters are listed in Tables 1 through 2.
2.2 Data from the Literature
In order to put the results from the 11 quasars in perspective, we have analysed their results along with the sample of 13 FR narrow-line radio galaxies that were studied in Kharb et al. (2008); the radio galaxy study had been carried out at similar frequencies and resolutions as the 11 quasars presented in this paper. Some properties of the radio galaxies have been re-estimated (like , 1.4 GHz flux density from NVSS) to compare them with the current quasar sample, and are tabulated in Table 6. We note that the combined quasar and radio galaxy sample studied here is statistically small and eclectic. Some statistical tests, therefore, suffer from small number statistics.
2.3 Physical Parameters Estimation Method
The spectral index is estimated from the spectral index images. The spectral index of the core is obtained by placing a box around the core using TVWIN and estimating the mean using IMSTAT. The hotspot spectral index is measured at the image position corresponding to the peak intensity in 1.4 and 5 GHz image, where the peak intensity is determined using TVMAXFIT.
The arm-length ratio is defined as the ratio of the longer core-hotspot distance to the shorter core-hotspot distance. The distance is measured from the peak flux position at the core to the peak flux position at the hotspot. The peak flux density was estimated using the procedure TVMAXFIT. The 5 GHz radio maps were used in the determination of the arm-length ratio.
The axial ratio is defined as the ratio of the length to the width of the lobe (Leahy et al., 1989). The core-hotspot distance is taken as the length of the lobe. Although lobe emission is not observed to be extending all the way to the core in the 1.4 GHz emission, lower frequency images (610 MHz GMRT observations, Vaddi et al. 2018 in preparation) show extended emission and thus it is safe to assume the lobe length to be the core-hotspot distance. The FWHM of the lobe is measured using the Gaussian-fitting task SLFIT to a one-dimensional slice across the lobe. The slice is chosen by-eye to be perpendicular to and at half the distance between the core and hotspot. To account for the finite resolution of the observing beam, we subtract the beam in quadrature from the measured Gaussian fit: where is the measured FWHM from the Gaussian fit, and is the average FWHM of the observing beam (), and is the FWHM of the Gaussian fit deconvolved with the beam. The lobe width is then taken as 2/ of the deconvolved Gaussian fit (Leahy et al., 1989). Correction for the radio galaxies is negligible since its lobe widths are large compared to the beam size (Daly et al., 2010).
The largest linear size (LLS) is defined as the distance between two hotspots (definition consistent with Kharb et al. (2008)). The distance is estimated by first obtaining the position of the hotspot peak intensity with TVMAXFIT and then measuring the distance between the hotspot peak intensity using IMDIST. The 5 GHz maps were used in the estimates.
The misalignment angle is the difference between the position angles of the core-hotspot segments on either side of the core. When there is more than one hotspot on one side of the core, the other side is used to define the core-hotspot line segment and misalignment angles to each of the two hotspots on the opposite side are obtained. Figure 27 in the appendix identifies different components that were estimated using AIPS task TVDIST.
2.4 Errors
Uncertainties in the flux measurements were estimated by combining all errors in quadrature. The errors considered are the calibration errors, rms, and random errors (). We assumed that the calibration errors are 10 per cent of the measured flux densities. RMS is the error reported in the tasks (JMFIT, IMSTAT) used for flux measurement. Random error of 5 per cent accounts for the different values we get when taking repeated measurements; an average of five measurements were obtained for a few sources and the per cent standard deviation from the average values was noted. This gives the uncertainties of the order of few mJy beam*-1*.
The error in the spectral index is estimated by adding in quadrature the noise map produced in the COMB task and a 5 per cent random error for , 15 per cent random error for the , and a 10 per cent random error for the .
The uncertainty on the measurement of the lobe width is obtained by adding in quadrature the 5 per cent random error on the measured width and the error in the Gaussian fit to the slice. The uncertainty of the deconvolved lobe width is obtained using the expression where is the uncertainty of the lobe width and is the uncertainty of the observing beam. Lobes whose widths are similar to the beam have larger uncertainty.
The uncertainties on angular size, arm-length ratios are estimated by repeated measurements (five times) on few sources and taking the standard deviation of these measurements. The ratios are calculated using the usual error propagation. These uncertainties are between 520 per cent.
The LLS is estimated from peak hotspot-to-hotspot distance which need not be the true LLS and true LLS is usually larger than the hotspot-hotspot distance. To quote the uncertainty in the LLS due to this, we estimate an average of several lobe edge-to-edge distance measurement obtained using TVDIST and the deviation of this estimate with the hotspot-hotspot distance is taken as the uncertainty. The uncertainty of the misalignment angle is estimated based on the uncertainties of the hotspot and core positions.
3 Results
In our quasar sample, 70 per cent of sources show one-sided radio jets, while the remaining 30 per cent sources do not clearly show a jet. In all but one source (3C 249.1 which shows a hybrid FR I/FR II morphology) two-sided hotspots located at the end of the radio lobe (i.e. FR II type morphology) are seen. For two sources (3C 109, 3C 336), jet-like features are seen in higher resolution images in the literature (Giovannini et al., 1987). In the radio galaxy sample of Kharb et al. (2008), 30 per cent of the sources show one-sided jet-like features (see Table 6 of (Kharb et al., 2008)). The 5 GHz images of the quasars show the presence of multiple hotspots in some sources. The southern lobe of 3C 263 and northern lobe of 3C 351 have off-axis extended emission i.e., the emission is extended on one side of the jet axis. See the Appendix for the contour images of the quasars at 1.4 GHz and 5 GHz.
We studied various properties of the cores, lobes, and hotspots of quasars and radio galaxies. We have quantified the observed trends with appropriate statistical tests. The significance of the relationship between properties is assessed using the non-parametric Spearman and Kendall rank tests. Both tests produced nearly identical results and so we have quoted here the results from the Spearman rank test. The p-value is a measure of the significance of the null hypothesis that two variables are uncorrelated. To determine if two distributions are the same, we have used the two sample KolmogorovSmirnov (KS) test. The null hypothesis to be tested in the KS test assumes that both the samples are drawn from the same population. A p-value of 0.05 and lower is considered to be significant to reject the null hypothesis. gives the per cent confidence level at which the null hypothesis is rejected.
3.1 Structure Asymmetry
Arm-length ratio is an asymmetry parameter and is defined as the ratio of the longer to the shorter arm-length, where arm-length is the core-to-hotspot distance. For quasars that show two hotspots on one side, we measured the arm-length ratio using both the hotspots. Figure 1 shows the arm-length ratio with respect to the redshift for the quasars and FR radio galaxies. The arm-length ratio distribution for quasars shows that there are both symmetric and asymmetric structures. A larger fraction of the sources in the radio galaxy sample are symmetric sources than in the quasar sample. Combining both the samples, we find that per cent of the sources have arm-length ratios lower than the combined average arm-length ratios (=1.49) suggesting symmetric structures. The KS test indicates that there is no substantial difference in arm-length ratios between the two types of AGNs and they may be drawn from the same population (). The arm-length ratios appear to be uncorrelated with redshift () (see also Kharb et al., 2008) suggesting that symmetric and asymmetric structures prevail at all redshifts.
The (brighter) jet side is observed in seven out of 10 quasars (3C 249.1 is excluded because of its hybrid structure). The remaining three sources have ambiguity in the jet feature (see Table 5). Of these seven jetted quasars, four quasars have their bright hotspot side on the same side as the shorter lobe. If we also include the “tentative” brighter jet sides from our images, 50 per cent (i.e. five out of ten quasars) have their brighter hotspot side on the same side as the shorter lobe. Simple relativistic beaming theories, however expect the longer arm to be the brighter arm (Ryle & Longair, 1967). Our results suggest that the brighter hotspot side is equally likely to be either on the shorter arm or on the longer arm . This asymmetry in brightness and morphology may be a result of asymmetric environments, rather than orientation-dependent Doppler effects. The environmental effects that constrained the jet from expanding might also make the hotspot appear brighter (Saikia et al., 1995; Gopal-Krishna & Wiita, 2004).
3.2 Hotspot Spectral Indices
The correlation of spectral index with redshift or radio luminosity is well established (Tielens et al., 1979; Blumenthal & Miley, 1979; Wellman et al., 1997; Dennett-Thorpe et al., 1999; Ishwara-Chandra & Saikia, 2000; Kharb et al., 2008). The tendency for the integrated spectral index to steepen with redshift is conventionally attributed to a “k-correction”. At higher redshift, the spectrum gets shifted to lower frequencies bringing the steep part of the spectrum closer to the observed frequency (Bolton, 1966; Blumenthal & Miley, 1979; Laing & Peacock, 1980; Gopal-Krishna, 1988). This interpretation for the is based on the observation that the SEDs of radio galaxies steepen at high frequencies. However, the k-correction mechanism has been disputed owing to the remarkable linearity and the lack of evidence for a spectral curvature in wide band radio spectrum of high redshift radio galaxies (Chambers et al., 1990; Klamer et al., 2006; Klamer et al., 2007). Gopal-Krishna et al. (2012) however argue that since the steepness of the spectral index anti-correlates with the spectral curvature (Mangalam & Gopal-Krishna, 1995), the spectral steepening becomes negligible for sources with steep spectra (, also reflected in Klamer et al. (2006)). k-correction being insufficient to explain the observed relation, alternative effects have been proposed. A density-dependent effect in which the jets are working against a denser environment at higher redshift can result in greater synchrotron losses at higher frequencies (Athreya & Kapahi, 1998; Klamer et al., 2006; Klamer et al., 2007; Ker et al., 2012). It has also been suggested that for high- sources, there can be an enhanced spectral aging due to increased inverse-Compton losses of the relativistic electrons against the cosmic microwave background (CMB), since the CMB energy density increases as (Gopal-Krishna et al., 1989; Krolik & Chen, 1991; Morabito & Harwood, 2018). High power sources have higher magnetic fields in the hotspots than in the less powerful sources. This can result in higher synchrotron losses, thus steepening the spectra (Meisenheimer et al., 1989; Gopal-Krishna & Wiita, 1990; Blundell et al., 1999). Yet another possibility is that the is an indirect manifestation of an intrinsic correlation between the luminosity and spectral index which when coupled to Malmquist bias translates to the correlation (Laing & Peacock, 1980; Chambers et al., 1990; Blundell et al., 1999). This, however, seems to be not the case in our sample as we discuss ahead.
Unlike previous studies where integrated spectral index was used, we examine the hotspot spectral index. Since hotspots are located far away from the galaxy, their study can provide physical insight into the environments of the radio source, and thereby test the effect of ambient density on source properties. Figure 2 shows the variation of spectral index of the hotspot with redshift in quasars and radio galaxies. We see a spectral steepening in the quasar and radio galaxies with increasing redshift consistent with earlier findings. The Spearman rank test indicates a statistically significant correlation with the correlation coefficient of -0.42 at a significance level greater than 99.99 per cent. A stronger dependence of spectral index on the luminosity is also observed in our sample and is shown in Figure 3. High radio power sources have steeper spectral indices and follow the relation at a significance level of 99.99 per cent. However, this correlation turns out to be not significant when the partial correlation 555The ppcor package in R programming language has been used in partial correlation statistics. coefficient (, ) was measured considering the effect of redshift. The observed correlation is likely due to the strong correlations between (correlation coefficient = 0.7) and (correlation coefficient = 0.5) at > 99.9% significance.
The average spectral index of the hotspots in quasars is 0.880.04 similar to the radio galaxies which has =0.870.02 . The KS test indicates that the distribution of quasars and radio galaxies are the same (). In 80 per cent of the quasars, the hotspot with the flatter spectral index is also the one with higher peak power at 1.4 GHz; this fraction is per cent in the case of the radio galaxies. In 70 per cent of the quasars (i.e., 7 out of 10 quasars excluding the hybrid quasar 3C349.1) where one-sided jets are observed, the jet-sided lobe has a significantly brighter hotspot than the non-jetted side in 6 out these 7 quasars. On the other hand, in the radio galaxies, only per cent of them show one-sided jet-like features and again the jet side points to the brightest hotspot. These observations (i.e. a similar distribution of hotspot spectral indices, jet-side pointing to the brightest hotspot in quasars and radio galaxies) imply that the brighter hotspot is likely, not due to Doppler boosting. The hotspots of quasars and radio galaxies considered here, do not have a bulk motion that is relativistic. The hotspot power is between 10-20 per cent of the total power at 1.4 GHz. The hotspot power correlates and increases with the total power in quasars and radio galaxies. The statistical significance of the relation for the hotspot () is stronger than the relation for core spectral index ().
3.3 Radio Core Prominence
The radio core prominence , is the ratio of the (beamed) radio core flux density to the (unbeamed) lobe flux density (Orr & Browne, 1982; Zirbel & Baum, 1995). It is a statistical indicator of beaming and thereby orientation (Padovani & Urry, 1992; Kharb & Shastri, 2004); larger values correspond to smaller viewing angles. The k-corrected radio core prominence is calculated using
[TABLE]
where is the core flux density at 5 GHz and is the total radio emission at 1.4 GHz. Figure 4 shows log with respect to redshift for quasars and radio galaxies. For the radio galaxies, the 5 GHz core flux density is estimated from the 8 GHz measurements assuming a spectral index of (average core spectral index in Kharb et al. (2008) sample), and the total flux density is taken at 1.4 GHz. It can be seen that quasars have larger values when compared to the radio galaxies. This is consistent with the orientation-based unification of quasars with radio galaxies; quasars are oriented at relatively smaller angles compared to the radio galaxies.
The median log for the quasars and radio galaxies is and , respectively, consistent with the RL unification scheme. For quasars, , and for the radio galaxies, . Kharb & Shastri (2004) (see their Figure 3) have estimated log values for a large number of AGN ranging from the plane-of-sky narrow-line radio galaxies to nearly pole-on blazars. Using the relations between , jet speed and orientation angles presented in Appendix A of Kharb & Shastri (2004), a minimum log () value of (see Kharb & Shastri, 2004), bulk jet Lorentz factor of 10 (Kapahi & Saikia, 1982), we deduce that the radio galaxies considered here lie at orientation angles between 341 and 742, while the quasars lie at orientation angles between 191 and 381. 85% of the radio galaxies have orientation angles greater than 45. Barthel (1989) suggested that the average for quasars, and for radio galaxies. These are, of course, crude representative estimates because is a statistical indicator of orientation and cannot be used to derive orientation angles for individual sources. It is important to note that while quasars are at relatively small angles compared to the radio galaxies, as indicated by their relative values, the overlapping spread in values indicates that the quasars are not oriented at small enough angles () to qualify as blazars.
Figure 5 shows log plotted against the arm-length ratio, Q. We find a weak anti-correlation in the combined sample with a correlation coefficient of at a significance of 95 per cent (). This marginal relation, which also has a large scatter, might suggest that the arm-length ratio is not strongly affected by projection but is an intrinsic asymmetry. We suggest that this weak relation is apparent and may be a result of small sample size, and and arm-length ratio may not be correlated. This is supported by the orientation independent statistics for the arm-length ratio in this sample (see Section 3.1).
3.4 Largest Linear Size
The jet orientation angle with respect to our line of sight influences the observed size of a radio source such that, . Sources with jets oriented close to the line of sight will have smaller apparent sizes than sources oriented close to the plane of the sky for sources with similar intrinsic sizes. In the present work, the angular sizes of the radio galaxies appear to be larger than the quasars. The average projected LLS of radio galaxies is 1.6 times larger than those of quasars. Following the Unification scheme, if we assume that quasars and radio galaxies are the same objects and have similar distributions in true linear sizes, then . Using an average and for quasars and radio galaxies, respectively, as indicated by their values, gives , fully consistent with .
The LLS ratio can be used to estimate the angle dividing quasars and radio galaxies () using the foreshortening calculations given by:
[TABLE]
For an LLS ratio of 1.6, this gives . These are in agreement with the torus opening angle estimated from the projected linear sizes of the radio sources. The average projected LLS of quasars is 24028 kpc, while the radio galaxies have an average projected LLS of 390 kpc. Using Eq (7) of Arshakian (2005), the torus opening angle for the quasars in our sample is 578. This yields an average orientation angle of 74 for radio galaxies and 39 for quasars. Encouragingly, this matches with the upper end of the average orientation angles estimated using Rc.
Wan & Daly (1998) and Barthel (1989) reported ratios of 1.25 and 1.8 respectively for 3C sources. Recently, using galaxies from the Bootes survey observed using LOFAR, Morabito et al. (2017) found the ratio of the average projected linear size of radio galaxies to the quasars to be , and after correcting for redshift evolution. The LLS estimates in this study are consistent with the LLS estimates in the literature, and minimal difference in the LLS ratio is likely due to the way the samples are selected. The selection criteria in our study (and Wan & Daly (1998)) favours quasars that have well-resolved extended lobes which in turn correspond to quasars at relatively larger viewing angles to the line of sight.
Overall, our results are consistent with the orientation-based unification scheme for quasars and radio galaxies. However, the range of orientation angles might not be large enough to show clear correlations with respect to orientation-indicators like or misalignment angles. We find the inverse dependence of the angular size with the redshift (widely known as the D-z relation) consistent with previous studies (Miley, 1968; Kapahi, 1985; Blundell et al., 1999; Ker et al., 2012; Morabito et al., 2017).
The dependence of the GHz spectral index of the hotspots on the projected core-hotspot distance is shown in Figure 6. We have plotted the spectral index of the hotspots of both the lobes for each source. We find a correlation between the spectral index and the projected core-hotspot distance of the quasars at a 96 per cent significance level. This apparent correlation, however, disappears once a partial correlation test is performed between and the core-hotspot distance neglecting the effects of redshift. This is indicative of the evolution of spectral index and size with redshift.
3.5 Lobe Axial Ratios
The axial ratios (lobe length-to-width ratio) for quasars are lower and statistically different from the radio galaxies in our study consistent with results obtained by Leahy et al. (1989). The intensity maps show more asymmetric and distorted lobes (such as extended diffuse emission to one side of the jet) in quasars (e.g. 3C 14, 3C 204, 3C 249.1, 3C 263, 3C 351) than in radio galaxies. The KS test suggests that quasars and radio galaxies have different distributions at 96 per cent significance (). The median axial ratios of quasars and radio galaxies are 3.0 and 4.2 respectively. The lobe widths of quasars are on an average only 1.3 times fatter (i.e. larger lobe width) than the radio galaxies and their median lobe widths are similar. Thus, lobe widths of quasars studied here are rather insensitive to projection effects, consistent with Wan & Daly (1998) who showed that projection has little effect on lobe width except at very small angles or for source widths measured close to the hotspot. This indicates that the ratio of lobe length to the width i.e. axial ratio is likely governed by the projection of the length. The average core-hotspot size ratio for similar lobe widths is 1.45. This likely suggests that the different axial ratios can be accounted for by projection effects alone. Leahy et al. (1989) however find that the differences in axial ratios cannot be explained by projection effects alone. Larger lobe structure distortions seen in quasars (also seen in our quasar sample) can partly explain the low axial ratios. Such lobe distortions can, for example, arise due to variations in the outflow angles, as discussed in Section 3.6. Any evolutionary relationship between quasars and radio galaxies leading to smaller linear sizes in quasars compared to radio galaxies, however, cannot be ruled out (e.g. Miley, 1971; Wardle & Miley, 1974).
The axial ratios show a marginal correlation with the spectral index of the lobes (correlation coefficient = 0.3, ) (see Figure 7) and an anti-correlation with the lobe flux density (=-0.4, ), both at a greater significance level i.e. fatter lobes are brighter and have steeper spectra. This is believed to be due to confinement of the lobes which enhances the radio emission through synchrotron losses while reducing the adiabatic expansion losses (e.g., Roland et al., 1985; Barthel & Arnaud, 1996). We further compare the axial ratios against the structure asymmetry using the arm-length ratio in Figure 8. Individually, the quasars and radio galaxies do not show any correlation, but the combined sample shows an anti-correlation at a significance level of 96 per cent (). This anti-correlation implies that sources that are narrower i.e. with larger axial ratios are also more symmetric in terms of their lobe sizes (i.e., have arm-length ratios closer to unity). The shorter or fatter sources are the more asymmetric. Such an anti-correlation was also observed in the combined radio galaxy sample in Kharb et al. (2008), who found that the shorter lobe was fatter and had a steeper spectral index.
3.6 Jet Axis Misalignment
Misalignment angle is the complement of the angle between the lines joining the core to the two hotspots. This is an asymmetry parameter and is sensitive to the orientation of the source to the line of sight. An intrinsic jet axis misalignment can appear enhanced for small angles to the line of sight, due to projection effects (Reynolds, 1980; Conway & Murphy, 1993). However, the misalignment angle is also found to be sensitive to environmental effects (Kharb et al., 2008). The misalignment angles for quasars are typically greater than and the average misalignment angle for quasars is about 1.8 times greater than radio galaxies (See Table 1), however, we cannot reject the hypothesis that both the samples have similar distribution at 28 per cent (). Although differences in misalignment angles can become amplified due to a geometrical projection effect, the difference seen in our sample is marginal because amplification is likely to be significant only in sources that are oriented close to the line of sight (Moore et al., 1981; Saikia et al., 1995).
Figure 9 shows the misalignment angle versus the radio core prominence, log . Statistical tests show no correlation for the quasars and radio galaxies taken together () or separately for quasars (). Kharb et al. (2008) had found a weak positive correlation for a large sample of radio galaxies. The radio core prominencemisalignment angle correlation has also been observed in quasars by Kapahi & Saikia (1982); Hough & Readhead (1989). This implies that either or the misalignment angles or both are not uniquely orientation-dependent. Environmental factors may be instrumental in this lack of a correlation. could be influenced by environmental asymmetries on parsec-scales, while misalignment angles could be influenced by environmental asymmetries on parsec as well as kpc-scales (Kharb et al., 2010).
Figure 10 shows the misalignment angles against axial ratios (longer and narrower sources have larger axial ratios). Quasars show smaller axial ratios and larger misalignment angles. This is consistent with the idea that environment plays an important role in tracing the jet and lobe structure. When the jet and the counter-jet plough into a medium with different densities, it will cause the jet to bend leading to greater misalignment angles (McCarthy et al., 1991). It is also possible that outflow angle is changing with time, as in the \qqdentist drill model (Scheuer, 1982) the jet axis random walks sending out the jet in different directions.
3.7 Minimum Internal Pressures
The internal pressure of the jet when compared with the external gas pressure of the ambient medium can be used to understand whether the ambient medium is sufficient to confine the outer parts of the radio structure. The non-thermal internal pressure in the radio jets contributed by the relativistic particles and magnetic field can be crudely estimated from the minimum energy arguments for synchrotron radiation. It assumes equipartition of energy density in the particles and in the magnetic field (e.g., Burbidge, 1956; Miley, 1980). This internal pressure will be a lower limit and is given by (following O’Dea & Owen (1987))
[TABLE]
where,
[TABLE]
is the magnetic field at minimum pressure condition, and
[TABLE]
is the energy of the particles at minimum pressure. is the ratio of proton to electron energy and is assumed to be 1, is a constant from synchrotron theory that depends on the spectral index, and the upper and lower cut-off frequencies, is the radio luminosity within the upper and lower cut-off frequencies. We use = 2.8 by assuming a spectral index of , 10 MHz and 100 GHz as the lower and upper cut-off frequencies, respectively. is the volume filling factor and is assumed to be equal to 1. To estimate the volume V (in cm*-3*), we make an assumption that the radio lobe is cylindrical with half the lobe width as the radius and core-to-hotspot distance as the height of the cylinder. The radio luminosity is calculated using (O’Dea & Owen, 1987)
[TABLE]
The average for quasars is dyne cm*-2*. We follow the same procedure and estimate minimum pressure for the radio galaxies studied here. The average dynes cm*-2* for radio galaxies. This is roughly a factor of 4 lower than the in quasars. Given that a greater fraction of radio galaxies are below redshift of 1, the lower of radio galaxies may likely be a result of the increasing with redshift as seen in Figure 11 and described below.
Figure 11 shows the minimum lobe pressure as a function of redshift. The lobe pressure correlates strongly with redshift with a correlation coefficient of 0.83 obtained at a significance greater than 99.99 per cent. The slope of the linear fit is 1.32. The plot suggests that high redshift sources have lobes with higher minimum non-thermal internal pressure (see also O’Dea et al., 2009). This could be an outcome of the increasing radio luminosity and decreasing lobe widths observed in this sample with redshift. This also results in increasing with redshift. Wan & Daly (1996) suggest that inner regions on radio lobes (close to the core) are in pressure equilibrium with their surroundings and the ambient gas densities are found to have generally no dependence on redshift (O’Dea et al., 2009). This implies that the actual source magnetic field strengths are lower than the minimum energy magnetic field strengths. Alternatively, high redshift sources are likely interacting with the ambient medium and have not yet reached pressure equilibrium with their surrounding. High redshift sources have denser ambient gas than at low redshift. This requires high pressures in the radio lobes to maintain pressure equilibrium with the ambient dense gas.
4 Low frequency emission in quasars and radio galaxies
TGSS is a radio survey carried out using the GMRT at 150 MHz and covers a declination range of to . It has the highest resolution of in the low frequency regime. The TGSS Alternative Data Release (ADR1, Intema et al., 2017) re-processed the TGSS images using a robust automated pipeline SPAM (Intema, 2014) which employs techniques such as direction dependent calibration of ionospheric phase errors and image-based flagging to deal with the widespread RFI at the low frequency. We use the TGSS ADR1 images to study the morphology of this quasar and radio galaxy sample at low frequencies and derive total energy transported by the jets during the lifetime of the AGN. We find that the classical FR type morphology seen in the 1.4 and 5 GHz images is no longer the same at 150 MHz. Several sources (see Figure 26) show signatures of wings and extended emission in directions perpendicular and at an oblique orientation to the primary lobes observed at GHz frequencies, reminiscent of X-shaped galaxies (Leahy & Williams, 1984; Leahy & Parma, 1992). Here we simply refer to all such structures as “distorted”. We note that these wings persist in higher resolution GMRT images at 610 MHz (Vaddi et al. 2018, in preparation), attesting to the reality of these features. These distorted radio morphologies may indicate jet axis re-alignment and/or multiple activity episodes in these AGNs.
4.1 Jet Kinetic Power
The synchrotron emitting plasma in the lobes can be used as an indicator of the amount of energy supplied by the jets to the lobes (Godfrey & Shabala, 2016). The time-averaged total bulk kinetic jet power, defined as the total energy transported by the jets over the lifetime of the source, is given by (Equation 12 of Willott et al. (1999))
[TABLE]
where is the total flux density at 150 MHz in W m*-2* Hz*-1*, is the spectral index estimated using the flux density at 150 MHz and 1.4 GHz, and is the luminosity distance. Equation 7 assumes that synchrotron radiation losses are negligible. The average quasar jet power is 6.40.07 1045 erg s*-1*. The KS test shows that the time-averaged jet kinetic power for quasars and radio galaxies may have similar distribution (). The jet power ranges from erg s*-1*, with 3C 470 radio galaxy having the highest jet kinetic power (refer Table 1, and Table 6). The jet power and redshift are correlated which is expected from a flux limited 3C sample. For the quasar sample, we notice that jet power is related to the mass of the supermassive black hole such that higher jet powers are associated with more massive black holes. However, a partial correlation test between jet power and the black hole mass while ignoring the effect of redshift does not show any correlation (=-0.3, =0.5) between jet power and black hole mass suggesting that the observed trend is likely a manifestation of Malmquist bias.
5 Possible Biases
5.1 Selection Bias
The Kharb et al. (2008) radio galaxies have radio power > 1028 WHz*-1* at 178 MHz. Their angular sizes are > 27" and the sources have a redshift range from 0.4 to 1.65. The quasars were selected to have BH mass determinations, at least one side of the source without a radio jet, and radio powers 31026h*-1* W Hz*-1* sr*-1* which for our chosen H0 gives 51026 WHz*-1* sr*-1*; the sources have redshift between 0.3 and 1.78. The black hole masses for the quasars were obtained from McLure et al. (2006), who studied the FWHM of the broad-line emission and typically found v 2000 km s*-1*. Based on these selection criteria, we can say the following about its effects on our results. The selected radio galaxies are NLRG that are mostly in the plane of the sky and lack the broad-lines for BH mass estimates through virial methods. This BH mass selection criteria for quasars, therefore, does not change the results. Since the quasars in our sample are selected such that spectral ageing can be performed, it automatically selects quasars that have larger lobe extents which in turn selects sources that are intrinsically larger or are at larger angles to the line of sight. This gives the quasars in our sample a larger average LLS compared to previous studies. This, however, does not change the general conclusion that quasars have statistically smaller projected sizes than the radio galaxies, in line with unification.
6 Summary and Conclusions
We have presented the results from our 1.4 and 5 GHz study with the VLA of 11 FR radio-loud quasars. The GHz spectral index images reveal spectral steepening along the lobes, moving away from the hotspots. These maps and results will be used in the spectral ageing analysis, the results of which will be presented in Paper II. We compare these quasars with FR radio galaxies to study the radio-loud unification scheme. We consider 13 narrow-line FR radio galaxy sample that span similar redshift, luminosity, and have matched frequency/resolution data. We see that radio-loud unification largely holds but environmental factors also play an important role. We have shown that:
The radio core prominence () parameter, which is a statistical indicator of beaming and thereby orientation, suggests that the radio galaxies considered in our study lie at orientation angles between and , while the quasars lie at orientation angles between and , broadly consistent with the radio-loud unification scheme. 2. 2.
About 90 per cent of the sources have brighter hotspots on the jetted side and the brighter hotspot is equally likely to be on the shorter or on the longer arm indicating that beaming effects may not be the primary contributor to the observed structural and brightness asymmetry. 3. 3.
The hotspot spectral index () distribution is similar for quasars and radio galaxies. The is strongly anti-correlated with redshift and it is argued to be due to higher synchrotron losses at higher redshifts due to density-dependent effects while luminosity-dependent or size-dependent effects are not seen. 4. 4.
The average projected LLS of radio galaxies is 1.60.3 times larger than those of quasars in our sample. This is comparable to the ratio of the sine of average orientation angles calculated from . This is consistent with the orientation based unification scheme and LLS is likely probing orientation. 5. 5.
The median axial ratios of quasars are roughly a factor of 1.4 lower than the radio galaxies. This may be explained by projection in length since the differences in lobe widths between quasars and radio galaxies are minimal. However, larger and significant lobe distortions seen in quasars compared to radio galaxies support the importance of intrinsic/environmental differences. It is possible that the direction of the quasar jets in our sample is changing with time, as in the \qqdentist drill model, while radio galaxies have stable long-term collimated jets. Alternatively, evolutionary relationship between quasars and radio galaxies, however, cannot be ruled out. 6. 6.
The misalignment angles for quasars are roughly 1.8 times greater than radio galaxies likely indicating an orientation effect. However, the radio core prominence which is also a statistical indicator of orientation does not show a correlation with the misalignment angle. Environmental factors may be instrumental in this lack of correlation. Rc could be influenced by small-scale environment asymmetries while the misalignment angle could be influenced by both small scale and large scale environment asymmetries. 7. 7.
The minimum lobe pressure is strongly correlated with redshift. High redshift sources have lobes with higher minimum internal pressure to maintain pressure equilibrium with the denser ambient gas. 8. 8.
Low frequency emission at 150 MHz with the GMRT reveal the lobes to move away from the classical FR morphologies; many sources reveal extended, winged morphologies that could suggest changes in jet orientation and/or restarted AGN activity in these quasars. 9. 9.
The time-averaged jet power estimated using the 150 MHz data is similar for quasars and radio galaxies.
In summary, several of the observed properties of quasars are consistent with the orientation-based radio-loud unification scheme. However, differences and relations between orientation indicators such as axial ratio, misalignment angle, and larger lobe distortions in quasars can be explained by intrinsic/environmental asymmetries.
Acknowledgements
We thank the referee for constructive and insightful suggestions that have significantly improved the quality of the manuscript. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. The work of Stefi Baum and Chris O’Dea was supported by NSERC (Natural Sciences and Engineering Research Council of Canada). The work of Ruth Daly and Trent Barbusca was supported by Penn State University. 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. We thank the staff of the GMRT who have made these observations possible. GMRT is run by the National Centre for Radio Astrophysics of the Tata Institute of Fundamental Research.
Appendix A Spectral Index and continuum Images of the Quasars
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Antonucci (1984) Antonucci R. R. J., 1984, Ap J , 278, 499 · doi ↗
- 2Antonucci (1993) Antonucci R., 1993, ARA&A , 31, 473 · doi ↗
- 3Antonucci & Barvainis (1990) Antonucci R., Barvainis R., 1990, Ap J , 363, L 17 · doi ↗
- 4Antonucci & Miller (1985) Antonucci R. R. J., Miller J. S., 1985, Ap J , 297, 621 · doi ↗
- 5Antonucci & Ulvestad (1985) Antonucci R. R. J., Ulvestad J. S., 1985, Ap J , 294, 158 · doi ↗
- 6Arshakian (2005) Arshakian T. G., 2005, A&A , 436, 817 · doi ↗
- 7Athreya & Kapahi (1998) Athreya R. M., Kapahi V. K., 1998, Journal of Astrophysics and Astronomy , 19, 63 · doi ↗
- 8Barthel (1989) Barthel P. D., 1989, Ap J , 336, 606 · doi ↗
