Changing-look phenomenon observed now in a growing number of active galaxies challenges our understanding of the accretion process close to a black hole. We propose a simple explanation for the sources where multiple semi-periodic outbursts are observed, and the sources are operating at a few per cent of the Eddington limit. The outburst are caused by the radiation pressure instability operating in the narrow ring between the standard gas-dominated outer disk and the hot optically thin inner Advection-Dominated Accretion Flow. The corresponding limit cycle is responsible for periodic outbursts, and the timescales are much shorter than the standard viscous timescale due to the narrowness of the unstable radial zone. Our toy model gives quantitative predictions and works well for multiple outbursts like those observed in NGC 1566, NGC 4151, NGC 5548 and GSN 069, although the shapes of the outbursts are not yet well modeled, and further development of the model is necessary.

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.

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.

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.

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.

We study the long-term variability in the optical monitoring database of Ark~120, a nearby radio-quiet active galactic nucleus (AGN) at a distance of 143 Mpc (z=0.03271). We compiled the historical archival photometric and spectroscopic data since 1974 and conducted a new two-year monitoring campaign in 2015-2017, resulting in a total temporal baseline over four decades. The long-term variations in the optical continuum exhibit a wave-like pattern and the Hbeta integrated flux series varies with a similar behavior. The broad Hbeta profiles have asymmetric double peaks, which change strongly with time and tend to merge into a single peak during some epochs. The period in the optical continuum determined from various period-search methods is about 20 yr and the estimated false alarm probability with null hypothesis simulations is about 1*10^-3. The overall variations of the broad Hbeta profiles also follow the same period. However, the present database only covers two cycles of the suggested period, which strongly encourages continued monitoring to track more cycles and confirm the periodicity. Nevertheless, in light of the possible periodicity and the complicated Hbeta profile, Ark~120 is one candidate of the nearest radio-quiet AGNs with possible periodic variability, and it is thereby a potential candidate host for a sub-parsec supermassive black hole binary.

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."

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.

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.

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.

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.

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.

Narrow-line Seyfert 1s (NLSy1s) are an ill-defined class. Work done over the past 20 years as well as recent analyses show a continuity in properties (e.g., Balmer line profiles, blueshifts of high-ionization lines) between sources with FWHM above and below 2000 km/s, the defining boundary of NLSy1s. This finding alone suggests that comparisons between samples of NLSy1s and rest of broad-line AGNs are most likely biased. NLSy1s can be properly contextualized by their location on the quasar main sequence originally defined by Sulentic et al 2000. At one end, NLSy1s encompass sources with strong FeII emission and associated with high Eddington ratio that hold the promise of becoming useful distance indicators; at the other end, at least some of them are sources with broad profiles seen face-on. Any rigid FWHM limit gives rise to some physical ambiguity, as the FWHM of low-ionization lines depends in a complex way on mass, Eddington ratio, orientation, and luminosity. In addition, if the scaling derived from luminosity and virial dynamics applies to the broad line regions, NLSy1s at luminosity higher than 1E47 erg/s become physically impossible. Therefore, in a broader context, a proper subdivision of two distinct classes of AGNs and quasars may be achieved by the distinction between Pop. A and B with boundary at = 4000 km/s in samples at z < 1, or on the basis of spectrophotometric properties which may ultimately be related to differences in accretion modes if high-luminosity quasars are considered.

Here we investigate the connection of broad emission line shapes and continuum light curve variability time scales of type-1 Active Galactic Nuclei (AGN). We developed a new model to describe optical broad emission lines as an accretion disk model of a line profile with additional ring emission. We connect ring radii with orbital time scales derived from optical light curves, and using Kepler's third law, we calculate mass of central supermassive black hole (SMBH). The obtained results for central black hole masses are in a good agreement with {other methods. This indicates that the variability time scales of AGN may not be stochastic, but rather connected to the orbital time scales which depend on the central SMBH mass.

Active galactic nuclei (AGN) are one of the most luminous objects in the Universe, emitting powerful continuum and line emission across all wavelength bands. They represent an important link in the investigations of the galaxy evolution and cosmology. The resolving of the AGN inner structure is still a difficult task with current instruments, therefore the spectroscopy and spectropolarimetry are crucial tools to investigate these objects and their components, such as the properties of the supermassive black hole, the broad line region, and the dusty torus. In this review, we present the results of the project "Astrophysical spectroscopy of extragalactic objects", from the observations, data processing and analysis, to the modelling of different regions in AGN.

The last 25 years saw a major step forward in the analysis of optical and UV spectroscopic data of large quasar samples. Multivariate statistical approaches have led to the definition of systematic trends in observational properties that are the basis of physical and dynamical modeling of quasar structure. We discuss the empirical correlates of the so-called "main sequence" associated with the quasar Eigenvector 1, its governing physical parameters and several implications on our view of the quasar structure, as well as some luminosity effects associated with the virialized component of the line emitting regions. We also briefly discuss quasars in a segment of the main sequence that includes the strongest FeII emitters. These sources show a small dispersion around a well-defined Eddington ratio value, a property which makes them potential Eddington standard candles.

Many active galactic nuclei (AGN) show strong variability of the optical continuum. Since the line flux, profile shapes and intensity ratios are changing, we analyze the variability patterns and possible periodicity of Type 1 AGN NGC 5548, using the Eigenvector 1 (EV1) diagram in different variability states, taking advantage of very long term monitoring campaign data. The preliminary results suggest that NGC 5548 - a highly variable object that over several decades has shown large amplitude continuum fluctuations and flaring behavior - remains Pop B. This means that the range in Eddington ratio, even when the source is in a bright state, remains consistent with the value of the low accreting Pop B. We inspected EV 1 parameters of a single object though long term monitoring, assuming an inclination and black hole mass to be constant during the observational time. Our results imply that the main driver for the variations along the EV 1 diagram could be dimensionless accretion rate. If so, then it appears that the source never crossed the boundary for structural changes, indicatively placed at $L/L_{Edd} \sim$ 0.2.

Here we present the evidence for periodicity of an optical emission detected in several AGN. Significant periodicity is found in light curves and radial velocity curves. We discuss possible mechanisms that could produce such periodic variability and their implications. The results are consistent with possible detection of the orbital motion in proximity of the AGN central supermassive black holes.

We present an analysis of 43 years (1972 to 2015) of spectroscopic observations of the Seyfert 1 galaxy NGC 5548. This includes 12 years of new unpublished observations (2003 to 2015). We compiled about 1600 H$\beta$ spectra and analyzed the long-term spectral variations of the 5100 \AA\ continuum and the H$\beta$ line. Our analysis is based on standard procedures including the Lomb-Scargle method, which is known to be rather limited to such heterogeneous data sets, and new method developed specifically for this project that is more robust and reveals a $\sim$5700 day periodicity in the continuum light curve, the H$\beta$ light curve and the radial velocity curve of the red wing of the H$\beta$ line. The data are consistent with orbital motion inside the broad emission line region of the source. We discuss several possible mechanisms that can explain this periodicity, including orbiting dusty and dust-free clouds, a binary black hole system, tidal disruption events, and the effect of an orbiting star periodically passing through an accretion disc.

The detection of gravitationally redshifted optical emission lines has been reported just for a few active galaxies. In this paper we give a short overview of studies that analyzed or exploited the detection of the gravitational redshift in optical AGN spectra. In addition, we tested the consistency of gravitational redshift as the physical origin of the redward shifts observed in their spectra using a sample of $\approx$ 50 Hamburg-ESO intermediate to high redshift quasars that are among the most luminous quasars known ($10^{47} \lesssim L \lesssim 10^{48}$ erg/s), and are expected to host very massive black holes. To this aim we modeled the line profile with accretion disk models around a black hole.

We investigated the efficiency of estimating characteristics of stellar populations (SP) and Active Galactic Nuclei (AGN) emission using ULySS code. To analyse simultaneously AGN and SP components in the integrated spectrum of Type 2 active galaxies, we modelled the featureless continuum (FC) and emission lines, and we used PEGASE.HR stellar population models provided by ULySS. In order to validate the method, we simulated over 7000 integrated spectra of Seyfert 2 galaxies. Spectra were generated using different characteristics of the featureless AGN continuum, signal-to-noise ratio (SNR), spectral ranges, properties of emission lines and single stellar population (SSP) model whose initial mass function (IMF) and abundance pattern is similar to the solar neighbourhood. Simulated spectra were fitted with ULySS to evaluate the ability of the method to extract SP and AGN properties. We found that the analysis with ULySS can efficiently restore the characteristics of SP in spectra of Seyfert 2 AGNs, where signal-to-noise ratio is higher than 20, and where SP contributes with more than 10% to the total flux. Degeneracies between AGN and SP parameters increase with increasing the AGN continuum fraction, which points out the importance of simultaneous fitting of the FC and SP contributions.

One of the most intriguing scenarios proposed to explain how active galactic nuclei are triggered involves the existence of a supermassive binary black hole system in their cores. Here we present an observational evidence for the first spectroscopically resolved sub-parsec orbit of a such system in the core of Seyfert galaxy NGC 4151. Using a method similar to those typically applied for spectroscopic binary stars we obtained radial velocity curves of the supermassive binary system, from which we calculated orbital elements and made estimates about the masses of components. Our analysis shows that periodic variations in the light and radial velocity curves can be accounted for an eccentric, sub-parsec Keplerian orbit of a 15.9-year period. The flux maximum in the lightcurve correspond to the approaching phase of a secondary component towards the observer. According to the obtained results we speculate that the periodic variations in the observed H{\alpha} line shape and flux are due to shock waves generated by the supersonic motion of the components through the surrounding medium. Given the large observational effort needed to reveal this spectroscopically resolved binary orbital motion we suggest that many such systems may exist in similar objects even if they are hard to find. Detecting more of them will provide us with insight into black hole mass growth process.