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 ∝ M Γ , Γ = −1 ± 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Ṁ ∝ 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 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; Γ → −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.

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.

The filamentary nature of accretion streams found around embedded sources suggests that protostellar disks experience heterogenous infall from the star-forming environment, consistent with the accretion behavior onto star-forming cores in top-down star-cluster formation simulations. This may produce disk substructures in the form of rings, gaps, and spirals that continue to be identified by high-resolution imaging surveys in both embedded Class 0/I and later Class II sources. We present a parameter study of anisotropic infall, informed by the properties of accretion flows onto protostellar cores in numerical simulations, and varying the relative specific angular momentum of incoming flows as well as their flow geometry. Our results show that anisotropic infall perturbs the disk and readily launches the Rossby wave instability. It forms vortices at the inner and outer edges of the infall zone where material is deposited. These vortices drive spiral waves and angular momentum transport, with some models able to drive stresses corresponding to a viscosity parameter on the order of α ∼ 10−2. The resulting azimuthal shear forms robust pressure bumps that act as barriers to radial drift of dust grains, as demonstrated by postprocessing calculations of drift-dominated dust evolution. We discuss how a self-consistent model of anisotropic infall can account for the formation of millimeter rings in the outer disk as well as producing compact dust disks, consistent with observations of embedded sources.

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 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 ∼ 0.05 pc (∼ 10 4 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

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.