skip to main content


Title: Modeling of the Effects of Stellar Feedback during Star Cluster Formation Using a Hybrid Gas and N-Body Method
Abstract

Understanding the formation of stellar clusters requires following the interplay between gas and newly formed stars accurately. We therefore couple the magnetohydrodynamics codeFLASHto theN-body codeph4and the stellar evolution codeSeBausing the Astrophysical Multipurpose Software Environment (AMUSE) to model stellar dynamics, evolution, and collisionalN-body dynamics and the formation of binary and higher-order multiple systems, while implementing stellar feedback in the form of radiation, stellar winds, and supernovae inFLASH. We here describe the algorithms used for each of these processes. We denote this integrated package Torch. We then use this novel numerical method to simulate the formation and early evolution of several examples of open clusters of ∼1000 stars formed from clouds with a mass range of 103Mto 105M. Analyzing the effects of stellar feedback on the gas and stars of the natal clusters, we find that in these examples, the stellar clusters are resilient to disruption, even in the presence of intense feedback. This can even slightly increase the amount of dense, Jeans unstable gas by sweeping up shells; thus, a stellar wind strong enough to trap its own H iiregion shows modest triggering of star formation. Our clusters are born moderately mass segregated, an effect enhanced by feedback, and retained after the ejection of their natal gas, in agreement with observations.

 
more » « less
Award ID(s):
1815461 1109395 1814772
NSF-PAR ID:
10485931
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
904
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 192
Size(s):
["Article No. 192"]
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 104Mgiant 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
  2. Abstract

    We describe the public release of the Cluster Monte Carlo (CMC) code, a parallel, star-by-starN-body code for modeling dense star clusters.CMCtreats collisional stellar dynamics using Hénon’s method, where the cumulative effect of many two-body encounters is statistically reproduced as a single effective encounter between nearest-neighbor particles on a relaxation timescale. The star-by-star approach allows for the inclusion of additional physics, including strong gravitational three- and four-body encounters, two-body tidal and gravitational-wave captures, mass loss in arbitrary galactic tidal fields, and stellar evolution for both single and binary stars. The public release ofCMCis pinned directly to theCOSMICpopulation synthesis code, allowing dynamical star cluster simulations and population synthesis studies to be performed using identical assumptions about the stellar physics and initial conditions. As a demonstration, we present two examples of star cluster modeling: first, we perform the largest (N= 108) star-by-starN-body simulation of a Plummer sphere evolving to core collapse, reproducing the expected self-similar density profile over more than 15 orders of magnitude; second, we generate realistic models for typical globular clusters, and we show that their dynamical evolution can produce significant numbers of black hole mergers with masses greater than those produced from isolated binary evolution (such as GW190521, a recently reported merger with component masses in the pulsational pair-instability mass gap).

     
    more » « less
  3. Abstract

    The globular cluster 47 Tucanae (47 Tuc) is one of the most massive star clusters in the Milky Way and is exceptionally rich in exotic stellar populations. For several decades it has been a favorite target of observers, and yet it is computationally very challenging to model because of its large number of stars (N≳ 106) and high density. Here we present detailed and self-consistent 47 Tuc models computed with theCluster Monte Carlocode (CMC). The models include all relevant dynamical interactions coupled to stellar and binary evolution, and reproduce various observations, including the surface brightness and velocity dispersion profiles, pulsar accelerations, and numbers of compact objects. We show that the present properties of 47 Tuc are best reproduced by adopting an initial stellar mass function that is both bottom-heavy and top-light relative to standard assumptions (as in, e.g., Kroupa 2001), and an initial Elson profile (Elson et al. 1987) that is overfilling the cluster’s tidal radius. We include new prescriptions inCMCfor the formation of binaries through giant star collisions and tidal captures, and we show that these mechanisms play a crucial role in the formation of neutron star binaries and millisecond pulsars in 47 Tuc; our best-fit model contains ∼50 millisecond pulsars, 70% of which are formed through giant collisions and tidal captures. Our models also suggest that 47 Tuc presently contains up to ∼200 stellar-mass black holes, ∼5 binary black holes, ∼15 low-mass X-ray binaries, and ∼300 cataclysmic variables.

     
    more » « less
  4. Abstract

    Gravitational lenses can magnify distant galaxies, allowing us to discover and characterize the stellar populations of intrinsically faint, quiescent galaxies that are otherwise extremely difficult to directly observe at high redshift from ground-based telescopes. Here, we present the spectral analysis of two lensed, quiescent galaxies atz≳ 1 discovered by theASTRO 3D Galaxy Evolution with Lensessurvey:AGEL1323 (M*∼ 1011.1M,z= 1.016,μ∼ 14.6) andAGEL0014 (M*∼ 1011.5M,z= 1.374,μ∼ 4.3). We measured the age, [Fe/H], and [Mg/Fe] of the two lensed galaxies using deep, rest-frame-optical spectra (S/N ≳40 Å−1) obtained on the Keck I telescope. The ages ofAGEL1323 andAGEL0014 are5.60.8+0.8Gyr and3.10.3+0.8Gyr, respectively, indicating that most of the stars in the galaxies were formed less than 2 Gyr after the Big Bang. Compared to nearby quiescent galaxies of similar masses, the lensed galaxies have lower [Fe/H] and [Mg/H]. Surprisingly, the two galaxies have comparable [Mg/Fe] to similar-mass galaxies at lower redshifts, despite their old ages. Using a simple analytic chemical evolution model connecting the instantaneously recycled element Mg with the mass-loading factors of outflows averaged over the entire star formation history, we found that the lensed galaxies may have experienced enhanced outflows during their star formation compared to lower-redshift galaxies, which may explain why they quenched early.

     
    more » « less
  5. 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 massiveM∼ 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 ≳105Mof dense clouds that have been pressurized toP∼ 109K cm−3by perhaps stellar radiation at within 10 pc. These should have large neutral columnsNHI> 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 ifNHI> 1024cm−2or if they sink farther toward the cluster center. Detecting strong [Niii]λλ1750,1752, we infer heavy nitrogen enrichmentlog(N/O)=0.210.11+0.10. This requires efficiently retaining ≳500Mof nitrogen in the high-pressure clouds from massive stars heavier than 60Mup 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