Cosmological simulations are reaching the resolution necessary to study ultra-faint dwarf galaxies. Observations indicate that in small populations, the stellar initial mass function (IMF) is not fully populated; rather, stars are sampled in a way that can be approximated as coming from an underlying probability density function. To ensure the accuracy of cosmological simulations in the ultra-faint regime, we present an improved treatment of the IMF. We implement a self-consistent, stochastically populated IMF in cosmological hydrodynamic simulations. We test our method using high-resolution simulations of a Milky Way halo, run to z = 6, yielding a sample of nearly 100 galaxies. We also use an isolated dwarf galaxy to investigate the resulting systematic differences in galaxy properties. We find that a stochastic IMF in simulations makes feedback burstier, strengthening feedback, and quenching star formation earlier in small dwarf galaxies. For galaxies in haloes with mass ≲ 108.5 M⊙, a stochastic IMF typically leads to lower stellar mass compared to a continuous IMF, sometimes by more than an order of magnitude. We show that existing methods of ensuring discrete supernovae incorrectly determine the mass of the star particle and its associated feedback. This leads to overcooling of surrounding gas, with at least ∼10 per cent higher star formation and ∼30 per cent higher cold gas content. Going forwards, to accurately model dwarf galaxies and compare to observations, it will be necessary to incorporate a stochastically populated IMF that samples the full spectrum of stellar masses.
- Award ID(s):
- 1906976
- NSF-PAR ID:
- 10287438
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 503
- Issue:
- 2
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 2380 to 2405
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
The condensation of baryons within a dark matter (DM) halo during galaxy formation should result in some contraction of the halo as the combined system settles into equilibrium. We quantify this effect on the cuspy primordial halos predicted by DM-only simulations for the baryon distributions observed in the galaxies of the SPARC database. We find that the DM halos of high surface brightness galaxies (with Σ eff ≳ 100 L ⊙ pc −2 at 3.6 μm) experience strong contraction. Halos become more cuspy as a result of compression: the inner DM density slope increases with the baryonic surface mass density. We iteratively fit rotation curves to find the balance between initial halo parameters (constrained by abundance matching), compression, and stellar mass-to-light ratio. The resulting fits often require lower stellar masses than expected for stellar populations, particularly in galaxies with bulges: stellar mass must be reduced to make room for the DM it compresses. This trade off between dark and luminous mass is reminiscent of the cusp-core problem in dwarf galaxies, but occurs in more massive systems: the present-epoch DM halos cannot follow from cuspy primordial halos unless (1) the stellar mass-to-light ratios are systematically smaller than expected from standard stellar population synthesis models, and/or (2) there is a net outward mass redistribution from the initial cusp, even in massive galaxies widely considered to be immune from such effects.more » « less
-
more » « less
ABSTRACT The mass distribution in massive elliptical galaxies encodes their evolutionary history, thus providing an avenue to constrain the baryonic astrophysics in their evolution. The power-law assumption for the radial mass profile in ellipticals has been sufficient to describe several observables to the noise level, including strong lensing and stellar dynamics. In this paper, we quantitatively constrained any deviation, or the lack thereof, from the power-law mass profile in massive ellipticals through joint lensing–dynamics analysis of a large statistical sample with 77 galaxy–galaxy lens systems. We performed an improved and uniform lens modelling of these systems from archival Hubble Space Telescope imaging using the automated lens modelling pipeline dolphin. We combined the lens model posteriors with the stellar dynamics to constrain the deviation from the power law after accounting for the line-of-sight lensing effects, a first for analyses on galaxy–galaxy lenses. We find that the Sloan Lens ACS Survey lens galaxies with a mean redshift of 0.2 are consistent with the power-law profile within 1.1σ (2.8σ) and the Strong Lensing Legacy Survey lens galaxies with a mean redshift of 0.6 are consistent within 0.8σ (2.1σ), for a spatially constant (Osipkov–Merritt) stellar anisotropy profile. We adopted the spatially constant anisotropy profile as our baseline choice based on previous dynamical observables of local ellipticals. However, spatially resolved stellar kinematics of lens galaxies are necessary to differentiate between the two anisotropy models. Future studies will use our lens models to constrain the mass distribution individually in the dark matter and baryonic components.
-
more » « less
ABSTRACT Velocity dispersion (σ) is a key driver for galaxy structure and evolution. We here present a comprehensive semi-empirical approach to compute σ via detailed Jeans modelling assuming both a constant and scale-dependent mass-to-light ratio M*/L. We compare with a large sample of local galaxies from MaNGA and find that both models can reproduce the Faber–Jackson (FJ) relation and the weak dependence of σ on bulge-to-total (B/T) ratio (for B/T ≳ 0.25). The dynamical-to-stellar mass ratio within R ≲ Re can be fully accounted for by a gradient in M*/L. We then build velocity dispersion evolutionary tracks σap[M*, z] (within an aperture) along the main progenitor dark matter haloes assigning stellar masses, effective radii, and Sérsic indices via a variety of abundance matching and empirically motivated relations. We find: (1) clear evidence for downsizing in σap[M*, z] along the progenitor tracks; (2) at fixed stellar mass σ ∝ (1 + z)0.2−0.3 depending on the presence or not of a gradient in M*/L. We extract σap[M*, z] from the TNG50 hydrodynamic simulation and find very similar results to our models with constant M*/L. The increasing dark matter fraction within Re tends to flatten the σap[M*, z] along the progenitors at z ≳ 1 in constant M*/L models, while σap[M*, z] have a steeper evolution in the presence of a stellar gradient. We then show that a combination of mergers and gas accretion is likely responsible for the constant or increasing σap[M*, z] with time. Finally, our σap[M*, z] are consistent with a nearly constant and steep Mbh − σ relation at z ≲ 2, with black hole masses derived from the LX − M* relation.
-
more » « less
ABSTRACT Gradients in the stellar populations (SP) of galaxies – e.g. in age, metallicity, stellar initial mass function (IMF) – can result in gradients in the stellar-mass-to-light ratio, M*/L. Such gradients imply that the distribution of the stellar mass and light is different. For old SPs, e.g. in early-type galaxies at z ∼ 0, the M*/L gradients are weak if driven by variations in age and metallicity, but significantly larger if driven by the IMF. A gradient which has larger M*/L in the centre increases the estimated total stellar mass (M*) and reduces the scale which contains half this mass (Re,*), compared to when the gradient is ignored. For the IMF gradients inferred from fitting MILES simple SP models to the H β, 〈Fe〉, [MgFe], and TiO2SDSS absorption lines measured in spatially resolved spectra of early-type galaxies in the MaNGA survey, the fractional change in Re,* can be significantly larger than that in M*, especially when the light is more centrally concentrated. The Re,*–M* correlation which results from accounting for IMF gradients is offset to smaller sizes by 0.3 dex compared to when these gradients are ignored. Comparisons with ‘quiescent’ galaxies at higher z must account for evolution in SP gradients (especially age and IMF) and in the light profile before drawing conclusions about how Re,* and M* evolve. The implied merging between higher z and the present is less contrived if Re,*/Re at z ∼ 0 is closer to our IMF-driven gradient calibration than to unity.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/mnras/stab536
