Semi-analytic models (SAMs) are a promising means of tracking the physical processes associated with galaxy formation, but many of their approximations have not been rigorously tested. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, we compare predictions from the FIRE-2 hydrodynamical "zoom-in" simulations to those from the Santa Cruz SAM run on the same halo merger trees, with an emphasis on the global mass flow cycle. Our study includes 13 halos spanning low-mass dwarfs (M_vir~10^10 M_sun at z=0), intermediate-mass dwarfs (M_vir~10^11 M_sun) and Milky Way-mass galaxies (M_vir~10^12 M_sun). The SAM and FIRE-2 predictions agree relatively well with each other in terms of stellar and interstellar mass, but differ dramatically on circumgalactic mass (the SAM is lower than FIRE-2 by ~3 orders of magnitude for dwarfs). Strikingly, the SAM predicts higher gas accretion rates for dwarfs compared to FIRE-2 by factors of ~10-100, and this is compensated for with higher mass outflow rates in the SAM. We argue that the most severe model discrepancies are caused by the lack of preventative stellar feedback and the assumptions for halo gas cooling and recycling in the SAM. As a first step towards resolving these model tensions, we present a simple yet promising new preventative stellar feedback model in which the energy carried by supernova-driven winds is allowed to heat some fraction of gas outside of halos to at least the virial temperature such that accretion is suppressed.

Intermediate and massive stars drive fast and powerful isotropic winds that interact with the winds of nearby stars in star clusters and the surrounding interstellar medium (ISM). Wind-ISM collisions generate astrospheres around these stars that contain hot $T\sim 10^7$ K gas that adiabatically expands. As individual bubbles expand and collide they become unstable, potentially driving turbulence in star clusters. In this paper we use hydrodynamic simulations to model a densely populated young star cluster within a homogeneous cloud to study stellar wind collisions with the surrounding ISM. We model a mass-segregated cluster of 20 B-type young main sequence stars with masses ranging from 3--17 $M_{\odot}$. We evolve the winds for $\sim$11 kyrs and show that wind-ISM collisions and over-lapping wind-blown bubbles around B-stars mixes the hot gas and ISM material generating Kolmogorov-like turbulence on small scales early in its evolution. We discuss how turbulence driven by stellar winds may impact the subsequent generation of star formation in the cluster

We examine the role of the interstellar magnetic field to modulate the orientation of turbulent flows within the Taurus molecular cloud using spatial gradients of thin velocity slices of 12CO and 13CO antenna temperatures. Our analysis accounts for the random errors of the gradients that arise from the thermal noise of the spectra. The orientations of the vectors normal to the antenna temperature gradient vectors are compared to the magnetic field orientations that are calculated from Planck 353~GHz polarization data. These relative orientations are parameterized with the projected Rayleigh statistic and mean resultant vector. For 12CO, 28 percent and 39 percent of the cloud area exhibit strongly parallel or strongly perpendicular relative orientations respectively. For the lower opacity 13CO emission, strongly parallel and strongly perpendicular orientations are found in 7 percent and 43 percent of the cloud area respectively. For both isotopologues, strongly parallel or perpendicular alignments are restricted to localized regions with low levels of turbulence. If the relative orientations serve as an observational proxy to the Alfvenic Mach number then our results imply local variations of the Alfvenic Mach number throughout the cloud.

Density fluctuations produced by supersonic turbulence are of great importance to astrophysical chemical models. A property of these density fluctuations is that the two point correlation function decreases with increasing scale separation. The relation between the density decorrelation length scale ($L_{\rm dec}$) and the turbulence driving scale ($L_{\rm drive}$) determines how turbulence affects the density and chemical structures in the interstellar medium (ISM), and is a key component for using observations of atomic and molecular tracers to constrain turbulence properties. We run a set of numerical simulations of supersonic magnetohydrodynamic turbulence, with different sonic Mach numbers ($\mathcal{M}_s=4.5, 7$) , and driven on varying scales (1/2.5, 1/5, 1/7) the box length. We derive the $L_{\rm dec}-L_{\rm drive}$ relation as a function of Mach number, driving scale, and the orientation of the line-of-sight (LOS) in respect to the magnetic-field. We find that the mean ratio $L_{\rm dec}/L_{\rm drive} = 0.19 \pm 0.10$, when averaged over snapshots, Mach numbers, driving lengths, and the three LOSs. For LOS parallel to the magnetic field the density structures are statistically smaller and the $L_{\rm dec}-L_{\rm drive}$ relation is tighter, with $L_{\rm dec}/L_{\rm drive} = 0.112 \pm 0.024$. We discuss our results in the context of using observations of chemical tracers to constrain the dominant turbulence driving scale in the ISM.

We have obtained a contiguous set of long-slit spectra of a shock wave in the Cygnus Loop to investigate its structure, which is far from the morphology predicted by 1D models. Proper motions from Hubble Space Telescope images combined with the known distance to the Cygnus Loop provide an accurate shock speed. Earlier analyses of shock spectra estimated the shock speed, postshock density, temperature, and elemental abundances. In this paper we determine several more shock parameters: a more accurate shock speed, ram pressure, density, compression ratio, dust destruction efficiency, magnetic field strength, and vorticity in the cooling region. From the derived shock properties we estimate the emissivities of synchrotron emission in the radio and pion decay emission in the gamma rays. Both are consistent with the observations if we assume simple adiabatic compression of ambient cosmic rays as in the van der Laan mechanism. We also find that, although the morphology is far from that predicted by 1D models and the line ratios vary dramatically from point to point, the average spectrum is matched reasonably well by 1D shock models with the shock speed derived from the measured proper motion. A subsequent paper will analyze the development of turbulence in the cooling zone behind the shock.

In this paper we present a detailed analysis of the Faraday depth (FD) spectrum and its clean components obtained through the application of the commonly used technique of Faraday rotation measure synthesis to analyze spectro-polarimetric data. In order to directly compare the Faraday depth spectrum with physical properties of a magneto-ionic medium, we generated synthetic broad-bandwidth spectro-polarimetric observations from magnetohydrodynamic (MHD) simulations of a transonic, isothermal, compressible turbulent medium. We find that correlated magnetic field structures give rise to a combination of spiky, localized peaks at certain FD values, and broad structures in the FD spectrum. Although the majority of these spiky FD structures appear narrow, giving an impression of a Faraday thin medium, we show that they arise from strong synchrotron emissivity at that FD. Strong emissivity at a FD can arise because of both strong spatially-local polarized synchrotron emissivity at a FD or accumulation of weaker emissions along the distance through a medium that have Faraday depths within half the width of the rotation measure spread function. Such a complex Faraday depth spectrum is a natural consequence of MHD turbulence when the lines of sight pass through a few turbulent cells. This therefore complicates the convention of attributing narrow FD peaks to presence of a Faraday rotating medium along the line of sight. Our work shows that it is difficult to extract the FD along a line of sight from the Faraday depth spectrum using standard methods for a turbulent medium in which synchrotron emission and Faraday rotation occur simultaneously.

Similar to their low-mass counterparts, massive stars likely form via the collapse of pre-stellar molecular cores. Recent observations suggest that most massive cores are subvirial (i.e., not supported by turbulence) and therefore are likely unstable to gravitational collapse. Here we perform radiation hydrodynamic simulations to follow the collapse of turbulent massive pre-stellar cores with subvirial and virialized initial conditions to explore how their dynamic state affects the formation of massive stars and core fragmentation into companion stars. We find that subvirial cores undergo rapid monolithic collapse resulting in higher accretion rates at early times as compared to the collapse of virialized cores that have the same physical properties. In contrast, we find that virialized cores undergo a slower, gradual collapse and significant turbulent fragmentation at early times resulting in numerous companion stars. In the absence of strong magnetic fields and protostellar outflows we find that the faster growth rate of massive stars that are born out of subvirial cores leads to an increase in the radiative heating of the core thereby further suppressing fragmentation at early times when turbulent fragmentation occurs for virialized cores. Regardless of initial condition, we find that the massive accretion disks that form around massive stars dominant the accretion flow onto the star at late times and eventually become gravitationally unstable and fragment to form companion stars at late times.

Supersonic turbulence results in strong density fluctuations in the interstellar medium (ISM), which have a profound effect on the chemical structure. Particularly useful probes of the diffuse ISM are the ArH$^+$, OH$^+$, H$_2$O$^+$ molecular ions, which are highly sensitive to fluctuations in the density and the H$_2$ abundance. We use isothermal magnetohydrodynamic (MHD) simulations of various sonic Mach numbers, $\mathcal{M}_s$, and density decorrelation scales, $y_{\rm dec}$, to model the turbulent density field. We post-process the simulations with chemical models and obtain the probability density functions (PDFs) for the H$_2$, ArH$^+$, OH$^+$ and H$_2$O$^+$ abundances. We find that the PDF dispersions increases with increasing $\mathcal{M}_s$ and $y_{\rm dec}$, as the magnitude of the density fluctuations increases, and as they become more coherent. Turbulence also affects the median abundances: when $\mathcal{M}_s$ and $y_{\rm dec}$ are high, low-density regions with low H$_2$ abundance become prevalent, resulting in an enhancement of ArH$^+$ compared to OH$^+$ and H$_2$O$^+$. We compare our models with Herschel observations. The large scatter in the observed abundances, as well as the high observed ArH$^+$/OH$^+$ and ArH$^+$/H$_2$O$^+$ ratios are naturally reproduced by our supersonic $(\mathcal{M}_s=4.5)$, large decorrelation scale $(y_{\rm dec}=0.8)$ model, supporting a scenario of a large-scale turbulence driving. The abundances also depend on the UV intensity, CR ionization rate, and the cloud column density, and the observed scatter may be influenced by fluctuations in these parameters.

Compressive isothermal turbulence is known to have a near lognormal density probability distribution function (PDF) with a width that scales with the sonic Mach number and nature of the turbulent driving (solenoidal vs compressive). However, the physical processes that mold the extreme high and low density structures in a turbulent medium can be different, with the densest structures being composed of strong shocks that evolve on shorter timescales than the low density fluid. The density PDF in a turbulent medium exhibits deviations from lognormal due to shocks, that increases with the sonic Mach number, which is often ignored in analytic models for turbulence and star formation. We develop a simple model for turbulence by treating it as a continuous Markov process, which explains both the density PDF and the transient timescales of structures as a function of density, using a framework developed in n Scannapieco & Safarzadeh (2018). Our analytic model depends on only a single parameter, the effective compressive sonic Mach number, and successfully describes the non-lognormal behavior seen in both 1D and 3D simulations of supersonic and subsonic compressive isothermal turbulence. The model quantifies the non-lognormal distribution of density structures in turbulent environments, and has application to star forming molecular clouds and star formation efficiencies.

The interstellar medium (ISM) of galaxies is composed of a turbulent magnetized plasma. In order to quantitatively measure relevant turbulent parameters of the ISM, a wide variety of statistical techniques and metrics have been developed that are often tested using numerical simulations and analytic formalism. These metrics are typically based on the Fourier power spectrum, which does not capture the Fourier phase information that carries the morphological characteristics of images. In this work we use density slices of magnetohydrodyanmic turbulence simulations to demonstrate that a modern tool, convolutional neural networks, can capture significant information encoded in the Fourier phases. We train the neural network to distinguish between two simulations with different levels of magnetization. We find that, even given a tiny slice of simulation data, a relatively simple network can distinguish sub-Alfv\'enic (strong magnetic field) and super-Alfv\'enic (weak magnetic field) turbulence >98% of the time, even when all spectral amplitude information is stripped from the images. In order to better understand how the neural network is picking out differences betweem the two classes of simulations we apply a neural network analysis method called "saliency maps". The saliency map analysis shows that sharp ridge-like features are a distinguishing morphological characteristic in such simulations. Our analysis provides a way forward for deeper understanding of the relationship between magnetohydrodyanmic turbulence and gas morphology and motivates further applications of neural networks for studies of turbulence. We make publicly available all data and software needed to reproduce our results.

The Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

The current disparity in computational knowledge is a critical hindrance to the diversity and success of the field. Recommendations are outlined for policies and funding models to enable the growth and retention of a new generation of computational researchers that reflect the demographics of the undergraduate population in Astronomy and Physics.

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

We explore how the velocity gradient technique (VGT) can be applied to absorbing media in the case of$^{13}$CO 2-1 emission. The VGT is a new way to trace magnetic fields in the plane of the sky using only spectroscopic observations. We apply the VGT to magnetohydrodynamic turbulence simulations that have been post-processed to include $^{13}$CO 2-1 emission and we calculate the velocity centroid gradients. We find that the velocity centroid gradients trace the projected magnetic field in media with different $^{13}$CO abundances, densities and optical depths. Our study opens up the possibility of using velocity centroid gradients to trace magnetic fields in molecular clouds using 13CO emission.

Precise measurement of the foreground synchrotron emission, which contaminates the faint polarized cosmic microwave background radiation (CMB), is a major challenge for the next-generation of CMB experiments. To address this, dedicated foreground measurement experiments are being undertaken at radio frequencies between 2 and 40 GHz. Foreground polarized synchrotron emission measurements are particularly challenging, primarily due to the complicated frequency dependence in the presence of Faraday rotation, and are best recovered through broad fractional-bandwidth polarization measurements at frequencies $\lesssim 5$ GHz. A unique opportunity for measuring the foreground polarized synchrotron emission will be provided by the 15-m SKA-MPG telescope operating in the frequency range 1.7 to 3.5~GHz (S-Band). Here, we present the scope of a Southern sky survey in S-Band at 1 degree angular resolution and explore its added advantage for application of powerful techniques, such as, Stokes $Q$, $U$ fitting and RM-synthesis. A full Southern-sky polarization survey with this telescope, when combined with other on-going efforts at slightly higher frequencies, will provide an excellent frequency coverage for modeling and extrapolating the foreground polarized synchrotron emission to CMB frequencies ($\gtrsim80$~GHz) with rms brightness temperature better than 10~nK per 1 degree$^2$. We find that this survey will be crucial for understanding the effects of Faraday depolarization, especially in low Galactic latitude regions. This will allow better foreground cleaning and thus will contribute significantly in further improving component separation analyses and increase usable sky area for cosmological analysis of the \textit{Planck} data, and the \textit{LiteBIRD} mission in the future.

We analytically calculate the star formation efficiency and dense self-gravitating gas fraction in the presence of magneto-gravo-turbulence using the model of Burkhart (2018), which employs a piecewise lognormal and powerlaw density Probability Distribution Function (PDF). We show that the PDF transition density from lognormal to powerlaw forms is a mathematically motivated critical density for star formation and can be physically related to the density where the Jeans length is comparable to the sonic length, i.e. the post-shock critical density for collapse. When the PDF transition density is taken as the critical density, the instantaneous star formation efficiency ($\epsilon_{\rm inst}$) and depletion time ($\tau_{\rm depl}$) can be calculated from the dense self-gravitating gas fraction represented as the fraction of gas in the PDF powerlaw tail. We minimize the number of free parameters in the analytic expressions for $\epsilon_{\rm inst}$ and $\tau_{\rm depl}$ by using the PDF transition density instead of a parameterized critical density for collapse and thus provide a more direct pathway for comparison with observations. We test the analytic predictions for the transition density and self-gravitating gas fraction against AREPO moving mesh gravoturbulent simulations and find good agreement. We predict that, when gravity dominates the density distribution in the star-forming gas, the star formation efficiency should be weakly anti-correlated with the sonic Mach number while the depletion time should increase with increasing sonic Mach number. The star formation efficiency and depletion time depend primarily on the dense self-gravitating gas fraction, which in turn depends on the interplay of gravity, turbulence and stellar feedback. Our model prediction is in agreement with recent observations, such as the M51 PdBI Arcsecond Whirlpool Survey (PAWS).

Supersonically induced gas objects (SIGOs) with little to no dark matter component are predicted to exist in patches of the Universe with non-negligible relative velocity between baryons and the dark matter at the time of recombination. Using {\sc arepo} hydrodynamic simulations we find that the gas densities inside these objects are high enough to allow stars to form. An estimate of the luminosity of the first star clusters formed within these SIGOs suggests that they may be observed at high redshift using future HST and JWST observations. Furthermore, our simulations indicate that SIGOs lie in a distinct place in the luminosity-radius parameter space, which can be used observationally to distinguish SIGOs from dark-matter hosting gas systems. Finally, as a proof-of-concept, we model star formation before reionization and evolve these systems to current times. We find that SIGOs occupy a similar part of the magnitude-radius parameter space as globular clusters. These results suggest that SIGOs may be linked with present-day metal-poor local globular clusters. Since the relative velocity between the baryons and dark matter is coherent over a few Mpc scales, we predict that if this is the dominant mechanism for the formation of globular clusters, their abundance should vary significantly over these scales.

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.

Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing.

The Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $\mu$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $\mu$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5\sigma)$, or will rule out with more than $5\sigma$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4\sigma$, constrain the effective number of light particle species with $\Delta N_{\rm eff}<0.06~(2\sigma)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $\sigma_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science.

We present the first application of the angle-dependent 3-Point Correlation Function (3PCF) to the density fields magnetohydrodynamic (MHD) turbulence simulations intended to model interstellar (ISM) turbulence. Previous work has demonstrated that the angle-averaged bispectrum, the 3PCF's Fourier-space analog, is sensitive to the sonic and Alfv\'enic Mach numbers of turbulence. Here we show that introducing angular information via multipole moments with respect to the triangle opening angle offers considerable additional discriminatory power on these parameters. We exploit a fast, order $N_{\rm g} \log N_{\rm g}$ ($N_{\rm g}$ the number of grid cells used for a Fourier Transform) 3PCF algorithm to study a suite of MHD turbulence simulations with 10 different combinations of sonic and Alfv\'enic Mach numbers over a range from sub to super-sonic and sub to super-Alfv\'{e}nic. The 3PCF algorithm's speed for the first time enables full quantification of the time-variation of our signal: we study 9 timeslices for each condition, demonstrating that the 3PCF is sufficiently time-stable to be used as an ISM diagnostic. In future, applying this framework to 3-D dust maps will enable better treatment of dust as a cosmological foreground as well as reveal conditions in the ISM that shape star formation.

Magnetohydrodynamic(MHD) turbulence displays velocity anisotropies which reflect the direction of the magnetic field. This anisotropy has led to the development of a number of statistical techniques for studying magnetic fields in the interstellar medium. In this paper, we review and compare three techniques that use radio position-position-velocity data for determining magnetic field strength and morphology : the correlation function anisotropy (CFA), Principal Component Analysis of Anisotropies (PCAA), and the more recent Velocity Gradient Technique (VGT). We compare these three techniques and suggest improvements to the CFA and PCAA techniques to increase their accuracy and versatility. In particular, we suggest and successfully implement a much faster way of calculating non-periodic correlation functions for the CFA. We discuss possible improvements to the current implementation of the PCAA. We show the advantages of the VGT in terms of magnetic field tracing and stress the complementary nature with the other two techniques.

Stars form in supersonic turbulent molecular clouds that are self-gravitating. We present an analytic determination of the star formation rate (SFR) in a gravoturbulent medium based on the density probability distribution function of molecular clouds having a piecewise lognormal and power law form. This is in contrast to previous analytic SFR models that are governed primarily by interstellar turbulence which sets purely lognormal density PDFs. In the gravoturbulent SFR model described herein, low density gas resides in the lognormal portion of the PDF. Gas becomes gravitationally unstable past a critical density ($\rho_{crit}$), and the PDF begins to forms a power law. As the collapse of the cloud proceeds, the transitional density ($\rho_t$) between the lognormal and power law portions of the PDF moves towards lower-density while the slope of the power law ($\alpha$) becomes increasingly shallow. The star formation rate per free-fall time is calculated via an integral over the lognormal from $\rho_{crit}$ to $\rho_t$ and an integral over the power law from $\rho_t$ to the maximum density. As $\alpha$ becomes shallower the SFR accelerates beyond the expected values calculated from a lognormal density PDF. We show that the star formation efficiency per free fall time in observations of local molecular cloud increases with shallower PDF power law slopes, in agreement with our model. Our model can explain why star formation is spatially and temporally variable within a cloud and why the depletion times observed in local and extragalactic giant molecular clouds vary. Both star-bursting and quiescent star-forming systems can be explained without the need to invoke extreme variations of turbulence in the local interstellar environment.

Supersonic isothermal turbulence establishes a network of transient dense shocks that sweep up material and have a density profile described by balance between ram pressure of the background fluid versus the magnetic and gas pressure gradient behind the shock. These rare, densest regions of a turbulent environment can become Jeans unstable and collapse to form pre-stellar cores. Using numerical simulations of magneto-gravo-turbulence, we describe the structural properties of dense shocks, which are the seeds of gravitational collapse, as a function of magnetic field strength. In the regime of a weak magnetic field, the collapse is isotropic. Strong magnetic field strengths lead to significant anisotropy in the shocked distribution and collapse occurs preferentially parallel to the field lines. Our work provides insight into analyzing the magnetic field topology and density structures of young protostellar collapse, which the theory presented here predicts are associated with large-scale strong shocks that persist for at least a free-fall time.

We present the first results from a new, high resolution, $^{12}$CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) molecular line survey of the Orion A cloud, hereafter referred to as the CARMA-NRO Orion Survey. CARMA observations have been combined with single-dish data from the Nobeyama 45m telescope to provide extended images at about 0.01 pc resolution, with a dynamic range of approximately 1200 in spatial scale. Here we describe the practical details of the data combination in uv space, including flux scale matching, the conversion of single dish data to visibilities, and joint deconvolution of single dish and interferometric data. A $\Delta$-variance analysis indicates that no artifacts are caused by combining data from the two instruments. Initial analysis of the data cubes, including moment maps, average spectra, channel maps, position-velocity diagrams, excitation temperature, column density, and line ratio maps provides evidence of complex and interesting structures such as filaments, bipolar outflows, shells, bubbles, and photo-eroded pillars. The implications for star formation processes are profound and follow-up scientific studies by the CARMA-NRO Orion team are now underway. We plan to make all the data products described here generally accessible; some are already available at https://dataverse.harvard.edu/dataverse/CARMA-NRO-Orion

We introduce a new model for the structure and evolution of the gas in galactic discs. In the model the gas is in vertical pressure and energy balance. Star formation feedback injects energy and momentum, and non-axisymmetric torques prevent the gas from becoming more than marginally gravitationally unstable. From these assumptions we derive the relationship between galaxies' bulk properties (gas surface density, stellar content, and rotation curve) and their star formation rates, gas velocity dispersions, and rates of radial inflow. We show that the turbulence in discs can be powered primarily by star formation feedback, radial transport, or a combination of the two. In contrast to models that omit either radial transport or star formation feedback, the predictions of this model yield excellent agreement with a wide range of observations, including the star formation law measured in both spatially resolved and unresolved data, the correlation between galaxies' star formation rates and velocity dispersions, and observed rates of radial inflow. The agreement holds across a wide range of galaxy mass and type, from local dwarfs to extreme starbursts to high-redshifts discs. We apply the model to galaxies on the star-forming main sequence, and show that it predicts a transition from mostly gravity-driven turbulence at high redshift to star formation-driven turbulence at low redshift. This transition, and the changes in mass transport rates that it produces, naturally explain why galaxy bulges tend to form at high redshift and discs at lower redshift, and why galaxies tend to quench inside-out.

We derive two-dimensional spatial power spectra of four distinct interstellar medium tracers, HI, $^{12}$CO($J$=1--0), $^{13}$CO($J$=1--0), and dust, in the Perseus molecular cloud, covering linear scales ranging from $\sim$0.1 pc to $\sim$90 pc. Among the four tracers, we find the steepest slopes of $-3.23\pm0.05$ and $-3.22\pm0.05$ for the uncorrected and opacity-corrected HI column density images. This result suggests that the HI in and around Perseus traces a non-gravitating, transonic medium on average, with a negligible effect from opacity. On the other hand, we measure the shallowest slope of $-2.72\pm0.12$ for the 2MASS dust extinction data and interpret this as the signature of a self-gravitating, supersonic medium. Possible variations in the dust-to-gas ratio likely do not change our conclusion. Finally, we derive slopes of $-3.08\pm0.08$ and $-2.88\pm0.07$ for the ${}^{12}$CO(1--0) and ${}^{13}$CO(1--0) integrated intensity images. Based on theoretical predictions for an optically thick medium, we interpret these slopes of roughly $-3$ as implying that both CO lines are susceptible to the opacity effect. While simple tests for the impact of CO formation and depletion indicate that the measured slopes of ${}^{12}$CO(1--0) and ${}^{13}$CO(1--0) are not likely affected by these chemical effects, our results generally suggest that chemically more complex and/or fully optically thick media may not be a reliable observational tracer for characterizing turbulence.

Diagnostics of polarized emission provide us with valuable information on the Galactic magnetic field and the state of turbulence in the interstellar medium, which cannot be obtained from synchrotron intensity alone. In Paper I (Herron et al. 2017b), we derived polarization diagnostics that are rotationally and translationally invariant in the $Q$-$U$ plane, similar to the polarization gradient. In this paper, we apply these diagnostics to simulations of ideal magnetohydrodynamic turbulence that have a range of sonic and Alfv\'enic Mach numbers. We generate synthetic images of Stokes $Q$ and $U$ for these simulations, for the cases where the turbulence is illuminated from behind by uniform polarized emission, and where the polarized emission originates from within the turbulent volume. From these simulated images we calculate the polarization diagnostics derived in Paper I, for different lines of sight relative to the mean magnetic field, and for a range of frequencies. For all of our simulations, we find that the polarization gradient is very similar to the generalized polarization gradient, and that both trace spatial variations in the magnetoionic medium for the case where emission originates within the turbulent volume, provided that the medium is not supersonic. We propose a method for distinguishing the cases of emission coming from behind or within a turbulent, Faraday rotating medium, and a method to partly map the rotation measure of the observed region. We also speculate on statistics of these diagnostics that may allow us to constrain the physical properties of an observed turbulent region.

The column density probability distribution function (N-PDF) of GMCs has been used as a diagnostic of star formation. Simulations and analytic predictions have suggested the N-PDF is composed of a low density lognormal component and a high density power-law component, tracing turbulence and gravitational collapse, respectively. In this paper, we study how various properties of the true 2D column density distribution create the shape, or "anatomy" of the PDF. We test our ideas and analytic approaches using both a real, observed, PDF based on Herschel observations of dust emission as well as a simulation that uses the ENZO code. Using a dendrogram analysis, we examine the three main components of the N-PDF: the lognormal component, the power-law component, and the transition point between these two components. We find that the power-law component of an N-PDF is the summation of N-PDFs of power-law substructures identified by the dendrogram algorithm. We also find that the analytic solution to the transition point between lognormal and power-law components proposed by Burkhart, Stalpes & Collins (2017) is applicable when tested on observations and simulations, within the uncertainties. We reconfirm and extend the results of Lombardi, Alves & Lada (2015), which stated that the lognormal component of the N-PDF is difficult to constrain due to the artificial choice of the map area. Based on the resulting anatomy of the N-PDF, we suggest avoiding analyzing the column density structures of a star forming region based solely on fits to the lognormal component of an N-PDF. We also suggest applying the N-PDF analysis in combination with the dendrogram algorithm, to obtain a more complete picture of the global and local environments and their effects on the density structures.

Magnetically active stars possess stellar winds whose interaction with planetary magnetic fields produces radio auroral emission. We examine the detectability of radio auroral emission from Proxima b, the closest known exosolar planet orbiting our nearest neighboring star, Proxima Centauri. Using the Radiometric Bode's Law, we estimate the radio flux produced by the interaction of Proxima Centauri's stellar wind and Proxima b's magnetosphere for different planetary magnetic field strengths. For plausible planetary masses, Proxima b produces 6-83 mJy of auroral radio flux at frequencies of 0.3-0.8 MHz for planetary magnetic field strengths of 1-3 B$_{\oplus}$. According to recent MHD models that vary the orbital parameters of the system, this emission is expected to be highly variable. This variability is due to large fluctuations in the size of Proxima b's magnetosphere as it crosses the equatorial streamer regions of the dense stellar wind and high dynamic pressure. Using the MHD model of Garraffo et al. 2016 for the variation of the magnetosphere radius during the orbit, we estimate that the observed radio flux can vary nearly by an order of magnitude over the 11.2 day period of Proxima b. The detailed amplitude variation depends on the stellar wind, orbital, and planetary magnetic field parameters. We discuss observing strategies for proposed future space-based observatories to reach frequencies below the ionospheric cut off ($\sim 10$ MHz) as would be required to detect the signal we investigate.

We report ALMA observations of polarized dust emission from the protostellar source Ser-emb 8 at a linear resolution of 140 au. Assuming models of dust-grain alignment hold, the observed polarization pattern gives a projected view of the magnetic field structure in this source. Contrary to expectations based on models of strongly magnetized star formation, the magnetic field in Ser-emb 8 does not exhibit an hourglass morphology. Combining the new ALMA data with previous observational studies, we can connect magnetic field structure from protostellar core (~80,000 au) to disk (~100 au) scales. We compare our observations with four magnetohydrodynamic gravo-turbulence simulations made with the AREPO code that have initial conditions ranging from super-Alfv\'enic (weakly magnetized) to sub-Alfv\'enic (strongly magnetized). These simulations achieve the spatial dynamic range necessary to resolve the collapse of protostars from the parsec scale of star-forming clouds down to the ~100 au scale probed by ALMA. Only in the very strongly magnetized simulation do we see both the preservation of the field direction from cloud to disk scales and an hourglass-shaped field at < 1000 au scales. We conduct an analysis of the relative orientation of the magnetic field and the density structure in both the Ser-emb 8 ALMA observations and the synthetic observations of the four AREPO simulations. We conclude that the Ser-emb 8 data are most similar to the weakly magnetized simulations, which exhibit random alignment, in contrast to the strongly magnetized simulation, where the magnetic field plays a role in shaping the density structure in the source. In the weak-field case, it is turbulence -- not the magnetic field -- that shapes the material that forms the protostar, highlighting the dominant role that turbulence can play across many orders of magnitude in spatial scale.

We study the effect of density fluctuations induced by turbulence on the HI/H$_2$ structure in photodissociation regions (PDRs) both analytically and numerically. We perform magnetohydrodynamic numerical simulations for both subsonic and supersonic turbulent gas, and chemical HI/H$_2$ balance calculations. We derive atomic-to-molecular density profiles and the HI column density probability density function (PDF) assuming chemical equilibrium. We find that while the HI/H$_2$ density profiles are strongly perturbed in turbulent gas, the mean HI column density is well approximated by the uniform-density analytic formula of Sternberg et al. (2014). The PDF width depends on (a) the radiation intensity to mean density ratio, (b) the sonic Mach number and (c) the turbulence decorrelation scale, or driving scale. We derive an analytic model for the HI PDF and demonstrate how our model, combined with 21 cm observations, can be used to constrain the Mach number and driving scale of turbulent gas. As an example, we apply our model to observations of HI in the Perseus molecular cloud. We show that a narrow observed HI PDF may imply small scale decorrelation, pointing to the potential importance of subcloud-scale turbulence driving.

Star formation in our Galaxy occurs in molecular clouds that are self-gravitating, highly turbulent, and magnetized. We study the conditions under which cloud cores inherit large-scale magnetic field morphologies and how the field is governed by cloud turbulence. We present four moving-mesh simulations of supersonic, turbulent, isothermal, self-gravitating gas with a range of magnetic mean-field strengths characterized by the Alfv\'enic Mach number $\mathcal{M}_{{\rm A}, 0}$, resolving pre-stellar core formation from parsec to a few AU scales. In our simulations with the turbulent kinetic energy density dominating over magnetic pressure ($\mathcal{M}_{{\rm A}, 0}>1$), we find that the collapse is approximately isotropic with $B\propto\rho^{2/3}$, core properties are similar regardless of initial mean-field strength, and the field direction on $100$ AU scales is uncorrelated with the mean field. However, in the case of a dominant large-scale magnetic field ($\mathcal{M}_{{\rm A}, 0}=0.35$), the collapse is anisotropic with $B\propto\rho^{1/2}$. This transition at $\mathcal{M}_{{\rm A}, 0}\sim1$ is not expected to be sharp, but clearly signifies two different paths for magnetic field evolution in star formation. Based on observations of different star forming regions, we conclude that star formation in the interstellar medium may occur in both regimes. Magnetic field correlation with the mean-field extends to smaller scales as $\mathcal{M}_{{\rm A}, 0}$ decreases, making future ALMA observations useful for constraining $\mathcal{M}_{{\rm A}, 0}$ of the interstellar medium.

We investigate the effects of AGN heating and the ultraviolet background on the low-redshift Lyman-$\alpha$ forest column density distribution (CDD) using the Illustris simulation. We show that Illustris reproduces observations at $z =0.1$ in the column density range $10^{12.5} - 10^{13.5}$cm$^{-2}$, relevant for the "photon underproduction crisis." We attribute this to the inclusion of AGN feedback, which changes the gas distribution so as to mimic the effect of extra photons, as well as the use of the Faucher-Gigu\`ere ultraviolet background, which is more ionizing at $z=0.1$ than the Haardt & Madau background previously considered. We show that the difference between simulations run with smoothed particle hydrodynamics and simulations using a moving mesh is small in this column density range but can be more significant at larger column densities. We further consider the effect of supernova feedback, Voigt profile fitting and finite resolution, all of which we show to have little influence on the CDD. Finally, we identify a discrepancy between our simulations and observations at column densities $10^{14} - 10^{16}$cm$^{-2}$, where Illustris produces too few absorbers, which suggests the AGN feedback model should be further refined. Since the "photon underproduction crisis" primarily affects lower column density systems, we conclude that AGN feedback and standard ionizing background models can resolve the crisis.

The Breakthrough Starshot initiative aims to launch a gram-scale spacecraft to a speed of $v\sim 0.2$c, capable of reaching the nearest star system, $\alpha$ Centauri, in about 20 years. However, a critical challenge for the initiative is the damage to the spacecraft by interstellar gas and dust during the journey. In this paper, we quantify the interaction of a relativistic spacecraft with gas and dust in the interstellar medium. For gas bombardment, we find that damage by track formation due to heavy elements is an important effect. We find that gas bombardment can potentially damage the surface of the spacecraft to a depth of $\sim 0.1$ mm for quartz material after traversing a gas column of $N_{\rm H}\sim 2\times 10^{18}\rm cm^{-2}$ along the path to $\alpha$ Centauri, whereas the effect is much weaker for graphite material. The effect of dust bombardment erodes the spacecraft surface and produces numerous craters due to explosive evaporation of surface atoms. For a spacecraft speed $v=0.2c$, we find that dust bombardment can erode a surface layer of $\sim 0.5$ mm thickness after the spacecraft has swept a column density of $N_{\rm H}\sim 3\times 10^{17}\rm cm^{-2}$, assuming the standard gas-to-dust ratio of the interstellar medium. Dust bombardment also damages the spacecraft surface by modifying the material structure through melting. We calculate the equilibrium surface temperature due to collisional heating by gas atoms as well as the temperature profile as a function of depth into the spacecraft. Our quantitative results suggest methods for damage control, and we highlight possibilities for shielding strategies and protection of the spacecraft.

Previous studies have shown that star formation depends on the driving of molecular cloud turbulence, and differences in the driving can produce an order of magnitude difference in the star formation rate. The turbulent driving is characterised by the parameter $\zeta$, with $\zeta=0$ for compressive, curl-free driving (e.g. accretion or supernova explosions), and $\zeta=1$ for solenoidal, divergence-free driving (e.g. Galactic shear). Here we develop a new method to measure $\zeta$ from observations of synchrotron emission from molecular clouds. We calculate statistics of mock synchrotron intensity images produced from magnetohydrodynamic simulations of molecular clouds, in which the driving was controlled to produce different values of $\zeta$. We find that the mean and standard deviation of the log-normalised synchrotron intensity are sensitive to $\zeta$, for values of $\zeta$ between $0$ (curl-free driving) and $0.5$ (naturally-mixed driving). We quantify the dependence of zeta on the direction of the magnetic field relative to the line of sight. We provide best-fit formulae for $\zeta$ in terms of the log-normalised mean and standard deviation of synchrotron intensity, with which $\zeta$ can be determined for molecular clouds that have similar Alfv\'enic Mach number to our simulations. These formulae are independent of the sonic Mach number. Signal-to-noise ratios larger than $5$, and angular resolutions smaller than $5\%$ of the cloud diameter, are required to apply these formulae. Although there are no firm detections of synchrotron emission from molecular clouds, by combining Green Bank Telescope and Very Large Array observations it should be possible to detect synchrotron emission from molecular clouds, thereby constraining the value of $\zeta$.

We characterize the column density probability distributions functions (PDFs) of the atomic hydrogen gas, HI, associated with seven Galactic molecular clouds (MCs). We use 21 cm observations from the Leiden/Argentine/Bonn Galactic HI Survey to derive column density maps and PDFs. We find that the peaks of the HI PDFs occur at column densities ranging from ~1-2$\times 10^{21}$ cm$^2$ (equivalently, ~0.5-1 mag). The PDFs are uniformly narrow, with a mean dispersion of $\sigma_{HI}\approx 10^{20}$ cm$^2$ (~0.1 mag). We also investigate the HI-to-H$_2$ transition towards the cloud complexes and estimate HI surface densities ranging from 7-16 $M_\odot$ pc$^{-2}$ at the transition. We propose that the HI PDF is a fitting tool for identifying the HI-to-H$_2$ transition column in Galactic MCs.

We derive an analytic expression for the transitional column density value ($s_t$) between the lognormal and power-law form of the probability distribution function (PDF) in star-forming molecular clouds. Our expression for $s_t$ depends on the mean column density, the variance of the lognormal portion of the PDF, and the slope of the power-law portion of the PDF. We show that $s_t$ can be related to physical quantities such as the sonic Mach number of the flow and the power-law index for a self-gravitating isothermal sphere. This implies that the transition point between the lognormal and power-law density/column density PDF represents the critical density where turbulent and thermal pressure balance, the so-called "post-shock density." We test our analytic prediction for the transition column density using dust PDF observations reported in the literature as well as numerical MHD simulations of self-gravitating supersonic turbulence with the Enzo code. We find excellent agreement between the analytic $s_t$ and the measured values from the numerical simulations and observations (to within 1.5 A$_V$). We discuss the utility of our expression for determining the properties of the PDF from unresolved low density material in dust observations, for estimating the post-shock density, and for determining the HI-H$_2$ transition in clouds.

Using data from the ALFALFA Arecibo HI survey of galaxies and the Virgo cluster X-ray pressure profiles from XMM-Newton, we investigate the possibility that starless dark HI clumps, also known as "dark galaxies", are supported by external pressure in the surrounding intercluster medium. We find that the starless HI clump masses, velocity dispersions and positions allow these clumps to be in pressure equilibrium with the X-ray gas near the virial radius of the Virgo cluster. We predict the sizes of these clumps to range from 1-10 kpc, in agreement with the range of sizes found for spatially resolved HI starless clumps outside of Virgo. Based on the predicted HI surface density of the Virgo sources, as well as a sample of other similar resolved ALFALFA HI dark clumps with follow up optical/radio observations, we predict that most of the HI dark clumps are on the cusp of forming stars. These HI sources therefore mark the transition between starless HI clouds and dwarf galaxies with stars.

We investigate how observations of synchrotron intensity fluctuations can be used to probe the sonic and Alfv\'enic Mach numbers of interstellar turbulence, based on mock observations performed on simulations of magnetohydrodynamic turbulence. We find that the structure function slope, and a diagnostic of anisotropy that we call the integrated quadrupole ratio modulus, both depend on the Alfv\'enic Mach number. However, these statistics also depend on the orientation of the mean magnetic field in the synchrotron emitting region relative to our line of sight, and this creates a degeneracy that cannot be broken by observations of synchrotron intensity alone. We conclude that the polarization of synchrotron emission could be analyzed to break this degeneracy, and suggest that this will be possible with the Square Kilometre Array.

Galaxies' interstellar media (ISM) are observed to be supersonically-turbulent, but the ultimate power source that drives turbulent motion remains uncertain. The two dominant models are that the turbulence is driven by star formation feedback and/or that it is produced by gravitational instability in the gas. Here we show that, while both models predict that the galaxies' ISM velocity dispersions will be positively correlated with their star formation rates, the forms of the correlation predicted by these two models are subtly but measurably different. A feedback-driven origin for the turbulence predicts a velocity dispersion that rises more sharply with star formation rate, and that does not depend on the gas fraction (i.e. $\dot{M}_* \propto \sigma^2$), while a gravity-driven model yields a shallower rise and a strong dependence on gas fraction (i.e. $\dot{M}_* \propto f_g^2 \sigma$). We compare the models to a collection of data on local and high-redshift galaxies culled from the literature, and show that the correlation expected for gravity-driven turbulence is a better match to the observations than a feedback-driven model. This suggests that gravity is the ultimate source of ISM turbulence, at least in the rapidly-star-forming, high velocity dispersion galaxies for which our test is most effective. We conclude by discussing the limitations of the present data set, and the prospects for future measurements to enable a more definitive test of the two models.

Studies of MHD turbulence often investigate the Fourier power spectrum to provide information on the nature of the turbulence cascade. However, the Fourier power spectrum only contains the Fourier amplitudes and rejects all information regarding the Fourier phases. Here we investigate the utility of two statistical diagnostics for recovering information on Fourier phases in ISM density data: the averaged amplitudes of the bispectrum and the phase coherence index (PCI), a new phase technique for the ISM. We create 3D density and 2D column density maps using a set of simulations of isothermal ideal MHD turbulence with a wide range of sonic and Alfv\'enic Mach numbers. We find that the bispectrum averaged along different angles with respect to either the $k_1$ or $k_2$ axis is primarily sensitive to the sonic Mach number while averaging the bispectral amplitudes over different annuli is sensitive to both the sonic and Alfv\'enic Mach numbers. The PCI of density suggests that the most correlated phases occur in supersonic sub-Alfv\'enic turbulence and also near the numerical dissipation regime. This suggests that non-linear interactions with correlated phases are strongest in shock dominated regions, in agreement with findings from the solar wind. Additionally, our results are particularly encouraging as they suggests the phase information contained in the bispectrum and PCI can be used to find parameters of turbulence in column density maps.

In this work we investigate the Principal Component Analysis (PCA) sensitivity to the velocity power spectrum in high opacity regimes of the interstellar medium (ISM). For our analysis we use synthetic Position-Position-Velocity (PPV) cubes of fractional Brownian motion (fBm) and magnetohydrodynamics (MHD) simulations, post processed to include radiative transfer effects from CO. We find that PCA analysis is very different from the tools based on the traditional power spectrum of PPV data cubes. Our major finding is that PCA is also sensitive to the phase information of PPV cubes and this allows PCA to detect the changes of the underlying velocity and density spectra at high opacities, where the spectral analysis of the maps provides the universal -3 spectrum in accordance with the predictions of Lazarian \& Pogosyan (2004) theory. This makes PCA potentially a valuable tool for studies of turbulence at high opacities provided that the proper gauging of the PCA index is made. The later, however, we found to be not easy, as the PCA results change in an irregular way for data with high sonic Mach numbers. This is in contrast to synthetic Brownian noise data used for velocity and density fields that show monotonic PCA behavior. We attribute this difference to the PCA's sensitivity to Fourier phase information.

The shape of the probability distribution function (PDF) of molecular clouds is an important ingredient for modern theories of star formation and turbulence. Recently, several studies have pointed out observational difficulties with constraining the low column density (i.e. Av <1) PDF using dust tracers. In order to constrain the shape and properties of the low column density probability distribution function, we investigate the PDF of multiphase atomic gas in the Perseus molecular cloud using opacity-corrected GALFA-HI data and compare the PDF shape and properties to the total gas PDF and the N(H2) PDF. We find that the shape of the PDF in the atomic medium of Perseus is well described by a lognormal distribution, and not by a power-law or bimodal distribution. The peak of the atomic gas PDF in and around Perseus lies at the HI-H2 transition column density for this cloud, past which the N(H2) PDF takes on a powerlaw form. We find that the PDF of the atomic gas is narrow and at column densities larger than the HI-H2 transition the HI rapidly depletes, suggesting that the HI PDF may be used to find the HI-H2 transition column density. We also calculate the sonic Mach number of the atomic gas by using HI absorption line data, which yields a median value of Ms=4.0 for the CNM, while the HI emission PDF, which traces both the WNM and CNM, has a width more consistent with transonic turbulence.

We present the results of the Velocity Coordinate Spectrum (VCS) technique to calculate the velocity power spectrum of turbulence in the Small Magellanic Cloud (SMC) in 21cm emission. We have obtained a velocity spectral index of -3.85 and an injection scale of 2.3 kpc. The spectral index is steeper than the Kolmogorov index which is expected for shock-dominated turbulence which is in agreement with past works on the SMC gas dynamics. The injection scale of 2.3 kpc suggests that tidal interactions with the Large Magellanic Cloud are the dominate driver of turbulence in this dwarf galaxy. This implies turbulence maybe driven by multiple mechanisms in galaxies in addition to supernova injection and that galaxy-galaxy interactions may play an important role.

We present a catalog of 1964 isolated, compact neutral hydrogen clouds from the Galactic Arecibo L-Band Feed Array Survey Data Release One (GALFA-HI DR1). The clouds were identified by a custom machine-vision algorithm utilizing Difference of Gaussian kernels to search for clouds smaller than 20'. The clouds have velocities typically between |VLSR| = 20-400 km/s, linewidths of 2.5-35 km/s, and column densities ranging from 1 - 35 x 10^18 cm^-2. The distances to the clouds in this catalog may cover several orders of magnitude, so the masses may range from less than a Solar mass for clouds within the Galactic disc, to greater than 10^4 Solar Masses for HVCs at the tip of the Magellanic Stream. To search for trends, we separate the catalog into five populations based on position, velocity, and linewidth: high velocity clouds (HVCs); galaxy candidates; cold low velocity clouds (LVCs); warm, low positive-velocity clouds in the third Galactic Quadrant; and the remaining warm LVCs. The observed HVCs are found to be associated with previously-identified HVC complexes. We do not observe a large population of isolated clouds at high velocities as some models predict. We see evidence for distinct histories at low velocities in detecting populations of clouds corotating with the Galactic disc and a set of clouds that is not corotating.

We study the observable signatures of self-gravitating MHD turbulence by applying the probability density functions (PDFs) and the spatial density power spectrum to synthetic column density maps. We find that there exists three characterizable stages of the evolution of the collapsing cloud which we term "early," "intermediate," and "advanced." At early times, i.e. $t<0.15t_{ff}$, the column density has a power spectral slope similar to nongravitating supersonic turbulence and a lognormal distribution. At an intermediate stage, i.e. $0.15t_{ff}< t \leq 0.35t_{ff}$, there exists signatures of the prestellar cores in the shallower PDF and power spectrum power law slopes. The column density PDF power law tails at these times have line of sight averaged slopes ranging from -2.5 to -1.5 with shallower values belonging to simulations with lower magnetic field strength. The density power spectrum slope becomes shallow and can be characterized by $P(k)=A_1k^{\beta_2}e^{-k/k_c}$, where $A_1$ describes the amplitude, $k^{\beta_2}$ describes the classical power law behavior and the scale $k_c$ characterizes the turn over from turbulence dominated to self-gravity dominated. At advanced stages of collapse, i.e. $\approx t>0.35t_{ff}$, the power spectral slope is positive valued, and a dramatic increase is observed in the PDF moments and the Tsallis incremental PDF parameters, which gives rise to deviations between PDF-sonic Mach number relations. Finally, we show that the imprint of gravity on the density power spectrum can be replicated in non-gravitating turbulence by introducing a delta-function with amplitude equivalent to the maximum valued point in a given self-gravitating map. We find that the turbulence power spectrum restored through spatial filtering of the high density material.

We investigate the nature of the Alfv\'enic turbulence cascade in two fluid MHD simulations in order to determine if turbulence is damped once the ion and neutral species become decoupled at a critical scale called the ambipolar diffusion scale (L$_{AD}$). Using mode decomposition to separate the three classical MHD modes, we study the second order structure functions of the Alfv\'en mode velocity field of both neutrals and ions in the reference frame of the local magnetic field. On scales greater than L$_{AD}$ we confirm that two fluid turbulence strongly resembles single fluid MHD turbulence. Our simulations show that the behavior of two fluid turbulence becomes more complex on scales less than L$_{AD}$. We find that Alfvenic turbulence can exist past L$_{AD}$ when the turbulence is globally super-Alfv\'enic, with the ions and neutrals forming separate cascades once decoupling has taken place. When turbulence is globally sub-Alfvenic and hence strongly anisotropic with a large separation between the parallel and perpendicular decoupling scales, turbulence is damped at L$_{AD}$. We also find that the power spectrum of the kinetic energy in the damped regime is consistent with a $k^{-4}$ scaling (in agreement with the predictions of Lazarian, Vishniac & Cho 2004).

The dynamics of the interstellar medium (ISM) are strongly affected by turbulence, which shows increased anisotropy in the presence of a magnetic field. We expand upon the Esquivel & Lazarian method to estimate the Alfven Mach number using the structure function anisotropy in velocity centroid data from position-position-velocity maps. We utilize 3D magnetohydrodynamic (MHD) simulations of fully developed turbulence, with a large range of sonic and Alfvenic Mach numbers, to produce synthetic observations of velocity centroids with observational characteristics such as thermal broadening, cloud boundaries, noise, and radiative transfer effects of carbon monoxide. In addition, we investigate how the resulting anisotropy-Alfven Mach number dependency found in Esquivel & Lazarian (2011) might change when taking the second moment of the position-position-velocity cube or when using different expressions to calculate the velocity centroids. We find that the degree of anisotropy is related primarily to the magnetic field strength (i.e. Alfven Mach number) and the line-of-sight orientation, with a secondary effect on sonic Mach number. If the line-of-sight is parallel to up to ~45 deg off of the mean field direction, the velocity centroid anisotropy is not prominent enough to distinguish different Alfvenic regimes. The observed anisotropy is not strongly affected by including radiative transfer, although future studies should include additional tests for opacity effects. These results open up the possibility of studying the magnetic nature of the ISM using statistical methods in addition to existing observational techniques.

Radio synchrotron polarization maps of the Galaxy can be used to infer the properties of interstellar turbulence in the diffuse warm ionized medium (WIM). In this paper, we investigate the spatial gradient of linearly polarized synchrotron emission ($|\nabla\textbf{P}|/|\textbf{P}|$) as a tracer of turbulence, the relationship of the gradient to the sonic Mach number of the WIM, and changes in morphology of the gradient as a function of Galactic position in the southern sky. We use data from the S-band Polarization All Sky Survey (S-PASS) to image the spatial gradient of the linearly polarized synchrotron emission ($|\nabla \textbf{P}|/|\textbf{P}|$) of the entire southern sky at $2.3$~GHz. The spatial gradient of linear polarization reveals rapid changes of the density and magnetic fluctuations in the WIM due to magnetic turbulence as a function of Galactic position; we make comparisons of these data to ideal MHD numerical simulations. In order to constrain the sonic Mach number ($M_{s}$), we apply a high order moments analysis to the observations and to the simulated diffuse, isothermal ISM with ideal magneto-hydrodynamic turbulence. We find that polarization gradient maps reveal elongated structures, which we associate with turbulence in the ISM. Our analysis corroborates the view of a turbulent WIM in a transonic regime $M_{s}\lesssim 2$. Filamentary structures with typical widths down to the angular resolution are seen and the observed morphologies match closely with numerical simulations and in some cases H$\alpha$ contours. The $|\nabla \textbf{P}|/|\textbf{P}|$ intensity is found to be approximately log-normal distributed. No systematic variations of the sonic Mach number are observed as a function of Galactic coordinates, which is consistent with turbulence in the WIM inferred by the analysis of H$\alpha$ data.

We study how the estimation of the sonic Mach number ($M_s$) from $^{13}$CO linewidths relates to the actual 3D sonic Mach number. For this purpose we analyze MHD simulations which include post-processing to take radiative transfer effects into account. As expected, we find very good agreement between the linewidth estimated sonic Mach number and the actual sonic Mach number of the simulations for optically thin tracers. However, we find that opacity broadening causes $M_s$ to be overestimated by a factor of ~ 1.16 - 1.3 when calculated from optically thick $^{13}$CO lines. We also find that there is a dependency on the magnetic field: super-Alfv\'enic turbulence shows increased line broadening as compared with sub-Alfv\'enic turbulence for all values of optical depth for supersonic turbulence. Our results have implications for the observationally derived sonic Mach number--density standard deviation ($\sigma_{\rho/<\rho>}$) relationship, $\sigma^2_{\rho/<\rho>}=b^2M_s^2$, and the related column density standard deviation ($\sigma_{N/<N>}$) sonic Mach number relationship. In particular, we find that the parameter b, as an indicator of solenoidal vs. compressive driving, will be underestimated as a result of opacity broadening. We compare the $\sigma_{N/<N>}-M_s$ relation derived from synthetic dust extinction maps and $^{13}$CO linewidths with recent observational studies and find that solenoidally driven MHD turbulence simulations have values of $\sigma_{N/<N>}$ which are lower than real molecular clouds. This may be due to the influence of self-gravity which should be included in simulations of molecular cloud dynamics.

We investigate turbulent properties of the non-star-forming, translucent molecular cloud, MBM16 by applying the statistical technique of a two-dimensional spatial power spectrum (SPS) on the neutral hydrogen (HI) observations obtained by the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) survey. The SPS, calculated over the range of spatial scales from 0.1 to 17 pc, is well represented with a single power-law function, with a slope ranging from -3.3 to -3.7 and being consistent over the velocity range of MBM16 for a fixed velocity channel thickness. However, the slope varies significantly with the velocity slice thickness, suggesting that both velocity and density contribute to HI intensity fluctuations. By using this variation we estimate the slope of 3D density fluctuations in MBM16 to be -3.7\pm0.2. This is significantly steeper than what has been found for HI in the Milky Way plane, the Small Magellanic Cloud, or the Magellanic Bridge, suggesting that interstellar turbulence in MBM16 is driven on scales >17 pc and that the lack of stellar feedback could be responsible for the steep power spectrum.

We present high resolution simulations of two-fluid (ion-neutral) MHD turbulence with resolutions as large as 512^3. The simulations are supersonic and mildly sub-Alfvenic, in keeping with the conditions present in molecular clouds. Such turbulence is thought to influence star formation processes in molecular clouds because typical cores form on length scales that are comparable to the dissipation scales of this turbulence in the ions. The simulations are motivated by the fact that recent studies of isophotologue lines in molecular clouds have found significant differences in the linewidth-size relationship for neutral and ion species. The goals of this paper are to explain those observations using simulations and analytic theory, present a new set of density-based diagnostics by drawing on similar diagnostics that have been obtained by studying single-fluid turbulence, and show that our two-fluid simulations play a vital role in reconciling alternative models of star formation. The velocity-dependent diagnostics display a complementarity with the density-dependent diagnostics. We find that the linewidth-size relationships show a prominent difference between ions and neutrals when the line of sight is orthogonal to the mean field. This is because the MHD waves in the ions differ from hydrodynamic waves in the neutrals. We also find that the density probability distribution functions (PDFs) show prominent differences between the ions and neutrals when the line of sight is parallel to the mean field. This is because the velocity fluctuations in the ions tend to produce column density fluctuations along field lines. When the magnetic field makes an angle to the line of sight, both observable differences should be visible. These diagnostics should be easy for observers to test. This analysis assumes optically thin lines and a mean magnetic field that is uniform in direction within a cloud.

The Fourier power spectrum is one of the most widely used statistical tools to analyze the nature of magnetohydrodynamic turbulence in the interstellar medium. Lazarian & Pogosyan (2004) predicted that the spectral slope should saturate to -3 for an optically thick medium and many observations exist in support of their prediction. However, there have not been any numerical studies to-date testing these results. We analyze the spatial power spectrum of MHD simulations with a wide range of sonic and Alfv\'enic Mach numbers, which include radiative transfer effects of the $^{13}$CO transition. We confirm numerically the predictions of Lazarian & Pogosyan (2004) that the spectral slope of line intensity maps of an optically thick medium saturates to -3. Furthermore, for very optically thin supersonic CO gas, where the density or CO abundance values are too low to excite emission in all but the densest shock compressed gas, we find that the spectral slope is shallower than expected from the column density. Finally, we find that mixed optically thin/thick CO gas, which has average optical depths on order of unity, shows mixed behavior: for super-Alfv\'enic turbulence, the integrated intensity power spectral slopes generally follow the same trend with sonic Mach number as the true column density power spectrum slopes. However, for sub-Alfv\'enic turbulence the spectral slopes are steeper with values near -3 which are similar to the very optically thick regime.

Magnetohydrodynamic turbulence is able to create hierarchical structures in the interstellar medium that are correlated on a wide range of scales via the energy cascade. We use hierarchical tree diagrams known as dendrograms to characterize structures in synthetic Position-Position-Velocity (PPV) emission cubes of optically thin isothermal magnetohydrodynamic turbulence. We show that the structures and degree of hierarchy observed in PPV space are related to the physics of the gas, i.e. self-gravity and the global sonic and Alfvenic Mach number. Simulations with higher Alfvenic Mach number, self-gravity and supersonic flows display enhanced hierarchical structure. We observed a strong sonic and Alfvenic dependency when we apply the the statistical moments (i.e. mean, variance, skewness, kurtosis) to the dendrogram distribution. Larger magnetic field and sonic Mach number correspond to larger values of the moments. Application of the dendrogram to 3D density cubes, also known as Position-Position-Position cubes (PPP), reveals that the dominant emission contours in PPP and PPV are related for supersonic gas but not for subsonic. We also explore the effects of smoothing, thermal broadening and velocity resolution on the dendrograms in order to make our study more applicable to observational data. These results all point to hierarchical tree diagrams as being a promising additional tool for studying ISM turbulence and star forming regions in the direction of obtaining information on the degree of self-gravity, the Mach numbers and the complicated relationship between PPV and PPP.

We study the effects of radiative transfer on the Probability Distribution Functions (PDFs) of simulations of magnetohydrodynamic turbulence in the widely studied $^{13}$CO 2-1 transition. We find that the integrated intensity maps generally follow a log-normal distribution, with the cases that have $\tau \approx 1$ best matching the PDF of the column density. We fit a 2D variance-sonic Mach number relationship to our logarithmic PDFs of the form $\sigma_{ln(\Sigma/\Sigma_0)}^2=A\times ln(1+b^2{\cal M}_s^2)$ and find that, for parameter $b=1/3$, parameter $A$ depends on the radiative transfer environment. We also explore the variance, skewness, and kurtosis of the linear PDFs finding that higher moments reflect both higher sonic Mach number and lower optical depth. Finally, we apply the Tsallis incremental PDF function and find that the fit parameters depend on both Mach numbers, but also are sensitive to the radiative transfer parameter space, with the $\tau \approx 1$ case best fitting the incremental PDF of the true column density. We conclude that, for PDFs of low optical depth cases, part of the gas is always sub-thermally excited so that the spread of the line intensities exceeds the spread of the underlying column densities and hence the PDFs do not reflect the true column density. Similarly, PDFs of optically thick cases are dominated by the velocity dispersion and therefore do not represent the true column density PDF. Thus, in the case of molecules like carbon monoxide, the dynamic range of intensities, structures observed and consequently, the observable PDFs, are less determined by turbulence and more-often determined by radiative transfer effects.

Although there are a wealth of column density tracers for both the molecular and diffuse interstellar medium, there are very few observational studies investigating the relationship between the density variance ($\sigma^2$) and the sonic Mach number (${\cal M}_s$). This is in part due to the fact that the $\sigma^2$-${\cal M}_s$ relationship is derived, via MHD simulations, for the 3D density variance only, which is not a direct observable. We investigate the utility of a 2D column density $\sigma_{\Sigma/\Sigma_0}^2$-${\cal M}_s$ relationship using solenoidally driven isothermal MHD simulations and find that the best fit follows closely the form of the 3D density $\sigma_{\rho/\rho_0}^2$-${\cal M}_s$ trend but includes a scaling parameter $A$ such that: $\sigma_{ln(\Sigma/\Sigma_0)}^2=A\times ln(1+b^2{\cal M}_s^2)$, where $A=0.11$ and $b=1/3$. This relation is consistent with the observational data reported for the Taurus and IC 5146 molecular clouds with $b=0.5$ and $A=0.16$ and $b=0.5$ and $A=0.12$, respectively. These results open up the possibility of using the 2D column density values of $\sigma^2$ for investigations of the relation between the sonic Mach number and the PDF variance in addition to existing PDF sonic Mach number relations.

Faraday rotation of linearly polarized radio signals provides a very sensitive probe of fluctuations in the interstellar magnetic field and ionized gas density resulting from magnetohydrodynamic (MHD) turbulence. We used a set of statistical tools to analyze images of the spatial gradient of linearly polarized radio emission ($|\nabla \textbf{P}|$) from the ISM for both observational data from a test image of the Southern Galactic Plane Survey (SGPS) and isothermal simulations of MHD turbulence. We compared the observational data with results of synthetic observations obtained with the simulations of 3D turbulence. Visually, in both data sets, a complex network of filamentary structures is seen. Our analysis shows that the filaments in the gradient can be produced by shocks as well as random fluctuations characterizing the non-differentiable field of MHD turbulence. The latter dominates for subsonic turbulence, while the former dominates for supersonic turbulence. In order to quantitatively characterize these differences we use the topology tool known as a genus curve as well as the moments of the image distribution. We find that higher values for the moments correspond to cases of $|\nabla \textbf{P}|$ with larger Mach numbers, but the strength of the dependency is connected to the telescope angular resolution. In regards to the topology, the supersonic filaments observed in $|\nabla \textbf{P}|$ have a positive genus shift, which indicates a "swisscheese" like topology, while the subsonic cases show a negative genus, indicating a "clump" like topology. In the case of the genus, the dependency on the telescope resolution is not as strong. The SGPS test region data has a distribution and morphology that matches subsonic to transsonic type turbulence, which independently confirms what is now expected for the WIM.

The interstellar medium of the Milky Way is multi-phase, magnetized and turbulent. Turbulence in the interstellar medium produces a global cascade of random gas motions, spanning scales ranging from 100 parsecs to 1000 kilometres. Fundamental parameters of interstellar turbulence such as the sonic Mach number (the speed of sound) have been difficult to determine because observations have lacked the sensitivity and resolution to directly image the small-scale structure associated with turbulent motion. Observations of linear polarization and Faraday rotation in radio emission from the Milky Way have identified unusual polarized structures that often have no counterparts in the total radiation intensity or at other wavelengths, and whose physical significance has been unclear. Here we report that the gradient of the Stokes vector (Q,U), where Q and U are parameters describing the polarization state of radiation, provides an image of magnetized turbulence in diffuse ionized gas, manifested as a complex filamentary web of discontinuities in gas density and magnetic field. Through comparison with simulations, we demonstrate that turbulence in the warm ionized medium has a relatively low sonic Mach number, M_s <~ 2. The development of statistical tools for the analysis of polarization gradients will allow accurate determinations of the Mach number, Reynolds number and magnetic field strength in interstellar turbulence over a wide range of conditions.

MHD Turbulence is a critical component of the current paradigms of star formation, particle transport, magnetic reconnection and evolution of the ISM. Progress on this difficult subject is made via numerical simulations and observational studies. However, due to limitations of resolution, scale discrepancies, and complexity of the observations, the best approach for connecting numerics to observations is not always obvious. Here we advocate for a approach that invokes statistical techniques to understand the underlying physics of turbulent astrophysical systems. The wealth of numerical and observational data calls for new statistical tools to be developed in order to study turbulence in the interstellar medium. We briefly review some of the recently developed statistics that focus on characterizing gas compressibility and magnetization and their uses to interstellar studies.

In an effort to characterize the Mach numbers of ISM magnetohydrodynamic (MHD) turbulence, we study the probability distribution functions (PDFs) of patial increments of density, velocity, and magnetic field for fourteen ideal isothermal MHD simulations at resolution 512^3. In particular, we fit the PDFs using the Tsallis function and study the dependency of fit parameters on the compressibility and magnetization of the gas. We find that the Tsallis function fits PDFs of MHD turbulence well, with fit parameters showing sensitivities to the sonic and Alfven Mach numbers. For 3D density, column density, and position-position-velocity (PPV) data we find that the amplitude and width of the PDFs shows a dependency on the sonic Mach number. We also find the width of the PDF is sensitive to global Alfvenic Mach number especially in cases where the sonic number is high. These dependencies are also found for mock observational cases, where cloud-like boundary conditions, smoothing, and noise are introduced. The ability of Tsallis statistics to characterize sonic and Alfvenic Mach numbers of simulated ISM turbulence point to it being a useful tool in the analysis of the observed ISM, especially when used simultaneously with other statistical techniques.

MHD Turbulence is a critical component of the current paradigms of star formation, particle transport, magnetic reconnection and evolution of the ISM, to name just a few. Progress on this difficult subject is made via numerical simulations and observational studies, however in order to connect these two, statistical methods are required. This calls for new statistical tools to be developed in order to study turbulence in the interstellar medium. Here we briefly review some of the recently developed statistics that focus on characterizing gas compressibility and magnetization and their uses to interstellar studies.

We investigate the nature and spatial variations of turbulence in the Small Magellanic Cloud (SMC) by applying several statistical methods on the neutral hydrogen (HI) column density image of the SMC and a database of isothermal numerical simulations. By using the 3rd and 4th statistical moments we derive the spatial distribution of the sonic Mach number (M_s) across the SMC. We find that about 90% of the HI in the SMC is subsonic or transonic. However, edges of the SMC `bar' have M_s=4 and may be tracing shearing or turbulent flows. Using numerical simulations we also investigate how the slope of the spatial power spectrum depends on both sonic and Alfven Mach numbers. This allows us to gauge the Alfven Mach number of the SMC and conclude that its gas pressure dominates over the magnetic pressure. The super-Alfvenic nature of the HI gas in the SMC is also highlighted by the bispectrum, a three-point correlation function which characterizes the level of non-Gaussianity in wave modes. We find that the bispectrum of the SMC HI column density displays similar large-scale correlations as numerical simulations, however it has localized enhancements of correlations. In addition, we find a break in correlations at a scale of 160 pc. This may be caused by numerous expanding shells of a similar size.

We present a number of statistical tools for obtaining studying turbulence in molecular clouds and diffuse interstellar medium. For our tests we used of three-dimensional 512 cube compressible MHD isothermal simulations performed for different sonic and Alfvenic Mach numbers. We introduce the bispectrum, a new tool for statistical studies of the interstellar medium which, unlike an ordinary power spectrum of turbulence, preserves the phase information of the stochastic field. We show that the bispectra of the 3D stochastic density field and of column densities, available from observations, are similar. We use the bispectrum technique to define the role of non-linear wave-wave interactions in the turbulent energy cascade. We also obtained the bispectrum function for density and column densities with varying magnetic field strength. Larger values of sonic Mach number result in increased correlations for modes with different wavenumbers. This effect becomes more evident with increasing magnetic field intensity. In addition to the bispectrum, we calculated the 3rd and 4th statistical moments of density and column density, namely, skewness and kurtosis, respectively. We found a strong dependence of skewness and kurtosis with the sonic Mach number. In particular, as this number increases, so does the asymmetry of the density distribution. We also studied the correlations of 2D column density with dispersion of velocities and magnetic field, as well as the correlations of 3D density with magentic and kinetic energy and Alfven Mach number for comparison. Our results show that column density is linearly correlated with magnetic field for high sonic Mach number. This trend is independent of the turbulent kinetic energy and can be used to characterize inhomogeneities of physical properties in low density clumps in the ISM.