Attention:The NSF Public Access Repository (PAR) system and access will be unavailable from 5:00 PM ET until 8:00 PM ET on Friday, September 11 due to maintenance. We apologize for the inconvenience.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


  • NSF Public Access
  • Search Results
  • RIGEL: Simulating dwarf galaxies at solar mass resolution with radiative transfer and feedback from individual massive stars
Title: RIGEL: Simulating dwarf galaxies at solar mass resolution with radiative transfer and feedback from individual massive stars
Context.Feedback from stars in the form of radiation, stellar winds, and supernovae is crucial to regulating the star formation activity of galaxies. Dwarf galaxies are especially susceptible to these processes, making them an ideal test bed for studying the effects of stellar feedback in detail. Recent numerical models have aimed to resolve the interstellar medium (ISM) in dwarf galaxies with a very high resolution of several solar masses. However, when it comes to modeling the radiative feedback from stars, many models opt for simplified approaches instead of explicitly solving radiative transfer (RT) because of the computational complexity involved. Aims.We introduce the Realistic ISM modeling in Galaxy Evolution and Lifecycles (RIGEL) model, a novel framework to self-consistently model the effects of stellar feedback in the multiphase ISM of dwarf galaxies with explicit RT on a star-by-star basis. Methods.The RIGEL model integrates detailed implementations of feedback from individual massive stars into the state-of-the-art radiation-hydrodynamics code,AREPO-RT. It forms individual massive stars from the resolved multiphase ISM by sampling the initial mass function and tracks their evolution individually. The lifetimes, photon production rates, mass-loss rates, and wind velocities of these stars are determined by their initial masses and metallicities based on a library that incorporates a variety of stellar models. The RT equations are solved explicitly in seven spectral bins accounting for the infrared to He IIionizing bands, using a moment-base scheme with the M1 closure relation. The thermochemistry model tracks the nonequilibrium H, He chemistry as well as the equilibrium abundance of C I, C II, O I, O II, and CO in the irradiated ISM to capture the thermodynamics of all ISM phases, from cold molecular gas to hot ionized gas. Results.We evaluated the performance of the RIGEL model using 1 M⊙resolution simulations of isolated dwarf galaxies. We found that the star formation rate (SFR) and interstellar radiation field (ISRF) show strong positive correlations with the metallicity of the galaxy. Photoionization and photoheating can reduce the SFR by an order of magnitude by removing the available cold, dense gas fuel for star formation. The presence of ISRF also significantly changes the thermal structure of the ISM. Radiative feedback occurs immediately after the birth of massive stars and rapidly disperses the molecular clouds within 1 Myr. As a consequence, radiative feedback reduces the age spread of star clusters to less than 2 Myr, prohibits the formation of massive star clusters, and shapes the cluster initial mass function to a steep power-law form with a slope of ∼ − 2. The mass-loading factor (measured atz = 1 kpc) of the fiducial galaxy has a median ofηM ∼ 50, while turning off radiative feedback reduces this factor by an order of magnitude. Conclusions.We demonstrate that RIGEL effectively captures the nonlinear coupling of early radiative feedback and supernova feedback in the multiphase ISM of dwarf galaxies. This novel framework enables the utilization of a comprehensive stellar feedback and ISM model in cosmological simulations of dwarf galaxies and various galactic environments spanning a wide dynamic range in both space and time.  more » « less
Award ID(s):
2108470
PAR ID:
10632676
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
EDP Sciences
Date Published:
Journal Name:
Astronomy & Astrophysics
Volume:
691
ISSN:
0004-6361
Page Range / eLocation ID:
A231
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ; ; ; ; ; et al (, Astronomy & Astrophysics)
    Context. The fundamental process of star formation in galaxies involves the intricate interplay between the fueling of star formation via molecular gas and the feedback from recently formed massive stars that can, in turn, hinder the conversion of gas into stars. This process, by which galaxies evolve, is also closely connected to the intrinsic properties of the interstellar medium (ISM), such as structure, density, pressure, and metallicity. Aims. To study the role that different molecular and atomic phases of the ISM play in star formation, and to characterize their physical conditions, we zoom into our nearest neighboring galaxy, the Large Magellanic Cloud (LMC; 50 kpc), the most convenient laboratory in which to study the effects of the lower metal abundance on the properties of the ISM. The LMC offers a view of the ISM and star formation conditions in a low-metallicity (Z~ 0.5 Z⊙) environment similar, in that regard, to the epoch of the peak of star formation in the earlier Universe (z~ 1.5). Following up on studies carried out at galactic scales in low-Z galaxies, we present an unprecedentedly detailed analysis of well-known star-forming regions (SFRs) at a spatial resolution of a few parsecs. Methods. We mapped a 610pc× 260pc region in the LMC molecular ridge in [C II]λ158 µm and the [O III]λ88 µm using the FIFI-LS instrument on the SOFIA telescope. We compared the data with the distribution of the CO(2−1) emission from ALMA, the modeled total infrared luminosity, and the Spitzer/MIPS 24 µm continuum and Hα. Results. We present new large maps of [CII] and [OIII] and perform a first comparison with CO(2−1) line and LTIR emission. We also provide a detailed description of the observing strategy with SOFIA/FIFI-LS and the data reduction process. Conclusions. We find that [CII] and [OIII] emission is associated with the SFRs in the molecular ridge, but also extends throughout the mapped region, and is not obviously associated with ongoing star formation. The CO emission is clumpier than the [C II] emission and we find plentiful [C II] present where there is little CO emission, possibly holding important implications for “CO-dark” gas. We find a clear trend of the L[C II]/LTIRratio decreasing with increasing LTIRin the full range. This suggests a strong link between the “[C II]-deficit” and the local physical conditions instead of global properties. 
    more » « less
  2. ; ; ; ; ; ; ; ; (, Astronomy & Astrophysics)
    The chemical abundances of stars in galaxies are a fossil record of the star formation and stellar evolution processes that regulate galaxy formation, including the stellar initial mass function, the fraction and timing of type Ia supernovae (SNeIa), and nucleosynthesis inside massive stars. In this paper, we systematically explore uncertainties associated with modeling chemical enrichment in dwarf galaxies. We repeatedly simulate a singleEDGE-INFERNOdwarf (M★ ≈ 105 M⊙), varying the chemical yields of massive stars, the timing and yields of SNeIa, and the intrinsic stochasticity that arises from sampling individual stars and galaxy formation chaoticity. All simulations are high-resolution (3.6 pc), cosmological zoom-in hydrodynamical simulations that track the stellar evolution of all individual stars with masses of > 0.5 M⊙. We find that SNeIa make significant contributions to the iron content of low-mass, reionization-limited galaxies, with possible variations in mean abundance ratios and [Fe/H] related to minor changes in their evolutionary timescales. In contrast, different massive star yields, accounting (or not) for stellar rotation, result in mean abundance variations comparable to those arising from stochasticity, with the possible exception of extremely rapidly rotating stars. Nonetheless, massive stars significantly affect the shape of abundance trends with [Fe/H], for example, through the existence (or not) of a bimodality in the [X/Fe]–[Fe/H] planes, particularly in [Al/Fe]. Finally, we find that the variance arising from random sampling severely limits the interpretation of single galaxies. Our analysis showcases the power of star-by-star cosmological models to unpick how both systematic uncertainties (e.g., assumptions in low-metallicity chemical enrichment) and statistical uncertainties (e.g., averaging over enough galaxies and stars within a galaxy) affect the interpretation of chemical observables in ultra-faint dwarf galaxies. 
    more » « less
  3. ; ; ; ; ; ; ; ; (, The Astrophysical Journal)
    Abstract Feedback from supermassive black holes is believed to be a critical driver of the observed color bimodality of galaxies above the Milky Way mass scale. Active galactic nuclei (AGN) feedback has been modeled in many galaxy formation simulations, but most implementations have involved simplified prescriptions or a coarse-grained interstellar medium (ISM). We present the first set of Feedback In Realistic Environments (FIRE)-3 cosmological zoom-in simulations with AGN feedback evolved toz∼ 0, examining the impact of AGN feedback on a set of galaxies with halos in the mass range 1012–1013M⊙. These simulations combine detailed stellar and ISM physics with multichannel AGN feedback including radiative feedback, mechanical outflows, and, in some simulations, cosmic rays (CRs). We find that massive (>L*) galaxies in these simulations can match local scaling relations including the stellar mass–halo mass relation and theMBH–σrelation; in the stronger model with CRs, they also match the size–mass relation and the Faber–Jackson relation. Many of the massive galaxies in the simulations with AGN feedback have quenched star formation and elliptical morphologies, in qualitative agreement with observations. In contrast, simulations at the massive end without AGN feedback produce galaxies that are too massive and form stars too rapidly, are order-of-magnitude too compact, and have velocity dispersions well above Faber–Jackson. Despite these successes, the AGN models analyzed do not produce uniformly realistic galaxies when the feedback parameters are held constant: While the stronger model produces the most realistic massive galaxies, it tends to overquench the lower-mass galaxies. This indicates that further refinements of the AGN modeling are needed. 
    more » « less
  4. ; ; ; ; ; ; ; ; ; ; et al (, Astronomy & Astrophysics)
    Characterizing the physical conditions (density, temperature, ionization state, metallicity, etc) of the interstellar medium is critical to improving our understanding of the formation and evolution of galaxies. In this work, we present a multi-line study of the interstellar medium in the host galaxy of a quasar atz ≈ 6.4, that is, when the universe was 840 Myr old. This galaxy is one of the most active and massive objects emerging from the dark ages and therefore represents a benchmark for models of the early formation of massive galaxies. We used the Atacama Large Millimeter Array to target an ensemble of tracers of ionized, neutral, and molecular gas, namely the following fine-structure lines: [O III] 88 μm, [N II] 122 μm, [C II] 158 μm, and [C I] 370 μm – as well as the rotational transitions of CO(7–6), CO(15–14), CO(16–15), and CO(19–18); OH 163.1 μm and 163.4 μm; along with H2O 3(0,3)–2(1,2), 3(3,1)–4(0,4), 3(3,1)–3(2,2), 4(0,4)–3(1,3), and 4(3,2)–4(2,3). All the targeted fine-structure lines were detected, along with half of the targeted molecular transitions. By combining the associated line luminosities with the constraints on the dust temperature from the underlying continuum emission and predictions from photoionization models of the interstellar medium, we find that the ionized phase accounts for about one-third of the total gaseous mass budget and is responsible for half of the total [C II] emission. This phase is characterized by a high density (n ∼ 180 cm−3) that typical of HII regions. The spectral energy distribution of the photoionizing radiation is comparable to that emitted by B-type stars. Star formation also appears to be driving the excitation of the molecular medium. We find marginal evidence for outflow-related shocks in the dense molecular phase, but not in other gas phases. This study showcases the power of multi-line investigations in unveiling the properties of the star-forming medium in galaxies at cosmic dawn. 
    more » « less
  5. ; ; ; ; ; ; ; ; ; ; et al (, Astronomy & Astrophysics)
    The mode of star formation that results in the formation of globular clusters and young massive clusters is difficult to constrain through observations. We present models of massive star cluster formation using the TORCHframework, which uses the Astrophysical MUltipurpose Software Environment (AMUSE) to couple distinct multi-physics codes that handle star formation, stellar evolution and dynamics, radiative transfer, and magnetohydrodynamics. We upgraded TORCHby implementing the N-body code PETAR, thereby enabling TORCHto handle massive clusters forming from 106M⊙clouds with ≥105individual stars. We present results from TORCHsimulations of star clusters forming from 104,  105, and 106M⊙turbulent spherical gas clouds (named M4, M5, M6) of radiusR= 11.7 pc. We find that star formation is highly efficient and becomes more so at a higher cloud mass and surface density. For M4, M5, and M6 with initial surface densities 2.325 × 101,2,3M⊙pc−2, after a free-fall time oftff= 6.7,2.1,0.67 Myr, we find that ∼30%, 40%, and 60% of the cloud mass has formed into stars, respectively. The end of simulation-integrated star formation efficiencies for M4, M5, and M6 areϵ⋆ = M⋆/Mcloud = 36%, 65%, and 85%. Observations of nearby clusters similar in mass and size to M4 have instantaneous star formation efficiencies ofϵinst ≤ 30%, which is slightly lower than the integrated star formation efficiency of M4. The M5 and M6 models represent a different regime of cluster formation that is more appropriate for the conditions in starburst galaxies and gas-rich galaxies at high redshift, and that leads to a significantly higher efficiency of star formation. We argue that young massive clusters build up through short efficient bursts of star formation in regions that are sufficiently dense (Σ ≥ 102M⊙pc−2) and massive (Mcloud≥ 105M⊙). In such environments, stellar feedback from winds and radiation is not strong enough to counteract the gravity from gas and stars until a majority of the gas has formed into stars. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Text Encoded Conversion
  • XML
  • Save / Share this Record
  • Send to Email