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: 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 Mresolution 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. 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. ABSTRACT We investigate the spatial structure and evolution of star formation and the interstellar medium (ISM) in interacting galaxies. We use an extensive suite of parsec-scale galaxy-merger simulations (stellar mass ratio = 2.5:1), which employs the ‘Feedback In Realistic Environments-2’ model (fire-2). This framework resolves star formation, feedback processes, and the multiphase structure of the ISM. We focus on the galaxy-pair stages of interaction. We find that close encounters substantially augment cool (H i) and cold-dense (H2) gas budgets, elevating the formation of new stars as a result. This enhancement is centrally concentrated for the secondary galaxy, and more radially extended for the primary. This behaviour is weakly dependent on orbital geometry. We also find that galaxies with elevated global star formation rate (SFR) experience intense nuclear SFR enhancement, driven by high levels of either star formation efficiency (SFE) or available cold-dense gas fuel. Galaxies with suppressed global SFR also contain a nuclear cold-dense gas reservoir, but low SFE levels diminish SFR in the central region. Concretely, in the majority of cases, SFR enhancement in the central kiloparsec is fuel-driven (55 per cent for the secondary, 71 per cent for the primary) – while central SFR suppression is efficiency-driven (91 per cent for the secondary, 97 per cent for the primary). Our numerical predictions underscore the need of substantially larger, and/or merger-dedicated, spatially resolved galaxy surveys – capable of examining vast and diverse samples of interacting systems – coupled with multiwavelength campaigns aimed to capture their internal ISM structure. 
    more » « less
  3. 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. 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. ABSTRACT The formation and evolution of galaxies have proved sensitive to the inclusion of stellar feedback, which is therefore crucial to any successful galaxy model. We present INFERNO, a new model for hydrodynamic simulations of galaxies, which incorporates resolved stellar objects with star-by-star calculations of when and where the injection of enriched material, momentum, and energy takes place. INFERNO treats early stellar kinematics to include phenomena such as walkaway and runaway stars. We employ this innovative model on simulations of a dwarf galaxy and demonstrate that our physically motivated stellar feedback model can drive vigorous galactic winds. This is quantified by mass and metal loading factors in the range of 10–100, and an energy loading factor close to unity. Outflows are established close to the disc, are highly multiphase, spanning almost 8 orders of magnitude in temperature, and with a clear dichotomy between mass ejected in cold, slow-moving (T ≲ 5 × 104 K, v < 100 km s−1) gas and energy ejected in hot, fast-moving (T > 106 K, v > 100 km s−1) gas. In contrast to massive disc galaxies, we find a surprisingly weak impact of the early stellar kinematics, with runaway stars having little to no effect on our results, despite exploding in diffuse gas outside the dense star-forming gas, as well as outside the galactic disc entirely. We demonstrate that this weak impact in dwarf galaxies stems from a combination of strong feedback and a porous interstellar medium, which obscure any unique signatures that runaway stars provide. 
    more » « less