X-ray spectroscopy of the quasar SDSS J1148+5251
We present the 78-ks Chandra observations of the quasar SDSS J1148+5251. The source is clearly detected in the energy range 0.3-7 keV with 42 counts (with a significance ). The X-ray spectrum is best-fitted by a power-law with photon index absorbed by a gas column density of . We measure an intrinsic luminosity at 2-10 keV and 10-40 keV equal to , comparable with luminous local and intermediate-redshift quasar properties. Moreover, the X-ray to optical power-law slope value () of J1148 is consistent with the one found in quasars with similar rest-frame 2500 Å luminosity (). Then we use Chandra data to test a physically motivated model that computes the intrinsic X-ray flux emitted by a quasar starting from the properties of the powering black hole and assuming that X-ray emission is attenuated by intervening, metal-rich () molecular clouds distributed on kpc scales in the host galaxy. Our analysis favors a black hole mass and a molecular hydrogen mass , in good agreement with estimates obtained from previous studies. We finally discuss strengths and limits of our analysis.
keywords:X-rays: galaxies; (galaxies:) quasars: super-massive black holes; galaxies: ISM; quasars: individual: SDSS J1148+5251
Studying the properties of quasars is important to understand the formation and evolution of super massive black holes (SMBHs) at the epoch of cosmic reionization and to investigate how they interact with the host galaxies, modifying their star formation history (e.g. Fan, 2012; Volonteri, 2012; Kormendy & Ho, 2013; Valiante et al., 2016). Hundreds of quasars have been discovered so far at (see recent numbers in Bañados et al., 2016), mostly through optical telescopes, and followed-up at different wavelengths, from the near-infrared to the radio regime (e.g. Gallerani et al., 2017). However, only tens of these high redshift quasars have been studied through their X-ray emission (Brandt et al., 2001; Mathur et al., 2002; Bechtold et al., 2003; Vignali et al., 2003b; Shemmer et al., 2006; Page et al., 2014; Moretti et al., 2014; Ai et al., 2016).
These studies were focused on X-ray luminous objects () and have found typical X-ray to optical power-law slopes that vary in the range . The X-ray spectra of these sources have been modelled by a power law with absorbed by low column densities of gas that are generally difficult to be robustly constrained.
Local and intermediate-redshift quasars of similar X-ray luminosities are typically characterized by (Just et al., 2007) and (Piconcelli et al., 2005). Thus, quasars do not show evident differences from their lower redshift counterparts. The only high- quasar that shows peculiar X-ray properties ( and ) is 0100+2802 at (Ai et al., 2016), powered by the most massive black hole () known so far at this redshift. However, these results are based on exploratory Chandra observations that consists of only 14 counts. In general, at these epochs, the number of X-ray observations is still sparse and often based on a limited number of counts ( 10-30), with the exception of ULAS J1120+0641 at (Page et al., 2014; Moretti et al., 2014) and SDSS J1030+0524 at (Farrah et al., 2004) that have been detected with 100 counts (Nanni et al., submitted to A&A).
In this work, we present X-ray observations of J1148+5251 (hereafter J1148), a quasar at , obtained with the Chandra telescope. Since its discovery (Fan et al., 2003; White et al., 2003), J1148 has been observed at several different wavelengths (e.g. Gallerani et al., 2008; Juarez et al., 2009; Gallerani et al., 2010; Riechers et al., 2009; Maiolino et al., 2012; Gallerani et al., 2014; Cicone et al., 2015; Stefan et al., 2015) and certainly represents the best studied case of quasars. In particular, from near-infrared observations of the MgII and CIV emission lines a black hole mass (Willott et al., 2003; Barth et al., 2003) has been derived. Moreover, several CO observations have been obtained for this object in the last years (Riechers et al., 2009; Gallerani et al., 2014; Stefan et al., 2015) suggesting a molecular hydrogen mass M and a size R kpc (Bertoldi et al., 2003; Walter et al., 2003; Walter et al., 2004).
In this paper we analyze the X-ray properties of J1148 and discuss whether X-ray observations, can be used to further constrain the molecular hydrogen mass contained in the host galaxy of this high-redshift quasar. The paper is organized as follows: in Sec. 3.1, we present the Chandra observations, the procedure adopted for data reduction and calibration, and the results from the X-ray spectral analysis; in Sec. 3.2 we discuss the serendipitous, tentative detection of X-rays at the location of an optically faint quasar, previously discovered by Mahabal et al. (2005); in Sec. 4, we describe a simple, physically motivated model for the X-ray emission from quasars and we use it to interpret our data; finally in Sec. 6 we discuss and summarize the results of this study.
2 Data reduction
The 78 ks Chandra observation (ObsId 17127) was performed in September 2, 2015. The target was observed with ACIS-S in Very Faint mode. In order to obtain the most updated data reduction and calibration files we reprocessed the data using the script chandra_repro included in CIAO v. 4.8 with CALDB v. 4.7.0. We created a S3 exposure corrected 0.5-7 keV image of the field of view (see Fig. 1) and run on it the source detection algorithm wavdetect with scales 1.0 2.0 and 4.0 pixels and false probability threshold
3.1 Spectral analysis
For the spectral analysis we used XSPEC v. 12.8.2 (Arnaud, 1996). Given the low number of counts (42) we grouped the spectrum at 1 net-count per bin and used C-stat (Cash, 1979) with direct background subtraction (Wstat in XSPEC; Wachter et al., 1979). Furthermore, in order to check the validity of our results against possible systematics, we also performed a joint fit with the source (not background subtracted) and background spectrum (grouped at 1 count per bin) using C-stat. We performed all the modellings in the energy range 0.3-5 keV. Indeed the source is detected up to keV. As J1148 is at , we are sampling a rest-frame energy range keV. The Galactic column density of , derived from Kalberla et al. (2005), was applied to all the models. In the following, errors are quoted at level and upper limits at .
We initially fitted the data with a simple power-law model (see Fig. 2; left panel) obtaining a best-fit slope of . This slope is slightly flatter than (but still consistent with) the typical value at found for quasars (; Piconcelli et al., 2005). We notice that the model slightly overestimates the data in the softest X-ray portion of the spectrum. For this reason, we added an intrinsic photoelectric component (ZWABS in XSPEC) to investigate the possible presence of absorbing material along the line of sight to the quasar. This resulted in a rather uncertain and steeper power-law photon index, , and a loosely constrained column density value for which we obtained an upper limit of . In order to obtain better constraints on the column density, we fix the slope to , consistently with the value found for the simple power-law model; the value is also within the range typically found for optically selected high- quasars (e.g. Vignali et al., 2003b). In this case we obtain a column density , still consistent with no absorption at 1.3.
We checked for the possible presence of the Compton reflection hump (George & Fabian, 1991) by adding to the simple power-law model a Compton reflection component, empirically parametrized in XSPEC by the PEXRAV model (Magdziarz & Zdziarski, 1995). This model calculates the reflection of the power-law primary AGN flux on an infinite slab of neutral material with infinite optical depth. The reflection strength is described by the R parameter which is defined as R=, where defines the solid angle subtended by the reflector. We set the energy cut-off of the input power-law spectrum in PEXRAV to 200 keV, the abundance to Solar and the inclination angle to . By fitting the observed spectrum with this model we obtain a photon index that is too steep (): the reflection component dominates almost all the spectral range at the expenses of the power-law, which is relevant at the softest energies. This result is unphysical and is due to the limited statistics. Therefore we fix again the power law slope to the canonical value, obtaining a . This value is in agreement with the low level of reflection observed in higher quality spectra of lower redshift luminous broad-line AGN (Vignali et al., 1999; Reeves & Turner, 2000; Page et al., 2005, Zappacosta et al. in prep.).
We checked our results also by modelling the total source+background extracted spectrum (i.e. without any background subtraction) jointly with the background spectrum using C-stat as fit statistics similarly to what has been done in ACIS-I, for instance, by Lanzuisi
et al. (2013). In order to do this we modelled the ACIS-S3 unfocused instrumental background spectrum (grouped to 1 count per bin) as a three-segment broken power law with slopes and the two energy break free to vary and two instrumental lines at 1.77 keV and 2.1 keV. The focused cosmic background was modelled with two thermal components with temperatures 0.07 keV and 0.2 keV with Solar metallicity (parametrizing the diffuse Galactic foregrounds; Lumb et al., 2002) and a power-law for the unresolved Cosmic X-ray Background with and normalization at 1 keV set to (De Luca &
By assuming the absorbed power-law model as fiducial best-fit parametrization, we estimated fluxes in the soft (0.5-2 keV) and hard (2-10 keV) X-ray band of respectively and . The intrinsic luminosities in the 2-10 and 10-40 keV energy bands are respectively of and . These luminosities are comparable to those typically measured in local and intermediate redshift quasars (i.e. Piconcelli et al., 2005, Zappacosta et al. in prep.).
Finally, we measure the X-ray to optical power-law slope defined as: , where and are the monochromatic flux densities at rest-frame 2 keV and 2500 Å, respectively. In particular, we use the composite quasar spectrum presented in Vanden Berk et al. (2001), to convert the broad-band z-filter measurements in , as in Vignali et al. (2003a). We find . Just et al. (2007) measured the in 32 luminous quasars at finding that anti-correlates with the their luminosity at Å (). By adopting such relation (see eq. 3 in Just et al., 2007), and considering that J1148 has , we infer , perfectly consistent with our estimate. This result implies that the of this luminous high- quasar is comparable with the ones found in lower redshift quasars.
3.2 Tentative X-ray emission from a quasar
The observed field of view encloses RD J11485253, a quasar discovered by Mahabal et al. (2005) at through optical observations in the R, z, and J bands. At the position of RD J11485253, four photons (three of these being contiguous) are detected in the observed-frame 0.5–7 keV (full) band in a circular region of 1.5″ radius. To understand whether this might be considered a detection and derive the corresponding significance level, we adopted three methods.
At first, similarly to Vignali et al. (2001), we extracted a 400400 pixel region centered on the optical position of the quasar, excluding the immediate vicinity of the quasar itself and masking few other detected sources. This region was covered with 30,000 circles of 1.5″ radius whose centers were randomly chosen through a Monte Carlo procedure. At this stage, to be conservative in the detection significance, we considered only the three contiguous counts (two of which in the hard, 2–7 keV, band) and found 161 and 518 extractions with at least three and two counts (out of the 30000 trials) in the full and hard band, respectively, corresponding to about 2.9 and 2.5 confidence level for a one-tailed Gaussian distribution in the two bands.
To verify this tentative detection, we considered the mean number of counts extracted in the full and hard band in the 30000 trials and assumed them as indicative of the background levels in the two bands, obtaining a probability of 2.6 and 2.2, respectively. Finally, we ran the source-detection tool wavdetect Freeman et al. (2002) on the full-band image to verify the results reported above. The source is detected at a false-positive threshold of 10. Although wavdetect is not properly a tool for photometry, we used its results in terms of number of counts. The source is detected with four counts, 0.7 of these due to the background. The Poisson probability of obtaining four total (sourcebackground) counts when 0.7 are expected in the extraction region is , corresponding to . From the source net counts and assuming Poissonian uncertainties (Gehrels, 1986) we estimate the 0.5-7 keV source count-rate to be . Assuming a power-law with canonical this corresponds to a 0.5-2 keV flux of . We then calculate the probability that the optical source and X-ray counterparts are associated by chance. Given the 0.5-2 keV fluxes we are probing and taking as reference the source counts by Lehmer et al. (2012), we expect to have in the ACIS-S3 field between 23 and 57 sources (accounting for uncertainty in flux). By assuming 1.5 radius for the X-ray counterpart we expect a random association probability in the range .
Summarizing, without assuming any prior (e.g., the optical position of the source, as in the Monte Carlo methods), the source is tentatively detected with a significance level of 2.6–2.9 in the full band and 2.2–2.5 in the hard band. The count rate in the 0.5-7 keV band () corresponds to an intrinsic luminosity of , assuming (). Though errors are large, this luminosity is typical of sources that are intermediate between bright Seyfert and faint quasars. Moreover, given the observed J-band magnitude () of RD J11485253, and extrapolating the flux to Å with the Vanden Berk et al. (2001) template (), we infer .
4 Modelling X-ray emission in J1148
In Sec. 3.1, we have seen that if we fix the spectral slope to the value that is commonly found at lower redshift () we can marginally constrain the column density () of the absorbing gas that is intervening along the line of sight towards us. Although our detection of the X-ray absorption is at a low significance level, here we speculate on its origin and assume that it is due to metals locked in molecular clouds, distributed on kpc scales (see the model by Gallerani et al., 2014, hereafter G14). In Sec. 6 we extensively discuss this assumption.
We compute the observed X-ray flux as , where is the intrinsic X-ray spectrum of J1148 and is the optical depth of neutral hydrogen encountered by X-ray photons in the host galaxy along the line of sight towards us. We model with a single power-law having slope . Given the bolometric luminosity , where is the Eddington ratio, and is the Eddington luminosity, the X-ray luminosity
where is the cloud column density given by , and is the number of clouds along a line of sight:
where is total number of clouds of mass , is the volume occupied by , and is the geometrical cross section of each cloud (). Thus, can be rewritten as:
To summarize, under the aforementioned assumptions, is completely specified by , , and for the intrinsic spectrum, and , and for the absorption model. A comprehensive analysis of the full parameter space is hampered by the limited number of counts of the Chandra spectrum combined with the large number of free parameters of our model. Moreover, some of the free parameters are degenerate, as with and . Nevertheless, J1148 represents the best studied case of quasars, thus providing a perfect laboratory to test the hypothesis that the X-ray photons absorbing material may be constituted by molecular clouds distributed on kpc scales. In what follows, we consider and as free parameters of the model, and we assume for the others reasonable values taken from the literature.
We assume that J1148 shines at the Eddington luminosity, i.e. . This assumption is supported by several works, e.g. Willott et al. (2003) and Schneider et al. (2015) for the specific case of J1148 and Wu et al. (2015) for the general case of high-z quasars (see their Fig. 4). For the intrinsic spectrum, we adopt (as assumed in our analysis in Sec. 3.1), and the bolometric correction by Lusso et al. (2012). We note that the value at the J1148 bolometric luminosity is extrapolated from results obtained at lower ; however, it is consistent within 1 with the work by Hopkins et al. (2007) (i.e. ), and with the results by Marconi et al. (2004) (i.e. ) for the same bolometric luminosity. More importantly, our choice is justified by the fact that if we take into account the black hole mass measured by Willott et al. (2003) and Barth et al. (2003), , the bolometric correction by Lusso et al. (2012) returns an X-ray luminosity () that is perfectly consistent with our measurements (see Sec. 3.1).
For what concerns the absorption model, high angular resolution (0.15”, namely 1 kpc at ) VLA observations of the CO(3-2) emission by Walter et al. (2004) have shown that the molecular gas is extended out to radii of . We adopt the metallicity value obtained from near-infrared observations of quasars (Juarez et al., 2009; Nagao et al., 2006), scaled for the solar metallicity (Asplund et al., 2009).
We constrain and by comparing our model with the X-ray data, shown in Fig. 3 as red filled circles. For this analysis, we re-bin the spectrum (see Fig. 2, left panel) in such a way that each bin contains at least 12-14 background-subtracted counts
The flux density obtained with the fiducial unabsorbed spectrum is shown in Fig. 3 with the black dashed line and the gray shaded region that reflects the uncertainties; the solid line and the cyan hatched region represent the best-fit model to the observed spectrum with the corresponding uncertainties. The values that minimize the are and , where the errors provide the uncertainties for a chi-square distribution with one degree of freedom. The mass range corresponds to a column density . Resulting values of and are perfectly consistent with previous black hole mass estimates by Willott et al. (2003) and Barth et al. (2003) and the mass inferred by Walter et al. (2003) from CO(3-2) observations.
Finally, we explore how the best-fit value of and would change if we vary the fiducial parameters (, , , ). We do not consider variations on since the results would be degenerate with . Moreover, for what concerns , we note that the measurements by Nagao et al. (2006) refers to the broad line regions (BLRs), small nuclear regions ( few pc) that may be characterized by larger metallicity values with respect to the one in the host galaxy (e.g. Valiante et al., 2011). Thus, we consider as lower limit for the metallicity of J1148. The results of this analysis are reported in Table 1. Both and shows the largest variations (up to a factor of 10) with .
These results represent an important consistency check of our model. In other words, we are saying that Chandra observations can be explained by assuming that the intrinsic X-ray spectrum of J1148 is attenuated by intervening, metal-rich () molecular clouds of total mass , distributed at kpc distance from the center of the galaxy, characterized by and density .
However, given the large uncertainties that are plaguing our analysis and the novelty of our approach some caution is needed. In the next section we thoroughly discuss both strengths and limits of this study.
Estimates of rely on the relation , where is the conversion factor, in units of , and is the CO(1-0) luminosity, in units of . For what concerns the conversion factor, in Galactic molecular clouds , while ultra-luminous infrared galaxies (ULIRGs, and nuclear star-burst galaxies have . However, this factor remains highly uncertain (see Bolatto
et al., 2013, for a review on this subject) and generally presents strong variations with the metallicity (e.g. Narayanan et al., 2012) and with the star formation rate (e.g. Clark &
Glover, 2015). Moreover, CO(1-0) data are generally poorly constrained at high- since the CO(1-0) frequency () falls outside the ALMA bands for . Thus, no quasar has been detected so far in this transition. In the J1148 case, for example, the estimate is based on the assumption of a constant brightness in the CO(3-2) and CO(1-0) transitions.
On the one hand, all these caveats make the possibility of measuring through alternative, independent methods extremely appealing. On the other hand, our result is valid under the assumption that the intervening gas absorbing X-ray photons is constituted by molecular clouds (MC) distributed on kpc scales, larger than those generally considered in more standard scenarios.
According to the unified-model (e.g. Antonucci, 1993; Urry & Padovani, 1995), the large variety of AGN types is just the result of varying orientation relative to the line of sight. In this scenario, objects that appear as optical TypeII AGN are observed along a line of sight covered by an axisymmetric dusty structure (the “torus”) with dimensions pc. TypeI AGNs are instead not obscured since the line of sight does not intercept the torus itself. Since its original formulation, the unified-model has been often revisited (see Netzer, 2015, for a recent review discussing the origin and properties of the central obscurer and the connection with its surrounding). For example, another possible origin for X-ray absorption concerns with gas that is located on broad line region (BLR, e.g. Peterson, 2006) scales ( pc) and constituted by high-density (), metal-rich () clouds (Risaliti et al., 2002). This scenario (hereafter called BLR-scenario) nicely accounts for the X-ray emission variability of AGN
The existence of all the aforementioned scenarios demonstrates that the origin of the X-ray absorption is far from being definitively understood and must occur on a wide range of scales. This is in line with the complex picture discussed by Elvis (2012) that takes into account the BLR, the torus obscuration, and the presence of outer ( kpc) absorbing regions connected to the host disk and/or dust lanes (see also the discussion on the torus-galaxy connection by Netzer, 2015). In particular, high- quasars are expected to form in massive, over-dense regions, with dark matter halo mass and virial radii kpc, that result from of a long history of lower mass systems merging (e.g. Valiante et al., 2014). More in general, since high-redshift galaxies are thought to be very compact and gas-rich systems, it is plausible that their interstellar medium may provide substantial contribution to the X-ray emission obscuration (e.g. Bournaud et al., 2011; Juneau et al., 2013; Gilli et al., 2014). Thus, the assumption of having intervening material distributed also on kpc scales sounds plausible.
We further note that the mass inferred from our analysis corresponds to a column density of neutral hydrogen . Since we are assuming a metallicity close to (or even higher than) the solar one, if the Galactic value (e.g. Güver &
Özel, 2009) would apply to J1148, this column density implies an . This value is much larger that the one inferred in this object from near infrared observations (, Gallerani
et al., 2010). This discrepancy is commonly found in several optically bright quasars (e.g. Maccacaro
et al., 1982; Granato et al., 1997; Maiolino
et al., 2001a).
In the BLR scenario, high gas column densities can provide a small since BLRs are confined within the dust sublimation radius ( pc in the case of J1148, from eq. 1 by Netzer, 2015) and are consequently composed of a dust-poor gas. Moreover, anomalous properties of dust grains could also explain an value smaller than the Galactic one Maiolino et al. (2001b).
In fact, for a given total dust mass and fixed dust composition, a grain size distribution flatter than the Galactic one
Also in the MC scenario there are lines of reasoning that can make large gas column densities consistent with small value. Ionized, atomic and molecular outflows are commonly observed towards local and high- quasars (e.g. Feruglio et al., 2010; Rupke & Veilleux, 2011; Carniani et al., 2015, 2016; Feruglio et al., 2015, see also Fabian (2012) for a review on this subject). In particular, observations of the CO(1-0) emission line in local ultra-luminous galaxies have found that a non-negligible fraction (from few% up to 30%) of the total mass is located in outflows extended on kpc scales (Cicone et al., 2014). For what concerns J1148, while observations of CO emission lines are still not sensitive enough for detecting molecular outflows, PdBI data of the [CII] 158 line have provided the evidence for outflowing gas in this high- quasar
To conclude, our study has the important implication that X-ray observations can provide measurements that are independent from millimeter observations of the CO emission lines. Nevertheless, given the marginal detection of the X-ray absorption resulting from our analysis and the large uncertainty on the origin of X-ray absorbers, our interpretation of the results remains speculative and need to be tested with higher quality X-ray data. Deeper Chandra observations are necessary to better constrain the X-ray power-law photon index, and thus the column density of the absorbing gas. By detecting twice the number of X-ray photons collected in this work the error on would be reduced by 30%. Moreover, given the novelty of the approach we have used, the same analysis should be applied to other quasars for which both and have been measured through more standard techniques (e.g. MgII, CIV and CO emission lines).
We have presented the X-ray observation of the quasar J1148+5251 () obtained with 78 ks of Chandra observing time. We have clearly detected the source with 42 net counts (with a significance of ) in the observed 0.3–7 keV energy band. We have modelled the X-ray spectrum with a power-law (photon index ) absorbed by a gas column density of . We have found that the X-ray properties of J1148 are comparable to those of lower redshift quasars: (i) we have measured intrinsic luminosities of and , respectively; (ii) we have checked for the presence of a reflection component over the simple power-law modelling finding a low degree of reflection from circum-nuclear material; (iii) we have measured an X-ray to optical power-law slope ().
In the same field of J1148, at the position of RD J11485253, the quasar discovered by Mahabal et al. (2005), we have detected three (two) photons in the full (hard) X-ray band. We have adopted three different methods to estimate the reliability of these detections and we have found a significance level of 2.8 (2.3), i.e. the source is only marginally detected by Chandra. Deeper observations are required to confirm the detection of this source. The tentatively inferred X-ray luminosity of RD J11485253 () is typical of sources that are intermediate between the Seyfert and quasar luminosity regime. Moreover, given its observed J band magnitude (), we estimate .
Finally, we have used Chandra observations to test a physically motivated model that computes the intrinsic X-ray flux emitted by a quasar starting from the properties of the powering black hole and assumes that X-ray emission is attenuated by intervening, metal-rich () molecular clouds distributed on kpc scales in the host galaxy. Our analysis favors a black hole mass and a molecular hydrogen mass , in good agreement with estimates obtained from previous studies, thus providing a solid consistency check to our model.
This work highlights the importance of the synergy between X-ray and millimeter data for studying the properties of high- quasars and their host galaxies.
RM acknowledges support from the ERC Advanced Grant 695671 “QUENCH” and from the Science and Technology Facilities Council (STFC).
- pagerange: X-ray spectroscopy of the quasar SDSS J1148+5251–References
- pubyear: 2016
- This value translates in the whole S3 chip to 1 statistical fluctuation identified as point-source.
- This kind of instrumental+cosmic background modeling is usually employed in faint point-source and diffuse source studies (Fiore et al., 2012; Humphrey et al., 2006)
- The subscript in quantities as , , refers to the band 2-10 keV, if not differently stated.
- The analysis can provide biased results when falls below 10-20 counts per bin Cash (1979); Humphrey et al. (2009). Nevertheless, the large uncertainties on the adopted fiducial parameters cover a range of values larger than a possible systematic bias.
- Variability studies of local AGN () indicate that more luminous sources vary with a lower amplitude, implying that the X-ray variability of J1148 is expected to be modest (e.g. Lawrence & Papadakis, 1993; Almaini et al., 2000; Paolillo et al., 2004). However, there are hints that quasars of the same X-ray luminosity are more variable at (e.g. Vagnetti et al., 2016; Manners et al., 2002). Nevertheless, this redshift dependence is only tentative and can be explained by selection effects (e.g. Lanzuisi et al., 2014; Shemmer et al., 2014). Given the high luminosity of J1148 and the uncertain redshift dependence of AGN variability, we have ignored in our analysis possible X-ray variability effects.
- The dust grain sizes distribution can be modelled through a power-law: , where is the grain diameter that varies in the range . The Galactic curve in the diffuse ISM can be described by a mixture of graphite and silicate grains having sizes distributed with , m and m (Mathis et al., 1977).
- We note that the [CII] emission is extended on scales ( kpc) larger than the CO emission. The origin of this extended [CII] emission is not clear, as discussed in details in Cicone et al. (2014). [CII] emission can both arise from dense () photo-dissociation regions (PDRs) and more diffuse () neutral/ionized gas, though the dominant contribution is expected to arise from PDRs (e.g. Vallini et al., 2013, 2015). From Table 1, we note that if we would include contribution from PDRs, namely we would consider an , we would get a total mass that is a factor larger than the one inferred by CO observations. In the case more sensitive observations would detect CO emission on the [CII] scales, our model could be easily updated to acount for larger masses. Viceversa, if the extended [CII] emission arises from diffuse gas, its contribution to the X-ray absorbing gas would be negligible (). For this reason, we have not taken into account possible contribution from gas extended on the [CII] emission scales.
- Ai Y., Dou L., Fan X., Wang F., Wu X.-B., Bian F., 2016, ApJ, 823, L37
- Almaini O., et al., 2000, MNRAS, 315, 325
- Antonucci R., 1993, ARA&A, 31, 473
- Arnaud K. A., 1996, in Jacoby G. H., Barnes J., eds, Astronomical Society of the Pacific Conference Series Vol. 101, Astronomical Data Analysis Software and Systems V. p. 17
- Asplund M., Grevesse N., Sauval A. J., Scott P., 2009, ARA&A, 47, 481
- Bañados E., et al., 2016, preprint, (arXiv:1608.03279)
- Barth A. J., Martini P., Nelson C. H., Ho L. C., 2003, ApJ, 594, L95
- Bechtold J., et al., 2003, ApJ, 588, 119
- Bertoldi F., et al., 2003, A&A, 409, L47
- Bianchi S., Piconcelli E., Chiaberge M., Bailón E. J., Matt G., Fiore F., 2009, ApJ, 695, 781
- Bolatto A. D., Wolfire M., Leroy A. K., 2013, ARA&A, 51, 207
- Bournaud F., Dekel A., Teyssier R., Cacciato M., Daddi E., Juneau S., Shankar F., 2011, ApJ, 741, L33
- Brandt W. N., Guainazzi M., Kaspi S., Fan X., Schneider D. P., Strauss M. A., Clavel J., Gunn J. E., 2001, AJ, 121, 591
- Cappi M., et al., 2016, A&A, 592, A27
- Carniani S., et al., 2015, A&A, 580, A102
- Carniani S., et al., 2016, A&A, 591, A28
- Cash W., 1979, ApJ, 228, 939
- Cicone C., et al., 2014, A&A, 562, A21
- Cicone C., et al., 2015, A&A, 574, A14
- Clark P. C., Glover S. C. O., 2015, MNRAS, 452, 2057
- De Luca A., Molendi S., 2004, A&A, 419, 837
- Dwek E., Foster S. M., Vancura O., 1996, ApJ, 457, 244
- Elvis M., 2012, Journal of Physics Conference Series, 372, 012032
- Fabian A. C., 2012, ARA&A, 50, 455
- Fan X., 2012, Research in Astronomy and Astrophysics, 12, 865
- Fan X., et al., 2003, AJ, 125, 1649
- Farrah D., Priddey R., Wilman R., Haehnelt M., McMahon R., 2004, ApJ, 611, L13
- Ferrara A., Scannapieco E., 2016, preprint, (arXiv:1609.02561)
- Feruglio C., Maiolino R., Piconcelli E., Menci N., Aussel H., Lamastra A., Fiore F., 2010, A&A, 518, L155
- Feruglio C., et al., 2015, A&A, 583, A99
- Fiore F., et al., 2012, A&A, 537, A16
- Freeman P. E., Kashyap V., Rosner R., Lamb D. Q., 2002, ApJS, 138, 185
- Gallerani S., Ferrara A., Fan X., Choudhury T. R., 2008, MNRAS, 386, 359
- Gallerani S., et al., 2010, A&A, 523, A85
- Gallerani S., Ferrara A., Neri R., Maiolino R., 2014, MNRAS, 445, 2848
- Gallerani S., Fan X., Maiolino R., Pacucci F., 2017, preprint, (arXiv:1702.06123)
- Gehrels N., 1986, ApJ, 303, 336
- George I. M., Fabian A. C., 1991, MNRAS, 249, 352
- Gilli R., et al., 2014, A&A, 562, A67
- Granato G. L., Danese L., Franceschini A., 1997, ApJ, 486, 147
- Güver T., Özel F., 2009, MNRAS, 400, 2050
- Hopkins P. F., Richards G. T., Hernquist L., 2007, ApJ, 654, 731
- Humphrey P. J., Buote D. A., Gastaldello F., Zappacosta L., Bullock J. S., Brighenti F., Mathews W. G., 2006, ApJ, 646, 899
- Humphrey P. J., Liu W., Buote D. A., 2009, ApJ, 693, 822
- Juarez Y., Maiolino R., Mujica R., Pedani M., Marinoni S., Nagao T., Marconi A., Oliva E., 2009, A&A, 494, L25
- Juneau S., et al., 2013, ApJ, 764, 176
- Just D. W., Brandt W. N., Shemmer O., Steffen A. T., Schneider D. P., Chartas G., Garmire G. P., 2007, ApJ, 665, 1004
- Kalberla P. M. W., Burton W. B., Hartmann D., Arnal E. M., Bajaja E., Morras R., Pöppel W. G. L., 2005, A&A, 440, 775
- Kormendy J., Ho L. C., 2013, ARA&A, 51, 511
- Lanzuisi G., et al., 2013, MNRAS, 431, 978
- Lanzuisi G., et al., 2014, ApJ, 781, 105
- Lawrence A., Papadakis I., 1993, ApJ, 414, L85
- Lehmer B. D., et al., 2012, ApJ, 752, 46
- Lumb D. H., Warwick R. S., Page M., De Luca A., 2002, A&A, 389, 93
- Lusso E., et al., 2012, MNRAS, 425, 623
- Maccacaro T., Perola G. C., Elvis M., 1982, ApJ, 257, 47
- Magdziarz P., Zdziarski A. A., 1995, MNRAS, 273, 837
- Mahabal A., Stern D., Bogosavljević M., Djorgovski S. G., Thompson D., 2005, ApJ, 634, L9
- Maiolino R., Marconi A., Salvati M., Risaliti G., Severgnini P., Oliva E., La Franca F., Vanzi L., 2001a, A&A, 365, 28
- Maiolino R., Marconi A., Oliva E., 2001b, A&A, 365, 37
- Maiolino R., et al., 2012, MNRAS, 425, L66
- Manners J., Almaini O., Lawrence A., 2002, MNRAS, 330, 390
- Marconi A., Risaliti G., Gilli R., Hunt L. K., Maiolino R., Salvati M., 2004, MNRAS, 351, 169
- Mathis J. S., Rumpl W., Nordsieck K. H., 1977, ApJ, 217, 425
- Mathur S., Wilkes B. J., Ghosh H., 2002, ApJ, 570, L5
- Matt G., 2000, A&A, 355, L31
- Moretti A., et al., 2014, A&A, 563, A46
- Morrison R., McCammon D., 1983, ApJ, 270, 119
- Nagao T., Marconi A., Maiolino R., 2006, A&A, 447, 157
- Narayanan D., Krumholz M. R., Ostriker E. C., Hernquist L., 2012, MNRAS, 421, 3127
- Netzer H., 2015, ARA&A, 53, 365
- Page K. L., Reeves J. N., O’Brien P. T., Turner M. J. L., 2005, MNRAS, 364, 195
- Page M. J., Simpson C., Mortlock D. J., Warren S. J., Hewett P. C., Venemans B. P., McMahon R. G., 2014, MNRAS, 440, L91
- Paolillo M., Schreier E. J., Giacconi R., Koekemoer A. M., Grogin N. A., 2004, ApJ, 611, 93
- Peterson B. M., 2006, in Alloin D., ed., Lecture Notes in Physics, Berlin Springer Verlag Vol. 693, Physics of Active Galactic Nuclei at all Scales. p. 77, doi:10.1007/3-540-34621-X˙3
- Piconcelli E., Jimenez-Bailón E., Guainazzi M., Schartel N., Rodríguez-Pascual P. M., Santos-Lleó M., 2005, A&A, 432, 15
- Reeves J. N., Turner M. J. L., 2000, MNRAS, 316, 234
- Riechers D. A., et al., 2009, ApJ, 703, 1338
- Risaliti G., Elvis M., Nicastro F., 2002, ApJ, 571, 234
- Risaliti G., Elvis M., Fabbiano G., Baldi A., Zezas A., 2005, ApJ, 623, L93
- Rupke D. S. N., Veilleux S., 2011, ApJ, 729, L27
- Schneider R., Bianchi S., Valiante R., Risaliti G., Salvadori S., 2015, A&A, 579, A60
- Shemmer O., et al., 2006, ApJ, 644, 86
- Shemmer O., et al., 2014, ApJ, 783, 116
- Stefan I. I., et al., 2015, MNRAS, 451, 1713
- Urry C. M., Padovani P., 1995, PASP, 107, 803
- Vagnetti F., Middei R., Antonucci M., Paolillo M., Serafinelli R., 2016, A&A, 593, A55
- Valiante R., Schneider R., Salvadori S., Bianchi S., 2011, MNRAS, 416, 1916
- Valiante R., Schneider R., Salvadori S., Gallerani S., 2014, MNRAS, 444, 2442
- Valiante R., Schneider R., Volonteri M., 2016, Publ. Astron. Soc. Australia, 33, e032
- Vallini L., Gallerani S., Ferrara A., Baek S., 2013, MNRAS, 433, 1567
- Vallini L., Gallerani S., Ferrara A., Pallottini A., Yue B., 2015, ApJ, 813, 36
- Vanden Berk D. E., et al., 2001, AJ, 122, 549
- Vignali C., Comastri A., Cappi M., Palumbo G. G. C., Matsuoka M., Kubo H., 1999, ApJ, 516, 582
- Vignali C., Brandt W. N., Fan X., Gunn J. E., Kaspi S., Schneider D. P., Strauss M. A., 2001, AJ, 122, 2143
- Vignali C., Brandt W. N., Schneider D. P., 2003a, AJ, 125, 433
- Vignali C., et al., 2003b, AJ, 125, 2876
- Volonteri M., 2012, Science, 337, 544
- Wachter K., Leach R., Kellogg E., 1979, ApJ, 230, 274
- Walter F., et al., 2003, Nature, 424, 406
- Walter F., Carilli C., Bertoldi F., Menten K., Cox P., Lo K. Y., Fan X., Strauss M. A., 2004, ApJ, 615, L17
- White R. L., Becker R. H., Fan X., Strauss M. A., 2003, AJ, 126, 1
- Willott C. J., McLure R. J., Jarvis M. J., 2003, ApJ, 587, L15
- Wu X.-B., et al., 2015, Nature, 518, 512
- Yaqoob T., 1997, ApJ, 479, 184