Globular clusters (GCs) are particularly efficient at forming millisecond pulsars. Among these pulsars, about half lack a companion star, a significantly higher fraction than in the Galactic field. This fraction increases further in some of the densest GCs, especially those that have undergone core collapse, suggesting that dynamical interaction processes play a key role. For the first time, we create
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.
more » « less- NSF-PAR ID:
- 10400092
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 521
- Issue:
- 1
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 124-147
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract N -body models that reproduce the ratio of single-to-binary pulsars in Milky Way–like GCs. We focus especially on NGC 6752, a typical core-collapsed cluster with many observed millisecond pulsars. Previous studies suggested that an increased rate of neutron star binary disruption in the densest clusters could explain the overabundance of single pulsars in these systems. Here, we demonstrate that binary disruption is ineffective and instead we propose that two additional dynamical processes play dominant roles: (1) tidal disruption of main-sequence stars by neutron stars and (2) gravitational collapse of heavy white dwarf binary merger remnants. Neutron stars formed through these processes may also be associated with fast radio bursts similar to those observed recently in an extragalactic GC. -
more » « less
ABSTRACT Galaxy mergers are known to host abundant young massive cluster (YMC) populations, whose formation mechanism is still not well-understood. Here, we present a high-resolution galaxy merger simulation with explicit star formation and stellar feedback prescriptions to investigate how mergers affect the properties of the interstellar medium and YMCs. Compared with a controlled simulation of an isolated galaxy, the mass fraction of dense and high-pressure gas is much higher in mergers. Consequently, the mass function of both molecular clouds and YMCs becomes shallower and extends to higher masses. Moreover, cluster formation efficiency is significantly enhanced and correlates positively with the star formation rate surface density and gas pressure. We track the orbits of YMCs and investigate the time evolution of tidal fields during the course of the merger. At an early stage of the merger, the tidal field strength correlates positively with YMC mass, λtid ∝ M0.71, which systematically affects the shape of the mass function and age distribution of the YMCs. At later times, most YMCs closely follow the orbits of their host galaxies, gradually sinking into the centre of the merger remnant due to dynamical friction, and are quickly dissolved via efficient tidal disruption. Interestingly, YMCs formed during the first passage, mostly in tidal tails and bridges, are distributed over a wide range of galactocentric radii, greatly increasing their survivability because of the much weaker tidal field in the outskirts of the merger system. These YMCs are promising candidates for globular clusters that survive to the present day.
-
more » « less
Abstract Globular clusters (GCs) provide valuable insight into the properties of their host galaxies’ dark matter halos. Using
N -body simulations incorporating semianalytic dynamical friction and GC−GC merger prescriptions, we study the evolution of GC radial distributions and mass functions in cuspy and cored dark matter halos. Modeling the dynamics of the GC-rich system in the dwarf galaxy UGC 7369, we find that friction-induced inspiral and subsequent mergers of massive GCs can naturally and robustly explain the mass segregation of the GCs and the existence of a nuclear star cluster (NSC). However, the multiple mergers required to form the NSC only take place when the dark matter halo is cuspy. In a cored halo, stalling of the dynamical friction within the core halts the inspiral of the GCs, and so the GC merger rate falls significantly, precluding the formation of an NSC. We therefore argue that the presence of an NSC requires a cusp in UGC 7369. More generally, we propose that the presence of an NSC and the corresponding alteration of the GC mass function due to mergers may be used as an indicator of a cuspy halo for galaxies in which we expect NSC formation to be merger dominated. These observables represent a simple, powerful complement to other inner halo density profile constraint techniques and should allow for straightforward extension to larger samples. -
more » « less
Abstract The recent discovery of two detached black hole–star (BH–star) binaries from Gaia’s third data release has sparked interest in understanding the formation mechanisms of these systems. We investigate the formation of these systems by dynamical processes in young star clusters (SCs) and via isolated binary (IB) evolution, using a combination of direct
N -body and population synthesis simulations. We find that dynamical formation in SCs is nearly 50 times more efficient per unit of star formation at producing BH–star binaries than IB evolution. We expand this analysis to the full Milky Way (MW) using a FIRE-2 hydrodynamical simulation of an MW-mass galaxy. Even assuming that only 10% of star formation goes into SCs, we find that approximately four out of every five BH–star systems are formed dynamically, and that the MW contains a total of ∼2 × 105BH–star systems. Many of these dynamically formed systems have longer orbital periods, greater eccentricities, and greater black hole masses than their isolated counterparts. For binaries older than 100 Myr, we show that any detectable system withe ≳ 0.5 orM BH≳ 10M ⊙canonly be formed through dynamical processes. Our MW model predicts between 64 and 215 such detections from the complete DR4 Gaia catalog, with the majority of systems being dynamically formed in massive and metal-rich SCs. Finally, we compare our populations to the recently discovered Gaia BH1 and Gaia BH2, and conclude that the dynamical scenario is the most favorable formation pathway for both systems. -
more » « less
ABSTRACT The velocity dispersion of globular clusters (GCs) around ultra-diffuse galaxies (UDGs) in the Virgo cluster spans a wide range, including cases where GC kinematics suggest haloes as massive as (or even more massive than) that of the Milky Way around these faint dwarfs. We analyse the catalogues of GCs derived in post-processing from the TNG50 cosmological simulation to study the GC system kinematics and abundance of simulated UDGs in galaxy groups and clusters. UDGs in this simulation reside exclusively in dwarf-mass haloes with M200 ≲ 1011.2 M⊙. When considering only GCs gravitationally bound to simulated UDGs, we find GCs properties that overlap well with several observational measurements for UDGs. In particular, no bias towards overly massive haloes is inferred from the study of bound GCs, confirming that GCs are good tracers of UDG halo mass. However, we find that contamination by intracluster GCs may, in some cases, substantially increase velocity dispersion estimates when performing projected mock observations of our sample. We caution that targets with less than 10 GC tracers are particularly prone to severe uncertainties. Measuring the stellar kinematics of the host galaxy should help confirm the unusually massive haloes suggested by GC kinematics around some UDGs.
