Next Article in Journal
Shaping Dynamical Casimir Photons
Next Article in Special Issue
Project MOMO: Multiwavelength Observations and Modeling of OJ 287
Previous Article in Journal
New Black Hole Solutions in N = 2 and N = 8 Gauged Supergravity
Previous Article in Special Issue
On the Heating of AGN Magnetospheres
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Jet-Induced Feedback in the [O III] Lines of Early Evolution Stage Active Galactic Nuclei

1
Finnish Centre for Astronomy with ESO, University of Turku, Vesilinnantie 5, FIN-20500 Turku, Finland
2
Aalto University Metsähovi Radio Observatory, Metsähovintie 114, FIN-02450 Kylmälä, Finland
3
European Space Agency (ESA), European Space Astronomy Centre (ESAC), Camino Bajo del Castillo s/n, Villanueva de la Cañada, 28692 Madrid, Spain
*
Author to whom correspondence should be addressed.
Universe 2021, 7(6), 188; https://doi.org/10.3390/universe7060188
Submission received: 22 January 2021 / Revised: 20 May 2021 / Accepted: 2 June 2021 / Published: 6 June 2021
(This article belongs to the Special Issue Panchromatic View of the Life-Cycle of AGN)

Abstract

:
It is well known that active galactic nuclei (AGN) show various forms of interaction with their host galaxy, in a number of phenomena generally called AGN feedback. In particular, the relativistic plasma jets launched by a fraction of AGN can strongly affect their environment. We present here a study of the [O III] λ λ 4959,5007 lines in a diverse sample of early evolution stage AGN–specifically narrow-line Seyfert 1 galaxies. Radio imaging observations of all of the sources enable a division to jetted and non-jetted sources, and exploiting this we show that the ionized gas properties are significantly influenced by the presence of the jets, as we often find the [O III] lines (blue-)shifted with respect to their restframe wavelength. We also show how the radio morphology and the radio spectral index do not seem to play a role in the origin of the [O III] shifts, thus suggesting that the source inclination is not relevant to the lines displacement. We do not find a strong relation between the [O III] line properties and the bolometric luminosity, suggesting that within our sample radiatively driven outflows do not seem to have a significant contribution to the [O III] line kinematics. We finally suggest that [O III] shifts may be a good proxy to identify the presence of relativistic jets. Additional studies, especially with integral-field spectroscopy, will provide a deeper insight into the relation between jets and their environment in early evolution stage AGN.

1. Introduction

Active galactic nuclei (AGN) play an essential role in the evolution of galaxies across cosmic time [1], as the central black hole (BH) and its host galaxy undergo co-evolution [2]. The nucleus can experience multiple activity phases throughout its life cycle, and during every phase there is interplay between the AGN and its host galaxy. The latter regulates the gas supply to the nucleus, enabling or hampering the triggering of the AGN, and inducing changes to its activity state. On the other hand, the AGN also affects its environment. The radiation produced by the process of accretion onto the BH can heat or ionize the gas up to several kiloparsecs (kpc) away (e.g., [3,4,5]), and by means of outflows or relativistic jets it can suppress or in some cases enhance [6] the star formation activity across the host galaxy [7]. All the examples of mutual interaction between an AGN and its host galaxy are collectively known as AGN feedback.
Feedback is still a relatively poorly understood aspect of galaxy evolution, and it is particularly interesting to study if and how this phenomenon occurs in AGN in early phases of their life cycle, that is, when they are “young”. Carrying out such research at high redshift, when all sources were young, is not an easy task. An alternative approach is to focus on sources that are the low redshift analogs to the first quasars. Narrow-line Seyfert 1 (NLS1) galaxies are in this sense ideal laboratories to study the early phases of AGN feedback since, like high-redshift quasars, they are characterized by fast growing black holes with high Eddington ratios, and also by high metallicity, especially prominent Fe II multiplets [8].
NLS1s are mostly non-jetted sources, but some of them do actually harbor powerful relativistic jets [9,10,11,12,13]. Typically hosted by spiral galaxies [14,15,16,17,18,19,20,21,22,23,24,25], due to their low-mass BHs and high Eddington ratios NLS1s are believed to be the progenitors of other AGN classes, such as broad-line Seyfert 1 galaxies and flat-spectrum radio quasars [8,26,27,28,29,30,31,32,33]. Since their classification is exclusively based on the width of the H β line [34], NLS1s can currently be found only up to z∼1 [35], although unclassified objects with similar properties are present also at higher redshift [36,37]. A feature which is often observed in the optical spectra of this class of AGN is a shift, usually toward lower wavelength, in the emission lines of [O III] λ λ 4959, 5007, corresponding to a velocity | v | > 150 km s 1 [38]. Such sources, in the literature, are called blue outliers [38,39,40]. They are mainly found among NLS1s but, more in general, they are observed almost exclusively in population A sources on the quasar main sequence (MS). The MS is the locus on the plane defined by the flux ratio between the Fe II multiplets and H β , known as R4570, and the FWHM(H β ), where all type 1 AGN lie. This sequence was identified by means of principal component analysis, and it seems to be mainly driven by the Eddington ratio [41]. However, some other factors seem to play a role in it, such as metallicity, and the inclination compared to the line of sight [42,43]. Sources on the MS can be divided into two main populations according to the full width at half maximum (FWHM) of their H β line. Population A sources are characterized by FWHM <4000 km s 1 , and they often exhibit prominent Fe II emission [27,36,40,44,45,46]. The origin of blue outliers is still unclear [47,48,49], and it could be connected to the high Eddington ratio observed in population A sources and NLS1s [50,51,52,53]. However, their presence is often associated with strong radio emission indicative of relativistic jets (hereafter B16, [48,54]). If this is true, blue outliers may be a remarkable example of feedback between kinematically young relativistic jets, as those present in NLS1s, and the ionized gas of their host galaxy.
In this work, we analyzed the properties of the [O III] lines, with particular attention to blue outliers, in a sample of NLS1s whose radio morphology has been observed with the Karl G. Jansky Very Large Array (JVLA), to investigate their relation with radio emission. In Section 2 we describe our sample selection and the data analysis process, in Section 3 we present our results, in Section 4 we discuss our findings, and in Section 5 we present a brief summary of this work. Throughout the paper, we adopt the standard ΛCDM cosmology, with a Hubble constant H 0 = 70 km s 1 Mpc 1 , and ΩΛ = 0.73 [55].

2. Sample Selection and Data Analysis

Our sample includes the 74 NLS1s that were observed with the JVLA at 5 GHz in its most extended A configuration (longest baseline 36 km, resolution 0.5”) in the survey described in (hereafter B18, [56]. We also decided to include into our sample seven additional NLS1s studied in the same JVLA configuration. These peculiar objects show mJy-level or lower flux density at 5 GHz, but they were observed flaring at Jy-level at 37 GHz, indicating that they harbor relativistic jets [57]. The reason for the weak radio emission at low frequencies may be a form of absorption, either free-free or synchrotron self-absorption [58].
We were able to retrieve an optical spectrum for 80 out of 81 sources (all but J0100-0200). In another object, J1348+2622, the [O III] lines are not visible in the optical spectrum due to its high redshift (z = 0.917), and it was therefore excluded from this analysis. For 62 of these NLS1s, a measurement of the [O III] line properties was already available in B16. For the remaining 17 objects, 13 spectra were obtained from the Sloan Digital Sky Survey (SDSS) Data Release 16 (DR16) [59], while the remaining four were observed with the 1.22 m telescope of the Asiago Astrophysical Observatory. The details of these new observations are reported in Table 1. The measurement process we used is the same as the one described in detail in B16. In the following we provide a brief summary.
The first challenge was obtaining a precise estimate of the redshift of each source. Ideally, stellar absorption lines would be the best proxy for the recessional velocity of a galaxy, but in type 1 AGN they are typically not visible. Emission lines are easily detectable, but not all of them are reliable redshift indicators. Broad lines and high-ionization lines often originate from outflows and appear shifted with respect to the galaxy frame, so the ideal solution is to use low-ionization lines, such as [S II] λ λ 6717, 6731, [O I] λ 6300, or [O II] λ 3727 [38,61]. It is worth noting that the [O II] line is an unresolved doublet in low-resolution spectra, therefore the other lines, if available, grant a better estimate of the redshift. However, due to the redshift, [O II] is often the only available reference. We measured the wavelength of the reference line by using a fitting procedure in Python. In this procedure there are two sources of systematic error. The first one is due to wavelength calibration errors. For SDSS its value is 2 km s 1 [62], while for the other spectra it is of the order of 20 km s 1 (B16). The other error is produced by the fit of the lines, and we computed it using a Monte Carlo method. We added a Gaussian noise, whose amplitude was proportional to the noise in the AGN continuum measured around λ 5100, to the line. We then repeated the fitting procedure 100 times with different noise, getting a median value of redshift and a standard deviation. On average, the latter is rather small, ∼2 km s 1   as it is obtained by a simple fit with one Gaussian of a relatively strong emission line.
After moving the spectra to restframe, we subtracted the AGN continuum using IRAF and, when necessary, we removed the Fe II multiplets using an online software 1 [63,64]. At this point, we decomposed each of the [O III] lines using two Gaussian components, one to reproduce the line core, and one for the blue-shifted wing of the line. The line core component originates from the ionized gas of the narrow-line region (NLR). Its kinematics is dominated by the gravitational potential of the bulge, and therefore it usually has the same recessional velocity as the host galaxy. The wing component, instead, probably comes from a turbulent outflow in the inner part of the NLR. Since in type 1 AGN we preferentially have a pole-on view of the central engine, the outflow appears as a blue-shifted (approaching) component in the line with a FWHM larger than that of the core due to its turbulence. Our aim was obtaining the position of the [O III] core component to derive the velocity of the ionized gas with respect to the galaxy frame, the position of the line wing to measure the velocity of the outflow, and the FWHM of both components to derive the gas velocity dispersion. To reduce the number of free parameters in the fitting procedure, we fixed the ratio between each component of the λ 5007 and λ 4959 lines to its theoretical ratio of 3.0 [65]. Furthermore, we forced the FWHM of the line core, the FWHM of the wing, and the wing velocity, expressed in km s 1 , to be the same in both lines. The errors were computed adopting the previously described Monte Carlo procedure. In conclusion, for each line we measured four parameters: the core shift with respect to the laboratory wavelength of [O III] (5006.843Å), the line core FWHM, the wing velocity with respect to the core component, and the wing FWHM. An example of a fit is shown in Figure 1. We successfully reproduced the [O III] lines in 76 sources, while the remaining three did not show any oxygen in their spectra (see Discussion). The results of our calculation are shown in Table 2. In the following, we also considered some additional physical properties of the sources such as their bolometric luminosity, radio luminosity, and radio loudness parameter. The bolometric luminosity was estimated using its scaling relation with the H β luminosity, as in [66]. The radio luminosities are derived from B18 and [58]. The radio loudness parameter is defined as the ratio between the radio flux density at 5 GHz and the flux density in the optical B band. When the ratio is larger than 10, the source is radio loud. Conversely, it is radio quiet [67]. To estimate it, we either used the values published from [13,66], or we calculated it using the B-band magnitude values reported in the NASA Extragalactic Database (NED) and our own measurements at 5 GHz.

3. Results

3.1. Correlations with the Whole Sample

The property we are most interested in is the shift of the whole [O III] line, that is, the blue outlier phenomenon. Our aim is to have a better understanding of what may be causing this shift. We tested the core velocity against four physical properties that may be associated with it, that is the actual [O III] line luminosity, the bolometric luminosity, the radio luminosity, and the radio loudness parameter. The first two, indeed, have been found to possibly affect the [O III] core shift especially at high redshift [68]. The radio luminosity and the radio loudness parameter, instead, may be related with the presence of relativistic jets (with some caveats, see e.g., [69,70]). The results of these tests are reported in Table 3. None of the tested quantities is strongly correlated to the velocity of the [O III] core. This is true also when the absolute value of the shift is considered, neglecting whether the line is blueshifted or redshifted. An example is shown in Figure 2. Since it is known that [O III] shift can also be radiatively driven in high luminosity sources [68], and a weak, but significant correlation is identified by the Spearman’s rank correlation coefficient between L bol and |v c |, we wanted to investigate whether an L bol threshold exists above which a source has a higher possibility of showing an [O III] shift or a higher FWHM([O III]). As seen in the left panel of Figure 2, however, no clear threshold seems to be present to produce an [O III] shift. In the FWHM distribution, the highest values of FWHM are found mostly, but not exclusively, among the highest luminosity sources, although the statistic is too poor to draw any meaningful conclusion.
It is also worth mentioning that the Spearman’s rank correlation coefficient has a rather low p-value when the radio loudness parameter and, even more noteworthy, the radio luminosity are tested against the absolute value of the core shift velocity, suggesting that these two weak correlations may be significant. The radio luminosity, however, contains some dependence on the redshift. Indeed, B18 already identified a correlation between the radio luminosity and the redshift, present also in our sample (r P = 0.51, p-value = 3 × 10 6 ), which is due to the sample selection. For this reason, we carried out a partial correlation analysis, using the redshift as a control parameter. In this case, the Spearman’s rank correlation coefficient becomes 0.26 (p-value = 0.03), so the correlation seems to disappear. We carried out the same analysis on the other correlations involving the radio luminosity shown in Table 3, but we found no statistically significant result.
The radio loudness parameter instead is redshift independent, and it also shows a weak but statistically significant correlation with the [O III] core shift. As we already mentioned, radio loudness is often used as a proxy to detect the presence of relativistic jets. This, therefore, suggests that sources harboring relativistic jets may tend to present [O III] core shift more often than those without a jet. To verify this hypothesis, we decided to split the sample of sources according to the presence or not of a relativistic jet and of other radio properties derived in B18, and to carry out some statistical tests to verify if and how these features affect the optical lines’ profile.

3.2. Radio Classification

The first radio-based classification, jetted (J) or non-jetted (NJ), is based on the data from B18, where each source was analyzed individually. In general, the classification criteria were the radio luminosity, which is typically higher than 10 40 ergs s 1 in jetted NLS1s, the spectral index, since a flat radio spectrum with α ν < 0.5 indicates the presence of relativistic beamed jets, and the multiwavelength properties of the sources. In previous studies as B16, each source could be classified only according to the radio-loudness parameter [67], defined as the ratio between the flux density at 5 GHz and in the optical B band. Sources with radio-loudness larger than 10 are classified as radio-loud, conversely as radio-quiet. Objects with no previous detection in radio surveys are often labeled as radio-silent. Usually, radio-loud sources harbor relativistic jets, while the origin of radio emission in radio-quiet objects is strongly debated (e.g., [71]). However, several sources do not obey this rule, therefore a classification criterion based on physical properties is needed [69]. Based on the radio observations, B18 established whether relativistic jets are present or not in each source. An example is J1038 + 4227, which was classified as radio-loud by Foschini et al. [13] (and so, estimated to be jetted), but was later found to not have relativistic jets, since the radio emission originates from prominent star forming activity [72].
The second subdivision is based on the radio morphology of our NLS1s. B18 classified them with a compactness criterion. They used the ratio R = S p / S i , that is the ratio between the peak flux density and the integrated flux density of the whole emitting region. Sources with R 0.95 were classified as compact (C), those with 0.75 R < 0.95 were labeled as intermediate (I), and finally those with R < 0.75 were classified as extended (E). While jetted sources were mostly (but no exclusively) compact, non-jetted sources typically had an extended morphology. The morphological parameters were calculated using the flux densities reported in B18 and [58]. For the sources without a detection in [58], we assumed a compact morphology. It is worth noting that since the radio observations were carried out in the most extended configuration of JVLA, some large-scale extended emission may had been resolved out. However, a sample of NLS1s were investigated at the same frequency by Chen et al. [73] using the JVLA in a more compact configuration. Their results seem to indicate that the morphology of jetted sources is typically compact, with no signs of large-scale radio lobes as those seen in other classes of jetted AGN. Furthermore, the spatial scale of these structures even in A configuration would be several kpc and possibly outside of the host galaxy.
The last classification was used only for jetted sources, and is based on the division to sources with flat (F) or steep (S) radio spectral index α ν ( F ν ν α ν ). As in other AGN, in jetted NLS1s the radio spectral index is a rather good indicator of the source inclination, that is the angle between the line of sight and the relativistic jet axis. In most AGN, the spectrum of the radio core appears flat due to the τ = 1 surface, where τ is the optical depth, moving closer to the BH as a function of frequency [74]. When our line of sight is close to the axis of the system, the emission of the radio core experiences relativistic beaming, that can result in the radio core emission to dominate the whole spectrum. If we are seeing the radio core at larger angles, the relativistic effects diminish and the contribution of extended, often optically thin, emission becomes more important. A properly formed, powerful jet is often characterized by several components, launched at different times, that each have a convex spectrum, that is synchrotron self-absorbed below the turnover frequency, ν t , and shows an inverted spectral index ( α ν −2.5), and optically thin with a spectral index of α ν 0.7 above ν t . As a component ages it loses energy and expands, and the peak moves to lower frequencies. A superposition of several components with different ages results in a seemingly flat, but variable, radio spectrum [75]. This is especially pronounced in sources whose jet axis is close to our line of sight, since the emission is enhanced by relativistic beaming, and in these sources the beamed jet and radio core emission dominates the whole radio spectrum. In sources that we observe at larger angles, the beaming effects become less significant, and whereas the radio core still appears flat, the overall radio spectrum becomes dominated by extended, optically thin, emission, with a steeper spectral index. The spectral classification used in this work is derived from B18 and Järvelä et al. (in prep.). All the classes are reported in Table 4.

3.3. Statistical Tests

We studied the properties of the [O III] lines by grouping them according to the radio properties as described in the previous section. For each distribution, we estimated its median value, its standard deviation, and its interquartile range (IQR). We remark that the standard deviation is not used here as an estimate of the error associated to the median, but instead as a proxy of how much each distribution is concentrated around its median value. Furthermore, to examine whether two distributions are statistically different or not, we used both the Kolmogorov–Smirnov (K-S) and Anderson–Darling (A-D) test. The A-D test is more powerful than the K-S in most cases [76], but we decided to adopt both tests to provide results comparable to the literature and B16. The null hypothesis of these tests is that the two distributions originate from the same population of sources. In the following, we will reject it only if p-value ≤ 0.05. Finally, to account for the measurement errors, in addition to the observed samples we created simulated samples of sources. In each distribution, to each measurement we added a Gaussian noise whose amplitude was equal to the measurement error, thus creating new samples with the same size as the observed ones. We repeated this operation 100 times and performed every time both statistical tests, to obtain a median p-value for both K-S and A-D test. The goal of these simulations is to verify the stability of the p-values we measured on the observed samples. Therefore, we also estimated the standard deviation of the distributions of p-values obtained from the simulations.

3.4. Jets and [O III]

3.4.1. Line Core Component

First of all, we examined the properties of the [O III] lines separately in jetted and non-jetted sources, as indicated in the second to last column of Table 4. We first compared the velocity of the core component of [O III] with respect to the galaxy rest frame in the two samples. We remark that, by definition, a source is considered a red or blue outlier when the velocity of the line core component is | v c | > 150 km s 1 , following the definition by Komossa et al. [38]. The distributions are shown in the left panel of Figure 3. Among our 28 non-jetted sources, we found only one blue outlier (∼4%) and no red outliers. In our sample of 48 jetted sources, instead, we found 15 blue outliers (∼31%) and 3 red outliers (∼6%).
The statistical properties of the two distributions are shown in Table 5. It is evident that while the median values of the two samples are similar, both the standard deviation and the IQR are very different. Both quantities are much larger in jetted sources, suggesting that the [O III] core velocity distribution in this subsample is significantly more spread out. Despite this difference, the K-S test indicates that the null hypothesis cannot be rejected, but the A-D test provides the opposite result (p-value = 0.01). We repeated the tests on the simulated samples, and the results did not change significantly. However, since our goal is to understand the nature of blue outliers, and they are concentrated in the distributions’ tail, the A-D test is more reliable than the K-S test. We can therefore conclude that the null hypothesis can be rejected, and that the core shift distributions do not originate from the same population of sources.
We then tested the distributions of FWHM of the [O III] core component in the two samples, shown in the right panel of Figure 3. As we did before, we tested the differences between the two samples, and also in this case the distribution of jetted sources has larger standard deviation and IQR. In this case, both statistical tests suggest that the null hypothesis can be rejected in the observed samples, and also the simulated ones. Therefore, we can conclude that the two distributions of [O III] core FWHM do not originate from the same population of sources.
We also studied whether a correlation is present among the core shift and the FWHM of the [O III] core component. The result is shown in Figure 4, while the values we obtained are reported in Table 6. The Pearson correlation coefficient for all sources is r = 0.70 , and its p-value, that is the probability that this result is due to a random distribution of sources, is 8 × 10 13 . However, as we have seen before, a core shift is present almost exclusively in jetted sources. Therefore, we decided to test the two samples separately. In non-jetted sources, no correlation is present ( r = 0.02 , p-value = 0.93), indicating that the FWHM of the [O III] core and its position are independent of each other. Conversely, in jetted sources the Pearson correlation coefficient is higher ( r = 0.73 ), and the p-value is 4 × 10 9 , indicating that the correlation is present only in this sample of objects. It is worth noting that this correlation is utterly dominated by the presence of the blue outliers. Indeed, when the Pearson coefficient is evaluated in this sample rejecting the blue outliers, its value becomes r = 0.28 (p-value = 0.13). Vice versa, when only the blue outliers are tested, the correlation remains, although with a lower p-value (r = 0.67 , p-value = 3 × 10 3 ). In conclusion, [O III] core velocity and FWHM do not seem to correlate in the bulk of NLS1s population, but they are related among blue outliers.

3.4.2. Wing Component

As previously mentioned, the wing component is likely associated to an outflow coming from the inner part of the NLR. A wing component could be identified during the decomposition process in all but ten sources, three non-jetted and seven jetted. We tested the distribution of wing velocity in jetted and non-jetted sources, shown in the left panel of Figure 5, and whose parameters are summarized in Table 5. As for the distributions of the [O III] core properties, the wing velocity distribution is somewhat peaked in non-jetted sources, with a large majority of sources having velocities between −400 and −200 km s 1 , while it is more spread out in the jetted sample.
However, the statistical tests seem to indicate that the null hypothesis cannot be rejected with the required confidence level. The same is true for the distributions of wing FWHM, displayed in the right panel of Figure 5. The presence of relativistic jets does not ultimately seem to affect the FWHM of the [O III] wing, which remains as turbulent in the jetted NLS1s as in the non-jetted NLS1s. However, it is worth noting that the standard deviation of the p-value derived in the simulations is rather large because of the relatively large measurement errors on the wing properties. Therefore, our data may not be enough to reach a robust conclusion.
Finally, to better understand if and how the physical quantities we measured are related to each other, we tested various correlations by evaluating the Pearson coefficients. The results are shown in Table 6. We first tested the presence of a correlation between the wing velocity and both [O III] core velocity and FWHM velocity, but we do not find any strong correlation, except for a trend between the wing velocity and the [O III] core FWHM in jetted sources (r = −0.50, p-value = 8 × 10 4 ). In a similar fashion, the FWHM of the wing component does not appear to correlate with any property of the [O III] core component, except again for the [O III] core FWHM in jetted sources (r = 0.54, p-value = 2 × 10 4 ). We finally tested the presence of a correlation between the wing velocity and its FWHM, and we find another possible trend among these two quantities in non-jetted sources (r = −0.51, p-value = 0.01).

3.5. Radio Morphology and [O III]

We tried to investigate whether the radio morphology is somehow connected to the [O III] properties. As before, we first analyzed the [O III] core velocity distributions, which are shown in the left side of Figure 6. Outliers do not seem to appear preferentially in any specific sample. Among the 26 C sources, we found one red outlier and six blue, among 27 I sources there is one red outlier and six blue, and finally in E sources we found one red outlier and four blue. To check whether the distributions were different from each other, we used the K-S and A-D tests both on the observed samples and on the simulated samples. All tests, even when accounting for the p-values distributions, confirm that the null hypothesis can never be rejected. The results of all tests, as well as the median velocity, its standard deviation, and the IQR are all reported in Table 7.
We also tested the distributions of FWHM of the [O III] core component across the different morphological samples. The results are shown in the right side of Figure 6. As for the [O III] core velocity, however, this parameter does not seem to be related to the radio morphology. All the tests we carried out, and that are summarized in Table 7, seem to point out that the distributions of FWHM are essentially the same across the three samples. Indeed, median values, standard deviation and IQR are all very similar.
There is only one minor detail possibly worth noting, that is the blue outliers in the E sample are those with velocity closest to the threshold. Indeed, the [O III] core velocity distribution of E sources has the lowest standard deviation and IQR. Since most sources with extended morphology do not have relativistic jets, this result could be a (weak) reflection of the difference between the distributions of jetted and non-jetted sources we found in the previous section. However, the samples we are examining do not allow us to draw any statistically robust conclusions.
The behavior of [O III] wings, as for the core component, does not seem to be related in any way to the radio morphology. We tested the wing velocity and FWHM among the three samples, and found that there is no evidence against the null hypothesis. Detailed results are shown in Table 7.

3.6. Effect of the Redshift

Before moving any further, we tried to quantify the impact that the redshift could have on some possibly redshift-dependent properties. The redshift distribution of jetted and non-jetted sources is indeed rather different, as shown in the left panel of Figure 7, with jetted sources detected up to higher redshift. As shown in B18, the same effect is present in the morphological samples, since sources with extended morphologies are usually not detected at high redshift due to their low surface brightness densities. We do not expect significant cosmological evolution of the sources within the redshift range of our samples, both jetted and non-jetted. Indeed, in the radio luminosity function of jetted NLS1s, studied by [30], no evolution is observed up to z = 0.6. Only five out of 48 of our jetted sources lie above this threshold, therefore we do not expect them to affect significantly our results. Furthermore, even if some galaxy evolution were present, it is unlikely to affect the interaction between the relativistic jets and the NLR, since the physical mechanisms at play remain the same.
We tried also to verify whether the distribution of [O III] luminosity is different among the jetted and non-jetted sources. Indeed, since jetted sources are located at higher redshift, there could be a selection bias toward more luminous sources. The two distributions are shown in the right panel of Figure 7, and they are actually different, as both K-S and A-D test show p-values 10 4 . However, we tried to verify whether the [O III] luminosity within the jetted sample is dependent on the redshift, and we found that only a weak correlation is present (r = 0.39, p-value = 0.01). It is also worth noting that, at z < 0.6, even this weak correlation disappears. Sources with bright [O III], therefore, are not necessarily located at higher redshift. This suggests that the different luminosity observed in jetted and non-jetted sources is not a selection effect, but a real physical difference between the two classes. This should not come as a surprise, since the presence of a correlation between the [O III] luminosity and the radio luminosity, which in turn is related to the presence of a relativistic jet, has been known since the seventies [77].
We investigate whether it is possible that the higher [O III] luminosity instead of the relativistic jets affects the presence of blue outliers. A higher line luminosity could originate because of a higher disk luminosity which, in turn, may induce outflows via radiation pressure, and disturbed and asymmetric profiles are often observed in the [O III] and C IV lines of high-luminosity sources [47,68,78]. However, this may not be the case for our sources. As found already by B16, and shown by the correlation coefficients of Table 3 and Table 6, there is no evidence that the [O III] core shift correlates with the [O III] line luminosity, neither when the whole sample nor when the jetted and non jetted subsamples are considered. We also investigated the redshift dependency of the bolometric luminosity. Naturally, a strong correlation is present (see Figure 8) with higher luminosity sources seen at higher redshift (Spearman’s r s = 0.67, p-value 3 × 10 11 ). However, as seen in Section 3.1, the bolometric luminosity also does not appear to correlate with the [O III] line properties. Therefore, the presence of the relativistic jets remains the most likely explanation for blue outliers and the other observed properties, at least for the sources included in our sample.

3.7. [O III] vs. Spectral Index

Testing the [O III] properties as a function of the spectral index will allow us to understand, at least at a first approximation, whether they are dependent on the source inclination or not.
The resulting distributions of jetted NLS1s are shown in the left panel of Figure 9. The number of blue and red outliers in the two samples is essentially the same. Out of 21 steep-spectrum sources, we found two red outliers and six blue, while out of 27 flat-spectrum sources we found one red outlier and nine blue. Neither the K-S nor the A-D test allows us to reject the null hypothesis that the two distributions of core shift velocities originate from the same sample of sources. The p-value is always above 0.9, both in our samples and in the simulated ones, indicating that the two distributions are very similar. Also median velocity, standard deviation and IQR are comparable among the two samples. We also analyzed the distributions of FHWM of the [O III] core component, shown on the right side of Figure 9. As for the core shift distributions, the null hypothesis can never be rejected. In a similar way, the properties of the [O III] wings among the two samples, both their velocity and FWHM, show that they originate from the same population of sources. All the values we measured are shown in Table 8.
It is worth noting that we decided to include the seven sources with absorbed radio jets in the flat-spectrum sample even if, formally, at low radio frequency, they either show a steep spectrum, or they are not detected. The reason for this is that their rapid variability observed at high radio frequency indicates that they harbor relativistic jets aligned with the line of sight, as seen in blazars. In fact, their spectral indexes are steep because the radio emission at low frequency is dominated by star formation, but their intrinsic spectrum, if corrected for absorption, has to be flat. However, our conclusions are independent of this assumption. Even when all of them are included in the steep-spectrum sample, the statistical tests still do not provide any evidence against the null hypothesis. This indicates that our conclusion is rather robust, and that the radio spectral index is not directly connected to the physical properties of the [O III] line.

4. Discussion

4.1. Origin of the Bulk Outflows

The [O III] lines of our sample could be modeled with two components in 67 out of 77 sources in which the line is detected (87%), a fraction which is in agreement with previous findings in the literature (e.g., [79,80]). In type 2 and intermediate type AGN, as well as in star forming galaxies, instead, the fraction of two-components [O III] lines is significantly lower [81,82]. This indicates that in type 1 AGN two separate kinematic components are directly visible. The first one is the [O III] core, which is usually not moving with respect to the galaxy restframe, and the second one is the wing, associated to an outflowing and more turbulent component closer to the central engine [83].
In some cases both the [O III] core and wing components are shifted toward bluer wavelengths, although outflows in the opposite direction can also be observed, likely depending on the geometry of the source. The presence of such wavelength shift is interpreted as a sign of a bulk outflowing motion in their NLR. The existence of these outflows in the NLR of AGN has been known for a long time, and it is particularly common in NLS1s and other population A sources [39,41,84]. In particular, Ganci et al. [49], analyzing high-quality optical spectra, found that line displacements are often found in sources with prominent Fe II multiplets, a feature which is typically observed in NLS1s. This spectral property has often been associated to the high Eddington ratio that is characteristic of population A sources, but past studies showed that these two quantities are not directly correlated (B16, [80,85]). An alternative explanation is that the origin of this phenomenon resides in the relativistic jets crossing the ISM and interacting with it [38]. Such phenomenon, indeed, has been observed not only in NLS1s, but in several galaxies harboring compact jets [86].
Our results seem to support the latter scenario. As seen in the left panel of Figure 3, the number of blue (and red) outliers observed in sources with relativistic jets is significantly higher than what is found in the non-jetted sample. This finding is also in agreement with the main result of B16. In that work, however, the adopted radio classification was based on the radio-loudness criterion, which has been repeatedly proven to be rather ineffective [33,57,58,69,70,72,73]. In this work, instead, we make use of a more rigorous classification based directly on the observed radio properties and morphology of the sources. The relativistic jets of NLS1s are usually rather compact [73,87,88,89], although some exceptions are known (e.g., [90,91,92,93,94,95]). This has led some authors to suggest that NLS1s may be connected with classes of young radio galaxies such as compact steep-spectrum (CSS) sources and gigahertz-peaked sources (GPS) [38,89,96,97,98,99,100,101,102]. Indeed, several radio galaxies show signs of line shifts and broadening induced by the presence of jets [40,86,103,104], and this seems particularly common in young radio galaxies, where the jet is directly interacting with the ISM [105,106]. In this framework, the detection of a large fraction of outliers among jetted NLS1s is not surprising.
As previously mentioned, Ganci et al. [49] found that line displacements are commonly, although not exclusively, present in sources with prominent Fe II emission. Our data do not allow to test this conclusion, because accurate Fe II measurements would require high signal-to-noise (S/N) ratio spectra which are currently unavailable. On the other hand, their sources are not prominent radio emitters according to the FIRST and NVSS radio surveys at 1.4 GHz [107,108]. Potentially, this could indicate that they do not harbor relativistic jets, and that two different physical mechanisms causing outflows can (co-)exist, a first one connected to the relativistic jet, and a second one connected to the presence of iron multiplets. In particular, the iron abundance may be connected to a high Eddington ratio, as postulated by the quasar main sequence [37]. This property can produce strong outflows that can ultimately produce the NLR bulk motion we observed. Indeed, Marziani et al. [68] observed that at intermediate redshift (0.9 < z < 2.2), blue outliers seem to be significantly more common than what is seen at low redshifts, with a whopping detection rate of ∼40%. Those quasars, however, have bolometric luminosity in the range 47 < log L Bol < 48.5 erg s 1 , that is one order of magnitude more than the brightest source included in our sample. We therefore speculate that, at such extreme luminosity, radiation pressure driven outflows are the most likely mechanism to produce blue outliers. This, of course, preferably happens at relatively high redshift, since in the local Universe such extremely bright objects are usually not observed. When, instead, the bolometric luminosity is smaller, it is possible that the formation of blue outliers is driven mostly by the presence of relativistic jets, although the radiation pressure channel cannot be fully ruled out, since we do observe outliers among non-jetted AGN.
It is worth noting that a weak radio emission at low frequencies does not necessarily rule out the presence of relativistic jets. As found by Berton et al. [58], the radio emission of some small-scale jets can be entirely absorbed via free-free or synchrotron self-absorption, thus making these features invisible in the low-frequency domain (<10 GHz). Their interaction with the ISM is instead well visible, since two out of seven of these peculiar objects have shifted [O III] lines (J1232+4957, J1509+6137). Therefore, we suggest that blue outliers, especially those showing the most prominent line displacements, could be potential targets in which an absorbed jet could be present. Only high-frequency radio observations will allow us to clarify this aspect.
As in Komossa et al. [48], we found a rather strong correlation between the FWHM and the shift of the [O III] core component (Figure 4). This seem to indicate that a bulk motion in the NLR (the core shift) is characterized by a larger intrinsic turbulence of the gas (the [O III] core FHWM). This result was observed not only in [O III] lines, but also in other NLR lines (such as [Fe VII], [Ne III], [Ne V], [Fe X]), and in the BLR line C IV λ 1549 [38,80,109,110]. In our analysis, this correlation is dominated by the outliers found in the jetted sources, again suggesting that this behavior can be produced by the relativistic jets crossing the NLR. Another effect of this interaction may be seen in the wing component, since its velocity distribution is different in jetted and non-jetted sources, with the former having a larger scatter than the latter.
Another interesting result of our analysis is the non-detection of the [O III] line in three of our sources, J0744+5149, J0945+1915, and J1421+2824. Except for J0744+5149, where the S/N ratio in the continuum is low (∼3), the other sources have relatively good spectra (S/N > 10). Therefore, the non-detection in these cases is not due to a poor data quality, but to an intrinsic weakness of the [O III] lines. All of these sources have relativistic jets with a compact morphology and a flat radio spectrum. This may not be a coincidence. The presence of relativistic jets can decrease the covering factor of the clouds where high ionization lines such as [O III] originate, and as a consequence reduce the equivalent width and luminosity of the lines [111,112]. Therefore, the non detection of these lines may be another sign of jet/NLR interaction.
The reason why all jetted sources are not blue outliers is still not fully understood. A possible explanation for this behavior, however, was provided in the simulations by Wagner et al. [113]. The interaction between the relativistic jets and the NLR can indeed occur only if the ratio between the total jet power and the Eddington luminosity of the source exceeds a threshold value (log(P j /L E d d ) > −4). When this condition is met, the velocity of the outflow produced by the relativistic jet of power P j exceeds the gas velocity dispersion derived from the relation between the black hole mass and the stellar velocity dispersion [114] of a galaxy for which the Eddington luminosity is L E d d . Therefore, the gas can escape the gravitational potential of the black hole. Unfortunately, our data do not allow us to test this hypothesis for all of our sources. However, for a few sources the jet power was calculated by Foschini et al. [13], and all of them respect this condition, as already shown by B16. In conclusion, we speculate that only the high-power jets launched by relatively low-mass-and hence low Eddington luminosity-sources, such as in NLS1s, can produce blue outliers.

4.2. [O III] and Radio Properties

A direct comparison between the radio morphology and the [O III] properties in high Eddington ratio AGN such as NLS1s has never been carried out, since the number of radio surveys aimed at studying the kpc-scale properties of these objects is very limited. The role of the radio morphology was studied by dividing the sources in compact, intermediate, and extended, based on the concentration of the flux density in the radio core. As found by B18, jetted sources often have compact or intermediate morphology, while non-jetted ones are mostly extended. However, our results seem to indicate that this is not connected to any of the [O III] properties we measured. The number of outliers found in the different morphological samples is basically the same, and also the other properties do not show preferences for any sample. This is somewhat surprising taking into account that jetted sources preferably have a compact morphology, and we therefore would have expected to see this trend reflected on the morphological distributions.
Interestingly, we also found that the [O III] properties in jetted sources are the same when the radio spectral index is flat or steep. As previously described, the spectral index can be used as a rough proxy to estimate the source inclination. Flat-spectrum objects are preferentially observed at low inclination with respect to the line of sight, while in steep-spectrum sources the relativistic jets are observed at larger angles. This suggests that the perturbations in the ionized gas kinematics can be seen regardless of the jet inclination. Since the structure of AGN has a cylindrical symmetry, also this result is slightly counterintuitive, as one would expect to observe outflows preferentially at low inclination [9,10], but it is not utterly surprising. If orientation were the sole responsible for blue outliers, we should see them in sources with the highest relativistic beaming and, as a consequence, the highest radio luminosity. However, this has never been the case, as clearly shown in the past [48,54]. Our work confirms this result. The Pearson correlation coefficient between the core shift velocity and the total radio luminosity at 5 GHz is 0.005, with a p-value of 0.97, providing no evidence that these two quantities are connected. Also no correlation was found in the jetted and non-jetted subsamples (see Table 6).
Determining the exact inclination of an AGN, especially non-jetted, is not an easy task. The most reliable technique relies on the geometry of the ionization cones, but it can be applied only for a limited number of sources [115]. Blue outliers, until now, were commonly considered low-inclination sources. Indeed, Boroson [116] suggested that they could be used, along with the red-shift of Fe II multiplets ([117], but see [118]), as inclination indicators. However, our findings seriously challenge this conclusion, since blue outliers seem to be uniformly distributed among sources with different inclinations.
The cloud acceleration process, therefore, seems to occur in every direction, and not only along the jet propagation axis. This behavior is similar to what was predicted by hydrodynamical simulations by Wagner and Bicknell [119] and Wagner et al. [113] 2. In their simulations, the jet can accelerate the clouds to velocities up to 10 3 km s 1 , which are comparable to the highest values we observed in our sample. Similarly, Cielo et al. [120] found that young jets can accelerate the gas clouds up to these velocities much more efficiently than what nuclear outflows induced by disk winds can do. Therefore, we can conclude that faster blue outliers are the most likely to originate from jet-ISM interaction.

5. Summary

In this paper, we analyzed the properties of the [O III] lines in a sample of low-mass high Eddington ratio AGN belonging to the NLS1 class, and compared them with their radio properties reported by a recent JVLA survey at 5 GHz. We found strong evidence that the presence of relativistic jets in these AGN can produce a bulk outflowing motion in the NLR, causing the [O III] lines to appear blue-shifted with respect to their restframe wavelength. Among the non-jetted AGN we found one outlier out of 28 sources, while amidst the jetted AGN we identified 18 outliers out of 48 sources. No strong correlation between the bolometric luminosity of the sources and the [O III] line properties was found, suggesting only a minor contribution of radiation pressure driven outflows in producing [O III] shifts. We therefore suggest that the most common origin for these so-called blue outliers is the interaction between the NLR gas clouds producing [O III] lines and the relativistic jets. The presence of blue outliers seems not to be connected either with the radio morphology of the source, or with the inclination of radio jets. This possibly indicates that the acceleration of the clouds produced by the propagation of the jets in the NLR medium does not have a preferred direction. The association between the blue outliers and the relativistic jets might indicate that the properties of the [O III] lines can be used as a powerful tool to identify the presence of jets whose low-frequency radio emission is absorbed, and that have gone undetected in the past radio surveys. Further studies, especially by means of spatially resolved integral-field spectroscopy, are needed to fully understand the interplay between the relativistic jets and the NLR clouds.

Author Contributions

Conceptualization, analysis and writing, M.B.; analysis and writing, E.J. All authors have read and agreed to the published version of this manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Most of the data used for this research are already publicly available in the literature. The non-public data can be made available upon reasonable request to the authors.

Acknowledgments

This work is based on observations made with the Galileo 1.22 m telescope of the Asiago Astrophysical Observatory operated by the Department of Physics and Astronomy “G. Galilei” of the University of Padova. 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 Admistration. Funding for the Sloan Digital Sky Survey has been provided by the Alfred P. Sloan Foundation, and the U.S. Department of Energy Office of Science. The SDSS web site is http://www.sdss.org (accessed on 5 June 2021). SDSS-III is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSS-III Collaboration including the University of Arizona, the Brazilian Participation Group, Brookhaven National Laboratory, Carnegie Mellon University, University of Florida, the French Participation Group, the German Participation Group, Harvard University, the Instituto de Astrofisica de Canarias, the Michigan State/Notre Dame/JINA Participation Group, Johns Hopkins University, Lawrence Berkeley National Laboratory, Max Planck Institute for Astrophysics, Max Planck Institute for Extraterrestrial Physics, New Mexico State University, University of Portsmouth, Princeton University, the Spanish Participation Group, University of Tokyo, University of Utah, Vanderbilt University, University of Virginia, University of Washington, and Yale University.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Most of the data used for this research are already publicly available in the literature. The non-public data can be made available upon reasonable request to the authors.

Notes

1.
http://servo.aob.rs/FeII_AGN/ (accessed on 5 June 2021).
2.
A very explanatory video can be found at https://www.youtube.com/watch?v=6-CO4hp6PYg (accessed on 5 June 2021).

References

  1. Taylor, P.; Kobayashi, C. The effects of AGN feedback on present-day galaxy properties in cosmological simulations. Mon. Not. R. Astron. Soc. 2015, 448, 1835–1846. [Google Scholar] [CrossRef]
  2. Kormendy, J.; Ho, L.C. Coevolution (Or Not) of Supermassive Black Holes and Host Galaxies. Annu. Rev. Astron. Astrophys. 2013, 51, 511–653. [Google Scholar] [CrossRef] [Green Version]
  3. Schmitt, H.R.; Ulvestad, J.S.; Antonucci, R.R.J.; Kinney, A.L. Jet Directions in Seyfert Galaxies: Radio Continuum Imaging Data. Astrophys. J. Suppl. Ser. 2001, 132, 199–209. [Google Scholar] [CrossRef]
  4. Congiu, E.; Contini, M.; Ciroi, S.; Cracco, V.; Di Mille, F.; Berton, M.; Frezzato, M.; La Mura, G.; Rafanelli, P. Extended Narrow-Line Region in Seyfert Galaxies. Front. Astron. Space Sci. 2017, 4, 27. [Google Scholar] [CrossRef]
  5. Venturi, G.; Nardini, E.; Marconi, A.; Carniani, S.; Mingozzi, M.; Cresci, G.; Mannucci, F.; Risaliti, G.; Maiolino, R.; Balmaverde, B.; et al. MAGNUM survey: A MUSE-Chandra resolved view on ionized outflows and photoionization in the Seyfert galaxy NGC1365. Astron. Astrophys. 2018, 619, A74. [Google Scholar] [CrossRef] [Green Version]
  6. Fragile, P.C.; Anninos, P.; Croft, S.; Lacy, M.; Witry, J.W.L. Numerical Simulations of a Jet-Cloud Collision and Starburst: Application to Minkowski’s Object. Astrophys. J. 2017, 850, 171. [Google Scholar] [CrossRef] [Green Version]
  7. Morganti, R. The many routes to AGN feedback. Front. Astron. Space Sci. 2017, 4, 42. [Google Scholar] [CrossRef] [Green Version]
  8. Mathur, S. Narrow-line Seyfert 1 galaxies and the evolution of galaxies and active galaxies. Mon. Not. R. Astron. Soc. 2000, 314, L17–L20. [Google Scholar] [CrossRef] [Green Version]
  9. Komossa, S.; Voges, W.; Xu, D.; Mathur, S.; Adorf, H.M.; Lemson, G.; Duschl, W.J.; Grupe, D. Radio-loud Narrow-Line Type 1 Quasars. Astrophys. J. 2006, 132, 531–545. [Google Scholar] [CrossRef] [Green Version]
  10. Yuan, W.; Zhou, H.Y.; Komossa, S.; Dong, X.B.; Wang, T.G.; Lu, H.L.; Bai, J.M. A Population of Radio-Loud Narrow-Line Seyfert 1 Galaxies with Blazar-Like Properties? Astrophys. J. 2008, 685, 801. [Google Scholar] [CrossRef] [Green Version]
  11. Abdo, A.A.; Ackermann, M.; Ajello, M.; Axelsson, M.; Baldini, L.; Ballet, J.; Barbiellini, G.; Bastieri, D.; Battelino, M.; Baughman, B.M.; et al. Fermi/Large Area Telescope Discovery of Gamma-Ray Emission from a Relativistic Jet in the Narrow-Line Quasar PMN J0948+0022. Astrophys. J. 2009, 699, 976–984. [Google Scholar] [CrossRef] [Green Version]
  12. Foschini, L. Evidence of powerful relativistic jets in narrow-line Seyfert 1 galaxies. arXiv 2011, arXiv:1105.0772. [Google Scholar]
  13. Foschini, L.; Berton, M.; Caccianiga, A.; Ciroi, S.; Cracco, V.; Peterson, B.M.; Angelakis, E.; Braito, V.; Fuhrmann, L.; Gallo, L.; et al. Properties of flat-spectrum radio-loud narrow-line Seyfert 1 galaxies. Astron. Astrophys. 2015, 575, A13. [Google Scholar] [CrossRef]
  14. Krongold, Y.; Dultzin-Hacyan, D.; Marziani, P. Host Galaxies and Circumgalactic Environment of “Narrow Line” Seyfert 1 Nuclei. Astrophys. J. 2001, 121, 702–709. [Google Scholar] [CrossRef] [Green Version]
  15. Crenshaw, D.M.; Kraemer, S.B.; Gabel, J.R. The Host Galaxies of Narrow-Line Seyfert 1 Galaxies: Evidence for Bar-Driven Fueling. Astrophys. J. 2003, 126, 1690–1698. [Google Scholar] [CrossRef]
  16. Deo, R.P.; Crenshaw, D.M.; Kraemer, S.B. The Host Galaxies of Narrow-Line Seyfert 1 Galaxies: Nuclear Dust Morphology and Starburst Rings. Astrophys. J. 2006, 132, 321–346. [Google Scholar] [CrossRef]
  17. Ohta, K.; Aoki, K.; Kawaguchi, T.; Kiuchi, G. A Bar Fuels a Supermassive Black Hole?: Host Galaxies of Narrow-Line Seyfert 1 Galaxies. Astrophys. J. Suppl. Ser. 2007, 169, 1–20. [Google Scholar] [CrossRef] [Green Version]
  18. Antón, S.; Browne, I.W.A.; Marchã, M.J. The colour of the narrow line Sy1-blazar 0324+3410. Astron. Astrophys. 2008, 490, 583–587. [Google Scholar] [CrossRef] [Green Version]
  19. Orban de Xivry, G.; Davies, R.; Schartmann, M.; Komossa, S.; Marconi, A.; Hicks, E.; Engel, H.; Tacconi, L. The role of secular evolution in the black hole growth of narrow-line Seyfert 1 galaxies. Mon. Not. R. Astron. Soc. 2011, 417, 2721–2736. [Google Scholar] [CrossRef] [Green Version]
  20. Mathur, S.; Fields, D.; Peterson, B.M.; Grupe, D. Supermassive Black Holes, Pseudobulges, and the Narrow-line Seyfert 1 Galaxies. Astrophys. J. 2012, 754, 146. [Google Scholar] [CrossRef] [Green Version]
  21. Kotilainen, J.K.; León-Tavares, J.; Olguín-Iglesias, A.; Baes, M.; Anórve, C.; Chavushyan, V.; Carrasco, L. Discovery of a Pseudobulge Galaxy Launching Powerful Relativistic Jets. Astrophys. J. 2016, 832, 157. [Google Scholar] [CrossRef] [Green Version]
  22. Olguín-Iglesias, A.; Kotilainen, J.K.; León Tavares, J.; Chavushyan, V.; Añorve, C. Evidence of bar-driven secular evolution in the gamma-ray narrow-line Seyfert 1 galaxy FBQS J164442.5+261913. Mon. Not. R. Astron. Soc. 2017, 467, 3712–3722. [Google Scholar] [CrossRef]
  23. Järvelä, E.; Lähteenmäki, A.; Berton, M. Near-infrared morphologies of the host galaxies of narrow-line Seyfert 1 galaxies. Astron. Astrophys. 2018, 619, A69. [Google Scholar] [CrossRef] [Green Version]
  24. Berton, M.; Congiu, E.; Ciroi, S.; Komossa, S.; Frezzato, M.; Di Mille, F.; Antón, S.; Antonucci, R.; Caccianiga, A.; Coppi, P.; et al. The Interacting Late-type Host Galaxy of the Radio-loud Narrow-line Seyfert 1 IRAS 20181-2244. Astrophys. J. 2019, 157, 48. [Google Scholar] [CrossRef] [Green Version]
  25. Olguín-Iglesias, A.; Kotilainen, J.; Chavushyan, V. The disc-like host galaxies of radio-loud Narrow-Line Seyfert 1s. arXiv 2019, arXiv:1912.02746. [Google Scholar] [CrossRef] [Green Version]
  26. Grupe, D. Statistical properties of narrow-line Seyfert 1 galaxies. New Astron. Rev. 2000, 44, 455–460. [Google Scholar] [CrossRef] [Green Version]
  27. Sulentic, J.W.; Zwitter, T.; Marziani, P.; Dultzin-Hacyan, D. Eigenvector 1: An Optimal Correlation Space for Active Galactic Nuclei. Astrophys. J. 2000, 536, L5–L9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Mathur, S.; Kuraszkiewicz, J.; Czerny, B. Evolution of active galaxies: Black-hole mass-bulge relations for narrow line objects. New Astron. 2001, 6, 321–329. [Google Scholar] [CrossRef] [Green Version]
  29. Peterson, B.M. Masses of Black Holes in Active Galactic Nuclei: Implications for Narrow-Line Seyfert 1 Galaxies. In Proceedings of the Science, Volume Narrow-Line Seyfert 1 Galaxies and Their Place in the Universe, Milan, Italy, 4–6 April 2011; p. 32. [Google Scholar]
  30. Berton, M.; Caccianiga, A.; Foschini, L.; Peterson, B.M.; Mathur, S.; Terreran, G.; Ciroi, S.; Congiu, E.; Cracco, V.; Frezzato, M.; et al. Compact steep-spectrum sources as the parent population of flat-spectrum radio-loud narrow-line Seyfert 1 galaxies. Astron. Astrophys. 2016, 591, A98. [Google Scholar] [CrossRef] [Green Version]
  31. Paliya, V.S. Gamma-ray emitting narrow-line Seyfert 1 galaxies: Past, present, and future. J. Astrophys. Astron. 2019, 40, 39. [Google Scholar] [CrossRef] [Green Version]
  32. Berton, M.; Björklund, I.; Lähteenmäki, A.; Congiu, E.; Järvelä, E.; Terreran, G.; La Mura, G. Line shapes in narrow-line Seyfert 1 galaxies: A tracer of physical properties? Contrib. Astron. Obs. Skaln. Pleso 2020, 50, 270–292. [Google Scholar] [CrossRef]
  33. Foschini, L. Jetted Narrow-Line Seyfert 1 Galaxies & Co.: Where Do We Stand? Universe 2020, 6, 136. [Google Scholar] [CrossRef]
  34. Osterbrock, D.E.; Pogge, R.W. The spectra of narrow-line Seyfert 1 galaxies. Astrophys. J. 1985, 297, 166–176. [Google Scholar] [CrossRef]
  35. Rakshit, S.; Stalin, C.S.; Chand, H.; Zhang, X.G. A Catalog of Narrow Line Seyfert 1 Galaxies from the Sloan Digital Sky Survey Data Release 12. Astrophys. J. Suppl. Ser. 2017, 229, 39. [Google Scholar] [CrossRef]
  36. Sulentic, J.W.; Marziani, P.; Zamanov, R.; Bachev, R.; Calvani, M.; Dultzin-Hacyan, D. Average Quasar Spectra in the Context of Eigenvector 1. Astrophys. J. 2002, 566, L71–L75. [Google Scholar] [CrossRef]
  37. Marziani, P.; Dultzin, D.; Sulentic, J.W.; Del Olmo, A.; Negrete, C.A.; Martínez-Aldama, M.L.; D’Onofrio, M.; Bon, E.; Bon, N.; Stirpe, G.M. A main sequence for quasars. Front. Astron. Space Sci. 2018, 5, 6. [Google Scholar] [CrossRef] [Green Version]
  38. Komossa, S.; Xu, D.; Zhou, H.; Storchi-Bergmann, T.; Binette, L. On the Nature of Seyfert Galaxies with High [O III] λ5007 Blueshifts. Astrophys. J. 2008, 680, 926–938. [Google Scholar] [CrossRef] [Green Version]
  39. Zamanov, R.; Marziani, P.; Sulentic, J.W.; Calvani, M.; Dultzin-Hacyan, D.; Bachev, R. Kinematic Linkage between the Broad- and Narrow-Line-emitting Gas in Active Galactic Nuclei. Astrophys. J. 2002, 576, L9–L13. [Google Scholar] [CrossRef]
  40. Marziani, P.; Zamanov, R.K.; Sulentic, J.W.; Calvani, M. Searching for the physical drivers of eigenvector 1: Influence of black hole mass and Eddington ratio. Mon. Not. R. Astron. Soc. 2003, 345, 1133–1144. [Google Scholar] [CrossRef] [Green Version]
  41. Boroson, T.A.; Green, R.F. The emission-line properties of low-redshift quasi-stellar objects. Astrophys. J. Suppl. Ser. 1992, 80, 109–135. [Google Scholar] [CrossRef]
  42. Shen, Y.; Ho, L.C. The diversity of quasars unified by accretion and orientation. Nature 2014, 513, 210–213. [Google Scholar] [CrossRef] [Green Version]
  43. Panda, S.; Marziani, P.; Czerny, B. The Quasar Main Sequence Explained by the Combination of Eddington Ratio, Metallicity, and Orientation. Astrophys. J. 2019, 882, 79. [Google Scholar] [CrossRef] [Green Version]
  44. Marziani, P.; Sulentic, J.W.; Zwitter, T.; Dultzin-Hacyan, D.; Calvani, M. Searching for the Physical Drivers of the Eigenvector 1 Correlation Space. Astrophys. J. 2001, 558, 553–560. [Google Scholar] [CrossRef]
  45. Marziani, P.; del Olmo, A.; D’Onofrio, M.; Dultzin, D.; Negrete, C.A.; Martínez-Aldama, M.L.; Bon, E.; Bon, N.; Stirpe, G.M. Narrow-line Seyfert 1s: What is wrong in a name? arXiv 2018, arXiv:1807.03003. [Google Scholar]
  46. Czerny, B.; Panda, S.; Śniegowska, M.; Kozlowski, S.; Nikolajuk, M.; Du, P.; You, B. Narrow Line Seyfert 1 galaxies in the context of Quasar Main Sequence. arXiv 2018, arXiv:1806.06741. [Google Scholar]
  47. Marziani, P.; Martínez Carballo, M.A.; Sulentic, J.W.; Del Olmo, A.; Stirpe, G.M.; Dultzin, D. The most powerful quasar outflows as revealed by the Civ λ1549 resonance line. Astrophys. Space Sci. 2016, 361, 29. [Google Scholar] [CrossRef] [Green Version]
  48. Komossa, S.; Xu, D.W.; Wagner, A.Y. Extreme gaseous outflows in radio-loud narrow-line Seyfert 1 galaxies. Mon. Not. R. Astron. Soc. 2018, 477, 5115–5126. [Google Scholar] [CrossRef]
  49. Ganci, V.; Marziani, P.; D’Onofrio, M.; del Olmo, A.; Bon, E.; Bon, N.; Negrete, C.A. Radio loudness along the quasar main sequence. Astron. Astrophys. 2019, 630, A110. [Google Scholar] [CrossRef]
  50. Boroson, T.A. Black Hole Mass and Eddington Ratio as Drivers for the Observable Properties of Radio-loud and Radio-quiet QSOs. Astrophys. J. 2002, 565, 78–85. [Google Scholar] [CrossRef]
  51. Grupe, D.; Mathur, S.; Komossa, S. Markarian 1239: A Highly Polarized Narrow-Line Seyfert 1 Galaxy with a Steep X-Ray Spectrum and Strong Ne IX Emission. Astrophys. J. 2004, 127, 3161–3167. [Google Scholar] [CrossRef]
  52. Grupe, D.; Komossa, S.; Leighly, K.M.; Page, K.L. The Simultaneous Optical-to-X-Ray Spectral Energy Distribution of Soft X-Ray Selected Active Galactic Nuclei Observed by Swift. Astrophys. J. Suppl. Ser. 2010, 187, 64–106. [Google Scholar] [CrossRef] [Green Version]
  53. Xu, D.; Komossa, S.; Zhou, H.; Lu, H.; Li, C.; Grupe, D.; Wang, J.; Yuan, W. Correlation Analysis of a Large Sample of Narrow-line Seyfert 1 Galaxies: Linking Central Engine and Host Properties. Astrophys. J. 2012, 143, 83. [Google Scholar] [CrossRef] [Green Version]
  54. Berton, M.; Foschini, L.; Ciroi, S.; Cracco, V.; La Mura, G.; Di Mille, F.; Rafanelli, P. [O III] line properties in two samples of radio-emitting narrow-line Seyfert 1 galaxies. Astron. Astrophys. 2016, 591, A88. [Google Scholar] [CrossRef]
  55. Komatsu, E.; Smith, K.M.; Dunkley, J.; Bennett, C.L.; Gold, B.; Hinshaw, G.; Jarosik, N.; Larson, D.; Nolta, M.R.; Page, L.; et al. Seven-year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpretation. Astrophys. J. Suppl. Ser. 2011, 192, 18. [Google Scholar] [CrossRef] [Green Version]
  56. Berton, M.; Congiu, E.; Järvelä, E.; Antonucci, R.; Kharb, P.; Lister, M.L.; Tarchi, A.; Caccianiga, A.; Chen, S.; Foschini, L.; et al. Radio-emitting narrow-line Seyfert 1 galaxies in the JVLA perspective. Astron. Astrophys. 2018, 614, A87. [Google Scholar] [CrossRef]
  57. Lähteenmäki, A.; Järvelä, E.; Ramakrishnan, V.; Tornikoski, M.; Tammi, J.; Vera, R.J.C.; Chamani, W. Radio jets and gamma-ray emission in radio-silent narrow-line Seyfert 1 galaxies. Astron. Astrophys. 2018, 614, L1. [Google Scholar] [CrossRef] [Green Version]
  58. Berton, M.; Järvelä, E.; Crepaldi, L.; Lähteenmäki, A.; Tornikoski, M.; Congiu, E.; Kharb, P.; Terreran, G.; Vietri, A. Absorbed relativistic jets in radio-quiet narrow-line Seyfert 1 galaxies. Astron. Astrophys. 2020, 636, A64. [Google Scholar] [CrossRef]
  59. Blanton, M.R.; Bershady, M.A.; Abolfathi, B.; Albareti, F.D.; Allende Prieto, C.; Almeida, A.; Alonso-García, J.; Anders, F.; Anderson, S.F.; Andrews, B.; et al. Sloan Digital Sky Survey IV: Mapping the Milky Way, Nearby Galaxies, and the Distant Universe. Astrophys. J. 2017, 154, 28. [Google Scholar] [CrossRef]
  60. Tripathi, S.; McGrath, K.M.; Gallo, L.C.; Grupe, D.; Komossa, S.; Berton, M.; Kriss, G.; Longinotti, A.L. Tracking the year-to-year variation in the spectral energy distribution of the narrow-line Seyfert 1 galaxy Mrk 335. Mon. Not. R. Astron. Soc. 2020, 499, 1266–1286. [Google Scholar] [CrossRef]
  61. Sulentic, J.W.; Bachev, R.; Marziani, P.; Negrete, C.A.; Dultzin, D. C IV λ1549 as an Eigenvector 1 Parameter for Active Galactic Nuclei. Astrophys. J. 2007, 666, 757–777. [Google Scholar] [CrossRef] [Green Version]
  62. Abazajian, K.N.; Adelman-McCarthy, J.K.; Agüeros, M.A.; Allam, S.S.; Allende Prieto, C.; An, D.; Anderson, K.S.J.; Anderson, S.F.; Annis, J.; Bahcall, N.A.; et al. The Seventh Data Release of the Sloan Digital Sky Survey. Astrophys. J. Suppl. Ser. 2009, 182, 543–558. [Google Scholar] [CrossRef]
  63. Kovačević, J.; Popović, L.Č.; Dimitrijević, M.S. Analysis of Optical Fe II Emission in a Sample of Active Galactic Nucleus Spectra. Astrophys. J. Suppl. Ser. 2010, 189, 15–36. [Google Scholar] [CrossRef] [Green Version]
  64. Shapovalova, A.I.; Popović, L.Č.; Burenkov, A.N.; Chavushyan, V.H.; Ilić, D.; Kovačević, A.; Kollatschny, W.; Kovačević, J.; Bochkarev, N.G.; Valdes, J.R.; et al. Spectral Optical Monitoring of the Narrow-line Seyfert 1 Galaxy Ark 564. Astrophys. J. Suppl. Ser. 2012, 202, 10. [Google Scholar] [CrossRef] [Green Version]
  65. Dimitrijević, M.S.; Popović, L.Č.; Kovačević, J.; Dačić, M.; Ilić, D. The flux ratio of the [OIII] λλ5007, 4959 lines in AGN: Comparison with theoretical calculations. Mon. Not. R. Astron. Soc. 2007, 374, 1181–1184. [Google Scholar] [CrossRef] [Green Version]
  66. Berton, M.; Foschini, L.; Ciroi, S.; Cracco, V.; La Mura, G.; Lister, M.L.; Mathur, S.; Peterson, B.M.; Richards, J.L.; Rafanelli, P. Parent population of flat-spectrum radio-loud narrow-Line Seyfert 1 galaxies. Astron. Astrophys. 2015, 578, A28. [Google Scholar] [CrossRef] [Green Version]
  67. Kellermann, K.I.; Sramek, R.; Schmidt, M.; Shaffer, D.B.; Green, R. VLA observations of objects in the Palomar Bright Quasar Survey. Astrophys. J. 1989, 98, 1195–1207. [Google Scholar] [CrossRef]
  68. Marziani, P.; Sulentic, J.W.; Stirpe, G.M.; Dultzin, D.; Del Olmo, A.; Martínez-Carballo, M.A. Blue outliers among intermediate redshift quasars. Astrophys. Space Sci. 2016, 361, 3. [Google Scholar] [CrossRef] [Green Version]
  69. Padovani, P. On the two main classes of active galactic nuclei. Nat. Astron. 2017, 1, 0194. [Google Scholar] [CrossRef] [Green Version]
  70. Järvelä, E.; Lähteenmäki, A.; Lietzen, H.; Poudel, A.; Heinämäki, P.; Einasto, M. Large-scale environments of narrow-line Seyfert 1 galaxies. Astron. Astrophys. 2017, 606, A9. [Google Scholar] [CrossRef] [Green Version]
  71. Panessa, F.; Baldi, R.D.; Laor, A.; Padovani, P.; Behar, E.; McHardy, I. The origin of radio emission from radio-quiet active galactic nuclei. Nat. Astron. 2019, 3, 387–396. [Google Scholar] [CrossRef]
  72. Caccianiga, A.; Antón, S.; Ballo, L.; Foschini, L.; Maccacaro, T.; Della Ceca, R.; Severgnini, P.; Marchã, M.J.; Mateos, S.; Sani, E. WISE colours and star formation in the host galaxies of radio-loud narrow-line Seyfert 1. Mon. Not. R. Astron. Soc. 2015, 451, 1795–1805. [Google Scholar] [CrossRef]
  73. Chen, S.; Järvelä, E.; Crepaldi, L.; Zhou, M.; Ciroi, S.; Berton, M.; Kharb, P.; Foschini, L.; Gu, M.; La Mura, G.; et al. Radio morphology of southern narrow-line Seyfert 1 galaxies with Very Large Array observations. Mon. Not. R. Astron. Soc. 2020, 498, 1278–1297. [Google Scholar] [CrossRef]
  74. Falcke, H.; Biermann, P.L. The jet-disk symbiosis. I. Radio to X-ray emission models for quasars. Astron. Astrophys. 1995, 293, 665–682. [Google Scholar]
  75. Blandford, R.D.; Königl, A. Relativistic jets as compact radio sources. Astrophys. J. 1979, 232, 34–48. [Google Scholar] [CrossRef]
  76. Hou, A.; Parker, L.C.; Harris, W.E.; Wilman, D.J. Statistical Tools for Classifying Galaxy Group Dynamics. Astrophys. J. 2009, 702, 1199–1210. [Google Scholar] [CrossRef] [Green Version]
  77. de Bruyn, A.G.; Wilson, A.S. The radio properties of Seyfert galaxies. Astron. Astrophys. 1978, 64, 433–444. [Google Scholar]
  78. Sulentic, J.W.; del Olmo, A.; Marziani, P.; Martínez-Carballo, M.A.; D’Onofrio, M.; Dultzin, D.; Perea, J.; Martínez-Aldama, M.L.; Negrete, C.A.; Stirpe, G.M.; et al. What does CIVλ1549 tell us about the physical driver of the Eigenvector quasar sequence? Astron. Astrophys. 2017, 608, A122. [Google Scholar] [CrossRef] [Green Version]
  79. Mullaney, J.R.; Alexander, D.M.; Fine, S.; Goulding, A.D.; Harrison, C.M.; Hickox, R.C. Narrow-line region gas kinematics of 24 264 optically selected AGN: The radio connection. Mon. Not. R. Astron. Soc. 2013, 433, 622–638. [Google Scholar] [CrossRef] [Green Version]
  80. Cracco, V.; Ciroi, S.; Berton, M.; Di Mille, F.; Foschini, L.; La Mura, G.; Rafanelli, P. A spectroscopic analysis of a sample of narrow-line Seyfert 1 galaxies selected from the Sloan Digital Sky Survey. Mon. Not. R. Astron. Soc. 2016, 462, 1256–1280. [Google Scholar] [CrossRef] [Green Version]
  81. Vaona, L.; Ciroi, S.; Di Mille, F.; Cracco, V.; La Mura, G.; Rafanelli, P. Spectral properties of the narrow-line region in Seyfert galaxies selected from the SDSS-DR7. Mon. Not. R. Astron. Soc. 2012, 427, 1266–1283. [Google Scholar] [CrossRef] [Green Version]
  82. Concas, A.; Popesso, P.; Brusa, M.; Mainieri, V.; Erfanianfar, G.; Morselli, L. Light breeze in the local Universe. Astron. Astrophys. 2017, 606, A36. [Google Scholar] [CrossRef] [Green Version]
  83. Bian, W.; Yuan, Q.; Zhao, Y. The blueshift of the [OIII] emission line in narrow-line Seyfert 1 galaxies. Mon. Not. R. Astron. Soc. 2005, 364, 187–194. [Google Scholar] [CrossRef] [Green Version]
  84. Negrete, C.A.; Dultzin, D.; Marziani, P.; Esparza, D.; Sulentic, J.W.; del Olmo, A.; Martínez-Aldama, M.L.; García López, A.; D’Onofrio, M.; Bon, N.; et al. Highly accreting quasars: The SDSS low-redshift catalog. Astron. Astrophys. 2018, 620, A118. [Google Scholar] [CrossRef] [Green Version]
  85. Aoki, K.; Kawaguchi, T.; Ohta, K. The Largest Blueshifts of the [O III] Emission Line in Two Narrow-Line Quasars. Astrophys. J. 2005, 618, 601–608. [Google Scholar] [CrossRef] [Green Version]
  86. Venturi, G.; Cresci, G.; Marconi, A.; Mingozzi, M.; Nardini, E.; Carniani, S.; Mannucci, F.; Marasco, A.; Maiolino, R.; Perna, M.; et al. MAGNUM survey: Compact jets causing large turmoil in galaxies— Enhanced line widths perpendicular to radio jets as tracers of jet-ISM interaction. arXiv 2020, arXiv:2011.04677. [Google Scholar]
  87. Gu, M.; Chen, Y.; Komossa, S.; Yuan, W.; Shen, Z.; Wajima, K.; Zhou, H.; Zensus, J.A. The Radio Properties of Radio-loud Narrow-line Seyfert 1 Galaxies on Parsec Scales. Astrophys. J. Suppl. Ser. 2015, 221, 3. [Google Scholar] [CrossRef] [Green Version]
  88. Berton, M. An updated view on the parent population of γ-ray emitting narrow-line Seyfert 1 galaxies. In Proceedings of the Science, Volume Revisiting narrow-line Seyfert 1 galaxies and their place in the Universe, Padova Botanical Garden, Italy, 9–13 April 2018; p. 26. [Google Scholar]
  89. Liao, M.; Gu, M. Investigation on young radio AGNs based on SDSS spectroscopy. Mon. Not. R. Astron. Soc. 2020, 491, 92–112. [Google Scholar] [CrossRef] [Green Version]
  90. Gliozzi, M.; Papadakis, I.E.; Grupe, D.; Brinkmann, W.P.; Raeth, C.; Kedziora-Chudczer, L. A Panchromatic View of PKS 0558-504: An Ideal Laboratory to Study the Disk-Jet Link. Astrophys. J. 2010, 717, 1243–1252. [Google Scholar] [CrossRef]
  91. Richards, J.L.; Lister, M.L. Kiloparsec-Scale Jets in Three Radio-Loud Narrow-Line Seyfert 1 Galaxies. Astrophys. J. 2015, 800, L8. [Google Scholar] [CrossRef] [Green Version]
  92. Congiu, E.; Berton, M.; Giroletti, M.; Antonucci, R.; Caccianiga, A.; Kharb, P.; Lister, M.L.; Foschini, L.; Ciroi, S.; Cracco, V.; et al. Kiloparsec-scale emission in the narrow-line Seyfert 1 galaxy Mrk 783. Astron. Astrophys. 2017, 603, A32. [Google Scholar] [CrossRef] [Green Version]
  93. Rakshit, S.; Stalin, C.S.; Hota, A.; Konar, C. Rare Finding of a 100 Kpc Large, Double-lobed Radio Galaxy Hosted in the Narrow-line Seyfert 1 Galaxy SDSS J103024.95+551622.7. Astrophys. J. 2018, 869, 173. [Google Scholar] [CrossRef] [Green Version]
  94. Gabányi, K.É.; Frey, S.; Paragi, Z.; Järvelä, E.; Morokuma, T.; An, T.; Tanaka, M.; Tar, I. The radio structure of the peculiar narrow-line Seyfert 1 galaxy candidate J1100+4421. Mon. Not. R. Astron. Soc. 2018, 473, 1554–1561. [Google Scholar] [CrossRef] [Green Version]
  95. Congiu, E.; Kharb, P.; Tarchi, A.; Berton, M.; Caccianiga, A.; Chen, S.; Crepaldi, L.; Di Mille, F.; Järvelä, E.; Jarvis, M.E.; et al. The radio structure of the narrow-line Seyfert 1 Mrk 783 with VLBA and e-MERLIN. Mon. Not. R. Astron. Soc. 2020, 499, 3149–3157. [Google Scholar] [CrossRef]
  96. Oshlack, A.Y.K.N.; Webster, R.L.; Whiting, M.T. A Very Radio Loud Narrow-Line Seyfert 1: PKS 2004-447. Astrophys. J. 2001, 558, 578–582. [Google Scholar] [CrossRef] [Green Version]
  97. Gallo, L.C. Investigating the nature of narrow-line Seyfert 1 galaxies with high-energy spectral complexity. Mon. Not. R. Astron. Soc. 2006, 368, 479–486. [Google Scholar] [CrossRef] [Green Version]
  98. Caccianiga, A.; Antón, S.; Ballo, L.; Dallacasa, D.; Della Ceca, R.; Fanali, R.; Foschini, L.; Hamilton, T.; Kraus, A.; Maccacaro, T.; et al. SDSS J143244.91+301435.3: A link between radio-loud narrow-line Seyfert 1 galaxies and compact steep-spectrum radio sources? Mon. Not. R. Astron. Soc. 2014, 441, 172–186. [Google Scholar] [CrossRef] [Green Version]
  99. Berton, M.; Foschini, L.; Caccianiga, A.; Ciroi, S.; Congiu, E.; Cracco, V.; Frezzato, M.; La Mura, G.; Rafanelli, P. An orientation-based unification of young jetted active galactic nuclei: The case of 3C 286. Front. Astron. Space Sci. 2017, 4, 8. [Google Scholar] [CrossRef] [Green Version]
  100. Caccianiga, A.; Dallacasa, D.; Antón, S.; Ballo, L.; Berton, M.; Mack, K.H.; Paulino-Afonso, A. SDSSJ143244.91+301435.3 at VLBI: A compact radio galaxy in a narrow-line Seyfert 1. Mon. Not. R. Astron. Soc. 2017, 464, 1474–1480. [Google Scholar] [CrossRef] [Green Version]
  101. O’Dea, C.P.; Saikia, D.J. Compact steep-spectrum and peaked-spectrum radio sources. arXiv 2020, arXiv:2009.02750. [Google Scholar] [CrossRef]
  102. Yao, S.; Komossa, S. Spectroscopic classification, variability and SED of the Fermi-detected CSS 3C 286: The radio-loudest NLS1 galaxy? arXiv 2020, arXiv:2012.01785. [Google Scholar]
  103. Tadhunter, C.; Wills, K.; Morganti, R.; Oosterloo, T.; Dickson, R. Emission-line outflows in PKS1549-79: The effects of the early stages of radio-source evolution? Mon. Not. R. Astron. Soc. 2001, 327, 227–232. [Google Scholar] [CrossRef] [Green Version]
  104. Holt, J.; Tadhunter, C.; Morganti, R.; Bellamy, M.; González Delgado, R.M.; Tzioumis, A.; Inskip, K.J. The co-evolution of the obscured quasar PKS 1549-79 and its host galaxy: Evidence for a high accretion rate and warm outflow. Mon. Not. R. Astron. Soc. 2006, 370, 1633–1650. [Google Scholar] [CrossRef]
  105. Morganti, R.; Oosterloo, T.; Oonk, J.B.R.; Frieswijk, W.; Tadhunter, C. The fast molecular outflow in the Seyfert galaxy IC 5063 as seen by ALMA. Astron. Astrophys. 2015, 580, A1. [Google Scholar] [CrossRef]
  106. Orienti, M. Radio properties of Compact Steep Spectrum and GHz-Peaked Spectrum radio sources. Astron. Nachrichten 2016, 337, 9. [Google Scholar] [CrossRef] [Green Version]
  107. Becker, R.H.; White, R.L.; Helfand, D.J. The FIRST Survey: Faint Images of the Radio Sky at Twenty Centimeters. Astrophys. J. 1995, 450, 559. [Google Scholar] [CrossRef]
  108. Condon, J.J.; Cotton, W.D.; Greisen, E.W.; Yin, Q.F.; Perley, R.A.; Taylor, G.B.; Broderick, J.J. The NRAO VLA Sky Survey. Astrophys. J. 1998, 115, 1693–1716. [Google Scholar] [CrossRef]
  109. Spoon, H.W.W.; Holt, J. Discovery of Strongly Blueshifted Mid-Infrared [Ne III] and [Ne V] Emission in ULIRGs. Astrophys. J. 2009, 702, L42–L46. [Google Scholar] [CrossRef] [Green Version]
  110. Gaskell, C.M.; Goosmann, R.W. Line Shifts, Broad-line Region Inflow, and the Feeding of Active Galactic Nuclei. Astrophys. J. 2013, 769, 30. [Google Scholar] [CrossRef] [Green Version]
  111. Cecil, G.; Dopita, M.A.; Groves, B.; Wilson, A.S.; Ferruit, P.; Pécontal, E.; Binette, L. Spatial Resolution of High-Velocity Filaments in the Narrow-Line Region of NGC 1068: Associated Absorbers Caught in Emission? Astrophys. J. 2002, 568, 627–638. [Google Scholar] [CrossRef] [Green Version]
  112. Ludwig, R.R.; Greene, J.E.; Barth, A.J.; Ho, L.C. Physical Properties of the Narrow-line Region of Low-mass Active Galaxies. Astrophys. J. 2012, 756, 51. [Google Scholar] [CrossRef] [Green Version]
  113. Wagner, A.Y.; Bicknell, G.V.; Umemura, M. Driving Outflows with Relativistic Jets and the Dependence of Active Galactic Nucleus Feedback Efficiency on Interstellar Medium Inhomogeneity. Astrophys. J. 2012, 757, 136. [Google Scholar] [CrossRef]
  114. Ferrarese, L.; Merritt, D. A Fundamental Relation between Supermassive Black Holes and Their Host Galaxies. Astrophys. J. 2000, 539, L9–L12. [Google Scholar] [CrossRef] [Green Version]
  115. Marin, F. Are there reliable methods to estimate the nuclear orientation of Seyfert galaxies? Mon. Not. R. Astron. Soc. 2016, 460, 3679–3705. [Google Scholar] [CrossRef] [Green Version]
  116. Boroson, T.A. A New Orientation Indicator for Radio-quiet Quasars. Astrophys. J. 2011, 735, L14. [Google Scholar] [CrossRef]
  117. Hu, C.; Wang, J.M.; Ho, L.C.; Chen, Y.M.; Zhang, H.T.; Bian, W.H.; Xue, S.J. A Systematic Analysis of Fe II Emission in Quasars: Evidence for Inflow to the Central Black Hole. Astrophys. J. 2008, 687, 78–96. [Google Scholar] [CrossRef] [Green Version]
  118. Bon, E.; Marziani, P.; Berton, M.; Bon, N.; Antonucci, R.; Gaskell, M.; Ferland, G. Fe II velocity shifts in optical spectra of type 1 AGN. arXiv 2018, arXiv:1812.05828. [Google Scholar]
  119. Wagner, A.Y.; Bicknell, G.V. Relativistic Jet Feedback in Evolving Galaxies. Astrophys. J. 2011, 728, 29. [Google Scholar] [CrossRef]
  120. Cielo, S.; Bieri, R.; Volonteri, M.; Wagner, A.Y.; Dubois, Y. AGN feedback compared: Jets versus radiation. Mon. Not. R. Astron. Soc. 2018, 477, 1336–1355. [Google Scholar] [CrossRef] [Green Version]
Figure 1. [O III] line fitting of the source J1114+3241. The black solid line represents the restframe spectrum after continuum and Fe II subtraction, the green dashed line represents the [O III] core components, the orange dashed line represents the wing components, and the best fit is shown as the red dashed line. The blue solid line shows the residuals, and the black horizontal dashed line corresponds to a zero-level flux.
Figure 1. [O III] line fitting of the source J1114+3241. The black solid line represents the restframe spectrum after continuum and Fe II subtraction, the green dashed line represents the [O III] core components, the orange dashed line represents the wing components, and the best fit is shown as the red dashed line. The blue solid line shows the residuals, and the black horizontal dashed line corresponds to a zero-level flux.
Universe 07 00188 g001
Figure 2. Left: [O III] core velocity shift in km s 1 as a function of the bolometric luminosity in ergs s 1 . The horizontal dashed lines indicate the limits for red or blue outliers at 150 km s 1 . Right: FWHM([O III]) core component in km s 1 as a function of the bolometric luminosity in ergs s 1 .
Figure 2. Left: [O III] core velocity shift in km s 1 as a function of the bolometric luminosity in ergs s 1 . The horizontal dashed lines indicate the limits for red or blue outliers at 150 km s 1 . Right: FWHM([O III]) core component in km s 1 as a function of the bolometric luminosity in ergs s 1 .
Universe 07 00188 g002
Figure 3. Left: distributions of the [O III] core velocity shift in km s 1 , in non-jetted (top panel) and jetted sources (bottom panel). The vertical dashed lines indicate the limits for red or blue outliers at 150 km s 1 . Binning of 25 km s 1 . Right: distributions of the FWHM([O III]) core component in km s 1 , in non-jetted (top panel) and jetted sources (bottom panel). Binning of 25 km s 1 .
Figure 3. Left: distributions of the [O III] core velocity shift in km s 1 , in non-jetted (top panel) and jetted sources (bottom panel). The vertical dashed lines indicate the limits for red or blue outliers at 150 km s 1 . Binning of 25 km s 1 . Right: distributions of the FWHM([O III]) core component in km s 1 , in non-jetted (top panel) and jetted sources (bottom panel). Binning of 25 km s 1 .
Universe 07 00188 g003
Figure 4. Correlation between the FWHM of the [O III] core component and its shift. Non-jetted sources are represented with the red squares, while jetted sources are indicated with the blue circles.
Figure 4. Correlation between the FWHM of the [O III] core component and its shift. Non-jetted sources are represented with the red squares, while jetted sources are indicated with the blue circles.
Universe 07 00188 g004
Figure 5. Left: distributions of the [O III] wing velocity shift in km s 1 , in non-jetted (top panel) and jetted sources (bottom panel). Binning of 25 km s 1 . Right: distributions of the FWHM([O III]) wing component in km s 1 , in non-jetted (top panel) and jetted sources (bottom panel). Binning of 50 km s 1 .
Figure 5. Left: distributions of the [O III] wing velocity shift in km s 1 , in non-jetted (top panel) and jetted sources (bottom panel). Binning of 25 km s 1 . Right: distributions of the FWHM([O III]) wing component in km s 1 , in non-jetted (top panel) and jetted sources (bottom panel). Binning of 50 km s 1 .
Universe 07 00188 g005
Figure 6. Left: distributions of the [O III] core velocity shift in km s 1 , in compact (top panel), intermediate (middle panel) and extended sources (bottom panel). The vertical dashed line indicate the limit for red or blue outliers at 150 km s 1 . Binning of 25 km s 1 . Right: distributions of the FWHM([O III]) core component in km s 1 in the same samples. Binning of 25 km s 1 .
Figure 6. Left: distributions of the [O III] core velocity shift in km s 1 , in compact (top panel), intermediate (middle panel) and extended sources (bottom panel). The vertical dashed line indicate the limit for red or blue outliers at 150 km s 1 . Binning of 25 km s 1 . Right: distributions of the FWHM([O III]) core component in km s 1 in the same samples. Binning of 25 km s 1 .
Universe 07 00188 g006
Figure 7. Left: Left: Redshift distribution of the jetted (top panel) and the non-jetted (bottom panel) sources. Binning of 0.025. Right: distributions of the logarithm of the [O III] luminosity (in erg s 1 ) for the jetted (top panel) and the non-jetted (bottom panel) sources. Binning of 0.25 dex.
Figure 7. Left: Left: Redshift distribution of the jetted (top panel) and the non-jetted (bottom panel) sources. Binning of 0.025. Right: distributions of the logarithm of the [O III] luminosity (in erg s 1 ) for the jetted (top panel) and the non-jetted (bottom panel) sources. Binning of 0.25 dex.
Universe 07 00188 g007
Figure 8. The bolometric luminosity of our sources as a function of their redshift. Non-jetted objects are shown in red, while jetted sources are plotted in blue.
Figure 8. The bolometric luminosity of our sources as a function of their redshift. Non-jetted objects are shown in red, while jetted sources are plotted in blue.
Universe 07 00188 g008
Figure 9. Left: distributions of the [O III] core velocity shift in km s 1 , in steep-spectrum (top panel) and flat-spectrum jetted sources (bottom panel). The vertical dashed line indicate the limit for red or blue outliers at 150 km s 1 . Binning of 25 km s 1 . Right: distributions of the FWHM([O III]) wing component in km s 1 in the same samples. Binning of 25 km s 1 .
Figure 9. Left: distributions of the [O III] core velocity shift in km s 1 , in steep-spectrum (top panel) and flat-spectrum jetted sources (bottom panel). The vertical dashed line indicate the limit for red or blue outliers at 150 km s 1 . Binning of 25 km s 1 . Right: distributions of the FWHM([O III]) wing component in km s 1 in the same samples. Binning of 25 km s 1 .
Universe 07 00188 g009
Table 1. Observational details of the spectra observed with the Asiago 1.22 m telescope. (1) Short name; (2) exposure time (seconds); (3) resolution; (4) observing date. The asterisk indicates that the spectra were observed on five separate dates, 6 December 2013, 21 September 2014, 12 August 2015, 15 November 2015, 2018 July 2018, and later combined [60]. This operation is possible because the interval between these spectra is shorter than the typical scale of the NLR.
Table 1. Observational details of the spectra observed with the Asiago 1.22 m telescope. (1) Short name; (2) exposure time (seconds); (3) resolution; (4) observing date. The asterisk indicates that the spectra were observed on five separate dates, 6 December 2013, 21 September 2014, 12 August 2015, 15 November 2015, 2018 July 2018, and later combined [60]. This operation is possible because the interval between these spectra is shorter than the typical scale of the NLR.
Short NameExp. TimeRDate
J0006+201225,5007002013–2018 *
J0347+0105480070016 November 2015
J0629−0545600070010 December 2015
J2242+294318,0007006 August 2015
Table 2. [O III] and radio properties of the sample. Columns: (1) short name; (2) logarithm of the [O III] luminosity (erg s 1 ); (3) wavelength of the [O III] core component (Å); (4) velocity of the [O III] core component (km s 1 ); (5) FWHM of the [O III] core component (km s 1 ); (6) wavelength of the [O III] wing component (Å); (7) velocity of the [O III] wing component (km s 1 ); (8) FWHM of the [O III] wing component (km s 1 ); (9) logarithm of the integrated radio luminosity at 5 GHz (erg s 1 ); (10) logarithm of the peak luminosity at 5 GHz (erg s 1 ); (11) logarithm of the diffuse luminosity at 5 GHz (erg s 1 ).
Table 2. [O III] and radio properties of the sample. Columns: (1) short name; (2) logarithm of the [O III] luminosity (erg s 1 ); (3) wavelength of the [O III] core component (Å); (4) velocity of the [O III] core component (km s 1 ); (5) FWHM of the [O III] core component (km s 1 ); (6) wavelength of the [O III] wing component (Å); (7) velocity of the [O III] wing component (km s 1 ); (8) FWHM of the [O III] wing component (km s 1 ); (9) logarithm of the integrated radio luminosity at 5 GHz (erg s 1 ); (10) logarithm of the peak luminosity at 5 GHz (erg s 1 ); (11) logarithm of the diffuse luminosity at 5 GHz (erg s 1 ).
Short Name logL [ O III ] λ c v c FWHM c λ w v w FWHM w logL int logL peak logL diff
J0006+201241.475005.65 ± 0.02−71.63 ± 1.03458.28 ± 4.775001.29 ± 0.11−332.20 ± 6.841387.51 ± 17.0238.4238.4136.88
J0146-004041.165006.67 ± 0.05−10.36 ± 2.99143.83 ± 4.195005.78 ± 0.22−53.34 ± 14.37575.32 ± 19.1637.9537.9436.40
J0347+010541.325004.23 ± 0.12−156.55 ± 7.33328.92 ± 26.555000.64 ± 0.15−371.31 ± 9.111074.27 ± 27.4239.1739.0138.66
J0629-054342.635007.91 ± 0.2563.91 ± 15.15621.47 ± 60.714999.09 ± 0.92−463.97 ± 55.211627.96 ± 52.2040.4340.3639.59
J0632+634040.285007.03 ± 0.2211.20 ± 13.17547.12 ± 4.195002.18 ± 3.92−290.47 ± 242.50999.76 ± 115.5637.7637.5437.39
J0706+390141.335003.61 ± 0.07−193.58 ± 4.19736.07 ± 2.990.000.000.0039.3339.2338.63
J0713+382042.395011.86 ± 0.14300.4 ± 8.38504.81 ± 12.574999.46 ± 0.58−741.81 ± 41.311494.70 ± 17.9639.8239.6839.27
J0744+5149N.D. 40.9840.9739.53
J0752+261741.005006.13 ± 0.06−42.69 ± 3.59227.02 ± 2.405001.64 ± 0.31−268.66 ± 20.36720.56 ± 8.9839.1938.1736.64
J0758+392042.675005.40 ± 0.07−86.40 ± 4.19644.41 ± 2.994993.81 ± 0.40−694.39 ± 26.35468.15 ± 17.9639.6839.3939.36
J0804+385342.195007.01 ± 0.0410.00 ± 2.40157.93 ± 1.205000.52 ± 0.15−388.22 ± 9.58760.04 ± 3.5939.6939.3339.43
J0806+724841.815001.90 ± 0.20−295.97 ± 11.98375.08 ± 8.984994.79 ± 1.49−426.03 ± 99.40733.25 ± 61.6740.1740.1538.93
J0814+560941.734997.51 ± 0.51−558.83 ± 30.54740.30 ± 14.970.000.000.0042.1942.1341.28
J0849+510841.745011.26 ± 0.18264.47 ± 10.78231.26 ± 8.985002.93 ± 1.09−498.77 ± 74.25703.64 ± 50.3043.2143.2041.22
J0850+462642.325008.52 ± 0.16100.41 ± 9.58248.18 ± 14.975005.76 ± 0.52−164.70 ± 38.921558.15 ± 28.7441.7741.7639.98
J0902+044342.205007.75 ± 0.2754.31 ± 16.17286.25 ± 26.354999.03 ± 0.66−521.95 ± 54.491231.01 ± 20.3642.7742.770.00
J0913+365841.435005.90 ± 0.08−56.46 ± 4.79329.96 ± 5.395000.78 ± 0.32−306.62 ± 22.15799.52 ± 7.7838.6638.6137.74
J0925+521741.995007.78 ± 0.5356.1 ± 31.73215.74 ± 1.800.000.000.0038.7238.4238.42
J0926+124441.145006.32 ± 0.61−31.32 ± 36.52337.01 ± 31.735005.00 ± 0.96−79.02 ± 92.21846.06 ± 32.9338.5238.4537.70
J0937+361541.765006.98 ± 0.148.20 ± 8.38404.70 ± 11.985002.80 ± 0.69−249.94 ± 48.50999.76 ± 25.1539.7139.5839.11
J0945+1915N.D. 40.9040.8839.52
J0948+502941.485006.74 ± 0.05−6.17 ± 2.99337.01 ± 1.805001.76 ± 0.17−298.60 ± 11.38868.62 ± 5.3938.9138.8937.52
J0952-013641.675005.96 ± 0.08−52.87 ± 4.79558.40 ± 4.194994.51 ± 0.52−685.82 ± 33.531583.53 ± 13.7739.0438.9938.02
J0957+243341.365005.64 ± 0.06−72.03 ± 3.59327.14 ± 2.405001.11 ± 0.27−271.66 ± 17.96958.86 ± 9.5838.3238.3136.73
J1031+423441.835009.07 ± 0.30133.35 ± 17.96231.26 ± 10.785005.45 ± 0.92−216.94 ± 65.86744.53 ± 10.1841.1441.1339.55
J1034+393841.035006.76 ± 0.04−4.97 ± 2.40228.44 ± 1.805002.54 ± 0.17−253.01 ± 10.78901.05 ± 4.7939.2439.1838.34
J1037+003641.605008.56 ± 0.42102.81 ± 25.15280.61 ± 7.784999.72 ± 0.58−529.11 ± 58.08356.75 ± 22.1542.1942.190.00
J1038+422741.625005.11 ± 0.12−103.77 ± 7.19314.45 ± 7.784999.87 ± 0.52−313.96 ± 36.521146.41 ± 28.1440.6740.1640.51
J1047+472542.895006.66 ± 0.07−10.96 ± 4.19280.61 ± 4.795004.82 ± 0.39−110.43 ± 25.75844.65 ± 24.5543.7743.5143.42
J1048+222241.795003.35 ± 0.31−209.15 ± 18.56389.19 ± 19.764994.49 ± 0.75−530.88 ± 61.071238.06 ± 16.7739.7539.7238.59
J1102+223942.384997.40 ± 0.36−565.41 ± 21.56879.90 ± 11.384984.90 ± 2.41−749.79 ± 164.061319.85 ± 83.2340.4740.4439.35
J1110+365341.735007.53 ± 0.2241.14 ± 13.17400.47 ± 16.770.000.000.0041.8941.8340.98
J1114+324142.134993.71 ± 0.08−786.31 ± 5.06971.78 ± 22.394979.31 ± 0.67−1648.58 ± 40.062101.97 ± 36.3241.0441.0339.35
J1121+535142.015007.09 ± 0.0614.79 ± 3.59305.99 ± 2.995003.61 ± 0.15−208.07 ± 10.78734.66 ± 4.7939.2539.1738.45
J1138+365341.685005.39 ± 0.17−87.00 ± 10.18234.08 ± 16.175004.59 ± 0.43−48.40 ± 34.73923.61 ± 28.7441.0241.0039.77
J1146+323642.005007.10 ± 0.1015.39 ± 5.99307.40 ± 4.195001.47 ± 0.40−336.86 ± 28.74620.44 ± 7.7841.9041.8939.94
J1159+283842.365007.75 ± 0.0454.32 ± 2.17397.77 ± 7.335002.81 ± 0.05−241.67 ± 2.801201.68 ± 5.5539.7639.7438.51
J1203+443139.945006.94 ± 0.455.81 ± 26.94318.68 ± 5.995002.02 ± 0.91−294.61 ± 59.88606.34 ± 81.4336.4635.9336.31
J1209+321741.745006.17 ± 0.19−40.30 ± 11.38410.34 ± 6.595002.46 ± 0.24−221.92 ± 23.951405.86 ± 12.5739.3139.2638.35
J1215+544241.785005.69 ± 0.16−69.04 ± 9.58384.96 ± 7.194997.34 ± 0.58−500.47 ± 42.51661.33 ± 15.5739.3039.2538.34
J1218+294841.265007.70 ± 0.3851.31 ± 22.75341.24 ± 2.405005.81 ± 0.47−113.22 ± 50.90774.14 ± 6.5938.4738.3437.87
J1227+321442.115007.84 ± 0.0559.70 ± 2.99242.54 ± 1.805004.35 ± 0.12−209.05 ± 8.38592.24 ± 3.5939.9439.9038.92
J1238+394241.985005.75 ± 0.15−65.45 ± 8.98380.73 ± 15.570.000.000.0042.1642.160.00
J1242+331741.655006.86 ± 0.041.02 ± 2.40307.40 ± 1.205003.30 ± 0.10−213.29 ± 6.59888.36 ± 2.9938.7438.5138.35
J1246+022241.125006.03 ± 0.07−48.68 ± 4.19310.22 ± 3.595002.24 ± 0.49−227.11 ± 31.73738.89 ± 13.7738.3138.1237.87
J1246+023841.645007.51 ± 0.3539.94 ± 20.96382.14 ± 20.364999.10 ± 1.96−503.49 ± 136.521245.11 ± 70.0641.6941.6740.20
J1302+162442.315003.65 ± 0.05−191.19 ± 2.99408.93 ± 1.200.000.000.0039.9139.2839.80
J1305+511643.145003.62 ± 1.36−192.98 ± 81.431012.45 ± 77.844993.35 ± 2.45−615.00 ± 226.331394.58 ± 43.1142.9142.7142.47
J1317+601041.745007.03 ± 0.0311.31 ± 1.55341.35 ± 6.695000.55 ± 0.34−376.82 ± 20.531435.93 ± 33.8839.4239.2838.86
J1333+414142.345005.35 ± 0.08−89.4 ± 4.79521.73 ± 3.595001.35 ± 0.11−239.37 ± 10.181745.69 ± 5.3939.8539.8438.21
J1337+242342.245005.74 ± 0.25−66.04 ± 14.97352.52 ± 16.174995.76 ± 1.29−597.81 ± 90.412240.64 ± 25.7540.2640.2339.18
J1346+312141.395005.55 ± 0.51−77.42 ± 30.54248.18 ± 55.095001.72 ± 2.60−229.05 ± 184.42920.79 ± 102.3940.6340.630.00
J1355+561242.225007.64 ± 0.0447.72 ± 2.40407.52 ± 1.205003.83 ± 0.10−228.43 ± 7.191013.86 ± 3.5939.6739.6538.15
J1358+265842.365007.86 ± 0.0560.89 ± 2.99448.41 ± 3.595003.64 ± 0.24−252.57 ± 16.171366.38 ± 5.3940.0539.9939.16
J1402+215941.495006.18 ± 0.05−39.70 ± 2.99291.89 ± 1.805003.45 ± 0.13−163.17 ± 8.98762.86 ± 5.3938.3438.2437.67
J1421+2824N.D. 42.2042.1840.66
J1443+472542.324999.40 ± 0.54−445.66 ± 32.331199.99 ± 33.530.000.000.0042.9842.9741.43
J1505+032641.775000.96 ± 0.99−352.25 ± 59.28597.88 ± 60.484995.52 ± 4.34−326.51 ± 317.351180.25 ± 131.1343.0943.090.00
J1536+543341.065007.23 ± 0.7523.17 ± 44.91455.46 ± 16.170.000.000.0037.9837.9237.13
J1537+494241.675006.48 ± 0.08−21.74 ± 4.79153.70 ± 6.595000.77 ± 0.37−341.82 ± 25.15837.60 ± 16.1740.0940.090.00
J1548+351142.405006.63 ± 0.07−12.75 ± 4.19265.10 ± 2.404999.95 ± 0.59−400.03 ± 37.72971.55 ± 27.5437.8537.2337.73
J1555+191140.825005.91 ± 0.43−55.86 ± 25.75349.70 ± 4.790.000.000.0038.1937.7538.00
J1559+350140.585005.48 ± 0.07−81.61 ± 4.19293.30 ± 3.594999.65 ± 0.38−349.15 ± 25.15806.57 ± 14.9742.4642.460.00
J1612+421942.315004.23 ± 0.14−156.46 ± 8.38741.71 ± 12.575005.96 ± 0.36103.66 ± 28.141913.50 ± 20.9639.8639.4539.65
J1629+400742.065009.08 ± 0.06133.94 ± 3.59211.51 ± 2.405004.52 ± 0.34−272.60 ± 22.15830.54 ± 19.1641.7141.7039.06
J1633+471841.655006.57 ± 0.05−16.35 ± 2.99255.23 ± 1.805001.91 ± 0.23−278.77 ± 14.971185.89 ± 13.1740.6740.6439.46
J1634+480941.965002.47 ± 0.44−261.84 ± 26.35389.19 ± 24.554993.44 ± 1.69−541.39 ± 125.141109.74 ± 56.2841.4241.4040.03
J1703+454041.555010.34 ± 1.23209.39 ± 73.65335.60 ± 49.105004.73 ± 3.96−335.68 ± 309.56645.82 ± 110.7740.2040.1639.04
J1709+234842.145008.43 ± 0.0995.02 ± 5.39401.88 ± 3.595003.06 ± 0.39−321.53 ± 26.94786.83 ± 9.5839.5539.550.00
J1713+352340.984995.59 ± 0.88−673.79 ± 52.691641.35 ± 47.900.000.000.0039.5439.5138.33
J2242+294341.345006.77 ± 0.01−4.46 ± 0.72226.27 ± 3.754997.42 ± 1.64−564.31 ± 98.631089.89 ± 279.0638.8838.6238.54
J2314+224342.315006.86 ± 0.811.02 ± 48.50602.11 ± 25.754989.64 ± 2.20−1031.42 ± 178.431579.30 ± 55.6940.4840.4239.58
J1029+555641.825006.96 ± 0.226.97 ± 13.39335.37 ± 37.375004.16 ± 1.65−160.73 ± 98.892004.90 ± 387.5438.4638.460.00
J1228+501741.835008.56 ± 0.03102.76 ± 1.71218.21 ± 5.785003.58 ± 0.27−195.21 ± 16.48978.24 ± 47.5239.0038.9138.27
J1232+495741.665003.32 ± 1.28−211.09 ± 76.44578.97 ± 156.314994.82 ± 4.45−719.77 ± 266.801142.04 ± 297.7738.2338.230.00
J1509+613740.945003.87 ± 0.12−178.29 ± 7.34320.40 ± 24.245005.87 ± 1.91−58.51 ± 114.531690.48 ± 189.3637.6837.680.00
J1510+554741.295006.85 ± 0.030.67 ± 1.56123.77 ± 8.755004.09 ± 0.29−164.99 ± 17.14547.56 ± 34.9137.6937.4837.27
J1522+393440.305007.43 ± 0.3635.10 ± 21.3321.12 ± 79.065001.67 ± 3.29−310.03 ± 197.071306.38 ± 324.2337.9537.8837.12
J1641+345441.455006.61 ± 0.07−14.01 ± 4.46448.44 ± 19.385006.19 ± 0.62−39.20 ± 37.301274.45 ± 136.6938.9238.7938.33
Table 3. Results of the correlations between the [O III] core shift and the physical properties of the sources. Columns: (1) First tested sample; (2) second tested sample; (3) Pearson correlation coefficient; (4) p-value of the Pearson coefficient; (5) Spearman’s rank correlation coefficient; (6) p-value of the Spearman coefficient.
Table 3. Results of the correlations between the [O III] core shift and the physical properties of the sources. Columns: (1) First tested sample; (2) second tested sample; (3) Pearson correlation coefficient; (4) p-value of the Pearson coefficient; (5) Spearman’s rank correlation coefficient; (6) p-value of the Spearman coefficient.
Q1Q2rPp-Value (P)rSp-Value (S)
L [ O III ] v c −0.050.680.100.38
L [ O III ] |v c |0.070.540.240.03
L bol v c −0.100.370.010.97
L bol |v c |0.150.180.280.01
L 5 GHz v c 0.010.970.010.92
L 5 GHz |v c |0.030.780.390.0004
RLv c 0.180.120.030.81
RL|v c |0.120.290.360.001
Table 4. Sample studied in this work. Columns: (1) short name; (2) NED alias; (3) right ascension; (4) declination; (5) redshift; (6) jet (Y stands for yes, N for no); (7) radio morphology (C stands for compact, I for intermediate, E for extended); (8) spectral index for jetted sources (F stands for flat spectrum, S stands for steep spectrum).
Table 4. Sample studied in this work. Columns: (1) short name; (2) NED alias; (3) right ascension; (4) declination; (5) redshift; (6) jet (Y stands for yes, N for no); (7) radio morphology (C stands for compact, I for intermediate, E for extended); (8) spectral index for jetted sources (F stands for flat spectrum, S stands for steep spectrum).
Short NameNED AliasR.A.Dec.zJetMorphSp.Ind.
J0006+2012Mrk 33500 06 19.50+20 12 10.00.026NC
J0146−00402MASX J01464481−004042601 46 44.80−00 40 43.00.083NC
J0347+0105IRAS 03450+005503 47 40.20+01 05 14.00.031NE
J0629−0545IRAS 06269−054306 29 24.70−05 45 30.00.117YIS
J0632+6340UGC 347806 32 47.20+63 40 25.00.013NE
J0706+3901FBQS J0706+390107 06 25.15+39 01 51.60.086YIS
J0713+3820FBQS J0713+382007 13 40.29+38 20 40.10.123YES
J0744+5149NVSS J074402+51491707 44 02.24+51 49 17.50.460YCF
J0752+2617FBQS J0752+261707 52 45.60+26 17 36.00.082NC
J0758+3920FBQS J075800.0+39202907 58 00.05+39 20 29.10.096YES
J0804+3853FBQS J0804+385308 04 09.24+38 53 48.70.212YES
J0806+7248RGB J0806+72808 06 38.96+72 48 20.40.098YIS
J0814+5609SDSS J081432.11+560956.608 14 32.13+56 09 56.60.510YIF
J0849+5108SBS 0846+51308 49 57.99+51 08 28.80.585YCF
J0850+4626SDSS J085001.17+462600.508 50 01.17+46 26 00.50.524YCF
J0902+0443SDSS J090227.16+044309.509 02 27.15+04 43 09.40.533YCF
J0913+3658RX J0913.2+365809 13 13.70+36 58 17.00.107NI
J0925+5217Mrk 11009 25 12.90+52 17 11.00.035NE
J0926+1244Mrk 70509 26 03.30+12 44 04.00.029NI
J0937+3615SDSS J093703.03+361537.209 37 03.01+36 15 37.30.180NE
J0945+1915SDSS J094529.23+191548.809 45 29.21+19 15 48.90.284YCF
J0948+5029Mrk 12409 48 42.60+50 29 31.00.056NC
J0952−0136Mrk 123909 52 19.10−01 36 43.00.020YIS
J0957+2433RX J0957.1+243309 57 07.20+24 33 16.00.082NC
J1029+5556SDSS J102906.69+555625.210 29 06.69+55 56 25.20.451YCF
J1031+4234SDSS J103123.73+423439.310 31 23.73+42 34 39.40.377YCS
J1034+3938KUG 1031+39810 34 38.60+39 38 28.00.042YIS
J1037+0036SDSS J103727.45+003635.610 37 27.45+00 36 35.80.595YCF
J1038+4227SDSS J103859.58+422742.310 38 59.59+42 27 42.00.220NE
J1047+4725SDSS J104732.68+472532.010 47 32.65+47 25 32.20.799YEF
J1048+2222SDSS J104816.57+222238.910 48 16.56+22 22 40.10.330YIS
J1102+2239SDSS J110223.38+223920.711 02 23.36+22 39 20.70.453YIS
J1110+3653SDSS J111005.03+365336.311 10 05.03+36 53 36.10.629YIF
J1114+3241B2 1111+3211 14 38.89+32 41 33.40.189YCF
J1121+5351SBS 1118+54111 21 08.60+53 51 21.00.103NI
J1138+3653SDSS J113824.54+365327.111 38 24.54+36 53 27.10.356YIS
J1146+3236SDSS J114654.28+323652.311 46 54.30+32 36 52.20.465YCF
J1159+2838SDSS J115917.32+283814.511 59 17.31+28 38 14.80.210YIS
J1203+4431NGC 405112 03 09.60+44 31 53.00.002NE
J1209+3217RX J1209.7+321712 09 45.20+32 17 01.00.144NI
J1215+5442SBS 1213+549A12 15 49.40+54 42 24.00.150NI
J1218+2948Mrk 76612 18 26.50+29 48 45.80.013NE
J1227+3214SDSS J122749.14+321458.912 27 49.15+32 14 59.00.136YIS
J1228+5017SDSS J122844.81+501751.212 28 44.81+50 17 51.20.262YIF
J1232+4957SDSS J123220.11+495721.812 32 20.11+49 57 21.80.262YCF
J1238+3942SDSS J123852.12+394227.812 27 49.15+39 42 27.60.622YCF
J1242+3317WAS 6112 42 10.60+33 17 03.00.044NE
J1246+0222PG 1244+02612 46 35.20+02 22 09.00.048NE
J1246+0238SDSS J124634.65+023809.012 46 34.68+02 38 09.00.367YCF
J1302+1624Mrk 78313 02 58.80+16 24 27.00.067YES
J1305+5116SDSS J130522.74+511640.213 05 22.75+51 16 40.30.788YEF
J1317+6010SBS 1315+60413 17 50.30+60 10 41.00.137NE
J1333+4141SDSS J133345.47+414127.713 33 45.47+41 41 28.20.225YCS
J1337+2423IRAS 13349+243813 37 18.70+24 23 03.00.108NI
J1346+3121SDSS J134634.97+312133.713 46 35.07+31 21 33.90.246YCF
J1355+5612SBS 1353+56413 55 16.50+56 12 45.00.122NC
J1358+2658SDSS J135845.38+265808.413 58 45.40+26 58 08.30.331YIS
J1402+2159RX J1402.5+215914 02 34.40+21 59 52.00.066NI
J1421+2824SDSS J142114.05+282452.814 21 14.07+28 24 52.20.539YCF
J1443+4725SDSS J144318.56+472556.714 43 18.56+47 25 56.70.705YCF
J1505+0326SDSS J150506.47+032630.815 05 06.47+03 26 30.80.408YCF
J1509+6137SDSS J150916.18+613716.715 09 16.18+61 37 16.70.201YCF
J1510+5547SDSS J151020.06+554722.015 10 20.06+55 47 22.00.150YEF
J1522+3934SDSS J152205.41+393441.315 22 05.41+39 34 41.30.077YIF
J1536+5433Mrk 48615 36 38.30+54 33 33.00.039NI
J1537+4942SDSS J153732.61+494247.515 37 32.60+49 42 48.00.280NC
J1548+3511SDSS J154817.92+351128.015 48 17.92+35 11 28.40.479YCF
J1555+1911Mrk 29115 55 07.90+19 11 33.00.035NE
J1559+3501Mrk 49315 59 09.60+35 01 47.00.031NE
J1612+4219SDSS J161259.83+421940.316 12 59.83+42 19 40.00.233YEF
J1629+4007SDSS J162901.30+400759.916 29 01.31+40 07 59.60.272YCF
J1633+4718SDSS J163323.58+471858.916 33 23.58+47 18 59.00.116YIF
J1634+4809SDSS J163401.94+480940.216 34 01.94+48 09 40.10.495YCF
J1641+3454SDSS J164100.10+345452.716 41 00.10+34 54 52.70.164YEF
J1703+4540SDSS J170330.38+454047.117 03 30.38+45 40 47.20.060YIS
J1709+2348SDSS J170907.80+234837.717 09 07.82+23 48 38.20.254YCS
J1713+3523FBQS J1713+352317 13 04.48+35 23 33.40.083YIS
J2242+2943Ark 56422 42 39.30+29 43 31.00.025NE
J2314+2243RX J2314.9+224323 14 55.70+22 43 25.00.169YIS
Table 5. Results of the statistical tests for the spectral index samples. Columns: (1) Tested sample; (2) median velocity; (3) standard deviation of the velocity; (4) interquartile range; (5) p-value of the K-S tests against the other sample, observed in top row, simulated in middle row, and standard deviation of the simulated p-value in the bottom row; (6) p-value of the A-D tests against the other sample, observed in top row, simulated in middle row, and standard deviation of the simulated p-value in the bottom row.
Table 5. Results of the statistical tests for the spectral index samples. Columns: (1) Tested sample; (2) median velocity; (3) standard deviation of the velocity; (4) interquartile range; (5) p-value of the K-S tests against the other sample, observed in top row, simulated in middle row, and standard deviation of the simulated p-value in the bottom row; (6) p-value of the A-D tests against the other sample, observed in top row, simulated in middle row, and standard deviation of the simulated p-value in the bottom row.
SampleMedian
(km s 1 )
Std
(km s 1 )
IQR
(km s 1 )
K-S obs
K-S sim
K-S std
A-D obs
A-D sim
A-D std
v c (NJ)−26.548.571.00.06
0.06
0.02
0.01
0.02
0.01
v c (J)−11.9226.4249.10.06
0.06
0.02
0.01
0.02
0.01
FWHM c (NJ)329.485.969.10.05
0.03
0.02
0.03
0.02
0.01
FWHM c (J)385.7289.8328.30.05
0.03
0.02
0.03
0.02
0.01
v w (NJ)−268.9127.798.80.06
0.10
0.09
0.05
0.07
0.05
v w (J)−335.7244.1298.60.06
0.10
0.09
0.05
0.07
0.05
FWHM w (NJ)868.6346.2311.40.03
0.06
0.05
0.07
0.07
0.04
FWHM w (J)1180.3427.8607.70.03
0.06
0.05
0.07
0.07
0.04
Table 6. Results of the correlations tested throughout this work. Columns: (1) First tested sample; (2) second tested sample; (3) Pearson correlation coefficient; (4) p-value.
Table 6. Results of the correlations tested throughout this work. Columns: (1) First tested sample; (2) second tested sample; (3) Pearson correlation coefficient; (4) p-value.
Q1Q2rp-Value
v c t o t FWHM c t o t −0.708 × 10 13
v c N J FWHM c N J 0.020.93
v c J FWHM c J −0.734 × 10 9
v w t o t v c t o t 0.010.92
v w N J v c N J 0.200.35
v w J v c J 0.240.13
v w t o t FWHM c t o t −0.090.43
v w N J FWHM c N J −0.060.79
v w J FWHM c J −0.508 × 10 4
FWHM w t o t v c t o t 0.020.85
FWHM w N J v c N J −0.210.32
FWHM w J v c J −0.360.02
FWHM w t o t FWHM c t o t −0.010.95
FWHM w N J FWHM c N J 0.340.09
FWHM w J FWHM c J 0.552 × 10 4
v w t o t FWHM w t o t −0.471 × 10 5
v w N J FWHM w N J −0.510.01
v w J FWHM w J −0.080.61
v c N J L i n t 5 G H z −0.300.12
v c J L i n t 5 G H z 0.030.84
v c N J L [ O I I I ] 0.220.26
v c J L [ O I I I ] −0.010.94
Table 7. Statistical tests results for the morphological samples. Columns: (1) Tested sample; (2) median velocity; (3) standard deviation of the velocity; (4) interquartile range; (5) p-value of the K-S tests against the corresponding C sample, observed in top row, simulated in middle row, and standard deviation of the simulated p-value in bottom row; (6) p-value of the A-D tests against the corresponding C sample, observed in top row, simulated in middle row, and standard deviation of the simulated p-value in bottom row; (7) K-S tests against the corresponding I sample, rows as before; (8) A-D tests against the corresponding I sample, rows as before; (9) K-S tests against the corresponding E sample, rows as before; (10) A-D tests against the corresponding E sample, rows as before.
Table 7. Statistical tests results for the morphological samples. Columns: (1) Tested sample; (2) median velocity; (3) standard deviation of the velocity; (4) interquartile range; (5) p-value of the K-S tests against the corresponding C sample, observed in top row, simulated in middle row, and standard deviation of the simulated p-value in bottom row; (6) p-value of the A-D tests against the corresponding C sample, observed in top row, simulated in middle row, and standard deviation of the simulated p-value in bottom row; (7) K-S tests against the corresponding I sample, rows as before; (8) A-D tests against the corresponding I sample, rows as before; (9) K-S tests against the corresponding E sample, rows as before; (10) A-D tests against the corresponding E sample, rows as before.
SampleMedian
(km s 1 )
Std
(km s 1 )
IQR
(km s 1 )
K-S C obs
K-S C sim
K-S C std
A-D C obs
A-D C sim
A-D C std
K-S I obs
K-S I sim
K-S I std
A-D I obs
A-D I sim
A-D I std
K-S E obs
K-S E sim
K-S E std
A-D E obs
A-D E sim
A-D E std
v c (C)−11.6210.4125.70.87
0.87
0.12
0.89
0.74
0.11
0.36
0.33
0.15
0.29
0.30
0.08
v c (I)−31.3207.2116.10.87
0.87
0.12
0.89
0.74
0.11
0.42
0.52
0.18
0.32
0.36
0.10
v c (E)−11.0101.097.00.36
0.33
0.15
0.29
0.30
0.08
0.42
0.52
0.18
0.32
0.36
0.10
FWHM c (C)323.8232.2147.30.36
0.20
0.18
0.15
0.15
0.09
0.63
0.60
0.15
0.65
0.61
0.12
FWHM c (I)385.0285.0193.40.36
0.20
0.18
0.15
0.15
0.09
0.64
0.43
0.20
0.56
0.57
0.12
FWHM c (E)341.2194.6176.30.63
0.60
0.15
0.65
0.61
0.12
0.64
0.43
0.20
0.56
0.57
0.12
v w (C)−309.8191.9272.20.95
0.69
0.23
0.92
0.70
0.21
0.68
0.44
0.28
0.54
0.57
0.25
v w (I)−302.4229.5276.60.95
0.69
0.23
0.92
0.70
0.21
0.59
0.59
0.24
0.60
0.55
0.24
v w (E)−270.3210.8199.10.68
0.44
0.28
0.54
0.57
0.25
0.59
0.59
0.24
0.60
0.55
0.24
FWHM w (C)965.2425.0472.40.65
0.89
0.14
0.81
0.93
0.17
1.00
0.99
0.05
1.00
1.00
0.10
FWHM w (I)1082.1404.8561.90.65
0.89
0.14
0.81
0.93
0.17
0.75
0.93
0.16
0.89
0.84
0.15
FWHM w (E)999.8383.0567.41.00
0.99
0.05
1.00
1.00
0.10
0.75
0.93
0.16
0.89
0.84
0.15
Table 8. Results of the statistical tests for the spectral index samples. Columns: (1) Tested sample; (2) median velocity; (3) standard deviation of the velocity; (4) interquartile range; (5) p-value of the K-S tests against the other sample, observed in top row, simulated in middle row, and standard deviation of the simulated p-value in bottom row; (6) p-value of the A-D tests against the other sample, observed in top row, simulated in middle row, and standard deviation of the simulated p-value in bottom row.
Table 8. Results of the statistical tests for the spectral index samples. Columns: (1) Tested sample; (2) median velocity; (3) standard deviation of the velocity; (4) interquartile range; (5) p-value of the K-S tests against the other sample, observed in top row, simulated in middle row, and standard deviation of the simulated p-value in bottom row; (6) p-value of the A-D tests against the other sample, observed in top row, simulated in middle row, and standard deviation of the simulated p-value in bottom row.
SampleMedian
(km s 1 )
Std
(km s 1 )
IQR
(km s 1 )
K-S obs
K-S sim
K-S std
A-D obs
A-D sim
A-D std
v c (S)−4.97226.6252.10.71
0.92
0.14
0.97
0.95
0.11
v c (F)−12.8225.8226.20.71
0.92
0.14
0.97
0.95
0.11
FWHM c (S)408.9311.6266.50.15
0.14
0.11
0.25
0.28
0.09
FWHM c (F)321.1268.0253.50.15
0.14
0.11
0.25
0.28
0.09
v w (S)−362.1246.5355.60.50
0.50
0.23
0.26
0.32
0.23
v w (F)−326.1229.6339.90.50
0.50
0.23
0.26
0.32
0.23
FWHM w (S)1062.6399.4678.30.68
0.68
0.05
0.67
0.69
0.10
FWHM w (F)1180.3445.3467.80.68
0.68
0.05
0.67
0.69
0.10
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Berton, M.; Järvelä, E. Jet-Induced Feedback in the [O III] Lines of Early Evolution Stage Active Galactic Nuclei. Universe 2021, 7, 188. https://doi.org/10.3390/universe7060188

AMA Style

Berton M, Järvelä E. Jet-Induced Feedback in the [O III] Lines of Early Evolution Stage Active Galactic Nuclei. Universe. 2021; 7(6):188. https://doi.org/10.3390/universe7060188

Chicago/Turabian Style

Berton, Marco, and Emilia Järvelä. 2021. "Jet-Induced Feedback in the [O III] Lines of Early Evolution Stage Active Galactic Nuclei" Universe 7, no. 6: 188. https://doi.org/10.3390/universe7060188

APA Style

Berton, M., & Järvelä, E. (2021). Jet-Induced Feedback in the [O III] Lines of Early Evolution Stage Active Galactic Nuclei. Universe, 7(6), 188. https://doi.org/10.3390/universe7060188

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop