F. De Paolis - G. Ingrosso - A. A. Nucita - D. Orlando
Dipartimento di Fisica, Università di Lecce and INFN, Sezione di Lecce, Via Arnesano CP 193, 73100 Lecce, Italy
Received 24 October 2001 / Accepted 21 October 2002
Abstract
Observations with the Extreme UltraViolet Explorer (EUVE)
satellite have shown that the inner region of the Virgo cluster
(centered in M 87 galaxy) has a strong Extreme UltraViolet (EUV)
emission (up to
)
in excess to the low-energy tail
expected from the hot, diffuse IntraCluster Medium (ICM).
Detailed observations of large scale radio emission and upper
limits for hard, non-thermal X-ray emission in the 2 - 10 keV
energy band have been also reported. Here we show that all
available observations can be accounted for by the existence of
two electron Populations (indicated as I and II) in the M 87
Galaxy. The mildly relativistic Population I is responsible for
the EUV excess emission via IC scattering of CBR and starlight
photons. Population II electrons (with higher energy) are instead
responsible for the radio emission through synchrotron mechanism.
The same electrons also give rise to hard non-thermal X-ray
emission (via IC scattering of CBR photons), but the resulting
power is always below the upper bounds placed by present
observations. The non-negligible energy budget of the two electron
populations with respect to that associated with thermal
electrons indicates that the M 87 galaxy is not today in a
quiescent (relaxed) phase. Nuclear activity and merging processes
could have made available this energy budget that today is
released in the form of relativistic electrons.
Key words: galaxies: clusters: individual: Virgo - radio continuum: galaxies - ultraviolet: galaxies
Observations with the EUVE satellite (Bowyer & Malina 1991) have provided evidence that a number of clusters of galaxies produce intense EUV emission, substantially in excess of the low-energy tail expected from the X-ray emission by the hot, diffuse ICM at temperature of a few keV. EUV excesses have been reported for Virgo (Lieu et al. 1996a), Coma (Lieu et al. 1996b), A1795 (Mittaz et al. 1998), A2199 galaxy clusters (Lieu et al. 1999), A4038 (Bowyer et al. 1998) and A4059 (Durret et al. 2002). However, although there is no doubt about the detection of the EUV excess emission for Virgo and Coma clusters, considerable controversy still exists for the other clusters since EUVE results are affected by the variation of the telescope sensitivity over the field of view and depend also on subtraction of background signals that do not come from the cluster (see Berghöfer & Bowyer 2002; Durret et al. 2002).
The initial explanation for the EUV excess emission was that it
is produced by a warm (
K) gas component of the ICM
(Lieu et al. 1996b). For Virgo and Coma clusters,
Bonamente et al. (2001) have proposed that HI
cold clouds might have, at least in part, been released in the
intergalactic space by ram pressure stripping. In this framework,
warm gas can be generated at the interface between the cold phase
(HI clouds) and the hot ICM, e.g. by the mixing layer
mechanism (Fabian 1997). The resulting scenario is then a
three phase gas model, where the hot and warm gas components are
responsible for the bulk of X-ray emission and the soft excess
emission, respectively, while the cold component accounts for
X-ray absorption (Buote 2001).
Models invoking the presence in galaxy clusters of the warm gas
component have problems since they refer to unrelaxed or
turbulent gas conditions: the longevity of such multiphase gas
would then require the presence of a heat source to carefully
balance the radiative cooling, but no obvious heating mechanisms
have been found (Fabian 1996). Indeed, it is difficult for
a warm intracluster gas to remain for a long time at 106 K
since at this temperature the cooling function for a gas with
solar chemical composition has a value
erg cm3 s-1 (Saito & Shigeyama 1999)
implying for the M 87 galaxy warm gas component with metallicity
(Gastaldello & Molendi
2002) and mean number density
cm-3 (Bonamente et al. 2001) a cooling
time of about 150 Myr.
Although strong evidence for cooling flows has been found in
low-resolution X-ray imaging and spectra, high-resolution
observations with the Reflection Grating Spectrometer on the
XMM-Newton satellite show inconsistencies with the simple cooling
flow models. The main problem is the lack of the emission lines
expected from gas cooling below 1-2 keV. Indeed,
observations with both the Hopkins Ultraviolet Telescope (Dixon et al. 1996)
and FUSE (Dixon et al.
2001a,2001b) found no significant far-UV line
emission from gas at 106 K
. Also the observation of a large number of
clusters by the XMM-Newton satellite (Tamura et al.
2001,2002) detected no warm gas component. All
these pieces of evidence seem overwhelming in the sense that a
thermal mechanism for the EUV excess can be ruled out.
Therefore, other mechanisms must be investigated as the source of the EUV excess emission in clusters of galaxies. Indeed, Inverse Compton (IC) scattering of Cosmic Background Radiation (CBR) photons by relativistic electrons present in the ICM was proposed to explain the EUV excess in the case of the Coma Cluster (Hwang 1997; Ensslin & Biermann 1998; Ensslin et al. 1999).
Here, we will examine in detail the case of Virgo cluster, referring in particular to the results of a re-analysis of archival EUVE data for the central region of the cluster performed by Berghöfer et al. (2000a).
In Fig. 2 of the above cited paper, the azimuthally averaged
radial intensity profile of the background-subtracted EUV emission
in the central region of the Virgo cluster (centered in the M 87
galaxy) is given, as well as the contribution of the low-energy
tail of the X-ray emitting hot ICM. Comparison between the two
curves shows the presence of a diffuse EUV excess which extends
up to a radius of
(corresponding to
65 kpc, for an assumed M 87 distance D = 17 Mpc), in the M 87 halo
region
.
The aim of the present paper is to show that IC scattering of CBR
and starlight photons by a population (hereafter denoted as
Population I) of mildly relativistic electrons with energy up to
a few hundred MeV is able to account for the EUV excess emission
observed from the central region of the Virgo cluster, centered
in M 87 galaxy. As we will see in Sect. 2, Population I
electrons have to be described by a power law energy spectrum
with an energy cut-off at
MeV (or slightly
above) in order to avoid both a hard, non thermal, X-ray excess
(due to IC scattering of CBR photons) with respect to the upper
limit of
erg cm-2 s-1 in the
2-10 keV energy band, set by observations (Reynolds et al.
1999) and over-production of radio waves by synchrotron
emission (Owen et al. 2000).
The existence in the M 87 galaxy of a highly relativistic electron
population (which we refer to as Population II) is well
demonstrated by the observations of large-scale radio emission in
the 10 MHz-100 GHz band (Owen et al. 2000).
This radio emission is generally interpreted as synchrotron
radiation produced by relativistic electrons interacting with the
halo magnetic field. To b exact, taking a value
G as typical for the M 87 halo magnetic field, one gets an
electron energy range
1 - 170 GeV. Population II electrons
have to be described with a power law energy spectrum with slope
in order to fit the observed radio power law
spectrum
.
Electrons belonging to the same
Population II, but with energy in the range
0.7 - 1.7 GeV, also
could give rise to a hard, non thermal X-ray emission in the 2 -10 keV region (through IC scattering of CBR photons), but the
resulting power is always negligible with respect to the upper
limit quoted by Reynolds et al. (1999).
Regarding the origin of Population I and II electrons, it is well
known that the main source of non-thermal energy in the M 87 galaxy
is the active nucleus and jet. The current bolometric luminosity
of the inner active region is 1042 erg s-1 -
mainly emitted in the radio band (Owen et al.
2000) - although the kinetic power in the jet (which
matches with the nuclear power for accretion at the Bondi rate)
is
1044 erg s-1 (Di Matteo et al. 2002).
This is believed to be a lower limit for the total power
available for particle acceleration, since the central engine
would have had a cyclic activity with time scale of the order of
100-200 Myr and today be in a quiescent phase (Owen et al. 2000; Corbin et al. 2002).
Therefore, simple lifetime considerations indicate that the
available current energy content is at least
1060 erg.
On the other hand, in the past, the M 87 galaxy could have
undergone a merging which would have made available a further
energy budget (as suggested for the Coma cluster by Blasi
2000). This energy budget could be stored by accelerated
protons that, being long-lived in the intergalactic medium (with
lifetimes of the order of a Hubble time) could be confined in the
galaxy and then released part of their energy to electrons by
Coulomb and/or hadronic interactions (see e.g. Ensslin
2002) or by shock waves in plasma turbulence.
The observed radio luminosity requires that the power emitted by
Population II electrons is 1042 erg s-1. Since the
cooling time of these electrons is dominated by synchrotron
losses with a time scale of
1015 s, one gets a
Population II energy budget of
1057 erg, much less than
the current available energy content. This implies that most of
the total energy made available by the M 87 nucleus could be
stored in Population I electrons. If these electrons are
responsible for the EUV excess emission by IC scattering of CBR
and starlight photons, their power is about 1044 erg s-1, which multiplied by the IC cooling time
1017 s, gives a Population I energy content close to that available.
As far as the electron acceleration mechanism is concerned, it is expected that Population II electrons have been directly accelerated in the active nucleus and jet with high Lorentz factors. Indeed, the two lobes in the observed radio map trace on a large-scale the jet injection, keeping initially the direction of collimated outflows from the core. Further diffusion by magnetic field inhomogeneities may give rise to a diffuse Population II component responsible for the overall radio emission. It is also known that the active nucleus and jet in the M 87 inner region are responsible for inflating bubbles of accelerated cosmic rays which rise through the cooling gas at roughly half the sound speed (Churazov et al. 2001). Shock acceleration at the bubble boundaries could give rise to Population I electrons.
Finally we note that the existence for Population I electrons of
an energy cutoff at
MeV may be justified by
considering that a time evolution for power-law spectrum injected
electrons is expected as a consequence of transport, diffusion
and cooling effects (see Petrosian 2001 for a
detailed analysis for the Coma cluster).
In Sect. 2 we calculate the EUV source function for the IC scattering by relativistic electrons on CBR and starlight photons and the expected EUV excess radial profile towards the M 87 galaxy. In Sect. 3 we present our model results and address some conclusions.
As stated in the Introduction, we assume that Population
I electrons account for
the large-scale EUV excess in M 87 galaxy through
IC scattering of CBR and starlight photons.
As far as the electron spatial distribution is concerned,
despite evidence of granularity in the EUV excess map of M 87
(Fig. 5 in Berghöfer et al. 2000a),
here we adopt for simplicity a continuous radial profile
as representative of an
azimuthally averaged radial distribution of the observed structures.
Therefore, the electron density distribution
as a function of
the electron energy
and the radial coordinate r,
in units of
,
is
In our model the target photons are supplied by the CBR and
Starlight (S) photons, with temperature
K and
K
,
respectively.
According to Eq. (4), Population I
electrons in the energy bands 125-250 MeV and 3-10 MeV
are involved in order to scatter the CBR and starlight photons,
respectively, into the EUV energy range 50-200 eV.
The source function for EUV photon production through IC scattering
is given by (Lang 1974), in units of
, as
As far as the average background photon density
in Eq. (6) is concerned,
in the case of CBR photons it is clearly constant
(
cm-3) within the M 87 galaxy,
whereas for starlight photons it is obtained by the relation
![]() |
(8) |
At this point we are ready to estimate the expected EUV excess radial
profile in the M 87 galaxy.
By integrating the source function given
in Eq. (6)
along the line of sight at impact parameter b(neglecting internal absorption in M 87),
we get the surface brightness profile (in units of
)
![]() |
(13) |
Following the formalism discussed in Sect. 2, we fit the available EUV excess data towards the M 87 galaxy. For this purpose, experimental points (given with the respective error bars in Fig. 1) are derived from those given in Fig. 2 in Berghöfer et al. (2000a), by subtracting the contribution of the low-energy tail due to the X-ray emitting ICM.
Preliminarily, we have fitted the M 87 observed EUV excess by
considering only the IC scattering of Population I electrons on
CBR photons (Model I, with 11 d.o.f.). In this way, we obtain a
poor fit with
by considering all the 15 observed
points (see the dotted line in Fig. 1)
.
![]() |
Figure 1:
The expected photon number count
![]() ![]() |
Open with DEXTER |
Consideration of the IC scattering by the same electron
population also on starlight photons, distributed accordingly to
Eqs. (7)-(11), (Model II, with 9 d.o.f.),
allows us to substantially increase the expected photon number
counts from the innermost part of the galaxy, therefore obtaining
a much better confidence level
,
considering
all the 15 points (see solid line in Fig. 1). This
represents a significant improvement of the fit to all the
observed points since the outcome of the F-test implies that
Model II has to be preferred over Model I within a
confidence level (the F-test value is 3.9).
Selected models with
give a range of values for
and
given by
kpc and
.
As far as the remaining two parameters
and KI- which are actually related by Eq. (3) - our acceptable
fits give
and, correspondingly,
eV-1 cm-3. The obtained
relative ratio between the EUV luminosity from IC scattering on
starlight and CBR photons is in the range of
.
The
parameter
cannot be further constrained by EUV excess
data due to the absence of spectral measurements.
Recent experimental results (Reynolds et al. 1999) set an
upper limit of
erg s-1 cm-2 for the
hard, non-thermal X-ray flux in the 2-10 keV energy band from the
M 87 galaxy. Moreover, radio observations in the 10 MHz-100 GHz
band give a total luminosity
erg s-1 (Herbig & Readhead 1992; Owen et al.
2000). If we do not introduce an energy cutoff for
Population I electrons (see Eq. (2)), we would find both a
hard, non-thermal X-ray excess emission (in the 2-10 keV energy
band, due to IC scattering of CBR photons) and radio wave excess
by synchrotron mechanism, above the upper bound allowed by
observations.
In the following, as stated in the Introduction, we assume that
the radio emission is accounted for by Population II electrons
described by the same distribution law as in Eq. (1) for
Population I electrons, but with different parameter values for
,
,
,
and without energy cutoff
.
Here we recall that the characteristic frequency and the
differential luminosity of the ultrarelativistic synchrotron
radiation are given by (see e.g. Lang 1974)
It is remarkable that the same set of values for Population II
electron parameters that fit radio observations give a hard,
non-thermal X-ray flux by IC scattering of CBR photons (by
electrons in the energy range 0.7-1.7 GeV) in accordance with
the upper bounds given by Reynolds et al. (1999) for the
M 87 galaxy. More recently, the XMM-Newton observatory has
been pointed towards M 87 galaxy with the aim of investigating if
non-thermal IC emission is present in the arms regions and if it
is linked to the power emitted by the radio jet. The upper limit
to the non-thermal X-ray emission has been found to be
erg cm-2 s-1 in the 0.5-8 keV band,
representing less than
of the flux from the thermal
components (Belsole et al. 2001). Our model results
give a hard non-thermal X-ray flux from the inner
region of M 87 galaxy close to the above upper bound.
Clearly, our model results depend on the chosen parameters for
population II electrons which are not, at present, tightly
determined by radio observations. For example, by increasing the
Population II electron core radius (or decreasing the value of
)
we can dilute the X-ray excess (and also the radio
emission) within a wider region, leaving the same total power.
A few comments about the energy budget involved for the two
electron Populations are in order. Observations show that the
current bolometric luminosity of the M 87 inner active region is
1042 erg s-1 - mainly emitted in the radio band
(Owen et al. 2000) - although the kinetic power
in the jet is
1044 erg s-1 (Di Matteo et al.
2002). Simple lifetime considerations then indicate
that the available electron energy content is at least
1060 erg.
However, the total power available for particle
acceleration may be larger since the central engine could have
had a cyclic activity with a time scale of the order of 100-200 Myr
and today be in a quiescent phase (Owen et al.
2000; Corbin et al. 2002). Moreover, in
the past the M 87 galaxy could have undergone a merging phenomenon
which would have made available a further energy budget,
initially stored by relativistic protons and later released
(partially) to electrons by Coulomb and/or hadronic interactions
or by shock waves in plasma turbulence (Ensslin
2002).
Coming back again to our model results, we find that the energy
budget associated with the two electron Populations is about
1060 erg (mostly in Population I electrons), which is always
less or approximately equal to the total energy associated with
thermal electrons.
This result represents an indication that the M 87 galaxy is not
today in a quiescent (relaxed) phase.
However, further more accurate spectroscopic observations in both the EUV and X-ray bands (especially towards the external region of the M 87 galaxy) are necessary in order to confirm the presence of the two electron Populations invoked in the present paper. Moreover, observations towards the other galaxy clusters for which the EUV emission is still uncertain may establish whether the existence of these two electron Populations is a common feature in galaxy clusters or if it is particular to the central, massive, cD galaxies like M 87.
Acknowledgements
We thank the anonymous referee for pointing out some important issues.