The current generation of galaxy simulations can resolve individual giant molecular clouds, the progenitors of dense star clusters. But the evolutionary fate of these young massive clusters, and whether they can become the old globular clusters (GCs) observed in many galaxies, is determined by a complex interplay of internal dynamical processes and external galactic effects. We present the first star-by-star N-body models of massive (N ∼ 105–107) star clusters formed in a FIRE-2 MHD simulation of a Milky Way-mass galaxy, with the relevant initial conditions and tidal forces extracted from the cosmological simulation. We select 895 (∼30 per cent) of the YMCs with >6 × 104 M⊙ from Grudić et al. 2022 and integrate them to z = 0 using the cluster Monte Carlo code, CMC. This procedure predicts a MW-like system with 148 GCs, predominantly formed during the early, bursty mode of star formation. Our GCs are younger, less massive, and more core-collapsed than clusters in the Milky Way or M31. This results from the assembly history and age-metallicity relationship of the host galaxy: Younger clusters are preferentially born in stronger tidal fields and initially retain fewer stellar-mass black holes, causing them to lose mass faster and reach core collapse sooner than older GCs. Our results suggest that the masses and core/half-light radii of GCs are shaped not only by internal dynamical processes, but also by the specific evolutionary history of their host galaxies. These results emphasize that N-body studies with realistic stellar physics are crucial to understanding the evolution and present-day properties of GC systems.
Milky Way globular clusters (GCs) display chemical enrichment in a phenomenon called multiple stellar populations (MSPs). While the enrichment mechanism is not fully understood, there is a correlation between a cluster’s mass and the fraction of enriched stars found therein. However, present-day GC masses are often smaller than their masses at the time of formation due to dynamical mass-loss. In this work, we explore the relationship between mass and MSPs using the stellar stream 300S. We present the chemical abundances of eight red giant branch member stars in 300S with high-resolution spectroscopy from Magellan/MIKE. We identify one enriched star characteristic of MSPs and no detectable metallicity dispersion, confirming that the progenitor of 300S was a GC. The fraction of enriched stars (12.5 per cent) observed in our 300S stars is less than the 50 per cent of stars found enriched in Milky Way GCs of comparable present-day mass (∼104.5 $\mathrm{\, {\rm M}_{\odot }}$). We calculate the mass of 300S’s progenitor and compare it to the initial masses of intact GCs, finding that 300S aligns well with the trend between the system mass at formation and enrichment. 300S’s progenitor may straddle the critical mass threshold for the formation of MSPs and can therefore serve as a benchmark for the stellar enrichment process. Additionally, we identify a CH star, with high abundances of s-process elements, probably accreted from a binary companion. The rarity of such binaries in intact GCs may imply stellar streams permit the survival of binaries that would otherwise be disrupted.
more » « less- NSF-PAR ID:
- 10496216
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 529
- Issue:
- 3
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 2413-2427
- Size(s):
- ["p. 2413-2427"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
Abstract We present stellar parameters and abundances of 13 elements for 18 very metal-poor (VMP; [Fe/H] < –2.0) stars, selected as extremely metal-poor (EMP; [Fe/H] < –3.0) candidates from the Sloan Digital Sky Survey and Large sky Area Multi-Object Fiber Spectroscopic Telescope survey. High-resolution spectroscopic observations were performed using GEMINI-N/GRACES. We find 10 EMP stars among our candidates, and we newly identify three carbon-enhanced metal-poor stars with [Ba/Fe] < 0. Although chemical abundances of our VMP/EMP stars generally follow the overall trend of other Galactic halo stars, there are a few exceptions. One Na-rich star ([Na/Fe] = +1.14) with low [Mg/Fe] suggests a possible chemical connection with second-generation stars in a globular cluster. The progenitor of an extremely Na-poor star ([Na/Fe] = –1.02) with high K- and Ni-abundance ratios may have undergone a distinct nucleosynthesis episode, associated with core-collapse supernovae (SNe) having a high explosion energy. We have also found a Mg-rich star ([Mg/Fe] = +0.73) with slightly enhanced Na and extremely low [Ba/Fe], indicating that its origin is not associated with neutron-capture events. On the other hand, the origin of the lowest Mg abundance ([Mg/Fe] = –0.61) star could be explained by accretion from a dwarf galaxy, or formation in a gas cloud largely polluted by SNe Ia. We have also explored the progenitor masses of our EMP stars by comparing their chemical-abundance patterns with those predicted by Population III SNe models, and find a mass range of 10–26 M ⊙ , suggesting that such stars were primarily responsible for the chemical enrichment of the early Milky Way.more » « less
-
more » « less
Abstract In the era of large-scale spectroscopic surveys in the Local Group, we can explore using chemical abundances of halo stars to study the star formation and chemical enrichment histories of the dwarf galaxy progenitors of the Milky Way (MW) and M31 stellar halos. In this paper, we investigate using the chemical abundance ratio distributions (CARDs) of seven stellar halos from the Latte suite of FIRE-2 simulations. We attempt to infer galaxies’ assembly histories by modeling the CARDs of the stellar halos of the Latte galaxies as a linear combination of
template CARDs from disrupted dwarfs, with different stellar massesM ⋆and quenching timest 100. We present a method for constructing these templates using present-day dwarf galaxies. For four of the seven Latte halos studied in this work, we recover the mass spectrum of accreted dwarfs to a precision of <10%. For the fraction of mass accreted as a function oft 100, we find the residuals of 20%–30% for five of the seven simulations. We discuss the failure modes of this method, which arise from the diversity of star formation and chemical enrichment histories that dwarf galaxies can take. These failure cases can be robustly identified by the high model residuals. Although the CARDs modeling method does not successfully infer the assembly histories in these cases, the CARDs of these disrupted dwarfs contain signatures of their unusual formation histories. Our results are promising for using CARDs to learn more about the histories of the progenitors of the MW and M31 stellar halos. -
ABSTRACT The contribution of dissolved globular clusters (GCs) to the stellar content of the Galactic halo is a key constraint on models for GC formation and destruction, and the mass assembly history of the Milky Way. Earlier results from APOGEE pointed to a large contribution of destroyed GCs to the stellar content of the inner halo, by as much as 25 ${{\ \rm per\ cent}}$, which is an order of magnitude larger than previous estimates for more distant regions of the halo. We set out to measure the ratio between nitrogen-rich (N-rich) and normal halo field stars, as a function of distance, by performing density modelling of halo field populations in APOGEE DR16. Our results show that at 1.5 kpc from the Galactic Centre, N-rich stars contribute a much higher 16.8$^{+10.0}_{-7.0}\, {{\ \rm per\ cent}}$ fraction to the total stellar halo mass budget than the 2.7$^{+1.0}_{-0.8}\, {{\ \rm per\ cent}}$ ratio contributed at 10 kpc. Under the assumption that N-rich stars are former GC members that now reside in the stellar halo field, and assuming the ratio between first and second population GC stars being 1:2, we estimate a total contribution from disrupted GC stars of the order of 27.5$^{+15.4}_{-11.5}\, {{\ \rm per\ cent}}$ at r = 1.5 kpc and 4.2$^{+1.5}_{-1.3}\, {{\ \rm per\ cent}}$ at r = 10 kpc. Furthermore, since our methodology requires fitting a density model to the stellar halo, we integrate such density within a spherical shell from 1.5 to 15 kpc in radius, and find a total stellar mass arising from dissolved and/or evaporated GCs of MGC,total = 9.6$^{+4.0}_{-2.6}\, \times$ 107 M⊙.more » « less
-
null (Ed.)ABSTRACT By means of 3D hydrodynamic simulations, we study how Type Ia supernovae (SNe) explosions affect the star formation history and the chemical properties of second-generation (SG) stars in globular clusters (GC). SG stars are assumed to form once first generation asymptotic giant branch (AGB) stars start releasing their ejecta; during this phase, external gas is accreted by the system and SNe Ia begin exploding, carving hot and tenuous bubbles. Given the large uncertainty on SNe Ia explosion times, we test two different values for the ‘delay time’. We run two different models for the external gas density: in the low-density scenario with short delay time, the explosions start at the beginning of the SG star formation, halting it in its earliest phases. The external gas hardly penetrates the system, therefore most SG stars present extreme helium abundances (Y > 0.33). The low-density model with delayed SN explosions has a more extended SG star formation epoch and includes SG stars with modest helium enrichment. On the contrary, the high-density model is weakly affected by SN explosions, with a final SG mass similar to the one obtained without SNe Ia. Most of the stars form from a mix of AGB ejecta and pristine gas and have a modest helium enrichment. We show that gas from SNe Ia may produce an iron spread of ∼0.14 dex, consistent with the spread found in about $20{{\ \rm per\ cent}}$ of Galactic GCs, suggesting that SNe Ia might have played a key role in the formation of this sub-sample of GCs.more » « less
