A&A 365, L74-L79 (2001)
D. M. Neumann1 -
M. Arnaud1 -
R. Gastaud2 -
N. Aghanim3 -
D. Lumb4 -
U. G. Briel5 -
W. T. Vestrand6 -
G. C. Stewart7 -
S. Molendi8 -
J. P. D. Mittaz9
Send offprint request: D. M. Neumann
1 - CEA/DSM/DAPNIA Saclay, Service d'Astrophysique,
L'Orme des Merisiers, Bât. 709, 91191 Gif-sur-Yvette, France
2 -
CEA/DSM/DAPNIA Saclay, SEI, 91191 Gif-sur-Yvette, France
3 -
IAS-CNRS, Universié Paris Sud, Bât. 121, 91405 Orsay Cedex, France
4 -
ESTEC, European Space & Technology Centre, Keplerlaan 1, Postbus 1,
2200 AG Noordwijk, The Netherlands
5 -
Max-Planck Institut für extraterrestrische Physik,
Giessenbachstr.,
85740 Garching,
Germany
6 -
NIS-2, MS D436,
Los Alamos National Laboratory,
Los Alamos, NM 87454, USA
7 -
Dept. of Physics and Astronomy, Leicester University,
Leicester LE1 7RH, UK
8 -
IFC/CNR, Via Bassini 15, 20133 Milano, Italy
9 -
University College London,
Mullard Space Science Laboratory,
Holmbury St. Mary,
Nr. Dorking,
Surrey,
RH5 6NT,
UK
Received 2 October 2000 / Accepted 17 October 2000
Abstract
We present here the first analysis of the XMM-Newton EPIC-MOS data of the
galaxy group around NGC 4839, which lies at a projected distance to
the Coma cluster center of 1.65-1
Mpc. In
our analysis, which includes imaging, spectro-imaging and spectroscopy
we find compelling evidence for the sub group being on its first infall onto
the Coma cluster. The complex temperature structure around NGC 4839 is
consistent with simulations of galaxies falling into a cluster environment. We
see indications of a bow shock and of ram
pressure stripping around NGC 4839. Furthermore our data reveal a
displacement between NGC 4839 and the center of the hot gas in the group of
about 3005-1 kpc. With a simple approximation we can explain this
displacement by the pressure force originating from the infall, which acts
much stronger on the group gas than on the galaxies.
Key words: galaxies: clusters: general - intergalactic medium, general - Cosmology: observations - large-scale structure of the Universe - X-rays: general
Author for correspondance: ddon@cea.fr
There is now compelling evidence from studies in different wavelength bands
that the Coma cluster has substructure. The search for substructure in the
Coma cluster starts more than a decade ago: Fitchett & Webster
(1987) find a double structure in the cluster center, and Mellier
et al. (1988) find a
second peak in the large scale galaxy distribution around NGC 4839, which
was further confirmed in a study by Merritt & Tremblay (1994).
While these
studies were based on optical data, Briel et al. (1992)
using X-ray data coming from the all-sky-survey of the ROSAT satellite, found
substructures coinciding with NGC 4911, and again, NGC 4839. Another
study by White et al. (1993) based on a deeper ROSAT
observation
revealed even more structure. In this paper we concentrate on the XMM-Newton
(Jansen et al. 2001)
observation of the sub group around NGC 4839, which has a projected distance
from the Coma cluster of 1.65-1 Mpc. In Sect. 2 we review briefly the
observation and data treatment of the EPIC data (see Turner et al.
2001).
This is followed by sections dealing
with imaging
and spectro-imaging. In Sect. 5 we perform a spectroscopic analysis in
different regions which is followed by discussion and conclusion.
This paper is one of three papers which present XMM-Newton data on the
Coma cluster (see also Arnaud et al. 2001a; Briel et al.
2001).
Throughout this paper, we assume H0=50 km s-1 Mpc-1 5-1 and
q0 = 1/2. At a redshift of z=0.0232,
corresponds to 38.9 kpc.
The subgroup associated with NGC 4839 was observed for 36 ksec with XMM-Newton in full frame mode with the medium filter. Due to remaining problems with the CTI (Charge Transfer Inefficiency) calibration in the extended full frame mode of the PN camera we concentrate in this paper on the data from the EPIC-MOS1 and EPIC-MOS2 cameras. The background intensity of the instruments varies considerably with time over several orders of magnitude. In order to optimize the signal-to-noise of the observation we reject time intervals with high background count rates. For this selection we binned the counts in the energy band 10-12 keV in time intervals of 100 s. In this energy band we hardly detect any astrophysical sources, due to the very small effective area of the mirrors. We reject time bins with more than 15 counts per 100 s. After this selection we obtain exposure times with 30.6 ksec for EPIC-MOS1 and 30.9 ksec for EPIC-MOS2.
Figures 1 and 2 display the images obtained in the band 0.5-2.0 keV and 2.0-5.0 keV, respectively. The band 5.0-10.0 keV is dominated by the particle background. As already observed with ROSAT (Briel et al. 1992; White et al. 1993; Dow & White 1995) we see a tail-like structure around NGC 4839 pointing away from the Coma cluster. The orientation of NGC 4839 in the optical coincides with the orientation of the tail. The X-ray emission of NGC 4839 seems somewhat disconnected to the diffuse X-ray emission coming from the group in the center of the pointing.
![]() |
Figure 1:
The 0.5-2.0 keV image of the EPIC-MOS1 and EPIC-MOS2
camera. To correct for
vignetting we corrected each individual photon with a weighting function
(see Arnaud et al. 2001b).
The image is gauss filtered with
![]() |
Open with DEXTER |
![]() |
Figure 2:
The 2.0-5.0 keV image of the EPIC-MOS1 and EPIC-MOS2 camera. To
correct for
vignetting we corrected each individual photon with a weighting function
(see Arnaud et al. 2001b).
The image is gauss filtered with
![]() |
Open with DEXTER |
In order to look for spectral variations we calculated a hardness ratio image of the data (see Fig. 3). To account for the background, especially the particle background, we used another Coma pointing with the central coordinates of RA = 13h 01m 37s Dec = 27d 19m 52s (J2000), which is situated South-East of the Coma cluster at roughly the same offset angle than the sub-group. (46 arcmin in comparison to 51 arcmin for the subgroup pointing). We orientate this background pointing in such a way that it accounts for the emission of Coma itself plus stray light in the pointing for NGC 4839. However, this rotation has no measurable effect on the hardness ratio image, the features with or without this rotation rest identical. This indicates that the background we measure is dominated by particles.
As we are not interested here in the emission of different point sources but only in spectral differences of the extended X-ray emission we use a wavelet algorithm (the MVM algorithm Rué & Bijaoui 1997) to filter out point-like structures for the hardness ratio images. Our hardness ratio image is calculated as hr = (image(2-5 keV) - image(0.5-2.0 keV))/(image(2-5 keV) + image(0.5-2.0 keV)).
We do not take into account any vignetting effects in the hardness ratio image since we know from calibrations that the vignetting functions are parallel to each other for energies below 5 keV.
The resulting image is shown in Fig. 3.
We can see that there are indeed spectral variations across the field-of-view.
Around NGC 4839 we see that the tail-like structure has a different colour
(is softer) than the regions North or South surrounding it
(the tail has a mean hardness ratio of
).
The region, which shows a harder spectrum in the South
(
)
seems also to have a
tail-like structure, which is parallel to the tail seen in the original image.
The core of NGC 4839 shows
.
The appearance of a region with hard spectrum in the South-East or North, at the border of the hardness ratio image is not statistically significant. We cannot state with certainty that these regions are hotter due to an interaction of the group with the intracluster medium (herafter ICM) of the Coma cluster. These ``hot'' regions are not unlikely to be an artifact which is created by low photon statistics together with remaining stray light.
![]() |
Figure 3: The background subtracted hardness ratio image in the bands 2.-5.0 keV and 0.5-2.0 keV. The image shows the data only up to a radius of 12 arcmin. Regions outside this area do not show sufficient statistics to give meaningful results. The contours are the results from a wavelet filtering of an image in the band 0.5-2.0 keV |
Open with DEXTER |
For our spectral analysis we selected different regions in the pointing, which show spectral differences in the hardness ratio image. In our spectral analysis we restrict ourselves to the energy range between 0.3-10.0 keV. The summary of the spectral fitting is given in Table 1. The position and extraction regions are displayed in Table 2. The spectra for some of the regions are shown in Fig. 5. For the spectral fitting we used XSPEC (Arnaud 1996) and the MEKAL code (Mewe et al. 1985; Mewe et al. 1986; Kaastra 1992; Liedahl et al. 1995). The galactic absorption was fixed to the value of the 21 cm line ( nH=9 1019 cm-2; Dickey & Lockman 1990). For the background spectrum determination we used the observations performed on the Coma background pointing (see above). In order to take into account instrumental variations across the field-of-view we selected for the background the same regions on the detector than for the source spectra. As the particle contribution is somewhat detector dependent, it is better to take the same regions on the detector rather than the regions corresponding to the Coma emission (the orientated regions, like for the hardness ratio image), since our background is highly dominated by particles and not by cluster emission/stray light. In order to estimate the influence of background variations we also used the observations of the Lockman hole as background. For EPIC-MOS1 we used the Lockman hole observations with the thin filter, which is not entirely consistent, since all Coma pointings were observed with the medium filter. However, this will give us an estimate on the intrinsic variations of the background. The comparison of the best fit temperature estimates are shown in Fig. 4. As one can see, the temperatures are very similar.
![]() |
Figure 4: The temperature estimates of the four regions with different background observations. The horizontal axis: Coma background, the vertical axis: Lockman hole background. Error bars are 90% confidence level |
Open with DEXTER |
To account for vignetting we used the approach described in Arnaud et al. (2001b).
![]() |
Figure 5: The spectra of region 2 (left) and region 3 (right). The black crosses show EPIC-MOS1 data, the red crosses show EPIC-MOS2 data |
Open with DEXTER |
Region | description | kT | metal | ![]() |
![]() |
Abund. | |||||
1 | core of | 3.12.43.7 | 0.560.280.96 | 0.82 | 1.781.891.69 |
NGC 4839 | |||||
2 | tail of | 4.84.06.0 | 0.700.311.6 | 0.99 | |
NGC 4839 | |||||
3 | region south | 8.05.613. | 0.010.02.8 | 0.69 | 1.511.431.62 |
of tail | |||||
4 | main body of | 4.44.04.7 | 0.360.260.50 | 0.87 | |
sub group |
region | center RA | center Dec. | radius |
(h:m:s) | (deg:arcmin:arcsec) | in arcsec | |
1 | 12:57:24 | 27:29:54 | 40 |
2 | 12:57:17 | 27:29:37 | 53 |
12:57:08 | 27:29:06 | 65 | |
3 | 12:57:23 | 27:27:50 | 68 |
12:57:15 | 27:26:49 | 63 | |
4 | 12:56:53 | 27:24:27 | 205 |
The origin of this X-ray emission is somewhat unclear. It might be the AGN
of NGC 4839, which is known as a radio source, or the emission might be thermal.
Fitting
a power law with the results as shown in Table 1 gives a reduced
in the order of 1.3. Fitting a thermal plasma improves the fit with
a reduced
of about 0.8-0.9.
This region, for which we only adopted a thermal emission model shows high temperature and metallicity. The temperature is too high for a galaxy, which ranges more about 1 keV and not around 4.5 keV, as observed. This might be explained by ram pressure stripping: the hot gas is blown out as the medium encounters the ICM of the Coma cluster. Due to the high velocity of the infall of NGC 4839 into the Coma cluster, shocks could occur, which heat the intergalactic medium. The fact that we see the tail pointing away from the center of the Coma cluster is a strong indication for this scenario (see also discussion).
As can be seen in Fig. 3, the region South of NGC 4839 seems to have a different spectrum with respect to center or tail of NGC 4839. We fit a thermal plasma and a power law to the spectrum of this region. The results, again are shown in Table 1. Both models for the emission are acceptable. The intergalactic medium around NGC 4839 could be heated by shocks caused from the infall of NGC 4839 onto the Coma cluster (see also the discussion section).
A power law emission could be linked to the particles responsible for the radio emission observed in NGC 4839. NGC 4839 is a radio source showing two opposite jets and a complex radio tail. This radio tail is extended in S-W direction and coincides nicely with region number 3 (see Venturi et al. 1990). If the X-ray emission in this region is linked to the radio tail it could come from Inverse Compton scattering of cosmic microwave background photons as they interact with the relativistic electrons originating from the radio jet. However, in this case it is not clear whether XMM has intrinsically enough sensitivity to resolve this IC emission.
Comparing the surface brightness of region 2 and region 3 we find a difference
in intensity of a factor of about 3.5 in the energy band 0.5-2.0 keV. If we
assume that the emission is thermal and thus the radiated energy
and furthermore that regions 2 and 3 are in pressure equilibrium we expect a
difference in surface brightness of the order of 2.5, since region 2 has about
kT=5 keV and region 3 has most likely a temperature of about kT=8 keV. The
remaining difference between the surface brightness intensities might, in this
case be due to differences in metallicities or extend in the line-of-sight.
Region 2 has a high metallicity and is in a temperature range, where
line emission is more dominant than at 8 keV.
Taking this all together (see also discussion) we think that it is more likely that the emission of region 3 is due to thermal bremsstrahlung rather than due to IC scattering.
This is the center of the extended X-ray emission of the subgroup. The fitted temperature of about kT=4.4 keV is consistent with ASCA results by Watanabe et al. (1999). The metallicity of this group seems to be typical with a best fit abundance of 0.3. The abundance seems to be significantly smaller than the one observed in the tail of NGC 4839 (region 2).
Another feature which is striking in our observations is the displacement
of the extended X-ray emission belonging to the sub group of NGC 4839 and NGC
4839 itself of about 8 arcmin (see for example Fig. 3). This
corresponds to roughly 300
kpc.
From the optical study by CD96
NGC 4839 seems to be right in the centroid of the galaxy distribution
belonging to the
sub group. It seems therefore likely that the galaxies belonging to the
sub group are on average closer to the Coma center than the hot gas belonging
to the sub group. Such a displacement between galaxies and gas is expected when
the sub group encounters the gas of the main cluster. The
ICM of the main cluster slows down the infall. As the galaxies in the
sub group show larger densities and therefore a smaller surface-to-mass ratio
than the gas, the galaxies are less affected by the ram pressure force than the
intra group gas. In the following we make a first order calculation
to
see whether this effect is large enough to explain the displacement of about
300
h50-1 kpc between NGC 4839 and the center of the gas. For a rough
calculation we estimate that the acceleration due to gravitational forces is
much stronger than the counteracting ram pressure force. We furthermore assume
a total mass of the Coma cluster of
(Hughes 1989; Briel et al. 1992; CD96),
and that the Coma cluster mass is a factor of ten higher than the mass of the
infalling group (CD96). In a very rough approximation we suppose that the Coma
mass
distribution is to first order point-like. This gives an infall velocity
of NGC 4839 of 2400 kms-1, supposing a distance to the Coma
cluster to r=1.6 Mpc (this is the projected distance see also
Dow & White 1995, or
CD96 - the most likely infall direction is close to the plane-of-sky). The
approximated infall velocity is not too
different to the calculated value of about
kms-1 by
CD96.
For the ram pressure we use
and for the density
distribution of
the ICM a beta-model (
). For
simplicity
we assume
and
kpc, which is in good
agreement with the ROSAT measurements by Briel et al. (1992).
For
we use their calculated value of
3 10-3h50-1/2
cm-3 and transform it into
.
For the surface of the gas encountering the ICM of Coma we suppose a circular
surface with a radius of
300h50-1 kpc (we take the radius of the
outermost contour in Fig. 3) and a gas mass of about
.
The force exerted by the ram pressure is
.
S
is
the
surface of the sub group gas, m its mass and a the deceleration.
For the slowing down we calculate
.
We suppose
,
which is
certainly true at
distances greater 1.6 Mpc. Filling in all parameters we get
kms-1, which means that the actual infall velocity is
diminished by 540 kms-1 for the hot gas of the sub group. Extrapolating a
difference in velocity between NGC 4839 and
the gas of 540 kms-1 (we assume that the ram pressure force is negligible for
the galaxy) to the time which the group needed to go from the virial radius
(about 3 Mpc) to the observed distance of 1.6 Mpc we see that we obtain a
difference in location of about 400 kpc assuming a constant infall velocity of
2000 kms-1 for NGC 4839. This value which was determined with many
approximations is of the same order of magnitude as the observed displacement
between gas and galaxy. The differences in location
between NGC 4839 and the group gas is thus likely to be due to the ram
pressure
exerced by the pressure of Coma's ICM. In addition,
more detailed modeling of
this infall is certainly very interesting, since it is very likely that
the displacement can be used to determine the density
or temperature of the ICM at large radii.
We obtain a temperature of the sub group of about 4.5 keV. Comparing this to
the central temperature of the Coma cluster of about 8-9 keV and applying
the M-T relation of
we obtain a mass ratio of
.
This ratio is very high in
comparison
to the
result of CD96, who find that the sub group has 5-10% the mass of the Coma
cluster. This discrepancy can be either explained by superposition of the
gas in the group with the hotter ICM of the Coma cluster or
the fact that the ICM of the
NGC 4839 group is heated up above its initial virial temperature. A plausible
scenario is that the subgroup is in a process of internal merging. - The
somewhat disrupted morphology of the extended emission seen
in Fig. 1
and Fig. 2 is an indication for that.
Such processes of internal merging
were already observed in hydrodynamic simulations (see for
example Schindler & Müller 1993). If the ICM of the sub group
were heated
up by a factor of two, which is reasonable, and its initial virial temperature
was 2 keV, we would obtain a mass ratio of
,
which is in good agreement with the values obtained by
CD96.
We performed a first analysis of the XMM-Newton data of the NGC 4839 group, which has a projected distance of the Coma cluster of 1.6 h50-1 Mpc. In our analysis we resolve for the first time a complex temperature structure around NGC 4839 which can be explained by a first infall of the galaxy onto the Coma cluster. Furthermore, we measure a displacement between NGC 4839, the dominant galaxy of the group and the extended X-ray emission coming from its ICM. In a simple approximation we can explain this displacement by a ram pressure force of Coma's ICM acting on the hot gas of the sub group as it falls onto the Coma center.
Acknowledgements
We would like to thank J. Ballet for support concerning the SAS software and J.-L. Sauvageot for providing the gain correction. We want to thank L. Feretti and A. Bykov for very useful discussions.