No Arabic abstract
We introduce a theory for the development of a transitional column density $Sigma_{rm TP}$ between the lognormal and the power-law forms of the probability distribution function (PDF) in a molecular cloud. Our turbulent magnetohydrodynamic simulations show that the value of $Sigma_{rm TP}$ increases as the strength of both the initial magnetic field and turbulence increases. We develop an analytic expression for $Sigma_{rm TP}$ based on the interplay of turbulence, a (strong) magnetic field, and gravity. The transition value $Sigma_{rm TP}$ scales with $mathcal{M}^2_{rm 0}$, the square of the initial sonic Mach number, and $beta_{0}$, the initial ratio of gas pressure to magnetic pressure. We fit the variation of $Sigma_{rm TP}$ among different model clouds as a function of $mathcal{M}^2_{rm 0} beta_{0}$, or equivalently the square of the initial Alfvenic Mach number $mathcal{M}^2_{rm A0}$. This implies that the transition value $Sigma_{rm TP}$ is an imprint of cloud initial conditions and is set by turbulent compression of a magnetic cloud. Physically, the value of $Sigma_{rm TP}$ denotes the boundary above which the mass-to-flux ratio becomes supercritical and gravity drives the evolution.
Recent observations of column densities in molecular clouds find lognormal distributions with power-law high-density tails. These results are often interpreted as indications that supersonic turbulence dominates the dynamics of the observed clouds. We calculate and present the column-density distributions of three clouds, modeled with very different techniques, none of which is dominated by supersonic turbulence. The first star-forming cloud is simulated using smoothed particle hydrodynamics (SPH); in this case gravity, opposed only by thermal-pressure forces, drives the evolution. The second cloud is magnetically subcritical with subsonic turbulence, simulated using nonideal MHD; in this case the evolution is due to gravitationally-driven ambipolar diffusion. The third cloud is isothermal, self-gravitating, and has a smooth density distribution analytically approximated with a uniform inner region and an r^-2 profile at larger radii. We show that in all three cases the column-density distributions are lognormal. Power-law tails develop only at late times (or, in the case of the smooth analytic profile, for strongly centrally concentrated configurations), when gravity dominates all opposing forces. It therefore follows that lognormal column-density distributions are generic features of diverse model clouds, and should not be interpreted as being a consequence of supersonic turbulence.
We introduce a new dual power law (DPL) probability distribution function for the mass distribution of stellar and substellar objects at birth, otherwise known as the initial mass function (IMF). The model contains both deterministic and stochastic elements, and provides a unified framework within which to view the formation of brown dwarfs and stars resulting from an accretion process that starts from extremely low mass seeds. It does not depend upon a top down scenario of collapsing (Jeans) masses or an initial lognormal or otherwise IMF-like distribution of seed masses. Like the modified lognormal power law (MLP) distribution, the DPL distribution has a power law at the high mass end, as a result of exponential growth of mass coupled with equally likely stopping of accretion at any time interval. Unlike the MLP, a power law decay also appears at the low mass end of the IMF. This feature is closely connected to the accretion stopping probability rising from an initially low value up to a high value. This might be associated with physical effects of ejections sometimes (i.e., rarely) stopping accretion at early times followed by outflow driven accretion stopping at later times, with the transition happening at a critical time (therefore mass). Comparing the DPL to empirical data, the critical mass is close to the substellar mass limit, suggesting that the onset of nuclear fusion plays an important role in the subsequent accretion history of a young stellar object.
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.
Simulations generally show that non-self-gravitating clouds have a lognormal column density ($Sigma$) probability distribution function (PDF), while self-gravitating clouds with active star formation develop a distinct power-law tail at high column density. Although the growth of the power law can be attributed to gravitational contraction leading to the formation of condensed cores, it is often debated if an observed lognormal shape is a direct consequence of supersonic turbulence alone, or even if it is really observed in molecular clouds. In this paper we run three-dimensional magnetohydrodynamic simulations including ambipolar diffusion with different initial conditions to see the effect of strong magnetic fields and nonlinear initial velocity perturbations on the evolution of the column density PDFs. Our simulations show that column density PDFs of clouds with supercritical mass-to-flux ratio, with either linear perturbations or nonlinear turbulence, quickly develop a power-law tail such that $dN/d log Sigma propto Sigma^{-alpha}$ with index $alpha simeq 2$. Interestingly, clouds with subcritical mass-to-flux ratio also proceed directly to a power-law PDF, but with a much steeper index $alpha simeq 4$. This is a result of gravitationally-driven ambipolar diffusion. However, for nonlinear perturbations with a turbulent spectrum ($v_{k}^{2} propto k^{-4}$), the column density PDFs of subcritical clouds do retain a lognormal shape for a major part of the cloud evolution, and only develop a distinct power-law tail with index $alpha simeq 2$ at greater column density when supercritical pockets are formed.
Both observational and theoretical research over the past decade has demonstrated that the probability distribution function (PDF) of the gas density in turbulent molecular clouds is a key ingredient for understanding star formation. It has recently been argued that the PDF of molecular clouds is a pure power-law distribution. It has been claimed that the log-normal part is ruled out when using only the part of the PDF up/down to which it is complete, that is where the column density contours are still closed. By using the results from high-resolution magnetohydrodynamical simulations of molecular cloud formation and evolution, we find that the column density PDF is indeed composed of a log-normal and, if including self-gravity, a power-law part. We show that insufficient sampling of a molecular cloud results in closed contours that cut off the log-normal part. In contrast, systematically increasing the field of view and sampling the entire cloud yields a completeness limit at the lower column densities, which also recovers the log-normal part. This demonstrates that the field of view must be sufficiently large for the PDF to be complete down to its log-normal part, which has important implications for predictions of star-formation activity based on the PDF.