Search

The Eigenvector 1 schema, or the main sequence of quasars, was introduced as an analogous scheme to the HR diagram that would allow us to understand the more complex, extended sources - active galactic nuclei (AGNs) that harbor accreting supermassive black holes. The study has spanned more than three decades and has advanced our knowledge of the diversity of Type-1 AGNs from both observational and theoretical aspects. The quasar main sequence, in its simplest form, is the plane between the FWHM of the broad H$\beta$ emission line and the strength of the optical FeII emission to the H$\beta$. While the former allows the estimation of the black hole mass, the latter enables direct measurement of the metal content and traces the accretion rate of the AGN. Together, they allow us to track the evolution of AGN in terms of the activity of the central nuclei, its effect on the line-emitting regions surrounding the AGN, and their diversity making them suitable distance indicators to study the expansion of our Universe. This mini-review aims to provide (i) a brief history leading up to the present day in the study of the quasar main sequence, (ii) introduce us to the many possibilities to study AGNs with the main sequence as a guiding tool, and (iii) highlight some recent, exciting lines of researches at the frontier of this ever-growing field., Comment: Mini Review Article based on the invited talk in Frontiers in Astronomy and Space Sciences: A Decade of Discovery and Advancement - 10th Anniversary Conference. This version matches the published version (19-09-2024)

The gamma-ray emitting narrow-line Seyfert 1 galaxies are a unique class of objects that launch powerful jets from relatively lower-mass black hole systems compared to the Blazars. However, the black hole masses estimated from the total flux spectrum suffer from the projection effect, making the mass measurement highly uncertain. The polarized spectrum provides a unique view of the central engine through scattered light. We performed spectro-polarimetric observations of the gamma-ray emitting narrow-line Seyfert 1 galaxy 1H0323+342 using SPOL/MMT. The degree of polarization and polarization angle is 0.122 $\pm$ 0.040 % and 142 $\pm$ 9 degrees, while the H$\alpha$ line is polarized at 0.265 $\pm$ 0.280 %. We decomposed the total flux spectrum and estimated broad H$\alpha$ FWHM of 1015 km/s. The polarized flux spectrum shows a broadening similar to the total flux spectrum, with a broadening ratio of 1.22. The Monte Carlo radiative transfer code `STOKES' applied to the data provides the best fit for a small viewing angle of 9-24 degrees and a small optical depth ratio between the polar and the equatorial scatters. A thick BLR with significant scale height can explain a similar broadening of the polarized spectrum compared to the total flux spectrum with a small viewing angle., Comment: Accepted for publication in MNRAS (Jul. 08th, 2024). 12 pages, 6 figures

Reverberation mapping accurately determines virial black hole masses only for redshifts $z <$ 0.2 by utilizing the relationship between the H$\beta$ broad-line region (BLR) size and the 5100 Angstroms continuum luminosity established with $\sim 200$ active galactic nuclei (AGN). For quasars at $z \sim 2-3$ determining the BLR size is time-consuming and limited by seasonal gaps, requiring e.g., $\sim$ 20 years of monitoring of the CIV emission lines. In this work, we demonstrate that an efficient alternative is to use a continuum size-luminosity relation, which can be obtained up to 150 times faster than BLR sizes using photometric reverberation mapping (PRM). We outline the method and its feasibility based on simulations and propose an observational strategy that can be carried out with meter-class telescopes. In particular, we focus on the ESO La Silla 2.2 meter telescope as it is suitable for an efficient PRM campaign. These observations will provide the scaling factor between the accretion disk and the BLR size (for CIV-1350 Angstroms), which is crucial for estimating the masses of black holes at higher redshifts ($z \gtrsim 2-3$)., Comment: 12 pages, 5 figures, accepted to the Astrophysical Journal Letters

We use the spectroscopic data collected by the Magellanic Quasars Survey (MQS) as well as the photometric V- and I-band data from the Optical Gravitational Lensing Experiment (OGLE) to measure the physical parameters for active galactic nuclei (AGNs) located behind the Magellanic Clouds. The flux-uncalibrated MQS spectra were obtained with the 4-m Anglo-Australian Telescope and the AAOmega spectroscope (R=1300) in a typical ~1.5-hour visit. They span a spectral range of 3700-8500 Angstroms and have S/N ratios in a range of 3-300. We report the discovery and observational properties of 161 AGNs in this footprint, which expands the total number of spectroscopically confirmed AGNs by MQS to 919. After converting the OGLE mean magnitudes to the monochromatic luminosities at 5100 Angstroms, 3000 Angstroms, and 1350 Angstroms, we reliably measured the black hole masses for 165 out of 919 AGNs. The remaining physical parameters we provide are the bolometric luminosities and the Eddington ratios. A fraction of these AGNs have been observed by the OGLE survey since 1997 (all of them since 2001), enabling studies of correlations between their variability and physical parameters., Comment: 15 pages, 8 figures, 8 tables, accepted for publication in The Astrophysical Journal Supplement Series

Constraining the physical conditions of the ionized media in the vicinity of an active supermassive black hole (SMBH) is crucial to understanding how these complex systems operate. Metal emission lines such as iron (Fe) are useful probes to trace the gaseous media's abundance, activity, and evolution in these accreting systems. Among these, the FeII emission has been the focus of many prior studies to investigate the energetics, kinematics, and composition of the broad-emission line region (BELR) from where these emission lines are produced. In this work, we present the first simultaneous FeII modeling in the optical and near-infrared (NIR) regions. We use CLOUDY photoionization code to simulate both spectral regions in the wavelength interval 4000-12000 Angstroms. We compare our model predictions with the observed line flux ratios for IZw1 - a prototypical strong FeII-emitting active galactic nuclei (AGN). This allows putting constraints on the BLR cloud density and metal content that is optimal for the production of the FeII emission, which can be extended to IZw1-like sources, by examining a broad parameter space. We demonstrate the salient and distinct features of the FeII pseudo-continuum in the optical and NIR, giving special attention to the effect of micro-turbulence on the intensity of the FeII emission., Comment: 17 pages, 6 figures, 1 table, accepted in MDPI Physics

The investigation of emission line regions within active galaxies (AGNs) has a rich and extensive history, now extending to the use of AGNs and quasars as "standardizable" cosmological indicators, shedding light on the evolution of our universe. As we enter the era of advanced observatories, such as the successful launch of JWST and the forthcoming Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST), the landscape of AGN exploration across cosmic epochs is poised for exciting advancements. In this work, we delve into recent developments in AGN variability research, anticipating the substantial influx of data facilitated by LSST. The article highlights recent strides made by the AGN Polish Consortium in their contributions to LSST. The piece emphasizes the role of quasars in cosmology, dissecting the intricacies of their calibration as standard candles. The primary focus centers on the relationship between the broad-line region size and luminosity, showcasing recent breakthroughs that enhance our comprehension of this correlation. These breakthroughs encompass a range of perspectives, including spectroscopic analyses, photoionization modeling, and collaborative investigations with other cosmological tools. The study further touches on select studies, underlining how the synergy of theoretical insights and advancements in observational capabilities has yielded deeper insights into these captivating cosmic entities., Comment: 34 pages, 5 figures, accepted for publication as a review in MDPI/Universe

Previously, we demonstrated that MgII and CIV reverberation-mapped quasars (RM QSOs) are standardizable and that the cosmological parameters inferred using the broad-line region radius-luminosity (R-L) relation are consistent with those determined from better-established cosmological probes. With more data expected from ongoing and future spectroscopic and photometric surveys, it is imperative to examine how new QSO data sets of varied quality, with their own specific luminosity and time-delay distributions, can be best used to determine more restrictive cosmological parameter constraints. In this study, we test the effect of adding 25 OzDES MgII RM QSOs as well as 25 lower-quality SDSS RM CIV QSOs, which increases the previous sample of RM QSOs by $\sim 36\%$. Although cosmological parameter constraints become tighter for some cosmological models after adding these new QSOs, the new combined data sets have increased differences between R-L parameter values obtained in different cosmological models and thus a lower standardizability for the larger MgII + CIV compilation. Different time-delay methodologies, particularly the ICCF and CREAM methods used for inferring time delays of SDSS RM QSOs, slightly affect cosmological and R-L relation parameter values, however, the effect is negligible for (smaller) compilations of robust time-delay detections. Our analysis indicates that increasing the sample size is not sufficient for tightening cosmological constraints and a quality cut is necessary to obtain a standardizable RM QSO sample., Comment: 25 pages, 14 figures, 6 tables; accepted by the MNRAS Main Journal

NGC 5548 has been hailed as an archetypical type-1 active galactic nuclei (AGN) and serves as a valuable laboratory to study the long-term variation of its broad-line region (BLR). In this work, we re-affirm our finding on the connection between the continuum variability in the optical regime and the corresponding H$\beta$ response to it, in order to realize the increase, albeit with a gradual saturation, in the H$\beta$ emitting luminosity with increasing AGN continuum. This effect is also known as the Pronik-Chuvaev effect after the authors who first demonstrated this effect using long-term monitoring of another well-studied Type-1 AGN - Mrk 6. We employ a homogeneous, multi-component spectral fitting procedure over a broad range of spectral epochs that is then used to create the continuum and H$\beta$ light curves. We focus on the epoch range 48636-49686 MJD, different from our previous analysis. We again notice a clear signature of shallowing in the trend between the H$\beta$ and the continuum luminosities. We attempt to recover this H$\beta$ emission trend as a response to a significant continuum flux increase using CLOUDY photoionization simulations and employ a suitable broad-band spectral energy distribution for this source. We explore the wide range in the physical parameters space for modelling the H$\beta$ emission from the BLR appropriate for this source. We employ a constant density, single cloud model approach in this study and successfully recover the observed shallowing of the H$\beta$ emission with respect to the rising AGN continuum and provide constraints on the local BLR densities and the location of the H$\beta$ emitting BLR which agrees with the H$\beta$ time-lags reported from the long-term reverberation mapping monitoring. On the contrary, we do not find a significant breathing effect in the location of the H$\beta$ line-emitting BLR for this epoch in NGC 5548., Comment: 5 pages, 7 figures, version matches the published version in Proceedings of the XLV Annual Meeting of the Brazilian Astronomical Society. arXiv admin note: substantial text overlap with arXiv:2111.05378

Broad lines in active galaxies are primarily emitted by photoionization processes that are driven by the incident continuum arising from a complex geometrical structure circumscribing the supermassive black hole. A model of the broad-band spectral energy distribution (SED) effective in ionizing the gas-rich broad line emitting region (BLR) is needed to understand the various radiative processes that eventually lead to the emission of emission lines from diverse physical conditions. Photoionization codes are a useful tool to investigate two aspects - the importance of the shape of the SED, and the physical conditions in the BLR. In this work, we focus on the anisotropy of continuum radiation from the very centre a direct consequence of the development of a funnel-like structure at regions very close to the black hole. Accounting for the diversity of Type-1 active galactic nuclei (AGNs) in the context of the main sequence of quasars permits us to locate the super Eddington sources along the sequence and constrain the physical conditions of their line-emitting BLR., Comment: 5 pages, 3 figures, version matches the published version in Proceedings of the XLV Annual Meeting of the Brazilian Astronomical Society. arXiv admin note: text overlap with arXiv:2210.15041

Double-peaked profiles associated with the broad-line region (BLR) of active galactic nuclei (AGNs) are regarded as the clearest evidence of the presence of an accretion disk. They are most commonly detected by means of optical spectroscopy in the Balmer lines and in the Mg II $\lambda$2798 ultraviolet line. Here, we report the first unambiguous detection of a double-peak broad emission line associated with the O I $\lambda$11297 emission line in the near-infrared (NIR) in the local Seyfert 1 galaxy III Zw 002. Additionally, we detect simultaneously in the spectrum the double-peak emission in the Pa$\alpha$ line and very likely in the He II $\lambda$10830. This is the first time that several broad double-peaked NIR emission lines have been detected simultaneously. The double-peaked profiles are fit using a disk-based model, with an additional Gaussian component attributed to nondisk clouds, which represents the classical BLR. Our results obtained from the fits reveal important parameters, such as disk inclination and geometry. From the double-peaked profile fits, we suggest that the BLR in III Zw 002 has a disk-like geometry, as it extends up to the outer edge of the BLR., Comment: 8 pages, 3 figures, 1 table, published in ApJ Letters

A landmark of accretion processes in active galactic nuclei (AGN) is the continuum originating from a complex structure, i.e., an accretion disk and a corona around a supermassive black hole. Modelling the broad-band spectral energy distribution (SED) effectively ionizing the gas-rich broad emission line region (BLR) is key to understanding the various radiative processes at play and their importance that eventually leads to the emission from diverse physical conditions. Photoionization codes are a useful tool to investigate two aspects, the importance of the shape of the spectral energy distribution, and the physical conditions in the broad emission line region. In this work, we critically review long-standing issues pertaining to the spectral energy distribution shape and the anisotropic continuum radiation from the central regions around the accreting supermassive black holes (a few 10-100 gravitational radii), with a focus on black holes accreting at high rates, possibly much above the Eddington limit. The anisotropic emission is a direct consequence of the development of a geometrically and optically thick structure at regions very close to the black hole due to a marked increase in the accretion rates. The analysis presented in this paper took advantage of the look at the diversity of the type-1 active galactic nuclei provided by the main sequence of quasars. The main sequence permitted us to assess the importance of the Eddington ratio and hence to locate the super Eddington sources in observational parameter space, as well as to constrain the distinctive physical conditions of their line-emitting BLR. This feat is posing the basis for the exploitation of quasars as cosmological distance indicators, hopefully allowing us to use the fascinating super Eddington quasars up to unprecedented distances., Comment: Identical to the published version in Frontiers in Astronomy and Space Sciences (09-05-2023). 31 pages, 8 figures, 2 tables, based on the invited talk presented at the SPIG - 31st Summer School and International Symposium on the Physics of Ionized Gases, held between 05th - 09th September 2022

As Setti & Woltjer noted back in 1973, quasars could be used to construct the Hubble diagram but the actual application was not that straightforward. It took years to implement the idea successfully. Most of the ways to use quasars for cosmology now require an advanced understanding of their structure, step by step. We briefly review this progress, with unavoidable personal bias, and concentrate on bright unobscured sources. We will mention the problem of the gas flow character close to the innermost stable circular orbit close to the black hole, discussed 50 years ago, which later led to the development of the slim disk scenario, but was recently revived in the context of Magnetically Arrested Disks (MAD) and Standard and Normal Evolution (SANE) disk models. We also discuss the hot/warm corona issue, which is still under discussion and complicates the analysis of the X-ray reflection. We present the scenario of the formation of the low ionization part of the Broad Line Region as a mostly failed wind powered by radiation pressure acting on dust (FRADO - Failed Radiatively Driven Dusty Outflow model). Next, we discuss the cosmological constraints that are currently achievable with quasars, mostly concentrating on light echo methods (continuum time delays and spectral line time delays with respect to the continuum) which are (or should be) incorporating the progress mentioned above. Finally, we briefly mention future prospects in this direction., Comment: Astrophysics and Space Science (in press)

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will detect an unprecedentedly large sample of actively accreting supermassive black holes with typical accretion disk (AD) sizes of a few light days. This brings us to face challenges in the reverberation mapping (RM) measurement of AD sizes in active galactic nuclei (AGNs) using interband continuum delays. We examine the effect of LSST cadence strategies on AD RM using our metric AGNTimeLagMetric. It accounts for redshift, cadence, the magnitude limit, and magnitude corrections for dust extinction. Running our metric on different LSST cadence strategies, we produce an atlas of the performance estimations for LSST photometric RM measurements. We provide an upper limit on the estimated number of quasars for which the AD time lag can be computed within 01000 sources in each Deep Drilling field (DDF, 10 sq. deg) in any filter, with the redshift distribution of these sources peaking at z~1. We find the LSST observation strategies with a good cadence (~ 5 days) and a long cumulative season (~9 yr), as proposed for LSST DDF, are favored for the AD size measurement. We create synthetic LSST light curves for the most suitable DDF cadences and determine RM time lags to demonstrate the impact of the best cadences based on the proposed metric., Comment: accepted for publication in the Astrophysical Journal Focus Issue on Rubin LSST cadence and survey strategy

Context. Continuum time delays from accretion disks in active galactic nuclei (AGN) has been proposed long time ago as a tool for measuring distances to the monitored sources. However, the method faces serious problems as a number of effects must be taken into account, including the contribution from the Broad Line Region (BLR). Aims. In this paper, we model the expected time delays when both the disk reprocessing of the incident X-ray flux and further reprocessing by the BLR are included, with the aim to see if the two effects can be disentangled. Methods. We use simple response function for the accretion disk, without relativistic effects, and we use a parametric description to account for the BLR contribution. We include only scattering of the disk emission by the BLR inter-cloud medium. We also use artificial lightcurves with 1-day sampling to check if the effects are likely to be seen in real data. Results. We show that the effect of the BLR scattering on the predicted time delay is very similar to the effect of rising height of the X-ray source, without any BLR contribution. This brings additional degeneracy if we want in the future to recover the parameters of the system from the observed time delays in a specific object. Both effects, however, modify the slope of the delay versus wavelength curve when plotted in log space which opens a way to obtaining bare disk time delay needed for cosmology. In addition, when the disk irradiation is strong, the modification of the predicted delay by the BLR scattering and by X-ray source height become considerably different. However, the required data quality must be high since our artificial curves with 1-day sampling give rather noisy results., Comment: 14 pages, 16 figures, and 1 table. Published in A&A. Comments are welcome

This paper is the first in a series of preparing and analyzing spectral and other properties for a database of already discovered Changing-Look active galactic nuclei (CL AGNs). Here, we focus on the spectral fitting and analysis of broad emission lines in a sample of 93 CL AGNs collected from the literature with existing SDSS/BOSS/eBOSS spectroscopy where the H$\beta$ emission line profile does not completely disappear in any epochs. Additionally, we have gathered older/newer spectral epochs from all the available SDSS data releases to make the database more complete. We use PyQSOFit and perform a homogeneous spectral decomposition of all of our SDSS spectra and tabulate the AGN continuum and emission line properties per epoch per source, chronologically. This further allows us to categorize the sources in our sample as \textit{Turn-On} or \textit{Turn-Off} and subsequently check for repeated occurrences of such phases. We then estimate the black hole mass (M$_{\rm BH}$) and the Eddington ratio (L$_{\rm bol}$/L$_{\rm Edd}$) per epoch per source where the required parameters are available and well-estimated. We realize the movement of the source in the M$_{\rm BH}$ versus L$_{\rm bol}$/L$_{\rm Edd}$ plane allowing us to check for systematic changes in the source's fundamental properties. We then track their transition along the optical plane of the Eigenvector 1 (EV1) schema and categorize sources that either stay within the same Population (A or B) or make an inter-population movement as a function of spectral epoch. We also test the Balmer decrement (H$\alpha$/H$\beta$) of a subset of our sample of CL AGNs as a function of time and AGN luminosity., Comment: 35 pages, 21 figures, 4 tables, accepted for publication in The Astrophysical Journal Supplement Series

Narrow Line Seyfert 1 (NLSy1) galaxies have been shown to have high Eddington ratios and relatively small black hole masses. The measurement of the black hole masses is based on the virial relation which is dependent on the distribution of the line-emitting gas and the viewing angle to the source. Spectropolarimetry enables us to probe the geometry of this line-emitting gas and allows us to estimate independently the viewing angle of the source by comparing the spectrum viewed under natural light and in the polarized light. We performed spectropolarimetric observations of three NLSy1 - Mrk 1044, SDSS J080101.41+184840.7, and IRAS 04416+1215 using the European Southern Observatory's Very Large Telescope. We use ESO Reflex workflow to perform standard data reduction and extract the natural and polarized spectra. We estimate the Stokes parameters and the viewing angles of the three sources. We model the Stokes parameters and infer the properties of the scattering media - located in the equatorial and polar regions, and simulate the spectra observed both in natural light and in polarized light using the polarization radiative transfer code STOKES. We confirm that all three sources are high Eddington ratio objects. We are successful in recovering the observed H$\alpha$ line profile both in the natural and polarized light using the STOKES modelling. We recover the polarization fractions of the order of 0.2-0.5% for the three sources. Our principal component analysis shows that the sample of the 25 sources including our sources, Fairall 9 from Jiang et al. (2021), and sources from Capetti et al. (2021) are mainly driven by the black hole mass and Eddington ratio. We re-affirm the connection of the strength of the optical FeII emission with the Eddington ratio, but the dependence on the viewing angle is moderate, more like a secondary effect., Comment: 25 pages, 21 figures, 11 tables, to be published in Astronomy & Astrophysics, version after proofs

The study of the broad-line region (BLR) using reverberation mapping has allowed us to establish an empirical relation between the size of this line emitting region and the continuum luminosity that drives the line emission (i.e. the R$_{\rm BLR}$-L$_{\rm 5100}$ relation). To realize its full potential, the intrinsic scatter in the R$_{\rm BLR}$-L$_{\rm 5100}$ relation needs to be understood better. The Eddington ratio plays a key role in addressing this problem. On the other hand, the Eigenvector 1 schema has helped to reveal an almost clear connection between the Eddington ratio and the strength of the optical FeII emission which has its origin from the BLR. This paper aims to reveal the connection between theoretical entities, like, the ionization parameter (U) and cloud mean density (n$_{\rm H}$) of the BLR, with physical observables obtained directly from the spectra, such as optical FeII strength (R$_{\rm FeII}$) that has immense potential to trace the accretion rate. We utilize the photoionization code CLOUDY and perform a suite of models to reveal the BLR in the U-n parameter space and estimate RFeII. We compare the SEDs for a prototypical Population A and Population B source, I Zw 1 and NGC 5548, respectively, in this study. The results from the photoionization modelling are combined with existing reverberation mapped sources with observed R$_{\rm FeII}$ estimates, allowing us to provide an analytical formulation to tie together the aforementioned quantities. We utilize the comparison of the modelled equivalent widths for the low-ionization emission lines to their observed values to identify the optimal (U,n$_{\rm H}$). The recovery of the correct physical conditions in the BLR suggests that the BLR `sees' a different, filtered ionizing continuum with only a very small fraction (~1-10%) that leads to the line emission in the dustless, low-ionization BLR., Comment: 33 pages, 12 figures, 4 tables, accepted for publication in `Frontiers in Astronomy and Space Sciences: Rising Stars' Issue

The issue of the difference between optical and UV properties of radio-quiet and radio-loud (relativistically "jetted") active galactic nuclei (AGN) is a long standing one, related to the fundamental question of why a minority of powerful AGN possess strong radio emission due to relativistic ejections. This paper examines a particular aspect: the singly-ionized iron emission in the spectral range 4400 -- 5600 A, where the prominent HI H$\beta$ and [OIII] 4959, 5007 lines are also observed. We present a detailed comparison of the relative intensity of Fe II multiplets in the spectral types of the quasar main sequence where most jetted sources are found, and afterwards discuss radio-loud narrow-line Seyfert 1 (NLSy1) nuclei with $\gamma$-ray detection and with prominent Fe II emission. An Fe II template based on I Zw 1 provides an accurate representation of the optical Fe II emission for RQ and, with some caveats, also for RL sources. CLOUDY photoionization simulations indicate that the observed spectral energy distribution can account for the modest Fe II emission observed in composite radio-loud spectra. However, spectral energy differences alone cannot account for the stronger Fe II emission observed in radio-quiet sources, for similar physical parameters. As for RL NLSy1s, they do not seem to behave like other RL sources, likely because of their different physical properties that could be ultimately associated with a higher Eddington ratio., Comment: 27 pages, to appear in Universe as a contribution to the special issue "Panchromatic View of the Life-Cycle of AGN" (M. Berton, Editor)

Broad-band spectra of active galaxies contain a wide range of information that help reveal the nature and activity of the central continuum source and their immediate surroundings. Understanding the evolution of metals in the spectra of Active Galactic Nuclei (AGN) and linking them with the various fundamental black hole (BH) parameters, for example, BH mass, the bolometric luminosity of the source, its accreting power, can help address the connection between the growth of the BH across cosmic time. We investigate the role of selected metallicity indicators utilizing the rich spectroscopic database of emission lines covering a wide range in redshift in the recent spectroscopic data release of Sloan Digital Sky Survey (SDSS). We make careful filtering of the parent sample to prepare a pair of high-quality, redshift-dependent sub-samples and present the first results of the analysis here. To validate our findings from the simple correlations, we execute and evaluate the performance of a linear dimensionality reduction technique - principal component analysis (PCA), over our sub-samples and present the projection maps highlighting the primary drivers of the observed correlations. The projection maps also allow us to isolate peculiar sources of potential interest., Comment: 3 pages, 2 figures, submitted to the Proceedings of the 40th Polish Astronomical Society Meeting, 13-17 September 2021, Szczecin, Poland

Using bonafide black hole (BH) mass estimates from reverberation mapping and the line ratio [Si VI]1.963$\mu$m/Br$\gamma_{\rm broad}$ as tracer of the Active Galactic Nuclei (AGN) ionizing continuum, we find a novel BH-mass scaling relation of the form log($M_{\rm BH}) = (6.40\pm 0.17) - (1.99\pm 0.37) \times$ log ([Si VI]/Br$\gamma_{\rm broad})$, over the BH mass interval, $10^6 - 10^8$ M$_{\odot}$. The current sample consists of 21 Type-1 AGNs with the overall dispersion in our scaling relation at 0.47 dex, one that emulates the well-established M-$\sigma$ relation which shows a dispersion $\sim$0.44 dex. The new scaling offers an economic, and physically motivated alternative for BH estimate using single epoch spectra, avoiding large telescope time (reverberation mapping) or absolute flux calibration (the continuum luminosity method). With the advent of big surveys in the infrared in the near future, we aim to reduce the scatter in the relation by supplementing with more sources., Comment: 3 pages, 1 figure, submitted to the Proceedings of the 40th Polish Astronomical Society Meeting, 13-17 September 2021, Szczecin, Poland

The CaFe Project involves the study of the properties of the low ionization emission lines (LILs) pertaining to the broad-line region (BLR) in active galaxies. These emission lines, especially the singly-ionized iron (Fe II) in the optical and the corresponding singly-ionized calcium (Ca II) in the near-infrared (NIR) are found to show a strong correlation in their emission strengths, i.e., with respect to the broad H$\beta$ emission line, the latter also belonging to the same category of LILs. We outline the progress made in the past years that has developed our understanding of the location and the efficient production of these emission lines. We have yet to realize the full potential of Ca II emission and its connection to the black hole and the BLR parameters which can be useful in - (1) the classification of Type-1 active galactic nuclei (AGNs) in the context of the main sequence of quasars, (2) to realize an updated radius-luminosity relation wherein the inclusion of the strength of this emission line with respect to H$\beta$ can be an effective tracer of the accretion rate of the AGN, and, (3) the close connection of Ca II to Fe II can allow us to use the ratio of the two species to quantify the chemical evolution in these active galaxies across cosmic time. In this paper, we use our current sample and utilize a non-linear dimensionality reduction technique - t-distributed Stochastic Neighbour Embedding (tSNE), to understand the clustering in our dataset based on direct observables., Comment: 6 pages, 3 figures, submitted to the Proceedings of the 40th Polish Astronomical Society Meeting, 13-17 September 2021, Szczecin, Poland

We report preliminary results on the analysis of the continuum and H$\beta$ light curves of the type-1 active galactic nucleus (AGN) NGC 5548. We notice a clear signature of shallowing in the trend between the H$\beta$ and the continuum luminosities. We attempt the recovery of this observed H$\beta$ emission trend as a response to large continuum flux increase using CLOUDY photoionization simulations. We explore a wide range in the physical parameters space for modelling the H$\beta$ emission from the broad-line region (BLR) appropriate for this source. We employ a constant density, single cloud model approach in this study and successfully recover the observed shallowing of the H$\beta$ emission with respect to the rising AGN continuum. With our modelling, we are able to provide constraints on the local BLR cloud density and recover the BLR distances (from the continuum source) consistent with the H$\beta$ reverberation mapping estimates. We further discuss the implications of the BLR covering factor and their sizes on recovering the observed trend., Comment: 10 pages, 7 figures, accepted version after proofs in Astronomische Nachrichten (Special Issue, 13th Serbian Conference on Spectral Line Shapes in Astrophysics)

The CaFe Project involves the study of the properties of the low ionization emission lines (LILs) pertaining to the broad-line region (BLR) in active galaxies. These emission lines, especially the singly-ionized iron (Fe II) in the optical and the corresponding singly-ionized calcium (Ca II) in the near-infrared (NIR) are found to show a strong correlation in their emission strengths, i.e. with respect to the broad H$\beta$ emission line, the latter also belonging to the same category of LILs. The origin of this correlation is attributed to the similarity in the physical conditions necessary to emit these lines - especially in terms of the strength of the ionization from the central continuum source and the local number density of available matter in these regions. In this paper, we focus on the issue of the spectral energy distribution (SED) characteristic to a prototypical Type-1 Narrow-line Seyfert galaxy (NLS1) - I Zw 1. We extract the continuum from quasi-simultaneous spectroscopic measurements ranging from the near-UV ($\sim$1200A) to the near-infrared ($\sim$24000A) to construct the SED and supplement it with archival X-ray measurements available for this source. Using the photoionization code CLOUDY, we assess and compare the contribution of the prominent "Big Blue Bump" seen in our SED versus the SED used in our previous work, wherein the latter was constructed from archival, multi-epoch photometric measurements. Following the prescription from our previous work, we constrain the physical parameter space to optimize the emission from these LILs and discuss the implication of the use of a "better" SED., Comment: 11 pages, 5 figures, accepted version after proofs, Acta Astrophysica Taurica proceedings for Crimean-2021 AGN Conference "Galaxies with Active Nuclei on Scales from Black Hole to Host Galaxy" dedicated to the E.A. Dibai's 90th anniversary (13-17 September 2021)

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 the monitoring of the AGN continuum and MgII broad line emission for the quasar HE 0413-4031 ($z=1.38$) based on the six-year monitoring by the South African Large Telescope (SALT). We managed to estimate a time-delay of $302.6^{+28.7}_{-33.1}$ days in the rest frame of the source using seven different methods: interpolated cross-correlation function (ICCF), discrete correlation function (DCF), $z$-transformed DCF, JAVELIN, two estimators of data regularity (Von Neumann, Bartels), and $\chi^2$ method. This time-delay is below the value expected from the standard radius-luminosity relation. However, based on the monochromatic luminosity of the source and the SED modelling, we interpret this departure as the shortening of the time-delay due to the higher accretion rate of the source, with the inferred Eddington ratio of $\sim 0.4$. The MgII line luminosity of HE 0413-4031 responds to the continuum variability as $L_{\rm line}\propto L_{\rm cont}^{0.43\pm 0.10}$, which is consistent with the light-travel distance of the location of MgII emission at $R_{\rm out} \sim 10^{18}\,{\rm cm}$. Using the data of 10 other quasars, we confirm the radius-luminosity relation for broad MgII line, which was previously determined for broad H$\beta$ line for lower-redshift sources. In addition, we detect a general departure of higher-accreting quasars from this relation in analogy to H$\beta$ sample. After the accretion-rate correction of the light-travel distance, the MgII-based radius-luminosity relation has a small scatter of only $0.10$ dex., Comment: 39 pages (23 pages - Main text, 16 pages - Appendix), 21 figures, 14 tables; accepted for publication in the Astrophysical Journal

Modelling the low ionization lines (LIL) in active galactic nuclei still faces problems in explaining the observed equivalent widths (EWs). We examine the optical Fe II and near-infrared Ca II triplet (CaT) emission strengths using the photoionization code CLOUDY. Using an incident continuum for I Zw 1 - a prototypical Type-1 narrow-line Seyfert galaxy, we can recover the line ratios for the optical Fe II (i.e. R$_{\rm{Fe II}}$) and the NIR CaT (i.e. R$_{\rm{CaT}}$) in agreement to the observed estimates. Although, the pairs of (U,$\rm{n_{H}}$) that reproduce the conforming line ratios, unfortunately, do not relate to agreeable line EWs. We thus propose that the LIL region of the BLR cloud doesn't see the same continuum seen by a distant observer that is emanated from the accretion disk, rather it sees a filtered version of the original continuum. The assumption of the filtered continuum as the source of BLR irradiation recovers realistic EWs for LIL species. However, our study finds that to account for the adequate R$_{\rm{Fe II}}$ (Fe II/H$\beta$ flux ratio) emission, the BLR needs to be selectively overabundant in iron. On the other hand, the R$_{\rm{CaT}}$ (CaT/H$\beta$ flux ratio) emission spans a broader range from solar to super-solar metallicities. In all these models the BLR cloud density is found to be consistent with our conclusions from prior works, i.e. $\rm{n_{H}} \sim 10^{12}$ cm$^{-3}$. An interesting result obtained here is the reduction in the value of the metallicity by up to a factor 10 for the R$_{\rm{Fe II}}$ cases when the microturbulence is invoked, suggesting that microturbulence can act as an apparent metallicity controller for the Fe II. On the contrary, the R$_{\rm{CaT}}$ cases are rather unaffected by the effect of microturbulence., Comment: 15 pages, 8 figures, accepted in Astronomy & Astrophysics

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

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)

In the near future, the arrival of new surveys promises a significant increment in the number of quasars (QSO) at large redshift ranges. 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 analysed using the photometric reverberation mapping technique. Oncoming LSST data will decrease the uncertainties in the dark energy determination using reverberating--mapped sources, particularly at high redshifts. 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 LSST. With the onset of the LSST era, we expect a huge rise in the overall quasar counts and redshift range covered ($z\leq7.0$), which will provide a better constraint of AGN properties with cosmological purposes., Comment: 3 pages, 1 figure, contribution to the Proceedings of the XXXIX Polish Astronomical Society Meeting (09-12 September 2019, Olsztyn, Poland)

Broad-line regions (BLR) are one of the main components that constitute the phenomenological picture of active galaxies near the vicinity of the accreting supermassive black holes. Both theoretical and observational studies have shown that the BLR is made of dense, ionized gas clumps that have a strong virialized distribution at parsec-scale distances from the nuclei. Using a theoretically motivated photoionized gas model, I constrain the ionisation parameter (U) and cloud density (n$_{\rm{H}}$) as a function of the strength of the FeII emission. Recent observations in the reverberation mapping studies have contested the standard radius-luminosity relation showing increased dispersion in the relation, in particular, after the inclusion of highly accreting quasars. I incorporate the departure coefficient that accounts for this dispersion. This departure term in terms of the dimensionless accretion rate ($\dot{\mathcal{M}}$) and Eddington ratio ($L\mathrm{_{bol}}/L\mathrm{_{Edd}}$), also includes the virial factor that accounts for the BLR geometry. I then combine the fundamental plane relation for the BLR to connect the FeII strength (R$_{\rm{FeII}}$) in terms of U and n$_{\rm{H}}$ using selected values for the shape of the broad H$\beta$ profile., Comment: 6 pages, 3 figures, contribution to the Proceedings of the XXXIX Polish Astronomical Society Meeting (09-12 September 2019, Olsztyn, Poland)

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

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"

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

Bright active galaxies show a range of properties but many of these properties are correlated which has led to the concept of the Quasar Main Sequence. We test whether our current understanding of the quasar structure allows to reproduce the pattern observed in the optical plane formed by the kinematic line width of H$\beta$ and the relative importance of the Fe II optical emission. We performed simulations of the H$\beta$ and Fe II production using the code CLOUDY and well justified assumptions about the broad band spectra, distance of the emission line region, and the cloud properties. We show that the presence of the warm corona is an important element of the broad band spectrum which decreases the dependence of the relative Fe II emissivity on the Eddington ratio, and allows to reproduce the rare cases of the particularly strong Fe II emitters. Results are sensitive to the adopted cloud distance, and strong Fe II emission can be obtain either by adopting strongly super-solar metallicity, or much shorter distance than traditionally obtained from reverberation mapping. We modeled in a similar way the UV plane defined by the Mg II line and Fe II UV pseudo-continuum, but here our approach is less successful, in general overproducing the Fe II strength. We found that the Fe II optical and UV emissivity depend in a different way on the turbulent velocity and metallicity, and the best extension of the model in order to cover both planes is to allow very large turbulent velocities in the Broad Line Region clouds., Comment: 15 pages, 12 figures, accepted for publication in ApJ

The work initially started as a test to retrace the Shen & Ho (2014) Quasar Main Sequence diagram where they claimed that the parameter RFeII, which defines the Eigenvector 1 (EV1) is driven by the Eddington ratio alone. We subsequently construct a refined (error and redshift limited) sample from the original Shen et al. 2011 QSO catalogue. Based on our hypothesis - the main driver of the Quasar Main Sequence is the maximum of the accretion disk temperature defined by the Big Blue Bump on the Spectral Energy Distribution. We select the four extreme sources that have RFeII $\geq$ 4.0 and use the SED modelling code CIGALE to fit the multi-band photometric data for these sources. We also perform detailed spectral fitting including the FeII pseudo-continuum to estimate and compare the value of RFeII for them. We show the metallicity dependent FeII strength in the context of this study., Comment: International Astronomical Union Proceedings: IAU Symposium 341 - PanModel2018: Challenges in Panchromatic Galaxy Modelling with next generation facilities

Most results of the reverberation monitoring of active galaxies showed a universal scaling of the time delay of the Hbeta emission region with the monochromatic flux at 5100 A, with very small dipersion. Such a scaling favored the dust-based formation mechanism of the Broad Line Region (BLR). Recent reverberation measurements showed that actually a significant fraction of objects exhibits horter lags than the previously found scaling. Here we demonstrate that these shorter lags can be explained by the old concept of scaling of the BLR size with the ionization parameter. Assuming a universal value of this parameter and universal value of the cloud density reproduces the distribution of observational points in the time delay vs. monochromatic flux plane, provided that a range of black hole spins is allowed. However, a confirmation of the new measurements for low/moderate Eddington ratio sources is strongly needed before the dust-based origin of the BLR can be excluded., Comment: Accepted for publication in ApJ

Active galaxies form a clear pattern in the optical plane showing the correlation between the Full Width Half Maximum (FWHM) of the H$\beta$ line and the ratio of the Equivalent Width (EW) of the optical FeII emission and the EW(H$\beta$). This pattern is frequently refereed to as the quasar main sequence. In this paper we study the UV plane showing FWHM of the MgII line against the ratio of EW of UV FeII emission to EW(MgII). We show that the UV plane trends are different, with the underlying strong correlation between FWHM(MgII) and EW(MgII). This correlation is entirely driven by the choice of the continuum used to measure EW(MgII). If we correct for the strong curvature in the UV continuum related to the character of emission expected from accretion disk, the behavior of the FWHM vs. EW for MgII becomes similar to the behaviour of H$\beta$. Such a similarity is expected since both lines belong to the Low Ionization group of emission lines and come from a similar region. Such a similarity and weak coupling between FWHM and EW of the corresponding line is then also consistent with the Broad Line Region model based on the presence of dust in the accretion disk atmosphere., Comment: 15 pages, 11 figures

The concept of the quasar main sequence is very attractive since it stresses correlations between various parameters and implies the underlying simplicity. In the optical plane defined by the width of the H\{beta} line and the ratio of the equivalent width of the Fe II to H\{beta} observed objects form a characteristic pattern. In this paper, we use a physically motivated model to explain the distribution of quasars in the optical plane. Continuum is modelled as an accretion disk with a hard X-ray power law uniquely tight to the disk at the basis of observational scaling, and the Broad Line Region distance is determined also from observational scaling. We perform the computations of the FeII and H\{beta} line production with the code CLOUDY. We have only six free parameters for an individual source: maximum temperature of the accretion disk, Eddington ratio, cloud density, cloud column density, microturbulence, and iron abundance, and only the last four remain as global parameters in our modelling of the whole sequence. Our theoretically computed points cover well the optical plane part populated with the observed quasars, particularly if we allow for super-Solar abundance of heavy elements. Explanation of the exceptionally strong Fe II emitter requires a stronger contribution from the dark sides of the clouds. Analyzing the way how our model covers the optical plane we conclude that there is no single simple driver behind the sequence, as neither the Eddington ratio nor broad band spectrum shape plays the dominant role. Also, the role of the viewing angle in providing the dispersion of the quasar main sequence is apparently not as strong as expected., Comment: 16 pages, 10 figures, accepted for publication in ApJ

The work initially started as a test to retrace the Shen & Ho (2014) Quasar Main Sequence diagram where they (and references therein) claimed that the parameter RFeII, which defines the Eigenvector 1 (EV1) is driven by the Eddington ratio alone. We subsequently construct a refined (error and redshift limited) sample from the original Shen et. al (2011) QSO catalogue. Based on our hypothesis - the main driver of the Quasar Main Sequence is the maximum of the accretion disk temperature defined by the Big Blue Bump on the Spectral Energy Distribution. We select five extreme sources that have RFeII >= 4.0 and use the SED modelling code CIGALE to fit the multi-band photometric data for these sources. Incorporating the prescription for broad emission-line species dependent virial factors, we derive the black hole masses for the entire Shen+11 catalogue and compare their distribution with the black hole masses in the original catalogue. Here, we show the results of our analysis for one of the 5 extreme sources, SDSSJ082358.30+213545.2. We also show the detailed modelling, including the Fe II pseudo-continuum to estimate and compare the value of RFeII for this object., Comment: Proceedings of Science - Revisiting narrow-line Seyfert 1 galaxies and their place in the Universe, 9-13 April 2018, Padova, Italy

Eigenvector 1 (EV1) was found to be the dominant component behind the significant correlations for the measured parameters in quasar spectra (Boroson & Green, 1992). The parameter R$_{\mathrm{FeII}}$, which strongly correlates to the EV1, is the ${\mathrm{FeII}}$ strength, defined to be the ratio of the equivalent width of ${\mathrm{FeII}}$ to the equivalent width of ${\mathrm{H\beta}}$. This allows to construct a quasar main sequence analogous to the stellar properties driven HR diagram (Sulentic et al. 2001). We try to find the main driver behind the EV1 among the basic (theoretically motivated) parameters of an active nucleus (Eddington ratio, black hole mass, accretion rate, spin, and viewing angle). Based on theoretical modeling using the photoionization code CLOUDY (Ferland et al. 2013), we test the hypothesis that the physical driver of EV1 is the maximum of the accretion disk temperature ($\mathrm{T_{BBB}}$), reflected in the shape of the spectral energy distribution (SED). We have assumed that both H$\mathrm{\beta}$ and Fe${\mathrm{II}}$ emission come from the Broad Line Region represented as a constant density cloud in a plane-parallel geometry. We test the effect of changing Eddington ratio on the $\mathrm{R_{FeII} - T_{BBB}}$ trends with varying mean hydrogen densities. We also test the effect of adding microturbulence that affect the line intensities on the overall $\mathrm{R_{FeII} - T_{BBB}}$ picture., Comment: 4 pages, 4 figures, to be published in the proceedings of the 3rd Cosmology School (Polish Astronomical Society), 10-23 July 2017, Krak\'ow, Poland

Quasars are among the most luminous sources characterized by their broad band spectra ranging from radio through optical to X-ray band, with numerous emission and absorption features. Using the Principal Component Analysis (PCA), Boroson & Green (1992) were able to show significant correlations between the measured parameters. Among the significant correlations projected, the leading component, related to Eigenvector 1 (EV1) was dominated by the anti-correlation between the Fe${\mathrm{II}}$ optical emission and [OIII] line where the EV1 alone contained 30% of the total variance. This introduced a way to define a quasar main sequence, in close analogy to the stellar main sequence in the Hertzsprung-Russel (HR) diagram (Sulentic et. al 2001). Which of the basic theoretically motivated parameters of an active nucleus (Eddington ratio, black hole mass, accretion rate, spin, and viewing angle) is the main driver behind the EV1 yet remains to be answered. We currently limit ourselves to the optical waveband, and concentrate on theoretical modelling the Fe${\mathrm{II}}$ to H$\mathrm{\beta}$ ratio, and test the hypothesis that the physical driver of EV1 is the maximum of the accretion disk temperature, reflected in the shape of the spectral energy distribution (SED). We performed computations of the H$\mathrm{\beta}$ and optical Fe${\mathrm{II}}$ for a broad range of SED peak position using CLOUDY photoionisation code. We assumed that both H$\mathrm{\beta}$ and Fe${\mathrm{II}}$ emission come from the Broad Line Region represented as a constant density cloud in a plane-parallel geometry. We compare the results for two different approaches: (1) considering a fixed bolometric luminosity for the SED; (2) considering $\mathrm{L_{bol}/L_{Edd}}$ = 1., Comment: 38th Polish Astronomical Society Assembly Conference Proceedings, Zielona G\'ora, Poland

Quasars are complex sources, characterized by broad band spectra from radio through optical to X-ray band, with numerous emission and absorption features. However, Boroson & Green (1992) used Principal Component Analysis (PCA), and with this analysis they were able to show significant correlations between the measured parameters. The leading component, related to Eigenvector 1 (EV1) was dominated by the anticorrelation between the Fe${\mathrm{II}}$ optical emission and [OIII] line and EV1 alone contained 30% of the total variance. It opened a way in defining a quasar main sequence, in close analogy to the stellar main sequence on the Hertzsprung-Russel (HR) diagram (Sulentic et al. 2001). The question still remains which of the basic theoretically motivated parameters of an active nucleus (Eddington ratio, black hole mass, accretion rate, spin, and viewing angle) is the main driver behind the EV1. Here we limit ourselves to the optical waveband, and concentrate on theoretical modelling the Fe${\mathrm{II}}$ to H$\mathrm{\beta}$ ratio, and we test the hypothesis that the physical driver of EV1 is the maximum of the accretion disk temperature, reflected in the shape of the spectral energy distribution (SED). We performed computations of the H$\mathrm{\beta}$ and optical Fe${\mathrm{II}}$ for a broad range of SED peak position using CLOUDY photoionisation code. We assumed that both H$\mathrm{\beta}$ and Fe${\mathrm{II}}$ emission come from the Broad Line Region represented as a constant density cloud in a plane-parallel geometry. We expected that a hotter disk continuum will lead to more efficient production of Fe${\mathrm{II}}$ but our computations show that the Fe${\mathrm{II}}$ to H$\mathrm{\beta}$ ratio actually drops with the rise of the disk temperature. Thus either hypothesis is incorrect, or approximations used in our paper for the description of the line emissivity is inadequate., Comment: 12 pages, 4 figures, Accepted for publication in the Journal Frontiers in Astronomy and Space Sciences