Physical properties of the star-forming clusters in NGC 6334

Aims. We aim to characterise certain physical properties of high-mass star-forming sites in the NGC 6334 molecular cloud, such as the core mass function (CMF), spatial distribution of cores, and mass segregation. Methods. We used the Atacama Large Millimeter/sub-millimeter Array (ALMA) to image the embedded clusters NGC 6334-I and NGC 6334-I(N) in the continuum emission at 87.6 GHz. We achieved a spatial resolution of 1300 au, enough to resolve different compact cores and fragments, and to study the properties of the clusters. Results. We detected 142 compact sources distributed over the whole surveyed area. The ALMA compact sources are clustered in different regions. We used different machine-learning algorithms to identify four main clusters: NGC 6334-I, NGC 6334-I(N), NGC 6334-I(NW), and NGC 6334-E. The typical separations between cluster members range from 4000 au to 12 000 au. These separations, together with the core masses (0.1–100 M⊙), are in agreement with the fragmentation being controlled by turbulence at scales of 0.1 pc. We find that the CMFs show an apparent excess of high-mass cores compared to the stellar initial mass function. We evaluated the effects of temperature and unresolved multiplicity on the derived slope of the CMF. Based on this, we conclude that the excess of high-mass cores might be spurious and due to inaccurate temperature determinations and/or resolution limitations. We searched for evidence of mass segregation in the clusters and we find that clusters NGC 6334-I and NGC 6334-I(N) show hints of segregation with the most massive cores located in the centre of the clusters. Conclusions. We searched for correlations between the physical properties of the four embedded clusters and their evolutionary stage (based on the presence of H II regions and infrared sources). NGC 6334-E appears as the most evolved cluster, already harbouring a well-developed H II region. NGC 6334-I is the second-most evolved cluster with an ultra-compact H II region. NGC 6334-I(N) contains the largest population of dust cores distributed in two filamentary structures and no dominant H II region. Finally, NGC 6334-I(NW) is a cluster of mainly low-mass dust cores with no clear signs of massive cores or H II regions. We find a larger separation between cluster members in the more evolved clusters favouring the role of gas expulsion and stellar ejection with evolution. The mass segregation, seen in the NGC 6334-I and NGC 6334-I(N) clusters, suggests a primordial origin for NGC 6334-I(N). In contrast, the segregation in NGC 6334-I might be due to dynamical effects. Finally, the lack of massive cores in the most evolved cluster suggests that the gas reservoir is already exhausted, while the less evolved clusters still have a large gas reservoir along with the presence of massive cores. In general, the fragmentation process of NGC 6334 at large scales (from filament to clump, i.e. at about 1 pc) is likely governed by turbulent pressure, while at smaller scales (scale of cores and sub-fragments, i.e. a few hundred au) thermal pressure starts to be more significant.

On the mass segregation of cores and stars

ABSTRACT Observations of pre- and proto-stellar cores in young star-forming regions show them to be mass segregated, i.e. the most massive cores are centrally concentrated, whereas pre-main-sequence stars in the same star-forming regions (and older regions) are not. We test whether this apparent contradiction can be explained by the massive cores fragmenting into stars of much lower mass, thereby washing out any signature of mass segregation in pre-main-sequence stars. Whilst our fragmentation model can reproduce the stellar initial mass function, we find that the resultant distribution of pre-main sequence stars is mass segregated to an even higher degree than that of the cores, because massive cores still produce massive stars if the number of fragments is reasonably low (between one and five). We therefore suggest that the reason cores are observed to be mass segregated and stars are not is likely due to dynamical evolution of the stars, which can move significant distances in star-forming regions after their formation.

Star formation activity and the spatial distribution and mass segregation of dense cores in the early phases of star formation

We examine the spatial distribution and mass segregation of dense molecular cloud cores in a number of nearby star forming regions (the region L1495 in Taurus, Aquila, Corona Australis, and W43) that span about four orders of magnitude in star formation activity. We used an approach based on the calculation of the minimum spanning tree, and for each region, we calculated the structure parameter 𝒬 and the mass segregation ratio ΛMSR measured for various numbers of the most massive cores. Our results indicate that the distribution of dense cores in young star forming regions is very substructured and that it is very likely that this substructure will be imprinted onto the nascent clusters that will emerge out of these clouds. With the exception of Taurus in which there is nearly no mass segregation, we observe mild-to-significant levels of mass segregation for the ensemble of the 6, 10, and 14 most massive cores in Aquila, Corona Australis, and W43, respectively. Our results suggest that the clouds’ star formation activity are linked to their structure, as traced by their population of dense cores. We also find that the fraction of massive cores that are the most mass segregated in each region correlates with the surface density of star formation in the clouds. The Taurus region with low star forming activity is associated with a highly hierarchical spatial distribution of the cores (low 𝒬 value) and the cores show no sign of being mass segregated. On the other extreme, the mini-starburst region W43-MM1 has a higher 𝒬 that is suggestive of a more centrally condensed structure. Additionally, it possesses a higher fraction of massive cores that are segregated by mass. While some limited evolutionary effects might be present, we largely attribute the correlation between the star formation activity of the clouds and their structure to a dependence on the physical conditions that have been imprinted on them by the large scale environment at the time they started to assemble.

Planetesimal fragmentation and giant planet formation

Context. Most planet formation models that incorporate planetesimal fragmentation consider a catastrophic impact energy threshold for basalts at a constant velocity of 3 km s−1 throughout the process of the formation of the planets. However, as planets grow, the relative velocities of the surrounding planetesimals increase from velocities of the order of meters per second to a few kilometers per second. In addition, beyond the ice line where giant planets are formed, planetesimals are expected to be composed roughly of 50% ices. Aims. We aim to study the role of planetesimal fragmentation on giant planet formation considering the planetesimal catastrophic impact energy threshold as a function of the planetesimal relative velocities and compositions. Methods. We improved our model of planetesimal fragmentation incorporating a functional form of the catastrophic impact energy threshold with the planetesimal relative velocities and compositions. We also improved in our model the accretion of small fragments produced by the fragmentation of planetesimals during the collisional cascade considering specific pebble accretion rates. Results. We find that a more accurate and realistic model for the calculation of the catastrophic impact energy threshold tends to slow down the formation of massive cores. Only for reduced grain opacity values at the envelope of the planet is the cross-over mass achieved before the disk timescale dissipation. Conclusions. While planetesimal fragmentation favors the quick formation of massive cores of 5–10 M⊕ the cross-over mass could be inhibited by planetesimal fragmentation. However, grain opacity reduction or pollution by the accreted planetesimals together with planetesimal fragmentation could explain the formation of giant planets with low-mass cores.

Gravity drives the evolution of infrared dark hubs: JVLA observations of SDC13

Context. Converging networks of interstellar filaments, that is hubs, have been recently linked to the formation of stellar clusters and massive stars. Understanding the relationship between the evolution of these systems and the formation of cores and stars inside them is at the heart of current star formation research. Aims. The goal is to study the kinematic and density structure of the SDC13 prototypical hub at high angular resolution to determine what drives its evolution and fragmentation. Methods. We have mapped SDC13, a ~1000 M⊙ infrared dark hub, in NH3(1,1) and NH3(2,2) emission lines, with both the Jansky Very Large Array and Green Bank Telescope. The high angular resolution achieved in the combined dataset allowed us to probe scales down to 0.07 pc. After fitting the ammonia lines, we computed the integrated intensities, centroid velocities and line widths, along with gas temperatures and H2 column densities. Results. The mass-per-unit-lengths of all four hub filaments are thermally super-critical, consistent with the presence of tens of gravitationally bound cores identified along them. These cores exhibit a regular separation of ~0.37 ± 0.16 pc suggesting gravitational instabilities running along these super-critical filaments are responsible for their fragmentation. The observed local increase of the dense gas velocity dispersion towards starless cores is believed to be a consequence of such fragmentation process. Using energy conservation arguments, we estimate that the gravitational to kinetic energy conversion efficiency in the SDC13 cores is ~35%. We see velocity gradient peaks towards ~63% of cores as expected during the early stages of filament fragmentation. Another clear observational signature is the presence of the most massive cores at the filaments’ junction, where the velocity dispersion is largest. We interpret this as the result of the hub morphology generating the largest acceleration gradients near the hub centre. Conclusions. We propose a scenario for the evolution of the SDC13 hub in which filaments first form as post-shock structures in a supersonic turbulent flow. As a result of the turbulent energy dissipation in the shock, the dense gas within the filaments is initially mostly sub-sonic. Then gravity takes over and starts shaping the evolution of the hub, both fragmenting filaments and pulling the gas towards the centre of the gravitational well. By doing so, gravitational energy is converted into kinetic energy in both local (cores) and global (hub centre) potential well minima. Furthermore, the generation of larger gravitational acceleration gradients at the filament junctions promotes the formation of more massive cores.