No Arabic abstract
Numerical simulations of star formation have found that a power-law mass function can develop at high masses. In a previous paper, we employed isothermal simulations which created large numbers of sinks over a large range in masses to show that the power law exponent of the mass function, $dN/dlog M propto M^{Gamma}$, asymptotically and accurately approaches $Gamma = -1.$ Simple analytic models show that such a power law can develop if the mass accretion rate $dot{M} propto M^2$, as in Bondi-Hoyle accretion; however, the sink mass accretion rates in the simulations show significant departures from this relation. In this paper we show that the expected accretion rate dependence is more closely realized provided the gravitating mass is taken to be the sum of the sink mass and the mass in the near environment. This reconciles the observed mass functions with the accretion rate dependencies, and demonstrates that power-law upper mass functions are essentially the result of gravitational focusing, a mechanism present in, for example, the competitive accretion model.
We introduce a new multi-power-law distribution for the Initial Mass Function (IMF) to explore its potential properties. It follows on prior work that introduced mechanisms accounting for mass accretion in star formation, developed within the framework of general evolution equations for the mass distribution of accreting and non-accreting (proto)stars. This paper uses the same fundamental framework to demonstrate that the interplay between a mass-dependent and a time-dependent step-like dropout rate from accretion leads to IMFs that exhibit multiple power laws for an exponential mass growth. While the mass-dependent accretion and its dropout is intrinsic to each star, the time-dependent dropout might be tied to a specific history such as the rapid consumption of nebular material by nearby stars or the sweeping away of some material by shock waves. The time-dependent dropout folded into the mass-dependent process of star formation is shown to have a significant influence on the IMFs.
We derive an analytical theory of the PDF of density fluctuations in supersonic turbulence in the presence of gravity in star-forming clouds. The theory is based on a rigorous derivation of a combination of the Navier-Stokes continuity equations for the fluid motions and the Poisson equation for the gravity. It extends upon previous approaches first by including gravity, second by considering the PDF as a dynamical system, not a stationary one. We derive the transport equations of the density PDF, characterize its evolution and determine the density threshold above which gravity strongly affects and eventually dominates the dynamics of turbulence. We demonstrate the occurence of {it two} power law tails in the PDF, with two characteristic exponents, corresponding to two different stages in the balance between turbulence and gravity. Another important result of this study is to provide a procedure to relate the observed {it column density} PDFs to the corresponding {it volume density} PDFs. This allows to infer, from the observation of column-densities, various physical parameters characterizing molecular clouds, notably the virial parameter. Furthermore, the theory offers the possibility to date the clouds in units of ${t}_{rm coll}$, the time since a statistically significant fraction of the cloud started to collapse. The theoretical results and diagnostics reproduce very well numerical simulations and observations of star-forming clouds. The theory provides a sound theoretical foundation and quantitative diagnostics to analyze observations or numerical simulations of star-forming regions and to characterize the evolution of the density PDF in various regions of molecular clouds.
We report the novel detection of complex high-column density tails in the probability distribution functions (PDFs) for three high-mass star-forming regions (CepOB3, MonR2, NGC6334), obtained from dust emission observed with Herschel. The low column density range can be fit with a lognormal distribution. A first power-law tail starts above an extinction (Av) of ~6-14. It has a slope of alpha=1.3-2 for the rho~r^-alpha profile for an equivalent density distribution (spherical or cylindrical geometry), and is thus consistent with free-fall gravitational collapse. Above Av~40, 60, and 140, we detect an excess that can be fitted by a flatter power law tail with alpha>2. It correlates with the central regions of the cloud (ridges/hubs) of size ~1 pc and densities above 10^4 cm^-3. This excess may be caused by physical processes that slow down collapse and reduce the flow of mass towards higher densities. Possible are: 1. rotation, which introduces an angular momentum barrier, 2. increasing optical depth and weaker cooling, 3. magnetic fields, 4. geometrical effects, and 5. protostellar feedback. The excess/second power-law tail is closely linked to high-mass star-formation though it does not imply a universal column density threshold for the formation of (high-mass) stars.
In the star formation process, the vital impact of environmental factors such as feedback from massive stars and stellar density on the form of the initial mass function (IMF) at low-mass end is yet to be understood. Hence a systematic, highly sensitive observational analysis of a sample of regions under diverse environmental conditions is essential. We analyse the IMF of eight young clusters ($<$5 Myr), namely IC1848-West, IC1848-East, NGC 1893, NGC 2244, NGC 2362, NGC 6611, Stock 8 and Cygnus OB2, which are located at the Galactocentric distance ($R_g$) range $sim$6-12 kpc along with nearby cluster IC348 using deep near-IR photometry and Gaia DR2. These clusters are embedded in massive stellar environments of radiation strength $log(L_{FUV}/L_{odot})$ $sim$2.6 to 6.8, $log(L_{EUV})$ $sim$42.2 to 50.85 photons/s, with stellar density in the range of $sim$170 - 1220 stars/pc$^2$. After structural analysis and field decontamination we obtain an unbiased, uniformly sensitive sample of pre-main-sequence members of the clusters down to brown-dwarf regime. The lognormal fit to the IMF of nine clusters gives the mean characteristic mass ($m_c$) and $sigma$ of 0.32$pm$0.02 $M_odot$ and 0.47$pm$0.02, respectively. We compare the IMF with that of low- and high-mass clusters across the Milky Way. We also check for any systematic variation with respect to the radiation field strength, stellar density as well with $R_g$. We conclude that there is no strong evidence for environmental effect in the underlying form of the IMF of these clusters.
We study the star formation (SF) law in 12 Galactic molecular clouds with ongoing high-mass star formation (HMSF) activity, as traced by the presence of a bright IRAS source and other HMSF tracers. We define the molecular cloud (MC) associated to each IRAS source using 13CO line emission, and count the young stellar objects (YSOs) within these clouds using GLIMPSE and MIPSGAL 24 micron Spitzer databases.The masses for high luminosity YSOs (Lbol>10~Lsun) are determined individually using Pre Main Sequence evolutionary tracks and the evolutionary stages of the sources, whereas a mean mass of 0.5 Msun was adopted to determine the masses in the low luminosity YSO population. The star formation rate surface density (sigsfr) corresponding to a gas surface density (siggas) in each MC is obtained by counting the number of the YSOs within successive contours of 13CO line emission. We find a break in the relation between sigsfr and siggas, with the relation being power-law (sigsfr ~ siggas^N) with the index N varying between 1.4 and 3.6 above the break. The siggas at the break is between 150-360 Msun/pc^2 for the sample clouds, which compares well with the threshold gas density found in recent studies of Galactic star-forming regions. Our clouds treated as a whole lie between the Kennicutt (1998) relation and the linear relation for Galactic and extra-galactic dense star-forming regions. We find a tendency for the high-mass YSOs to be found preferentially in dense regions at densities higher than 1200 Msun/pc^2 (~0.25 g/cm^2).