The emergence of a complex, large-scale organisation of cosmic matter into the Cosmic Web is a beautiful exemplification of how complexity can be produced by simple initial conditions and simple physical laws. In the epoch of Big Data in astrophysics, connecting the stunning variety of multi-messenger observations to the complex interplay of fundamental physical processes is an open challenge. In this contribution, I discuss a few relevant applications of Information Theory to the task of objectively measuring the complexity of modern numerical simulations of the Universe. When applied to cosmological simulations, complexity analysis makes it possible to measure the total information necessary to model the cosmic web. It also allow us to monitor which physical processes are mostly responsible for the emergence of complex dynamical behaviour across cosmic epochs and environments, and possibly to improve mesh refinement strategies in the future.

Structure formation shocks are believed to be the largest accelerators of cosmic rays in the Universe. However, little is still known about their efficiency in accelerating relativistic electrons and protons as a function of their magnetization properties, i.e. of their magnetic field strength and topology. In this work, we analyzed both uniform and adaptive mesh resolution simulations of large-scale structures with the magnetohydrodynamical grid code Enzo, studying the dependence of shock obliquity with different realistic scenarios of cosmic magnetism. We found that shock obliquities are more often perpendicular than what would be expected from a random three-dimensional distribution of vectors, and that this effect is particularly prominent in the proximity of filaments, due to the action of local shear motions. By coupling these results to recent works from particle-in-cell simulations, we estimated the flux of cosmic-ray protons in galaxy clusters, and showed that in principle the riddle of the missed detection of hadronic gamma-ray emission by the Fermi-LAT can be explained if only quasi-parallel shocks accelerate protons. On the other hand, for most of the cosmic web the acceleration of cosmic-ray electrons is still allowed, due to the abundance of quasi-perpendicular shocks. We discuss quantitative differences between the analyzed models of magnetization of cosmic structures, which become more significant at low cosmic overdensities.

The Square Kilometre Array (SKA) will answer fundamental questions about the origin, evolution, properties, and influence of magnetic fields throughout the Universe. Magnetic fields can illuminate and influence phenomena as diverse as star formation, galactic dynamics, fast radio bursts, active galactic nuclei, large-scale structure, and Dark Matter annihilation. Preparations for the SKA are swiftly continuing worldwide, and the community is making tremendous observational progress in the field of cosmic magnetism using data from a powerful international suite of SKA pathfinder and precursor telescopes. In this contribution, we revisit community plans for magnetism research using the SKA, in the light of these recent rapid developments. We focus in particular on the impact that new radio telescope instrumentation is generating, thus advancing our understanding of key SKA magnetism science areas, as well as the new techniques that are required for processing and interpreting the data. We discuss these recent developments in the context of the ultimate scientific goals for the SKA era.

Measuring the properties of extragalactic magnetic fields through the effect of Faraday rotation provides a means to understand the origin and evolution of cosmic magnetism. Here we use data from the LOFAR Two-Metre Sky Survey (LoTSS) to calculate the Faraday rotation measure (RM) of close pairs of extragalactic radio sources. By considering the RM difference ($\Delta$RM) between physical pairs (e.g. double-lobed radio galaxies) and non-physical pairs (i.e. close projected sources on the sky), we statistically isolate the contribution of extragalactic magnetic fields to $\Delta$RM along the line of sight between non-physical pairs. From our analysis, we find no significant difference between the $\Delta$RM distributions of the physical and non-physical pairs, limiting the excess Faraday rotation contribution to $< 1.9$ rad/m$^2$ ($\sim$$95\%$ confidence). We use this limit with a simple model of an inhomogeneous universe to place an upper limit of 4 nG on the cosmological co-moving magnetic field strength on Mpc scales. We also compare the RM data with a more realistic suite of cosmological MHD simulations, that explore different magnetogenesis scenarios. Both magnetization of the large scale structure by astrophysical processes such as galactic and AGN outflows, and simple primordial scenarios with seed magnetic field strengths $< 0.5$ nG cannot be rejected by the current data; while stronger primordial fields or models with dynamo amplification in filaments are disfavoured.

The origin of radio relics is usually explained via diffusive shock acceleration (DSA) or re-acceleration of electrons at/from merger shocks in galaxy clusters. The case of acceleration is challenged by the low predicted efficiency of low-Mach number merger shocks, unable to explain the power observed in most radio relics. In this Letter we present the discovery of a new giant radio relic around the galaxy cluster Abell 2249 ($z=0.0838$) using LOFAR. It is special since it has the lowest surface brightness of all known radio relics. We study its radio and X-ray properties combinig LOFAR data with uGMRT, JVLA and XMM. This object has a total power of $L_{1.4\rm GHz}=4.1\pm 0.8 \times 10^{23}$ W Hz$^{-1}$ and integrated spectral index $\alpha = 1.15\pm 0.23$. We infer for this radio relic a lower bound on the magnetisation of $B\geq 0.4\, \mu$G, a shock Mach number of $\mathcal{M}\approx 3.79$, and a low acceleration efficiency consistent with DSA. This result suggests that a missing population of relics may become visible thanks to the unprecedented sensitivity of the new generation of radio telescopes.

Observations of large-scale radio emissions prove the existence of shock accelerated cosmic-ray electrons in galaxy clusters, while the lack of detected $\gamma$-rays limits the acceleration of cosmic-ray protons in galaxy clusters. This challenges our understanding of how diffusive shock acceleration works. In this work, we couple the most updated recipes for shock acceleration in the intracluster medium to state-of-the-art magneto-hydrodynamical simulations of massive galaxy clusters. Furthermore, we use passive tracer particles to follow the evolution of accelerated cosmic-rays. We show that when the interplay between magnetic field topology and the feedback from accelerated cosmic rays is taken into account, the latest developments of particle acceleration theory give results which are compatible with observational constraints.

We used magneto-hydrodynamical cosmological simulations to investigate the cross-correlation between different observables (i.e. X-ray emission, Sunyaev-Zeldovich signal at 21 cm, HI temperature decrement, diffuse synchrotron emission and Faraday Rotation) as a probe of the diffuse matter distribution in the cosmic web. We adopt an uniform and simplistic approach to produce synthetic observations at various wavelengths, and we compare the detection chances of different combinations of observables correlated with each other and with the underlying galaxy distribution in the volume. With presently available surveys of galaxies and existing instruments, the best chances to detect the diffuse gas in the cosmic web outside of halos is by cross-correlating the distribution of galaxies with Sunyaev-Zeldovich observations. We also find that the cross-correlation between the galaxy network and the radio emission or the Faraday Rotation can already be used to limit the amplitude of extragalactic magnetic fields, well outside of the cluster volume usually explored by existing radio observations, and to probe the origin of cosmic magnetism with the future generation of radio surveys.

The degree of turbulent pressure support by residual gas motions in galaxy clusters is not well known. Mass modelling of combined X-ray and Sunyaev Zel'dovich observations provides an estimate of turbulent pressure support in the outer regions of several galaxy clusters. Here, we test two different filtering techniques to disentangle bulk from turbulent motions in non-radiative high-resolution cosmological simulations of galaxy clusters using the cosmological hydro code ENZO. We find that the radial behavior of the ratio of non-thermal pressure to total gas pressure as a function of cluster-centric distance can be described by a simple polynomial function. The typical non-thermal pressure support in the centre of clusters is $\sim$5%, increasing to $\sim$15% in the outskirts, in line with the pressure excess found in recent X-ray observations. While the complex dynamics of the ICM makes it impossible to reconstruct a simple correlation between turbulent motions and hydrostatic bias, we find that a relation between them can be established using the median properties of a sample of objects. Moreover, we estimate the contribution of radial accelerations to the non-thermal pressure support and conclude that it decreases moving outwards from 40% (in the core) to 15% (in the cluster's outskirts). Adding this contribution to one provided by turbulence, we show that it might account for the entire observed hydrostatic bias in the innermost regions of the clusters, and for less than 80% of it at $r > 0.8 r_{200, m}$.

Fast radio bursts remain one of the most enigmatic astrophysical sources. Observations have significantly progressed over the last few years, thanks to the capabilities of new radio telescopes and the refurbishment of existing ones. Here we describe the upgrade of the Northern Cross radio telescope, operating in the 400-416~MHz frequency band, with the ultimate goal of turning the array into a dedicated instrument to survey the sky for fast radio bursts. We present test observations of the pulsar B0329+54 to characterize the system performance and forecast detectability. Observations with the system currently in place are still limited by modest sky coverage ($\sim 9.4$~deg$^2$) and biased by smearing of high dispersion measure events within each frequency channels. In its final, upgraded configuration, however, the telescope will be able to carry out unbiased fast radio burst surveys over a $\sim 350$~deg$^2$ instantaneous field of view up to $z \sim 5$, with a (nearly constant) $\sim 760 \, (\tau/{\rm ms})^{-0.5}$~mJy rms sensitivity.

Radio relics are diffuse extended synchrotron sources that originate from shock fronts induced by galaxy cluster mergers. The actual particle acceleration mechanism at the shock fronts is still under debate. The galaxy cluster 1RXS J0603.3+4214 hosts one of the most intriguing examples of radio relics, known as the Toothbrush. We present new wideband radio continuum observations made with uGMRT and VLA. Our new observations, in combination with previously published data, allowed us to carry out a detailed high spatial resolution spectral and curvature analysis of the known diffuse radio emission sources, over a broad range of frequencies. The integrated spectrum of the Toothbrush follows closely a power law over close to 2 decades in frequency, with a spectral index of $-1.16\pm0.02$. We do not find any evidence of spectral steepening below 8 GHz. The subregions of the Toothbrush show an identical spectral slopes, suggesting that observed spectral index is rather set by the distribution of Mach numbers which may have a similar shape at different parts of the shock front. Indeed, numerical simulations show an intriguing similar spectral index, indicating that the radio spectrum is dominated by the average over the inhomogeneities within the shock, with most of the emission coming from the tail of the Mach number distribution. In contrast to the Toothbrush, the spectrum of the fainter relics show a high frequency steepening. The integrated spectrum of the halo follows a power law from 150 MHz to 3 GHz with a spectral index of $-1.16\pm0.04$. We do not find any evidence for spectral curvature, not even in subareas of the halo. This suggest a homogeneous acceleration throughout the cluster volume. Between the brush region of the Toothbrush and the halo, the color-color analysis revealed emission that was consistent with an overlap between the two different spectral regions.

Radio relics are vast synchrotron sources that sit on the outskirts of merging galaxy clusters. In this work we model their formation using a Press-Schechter formalism to simulate merger rates, analytical models for the intracluster medium and the shock dynamics, as well as a simple model for the cosmic-ray electrons at the merger shocks. We show that the statistical properties of the population of radio relics are strongly dependent on key physical parameters, such as the acceleration efficiency, the magnetic field strength at the relic, the geometry of the relic and the duration of the electron acceleration at merger shocks. It turns out that the flux distribution as well as the power-mass relation can constrain key parameters of the intracluster medium. With the advent of new large-area radio surveys, statistical analyses of radio relics will complement what we have learned from observations of individual objects.

Radio observations at low frequencies with the Low Frequency Array (LOFAR) start discovering gigantic radio bridges connecting pairs of massive galaxy clusters. These observations probe unexplored mechanisms of in situ particle acceleration that operate on volumes of several Mpc$^3$. Numerical simulations suggest that such bridges are dynamically complex and that weak shocks and super-Alfv\'{e}nic turbulence can be driven across the entire volume of these regions. In this Letter we explore, for the first time, the role of second order Fermi mechanisms for the reacceleration of relativistic electrons interacting with turbulence in these peculiar regions. We assume the turbulent energy flux measured in simulations and adopt a scenario in which relativistic particles scatter with magnetic field lines diffusing in super-Alfv\'enic turbulence and magnetic fields are amplified by the same turbulence. We show that steep spectrum and volume filling synchrotron emission can be generated in the entire intra-cluster bridge region thus providing a natural explanation for radio bridges. Consequently, radio observations have the potential to probe the dissipation of energy on scales larger than galaxy clusters and second order Fermi mechanisms operating in physical regimes that are still poorly explored. This has a potential impact on several branches of astrophysics and cosmology.

The growth of large-scale cosmic structure is a beautiful exemplification of how complexity can emerge in our Universe, starting from simple initial conditions and simple physical laws. Using {\enzo} cosmological numerical simulations, I applied tools from Information Theory (namely, "statistical complexity") to quantify the amount of complexity in the simulated cosmic volume, as a function of cosmic epoch and environment. This analysis can quantify how much difficult to predict, at least in a statistical sense, is the evolution of the thermal, kinetic and magnetic energy of the dominant component of ordinary matter in the Universe (the intragalactic medium plasma). The most complex environment in the simulated cosmic web is generally found to be the periphery of large-scale structures (e.g. galaxy clusters and filaments), where the complexity is on average $\sim 10-10^2$ times larger than in more rarefied regions, even if the latter dominate the volume-integrated complexity of the simulated Universe. If the energy evolution of gas in the cosmic web is measured on a $\approx 100 $ $\rm kpc/h$ resolution and over a $\approx 200$ $\rm Myr$ timescale, its total complexity is the range of $\sim 10^{16}-10^{17} \rm ~bits$, with little dependence on the assumed gas physics, cosmology or cosmic variance.

The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy and spacetime.

We reconstruct the 3D structure of magnetic fields, which were seeded by density perturbations during the radiation dominated epoch of the Universe and later on were evolved by structure formation. To achieve this goal, we rely on three dimensional initial density fields inferred from the 2M++ galaxy compilation via the Bayesian $\texttt{BORG}$ algorithm. Using those, we estimate the magnetogenesis by the so called Harrison mechanism. This effect produced magnetic fields exploiting the different photon drag on electrons and ions in vortical motions, which are exited due to second order perturbation effects in the Early Universe. Subsequently we study the evolution of these seed fields through the non-linear cosmic structure formation by virtue of a MHD simulation to obtain a 3D estimate for the structure of this primordial magnetic field component today. At recombination we obtain a reliable lower limit on the large scale magnetic field strength around $10^{-23} \mathrm{G}$, with a power spectrum peaking at about $ 2\, \mathrm{Mpc}^{-1}h$ in comoving scales. At present we expect this evolved primordial field to have strengthts above $\approx 10^{-27}\, \mathrm{G}$ and $\approx 10^{-29}\, \mathrm{G}$ in clusters of galaxies and voids, respectively. We also calculate the corresponding Faraday rotation measure map and show the magnetic field morphology and strength for specific objects of the Local Universe.

Radio emission in the form of giant radio relics is observed at the periphery of galaxy clusters. This non-thermal emission is an important tracer for cosmic-ray electrons and intracluster magnetic fields. One striking observational feature of these objects is their high degree of polarisation which provides information on the magnetic fields at the relics' positions. In this contribution, we test if state-of-the-art high resolution cosmological simulations are able to reproduce the polarisation features of radio relics. Therefore, we present a new analysis of high-resolution cosmological simulations to study the polarisation properties of radio relics in detail. In order to compare our results with current and future radio observations, we create mock radio observations of the diffuse polarised emission from a massive galaxy clusters using six different projections, for different observing frequencies and for different telescopes. Our simulations suggest that, due to the effect of Faraday rotation, it is extremely difficult to relate the morphology of the polarised emission for observing frequencies below $1.4 \ \mathrm{GHz}$ to the real magnetic field structure in relics. We can reproduce the observed degree of polarisation and also several small-scale structures observed in real radio relics, but further work would be needed to reproduce some large-scale spectacular features as observed in real radio relics, such as the "Sausage" and the "Toothbrush" relics.

With ENZO simulations run on the J\"ulich supercomputers, we have investigated the evolution of magnetic fields in the largest cosmic structures (namely galaxy clusters and filaments connecting them) with unprecedented dynamical range. These simulations revealed the full development of the small-scale dynamo in Eulerian cosmological magneto-hydrodynamical simulations. The turbulent motions developed during the formation of clusters are energetic enough to foster the growth of magnetic fields by several orders of magnitude, starting from weak magnetic fields up strengths of $\sim \rm \mu G$ as observed. Furthermore, shock waves are launched during cluster formation and they are able to accelerate cosmic-ray electrons, that emit in the radio wavelengths. Radio observations of this emission provide information on the local magnetic field strength. We have incorporated, for the first time, the cooling of cosmic-ray electrons when modelling this emission. In this contribution, we present our advances in modelling these physical processes. Here, we mostly focus on the most interesting object in our sample of galaxy clusters, which shows the complexity of magnetic fields and the potential of existing and future multi-wavelengths observations in the study of the weakly collisional plasma on $\sim$ Megaparsecs scales.

Radio relics are sites of electron (re)acceleration in merging galaxy clusters but the mechanism of acceleration and the topology of the magnetic field in and near relics are yet to be understood. We are carrying out an observational campaign on double relic galaxy clusters starting with RXC J1314.4-2515. With $Jansky Very Large Array$ multi-configuration observations in the frequency range 1-4 GHz, we perform both spectral and polarization analyses, using the Rotation Measure synthesis technique. We use archival $XMM-Newton$ observations to constrain the properties of the shocked region. We discover a possible connection between the activity of a radio galaxy and the emission of the eastern radio relic. In the northern elongated arc of the western radio relic, we detect polarized emission with an average polarization fraction of $31 \ \%$ at 3 GHz and we derive the Mach number of the underlying X-ray shock. Our observations reveal low levels of fractional polarization and Faraday-complex structures in the southern region of the relic, which point to the presence of thermal gas and filamentary magnetic field morphology inside the radio emitting volume. We measured largely different Rotation Measure dispersion from the two relics. Finally, we use cosmological magneto-hydrodynamical simulations to constrain the magnetic field, viewing angle, and to derive the acceleration efficiency of the shock. We find that the polarization properties of RXC J1314.4-2515 are consistent with a radio relic observed at $70^{\circ}$ with respect to the line of sight and that efficient re-acceleration of fossil electrons has taken place.

The majority of the ordinary matter in the local Universe has been heated by strong structure formation shocks and resides in a largely unexplored hot, diffuse, X-ray emitting plasma that permeates the halos of galaxies, galaxy groups and clusters, and the cosmic web. We propose a next-generation "Cosmic Web Explorer" that will permit a complete and exhaustive understanding of these unseen baryons. This will be the first mission capable to reach the accretion shocks located several times farther than the virial radii of galaxy clusters, and reveal the out-of-equilibrium parts of the intra-cluster medium which are live witnesses to the physics of cosmic accretion. It will also enable a view of the thermal, kinematic, and chemical composition of the circum-galactic medium in galaxies with masses similar to the Milky Way, at the same level of detail that $Athena$ will unravel for the virialized regions of massive galaxy clusters, delivering a transformative understanding of the evolution of those galaxies in which most of the stars and metals in the Universe were formed. Finally, the proposed X-ray satellite will connect the dots of the large-scale structure by mapping, at high spectral resolution, as much as 100% of the diffuse gas hotter than $10^6$ K that fills the filaments of the cosmic web at low redshifts, down to an over-density of 1, both in emission and in absorption against the ubiquitous cosmic X-ray background, surveying up to half of the sky over 5 years in orbit. This requires a large effective area (~10 m$^2$ at 1 keV) over a large field of view (~1 deg$^2$), a megapixel cryogenic microcalorimeter array providing integral field spectroscopy with a resolving power $E/\Delta E$ = 2000 at 0.6 keV and a spatial resolution of 5" in the soft X-ray band, and a Low-Earth orbit that ensures a stable instrumental background and high sensitivity to faint, extended emission.

We investigate the possibility of measuring intergalactic magnetic fields using the dispersion measures and rotation measures of fast radio bursts. With Bayesian methods, we produce probability density functions for values of these measures. We distinguish between contributions from the intergalactic medium, the host galaxy and the local environment of the progenitor. To this end, we use constrained, magnetohydrodynamic simulations of the local Universe to compute lines-of-sight integrals from the position of the Milky Way. In particular, we differentiate between predominantly astrophysical and primordial origins of magnetic fields in the intergalactic medium. We test different possible types of host galaxies and probe different distribution functions of fast radio burst progenitor locations inside the host galaxy. Under the assumption that fast radio bursts are produced by magnetars, we use analytic predictions to account for the contribution of the local environment. We find that less than 100 fast radio bursts from magnetars in stellar-wind environments hosted by starburst dwarf galaxies at redshift $z \gtrsim 0.5$ suffice to discriminate between predominantly primordial and astrophysical origins of intergalactic magnetic fields. However, this requires the contribution of the Milky Way to be removed with a precision of $\approx 1 \rm~rad~m^{-2}$. We show the potential existence of a subset of fast radio bursts whose rotation measure carry information on the strength of the intergalactic magnetic field and its origins.

Galaxy clusters are the most massive gravitationally bound structures in the Universe. They grow by accreting smaller structures in a merging process that produces shocks and turbulence in the intra-cluster gas. We observed a ridge of radio emission connecting the merging galaxy clusters Abell 0399 and Abell 0401 with the Low Frequency Array (LOFAR) at 140 MHz. This emission requires a population of relativistic electrons and a magnetic field located in a filament between the two galaxy clusters. We performed simulations to show that a volume-filling distribution of weak shocks may re-accelerate a pre-existing population of relativistic particles, producing emission at radio wavelengths that illuminates the magnetic ridge.

Detecting the thermal and non-thermal emission from the shocked cosmic gas surrounding large-scale structures represents a challenge for observations, as well as a unique window into the physics of the warm-hot intergalactic medium. In this work, we present synthetic radio and X-ray surveys of large cosmological simulations in order to assess the chances of jointly detecting the cosmic web in both frequency ranges. We then propose best observing strategies tailored for existing (LOFAR, MWA and XMM) or future instruments (SKA-LOW and SKA-MID, ATHENA and eROSITA). We find that the most promising targets are the extreme peripheries of galaxy clusters in an early merging stage, where the merger causes the fast compression of warm-hot gas onto the virial region. By taking advantage of a detection in the radio band, future deep X-ray observations will probe this gas in emission, and help us to study plasma conditions in the dynamic warm-hot intergalactic medium with unprecedented detail.

We investigate the evolution of magnetic fields in galaxy clusters starting from constant primordial fields using highly resolved ($\approx \rm 4 ~kpc$) cosmological MHD simulations. The magnetic fields in our sample exhibit amplification via a small-scale dynamo and compression during structure formation. In particular, we study how the spectral properties of magnetic fields are affected by mergers, and we relate the measured magnetic energy spectra to the dynamical evolution of the intracluster medium. The magnetic energy grows by a factor of $\sim$ 40-50 in a time-span of $\sim 9$ Gyr and equipartition between kinetic and magnetic energy occurs on a range of scales ($< 160 \rm ~kpc$ at all epochs) depending on the turbulence state of the system. We also find that, in general, the outer scale of the magnetic field and the MHD scale are not simply correlated in time. The effect of major mergers is to shift the peak magnetic spectra to it smaller scales, whereas the magnetic amplification only starts after $\lesssim$ 1 Gyr. In contrast, continuous minor mergers promote the steady growth of the magnetic field. We discuss the implications of these findings in the interpretation of future radio observations of galaxy clusters.

In this paper, we exploit a large suite of {\enzo} cosmological magneto-hydrodynamical simulations adopting uniform mesh resolution, to investigate the properties of cosmic filaments under different baryonic physics and magnetogenesis scenarios. We exploit a isovolume based algorithm to identify filaments and determine their attributes from the continuous distribution of gas mass density in the simulated volumes. The global (e.g. mass, size, mean temperature and magnetic field strength, enclosed baryon fraction) and internal (e.g. density, temperature, velocity and magnetic field profiles) properties of filaments in our volume are calculated across almost four orders of magnitude in mass. The inclusion of variations in non-gravitational physical processes (radiative cooling, star formation, feedback from star forming regions and active galactic nuclei) as well as in the seeding scenarios for magnetic fields (early magnetisation by primordial process vs later seeding by galaxies) allows us to study both the large-scale thermodynamics and the magnetic properties of the Warm-Hot Intergalactic Medium (WHIM) with an unprecedented detail. We show how the impact of non-gravitational physics on the global thermodynamical properties of filaments is modest, with the exception of the densest gas environment surrounding galaxies in filaments. Conversely, the magnetic properties of the WHIM in filament are found to dramatically vary as different seeding scenarios are considered. We study the correlation between the properties of galaxy-sized halos and their host filaments, as well as between the halos and the local WHIM in which they lie. Significant general statistical trends are reported.

In recent years, the outskirts of galaxy clusters have emerged as one of the new frontiers and unique laboratories for studying the growth of large scale structure in the universe. Modern cosmological hydrodynamical simulations make firm and testable predictions of the thermodynamic and chemical evolution of the X-ray emitting intracluster medium. However, recent X-ray and Sunyaev-Zeldovich effect observations have revealed enigmatic disagreements with theoretical predictions, which have motivated deeper investigations of a plethora of astrophysical processes operating in the virialization region in the cluster outskirts. Much of the physics of cluster outskirts is fundamentally different from that of cluster cores, which has been the main focus of X-ray cluster science over the past several decades. A next-generation X-ray telescope, equipped with sub-arcsecond spatial resolution over a large field of view along with a low and stable instrumental background, is required in order to reveal the full story of the growth of galaxy clusters and the cosmic web and their applications for cosmology.

We simulate the propagation of cosmic rays at ultra-high energies, $\gtrsim 10^{18}$ eV, in models of extragalactic magnetic fields in constrained simulations of the local Universe. We investigate the impact of different magneto-genesis scenarios, both, primordial and astrophysical, on the propagation of cosmic rays. Our study shows that different scenarios of magneto-genesis do not have a large impact on the anisotropy measurements. The distribution of nearby sources causes anisotropy at very high energies, independent of the magnetic field model. We compare our results to the dipole signal measured by the Pierre Auger Observatory. All our models could reproduce the observed dipole amplitude with a pure iron injection composition. This is due to clustering of secondary nuclei in direction of nearby sources of heavy nuclei. A light injection composition is disfavoured by the non-observation of anisotropy at energies of 4 - 8 EeV.

Galaxy clusters grow by gas accretion, mostly from mergers of substructures, which release powerful shock waves into cosmic plasmas and convert a fraction of kinetic energy into thermal energy, amplification of magnetic fields and into the acceleration of energetic particles. The modeling of the radio signature of cosmic shocks, combined with the lack of detected gamma-rays from cosmic ray (CR) protons, poses challenges to our understanding of how cosmic rays get accelerated and stored in the intracluster medium. Here we review the injection of CRs by cosmic shocks of different strengths, combining the detailed "microscopic" view of collisionless processes governing the creation of non-thermal distributions of electrons and protons in cluster shocks (based on analytic theory and particle-in-cell simulations), with the "macroscopic" view of the large-scale distribution of cosmic rays, suggested by modern cosmological simulations. We discuss time dependent non-linear kinetic models of particle acceleration by multiple internal shocks with large scale compressible motions of plasma. The models produce soft CR spectra containing a noticeable energy density in the super-thermal protons of energies below a few GeV which are difficult to constrain by Fermi observations. We consider the effect of plasma composition on CR injection and super-thermal particle population in the hot intracluster matter which can be constrained by fine high resolution X-ray spectroscopy of Fe ions.

The detection of filaments in the cosmic web will be crucial to distinguish between the possible magnetogenesis scenarios and future large polarization surveys will be able to shed light on their magnetization level. In this work, we use numerical simulations of galaxy clusters to investigate their possible detection. We compute the Faraday Rotation signal in intracluster filaments and compare it to its surrounding environment. We find that the expected big improvement in sensitivity with the SKA-MID will in principle allow the detection of a large fraction of filaments surrounding galaxy clusters. However, the contamination of the intrinsic Faraday Rotation of background polarized sources will represent a big limiter to the number of objects that can be significantly detected. We discuss possible strategies to minimize this effect and increase the chances of detection of the cosmic web with the large statistics expected from future surveys.

Centrally located diffuse radio emission has been observed in both merging and non-merging galaxy clusters. Depending on their morphology and size, we distinguish between giant radio haloes, which occur predominantly in merging clusters, and mini haloes, which are found in non-merging, cool-core clusters. Low-frequency sensitive observations are required to assess whether the emission discovered in these few cases is common in galaxy clusters or not. With this aim, we carried out a campaign of observations with the LOw Frequency ARay (LOFAR) in the frequency range 120 - 168 MHz of nine massive clusters selected from the \textit{Planck} SZ catalogue, which had no sign of major mergers. In this paper, we discuss the results of the observations that have led to the largest cluster sample studied within the LOFAR Two-metre Sky Survey, and we present Chandra X-ray data used to investigate the dynamical state of the clusters, verifying that the clusters are currently not undergoing major mergers, and to search for traces of minor or off-axis mergers. We discover large-scale steep-spectrum emission around mini haloes in the cool-core clusters PSZ1G139 and RXJ1720, which is not observed around the mini halo in the non-cool-core cluster A1413. We also discover a new 570 kpc-halo in the non-cool-core cluster RXCJ0142. We derived new upper limits to the radio power for clusters in which no diffuse radio emission was found, and we discuss the implication of our results to constrain the cosmic-ray energy budget in the ICM. We conclude that radio emission in non-merging massive clusters is not common at the sensitivity level reached by our observations and that no clear connection with the cluster dynamical state is observed. Our results might indicate that the sloshing of a dense cool core could trigger particle acceleration on larger scales and generate steep-spectrum radio emission.

Cosmological simulations predict that an intergalactic magnetic field (IGMF) pervades the large scale structure (LSS) of the Universe. Measuring the IGMF is important to determine its origin (i.e. primordial or otherwise). Using data from the LOFAR Two Metre Sky Survey (LoTSS), we present the Faraday rotation measure (RM) and depolarisation properties of the giant radio galaxy J1235+5317, at a redshift of $z = 0.34$ and 3.38 Mpc in size. We find a mean RM difference between the lobes of $2.5\pm0.1$ rad/m$^2$ , in addition to small scale RM variations of ~0.1 rad/m$^2$ . From a catalogue of LSS filaments based on optical spectroscopic observations in the local universe, we find an excess of filaments intersecting the line of sight to only one of the lobes. Associating the entire RM difference to these LSS filaments leads to a gas density-weighted IGMF strength of ~0.3 {\mu}G. However, direct comparison with cosmological simulations of the RM contribution from LSS filaments gives a low probability (~5%) for an RM contribution as large as 2.5 rad/m$^2$ , for the case of IGMF strengths of 10 to 50 nG. It is likely that variations in the RM from the Milky Way (on 11' scales) contribute significantly to the mean RM difference, and a denser RM grid is required to better constrain this contribution. In general, this work demonstrates the potential of the LOFAR telescope to probe the weak signature of the IGMF. Future studies, with thousands of sources with high accuracy RMs from LoTSS, will enable more stringent constraints on the nature of the IGMF.

The hot plasma in galaxy clusters is expected to be heated to high temperatures through shocks and adiabatic compression. The thermodynamical properties of the gas encode information on the processes leading to the thermalization of the gas in the cluster's potential well as well as non-gravitational processes such as gas cooling, AGN feedback and kinetic energy. In this work we present the radial profiles of the thermodynamic properties of the intracluster medium (ICM) out to the virial radius for a sample of 12 galaxy clusters selected from the Planck all-sky survey. We determine the universal profiles of gas density, temperature, pressure, and entropy over more than two decades in radius. We exploit jointly X-ray information from XMM and Sunyaev-Zel'dovich constraints from Planck to recover thermodynamic properties out to 2 R500. We provide average functional forms for the radial dependence of the main quantities and quantify the slope and intrinsic scatter of the population as a function of radius. We find that gas density and pressure profiles steepen steadily with radius, in excellent agreement with previous observational results. Entropy profiles beyond R500 closely follow the predictions for the gravitational collapse of structures. The scatter in all thermodynamical quantities reaches a minimum in the range [0.2-0.8] R500 and increases outwards. Somewhat surprisingly, we find that pressure is substantially more scattered than temperature and density. Our results indicate that once accreting substructures are properly excised, the properties of the ICM beyond the cooling region R > 0.3 R500) follow remarkably well the predictions of simple gravitational collapse and require little non-gravitational corrections.

We review the present theoretical and numerical understanding of magnetic field amplification in cosmic large-scale structure, on length scales of galaxy clusters and beyond. Structure formation drives compression and turbulence, which amplify tiny magnetic seed fields to the microGauss values that are observed in the intracluster medium. This process is intimately connected to the properties of turbulence and the microphysics of the intra-cluster medium. Additional roles are played by merger induced shocks that sweep through the intra-cluster medium and motions induced by sloshing cool cores. The accurate simulation of magnetic field amplification in clusters still poses a serious challenge for simulations of cosmological structure formation. We review the current literature on cosmological simulations that include magnetic fields and outline theoretical as well as numerical challenges.

We introduce a new sample of galaxy clusters obtained from a cosmological simulation covering an unprecedented dynamical range. All the clusters in our sample show a clear signature of small-scale dynamo amplification. We show that it is possible to use dynamo theory for studying the magnetic spectrum in the intracluster medium. We study if the intrinsic variations on the spectra depend on the dynamical history of each cluster or on some host cluster properties.

Galaxy clusters are the endpoints of structure formation and are continuously growing through the merging and accretion of smaller structures. Numerical simulations predict that a fraction of their energy content is not yet thermalized, mainly in the form of kinetic motions (turbulence, bulk motions). Measuring the level of non-thermal pressure support is necessary to understand the processes leading to the virialization of the gas within the potential well of the main halo and to calibrate the biases in hydrostatic mass estimates. We present high-quality measurements of hydrostatic masses and intracluster gas fraction out to the virial radius for a sample of 12 nearby clusters with available XMM-Newton and Planck data. We compare our hydrostatic gas fractions with the expected universal gas fraction to constrain the level of non-thermal pressure support. We find that hydrostatic masses require little correction and infer a median non-thermal pressure fraction of $\sim6\%$ and $\sim10\%$ at $R_{500}$ and $R_{200}$, respectively. Our values are lower than the expectations of hydrodynamical simulations, possibly implying a faster thermalization of the gas. If instead we use the mass calibration adopted by the Planck team, we find that the gas fraction of massive local systems implies a mass bias $1-b=0.85\pm0.05$ for SZ-derived masses, with some evidence for a mass-dependent bias. Conversely, the high bias required to match Planck CMB and cluster count cosmology is excluded by the data at high significance, unless the most massive halos are missing a substantial fraction of their baryons.

In this paper we introduce a reliable, fully automated and fast algorithm to detect extended extragalactic radio sources (cluster of galaxies, filaments) in existing and forthcoming surveys (like LOFAR and SKA). The proposed solution is based on the adoption of a Deep Learning approach, more specifically a Convolutional Neural Network, that proved to perform outstandingly in the processing, recognition and classification of images. The challenge, in the case of radio interferometric data, is the presence of noise and the lack of a sufficiently large number of labeled images for the training. We have specifically addressed these problems and the resulting software, COSMODEEP proved to be an accurate, efficient and effective solution for detecting very faint sources in the simulated radio images. We present the comparison with standard source finding techniques, and discuss advantages and limitations of our new approach.

Due to their late formation in cosmic history, clusters of galaxies are not fully in hydrostatic equilibrium and the gravitational pull of their mass at a given radius is expected not to be entirely balanced by the thermal gas pressure. Turbulence may supply additional pressure, and recent (X-ray and SZ) hydrostatic mass reconstructions claim a pressure support of $\sim 5-15\%$ of the total pressure at $R_{\rm 200}$. In this work we show that, after carefully disentangling bulk from small-scale turbulent motions in high-resolution simulations of galaxy clusters, we can constrain which fraction of the gas kinetic energy effectively provides pressure support in the cluster's gravitational potential. While the ubiquitous presence of radial inflows in the cluster can lead to significant bias in the estimate of the non-thermal pressure support, we report that only a part of this energy effectively acts as a source of pressure, providing a support of the order of $\sim 10\%$ of the total pressure at $R_{\rm 200}$.

The Advanced Telescope for High Energy Astrophysics (Athena) is the X-ray observatory large mission selected by the European Space Agency (ESA), within its Cosmic Vision 2015-2025 programme, to address the "Hot and Energetic Universe" scientific theme (Nandra et al. 2013), and it is provisionally due for launch in the early 2030s. The Square Kilometer Array (SKA) is the next generation radio observatory and consists of two telescopes, one comprised of dishes operating at mid frequencies (SKA1-MID) and located in South Africa, and the other comprised of Log-Periodic antennas operating at low radio frequencies (SKA1-LOW), which will be located in Australia (Braun et al. 2017). The scientific commissioning of the radio telescope is planned to begin in 2021-2022. The SKA-Athena Synergy Team (SAST) has been tasked to single out the potential scientific synergies between Athena and SKA. The astrophysical community was involved in this exercise primarily through a dedicated SKA-Athena Synergy Workshop, which took place on April 24-25, 2017 at SKAO, Jodrell Bank, Manchester. The final result of the synergy exercise, this White Paper, describes in detail a number of scientific opportunities that will be opened up by the combination of Athena and SKA, these include: 1. the Cosmic Dawn; 2. the Evolution of black holes and galaxies; 3. Active galaxy feedback in galaxy clusters; 4. Non-thermal phenomena in galaxy clusters; 5. Detecting the cosmic web; 6. Black-hole accretion physics and astrophysical transients; 7. Galactic astronomy: stars, planets, pulsars and supernovae.

The joint analysis of the Dispersion and Faraday Rotation Measure from distant, polarised Fast Radio Bursts may be used to put constraints on the origin and distribution of extragalactic magnetic fields on cosmological scales. While the combination of Dispersion and Faraday Rotation Measure can in principle give the average magnetic fields along the line-of-sight, in practice this method must be used with care because it strongly depends on the assumed magnetisation model on large cosmological scales. Our simulations show that the observation of Rotation Measures with $\geq 1-10 ~\rm rad/m^2$ in $\sim 10^2$ Fast Radio Bursts will be able to discriminate between extreme scenarios for the origin of cosmic magnetic fields, independent of the exact distribution of sources with redshift. This represent a strong case for incoming (e.g. ALERT, CHIME) and future (e.g. with the Square Kilometer Array) radio polarisation surveys of the sky.

We present results from LOFAR and GMRT observations of the galaxy cluster MACS$\,$J0717.5$+$3745. The cluster is undergoing a violent merger involving at least four sub-clusters, and it is known to host a radio halo. LOFAR observations reveal new sources of radio emission in the Intra-Cluster Medium: (i) a radio bridge that connects the cluster to a head-tail radio galaxy located along a filament of galaxies falling into the main cluster, (ii) a 1.9 Mpc radio arc, that is located North West of the main mass component, (iii) radio emission along the X-ray bar, that traces the gas in the X-rays South West of the cluster centre. We use deep GMRT observations at 608 MHz to constrain the spectral indices of these new radio sources, and of the emission that was already studied in the literature at higher frequency. We find that the spectrum of the radio halo and of the relic at LOFAR frequency follows the same power law as observed at higher frequencies. The radio bridge, the radio arc, and the radio bar all have steep spectra, which can be used to constrain the particle acceleration mechanisms. We argue that the radio bridge could be caused by the re-acceleration of electrons by shock waves that are injected along the filament during the cluster mass assembly. Despite the sensitivity reached by our observations, the emission from the radio halo does not trace the emission of the gas revealed by X-ray observations. We argue that this could be due to the difference in the ratio of kinetic over thermal energy of the intra-cluster gas, suggested by X-ray observations.

We report the detection of diffuse radio emission which might be connected to a large-scale filament of the cosmic web covering a 8deg x 8deg area in the sky, likely associated with a z~0.1 over-density traced by nine massive galaxy clusters. In this work, we present radio observations of this region taken with the Sardinia Radio Telescope. Two of the clusters in the field host a powerful radio halo sustained by violent ongoing mergers and provide direct proof of intra-cluster magnetic fields. In order to investigate the presence of large-scale diffuse radio synchrotron emission in and beyond the galaxy clusters in this complex system, we combined the data taken at 1.4 GHz obtained with the Sardinia Radio Telescope with higher resolution data taken with the NRAO VLA Sky Survey. We found 28 candidate new sources with a size larger and X-ray emission fainter than known diffuse large-scale synchrotron cluster sources for a given radio power. This new population is potentially the tip of the iceberg of a class of diffuse large-scale synchrotron sources associated with the filaments of the cosmic web. In addition, we found in the field a candidate new giant radio galaxy.

We simulate the propagation of cosmic rays at ultra-high energies, $\gtrsim 10^{18}$ eV, in models of extragalactic magnetic fields in constrained simulations of the local Universe. We use constrained initial conditions with the cosmological magnetohydrodynamics code {\sc ENZO}. The resulting models of the distribution of magnetic fields in the local Universe are used in the \crpropa code to simulate the propagation of ultra-high energy cosmic rays. We investigate the impact of six different magneto-genesis scenarios, both primordial and astrophysical, on the propagation of cosmic rays over cosmological distances. Moreover, we study the influence of different source distributions around the Milky Way. Our study shows that different scenarios of magneto-genesis do not have a large impact on the anisotropy measurements of ultra-high energy cosmic rays. However, at high energies above the GZK-limit, there is anisotropy caused by the distribution of nearby sources, independent of the magnetic field model. This provides a chance to identify cosmic ray sources with future full-sky measurements and high number statistics at the highest energies. Finally, we compare our results to the dipole signal measured by the Pierre Auger Observatory. All our source models and magnetic field models could reproduce the observed dipole amplitude with a pure iron injection composition. Our results indicate that the dipole is observed due to clustering of secondary nuclei in direction of nearby sources of heavy nuclei. A light injection composition is disfavoured by the non-observation of anisotropy at energies of $4-8 \rm\ EeV$.

Faraday rotation and synchrotron emission from extragalactic radio sources give evidence for the presence of magnetic fields extending over ~Mpc scales. However, the origin of these fields remains elusive. With new high-resolution grid simulations we studied the growth of magnetic fields in a massive galaxy cluster that in several aspects is similar to the Coma cluster. We investigated models in which magnetic fields originate from primordial seed fields with comoving strengths of 0.1 nG at redshift z=30. The simulations show evidence of significant magnetic field amplification. At the best spatial resolution (3.95 kpc), we are able to resolve the scale where magnetic tension balances the bending of magnetic lines by turbulence. This allows us to observe the final growth stage of the small-scale dynamo. To our knowledge this is the first time that this is seen in cosmological simulations of the intracluster medium. Our mock observations of Faraday Rotation provide a good match to observations of the Coma cluster. However, the distribution of magnetic fields shows strong departures from a simple Maxwellian distribution, suggesting that the three-dimensional structure of magnetic fields in real clusters may be significantly different than what is usually assumed when inferring magnetic field values from rotation measure observations.

The origin of extragalactic magnetic fields is still poorly understood. Based on a dedicated suite of cosmological magneto-hydrodynamical simulations with the ENZO code we have performed a survey of different models that may have caused present-day magnetic fields in galaxies and galaxy clusters. The outcomes of these models differ in cluster outskirts, filaments, sheets and voids and we use these simulations to find observational signatures of magnetogenesis. With these simulations, we predict the signal of extragalactic magnetic fields in radio observations of synchrotron emission from the cosmic web, in Faraday Rotation, in the propagation of Ultra High Energy Cosmic Rays, in the polarized signal from Fast Radio Bursts at cosmological distance and in spectra of distant blazars. In general, primordial scenarios in which present-day magnetic fields originate from the amplification of weak (<nG) uniform seed fields result more homogeneous and relatively easier to observe magnetic fields than than astrophysical scenarios, in which present-day fields are the product of feedback processes triggered by stars and active galaxies. In the near future the best evidence for the origin of cosmic magnetic fields will most likely come from a combination of synchrotron emission and Faraday Rotation observed at the periphery of large-scale structures.

Large-scale extragalactic magnetic fields may induce conversions between very-high-energy photons and axionlike particles (ALPs), thereby shielding the photons from absorption on the extragalactic background light. However, in simplified "cell" models, used so far to represent extragalactic magnetic fields, this mechanism would be strongly suppressed by current astrophysical bounds. Here we consider a recent model of extragalactic magnetic fields obtained from large-scale cosmological simulations. Such simulated magnetic fields would have large enhancement in the filaments of matter. As a result, photon-ALP conversions would produce a significant spectral hardening for cosmic TeV photons. This effect would be probed with the upcoming Cherenkov Telescope Array detector. This possible detection would give a unique chance to perform a tomography of the magnetized cosmic web with ALPs.

We observed the galaxy cluster CIZA J2242.8+5301 with the Sardinia Radio Telescope to provide new constraints on its spectral properties at high frequency. We conducted observations in three frequency bands centred at 1.4 GHz, 6.6 GHz and 19 GHz, resulting in beam resolutions of 14$^{\prime}$, 2.9$^{\prime}$ and 1$^{\prime}$ respectively. These single-dish data were also combined with archival interferometric observations at 1.4 and 1.7 GHz. From the combined images, we measured a flux density of ${\rm S_{1.4GHz}=(158.3\pm9.6)\,mJy}$ for the central radio halo and ${\rm S_{1.4GHz}=(126\pm8)\,mJy}$ and ${\rm S_{1.4GHz}=(11.7\pm0.7)\,mJy}$ for the northern and the southern relic respectively. After the spectral modelling of the discrete sources, we measured at 6.6 GHz ${\rm S_{6.6GHz}=(17.1\pm1.2)\,mJy}$ and ${\rm S_{6.6GHz}=(0.6\pm0.3)\,mJy}$ for the northern and southern relic respectively. Assuming simple diffusive shock acceleration, we interpret measurements of the northern relic with a continuous injection model represented by a broken power-law. This yields an injection spectral index ${\rm \alpha_{inj}=0.7\pm0.1}$ and a Mach number ${\rm M=3.3\pm0.9}$, consistent with recent X-ray estimates. Unlike other studies of the same object, no significant steepening of the relic radio emission is seen in data up to 8.35 GHz. By fitting the southern relic spectrum with a simple power-law (${\rm S_{\nu}\propto\nu^{-\alpha}}$) we obtained a spectral index ${\rm \alpha\approx1.9}$ corresponding to a Mach number (${\rm M\approx1.8}$) in agreement with X-ray estimates. Finally, we evaluated the rotation measure of the northern relic at 6.6 GHz. These results provide new insights on the magnetic structure of the relic, but further observations are needed to clarify the nature of the observed Faraday rotation.

Turbulence generated by large-scale motions during structure formation affects the evolution of the thermal and non-thermal components of the intracluster medium. As enstrophy is a measure of the magnitude of vorticity, we study the generation and evolution of turbulence by analysing the Lagrangian history of enstrophy. For this purpose we combine cosmological simulations carried out with the ENZO-code with our Lagrangian post-processing tool CRaTer. This way we are able to quantify the individual source terms of enstrophy in the course of the accretion of groups onto galaxy clusters. Here we focus on the redshift range from $z=1$ to $z=0$. Finally, we measure the rate of dissipation of turbulence and estimate the resulting amplification of intracluster magnetic fields. We find that compressive and baroclinic motions are the main sources of enstrophy, while stretching motions and dissipation affect most of the ensuing enstrophy evolution. The rate of turbulent dissipation is able to sustain the amplification of intracluster magnetic fields to observed levels.

Giant radio halos are Mpc-scale diffuse radio sources associated with the central regions of galaxy clusters. The most promising scenario to explain the origin of these sources is that of turbulent re-acceleration, in which MeV electrons injected throughout the formation history of galaxy clusters are accelerated to higher energies by turbulent motions mostly induced by cluster mergers. In this Letter, we use the amplitude of density fluctuations in the intracluster medium as a proxy for the turbulent velocity and apply this technique to a sample of 51 clusters with available radio data. Our results indicate a segregation in the turbulent velocity of radio halo and radio quiet clusters, with the turbulent velocity of the former being on average higher by about a factor of two. The velocity dispersion recovered with this technique correlates with the measured radio power through the relation $P_{\rm radio}\propto\sigma_v^{3.3\pm0.7}$, which implies that the radio power is nearly proportional to the turbulent energy rate. Our results provide an observational confirmation of a key prediction of the turbulent re-acceleration model and possibly shed light on the origin of radio halos.

We study the intra-cluster magnetic field in the poor galaxy cluster Abell 194 by complementing radio data, at different frequencies, with data in the optical and X-ray bands. We analyze new total intensity and polarization observations of Abell 194 obtained with the Sardinia Radio Telescope (SRT). We use the SRT data in combination with archival Very Large Array observations to derive both the spectral aging and Rotation Measure (RM) images of the radio galaxies 3C40A and 3C40B embedded in Abell 194. The optical analysis indicates that Abell 194 does not show a major and recent cluster merger, but rather agrees with a scenario of accretion of small groups. Under the minimum energy assumption, the lifetimes of synchrotron electrons in 3C40B measured from the spectral break are found to be 157 Myrs. The break frequency image and the electron density profile inferred from the X-ray emission are used in combination with the RM data to constrain the intra-cluster magnetic field power spectrum. By assuming a Kolmogorov power law power spectrum, we find that the RM data in Abell 194 are well described by a magnetic field with a maximum scale of fluctuations of Lambda_max=64 kpc and a central magnetic field strength of <B0>=1.5 microG. Further out, the field decreases with the radius following the gas density to the power of eta=1.1. Comparing Abell 194 with a small sample of galaxy clusters, there is a hint of a trend between central electron densities and magnetic field strengths.

Observations of diffuse radio emission in galaxy clusters indicate that cosmic-ray electrons are accelerated on $\sim$Mpc scales. However, protons appear to be accelerated less efficiently since their associated hadronic $\gamma$-ray emission has not yet been detected. Inspired by recent particle-in-cell simulations, we study the cosmic-ray production and its signatures under the hypothesis that the efficiency of shock acceleration depends on the Mach number and on the shock obliquity. For this purpose, we combine ENZO cosmological magneto-hydrodynamical simulations with a Lagrangian tracer code to follow the properties of the cosmic rays. Our simulations suggest that the distribution of obliquities in galaxy clusters is random to first order. Quasi-perpendicular shocks are able to accelerate cosmic-ray electrons to the energies needed to produce observable radio emission. However, the $\gamma$-ray emission is lowered by a factor of a few, $\sim$3, if cosmic-ray protons are only accelerated by quasi-parallel shocks, reducing (yet not entirely solving) the tension with the non-detection of hadronic $\gamma$-ray emission by the \textit{Fermi}-satellite.

The emergence of cosmic structure is commonly considered one of the most complex phenomena in Nature. However, this complexity has never been defined nor measured in a quantitative and objective way. In this work we propose a method to measure the information content of cosmic structure and to quantify the complexity that emerges from it, based on Information Theory. The emergence of complex evolutionary patterns is studied with a statistical symbolic analysis of the datastream produced by state-of-the-art cosmological simulations of forming galaxy clusters. This powerful approach allows us to measure how many bits of information are necessary to predict the evolution of energy fields in a statistical way, and it offers a simple way to quantify when, where and how the cosmic gas behaves in complex ways. The most complex behaviors are found in the peripheral regions of galaxy clusters, where supersonic flows drive shocks and large energy fluctuations over a few tens of million years. Describing the evolution of magnetic energy requires at least a twice as large amount of bits than for the other energy fields. When radiative cooling and feedback from galaxy formation are considered, the cosmic gas is overall found to double its degree of complexity. In the future, Cosmic Information Theory can significantly increase our understanding of the emergence of cosmic structure as it represents an innovative framework to design and analyze complex simulations of the Universe in a simple, yet powerful way.

Galaxy cluster merger shocks are the main agent for the thermalization of the intracluster medium and the energization of cosmic ray particles in it. Shock propagation changes the state of the tenuous intracluster plasma, and the corresponding signal variations are measurable with the current generation of X-ray and Sunyaev-Zel'dovich (SZ) effect instruments. Additionally, non-thermal electrons (re-)energized by the shocks sometimes give rise to extended and luminous synchrotron sources known as radio relics, which are prominent indicators of shocks propagating roughly in the plane of the sky. In this short review, we discuss how the joint modeling of the non-thermal and thermal signal variations across radio relic shock fronts is helping to advance our knowledge of the gas thermodynamical properties and magnetic field strengths in the cluster outskirts. We describe the first use of the SZ effect to measure the Mach numbers of relic shocks, for both the nearest (Coma) and the farthest (El Gordo) clusters with known radio relics.

We present ALMA measurements of a merger shock using the thermal Sunyaev-Zel'dovich (SZ) effect signal, at the location of a radio relic in the famous El Gordo galaxy cluster at $z \approx 0.9$. Multi-wavelength analysis in combination with the archival Chandra data and a high-resolution radio image provides a consistent picture of the thermal and non-thermal signal variation across the shock front and helps to put robust constraints on the shock Mach number as well as the relic magnetic field. We employ a Bayesian analysis technique for modeling the SZ and X-ray data self-consistently, illustrating respective parameter degeneracies. Combined results indicate a shock with Mach number ${\cal M} = 2.4^{+1.3}_{-0.6}$, which in turn suggests a high value of the magnetic field (of the order of $4-10 ~\mu$G) to account for the observed relic width at 2 GHz. At roughly half the current age of the universe, this is the highest-redshift direct detection of a cluster shock to date, and one of the first instances of an ALMA-SZ observation in a galaxy cluster. It shows the tremendous potential for future ALMA-SZ observations to detect merger shocks and other cluster substructures out to the highest redshifts.

1) Background: the budget of non-thermal energy in galaxy clusters is not well constrained, owing to the observational and theoretical difficulties in studying these diluted plasmas on large scales. 2) Method: we use recent cosmological simulations with complex physics in order to connect the emergence of non-thermal energy to the underlying evolution of gas and dark matter. 3) Results: the impact of non-thermal energy (e.g. cosmic rays, magnetic fields and turbulent motions) is found to increase in the outer region of galaxy clusters. Within numerical and theoretical uncertainties, turbulent motions dominate the budget of non-thermal energy in most of the cosmic volume. 4) Conclusion: assessing the distribution non-thermal energy in galaxy clusters is crucial to perform high-precision cosmology in the future. Constraining the level of non-thermal energy in cluster outskirts will improve our understanding of the acceleration of relativistic particles by cosmic shocks and of the origin of extragalactic magnetic fields.

We review the observational and theoretical constraints on extragalactic magnetic fields across cosmic environment. In the next decade, the combination of sophisticated numerical simulations and various observational probes might succeed in constraining the still elusive origin of magnetic fields on the largest scales in the Universe.

Weak shocks in the intracluster medium may accelerate cosmic-ray protons and cosmic-ray electrons differently depending on the angle between the upstream magnetic field and the shock normal. In this work, we investigate how shock obliquity affects the production of cosmic rays in high-resolution simulations of galaxy clusters. For this purpose, we performed a magneto-hydrodynamical simulation of a galaxy cluster using the mesh refinement code \enzo. We use Lagrangian tracers to follow the properties of the thermal gas, the cosmic rays and the magnetic fields over time. We tested a number of different acceleration scenarios by varying the obliquity-dependent acceleration efficiencies of protons and electrons, and by examining the resulting hadronic $\gamma$-ray and radio emission. We find that the radio emission does not change significantly if only quasi-perpendicular shocks are able to accelerate cosmic-ray electrons. Our analysis suggests that radio emitting electrons found in relics have been typically shocked many times before $z=0$. On the other hand, the hadronic $\gamma$-ray emission from clusters is found to decrease significantly if only quasi-parallel shocks are allowed to accelerate cosmic-ray protons. This might reduce the tension with the low upper limits on $\gamma$-ray emission from clusters set by the \textit{Fermi}-satellite.

Turbulence is a key ingredient for the evolution of the intracluster medium, whose properties can be predicted with high resolution numerical simulations. We present initial results on the generation of solenoidal and compressive turbulence in the intracluster medium during the formation of a small-size cluster using highly resolved, non-radiative cosmological simulations, with a refined monitoring in time. In this first of a series of papers, we closely look at one simulated cluster whose formation was distinguished by a merger around $z \sim 0.3$. We separate laminar gas motions, turbulence and shocks with dedicated filtering strategies and distinguish the solenoidal and compressive components of the gas flows using Hodge-Helmholtz decomposition. Solenoidal turbulence dominates the dissipation of turbulent motions ($\sim 95\%$) in the central cluster volume at all epochs. The dissipation via compressive modes is found to be more important ($\sim 30 \%$ of the total) only at large radii ($\geq 0.5 ~r_{\rm vir}$) and close to merger events. We show that enstrophy (vorticity squared) is good proxy of solenoidal turbulence. All terms ruling the evolution of enstrophy (i.e. baroclinic, compressive, stretching and advective terms) are found to be significant, but in amounts that vary with time and location. Two important trends for the growth of enstrophy in our simulation are identified: first, enstrophy is continuously accreted into the cluster from the outside, and most of that accreted enstrophy is generated near the outer accretion shocks by baroclinic and compressive processes. Second, in the cluster interior vortex stretching is dominant, although the other terms also contribute substantially.

We use the CRPropa code to simulate the propagation of ultra high energy cosmic rays (with energy $\geq 10^{18} \rm eV$ and pure proton composition) through extragalactic magnetic fields that have been simulated with the cosmological ENZO code.We test both primordial and astrophysical magnetogenesis scenarios in order to investigate the impact of different magnetic field strengths in clusters, filaments and voids on the deflection of cosmic rays propagating across cosmological distances. We also study the effect of different source distributions of cosmic rays around simulated Milky-Way like observers. Our analysis shows that the arrival spectra and anisotropy of events are rather insensitive to the distribution of extragalactic magnetic fields, while they are more affected by the clustering of sources within a $\sim 50$ Mpc distance to observers. Finally, we find that in order to reproduce the observed degree of isotropy of cosmic rays at $\sim $ EeV energies, the average magnetic fields in cosmic voids must be $\sim 0.1 \rm \ nG$, providing limits on the strength of primordial seed fields.

Radio relics are Mpc-scale diffuse radio sources at the peripheries of galaxy clusters which are thought to trace outgoing merger shocks. We present XMM-Newton and Suzaku observations of the galaxy cluster Abell 2744 (z=0.306), which reveal the presence of a shock front 1.5 Mpc East of the cluster core. The surface-brightness jump coincides with the position of a known radio relic. Although the surface-brightness jump indicates a weak shock with a Mach number $\mathcal{M}=1.7_{-0.3}^{+0.5}$, the plasma in the post-shock region has been heated to a very high temperature ($\sim13$ keV) by the passage of the shock wave. The low acceleration efficiency expected from such a weak shock suggests that mildly relativistic electrons have been re-accelerated by the passage of the shock front.

Despite containing about a half of the total matter in the Universe, at most wavelengths the filamentary structure of the cosmic web is difficult to observe. In this work, we use large unigrid cosmological simulations to investigate how the geometrical, thermodynamical and magnetic properties of cosmological filaments vary with mass and redshift (z $\leq 1$). We find that the average temperature, length, volume and magnetic field of filaments are tightly log-log correlated with the underlying total gravitational mass. This reflects the role of self-gravity in shaping their properties and enables statistical predictions of their observational properties based on their mass. We also focus on the properties of the simulated population of galaxy-sized halos within filaments, and compare their properties to the results obtained from the spectroscopic GAMA survey. Simulated and observed filaments with the same length are found to contain an equal number of galaxies, with very similar distribution of halo masses. The total number of galaxies within each filament and the total/average stellar mass in galaxies can now be used to predict also the large-scale properties of the gas in the host filaments across tens or hundreds of Mpc in scale. These results are the first steps towards the future use of galaxy catalogues in order to select the best targets for observations of the warm-hot intergalactic medium.

(Abridged) Radio relics in galaxy clusters are believed to be associated with powerful shock fronts that originate during cluster mergers, and are a testbed for the acceleration of relativistic particles in the intracluster medium. Recently, radio relic observations have pushed into the cm-wavelength domain (1-30 GHz) where a break from the standard synchrotron power-law spectrum has been found, most noticeably in the famous 'Sausage' relic. In this paper, we point to an important effect that has been ignored or considered insignificant while interpreting these new high-frequency radio data, namely the contamination due to the Sunyaev-Zel'dovich (SZ) effect that changes the observed synchrotron flux. Even though the radio relics reside in the cluster outskirts, the shock-driven pressure boost increases the SZ signal locally by roughly an order of magnitude. The resulting flux contamination for some well-known relics are non-negligible already at 10 GHz, and at 30 GHz the observed synchrotron fluxes can be diminished by a factor of several from their true values. Interferometric observations are not immune to this contamination, since the change in the SZ signal occurs roughly at the same length scale as the synchrotron emission, although there the flux loss is less severe than single-dish observations. We present a simple analytical approximation for the synchrotron-to-SZ flux ratio, based on a theoretical radio relic model that connects the non-thermal emission to the thermal gas properties, and show that by measuring this ratio one can potentially estimate the relic magnetic fields or the particle acceleration efficiency.

We study the acceleration of cosmic rays by collisionless structure formation shocks with ENZO grid simulations. Data from the FERMI satellite enable the use of galaxy clusters as a testbed for particle acceleration models. Based on advanced cosmological simulations that include different prescriptions for gas and cosmic rays physics, we use the predicted {\gamma}-ray emission to constrain the shock acceleration efficiency. We infer that the efficiency must be on average <0.1% for cosmic shocks, particularly for the 2<M<5 merger shocks that are mostly responsible for the thermalisation of the intracluster medium. These results emerge, both, from non-radiative and radiative runs including feedback from active galactic nuclei, as well as from zoomed resimulations of a cluster resembling MACSJ1752.0+0440. The limit on the acceleration efficiency we report is lower than what has been assumed in the literature so far. Combined with the information from radio emission in clusters, it appears that a revision of the present understanding of shock acceleration in the ICM is unavoidable.

We study the challenges to detect the cosmic web at radio wavelengths with state-of-the-art cosmological simulations of extragalactic magnetic fields. The incoming generation of radio surveys operating at low frequency, like LOFAR, SKA-LOW and MWA will have the best chance to detect the large-scale, low surface brightness emission from the shocked cosmic web. The detected radio emission will enable to constrain the average magnetisation level of the gas in filaments and the acceleration efficiency of electrons by strong shocks. In case of detections, through statistical modelling (e.g. correlation functions) it will be possible to discriminate among competing scenarios for the magnetisation of large-scale structures (i.e. astrophysical versus primordial scenarios), making radio surveys an important probe of cosmic magnetogenesis.

We measure the thermodynamic properties of cluster outer regions to provide constraints on the processes that rule the formation of large scale structures. We derived the thermodynamic properties of the intracluster gas (temperature, entropy) by combining the SZ thermal pressure from Planck and the X-ray gas density from ROSAT. This method allowed us to reconstruct for the first time temperature and entropy profiles out to the virial radius and beyond in a large sample of objects. At variance with several recent Suzaku studies, we find that the entropy rises steadily with radius, albeit at at a somewhat lower rate than predicted by self-similar expectations. We note significant differences between relaxed, cool-core systems and unrelaxed clusters in the outer regions. Relaxed systems appear to follow the self-similar expectations more closely than perturbed objects. Our results indicate that the well-known entropy excess observed in cluster cores extends well beyond the central regions. When correcting for the gas depletion, the observed entropy profiles agree with the prediction from gravitational collapse only, especially for cool-core clusters.

Several recent studies used the hot gas fraction of galaxy clusters as a standard ruler to constrain dark energy, which provides competitive results compared to other techniques. This method, however, relies on the assumption that the baryon fraction in clusters agrees with the cosmic value Omega_b/Omega_m, and does not differ from one system to another. We test this hypothesis by measuring the gas mass fraction over the entire cluster volume in a sample of local clusters. Combining the SZ thermal pressure from Planck and the X-ray gas density from ROSAT, we measured for the first time the average gas fraction (fgas) out to the virial radius and beyond in a large sample of clusters. We also obtained azimuthally-averaged measurements of the gas fraction for 18 individual systems, which we used to compute the scatter of fgas around the mean value at different radii and its dependence on the cluster's temperature. The gas mass fraction increases with radius and reaches the cosmic baryon fraction close to R200. At R200, we measure fgas,200=0.176+/-0.009. We find significant differences between the baryon fraction of relaxed, cool-core (CC) systems and unrelaxed, non-cool core (NCC) clusters in the outer regions. In average, the gas fraction in NCC clusters slightly exceeds the cosmic baryon fraction, while in CC systems the gas fraction converges to the expected value when accounting for the stellar content, without any evidence for variations from one system to another. We find that fgas estimates in NCC systems slightly disagree with the cosmic value approaching R200. This result could be explained either by a violation of the assumption of hydrostatic equilibrium or by an inhomogeneous distribution of the gas mass. Conversely, cool-core clusters are found to provide reliable constraints on fgas at overdensities >200, which makes them suitable for cosmological studies.

We present the analysis of a local (z = 0.04 - 0.2) sample of 31 galaxy clusters with the aim of measuring the density of the X-ray emitting gas in cluster outskirts. We compare our results with numerical simulations to set constraints on the azimuthal symmetry and gas clumping in the outer regions of galaxy clusters. We exploit the large field-of-view and low instrumental background of ROSAT/PSPC to trace the density of the intracluster gas out to the virial radius. We perform a stacking of the density profiles to detect a signal beyond r200 and measure the typical density and scatter in cluster outskirts. We also compute the azimuthal scatter of the profiles with respect to the mean value to look for deviations from spherical symmetry. Finally, we compare our average density and scatter profiles with the results of numerical simulations. As opposed to some recent Suzaku results, and confirming previous evidence from ROSAT and Chandra, we observe a steepening of the density profiles beyond \sim r500. Comparing our density profiles with simulations, we find that non-radiative runs predict too steep density profiles, whereas runs including additional physics and/or treating gas clumping are in better agreement with the observed gas distribution. We report for the first time the high-confidence detection of a systematic difference between cool-core and non-cool core clusters beyond \sim 0.3r200, which we explain by a different distribution of the gas in the two classes. Beyond \sim r500, galaxy clusters deviate significantly from spherical symmetry, with only little differences between relaxed and disturbed systems. We find good agreement between the observed and predicted scatter profiles, but only when the 1% densest clumps are filtered out in the simulations. [Abridged]

Radio halos are synchrotron radio sources detected in some massive galaxy clusters. Their Mpc-size indicates that (re)acceleration processes are taking place in the host cluster. X-ray catalogues of galaxy clusters have been used in the past to search for radio halos and to understand their connection with cluster-cluster mergers and with the thermal component of the intra-cluster medium. More recently, the Sunyaev-Zel'dovich effect has been proven to be a better route to search for massive clusters in a wider redshift range. With the aim of discovering new radio halos and understanding their connection with cluster-cluster mergers, we have selected from the Planck Early source catalog the most massive clusters, and we have observed with the Giant Metrewave Radio Telescope at 323 MHz those objects for which deep observations were not available. We have discovered new peculiar radio emission in three of the observed clusters finding: (i) a radio halo in the cluster RXCJ0949.8+1708; (ii) extended emission in Abell 1443 that we classify as a radio halo plus a radio relic, with a bright filament embedded in the radio halo; (iii) low-power radio emission is found in CIZA J1938.3+5409 which is ten times below the radio - X-ray correlation, and represents the first direct detection of the radio emission in the "upper-limit" region of the radio - X-ray diagram. We discuss the properties of these new radio halos in the framework of theoretical models for the radio emission.

Using a new parallel algorithm implemented within the VisIt framework, we analysed large cosmological grid simulations to study the properties of baryons in filaments. The procedure allows us to build large catalogues with up to $\sim 3 \cdot 10^4$ filaments per simulated volume and to investigate the properties of cosmic filaments for very large volumes at high resolution (up to $300^3 ~\rm Mpc^3$ simulated with $2048^3$ cells). We determined scaling relations for the mass, volume, length and temperature of filaments and compared them to those of galaxy clusters. The longest filaments have a total length of about $200 ~\rm Mpc$ with a mass of several $10^{15} M_{\odot}$. We also investigated the effects of different gas physics. Radiative cooling significantly modifies the thermal properties of the warm-hot-intergalactic medium of filaments, mainly by lowering their mean temperature via line cooling. On the other hand, powerful feedback from active galactic nuclei in surrounding halos can heat up the gas in filaments. The impact of shock-accelerated cosmic rays from diffusive shock acceleration on filaments is small and the ratio of between cosmic ray and gas pressure within filaments is of the order of $\sim 10-20$ percent.

The cosmic web contains a large fraction of the total gas mass in the universe but is difficult to detect at most wavelengths. Synchrotron emission from shock-accelerated electrons may offer the chance of imaging the cosmic web at radio wavelengths. In this work we use 3D cosmological ENZO-MHD simulations (combined with a post-processing renormalisation of the magnetic field to bracket for missing physical ingredients and resolution effects) to produce models of the radio emission from the cosmic web. In post-processing we study the capabilities of 13 large radio surveys to detect this emission. We find that surveys by LOFAR, SKA1-LOW and MWA have a chance of detecting the cosmic web, provided that the magnetisation level of the tenuous medium in filaments is of the order of 1% of the thermal gas energy.

Radio relics in galaxy clusters are associated with powerful shocks that (re)accelerate relativistic electrons. It is widely believed that the acceleration proceeds via diffusive shock acceleration. In the framework of thermal leakage, the ratio of the energy in relativistic electrons to the energy in relativistic protons should should be smaller than $K_{\rm e/p} \sim 10^{-2}$. The relativistic protons interact with the thermal gas to produce $\gamma$-rays in hadronic interactions. Combining observations of radio relics with upper limits from $\gamma$-ray observatories can constrain the ratio $K_{\rm e/p}$. In this work we selected 10 galaxy clusters that contain double radio relics, and derive new upper limits from the stacking of $\gamma$-ray observations by FERMI. We modelled the propagation of shocks using a semi-analytical model, where we assumed a simple geometry for shocks and that cosmic ray protons are trapped in the intracluster medium. Our analysis shows that diffusive shock acceleration has difficulties in matching simultaneously the observed radio emission and the constraints imposed by FERMI, unless the magnetic field in relics is unrealistically large ($\gg 10 ~\rm \mu G$). In all investigated cases (also including realistic variations of our basic model and the effect of re-acceleration) the mean emission of the sample is of the order of the stacking limit by FERMI, or larger. These findings put tension on the commonly adopted model for the powering of radio relics, and imply that the relative acceleration efficiency of electrons and protons is at odds with predictions of diffusive shock acceleration, requiring $K_{\rm e/p} \geq 10-10^{-2}$.

The detection of the diffuse gas component of the cosmic web remains a formidable challenge. In this work we study synchrotron emission from the cosmic web with simulated SKA1 observations, which can represent an fundamental probe of the warm-hot intergalactic medium. We investigate radio emission originated by relativistic electrons accelerated by shocks surrounding cosmic filaments, assuming diffusive shock acceleration and as a function of the (unknown) large-scale magnetic fields. The detection of the brightest parts of large ($>10 \rm Mpc$) filaments of the cosmic web should be within reach of the SKA1-LOW, if the magnetic field is at the level of a $\sim 10$ percent equipartition with the thermal gas, corresponding to $\sim 0.1 \mu G$ for the most massive filaments in simulations. In the course of a 2-years survey with SKA1-LOW, this will enable a first detection of the "tip of the iceberg" of the radio cosmic web, and allow for the use of the SKA as a powerful tool to study the origin of cosmic magnetism in large-scale structures. On the other hand, the SKA1-MID and SKA1-SUR seem less suited for this science case at low redshift ($z \leq 0.4$), owing to the missing short baselines and the consequent lack of signal from the large-scale brightness fluctuations associated with the filaments. In this case only very long exposures ($\sim 1000$ hr) may enable the detection of $\sim 1-2$ filament for field of view in the SKA1-SUR PAF Band1.

We investigate the possibility for the SKA to detect and study the magnetic fields in galaxy clusters and in the less dense environments surrounding them using Faraday Rotation Measures. To this end, we produce 3-dimensional magnetic field models for galaxy clusters of different masses and in different stages of their evolution, and derive mock rotation measure observations of background radiogalaxies. According to our results, already in phase I, we will be able to infer the magnetic field properties in galaxy clusters as a function of the cluster mass, down to $10^{13}$ solar-masses. Moreover, using cosmological simulations to model the gas density, we have computed the expected rotation measure through shock-fronts that occur in the intra-cluster medium during cluster mergers. The enhancement in the rotation measure due to the density jump will permit to constraint the magnetic field strength and structure after the shock passage. SKA observations of polarised sources located behind galaxy clusters will answer several questions about the magnetic field strength and structure in galaxy clusters, and its evolution with cosmic time.

Galaxy clusters are known to host a variety of extended radio sources: tailed radio galaxies whose shape is modelled by the interaction with the intra-cluster medium (ICM); radio bubbles filling cavities in the ICM distribution and rising buoyantly through the thermal gas; diffuse giant radio sources ("halos" and "relics") revealing the presence of relativistic electrons and magnetic fields in the intra-cluster volume. It is currently the subject of an active debate how the non-thermal components that we observe at radio wavelengths affect the physical properties of the ICM and depend on the dynamical state of galaxy clusters. In this work we start our SKA1 feasibility study of the "radio cluster zoo" through simulations of a typical radio-loud cluster, hosting several bright tailed radio galaxies and a diffuse radio halo. Realistic simulations of SKA1 observations are obtained through the MeqTrees software. A new deconvolution algorithm, based on sparse representations and optimised for the detection of faint diffuse astronomical sources, is tested and compared to the classical CLEAN method.

The reconstruction of galaxy cluster's gas density profiles is usually performed by assuming spherical symmetry and averaging the observed X-ray emission in circular annuli. In the case of a very inhomogeneous and asymmetric gas distribution, this method has been shown to return biased results in numerical simulations because of the $n^2$ dependence of the X-ray emissivity. We propose a method to recover the true density profiles in the presence of inhomogeneities, based on the derivation of the azimuthal median of the surface brightness in concentric annuli. We demonstrate the performance of this method with numerical simulations, and apply it to a sample of 31 galaxy clusters in the redshift range 0.04-0.2 observed with ROSAT/PSPC. The clumping factors recovered by comparing the mean and the median are mild and show a slight trend of increasing bias with radius. For $R<R_{500}$, we measure a clumping factor $\sqrt{C}<1.1$, which indicates that the thermodynamic properties and hydrostatic masses measured in this radial range are only mildly affected by this effect. Comparing our results with three sets of hydrodynamical numerical simulations, we found that non-radiative simulations significantly overestimate the level of inhomogeneities in the ICM, while the runs including cooling, star formation, and AGN feedback reproduce the observed trends closely. Our results indicate that most of the accretion of X-ray emitting gas is taking place in the diffuse, large-scale accretion patterns rather than in compact structures.

A number of scientific numerical codes can currently exploit GPUs with remarkable performance. In astrophysics, Enzo and Ramses are prime examples of such applications. The two codes have been ported to GPUs adopting different strategies and programming models, Enzo adopting CUDA and Ramses using OpenACC. We describe here the different solutions used for the GPU implementation of both cases. Performance benchmarks will be presented for Ramses. The results of the usage of the more mature GPU version of Enzo, adopted for a scientific project within the CHRONOS programme, will be summarised.

We compare two different codes for simulations of cosmological structure formation to investigate the sensitivity of hydrodynamical instabilities to numerics, in particular, the hydro solver and the application of adaptive mesh refinement (AMR). As a simple test problem, we consider an initially spherical gas cloud in a wind, which is an idealized model for the merger of a subcluster or galaxy with a big cluster. Based on an entropy criterion, we calculate the mass stripping from the subcluster as a function of time. Moreover, the turbulent velocity field is analyzed with a multi-scale filtering technique. We find remarkable differences between the commonly used PPM solver with directional splitting in the Enzo code and an unsplit variant of PPM in the Nyx code, which demonstrates that different codes can converge to systematically different solutions even when using uniform grids. For the test case of an unbound cloud, AMR simulations reproduce uniform-grid results for the mass stripping quite well, although the flow realizations can differ substantially. If the cloud is bound by a static gravitational potential, however, we find strong sensitivity to spurious fluctuations which are induced at the cutoff radius of the potential and amplified by the bow shock. This gives rise to substantial deviations between uniform-grid and AMR runs performed with Enzo, while the mass stripping in Nyx simulations of the subcluster is nearly independent of numerical resolution and AMR. Although many factors related to numerics are involved, our study indicates that unsplit solvers with advanced flux limiters help to reduce grid effects and to keep numerical noise under control, which is important for hydrodynamical instabilities and turbulent flows.

The amplification of primordial magnetic fields via a small-scale turbulent dynamo during structure formation might be able to explain the observed magnetic fields in galaxy clusters. The magnetisation of more tenuous large-scale structures such as cosmic filaments is more uncertain, as it is challenging for numerical simulations to achieve the required dynamical range. In this work, we present magneto-hydrodynamical cosmological simulations on large uniform grids to study the amplification of primordial seed fields in the intracluster medium (ICM) and in the warm-hot-intergalactic medium (WHIM). In the ICM, we confirm that turbulence caused by structure formation can produce a significant dynamo amplification, even if the amplification is smaller than what is reported in other papers. In the WHIM inside filaments, we do not observe significant dynamo amplification, even though we achieve Reynolds numbers of $R_{\rm e} \sim 200-300$. The maximal amplification for large filaments is of the order of $\sim 100$ for the magnetic energy, corresponding to a typical field of a few $\sim \rm nG$ starting from a primordial weak field of $10^{-10}$ G (comoving). In order to start a small-scale dynamo, we found that a minimum of $\sim 10^2$ resolution elements across the virial radius of galaxy clusters was necessary. In filaments we could not find a minimum resolution to set off a dynamo. This stems from the inefficiency of supersonic motions in the WHIM in triggering solenoidal modes and small-scale twisting of magnetic field structures. Magnetic fields this small will make it hard to detect filaments in radio observations.

Radio relics are diffuse synchrotron sources in galaxy clusters that are believed to trace large-scale shock waves. We have discovered a new double radio relic system in PSZ1 G096.89+24.17 (z=0.3) and have carried out a full-polarization radio observation using the Westerbork Synthesis Radio Telescope (WSRT) at 1.4 GHz. The observation revealed the presence of two relics located on the two diametrically opposite sides of the cluster and hints of a central radio halo. The linear sizes of the relics are $\sim 0.9$ and $\sim 1.4$ Mpc. We carried out an analysis of all known double radio relics by using radio, X-ray and Sunyaev-Zeldovich (SZ) data. We find that the radio luminosity of double relics is a steep function of the cluster mass, with $L_R \propto M^{2.83\pm0.39}$. If we include single radio relics, this relation is maintained. This dependence has implications for the origin of magnetic fields at the relic's locations.

Radio relics in galaxy clusters are thought to be associated with powerful shock waves that accelerate particles via diffusive shock acceleration (DSA). Among the particles accelerated by DSA, relativistic protons should outnumber electrons by a large factor. While the relativistic electrons emit synchrotron emission detectable in the radio band, the protons interact with the thermal gas to produce gamma-rays in hadronic interactions. Using simple models for the propagation of shock waves through clusters, the distribution of thermal gas and the efficiency of DSA, we find that the resulting hadronic $\gamma$-ray emission lies very close or above the upper limits from the FERMI data on nearby clusters. This suggests that the relative acceleration efficiency of electrons and protons is at odds with predictions from DSA. The inclusion of re-accelerated "fossil" particles does not seem to solve the problem. Our study highlights a possible tension of the commonly assumed scenario for the formation of radio relics and we discuss possible solutions to the problem.

Non-thermal (relativistic) particles are injected into the cosmos by structure formation shock waves, active galactic nuclei and stellar explosions. We present a suite of unigrid cosmological simulations (up to $2048^3$) using a two-fluid model in the grid code ENZO. The simulations include the dynamical effects of cosmic-ray (CR) protons and cover a range of theoretically motivated acceleration efficiencies. For the bulk of the cosmic volume the modelling of CR processes is rather stable with respect to resolution, provided that a minimum (cell) resolution of $\approx 100 ~\rm kpc/h$ is employed. However, the results for the innermost cluster regions depend on the assumptions for the baryonic physics. Inside clusters, non-radiative runs at high resolution tend to produce an energy density of CRs that are below available upper limits from the FERMI satellite, while the radiative runs are found to produce a higher budget of CRs. We show that weak ($M \leq 3-5$) shocks and shock-reacceleration are crucial to set the level of CRs in the innermost region of clusters, while in the outer regions the level of CR energy is mainly set via direct injection by stronger shocks, and is less sensitive to cooling and feedback from active galactic nuclei and supernovae.

As the nodes of the cosmic web, clusters of galaxies trace the large-scale distribution of matter in the Universe. They are thus privileged sites in which to investigate the complex physics of structure formation. However, the complete story of how these structures grow, and how they dissipate the gravitational and non-thermal components of their energy budget over cosmic time, is still beyond our grasp. Fundamental questions such as How do hot diffuse baryons accrete and dynamically evolve in dark matter potentials? How and when was the energy that we observe in the ICM generated and distributed? Where and when are heavy elements produced and how are they circulated? are still unanswered. Most of the cluster baryons exists in the form of a diffuse, hot, metal-enriched plasma that radiates primarily in the X-ray band (the intracluster medium, ICM), allowing the X-ray observations of the evolving cluster population to provide a unique opportunity to address these topics. Athena+ with its large collecting area and unprecedented combination of high spectral and angular resolution offers the only way to make major advances in answering these questions. Athena+ will show how the baryonic gas evolves in the dark matter potential wells by studying the motions and turbulence in the ICM. Athena+ will be able to resolve the accreting region both spatially and spectroscopically, probing the true nature and physical state of the X-ray emitting plasma. Athena+ has the capabilities to permit a definitive understanding of the formation and evolution of large-scale cosmic structure through the study of the cluster population.

The aim of this work is to probe the magnetic field properties in relics and infall regions of galaxy clusters using Faraday Rotation Measures. We present Very Large Array multi-frequency observations of seven sources in the region South-West of the Coma cluster, where the infalling group NGC4839 and the relic 1253+275 are located. The Faraday Rotation Measure maps for the observed sources are derived and analysed to study the magnetic field in the South-West region of Coma. We discuss how to interpret the data by comparing observed and mock rotation measures maps that are produced simulating different 3-dimensional magnetic field models. The magnetic field model that gives the best fit to the Coma central region underestimates the rotation measure in the South-West region by a factor ~6, and no significant jump in the rotation measure data is found at the position of the relic. We explore different possibilities to reconcile observed and mock rotation measure trends, and conclude that an amplification of the magnetic field along the South-West sector is the most plausible solution. Our data together with recent X-ray estimates of the gas density obtained with Suzaku suggest that a magnetic field amplification by a factor ~3 is required throughout the entire South-West region in order to reconcile real and mock rotation measures trends. The magnetic field in the relic region is inferred to be ~2 microG, consistent with Inverse Compton limits.

The spatial distribution of gas matter inside galaxy clusters is not completely smooth, but may host gas clumps associated with substructures. These overdense gas substructures are generally a source of unresolved bias of X-ray observations towards high density gas, but their bright luminosity peaks may be resolved sources within the ICM, that deep X-ray exposures may be (already) capable to detect. In this paper we aim at investigating both features, using a set of high-resolution cosmological simulations with ENZO. First, we monitor how the bias by unresolved gas clumping may yield incorrect estimates of global cluster parameters and affects the measurements of baryon fractions by X-ray observations. We find that based on X-ray observations of narrow radial strips, it is difficult to recover the real baryon fraction to better than 10 - 20 percent uncertainty. Second, we investigated the possibility of observing bright X-ray clumps in the nearby Universe (z<=0.3). We produced simple mock X-ray observations for several instruments (XMM, Suzaku and ROSAT) and extracted the statistics of potentially detectable bright clumps. Some of the brightest clumps predicted by simulations may already have been already detected in X- ray images with a large field of view. However, their small projected size makes it difficult to prove their existence based on X-ray morphology only. Preheating, AGN feedback and cosmic rays are found to have little impact on the statistical properties of gas clumps.

We investigate the observable effects of feedback from Active Galactic Nuclei (AGN) on non-thermal components of the intracluster medium (ICM). We have modelled feedback from AGN in cosmological simulations with the adaptive mesh refinement code ENZO, investigating three types of feedback that are sometimes called quasar, jet and radio mode. Using a small set of galaxy clusters simulated at high resolution, we model the injection and evolution of Cosmic Rays, as well as their effects on the thermal plasma. By comparing, both, the profiles of thermal gas to observed profiles from the ACCEPT sample, and the secondary gamma-ray emission to the available upper limits from FERMI, we discuss how the combined analysis of these two observables can constrain the energetics and mechanisms of feedback models in clusters. Those modes of AGN feedback that provide a good match to X-ray observations, yield a gamma-ray luminosity resulting from secondary cosmic rays that is about below the available upper limits from FERMI. Moreover, we investigate the injection of turbulent motions into the ICM from AGN, and the detectability of these motions via the analysis of line broadening of the Fe XXIII line. In the near future, deeper observations/upper-limits of non-thermal emissions from galaxy clusters will yield stringent constraints on the energetics and modes of AGN feedback, even at early cosmic epochs.

We have performed 323 MHz observations with the Giant Metrewave Radio Telescope of the most promising candidates selected from the MACS catalog. The aim of the work is to extend our knowledge of the radio halo and relic populations to z>0.3, the epoch in which massive clusters formed. In MACSJ1149.5+2223 and MACSJ1752.1+4440, we discovered two double-relic systems with a radio halo, and in MACSJ0553.4-3342 we found a radio halo. Archival Very Large Array observations and Westerbork Synthesis Radio Telescope observations have been used to study the polarization and spectral index properties. The radio halo in MACSJ1149.5+2223 has the steepest spectrum ever found so far in these objects (alpha > 2). The double relics in MACSJ1149.5+2223 are peculiar in their position that is misaligned with the main merger axis. The relics are polarized up to 30% and 40% in MACSJ1149.5+2223 and MACSJ1752.040+44, respectively. In both cases, the magnetic field is roughly aligned with the relics' main axes. The spectra in the relics in MACSJ1752.040+44 steepen towards the cluster centre, in agreement with model expectations. X-ray data on MACSJ0553.4-3342 suggests that this cluster is undergoing a major merger, with the merger axis close to the plane of the sky. The cores of the disrupted clusters have just passed each other, but no radio relic is detected in this system. If turbulence is responsible for the radio emission, we argue that it must develop before the core passage. A comparison of double relic plus halo system with cosmological simulations allows a simultaneous estimate of the acceleration efficiencies at shocks (to produce relics) and of turbulence (to produce the halo).

We have designed a simple multi-scale method that identifies turbulent motions in hydrodynamical grid simulations. The method does not assmume ant a-priori coherence scale to distinguish laminar and turbulent flows. Instead, the local mean velocity field around each cell is reconstructed with a multi-scale filtering technique, yielding the maximum scale of turbulent eddies by means of iterations. The method is robust, fast and easily applicable to any grid simulation. We present here the application of this technique to the study of spatial and spectral properties of turbulence in the intra cluster medium, measuring turbulent diffusion and anisotropy of the turbulent velocity field for a variety of driving mechanisms: a) accretion of matter in galaxy clusters (simulated with ENZO); b) sloshing motions around cool-cores (simulated with FLASH); c) jet outflows from active galactic nuclei, AGN (simulated with FLASH). The turbulent velocities driven by matter accretion in galaxy clusters are mostly tangential in the inner regions (inside the cluster virial radius) and isotropic in regions close to the virial radius. The same is found for turbulence excited by cool core sloshing, while the jet outflowing from AGN drives mostly radial turbulence motions near its sonic point and beyond. Turbulence leads to a diffusivity in the range =10^29-10^30 cm^2/s in the intra cluster medium. On average, the energetically dominant mechanism of turbulence driving in the intra cluster medium is represented by accretion of matter and major mergers during clusters evolution.

We present a numerical scheme, implemented in the cosmological adaptive mesh refinement code ENZO, to model the injection of Cosmic Ray (CR) particles at shocks, their advection and their dynamical feedback on thermal baryonic gas. We give a description of the algorithms and show their tests against analytical and idealized one-dimensional problems. Our implementation is able to track the injection of CR energy, the spatial advection of CR energy and its feedback on the thermal gas in run-time. This method is applied to study CR acceleration and evolution in cosmological volumes, with both fixed and variable mesh resolution. We compare the properties of galaxy clusters with and without CRs, for a sample of high-resolution clusters with different dynamical states. At variance with similar simulations based on Smoothed Particles Hydrodynamics, we report that the inclusion of CR feedback in our method decreases the central gas density in clusters, thus reducing the X-ray and Sunyaev-Zeldovich effect from the clusters centre, while enhancing the gas density and its related observables near the virial radius.

We present the numerical implementation of run-time injection of Cosmic Rays energy, their spatial advection and their dynamical feedack on baryonic gas in the cosmological grid code ENZO. We discuss the results of its application to large-scale simulations showing that the CR energy inside clusters of galaxies is small compared to the gas energy (less than a few percent), while the ratio is larger near the accretion regions of clusters and filaments (about 0.1-0.3). CR feedback has a small, but significant impact on the X-ray emission and Sunyaev-Zeldovich effect from clusters.

In this paper we address the question why cluster radio relics that are connected to shock acceleration, so-called radio gischt, have preferentially been found in the outskirts of galaxy clusters. By identifying merger shock waves in cosmological grid simulations, we explore several prescriptions for relating the energy dissipated in shocks to the energy emitted in the radio band. None of the investigated models produce detectable radio relics within 100-200 kpc from the cluster centre. All models cause > 50 per cent of the detectable relic emission at projected distances > 800 kpc. Central radio relics caused by shocks that propagate along the line-of-sight are rare events for simple geometrical reasons, and they have a low surface brightness making them elusive for current instruments. Our simulations show that the radial distribution of observed relics can be explained by the radial trend of dissipated kinetic energy in shocks, that increases with distance from the cluster centre up until half of the virial radius.

[...] We present results for the statistics of thermal gas and the shock wave properties for a large volume simulated with three different cosmological numerical codes: the Eulerian total variations diminishing code TVD, the Eulerian piecewise parabolic method-based code ENZO, and the Lagrangian smoothed-particle hydrodynamics code GADGET. Starting from a shared set of initial conditions, we present convergence tests for a cosmological volume of side-length 100 Mpc/h, studying in detail the morphological and statistical properties of the thermal gas as a function of mass and spatial resolution in all codes. By applying shock finding methods to each code, we measure the statistics of shock waves and the related cosmic ray acceleration efficiencies, within the sample of simulations and for the results of the different approaches. We discuss the regimes of uncertainties and disagreement among codes, with a particular focus on the results at the scale of galaxy clusters. Even if the bulk of thermal and shock properties are reasonably in agreement among the three codes, yet some significant differences exist (especially between Eulerian methods and smoothed particle hydrodynamics). In particular, we report: a) differences of huge factors (10-100) in the values of average gas density, temperature, entropy, Mach number and shock thermal energy flux in the most rarefied regions of the simulations between grid and SPH methods; b) the hint of an entropy core inside clusters simulated in grid codes; c) significantly different phase diagrams of shocked cells in grid codes compared to SPH; d) sizable differences in the morphologies of accretion shocks between grid and SPH methods.

Aims: We discuss the dependence of shocks, cosmic rays acceleration and turbulence on the dynamical state of the host clusters. Method: We perform cosmological simulations with the grid code ENZO 1.5, with a mesh refinement scheme tailored to follow at high resolution shocks and turbulence developed in the clusters volume. Results: Sizable differences are found when some important properties, connected to non-thermal activity in clusters, are compared for post-merger, merging and relaxing systems.

We study the properties of chaotic motions in the intra cluster medium using a set of 20 galaxy clusters simulated with large dynamical range, using the Adaptive Mesh Refinement code ENZO (e.g. Norman et al.2007). The adopted setup allows us to study the spectral and spatial properties of turbulent motions in galaxy clusters with unprecedented detail, achieving an maximum available Reynolds number of the order of R=500-1000 for the largest eddies. The correlations between the energy of these motions in the Intra Cluster Medium and the dynamical state of the host systems are studied, and the statistical properties of turbulent motions and their evolution with time support that major merger events are responsible for the injection of the bulk of turbulent kinetic energy inside cluster. Turbulence is found to account for a 20-30 per cent of the thermal energy in merging clusters, while it accounts for a 5 per cent in relaxed clusters. A comparison of the energies of turbulence and motions in our simulated clusters with present upper-limits in real nearby clusters, recently derived with XMM-Newton (Sanders et al.2010), is provided. When the same spatial scales of turbulent motions are compared, the data from simulations result well within the range presently allowed by observations. Finally, we comment on the possibility that turbulence may accelerate relativistic particles leading to the formation of giant radio halos in turbulent (merging) clusters. Based on our simulations we confirm previous semi-analytical studies that suggest that the fraction of turbulent clusters is consistent with that of clusters hosting radio halos.

We study the azimuthal scatter in the radial profiles of X-ray luminous galaxy clusters, with two sets of high-resolution cosmological re-simulations obtained with the codes ENZO and GADGET2. The average gas profiles are computed for different angular sectors of the cluster projected volume, and compared with the mean cluster profiles at each radius from the center. We report that in general the level of azimuthal scatter is found to be about 10 per cent for gas density, temperature and entropy inside R200, and about 25 per cent for X-ray luminosity for the same volume. These values generally doubles going to 2 R200 from the cluster center, and are generally found to be higher (by about 20-40 percent) in the case of perturbed systems. A comparison with results from recent SUZAKU observations is discussed, showing the possibility to simply interpret the large azimuthal scatter of observables in light of our simulated results.

We investigated the numerical and physical reasons leading to a flat distribution of low gas entropy in the core region of galaxy clusters, as commonly found in grid cosmological simulations. To this end, we run a set of 30 high resolution re-simulations of a 3 x 10^14 M_sol/h cluster of galaxies with the AMR code ENZO, exploring and investigating the details involved in the production of entropy in simulated galaxy clusters. The occurrence of the flat entropy core is found to be mainly due to hydro-dynamical processes resolved in the code and that additional spurious effects of numerical origin (e.g. artificial heating due to softening effects or N-body noise) can affect the size and level of the entropy core only in a minor way. We show that the entropy profile of non-radiative simulations is produced by a mechanism of "sorting in entropy" which takes place with regularity during the cluster evolution. Using gas tracers we prove that the flat entropy core is caused by physical mixing of gravity-driven subsonic motions within the shallow inner cluster potential. Re-simulations were also produced for the same cluster object with the addition of radiative cooling, uniform pre-heating at high redshift (z=10) and late (z<1) thermal energy feedback from AGN activity in the cluster, in order to assess the effects of such mechanisms on the final entropy profile of the cluster. We report on the infeasibility of balancing the catastrophic cooling and recovering a flat entropy profile with the investigated trials for AGN activity alone, while for a sub-set of pre-heating models, or AGN feedback plus pre-heating models, a flat entropy distribution similar to non-radiative runs can be obtained with a viable energy requirement, and in good consistency with X-ray observations.

On May 27th 2010, the Italian astronomical community learned with concern that the National Institute for Astrophysics (INAF) was going to be suppressed, and that its employees were going to be transferred to the National Research Council (CNR). It was not clear if this applied to all employees (i.e. also to researchers hired on short-term contracts), and how this was going to happen in practice. In this letter, we give a brief historical overview of INAF and present a short chronicle of the few eventful days that followed. Starting from this example, we then comment on the current situation and prospects of astronomical research in Italy.

We present a sample of 20 massive galaxy clusters with total virial masses in the range of 6 10^14 M_sol<M(vir)< 2 10^15M_sol, re-simulated with a customized version of the 1.5. ENZO code employing Adaptive Mesh Refinement. This technique allowed us to obtain unprecedented high spatial resolution (25kpc/h) up to the distance of 3 virial radii from the clusters center, and makes it possible to focus with the same level of detail on the physical properties of the innermost and of the outermost cluster regions, providing new clues on the role of shock waves and turbulent motions in the ICM, across a wide range of scales. In this paper, a first exploratory study of this data set is presented. We report on the thermal properties of galaxy clusters at z=0. Integrated and morphological properties of gas density, gas temperature, gas entropy and baryon fraction distributions are discussed, and compared with existing outcomes both from the observational and from the numerical literature. Our cluster sample shows an overall good consistency with the results obtained adopting other numerical techniques (e.g. Smoothed Particles Hydrodynamics), yet it provides a more accurate representation of the accretion patterns far outside the cluster cores. We also reconstruct the properties of shock waves within the sample by means of a velocity-based approach, and we study Mach numbers and energy distributions for the various dynamical states in clusters, giving estimates for the injection of Cosmic Rays particles at shocks. The present sample is rather unique in the panorama of cosmological simulations of massive galaxy clusters, due to its dynamical range, statistics of objects and number of time outputs. For this reason, we deploy a public repository of the available data, accessible via web portal at http://data.cineca.it.

We present a study of the mixing properties of the simulated intra cluster Medium, using tracers particles that are advected by the gas flow during the evolution of cosmic structures. Using a sample of seven galaxy clusters (with masses in the range of M=2-3 10^14Msol/h) simulated with a peak resolution of 25kpc/h up to the distance of two virial radii from their centers, we investigate the application of tracers to some important problems concerning the mixing of the ICM. The transport properties of the evolving ICM are studied through the analysis of pair dispersion statistics and mixing distributions. As an application, we focus on the transport of metals in the ICM. We adopt simple scenarios for the injection of metal tracers in the ICM, and find remarkable differences of metallicity profiles in relaxed and merger systems, also through the analysis of simulated emission from Doppler-shifted Fe XXIII lines.

We have implemented an Adaptive Mesh Refinement criterion explicitly designed to increase spatial resolution around discontinuities in the velocity field in ENZO cosmological simulations. With this technique, shocks and turbulent eddies developed during the hierarchical assembly of galaxy clusters are followed with unprecedented spatial resolution, even at large distances from the clusters center. By measuring the spectral properties of the gas velocity field, its time evolution and the properties of shocks for a reference galaxy cluster, we investigate the connection between accretion processes and the onset of chaotic motions in the simulated Inter Galactic Medium over a wide range of scales

We analyze the properties of Large Scale Shocks in a cosmological volume of size 103Mpc/h simulated with the public 1.0.1 release of the ENZO code. Different methods to identify and characterize shocks in post processing are discussed together with their uncertainties. Re-ionization affects the properties of shocks in simulations, and we propose a fitting procedure to model accurately the effect of re-ionization in non--radiative simulations, with a post--processing procedure. We investigate the properties of shocks by means of a procedure which uses jumps in the velocity variables across the cells in the simulations and this allows us to have a viable description of shocks also in relatively under-dense cosmic regions. We derive the distributions of the number of shocks and of the energy dissipated at these shocks as a function of their Mach number, and discuss the evolution of these distributions with cosmological time and across different cosmic environments (clusters, outskirts, filaments, voids). In line with some previous numerical studies (e.g. Ryu et al.2003, Pfrommer et al.2006) relatively weak shocks are found to dominate the process of energy dissipation in the simulated cosmic volume, although we find a larger ratio between weak and strong shocks with respect to previous studies. We estimate the rate of injection of CR at Large Scale Shocks by adopting injection efficiencies taken from previous numerical calculations. The flux dissipated in the form of CR within the whole simulated volume at the present epoch is about 20 per cent of the thermal energy dissipated at shocks; this fraction is smaller within galaxy clusters.

Aims: morpholgies, number and energy distributions of Cosmological Shock Waves from a set of ENZO cosmological simulations are produced, along with a study of the connection with Cosmic Rays processes in different environments. Method: we perform cosmological simulations with the public release of the PPM code ENZO, adopt a simple and physically motivated numerical setup to follow the evolution of cosmic structures at the resolution of 125kpc per cell, and characterise shocks with a new post processing scheme. Results: we estimate the efficency of the acceleration of Cosmic Ray particles and present the first comparison of our results with existing limits from observations of galaxy clusters.

We present a study of the turbulent velocity fields in the Intra Cluster Medium of a sample of 21 galaxy clusters simulated by the SPH--code Gadget2, using a new numerical scheme where the artificial viscosity is suppressed outside shocks. The turbulent motions in the ICM of our simulated clusters are detected with a novel method devised to better disentangle laminar bulk motions from chaotic ones. We focus on the scaling law between the turbulent energy content of the gas particles and the total mass, and find that the energy in the form of turbulence scales approximatively with the thermal energy of clusters. We follow the evolution with time of the scaling laws and discuss the physical origin of the observed trends. The simulated data are in agreement with independent semi--analytical calculations, and the combination between the two methods allows to constrain the scaling law over more than two decades in cluster mass.

Smoothed particle hydrodynamics (SPH) employs an artificial viscosity to properly capture hydrodynamical shock waves. In its original formulation, the resulting numerical viscosity is large enough to suppress structure in the velocity field on scales well above the nominal resolution limit, and to damp the generation of turbulence by fluid instabilities. This could artificially suppress random gas motions in the intracluster medium (ICM), which are driven by infalling structures during the hierarchical structure formation process. We show that this is indeed the case by analysing results obtained with an SPH formulation where an individual, time-variable viscosity is used for each particle (Monaghan 1997). Using test calculations involving strong shocks, we demonstrate that this scheme captures shocks as well as the original formulation of SPH, but, in regions away from shocks, the numerical viscosity is much smaller. In a set of nine high-resolution simulations of cosmological galaxy cluster formation, we find that this low--viscosity formulation of SPH produces substantially higher levels of turbulent gas motions in the ICM, reaching a kinetic energy content in random gas motions (measured within a 1Mpc cube) of up to 5%-30% of the thermal energy content, depending on cluster mass. This has also significant effects on radial gas profile. We find a central flattening of the entropy profile and a reduction of the central gas density in the low--viscosity scheme. Interestingly, this tends to reduce the differences seen in SPH and adaptive mesh refinement simulations of cluster formation. Finally, invoking a model for particle acceleration by MHD waves driven by turbulence, we find efficient electron acceleration to power diffuse radio emission.