Using the latest cosmological hydrodynamic N-body simulations of groups and clusters, we study how location in phase-space coordinates at z=0 can provide information on environmental effects acting in clusters. We confirm the results of previous authors showing that galaxies tend to follow a typical path in phase-space as they settle into the cluster potential. As such, different regions of phase-space can be associated with different times since first infalling into the cluster. However, in addition, we see a clear trend between total mass loss due to cluster tides, and time since infall. Thus we find location in phase-space provides information on both infall time, and tidal mass loss. We find the predictive power of phase-space diagrams remains even when projected quantities are used (i.e. line-of-sight velocities, and projected distances from the cluster). We provide figures that can be directly compared with observed samples of cluster galaxies and we also provide the data used to make them as supplementary data, in order to encourage the use of phase-space diagrams as a tool to understand cluster environmental effects. We find that our results depend very weakly on galaxy mass or host mass, so the predictions in our phase-space diagrams can be applied to groups or clusters alike, or to galaxy populations from dwarfs up to giants.

Using high resolution hydrodynamical cosmological simulations, we conduct a comprehensive study of how tidal stripping removes dark matter and stars from galaxies. We find that dark matter is always stripped far more significantly than the stars -galaxies that lose ∼80% of their dark matter, typically lose only 10% of their stars. This is because the dark matter halo is initially much more extended than the stars. As such, we find the stellar-to-halo size-ratio (measured using r eff /r vir ) is a key parameter controlling the relative amounts of dark matter and stellar stripping. We use simple fitting formulae to measure the relation between the fraction of bound dark matter and fraction of bound stars. We measure a negligible dependence on cluster mass or galaxy mass. Therefore these formulae have general applicability in cosmological simulations, and are ideal to improve stellar stripping recipes in semi-analytical models, and/or to estimate the impact that tidal stripping would have on galaxies when only their halo mass evolution is known.

Hydrodynamical cosmological simulations are increasing their level of realism by considering more physical processes and having greater resolution or larger statistics. However, usually either the statistical power of such simulations or the resolution reached within galaxies are sacrificed. Here, we introduce the NEWHORIZON project in which we simulate at high resolution a zoom-in region of ∼(16 Mpc)3 that is larger than a standard zoom-in region around a single halo and is embedded in a larger box. A resolution of up to 34 pc, which is typical of individual zoom-in, up-to-date resimulated halos, is reached within galaxies; this allows the simulation to capture the multi-phase nature of the interstellar medium and the clumpy nature of the star formation process in galaxies. In this introductory paper, we present several key fundamental properties of galaxies and their black holes, including the galaxy mass function, cosmic star formation rate, galactic metallicities, the Kennicutt–Schmidt relation, the stellar-to-halo mass relation, galaxy sizes, stellar kinematics and morphology, gas content within galaxies and its kinematics, and the black hole mass and spin properties over time. The various scaling relations are broadly reproduced by NEWHORIZON with some differences with the standard observables. Owing to its exquisite spatial resolution, NEWHORIZON captures the inefficient process of star formation in galaxies, which evolve over time from being more turbulent, gas rich, and star bursting at high redshift. These high-redshift galaxies are also more compact, and they are more elliptical and clumpier until the level of internal gas turbulence decays enough to allow for the formation of discs. The NEWHORIZON simulation gives access to a broad range of galaxy formation and evolution physics at low-to-intermediate stellar masses, which is a regime that will become accessible in the near future through surveys such as the LSST.

We used the time since infall (TSI) of galaxies, obtained from the Yonsei Zoom-in Cluster Simulation, and the star formation rate (SFR) from the Sloan Digital Sky Survey (SDSS) Data Release 10 to study how quickly star formation of disk galaxies is quenched in cluster environments. We first confirm that both simulated and observed galaxies are consistently distributed in phase space. We then hypothesize that the TSI and SFR are causally connected; thus, both the TSI and SFR of galaxies at each position of phase space can be associated through abundance matching. Using a flexible model, we derive the star formation history (SFH) of cluster galaxies that best reproduces the relationship between the TSI and SFR at z ∼ 0.08. According to this SFH, we find that the galaxies with M * > 10 9.5 M generally follow the so-called "delayed-then-rapid" quenching pattern. Our main results are as following: (i) Part of the quenching takes place outside clusters through mass quenching and pre-processing. The e-folding timescale of this "ex-situ quenching phase" is roughly 3 Gyr with a strong inverse mass dependence. (ii) The pace of quenching is maintained roughly for 2 Gyr ("delay time") during the first crossing time into the cluster. During the delay time, quenching remains gentle probably because gas loss happens primarily on hot and neutral gases. (iii) Quenching becomes more dramatic (e-folding timescale of roughly 1 Gyr) after delay time, probably because ram pressure stripping is strongest near the cluster center. Counter-intuitively, more massive galaxies show shorter quenching timescales mainly because they enter their clusters with lower gas fractions due to ex-situ quenching.

The environmental effect is commonly used to explain the excess of gas-poor galaxies in galaxy clusters. Meanwhile, the presence of gas-poor galaxies at cluster outskirts, where galaxies have not spent enough time to feel the cluster environmental effect, hints at the presence of preprocessing. Using cosmological hydrodynamic simulations on 16 clusters, we investigate the mechanisms of gas depletion of galaxies found inside clusters. The gas-depletion mechanisms can be categorized into three channels based on where and when they took place. First, 34% of our galaxies are gas poor before entering clusters (“preprocessing”). They are mainly satellites that have undergone the environmental effect inside group halos. Second, 43% of the sample quickly became gas deficient in clusters before the first pericentric pass (“fast cluster processing”). Some of them were group satellites that are low in gas at the time of cluster entry compared to the galaxies directly coming from the field. Even the galaxies with large gas fractions take this channel if they fall into massive clusters (≳1014.5 M ⊙) or approach cluster centers through radial orbits. Third, 24% of our sample retain gas even after their first pericentric pass (“slow cluster processing”) as they fall into the less massive clusters or have circular orbits. The relative importance of each channel varies with a cluster’s mass, while the exact degree of significance is subject to large uncertainties. Group preprocessing accounts for one-third of the total gas depletion, but it also determines the gas fraction of galaxies at their cluster entry, which in turn determines whether a galaxy should take the fast or slow cluster processing.