A&A 442, 469-478 (2005)
DOI: 10.1051/0004-6361:20053558

XMM-Newton spectral properties of the ultraluminous IRAS galaxy Mrk 273

I. Balestra - Th. Boller - L. Gallo - D. Lutz - S. Hess

Max-Planck-Institut für extraterrestrische Physik, Postfach 1312, 85741 Garching, Germany

Received 2 June 2005 / Accepted 12 July 2005

Abstract
We present a 23 ks XMM-Newton observation of the Ultraluminous Infrared Galaxy (ULIRG) Mrk 273. The hard X-ray spectrum can be modeled by a highly absorbed ($\sim$ $7\times10^{23}$ cm-2) power law plus an Fe K$\alpha $ emission line. The iron line (detected at more than 99% c.l.) is broad ( $\sigma=0.26^{+0.37}_{-0.17}$ keV), suggesting possible superposition of a neutral iron line at 6.4 keV, and a blend of ionized iron lines from Fe XXV and Fe XXVI. Given the relatively short exposure, the three line components cannot be individually resolved with high statistical significance: the neutral component is detected at ${\sim}2.5\sigma$ and the Fe XXV line at $\sim$$2\sigma$ c.l., while for the Fe XXVI line we can only estimate an upper limit. The broad band spectrum requires, in addition to a highly absorbed power law, at least three collisionally ionized plasma components, which may be associated with star-forming regions. The temperatures of the three plasmas are about 0.3, 0.8 and 6 keV, where the highest of the three is sufficient to produce ionized iron emission lines. An alternative interpretation for the origin of the soft emission might also be reflection off photoionized gas, as has been observed in a number of nearby Compton-thick Seyfert 2 galaxies (e.g. NGC 1068, Circinus, Mrk 3, NGC 4945). A hot gas, photoionized by the primary, continuum can also produce ionized iron lines. Unfortunately, given the limited statistics and the lack of high resolution spectroscopy, it is not possible to distinguish between the two models investigated. We further compare the XMM-Newton findings with the Chandra data obtaining consistent spectral results. The absorption-corrected hard X-ray luminosity of Mrk 273 is $L_{2-10~{\rm keV}}\sim7\times10^{42}$ erg s-1, corresponding to $\sim$0.2% of the far-IR luminosity, similar to typical values found in pure starbursts. The thermal contribution to the soft X-ray luminosity is approximately $0.2{-}0.7\times10^{42}$ erg s-1, comparable to those found in NGC 6240 and other starburst-dominated ULIRGs. We also analyze the XMM-Newton spectrum of Mrk 273x, an unabsorbed Seyfert 2 galaxy at redshift z=0.458, which lies in the field of view of Mrk 273.

Key words: galaxies: individual: Mrk 273 - galaxies: individual: Mrk 273x - galaxies: Seyfert - X-rays: galaxies

1 Introduction

The InfraRed Astronomical Satellite (IRAS) has detected a very large number of galaxies in the local universe (z<0.3) which exhibit extraordinarily high infrared luminosity ( $ L_{\rm IR}\ge 10^{12}~L_{\odot}$ for $ H_0 = 75~{\rm km~s}^{-1}~{\rm Mpc}^{-1}$; see Sanders & Mirabel 1996, for a review). The number density of these so called Ultraluminous InfraRed Galaxies (ULIRGs) exceeds that of optically selected Seyfert galaxies and QSOs with comparable bolometric luminosities (Sanders et al. 1988a; Soifer et al. 1987; Sanders et al. 1988b; Sanders & Mirabel 1996) by a factor of ${\sim}1.5{-}2$ (Sanders et al. 1999). The bulk of luminosity in these sources is infrared (IR) emission from warm dust. Only two dust-heating mechanisms are capable of producing such an extraordinary IR luminosity: one involves the presence of a strong starburst region, and the other a dust-enshrouded AGN. However it is still not clear what the relative contribution of each component to their bolometric luminosity should be.

Spectroscopic surveys of samples of ULIRGs carried out with the Infrared Space Observatory (ISO) have revealed that about 80% of such objects are predominantly powered by star formation, but the fraction of AGN-powered objects increases with luminosity (Genzel et al. 1998; Lutz et al. 1998).

X-ray observations are a fundamental tool to probe the highly obscured innermost regions of these objects, and therefore potentially unveil the physical processes at work. In fact a considerable portion of ULIRGs has been found to contain a hard X-ray source, highly absorbed by a molecular torus, which indicates the presence of a hidden AGN (Mitsuda 1995; Ptak et al. 2003; Braito et al. 2003; Franceschini et al. 2003; Brandt et al. 1997; Vignati et al. 1999; Kii et al. 1997).

Mrk 273 is a well studied ULIRG at a redshift z = 0.03778 with a Seyfert 2 nucleus (Sanders et al. 1988a; Koski 1978). It presents evidences for strong star formation (Goldader et al. 1995), such as strong Polycylic Aromatic Hydrocarbons (PAH) features in the near-IR (Genzel et al. 1998) and extended ($\sim$370 pc) radio emission, punctuated by a number of compact sources identifiable as supernovae or supernova remnants (Carilli & Taylor 2000; Bondi et al. 2005). Like many ULIRGs, Mrk 273 presents clear indications of an ongoing merging process such as a long tidal tail and a double nucleus (Carilli & Taylor 2000; Soifer et al. 2000; Knapen et al. 1997).

ASCA X-ray observations of Mrk 273 showed evidence of a highly absorbed ( ${\sim}4\times10^{23}$ cm-2) hard component above about 3 keV together with soft thermal emission (Iwasawa 1999; Turner et al. 1998,1997). These findings were indeed confirmed by Chandra, which revealed a compact hard X-ray nucleus inside a much more extended soft halo (Xia et al. 2002, hereafter Paper I).

In this paper we will discuss the analysis of the XMM-Newton data of Mrk 273 and compare them with the previous Chandra observation (Paper I). The data analysis and details of the XMM-Newton observation are described in Sect. 2, together with our updated reduction of the Chandra data. In Sect. 3 the spectral fitting results are presented. A comparison with the Chandra observation is given in Sect. 4. Our results are then discussed in Sect. 5 and a summary of the XMM-Newton results on Mrk 273 is given in Sect. 6. In Appendix A we discuss the XMM-Newton spectrum of Mrk 273x, an unabsorbed Seyfert 2 observed serendipitously in the field of view of Mrk 273.

Throughout this paper we use the following cosmological parameters: H0=70 km s-1 Mpc-1, $\Lambda_0=0.7$ and q0=0. Within the adopted cosmology 1'' corresponds to 0.75 kpc at the redshift of Mrk 273.

2 Observations and data reduction

2.1 XMM-Newton

Mrk 273 was observed with XMM-Newton (Jansen et al. 2001) for 23 ks on 2002 May 07 (revolution 0441). During this time the EPIC-PN (Strüder et al. 2001) and MOS (MOS1 and MOS2; Turner et al. 2001) cameras, as well as the Optical Monitor (OM; Mason et al. 2001) and the Reflection Grating Spectrometers (RGS1 and RGS2; den Herder et al. 2001) collected data. The EPIC-PN and MOS cameras were operated in full-frame mode and utilized the thick filter.

The Observation Data Files were processed to produce calibrated event lists using the XMM-Newton Science Analysis System (SAS v6.1.0). Unwanted hot, dead, or flickering pixels were removed as were events due to electronic noise. Event energies were corrected for charge-transfer losses, and EPIC response matrices were generated using the SAS tasks ARFGEN and RMFGEN. Light curves were extracted from these event lists to search for periods of high background flaring. Background flaring was negligible. The total good exposure times selected for the PN and MOS were 18 and 22 ks, respectively.

The source plus background photons were extracted from a circular region with a radius of 50'', and the background was selected from an off-source region with a radius of 70'' and appropriately scaled to the source region. Single and double events were selected for the PN detector, and single-quadruple events were selected for the MOS. The resulting PHA files were grouped with a minimum of 20 counts per bin. Pile-up effects were determined to be negligible.

The RGS were operated in standard Spectro+Q mode. The first-order RGS spectra were extracted using the SAS task RGSPROC, and the response matrices were generated using RGSRMFGEN. Unfortunately, due to the low signal-to-noise the RGS data were not useful.

The OM was operated in imaging mode for the entire observation. Fourteen images were taken in three filters: 5 in UVW1 (245-320 nm), 9 in UVM2 (205-245 nm), and 1 in UVW2 (180-225 nm). The average exposure time (800 s) was short for such a faint source and magnitudes could not be calculated for all UVM2 images. The average apparent magnitude in each filter was $UVW1=16.40\pm0.07$, $UVM2=16.51\pm0.08$, and $UVW2=16.14\pm0.07$, somewhat larger than the ground-based measured U (309-373 nm) magnitude $m_U=15.1\pm0.2$ for the whole galaxy (Surace & Sanders 2000), but not inconsistent when accounting for the different aperture sizes used.

2.2 Chandra

Chandra observed Mrk 273 on April 2000 with the Advanced CCD Imaging Spectrometer (ACIS-S3 Back Illuminated chip) for about 47 ks. The data relevant to this observation were already published in Paper I, where they were processed following the Standard Data Processing available at that time (R4CU5UPD13.2, January 2001).

Here, data are reprocessed using much more recent versions of the Chandra Interactive Analysis of Observations software (CIAO 3.2) and Chandra Calibration Database (CALDB 3.0.0). The most effective improvements introduced concern the more accurate background subtraction and more careful computation of the effective areas and response matrices.

We start processing data from the level = 1 event file. We apply the recently released, time-dependent gain correction[*], which is necessary to adjust the "effective gains'', which have been drifting with time due to increasing charge transfer inefficiency. Since the observation was taken in the VFAINT mode, we run the tool acis_process_events to flag probable background events using all the information of the pulse heights in a $5\!\times \!5$ event island to help distinguishing between good X-ray events and bad events that are most likely associated with cosmic rays. With this procedure, the ACIS particle background can be reduced significantly compared to the standard grade selection[*]. Real X-ray photons are practically not affected by such cleaning (only about 2% of them are rejected, independently of the energy band, provided there is no pileup).

The data are filtered to include only the standard event grades 0, 2, 3, 4 and 6. We finally filter time intervals with high background by performing a 3$\sigma$ clipping of the background level using the script analyze_ltcrv[*]. This yield an effective exposure time of about 43 ks for the ACIS-S3 chip in the energy range 0.3-10 keV.

As in Paper I, we separate the spectral analysis of the nuclear emission from the extended soft X-ray halo. The spectrum of the nuclear region of the source is extracted from a circular region of radius 10'' centered at the peak of the hard X-ray emission, while, for the extended soft halo, we choose the region lying between the inner circle and the ellipse shown in Fig. 1 in order to have the best possible signal-to-noise ratio.


  \begin{figure}
\par\includegraphics[width=7.8cm,clip]{3558fig1.ps}\end{figure} Figure 1: Adaptively smoothed 0.3-10 keV Chandra image of Mrk 273 with a smoothing scale of 2 pixels ($\approx $1''). The circle has a radius of 10'' and it is centered at the peak of the hard X-ray emission. The ellipse has a size of $42'' \times 25''$ and a position angle of $13\hbox {$^\circ $ }$. The bright source visible in the upper left corner is Mrk 273x.
Open with DEXTER

We use the events included in each of the extraction regions defined above to produce a spectrum (PHA) file, which we grouped with a minimum of 20 counts per bin, unless otherwise specified. The background is obtained from a source free circular region of radius 50'' located on the same chip. The background file is then scaled to the source file by the ratio of the geometrical area. The response matrices and the ancillary response matrices of each spectrum are computed with acisspec for the same regions from which the spectra are extracted.

The X-ray spectral analysis has been performed using XSPEC version 11.3.1.

The reported errors on model parameters are at the 90% confidence level for one interesting parameter ( $\Delta \chi^2=2.71$).

Table 1: XMM- Newton and Chandra net number of counts (background subtracted) collected by each instrument in the specifiedextraction region and spectral range.

3 Spectral analysis

In Table 1 we give a comparison between the XMM-Newton and Chandra data sets. The Chandra ACIS-S spectrum of the nuclear region (inner 10'') of Mrk 273 is dominated by the background above 8 keV, while the spectrum of the extended soft X-ray halo is source dominated only between 0.5 and 0.9 keV. Both the EPIC-PN and the two MOS spectra are source dominated between 0.3 and 10 keV, therefore the entire energy range can be analyzed with relatively good signal-to-noise. XMM-Newton therefore provides some further constraints on the shape of the continuum between 8 and 10 keV.

The extracted PN and ACIS-S light curves do not show any significant (within 30%) flux or spectral variability.

Throughout our spectral analysis we fix the abundances of the elements to the solar values (Anders & Grevesse 1989) and we model Galactic absorption with the photoelectric absorption model PHABS in XSPEC, where we fixed the absorbing column density to $N{\rm _{g}}=1.09\times10^{20}$ cm-2 (HEASARC W3nH[*]). Line energies are given in the rest frame of the source.

3.1 The complex Fe K$\alpha $ line

In order to study the properties of the Fe K$\alpha $ line we concentrate our spectral analysis only on the hard X-ray (4-10 keV) spectrum of the EPIC-PN. In this energy range the continuum can not be reproduced by a simple unabsorbed power law, since this model gives a negative spectral index and a poor fit. Allowing for extra absorption at the redshift of the source (ZPHABS in XSPEC), we find a much better fit result ( $\chi^2_r\simeq1.2$), but, given the limited statistics available, if both $\Gamma $ and $N{\rm _{H}}$ are left free to vary, we find very large uncertainties (e.g. $\Gamma=2.4^{+1.2}_{-0.7}$). Fixing the spectral index to the "canonical'' value $\Gamma=1.9$ allows us to measure a column density $N{\rm _{H}}=(6.2\pm1.5)\times10^{23}$ cm-2 intrinsic to the source. As displayed in Fig. 2, the fit is still not completely satisfactory ( $\chi^2/{\rm d.o.f}=24.5/20$), mainly because of the presence of a feature in the spectrum at about 6.4 keV, which can be identified with an Fe K$\alpha $ emission line. A significant improvement in the fit is obtained when a single broad Gaussian is added to the model ( $\chi^2/{\rm d.o.f.}=15.2/17$). The line centroid is at $6.51\pm0.19$ keV and the line is found to be broad ( $\sigma=0.26^{+0.37}_{-0.17}$ keV), with an equivalent width EW=560+490-330 eV.

Such a broad iron line ( ${\it FWHM}~\sim29~000$ km s-1) could, in principle, originate from the accretion disc, since the primary radiation is obscured only up to about 4-5 keV and therefore should only partly affect the photons at the iron line energy. We try to model the iron line profile with the DISKLINE model in XSPEC for a Schwarzschild black hole, which parameterizes the radial emissivity of the disc as a power law (i.e. ${\propto} r^{-q}$). We fix the line rest energy to 6.4 keV, the emissivity of the disc q to -2, the inner and outer radii respectively to 6 and to 1000 $r{\rm _{g}}$ and the inclination angle to $30^{\circ}$. We find an equivalent width $EW=520\pm300$ eV and a fit result statistically equivalent to the one obtained with a single broad Gaussian ( $\chi^2/{\rm d.o.f}=17.1/19$). Unfortunately, the statistics are too poor to derive any physical properties of the disc.

The hypothesis of a broad iron line originating from the disc cannot be ruled out, but the high column density, which is obscuring the AGN up to about 4 keV, certainly makes the detection of the redshifted wing of a relativistic iron line profile much less feasible. A more likely explanation for the prominent iron line feature observed may be given in terms of a complex profile due to the superposition of, at least, two unresolved line components: a neutral iron line at 6.4 keV, probably originating from a Compton-thick torus, and an ionized (Fe XXV) iron line at 6.7 keV. Indeed, the latter could be produced either by some diffuse collisionally ionized plasma at temperature $kT\geq5$ keV, as observed in at least another two ULIRGs (namely NGC 6240, Boller et al. 2003; Netzer et al. 2005 and Arp 220, Iwasawa et al. 2005), or by reflection off some photoionized gas surrounding the nucleus, as often seen in Compton-thick Seyfert galaxies (e.g. Circinus, Sambruna et al. 2001; NGC 1068, Kinkhabwala et al. 2002; NGC 4945, Done et al. 2003; Mrk 3, Bianchi et al. 2005).

When the line profile is modeled with two narrow ($\sigma=0$ eV) Gaussian lines with energies fixed at 6.4 and 6.7 keV the result is statistically equivalent to a single broad line model ( $\chi^2/{\rm d.o.f.}=15.5/18$). In this case we measure an equivalent width EW=165+120-110 eV and $EW=120\pm95$ eV for the neutral and for the ionized component respectively. Finally, we find no improvement in the fit for the addition of a third ionized iron line at 6.97 keV from Fe XXVI. However upper limits on the flux ( ${<}0.18\times10^{-13}$ erg cm-2 s-1) and the equivalent width (<85 eV) of this line can be derived.

In the following we discuss two distinct models capable of reproducing the broad band spectrum as well as the possible ionized iron lines: one through thermal emission, therefore involving the presence of a hot gas associated with a starburst, and the other one through reflection of the primary continuum by some optically thin photoionized gas.


  \begin{figure}
\par\includegraphics[angle=-90,width=8.2cm,clip]{3558fig2.ps}\end{figure} Figure 2: The 4.0-10 keV EPIC-PN spectrum fitted with a simple absorbed power law model ($\Gamma $ is fixed to 1.9). Residuals are found at Fe K$\alpha $ emission line energy (at about 6.5 keV).
Open with DEXTER

3.2 The thermal emission model

At first we concentrate our analysis on the EPIC-PN spectrum. The better statistics, especially at high energies, permit more robust constraints of the model parameters. Once the best fit for the PN is found, we perform a combined fit of the PN and the two MOS, applying the same model, in order to further improve the statistics.

When the 0.3-10 keV EPIC-PN spectrum is modeled with a simple absorbed (Galactic plus intrinsic) power law, we obtain a statistically unacceptable fit ( $\chi^2_r>4$). Strong residuals remain primarily above 4 keV, but also between 0.6 and 1 keV (Fig. 3).

As a starting point for our spectral modelling we adopt the best fit absorbed power law plus Gaussian model for the 4-10 keV band (see Sect. 3.1) and extrapolate to 0.3 keV. A very prominent soft excess is clearly noticeable. Since we intend to investigate both the origin of the soft excess and of the possible ionized iron, in the following we fix the parameters of the Gaussian profile to E=6.4 keV and $\sigma=0$ eV. This will be considered as our baseline model.


  \begin{figure}
\par\includegraphics[angle=-90,width=8.38cm,clip]{3558fig3.ps}\end{figure} Figure 3: A simple absorbed power law model applied to the EPIC-PN spectrum of Mrk 273. Strong residuals (data/model) are clearly present above 4 keV and also between 0.6 and 1 keV.
Open with DEXTER

Table 2: XMM- Newton EPIC-PN spectral fitting results when the thermal emission model is applied (2 APEC: two temperature model; 3 APEC: three temperature model). The last column displays the result of the combined PN+MOS fit with the 3 APEC model. Fluxes are given in units of 10-13 erg cm-2 s-1, luminosities in units of 1043 erg s-1; both are corrected for absorption. * denotes a fixed parameter.

Diffuse, optically thin, thermal plasma emission is often observed in Starburst galaxies and ULIRGs. It is typically associated with star-forming regions and shock-heated gas in merging systems. We try to model the soft, extended, X-ray emission with a thermal spectrum from a collisionally ionized plasma, using the APEC model (Smith et al. 2001) in XSPEC. A single temperature plasma component gives a considerable improvement to the fit with respect to the baseline model ( $\chi^2/{\rm d.o.f.}=259/82$), but can not simultaneously account for the observed peak at about 0.7-0.9 keV and the possible ionized iron lines.

Therefore, we add to the baseline model two collisionally ionized plasma components of different temperatures. The free parameters are now the temperature and the normalization of each plasma component, the column density of the absorber at the redshift of the source, the normalization of the power law and the normalization of the Gaussian line. The best fit parameters for this model are given in Table 2. A plasma temperature of $\sim$0.7 keV is needed in order to model the peak in the emission at low energies, while a much hotter ($\sim$5 keV) plasma is required to fit part of the high energy continuum and also to give rise to the ionized iron line at 6.7 keV. The second, hotter, plasma is required with a high statistical significance ( $\Delta\chi^2\simeq175$ for the addition of two free parameters). However the fit is still not completely satisfactory ( $\chi^2/{\rm d.o.f.}=84.2/80$) due to the presence of some excess below about 0.6 keV.

The addition of a third, lower temperature, plasma component provides a further improvement to the fit ( $\Delta\chi^2\simeq20$ for the addition of two free parameters). In this case we find an undoubtedly acceptable result ( $\chi^2/{\rm d.o.f.}=63.8/78$). The lowest temperature component is required with a probability of more than 99.99% according to an F-test. Its temperature ( kT=0.26+0.08-0.06 keV) is consistent with the one derived from the spectrum of the extended hot gas halo ( $kT=0.38\pm0.04$ keV) observed in the Chandra image (see Sect. 4.2).

The best fit to the 0.3-10 keV PN spectrum is therefore obtained using three thermal plasma components with different temperatures, a highly absorbed power law and a neutral Fe K$\alpha $ line (Fig. 4). According to this model the emission below about 4 keV would be purely thermal. We will discuss the possible implications of these new findings in Sect. 5.


  \begin{figure}
\par\includegraphics[angle=-90,width=8.3cm,clip]{3558fi4a.ps}\par\includegraphics[angle=-90,width=8.3cm,clip]{3558fi4b.ps}\par\end{figure} Figure 4: ( Upper panel) The EPIC-PN spectrum of Mrk 273 fitted with three different temperature APEC models, a highly absorbed power law and a neutral Fe K$\alpha $ line. Each model component is displayed. ( Lower panel) The EPIC-PN and the two MOS combined fit with the same model. The closed and open circles correspond to the two MOS data sets. The best fit parameters are given in Table 2.
Open with DEXTER

Simultaneously fitting the PN and the two MOS data results in an equivalent fit, as shown in Fig. 4. All the parameters are consistent with those obtained using the PN only (see Table 2). Combining the PN and the MOS, we find a temperature kT=0.31+0.09-0.07 keV for the coolest thermal component, which is in very good agreement with the one measured by Chandra for the hot extended halo.

There is no need for further absorption at the redshift of the source for any of the three thermal components. When left free to vary the absorbing column density of each component is consistent with the Galactic value and does not improve the fit ( $\Delta\chi^2<1$ for the addition of one free parameter).

3.3 The photoionized gas model

An alternative interpretation for the origin of the soft X-ray continuum and the ionized iron line observed in the XMM-Newton spectrum is scattering of the primary radiation by some photoionized, Compton-thin, hot gas surrounding the AGN. Ionized Fe K$\alpha $ emission lines from H- and He-like iron ions are expected to be produced through recombination and resonant scattering (Bianchi & Matt 2002).

The analysis of the Chandra image in Paper I revealed that the emission below 1 keV is extended over a region of radius $\sim$7.5 kpc surrounding the nucleus and, partly, over a much larger region ($\sim$45 kpc) encompassing the tidal tail. A pure reflection model cannot account for such extended X-ray emission, which must be of thermal origin. If we model the 0.3-10 keV PN continuum with a highly absorbed (primary) plus a less absorbed (reflected) power law having the same spectral index, we find an unacceptable fit ( $\chi^2_r\simeq2.6$). Therefore we added to the model two narrow Gaussian lines, at 6.4 keV from neutral iron and at 6.7 keV from Fe XXV, plus a hot collisionally ionized plasma component. In this case the fit result is acceptable ( $\chi^2/{\rm d.o.f}=66.9/79$). However, there is no need for an additional lower temperature thermal plasma component, which is clearly detected by Chandra.

Table 3: XMM- Newton EPIC-PN spectral fitting results using the photoionized gas model described in the text. Fluxes are given in units of 10-13 erg cm-2 s-1, luminosities in units of 1043 erg s-1; both are corrected for absorption. * denotes a fixed parameter.

This apparent discrepancy between the two data sets may be reconciled simply by taking into account the absorption of the secondary power law, which is apparently negligible in the XMM-Newton spectrum, but ${\sim}7\times10^{20}$ cm-2 in the Chandra spectrum of the "nuclear'' region (see Sect. 4.1). Therefore we fix the column density of the material obscuring the reflected power law to the value measured by Chandra. The free parameters of the model are: the temperature and the normalization of the thermal component; the column density of the absorber at the redshift of the source, which is obscuring only the primary power law; the normalizations of the two power law components; and the normalizations of the two Gaussian lines.

When this model is applied to the PN spectrum the fit result is still not completely satisfactory ( $\chi^2/{\rm d.o.f}=84.9/80$), mostly due to some excess still present below 0.6 keV. The best fit parameters for this model are given in Table 3.

The addition of a second, lower temperature, plasma component results in a considerable improvement to the fit ( $\Delta\chi^2\simeq20$ for the addition of two free parameters, $\chi^2/{\rm d.o.f.}=64.0/78$). Figure 5 shows the best fit obtained for this model. The cooler thermal plasma component is required with a probability of more than 99.99% according to an F-test. The temperature found kT=0.26+0.08-0.03 keV is consistent, within the uncertainties, with that measured by Chandra for the extended hot gas halo.

The flux of the reflected power law amounts to 2-6% of the primary component. From this ratio we can estimate the column density of the reflecting material. Assuming a covering factor of 0.5, the column density of the photoionized gas would be approximately 1023 cm-2. The equivalent widths of the ionized iron lines expected to be produced by a hot gas with such a density are consistent with the measured values within the uncertainties.

The best fit model obtained here is similar to the one presented in Paper I, with the only noticeable differences being a slightly higher absorbing column density obscuring the primary radiation and a slightly lower temperature of the cooler thermal component.

4 Comparison with Chandra

For a detailed spatial analysis of the Chandra X-ray image of Mrk 273 we refer the reader to Paper I. Here, as in Paper I, we extract separately the spectra of the nuclear region and that of the extendend soft halo from the two regions shown in Fig. 1, in order to have the best possible signal-to-noise ratio.

4.1 The inner 10'' region

We apply to the Chandra spectrum of the inner region ( $10''\sim7.5$ kpc surrounding the AGN) both models successfully tested in Sect. 3 for the XMM-Newton data.

The thermal emission model, which was not tested in Paper I, when applied to the Chandra spectrum of the "nuclear'' region, also gives a good fit ( $\chi^2/{\rm d.o.f}=65.0/75$). The best fit parameters are given in Table 4. In this case, only two collisionally ionized plasma components are required, as shown in Fig. 6.


  \begin{figure}
\par\includegraphics[angle=-90,width=8.2cm,clip]{3558fig5.ps}\end{figure} Figure 5: The EPIC-PN spectrum fitted with two thermal plasma components, a highly absorbed and a less absorbed power law with the same spectral index, plus a narrow Fe K$\alpha $ line at 6.4 keV. The best fit parameters for this model are given in Table 3.
Open with DEXTER

The photoionized gas model, which we discuss here, is similar to the one presented in Paper I. It consists of a highly absorbed (primary) plus a less absorbed (secondary) power law having the same spectral index, two Fe K$\alpha $ lines at 6.4 and 6.7 keV and a single thermal plasma component. This model also provides a good fit ( $\chi^2/{\rm d.o.f}=62.0/75$, see Fig. 7). A column density of $N{\rm _{H}}=6.8^{+2.3}_{-2.8}\times10^{20}$ cm-2 obscuring the secondary power law is needed, in agreement with Paper I and the XMM-Newton results.

A narrow, unresolved, Fe K$\alpha $ line is detected at about 3$\sigma$ significance in the ACIS-S spectrum. Its centroid is found at $6.34\pm0.04$ keV, which is slightly lower than the expected value of 6.4 keV from neutral iron, probably due to uncertainties in the calibration. No ionized iron lines are required by the data; however, we derive upper limits on the fluxes of the Fe XXV and Fe XXVI lines (respectively ${<}0.37\times10^{-13}$ and ${<}0.39\times10^{-13}$ erg cm-2 s-1), which are consistent with those measured by XMM-Newton.

The best fit parameters for both models (summarized in Table 4) are consistent, within the uncertainties, with those found by XMM-Newton. The only exception is in the value of the column density obscuring the primary radiation, which is slightly different between the two observations ($41\pm6$ and $69^{+16}_{-19}\times10^{22}$ cm-2 for the Chandra and the XMM-Newton observation respectively).

Table 4: Chandra ACIS-S spectral fitting results when the thermal emission model (Model 1) and the photoionized gas model (Model 2) are applied to the spectrum of the inner 10''. The last column refers to the spectrum of the extended hot gas halo instead. Fluxes are given in units of 10-13 erg cm-2 s-1, luminosities in units of 1043 erg s-1; both are corrected for absorption. * denotes a fixed parameter.


  \begin{figure}
\par\includegraphics[angle=270,width=8.2cm,clip]{3558fig6.ps}\par\end{figure} Figure 6: The 0.3-8 keV ACIS-S spectrum of the inner 10'' region fitted with two thermal plasma components, a highly absorbed power law and a narrow Fe K$\alpha $ line at 6.4 keV. The best fit parameters are given in Table 4.
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=270,width=8.3cm,clip]{3558fig7.ps}\par\end{figure} Figure 7: The 0.3-8 keV ACIS-S spectrum of the inner 10'' region fitted with a single thermal plasma component, a highly absorbed and a less absorbed power law with the same spectral index, plus a narrow Fe K$\alpha $ line at 6.4 keV. The best fit parameters are given in Table 4.
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=270,width=8.1cm,clip]{3558fig8.ps}\end{figure} Figure 8: The 0.5-0.9 keV unbinned spectrum of the extended hot gas halo. A single APEC model with a temperature of $0.38\pm 0.04$ keV provides a satisfactory fit.
Open with DEXTER

4.2 The hot gas halo

The ACIS-S spectrum of the extended hot gas halo, which is extracted from the region defined in Sect. 2.2, is source dominated only in the 0.5-0.9 keV band. Given the limited amount of photons (315) detected in this band, we prefer to use Cash statistics on the unbinned spectrum. As displayed in Fig. 8, the spectrum is consistent with a single temperature plasma component with $kT=0.38\pm0.04$ keV. We performed a Monte Carlo calculation to test the goodness of the fit and we found that 63% of the 104 simulated spectra have a fit statistics less than that for the data. A fraction of approximately 50% is expected, if the observed spectrum was produced by the model. If the $\chi^2$ statistics is applied to the spectrum grouped to have a minimum of 20 counts per bin, fixing the model parameters to those previously found (see Table 4), we obtain $\chi^2/{\rm d.o.f}=15/11$.

In this analysis we adopted a more conservative approach in subtracting the background in comparison to Paper I. This resulted in a slightly lower temperature for the extended hot gas halo with respect to the previous analysis, but in very good agreement with the XMM-Newton results.

5 Discussion

A broad Fe K$\alpha $ emission line with $\sigma=0.26^{+0.37}_{-0.17}$ keV is detected at high significance in the broad band (0.3-10 keV) EPIC-PN spectrum ( $\Delta\chi^2=11.4$ for three additional free parameters, corresponding to more than 99% confidence level). The hypothesis of a broad iron line originating from the disc is plausible, but the high column density, obscuring the AGN up to about 4 keV, makes the detection of the redshifted wing of a relativistic iron line profile more difficult. The most likely explanation for such a broad line profile is the superposition of a narrow, unresolved, neutral iron line at 6.4 keV, probably associated with the AGN torus, and a blend of Fe XXV and Fe XXVI lines, associated with a hot plasma, in which this transitions are either thermally excited or induced by photoionization. Unfortunately, the modest statistics do not allow us to individually resolve the three possible lines. However the fluxes and equivalent widths measured by Chandra and XMM-Newton are consistent with each other.

There is no evidence for extra absorbing material obscuring any of the thermal plasma components. This could be due to a non-spherical geometry of the gas surrounding the nuclear region, given the complex morphology of this merger, or it could be simply due to the low statistics. A slightly different column density is found to be obscuring the primary radiation between the two observations. Again this could be, associated with the complex morphology of the merger, since, for instance, the projected separation between the two nuclei is $\sim$1'', or it could be an effect due to the low statistics available at high energies.

According to the thermal emission model, the temperatures of the three plasma components found in Mrk 273 are remarkably similar to those observed in other ULIRGs (e.g. NGC 6240 and Arp 220) and in the local starburst galaxy NGC 253 (Pietsch et al. 2001). The thermal emission at ${\sim}0.3$ keV, as clearly shown from the Chandra image, extends on a large scale ($\sim$45 kpc) and embraces the long tidal tail of the merger. Therefore, it must be associated with hot gas distributed in the halo. Interestingly, thermal emission at ${\sim}0.7$ keV seems to be ubiquitous in the spectra of ULIRGs, probably being associated with a nuclear or circumnuclear starburst (Franceschini et al. 2003). The presence of a high temperature ($\gtrsim$5 keV) thermal component is less frequently observed, but not so unusual, being also detected in the XMM-Newton spectra of NGC 6240 (Boller et al. 2003) and Arp 220 (Iwasawa et al. 2005). The ratio of the 2-10 keV unabsorbed X-ray luminosity of the highest temperature thermal component to the bolometric luminosity in Mrk 273 is $L_{2{-}10~{\rm keV}}/L_{\rm FIR}\simeq7\times10^{-5}$ ( $L_{\rm FIR}=L(8{-}1000~\mu{\rm m}\sim L_{\rm bol}$), a factor of 20 smaller than in NGC 6240 and a factor of 5 larger than in Arp 220.

On the other hand, the second model investigated suggests that a secondary hard X-ray component originates from the reflection of the primary continuum on some photoionized circumnuclear material. This is often observed in Compton-thick Seyfert 2 galaxies (Done et al. 2003; Bianchi et al. 2005; Sambruna et al. 2001; Kinkhabwala et al. 2002), which, in some cases, also show evidence for starburst activity (e.g. Circinus, Ruiz et al. 2001; NGC 1068, Smith & Wilson 2003). As was observed in Paper I, two bright clumps of soft X-ray emission coincide with [OIII] emission. This evidence seems to support the hypothesis of a photoionized gas surrounding the nucleus, since in several Compton-thick Seyfert galaxies, the extended soft X-ray emission expected to be produced by reflection from that photoionized gas has been observed to coincide with the [OIII] emission and the ionization cone (e.g. Sako et al. 2000; Young et al. 2001).

The absorption-corrected hard X-ray luminosity of Mrk 273 is $L_{2-10~{\rm keV}}\sim7\times10^{42}$ erg s-1, which corresponds to ${\sim}2\times10^{-3}$ of the far-IR luminosity ( $L_{\rm FIR}\sim10^{12}~L_{\odot}$, e.g. Genzel et al. 1998). Such a ratio approaches values of pure starbursts and is orders of magnitude smaller than in Compton-thin AGNs (Risaliti et al. 2000). According to the thermal emission model, the absorption-corrected 0.3-10 keV luminosity due to the three thermal components would be approximately $0.7\times10^{42}$ erg s-1. Comparable soft X-ray luminosities have been measured for NGC 6240 (Netzer et al. 2005) and other ULIRGs in which the starburst contribution to the bolometric luminosity is dominant (Franceschini et al. 2003). From the photoionized gas model the resulting thermal contribution to the X-ray luminosity is about a factor of 3 smaller. For Mrk 273 the unobscured AGN X-ray luminosity is a more modest fraction of the bolometric luminosity than in typical AGN. This suggests additional strong star formation which will naturally explain the thermal X-ray components invoked in one of our two scenarios. The presence of strong star formation agrees with the detection from mid-infrared spectroscopy (Genzel et al. 1998) of both very powerful star formation and an AGN in Mrk 273.

Besides the two models discussed above, other interpretations are possible. X-ray emission from luminous radio supernovae like SN1995J (Fox et al. 2000), or the integrated contribution from low- and high-mass X-ray binaries, as observed from the spectra of local starburst galaxies (Persic & Rephaeli 2002), may both give rise to a very hot ($kT\sim5$ keV) thermal component, or a $\Gamma\sim2$ power law. Following Kennicutt (1998), the star formation rate estimated assuming that the entire $L_{\rm FIR}$ is due to starburst would be $\sim$200 $M_{\odot}$ yr-1. Therefore, the inferred supernova rate is $\sim$2 yr-1, which is in agreement with the value estimated from radio observations Bondi et al. (2005). On the other hand, given a 200 $M_{\odot}$ yr-1 star formation rate, the estimated emission from X-ray binaries, following Franceschini et al. (2003), would be $L_{2-10~{\rm keV}}{\sim}2\times10^{41}$ erg s-1, which is very close to the observed value ( ${\sim}3\times10^{41}$ erg s-1). X-ray spectra from supernovae have shown evidence for a strong Fe K line at 6.7 keV (e.g. SN1986J, Houck et al. 1998), therefore could also account for this feature, while emission from X-ray binaries cannot.

6 Conclusions

We analyzed the X-ray spectral properties of Mrk 273 combining the high throughput of XMM-Newton with the high spatial resolution of Chandra.

From the XMM-Newton spectrum we found a broad Fe K$\alpha $ emission line with a high statistical significance ($\gtrsim$99% c.l.). We suggested the superposition of multiple unresolved iron line features: one from neutral iron at 6.4 keV ($EW\sim170$ eV) and one from a blend of Fe XXV at 6.7 keV ($EW\sim120$ eV) and Fe XXVI at 6.97 keV (EW<85 eV).

We tested two different models capable of reproducing the possible ionized iron lines detected: one through thermal emission, therefore involving the presence of a hot gas associated with a starburst, and the other one through reflection of the primary continuum by optically thin photoionized gas. Unfortunately, given the available statistics, we cannot rule out either of the two possible explanations. Some indications that a reflection componet is indeed present in the nuclear region of Mrk 273 may come from polarization measures, which are already available in the IR (Siebenmorgen & Efstathiou 2001), but still not in the X-rays. High resolution spectroscopy could be helpful to detect some of the forbidden lines or the Radiative Recombination Continua (RRC), which are considered as typical signatures of photoionization (Porquet & Dubau 2000; Liedahl 1999; Liedahl & Paerels 1996).

Acknowledgements
We thank Günther Hasinger and Andrea Comastri for useful discussions, Stefano Bianchi and Paolo Tozzi for their help on the reduction and spectral analysis of the XMM-Newton and Chandra data and the anonymous referee for valuable comments and suggestions. This paper is based on observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly founded by ESA Member States and NASA.

   
Appendix A: The unabsorbed Seyfert 2 Mrk 273x

Over recent years, much attention has been directed to the study of unabsorbed Seyfert 2 galaxies (e.g. Ptak et al. 1996; Panessa & Bassani 2002; Gallo et al. 2005). These galaxies are of interest for two reason: they may account for as much as 30% of the Seyfert 2 population (Panessa & Bassani 2002), and they are a challenge to explain with the existing Seyfert unification model.

Mrk 273x (z = 0.458) is one of the better known unabsorbed Seyfert 2 galaxies, perhaps because it falls in the field-of-view of most observations of Mrk 273. In this appendix, we discuss the properties of Mrk 273x as seen during the XMM-Newton observation of Mrk 273.

In the EPIC-PN image, Mrk 273x falls directly in a chip gap; thus we used the combined MOS 1/2 data for spectral analysis. The combined spectrum is source dominated in the 0.3-10 keV band ( 0.44-14.6 keV in the rest-frame). The total number of source plus background counts is 526, whereas the total number of background counts is 65 (scaled to source cell size).

The best-fit continuum model of Paper I to the Chandra data of Mrk 273x was a power law plus a weak thermal component. An emission line at 1.82 keV, which the authors suggested could be due to Si XVI, was significantly detected above the continuum. The XMM-Newton data are consistent with a power law ( $\Gamma = 1.49 \pm 0.12$) modified by Galactic absorption ( $\chi^2/{\rm d.o.f.}=60.2/82$; Fig. A.1).


  \begin{figure}
\par\includegraphics[angle=270,width=8.1cm,clip]{3558fig9.ps}\end{figure} Figure A.1: The combined MOS 1 and MOS2 0.3-10 keV data fitted with a power law modified by Galactic absorption.

The 90% confidence upper limit on the intrinsic absorption is ${<}4.5\times10^{20}$ cm-2, at least a factor of three less than what was required in the Chandra model.

Multi-component continuum models (e.g. blackbody plus power law, double power law, or a broken power law) also fit the data well, but they are not significant improvements over the single power law. The addition of a $\sim$1.82 keV emission feature is not required by the data ( $\Delta\chi^2=1$ for two additional free parameters; ${\it EW}<33$ eV). Positive residuals at about 4.3 keV can be modelled with a narrow Gaussian profile with a rest-frame energy of E=6.29+0.14-0.34 keV. The addition of the Gaussian profile is an improvement to the residuals, but it is not statistically significant ( $\Delta\chi^2=3$ for three additional free parameters).

Based on the single power law fit, the unabsorbed 0.5-10 keV flux is ( $1.43^{+0.22}_{-0.18})\times10^{-13}$ erg s-1 cm-2( $1.18\times10^{-13}$ erg s-1 cm-2 in the 2-10 keV band). The intrinsic 0.3-10 keV unabsorbed luminosity is $9.2\times10^{43}$ erg s-1. The broad band flux measured during this XMM-Newton observation is comparable to the flux observed during the Chandra observation two years earlier. The $\sim$22 ks light curve of Mrk 273x shows no variability and is perfectly fitted by a constant.

References

 

Copyright ESO 2005