Source author record

F. Vazza

F. Vazza appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

42works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

42 published item(s)

preprint2025arXiv

On the interpretation of XRISM X-ray measurements of turbulence in the intracluster medium: a comparison with cosmological simulations

We investigate whether the properties of turbulent gas motions recently measured via X-ray spectroscopy in the Coma cluster of galaxies by XRISM are in tension with the turbulent picture established by current numerical cosmological simulations. We use a high-resolution simulation of a Coma-like cluster and show that the simulation yields velocity structure functions and X-ray line-widths that are compatible with those measured by the XRISM observations of Coma. In particular, it has been previously suggested that a much steeper turbulence spectrum than the Kolmogorov would be needed to explain the XRISM observations under a homogeneous, cluster volume-filling turbulence model. Our results show that this tension is overcome thanks to the more complicated turbulent picture in cosmological simulations, that indeed shows a patchy distribution of turbulent regions in galaxy clusters, with a spectrum that is generally consistent with a Kolmogorov power-law over a fairly wide range of scales. More generally, our study highlights the fact that the interpretation of XRISM data of galaxy clusters depends on the turbulence model used and the importance of combining data and advanced simulations in the future steps.

preprint2022arXiv

Deep low-frequency radio observations of Abell 2256 I: The filamentary radio relic

We present deep and high fidelity images of the merging galaxy cluster Abell 2256 at low frequencies, using the upgraded Giant Metrewave Radio Telescope (uGMRT) and LOw-Frequency ARray (LOFAR). This cluster hosts one of the most prominent known relics, with a remarkably spectacular network of filamentary substructures. The new uGMRT (300-850 MHz) and LOFAR (120-169 MHz) observations, combined with the archival Karl G. Jansky Very Large Array (VLA; 1-4 GHz) data, allowed us to carry out the first spatially resolved spectral analysis of the exceptional relic emission down to 6 arcsec resolution over a broad range of frequencies. Our new sensitive radio images confirm the presence of complex filaments of magnetized relativistic plasma also at low frequencies. We find that the integrated spectrum of the relic is consistent with a single power law, without any sign of spectral steepening, at least below 3 GHz. Unlike previous claims, the relic shows an integrated spectral index of $-1.07\pm0.02$ between 144 MHz and 3 GHz, which is consistent with the (quasi)stationary shock approximation. The spatially resolved spectral analysis suggests that the relic surface very likely traces the complex shock front, with a broad distribution of Mach numbers propagating through a turbulent and dynamically active intracluster medium. Our results show that the northern part of the relic is seen edge-on and the southern part close to face-on. We suggest that the complex filaments are regions where higher Mach numbers dominate the (re-)acceleration of electrons that are responsible for the observed radio emission.

preprint2022arXiv

Life cycle of cosmic-ray electrons in the intracluster medium

We simulate the evolution of relativistic electrons injected into the medium of a small galaxy cluster by a central radio galaxy, studying how the initial jet power affects the dispersal and the emission properties of radio plasma. By coupling passive tracer particles to adaptive-mesh cosmological MHD simulations, we study how cosmic-ray electrons are dispersed as a function of the input jet power. We also investigate how the latter affects the thermal and non-thermal properties of the intracluster medium, with differences discernible up to $\sim$ Gyr after the start of the jet. We evolved the energy spectra of cosmic-ray electrons, subject to energy losses that are dominated by synchrotron and inverse Compton emission as well as energy gains via re-acceleration by shock waves and turbulence. We find that in the absence of major mergers the amount of re-acceleration experienced by cosmic-ray electrons is not enough to produce long-lived detectable radio emissions. However, for all simulations the role of re-acceleration processes is crucial to maintain a significant and volume-filling reservoir of fossil electrons ($γ\sim 10^3$) for several Gyrs after the first injection by jets. This is important to possibly explain recent discoveries of cluster-wide emission and other radio phenomena in galaxy clusters.

preprint2022arXiv

Mapping "out-of-the-box" the properties of the baryons in massive halos

We study the distributions of the baryons in massive halos ($M_{vir} > 10^{13} \ h^{-1}M_{\odot}$) in the $Magneticum$ suite of Smoothed Particle Hydrodynamical cosmological simulations, out to the unprecedented radial extent of $10 R_{500,\mathrm c}$. We confirm that, under the action of non-gravitational physical phenomena, the baryon mass fraction is lower in the inner regions ($<R_{500,\mathrm c}$) of increasingly less massive halos, and rises moving outwards, with values that spans from 51% (87%) in the regions around $R_{500,\mathrm c}$ to 95% (100%) at $10R_{500,\mathrm c}$ of the cosmological value in the systems with the lowest (highest; $M_{vir} \sim 5 \times 10^{14} \ h^{-1}M_{\odot}$) masses. The galaxy groups almost match the gas (and baryon) fraction measured in the most massive halos only at very large radii ($r>6 R_{500,\mathrm c}$), where the baryon depletion factor $Y_{\rm bar} = f_{\rm bar} / (Ω_{\rm b}/Ω_{\rm m})$ approaches the value of unity, expected for "closed-box" systems. We find that both the radial and mass dependency of the baryon, gas, and hot depletion factors are predictable and follow a simple functional form. The star mass fraction is higher in less massive systems, decreases systematically with increasing radii, and reaches a constant value of $Y_{\rm star} \approx 0.09$, where also the gas metallicity is constant, regardless of the host halo mass, as a result of the early ($z>2$) enrichment process.

preprint2022arXiv

The Coma cluster at LOFAR frequencies II: the halo, relic, and a new accretion relic

We present LOw Frequency ARray observations of the Coma cluster field at 144\,MHz. The cluster hosts one of the most famous radio halos, a relic, and a low surface-brightness bridge. We detect new features that allow us to make a step forward in the understanding of particle acceleration in clusters. The radio halo extends for more than 2 Mpc, which is the largest extent ever reported. To the North-East of the cluster, beyond the Coma virial radius, we discover an arc-like radio source that could trace particles accelerated by an accretion shock. To the West of the halo, coincident with a shock detected in the X-rays, we confirm the presence of a radio front, with different spectral properties with respect to the rest of the halo. We detect a radial steepening of the radio halo spectral index between 144 MHz and 342 MHz, at $\sim 30^{\prime}$ from the cluster centre, that may indicate a non constant re-acceleration time throughout the volume. We also detect a mild steepening of the spectral index towards the cluster centre. For the first time, a radial change in the slope of the radio-X-ray correlation is found, and we show that such a change could indicate an increasing fraction of cosmic ray versus thermal energy density in the cluster outskirts. Finally, we investigate the origin of the emission between the relic and the source NGC 4789, and we argue that NGC4789 could have crossed the shock originating the radio emission visible between its tail and the relic.

preprint2022arXiv

The galaxy group NGC 507: newly detected AGN remnant plasma transported by sloshing

Jets from active galactic nuclei (AGN) are known to recurrently enrich their surrounding medium with mildly-relativistic particles and magnetic fields. Here, we present a detailed multi-frequency analysis of the nearby (z=0.01646) galaxy group NGC 507. In particular, we present new high-sensitivity and high spatial resolution radio images in the frequency range 144-675 MHz obtained using LOFAR and uGMRT observations. These reveal the presence of previously undetected diffuse radio emission with complex, filamentary morphology, likely related to a previous outburst of the central galaxy. Based on spectral ageing considerations, we derived that the plasma was first injected by the AGN 240-380 Myr ago and is now cooling. Our analysis of deep archival XMM-Newton data confirms that the system is dynamically disturbed, as previously suggested. We detect two discontinuities in the X-ray surface brightness distribution (in East and South direction) tracing a spiral pattern, which we interpret as cold fronts produced by sloshing motions. The remarkable spatial coincidence observed between the newly-detected arc-like radio filament and the southern concave X-ray discontinuity strongly suggests that the remnant plasma has been displaced by the sloshing motions on large scales. Overall, NGC 507 represents one of the clearest examples known to date in which a direct interaction between old AGN remnant plasma and the external medium is observed in a galaxy group. Our results are consistent with simulations, which suggest that filamentary emission can be created by the cluster/group weather disrupting AGN lobes and spreading their relativistic content into the surrounding medium.

preprint2021arXiv

The intra-cluster magnetic field in the double relic galaxy cluster Abell 2345

Magnetic fields are ubiquitous in galaxy clusters, yet their radial profile, power spectrum, and connection to host cluster properties are poorly known. Merging galaxy clusters hosting diffuse polarized emission in the form of radio relics offer a unique possibility to study the magnetic fields in these complex systems. In this paper, we investigate the intra-cluster magnetic field in Abell 2345. This cluster hosts two radio relics that we detected in polarization with 1-2 GHz JVLA observations. X-ray XMM-Newton images show a very disturbed morphology. We derived the Rotation Measure (RM) of five polarized sources within $\sim$ 1 Mpc from the cluster center applying the RM synthesis. Both, the average RM and the RM dispersion radial profiles probe the presence of intra-cluster magnetic fields. Using the thermal electron density profile derived from X-ray analysis and simulating a 3D magnetic field with fluctuations following a power spectrum derived from magneto-hydrodynamical cosmological simulations, we build mock RM images of the cluster. We constrained the magnetic field profile in the eastern radio relic sector by comparing simulated and observed RM images. We find that, within the framework of our model, the data require a magnetic field scaling with thermal electron density as $B(r)\propto n_e(r)$. The best model has a central magnetic field (within a 200 kpc radius) of $2.8\pm0.1$ $μ$G. The average magnetic field at the position of the eastern relic is $\sim$0.3 $μ$G, a factor 2.7 lower than the equipartition estimate.

preprint2020arXiv

Analytical model for cluster radio relics

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.

preprint2020arXiv

New constraints on the magnetization of the cosmic web using LOFAR Faraday rotation observations

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 ($Δ$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 $Δ$RM along the line of sight between non-physical pairs. From our analysis, we find no significant difference between the $Δ$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.

preprint2020arXiv

New mysteries and challenges from the Toothbrush relic: wideband observations from 550 MHz to 8 GHz

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.

preprint2020arXiv

Physical insights from the spectrum of the radio halo in MACS J0717.5+3745

We present new LOFAR observations of the massive merging galaxy cluster MACS J0717.5+3745. The cluster hosts the most powerful radio halo known to date. These new observations, in combination with published uGMRT (300$-$850 MHz) and VLA (1$-$6.5 GHz) data, reveal that the halo is more extended than previously thought, with a largest linear size of $\sim2.2 \rm Mpc$. The halo shows a steep spectrum ($α_{144\,\text{MHz}}^{1.5\,\text{GHz}}\sim-1.4$) and a steepening ($α_{1.5 \text{GHz}}^{5.5 \text{GHz}}\sim-1.9$) above 1.5 GHz. We find a strong scattering in spectral index maps on scales of 50$-$100 kpc. We suggest that such a strong scattering may be a consequence of the regime where inverse Compton dominate the energy losses of electrons. The spectral index becomes steeper and shows an increased curvature in the outermost regions of the halo. We combined the radio data with \textit{Chandra} observations to investigate the connection between the thermal and non-thermal components of the intracluster medium (ICM). Despite a significant substructure in the halo emission, the radio brightness correlates strongly with the X-ray brightness at all observed frequencies. The radio-versus-X-ray brightness correlation slope steepens at a higher radio frequency (from $b_{144 \text{MHz}}=0.67\pm0.05$ to $b_{3.0 \text{GHz}}=0.98\pm0.09$) and the spectral index shows a significant anti correlation with the X-ray brightness. Both pieces of evidence further support a spectral steepening in the external regions. The compelling evidence for a steep spectral index, the existence of a spectral break above 1.5 GHz, and the dependence of radio and X-ray surface brightness correlation on frequency are interpreted in the context of turbulent reacceleration models. Under this scenario, our results allowed us to constrain that the turbulent kinetic pressure of the ICM is up to 10%.

preprint2020arXiv

Second order Fermi reacceleration mechanisms and large scale synchrotron radio emission in intra-cluster bridges

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é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énic 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.

preprint2020arXiv

Turbulent pressure support and hydrostatic mass-bias in the intracluster medium

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}$.

preprint2020arXiv

Understanding the radio relic emission in the galaxy cluster MACS J0717.5+3745: spectral analysis

Radio relics are diffuse, extended synchrotron sources that originate from shock fronts generated during cluster mergers. The massive merging galaxy cluster MACS J0717.5+3745 hosts one of the more complex relics known to date. We present upgraded Giant Metrewave Radio Telescope band 3 (300-500 MHz) and band 4 (550-850 MHz) observations. These new observations, combined with published VLA and the new LOFAR HBA data, allow us to carry out a detailed, high spatial resolution spectral analysis of the relic over a broad range of frequencies. The integrated spectrum of the relic closely follows a power-law between 144 MHz and 5.5 GHz with a mean spectral slope $α=-1.16\pm0.03$. Despite its complex morphology, the subregions of the relic and the other isolated filaments also follow power-law behaviors, and show similar spectral slopes. Assuming Diffusive Shock Acceleration, we estimate a dominant Mach number of $\sim 3.7$ for the shocks that make up the relic. Comparison with recent numerical simulations suggests that in the case of radio relics, the slopes of the integrated radio spectra are determined by the Mach number of the accelerating shock, with $α$ nearly constant, namely between $-1.13$ and $-1.17$, for Mach numbers $3.5 - 4.0$. The spectral shapes inferred from spatially resolved regions show curvature, we speculate that the relic is inclined along the line-of-sight. The locus of points in the simulated color-color plots changes significantly with the relic viewing angle. We conclude that projection effects and inhomogeneities in the shock Mach number dominate the observed spectral properties of the relic in this complex system. Based on the new observations we raise the possibility that the relic and a narrow-angle-tailed radio galaxy are two different structures projected along the same line-of-sight.

preprint2019arXiv

Fundamental Physics with the Square Kilometre Array

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.

preprint2019arXiv

How complex is the cosmic web?

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.

preprint2016arXiv

Constraining the efficiency of cosmic ray acceleration by cluster shocks

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 γ-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.

preprint2016arXiv

Detecting the cosmic web with radio surveys

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.

preprint2016arXiv

Turbulence and Vorticity in Galaxy Clusters Generated by Structure Formation

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.

preprint2015arXiv

Constraining the electron and proton acceleration efficiencies in merger shocks in galaxy clusters

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 $γ$-rays in hadronic interactions. Combining observations of radio relics with upper limits from $γ$-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 $γ$-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 μ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}$.

preprint2015arXiv

Forecasts for the detection of the magnetised cosmic web from cosmological simulations

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.

preprint2015arXiv

Radio halos in SZ-selected clusters of galaxies: the making of a halo?

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.

preprint2015arXiv

The gas distribution in the outer regions of galaxy clusters

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]

preprint2015arXiv

Unravelling the origin of large-scale magnetic fields in galaxy clusters and beyond through Faraday Rotation Measures with the SKA

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.

preprint2014arXiv

A new double radio relic in PSZ1 G096.89+24.17 and a radio relic mass-luminosity relation

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.

preprint2014arXiv

Influence of adaptive mesh refinement and the hydro solver on shear-induced mass stripping in a minor-merger scenario

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.

preprint2014arXiv

On the amplification of magnetic fields in cosmic filaments and galaxy clusters

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.

preprint2013arXiv

Measurements and simulation of Faraday rotation across the Coma radio relic

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.

preprint2013arXiv

The Hot and Energetic Universe: The astrophysics of galaxy groups and clusters

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.

preprint2012arXiv

Discovery of radio halos and double-relics in distant MACS galaxy clusters: clues to the efficiency of particle acceleration

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).

preprint2012arXiv

Modelling injection and feedback of Cosmic Rays in grid-based cosmological simulations: effects on cluster outskirts

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.

preprint2012arXiv

Properties of gas clumps and gas clumping factor in the intra cluster medium

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.

preprint2012arXiv

Thermal and non-thermal traces of AGN feedback: results from cosmological AMR simulations

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.

preprint2012arXiv

Turbulence in the ICM from mergers, cool-core sloshing and jets: results from a new multi-scale filtering approach

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.

preprint2011arXiv

A Comparison of Cosmological Codes: Properties of Thermal Gas and Shock Waves in Large Scale Structures

[...] 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.

preprint2011arXiv

Massive and Refined. II. The statistical properties of turbulent motions in massive galaxy clusters with high spatial resolution

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.

preprint2011arXiv

The injection and feedback of Cosmic Rays in large-scale structures

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.

preprint2011arXiv

Why are central radio relics so rare?

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.

preprint2010arXiv

The entropy core in galaxy clusters: numerical and physical effects in cosmological grid simulations

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.

preprint2010arXiv

The Mixing and Transport Properties of the Intra Cluster Medium: a numerical study using tracers particles

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.

preprint2010arXiv

The scatter in the radial profiles of X-ray luminous galaxy clusters as diagnostic of the thermodynamical state of the ICM

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.

preprint2009arXiv

Turbulent Motions and Shocks Waves in Galaxy Clusters simulated with AMR

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