Issue |
A&A
Volume 520, September-October 2010
|
|
---|---|---|
Article Number | L4 | |
Number of page(s) | 5 | |
Section | Letters | |
DOI | https://doi.org/10.1051/0004-6361/201015486 | |
Published online | 27 September 2010 |
LETTER TO THE EDITOR
Molecular lines as tracers of Compton-thick AGN?![[*]](/icons/foot_motif.png)
I. Georgantopoulos1,2 - E. Rovilos3 - A. Akylas2 - E. Xilouris2
1 - INAF-Osservatorio Astronomico di Bologna, via Ranzani 1, 40127
Bologna, Italy
2 -
Institute of Astronomy and Astrophysics, National Observatory of
Athens, I. Metaxa & V. Pavlou str., Palaia Penteli, 15236 Athens,
Greece
3 -
Max-Planck-Institut für Extraterrestrische Physik,
Giessenbachstraße, 85748 Garching, Germany
Received 28 July 2010 / Accepted 2 September 2010
Abstract
Recent sub-mm CO molecular line
observations of nearby ultra-luminous IRAS galaxies, (U)LIRGs with
,
have shown that exceptionally
large gas column densities (
)
can
be present across some of the very dense gaseous disks that are typically
found in these objects. In particular, Papadopoulos et al. proposed
a diagnostic for finding such obscured sources using CO and HCN molecular lines.
Given that these high column densities are expected to absorb any
X-ray luminous active galactic nuclei (AGN), yielding Compton-thick sources, we set out to
explore whether this can be discerned using X-ray observations. More
specifically we examine X-ray spectral observations of 14 sources in
their sample, using public Chandra observations (0.5-10 keV) for eleven sources
as well as BeppoSAX results (2-100 keV) from the literature for another three sources.
Our goal is to find candidate Compton-thick AGN
and to check whether the molecular line selection criterion is
successful in selecting these systems.
X-ray spectroscopy reveals four candidate Compton-thick AGN, of which two fall within the
high obscuration region in the molecular line ratio diagnostics. Of
the remaining five sources falling into the ``high dust obscuration''
box, one (Mrk273) is highly obscured
(
), while in the other four the
X-ray emission is most probably associated with star-forming
processes rather than an AGN on the basis of their X-ray and
mid-infrared properties. Overall, we argue that although this method
cannot of course recover all Compton-thick AGN, there are no examples of
X-ray luminous AGN inside that region that have low obscuration, suggesting
that this method
is efficient in finding heavily obscured AGN in dust-enshrouded
star-forming galaxies. The above results bear important implications for
future joint ALMA and
X-ray observations for the detection of Compton-thick AGN.
Key words: X-rays: general - X-rays: diffuse background - X-rays: galaxies - galaxies: active - submillimeter: galaxies - infrared: galaxies
1 Introduction
X-ray surveys provide the most efficient tool for identifying active galactic nuclei (AGN) (Brandt & Hasinger 2005). This is because the hard X-rays can penetrate large amounts of dust and gas and thus can reveal hidden AGN, which would be very difficult to detect at optical wavelengths. However, even the hardest X-rays cannot easily escape from very large absorbing column densities, >1024



These extreme column densities could be relatively easily explored in the IR or
sub-mm part of the spectrum through their thermal or molecular emission.
For example, sub-mm observations have drawn much attention recently as a means
of discovering heavily obscured AGN (e.g. Alexander et al. 2005). The discovery
of significant dust optical depths even at submillimeter wavelengths
(in the case of Arp220; Sakamoto et al. 2008) reveals that such highly
obscured systems may not be uncommon. Recently, Papadopoulos et al. (2010)
(hereafter P10) made the interesting suggestion that observations of suitable
molecular lines can reveal cases of very high column densities (and thus
possible obscuration) present in some reservoirs of molecular gas in LIRGs. The
average molecular gas densities in some of these systems can be so high
that if the gas is distributed in a gaseous disk, it can yield extremely high
column densities and dust optical depths that suppress the high-J (6-5) CO line
emission to levels well below those expected in starburst galaxies such as
(ultra)-luminous IR galaxies, ULIRGs. In these cases the line of sight towards
the nucleus would be severely obscured making
this source Compton-thick, if an AGN were embedded in the
(U)LIRG. A molecular line diagnostic to identify these
severely obscured sources was proposed by P10. This helps to distinguish
between the degenerate possibilities of high dust obscuration and large
reservoirs of non star-forming cool and diffuse gas, whose high-J (6-5) CO
lines are genuinely faint rather than immersed to a strong dust continuum.
They suggest that high fractions of very dense molecular gas (with
)
can be traced by a high
HCN(1-0)/CO(1-0) ratio (R>0.1).
Then a low (R<0.3) CO(6-5)/(3-2) ratio
can only be the result of high dust optical depths (
).
Consequently there should be a very
high probability of Compton-thick AGN for (U)LIRG hosts lying inside their
aforementioned selection ``window''. On the other hand, we note that gaseous
disks (in the center of which the AGN are presumably embedded) with low average
column densities can often yield high column densities along the particular
line of sight towards the AGN as a result of their highly turbulent nature
(Wada et al. 2009). Thus Compton-thick sources can be expected also outside the
dust-obscuration selection window, but the vast majority of sources inside it
should be Compton-thick. In other words, if molecular gas disks have on
average high column densities, as the supression of CO J=6-5 emission from
their star-forming gas would imply, this will almost certainly be the case for
the special line of sight towards the AGN as well, but the reverse does not
necessarily hold. Thus the selection window of highly dust-obscured
star-forming systems advocated by P10 ought to select almost exclusively
Compton-thick AGN, but not all of them.
In this paper we are examining the X-ray properties of the systems in P10,
using both literature spectra as well as Chandra archival observations.
Our goal is to test whether AGN residing in the P10 box have a higher
probability for being classified as Compton-thick sources in X-rays.
We adopt
,
,
throughout the paper.
2 Data and analysis
2.1 The sample
Papadopoulos et al. (2010) present molecular line spectra for 20 (U)LIRGs. For six objects there are no HCN observations available so that the HCN/CO vs. CO(6-5)/CO(3-2) diagnostic cannot be applied. For the remaining 14 objects there are X-ray spectra available either in the literature or in publicly available Chandra observations. Although in many cases the Chandra observations may have been presented elsewhere, we chose to re-analyze the data in order to provide a uniform analysis of the sample. To our knowledge, in six cases X-ray spectral fits to the Chandra data are presented here for the first time: VII Zw 031, IRAS 10565+2448, IRAS 12112+0305, Arp 193, IRAS 17208-0014, IRAS 23365+3604. However, X-ray luminosities or hardness ratios derived from Chandra observations have been presented before in the literature for some of those (see Table 1). In the cases of Mrk 231, NGC 6240 and NGC 7469, where there are high signal-to-noise spectra in the 2-100 keV band reported in the literature from BeppoSAX, we prefer to present those data instead, as high energy observations can be much more efficient in detecting Compton-thick sources.
2.2 Spectral fits
We fit the Chandra data in the 0.5-10 keV band using XSPEC v.12.6 (Arnaud 1996). Initially, we fit an absorbed power-law model to the spectra. For more complex spectra, and judging on the basis of the data residuals from the model fit, we additionally fit a thermal component (Raymond-Smith) to account for X-ray emission from star-forming processes, and/or an FeK

We use the following criteria to classify a source as Compton-thick: a) the detection of an
absorption turnover, b) a flat X-ray spectrum
with
or flatter, suggestive of a reflection-dominted spectrum,
or finally c) a high EW (>
)
FeK
line, indicating a strongly
absorbed continuum
(see discussion in Georgantopoulos et al. 2009; Akylas & Georgantopoulos 2009; Matt et al. 2004, for justification and explanation of these diagnostics).
On the basis of the above criteria, we consider four Compton-thick sources:
IRAS 09320+6134, Mrk 231, Arp 220, and NGC 6240. For
Mrk 231 and NGC 6240 the absorption turnovers have been directly
detected by BeppoSAX (Braito et al. 2004 and Vignati et al. 1999
respectively) and therefore these can be considered as bona-fide Compton-thick sources.
In the case of Arp 220, the spectrum
is flat (
), but no Fe line can be detected in the Chandra data. However, in
the XMM-Newton data (Iwasawa et al. 2005) a high EW (1.9 keV) FeK
line at 6.7 keV
has been clearly detected, rendering this source a highly probable Compton-thick source.
The case of IRAS 09320+6134 (UGC 5101) is the most ambiguous one;
we detect a flat spectrum (
), albeit with a large uncertainty, together with a high
EW (
)
FeK
line.
The combination of a flat spectrum with a high EW FeK
line
rather points towards a Compton-thick AGN.
The EW derived here is consistent with the EW obtained by Ptak et al. (2003).
Imanishi et al. (2003) analyze both the XMM-Newton and the Chandra data and find
a considerably lower EW (
).
González-Martin et al. (2009)
claim the detection of a mildly Compton-thick AGN with a column density of
.
Finally, for IRAS 23365+3604, the spectrum is flat, but with a large uncertainty. Moreover, there is no
evidence for a high-EW FeK
line, therefore the possibility that this
source is Compton-thick is weaker.
2.3 Galaxies or AGN ?
From Table 1 we see that a large number of sources (with the exception of Mrk 273, IRAS 05189-2524 and the four Compton-thick AGN) are unobscured with column densities

In order to further investigate the energy mechanism of the sources that lie
in the window proposed by P10, we use the mid-IR diagnostic diagram of
Laurent et al. (2000). This traces the excess of the 5.5 m continuum
emission with respect to the 15
m continuum and 6.2
m PAH emission;
essentially it is a tracer of hot dust from the AGN, which has a higher
5.5
m/15
m flux ratio than the cooler dust associated with
H II regions, and low 6.2
m PAH flux. The sources within
the P10 box, namely Mrk 273,
IRAS 09320+6134, IRAS 10565+2448, Arp 220, IRAS 17208-0014, IRAS 12112+0305,
IRAS 23365+3604, show a fractional AGN contribution of 77%, 80%, 5%, 60%,
14%, 21% and 22% respectively,
(Armus et al. 2007, and V. Charmandaris private communication).
For comparison, we note that the two Compton-thick AGN
outside the box (Mrk 231 and NGC 6240) show AGN contributions of 95% and 50%
respectively (Armus et al. 2007).
Table 2: AGN diagnostics.
In Table 2 we summarize the results of the AGN diagnostics
used to examine the nature of the sources in the P10 box. There is an agreement
between the mid-IR and the X-ray diagnostics which find that 3/7 sources host
an AGN, if we use a 50% as a threshold for the Laurent et al. (2000) diagnostic.
We note that one source (IRAS 12112+0305)
has been observed with the VLA (Nagar et al. 2003). This source, which is
optically classified as a LINER, appears to show a radio core at 15 GHz with
a radio power of
.
Therefore, the possibility that this source hosts a Compton-thick AGN
cannot be securely excluded.
![]() |
Figure 1:
X-ray (2-10 keV) vs. FIR (4-100 |
Open with DEXTER |
![]() |
Figure 2: Diagnostic diagram of Papadopoulos et al. (2010). High HCN(1-0)/CO(1-0) together with low CO(6-5)/(3-2) ratios (region inside the box) mark the region where high average dust absorption likely surpresses the high-J (in this case J=6-5) CO emission that is emitted by a heavily dust-enshrouded starburst with large supplies of dense gas (and hence the large HCN/CO J=1-0 ratio). Filled circles denote the excess X-ray emission sources (i.e. likely AGN) according Fig. 1. Open circles denote sources falling within the Ranalli et al. (2003) relation for normal galaxies (all four CO(6-5)/(3-2) upper limits fall again within the same relation). Crosses denote the four Compton-thick sources according to X-ray spectroscopy. |
Open with DEXTER |
3 Discussion and conclusions
In Fig. 2 we present the diagnostic diagram of P10. We see
that out of the four Compton-thick sources (marked with crosses), two
(IRAS 09320+6134 and Arp 220) fall within the ``high dust obscuration''
window. The other one (NGC 6240) lies very close to the box with a
marginally lower HCN(1-0)/CO(1-0) ratio. The fourth source (Mrk 231)
lies away from the box showing highly excited CO lines and a high
HCN/CO J=1-0 line
ratio. This presents no fundamental problem in the Compton-thick
selection using molecular lines, because Compton-thick AGN are also
expected to lie outside the ``window'' of the diagnostic diagram
of P10. Moreover, this source displays peculiar characteristics
compared to the other three Compton-thick sources: although it is a
bona-fide Compton-thick source, it is a broad absorption line (BAL)
Quasi-Stellar Object (QSO). BAL is a special class of QSOs that show
highly ionized gas
flowing away from the central source with speeds in some cases
exceeding 30 000
.
Most BAL QSOs present high obscuring
columns in X-rays, but their broad-line-regions are still viewed
relatively unobscured by dust, probably meaning that their dense
obscuring medium is dust-free (e.g. Véron-Cetty & Véron 2000). It is
then reasonable that in the case of Mrk 231 the high
CO(6-5)/CO(3-2)
ratio implies that the CO emission is not suppressed by dust, despite
the presence of a high X-ray obscuring column density.
In addition, we need to know the
nature of the remaining sources within the P10 ``high dust obscuration''
box. Apart from the two Compton-thick AGN there are five more
sources presenting characteristically low CO(6-5)/CO(3-2) ratios in
combination with high HCN/CO J=1-0 ones. These are Mrk 273,
IRAS 17208-0014, IRAS 10565+2448, IRAS 23365+3604, and IRAS 12112+0305.
Among them only Mrk 273 has a known X-ray luminous, highly
obscured (
), but not
Compton-thick AGN.
Thus the total mean column density of a gaseous disk around the AGN
that would be responsible for its X-ray obscuration would be
.
On the other hand mean
disk column densities as low as
are
compatible with the observed suppression of the CO 6-5/3-2 ratio in
this object, which are
6.5 times higher than the one deduced by the X-ray data. This discrepancy remains considerable,
although it is not very large given the totally independent means used to obtain
(molecular lines vs. X-rays).
However, we note that the substantial difficulty of measuring CO
J=6-5 fluxes from the Earth with a large telescope like the JCMT
will more often yield underestimates rather than overestimates (because
pointing errors become substantial with respect to the narrow beams
used). This could tend to make molecular line-deduced and
X-ray-deduced hydrogen column densities more compatible in cases
like Mrk 273.
The X-ray emission of the remaining four sources in the box is
probably dominated by star-forming processes according to their X-ray
spectra, the
relation and the MIR diagnostics of
Laurent et al. (2000). This could explain why the sources do not appear as
Compton-thick in X-ray wavelengths, because the origin of the X-rays is
not a potentially deeply buried AGN but stellar. Thus these sources
are not useful in deducing mean column densities of gaseous disks
around AGN (and could be indeed found inside the selection
``window'' of P10 without invalidating its utility).
As a consistency check it is interesting to discuss the Si
properties of the objects within the P10 box, because it is
clearly impossible that these have substantial dust optical depths in short
sub-mm wavelengths as suggested by P10, while at the same time
they exhibit modest or low
Si (9.7
)
optical depths (e.g. Shi et al. 2006). Therefore
the highly obscured objects according to P10 should present
substantial Si optical depths. There are Spitzer-IRS spectroscopic
observations for four sources in the P10 box: Mrk 273, IRAS 09320+6134,
Arp 220 and IRAS 12112+0305 (Armus et al. 2007). Their Si (
)
optical depths are high: 1.8
0.4, 1.6
0.3, 3.3
0.2 and
1.3
0.3 respectively. This clearly suggests an association of all
the sources in the P10 box with optically-thick systems at
.
In conclusion, our results show that the diagnostic proposed by P10 offers a potentially useful tool for the detection of highly obscured X-ray luminous AGN. This bears important implications for future facilities such as ALMA. However, our work identified two caveats. First, the AGN inside the P10 box are heavily obscured but not necessarily Compton-thick, as the case of Mrk273 implies. Second, there can be large contaminations by powerful star-forming systems. Additional X-ray spectral observations might be necessary in order to sift those out. More molecular line data, ranging from low to high-J CO lines, as well as low-frequency (and thus unaffected by dust exinction) density indicators, such as low-to-modest-J HCN and CS lines of nearby (U)LIRGs are needed (along with X-ray data) to verify(or not) this picture in local (U)LIRGs that are known to have X-ray luminous AGN. The now superbly operating Herschel Space Observatory which can detect numerous high-J CO lines in local (U)LIRGs (e.g. van der Werf et al. 2010) in combination with ground based low-J CO and HCN, CS observations presents an excellent opportunity to do so. The large (U)LIRG sample that is now observed as part of the Herschel Key Project HerCULES, in combination with the extensive ground-based low-frequency CO and HCN survey will be ideal for such studies.
AcknowledgementsWe thank the anonymous referee for many useful comments and suggestions, which helped to substantially improve this paper. I.G. acknowledges the receipt of a Marie Curie fellowship FP7-PEOPLE-IEF-2008 Prop. 235285. We thank P. P. Papadopoulos for reading the manuscript and providing useful comments and suggestions. The Chandra data were taken from the Chandra Data Archive at the Chandra X-ray Center.
References
- Akylas, A., & Georgantopoulos, I. 2009, A&A, 500, 999 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Alexander, D. M., Smail, I., Bauer, F. E., et al. 2005, Nature, 434, 738 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Armus, L., Charmandaris, V., Bernard-Salas, J., et al. 2007, ApJ, 656, 148 [NASA ADS] [CrossRef] [Google Scholar]
- Arnaud, K. A. 1996, Astronomical Data Analysis Software and Systems V, ed. G. Jacoby, & J. Barnes, ASP Conf. Ser., 101, 17 [Google Scholar]
- Braito, V., Della Ceca, R., Piconcelli, E., et al. 2004, A&A, 420, 79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Brandt, W. N., & Hasinger, G. 2005, ARA&A, 43, 827 [Google Scholar]
- Cash, W. 1979, ApJ, 228, 939 [NASA ADS] [CrossRef] [Google Scholar]
- Clements, D. L., McDowell, J. C., Shaked, S., et al. 2002, ApJ, 581, 974 [NASA ADS] [CrossRef] [Google Scholar]
- Dadina, M. 2007, A&A, 461, 1209 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Georgantopoulos, I., Akylas, A., Georgakakis, A., & Rowan-Robinson, M. 2009, A&A, 507, 747 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- González-Martin, O., Masegosa, J., Marquez, M., Guainazzi, M., & Jimnez-Bailn, E. 2009, A&A, 506, 1107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Imanishi, M., Terashima, Y., Anabuki, N., & Nakagawa, T. 2003, ApJ, 596, L167 [NASA ADS] [CrossRef] [Google Scholar]
- Iwasawa, K., Sanders, D. B., Evans, A. S., et al. 2005, MNRAS, 357, 565 [NASA ADS] [CrossRef] [Google Scholar]
- Iwasawa, K., Sanders, D. B., Evans, A. S., et al. 2009, ApJ, 695, L103 [NASA ADS] [CrossRef] [Google Scholar]
- Laurent, O., Mirabel, I. F., Charmandaris, V., et al. 2000, A&A, 359, 887 [NASA ADS] [Google Scholar]
- Matt, G., Bianchi, S., Guainazzi, M., & Molendi, S. 2004, A&A, 421, 473 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Murphy, K. D., & Yaqoob, T. 2009, MNRAS, 397, 1549 [NASA ADS] [CrossRef] [Google Scholar]
- Nagar, N. M., Wilson, A. S., Falcke, H., Veilleux, S., & Maiolino, R. 2003, A&A, 409, 115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Papadopoulos, P. P., van der Werf, P., Isaak, K., & Xilouris, E. M. 2010, ApJ, 715, 775 [NASA ADS] [CrossRef] [Google Scholar]
- Ptak, A., Heckman, T., Levenson, N. A., Weaver, K., & Strickland, D. 2003, ApJ, 592, 782 [NASA ADS] [CrossRef] [Google Scholar]
- Ranalli, P., Comastri, A., & Setti, G. 2003, A&A, 399, 39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sakamoto, K., Wang, J., Wiedner, M. C., et al. 2008, ApJ, 684, 957 [NASA ADS] [CrossRef] [Google Scholar]
- Shi, Y., Rieke, G. H., Hines, D. C., et al. 2006, ApJ, 653, 127 [NASA ADS] [CrossRef] [Google Scholar]
- Teng, S. H., Wilson, A. S., Veilleux, S., et al. 2005, ApJ, 633, 664 [NASA ADS] [CrossRef] [Google Scholar]
- van der Werf, P. P., Isaak, K. G., Meijerink, R., et al. 2010, A&A, 518, L42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Veilleux, S., Sanders, D. B., & Kim, D. C. 1999, ApJ, 522, 139 [NASA ADS] [CrossRef] [Google Scholar]
- Véron-Cetty, M. P., & Véron, P. 2000, A&ARv, 10, 81 [NASA ADS] [CrossRef] [Google Scholar]
- Vignati, P., Molendi, S., Matt, G., et al. 1999, A&A, 349, L57 [NASA ADS] [Google Scholar]
- Wada, K., Papadopoulos, P. P., & Spaans, M. 2009, ApJ, 702, 63 [NASA ADS] [CrossRef] [Google Scholar]
- Xia, X. Y., Xue, S. J., Mao, S., et al. 2002, ApJ, 564, 196 [NASA ADS] [CrossRef] [Google Scholar]
Online Material
Table 1: X-ray properties.
Footnotes
- ... AGN?
- Table 1 is only available at http://www.aanda.org
All Tables
Table 2: AGN diagnostics.
Table 1: X-ray properties.
All Figures
![]() |
Figure 1:
X-ray (2-10 keV) vs. FIR (4-100 |
Open with DEXTER | |
In the text |
![]() |
Figure 2: Diagnostic diagram of Papadopoulos et al. (2010). High HCN(1-0)/CO(1-0) together with low CO(6-5)/(3-2) ratios (region inside the box) mark the region where high average dust absorption likely surpresses the high-J (in this case J=6-5) CO emission that is emitted by a heavily dust-enshrouded starburst with large supplies of dense gas (and hence the large HCN/CO J=1-0 ratio). Filled circles denote the excess X-ray emission sources (i.e. likely AGN) according Fig. 1. Open circles denote sources falling within the Ranalli et al. (2003) relation for normal galaxies (all four CO(6-5)/(3-2) upper limits fall again within the same relation). Crosses denote the four Compton-thick sources according to X-ray spectroscopy. |
Open with DEXTER | |
In the text |
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.