No Arabic abstract
Dark matter-only simulations are able to produce the cosmic structure of a $Lambda$CDM universe, at a much lower computational cost than more physically motivated hydrodynamical simulations. However, it is not clear how well smaller substructure is reproduced by dark matter-only simulations. To investigate this, we directly compare the substructure of galaxy clusters and of surrounding galaxy groups in hydrodynamical and dark matter-only simulations. We utilise TheThreeHundred project, a suite of 324 simulations of galaxy clusters that have been simulated with hydrodynamics, and in dark matter-only. We find that dark matter-only simulations underestimate the number density of galaxies in the centres of groups and clusters relative to hydrodynamical simulations, and that this effect is stronger in denser regions. We also look at the phase space of infalling galaxy groups, to show that dark matter-only simulations underpredict the number density of galaxies in the centres of these groups by about a factor of four. This implies that the structure and evolution of infalling groups may be different to that predicted by dark matter-only simulations. Finally, we discuss potential causes for this underestimation, considering both physical effects, and numerical differences in the analysis.
In the outer regions of a galaxy cluster, galaxies may be either falling into the cluster for the first time, or have already passed through the cluster centre at some point in their past. To investigate these two distinct populations, we utilise TheThreeHundred project, a suite of 324 hydrodynamical resimulations of galaxy clusters. In particular, we study the backsplash population of galaxies; those that have passed within $R_{200}$ of the cluster centre at some time in their history, but are now outside of this radius. We find that, on average, over half of all galaxies between $R_{200}$ and $2R_{200}$ from their host at $z=0$ are backsplash galaxies, but that this fraction is dependent on the dynamical state of a cluster, as dynamically relaxed clusters have a greater backsplash fraction. We also find that this population is mostly developed at recent times ($zleq0.4$), and is dependent on the recent history of a cluster. Finally, we show that the dynamical state of a given cluster, and thus the fraction of backsplash galaxies in its outskirts, can be predicted based on observational properties of the cluster.
Recent X-ray observations of galaxy clusters show that the distribution of intra-cluster medium (ICM) metallicity is remarkably uniform in space and time. In this paper, we analyse a large sample of simulated objects, from poor groups to rich clusters, to study the dependence of the metallicity and related quantities on the mass of the systems. The simulations are performed with an improved version of the Smoothed-Particle-Hydrodynamics texttt{GADGET-3} code and consider various astrophysical processes including radiative cooling, metal enrichment and feedback from stars and active galactic nuclei (AGN). The scaling between the metallicity and the temperature obtained in the simulations agrees well in trend and evolution with the observational results obtained from two data samples characterised by a wide range of masses and a large redshift coverage. We find that the iron abundance in the cluster core ($r<0.1R_{500}$) does not correlate with the temperature nor presents a significant evolution. The scale invariance is confirmed when the metallicity is related directly to the total mass. The slope of the best-fitting relations is shallow ($betasim-0.1$) in the innermost regions ($r<0.5R_{500}$) and consistent with zero outside. We investigate the impact of the AGN feedback and find that it plays a key role in producing a constant value of the outskirts metallicity from groups to clusters. This finding additionally supports the picture of early enrichment.
We analyse the gas content evolution of infalling haloes in cluster environments from THE THREE HUNDRED project, a collection of 324 numerically modelled galaxy clusters. The haloes in our sample were selected within $5R_{200}$ of the main cluster halo at $z=0$ and have total halo mass $M_{200}geq10^{11} h^{-1} M_{odot}$. We track their main progenitors and study their gas evolution since their crossing into the infall region, which we define as $1-4R_{200}$. Studying the radial trends of our populations using both the full phase space information and a line-of-sight projection, we confirm the Arthur et al. (2019) result and identify a characteristic radius around $1.7R_{200}$ in 3D and at $R_{200}$ in projection at which infalling haloes lose nearly all of the gas prior their infall. Splitting the trends by subhalo status we show that subhaloes residing in group-mass and low-mass host haloes in the infall region follow similar radial gas-loss trends as their hosts, whereas subhaloes of cluster-mass host haloes are stripped of their gas much further out. Our results show that infalling objects suffer significant gaseous disruption that correlates with time-since-infall, cluster-centric distance and host mass, and that the gaseous disruption they experience is a combination of subhalo pre-processing and object gas depletion at a radius which behaves like an accretion shock.
We address the issue of numerical convergence in cosmological smoothed particle hydrodynamics simulations using a suite of runs drawn from the EAGLE project. Our simulations adopt subgrid models that produce realistic galaxy populations at a fiducial mass and force resolution, but systematically vary the latter in order to study their impact on galaxy properties. We provide several analytic criteria that help guide the selection of gravitational softening for hydrodynamical simulations, and present results from runs that both adhere to and deviate from them. Unlike dark matter-only simulations, hydrodynamical simulations exhibit a strong sensitivity to gravitational softening, and care must be taken when selecting numerical parameters. Our results--which focus mainly on star formation histories, galaxy stellar mass functions and sizes--illuminate three main considerations. First, softening imposes a minimum resolved escape speed, $v_epsilon$, due to the binding energy between gas particles. Runs that adopt such small softening lengths that $v_epsilon gt 10,{rm km s^{-1}}$ (the sound speed in ionised $sim 10^4,{rm K}$ gas) suffer from reduced effects of photo-heating. Second, feedback from stars or active galactic nuclei may suffer from numerical over-cooling if the gravitational softening length is chosen below a critical value, $epsilon_{rm eFB}$. Third, we note that small softening lengths exacerbate the segregation of stars and dark matter particles in halo centres, often leading to the counter-intuitive result that galaxy sizes {em increase} as softening is reduced. The structure of dark matter haloes in hydrodynamical runs respond to softening in a way that reflects the sensitivity of their galaxy populations to numerical parameters.
We searched for isolated dark matter deprived galaxies within several state-of-the-art hydrodynamical simulations: Illustris, IllustrisTNG, EAGLE, and Horizon-AGN and found a handful of promising objects in all except Horizon-AGN. While our initial goal was to study their properties and evolution, we quickly noticed that all of them were located at the edge of their respective simulation boxes. After carefully investigating these objects using the full particle data, we concluded that they are not merely caused by a problem with the algorithm identifying bound structures. We provide strong evidence that these oddballs were created from regular galaxies that get torn apart due to unphysical processes when crossing the edge of the simulation box. We show that these objects are smoking guns indicating an issue with the implementation of the periodic boundary conditions of the particle data in Illustris, IllustrisTNG, and EAGLE, which was eventually traced down to be a minor bug occurring for a very rare set of conditions.