scholarly journals The impact of stellar feedback on high-z galaxy populations

2015 ◽  
Vol 11 (S319) ◽  
pp. 26-26
Author(s):  
Michaela Hirschmann ◽  
Gabriella De Lucia

AbstractOne major deficiency of state-of-the-art galaxy formation models consists in their inability of capturing the observed galaxy downsizing trend significantly over-estimating the number density of low-mass galaxies, in particular at high redshifts. Employing an enhanced galaxy formation model with a full chemical enrichment scheme (DeLucia et al., 2014), we present an improved model for stellar feedback (based on parametrizations from cosmological zoom simulations), in which strong gas outflows occur due to bursty star formation at high z, while star formation is mainly “quiescent” not causing any significant outflows anymore at low z. Due to the stronger gas outflows at high z, early star formation is strongly delayed towards later times. This helps to sufficiently detach the evolution of galaxy growth from the hiearchical dark matter assembly resulting in a fairly good agreement with the evolution of the observed stellar mass function (SMF, see Fig. 1). With our new feedback scheme, we can also successfully reproduce many other observational constraints, such as the metallicity content, the cold gas fractions or the quiescent galaxy fractions at both low and high redshifts. The resulting new-generation galaxy catalogues (Hirschmann et al., in prep) based on that model are expected to significantly contribute to the interpretation of current and up-coming large-scale surveys (HST, JWST, Euclid). This will, in turn, provide a rapid verification and refinement of our modeling.

2019 ◽  
Vol 627 ◽  
pp. A131 ◽  
Author(s):  
M. Cousin ◽  
P. Guillard ◽  
M. D. Lehnert

Context. Star formation in galaxies is inefficient, and understanding how star formation is regulated in galaxies is one of the most fundamental challenges of contemporary astrophysics. Radiative cooling, feedback from supernovae and active galactic nuclei (AGN), and large-scale dynamics and dissipation of turbulent energy act over various time and spatial scales and all regulate star formation in a complex gas cycle. Aims. This paper presents the physics implemented in a new semi-analytical model of galaxy formation and evolution called the Galaxy Assembler from dark-matter Simulation (G.A.S.). Methods. The fundamental underpinning of our new model is the development of a multiphase interstellar medium (ISM) in which energy produced by supernovae and AGN maintains an equilibrium between a diffuse, hot, and stable gas and a cooler, clumpy, and low-volume filling factor gas. The hot gas is susceptible to thermal and dynamical instabilities. We include a description of how turbulence leads to the formation of giant molecular clouds through an inertial turbulent energy cascade, assuming a constant kinetic energy transfer per unit volume. We explicitly modelled the evolution of the velocity dispersion at different scales of the cascade and accounted for thermal instabilities in the hot halo gas. Thermal instabilities effectively reduce the impact of radiative cooling and moderates accretion rates onto galaxies, and in particular, for those residing in massive haloes. Results. We show that rapid and multiple exchanges between diffuse and unstable gas phases strongly regulates star formation rates in galaxies because only a small fraction of the unstable gas is forming stars. We checked that the characteristic timescales describing the gas cycle, gas depletion timescale, and star-forming laws at different scales are in good agreement with observations. For high-mass haloes and galaxies, cooling is naturally regulated by the growth of thermal instabilities, so we do not need to implement strong AGN feedback in this model. Our results are also in good agreement with the observed stellar mass function from z ≃ 6.0 to z ≃ 0.5. Conclusion. Our model offers the flexibility to test the impact of various physical processes on the regulation of star formation on a representative population of galaxies across cosmic times. Thermal instabilities and the cascade of turbulent energy in the dense gas phase introduce a delay between gas accretion and star formation, which keeps galaxy growth inefficient in the early Universe. The main results presented in this paper, such as stellar mass functions, are available in the GALAKSIENN library.


2020 ◽  
Vol 499 (2) ◽  
pp. 2401-2415
Author(s):  
A C Trapp ◽  
Steven R Furlanetto

ABSTRACT Cosmic variance is the intrinsic scatter in the number density of galaxies due to fluctuations in the large-scale dark matter density field. In this work, we present a simple analytic model of cosmic variance in the high-redshift Universe (z ∼ 5–15). We assume that galaxies grow according to the evolution of the halo mass function, which we allow to vary with large-scale environment. Our model produces a reasonable match to the observed ultraviolet (UV) luminosity functions in this era by regulating star formation through stellar feedback and assuming that the UV luminosity function is dominated by recent star formation. We find that cosmic variance in the UV luminosity function is dominated by the variance in the underlying dark matter halo population, and not by differences in halo accretion or the specifics of our stellar feedback model. We also find that cosmic variance dominates over Poisson noise for future high-z surveys except for the brightest sources or at very high redshifts (z ≳ 12). We provide a linear approximation of cosmic variance for a variety of redshifts, magnitudes, and survey areas through the public python package galcv. Finally, we introduce a new method for incorporating priors on cosmic variance into estimates of the galaxy luminosity function and demonstrate that it significantly improves constraints on that important observable.


2016 ◽  
Vol 11 (S321) ◽  
pp. 99-101 ◽  
Author(s):  
Gerhard Hensler ◽  
Patrick Steyrleithner ◽  
Simone Recchi

AbstractDue to their low masses dwarf galaxies experience low star-formation rates resulting in stellar cluster masses insufficient to fill the initial mass function (IMF) to the uppermost mass. Numerical simulations usually do not account for the completeness of the IMF, but treat a filed IMF by numbers, masses, and stellar feedback by fractions. To ensure that only entire stars are formed, we consider an IMF filled from the lower-mass regime and truncated where at least one entire massive star is formed.By 3D simulations we investigate the effects of two possible IMFs on the evolution of dwarf galaxies: filled vs. truncated IMF. For the truncated IMF the star-formation self-regulation is suppressed, while the energy release by typeII supernovae is larger, both compared to the filled IMF. Moreover, the abundance ratios of particular elements yielded from massive and intermediate-mass stars differ significantly between the two IMF distributions.


2020 ◽  
Vol 494 (1) ◽  
pp. 1002-1017 ◽  
Author(s):  
L Y Aaron Yung ◽  
Rachel S Somerville ◽  
Gergö Popping ◽  
Steven L Finkelstein

ABSTRACT The James Webb Space Telescope is expected to enable transformational progress in studying galaxy populations in the very early Universe, during the epoch of reionization. A critical parameter for understanding the sources that reionized the Universe is the Lyman-continuum production efficiency, ξion, defined as the rate of production of ionizing photons divided by the intrinsic UV luminosity. In this work, we combine self-consistent star formation and chemical enrichment histories predicted by semi-analytic models of galaxy formation with stellar population synthesis (SPS) models to predict the expected dependence of ξion on galaxy properties and cosmic epoch from z = 4–10. We then explore the sensitivity of the production rate of ionizing photons, $\dot{N}_\text{ion}$, to the choice of SPS model and the treatment of stellar feedback in our galaxy formation model. We compare our results to those of other simulations, constraints from empirical models, and observations. We find that adopting SPS models that include binary stars predict about a factor of 2 more ionizing radiation than models that only assume single stellar populations. We find that UV-faint, low-mass galaxies have values of ξion about 0.25 dex higher than those of more massive galaxies, but find weak evolution with cosmic time, about 0.2 dex from z ∼ 12–4 at fixed rest-UV luminosity. We provide predictions of $\dot{N}_\text{ion}$ as a function of Mh and a number of other galaxy properties. All results presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.


2019 ◽  
Vol 490 (4) ◽  
pp. 5375-5389 ◽  
Author(s):  
I Koutsouridou ◽  
A Cattaneo

ABSTRACT The difference in stellar metallicity between red and blue galaxies with the same mass constrains the time-scale over which red galaxies ceased to form stars. Here we investigate this constraint with the galics 2.0 semi-analytic model of galaxy formation. The advantage of this approach is that the time of pericentric passages for satellite galaxies and the mass-loading factor for galactic winds are not free parameters of the chemical evolution model. The former is determined by the N-body simulation used to construct the merger trees, the latter by the requirement that galics 2.0 should reproduce the stellar mass function of galaxies. When we compare our theoretical predictions with observations, we find that galics 2.0 can reproduce the observed metallicity difference only if quenching is preceded by a burst of star formation, which contributes to the chemical enrichment of the stellar population. Physically, this burst can be explained as tidally induced star formation or as an effect of ram pressure, which not only strips gas from galaxies but also compresses it, accelerating its conversion into stars.


2019 ◽  
Vol 489 (1) ◽  
pp. 487-496 ◽  
Author(s):  
Boyan K Stoychev ◽  
Keri L Dixon ◽  
Andrea V Macciò ◽  
Marvin Blank ◽  
Aaron A Dutton

ABSTRACT We use 38 high-resolution simulations of galaxy formation between redshift 10 and 5 to study the impact of a 3 keV warm dark matter (WDM) candidate on the high-redshift Universe. We focus our attention on the stellar mass function and the global star formation rate and consider the consequences for reionization, namely the neutral hydrogen fraction evolution and the electron scattering optical depth. We find that three different effects contribute to differentiate warm and cold dark matter (CDM) predictions: WDM suppresses the number of haloes with mass less than few 109 M⊙; at a fixed halo mass, WDM produces fewer stars than CDM, and finally at halo masses below 109 M⊙, WDM has a larger fraction of dark haloes than CDM post-reionization. These three effects combine to produce a lower stellar mass function in WDM for galaxies with stellar masses at and below 107 M⊙. For z > 7, the global star formation density is lower by a factor of two in the WDM scenario, and for a fixed escape fraction, the fraction of neutral hydrogen is higher by 0.3 at z ∼ 6. This latter quantity can be partially reconciled with CDM and observations only by increasing the escape fraction from 23 per cent to 34 per cent. Overall, our study shows that galaxy formation simulations at high redshift are a key tool to differentiate between dark matter candidates given a model for baryonic physics.


2021 ◽  
Vol 502 (4) ◽  
pp. 5417-5437
Author(s):  
Matthew C Smith

ABSTRACT Galaxy formation simulations frequently use initial mass function (IMF) averaged feedback prescriptions, where star particles are assumed to represent single stellar populations that fully sample the IMF. This approximation breaks down at high mass resolution, where stochastic variations in stellar populations become important. We discuss various schemes to populate star particles with stellar masses explicitly sampled from the IMF. We use Monte Carlo numerical experiments to examine the ability of the schemes to reproduce an input IMF in an unbiased manner while conserving mass. We present our preferred scheme which can easily be added to pre-existing star formation prescriptions. We then carry out a series of high-resolution isolated simulations of dwarf galaxies with supernovae (SNe), photoionization, and photoelectric heating to compare the differences between using IMF averaged feedback and explicitly sampling the IMF. We find that if SNe are the only form of feedback, triggering individual SNe from IMF averaged rates gives identical results to IMF sampling. However, we find that photoionization is more effective at regulating star formation when IMF averaged rates are used, creating more, smaller H ii regions than the rare, bright sources produced by IMF sampling. We note that the increased efficiency of the IMF averaged feedback versus IMF sampling is not necessarily a general trend and may be reversed depending on feedback channel, resolution and other details. However, IMF sampling is always the more physically motivated approach. We conservatively suggest that it should be used for star particles less massive than $\sim 500\, \mathrm{M_\odot }$.


2020 ◽  
Vol 497 (1) ◽  
pp. 279-301
Author(s):  
A Cattaneo ◽  
I Koutsouridou ◽  
E Tollet ◽  
J Devriendt ◽  
Y Dubois

ABSTRACT Dekel & Birnboim proposed that the mass-scale that separates late-type and early-type galaxies is linked to the critical halo mass $M_{\rm vir}^{\rm crit}$ for the propagation of a stable shock and showed that they could reproduce the observed bimodality scale for plausible values of the metallicity of the accreted gas Zaccr and the shock radius rs. Here, we take their analysis one step further and present a new semianalytic model that computes rs from first principles. This advancement allows us to compute $M_{\rm vir}^{\rm crit}$ individually for each halo. Separating cold-mode and hot-mode accretion has little effect on the final galaxy masses if feedback does not preferentially couple to the hot gas. We also present an improved model for stellar feedback where ${\sim }70{{\ \rm per\ cent}}$ of the wind mass is in a cold galactic fountain with a shorter reaccretion time-scale at high masses. The latter is the key mechanism that allows us to reproduce the low-mass end of the mass function of galaxies over the entire redshift range 0 < z < 2.5. Cooling must be mitigated to avoid overpredicting the number density of galaxies with stellar mass $M_{\rm stars}\gt 10^{11}\, {\rm M}_\odot$ but is important to form intermediate-mass galaxies. At $M_{\rm vir}\gt 3\times 10^{11}\, {\rm M}_\odot$, cold accretion is more important at high z, where gas is accreted from smaller solid angles, but this is not true at lower masses because high-z filaments have lower metallicities. Our predictions are consistent with the observed metallicity evolution of the intergalactic medium at 0 < z < 5.


Author(s):  
Mélanie Habouzit ◽  
Yuan Li ◽  
Rachel S Somerville ◽  
Shy Genel ◽  
Annalisa Pillepich ◽  
...  

Abstract The past decade has seen significant progress in understanding galaxy formation and evolution using large-scale cosmological simulations. While these simulations produce galaxies in overall good agreement with observations, they employ different sub-grid models for galaxies and supermassive black holes (BHs). We investigate the impact of the sub-grid models on the BH mass properties of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations, focusing on the MBH − M⋆ relation and the BH mass function. All simulations predict tight MBH − M⋆ relations, and struggle to produce BHs of $M_{\rm BH}\leqslant 10^{7.5}\, \rm M_{\odot }$ in galaxies of $M_{\star }\sim 10^{10.5}-10^{11.5}\, \rm M_{\odot }$. While the time evolution of the mean MBH − M⋆ relation is mild ($\rm \Delta M_{\rm BH}\leqslant 1\, dex$ for 0 ≤ z ≤ 5) for all the simulations, its linearity (shape) and normalization varies from simulation to simulation. The strength of SN feedback has a large impact on the linearity and time evolution for $M_{\star }\leqslant 10^{10.5}\, \rm M_{\odot }$. We find that the low-mass end is a good discriminant of the simulation models, and highlights the need for new observational constraints. At the high-mass end, strong AGN feedback can suppress the time evolution of the relation normalization. Compared with observations of the local Universe, we find an excess of BHs with $M_{\rm BH}\geqslant 10^{9}\, \rm M_{\odot }$ in most of the simulations. The BH mass function is dominated by efficiently accreting BHs (log10 fEdd ≥ −2) at high redshifts, and transitions progressively from the high-mass to the low-mass end to be governed by inactive BHs. The transition time and the contribution of active BHs are different among the simulations, and can be used to evaluate models against observations.


2019 ◽  
Vol 488 (1) ◽  
pp. 419-437 ◽  
Author(s):  
Xiaohan Wu ◽  
Rahul Kannan ◽  
Federico Marinacci ◽  
Mark Vogelsberger ◽  
Lars Hernquist

Abstract We present self-consistent radiation hydrodynamic simulations of hydrogen reionization performed with arepo-rt complemented by a state-of-the-art galaxy formation model. We examine how photoheating feedback, due to reionization, shapes the galaxies properties. Our fiducial model completes reionization by z ≈ 6 and matches observations of the Ly α forest, the cosmic microwave background electron scattering optical depth, the high-redshift ultraviolet (UV) luminosity function, and stellar mass function. Contrary to previous works, photoheating suppresses star formation rates by more than $50{{\ \rm per\ cent}}$ only in haloes less massive than ∼108.4 M⊙ (∼108.8 M⊙) at z = 6 (z = 5), suggesting inefficient photoheating feedback from photons within galaxies. The use of a uniform UV background that heats up the gas at z ≈ 10.7 generates an earlier onset of suppression of star formation compared to our fiducial model. This discrepancy can be mitigated by adopting a UV background model with a more realistic reionization history. In the absence of stellar feedback, photoheating alone is only able to quench haloes less massive than ∼109 M⊙ at z ≳ 5, implying that photoheating feedback is sub-dominant in regulating star formation. In addition, stellar feedback, implemented as a non-local galactic wind scheme in the simulations, weakens the strength of photoheating feedback by reducing the amount of stellar sources. Most importantly, photoheating does not leave observable imprints in the UV luminosity function, stellar mass function, or the cosmic star formation rate density. The feasibility of using these observables to detect imprints of reionization therefore requires further investigation.


Sign in / Sign up

Export Citation Format

Share Document