All molecular clouds are observed to be turbulent, but the origin, means of sustenance, and evolution of the turbulence remain debated. One possibility is that stellar feedback injects enough energy into the cloud to drive observed motions on parsec scales. Recent numerical studies of molecular clouds have found that feedback from stars, such as protostellar outflows and winds, injects energy and impacts turbulence. We expand upon these studies by analyzing magnetohydrodynamic simulations of winds interacting with molecular clouds which vary the stellar mass-loss rates and magnetic field strength. We generate synthetic ¹²CO(1-0) maps assuming that the simulations are at the distance of the nearby Perseus molecular cloud. By comparing the outputs from different initial conditions and evolutionary times, we identify differences in the synthetic observations and characterize these using common astrostatistics. We quantify the different statistical responses using a variety of metrics proposed in the literature. We find that multiple astrostatistics, such as principle component analysis, velocity component spectrum, and dendrograms, are sensitive to changes in stellar mass-loss rates and/or magnetic field strength. This demonstrates that stellar feedback influences molecular cloud turbulence and can be identified and quantified observationally using such statistics.

INTRODUCTION Young protostars are observed to launch energetic collimated bipolar mass outflows . These protostellar outflows play a fundamental role in the star formation process on a variety of scales. On sub-pc scales they entrain and unbind core gas, thus setting the efficiency at which dense gas turns into stars . Interaction between outflows and infalling material may regulate protostellar accretion and, ultimately, terminate it . On sub-pc up to cloud scales, outflows inject substantial energy into their surroundings, potentially providing a means of sustaining cloud turbulence over multiple dynamical times. The origin of outflows is attributed to the presence of magnetic fields, and a variety of different models have been proposed to explain the launching mechanism \citep[e.g.,][]{arce07}. Of these, the “disk-wind" model , in which the gas is centrifugally accelerated from the accretion disk surface, and the “X-wind" model , in which gas is accelerated along tightly wound field lines, are most commonly invoked to explain observed outflow signatures. However, investigating the launching mechanism is challenging because launching occurs on scales of a few stellar radii and during times when the protostar is heavily extincted by its natal gas. Consequently, separating outflow gas from accreting core gas, discriminating between models, and determining fundamental outflow properties are nontrivial. Three main approaches have been applied to studying outflows. First, single-dish molecular line observations have been successful in mapping the extent of outflows and their kinematics on core to cloud scales \citep[][]{bourke97,arce10,dunham14}. However, outflow gas with velocities comparable to the cloud turbulent velocity can only be extracted with additional assumptions and modeling \citep[e.g.,][]{arce01b,dunham14}, which are difficult to apply to confused, clustered star forming environments . Second, interferometry provide a means of mapping outflows down to 1,000 AU scales scales , and the Atacama Large Millimeter/submilllimeter Antenna (ALMA) is extending these limits down to sub-AU scales . However, interferometry is not suitable for producing large high-resolution maps and it resolves out larger scale structure. Consequently, it is difficult to assemble a complete and multi-scale picture of outflow properties with these observations. Finally, numerical simulations provide a complementary approach that supplies three-dimensional predictions for launching, entrainment and energy injection . The most promising avenue for understanding outflows lies at the intersection of numerical modeling and observations. By performing synthetic observations to model molecular and atomic lines, continuum, and observational effects, simulations can be mapped into the observational domain where they can be compared directly to observations \citep[e.g.,][]{Offner11,Offner12b,Mairs13}. Such direct comparisons are important for assessing the “reality" of the simulations, to interpret observational data and to assess observational uncertainties . In addition to observational instrument limitations, chemistry and radiative transfer introduce additional uncertainties that are difficult to quantify without realistic models . Synthetic observations have previously been performed in the context of understanding outflow opening angles , observed morphology , and impact on spectral energy distributions . The immanent completion of ALMA provides further motivation for predictive synthetic observations. Although ALMA will have unprecedented sensitivity and resolution compared to existing instruments, by nature interferometry resolves out large-scale structure and different configurations will be sensitive to different scales. Atmospheric noise and total observing time may also effect the fidelity of the data. Previous synthetic observations performed by suggest that the superior resolution of full ALMA and the Atacama Compact Array (ACA) will be able to resolve core structure and fragmentation prior to binary formation. predicts that ALMA will be able to resolve complex outflow velocity structure and helical structure in molecular emission. In this paper we seek to quantify the accuracy of different ALMA configurations in recovering fundamental gas properties such as mass, line-of-sight momentum, and energy. We use the casa software package to synthetically observe protostellar outflows in the radiation-hydrodynamic simulations of . By modeling the emission at different times, inclinations, molecular lines, and observing configurations we evaluate how well physical quantities can be measured in the star formation process. In section §[Methods] we describe our methods for modeling and observing outflows. In section §[results] we evaluate the effects of different observational parameters on bulk quantities. We discuss results and summarize conclusions in §[Conclusions].

ABSTRACT. Broadly, we seek to understand the role of pressure in star forming molecular clouds. We examine molecular line data of the Perseus region from the COMPLETE survey alongside radiative transfer-processed ‘observations’ of the turbulent simulations of S. Offner to try to (1) understand to what extent we can actually measure pressure through observations, and (2) study how pressure changes within a cloud’s substructure.