An Investigation into the Effects of Luminosity on the Mid‐Infrared Spectral Energy Distributions of Radio‐quiet Quasars

We explore the mid-infrared (mid-IR) through ultraviolet (UV) spectral energy distributions (SEDs) of 119,652 luminous broad-lined quasars with 0.064<z<5.46 using mid-IR data from Spitzer and WISE, near-infrared data from Two Micron All Sky Survey and UKIDSS, optical data from Sloan Digital Sky Survey, and UV data from Galaxy Evolution Explorer. The mean SED requires a bolometric correction (relative to 2500A) of BC=2.75+-0.40 using the integrated light from 1um-2keV, and we further explore the range of bolometric corrections exhibited by individual objects. In addition, we investigate the dependence of the mean SED on various parameters, particularly the UV luminosity for quasars with 0.5<z<3 and the properties of the UV emission lines for quasars with z>1.6; the latter is a possible indicator of the strength of the accretion disk wind, which is expected to be SED dependent. Luminosity-dependent mean SEDs show that, relative to the high-luminosity SED, low-luminosity SEDs exhibit a harder (bluer) far-UV spectral slope, a redder optical continuum, and less hot dust. Mean SEDs constructed instead as a function of UV emission line properties reveal changes that are consistent with known Principal Component Analysis (PCA) trends. A potentially important contribution to the bolometric correction is the unseen extream-UV (EUV) continuum. Our work suggests that lower-luminosity quasars and/or quasars with disk-dominated broad emission lines may require an extra continuum component in the EUV that is not present (or much weaker) in high-luminosity quasars with strong accretion disk winds. As such, we consider four possible models and explore the resulting bolometric corrections. Understanding these various SED-dependent effects will be important for accurate determination of quasar accretion rates.

We present a new method to calculate bolometric luminosities for unobscured, type 1 quasars with multiband photometric data. Bolometric luminosity is a fundamental property for understanding quasars, and it is commonly estimated from monochromatic luminosities using bolometric corrections that often neglect quasar spectral energy distribution (SED) diversity. We take advantage of the fact that most quasars now have multiband observations from UV to mid-IR and construct SEDs for a well-defined sample of Sloan Digital Sky Survey (SDSS) quasars at 0.5 ≤ z ≤ 2. Based on this fiducial sample, we explore quasar SEDs, their diversity, and their relations with bolometric luminosities. We then use unsupervised neural network self-organizing maps (SOM) to describe the SED diversity and compute the bolometric luminosities with a fully trained SOM model. This method reduces systematical uncertainties compared to the traditional method. In addition, we update the multilinear regression relations between bolometric luminosity and monochromatic luminosities at rest frame 1450 Å, 3000 Å, and 5100 Å. Our method is applicable to large quasar samples with a wide range of luminosity and redshift. We have applied it to the SDSS Data Release 16 quasars. We have also made our code publicly available.

Aims. We aim to characterize the X-ray-to-optical/near-infrared(NIR) broad-band emission of luminous quasars (QSOs) in the first gigayear (Gyr) of cosmic evolution in order to decipher whether or not they exhibit differences compared to the lower-z QSO population. Our goal is also to provide a reliable and uniform catalog of derivable properties for these objects (from fitting their spectral energy distribution), such as bolometric and monochromatic luminosities, Eddington ratios, dust extinction, and the strength of the hot dust emission. Methods. We gathered all available photometry –from XMM-Newton proprietary data in X-rays to rest-frame NIR wavelengths– for the 18 QSOs in the HYPERION samples (6.0 ≤ z ≤ 7.5). For sources lacking uniform NIR coverage, we conducted NIR observations in the J, H, and K bands. To increase the statistical robustness of our analysis across the UV-to-NIR region, we add 36 additional sources to our sample from the E-XQR-30 sample with 5.7 ≲ z ≲ 6.6. We characterized the X-ray/UV emission of each QSO using average SEDs from luminous Type 1 sources and calculated bolometric and monochromatic luminosities. Finally, we constructed a mean SED extending from the X-rays to the NIR bands. Results. We find that the UV-optical emission of these QSOs can be modeled with templates of z ∼ 2 luminous QSOs. We observe that the bolometric luminosities derived while adopting some bolometric corrections at 3000 Å (BC3000 Å) largely used in the literature are slightly overestimated, by 0.13 dex, as they also include reprocessed IR emission. We estimate a revised value of BC3000 Å = 3.3, which can be used to derive Lbol in z ≥ 6 QSOs. We provide a subsample of 11 QSOs with rest-frame NIR photometry; these show a broad range of hot dust emission strength, with two sources exhibiting low levels of emission. Despite potential observational biases arising from nonuniform photometric coverage and selection biases, we produce an X-ray-to-NIR mean SED for QSOs at z ≳ 6 that is a good match to templates of lower-redshift, luminous QSOs up to the UV–optical range, with a slightly enhanced contribution from hot dust in the NIR.

We measure the bolometric luminosity of a complete and unbiased 12 μm-selected sample of active galactic nuclei (AGN) in the local Universe. For each galaxy, we used a 10-band radio-to-X-ray spectral energy distribution (SED) to isolate the genuine AGN continuum in each band, including subarcsecond measurements where available, and correcting those contaminated by the host galaxy. We derive the median SED of Seyfert type 1 AGN, Seyferts with hidden broad lines (HBLs), Seyferts of type 2, and LINER nuclei in our sample. The median Seyfert 1 SED shows the characteristic blue bump feature in the UV, but nevertheless, the largest contribution to the bolometric luminosity comes from the IR and X-ray continua. The median SEDs of both HBL and type 2 AGN are affected by starlight contamination in the optical/UV. The median SED of HBL AGN is consistent with that of Seyfert 1s, when an extinction of A V ∼ 1.2 mag is applied. The comprehensive SEDs allowed us to measure accurate bolometric luminosities and derive robust bolometric corrections for the different tracers. The 12 μm and K-band nuclear luminosities have good linear correlations with the bolometric luminosity, similar to those in the X-rays. We derive bolometric corrections for either continuum bands (K band, 12 μm, 2–10 keV, and 14–195 keV) or narrow emission lines (mid-IR high-ionization lines of [O iv] and [Ne v] and optical [O iii] 5007 Å) as well as for combinations of IR continuum and line emission. A combination of continuum plus line emission accurately predicts the bolometric luminosity up to quasar luminosities (∼1046 erg s−1).

Context. Hyperluminous quasi-stellar objects (QSOs) are ideal laboratories to investigate active galactic nucleus (AGN) feedback mechanisms. Their formidable energy release causes powerful winds at all scales, and thus the maximum feedback is expected. Aims. Our aim is to derive the mean spectral energy distribution (SED) of a sample of 85 WISE-SDSS selected hyperluminous (WISSH) quasars. Since the SED provides a direct way to investigate the AGN structure, our goal is to understand if quasars at the bright end of the luminosity function have peculiar properties compared to the bulk of the QSO population. Methods. We collected all the available photometry, from X-rays to the far-infrared (FIR); each WISSH quasar is observed in at least 12 different bands. We then built a mean intrinsic SED after correcting for the dust extinction, absorption and emission lines, and intergalactic medium absorption. We also derived bolometric, IR band, and monochromatic luminosities together with bolometric corrections at λ = 5100 Å and 3 μm. We define a new relation for the 3 μm bolometric correction. Results. We find that the mean SED of hyperluminous WISSH QSOs shows some differences compared to that of less luminous sources (i.e., a lower X-ray emission and a near- and mid-IR excess which can be explained assuming a larger dust contribution. WISSH QSOs have stronger emission from both warm (T ∼ 500 − 600 K) and very hot (T ≥ 1000 K) dust, the latter being responsible for shifting the typical dip of the AGN SED from 1.3 μm to 1.1 μm. We also derived the mean SEDs of two subsamples created based on their spectral features (presence of broad absorption lines and equivalent width of CIV line). We confirm that broad absorption lines (BALs) are X-ray weak and that they have a reddened UV-optical continuum. We also find that BALs tend to have stronger emission from the hot dust component. For sources with a weaker CIV line, our main result is the confirmation of their lower X-ray emission. By populating the LIR vs. z diagram proposed by Symeonidis &amp; Page (MNRAS, 503, 3992), we found that ∼90% of WISSH QSOs with z ≥ 3.5 have their FIR emission dominated by star-forming activity. Conclusions. This analysis suggests that hyperluminous QSOs have a peculiar SED compared to less luminous objects. It is therefore critical to use SED templates constructed exclusively from very bright quasar samples (such as this one) when dealing with particularly luminous sources, such as high-redshift QSOs.

We present infrared bolometric luminosity corrections derived from the detailed spectral energy distributions of 62 bright quasars of low- to moderate-redshift (z=0.03-1.4). At 1.5, 2, 3, 7, 12, 15, and 24 microns we provide bolometric corrections of the mathematical forms L_iso=\zeta \lambda L_\lambda and log(L_iso)=A+B log(\lambda L_\lambda). Bolometric corrections for radio-loud and radio-quiet objects are consistent within 95% confidence intervals, so we do not separate them. Bolometric luminosities estimated using these corrections are typically smaller than those derived from some commonly used in the literature. We investigate the possibility of a luminosity dependent bolometric correction and find that, while the data are consistent with such a correction, the dispersion is too large and the luminosity range too small to warrant such a detailed interpretation. Bolometric corrections at 1.5 $\mu$m are appropriate for objects with properties that fall in the range log(L_bol)=45.4-47.3 and bolometric corrections at all other wavelengths are appropriate for objects with properties that fall in the range log(L_bol)=45.1-47.0.

Context. Astronomers have yet to establish whether high-mass protostars form from high-mass prestellar cores, similar to their lower-mass counterparts, or from lower-mass fragments at the heart of a pre-protostellar cluster undergoing large-scale collapse. Part of the uncertainty is due to a shortage of envelope structure data on protostars of a few tens of solar masses, where we expect to see a transition from intermediate-mass star formation to the high-mass process. Aims. We sought to derive the masses, luminosities, and envelope density profiles for eight sources in Cygnus-X, whose mass estimates in the literature placed them in the sampling gap. Combining these sources with similarly evolved sources in the literature enabled us to perform a meta-analysis of protostellar envelope parameters over six decades in source luminosity. Methods. We performed spectral energy distribution fitting on archival broadband photometric continuum data from 1.2 to 850 μm to derive bolometric luminosities for our eight sources plus initial mass and radius estimates for modelling density and temperature profiles with the radiative-transfer package Transphere. Results. The envelope masses, densities at 1000 AU, outer envelope radii, and density power law indices as functions of bolometric luminosity all follow established trends in the literature spanning six decades in luminosity. Most of our sources occupy an intermediate to moderately high range of masses and luminosities, which helps to more firmly establish the continuity between low- and high-mass star formation mechanisms. Our density power law indices are consistent with observed values in the literature, which show no discernible trends with luminosity, and have a mean p = −1.4 ± 0.4. However, our sub-sample, with a mean power law index of −1.1 ± 0.3, is slightly flatter than would be expected for spherical envelopes in free fall (p = −1.5). Conclusions. We attribute flattened density profiles for our eight sources to one or more of the following: ongoing accretion from their natal filaments, convolution of sources with neighbours or the larger filament, spherical averaging of asymmetric features (for example fragments), or inflation of the envelope by a moderate far-ultraviolet field. Finally, we show that the trends in all of the envelope parameters for high-mass protostars are statistically indistinguishable from trends in the same variables for low-mass protostars.

We present one of the largest multiwavelength studies of simultaneous optical-to-X-ray spectral energy distributions (SEDs) of unobscured (NH &lt; 1022 cm−2) active galactic nuclei (AGN) in the local Universe. Using a representative sample of hard-X-ray-selected AGN from the 70-month Swift/BAT catalog, with optical/UV photometric data from Swift/UVOT and X-ray spectral data from Swift/XRT, we constructed broadband SEDs of 236 nearby AGN (0.001 &lt; z &lt; 0.3). We employed GALFIT to estimate host galaxy contamination in the optical/UV and determine the intrinsic AGN fluxes. We used an absorbed power law with a reflection component to model the X-ray spectra and a dust-reddened multi-temperature blackbody to fit the optical/UV SED. We calculated intrinsic luminosities at multiple wavelengths, total bolometric luminosities (Lbol), optical-to-X-ray spectral indices (αox), and multiple bolometric corrections (κλ) in the optical, UV, and X-rays. We used black hole masses obtained by reverberation mapping and the virial method to estimate Eddington ratios (λEdd) for all our AGN. We confirm the tight correlation (scatter = 0.45 dex) between UV (2500 Å) and X-ray (2 keV) luminosity for our sample. We observe a significant decrease in αox with Lbol and λEdd, suggesting that brighter sources emit more UV photons per X-rays. We report a second-order regression relation (scatter = 0.15 dex) between the 2–10 keV bolometric correction (κ2 − 10) and αox, which is useful to compute Lbol in the absence of multiband SEDs. We also investigate the dependence of optical/UV bolometric corrections on the physical properties of AGN and obtain a significant increase in the UV bolometric corrections (κW2 and κM2) with Lbol and λEdd, unlike those in the optical (κV and κB), which are constant across five orders of Lbol and λEdd. We obtain significant dispersions (∼0.1–1 dex) in all bolometric corrections, and hence recommend using appropriate relations with observed quantities while including the reported scatter, instead of their median values.

Using a realistic model for line emission from the broad emission line regions of quasars, we are able to reproduce the previously observed correlations of emission-line ratios with the shape of the spectral energy distribution (SED). In agreement with previous studies, we find that the primary driving force behind the Baldwin effect (Wλ ∝ Lβ, β < 0) is a global change in the SED with quasar luminosity, in that more luminous quasars must have characteristically softer ionizing continua. This is completely consistent with observations that show (1) a correlation between Luv, αox, and αuvx; (2) correlations of SED shape-sensitive line ratios with αox, αuvx, and Luv; and (3) correlations between line equivalent widths and αox, αuvx, and Luv. However, to explain the complete lack of a correlation in the Wλ(N V)-Luv diagram, we propose that the more luminous quasars have characteristically larger gas metallicities (Z). As a secondary element, nitrogen's rapidly increasing abundance with increasing Z compensates for the losses in Wλ(N V) emitted by gas illuminated by softer continua in higher luminosity quasars. A characteristic relationship between Z and L has an impact on the Wλ-Luv relations for other lines as well. For a fixed SED, an increasing gas metallicity reduces the Wλ of the stronger metal lines (the gas cools), as well as that of Lyα and especially He II (because of the increasing metal opacity), while the weaker lines (e.g., C III] λ1909) generally respond positively. The interplay between the effects of a changing SED and Z with L results in the observed luminosity-dependent spectral variations. All of the resulting dependences on Luv are within the range of the observed slopes.

We employ contemporaneous optical, ultraviolet (UV) and X-ray observations from the XMM–Newton European Photon Imaging Camera (EPIC-pn) and Optical Monitor (OM) archives to present, for the first time, simultaneous spectral energy distributions (SEDs) for the majority of the Peterson et al. reverberation mapped sample of active galactic nuclei (AGN). The raw data were reduced using the latest pipelines and are all analysed consistently. The virial mass estimates from Peterson et al. allow us to calculate Eddington ratios λEdd for the sample using the bolometric accretion luminosities determined directly from the SEDs. We calculate hard X-ray bolometric corrections κ2–10 keV for the sample and confirm a trend for increasing bolometric correction with Eddington ratio proposed in previous studies. Our comparison with previous work on these objects suggests that the OM bandpass may be less susceptible to intrinsic reddening than the far-UV peak of the thermal disc spectrum in AGN, yielding larger bolometric corrections than previous work: κ2–10 keV≈ 15–30 for λEdd≲ 0.1, κ2–10 keV≈ 20–70 for 0.1 ≲λEdd≲ 0.2 and κ2–10 keV≈ 70–150 for λEdd≳ 0.2, but part of this increase could be attributed to spectral complexity preventing accurate recovery of the intrinsic luminosity in some sources. Long-term optical–UV variability contributes a second order, but significant change to the total bolometric luminosity when comparing multiple observations for individual objects. We also consider the effect of a recently proposed correction for radiation pressure when determining black hole masses with reverberation mapping, and find that the revised mass estimates do not significantly alter the range of bolometric corrections seen but may yield a narrower distribution of Eddington ratios.

Correlations of hot dust emission with outflow properties are investigated, based on a large z ∼ 2 non-broad absorption line quasar sample built from the Wide-field Infrared Survey and the Sloan Digital Sky Survey data releases. We use the near-infrared slope and the infrared to UV luminosity ratio to indicate the hot dust emission relative to the emission from the accretion disk. In our luminous quasars, these hot dust emission indicators are almost independent of the fundamental parameters, such as luminosity, Eddington ratio and black hole mass, but moderately dependent on the blueshift and asymmetry index (BAI) and FWHM of C iv lines. Interestingly, the latter two correlations dramatically strengthen with increasing Eddington ratio. We suggest that, in high Eddington ratio quasars, C iv regions are dominated by outflows so the BAI and FWHM (C iv) can reliably reflect the general properties and velocity of outflows, respectively. In low Eddington ratio quasars, on the other hand, C iv lines are primarily emitted by virialized gas so the BAI and FWHM (C iv) become less sensitive to outflows. Therefore, the correlations for the highest Eddington ratio quasars are more likely to represent the true dependence of hot dust emission on outflows and the correlations for the entire sample are significantly diluted by the low Eddington ratio quasars. Our results show that an outflow with a large BAI or velocity can double the hot dust emission on average. We suggest that outflows either contain hot dust in themselves or interact with the dusty interstellar medium or torus.

Thermal mid-infrared emission of quasars requires an obscuring structure that can be modeled as a magneto-hydrodynamic wind in which radiation pressure on dust shapes the outflow. We have taken the dusty wind models presented by Keating and collaborators that generated quasar mid-infrared spectral energy distributions (SEDs), and explored their properties (such as geometry, opening angle, and ionic column densities) as a function of Eddington ratio and X-ray weakness. In addition, we present new models with a range of magnetic field strengths and column densities of the dust-free shielding gas interior to the dusty wind. We find this family of models -- with input parameters tuned to accurately match the observed mid-IR power in quasar SEDs -- provides reasonable values of the Type 1 fraction of quasars and the column densities of warm absorber gas, though it does not explain a purely luminosity-dependent covering fraction for either. Furthermore, we provide predictions of the cumulative distribution of E(B-V) values of quasars from extinction by the wind and the shape of the wind as imaged in the mid-infrared. Within the framework of this model, we predict that the strength of the near-infrared bump from hot dust emission will be correlated primarily with L/L_Edd rather than luminosity alone, with scatter induced by the distribution of magnetic field strengths. The empirical successes and shortcomings of these models warrant further investigations into the composition and behaviour of dust and the nature of magnetic fields in the vicinity of actively accreting supermassive black holes.

We present observed mid-infrared and optical colors and composite spectral energy distributions (SEDs) of type 1 (broad-line) and 2 (narrow-line) quasars selected from Sloan Digital Sky Survey (SDSS) spectroscopy. A significant fraction of powerful quasars are obscured by dust and are difficult to detect in optical photometric or spectroscopic surveys. However, these may be more easily identified on the basis of mid-infrared (MIR) colors and SEDs. Using samples of SDSS type 1 and 2 matched in redshift and [O iii] luminosity, we produce composite rest-frame 0.2–15 μm SEDs based on SDSS, UKIDSS, and Wide-field Infrared Survey Explorer photometry and perform model fits using simple galaxy and quasar SED templates. The SEDs of type 1 and 2 quasars are remarkably similar, with the differences explained primarily by the extinction of the quasar component in the type 2 systems. For both types of quasar, the flux of the active galactic nucleus (AGN) relative to the host galaxy increases with AGN luminosity ( L [ O III ] ) and redder observed MIR color, but we find only weak dependencies of the composite SEDs on mechanical jet power as determined through radio luminosity. We conclude that luminous quasars can be effectively selected using simple MIR color criteria similar to those identified previously ( W 1 – W 2 &gt; 0.7 ; Vega), although these criteria miss many heavily obscured objects. Obscured quasars can be further identified based on optical–IR colors (for example, ( u – W 3 [ AB ] ) &gt; 1.4 ( W 1 – W 2 [ Vega ] ) + 3.2 ). These results illustrate the power of large statistical studies of obscured quasars selected on the basis of MIR and optical photometry.

We present Spitzer infrared spectra and ultra-violet to mid-infrared spectral energy distributions (SEDs) of 25 luminous type 1 quasars at z \sim 2. In general, the spectra show a bump peaking around 3 {\mu}m, and the 10 {\mu}m silicate emission feature. The 3 {\mu}m emission is identified with hot dust emission at its sublimation temperature. We explore two approaches to modeling the SED: (i) using the Clumpy model SED from Nenkova et al. (2008a), and (ii) the Clumpy model SED, and an additional blackbody component to represent the 3 {\mu}m emission. In the first case, a parameter search of \sim 1.25 million Clumpy models shows: (i) if we ignore the UV-to-near-IR SED, models fit the 2-8 {\mu}m region well, but not the 10 {\mu}m feature; (ii) if we include the UV-to-near-IR SED in the fit, models do not fit the 2-8 {\mu}m region. The observed 10 {\mu}m features are broader and shallower than those in the best-fit models in the first approach. In the second case, the shape of the 10 {\mu}m feature is better reproduced by the Clumpy models. The additional blackbody contribution in the 2-8 {\mu}m range allows Clumpy models dominated by cooler temperatures (T < 800K) to better fit the 8-12{\mu}m SED. A centrally concentrated distribution of a small number of torus clouds is required in the first case, while in the second case the clouds are more spread out radially. The temperature of the blackbody component is ~ 1200 K as expected for graphite grains.

Context. The bolometric luminosity is one of the most fundamental stellar parameters that is generally not directly observable. When integrating whole spectral energy distributions, bolometric corrections are widely used. These corrections are typically specified for a well defined object type, a measured photometric colour, and the colour index range for which it is calibrated. Aims. We provide this bolometric correction for the near-infrared colours of red late giants. Methods. We compare bolometric luminosities derived from fits to spectral energy distributions covering the visual and near, mid, and far infrared to different near-infrared colour indices to define subgroups within our stellar sample and finally obtain specific bolometric corrections. Results. For well-defined subgroups of different near-infrared colours and atmospheric chemistry, we present four distinct bolometric colour relations and compare these to earlier results from the literature.