Spectra Of Active Galactic Nuclei Research Articles

The VMC Survey – LI. Classifying extragalactic sources using a probabilistic random forest supervised machine learning algorithm

Abstract We used a supervised machine learning algorithm (probabilistic random forest) to classify ∼130 million sources in the VISTA Survey of the Magellanic Clouds (VMC). We used multi-wavelength photometry from optical to far-infrared as features to be trained on, and spectra of Active Galactic Nuclei (AGN), galaxies and a range of stellar classes including from new observations with the Southern African Large Telescope (SALT) and SAAO 1.9m telescope. We also retain a label for sources that remain unknown. This yielded average classifier accuracies of ∼79% (SMC) and ∼87% (LMC). Restricting to the 56,696,719 sources with class probabilities (Pclass) > 80% yields accuracies of ∼90% (SMC) and ∼98% (LMC). After removing sources classed as ‘Unknown’, we classify a total of 707,939 (SMC) and 397,899 (LMC) sources, including >77,600 extragalactic sources behind the Magellanic Clouds. The extragalactic sources are distributed evenly across the field, whereas the Magellanic sources concentrate at the centres of the Clouds, and both concentrate in optical/IR colour–colour/magnitude diagrams as expected. We also test these classifications using independent datasets, finding that, as expected, the majority of X-ray sources are classified as AGN (554/883) and the majority of radio sources are classed as AGN (1756/2694) or galaxies (659/2694), where the relative AGN–galaxy proportions vary substantially with radio flux density. We have found: >49,500 hitherto unknown AGN candidates, likely including more AGN dust dominated sources which are in a critical phase of their evolution; >26,500 new galaxy candidates and >2800 new Young Stellar Object (YSO) candidates.

Modelling absorption and emission profiles from accretion disc winds with WINE

Context. Fast and massive winds are ubiquitously observed in the UV and X-ray spectra of active galactic nuclei (AGNs) and other accretion-powered sources. Several theoretical and observational pieces of evidence suggest they are launched at accretion disc scales, carrying significant mass and angular momentum. Thanks to such high-energy output, they may play an important role in transferring the energy released by accretion to the surrounding environment. In the case of AGNs, this process can help to set the so-called co-evolution between an AGN and its host galaxy, which mutually regulates their growth across cosmic time. To precisely assess the effective role of UV and X-ray winds at accretion disc scales, it is necessary to accurately measure their properties, including mass and energy rates. However, this is a challenging task, due to both the limited signal-to-noise ratio of available observations and the limitations of the models currently used in the spectral analysis. Aims. We aim to maximise the scientific return of current and future observations by improving the theoretical modelling of these winds through our Winds in the Ionised Nuclear Environment (WINE) model. WINE is a spectroscopic model specifically designed for disc winds in AGNs and compact accreting sources, which couples photoionisation and radiative transfer with special relativistic effects and a three-dimensional model of the emission profiles. Methods. We explore with WINE the main spectral features associated with the disc winds in AGNs, with a particular emphasis on the detectability of the wind emission in the total transmitted spectrum. We explore the impact of the wind ionisation, column density, velocity field, and geometry in shaping the emission profiles. We simulated observations with the X-ray microcalorimeter Resolve on board the recently launched XRISM satellite and the X-IFU on board the future Athena mission. This allows us to assess the capabilities of these telescopes in the study of disc winds in X-ray spectra of AGNs for the typical physical properties and exposure times of the sources included in the XRISM performance verification phase. Results. The wind kinematic and geometry (together with the ionisation and column density) deeply affect both shape and strength of the wind spectral features. Thanks to this, both Resolve and, on a longer timescale, X-IFU will be able to accurately constrain the main properties of disc winds over a broad range of ionisation, column densities, and covering factors. We also investigate the impact of the spectral energy distribution (SED) on the resulting appearance of the wind. Our findings reveal a dramatic difference in the gas opacity when using a soft, Narrow Line Seyfert 1-like SED compared to a canonical powerlaw SED with a spectral index of Γ ≈ 2.

Separating the Spectral Counterparts in NGC 1275/Perseus Cluster in X-Rays

We present a model-independent method for separating the spectral counterparts of the active galactic nucleus (AGN) NGC 1275 from the surrounding emission of the Perseus cluster, as observed by Suzaku/XIS cameras. The Perseus cluster emission extends to higher energies than typically observed in AGN environments, reaching up to 9–10 keV. This necessitates precise separation of AGN and cluster spectra. To circumvent the degeneracy arising from numerous spectral fitting parameters, including elemental abundances, thermal and Compton emissions from the nucleus, and spectral parameters of the jet synchrotron self-Compton/inverse Compton emissions, we avoid traditional spectral fitting methods. Instead, we leverage spatial resolution and employ a double background subtraction approach. We apply this procedure to the complete set of Suzaku/XIS observational data for NGC 1275, resulting in cleaned spectra and a light curve of the AGN emission in this system. To demonstrate the applicability of our method, we also utilize the available XMM-Newton/EPIC data.

Systematic collapse of the accretion disc across the supermassive black hole population

ABSTRACT The structure of the accretion flow on to supermassive black holes is not well understood. Standard disc models match to zeroth-order in predicting substantial energy dissipation within optically thick material producing a characteristic strong blue/ultraviolet continuum. However, they fail at reproducing more detailed comparisons to the observed spectral shapes along with their observed variability. Based on stellar mass black holes within our Galaxy, accretion discs should undergo a transition into an X-ray hot, radiatively inefficient flow, below a (mass scaled) luminosity of $\sim 0.02\, L_{\rm {Edd}}$. While this has been seen in limited samples of nearby low-luminosity active galactic nuclei (AGN) and few rare changing-look AGN, it is not at all clear whether this transition is present in the wider AGN population across cosmic time. A key issue is the difficulty in disentangling a change in spectral state from increased dust obscuration and/or host galaxy contamination, effectively drowning out the AGN emission. Here, we use the new eROSITA eFEDS Survey to identify unobscured AGN from their X-ray emission, matched to excellent optical imaging from Subaru’s Hyper Suprime-Cam; allowing the subtraction of the host galaxy contamination. The resulting, uncontaminated, AGN spectra reveal a smooth transition from a strongly disc-dominated state in bright AGN, to the collapse of the disc into an inefficient X-ray plasma in the low-luminosity AGN, with the transition occurring at $\sim 0.02\, L_{\rm {Edd}}$; revealing fundamental aspects of accretion physics in AGN.

X-ray view of dissipative warm corona in active galactic nuclei

Context. Detections of the X-ray spectra of active galactic nuclei (AGNs) typically reveal a noticeable excess of soft X-rays beyond the extrapolation of the power-law trend observed between 2 and 10 keV. However, the cause of this surplus remains unclear. In the scenario of soft Comptonization, observations suggest a warm corona temperature, ranging from 0.1 to 1 keV, and an optical depth of approximately 10 to 30. Furthermore, according to radiative constraints derived from spectral analyses employing Comptonization models, it has been suggested that most of the accretion power is released within the warm corona. At the same time, the disk beneath it is largely non-dissipative, mainly emitting the reprocessed radiation from the corona. Aims. We tested the dissipative warm corona model using the radiative transfer code TITAN-NOAR on a sample of 82 XMM-Newton EPIC-pn observations of 21 AGNs. Based on a spectral modeling of the X-ray data, we aim to estimate the total amount of internal heating inside the warm corona on top of the accretion disk. Methods. By modeling the 0.3–10 keV EPIC-pn spectra with the TITAN-NOAR model component, we estimated the internal heating and optical depth of the warm corona and checked their correlations with such global parameters as the hot corona spectral index, black hole mass, and accretion rate. From the model normalization, we computed the radial extent of the warm corona on top of the cold accretion disk. Results. Our model infers the presence of dissipative warm corona, with optical depths distributed across the range of ∼6–30 and a total internal heating in the range of ∼1–29 × 10−23 erg s−1 cm3. We do not detect any variation between these properties and global properties, such as the black hole mass and accretion rate. The extent of the warm corona is spread across a broad range, from 7–408 gravitational radii, and we find that the warm corona is more extended for higher accretion rates. Conclusions. Soft excess emission is ubiquitous across a wide mass range and accretion rate in AGNs. We confirm that the warm corona responsible for producing the soft X-ray excess is highly dissipative with greater optical depths being associated with lower internal heating and vice versa. The cold standard accretion disk regulates the extent of the warm corona.

The X-ray variability of active galactic nuclei: Power spectrum and variance analysis of the Swift/BAT light curves

Aims. We study the X-ray power spectrum of active galactic nuclei (AGN) in order to investigate whether Seyfert I and II power spectra are similar or not and whether AGN variability depends on the mass and accretion rate of black holes as well as to compare the power spectra of AGN with the power spectra of Galactic X-ray black hole binaries. Method. We used 14–195 keV band light curves from the 157-month Swift/BAT hard X-ray survey, and we computed the mean power spectrum and excess variance of AGN in narrow black hole mass and AGN luminosity bins. We fitted a power-law model to the AGN power spectra, and we investigated whether the power spectrum parameters and the excess variance depend on the black hole mass, luminosity, and accretion rate of AGN. Results. We found the Seyfert I and Seyfert II power spectra to be identical, in agreement with AGN unification models. The mean AGN X-ray power spectrum has the same power-law like shape, with a slope of −1 in all AGN irrespective of their luminosity and black hole mass. We did not detect any flattening to a slope of zero at frequencies as low as 10−9 Hz. We detected an anti-correlation between the power spectral density function (PSD) amplitude and the accretion rate, similar to what has been seen in the past in the 2–10 keV band. This implies that the variability amplitude in AGN decreases with an increasing accretion rate. The universal AGN power spectrum is consistent with the mean 2–9 keV band Cyg X-1 power spectrum in its soft state. We detected a small difference in amplitude, but this is probably due to the difference in energy. Conclusions. The mean low-frequency AGN X-ray power spectrum is consistent with the extension of the mean 0.01–25 Hz Cyg X-1 power spectrum in its soft state to lower frequencies. We cannot prove that the mean AGN PSD is analogous to the mean Cyg X-1 PSD in its soft state, as we do not know the location of the high-frequency break in the hard X-ray AGN PSDs. However, if this is the case, then the accretion disc in AGN probably extends to the radius of the innermost circular stable orbit (as is probably the case with the black hole binaries in their soft state). The X-ray corona will then be located on top, illuminating the disc and producing the X-ray reflection and disc reverberation phenomena commonly observed in these objects. Furthermore, the agreement between the PSD amplitude in AGN and the Cyg X-1 (either in the soft or the hard state) over many decades in frequency indicates that the X-ray variability process is probably the same in all accreting objects, irrespective of the mass of the compact object. We plan to investigate this issue further in the near future.

Euclidpreparation

The statistical distribution and evolution of key properties of active galactic nuclei (AGN), such as their accretion rate, mass, and spin, remains a subject of open debate in astrophysics. The ESAEuclidspace mission, launched on July 1 2023, promises a breakthrough in this field. We create detailed mock catalogues of AGN spectra from the rest-frame near-infrared down to the ultraviolet – including emission lines – to simulate whatEuclidwill observe for both obscured (type 2) and unobscured (type 1) AGN. We concentrate on the red grisms of the NISP instrument, which will be used for the wide-field survey, opening a new window for spectroscopic AGN studies in the near-infrared. We quantify the efficiency in the redshift determination as well as in retrieving the emission line flux of the Hα+[N II] complex, asEuclidis mainly focused on this emission line, given that it is expected to be the brightest one in the probed redshift range. Spectroscopic redshifts are measured for 83% of the simulated AGN in the interval where the Hαis visible (i.e. 0.89 < z < 1.83 at a line flux of > 2 × 10−16erg s−1cm−2, encompassing the peak of AGN activity atz ≃ 1 − 1.5) within the spectral coverage of the red grism. Outside this redshift range, the measurement efficiency decreases significantly. Overall, a spectroscopic redshift iscorrectly determined for about 90% of type 2 AGN down to an emission line flux of roughly 3 × 10−16erg s−1cm−2, and for type 1 AGN down to 8.5 × 10−16erg s−1cm−2. Recovered values for black hole mass show a small offset with respect to the input values by about 10%, but the agreement is good overall. With such a high spectroscopic coverage atz < 2, we will be able to measure AGN demography, scaling relations, and clustering from the epoch of the peak of AGN activity down to the present-day Universe for hundreds of thousands of AGN with homogeneous spectroscopic information.

TELAMON: Effelsberg monitoring of AGN jets with very-high-energy astroparticle emission

Aims.We introduce the TELAMON program which is using the Effelsberg 100-m telescope to monitor the radio spectra of active galactic nuclei (AGN) under scrutiny in astroparticle physics, specifically TeV blazars and candidate neutrino-associated AGN. Here, we present and characterize our main sample of TeV-detected blazars.Methods.We analyzed the data sample from the first ∼2.5 yr of observations between August 2020 and February 2023 in the range from 14 GHz to 45 GHz. During this pilot phase, we observed all 59 TeV-detected blazars in the Northern Hemisphere (i.e., Dec > 0°) known at the time of observation. We discuss the basic data reduction and calibration procedures used for all TELAMON data and introduce a sub-band averaging method used to calculate average light curves for the sources in our sample.Results.The TeV-selected sources in our sample exhibit a median flux density of 0.12 Jy at 20 mm, 0.20 Jy at 14 mm, and 0.60 Jy at 7 mm. The spectrum for most of the sources is consistent with a flat radio spectrum and we found a median spectral index (S(ν)∝να) ofα = −0.11. Our results on flux density and spectral index are consistent with previous studies of TeV-selected blazars. Compared to the GeV-selected F-GAMMA sample, TELAMON sources are significantly fainter in the radio band. This is consistent with the double-humped spectrum of blazars being shifted towards higher frequencies for TeV-emitters (in particular for high-synchrotron peaked BL Lac type objects), which results in a lower radio flux density. The spectral index distribution of our TeV-selected blazar sample is not significantly different from the GeV-selected F-GAMMA sample. Moreover, we present a strategy to track the light curve evolution of sources in our sample for future variability and correlation analysis.

UV FeII emission model of HE 0413−4031 and its relation to broad-line time delays

Context. FeII emission is a well-known contributor to the UV spectra of active galactic nuclei and the modeling of this part may affect the results obtained for the MgII λ2800 emission, which is one of the lines used for black hole mass measurements and cosmological applications. Aims. We tested different FeII emission models when modeling the UV emission of the intermediate-redshift quasar HE 0413−4031 to see how the use of a specific template affects the MgII λ2800 line properties and the measurement of the MgII λ2800 and UV FeII time delays with respect to the continuum. Methods. We used the 11-year monitoring of the selected quasar HE 0413−4031 with the South African Large Telescope (SALT), and we supplemented this monitoring with the near-IR spectrum taken with the SOAR telescope, which gave access to the Hβλ4861 and [OIII] λλ4959, 5007 emission lines at the rest frame and allowed for a precise measurement of the redshift. Results. A new redshift determination (z = 1.39117 ± 0.00017) using [OIII] λλ4959, 5007 gave a very different value than the previous determination based only on the UV FeII pseudocontinuum (z = 1.3764). It favors a different decomposition of the spectrum into MgII and UV FeII emissions. The line characteristics and the time delay of the MgII emission (224−23+21 days) are not significantly affected. However, in comparison with the previous analysis, the rest-frame UV FeII time delay (251−7+9 days) is consistent with the inferred UV FeII line full width at half maximum (FWHM) of 4200 km s−1 that is only slightly smaller than the MgII line FWHM. Hence the FeII-emitting material is more distant than the MgII-emitting gas in HE 0413−4031 by ∼0.023 pc (4700 AU). The inferred velocity shift of both MgII and UV FeII lines with respect to the systemic redshift is now rather low, below 300 km s−1. In addition, we constructed an updated MgII radius-luminosity (R − L) relation from 194 sources, which is more than double the previous sample. The MgII R − L relation is flatter than the UV FeII, optical FeII, and HβR − L relations. While the new decomposition of the spectrum is satisfactory, we see a need to create better FeII templates using the newest version of the code CLOUDY.

Model of a “Warm Corona” as the origin of the soft X-ray excess of active galactic nuclei

Abstract The soft X-ray excess in the spectra of active galactic nuclei is characterized by similar electron temperatures of 0.1–0.3 keV and similar photon indices around 2.2–3, if fitted with inverse Comptonization. It remains a puzzle why both values are not sensitive to the black hole mass nor the accretion rate. Supposing that the scattering-dominated surface layer of an accretion disk can act as a warm corona, we construct a vertical one-zone model to understand what determines its temperature. By solving the equations of (1) the condition for the effective optical depth, (2) the energy balance, and (3) the dominance of the Compton cooling over the bound–free cooling, we could reproduce the basic observational features of the soft excess, provided that anomalous heating (excess heating other than what is expected by local energy dissipation) takes place in the warm corona in agreement with similar studies done so far. The similar temperatures can be understood, since both the anomalous heating and Compton cooling rates are proportional to the dissipation rate of the accretion energy, while similar photon indices are a natural consequence of the fact that observed photons are finally emitted from the layer of Compton y ∼ 1. The soft excess is not observed in black hole binaries, since disk temperatures are too high for the Compton scattering to work as cooling. The derived temperatures are somewhat of an underestimation, however. This may indicate a necessity of multi-zone corona structure. The stability of the warm corona and its consequences are briefly discussed.

Exploring the AGN Accretion Disks Using Continuum Reverberation Mapping

In the innermost regions of Active Galactic Nuclei (AGN), matter is understood to be flowing onto the Supermassive black hole (SMBH), which forms an accretion disk. This disk is responsible for the optical/UV continuum emission observed in the spectra of AGN. Reverberation Mapping of the accretion disk using multiple bands can yield the structure of the disk. The emission is expected to be of the black body type peaking at different wavelengths. Hence, depending on the temperature of the disk, continuous, simultaneous monitoring in multiple wavelength ranges to cover hotter inner regions and cooler outer regions can yield the structure and temperature profile of the accretion disk itself. In this study, we present initial results from our accretion disk reverberation mapping campaign targeting AGN with Super High Eddington Accreting Black Holes (SEAMBH). Our analysis on one of the sources, IRAS 04416+1215, based on the broadband observations using the Growth India telescope (GIT), reveals that the size of the accretion disk for this source, calculated by cross-correlating the continuum light curves is larger than expected from the theoretical model. We fit the light curves directly using the thin disk model available in JAVELIN and find that the disk sizes are approximately four times larger than expected from the Shakura Sunyaev (SS) disk model. Further studies are needed to understand better the structure and physics of AGN accretion disks and their role in the evolution of galaxies.

Time Evolving Photo Ionisation Device (TEPID): A novel code for out-of-equilibrium gas ionisation

Context. Photoionisation is one of the main mechanisms at work in the gaseous environment of bright astrophysical sources. A great deal of information on the gas physics, chemistry and kinematics, and on the ionising source itself, can be gathered through optical to X-ray spectroscopy. While several public time equilibrium photoionisation codes are readily available and can be used to infer average gas properties at equilibrium, time-evolving photoionisation models have only very recently started to become available. They are needed when the ionising source varies faster than the typical gas equilibration timescale. Using equilibrium models to analyse spectra of non-equilibrium photoionised gas may lead to inaccurate results, and prevents a solid assessment of gas density, physics, and geometry. Aims. Our main objective is to present and make available the Time-Evolving PhotoIonisation Device (TEPID), a new code that self-consistently solves time evolving photoionisation equations (both thermal and ionisation balance) and accurately follows the response of the gas to changes in the ionising source. Methods. TEPID self-consistently follows the gas temperature and ionisation in time by including all the main ionisation/recombination and heating/cooling mechanisms. The code takes in input the ionising light curve and spectral energy distribution and solves the time-evolving equations as a function of gas electron density and of time. The running time is intelligently optimised by an internal algorithm that initially scans the input light curve to set a time-dependent integration frequency. The code is built in a modular way, can be applied to a variety of astrophysical scenarios and produces time-resolved gas absorption spectra to fit the data. Results. To describe the structure and main features of the code, we present two applications of TEPID to two dramatically different astrophysical scenarios: the typical ionised absorbers observed in the X-ray spectra of active galactic nuclei (e.g. warm absorbers and ultra-fast outflows), and the circumburst environment of a gamma-ray burst. For both cases we show how the gas energy and ionisation balances vary as a function of time, gas density and distance from the ionising source. We show that time-evolving photoionisation leads to unique ionisation patterns that cannot be reproduced by stationary photoionisation codes when the gas is out of equilibrium. This demonstrates the need for codes such as TEPID in view of the unprecedented capabilities that will be offered by the upcoming high-resolution X-ray spectrometers on board missions like XRISM or Athena.

Influence of the optical Fe II quasi-continuum on measuring the spectral parameters of active galactic nuclei

Aims. We explore the influence of optical Fe II quasi-continuum on the measured spectral parameters in the λλ4150−5500 Å range for the spectra of Type 1 active galactic nuclei (AGNs). Methods. We assume that the broad line region is composed of two sub-regions: the very broad line region (VBLR) and the intermediate line region (ILR). We constructed a large set of synthetic AGN spectra by taking different portions of the VBLR and ILR contributions, where initially the VBLR and ILR model spectra were constructed on the basis of prototypes of two observed spectra with dominant VBLR (i.e. ILR) emission. To investigate the influence of the optical Fe II quasi-continuum on the AGN measured spectral parameters, we fit the power-law continuum and emission lines in a set of model spectra, as commonly done for observed AGN spectra. We then compared the spectral parameters obtained after the fitting procedure with those of the model. Results. We find that the optical Fe II quasi-continuum can be very strong in the case of spectra with strong and very broad Fe II lines and it is difficult to fully separate it from the power-law continuum. This gives the effect of a slightly underestimated Hβ width and underestimated fluxes of the Hβ and Fe II lines, while the continuum flux is then slightly overestimated. The most affected spectral parameters are the line equivalent widths (EWs), especially EW Fe II, which may be strongly underestimated. We discuss the possible underlying physics in the quasar main sequence, as implied by the results of our spectral modelling. We find that the set of AGN model spectra assuming different ILR and VBLR contributions can aptly reproduce the quasar main sequence, that is, the full width at half maximum Hβ versus Fe II/Hβ anti-correlation, where both parameters in this anti-correlation are strongly dependent on the ILR and VBLR contribution rate.

Testing Physical Scenarios for the Reflection Features of Type-1 AGNs using XMM-Newton and NuSTAR Simultaneous Observations

Above ∼3 keV, the X-ray spectrum of active galactic nuclei (AGNs) is characterized by the intrinsic continuum and Compton reflection features. For type-1 AGNs, several regions may contribute to the reflection. To investigate the nature of the reflecting medium, we perform a systematic analysis of the reflector using XMM-Newton and Nuclear Spectroscopic Telescope Array observations of a sample of 22 type-1 AGNs. We create a baseline model that includes Galactic absorption and an intrinsically absorbed power law, plus a reflection model. We test a set of nine reflection models in a subsample of five objects. Based on these results, we select three models to be tested on the entire sample, accounting for distinct physical scenarios: neutral/distant reflection, ionized/relativistic reflection, and neutral/distant+ionized/relativistic reflection, namely, a hybrid model. We find that 18 sources require the reflection component to fit their spectra. Among them, 67% prefer the hybrid model. Neutral and ionized models are equally preferred by three sources. We conclude that both the neutral/distant reflector most probably associated with the inner edges of the torus and the ionized/relativistic reflector associated with the accretion disk are required to describe the reflection in type-1 AGNs.

Star Formation in Self-gravitating Disks in Active Galactic Nuclei. III. Efficient Production of Iron and Infrared Spectral Energy Distributions

Strong iron lines are a common feature of the optical spectra of active galactic nuclei (AGNs) and quasars from z ∼ 6−7 to the local universe, and [Fe/Mg] ratios do not show cosmic evolution. During active episodes, accretion disks surrounding supermassive black holes (SMBHs) inevitably form stars in the self-gravitating part, and these stars accrete with high accretion rates. In this paper, we investigate the population evolution of accretion-modified stars (AMSs) to produce iron and magnesium in AGNs. The AMSs, as a new type of star, are allowed to have any metallicity but without significant loss from stellar winds, since the winds are choked by the dense medium of the disks and return to the core stars. Mass functions of the AMS population show a pile-up or cutoff pile-up shape in top-heavy or top-dominant forms if the stellar winds are strong, consistent with the narrow range of supernovae (SNe) explosions driven by the known pair-instability. This provides an efficient way to produce metals. Meanwhile, SN explosions support an inflated disk as a dusty torus. Furthermore, the evolving top-heavy initial mass functions lead to bright luminosity in infrared bands in dusty regions. This contributes a new component in infrared bands, which is independent of the emissions from the central part of accretion disks, appearing as a long-term trending of the NIR continuum compared to optical variations. Moreover, the model can be further tested through reverberation mapping of emission lines, including LIGO/LISA detections of gravitational waves and signatures from spatially resolved observations of GRAVITY+/VLTI.

Fantastic Fits with fantasy of Active Galactic Nuclei Spectra: Exploring the Fe ii Emission near the Hα Line

In this study, a refined approach for multicomponent fitting of active galactic nuclei (AGNs) spectra is presented utilizing the newly developed Python code Fully Automated pythoN Tool for AGN Spectra analYsis (fantasy). AGN spectra are modeled by simultaneously considering the underlying broken power-law continuum, predefined emission line lists, and an Fe ii model, which is here extended to cover the wavelength range 3700–11000 Å. The Fe ii model, founded solely on atomic data, effectively describes the extensive emission of the complex iron ion in the vicinity of the Hγ and Hβ lines, as well as near the Hα line, which was previously rarely studied. The proposed spectral fitting approach is tested on a sample of high-quality AGN spectra from the Sloan Digital Sky Survey Data Release 17. The results indicate that when Fe ii emission is present near Hβ, it is also detected redward from Hα, potentially contaminating the broad Hα line wings and thus affecting the measurements of its flux and width. The production of Fe ii emission is found to be strongly correlated with Eddington luminosity and appears to be controlled by a similar mechanism as the hydrogen Balmer lines. The study highlights the benefits of fitting type 1 AGN spectra with the fantasy code, pointing that it may be used as a robust tool for analyzing a large number of AGN spectra in the coming spectral surveys.

What Can We Learn about Compton-thin AGN Tori from Their X-Ray Spectra?<sup>1</sup>

We have developed a Monte Carlo code for simulation of X-ray spectra of active galactic nuclei (AGN) based on a model of a clumpy obscuring torus. Using this code, we investigate the diagnostic power of X-ray spectroscopy of obscured AGN with respect to the physical properties and orientation of the torus, namely: the average column density, \(\langle N_{\textrm{H}}\rangle\), the line-of-sight column density, \(N_{\textrm{H}}\), the abundance of iron, \(A_{\textrm{Fe}}\), the clumpiness (i.e. the average number of gas clouds along the line of sight), \(\langle N\rangle\), and the viewing angle, \(\alpha\). In this first paper of a series, we consider the Compton-thin case, where both \(\langle N_{\textrm{H}}\rangle\) and \(N_{\textrm{H}}\) do not exceed \(10^{24}\) cm\({}^{-2}\). To enable quantitative comparison of the simulated spectra, we introduce five measurable spectral characteristics: the low-energy hardness ratio (ratio of the continuum fluxes in the 7–11 keV and 2–7 keV energy bands), the high-energy hardness ratio (ratio of the continuum fluxes in the 10–100 keV and 2–10 keV energy bands), the depth of the iron K absorption edge, the equivalent width of the Fe K\(\alpha\) line, and the fraction of the Fe K\(\alpha\) flux contained in the Compton shoulder. We demonstrate that by means of X-ray spectroscopy it is possible to tightly constrain \(\langle N_{\textrm{H}}\rangle\), \(N_{\textrm{H}}\) and \(A_{\textrm{Fe}}\) in the Compton-thin regime, while there is degeneracy between clumpiness and viewing direction.

The universal shape of the X-ray variability power spectrum of AGN up to z ∼ 3

Aims. We study the ensemble X-ray variability properties of active galactic nuclei (AGN) over large ranges of timescale (20 ks ≤ T ≤ 14 yr), redshift (0 ≤ z ≲ 3), luminosity (1040 erg s−1 ≤ LX ≤ 1046 erg s−1), and black hole (BH) mass (106 ≤ M⊙ ≤ 109). Methods. We propose the use of the variance-frequency diagram as a viable alternative to the study of the power spectral density (PSD), which is not yet accessible for distant, faint, and/or sparsely sampled AGN. Results. We show that the data collected from archival observations and previous literature studies are fully consistent with a universal PSD form, which does not show any evidence for systematic evolution of shape or amplitude with redshift or luminosity, even if there may be differences between individual AGN at a given redshift or luminosity. We find new evidence that the PSD bend frequency depends on BH mass and possibly on accretion rate. We finally discuss the implications for current and future AGN population and cosmological studies.

Ray-tracing Simulations and Spectral Models of X-Ray Radiation in Dusty Media

Dust can play an important role in shaping the X-ray spectra and images of astrophysical sources. In this work we report on the implementation of dust in the ray-tracing platform RefleX. We illustrate the different effects associated with the interaction between X-ray photons and dust grains, such as dust scattering, near-edge X-ray absorption fine structures, and shielding. We show how the cross sections of the photon–gas interaction change depending on the fraction of metals in dust grains (i.e., the dust depletion factor). We compare RefleX simulations to the most widely used absorption model that includes dust and show how X-ray spectra are affected by the presence of dust in the absorbing/reprocessing medium for different geometries. We also show how RefleX can be used to reproduce the dust scattering halos observed in Galactic sources, and we release the first torus X-ray spectral model that considers dust absorption and scattering (RXTorusD), to reproduce the spectra of active galactic nuclei (AGNs). RXTorusD also considers other physical processes that are not included in the most widely used AGN torus models, such as Rayleigh scattering and scattering on molecular gas, which can lead to remarkable differences in the predicted X-ray spectra for the same set of geometrical and physical parameters.

An empirical relation to estimate host galaxy stellar light from AGN spectra

ABSTRACT Measurement of black hole mass for low-z (z≤ 0.8) Active Galactic Nuclei (AGNs) is difficult due to the strong contribution from host galaxy stellar light necessitating detailed spectral decomposition to estimate the AGN luminosity. Here, we present an empirical relation to estimate host galaxy stellar luminosity from the optical spectra of AGNs at z ≤ 0.8. The spectral data were selected from the fourteenth data release of the Sloan Digital Sky Survey (SDSS-DR14) quasar catalogue having a signal-to-noise ratio at 5100 Å (SNR5100) >10 containing 11 415 quasars. The median total luminosity (log (Ltotal/[erg s−1])), stellar luminosity (log (Lstar/[erg s−1])), and AGN continuum luminosity ((log Lcont/[erg s−1])) in our sample are 44.52, 44.06, and 44.30, respectively. We fit the AGN power-law continuum, host galaxy, and iron blend contribution, simultaneously over the entire available spectrum. We found the host galaxy fraction to anticorrelate with continuum luminosity and can be well-represented by a polynomial function, which can be used to correct the stellar light contribution from AGN spectra. We also found anticorrelation between host galaxy fraction and iron strength, Eddington ratio, and redshift. The empirical relation gives comparable results of host-fraction with the image decomposition method.