Source author record

Franco Vazza

Franco 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

37works
6topics
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

37 published item(s)

preprint2022arXiv

Filamentary Baryons and Where to Find Them: A forecast of synchrotron radiation from merger and accretion shocks in the local Cosmic Web

We generate probabilistic predictions of the low-redshift ($z < 0.2$) synchrotron Cosmic Web for half of the Northern Sky. In particular, we predict the contribution to the specific intensity function at $ν_\mathrm{obs} = 150\ \mathrm{MHz}$ from merger shocks in clusters and accretion shocks in filaments, both of which arise during large-scale structure formation. We assume a primordial magnetogenesis scenario, but our method is general enough to allow for an exploration of alternative magnetogenesis scenarios - and even alternative radiation mechanisms and spectral windows. In the future, by comparing different predictions, one could infer the most plausible physical model from a synchrotron Cosmic Web detection. Our method combines Bayesian large-scale structure reconstructions, snapshots of an MHD cosmological simulation, a Gaussian random field, and a ray tracing approach. The results help to select targets for deep observations and can be used in actual detection experiments. We highlight predictions for the Hercules Cluster, the Coma Cluster, Abell 2199, Abell 2255, the Lockman Hole and the Ursa Major Supercluster. At degree-scale resolution, the median specific intensity reaches $m_{I_ν} \sim 10^{-1}\ ξ_e\ \mathrm{Jy\ deg^{-2}}$, where $ξ_e$ is the electron acceleration efficiency.

preprint2022arXiv

Stacking the Synchrotron Cosmic Web with FIGARO

Recently Vernstrom et al. (2021) claimed the first definitive detection of the synchrotron cosmic web, obtained by `stacking' hundreds of thousands of pairs of close-proximity clusters in low-frequency radio observations and looking for a residual excess signal spanning the intracluster bridge. A reproduction study by Hodgson et al. (2022), using both the original radio data as well as new observations with the Murchison Widefield Array, failed to confirm these findings. Whilst the detection remains unsure, we here turn to stacking a simulated radio sky to understand what kind of excess radio signal is predicted by our current best cosmological models of the synchrotron cosmic web. We use the FIlaments & GAlactic RadiO (Hodgson et al., 2021a) simulation, which models both the synchrotron cosmic web as well as various subtypes of active galactic nucleii and star forming galaxies. Being a simulation, we have perfect knowledge of the location of clusters and galaxy groups which we use in our own stacking experiment. Whilst we do find an excess radio signature in our stacks that is attributable to the synchrotron cosmic web, its distribution is very different to that found by Vernstrom et al. (2021). Instead, we observe the appearance of excess emission on the immediate interiors of cluster pairs as a result of asymmetric, `radio relic'-like shocks surrounding cluster cores, whilst the excess emission spanning the intracluster region -- attributable to filaments proper -- is two orders of magnitude lower and undetectable in our experiment even under ideal conditions.

preprint2022arXiv

The redshift evolution of extragalactic magnetic fields

Faraday rotation studies of distant radio sources can constrain the evolution and the origin of cosmic magnetism. We use data from the LOFAR Two Metre Sky Survey: Data Release 2 (LoTSS DR2) to study the dependence of the Faraday rotation measure (RM) on redshift. By focusing on radio sources that are close in terms of their projection on the sky, but physically unrelated (random pairs), we measure the RM difference, $Δ$RM, between the two sources. Thus, we isolate the extragalactic contribution to $Δ$RM from other contributions. We present a statistical analysis of the resulting sample of random pairs and find a median absolute RM difference |$Δ$RM| $ = (1.79 \pm 0.09)$ rad/m$^{2}$ , with |$Δ$RM| uncorrelated both with respect to the redshift difference of the pair and the redshift of the nearer source, and a median excess of random pairs over physical pairs of $(1.65 \pm 0.10)$ rad/m$^{2}$. We seek to reproduce this result with Monte Carlo simulations assuming a non vanishing seed cosmological magnetic field and a redshift evolution of the comoving magnetic field strength that varies as $1/(1 + z)^γ$. We find the best fitting results $B_0 \equiv B_{\rm comoving}(z = 0) \lesssim (2.0 \pm 0.2)$ nG and $γ\lesssim 4.5 \pm 0.2$ that we conservatively quote as upper limits due to an unmodelled but non vanishing contribution of local environments to the RM difference. A comparison with cosmological simulations shows our results to be incompatible with primordial magnetogenesis scenarios with uniform seed fields of order nG.

preprint2021arXiv

Exploring the relation between turbulent velocity and density fluctuations in the stratified intracluster medium

The dynamics of the intracluster medium (ICM) is affected by turbulence driven by several processes, such as mergers, accretion and feedback from active galactic nuclei. X-ray surface brightness fluctuations have been used to constrain turbulence in galaxy clusters. Here, we use simulations to further investigate the relation between gas density and turbulent velocity fluctuations, with a focus on the effect of the stratification of the ICM. In this work, we studied the turbulence driven by hierarchical accretion by analysing a sample of galaxy clusters simulated with the cosmological code ENZO. We used a fixed scale filtering approach to disentangle laminar from turbulent flows. In dynamically perturbed galaxy clusters, we found a relation between the root mean square of density and velocity fluctuations, albeit with a different slope than previously reported. The Richardson number is a parameter that represents the ratio between turbulence and buoyancy, and we found that this variable has a strong dependence on the filtering scale. However, we could not detect any strong relation between the Richardson number and the logarithmic density fluctuations, in contrast to results by recent and more idealised simulations. In particular, we find a strong effect from radial accretion, which appears to be the main driver for the gas fluctuations. The ubiquitous radial bias in the dynamics of the ICM suggests that homogeneity and isotropy are not always valid assumptions, even if the turbulent spectra follow Kolmogorov's scaling. Finally, we find that the slope of the velocity and density spectra are independent of cluster-centric radii.

preprint2021arXiv

FIGARO Simulation: FIlaments & GAlactic RadiO Simulation

We produce the first low to mid frequency radio simulation that incorporates both traditional extragalactic radio sources as well as synchrotron cosmic web emission. The FIlaments \& GAlactic RadiO (FIGARO) simulation includes ten unique \SI{4x4}{\degree} fields, incorporating active galactic nucleii (AGNs), star forming galaxies (SFGs) and synchrotron cosmic web emission out to a redshift of $z = 0.8$ and over the frequency range 100-1400 MHz. To do this, the simulation brings together a recent $100^3$ Mpc$^3$ magneto-hydrodynamic simulation (Vazza et al., 2019), calibrated to match observed radio relic population statistics, alongside updated `T-RECS' code for simulating extragalactic radio sources (Bonaldi et al., 2019). Uniquely, the AGNs and SFGs are populated and positioned in accordance with the underlying matter density of the cosmological simulation. In this way, the simulation provides an accurate understanding of the apparent morphology, angular scales, and brightness of the cosmic web as well as -- crucially -- the clustering properties of the cosmic web with respect to the embedded extragalactic radio population. We find that the synchrotron cosmic web does not closely trace the underlying mass distribution of the cosmic web, but is instead dominated by shocked shells of emission surrounding dark matter halos and resembles a large, undetected population of radio relics. We also show that, with accurate kernels, the cosmic web radio emission is clearly detectable by cross-correlation techniques and this signal is separable from the embedded extragalactic radio population. We offer the simulation as a public resource towards the development of techniques for detecting and measuring the synchrotron cosmic web.

preprint2021arXiv

Joint inference on the redshift distribution of fast radio burst and on the intergalactic baryon content

Context: Fast radio bursts are transient radio pulses of extragalactic origin. Their dispersion measure is indicative of the baryon content in the ionized intergalactic medium between the source and the observer. However, inference using unlocalized fast radio bursts is degenerate to the distribution of redshifts of host galaxies. Method: We perform a joint inference of the intergalactic baryon content and the fast radio burst redshift distribution with the use of Bayesian statistics by comparing the likelihood of different models to reproduce the observed statistics in order to infer the most likely models. In addition to two models of the intergalactic medium, we consider contributions from the local environment of the source, assumed to be a magnetar, as well as a representative ensemble of host and intervening galaxies. Results: Assuming that the missing baryons reside in the ionized intergalactic medium, our results suggest that the redshift distribution of observed fast radio bursts peaks at $z \lesssim 0.6$. However, conclusions from different instruments regarding the intergalactic baryon content diverge and thus require additional changes to the observed distribution of host redshifts, beyond those caused by telescope selection effects.

preprint2021arXiv

Ultra-Steep Spectrum Radio `Jellyfish' Uncovered in Abell 2877

We report on the discovery of a mysterious ultra-steep spectrum (USS) synchrotron source in the galaxy cluster Abell 2877. We have observed the source with the Murchison Widefield Array at five frequencies across 72-231 MHz and have found the source to exhibit strong spectral curvature over this range as well the steepest known spectra of a synchrotron cluster source, with a spectral index across the central three frequency bands of $α= -5.97^{+0.40}_{-0.48}$. Higher frequency radio observations, including a deep observation with the Australia Telescope Compact Array, fail to detect any of the extended diffuse emission. The source is approximately 370 kpc wide and bears an uncanny resemblance to a jellyfish with two peaks of emission and long tentacles descending south towards the cluster centre. Whilst the `USS Jellyfish' defies easy classification, we here propose that the phenomenon is caused by the reacceleration and compression of multiple aged electron populations from historic active galactic nucleus (AGN) activity, so-called `radio phoenix', by an as yet undetected weak cluster-scale mechanism. The USS Jellyfish adds to a growing number of radio phoenix in cool-core clusters with unknown reacceleration mechanisms; as the first example of a polyphoenix, however, this implies the mechanism is on the scale of the cluster itself. Indeed, we show that in simulations, emission akin to the USS Jellyfish can be produced as a short-lived, transient phase in the evolution of multiple interacting AGN remnants when subject to weak external shocks.

preprint2020arXiv

Discovering the most elusive radio relic in the sky: Diffuse Shock Acceleration caught in the act?

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 $α= 1.15\pm 0.23$. We infer for this radio relic a lower bound on the magnetisation of $B\geq 0.4\, μ$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.

preprint2020arXiv

Limiting the shock acceleration of cosmic-ray protons in the ICM

Observations of large-scale radio emissions prove the existence of shock accelerated cosmic-ray electrons in galaxy clusters, while the lack of detected $γ$-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.

preprint2020arXiv

Magnetism Science with the Square Kilometre Array

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.

preprint2020arXiv

Multi wavelength cross-correlation analysis of the simulated cosmic web

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.

preprint2020arXiv

Redshift estimates for fast radio bursts and implications on intergalactic magnetic fields

Context: Fast Radio Bursts are transient radio pulses from presumably compact stellar sources of extragalactic origin. With new telescopes detecting multiple events per day, statistical methods are required in order to interpret observations and make inferences regarding astrophysical and cosmological questions. Purpose: We present a method that uses probability estimates of fast radio burst observables to obtain likelihood estimates for the underlying models. Method: Considering models for all regions along the line-of-sight, including intervening galaxies, we perform Monte-Carlo simulations to estimate the distribution of the dispersion measure, rotation measure and temporal broadening. Using Bayesian statistics, we compare these predictions to observations of Fast Radio Bursts. Results: By applying Bayes theorem, we obtain lower limits on the redshift of Fast Radio Bursts with extragalactic DM $\gtrsim 400$ pc cm$^{-3}$. We find that intervening galaxies cannot account for all highly scattered Fast Radio Bursts in FRBcat, thus requiring a denser and more turbulent environment than a SGR 1935+2154-like magnetar. We show that a sample of $\gtrsim 10^3$ unlocalized Fast Radio Bursts with associated extragalactic RM $\geq 1$ rad m$^{-2}$ can improve current upper limits on the strength of intergalactic magnetic fields.

preprint2020arXiv

Shock waves in the magnetized cosmic web: the role of obliquity and cosmic-ray acceleration

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.

preprint2020arXiv

Spectro-polarimetric observations of the CIZA J2242.8+5301 northern radio relic: no evidence of high-frequency steepening

Observations of radio relics at very high frequency (>10 GHz) can help to understand how particles age and are (re-)accelerated in galaxy cluster outskirts and how magnetic fields are amplified in these environments. In this work, we present new single-dish 18.6 GHz Sardinia Radio Telescope and 14.25 GHz Effelsberg observations of the well known northern radio relic of CIZA J2242.8+5301. We detected the relic which shows a length of $\sim$1.8 Mpc and a flux density equal to $\rm S_{14.25\,GHz}=(9.5\pm3.9)\,mJy$ and $\rm S_{18.6\,GHz}=(7.67\pm0.90)\,mJy$ at 14.25 GHz and 18.6 GHz respectively. The resulting best-fit model of the relic spectrum from 145 MHz to 18.6 GHz is a power-law spectrum with spectral index $α=1.12\pm0.03$: no evidence of steepening has been found in the new data presented in this work. For the first time, polarisation properties have been derived at 18.6 GHz, revealing an averaged polarisation fraction of $\sim40\%$ and a magnetic field aligned with the 'filaments' or 'sheets' of the relic.

preprint2020arXiv

The Northern Cross Fast Radio Burst project. I. Overview and pilot observations at 408 MHz

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 \, (τ/{\rm ms})^{-0.5}$~mJy rms sensitivity.

preprint2019arXiv

Polarisation of Radio Relics in Galaxy Clusters

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.

preprint2016arXiv

A shock front at the radio relic of Abell 2744

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.

preprint2016arXiv

ALMA-SZ Detection of a Galaxy Cluster Merger Shock at Half the Age of the Universe

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 ~μ$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.

preprint2016arXiv

Evolution of cosmic filaments and of their galaxy population from MHD cosmological simulations

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.

preprint2016arXiv

Galaxy cluster outskirts from the thermal SZ and non-thermal synchrotron link

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.

preprint2016arXiv

On the complexity and the information content of cosmic structures

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.

preprint2016arXiv

On the non-thermal energy content of cosmic structures

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.

preprint2016arXiv

Propagation of Ultra High Energy Cosmic Rays in Extragalactic Magnetic Fields: A view from cosmological simulations

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.

preprint2016arXiv

Studying the Effect of Shock Obliquity on the $γ$-ray and Diffuse Radio Emission in Galaxy Clusters

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 $γ$-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 $γ$-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 $γ$-ray emission by the \textit{Fermi}-satellite.

preprint2016arXiv

The impact of the SZ effect on cm-wavelength (1-30 GHz) observation of galaxy cluster radio relics

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

preprint2015arXiv

Filaments of the radio cosmic web: opportunities and challenges for SKA

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

preprint2015arXiv

Properties of Cosmological Filaments extracted from Eulerian Simulations

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.

preprint2015arXiv

The X-ray/SZ view of the virial region. I. Thermodynamic properties

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.

preprint2015arXiv

The X-ray/SZ view of the virial region. II. Gas mass fraction

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.

preprint2014arXiv

Do radio relics challenge diffusive shock acceleration?

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

preprint2014arXiv

Non-thermal emission from galaxy clusters: feasibility study with SKA1

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.

preprint2014arXiv

Numerical cosmology on the GPU with Enzo and Ramses

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.

preprint2014arXiv

Simulations of cosmic rays in large-scale structures: numerical and physical effects

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.

preprint2011arXiv

Merger and non-merger galaxy clusters in cosmological AMR simulations

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.

preprint2010arXiv

A decline and fall in the future of Italian Astronomy?

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.

preprint2010arXiv

Massive and refined: a sample of large galaxy clusters simulated at high resolution. I:Thermal gas and shock waves properties

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.