Title: The Life and Times of Star-forming Cores: An Analysis of Dense Gas in the STARFORGE Simulations
Dense gas in molecular clouds is an important signature of ongoing and future star formation. We identify and track dense cores in the Starforge simulations, following the core evolution from birth through dispersal by stellar feedback for typical Milky Way cloud conditions. Only ∼8% of cores host protostars, and most disperse before forming stars. The median starless and protostellar core lifetimes are ∼0.5–0.6 Myr and ∼0.8–1.1 Myr, respectively, where the protostellar phase lasts 0.1 Myr. While core evolution is stochastic, we find that virial ratios and line widths decline in prestellar cores, coincident with turbulent decay. Collapse occurs over ∼0.1 Myr, once the central density exceeds ≳106cm−3. Starless cores, only, follow line-width–size and mass–size relations,σ∝R0.3andM∝R1. The core median mass, radius, and velocity dispersion scale weakly with the cloud magnetic field strength. We cluster the core properties and find that protostellar cores have >80% likelihood of belonging to three particular groups that are characterized by high central densities, compact radii, and lower virial parameters. Overall, core evolution appears to be universally set by the interplay of gravity and magnetized turbulence, while stellar feedback dictates protostellar core properties and sets the protostellar phase lifetime.  more » « less
Award ID(s):
2107942 2421782
PAR ID:
10625498
Author(s) / Creator(s):
; ;
Publisher / Repository:
ApJ
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
982
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
138
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Star formation occurs within dense regions of giant molecular clouds (GMCs); however, exactly how gas collects and evolves to form individual stars and what role dense cores play remains unclear. We use the Lagrangian cell information in the STARFORGE simulation suite to track star-forming gas in three GMCs with varying magnetic field strengths. We find that once a protostar forms, the lifetime of the unaccreted gas correlates with the final stellar mass, where low-mass stars (M* < 0.5M⊙) accrete for 0.5–0.6 Myr from a relatively local reservoir of gas and high-mass stars (M* > 2M⊙) accrete over 3.3–4.7 Myr from a much larger volume. Although the protostellar accretion time increases weakly with magnetic field strength, the accreting gas radii, velocity dispersions, virial parameters, and magnetic energy ratios are largely insensitive to the global cloud properties. At the time of protostar formation, the unaccreted gas exhibits linewidth-size and mass-size relations characteristic of turbulently regulated, isothermal dense cores, followingσv ∝ R1.0−1.1andM ∝ R0.47−0.55, respectively. Low- and intermediate-mass stars undergo relatively continuous accretion, and their accretion histories are well-fit by isothermal sphere, turbulent core, or competitive accretion models, where no one model fits all masses. However, many high-mass stars experience intermittent accretion, and their accretion histories are not well-fit by any of these models. While the distribution of accreting gas is more extended than typically defined dense cores, the physical properties and structure of the star-forming gas resemble those of observed cores and are largely regulated by turbulence and feedback. 
    more » « less
  2. Stars form in dense cores within molecular clouds, and newly formed stars influence their natal environments. How stellar feedback impacts core properties and evolution has been the subject of extensive investigation. We performed a hierarchical clustering (dendrogram) analysis of a STARFORGE (STAR FORmation in Gaseous Environments) simulation, modelling a giant molecular cloud to identify gas overdensities (cores) and study changes in their radius, mass, velocity dispersion, and virial parameter with respect to stellar feedback. We binned these cores on the basis of the fraction of gas affected by protostellar outflows, stellar winds, and supernovae and analysed the property distributions for each feedback bin. We find that cores that experience more feedback influence are smaller. Feedback notably enhances the velocity dispersion and virial parameter of the cores, more so than it reduces their radius. This is also evident in the linewidth–size relation, according to which cores in higher-feedback bins exhibit higher velocities than their similarly sized pristine counterparts. We conclude that stellar feedback mechanisms, which impart momentum to the molecular cloud, simultaneously compress and disperse the dense molecular gas. 
    more » « less
  3. ABSTRACT We study the formation, evolution, and collapse of dense cores by tracking structures in a magnetohydrodynamic simulation of a star-forming cloud. We identify cores using the dendrogram algorithm and utilize machine learning techniques, including Neural Gas prototype learning and Fuzzy c-means clustering to analyse the density and velocity dispersion profiles of cores together with six bulk properties. We produce a 2-d visualization using a Uniform Manifold Approximation and Projection (UMAP), which facilitates the connection between physical properties and three partially-overlapping phases: i) unbound turbulent structures (Phase I), ii) coherent cores that have low turbulence (Phase II), and iii) bound cores, many of which become protostellar (Phase III). Within Phase II, we identify a population of long-lived coherent cores that reach a quasi-equilibrium state. Most prestellar cores form in Phase II and become protostellar after evolving into Phase III. Due to the turbulent cloud environment, the initial core properties do not uniquely predict the eventual evolution, i.e. core evolution is stochastic, and cores follow no one evolutionary path. The phase lifetimes are 1.0 ± 0.1 × 105 yr, 1.3 ± 0.2 × 105 yr, and 1.8 ± 0.3 × 105 yr for Phase I, II, and III, respectively. We compare our results to NH3 observations of dense cores. Known coherent cores predominantly map into Phase II, while most turbulent pressure-confined cores map to Phase I or III. We predict that a significant fraction of observed starless cores have unresolved coherent regions and that ≳20 per cent of observed starless cores will not form stars. Measurements of core radial profiles in addition to the usual bulk properties will enable more accurate predictions of core evolution. 
    more » « less
  4. Abstract We present a population study of Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 6 observations of the 100 most gravitationally unstable dense cores in Aquila using a simple mass versus size analysis. We identify 66 continuum sources from ALMA 12 m observations at 106 GHz and through comparisons with known protostellar catalogs; two of these detected dense cores appear to be completely starless, without any accompanying/nearby protostar detections. Additionally, we find nine other starless ALMA 12 m detections within protostellar cores that have fragmented into a mixture of starless and protostellar substructures. We test the turbulent core collapse model by conducting synthetic observations of turbulent magnetohydrodynamical simulations of collapsing starless cores in order to predict how many starless cores should be detected given their central density and density profile. The simulations predict at least one (1.19) detection, consistent with our two detections of ALMA 12 m emission within completely starless cores. We also use a combination of ALMA Compact Array Cycle 4 observations and the Herschel Gould Belt Survey data to analyze how mass is distributed on three distinct spatial scales, in order to understand how turbulence shapes the evolution of substructure development as dense cores collapse to form new star systems. We find an increase in multiplicity at the smallest scales when the parent larger-scale structure also has a higher degree of fragmentation. 
    more » « less
  5. Abstract The mass distribution of dense cores is a potential key to understanding the process of star formation. Applying dendrogram analysis to the CARMA-NRO Orion C 18 O ( J = 1–0) data, we identify 2342 dense cores, about 22% of which have virial ratios smaller than 2 and can be classified as gravitationally bound cores. The derived core mass function (CMF) for bound starless cores that are not associate with protostars has a slope similar to Salpeter’s initial mass function (IMF) for the mass range above 1 M ⊙ , with a peak at ∼0.1 M ⊙ . We divide the cloud into four parts based on decl., OMC-1/2/3, OMC-4/5, L1641N/V380 Ori, and L1641C, and derive the CMFs in these regions. We find that starless cores with masses greater than 10 M ⊙ exist only in OMC-1/2/3, whereas the CMFs in OMC-4/5, L1641N, and L1641C are truncated at around 5–10 M ⊙ . From the number ratio of bound starless cores and Class II objects in each subregion, the lifetime of bound starless cores is estimated to be 5–30 freefall times, consistent with previous studies for other regions. In addition, we discuss core growth by mass accretion from the surrounding cloud material to explain the coincidence of peak masses between IMFs and CMFs. The mass accretion rate required for doubling the core mass within a core lifetime is larger than that of Bondi–Hoyle accretion by a factor of order 2. This implies that more dynamical accretion processes are required to grow cores. 
    more » « less