An official website of the United States government
Abstract Feedback from massive stars plays an important role in the formation of star clusters. Whether a very massive star is born early or late in the cluster formation timeline has profound implications for the star cluster formation and assembly processes. We carry out a controlled experiment to characterize the effects of early-forming massive stars on star cluster formation. We use the star formation software suiteTorch, combining self-gravitating magnetohydrodynamics, ray-tracing radiative transfer,N-body dynamics, and stellar feedback, to model four initially identical 104M⊙giant molecular clouds with a Gaussian density profile peaking at 521.5 cm−3. Using theTorchsoftware suite through theAMUSEframework, we modify three of the models, to ensure that the first star that forms is very massive (50, 70, and 100M⊙). Early-forming massive stars disrupt the natal gas structure, resulting in fast evacuation of the gas from the star-forming region. The star formation rate is suppressed, reducing the total mass of the stars formed. Our fiducial control model, without an early massive star, has a larger star formation rate and total efficiency by up to a factor of 3, and a higher average star formation efficiency per freefall time by up to a factor of 7. Early-forming massive stars promote the buildup of spatially separate and gravitationally unbound subclusters, while the control model forms a single massive cluster.
more »
« less
Triggered Population III star formation: the effect of H2 self-shielding
ABSTRACT The multiplicity of metal-free (Population III) stars may influence their feedback efficiency within their host dark matter haloes, affecting subsequent metal enrichment and the transition to galaxy formation. Radiative feedback from massive stars can trigger nearby star formation in dense self-shielded clouds. In model radiation self-shielding, the H2 column density must be accurately computed. In this study, we compare two local approximations based on the density gradient and Jeans length with a direct integration of column density along rays. After the primary massive star forms, we find that no secondary stars form for both the direct integration and density gradient approaches. The approximate method reduces the computation time by a factor of 2. The Jeans length approximation overestimates the H2 column density by a factor of 10, leading to five numerically enhanced self-shielded, star-forming clumps. We conclude that the density gradient approximation is sufficiently accurate for larger volume galaxy simulations, although one must still caution that the approximation cannot fully reproduce the result of direct integration.
more »
« less
- Award ID(s):
- 1835213
- PAR ID:
- 10398587
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 520
- Issue:
- 4
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 5077-5089
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT ‘Runaway stars’ might play a role in driving galactic outflows and enriching the circumgalactic medium with metals. To study this effect, we carry out high-resolution dwarf galaxy simulations that include velocity ‘kicks’ to massive stars above eigth solar masses. We consider two scenarios, one that adopts a power law velocity distribution for kick velocities, resulting in more stars with high-velocity kicks, and a more moderate scenario with a Maxwellian velocity distribution. We explicitly resolve the multiphase interstellar medium (ISM) and include non-equilibrium cooling and chemistry. We sample individual massive stars from an IMF and follow their radiation input and SN feedback (core-collapse) channel at the end of their lifetime. In the simulations with runaway stars, we add additional (natal) velocity kicks that mimic two- and three-body interactions that cannot be fully resolved in our simulations. We find that including runaway or ‘walkaway’ star scenarios impacts mass, metal, momentum, and energy outflows as well as the corresponding loading factors. The effect on the mass loading factor is small, but we find an increase in the metal loading by a factor of 1.5 to 2. The momentum loading increases by a factor of 1.5–2. The energy loading increases by roughly a factor of 5 when runaway stars are included. Additionally, we find that the overall level of star formation is increased in the models that include runaway stars. We conclude that the inclusion of runaway stars could have an impact on the global star formation and subsequent outflow properties of dwarf galaxies.more » « less
-
null (Ed.)ABSTRACT The formation of Population III (Pop III) stars is a critical step in the evolution of the early Universe. To understand how these stars affected their metal-enriched descendants, the details of how, why and where Pop III formation takes place needs to be determined. One of the processes that is assumed to greatly affect the formation of Pop III stars is the presence of a Lyman–Werner (LW) radiation background, that destroys H2, a necessary coolant in the creation of Pop III stars. Self-shielding can alleviate the effect the LW background has on the H2 within haloes. In this work, we perform a cosmological simulation to study the birthplaces of Pop III stars, using the adaptive mesh refinement code enzo. We investigate the distribution of host halo masses and its relationship to the LW background intensity. Compared to previous work, haloes form Pop III stars at much lower masses, up to a factor of a few, due to the inclusion of H2 self-shielding. We see no relationship between the LW intensity and host halo mass. Most haloes form multiple Pop III stars, with a median number of four, up to a maximum of 16, at the instance of Pop III formation. Our results suggest that Pop III star formation may be less affected by LW radiation feedback than previously thought and that Pop III multiple systems are common.more » « less
-
ABSTRACT The feedback loop between the galaxies producing the background radiation field for reionization and their growth is crucial, particularly for low-mass haloes. Despite this, the vast majority of galaxy formation studies employ a spatially uniform, time-varying reionizing background, with the majority of reionization studies employing galaxy formation models only required to work at high redshift. This paper uses the well-studied TNG galaxy formation model, calibrated at low redshift, coupled to the arepo-rt code, to self-consistently solve the coupled problems of galaxy evolution and reionization, evaluating the impact of patchy (and slow) reionization on early galaxies. thesan-hr is an extension of the thesan project to higher resolution (a factor of 50 increase, with a baryonic mass of mb ≈ 104 M⊙), to additionally enable the study of ‘mini-haloes’ with virial temperatures Tvir < 104 K. Comparing the self-consistent model to a uniform UV background, we show that galaxies in thesan-hr are predicted to be larger in physical extent (by a factor ∼2), less metal enriched (by ∼0.2 dex), and less abundant (by a factor ∼10 at M1500 = − 10) by z = 5. We show that differences in star formation and enrichment patterns lead to significantly different predictions for star formation in low mass haloes, low-metallicity star formation, and even the occupation fraction of haloes. We posit that cosmological galaxy formation simulations aiming to study early galaxy formation (z ≳ 3) must employ a spatially inhomogeneous UV background to accurately reproduce galaxy properties.more » « less
-
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
