Spectral Energy Distribution Decomposition Research Articles

Radio-to-submillimetre spectral energy distributions of NGC 1365

ABSTRACT We analyse the radio-to-submillimetre spectral energy distribution (SED) for the central pseudo-bulge of NGC 1365 using archival data from the Atacama Large Millimeter/submillimeter Array (ALMA) and the Very Large Array. This analysis shows that free–free emission dominates the continuum emission at 50–120 GHz and produces about 75 per cent of the 103 GHz continuum emission. However, the fraction of 103 GHz continuum emission originating from free–free emission varies significantly among different subregions in the pseudo-bulge, particularly for an outflow from the active galactic nuclei (AGN) on the eastern pseudo-bulge where the synchrotron emission produces half of the 103 GHz continuum emission. Free–free emission also dominates at 103 GHz within the central 400 pc diameter region, but this emission is associated with the AGN rather than star formation. The star formation rate (SFR) within the pseudo-bulge derived from the ALMA free–free emission is 8.9 ± 1.1 M⊙ yr−1. This is comparable to the SFR from the mid-infrared emission but higher than the SFR from the extinction-corrected H α line emission, mainly because the pseudo-bulge is heavily dust obscured. The 1.5 GHz emission yields a comparable SFR for the pseudo-bulge but may have lower SFRs within subregions of the pseudo-bulge because of the diffusion outside of these regions of the electrons producing the synchrotron radiation. We propose that applying a correction factor of 75 per cent to the 80–110 GHz continuum emission could provide valuable estimates of the free–free emission without performing any SED decomposition, which could derive extinction-free SFRs within 20 per cent accuracy.

Dust-enshrouded AGNs Can Dominate Host-galaxy-scale Cold Dust Emission

Abstract It is widely assumed that long-wavelength infrared (IR) emission from cold dust (T ∼ 20–40 K) is a reliable tracer of star formation even in the presence of a bright active galactic nucleus (AGN). Based on radiative transfer (RT) models of clumpy AGN tori, hot dust emission from the torus contributes negligibly to the galaxy spectral energy distribution (SED) at λ ≳ 100 μm. However, these models do not include AGN heating of host-galaxy-scale diffuse dust, which may have far-IR (FIR) colors comparable to cold diffuse dust heated by stars. To quantify the contribution of AGN heating to host-galaxy-scale cold dust emission at λ ≳ 100 μm, we perform dust RT calculations on a simulated galaxy merger both including and excluding the bright AGN that it hosts. By differencing the SEDs yielded by RT calculations with and without AGNs that are otherwise identical, we quantify the FIR cold dust emission arising solely from reprocessed AGN photons. In extreme cases, AGN-heated host-galaxy-scale dust can increase galaxy-integrated FIR flux densities by factors of 2–4; star formation rates calculated from the FIR luminosity assuming no AGN contribution can overestimate the true value by comparable factors. Because the FIR colors of such systems are similar to those of purely star-forming galaxies and redder than torus models, broadband SED decomposition may be insufficient for disentangling the contributions of stars and heavily dust-enshrouded AGNs in the most IR-luminous galaxies. We demonstrate how kiloparsec-scale resolved observations can be used to identify deeply dust-enshrouded AGNs with cool FIR colors when spectroscopic and/or X-ray detection methods are unavailable.

How Does the Polar Dust Affect the Correlation between Dust Covering Factor and Eddington Ratio in Type 1 Quasars Selected from the Sloan Digital Sky Survey Data Release 16?

Abstract We revisit the dependence of the covering factor (CF) of dust torus on physical properties of active galactic nuclei (AGNs) by taking into account an AGN polar dust emission. The CF is converted from a ratio of infrared (IR) luminosity contributed from AGN dust torus ( ) and AGN bolometric luminosity (L bol), by assuming a nonlinear relation between luminosity ratio and intrinsic CF. We select 37,181 type 1 quasars at z < 0.7 from the Sloan Digital Sky Survey Data Release 16 quasar catalog. Their L bol, black hole mass (M BH), and Eddington ratio (λ Edd) are derived by spectral fitting with QSFit. We conduct spectral energy distribution decomposition by using X-CIGALE with a clumpy torus and polar dust model to estimate without being affected by the contribution of stellar and AGN polar dust to IR emission. For 5752 quasars whose physical quantities are securely determined, we perform a correlation analysis on CF and (i) L bol, (ii) M BH, and (iii) λ Edd. As a result, anticorrelations for CF–L bol, CF–M BH, and CF–λ Edd are confirmed. We find that incorporating the AGN polar dust emission makes those anticorrelations stronger compared to those without considering it. This indicates that polar dust wind probably driven by AGN radiative pressure is one of the key components to regulate obscuring material of AGNs.

X-ray flux in SED modelling: An application of X-CIGALE in the XMM-XXL field

X-CIGALE is built on the spectral energy distribution (SED) code of CIGALE and implements important new features: the code accounts for obscuring material in the poles of AGNs and has the ability to fit X-ray fluxes. In this work, we used ∼2500 spectroscopic X-ray AGNs from the XMM-XXL-north field and examined the improvements the new features bring to the SED modelling analysis. Based on our results, X-CIGALE successfully connects the X-ray with the UV luminosity in the whole range spanned by our sample (log LX(2−10 keV) = (42−46) erg s−1). The addition of the new features generally improves the efficiency of X-CIGALE in the estimation and characterisation of the AGN component. Classification as type 1 or type 2 based on their inclination angle is improved, especially at redshifts lower than 1. The statistical significance of AGN fraction, fracAGN, measurements is increased, in particular for luminous X-ray sources (LX > 1045 erg s−1). These conclusions hold under the condition that (mid-) IR photometry is available in the SED fitting process. The addition of polar dust increases the AGN fraction and the efficiency of the SED decomposition to detect AGNs among X-ray selected sources. X-CIGALE estimates a strong AGN (fracAGN > 0.3) in more than 90% of the IR-selected AGNs and 75% of X-ray-detected AGNs not selected by IR colour criteria. The latter drops to ∼50% when polar dust is not included. The ability of X-CIGALE to include X-ray information in the SED fitting process can be instrumental in the optimal exploitation of the wealth of data that current (eROSITA) and future (ATHENA) missions will provide us.

The role of SPICA-like missions and the Origins Space Telescope in the quest for heavily obscured AGN and synergies with Athena

Abstract In the black hole (BH)–galaxy co-evolution framework, most of the star formation (SF) and the BH accretion are expected to take place in highly obscured conditions. The large amount of gas and dust absorbs most of the UV-to-soft-X radiation and re-emits it at longer wavelengths, mostly in the IR. Thus, obscured active galactic nuclei (AGN) are very difficult to identify in optical or X-ray bands but shine bright in the IR. Moreover, X-ray background (XRB) synthesis models predict that a large fraction of the yet-unresolved XRB is due to the most obscured (Compton thick, CT: N $_{\text{H}}\ge 10^{24} \,\mathrm{cm}^{-2}$ ) of these AGN. In this work, we investigate the synergies between putative IR missions [using SPace Infrared telescope for Cosmology and Astrophysics (SPICA), proposed for European Space Agency (ESA)/M5 but withdrawn in 2020 October, and Origins Space Telescope, OST, as ‘templates’] and the X-ray mission Athena (Advanced Telescope for High ENergy Astrophysics), which should fly in early 2030s, in detecting and characterising AGN, with a particular focus on the most obscured ones. Using an XRB synthesis model, we estimated the number of AGN and the number of those which will be detected in the X-rays by Athena. For each AGN, we associated an optical-to-Far InfraRed (FIR) spectral energy distribution (SED) from observed AGN with both X-ray data and SED decomposition and used these SEDs to check if the AGN will be detected by SPICA-like or OST at IR wavelengths. We expect that, with the deepest Athena and SPICA-like (or OST) surveys, we will be able to photometrically detect in the IR more than 90% of all the AGN (down to $L_{2-10\text{keV}} \sim 10^{42}\,\mathrm{erg\ s}^{-1}$ and up to $z \sim 10$ ) predicted by XRB synthesis modeling, and we will detect at least half of them in the X-rays. The spectroscopic capabilities of the OST can provide ${\approx}51\,000$ and ${\approx}3\,400$ AGN spectra with $R= 300$ at 25–588 $\unicode[Times]{x03BC}$ m in the wide and deep surveys, respectively, the last one up to $z\approx 4$ . Athena will be extremely powerful in detecting and discerning moderate- and high-luminosity AGN, allowing us to properly select AGN even when the mid-IR torus emission is ‘hidden’ by the host galaxy contribution. We will constrain the intrinsic luminosity and the amount of obscuration for $\sim\!20\%$ of all the AGN (and $\sim\!50\%$ of those with $L_{2-10\text{keV}} > 3.2 \times 10^{43}\,\mathrm{erg\ s}^{-1}$ ) using the X-ray spectra provided by Athena WFI. We find that the most obscured and elusive CT-AGN will be exquisitely sampled by SPICA-like mission or OST and that Athena will allow a fine characterisation of the most luminous ones. This will provide a significant step forward in the process of placing stronger constraints on the yet-unresolved XRB and investigating the BH accretion rate evolution up to very high redshift ( $z \ge 4$ ).

An obscured AGN population hidden in the VIPERS galaxies: identification through spectral energy distribution decomposition

ABSTRACT The detection of X-ray emission constitutes a reliable and efficient tool for the selection of active galactic nuclei (AGNs), although it may be biased against the most heavily absorbed AGNs. Simple mid-infrared (IR) broad-band selection criteria identify a large number of luminous and absorbed AGNs, yet again host contamination could lead to non-uniform and incomplete samples. Spectral energy distribution (SED) decomposition is able to decouple the emission from the AGN versus that from star-forming regions, revealing weaker AGN components. We aim to identify the obscured AGN population in the VIMOS Public Extragalactic Redshift Survey in the Canada–France–Hawaii Telescope Legacy Survey W1 field through SED modelling. We construct SEDs for 6860 sources and identify 160 AGNs at a high confidence level using a Bayesian approach. Using optical spectroscopy, we confirm the nature of ∼85 per cent of the AGNs. Our AGN sample is highly complete (∼92 per cent) compared to mid-IR colour-selected AGNs, including a significant number of galaxy-dominated systems with lower luminosities. In addition to the lack of X-ray emission (80 per cent), the SED fitting results suggest that the majority of the sources are obscured. We use a number of diagnostic criteria in the optical, IR, and X-ray regimes to verify these results. Interestingly, only 35 per cent of the most luminous mid-IR-selected AGNs have X-ray counterparts suggesting strong absorption. Our work emphasizes the importance of using SED decomposition techniques to select a population of type II AGNs, which may remain undetected by either X-ray or IR colour surveys.

AGNs Are Not That Cool: Revisiting the Intrinsic AGN Far-infrared Spectral Energy Distribution

Abstract We investigate the intrinsic spectral energy distribution (SED) of active galactic nuclei (AGNs) at infrared (IR) bands with 42 z < 0.5 optically luminous Palomar Green survey quasars through SED decomposition. We decompose the SEDs of the 42 quasars by combining an AGN IR template library that covers a wide range of the AGN parameter space with three commonly used galaxy template libraries. We determine the median AGN SED from the best-fitting results. The far-IR (FIR) contribution of our median AGN SED is significantly smaller than that of Symeonidis et al., but roughly consistent with that of Lyu & Rieke. The AGN IR SED becomes cooler with increasing bolometric luminosity, which might be due to the fact that more luminous AGNs might have stronger radiative feedback to change torus structures and/or their tori might have higher metallicities. Our conclusions do not depend on the choice of galaxy template libraries. However, as the predicted polycyclic aromatic hydrocarbon (PAH) emission line flux is galaxy template dependent, cautions should be taken on deriving galaxy FIR contribution from PAH fluxes.

No Significant Evolution of Relations between Black Hole Mass and Galaxy Total Stellar Mass Up to z ∼ 2.5

Abstract We investigate the cosmic evolution of the ratio between black hole (BH) mass (M BH) and host galaxy total stellar mass (M stellar) out to z ∼ 2.5 for a sample of 100 X-ray-selected moderate-luminosity, broad-line active galactic nuclei (AGNs) in the Chandra-COSMOS Legacy Survey. By taking advantage of the deep multiwavelength photometry and spectroscopy in the COSMOS field, we measure in a uniform way the galaxy total stellar mass using an spectral energy distribution decomposition technique and the BH mass based on broad emission line measurements and single-epoch virial estimates. Our sample of AGN host galaxies has total stellar masses of 1010−12 M ⊙, and BH masses of 107.0–9.5 M ⊙. Combining our sample with the relatively bright AGN samples from the literature, we find no significant evolution of the M BH–M stellar relation with the BH-to-host total stellar mass ratio of M BH/M stellar ∼ 0.3% at all redshifts probed. We conclude that the average BH-to-host stellar mass ratio appears to be consistent with the local value within the uncertainties, suggesting a lack of evolution of the M BH–M stellar relation up to z ∼ 2.5.

X-ray emission of z > 2.5 active galactic nuclei can be obscured by their host galaxies

We present a multiwavelength study of seven active galactic nuclei (AGN) at spectroscopic redshift >2.5 in the 7 Ms Chandra Deep Field South that were selected for their good far-infrared (FIR) and submillimeter (submm) detections. Our aim is to investigate the possibility that the obscuration observed in the X-rays can be produced by the interstellar medium (ISM) of the host galaxy. Based on the 7 Ms Chandra spectra, we measured obscuring column densities NH, X in excess of 7 × 1022 cm−2 and intrinsic X-ray luminosities LX > 1044 erg s−1 for our targets, as well as equivalent widths for the Fe Kα emission line EWrest ≳ 0.5−1 keV. We built the UV-to-FIR spectral energy distributions (SEDs) by using broadband photometry from the CANDELS and Herschel catalogs. By means of an SED decomposition technique, we derived stellar masses (M* ∼ 1011 M⊙), IR luminosities (LIR > 1012 L⊙), star formation rates (SFR ∼ 190−1680 M⊙ yr−1) and AGN bolometric luminosities (Lbol ∼ 1046 erg s−1) for our sample. We used an empirically calibrated relation between gas masses and FIR/submm luminosities and derived Mgas ∼ 0.8−5.4 × 1010 M⊙. High-resolution (0.3−0.7″) ALMA data (when available, CANDELS data otherwise) were used to estimate the galaxy size and hence the volume enclosing most of the ISM under simple geometrical assumptions. These measurements were then combined to derive the column density associated with the ISM of the host, which is on the order of NH, ISM ∼ 1023−24 cm−2. The comparison between the ISM column densities and those measured from the X-ray spectral analysis shows that they are similar. This suggests that at least at high redshift, significant absorption on kiloparsec scales by the dense ISM in the host likely adds to or substitutes that produced by circumnuclear gas on parsec scales (i.e., the torus of unified models). The lack of unobscured AGN among our ISM-rich targets supports this scenario.

BAT AGN Spectroscopic Survey. XI. The Covering Factor of Dust and Gas in Swift/BAT Active Galactic Nuclei

Abstract We quantify the luminosity contribution of active galactic nuclei (AGNs) to the 12 μm, mid-infrared (MIR; 5–38 μm), and total IR (5–1000 μm) emission in the local AGNs detected in the all-sky 70 month Swift/Burst Alert Telescope (BAT) ultrahard X-ray survey. We decompose the IR spectral energy distributions (SEDs) of 587 objects into the AGN and starburst components using templates for an AGN torus and a star-forming galaxy. This enables us to recover the emission from the AGN torus including the low-luminosity end, down to , which typically has significant host galaxy contamination. The sample demonstrates that the luminosity contribution of the AGN to the 12 μm, the MIR, and the total IR bands is an increasing function of the 14–150 keV luminosity. We also find that for the most extreme cases, the IR pure-AGN emission from the torus can extend up to 90 μm. The total IR AGN luminosity obtained through the IR SED decomposition enables us to estimate the fraction of the sky obscured by dust, i.e., the dust covering factor. We demonstrate that the median dust covering factor is always smaller than the median X-ray obscuration fraction above an AGN bolometric luminosity of . Considering that the X-ray obscuration fraction is equivalent to the covering factor coming from both the dust and gas, this indicates that an additional neutral gas component, along with the dusty torus, is responsible for the absorption of X-ray emission.

Dust-deficient Palomar-Green Quasars and the Diversity of AGN Intrinsic IR Emission

Abstract To elucidate the intrinsic broadband infrared (IR) emission properties of active galactic nuclei (AGNs), we analyze the spectral energy distributions (SEDs) of 87 z ≲ 0.5 Palomar-Green (PG) quasars. While the Elvis AGN template with a moderate far-IR correction can reasonably match the SEDs of the AGN components in ∼60% of the sample (and is superior to alternatives such as that by Assef), it fails on two quasar populations: (1) hot-dust-deficient (HDD) quasars that show very weak emission thoroughly from the near-IR to the far-IR, and (2) warm-dust-deficient (WDD) quasars that have similar hot dust emission as normal quasars but are relatively faint in the mid- and far-IR. After building composite AGN templates for these dust-deficient quasars, we successfully fit the 0.3–500 μm SEDs of the PG sample with the appropriate AGN template, an infrared template of a star-forming galaxy, and a host galaxy stellar template. 20 HDD and 12 WDD quasars are identified from the SED decomposition, including seven ambiguous cases. Compared with normal quasars, the HDD quasars have AGNs with relatively low Eddington ratios and the fraction of WDD quasars increases with AGN luminosity. Moreover, both the HDD and WDD quasar populations show relatively stronger mid-IR silicate emission. Virtually identical SED properties are also found in some quasars from z = 0.5 to 6. We propose a conceptual model to demonstrate that the observed dust deficiency of quasars can result from a change of structures of the circumnuclear tori that can occur at any cosmic epoch.

Tracing black hole accretion with SED decomposition and IR lines: from local galaxies to the high-zUniverse

We present new estimates of AGN accretion and star-formation luminosity in galaxies obtained for the local 12-$\mu$m sample of Seyfert galaxies (12MGS), by performing a detailed broad-band spectral energy distribution (SED) decomposition including the emission of stars, dust heated by star formation and a possible AGN dusty torus. Thanks to the availability of data from the X-rays to the sub-millimetre, we constrain and test the contribution of the stellar, AGN and star-formation components to the SEDs. The availability of Spitzer-IRS low resolution mid-infrared (mid-IR) spectra is crucial to constrain the dusty torus component at its peak wavelengths. The results of SED-fitting are also tested against the available information in other bands: the reconstructed AGN bolometric luminosity is compared to those derived from X-rays and from the high excitation IR lines tracing AGN activity like [Ne V] and [O IV]. The IR luminosity due to star-formation (SF) and the intrinsic AGN bolometric luminosity are shown to be strongly related to the IR line luminosity. Variations of these relations with different AGN fractions are investigated, showing that the relation dispersions are mainly due to different AGN relative contribution within the galaxy. Extrapolating these local relations between line and SF or AGN luminosities to higher redshifts, by means of recent Herschel galaxy evolution results, we then obtain mid- and far-IR line luminosity functions useful to estimate how many star-forming galaxies and AGN we expect to detect in the different lines at different redshifts and luminosities with future IR facilities (e.g., JWST, SPICA).

Revisiting the relationship between 6 μm and 2–10 keV continuum luminosities of AGN

We have determined the relation between the AGN luminosities at rest-frame 6 {\mu}m associated to the dusty torus emission and at 2-10 keV energies using a complete, X-ray flux limited sample of 232 AGN drawn from the Bright Ultra-hard XMM-Newton Survey. The objects have intrinsic X-ray luminosities between 10^42 and 10^46 erg/s and redshifts from 0.05 to 2.8. The rest-frame 6 {\mu}m luminosities were computed using data from the Wide-Field Infrared Survey Explorer and are based on a spectral energy distribution decomposition into AGN and galaxy emission. The best-fit relationship for the full sample is consistent with being linear, L_6 {\mu}m $\propto$ L_2-10 keV^0.99$\pm$0.032, with intrinsic scatter, ~0.35 dex in log L_6 {\mu}m. The L_6 {\mu}m/L_2-10 keV luminosity ratio is largely independent on the line-of-sight X-ray absorption. Assuming a constant X-ray bolometric correction, the fraction of AGN bolometric luminosity reprocessed in the mid-IR decreases weakly, if at all, with the AGN luminosity, a finding at odds with simple receding torus models. Type 2 AGN have redder mid-IR continua at rest-frame wavelengths <12 {\mu}m and are overall ~1.3-2 times fainter at 6 {\mu}m than type 1 AGN at a given X-ray luminosity. Regardless of whether type 1 and type 2 AGN have the same or different nuclear dusty toroidal structures, our results imply that the AGN emission at rest-frame 6 {\mu}m is not isotropic due to self-absorption in the dusty torus, as predicted by AGN torus models. Thus, AGN surveys at rest-frame 6 {\mu}m are subject to modest dust obscuration biases.

ALMA OBSERVATIONS OF A CANDIDATE MOLECULAR OUTFLOW IN AN OBSCURED QUASAR

We present Atacama Large Millimeter/Submillimeter Array (ALMA) CO (1-0) and CO (3-2) observations of SDSS J135646.10+102609.0, an obscured quasar and ultra-luminous infrared galaxy (ULIRG) with two merging nuclei and a known 20-kpc-scale ionized outflow. The total molecular gas mass is M_{mol} ~ 9^{+19}_{-6} x 10^8 Msun, mostly distributed in a compact rotating disk at the primary nucleus (M_{mol} ~ 3 x 10^8 Msun) and an extended tidal arm (M_{mol} ~ 5 x 10^8 Msun). The tidal arm is one of the most massive molecular tidal features known; we suggest that it is due to the lower chance of shock dissociation in this elliptical/disk galaxy merger. In the spatially resolved CO (3-2) data, we find a compact (r ~ 0.3 kpc) high velocity (v ~ 500 km/s) red-shifted feature in addition to the rotation at the N nucleus. We propose a molecular outflow as the most likely explanation for the high velocity gas. The outflowing mass of M_{mol} ~ 7 x 10^7 Msun and the short dynamical time of t_{dyn} ~ 0.6 Myr yield a very high outflow rate of \dot{M}_{mol} ~ 350 Msun/yr and can deplete the gas in a million years. We find a low star formation rate (< 16 Msun/yr from the molecular content and < 21 Msun/yr from the far-infrared spectral energy distribution decomposition) that is inadequate to supply the kinetic luminosity of the outflow (\dot{E} ~ 3 x 10^43 erg/s). Therefore, the active galactic nucleus, with a bolometric luminosity of 10^46 erg/s, likely powers the outflow. The momentum boost rate of the outflow (\dot{p}/(Lbol/c) ~ 3) is lower than typical molecular outflows associated with AGN, which may be related to its compactness. The molecular and ionized outflows are likely two distinct bursts induced by episodic AGN activity that varies on a time scale of 10^7 yr.

SED DECOMPOSITION OF INFRARED-LUMINOUS GALAXIES

We select infrared-luminous galaxies by cross-matching the SDSS spectroscopic sample of galaxies with the WISE all-sky survey catalog. Based on photometric data points covering from SDSS u-band to WISE , their spectral energy distributions (SEDs) are separated into AGN, elliptical, spiral, and irregular galaxy components. The derived luminosities of spiral galaxy and AGN are well correlated with and [OIII] line luminosities, respectively. Most galaxies are dominated by young stellar populations even for optical AGNs, but at least 10% of optical non-AGNs appear to harbor buried AGNs. The AGN contribution increases dramatically with the total luminosity. These results show that the SED decomposition is successful and is useful to understand the true nature of dusty galaxies.