Host Galaxy Components Research Articles

Prior-informed Active Galactic Nucleus Host Spectral Decomposition Using PyQSOFit

We introduce an improved method for decomposing the emission of active galactic nuclei (AGN) and their host galaxies using templates from principal component analysis (PCA). This approach integrates prior information from PCA with a penalized pixel fitting mechanism that improves the precision and effectiveness of the decomposition process. Specifically, we have reduced the degeneracy and overfitting in AGN host decomposition, particularly for those with low signal-to-noise ratios (SNRs), where traditional methods tend to fail. By applying our method to 76,565 Sloan Digital Sky Survey Data Release 16 quasars with z < 0.8, we achieve a success rate of ≈94%, thus establishing the largest host-decomposed spectral catalog of quasars to date. Our fitting results consider the impact of the host galaxy on the overestimation of the AGN luminosity and black hole mass (M BH). Furthermore, we obtained stellar velocity dispersion (σ ⋆) measurements for 4137 quasars. The slope of the M BH−σ ⋆ relation in this subsample is generally consistent with previous quasar studies beyond the local Universe. Our method provides a robust and efficient approach to disentangle the AGN and host galaxy components across a wide range of SNRs and redshifts.

SMART: Spectral energy distributions Markov chain Analysis with Radiative Transfer models

Abstract In this paper we present the publicly available open-source spectral energy distribution (SED) fitting code SMART (Spectral energy distributions Markov chain Analysis with Radiative Transfer models). Implementing a Bayesian Markov chain Monte Carlo (MCMC) method, SMART fits the ultraviolet to millimetre SEDs of galaxies exclusively with radiative transfer models that currently constitute four types of pre-computed libraries, which describe the starburst, active galactic nucleus (AGN) torus, host galaxy and polar dust components. An important novelty of SMART is that, although it fits SEDs exclusively with radiative transfer models, it takes comparable time to popular energy balance methods to run. Here we describe the key features of SMART and test it by fitting the multi-wavelength SEDs of the 42 local ultraluminous infrared galaxies (ULIRGs) that constitute the HERschel Ultraluminous Infrared Galaxy Survey (HERUS) sample. The Spitzer spectroscopy data of the HERUS ULIRGs are included in the fitting at a spectral resolution, which is matched to that of the radiative transfer models. We also present other results that highlight the performance and versatility of SMART. SMART promises to be a useful tool for studying galaxy evolution in the James Webb Space Telescope (JWST) era. SMART is developed in Python and is available at https://github.com/ch-var/SMART.git.

Coevolution and Nuclear Structure in the Dwarf Galaxy POX 52 Studied by Multiwavelength Data from Radio to X-Ray

The nearby dwarf galaxy POX 52 at z = 0.021 hosts an active galactic nucleus (AGN) with a black hole (BH) mass of M BH ∼ 105–6 M ⊙ and an Eddington ratio of ∼0.1–1. This object provides the rare opportunity to study both AGN and host-galaxy properties in a low-mass highly accreting system. To do so, we collected its multiwavelength data from X-ray to radio. First, we construct a spectral energy distribution, and by fitting it with AGN and host-galaxy components, we constrain AGN-disk and dust-torus components. Then, while considering the AGN-disk emission, we decompose optical Hubble Space Telescope images. As a result, it is found that a classical bulge component is probably present, and its mass (M bulge) is consistent with an expected value from a local relation. Lastly, we analyze new quasi-simultaneous X-ray (0.2–30 keV) data obtained by the Nuclear Spectroscopic Telescope Array and XMM-Newton. The X-ray spectrum can be reproduced by multicolor blackbody, warm and hot coronae, and disk and torus reflection components. Based on this, the spin is estimated to be a spin = 0.998−0.814, which could suggest that most of the current BH mass was achieved by prolonged mass accretion. Given the presence of the bulge, POX 52 would have undergone a galaxy merger, while the M BH–M bulge relation and the inferred prolonged accretion could suggest that AGN feedback occurred. Regarding the AGN structure, the spectral slope of the hot corona, its relative strength to the bolometric emission, and the torus structure are found to be consistent with Eddington-ratio dependencies found for nearby AGNs.

Testing super-eddington accretion on to a supermassive black hole: reverberation mapping of PG 1119+120

ABSTRACT We measure the black hole mass and investigate the accretion flow around the local (z = 0.0502) quasar PG 1119+120. Spectroscopic monitoring with Calar Alto provides H β lags and linewidths from which we estimate a black hole mass of log (M•/M⊙) = 7.0, uncertain by ∼0.4 dex. High cadence photometric monitoring over 2 yr with the Las Cumbres Observatory provides light curves in seven optical bands suitable for intensive continuum reverberation mapping. We identify variability on two time-scales. Slower variations on a 100-d time-scale exhibit excess flux and increased lag in the u′ band and are thus attributable to diffuse bound-free continuum emission from the broad-line region. Faster variations that we attribute to accretion disc reprocessing lack a u′-band excess and have flux and delay spectra consistent with either τ ∝ λ4/3, as expected for a temperature structure of T(R) ∝ R−3/4 for a thin accretion disc, or τ ∝ λ2 expected for a slim disc. Decomposing the flux into variable (disc) and constant (host galaxy) components, we find the disc SED to be flatter than expected with $f_{\nu } \sim \rm {const}$. Modelling the SED predicts an Eddington ratio of λEdd &amp;gt; 1, where the flat spectrum can be reproduced by a slim disc with little dust extinction or a thin disc that requires more dust extinction. While this accretion is super-Eddington, the geometry is still unclear; however, a slim disc is expected due to the high radiation pressure at these accretion rates, and is entirely consistent with our observations.

Molecular Gas Heating, Star Formation Rate Relations, and AGN Feedback in Infrared-Luminous Galaxy Mergers

We examine the origin of molecular gas heating in a sample of 42 infrared-luminous galaxies at z&lt;0.3 by combining two sets of archival data: first, integrated CO line luminosities in the 1–0 and 5–4 through 13–12 transitions; second, results from radiative transfer modelling that decompose their bolometric emission into starburst, AGN, and host galaxy components. We find that the CO 1–0 and 5–4 through 9–8 lines primarily arise via radiative heating in the starburst and the host galaxy. In contrast, the CO 10–9 through 13–12 lines may arise primarily in the starburst and AGN, with an increasing contribution from mechanical heating and shocks. For the sample as a whole, we find no evidence that AGN luminosity affects the heating of molecular gas by star formation. However, for starbursts with low initial optical depths, a more luminous AGN may reduce the efficiency of starburst heating of the CO 5–4 and above lines, consistent with negative AGN feedback.

Dust and the intrinsic spectral index of quasar variations: hints of finite stress at the innermost stable circular orbit

ABSTRACT We present a study of 9 242 spectroscopically confirmed quasars with multiepoch ugriz photometry from the SDSS Southern Survey. By fitting a separable linear model to each quasar’s spectral variations, we decompose their five-band spectral energy distributions into variable (disc) and non-variable (host galaxy) components. In modelling the disc spectra, we include attenuation by dust on the line of sight through the host galaxy to its nucleus. We consider five commonly used attenuation laws, and find that the best description is by dust similar to that of the Small Magellanic Cloud, inferring a lack of carbonaceous grains from the relatively weak 2175-Å absorption feature. We go on to construct a composite spectrum for the quasar variations spanning 700–8000 Å. By varying the assumed power-law Lν ∝ να spectral slope, we find a best-fitting value α = 0.71 ± 0.02, excluding at high confidence the canonical Lν ∝ ν1/3 prediction for a steady-state accretion disc with a T ∝ r−3/4 temperature profile. The bluer spectral index of the observed quasar variations instead supports the model of Agol &amp; Krolik, and Mummery &amp; Balbus, in which a steeper temperature profile, T ∝ r−7/8, develops as a result of finite magnetically induced stress at the innermost stable circular orbit extracting energy and angular momentum from the black hole spin.

Subaru Hyper Suprime-Cam view of quasar host galaxies at z &amp;lt; 1

Abstract Active galactic nuclei (AGNs) are key in understanding the coevolution of galaxies and supermassive black holes (SMBHs). AGN activity is thought to affect the properties of their host galaxies via a process called “AGN feedback,” which drives the coevolution. From a parent sample of 1151 z &amp;lt; 1 type-1 quasars from the Sloan Digital Sky Survey quasar catalog, we detected the host galaxies of 862 of them in the high-quality grizy images of the Subaru Hyper Suprime-Cam survey. The unprecedented combination of the survey area and depth allowed us to perform a statistical analysis of the quasar host galaxies, with small sample variance. We fitted the radial image profile of each quasar as a linear combination of the point spread function and the Sérsic function, decomposing the images into the quasar nucleus and the host galaxy components. We found that the host galaxies are massive, with stellar mass Mstar ≳ 1010 $M_\odot$, and are mainly located on the green valley. This trend is consistent with a scenario in which star formation in the host galaxies is suppressed by AGN feedback, that is, AGN activity may be responsible for the transition of these galaxies from the blue cloud to the red sequence. We also investigated the SMBH mass to stellar mass relation of the z &amp;lt; 1 quasars, and found a consistent slope with the local relation, while the SMBHs may be slightly undermassive. However, the above results are subject to our sample selection, which biases against host galaxies with low masses and/or large quasar-to-host flux ratios.

Ionization Mechanisms in Quasar Outflows

Abstract The various ionization mechanisms at play in active galactic nuclei (AGN) and quasars have been well studied, but relatively little has been done to separately investigate the contributions of these ionization mechanisms within the host galaxy and outflowing components. Using Gemini integral field spectroscopy data, we study the ionization properties of these two components in four nearby (z ≲ 0.2) radio-quiet Type 1 quasars. Emission-line ratios and widths are employed to identify the dominant ionization mechanisms for the host and outflow components in each object. We find that photoionization by the AGN often dominates the ionization of both gaseous components in these systems. In three cases, the outflowing gas is more highly ionized than the gas in the host, indicating that it is more strongly exposed to the ionizing radiation field of the AGN. In two objects, a positive correlation between the line widths and line ratios in the outflowing gas component indicates that shocks with velocities of order 100–500 km s−1 may also be contributing to the ionization and heating of the outflowing gas component. The line ratios in the outflowing gas of one of these two objects also suggest a significant contribution from photoionization by hot, young stars in the portion of the outflow that is closest to star-forming regions in the host galaxy component. The data thus favor photoionization by hot stars in the host galaxy rather than stars formed in the outflow itself.

Minute-scale Rapid Variability of Mrk 501 by Multi-band Photometric Monitoring from 2010 to 2017

Mrk 501 is one of our long-term photometric monitoring objects with the optical telescopes at Yunnan Astronomical Observatory (YAO). In this work, observations were performed from 2010 March 11 to 2017 June 7 in the BVRI bands. We detected minute-scale rapid variabilities in the optical multi-band light curves. Durations of these rapid variabilities are from 5 to 77 minutes. The minimal decay timescale is about 3 minutes (2.330 ± 1.572 minutes on 2012 May 22 and 3.172 ± 1.626 minutes on 2017 May 20). After subtracting the host-galaxy component, The minimal decay timescale was about 8 minutes on 2017 June 3, and the corresponding duration is about 10 minutes. This fast variability requires the strength of the co-moving magnetic field in the jet larger than 38.60 G (19.69 G for corrected results, with Doppler factor δD ∼ 3). The discrete correlation function (DCF) analysis between B and I band short-term light curves show significant correlated variability; however, no significant time lags are detected. The color–magnitude variation shows significant negative correlation on intranight timescales, which is due to the superposition of multiple components.

AGN mass estimates in large spectroscopic surveys: the effect of host galaxy light

Virial–based methods for estimating active supermassive black hole masses are now commonly used on extremely large spectroscopic quasar catalogs. Most spectral analyses, though, do not pay enough attention to the detailed continuum decomposition. To understand how this affects virial mass estimates, we test the influence of host galaxy light on them, along with a Balmer continuum component. A detailed fit with the new spectroscopic analysis software QSFIT demonstrates that the presence or absence of continuum components does not significantly affect the virial-based results for our sample. Taking a host galaxy component into consideration or not, instead, affects the emission line fitting in a more pronounced way at lower redshifts, where in fact we observe dimmer quasars and more visible host galaxies.

A Catalog of Post-starburst Quasars from Sloan Digital Sky Survey Data Release 7

Abstract We present a catalog of nearby (z ≤ 0.5) quasars with significant features of post-starburst stellar populations in their optical spectra: so-called post-starburst quasars, or PSQs. After carefully decomposing spectra from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) Quasar Catalog into quasar and host-galaxy components, we derive a sample of 208 PSQs. Their host-galaxy components have strong Hδ absorption (EW ≥ 6 Å) indicating a significant contribution of an intermediate-aged stellar population formed in a burst of star formation within the past 1 Gyr, which makes them potentially useful for studying the co-evolution of supermassive black holes and their host galaxies.

A Hubble Space Telescope imaging study of four FeLoBAL quasar host galaxies

We study the host galaxies of four Iron Low-Ionization Broad Absorption-line Quasars (FeLoBALs) using Hubble Space Telescope imaging data, investigating the possibility that they represent a transition between an obscured AGN and an ordinary optical quasar. In this scenario, the FeLoBALs represent the early stage of merger-triggered accretion, in which case their host galaxies are expected to show signs of an ongoing or recent merger. Using PSF subtraction techniques, we decompose the images into host galaxy and AGN components at rest-frame ultraviolet and optical wavelengths. The ultraviolet is sensitive to young stars, while the optical probes stellar mass. In the ultraviolet we image at the BAL absorption trough wavelengths so as to decrease the contrast between the quasar and host galaxy emission. We securely detect an extended source for two of the four FeLoBALs in the rest-frame optical; a third host galaxy is marginally detected. In the rest-frame UV we detect no host emission; this constrains the level of unobscured star formation. Thus, the host galaxies have observed properties that are consistent with those of non-BAL quasars with the same nuclear luminosity, i.e., quiescent or moderately starforming elliptical galaxies. However, we cannot exclude starbursting hosts that have the stellar UV emission obscured by modest amounts of dust reddening. Thus, our findings also allow the merger-induced young quasar scenario. For three objects, we identify possible close companion galaxies that may be gravitationally interacting with the quasar hosts.

MULTI-COLOR OPTICAL MONITORING OF MRK 501 FROM 2010 TO 2015

ABSTRACT We have monitored the BL Lac object Mrk 501 in the optical V, R, and I bands from 2010 to 2015. For Mrk 501, the presence of a strong host galaxy component can affect the results of photometry. After subtracting the host galaxy contributions, the source shows intraday and long-term variabilities for optical flux and color indices. The average variability amplitudes of the V, R, and I bands are 22.05 &percnt; , 22.25 &percnt; , and 23.82 &percnt; ?> , respectively, and the value of the duty cycle is 14.87%. A minimal variability timescale of 106 minutes is detected. No significant time lag between the V and I bands is found on one night. The bluer-when-brighter (BWB) trend is dominant for Mrk 501 on intermediate, short, and intraday timescales, which supports the shock-in-jet model. For the long timescale, Mrk 501, in different states, can have different BWB trends. The corresponding results for non-correcting host galaxy contributions are also presented.

A CONNECTION BETWEEN OBSCURATION AND STAR FORMATION IN LUMINOUS QUASARS

We present a measurement of the star formation properties of a uniform sample of mid-IR selected, unobscured and obscured quasars (QSO1s and QSO2s) in the Bo\"otes survey region. We use an spectral energy distribution (SED) analysis for photometric data spanning optical to far-IR wavelengths to decompose AGN and host galaxy components. We find that when compared to a matched sample of QSO1s, the QSO2s have higher far-IR detection fractions, far-IR fluxes and infrared star formation luminosities ($L_{\rm IR}^{\rm SF}$) by a factor of $\sim2$. Correspondingly, we show that the AGN obscured fraction rises from 0.3 to 0.7 between $4-40\times10^{11}L_\odot$. We also find evidence associating the absorption in the X-ray emission with the presence of far-IR emitting dust. Overall, these results are consistent with galaxy evolution models in which quasar obscurations can be associated with a dust-enshrouded starburst galaxies.

The ensemble optical variability of type-1 AGN in the Sloan Digital Sky Survey Data Release 7

We use a sample of over 5000 active galactic nuclei (AGN) with extended morphologies at z < 0.8 from the Sloan Digital Sky Survey (SDSS) to study the ensemble optical variability as a function of rest-frame time lag and AGN luminosity with the aim of investigating these parameter relationships at lower luminosities than previously studied. We compare photometry from imaging data with spectrophotometry obtained weeks to years later in the Sloan g, r, and i bands. We employ quasar and galaxy eigenspectra fitting to separate the AGN and host galaxy components. A strong correlation between the variability amplitude and rest-frame time lag is observed, in agreement with quasar structure functions but extending to AGN several magnitudes fainter than previously studied. The structure function slopes for our fainter AGN sample are slightly shallower than those found in quasars studies. An anticorrelation with luminosity is clearly detected, with lower luminosity AGN displaying greater variability amplitudes. We demonstrate for the first time that this anticorrelation extends to AGN as faint as MAGNi ∼ -18.5, with a slight trend towards shallower slopes at luminosities fainter than MAGNi ∼ -20.2.

A TALE OF TWO FEEDBACKS: STAR FORMATION IN THE HOST GALAXIES OF RADIO AGNs

Several lines of argument support the existence of a link between activity at the nuclei of galaxies, in the form of an accreting supermassive black hole, and star-formation activity in these galaxies. The exact nature of this link is still under debate. Radio jets have long been argued to be an ideal mechanism that allows AGN to interact with their host galaxy and regulate star-formation. In this context, we are using a sample of radio sources in the North Ecliptic Pole (NEP) field to study the nature of the putative link between AGN activity and star-formation. This is done by means of spectral energy distribution (SED) fitting. We use the excellent spectral coverage of the AKARI infrared space telescope together with the rich ancillary data available in the NEP to build SEDs extending from UV to far-IR wavelengths. Through SED fitting we constrain both the AGN and host galaxy components. We find a significant AGN component in our sample of relatively faint radio-sources ($<$mJy), that increases in power with increasing radio-luminosity. At the highest radio-luminosities, the presence of powerful jets dominates the radio emission of these sources. A positive correlation is found between the luminosity of the AGN component and that of star-formation in the host galaxy, independent of the radio luminosity. By contrast, for a given redshift and AGN luminosity, we find that increasing radio-luminosity leads to a decrease in the specific star-formation rate. The most radio-loud AGN are found to lie on the main sequence of star-formation for their respective redshifts. For the first time, such a two-sided feedback process is seen in the same sample. We conclude that radio jets do suppress star-formation in their host galaxies but appear not to totally quench it. Our results therefore support the maintenance nature of "radio-mode" feedback from radio-AGN jets.

Through the kaleidoscope: star formation the host galaxies of radio-AGN

AbstractThere exist strong evidence supporting the co-evolution of central supermassive black holes and their host galaxies; however it is still under debate how such a relation comes about and whether it is relevant for all or only a subset of galaxies. A rich multi-wavelength dataset is available for the North Ecliptic Pole field, most notably surveyed by the AKARI infrared space telescope. We investigate the star-formation properties of the host galaxies of radio-AGN together with the radio feedback mechanism, potentially responsible for the eventual quenching of star formation. Using broadband SED modelling, the nuclear and host galaxy components of these sources are studied as a function of their radio luminosity. Here we present results concerning the AGN content of the radio sources in this field, while offering evidence supporting a “maintenance” type of feedback from powerful radio-jets.

LOW-RESOLUTION SPECTRAL TEMPLATES FOR ACTIVE GALACTIC NUCLEI AND GALAXIES FROM 0.03 TO 30 μm

We present a set of low resolution empirical SED templates for AGNs and galaxies in the wavelength range from 0.03 to 30 microns based on the multi-wavelength photometric observations of the NOAO Deep-Wide Field Survey Bootes field and the spectroscopic observations of the AGN and Galaxy Evolution Survey. Our training sample is comprised of 14448 galaxies in the redshift range 0<~z<~1 and 5347 likely AGNs in the range 0<~z<~5.58. We use our templates to determine photometric redshifts for galaxies and AGNs. While they are relatively accurate for galaxies, their accuracies for AGNs are a strong function of the luminosity ratio between the AGN and galaxy components. Somewhat surprisingly, the relative luminosities of the AGN and its host are well determined even when the photometric redshift is significantly in error. We also use our templates to study the mid-IR AGN selection criteria developed by Stern et al.(2005) and Lacy et al.(2004). We find that the Stern et al.(2005) criteria suffers from significant incompleteness when there is a strong host galaxy component and at z =~ 4.5, when the broad Halpha emission line is redshifted into the [3.6] band, but that it is little contaminated by low and intermediate redshift galaxies. The Lacy et al.(2004) criterion is not affected by incompleteness at z =~ 4.5 and is somewhat less affected by strong galaxy host components, but is heavily contaminated by low redshift star forming galaxies. Finally, we use our templates to predict the color-color distribution of sources in the upcoming WISE mission and define a color criterion to select AGNs analogous to those developed for IRAC photometry. We estimate that in between 640,000 and 1,700,000 AGNs will be identified by these criteria, but will have serious completeness problems for z >~ 3.4.

GEMINI 3D spectroscopy of BAL + IR + Fe ii QSOs - I. Decoupling the BAL, QSO, starburst, NLR, supergiant bubbles and galactic wind in Mrk 231

In this work we present the first results of a study of BAL QSOs (at low and high redshift), based on very deep Gemini GMOS integral field spectroscopy. In particular, the results obtained for the nearest BAL IR QSO Mrk 231 are presented. Very deep three-dimensional spectra and maps clearly show that the BAL systems I and II are extended (reaching 1.4-1.6" = 1.2-1.3 kpc, from the nucleus) and clearly elongated at the position angle close to the radio jet PA. Which suggest that the BAL systems I and II are both associated with the radio jet, and supporting the bipolar jet-wind model for some BALs. For the nuclear region of Mrk 231, the QSO and host-galaxy components were modelled, using a new technique of decoupling 3D spectra. From this study, the following main results were found: (i) in the pure host galaxy spectrum an extreme nuclear starburst component was clearly observed, mainly as a very strong increase in the flux, at the blue wavelengths; (ii) the BAL system I is observed in the spectrum of the host galaxy; (iii) in the clean/pure QSO emission spectrum, only broad lines were detected. 3D GMOS individual spectra (specially in the IR Ca II triplet) and maps confirm the presence of an extreme and young nuclear starburst (8 < age < 15 Myr), which was detected mainly in a ring or toroid with a radius r = 0.3" = 200 pc, around the very nucleus. The physical properties of the four expanding nuclear bubbles were analysed, using the GMOS 3D spectra and maps. These results suggest that an important part of the nuclear NLR is generated by the OF process and the associated low velocity ionizing shocks.

An extreme EXO: a type 2 QSO atz = 1.87

We aim to understand the multi-wavelength properties of 2XMM J123204+215255, the source with the most extreme X-ray-to-optical flux ratio amongst a sample of bright X-ray selected EXOs drawn from a cross-correlation of the 2XMMp catalogue with the SDSS-DR5 catalogue. We use 2XMMp X-ray data, SDSS-DR5, NOT and UKIRT optical/NIR photometric data and Subaru MOIRCS IR spectroscopy to study the properties of 2XMM J123204+215255. We created a model SED including an obscured QSO and the host galaxy component to constrain the optical/IR extinction and the relative contribution of the AGN and the galaxy to the total emission. 2XMM J123204+215255 is a bright X-ray source with f_X~10^{-12} erg cm^{-2} s^{-1} (2-10 keV energy band) which has no detection down to a magnitude i' > 25.2. NIR imaging reveals a faint K-band counterpart and NIR spectroscopy shows a single broad (FWHM=5300 km/s) emission line, which is almost certainly H-alpha at z=1.87. The X-ray spectrum shows evidence of significant absorption (N_H > 10^{23} cm^{-2}), typical of type 2 AGN, but the broad H-alpha emission suggests a type 1 AGN classification. The very red optical/NIR colours (i'-K > 5.3) strongly suggest significant reddening however. We find that simple modelling can successfully reproduce the NIR continuum and strongly constrain the intrinsic nuclear optical/IR extinction to A_V~4, which turns out to be much smaller than the expected from the X-ray absorption (assuming Galactic gas-to-dust ratio).