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: Bursty Star Formation in Dwarfs is Sensitive to Numerical Choices in Supernova Feedback Models
Abstract Simulations of galaxy formation are mostly unable to resolve the energy-conserving phase of individual supernova events, having to resort to subgrid models to distribute the energy and momentum resulting from stellar feedback. However, the properties of these simulated galaxies, including the morphology, stellar mass formed, and the burstiness of the star formation history, are highly sensitive to the numerical choices adopted in these subgrid models. Using the SMUGGLE stellar feedback model, we carry out idealized simulations of anMvir∼ 1010Mdwarf galaxy, a regime where most simulation codes predict significant burstiness in star formation, resulting in strong gas flows that lead to the formation of dark matter cores. We find that by varying only the directional distribution of momentum imparted from supernovae to the surrounding gas, while holding the total momentum per supernova constant, bursty star formation may be amplified or completely suppressed, and the total stellar mass formed can vary by as much as a factor of ∼3. In particular, when momentum is primarily directed perpendicular to the gas disk, less bursty and lower overall star formation rates result, yielding less gas turbulence, more disky morphologies, and a retention of cuspy dark matter density profiles. An improved understanding of the nonlinear coupling of stellar feedback into inhomogeneous gaseous media is thus needed to make robust predictions for stellar morphologies and dark matter core formation in dwarfs independent of uncertain numerical choices in the baryonic treatment.  more » « less
Award ID(s):
1945310 2107993
PAR ID:
10587239
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Publisher / Repository:
Zhang et al. 2024
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
975
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
229
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Due to their inability to self-regulate, ultrafaint dwarfs are sensitive to prescriptions in subgrid physics models that converge and regulate at higher masses. We use high-resolution cosmological simulations to compare the effect of bursty star formation histories (SFHs) on dwarf galaxy structure for two different subgrid supernova (SN) feedback models, superbubble and blastwave, in dwarf galaxies with stellar masses from 5000 <M*/M< 109. We find that in the “MARVEL-ous Dwarfs” suite both feedback models produce cored galaxies and reproduce observed scaling relations for luminosity, mass, and size. Our sample accurately predicts the average stellar metallicity at higher masses, however low-mass dwarfs are metal poor relative to observed galaxies in the Local Group. We show that continuous bursty star formation and the resulting stellar feedback are able to create dark matter (DM) cores in the higher dwarf galaxy mass regime, while the majority of ultrafaint and classical dwarfs retain cuspy central DM density profiles. We find that the effective core formation peaks atM*/Mhalo≃ 5 × 10−3for both feedback models. Both subgrid SN models yield bursty SFHs at higher masses; however, galaxies simulated with superbubble feedback reach maximum mean burstiness values at lower stellar mass fractions relative to blastwave feedback. As a result, core formation may be better predicted by stellar mass fraction than the burstiness of SFHs. 
    more » « less
  2. ABSTRACT We present a suite of 16 high-resolution hydrodynamic simulations of an isolated dwarf galaxy (gaseous and stellar disc plus a stellar bulge) within an initially cuspy dark matter (DM) halo, including self-interactions between the DM particles; as well as stochastic star formation and subsequent supernova feedback (SNF), implemented using the stellar feedback model SMUGGLE. The simulations start from identical initial conditions, and we regulate the strength of DM self-interactions and SNF by systematically varying the self-interacting DM (SIDM) momentum transfer cross-section and the gas density threshold for star formation. The DM halo forms a constant density core of similar size and shape for several combinations of those two parameters. Haloes with cores that are formed due to SIDM (adiabatic cusp-core transformation) have velocity dispersion profiles that are closer to isothermal than those of haloes with cores that are formed due to SNF in simulations with bursty star formation (impulsive cusp-core transformation). Impulsive SNF can generate positive stellar age gradients and increase random motion in the gas at the centre of the galaxy. Simulated galaxies in haloes with cores that were formed adiabatically are spatially more extended, with stellar metallicity gradients that are shallower (at late times) than those of galaxies in other simulations. Such observable properties of the gas and the stars, which indicate either an adiabatic or an impulsive evolution of the gravitational potential, may be used to determine whether observed cores in DM haloes are formed through DM self-interactions or in response to impulsive SNF. 
    more » « less
  3. Abstract Bursty star formation—a key prediction for high-redshift galaxies from cosmological simulations explicitly resolving stellar feedback in the interstellar medium—has recently been observed to prevail among galaxies at redshiftz≳ 6. Line intensity mapping (LIM) of the 158μm [Cii] line as a star formation rate (SFR) indicator offers unique opportunities to tomographically constrain cosmic star formation at high redshift, in a way complementary to observations of individually detected galaxies. To understand the effects of bursty star formation on [Cii] LIM, which have remained unexplored in previous studies, we present an analytic modeling framework for high-zgalaxy formation and [Cii] LIM signals that accounts for bursty star formation histories induced by delayed supernova feedback. We use it to explore and characterize how bursty star formation can impact and thus complicate the interpretation of the [Cii] luminosity function and power spectrum. Our simple analytic model indicates that bursty star formation mainly affects low-mass galaxies by boosting their average SFR and [Cii] luminosity, and in the [Cii] power spectrum it can create a substantial excess in the large-scale clustering term. This distortion results in a power spectrum shape that cannot be explained by invoking a mass-independent logarithmic scatter. We conclude that burstiness must be accounted for when modeling and analyzing [Cii] data sets from the early Universe, and that in the extreme, the signature of burstiness may be detectable with first-generation experiments such as TIME, CONCERTO, and CCAT-DSS. 
    more » « less
  4. The formation of compact high-redshift star-forming clumps, along with the physical processes driving their evolution and their potential connection to present-day globular clusters (GCs), are key open questions in studies of galaxy formation. In this work, we aim to shed light on these aspects using the SImulating the Environment where Globular clusters Emerged (SIEGE) project, a suite of cosmological zoom-in simulations with subparsec resolution that is specifically designed to investigate the physical conditions behind the origin of compact stellar systems in high-redshift environments. The simulations analyzed in this study are focused on a dwarf galaxy with a virial mass of a few 109Matz= 6.14, where the spatial resolution reaches 0.3 pc h−1. Individual stars are formed directly by sampling the initial mass function, with a 100% star formation efficiency. This setup is designed to explore the impact of a high star formation efficiency under high-redshift conditions. The simulation reveals the emergence of numerous stellar clumps with sizes of 1–3 pc, stellar surface densities up to almost 104Mpc−2, and masses predominantly spanning 103Mto several 104M, with a few reaching 105Mand up to 106M. All clumps form during intense, short bursts of star formation lasting less than a megayear, without noticeable signs of second peaks of star formation or accretion, often with negligible dark matter content (i.e., dark-to-stellar mass ratios below 1 within three times their effective radii). We measured a clear correlation between mass and size, with a clump mass function described by a power law with a slope of −2. Star formation conditions in the simulation reveal a behaviour that is similar to that of a feedback-free starburst scenario, where dense clumps form due to inefficient stellar feedback over small timescales. Notably, some clumps exhibit properties that closely resemble those of present-day globular clusters, highlighting their potential evolutionary connection. 
    more » « less
  5. Context.The detection of supermassive black holes (SMBHs) in high-redshift luminous quasars may require a phase of rapid accretion, and as a precondition, substantial gas influx toward seed black holes (BHs) from kiloparsec or parsec scales. Our previous research demonstrated the plausibility of such gas supply for BH seeds within star-forming giant molecular clouds (GMCs) with high surface density (∼104 M pc−2), facilitating “hyper-Eddington” accretion via efficient feeding by dense clumps, which are driven by turbulence and stellar feedback. Aims.This article presents an investigation of the impacts of feedback from accreting BHs on this process, including radiation, mechanical jets, and highly relativistic cosmic rays. Methods.We ran a suite of numerical simulations to explore diverse parameter spaces of BH feedback, including the subgrid accretion model, feedback energy efficiency, mass loading factor, and initial metallicity. Results.Using radiative feedback models inferred from the slim disk, we find that hyper-Eddington accretion is still achievable, yielding BH bolometric luminosities of as high as 1041 − 1044 erg/s, depending on the GMC properties and specific feedback model assumed. We find that the maximum possible mass growth of seed BHs (ΔMmaxBH) is regulated by the momentum-deposition rate from BH feedback,ṗfeedback/(ṀBHc), which leads to an analytic scaling that agrees well with simulations. This scenario predicts the rapid formation of ∼104Mintermediate-massive BHs (IMBHs) from stellar-mass BHs within ∼1 Myr. Furthermore, we examine the impacts of subgrid accretion models and how BH feedback may influence star formation within these cloud complexes. 
    more » « less