Open access peer-reviewed chapter

Evolution of Radio Source Components and the Quasar/Galaxy Unification Scheme

Written By

Costecia Ifeoma Onah, Augustine A. Ubachukwu and Finbarr C. Odo

Submitted: 23 March 2022 Reviewed: 04 July 2022 Published: 06 August 2022

DOI: 10.5772/intechopen.106244

From the Edited Volume

Astronomy and Planetary Science - From Cryovolcanism to Black Holes and Galactic Evolution

Edited by Yann-Henri Chemin

Chapter metrics overview

106 Chapter Downloads

View Full Metrics

Abstract

In this work, a theoretical model is developed for explanation of temporal evolution of extragalactic radio sources via beaming, orientation effects and asymmetries. Equation of the form D≈P±q1+z−m is used to account for the D ∼ P/z relation. Also, D≈D01+z−1+z1+z2 accounted properly for Ω0=1 cosmology than the Ω0=0 counterpart in linear size versus redshift of radio sources. Similarly, D=Dc1∓lnPPc1/2 model explained redshift-luminosity relationship of extragalactic radio sources. The results from the regression analyses are q = +0.003 (r = 0.04) for sources with z < 1 and q = −1.59 (r = −0.6) for all z≥1 sources. A critical linear size, Dc of 316kpc which matches the maximum theoretical linear size, Dmax of 0.15D0 at a critical redshift zc∼1 and a critical luminosity Pc=26.33WHz−1 are obtained. The indication of all these results is that the linear size of radio sources evolves up to a certain limit in D–P plane and thereafter decreases with increasing luminosity as predicted in this work.

Keywords

  • AGNs
  • ESS quasars
  • ESS galaxies
  • evolutions: temporal
  • cosmological and general

1. Introduction

1.1 The active galactic nucleus (AGN)

The active galactic nucleus (AGN) is the existence of energetic phenomena in the nuclei, or core regions, of galaxies that cannot be clearly and directly explained to the interactions between the stars and interstellar medium. The source of radiation of AGN is definitely believed to emanate from the gravitational potential of gas from a supermassive black hole accreting mass at the core of the host galaxy. This emitted radiation by AGN is throughout the electromagnetic spectrum. The energy radiation in AGN is non-thermal unlike spectra of stars rather it is primarily as a result of the process of synchrotron radiation. In this scenario, power-law spectra, Sννα and high degree of linear polarization assigned to AGN object stand as evidence of synchrotron theory. AGNs have typical luminosity in the range of 1033 to 1040WHz1 [1].

Generally, an AGN possesses peculiar properties like, intense bright and point-like nucleus, radio cores with compact flat spectrum, highly ionized gas relativistically beaming out, variable fluxes observed on a wide timescale range from minutes upwards, extremely high luminosity from the range of 1061014 solar luminosity, narrow, broad, and sometimes without lines of emission, extended radio jets and lobes.

1.2 Classification of AGN

The naming of AGN into classes and subclasses is mostly based on the exhibition of their properties or morphologies. Some AGN classes may be similar as a result of their evolution; some may be due to observed variability of luminosities [2]. Some AGN are classified based on their viewing angle as seen by the observer, which depends mostly on the obscuration of torus [3]. This brings about an unification scheme in AGNs due to the relativistic beaming and orientation effects.

Optically, AGN may be classified as type I or type II based on the optical spectral line emission. The classification of AGN using their response to radio-loudness is preferable because of the clearer observation obtained from the radio band compared to any bands’ counterpart in the electromagnetic window. In this scenario, AGN can be classified as radio-loud and radio-quiet depending on their radio brightness. When they have ratios of radio (5GHz) to optical (B-band) flux F5/FB10, they are called radio-loud, otherwise known as radio-quiet if F5/FB10 [2, 4]. Since this work is entirely based on the observed radio properties of various AGN samples plus the clarification and consideration of radio window observation, the classification of AGNs based on radio-loudness is represented in Figure 1.

Figure 1.

Classification of AGNs based on radio-loudness.

1.3 Radio-loud AGNs

Radio-loud sources emit more energies in the radio band than in the optical waveband, and hence possess F5/FB fluxes ≥ 10 [2]. There are about 15–20% of radio-loud AGN. They are mostly elliptical galaxies in accordance with Hubble turning-fork proposition (population II)—meaning that they contain mostly of old stars with little interstellar gas. These objects are further classified into two: high luminosity objects and low luminosity objects.

1.4 High luminosity objects

They have total radio luminosity, P1781035WHz1, with a highly ionized accretion disk [5]. These comprised of radio-loud quasars, Compact Steep-Spectrum Sources (CSSs) and Fanaroff-Riley class II radio galaxies. The works are mostly on these classes of AGN.

1.4.1 Fanaroff-Riley class II radio galaxies

These are known for their high luminous intensity with the possession of extended double lobes/jets where one side is Doppler enhanced. They also have smooth jets as a result of highly supersonic flows. The jets are also edge-brightened and terminate in hotspots [5, 6, 7]. They polarize linearly with the electric field vector being perpendicularly to the jets. In a unification scheme, [8] suggested that FR II sources are misaligned counterparts of core-dominated quasars.

1.4.2 Quasars

Quasars are classified further into core-dominated and lobe-dominated quasars. These are core-dominated if the radio emissions emanate mostly from the core, otherwise it is a lobe-dominated one.

The core-dominated ones possess properties like flat radio spectra with spectra index, α0.5,sννα as a result of synchrotron self-absorption mechanism, cores with extremely brightness, broad emission lines and one-sided jets/lobes. These types of quasars dominated the survey at high frequencies and high redshifts. These classes of quasars show more asymmetry than the lobe-dominated counterpart [9].

On the other hand, lobe-dominated quasars, unlike the core-dominated, have two extended lobes straddling a weak compact core. They are also high luminosity sources with total luminosity, P1781035WHz1. They are characterized by steep radio spectra (α>0.5). They show spectra with broad emission lines; hence, they are referred to as broad-line region sources (BLRS). They also have higher redshifts when compared with radio galaxies [9].

1.4.3 Compact Steep-Spectrum Sources (CSSs)

The CSS sources are characterized by sharp peaks exhibited in their radio spectra. As their names imply, they are compact bright radio sources with a population up to 30% [10]. They are also called “youth” scenario [10] being believed to be the younger phases of powerful large-scale extragalactic radio sources. They have a small radio size, D15kpc, with a steep radio spectrum, α0.5, a very high radio luminosity, logP>1026WHz1 at frequency, ν = 2.7 GHz [10]. CSS radio sources exhibit brightness temperature up to 1010K [10]. Their radio structure is symmetric with low radio polarity and large Faraday rotation measures.

They are CSS radio galaxies if they have double lobes with weak jets and cores emitting weakly, otherwise CSS quasars if they exhibited brighter cores and jets [11]. Majority of CSS radio jets are one-sided and superluminal [11]. From the theory of orientation-based unification scheme, the morphological difference between the CSS radio galaxies and CSS quasars is that objects seen close to the line of sight of the observer are referred to CSS quasars otherwise radio galaxies [2].

1.5 Low luminosity radio sources

Low luminosity radio sources have total radio luminosity, P178<1026WHz1 with less ionized accretion disks [12]. Examples of these sources include FR I radio galaxies and BL Lacertae objects.

1.5.1 Fanaroff and Riley class I (FR I) radio galaxies

According to [12], FR I radio galaxies are characterized by extended double-lobed with low frequency, ν ∼ 400 MHz. FR I has RFR<0.5 as the source faints away from the nucleus, while FR II with brightness further away from the nucleus has RFR0.5. RFR ratio is the ratio of the distance between the regions of highest surface brightness to the lowest brightness contour of the central galaxy.

Moreover, FR I sources are symmetric with smooth and continuous jets which begin as one-sided nearer the core and two-sided at a few kilo parsecs away. FR I sources are located in rich clusters that highly emit x-ray gas. The x-ray gas sweeps back and distorts the FR I radio structure as it moves across the interstellar cluster, hence giving FR I object narrow-angle-tail or wide-angle-tail according to the strength of the ram pressure of the gas [13, 14].

1.5.2 BL Lacertae sources (BL Lacs)

BL Lacs objects are the most violent AGN known. They have properties like very weak or sometimes no radio emission lines, compact radio core, rapid and high peak variable fluxes, superluminal flows. They are also known as blazars just like optically violently variable (OVV) quasars.

1.6 Radio quiet objects

These are objects that emit more in the optical window than they do in the radio band. AGN sources are said to be radio quiet if they emit more of their energy in the optical waveband than in the radio waveband. They have properties like F5/FB fluxes ≤10, low luminosity at 6 cm less than 1026WHz1, short jets, few relativistic particles and total weakness of the radio sources [15]. Examples of these objects are radio quiet quasars and seyferts (seyfert I and seyfert II). In this work, it is centered more on high luminosity radio-loud objects like FR II radio galaxy and quasar.

1.7 Features observed in AGNs

The observational morphological features of a radio source are radio core, jets, lobes and hotspots, though every source may not exhibit all these features.

1.7.1 Radio core

This is the core engine where the energetic radio emission mechanism in EGRSs is assumed to originate. This core is divided into steep spectrum cores and ultra-compact flat spectrum [16]. The steep spectrum cores are characterized in some radio galaxies by having more extended radio cores of few kilo parsecs in size and steeper spectra (α0.5), example as found in seyferts and some spiral galaxies. On the other hand, the ultra-compact cores’ counterparts possess properties like 1kpc in sizes and flat radio spectra (α<0.5), signifying synchrotron self-absorption that arises as a result of the re-absorption of some radiate relativistic electrons within the radio source. Quasars objects’ core emissions appear to be more powerful than that of the radio galaxies.

1.7.2 Jets

These are the conduits through which the high-energy particles are transporting from the cores to the other extended radio structures. A typical radio jet is expected to be at least four times as long as its width, in line with the radio core and separable from other features at high resolution [17]. The radio jets which can be one or two-sided exhibit properties like emissions with steep spectra of spectral index, α ̴ 0.5–0.9.

1.7.3 Radio lobes

These are one of the extended structures of radio sources that cover a range of hundred kilo parsec to a few mega parsecs. These are characterized by possession of non-thermal steep spectra of spectral indices, α0.5, with a high degree of polarization at high frequency. They exhibit morphological features like tail, plumes, bridges and haloes. Tails are structures formed as a result of deflected plasma interacting with external medium, while plumes are extended regions with low luminosity that faints away from the whole source. Bridges occur in the inner lobe regions of radio galaxies, while the haloes are low surface brightness structures containing old aged plasma [18].

1.7.4 Hotspots

Hotspots are the brightest region of the extended radio structure formed at the end of the lobe where kinetic power of the jets is converted into random motion within the relativistic plasma and strengthened magnetic fields [19]. They have a linear size of 1kpc and steep spectra, α ̴ 0.5–1 slightly flatter than that of the surrounding diffuse emission [20]. See Figure 2.

Figure 2.

Unified structure of an typical extragalactic radio source [21].

It has been established that the appearance of EGRSs is substantially modified by relativistic beaming and orientation of the radio axes with the line of sight, leading to asymmetries in the observed radio structures. Similarly, radio sources are known to undergo some form of cosmological as well as temporal evolution. However, the amount of relativistic beaming and the nature of the evolution present in different classes and subclasses of the EGRSs are still a subject of intensive research. In particular, different source samples show a wide range of the amount and nature of temporal evolution as reported in literature. Hence, the aim of this work is to analytically examine the observed radio properties of different samples of EGRSs for radio source structural asymmetry, use relativistic beaming and source orientation model to explain any observed structural asymmetry in the radio sources and finally develop a model that can unambiguously explain the temporal evolution in extragalactic radio sources.

Advertisement

2. Theory of relationships

In this approach, relations between various parameters of the radio sources would be derived based on established laws and theories under certain assumptions. The theories of relationships are outlined.

2.1 Theory of temporal evolution in radio sources

The standard relationship for an ideal temporal evolution model of extragalactic radio sources can be studied using the relationship between the two key parameters known as the observed linear size (D) and the spectral luminosity (Pv) in a general power-law function as [22]:

DD0Pv±qE1

where D0 (a constant) represents D at z ∼ 0 while the slope (q) is the temporal evolution parameter.

Meanwhile, from the Friedmann-Robertson-Walker universe, the radio angular size–redshift θz relation is described [23, 24] as

θ=D1+z2dL,E2

where D can be represented using the unit of kpc and dL is the luminosity distance expressed [24, 25] as

dL=2cH0Ω02Ω0z+Ω02Ω0z+1121,E3

and c is the speed of light. Moreover, the radio luminosity (P) of a radio source at redshift (z) can be defined as a function of the spectral flux density (Sv) at observing frequency (ν) according to the [26] as:

P=P01+zβE4

where α stands for spectral index. Thus, from [27] work, the relationship between P and z can be approximated to a simple power-law function of the form:

P=P01+zβE5

where β is the slope of the P–z which is supposed to be constant over all values of z for a given sample of sources. In this scenario, a significant correlation for all values of z in the line with β>0 is expected in a P–z data of different samples of extragalactic radio sources, due to luminosity selection effects in flux density limited samples [27].

Besides, researchers [23, 24, 28] suggest that linear sizes of extragalactic radio sources evolve with cosmological epoch in the form:

DD01+zkE6

where D0 is the normalized linear size which depends on the assumed cosmology and k represents the evolution parameter which could be a function of both cosmological evolution and luminosity selection effects.

Alternatively, with different values of Ω0, the observed θz data of extragalactic radio sources deviate significantly from the standard Friedmann models [29], as a result of an entanglement of two effects, namely linear size evolution [29] and luminosity selection effect [30]. Hence, the linear size evolution of extragalactic radio sources can therefore be expressed as a function of both redshift and luminosity in the form of [24, 30, 31]

DPzP±q1+zmE7

where m is the residual cosmological evolution parameter defined [24, 27] as:

m=k±E8

When the effect of luminosity is controlled and this depends on the product of temporal evolution (q) and luminosity selection effect (β).

It is certain from above relation that if luminosity selection effect is above cosmological evolution, then, when the effect of luminosity is corrected, m < 0, so that Eq. (8) becomes

k<0E9

and q>kβ. For such sources, a significant positive D–P correlation is envisaged, suggesting a linear relationship between the two parameters. On the other hand, if there is a residual linear size evolution at certain values of k, q < 0 and m > 0, then q<kβ, indicating a significant D–P anti-correlation and also implies that the luminosity decreases with the expansion of the sources. Hence, the values of k are expected to vary in a given sample. In this scenario, the temporal evolution model for an assumed cosmology can be constrained using the D–z plane.

It is also believed that D should not increase with z for all z. Although, considering the fact that radio sources are not just rigid rods and the θz plane depends on the assumed cosmology, which is also, characterized by the value of Ω0, in principle, this may not always be the case. Hence, the P–D track of all inclusive radio source samples could be defined [24, 29] as:

P=Pmaxexp±DDc12E10

where Pmax is the maximum luminosity and Dc is the critical linear size at which the Pmax is emitted by a radio source. Rearranging Eq. (10) gives

D=Dc1lnPPc1/2E11

Eq. (11) comprises of two separable components corresponding to k<0 and k>0+. The D–P relation (c.f. Eq. (11)) is shown in Figure 7.

However, using Eq. (3) in (2), the linear size can be defined as:

D=D0Ω0z+Ω02Ω0z+11/211+z2,E12

where D0=2H0Ω02 is the intrinsic linear size and a constant. Hence, for inflationary universe, Ω0=1, [32], the linear size of a radio source depends on z in the form [27]:

DD01+z1+z1+z2E13

Eq. (13) indicates two components of z: z < 1 and z > 1, where D increases with increasing z in the first so that k > 0, while the revised is the case in the last component for all D and k. On the other hand, assuming Ω0=0, on supposition of a low-density universe, for which dL is given as [33]:

dL=2cH01+z21,E14

Eq. (2) yields,

θ=DH01+z22c1+z21.E15

Hence, D can be expressed in relation with z as

D=D0111+z2E16

where the intrinsic radio size, D0=θ2cH0. The variation of D with z for both cosmological models is shown in Figure 3.

Figure 3.

Variation of D with z for Ω0=1 (a) and Ω0=0 (b) cosmologies.

It is obvious from Figure 3 that there is increase in D as z increases up to a certain maximum point known as a critical value of DC0.146D0 (kpc) at a redshift maximum called critical zc1, and decreases after for Ω0=1. This implies that for parameters, z < 1, k > 0 and q > 0, there is a positive temporal evolution, while for z >1, k < 0 and q < 0 a negative temporal evolution is envisaged which is in good agreement with the predictions of Eq. (13).

On the other hand, there is increase in D up to a critical point Dc0.8D0 after which it remains constant for Ω0=0. It is now clear that in the two cosmological models, there is indication of Dc at which luminosity is maximum, suggesting that for any assumed cosmological model, the Dc value obtained from D–P turnover can still be found using the D–z plane of the same data. Hence, in any assumed cosmological model, there will be an expected range of the Dc value of 0.15D0 to 0.8D0 bounded by the two cosmologies. In this scenario, temporal evolution in the current sample of EGRSs would be modeled in terms of the current inflationary one with Ω0=1.

2.1.1 Relativistic beaming based on orientation and radio source asymmetries

The standard relation for explaining an ideal relativistic beaming and orientation effects for extragalactic sources is often carried using a key parameter known as core-dominance, (R) defined [34] as:

R=PCPE=P5GHzCP1.4GHzE1.45αE1+zαEE17

where PC represents core-luminosity at 5 GHz, PE is the extended/lobe luminosity at 1.4 GHz and αE is the lobe spectral index. However, if relativistic beaming effect at small orientation angle, then R can be expressed in terms of the jet speed (β) and inclination angle (ϕ) in the form [35]:

R=PCPE=RT21βcosϕn+α+1+βcosϕn+αE18

where RT is the value of R at ϕ=90° and n is a parameter that depends on the assumed flow model of the radio jet. For radiating plasma with continuous jet n = 2, otherwise n = 3 if the jet consists of blobs.

An obvious outcome of relativistic beaming and orientation effects in AGNs is the wide range of asymmetry observed in their radio structures. The radio source asymmetry can be explained using the arm-length ratio (Q), defined as the ratio of the distance, from the central engine, of a plasma element emitting radio waves on the approaching jet side (dapp) to that on the receding jet side (drec), [36, 37] given as:

Q=dappdrec=1+βcosϕ1βcosϕE19

On the other hand, [38, 39] suggested that,

x=Q1Q+1E20

where, x represents the index of the asymmetry.

This x parameter which believed to have better relationship with orientation when compared to Q can further be defined in terms of the viewing angle as [40, 41]:

x=βcosϕE21

The relativistic beaming in AGN at small angles to the line-of-sight is fundamentally characterized by beaming enhancement factor (δ) expressed [40, 42] as:

δ=γ11βcosϕ1=E22

where, γ is the bulk Lorentz factor of the jet [8, 42, 43] defined as:

γ=11β212E23

Also, assuming ϕ=00 in Eq. (18) and analyzing further, [44, 45] gives

RmaxRTγ22γ212RTγ4E24

While at α=0 and β1, for high luminosity radio-loud AGN sources emitting at small angle to the line-of-sight of the observer, (to a first approximation) Eq. (18) reduces to

cosϕm=12RmRT1nE25

where Rm is the mean value of the R-distribution and ϕm is the mean observation angle.

Thus, it can be shown from Eq. (23) that the asymmetry parameter x can be expressed in terms of the beaming enhancement factor as:

x=δγ1δγE26

Eq. (26) implies that there is an association of relativistic beaming and radio source asymmetry. In asymmetric sources with Q > 1, x–D anti-correlation is expected if geometric projection at small viewing angles is responsible for the observed asymmetry. Following [40, 41], we assume a linear x–D relation of the form:

x=xmaxλD,E27

where xmax represents the maximum x for a sample at D0ϕc0 to the line-of-sight [44] and λ is the slope. Thus, if β1 for the relativistic jets, analysis of Eq. (20) down to Eq. (28) for optimum beaming gives [40, 45]:

ϕcsin11γcos1xm,E28

Hence, if relativistic beaming at small viewing angles is responsible for the observed structural asymmetry, the critical viewing angle ϕc as well as the Lorentz factor γ can be obtained using the x–D data.

2.2 Statistical analyses and results

2.2.1 The source samples

The data used in the present analysis were drawn from a well-defined source sample of [46] compilation which contains required information on the two objects of interest–ESS quasars and ESS radio galaxies, [34] compilation of 542 extragalactic radio sources and the deep VLA sample of FSRQs compiled by [47]. In these samples, there is wide dispersion in the distributions of the observed parameters.

2.2.2 Distributions of observed radio parameters

The distributions of the linear size, D, of the extended steep-spectrum sources (ESSs) in logarithm scales are represented in Figure 4a. The graph shows D-values of 1896.3, 17.1 and 1879 kpc for the maximum, minimum and range respectively for ESS quasars, while D-values of 5853.30, 29.80 and 5823.5 kpc were obtained for maximum, minimum and range respectively for ESS galaxies. The entire sample yields D-values of 5853.3, 17.1 and 5836.2 kpc for maximum, minimum and range respectively. Further analyses yield median D-values of 148.90, 323.59 and 201.90 kpc for ESS quasars, ESS galaxies and entire sample respectively. The mean D-values of 144.54 ± 7.52 kpc and 288.40 ± 33.66 kpc were obtained respectively for ESS quasars and radio galaxies.

Figure 4.

Distribution of (a) D, (b) z and (c) P respectively for ESS quasars (lines) and ESS galaxies (plane).

The distributions of the redshift, z of the ESSs in logarithm scales are represented in Figure 4b. The graph shows z-values of 2.88, 0.05 and 2.83 for the maximum, minimum and range respectively for ESS quasars, while z-values of 3.22, 0.006 and 3.21 were obtained for maximum, minimum and range respectively for ESS galaxies. The entire sample yields z-values of 3.22, 0.006 and 3.21 for maximum, minimum and range respectively. Further analyses yield median z-values of 1.89, 1.17 and 1.62 for ESS quasars, ESS galaxies and entire sample respectively. The mean z-values of 1.95 ± 0.01 and 1.30 ± 0.02 were obtained respectively for ESS quasars and radio galaxies.

The distributions of the radio luminosity, P of the ESSs in logarithm scales are represented in Figure 4c. The graph shows P-values of 27.90, 25.71 and 2.18 WHz1 for the maximum, minimum and range respectively for ESS quasars, while P-values of 28.01, 24.97 and 3.04 WHz1 were obtained for maximum, minimum and range respectively for ESS galaxies. The entire sample yields P-values of 28.01, 24.97 and 3.04 WHz1 for maximum, minimum and range respectively. Further analyses yield median P-values of 27.04, 26.39 and 26.84 WHz1 respectively for ESS quasars, ESS galaxies and the entire data. The mean P-values of 27.01 ± 0.01 WHz1 and 26.39 ± 0.02 WHz1 were obtained respectively for ESS quasars and radio galaxies Figure 4.

2.2.3 D: P/z correlation

Figure 5 represents the scatter plot of linear size, D against the redshift, z. The median value data in seven redshift bins is superimposed on the plot. Critical investigation of the plot shows that on average, the linear size increases with increasing redshift up to a value logDc = 2.5 kpc (Dc = 316.23 kpc) at zc=1, after which it decreases with increasing redshift [24]. This is in agreement with the prediction made in Figure 3(a). Hence, the present data obviously proved consistency with the inflationary model of the universe Ω0=1. The median values give significant trends with correlation coefficients of +0.95 and − 0.90 for zc=1 and zc1 respectively. Results of the regression analyses of the D–P/z data for zc=1 and zc1 are summarized in Table 1.

Figure 5.

Scatterplot of logD (kpc) against z for entire sources (circle) with median values (square) superimposed [24].

ParametersD–PD–z
D0qrD0kr
z<12.370.0030.0042.55−1.09−0.12
z12.67−1.59−0.6−1.592.67−0.5
P<Pc4.910.170.212.510.81−0.05
PPc16.48−0.32−0.72.63−1.50−0.6

Table 1.

D–P/z regression analyses results for both z < 1 and z1/P<Pc and PPc [24].

In modelling the temporal evolution of the sample, Figure 6 represents the scatter plot of projected linear size (D) and the radio luminosity (P). Similarly, the median value data in eight uniform luminosity bins are superimposed on the plot. There is an obvious trend indicating that the linear size first increases

Figure 6.

Scatterplot of logD(kpc) against logP (WHz1) for both sources with PPc (circle), PPc (square) and median (triangle) values superimposed [24].

with increasing luminosity up to a certain value and thereafter decreases with increasing luminosity. The median value data showed very significant trends. This suggests that the turnover occurs at a critical point of luminosity, logPc=26.33WHz1 and logDc = 2.51 kpc (316.23 kpc) [24]. A summary of the results of these regression analyses of the D–P and D–z data for sources with P ≤ Pc and P > Pc is presented in Table 1.

The results in Table 1 did not show any obvious trend in D–P relation for sources with P<Pc. However, for sources with PPc, there is a fairly strong correlation. Therefore, the weak trend found in the region below z = 1 and P<Pc is believed to be due to the effects of luminosity selection. In this scenario, the low redshift samples, z < 1 have more impacts on average luminosity-redshift plane than the high redshift, z1 counterparts in any flux density limited samples.

2.2.4 Luminosity selection effect on temporal evolution model

Figure 7 represents the P–z scatter plot of the sample. There is a turnover at critical point logPc=26.33 W/Hz and z ∼ 0.05 in the P–z plane. This point of P is presumably the value of luminosity in Figure 6 that will correspond to Dc at zc1. However, the critical luminosity, Pc found in the P–z plane is inconsistent with sources at z = 1. Hence, the luminosity-redshift relation of the current data did not assume Ω0=1 model but Ω0=0 cosmological model (low-density universe) [24].

Figure 7.

Scatterplot of logP (W/Hz) against log (1 + z) for all sample.

In Figure 6, the effect of Eqs. (10) and (11) in the light of the temporal evolution model was considered. Hence, the data are grouped into two; P < Pc and PPc. The results of the regression analyses are shown in Table 1. These opine that the linear size, D of radio sources increases up to critical radio luminosity, logPc=26.33WHz1 and decreases thereafter. This is in agreement with Eq. (10), hence suggesting that eqn. (10) is approximately correct to zero order. This is applicable when critical linear size Dc is used. Dc=316kpc is obtained from the D - z plane at the turning point. This Dc corresponds to theoretical Dc=0.14D0 at zc=1 for Ω0=1 in Figure 3a and logPc = 26.33 WHz1 in Figure 7. The indication of this is that D0 (∼ 2100 kpc in the observed data) is approximately the linear size at the earliest epoch (at z ∼ 0.02 in the observed data of the present sample). The correlation coefficient, r ∼ +0.4, −0.5 and − 0.9 respectively for sources with P < Pc, P ≥ Pc and the median values respectively were obtained. These suggest that there are positive and negative correlations in the D–P track at P < Pc and PPc respectively, indicating that the Dc=316kpc of the observed samples is consistent with Dc = 0.14D0 just in accordance with the theory for the inflationary model of the universe, Ω0=1 only, at zc=1 and Pc=26.33WHz1, proving Eq. (11) to be perfectly correct. Hence, temporal evolution in extragalactic radio sources can be explained.

2.2.5 The x – D relationship

There is a fairly significant trend in the x–D plot, which is obvious at the upper envelope function. This yields: x = 0.35–0.0006 D with a correlation coefficient r ∼ − 0.5, chance probability ρ ∼ 10−10 and critical viewing angle, ϕc70o, which corresponds to γ=1.1. The upper envelope function gave correlation coefficient, r ∼ 0.9, ϕc48o and γ=1.3. The analyses for separate objects, radio galaxy and quasar sub-samples, of upper envelope functions yield ϕc59o; γ = 1.2 and ϕc33o; γ=1.8, for radio galaxies and quasars, respectively. These results correspond to angular separation ϕsep of 26o, based on the upper envelope functions [40]. The results are shown in Table 2 and Figure 8.

Source parametersRadio galaxiesQuasarsGalaxies + quasars
xm0.160.250.20
λ+0.000056−0.0002+0.000002
r+0.3−0.3−0.5
Sig. (%)5.05.05.0

Table 2.

x–D regression analyses results [40].

Figure 8.

Scatterplot of the fractional separation difference (x) as against projected linear size (D) for quasars (•) and radio galaxies (∆) [40].

The results opine that the relativistic beaming and source orientation effects are the major cause of large-scale structural asymmetries observed among powerful extragalactic radio sources, which are more obvious in core-dominated sources with large core-to-lobe luminosity ratio.

2.3 General discussions and conclusions

2.3.1 Discussions

The study of the redshift effect dependence of radio size in extragalactic radio sources has been of great interest since inception of the universe. This is obviously important in cosmological studies more especially in testing of the evolution in the extragalactic radio sources. In this scenario, a theoretical model that can best interpret the D–P track is developed. According to [48] the analyzes obtained, the high redshift radio sources show high bending angles, distorted and smaller structures than the low redshift radio sources. According to the researchers, the effects on the radio source evolution depend mostly on intracluster/ or interstellar medium. A plausible model was first developed by [49] for strong linear size evolution for radio galaxies. This model explained that the typical radio sources, both giant galaxies and sub-galactic quasars, evolve at high redshifts, z.

The temporal evolution of radio linear size could be interpreted in the light of observed D–z plane for radio source samples. Hence, the D–z correlation helps in constraining cosmological/temporal models. According to [50], the amount of linear size evolution required to interpret the observed θz data can be in the range of k = 2.0 and 1.5 to k = 1.2 and 0.75 for Ω=1 and 0 respectively. A clear investigation from the median values in Figure 5 shows that the D–z curve of the observed data is in good agreement with that of the theoretical D–z plane of Figure 3. In other word, for Ω=1, there is an increase in the linear size of radio sources up to a certain point known as critical linear size, Dc which corresponds to critical luminosity, Pc at critical redshift, zc1 and thereafter decreases. This steep change in P–z relation of extragalactic radio sources may occur around z = 0.3 [27, 33] or z = 1 (see also [46], Fig. 2) [51]. This is due to the luminosity selection effect occurred as a result of Malmquist bias in flux density limited samples [27]. In this effect, the present work adopted samples based on zc=1 as predicted earlier with the theory. This raise and fall of the radio sources in the D–z plane are consistent with the observational results showing that radio galaxies and quasars exhibit different D–P correlations.

In this scenario, the median values of D–P of the entire data superimposed in Figure 6 clearly show that the observed D–P data of the sample are consistent with the theoretical curves (c.f. Figure 3a). Hence, this suggests that the P–D data of the present sample could be used to constrain the temporal evolution model.

A cursory look in Table 1 shows that there is steep change in q as predicted above in the subsample around z ≥ 1 with fairly trends with a significant correlation in the median values obtained from the eight appropriate luminosity bins of the entire sample. The results show that the apparent lack of D–P correlation in the samples below z = 1 was due to the luminosity selection effects in the observed data. In other words, in flux density limited samples, the low redshift (z < 1) sources exhibit more dependence on luminosity as a function of redshift than the high redshift (z > 1) source counterparts. Below z = 1, the P–z weak correlation is expected, while Beyond z = 1, the luminosities have a much milder dependence on z. The regression analyses yield q = +0.002 for all samples around z ≤ 1 and q = −1.59 at z = ≥ 1. Hence, at z < 1, there is a positive temporal evolution model while at z ≥ 1, negative temporal evolution is obtained for zc1 and Ω0=1 as predicted. This shows that the temporal evolution model can be well understood using Ω0=1 than Ω0=0. In this scenario, the temporal evolution is positive when z < 1, k > 0 and q > 0, and otherwise for z > 1, k < 0 and q < 0. Thus, the linear size increases as a function of radio luminosity up to a maximum value called critical D-value, Dc=316kpc in consistence with the maximum theoretical linear size, Dmax=0.14D0 at zc1 and Pc=26.33WHz1 and thereafter decreases with increasing luminosity, as predicted earlier in this work.

Several authors have argued that the unified scheme will vanish out if the derived positive D–P dependence for radio galaxies (DPq for q ∼ 0.3) contrasts with the negative D–P dependence in quasars [52]. Ubachukwu and Onuora [50] suggested an inverse correlation of D–P of the form DP0.52. They found out that radio galaxies locate at low redshifts with q ∼ 0.3 [31], while radio quasars located at high redshifts with q ∼ −0.5 [48]. The D–P turnover in extragalactic radio sources occurs around Dc = 100kpc [29]. It was argued [24] that at D ∼ 1Mpc and critical luminosity, Pc26WHz1, extragalactic radio sources evolve and thereafter decrease.

In x–D regression analyses, there is an apparent lack of x–D anti-correlation in the entire samples and radio galaxies subsample. There was a fairly x–D anti-correlation in quasars subsample which was obvious in the upper envelope function with r ∼ −0.5, 0.30 and − 0.3 for the entire sample, radio galaxies and quasars respectively. Hence, the results suggest that orientation effects may be more necessary in explaining the properties of the radio quasars population than that of the radio galaxies populations in which the expected x–D anti-correlation was not observed and also implies that the two objects are the same, the difference is just the angle at which the observer took in viewing them.

Advertisement

3. Conclusions

The following major conclusions were drawn from the work. A satisfactory theoretical model was obtained to explain the observed D–P track, with sources at z ≥ 1 showing stronger anti-correlation than those around z < 1. A problem was developed for these objects to match-up with D–z positive correlation as predicted in this paper at z < 1. Finally, the distributions of radio galaxies and quasars on the x–D plane are consistent with the beaming scenario, with quasars being more consistent with relativistic beaming and source orientation model than radio galaxies.

References

  1. 1. Wang TL. Black hole mass and velocity dispersion of narrow line region in active galactic nuclei and narrow line Seyfert 1 galaxies. Astronomy & Astrophysics. 2001;377:52. DOI: 10.1051/0004-6361:20011071
  2. 2. Urry CM, Padovani P. Unified schemes for radio-loud active galactic nuclei. Publications of the Astronomical Society of the Pacific. 1995;107:803
  3. 3. Antonucci RRJ. Unified models for active galactic nuclei and quasars. Annual Review of Astronomy and Astrophysics. 1993;31:473. DOI: 10.1146/annurev.aa.31.090193.002353
  4. 4. Kellermann KI, Sramek R, Schmidt m, Shaffer DB, and Green R. VLA observations of objects in the palomar bright quasar survey. Astronomy Journal. 1989;98:1195
  5. 5. Kinwah W, Mark C, Gavin R, Kazuhiro S. An electrically powered binary star? Monthly Notices of the Royal Astronomical Society. 2002;331:221
  6. 6. Garrington ST, Conway RG. The interpretation of asymmetric depolarization in extragalactic radio sources. Monthly Notices of the Royal Astronomical Society. 1991;250:198
  7. 7. Zirbel E, Baum SA. On the FR I/FR II Dichotomy in powerful radio sources: Analysis of their emission-line and radio luminosities. Astrophysical Journal. 1995;448:521
  8. 8. Ubachukwu AA, Chukwude AE. On the relativistic beaming and orientation effects in core-dominated quasars. Journal of Astrophysics & Astronomy. 2002;23:253. DOI: 10.1007/BF02702285
  9. 9. Arshakian TG, Longair MS. An asymmetric relativistic model for classical double radio sources. Monthly Notices of the Royal Astronomical Society. 2000;311:846
  10. 10. O’Dea CP. The compact steep-spectrum and gigahertz peaked-spectrum radio sources. Publications of the Astronomical Society of the Pacific. 1998;110:493
  11. 11. Akujor CE, Spencer RE, Zhang FJ, Davis RJ, Browne WA, Fanti C. MERLIN observations of steep-spectrum radio sources at 6 cm. Monthly Notices of the Royal Astronomical Society. 1991;250:215
  12. 12. Fanaroff BL, Riley JM. The morphology of extragalactic radio sources of high and low luminosity. Monthly Notices of the Royal Astronomical Society. 1974;167:31
  13. 13. Owen FN, Laing RA. CCD surface photometry of radio galaxies -I. FR class I and II sources. Monthly Notices of the Royal Astronomical Society. 1989;238:357
  14. 14. Prestage RM, Peacock A. The cluster environments of powerful radio galaxies. Monthly Notices of the Royal Astronomical Society. 1988;230:131
  15. 15. Ferrari A. Modeling extragalactic jets. Annual Review of Astronomy and Astrophysics. 1998;36:539
  16. 16. Miley GK. The structure of extended extragalactic radio sources. Annual Review of Astronomy and Astrophysics. 1980;18:165
  17. 17. Perley RA, Bridle AH, Willis AG. High-resolution VLA observations of the radio jet in NGC 6251. Astrophysical Journal (Supplement). 1984;54:291
  18. 18. Kharb P, O’Dea CP, Baum SA, Daly RA, Mory MP, Donahue M, et al. A study of 13 powerful classical double radio galaxies. Astrophysical Journal Supplements. 2008;174:74
  19. 19. Bridle AH, Hough DH, Lonsdale CJ, Burns JO, Laing RA. Deep VLA imaging of twelve extended 3CR Quasars. Astronomical Journal. 1994;108:776
  20. 20. Longsdale CJ. In: Meisenheimer R, Roser HJ, editors. Hotspots in Radio Galaxies. Berlin: Springer; 1989. p. 27
  21. 21. Burke BF, Graham-Smith F. An Introduction to Radio Astronomy. United Kingdom: Cambridge University Press; 1997
  22. 22. Aird J, Wandra K, Laird ES, Georgakakis A, Ashby M. The evolution of the hard X-ray luminosity function of AGN. Monthly Notices of the Royal Astronomical Society. 2010;401:2531
  23. 23. Miley GK. Variation of the angular sizes of quasars with red-shift. Nature. 1968;218:933
  24. 24. Onah CI, Ubachukwu AA, Odo FC, Onuchukwu CC. Evolution of extragalactic radio sources and quasar/galaxy unification. Revista Mexicana de Astronomia y Astrofisica. 2018;54:271
  25. 25. Mattig W. Astronomische Nachrichten. 1959;285:1
  26. 26. Onuora LI, Okoye SE. The variation of radio luminosity with epoch and its effect on the angular diameter-redshift relation. Astrophysical Journal. 1983;270:360
  27. 27. Ubachukwu AA, Ogwo JN. Luminosity selection effects and linear size evolution in the quasar/galaxy unification scheme. American Journal of Physics. 1998;51:143
  28. 28. Gall C, Anderson AC, Hjorth J. Genesis and evolution of dust in galaxies in the early Universe I. Modelling dust evolution in starburst galaxies. Astronomy & Astrophysics. 2011;528:A13. DOI: 10.1051/0004-6361/201015286
  29. 29. Okoye SE, Onuora LI. Angular diameter-redshift relations for extragalactic radio sources. Astrophysical Journal. 1982;260:37
  30. 30. Ubachukwu AA. The implications of intrinsic luminosity evolution on the value of the density parameter and the evolution of radio sizes of radio galaxies and quasars. Astrophysics and Space Science. 1995;228:195
  31. 31. Kapahi VK. Redshift and Luminosity dependence of the linear sizes of powerful radio galaxies. Astronomical Journal. 1989;97:1
  32. 32. Peebles PJ. The theoretical aspects of the nebular redshift, 33 years later. Publications of the Astronomical Society of the Pacific. 1988;100:670
  33. 33. Alhassan JA, Ubachukwu AA, Odo FC. On the absence of core luminosity-core-dominance parameter (P C - R) correlation in radio galaxies and BL lacs. Journal of Astrophysics & Astronomy. 2013;34:61
  34. 34. Fan JH and Zhang JS. JE The core dominance parameter of extragalactic radio sources. Astronomy & Astrophysics. 2003;407:899
  35. 35. Bell MB and Comeau. Evidence that the bursting component of the X-ray radiation from 3C 111 originates in the PC-scale jet. Astrophysics and Space Science. 2012;339:143
  36. 36. Onuchukwu CC, Ubachukwu AA, Onah CI. On the expansion speed of highly asymmetric extragalactic radio sources. Nigerian Journal of Space Research. 2014;12:39
  37. 37. Onah CI, Ubachukwu AA, Odo FC. Environmental effects on structural asymmetry in extragalactic radio sources and quasar/galaxy unification. European Physical Journal - Plus. 2021;136:576
  38. 38. Ubachukwu AA. Relativistic beaming and orientation effects in lobe-dominated quasars. Astrophysics and Space Science. 1998;257:23
  39. 39. Banhatti DG. Expansion speeds in extended extragalactic double radio sources from angular structure. Astronomy & Astrophysics. 1980;84:112
  40. 40. Onah CI, Ubachukwu AA, Odo FC. Asymmetries in powerful radio sources and the quasar/galaxy unification. Journal of Astrophysics and Astronomy. 2014;35:619. DOI: 10.1007/s12036-014-9311-z
  41. 41. Onuchukwu CC, Ubachukwu AA. Structural asymmetries, relativistic beaming and orientation effects in lobe-dominated quasars. Astrophysics and Space Science. 2013;344:211. DOI: 10.1007/s10509-012-1325-x
  42. 42. Bai JM, Lee MG. New evidence for the unified scheme of bl lacertae objects and fr i radio galaxies. American Journal of Physics. 2001;548:244. DOI: 10.1086/318695
  43. 43. Odo FC, Ubachukwu AA, Chukwude AE. On FRI-BL lac unification. African Skies/Cieux Africains. 2011;15:15-17
  44. 44. Odo FC, Ubachukwu AA. Linear size-extended radio luminosity (D- PE) correlation in BL Lacertae objects: evidence for large scale beaming? Astrophysics and Space Science. 2013;347:357. DOI: 10.1007/s10509-013-1528-9
  45. 45. Vermeulem RC, Cohen MH. Superluminal motion statistics and cosmology, Astrophysical Journal. 1994;279:93
  46. 46. Nilsson K. Kinematical models of double radio sources and the unified scheme II. The database. Astronomy & Astrophysics. 1998;132:31
  47. 47. Landt H, Padovani P, Perlman ES. VLA observations of a new population of blazars. Astrophysical Journal. 2006;637:183. DOI: 10.1086/498261
  48. 48. Barthel PD, Miley GK. Evolution of radio structure in quasars: A new probe of protogalaxies? Nature. 1988;333:318. DOI: 10.1038/333319a0
  49. 49. Gopal-Krishna and Wiita PJ. The expansion and cosmological evolution of powerful radio sources. Monthly Notices of the Royal Astronomical Society. 1987;226:531-542. DOI: 10.1093/mnras/226.3.531
  50. 50. Ubachukwu AA, Onuora LI. Radio source orientation and the cosmological interpretation of the angular size–redshift relation. Astrophysics and Space Science. 1993;209:169-180
  51. 51. Kollgard RI, Wardle JF, Roberts D, Gabuzda D. Astrophysical Journal. 1992;104:1687
  52. 52. Singal AK. In: Ekers R et al., editors. Extragalactic Radio Sources, IAU Symp. 175. Kluwer: Dordrecht; 1996. p. 563

Written By

Costecia Ifeoma Onah, Augustine A. Ubachukwu and Finbarr C. Odo

Submitted: 23 March 2022 Reviewed: 04 July 2022 Published: 06 August 2022