Diffuse Radio Emission from Galaxy Clusters

TL;DR

This paper reviews the properties and recent observational findings of diffuse radio sources in galaxy clusters, discussing their classifications, physical processes, and implications for intracluster medium physics, including cosmic web emission detection.

Contribution

It provides a comprehensive overview of diffuse cluster radio sources, including recent observational results and an updated catalog, highlighting advances in understanding ICM magnetic fields and acceleration processes.

Findings

Diffuse radio sources reveal cosmic rays and magnetic fields in ICM.

Classification of sources into halos, relics, and revived fossil plasma.

Recent observational data support models of shock acceleration and fossil plasma revival.

Abstract

In a growing number of galaxy clusters diffuse extended radio sources have been found. These sources are not directly associated with individual cluster galaxies. The radio emission reveal the presence of cosmic rays and magnetic fields in the intracluster medium (ICM). We classify diffuse cluster radio sources into radio halos, cluster radio shocks (relics), and revived AGN fossil plasma sources. Radio halo sources can be further divided into giant halos, mini-halos, and possible `intermediate' sources. Halos are generally positioned at cluster center and their brightness approximately follows the distribution of the thermal ICM. Cluster radio shocks (relics) are polarized sources mostly found in the cluster's periphery. They trace merger induced shock waves. Revived fossil plasma sources are characterized by their radio steep-spectra and often irregular morphologies. In this review we give an overview of the properties of diffuse cluster radio sources, with an emphasis on recent observational results. We discuss the resulting implications for the underlying physical acceleration processes that operate in the ICM, the role of relativistic fossil plasma, and the properties of ICM shocks and magnetic fields. We also compile an updated list of diffuse cluster radio sources which will be available on-line http://galaxyclusters.com. We end this review with a discussion on the detection of diffuse radio emission from the cosmic web.

Figures

Tables

Table 1. Table 1: A table of clusters which show evidence for shock waves in X-ray observations and that coincide with the location of cluster radio shocks (Akamatsu et al. in prep.)

Name	$T$ jump	$\rho$ jump	Spec index	X-ray ref	Radio ref
Coma SW	✔		✔	Akamatsu et al. (2013)	Giovannini et al. (1991)
				Ogrean & Brüggen (2013)	Thierbach et al. (2003)
Abell 115	✔	✔	✔	Botteon et al. (2016a)	Govoni et al. (2001b)
Abell 754	✔	✔	✔	Krivonos et al. (2003)	Kassim et al. (2001)
				Macario et al. (2011)	Kale & Dwarakanath (2009)
Abell 1240		✔	✔	Hoang et al. (2018b)	Kempner & Sarazin (2001)
					Bonafede et al. (2009a)
					Hoang et al. (2018b)
Abell 3667 NW	✔	✔	✔	Finoguenov et al. (2010)	Röttgering et al. (1997)
				Akamatsu et al. (2012a)	Johnston-Hollitt (2003)
				Sarazin et al. (2016)	Hindson et al. (2014)
					Johnston-Hollitt & Pratley (2017)
Abell 3667 SE	✔	✔	✔	Akamatsu & Kawahara (2013)	Johnston-Hollitt (2003)
				Storm et al. (2018)	Hindson et al. (2014)
					Riseley et al. (2015)
Abell 3376 W	✔	✔	✔	Akamatsu et al. (2012b)	Bagchi et al. (2006)
				Urdampilleta et al. (2018)	Kale et al. (2012)
					George et al. (2015)
Abell 3376 E	✔	✔	✔	Urdampilleta et al. (2018)	Kale et al. (2012)
Abell 2255	✔	✔	✔	Akamatsu et al. (2017)	Pizzo & de Bruyn (2009)
Abell 2256	✔		✔	Trasatti et al. (2015)	Röttgering et al. (1994)
					Clarke & Enßlin (2006)
					van Weeren et al. (2012c)
Abell 2744	✔	✔	✔	Eckert et al. (2016)
				Hattori et al. (2017)	Pearce et al. (2017)
Sausage N	✔		✔	Akamatsu & Kawahara (2013)	van Weeren et al. (2010)
				Ogrean et al. (2013a)	Stroe et al. (2013, 2014b, 2016)
				Ogrean et al. (2014a)	Hoang et al. (2017)
				Akamatsu et al. (2015)	Loi et al. 2017 Loi et al. (2017)
					Kierdorf et al. (2017)
					Di Gennaro et al. (2018)
Sausage S	✔		✔	Akamatsu et al. (2015)	Hoang et al. (2017)
					Stroe et al. (2013)
					Di Gennaro et al. (2018)
Toothbrush	✔	✔	✔	Ogrean et al. (2013b)	Rajpurohit et al. (2018)
				Itahana et al. (2015)	van Weeren et al. (2012b)
				van Weeren et al. (2016)	van Weeren et al. (2016)
					Stroe et al. (2016)
El Gordo	✔	✔	✔	Botteon et al. (2016b)	Lindner et al. (2014)
				Basu et al. (2016a)
ZwCl 2341.1+0000		✔	✔	Ogrean et al. (2014b)	Bagchi et al. (2002)
					van Weeren et al. (2009d)
					Giovannini et al. (2010)
					Benson et al. (2017)
Bullet reverse	✔	✔	✔	Shimwell et al. (2015)	Shimwell et al. (2015)
Abell 2146	✔	✔		Russell et al. (2010, 2011, 2012)	Hlavacek-Larrondo et al. (2018)
Abell 521	✔	✔	✔	Bourdin et al. (2013)	Brunetti et al. (2008)
					Giacintucci et al. (2008)
RXJ1314.4-2515	✔	✔		Mazzotta et al. (2011)	Feretti et al. (2005)
					Venturi et al. (2007)

Table 2. Table 2: Clusters with diffuse radio emission

cluster	z	classification	notes	references
Abell S753	0.014	RS/F		Goss et al. (1987); Subrahmanyan et al. (2003)
Perseus cluster	0.018	MH		Miley & Perola (1975); Noordam & de Bruyn (1982)
				Pedlar et al. (1990); Burns et al. (1992); Sijbring (1993)
				Sijbring & de Bruyn (1998); Brentjens (2011)
				Gendron-Marsolais et al. (2017)
Abell 1367	0.022	RS		Gavazzi (1978); Hanisch (1980); Ballarati et al. (1981)
				Gavazzi & Trinchieri (1983); Gavazzi & Jaffe (1987)
				Farnsworth et al. (2013)
Coma clustera	0.023	H, RS		Large et al. (1959); Willson (1970); Jaffe & Rudnick (1979)
				Ballarati et al. (1981); Andernach et al. (1984)
				Giovannini et al. (1985); Schlickeiser et al. (1987)
				Giovannini et al. (1991, 1993)
				Deiss et al. (1997); Thierbach et al. (2003); Pizzo (2010)
				Brown & Rudnick (2011)
Ophiuchus cluster	0.028	MH, F		Murgia et al. (2009); Govoni et al. (2009); Murgia et al. (2010a)
				Werner et al. (2016)
Abell 4038	0.030	F		Slee & Reynolds (1984); Slee et al. (1994); Slee & Roy (1998)
				Slee et al. (2001); Kale & Dwarakanath (2012); Kale et al. (2018)
				Duchesne et al. (2017)
2A 0335+096b	0.036	MH		Sarazin et al. (1995)
Abell 1314	0.034	F		Rudnick & Lemmerman (2009); Wilber et al. (2018a)
Abell 548b	0.045	RS/F		Giovannini et al. (1999, 2006); Feretti et al. (2006)
				George et al. (2017)
Abell 168	0.045	mRS		Dwarakanath et al. (2018)
Abell 3376	0.046	dRS		Bagchi et al. (2006); Kale et al. (2012); George et al. (2015)
Abell 1213	0.047	cH		Giovannini et al. (2009)
Abell 3562	0.049	H		Venturi et al. (2000, 2003); Giacintucci et al. (2005)
Abell 754	0.054	H, RS/F		Wielebinski et al. (1977); Mills et al. (1978); Harris et al. (1980a)
				Giovannini et al. (1999); Kassim et al. (2001)
				Bacchi et al. (2003); Kale & Dwarakanath (2009)
				Macario et al. (2011)
Abell 85	0.055	F		Slee & Reynolds (1984); Slee et al. (1994)
				Giovannini & Feretti (2000); Duchesne et al. (2017)
Abell 2626	0.055	U	previously classified as MH	Gitti et al. (2004); Gitti (2013); Ignesti et al. (2017)
				Kale & Gitti (2017)
Abell 3667	0.056	dRS	non-confirmed MH	Schilizzi & McAdam (1975); Goss et al. (1982)
				Jones & McAdam (1992); Röttgering et al. (1997)
				Johnston-Hollitt (2003); Riseley et al. (2015)
				Hindson et al. (2014); Johnston-Hollitt & Pratley (2017)
Abell 2319	0.056	H		Harris & Miley (1978); Feretti et al. (1997b)
				Giovannini et al. (1999); Kempner & Sarazin (2001)
				Farnsworth et al. (2013); Storm et al. (2015)
Abell 133	0.057	cF		Slee & Reynolds (1984); Slee et al. (2001, 1994, 2001)
				Giovannini et al. (1999); Randall et al. (2010)
				Duchesne et al. (2017)
Abell 2256	0.058	H, RS, mF		Bridle & Fomalont (1976); Masson & Mayer (1978)
				Bridle et al. (1979); Röttgering et al. (1994); Kim (1999)
				Giovannini et al. (1999); Kempner & Sarazin (2001)
				Clarke & Enßlin (2006); Brentjens (2008)
				van Weeren et al. (2009a); Rudnick & Lemmerman (2009)
				Kale & Dwarakanath (2010); van Weeren et al. (2012c)
				Owen et al. (2014); Trasatti et al. (2015); Ozawa et al. (2015)
RXC J0225.1–2928	0.060	cRS		Shakouri et al. (2016)
Abell 1795	0.062	cMH		Giacintucci et al. (2014b)
CIZA J0649.3+1801	0.064	cRS		van Weeren et al. (2011b)
CL 0217$+$70	0.065c	dRS, H		Brown et al. (2011a)
RXC J1053.7+5452	0.070	RS		Rudnick & Lemmerman (2009); van Weeren et al. (2011b)
Abell 399	0.072	H	close to Abell 401	Murgia et al. (2010b)
Abell 2065	0.072	cH		Farnsworth et al. (2013)
Abell 401	0.074	H	close to Abell 399	Harris et al. (1980b); Roland et al. (1981); Giovannini et al. (1999)
				Giovannini & Feretti (2000); Bacchi et al. (2003)
				Murgia et al. (2010b)
Abell 2067	0.074	cR		Farnsworth et al. (2013)
ZwCl 1742.1+3306	0.076	cMH		Giacintucci et al. (2014b)
Abell 2029	0.077	MH		Murgia et al. (2009); Govoni et al. (2009)
Abell 2061	0.078	RS, cH		Kempner & Sarazin (2001); Rudnick & Lemmerman (2009)
				van Weeren et al. (2011b); Farnsworth et al. (2013)
Abell 2255	0.081	H, cRS, mF	cUSS	Jaffe & Rudnick (1979); Harris et al. (1980b); Burns et al. (1995)
				Feretti et al. (1997a); Giovannini et al. (1999)
				Kempner & Sarazin (2001); Govoni et al. (2005)
				Pizzo et al. (2008); Pizzo & de Bruyn (2009)
				Rudnick & Lemmerman (2009); Pizzo et al. (2011)
Abell 2556	0.087	cF	cUSS	Duchesne et al. (2017)
Abell 478	0.088	MH		Giacintucci et al. (2014b); Savini et al. (2018a)
Abell 2142	0.089	H, cMH	H-MH	Giovannini et al. (1999); Giovannini & Feretti (2000)
				Farnsworth et al. (2013); Venturi et al. (2017)
Abell 725	0.092	cRS		Kempner & Sarazin (2001)
Abell 3365	0.093	dcRS	at least one RS present	van Weeren et al. (2011b)
Abell 13	0.094	F		Slee & Reynolds (1984); Slee et al. (1994)
				Giovannini et al. (1999)
				Slee et al. (2001); George et al. (2017)
				Duchesne et al. (2017)
Abell 610	0.095	RS		Giovannini & Feretti (2000)
Abell 2048	0.097	F		van Weeren et al. (2009b, 2011d)
PKS 0745–191	0.103	cMH		Baum & O’Dea (1991); Gitti et al. (2004)
				Venturi et al. (2007)
ZwCl 0008.8+5215	0.103	dRS		van Weeren et al. (2011c); Kierdorf et al. (2017)
				Golovich et al. (2017a)
Abell 523	0.104	cH/cRS		Giovannini et al. (2011); van Weeren et al. (2011b)
				Girardi et al. (2016)
Abell 2798	0.105	RS		Duchesne et al. (2017)
CIZA J0107.7+5408	0.107	cDR, cH, cF	diffuse emission present	van Weeren et al. (2011b); Randall et al. (2016)
Abell 2751	0.107	RS		Duchesne et al. (2017)
Abell S0084	0.108	cH		Duchesne et al. (2017)
Abell 2443	0.108	F/RS	complex	Cohen & Clarke (2011)
Abell 2811	0.108	H		Duchesne et al. (2017)
Abell 2034	0.113	cF, cRS, cH	diffuse emission present	Kempner & Sarazin (2001); Giovannini et al. (2009)
				van Weeren et al. (2011b); Rudnick & Lemmerman (2009)
				Shimwell et al. (2016)
Abell 2721	0.114	cRS		Duchesne et al. (2017)
Abell 2069	0.115	H		Farnsworth et al. (2013); Drabent et al. (2015)
Abell 2496	0.122	cRS		Duchesne et al. (2017)
Abell 1664	0.128	RS/F		Giovannini et al. (1999); Govoni et al. (2001b)
				Kale & Dwarakanath (2012)
Abell 1033	0.130	F, GReET		Rudnick & Lemmerman (2009); de Gasperin et al. (2015b)
				de Gasperin et al. (2017a)
Abell 1132	0.137	H	cUSS	Rudnick & Lemmerman (2009); Wilber et al. (2018b)
Abell 1068	0.138	cMH		Govoni et al. (2009)
Abell 22	0.142	cH/cRS		Duchesne et al. (2017)
Abell 1413	0.143	MH		Govoni et al. (2009); Savini et al. (2018a)
24P73	0.150d	F		van Weeren et al. (2009b, 2011d)
Abell 3888	0.151	H		Shakouri et al. (2016)
Abell 2204	0.152	MH		Giacintucci et al. (2014b)
Abell 545	0.154	H		Giovannini et al. (1999); Bacchi et al. (2003)
Abell 1240	0.159	dRS		Kempner & Sarazin (2001); Bonafede et al. (2009a)
				Hoang et al. (2018b)
WHL J143150.1+133205	0.160	F		van Weeren et al. (2009b, 2011d); Ogrean et al. (2011)
				Shulevski et al. (2015)
Abell 3411–3412	0.162	H, mRS		van Weeren et al. (2013); Giovannini et al. (2013)
				van Weeren et al. (2017b)
RXC J1720.1+2637	0.164	MH	H-MH	Mazzotta & Giacintucci (2008); Giacintucci et al. (2014a)
				Savini et al. (2018a)
Abell 2294	0.169	cH		Owen et al. (1999); Giovannini et al. (2009)
Abell 907	0.167	MH		Giacintucci et al. (2017)
Abell 1914	0.171	F, cH		Roland et al. (1985); Giovannini et al. (1999)
				Kempner & Sarazin (2001); Bacchi et al. (2003)
				Rudnick & Lemmerman (2009); Mandal et al. (2018)
Abell 2218	0.171	H		Giovannini et al. (1999); Giovannini & Feretti (2000)
				Kempner & Sarazin (2001); Rudnick & Lemmerman (2009)
Abell 2073	0.172	cR		Farnsworth et al. (2013)
Abell 2693	0.173	cH		Duchesne et al. (2017)
ZwCl 0634.1+4750	0.174	H		Cuciti et al. (2018)
Abell 2680	0.177	cH		Duchesne et al. (2017)
Abell 2345	0.177	dRS		Giovannini et al. (1999); Bonafede et al. (2009a)
				George et al. (2017)
Abell 1931	0.178	cF		Brüggen et al. (2018)
Abell  2254	0.178	H		Giovannini et al. (1999); Govoni et al. (2001b)
				George et al. (2017)
Abell 1612	0.179	RS/F		van Weeren et al. (2011b); Kierdorf et al. (2017)
Abell 665	0.182	H		Giovannini et al. (1999); Giovannini & Feretti (2000)
				Kempner & Sarazin (2001); Feretti et al. (2004b)
				Vacca et al. (2010); Rudnick & Lemmerman (2009)
Abell 1689	0.183	H		Vacca et al. (2011)
CIZA J2242.8+5301e	0.192	F, mdRS, H		van Weeren et al. (2010); Stroe et al. (2013, 2016)
				Kierdorf et al. (2017); Loi et al. (2017)
				Di Gennaro et al. (2018)
Abell 115	0.197	RS		Giovannini et al. (1999); Govoni et al. (2001b)
				Botteon et al. (2016a); Hallman et al. (2018)
Abell 1451	0.199	H		Cuciti et al. (2018)
Abell 3527-bis	0.200	RS		de Gasperin et al. (2017b)
Abell 2163	0.203	H, RS		Herbig & Birkinshaw (1994); Giovannini et al. (1999)
				Feretti et al. (2001, 2004b); George et al. (2017)
				Thölken et al. (2018)
Abell 520	0.203	H, cmRS		Giovannini et al. (1999); Govoni et al. (2001b)
				Bacchi et al. (2003); Vacca et al. (2014)
				Wang et al. (2018); Hoang et al. (2018c)
Abell 910	0.206	RS		Govoni et al. (2012)
Abell 209	0.206	H		Giovannini et al. (1999, 2006); Venturi et al. (2007)
				Giovannini et al. (2009); Venturi et al. (2013)
RXC J1504.1–0248	0.215	MH		Giacintucci et al. (2011b)
Abell 773	0.217	H		Giovannini et al. (1999); Kempner & Sarazin (2001)
				Govoni et al. (2001b); Rudnick & Lemmerman (2009)
PLCK G200.9–28.2	0.220	RS		Kale et al. (2017)
Abell 800	0.222	H		Govoni et al. (2012)
RXC J1514.9–1523	0.223	H	USS	Giacintucci et al. (2011a)
Abell 2261	0.224	H		Venturi et al. (2008); Sommer et al. (2017)
				Burke-Spolaor et al. (2017); Savini et al. (2018a)
[VRI2012] Toothbrushf	0.225	mRS, cF, H		van Weeren et al. (2012b); Stroe et al. (2016)
				van Weeren et al. (2016); Kierdorf et al. (2017)
				Rajpurohit et al. (2018)
Abell 2667	0.226	MH		Giacintucci et al. (2017)
Abell 1682	0.226	cH, cRS, cmFS	diffuse emission present	Venturi et al. (2008); Macario et al. (2013)
				Venturi et al. (2013); Rudnick & Lemmerman (2009)
Abell 2219	0.228	H		Giovannini et al. (1999); Kempner & Sarazin (2001)
				Bacchi et al. (2003); Orrú et al. (2007)
RXC J1234.2+0947g	0.229	cH, RS		Kale et al. (2015)
Abell 141	0.230	H	USS	Duchesne et al. (2017)
Abell 2146	0.232	cH, cDRS	diffuse emission present	Hlavacek-Larrondo et al. (2018)
				Hoang et al. (2018a)
Abell 746	0.232	H, RS		van Weeren et al. (2011b)
Abell 2390	0.233	…	non-confirmed H	Giovannini et al. (1999); Bacchi et al. (2003)
				Sommer et al. (2017); Savini et al. (2018a)
Abell 33	0.234	cRS		Duchesne et al. (2017)
RX J2129.6+0005	0.235	MH		Kale et al. (2015)
Abell S780	0.236	MH		Kale et al. (2015)
RXC J1314.4–2515	0.247	H, dRS		Feretti et al. (2005); Venturi et al. (2007, 2013)
				George et al. (2017)
Abell 521	0.248	H, RS	USS	Ferrari et al. (2006); Venturi et al. (2007)
				Dallacasa et al. (2009); Giovannini et al. (2009)
				Brunetti et al. (2008); Giacintucci et al. (2008)
				Venturi et al. (2013); Macario et al. (2013)
RXC J2351.0–1954	0.248	cdRS, cH		Duchesne et al. (2017)
Abell 3444	0.253	MH		Venturi et al. (2007); Giovannini et al. (2009)
				Kale et al. (2015)
Abell 1835	0.253	MH		Murgia et al. (2009); Govoni et al. (2009)
Abell 1550	0.254	H		Govoni et al. (2012)
ZwCl 1454.8+2233h	0.258	MH		Venturi et al. (2008)
CIZAJ1938.3+5409	0.260	cH		Bonafede et al. (2015a)
PSZ1 G139.61+24.20	0.267	MH	H-MH	Giacintucci et al. (2017); Savini et al. (2018b)
Abell 1443	0.269	cH, cRS		Bonafede et al. (2015a)
PSZ1 G171.96–40.64	0.270	H	cUSS	Giacintucci et al. (2013)
ZwCl 2341.1+0000	0.270	dRS, cH		Bagchi et al. (2002); van Weeren et al. (2009d)
				Giovannini et al. (2010); Benson et al. (2017)
Abell 1758Ni	0.280	H	close to Abell 1758S	Giovannini et al. (1999, 2006, 2009)
				Kempner & Sarazin (2001)
				Rudnick & Lemmerman (2009); Venturi et al. (2013)
				Botteon et al. (2018a)
Abell 1758Sj	0.280	H, cRS	close to Abell 1758N	Botteon et al. (2018a)
Abell 697	0.282	H	USS	Kempner & Sarazin (2001); Venturi et al. (2008)
				Giovannini et al. (2009)
				Rudnick & Lemmerman (2009)
				Macario et al. (2010); van Weeren et al. (2011b)
				Venturi et al. (2013); Macario et al. (2013)
ZwCl 1021.0+0426k	0.291	MH		Giacintucci et al. (2014b); Kale et al. (2015)
RXC J1501.3+4220	0.292	H		Wilber et al. (2018a)
Bullet clusterl	0.296	H, RS		Liang et al. (2000); Shimwell et al. (2014, 2015)
Abell 781	0.298	cRS/cF	non-confirmed H	Venturi et al. (2008); Rudnick & Lemmerman (2009)
				Govoni et al. (2011); Venturi et al. (2011, 2013)
				Botteon et al. (2018b)
PSZ1 G096.89+24.17m	0.300	dRS, cH		de Gasperin et al. (2014)
SPT-CL J0245-5302	0.300	cdRS		Zheng et al. (2018)
Abell 2552	0.305	cH		Kale et al. (2015)
Abell 2744	0.308	H, mRS		Giovannini et al. (1999); Govoni et al. (2001b)
				Bacchi et al. (2003); Orrú et al. (2007)
				Venturi et al. (2013); Pearce et al. (2017)
				George et al. (2017); Duchesne et al. (2017)
Abell 1300	0.308	H, RS, cRS	USS	Venturi et al. (2013); Reid et al. (1999)
				Giovannini et al. (1999); Parekh et al. (2017)
RXC J2003.5–2323	0.317	H		Venturi et al. (2007); Giacintucci et al. (2009a)
				Venturi et al. (2013)
Abell 1995	0.318	H		Owen et al. (1999); Giovannini et al. (2009)
MACS J0257.6–2209	0.322	cH		Giacintucci et al. (2017)
Abell 1351	0.322	H		Giacintucci et al. (2009b); Giovannini et al. (2009)
WHL J091541.3+251206	0.324	cF		van Weeren et al. (2009b, 2011d)
PSZ1 G094.00+27.41n	0.332	H		Bonafede et al. (2014b); Kale & Parekh (2016)
PSZ1 G108.18–11.53	0.335	dRS, H		de Gasperin et al. (2015a)
MACS J0520.7–1328	0.336	cRS	close to 1WGA J0521.0–1333	Macario et al. (2014)
1WGA J0521.0–1333	0.340	cRS, cH	diffuse emission present	Macario et al. (2014)
			close to MACS J0520.7-1328
RXC J1115.8+0129	0.345	cMH		Pandey-Pommier et al. (2016)
RBS 797	0.345	MH		Gitti et al. (2006); Doria et al. (2012)
MACS J1931.8–2634	0.352	U		Giacintucci et al. (2014b)
MACS J0308.9+2645	0.356	cH		Parekh et al. (2017)
Abell S1121	0.358	H		Duchesne et al. (2017)
RXC J0256.5+0006o	0.360	H		Knowles et al. (2016, 2018)
RX J1532.9+3021	0.363	MH		Kale et al. (2013); Hlavacek-Larrondo et al. (2013)
				Giacintucci et al. (2014b)
MACS J1752.0+4440	0.366	dRS, H		van Weeren et al. (2012a); Bonafede et al. (2012)
ZwCl 1447.2+2619p	0.370	H, RS		Giovannini et al. (2009); Govoni et al. (2012)
ZwCl 0847.2+3617q	0.373	cU		Rudnick & Lemmerman (2009); Kale et al. (2015)
RXC J0949.8+1707r	0.383	cH		Venturi et al. (2008, 2013)
MACS J0949.8+1708	0.383	H		Bonafede et al. (2015a)
PLCK G287.0+32.9	0.390	H, RS, F		Bagchi et al. (2011); Bonafede et al. (2014a)
				George et al. (2017)
PLCKESZ G284.99–23.70	0.390	H		Martinez Aviles et al. (2016, 2018)
RX J1720.2+3536s	0.391	cMH		Giacintucci et al. (2017)202020Venturi et al. in prep
GMBCG J357.91841–08.97978t	0.394	cH		Duchesne et al. (2017)
MACS J0416.1–2403	0.396	H	cUSS	Ogrean et al. (2015); Pandey-Pommier et al. (2015)
MACS J0159.8–0849	0.405	MH		Giacintucci et al. (2014b, 2017)
MACSJ0553.4–3342	0.407	H		Bonafede et al. (2012)
Abell 851	0.407	cH		Owen et al. (1999); Giovannini et al. (2009)
PSZ1 G262.72–40.92	0.421	U		Martinez Aviles et al. (2018)
MACS J0417.5–1154	0.443	H, RS		Dwarakanath et al. (2011); Parekh et al. (2017)
				Sandhu et al. (2018)
MACS J2243.3–0935	0.447	H, dRS		Parekh et al. (2017); Cantwell et al. (2016)
				Duchesne et al. (2017)
MACS J0329.6–0211	0.450	MH		Giacintucci et al. (2014b, 2017)
RX J1347.5–1145	0.452	MH, cRS		Gitti et al. (2007); Ferrari et al. (2011)
PLCK G004.5-19.5	0.516	H, RS		Albert et al. (2017)
ACT-CL J0014.9–0057	0.533	cRS		Knowles et al. (2018)
CL 0016+16u	0.541	H		Moffet & Birkinshaw (1989)
				Giovannini & Feretti (2000)
MACS J1149.5+2223	0.544	RS, cRS, H	USS	Bonafede et al. (2012)
ACT-CL J0045.2–0152	0.545	U		Knowles et al. (2018)
MACS J0717.5+3745	0.546	H, mRS		van Weeren et al. (2009b); Bonafede et al. (2009b)
				Pandey-Pommier et al. (2013); Bonafede et al. (2018)
				Bonafede et al. (2018)
MACS J0025.4–1222	0.586	cdRS		Riseley et al. (2017)
Phoenix cluster	0.596	MH		van Weeren et al. (2014)
ACT-CL J0022.2–0036	0.805	cMH		Knowles et al. (2018)
ACT-CL J0102–4915v	0.870	H, dRS		Menanteau et al. (2012); Lindner et al. (2014)
				Botteon et al. (2016b)
RXC J2351.0–1954	…	cH		Duchesne et al. (2017)

Equations

$$n(E)\mathrm{d}E\propto E^{-p}\mathrm{d}E\text{mbox}.$$

$$p=1-2\alpha \text{mbox}.$$

$$t_{\mathrm{age}}\text{mbox}[yr]\approx 3.2\times 10^{10}\frac{B^{1/2}}{B^2+B_{\mathrm{CMB}}^2}{\left[(1+z)\nu \right]}^{-1/2}\text{mbox},$$

$$\varphi (\mathbf{r})=0.81{\int }_{\mathrm{source}}^{\mathrm{telescope}}{n}_{\mathrm{e}}\mathbf{B}\cdot d\mathbf{r}\text{mbox}[\text{mbox}{radm}^{-2}]\text{mbox},$$

$$\mathrm{RM}=\frac{d\chi ({\lambda }^2)}{d{\lambda }^2}\text{mbox},$$

$$\chi ={\chi }_0+\varphi {\lambda }^2\text{mbox}.$$

$$R_{\mathrm{obs}}\equiv \frac{f_{\mathrm{IC}}(kT)}{f_{\mathrm{sync}}(\nu )}$$

$$\times {\left(\frac{kT}{20\mathrm{keV}}\right)}^{-\Gamma }{\left(\frac{\nu }{\mathrm{GHz}}\right)}^{\Gamma -1}$$

$$\times {\left(\frac{T_{\mathrm{CMB}}}{2.8\mathrm{K}}\right)}^{\Gamma +2}{\left(\frac{B}{\mu \mathrm{Gauss}}\right)}^{-\Gamma }c(p),$$

$$B=\left(\frac{20\mathrm{keV}}{kT}\right){\left(\frac{\nu }{\mathrm{GHz}}\right)}^{(p-1)/(p+1)}{e}^{\frac{2.84(p-r)}{p+1}}\mu \mathrm{Gauss},$$

$$r=0.7\ln \left[\frac{R_{\mathrm{obs}}(kT,\nu )}{1.11\times 10^{-8}}\right].$$

$${\nu }_{\mathrm{b}}\propto M^{4/3}\text{mbox},$$

$${\nu }_{\mathrm{b}}\propto {\left(1+\Delta M/M\right)}^3\text{mbox},$$

$$I(r)=I_0{e}^{-r/r_e}\text{mbox},$$

$$\frac{T_2}{T_1}=\frac{5{\mathcal{M}}^4+14{\mathcal{M}}^2-3}{16{\mathcal{M}}^2}$$

$$\frac{{\rho }_2}{{\rho }_1}=C=\frac{4{\mathcal{M}}^2}{{\mathcal{M}}^2+3}\text{mbox},$$

$${\mathcal{M}}_{\mathrm{radio}}=\sqrt{\frac{2{\alpha }_{\mathrm{inj}}-3}{2{\alpha }_{\mathrm{inj}}+1}}\text{mbox}.$$

Full text

1 Introduction

Galaxy clusters are the largest virialized objects in our Universe, with masses up to $\sim 10^{15}$ M*⊙*. Elongated filaments of galaxies, located between clusters, form even larger unbound structures, making up the cosmic web. Galaxy clusters are located at the nodes of filaments, like “spiders” in the cosmic web.

Clusters contain up to several thousands of galaxies. However, the galaxies comprise only a few percent of a cluster’s total mass. Most of the baryonic mass of clusters is contained in a hot ($10^{7}$–$10^{8}$ K) ionized intracluster medium (ICM), held together by the clusters’ gravitational pull. This dilute magnetized plasma ($\sim 10^{-3}$ particles cm*-3*) emits thermal Bremsstrahlung at X-ray wavelengths, permeating the cluster’s volume (e.g., Mitchell et al., 1976; Serlemitsos et al., 1977; Forman & Jones, 1982), see Figure 1. The ICM makes up $\sim$15% of a cluster’s mass budget. Most of the mass, $\sim 80\%$, is in the form of dark matter (e.g., Blumenthal et al., 1984; White & Fabian, 1995; Jones & Forman, 1999; Arnaud & Evrard, 1999; Sanderson et al., 2003; Vikhlinin et al., 2006).

Elongated filaments of galaxies span the regions between clusters. The so-called warm-hot intergalactic medium (WHIM) pervades these galaxy filaments (Cen & Ostriker, 1999). Compared to the ICM, the intergalactic medium of galaxy filaments (WHIM) has a significantly lower density ($\lesssim 10^{-4}$ particles cm*-3*) and cooler temperature ($10^{5}$–$10^{7}$ K). About half of the Universe’s baryons reside in this WHIM (e.g., Cen & Ostriker, 1999; Davé et al., 2001; Eckert et al., 2015). Galaxy filaments are expected to be surrounded by strong accretion shocks, where the plasma is first shock-heated (Sunyaev & Zeldovich, 1972). However, studying the WHIM and associated shocks is difficult due to a lack of sensitive observational tools. Galaxy clusters form by accretion from the WHIM and through a sequence of mergers of clusters and groups (e.g., Peebles & Yu, 1970; Press & Schechter, 1974; Voit, 2005; Kravtsov & Borgani, 2012). Cluster mergers are very energetic events, releasing energies up to $\sim 10^{64}$ ergs on a few Gyr timescale. This energy is dissipated through low-Mach number shocks and turbulence, heating the ICM (e.g., Markevitch & Vikhlinin, 2007). Clusters can thus be divided as either “relaxed” (undisturbed) or “merging” (disturbed) systems, depending on their dynamical (merging) state.

Galaxy clusters often host a number of active galactic nuclei (AGN) that emit radio synchrotron emission (i.e., radio galaxies) (e.g., De Young, 1984; de Young, 2002; Tadhunter, 2016). The sizes of these sources range from a few kpc to about $\sim$1 Mpc, extending well beyond the host galaxy. A major difference with radio galaxies that are located outside clusters (and groups) is that the jets and lobes of cluster radio galaxies often show signs of interaction with the ICM (e.g., Miley, 1980; Burns, 1998; Johnston-Hollitt et al., 2015a). These interactions result in morphologies that range from wide-angle (WAT), narrow angle (NAT), to “head-tail” radio sources.

Gas in the central regions of many relaxed clusters has a radiative cooling time that is much shorter than the Hubble time. In the absence of a heating source, a cooling flow is expected to develop, whereby the temperature in the central region of the cluster drops and gas flows inwards (e.g., Fabian, 1994; Peterson & Fabian, 2006; Fabian, 2012; McNamara & Nulsen, 2012). X-ray observations do show these temperature drops in some cluster cores (“cool core” clusters), but there is much less cool gas than what would be expected from the short radiative cooling time (Kaastra et al., 2001; Peterson et al., 2001, 2003). Therefore, some source of heating must balance the radiative losses. Radio galaxies, associated with the brightest cluster galaxy (BCG), have been identified as the main source of energy input into the ICM. X-ray observations show numerous cavities in cool core clusters, coincident with the lobes of the central radio galaxy. Here the radio plasma has displaced the X-ray emitting gas, creating a low-density bubble which rises buoyantly and expands, distributing energy to the surrounding ICM (e.g., Churazov et al., 2002). This process is commonly referred to as “radio-mode” feedback, although it is still being debated what the precise mechanism is that transfers the energy to the ICM.

1.1 Extended synchrotron radio emission from galaxy clusters

Radio observations have shown that the ICM can also contain a non-thermal component of cosmic rays (CR, see Figure 1) which is not directly associated with cluster radio galaxies (e.g., Large et al., 1959; Willson, 1970). These GeV CR electrons (i.e., Lorentz factors of $\gamma>10^{3}$) emit synchrotron radiation in the presence of $\sim\mu$Gauss ICM magnetic fields. During the last decade significant progress has been made in our understanding of this non-thermal component, through observations, theoretical, and numerical work. There is now compelling evidence that ICM shocks waves, and likely also turbulence, are able to (re-)accelerate particle to relativistic energies creating this non-thermal CR component of the ICM.

The presence of extended synchrotron emission also indicates the existence of large-scale ICM magnetic fields with a strength of the order of 0.1–10 $\mu$Gauss (e.g., Brüggen et al., 2012; Clarke et al., 2001; Johnston-Hollitt et al., 2015b). Cluster magnetic fields play an important role in particle acceleration processes. Additionally, magnetic fields inhibit transport processes like heat conduction, spatial mixing of gas, and the propagation of cosmic rays (e.g., Pfrommer et al., 2017; Ruszkowski & Oh, 2010). However, few details are known about the precise properties of these fields since they are difficult to measure (e.g., Govoni & Feretti, 2004).

The synchrotron emitting CR electrons should scatter photons from the cosmic microwave background (CMB) to X-ray energies, resulting in a hard tail on top of the thermal X-ray spectrum of clusters (Rephaeli, 1979; Rephaeli et al., 1994; Sarazin & Kempner, 2000). So far, no conclusive detection of this inverse-Compton (IC) radiation has been made (e.g., Fusco-Femiano et al., 2000, 2001; Rephaeli & Gruber, 2004; Rossetti & Molendi, 2004; Fusco-Femiano, 2004; Rephaeli et al., 2008; Eckert et al., 2008; Wik et al., 2009, 2014). However, even a non-detection of IC X-ray emission, in combination with radio observations, is useful to set lower limits on the ICM magnetic fields strength (e.g., Sugawara et al., 2009; Finoguenov et al., 2010; Itahana et al., 2015). Similarly, CR protons can interact hadronically with the protons of the ICM and generate pions that can then decay into gamma-rays (c.f., Dennison, 1980; Blasi & Colafrancesco, 1999; Blasi et al., 2007). Gamma-ray observations are particularly important to understand the dynamical role of CR protons in clusters, and the role of secondary electrons, also coming from pion decays, in generating the extended radio emission.

1.2 This review

Galaxy clusters provide a unique environment to study the physics of particle acceleration in collisionless, high-$\beta$, turbulent plasmas, where $\beta$ is the ratio of the thermal pressure to the magnetic pressure111$\beta=\frac{8\pi nT}{B^{2}}\sim 100$ for the ICM, taking $T=5$ keV, $B=3$ $\mu$Gauss, and $n=5\times 10^{-3}$ cm*-3*, and at low Mach numbers shocks. Furthermore, diffuse radio emission from clusters can be used as a signpost of ICM shocks and turbulence, which are often difficult to detect and characterize at other wavelengths. Since shocks and turbulence trace the dynamical state of the ICM, radio observations also provide us with a probe of the cluster’s evolutionary stage, important for our understanding of structure formation in the Universe. Finally, diffuse radio emission can be used as a complementary method to discover clusters that were missed by X-ray, SZ, or optical surveys (Brown et al., 2011a; van Weeren et al., 2012b; Macario et al., 2014; de Gasperin et al., 2017b).

In this paper we review the observational properties of diffuse extended cluster radio emission. Previous observational reviews on this subject were presented by Feretti (2002); Giovannini & Feretti (2002); Feretti (2003); Ferrari et al. (2008); Feretti et al. (2012). Here we provide an update, encompassing recent results that have helped to improve our understanding of these sources. For a more theoretical review we refer the reader to Brunetti & Jones (2014). Observational progress in this field has been made through a combination of high-resolution multi-frequency studies, the availability of deep low-frequency observations, an increasing number of polarimetric studies, the compilation of larger cluster samples with deep radio data, and high-frequency detections. The joint analyses of radio data and observations at other wavelengths, in particular in the X-ray and Gamma-ray bands, has also played an important role.

The outline of this paper is as follows. In Section 2 we briefly discuss synchrotron radiation and particle acceleration mechanisms. The classification of diffuse cluster radio sources is discussed in Section 3. A review of cluster magnetic fields is given in Section 4. Overviews of radio halos, including mini-halos, and cluster radio shocks and revived fossil plasma sources are presented in Sections 5 and 6. In Section 7 we end this review with a discussion on the detection of diffuse radio emission outside cluster environments.

2 Synchrotron radiation and radio spectra

In this section we briefly discuss some relevant theory about the synchrotron spectra of CR electrons. For a more detailed treatment of synchrotron radiation we refer the reader to the references provided in Feretti et al. (2012). A standard assumption is that the ICM CR population can be described by a power law energy ($E$) distribution

[TABLE]

The index of the energy (or momentum) distribution $p$ is directly related to the radio spectral index222$F_{\nu}\propto\nu^{\alpha}$, where $\alpha$ is the spectral index

[TABLE]

Diffuse cluster radio emission typically has a steep spectral index , i.e., $\alpha\lesssim-1$. The spectral shape is related to the physics of the acceleration mechanism and the electron synchrotron and IC energy losses. The characteristic lifetime ($t_{\mathrm{age}}$) of the synchrotron emitting electrons ($\gamma\sim 10^{4}$; GeV energy) due to these energy losses is

[TABLE]

where $B$ the magnetic field strength, $z$ the source redshift, $B_{\mathrm{CMB}}$ the equivalent magnetic field strength of the CMB ($B_{\mathrm{CMB}}\mbox{ }[\mu\rm{Gauss}]\approx 3.25\left(1+z\right)^{2}$), and $\nu$ the observing frequency in MHz. In clusters, we have $t_{\mathrm{age}}\lesssim 10^{8}$ yrs. The typical diffusion length-scale in the ICM of a GeV electron, using the Bohm approximation, is of the order of 10 pc (e.g., Bagchi et al., 2002). Plasma motions can increase the distance over which GeV electrons travel, but this distance is still expected to remain well below a Mpc. This means that Mpc-scale diffuse radio sources cannot trace CR electrons that are accelerated at a single location in the ICM. Instead, they need to be (re-)accelerated or produced in-situ (Jaffe, 1977), providing important constraints on the possible acceleration/production mechanisms.

Due to the energy losses, the initial power-law spectrum steepens beyond a break frequency, whose position is related to the time since acceleration. The power-law spectrum is commonly refereed to as the injection spectrum, characterized by an injection spectral index ($\alpha_{\rm{inj}}$). For the JP (Jaffe-Perola) synchrotron spectrum (Jaffe & Perola, 1973), one assumes that there is a continuous isotropization of the electron pitch angles (i.e., angle between the magnetic field and the electron velocity) on a timescale that is shorter than $t_{\mathrm{age}}$. A JP spectrum describes a synchrotron spectrum from a single burst of acceleration and then aging. The KP (Kardashev-Pacholczyk) model (Kardashev, 1962; Pacholczyk, 1970) also represents such a spectrum, but without the isotropization of the pitches angles. A collection of spectral shapes is displayed in Figure 2.

Since it is usually difficult to spatially isolate electrons that all have the same spectral age, there are also composite models. These models sum JP (or KP) spectra with different amounts of spectral aging. The CI (continuous injection) composite model (Pacholczyk, 1970) describes the integrated spectrum of a source with continuous particle injection. For the KGJP/KGKP (Komissarov-Gubanov) model (Komissarov & Gubanov, 1994), the particles are only injected for a finite amount of time before the injection in the source stops.

2.1 Particle acceleration mechanisms

There are several physical mechanisms to accelerate particles in the ICM and produce the synchrotron emitting CR electrons. We briefly give an overview of these processes below. Further details will be discussed in later sections where relevant.

• •

First order Fermi acceleration (Fermi-I): This process, also known as diffusive shock acceleration (DSA), plays an important role in various astrophysical environments (Krymskii, 1977; Axford et al., 1977; Bell, 1978a, b; Blandford & Ostriker, 1978; Drury, 1983; Blandford & Eichler, 1987; Jones & Ellison, 1991; Malkov & O’C Drury, 2001). For DSA, particles are accelerated at a shock with the acceleration taking place diffusively. In this process, particles cross back and forward across the shock front as they scatter from magnetic inhomogeneities in the shock down and upstream region. At each crossing, particles gain additional energy, forming a power-law energy distribution of CR.

• •

Second order Fermi acceleration (Fermi-II): This is a stochastic process where particles scatter from magnetic inhomogeneities, for example from magneto-hydrodynamical (MHD) turbulence (Schlickeiser et al., 1987; Schlickeiser & Achatz, 1993; Brunetti et al., 2001; Petrosian, 2001). Particles can either gain or loose energy when scattering. When the motions are random, the probability for a head-on collision, where energy is gained, is slightly larger. Because of its random nature, second order Fermi acceleration is an inefficient process.

• •

Adiabatic compression: A shock wave can adiabatically compress a bubble/lobe/cocoon of (old) relativistic radio plasma from an AGN. Due to the compression, the CR electrons in the cocoon regain energy boosting the radio synchrotron emission (Enßlin & Gopal-Krishna, 2001; Enßlin & Brüggen, 2002).

• •

Secondary models: Another mechanism to produce CR electrons is via a secondary process, meaning that the CR electrons are produced as secondary particles (decay products). In the hadronic model, collisions between relativistic protons and the thermal ions produce secondary CR electrons (Dennison, 1980; Blasi & Colafrancesco, 1999; Dolag & Enßlin, 2000; Miniati et al., 2001a; Keshet & Loeb, 2010; Donnert et al., 2010; Enßlin et al., 2011). Since CR protons have a very long lifetime compared to CR electrons, they will accumulate over the lifetime of a cluster once they are accelerated. Possible mechanisms to produce CR protons are first order Fermi acceleration at shocks, AGN activity, and galactic outflows (supernovae, winds).

3 Classification

Diffuse cluster radio sources have historically been divided into three main classes, relics, halos, and mini-halos (Feretti & Giovannini, 1996). In addition, radio filaments were proposed to trace the large-scale filaments of the cosmic web, outside of clusters. Note that the term filament has also sometimes been used to describe radio relics (or relic-type structures) in clusters. We will discuss radio emission outside the cluster environment in Section 7.

Radio halos are centrally located diffuse sources in merging clusters. They do not have any optical counterparts. Mini-halos have smaller sizes and are located in relaxed cool core clusters which also host a powerful radio galaxy associated with the BCG. Radio relics have been defined as extended sources that show high levels of polarization ($\gtrsim 10\%$ at GHz frequencies) and are located in the cluster periphery. Similar to radio halos, they not show optical counterparts. Relics were further subdivided (Kempner et al., 2004) into large Radio Gischt, large Mpc-size sources that trace particles accelerated at shocks via Fermi-I processes; Radio Phoenices, AGN fossil plasma compressed and revived by merger shocks; and AGN Relics, fossil radio plasma that is passively evolving from an AGN that has been switched off. For radio relics, the boundaries between the different categories is not always very obvious and the term relics itself is somewhat unfortunate because large relics could be “young” sources with on-going (re-)acceleration.

Here we propose to classify cluster emission into three broad classes:

• $\bullet$

Radio halos are extended sources that roughly follow the ICM baryonic mass distribution. This class includes giant radio halos and mini-halos, see Figure 3. This class would also contain possible “intermediate” or “hybrid” radio halos, with properties falling somewhere in between those of classical giant radio halos and mini-halos. Another property of the halo class is that these sources are not localized, in the sense that particle (re-)acceleration/production occurs throughout a significant volume of the cluster and is not associated with a particular shock which location can be pint-pointed. In terms of a physical interpretation, these “global” sources should trace Fermi-II processes and/or secondary electrons.

• $\bullet$

Cluster radio shocks (radio relics) are extended diffuse sources tracing particles that are (re-)accelerated at ICM shock waves (Figure 3). They have commonly been referred to as radio relics. This radio shock classification is somewhat similar to the that of Gischt, but it does not necessarily require DSA or Fermi-I type acceleration. In that sense, cluster radio shocks are an observationally defined class, unrelated to the details of the actual acceleration mechanism. However, based on our current understanding of these sources, we do anticipate that in most cases cluster radio shocks are associated with Fermi-I acceleration processes. It is not required that cluster radio shocks are located in the cluster periphery, although for large cluster radio shocks that will typically be the case. Due to their nature, the large majority of these sources are expected to show a high degree of polarization. Sources previously classified as large radio relics, Gischt, and double relics, fall in the cluster radio shock category. Unlike radio halos, cluster radio shocks can be associated to a specific cluster region where a shock wave is present, or where a shock wave recently passed. A drawback of the radio shock classification is that the detection of shocks in the ICM is observationally challenging. Therefore, the classification will remain uncertain for some sources. However, for a number of sources the presence of a shock at their location has been confirmed by X-ray observations (see Section 6.1.5) which we argue warrants the creation of a radio shock class. In this review we will use the term radio shock for sources previously classified as large radio relics, Gischt, and double relics. It is important to keep in mind that for a number of sources the presence of a shock remains to be confirmed.

• $\bullet$

Revived AGN fossil plasma sources, phoenices, and GReET In this class we group sources that trace AGN radio plasma that has somehow been re-energized through processes in the ICM, unrelated to the radio galaxy itself. Low-frequency observations are starting to reveal more and more of these type of sources. However, their precise origin and connection to cluster radio shocks and possibly also halos is still uncertain. The main observational property that the sources have in common is the AGN origin of the plasma and their ultra-steep radio spectra due to their losses. For this review we decided to keep the radio phoenix classification (Kempner et al., 2004). Often these phoenices display irregular filamentary morphologies. They have relatively small sizes of at most several hundreds of kpc. Gently re-energized tails (GReETs; de Gasperin et al., 2017a) are tails of radio galaxies that are somehow revived, showing unexpected spectral flattening, opposite from the general steepening trend caused by electron energy losses. With the new and upgraded low-frequency radio telescopes that have become operational, we expect that the nature of these revived fossil plasma sources will become more clear over the next decade.

Fossil radio plasma plays and important role in some of the models for the origin of radio halos and cluster radio shocks. In these models fossil plasma is re-accelerated via first and second order Fermi processes. This implies that when clusters are observed at low enough frequencies, both halos and cluster radio shocks will blend with regions of old AGN radio plasma, complicating the classification.

The classification can also be hindered by projection effects. For example, a cluster radio shock observed in front of the cluster center might mimic halo-type emission if the signal to noise of the image is not very high. However, these are observation related difficulties, which can in principle be resolved with better data. On the website http://galaxyclusters.com we provide an up to date list of the currently known diffuse cluster radio sources and their classification. An up-to-date list of clusters with (candidate) diffuse radio emission at the time of writing (September 2018) is shown in Table LABEL:tab:clusterlist.

4 Cluster magnetic fields

4.1 Global

Magnetic fields permeate galaxy clusters and the intergalactic medium on Mpc-scales. These fields play key roles in particle acceleration and on the process of large scale structure formation, having effects on turbulence, cloud collapse, large-scale motions, heat and momentum transport, convection, viscous dissipation, etc. In particular, cluster magnetic fields inhibit transport processes like heat conduction, spatial mixing of gas, and propagation of cosmic rays. The origin of the fields that are currently observed remains largely uncertain. A commonly accepted hypothesis is that they result from the amplification of much weaker pre-existing seed fields via shock/compression and/or turbulence/dynamo amplification during merger events and structure formation, and different magnetic field scales survive as the result of turbulent motions (e.g., Kahniashvili et al., 2013). The origin of seed fields is is unknown. They could be either primordial, i.e., generated in the early Universe prior to recombination, or produced locally at later epochs of the Universe, in early stars and/or (proto)galaxies, and then injected in the interstellar and intergalactic medium (Rees, 2006). For a review about magnetic field amplification in clusters we refer the reader to Donnert et al. (2018).

Magnetic fields are difficult to measure. Some estimates have relied on the idea that the energies in cosmic rays and magnetic fields in the radio emitting regions are the same (“equipartition”; Beck & Krause, 2005). In this way, magnetic field values in the range 0.1–10 $\mu$Gauss are obtained. However, this method is inherently uncertain due to the many assumptions that are required. Cosmological simulations of clusters predict $\mu$Gauss-level magnetic field strengths in the cluster centers and a decrease of the magnetic field strength with radius in the outer regions (Dolag et al., 1999, 2001, 2002; Vazza et al., 2014, 2018). These values are roughly consistent with equipartition magnetic field strengths estimates of the order of a $\mu$Gauss.

The most promising technique to derive a more detailed view of the magnetic fields in clusters is via the analysis of the Faraday rotation of radio galaxies located inside and behind the cluster (e.g., Clarke, 2004; Govoni & Feretti, 2004). Faraday rotation changes the intrinsic polarization angle ($\chi_{0}$). The Faraday depth ($\phi$) is related to the properties of the plasma that cause the Faraday rotation (Burn, 1966; Brentjens & de Bruyn, 2005)

[TABLE]

where $n_{\mathrm{e}}$ is the electron density in units of cm*-3*, $\mathbf{B}$ the magnetic field in units of $\mu$Gauss, and $d\mathbf{r}$ is an infinitesimal path length in along the line of sight in units of parsec. The rotation measure (RM) is defined as

[TABLE]

where $\lambda$ is the observing wavelength. The Faraday depth equals the RM if there is only one source along the line of sight (and there is no internal Faraday rotation). This means that the RM does not depend on the observing wavelength. Also, all polarized emission comes from a single Faraday depth $\phi$ and the measured polarization angle ($\chi$) is given by

[TABLE]

From RM measurements, the strength and structure of cluster magnetic fields can be constrained by semi-analytical approaches, numerical techniques or RM synthesis (Brentjens & de Bruyn, 2005). To this aim, a spherically symmetric model ($\beta$-model) is generally assumed for the thermal gas. Moreover, one needs to assume that the interaction between the ICM and the radio galaxy plasma does not affect the measured RM. It is still being debated to what extent this assumption holds. Deviations of the Faraday rotation from the simple $\lambda^{2}$–law (Equation 6) have been detected (e.g., Bonafede et al., 2009b), likely implying either that the magnetized screen is non–uniform and/or that the ICM thermal plasma is mixed with the relativistic plasma.

4.1.1 Results from RM studies

The presence of magnetic field in clusters is demonstrated by statistical studies. The comparison between the RMs of polarized extragalactic radio sources in the line of sight of galaxy clusters and RM measurements made outside of the projected cluster regions shows excess of the standard deviations of RM values in the cluster areas (c.f., Clarke et al., 2001; Böhringer et al., 2016), see Figure 4. This is consistent with ubiquitous cluster magnetic fields of a few $\mu$Gauss strength, coherent cells of about 10 kpc, and a magnetic field energy density of a few per mille of the thermal energy density.

Information about the magnetic field in individual clusters through RM studies has been obtained so far for about 30 objects, including both merging and relaxed clusters. The best studied cluster is Coma, whose magnetic field has been obtained with RM information on 7 radio galaxies in the cluster central region (Bonafede et al., 2010), and 7 additional radio galaxies in the peripheral Coma southwest region, where the NGC 4839 infalling group and the cluster radio shock are located (Bonafede et al., 2013). A single-cell model is not appropriate to describe the observed data, which are generally consistent with a turbulent field following a Kolmogorov power-law spectrum. From energy considerations, i.e., to avoid that the magnetic pressure exceeds the thermal pressure in the outer cluster regions, it is inferred that the magnetic field profile scales with the gas density $n_{th}$ as $B\propto n_{th}^{\eta}$. The value of the index $\eta$ reflects the magnetic field formation and amplification. It is expected that $\eta$=2/3 in the case of adiabatic compression during a spherical collapse due to gravity. In this case, the field lines are frozen into the plasma and compression of the plasma results in compression of the flux lines (as a consequence of magnetic flux conservation). A value $\eta$=1/2 is instead expected if the energy in the magnetic field scales as the energy in the thermal plasma. Other values of $\eta$ may be obtained by specific combinations of compression orientation and magnetic field orientation.

The Coma cluster magnetic field is well represented by a Kolmogorov power spectrum with minimum scale of $\sim$2 kpc and maximum scale of $\sim$34 kpc. The central field strength is 4.7 $\mu$Gauss and the radial slope is $\propto n_{th}^{0.7}$ (Bonafede et al., 2010), see Figure 5. The magnetic field of the southwest peripheral region is found to be $\sim$2 $\mu$Gauss, i.e., higher than that derived from the extrapolation of the radial profile obtained for the cluster center; a boost of magnetic field of $\sim$ a factor of 3 is required. The magnetic field amplification does not appear to be limited to the cluster radio shock region, but it must occur throughout the whole southwestern cluster sector, including the NGC 4839 group (Bonafede et al., 2013).

In the clusters analyzed so far, it is derived that cool core clusters have central magnetic field intensities of the order of a few 10 $\mu$Gauss, while merging clusters are characterized by intensities of a few $\mu$Gauss. The fields are turbulent, with spatial scales in the range 5–500 kpc, and coherence lengths of a few 10 kpc. The values of the profile index $\eta$ are in the range 0.4–1, therefore no firm conclusion can be drawn on the radial trend of the magnetic field. Recently, Govoni et al. (2017) found a correlation between the central electron density and mean central magnetic field strength ($\eta$=0.47) using data for 9 clusters. No correlation seems to be present between the mean central magnetic field and the cluster temperature. In conclusion, good information about the central magnetic field intensity in clusters has been obtained, whereas the magnetic field structure (profile, coherence scale, minimum and maximum scales, power spectrum, link to cluster properties) is still poorly known.

4.1.2 Statistical studies from fractional polarization

From the analysis of the fractional polarization of radio sources in a sample of X-ray luminous clusters from the NVSS, a clear trend of the fractional polarization increasing with the distance from the cluster center has been derived (Bonafede et al., 2011). The low fractional polarization in sources closer to the cluster center is interpreted as the result of higher beam depolarization, occurring in the ICM because of fluctuations within the observing beam and higher magnetic field and gas densities in these regions. Results are consistent with fields of a few $\mu$Gauss, regardless of the presence or not of radio halos. A marginally significant difference between relaxed and merging clusters has been found.

4.1.3 Lower limits from IC emission

CR electrons present in the ICM should scatter photons from the CMB, creating a hard power-law of X-ray emission, on top of the thermal Bremsstrahlung from the ICM (Rephaeli, 1979; Rephaeli et al., 1994; Sarazin & Kempner, 2000). Despite several claims made over the last decades, it seems that there is no conclusive evidence yet for this IC emission from the diffuse CR component of the ICM (e.g., Fusco-Femiano et al., 2000, 2001; Rephaeli & Gruber, 2004; Rossetti & Molendi, 2004; Fusco-Femiano, 2004; Rephaeli et al., 2008; Eckert et al., 2008; Wik et al., 2009, 2014; Ajello et al., 2009; Molendi & Gastaldello, 2009; Kawaharada et al., 2010; Wik et al., 2012; Gastaldello et al., 2015). The difficultly associated with the detection of IC emission is related to the requirement of accurately modeling the contributions of the instrumental and astronomical backgrounds.

Following Petrosian (2001); Randall et al. (2016), the monochromatic IC X-ray and synchrotron radio flux ratio ($R_{\rm obs}$) can be written as

[TABLE]

where $\Gamma=(p+1)/2$, $p$ is the power-law slope of the electron energy distribution $N(E)\propto E^{-p}$ (see Equation 2 for the relation between radio spectral index $\alpha$ and $p$), $f_{\rm IC}(kT)$ is the IC flux density at energy $kT$, $f_{\rm sync}(\nu)$ is the synchrotron flux density at frequency $\nu$, $T_{\rm CMB}$ is the CMB temperature at the cluster’s redshift, and $c(p)$ is a normalization factor that is a function of $p$. For typical values of $p$, $10<c(p)<1000$, see Rybicki & Lightman (1979). The function $c(p)$, for values of $2\lesssim p\lesssim 5$ can be approximated as $c(p)\approx e^{1.42p-0.51}$. With Equation 7 and this approximation the expression for the magnetic field strength becomes

[TABLE]

In the above derivations a power-law distribution of electrons down to low energies is assumed. If this assumption does not hold (e.g., Bartels et al., 2015), for example because there is flattening of the spectrum at low frequencies, the magnetic field values will be overestimated.

By deriving upper limits on the IC X-ray emission and combining that with radio flux density measurements of radio halos, lower limits on the global ICM magnetic field strength can be computed. For radio halos, it is generally challenging to obtain stringent lower limits. The reason is that radio halos are typically faint. In addition, the IC emission is co-spatial with the thermal ICM, making it harder to separate the components. Furthermore, bright radio galaxies located in the cluster center can also produce non-thermal X-ray emission. The obtained lower magnetic field strength limits are therefore less constraining than the ones obtained for radio shocks (see Section 4.2). The lower limits that have been computed for radio halo hosting clusters range around $0.1-0.5$ $\mu$Gauss. For example, for the Coma cluster Rossetti & Molendi (2004) found $B>0.2-0.4$ $\mu$Gauss and Wik et al. (2009) reported $B>0.15$ $\mu$Gauss. For the Bullet cluster a limit of $B>0.2$ $\mu$Gauss was determined (Wik et al., 2014). Magnetic field strength limits for the cluster Abell 2163 are $B>0.2$ $\mu$Gauss and $B>0.1$ $\mu$Gauss (Sugawara et al., 2009; Ota et al., 2014). A recent overview of constraints on the volume-average magnetic field for radio halo and relic hosting clusters is given by Bartels et al. (2015).

4.2 Magnetic fields at cluster radio shocks

Similar to radio halos, measurements of IC X-ray emission can be used to determine magnetic field strength at the location of cluster radio shocks (Rephaeli, 1979; Rephaeli et al., 1994; Sarazin & Kempner, 2000; Randall et al., 2016), but so far no undisputed detections have been made. With deep X-ray observations, mostly from the XMM-Newton and Suzaku satellites, interesting lower limits on the magnetic field strength have been determined. Finoguenov et al. (2010) placed a lower limit of $3$ $\mu$Gauss on the northwest cluster radio shock region in Abell 3667, consistent with an earlier reported lower limit of $1.6$ $\mu$Gauss by Nakazawa et al. (2009). Itahana et al. (2015) reported a lower limit of $1.6$ $\mu$Gauss for the Toothbrush Cluster. For the radio shock in the cluster RXC J1053.7+5453, the lower limits was found to be $0.7$ $\mu$Gauss (Itahana et al., 2017).

Another method to constrain the magnetic field strength at the location of cluster radio shocks is to use the source’s width. Here the assumption is that the source’s width is determined the characteristic timescale of electron energy losses (synchrotron and IC) and the shock downstream velocity. Using this method, values of either $\sim$1 or $\sim$5 $\mu$Gauss were found for the Sausage Cluster (van Weeren et al., 2010). However, recent work by Rajpurohit et al. (2018) suggests that there are more factors affecting the downstream radio brightness profiles making the interpretation more complicated, for example, due to the presence of filamentary structures in the radio shock and a distribution of magnetic fields strengths (see also Di Gennaro et al., 2018). Taking some of these complications into account, Rajpurohit et al. (2018) concluded that the magnetic field strength is less than 5 $\mu$Gauss for the Toothbrush cluster.

4.3 Future prospects

Surveys at frequencies of $\gtrsim 1$ GHz, such ongoing VLA Sky Survey at 2–4 GHz (VLASS; Lacy et al., 2016; Myers et al., 2016), and future surveys carried out with MeerKat (Booth et al., 2009; Jonas, 2009), ASKAP (Norris et al., 2011; Gaensler et al., 2010), and WSRT-APERTIF (Verheijen et al., 2008; Adams et al., 2018) will provide larger samples of polarized radio sources that can be utilized for ICM magnetic field studies. In the more distant future, the SKA will provide even larger samples. This will enable the detailed characterization of magnetic fields in some individual (nearby) clusters, employing background and cluster sources (Krause et al., 2009; Bonafede et al., 2015b; Johnston-Hollitt et al., 2015b; Roy et al., 2016).

Another important avenue to further pursue are hard X-ray observations to directly measure the IC emission from the CRe in the ICM (e.g., Bartels et al., 2015). This will enable direct measurements of the ICM magnetic field strength at the location of radio shocks and halos.

5 Radio halos

5.1 Giant radio halos

Radio halos are diffuse extended sources that roughly follow the brightness distribution of the ICM. Giant Mpc-size radio halos are mostly found in massive dynamically disturbed clusters (Giovannini et al., 1999; Buote, 2001; Cassano et al., 2010b). The prototypical example is the radio halo found in the Coma cluster (e.g., Large et al., 1959; Willson, 1970; Giovannini et al., 1993; Thierbach et al., 2003; Brown & Rudnick, 2011). In Table LABEL:tab:clusterlist we list the currently known giant radio halos and candidates. Some examples of clusters hosting giant radio halos are shown in Figure 6.

Giant radio halos have typical sizes of about 1–2 Mpc. The most distant radio halo is found in El Gordo at $z=0.87$ (Menanteau et al., 2012; Lindner et al., 2014; Botteon et al., 2016b). The 1.4 GHz radio powers of observed halos range between about $10^{23}$ and $10^{26}$ W Hz*-1*, with the most powerful radio halo ($P_{\rm 1.4GHz}=1.6\times 10^{26}$ W Hz*-1*) being present in the quadruple merging cluster MACS J0717.5+3745 (Bonafede et al., 2009b; van Weeren et al., 2009c). The radio halo with the lowest power known to date ($P_{\rm 1.4GHz}=3.1\times 10^{23}$ W Hz*-1*) is found in ZwCl 0634.1+4747 (Cuciti et al., 2018). Other noteworthy examples are the double radio halos in the pre-merging cluster pairs Abell 399–401 (Murgia et al., 2010b) and Abell 1758N–1758S (Botteon et al., 2018a).

Currently there are about 65 confirmed radio halos. Initially, most halos were found via the NVSS333NRAO VLA Sky Survey (Condon et al., 1998) and WENSS444Westerbork Northern Sky Survey (Rengelink et al., 1997) surveys (e.g., Giovannini et al., 1999; Kempner & Sarazin, 2001; Rudnick & Lemmerman, 2009; van Weeren et al., 2011b; George et al., 2017). More recently, halos have been uncovered with targeted GMRT campaigns555Giant Metrewave Radio Telescope (Venturi et al., 2008, 2007; Kale et al., 2013, 2015; Knowles et al., 2018), and via low-frequency surveys such as GLEAM666GaLactic and Extragalactic All-sky MWA Survey (Wayth et al., 2015; Hurley-Walker et al., 2017) and LoTSS777The LOFAR Two-metre Sky Survey (Shimwell et al., 2017, 2018). In addition, radio halo searches have been carried out with the VLA888Very Large Array, ATCA999Australia Telescope Compact Array, MWA101010Murchison Widefield Array, KAT-7111111Seven-dish MeerKAT precursor array, and LOFAR121212LOw-Frequency ARray (Giovannini et al., 2009; Shakouri et al., 2016; Martinez Aviles et al., 2016, 2018; Bernardi et al., 2016; Cuciti et al., 2018; Wilber et al., 2018a; Savini et al., 2018a).

5.1.1 Morphology

Radio halos typically have a smooth and regular morphology with the radio emission approximately following the distribution of the thermal ICM. This is supported by quantitative studies which find a point-to-point correlation between the radio and X-ray brightness distributions (Govoni et al., 2001a; Feretti et al., 2001; Giacintucci et al., 2005; Brown & Rudnick, 2011; Rajpurohit et al., 2018)), although there are some exceptions. One example is the Bullet cluster, where no clear correlation is found (Shimwell et al., 2014).

A few radio halos with more irregular shapes have been uncovered (e.g., Giacintucci et al., 2009b; Giovannini et al., 2009, 2011). One striking example is MACS J0717.5+3745, where a significant amount of small scale structure is present within the radio halo (van Weeren et al., 2017a). Although, it is not yet clear whether these structures really belong to the radio halo or if they are projected on top of it. Two other peculiar cases are the “over-luminous” halos in the low luminosity X-ray cluster Abell 1213 (Giacintucci et al., 2009b) and 0217+70 (Brown et al., 2011a). Giovannini et al. (2011) discussed the interesting possibility that over-luminous halos represent a new class. However, better data is required to further investigate this possibility since none of these “peculiar” halos have been studied in great detail, making the classification and interpretation more uncertain. For example, the peculiar “halo” in A523 has also been classified as a possible radio shock by van Weeren et al. (2011b).

5.1.2 Radio spectra

The spectral properties of radio halos can provide important information about their origin. Therefore, considerable amount of work has gone into measuring the spectral properties of halos.

A complication is that reliable flux density measurements of extended low signal to noise ratio sources are often not trivial to obtain. Reported uncertainties on flux density measurements in the literature often take into account the (1) map noise, assuming the noise is Gaussian distributed and not varying spatially across the radio halo, (2) flux-scale uncertainty, usually somewhere between 2 and 20%, and (3) uncertainty in the subtraction of flux from discrete sources embedded in the diffuse emission. Correctly assessing latter effect can be hard, in particular at low frequencies when extended emission from radio galaxies (i.e., their tails and lobes) becomes more prominent and partly blends with the halo emission. Errors from incomplete uv-coverage and deconvolution are usually not included in the uncertainties. However, in principle they can be determined but this requires some amount of work. The uncertainties related to calibration errors, for example coming from model incompleteness or ionosphere, are often not fully taken into account. Calibration errors affect discrete source subtraction, the map noise distribution, deconvolution, and can lead to flux “absorption”. For the above reasons, the reported uncertainties on radio halo flux-density measurements and spectral index maps in the literature can usually be thought of as lower limits on the true uncertainty.

5.1.3 Integrated spectra

Most radio halos have integrated spectral indices in the range $-1.4<\alpha<-1.1$ (e.g., Giovannini et al., 2009).

The spectral information of most radio halos is based on measurements at just two frequencies. Recently, two systematic campaigns have been carried out with the GMRT to follow-up clusters at lower frequencies to obtain spectra (Macario et al., 2013; Venturi et al., 2013). Flux density measurements at more than three frequencies that also cover a large spectral baseline are rare. Therefore, deviations from power-law spectral shapes are difficult to detect. The best example of a radio halo with an observed spectral steepening, displayed in Figure 7, is the Coma cluster (Thierbach et al., 2003). Importantly, it has also been shown that most of this steepening is not due to the Sunyaev-Zel’dovich effect (SZ) decrement (Brunetti et al., 2013). Other halos with well sampled spectra include the Toothbrush and Bullet cluster which show power-law spectral shapes (Liang et al., 2000; van Weeren et al., 2012b; Shimwell et al., 2014).

There is some evidence that the integrated spectra of radio halos show a correlation with the global ICM temperature of clusters, where hotter clusters host halos with flatter spectra (Feretti et al., 2004a; Giovannini et al., 2009). However, Kale & Dwarakanath (2010) pointed out that comparing the average values of ICM temperatures and of spectral indices can give inconclusive results.

5.1.4 Resolved spectra

The first detailed study of the spatial distribution of the radio spectral index across a radio halo was carried out by Giovannini et al. (1993). They found a smooth spectral index distribution for the Coma cluster radio halo, with evidence for radial spectral steepening. For Abell 665 and Abell 2163 hints of radial spectral steepening where also found in undisturbed cluster regions (Feretti et al., 2004b). A caveat of these studies is that they were not done with matched uv-coverage, which could lead to errors in the derived spectral index distributions. Some other studies of radio halo spectral index distributions are Giacintucci et al. (2005); Orrú et al. (2007); Pizzo & de Bruyn (2009); Kale & Dwarakanath (2010); Shimwell et al. (2014); Pearce et al. (2017). Two examples radio halo spectral index maps, for the massive merging clusters Abell 2744 and the Toothbrush, are shown in Figure 8. It shows that the spectral index is rather uniform across these radio halos.

A spatial correlation between radio spectral index and ICM temperature ($T$) for Abell 2744 was reported by Orrú et al. (2007), with flatter spectral index regions corresponding to higher temperatures. However, using deeper VLA and Chandra data this result was not confirmed (Pearce et al., 2017). Similarly, no clear evidence for such a correlation was founding in Abell 520 (Vacca et al., 2014), the Toothbrush Cluster (van Weeren et al., 2016), the Bullet cluster (Shimwell et al., 2014), and Abell 2256 (Kale & Dwarakanath, 2010). The current results therefore indicate there is no strong $T-\alpha$ correlation present, although more studies are necessary. It has been noted that even in the presence of an underlying $T-\alpha$ correlation, projection effects might also significantly reduce the detectability (Kale & Dwarakanath, 2010).

5.1.5 Ultra-steep spectrum radio halos

Some halos have been found that have ultra-steep spectra, up to $\alpha\sim-2$. Radio halos with $\lesssim-1.6$ have been called ultra-steep spectrum radio halos (USSRH). The existence of USSRH is expected if the integrated spectra of radio halos include a cutoff. When we measure the spectral index close to the cutoff frequency ($\nu_{\rm b}$) it becomes very steep. Any radio halo can thus appear as an USSRH as along as we observe it close to (or beyond) the cutoff frequency. It is expected that only the most luminous radio halos, corresponding to the most energetic merger events, have cutoff frequencies of $\gtrsim 1$ GHz. In the turbulent re-acceleration model, the location of the cutoff frequency approximately scales as (Cassano et al., 2010a),

[TABLE]

where $M$ is the mass of the main cluster. In connection with major merger events

[TABLE]

where $\Delta M$ the mass the merging subcluster. Because of these scalings, it is expected that more USSRH radio halos, corresponding to less energetic merger events, can be uncovered with sensitive observations at low frequencies.

The prime example of a USSRH is found in Abell 521 (Brunetti et al., 2008; Dallacasa et al., 2009), Other clusters with USSRH or candidate USSRH are Abell 697 (Macario et al., 2010; van Weeren et al., 2011b; Macario et al., 2013), Abell 2256 (Brentjens, 2008), Abell 2255 (Feretti et al., 1997a; Pizzo & de Bruyn, 2009), Abell 1132 (Wilber et al., 2018b), MACS J0416.1–2403 (Pandey-Pommier et al., 2015), MACS J1149.5+2223 (Bonafede et al., 2012), Abell 1300 (Reid et al., 1999; Venturi et al., 2013), and PSZ1 G171.96–40.64 (Giacintucci et al., 2013). It should be noted that a number of these USSRH still need to be confirmed. The reason is that reliable spectral index measurements are difficult to obtain because of differences in uv-coverage, sensitivity, resolution, and absolute flux calibration. This situation will improve with the new and upgraded radio telescopes that have become operational, in particular at low frequencies. One example of a candidate radio halo with an ultra-steep spectrum was Abell 1914 (Bacchi et al., 2003). Recent LOFAR and GMRT observations suggest that the most of the diffuse emission in this cluster does not come from a halo but instead from a radio phoenix (Mandal et al., 2018).

5.1.6 Polarization

Radio halos are found to be generally unpolarized. This likely is caused by the limited angular resolution of current observations, resulting in beam depolarization. This effect is significant when the beam size becomes larger than the angular scale of coherent magnetic field regions. Even at high-angular resolution, magnetic field reversals and resulting Faraday rotation will reduce the amount of observed polarized flux.

For three clusters, Abell 2255, MACS J0717.5+3745, and Abell 523 significant polarization has been reported (Govoni et al., 2005; Bonafede et al., 2009b; Girardi et al., 2016), but it is not yet fully clear whether this emission is truly from the radio halos, or from polarized cluster radio shocks projected on-top or near the radio halo emission (Pizzo et al., 2011; van Weeren et al., 2017a).

Govoni et al. (2013) modeled the radio halo polarization signal at 1.4 GHz and inferred that radio halos should be intrinsically polarized. The fractional polarization at the cluster centers is about 15–35%, varying from cluster to cluster, and increasing with radial distance. However, the polarized signal is generally undetectable if it is observed with the low sensitivity and resolution of current radio interferometers. The Govoni et al. (2013) results are based on MHD simulations by Xu et al. (2011, 2012) which are probably not accurate enough yet to resolve the full dynamo amplification. Whether this will affect the predicted fractional polarization levels is not yet clear, see Donnert et al. (2018). If the polarization properties of radio halos can be obtained from future observations it would provide very valuable information on the ICM magnetic field structure.

5.1.7 Samples and scaling relations, merger connection

Statistical studies of how the radio halo properties relate to the ICM provide important information on the origin of the non-thermal CR component.

It is well known (e.g., Liang et al., 2000; Enßlin & Röttgering, 2002; Feretti, 2003; Yuan et al., 2015) that the radio power (luminosity) of giant halos correlates with the cluster X-ray luminosity ($L_{\rm{X}}$), and thus cluster mass. For observational reasons, the radio power at 1.4 GHz ($P_{\rm{1.4GHz}}$) is commonly used to study scaling relations. The X-ray luminosity is often reported in the 0.1–2.4 keV ROSAT band. Figure 9 shows a compilation of radio halos and upper limits on a mass-$P_{\rm{1.4GHz}}$ and $L_{\rm{X}}$-$P_{\rm{1.4GHz}}$ diagram. Detailed investigations of the scaling relations between radio power and X-ray luminosity (or mass), based on the turbulent re-acceleration model, were performed by Cassano et al. (2006, 2007, 2008a). These models were also used to predict the resulting statistics for upcoming radio surveys (Cassano et al., 2010a; Cassano, 2010; Cassano et al., 2012). More recently, the integrated Sunyaev-Zel’dovich Effect signal (i.e., the Compton $Y_{\rm{SZ}}$ parameter) has been used as a proxy of cluster mass (Basu, 2012; Cassano et al., 2013; Sommer & Basu, 2014). The advantage from using this proxy stems from the fact that $Y_{\rm{SZ}}$ should be less affected by the dynamical state of a cluster, providing less scatter compared to $L_{\rm{X}}$ (e.g., Motl et al., 2005; Wik et al., 2008).

To determine radio halo power or upper limits for statistical studies, it is important to derive these quantities in a homogeneous way and minimize the dependence on map noise or uv-coverage. This argues against using a certain contour level, often $3\sigma_{\rm rms}$ has been used, to define the radio halo flux density integration area. Assumptions have to be made on the brightness distribution to determine upper limits for non-detections (Brunetti et al., 2007; Murgia et al., 2009; Russell et al., 2011). For example, Bonafede et al. (2017) used an exponential radial profile of the form

[TABLE]

with added brightness fluctuations, with the characteristic sizes ($r_{e}$, e-folding radius) determined from previously found correlations between power and size (Cassano et al., 2007; Murgia et al., 2009). In addition, ellipsoidal profiles were employed for clusters with very elongated X-ray brightness distributions. The effects of uv-coverage, visibility weighting, mosaicking (for observations that combine several pointings), and deconvolution can be quantified by injection of mock radio halos into the uv-data (Brunetti et al., 2007; Johnston-Hollitt & Pratley, 2017).

Radio halos are rather common in massive clusters. An early study by Giovannini et al. (1999) showed that about 6%–9% of $L_{\rm{X}}<5\times 10^{44}$ erg s*-1* clusters host halos at the limit of the NVSS survey, while this number increases to 27%–44% above this luminosity. Extensive work, mainly using the GMRT, provided further improvements on the statistics, showing that the occurrence fraction for clusters with $L_{\rm{X}}>5\times 10^{44}$ erg s*-1* is about 30% (Venturi et al., 2007, 2008; Cassano et al., 2013; Kale et al., 2015). For a mass-selected sample ($M>6\times 10^{14}$ M*⊙), Cuciti et al. (2015) found evidence for a drop in the halo occurrence fraction for lower mass clusters. For clusters with $M>8\times 10^{14}$ M⊙* this fraction is $\approx 60\%-80\%$, dropping to $\approx 20\%-30\%$ below this mass.

An important result from observations is that giant radio halos are predominately found in merger clusters, as indicated by a disturbed ICM and/or other indicators of the cluster’s dynamical state, e.g., the velocity distribution of cluster member galaxies, presence of multiple BCGs, and galaxy distribution. Early work already established evidence that radio halos were related to cluster merger events as determined from X-ray observations (e.g., Feretti et al., 2000; Buote, 2001; Schuecker et al., 2001, 2002; Feretti, 2002; Giovannini & Feretti, 2002; Böringer & Schuecker, 2002). This conclusion is also supported by optical studies (Ferrari et al., 2003; Boschin et al., 2004, 2006; Girardi et al., 2006; Barrena et al., 2007a; Girardi et al., 2008; Boschin et al., 2008, 2009; Girardi et al., 2010, 2011; Boschin et al., 2012b, a; Barrena et al., 2014; Girardi et al., 2016; Golovich et al., 2016). A common method is to use the cluster’s X-ray morphology as an indicator of the cluster’s dynamical state, such as the centroid shift, power ratio, and concentration parameter (Buote, 2001; Cassano et al., 2010b). Almost all giant ($\gtrsim 1$ Mpc) radio halos so far have been found in dynamically disturbed clusters. Recent studies also confirm this general picture (Cassano et al., 2013; Kale et al., 2015; Cuciti et al., 2015), but see Section 5.2.3 for some exceptions.

Further support for the relation between cluster mergers and the presence of radio halos was presented by Brunetti et al. (2009). They found that there is a radio bi-modality between merging and relaxed clusters. Merging clusters host radio halos, with the radio power increasing with $L_{\rm X}$. Relaxed clusters do not show the presence of halos, with upper limits located well below the expected correlation. Similarly, Rossetti et al. (2011); Brown et al. (2011b) find that the occurrence of halos is related to the cluster’s evolutionary stage. Early work by Basu (2012) reported a lack of a radio bimodality in the Y–P plane. However, this was not confirmed by Cassano et al. (2013). On the other hand, X-ray selected cluster samples are biased towards selecting cool core clusters, which generally do not host giant radio halos, and hence the occurrence fraction of radio halos in SZ-selected samples is expected to be higher (Sommer & Basu, 2014; Andrade-Santos et al., 2017). Recently, Cuciti et al. (2018) found two radio halos that occupy the region below the mass-$P_{\rm{1.4GHz}}$ correlation. These two underluminous radio halos do not have steep spectra and could be generated during minor mergers where turbulence has been dissipated in smaller volumes, or be “off-state” radio halos originating from hadronic collisions in the ICM.

Some merging clusters that host cluster double radio shocks (see Section 6.1.2), do not show the presence of a radio halo (Bonafede et al., 2017). This absence of a radio halo might be related to early or late phase mergers, and the timescale of halo formation and disappearance. Although, these results are not yet statistically significant given the small sample size.

Cassano et al. (2016) investigated whether giant radio halos can probe the merging rate of galaxy clusters. They suggested that merger events generating radio halos are characterized by larger mass ratios. Another possible explanation is that radio halos may be generated in all mergers but their lifetime is shorter than the timescale of the merger-induced disturbance. The lack of radio halos in some merging clusters can also be caused by the lack of sufficiently deep observations. One prime example is Abell 2146 (Russell et al., 2011) where no diffuse emission was found in GMRT observations. However, recent deep VLA and LOFAR observations revealed the presence of a radio halo in this cluster (Hlavacek-Larrondo et al., 2018; Hoang et al., 2018a).

5.1.8 Origin of radio halos

The origin of radio halos have been historically debated between two models: the hadronic and turbulent re-acceleration models. In the hadronic model, radio emitting electrons are produced in the hadronic interaction between CR protons and ICM protons (Dennison, 1980; Blasi & Colafrancesco, 1999; Dolag & Enßlin, 2000; Miniati et al., 2001a; Pfrommer et al., 2008; Keshet & Loeb, 2010; Enßlin et al., 2011). In the re-acceleration model, a population of seed electrons (e.g., Pinzke et al., 2017) is re-accelerated during powerful states of ICM turbulence (Brunetti et al., 2001; Petrosian, 2001; Donnert et al., 2013; Donnert & Brunetti, 2014), as a consequence of a cluster merger event. While indirect arguments against the hadronic model can be drawn from the integrated radio spectral (Brunetti et al., 2008) and spatial characteristics of halos, and from radio–X-ray scaling relations (for a review see Brunetti & Jones, 2014), only gamma-ray observations, which will be discussed in more detail below (Section 5.1.9), of the Coma cluster directly determined that radio halos cannot be of hadronic origin. The spatial distribution of spectral indices across radio halos, which can go from being very uniform to more patchy, might provide further tests for turbulent re-acceleration model. Furthermore, additional high-frequency ($\gtrsim 5$ GHz) observations of known radio halos would enable a search for possible spectral cutoffs. Such cutoffs are expected in the framework of the turbulent re-acceleration model, but have so far rarely been observed (see Sections 5.1.3 and 5.1.5). Such measurements would be quite challenging though, requiring single dish observations to avoid resolving out diffuse emission.

Nowadays, turbulent re-acceleration is thought to be the main mechanism responsible for generating radio halos, even if other mechanisms as magnetic reconnection have been proposed (e.g., Brunetti & Lazarian, 2016). However, one of the main open questions for the re-acceleration model is the source of the seed electrons. There are several possibilities, with secondary electrons coming from proton-proton interactions being an obvious candidate (Brunetti & Blasi, 2005; Brunetti & Lazarian, 2011). The seed electrons could also have been previously accelerated at cluster merger and accretion shocks. A third possibility is that the seed electrons are related to galaxy outflows and AGN activity. The latter, in particular, is becoming more and more evident thanks to the recent low-frequency observation of re-energized tails (de Gasperin et al., 2017a, see Section 6.3) and fossil plasma sources (e.g., Shimwell et al., 2016). While it is difficult to determine the possible contribution of these primary sources of seed electrons, gamma-ray observations can be used to study the contribution of secondary electrons. Another important open question in this context is the connection with the generation mechanism for mini-halos that will be discussed in Section 5.2.3.

Eckert et al. (2017) used the amplitude of density fluctuations in the ICM as a proxy for the turbulent velocity. Importantly, they inferred that radio halo hosting clusters have one average and a factor of two higher turbulent velocities. However, this indirect method relies on number of assumptions making the result somewhat open to interpretation. Direct measurements of ICM turbulence have so far only been performed for the Perseus cluster with the Hitomi satellite (Hitomi Collaboration et al., 2016, 2018), finding a line-of-sight velocity dispersion of $164\pm 10$ km s*-1*. Future measurements with XRISM (X-ray Imaging and Spectroscopy Mission) and Athena (Nandra et al., 2013; Barret et al., 2016) of the turbulent motions in halo and non-halo hosting clusters will provide crucial tests for the turbulent re-acceleration model.

5.1.9 Gamma-ray upper limits

Gamma-rays in clusters of galaxies are expected from neutral pion decays coming from proton-proton interactions (for more details see Reimer, 2004; Blasi et al., 2007; Pinzke et al., 2011). As mentioned earlier, CR protons can be injected in clusters by structure formation shocks and galaxy outflows, and can accumulate there for cosmological times. The quest for the detection of these gamma-rays have been going on for about two decades now (Reimer et al., 2003; Reimer & Sreekumar, 2004; Aharonian et al., 2009; Ackermann et al., 2010; Aleksić et al., 2010; Arlen et al., 2012; Huber et al., 2012, 2013; Ackermann et al., 2014; Zandanel & Ando, 2014; Prokhorov & Churazov, 2014; Griffin et al., 2014; Ackermann et al., 2016; Liang et al., 2016; Branchini et al., 2017). Unfortunately, the detection of diffuse gamma-ray emission connected with the ICM has been so far elusive. There is no conclusive evidence for an observation yet.

Nevertheless, gamma-ray observations have been very important in the last few years for three reasons: to put a direct limit on the CR content in clusters, to test the hadronic nature of radio halos and mini-halos, and to test the contribution of secondary electrons in re-acceleration models. The number of works on this topic are numerous, thanks to the observations of imaging atmospheric Cherenkov telescopes and of gamma-ray satellites, and the most relevant ones have been cited in the previous paragraph.

Of particular importance for this review are the observations of Coma and Perseus clusters (results for the Perseus cluster will be discussed in Section 5.2.4), and of larger combined samples of nearby massive and X-ray luminous clusters. The combined likelihood analysis of the Fermi-Large Area Telescope (LAT; Atwood et al., 2009) satellite of 50 HIFLUGCS clusters have been a milestone in constraining the amount of CR protons in merging clusters to be below a few percent (Ackermann et al., 2014). However, the most constraining object is the Coma cluster due to its high mass, closeness and radio-halo brightness. In fact, thanks to the Fermi-LAT observations, we are now able to exclude the hadronic origin of the prototypical radio halo of Coma independently from the exact magnetic field value in the cluster (Brunetti et al., 2012, 2017), a long standing issue in the field (e.g., Jeltema & Profumo, 2011). In particular, the CR-to-thermal energy in Coma is limited to be $\lesssim 10$%, almost independently (within a factor or two) from the specific model considered, i.e., re-acceleration or hadronic, and from the magnetic field (Brunetti et al., 2017). Additionally, the Fermi-LAT observations of Coma are starting to test re-acceleration models. These first gamma-ray constraints on re-acceleration are obtained under the assumption that only CR protons and their secondaries are present in the ICM (Brunetti et al., 2017). While we obviously know that this is not the case (see the discussion in the previous Sec. 5.1.8), it is possible that CR protons and their secondaries give the dominant seed contribution.

5.1.10 Radio halo-shock edges

In a handful of clusters the radio halo emission seems to be bounded by cluster shock fronts (Markevitch et al., 2005; Brown & Rudnick, 2011; Markevitch, 2010; Planck Collaboration et al., 2013; Vacca et al., 2014; Shimwell et al., 2014; van Weeren et al., 2016). Two of these examples of “halo-shock edges” are shown in Figure 10. The nature of these sharp edges is still unclear.

It is possible that some of the “halo” emission near these shocks comes from CR electrons compressed at the shock. Alternatively, these edges are cluster radio shocks where electrons are (re-) accelerated. When these electrons move further downstream they will be re-accelerated again, but now by turbulence generated by the merger. Then, depending on the observing frequency, magnetic field strength (which sets the cooling time), and timescale for the turbulent cascade and re-acceleration, the radio shock and halo emission might blend forming these apparent halo-shock edges.

On the other hand, so far no polarized emission has been observed at these halo-shock edges (Shimwell et al., 2014) which would indicate compression. Also, no clear strong downstream spectral gradients due to electron energy losses have been found so far (e.g., van Weeren et al., 2016; Rajpurohit et al., 2018; Hoang et al., 2018c). If the synchrotron emission purely comes from a second order Fermi process at these edges, it would imply that there is sufficient post-shock MHD turbulence immediately after the shock (see for example Fujita et al., 2015). However, if this turbulence is generated by the shock passage downstream there might be insufficient time for this turbulence to decay to the smaller scales that are relevant for particle acceleration. To fully understand the nature of halo-shock edges, future high-resolution spectral and polarimetric observations will be crucial.

5.2 Mini-halos

Radio mini-halos have sizes of $\sim$100–500 kpc and are found in relaxed cool core clusters, with the radio emission surrounding the central radio loud BCG (for a recent overview of mini-halos see Gitti et al., 2015). The sizes of mini-halos are comparable to that of the central cluster cooling regions. The prototypical mini-halo is the one found in the Perseus cluster (Miley & Perola, 1975; Noordam & de Bruyn, 1982; Pedlar et al., 1990; Burns et al., 1992; Sijbring, 1993; Sijbring & de Bruyn, 1998), see Figures 11 and 12. Although smaller than radio halos, radio mini-halos also require in-situ acceleration given the short lifetime of synchrotron emitting electrons. The radio emission from mini-halos does therefore not directly originate from the central ANG, unlike the radio lobes that coincide with X-ray cavities in the ICM.

Radio mini-halos have 1.4 GHz radio powers in the range of $10^{23}-10^{25}$ W Hz*-1*. The most luminous mini-halos known are located in the clusters PKS 0745–191 (Baum & O’Dea, 1991) and RX J1347.5–1145 (Gitti et al., 2007), although the classification of the radio emission in PKS 0745–191 as a mini-halo is uncertain (Gitti et al., 2004; Venturi et al., 2007). The most distant mini-halo is found in the Phoenix Cluster (van Weeren et al., 2014), although very recently a possible mini-halo in ACT-CL J0022.2–0036 at $z=0.8050$ has been reported by Knowles et al. (2018).

Compared to giant radio halos, the synchrotron volume emissivities of mini-halos are generally higher (Cassano et al., 2008b; Murgia et al., 2009). Murgia et al. (2009) fitted exponential azimuthal surface brightness profiles (see Equation 11) and showed that mini-halos have smaller e-folding radii ($r_{e}$) compared to giant halos, as expected from their smaller sizes with the emission being mostly confined to the X-ray cooling region.

Since the mini-halo emission surround the central radio galaxy, whose lobes often have excavated cavities in the X-ray emitting gas, the separation between AGN lobes and mini-halos can be difficult, in particular in the absence of high-resolution images. Radio emission that directly surrounds the central AGN (less than a few dozens of kpc), does not necessarily require in-situ re-acceleration. This emission has also been classified as ‘core-halo’ sources. The separation between core-halo sources, amorphous lobe-like structures, and mini-halos is often not clear (Baum & O’Dea, 1991; Mazzotta & Giacintucci, 2008). In addition, the central radio galaxies are sometimes very bright, requiring high-dynamic range imaging to bring out the low-surface brightness mini-halos. The classification as a mini-halo is also difficult without X-ray data (e.g., Bagchi et al., 2009). Because of these observational limitations, there is currently a rather strong observational selection bias. For that reason many fainter radio mini-halos could be missing since they fall below the detection limit of current telescopes. Despite these observational difficulties the number of known mini-halo has steadily been increasing (Gitti et al., 2006; Doria et al., 2012; Giacintucci et al., 2011b, 2014b, 2017). In Table LABEL:tab:clusterlist we list the currently known radio mini-halos and candidates.

An example of a source that is difficult to classify is the one found in the central parts of the cluster Abell 2626. This source was initially named as a mini-halo by Gitti et al. (2004). More detailed studies (Gitti, 2013; Ignesti et al., 2017; Kale & Gitti, 2017) reveal a complex “kite-like” radio structure, complicating the interpretation and classification. The cluster RX J1347.5–1145 presents another interesting case. It was found to host a luminous radio mini-halo (Gitti et al., 2007) with an elongation to the south-east. This elongation seems to correspond to a region of shock heated gas induced by a merger event, also detected in the SZ (Komatsu et al., 2001; Kitayama et al., 2004; Mason et al., 2010; Korngut et al., 2011; Johnson et al., 2012). This suggests that the south-east emission is not directly related to the central mini-halo, but rather is a separate source (Ferrari et al., 2011) which could be classified as a cluster radio shock.

Few detailed high-quality resolved images of mini-halos exist. This makes it hard to study the morphology of mini-halos in detail. Interestingly, Mazzotta & Giacintucci (2008) found that mini-halos are often confined by the cold fronts of cool core clusters (but see Section 5.2.3). The most detailed morphological information is available for the Perseus cluster mini-halo. Gendron-Marsolais et al. (2017) presented 230–470 MHz images which revealed filamentary structures in this mini-halo, extending in various directions (Figure 12). Hints of these structures are already visible at 1.4 GHz (Sijbring et al., 1989). These structures could be related to variations in the ICM magnetic field strength, localized sites of particle re-acceleration, or a non-uniform distribution of fossil electrons. The Perseus cluster mini-halo emission also follows some of the structures observed in X-ray images. Most of the mini-halo emission is contained within a cold front. However, some faint emission extends (“leaks”) beyond the cold front. Similarly, the RX J1720.1+2638 mini-halo also displays substructure suggesting that when observed at high resolution and signal-to-noise mini-halos are not fully diffuse.

Spectral indices of radio mini-halos are similar to giant radio halos, although few detailed studies exist. The integrated spectrum for the Perseus mini-halo is consistent with a power-law shape (Sijbring, 1993). A hint of spectral steepening above 1.4 GHz is found for RX J1532.9+3021 (Hlavacek-Larrondo et al., 2013; Giacintucci et al., 2014b). An indication of radial spectral steepening for the Ophiuchus cluster (Govoni et al., 2009; Pérez-Torres et al., 2009) was reported by Murgia et al. (2010a). The most detailed spectral study so far has been carried out on RX J1720.1+2638 (Mazzotta & Giacintucci, 2008; Giacintucci et al., 2014a). This mini-halo shows a spiral-shaped tail, with spectral steepening along the tail. Possible steepening of the integrated spectrum for RX J1720.1+2638 at high frequencies has also been reported (Giacintucci et al., 2014a). So far no targeted polarization studies of mini-halos have been performed.

5.2.1 Statistics

Giacintucci et al. (2014b) found no clear correlation between the mini-halo radio power and cluster mass, unlike giant radio halos. However, Cassano et al. (2008b); Kale et al. (2013); Gitti et al. (2015) did report evidence for a correlation between radio power and X-ray luminosity. The slope of the correlation was found to be similar to that of giant radio halos (Gitti et al., 2015). Larger samples are required to obtain better statistics and confirm the found correlations, or lack thereof.

Giacintucci et al. (2017) determined the occurrence of radio mini halos in a sample of 58 clusters with M${}_{500}>6\times 10^{14}$ M⊙. They found that 80% of the cool core clusters hosted mini-halos. Therefore, mini-halos are common phenomenon in such systems. No mini-halos were found in non-cool core systems. In addition, tentative evidence was found for a drop in the occurrence rate for lower cluster masses. Kale et al. (2013) found a mini-halo occurrence rate of about 50% in the Extended GMRT Radio Halo Survey ($L_{\rm{X,0.1-2.4keV}}>5\times 10^{44}$ erg s*-1*, $0.2<z<0.4$), also indicating mini-halos are rather common.

5.2.2 Origin of radio mini-halos

Similar to giant radio halos, hadronic (e.g., Pfrommer & Enßlin, 2004) or turbulent re-acceleration models (Gitti et al., 2002) been invoked to explain the presence of the CR synchrotron emitting electrons. Unlike giant radio halos, where the turbulence is induced by major cluster mergers, mini-halos would trace turbulence in the cluster cores generated by gas sloshing (ZuHone et al., 2013, 2015). The central AGN is a likely candidate for the source of the fossil electrons that are re-accelerated (e.g., Fujita et al., 2007). The confinement of mini-halos by cold fronts (Mazzotta & Giacintucci, 2008) support a scenario where turbulence induced by gas sloshing motions re-accelerates particles. Simulations by Fujita & Ohira (2013); ZuHone et al. (2013, 2015) provided further support for this scenario, reproducing some of the observed morphology, where the emission is bounded by cold fronts.

The radio spectral properties of mini-halos provide another discriminator for the origin of the CR electrons. If the electrons are re-accelerated by magnetohydrodynamical turbulence, the integrated spectra of mini-halos should display a spectral break caused by a cutoff in the electron energy distribution. Due to the limited number of spectral studies available, no clear conclusion can be drawn on the general occurrence of spectral breaks in mini-halo spectra.

5.2.3 Unification

Despite of the their differences, it is possible that mini-halos and giant halos in clusters are physically related to each other. For example, cluster merger events could transport CR from cluster cores to larger-scales where they are re-accelerated again (see Brunetti & Jones, 2014). This could lead to “intermediate” cases where mini-halos could evolve into giant radio halos and vice-versa. This could either be a transition been turbulent re-acceleration due to core sloshing and merger induced turbulent re-acceleration. Or alternatively, a transition between hadronic mini-halos and merger induced turbulent re-acceleration (Zandanel et al., 2014). Recent observations have provided evidence for such scenarios, finding (mini-)halos with unusual properties.

Bonafede et al. (2014b) discovered a large 1.1 Mpc radio halo in CL1821+643 which contains a strong cool core. If this halo is caused by a merger event, the cluster is in a stage where the merger has not (yet) been able to disrupt the cool core as also noted by Kale & Parekh (2016). For example, because the merger is an off-axis event, or the merger is still in an early stage. CL1821+643 could therefore be a transitional object, where a mini-halo is switching off and a giant radio halo is just being formed. Similarly, Sommer et al. (2017); Savini et al. (2018a) reported the presence of a $\sim$1 Mpc radio halo in the semi-relaxed cluster Abell 2261131313The classification of Abell 2390 as giant radio halo by Sommer et al. (2017) was not confirmed by Savini et al. (2018a) which suggested the emission belongs to a double lobe radio galaxy., questioning the assumption that giant radio halos only occur in clusters undergoing major mergers.

Another peculiar case is the sloshing, minor-merger cluster Abell 2142. Early work already hinted at the presence of diffuse emission Harris et al. (1977); Giovannini et al. (1999); Giovannini & Feretti (2000) in this cluster. This was confirmed by Farnsworth et al. (2013) which showed a 2 Mpc radio halo. Venturi et al. (2017) found that the radio halo consists of two components. The inner component has a higher surface brightness, with properties similar to that of a mini-halo. The outer larger component has a steeper spectrum. They proposed that the inner component is powered by central sloshing turbulence. The outer component might probe turbulent re-acceleration induced by a less energetic merger event. Alternatively, the different components are the result from a transition between hadronic and turbulent re-acceleration processes.

The cluster PSZ1 G139.61+24.20 ($z=0.267$) was listed as a candidate mini-halo by Giacintucci et al. (2017). Savini et al. (2018b) presented the discovery of steep-spectrum emission extending beyond the cool core region of the cluster with LOFAR. They argued that the emission outside the core is produced by turbulent re-acceleration from a minor merger event that has not disrupted the cool core. If this scenario is correct, it indicates that both a giant radio halo and mini halo could co-exist. A very similar situation has recently been found in the cluster RX J1720.1+2638. Here, Savini et al. (2018a) discovered extended faint diffuse steep spectrum emission beyond the cold front and mini-halo region (Mazzotta & Giacintucci, 2008; Giacintucci et al., 2014a).

5.2.4 Gamma-ray upper limits

The most important gamma-ray limits on mini-halos come from the observations of the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) telescopes of Perseus (Aleksić et al., 2010, 2012; Ahnen et al., 2016), and from the combined likelihood analysis of HIFLUGCS clusters with the Fermi-LAT satellite data (Ackermann et al., 2014). As is the case of Coma for merging clusters, Perseus is the most constraining object when it comes to mini-halos because of its high mass, closeness, and mini-halo brightness. Perseus host two gamma-ray bright AGNs - the central radio galaxy NGC 1275 and IC 310 - detected both by Fermi (Abdo et al., 2009) and by MAGIC (Aleksić et al., 2012; Ahnen et al., 2016). The poor angular resolution of Fermi at low ($<$10 GeV) energies makes it difficult to target the possible diffuse gamma-ray emission in Perseus, and makes the MAGIC Perseus observations the most constraining for relaxed cool core clusters hosting mini-halos.

Differently from the case of the Coma radio halo, the gamma-ray upper limits on Perseus do not yet allow to exclude the hadronic origin of its mini-halo. The CR energy density in Perseus is constrained to be below about 1–10% of the thermal energy density, with the exact number depending on the assumptions made regarding the CR-spectral and spatial distribution, e.g., the steeper the spectrum and/or the flatter spatial (radial) distribution, the looser the constraints become. This strong dependence of the constraints on the CR content in clusters on the proton spectral and spatial distributions should be kept in mind when quoting these limits.

Assuming the mini-halo emission is hadronic, the gamma-ray upper limit can be turned in a lower limit on the magnetic field needed to generate the radio emission with secondary electrons. This is similar to what has been done for the Coma radio halo where the magnetic field needed for the hadronic interpretation would be extremely high with an energy density of 1/3 or more of the thermal energy density (Brunetti et al., 2017). In the case of Perseus, current gamma-ray limits imply central magnetic fields above $\sim 5$ $\mu$Gauss, still well below the $\sim 25$ $\mu$Gauss inferred from Faraday rotation measurements (Taylor et al., 2006).

5.3 Future Gamma-ray studies

Future gamma-ray observations of clusters of galaxies will be fundamental for this field as only thanks to gamma-rays the exact amount of CR protons can be directly studied and the degeneracy between secondary and primary sources of electrons in radio-halo models can be addressed. In particular, future observations of the Perseus cluster – as envisioned in the key science projects of the in-construction Cherenkov Telescope Array (CTA; Cherenkov Telescope Array Consortium et al., 2017) – will allow to eventually test the hadronic interpretation of mini-halos, and, more importantly, to limit the CR energy density to below about 2% independently from the assumptions on the CR-proton spectral and spatial distribution. Such low limits will also allow to test the role of AGNs, particularly, the protons confinement in AGN bubbles and how protons are transported from the central AGNs to cluster peripheries. Paramount for an “order-of-magnitude” jump in constraining power, also for gamma-ray observations of cluster radio shocks, and hopefully to aim for several detections, will be the satellites proposed as successors of Fermi (Knödlseder, 2016; McConnell, 2016). Note, however, that if nature is “kind” and the electrons generating the radio halo of Coma are re-accelerated secondaries, continued Fermi observations could reach a detection in the near future (Brunetti et al., 2017).

5.4 Upcoming large cluster samples

With new deep low-frequency radio surveys covering a significant fraction of the sky (such as LoTSS, Shimwell et al., 2017, 2018) many new radio (mini-)halos are expected to be discovered. In particular those with steep radio spectra. New surveys are also planned at GHz frequencies (Norris et al., 2011, 2013) which should also uncover additional diffuse cluster radio sources (Cassano et al., 2012).

With the improved statistics offered by larger samples, the properties and occurrence rates as a function of cluster mass, dynamical state, and other global cluster properties can be investigated in detail. These samples should also contain a population of ultra-steep spectrum radio halos that are predicted in the framework of the turbulent re-acceleration model (Cassano et al., 2010a). Furthermore, large samples might shed more light on (i) the possible connections between halos and mini-halos and (ii) the evolution of diffuse cluster radio sources over cosmic time, from $z\sim$1 to the present epoch. For example, changes in the occurrence rate are expected due to the increase of Inverse Compton losses with redshift, change of the cluster merger rate, evolution of cluster magnetic fields.

6 Cluster radio shocks (relics) and revived fossil plasma sources

Apart from radio halos, we broadly divide diffuse cluster sources into cluster radio shocks and revived fossil plasma sources. The distinction between radio shocks and fossil plasma sources is not always straightforward, since it requires the detection of shocks via SZ or X-ray measurements and the availability of radio spectra. Our adopted classification is similar to that of Kempner et al. (2004)141414We do not consider dying radio galaxies here that have not interacted with the ICM, see Murgia et al. (2011). who defined radio Gischt and Phoenix. Given that there is now compelling evidence that Gischt trace shock waves (e.g., Finoguenov et al., 2010), we propose to simply call these cluster radio shocks. This still leaves open the questions of the underlying (re-)acceleration mechanism that produces the synchrotron emitting CR at these shocks.

Radio shocks and fossil sources are detected in clusters covering a wide range in mass, unlike radio halos which are almost exclusively found in massive systems. Some examples of radio shocks and revived fossil sources in lower mass clusters are discussed in Subrahmanyan et al. (2003); Kale et al. (2017); de Gasperin et al. (2017b); Brüggen et al. (2018); Dwarakanath et al. (2018). Phoenices and other revived AGN fossil sources (such as GreEt) are characterized by their steep radio spectra and presence of high frequency spectral breaks. These sources will be discussed in Section 6.3.

Similar to giant radio halos and mini-halos, there are “hybrid” or “intermediate” sources which share some properties between these two categories. For example AGN fossil plasma that is re-accelerated at a large cluster merger shock (e.g., in Abell 3411-3412; van Weeren et al., 2017b).

6.1 Cluster radio shocks (relics)

Cluster radio shocks are mostly found in the outskirts of galaxy clusters, see Figure 13. Unlike radio halos, they have elongated shapes. In addition, radio shocks are strongly polarized at frequencies $\gtrsim 1$ GHz, with polarization fractions of $\gtrsim 20\%$ (Enßlin et al., 1998), see Figures 17 and  18.

The first identified cluster radio shock was the source 1253+275 in the Coma cluster (Jaffe & Rudnick, 1979; Ballarati et al., 1981). This radio source has been studied in considerable detail early on by Giovannini et al. (1985, 1991). Recently, evidence for a shock at this location has also been obtained (Ogrean & Brüggen, 2013; Akamatsu et al., 2013), see also Section 6.1.5. A couple of other cluster radio shocks that were studied after the discovery of 1253+275 were the ones found in Abell 2256 (e.g., Bridle & Fomalont, 1976; Röttgering et al., 1994) and Abell 3667 (e.g., Röttgering et al., 1997). The number of detected radio shocks increased significantly with the availability of the NVSS and WENSS surveys (Giovannini et al., 1999; Kempner & Sarazin, 2001). A list of cluster radio shocks in given in Table LABEL:tab:clusterlist. The most powerful cluster radio shock is found in MACSJ0717.5+3745 (Bonafede et al., 2009b; van Weeren et al., 2009c). Interestingly, this cluster also hosts the most powerful radio halo. The most distant radio shocks are located in “El Gordo” at $z=0.87$ (Menanteau et al., 2012; Lindner et al., 2014; Botteon et al., 2016b).

In an idealized binary merger, ‘equatorial’ shocks form first and move outwards in the equatorial plane, see Figure 14). After the dark matter core passage, two ‘merger’ shocks lunch into the opposite directions along the merger axis, which can explain the formation of cluster double radio shocks in observed merging clusters (e.g., van Weeren et al., 2011a, c; Molnar & Broadhurst, 2017), see also Section 6.1.2. Vazza et al. (2012) investigated why cluster radio shocks are mostly found in the periphery of clusters using simulations. They showed that the radial distribution of observed radio shocks can be explained by the radial trend of dissipated kinetic energy in shocks, which increases with cluster centric distance up until half of the virial radius. Analyzing the properties of shocks associated with synthetic merging clusters in structure formation simulations, Ha et al. (2018) found that the CR production peaks at $\sim$1 Gyr after the core passage, with the shock-kinetic-energy-weighted Mach number $\left<\mathcal{M}_{s}\right>_{\phi}\simeq 2-3$ and the CR-flux-weighted Mach number $\left<\mathcal{M}_{s}\right>_{\rm CR}\simeq 3-4$. Simulations by Skillman et al. (2011); Vazza et al. (2012); Bonafede et al. (2012); Skillman et al. (2013); Vazza et al. (2016); Nuza et al. (2017); Wittor et al. (2017) also produce large-scale radio shock morphologies that provide a reasonable match to what is found in observations.

Some examples of studies showing the connection between radio shocks and cluster mergers using optical spectroscopy and imaging are Barrena et al. (2007b); Boschin et al. (2010); Barrena et al. (2009); Boschin et al. (2013); Dawson et al. (2015); Golovich et al. (2017a, 2016); Benson et al. (2017). This connection is also corroborated by weak lensing studies that reveal multiple mass peaks in some radio shock hosting clusters (e.g., Okabe et al., 2015; Jee et al., 2016, 2015). The most comprehensive analysis of a sample of 29 radio shock hosting clusters was performed by Golovich et al. (2017b, 2018). They found that the merger axes of radio shock hosting clusters are generally in or near the plane of the sky. This indicates that there are selection biases for finding cluster radio shocks based on the viewing angle. Due to this selection effect, many radio shocks with less favorable orientations are probably missing in current samples.

Cluster radio shocks seems to be less common than radio halos or mini-halos, the occurrence of radio shocks was found to be about about $5\%\pm 3\%$ by Kale et al. (2015). However, unlike radio halos or mini-halos, the merger axis orientation probably plays an important role in detecting these sources, as mentioned.

Some giant cluster radio shocks such as the Sausage and the Toothbrush are thought to be associated with major mergers with a subclump mass ratio $\lesssim 3$ (Okabe et al., 2015; Jee et al., 2015, 2016), while the cluster ZwCl 0008.8+5215 with a double radio shock and PLCK G287.0+32.9 with multiple radio shocks are merging systems with a mass ratio $\gtrsim 5$ (Golovich et al., 2017a; Finner et al., 2017).

In a few clusters the emission from the cluster radio shocks is attached or overlaps with that of the radio halo (e.g., Dallacasa et al., 2009; van Weeren et al., 2016). The nature of these “bridges” between halos and cluster radio shocks is still unclear. In some cases, the radio halo emission covers the entire region between double radio shocks (Bonafede et al., 2012; Hoang et al., 2017; Di Gennaro et al., 2018). One possibility is that we observe a transition from first order Fermi (re-)acceleration to second order re-acceleration by turbulence that develops in the post shock region.

6.1.1 Morphology and sizes

Cluster radio shocks typically have elongated shapes, examples are the sources found in the Coma cluster (Giovannini et al., 1991), CIZA J2242.8+5301 (van Weeren et al., 2010), Abell 3667 (Röttgering et al., 1997; Johnston-Hollitt, 2003), Abell 115 (Govoni et al., 2001b), and Abell 168 (Dwarakanath et al., 2018). These elongated shapes are expected for sources that trace shock waves in the cluster outskirts and are seen close to edge-on. Examples of radio shocks that are less elongated are found in Abell 2256 (e.g., Clarke & Enßlin, 2006) and ZwCl 2341.1+0000 (Bagchi et al., 2002; van Weeren et al., 2009d). Cluster radio shocks have sizes that roughly range between 0.5 to 2 Mpc, see Figure 15. Most large radio shocks that are found in the cluster outskirts show asymmetric transverse brightness profiles, with a sharp edge on the side away from the cluster center. On the side of the cluster center, the emission fades more gradually, see Figures 17 and  25.

Deep high-resolution observations of large elongated radio shocks have also revealed a significant amount of filamentary substructures, see Figures 24 and 25. Large radio shocks that display these filamentary structures are found in Abell 2256, CIZA J2242.8+5301, Toothbrush, MACS J0717.5+3745, Abell 3376, and Abell 3667. The nature of the filamentary structures is not fully understood. One possibility is that they trace changes in the magnetic field. Alternatively, they reflect the complex shape of the shock surfaces. The filamentary morphology of cluster radio shocks seems to be ubiquitous because all radio shocks that have been studied with good signal to noise and at high resolution display them.

6.1.2 Cluster double radio shocks

A particular interesting class of cluster radio shocks are so-called “double shocks”. Here two large elongated convex radio shocks are found diametrically with respect to the cluster center, see Figure 14. The radio shocks are oriented perpendicular with the respect to the elongated ICM distribution (and merger axis) of the cluster, see Figures 16 and 19. Double radio shocks are an important subclass of radio shocks as the cluster merger scenario can be relatively well constrained. In addition, these system seems to be observed close to edge-on. Note that we reserve the classification of a double radio shock for a pair of shock waves that were generated at the same time during core passage. So the presence of two radio shocks in a cluster alone is not a sufficient condition to classify it as a double radio shocks.

About a dozen well-defined double radio shock systems are known, see Table LABEL:tab:clusterlist. The first cluster double radio shock was found in Abell 3667 (Röttgering et al., 1997). It was realized by Roettiger et al. (1999); Johnston-Hollitt et al. (1999) that these radio sources could have resulted from particles accelerated at shocks from a binary merger event. The presence of a shock in the ICM at the location of the northwestern radio source in Abell 3667 was confirmed via X-ray observations by Finoguenov et al. (2010). The second double radio shock system was discovered by Bagchi et al. (2006) in Abell 3376. Other well studied cluster double radio shocks are the ones in CIZA J2242.8+5301 (van Weeren et al., 2010), ZwCl 0000.8+5215 (van Weeren et al., 2011c), MACS J1752.0+4440 (van Weeren et al., 2012a; Bonafede et al., 2012), PSZ1 G108.18-11.53 (de Gasperin et al., 2015a), and Abell 1240 (Kempner & Sarazin, 2001; Bonafede et al., 2009a).

6.1.3 Radio spectra

The integrated radio spectra of cluster radio shocks display power-law shapes (but see Section 6.1.8), with spectral indices ranging from about $-1.0$ to $-1.5$ (e.g., Bonafede et al., 2012; Feretti et al., 2012; de Gasperin et al., 2014). One notable exception of a flatter integrated spectrum, with good data available, is Abell 2256 where the spectral index is about $-0.8$ (Brentjens, 2008; van Weeren et al., 2012c; Trasatti et al., 2015). This flat spectral index is difficult to reconcile with particle acceleration models and electron energy losses, see van Weeren et al. (2012c) for a discussion. Another example appeared to be ZwCl 2341.1+0000 (van Weeren et al., 2009d) but more recent observations indicate that the spectral index is within the normally observed range (Giovannini et al., 2010; Benson et al., 2017).

Cluster radio shocks often show a clear spectral index gradient across their width, see Figures 17 and 19. The region with the flattest spectral index is located on the side away from the cluster center. Towards the cluster center the spectral index steepens. This steepening is thought to be caused by synchrotron and IC losses in the shock downstream region. The majority of well-studied cluster radio shocks, both single shocks (see Figure 8) and double shocks, show this behavior.

6.1.4 Polarization

Cluster radio shocks are amongst the most polarized sources in the extragalactic sky. Very elongated radio shocks usually show the highest polarization fraction, which is expected if they trace edge-on shock waves (Enßlin et al., 1998). For example, CIZA J2242.8+5301 shows polarization fractions of $\sim$50% or more at GHz frequencies for some parts of the radio shock (van Weeren et al., 2010), see Figure 17.

For large cluster radio shocks the intrinsic polarization angles, corrected for the effect of Faraday Rotation, are found to be well aligned. The polarization magnetic field vectors are oriented within the plane of the radio shock (e.g., Bonafede et al., 2009a; van Weeren et al., 2010; Bonafede et al., 2012; Pearce et al., 2017, see also Figures 17 and 18). Only a few Faraday rotation studies have been performed so far of radio shocks. They indicate that for radio shocks projected at large cluster centric radii the Faraday Rotation is mostly caused by the galactic foreground. Faraday Rotation caused by the cluster can be seen for (parts of) radio shocks at smaller cluster centric radii (Bonafede et al., 2009b; Pizzo et al., 2011; van Weeren et al., 2012b; Owen et al., 2014). From the limited studies available, it seems that large cluster radio shocks strongly depolarize at frequencies $\lesssim 1$ GHz (Brentjens, 2008; Pizzo et al., 2011; Ozawa et al., 2015). Therefore, high-frequency observations (above $\gtrsim 2$ GHz) are best suited to probe the intrinsic polarization properties of radio shocks. For example, the fractional linear polarization in for main ‘Sausage’ and ‘Toothbrush’ radio shocks is on average 40% at 5–10 GHz, reaching 70% in localized areas (Kierdorf et al., 2017; Loi et al., 2017).

6.1.5 Comparison between radio and X-ray observations of ICM shocks

Because of their shapes, locations, and spectral and polarimetric properties, cluster radio shocks are considered to trace particles accelerated at shocks. These shocks can be generated by cluster merger activity or accretion flows from surrounding large-scale structures (e.g., Enßlin et al., 1998). If this assumption is correct, shock waves should coexist at the location of radio shocks. From X-ray observations, the intensity of shock structure can be estimated from the Rankine-Hugoniot jump condition (Landau & Lifshitz, 1959). Assuming a ratio of specific heats as $\gamma=5/3$, we have

[TABLE]

[TABLE]

where the subscripts 1 and 2 refer to the pre- and post-shock ICM density ($\rho$) or temperature ($T$), respectively. The ratios of ICM properties as a function of the shock strength ($\cal M$) are shown in Figure 20. On the other hand, based on the assumption of simple DSA theory, the Mach number can be also estimated from the radio injection spectral index ($\alpha_{\rm inj}$) via

[TABLE]

In principle, both X-ray and radio approaches are independent methods to characterize the shock strength, meaning shock strengths inferred from these different wavelength regimes should match each other, if underlaying assumptions are correct. Therefore, the comparison of the shock properties inferred from X-ray and radio data is an important tool to investigate shock related ICM physics. Until recently, observational information of radio shocks at X-ray wavelengths were limited because radio shocks are typically located in the cluster periphery, where the ICM X-ray emission is very faint. This makes it challenging to characterize the X-ray shock properties.

The first detection of a shock wave, co-located with a cluster radio shock (relic), was in the nearby merging cluster Abell 3667 using XMM-Newton observations. Finoguenov et al. (2010) found a sharp X-ray surface brightness discontinuity at the outer edge of the radio shock, and a significant drop in the ICM temperature at the same location. These discontinuities are consistent with a $\mathcal{M}\sim$2 shock. These results have been confirmed by Akamatsu et al. (2012b); Sarazin et al. (2016).

The Suzaku satellite, being in a low orbit within Earth’s magnetopause, provided a significantly lower and stable particle background compared to previous X-ray observatories (XMM-Newton and Chandra). The low and well calibrated instrumental background of Suzaku made it well-suited to study the faint cluster periphery. The first systematic Suzaku investigation of cluster radio shocks was performed by Akamatsu & Kawahara (2013). Since the first detection of the shock wave associated with Abell 3667, there are about 20 X-ray detected shocks corresponding to radio shocks. An overview of radio shocks with X-ray detections is shown in Table 1.

At radio wavelength, there are also observational challenges to derive shock properties. One particular difficulty is to measure $\alpha_{\rm inj}$. The integrated spectral index of a radio shock reflects a balance between acceleration and energy losses. As a result, the index of the integrated spectrum is 0.5 steeper compared to $\alpha_{\rm inj}$. This relation ($\alpha_{\rm int}=\alpha_{\rm inj}+0.5$ (Kardashev, 1962)) is however somewhat simplistic, since the shock properties do evolve over time, see also Section 6.2. Alternatively, with spatially resolved spectral index maps one can obtain more reliable measurements of $\alpha_{\rm inj}$, avoiding some of the problems with energy losses. Here one needs to measure the spectral index as close as possible to the shock location. However, even in this case some mixing of different electron energy populations will occur, depending on the spatial resolution, shape of the shock surface, and projection effects.

For the northern radio shock in CIZA J2242.8+5301 a number of detailed comparison between the radio and X-ray derived Mach numbers have been performed. van Weeren et al. (2010) reported a radio injection spectral index of $-0.60\pm 0.05$ resulting in $\mathcal{M}=4.6^{+1.3}_{-0.9}$ (68% confidence range). In the X-rays, Akamatsu & Kawahara (2013); Akamatsu et al. (2015) reported a temperature increase across the radio shock with an amplitude of a factor $\sim$3 resulting giving ${\cal M}=2.7^{+0.7}_{-0.4}$ (including systematics due to the background estimation). This kind of tension, ${\cal M}_{\rm radio}>{\cal M}_{\rm X}$, has been found for other radio shocks, see Figure 22. If this discrepancy is indeed real, this may point to problems in the DSA scenario for shocks in clusters. To explain the observational results, several solutions have been proposed.

For example, it is possible that the X-ray derived Mach numbers are somewhat underestimated due to unfavorable viewing angles and the complexity of the shock surface. In addition, the shock acceleration efficiency is a thought to be a strong function of shock Mach number (Hoeft & Brüggen, 2007). Therefore the CR-energy-weighted Mach number is expected to be higher than the kinetic-energy-weighted Mach number (Ha et al., 2018). Thus radio measured Mach numbers will be biased towards parts of the shock with the highest Mach numbers. Difficulties and possible biases with radio based measurements are discussed in Stroe et al. (2014a); van Weeren et al. (2016); Hoang et al. (2017). The re-acceleration of fossil plasma has also been invoked, see Section 6.2. Akamatsu et al. (2017) investigated possible systematic errors associated with X-ray observations. We refer the reader to Sect 4.3. in their paper for more details.

Future X-ray satellites, such as Athena (Nandra et al., 2013), will provide precise measurements of cluster merger shocks. This will shed further light on the apparent discrepancy between the Mach numbers derived from radio and X-ray observations. With the improved collecting area with respect to current satellites, the shock properties in faint cluster outskirt can also be determined.

6.1.6 SZ observations

The thermal ICM electrons in galaxy clusters interact with CMB photons through inverse Compton scattering, resulting in the so-called SZ effect (Sunyaev & Zeldovich, 1970). The SZ effect provides a complementary way of studying the ICM and, because of its redshift-independent nature, is particularly powerful at high-redshift where the X-ray surface brightness suffers from significant cosmological dimming. Low-resolution studies measuring the bulk SZ signal have been very successful at selecting large samples of both relaxed and disturbed clusters up to $z\sim$1.5 (e.g., Planck Collaboration et al., 2016; Bleem et al., 2015).

In the last decade, efforts at the very highest radio frequencies (above 90 GHz) have focused on measuring SZ at high spatial resolution with the aim of detecting small scale features in the ICM, such as shocks in merging clusters. Great strides have been made possible by the introduction of high-resolution, large field-of-view instruments such as MUSTANG-2151515http://www.gb.nrao.edu/mustang/ installed on the 100-m Green Bank Telescope ($9^{\prime\prime}$ resolution at 90 GHz) and NIKA/NIKA2161616http://ipag.osug.fr/nika2/Welcome.html on the 30-m IRAM telescope (reaching $10^{\prime\prime}$–$20^{\prime\prime}$ resolution at 150 and 260 GHz). The power of these instruments has already been demonstrated through high-resolution SZ images showing substructure in merging clusters and in the cores of relaxed clusters (Adam et al., 2017; Romero et al., 2018; Adam et al., 2018). For the nearby Coma cluster the Planck satellite has provided resolved SZ images (Planck Collaboration et al., 2013), including the likely detection of two $\mathcal{M}\sim 2$ shocks in the cluster periphery.

Following pioneering work detecting a weak shock in MACS0744+3927 (Korngut et al., 2011), more recent observational work with the Atacama Large Millimeter/submillimeter Array (ALMA)171717http://www.eso.org/public/teles-instr/alma/receiver-bands/ has enabled a direct detection and measurement of a cluster merger shock in ‘El Gordo’ (Basu et al., 2016a). These observations demonstrate great potential for future SZ determinations of shock properties (particularly the Mach number), especially at large cluster-centric distances and high-redshift, where X-ray measurements of the ICM properties become challenging.

6.1.7 Gamma-rays from cluster radio shocks

Apart from (re-)accelerating electrons, shocks should also accelerate protons. For DSA, the number of accelerated protons should be much larger than electrons. Similar to the secondary model for radio halos, these CR protons should collide with the thermal ICM and produce gamma-rays via hadronic interactions.

It has been noted by Vazza & Brüggen (2014); Vazza et al. (2015a, 2016) that the expected gamma-ray emission for DSA shock acceleration at radio shocks is in tension with gamma-ray upper limits for some clusters. This indicates that the relative acceleration efficiency of electrons and protons is at odds with predictions from DSA. Adding the re-acceleration of fossil particles to this prediction does not change this conclusion. One possible explanation for the lack CR protons is that the magnetic field at radio shocks is predominantly perpendicular to the shock normal. Work by Caprioli & Spitkovsky (2014) indicates that the acceleration efficiency of protons is strongly suppressed at such shocks. Simulations by Wittor et al. (2017) indicate this could reduce the tension with the low gamma-ray upper limits.

Recently, claims of gamma-ray emission from the virial shocks around the Coma cluster (Keshet et al., 2017), as well as from a stacking of other clusters (Reiss et al., 2017), have been put forward. We underline, however, that so far these claims have been not been confirmed (Arlen et al., 2012; Zandanel & Ando, 2014; Prokhorov, 2014; Ackermann et al., 2016).

6.1.8 High-frequency studies of radio shocks

Owing to their their steep spectra, radio shocks have been classically observed at relatively low frequencies ($<2$ GHz). In this Section we review the current state-of-the-art high-frequency observations of radio shocks, by focusing on observations above 5 GHz. High-frequency observations pose particular challenges: (i) radio shocks have steep-spectra making them very faint at high-frequencies; (ii) radio interferometers typically have small fields of view at high frequency and thus have difficulty in detecting extended diffuse sources. Until 2014, the highest frequency detection of a radio shocks were in the clusters Abell 521 and MACS J0717.5+3745 at 5 GHz (Giacintucci et al., 2008; Bonafede et al., 2009b). The interest in high-frequency observations of cluster radio shocks, and the number of detections, has grown over the past few years. This interest has been motivated by the study of the injected electrons and their aging mechanism (as discussed for example by Kang, 2016a; Donnert et al., 2016; Kang & Ryu, 2016).

Instruments that helped make progress at high frequencies include interferometers, such as the Arcminute Microkelvin Imager (AMI, 16 GHz), the Combined Array for Research in Millimeter-wave Astronomy (CARMA, 30 GHz) and the VLA (4–10 GHz), and single dish antennas such as Effelsberg (up to 10 GHz) and the Sardinia Radio Telescope (SRT, up to 19 GHz), see Figure 23.

At the moment of writing, six clusters benefit from radio shock detections above 5 GHz: the main radio shocks in the ‘Sausage’ and the ‘Toothbrush’ clusters (both up to 30 GHz, Stroe et al., 2014b, 2016; Kierdorf et al., 2017; Loi et al., 2017), Abell 2256 (at 5 GHz, Trasatti et al., 2015) , the ‘Bullet’ cluster radio shock (5.5 and 9 GHz, Malu et al., 2016), ZwCl 0008.8+5215 and Abell 1612 (at 5 and 8 GHz, Kierdorf et al., 2017). In combination with low frequency measurements, integrated cluster radio shock spectra spanning over 3 orders of magnitude in frequency have been produced, for example, covering the range from 74/150 MHz to 30 GHz, as is the case for the ‘Sausage’ and the ‘Toothbrush’ radio shocks (van Weeren et al., 2012b; Stroe et al., 2016).

Interferometric observations from 150 MHz to 30 GHz have revealed a possible steepening of the integrated radio shock spectra beyond 2–5 GHz (Stroe et al., 2014b, 2016; Trasatti et al., 2015), which challenges the radio shock formation model involving DSA acceleration at planar shocks. However, studies combining high-frequency single-dish observations with low-frequency interferometric observations (Kierdorf et al., 2017; Loi et al., 2017) do not corroborate this finding (for more details on the caveats of both methods, see below). The mismatch between observations and theory has sparked a discussion as to what is causing the decrement in the flux density of cluster radio shocks at high frequencies (see also Section 6.2). One possibility is that the decrement is not intrinsic to the CR electron distribution at the shock, but is caused by the SZ effect. At 10–30 GHz, the SZ effect is expected to result in a decrement in flux density. Even though the radio shocks are typically located $1-1.5$ Mpc away from the cluster center, authors have argued that the sharp pressure discontinuity from the shock could explain $\sim 20-50\%$ of the decrement (for more typical examples such as the ‘Sausage’, ‘Toothbrush’, or Coma cluster), even up to 100% at the highest frequencies for extreme cases, such as the ‘El Gordo’ or Abell 2256 clusters (depending on the shock geometry, Erler et al., 2015; Basu et al., 2016b).

Various alternatives to the simple shock acceleration model have also been proposed. By contrast to acceleration at time invariant shocks, which results in power-law integrated spectra, curved spectra could be a natural result of spherically-expanding ICM shocks (Kang, 2015b, a). The simple radio shock formation model assumes that the associated shock wave injects thermal electrons. A scenario where the shock predominantly injects non-thermal fossil electrons, pre-accelerated by previous AGN activity, could also reproduce the observed curved radio spectra (Kang & Ryu, 2015). The downstream steepening, as well as the steepening of the integrated spectrum, can be recovered if there is non-uniform magnetic field in the downstream area of the shock (Donnert et al., 2016) or if the electrons, after shock acceleration, are further re-accelerated by turbulence (Fujita et al., 2015; Kang, 2017). Tailored DSA simulations aimed at reproducing the observed parameters of radio shocks with good spectral coverage are now also becoming available (Kang & Ryu, 2015; Kang, 2016b; Kang et al., 2017).

Limitations and caveats:

The combination of low-frequency and high-frequency data to produce wide-frequency spectra can be complicated. Different approaches have been taken in the literature to achieve this: (i) using common baselines at all frequencies trying to ensure recovery of the same spatial scales (Stroe et al., 2014b, 2016) or (ii) combining low-frequency datasets with the best available coverage at short baselines with single dish measurements (Trasatti et al., 2015; Kierdorf et al., 2017; Loi et al., 2017). However, both these approaches come with caveats. Unlike low-frequencies, at high-frequencies, current interferometers do not have very good short-baseline coverage and therefore have trouble recovering extended emission. When using common baseline coverage, the data at high frequency can become too noisy and the spectral behavior of the extended flux is lost. Using low-frequency interferometric data together with single dish measurements has the intrinsic problem of resolving out flux in the interferometric data. In addition, the removal of flux from compact sources that contaminate the single dish measurements is not always straightforward. Current cluster radio shock observations at $10-90$ GHz are limited by their resolution. Interferometric observations in the literature can push down to half arcmin resolution at 10–30 GHz (Stroe et al., 2014b, 2016). By contrast, single dish measurements have the advantage of measuring the total power, but the resolution they can achieve is relatively poor. The largest single dish telescopes (such as the 100-m Effelsberg) can reach $20^{\prime\prime}$ resolution at the highest frequency, but can only achieve one to a few arcmin resolution at 10–30 GHz (Trasatti et al., 2015; Stroe et al., 2016; Kierdorf et al., 2017; Loi et al., 2017).

Future prospects:

The number of radio shock detections above 5 GHz is expected to steadily rise in the following years with observations coming from current instruments, such as the VLA and single-dish telescopes. A number of new facilities are coming or will shortly come online, which will have a significant impact on the study of radio shocks at high-frequencies. Particularly, instruments mounted on large single dish telescopes, such as MUSTANG-2 and NIKA2, will enable SZ studies at high resolution and thus pave the way for joint SZ and X-ray studies of shocks.

Upcoming interferometers will enable the study of the diffuse synchrotron emission from radio shocks at never-before achieved resolution. The low-bands of the ALMA, will provide $5^{\prime\prime}$–$15^{\prime\prime}$ resolution over the 35–50 GHz (Band 1) and $65-90$ GHz (Band 2) range in its most compact configuration. Particularly interesting will be the combination of the ALMA 12-m array with the Atacama Compact Array (a compact configuration of 7-m dishes), which is expected to provide a good compromise in terms of mapping of large scale structures and resolution. In the 2020s, the Square Kilometre Array (SKA)181818https://www.skatelescope.org/ will have observing capabilities up to 10 GHz providing exquisite low surface brightness sensitivity at high resolution (up to 2 milliarcsec at 10 GHz).

6.1.9 Scaling relations

Similar to radio halos, a correlation is found between cluster X-ray luminosity and radio power of cluster radio shocks (Feretti et al., 2012). This correlation likely reflects an underlying correlation between mass an radio power, with $P\propto M^{2.8}$ (de Gasperin et al., 2014). In addition, there is a correlation between the largest linear size (LLS) and distance from the cluster center of the radio shock (van Weeren et al., 2009b; Bonafede et al., 2012; de Gasperin et al., 2014). This is in line with the prediction that in the periphery of clusters the shock surfaces are larger. There is no clear evidence for a correlation between LLS and radio spectral index. Previously, the existence such a correlation has been reported by van Weeren et al. (2009b). However, this LSS–$\alpha$ correlation was produced by the radio phoenices present in the van Weeren et al. (2009b) sample, because radio phoenices generally have smaller LLS and steeper spectra than radio shocks. Nuza et al. (2012); Araya-Melo et al. (2012); Nuza et al. (2017) investigated whether simulations can reproduce the luminosity function, shapes, and LLS distribution of radio shocks. They found reasonable agreement with the properties of radio shocks detectable in the NVSS survey.

6.2 Cosmic ray acceleration modeling at cluster shocks

The acceleration of CR electrons at ICM shocks depends critically on the injection of background electrons into the Fermi-I process and the self-generation of plasma waves that scatter electrons both upstream and downstream of the shock. Background electrons need to be pre-accelerated above the injection momentum, $p_{\rm inj}\sim 130-200~{}p_{\rm th,e}$ ($p_{\rm th,e}=\sqrt{2m_{\rm e}kT_{2}}$), in order to cross the shock transition layer whose width is of the order of gyro radii of thermal protons. Particle-in-cell (PIC) simulations of low Mach number shocks in high beta plasma by Guo et al. (2014a, b) demonstrated that incoming electrons are specularly reflected at the shock ramp by magnetic mirrors and gain energy via multiple cycles of shock drift acceleration, resulting in a suprathermal power-law population of electrons. Necessary scattering waves are self-excited by the firehorse instability. This process is most efficient at quasi-perpendicular shocks where the mean background magnetic field is nearly perpendicular to the shock flow direction. However, the full Fermi-I acceleration that involves scattering of electrons in both upstream and downstream regions of the shock has yet to be studied by PIC simulations.

The merger-shock DSA models for cluster radio shocks have to adopt a set of shock parameters including the pre-shock temperature, $kT_{1}$, sonic Mach number, $\mathcal{M}_{s}$, post-shock magnetic field strength profile, $B_{2}(r)$, and optionally a turbulent acceleration timescale, $\tau_{\rm acc}$, and assume a specific viewing geometry often parameterized with extension angles, $\psi$’s. In addition, in the re-acceleration model, one assumes a fossil electron population with energy spectrum, $N_{\rm fossil}=N_{e}(r)E^{-p}\exp[-(E/E_{c})^{2}]$, in a large volume over $\sim$1 Mpc scale. The power-law slope, $p$, and the energy cutoff, $E_{c}$, can be adjusted to reproduce radio observations.

In particular, DSA models have been successful in reproducing some of observed properties of giant radio shocks such as the thin elongated morphologies, radio flux ($S_{\nu}$) and spectral index ($\alpha_{\nu}$) profiles, and integrated radio spectra ($J_{\nu}$) (Kang et al., 2012). In the case where the radio-inferred Mach number, $\mathcal{M}_{\rm radio}=[(2\alpha_{\rm sh}-3)/(2\alpha_{\rm sh}+1)]^{1/2}\approx 3-4$, is greater than the X-ray-inferred Mach number, $\mathcal{M}_{\rm X}\approx 1.5-3$, the re-acceleration of fossil electrons with a flat spectrum could explain the observed discrepancy (Pinzke et al., 2013; Kang & Ryu, 2016). On the other hand, Zimbardo & Perri (2017, 2018) suggested superdiffusive shock acceleration (SSA) as an alternative explanation for the Mach number discrepancy. SSA is based on superdiffusive transport of energetic particles due to a non-Gaussian (Lévy) random walk. It may lead to CR energy spectra flatter than expected from DSA with normal diffusion.

In addition, the Fermi-II acceleration by post-shock turbulence via transit-time-damping resonance has been invoked to explain the broad downstream steepening of the radio spectrum behind the observed radio shocks (Brunetti & Lazarian, 2007; Kang et al., 2017; Kang, 2017). Thus, the model parameters need to be fine-tuned by comparing theoretical predictions against observations of radio shocks, especially, $S_{\nu}(R)$, $\alpha_{\nu}(R)$, and $J_{\nu}$, at the least.

Figures 26 and 27 demonstrate that such shock re-acceleration models could reproduce the radio observations of the Sausage and the Toothbrush radio shocks. In Figure 26, the shock Mach number is $\mathcal{M}_{\rm s}\approx 3.2$ at $t_{\rm age}=211$ Myr for the Sausage radio shock and $\mathcal{M}_{\rm s}\approx 3.0$ at $t_{\rm age}=144$ Myr for the Toothbrush radio shock. The turbulence acceleration time scale $\tau_{\rm acc}\approx 10^{8}$ yrs is adopted, and the DSA models with (black lines) and without (red lines) post-shock turbulent acceleration are compared. Note that the model shock Mach number $\mathcal{M}_{s}$ is slightly higher than $\mathcal{M}_{\rm X}\approx 2.7$ for the Sausage (Akamatsu et al., 2015), while $\mathcal{M}_{s}$ is much higher than $\mathcal{M}_{\rm X}\approx 1.2-1.5$ for the Toothbrush (van Weeren et al., 2016). Although alternative models with $\mathcal{M}_{s}$ closer to $\mathcal{M}_{\rm X}$ could be adopted to explain both radio and X-ray observations, fossil electrons with hard spectra ($s\approx 1-2\alpha_{\rm obs}$) should be present over a large volume in the ICM. But that seems unrealistic because of the fast cooling time scales of GeV electrons. For the Sausage, the merger scenario itself also remains unclear. Numerical models have trouble to both produce the X-ray properties (like post-shock temperature) and low Mach numbers (Donnert et al., 2017).

As mentioned in Section 6.1.8, the steepening above 2 GHz in the integrated spectrum of the Sausage radio shock has yet to be understood (Stroe et al., 2016). Re-acceleration of electrons by post-shock Alfvénic turbulence (Fujita et al., 2015) and magnetic field amplification behind the shock (Donnert et al., 2016) have been suggested to explain such steepening at high frequencies. In the model shown in Figure 27, the shock sweeps through a finite region of fossil electrons, resulting in continuous softening of $J_{\nu}$ in time (Kang & Ryu, 2016; Kang, 2017). The spectrum at 211 Myr (black solid line) shows the best match with the observed radio data.

6.3 Fossil plasma and CRe re-energization

The study of mildly relativistic AGN fossil plasma throughout clusters is and important topic since, as discussed, old populations of relativistic electrons have been invoked as seed particles for the formation of radio halos and cluster radio shocks. They also retrace past AGN activity and constitute a source of non-thermal pressure in the ICM. Examples of radio phoenices and revived fossil plasma sources are shown in Figure 28.

6.3.1 Radio phoenices and revived fossil plasma

The currently favored scenario is that phoenices trace old radio plasma from past episodes of AGN activity. When a shock compresses this old plasma, the resulting increase in the momentum of the relativistic electrons and the magnetic field strength can produce a source characterized by a steep and curved radio spectrum (Enßlin & Gopal-Krishna, 2001). Simulations also predict that these sources should often have complex morphologies (Enßlin & Brüggen, 2002). It should be noted that so far direct observational evidence for a connection between shocks waves and phoenices is still missing. Therefore, the formation scenario for these revived fossil plasma sources remains somewhat uncertain.

Compared to cluster radio shocks, revived fossil plasma sources and phoenices are on average found at smaller cluster centric distances (Feretti et al., 2012), have smaller sizes ($\lesssim$300–400 kpc, see Figure 15, and have lower radio powers. These revived fossil sources have a range of morphologies, from roundish shapes (e.g., Abell 1664, Giovannini et al., 1999; Govoni et al., 2001b; Kale & Dwarakanath, 2012)) to elongated and filamentary (e.g., Abell 13, Abell 85, Abell 2048, Abell 4038, Abell 2443, Abell 1033191919Not to be confused with the GReET discussed in Section 6.3.3., Abell 1914, Abell 1931, and the Ophiuchus cluster, Slee et al., 1983, 2001; van Weeren et al., 2011d; de Gasperin et al., 2015b; van Weeren et al., 2009a; Werner et al., 2016; Murgia et al., 2010a; Brüggen et al., 2018; Mandal et al., 2018). The elongated and filamentary morphologies, see Figures 30, are the most common (e.g., Slee & Roy, 1998; Slee et al., 2001). Some of these objects are found in cool core clusters such as Abell 85, Abell 1664, and Abell 4038, unlike cluster radio shocks. This indicates that major merger events are not required for their formation.

Radio phoenices and revived fossil sources have integrated spectra that are typically steeper than $-1.5$. In many instances the spectra are curved (Cohen & Clarke, 2011; Slee et al., 2001; van Weeren et al., 2009b), showing high-frequency spectral steepening, see Figure 29 for an example. The spectral index distribution across these sources is irregular without clear common trends (van Weeren et al., 2011d; Cohen & Clarke, 2011; Kale & Dwarakanath, 2012).

Polarized emission from radio phoenices and revived fossil sources has also been detected. The polarization fractions are generally lower than for cluster radio shocks and show larger variations (e.g., Slee et al., 2001). However, it should be remarked that only a few polarization studies have been performed so far of these sources.

6.3.2 Re-acceleration and fossil plasma

As discussed before, DSA shock models proposed for CRe acceleration have found that the acceleration efficiency is often low when electrons are accelerated directly from the thermal pool. This low efficiency is hard to reconcile with the observed brightness and radio spectrum of some cluster radio shocks which suggest a higher acceleration efficiency (e.g., Kang & Ryu, 2011; Vazza & Brüggen, 2014). AGN activity continuously supplies fresh CRs in the ICM creating bright radio galaxies. Due to synchrotron losses, these CRs are visible only for few tens of Myr at Gigahertz frequencies. Although direct observations are prohibitive, a certain amount of CRe with $\gamma\sim 100$ should be present mixed with the ICM Sarazin (1999); Petrosian (2001); Pinzke et al. (2013). Therefore, CR electrons might be re-accelerated from this seed population in the ICM (Enßlin et al., 1998; Markevitch et al., 2005; Kang & Ryu, 2011; Kang et al., 2012), mitigating some of the DSA requirements (see also Section 6.2). An underlying assumption here is that the jets and lobes of radio galaxies are lepton-dominated (e.g., Vazza et al., 2016). Otherwise many CR protons would be re-accelerated, possibly causing problems with the Fermi gamma-ray upper limits (Section 6.1.7).

Instead of fossil ($\gamma\sim 100$) CRe, more energetic CRe from the lobes of a currently active radio galaxy could also re-accelerated (Kang & Ryu, 2016). A few observational pieces of evidence for this scenario were recently reported. In PLCK G287.0+32.9, two large radio shocks have been discovered (Bagchi et al., 2011; Bonafede et al., 2014a). One of the two radio shocks appears to be connected to the lobes of a radio galaxy. However, no optical counterpart for the radio galaxy could be located and the radio spectral index across the source remains difficult to interpret. In the Bullet cluster (1E 0657–55.8), a 930 kpc long radio shock is located opposite to the bullet direction (Shimwell et al., 2015). In this radio shock, a region of 330 kpc has a much higher surface brightness. This might haven be caused by a pre-existing population of CRe of AGN origin. The best example of CRe of AGN origin re-accelerated by a merger shock comes from Abell 3411-3412 (van Weeren et al., 2017b; Johnston-Hollitt, 2017). In this mering system a morphological connection between a radio galaxy and a radio shock is evident. Both polarization and spectral features are in agreement with particle re-acceleration. Furthermore, X-ray data show the presence of a surface brightness discontinuity at the radio shock’s outer edge. However, in the great majority of cases, the presence of a source of CR electrons near the radio shock is missing, leaving unanswered the question: are pre-energized CRe necessary to power all radio shocks? A similar problem is present with radio halos that also require an initial reservoir of mildly energetic CRe to re-energize (Brunetti & Jones, 2014).

With the increase in resolution, sensitivity, and sky coverage of low-frequency telescopes more steep spectrum fossil sources are being discovered. This should shed more light on the connection between diffuse cluster radio sources and AGN fossil plasma in the near future. It has already become clear that galaxy cluster host sources with such steep spectra that they are completely missed at GHz frequencies. Several of these examples have now been uncovered with LOFAR, such as in Abell 1033 (see Section 6.3.3), Abell  1931, and Abell 2034. Recently, the MWA has also found a significant number of new fossil plasma sources and candidates (Duchesne et al., 2017).

6.3.3 GReET

From observations of extended radio sources in the galaxy cluster Abell 1033 (see Figure 31; de Gasperin et al., 2017a), the presence of a possible new mechanism to energize old radio plasma was inferred. In this cluster a WAT source fades into a pair of fairly thin filaments within which the emission starts to brighten and the synchrotron spectrum flattens again. This process of re-energisation is so gentle that it barely balances the radiative losses of cosmic rays, with a particle acceleration time-scale comparable to the radiative loss time-scale of the electrons emitting at $<$100 MHz. This source has been labeled “GReET” (gently re-energized tail).

A proposed physical explanation for the re-energisation mechanism is that Rayleigh-Taylor and Kelvin-Helmholtz instabilities in the tails generate turbulent waves that re-accelerate electrons via second order Fermi mechanisms. The challenge is to understand how the re-acceleration rate is maintained quasi-constant in the tail over a long time-scale. A proposed solution is to assume that turbulence is continuously forced in the tail by the interaction between perturbations in the surrounding medium with the tail itself (de Gasperin et al., 2017a). These perturbations are driven in the medium by the cluster dynamics for a time-scale and on spatial-scales that are larger/comparable to that of the GReET.

If this gentle re-energizing process observed in Abell 1033 is common in tails of radio galaxies in galaxy clusters, then electrons released by radio galaxies in the ICM could live as long as seen in the case of Abell 1033 ($>0.5$ Gyr) and they would be able to accumulate in larger quantities and with higher energies. This could produce a seed population of energetic particles for merger-induced re-acceleration mechanisms, such as turbulence and shocks, that were proposed to explain cluster-scale radio sources. Two other possible GReETs are present in ZwCl 0634.1+4750 (Cuciti et al., 2018) and in Abell 1314 (Wilber et al., 2018a). In both cases a tailed radio galaxy shows an increase in surface brightness along its tail and an unexpected flattening in the spectral index. Because very few examples of GReETs are known, the precise nature of GReETs and their existence as a distinct class of objects remains to be confirmed.

6.3.4 Future prospects

A vast phenomenology of re-energized plasma of AGN origin has recently been emerging, and it attests to the different mechanisms at play: compression (radio phoenices; Enßlin & Gopal-Krishna, 2001), Fermi-I shock re-acceleration (cluster radio shocks; van Weeren et al., 2017b), turbulence ((mini-)halos; ZuHone et al., 2013) or complex plasma interactions (GReETs; de Gasperin et al., 2017a). In most of these cases the re-energization is mild and the radio spectrum is steep, implying that conventional GHz-frequency telescopes overlooked the great majority of these phenomena.

With the current low-frequency telescopes LOFAR, MWA, and the uGMRT, and the future SKA-low, many more revived fossil plasma sources are to be discovered. This should help to better understand the variety of sources present and their spatial distribution in the ICM. Future low-frequency observations should also reveal more connections with cluster radio shocks and possibly with radio halos. These connections can then be studied in more detail. Particularly interesting will be to push these observations towards to lowest frequencies possible ($\lesssim 10-50$ MHz) as current MWA, LOFAR, and uGMRT observations in the 100–300 MHz range probably only probe the tip of the iceberg.

7 Diffuse radio emission outside clusters

Elongated filaments of galaxies span the regions between clusters. Compared to the intracluster medium, the intergalactic medium of galaxy filaments has a significantly lower density ($\lesssim 10^{-4}$ particles cm*-3*) and cooler temperature ($10^{5-7}$ K). About half of the Universe’s baryons reside in this WHIM (Cen & Ostriker, 1999). Galaxy filaments are expected to be surrounded by strong accretion shocks, where the plasma is first shock-heated (Sunyaev & Zeldovich, 1972; Miniati et al., 2000). However, studying the WHIM and associated shocks is notoriously difficult due to a lack of sensitive observational tools.

Owing to their high-Mach numbers ($\mathcal{M}\gtrsim 10$), WHIM accretion shocks should accelerate particles (Miniati et al., 2001b; Ryu et al., 2003; Keshet et al., 2003), similar to what happens in clusters. Radio studies of the WHIM would offer a unique diagnostic tool to determine the strength of the WHIM magnetic field and pinpoint the location of the accretion shocks. However, the detection of the very faint radio emission from these shocks around galaxy filaments is difficult. With larger catalogs of polarized sources, deep continuum images, and simulations, some progress has already been made in deriving the properties of magnetic fields beyond clusters in the cosmic web (e.g., Xu et al., 2006; Brown et al., 2017; Vernstrom et al., 2017; Vazza et al., 2015b, 2017; Gheller et al., 2016). Deep observations with the SKA and other radio telescopes might have the sensitivity to detect radio emission beyond cluster environments (Araya-Melo et al., 2012; Vazza et al., 2015b) and in low mass systems such as groups. A challenge will be to properly classify such emission, since deep observations will also pick up extended low-surface brightness emission associated with (old) radio galaxies.

Despite the expected challenges, some studies have already reported possible extended synchrotron sources from poor clusters and group environments (Nikiel-Wroczyński et al., 2017; Brown & Rudnick, 2009). Also, a candidate for a radio filament was found 5 Mpc away from the cluster Abell 786 (Dewdney et al., 1991; Harris et al., 1993; Giovannini & Feretti, 2000). However, a more recent study suggests that the emission belongs to an old giant radio galaxy (Kale & Dwarakanath, 2012). Radio emission surrounding the ZwCl 2341.1+0000 cluster was reported by Bagchi et al. (2002). Later studies indicate that the radio radio emission is probably associated with a cluster merger event (Giovannini et al., 2010; van Weeren et al., 2009d). Other possible cases of radio emission around clusters associated with accretion shocks (and not from merger events) are the extended radio emission located near MACS J0520.7–1328 (Macario et al., 2014), Abell 3444, (Giovannini et al., 2009), Abell 2255 (Pizzo et al., 2008; Rudnick & Lemmerman, 2009), Abell 1758N–1758S (Botteon et al., 2018a), and MACS J0717.5+3745 (Bonafede et al., 2018). Additionally, candidate radio emission connected to large-scale filaments was recently reported based on the SRT observations (Vacca et al., 2018). However, in all the above mentioned cases the nature of these radio sources still needs to be confirmed, requiring additional observations to shed more light on their origin.

H = giant radio halo; MH = radio-mini-halo; F = revived fossil plasma source; RS = cluster radio shock (relic); U = unclassified; d = double; m = multiple; c = candidate; USS = ultra-steep spectrum; H-MH = “intermediate” or “hybrid” halo–mini-halo

aAbell 1656; bRXC J0338.6+0958; cuncertain; duncertain; $e$Sausage cluster; fRX J0603.3+4214; gZ5247; hZ7160, MS 1455.0+2232; iAbell 1758a ;jAbell 1758b; kZ3146; l1E 0657–5655; mZwCl 1856.8+6616; nCL 1821+643; oACT-CL J0256.5+0006; pCL 14 46$+$26, CL 1447+26; q Z1953; rZ2661; s Z8201; tWHL J235151.0–0.085929; uMACS J0018.5+1626, RXC J0018.5+1626; vEl Gordo

Acknowledgements.

We thank the reviewer for the constructive feedback and suggestions. We thank Francesca Loi, Simona Giacintucci, Annalisa Bonafede, Felipe Andrade-Santos, Marie-Lou Gendron-Marsolais, Frazer Owen, Shea Brown, Lawrence Rudnick, Nobert Werner, Tim Shimwell, and Hans Böhringer for sharing their images and data. We thank Kathrin Böckmann for helping to compile the cluster catalogs. We thank Amanda Wilber, Federica Savini, and Julius Donnert for feedback on the manuscript. We acknowledge the help of Soumyajit Mandal, Huib Intema, and Dawoon Jung for making some of the displayed GMRT images with SPAM (Intema et al., 2009; Intema, 2014).

XMM-Newton images from the Perseus cluster, Abell 2256, and Coma cluster are from Steve Snowden NASA/GSFC. RJvW acknowledges support from the ERC Advanced Investigator programme NewClusters 321271 and the VIDI research programme with project number 639.042.729, which is financed by the Netherlands Organisation for Scientific Research (NWO). FdG is supported by the VENI research programme with project number 639.041.542, which is financed by the Netherlands Organisation for Scientific Research (NWO). HA acknowledges support from the VENI research programme, which is financed by NWO. HK is supported by the National Research Foundation (NRF) of Korea through grants 2016R1A5A1013277 and 2017R1D1A1A09000567. This research made use of APLpy, an open-source plotting package for Python (Robitaille & Bressert, 2012). This research made use of Astropy, a community-developed core Python package for Astronomy (Astropy Collaboration et al., 2013).