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.


This content will become publicly available on November 1, 2025

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
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. 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
  5. Abstract We present a new suite of numerical simulations of the star-forming interstellar medium (ISM) in galactic disks using the TIGRESS-NCR framework. Distinctive aspects of our simulation suite are (1) sophisticated and comprehensive numerical treatments of essential physical processes including magnetohydrodynamics, self-gravity, and galactic differential rotation, as well as photochemistry, cooling, and heating coupled with direct ray-tracing UV radiation transfer and resolved supernova feedback and (2) wide parameter coverage including the variation in metallicity over Z Z / Z 0.1 - 3 , gas surface density Σgas∼ 5–150Mpc−2, and stellar surface density Σstar∼ 1–50Mpc−2. The range of emergent star formation rate surface density is ΣSFR∼ 10−4–0.5Mkpc−2yr−1, and ISM total midplane pressure isPtot/kB= 103–106cm−3K, withPtotequal to the ISM weight W . For given Σgasand Σstar, we find Σ SFR Z 0.3 . We provide an interpretation based on the pressure-regulated feedback-modulated (PRFM) star formation theory. The total midplane pressure consists of thermal, turbulent, and magnetic stresses. We characterize feedback modulation in terms of the yield ϒ, defined as the ratio of each stress to ΣSFR. The thermal feedback yield varies sensitively with both weight and metallicity as ϒ th W 0.46 Z 0.53 , while the combined turbulent and magnetic feedback yield shows weaker dependence ϒ turb + mag W 0.22 Z 0.18 . The reduction in ΣSFRat low metallicity is due mainly to enhanced thermal feedback yield, resulting from reduced attenuation of UV radiation. With the metallicity-dependent calibrations we provide, PRFM theory can be used for a new subgrid star formation prescription in cosmological simulations where the ISM is unresolved. 
    more » « less