EDP Sciences
Free Access
Issue
A&A
Volume 504, Number 1, September II 2009
Page(s) 53 - 59
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/200811516
Published online 09 July 2009

A Suzaku observation of the ULIRG IRAS19254-7245: discerning the AGN component

V. Braito1,2 - J. N. Reeves2,3 - R. Della Ceca4 - A. Ptak2,5 - G. Risaliti6,7 - T. Yaqoob2,5

1 - Department of Physics and Astronomy, Leicester University, Leicester LE1 7RH, UK
2 - Department of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA
3 - Astrophysics Group, School of Physical and Geographical Sciences, Keele University, Keele, Staffordshire ST5 5BG, UK
4 - INAF - Osservatorio Astronomico di Brera, via Brera 28, 20121 Milan, Italy
5 - Astrophysics Science Division, Code 662, NASA/Goddard Space Flight Center, Greenbelt, MD 20771, USA
6 - INAF - Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 50125 Florence, Italy
7 - Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA

Received 12 December 2008 / Accepted 27 April 2009

Abstract
We discuss a long Suzaku observation of IRAS 19254-7245 (also known as the Superantennae), one of the brightest and well studied Ultra Luminous Infrared Galaxies in the local Universe.
This long observation provided the first detection of IRAS 19254-7245 above 10 keV, and a measurement of a 15-30 keV flux of $\sim$ $5\times
10^{-12}$ erg cm-2 s-1. The detection above 10 keV has allowed us to determine, for the first time, the intrinsic luminosity of the AGN hosted in IRAS 19254-7245, which is strongly absorbed ( $N_{\rm H}\sim 3\times10^{24}$ cm-2) and has an intrinsic luminosity (L(2-10 keV) $\sim$ $3\times
10^{44}$ erg s-1) in the QSO regime.
The 2-10 keV spectrum of IRAS 19254-7245 is remarkably hard ( $\Gamma\sim 1.2$), and presents a strong iron line ( $EW \sim 0.7$ keV), clearly suggesting that below 10 keV we observe only reprocessed radiation. Since the energy of the Fe K emission is found to be $\sim$6.7 keV, consistent with He-like Fe, and its EW is too high to be explained by a starburst-dominated scenario, we suggest that the 2-10 keV emission of IRAS 19254-7245 is dominated by reflection/scattering from highly ionized matter. Within this latter scenario, we found that the photon index of the illuminating source is $\Gamma=1.87^{+0.11}_{-0.28}$, in excellent agreement with the mean value found for radio quiet unobscured AGN.

Key words: galaxies: active - galaxies: individual: IRAS 19254-7245- galaxies: Seyfert - X-rays: galaxies

1 Introduction

Ultra luminous infrared galaxies (hereafter ULIRGs) are an enigmatic class of sources which emit most of their energy in the far-infrared (FIR, 8-1000  $\mu
{\rm m}$) domain (Sanders & Mirabel 1996), with luminosities above $\sim$1012 $L_\odot$, (i.e., comparable to QSO luminosities). The importance of understanding the physical processes at work in ULIRGs is strengthened by the observational evidence that they are generally advanced mergers of gas-rich galaxies; these events are now considered to be at the origin of some of the massive elliptical and S0 galaxies (Hopkins et al. 2006,2005; Springel et al. 2005) and the QSO stage could be a phase in the evolution of these systems. However, understanding their physical nature is complicated by the large amount of obscuration from dust present in these sources, which makes it difficult to directly observe the nuclear source.

X-ray observations of ULIRGs performed with XMM-Newton (Franceschini et al. 2003; Braito et al. 2003), Chandra (Teng et al. 2005; Ptak et al. 2003), and more recently Suzaku (Teng et al. 2009) have proved to be a fundamental tool in distinguishing between the contributions of starburst and AGN activity and to investigate the presence of hidden AGNs in these sources. These observations have shown that ULIRGs are intrinsically faint X-ray sources, their observed X-ray luminosities being typically L(2- $10~{\rm keV})\lesssim
10^{42}$-1043 erg s-1. The X-ray spectra of ULIRGs are complex and exhibit characteristics of both starburst and AGN activity, confirming the composite nature of ULIRGs. These studies have also shown that more than half of the local brightest AGN-ULIRGs (5/8) host an obscured AGN, with three being Comption Thick ( $N_{\rm H} >
10^{24}$ cm-2; NGC6240, Vignati et al. 1999; Mrk 231, Braito et al. 2004; UGC5101, Imanishi et al. 2003). Observations above 10 keV are thus fundamental for measuring the intrinsic X-ray luminosity of obscured AGN hosted in ULIRGs and its contribution to their high observed FIR emission.

IRAS 19254-7245 (also known as the SuperAntennae) belongs to a flux limited sample at $60~\mu$m consisting of the 15 brightest nearby ULIRGs (Genzel et al. 1998). Located at z=0.062, IRAS 19254-7245 has an infrared luminosity of $L_{\rm 8-1000~\:\mu\mbox{\scriptsize m}}=1.1\times 10^{12}~L_\odot$ corresponding to a bolometric luminosity of $L_{{\rm bol}}\sim 4\times 10^{45}$ erg s-1. Like most ULIRGs, IRAS 19254-7245 is a merger system of two gas-rich spiral galaxies.

The southern nucleus, optically classified as a Seyfert 2, is one of the brightest nearby ULIRGs, which has been found to host both a powerful starburst and an obscured AGN, while there is no evidence of AGN activity in the northern nucleus.

A previous X-ray observation of IRAS 19254-7245 performed with XMM-Newton suggested that this ULIRG harbors a heavily obscured and high-luminosity AGN. The hard power-law continuum above 2 keV (photon index $\Gamma
=1.3$) and the detection of a strong Fe-K$\alpha$ emission line at $6.5\pm 0.1$keV ( $EW \sim 1.4$ keV) were indeed highly indicative of a Compton-thick source (Braito et al. 2003). Since the two nuclei are located only $\sim$9 arcsec apart from each other, XMM-Newton was unable to resolve them; however, the centroid of the hard X-ray emission was found to be spatially coincident with the southern nucleus. A Chandra observation that could resolve this issue has not yet been performed.

The best-fit model obtained for the 0.5-10 keV X-ray emission detected with XMM-Newton consisted of by a strong soft thermal component associated with starburst emission and a hard X-ray component associated with AGN activity. This AGN component was parametrized by a Compton thick AGN model, and consisted of a pure Compton-reflected continuum (with $\Gamma\sim 1.8$), a scattered power law component and a strong Fe emission line. The observed 2-10 keV luminosity of the AGN was found to be $\sim$ $4\times 10^{42}$ erg s-1. Due to the limited energy bandpass of XMM-Newton, this observation did not allow us to directly see the intrinsic continuum, and thus to measure the absorbing column density and the intrinsic X-ray luminosity of IRAS 19254-7245.

Here we present the analysis of a deep Suzaku observation ($\sim$150 ks) of this system, which allowed us for the first time to constrain the intrinsic power of the AGN hosted by IRAS 19254-7245, as well as to investigate in detail the properties of the Fe line complex. In Sect. 2 we present the Suzaku data analysis and results, while in Sect. 3 we discuss the overall scenario for X-ray emission of IRAS 19254-7245. Throughout this paper, the current standard cosmology is assumed with H0=73 km s-1 Mpc-1, $\Omega_{\rm {M}}=0.27$, and $\Omega_\lambda=0.73$.

2 Observations and data reduction

Suzaku (Mitsuda et al. 2007) is the fifth Japanese X-ray satellite, which carries on board four sets of X-ray mirrors, with a X-ray CCD (XIS; three front-illuminated, FI, and one back-illuminated, BI Koyama et al. 2007) at their focal plane, and a non imaging hard X-ray detector (HXD, Takahashi et al. 2007). The latter consist of 2 main instruments: the Si PIN photodiodes and the GSO scintillator counter. Altogether, the XIS and the HXD-PIN cover the 0.5-10 keV and 12-70 keV bands, respectively.

Suzaku observed IRAS19254-7245 for a total exposure time of about 150 ksec; the observation was performed at the beginning of November 2005, when all 4 XIS were still working[*].

Cleaned event files from the version 2 of the Suzaku pipeline processing were used with the standard screening[*]. The net exposure times are 97.9 ks for each of the XIS and 142.1 ks for the HXD-PIN. The XIS source spectra were extracted from a circular region of 2.9' radius (which correspond to an energy encircled fraction of 90%[*]) centered on the source. Background spectra were extracted from two circular regions of 2.4' radius offset from the source and the calibration sources. Using the latest calibration files available, the XIS response (rmfs) and ancillary response (arfs) files were produced with the ftools tasks xisrmfgen and xissimarfgen, respectively. The net 0.5-10 keV count rates were: $(1.67\pm
0.07)\times 10^{-2}$ cts/s, $(1.57\pm 0.06)\times 10^{-2}$ cts/s, $(1.55\pm 0.06)\times
10^{-2}$ cts/s and $(1.84\pm 0.09)\times 10^{-2}$ cts/s for the XIS0, XIS2, XIS3, and XIS1 respectively. The source spectra from the three FI CCDs were then combined, while the BI (the XIS1) spectrum was kept separate and fitted simultaneously. The net XIS source spectra were then binned with a minimum S/N of 4 in each energy bin and $\chi^2$ statistics have been used.

2.1 HXD-PIN data reduction

For the HXD-PIN data reduction and analysis, we followed the latest Suzaku data reduction guide (the ABC guide Version 2)[*]. For the analysis, we used the rev2 data, which include all 4 cluster units, and the backgrounds available (Fukazawa et al. 2009), which account for the instrumental background (NXB; Takahashi et al. 2007; Kokubun et al. 2007). We then simulated a spectrum for the cosmic X-ray background counts (Gruber et al. 1999; Boldt 1987) and added this to the instrumental background.

At the time of writing, two different instrumental background files have been released (background A or ``quick'' background, and the background D or ``tuned'' background; Mizuno et al. 2008[*], and Fukazawa et al. 2009). We tested both the instrumental backgrounds and we included a $\pm 10\%$ uncertainty in the level of the CXB. The inspection of the IRAS 19254-7245 net spectrum shows that the source is detected in the 15-30 keV band with both the two background files. The net count rate in the 15-30 keV using background A and D are respectively $1.59\pm0.14\times10^{-2}$ cts s-1 and $1.48\pm0.14\times10^{-2}$ cts s-1, and the corresponding background count rates are $0.25\pm 0.004$ cts s-1 and $0.26\pm 0.003$ cts s-1.

We decided to use the latest release (background D), which is affected by lower systematic uncertainties (of about 1.3% at 1$\sigma$), which correspond to about half of the first release of the NXB. Using this background, IRAS 19254-7245 is detected in the 15-30 keV band at $\sim$5.5% above the background (a total of $\sim$2000 net counts were collected), corresponding to a signal-to noise ratio $S/N\simeq 10.8$. The dominant component in the background is the instrumental one with a count rate of $0.24\pm 0.001$ cts s-1, while the CXB count rate ranges from $1.4 \times 10^{-2}$to $1.6 \times 10^{-2}$ when we include the $\pm 10\%$ uncertainty in its level. If we assume a 10% higher CXB level, the source is still detected at 5.0% above the background (its mean count rate in the 15-30 keV is $1.35\pm0.14\times10^{-2}$ cts s-1) with $S/N\simeq 9.7$, thus the detection of IRAS 19254-7245 is not dependent on the CXB absolute level. As a further check of the level of the CXB, we analysed the Suzaku observation of the Lockman Hole performed in May 2007. We performed an identical analysis of the CXB HXD-PIN observation as we completed for IRAS 19254-7245. The flux of the CXB measured with the Lockman Hole observation was F(15- $50~{\rm keV})=
1.1\pm0.2\times 10^{-11}$ erg cm-2 s-1 (corresponding to a flux density of $3\times 10^{-11}$ erg cm-2 s-1 deg-2), in agreement with the flux of the simulated CXB (F(15- $50~{\rm keV})=1.0\pm0.2\times
10^{-11}$ erg cm-2 s-1) and the flux measured with both BeppoSAX (Frontera et al. 2007) and Swift (Moretti et al. 2009).

Since the HXD-PIN is a non-imaging detector, and taking into account the large field of view of the instrument ( $0.56~{\rm deg} \times 0.56~{\rm deg}$), we checked that the detection is not caused by to another X-ray source. In particular, we searched the NED database for known AGN in the HXD field of view and inspected the available XMM-Newton observation. We found that two X-ray sources with a 2-10 keV flux comparable to IRAS 19254-7245's emission were detected with XMM-Newton. The two sources are both AGNs, belonging to the XMM Bright Survey sample (XBS J193138.9-725115 and XBS J193248.8-723355;; Caccianiga et al. 2008). XBS J193138.9-725115 is a Type 1 AGN (z=0.701) and its XMM-Newton spectrum is well modeled by a single unabsorbed power law component ( $\Gamma\sim 2$) with no evidence of absorption. XBS J193248.8-723355 is a Seyfert 2 at z=0.287; the X-ray emission of this source agrees with its classification as a Compton-Thin Seyfert; indeed a low energy cut-off is present in the XMM-Newton spectra corresponding to $N_{\rm H}\sim 10^{22}$cm-2 and there is no evidence that this source could be Compton thick. The predicted 15-30 keV emission from these sources (derived from the analysis of the XMM-Newton and Suzaku data) is less than $\sim$10-13 erg cm-2 s-1, which is below the HXD-PIN sensitivity and a factor of $\sim$50 below the measured 15-30 keV flux.

For the spectral analysis, we rebinned the HXD-PIN spectrum of IRAS 19254-7245 to have a signal-to-noise ratio of 5 in each energy bin. To derive a first estimate of the 15-30 keV flux and luminosity of IRAS 19254-7245, we fitted the HXD-PIN spectrum assuming a power law model with $\Gamma=1.9$ (i.e., a standard AGN value; Reeves & Turner 2000; Page et al. 2004; Caccianiga et al. 2004). Taking into account the systematic uncertainties in the NXB model, with this simple model we obtained F(15- $30~{\rm keV})\sim 5.2\pm 1.1\times 10^{-12}$ erg cm-2 s-1and L(15- $30~{\rm keV})\sim
4.7\times 10^{43}$ erg s-1. The extrapolation of this model down to 2 keV predicts an intrinsic luminosity (which does not include a correction for Compton scattering) of the AGN in IRAS 19254-7245 of L(2- $10~{\rm keV})\sim 9.5\times 10^{43}$ erg s -1.

2.2 The broad-band spectrum

Overall the Suzaku observation confirms the XMM-Newton results. A good fit to the 0.5-10 keV Suzaku data is obtained with a model consinting of: a thermal emission component ( $kT=0.64\pm 0.10$ keV and abundance $Z=Z_\odot$, likely associated with starburst activity), a strong hard power-law component ( $\Gamma=1.2\pm 0.1$, likely associated with AGN emission), and a strong iron K emission line ( EW= 710+190-170 eV, with respect to the observed continuum). The flux and observed luminosity F(2-10 keV) $\sim 2.9\times
10^{-13}$ erg cm-2 s-1, and L(2-10 keV) $\sim$ $2.5\times 10^{42}$ erg s-1are found to be consistent with the values measured with XMM-Newton. The de-absorbed luminosity of the starburst component is L(0.5- $2~{\rm keV}) \sim 4\times 10^{41}$ erg s-1  also in agreement with the luminosity measured by XMM-Newton.

We then tested the best-fit model obtained for the XMM-Newton spectrum. In this model, the soft X-ray emission is still modeled with a thermal component, while the hard ( $\Gamma\sim 1.2$) power law component is replaced by a pure Compton-reflected continuum (the PEXRAV model in Xspec, Magdziarz & Zdziarski 1995, with an intrinsic $\Gamma =1.8$) combined with a moderately absorbed ( $N_{\rm H}\sim 5\times 10^{21}$ cm-2) power law component with the same $\Gamma$, representing the possible scattered emission. The parameters of the reflection component are: an inclination angle $i=45^\circ$, abundance $Z=Z_\odot$ (using the abundances of Wilms et al. 2000), and a reflection fraction (defined by the subtending solid angle of the reflector $R=\Omega/2\pi$) R fixed to be 1. The normalization of this component was allowed to vary. This model is a good fit to the XIS data alone ( $\chi^2/{\rm d.o.f.}=159/144$, see Fig. 1) and the fluxes and observed luminosity of the AGN components are consistent with the values measured with the XMM-Newton observation. However, this model clearly underpredicts the counts detected above 10 keV (see Fig. 1, green data points). When we include in the fit the HXD-PIN data, fixing the cross-normalization between the XIS and the PIN to 1.16 (Manabu et al. 2007; Maeda et al. 2008[*]), the model is statistically unacceptable ( $\chi^2/{\rm d.o.f.}=231/147$) and even allowing for a harder photon index ( $\Gamma\sim 1.3$) it is not a good representation of the 0.5-30 keV emission ( $\chi^2/{\rm d.o.f.}=220/146$).

 \begin{figure}
\par\includegraphics[width=6.15cm,clip,angle=-90]{11516fg1.eps}
\par\end{figure} Figure 1:

0.5-30 keV Suzaku XIS-FI (black data points), XIS-BI (red), and HXD-PIN (green) data and their ratio to the 0.5-10 keV best-fit model. The model is composed by: a thermal emission component dominating the 0.5-2 keV emission, a Compton-reflected continuum (R=1, $\Gamma =1.8$) and a strong Fe emission line at 6.7 keV ( $EW\sim 700$ eV). A clear excess is present above 10 keV, which we attribute to the intrinsic X-ray emission of IRAS 19254-7245 transmitted through the high column density absorber.

Open with DEXTER

To account for the excess detected above 10 keV, we added to the model a second heavily absorbed power law component. Since the HXD/PIN residuals suggest the presence of a high column density absorber, we used for this component the model by Yaqoob (1997) ( PLCABS in Xspec), which correctly takes into account Compton down-scattering. For high column densities the observed X-ray continuum may also be suppressed by Compton down-scattering, and the intrinsic luminosity must be corrected by a factor $e^{\tau}$, where $\tau=N_{\rm H}\sigma_\tau$ and $\sigma_\tau = 6.65 \times 10^{-25}$ cm2 is the Thomson cross-section.

Table 1:   Results of the spectral fit.

 \begin{figure}
\par\includegraphics[width=6.5cm,angle=-90,clip]{11516fg2.eps}\par\includegraphics[width=6.5cm,angle=-90,clip]{11516fg3.eps}
\end{figure} Figure 2:

Upper panel: Suzaku XIS and HXD spectra of IRAS 19254-7245 where the underlying AGN continuum is modelled with a neutral Compton-reflected component and power-law component seen in transmission through a high column density absorber ( $N_{\rm H}\sim 4\times 10^{24}$ cm-2). Lower panel: same as above but with the 2-10 keV AGN emission modelled with an ionized reflected component (the high energy component is modelled as in the upper panel with a power-law component seen in transmission through the high column density absorber).

Open with DEXTER

This model provides a good fit to the 0.5-30 keV spectrum ( $\chi^2/\rm {d.o.f.}=164/145$; see Fig. 2 upper panel and Table 1 model A). However, we found a low value of the reflection fraction with respect to this primary absorbed power law component (R<0.1). This suggests that the line and the hard 2-10 keV spectrum are unlikely to be produced by reflection off cold material. The broad band continuum could also be reproduced by a model without the reflected component (see Table 1 model B). Statistically, this model provides a slightly worse fit ( $\chi^2/\rm {d.o.f.}=181/146$) than the previous one, but it is unable to account for the hardness of the 2-10 keV emission. In particular, clear residuals are present in the 5-10 keV band where the reflected component dominates (see Fig. 2 top panel). Finally, if we allowed the photon index to vary, we found that although the fit improves we again need an unusually hard photon index of the power law component ( $\Gamma=1.20^{+0.11}_{-0.04}$, $N_{\rm H}=4\pm 1\times 10^{24}$ cm-2; see Table 1 model C).

One of the main results of this observation is that, although we can confirm the presence of a strong Fe line as detected with XMM-Newton, the centroid of this line is now at $ 6.67\pm 0.05$ keV in the rest frame ( $EW \sim 0.7$ keV; for the model without the reflected component), consistent with He-like Fe. Furthermore, the line appears to be marginally broad ( $\sigma=0.12\pm 0.06$ keV). The inclusion of the line in the model improves the fit by $\Delta\chi^2= 51$ for 3 degree of freedom[*]. However, if we constrain the line to be unresolved the fit is worse only by $\Delta\chi^2=5$. To check the energy and the intrinsic width of the Fe line detected in the XIS, we examined the spectra of the 55Fe calibration source lines, which are located at two corners of each XIS. The calibration source produces lines from Mn $
{\rm K}\alpha_{1}$ at 5.899 keV and Mn $ {\rm K}\alpha_{2}$ at 5.888 keV. From the spectrum of the calibration source, we found that the line energy is shifted redwards by about 25 eV, while the residual width is $\sigma \sim 50 $ eV[*], confirming that the broadening of the line is intrinsic to the source and not instrumental[*]. After the subtraction in quadrature of this residual width, we obtained $ \sigma_{{\rm int}}^{2} =110\pm60 $ eV). However, taking into account the present count statistics of the data, this broadening could be due to the presence of other line components, which are not resolved. In particular, the line profile can be explained with three unresolved Gaussian lines (at 6.4 keV, 6.7 keV, and 6.96 keV; see Fig. 3). Statistically, this model provides a similarly good fit as models with a single broad line ( $\chi^2/\rm {d.o.f.}=161/145$), the strongest line being the 6.7 keV line ( $EW\sim 400$eV). Although the other two lines are not statistically required, a weak (EW< 200 eV) emission line could be present at the energy of the neutral Fe K$\alpha$, while the 90% upper limit to the 6.96 keV line is 150 eV.

 \begin{figure}
\par\includegraphics[width=6.7cm,angle=-90,clip]{11516fg4.eps} %
\par\end{figure} Figure 3:

Residuals of the data/model of IRAS 19254-7245 XIS data at the Fe band, when no iron line is included in the model. The three vertical lines indicate the energy centroids of the three possible components of the Fe line complex (6.4 keV, 6.7 keV and 6.97 keV). The energy scale is in the rest frame.

Open with DEXTER

Since the energy centroid of the line detected by Suzaku appears to disagree with the XMM-Newton results and the EW appears to be lower, we reanalyzed the XMM-Newton data and compared them with the Suzaku results. The exposure time in the XMM-Newton observation was only 20 ks, and after taking into account the errors in the flux and line continuum, we found that the two lines are consistent with each other. Furthermore, the energy centroids are consistent to within the errors ( $E_{\rm XMM}=6.5\pm0.1$ keV; $E_{{\rm Suzaku}}=6.66\pm 0.05$ keV). Finally, a possible blending of 3 lines was also present in the XMM-Newton data, but again the low exposure time of this observation did not allow a more detailed analysis of the Fe line profile.

From a statistical point of view, all these models are a good representation of the 0.5-30 keV emission, but they are unable to account for the hardness of the continuum. Furthermore, the line energy of the strongest emission line is at odds with a scenario where the 2-10 keV emission is dominated by reflection/scattering off cold material as assumed by the continuum model tested above. One possibility is that the 6.7 keV line is due to reflection from highly ionized matter; we thus replaced the cold reflected ( PEXRAV) power law component with an ionized reflected component, as described by the Ross & Fabian (2005) table (otherwise known as the REFLION model). This model allows different values for the ionization parameter of the reflecting material and it also includes the Fe K emission line, as well as emission lines from other elements in addition to the reflected continuum. We fixed the Fe abundance to be solar and we included a lower column density in front of the reflected component. The photon index of the illuminating X-ray source was kept as a free parameter. A good fit (see Fig. 2 lower panel and Table 1 model D) to the 0.5-30 keV emission was obtained with a ionization parameter of $\xi=
1000^{+170}_{-430}$ erg cm s-1, where the value of the ionization was determined mainly by the strength and energy of the Fe line; at this ionization level, Fe K emission is almost entirely due to Fe xxv. The reflected component is modified by a lower column density absorber with $N_{\rm H}\sim 10^{22}$ cm-2, which is probably on a larger scale than the inner high column density absorber. We emphasize that a second highly absorbed ( $N_{\rm H}\sim 3\times10^{24}$ cm-2) power-law component is still required to account for the HXD-PIN emission and the intrinsic 2-10 keV luminosity is $\sim$ $3\times
10^{44}$ erg s-1. The photon index of the illuminating source becomes $\Gamma=1.87^{+0.11}_{-0.28}$, consistent with the mean value measured in unobscured radio quiet AGN ( $\Gamma_{\rm mean}=1.9$; Reeves & Turner 2000).

In summary, this model is now able to reproduce in a consistent way all the main characteristic of the the broad-band X-ray emission of IRAS 19254-7245 and in particular the Fe emission line detected at $\sim$6.7 keV and the flatness of the observed continuum. Finally, it is worth noting that, independently of the assumed model for the 0.5-10 keV emission, we always need to include an absorbed power law component to account for the HXD-PIN data with a $N_{\rm H}\sim$ 3-4  $\times 10^{24}$cm-2, and the derived 2-10 keV intrinsic luminosity is always above L(2-10 keV) $\sim 10^{44}$ erg s-1, ranging from $2\times 10^{44}$ erg s-1to $4\times 10^{44}$ erg s-1 (see Table 1).

3 Discussion and conclusions

The detection of the Fe K emission line at 6.7 keV instead of the 6.4 keV emission line expected from neutral iron could suggest that this line is associated with strong starburst activity and that the emission below 10 keV is not caused by the AGN, but rather by a hot thermal plasma as expected in a starburst-dominated scenario. From a statistical point of view, we can obtain a good fit ( $\chi^2/d.o.f.=147/144$) by replacing the AGN reflected emission with a thermal component. This model gives a best-fit temperature of kT=8.1+1.2-1.3 keV, $N_{\rm H}\sim 6\times 10^{21}$ cm-2, twice solar abundance ( $Z\sim
2.1~Z_\odot$), and a luminosity of L(2- $10\; \rm {keV})\sim 2\times 10^{42}$ erg s-1.

A possible origin of the high temperature ($kT\sim 8$ keV) plasma could be the presence of several supernovae (SNe); their X-ray emission could in principle explain the high temperature as well as the presence of the strong 6.7 keV line (Persic & Rephaeli 2002). From the FIR luminosity, we estimated a star formation rate (SFR) for IRAS 19254-7245 of $\sim$200 $M_\odot$ yr-1 and a supernova (SN) rate of $\sim$2 SNe yr-1(Mannucci et al. 2003). However, even assuming that the SNe are at the higher end of the expected range of X-ray luminosity ( $L_{\rm x}=10^{40}$-1041 erg s-1), a factor of 10 more SNe than the predicted rate would be required to maintain the observed hard X-ray emission. Furthermore, X-ray observations of starburst (SB) galaxies showed that the major contributor to the 2-10 keV emission is the integrated emission from high-mass X-ray binaries (HMXB) and that there is a linear relation between the SFR and the 2-10 keV emission from HMXB. In the case of IRAS 19254-7245, if all the FIR luminosity is produced by the SB activity, then the predicted 2-10 keV luminosity of the HMXB is $\sim$1042 erg s-1 (Persic et al. 2004; Grimm et al. 2003; Ranalli et al. 2003). In this scenario, we would thus require that, in contrast to other SB, the hot diffuse emission has a luminosity similar to the contribution from HMXB. In summary, although the starburst model can well reproduce the line intensity and the overall shape of the 2-10 keV continuum, the luminosity of this thermal component (L(2- $10~{\rm keV})\sim 2 \times 10^{42}$  rg s-1) is likely to be too high to be consistent with a pure starburst scenario.

Finally, we note that if we assume that the emission detected above 2 keV is dominated by the emission of unresolved HMXB we still cannot explain the 6.7 keV emission line. We would expect a lower EW of the Fe line ( $EW\sim 0.3$ keV; Persic & Rephaeli 2002; White et al. 1983), which is inconsistent with the high value observed in IRAS 19254-7245. One possibility is that only a fraction of the Fe line at 6.7 keV originates in a high temperature plasma. We thus included in our best-fit model a second high temperature emission fixing the abundance to the solar value. This high temperature thermal emission can account for $\sim$$30\%$ of the flux of the line at 6.7 keV. Although the inclusion of this emission can account for a fraction of the Fe line at 6.7 keV, we still need a strong ionized reflection component to account at the same time for the continuum shape and the Fe line.

On the other hand, if we attribute the hard X-ray emission to the presence of the AGN the flatness of its observed continuum together with its strong Fe emission lines are usually considered evidence of a Compton-Thick source, where no direct emission is seen below 10 keV and the shape is produced by reflection off cold matter. However, in contrast to what observed for other Compton thick sources, here the Fe K emission is predominately in the He-like state with little or no 6.4 keV line.

It is often assumed that a hallmark of Compton-thick X-ray sources is a prominent Fe K$\alpha$emission line at $\sim$6.4 keV, from neutral matter, with an EW that can exceed a keV. However, we show below that in practice, the Fe K$\alpha$ line EW can be reduced significantly and become undetectable for large column densities of the order of $10^{25} \
\rm cm^{-2}$ or more. The observed EW of the Fe K line is at its largest when only the reflected continuum associated with the production of the line is observed. However, if any additional continuum is observed and if its magnitude is comparable to, or greater, than the scattered/reflected continuum, the Fe K line will be ``diluted'' and the EW reduced. This is because for lines of sight that intercept the reprocessor, the absolute luminosity of the 6.4 keV Fe K line and its associated reflected continuum can be much less than the continuum from the Compton-thin scattering zone, so that the EW of the line measured against the total observed continuum can be much less than that measured against the continuum emerging from the Compton-thick reflector.

In this scenario, the ionized Fe emission line could be produced in the same Compton-thin zone that scatters a fraction of the intrinsic continuum into the line-of-sight. If this scattered continuum is the dominant observed continuum (as would be the case if the direct line-of-sight is obscured by Compton-thick matter), then the EW of the emission lines from ionized Fe could be large, of the order of hundreds of eV or more. The details depend on several factors, principally the ionization parameter and column density of the line-emitting region, as well as the shape of the ionizing continuum (see for example, detailed calculations of the EW of the Fe  XXV and Fe  XXVI emission lines in Bianchi & Matt 2002). Suppose that the Compton-thick reprocessor subtends a solid angle $\Delta \Omega/4\pi$at the X-ray source and that the system is observed along a line-of-sight that does not provide a direct view of the X-ray source and has the greatest column density, NH, (to maximize the Fe K line EW). Further suppose that the space between the reprocessor is filled by a warm, optically-thin scattering zone with Thomson depth $\tau_{\rm thin}$ (with $\tau_{\rm thin} \ll 1 $), subtending a solid angle $1-(\Delta \Omega/4\pi)$at the X-ray source. Then, a fraction $f\equiv \tau_{\rm thin}[1-(\Delta \Omega/4\pi)]$ of the intrinsic X-ray continuum is scattered into the observer's line-of-sight and will reduce the EWof the Fe K line if it dominates over the zeroth-order continuum that is observed directly through the Compton-thick absorber. Specifically, for column densities greater than a few $\times
10^{24} \ \rm cm^{-2}$, even a small value of f can significantly reduce the EW of the Fe K line. For example, for a column density of $10^{25} \
\rm cm^{-2}$, a scattering fraction of f>0.01 will reduce the EW of an Fe K line by more than an order of magnitude, so that an EW of 1 keV would be reduced to less than $\sim$100 eV, and it could render the line undetectable (see Ghisellini et al. 1994).

Furthermore, the intrinsic EW (i.e., prior to the dilution effect) of Fe K lines depends on several factors, not only the column density of the absorber but also the geometry of the absorber (e.g., the half opening angle of the putative torus; see Matt et al. 1996; Ikeda et al. 2009; Ghisellini et al. 1994). For example, for an half opening angle of $30^{\circ}$ and our estimate of the column density of the neutral absorber, the intrinsic EW of the 6.4 keV Fe line can span the range from 1 to 4 keV (Ghisellini et al. 1994).

The inferred column density of the Compton-thick reprocessor implies that a scattering fraction of only $\sim$$0.1\%$ in the optically-thin zone is required to begin to dilute the 6.4 keV Fe K line and a scattering fraction of a few percent is sufficient to reduce the EW of the line well below 100 eV, consistent with the upper limit to the EW and the $\sim$2% scattering fraction (measured with respect to the de-absorbed primary power law component) as measured with the Suzaku data.

A second possible geometry is that we have a direct view of the inner surface of the Compton-thick reprocessor, but the outer part of this reprocessor is ionized. In this case, if the remaining part of the reprocessor is Compton-thin, the EWof the 6.4 keV Fe K$\alpha$ line will be reduced and the emission detected below 10 keV is the reflected emission from this inner ionized surface of the torus, which will also produce a strong 6.7 keV line. Thus, we see that the lack of a large EW neutral Fe K in IRAS 19254-7245 is not unexpected. The spectrum below 10 keV is then dominated by this optically-thin scattered continuum and the dominance of the emission line from ionized Fe is consistent with this picture.

We note that this is not a unique case of a detection of a strong 6.7 keV line in a luminous infrared galaxy. Other examples are Arp299 (Ballo et al. 2004), Arp220 (Teng et al. 2009; Iwasawa et al. 2005), and IRAS 00182-7112 (Nandra & Iwasawa 2007). For all these sources although the optical spectra show no clear signature of AGN activity, their X-ray emission and the 6.7 keV line can be explained by the presence of an AGN and an ionized reflector as in the case of IRAS 19254-7245. However, while in the case of Arp 299 and Arp 220, the X-ray luminosity is not indicative of the major contributor to the bolometric luminosity being a high-luminosity AGN, in the case of IRAS 00182-7112 the X-ray luminosity is too high to be accounted for by the strong starforming activity (L(2- 10)> 1044 erg s-1) as for IRAS 19254-7245. For all these sources, the presence of a strong ionized Fe line, with little or no 6.4 keV line, could be reconciled with the picture of a heavily obscured AGN assuming that we do not have a direct view of the reflected continuum; this is because the continuum is diluted by the scattering from the ionized matter, which is associated with the production of the line, or if the surface of the putative Compton-thick reprocessor is highly ionized.

Overall to account for the X-ray emission of IRAS 19254-7245 above 2 keV, we need two absorbing/reflecting media: one neutral and Compton-thick and probably seen in transmission, and one ionized and seen probably in reflection. This latter is responsible for the flat X-ray spectrum emerging below 10 keV, for the He-like Fe K line and probably for the dilution of the 6.4 keV Fe line produced in the neutral Compton-thick absorber.

Despite the various possible models for the 0.5-10 keV emission, we always need a neutral high column density absorber to account for the emission emerging above 10 keV. Once we have corrected for the amount of absorption, the intrinsic 2-10 keV luminosity of the primary AGN component is $\sim$1044 erg s-1. This high column density absorber may also be the one responsible for the deep hydrocarbon absorption detected in the L-band spectrum (Risaliti et al. 2003) at 3.4  $\mu
{\rm m}$. Finally, independently of the assumed model for the 2-35 keV emission (i.e., ionized reflection or scattered power-law component), we always require the soft thermal component to account for the 0.5-2 keV emission. As already found for other ULIRGs, this thermal component has a temperature $kT\sim 0.7$ keV, and the luminosity of this component is L(0.5-2 keV)  $\sim 4\times 10^{41}$ erg s-1. Although, we cannot exclude a possible contribution from the ionized reflector, the measured soft X-ray luminosity agrees with that expected from the FIR luminosity and the SFR of IRAS 19254-7245 (Persic et al. 2004; Ranalli et al. 2003).

The total X-ray luminosity estimated from the intrinsic component at E>10 keV can be converted into a bolometric luminosity and compared with the total infrared emission. We adopted the $\alpha_{OX}$-luminosity correlation of Steffen et al. (2006) to estimate the 2500 Å luminosity, and the Elvis et al. (1994) quasar spectral energy distribution to estimate the total luminosity of the AGN component. Assuming $L_{2-10}\sim10^{44}$ erg s-1, we obtained $L_{\rm BOL}(AGN)\sim2\times10^{45}$ erg s-1, i.e., about 50% of the infrared luminosity. A different way of estimating the AGN luminosity is from its emission in the mid-infrared, the only other band where the continuum emission of the AGN is not completely absorbed. From an analysis of the Spitzer-IRS spectrum, Nardini et al. (2008) estimated an AGN contribution to the bolometric luminosity of IRAS 19254-7245 of $\sim$25%. Considering the uncertainties in the bolometric corrections in both the X-ray and the mid-infrared bands, the two estimates can be considered to be in rough agreement. If the difference is assumed to be real, this could be an indication of the non-complete covering factor of the AGN circumnuclear absorber: indeed, the estimate from the mid-infrared spectrum is done assuming a complete reprocessing of the intrinsic AGN emission in the infrared. However, the optical classification of IRAS 19254-7245 as a Seyfert 2 suggests that the obscuration of the nuclear source is not complete. The estimates from the X-ray and infrared spectra could then be perfectly reconciled by assuming a covering factor of the obscuring material of about 50%.

In conclusion, this deep Suzaku observation has allowed us to measure for the first time the hard X-ray emission of IRAS 19254-7245 and infer that its intrinsic 2-10 keV luminosity is about few times 1044 erg s-1. We have found evidence that the AGN hosted in IRAS 19254-7245 is highly obscured, with a measured column density of the neutral absorber of $N_{\rm H}\sim 3\times10^{24}$ cm-2. We confirm the presence of a strong iron K emission line with an $EW \sim 0.7$ keV. The energy of iron K emission line is found to be consistent with Fe XXV. We propose that the X-ray emission detected below 10 keV can be ascribed to scattered/reflected emission from highly ionized matter, which could be identified with the warm Compton-thin gas that fills the space between the neutral Compton-thick reprocessor.

Acknowledgements
V.B. acknowledge support from the UK STFC research council. R.D.C. acknowledge financial support from the ASI (Agenzia Spaziale Italiana) grant I/088/06/0. Support for this work was provided by the National Aeronautics and Space Administration through the NASA grant NNG04GB78A. We thank the anonymous referee for his/her useful comments, which have improved this paper.

References

Footnotes

... working[*]
Later that month Suzaku XIS2 failed. No charge injection (see: http://suzaku.gsfc.nasa.gov/docs/suzaku/analysis/sci.html) was applied at the time of the observation so the nominal energy resolution of the XIS at 6 keV was degraded with respect to the prelaunch one
... screening[*]
The screening filter all events within the South Atlantic Anomaly (SAA) as well as with an Earth elevation angle (ELV) < $5\ensuremath {{}^{\circ }} $ and Earth day-time elevation angles (DYE_ELV) less than $ 20\ensuremath {{}^{\circ }} $. Furthermore also data within 256 s of the SAA were excluded from the XIS and within 500s of the SAA for the HXD. Cut-off rigidity (COR) criteria of $ > 8 ~{\rm GV}$ for the HXD data and > $6 ~{\rm GV}$for the XIS were used.
...%[*]
see ftp://legacy.gsfc.nasa.gov/suzaku/doc/xrt/suzakumemo-2008-04.pdf
... Version 2)[*]
http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/
... 2008[*]
http://www.astro.isas.jaxa.jp/suzaku/doc/suzakumemo/suzakumemo-2008-03.pdf
... 2008[*]
http://www.astro.isas.jaxa.jp/suzaku/doc/suzakumemo/suzakumemo-2007-11.pdf; http://www.astro.isas.jaxa.jp/suzaku/doc/suzakumemo/suzakumemo-2008-06.pdf
... freedom[*]
for the scenario with the reflected component the fit improves by $\Delta\chi^2= 40$.
...$\sigma \sim 50 $ eV[*]
This residual width is due to the degradation of the XIS after the launch and prior to the correction with the charge injection.
... instrumental[*]
The intrinsic width of the Fe line in the spectrum of IRAS 19254-7245 can be $ \sigma_{{\rm int}}^{2} = \sigma_{{\rm meas}}^{2} -
\sigma_{{\rm lamp}}^{2}$ (where $ \sigma_{{\rm meas}}$ is the measured width and $\sigma_{\rm lamp}$ is the width of the calibration lines)

All Tables

Table 1:   Results of the spectral fit.

All Figures

  \begin{figure}
\par\includegraphics[width=6.15cm,clip,angle=-90]{11516fg1.eps}
\par\end{figure} Figure 1:

0.5-30 keV Suzaku XIS-FI (black data points), XIS-BI (red), and HXD-PIN (green) data and their ratio to the 0.5-10 keV best-fit model. The model is composed by: a thermal emission component dominating the 0.5-2 keV emission, a Compton-reflected continuum (R=1, $\Gamma =1.8$) and a strong Fe emission line at 6.7 keV ( $EW\sim 700$ eV). A clear excess is present above 10 keV, which we attribute to the intrinsic X-ray emission of IRAS 19254-7245 transmitted through the high column density absorber.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=6.5cm,angle=-90,clip]{11516fg2.eps}\par\includegraphics[width=6.5cm,angle=-90,clip]{11516fg3.eps}
\end{figure} Figure 2:

Upper panel: Suzaku XIS and HXD spectra of IRAS 19254-7245 where the underlying AGN continuum is modelled with a neutral Compton-reflected component and power-law component seen in transmission through a high column density absorber ( $N_{\rm H}\sim 4\times 10^{24}$ cm-2). Lower panel: same as above but with the 2-10 keV AGN emission modelled with an ionized reflected component (the high energy component is modelled as in the upper panel with a power-law component seen in transmission through the high column density absorber).

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=6.7cm,angle=-90,clip]{11516fg4.eps} %
\par\end{figure} Figure 3:

Residuals of the data/model of IRAS 19254-7245 XIS data at the Fe band, when no iron line is included in the model. The three vertical lines indicate the energy centroids of the three possible components of the Fe line complex (6.4 keV, 6.7 keV and 6.97 keV). The energy scale is in the rest frame.

Open with DEXTER
In the text


Copyright ESO 2009

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.