Free Access
Issue
A&A
Volume 567, July 2014
Article Number A102
Number of page(s) 10
Section Cosmology (including clusters of galaxies)
DOI https://doi.org/10.1051/0004-6361/201423928
Published online 22 July 2014

© ESO, 2014

1. Introduction

The deep potential wells of clusters of galaxies are permeated by a hot (kT ~ 1–10 keV), optically thin and spatially diffused plasma, the intra-cluster medium (ICM), which is also enriched by heavy elements. It follows that observations in the X-ray band provide a powerful method for studying and classifying these systems. Galaxy clusters can be divided into two main classes, cool-core and non cool-core systems, on the basis of the thermodynamical and chemical properties of their central regions (i.e. within radii 0.15 r500). Cool-core clusters show an almost regular, elliptical morphology on large scales (e.g. Nurgaliev et al. 2013), and their innermost region displays several characteristic properties: the surface brightness and the iron (Fe) abundance are both strongly peaked at the centre (where central cooling times are significantly shorter than the Hubble time) (e.g. Santos et al. 2008; Sanderson & Ponman 2010), the temperature decreases to about one third of the external value (e.g. Hudson et al. 2010; Panagoulia et al. 2014), and the central entropy is 30 keV cm2 (e.g. Cavagnolo et al. 2009; Pratt et al. 2010; McDonald et al. 2013). Conversely, non-cool-core clusters do not show prominent surface-brightness peaks and temperature drops in the cores, and very often, they feature the irregular X-ray morphologies and galaxy distributions indicative of recent merging events.

The study of the spatial distribution of iron in the ICM provides important information that can be used to deepen our understanding of the chemical history and evolution of galaxy clusters in general and of cool-core clusters in particular. Outside the core, the Fe abundance of the ICM is on average one-fourth solar for both cool-core and non-cool-core clusters (De Grandi & Molendi 2001, hereafter DM01; Leccardi et al. 2010; Matsushita 2011), and it shows only a moderate evolution up to a redshift of ~1.3 (Balestra et al. 2007; Maughan et al. 2008; Baldi et al. 2012).

The abundance excess found at the centre of cool core clusters is always associated with a giant, early-type galaxy, the brightest cluster galaxy (BCG), often showing a very extended diffuse stellar envelope. The favoured scenario for the formation of the BCGs is through mergers of galaxies falling in along primordial filaments early in the history of the Universe (e.g. Dubinski 1998; De Lucia & Blaizot 2007; Whiley et al. 2008; Collins et al. 2009; Hung & Ebeling 2012; Ruszkowski & Springel 2009; Lin et al. 2013; Lidman et al. 2013). The presence of a BCG makes this galaxy the prime candidate for producing the peaked abundance profiles (e.g. Fukazawa et al. 2000; Böhringer et al. 2004; De Grandi et al. 2004; Leccardi & Molendi 2008a). De Grandi et al. (2004) estimated that BCGs in low-redshift clusters (i.e. z < 0.1) are able to produce the Fe mass observed in the cluster core during their star formation histories. Furthermore, the observations show that the stellar light profiles of the BCGs are less extended than the Fe abundance profiles, suggesting that there must be a mixing mechanism that moved the metals from their production sites (DM01; Churazov et al. 2003; Rebusco et al. 2005; Graham et al. 2006).

Very little is known about the formation time and subsequent evolution of these peaked Fe abundance profiles in cool-core clusters. Baldi et al. (2012) present the first attempt to determine, with XMM-Newton and spatially resolved data, the evolution of the abundance at different positions in the clusters and with redshift (up to z ~ 1.4) and find no evidence of any evolution in the Fe content in the cores, although with low statistical significance and without distinguishing cool cores from non-cool cores. The major sources of Fe in the Universe are type-Ia supernovae (SNIa), whereas other elements (alpha elements) are produced by core-collapsed SN (SNcc) or by a mix of the two SN types. A measure of the alpha/Fe ratios in the ICM gives a first-order approximation of the relative importance of the SN types to the enrichment (for a review see Werner et al. 2008; Borgani et al. 2008). Early claims of an increasing type-Ia supernovae (SNIa) contribution, with respect to the SN core-collapsed (SNcc) to the ICM enrichment towards cluster centres (Finoguenov et al. 2000; Rasmussen & Ponman 2007), contributed to establishing a scenario where the Fe abundance in cores increased steadily from the time of formation to z = 0 (e.g. Fukazawa et al. 2000; Böhringer et al. 2004). More recent observations from XMM-Newton, Chandra, and Suzaku (e.g. Tamura et al. 2004; Million et al. 2011; Sato et al. 2007; Matsushita et al. 2013), which found a constant ratio between SN types at any radius in cool-core clusters contradicting the early results, requires to reconsider this picture.

In this paper we provide a detailed investigation of iron content and spatial distribution in the X-ray-selected cluster WARPJ1415.1+3612 (hereafter WARPJ1415) at redshift z = 1.03 (Perlman et al. 2002; Fritz et al. 2009), using a deep Chandra observation recently explored by Santos et al. (2012, hereafter S12). This cluster is a high-redshift, strong cool-core cluster (S12, Santos et al. 2010), which is strikingly similar in its thermodynamic characteristics to nearby cool-core clusters. In particular we study the connection between its Fe content and its central galaxy, which was observed with the Subaru and Spitzer telescopes. Our aim is to compare the properties of this distant cluster with those of local (i.e. z < 0.1) cool-core systems to try to constrain the formation history of their iron peaks (DM01; De Grandi et al. 2004). We extended the analysis in DM01 with XMM-Newton data, selecting from the XMM-Newton public archive a small sample of nearby clusters for which we also have near-infrared (NIR) observations of the central BCGs, in order to compare their chemical properties with those of the distant cluster WARPJ1415.

The structure of this paper is as follows. In Sect. 2, we present the X-ray and NIR datasets and the data analysis procedures used to study WARPJ1415 and the nearby cool-core clusters. In Sect. 3, we measure the Fe content in WARPJ1415 and in the local clusters and compare them. We estimate the enrichment times and the supernovae rates needed to produce the observed Fe mass in WARPJ1415. Section 4 is devoted to comparing the BCG’s star-light profiles and the Fe abundance distributions in the ICM in our distant and local clusters and to discussing the ensuing implications for the enrichment histories in cool-core cluster cores. Finally, we summarize our results in Sect. 5.

Throughout this paper we assume the concordance flat ΛCDM cosmology with H0 = 71 km s-1 Mpc-1, Ωm = 0.27, and ΩΛ = 0.73, that for WARPJ1415 at z = 1.03 gives a scale of 8.099 kpc/′′. Iron abundances are measured relative to the solar photospheric values of Asplund et al. (2009), where the solar Fe abundance is ZFe,⊙ = 3.16 × 10-5.

Quoted errors are at the 1σ level, unless otherwise stated.

2. Data and analysis

2.1. WARPJ1415

2.1.1. X-ray data

The galaxy cluster WARPJ1415 was detected in the Wide Angle ROSAT Pointed Survey (Jones et al. 1998; Perlman et al. 2002) and then observed by Chandra ACIS-S (280 ks) and ACIS-I (90 ks). S12 performed a detailed spatially resolved spectral analysis using both Chandra datasets for a total of 370 ks exposure time. Their analysis so far represents the most detailed X-ray analysis of any cool-core cluster at such high redshift, z ~ 1. Details on the data reduction with the CIAO 4.3 software package with CALDB 4.4.5 (Fruscione et al. 2006) are given in Sect. 2 of S12. In S12 the authors extracted ACIS-I and ACIS-S spectra from nine rings at different radii up to ~400 kpc and performed the spectral analysis on the combined spectra (apart for the last ring, 300–400 kpc, where there was no useful signal in the ACIS-I data) by assuming a single-temperature model. The resulting projected and de-projected temperature and Fe abundance profiles are shown in their Fig. 2. (Details on spectral analysis and de-projection technique are in S12.)

The high X-ray count statistic of WARPJ1415 allows the most detailed analysis of the Fe abundance distribution in the core of a cluster at z ~ 1 to date. We performed a new spectral analysis with the XSPEC package (v12.6, Arnaud 1996) starting from the same spectra, response matrices, and the ancillary response matrices computed by S12, but using the apec spectral model instead of the mekal one. We decided to use apec because this model is continuously updated within XSPEC and we want to be consistent with the analysis of the XMM-Newton data of the local clusters (see Sect. 2.2.1). We first checked that our re-analysis agrees with those published in S12. To this aim we combined the ACIS-S and ACIS-I spectra and assumed a single-temperature apec model over the energy range 0.5–7.0 keV and a fixed local absorption equal to the Galactic neutral hydrogen column density (NH) measured at the cluster position (Kalberla et al. 2005) of NH = 1.05 × 1020 cm-2. The gas temperature, metal abundance, and normalization were left unconstrained in the fitting procedure. The resulting projected and de-projected temperature and Fe abundance profiles are consistent within the uncertainties with those published in S12.

We tested whether the Fe abundance might be affected by biases induced by the presence of multi-phase gas in the cluster core. We do not expect the abundance measurement to be strongly affected by the Fe-bias (Buote 2000a,b) or by the inverse Fe-bias (Rasia et al. 2008; Simionescu et al. 2009; Gastaldello et al. 2010) since the measured emission-weighted temperature in the innermost bin is >4 keV, and the Fe L-shell complex (energy at rest ~0.7–1.5 keV) for this cluster at z = 1 only marginally overlaps the energy range considered for the spectral analysis (~0.35–0.7 keV observer frame for a z = 1 cluster). We checked the stability of the Fe abundance measurement by fitting the combined ACIS-S (495 cts) and ACIS-I (94 cts) spectrum of the innermost bin, a circle with a 20 kpc radius, with a two-temperature apec model, where the temperature of the first component is fixed and the temperature of the second component is left free to vary. The values assigned to the first temperature are reported in Table 1. As in S12 we use the Cash statistics, which is preferable for low S/N spectra (Nousek & Shue 1989). From the results of the best-fit shown in Table 1 it is clear that, if present, a second component provides a modest contribution to the total emission and that the value of the measured iron abundance is essentially unaffected by the multi-temperature structure.

Table 1

Best-fitting values for a 2-temperature apec model for the ACIS (S and I) spectrum of the innermost ring.

2.1.2. Near-infrared data

The BCG of WARPJ1415 is a large and massive elliptical galaxy with luminosity in the B band of LB = 3 × 1011LB,. We retrieved the raw broadband images taken with Subaru/Suprime (BVRiz bands; Miyazaki et al. 2002) and Subaru/MOIRCS (J,Ks bands; Suzuki et al. 2008), from the Subaru archive, SMOKA1. MOIRCS provides wide-field imaging in a 4′ × 7′ field of view covered by two 2048 × 2048 arrays with the spatial resolution of 0.117 arcsec/pixel. The radial stellar profile in the NIR Ks band is extracted from MOIRCS data using the IRAF ELLIPSE routine, a procedure that measures best-fit isophotes from the data. The sky-background-subtracted profile shown in Fig. 1 is traced out to a radius of ~80 kpc; however, beyond a radius of ~30 kpc, the profile rapidly dilutes into the background.

thumbnail Fig. 1

Subaru/MOIRCS Ks surface brightness profile of the BCG in WARPJ1415.

Open with DEXTER

Using our own multi-colour catalogue from broadband optical, infrared, and IRAC2 photometry, we estimated the stellar mass of the BCG by fitting the galaxy spectral energy distribution (SED) with LePhare (Arnouts & Ilbert 2011) obtaining M = 9.3 × 1011 M for a Chabrier IMF (Chabrier 2003). With the SED fitting we measured an apparent, rest frame Ks magnitude of 18.69 AB. To estimate the Ks rest frame magnitude, we used all photometric points computed within 1.5′′ radius aperture, then corrected to 2.3′′ (18 kpc rest frame) aperture assuming no colour dependence, with respect to the z band. This step was adopted to avoid BCG flux contamination in IRAC bands from nearby sources. The estimated absolute Ks rest frame magnitude is − 24.73 ± 0.20 AB.

Table 2

The local sample of galaxy clusters collected as a reference sample in this work.

2.2. Local clusters of galaxies

From its thermodynamical properties, the distant cluster WARPJ1415 is very similar to cool-core clusters in the nearby Universe. In this work, our aim is to study the chemical properties of WARPJ1415 by comparing it with those of local systems. To this purpose we selected, from the XMM-Newton public archive, six strong cool-core clusters at redshifts lower than 0.09 that have medium–high total masses (i.e. M500 ≳ 3.5 × 1014 M, Chen et al. 2007), similar to the total mass of WARPJ1415 (M500 = 2.4 × 1014 M, S12). In particular, we included in our sample A3526 (the Centaurus cluster) since its iron abundance profile is very peaked and reaches over-solar Fe abundance in the center, as in the case of WARPJ1415. We point out that the spatial resolution for these low-redshift clusters sample observed with XMM-Newton is very similar to the spatial resolution for a cluster observed with Chandra at z ~ 1. This follows from the fact that the ratio of the angular distances computed at z = 0.05 and z = 1 is similar to the ratio between the point spread functions of the XMM-Newton/EPIC and Chandra/ACIS instruments.

The local clusters are listed in Table 2.

2.2.1. X-ray data

We analysed XMM-Newton archival data for the clusters in Table 2, following the data analysis described in details in our previous work De Grandi & Molendi (2009). We reprocessed the observation data files (ODF) using the Science Analysis System (SAS) version 10.0.0. The effective exposure times of the observations obtained after soft proton cleaning are reported in Table 2. Using the soft-proton, cleaned event files we extracted spectra from annular regions centred on the cluster emission peaks. We subtracted the background using blank-sky fields for EPIC MOS and pn produced by Leccardi & Molendi (2008b). All the spectral fits were performed with the XSPEC package using a one-temperature thermal model (apec) and a two-temperature thermal model (apec+apec). The Galactic hydrogen column is accounted for with a wabs absorption model component assuming a Galactic HI value by Kalberla et al. (2005). Redshifts, temperatures, iron abundances, and normalizations of the two models were left free to vary. Whenever the two-temperature fit provided significant improvements over the single-temperature model (using an F-test), we used the Fe abundance derived from the two-temperature model.

Table 3

Ks band observations parameters of the BCGs in the local clusters.

2.2.2. Near-infrared data

When deriving the NIR surface brightness profiles of the BCGs in the local cool-core clusters, we used Ks band data that were acquired using SofI (Moorwood et al. 1998) mounted on NTT telescope at the ESO/La Silla observatory. Since the targeted galaxies have an angular size ranging from 2 to 4 arcmin, the observations were done with the large field mode, corresponding to a field of view of 5 × 5 arcmin, with a pixel scale of 0.288 arsec/pix.

We performed a standard data reduction of the raw frames and calibration standard star fields using IRAF procedures. The main steps of the processing are dark subtraction, flat-field correction, and sky subtraction. Given the large angular size of the galaxies, sky subtraction is achieved through an on-source/off-source imaging strategy, where one frame targeted the science object and the next observed a blank field. Photometric calibration standards (Persson et al. 1998) were acquired several times during the observation run. The scatter in the zero points ranges from 0.02 to 0.14 mag (see Table 3).

thumbnail Fig. 2

Ks band surface brightness profiles of the BCGs in the local clusters listed in Table 2.

Open with DEXTER

The Ks band background subtracted surface brightness profiles for the BCGs in the local cluster sample are shown in Fig. 2. These Ks light profiles are integrated up to the background level so they encompass faint light contribution, such as haloes or intra-cluster light, up to measured radii. For A3526, the CCD frame border is at 1.4′ (~17 kpc) from the centre of the BCG, and therefore the decreasing of the profile at larger radii is a border effect. For this cluster we limited our analysis to radii up to 17 kpc. In all other cases, we verified that the CCD borders are well outside the range considered in these plots.

We performed a 1D fit of the BCG light profiles shown in Fig. 2 using a Sersic model and integrated the best-fit model up to a radius of 18 kpc to obtain the aperture Ks magnitudes of the galaxies (except for A3526 where we used the magnitude estimated by Hoffer et al. 2012). The absolute Ks magnitudes (AB reference system) are reported in Table 4, and errors on the magnitudes are ±0.10 mag.

2.3. A note on the near-infrared band

To probe the mass distribution of the BCG in detail and to measure possible faint contribution from intra-cluster light, we have based our analysis on deep Ks band data rather than other filter passbands. Indeed, NIR luminosities and, in particular, K-band magnitudes, represent a more direct link with the stellar mass content of galaxies than the visual ones (Kauffmann & Charlot 1998). The light in the visual bands is dominated by the brightest main-sequence stars, whose luminosity decreases as stellar populations age. In contrast, a sort of conspiracy in the post main sequence evolution of stars maintains the K band luminosity of a stellar population roughly constant when it ages (Charlot 1996). The K-band magnitudes vs. X-ray-derived physical quantities are expected to correlate with a smaller scatter with respect to the optical ones, and this will allow us to put stronger constraints on the link between the BCGs and the ICM properties.

The rest frame Ks band at z ~ 1 falls on the 4.5 μm channel of Spitzer/IRAC (Infrared Array Camera). However, IRAC does not have enough spatial resolution to effectively estimate surface brightness profile for the BCG of WARPJ1415 at z = 1.03. We thus applied a K correction to the measured Ks-band data, based on a fit to the galaxy SED, to compare the BCG magnitude with the local sample. We have not considered any size or luminosity evolution for the distant BCG.

The Ks profile of the distant BCG, shown in Fig.1, samples a slightly different rest frame with respect to the local BCGs, but still in the NIR band (i.e. in the J band corresponding to ~1.2 μm), where the contamination from stellar evolution of the brightest main-sequence stars is still negligible.

Table 4

ICM Fe mass computed within overdensity Δ = 2500, Fe mass excess, and magnitude within an aperture of 18 kpc for the distant cluster WARPJ1415 and for the local XMM-Newton cluster sample.

3. Iron mass and iron mass excess

3.1. The iron mass

The iron mass in the ICM (in solar units), enclosed within a sphere of radius r from the cluster centre, is given by the following integral: (1)where AFe is the atomic weight of iron, mH the atomic unit mass, and ZFe = nFe/nH is the de-projected Fe abundance in solar units. (nFe and nH are the iron and hydrogen densities by number.) We estimate the iron mass in WARPJ1415 within r2500 (i.e. kpc, see Table 1 in S12) linearly interpolating the cumulative iron mass profile computed with Eq. (1).

The errors on this mass are dominated by the systematics due to the large errors in the Fe abundance profile of WARPJ1415 for radii 0.2r200 ~ 180 kpc (S12), . To estimate these systematic errors we compute the Fe mass with Eq. (1) twice, fixing to two extreme values: the best-fitting value of WARPJ1415 de-projected iron abundance profile at r> 0.2r200 and the iron abundance outside the core measured in nearby cool-core clusters (DM01; Leccardi et al. 2010). The former value (0.21 ± 0.16) represents a lower limit for the Fe abundance outside the core in WARPJ1415, and it takes the conservative possibility of a positive evolution of from z = 1 to z = 0 into account. The latter value (0.34 ± 0.01, from Leccardi et al. 2010, converted into Asplund et al. 2009 solar units) represents an upper limit to the Fe abundance in the cluster outskirts or, in other words, assumes that at redshift 1 is the same as the one at z < 0.3 (Balestra et al. 2007; Leccardi et al. 2010). We finally compute the mean between these upper and lower Fe masses to derive the Fe mass within r2500: MFe,2500 = (7.57 ± 1.26) × 109 M, with the systematic errors computed as half the difference between the two mass limits.

We de-project the XMM-Newton profiles of the local clusters in Table 2 following the method described in Ettori (2002) and Ghizzardi et al. (2004) and estimate the Fe masses of these clusters within r2500 with Eq. (1). We use the relation r2500 ~ 0.3r200 (e.g. Ettori et al. 2010, see their Appendix) with r200 taken from Reiprich & Böhringer (2002). The Fe masses within r2500 are given in Table 4. From this table and the results in DM01 (for clusters in common between the two samples we take the new XMM-Newton values), we find that the Fe mass distribution for cool-core clusters in the local Universe ranges between ~0.4–6 ×1010 M and that the Fe mass obtained for WARPJ1415 is consistent with the Fe masses measured in nearby clusters.

3.2. The iron mass excess

We use Eq. (1) to also measure the iron mass excess, , in WARPJ1415 and in the local clusters. This mass is computed by substituting in Eq. (1) ZFe with , which is defined as the iron abundance associated with the abundance peak in the cluster core after subtracting the average Fe abundance measured at radii larger than 0.2r200: .

We compute the Fe mass excess and its errors for WARPJ1415 taking the two limits discussed previously for into account (i.e., 0.21 and 0.34), obtaining ±0.40) × 109 M3. This measure is consistent with the Fe mass excess distribution for nearby galaxy clusters: ~1.0–10 × 109 M (see Table 4 and DM01). For WARPJ1415 we find a ratio /MFe of 0.22 ± 0.06, which again agrees within the uncertainties with the local cool-core clusters mean value, 0.20 ± 0.05.

By comparing the Ks absolute magnitude of the BCG in WARPJ1415 with those of the local sample shown in Table 4, we find that their values are comparable. This, and the similarity between the Fe mass excess in this and in local systems (see also De Grandi et al. 2004), strengthens our hypothesis that the iron present in the excess has also been produced in this distant cluster by the stars that form the BCG at z = 1.

3.3. Iron abundance excess formation time scale

The high redshift of WARPJ1415 allows us to put some constraints on the formation time scale of the Fe peak. We use the simple formalization given by Böhringer et al. (2004), (2)where is the luminosity in the B-band of the BCG, SNR is the SN type Ia rate at the redshift of the cluster measured in SNU (1 SNU = 10-12 SNe yr-1), and ηFe is the iron yield per SN. In Eq. (2) the Fe mass excess is assumed to be entirely produced by Type Ia SNe neglecting additional minor contribution such as stellar mass loss or newly formed stars in the BCG. The rate of SNIa has been recently measured to be SNU in galaxy clusters at 0.5 < z < 1.46 (Barbary et al. 2012). This value suggests that there is little or no evolution in the cluster SNe rate from z ~ 1 to the present time. (For a compilation of recent cluster rate measurements see Dilday et al. 2010; Maoz & Mannucci 2012.) With a standard yield of MFe = 0.7 M per SNIa (Iwamoto et al. 1999), we find tenrich = 16.2 Gyr, which is significantly greater than the age of the Universe at the redshift of WARPJ1415 (i.e. tH(z = 1) = 5.9 Gyr). We conclude that current estimates of the Fe production rates at high redshift fail to account for the presence of the high Fe abundance peak at z = 1 in WARPJ1415.

3.4. Constraints on the SN rate in BCGs at z > 1

As shown in the previous section, a constant SNR is inconsistent with our results, it is also inconsistent with other observational evidence and with theoretical work (e.g. Maoz & Mannucci 2012). Maoz et al. (2010) argue that the SN delay time distribution (DTD) can be described as a power law with index close to − 1. To test that our data are consistent with such a model, we estimated the total number of SN that have exploded in the BCG by assuming a SNR of the form: SNR = SNRo·(t/to)-1, where to is the time that has passed between the formation of the BCG and our observation of WARPJ1415; SNRo is the SN rate at to, which we derive by multiplying the stellar mass of our BCG, M, with the SNR per unit mass (SNuM), i.e. SNRo = SNuM·M, where we adopt the SNuM derived by Barbary et al. (2012) for high-redshift red-sequence early type galaxies (see their Table 7); and M is taken from Santos et al. (2012). This leads to a value of SNRo = 0.474 yr-1. The total number of SNIa, SNtot, is obtained by integrating the SNR from a minimum delay time, tmin, to to, which leads to SNtot = SNRoto ln(to/tmin). Finally, the total mass of Fe produced is obtained by multiplying the total expected SN, SNtot, by the yield per SN, ηFe. Other authors (e.g. Greggio 2005) have proposed that at early times the SNR might fall off less rapidly, with a slope of −1/2. We have modelled the SNR in this scenario with a broken power-law model with index −1/2 from tmin to a break time, tb, and − 1 from tb to to. In this case the total number of SN, SNtot, turns out to be SNtot = SNRoto {2[1−(tmin/tb)1/2] + ln(to/tb)}.

Before comparing the Fe masses that can be estimated from different choices of the model parameters, we note that these predicted masses should be compared with the total Fe mass, which is given by the sum of the Fe excess mass observed in X-rays and the Fe mass locked up in stars in the BCG. The latter can be estimated from the total stellar mass assuming a typical metal abundance of ZFe = 1.2, (e.g. Maoz et al. 2010, and refs. therein). By performing this calculation, we find an Fe mass locked up in stars in the BCG of 2.1 × 109 M that, as expected, is close to the Fe excess and a total estimated Fe mass of 3.9 × 109 M, which can now be compared with expectations from the different DTDs.

We then start by considering the case of a DTD shaped like a broken power law with tmin = 0.04 Gyr and tb = 0.4 Gyr and consider the two cases of BCG formation at z = 2 and z = 3 corresponding to to of 2.2 and 3.6 Gyr, respectively. We find a total expected Fe mass of 1.8 × 109 M and 3.0 × 109 M, respectively, but both these values fall short of the measured Fe mass in the case of a late BCG formation, to = 2.2 Gyr, by more than a factor of 2 and in the case of an early formation, to = 3.6, by about 30%. In the case of a power-law-shaped DTD with tmin = 0.04 Gyr, we find a total expected Fe mass of 2.3 × 109 M for to = 2.2 Gyr and 3.8 × 109 M for to = 3.6 Gyr, respectively. For the late BCG formation scenario, we are still below, by about a factor of 2, the observed Fe mass, while for the early BCG formation scenario we are in good agreement.

Summarizing, our calculations suggest that a late BCG formation is unlikely, while an early BCG formation with a relatively flat SNR at early times, the broken power-law model, is below but possibly not inconsistent with the measured value. The best agreement is found with a combination of early BCG formation and steep SNR extending down to very short delay times, namely 0.04 Gyr.

At the end of this section we can draw some conclusions. The galaxy cluster WARPJ1415 at z = 1.03, from its X-ray and NIR properties, is virtually indistinguishable from a strong cool-core cluster at z ≲ 0.1. Indeed, its Fe mass and Fe mass excess, its luminosity, and the stellar mass of its central BCG all agree with what it is observed in local cool-core clusters.

The presence of a well-formed Fe peak at the centre of this distant cool-core cluster, all similar to the ones found in nearby clusters, places strong constraints on the formation time of the peak itself. As shown previously, the enrichment time of the Fe abundance peak is ~4 Gyr, and to produce all the Fe in the central excess and in the stars of the BCG, we require SNIa rates that are much higher than those observed in both low-redshift or high-redshift galaxy clusters.

The most straightforward implication of our results is that the bulk of the enrichment occurred at very early times, within a few 100 Myr of the formation of the BCG, possibly during the BCG assembly process. Indeed, in the case of the power-law DTD, which best reproduces our data, half of the SN explode within 0.4 Gyr of the formation of the BCG. This is consistent with observational evidence showing that the star formation rate (SFR) during this time could have been extremely high. A good example is the so-called Spiderweb galaxy at z = 2 with an observed SFR of 1390 ± 150 M yr-1 (Seymour et al. 2012), which converts into an equally huge SNR (e.g. Calura et al. 2007; Maoz et al. 2012). Other observations show that ~50–60 per cent of the iron released by SNIa explosions occurs within ~1 Gyr of a burst of star formation (e.g. Maoz et al. 2010, 2012, and references therein for the most recent measure of the delay-time distribution of SNe Ia as a function of time). Both observation (e.g. Daddi et al. 2004; Rosati et al. 2009; Collins et al. 2009; Stott et al. 2010) and theory (e.g. Dubinski 1998; Conroy et al. 2007a; De Lucia & Blaizot 2007), point out that BCGs form at early times (z > 2) relatively quickly via hierarchical merging (or cannibalism) of small groups of galaxies falling along the filamentary large-scale structure prior to cluster virialization. The growth of the BCGs continue to be fed by mass flows up to z = 1.5−2, when the stellar population of the BCG start to evolve mostly passively (e.g. Brodwin et al. 2013). Nearly all BCGs in the centres of local cool-core clusters still show little ongoing star formation, with SFRs that correlate with the ICM cooling rate (e.g. Rafferty et al. 2008; O’Dea et al. 2010; Donahue et al. 2010; McDonald et al. 2011); however, these rates contribute to a negligible Fe production. The SFR in WARPJ1415 of ~2–8 M/yr-1 (S12), as inferred from the [OII] emission line in the BCG spectrum, is comparable to typical SFRs measured in BCG in cool-core clusters at low redshift (i.e., 1–10 M/yr-1O’Dea et al. 2010), indicating that the major epoch of metal production in this cluster consistently happened at redshifts higher than 1.

All these observations suggest that the bulk of the Fe abundance excess in WARPJ1415 was likely produced by the prompt release of the SNIa shortly after its BCG assembly at z ≫ 1 and relatively short delay times.

thumbnail Fig. 3

Comparison between the observed (black dots) and expected (red lines) Fe abundance excess profiles for the galaxy clusters in the XMM-Newton local sample.

Open with DEXTER

4. BCG stellar light and Fe abundance profiles

In the previous section we found that the Fe mass excess in WARPJ1415 was most probably produced by the stars forming the BCG at redshifts 1.5–2. Now we want to compare the stellar light distribution of the BCG with the Fe abundance excess profile.

To this aim we compute the Fe abundance profile expected when the Fe distribution traces the light distribution of the BCG. Following DM01, we find that the projected abundance ZFe,proj(bmin,bmax) measured within a bin with bounding radii bmin, bmax is related to the de-projected Fe abundance distribution, ZFe(r), by the following equation: (3)where nH(r) is the hydrogen density (by number), and ZFe(r) = nFe(r) /nH(r) is the abundance excess profile with nFe(r) equal to the Fe density (by number). We assume that nFe(r) follows the stellar light distribution of the BCG, l(r), i.e. nFe(r) ∝ l(r), and then compute the projected Fe abundance profile expected, ZFe,proj, under this hypothesis (ZFe(r) ∝ l(r) /nH(r)) with Eq. (3). We, finally, compare this expected Fe excess profile with the observed one.

DM01 performed a similar comparison on a sample of local clusters observed with BeppoSAX, showing that the observed Fe abundance profile in the Perseus cluster was much broader than the predicted one. This result was later confirmed by several authors working on nearby clusters observed with the XMM-Newton satellite (e.g. Churazov et al. 2003; Rebusco et al. 2005; David & Nulsen 2008). Under our assumption that the iron in the excess was formed by the stars in the BCG, this difference indicated that the Fe produced in the BCG spread out several tens of kpc from the place of origin.

We redid this exercise for the local cool-core clusters listed in Table 2 and for WARPJ1415 at z = 1. The results are plotted in Fig. 3 for the local clusters and in Fig. 4 for WARPJ1415: in each panel the data points are the observed Fe abundance excess, whereas the continuous lines are the expected Fe abundance excess profiles. The profiles are arbitrarily normalized to the innermost bin because in this context, we are only interested in the comparison between the shapes of the profiles. From Fig. 3 it is clear that for all the local clusters, the observed profiles are substantially broader than the expected ones.

In Fig. 4, we plot two expected profiles for WARPJ1415 to include the uncertainties in the Ks band BCG’s profile (see Fig. 1). These uncertainties have been included by computing the maximum and minimum expected Fe abundance excess profiles with Eq. (3) using the Ks band profile twice, the first time adding its 1σ statistical errors to the Ks band profile and the second time, subtracting it (Fig. 4). Since the errors in the cumulative Ks band profiles are not independent, these systematic errors may be overestimated. We have also experimented with different spatial binnings of the X-ray Fe abundance profile and find consistent results. We conclude that, despite the larger scatter produced by the X-ray and NIR measurements in WARPJ1415, its observed Fe abundance excess profile is closer, or at most it differs marginally from the predicted one. This is at variance with what it is observed in low-redshift cool-core clusters, where the Fe abundance excess profiles are always significantly broader than the expected ones (Fig. 3).

Under the reasonable assumption that WARPJ1415 will evolve into a present-day cool-core cluster, the smaller difference between the observed and expected profiles (Fig. 4) may indicate that the processes responsible for the broadening of the Fe excess profile in local systems have not yet had time to produce significant modifications in this distant cluster.

The time available for widening the Fe abundance peak excess is 7.9 Gyr, i.e., the look-back time at z = 1.03. On the basis of this observation alone, it is not possible to know whether there is a dominant mechanism that broadens the Fe profile or whether this widening is the result of the combined action of more mechanisms. Neither it is possible to hypothesize whether the Fe broadening is the result of many smaller events during a long period of time or if it is the result of a single event happening in a relatively short time period. What is apparent is that the time available to broaden the Fe peak is very long.

The mechanisms that could mix the ICM and produce the observed Fe profiles are several (see Schindler & Diaferio 2008, and references therein). One possible mechanism is turbulent mixing of the ICM due to AGN activity (e.g. Rebusco et al. 2005; Simionescu et al. 2009, 2010; Gitti et al. 2011; O’Sullivan et al. 2011; McNamara et al. 2014). Kirkpatrick et al. (2011) found that clusters where the BCGs have experienced recent AGN activity (observed in the form of bright radio emission, cavities, and shock fronts embedded in the ICM) have heavy elements distributed anisotropically and aligned with the large-scale radio and cavity axes, indicating that they are being transported from the halo of the BCG into the ICM along the large-scale outflows driven by the radio jets. The amount of transported gas is substantial, and it is consistent with the results of simulations showing that AGN outflows are able to advect ambient, iron-rich material from the core to much larger radii of hundreds of kpc (Gaspari et al. 2011b,a; Planelles et al. 2014).

thumbnail Fig. 4

Comparison between the observed (black dots) and expected (lines) Fe abundance excess profiles for WARPJ1415 at z = 1.03. The red and blue lines are the expected profiles computed taking the errors in the observed BCG’s Ks profile shown in Fig. 1 into account.

Open with DEXTER

Another mechanism at work could be the ICM-sloshing induced by minor mergers in cool-core clusters (Ascasibar & Markevitch 2006; Roediger et al. 2011), a phenomenon that is now more frequently discovered in nearby cool-core clusters (Simionescu et al. 2010; Ghizzardi et al. 2010). Indeed, in Ghizzardi et al. (2013) the observations of the central regions of A496 show that the gas-sloshing is able to shift the metals up to several tens of kpc out from the BCG (see also O’Sullivan et al. 2014).

Although still debated, the stellar mass of the BCGs appears to increase by a factor ~2–3 from redshift 1 to 0 (Lin et al. 2013; Laporte et al. 2013; Lidman et al. 2012, see however Collins et al. 2009; Stott et al. 2010). This mass accretion occurs through mergers due to fast and close encounters of small groups of galaxies onto the BCGs (e.g. Lidman et al. 2013; Tonini et al. 2012). The fate of the merging galaxy stars, which are collisionless, is to accrete onto the central giant galaxy increasing the stellar mass and the intra-cluster light surrounding such a galaxy (e.g. Conroy et al. 2007b; Murante et al. 2007; Henriques & Thomas 2010). Indeed, once the streams of stars have been formed, they begin to decay since they are disrupted by the tidal field of the cluster within timescales that are 1.5 times their dynamical time in the cluster (tdyn ≲ 1 Gyr, Rudick et al. 2009). Less obvious is the fate of the enriched ISM and of the infalling galaxies, which behave as a collisional fluid. This gas is probably tidally or ram-pressure stripped (Kapferer et al. 2007) from its disrupted galaxy and mixed with the ICM at a greater distance from the BCG, with respect to the stars, in the end possibly contributing to the widening of the Fe peak between redshift 1 and 0. The relative importance of this process should be investigated in depth both observationally and with numerical simulations (see Vazza et al. 2010).

5. Summary

Spatially resolved studies of the iron content in high-redshift clusters are very important because they provide much needed information on the enrichment process in these systems. Unfortunately, thus far, only a handful of high redshift cluster have been observed at the necessary sensitivity level, amongst them, the X-ray luminous cool-core cluster WARPJ1415 at z = 1.03, which is the subject of this paper. Our main findings can be summarized as follows.

  • 1.

    By comparing the central Fe mass excess from WARPJ1415with those of local cool-core systems, we find that theFe abundance excess produced by the BCG wasalready in place at z = 1.

  • 2.

    The large amount of metals observed at z = 1 implies that the enrichment process must have occurred early on. More specifically, since SN rates at z = 1 would imply enrichment timescales that are longer than the age of the Universe, our measures point to an early and intense period of star formation most likely associated with the formation of the BCG.

  • 3.

    The Fe mass excess in WARPJ1415 agrees with expectations from the combination of early BCG formation (z ~ 3) and steeply increasing SNR in the past extending down to very short, i.e. 0.04 Gyr, delay times.

  • 4.

    While for local cool-core clusters the Fe distribution is broader than the NIR light distribution of the BCG, in WARPJ1415 the two distributions are consistent, indicating that the process responsible for broadening the Fe distribution in local systems has not yet started in this distant cluster.


1

Based in part on data collected at Subaru Telescope and obtained from the SMOKA, Subaru-Mitaka-Okayama-Kiso Archive, which is operated by the Astronomy Data Centre, National Astronomical Observatory of Japan (Baba et al. 2002).

2

This work is based in part on observations made with the Spitzer Space Telescope obtained from the NASA/IPAC Infrared Science Archive, both of which are operated by the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration.

3

In agreement with S12, who find a Fe mass excess of ×109 M, the small difference is due both to our slightly different spectral analysis and to the integration routine used to calculate the Fe mass excess.

Acknowledgments

The authors thank S. Borgani, S. Ghizzardi, and F. Mannucci for useful comments that helped to improve this work. S.D., P.T., S.M., M.R., and P.R. acknowledge the support by INAF PRIN 2012 (P.I. S. Molendi), A.F. acknowledges the support by INAF PRIN 2010 (P.I. L. Guzzo). J.S. has received funding from the European Programme FP7/2007-2013 under grant agreement 267251 “Astronomy Fellowships in Italy” (Astrofit). This research has made use of the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, of the NASA’s High Energy Astrophysics Science Archive Research Center (HEASARC), and, of the XMM-Newton archive.

References

All Tables

Table 1

Best-fitting values for a 2-temperature apec model for the ACIS (S and I) spectrum of the innermost ring.

Table 2

The local sample of galaxy clusters collected as a reference sample in this work.

Table 3

Ks band observations parameters of the BCGs in the local clusters.

Table 4

ICM Fe mass computed within overdensity Δ = 2500, Fe mass excess, and magnitude within an aperture of 18 kpc for the distant cluster WARPJ1415 and for the local XMM-Newton cluster sample.

All Figures

thumbnail Fig. 1

Subaru/MOIRCS Ks surface brightness profile of the BCG in WARPJ1415.

Open with DEXTER
In the text
thumbnail Fig. 2

Ks band surface brightness profiles of the BCGs in the local clusters listed in Table 2.

Open with DEXTER
In the text
thumbnail Fig. 3

Comparison between the observed (black dots) and expected (red lines) Fe abundance excess profiles for the galaxy clusters in the XMM-Newton local sample.

Open with DEXTER
In the text
thumbnail Fig. 4

Comparison between the observed (black dots) and expected (lines) Fe abundance excess profiles for WARPJ1415 at z = 1.03. The red and blue lines are the expected profiles computed taking the errors in the observed BCG’s Ks profile shown in Fig. 1 into account.

Open with DEXTER
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.