A&A 420, 853-872 (2004)
DOI: 10.1051/0004-6361:20035687

The XMM-NEWTON $\Omega $ project

I. The X-ray luminosity-temperature relation at ${\textit{z}>0.4}$[*]

D. H. Lumb1 - J. G. Bartlett2,3 - A. K. Romer4,9 - A. Blanchard5 - D. J. Burke6 - C. A. Collins7 - R. C. Nichol4 - M. Giard8 - P. B. Marty5 - J. Nevalainen6 - R. Sadat5 - S. C. Vauclair5

1 - Science Payloads and Advanced Concepts Office, European Space Agency, ESTEC, 2200AG Noordwijk, The Netherlands
2 - APC - Université Paris 7/PCC - Collège de France, 11 Pl. Marcelin Berthelot, 75231 Paris Cedex 05, France
3 - Centre de Données astronomiques de Strasbourg, 11 rue de l'Université, 67000 Strasbourg, France
4 - Physics Department, Carnegie Mellon University, Pittsburgh, PA 15213, USA
5 - Laboratoire d'Astrophysique, OMP, CNRS, UPS, 14 Av. Ed. Belin, 31400 Toulouse, France
6 - Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA
7 - Astrophysics Research Institute, Liverpool John Moores University, Twelve Quays House, Egerton Wharf, Birkenhead CH41 1LD, UK
8 - Centre d'Étude Spatiale des Rayonnements, 9 avenue du Colonel Roche, BP 4346, 31028 Toulouse, France
9 - Astronomy Centre, Dept. of Physics and Astronomy, University of Sussex, Brighton, UK

Received 14 November 2003 / Accepted 17 March 2004

We describe XMM-Newton Guaranteed Time observations of a sample of eight high redshift (0.45<z<0.62) clusters. The goal of these observations was to measure the luminosity and the temperature of the clusters to a precision of $\sim $10%, leading to constraints on the possible evolution of the luminosity-temperature ( Lx-Tx) relation, and ultimately on the values of the matter density, $\Omega_{\rm M}$, and, to a lesser extent, the cosmological constant $\Omega_{\Lambda}$. The clusters were drawn from the SHARC and 160 Square Degree (160SD) ROSAT surveys and span a bolometric (0.0-20 keV) luminosity range of 2.0 to $14.4 \times 10^{44}$ erg s-1 (H0=50, $\Omega_{\rm M}=1$, $\Omega_{\Lambda}=0$). Here we describe our data analysis techniques and present, for the first time with XMM-Newton, a Lx-Tx relation. For each of the eight clusters in the sample, we have measured total ( $r<r_{\rm virial}$) bolometric luminosities, performed $\beta $-model fits to the radial surface profiles and made spectral fits to a single temperature isothermal model. We describe data analysis techniques that pay particular attention to background mitigation. We have also estimated temperatures and luminosities for two known clusters (Abell 2246 and RX J1325.0-3814), and one new high redshift cluster candidate (XMMU J084701.8+345117), that were detected off-axis. Characterizing the Lx-Txrelation as $L_{x} = L_{6} (\frac{T}{6~\rm keV})^{\alpha}$, we find $L_{6}=15.9 ^{+7.6}_{-5.2} \times
10^{44}$ erg s-1  and $\alpha=2.7\pm0.4$ for an $\Omega_{\Lambda}=0.0, \Omega_{\rm M}=1.0$, H0=50 km s-1 Mpc-1 cosmology at a typical redshift $z \sim 0.55$. Comparing with the low redshift study by Markevitch (1998), we find $\alpha$ to be in agreement, and assuming Lx-Txto evolve as (1+z)A , we find $A=0.68\pm0.26$for the same cosmology and A=1.52+0.26-0.27 for an $\Omega _{\Lambda }=0.7, \Omega _{\rm M}=0.3$ cosmology. Our A values are very similar to those found previously by Vikhlinin et al. (2002) using a compilation of Chandra observations of 0.39<z<1.26 clusters. We conclude that there is now evidence from both XMM-Newton and Chandra for an evolutionary trend in the Lx-Tx relation. This evolution is significantly below the level expected from the predictions of the self-similar model for an $\Omega_{\Lambda}=0.0, \Omega_{\rm M}=1.0$, cosmology, but consistent with self-similar model in an $\Omega _{\Lambda }=0.7, \Omega _{\rm M}=0.3$ cosmology. Our observations lend support to the robustness and completeness of the SHARC and 160SD surveys.

Key words: X-rays: galaxies: clusters

1 Introduction

1.1 Motivation

As the most massive gravitationally bound objects in the universe, galaxy clusters are particularly sensitive to the evolution of the density perturbations responsible for their formation. Cluster abundance as a function of mass and redshift is dictated by the mass function (Press & Schechter 1974; Jenkins et al. 2001), which gives the comoving space density of collapsed objects as a function of mass and redshift. In standard models, the mass function falls off as a Gaussian at the high mass end, reflecting the Gaussian nature of the density perturbations. Galaxy cluster abundance is therefore exquistely sensitive to the amplitude of the perturbations and its evolution. Since this evolution is controlled by the underlying cosmological background, observations of cluster abundance offer an effective way to constrain certain cosmological parameters, such as the density parameter (Oukbir & Blanchard 1992; Blanchard & Bartlett 1998) or the dark energy parameter and equation-of-state (Majumdar & Mohr 2003; Wang & Steinhardt 1998). The present day cluster abundance is degenerate in the matter density and amplitude of the power spectrum; evolution of the abundance breaks this degeneracy. Constraints obtained in this manner are complementary to others that essentially rely on determinations of cosmological distances; for example, those found by observations of supernovae type Ia or by measurements of cosmic microwave background anisotropies.

Cluster mass is, however, difficult to measure directly, and in practice one seeks a direct observable that is closely related to virial mass. Lensing surveys would seem the most suited to the task, as the effects of lensing are of course directly related to mass (albeit projected along the line-of-sight). However, cluster mass estimates from weak and strong lensing remain controversial as they suffer from several systematic uncertainties (e.g. projection effects). Among X-ray observables, intracluster gas temperature is expected to be tightly correlated with virial mass, an expectation borne out by simple hydrostatic considerations as well as by numerical simulations (Evrard et al. 1996; Bryan & Norman 1998). X-ray luminosity, on the other hand, is a much less robust mass indicator, despite being significantly easier to measure, because it depends on the density profile of the intracluster gas, the physics of which is currently difficult to model.

\par\includegraphics[width=8.2cm,clip]{0687fig1.ps}\end{figure} Figure 1: The predicted X-ray temperature function at redshifts z=0.05 (black), 0.33 (light grey) and 0.5 (dark grey) for a flat high density model (left; $\Omega _{\Lambda }=0.12$) and a concordance model (right; $\Omega _{\Lambda }=0.7$), both fitted to the local ( $z \sim 0.05$) temperature function.
Open with DEXTER

With a calibrated Tx-M (temperature-mass) relation, the mass function can be translated into an observable cluster temperature distribution function ${\rm d}n/{\rm d}T$. The exact T-M relation to use is of course a key ingredient, one that may be addressed, for example, using numerical simulations, or directly from detailed observations that determine both cluster mass and temperature (e.g. Nevalainen et al. 2000). Figure 1 compares temperature function predictions for a high matter density and a flat model, both normalized to the present-day, observed ${\rm d}n/{\rm d}T$ (e.g. Henry & Arnaud 1991; Edge et al. 1990; Blanchard et al. 2000; Ikebe et al. 2002). Evolution toward higher redshift is strikingly different in two cosmological models differing by their matter content alone and illustrating the above argument.

Despite this obvious promise, it has been difficult to use ${\rm d}n/{\rm d}T(z)$to explore the cosmological model because direct measurements of Tx at z>0 require long integration with satellite observations. Prior to the launch of Chandra and XMM-Newton, the number of measured cluster temperatures at high redshifts was very small. This is attested to by the fact that the most distant determination of ${\rm d}n/{\rm d}T$ (Henry 2000) to date was based on a sample of only fourteen z>0.3 clusters (mean $\sim $0.38) from the Einstein Medium Sensitivity Survey (EMSS; Gioia et al. 1990).

An alternative, but related approach is to apply a luminosity-temperature ( Lx-Tx) relation to a flux-limited sample, thereby obtaining either the temperature function, or a redshift distribution at given temperature (Oukbir & Blanchard 1997; Sadat et al. 1998; Reichart et al. 1999; Borgani et al. 1999; Borgani et al. 2001). The advantage of this approach is that it does not require detailed X-ray spectroscopy of all the clusters in a flux-limited sample to establish the Lx-Tx relation and its evolution. There are now several flux limited catalogs of medium to high redshift clusters, in addition to the EMSS, to which one can apply this technique. These include those based on data from the ROSAT All Sky Survey (e.g. Ebeling et al. 2001; Zhang et al. 2003), NEP (e.g. Gioia et al. 2001; Henry et al. 2001); and from the ROSAT pointing archive (e.g. Romer et al. 2000; Burke et al. 2003; Vikhlinin et al. 1998; Mullis et al. 2003; Scharf et al. 1997; Rosati et al. 1995).

1.2 The XMM-Newton $\Omega $-project

The goal of the XMM-Newton $\Omega $-project (Bartlett et al. 2001) is to increase the number of high quality X-ray cluster temperature measurements at z>0.3, and thus enhance the scientific yield from the various ROSAT surveys. We describe below XMM-Newton observations of eight ROSAT clusters at z>0.4 (median of z=0.54) performed as part of the Guaranteed Time programme. Related observations conducted during the open time programmes of seven 0.3<z<0.4 clusters are discussed elsewhere (Majerowicz et al. 2002b; Majerowicz et al. 2004). Cosmological interpretation is presented in Vauclair et al. (2003). In Sects. 2 and 3 we describe the observations and data analysis techniques. In Sect. 4 we discuss each of the eight clusters in turn. In Sect. 5 we present our Lx-Tx relation and compare it to previous work.

Table 1: Summary of pointing directions for each observation, the overall scheduled durations, and redshift data (s Southern SHARC, Burke et al. 2003; b Bright SHARC, Romer et al. 2000; v 160SD Vikhlinin et al. 1998).

This paper serves also to introduce the Project, and represents an opportunity to provide detailed descriptions of analysis techniques used in the Project and which can be be used generally for XMM-Newton observations of cluster targets. Except where explicitly stated elsewhere we use a $\Omega_{\rm M}=1$, $\Omega_{\Lambda}=0$, q0=0.5 model with H0=50 km s-1 Mpc-1 that has been most frequently used in the past as the parameter set for X-ray cluster studies. When necessary for examining cosmological implications we correct our results to a concordance model.

2 Observation programme

2.1 The sample

Table 1 summarizes the locations, date of observations and other details of the eight clusters in this program. Seven of the eight clusters were drawn from the SHARC Surveys, four from the Southern SHARC (Burke et al. 2003) sample and three from the Bright SHARC (Romer et al. 2000) sample. The SHARC cluster samples are based on searches for clusters serendipitously detected in ROSAT PSPC observations. They complement each other in that they cover, respectively, 17.7 degrees2 to a flux limit of ${\simeq} 3.9\times 10^{-14}$ erg s-1 cm-2(Collins et al. 1997) and 178.6 degrees2 to a flux limit of ${\simeq}
1.4\times 10^{-13}$ erg s-1 cm-2 (Nichol et al. 1999). This strategy has yielded a combined cluster catalogue that straddles L*over the redshift range 0.2 < z < 0.8 and shows a consistent picture of a non-evolving cluster luminosity function (Collins et al. 1997; Burke et al. 1997) except, possibly, at luminosities greater than $L_{x}=5\times 10^{44}$ erg s-1 (Nichol et al. 1999; Adami et al. 2000). Two aspects of the SHARC surveys, makes them particularly well suited to the XMM-Newton $\Omega $ project; they have been subjected both to detailed optical follow-up (Romer et al. 2000; Burke et al. 2003) and to extensive selection function simulations (Adami et al. 2000; Burke et al. 2003). The eighth cluster in our observation programme was taken from the 160 Square Degree ROSAT Survey (160SD hereafter, Vikhlinin et al. 1998; Mullis et al. 2003). There is considerable overlap, both in terms of methodology and cluster members, between the 160SD and the SHARC surveys (3 of the 7 SHARC clusters are also members of the 160SD catalogue), however, this particular cluster (RX J0847.2) was not a member of either SHARC sample because its host PSPC observation did not meet the SHARC exposure time criterion. The selection of clusters sampled in this programme was driven solely by a requirement to observe all $z\geq0.5$ targets from the SHARC surveys. Following the visibility and observability constraints of the XMM-Newton Guaranteed Time programme we added two clusters with slightly lower redshift from the SHARC sample in addition to the one from the 160SD catalogue.

2.2 XMM-Newton

XMM-Newton (Jansen et al. 2001) comprises 3 co-aligned telescopes, each with effective area at 1.5 keV of $\sim $1500 cm2, and Full Width Half Maximum (FWHM) angular resolution of $\sim $5 arcsec. This combination of the highest ever focused X-ray collection area, and the ability to resolve clusters at all redshifts, makes XMM-Newton the best suited observatory for this programme. The 3 telescopes each have a focal plane CCD imaging spectrometer camera provided by the EPIC consortium. Two also have a reflection grating array, which splits off half the light, to provide simultaneous high resolution dispersive spectra. These two telescopes are equipped with EPIC MOS cameras (Turner et al. 2001), which are conventional CMOS CCD-based imagers enhanced for X-ray sensitivity. The third employs the EPIC PN camera (Strüder et al. 2001) which is based on a pn-junction multi-linear readout CCD. The EPIC cameras offer a field of view (FOV) of $\sim $30 arc-min diameter, and an energy resolution of typically 100 eV (FWHM) in the range $\sim $0.2-10 keV. The 7 CCDs in each of the MOS cameras are about 10 arc-min square each. The central chip encompasses the whole of our z>0.4 clusters. The 12 CCDs in the PN are about $4 \times 13$ arc-min. Even if correctly centered on the boresight PN CCD, some portion of the target cluster emission may spill onto neighbouring CCDs, and across dead zones between CCDs. In each camera, an aluminized optical blocking filter was deployed, the thickness of which was chosen to suit the expected brightness of nearby serendipitous objects in each field. Seven of the eight clusters were observed in the Full Frame Imaging mode, as appropriate for weak extended targets, but one (RX J1325.5) was observed when the MOS cameras were in "Window'' mode, due to the presence of an unusually bright nearby point source.

The data were processed through the XMM-Newton Science Analysis Sub-System (SAS; Watson et al. 2001) version 5.3, in order to register photons from detector to sky co-ordinates, to correct energy data for gain and charge transfer losses, and remove instrument noise artifacts. This provided calibrated event lists as a starting point for our detailed data reduction.

3 Data reduction

3.1 Rate filtering

Part of the XMM-Newton orbit lies within the magnetosphere, and consequently the spacecraft can encounter clouds of protons accelerated by magnetic reconnection. When these particles scatter through the mirror system they are concentrated onto the focal plane, and an enhanced background rate can occur. These intervals were identified by forming histograms of events with energy $\geq$10 keV, located in single pixels, in time bins of 50 (100) s in the PN (MOS) camera(s). Due to variations in the baseline raw cosmic ray rate experienced through the mission to date (probably due to modulation by the solar activity and/or the seasonal variation of satellite apogee direction), we prefer not filter at fixed background count rates, rather for each exposure we define ${\pm } 3 \sigma $ limits after ignoring the highest count rate periods (Figs. 2 and 3). Table 2 summarizes the total on-axis Good Time Interval (GTI) exposure times extracted for each observation using this procedure.

\par\includegraphics[width=8.1cm,clip]{0687fig2.ps}\end{figure} Figure 2: Count rate in PN ( top) and MOS ( bottom) cameras in the RX J1120.1 observation after selection for single pixel events above 10 keV. Typical rates $\sim $15 cts/50 s bin (PN) and $\sim $8 cts/100 s bin (MOS)
Open with DEXTER

\par {\includegraphics[width=7.3cm,clip]{0687fig3.ps} }
\par\end{figure} Figure 3: Histogram of the count rate sequence for PN ( top) and MOS ( bottom) cameras in the RX J1120.1 observation after selection for single pixel events above 10 keV. The dashed lines indicate the upper and lower ( ${\pm } 3 \sigma $) bounds used for the Good Time Interval filtering.
Open with DEXTER

Table 2: Summary of the usable exposure durations obtained for the different clusters after accounting for the GTI filtering, telemetry losses etc. In some observations data from only one camera was available. For operational reasons the PN camera commences exposures somewhat later than the MOS cameras. Losses due to proton flares, telemetry drops etc. are given in Col. 4 as a percentage of the total exposure time.

For clusters RX J1354.2-0222 and RX J1701.3+6414 the high prevalence of proton flares mitigated against the automated filtering procedure, and a manual analysis of the the "curve of growth'' for signal:noise in the cluster target region was employed.

3.2 Vignetting

The clusters in our sample are expected to be extended on a scale of 2-3 arc-min, and in order to assess accurately the surface brightness distribution, a small correction for energy dependent telescope vignetting must be applied. This was performed by the SAS task EVIGWEIGHT which assigns a weighting value to each detected photon, accounting for telescope and CCD efficiency variations (the latter being negligible by comparison). Thereafter, spectrum extraction and image product generation is automatically weighted for the vignetting correction, and response matrices for the on-axis location can be used[*].

An exposure map was also generated for each camera observation. These maps accounted for spatial variations, such as dead pixels, noisy readout columns and chip gaps, not included in the vignetting correction. The XMM-Newton satellite is usually very stable during observations, so the inter-CCD gaps remain approximately fixed in sky projection. It is therefore important to correct for flux lost from extended cluster regions using the exposure maps. This is especially true for the PN images, because the PN CCDs are smaller and the gaps closer to the target than those in the MOS cameras.

3.3 Background subtraction techniques

Background subtraction is a very important step in our analysis. Even after the proton flare removal described above (Sect. 3.1), and after masking out point sources, our cluster observations will be contaminated, to varying degrees, by three different background signals; (i) the cosmic X-ray background, (ii) a particle background induced by incident cosmic rays and (iii) some residual soft proton contamination. We will refer to these as the cosmic, particle and proton backgrounds respectively, hereafter. There are two approaches to background subtraction, either one can use the source observation itself (this is known as in-field subtraction), or one can use background template files. The templates are generated by combining several deep, blank field, observations and are especially useful when the source of interest covers a large fraction of the field of view (e.g. a nearby cluster). We have used a combination of both techniques in the analysis below.

For each of the eight clusters in our study, we applied the same rate filtering, vignetting correction and detector-to-sky conversion used for the reduction of the cluster observation, to the background template files provided by the XMM-SOC (Lumb et al. 2002). This allowed us to extract identical physical detector regions for both the cluster and the background using the SAS task ATTCALC. We then re-normalized the background template to account for any differences in the particle background count rates between the cluster and background template observations. This step was necessary because the particle background varies with observation epoch (Sect. 3.3.1). To make this correction, we took advantage of the fact that high energy particles are not focused by the telescope optics, so the particle background count rate can be measured from areas outside the telescope FOV.

We used re-normalized background templates (energy range 0.3-4.5 keV) during the spatial analysis (Sect. 3.5). We also used background templates (energy range 0.3-10 keV) during the spectral analysis of two of our two brightest clusters (RX J1120.1 and RX J1334.3). For this, we employed the so-called "double subtraction'' technique (Arnaud et al. 2001), which involves making an additional correction to compensate for the fact that the Galactic Halo and the Local Hot Bubble component of the diffuse cosmic X-ray background varies significantly across the sky. To determine this correction, we compared off-axis, source free, regions in both the cluster fields and their corresponding background templates. We found the correction to be small, due to the lack of soft X-ray emission features at the high Galactic latitude locations of RX J1120.1 and RX J1334.3. We have used the results from the double subtraction analysis of RX J1120.1 and RX J1334.3 to validate the results from the in-field background subtraction, see below.

For a variety of reasons, e.g. the choice of filter[*], low signal to noise, non-zero off-axis angle etc., we were not able to apply the double subtraction technique to all our clusters. Instead we have to rely on the in-field technique. The radial dependence of the vignetting means that this technique is most successful for point sources, however it should still work well for our clusters since, at z>0.4, they are barely more extended than the instrument point spread function (PSF). We illustrate this by considering a cluster described by a $\beta $-model with $\beta \sim0.67$ and a core radius of 25'' ( ${\sim} 185~ h^{-1}_{50}$ kpc typical values for the clusters in our sample, Table 5). At a radius enclosing 75% of the counts from such a cluster, the vignetting is different by less than 1% from the aim-point value. Therefore we posit that an accurate representation of the cosmic X-ray background spectrum at the cluster position can be derived using a nearby in-field background aperture. For the spectral analysis of our clusters we generally used an annular background aperture centered on the source. The annuli were chosen on a case by case basis, but were typically 2-3' wide with an inner radius no less than 3' from the cluster center (to ensure that no flux from outlying regions of the cluster was erroneously removed). We note that in the case of RX J1325.5, we could not use an annulus, due to the proximity of a very bright point source (Sect. 4.6). In Sects. 4.1 and 4.2, we compare the results of the spectral fitting for the clusters where both in field and double subtraction could be applied. As we did not find significant differences, we conclude that the in field background subtraction technique should be valid for our sample.

\par {\includegraphics[width=8.5cm,clip]{0687fig4.ps} }
\par \end{figure} Figure 4: Comparison of the relative scaling of particle (triangle) and proton (square) background rates ( after GTI filtering) in the eight cluster observations, with respect to the components of the background template files. Data averaged for all cameras in observation.
Open with DEXTER

3.3.1 Impact of residual proton contamination

The ratio of the proton to particle backgrounds is not constant; even after proton flare mitigation (Sect. 3.1), the proton background may have a significantly different counting rate in the cluster observation compared to background template files (Markevitch 2002). This is illustrated in Fig. 4. The triangles indicate the ratio of the particle background in the cluster observation to that in the background template file (Col. 2 in Table 3). The squares show the equivalent information for the proton background (Col. 3 in Table 3). For this comparison, the particle background was estimated from the $\geq$10 keV count rate in CCD areas outside the telescope FOV, whereas the proton background was estimated by subtracting that value from the $\geq$10 keV count rate within the FOV. From the Fig. 4, it is clear that the proton background varies with epoch and is anti-correlated with the particle background. This anti-correlation can be explained if

Although this anti-correlation tends to maintain a more constant total background rate, the result of this anti-correlation is that we cannot be certain the background template re-normalizations described above, compensate correctly for the temporal variations in the proton background. This suggests that the in-field background subtraction method may be preferable to double subtraction technique for spectral analysis. However, the in-field technique is complicated by the fact that the scattering of low energy protons at the mirrors occurs over angles somewhat larger than for the X-ray reflection (Rasmussen et al. 1999). Therefore, when using task EVIGWEIGHT (Sect. 3.2), there is the potential to over-weight the proton background. The proton background spectrum is hard, so this over-weighting could result in an artificial softening of the cluster spectrum and hence a lowering the estimated cluster temperature, which might in turn mimic evolution in the measured Lx-Tx relation (see Sect. 5.1).

Table 3: The scaling ratio of the >10 keV count rates outside (Col. 2) and inside (Col. 3) the telescope FOV in the background template files compared to the cluster files. The former ratio is an indication of the particle (cosmic ray) background rate, while the latter indicates the level of the soft proton background rate. These ratios are plotted as a function of Julian date in Fig. 4.

To examine the likely impact of systematic errors in the proton vignetting correction, we used the EVIGWEIGHT vignetting corrected background template files to create an image that was essentially free from proton contamination. We did so by applying very strict count-rate filtering. We then made a comparison image by applying less conservative count-rate filtering to the background template files. For this we used the filtering criteria derived from the RX J1701.3 observation (which suffered from unusually high proton contamination); the comparison is presented in Fig. 5. The two count rate images we derived, should be identical within the noise, except with regard to the proton background. By dividing one by the other, and then fitting a radial surface brightness profile, we should therefore be able to get an impression of how the proton background is vignetted. The results of this test (in the 0.5-7.0 keV band to emphasize an effect of hard proton spectrum) are shown in Fig. 5. We also generated the equivalent figure using the RX J1354.2 filtering criteria with almost identical results. From this investigation, we conclude that the proton background is over-weighted by only a few percent by EVIGWEIGHT even at the edges of the field of view in observations with significant proton contamination. We therefore chose to ignore the over-weighting of the proton background in our analysis, except in the cases of RX J1701.3 and RX J1354.2. For these, we applied a small ($\sim $2%) scaling to the background spectrum derived from an in-field annulus around the cluster[*].

\par {\includegraphics[width=8cm,clip]{0687fig5.ps} }
\end{figure} Figure 5: Ratio of the surface brightness distributions in blank fields, showing proton contamination in RX J1701.3+6414 field.
Open with DEXTER

3.4 Image products

Images for each camera were compiled in the 0.3-4.5 keV energy band, this maximizes source count rates for cluster spectra that are characterized by a $\sim $4 keV temperature, and furthermore provides a guard against corrections at the softest and hardest energy bands that are subject to the largest potential Galactic and cosmic ray subtraction inaccuracy. A spatial binning of 4.3 arc-sec per pixel, was employed, which slightly oversamples the mirror FWHM. The 3 separate count-rate images from the EPIC cameras were exposure corrected and co-added. In each image field, we use the SAS task EBOXDETECT to identify point sources via a sliding box detection algorithm (see Table A.1). In Figs. 6 and 7, we show the inner portion of the eight cluster images after a Gaussian smoothing of ${\sim} 4''$ (using the FTOOL  fGauss). It is clear from these images that the clusters are well resolved and encompass a variety of morphologies. Point sources masked during the spatial and spectral analysis (Sect. 4) are indicated by small circles. The larger circles denote the apertures used to generate the cluster spectrum.
\par {\includegraphics[width=17.5cm]{0687f6.eps} }\end{figure} Figure 6: a) Central ${\sim } 18 \times 18$ arc-min of the RX J1120.1 field - smoothed by a $4'' \sigma $ Gaussian. The large circle represents the spectral extraction radius, and the smaller circles locate the point sources which were excised. Similar images for RX J1334.3 b), RX J0337.7 c) and RX J0505.3  d). The RX J1334.3 image $15\times 15$ arc-min; the RX J0337.7 and RX J0505.3 images are $14\times 14$ arc-min.
Open with DEXTER

\par {\includegraphics[width=17.5cm,clip]{0687fig7.ps} }\end{figure} Figure 7: a) Central MOS CCD image (10 arc-min) of the RX J0847.2 field - smoothed by a 4''Gaussian. The large circle represents the spectral extraction radius, and the smaller circles locate the point sources which were excised. Similar images for RX J1325.5 b), RX J1701.3  c) and RX J1354.2 d). A newly discovered candidate distant cluster (XMMU J084701.8+345117) is labeled to the upper right of RX J0847.2. The RX J1325.5 and RX J1325.0 Southern SHARC clusters fall in an observation centered on IRAS 13224-3809. The long steak connected to IRAS 13224-3809 comes from the "out-of-time" events collected during readout of the PN camera. In the RX J1701.3 field, the nearby Abell 2246 and QSO objects are labeled. The RX J1325.5 image is ${\sim } 18 \times 18$ arcmin; the RX J1701.3 and RX J1354.2 images are ${\sim } 14\times 14$ arcmin.
Open with DEXTER

3.5 Radial surface brightness profiles

Before generating the profile, we defined the centroid of the cluster brightness distribution using a 2-d Gaussian fit around the core of the raw cluster image (0.3-4.5 keV). These centroids are given in Table 5, row 1. Next, a background correction was applied by subtracting, pixel by pixel, the corresponding re-normalized background template (Sect. 3.3). The radial bins were chosen so that the background-subtracted counts per bin, in the co-added profile, were at ${\geq} 3\sigma$ significance. A $\beta $-model convolved with the XMM-Newton telescope PSF appropriate at the position of the cluster centroid was then fit to these profiles, using simple $\chi^{2}$-minimization (background was fixed). For the PSF convolution, we used calibration file XRTn_XPSF_0004.CCF in SAS medium accuracy mode, available from the XMM-SOC calibration ftp site. The PSF was constructed by co-adding the contribution from PSFs at different energies using a weighting scheme appropriate for a 4 keV thermal spectrum. Fortuitously the XMM-Newton telescope PSF is a rather weak function of X-ray energy, so that any deviations from this default spectrum would have a negligible impact on our fits. We were able to use the on-axis PSF for the convolution, except in the case of RX J1325.5. For this cluster, which was observed off-axis, the appropriate off-axis PSF model was generated. We note that the PSF correction was applied separately for each camera. Fitted values for $\beta $ and the core radius $r\rm _c$ are given in Table 5 rows 8 and 9. To convert the core radii from angles to distances for this table, we assumed a spatially flat cosmology with $\Omega_{\rm M}=1$ and H0 = 50 km s-1 Mpc-1. The eight radial profiles and their respective best fit PSF convolved $\beta $ models are presented in Figs. 19 and 20.

3.6 Spectral extraction

Cluster spectra were extracted from circular regions centred on the cluster centroids, with radii ranging from 90 to 145 arcsec (typically 120'', see Table 5, row 13). From these regions we extracted data from the event list using the recommended pixel selections. For the PN camera we used single and double pixel events (event patterns 0-4) with all the recommended selection flags applied. For the MOS cameras, we used the 0-12 event patterns. We account for background contamination in the spectra using either the double and/or in-field subtraction methods (Sect. 3.3).

Spectral fitting was implemented with the XSPEC v11 package (Arnaud 1996), using the MEKAL models (Mewe et al. 1986) for a thermal spectrum, modified with interstellar absorption (McCammon & Sanders 1990) appropriate for the Galactic column density (Dickey & Lockman 1990). To facilitate fitting via Chi-square minimization, the spectral files were re-binned to ensure at least 25 counts per bin in order to approximate Gaussian statistics. The energy range was 0.3-10 keV. In general, the following on-axis response files were used during the fits, epn_ff20_sdY9_filter.rmf for the PN camera and m1(2)_filterv9q20t5r6_all_15.rsp for the MOS1(2) cameras. Exceptions to this, e.g. during the analysis of RX J1325.5, are noted below. Fitted spectral parameters are quoted with 1$\sigma $ confidence limits on one interesting parameter. The spectral fits were typically performed simultaneously on the PN, MOS1(2) spectra. Again exceptions to this, e.g. during the analysis of RX J0847.2, are noted in the relevant sections.

Following Markevitch (1998), we have also investigated the impact of cooling cores on our fitted temperatures, and hence the measured Lx-Tx relation, by performing spectral fits after excising the central region. For all clusters except RX J1701 and RX J1325, a 50  h-150 kpc region was excised. For RX J1701 & 1325 a 120 kpc region was excised (see Sect. 4.7). The results from these fits are given at the beginning of Table 5. See Sect. 5.2.1 for a discussion.

3.7 Deriving bolometric luminosities

For the luminosity calculations, we adopted a physically meaningful circular aperture, with a virial radius, $r_{\rm v}$. We derived $r_{\rm v}$according the $T-r_{\rm v}$ relation of Evrard et al. (1996) for each cluster using the best fit temperatures from the spectral fits. We note that uncertainties both in the measured temperature and in the $T-r_{\rm v}$relation will introduce a systematic error in the derived Lxvalues. However, this should be insignificant, as very little cluster flux falls at radii close to the virial radius.

We integrated the background-subtracted counts inside $r_{\rm v}$ and then corrected for any cluster flux lost in areas under masked point sources or inter-chip gaps. The correction factors are given in row 12 of Table 5. These were derived in a model independent fashion; the count-rate was accumulated in 1 pixel wide annuli centered on the cluster. If any pixels within a specific annulus fell in chip gaps or in the point source masked, the annulus count-rate was scaled up accordingly. Next, we used XSPEC, together with the appropriate EPIC response functions, to determine the absorbed 0.3-4.5 keV flux within $r_{\rm v}$, using the best fit absorbed MEKAL model from our spectral analysis so that it yielded the measured $r<r_{\rm v}$ count rate (the MEKAL normalization factors are given in row 10 of Table 5). For comparison with equivalent ROSAT data such as SHARC or 160SD these fluxes (converted to the normal ROSAT 0.5-2 keV band) are listed in row 6 of Table 5. We then set the hydrogen column density to zero and re-determined the flux inside a pseudo bolometric band of 0.0-20 keV. From the unabsorbed bolometric flux, it is trivial to determine the bolometric luminosity for any given cosmological model. To aid comparisons with previous work, we give in, row 3 of Table 5, the Lxvalues appropriate for a $\Omega_{\rm M}=1$, $\Omega_{\Lambda}=0$, q0=0.5 model with H0=50 km s-1 Mpc-1 that has been most frequently used in the past as the parameter set for X-ray cluster studies.

Also following Markevitch (1998), we estimated the cooling flow corrected bolometric luminosities by summing counts in a ( $r_{\min}
\leq r \leq r_{\rm v}$) ring. For clusters other than RX J1701 and RX J1325, $r_{\min}$ = 50  h-150 kpc, while for RX J1701 & 1325 $r_{\min}=120$  h-150 kpc (see Sect. 4.7). To account for the non-cooling flow flux falling at $r\leq r_{\min}$, we applied a cluster specific renormalisation parameter calculated using the best fit $\beta $-model. We note that our correction approach differs from that of Markevitch (1998), therein a single multiplicative factor of 1.06 was used.

4 Individual clusters

We now discuss our analysis for each cluster in turn. We have ordered our discussion according to the Julian date on which they were observed (Table 1).

4.1 RX J1120.1+4318

Figure 6a shows the vignetting corrected, background subtracted and co-added (PN+MOS1+MOS2) 0.3-4.5 keV image of RX J1120.1+4318. Point sources removed during analysis are circled and their positions listed in Table A.1. Figure 19 shows the corresponding radial surface brightness distribution and best fit $\beta $ model; $\beta=0.77\pm0.03$, $\theta_{\rm c}=27.4\pm1.2$ arcsec ( $r_{\rm c}=209^{+9}_{-8}~ h^{-1}_{50}$ kpc, Table 5).

From a spectral extraction region with a radius of 145'' and using the in-field background subtraction technique, we measured the following temperature, metal abundance and redshift by fitting to the PN, MOS1 and MOS2 data simultaneously; Tx= 5.45 +0.26-0.35 keV, $Z=0.47 \pm 0.09$, $z=0.60 \pm 0.08$. We note the consistency between the redshift measured from the X-ray spectrum with the optically determined value (z=0.60; Romer et al. 2000). The overall $\chi_{\nu}$ of the spectral fit was 390/360 degrees of freedom. During the fit, the hydrogen column density was fixed at the Galactic value ($N_{\rm H}$   $=2.1\times 10^{20}$ atoms cm-2), but we note that, when left as a free parameter, its best fit value was very similar ($N_{\rm H}$   $=2.2 \pm 0.4 \times 10^{20}$ atoms cm-2).

\par {\includegraphics[angle=270,width=8.8cm,clip]{0687fig8.ps} }\end{figure} Figure 8: Spectral fit and residuals for RX J1120.1; black - PN, dark and light grey - MOS1 & 2 (double background subtraction).
Open with DEXTER

\par {\includegraphics[width=8.5cm,clip]{0687fig9.ps} }
\end{figure} Figure 9: The spatially resolved temperature profile of RX J1120.
Open with DEXTER

Of the eight clusters in our sample, this object has the greatest signal-to-noise ratio. We have taken advantage of this to perform additional spectral analyses. First we have been able to use the double subtraction technique to investigate how the adopted background subtraction method impacts the spectral fits. Doing so, we derive Tx=5.6+0.25-0.3 keV, $Z=0.43 \pm 0.06$, $z=0.605 \pm
0.08$ (Fig. 8). The overall $\chi_{\nu}$ of the spectral fit was 330/307 degrees of freedom. During the fit, the hydrogen column density was again fixed at the Galactic value ($N_{\rm H}$   $=2.1\times 10^{20}$ atoms cm-2), but when left as a free parameter, its best fit value was ($N_{\rm H}$   $=2.2 ^{+0.2}_{-0.4}\times
10^{20}$ atoms cm-2). It is clear, therefore, that the spectral fits are not significantly changed by the choice of background subtraction technique (see also Sect. 4.2).

We have also investigated how the spectral fits differ when we treat the PN and MOS data separately. We find them to be in good agreement; $T_{{x},\rm PN}=5.30 \pm 0.6$ keV and $T_{{x},\rm MOS}=5.7
\pm 0.8$ keV. This is encouraging since, in two other cases (RX J1325.5 and RX J0847.2), we do not have access to data from all three cameras. We have also been able to subdivide the spectral extraction region into three radial bins, and determined a crude temperature profile (Fig. 9). The profile is essentially flat, indicating that there is not a "cooling flow'' region at the cluster core (this conclusion is supported by the absence of a central spike in the surface brightness profile). We note that for the PN/MOS comparison and for the radial profile, we used the double subtraction technique to account for the background.

4.1.1 Comparison with Arnaud et al. (2002)

The RX J1120.1 observation described here have been previously analyzed and interpreted in an earlier paper in this series (Arnaud et al. 2002). The analysis procedures developed in that paper have formed the basis of the analysis of the eight clusters described here. For example, the vignetting correction technique used in Arnaud et al. (2002) has been subsequently implemented in SAS as the EVIGWEIGHT routine mentioned above. However, there are certain differences in our respective data analysis techniques, even with regard to the RX J1120.1 observation. For example, we use updated calibration information and an updated SAS processing version. We also used a different method (3-sigma clipping) for the GTI filtering. We have also adopted the in-field background subtraction method as the standard for our cluster analysis (in Arnaud et al. 2002, only the double subtraction technique was used). Finally, we note our use of the both single and double PN events, compared to the selection of single events by Arnaud et al. (2002). In general the revised analysis of RX J1120.1 has yielded very similar results to Arnaud et al. (2002) and we concur that, given its isothermal temperature profile and the absence of significant substructure, RX J1120.1 is consistent with being a relaxed cluster. It is noteworthy, however, that in Arnaud et al. (2002) the mean temperature values determined using only MOS data differed from those determined using only PN data by more than 1 keV; $T_{{x}, \rm MOS}=5.8^{+1.0}_{-0.7}$ keV and $T_{{x},\rm PN}=4.5^{+0.8}_{-0.5}$keV. By contrast, we measured temperatures that differed by only 8% (see above). This improvement is most likely attributable to the improved calibration that has become available since Arnaud et al. (2002) was published. For completeness we compare in Table 4 the values derived for various fitted parameters in the two analyses. For consistency with Arnaud et al. (2002), we quote the mean temperature derived using the double subtraction technique (in Table 5 we quote the value from based on the in field subtraction technique).

Table 4: Comparison of parameter fits to the RX J1120.1 observation by Arnaud et al. (2002) (first line) and this work (second line). Here Lx refers to the bolometric luminosity.

4.2 RX J1334.3+5030

Figure 6b shows the image of the RX J1334.3+5030 field. Point sources removed during analysis are circled and their positions listed in Table A.1. Figure 19 shows the corresponding radial surface brightness distribution and best fit $\beta $ model; $\beta=0.66~\pm~0.02$, $\theta_{\rm c}=20 \pm1 $ arcsec ( $r_{\rm c}=154\pm10~ h^{-1}_{50}$ kpc, Table 5).

From a spectral extraction region with a radius of 120'' and using the in-field background subtraction technique, we measured the following temperature and metal abundance; Tx=5.20+0.26-0.28 keV, $Z=0.15 \pm 0.08$ (Fig. 10). The overall $\chi_{\nu}$ of the spectral fit was 464/473 degrees of freedom. The hydrogen column was fixed at the Galactic value ($N_{\rm H}$   $=1.05\times 10^{20}$ atoms cm-2), but we note that, when left as a free parameter, its best fit value was very similar ($N_{\rm H}$   $=0.8\pm 0.5 \times 10^{20}$ atoms cm-2). Likewise, the redshift was fixed at its optically derived value of z=0.62 (Romer et al. 2000), but when left as a free parameter, its best fit value was $0.63 \pm0.02$. The derived values for the bolometric luminosity and absorbed flux inside the viral radius are given in Table 5. We were also able to apply the double subtraction technique to these data and derived a consistent mean temperature value; $T_{x}=5.05\pm0.3$ keV.

\par {\includegraphics[angle=270,width=8.8cm,clip]{0687fig10.ps} }\end{figure} Figure 10: Spectral fit and residuals for RX J1334.3; black - PN, dark & light gray - MOS1 & 2 (in-field background subtraction).
Open with DEXTER

4.3 RX J0337.7-2522

Figure 6c shows the vignetting corrected, background subtracted image of the RX J0337.7-2522 field. Two excised sources within the spectral extraction radius, XMMU J033742.9-252208 & XMMU J033745.9-252206, are tentatively associated with two stellar objects on the USNO Catalogue. The former is U0600-01432100, also identified as blue stellar object PHL4470, the latter is identified as U0600-01432383 (a 16.5 r magnitude object). Figure 19 shows the radial surface brightness distribution and best fit $\beta $ model; $\beta=0.76^{+0.08}_{-0.04}$ and $r_{\rm c}=19.4\pm2.5$ arcsec. ( $r_{\rm c}=145\pm18~ h^{-1}_{50}$ kpc, Table 5).

\par {\includegraphics[angle=270,width=8.8cm,clip]{0687fig11.ps} }\end{figure} Figure 11: Spectral fit and residuals for RX J0337.7; black - PN, dark & light grey - MOS1 & 2.
Open with DEXTER

From a spectral extraction region with a radius of 120'' and using the in-field background subtraction technique, we measured the following temperature, metal abundance and redshift; $T_{x}=2.6 \pm~0.35$ keV, $Z=0.38\pm 0.09$, $z=0.57\pm0.3$ (Fig. 11). The overall $\chi_{\nu}$ of the spectral fit was 193/214 degrees of freedom. We note the consistency between the redshift measured from the X-ray spectrum with the optically determined value (z=0.577; Burke et al. 2003). The hydrogen column density was fixed at the Galactic value ($N_{\rm H}$   $=0.99\times 10^{20}$ atoms cm-2), when left as a free parameter, its best fit value was very similar ($N_{\rm H}$   $=8.7^{+6.4}_{-5.4} \times 10^{19}$ atoms cm-2). The derived values for the bolometric luminosity and absorbed flux inside the viral radius are given in Table 5.

4.4 RX J0505.3-2849

The RX J0505.3-2849 field is shown in Fig. 6d. Four point sources removed during analysis are circled and their positions listed in Table A.1. Tentative identifications of these sources with objects in the UK APM survey (Maddox et al. 1990) suggest an absolute astrometric accuracy for our observations of $\sim $2 arcsec. Figure 19 shows the corresponding radial surface brightness distribution and best fit $\beta $ model; $\beta=0.66^{+0.05}_{-0.04}$ and $r_{\rm c}=22.8\pm2.4$ arcsec ( $r_{\rm c}=164\pm17~ h^{-1}_{50}$ kpc, Table 5).

From a spectral extraction region with a radius of 120'' and using the in-field background subtraction technique, we measured; $T_{x}=2.5~\pm0.3$ keV, $Z=0.17 \pm
0.08$ (Fig. 12). The overall  $\chi_{\nu}$ of the spectral fit was 279/248 degrees of freedom. The fixed Galactic hydrogen value was $N_{\rm H}$   $=1.5\times 10^{20}$ atoms cm-2), compared with a free fit parameter, of $N_{\rm H}$   $=1.1\pm 0.6 \times 10^{20}$ atoms cm-2). The redshift was fixed at its optically derived value of z=0.51 (Burke et al. 2003), but when left as a free parameter, its best fit value was $z=0.53 \pm0.04$. As a test, we have also performed a spectral fit without excluding the four point sources. We find that, within the errors, the fitted temperature was unchanged.

\end{figure} Figure 12: Spectral fit and residuals for RX J0505.3; black - PN, dark & light grey - MOS.
Open with DEXTER

4.5 RX J0847.2+3449 (and XMMU J084701.8+345117)

Figure 7a shows the image of RX J0847.2+3449. Of the eight clusters in our sample, RX J0847.2+3449 is noteworthy as being the only one that is not a SHARC cluster; it was selected instead from the 160SD (Vikhlinin et al. 1998). RX J0847.2 was observed during a period when the PN camera was temporarily disabled by a hardware fault, and the requested exposure duration was partly compensated by an extension of the MOS observation. To improve the contrast of the cluster against the nearby bright QSO (PG 0844+34), only the inner part of the image, that covering central CCD of the MOS cameras, is shown in Fig. 7a. We note that the mirror scattering is sufficiently low that no flux from PG 0844+34 impedes our analysis. The bright point source $\sim $3 arcmin west of the RX J0847.2 cluster in Fig. 7a was identified with a second, fainter, quasar, QSO 0843+349 at z=1.57. The image also shows evidence for a nearby fainter extended object to the NW, which we identify as a cluster candidate with provisional designation XMMU J084701.8+345117. Vikhlinin (private communication) claims to have found a concentration of faint galaxies in this region, but to date no galaxy redshifts are available. Neither QSO 0843+349 nor XMMU J084701.8+345117 fall inside our spectral extraction region. The source detection software picked out four sources in that region, none of which have counterparts in the NED catalogs. These sources were excluded from the spatial and spectral analysis, they are circled in Fig. 7 and their positions are listed in Table A.1.

Figure 20 shows the radial surface brightness distribution for RX J0847.2 and the best fit $\beta $ model; $\beta=0.81\pm 0.07$ and $r_{\rm c}=42^{+4}_{-4}$ arcsec. ( $r_{\rm c}=307\pm30~ h^{-1}_{50}$ kpc, Table 5). From a spectral extraction region with a radius of 120'' and using the in-field background subtraction technique, we measured Tx=3.62+0.58-0.51 keV, $Z=0.30 \pm 0.28$ (Fig. 13). The overall $\chi_{\nu}$ of the spectral fit was 152/180 degrees of freedom. The fixed hydrogen column density was $N_{\rm H}$   $=3.2\times 10^{20}$ atoms cm-2, as a free parameter, its best fit value was $N_{\rm H}$   $=2.8\pm 0.14 \times
10^{20}$ atoms cm-2. The redshift was fixed at its optically derived value of z=0.56but its best fit value was $z=0.54 \pm0.04$. As a test, we have also performed a spectral fit without excluding the four point sources. Doing so changed the fitted temperature by 0.1 keV and the measured flux by $\leq$10%.

We have also estimated the temperature of the second cluster in the field, XMMU J084701.8+345117. For this we used a spectral extraction radius of 80 arcsec and the in field background subtraction technique. We do not have an independent estimate of the cluster redshift, and it cannot be constrained by the X-ray spectrum, so we have assumed that the cluster lies at the same redshift as RX J0847.2 (z=0.56). We have also fixed the metal abundance to be Z=0.3 and the hydrogen column density to the same Galactic value. We measure a cluster temperature of Tx=1.8+1.0-0.4 keV and a corresponding ( $r<r_{\rm v}$, 0.5-2 keV) flux and bolometric luminosity of $1.3\times$ 10- 14 erg cm-2 s-1 and $1\times~10^{43}$ erg s-1 respectively. These results are consistent with measured low red-shift Lx-Tx relations.

\par {\includegraphics[angle=270,width=8.8cm,clip]{0687fig13.ps} }\end{figure} Figure 13: Spectral fit and residuals for RX J0847.2 (MOS data only).
Open with DEXTER

4.6 RX J1325.5-3826 (and RX J1325.0-3814)

The data for this cluster were obtained by kind permission of Guaranteed Time holder M Watson, from an observation of the nearby object IRAS 13224-3809. Unfortunately for our own purposes, this bright source had been observed in the MOS cameras in a "Window'' mode of readout on the central CCDs in order to minimize effects of photon pile-up, and consequently the area of focal plane containing the SHARC cluster was not exposed. Therefore, our spatial and spectral analysis is restricted to the PN camera data only. Figure 7b shows the vignetting corrected, background subtracted 0.3-4.5 keV PN image of the region surrounding IRAS 13224-3809. Clusters RX J1325.5-3826 and RX J1325.0-3814 (see below) are labeled, as is the streak corresponding to the "out of time'' events from IRAS 13224-3809. IRAS 13224-3809 is bright enough that the "spider'' like structure in the PSF is visible around it. Given the proximity of the cluster to the IRAS source, the surface brightness profile for RX J1325.5-3826 (Fig. 20) was extracted from manually selected clean regions that avoided the IRAS 13224-3809 PSF.

For the $\beta $ model fitting, a PSF appropriate for the PN at the off-axis angle of the cluster was used (for all the other clusters in our sample, we had used the on-axis PSF). The best fit parameters in the spatial fit are $\beta=0.64^{+0.09}_{-0.07}$ and $\theta_{\rm c}=17.3 \pm 3.0$ arcsec ( $r_{\rm c}=115 \pm 20~ h^{-1}_{50}$ kpc, Table 5). There is evidence for possible excess brightness in the inner bins of the radial profile. This may be in an indication of a central cooling core, but we do not have adequate signal to noise to confirm this spectroscopically. Following the procedure adopted for RX J1701.3 (see below) we also performed a $\beta $ fit after exclusion of the core 120 kpc and found a value for $\beta $ of $0.71\pm0.05$, given a fixed $r_{\rm c}$ of 200  h50-1 kpc (we fixed this value arbitrarily to match typical cluster values).

We used a 90'' radius extraction region to fit the spectrum from this cluster. This is a smaller region than was used for the other clusters (120'' or 145'', Table 5) because of the proximity of IRAS 13224-3809. Even though the proton background was low enough in this exposure to permit spectral measurements using the double subtraction technique, we decided to use instead only the in-field background subtraction technique. This was because we wanted to mimic the off-axis angle dependence of the point source contamination. Rather than using an annulus around the cluster to determine the background, we used two source and streak free circular regions at the same off-axis angle as the cluster. From the background subtracted PN spectrum, we measured the following temperature, metal abundance and redshift; Tx=4.15+0.4-0.3 keV, Z=0.31 +0.19-0.17, $z=0.44\pm{0.01}$(Fig. 14). We note the consistency between the redshift measured from the X-ray spectrum with the optically determined value (z=0.445; Burke et al. 2003). The overall $\chi_{\nu}$ of the spectral fit was 229/242 degrees of freedom. The fixed and free hydrogen column density values were $N_{\rm H}$   $=4.8\times 10^{20}$ atoms cm-2 and 4.5 $^{+0.9}_{-0.8}
\times 10^{20}$ atoms cm-2 respectively.

\par {\includegraphics[angle=270,width=8.8cm,clip]{0687fig14.ps} }\end{figure} Figure 14: Spectral fit and residuals for RX J1325.5 (PN data only).
Open with DEXTER

4.6.1 RX J1325.0-3814
For completeness, we note that another, fainter, Southern SHARC cluster lies within the same field of view. The PN image of this cluster, RX J1325.5-3814, can be seen in Fig. 7b. The cluster is also visible in the MOS data (not shown) because it lies outside the region blanked by the "Window'' mode. The cluster is too faint and too far off axis (i.e. where the PSF is poorly defined) to permit a detailed spatial analysis, however we were able to make a spectral analysis. For this we generated off-axis response matrices for each camera using the SAS ARFGEN package. From a spectral extraction region with a radius of 90'', and using the in-field background subtraction technique, we measured the following temperature and redshift ($N_{\rm H}$ fixed at the same galactic value of $4.8\times 10^{20}$ atoms cm-2 and abundance at 0.3) $T_{x}=3.2
\pm0.4$ keV and a redshift of $z=0.29 \pm0.02$. We note the consistency between the redshift measured from the X-ray spectrum with the optically determined value (z=0.296; Burke et al. 2003).

The cluster is observed close to the edge of the outer CCDs, preventing the use of an extraction region large enough to correspond with the virial radius. Within a radius $\sim $500  h50-1 kpc, the measured flux and bolometric luminosity are 6.1 ($\pm$0.3) 10-14 erg cm-2 s-1 (0.5-2 keV) and $L_{x}=1.1 \pm 0.23\times10^{44}$ erg s-1 respectively. This combination of Tx and Lx are not consistent with the measured Lx-Tx relation for nearby clusters. For example, for Tx=3.2 keV, one would expect a bolometric luminosity of $L_{x}= 2 \times
10^{44}$ erg s-1 based on the Markevitch (1998) relation, and is probably due mainly to the missing flux outside our spectral extraction radius.

4.7 RX J1701.3+6414 (and Abell 2246)

This cluster was observed during an exposure scheduled near the end of an orbit, where the spacecraft was approaching the particle belts. Not only were the spacecraft operations terminated prematurely, but most of the exposure was dominated by high soft proton background from the edge of the particle belts. Fortunately this is one of the brightest clusters in our sample and we were still able to produce acceptable quality images and spectra.

Figure 7c shows the image of the RX J1701.3+6414 field. Two nearby sources are marked in the figure; a QSO (HS1700+6416) and a known cluster (A2246).

The XMM-Newton data indicate an anomalously low $r_{\rm c}$ value in the $\beta $-fit, and this is supported by CHANDRA data (Vikhlinin et al. 2002) where this cluster was observed with a suspected cooling flow central brightness enhancement. We therefore excluded the central bins (120  h50-1 kpc) and fixed the $r_{\rm c}$ to the CHANDRA value of 0.5 arcmin (204 kpc). Figure 20 shows the corresponding radial surface brightness distribution and resulting best fit $\beta $ model; $\beta=0.64 \pm{0.05}$, $\theta_{\rm c}=30$ arcsec. This value for $\beta $ is consistent with the CHANDRA measurement ( $0.62 \pm 0.03$).

From a spectral extraction region with a radius of 120'', we measured the following temperature, metal abundance and redshift (without core excision); Tx=4.5+1.5-1.0 keV, $Z=0.24 \pm 0.2$, $z=0.44\pm0.02$ (Fig. 15). The overall $\chi_{\nu}$ of the spectral fit was 58/53 degrees of freedom. The hydrogen column density was fixed at the Galactic value ($N_{\rm H}$   $=2.6\times 10^{20}$ atoms cm-2), and as a free parameter, its best fit value was very similar ($N_{\rm H}$   $=2.6\pm 1.5 \times 10^{20}$ atoms cm-2).

After excluding the 120 kpc cooling flow region we obtain $T_{\rm
x}=4.8^{+1.9}_{-1.3}$ keV. The derived values for the bolometric luminosity and absorbed flux inside the viral radius are given in Table 5. We were unable to apply the double subtraction technique to these data because of the enhanced proton background during the observation (see Fig. 5). The poorer quality of data in this observation limits the quality of determination of T and L. Agreement with the CHANDRA data is acceptable for the temperature, but not the luminosity ($5.8\pm0.5$ keV and $15.7\times 10^{44}$ erg/s respectively).

\par {\includegraphics[angle=270,width=8.8cm,clip]{0687fig15.ps} }
\par\end{figure} Figure 15: Spectral fit and residuals for RX J1701.3; black - PN, dark & light grey - MOS.
Open with DEXTER

4.7.1 Abell 2246
For completeness we also offer our interpretation of the Abell 2246 spectrum. From a spectral extraction region with a radius of $90 \hbox{$^{\prime\prime}$ }$, and using the in-field background subtraction technique, we measured the following temperature, metal abundance and redshift by fitting to the PN, MOS1 and MOS2 data simultaneously; Tx=2.7 +0.6-0.5keV, $Z=0.32\pm0.13$, $z=0.22\pm0.04$. We note the consistency between the redshift measured from the X-ray spectrum with the optically determined value (z=0.225; Struble & Rood 1991). Within a radius $\sim $415  h50-1 kpc, the flux and bolometric luminosity are $3.4\pm0.2\times10^{-13}$ erg cm-2 s-1 (0.5-2 keV) and $2.1\times10^{44}$ erg s-1 respectively. Despite the rather restricted radius from which the luminosity is determined, the local Lx-Txrelation would predict a slightly lower luminosity than this estimate.

4.8 RX J1354.2-0222

The RX J1354.2-0222 field is shown in Fig. 7d. A point source that was removed during analysis is circled and its position listed in Table A.1. Figure 20 shows the corresponding radial surface brightness distribution and best fit $\beta $ model that gives $\beta=0.68 \pm0.06$ and $r_{\rm c}=33.6^{+4.9}_{-3.6}$ arcsec ( $r_{\rm c}=248 ^{+36}_{-26}~h^{-1}_{50}$ kpc, Table 5).

From a spectral extraction region with a radius of 120'', and using the in-field background subtraction technique, we measured; Tx=3.66+0.6-0.5 keV, $Z=0.25
\pm 0.14$, (Fig. 16). The overall $\chi_{\nu}$ of the spectral fit was 120/147 degrees of freedom. The fixed Galactic hydrogen column value was $N_{\rm H}$   $=3.4\times 10^{20}$ atoms cm-2, but $N_{\rm H}$   $=3.2\pm 1.2 \times 10^{20}$ atoms cm-2) when free. The redshift best fit value was $0.53
\pm0.04$. Similarly to RX J1701.3, we were unable to apply the double-background subtraction technique to these data because of the enhanced proton background (see Fig. 5).

\par {\includegraphics[angle=270,width=8.8cm,clip]{0687fig16.ps} }
\end{figure} Figure 16: Spectral fit and residuals for RX J1354.2 black - PN, dark & light grey - MOS.
Open with DEXTER

5 Discussion

5.1 Luminosity-temperature relations

Table 5 summarizes our analysis of the eight cluster observations described above. The mean elemental abundance ( $Z=0.28\pm
0.08$), core radius ( $r_{\rm c}=195\pm50~ h_{50}^{-1}$ kpc) and $\beta $ ( $\beta=0.70\pm0.05$), are consistent with typical values measured at low redshift (e.g. Fukazawa et al. 1998; Jones & Forman 1999; Mohr et al. 1999). The clusters were derived from well understood surveys (SHARC and 160SD) and thus should be representative of the cluster population as a whole at these redshifts. We have used our observations to determine the Lx-Tx relation at 0.45<z<0.62 and to investigate evolutionary effects. We plot luminosity versus temperature, with one sigma errors, for the 8 cluster targets in Fig. 17, using the values $T^{\rm cf}$, $L^{\rm cf}$ of Table 5. The errors are clearly dominated by those of the temperature measurements.

Characterizing the Lx-Tx relation as

\begin{displaymath}L_{x} = L_{6} \left(\frac{T}{6 ~\rm keV} \right)^{\alpha} %
\end{displaymath} (1)

and assuming an EdS cosmology, we find $\alpha=2.7\pm0.4$, and $L_{6}=15.9 ^{+7.9}_{-5.5} \times 10^{44}$ erg s-1 using the bisector variant of the BCES fitting package (Akritas & Bershady 1996). This relation is shown as the solid line in Fig. 17. The slope of the relation is similar to most previous measurements, e.g. for an EdS cosmology, $\alpha=2.64\pm0.27$(Markevitch 1998); $\alpha=2.33\pm0.43$ (Allen & Fabian 1998); $\alpha=2.88\pm0.15$ (Arnaud & Evrard 1999); $\alpha=2.47\pm0.14$(Ikebe et al. 2002); $\alpha=2.82\pm0.32$ (Novicki et al. 2002).

\par {\includegraphics[width=7.8cm,clip]{0687fig17.ps} }\end{figure} Figure 17: Lx-Tx relation for an EdS model (q0=0.5) and H0=50 km s-1 Mpc-1. The solid line is the best fit for our sample after cooling flow excision, the dotted line is the relation from Markevitch (1998) for low red-shift clusters.
Open with DEXTER

\par {\includegraphics[width=8.3cm,clip]{0687fig18.ps} }\par\end{figure} Figure 18: Lx-Tx relation for H0=70, $\Omega _{\Lambda }=0.7, \Omega _{\rm M}=0.3$. The dotted horizontal lines are the CHANDRA sample reported by Vikhlinin et al. (2002), and the dashed diagonal line the Markevitch (1998) Lx-Txrelation corrected to this cosmology. The solid crosses indicate the results from the clusters in our study, corrected for this cosmological model, and the diagonal solid line is the best fit for our data.
Open with DEXTER

To characterize the possible evolution with z, we assume

\begin{displaymath}L_{x} = L_{6}\left(\frac{T}{6 ~\rm keV} \right)^{\alpha} (1+z)^{A }.
\end{displaymath} (2)

Adopting the low redshift Lx-Tx relation from Markevitch (1998) ( $\alpha=2.64\pm0.27$, $L_{6}=12.44\pm 1.08 \times 10^{44}$ erg s-1) we performed a $\chi^{2}$ minimization on our data points corrected by a (1+z)A term. The Markevitch (1998) relation was derived from ROSAT and ASCA observations of 30 clusters at 0.04<z<0.1 ( $\bar z=0.05$) and is plotted as a dotted line in Fig. 17. On average, the XMM-Newton data points lie away from the Markevitch (1998) relation, suggesting an evolutionary effect whereby clusters of the same temperature were more luminous in the past. The best fit value for A is 0.68$\;\pm\;$0.26 ( $\chi_{\rm red}=1.2$). We have also determined the BCES best fit for our sample when all the points are corrected with this evolution to a single redshift (low redshift z=0.05 of Markevitch) bin, and find $\alpha=2.71~\pm~0.22$, and $L_{6}=12.5^{+4.9}_{-3.5} \times 10^{44}$ erg s-1. We note that we have chosen to compare with Markevitch (1998), rather than with the recent study of 82 clusters with ASCA temperatures by Ikebe et al. (2002), because the Markevitch (1998) method to derive cluster temperatures is closer to our own (Sect. 5.2.3), thus we try to minimise any systematic effects from different data treatment. In Fig. 18 we re-plot our data and the Markevitch (1998) relation after converting them to a common cosmology of H0=70 km s-1 Mpc-1, $\Omega _{\Lambda }=0.7$, $\Omega_{\rm M}=0.3$. In this cosmology we find A =1.52+0.26-0.27.

\end{figure} Figure 19: Surface brightness distribution of clusters, compared with best $\beta $-model fit (dotted line). Reading clockwise from top left: RX J1120.1, RX J1334.3, RX J0337.7 and RX J0505.3.
Open with DEXTER

For comparison, we include the data points obtained by Vikhlinin et al. (2002) from a compilation of Chandra observations of high-z clusters ( 0.39<z<1.26) in Fig. 18 (after adjusting to H0=70 km s-1 Mpc-1). From this figure, it is clear that both the XMM-Newton and Chandra points lie away from the Markevitch (1998) relation and that there is no obvious systematic offset between the XMM and Chandra data (see Sect. 5.2.3). Vikhlinin et al. (2002) also used Markevitch (1998) as a low redshift benchmark to investigate Lx-Tx evolution. From their data, Vikhlinin et al. (2002) found $A=0.6\pm0.3$ for $\Omega_{\Lambda}=0$, $\Omega_{\rm M}=1$ and $A=1.5\pm0.3$ for a $\Omega _{\Lambda }=0.7, \Omega _{\rm M}=0.3$ cosmology.
\end{figure} Figure 20: Surface brightness distribution of clusters, compared with best $\beta $-model fit (dotted line). Reading clockwise from top left: RX J0847.2, RX J1325.5, RX J1701.1 and RX J1354.2.
Open with DEXTER

The value of A has important implications for our understanding of structure formation and cluster evolution. In a self-similar model (Kaiser 1981) cluster X-ray properties are driven by gravitational processes, such as shock heating, and a value of $A\sim1.5$ would be expected in most cosmological models. A self-similar EdS model predicts exactly 1.5, and the value is only slightly different in a low density concordance model, the correction due to the redshift dependence of the virial density (Bryan & Norman 1998) being small. In a model where the cluster X-ray properties are influenced by non-gravitational processes, such as energy injection by AGN's or supernovae, lower values of A are predicted (e.g. Tozzi & Norman 2001). The evolution deduced from our observations assuming an EdS model is significantly below the predicted value. Adopting a low density concordance model, on the other hand, leads to a value consistent with predictions. We conclude that there is now evidence from both XMM-Newton and Chandra for an evolutionary trend in the Lx-Tx relation.

Previous studies of Lx-Txevolution based on either ASCA (e.g. Mushotzky & Scharf 1997; Allen & Fabian 1998; Novicki et al. 2002; Sadat et al. 1998), ROSAT PSPC (e.g. Fairley et al. 2000) or Chandra (e.g. Borgani et al. 2001; Holden et al. 2002) temperature measurements found results consistent with no evolution (Sadat et al. 1998). When these studies quote values for A, these values are generally smaller than those measured by us or Vikhlinin et al. (2002), but are usually still at least one sigma away from A =0. Sadat et al. (1998) found a positive evolution $A =0.5\pm 0.3$ for $\Omega_{\rm M}=1,\Omega_{\Lambda}=0$ and Novicki et al. (2002) $A =1.1\pm1.1$ for $\Omega_{\rm M}=1,\Omega_{\Lambda}=0$, and $A=2.1\pm1.0$ for $\Omega _{\Lambda }=0.7$, $\Omega_{\rm M}=0.3$ respectively. It is clear that we are beginning to probe evidence that Lx-Tx evolution, although many more clusters need to be studied, and systematic biases (see below) examined in detail, before solid conclusions regarding structure formation models can be drawn.

5.2 Data treatment and systematics

We have made every effort to ensure that the Lx-Tx relation presented above is robust. We have adopted a uniform approach to the spatial and spectral analysis of the eight clusters in the sample. We have excised regions with clear signs of point source contamination before extracting spectra and surface brightness profiles. We have not used $\beta $-models to calculate total cluster count rates, but have rather summed up the counts within a viral radius. We have tried to ensure that our approach to background subtraction does not bias the measured cluster parameters (see below). We have used measured temperatures to make conversions between count rates and fluxes and, where possible, we have used data from all three EPIC cameras to derive Lx and Tx values. For RX J0847.2 and RX J1325.5, data were not available from all three cameras. However, we are confident that the derived quantities for these two clusters are robust, as we have shown, using the RX J1120.1 observation, that PN-only or MOS-only spectra are in excellent agreement.

The background subtraction has been complicated by the extended nature of the targets and the high proton flare contamination in some of the observations. We used background template files during the spatial analysis. These files covered a limited energy range (0.3-4.5 keV), to minimize any calibration uncertainties. They were filtered using the same rate filtering criteria as their respective cluster observation and were normalized to it using the particle background rates in the CCD regions outside the FOV. For various reasons (see Sect. 3.3), we did not use the background templates for the spectral analysis. Instead, we adopted an in-field background subtraction technique. With this method comes the concern that the proton background can be over-vignetted. We have checked for this by applying a different method, the double subtraction technique, to the RX J1120.1 and RX J1334.3 observations. We found that the results from both techniques agree within the statistical errors. For the two observations with the worst proton flare contamination, RX J1354.2 and RX J1701.3, we applied a 2% scaling to the in-field background spectrum to compensate for potential over-vignetting of the proton background. In Fig. 5 we demonstrate that this over-vignetting should not present a problem for the other clusters in our sample.

Despite all these quality controls, we still cannot rule out the possibility of there being some systematic bias in our results. We discuss possible sources of such bias and the potential impact on our claimed Lx-Tx evolution below.

5.2.1 Impact of cooling core clusters

Cooling cores have the joint effect of boosting the cluster luminosity and of reducing the measured global cluster temperature. These effects have been shown to introduce scatter in the Lx-Tx relation and to modify its slope (Markevitch 1998; Allen & Fabian 1998). Various approaches have been taken to alleviate the effects of cooling cores on measured Lx-Txrelations, e.g. by excluding cooling core clusters from the analysis, Arnaud & Evrard (1999), by fitting two temperature models (Ikebe et al. 2002), by masking the central regions from obvious cooling core clusters (Vikhlinin et al. 2002), or by masking the central regions from all clusters (Markevitch 1998). Based on the images (Figs. 6 and 7) and radial profiles (Figs. 19 and 20 ) of the eight clusters in our sample, we do not expect cooling cores to have much, if any, impact on the measured Lx and Tx values; only RX J1701.3 and RX J1325.5 show evidence of a central surface brightness excess. We have investigated this by comparing the results before and after excising a circular region[*] from the cluster cores.

The results of the re-analysis can be found in Table 5, (rows 4 and 5). Temperatures changed very little - always within the statistical error - after we excised the central r=50  h50-1kpc region. In two cases, (RX J1120.1 and RX J1701.3), the luminosity changes by more than $1\sigma$, but still by only 5% and 9% respectively. Without the excision the best fit Lx-Txrelation is $L_{6}=17.6 ^{+5.3}_{4.1}\times
10^{44}$ erg s-1  and $\alpha=2.66\pm0.25$ for an $\Omega_{\Lambda}=0.0, \Omega_{\rm M}=1.0$, H0=50 km s-1 Mpc-1 cosmology at a typical redshift $z \sim 0.55$. In summary, we find the Lx-Tx relations before and after the core excision are consistent within the errors. However, we caution that some residual cooling flow contamination may still remain after this excision; $50 ~h_{50}^{-1}~\rm kpc \simeq 14''$, corresponding to an encircled energy fraction of a little under 50% of the on-axis mirror PSFs.

5.2.2 Incompleteness and flux errors in the ROSAT catalogs

Our Lx-Tx evolution result implies that high redshift clusters in our sample are more luminous than clusters of the same temperature at lower redshift. We interpret this as evidence for a general evolutionary trend in the Lx-Txrelation, but it might also reflect an underlying selection bias in the cluster catalogs from which the sample was drawn, e.g. if the SHARC and 160SD surveys were biased toward merging systems. It is possible that the luminosity and temperature of clusters can change dramatically during the process of a major merging event (Ricker & Sarazin 2001). Based on the XMM-Newton images of the 8 clusters studied here (Figs. 6 and 7), such a bias does not appear to exist; only RX J1354 has evidence for a possible double component, and excluding this from the best fit Lx-Txrelation changes the luminosity normalisation only $\sim $2%. As long as the angular resolution of observations allows sufficient discrimination against obvious mergers, then presumably an ensemble sample of local and distant clusters are similarly affected by merger boosting of luminosity.

Moreover, we find no evidence for any incompleteness in the SHARC and 160SD surveys, at least in the small area covered by these observations. No new clusters were found in the fields surrounding the seven SHARC clusters in our sample. In the field of the 160SD cluster, RX J0847.2, we have discovered a new X-ray cluster (XMMU J084701.8+345117). However this cluster was too faint to have been included in the original 160SD survey. We estimate the ( $r<r_{\rm v}$) flux of this object to be $1.3\times 10^{-14}$ erg cm-2 s-1 (based on a fitted temperature of Tx=1.8 keV and an assumed redshift of z=0.56), compared to the 160SD flux limit of $3.7\times 10^{-14}$ erg cm-2 s-1, so this cluster is too faint to have been included in the original 160SD survey. In summary we do not find evidence for any intrinsic bias in the SHARC or 160SD survey selection functions that could mimic Lx-Tx evolution.

Another possible way to mimic the claimed evolution would be to systematically overestimate the Lx values. We have investigated this by comparing our flux measurements to those published in several works (Romer et al. 2000; Burke et al. 2003; Vikhlinin et al. 1998). We report here that a systematic offset does appear to exist, but in the opposite direction to that needed to mimic evolution. The XMM-Newton fluxes quoted in Table 5 are all lower than previously determined values, typically by 20% but as much as by 40% in the case of RX J0505.3. We note that this discrepancy persists regardless of the ROSAT catalog from which the cluster was selected. The average discrepancy between the Southern SHARC (4 clusters), Bright SHARC (3 clusters) and 160SD (3 clusters) fluxes and the XMM-Newton fluxes is 22%, 21% and 16% respectively.

Improvements in point source rejection, count rate to flux conversions and the signal to noise have meant that the typical statistical error on the XMM-Newton determined fluxes in Table 5 is ${\simeq}
3\%$, compared to ${\simeq} 10\%$ for the ROSAT determinations. Therefore, even accounting for a possible offset in the absolute cross calibration of the two observatories, it is clear that a systematic flux offset does exist. There are likely to be several reasons for this offset. Of these, the use of an incorrect spatial model, to convert between an aperture flux and a total flux, is probably the most significant. As shown in Romer et al. (2000), the use of a generic ( $\beta=0.67$, $r_{\rm c}=250~ h_{50}^{-1}$ kpc), rather than best fit, $\beta $-model resulted in a typical enhancement of 10% in the total cluster flux. In our XMM-Newton analysis, we do not use the $\beta $-model fits to determine the total flux. Instead we directly sum up the counts within a virial radius. For some clusters this difference in approach can explain the entire ROSAT to XMM-Newton discrepancy. For example, in the case of RX J1334.3, our estimate of $r_{\rm v}$ is almost identical to the r80 radius within which 80% of the total flux from a generic $\beta $-model would fall. Romer et al. (2000) divided the ROSAT counts within r80 by 0.8 to estimate the total count rate. For RX J1334.3 this division would artificially boost the flux by 25% which more than accounts for the 22% the discrepancy between the XMM-Newton and Bright SHARC values.

For the four Southern SHARC clusters in our sample, we have recalculated the total ROSAT PSPC fluxes using the $\beta $ and $r_{\rm c}$values given in Table 5 (in Burke et al. 2003 generic values were used; $\beta=0.67$, $r_{\rm c}=250~ h_{50}^{-1}$ kpc). By doing so, we improve the consistency between the ROSAT and XMM-Newton fluxes to within 5% for RX J1354.2 and RX J1325.5 and within 15% for RX J0337.7. For RX J0505.3 there remains a large, 35%, discrepancy which requires further investigation.

For the three Bright SHARC clusters in our sample, we have gone back to the original Romer et al. (2000) data and measured count rates inside $r_{\rm v}$. Doing so reduces the ROSAT determined flux by ${\simeq}
30\%$ in the case of both RX J1120.1 and RX J1334.3, i.e. so that they are ${\simeq} 10\%$ lower than the values quoted in Table 5. In the case of RX J1701.3, the measured ROSAT flux actually increased slightly when measured inside $r_{\rm v}$ rather than r80. On closer examination of the ROSAT data, it was apparent that the nearby QSO and cluster (see Fig. 7) were contaminating both the source and background apertures in the PSPC image. When the source and background were accumulated only from the (source free) region to the West of the cluster, the ROSAT flux within $r_{\rm v}$ dropped to within ${\simeq} 10\%$ of the XMM-Newton determined value. These examples demonstrate the importance of using XMM-Newton to re-calculate fluxes and luminosities for high redshift clusters detected at low signal to noise by ROSAT. We stress that, despite these flux uncertainties, the SHARC and 160SD catalogs are still fair representations of the underlying cluster population and can still be used to probe the Lx-Tx relation. It is also important to note that the these flux uncertainties do not apply to the ROSAT observations used by Markevitch (1998) to determine Lx values for the 30 low redshift clusters in his Lx-Tx analysis. Those observations have exquisite signal to noise, so that the Lxvalues derived from them will be limited only by the absolute calibration of the instrument.

The fact that we find a Lx-Txrelation consistent with the low redshift value argues against a bias in our selection towards the more luminous tail of the cluster population; for either

Nevertheless the reader might be cautioned that should there remain a common bias for all cluster samples, despite these cross-checks, it is possible the clusters selected represent the brightest portion of the intrinsic high-z sample, and thus remain on the high luminosity end of the true population, hence mimicking the luminosity evolution we see. Further validation of the ROSAT data with CHANDRA and XMM-Newton should help to close this issue.

5.2.3 Comparison with Markevitch (1998) and Vikhlinin et al. (2002)

We come to very much the same conclusions regarding Lx-Tx evolution as did Vikhlinin et al. (2002). Given that these conclusions provide evidence for significant evolution, contrary to previous analyses, it is important to investigate whether this concordance is genuine or coincidental. Both studies rely on Markevitch (1998) for the low redshift comparison, so let us first consider how our XMM-Newton clusters might appear either too luminous, or too cool, compared to the Markevitch (1998) sample:

Markevitch (1998) measured cluster fluxes by summing the count rate inside circular apertures. However, rather than using a cluster specific virial radius, he used a fixed metric aperture of r=1 h50-1 Mpc and ignored any flux that might lie outside this radius. However, this would be a very small effect; Vikhlinin et al. (2002) estimate that the no more than 4% of the total flux would be missed by adopting a fixed r=1 h50-1 Mpc aperture.
When Markevitch (1998) excised the central r=50  h50-1 kpc region from all the clusters in his low redshift sample before measuring the flux, he applied a multiplication factor of 1.06 to account for flux falling inside the excised region. We carried out an excision and flux lost correction on our eight clusters with a cluster-specific surface brightness model (Sect. 5.2.1, and Table 5), we found that for six of the eight clusters, the measured Lx values dropped slightly. The difference in our techniques could systematically affect cluster brightness, compared with the low redshift counterparts, but only at the few percent level.
Uncertainties in the cross calibration between XMM-Newton and the instruments used for the Markevitch (1998) analysis (ROSAT and ASCA) may mimic evolution in the Lx-Tx relation. At this time it is not possible to rule out that possibility; to date, a full comparison of XMM-Newton and ASCA determined cluster temperatures has not been carried out.
Ikebe et al. (2002) have used a different approach to Markevitch (1998) to analyse ASCA observations of low redshift clusters. They measure temperatures that are on average lower than those measured by Markevitch (1998), with the trend becoming more pronounced as Tx increases. However, this should not impact our conclusions regarding Lx-Tx evolution, given that we used a similar technique to Markevitch (1998) to measure Tx.
Some estimate of the effect of external systematic error could be taken by considering the value of A in fitting to a different low-z Lx-Txsample. We have done so using the Arnaud & Evrard (1999) relation and find $A=0.95 \pm 0.2$ (again in the Einstein de Sitter cosmology). Therefore, despite a possible contamination in that sample from cooling flows we still see evidence for evolution.
We note that the first two concerns listed above do not apply to the Vikhlinin et al. (2002) study, as they deliberately adopted the same methodology as Markevitch (1998) and the last three concerns are not likely to be important either; Vikhlinin et al. (2002) have carried out a cross comparison of Chandra and ASCA Tx values for $\simeq
$20 low redshift clusters and find them to be in agreement at better than the 5% level, with no systematic offset. Despite this, we should still consider the possibility that both studies (ours and Vikhlinin et al. 2002) might be pointing to the wrong value of A. At present we cannot quantify how factors such as uncertainties in the cross-calibration of Chandra and XMM-Newton, especially at low energies, might impact our conclusions. We note in particular that Vikhlinin et al. (2002) applied an empirical factor of 0.93 to improve the cross-calibration of the ACIS front- and back-illuminated CCDs below 1.8 keV. Other differences between the observatories may also be important. For example, it will be easier to account for point source and cooling core contamination using Chandra data, because of the improved spatial resolution (see for example the case of RX J1701.3). Alternatively, it should be easier to correctly account for the various particle (cosmic ray and proton) backgrounds using XMM-Newton data, because its CCD's cover a larger area both inside and outside the FOV. In summary, although we cannot rule out the possibility that a combination of factors have lead us to measure an artificially large value for the A parameter, we believe that our XMM measurements provide clear evidence for evolution in the Lx-Tx relation.

We thank M. Arnaud for very useful comments concerning the analysis and interpretation of the data. We also thank M. Markevitch and A. Vikhlinin for their assistance. M. Watson is acknowledged for allowing us to use data from the GT observation of IRAS 13224-3809. The EPIC instrument team are thanked for their continued help in improving the calibration knowledge and help in understanding various effects related to the instrument background. We thank the referee, F. Castander for the careful review and comments that helped to improve the interpretations. D.J.B. acknowledges the support of NASA contract NAS8-39073 (CXC). A.K.R. and R.C.N. acknowledge support from the NASA-LTSA program, contract NAG5-11634 and the hospitality of the Durham University Physics department during the summer of 2002. These data were obtained from observations made by the XMM-Newton Observatory. This is an ESA science mission with instruments and contributions directly funded by ESA Member States and the USA (NASA). The XMM-Newton $\Omega $ Project GT data were provided from a combination of data-rights holders, including the CDS Strasbourg (XMM Science Survey Centre), CEA-Saclay (XMM EPIC PI Team) and the XMM-Newton SOC. AKR acknowledges the Aspen Centre for Physics for hospitality in 2003.


Online Material

Table 5: Summary of cluster parameters for EdS model H0=50, q0=0.5. Spectral fitting errors, L, $\beta $ and $r_{\rm c}$ are 1$\sigma $ on one parameter. Fluxes are the measured, absorbed fluxes in ROSAT band. Fractr is the fractional correction made from the spectral extraction region, to the total counts within the $r_{\rm v}$ radius after accounting for point source excision, inter-chip gaps loss etc.. The 3rd & 4th rows summarise the data for the case when the core region has been excised from the cluster core. (Except for RX J1701.3, where following the CHANDRA data, we excise a cooling flow enhancement to 120 kpc).

Appendix A: Excised point sources

Table A.1: The identifications (XMM-Newton informal ID, nominal RA and Dec ) of point sources that were excised from the spectral and imaging analysis.

Copyright ESO 2004