The density structure of the interstellar medium (ISM) determines where stars form and release energy, momentum, and heavy elements, driving galaxy evolution. Density variations are seeded and amplified by gas motion, but the exact nature of this motion is unknown across spatial scale and galactic environment. Although dense star-forming gas likely emerges from a combination of instabilities, convergent flows, and turbulence, establishing the precise origin is challenging because it requires quantifying gas motion over many orders of magnitude in spatial scale. Here we measure the motion of molecular gas in the Milky Way and in nearby galaxy NGC 4321, assembling observations that span an unprecedented spatial dynamic range ($10^{-1}{-}10^3$ pc). We detect ubiquitous velocity fluctuations across all spatial scales and galactic environments. Statistical analysis of these fluctuations indicates how star-forming gas is assembled. We discover oscillatory gas flows with wavelengths ranging from $0.3{-}400$ pc. These flows are coupled to regularly-spaced density enhancements that likely form via gravitational instabilities. We also identify stochastic and scale-free velocity and density fluctuations, consistent with the structure generated in turbulent flows. Our results demonstrate that ISM structure cannot be considered in isolation. Instead, its formation and evolution is controlled by nested, interdependent flows of matter covering many orders of magnitude in spatial scale.

The dynamical production of low-mass X-ray binaries and brighter cataclysmic variables (CVs) in dense globular clusters is well-established. We investigate how the X-ray emissivity of fainter X-ray binaries (principally CVs and coronally active binaries) varies between different environments. We compile calculations (largely from the literature) of the X-ray emissivity of old stellar populations, including open and globular clusters and several galaxies. We investigate three literature claims of unusual X-ray sources in low-density stellar populations. We show that a suggested quiescent neutron star in the open cluster NGC 6819 is a foreground M dwarf. We show that the suggested diffuse X-ray emission from an old nova shell in the globular cluster NGC 6366 is actually a background galaxy cluster. And we show that a suggested population of quiescent X-ray binaries in the Sculptor Dwarf Galaxy is mostly (perhaps entirely) background galaxies. We find that above densities of $10^4$ M$_{\odot}$/pc$^3$, the X-ray emissivity of globular clusters increases, due to dynamical production of X-ray emitting systems. Below this density, globular clusters have lower X-ray emissivity than the other populations, and we do not see a strong dependence of X-ray emissivity due to density effects. We find significant correlations between X-ray emissivity and binary fraction, metallicity, and density. Sampling these fits via bootstrap techniques gives less significant correlations, but confirms the effect of metallicity on low-density populations, and that of density on the full globular cluster sample.

We analyze the 1D spatial power spectra of dust surface density and mid to far-infrared emission at $24-500\,\mu$m in the LMC, SMC, M31, and M33. By forward-modelling the point-spread-function (PSF) on the power spectrum, we find that nearly all power spectra have a single power-law and point source component. A broken power-law model is only favoured for the LMC 24 $\mu$m MIPS power spectrum and is due to intense dust heating in 30 Doradus. We also test for local power spectrum variations by splitting the LMC and SMC maps into $820$ pc boxes. We find significant variations in the power-law index with no strong evidence for breaks. The lack of a ubiquitous break suggests that the spatial power spectrum does not constrain the disc scale height. This contradicts claims of a break where the turbulent motion changes from 3D to 2D. The power spectrum indices in the LMC, SMC, and M31 are similar ($2.0-2.5$). M33 has a flatter power spectrum ($1.3$), similar to more distant spiral galaxies with a centrally-concentrated H$_{2}$ distribution. We compare the power spectra of HI, CO, and dust in M31 and M33 and find that HI power spectra are consistently flatter than CO power spectra. These results cast doubt on the idea that the spatial power spectrum traces large scale turbulent motion in nearby galaxies. Instead, we find that the spatial power spectrum is influenced by (1) the PSF on scales below $\sim3$ times the FWHM, (2) bright compact regions (30 Doradus), and (3) the global morphology of the tracer (an exponential CO disc).

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 investigate how the spectral properties of atomic (HI) and molecular (H$_2$) gas, traced by CO(2-1), are related in M33 on $80$ pc scales. We find the HI and CO(2-1) velocity at peak intensity to be highly correlated, consistent with previous studies. By stacking spectra aligned to the velocity of HI peak intensity, we find that the CO line width ($\sigma_{\rm HWHM}=4.6\pm0.9$ km s$^{-1}$; $\sigma_{\rm HWHM}$ is the effective Gaussian width) is consistently smaller than the HI line width ($\sigma_{\rm HWHM}=6.6\pm0.1$ km s$^{-1}$), with a ratio of ${\sim}0.7$, in agreement with Druard et al. (2014). The ratio of the line widths remains less than unity when the data are smoothed to a coarser spatial resolution. In other nearby galaxies, this line width ratio is close to unity which has been used as evidence for a thick, diffuse molecular disk that is distinct from the thin molecular disk dominated by molecular clouds. The smaller line width ratio found here suggests that M33 has a marginal thick molecular disk. From modelling individual lines-of-sight, we recover a strong correlation between HI and CO line widths when only the HI located closest to the CO component is considered. The median line width ratio of the line-of-sight line widths is $0.56\pm0.01$. There is substantial scatter in the HI--CO(2-1) line width relation, larger than the uncertainties, that results from regional variations on $<500$ pc scales, and there is no significant trend in the line widths, or their ratios, with galactocentric radius. These regional line width variations may be a useful probe of changes in the local cloud environment or the evolutionary state of molecular clouds.

We present new L-band ($1\mbox{--}2$ GHz) observations of the nearby spiral galaxy M33 with 80 pc resolution obtained with the Karl G. Jansky Very Large Array. The HI observations, combined with HI measurements from the Green Bank Telescope, improve the spectral resolution and sensitivity ($2.8$ K rms noise in a $0.2$ km s$^{-1}$ channel) compared to previous observations. We find individual profiles are usually non-Gaussian, harbouring line wings, multiple components, and asymmetries. Given this spectral complexity, we quantify the motions in the atomic ISM through moment analysis of the spectra and fits to aligned, stacked profiles. The measured value of the HI line width depends strongly on the method used, with the velocity stacked profiles aligned to the peak velocity giving the minimum value of $\sigma = 7$ km s$^{-1}$ and all other methods giving higher values ($\sigma\sim10$ km s$^{-1}$). All measurements of the line width show a shallow radial trend, with $\sigma$ decreasing by $\sim2$ km s$^{-1}$ from $R_{\rm gal}=0$ to $R_{\rm gal}=8$ kpc. We consider a number of energy sources that might maintain the line width against turbulent dissipation, but no single source is adequate. We find excess emission relative to a Gaussian in the stacked profile line wings, ranging from 9% to 26% depending on how the spectra are aligned. By splitting the line wings into symmetric and asymmetric components, we find that the lagging rotational disk accounts for one-third of the line wing flux. We also find emission far from the rotation axis of the galaxy in multiple discrete HI clouds, including a filament with a projected length of $\sim8$ kpc.

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 report one OH 1665-MHz detection in new observations with the NSF's Karl G. Jansky Very Large Array. To the authors' knowledge, this is the first OH maser detection in M33. The source occurs at a position of $\alpha=23.50069(9)^{\circ}$ and $\delta = 30.67984(9)^{\circ}$ (J2000) and is unresolved at the beam size ($9."7$). Two H$\alpha$ sources from the Moody et al. (2017) catalogue --- C1-1 and C1-2 --- are consistent within the beam size. We find a peak flux density of 24.1 mJy in a 1.5 km s$^{-1}$ channel. At the distance of M33, the source's specific luminosity is consistent with the brightest Galactic OH 1665-MHz masers.

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 have modified the iterative procedure introduced by Lin et al. (2016), to systematically combine the submm images taken from ground based (e.g., CSO, JCMT, APEX) and space (e.g., Herschel, Planck) telescopes. We applied the updated procedure to observations of three well studied Infrared Dark Clouds (IRDCs): G11.11-0.12, G14.225-0.506 and G28.34+0.06, and then performed single-component, modified black-body fits to derive $\sim$10$"$ resolution dust temperature and column density maps. The derived column density maps show that these three IRDCs exhibit complex filamentary structures embedding with rich clumps/cores. We compared the column density probability distribution functions (N-PDFs) and two-point correlation (2PT) functions of the column density field between these IRDCs with several OB cluster-forming regions. Based on the observed correlation and measurements, and complementary hydrodynamical simulations for a 10$^{4}$ $\rm M_{\odot}$ molecular cloud, we hypothesize that cloud evolution can be better characterized by the evolution of the (column) density distribution function and the relative power of dense structures as a function of spatial scales, rather than merely based on the presence of star-forming activity. Based on the small analyzed sample, we propose four evolutionary stages, namely: {\it cloud integration, stellar assembly, cloud pre-dispersal and dispersed-cloud.} The initial {\it cloud integration} stage and the final {\it dispersed cloud} stage may be distinguished from the two intermediate stages by a steeper than $-$4 power-law index of the N-PDF. The {\it cloud integration} stage and the subsequent {\it stellar assembly} stage are further distinguished from each other by the larger luminosity-to-mass ratio ($>$40 $\rm L_{\odot}/M_{\odot}$) of the latter.

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 have developed an iterative procedure to systematically combine the millimeter and submillimeter images of OB cluster-forming molecular clouds, which were taken by ground based (CSO, JCMT, APEX, IRAM-30m) and space telescopes (Herschel, Planck). For the seven luminous ($L$$>$10$^{6}$ $L_{\odot}$) Galactic OB cluster-forming molecular clouds selected for our analyses, namely W49A, W43-Main, W43-South, W33, G10.6-0.4, G10.2-0.3, G10.3-0.1, we have performed single-component, modified black-body fits to each pixel of the combined (sub)millimeter images, and the Herschel PACS and SPIRE images at shorter wavelengths. The $\sim$10$"$ resolution dust column density and temperature maps of these sources revealed dramatically different morphologies, indicating very different modes of OB cluster-formation, or parent molecular cloud structures in different evolutionary stages. The molecular clouds W49A, W33, and G10.6-0.4 show centrally concentrated massive molecular clumps that are connected with approximately radially orientated molecular gas filaments. The W43-Main and W43-South molecular cloud complexes, which are located at the intersection of the Galactic near 3-kpc (or Scutum) arm and the Galactic bar, show a widely scattered distribution of dense molecular clumps/cores over the observed $\sim$10 pc spatial scale. The relatively evolved sources G10.2-0.3 and G10.3-0.1 appear to be affected by stellar feedback, and show a complicated cloud morphology embedded with abundant dense molecular clumps/cores. We find that with the high angular resolution we achieved, our visual classification of cloud morphology can be linked to the systematically derived statistical quantities (i.e., the enclosed mass profile, the column density probability distribution function, the two-point correlation function of column density, and the probability distribution function of clump/core separations).

We present a new algorithm for detecting filamentary structure FilFinder. The algorithm uses the techniques of mathematical morphology for filament identification, presenting a complementary approach to current algorithms which use matched filtering or critical manifolds. Unlike other methods, FilFinder identifies filaments over a wide dynamic range in brightness. We apply the new algorithm to far infrared imaging data of dust emission released by the Herschel Gould Belt Survey team. Our preliminary analysis characterizes both filaments and fainter striations. We find a typical filament width of 0.09 pc across the sample, but the brightness varies from cloud to cloud. Several regions show a bimodal filament brightness distribution, with the bright mode (filaments) being an order of magnitude brighter than the faint mode (striations). Using the Rolling Hough Transform, we characterize the orientations of the striations in the data, finding preferred directions that agree with magnetic field direction where data are available. There is a suggestive but noisy correlation between typical filament brightness and literature values of the star formation rates for clouds in the Gould Belt.

We report observations using the Swift/XRT, NuSTAR, and Chandra X-ray telescopes of the transient X-ray source CXOGC J174540.0-290005, during its 2013 outburst. Due to its location in the field of multiple observing campaigns targeting Sgr A*, this is one of the best-studied outbursts of a very faint X-ray binary (VFXB; peak $L_X<10^{36}$ erg/s) yet recorded, with detections in 173 ks of X-ray observations over 50 days. VFXBs are of particular interest, due to their unusually low outburst luminosities and time-averaged mass transfer rates, which are hard to explain within standard accretion physics and binary evolution. The 2013 outburst of CXOGC J174540.0-290005 peaked at $L_X$(2-10 keV)=$5.0\times10^{35}$ erg/s, and all data above $10^{34}$ ergs/s were well-fit by an absorbed power-law of photon index $\sim1.7$, extending from 2 keV out to $\sim$70 keV. We discuss the implications of these observations for the accretion state of CXOGC J174540.0-290005.

The FourStar infrared camera is a 1.0-2.5 micron (JHKs) near infrared camera for the Magellan Baade 6.5m telescope at Las Campanas Observatory (Chile). It is being built by Carnegie Observatories and the Instrument Development Group at Johns Hopkins and is scheduled for completion in 2009. The instrument uses four Teledyne HAWAII-2RG arrays that produce a 10.9 x 10.9 arcmin field of view. The outstanding seeing at the Las Campanas site coupled with FourStar's high sensitivity and large field of view will enable many new survey and targeted science programs.

PANIC (Persson's Auxiliary Nasmyth Infrared Camera) is a near-infrared camera designed to operate at any one of the f/11 folded ports of the 6.5m Magellan telescopes at Las Campanas Observatory, Chile. The instrument is built around a simple, all-refractive design that reimages the Magellan focal plane to a plate scale of 0.125'' pixel^{-1} onto a Rockwell 1024x1024 HgCdTe detector. The design goals for PANIC included excellent image quality to sample the superb seeing measured with the Magellan telescopes, high throughput, a relatively short construction time, and low cost. PANIC has now been in regular operation for over one year and has proved to be highly reliable and produce excellent images. The best recorded image quality has been ~0.2'' FWHM.