skip to main content


Title: Turbulence, coherence, and collapse: Three phases for core evolution
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
Award ID(s):
2142300 2008101 2107942 1908419
NSF-PAR ID:
10403473
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
517
Issue:
1
ISSN:
0035-8711
Page Range / eLocation ID:
885 to 909
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Context. We started a multi-scale analysis of star formation in G202.3+2.5, an intertwined filamentary sub-region of the Monoceros OB1 molecular complex, in order to provide observational constraints on current theories and models that attempt to explain star formation globally. In the first paper (Paper I), we examined the distributions of dense cores and protostars and found enhanced star formation activity in the junction region of the filaments. Aims. In this second paper, we aim to unveil the connections between the core and filament evolutions, and between the filament dynamics and the global evolution of the cloud. Methods. We characterise the gas dynamics and energy balance in different parts of G202.3+2.5 using infrared observations from the Herschel and WISE telescopes and molecular tracers observed with the IRAM 30-m and TRAO 14-m telescopes. The velocity field of the cloud is examined and velocity-coherent structures are identified, characterised, and put in perspective with the cloud environment. Results. Two main velocity components are revealed, well separated in radial velocities in the north and merged around the location of intense N 2 H + emission in the centre of G202.3+2.5 where Paper I found the peak of star formation activity. We show that the relative position of the two components along the sightline, and the velocity gradient of the N 2 H + emission imply that the components have been undergoing collision for ~10 5 yr, although it remains unclear whether the gas moves mainly along or across the filament axes. The dense gas where N 2 H + is detected is interpreted as the compressed region between the two filaments, which corresponds to a high mass inflow rate of ~1 × 10 −3 M ⊙ yr −1 and possibly leads to a significant increase in its star formation efficiency. We identify a protostellar source in the junction region that possibly powers two crossed intermittent outflows. We show that the H  II region around the nearby cluster NCG 2264 is still expanding and its role in the collision is examined. However, we cannot rule out the idea that the collision arises mostly from the global collapse of the cloud. Conclusions. The (sub-)filament-scale observables examined in this paper reveal a collision between G202.3+2.5 sub-structures and its probable role in feeding the cores in the junction region. To shed more light on this link between core and filament evolutions, one must characterise the cloud morphology, its fragmentation, and magnetic field, all at high resolution. We consider the role of the environment in this paper, but a larger-scale study of this region is now necessary to investigate the scenario of a global cloud collapse. 
    more » « less
  2. ABSTRACT

    Young massive clusters (YMCs) are compact (≲1 pc), high-mass (>104 M⊙) stellar systems of significant scientific interest. Due to their rarity and rapid formation, we have very few examples of YMC progenitor gas clouds before star formation has begun. As a result, the initial conditions required for YMC formation are uncertain. We present high resolution (0.13 arcsec, ∼1000 au) ALMA observations and Mopra single-dish data, showing that Galactic Centre dust ridge ‘Cloud d’ (G0.412 + 0.052, mass = 7.6 × 104 M⊙, radius = 3.2 pc) has the potential to become an Arches-like YMC (104 M⊙, r ∼ 1 pc), but is not yet forming stars. This would mean it is the youngest known pre-star-forming massive cluster and therefore could be an ideal laboratory for studying the initial conditions of YMC formation. We find 96 sources in the dust continuum, with masses ≲3 M⊙ and radii of ∼103 au. The source masses and separations are more consistent with thermal rather than turbulent fragmentation. It is not possible to unambiguously determine the dynamical state of most of the sources, as the uncertainty on virial parameter estimates is large. We find evidence for large-scale (∼1 pc) converging gas flows, which could cause the cloud to grow rapidly, gaining 104 M⊙ within 105 yr. The highest density gas is found at the convergent point of the large-scale flows. We expect this cloud to form many high-mass stars, but find no high-mass starless cores. If the sources represent the initial conditions for star formation, the resulting initial mass function will be bottom heavy.

     
    more » « less
  3. Abstract Magnetic fields have an important role in the evolution of interstellar medium and star formation 1,2 . As the only direct probe of interstellar field strength, credible Zeeman measurements remain sparse owing to the lack of suitable Zeeman probes, particularly for cold, molecular gas 3 . Here we report the detection of a magnetic field of +3.8 ± 0.3 microgauss through the H  I narrow self-absorption (HINSA) 4,5 towards L1544 6,7 —a well-studied prototypical prestellar core in an early transition between starless and protostellar phases 8–10 characterized by a high central number density 11 and a low central temperature 12 . A combined analysis of the Zeeman measurements of quasar H  I absorption, H  I emission, OH emission and HINSA reveals a coherent magnetic field from the atomic cold neutral medium (CNM) to the molecular envelope. The molecular envelope traced by the HINSA is found to be magnetically supercritical, with a field strength comparable to that of the surrounding diffuse, magnetically subcritical CNM despite a large increase in density. The reduction of the magnetic flux relative to the mass, which is necessary for star formation, thus seems to have already happened during the transition from the diffuse CNM to the molecular gas traced by the HINSA. This is earlier than envisioned in the classical picture where magnetically supercritical cores capable of collapsing into stars form out of magnetically subcritical envelopes 13,14 . 
    more » « less
  4. null (Ed.)
    ABSTRACT We discuss a theoretical model for the early evolution of massive star clusters and confront it with the ALMA, radio, and infrared observations of the young stellar cluster highly obscured by the molecular cloud D1 in the nearby dwarf spheroidal galaxy NGC 5253. We show that a large turbulent pressure in the central zones of D1 cluster may cause individual wind-blown bubbles to reach pressure confinement before encountering their neighbours. In this case, stellar winds energy is added to the hot shocked wind pockets of gas around individual massive stars that leads them to meet and produce a cluster wind in time-scales less than 105 yr. In order to inhibit the possibility of cloud dispersal, or the early negative star formation feedback, one should account for mass loading that may come, for example, from pre-main-sequence (PMS) low-mass stars through photoevaporation of their protostellar discs. Mass loading at a rate in excess of 8 × 10−9 M⊙ yr−1 per each PMS star is required to extend the hidden star cluster phase in this particular cluster. In this regime, the parental cloud remains relatively unperturbed, while pockets of molecular, photoionized and hot gas coexist within the star-forming region. Nevertheless, the most likely scenario for cloud D1 and its embedded cluster is that the hot shocked winds around individual massive stars should merge at an age of a few million of years when the PMS star protostellar discs vanish and mass loading ceases that allows a cluster to form a global wind. 
    more » « less
  5. null (Ed.)
    Context. Stars form in cold dense cores showing subsonic velocity dispersions. The parental molecular clouds display higher temperatures and supersonic velocity dispersions. The transition from core to cloud has been observed in velocity dispersion, but temperature and abundance variations are unknown. Aims. We aim to measure the temperature and velocity dispersion across cores and ambient cloud in a single tracer to study the transition between the two regions. Methods. We use NH 3 (1,1) and (2,2) maps in L1688 from the Green Bank Ammonia Survey, smoothed to 1′, and determine the physical properties by fitting the spectra. We identify the coherent cores and study the changes in temperature and velocity dispersion from the cores to the surrounding cloud. Results. We obtain a kinetic temperature map extending beyond dense cores and tracing the cloud, improving from previous maps tracing mostly the cores. The cloud is 4–6 K warmer than the cores, and shows a larger velocity dispersion (Δ σ v = 0.15–0.25 km s −1 ). Comparing to Herschel -based dust temperatures, we find that cores show kinetic temperatures that are ≈1.8 K lower than the dust temperature, while the gas temperature is higher than the dust temperature in the cloud. We find an average p-NH 3 fractional abundance (with respect to H 2 ) of (4.2 ± 0.2) × 10 −9 towards the coherent cores, and (1.4 ± 0.1) × 10 −9 outside the core boundaries. Using stacked spectra, we detect two components, one narrow and one broad, towards cores and their neighbourhoods. We find the turbulence in the narrow component to be correlated with the size of the structure (Pearson- r = 0.54). With these unresolved regional measurements, we obtain a turbulence–size relation of σ v,NT ∝ r 0.5 , which is similar to previous findings using multiple tracers. Conclusions. We discover that the subsonic component extends up to 0.15 pc beyond the typical coherent boundaries, unveiling larger extents of the coherent cores and showing gradual transition to coherence over ~0.2 pc. 
    more » « less