Showing papers by "Klaus Dolag published in 2020"

First Simulations of Axion Minicluster Halos.

Abstract: We study the gravitational collapse of axion dark matter fluctuations in the postinflationary scenario, so-called axion miniclusters, with $N$-body simulations. Largely confirming theoretical expectations, overdensities begin to collapse in the radiation-dominated epoch and form an early distribution of miniclusters with masses up to ${10}^{\ensuremath{-}12}\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$. After matter-radiation equality, ongoing mergers give rise to a steep power-law distribution of minicluster halo masses. The density profiles of well-resolved halos are Navarro--Frenk--White-like to good approximation. The fraction of axion dark matter in these bound structures is $\ensuremath{\sim}0.75$ at redshift $z=100$.

The Three Hundred Project: Correcting for the hydrostatic-equilibrium mass bias in X-ray and SZ surveys

Abstract: Accurate and precise measurement of the masses of galaxy clusters is key to deriving robust constraints on cosmological parameters. However, increasing evidence from observations confirms that X-ray masses obtained under the assumption of hydrostatic equilibrium might be underestimated, as previously predicted by cosmological simulations. We analyze more than 300 simulated massive clusters from the Three Hundred Project, and investigate the connection between mass bias and several diagnostics extracted from synthetic X-ray images of these simulated clusters. We find that the azimuthal scatter measured in 12 sectors of the X-ray flux maps is a statistically significant indication of the presence of an intrinsic (i.e., 3D) clumpy gas distribution. We verify that a robust correction to the hydrostatic mass bias can be inferred when estimates of the gas inhomogeneity from X-ray maps (such as the azimuthal scatter or the gas ellipticity) are combined with the asymptotic external slope of the gas density or pressure profiles, which can be respectively derived from X-ray and millimeter (Sunyaev-Zeldovich effect) observations. We also obtain that mass measurements based on either gas density and temperature or gas density and pressure result in similar distributions of the mass bias. In both cases, we provide corrections that help reduce both the dispersion and skewness of the mass bias distribution. These are effective even when irregular clusters are included leading to interesting implications for the modeling and correction of hydrostatic mass bias in cosmological analyses of current and future X-ray and SZ cluster surveys.

On the impact of baryons on the halo mass function, bias, and cluster cosmology

Abstract: Luminous matter produces very energetic events, such as active galactic nuclei and supernova explosions, that significantly affect the internal regions of galaxy clusters. Although the current uncertainty in the effect of baryonic physics on cluster statistics is subdominant as compared to other systematics, the picture is likely to change soon as the amount of high-quality data is growing fast, urging the community to keep theoretical systematic uncertainties below the ever-growing statistical precision. In this paper, we study the effect of baryons on galaxy clusters, and their impact on the cosmological applications of clusters, using the Magneticum suite of cosmological hydrodynamical simulations. We show that the impact of baryons on the halo mass function can be recast in terms on a variation of the mass of the halos simulated with pure N-body, when baryonic effects are included. The halo mass function and halo bias are only indirectly affected. Finally, we demonstrate that neglecting baryonic effects on halos mass function and bias would significantly alter the inference of cosmological parameters from high-sensitivity next-generations surveys of galaxy clusters.

The DIANOGA simulations of galaxy clusters: characterizing star formation in proto-clusters

Abstract: We studied the star formation rate (SFR) in cosmological hydrodynamical simulations of galaxy (proto-)clusters in the redshift range $0

Encounters of merger and accretion shocks in galaxy clusters and their effects on intracluster medium

Abstract: Several types/classes of shocks naturally arise during formation and evolution of galaxy clusters. One such class is represented by accretion shocks, associated with deceleration of infalling baryons. Such shocks, characterized by a very high Mach number, are present even in 1D models of cluster evolution. Another class is composed of "runaway merger shocks", which appear when a merger shock, driven by a sufficiently massive infalling subcluster, propagates away from the main-cluster center. We argue that, when the merger shock overtakes the accretion shock, a new long-living shock is formed that propagates to large distances from the main cluster (well beyond its virial radius) affecting the cold gas around the cluster. We refer to these structures as Merger-accelerated Accretion shocks (MA-shocks) in this paper. We show examples of such MA-shocks in 1D and 3D simulations and discuss their characteristic properties. In particular, (1) MA-shocks shape the boundary separating the hot intracluster medium (ICM) from the unshocked gas, giving this boundary a "flower-like" morphology. In 3D, MA-shocks occupy space between the dense accreting filaments. (2) Evolution of MA-shocks highly depends on the Mach number of the runaway merger shock and the mass accretion rate parameter of the cluster. (3) MA-shocks may lead to the misalignment of the ICM boundary and the splashback radius.

The DIANOGA simulations of galaxy clusters: characterising star formation in protoclusters

Abstract: Aims. We studied the star formation rate (SFR) in cosmological hydrodynamical simulations of galaxy (proto-)clusters in the redshift range 0 We analyse a set of zoom-in cosmological hydrodynamical simulations centred on 12 clusters. The simulations are carried out with the GADGET-3 Tree-PM smoothed-particle hydro-dynamics code which includes various subgrid models to treat unresolved baryonic physics, including AGN feedback.Results. Simulations do not reproduce the high values of SFR observed within protocluster cores, where the values of SFR are underpredicted by a factor ≳4 both at z ∼ 2 and z ∼ 4. The difference arises as simulations are unable to reproduce the observed starburst population and is greater at z ∼ 2 because simulations underpredict the normalisation of the main sequence (MS) of star forming galaxies (i.e. the correlation between stellar mass and SFR) by a factor of ∼3. As the low normalisation of the MS seems to be driven by an underestimated gas fraction, it remains unclear whether numerical simulations miss starburst galaxies due to overly underpredicted gas fractions or overly low star formation efficiencies. Our results are stable against varying several parameters of the star formation subgrid model and do not depend on the details of AGN feedback.Conclusions. The subgrid model for star formation, introduced to reproduce the self-regulated evolution of quiescent galaxies, is not suitable to describe violent events like high-redshift starbursts. We find that this conclusion holds, independently of the parameter choice for the star formation and AGN models. The increasing number of multi-wavelength high-redshift observations will help to improve the current star formation model, which is needed to fully recover the observed star formation history of galaxy clusters.

Stellar property statistics of massive haloes from cosmological hydrodynamics simulations: common kernel shapes

Abstract: We study stellar property statistics, including satellite galaxy occupation, of massive halo populations realized by three cosmological hydrodynamics simulations: BAHAMAS + MACSIS, TNG300 of the IllustrisTNG suite, and Magneticum Pathfinder. The simulations incorporate independent sub-grid methods for astrophysical processes with spatial resolutions ranging from $1.5$ to $6$ kpc, and each generates samples of $1000$ or more halos with $M_{\rm halo}> 10^{13.5} M_{\odot}$ at redshift $z=0$. Applying localized, linear regression (LLR), we extract halo mass-conditioned statistics (normalizations, slopes, and intrinsic covariance) for a three-element stellar property vector consisting of: i) $N_{sat}$, the number of satellite galaxies with stellar mass, $M_{\star, \rm sat} > 10^{10} M_{\odot}$ within radius $R_{200c}$ of the halo; ii) $M_{\star,\rm tot}$, the total stellar mass within that radius, and; iii) $M_{\star,\rm BCG}$, the gravitationally-bound stellar mass of the central galaxy within a $100 \, \rm kpc$ radius. Scaling parameters for the three properties with halo mass show mild differences among the simulations, in part due to numerical resolution, but there is qualitative agreement on property correlations, with halos having smaller than average central galaxies tending to also have smaller total stellar mass and a larger number of satellite galaxies. Marginalizing over total halo mass, we find the satellite galaxy kernel, $p(\ln N_{sat}\,|\,M_{\rm halo},z)$ to be consistently skewed left, with skewness parameter $\gamma = -0.91 \pm 0.02$, while that of $\ln M_{\star,\rm tot}$ is closer to log-normal, in all three simulations. The highest resolution simulations find $\gamma \simeq -0.8$ for the $z=0$ shape of $p(\ln M_{\star,\rm BCG}\,|\,M_{\rm halo},z)$ and also that the fractional scatter in total stellar mass is below $10\%$ in halos more massive than $10^{14.3} M_{\odot}$.

Kinematics of simulated galaxies II: Probing the stellar kinematics of galaxies out to large radii

Abstract: We investigate the stellar kinematics of a sample of galaxies extracted from the hydrodynamic cosmological Magneticum Pathfinder simulations out to $5$ half-mass radii. We construct differential radial stellar spin profiles quantified by the observationally widely used $\lambda_\mathrm{R}$ and the closely related $(V/\sigma)$ parameters. We find three characteristic profile shapes: profiles exhibiting a (i) peak within $2.5$ half-mass radii and a subsequent decrease (ii) continuous increase that plateaus at larger radii typically with a high amplitude (iii) completely flat behaviour typically with low amplitude, in agreement with observations. This shows that the kinematic state of the stellar component can vary significantly with radius, suggesting a distinct interplay between in-situ star formation and ex-situ accretion of stars. Following the evolution of our sample through time, we provide evidence that the accretion history of galaxies with decreasing profiles is dominated by the anisotropic accretion of low mass satellites that get disrupted beyond $ \sim 2.0$ half-mass radii, building up a stellar halo with non-ordered motion while maintaining the central rotation already present at $z=2$. In fact, at $z=2$ decreasing profiles are the predominant profile class. Hence, we can predict a distinct formation pathway for galaxies with a decreasing profile and show that the centre resembles an old embedded disk. Furthermore, we show that the radius of the kinematic transition provides a good estimation for the transition radius from in-situ stars in the centre to accreted stars in the halo.

Cosmology dependence of halo masses and concentrations in hydrodynamic simulations

Abstract: We employ a set of Magneticum cosmological hydrodynamic simulations that span over $15$ different cosmologies, and extract masses and concentrations of all well-resolved haloes between $z=0-1$ for critical over-densities $\Delta_\texttt{vir}, \Delta_{200c}, \Delta_{500c}, \Delta_{2500c}$ and mean overdensity $\Delta_{200m}.$ We provide the first mass-concentration (Mc) relation and sparsity relation (i.e. $M_{\Delta1} - M_{\Delta2}$ mass conversion) of hydrodynamic simulations that is modelled by mass, redshift and cosmological parameters $\Omega_m, \Omega_b, \sigma_8, h_0$ as a tool for observational studies. We also quantify the impact that the Mc relation scatter and the assumption of NFW density profiles have on the uncertainty of the sparsity relation. We find that converting masses with the aid of a Mc relation carries an additional fractional scatter ($\approx 4\%$) originated from deviations from the assumed NFW density profile. For this reason we provide a direct mass-mass conversion relation fit that depends on redshift and cosmological parameters. We release the package hydro\_mc, a python tool that perform all kind of conversions presented in this paper.

The Cluster HEritage project with XMM-Newton: Mass Assembly and Thermodynamics at the Endpoint of structure formation. I. Programme overview

Abstract: The Cluster HEritage project with XMM-Newton - Mass Assembly and Thermodynamics at the Endpoint of structure formation (CHEX-MATE) is a three mega-second Multi-Year Heritage Programme to obtain X-ray observations of a minimally-biased, signal-to-noise limited sample of 118 galaxy clusters detected by Planck through the Sunyaev-Zeldovich effect. The programme, described in detail in this paper, aims to study the ultimate products of structure formation in time and mass. It is composed of a census of the most recent objects to have formed (Tier-1: 0.05 7.25 x 10e14 M_sun). The programme will yield an accurate vision of the statistical properties of the underlying population, measure how the gas properties are shaped by collapse into the dark matter halo, uncover the provenance of non-gravitational heating, and resolve the major uncertainties in mass determination that limit the use of clusters for cosmological parameter estimation. We will acquire X-ray exposures of uniform depth, designed to obtain individual mass measurements accurate to 15-20% under the hydrostatic assumption. We present the project motivations, describe the programme definition, and detail the ongoing multi-wavelength observational (lensing, SZ, radio) and theoretical effort that is being deployed in support of the project.

On the origin of magnetic driven winds and the structure of the galactic dynamo in isolated galaxies

Abstract: We investigate the build-up of the galactic dynamo and subsequently the origin of a magnetic driven outflow. We use a setup of an isolated disc galaxy with a realistic circum-galactic medium (CGM). We find good agreement of the galactic dynamo with theoretical and observational predictions from the radial and toroidal components of the magnetic field as function of radius and disc scale height. We find several field reversals indicating dipole structure at early times and quadrupole structure at late times. Together with the magnetic pitch angle and the dynamo control parameters Rα, Rω and D we present strong evidence for an α2-Ω dynamo. The formation of a bar in the centre leads to further amplification of the magnetic field via adiabatic compression which subsequently drives an outflow. Due to the Parker-Instability the magnetic field lines rise to the edge of the disc, break out and expand freely in the CGM driven by the magnetic pressure. Finally, we investigate the correlation between magnetic field and star formation rate. Globally, we find that the magnetic field is increasing as function of the star formation rate surface density with a slope between 0.3 and 0.45 in good agreement with predictions from theory and observations. Locally, we find that the magnetic field can decrease while star formation increases. We find that this effect is correlated with the diffusion of magnetic field from the spiral arms to the inter-arm regions which we explicitly include by solving the induction equation and accounting for non-linear terms.

Cosmology dependence of galaxy cluster scaling relations

Abstract: The abundance of galaxy clusters as a function of mass and redshift is a well known powerful cosmological probe, which relies on underlying modelling assumptions on the mass-observable relations (MOR). Some of the MOR parameters can be constrained directly from multi-wavelength observations, as the normalization at some reference cosmology, the mass-slope, the redshift evolution and the intrinsic scatter. However, the cosmology dependence of MORs cannot be tested with multi-wavelength observations alone. We use {\tt Magneticum} simulations to explore the cosmology dependence of galaxy cluster scaling relations. We run fifteen hydro-dynamical cosmological simulations varying $\Omega_m$, $\Omega_b$, $h_0$ and $\sigma_8$ (around a reference cosmological model). The MORs considered are gas mass, baryonic mass, gas temperature, $Y$ and velocity dispersion as a function of virial mass. We verify that the mass and redshift slopes and the intrinsic scatter of the MORs are nearly independent of cosmology with variations significantly smaller than current observational uncertainties. We show that the gas mass and baryonic mass sensitively depends only on the baryon fraction, velocity dispersion and gas temperature on $h_0$, and $Y$ on both baryon fraction and $h_0$. We investigate the cosmological implications of our MOR parameterization on a mock catalog created for an idealized eROSITA-like experiment. We show that our parametrization introduces a strong degeneracy between the cosmological parameters and the normalization of the MOR. Finally, the parameter constraints derived at different overdensity ($\Delta_{500c}$), for X-ray bolometric gas luminosity, and for different subgrid physics prescriptions are shown in the appendix.

Collision of merger and accretion shocks: formation of Mpc-scale contact discontinuity in the Perseus cluster

Abstract: Two Mpc-size contact discontinuities have recently been identified in the XMM-Newton and Suzaku X-ray observations of the outskirts of the Perseus cluster (Walker et al 2020) These structures have been tentatively interpreted as "sloshing cold fronts", which are customarily associated with differential motions of the cluster gas, perturbed by a merger In this study we consider an alternative scenario, namely, that the most prominent discontinuity near the cluster virial radius is the result of the collision between the accretion shock and the "runaway" merger shock We also discuss the possible origin of the second discontinuity at ~12 Mpc

Constraining the origin and models of chemical enrichment in galaxy clusters using the Athena X-IFU

Abstract: Chemical enrichment of the Universe at all scales is related to stellar winds and explosive supernovae phenomena. Metals produced by stars and later spread throughout the intracluster medium (ICM) at the megaparsec scale become a fossil record of the chemical enrichment of the Universe and of the dynamical and feedback mechanisms determining their circulation. As demonstrated by the results of the soft X-ray spectrometer onboard Hitomi, high-resolution X-ray spectroscopy is the path to differentiating among the models that consider different metal-production mechanisms, predict the outcoming yields, and are a function of the nature, mass, and/or initial metallicity of their stellar progenitor. Transformational results shall be achieved through improvements in the energy resolution and effective area of X-ray observatories, allowing them to detect rarer metals (e.g. Na, Al) and constrain yet-uncertain abundances (e.g. C, Ne, Ca, Ni). The X-ray Integral Field Unit (X-IFU) instrument onboard the next-generation European X-ray observatory Athena is expected to deliver such breakthroughs. Starting from 100 ks of synthetic observations of 12 abundance ratios in the ICM of four simulated clusters, we demonstrate that the X-IFU will be capable of recovering the input chemical enrichment models at both low (z = 0.1) and high (z = 1) redshifts, while statistically excluding more than 99.5% of all the other tested combinations of models. By fixing the enrichment models which provide the best fit to the simulated data, we also show that the X-IFU will constrain the slope of the stellar initial mass function within ∼12%. These constraints will be key ingredients in our understanding of the chemical enrichment of the Universe and its evolution.

Characterizing the Uncertainty in Cluster Magnetic Fields Derived from Rotation Measures

Abstract: Magnetic fields play vital roles in intracluster media (ICMs), but estimating their strengths and distributions from observations is a major challenge. Faraday rotation measures (RMs) are widely applied to this task, so it is critical to understand inherent uncertainties in RM analysis. In this paper, we seek to characterize those uncertainties given the types of information available today, independent of the specific technique used. We conduct synthetic RM observations through the ICM of a galaxy cluster drawn from an MHD cosmological simulation in which the magnetic field is known. We analyze the synthetic RM observations using an analytical formalism based on commonly used model assumptions allowing us to relate model physical variables to outcome uncertainties. Despite the simplicity of some assumptions, and unknown physical parameters, we are able to extract an approximate magnitude of the central magnetic field within an apparently irreducible uncertain factor approximately 3. Principal, largely irreducible, uncertainties come from the unknown depth along the line of sight of embedded polarized sources, the lack of robust coherence lengths from area-constrained polarization sampling, and the unknown scaling between ICM electron density and magnetic field strength. The RM-estimated central magnetic field strengths span more than an order of magnitude including the full range of synthetic experiments.

Constraining the origin and models of chemical enrichment in galaxy clusters using the Athena X-IFU

Abstract: The chemical enrichment of the Universe at all scales is related to stellar winds and explosive supernovae phenomena. Metals produced by stars and later spread at the mega-parsec scale through the intra-cluster medium (ICM) become a fossil record of the chemical enrichment of the Universe and of the dynamical and feedback mechanisms determining their circulation. As demonstrated by the results of the soft X-ray spectrometer onboard Hitomi, high resolution X-ray spectroscopy is the path to to differentiate among the models that consider different metal production mechanisms, predict the outcoming yields, and are function of the nature, mass, and/or initial metallicity of their stellar progenitor. Transformational results shall be achieved through improvements in the energy resolution and effective area of X-ray observatories to detect rare metals (e.g. Na, Al) and constrain yet uncertain abundances (e.g. C, Ne, Ca, Ni). The X-ray Integral Field Unit (X-IFU) instrument onboard the next-generation European X-ray observatory Athena is expected to deliver such breakthroughs. Starting from 100 ks of synthetic observations of 12 abundance ratios in the ICM of four simulated clusters, we demonstrate that the X-IFU will be capable of recovering the input chemical enrichment models at both low ($z = 0.1$) and high ($z = 1$) redshifts, while statistically excluding more than 99.5% of all the other tested combinations of models. By fixing the enrichment models which provide the best fit to the simulated data, we also show that the X-IFU will constrain the slope of the stellar initial mass function within $\sim$12%. These constraints will be key ingredients in our understanding of the chemical enrichment of the Universe and its evolution.

Impact of baryons in cosmic shear analyses with tomographic aperture mass statistics

Abstract: NonGaussian cosmic shear statistics based on weak-lensing aperture mass ($M_{\rm ap}$) maps can outperform the classical shear two-point correlation function ($\gamma$-2PCF) in terms of cosmological constraining power. However, reaching the full potential of these new estimators requires accurate modeling of the physics of baryons as the extra nonGaussian information mostly resides at small scales. We present one such modeling based on the Magneticum hydrodynamical simulation for the KiDS-450 and DES-Y1 surveys and a Euclid-like survey. We compute the bias due to baryons on the lensing PDF and the distribution of peaks and voids in $M_{\rm ap}$ maps and propagate it to the cosmological forecasts on the structure growth parameter $S_8$, the matter density parameter $\Omega_{\rm m}$, and the dark energy equation of state $w_0$ using the SLICS and cosmo-SLICS sets of dark-matter-only simulations. We report a negative bias of a few percent on $S_8$ and $\Omega_{\rm m}$ and also measure a positive bias of the same level on $w_0$ when including a tomographic decomposition. These biases reach $\sim 5$% when combining $M_{\rm ap}$ statistics with the $\gamma$-2PCF as these estimators show similar dependency on the AGN feedback. We verify that these biases constitute a less than $1\sigma$ shift on the probed cosmological parameters for current cosmic shear surveys. However, baryons need to be accounted for at the percentage level for future Stage IV surveys and we propose to include the uncertainty on the AGN feedback amplitude by marginalizing over this parameter using multiple simulations such as those presented in this paper. Finally, we explore the possibility of mitigating the impact of baryons by filtering the $M_{\rm ap}$ map but find that this process would require to suppress the small-scale information to a point where the constraints would no longer be competitive.

Characterizing the Uncertainty in Cluster Magnetic Fields derived from Rotation Measures

Abstract: Magnetic fields play vital roles in intracluster media (ICMs), but estimating their strengths and distributions from observations is a major challenge. Faraday rotation measures (RMs) are widely applied to this task, so it is critical to understand inherent uncertainties in RM analysis. In this paper, we seek to characterize those uncertainties given the types of information available today, independent of the specific technique used. We conduct synthetic RM observations through the ICM of a galaxy cluster drawn from an MHD cosmological simulation in which the magnetic field is known. We analyze the synthetic RM observations using an analytical formalism based on commonly used model assumptions allowing us to relate model physical variables to outcome uncertainties. Despite the simplicity of some assumptions, and unknown physical parameters, we are able to extract an approximate magnitude of the central magnetic field within an apparently irreducible uncertain factor approximately 3. Principal, largely irreducible, uncertainties come from the unknown depth along the line of sight of embedded polarized sources, the lack of robust coherence lengths from area-constrained polarization sampling, and the unknown scaling between ICM electron density and magnetic field strength. The RM-estimated central magnetic field strengths span more than an order of magnitude including the full range of synthetic experiments.

On the phase-space structure of galaxy clusters from cosmological simulations

Abstract: Cosmological N-body simulations represent an excellent tool to study the formation and evolution of dark matter (DM) halos and the mechanisms that have originated the universal profile at the largest mass scales in the Universe. In particular, the combination of the velocity dispersion $\sigma_\mathrm{v}$ with the density $\rho$ can be used to define the pseudo-entropy $S(r)=\sigma_\mathrm{v}^2/\rho^{\,2/3}$, whose profile is well-described by a simple power-law $S\propto\,r^{\,\alpha}$. We analyze a set of cosmological hydrodynamical re-simulations of massive galaxy clusters and study the pseudo-entropy profiles as traced by different collisionless components in simulated galaxy clusters: DM, stars, and substructures. We analyze four sets of simulations, exploring different resolution and physics (N-body and full hydrodynamical simulations) to investigate convergence and the impact of baryons. We find that baryons significantly affect the inner region of pseudo-entropy profiles as traced by substructures, while DM particles profiles are characterized by an almost universal behavior, thus suggesting that the level of pseudo-entropy could represent a potential low-scatter mass-proxy. We compare observed and simulated pseudo-entropy profiles and find good agreement in both normalization and slope. We demonstrate, however, that the method used to derive observed pseudo-entropy profiles could introduce biases and underestimate the impact of mergers. Finally, we investigate the pseudo-entropy traced by the stars focusing our interest in the dynamical distinction between intracluster light (ICL) and the stars bound to the brightest cluster galaxy (BCG): the combination of these two pseudo-entropy profiles is well-described by a single power-law out to almost the entire cluster virial radius.

Properties of the CGM and IGM: constraints on galaxy formation models from the Sunyaev-Zel'dovich effect

Abstract: We present a comparison of the physical properties of the gas in the circumgalactic (CGM) and intergalactic (IGM) media at $z\sim0$ between observations and four cosmological hydrodynamical simulations: Illustris, TNG300 of the IllustrisTNG project, EAGLE, and one of the Magneticum simulations. For the observational data, we use the gas properties that are inferred from cross-correlating the Sunyaev-Zel'dovich effect (SZE) from the {\it Planck} CMB maps with the haloes and the large-scale structure reconstructed from Sloan Digital Sky Survey data. Both the observational and simulation results indicate that the integrated gas pressure in haloes deviates from the self-similar case, showing that feedback impacts haloes with $M_{500}\sim 10^{12-13}\,{\rm M_\odot}$. The simulations predict that more than half the baryons are displaced from haloes, while the gas fraction inferred from our observational data roughly equals the cosmic baryon fraction throughout the $M_{500}\sim 10^{12-14.5}\,{\rm M_\odot}$ halo mass range. All simulations tested here predict that the mean gas temperature in haloes is about the virial temperature, while that inferred from the SZE is up to one order of magnitude lower than that from the simulations (and also from X-ray observations). While a remarkable agreement is found for the average properties of the IGM between the observation and some simulations, we show that their dependence on the large-scale tidal field can break the degeneracy between models that show similar predictions otherwise. Finally, we show that the gas pressure and the electron density profiles from simulations are not well described by a generalized NFW (GNFW) profile. Instead, we present a new model with a mass-dependent shape that fits the profiles accurately.

Gadget3 on GPUs with OpenACC.

Abstract: We present preliminary results of a GPU porting of all main Gadget3 modules (gravity computation, SPH density computation, SPH hydrodynamic force, and thermal conduction) using OpenACC directives. Here we assign one GPU to each MPI rank and exploit both the host and accellerator capabilities by overlapping computations on the CPUs and GPUs: while GPUs asynchronously compute interactions between particles within their MPI ranks, CPUs perform tree-walks and MPI communications of neighbouring particles. We profile various portions of the code to understand the origin of our speedup, where we find that a peak speedup is not achieved because of time-steps with few active particles. We run a hydrodynamic cosmological simulation from the Magneticum project, with $2\cdot10^{7}$ particles, where we find a final total speedup of $\approx 2.$ We also present the results of an encouraging scaling test of a preliminary gravity-only OpenACC porting, run in the context of the EuroHack17 event, where the prototype of the porting proved to keep a constant speedup up to $1024$ GPUs.

Gadget3 on GPUs with OpenACC

Abstract: We present preliminary results of a GPU porting of all main Gadget3 modules (gravity computation, SPH density computation, SPH hydrodynamic force, and thermal conduction) using OpenACC directives. Here we assign one GPU to each MPI rank and exploit both the host and accellerator capabilities by overlapping computations on the CPUs and GPUs: while GPUs asynchronously compute interactions between particles within their MPI ranks, CPUs perform tree-walks and MPI communications of neighbouring particles. We profile various portions of the code to understand the origin of our speedup, where we find that a peak speedup is not achieved because of time-steps with few active particles. We run a hydrodynamic cosmological simulation from the Magneticum project, with $2\cdot10^{7}$ particles, where we find a final total speedup of $\approx 2.$ We also present the results of an encouraging scaling test of a preliminary gravity-only OpenACC porting, run in the context of the EuroHack17 event, where the prototype of the porting proved to keep a constant speedup up to $1024$ GPUs.

Propagation of cosmic rays and their secondaries in the intracluster medium

Abstract: We present results of the propagation of high-energy cosmic rays (CRs) and their secondaries in the intracluster medium (ICM) To this end, we employ three-dimensional cosmological magnetohydrodynamical simulations of the turbulent intergalactic medium to explore the propagation of CRs with energies between 1014 and 1019 eV We study the interaction of test particles with this environment considering all relevant electromagnetic, photohadronic, photonuclear, and hadronuclear processes Finally, we discuss the consequences of the confinement of high-energy CRs in clusters for the production of gamma rays and neutrinos

Driving galactic winds with magnetic fields at low and high redshift

Abstract: Although playing a key role for our understanding of the evolution of galaxies, the exact way how observed galactic outflows are driven is still far from being understood and therefore our understanding of associated feedback mechanisms that control the evolution of galaxies is still plagued by many enigmas. In this work we present a simple toy model that can provide insight on how non-axis-symmetric instabilities in galaxies (bars, spiral-arms, warps) can lead to local exponential magnetic field growth by a radial flows beyond the equipartition value by at least two orders of magnitude on a time-scale of a few $100$ Myr. Our predictions show that the process can lead to galactic outflows in barred spiral galaxies with a mass loading factor $\eta \approx 0.1$, in agreement with our numerical simulations. Moreover, our outflow mechanism could contribute to an understanding of the large fraction of bared spiral galaxies that show signs of galactic outflows in the CHANG-ES survey. Extending our model shows the importance of such processes in high redshift galaxies by assuming equipartition between magnetic energy and turbulent energy. Simple estimates for the star formation rate (SFR) in our model together with cross-correlated masses from the star-forming main-sequence at redshifts $z\sim2$ allow us to estimate the outflow rate and mass loading factors by non-axis-symmetric instabilities and a subsequent radial inflow dynamo, giving mass loading factors of $\eta \approx 0.1$ for galaxies in the range of M$_{\star}=10^9 - 10^{12}$ M$_{\odot}$, in good agreement with recent results of Sinfoni and KMOS$^{3\mathrm{D}}$.