A uniform metal distribution in the intergalactic medium of the Perseus cluster of galaxies
Norbert Werner1,2, Ondrej Urban1,2,3, Aurora Simionescu1,2,4, Steven W. Allen1,2,3
1Kavli
Institute for Particle Astrophysics and Cosmology, Stanford University, 452 Lomita Mall,
Stanford, CA 94305-4085, USA
2Department of Physics, Stanford University, 382 Via Pueblo Mall, Stanford, CA 94305-4060,
USA
3SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA 94025, USA
4Institute of Space and Astronautical Science (ISAS), JAXA, 3-1-1 Yoshinodai, Chuo-ku,
Sagamihara, Kanagawa 252-5210, Japan
Most of the metals (elements heavier than helium) ever produced by stars in the member
galaxies of galaxy clusters currently reside within the hot, X-ray emitting intra-cluster gas.
Observations of X-ray line emission from this intergalactic medium have suggested a
relatively small cluster-to-cluster scatter outside of the cluster centers1,2 and enrichment
with iron out to large radii3,4,5, leading to the idea that the metal enrichment occurred early
in the history of the Universe3. Models with early enrichment predict a uniform metal
distribution at large radii in clusters, while late-time enrichment, favored by some previous
studies6,7, is expected to introduce significant spatial variations of the metallicity. To
discriminate clearly between these competing models, it is essential to test for potential
inhomogeneities by measuring the abundances out to large radii along multiple directions
in clusters, which has not hitherto been done. Here we report a remarkably uniform
measured iron abundance, as a function of radius and azimuth, that is statistically
consistent with a constant value of ZFe = 0.306 ± 0.012 Solar out to the edge of the nearby
Perseus Cluster. This homogeneous distribution requires that most of the metal enrichment
of the intergalactic medium occurred before the cluster formed, likely over 10 billion years
ago, during the period of maximal star formation and black hole activity.
Between 2009 and 2011, we obtained a total of 84 observations of the Perseus Cluster with the
Suzaku X-ray satellite, as a Key Project for that mission. The pointings covered eight azimuthal
directions from the cluster center out to an offset angle of 2o, with a total exposure time of over 1
million seconds. We have analyzed the data from all three functioning X-ray Imaging
Spectrometers, extracting spectra from annuli centered on the cluster center. We modeled the
spectra from each of the 76 independent regions as a single-temperature thermal plasma in
collisional ionization equilibrium, with the temperature, iron abundance, and spectrum
normalization included as free parameters8.
The measured radial and azimuthal variation of the iron abundance of the hot intra-cluster
medium (ICM) out to the edge of the Perseus Cluster along the eight different directions is
presented in Fig. 1. We define this `edge’ as r200, the radius within which the mean enclosed mass
density of the cluster is 200 times the critical density of the Universe at the cluster redshift (for
the Perseus Cluster r200 = 1.8 Mpc, corresponding to 82 arcmin4). We have tested the statistical
uniformity of the measured iron abundance at radii r > 400 kpc (i.e. beyond the central
metallicity peak associated with the brightest cluster galaxy4,9) by modeling the data from all
radii and azimuths with a constant abundance. The measured chi-square value of 65.8 for 75
degrees of freedom is consistent with the null hypothesis of a constant iron abundance with a
value of ZFe = 0.306±0.012 Solar (68% confidence limit; we adopt a Solar iron abundance of
3.24×10−5 relative to hydrogen by number10). This abundance is also consistent with the mean
value measured at intermediate radii (~0.2-0.5r200) for a sample of 48 clusters previously
observed with XMM-Newton2. The iron abundance values measured at large radii in the less
massive Virgo Cluster are lower, although these are likely biased by the multi-temperature
structure of the gas in that system5. The azimuthally averaged iron abundance profile for the
Perseus Cluster is shown in Fig. 2.
Enrichment scenarios in which the ICM is predominantly enriched with metals after the
formation of clusters predict a non-uniform metal distribution, which is inconsistent with the
observations presented here. Ram-pressure stripping11 would result in a negative gradient in the
metal abundance profiles out to large radii, as well as significant azimuthal variations, with
higher metal abundances expected along directions connecting to surrounding large scale
structure filaments12. Furthermore, metals ejected by galaxies at later times should roughly
follow the distribution of galaxies, producing an approximately constant metal mass to light
ratio, which is also inconsistent with observations13. Although large scale sloshing motions14 may
be able to mix metals to some degree, the strong gradients in the entropy distributions of cluster
atmospheres15, including Perseus4,8, make the ICM convectively stable, prohibiting the efficient
mixing of metals that are initially non-uniformly distributed across large radial ranges.
The observed uniform iron abundance distribution at large radii in the Perseus Cluster
requires early enrichment of the intergalactic gas in the proto-cluster environment. This
enrichment was most likely driven primarily by galactic winds16 - energetic outflows of metal
enriched gas - which are expected to be strongest at epochs around the peak of star formation17
and active galactic nuclei (AGN) activity18 (redshifts z ~ 2-3, or lookback times of 10-12 billion
years). The combined energy of supernova explosions16 and AGN19 must have been strong
enough to expel most of the metals from the galaxies at early times, and enrich and mix the
intergalactic gas. This gas was later accreted by clusters and virialized (increasing its entropy
through shock heating) to form the present ICM.
The total iron mass in the ICM of the Perseus Cluster, calculated from the measured iron
abundance and ICM density profiles8, is about 50 billion Solar masses, with about 60% residing
beyond 0.5r200. The dominant fraction of this iron was likely supplied by type Ia supernovae
(SNIa), which are thought to have produced 60-90% of iron in the Perseus Cluster13 depending
on the assumed supernova yields. Based on our measurements, we estimate that at least 40
billion SNIa contributed to the chemical enrichment of the proto-cluster environment that later
formed the Perseus Cluster. For such a large iron mass to be expelled by the winds from the
galaxies at early times, a significant fraction of these SNIa must have exploded shortly after the
epoch of peak star-formation. This is consistent with recent findings based on SNIa delay time
distributions20,21, which imply that a large fraction of SNIa explode less than ~5×108 years after
the formation of the progenitor binary system (prompt SNIa). SNIa with longer delay times will
continue contributing to the enrichment of the ICM after clusters virialize, and are expected to be
partly responsible for the metallicity peaks surrounding the brightest cluster galaxies in the
centers of many clusters22.
A unique prediction of the early enrichment scenario is that essentially all galaxy clusters
with masses comparable to the Perseus Cluster should have homogeneous iron abundance
distributions of about one-third Solar at large radii. Also, contrary to some initial findings6,7,
there should be no substantial redshift evolution in the ICM metallicity outside of the central
regions of clusters, out to z ~ 2. The presence and strength of such evolution are a matter of
ongoing debate2,23,24. The observed large iron abundance of the high-entropy gas in the outskirts
of the Perseus Cluster is also consistent with the idea that the highest-energy cosmic rays (above
a few 1018 eV) may be primarily iron nuclei25 accelerated by cluster formation shocks26.
Additionally, because much of the metal rich ICM seen at large radii has been accreted from the
surrounding large scale structure filaments, our model predicts that the tenuous warm-hot
intergalactic medium (WHIM) permeating the cosmic web, in which up to half of the baryons in
the Universe currently reside27,28, is likely to be substantially enriched in metals. The same is true
for the hot circumgalactic gas accreted by galaxies from the intergalactic-medium.
References:
1. Matsushita, K. Radial profiles of Fe abundance in the intracluster medium of nearby clusters
observed with XMM-Newton. Astron. Astrophys. 527, 134-146 (2011)
2. Leccardi, A. & Molendi, S. Radial metallicity profiles for a large sample of galaxy clusters
observed with XMM-Newton. Astron. Astrophys. 487, 461-466 (2008)
3. Fujita, Y., et al. High Metallicity of the X-Ray Gas Up to the Virial Radius of a Binary Cluster
of Galaxies: Evidence of Galactic Superwinds at High-Redshift Publ. Astron. Soc. Jpn. 60,
343-349 (2008)
4. Simionescu, A., et al. Baryons at the Edge of the X-ray–Brightest Galaxy Cluster. Sci. 331,
1576-1579 (2011)
5. Urban, O., et al. X-ray spectroscopy of the Virgo Cluster out to the virial radius. Mon. Not. R.
Astron. Soc. 414, 2101-2111 (2011)
6. Balestra, I., et al. Tracing the evolution in the iron content of the intra-cluster medium. Astron.
Astrophys. 462, 429-442 (2007)
7. Maughan, B.J., Jones, C., Forman, W. & Van Speybroeck, L. Images, Structural Properties,
and Metal Abundances of Galaxy Clusters Observed with Chandra ACIS-I at 0.1 < z < 1.3.
Astrophys. J. 174, 117-135 (2008)
8. Urban, O., et al. Azimuthally Resolved X-Ray Spectroscopy of the Perseus Cluster Out to Its
Edge. Mon. Not. R. Astron. Soc. submitted [arXiv: 1307.3592] (2013)
9. De Grandi, S., Ettori, S., Longhetti, M. & Molendi, S. On the iron content in rich nearby
clusters of galaxies. Astron. Astrophys. 419, 7-18 (2004)
10. Feldman, U. Elemental abundances in the upper solar atmosphere. Physica Scripta 46,
202-220 (1992)
11. Gunn, J. E. & Gott, J. R. III On the Infall of Matter Into Clusters of Galaxies and Some
Effects on Their Evolution. Astrophys. J. 176, 1-19 (1972)
12. Domainko, W., et al. Enrichment of the ICM of galaxy clusters due to ram-pressure stripping.
Astron. Astrophys. 452, 795-802 (2006)
13. Matsushita, K., Sakuma, E., Sasaki, T., Sato, K. & Simionescu, A. Metal-mass-to-light ratios
of the Perseus Cluster out to the virial radius. Astrophys. J. 764, 147-156 (2013)
14. Simionescu, A., et al. Large-scale Motions in the Perseus Galaxy Cluster. Astrophys. J. 757,
182-187 (2012)
15. Voit, G. M., Kay, S. T. & Bryan G. L. The baseline intracluster entropy profile from
gravitational structure formation. Mon. Not. R. Astron. Soc. 364, 909-916 (2005)
16. De Young, D.S. On the origin and evolution of iron-enriched gas in clusters of galaxies.
Astrophys. J. 223, 47-55 (1978)
17. Madau, P., et al. High-redshift galaxies in the Hubble Deep Field: colour selection and star
formation history to z ~ 4. Mon. Not. R. Astron. Soc. 283, 1388-1404 (1996)
18. Brandt, W. N. & Hasinger, G. Deep Extragalactic X-Ray Surveys, Annu. Rev. Astron.
Astrophys. 43, 827-859 (2005)
19. Fabjan, D., et al. Simulating the effect of active galactic nuclei feedback on the metal
enrichment of galaxy clusters, Mon. Not. R. Astron. Soc. 401, 1670-1690 (2010)
20. Mannucci, F., Della Valle, M. & Panagia, N. Two populations of progenitors for Type Ia
supernovae? Mon. Not. R. Astron. Soc. 370, 773-783 (2006)
21. Maoz, D., Mannucci, F. & Brandt, T. D. The delay-time distribution of Type Ia supernovae
from Sloan II. Mon. Not. R. Astron. Soc. 426, 3282-3294 (2012)
22. Böhringer, H., Matsushita, K., Churazov, E., Finoguenov, A. & Ikebe, Y. Implications of the
central metal abundance peak in cooling core clusters of galaxies. Astron. Astrophys. 416, L21L25 (2004)
23. Baldi, A., et al. An XMM-Newton spatially-resolved study of metal abundance evolution in
distant galaxy clusters. Astron. Astrophys. 537, 142-153
24. Andreon, S. The enrichment history of the intracluster medium: a Bayesian approach. Astron.
Astrophys. 546, 6-17 (2012)
25. Abraham, J., et al. Measurement of the Depth of Maximum of Extensive Air Showers above
1018  eV. Phys. Rev. Lett. 104, 091101 (2010)
26. Norman, C. A., Melrose, D. B. & Achterberg, A., The Origin of Cosmic Rays above 1018.5
eV. Astrophys. J. 454, 60-68 (1995)
27. Cen, R. & Ostriker, J. P. Where Are the Baryons? Astrophys. J. 514, 1-6 (1999)
28. Davé, R., et al. Baryons in the Warm-Hot Intergalactic Medium. Astrophys. J. 552, 473-483
(2001)
29. Anders, E. & Grevesse, N. Abundances of the elements - Meteoritic and solar. Geochim.
Cosmochim. Acta 53, 197-214 (1989)
Acknowledgements The authors are grateful for discussions with other members of the Perseus
Cluster Suzaku Key project collaboration, as well as with Yu Lu, Paul Simeon, and Roger
Blandford. This work was supported by the Suzaku grants NNX09AV64G, NNX10AR48G,
NASA ADAP grant NNX12AE05G, and by the U.S. Department of Energy under contract
number DE-AC02-76SF00515.
Author Contributions NW led the writing of the manuscript. OU reduced and analyzed the
data. AS and NW contributed to the data analysis. SWA contributed to the writing of the
manuscript and is the principal investigator of the Suzaku Key Project data. All authors discussed
all results, developed the interpretation, and commented on the manuscript.
Competing Interests The authors declare that they have no competing financial interests.
Correspondence Correspondence and requests for materials should be addressed to NW (email:
[email protected]).
Figure 1: Iron abundance profiles expressed in Solar units10 measured along eight different
directions in the Perseus Cluster. At the temperature of the Perseus Cluster, the iron abundance
measurement is driven by the Fe–K lines. Compared to the Fe-L lines, these lines are
significantly less prone to systematic biases associated with the presence of multi-temperature
gas. The spectral extraction regions span a relatively small range in radius and azimuth, where
the plasma can be well approximated as single-temperature. There is no evidence for a bias in the
temperature measurements due to gas clumping8. We have verified that spectral fits with the FeL line complex (0.8-1.5 keV energy range) excluded give consistent results to those in the full
0.7-7 keV band. Collisional ionization and electron-ion equilibrium are likely to hold out to r200
in the Perseus Cluster, where the equilibration timescales are only ~3.5×108 yr and ~7.3×108 yr,
respectively. Within a radius r ~ 400 kpc (~0.2r200; not shown here) metal enrichment is strongly
influenced by the brightest cluster galaxy and the metallicity is centrally peaked4,9. Along the
eastern direction, large scale sloshing motions are present and are uplifting the ICM from smaller
radii out to r ~ 650 kpc14. Thus we also do not present results for r < 650 kpc in this direction.
Beyond these radii, the iron abundance measurements are consistent with a radially and
azimuthally constant value of ZFe = 0.306 ± 0.012 Solar, indicated by the dotted lines. Expressed
in units of an older, and historically more commonly used Solar iron abundance value29, our best
fit constant iron abundance is ZFe = 0.212 ± 0.008 Solar. The plotted error bars are 68%
confidence intervals.
Figure 2: Azimuthally averaged iron abundance profile for the Perseus Cluster. The profile
was determined by simultaneous fitting of the spectra at a given radius for all eight azimuthal
directions. While the iron abundance was assumed to be the same in all directions in the fits, the
temperatures and spectral normalizations were allowed to vary independently along the different
azimuths. The dotted line shows the best fit constant iron abundance value of ZFe = 0.306 Solar.
The plotted error bars are 68% confidence intervals. Our iron abundance values are lower than
the value measured in the 650-1100 kpc range with XMM-Newton13; however, the XMMNewton measurement suffers from significant systematic uncertainties due to the high and
variable particle background of that instrument compared to Suzaku.