We study the formation, evolution and collapse of dense cores by tracking density structures in a magnetohydrodynamic (MHD) simulation. We identify cores using the dendrogram algorithm and utilize machine learning techniques, including principal component analysis (PCA) and the k-means clustering algorithm to analyze the full density and velocity dispersion profiles of these cores. We find that there exists an evolutionary sequence consisting of three distinct phases: i) the formation of turbulent density structures (Phase I), ii) the dissipation of turbulence and the formation of coherent cores (Phase II), and iii) the transition to protostellar cores through gravitational collapse (Phase III). In dynamically evolving molecular clouds, the existence of these three phases corresponds to the coexistence of three populations of cores with distinct physical properties. The prestellar and protostellar cores frequently analyzed in previous studies of observations and simulations belong to the last phase in this evolutionary picture. We derive typical lifetimes of 1.4$\pm$1.0$\times$10$^5$ yr, 3.3$\pm$1.4$\times$10$^5$ yr and 3.3$\pm$1.4$\times$10$^5$ yr, respectively for Phase I, II and III. We find that cores can form from both converging flows and filament fragmentation and that cores may form both inside and outside the filaments. We then compare our results to previous observations of coherent cores and provide suggestions for future observations to study cores belonging to the three phases.

In this chapter, we review some historical understanding and recent advances on the Initial Mass Function (IMF) and the Core Mass Function (CMF), both in terms of observations and theories. We focus mostly on star formation in clustered environment since this is suggested by observations to be the dominant mode of star formation. The statistical properties and the fragmentation behaviour of turbulent gas is discussed, and we also discuss the formation of binaries and small multiple systems.

Star formation is a multi-scale, multi-physics problem ranging from the size scale of molecular clouds ($\sim$10s pc) down to the size scales of dense prestellar cores ($\sim$0.1 pc) that are the birth sites of stars. Several physical processes like turbulence, magnetic fields and stellar feedback, such as radiation pressure and outflows, are more or less important for different stellar masses and size scales. During the last decade a variety of technological and computing advances have transformed our understanding of star formation through the use of multi-wavelength observations, large scale observational surveys, and multi-physics multi-dimensional numerical simulations. Additionally, the use of synthetic observations of simulations have provided a useful tool to interpret observational data and evaluate the importance of various physical processes on different scales in star formation. Here, we review these recent advancements in both high- ($M \gtrsim 8 \, M_{\rm \odot}$) and low-mass star formation.

Understanding the evolution of self-gravitating, isothermal, magnetized gas is crucial for star formation, as these physical processes have been postulated to set the initial mass function (IMF). We present a suite of isothermal magnetohydrodynamic (MHD) simulations using the GIZMO code, that resolve the formation of individual stars in giant molecular clouds (GMCs), spanning a range of Mach numbers found in observed GMCs. As in past works, the mean and median stellar masses are sensitive to numerical resolution, because they are sensitive to low-mass stars that contribute a vanishing fraction of the overall stellar mass. The {\em mass-weighted} median stellar mass $M_\mathrm{50}$ becomes insensitive to resolution once turbulent fragmentation is well-resolved. Without imposing Larson-like scaling laws, our simulations find $M_\mathrm{50} \propto M_\mathrm{0} \mathcal{M}^{-3} \alpha_\mathrm{turb} \mathrm{SFE}^{1/3}$ for GMC mass $M_\mathrm{0}$, sonic Mach number $\mathcal{M}$, virial parameter $\alpha_\mathrm{turb}$, and star formation efficiency $\mathrm{SFE}=M_\mathrm{\star}/M_\mathrm{0}$. This fit agrees well with previous IMF results from the RAMSES, ORION2, and SphNG codes. Although $M_\mathrm{50}$ has no significant dependence on the magnetic field strength at the cloud scale, MHD is necessary to prevent a fragmentation cascade that results in non-convergent stellar masses. For initial conditions and SFE similar to star-forming GMCs in our Galaxy, we predict $M_\mathrm{50}$ to be $>20 M_{\odot}$, an order of magnitude larger than observed ($\sim 2 M_\odot$), together with an excess of brown dwarfs. Moreover, $M_\mathrm{50}$ is sensitive to initial cloud properties and evolves strongly in time within a given cloud, predicting much larger IMF variations than are observationally allowed. We conclude that physics beyond MHD turbulence and gravity are necessary ingredients for the IMF.

Star clusters form in dense, hierarchically collapsing gas clouds. Bulk kinetic energy is transformed to turbulence with stars forming from cores fed by filaments. In the most compact regions, stellar feedback is least effective in removing the gas and stars may form very efficiently. These are also the regions where, in high-mass clusters, ejecta from some kind of high-mass stars are effectively captured during the formation phase of some of the low mass stars and effectively channeled into the latter to form multiple populations. Star formation epochs in star clusters are generally set by gas flows that determine the abundance of gas in the cluster. We argue that there is likely only one star formation epoch after which clusters remain essentially clear of gas by cluster winds. Collisional dynamics is important in this phase leading to core collapse, expansion and eventual dispersion of every cluster. We review recent developments in the field with a focus on theoretical work.

Star formation is a complex multi-scale phenomenon that is of significant importance for astrophysics in general. Stars and star formation are key pillars in observational astronomy from local star forming regions in the Milky Way up to high-redshift galaxies. From a theoretical perspective, star formation and feedback processes (radiation, winds, and supernovae) play a pivotal role in advancing our understanding of the physical processes at work, both individually and of their interactions. In this review we will give an overview of the main processes that are important for the understanding of star formation. We start with an observationally motivated view on star formation from a global perspective and outline the general paradigm of the life-cycle of molecular clouds, in which star formation is the key process to close the cycle. After that we focus on the thermal and chemical aspects in star forming regions, discuss turbulence and magnetic fields as well as gravitational forces. Finally, we review the most important stellar feedback mechanisms.

We explore the relation between the stellar mass surface density and the mass surface density of molecular hydrogen gas in twelve nearby molecular clouds that are located at $<$1.5 kpc distance. The sample clouds span an order of magnitude range in mass, size, and star formation rates. We use thermal dust emission from $Herschel$ maps to probe the gas surface density and the young stellar objects from the most recent $Spitzer$ Extended Solar Neighborhood Archive (SESNA) catalog to probe the stellar surface density. Using a star-sampled nearest neighbor technique to probe the star-gas surface density correlations at the scale of a few parsecs, we find that the stellar mass surface density varies as a power-law of the gas mass surface density, with a power-law index of $\sim$2 in all the clouds. The consistent power-law index implies that star formation efficiency is directly correlated with gas column density, and no gas column density threshold for star formation is observed. We compare the observed correlations with the predictions from an analytical model of thermal fragmentation, and with the synthetic observations of a recent hydrodynamic simulation of a turbulent star-forming molecular cloud. We find that the observed correlations are consistent for some clouds with the thermal fragmentation model and can be reproduced using the hydrodynamic simulations.

We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particularly efficient. Methanol and sulphur monoxide are most abundant in a narrow zone that follows the eastern side of the core. This side is sheltered from the stronger external radiation field coming from the west. We show that photodissociation on the illuminated side can give rise to an asymmetric methanol distribution, but that the stark contrast observed in H-MM1 is hard to explain without assuming enhanced desorption on the shaded side. The region of the brightest emission has a wavy structure that rolls up at one end. This is the signature of Kelvin-Helmholtz instability occurring in sheared flows. We suggest that in this zone, methanol and sulphur are released as a result of grain-grain collisions induced by shear vorticity.

We investigate the time evolution of dense cores identified in molecular cloud simulations using dendrograms, which are a common tool to identify hierarchical structure in simulations and observations of star formation. We develop an algorithm to link dendrogram structures through time using the three-dimensional density field from magnetohydrodynamical simulations, thus creating histories for all dense cores in the domain. We find that the population-wide distributions of core properties are relatively invariant in time, and quantities like the core mass function match with observations. Despite this consistency, an individual core may undergo large (>40%), stochastic variations due to the redefinition of the dendrogram structure between timesteps. This variation occurs independent of environment and stellar content. We identify a population of short-lived (<200 kyr) overdensities masquerading as dense cores that may comprise ~20% of any time snapshot. Finally, we note the importance of considering the full history of cores when interpreting the origin of the initial mass function; we find that, especially for systems containing multiple stars, the core mass defined by a dendrogram leaf in a snapshot is typically less than the final system stellar mass. This work reinforces that there is no time-stable density contour that defines a star-forming core. The dendrogram itself can induce significant structure variation between timesteps due to small changes in the density field. Thus, one must use caution when comparing dendrograms of regions with different ages or environment properties because differences in dendrogram structure may not come solely from the physical evolution of dense cores.

The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular clouds. We find that the results on relative alignment between cores and background magnetic field in synthetic observations slightly disagree with those measured in fully 3D simulation data, which is partly because cores identified in projected 2D maps tend to coexist within filamentary structures, while 3D cores are generally more rounded. In addition, we examine the progression of magnetic field from pc- to core-scale in the simulation, which is consistent with the anisotropic core formation model that gas preferably flow along the magnetic field toward dense cores. When comparing the observed cores identified from the GBT Ammonia Survey (GAS) and Planck polarization-inferred magnetic field orientations, we find that the relative core-field alignment has a regional dependence among different clouds. More specifically, we find that dense cores in the Taurus molecular cloud tend to align perpendicular to the background magnetic field, while those in Perseus and Ophiuchus tend to have random (Perseus) or slightly parallel (Ophiuchus) orientations with respect to the field. We argue that this feature of relative core-field orientation could be used to probe the relative significance of the magnetic field within the cloud.

Recent observations of global velocity gradients across and along molecular filaments have been interpreted as signs of gas accreting onto and along these filaments, potentially feeding star-forming cores and proto-clusters. The behavior of velocity gradients in filaments, however, has not been studied in detail, particularly on small scales (< 0.1 pc). In this paper, we present MUFASA, an efficient, robust, and automatic method to fit ammonia lines with multiple velocity components, generalizable to other molecular species. We also present CRISPy, a Python package to identify filament spines in 3D images (e.g., position-position-velocity cubes), along with a complementary technique to sort fitted velocity components into velocity-coherent filaments. In NGC 1333, we find a wealth of velocity gradient structures on a beam-resolved scale of ~0.05 pc. Interestingly, these local velocity gradients are not randomly oriented with respect to filament spines and their perpendicular, i.e., radial, component decreases in magnitude towards the spine for many filaments. Together with remarkably constant velocity gradients on larger scales along many filaments, these results suggest a scenario in which gas falling onto filaments is progressively damped and redirected to flow along these filaments.

In recent years, exciting developments have taken place in the identification of the role of cosmic rays in star-forming environments. Observations from radio to infrared wavelengths and theoretical modelling have shown that low-energy cosmic rays (<1 TeV) play a fundamental role in shaping the chemical richness of the interstellar medium, determining the dynamical evolution of molecular clouds. In this review we summarise in a coherent picture the main results obtained by observations and by theoretical models of propagation and generation of cosmic rays, from the smallest scales of protostars and circumstellar discs, to young stellar clusters, up to Galactic and extragalactic scales. We also discuss the new fields that will be explored in the near future thanks to new generation instruments, such as: CTA, for the $\gamma$-ray emission from high-mass protostars; SKA and precursors, for the synchrotron emission at different scales; and ELT/HIRES, JWST, and ARIEL, for the impact of cosmic rays on exoplanetary atmospheres and habitability.

Giant molecular clouds (GMCs) are well-studied in the local Universe, however, exactly how their properties vary during galaxy evolution is poorly understood due to challenging resolution requirements, both observational and computational. We present the first time-dependent analysis of giant molecular clouds in a Milky Way-like galaxy and an LMC-like dwarf galaxy of the FIRE-2 (Feedback In Realistic Environments) simulation suite, which have sufficient resolution to predict the bulk properties of GMCs in cosmological galaxy formation self-consistently. We show explicitly that the majority of star formation outside the galactic center occurs within self-gravitating gas structures that have properties consistent with observed bound GMCs. We find that the typical cloud bulk properties such as mass and surface density do not vary more than a factor of 2 in any systematic way after the first Gyr of cosmic evolution within a given galaxy from its progenitor. While the median properties are constant, the tails of the distributions can briefly undergo drastic changes, which can produce very massive and dense self-gravitating gas clouds. Once the galaxy forms, we identify only two systematic trends in bulk properties over cosmic time: a steady increase in metallicity produced by previous stellar populations and a weak decrease in bulk cloud temperatures. With the exception of metallicity we find no significant differences in cloud properties between the Milky Way-like and dwarf galaxies. These results have important implications for cosmological star and star cluster formation and put especially strong constraints on theories relating the stellar initial mass function to cloud properties.

We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and Flat Spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and Flat Spectrum protostars are 44.9$^{+5.8}_{-3.4}$, 37.0$^{+4.9}_{-3.0}$, and 28.5$^{+3.7}_{-2.3}$ au, respectively, and the mean protostellar dust disk masses are 25.9$^{+7.7}_{-4.0}$, 14.9$^{+3.8}_{-2.2}$, 11.6$^{+3.5}_{-1.9}$ Earth masses, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% out of 379 detected 0.87 mm continuum sources plus 42 non-detections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than that of Class II disks by a factor of $>$4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase.

We adopt the deep learning method CASI (Convolutional Approach to Shell Identification) and extend it to 3D (CASI-3D) to identify signatures of stellar feedback in molecular line spectra, such as 13CO. We adopt magneto-hydrodynamics simulations that study the impact of stellar winds in a turbulent molecular cloud as an input to generate synthetic observations. We apply the 3D radiation transfer code radmc-3d to model 13CO (J=1-0) line emission from the simulated clouds. We train two CASI-3d models: ME1 predicts only the position of feedback, while MF predicts the fraction of the mass coming from feedback in each voxel. We adopt 75% of the synthetic observations as the training set and assess the accuracy of the two models with the remaining data. We demonstrate that model ME1 identifies bubbles in simulated data with 95% accuracy, and model MF predicts the bubble mass within 4% of the true value. We use bubbles previously visually identified in Taurus in 13CO to validate the models and show both perform well on the highest confidence bubbles. We apply our two models on the full 98 deg2 FCRAO 13CO survey of the Taurus cloud. Models ME1 and MF predict feedback gas mass of 2894 M and 302 M, respectively. When including a correction factor for missing energy due to the limited velocity range of the 13CO data cube, model ME1 predicts feedback kinetic energies of 4.0*1e46 ergs and 1.5*1e47 ergs with/without subtracting the cloud velocity gradient. Model MF predicts feedback kinetic energy of 9.6*1e45 ergs and 2.8*1e46 ergs with/without subtracting the cloud velocity gradient. Model ME1 predicts bubble locations and properties consistent with previous visually identified bubbles. However, model MF demonstrates that feedback properties computed based on visual identifications are significantly overestimated due to line of sight confusion and contamination from background and foreground gas.

We present and release the full dataset for the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. This survey used the Submillimeter Array (SMA) to image the 74 known protostars within the Perseus molecular cloud. The SMA was used in two array configurations to capture outflows for scales $>$30$^{\prime\prime}$ ($>$9000 au) and to probe scales down to $\sim$1$^{\prime\prime}$ ($\sim$300 au). The protostars were observed with the 1.3 mm and 850 $\mu$m receivers simultaneously to detect continuum at both wavelengths and molecular line emission from CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1), N$_2$D$^+$(3-2), CO(3-2), HCO$^+$(4-3), and H$^{13}$CO$^+$(4-3). Some of the observations also used the SMA's recently upgraded correlator, SWARM, whose broader bandwidth allowed for several more spectral lines to be observed (e.g., SO, H$_2$CO, DCO$^+$, DCN, CS, CN). Of the main continuum and spectral tracers observed, 84% of the images and cubes had emission detected. The median C$^{18}$O(2-1) linewidth is $\sim$1.0 km s$^{-1}$, which is slightly higher than those measured with single-dish telescopes at scales of 3000-20000 au. Of the 74 targets, six are suggested to be first hydrostatic core candidates, and we suggest that L1451-mm is the best candidate. We question a previous continuum detection toward L1448 IRS2E. In the SVS13 system, SVS13A certainly appears to be the most evolved source, while SVS13C appears to be hotter and more evolved than SVS13B. The MASSES survey is the largest publicly available interferometric continuum and spectral line protostellar survey to date, and is largely unbiased as it only targets protostars in Perseus. All visibility ($uv$) data and imaged data are publicly available at https://dataverse.harvard.edu/dataverse/full_MASSES/.

Stars rarely form in isolation. Nearly half of the stars in the Milky Way have a companion, and this fraction increases in star-forming regions. However, why some dense cores and filaments form bound pairs while others form single stars remains unclear. We present a set of three-dimensional, gravo-magnetohydrodynamic simulations of turbulent star-forming clouds, aimed at understanding the formation and evolution of multiple-star systems formed through large scale (>~$10^3$ AU) turbulent fragmentation. We investigate three global magnetic field strengths, with global mass-to-flux ratios of $\mu_\phi$=2, 8, and 32. The initial separations of protostars in multiples depends on the global magnetic field strength, with stronger magnetic fields (e.g., $\mu_\phi$=2) suppressing fragmentation on smaller scales. The overall multiplicity fraction (MF) is between 0.4-0.6 for our strong and intermediate magnetic field strengths, which is in agreement with observations. The weak field case has a lower fraction. The MF is relatively constant throughout the simulations, even though stellar densities increase as collapse continues. While the MF rarely exceeds 60% in all three simulations, over 80% of all protostars are part of a binary system at some point. We additionally find that the distribution of binary spin mis-alignment angles is consistent with a randomized distribution. In all three simulations, several binaries originate with wide separations and dynamically evolve to <~ $10^2$ AU separations. We show that a simple model of mass accretion and dynamical friction with the gas can explain this orbital evolution.

We present an analysis of the internal velocity structures of the newly identified sub-0.1 pc coherent structures, droplets, in L1688 and B18. By fitting 2D linear velocity fields to the observed maps of velocity centroids, we determine the magnitudes of linear velocity gradients and examine the potential rotational motions that could lead to the observed velocity gradients. The results show that the droplets follow the same power-law relation between the velocity gradient and size found for larger-scale dense cores. Assuming that rotational motion giving rise to the observed velocity gradient in each core is a solid-body rotation of a rotating body with a uniform density, we derive the "net rotational motions" of the droplets. We find a ratio between rotational and gravitational energies, $\beta$, of $\sim 0.046$ for the droplets, and when including both droplets and larger-scale dense cores, we find $\beta \sim 0.039$. We then examine the alignment between the velocity gradient and the major axis of each droplet, using methods adapted from the histogram of relative orientations (HRO) introduced by Soler et al. (2013). We find no definitive correlation between the directions of velocity gradients and the elongations of the cores. Lastly, we discuss physical processes other than rotation that may give rise to the observed velocity field.

We present ALMA (0.87~mm) and VLA (9~mm) observations toward OMC2-FIR4 and OMC2-FIR3 within the Orion integral-shaped filament that are thought to be the nearest regions of intermediate mass star formation. We characterize the continuum sources within these regions on $\sim$40~AU (0\farcs1) scales and associated molecular line emission at a factor of $\sim$30 better resolution than previous observations at similar wavelengths. We identify six compact continuum sources within OMC2-FIR4, four in OMC2-FIR3, and one additional source just outside OMC2-FIR4. This continuum emission is tracing the inner envelope and/or disk emission on less than 100~AU scales. HOPS-108 is the only protostar in OMC2-FIR4 that exhibits emission from high-excitation transitions of complex organic molecules (e.g., methanol and other lines) coincident with the continuum emission. HOPS-370 in OMC2-FIR3 with L~$\sim$~360~\lsun, also exhibits emission from high-excitation methanol and other lines. The methanol emission toward these two protostars is indicative of temperatures high enough to thermally evaporate methanol from icy dust grains; overall these protostars have characteristics similar to hot corinos. We do not identify a clear outflow from HOPS-108 in \twco, but find evidence of interaction between the outflow/jet from HOPS-370 and the OMC2-FIR4 region. The multitude of observational constraints indicate that HOPS-108 is likely a low to intermediate-mass protostar in its main mass accretion phase and it is the most luminous protostar in OMC2-FIR4. The high resolution data presented here are essential for disentangling the embedded protostars from their surrounding dusty environments and characterizing them.

We present the first results of high-spectral resolution (0.023 km/s) N$_2$H$^+$ observations of dense gas dynamics at core scales (~0.01 pc) using the recently commissioned Argus instrument on the Green Bank Telescope (GBT). While the fitted linear velocity gradients across the cores measured in our targets nicely agree with the well-known power-law correlation between the specific angular momentum and core size, it is unclear if the observed gradients represent core-scale rotation. In addition, our Argus data reveal detailed and intriguing gas structures in position-velocity (PV) space for all 5 targets studied in this project, which could suggest that the velocity gradients previously observed in many dense cores actually originate from large-scale turbulence or convergent flow compression instead of rigid-body rotation. We also note that there are targets in this study with their star-forming disks nearly perpendicular to the local velocity gradients, which, assuming the velocity gradient represents the direction of rotation, is opposite to what is described by the classical theory of star formation. This provides important insight on the transport of angular momentum within star-forming cores, which is a critical topic on studying protostellar disk formation.

We use the gas-grain chemistry code UCLCHEM to explore the impact of cosmic-ray feedback on the chemistry of circumstellar disks. We model the attenuation and energy losses of the cosmic-rays as they propagate outwards from the star and also consider ionization due to stellar radiation and radionuclides. For accretion rates typical of young stars, $\dot M_* \sim 10^{-9}-10^{-6}$ M_\odot yr$^{-1}$, we show that cosmic rays accelerated by the stellar accretion shock produce a cosmic-ray ionization rate at the disk surface $\zeta \gtrsim 10^{-15}$ s$^{-1}$, at least an order of magnitude higher than the ionization rate associated with the Galactic cosmic-ray background. The incident cosmic-ray flux enhances the disk ionization at intermediate to high surface densities ($\Sigma > 10$ g cm$^{-2}$) particularly within 10 au of the star. We find the dominant ions are C$^+$, S$^+$ and Mg$^+$ in the disk surface layers, while the H$_3^+$ ion dominates at surface densities above 1.0 g cm$^{-2}$. We predict the radii and column densities at which the magneto-rotational instability (MRI) is active in T Tauri disks and show that ionization by cosmic-ray feedback extends the MRI-active region towards the disk mid-plane. However, the MRI is only active at the mid-plane of a minimum mass solar nebula disk if cosmic-rays propagate diffusively ($\zeta \propto r^{-1}$) away from the star. The relationship between accretion, which accelerates cosmic rays, the dense accretion columns, which attenuate cosmic rays, and the MRI, which facilitates accretion, create a cosmic-ray feedback loop that mediates accretion and may produce luminosity variability.

We present initial results from the K-band focal plane array Examinations of Young STellar Object Natal Environments (KEYSTONE) survey, a large project on the 100-m Green Bank Telescope mapping ammonia emission across eleven giant molecular clouds at distances of $0.9-3.0$ kpc (Cygnus X North, Cygnus X South, M16, M17, MonR1, MonR2, NGC2264, NGC7538, Rosette, W3, and W48). This data release includes the NH$_3$ (1,1) and (2,2) maps for each cloud, which are modeled to produce maps of kinetic temperature, centroid velocity, velocity dispersion, and ammonia column density. Median cloud kinetic temperatures range from $11.4\pm2.2$ K in the coldest cloud (MonR1) to $23.0\pm6.5$ K in the warmest cloud (M17). Using dendrograms on the NH$_3$ (1,1) integrated intensity maps, we identify 856 dense gas clumps across the eleven clouds. Depending on the cloud observed, $40-100\%$ of the clumps are aligned spatially with filaments identified in H$_2$ column density maps derived from SED-fitting of dust continuum emission. A virial analysis reveals that 523 of the 835 clumps ($\sim63\%$) with mass estimates are bound by gravity alone. We find no significant difference between the virial parameter distributions for clumps aligned with the dust-continuum filaments and those unaligned with filaments. In some clouds, however, hubs or ridges of dense gas with unusually high mass and low virial parameters are located within a single filament or at the intersection of multiple filaments. These hubs and ridges tend to host water maser emission, multiple 70$\mu$m-detected protostars, and have masses and radii above an empirical threshold for forming massive stars.

We utilize a modified astrochemistry code which includes cosmic ray attenuation in-situ to quantify the impact of different cosmic-ray models on the CO-to-H$_2$ and CI-to-H$_2$ conversion factors, $X_{\rm CO}$ and $X_{\rm CI}$, respectively. We consider the impact of cosmic rays accelerated by accretion shocks, and show that clouds with star formation efficiencies greater than 2\% have $X_{\rm CO} = (2.5 \pm 1)\times10^{20}$ cm$^{-2}$(K km s$^{-1}$)$^{-1}$, consistent with Milky Way observations. We find that changing the cosmic ray ionization rate from external sources from the canonical $\zeta \approx 10^{-17}$ to $\zeta \approx 10^{-16}$ s$^{-1}$, which better represents observations in diffuse gas, reduces $X_{\rm CO}$ by 0.2 dex for clusters with surface densities below 3 g cm$^{-2}$. We show that embedded sources regulate $X_{\rm CO}$ and decrease its variance across a wide range of surface densities and star formation efficiencies. Our models reproduce the trends of a decreased $X_{\rm CO}$ in extreme cosmic ray environments. $X_{\rm CI}$ has been proposed as an alternative to $X_{\rm CO}$ due to its brightness at high redshifts. The inclusion of internal cosmic ray sources leads to 1.2 dex dispersion in $X_{\rm CI}$ ranging from $2\times10^{20} < X_{\rm CI} < 4\times10^{21}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$. We show that $X_{\rm CI}$ is highly sensitive to the underlying cosmic ray model.

During star formation, the accretion disk drives fast MHD winds which usually contain two components, a collimated jet and a radially distributed wide-angle wind. These winds entrain the surrounding ambient gas producing molecular outflows. We report recent observation of 12CO (2-1) emission of the HH 46/47 molecular outflow by the Atacama Large Millimeter/sub-millimeter Array, in which we identify multiple wide-angle outflowing shell structures in both the blue and red-shifted outflow lobes. These shells are highly coherent in position-position-velocity space, extending to >40-50 km/s in velocity and 10^4 au in space with well defined morphology and kinematics. We suggest these outflowing shells are the result of the entrainment of ambient gas by a series of outbursts from an intermittent wide-angle wind. Episodic outbursts in collimated jets are commonly observed, yet detection of a similar behavior in wide-angle winds has been elusive. Here we show clear evidence that the wide-angle component of the HH 46/47 protostellar outflows experiences similar variability seen in the collimated component.

We utilize techniques from deep learning to identify signatures of stellar feedback in simulated molecular clouds. Specifically, we implement a deep neural network with an architecture similar to U-Net and apply it to the problem of identifying wind-driven shells and bubbles using data from magneto-hydrodynamic simulations of turbulent molecular clouds with embedded stellar sources. The network is applied to two tasks, dense regression and segmentation, on two varieties of data, simulated density and synthetic 12 CO observations. Our Convolutional Approach for Shell Identification (CASI) is able to obtain a true positive rate greater than 90\%, while maintaining a false positive rate of 1\%, on two segmentation tasks and also performs well on related regression tasks. The source code for CASI is available on GitLab.

Citizen science has helped astronomers comb through large data sets to identify patterns and objects that are not easily found through automated processes. The Milky Way Project (MWP), a citizen science initiative on the Zooniverse platform, presents internet users with infrared (IR) images from Spitzer Space Telescope Galactic plane surveys. MWP volunteers make classification drawings on the images to identify targeted classes of astronomical objects. We present the MWP second data release (DR2) and an updated data reduction pipeline written in Python. We aggregate ${\sim}3$ million classifications made by MWP volunteers during the years 2012-2017 to produce the DR2 catalogue, which contains 2600 IR bubbles and 599 candidate bow-shock driving stars. The reliability of bubble identifications, as assessed by comparison to visual identifications by trained experts and scoring by a machine-learning algorithm, is found to be a significant improvement over DR1. We assess the reliability of IR bow shocks via comparison to expert identifications and the colours of candidate bow-shock driving stars in the 2MASS point-source catalogue. We hence identify highly-reliable subsets of 1394 DR2 bubbles and 453 bow-shock driving stars. Uncertainties on object coordinates and bubble size/shape parameters are included in the DR2 catalog. Compared with DR1, the DR2 bubbles catalogue provides more accurate shapes and sizes. The DR2 catalogue identifies 311 new bow shock driving star candidates, including three associated with the giant HII regions NGC 3603 and RCW 49.

We present the observation and analysis of newly discovered coherent structures in the L1688 region of Ophiuchus and the B18 region of Taurus. Using data from the Green Bank Ammonia Survey (GAS), we identify regions of high density and near-constant, almost-thermal, velocity dispersion. Eighteen coherent structures are revealed, twelve in L1688 and six in B18, each of which shows a sharp "transition to coherence" in velocity dispersion around its periphery. The identification of these structures provides a chance to study the coherent structures in molecular clouds statistically. The identified coherent structures have a typical radius of 0.04 pc and a typical mass of 0.4 Msun, generally smaller than previously known coherent cores identified by Goodman et al. (1998), Caselli et al. (2002), and Pineda et al. (2010). We call these structures "droplets." We find that unlike previously known coherent cores, these structures are not virially bound by self-gravity and are instead predominantly confined by ambient pressure. The droplets have density profiles shallower than a critical Bonnor-Ebert sphere, and they have a velocity (VLSR) distribution consistent with the dense gas motions traced by NH3 emission. These results point to a potential formation mechanism through pressure compression and turbulent processes in the dense gas. We present a comparison with a magnetohydrodynamic simulation of a star-forming region, and we speculate on the relationship of droplets with larger, gravitationally bound coherent cores, as well as on the role that droplets and other coherent structures play in the star formation process.

We present astrochemical photo-dissociation region models in which cosmic ray attenuation has been fully coupled to the chemical evolution of the gas. We model the astrochemical impact of cosmic rays, including those accelerated by protostellar accretion shocks, on molecular clouds hosting protoclusters. Our models with embedded protostars reproduce observed ionization rates. We study the imprint of cosmic ray attenuation on ions for models with different surface cosmic ray spectra and different star formation efficiencies. We find that abundances, particularly ions, are sensitive to the treatment of cosmic rays. We show the column densities of ions are under predicted by the `classic' treatment of cosmic rays by an order of magnitude. We also test two common chemistry approximations used to infer ionization rates. We conclude that the approximation based on the H$_3^+$ abundance under predicts the ionization rate except in regions where the cosmic rays dominate the chemistry. Our models suggest the chemistry in dense gas will be significantly impacted by the increased ionization rates, leading to a reduction in molecules such as NH$_3$ and causing H$_2$-rich gas to become [C II] bright.

We present TurbuStat (v1.0): a Python package for computing turbulence statistics in spectral-line data cubes. TurbuStat includes implementations of fourteen methods for recovering turbulent properties from observational data. Additional features of the software include: distance metrics for comparing two data sets; a segmented linear model for fitting lines with a break-point; a two-dimensional elliptical power-law model; multi-core fast-fourier-transform support; a suite for producing simulated observations of fractional Brownian Motion fields, including two-dimensional images and optically-thin HI data cubes; and functions for creating realistic world coordinate system information for synthetic observations. This paper summarizes the TurbuStat package and provides representative examples using several different methods. TurbuStat is an open-source package and we welcome community feedback and contributions.

Significant advances have been made over the past decade in the characterization of multiple protostar systems, enabled by the Karl G. Jansky Very Large Array (VLA), high-resolution infrared observations with the Hubble Space Telescope, and ground-based facilities. To further understand the mechanism(s) of multiple star formation, a combination of statistics, high-angular resolution radio/millimeter continuum imaging, characterization of kinematic structure, magnetic fields via polarimetry, and comparison with numerical simulations are needed. Thus, understanding the origin of stellar multiplicity in different regimes of companion separation will soon be within reach. However, to overcome challenges that studies in this field are now confronted with, a range of new capabilities are required: a new millimeter/centimeter wave facility with 10 mas resolution at {\lambda}=1 cm, space-based near to far-infrared observatories, continued development of low to high resolution spectroscopy on 3m to 10m class telescopes, and an ELT-class telescope with near to mid-infrared imaging/spectroscopic capability.

Knowledge of protostellar evolution has been revolutionized with the advent of surveys at near-infrared to submillimeter wavelengths. This has enabled the bolometric luminosities and bolometric temperatures (traditional protostellar evolution diagnostics) to be measured for large numbers of protostars. However, further progress is difficult without knowing the masses of the central protostars. Protostar masses can be most accurately determined via molecular line kinematics from millimeter interferometers (i.e., ALMA). Theoretical investigations have predicted the protostellar mass function (PMF) for various protostellar mass accretion models, and it is now imperative to observationally constrain its functional form. While ALMA has enabled protostellar mass measurements, samples approaching 100 sources are necessary to constrain the functional form of the PMF, and upgrades to ALMA and/or a new mm/cm facility will increase the feasibility of measuring such a large number of protostar masses. The masses of protostars will enable their stellar structure (radius and intrinsic luminosity), evolution, and accretion histories to be better understood. This is made more robust when effective temperatures and accretion rates can be measured via ground/space-based near to mid-infrared spectroscopy. Furthermore, access to supercomputing facilities is essential to fit the protostar masses via radiative transfer modeling and updated theoretical/numerical modeling of stellar structure may also be required.

Background: low-mass stars are the dominant product of the star formation process, and they trace star formation over the full range of environments, from isolated globules to clusters in the central molecular zone. In the past two decades, our understanding of the spatial distribution and properties of young low-mass stars and protostars has been revolutionized by sensitive space-based observations at X-ray and IR wavelengths. By surveying spatial scales from clusters to molecular clouds, these data provide robust measurements of key star formation properties. Goal: with their large numbers and their presence in diverse environments, censuses of low mass stars and protostars can be used to measure the dependence of star formation on environmental properties, such as the density and temperature of the natal gas, strengths of the magnetic and radiation fields, and the density of stars. Here we summarize how such censuses can answer three basic questions: i.) how is the star formation rate influenced by environment, ii.) does the IMF vary with environment, and iii.) how does the environment shape the formation of bound clusters? Answering these questions is an important step toward understanding star and cluster formation across the extreme range of environments found in the Universe. Requirements: sensitivity and angular resolution improvements will allow us to study the full range of environments found in the Milky Way. High spatial dynamic range (< 1arcsec to > 1degree scales) imaging with space-based telescopes at X-ray, mid-IR, and far-IR and ground-based facilities at near-IR and sub-mm wavelengths are needed to identify and characterize young stars.

The majority of the ultimate main-sequence mass of a star is assembled in the protostellar phase, where a forming star is embedded in an infalling envelope and encircled by a protoplanetary disk. Studying mass accretion in protostars is thus a key to understanding how stars gain their mass and ultimately how their disks and planets form and evolve. At this early stage, the dense envelope reprocesses most of the luminosity generated by accretion to far-infrared and submillimeter wavelengths. Time-domain photometry at these wavelengths is needed to probe the physics of accretion onto protostars, but variability studies have so far been limited, in large part because of the difficulty in accessing these wavelengths from the ground. We discuss the scientific progress that would be enabled with far-infrared and submillimeter programs to probe protostellar variability in the nearest kiloparsec.

Understanding the collapse of clouds and the formation of protoplanetary disks is essential to understanding the formation of stars and planets. Infall and accretion, the mass-aggregation processes that occur at envelope and disk scales, drive the dynamical evolution of protostars. While the observations of protostars at different stages constrain their evolutionary tracks, the impact of variability due to accretion bursts on dynamical and chemical evolution of the source is largely unknown. The lasting effects on protostellar envelopes and disks are tracked through multi-wavelength and time domain observational campaigns, requiring deep X-ray, infrared, and radio imaging and spectroscopy, at a sufficient level of spatial detail to distinguish contributions from the various substructures (i.e., envelope from disk from star from outflow). Protostellar models derived from these campaigns will illuminate the initial chemical state of protoplanetary disks during the epoch of giant planet formation. Insight from individual star formation in the Milky Way is also necessary to understand star formation rates in extragalactic sources. This cannot be achieved with ground-based observatories and is not covered by currently approved instrumentation. Requirements: High (v < 10 km/s for survey; v < 1 km/s for followup) spectral resolution capabilities with relatively rapid response times in the IR (3-500 um), X-ray (0.1-10 keV), and radio (cm) are critical to follow the course of accretion and outflow during an outburst. Complementary, AU-scale radio observations are needed to probe the disk accretion zone, and 10 AU-scale to probe chemical and kinematic structures of the disk-forming regions, and track changes in the dust, ice, and gas within protostellar envelopes.

We perform a virial analysis of starless dense cores in three nearby star-forming regions : L1688 in Ophiuchus, NGC 1333 in Perseus, and B18 in Taurus. Our analysis takes advantage of comprehensive kinematic information for the dense gas in all of these regions made publicly available through the Green Bank Ammonia Survey Data Release 1, which used to estimate internal support against collapse. We combine this information with ancillary data used to estimate other important properties of the cores, including continuum data from the James Clerk Maxwell Telescope Gould Belt Survey for core identification, core masses, and core sizes. Additionally, we used \textit{Planck} and \textit{Herschel}-based column density maps for external cloud weight pressure, and Five College Radio Astronomy Observatory $^{13}$CO observations for external turbulent pressure. Our self-consistent analysis suggests that many dense cores in all three star-forming regions are not bound by gravity alone, but rather require additional pressure confinement to remain bound. Unlike a recent, similar study in Orion~A, we find that turbulent pressure represents a significant portion of the external pressure budget. Our broad conclusion emphasizing the importance of pressure confinement in dense core evolution, however, agrees with earlier work.

Context: The study of dust emission at millimeter wavelengths is important to shed light on the dust properties and physical structure of pre-stellar cores, the initial conditions in the process of star and planet formation. Aims: Using two new continuum facilities, AzTEC at the LMT and MUSTANG-2 at the GBO, we aim to detect changes in the optical properties of dust grains as a function of radius for the well-known pre-stellar core L1544. Methods: We determine the emission profiles at 1.1 and 3.3 mm and examine whether they can be reproduced in terms of the current best physical models for L1544. We also make use of various tools to determine the radial distributions of the density, temperature, and the dust opacity in a self-consistent manner. Results: We find that our observations cannot be reproduced without invoking opacity variations. With the new data, new temperature and density profiles, as well as opacity variations across the core, have been derived. The opacity changes are consistent with the expected variations between uncoagulated bare grains, toward the outer regions of the core, and grains with thick ice mantles, toward the core center. A simple analytical grain growth model predicts the presence of grains of ~3-4 um within the central 2000 au for the new density profile.

We present combined observations of the NH3 (J,K) = (1,1) and (2,2) inversion transitions towards OMC1 in Orion A obtained by the Karl G. Jansky Very Large Array (VLA) and the 100 m Robert C. Byrd Green Bank Telescope (GBT). With an angular resolution of 6" (0.01 pc), these observations reveal with unprecedented detail the complex filamentary structure extending north of the active Orion BN/KL region in a field covering 6' x 7'. We find a 0.012 pc wide filament within OMC1, with an aspect ratio of ~37:1, that was missed in previous studies. Its orientation is directly compared to the relative orientation of the magnetic field from the James Clerk Maxwell Telescope BISTRO survey in Orion A. We find a small deviation of ~11 deg between the mean orientation of the filament and the magnetic field, suggesting that they are almost parallel to one another. The filament's column density is estimated to be 2-3 orders of magnitude larger than the filaments studied with Herschel and is possibly self-gravitating given the low values of turbulence found. We further produce maps of the gas kinematics by forward modeling the hyperfine structure of the NH3 (J,K) = (1,1) and (2,2) lines. The resulting distribution of velocity dispersions peaks at ~0.5 km/s, close to the subsonic regime of the gas. This value is about 0.2 km/s smaller than previously measured in single-dish observations of the same region, suggesting that higher angular and spectral resolution observations will identify even lower velocity dispersions that might reach the subsonic turbulence regime in dense gas filaments.

We present a 1.1~mm census of dense cores in the Mon~R2 Giant Molecular Cloud with the AzTEC instrument on the Large Millimeter Telescope (LMT). We detect 295 cores (209 starless, and 86 with protostars) in a two square degree shallow survey. We also carry out a deep follow-up survey of 9 regions with low to intermediate ($3<A_V<7$) gas column densities and detect 60 new cores in the deeper survey which allows us to derive a completeness limit. After performing corrections for low signal-to-noise cores, we find a median core mass of $\sim 2.1 \text{M}_{\odot}$ and a median size of $ 0.08$~pc. $46\%$ of the cores (141) have masses exceeding the local Bonor-Ebert mass for cores with T=12K, suggesting that in the absence of supporting non-thermal pressure, these regions are unstable to gravitational collapse. We present the core mass function (CMF) for various subdivisions of the core sample. We find that cores with masses $>$10~$M_{\odot}$ are exclusively found in regions with high core number densities and that the CMF of the starless cores has an excess of low-mass cores ($<$5~$M_{\odot}$) compared to the CMF of protostellar cores. We report a power law correlation of index $1.99 \pm 0.03$ between local core mass density and gas column density (as traced by Herschel) over a wide range of size scales (0.3-5~pc). This power law is consistent with that predicted for thermal fragmentation of a self-gravitating sheet.

A fundamental property of molecular clouds is that they are turbulent, but how this turbulence is generated and maintained is unknown. One possibility is that stars forming within the cloud regenerate turbulence via their outflows, winds and radiation ("feedback"). Disentangling motions created by feedback from the initial cloud turbulence is challenging, however. Here we confront the relationship between stellar feedback and turbulence by identifying and separating the local and global impact of stellar winds. We analyze magnetohydrodyanamic simulations in which we track wind material as it interacts with the ambient cloud. By distinguishing between launched material, gas entrained by the wind and pristine gas we show energy is transferred away from the sources via magnetic waves excited by the expanding wind shells. This action at a distance enhances the fraction of stirring motion compared to compressing motion and produces a flatter velocity power spectrum. We conclude stellar feedback accounts for significant energy transfer within molecular clouds, an impact enhanced by magnetic waves, which have previously been neglected by observations. Altogether, stellar feedback can partially offset global turbulence dissipation.

We present the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey, which uses the Submillimeter Array (SMA) interferometer to map the continuum and molecular lines for all 74 known Class 0/I protostellar systems in the Perseus molecular cloud. The primary goal of the survey is to observe an unbiased sample of young protostars in a single molecular cloud so that we can characterize the evolution of protostars. This paper releases the MASSES 1.3 mm data from the subcompact configuration ($\sim$4$^{\prime\prime}$ or $\sim$1000 au resolution), which is the SMA's most compact array configuration. We release both $uv$ visibility data and imaged data for the spectral lines CO(2-1),$^{13}$CO(2-1), C$^{18}$O(2-1), and N$_2$D$^+$(3-2), as well as for the 1.3 mm continuum. We identify the tracers that are detected toward each source. We also show example images of continuum and CO(2-1) outflows, analyze C$^{18}$O(2-1) spectra, and present data from the SVS 13 star-forming region. The calculated envelope masses from the continuum show a decreasing trend with bolometric temperature (a proxy for age). Typical C$^{18}$O(2-1) linewidths are 1.45 km s$^{-1}$, which is higher than the C$^{18}$O linewidths detected toward Perseus filaments and cores. We find that N$_2$D$^+$(3-2) is significantly more likely to be detected toward younger protostars. We show that the protostars in SVS 13 are contained within filamentary structures as traced by C$^{18}$O(2-1) and N$_2$D$^+$(3-2). We also present the locations of SVS 13A's high velocity (absolute line-of-sight velocities $>$150 km s$^{-1}$) red and blue outflow components. Data can be downloaded from https://dataverse.harvard.edu/dataverse/MASSES .

We construct a model for cosmic ray acceleration from protostellar accretion shocks and calculate the resulting cosmic ray ionization rate within star-forming molecular clouds. We couple a protostar cluster model with an analytic accretion shock model to calculate the cosmic ray acceleration from protostellar surfaces. We present the cosmic ray flux spectrum from keV to GeV energies for a typical low-mass protostar. We find that at the shock surface the spectrum follows a power-law trend across 6 orders of magnitude in energy. After attenuation, the spectrum at high energies steepens, while at low energies it is relatively flat. We calculate the cosmic ray pressure and cosmic ray ionization rate from relativistic protons at the protostellar surface and at the edge of the core. We present the cosmic ray ionization rate for individual protostars as a function of their instantaneous mass and final mass. The protostellar cosmic ray ionization rate is $\zeta \approx 0.01 - 1$ s$^{-1}$ at the accretion shock surface. However, at the edge of the core, the cosmic ray ionization rate drops substantially to between $\zeta \approx 10^{-20}$ to $10^{-17}$ s$^{-1}$. There is a large spatial gradient in the cosmic ray ionization rate, such that inner regions may experience cosmic ray ionization rates larger than the often assumed fiducial rate, $\zeta = 3\times10^{-17}$ s$^{-1}$. Finally, we calculate the cosmic ray ionization rate for protostellar clusters over 5 orders of magnitude of cluster size. We find that clusters with more than approximately 200 protostars produce a higher cosmic ray ionization rate within their natal cloud than the fiducial galactic value.

We analyze hydrodynamic simulations of turbulent, star-forming molecular clouds that are post-processed with the photo-dissociation region astrochemistry code 3D-PDR. We investigate the sensitivity of 15 commonly applied turbulence statistics to post-processing assumptions, namely variations in gas temperature, abundance and external radiation field. We produce synthetic $^{12}$CO(1-0) and CI($^{3}$P$_{1}$-$^{3}$P$_{0}$) observations and examine how the variations influence the resulting emission distributions. To characterize differences between the datasets, we perform statistical measurements, identify diagnostics sensitive to our chemistry parameters, and quantify the statistic responses by using a variety of distance metrics. We find that multiple turbulent statistics are sensitive not only to the chemical complexity but also to the strength of the background radiation field. The statistics with meaningful responses include principal component analysis, spatial power spectrum and bicoherence. A few of the statistics, such as the velocity coordinate spectrum, are primarily sensitive to the type of tracer being utilized, while others, like the delta-variance, strongly respond to the background radiation field. Collectively, these findings indicate that more realistic chemistry impacts the responses of turbulent statistics and is necessary for accurate statistical comparisons between models and observed molecular clouds.

We construct a semi-analytic model to study the effect of far-ultraviolet (FUV) radiation on gas chemistry from embedded protostars. We use the Protostellar Luminosity Function (PLF) formalism of Offner & McKee (2011) to calculate the total, FUV, and ionizing cluster luminosity for various protostellar accretion histories and cluster sizes. We compare the model predictions with surveys of Gould Belt star-forming regions and find the Tapered Turbulent Core model matches best the mean luminosities and the spread in the data. We combine the cluster model with the photo-dissociation region astrochemistry code, {\sc 3d-pdr}, to compute the impact of the FUV luminosity from embedded protostars on the CO to H$_2$ conversion factor, $X_{\rm CO}$, as a function of cluster size, gas mass and star formation efficiency. We find that $X_{\rm CO}$ has a weak dependence on the FUV radiation from embedded sources for large clusters due to high cloud optical depths. In smaller and more efficient clusters the embedded FUV increases X$_{\rm CO}$ to levels consistent with the average Milky Way values. The internal physical and chemical structure of the cloud are significantly altered, and $X_{\rm CO}$ depends strongly on the protostellar cluster mass for small efficient clouds.

We present a study of hierarchical structure in the Perseus molecular cloud, from the scale of the entire cloud ($\gtrsim$10 pc) to smaller clumps ($\sim$1 pc), cores ($\sim$0.05-0.1 pc), envelopes ($\sim$300-3000 AU) and protostellar objects ($\sim$15 AU). We use new observations from the Submillimeter Array (SMA) large project "Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES)" to probe the envelopes, and recent single-dish and interferometric observations from the literature for the remaining scales. This is the first study to analyze hierarchical structure over five scales in the same cloud complex. We compare the number of fragments with the number of Jeans masses in each scale to calculate the Jeans efficiency, or the ratio of observed to expected number of fragments. The velocity dispersion is assumed to arise either from purely thermal motions, or from combined thermal and non-thermal motions inferred from observed spectral line widths. For each scale, thermal Jeans fragmentation predicts more fragments than observed, corresponding to inefficient thermal Jeans fragmentation. For the smallest scale, thermal plus non-thermal Jeans fragmentation also predicts too many protostellar objects. However at each of the larger scales thermal plus non-thermal Jeans fragmentation predicts fewer than one fragment, corresponding to no fragmentation into envelopes, cores, and clumps. Over all scales, the results are inconsistent with complete Jeans fragmentation based on either thermal or thermal plus non-thermal motions. They are more nearly consistent with inefficient thermal Jeans fragmentation, where the thermal Jeans efficiency increases from the largest to the smallest scale.

Stellar feedback created by radiation and winds from massive stars plays a significant role in both physical and chemical evolution of molecular clouds. This energy and momentum leaves an identifiable signature ("bubbles") that affect the dynamics and structure of the cloud. Most bubble searches are performed "by-eye", which are usually time-consuming, subjective and difficult to calibrate. Automatic classifications based on machine learning make it possible to perform systematic, quantifiable and repeatable searches for bubbles. We employ a previously developed machine learning algorithm, Brut, and quantitatively evaluate its performance in identifying bubbles using synthetic dust observations. We adopt magneto-hydrodynamics simulations, which model stellar winds launching within turbulent molecular clouds, as an input to generate synthetic images. We use a publicly available three-dimensional dust continuum Monte-Carlo radiative transfer code, hyperion, to generate synthetic images of bubbles in three Spitzer bands (4.5 um, 8 um and 24 um). We designate half of our synthetic bubbles as a training set, which we use to train Brut along with citizen-science data from the Milky Way Project. We then assess Brut's accuracy using the remaining synthetic observations. We find that after retraining Brut's performance increases significantly, and it is able to identify yellow bubbles, which are likely associated with B-type stars. Brut continues to perform well on previously identified high-score bubbles, and over 10% of the Milky Way Project bubbles are reclassified as high-confidence bubbles, which were previously marginal or ambiguous detections in the Milky Way Project data. We also investigate the size of the training set, dust model, evolution stage and background noise on bubble identification.

Understanding the early stages of star formation is a research field of ongoing development, both theoretically and observationally. In this context, molecular data have been continuously providing observational constraints on the gas dynamics at different excitation conditions and depths in the sources. We have investigated the Barnard 59 core, the only active site of star formation in the Pipe Nebula, to achieve a comprehensive view of the kinematic properties of the source. These information were derived by simultaneously fitting ammonia inversion transition lines (1,1) and (2,2). Our analysis unveils the imprint of protostellar feedback, such as increasing line widths, temperature and turbulent motions in our molecular data. Combined with complementary observations of dust thermal emission, we estimate that the core is gravitationally bound following a virial analysis. If the core is not contracting, another source of internal pressure, most likely the magnetic field, is supporting it against gravitational collapse and limits its star formation efficiency.

We use Green Bank Ammonia Survey observations of NH$_3$ (1,1) and (2,2) emission with 32'' FWHM resolution from a ~ 10 pc$^{2}$ portion of the Cepheus-L1251 molecular cloud to identify hierarchical dense gas structures. Our dendrogram analysis of the NH$_3$ data results in 22 top-level structures, which reside within 13 lower-level, parent structures. The structures are compact (0.01 pc $\lesssim R_{eff} \lesssim$ 0.1 pc) and are spatially correlated with the highest H$_2$ column density portions of the cloud. We also compare the ammonia data to a catalog of dense cores identified by higher-resolution (18.2'' FWHM) Herschel Space Observatory observations of dust continuum emission from Cepheus-L1251. Maps of kinetic gas temperature, velocity dispersion, and NH$_3$ column density, derived from detailed modeling of the NH$_3$ data, are used to investigate the stability and chemistry of the ammonia-identified and Herschel-identified structures. We show that the dust and dense gas in the structures have similar temperatures, with median $T_{dust}$ and $T_K$ measurements of 11.7 $\pm$ 1.1 K and 10.3 $\pm$ 2.0 K, respectively. Based on a virial analysis, we find that the ammonia-identified structures are gravitationally dominated, yet may be in or near a state of virial equilibrium. Meanwhile, the majority of the Herschel-identified dense cores appear to be not bound by their own gravity and instead confined by external pressure. CCS $(2_0-1_0)$ and HC$_5$N $(9-8)$ emission from the region reveal broader line widths and centroid velocity offsets when compared to the NH$_3$ (1,1) emission in some cases, likely due to these carbon-based molecules tracing the turbulent outer layers of the dense cores.

We use CARMA 3mm continuum and molecular lines (NH2D, N2H+, HCO+, HCN and CS) at ~1000 au resolution to characterize the structure and kinematics of the envelope surrounding the deeply embedded first core candidate Per-Bolo 58. The line profile of the observed species shows two distinct peaks separated by 0.4-0.6 km/s, most likely arising from two different optically thin velocity components rather than the product of self-absorption in an optically thick line. The two velocity components, each with a mass of ~0.5-0.6 Msun, overlap spatially at the position of the continuum emission, and produce a general gradient along the outflow direction. We investigate whether these observations are consistent with infall in a turbulent and magnetized envelope. We compare the morphology and spectra of the N2H+(1-0) with synthetic observations of an MHD simulation that considers the collapse of an isolated core that is initially perturbed with a turbulent field. The proposed model matches the data in the production of two velocity components, traced by the isolated hyperfine line of the N2H+(1-0) spectra and shows a general agreement in morphology and velocity field. We also use large maps of the region to compare the kinematics of the core with that of the surrounding large-scale filamentary structure and find that accretion from the large-scale filament could also explain the complex kinematics exhibited by this young dense core.

The star-forming efficiency of dense gas is thought to be set within cores by outflow and radiative feedback. We use magneto-hydrodynamic simulations to investigate the relation between protostellar outflow evolution, turbulence and star formation efficiency. We model the collapse and evolution of isolated dense cores for $\gtrsim$ 0.5 Myr including the effects of turbulence, radiation transfer, and both radiation and outflow feedback from forming protostars. We show that outflows drive and maintain turbulence in the core environment even with strong initial fields. The star-formation efficiency decreases with increasing field strength, and the final efficiencies are $15-40$\%. The Stage 0 lifetime, during which the protostellar mass is less than the dense envelope, increases proportionally with the initial magnetic field strength and ranges from $\sim 0.1-0.4$ Myr. The average accretion rate is well-represented by a tapered turbulent core model, which is a function of the final protostellar mass and is independent of the magnetic field strength. By tagging material launched in the outflow, we demonstrate that the outflow entrains about $3$ times the actual launched gas mass, a ratio that remains roughly constant in time regardless of the initial magnetic field strength. However, turbulent driving increases for stronger fields since momentum is more efficiently imparted to non-outflow material. The protostellar outflow momentum is highest during the first 0.1 Myr and declines thereafter by a factor of $\gtrsim 10$ as the accretion rate diminishes.

We use gas temperature and velocity dispersion data from the Green Bank Ammonia Survey and core masses and sizes from the James Clerk Maxwell Telescope Gould Belt Survey to estimate the virial states of dense cores within the Orion A molecular cloud. Surprisingly, we find that almost none of the dense cores are sufficiently massive to be bound when considering only the balance between self-gravity and the thermal and non-thermal motions present in the dense gas. Including the additional pressure binding imposed by the weight of the ambient molecular cloud material and additional smaller pressure terms, however, suggests that most of the dense cores are pressure confined.

We present new Submillimeter Array (SMA) observations of CO(2-1) outflows toward young, embedded protostars in the Perseus molecular cloud as part of the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. For 57 Perseus protostars, we characterize the orientation of the outflow angles and compare them with the orientation of the local filaments as derived from $Herschel$ observations. We find that the relative angles between outflows and filaments are inconsistent with purely parallel or purely perpendicular distributions. Instead, the observed distribution of outflow-filament angles are more consistent with either randomly aligned angles or a mix of projected parallel and perpendicular angles. A mix of parallel and perpendicular angles requires perpendicular alignment to be more common by a factor of $\sim$3. Our results show that the observed distributions probably hold regardless of the protostar's multiplicity, age, or the host core's opacity. These observations indicate that the angular momentum axis of a protostar may be independent of the large-scale structure. We discuss the significance of independent protostellar rotation axes in the general picture of filament-based star formation.

We present a statistical framework to compare spectral-line data cubes of molecular clouds and use the framework to perform an analysis of various statistical tools developed from methods proposed in the literature. We test whether our methods are sensitive to changes in the underlying physical properties of the clouds or whether their behaviour is governed by random fluctuations. We perform a set of 32 self-gravitating magnetohydrodynamic simulations that test all combinations of five physical parameters -- Mach number, plasma parameter, virial parameter, driving scales, and solenoidal driving fraction -- each of which can be set to a low or high value. We create mock observational data sets of ${\rm ^{13}CO}$(1-0) emission from each simulation. We compare these mock data to a those generated from a set of baseline simulations using pseudo-distance metrics based on 18 different statistical techniques that have previously been used to study molecular clouds. We analyze these results using methods from the statistical field of experimental design and find that several of the statistics can reliably track changes in the underlying physics. Our analysis shows that the interactions between parameters are often among the most significant effects. A small fraction of statistics are also sensitive to changes in magnetic field properties. We use this framework to compare the set of simulations to observations of three nearby star-forming regions: NGC 1333, Oph A, and IC 348. We find that no one simulation agrees significantly better with the observations, although there is evidence that the high Mach number simulations are more consistent with the observations.

We present an overview of the first data release (DR1) and first-look science from the Green Bank Ammonia Survey (GAS). GAS is a Large Program at the Green Bank Telescope to map all Gould Belt star-forming regions with $A_V \gtrsim 7$ mag visible from the northern hemisphere in emission from NH$_3$ and other key molecular tracers. This first release includes the data for four regions in Gould Belt clouds: B18 in Taurus, NGC 1333 in Perseus, L1688 in Ophiuchus, and Orion A North in Orion. We compare the NH$_3$ emission to dust continuum emission from Herschel, and find that the two tracers correspond closely. NH$_3$ is present in over 60\% of lines-of-sight with $A_V \gtrsim 7$ mag in three of the four DR1 regions, in agreement with expectations from previous observations. The sole exception is B18, where NH$_3$ is detected toward ~ 40\% of lines-of-sight with $A_V \gtrsim 7$ mag. Moreover, we find that the NH$_3$ emission is generally extended beyond the typical 0.1 pc length scales of dense cores. We produce maps of the gas kinematics, temperature, and NH$_3$ column densities through forward modeling of the hyperfine structure of the NH$_3$ (1,1) and (2,2) lines. We show that the NH$_3$ velocity dispersion, ${\sigma}_v$, and gas kinetic temperature, $T_K$, vary systematically between the regions included in this release, with an increase in both the mean value and spread of ${\sigma}_v$ and $T_K$ with increasing star formation activity. The data presented in this paper are publicly available.

Context: Understanding how accretion proceeds is a key question of star formation, with important implications for both the physical and chemical evolution of young stellar objects. In particular, very little is known about the accretion variability in the earliest stages of star formation. Aims: To characterise protostellar accretion histories towards individual sources by utilising sublimation and freeze-out chemistry of CO. Methods: A sample of 24 embedded protostars are observed with the Submillimeter Array (SMA) in context of the large program "Mass Assembly of Stellar Systems and their Evolution with the SMA" (MASSES). The size of the C$^{18}$O emitting region, where CO has sublimated into the gas-phase, is measured towards each source and compared to the expected size of the region given the current luminosity. The SMA observations also include 1.3 mm continuum data, which are used to investigate whether a link can be established between accretion bursts and massive circumstellar disks. Results: Depending on the adopted sublimation temperature of the CO ice, between 20% and 50% of the sources in the sample show extended C$^{18}$O emission indicating that the gas was warm enough in the past that CO sublimated and is currently in the process of refreezing; something which we attribute to a recent accretion burst. Given the fraction of sources with extended C$^{18}$O emission, we estimate an average interval between bursts of 20000-50000 yr, which is consistent with previous estimates. No clear link can be established between the presence of circumstellar disks and accretion bursts, however the three closest known binaries in the sample (projected separations <20 AU) all show evidence of a past accretion burst, indicating that close binary interactions may also play a role in inducing accretion variability.

Compact substructure is expected to arise in a starless core as mass becomes concentrated in the central region likely to form a protostar. Additionally, multiple peaks may form if fragmentation occurs. We present ALMA Cycle 2 observations of 60 starless and protostellar cores in the Ophiuchus molecular cloud. We detect eight compact substructures which are >15 arcsec from the nearest Spitzer YSO. Only one of these has strong evidence for being truly starless after considering ancillary data, e.g., from Herschel and X-ray telescopes. An additional extended emission structure has tentative evidence for starlessness. The number of our detections is consistent with estimates from a combination of synthetic observations of numerical simulations and analytical arguments. This result suggests that a similar ALMA study in the Chamaeleon I cloud, which detected no compact substructure in starless cores, may be due to the peculiar evolutionary state of cores in that cloud.

The role of magnetic fields in the early stages of star formation is not well constrained. In order to discriminate between different star formation models, we analyze 3D magnetohydrodynamic simulations of low-mass cores and explore the correlation between magnetic field orientation and outflow orientation over time. We produce synthetic observations of dust polarization at resolutions comparable to millimeter-wave dust polarization maps observed by CARMA and compare these with 2D visualizations of projected magnetic field and column density. Cumulative distribution functions of the projected angle between the magnetic field and outflow show different degrees of alignment in simulations with differing mass-to-flux ratios. The distribution function for the less magnetized core agrees with observations finding random alignment between outflow and field orientations, while the more magnetized core exhibits stronger alignment. We find that fractional polarization increases when the system is viewed such that the magnetic field is close to the plane of the sky, and the values of fractional polarization are consistent with observational measurements. The simulation outflow, which reflects the underlying angular momentum of the accreted gas, changes direction significantly over the first $\sim0.1$ Myr of evolution. This movement could lead to the observed random alignment between outflows and the magnetic fields in protostellar cores.

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 molecular clouds, including stellar winds, with a range of stellar mass-loss rates and magnetic field strengths. We generate synthetic $^{12}$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, including principal component analysis, the spectral correlation function and the velocity coordinate spectrum, are sensitive to changes in stellar mass-loss rates and/or time evolution. A few statistics, including the Cramer statistic and velocity coordinate spectrum, are sensitive to the magnetic field strength. These findings demonstrate that stellar feedback influences molecular cloud turbulence and can be identified and quantified observationally using such statistics.

We present ALMA Cycle 1 observations of the HH46/47 molecular outflow using combined 12m array and ACA observations. The improved angular resolution and sensitivity of our multi-line maps reveal structures that help us study the entrainment process in much more detail and allow us to obtain more precise estimates of outflow properties than previous observations. We use 13CO(1-0) and C18O(1-0) emission to correct for the 12CO(1-0) optical depth to accurately estimate the outflow mass, momentum and kinetic energy. This correction increases the estimates of the mass, momentum and kinetic energy by factors of about 9, 5 and 2, respectively, with respect to estimates assuming optically thin emission. The new 13CO and C18O data also allow us to trace denser and slower outflow material than that traced by the 12CO maps, and they reveal an outflow cavity wall at very low velocities (as low as 0.2km/s with respect to the cores central velocity). Adding with the slower material traced only by 13CO and C18O, there is another factor of 3 increase in the mass estimate and 50% increase in the momentum estimate. The estimated outflow properties indicate that the outflow is capable of dispersing the parent core within the typical lifetime of the embedded phase of a low-mass protostar, and that it is responsible for a core-to-star efficiency of 1/4 to 1/3. We find that the outflow cavity wall is composed of multiple shells associated with a series of jet bow-shock events. Within about 3000AU of the protostar the 13CO and C18O emission trace a circumstellar envelope with both rotation and infall motions, which we compare with a simple analytic model. The CS(2-1) emission reveals tentative evidence of a slowly-moving rotating outflow, which we suggest is entrained not only poloidally but also toroidally by a disk wind that is launched from relatively large radii from the source.

The protostellar outflows of wide-separation forming binaries frequently appear misaligned. We use magneto-hydrodynamic simulations to investigate the alignment of protostellar spin and molecular outflows for forming binary pairs. We show that the protostellar pairs, which form from turbulent fragmentation within a single parent core, have randomly oriented angular momenta. Although the pairs migrate to closer separations, their spins remain partially misaligned. We produce $^{12}$CO(2-1) synthetic observations of the simulations and characterize the outflow orientation in the emission maps. The CO-identified outflows exhibit a similar random distribution and are also statistically consistent with the observed distribution of molecular outflows. We conclude that observed misalignment provides a clear signature of binary formation via turbulent fragmentation. The persistence of misaligned outflows and stellar spins following dynamical evolution may provide a signature of binary origins for more evolved multiple star systems.

We present an Atacama Large Millimeter/submillimeter Array (ALMA) 106 GHz (Band 3) continuum survey of the complete population of dense cores in the Chamaeleon I molecular cloud. We detect a total of 24 continuum sources in 19 different target fields. All previously known Class 0 and Class I protostars in Chamaeleon I are detected, whereas all of the 56 starless cores in our sample are undetected. We show that the Spitzer+Herschel census of protostars in Chamaeleon I is complete, with the rate at which protostellar cores have been misclassified as starless cores calculated as <1/56, or < 2%. We use synthetic observations to show that starless cores collapsing following the turbulent fragmentation scenario are detectable by our ALMA observations when their central densities exceed ~10^8 cm^-3, with the exact density dependent on the viewing geometry. Bonnor-Ebert spheres, on the other hand, remain undetected to central densities at least as high as 10^10 cm^-3. Our starless core non-detections are used to infer that either the star formation rate is declining in Chamaeleon I and most of the starless cores are not collapsing, matching the findings of previous studies, or that the evolution of starless cores are more accurately described by models that develop less substructure than predicted by the turbulent fragmentation scenario, such as Bonnor-Ebert spheres. We outline future work necessary to distinguish between these two possibilities.

We investigate the alignment between outflow axes in nine of the youngest binary/multiple systems in the Perseus Molecular Cloud. These systems have typical member spacing larger than 1000 AU. For outflow identification, we use 12CO(2-1) and 12CO(3-2) data from a large survey with the Sub- millimeter Array: Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES). The distribution of outflow orientations in the binary pairs is consistent with random or preferentially anti-aligned distributions, demonstrating that these outflows are misaligned. This result suggests that these systems are possibly formed in environments where the distribution of angular momentum is complex and disordered, and these systems do not come from the same co-rotating structures or from an initial cloud with aligned vectors of angular momentum.

We present continuum and molecular line observations at 230 GHz and 345 GHz from the Sub-millimeter Array (SMA) toward three protostars in the Perseus L1448N region. The data are from the large project "Mass Assembly of Stellar Systems and their Evolution with the SMA" (MASSES). Three dust continuum sources, Source B, Source NW, and Source A, are detected at both frequencies. These sources have corresponding emission peaks in C18O (J=2-1), 13CO (J=2-1), and HCO+ (J=4-3), and have offsets with N2D+ (J=3-2) peaks. High angular resolution data from a complimentary continuum survey with the Karl G. Jansky Very Large Array show that Source B is associated with three 8 mm continuum objects, Source NW with two, and Source A remains single. These results suggest that multiplicity in L1448N exists at different spatial scales from a few thousand AU to < 100 AU. Velocity gradients in each source obtained from two-dimensional fits to the SMA C18O emission are found to be perpendicular to within 20 degrees of the outflow directions as revealed by 12CO (J=2-1). We have observed that Sources B and NW with multiplicity have higher densities than Source A without multiplicity. This suggests that thermal Jeans fragmentation can be relevant in the fragmentation process. However, we have not observed a difference in the ratio between rotational and gravitational energy between sources with and without multiplicity. We also have not observed a trend between non-thermal velocity dispersions and the level of fragmentation. Our study has provided the first direct and comprehensive comparison between multiplicity and core properties in low-mass protostars, although based on small number statistics.

The stellar initial mass function (IMF) is a fundamental astrophysical quantity that impacts a wide range of astrophysical problems from heavy element distribution to galactic evolution to planetary system formation. However, the origin and universality of the IMF are hotly debated both observationally and theoretically. I review recent observations of the IMF across a variety of environments. These suggest the IMF is surprisingly invariant between star-forming regions, star clusters, and spiral galaxies but that it may also vary under extreme conditions, including within the Galactic center and early type galaxies.

Observations of nearby molecular clouds detect "shells", which are likely caused by winds from young main sequence stars. However, the progenitors of these observed features are not well characterized and the mass-loss rates inferred from the gas kinematics are several orders of magnitude greater than those predicted by atomic line-driven stellar wind models. We use magnetohydrodynamic simulations to model winds launching within turbulent molecular clouds and explore the impact of wind properties on cloud morphology and turbulence. We find that winds do not produce clear features in turbulent statistics such as the Fourier spectra of density and momentum but do impact the Fourier velocity spectrum. The density and velocity distribution functions, especially as probed by CO spectral lines, strongly indicate the presence and influence of winds. We show that stellar mass-loss rates for individual stars must be $\dot m_w \gtrsim 10^{-7}$ Msun yr$^{-1}$, similar to those estimated from observations, to reproduce shell properties. Consequently, we conclude that B and A-type main sequence stars have mass-loss rates several orders of magnitude larger that those predicted by models or that young stars are more variable than expected due to magnetic activity or accretion.

The core velocity dispersion (CVD) is a potentially useful tool for studying the turbulent velocity field of molecular clouds. CVD is based on centroid velocities of dense gas clumps, thus is less prone to density fluctuation and reflects more directly the cloud velocity field. Prior work demonstrated that the Taurus molecular cloud CVD resembles the well-known Larson's linewidth-size relation of molecular clouds. In this work, we studied the dependence of the CVD on the line-of-sight thickness of molecular clouds, a quantity which cannot be measured by direct means. We produced a simple statistical model of cores within clouds and analyzed the CVD of a variety of hydrodynamical simulations. We show that the relation between the CVD and the 2D projected separation of cores ($L_{2D}$) is sensitive to the cloud thickness. When the cloud is thin, the index of CVD-$L_{2D}$ relation ($\gamma$ in the relation CVD$\sim L_{2D}^{\gamma}$) reflects the underlying energy spectrum ($E(k)\sim k^{-\beta}$) in that $\gamma\sim(\beta-1)/2$. The CVD-$L_{2D}$ relation becomes flatter ($\gamma\to 0$) for thicker clouds. We used this result to constrain the thicknesses of Taurus, Perseus, and Ophiuchus. We conclude that Taurus has a ratio of cloud depth to cloud length smaller than about 1/10-1/8, i.e. it is a sheet. A simple geometric model fit to the linewidth-size relation indicates that the Taurus cloud has a $\sim 0.7$ pc line-of-sight dimension. In contrast, Perseus and Ophiuchus are thicker and have ratios of cloud depth to cloud length larger than about 1/10-1/8.

Protostellar outflows provide a means to probe the accretion process of forming stars and their ability to inject energy into their surroundings. However, conclusions based on outflow observations depend upon the degree of accuracy with which their properties can be estimated. We examine the quality of Atacama Large Millimeter/submillimeter Array (ALMA) observations of protostellar outflows by producing synthetic $^{12}$CO(1-0) and $^{13}$CO(1-0) observations of numerical simulations. We use various ALMA configurations, observational parameters, and outflow inclinations to assess how accurately different assumptions and setups can recover underlying properties. We find that more compact arrays and longer observing times can improve the mass and momentum recovery by a factor of two. During the first $\sim$0.3 Myr of evolution, $^{12}$CO(1-0) is optically thick, even for velocities $|v|\ge 1$ km s$^{-1}$, and outflow mass is severely underestimated without an optical depth correction. Likewise, $^{13}$CO(1-0) is optically thick during the first $\simeq 0.1$ Myr. However, underestimation due to shorter observing time, missing flux, and optical depth are partially offset by the assumption of local thermodynamic equilibrium and higher excitation temperatures. Overall, we expect that full ALMA $^{13}$CO(1-0) observations of protostellar sources within 500 pc with observing times $\gtrsim 1$ hrs and assumed excitation temperatures of $T<20$K will reliably measure mass and line-of-sight momentum to within 20\%.

We investigate the spectral correlations between different species used to observe molecular clouds. We use hydrodynamic simulations and a full chemical network to study the abundances of over 150 species in typical Milky Way molecular clouds. We perform synthetic observations in order to produce emission maps of a subset of these tracers. We study the effects of different lines of sight and spatial resolution on the emission distribution and perform a robust quantitative comparison of the species to each other. We use the Spectral Correlation Function (SCF), which quantifies the root mean squared difference between spectra separated by some length scale, to characterize the structure of the simulated cloud in position-position-velocity (PPV) space. We predict the observed SCF for a broad range of observational tracers, and thus, identify homologous species. In particular, we show that the pairs C and CO, C$^{+}$ and CN, NH$_3$ and H$_2$CS have very similar SCFs. We measure the SCF slope variation as a function of beam size for all species and demonstrate that the beam size has a distinct effect on different species emission. However, for beams of up to 10", placing the cloud at 1 kpc, the change is not large enough to move the SCF slopes into different regions of parameter space. The results from this study provide observational guidance for choosing the best tracer to probe various cloud length scales.

We have used deep near-infrared observations with adaptive optics to discover a distributed population of low-mass protostars within the filamentary Infrared Dark Cloud G34.43+00.24. We use maps of dust emission at multiple wavelengths to determine the column density structure of the cloud. In combination with an empirically-verified model of the magnitude distribution of background stars, this column density map allows us to reliably determine overdensities of red sources that are due to embedded protostars in the cloud. We also identify protostars through their extended emission in K-band which comes from excited H2 in protostellar outflows or reflection nebulosity. We find a population of distributed low-mass protostars, suggesting that low-mass protostars may form earlier than, or contemporaneously with, high-mass protostars in such a filament. The low-mass protostellar population may also produce the narrow linewidth SiO emission observed in some clouds without high-mass protostars. Finally, we use a molecular line map of the cloud to determine the virial parameter per unit length along the filament and find that the highest mass protostars form in the most bound portion of the filament, as suggested by theoretical models.

We review current theories for the origin of the Stellar Initial Mass Function (IMF) with particular focus on the extent to which the IMF can be considered universal across various environments. To place the issue in an observational context, we summarize the techniques used to determine the IMF for different stellar populations, the uncertainties affecting the results, and the evidence for systematic departures from universality under extreme circumstances. We next consider theories for the formation of prestellar cores by turbulent fragmentation and the possible impact of various thermal, hydrodynamic and magneto-hydrodynamic instabilities. We address the conversion of prestellar cores into stars and evaluate the roles played by different processes: competitive accretion, dynamical fragmentation, ejection and starvation, filament fragmentation and filamentary accretion flows, disk formation and fragmentation, critical scales imposed by thermodynamics, and magnetic braking. We present explanations for the characteristic shapes of the Present-Day Prestellar Core Mass Function and the IMF and consider what significance can be attached to their apparent similarity. Substantial computational advances have occurred in recent years, and we review the numerical simulations that have been performed to predict the IMF directly and discuss the influence of dynamics, time-dependent phenomena, and initial conditions.

We propose a new approach to comparing simulated observations that enables us to determine the significance of the underlying physical effects. We utilize the methodology of experimental design, a subfield of statistical analysis, to establish a framework for comparing simulated position-position-velocity data cubes to each other. We propose three similarity metrics based on methods described in the literature: principal component analysis, the spectral correlation function, and the Cramer multi-variate two sample similarity statistic. Using these metrics, we intercompare a suite of mock observational data of molecular clouds generated from magnetohydrodynamic simulations with varying physical conditions. Using this framework, we show that all three metrics are sensitive to changing Mach number and temperature in the simulation sets, but cannot detect changes in magnetic field strength and initial velocity spectrum. We highlight the shortcomings of one-factor-at-a-time designs commonly used in astrophysics and propose fractional factorial designs as a means to rigorously examine the effects of changing physical properties while minimizing the investment of computational resources.

We explore the utility of CI as an alternative high-fidelity gas mass tracer for Galactic molecular clouds. We evaluate the X$_{\rm CI}$-factor for the 609 $\mu$m carbon line, the analog of the CO X-factor, which is the ratio of the H$_2$ column density to the integrated $^{12}$CO(1-0) line intensity. We use 3D-PDR to post-process hydrodynamic simulations of turbulent, star-forming clouds. We compare the emission of CI and CO for model clouds irradiated by 1 and 10 times the average background and demonstrate that CI is a comparable or superior tracer of the molecular gas distribution for column densities up to $6 \times 10^{23}$ cm$^{-2}$. Our results hold for both reduced and full chemical networks. For our fiducial Galactic cloud we derive an average $X_{\rm CO}$ of $3.0\times 10^{20}$ cm$^{-2}$K$^{-1}$km$^{-1}$s and $X_{\rm CI}$ of $1.1\times 10^{21}$ cm$^{-2}$K$^{-1}$km$^{-1}$s.

Recent observations show that small, young, stellar groupings of ~10 to 40 members tend of have a centrally-located most massive member, reminiscent of mass segregation seen in large clustered systems. Here, we analyze hydrodynamic simulations which form small clusters and analyze their properties in a manner identical to the observations. We find that the simulated clusters possess similar properties to the observed clusters, including a tendency to exhibit mass segregation. In the simulations, the central location of the most massive member is not due to dynamical evolution, since there is little interaction between the cluster members. Instead, the most massive cluster member appears to form at the center. We also find that the more massive stars in the cluster form at slightly earlier times.

We investigate protostellar outflow evolution, gas entrainment, and star formation efficiency using radiation-hydrodynamic simulations of isolated, turbulent low-mass cores. We adopt an X-wind launching model, in which the outflow rate is coupled to the instantaneous protostellar accretion rate and evolution. We vary the outflow collimation angle from $\theta$=0.01-0.1 and find that even well collimated outflows effectively sweep up and entrain significant core mass. The Stage 0 lifetime ranges from 0.14-0.19 Myr, which is similar to the observed Class 0 lifetime. The star formation efficiency of the cores spans 0.41-0.51. In all cases, the outflows drive strong turbulence in the surrounding material. Although the initial core turbulence is purely solenoidal by construction, the simulations converge to approximate equipartition between solenoidal and compressive motions due to a combination of outflow driving and collapse. When compared to a simulation of a cluster of protostars, which is not gravitationally centrally condensed, we find that the outflows drive motions that are mainly solenoidal. The final turbulent velocity dispersion is about twice the initial value of the cores, indicating that an individual outflow is easily able to replenish turbulent motions on sub-parsec scales. We post-process the simulations to produce synthetic molecular line emission maps of $^{12}$CO, $^{13}$CO, and C$^{18}$O and evaluate how well these tracers reproduce the underlying mass and velocity structure.

Interpreting the nature of starless cores has been a prominent goal in star formation for many years. In order to characterise the evolutionary stages of these objects, we perform synthetic observations of a numerical simulation of a turbulent molecular cloud. We find that nearly all cores that we detect are associated with filaments and eventually form protostars. We conclude that observed starless cores which appear Jeans unstable are only marginally larger than their respective Jeans masses (within a factor of 3). We note single dish observations such as those performed with the JCMT appear to miss significant core structure on small scales due to beam averaging. Finally, we predict that interferometric observations with ALMA Cycle 1 will resolve the important small scale structure, which has so far been missed by mm-wavelength observations.

Stars form from the gravitational collapse of dense molecular cloud cores. In the protostellar phase, mass accretes from the core onto a protostar, likely through an accretion disk, and it is during this phase that the initial masses of stars and the initial conditions for planet formation are set. Over the past decade, new observational capabilities provided by the Spitzer Space Telescope and Herschel Space Observatory have enabled wide-field surveys of entire star-forming clouds with unprecedented sensitivity, resolution, and infrared wavelength coverage. We review resulting advances in the field, focusing both on the observations themselves and the constraints they place on theoretical models of star formation and protostellar evolution. We also emphasize open questions and outline new directions needed to further advance the field.

The physical properties of molecular clouds are often measured using spectral-line observations, which provide the only probes of the clouds' velocity structure. It is hard, though, to assess whether and to what extent intensity features in position-position-velocity (PPV) space correspond to "real" density structures in position-position-position (PPP) space. In this paper, we create synthetic molecular cloud spectral-line maps of simulated molecular clouds, and present a new technique for measuring the reality of individual PPV structures. Our procedure projects density structures identified in PPP space into corresponding intensity structures in PPV space and then measures the geometric overlap of the projected structures with structures identified from the synthetic observation. The fractional overlap between a PPP and PPV structure quantifies how well the synthetic observation recovers information about the 3D structure. Applying this machinery to a set of synthetic observations of CO isotopes, we measure how well spectral-line measurements recover mass, size, velocity dispersion, and virial parameter for a simulated star-forming region. By disabling various steps of our analysis, we investigate how much opacity, chemistry, and gravity affect measurements of physical properties extracted from PPV cubes. For the simulations used here, our results suggest that superposition induces a ~40% uncertainty in masses, sizes, and velocity dispersions derived from 13CO. The virial parameter is most affected by superposition, such that estimates of the virial parameter derived from PPV and PPP information typically disagree by a factor of ~2. This uncertainty makes it particularly difficult to judge whether gravitational or kinetic energy dominate a given region, since the majority of virial parameter measurements fall within a factor of 2 of the equipartition level alpha ~ 2.

We present APEX observations of C17O(2-1), N2H+(3-2), and N2D+(3-2) towards the subfragments inside the prestellar core SMM 6 in Orion B9. We combined these spectral line data with our previous SABOCA 350-{\mu}m dust continuum map of the source. The subfragments are characterised by subsonic internal non-thermal motions ({\sigma}NT~0.5cs), and most of them appear to be gravitationally bound. The dispersion of the N2H+ velocity centroids among the condensations is very low (0.02 km/s). The CO depletion factors we derive, fD=0.8+/-0.4 - 3.6+/-1.5, do not suggest any significant CO freeze-out but this may be due to the canonical CO abundance we adopt. The fractional abundances of N2H+ and N2D+ with respect to H2 are found to be ~0.9-2.3x10^-9 and ~4.9-9.9x10^-10, respectively. The deuterium fractionation of N2H+, or the N2D+/N2H+ column density ratio, lies in the range 0.30+/-0.07 - 0.43+/-0.09. The detected substructure inside SMM 6 is likely the result of cylindrical Jeans-type gravitational fragmentation. We estimate the timescale for this fragmentation to be ~1.8x10^5 yr. The condensations are unlikely to be able to interact with one another and coalesce before local gravitational collapse ensues. Moreover, significant mass growth of the condensations via competitive-like accretion from the parent core seems unfeasible. The high level of molecular deuteration in the condensations suggests that gas-phase CO should be strongly depleted. It also points towards an advanced stage of chemical evolution. The subfragments of SMM 6 might therefore be near the onset of gravitational collapse, but whether they can form protostellar or substellar objects (brown dwarfs) depends on the local star formation efficiency and remains to be clarified.

We use 3D-PDR, a three-dimensional astrochemistry code for modeling photodissociation regions (PDRs), to post-process hydrodynamic simulations of turbulent, star-forming clouds. We focus on the transition from atomic to molecular gas, with specific attention to the formation and distribution of H, C+, C, H2 and CO. First, we demonstrate that the details of the cloud chemistry and our conclusions are insensitive to the simulation spatial resolution, to the resolution at the cloud edge, and to the ray angular resolution. We then investigate the effect of geometry and simulation parameters on chemical abundances and find weak dependence on cloud morphology as dictated by gravity and turbulent Mach number. For a uniform external radiation field, we find similar distributions to those derived using a one-dimensional PDR code. However, we demonstrate that a three-dimensional treatment is necessary for a spatially varying external field, and we caution against using one-dimensional treatments for non-symmetric problems. We compare our results with the work of Glover et al. (2010), who self-consistently followed the time evolution of molecule formation in hydrodynamic simulations using a reduced chemical network. In general, we find good agreement with this in situ approach for C and CO abundances. However, the temperature and H2 abundances are discrepant in the boundary regions (Av < 5), which is due to the different number of rays used by the two approaches.

We present a 350 micron APEX/SABOCA map of IRAS 05399-0121/SMM 1, which is a dense double-core system in Orion B9. We combined these data with our previous LABOCA 870-micron data. The spatial resolution of the new SABOCA image, ~3400 AU, is about 2.6 times better than provided by LABOCA. We also make use of Spitzer infrared observations to characterise the star-formation activity in the source. The source is filamentary and remains a double-core system on the 3400 AU scale probed here, where the projected separation between IRAS 05399 and SMM 1 is 0.14 pc. The broadband spectral energy distribution of IRAS 05399 suggests that it is near the Stage 0/I borderline. A visual inspection of the Spitzer/IRAC images provides hints of a quadrupolar-like jet morphology around IRAS 05399, supporting the possibility that it is a binary system. The temperature map reveals warm spots towards IRAS 05399 and the southeastern tip of the source. These features are likely to be imprints of protostellar or shock heating, while external heating could be provided by the nearby high-mass star-forming region NGC 2024. A simple analysis suggests that the density profile at the position of SMM 1 has the form ~r^-(2.3_{-0.9}^{+2.2}). The source splitting into two subcores along the long axis can be explained by cylindrical Jeans-type fragmentation but the steepness of the density profile is shallower than what is expected for an isothermal cylinder. The difference between the evolutionary stages of IRAS 05399 (protostellar) and SMM 1 (starless) suggests that the former has experienced a phase of rapid mass accretion, supported by the very long outflow it drives. The protostellar jet from IRAS 05399 might have influenced the nearby core SMM 1.

The properties of unresolved protostars and their local environment (e.g., disk, envelope and outflow characteristics) are frequently inferred from spectral energy distributions (SEDs) through comparison with idealized model SEDs. However, if it is not possible to image a source and its environment directly, it is difficult to constrain and evaluate the accuracy of these derived properties. In this proceeding, I present a brief overview of the reliability of SED modeling by analyzing dust continuum synthetic observations of realistic simulations.

The properties of unresolved protostars and their local environment are frequently inferred from spectral energy distributions (SEDs) using radiative transfer modeling. We perform synthetic observations of realistic star formation simulations to evaluate the accuracy of properties inferred from fitting model SEDs to observations. We use ORION, an adaptive mesh refinement (AMR) three-dimensional gravito-radiation-hydrodynamics code, to simulate low-mass star formation in a turbulent molecular cloud including the effects of protostellar outflows. To obtain the dust temperature distribution and SEDs of the forming protostars, we post-process the simulations using HYPERION, a state-of-the-art Monte-Carlo radiative transfer code. We find that the ORION and HYPERION dust temperatures typically agree within a factor of two. We compare synthetic SEDs of embedded protostars for a range of evolutionary times, simulation resolutions, aperture sizes, and viewing angles. We demonstrate that complex, asymmetric gas morphology leads to a variety of classifications for individual objects as a function of viewing angle. We derive best-fit source parameters for each SED through comparison with a pre-computed grid of radiative transfer models. While the SED models correctly identify the evolutionary stage of the synthetic sources as embedded protostars, we show that the disk and stellar parameters can be very discrepant from the simulated values. Parameters such as the stellar accretion rate, stellar mass, and disk mass show better agreement, but can still deviate significantly, and the agreement may in some cases be artificially good due to the limited range of parameters in the set of model SEDs. Lack of correlation between the model and simulation properties in many individual instances cautions against over-interpreting properties inferred from SEDs for unresolved protostellar sources. (Abridged)

Determining the initial stellar multiplicity is a challenging problem since protostars are faint and deeply embedded at early times; once formed, multiple protostellar systems may significantly dynamically evolve before they are optically revealed. Interferometers such as CARMA and ALMA make it possible to probe the scales at which turbulent fragmentation occurs in dust continuum emission, potentially constraining early stellar multiplicity. In this Letter we present synthetic observations of starless and protostellar cores undergoing fragmentation on scales of a few thousand AU to produce wide binary systems. We show that interferometric observations of starless cores by CARMA should be predominantly featureless at early stages, although wide protostellar companions should be apparent. The enhanced capabilities of ALMA improve the detection of core morphology so that it may be possible to detect substructure at earlier times. In either case, spatial filtering from interferometry reduces the observed core substructure and often eradicates traces of existing filamentary morphology on scales down to 0.025 pc. However, some missing structure may be recaptured by combining data from the ALMA full science and Atacama compact arrays.

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.

We revisit the problem of low-mass pre-main-sequence (PMS) stellar evolution and its observational consequences for where stars fall on the Hertzsprung-Russell diagram (HRD). In contrast to most previous work, our models follow stars as they grow from small masses via accretion, and we perform a systematic study of how the stars' HRD evolution is influenced by their initial radius, by the radiative properties of the accretion flow, and by the accretion history, using both simple idealized accretion histories and histories taken from numerical simulations of star cluster formation. We compare our numerical results to both non-accreting isochrones and to the positions of observed stars in the HRD, with a goal of determining whether both the absolute ages and the age dispersions inferred from non-accreting isochrones are reliable. We show that non-accreting isochrones can sometimes overestimate stellar ages for more massive stars (those with effective temperatures above \sim 3500 K), thereby explaining why non-accreting isochrones often suggest a systematic age difference between more and less massive stars in the same cluster. However, we also find the only way to produce a similar overestimate for the ages of cooler stars is by selecting parameters that are strongly inconsistent with both theoretical models and observational constraints. We therefore conclude that inferred ages and age spreads for cool stars are reliable, at least to the extent that the observed bolometric luminosities and temperatures are accurate. Finally, we note that the time-dependence of the mass accretion rate has remarkably little effect on low-mass stars' evolution on the HRD, and that such time-dependence may be neglected for all stars except those with effective temperatures above \sim 4000 K.

The protostellar luminosity function (PLF) is the present-day luminosity function of the protostars in a region of star formation. It is determined using the protostellar mass function (PMF) in combination with a stellar evolutionary model that provides the luminosity as a function of instantaneous and final stellar mass. As in McKee & Offner (2010), we consider three main accretion models: the Isothermal Sphere model, the Turbulent Core model, and an approximation of the Competitive Accretion model. We also consider the effect of an accretion rate that tapers off linearly in time and an accelerating star formation rate. For each model, we characterize the luminosity distribution using the mean, median, maximum, ratio of the median to the mean, standard deviation of the logarithm of the luminosity, and the fraction of very low luminosity objects. We compare the models with bolometric luminosities observed in local star forming regions and find that models with an approximately constant accretion time, such as the Turbulent Core and Competitive Accretion models, appear to agree better with observation than those with a constant accretion rate, such as the Isothermal Sphere model. We show that observations of the mean protostellar luminosity in these nearby regions of low-mass star formation suggest a mean star formation time of 0.3$\pm$0.1 Myr. Such a timescale, together with some accretion that occurs non-radiatively and some that occurs in high-accretion, episodic bursts, resolves the classical "luminosity problem" in low-mass star formation, in which observed protostellar luminosities are significantly less than predicted. An accelerating star formation rate is one possible way of reconciling the observed star formation time and mean luminosity.

Low-mass protostars are less luminous than expected. This luminosity problem is important because the observations appear to be inconsistent with some of the basic premises of star formation theory. Two possible solutions are that stars form slowly, which is supported by recent data, and/or that protostellar accretion is episodic; current data suggest that the latter accounts for less than half the missing luminosity. The solution to the luminosity problem bears directly on the fundamental problem of the time required to form a low-mass star. The protostellar mass and luminosity functions provide powerful tools both for addressing the luminosity problem and for testing theories of star formation. Results are presented for the collapse of singular isothermal spheres, for the collapse of turbulent cores, and for competitive accretion.

Forming stars emit a significant amount of radiation into their natal environment. While the importance of radiation feedback from high-mass stars is widely accepted, radiation has generally been ignored in simulations of low-mass star formation. I use ORION, an adaptive mesh refinement (AMR) three-dimensional gravito-radiation-hydrodynamics code, to model low-mass star formation in a turbulent molecular cloud. I demonstrate that including radiation feedback has a profound effect on fragmentation and protostellar multiplicity. Although heating is mainly confined within the core envelope, it is sufficient to suppress disk fragmentation that would otherwise result in low-mass companions or brown dwarfs. As a consequence, turbulent fragmentation, not disk fragmentation, is likely the origin of low-mass binaries.

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 protostellar mass function (PMF) is the Present-Day Mass Function of the protostars in a region of star formation. It is determined by the initial mass function weighted by the accretion time. The PMF thus depends on the accretion history of protostars and in principle provides a powerful tool for observationally distinguishing different protostellar accretion models. We consider three basic models here: the Isothermal Sphere model (Shu 1977), the Turbulent Core model (McKee & Tan 2003), and an approximate representation of the Competitive Accretion model (Bonnell et al. 1997, 2001a). We also consider modified versions of these accretion models, in which the accretion rate tapers off linearly in time. Finally, we allow for an overall acceleration in the rate of star formation. At present, it is not possible to directly determine the PMF since protostellar masses are not currently measurable. We carry out an approximate comparison of predicted PMFs with observation by using the theory to infer the conditions in the ambient medium in several star-forming regions. Tapered and accelerating models generally agree better with observed star-formation times than models without tapering or acceleration, but uncertainties in the accretion models and in the observations do not allow one to rule out any of the proposed models at present. The PMF is essential for the calculation of the Protostellar Luminosity Function, however, and this enables stronger conclusions to be drawn (Offner & McKee 2010).

The kinematics of newly-formed star clusters are interesting both as a probe of the state of the gas clouds from which the stars form, and because they influence planet formation, stellar mass segregation, cluster disruption, and other processes controlled in part by dynamical interactions in young clusters. However, to date there have been no attempts to use simulations of star cluster formation to investigate how the kinematics of young stars change in response to variations in the properties of their parent molecular clouds. In this letter we report the results of turbulent self-gravitating simulations of cluster formation in which we consider both clouds in virial balance and those undergoing global collapse. We find that stars in these simulations generally have velocity dispersions smaller than that of the gas by a factor of ~ 5, independent of the dynamical state of the parent cloud, so that subvirial stellar velocity dispersions arise naturally even in virialized molecular clouds. The simulated clusters also show large-scale stellar velocity gradients of ~0.2-2 km s$^{-1}$ pc$^{-1}$ and strong correlations between the centroid velocities of stars and gas, both of which are observed in young clusters. We conclude that star clusters should display subvirial velocity dispersions, large-scale velocity gradients, and strong gas-star velocity correlations regardless of whether their parent clouds are in virial balance, and, conversely, that observations of these features cannot be used to infer the dynamical state of the parent gas clouds.

Forming stars emit a substantial amount of radiation into their natal environment. We use ORION, an adaptive mesh refinement (AMR) three-dimensional gravito-radiation-hydrodynamics code, to simulate low-mass star formation in a turbulent molecular cloud. We compare the distribution of stellar masses, accretion rates, and temperatures in the cases with and without radiative transfer, and we demonstrate that radiative feedback has a profound effect on accretion, multiplicity, and mass by reducing the number of stars formed and the total rate at which gas turns into stars. We also show, that once star formation reaches a steady state, protostellar radiation is by far the dominant source of energy in the simulation, exceeding viscous dissipation and compressional heating by at least an order of magnitude. Calculations that omit radiative feedback from protstars significantly underestimate the gas temperature and the strength of this effect. Although heating from protostars is mainly confined to the protostellar cores, we find that it is sufficient to suppress disk fragmentation that would otherwise result in very low-mass companions or brown dwarfs. We demonstrate that the mean protostellar accretion rate increases with the final stellar mass so that the star formation time is only a weak function of mass.

While both society and astronomy have evolved greatly over the past fifty years, the academic institutions and incentives that shape our field have remained largely stagnant. As a result, the astronomical community is faced with several major challenges, including: (1) the training that we provide does not align with the skills that future astronomers will need, (2) the postdoctoral phase is becoming increasingly demanding and demoralizing, and (3) our jobs are increasingly unfriendly to families with children. Solving these problems will require conscious engineering of our profession. Fortunately, this Decadal Review offers the opportunity to revise outmoded practices to be more effective and equitable. The highest priority of the Subcommittee on the State of the Profession should be to recommend specific, funded activities that will ensure the field meets the challenges we describe.

Massive stars produce so much light that the radiation pressure they exert on the gas and dust around them is stronger than their gravitational attraction, a condition that has long been expected to prevent them from growing by accretion. We present three-dimensional radiation-hydrodynamic simulations of the collapse of a massive prestellar core and find that radiation pressure does not halt accretion. Instead, gravitational and Rayleigh-Taylor instabilities channel gas onto the star system through non-axisymmetric disks and filaments that self-shield against radiation, while allowing radiation to escape through optically-thin bubbles. Gravitational instabilities cause the disk to fragment and form a massive companion to the primary star. Radiation pressure does not limit stellar masses, but the instabilities that allow accretion to continue lead to small multiple systems.

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.

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$(2\to 1)$, NH$_3(1,1)$, and N$_2$H$^+(1\to 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.

Molecular clouds are observed to be turbulent, but the origin of this turbulence is not well understood. As a result, there are two different approaches to simulating molecular clouds, one in which the turbulence is allowed to decay after it is initialized, and one in which it is driven. We use the adaptive mesh refinement (AMR) code, Orion, to perform high-resolution simulations of molecular cloud cores and protostars in environments with both driven and decaying turbulence. We include self-gravity, use a barotropic equation of state, and represent regions exceeding the maximum grid resolution with sink particles. We analyze the properties of bound cores such as size, shape, linewidth, and rotational energy, and we find reasonable agreement with observation. At high resolution, the different rates of core accretion in the two cases have a significant effect on protostellar system development. Clumps forming in a decaying turbulence environment produce high-multiplicity protostellar systems with Toomre-Q unstable disks that exhibit characteristics of the competitive accretion model for star formation. In contrast, cores forming in the context of continuously driven turbulence and virial equilibrium form smaller protostellar systems with fewer low-mass members. Our simulations of driven and decaying turbulence show some statistically significant differences, particularly in the production of brown dwarfs and core rotation, but the uncertainties are large enough that we are not able to conclude whether observations favor one or the other.

We present a scheme to solve the nonlinear multigroup radiation diffusion (MGD) equations. The method is incorporated into a massively parallel, multidimensional, Eulerian radiation-hydrodynamic code with adaptive mesh refinement (AMR). The patch-based AMR algorithm refines in both space and time creating a hierarchy of levels, coarsest to finest. The physics modules are time-advanced using operator splitting. On each level, separate level-solve packages advance the modules. Our multigroup level-solve adapts an implicit procedure which leads to a two-step iterative scheme that alternates between elliptic solves for each group with intra-cell group coupling. For robustness, we introduce pseudo transient continuation (PTC). We analyze the magnitude of the PTC parameter to ensure positivity of the resulting linear system, diagonal dominance and convergence of the two-step scheme. For AMR, a level defines a subdomain for refinement. For diffusive processes such as MGD, the refined level uses Dirichet boundary data at the coarse-fine interface and the data is derived from the coarse level solution. After advancing on the fine level, an additional procedure, the sync-solve (SS), is required in order to enforce conservation. The MGD SS reduces to an elliptic solve on a combined grid for a system of G equations, where G is the number of groups. We adapt the partial temperature scheme for the SS; hence, we reuse the infrastructure developed for scalar equations. Results are presented. (Abridged)