Building on our previous hydrodynamic study of the angular momenta of cloud cores formed during gravitational collapse of star-forming molecular gas in our previous work, we now examine core properties assuming ideal magnetohydrodynamics (MHD). Using the same sink-patch implementation for the \emph{Athena} MHD code, we characterize the statistical properties of cores, including the mass accretion rates, specific angular momenta, and alignments between the magnetic field and the spin axis of the core on the $0.1 \ \mathrm{pc}$ scale. Our simulations, which reproduce the observed relation between magnetic field strength and gas density, show that magnetic fields can help collimate low density flows and help seed the locations of filamentary structures. Consistent with our previous purely hydrodynamic simulations, stars (sinks) form within the heterogeneous environments of filaments, such that accretion onto cores is highly episodic leading to short-term variability but no long-term monotonic growth of the specific angular momenta. With statistical characterization of protostellar cores properties and behaviors, we aim to provide a starting point for building more realistic and self-consistent disk formation models, helping to address whether magnetic fields can prevent the development of (large) circumstellar disks in the ideal MHD limit.

We present the results of a suite of numerical simulations designed to explore the origin of the angular momenta of protostellar cores. Using the hydrodynamic grid code \emph{Athena} with a sink implementation, we follow the formation of protostellar cores and protostars (sinks) from the subvirial collapse of molecular clouds on larger scales to investigate the range and relative distribution of core properties. We find that the core angular momenta are relatively unaffected by large-scale rotation of the parent cloud; instead, we infer that angular momenta are mainly imparted by torques between neighboring mass concentrations and exhibit a log-normal distribution. Our current simulation results are limited to size scales $\sim 0.05$~pc ($\sim 10^4 \rm AU$), but serve as first steps toward the ultimate goal of providing initial conditions for higher-resolution studies of core collapse to form protoplanetary disks.

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/d\log 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 use numerical simulations of molecular cloud formation in the colliding flow scenario to investigate the reliability of observational estimates of the angular momenta of early-state, low-mass protostellar cores. We show that, with suitable corrections for projection factors, molecular line observations of velocity gradients in NH$_3$ can be used to provide reasonable estimates of core angular momenta within a factor of two to three, with a few large underestimates due to unfavorable viewing angles. Our results differ from previous investigations which suggested that observations might overestimate true angular momenta by as much as an order of magnitude; the difference is probably due to the much smoother velocity field on small scales in our simulations, which result from allowing turbulence to decay and gravitational infall to dominate. The results emphasize the importance of understanding the nature of 'turbulent' velocities, with implications for the formation of protostellar disks during core collapse.

The origin of the observed morphological and kinematic substructure of young star forming regions is a matter of debate. We offer a new analysis of data from simulations of globally gravitationally collapsing clouds of progenitor gas to answer questions about sub-structured star formation in the context of cold collapse. As a specific example, we compare our models to recent radial velocity survey data from the IN-SYNC survey of Orion and new observations of dense gas kinematics, and offer possible interpretations of kinematic and morphological signatures in the region. In the context of our model, we find the frequently-observed hub-filament morphology of the gas naturally arises during gravitational evolution, as well as the dynamically-distinct kinematic substructure of stars. We emphasize that the global and not just the local gravitational potential plays an important role in determining the dynamics of both clusters and filaments.

We discuss the possibility that gravitational focusing, is responsible for the power-law mass function of star clusters $N(\log M) \propto M^{-1}$. This power law can be produced asymptotically when the mass accretion rate of an object depends upon the mass of the accreting body as $\dot{M} \propto M^2$. While Bondi-Hoyle-Littleton accretion formally produces this dependence on mass in a uniform medium, realistic environments are much more complicated. However, numerical simulations in SPH allowing for sink formation yield such an asymptotic power-law mass function. We perform pure N-body simulations to isolate the effects of gravity from those of gas physics and to show that clusters naturally result with the power-law mass distribution. We also consider the physical conditions necessary to produce clusters on appropriate timescales. Our results help support the idea that gravitationally-dominated accretion is the most likely mechanism for producing the cluster mass function.

Sub-virial gravitational collapse is one mechanism by which star clusters may form. Here we investigate whether this mechanism can be inferred from observations of young clusters. To address this question, we have computed SPH simulations of the initial formation and evolution of a dynamically young star cluster through cold (sub-virial) collapse, starting with an ellipsoidal, turbulently seeded distribution of gas, and forming sink particles representing (proto)stars. While the initial density distributions of the clouds do not have large initial mass concentrations, gravitational focusing due to the global morphology leads to cluster formation. We use the resulting structures to extract observable morphological and kinematic signatures for the case of sub-virial collapse. We find that the signatures of the initial conditions can be erased rapidly as the gas and stars collapse, suggesting that kinematic observations need to be made either early in cluster formation and/or at larger scales, away from the growing cluster core. Our results emphasize that a dynamically young system is inherently evolving on short timescales, so that it can be highly misleading to use current-epoch conditions to study aspects such as star formation rates as a function of local density. Our simulations serve as a starting point for further studies of collapse including other factors such as magnetic fields and stellar feedback.

We report on a series of numerical simulations of gas clouds with self-gravity forming sink particles, adopting an isothermal equation of state to isolate the effects of gravity from thermal physics on the resulting sink mass distributions. Simulations starting with supersonic velocity fluctuations develop sink mass functions with a high-mass power-law tail $dN/d\log M \propto M^{\Gamma}$, $\Gamma = -1 \pm 0.1$, independent of the initial Mach number of the velocity field. Similar results but with weaker statistical significance hold for a simulation starting with initial density fluctuations. This mass function power-law dependence agrees with the asymptotic limit found by Zinnecker assuming Bondi-Hoyle-Littleton (BHL) accretion, even though the mass accretion rates of individual sinks show significant departures from the predicted $\mdot \propto M^2$ behavior. While BHL accretion is not strictly applicable due to the complexity of the environment, we argue that the final mass functions are the result of a {\em relative} $M^2$ dependence resulting from gravitationally-focused accretion. Our simulations may show the power-law mass function particularly clearly compared with others because our adoption of an isothermal equation of state limits the effects of thermal physics in producing a broad initial fragmentation spectrum; $\Gamma \rightarrow -1$ is an asymptotic limit found only when sink masses grow well beyond their initial values. While we have purposely eliminated many additional physical processes (radiative transfer, feedback) which can affect the stellar mass function, our results emphasize the importance of gravitational focusing for massive star formation.