No Arabic abstract
Star formation in a filamentary infrared dark cloud (IRDC) is simulated over a dynamic range of 4.2 pc to 28 au for a period of $3.5times 10^5$ yr, including magnetic fields and both radiative and outflow feedback from the protostars. At the end of the simulation, the star formation efficiency is 4.3 per cent and the star formation rate per free fall time is $epsilon_{rm ff}simeq 0.04$, within the range of observed values (Krumholz et al. 2012a). The total stellar mass increases as $sim,t^2$, whereas the number of protostars increases as $sim,t^{1.5}$. We find that the density profile around most of the simulated protostars is $sim,rhopropto r^{-1.5}$, as predicted by Murray & Chang (2015). At the end of the simulation, the protostellar mass function approaches the Chabrier (2005) stellar initial mass function. We infer that the time to form a star of median mass $0.2,M_odot$ is about $1.4times 10^5$~yr from the median mass accretion rate. We find good agreement among the protostellar luminosities observed in the large sample of Dunham et al. (2013), our simulation, and a theoretical estimate, and conclude that the classical protostellar luminosity problem Kenyon et al. (1990) is resolved. The multiplicity of the stellar systems in the simulation agrees to within a factor 2 of observations of Class I young stellar objects; most of the simulated multiple systems are unbound. Bipolar protostellar outflows are launched using a sub-grid model, and extend up to 1 pc from their host star. The mass-velocity relation of the simulated outflows is consistent with both observation and theory.
We perform ideal MHD high resolution AMR simulations with driven turbulence and self-gravity and find that long filamentary molecular clouds are formed at the converging locations of large-scale turbulence flows and the filaments are bounded by gravity. The magnetic field helps shape and reinforce the long filamentary structures. The main filamentary cloud has a length of ~4.4 pc. Instead of a monolithic cylindrical structure, the main cloud is shown to be a collection of fiber/web-like sub-structures similar to filamentary clouds such as L1495. Unless the line-of-sight is close to the mean field direction, the large-scale magnetic field and striations in the simulation are found roughly perpendicular to the long axis of the main cloud, similar to 1495. This provides strong support for a large-scale moderately strong magnetic field surrounding L1495. We find that the projection effect from observations can lead to incorrect interpretations of the true three-dimensional physical shape, size, and velocity structure of the clouds. Helical magnetic field structures found around filamentary clouds that are interpreted from Zeeman observations can be explained by a simple bending of the magnetic field that pierces through the cloud. We demonstrate that two dark clouds form a T-shape configuration which are strikingly similar to the Infrared dark cloud SDC13 leading to the interpretation that SDC13 results from a collision of two long filamentary clouds. We show that a moderately strong magnetic field (M_A ~ 1) is crucial for maintaining a long and slender filamentary cloud for a long period of time ~0.5 million years.
To study the early phases of massive star formation, we present ALMA observations of SiO(5-4) emission and VLA observations of 6 cm continuum emission towards 32 Infrared Dark Cloud (IRDC) clumps, spatially resolved down to $lesssim 0.05$ pc. Out of the 32 clumps, we detect SiO emission in 20 clumps, and in 11 of them the SiO emission is relatively strong and likely tracing protostellar outflows. Some SiO outflows are collimated, while others are less ordered. For the six strongest SiO outflows, we estimate basic outflow properties. In our entire sample, where there is SiO emission, we find 1.3 mm continuum and infrared emission nearby, but not vice versa. We build the spectral energy distributions (SEDs) of cores with 1.3 mm continuum emission and fit them with radiative transfer (RT) models. The low luminosities and stellar masses returned by SED fitting suggest these are early stage protostars. We see a slight trend of increasing SiO line luminosity with bolometric luminosity, which suggests more powerful shocks in the vicinity of more massive YSOs. We do not see a clear relation between the SiO luminosity and the evolutionary stage indicated by $L/M$. We conclude that as a protostar approaches a bolometric luminosity of $sim 10^2 : L_{odot}$, the shocks in the outflow are generally strong enough to form SiO emission. The VLA 6 cm observations toward the 15 clumps with the strongest SiO emission detect emission in four clumps, which is likely shock ionized jets associated with the more massive ones of these protostellar cores.
We present molecular line observations, made with angular resolutions of ~20, toward the filamentary infrared dark cloud G34.43+0.24 using the APEX [CO(3-2), 13CO(3-2), C18O(3-2) and CS(7-6) transitions], Nobeyama 45 m [CS(2-1), SiO(2-1), C34S(2-1), HCO+(1-0), H13CO+(1-0) and CH3OH(2-1) transitions], and SEST [CS(2-1) and C18O(2-1) transitions] telescopes. We find that the spatial distribution of the molecular emission is similar to that of the dust continuum emission observed with 11 resolution showing a filamentary structure and four cores. The cores have local thermodynamic equilibrium masses ranging from 3.3x10^2 - 1.5x10^3 solar masses and virial masses from 1.1x10^3 - 1.5x10^3 solar masses, molecular hydrogen densities between 1.8x10^4 and 3.9x10^5 cm^{-3}, and column densities >2.0x10^{22} cm^{-2}; values characteristics of massive star forming cores. The 13CO(3-2) profile observed toward the most massive core reveals a blue profile indicating that the core is undergoing large-scale inward motion with an average infall velocity of 1.3 km/s and a mass infall rate of 1.8x10^{-3} solar masses per year. We report the discovery of a molecular outflow toward the northernmost core thought to be in a very early stage of evolution. We also detect the presence of high velocity gas toward each of the other three cores, giving support to the hypothesis that the excess 4.5 $mu$ emission (green fuzzies) detected toward these cores is due to shocked gas. The molecular outflows are massive and energetic, with masses ranging from 25 -- 80 solar masses, momentum 2.3 - 6.9x10^2 Msun km/s, and kinetic energies 1.1 - 3.6x10^3 Msun km^2 s^{-2}; indicating that they are driven by luminous, high-mass young stellar objects.
We present a large suite of MHD simulations of turbulent, star-forming giant molecular clouds(GMCs) with stellar feedback, extending previous work by simulating 10 different random realizations for each point in the parameter space of cloud mass and size. It is found that oncethe clouds disperse due to stellar feedback, both self-gravitating star clusters and unbound stars generally remain, which arise from the same underlying continuum of substructured stellar density, ie. the hierarchical cluster formation scenario. The fraction of stars that are born within gravitationally-bound star clusters is related to the overall cloud star formation efficiency set by stellar feedback, but has significant scatter due to stochastic variations in the small-scale details of the star-forming gas flow. We use our numerical results to calibrate a model for mapping the bulk properties (mass, size, and metallicity) of self-gravitating GMCs onto the star cluster populations they form, expressed statistically in terms of cloud-level distributions. Synthesizing cluster catalogues from an observed GMC catalogue in M83, we find that this model predicts initial star cluster masses and sizes that are in good agreement with observations, using only standard IMF and stellar evolution models as inputs for feedback. Within our model, the ratio of the strength of gravity to stellar feedback is the key parameter setting the masses of star clusters, and of the various feedback channels direct stellar radiation(photon momentum and photoionization) is the most important on GMC scales.
Infrared dark clouds (IRDCs) are thought to be potential hosts of the elusive early phases of high-mass star formation. Here we conduct an in-depth kinematic analysis of one such IRDC, G034.43+00.24 (Cloud F), using high sensitivity and high spectral resolution IRAM-30m N$_2$H$^+$ ($1-0$) and C$^{18}$O ($1-0$) observations. To disentangle the complex velocity structure within this cloud we use Gaussian decomposition and hierarchical clustering algorithms. We find that four distinct coherent velocity components are present within Cloud F. The properties of these components are compared to those found in a similar IRDC, G035.39-00.33 (Cloud H). We find that the components in both clouds have: high densities (inferred by their identification in N$_2$H$^+$), trans-to-supersonic non-thermal velocity dispersions with Mach numbers of $sim$ $1.5-4$, a separation in velocity of $sim$3 km s$^{-1}$, and a mean red-shift of $sim$ 0.3 km s$^{-1}$ between the N$_2$H$^+$ (dense gas) and C$^{18}$O emission (envelope gas). The latter of these could suggest that these clouds share a common formation scenario. We investigate the kinematics of the larger-scale Cloud F structures, using lower-density-tracing $^{13}$CO ($1-0$) observations. A good correspondence is found between the components identified in the IRAM-30m observations and the most prominent component in the$^{13}$CO data. We find that the IRDC Cloud F is only a small part of a much larger structure, which appears to be an inter-arm filament of the Milky Way.