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.Comment: 24 pages, 20 figures. MNRAS accepted. A version of the paper with higher quality figures can be found at this url: http://www.ira.inaf.it/~vazza/papers/feedback_vazza.pd

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.Comment: MNRAS accepted, in press. 18 pages, 18 Figures. New version to match with the one published in MNRAS. Updated publication list and footnote added to the title as obituary notic

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.Comment: 6 pages, 4 figures. MNRAS accepte

Convolutional deep denoising autoencoders for radio astronomical images

We apply a Machine Learning technique known as Convolutional Denoising Autoencoder to denoise synthetic images of state-of-the-art radio telescopes, with the goal of detecting the faint, diffused radio sources predicted to characterize the radio cosmic web. In our application, denoising is intended to address both the reduction of random instrumental noise and the minimization of additional spurious artefacts like the sidelobes, resulting from the aperture synthesis technique. The effectiveness and the accuracy of the method are analysed for different kinds of corrupted input images, together with its computational perfoance. Specific attention has been devoted to create realistic mock observations for the training, exploiting the outcomes of cosmological numerical silations, to generate images corresponding to LOFAR HBA 8 h observations at 150 MHz. Our autoencoder can effectively denoise complex images identifying and extracting faint objects at the limits of the instrumental sensitivity. The method can efficiently scale on large data sets, exploiting high-perfoance computing solutions, in a fully automated way (i.e. no human supervision is required after training). It can accurately perfo image segmentation, identifying low brightness outskirts of diffused sources, proving to be a viable solution for detecting challenging extended objects hidden in noisy radio observations

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.Comment: 17 pages, 11 figures. MNRAS accepte

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.Comment: 19 pages, 20 figures. Astronomy and Astrophysics, in pres

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.Comment: 23 pages, 18 figures, accepted for publication by Astronomy and Computin

Detecting shocked intergalactic gas with X-ray and radio observations

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

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.Comment: 19 pages, 18 figures. A&A accepted, in press. The public repository of radio maps for the full volumes studied in this work is available at http://www.hs.uni-hamburg.de/DE/Ins/Per/Vazza/projects/Public_data.htm