Search

The last few years have seen the confirmation of several trends associated with the quasar main sequence. The idea of a main sequence for quasars is relatively recent, and its full potential for the observational classification and contextualization of quasar properties has yet to be fully exploited. The main sequence drivers are discussed in terms of the properties of extreme objects. We briefly summarize developments that constrain the viewing angle of the accretion disk in a particular class of quasars (extreme Population B, radiating at low Eddington ratio), as well as inferences on the chemical composition of the broad line emitting gas, and on the nature of radio emission along the quasar main sequence., Comment: Paper presented at the XIII SCSLSA, to appear in Astronomische Nachrichten

Improving our understanding of the nuclear properties of high-Eddington ratio ($\lambda_\mathrm{Edd}$) active galactic nuclei (AGN) is necessary since the bulk of X-ray spectroscopic studies have been focused on low-Eddington AGN. We present here the X-ray spectral analysis of 14 radio-quiet, $\lambda_\mathrm{Edd}\gtrsim1$ AGN at $0.4\leq z \leq 0.75$, observed with XMM-Newton. Optical/UV data from simultaneous Optical Monitor observations have been also considered. These AGN have been selected to have relatively high values of black hole mass ($M_\mathrm{BH}\sim10^{8-8.5}M_\odot$) and bolometric luminosity ($L_\mathrm{bol} \sim 10^{46}$ erg s$^{-1}$), in order to complement previous studies of high-Eddington AGN at lower $M_\mathrm{BH}$ and $L_\mathrm{bol}$. We studied the relation between $\lambda_\mathrm{Edd}$ and other key X-ray spectral parameters, such as the photon index of the power-law continuum $\Gamma$, the X-ray bolometric correction $k_\mathrm{bol,X}$ and $\alpha_\mathrm{ox}$. Despite the homogeneous optical and SMBH accretion properties, the X-ray properties of these high-Eddington AGN are quite heterogeneous. We measured values of $\Gamma$ comprised between 1.3 and 2.5, at odds with the expectations based on previously reported $\Gamma-\lambda_\mathrm{Edd}$ relations, by which $\Gamma\geq2$ would be an ubiquitous hallmark of AGN with $\lambda_\mathrm{Edd}\sim1$. We found that $\sim30\%$ of the sources are X-ray weak, with an X-ray emission about a factor of $\sim10-80$ fainter than that of typical AGN at similar UV luminosities. The X-ray weakness seems to be intrinsic and not due to intervening obscuration. This may indicate that high-Eddington AGN commonly undergo periods of intrinsic X-ray weakness. Furthermore, results from a follow-up monitoring with Swift of one of these X-ray weak AGN suggest that these periods can last for several years., Comment: 17 pages, 8 figures (with 2 more figures in Appendix A and B). Accepted for publication in A&A

We review the properties of the established Scaling Relations (SRs) of galaxies and active galactic nuclei (AGN), focusing on their origin and expected evolution back in time, providing a short history of the most important progresses obtained up to now and discussing the possible future studies. We also try to connect the observed SRs with the physical mechanisms behind them, examining to what extent current models reproduce the observational data. The emerging picture clarifies the complexity intrinsic to the galaxy formation and evolution process as well as the basic uncertainties still affecting our knowledge of the AGN phenomenon. At the same time, however, it suggests that the detailed analysis of the SRs can profitably contribute to our understanding of galaxies and AGN., Comment: 116, 15 figures

The quasar Main Sequence (MS) appears to be an incredibly powerful tool to organize the diversity in large samples of type-1 quasars but the most important physical parameters governing it are still unclear. Here we investigate the origin of the broadening and of a defining feature of Population B sources: a strong redward asymmetry of the Balmer emission lines. We focus on a prototypical source, Fairall 9. Spectropolarimetric data of the Fairall 9 broad H$\beta$ and H$\alpha$ profiles allowed for a view of the geometric and dynamical complexity of the line emitting regions. Measurements (1) provided evidence of rotational motion; (2) were helpful to test the presence of polar and equatorial scatterers, and their association with non-virial motions., Comment: 22 pages, 3 tables and 14 figures. Accepted for publication in MNRAS

PKS\,1351$-$018 at a redshift of $z=3.71$ is one of the most luminous, steady synchrotron sources with a luminosity $> 10^{47}$\,erg~s$^{-1}$. The synchrotron luminosity does not seem to vary by more than $\sim 25\%$ over 35 years. In order to appreciate this remarkable behavior, if it were at $z=0.5$, it would have a flux density at 15 GHz in a range of $110 - 137$\,Jy over 11 yrs. In spite of this steady behavior, two strong $\gamma$-ray flares $\lesssim 10^{49}$\,erg~s$^{-1}$ were detected in 2011 and 2016. There is blazar-like behavior coexisting with the steady behavior. This study is aimed at elucidating the dual nature of this source. We find that the radio source is extremely compact with a bright core and a steep spectrum secondary component, 12\,mas away, that appears to be constant in position and flux density in six epochs from 1995 to 2018. We estimate that a jet with a time averaged power of $(5.2 \pm 3.2) \times 10^{45}$\,erg~s$^{-1}$ terminates in this lobe that is advancing $\gtrsim 0.9 c$ at a deprojected distance of $1-3$\,kpc from the central engine. This is the rare case of a young ($\sim 6000$\,yr), very powerful radio source that is viewed a few degrees from the jet axis. We find evidence of a high velocity (4000\,km~s$^{-1}$), high ionization wind emanating form a luminous quasar. The young radio jet appears to experience modest bending as it navigates through the intense quasar environment., Comment: To appear in ApJ

Narrow-line Seyfert 1 (NLS1) galaxies are a class of active galactic nuclei (AGN) that, in some cases, can harbor powerful relativistic jets. One of them, PKS 2004-447, shows gamma-ray emission, and underwent its first recorded multifrequency flare in 2019. However, past studies revealed that in radio this source can be classified as a compact steep-spectrum source (CSS), suggesting that, unlike other gamma-ray sources, the relativistic jets of PKS 2004-447 have a large inclination with respect to the line of sight. We present here a set of spectroscopic observations of this object, aimed at carefully measuring its black hole mass and Eddington ratio, determining the properties of its emission lines, and characterizing its long term variability. We find that the black hole mass is $(1.5\pm0.2)\times10^7$ M$_\odot$, and the Eddington ratio is 0.08. Both values are within the typical range of NLS1s. The spectra also suggest that the 2019 flare was caused mainly by the relativistic jet, while the accretion disk played a minor role during the event. In conclusion, we confirm that PKS 2004-447 is one of the rare examples of gamma-ray emitting CSS/NLS1s hybrid, and that these two classes of objects are likely connected in the framework of AGN evolution., Comment: 14 pages, 8 figures, 4 tables, resubmitted to A&A after minor comments from the referee

The classification of the optical spectra of active galactic nuclei (AGN) into different types is well founded on AGN physics, but it involves some degree of human oversight and cannot be reliably scaled to large data sets. Machine learning (ML) tackles such a classification problem in a fast and reproducible way, but is often perceived as a black box. However, ML interpretability and explainability are active research areas in computer science, increasingly providing us with tools to alleviate this issue. We applied ML interpretability tools to a classifier trained to predict AGN type from spectra, to demonstrate the use of such tools in this context. We trained a support-vector machine on 3346 high-quality, low redshift AGN spectra from SDSS DR15 with an existing reliable classification as type 1, type 2, or intermediate type. On a selection of test-set spectra, we computed the gradient of the predicted class probability and we built saliency maps. We also visualized the high-dimensional space of AGN spectra using t-distributed stochastic neighbor embedding (t-SNE), showing where the spectra for which we computed a saliency map are located. Regions that affect the predicted AGN type often coincide with physically relevant features, such as spectral lines. t-SNE visualization shows good separability of type 1 and type 2 spectra, while intermediate-type spectra either lie in-between as expected or appear mixed with type 2 spectra. Saliency maps show why a given AGN type was predicted by our classifier, resulting in a physical interpretation in terms of regions of the spectrum that affected its decision, making it no longer a black box. These regions coincide with those used by human experts such as relevant spectral lines, and are even used in a similar way, with the classifier e.g. effectively measuring the width of a line by weighing its center and its tails oppositely., Comment: Under review by A&A

We study the limit behaviour of singularly-perturbed elliptic functionals of the form \[ \mathcal F_k(u,v)=\int_A v^2\,f_k(x,\nabla u)\.dx+\frac{1}{\varepsilon_k}\int_A g_k(x,v,\varepsilon_k\nabla v)\.dx\,, \] where $u$ is a vector-valued Sobolev function, $v \in [0,1]$ a phase-field variable, and $\varepsilon_k>0$ a singular-perturbation parameter, i.e., $\varepsilon_k \to 0$, as $k\to +\infty$. Under mild assumptions on the integrands $f_k$ and $g_k$, we show that if $f_k$ grows superlinearly in the gradient-variable, then the functionals $\mathcal F_k$ $\Gamma$-converge (up to subsequences) to a brittle energy-functional, i.e., to a free-discontinuity functional whose surface integrand does not depend on the jump-amplitude of $u$. This result is achieved by providing explicit asymptotic formulas for the bulk and surface integrands which show, in particular, that volume and surface term in $\mathcal F_k$ decouple in the limit. The abstract $\Gamma$-convergence analysis is complemented by a stochastic homogenisation result for stationary random integrands.

In this second paper in the series, we carefully analyze the observational properties of the optical FeII and NIR CaII triplet in Active Galactic Nuclei, as well as the luminosity, black hole mass, and Eddington ratio in order to define the driving mechanism behind the properties of our sample. The CaII shows an inverse Baldwin effect, bringing out the particular behavior of this ion with respect to the other low-ionization lines such as H$\beta$. We performed a Principal Component Analysis, where 81.2% of the variance can be explained by the first three principal components drawn from the FWHMs, luminosity, and equivalent widths. The first principal component (PC1) is primarily driven by the combination of black hole mass and luminosity with a significance over 99.9%, which in turn is reflected in the strong correlation of the PC1 with the Eddington ratio. The observational correlations are better represented by the Eddington ratio, thus it could be the primary mechanism behind the strong correlations observed in the CaII-FeII sample. Since calcium belongs to the $\alpha$-elements, the FeII/CaII flux ratio can be used as a chemical clock for determining the metal content in AGN and trace the evolution of the host galaxies. We confirm the de-enhancement of the ratio FeII/CaII by the Eddington ratio, suggesting a metal enrichment of the BLR in intermediate-$z$ with respect to low-$z$ objects. A larger sample, particularly at $z$>2, is needed to confirm the present results., Comment: 35 pages, 16 figures, 4 tables. Accepted for publication in ApJ

Recent measurements of the parameters of the Concordance Cosmology Model ($\Lambda$CDM) done in the low-redshift Universe with Supernovae Ia/Cepheids, and in the distant Universe done with Cosmic Microwave Background (CMB) imply different values for the Hubble constant (67.4 $\pm$ 0.5 km s$^{-1}$ Mpc$^{-1}$ from Planck vs 74.03 $\pm$ 1.42 km s$^{-1}$ Mpc$^{-1}$, Riess et al. 2019). This Hubble constant tension implies that either the systematic errors are underestimated, or the $\Lambda$CDM does not represent well the observed expansion of the Universe. Since quasars - active galactic nuclei - can be observed in the nearby Universe up to redshift z $\sim$ 7.5, they are suitable to estimate the cosmological properties in a large redshift range. Our group develops two methods based on the observations of quasars in the late Universe up to redshift z$\sim $4.5, with the objective to determine the expansion rate of the Universe. These methods do not yet provide an independent measurement of the Hubble constant since they do not have firm absolute calibration but they allow to test the $\Lambda$CDM model, and so far no departures from this model were found., Comment: 6 pages, 2 figures, contribution to the Proceedings of the 100 years of the Polish Physical Society Congress, 16-18 October 2020, Warsaw, Poland. http://info.ifpan.edu.pl/APP/

We present an analysis of UV spectra of 13 quasars believed to belong to extreme Population A (xA) quasars, aimed at the estimation of the chemical abundances of the broad line emitting gas. Metallicity estimates for the broad line emitting gas of quasars are subject to a number of caveats, although present data suggest the possibility of an increase along the quasar main sequence along with prominence of optical Fe II emission. Extreme Population A sources with the strongest Fe II emission offer several advantages with respect to the quasar general population, as their optical and UV emission lines can be interpreted as the sum of a low-ionization component roughly at quasar rest frame (from virialized gas), plus a blueshifted excess (a disk wind), in different physical conditions. Specifically, in terms of ionization parameter, cloud density, metallicity and column density. Capitalizing on these results, we analyze the component at rest frame and the blueshifted one, exploiting the dependence (of several intensity line ratios on metallicity $Z$). We find that the validity of intensity line ratios as metallicity indicators depends on the physical conditions. We apply the measured diagnostic ratios to estimate the physical properties of sources such as density, ionization, and metallicity of the gas. Our results confirm that the two regions (the low-ionization component and the blue-shifted excess) of different dynamical conditions also show different physical conditions and suggest metallicity values that are high, and probably the highest along the quasar main sequence, with $Z \gtrsim 10 Z_{\odot}$. We found some evidence of an overabundance of Aluminium with respect to Carbon, possibly due to selective enrichment of the broad line emitting gas by supernova ejecta., Comment: 39 pages, will appear in ApJ

We present a study of quasars with very redward asymmetric (RA) ultraviolet (UV) broad emission lines (BELs). An excess of redshifted emission has been previously shown to occur in the BELs of radio loud quasars and is most extreme in certain blazars. Paradoxically, blazars are objects that are characterized by a highly relativistic blue-shifted outflow towards Earth. We show that the red emitting gas resides in a very broad component (VBC) that is typical of Population B quasars that are defined by a wide H$\beta$ BEL profile. Empirically, we find that RA BEL blazars have both low Eddington rates ($\lesssim1\%$) and an inordinately large (order unity) ratio of long term time averaged jet power to accretion luminosity. The latter circumstance has been previously shown to be associated with a depressed extreme UV ionizing continuum. Both properties conspire to produce a low flux of ionizing photons, two orders of magnitude less than typical Population B quasars. We use CLOUDY models to demonstrate that a weak ionizing flux is required for gas near the central black hole to be optimally ionized to radiate BELs with high efficiency (most quasars over-ionize nearby gas, resulting in low radiative efficiency). The large gravitational redshift and transverse Doppler shift results in a VBC that is redshifted by $\sim 2000 -5000$~km~s$^{-1}$ with a correspondingly large line width. The RA BELs result from an enhanced efficiency (relative to typical Population B quasars) to produce a luminous, redshifted VBC near the central black hole., Comment: To appear in ApJ

3C 82 at a redshift of 2.87 is the most distant 3C (Third Cambridge Catalogue) quasar. Thus, it is a strong candidate to have the most luminous radio lobes in the Universe. 3C 82 belongs to the class of compact steep spectrum radio sources. We use single dish and interferometric radio observations in order to model the plasma state of these powerful radio lobes. It is estimated that the long-term time-averaged jet power required to fill these lobes with leptonic plasma is $\overline{Q} \approx 2.66 \pm 1.33 \times 10^{47} \rm{ergs/sec}$, among the largest time averaged jet powers from a quasar. Positing protonic lobes is not tenable since they would require two orders of magnitude more mass transport to the lobes than was accreted to the central black hole during their formation. The first high signal to noise optical spectroscopic observation obtained of this object indicates that there is a powerful high ionization broad line wind with a kinetic power $\sim 10^{45} \rm{ergs/sec}$ and a velocity $\sim 0.01$c. We also estimate from the broad lines in 2018 and the UV continuum in three epochs spread out over three decades that the accretion flow bolometric luminosity is $L_{\rm{bol}} \approx 3.2-5.8 \times 10^{46} \rm{ergs/sec}$. The ratio of $\overline{Q}/L_{\rm{bol}}\approx 6.91 \pm 3.41$, is perhaps the largest of any known quasar. Extremely powerful jets tend to strongly suppress powerful winds of ionized baryonic matter. Consequently, 3C 82 provides a unique laboratory for studying the dynamical limits of the central engine of outflow initiation in quasars., Comment: To appear in ApJ

A statistical study of intermediate Palomar Transient Factory supernovae (SNe) in Type 1 AGN has shown a major deficit of supernovae around Type 1 AGN host galaxies, with respect to Type 2 AGN hosts. The aim of this work is to test whether there is any preference for Type 1 AGN to host SN of a specific kind. Through the analysis of SN occurrence and their type (thermonuclear vs core-collapse), we can directly link the type of stars producing the SN events, thus this is an indirect way to study host galaxies in Type 1 AGN. We examine the detection fractions of SNe, the host galaxies and compare the sample properties to typical host galaxies in the Open Supernova Catalog (OSC; Guillochon et al. 2017). The majority of the host galaxies in the AGN sample are late-type, similar to typical galaxies hosting SN within the OSC. The findings are supportive of a deficiency of SNe near Type 1 AGN, although we cannot with certainty assess the overall detection fractions of SNe in Type 1 AGN relative to other SN host galaxies. We can state that Type 1 AGN has equal detection fractions of thermonuclear vs core-collapse SNe. However, we note the possibility of a higher detection rate of core-collapse supernovae in Type-1 AGN with insecure AGN classifications., Comment: 12 pages, 5 figures, 3 tables. Accepted for publication in MNRAS

We analyze the optical properties of Radio-Loud quasars along the Main Sequence (MS) of quasars. A sample of 355 quasars selected on the basis of radio detection was obtained by cross-matching the FIRST survey at 20cm and the SDSS DR12 spectroscopic survey. We consider the nature of powerful emission at the high-FeII end of the MS. At variance with the classical radio-loud sources which are located in the Population B domain of the MS optical plane, we found evidence indicating a thermal origin of the radio emission of the highly accreting quasars of Population A., Comment: 3 pages, 2 figures; Proceedings paper of the IAU symposium 356 "Nuclear Activity in Galaxies Across Cosmic Time" (Ethiopia) accepted to be published under the Cambridge University Press, eds. M. Povic, P. Marziani, J. Masegosa, H. Netzer, S. H. Negu, and S. B. Tessema

Optical Fe II emission is a strong feature in quasar spectra originating in the broad-line region (BLR). The difficulty in understanding the complex Fe II pseudo-continuum has led us to search for other reliable, simpler ionic species such as Ca II. In this first part of the series, we confirm the strong correlation between the strengths of two emission features, the optical Fe II and the NIR Ca II, both from observations and photoionization modelling. With the inclusion of an up-to-date compilation of observations with both optical Fe II and NIR Ca II measurements, we span a wider and more extended parameter space and confirm the common origin of these two spectral features with our photoionization models using CLOUDY. Taking into account the effect of dust into our modelling, we constrain the BLR parameter space (primarily, in terms of the ionization parameter and local cloud density) as a function of the strengths of Fe II and Ca II emission., Comment: 21 pages, 12 figures, 1 table; accepted for publication in ApJ

In this paper we derive a line tension model for dislocations in 3d starting from a geometrically nonlinear elastic energy with quadratic growth. In the asymptotic analysis, as the amplitude of the Burgers vectors (proportional to the lattice spacing) tends to zero, we show that the elastic energy linearises and the line tension energy density, up to an overall constant rotation, is identified by the linearised cell problem formula given in [17]., Comment: 52 pages, 2 figures

We present near-infrared spectroscopy of the NLS1 galaxy PHL1092 (z=0.394), the strongest FeII emitter ever reported, combined with optical and UV data. We modeled the continuum and the broad emission lines using a power-law plus a black body function and Lorentzian functions, respectively. The strength of the FeII emission was estimated using the latest FeII templates in the literature. We re-estimate the ratio between the FeII complex centered at 4570Ang and the broad component of H-Beta, R_FeII, obtaining a value of 2.58, nearly half of that previously reported (R_FeII=6.2), but still placing PHL1092 among extreme FeII emitters. The FWHM found for low ionization lines are very similar (FWHM~1200km/s), but significantly narrower than those of the Hydrogen lines (FWHM(H-Beta)~1900km/s). Our results suggest that the FeII emission in PHL1092 follows the same trend as in normal FeII emitters, with FeII being formed in the outer portion of the BLR and co-spatial with CaII, and OI, while H-Beta is formed closer to the central source. The flux ratio between the UV lines suggest high densities, log(n_H)~13.0 cm^{-3}, and a low ionization parameter, log(U)~-3.5. The flux excess found in the FeII bump at 9200Ang after the subtraction of the NIR FeII template and its comparison with optical FeII emission suggests that the above physical conditions optimize the efficiency of the ly-Alpha fluorescence process, which was found to be the main excitation mechanism in the FeII production. We discuss the role of PHL1092 in the Eigenvector 1 context., Comment: Accepted for publication at MNRAS

Recent work has shown that it is possible to systematize quasars (QSOs) spectral diversity in 4DE1 parameter space. The spectra contained in most of the surveys have low signal to noise ratio which fed the impression that all QSO\'{s} are spectroscopically similar. Exploration of 4DE1 parameter space gave rise to the concept of two populations of QSOs that present important spectroscopic differences. We aim to quantify broad emission line differences between radio quiet and radio loud sources by exploiting more complete samples of QSO with spectral coverage in Hbeta, MgII and CIV emission lines. We used a high redshift sample (0.35 < z < 1) of strong radio emitter QSOs observations from Calar Alto Observatory in Spain., Comment: Proceedings paper of the IAU symposium "Nuclear Activity in Galaxies Across Cosmic Time" (Ethiopia) accepted for publishing under the Cambridge University Press, eds. M. Povic, P. Marziani, J. Masegosa, H. Netzer, S. H. Negu, and S. B. Tessema

The main sequence offers a method for the systematization of quasar spectral properties. Extreme FeII emitters (or extreme Population A, xA) are believed to be sources accreting matter at very high rates. They are easily identifiable along the quasar main sequence, in large spectroscopic surveys over a broad redshift range. The very high accretion rate makes it possible that massive black holes hosted in xA quasars radiate at a stable, extreme luminosity-to-mass ratio. After reviewing the basic interpretation of the main sequence, we report on the possibility of identifying virial broadening estimators from low-ionization line widths, and provide evidence of the conceptual validity of redshift-independent luminosities based on virial broadening for a known luminosity-to-mass ratio., Comment: 5 pages, 2 figures, contribution presented at the IAU Symposium 356, Nuclear Activity in Galaxies Across Cosmic Time, held in Addis Ababa, Ethiopia, on Oct. 7-11, 2019. M. Povic et al. (Eds.)

Quasars accreting matter at very high rates (known as extreme Population A [xA] quasars, possibly associated with super-Eddington accreting massive black holes) may provide a new class of distance indicators covering cosmic epochs from present day up to less than 1 Gyr from the Big Bang. At a more fundamental level, xA quasars are of special interest in studies of the physics of AGNs and host galaxy evolution. However, their observational properties are largely unknown. xA quasars can be identified in relatively large numbers from major optical surveys over a broad range of redshifts, thanks to selection criteria defined from the systematic changes along the quasars main sequence. It has been possible to build a sample of about 250 quasars at low and intermediate redshift, and larger samples can be easily selected from the latest data releases of the Sloan Digital Sky Survey. A large sample can clarify the main properties of xA quasars which are expected - unlike the general population of quasars - to radiate at an extreme, well defined Eddington ratio with small scatter. As a result of the small scatter in Eddington ratio shown by xA quasars, we propose a method to derive the main cosmological parameters based on redshift-independent "virial luminosity" estimates from measurements of emission line widths, roughly equivalent to the luminosity estimates based from line width in early and late type galaxies. A major issue related to the cosmological application of the xA quasar luminosity estimates from line widths is the identification of proper emission lines whose broadening is predominantly virial over a wide range of redshift and luminosity. We report on preliminary developments using the AlIII 1860 intermediate ionization line and the Hydrogen Balmer line H-beta as virial broadening estimators, and we briefly discuss the perspective of the method based on xA quasars., Comment: 21 pages, 3 figures

The quasar class of extreme Population A (xA) (also known as super-Eddington accreting massive black holes, SEAMBHs) has been hailed as potential distance indicators for cosmology. The aim of this paper is to define tight criteria for their proper identification starting from the main selection criterion $R_{FeII} > 1 $, and to identify potential intruders not meeting the selection criteria, but nonetheless selected as xA because of the coarseness of automatic searches. Inclusion of the spurious xA sources may dramatically increase the dispersion in the Hubble diagram of quasars obtained from virial luminosity estimates. We studied a sample of 32 low-$z$ quasars originally selected from the SDSS DR7 as xA. All of them show moderate-to-strong FeII emission and the wide majority strong absorption features in their spectra are typical of fairly evolved stellar populations. We performed a simultaneous fit of a host galaxy spectrum, AGN continuum, FeII template and emission lines to spectra, using the fitting technique based on ULySS, full spectrum fitting package. For sources in our sample (of spectral types corresponding to relatively low Eddington ratio), we found an overall consistency between narrow components of H$\beta$ and [OIII]$\lambda\lambda$4959, 5007 line shifts and the mean stellar velocity obtained from the host galaxy fit (within $\lesssim |60|$ km/s). Only one source in our sample qualify as xA source. We found high fraction of host galaxy spectrum (in half of the sample even higher then 40\%). When absorption lines are prominent, and the fraction of the host galaxy is high, SSP is mimicking FeII, and that may result in a mistaken identification of FeII spectral features. We have identified several stellar absorption lines that, along with the continuum shape, may lead to an overestimate of $R_{FeII}$, and therefore to the misclassification of sources as xA sources., Comment: 28 pages, 16 figures, accepted for publication in A&A

The interpretation of the main sequence of quasars has become a frontier subject in the last years. The consider the effect of a highly flattened, axially symmetric geometry for the broad-line region (BLR) on the parameters related to the distribution of quasars along their Main Sequence. We utilize the photoionization code CLOUDY to model the BLR, assuming `un-constant' virial factor with a strong dependence on the viewing angle. We show the preliminary results of the analysis to highlight the co-dependence of the Eigenvector 1 parameter, R$_{\mathrm{FeII}}$ on the broad H$\beta$ FWHM (i.e. the line dispersion) and the inclination angle ($\theta$), assuming fixed values for the Eddington ratio ($\mathrm{L_{bol}}/\mathrm{L_{Edd}}$), black hole mass ($\mathrm{M_{BH}}$) and spectral energy distribution (SED) shape. We consider four cases with changing cloud density (n$\rm{_{H}}$) and composition. Understanding the FeII emitting region is crucial as this knowledge can be extended to the use of quasars as distance indicators for Cosmology., Comment: 5 pages, 2 figures, Proceedings of the IAU356: Nuclear Activity in Galaxies across Cosmic Time, Addis Ababa, Ethiopia (7-11 October, 2019)

We address the effect of orientation of the accretion disk plane and the geometry of the broad-line region (BLR) in the context of understanding the distribution of quasars along their Main Sequence. We utilize the photoionization code CLOUDY to model the BLR, incorporating the `un-constant' virial factor. We show the preliminary results of the analysis to highlight the co-dependence of the Eigenvector 1 parameter, $\mathrm{R_{FeII}}$ on the broad H$\beta$ FWHM (i.e. the line dispersion) and the inclination angle ($\theta$), assuming fixed values for the Eddington ratio ($\mathrm{L_{bol}}/\mathrm{L_{Edd}}$), black hole mass ($\mathrm{M_{BH}}$), spectral energy distribution (SED) shape, cloud density (n$\rm{_{H}}$) and composition., Comment: 3 pages, 1 figure, Proceedings of the IAU356: Nuclear Activity in Galaxies across Cosmic Time, Addis Ababa, Ethiopia (7-11 October, 2019)

Quasars accreting matter at very high rates (known as extreme Population A [xA]) may provide a new class of distance indicators covering cosmic epochs from present day up to less than 1 Gyr from the Big Bang. We report on the developments of a method that is based on "virial luminosity" estimates from measurements of emission line widths of xA quasars. The approach is conceptually equivalent to the virial estimates based on early and late type galaxies. The main issues related to the cosmological application of luminosity estimates from xA quasar line widths are the identification of proper emission lines whose broadening is predominantly virial over a wide range of luminosity, and the assessment of the effect of the emitting region orientation with respect to the line of sight. We report on recent developments concerning the use of the AlIII 1860 intermediate ionisation line and of the Hydrogen Balmer line H$\beta$ as "virial broadening estimators.", Comment: Submitted as part of the special issue "Spectral Line Shapes in Astrophysics and Related Topics" to appear in the Contributions of the Astronomical Observatory of Skalnate' Pleso; 13 pages, 5 figures

Context. This is the third study of a series dedicated to the observed parallelism of properties between Galaxy Clusters and Groups(GCGs) and early-type galaxies (ETGs). Aims. Here we investigate the physical origin of the Mass-Radius Relation (MRR). Methods. Having collected literature data on masses and radii for objects going from Globular Clusters (GCs) to ETGs and GCGs, we set up the MR-plane and compare the observed distribution with the MRR predicted by theoretical models both for the monolithic and hierarchical scenarios. Results. We argue that the distributions of stellar systems in the MR-plane is due to complementary mechanisms: (i) on one hand, as shown in paper II, the relation of the virial equilibrium does intersect with a relation that provides the total luminosity as a function of the star formation history; (ii) on the other hand, the locus predicted for the collapse of systems should be convolved with the statistical expectation for the maximum mass of the halos at each comsic epoch. This second aspect provides a natural boundary limit explaining either the curved distribution observed in the MR-plane and the existence of a zone of avoidance. Conclusions. The distribution of stellar systems in the MR-plane is the result of two combined evolution, that of the stellar component and that of the halo component., Comment: 21 pages, 10 figures, 8 tables

We address the effect of the viewing angle of the accretion disk plane and the geometry of the broad-line region (BLR) with the goal of interpreting the distribution of quasars along the main sequence (MS). We utilize photoionization code CLOUDY to model the BLR FeII emission, incorporating the grossly underestimated role of the form factor (f). We recover the dependence of the strength of the FeII emission in the optical (R$_{\rm{FeII}}$) on L$_{\rm{bol}}$/L$_{\rm{Edd}}$ ratio and related observational trends - as a function of the spectral energy distribution (SED) shape, cloud density, composition and intra-cloud dynamics, assumed following prior observational constraints. With this approach, we are now able to explain the diversity of quasars and the change of the quasar properties along the Main Sequence (MS). Our approach also explains the rarity of the highest FeII emitters known as the extreme xA sources and can be used as a predictive tool in future reverberation mapping studies of Type-1 AGNs. This approach further justifies the use of quasars as `cosmological probes'., Comment: 15 pages, 5 figures; Contributions of the Astronomical Observatory Skalnat\'e Pleso (CAOSP) as a Special Issue "Spectral Line Shapes in Astrophysics and Related Topics"

Following the established view of the AGNs inner workings, an AGN is radio-loud (RL) if associated with relativistic ejections emitting a radio synchrotron spectrum (i.e., a jetted AGN). If large samples of optically-selected quasars are considered, AGNs are identified as RL if their Kellermann's radio loudness ratio RK > 10. Our aims are to characterize the optical properties of different classes based on radio-loudness within the quasar main sequence (MS) and to test whether the condition RK > 10 is sufficient for the identification of RL AGNs. A sample of 355 quasars was selected by cross-correlating the FIRST survey with the SDSS DR14 quasar catalog. We classified the optical spectra according to their spectral types along the quasars MS. For each spectral type, we distinguished compact and extended morphology, and three classes of radio-loudness: detected (specific flux ratio in the g band and at 1.4GHz, RK' < 10, RD), intermediate (10 < RK' < 70, RI), and radio loud (RK' > 70). The analysis revealed systematic differences between RD, RI, and RL in each spectral type along the MS. We show that spectral bins that contain the extreme Population A sources have radio power compatible with emission by mechanisms ultimately due to star formation processes. RL sources of Population B are characteristically jetted. Their broad H-beta profiles can be interpreted as due to a binary broad-line region. We suggest that RL Population B sources should be preferential targets for the search of black hole binaries, and present a sample of binary black hole AGN candidates. The validity of the Kellermann's criterion may be dependent on the source location along the quasar MS. The consideration of the MS trends allowed to distinguish between sources whose radio emission mechanisms is jetted from the ones where the mechanism is likely to be fundamentally different.

Context. Galaxy evolution and the effect of environment are most often studied using scaling relations or some regression analyses around some given property. These approaches however do not take into account the complexity of the physics of the galaxies and their diversification. Aims. We here investigate the effect of cluster environment on the evolution of galaxies through multivariate unsupervised classification and phylogenetic analyses applied to two relatively large samples from the WINGS survey, one of cluster members and one of field galaxies (2624 and 1476 objects respectively). Methods. These samples are the largest ones ever analysed with a phylogenetic approach in astrophysics. To be able to use the Maximum Parsimony (cladistics) method, we first performed a pre-clustering in 300 clusters with a hierarchical clustering technique, before applying it to these pre-clusters. All these computations used seven parameters: B-V, log(Re), nV , $\mu$e , H$\beta$ , D4000 , log(M *). Results. We have obtained a tree for the combined samples and do not find different evolutionary paths for cluster and field galaxies. However, the cluster galaxies seem to have accelerated evolution in the sense they are statistically more diversified from a primitive common ancestor. The separate analyses show a hint for a slightly more regular evolution of the variables for the cluster galaxies, which may indicate they are more homogeneous as compared to field galaxies in the sense that the groups of the latter appear to have more specific properties. On the tree for the cluster galaxies, there is a separate branch which gathers rejunevated or stripped-off groups of galaxies. This branch is clearly visible on the colour-magnitude diagram, going back from the red sequence towards the blue one. On this diagram, the distribution and the evolutionary paths of galaxies are strikingly different for the two samples. Globally, we do not find any dominant variable able to explain either the groups or the tree structures. Rather, co-evolution appears everywhere, and could depend itself on environment or mass. Conclusions. This study is another demonstration that unsupervised machine learning is able to go beyond the simple scaling relations by taking into account several properties together. The phylogenetic approach is invaluable to trace the evolutionary scenarii and project them onto any biavariate diagram without any a priori modelling. Our WINGS galaxies are all at low redshift, and we now need to go to higher redshfits to find more primitve galaxies and complete the map of the evolutionary paths of present day galaxies., Comment: Astronomy and Astrophysics - A&A, EDP Sciences, In press

Context. This is the second work dedicated to the observed parallelism between galaxy clusters and early-type galaxies. The focus is on the distribution of these systems in the scaling relations (SRs) observed when effective radii, effective surface brightness, total luminosities and velocity dispersions are mutually correlated. Aims. Using the data of the Illustris simulation we try to speculate on the origin of the observed SRs. Methods. We compare the observational SRs extracted from the database of the WIde-field Nearby Galaxy-cluster Survey (WINGS) with the relevant parameters coming from the Illustris simulations. Then we use the simulated data at different redshift to infer the evolution of the SRs. Results. The comparison demonstrate that galaxy clusters (GCs) at z~0 follow the same log(L)-log(sigma) relation of early-type galaxies (ETGs) and that both in the log(Ie)-log(Re) and log(Re)-log(M*) planes the distribution of GCs is along the sequence defined by the brightest and massive early-type galaxies (BCGs). The Illustris simulation reproduces the tails of the massive galaxies visible both in the log(Ie)-log(Re) and log(Re)-log(M*) planes, but fail to give the correct estimate of the effective radii of the dwarf galaxies that appear too large and those of GCs that are too small. The evolution of the SRs up to z=4 permits to reveal the complex evolutionary paths of galaxies in the SRs and indicate that the line marking the Zone of Exclusion (ZoE), visible both in the log(Ie)-log(Re) and log(Re)-log(M*) planes, is the trend followed by virialized and passively evolving systems. Conclusions. We speculate that the observed SRs originate from the intersection of the virial theorem and a relation L=L_0 x sigma^beta where the luminosities depend on the star formation history., Comment: 22 pages, 14 figures, 4 tables

We address the effect of orientation of the accretion disk plane and the geometry of the broad-line region (BLR) as part of an effort to understand the distribution of quasars in the optical plane of the quasar main sequence. We utilize the photoionization code CLOUDY to model the BLR incorporating the grossly underestimated form factor ($f$). Treating the aspect of viewing angle appropriately, we confirm the dependence of the $R_{\mathrm{FeII}}$ sequence on Eddington ratio and on the related observational trends - as a function of the SED shape, cloud density and composition, verified from prior observations. Sources with $R_{\mathrm{FeII}}$ in the range 1 -- 2 (about 10\% of all quasars, the so-called extreme Population A [xA] quasars) are explained as sources of high, and possibly extreme Eddington ratio along the $R_{\mathrm{FeII}}$ sequence. This result has important implication for the exploitation of xA sources as distance indicators for Cosmology. $\mathrm{FeII}$ emitters with $R_{\mathrm{FeII}} > 2$ are very rare (<1\% of all type 1 quasars). Our approach also explains the rarity of these highest $\mathrm{FeII}$ emitters as extreme xA sources and constrains the viewing angle ranges with increasing H$\beta$ FWHM., Comment: 9 pages, 4 figures, 1 table; accepted for publication in ApJ

The inter-line comparison between high- and low-ionization emission lines has yielded a wealth of information on the quasar broad line region (BLR) structure and dynamics, including perhaps the earliest unambiguous evidence in favor of a disk + wind structure in radio-quiet quasars. We carried out an analysis of the CIV 1549 and Hbeta line profiles of 28 Hamburg-ESO high luminosity quasars and of 48 low-z, low luminosity sources in order to test whether the high-ionization line CIV 1549 width could be correlated with Hbeta and be used as a virial broadening estimator. We analyze intermediate- to high-S/N, moderate resolution optical and NIR spectra covering the redshifted CIV and H$\beta$ over a broad range of luminosity log L ~ 44 - 48.5 [erg/s] and redshift (0 - 3), following an approach based on the quasar main sequence. The present analysis indicates that the line width of CIV 1549 is not immediately offering a virial broadening estimator equivalent to H$\beta$. At the same time a virialized part of the BLR appears to be preserved even at the highest luminosities. We suggest a correction to FWHM(CIV) for Eddington ratio (using the CIV blueshift as a proxy) and luminosity effects that can be applied over more than four dex in luminosity. Great care should be used in estimating high-L black hole masses from CIV 1549 line width. However, once corrected FWHM(CIV) values are used, a CIV-based scaling law can yield unbiased MBH values with respect to the ones based on H$\beta$ with sample standard deviation ~ 0.3 dex., Comment: 43 pages, 15 Figures, submitted to A&A

We have analyzed the parallelism between the properties of galaxy clusters and early-type galaxies (ETGs) by looking at the similarity between their light profiles. We find that the equivalent luminosity profiles of all these systems in the \vfilt\ band, once normalized to the effective radius \re\ and shifted in surface brightness, can be fitted by the S\'ersic's law \Sers\ and superposed with a small scatter ($\le0.3$ mag). By grouping objects in different classes of luminosity, the average profile of each class slightly deviates from the other only in the inner and outer regions (outside $0.1\leq r/R_e\leq 3$), but the range of values of $n$ remains ample for the members of each class, indicating that objects with similar luminosity have quite different shapes. The "Illustris" simulation reproduces quite well the luminosity profiles of ETGs, with the exception of in the inner and outer regions where feedback from supernovae and active galactic nuclei, wet and dry mergers, are at work. The total mass and luminosity of galaxy clusters as well as their light profiles are not well reproduced. By exploiting simulations we have followed the variation of the effective half-light and half-mass radius of ETGs up to $z=0.8$, noting that progenitors are not necessarily smaller in size than current objects. We have also analyzed the projected dark+baryonic and dark-only mass profiles discovering that after a normalization to the half-mass radius, they can be well superposed and fitted by the S\'ersic's law., Comment: 25 pages, 19 figures

Quasars accreting matter at very high rates (known as extreme Population A [xA] or super-Eddington accreting massive black holes) provide a new class of distance indicators covering cosmic epochs from the present-day Universe up to less than 1 Gyr from the Big Bang. The very high accretion rate makes it possible that massive black holes hosted in xA quasars radiate at a stable, extreme luminosity-to-mass ratio. This in turns translates into stable physical and dynamical conditions of the mildly ionized gas in the quasar low-ionization line emitting region. In this contribution, we analyze the main optical and UV spectral properties of extreme Population A quasars that make them easily identifiable in large spectroscopic surveys at low-z (z < 1) and intermediate-z (2 < z < 2.6), and the physical conditions that are derived for the formation of their emission lines. Ultimately, the analysis supports the possibility of identifying a virial broadening estimator from low-ionization line widths, and the conceptual validity of the redshift-independent luminosity estimates based on virial broadening for a known luminosity-to-mass ratio., Comment: 13 pages, 5 figures. Invited lecture at SPIG 2018, Belgrade. To appear in Atom

Here we present a critical view of practical problems in the analysis of optical Fe II emission lines of Type 1 AGN spectra. Besides the very complex and unclear physical interpretation of the Fe II template shape, there are other issues that might affect the results of Fe II contribution as well (like the S/N, the AGN continuum component modelling, the complex shapes of other broad and narrow emission lines that are blending in the same part of the spectrum, the galactic host stellar component that could lead to mimicking of the Fe II template shape on some parts of spectrum, etc.). In this paper we concentrate particularly on the claims of possibility that in some objects Fe II could be strongly shifted to the red. We examine the effects in the fitting procedure that could artificially lead to such results., Comment: 8 pages, 2 figures, 1 table, presented at conference "Revisiting narrow-line Seyfert 1 galaxies and their place in the Universe"; PoS(NLS1-2018)007

PKS 0438$-$436 at a redshift of $z=2.856$ has been previously recognized as possessing perhaps the most luminous known synchrotron jet. Little is known about this source since the maximum elevation above the horizon is low for the Very Large Array (VLA). We present the first VLA radio image that detects the radio lobes. We use both the 151 MHz luminosity, as a surrogate for the isotropic radio lobe luminosity, and the lobe flux density from the radio image to estimate a long term, time averaged, jet power, $\overline{Q} =1.5\pm 0.7 \times 10^{47} \rm{ergs~s^{-1}}$. We analyze two deep optical spectra with strong broad emission lines and estimate the thermal bolometric luminosity of the accretion flow, $L_{\rm{bol}} = 6.7 \pm 3.0 \times 10^{46} \rm{ergs~s^{-1}}$. The ratio, $\overline{Q}/L_{\rm{bol}} = 3.3 \pm 2.6 $, is at the limit of this empirical metric of jet dominance seen in radio loud quasars and this is the most luminous accretion flow to have this limiting behavior. Despite being a very luminous blazar, it previously had no $\gamma$-ray detections (EGRET, AGILE or FERMI) until December 11 - 13 2016 (54 hours) when FERMI detected a flare that we analyze here. The isotropic apparent luminosity from 100 MeV - 100 GeV rivals the most luminous detected blazar flares (averaged over 18 hours), $\sim 5-6 \times 10^{49} \rm{ergs~s^{-1}}$. The $\gamma$-ray luminosity varies over time by two orders of magnitude, highlighting the extreme role of Doppler abberation and geometric alignment in producing the inverse Compton emission., Comment: To appear in ApJ. Revision required to replace the last figure file with the current version

NGC 1275 is one of the most conspicuous active galactic nuclei (AGN) in the local Universe. The radio jet currently emits a flux density of $\sim 10$ Jy at $\sim 1$ mm wavelengths, down from the historic high of $\sim 65$ Jy in 1980. Yet, the nature of the AGN in NGC 1275 is still controversial. It has been debated whether this is a broad emission line (BEL) Seyfert galaxy, an obscured Seyfert galaxy, a narrow line radio galaxy or a BL-Lac object. We clearly demonstrate a persistent H$\beta$ BEL over the last 35 years with a full width half maximum (FWHM) of 4150 - 6000 km/s. We also find a prominent P$\alpha$ BEL (FWHM $\approx 4770 $ km/s) and a weak CIV BEL (FWHM $\approx 4000 $ km/s), H$\beta$/CIV $\approx 2$. A far UV HST observation during suppressed jet activity reveals a low luminosity continuum. The H$\beta$ BEL luminosity is typical of broad line Seyfert galaxies with similar far UV luminosity. X-ray observations indicate a softer ionizing continuum than expected for a broad line Seyfert galaxy with similar far UV luminosity. This is opposite of the expectation of advection dominated accretion. The AGN continuum appears to be thermal emission from a low luminosity, optically thick, accretion flow with a low Eddington ratio, $\sim 0.0001$. The soft, weak ionizing continuum is consistent with the relatively weak CIV BEL. Evidence that the BEL luminosity is correlated with the jet mm wave luminosity is presented. Apparently, the accretion rate regulates jet power., Comment: To appear in ApJ

The most highly accreting quasars (xA) are of special interest in studies of the physics of AGNs and host galaxy evolution. Quasars accreting at high rates (L/LEdd $\sim$ 1) hold promise for use as 'standard candles': distance indicators detectable at very high redshift. However, their observational properties are largely unknown. A large sample of xA can clarify the main properties of quasars radiating near L/LEdd $\sim$ 1 in the Hb spectral range for redshift $\lt$ 0.8. We use selection criteria derived from 4DE1 studies to identify and analyze spectra for a sample of 334 candidate sources identified from the SDSS DR7 database. The source spectra were chosen to show a ratio RFeII between the FeII emission blend at 4570 \AA\ and Hb, RFeII $\gt$ 1. Composite spectra were analyzed for systematic trends as a function of FeII strength, line width, and [OIII] strength. We introduced tighter constraints on the signal-to-noise ratio (S/N) and RFeII values that allowed us to isolate sources most likely to be extreme accretors. We provide a database of detailed measurements. Analysis of the data allows us to confirm that Hb shows a Lorentzian function with a FWHM of Hb $\lt$ 4000 km s-1. Systematic [OIII] blue shifts, as well as a blueshifted component in Hb are revealed. We interpret the blueshifts as related to the signature of outflowing gas from the quasar central engine. The FWHM of Hb is still affected by the blueshifted emission, however, the effect is non-negligible if the FWHM Hb is used as a 'virial broadening estimator'. A effect of the viewing angle on Hb broadening, deriving a correction for those sources that shows disagreement between virial and concordance cosmology luminosity values. The scatter between concordance cosmology and virial luminosity estimates can be reduced if a correction for orientation effects is included in the FWHM Hb value., Comment: Tables 4 and 5 as well as Figure 2 are online only