No Arabic abstract
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.
The most accurate measurements of magnetic fields in star-forming gas are based on the Zeeman observations analyzed by Crutcher et al. (2010). We show that their finding that the 3D magnetic field scales approximately as density$^{0.65}$ can also be obtained from analysis of the observed line-of-sight fields. We present two large-scale AMR MHD simulations of several thousand $M_odot$ of turbulent, isothermal, self-gravitating gas, one with a strong initial magnetic field (Alfven Mach number $M_{A,0}= 1$) and one with a weak initial field ($M_{A,0}=10$). We construct samples of the 100 most massive clumps in each simulation and show that they exhibit a power-law relation between field strength and density in excellent agreement with the observed one. Our results imply that the average field in molecular clumps in the interstellar medium is $<B_{tot}> sim 42 n_{H,4}^{0.65} mu$G. Furthermore, the median value of the ratio of the line-of-sight field to density$^{0.65}$ in the simulations is within a factor of about (1.3, 1.7) of the observed value for the strong and weak field cases, respectively. The median value of the mass-to-flux ratio, normalized to the critical value, is 70% of the line-of-sight value. This is larger than the 50% usually cited for spherical clouds because the actual mass-to-flux ratio depends on the volume-weighted field, whereas the observed one depends on the mass-weighted field. Our results indicate that the typical molecular clump in the ISM is significantly supercritical (~ factor of 3). The results of our strong-field model are in very good quantitative agreement with the observations of Li et al. (2009), which show a strong correlation in field orientation between small and large scales. Because there is a negligible correlation in the weak-field model, we conclude that molecular clouds form from strongly magnetized (although magnetically supercritical) gas.
Recent surveys of dust continuum emission at sub-mm wavelengths have shown that filamentary molecular clouds are ubiquitous along the Galactic plane. These structures are inhomogeneous, with over-densities that are sometimes associated with infrared emission and active of star formation. To investigate the connection between filaments and star formation, requires an understanding of the processes that lead to the fragmentation of filaments and a determination of the physical properties of the over-densities (clumps). In this paper, we present a multi-wavelength study of five filamentary molecular clouds, containing several clumps in different evolutionary stages of star formation. We analyse the fragmentation of the filaments and derive the physical properties of their clumps. We find that the clumps in all filaments have a characteristic spacing consistent with the prediction of the `sausage instability theory, regardless of the complex morphology of the filaments or their evolutionary stage. We also find that most clumps have sufficient mass and density to form high-mass stars, supporting the idea that high-mass stars and clusters form within filaments.
The structure of molecular clouds (MCs) holds important clues on the physical processes that lead to their formation and subsequent evolution. While it is well established that turbulence imprints a self-similar structure to the clouds, other processes, such as gravity and stellar feedback, can break their scale-free nature. The break of self-similarity can manifest itself in the existence of characteristic scales that stand out from the underlying structure generated by turbulent motions. We investigate the structure of the Cygnus-X North and the Polaris MCs which represent two extremes in terms of their star formation activity. We characterize the structure of the clouds using the delta-variance ($Delta$-variance) spectrum. In Polaris, the structure of the cloud is self-similar over more than one order of magnitude in spatial scales. In contrast, the $Delta$-variance spectrum of Cygnus-X exhibits an excess and a plateau on physical scales of ~0.5-1.2 pc. In order to explain the observations for Cygnus-X, we use synthetic maps in which we overlay populations of discrete structures on top of a fractal Brownian motion (fBm) image. The properties of these structures such as their major axis sizes, aspect ratios, and column density contrasts are randomly drawn from parameterized distribution functions. We show that it is possible to reproduce a $Delta$-variance spectrum that resembles the one of the Cygnus-X cloud. We also use a reverse engineering approach in which we extract the compact structures in the Cygnus-X cloud and re-inject them on an fBm map. The calculated $Delta$-variance using this approach deviates from the observations and is an indication that the range of characteristic scales observed in Cygnus-X is not only due to the existence of compact sources, but is a signature of the whole population of structures, including more extended and elongated structures
We show that the inter-cloud Larson scaling relation between mean volume density and size $rhopropto R^{-1}$, which in turn implies that mass $Mpropto R^2$, or that the column density $N$ is constant, is an artifact of the observational methods used. Specifically, setting the column density threshold near or above the peak of the column density probability distribution function Npdf ($Nsim 10^{21}$ cmalamenos 2) produces the Larson scaling as long as the Npdf decreases rapidly at higher column densities. We argue that the physical reasons behind local clouds to have this behavior are that (1) this peak column density is near the value required to shield CO from photodissociation in the solar neighborhood, and (2) gas at higher column densities is rare because it is susceptible to gravitational collapse into much smaller structures in specific small regions of the cloud. Similarly, we also use previous results to show that if instead a threshold is set for the volume density, the density will appear to be constant, implying thus that $M propto R^3$. Thus, the Larson scaling relation does not provide much information on the structure of molecular clouds, and does not imply either that clouds are in Virial equilibrium, or have a universal structure. We also show that the slope of the $M-R$ curve for a single cloud, which transitions from near-to-flat values for large radii to $alpha=2$ as a limiting case for small radii, depends on the properties of the Npdf.
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.