Radiation-Hydrodynamic Simulations of Protostellar Outflows: Synthetic Observations and Data Comparisons

We present results from three-dimensional, self-gravitating, radiation-hydrodynamic simulations of low-mass protostellar outflows. We construct synthetic observations in 12CO in order to compare with observed outflows and evaluate the effects of beam resolution and outflow orientation on inferred outflow properties. To facilitate the comparison, we develop a quantitative prescription for measuring outflow opening angles. Using this prescription, we demonstrate that, in both simulations and synthetic observations, outflow opening angles broaden with time similarly to observed outflows. However, the interaction between the outflowing gas and the turbulent core envelope produces significant asymmetry between the red and blue shifted outflow lobes. We find that applying a velocity cutoff may result in outflow masses that are underestimated by a factor 5 or more, and masses derived from optically thick CO emission further underpredict the mass of the high-velocity gas by a factor of 5-10. Derived excitation temperatures indicate that outflowing gas is hotter than the ambient gas with temperature rising over time, which is in agreement with the simulation gas temperatures. However, excitation temperatures are otherwise not well correlated with the actual gas temperature.

The Formation of Low-Mass Binary Star Systems Via Turbulent Fragmentation

We characterize the infall rate onto protostellar systems forming in self-gravitating radiation-hydrodynamic simulations. Using two dimensionless parameters to determine disks susceptability to gravitational fragmentation, we infer limits on protostellar system multiplicity and the mechanism of binary formation. We show that these parameters give robust predictions even in the case of marginally resolved protostellar disks. We find that protostellar systems with radiation feedback predominately form binaries via turbulent fragmentation, not disk instability, and we predict turbulent fragmentation is the dominant channel for binary formation for low-mass stars. We clearly demonstrate that systems forming in simulations including radiative feedback have fundamentally different parameters than those in purely hydrodynamic simulations.

The Shapes of Molecular Cloud Cores in Simulations and Observations

In this study, we investigate the shapes of starless and protostellar cores using hydrodynamic, self-gravitating adaptive mesh refinement simulations of turbulent molecular clouds. We simulate observations of these cores in dust emission, including realistic noise and telescope resolution, and compare to the observed core shapes measured in Orion by Nutter & Ward-Thompson (2007). The simulations and the observations have generally high statistical similarity, with particularly good agreement between simulations and Orion B. Although protostellar cores tend to have semi-major axis to semi-minor axis ratios closer to one, the distribution of axis ratios for starless and protostellar cores are not significantly different for either the actual observations of Orion or the simulated observations. Because of the high level of agreement between the non-magnetic hydrodynamic simulations and observation, contrary to a number of previous authors, one cannot infer the presence of magnetic fields from core shape distributions.

The Kinematics of Molecular Cloud Cores in the Presence of Driven and Decaying Turbulence: Comparisons with Observations

In this study we investigate the formation and properties of prestellar and protostellar cores using hydrodynamic, self-gravitating Adaptive Mesh Refinement simulations, comparing the cases where turbulence is continually driven and where it is allowed to decay. We model observations of these cores in the C$^{18}$O$(2to 1)$, NH$_3(1,1)$, and N$_2$H$^+(1to 0)$ lines, and from the simulated observations we measure the linewidths of individual cores, the linewidths of the surrounding gas, and the motions of the cores relative to one another. Some of these distributions are significantly different in the driven and decaying runs, making them potential diagnostics for determining whether the turbulence in observed star-forming clouds is driven or decaying. Comparing our simulations with observed cores in the Perseus and $rho$ Ophiuchus clouds shows reasonably good agreement between the observed and simulated core-to-core velocity dispersions for both the driven and decaying cases. However, we find that the linewidths through protostellar cores in both simulations are too large compared to the observations. The disagreement is noticably worse for the decaying simulation, in which cores show highly supersonic infall signatures in their centers that decrease toward their edges, a pattern not seen in the observed regions. This result gives some support to the use of driven turbulence for modeling regions of star formation, but reaching a firm conclusion on the relative merits of driven or decaying turbulence will require more complete data on a larger sample of clouds as well as simulations that include magnetic fields, outflows, and thermal feedback from the protostars.