Optical Fe II emission is a strong feature of quasar spectra originating in the broad-line regions (BLRs). The difficulty in understanding the complex Fe II pseudo-continuum has led us in search of 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 species with our photoionization models using CLOUDY. Taking into account the effect of dust into our modelling, we retrieve the log U - log n parameter space as a function of the strengths of Fe II and Ca II emission. Incorporating a broad range of metallicity and extended cloud columns we reinstate this correlation and explain the emission from the strong Fe II emitting Narrow-Line Seyfert 1 galaxies (NLS1s).

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.

Reverberation mapping technique is an important milestone in AGN demographics, kinematics and the structure of the Broad Line Region (BLR) based on the time-delay response between the continuum and emission line. The time delay is directly related to the size of BLR which is related to the continuum luminosity of the source, producing the well-known Radius-Luminosity (RL) relation. The majority of sources have been monitored for their H$\beta$ emission line in low redshift sources ($z<0.1$), while there are some attempts using the MgII line for higher redshifts. We present a recent MgII monitoring for the quasar CTS C30.10 ($z=0.90$) observed with the 10-meter SALT, for which RL scaling based on MgII holds within measurement and time-delay uncertainties. One of the most important advantages of reverberation mapping technique is the independent determination to the distant source, and considering the large range of redshifts and luminosities found in AGNs, their application to cosmology is promising. However, recently it has been found that highly accreting sources show the time delays shorter than expected from RL relation. We propose correction for this effect using a sample of 117 H$\beta$ reverberating-mapped AGNs with $0.02<z<0.9$, recovering the low scatter along with the relation. We are able to determine the cosmological constants, $\Omega_m$ and $\Omega_\Lambda$ within 2$\sigma$ confidence level. We present the first steps in modelling of light curves for H$\beta$ and MgII and discuss the quasar selection in the context of photometric reverberation mapping with Large Synoptic Survey Telescope (LSST). With the onset of the LSST era, we expect a huge rise in the overall quasar counts and redshift range covered ($z\lesssim 7.0$), which will provide a better constraint of AGN properties with cosmological purposes.

In the near future, new surveys promise a significant increase in the number of quasars (QSO) at large redshifts. This will help to constrain the dark energy models using quasars. The Large Synoptic Survey Telescope (LSST) will cover over 10 million QSO in six photometric bands during its 10-year run. QSO will be monitored and subsequently analyzed using the photometric reverberation mapping (RM) technique. In low-redshift quasars, the combination of reverberation-mapped and spectroscopic results have provided important progress. However, there are still some facts which have to be taken into account for future results. It has been found that super-Eddington sources show time delays shorter than the expected from the well-known Radius-Luminosity ($R\mathrm{_{H\beta}}-L_{5100}$) relation. Using a sample of 117 H$\beta$ reverberation-mapped AGN with $0.02<z<0.9$, we propose a correction by the accretion rate effect recovering the classical $R\mathrm{_{H\beta}}-L_{5100}$ relation. We determined the cosmological constants, which are in agreement with $\Lambda$-Cold Dark Matter model within 2$\sigma$ confidence level, which is still not suitable for testing possible departures from the standard model. Upcoming LSST data will decrease the uncertainties in the dark energy determination using reverberation-mapped sources, particularly at high redshifts. We show the first steps in the modeling of the expected light curves for H$\beta$ and MgII.

Quasars have been proposed as a new class of standard candles analogous to Supernovae, as their large redshift range and high luminosities make them excellent candidates. Reverberation mapping (RM) method enables to estimate the distance to the source from the time delay measurement of the emission lines with respect to the continuum, since the time delay depends on the absolute luminosity of the source. This radius-luminosity (RL) relation exhibits a low scatter and it was thus proposed for cosmological purposes. However, in recent years the increase in the studied sample, and in particular the inclusion of highly accreting QSO has increased the dispersion in the RL relation, with many objects showing time delays shorter than the expected. Using H$\beta$ RM measurements for 117 sources with $0.002<z<0.9$ and $41.5<\mathrm{log}~L_{5100}<45.9$, we find a correction for the time delay based on the dimensionless accretion rate considering the virial factor anti-correlated with the FWHM of H$\beta$. This correction decreases the scattering in the accretion parameters compared with the typical value used, which is directly reflected by suppressing the RL relation dispersion. We also confirm the anti-correlation between the excess of variability and the accretion parameters. With this correction we are able to build the Hubble diagram and estimate the cosmological constants $\Omega_m$ and $\Omega_\Lambda$, which are consistent with the standard cosmological model at 2$\sigma$ confidence level. Therefore, RM results can be used to constrain cosmological models in the future.

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.

We present a multiwavelength spectral study of the quasar HE0359-3959, which has been identified as an extreme radiating source at intermediate redshift (z=1.5209). Along the spectral range, the different ionic species give information about the substructures in the broad line region. The presence of a powerful outflow with an extreme blueshifted velocity of $\sim$-6000$\pm$500 km s$^{-1}$ is shown in the CIV{\lambda}1549 emission line. A prominent blueshifted component is also associated with the 1900{\AA} blend, resembling the one observed in CIV{\lambda}1549. We detect a strong contribution of very the low-ionization lines, FeII and Near-Infrared CaII triplet. We find that the physical conditions for the low, intermediate and high-ionization emission lines are different, which indicate that the emission lines are emitted in different zones of the broad line region. The asymmetries shown by the profiles reveal different forces over emitter zones. The high-ionization region is strongly dominated by radiation forces, which also affect the low and intermediate-ionization emitter region, commonly governed by virial motions. These results support the idea that highly radiating sources host a slim disk.

Broad emission lines in quasars enable us to "resolve" structure and kinematics of the broad line emitting region (BLR) thought to in- volve an accretion disk feeding a supermassive black hole. Interpretation of broad line measures within the 4DE1 formalism simplifies the apparent confusion among such data by contrasting and unifying properties of so-called high and low accreting Population A and B sources. H{\beta} serves as an estimator of black hole mass, Eddington ratio and source rest frame, the latter a valuable input for Civ{\lambda}1549 studies which allow us to isolate the blueshifted wind component. Optical and HST-UV spectra yield H{\beta} and Civ{\lambda}1549 spectra for low-luminosity sources while VLT-ISAAC and FORS and TNG-LRS provide spectra for high Luminosity sources. New high S/N data for Civ in high-luminosity quasars are presented here for comparison with the other previously published data. Comparison of H{\beta} and Civ{\lambda}1549 profile widths/shifts indicates that much of the emission from the two lines arise in regions with different structure and kinematics. Covering a wide range of luminosity and redshift shows evidence for a correlation between Civ{\lambda}1549 blueshift and source Eddington ratio, with a weaker trend with source luminosity (similar amplitude outflows are seen over 4 of the 5 dex luminosity range in our combined samples). At low luminosity (z < 0.7) only Population A sources show evidence for a significant outflow while at high luminosity the outflow signature begins to appear in Population B quasars as well.

We present an unprecedented spectroscopic survey of the CaII triplet + OI for a sample of 14 luminous ($-$26 $\gtrsim$ M$_V$ $\gtrsim$ $-$29), intermediate redshift (0.85 $\lesssim$ $z$ $\lesssim$ 1.65) quasars. The ISAAC spectrometer at ESO VLT allowed us to cover the CaII NIR spectral region redshifted into the H and K windows. We describe in detail our data analysis which enabled us to detect CaII triplet emission in all 14 sources (with the possible exception of HE0048-2804) and to retrieve accurate line widths and fluxes of the triplet and OI $\lambda$8446. The new measurements show trends consistent with previous lower $z$ observations, indicating that CaII and optical FeII emission are probably closely related. The ratio between the CaII triplet and the optical FeII blend at $\lambda$4570 $\AA$ is apparently systematically larger in our intermediate redshift sample relative to a low-$z$ control sample. Even if this result needs a larger sample for adequate interpretation, higher CaII/optical FeII should be associated with recent episodes of star formation in the intermediate redshift quasars and, at least in part, explain an apparent correlation of CaII triplet equivalent width with $z$ and $L$. The CaII triplet measures yield significant constraints on the emitting region density and ionization parameter, implying CaII triplet emission from log n$_H$ $\gtrsim$ 11 [cm$^{-3}$] and ionization parameter log $U$ $\lesssim$ 1.5. Line width and intensity ratios suggest properties consistent with emission from the outer part of a high density broad line region (a line emitting accretion disk?).

Highly accreting quasars are characterized by distinguishing properties in the 4D eigenvector 1 parameter space that make them easily recognizable over a broad range range of redshift and luminosity. The 4D eigenvector 1 approach allows us to define selection criteria that go beyond the restriction to Narrow Line Seyfert 1s identified at low redshift. These criteria are probably able to isolate sources with a defined physical structure i.e., a geometrically thick, optically thick advection-dominated accretion disk (a "slim" disk). We stress that the importance of highly accreting quasars goes beyond the understanding of the details of their physics: their Eddington ratio is expected to saturate toward values of order unity, making them possible cosmological probes.

In order to investigate where and how low ionization lines are emitted in quasars we are studying a new collection of spectra of the CaII triplet at $\lambda$8498, $\lambda$8542, $\lambda$8662 observed with the Very Large Telescope (VLT) using the Infrared Spectrometer And Array Camera (ISAAC). Our sample involves luminous quasars at intermediate redshift for which CaII observations are almost nonexistent. We fit the CaII triplet and the OI $\lambda$8446 line using the H$\beta$ profile as a model. We derive constraints on the line emitting region from the relative strength of the CaII triplet, OI $\lambda$8446 and H$\beta$.