A&A 423, 49-56 (2004)
DOI: 10.1051/0004-6361:20047170
J. de Plaa1,2 - J. S. Kaastra1 - T. Tamura1,3 - E. Pointecouteau4 - M. Mendez1 - J. R. Peterson5
1 - SRON National Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands
2 - Astronomical Institute, Utrecht University, PO Box 80000, 3508 TA Utrecht, The Netherlands
3 - Institute of Space and Astronautical Science, JAXA, 3-1-1 Yoshinodai, Sagamihara, Kanagawa 229-8510, Japan
4 - CEA/DSM/DAPNIA Saclay, Service d'Astrophysique, L'Orme des Merisiers, Bât. 709, 91191 Gif-sur-Yvette, France
5 - KIPAC, Stanford University, PO Box 90450, MS 29, Stanford, CA 94039, USA
Received 30 January 2004 / Accepted 7 May 2004
Abstract
We report the results from a spatially resolved spectroscopy study with XMM-Newton on the relaxed cluster of galaxies Abell 478.
From the EPIC data we extract a temperature profile and radial abundance profiles for Ne, Mg, Si, S, Ca, Fe and Ni. The abundance profiles
follow the same trends as observed in other clusters.
The spectra of the core of the cluster can be best fitted with a multi-temperature model. We argue that this
multi-temperature behavior is mostly due to projection effects, because of the strong temperature gradient in the core.
Contributions from other effects, for example, intrinsic temperature stratification cannot be fully excluded.
For the first time we measure an underabundance of oxygen in the Galactic absorption component
toward a cluster. The measured oxygen abundance in this absorber is about 0.5 times the solar oxygen abundance as determined by Anders & Grevesse (1989).
Key words: galaxies: clusters: general - galaxies: clusters: individual: Abell 478 - galaxies: cooling flows - galaxies: intergalactic medium - X-rays: galaxies: clusters - ISM: abundances
The evolution and hydrodynamical structure of the hot diffuse X-ray emitting gas in clusters of galaxies is still not well understood. Several clusters show excess X-ray emission in their cores. The so-called cooling flow models (for a review see Fabian 1994) propose radiative cooling of gas in the core to explain the observed surface brightness and temperature profiles. The pressure decrease in the core due to the cooling gas causes a net inflow toward the center of the cluster, hence the name cooling flow. XMM-Newton observations show that the amount of cool gas in the core of cooling-flow clusters is much lower than predicted (Tamura et al. 2001b; Peterson et al. 2003; Kaastra et al. 2001,2004; Peterson et al. 2001). Therefore, the standard cooling-flow model needs to be adjusted. Several mechanisms to explain the observations have been proposed (see Peterson et al. 2001, for a list of possible explanations).
The Abell 478 cluster of galaxies is a good example of a highly relaxed cluster. Earlier X-ray observations with EXOSAT (Edge & Stewart 1991), Ginga and Einstein (Johnstone et al. 1992), ROSAT (White et al. 1994; Allen et al. 1993), ASCA (White 2000; Markevitch et al. 1998) and Chandra (Sun et al. 2003) show a significant excess of X-ray emission in the core, suggesting the presence of a massive cooling flow. The observations performed by missions with sufficient spatial resolution (e.g. ROSAT, ASCA and Chandra) are consistent with a radial temperature decrease toward the core. In a recent Chandra observation an X-ray cavity was discovered in the core which seems to be associated with a lobed radio source (Sun et al. 2003).
In this paper we study the properties of Abell 478 using high-resolution and spatially-resolved spectra obtained with the European Photon Imaging Camera (EPIC, Turner et al. 2001) and the Reflection Grating Spectrometer (RGS, den Herder et al. 2001) aboard XMM-Newton (Jansen et al. 2001). We focus on the physics and metal abundances in the core of Abell 478. In an associated paper Pointecouteau et al. (2004) describe the large scale temperature, gas and dark matter distribution in this cluster using a pointed and an offset XMM-Newton EPIC observation.
Throughout this paper we use H0 = 50 km s-1 Mpc-1 and q0 = 0.5. Using this cosmology 1
is 157 kpc at
the cluster redshift of 0.0881.
The observation of Abell 478 was performed as part of the Guaranteed Time program on February 15, 2002 and had a total duration of 126 ks. Both EPIC MOS cameras were operated in Full Frame mode and the EPIC pn camera in Extended Full Frame mode. For all EPIC cameras the thin filter was used.
For the analysis we use the event lists produced by the XMM-Newton SOC and the 5.4.1 version of the XMM Science Analysis System (SAS).
To correct for enhanced and variable background due to soft proton flares we use the events with energies larger than 10 keV from
pn and CCD 9 of RGS. The data are filtered using upper and lower count-rate thresholds. The
thresholds are calculated in a way analogous to the method used in Pratt & Arnaud (2002). We first make a histogram of the
light curve above 10 keV with a binsize of 100 s and pattern equal to 0, which we fit it with a Poissonian function. The count rate
thresholds are subsequently fixed to
,
where N is equal to the mean number of counts within a bin.
Because the pn is most sensitive to solar flares, we merge the Good Time Intervals (GTI) derived from the pn data also with the
preliminary GTI for MOS. Therefore pn and MOS have roughly the same exposure time. The maximum count rates and the resulting useful
exposure times are shown in Table 1. This simple variable background procedure aims at obtaining maximum statistics
in the core region, where the background is relatively less important, but excludes contributions from the large solar flares.
Because the cluster signal in the 0-4
region dominates over the background, it is not necessary to perform a double
background subtraction nor a multi-band flare search like described in Pointecouteau et al. (2004).
The countrate thresholds we chose are well within the advised threshold of
<1.0 cts s-1 (PN) provided by the XMM-Newton SOC (Kirsch 2003).
Table 1: Maximum count rates and resulting useful exposure times after screening for solar flares.
Because the cluster is slightly elongated along the north-east to south-west axis, we extract
the source and background spectra from elliptical annuli allowing patterns
12 for MOS and
4 for pn.
We use a ratio of 1.2 between the major and minor axis of the ellipse, which we obtain from an empirical 2D
-model fit to the EPIC image. Although the data show some minor deviations from this 2D model, the value is in
the range of 1.2-1.4 (ROSAT, White et al. 1994) and compatible with 1.22 derived by Pointecouteau et al. (2004) who uses the
same method and data.
Further on in this paper we characterize the extraction region by its semi major axis.
We extract the EPIC background spectra from the blank sky event files provided by Read & Ponman (2003).
Relatively bright point sources in the EPIC field of view have been excluded both from the Abell 478
and background datasets.
To gain statistics in each spatial bin we extract the spectra from annuli with a width of more than 30
.
This way, we are less sensitive to the energy dependent shape of the PSF and therefore we neglect this effect
in the rest of the analysis. The downside of choosing large binsizes is that mixing of temperatures within the
bin also increases with binsize if there are strong temperature gradients present. This effect may become
important during multi-temperature fitting.
Although we neglect the influence of the energy dependent part of the PSF, a substantial broadening effect
in the core remains. To correct the normalization of each annulus for PSF effects we calculate correction
factors by fitting the surface brightness profile from the Advanced Camera Imaging System (ACIS) on board
Chandra (Sun et al. 2003) with a double King profile. We normalize both the ACIS fit
and the EPIC surface brightness profile and calculate the ratios between the two for each bin. These
factors range from 0.73 in the core up to 1.03 at 4
.
Then, we multiply the effective area of each annulus with this factor and obtain a PSF corrected response.
Table 2: Systematic errors included in the EPIC spectral fits. (Adapted from Kaastra et al. 2004.)
For the spectral analysis we use the SPEX package (Kaastra et al. 2003) and fit the spectra over the 0.4-10 keV range. We include systematic errors due to uncertainties in the calibration and the background in the spectral fit. The applied systematic errors are shown in Table 2; they have been adapted from Kaastra et al. (2004). These systematic errors enclose possible variations of the Cosmic X-ray Background (CXB), background normalization and errors in the calibration.
After a preliminary analysis, the pn appears to have a gain problem, resulting in a significantly lower value for the
redshift of the source than measured by the MOS instrument.
From MOS we obtain a redshift of
,
which is consistent with the optical value of
(Zabludoff et al. 1990). But, our fits of the pn spectra result in a redshift of
.
To correct for this, we measure the centroid energies of the aluminum, nickel,
copper and zinc background lines. A linear fit (
)
through the measured and expected energies for these lines yields
and
.
After correcting
the energy values in the original eventfile according to this relation, we fit the resulting spectra again.
The corrected spectra show large discrepancies at low energies. Therefore, we conclude that the observed gain problem is
probably not linear, in line with the argumentation of Pointecouteau et al. (2004).
Unfortunately, the number and distribution of background lines is insufficient to derive the exact shape of the gain
correction. For the rest of the analysis we therefore use the original event file, fit MOS and pn
separately and let the redshift free in the fit.
We extract the RGS spectra with SAS version 5.4.1 following the same method as described in Tamura et al. (2001a).
We select the events from a rectangular area on the CCD strip with a full width of 1
in the cross-dispersion
direction and centered on the core of the source. We extract the background spectra from several blank field observations (Tamura & den Herder 2003).
Because the observation was performed before the cooling of the RGS instruments, the spectra contain a relatively
high number of bad columns due to warm pixels. These bad columns appear as absorption-like features in the spectra.
Unfortunately, also CCD gaps are present in the wavelength ranges where line
emission is expected. Together with the relatively high thermal Bremsstrahlung component due to the high temperature
of the cluster, the lines are difficult to resolve. Therefore, we also include the second order spectrum in the fit to
increase the statistics. All RGS spectra are fitted over the 8-35 Å range.
Because the RGS gratings operate without a slit, the resulting spectrum of an extended source is the sum of
all spectra in the (in our case)
field of view, convolved with the PSF
(see for a complete discussion about grating responses Davis 2001).
Extended line-emission appears to be broadened depending on the spatial extent of the source along the dispersion direction.
In order to describe the data properly, the spectral fits need to account for this effect. In practice, this is accomplished by
convolving the spectral models with the surface brightness profile of the source along the dispersion direction (Tamura et al. 2004).
For that purpose we extract the cluster intensity profile from MOS1 along the dispersion direction of RGS,
which we convolve with the RGS response during spectral fitting.
Because the radial profile of an ion can be different from the mean profile, this method is not ideal. We
let the scale of the width and the position of the profile free in the fit to match the profiles of the main
emission lines.
![]() |
Figure 1: Single-temperature fit results. EPIC MOS (stars), pn (diamonds) and RGS (plus-sign) are shown. The dashed line shows an empirical powerlaw fit to the Fe abundance distribution for comparison. |
| Open with DEXTER | |
Table 3:
Averaged EPIC results for Abell 478 for the single temperature (1) and wdem (2) fits.
Errors are given at the 1
confidence level and the <-sign denotes a 2
upper limit.
We leave the abundance of oxygen in the absorption component free to fit the oxygen edge near 0.5 keV.
This parameter is called
.
The abundances are given relative to solar abundances measured by Anders & Grevesse (1989).
We fit the spectra both with a single temperature collisionally ionized plasma model (MEKAL) and a differential emission measure (DEM) model
called wdem. In this particular model the emission measure (Y) is distributed as a function of temperature (T)
as shown in Eq. (1) adapted from Kaastra et al. (2004):
To investigate the role of projection effects, we fit the spectra of the core also with an extra temperature component of 6.5 keV. This way we emulate the spectral contribution of the outer parts of the cluster.
We notice from a preliminary analysis that the oxygen edge near 0.5 keV is not well fitted. Therefore, we use
an absorption model component with variable element abundances and leave the oxygen abundance free. We call
this particular oxygen abundance
further on in the paper.
The results of the EPIC single-temperature fits are presented in Fig. 1 and Table 3.
The upper left plot in Fig. 1 shows that the absorption is not a constant across the cluster.
The temperature profile shows a steep decrease toward the core with a hint for a bend outside the 2.0
radius. The abundances
of neon, sulfur and iron are consistent with an increase toward the center. However, neon and sulfur are also consistent with a flat
distribution, like the other abundances. Because of low statistics we had to discard the
8-10 keV band in the 3-4
annulus.
Figure 2 shows the residuals of the combined MOS spectra of the 0-0.5
region, fitted with a single
temperature component for which the abundances have been set to zero afterward. In this way we show the lines that are
detected in the EPIC spectrum. Except for neon, the abundance of all elements can be fitted
independently. Since the main neon line is close to the Fe-L complex blend at
1 keV, it is difficult to constrain
the neon abundance with EPIC. Because the Fe-K line is the strongest line, this line mainly determines the iron abundance.
The strength of the Fe-L complex is mostly dependent on the temperature structure.
In the core of the cluster the
value from the single-temperature fit is higher than can be expected from
the high signal-to-noise in that area. We therefore also fit the core with the multi-temperature wdem
component defined in Sect. 2.3. These results are also listed in Table 3. The wdem abundances
are all consistent with the single temperature abundances, apart from neon, which is blended by Fe-L. From the
values
the multi-temperature model appears to be preferable over the single-temperature model in the 0-0.5
region.
Then, we add a single temperature component fixed at 6.5 keV to emulate the
projection effect from the hot gas in front of the core. The value of 6.5 keV is roughly the mean temperature outside
the 1
radius. In Table 4 we show the results of this fit. The abundances are again consistent
apart from neon. We notice that the fitted temperatures for the core are now significantly cooler than without the hot component.
The value for 1/
in the wdem component is 0.38
0.07, lower than the value we previously found,
1.12
+0.12-0.3.
Table 4: Averaged results from the single-temperature (3) and wdem (4) fits with an extra 6.5 keV temperature component.
We obtain a reasonable fit to the RGS spectra of the core (0-0.5
)
with the single-temperature
model. The
/d.o.f. is 887/754. Unfortunately, the O VIII line complex is partly within
a CCD bad-column on RGS1 and falls in the dead area of RGS2, which makes the derived oxygen abundance more uncertain.
The Fe-L complex and Ne line near 12 Å are detected and resolved.
Because the EPIC spectra from the core region are best fitted with a multi-temperature model, we also fit the
RGS spectra with the wdem component. This results in a slightly better
/d.o.f. of 867/753. The results
of this fit and the single-temperature model are shown in Table 5 and Fig. 3.
Because the temperature determination with RGS depends mostly on the line emission, the weakness of the lines in
this spectrum causes the temperature to be less well constrained. Furthermore, the energy range of the RGS is
much smaller than EPIC and concentrated on the soft energy part of the spectrum, making it more sensitive
to the cool component of the spectrum. Therefore, the fit is less sensitive to the multi-temperature distribution. Because
of the increased freedom of the fit due to the use of multi-temperature components, the wdem fit is marginally
better than the single-temperature fit. We have also attempted to fit
the RGS spectrum with an extra hot component of 6.5 keV which we added to both a single-temperature and wdem model.
Because the temperatures in this model could not be constrained, we have not pursued this model any further.
![]() |
Figure 2: Residuals of the averaged MOS spectrum of the core of Abell 478 with the abundances put to zero. |
| Open with DEXTER | |
![]() |
Figure 3:
RGS averaged spectrum (0-0.5
|
| Open with DEXTER | |
We analyze the high-resolution spectra of the Abell 478 cluster of galaxies and we derive, for the first time, radial abundance profiles for several elements in this hot cluster. The exposure time of 126 ks allows us to resolve the main emission lines. By fitting these lines and the continuum emission we are able to put constraints on the temperature structure in the core and the abundance distribution of Ne, Mg, Si, S, Ca, Fe and Ni. Furthermore, by fitting the oxygen edge near 0.5 keV we measure an underabundance of oxygen in the absorption component.
![]() |
Figure 4:
Error ellipses for
|
| Open with DEXTER | |
The temperature profile obtained from the single-temperature fits (see Fig. 1 and
Table 3) of the core (<2
)
are
consistent with the temperature profile measured with Chandra (Sun et al. 2003) and ROSAT (White et al. 1994; Allen et al. 1993).
It confirms that Abell 478 has one of the steepest temperature gradients observed so far.
Table 5:
Results from the RGS spectral fits of the core region (0-0.5
)
of Abell 478.
The RGS spectra were fitted using the wdem model component and a single-temperature model.
Errors are given at the 1
confidence level and the <-sign denotes a 2
upper limit.
We left the abundance of oxygen in the absorption component free to fit the oxygen edge near 0.5 keV.
This parameter is called
.
In the core region, within 0.5
,
a multi-temperature (wdem) model fits better than the
single-temperature model. The value of 1/
from RGS is high compared to the cluster samples studied
by e.g. Peterson et al. (2003) and Kaastra et al. (2004). These studies find values for 1/
which are
roughly of the order of 0.5. If we add an extra temperature component of 6.5 keV to account for the emission from the outer
regions in front of the core, then the difference between the single-temperature and wdem component vanishes (see Table 4).
This suggests that the projection effect of the hot cluster material along the line of sight strongly affects the 1/
value
in the fits.
The multi-temperature behavior is not necessarily explained by contamination from the outer parts of the cluster.
Although the X-ray cavities reported by Sun et al. (2003) are not spatially resolved with
XMM-Newton, the value of 1/
and the cool area associated with the X-ray cavity may be related.
Together with the high resolution study of the temperature profile (Pointecouteau et al. 2004), which
reveals a steep temperature gradient in the core region, these effects could also partially account for the multi-temperature
behavior observed in the central bin. Unfortunately, the spatial resolution of XMM-Newton does not allow
us to draw conclusions from these data on the complex structure in the core.
From the single-temperature fits we measure a central increase of the iron abundance from
0.3 at about 4
to
0.4 in the center which is consistent with the mean abundance profile of Abell 478 derived by White (2000) using ASCA.
The central abundance is also consistent with earlier EXOSAT measurements by Edge & Stewart (1991).
The significant central increase of Fe is similar to the profiles of other hot (>6 keV) clusters (Tamura et al. 2004).
Also neon and sulfur seem to follow the same trend in our results. However, due to the large error bars these profiles could still
be consistent with a flat distribution. The central increase of sulfur was also found in the cluster sample of Tamura et al. (2004).
Neon is not resolved in EPIC, because of the lower spectral resolution and blending with lines from the iron L complex.
Therefore, the systematic uncertainty on the abundance of neon is quite large. Moreover, the single-temperature
model tries to fit the Fe-L complex by enhancing the Ne abundance, while the reason for the Fe-L to be enhanced is the
presence of colder material. Because RGS resolves the neon line, its value for the neon abundance is much more robust.
Neon and oxygen are thought to have the same origin. Therefore, we expect their values and spatial behavior
to be comparable, which excludes the very high numbers for neon from the single-temperature EPIC fits.
The nickel abundance we find is likely to be overestimated due to an error in the nickel line
energies (see also Gastaldello & Molendi 2004).
The oxygen abundance could not be constrained in our EPIC fits, but despite the occurrence of bad columns in the line
resolved by the RGS we could get a value for the oxygen abundance consistent with the EPIC upper limits.
We can compare our result with the predicted O/Fe ratios from theoretical models for supernovae Ia and II.
The O/Fe ratio obtained from the RGS wdem fit is
,
which is lower than the average value
of
for hot clusters reported by Tamura et al. (2004), but not significantly different because of the large
spread. The theoretical models predict an O/Fe ratio of <0.05 and 1.5-4 for supernova types Ia and II, respectively.
Our value of
is in between the two predictions, like in the other hot clusters.
In Tables 3 and 5 the values for
and
appear to be
different between EPIC and RGS. This is partly due to the fact that the EPIC/RGS cross-calibration near the oxygen edge
is still not optimal. The EPIC results are also slightly affected by correlations between
,
and the oxygen abundance,
because of blending due to the lower CCD resolution.
The
anti-correlates with absorption, but because fixing the edge results
in large residuals at lower energies, we decided to let the oxygen abundance in the absorption component free.
Figure 4 shows the error ellipses for the oxygen abundance against
.
From the plot it is apparent that
and
are anti-correlated.
If we fix the
abundance to 1.0, the
values we find are consistent with those from
Pointecouteau et al. (2004) and Sun et al. (2003). This plot also shows that
is significantly different from solar. Weisskopf et al. (2004) report from their Chandra observation of the
Crab nebula that, using the abundances from Anders & Grevesse (1989), they find an underabundance of oxygen in the Galactic absorption component
of
.
Takei et al. (2003) find values in the range from
to 0.74
0.14 using a Chandra
observation of Cyg X-2. Within the calibration uncertainties our results are consistent with these numbers.
It is likely that the observed underabundance is an artifact of the solar abundances we use. Allende Prieto et al. (2001)
reported a new solar photospheric abundance of log
dex. This is
0.58
0.08 times smaller than the value of Anders & Grevesse (1989) and in between our values for
.
Pointecouteau et al. (2004) show that the origin of the absorption in the direction of Abell 478 is
most likely Galactic and correlated with infra-red data.
Although
is consistent with the solar oxygen abundance of Allende Prieto et al. (2001)
we cannot fully exclude the existence of a real underabundance of oxygen in the absorbing material along the line of sight.
Finally, we find that the average redshift measured in the pn detector is 0.0775
0.0002, compared to 0.0881
0.0009
measured in the optical by Zabludoff et al. (1990). This is probably a gain-related problem which was not
corrected for during the XMM-Newton SOC data processing. Problems with the redshift determinations are not unique.
Discrepancies between X-ray measured redshifts and optical redshifts have been reported earlier
by Zhang et al. (2004) and Takahashi & Yamashita (2003).
We analyze the high-resolution XMM-Newton spectra of the Abell 478 cluster of galaxies and conclude that:
Acknowledgements
We would like to thank Monique Arnaud for providing useful comments and for the fruitful discussions we had while completing this paper. This work is based on observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA member states and the USA (NASA). The Space Research Organization of The Netherlands (SRON) is supported financially by NWO, the Netherlands Organization for Scientific Research.