No Arabic abstract
We undertake a systematic analysis of the early (< 0.5 Myr) evolution of clustering and the stellar initial mass function in turbulent fragmentation simulations. These large scale simulations for the first time offer the opportunity for a statistical analysis of IMF variations and correlations between stellar properties and cluster richness. The typical evolutionary scenario involves star formation in small-n clusters which then progressively merge; the first stars to form are seeds of massive stars and achieve a headstart in mass acquisition. These massive seeds end up in the cores of clusters and a large fraction of new stars of lower mass is formed in the outer parts of the clusters. The resulting clusters are therefore mass segregated at an age of 0.5 Myr, although the signature of mass segregation is weakened during mergers. We find that the resulting IMF has a smaller exponent (alpha=1.8-2.2) than the Salpeter value (alpha=2.35). The IMFs in subclusters are truncated at masses only somewhat larger than the most massive stars (which depends on the richness of the cluster) and an universal upper mass limit of 150 Msun is ruled out. We also find that the simulations show signs of the IGIMF effect proposed by Weidner & Kroupa, where the frequency of massive stars is suppressed in the integrated IMF compared to the IMF in individual clusters. We identify clusters through the use of a minimum spanning tree algorithm which allows easy comparison between observational survey data and the predictions of turbulent fragmentation models. In particular we present quantitative predictions regarding properties such as cluster morphology, degree of mass segregation, upper slope of the IMF and the relation between cluster richness and maximum stellar mass. [abridged]
We use two hydrodynamical simulations (with and without photoionising feedback) of the self-consistent evolution of molecular clouds (MCs) undergoing global hierarchical collapse (GHC), to study the effect of the feedback on the structural and kinematic properties of the gas and the stellar clusters formed in the clouds. During this early stage, the evolution of the two simulations is very similar (implying that the feedback from low mass stars does not affect the cloud-scale evolution significantly) and the star-forming region accretes faster than it can convert gas to stars, causing the instantaneous measured star formation efficiency (SFE) to remain low even in the absence of significant feedback. Afterwards, the ionising feedback first destroys the filamentary supply to star-forming hubs and ultimately removes the gas from it, thus first reducing the star formation (SF) and finally halting it. The ionising feedback also affects the initial kinematics and spatial distribution of the forming stars, because the gas being dispersed continues to form stars, which inherit its motion. In the non-feedback simulation, the groups remain highly compact and do not mix, while in the run with feedback, the gas dispersal causes each group to expand, and the cluster expansion thus consists of the combined expansion of the groups. Most secondary star-forming sites around the main hub are also present in the non-feedback run, implying a primordial rather than triggered nature. We do find one example of a peripheral star-forming site that appears only in the feedback run, thus having a triggered origin. However, this appears to be the exception rather than the rule, although this may be an artifact of our simplified radiative transfer scheme.
Magnetic and energetic properties are presented for 17 dense cores within a few hundred pc of the Sun. Their plane-of-sky field strengths are estimated from the dispersion of polarization directions, following Davis, Chandrasekhar and Fermi (DCF). Their ratio of mass to magnetic critical mass is 0.5-3, indicating nearly critical field strengths. The field strength B_pos is correlated with column density N as B_pos~N^p, where p=1.05+-0.08, and with density n as B_pos~n^q, where q=0.66+-0.05. These magnetic properties are consistent with those derived from Zeeman studies (Crutcher et al. 2010), with less scatter. Relations between virial mass M_V, magnetic critical mass M_B, and Alfven amplitude sigma_B/B match the observed range of M/M_B for cores observed to be nearly virial, with M/M_V=0.5-2, with moderate Alfven amplitudes, and with sigma_B/B=0.1-0.4. The B-N and B-n correlations in the DCF and Zeeman samples can be explained when such bound, Alfvenic, and nearly-critical cores have central concentration and spheroidal shape. For these properties, B~N because M/M_B is nearly constant compared to the range of N, and B~n^(2/3) because M^(1/3) is nearly constant compared to the range of n^(2/3). The observed core fields which follow B~n^(2/3) need not be much weaker than gravity, in contrast to core fields which follow B~n^(2/3) due to spherical contraction at constant mass (Mestel 1966). Instead, the nearly critical values of M/M_B suggest that the observed core fields are nearly as strong as possible, among values which allow gravitational contraction.
In this work, we aim to characterise high-mass clumps with infall motions. We selected 327 clumps from the Millimetre Astronomy Legacy Team 90-GHz (MALT90) survey, and identified 100 infall candidates. Combined with the results of He et al. (2015), we obtained a sample of 732 high-mass clumps, including 231 massive infall candidates and 501 clumps where infall is not detected. Objects in our sample were classified as pre-stellar, proto-stellar, HII or photo-dissociation region (PDR). The detection rates of the infall candidates in the pre-stellar, proto-stellar, HII and PDR stages are 41.2%, 36.6%, 30.6% and 12.7%, respectively. The infall candidates have a higher H$_{2}$ column density and volume density compared with the clumps where infall is not detected at every stage. For the infall candidates, the median values of the infall rates at the pre-stellar, proto-stellar, HII and PDR stages are 2.6$times$10$^{-3}$, 7.0$times$10$^{-3}$, 6.5$times$10$^{-3}$ and 5.5$times$10$^{-3}$ M$_odot$ yr$^{-1}$, respectively. These values indicate that infall candidates at later evolutionary stages are still accumulating material efficiently. It is interesting to find that both infall candidates and clumps where infall is not detected show a clear trend of increasing mass from the pre-stellar to proto-stellar, and to the HII stages. The power indices of the clump mass function (ClMF) are 2.04$pm$0.16 and 2.17$pm$0.31 for the infall candidates and clumps where infall is not detected, respectively, which agree well with the power index of the stellar initial mass function (2.35) and the cold Planck cores (2.0).
The full stellar population of NGC 6334, one of the most spectacular regions of massive star formation in the nearby Galaxy, have not been well-sampled in past studies. We analyze here a mosaic of two Chandra X-ray Observatory images of the region using sensitive data analysis methods, giving a list of 1607 faint X-ray sources with arcsecond positions and approximate line-of-sight absorption. About 95 percent of these are expected to be cluster members, most lower mass pre-main sequence stars. Extrapolating to low X-ray levels, the total stellar population is estimated to be 20-30,000 pre-main sequence stars. The X-ray sources show a complicated spatial pattern with about 10 distinct star clusters. The heavily-obscured clusters are mostly associated with previously known far-infrared sources and radio HII regions. The lightly-obscured clusters are mostly newly identified in the X-ray images. Dozens of likely OB stars are found, both in clusters and dispersed throughout the region, suggesting that star formation in the complex has proceeded over millions of years. A number of extraordinarily heavily absorbed X-ray sources are associated with the active regions of star formation.
NGC 253 hosts the nearest nuclear starburst. Previous observations show a region rich in molecular gas, with dense clouds associated with recent star formation. We used ALMA to image the 350 GHz dust continuum and molecular line emission from this region at 2 pc resolution. Our observations reveal ~14 bright, compact (~2-3 pc FWHM) knots of dust emission. Most of these sources are likely to be forming super star clusters (SSCs) based on their inferred dynamical and gas masses, association with 36 GHz radio continuum emission, and coincidence with line emission tracing dense, excited gas. One source coincides with a known SSC, but the rest remain invisible in Hubble near-infrared (IR) imaging. Our observations imply that gas still constitutes a large fraction of the overall mass in these sources. Their high brightness temperature at 350 GHz also implies a large optical depth near the peak of the IR spectral energy distribution. As a result, these sources may have large IR photospheres and the IR radiation force likely exceeds L/c. Still, their moderate observed velocity dispersions suggest that feedback from radiation, winds, and supernovae are not yet disrupting most sources. This mode of star formation appears to produce a large fraction of stars in the burst. We argue for a scenario in which this phase lasts ~1 Myr, after which the clusters shed their natal cocoons but continue to produce ionizing photons. The strong feedback that drives the observed cold gas and X-ray outflows likely occurs after the clusters emerge from this early phase.