skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Mergers, starbursts, and quenching in the simba simulation
ABSTRACT We use the simba cosmological galaxy formation simulation to investigate the relationship between major mergers ($$\lesssim$$4:1), starbursts, and galaxy quenching. Mergers are identified via sudden jumps in stellar mass M* well above that expected from in situ star formation, while quenching is defined as going from specific star formation rate (sSFR) $$\gt t_{\rm H}^{-1}$$ to $$\lt 0.2t_{\rm H}^{-1}$$, where tH is the Hubble time. At z ≈ 0–3, mergers show ∼2–3× higher SFR than a mass-matched sample of star-forming galaxies, but globally represent $$\lesssim 1{{\ \rm per\ cent}}$$ of the cosmic SF budget. At low masses, the increase in SFR in mergers is mostly attributed to an increase in the H2 content, but for $$M_*\gtrsim 10^{10.5} \,\mathrm{ M}_{\odot }$$ mergers also show an elevated star formation efficiency suggesting denser gas within merging galaxies. The merger rate for star-forming galaxies shows a rapid increase with redshift, ∝(1 + z)3.5, but the quenching rate evolves much more slowly, ∝(1 + z)0.9; there are insufficient mergers to explain the quenching rate at $$z\lesssim 1.5$$. simba first quenches galaxies at $$z\gtrsim 3$$, with a number density in good agreement with observations. The quenching time-scales τq are strongly bimodal, with ‘slow’ quenchings (τq ∼ 0.1tH) dominating overall, but ‘fast’ quenchings (τq ∼ 0.01tH) dominating in M* ∼ 1010–1010.5 M$$\odot$$ galaxies, likely induced by simba’s jet-mode black hole feedback. The delay time distribution between mergers and quenching events suggests no physical connection to either fast or slow quenching. Hence, simba predicts that major mergers induce starbursts, but are unrelated to quenching in either fast or slow mode.  more » « less
Award ID(s):
1908137 1715206
PAR ID:
10121814
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
490
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
p. 2139-2154
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT We revisit the question of ‘hot mode’ versus ‘cold mode’ accretion on to galaxies using steady-state cooling flow solutions and idealized 3D hydrodynamic simulations. We demonstrate that for the hot accretion mode to exist, the cooling time is required to be longer than the free-fall time near the radius where the gas is rotationally supported, Rcirc, i.e. the existence of the hot mode depends on physical conditions at the galaxy scale rather than on physical conditions at the halo scale. When allowing for the depletion of the halo baryon fraction relative to the cosmic mean, the longer cooling times imply that a virialized gaseous halo may form in halo masses below the threshold of $$\sim 10^{12}\, {\rm M_{\odot }}$$ derived for baryon-complete haloes. We show that for any halo mass there is a maximum accretion rate for which the gas is virialized throughout the halo and can accrete via the hot mode of $${\dot{M}}_{\rm crit}\approx 0.7(v_{\rm c}/100\, \rm km\ s^{-1})^{5.4}(R_{\rm circ}/10\, {\rm kpc})(Z/\, {\rm Z_{\odot }})^{-0.9}\, {\rm M_{\odot }}\, {\rm yr}^{-1}$$, where Z and vc are the metallicity and circular velocity measured at Rcirc. For accretion rates $$\gtrsim {\dot{M}}_{\rm crit}$$ the volume-filling gas phase can in principle be ‘transonic’ – virialized in the outer halo but cool and free-falling near the galaxy. We compare $${\dot{M}}_{\rm crit}$$ to the average star formation rate (SFR) in haloes at 0 < z < 10 implied by the stellar-mass–halo-mass relation. For a plausible metallicity evolution with redshift, we find that $${\rm SFR}\lesssim {\dot{M}}_{\rm crit}$$ at most masses and redshifts, suggesting that the SFR of galaxies could be primarily sustained by the hot mode in halo masses well below the classic threshold of $$\sim 10^{12}\, {\rm M_{\odot }}$$. 
    more » « less
  2. ABSTRACT The star formation and gas content of satellite galaxies around the Milky Way (MW) and Andromeda (M31) are depleted relative to more isolated galaxies in the Local Group (LG) at fixed stellar mass. We explore the environmental regulation of gas content and quenching of star formation in z = 0 galaxies at $$M_{*}=10^{5\!-\!10}\, \rm {M}_{\odot }$$ around 14 MW-mass hosts from the Feedback In Realistic Environments 2 (FIRE-2) simulations. Lower mass satellites ($$M_{*}\lesssim 10^7\, \rm {M}_{\odot }$$) are mostly quiescent and higher mass satellites ($$M_{*}\gtrsim 10^8\, \rm {M}_{\odot }$$) are mostly star forming, with intermediate-mass satellites ($$M_{*}\approx 10^{7\!-\!8}\, \rm {M}_{\odot }$$) split roughly equally between quiescent and star forming. Hosts with more gas in their circumgalactic medium have a higher quiescent fraction of massive satellites ($$M_{*}=10^{8\!-\!9}\, \rm {M}_{\odot }$$). We find no significant dependence on isolated versus paired (LG-like) host environments, and the quiescent fractions of satellites around MW-mass and Large Magellanic Cloud (LMC)-mass hosts from the FIRE-2 simulations are remarkably similar. Environmental effects that lead to quenching can also occur as pre-processing in low-mass groups prior to MW infall. Lower mass satellites typically quenched before MW infall as central galaxies or rapidly during infall into a low-mass group or a MW-mass galaxy. Most intermediate- to high-mass quiescent satellites have experienced ≥1–2 pericentre passages (≈2.5–5 Gyr) within a MW-mass halo. Most galaxies with $$M_{*}\gtrsim 10^{6.5}\, \rm {M}_{\odot }$$ did not quench before falling into a host, indicating a possible upper mass limit for isolated quenching. The simulations reproduce the average trend in the LG quiescent fraction across the full range of satellite stellar masses. Though the simulations are consistent with the Satellites Around Galactic Analogs (SAGA) survey’s quiescent fraction at $$M_{*}\gtrsim 10^8\, \rm {M}_{\odot }$$, they do not generally reproduce SAGA’s turnover at lower masses. 
    more » « less
  3. ABSTRACT The shape of the low-mass (faint) end of the galaxy stellar mass function (SMF) or ultraviolet luminosity function (UVLF) at $$z \gtrsim 6$$ is an open question for understanding which galaxies primarily drove cosmic reionization. Resolved photometry of Local Group low-mass galaxies allows us to reconstruct their star formation histories, stellar masses, and UV luminosities at early times, and this fossil record provides a powerful ‘near-far’ technique for studying the reionization-era SMF/UVLF, probing orders of magnitude lower in mass than direct HST/JWST observations. Using 882 low-mass ($$M_{\rm star}\lesssim 10^{9}\, \rm {M_\odot }$$) galaxies across 11 Milky Way (MW)- and Local Group-analogue environments from the FIRE-2 cosmological baryonic zoom-in simulations, we characterize their progenitors at $$z=6\!-\!9$$, the mergers/disruption of those progenitors over time, and how well their present-day fossil record traces the high-redshift SMF. A present-day galaxy with $$M_{\rm star}\sim 10^5\, \rm {M_\odot }$$ ($$\sim 10^9\, \rm {M_\odot }$$) had $$\approx 1$$ ($$\approx 30$$) progenitors at $$z\approx 7$$, and its main progenitor comprised $$\approx 100~{{\ \rm per\ cent}}$$ ($$\approx 10~{{\ \rm per\ cent}}$$) of the total stellar mass of all its progenitors at $$z\approx 7$$. We show that although only $$\sim 15~{{\ \rm per\ cent}}$$ of the early population of low-mass galaxies survives to present day, the fossil record of surviving Local Group galaxies accurately traces the low-mass slope of the SMF at $$z \sim 6 \!-\! 9$$. We find no obvious mass dependence to the mergers and accretion, and show that applying this reconstruction technique to just low-mass galaxies at $z = 0$ and not the MW/M31 hosts correctly recovers the slope of the SMF down to $$M_{\rm star} \sim 10^{4.5}\, \rm {{\rm M}_{\odot }}$$ at $$z \gtrsim 6$$. Thus, we validate the ‘near-far’ approach as an unbiased tool for probing low-mass reionization-era galaxies. 
    more » « less
  4. ABSTRACT Motivated by the early excess of bright galaxies seen by JWST, we run zoom-in cosmological simulations of a massive galaxy at Cosmic Dawn, in a halo of $$10^{11} {\rm M}_\odot$$ at $z = 9$, using the hydro-gravitational code ramses at an effective resolution $$\sim 10~{\rm pc}$$. We investigate physical mechanisms that enhance the star formation efficiencies (SFEs) at the high gas densities of the star-forming regions in this galaxy ($$\sim 3\times 10^3~{\rm cm^{-3}}$$, $$\sim 10^4~{\rm M}_\odot \,{\rm pc^{-2}}$$). Our fiducial star formation recipe uses a physically motivated, turbulence-based, multi-freefall model, avoiding ad hoc extrapolation from lower redshifts. By $z = 9$, our simulated galaxy is a clumpy, thick, rotating disc with a high stellar mass $$\sim 3\times 10^9~{\rm M}_\odot$$ and high star formation rate $$\sim 50~{\rm M}_\odot \,{\rm yr^{-1}}$$. The high gas density makes supernova (SN) feedback less efficient, producing a high local SFE $$\gtrsim 10~{{\ \rm per\ cent}}$$. The global SFE is set by feedback-driven outflows and only weakly correlated with the local SFE. Photoionization heating makes SN feedback more efficient, but the integrated SFE always remains high. Intense accretion at Cosmic Dawn seeds turbulence that reduces local SFE, but this only weakly affects the global SFE. The star formation histories of our simulated galaxies are similar to observed massive galaxies at Cosmic Dawn, despite our limited resolution. We set the stage for future simulations which treat radiation self-consistently and use a higher effective resolution $$\sim 1~{\rm pc}$$ that captures the physics of star-forming clouds. 
    more » « less
  5. ABSTRACT We study specific star formation rate (sSFR) and gas profiles of star-forming (SF) and green valley (GV) galaxies in the simba cosmological hydrodynamic simulation. SF galaxy half-light radii (Rhalf) at z = 0 and their evolution (∝(1 + z)−0.78) agree with observations. Passive galaxy Rhalf agree with observations at high redshift, but by z = 0 are too large, owing to numerical heating. We compare simbaz = 0 sSFR radial profiles for SF and GV galaxies to observations. simba shows strong central depressions in star formation rate (SFR), sSFR, and gas fraction in GV galaxies and massive SF systems, qualitatively as observed, owing to black hole X-ray feedback, which pushes central gas outwards. Turning off X-ray feedback leads to centrally peaked sSFR profiles as in other simulations. In conflict with observations, simba yields GV galaxies with strongly dropping sSFR profiles beyond ≳Rhalf, regardless of active galactic nucleus feedback. The central depression owes to lowering molecular gas content; the drop in the outskirts owes to reduced star formation efficiency. simba’s satellites have higher central sSFR and lower outskirts sSFR than centrals, in qualitative agreement with observations. At z = 2, simba does not show central depressions in massive SF galaxies, suggesting simba’s X-ray feedback should be more active at high-z. High-resolution tests indicate central sSFR suppression is not sensitive to numerical resolution. Reproducing the central sSFR depression in z = 0 GV galaxies represents a unique success of simba. The remaining discrepancies highlight the importance of SFR and gas profiles in constraining quenching mechanisms. 
    more » « less