We present Atacama Large Millimeter Array CO(3$-$2) and HCO$^+$(4$-$3) observations covering the central $1\rlap{.}'5$$\times$$1\rlap{.}'5$ region of the Orion Nebula Cluster (ONC). The unprecedented level of sensitivity ($\sim$0.1 mJy beam$^{-1}$) and angular resolution ($\sim$$0\rlap{.}''09 \approx 35$ AU) of these line observations enable us to search for gas-disk detections towards the known positions of submillimeter-detected dust disks in this region. We detect 23 disks in gas: 17 in CO(3$-$2), 17 in HCO$^+$(4$-$3), and 11 in both lines. Depending on where the sources are located in the ONC, we see the line detections in emission, in absorption against the warm background, or in both emission and absorption. We spectrally resolve the gas with $0.5$ km s$^{-1}$ channels, and find that the kinematics of most sources are consistent with Keplerian rotation. We measure the distribution of gas-disk sizes and find typical radii of $\sim$50-200 AU. As such, gas disks in the ONC are compact in comparison with the gas disks seen in low-density star-forming regions. Gas sizes are universally larger than the dust sizes. However, the gas and dust sizes are not strongly correlated. We find a positive correlation between gas size and distance from the massive star $\theta^1$ Ori C, indicating that disks in the ONC are influenced by photoionization. Finally, we use the observed kinematics of the detected gas lines to model Keplerian rotation and infer the masses of the central pre-main-sequence stars. Our dynamically-derived stellar masses are not consistent with the spectroscopically-derived masses, and we discuss possible reasons for this discrepancy.

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.

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 present ALMA 850 $\mu$m continuum observations of the Orion Nebula Cluster that provide the highest angular resolution ($\sim 0\rlap{.}''1 \approx 40$ AU) and deepest sensitivity ($\sim 0.1$ mJy) of the region to date. We mosaicked a field containing $\sim 225$ optical or near-IR-identified young stars, $\sim 60$ of which are also optically-identified "proplyds". We detect continuum emission at 850 $\mu$m towards $\sim 80$% of the proplyd sample, and $\sim 50$% of the larger sample of previously-identified cluster members. Detected objects have fluxes of $\sim 0.5$-80 mJy. We remove sub-mm flux due to free-free emission in some objects, leaving a sample of sources detected in dust emission. Under standard assumptions of isothermal, optically thin disks, sub-mm fluxes correspond to dust masses of $\sim 0.5$ to 80 Earth masses. We measure the distribution of disk sizes, and find that disks in this region are particularly compact. Such compact disks are likely to be significantly optically thick. The distributions of sub-mm flux and inferred disk size indicate smaller, lower-flux disks than in lower-density star-forming regions of similar age. Measured disk flux is correlated weakly with stellar mass, contrary to studies in other star forming regions that found steeper correlations. We find a correlation between disk flux and distance from the massive star $\theta^1$ Ori C, suggesting that disk properties in this region are influenced strongly by the rich cluster environment.

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.