Hydrodynamical Simulations Of Galaxy Formation Research Articles

The signature of galaxy formation models in the power spectrum of the hydrogen 21 cm line during reionization

ABSTRACT Observations of the 21 cm line of hydrogen are poised to revolutionize our knowledge of reionization and the first galaxies. However, harnessing such information requires robust and comprehensive theoretical modelling. We study the non-linear effects of hydrodynamics and astrophysical feedback processes, including stellar and AGN feedback, on the 21 cm signal by post-processing three existing cosmological hydrodynamical simulations of galaxy formation: Illustris, IllustrisTNG, and Eagle. Despite their different underlying galaxy-formation models, the simulations return similar predictions for the global 21 cm brightness temperature and its power spectrum. At fixed redshift, most differences are attributable to alternative reionization histories, in turn driven by differences in the buildup of stellar sources of radiation. However, several astrophysical processes imprint signatures in the 21 cm power spectrum at two key scales. First, we find significant small-scale ($k \ge 10\, \rm {Mpc}^{-1}$) differences between Illustris and IllustrisTNG, where higher velocity winds generated by supernova feedback soften density peaks, leading to lower 21 cm power in TNG. Thus, constraints at these scales could rule out extreme feedback models. Secondly, we find more 21 cm power at intermediate scales ($k \approx 0.8\, \rm {Mpc}^{-1}$) in Eagle due to ionization differences driven by highly effective stellar feedback, resulting in lower star formation, older and redder stellar populations, and lower ionizing luminosities for $M_h \gt 10^9 \, \rm M_\odot$. Different source models can manifest similarly in the 21 cm power spectrum, leading to often ignored degeneracies. These subtle features could allow future observations of the 21 cm signal, in conjunction with other observables, to constrain theoretical models for galactic feedback at high redshift.

Hydrodynamical simulations of galaxy formation with non-Gaussian initial conditions

Collisionless simulations of structure formation with significant local primordial non-Gaussianities at Mpc scales have shown that a non-Gaussian tail favouring underdensities, with a negative f NL parameter, can significantly change the merging history of galaxy-sized dark matter halos, which then typically assemble later than in vanilla ΛCDM. Moreover, such a small-scale negative f NL could have interesting consequences for the cosmological S 8 tension. Here, we complement our previous work on collisionless simulations with new hydrodynamical simulations of galaxy formation in boxes of 30 Mpc/h, using the RAMSES code. In particular, we show that all feedback prescriptions being otherwise identical, simulations with a negative f NL ∼ -1000 on small scales, hence forming galaxies a bit later than in vanilla ΛCDM, allow to form simulated galaxies with more disky kinematics than in the vanilla case. Therefore, such small-scale primordial non-Gaussianities could potentially help alleviate, simultaneously, tensions in cosmology and galaxy formation. These hydrodynamical simulations on small scales will need to be complemented with larger box simulations with scale-dependent non-Gaussianities, to statistically confirm these trends and explore their observational consequences in further detail.

Computationally efficient reionization in a large hydrodynamic galaxy formation simulation

ABSTRACT Accuracy in the topology and statistics of a simulated epoch of reionization (EoR) are vital to draw connections between observations and physical processes. While full radiative transfer models produce the most accurate reionization models, they are highly computationally expensive, and are infeasible for the largest cosmological simulations. Instead, large simulations often include EoR models that are pre-computed via the initial density field, or post-processed where feedback effects are ignored. We introduce astrid-es, a resimulation of the astrid epoch of reionization 20 > z > 5.5 which includes an on-the-fly excursion-set reionization algorithm. astrid-es produces more accurate reionization histories without significantly impacting the computational time. This model directly utilizes the star particles produced in the simulation to calculate the EoR history and includes an ultraviolet (UV) background which heats the gas particles after their reionization. We contrast the reionization topology and statistics in astrid-es with the previously employed parametric reionization model, finding that in astrid-es, ionized regions are more correlated with galaxies, and the 21cm power spectrum shows an increase in large-scale power. We calculate the relation between the size of H ii regions and the UV luminosity of the brightest galaxy within them. Prior to the overlap phase, we find a power-law fit of log(R) = −0.314MUV − 2.550log(1 + z) + 7.408 with a standard deviation σR < 0.15 dex across all mass bins. We also examine the properties of haloes throughout reionization, finding that while the properties of haloes in the simulation are correlated with the redshift of reionization, they are not greatly affected by reionization itself.

Probing z ≳ six massive black holes with gravitational waves

ABSTRACT We investigate the coalescence of massive black hole ($M_{\rm BH}\gtrsim 10^{6}~\rm {\rm M}_{\odot }$) binaries (MBHBs) at 6 < z < 10 by adopting a suite of cosmological hydrodynamical simulations of galaxy formation, zoomed-in on biased (>3σ) overdense regions (Mh ∼ 1012 M⊙ dark matter haloes at z = 6) of the Universe. We first analyse the impact of different resolutions and AGN feedback prescriptions on the merger rate, assuming instantaneous mergers. Then, we compute the halo bias correction factor due to the overdense simulated region. Our simulations predict merger rates that range between 3 and 15 $\rm yr^{-1}$ at z ∼6, depending on the run considered, and after correcting for a bias factor of ∼20−30. For our fiducial model, we further consider the effect of delay in the MBHB coalescence due to dynamical friction. We find that 83 per cent of MBHBs will merge within the Hubble time, and 21 per cent within 1 Gyr, namely the age of the Universe at z > 6. We finally compute the expected properties of the gravitational wave (GW) signals and find the fraction of LISA detectable events with high signal-to-noise ratio (SNR > 5) to range between 66 per cent and 69 per cent. However, identifying the electro-magnetic counterpart of these events remains challenging due to the poor LISA sky localization that, for the loudest signals ($\mathcal {M}_c\sim 10^6~{{\rm M}_{\odot }}$ at z = 6), is around 10 $\rm deg^2$.

The galaxy size to halo spin relation of disc galaxies in cosmological hydrodynamical simulations

ABSTRACT In the standard disc galaxy formation model, the sizes of galactic discs are tightly related to the spin parameters λ of their dark matter haloes. The model has been wildly adopted by various semi-analytical galaxy formation models which have been extremely successful to interpret a large body of observational data. However, the size–λ correlation was rarely seen in most modern hydrodynamical simulations of galaxy formation. In this short paper, we make use of 4 sets of large hydrodynamical simulations to explore the size–spin parameter relation with a large sample of simulated disc galaxies and compare it with a popular disc galaxy formation model of Mo et al. (1998). Intriguingly, galactic sizes correlate with spin parameters of their dark matter haloes in the simulations developed by the IllustrisTNG collaborations, albeit the relation does not always agree with prediction of MMW98 model overall stellar mass range we examined. There is also a size–spin correlation for the Milky Way analogies in the EAGLE simulations, while it is relatively weaker than that of the IllustrisTNG counterparts. For the dwarfs in the simulations from the EAGLE collaboration, there is NULL correlation. We conclude that either the detailed subgrid physics or hydrodynamics solvers account for the size-spin parameter relation, which will be explored in our future work.

Accurate predictions from small boxes: variance suppression via the Zel'dovich approximation

Simulations have become an indispensable tool for accurate modelling of observables measured in galaxy surveys, but can be expensive if very large dynamic range in scale is required. We describe how to combine Lagrangian perturbation theory models with N-body simulations to reduce the effects of finite computational volume in the prediction of ensemble average properties in the simulations within the context of control variates. In particular we use the fact that Zel'dovich displacements, computed during initial condition generation for any simulation, correlate strongly with the final density field. Since all the correlators of biased tracers can be computed with arbitrary precision for these displacements, pairing the Zel'dovich `simulation' with the N-body realization allows hundredfold reductions in sample variance for power spectrum or correlation function estimation. Zel'dovich control variates can accurately extend matter or tracer field emulators to larger scales than previously possible, as well as improving measurements of statistics in simulations which are inherently limited to small volumes, such as hydrodynamical simulations of galaxy formation and reionization.

Metal Lines Associated with the Lyα Forest from eBOSS Data

We investigate the metal species associated with the Lyα forest in eBOSS quasar spectra. Metal absorption lines are revealed in stacked spectra from cross-correlating the selected Lyα absorbers in the forest and the flux fluctuation field. Up to 13 metal species are identified as being associated with relatively strong Lyα absorbers (those with flux fluctuations − 1.0 < δ Lyα < − 0.6 and with a neutral hydrogen column density of ∼ 1015−16 cm−2) over the absorber redshift range of 2 < z abs < 4. The column densities of these species decrease toward higher redshift and for weaker Lyα absorbers. From modeling the column densities of various species, we find that the column density pattern suggests contributions from multiple gas components, both in the circumgalactic medium (CGM) and the intergalactic medium (IGM). While the low-ionization species (e.g., C ii, Si ii, and Mg ii) can be explained by high-density, cool gas (T ∼ 104 K) from the CGM, the high-ionization species may reside in low-density or high-temperature gas in the IGM. The measurements provide inputs for modeling the metal contamination in the Lyα forest baryon acoustic oscillation measurements. Comparisons with metal absorptions in high-resolution quasar spectra and hydrodynamic galaxy formation simulations can further elucidate the physical conditions of these Lyα absorbers.

Redshift and stellar mass dependence of intrinsic shapes of disc-dominated galaxies from COSMOS observations below z = 1.0

ABSTRACT The high abundance of disc galaxies without a large central bulge challenges predictions of current hydrodynamic simulations of galaxy formation. We aim to shed light on the formation of these objects by studying the redshift and mass dependence of their intrinsic 3D shape distributions in the COSMOS galaxy survey below redshift z = 1.0. This distribution is inferred from the observed distribution of 2D shapes, using a reconstruction method which we test using hydrodynamic simulations. Our tests reveal a moderate bias for the inferred average disc circularity and relative thickness, but a large bias on the dispersion of these quantities. Applying the reconstruction method on COSMOS data, we find variations of the average disc circularity and relative thickness with redshift of around ∼1 per cent and ∼10 per cent, respectively, which is comparable to the error estimates on these quantities. The average relative disc thickness shows a significant mass dependence which can be accounted for by the scaling of disc radius with galaxy mass. We conclude that our data provides no evidence for a strong dependence of the average circularity and absolute thickness of disc-dominated galaxies on redshift and mass that is significant with respect to the statistical uncertainties in our analysis. These findings are expected in the absence of disruptive merging or feedback events that would affect galaxy shapes. They hence support a scenario where present-day discs form early ( z &amp;gt; 1.0) and subsequently undergo a tranquil evolution in isolation. However, more data and a better understanding of systematics are needed to reaffirm our results.

The BPT Diagram in Cosmological Galaxy Formation Simulations: Understanding the Physics Driving Offsets at High Redshift

Abstract The Baldwin, Philips, &amp; Terlevich diagram of [O iii]/Hβ versus [N ii]/Hα (hereafter N2-BPT) has long been used as a tool for classifying galaxies based on the dominant source of ionizing radiation. Recent observations have demonstrated that galaxies at z ∼ 2 reside offset from local galaxies in the N2-BPT space. In this paper, we conduct a series of controlled numerical experiments to understand the potential physical processes driving this offset. We model nebular line emission in a large sample of galaxies, taken from the simba cosmological hydrodynamic galaxy formation simulation, using the cloudy photoionization code to compute the nebular line luminosities from H ii regions. We find that the observed shift toward higher [O iii]/Hβ and [N ii]/Hα values at high redshift arises from sample selection: when we consider only the most massive galaxies M * ∼ 1010–11 M ⊙, the offset naturally appears, due to their high metallicities. We predict that deeper observations that probe lower-mass galaxies will reveal galaxies that lie on a locus comparable to z ∼ 0 observations. Even when accounting for samples-selection effects, we find that there is a subtle mismatch between simulations and observations. To resolve this discrepancy, we investigate the impact of varying ionization parameters, H ii region densities, gas-phase abundance patterns, and increasing radiation field hardness on N2-BPT diagrams. We find that either decreasing the ionization parameter or increasing the N/O ratio of galaxies at fixed O/H can move galaxies along a self-similar arc in N2-BPT space that is occupied by high-redshift galaxies.

How gas flows shape the stellar–halo mass relation in the eagle simulation

ABSTRACT The difference in shape between the observed galaxy stellar mass function and the predicted dark matter halo mass function is generally explained primarily by feedback processes. Feedback can shape the stellar–halo mass (SHM) relation by driving gas out of galaxies, by modulating the first-time infall of gas on to galaxies (i.e. preventative feedback), and by instigating fountain flows of recycled wind material. We present and apply a method to disentangle these effects for hydrodynamical simulations of galaxy formation. We build a model of linear coupled differential equations that by construction reproduces the flows of gas on to and out of galaxies and haloes in the eaglecosmological simulation. By varying individual terms in this model, we isolate the relative effects of star formation, ejection via outflow, first-time inflow, and wind recycling on the SHM relation. We find that for halo masses $M_{200} \lt 10^{12} \, \mathrm{M_\odot }$ the SHM relation is shaped primarily by a combination of ejection from galaxies and haloes, while for larger M200 preventative feedback is also important. The effects of recycling and the efficiency of star formation are small. We show that if, instead of M200, we use the cumulative mass of dark matter that fell in for the first time, the evolution of the SHM relation nearly vanishes. This suggests that the evolution is due to the definition of halo mass rather than to an evolving physical efficiency of galaxy formation. Finally, we demonstrate that the mass in the circumgalactic medium is much more sensitive to gas flows, especially recycling, than is the case for stars and the interstellar medium.

The Diversity of Environments around Luminous Quasars at Redshift z ∼ 6

Abstract Significant clustering around the rarest luminous quasars is a feature predicted by dark matter theory combined with number density matching arguments. However, this expectation is not reflected by observations of quasars residing in a diverse range of environments. Here we assess the tension in the diverse clustering of visible i-band dropout galaxies around luminous z ∼ 6 quasars. Our approach uses a simple empirical method to derive the median luminosity–to–halo mass relation, L c (M h ), for both quasars and galaxies under the assumption of lognormal luminosity scatter, Σ Q and Σ G . We show that higher Σ Q reduces the average halo mass hosting a quasar of a given luminosity, thus introducing at least a partial reversion to the mean in the number count distribution of nearby Lyman-break galaxies. We generate a large sample of mock Hubble Space Telescope fields of view centered across rare z ∼ 6 quasars by resampling pencil beams traced through the dark matter component of the BlueTides cosmological simulation. We find that diverse quasar environments are expected for Σ Q &gt; 0.4, consistent with numerous observations and theoretical studies. However, we note that the average number of galaxies around the central quasar is primarily driven by galaxy evolutionary processes in neighboring halos, as embodied by our parameter Σ G , instead of a difference in the large-scale structure around the central quasar host, embodied by Σ Q . We conclude that models with Σ G &gt; 0.3 are consistent with current observational constraints on high-z quasars, and that such a value is comparable to the scatter estimated from hydrodynamical simulations of galaxy formation.

Galaxy formation in the brane world I: overview and first results

ABSTRACT We carry out ‘full-physics’ hydrodynamical simulations of galaxy formation in the normal-branch Dvali–Gabadadze–Porrati (nDGP) braneworld model using a new modified version of the arepo code and the IllustrisTNG galaxy formation model. We simulate two nDGP models (N5 and N1) that represent, respectively, weak and moderate departures from general relativity (GR), in boxes of sizes $62$ and $25\, h^{-1}{\rm Mpc}$ using 2 × 5123 dark matter particles and initial gas cells. This allows us to explore, for the first time, the impact of baryonic physics on galactic scales in braneworld models of modified gravity and to make predictions on the stellar content of dark matter haloes and galaxy evolution through cosmic time in these models. We find significant differences between the GR and nDGP models in the power spectra and correlation functions of gas, stars and dark matter of up to ∼25 per cent on large scales. Similar to their impact in the standard cosmological model (Λ cold dark matter), baryonic effects can have a significant influence over the clustering of the overall matter distribution, with a sign that depends on scale. Studying the degeneracy between modified gravity and galactic feedback in these models, we find that these two physical effects on matter clustering can be cleanly disentangled, allowing for a method to accurately predict the matter power spectrum with baryonic effects included, without having to run hydrodynamical simulations. Depending on the braneworld model, we find differences compared with GR of up to ∼15 per cent in galaxy properties such as the stellar-to-halo-mass ratio, galaxy stellar mass function, gas fraction, and star formation rate density. The amplitude of the fifth force is reduced by the presence of baryons in the very inner part of haloes, but this reduction quickly becomes negligible above ∼0.1 times the halo radius.

A giant central red disk galaxy at redshift z = 0.76: Challenge to theories of galaxy formation

We report a giant red central disk galaxy in the XMM-LSS north region. The region is covered with a rich variety of multiband photometric and spectroscopic observations. Using the photometric data of the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) and spectroscopic observation of the Baryon Oscillation Spectroscopic Survey (BOSS), we find that the galaxy has a stellar mass of $\sim 10^{11.6}$ times of the solar mass $M_\odot$. The galaxy has a red color and has an old stellar population, and thus its star formation has stopped. With the photometric image data of Hyper Suprime-Cam (HSC) Subaru Strategic Program, we demonstrate that its luminosity profile is perfectly described by a S\'ersic form with $n=1.22$ indicating disk morphology. We also analyze its environment based on the VIMOS Public Extragalactic Redshift Survey (VIPERS) photometric catalog, and find that its close neighbors are all less massive, indicating that our observed galaxy is sitting at the center of its host halo. Existence of the giant red central disk galaxy seriously challenges the current standard paradigm of galaxy formation, as there is no known physical mechanism to explain the quenching of its star formation. This conclusion is supported by state-of-the-art hydrodynamical simulations of galaxy formation.

Simulating dust grain-radiation coupling on a moving mesh

ABSTRACT We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code arepo-rt, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in arepo-rt, which treats UV, optical photons as single scattered and IR photons as multi scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.

How Well Can We Measure the Stellar Mass of a Galaxy: The Impact of the Assumed Star Formation History Model in SED Fitting

Abstract The primary method for inferring the stellar mass (M *) of a galaxy is through spectral energy distribution (SED) modeling. However, the technique rests on assumptions such as the galaxy star formation history (SFH) and dust attenuation law that can severely impact the accuracy of derived physical properties from SED modeling. Here we examine the effect that the assumed SFH has on the stellar properties inferred from SED fitting by ground-truthing them against mock observations of high-resolution cosmological hydrodynamic galaxy formation simulations. Classically, SFHs are modeled with simplified parameterized functional forms, but these forms are unlikely to capture the true diversity of galaxy SFHs and may impose systematic biases with underreported uncertainties on results. We demonstrate that flexible nonparametric SFHs outperform traditional parametric forms in capturing variations in galaxy SFHs and, as a result, lead to significantly improved stellar masses in SED fitting. We find a decrease in the average bias of 0.4 dex with a delayed-τ model to a bias under 0.1 dex for the nonparametric model, though this is heavily dependent on the choice of prior for the nonparametric model. Similarly, using nonparametric SFHs in SED fitting results in increased accuracy in recovered galaxy star formation rates and stellar ages.

Indirectly Measuring Stellar Velocity Dispersions in High-redshift Disk Galaxies

Abstract In the Milky Way, older stars tend to have higher vertical velocity dispersions and comprise a thicker component of the stellar disk. Historically, this has been attributed to the gradual vertical heating of stellar orbits. A competing explanation suggests that stars may have formed in thicker disks at earlier times, inheriting the turbulent velocities of the gas in which they formed. As the gas velocity dispersion in disks declines with time, stars should form in thinner disks. Both explanations are consistent with today’s Milky Way. We present a novel technique to break this degeneracy by indirectly measuring the vertical stellar velocity dispersion of high redshift galaxies using Hubble Space Telescope imaging. We assess this method using a high-resolution hydrodynamic galaxy formation simulation. By comparing the true vertical velocity dispersion of the simulated disk against what is inferred through its synthetic images, we demonstrate the potential of this approach.

Evaluating hydrodynamical simulations with green valley galaxies

ABSTRACT We test cosmological hydrodynamical simulations of galaxy formation regarding the properties of the blue cloud (BC), green valley (GV), and red sequence (RS), as measured on the 4000Å break strength versus stellar mass plane at z = 0.1. We analyse the RefL0100N1504 run of EAGLE and the TNG100 run of IllustrisTNG project, by comparing them with the Sloan Digital Sky Survey (SDSS), while taking into account selection bias. Our analysis focuses on the GV, within stellar mass $\log \, \mathrm{{\it M}_\star /M_{\odot }} \simeq 10\!-\!11$, selected from the bimodal distribution of galaxies on the Dn(4000) versus stellar mass plane, following Angthopo et al. methodology. Both simulations match the fraction of AGN in the GV. However, they overproduce quiescent GV galaxies with respect to observations, with IllustrisTNG yielding a higher fraction of quiescent GV galaxies than EAGLE. In both, GV galaxies have older luminosity-weighted ages with respect to the SDSS, while a better match is found for mass-weighted ages. We find EAGLE GV galaxies quench their star formation early, but undergo later episodes of star formation, matching observations. In contrast, IllustrisTNG GV galaxies have a more extended star formation history, and quench more effectively at later cosmic times, producing the excess of quenched galaxies in GV compared with SDSS, based on the 4000Å break strength. These results suggest the AGN feedback subgrid physics, more specifically, the threshold halo mass for black hole input and the black hole seed mass, could be the primary cause of the overproduction of quiescent galaxies found with respect to the observational constraints.

The detailed structure and the onset of galaxy formation in low-mass gaseous dark matter haloes

ABSTRACT We present a model for the formation of the first galaxies before and after the reionization of hydrogen in the early universe. In this model, galaxy formation can only take place in dark matter haloes whose mass exceeds a redshift-dependent critical value, which, before reionization, is equal (in the simplest case) to the mass at which atomic hydrogen cooling becomes effective and, after reionization, is equal to the mass above which gas cannot remain in hydrostatic equilibrium. We define the Halo Occupation Fraction (HOF) as the fraction of haloes that host a luminous galaxy as a function of halo mass. The HOF is established by the interplay between the evolution of the critical mass and the assembly history of haloes and depends on three factors: the minimum halo mass for galaxy formation before reionization, the redshift of reionization, and the intensity of the (evolving) external photoheating rate. Our fiducial model predicts a cutoff in the galaxy mass function at a present-day halo mass, $M_{200} \sim 3\times 10^{8} \, \mathrm{M}_{\odot }$; 100 per cent occupation at $M_{200} \gt 5\times 10^9 \, \mathrm{M}_{\odot }$; and a population of starless gaseous haloes of present-day mass in the range 106 ≲ M200/M⊙ ≲ 5 × 109, in which the gas is in thermal equilibrium with the ultraviolet background radiation and in hydrostatic equilibrium in the gravitational potential of the halo. The transition between HOF = 0 and HOF = 1 reflects the stochastic nature of halo mass growth. We explore how these characteristic masses vary with model assumptions and parameter values. The results of our model are in excellent agreement with cosmological hydrodynamic simulations of galaxy formation.

A new model for including galactic winds in simulations of galaxy formation – I. Introducing the Physically Evolved Winds (PhEW) model

ABSTRACT The propagation and evolution of cold galactic winds in galactic haloes is crucial to galaxy formation models. However, modelling of this process in hydrodynamic simulations of galaxy formation is oversimplified owing to a lack of numerical resolution and often neglects critical physical processes such as hydrodynamic instabilities and thermal conduction. We propose an analytic model, Physically Evolved Winds, that calculates the evolution of individual clouds moving supersonically through a uniform ambient medium. Our model reproduces predictions from very high resolution cloud-crushing simulations that include isotropic thermal conduction over a wide range of physical conditions. We discuss the implementation of this model into cosmological hydrodynamic simulations of galaxy formation as a subgrid prescription to model galactic winds more robustly both physically and numerically.

Modelling the tightest relation between galaxy properties and dark matter halo properties from hydrodynamical simulations of galaxy formation

ABSTRACT We investigate how a property of a galaxy correlates most tightly with a property of its host dark matter halo, using state-of-the-art hydrodynamical simulations of galaxy formation: EAGLE, Illustris, and IllustrisTNG. Unlike most of the previous work, our analyses focus on all types of galaxies, including both central and satellite galaxies. We find that the stellar mass of a galaxy at the epoch of the peak circular velocity with an evolution correction gives the tightest such correlation to the peak circular velocity Vpeak of the galaxy’s underling dark matter halo. The evolution of galaxy stellar mass reduces rather than increases scatter in such a relation. We also find that one major source of scatter comes from star stripping due to the strong interactions between galaxies. Even though, we show that the size of scatter predicted by hydrodynamical simulations has a negligible impact on the clustering of dense Vpeak-selected subhalo from simulations, which suggests that even the simplest subhalo abundance matching (SHAM), without scatter and any additional free parameter, can provide a robust prediction of galaxy clustering that can agree impressively well with the observations from the Sloan Digital Sky Survey (SDSS) main galaxy survey.