No Arabic abstract
Recent observations of molecular clouds show that dense filaments are the sites of present-day star formation. Thus, it is necessary to understand the filament formation process because these filaments provide the initial condition for star formation. Theoretical research suggests that shock waves in molecular clouds trigger filament formation. Since several different mechanisms have been proposed for filament formation, the formation mechanism of the observed star-forming filaments requires clarification. In the present study, we perform a series of isothermal magnetohydrodynamics simulations of filament formation. We focus on the influences of shock velocity and turbulence on the formation mechanism and identified three different mechanisms for the filament formation. The results indicate that when the shock is fast, at shock velocity v_sh = 7 km/s, the gas flows driven by the curved shock wave create filaments irrespective of the presence of turbulence and self-gravity. However, at a slow shock velocity v_sh = 2.5 km/s, the compressive flow component involved in the initial turbulence induces filament formation. When both the shock velocities and turbulence are low, the self-gravity in the shock-compressed sheet becomes important for filament formation. Moreover, we analyzed the line-mass distribution of the filaments and showed that strong shock waves can naturally create high-line-mass filaments such as those observed in the massive star-forming regions in a short time. We conclude that the dominant filament formation mode changes with the velocity of the shock wave triggering the filament formation.
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.
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 test some ideas for star formation relations against data on local molecular clouds. On a cloud by cloud basis, the relation between the surface density of star formation rate and surface density of gas divided by a free-fall time, calculated from the mean cloud density, shows no significant correlation. If a crossing time is substituted for the free-fall time, there is even less correlation. Within a cloud, the star formation rate volume and surface densities increase rapidly with the corresponding gas densities, faster than predicted by models using the free-fall time defined from the local density. A model in which the star formation rate depends linearly on the mass of gas above a visual extinction of 8 mag describes the data on these clouds, with very low dispersion. The data on regions of very massive star formation, with improved star formation rates based on free-free emission from ionized gas, also agree with this linear relation.
Dust is the usual minor component of the interstellar medium. Its dynamic role in the contraction of the diffuse gas into molecular clouds is commonly assumed to be negligible because of the small mass fraction, $f simeq 0.01$. However, as shown in this study, the collective motion of dust grains with respect to the gas may considerably contribute to the destabilisation of the medium on scales $lambda lesssim lambda_J$, where $lambda_J$ is the Jeans length-scale. The linear perturbations of the uniform self-gravitating gas at rest are marginally stable at $lambda simeq lambda_J$, but as soon as the drift of grains is taken into account, they begin growing at a rate approximately equal to $(f tau)^{1/3} t^{-1}_{ff}$, where $tau$ is the stopping time of grains expressed in units of the free fall time of the cloud, $t_{ff}$. The physical mechanism responsible for such a weak dependence of the growth rate on $f$ is the resonance of heavy sound waves stopped by the self-gravity of gas with weak gravitational attraction caused by perturbations of the dust fraction. Once there is stationary subsonic bulk drift of the dust, the growing gas-dust perturbations at $lambda < lambda_J$ become waves propagating with the drift velocity projected onto the wavevector. Their growth has a resonant nature as well and the growth rate is substantially larger than that of the recently discovered resonant instability of gas-dust mixture in the absence of self-gravity. The new instabilities can facilitate gravitational contraction of cold interstellar gas into clouds and additionally produce dusty domains of sub-Jeans size at different stages of molecular cloud formation and evolution.
Previous studies have shown that star formation depends on the driving of molecular cloud turbulence, and differences in the driving can produce an order of magnitude difference in the star formation rate. The turbulent driving is characterised by the parameter $zeta$, with $zeta=0$ for compressive, curl-free driving (e.g. accretion or supernova explosions), and $zeta=1$ for solenoidal, divergence-free driving (e.g. Galactic shear). Here we develop a new method to measure $zeta$ from observations of synchrotron emission from molecular clouds. We calculate statistics of mock synchrotron intensity images produced from magnetohydrodynamic simulations of molecular clouds, in which the driving was controlled to produce different values of $zeta$. We find that the mean and standard deviation of the log-normalised synchrotron intensity are sensitive to $zeta$, for values of $zeta$ between $0$ (curl-free driving) and $0.5$ (naturally-mixed driving). We quantify the dependence of zeta on the direction of the magnetic field relative to the line of sight. We provide best-fit formulae for $zeta$ in terms of the log-normalised mean and standard deviation of synchrotron intensity, with which $zeta$ can be determined for molecular clouds that have similar Alfvenic Mach number to our simulations. These formulae are independent of the sonic Mach number. Signal-to-noise ratios larger than $5$, and angular resolutions smaller than $5%$ of the cloud diameter, are required to apply these formulae. Although there are no firm detections of synchrotron emission from molecular clouds, by combining Green Bank Telescope and Very Large Array observations it should be possible to detect synchrotron emission from molecular clouds, thereby constraining the value of $zeta$.