Stellar clusters are critical constituents within galaxies: They are the result of highest-density star formation, and through their spatially and temporally correlated feedback, they regulate their host galaxy evolution. We present a novel numerical method to model star clusters as individual units of star formation using sink particles. In our method, star clusters grow via gas accretion and via merging with less massive clusters. We describe the implementation in the radiation hydrodynamics code
This content will become publicly available on October 1, 2025
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 T
- PAR ID:
- 10547502
- Publisher / Repository:
- EDP Sciences
- Date Published:
- Journal Name:
- Astronomy & Astrophysics
- Volume:
- 690
- ISSN:
- 0004-6361
- Page Range / eLocation ID:
- A94
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract GIZMO and run a large grid of marginally bound, turbulent clouds of 107M ⊙to explore the effect of modeling ingredients on the evolution of the clouds and the star clusters. We find both gas accretion and mergers to be critical processes to form star clusters of masses up to ∼105–106M ⊙, while ionizing radiation is the main feedback mechanism regulating the growth of star clusters. The majority of our star clusters assemble their mass in 0.3–2.6 Myr, and the most massive ones take ∼10 Myr. By removing high-density gas by accretion, our sink-based cluster formation prescription allows the newly formed star clusters to inject their stellar feedback in less dense environments. This makes feedback more efficient at ionizing and disrupting the cloud than if we were to use a standard star formation approach, indicating that our numerical method is the missing critical step to model the interplay between star clusters and their host galaxies. -
more » « less
Abstract We study the evolution of populations of binary stars within massive cluster-forming regions. We simulate the formation of young massive star clusters within giant molecular clouds with masses ranging from 2 × 104to 3.2 × 105
M ⊙. We use Torch , which couples stellar dynamics, magnetohydrodynamics, star and binary formation, stellar evolution, and stellar feedback through the Amuse framework. We find that the binary fraction decreases during cluster formation at all molecular cloud masses. The binaries’ orbital properties also change, with stronger and quicker changes in denser, more massive clouds. Most of the changes we see can be attributed to the disruption of binaries wider than 100 au, although the close binary fraction also decreases in the densest cluster-forming region. The binary fraction for O stars remains above 90%, but exchanges and dynamical hardening are ubiquitous, indicating that O stars undergo frequent few-body interactions early during the cluster formation process. Changes to the populations of binaries are a by-product of hierarchical cluster assembly: most changes to the binary population take place when the star formation rate is high, and there are frequent mergers between subclusters in the cluster-forming region. A universal primordial binary distribution based on observed inner companions in the Galactic field is consistent with the binary populations of young clusters with resolved stellar populations, and the scatter between clusters of similar masses could be explained by differences in their formation history. -
Nitrogen-enriched, Highly Pressurized Nebular Clouds Surrounding a Super Star Cluster at Cosmic Noonmore » « less
Abstract Strong lensing offers a precious opportunity for studying the formation and early evolution of super star clusters that are rare in our cosmic backyard. The Sunburst Arc, a lensed Cosmic Noon galaxy, hosts a young super star cluster with escaping Lyman continuum radiation. Analyzing archival Hubble Space Telescope images and emission line data from Very Large Telescope/MUSE and X-shooter, we construct a physical model for the cluster and its surrounding photoionized nebula. We confirm that the cluster is ≲4 Myr old, is extremely massive
M ⋆∼ 107M ⊙, and yet has a central component as compact as several parsecs, and we find a gas-phase metallicityZ = (0.22 ± 0.03)Z ⊙. The cluster is surrounded by ≳105M ⊙of dense clouds that have been pressurized toP ∼ 109K cm−3by perhaps stellar radiation at within 10 pc. These should have large neutral columnsN HI> 1022.8cm−2to survive rapid ejection by radiation pressure. The clouds are likely dusty as they show gas-phase depletion of silicon, and may be conducive to secondary star formation ifN HI> 1024cm−2or if they sink farther toward the cluster center. Detecting strong [Niii ]λ λ 1750,1752, we infer heavy nitrogen enrichmentM ⊙of nitrogen in the high-pressure clouds from massive stars heavier than 60M ⊙up to 4 Myr. We suggest a physical origin of the high-pressure clouds from partial or complete condensation of slow massive star ejecta, which may have an important implication for the puzzle of multiple stellar populations in globular clusters. -
more » « less
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 suite
Torch , 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 theTorch software suite through theAMUSE framework, 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
Abstract We introduce the Phoenix Simulations, a suite of highly resolved cosmological simulations featuring hydrodynamics, primordial gas chemistry, primordial and enriched star formation and feedback, UV radiative transfer, and saved outputs with Δ
t = 200 kyr. We observe 73,523 individual primordial stars within 3313 distinct regions forming 2110 second-generation enriched star clusters byz ≥ 12 within a combined 177.25 Mpc3volume across three simulations. The regions that lead to enriched star formation can contain ≳150 primordial stars, with 80% of regions having experienced combinations of primordial Type II, hypernovae, and/or pair-instability supernovae. Primordial supernovae enriched 0.8% of the volume, with 2% of enriched gas enriched by later-generation stars. We determine the extent of a primordial stellar region by its metal-rich or ionized hydrogen surrounding cloud; the metal-rich and ionized regions have time-dependent average radiir ≲ 3kpc. 7 and 17% of regions haver > 7 kpc for metal-rich and ionized radii, respectively. We find that the metallicity distribution function of second-generation stars overlaps that of subsequent Population II star formation, spanning metal-deficient (∼7.94 × 10−8Z ⊙) to supersolar (∼3.71Z ⊙), and that 30.5% of second-generation stars haveZ > 10−2Z ⊙. We find that the metallicity of second-generation stars depends on progenitor configuration, with metals from pair-instability supernovae contributing to the most metal-rich clusters; these clusters form promptly after the supernova event. Finally, we create an interpretable regression model to predict the radius of the metal-rich influence of Population III star systems within the first 7–18 Myr after the first Population III stars form in the region.
