Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userHow to empirically model star formation in dark matter halos: I. Inferences about central galaxies from numerical simulations
145
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We use TNG and EAGLE hydrodynamic simulations to investigate the central galaxy - dark matter halo relations that are needed for a halo-based empirical model of star formation in galaxies. Using a linear dimension reduction algorithm and a model ensemble method, we find that for both star-forming and quenched galaxies, the star formation history (SFH) is tightly related to the halo mass assembly history (MAH). The quenching of a low-mass galaxy is mainly due to the infall-ejection process related to a nearby massive halo, while the quenching of a high-mass galaxy is closely related to the formation of a massive progenitor in its host halo. The classification of star-forming and quenched populations based solely on halo properties contains contamination produced by sample imbalance and overlapping distributions of the two populations. Guided by the results from hydrodynamic simulations, we build an empirical model to predict the SFH of central galaxies based on the MAH of their host halos, and we model the star-forming and quenched populations separately. Our model is based on the idea of adopting star formation templates from hydrodynamic simulations to reduce model complexity. We use various tests to demonstrate that the model can recover star formation histories of individual galaxies, and can statistically reproduce the galaxy bimodal distribution, stellar mass - halo mass and star formation rate - halo mass relations from low to high redshift, and assembly bias. Our study provides a framework of using hydrodynamic simulations to discover, and to motivate the use of, key ingredients to model galaxy formation using halo properties.
rate research
Read More
Using the self-consistent modeling of the conditional stellar mass functions across cosmic time by Yang et al. (2012), we make model predictions for the star formation histories (SFHs) of {it central} galaxies in halos of different masses. The model requires the following two key ingredients: (i) mass assembly histories of central and satellite galaxies, and (ii) local observational constraints of the star formation rates of central galaxies as function of halo mass. We obtain a universal fitting formula that describes the (median) SFH of central galaxies as function of halo mass, galaxy stellar mass and redshift. We use this model to make predictions for various aspects of the star formation rates of central galaxies across cosmic time. Our main findings are the following. (1) The specific star formation rate (SSFR) at high $z$ increases rapidly with increasing redshift [$propto (1+z)^{2.5}$] for halos of a given mass and only slowly with halo mass ($propto M_h^{0.12}$) at a given $z$, in almost perfect agreement with the specific mass accretion rate of dark matter halos. (2) The ratio between the star formation rate (SFR) in the main-branch progenitor and the final stellar mass of a galaxy peaks roughly at a constant value, $sim 10^{-9.3} h^2 {rm yr}^{-1}$, independent of halo mass or the final stellar mass of the galaxy. However, the redshift at which the SFR peaks increases rapidly with halo mass. (3) More than half of the stars in the present-day Universe were formed in halos with $10^{11.1}msunh < M_h < 10^{12.3}msunh$ in the redshift range $0.4 < z < 1.9$. (4) ... [abridged]
We study the implementation of mechanical feedback from supernovae (SNe) and stellar mass loss in galaxy simulations, within the Feedback In Realistic Environments (FIRE) project. We present the FIRE-2 algorithm for coupling mechanical feedback, which can be applied to any hydrodynamics method (e.g. fixed-grid, moving-mesh, and mesh-less methods), and black hole as well as stellar feedback. This algorithm ensures manifest conservation of mass, energy, and momentum, and avoids imprinting preferred directions on the ejecta. We show that it is critical to incorporate both momentum and thermal energy of mechanical ejecta in a self-consistent manner, accounting for SNe cooling radii when they are not resolved. Using idealized simulations of single SN explosions, we show that the FIRE-2 algorithm, independent of resolution, reproduces converged solutions in both energy and momentum. In contrast, common fully-thermal (energy-dump) or fully-kinetic (particle-kicking) schemes in the literature depend strongly on resolution: when applied at mass resolution >100 solar masses, they diverge by orders-of-magnitude from the converged solution. In galaxy-formation simulations, this divergence leads to orders-of-magnitude differences in galaxy properties, unless those models are adjusted in a resolution-dependent way. We show that all models that individually time-resolve SNe converge to the FIRE-2 solution at sufficiently high resolution. However, in both idealized single-SN simulations and cosmological galaxy-formation simulations, the FIRE-2 algorithm converges much faster than other sub-grid models without re-tuning parameters.
We present a robust method to constrain average galaxy star formation rates, star formation histories, and the intracluster light as a function of halo mass. Our results are consistent with observed galaxy stellar mass functions, specific star formation rates, and cosmic star formation rates from z=0 to z=8. We consider the effects of a wide range of uncertainties on our results, including those affecting stellar masses, star formation rates, and the halo mass function at the heart of our analysis. As they are relevant to our method, we also present new calibrations of the dark matter halo mass function, halo mass accretion histories, and halo-subhalo merger rates out to z=8. We also provide new compilations of cosmic and specific star formation rates; more recent measurements are now consistent with the buildup of the cosmic stellar mass density at all redshifts. Implications of our work include: halos near 10^12 Msun are the most efficient at forming stars at all redshifts, the baryon conversion efficiency of massive halos drops markedly after z ~ 2.5 (consistent with theories of cold-mode accretion), the ICL for massive galaxies is expected to be significant out to at least z ~ 1-1.5, and dwarf galaxies at low redshifts have higher stellar mass to halo mass ratios than previous expectations and form later than in most theoretical models. Finally, we provide new fitting formulae for star formation histories that are more accurate than the standard declining tau model. Our approach places a wide variety of observations relating to the star formation history of galaxies into a self-consistent framework based on the modern understanding of structure formation in LCDM. Constraints on the stellar mass-halo mass relationship and star formation rates are available for download at http://www.peterbehroozi.com/data.html .
Multicomponent dark matter with self-interactions, which allows for inter-
Self-gravitating astronomical objects often show a central plateau in the density profile (core) whose physical origin is hotly debated. Cores are theoretically expected in N-body systems of maximum entropy, however, they are not present in the canonical N-body numerical simulations of cold dark matter (CDM). Our work shows that despite this apparent contradiction between theory and numerical simulations, they are fully consistent. Simply put, cores are characteristic of systems in thermodynamic equilibrium, but thermalizing collisions are purposely suppressed in CDM simulations. When collisions are allowed, N-body numerical simulations develop cored density profiles, in perfect agreement with the theoretical expectation. We compare theory and two types of numerical simulations: (1) when DM particles are self-interacting (SIDM) with enough cross-section, then the effective two-body relaxation timescale becomes shorter than the Hubble time resulting in cored DM haloes. The haloes thus obtained, with masses from dwarf galaxies to galaxy clusters, collapse to a single shape after normalization, and this shape agrees with the polytropic density profile theoretically expected. (2) The inner radii in canonical N-body numerical simulations are always discarded because the use of finite-mass DM particles artificially increases the two-body collision rate. We show that the discarded radii develop cores that are larger than the employed numerical softening and have polytropic shapes independently of halo mass. Our work suggests that the presence of cores in simulated (or observed) density profiles can used as evidence for systems in thermodynamic equilibrium.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
