Iron abundance distribution in the hot gas of merging galaxy clusters

TL;DR

This study uses XMM-Newton observations to analyze the distribution of iron and thermodynamic properties in merging galaxy clusters, revealing insights into their chemical enrichment and the effects of mergers on intracluster medium properties.

Contribution

It provides a comprehensive analysis of the chemical and thermodynamic properties of the ICM in merging clusters, highlighting the persistence of metal-rich cool cores and supporting early enrichment scenarios.

Findings

Fe abundance profiles are similar to non-cool-core clusters at large radii.

Remnants of cool cores can survive major mergers.

No clear correlation between Fe abundance and pseudo-entropy at outer radii.

Abstract

We present XMM-Newton/EPIC observations of six merging galaxy clusters and study the distributions of their temperature, iron (Fe) abundance and pseudo-entropy along the merging axis. For the first time, we focus simultaneously, and in a comprehensive way, on the chemical and thermodynamic properties of the freshly collided intracluster medium (ICM). The Fe distribution of these clusters along the merging axis is found to be in good agreement with the azimuthally-averaged Fe abundance profile in typical non-cool-core clusters out to $r_{500}$. In addition to showing a moderate central abundance peak, though less pronounced than in relaxed systems, the Fe abundance flattens at large radii towards $\sim$0.2-0.3 $Z_\odot$. Although this shallow metal distribution is in line with the idea that disturbed, non-cool-core clusters originate from the merging of relaxed, cool-core clusters, we find that in some cases, remnants of metal-rich and low entropy cool cores can persist after major mergers. While we obtain a mild anti-correlation between the Fe abundance and the pseudo-entropy in the (lower entropy, $K$ = 200-500 keV cm$^2$) inner regions, no clear correlation is found at (higher entropy, $K$ = 500-2300 keV cm$^2$) outer radii. The apparent spatial abundance uniformity that we find at large radii is difficult to explain through an efficient mixing of freshly injected metals, particularly in systems for which the time since the merger is short. Instead, our results provide important additional evidence in favour of the early enrichment scenario - in which the bulk of the metals are released outside galaxies at $z$ > 2-3 - and extend it from cool-core and (moderate) non-cool-core clusters to a few of the most disturbed merging clusters as well. These results constitute a first step towards a deeper understanding of the chemical history of merging clusters.

Tables

Table 1. Table 1: Observations and exposure times.

Name	Label	Sequence ID	Position (J2000)	Observation	MOS1	MOS2	pn
			(RA, Dec.)	starting date	Net (ks)	Net (ks)	Net (ks)
CIZA J2242.8+5301	CIZA2242	0654030101	(22:43:01.99, 53:07:30.0)	2010-12-13	71.5	71.9	68.4
1RXS J0603.3+4214	1RXSJ0603	0675060101	(06:03:13.39, 42:12:31.0)	2011-10-03	75.4	75.5	71.4
Abell 3376 East	A3376	0151900101	(06:02:08.59, -39:57:18.0)	2003-04-01	27.9	27.8	25.4
Abell 3376 West		0504140101	(06:00:51.14, -40:00:31.6)	2007-08-24	40.8	42.7	36.4
Abell 3667	A3667	0206850101	(20:13:04.79, -56:53:60.0)	2004-05-03	59.7	60.7	57.3
Abell 3667 offset		0653050201	(20:12:31.46, -56:21:39.3)	2010-09-21	28.9	28.6	24.7
Abell 665	A665	0109890501	(08:30:57.99, 65:50:20.0)	2001-09-23	80.3	81.4	–
Abell 2256	A2256	0401610101	(17:03:02.55, 78:45:00.0)	2006-08-04	39.7	42.6	6.6
		0141380101	(17:03:02.55, 78:45:00.0)	2003-04-27	9.8	10.1	9.5

Table 2. Table 2: Cluster sample description.

Label	$z$	$N_H$a	Radio features	Mass ratio	Dyn. stageb		Scale f.c
		(1020 cm-2)			(Gyr)	(keV)	(kpc/$\prime$)
CIZA2242	0.192	25.80–55.70	Doble relic, faint halo[1]	1–2:1[2,3]	$\sim$0.6[4]	7.9[4]	193.0
1RXSJ0603	0.225	21.50	Triple relic, halo[5]	3:1:0.1:0.1[6]	$\sim$2.0[7]	7.8[8]	218.0
A3376	0.046	4.84	Double relic[9]	3–6:1[10]	$\sim$0.5–0.6[10]	4.2[11]	54.0
A3667	0.055	4.44	Double relic, mini-halo[12,13]	5:1[14]	$\sim$1.0[14]	6.3[15]	64.5
A665	0.182	4.31	Halo[16]	1–2:1[17]	$\sim$1–2[17]	8.3[18]	184.5
A2256	0.058	4.24	Relic, halo[19,20]	3:1:0.3[21]	$\sim$–0.2[22]	6.4[23]	67.5

Table 3. Table 3: Averaged abundance profile.

Radius	$Z$	${\sigma }_{\text{stat}}$	${\sigma }_{\text{scatter}}$
(r/r500)	($Z_{\odot }$)
0.00–0.10	0.437	0.023	0.057
0.10–0.20	0.422	0.024	0.061
0.20–0.30	0.392	0.025	0.085
0.30–0.40	0.334	0.026	0.069
0.40–0.50	0.281	0.027	0.049
0.50–0.70	0.254	0.026	0.060
0.70–1.00	0.236	0.033	0.104

Table 4. Table 4: Best-fit parameters of abundance vs pseudo-entropy distributions.

K (keV cm2)	$\alpha$ (Z⊙ keV-1 cm-2)	Z0 ($Z_{\odot }$)
200–2300	0.32$\pm$0.06	0.51 $\pm$ 0.03
200–500	0.61$\pm$0.17	0.61 $\pm$ 0.05
500–2300	0.07$\pm$0.07	0.32 $\pm$ 0.05
10–1000a	1.49$\pm$0.51	0.86 $\pm$ 0.18

Table 5. Table 5: Averaged abundance distribution.

Radius	$Z$	${\sigma }_{\text{stat}}$	${\sigma }_{\text{scatter}}$
(r/r500)	($Z_{\odot }$)
0.00–0.10	0.434	0.024	0.071
0.10–0.20	0.391	0.025	0.060
0.20–0.30	0.356	0.024	0.070
0.30–0.40	0.316	0.025	0.059
0.40–0.50	0.279	0.027	0.048
0.50–0.70	0.255	0.026	0.060
0.70–1.00	0.236	0.033	0.104

Table 6. Table 6: Background components derived from the CIZA2242 observation.

	Norma	kT (keV)	$\Gamma$
MWH	57.55$\pm 0.28$	0.27 (fixed)	–
LHB	258.33$\pm 0.68$	0.08 (fixed)	–
CXB	6.75 (fixed)	–	1.41 (fixed)
HF	6.24$\pm 2.56$	0.7 (fixed)	–
SP MOS1	13.25$\pm 1.98$	–	1.40$\pm 0.19$
SP MOS2	18.10$\pm 1.77$	–	1.40$\pm 0.19$
SP PN	58.86$\pm 11.51$	–	1.20$\pm 0.15$

Table 7. Table 7: Background components derived from the 1RXSJ0603 observation.

	Norma	kT (keV)	$\Gamma$
MWH	2.58$\pm 0.07$	$0.21\pm 0.02$	–
LHB	13.40$\pm 0.06$	0.08 (fixed)	–
CXB	1.57 (fixed)	–	1.41 (fixed)
SP MOS1	5.22$\pm 0.46$	–	1.40$\pm 0.01$
SP MOS2	6.34$\pm 1.01$	–	1.40$\pm 0.01$
SP PN	42.71$\pm 0.90$	–	1.40$\pm 0.04$

Table 8. Table 8: Best-fit parameter values of the total background estimated for A3376 E observation.

	Norma	kT (keV)	$\Gamma$
MWH	2.31 (fixed)	0.27 (fixed)	–
LHB	86.00 (fixed)	0.08 (fixed)	–
CXB	3.90 (fixed)	–	1.41 (fixed)
SP MOS1	2.18$\pm 0.37$	–	1.40$\pm 0.04$
SP MOS2	1.30$\pm 0.37$	–	1.40$\pm 0.04$
SP PN	2.77$\pm 0.66$	–	1.09$\pm 0.12$

Table 9. Table 9: Best-fit parameter values of the total background estimated for A3376 W observation.

	Norma	kT (keV)	$\Gamma$
MWH	2.31 (fixed)	0.27 (fixed)	–
LHB	86.00 (fixed)	0.08 (fixed)	–
CXB	3.90 (fixed)	–	1.41 (fixed)
SP MOS1	2.87$\pm 0.40$	–	1.40$\pm 0.01$
SP MOS2	23.65$\pm 0.45$	–	1.40$\pm 0.01$
SP PN	7.57$\pm 0.57$	–	1.21$\pm 0.04$

Table 10. Table 10: Background components derived from the A3667 observation.

	Norma	kT (keV)	$\Gamma$
MWH	1.49$\pm 0.07$	$0.26\pm 0.01$	–
LHB	26.21$\pm 1.56$	0.08 (fixed)	–
CXB	0.74 (fixed)	–	1.41 (fixed)
SP MOS1	12.26$\pm 1.82$	–	1.40$\pm 0.09$
SP MOS2	7.55$\pm 3.12$	–	1.40$\pm 0.09$
SP PN	0.40$\pm 0.06$	–	0.35$\pm 0.07$

Table 11. Table 11: Background components derived from the A665 observation.

	Norma	kT (keV)	$\Gamma$
MWH	0.89$\pm 0.06$	$0.19\pm 0.03$	–
LHB	15.00$\pm 3.54$	0.08 (fixed)	–
CXB	0.41 (fixed)	–	1.41 (fixed)
SP MOS1	22.12$\pm 1.17$	–	1.40$\pm 0.04$
SP MOS2	9.88$\pm 0.48$	–	1.40$\pm 0.04$

Table 12. Table 12: Background components derived from the A2256 pointing 1 observation.

	Norma	kT (keV)	$\Gamma$
MWH	4.77$\pm 0.18$	$0.31\pm 0.04$	–
LHB	54.50$\pm 1.70$	0.08 (fixed)	–
CXB	8.17 (fixed)	–	1.41 (fixed)
SP MOS1	9.66$\pm 0.61$	–	0.91$\pm 0.02$
SP MOS2	9.99$\pm 0.72$	–	0.91$\pm 0.02$
SP PN	20.94$\pm 3.51$	–	1.08$\pm 0.08$

Table 13. Table 14: Best-fit parameters for CIZA2242 regions shown in Fig. 3 .

Region	Radius	kT	$Z$	Norm	NH	C-stat/d.o.f.
	(′)	(keV)	($Z_{\odot }$)	(1072 m-3 arcmin-2)	(1021 cm-2)
1	0.0 $\pm 1.5$	9.62 $\pm 0.54$	0.40 $\pm 0.09$	1.64 $\pm 0.04$	3.60 $\pm 0.10$	603/574
2	2.0 $\pm 0.5$	8.21$\pm 0.55$	0.19 $\pm 0.09$	1.98 $\pm 0.06$	3.79 $\pm 0.13$	553/539
3	3.0 $\pm 0.5$	9.14 $\pm 0.40$	0.26 $\pm 0.08$	1.63 $\pm 0.04$	3.85 $\pm 0.11$	570/560
4	4.3 $\pm 0.8$	8.07 $\pm 0.42$	0.29 $\pm 0.06$	1.46 $\pm 0.05$	3.85 $\pm 0.07$	624/606
5	6.0 $\pm 1.0$	8.68 $\pm 0.45$	0.40 $\pm 0.06$	0.95 $\pm 0.03$	4.27 $\pm 0.09$	653/619
5a (70–97$°$)	6.0 $\pm 1.0$	7.41$\pm$0.56	0.58$\pm$0.11	0.79 $\pm 0.03$	4.10 $\pm 0.13$	561/546
5b (97–117$°$)	6.0 $\pm 1.0$	8.69$\pm$0.67	0.29$\pm$0.09	1.44 $\pm 0.04$	4.32 $\pm 0.14$	525/559
5c (117–140$°$)	6.0 $\pm 1.0$	9.61$\pm$1.03	0.18$\pm$0.10	0.83 $\pm 0.03$	5.57 $\pm 0.07$	566/551
6	9.0 $\pm 2.0$	10.29 $\pm 0.79$	0.29 $\pm 0.10$	0.25 $\pm 0.01$	5.13 $\pm 0.07$	662/643
BCGN	6.0 $\pm 1.5$	8.16 $\pm 0.41$	0.26 $\pm 0.06$	1.37 $\pm 0.03$	4.33 $\pm 0.10$	645/591

Table 14. Table 15: Best-fit parameters for 1RXSJ0603 regions shown in Fig. 5 .

Region	Radius	kT	$Z$	Norm	C-stat/d.o.f.
	(′)	(keV)	($Z_{\odot }$)	(1072 m-3 arcmin-2)
1	-1.9 $\pm 0.4$	10.07 $\pm 0.46$	0.31 $\pm 0.09$	1.62 $\pm 0.09$	524/590
2	-1.1$\pm 0.4$	9.80$\pm 0.41$	0.38 $\pm 0.05$	4.16 $\pm 0.07$	606/625
3 (CORES)	0.0 $\pm 0.7$	8.43 $\pm 0.27$	0.51 $\pm 0.06$	7.02 $\pm 0.13$	613/572
4	2.0 $\pm 0.5$	8.92 $\pm 0.31$	0.37 $\pm 0.04$	8.85 $\pm 0.13$	640/647
5	3.0 $\pm 0.5$	8.57 $\pm 0.18$	0.24 $\pm 0.04$	9.20 $\pm 0.08$	717/652
6	4.0 $\pm 0.5$	7.86 $\pm 0.15$	0.19 $\pm 0.04$	8.10 $\pm 0.10$	783/652
7	5.0 $\pm 0.5$	7.42 $\pm 0.20$	0.14 $\pm 0.04$	5.08 $\pm 0.05$	729/631
8	6.0 $\pm 0.5$	6.92 $\pm 0.29$	0.20 $\pm 0.07$	2.43 $\pm 0.03$	619/602
9	10.0 $\pm 2.0$	2.93 $\pm 0.22$	0.16 $\pm 0.06$	0.36 $\pm 0.02$	974/661
COREN	3.0 $\pm 0.4$	8.85 $\pm 0.36$	0.37 $\pm 0.08$	5.05 $\pm 0.01$	606/562

Table 15. Table 16: Best-fit parameters for A3376 regions shown in Fig. 7 .

Region	Radius	kT	$Z$	Norm	C-stat/d.o.f.
	(′)	(keV)	($Z_{\odot }$)	(1070 m-3 arcmin-2)
1 (BCG2)	0.0 $\pm 2.0$	4.17 $\pm 0.13$	0.51 $\pm 0.05$	6.88 $\pm 0.12$	553/607
2	3.5 $\pm 1.5$	3.89$\pm 0.12$	0.44 $\pm 0.06$	7.14 $\pm 0.18$	642/579
3	6.5 $\pm 1.5$	4.00 $\pm 0.14$	0.37 $\pm 0.05$	5.61 $\pm 0.13$	647/603
4	9.5 $\pm 1.5$	3.79 $\pm 0.10$	0.31 $\pm 0.05$	4.54 $\pm 0.68$	1250/1185
5	12.5 $\pm 1.5$	3.88 $\pm 0.16$	0.30 $\pm 0.04$	4.29 $\pm 0.58$	1276/1195
6	16.0 $\pm 2.0$	4.21 $\pm 0.19$	0.29 $\pm 0.06$	1.63 $\pm 0.03$	617/643
7	20.0 $\pm 2.0$	5.36 $\pm 0.84$	0.15 $\pm 0.06$	0.68 $\pm 0.10$	657/599
BCG1	18.2 $\pm 2.0$	4.49 $\pm 0.35$	0.21 $\pm 0.12$	1.35 $\pm 0.01$	512/544

Table 16. Table 17: Best-fit parameters for A3667 regions shown in Fig. 9 .

Region	Radius	kT	$Z$	Norm	C-stat/d.o.f.
	(′)	(keV)	($Z_{\odot }$)	(1071 m-3 arcmin-2)
1	-7.6 $\pm 0.9$	5.46 $\pm 0.11$	0.43 $\pm 0.04$	1.24 $\pm 0.01$	802/664
2	-5.9 $\pm 0.9$	4.56$\pm 0.06$	0.52 $\pm 0.03$	2.57 $\pm 0.02$	843/682
3	-4.1 $\pm 0.9$	5.04 $\pm 0.07$	0.48 $\pm 0.03$	4.16 $\pm 0.02$	761/685
4	-2.4 $\pm 0.9$	5.54 $\pm 0.09$	0.45 $\pm 0.03$	6.12 $\pm 0.47$	756/651
5	0.0 $\pm 1.5$	6.79 $\pm 0.17$	0.38 $\pm 0.04$	4.29 $\pm 0.05$	760/652
6	2.4 $\pm 0.9$	7.05 $\pm 0.23$	0.39 $\pm 0.06$	3.54 $\pm 0.05$	675/618
7	4.1 $\pm 0.9$	7.29 $\pm 0.24$	0.30 $\pm 0.08$	2.35 $\pm 0.04$	729/631
8	6.8 $\pm 1.8$	6.79 $\pm 0.20$	0.27 $\pm 0.06$	1.48 $\pm 0.02$	733/673

Table 17. Table 18: Best-fit parameters for A665 regions shown in Fig. 11 .

Region	Radius	kT	$Z$	Norm	C-stat/d.o.f.
	(′)	(keV)	($Z_{\odot }$)	(1072 m-3 arcmin-2)
1	-3.5 $\pm 1.0$	6.73 $\pm 0.52$	0.29 $\pm 0.15$	0.26 $\pm 0.01$	393/400
2	-2.0 $\pm 0.5$	8.43 $\pm 0.85$	0.41 $\pm 0.17$	0.67 $\pm 0.02$	373/362
3	-1.0 $\pm 0.5$	7.30 $\pm 0.38$	0.57 $\pm 0.12$	2.51 $\pm 0.05$	413/372
4	0.0 $\pm 0.5$	7.80 $\pm 0.37$	0.40 $\pm 0.07$	11.79 $\pm 0.21$	431/421
5	1.0 $\pm 0.5$	7.39 $\pm 0.27$	0.28 $\pm 0.06$	7.17 $\pm 0.09$	452/405
6	2.0 $\pm 0.5$	7.36 $\pm 0.27$	0.33 $\pm 0.06$	3.35 $\pm 0.03$	476/405
7	3.0 $\pm 0.5$	8.72 $\pm 0.44$	0.23 $\pm 0.09$	1.66 $\pm 0.03$	412/395
8	4.0 $\pm 0.5$	8.32$\pm 0.59$	0.24 $\pm 0.12$	0.84 $\pm 0.02$	393/386
9	5.0 $\pm 0.5$	8.56 $\pm 0.79$	0.23 $\pm 0.15$	0.51 $\pm 0.01$	375/374

Table 18. Table 19: Best-fit parameters for A2256 regions shown in Fig. 13 .

Region	Radius	kT	$Z$	Norm	C-stat/d.o.f.
	(′)	(keV)	($Z_{\odot }$)	(1071 m-3 arcmin-2)
1	-4.0 $\pm 0.5$	9.03 $\pm 0.54$	0.28 $\pm 0.15$	1.86 $\pm 0.06$	1203/1095
2	-3.0 $\pm 0.5$	6.97 $\pm 0.27$	0.36 $\pm 0.09$	3.25 $\pm 0.08$	1161/1101
3	-2.0 $\pm 0.5$	6.56 $\pm 0.25$	0.37 $\pm 0.09$	3.83 $\pm 0.09$	1149/1081
4 (CORE1)	0.0 $\pm 1.5$	6.89 $\pm 0.15$	0.44 $\pm 0.06$	4.78 $\pm 0.07$	1218/1192
4	2.3 $\pm 0.8$	6.60 $\pm 0.20$	0.48 $\pm 0.07$	5.33 $\pm 0.09$	1077/1073
5	3.8 $\pm 0.8$	5.09$\pm 0.13$	0.43 $\pm 0.04$	5.41 $\pm 0.09$	1308/1145
6	5.3 $\pm 0.8$	4.24 $\pm 0.14$	0.49 $\pm 0.08$	0.70 $\pm 0.02$	1145/1024
CORE2	-3.3 $\pm 0.5$	4.70 $\pm 0.12$	0.56 $\pm 0.05$	5.79 $\pm 0.09$	1259/1096

Equations

$$r_{200}=2.77h_{70}^{-1}(⟨kT⟩/10\text{leavevmode}\mathrm{keV}{)}^{1/2}/\text{E(z)}\text{leavevmode}\mathrm{Mpc},$$

$$Z_{\text{ref}}(k)=\Bigg{(}\sum_{j=1}^{N}\sum_{i=1}^{M}w_{i,j,k}\frac{Z(i)_{j}}{\sigma^{2}_{Z(i)_{j}}}\Bigg{)}/\Bigg{(}\sum_{j=1}^{N}\sum_{i=1}^{M}w_{i,j,k}{\frac{1}{\sigma^{2}_{Z(i)_{j}}}}\Bigg{)},$$

$${\sigma }_{\text{stat}}(k)=\frac{1}{\sqrt{{\sum }_{j=1}^N{\sum }_{i=1}^M{w}_{i,j,k}\frac{1}{{\sigma }_{Z(i{)}_j}^2}}},$$

$${\sigma }_{\text{scatter}}(k)=\frac{\sqrt{{\sum }_{j=1}^N{\sum }_{i=1}^M{w}_{i,j,k}{\left(\frac{Z(i{)}_j-Z_{\text{ref}}(k)}{{\sigma }_{Z(i{)}_j}}\right)}^2}}{\sqrt{{\sum }_{j=1}^N{\sum }_{i=1}^M{w}_{i,j,k}\frac{1}{{\sigma }_{Z(i{)}_j}^2}}},$$

Full text

1 Introduction

In order to become the largest gravitationally bound structures seen in today’s Universe, galaxy clusters grow in a hierarchical way, via not only accretion but also merging of surrounding (sub) haloes. The latter case is particularly interesting as major cluster mergers are the most energetic events since the Big Bang. During the merging process, the hot intra-cluster medium (ICM) in galaxy clusters is violently compressed and heated, and large amounts of thermal and non-thermal energy are released, giving rise to shock fronts and turbulence (Markevitch & Vikhlinin 2007). Shocks can propagate into the ICM (re)accelerating electrons (Bell 1987; Blandford & Eichler 1987), which may produce elongated and polarized radio structures known as radio relics via synchrotron radio emission (for a review, see Ferrari et al. 2008; Brunetti & Jones 2014). On the other hand, particle (re)acceleration by means of the turbulence in the ICM could generate unpolarized cluster-wide sources known as radio haloes (for a review, see Brunetti et al. 2001; Feretti et al. 2012; van Weeren et al. 2019).

The ICM is also known to be rich in metals, which originate mainly from core-collapse (SNcc) and type Ia (SNIa) supernovae, which continuously release their nucleosynthesis products since the epoch of major star formation, about 10 billion years ago ($z$ $\sim$ 2–3) (Hopkins & Beacom 2006; Madau & Dickinson 2014). SNcc contribute primarily to the synthesis of light metals (O, Ne, Mg, Si and S). On the other hand, heavier metals (such as Ar, Ca, Mn, Fe and Ni) together with smaller relative contribution of Si and S are produced by SNIa. Finally, C and N are mostly released by low-mass stars on the asymptotic giant branch (AGB). Once the metals are ejected to the ICM by SN or AGB stars, they are later transported, mixed and redistributed via other processes, such as galactic winds (Kapferer et al. 2006, 2007; Baumgartner & Breitschwerdt 2009); ram-pressure stripping (Schindler et al. 2005; Kapferer et al. 2007), active galactic nucleus (AGN) outflows (Simionescu et al. 2008, 2009; Kapferer et al. 2009) or gas sloshing (Simionescu et al. 2010; Ghizzardi et al. 2014), among other processes. Detailed reviews on the observation of metals in the ICM include, e.g. Werner et al. (2008); Böhringer & Werner (2010); de Plaa (2013) and Mernier et al. (2018b).

In the last years, the iron (Fe) abundance distribution in galaxy clusters has been extensively studied using the different X-ray observatories (De Grandi & Molendi 2001; De Grandi et al. 2004; Baldi et al. 2007; Leccardi & Molendi 2008; Maughan et al. 2008; Matsushita 2011; Mernier et al. 2016, 2017; Mantz et al. 2017; Simionescu et al. 2017, 2019). The Fe radial distribution is different for cool-core (CC, i.e. dynamically relaxed) clusters and non-cool-core (NCC, i.e. dynamically disturbed) clusters. In the case of CC clusters, the Fe distribution shows a peak in the core, which decreases up to $\sim$0.3$r_{500}$ and flattens for larger radii. NCC clusters, however, have a flat distribution in the core and follow the same universal distribution as CC for outer radii. These abundance distributions do not present evidence of evolution up to $z$ $\sim$ 1 (McDonald et al. 2016; Mantz et al. 2017; Ettori et al. 2017). Moreover, recent X-ray observations with Suzaku (Fujita et al. 2008; Werner et al. 2013; Urban et al. 2017; Ezer et al. 2017; Simionescu et al. 2017) reveal an uniform radial distribution around $\sim$0.2–0.3 $Z_{\odot}$ in the outskirts of the galaxy clusters, which suggests an early enrichment of the ICM (i.e. $z$ $\sim$ 2–3). Numerical simulations later showed that active galactic nuclei (AGN) feedback at early times is needed to explain these observations (e.g. Biffi et al. 2017, 2018a, 2018b).

Whereas metals in the ICM have been extensively and systematically studied in CC and, to some extent, NCC clusters, the ICM enrichment in the extreme case of merging clusters has been poorly explored so far (except in a few specific cases; e.g. A3376, Bagchi et al. 2006; A3667, Lovisari et al. 2009; Laganá et al. 2019). Knowing in detail how metals are distributed across merging clusters, however, is a key ingredient in our understanding of the chemical evolution of large scale structures. On one hand metals act as passive tracers and may reveal valuable information on the dynamical history of the recently merged ICM. On the other hand, merging shocks could in principle trigger star formation in neighbouring galaxies (Stroe et al. 2015a; Sobral et al. 2015), hence potentially contributing to the ICM enrichment to some extent (Elkholy et al. 2015). Metallicities in the shocked gas of such disturbed clusters, and their connection with their thermodynamical properties and star formation of their galaxies remain widely unexplored.

In this work, we present for the first time a study specifically devoted to the chemical state of merging clusters. For this purpose, we derive the spatial distribution of temperature, abundance (Fe), and pseudo-entropy from specifically selected regions of six merging galaxy clusters: CIZA J2242.8+5301, 1RXS J0603.3+4214, Abell 3376, Abell 3667, Abell 665 and Abell 2256. The principal characteristics of these clusters are summarized in Table 2. The first four clusters of the list are essentially major mergers hosting a double radio relic, while Abell 665 hosts a giant radio halo and Abell 2256 a prominent single radio relic together with a radio halo. The selection criteria of these clusters have been: (i) XMM-Newton observation availability, (ii) net exposure time ¿ 25 ks (iii) bright and disturbed central ICM and (iv) probing a variety of diffuse radio features: we include merging clusters with double radio relics, a single radio relic and strong radio halo. This is not intended as a complete sample but rather as an exploratory first step towards future studies on metals in the ICM of merging clusters from larger observational samples. For this study, we report our results with respect to the protosolar abundances ($Z_{\odot}$) reported by Lodders et al. (2009). We assume the cosmological parameters $H_{0}$ = 70 km/s/Mpc, $\Omega_{\rm{M}}$ = 0.27 and $\Omega_{\Lambda}$= 0.73. All errors are given at 1$\sigma$ (68$\%$) confidence level unless otherwise stated and the spectral analysis uses the modified Cash statistics (Cash 1979; Kaastra 2017a).

2 Observations and data reduction

Table 1 summarizes the XMM-Newton observations, the different pointings and the net exposure time used for the cluster survey analysis. All the data are reduced with the XMM-Newton Science Analysis System (SAS) v17.0.0, with the calibration files dated by June 2018. We use the observations from the MOS and pn detectors of the EPIC instrument and reduce them first with the standard pipeline commands emproc and epproc. Next, we filter the soft-proton (SP) flares by building Good Time Interval (GTI) files following the method detailed in Appendix A.1 of Mernier et al. (2015). This method consists of extracting light curves at the high energy range of spectrum (10–12 keV for MOS and 12–14 keV for pn) in bins of 100s. Then, the mean count rate, $\mu$, and the standard deviation, $\sigma$, are calculated in order to apply a threshold of $\mu\leavevmode\nobreak\ \pm\leavevmode\nobreak\ 2\sigma$ to the generated distribution. Therefore, all the time bins with a number of counts outside the interval $\mu\leavevmode\nobreak\ \pm\leavevmode\nobreak\ 2\sigma$ are rejected. Afterwards, we apply the same process at the low energy range 0.3–2 keV with bins of 100s (Lumb et al. 2002). Only events satisfying the 2$\sigma$ threshold in the two energy ranges mentioned above are further considered. The point source identification in the complete field of view (FOV) is done with the SAS task edetect_chain applied in four different bands (0.3–2.0 keV, 2.0–4.5 keV, 4.5–7.5 keV, and 7.5–12 keV). After this, we combine the detections in these four bands and exclude the point sources with a circular region of 10$\arcsec$ radius, except some larger sources where an appropriate radius has been applied. The 10$\arcsec$ radius size, as explained in Appendix A.2 of Mernier et al. (2015), is optimized for avoiding the polluting flux of the point sources and not reducing considerably the emission of the ICM. We keep the single, double, triple and quadruple events in MOS (pattern$\leq$12) and only the single pixel events in pn (pattern==0). We also correct the out-of-time events from the pn detector. 111https://www.cosmos.esa.int/web/xmm-newton/sas-thread-epic-oot

3 Spectral analysis

3.1 Spectral analysis approach

In our spectral analysis of the clusters, we assume that the observed spectra include the following components (see EPIC MOS2 spectrum in Fig. 1): an optically thin thermal plasma emission from the ICM in collisional ionization equilibrium (CIE), the local hot bubble (LHB), the Milky Way halo (MWH), the cosmic X-ray background (CXB) and the non X-ray background (NXB), consisting of hard particle (HP) background and residual soft-proton (SP) component (see Sect. 3.2). An additional hot foreground (HF) component is included for CIZA J2242.8+5301 as suggested by Ogrean et al. (2013a) and Akamatsu et al. (2015). Further details about these components (and how we include them in the fits) are given in Sect. 3.2. The emission models are corrected for the cosmological redshift and absorbed by the galactic interstellar medium. The hydrogen column density value adopted for all the clusters is the weighted neutral hydrogen column density, $N_{\rm{HI}}$, estimated using the method of Willingale et al. (2013)222http://www.swift.ac.uk/analysis/nhtot/. Previous work (e.g. Mernier et al. 2016; de Plaa 2017) noted that the total $N_{\rm{H}}$ value (neutral + molecular hydrogen) may sometimes provide inaccurate fits whereas the weighted neutral hydrogen $N_{\rm{HI}}$ values are generally closer to the best-fit $N_{\rm{H}}$ values. Specifically, we have checked that leaving free the $N_{\rm{H}}$ in our fits provides values that do not deviate more than $\sim$10$\%$ from the above estimates. The only exception is CIZA J2242.8+530, for which we need to leave the value of $N_{\rm{H}}$ free within the range 0.8 x $N_{\rm{HI}}\leavevmode\nobreak\ \leq\leavevmode\nobreak\ N_{\rm{H}}\leavevmode\nobreak\ \leq$ 1.2 x $N_{\rm{H,tot}}$, where $N_{\rm{H,tot}}$ is the weighted total (neutral and molecular) hydrogen column density (Willingale et al. 2013).

In our spectral analysis, we use SPEX333https://www.sron.nl/astrophysics-spex (Kaastra et al. 1996, 2017b) version 3.04.00 with SPEXACT (SPEX Atomic Code and Tables) version 3.04.00. We carry out the spectral fitting in different regions as detailed in the sections below. The redistribution matrix file (RMF) and the ancillary response file (ARF) are processed using the SAS tasks rmfgen and arfgen, respectively. The spectra of the MOS (energy range from 0.5 to 10 keV) and pn (0.6–10 keV) detectors are fitted simultaneously and binned using the method of optimal binning described in Kaastra & Bleeker (2016). We assume a single temperature structure in the Collisional Ionization Equilibrium (cie) component for all the clusters (this assumption is justified in Sect. 3.3). Because the high temperature of merging clusters does not allow to measure other metals with reasonable accuracy, their abundances are coupled to Fe. The free parameters considered in this study are the temperature kT, the metal abundance Z and the normalization Norm. In the case of CIZA J2242.8+5301, the column density $N_{\rm{H}}$ is also free to vary.

3.2 Estimation of background spectra

A proper estimation of the background components is important especially in the regions where the emission of the cluster is low. This is particularly relevant outside of the core, where the spectrum can be background dominated. In this work, we model the background components directly from our spectra (Sect. 3.1) using the following list of components:

• •

The LHB, modelled by a non-absorbed $cie$ component with proto-solar abundance. The temperature is fixed to 0.08 keV.

• •

The MWH, modelled by an absorbed $cie$ component with proto-solar abundance. The temperature is free to vary except for CIZA J2242.8+5301 and Abell 3376, where it is fixed to 0.27 keV, following the values of Akamatsu et al. (2015) and Urdampilleta et al. (2018), respectively.

• •

The CXB, modelled by an absorbed power law with a fixed $\Gamma$ = 1.41 (De Luca & Molendi 2004). The normalization is a fixed value for each cluster assuming the CXB flux of 8.07 x 10*-15* W m*-2* deg*-2*, taking into account a detection limit of Sc = 3.83 x 10*-15* W m*-2* and calculated with the CXBTools (de Plaa 2017). See details in Appendix B.2 of Mernier et al. (2015).444Note that the flux of the faintest resolved and excluded point source, and hence the remaining unresolved CXB flux, depend on the exposure time. In our case, most of the observations have similar a exposure except two: 0151900101 and 0653050201. Even for these data sets, the SP power law norm is free to vary and will compensate for any residuals in the CXB subtraction.

• •

The HF, modelled by an absorbed $cie$ component with proto-solar abundance. The temperature is fixed to 0.7 keV.

• •

The HP background, modelled by a broken power law (unfolded by the effective area) together with several instrumental Gaussian profiles (eight and nine for MOS and pn, respectively). The photon indices of the power laws and the centroid energies of the profiles have been taken from Table B.1 and Table B.2 of Mernier et al. (2015), see their Appendix B.1 for more details. The normalizations are left free to vary.

• •

The SP background, modelled by a power law (unfolded by the effective area). The power law index of the SP component is left free to vary between 0.1 and 1.4.

The estimation of the sky background components (LHB, MHW and HF) is obtained from the offset regions (see green sectors in Fig. 3, 5, 9, 11 and 13) where a negligible (yet still modelled) thermal emission of the cluster is expected and the background is dominant. For Abell 3376 the values of Urdampilleta et al. (2018) are adopted. The parameters of the SP components are obtained fitting the total FOV (r\leavevmode\nobreak\ =\leavevmode\nobreak\15$\arcmin$) EPIC spectra, where the cie (i.e. ICM) parameters are also left free. Finally, the normalization of every background component (except the HP background) has been fixed, rescaled on the sky area of each region and corrected for vignetting if necessary. The best-fit parameters for each cluster are listed in Appendix A.

3.3 Systematic uncertainties

We account for the effect of systematic uncertainties related to the $\pm$10$\%$ variation in the normalization of the sky background and non X-ray background components (Mernier et al. 2017). The contribution of these systematic uncertainties are further included in our analysis and results.

Mernier et al. (2017) and the posterior review of Mernier et al. (2018b) list other potential systematic uncertainties, which could affect the abundance measurements of this work. These include MOS-pn discrepancies due to residual EPIC cross-calibration, projection effects, atomic code and thermal structure uncertainties. The former has limited impact on the relative Fe abundance distribution as explained by Mernier et al. (2017). Although recent work of Liu et al. (2018) suggests that projection effects may have non-negligible observational biases on the observed abundance evolution in the core of CC clusters, the lack of strong emission gradients in merging clusters makes us confident that such a bias is limited in our case. Moreover, deprojecting merging cluster emission is very challenging, as by definition these systems are far from being spherically symmetric.

Uncertainties regarding the atomic codes are limited by the fact that we are using the update of the code SPEXACT and SPEX, which incorporate more precise modelling of the atomic processes and extensive compilation of transitions. However, uncertainties due to a not complete atomic database are still present (Hitomi Collaboration et al. 2018) and their effects can be non-negligible (Mernier et al. 2018a).

For the thermal structure, we assume a single-temperature (1T) model. This is a simplified approximation; other more complex models reproducing Gaussian (gdem, de Plaa et al. 2006) or power-law (wdem, Kaastra et al. 2004) temperature distributions could represent better the inhomogeneities of the ICM. To quantify these possible effects, we have fitted the spectra assuming a simplified version of the gdem model using a 3T model with one main temperature free and the other two coupled with a factor of 0.64 and 1.56, respectively. We have realized the fits in the SPEXACT version 2.07.00 in order to significantly reduce the computing time. This version includes simpler atomic code and tables, which allows to fit the spectra faster, but with a similar accuracy as the most updated version (for more details, see Mernier et al. 2018a). We verify that the best-fit temperature of the main component in our 3T model remains within $\pm$ 5$\%$ of the best-fit temperature obtained when assuming a 1T model. The relative emission measure of the other two components are less than 1$\%$ of the best-fit emission measure of the 1T model, and the variation in the abundance is less than 10$\%$. This 1T-like behaviour is not surprising as, unlike CC clusters, merging clusters do not exhibit a strong temperature gradient, and at these high temperatures the Fe abundance is almost entirely determined from the Fe-K line (which is unaffected by the Fe-bias or the inverse Fe-bias, e.g. Buote 2000, Simionescu et al. 2009). For this reason and for simplicity we decide to maintain the 1T model and we will add these uncertainties to the other background systematic uncertainties as the thermal structure contribution.

4 Spectral analysis: individual clusters

4.1 CIZA J2242.8+5301

CIZA J2242.8+5301 (hereafter CIZA2242, also known as the ”Sausage” cluster) is a massive ($\sim$1015M⊙) and two-body (Jee et al. 2015) merging galaxy cluster at a redshift of $z$ = 0.192. It was discovered in the second CIZA sample from the ROSAT All-Sky Survey and identified as a major-merger by Kocevski et al. (2007). van Weeren et al. (2010) studied CIZA2242 for the first time using WSRT, GMRT and VLA. They observed a double radio relic feature located at $\sim$1.5 Mpc from the cluster center. XMM-Newton and Chandra observations (Ogrean et al. 2013a, 2014) reveal a strongly disturbed and complex merger. In fact, the ICM presents an elongated X-ray morphology with a bullet shape bright edge in the south (S) and more irregular edge in the north (N), suggesting that the infalling direction is N to S. Akamatsu et al. (2015) estimate the time after the core passage to be $\sim$0.6 Gyr, in good agreement with the Sunyaev–-Zel’dovich (SZ) analysis of Rumsey et al. (2017).

Finally, Stroe et al. (2014a); Stroe & Sobral (2015b); Sobral et al. (2015); Stroe et al. (2017) confirm an overdensity of Hα emitters, mainly cluster star-forming galaxies and AGNs, concentrated close to the subcores and northern post-shock region. The star-forming galaxies are found to be highly metal rich and to have strong outflows. The aboved mentionned authors propose that such an enhancement of the star formation ratio could be a shock-induced effect.

4.1.1 Results

We fit the cluster emission and the background spectra (see Table 6) in the sectors shown in Fig. 5, having excluded the point sources. We leave $N_{\rm{H}}$ free to vary between 0.8 x $N_{\rm{HI}}$ $\leq$  $N_{\rm{H}}$ $\leq$ 1.2 x $N_{\rm{HTotal}}$, where $N_{\rm{HI}}$ = 3.22 $\times$ 1021 cm*-2* and $N_{\rm{HTotal}}$ = 4.64 $\times$ 1021 cm*-2* (Willingale et al. 2013), as explained in Section 3.1. The two circular regions are centred close to the two brightest cluster galaxies (BCGs) or subcores of the bimodal merger system described by Dawson et al. (2015) (BCGN: $22^{\rm{h}}42^{\rm{m}}50\aas@@fstack{s}00$, $53\degr 05\arcmin 06\aas@@fstack{\prime\prime}00$ ; BCGS: $22^{\rm{h}}42^{\rm{m}}39\aas@@fstack{s}00$, $52\degr 58\arcmin 35\aas@@fstack{\prime\prime}00$). We use as centroid the southern ”core”, because it contains the X-ray peak emission and is close to the southern BCG.

The best-fit parameters are shown in Table 14, while the temperature and abundance distributions along the merging axis are plotted in Fig. 3. By convention, and for all the further relevant figures in this paper, the regions lying behind the projected main BCG trajectory are plotted with positive values on the x-axis. The temperature distribution presents a variation in the range between $\sim$8–10 keV within the inner regions of the X-ray emission. This temperature is higher in the last annular region, coincident with the post-shock region as described in Ogrean et al. (2014) and Akamatsu et al. (2015). These results are in good agreement within the uncertainties with the temperature map obtained by Ogrean et al. (2013a) and the radial profile of Akamatsu et al. (2015).

Between the two BCGs (sectors 2 to 4) the Fe abundance is measured at its lowest value, between $\sim$0.2–0.3 $Z_{\odot}$. Moreover, the abundance profile shows an apparent enhancement with a maximum value of $\sim$0.4 $Z_{\odot}$ in sectors 1 and 5, where the two BGCs are located. The northern BCG region evinces a lower abundance than the entire sector 5. If we divide sector 5 in three regions, we obtain a gradient in the abundance and slight increase in the temperature towards the east, see Table 14. BCGN coincides mostly with region 5b and the enhancement in abundance in sector 5 is mainly driven by region 5a.

4.2 1RXS J0603.3+4214

1RXS J0603.3+4214 (hereafter 1RXSJ0603), also known as the ”Toothbrush” cluster, is a bright, rich and massive cluster discovered by van Weeren et al. (2012a) and located at $z$ = 0.225. It hosts three radio relics and a giant ($\sim$2 Mpc) radio halo elongated along the main merger axis S–N, which follows the X-ray emission morphology excluding the southernmost part (Ogrean et al. 2013b; Rajpurohit et al. 2018). Stroe et al. (2014b) found no enhancement of Hα emitters near the northern radio relic, with a star formation rate (SFR) consistent with that of blank field galaxies at $z$ = 0.2. They concluded that the lack of Hα emitters could be possible because 1RXSJ0603 is a well evolved system, approximately $\sim$2 Gyr after the core passage.

4.2.1 Results

As shown in Fig. 5, we use various sectors to fit the cluster emission and the background spectra (see Table 7), excluding the point sources and using $N_{\rm{HI}}$ = 2.15 $\times$ 1021 cm*-2* (Willingale et al. 2013). We follow the X-ray emission S–N shape, coincident with the radio halo elongation with boxes 4–8, as well as the pre-shock region 9. We dedicate two circular sectors to the subcores of the main two subclusters, locating the brightest one in X-rays at the south (CORES: $6^{\rm{h}}03^{\rm{m}}13\aas@@fstack{s}39$, $42\degr 12\arcmin 31\aas@@fstack{\prime\prime}00$) and the second at the north (COREN: $6^{\rm{h}}03^{\rm{m}}19\aas@@fstack{s}17$, $42\degr 14\arcmin 24\aas@@fstack{\prime\prime}00$). We use CORES as centroid for our temperature and abundance distributions along the merging axis (see Fig. 5).

The best-fit parameters are presented in Table 15 and are in good agreement with the previous X-ray studies by Ogrean et al. (2013b); Itahana et al. (2015); and van Weeren et al. (2016). The subcore temperatures ($kT$ = 8.43 $\pm$ 0.27 keV for CORES and $kT$ = 8.85 $\pm$ 0.36 keV for COREN) are consistent with van Weeren et al. (2016). Interestingly, we find a peak in Fe abundance located in the inner core, just behind the cold front (CF) (see Fig. 5; its location is approximated based on van Weeren et al. 2016). In the south, we see a slight temperature gradient towards the southern shock front, anticorrelated with a decrease in Fe abundance. The COREN shows the same abundance as the adjacent inner region, but a higher value than in the outer sectors (6–9). It is possible that the gas belonging to this subcluster has been stripped during the merger and is mixing in that region. In the north, a gradual smooth decrease of temperature can be observed up to the post-shock region, followed by a significant drop after the radio relic. Meanwhile, the Fe abundance decreases in sector 5 and flattens later around $\sim$0.2 $Z_{\odot}$ for the outer sectors.

4.3 A3376

Abell 3376 (hereafter A3376) is a bright and nearby ($z$ = 0.046) merging galaxy cluster located in the southern hemisphere. The N-body hydrodynamical simulations of Machado & Lima Neto (2013) suggest a two body merger scenario with a mass ratio of 3–6:1. The more massive and diffuse western subcluster core (BCG1: $6^{\rm{h}}00^{\rm{m}}41\aas@@fstack{s}10$, $-40\degr 02\arcmin 40\aas@@fstack{\prime\prime}00$) has been disrupted by a dense and compact eastern subcluster (BCG2: $6^{\rm{h}}02^{\rm{m}}09\aas@@fstack{s}70$, $-39\degr 57\arcmin 05\aas@@fstack{\prime\prime}00$). Recent studies suggest that the core-passing took place $\sim$0.5–0.6 Gyr ago in the W–E direction (Durret et al. 2013; Monteiro-Oliveira et al. 2017; Urdampilleta et al. 2018). Due to the merger, the inner gas of A3376 follows a cometary tail X-ray morphology as shown in Fig. 7. A3376 hosts two giant ($\sim$Mpc) arc-shaped radio relics (Bagchi et al. 2006) in the periphery of the cluster. X-ray analyses with the Suzaku satellite (Akamatsu et al. 2012; Urdampilleta et al. 2018) confirm the presence of a X-ray shock front located at the western radio relic, another shock probably associated with the eastern ”notch” (Paul et al. 2011; Kale et al. 2012) and a cold front at 3$\arcmin$ from the center.

4.3.1 Results

In the spectral analysis of A3376 we use two different observations as listed in Table 1. Specifically, ObsI:0151900101 is used for regions 1, 2 and 3, while ObsID:0504140101 is used for regions 6 and 7. In addition, both these observations are used for the overlapping regions 4 and 5. The background modelling parameters for each of the observations are shown in Table 8 and Table 9. We fit the spectra of the annular sectors centred on BCG2, located close to the X-ray peak emission at the eastern subcluster. We assume the column density as $N_{\rm{HI}}$ = 4.84 $\times$ 1020 cm*-2* (Willingale et al. 2013). The best-fit parameters (see Table 16) are in good agreement with previous studies with the Suzaku satellite (Akamatsu et al. 2012; Urdampilleta et al. 2018).

The temperature distribution of A3376 (Fig. 7) shows a central region with an average temperature $\sim$4 keV out to the outer radii where an increase of the temperature is seen. It is associated to the post merger region behind the western X-ray shock as described by Urdampilleta et al. (2018). The Fe abundance profile shows a peak of $\sim$0.5 $Z_{\odot}$ in the core, coincident with the BCG2. At larger distances, the abundance decreases smoothly to reach $\sim$0.3 $Z_{\odot}$ in regions 4, 5, 6, including BCG1 (dashed annular sector in Fig. 7) and in the surrounding diffuse gas. The lower abundance value appears in the outermost region of the cluster before the X-ray shock.

4.4 A3667

Abell 3667 (hereafter A3667) is a widely studied bright, low redshift ($z$ = 0.0553) and bimodal merging galaxy cluster (Markevitch et al. 1999). As a consequence of a violent merger along the northwest (NW) direction, A3667 hosts two curved radio relics to the NW and SE (Röttgering et al. 1997; Johnston-Hollitt 2003; Johnston-Hollitt & Pratley 2017; Carretti et al. 2013; Hindson et al. 2014; Riseley et al. 2015). The inner central region shows an elongated morphology of high-abundance gas in SE-NW direction (Mazzotta et al. 2002; Briel et al. 2004; Lovisari et al. 2009) according to the subcluster infall direction. A3667 also includes a prominent mushroom-like cold front close to the center at SE (Vikhlinin et al. 2001; Briel et al. 2004; Owers et al. 2009; Datta et al. 2014; Ichinohe et al. 2017).

4.4.1 Results

We fit the cluster emission and the background spectra (see Table 10) in the elliptical regions as shown in Fig. 9, excluding the point sources and using $N_{\rm{HI}}$ = 4.44 $\times$ 1020 cm*-2* (Willingale et al. 2013). We use as centroid the main BCG ($20^{\rm{h}}12^{\rm{m}}27\aas@@fstack{s}43$, $-56\degr 49\arcmin 35\aas@@fstack{\prime\prime}85$), which is close to the X-ray emission peak ($20^{\rm{h}}12^{\rm{m}}27\aas@@fstack{s}$, $-56\degr 50\arcmin 11\aas@@fstack{\prime\prime}00$). We distribute the elliptical sectors along the merging axis assuming the BCG as the center and following the X-ray morphology described by Briel et al. (2004).

The best-fit parameters are summarized in Table 18 and represented in the distributions along the merging axis in Fig 9. The results are in good agreement with previous studies (Briel et al. 2004; Lovisari et al. 2009; Akamatsu et al. 2013a; Datta et al. 2014; Ichinohe et al. 2017). The distributions show a clear temperature decrease towards the cold front at SE. It has an opposite trend with the Fe abundance, which reaches a peak just behind the cold front. The value presented here is lower than the maximum found by Briel et al. (2004) and Lovisari et al. (2009) of $\sim$0.7 $Z_{\odot}$, probably because our value includes larger azimuthal angles with lower abundance. However, if we select a smaller region in sector 2 (yellow dashed lines in Fig. 9), the temperature drops to 3.77 $\pm$ 0.09 keV and the abundance increases up to 0.64 $\pm$ 0.06 $Z_{\odot}$. As pointed out by Lovisari et al. (2009), this behaviour agrees with the simulations of Heinz et al. (2003), which suggest a possible displacement of the central rich and cold gas of the subcluster center well to the front during the infall, due to the internal dynamics of the cold front. The sectors in the NW direction (6 to 9) show a temperature distribution close to $\sim$7 keV, i.e. the average temperature of the cluster (Markevitch et al. 1999), and consistent with the hot tail described by Briel et al. (2004). The abundance in the central regions (2 to 4) is almost uniform, around $\sim$0.4 $Z_{\odot}$, but higher than the main cluster value ($\sim$0.2–0.3 $Z_{\odot}$). Finally, the abundance is again flat in the outermost regions (7 to 8) close to $\sim$0.3 $Z_{\odot}$.

4.5 A665

The two-body merging cluster Abell 665 (hereafter A665) was discovered as a non-relaxed cluster in an optical observation of Geller & Beers (1982). The authors found the first evidence of an elongated morphology and non-relaxed galaxy cluster. A665 is a hot cluster with an average temperature of $\sim$8.3 keV (Hughes & Tanaka 1992) and located at $z$ = 0.182 (Gomez et al. 2000). The N-body hydrodynamical simulations and dynamic studies of Gomez et al. (2000) suggest a subcluster merger in the NW–SE direction with two similar masses. The bright core of the infalling subcluster is moving south, crossing a more diffuse halo, before being stripped by ram pressure. A665 hosts a giant radio halo (Jones & Saunders 1996; Feretti et al. 2004; Vacca et al. 2010), which follows the X-ray emission elongation. The X-ray observations of A665 (Markevitch & Vikhlinin 2001; Govoni et al. 2004; Dasadia et al. 2016a) have revealed only one X-ray peak ($8^{\rm{h}}30^{\rm{m}}59\aas@@fstack{s}8$, $65\degr 50\arcmin 31\aas@@fstack{\prime\prime}3$), close to the BCG, which indicates that only the infalling cluster core has survived the merger. The remnant core presents a nearby cold front in the SE, followed by a hot gas ($\sim$8–15 keV) thought to be a post-merger region of a shock (Markevitch & Vikhlinin 2001; Govoni et al. 2004).

4.5.1 Results

We use spectra from MOS, the only detector available for these observations, to fit the cluster emission together with the modelled background (Table 11). We analyse the circular sectors centred on the BCG ($8^{\rm{h}}30^{\rm{m}}57\aas@@fstack{s}.6$, $65\degr 50\arcmin 29\aas@@fstack{\prime\prime}6$) as shown in Fig. 11, using $N_{\rm{HI}}$ = 4.31 $\times$ 1020 cm*-2* (Willingale et al. 2013). The best-fit parameters are presented in the Table 19 and plotted in Fig. 11.

The temperature and Fe abundance distributions along the merging axis show the presence of the two cold fronts as described by previous X-ray studies (Markevitch & Vikhlinin 2001; Govoni et al. 2004; Dasadia et al. 2016a). We see a slight temperature decrease across them towards the center. The abundance peak is found just within the south cold front and reflects the same case described by Heinz et al. (2003) for A3667. The diffuse gas in the north presents a lower abundance and a more uniform distribution ($\sim$0.2–0.3 $Z_{\odot}$) than the gas in the core and just behind the cold front in the south. This could indicate that the gas belongs to the more diffuse cluster, which had been already disrupted by the core-crossing. The two cold fronts seem to delimit a central region with a lower temperature ($\sim$7 keV) than the surrounding gas. We found as well the cooler gas region ($kT$ = 7.4 $\pm$ 0.2, $Z$ = 0.33$\pm$ 0.06 $Z_{\odot}$) in the north observed by Markevitch & Vikhlinin (2001) and shown in Fig. 11 as a magenta box. This cooler region seems to belong to the infalling cluster coming from the NW direction as described by Markevitch & Vikhlinin (2001) and Govoni et al. (2004).

4.6 A2256

Abell 2256 (hereafter A2256) is a nearby ($z$ = 0.058) rich galaxy cluster with bright X-ray emission. Optical analyses (Fabricant et al. 1989; Berrington et al. 2002; Miller et al. 2003) revealed a complex dynamical state, which consists of at least a triple merging with a mass ratio of 3:1:0.3. Berrington et al. (2002) propose that the system includes a main cluster with an infalling subcluster NW at $\sim$0.2 Gyr from the core passage (Roettiger et al. 1995) and another infalling ”group” northeast (NE), possibly associated to a previous merger. Briel et al. (1991) detected two separate X-ray peaks corresponding to the main cluster and the infalling subcluster, lately confirmed by Sun et al. (2002) and Bourdin & Mazzotta (2008). They also found bimodality in the temperature structure along the cluster elongation axis, being the low temperature component ($\sim$4.5 keV) associated with the subcluster and the higher temperature component ($\sim$7–8 keV) with the main cluster. In addition, Tamura et al. (2011) found bulk motion of gas of the cooler temperature component.

4.6.1 Results

We use two observations (see Table 1) for the analysis of A2256. We divide the ICM emission region in circular sectors along the cluster major elongation axis centered in the X-ray peak of the main cluster (CORE1: $17^{\rm{h}}04^{\rm{m}}13\aas@@fstack{s}82$, $78\degr 38\arcmin 15\aas@@fstack{\prime\prime}0$, Bourdin & Mazzotta 2008). We fit the cluster emission and the background spectra (see Table 13) using $N_{\rm{HI}}$ = 4.24 $\times$ 1020 cm*-2* (Willingale et al. 2013). The best-fit parameters are presented in the Table 19 and plotted in Fig. 13.

The temperature profile shows the bimodal temperature structure already described by Sun et al. (2002) and Bourdin & Mazzotta (2008), as well as the hot ($\sim$9 keV) temperature component south of the cluster center. It shows a temperature discontinuity from the adjacent region at $\sim$7 keV, suggested to be caused by the presence of a cold front (Bourdin & Mazzotta 2008). However, and similarly to Bourdin & Mazzotta (2008), we do not find evidence for the presence of the ”shoulder” mentioned by Sun et al. (2002) east of the CORE1. The Fe abundance distribution is in agreement with the values obtained by Sun et al. (2002). We obtain a quasi uniform abundance distribution (0.3–0.5 $Z_{\odot}$) with a slight decrease at the south. The second X-ray peak, (CORE2: $17^{\rm{h}}03^{\rm{m}}07\aas@@fstack{s}43$, $78\degr 39\arcmin 07\aas@@fstack{\prime\prime}8$, Sun et al. 2002), associated to the infalling subcluster, presents the highest abundance value of $\sim$0.6 $Z_{\odot}$.

5 Discussion

5.1 Averaged abundance distribution

The individual Fe abundance distributions of the six merging galaxy clusters along their merging axis are compiled in Fig. 14 (left). We remind the reader that such distributions are not azimuthally averaged. Instead, in this work we have chosen to trace the central ICM elongation along the merging axis, essentially because the surface brightness allows better constraints on the Fe abundance and because no spectacular metal redistribution is expected to take place perpendicular to the merging axis. In the case of CIZA2242, 1RXSJ0603 and A3376 the interest region clearly extends the core (where the BCG resides) to the outwards direction. For A3667, A665 and A2256, sectors along two opposite directions from the core are used. In Fig. 14, all the clustercentric radii have been rescaled to fractions of $r_{500}$555$r_{\Delta}$ is formally defined as the radius within which the density is $\Delta$ times the critical density of the Universe at the redshift of the source., adopting $r_{500}$ and $r_{200}$ values found in the literature and assuming the conversion r_{500}\leavevmode\nobreak\ \simeq\leavevmode\nobreak\0.65$r_{200}$ of Reiprich et al. (2013). For our cosmology and redshifts we use the approximation of Henry et al. (2009):

[TABLE]

where $E(z)$=$(\Omega_{\rm{M}}(1+z)^{3}+1-\Omega_{\rm{M}})^{1/2}$. The mean temperature values, $\left<kT\right>$, are adopted from references as described in Table 2. Our averaged Fe abundance distribution, stacked from all the above measurements, is shown in Fig. 14 (right), and the numerical values are reported in Table 3.

We use the stacking method of Mernier et al. (2017) to determine this distribution. The averaged abundance distribution, $Z_{\text{ref}}(k)$, as a function of the $k$th reference bin is defined as:

[TABLE]

where $Z(i)_{j}$ is the individual abundance of the $j$th cluster at its $i$th region; $\sigma_{Z(i)_{j}}$ is its statistical error, $N$ is the number of clusters, $M$ is the number of regions analysed in each cluster and $w_{i,j,k}$ is the weighting factor. This factor represents the geometrical overlap of the $i$th region with the $k$th reference bin in cluster $j$. It varies from 0 to 1. The stacked statistical error is calculated as:

[TABLE]

We also obtained the scatter of the measurements for each $k$th reference bin as:

[TABLE]

It clearly appears from Fig. 14 that abundance distributions of these merging clusters are not flat. Instead, we note a decrease from the core up to $\sim$0.4$r_{500}$, followed by a flattening at larger radii. Figure 14 (right) compares our averaged distribution with the relaxed, azimuthally-averaged CC systems of Mernier et al. (2017) with $kT$ ¿ 1.7 keV and the disturbed, NCC systems of Lovisari & Reiprich (2019). In both cases we have included the statistical errors and the scatter. Remarkably, our values are in good agreement with the (azimuthally-averaged) ”NCC” profile of Lovisari & Reiprich (2019), which shows lower abundance values in the core than relaxed clusters and a shallower decrease towards the outer radii. The averaged Fe distribution values shown in Table 3 are therefore consistent with Lovisari & Reiprich (2019) results. We have additionally checked the averaged abundance distribution along the BCG displacement direction (assuming only the positive radii). In this case, the agreement with the NCC systems profile of Lovisari & Reiprich (2019) improves. For more details see Appendix A.

The lower abundance value in the core for disturbed systems is suggested by Lovisari & Reiprich (2019) to be caused by major mergers, where in some cases the core could be disrupted. As a result, the abundance peak in the center of the (CC) progenitors is spread out and distributed across larger radii. At first glance, the similarity of the Fe distribution observed between NCC clusters and these merging clusters tends to support this scenario. However, the steeper Fe peaks observed in the core of 1RXSJ0603 and in the vicinity of the core of A665 (perhaps offset from the BCG because of important sloshing motions) suggest that remnants of previous CC clusters can partly persist to (or, at least, not be entirely disrupted by) even the most major mergers in some cases (see Sect. 5.2). While this moderate central abundance increase could be associated with either (relatively old) core passage mergers, or on the contrary, early stage mergers, we note that our sample contains both cases. In fact, most of them are massive clusters, which have undergone a major merger with a core passage more than 0.6 Gyr ago. However, A2256 is a still on going merger, before core passage, although it is thought that the main cluster center is already affected by an older merger. Future larger samples, containing an equal proportion of early and late stage mergers, will help to relate the (re)-distribution of metals with merger history.

5.2 Abundance vs. pseudo-entropy relation

In order to study the thermal history of the ICM, we have calculated the pseudo-entropy, $K\equiv kT\times n^{-2/3}_{e}$, as a function of the region distance for each merging cluster. These pseudo-entropy distributions along the merging axis are shown together with the Fe abundance profiles in Fig. 15. The temperature of each region, $kT$, is directly obtained from the spectral analysis, and the electron density, $n_{e}$, can be inferred from the emission measure, $Norm$, as $Norm$ $\propto$ 1.2$n^{2}_{e}V$, where $V$ is the emitting volume projected onto the line-of-sight (LOS). To estimate $n_{e}$, we assume a spherical LOS in which only the sphere between the maximum and minimum radius contributes to the emission (Henry et al. 2004; Mahdavi et al. 2005). We assume the emission volume $V$ to be $V=2SL/3$ for the circular and elliptical sectors, and $V=\pi SL/2$ for rectangular ones. In each case, $L=2\sqrt{(R_{max}^{2}-R_{min}^{2})}$ and $S$ is the area of the sectors in the plane of the sky.

Previous studies (De Grandi et al. 2004; Leccardi et al. 2010; Ghizzardi et al. 2014; Liu et al. 2018) have demonstrated, especially for relaxed (CC) clusters, the existence of a negative correlation between the Fe abundance and the gas entropy. Low-entropy cores, where the cooling of the ICM might take place, are in fact associated with central abundance peaks. However, disturbed (NCC) clusters are not well characterised by this behaviour as mentioned by Leccardi et al. (2010) and Rossetti & Molendi (2010).

Pseudo-entropy quantifies the internal energy variation due to thermal processes such as, e.g. cooling, turbulent dissipation or shock heating. Regions under the influence of the latter, in particular, can be easily identified because they exhibit a positive gradient proportional to $r^{1.1}$ (Tozzi & Norman 2001). In particular for our merging cluster selection, we observe a similar trend in the post-shock regions, see Fig. 15 and Fig. 16 (left), although some systems (e.g. 1RXSJ0603, A3376) significantly deviate from this relation. If we fit simultaneously all our pseudo-entropy distributions with a power-law ($\sim r^{\alpha}$), we obtain a slope of $\alpha$ = 1.34 $\pm$ 0.63 (Fig. 16 left). In the case of 1RXSJ0603, we see a flat profile between $\sim$0.2–0.4 $r_{500}$, coincident with the COREN location, possibly indicating a sign of the subcluster core remnant after the merger event. A similar flattening is observed in A3376, although its origin may be rather related to the central low-entropy gas stripped during the merger, thereby forming the cometary shape ICM.

Most of NCC clusters are thought to have undergone a CC phase, before the (cool) core being disrupted completely or partially by major mergers (Leccardi et al. 2010; Rossetti & Molendi 2010; ZuHone 2011). The on-going interactions heat and increase the entropy in the core. Therefore, NCC tend to have high-entropy cores and smaller entropy gradients. Some of these clusters present regions with lower entropy that have associated an abundance excess, in most of the cases associated with the BCG or a giant elliptical galaxy, and known as ”CC-remnants” (Rossetti & Molendi 2010). It is suggested that these regions are the reminiscent of CC after the merging activity. Another consequence of merger events is the decrease of the abundance excess in the core, which can be erased by the shock-heating and gas mixing processes.

Figure 15 shows simultaneously the pseudo-entropy and Fe abundance distributions scaled by $r_{500}$ for our six merging galaxy clusters. In the case of CIZA2242 and A3376, the entropy is slightly flat in the core and the profile increases smoothly (though with a larger gradient) out to the recently shocked region. Interestingly, 1RXSJ0603 exhibits a relative low-entropy core associated with the abundance peak. One natural interpretation is that this region is the remainder of cool core after the merging activity, which agrees with ZuHone (2011) simulations taking into account the mass ratio of 3:1 and it is found in 12 clusters of Rossetti & Molendi (2010). An alternative scenario could be that 1RXSJ0603, as a well evolved merger ($\sim$2 Gyr since the shock producing the radio relic; e.g. Brüggen et al. (2012), Stroe et al. (2015a)), has started to recover its potential well and to re-build its entropy (and metal) stratification after the merger. For A3667 and A2256, we identify possible CC-remnants characterised by somewhat lower entropy and higher metallicity values (Rossetti & Molendi 2010). These remnants clearly appear after we include the measurements of the truncated yellow region after the A3667 cold front (see Fig. 9) and the CORE2 values of A2256 (see Fig. 13), in dashed lines in Fig. 15. A2256 also contains a second minimum corresponding to the CORE1. As explained in Leccardi et al. (2010), two low-entropy clumps can be found if we are in the presence of an early stage merger, as in this case, or an off-axis merger event. A665 also shows a relative low-entropy drop in the core, but in this case the abundance peak is displaced towards the cold front, as explained in Section 5.1, possibly due to the effect of the internal dynamics of the cold front during the infall of the subcluster. The above examples, and in particular the case of 1RXSJ0603, strongly suggest that in some cases, signatures of cool cores can remain even after the most major ICM merger events.

Furthermore, we try to determine a correlation between Fe abundance and pseudo-entropy for our merging clusters. Figure 16 (right) shows the distribution of $Z$-$K$. Generally speaking, our measured entropies are clearly higher than in the distribution of Liu et al. (2018), with our lowest value being $\sim$200 keV cm2 and the maximum abundance value being $\sim$0.5 $Z_{\odot}$. In fact, these values are characteristic of merging clusters (Leccardi et al. 2010; ZuHone 2011). We fit the distribution with a linear function $Z\sim Z_{0}-\alpha\,K$/1000. Our best-fit parameters are $Z_{0}=0.51\pm 0.03$ and $\alpha=0.32\pm 0.06$. This correlation shows a flatter slope than obtained by Liu et al. (2018) for relaxed (CC) galaxy clusters ($Z_{0}=0.86\pm 0.18$ and $\alpha=1.49\pm 0.51$). Our results suggest that the abundance and pseudo-entropy relation are not always tightly related. Instead, it is likely that, when a merger occurs, the entropy gets boosted in post-shock regions more rapidly than metals diffuse out. Interestingly, when we consider the values of pseudo-entropy only below $\sim$500 keV cm2 – which were presumably less affected by shock-heating, or which start to relax again (see vertical black dashed line in Fig. 16), the slope of the relation considerably increases (Table 4). This might suggest that in the lowest entropy regions ($\sim$200–500 keV cm2) of these merging clusters (NCC) – usually associated with the core remnant (see Fig. 15), the metal-rich gas (originally still coming from the BCG) gets mixed and redistributed at a moderate rate compared to the entropy changes during the merger event (ZuHone 2011). Alternatively, chemical histories may differ between relaxed and merging systems. Although beyond the scope of this paper, in some cases we cannot rule out the possibility of an extra Fe enrichment, produced by a boost of star formation or SNIa explosions, triggered by major mergers (e.g. Stroe et al. (2015a); Sobral et al. (2015)).

Finally, we do not find evidence of a correlation in the higher (¿500 keV cm2) entropy range – i.e. beyond $\sim$0.4–0.6 $r_{500}$. In these outer regions, the Fe abundance remains rather uniform around 0.2–0.3 Z⊙, with no notable trend with cluster mass, radius or azimuth (Fig. 14). Combined with the previous observations of a uniform $\sim$0.3 Z⊙ Fe abundance in the outskirts of (mostly) CC clusters (Fujita et al. 2008; Werner et al. 2013; Thölken et al. 2016; Urban et al. 2017; Ezer et al. 2017), we conclude that Fe distribution in cluster outskirts of these merging clusters depends very weakly (or not at all) on the dynamical state of the cluster. This has important consequences regarding the chemical history of large-scale structures. In fact, although these previous measurements in relaxed systems strongly suggested that the most of the ICM enrichment took place in early stages of cluster formation (i.e. at $z$ ¿ 2–3; for recent reviews, see Mernier et al. 2018b; Biffi et al. 2018b), the alternative scenario of late enrichment coupled with a very efficient mixing of metals in the ICM, could not be totally excluded. Among the entire cluster population, merging clusters are systems in which metals had the least time to mix and diffuse out since their last merger event. Consequently, if the bulk of metals were mixed efficiently in the entire cluster volume (after being injected into the ICM at late stages), one would expect to find important variations of the Fe abundance at large radii in such dynamically disturbed systems, with values considerably higher than 0.2–0.3 Z⊙ in some specific outer regions. Instead, the remarkable similarity of the Fe ndistribution in the outskirts of CC and these six merging clusters, and the above indication that in inner regions metal mixing lags somewhat behind entropy changes, constitute valuable additional evidence in favour of the ”early enrichment” scenario.

6 Summary

In this work, we have analysed for the first time the Fe distribution along the merging axis of six merging galaxy clusters ($0.05<z<0.23$) observed with XMM-Newton/EPIC. We have obtained and specifically studied the temperature, abundance and pseudo-entropy distributions across the ICM of these disturbed systems. Specifically, we have provided a first estimate of the averaged Fe distribution out to $r_{500}$ – as well as the relation between Fe abundance and pseudo-entropy – in six merging clusters. The main results can be summarised as follows.

• $\bullet$

The Fe distribution is far from being universal, as it strongly varies from one merging cluster to another. While in some cases (e.g. CIZA2242, A2256) no Fe peak nor pseudo-entropy drop is detected in central regions, some of our observations (e.g. 1RXSJ0603, A665) strongly suggest that remnants of previous (metal-rich) cool-cores can persist after major mergers.

• $\bullet$

Our averaged Fe distribution in these six merging clusters is in excellent agreement with that in (non-merging) ”non-cool-core” clusters presented by Lovisari & Reiprich (2019). The average profile clearly shows a moderate abundance central excess, thus present not only in CC systems but even in many spectacular mergers. Compared to the former, however, the central ($\lesssim 0.3r_{500}$) Fe profile in these merging clusters is significantly less steep. This can be interpreted as the metal-rich gas of the CC progenitors being mixed and spread across moderate radii by merging processes.

• $\bullet$

The minimum pseudo-entropy value obtained for these six merging clusters is $\sim$200 keV cm2. This value is characteristic of typical unrelaxed systems, whose internal activity is dominated by turbulence, large-scale motions and/or AGN feedback. Whereas some clusters exhibit a pseudo-entropy distribution that tends to follow the behaviour $\propto$ $r^{1.1}$ described by Tozzi & Norman (2001), some systems (e.g. 1RXSJ0603, A3376) deviate from this relation due to stripped cold gas from disrupted cores.

• $\bullet$

We find a mild negative correlation between abundance and pseudo-entropy for the inner core region ($K$ = 200-500 keV cm2), shallower than the anti-correlation obtained by Liu et al. (2018) for relaxed clusters. This might suggest that mergers affect the gas entropy more rapidly than they re-distribute metals (via gas mixing).

• $\bullet$

In the outermost, high-entropy regions ($K$ = 500-2300 keV cm2) studied here, we have no evidence of a correlation between metallicity and pseudo-entropy and the abundance is consistent with the uniform abundance value of 0.2–0.3 $Z_{\odot}$. This demonstrates that the lack of dependence of Fe abundance at such high entropies (and large radii) applies not only to the mass, radius and azimuthal characteristics of clusters, but also to the dynamical stage. In other words, our results provide important additional evidence towards the scenario of an early enrichment of the inter galactic medium ($z$ ¿2–3).

Acknowledgements.

The authors thank Dr. R. van Weeren for the LOFAR radio data, Dr. R. Kale, Dr. T. Shimwell and Dr. V. Vacca for providing the VLA radio data and Dr. A. Botteon for the SUMMSS data. S. Kara and. E. N. Ercan would like to thank TUBITAK for financial support through a 1002 project under the code 118F023. SRON is supported financially by NWO, the Netherlands Organization for Scientific Research. This work is based on observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA member states and the USA (NASA).

————————————————–

Appendix A Averaged abundance distribution along BCG displacement direction

In addition to stacking our measurements over the entire merger axis as explained in Section 5.1, we can alternatively choose to trace the ICM only from the BCG region (i.e. typically including the cluster X-ray peak), down to the freshly collided gas, i.e. along its cometary structure (defined in Figs. 3, 5, 7, 9, 11 and 13 as positive values on the x-axis). These regions better trace the disturbed and mixed plasma while avoiding special features such as cold fronts ahead of the merger (e.g. A3667, Fig. 9); therefore the effects of the merging activity in the ICM are expected to be more significant. For CIZA2242, 1RXSJ0603 and A3376 this covers the complete elongation or cometary tail present in the ICM. For A3367 and A665 this includes the sectors 5–8 and 4–9, respectively. In the case of A2256, we consider the sectors 4–7, where the ICM of the main subclusters is predominantly affected by the presence of the infalling subclsuter. The individual Fe abundance distribution of the six merging galaxy clusters included in this analysis are plotted in Fig. 17 (top), the averaged values obtained with the stacking method described in Section 5.1 are shown in Fig. 17 (bottom) and Table 5.

Interestingly, the Fe averaged distribution is now in excellent agreement with the Fe average profile of NCC systems by Lovisari & Reiprich (2019). The major differences with respect to the distribution shown in Section 5.1 appear in the 0.1–0.2 $r_{500}$ (8$\%$ lower) and 0.2–0.3 $r_{500}$ (10$\%$ lower) radii. These lower values are obtained because the displaced Fe peak outside cluster center (see Section 5.1 and 5.2 for more details) of A3667 and A665 are not included.

Appendix B EPIC background modelling for the cluster sample

The following tables summarise the EPIC background model parameters used of each one of the galaxy clusters using the method described in Section 3.2. The norms are given in standard SPEX definition, which are intrinsically dependent of the distance to the source.

Appendix C Temperature and abundance analysis for the cluster sample

The following tables summarise the best-fit parameters for the temperature and abundance distribution along the merging axis.