Low-Surface-Brightness Galaxies are missing in the observed Stellar Mass Function

We investigate the impact of the surface brightness (SB) limit on the galaxy stellar mass functions (GSMFs) using mock surveys generated from the Horizon Run 5 (HR5) simulation. We compare the stellar-to-halo-mass relation, GSMF, and size-stellar mass relation of the HR5 galaxies with empirical data and other cosmological simulations. The mean SB of simulated galaxies are computed using their effective radii, luminosities, and colors. To examine the cosmic SB dimming effect, we compute $k$-corrections from the spectral energy distributions of individual simulated galaxy at each redshift, apply the $k$-corrections to the galaxies, and conduct mock surveys based on the various SB limits. We find that the GSMFs are significantly affected by the SB limits at a low-mass end. This approach can ease the discrepancy between the GSMFs obtained from simulations and observations at $0.625\le z\le 2$. We also find that a redshift survey with a SB selection limit of \left^e = 28 mag arcsec${}^{-2}$ will miss 20% of galaxies with $M_\star^g=10^{9}~{\rm M_\odot}$ at $z=0.625$. The missing fraction of low-surface-brightness galaxies increases to 50%, 70%, and 98% at $z=0.9$, 1.1, and 1.9, respectively, at the SB limit.Comment: 27 pages, 30 figures, accepted for publication in Ap

Identification of Galaxy Protoclusters Based on the Spherical Top-hat Collapse Theory

We propose a new method for finding galaxy protoclusters that is motivated by structure formation theory and also directly applicable to observations. We adopt the conventional definition that a protocluster is a galaxy group whose virial mass $M_{\rm vir} < M_{\rm cl}$ at its epoch, where $M_{\rm cl}=10^{14}\,M_{\odot}$, but would exceed that limit when it evolves to $z=0$. We use the critical overdensity for complete collapse at $z = 0$ predicted by the spherical top-hat collapse model to find the radius and total mass of the regions that would collapse at $z=0$. If the mass of a region centered at a massive galaxy exceeds $M_{\rm cl}$, the galaxy is at the center of a protocluster. We define the outer boundary of protocluster as the zero-velocity surface at the turnaround radius so that the member galaxies are those sharing the same protocluster environment and showing some conformity in physical properties. We use the cosmological hydrodynamical simulation Horizon Run 5 (HR5) to calibrate this prescription and demonstrate its performance. We find that the protocluster identification method suggested in this study is quite successful. Its application to the high-redshift HR5 galaxies shows a tight correlation between the mass within the protocluster regions identified according to the spherical collapse model and the final mass to be found within the clusters at $z=0$, meaning that the regions can be regarded as the bona fide protoclusters with high reliability. We also confirm that the redshift-space distortion does not significantly affect the performance of the protocluster identification scheme.Comment: 23 pages, 21 figures. Accepted for publication in Ap

The Horizon Run 5 Cosmological Hydrodynamic Simulation: Probing Galaxy Formation from Kilo- to Giga-parsec Scales

Horizon Run 5 (HR5) is a cosmological hydrodynamical simulation which captures the properties of the Universe on a Gpc scale while achieving a resolution of 1kpc. Inside the simulation box we zoom-in on a high-resolution cuboid region with a volume of 1049×114×114cMpc3.The sub-grid physics chosen to model galaxy formation includes radiative heating/cooling, UV background, star formation, supernova feedback, chemical evolution tracking the enrichment of oxygen and iron, the growth of supermassive black holes and feedback from active galactic nuclei (AGN) in the form of a dual jet-heating mode. For this simulation we implemented a hybrid MPI-OMP version of RAMSES, specifically targeted for modern many-core many thread parallel architectures. In addition to the traditional simulation snapshots, light-cone data was generated on the fly. For the post-processing, we extended the Friends-of-Friend (FoF) algorithm and developed a new galaxy finder PGalF to analyse the outputs of HR5. The simulation successfully reproduces observations, such as the cosmic star formation history and connectivity of galaxy distribution, We identify cosmological structures at a wide range of scales, from filaments with a length of several cMpc, to voids with a radius of ~100 cMpc. The simulation also indicates that hydrodynamical effects on small scales impact galaxy clustering up to very large scales near and beyond the baryonic acoustic oscillation (BAO) scale. Hence, caution should be taken when using that scale as a cosmic standard ruler: one needs to carefully understand the corresponding biases. The simulation is expected to be an invaluable asset for the interpretation of upcoming deep surveys of the Universe

High-resolution three-dimensional simulations of gas removal from ultrafaint dwarf galaxies

Context. The faintest Local Group galaxies found lurking in and around the Milky Way halo provide a unique test bed for theories of structure formation and evolution on small scales. Deep Subaru and Hubble Space Telescope photometry demonstrates that the stellar populations of these galaxies are old and that the star formation activity did not last longer than 2 Gyr in these systems. A few mechanisms that may lead to such a rapid quenching have been investigated by means of hydrodynamic simulations, but these have not provided any final assessment so far. Aims. This is the first in a series of papers aimed at analyzing the roles of stellar feedback, ram pressure stripping, host-satellite tidal interactions, and reionization in cleaning the lowest mass Milky Way companions of their cold gas using high-resolution, three-dimensional hydrodynamic simulations. Methods. We simulated an isolated ultrafaint dwarf galaxy loosely modeled after Boötes I, and examined whether or not stellar feedback alone could drive a substantial fraction of the ambient gas out from the shallow potential well. Results. In contrast to simple analytical estimates, but in agreement with previous hydrodynamical studies, we find that most of the cold gas reservoir is retained. Conversely, a significant amount of the metal-enriched stellar ejecta crosses the boundaries of the computational box with velocities exceeding the local escape velocity and is, thus, likely lost from the system. Conclusions. Although the total energy output from multiple supernova explosions exceeds the binding energy of the gas, no galactic-scale outflow develops in our simulations and as such, most of the ambient medium remains trapped within the weak potential well of the model galaxy. It seems thus unavoidable that to explain the dearth of gas in ultrafaint dwarf galaxies, we will have to resort to environmental effects. This will be the subject of a forthcoming paper

Low-surface-brightness Galaxies are Missing in the Observed Stellar Mass Function

We investigate the impact of the surface-brightness (SB) limit on the galaxy stellar mass functions (GSMFs) using galaxy catalogs generated from the Horizon Run 5 ( HR5 ) simulation. We compare the stellar-to-halo-mass relation, GSMF, and size–stellar mass relation of the HR5 galaxies with observational data and other cosmological simulations. The mean SB of simulated galaxies are computed using their effective radii, luminosities, and colors. To examine the cosmic SB dimming effect, we compute k -corrections from the spectral energy distributions of individual simulated galaxy at each redshift, apply the k -corrections to the galaxies, and conduct mock surveys based on the various SB limits. We find that the GSMFs are significantly affected by the SB limits at the low-mass end. This approach can ease the discrepancy between the GSMFs obtained from simulations and observations at 0.6 ≲ z ≤ 2. We also find that a redshift survey with an SB selection limit of ${\left\langle {\mu }_{r}\right\rangle }^{e}=$ 25 mag arcsec ^−2 will miss 20% of galaxies with ${M}_{\star }^{g}={10}^{9}\,{M}_{\odot }$ at z = 0.625. The missing fraction of low-surface-brightness galaxies increases to 35%, 55%, and 80% at z = 0.9, 1.1, and 1.9, respectively, at the same SB limit

Formation and Morphology of the First Galaxies in the Cosmic Morning

International audienceWe investigate the formation and morphological evolution of the first galaxies in the cosmic morning (10 ≳ z ≳ 4) using the Horizon Run 5 (HR5) simulation. For galaxies above the stellar mass ${M}_{\star ,\min }=2\times {10}^{9}\,{M}_{\odot }$ , we classify them into disk, spheroid, and irregular types according to their asymmetry and stellar-mass morphology. We find that about two-thirds of the galaxies have a Sérsic index 10 M ⊙. Almost all the first galaxies with ${M}_{\star } {M}_{\star ,\min }$ at z > 4 form at initial peaks of the matter-density field. Large-scale structures in the universe emerge and grow like cosmic rhizomes as the underlying matter-density fluctuations grow and form associations of galaxies in rare overdense regions and the realm of the galactic world is stretched into relatively lower-density regions along evolving filaments. The cosmic web of galaxies forms at lower redshifts when most rhizomes globally percolate. The primordial angular momentum produced by the induced tidal torques on protogalactic regions is correlated with the internal kinematics of galaxies and tightly aligned with the angular momentum of the total galaxy mass. The large-scale tidal field imprinted in the initial conditions seems responsible for the dominance of disk morphology and for the tendency of galaxies to reacquire a disk postdistortion