Physical properties of the interstellar medium (ISM) at sub-galactic ($\sim$kpc) scales play an indispensable role in controlling the ability of gas to form stars. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, in this paper, we use the TNG50 cosmological simulation of the IllustrisTNG project to explore the physical parameter space comprised of 8 resolved ISM properties in star-forming regions across a wide range of host galaxy mass (M$_\star = 10^{7-11}$ M$_\odot$) and redshift ($0 \leq z \leq 3$). We deconstruct our simulated galaxies into regions of 1 kpc in size, and measure the gas/stellar surface densities, gas metallicity, vertical stellar velocity dispersion, epicyclic frequency and dark-matter volumetric density representative of each region, and study them in the context of their star formation activity and galactic environment (radial galactocentric location). By examining the star formation rate weighted distributions of these properties, we show that stars primarily form in two distinct regimes, which are brought about by an underlying bi-component radial star formation rate profile in galaxies. We examine how the relative prominence of these two regimes depends on host galaxy mass and cosmic time, and compare our findings with observations from the MaNGA IFU survey. Further, using principal component analysis, we characterise the aforementioned parameter space to reveal a high-degree of multicollinearity in relationships between ISM properties that drive the distribution of star formation at kiloparsec scales. Based on this, we find that a reduced three dimensional representation underpinned in essence by a multi-variate radius relationship is sufficient to capture most of the variance in the original 8D space.

A cool phase of the interstellar medium has been observed in many giant elliptical galaxies, but its origin remains unclear. We propose that uneven heating from Type Ia supernovae (SNe Ia), together with radiative cooling, can lead to the formation of the cool phase. The basic idea is that since SNe Ia explode randomly, gas parcels which are not directly heated by SN shocks will cool, forming multiphase gas. We run a series of idealized high-resolution numerical simulations, and find that cool gas develops even when the overall SNe heating rate $H$ exceeds the cooling rate $C$ by a factor as large as 1.4. We also find that the time for multiphase gas development depends on the gas temperature. When the medium has a temperature $T = 3\times 10^6$ K, the cool phase forms within one cooling time \tc; however, the cool phase formation is delayed to a few times \tc\ for higher temperatures. The main reason for the delay is turbulent mixing. Cool gas formed this way would naturally have a metallicity lower than that of the hot medium. For constant $H/C$, there is more turbulent mixing for higher temperature gas. We note that this mechanism of producing cool gas cannot be captured in cosmological simulations, which usually fail to resolve individual SN remnants.

Type Ia supernovae (SNe Ia) provide unique and important feedback in quiescent galaxies, but their impact has been underappreciated. In this paper, we analyze a series of high-resolution simulations to examine the energetics and turbulence of the medium under SNe Ia. We find that when SN remnants are resolved, their effects differ distinctly from a volumetric heating term, as is commonly assumed in unresolved simulations. First, the net heating is significantly higher than expected, by 30$\pm$10\% per cooling time. This is because a large fraction of the medium is pushed into lower densities which cool inefficiently. Second, the medium is turbulent; the root-mean-squared (RMS) velocity of the gas to 20-50 km s$^{-1}$ on a driving scale of tens of parsec. The velocity field of the medium is dominated by compressional modes, which are larger than the solenoidal components by a factor of 3-7. Third, the hot gas has a very broad density distribution. The ratio between the density fluctuations and the RMS Mach number, parameterized as $b$, is 2-20. This is in contrast to previous simulations of turbulent media, which have found $b\lesssim$ 1. The reason for the difference is mainly caused by the \textit{localized} heating of SNe Ia, which creates a large density contrast. Last, the typical length scale of a density fluctuation grows with time, forming increasingly larger bubbles and filamentary ridges. These underlying density fluctuations need to be included when X-ray observations are interpreted.

Large scale outflows in star-forming galaxies are observed to be ubiquitous, and are a key aspect of theoretical modeling of galactic evolution in a cosmological context, the focus of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project. Gas blown out from galactic disks, similar to gas within galaxies, consists of multiple phases with large contrasts of density, temperature, and other properties. To study multiphase outflows as emergent phenomena, we run a suite of ~pc-resolution local galactic disk simulations using the TIGRESS framework. Explicit modeling of the interstellar medium (ISM), including star formation and self-consistent radiative heating plus supernova feedback, regulates ISM properties and drives the outflow. We investigate the scaling of outflow mass, momentum, energy, and metal loading factors with galactic disk properties, including star formation rate (SFR) surface density (\Sigma_SFR~10^{-4}-1 M_sun/kpc^2/yr), gas surface density (~1-100 M_sun/pc^2), and total midplane pressure (or weight) (~10^3-10^6 k_B cm^{-3} K). The main components of outflowing gas are mass-delivering cool gas (T~10^4 K) and energy/metal-delivering hot gas (T~10^6 K). Cool mass outflow rates measured at outflow launch points (one or two scale heights) are 1-100 times the SFR (decreasing with \Sigma_SFR), although in massive galaxies most mass falls back due to insufficient outflow velocity. The hot galactic outflow carries mass comparable to 10% of the SFR, together with 10-20% of the energy and 30-60% of the metal mass injected by SN feedback. The characteristic outflow velocities of both phases scale very weakly with SFR, as v_out \propto \Sigma_SFR^{0.1~0.2}, consistent with observations. Importantly, our analysis demonstrates that in any physically-motivated cosmological wind model, it is crucial to include at least two distinct thermal wind components.

It has been hypothesized that photons from young, massive star clusters are responsible for maintaining the ionization of diffuse warm ionized gas seen in both the Milky Way and other disk galaxies. For a theoretical investigation of the warm ionized medium (WIM), it is crucial to solve radiation transfer equations where the ISM and clusters are modeled self-consistently. To this end, we employ a Solar neighborhood model of TIGRESS, a magnetohydrodynamic simulation of the multiphase, star-forming ISM, and post-process the simulation with an adaptive ray tracing method to transfer UV radiation from star clusters. We find that the WIM volume filling factor is highly variable, and sensitive to the rate of ionizing photon production and ISM structure. The mean WIM volume filling factor rises to ~0.15 at |z|~1 kpc. Approximately half of ionizing photons are absorbed by gas and half by dust; the cumulative ionizing photon escape fraction is 1.1%. Our time-averaged synthetic H$\alpha$ line profile matches WHAM observations on the redshifted (outflowing) side, but has insufficient intensity on the blueshifted side. Our simulation matches the Dickey-Lockman neutral density profile well, but only a small fraction of snapshots have high-altitude WIM density consistent with Reynolds Layer estimates. We compute a clumping correction factor C = <n_e>/sqrt<n_e^2>~0.2 that is remarkably constant with distance from the midplane and time; this can be used to improve estimates of ionized gas mass and mean electron density from observed H$\alpha$ surface brightness profiles in edge-on galaxies.

Spiral arms greatly affect gas flows and star formation in disk galaxies. We use local three-dimensional simulations of the vertically-stratified, self-gravitating, differentially-rotating, interstellar medium (ISM) subject to a stellar spiral potential to study the effects of spiral arms on star formation and formation of arm spurs/feathers. We adopt the TIGRESS framework of Kim & Ostriker (2017) to handle radiative heating and cooling, star formation, and ensuing supernova (SN) feedback. We find that more than 90% of star formation takes place in spiral arms, but the global star formation rate (SFR) in models with spiral arms is enhanced by less than a factor of 2 compared to the no-arm counterpart. This results from a quasi-linear relationship between the SFR surface density Sigma_SFR and the gas surface density Sigma, and supports the picture that spiral arms do not trigger star formation but rather concentrate star-forming regions. Correlated SN feedback produces gaseous spurs/feathers downstream from arms in both magnetized and unmagnetized models. These spurs/feathers are short-lived and have magnetic fields parallel to their length, in contrast to the longer-lived features with perpendicular magnetic fields induced by gravitational instability. SN feedback drives the turbulent component of magnetic fields, with the total magnetic field strength sublinearly proportional to Sigma. The total midplane pressure varies by a factor of ~10 between arm and interarm regions but agrees locally with the total vertical ISM weight, while Sigma_SFR is locally consistent with the prediction of pressure-regulated, feedback-modulated theory.

In typical environments of star-forming clouds, converging supersonic turbulence generates shock-compressed regions, and can create strongly-magnetized sheet-like layers. Numerical MHD simulations show that within these post-shock layers, dense filaments and embedded self-gravitating cores form via gathering material along the magnetic field lines. As a result of the preferred-direction mass collection, a velocity gradient perpendicular to the filament major axis is a common feature seen in simulations. We show that this prediction is in good agreement with recent observations from the CARMA Large Area Star Formation Survey (CLASSy), from which we identified several filaments with prominent velocity gradients perpendicular to their major axes. Highlighting a filament from the northwest part of Serpens South, we provide both qualitative and quantitative comparisons between simulation results and observational data. In particular, we show that the dimensionless ratio $C_v \equiv {\Delta v_h}^2/(GM/L)$, where $\Delta v_h$ is half of the observed perpendicular velocity difference across a filament, and $M/L$ is the filament's mass per unit length, can distinguish between filaments formed purely due to turbulent compression and those formed due to gravity-induced accretion. We conclude that the perpendicular velocity gradient observed in the Serpens South northwest filament can be caused by gravity-induced anisotropic accretion of material from a flattened layer. Using synthetic observations of our simulated filaments, we also propose that a density-selection effect may explain observed subfilaments (one filament breaking into two components in velocity space) as reported in Dhabal et al. (2018).

A common situation in galactic and intergalactic gas involves cold dense gas in motion relative to hot diffuse gas. Kelvin-Helmholtz instability creates a turbulent mixing layer and populates the intermediate-temperature phase, which often cools rapidly. The energy lost to cooling is balanced by the advection of hot high enthalpy gas into the mixing layer, resulting in growth and acceleration of the cold phase. This process may play a major role in determining the interstellar medium and circumgalactic medium phase structure, and accelerating cold gas in galactic winds and cosmic filaments. Cooling in these mixing layers occurs in a thin corrugated sheet, which we argue has an area with fractal dimension $D=5/2$ and a thickness that adjusts to match the hot phase mixing time to the cooling time. These cooling sheet properties form the basis of a new model for how the cooling rate and hot gas inflow velocity depend on the size $L$, cooling time $t_{\rm cool}$, relative velocity $v_{\rm rel}$, and density contrast $\rho_{\rm cold}/\rho_{\rm hot}$ of the system. Entrainment is expected to be enhanced in environments with short $t_{\rm cool}$, large $v_{\rm rel}$, and large $\rho_{\rm cold}/\rho_{\rm hot}$. Using a large suite of three dimensional hydrodynamic simulations, we demonstrate that this fractal cooling layer model accurately captures the energetics and evolution of turbulent interfaces and can therefore be used as a foundation for understanding multiphase mixing with strong radiative cooling.

We present the fourth of the Cholla Galactic OutfLow Simulations suite (CGOLS). Using a physically-motivated prescription for clustered supernova feedback, we successfully drive a multiphase outflow from a disk galaxy. The high resolution ($< 5\,\mathrm{pc}$) across a relatively large domain ($20\,\mathrm{kpc}$) allows us to capture the hydrodynamic mixing and dynamical interactions between the hot and cool ($T \sim 10^4\,\mathrm{K}$) phases in the outflow, which in turn leads to direct evidence of a qualitatively new mechanism for cool gas acceleration in galactic winds. We show that mixing of momentum from the hot phase to the cool phase accelerates the cool gas to $800\,\mathrm{km}\,\mathrm{s}^{-1}$ on kpc scales, with properties inconsistent with the physical models of ram pressure acceleration or with bulk cooling from the hot phase. The mixing process also affects the hot phase, modifying its radial profiles of temperature, density, and velocity from the expectations of radial supersonic flow. This mechanism provides a physical explanation for the high velocity, blue shifted, low ionization absorption lines often observed in the spectra of starburst and high redshift galaxies.

We compare the observed turbulent pressure in molecular gas, $P_\mathrm{turb}$, to the required pressure for the interstellar gas to stay in equilibrium in the gravitational potential of a galaxy, $P_\mathrm{DE}$. To do this, we combine arcsecond resolution CO data from PHANGS-ALMA with multi-wavelength data that traces the atomic gas, stellar structure, and star formation rate (SFR) for 28 nearby star-forming galaxies. We find that $P_\mathrm{turb}$ correlates with, but almost always exceeds the estimated $P_\mathrm{DE}$ on kiloparsec scales. This indicates that the molecular gas is over-pressurized relative to the large-scale environment. We show that this over-pressurization can be explained by the clumpy nature of molecular gas; a revised estimate of $P_\mathrm{DE}$ on cloud scales, which accounts for molecular gas self-gravity, external gravity, and ambient pressure, agrees well with the observed $P_\mathrm{turb}$ in galaxy disks. We also find that molecular gas with cloud-scale ${P_\mathrm{turb}}\approx{P_\mathrm{DE}}\gtrsim{10^5\,k_\mathrm{B}\,\mathrm{K\,cm^{-3}}}$ in our sample is more likely to be self-gravitating, whereas gas at lower pressure appears more influenced by ambient pressure and/or external gravity. Furthermore, we show that the ratio between $P_\mathrm{turb}$ and the observed SFR surface density, $\Sigma_\mathrm{SFR}$, is compatible with stellar feedback-driven momentum injection in most cases, while a subset of the regions may show evidence of turbulence driven by additional sources. The correlation between $\Sigma_\mathrm{SFR}$ and kpc-scale $P_\mathrm{DE}$ in galaxy disks is consistent with the expectation from self-regulated star formation models. Finally, we confirm the empirical correlation between molecular-to-atomic gas ratio and kpc-scale $P_\mathrm{DE}$ reported in previous works.

Galactic outflows produced by stellar feedback are known to be multiphase in nature. Both observations and simulations indicate that the material within several kpc of galactic disk mid-planes consists of warm clouds embedded within a hot wind. A theoretical understanding of the outflow phenomenon, including both winds and fountain flows, requires study of the interactions among thermal phases. We develop a method to quantify these interactions via measurements of mass, momentum, and energy flux exchanges using temporally and spatially averaged quantities and conservation laws. We apply this method to a star-forming ISM MHD simulation based on the TIGRESS framework, for Solar neighbourhood conditions. To evaluate the extent of interactions among the phases, we first examine the validity of the ``ballistic model,'' which predicts trajectories of the warm phase ($5050\,\rm{K}<T<2\times10^4\,\rm{K}$) treated as non-interacting clouds. This model is successful at intermediate vertical velocities ($ 50$ km s$^{-1}$ $\lesssim |v_z| \lesssim 100 $ km s$^{-1}$), but at higher velocities we observe an excess in simulated warm outflow compared to the ballistic model. This discrepancy cannot be fully accounted for by cooling of high-velocity intermediate-temperature ($2\times10^4\,\rm{K}<T<5\times10^5\,\rm{K}$) gas. By examining the fluxes of mass, momentum and energy, we conclude that warm phase gains mass via cooling of the intermediate phase, while momentum transfer occurs from the hot ($T>5\times10^5\,\rm{K}$) to the warm phase. The large energy flux from the hot outflow that is transferred to the warm and intermediate phases is quickly radiated away. A simple interaction model implies an effective warm cloud size in the fountain flow of a few 100~pc, showing that warm-hot flux exchange mainly involves a few large clouds rather than many small ones.

We apply gravity-based and density-based methods to identify clouds in numerical simulations of the star-forming, three-phase interstellar medium (ISM), and compare their properties and their global correlation with the star formation rate over time. The gravity-based method identifies bound objects, which have masses M ~ 10^3 - 10^4 M_solar at densities n_H ~ 100 cm^-3, and traditional virial parameters alpha_v ~ 0.5 - 5. For clouds defined by a density threshold n_H,min , the average virial parameter decreases, and the fraction of material that is genuinely bound increases, at higher n_H,min. Surprisingly, these clouds can be unbound even when alpha_v < 2, and high mass clouds (10^4 - 10^6 M_solar) are generally unbound. This suggests that the traditional alpha_v is at best an approximate measure of boundedness in the ISM. All clouds have internal turbulent motions increasing with size as sigma ~ 1 km/s(R/ pc)^1/2, similar to observed relations. Bound structures comprise a small fraction of the total simulation mass, with star formation efficiency per free-fall time epsilon_ff ~ 0.4. For n_H,min = 10 - 100 cm^-3, epsilon_ff ~ 0.03 - 0.3, increasing with density. Temporal correlation analysis between SFR(t) and aggregate mass M(n_H,min;t) at varying n_H,min shows that time delays to star formation are t_delay ~ t_ff(n_H,min). Correlation between SFR(t) and M(nH,min;t) systematically tightens at higher n_H,min. Considering moderate-density gas, selecting against high virial parameter clouds improves correlation with SFR, consistent with previous work. Even at high n_H,min, the temporal dispersion in (SFR-epsilon_ff M/t_ff )/<SFR> is ~ 50%, due to the large-amplitude variations and inherent stochasticity of the system.

We use spherically symmetric hydrodynamic simulations to study the dynamical evolution and internal structure of superbubbles (SBs) driven by clustered supernovae (SNe), focusing on the effects of thermal conduction and cooling in the interface between the hot bubble interior and cooled shell. Our simulations employ an effective diffusivity to account for turbulent mixing from nonlinear instabilities that are not captured in 1D. The conductive heat flux into the shell is balanced by a combination of cooling in the interface and evaporation of shell gas into the bubble interior. This evaporation increases the density, and decreases the temperature, of the SB interior by more than an order of magnitude relative to simulations without conduction. However, most of the energy conducted into the interface is immediately lost to cooling, reducing the evaporative mass flux required to balance conduction. As a result, the evaporation rate is typically a factor of $\sim$3-30 lower than predicted by the classical similarity solution of Weaver et al. (1977), which neglects cooling. Blast waves from the first $\sim$30 SNe remain supersonic in the SB interior because reduced evaporation from the interface lowers the mass they sweep up in the hot interior. Updating the Weaver solution to include cooling, we construct a new analytic model to predict the cooling rate, evaporation rate, and temporal evolution of SBs. The cooling rate, and hence the hot gas mass, momentum, and energy delivered by SBs, is set by the ambient ISM density and the efficiency of nonlinear mixing at the bubble/shell interface.

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

Using a suite of radiation hydrodynamic simulations of star cluster formation in turbulent clouds, we study the escape fraction of ionizing (Lyman continuum) and non-ionizing (FUV) radiation for a wide range of cloud masses and sizes. The escape fraction increases as H II regions evolve and reaches unity within a few dynamical times. The cumulative escape fraction before the onset of the first supernova explosion is in the range 0.05-0.58; this is lower for higher initial cloud surface density, and higher for less massive and more compact clouds due to rapid destruction. Once H II regions break out of their local environment, both ionizing and non-ionizing photons escape from clouds through fully ionized, low-density sightlines. Consequently, dust becomes the dominant absorber of ionizing radiation at late times and the escape fraction of non-ionizing radiation is only slightly larger than that of ionizing radiation. The escape fraction is determined primarily by the mean $\langle \tau\rangle$ and width $\sigma$ of the optical-depth distribution in the large-scale cloud, increasing for smaller $\langle \tau\rangle$ and/or larger $\sigma$. The escape fraction exceeds (sometimes by three orders of magnitude) the naive estimate $e^{-\langle \tau\rangle}$ due to non-zero $\sigma$ induced by turbulence. We present two simple methods to estimate, within $\sim20\%$, the escape fraction of non-ionizing radiation using the observed dust optical depth in clouds projected on the plane of sky. We discuss implications of our results for observations, including inference of star formation rates in individual molecular clouds, and accounting for diffuse ionized gas on galactic scales.

We investigate the prevalence, properties, and kinematics of extraplanar diffuse ionized gas (eDIG) in a sample of 25 edge-on galaxies selected from the CALIFA survey. We measure ionized gas scale heights from ${\rm H\alpha}$ and find that 90% have measurable scale heights with a median of $0.8^{+0.7}_{-0.4}$ kpc. From the ${\rm H\alpha}$ kinematics, we find that 60% of galaxies show a decrease in the rotation velocity as a function of height above the midplane. This lag is characteristic of eDIG, and we measure a median lag of 21 km s$^{-1}$ kpc$^{-1}$ which is comparable to lags measured in the literature. We also investigate variations in the lag with radius. $\rm H{\small I}$ lags have been reported to systematically decrease with galactocentric radius. We find both increasing and decreasing ionized gas lags with radius, as well as a large number of galaxies consistent with no radial lag variation, and investigate these results in the context of internal and external origins for the lagging ionized gas. We confirm that the ${\rm [S{\small II}]}$/${\rm H\alpha}$ and ${\rm [N{\small II}]}$/${\rm H\alpha}$ line ratios increase with height above the midplane as is characteristic of eDIG. The ionization of the eDIG is dominated by star-forming complexes (leaky ${\rm H{\small II}}$ regions). We conclude that the lagging ionized gas is turbulent ejected gas likely resulting from star formation activity in the disk as opposed to gas in the stellar thick disk or bulge. This is further evidence for the eDIG being a product of stellar feedback and for the pervasiveness of this WIM-like phase in many local star-forming galaxies.

Chemistry plays an important role in the interstellar medium (ISM), regulating heating and cooling of the gas, and determining abundances of molecular species that trace gas properties in observations. Although solving the time-dependent equations is necessary for accurate abundances and temperature in the dynamic ISM, a full chemical network is too computationally expensive to incorporate in numerical simulations. In this paper, we propose a new simplified chemical network for hydrogen and carbon chemistry in the atomic and molecular ISM. We compare results from our chemical network in detail with results from a full photo-dissociation region (PDR) code, and also with the Nelson & Langer (1999) (NL99) network previously adopted in the simulation literature. We show that our chemical network gives similar results to the PDR code in the equilibrium abundances of all species over a wide range of densities, temperature, and metallicities, whereas the NL99 network shows significant disagreement. Applying our network in 1D models, we find that the $\mathrm{CO}$-dominated regime delimits the coldest gas and that the corresponding temperature tracks the cosmic ray ionization rate in molecular clouds. We provide a simple fit for the locus of $\mathrm{CO}$ dominated regions as a function of gas density and column. We also compare with observations of diffuse and translucent clouds. We find that the $\mathrm{CO}$, $\mathrm{CHx}$ and $\mathrm{OHx}$ abundances are consistent with equilibrium predictions for densities $n=100-1000~\mathrm{cm^{-3}}$, but the predicted equilibrium $\mathrm{C}$ abundance is higher than observations, signaling the potential importance of non-equilibrium/dynamical effects.

The gyro-resonant cosmic-ray (CR) streaming instability is believed to play a crucial role in CR transport, leading to growth of Alfv\'en waves at small scales that scatter CRs, and impacts the interaction of CRs with the ISM on large scales. However, extreme scale separation ($\lambda \ll \rm pc$), low cosmic ray number density ($n_{\rm CR}/n_{\rm ISM} \sim 10^{-9}$), and weak CR anisotropy ($\sim v_A/c$) pose strong challenges for proper numerical studies of this instability on a microphysical level. Employing the recently developed magnetohydrodynamic-particle-in-cell (MHD-PIC) method, which has unique advantages to alleviate these issues, we conduct one-dimensional simulations that quantitatively demonstrate the growth and saturation of the instability in the parameter regime consistent with realistic CR streaming in the large-scale ISM. Our implementation of the $\delta f$ method dramatically reduces Poisson noise and enables us to accurately capture wave growth over a broad spectrum, equally shared between left and right handed Alfv\'en modes. We are also able to accurately follow the quasi-linear diffusion of CRs subsequent to wave growth, which is achieved by employing phase randomization across periodic boundaries. Full isotropization of the CRs in the wave frame requires pitch angles of most CRs to efficiently cross $90^\circ$, and can be captured in simulations with relatively high wave amplitude and/or high spatial resolution. We attribute this crossing to non-linear wave-particle interaction (rather than mirror reflection) by investigating individual CR trajectories. We anticipate our methodology will open up opportunities for future investigations that incorporate additional physics.

We present an accurate and efficient method to calculate the gravitational potential of an isolated system in three-dimensional Cartesian and cylindrical coordinates subject to vacuum (open) boundary conditions. Our method consists of two parts: an interior solver and a boundary solver. The interior solver adopts an eigenfunction expansion method together with a tridiagonal matrix solver to solve the Poisson equation subject to the zero boundary condition. The boundary solver employs James's method to calculate the boundary potential due to the screening charges required to keep the zero boundary condition for the interior solver. A full computation of gravitational potential requires running the interior solver twice and the boundary solver once. We develop a method to compute the discrete Green's function in cylindrical coordinates, which is an integral part of the James algorithm to maintain second-order accuracy. We implement our method in the {\tt Athena++} magnetohydrodynamics code, and perform various tests to check that our solver is second-order accurate and exhibits good parallel performance.

Recent observations have revealed that the remnants of stellar-coalescence transients are bipolar. This raises the questions of how these bipolar morphologies arise and what they teach us about the mechanisms of mass ejection during stellar mergers and common-envelope phases. In this paper, we analyze hydrodynamic simulations of the lead-in to binary coalescence, a phase of unstable Roche lobe overflow that takes the binary from the Roche limit separation to the engulfment of the more compact accretor within the envelope of the extended donor. As mass transfer runs away at increasing rates, gas trails away from the binary. Contrary to previous expectations, early mass loss from the system remains bound to the binary and forms a circumbinary torus. Later ejecta, generated as the accretor grazes the surface of the donor, have very different morphologies and are unbound. These two components of mass loss from the binary interact as later, higher-velocity ejecta collide with the circumbinary torus formed by earlier mass loss. Unbound ejecta are redirected toward the poles, and escaping material creates a bipolar outflow. Our findings show that the transition from bound to unbound ejecta from coalescing binaries can explain the bipolar nature of their remnants, with implications for our understanding of the origin of bipolar remnants of stellar-coalescence transients and, perhaps, some preplanetary nebulae.

NGC 253 hosts the nearest nuclear starburst. Previous observations show a region rich in molecular gas, with dense clouds associated with recent star formation. We used ALMA to image the 350 GHz dust continuum and molecular line emission from this region at 2 pc resolution. Our observations reveal ~14 bright, compact (~2-3 pc FWHM) knots of dust emission. Most of these sources are likely to be forming super star clusters (SSCs) based on their inferred dynamical and gas masses, association with 36 GHz radio continuum emission, and coincidence with line emission tracing dense, excited gas. One source coincides with a known SSC, but the rest remain invisible in Hubble near-infrared (IR) imaging. Our observations imply that gas still constitutes a large fraction of the overall mass in these sources. Their high brightness temperature at 350 GHz also implies a large optical depth near the peak of the IR spectral energy distribution. As a result, these sources may have large IR photospheres and the IR radiation force likely exceeds L/c. Still, their moderate observed velocity dispersions suggest that feedback from radiation, winds, and supernovae are not yet disrupting most sources. This mode of star formation appears to produce a large fraction of stars in the burst. We argue for a scenario in which this phase lasts ~1 Myr, after which the clusters shed their natal cocoons but continue to produce ionizing photons. The strong feedback that drives the observed cold gas and X-ray outflows likely occurs after the clusters emerge from this early phase.

We utilize the more than 100 gravitationally-bound dense cores formed in our three-dimensional, turbulent MHD simulations reported in Chen & Ostriker (2015) to analyze structural, kinematic, and magnetic properties of prestellar cores. Our statistical results disagree with the classical theory of star formation, in which cores evolve to be oblate with magnetic field parallel to the minor axes. Instead, we find that cores are generally triaxial, although the core-scale magnetic field is still preferentially most parallel to the core's minor axis and most perpendicular to the major axis. The internal and external magnetic field directions are correlated, but the direction of integrated core angular momentum is misaligned with the core's magnetic field, consistent with recent observations. The ratio of rotational/total kinetic and rotational/gravitational energies are independent of core size and consistent in magnitude with observations. The specific angular momentum also follows the observed relationship $L/M \propto R^{3/2}$, indicating rotation is acquired from ambient turbulence. With typical $E_\mathrm{rot}/E_K \sim 0.1$, rotation is not the dominant motion when cores collapse.

Luminous red nova transients, presumably from stellar coalescence, exhibit long-term precursor emission over hundreds of binary orbits, leading to impulsive outbursts with durations similar to a single orbital period. In an effort to understand these signatures, we present and analyze a hydrodynamic model of unstable mass transfer from a giant-star donor onto a more compact accretor in a binary system. Our simulation begins with mass transfer at the Roche limit separation and traces a phase of runaway decay leading up to the plunge of the accretor within the envelope of the donor. We characterize the fluxes of mass and angular momentum through the system and show that the orbital evolution can be reconstructed from measurements of these quantities. The morphology of outflow from the binary changes significantly as the binary orbit tightens. At wide separations, a thin stream of relatively high-entropy gas trails from the outer Lagrange points. As the orbit tightens, the orbital motion desynchronizes from the donor's rotation, and low-entropy ejecta trace a broad fan of largely ballistic trajectories. An order-of-magnitude increase in mass ejection rate accompanies the plunge of the accretor with the envelope of the donor. We argue that this transition marks the precursor-to-outburst transition observed in stellar coalescence transients.

We estimate the star formation efficiency per gravitational free fall time, $\epsilon_{\rm ff}$, from observations of nearby galaxies with resolution matched to the typical size of a Giant Molecular Cloud. This quantity, $\epsilon_{\rm ff}$, is theoretically important but so far has only been measured for Milky Way clouds or inferred indirectly in a few other galaxies. Using new, high resolution CO imaging from the PHANGS-ALMA survey, we estimate the gravitational free-fall time at 60 to 120 pc resolution, and contrast this with the local molecular gas depletion time to estimate $\epsilon_{\rm ff}$. Assuming a constant thickness of the molecular gas layer ($H = 100$ pc) across the whole sample, the median value of $\epsilon_{\rm ff}$ in our sample is $0.7\%$. We find a mild scale-dependence, with higher $\epsilon_{\rm ff}$ measured at coarser resolution. Individual galaxies show different values of $\epsilon_{\rm ff}$, with the median $\epsilon_{\rm ff}$ ranging from $0.3\%$ to $2.6\%$. We find the highest $\epsilon_{\rm ff}$ in our lowest mass targets, reflecting both long free-fall times and short depletion times, though we caution that both measurements are subject to biases in low mass galaxies. We estimate the key systematic uncertainties, and show the dominant uncertainty to be the estimated line-of-sight depth through the molecular gas layer and the choice of star formation tracers.

UV radiation feedback from young massive stars plays a key role in the evolution of giant molecular clouds (GMCs) by photoevaporating and ejecting the surrounding gas. We conduct a suite of radiation hydrodynamic simulations of star cluster formation in marginally-bound, turbulent GMCs, focusing on the effects of photoionization and radiation pressure on regulating the net star formation efficiency (SFE) and cloud lifetime. We find that the net SFE depends primarily on the initial gas surface density, $\Sigma_0$, such that the SFE increases from 4% to 51% as $\Sigma_0$ increases from $13\,M_{\odot}\,{\rm pc}^{-2}$ to $1300\,M_{\odot}\,{\rm pc}^{-2}$. Cloud destruction occurs within $2$-$10\,{\rm Myr}$ after the onset of radiation feedback, or within $0.6$-$4.1$ freefall times (increasing with $\Sigma_0$). Photoevaporation dominates the mass loss in massive, low surface-density clouds, but because most photons are absorbed in an ionization-bounded Str\"{o}mgren volume the photoevaporated gas fraction is proportional to the square root of the SFE. The measured momentum injection due to thermal and radiation pressure forces is proportional to $\Sigma_0^{-0.74}$, and the ejection of neutrals substantially contributes to the disruption of low-mass and/or high-surface density clouds. We present semi-analytic models for cloud dispersal mediated by photoevaporation and by dynamical mass ejection, and show that the predicted net SFE and mass loss efficiencies are consistent with the results of our numerical simulations.

We use new ALMA observations to investigate the connection between dense gas fraction, star formation rate, and local environment across the inner region of four local galaxies showing a wide range of molecular gas depletion times. We map HCN (1-0), HCO$^+$ (1-0), CS (2-1), $^{13}$CO (1-0), and C$^{18}$O (1-0) across the inner few kpc of each target. We combine these data with short spacing information from the IRAM large program EMPIRE, archival CO maps, tracers of stellar structure and recent star formation, and recent HCN surveys by Bigiel et al. and Usero et al. We test the degree to which changes in the dense gas fraction drive changes in the SFR. $I_{HCN}/I_{CO}$ (tracing the dense gas fraction) correlates strongly with $I_{CO}$ (tracing molecular gas surface density), stellar surface density, and dynamical equilibrium pressure, $P_{DE}$. Therefore, $I_{HCN}/I_{CO}$ becomes very low and HCN becomes very faint at large galactocentric radii, where ratios as low as $I_{HCN}/I_{CO} \sim 0.01$ become common. The apparent ability of dense gas to form stars, $\Sigma_{SFR}/\Sigma_{dense}$ (where $\Sigma_{dense}$ is traced by the HCN intensity and the star formation rate is traced by a combination of H$\alpha$ and 24$\mu$m emission), also depends on environment. $\Sigma_{SFR}/\Sigma_{dense}$ decreases in regions of high gas surface density, high stellar surface density, and high $P_{DE}$. Statistically, these correlations between environment and both $\Sigma_{SFR}/\Sigma_{dense}$ and $I_{HCN}/I_{CO}$ are stronger than that between apparent dense gas fraction ($I_{HCN}/I_{CO}$) and the apparent molecular gas star formation efficiency $\Sigma_{SFR}/\Sigma_{mol}$. We show that these results are not specific to HCN.

CO(J=1-0) line emission is a widely used observational tracer of molecular gas, rendering essential the X_CO factor, which is applied to convert CO luminosity to H_2 mass. We use numerical simulations to study how X_CO depends on numerical resolution, non-steady-state chemistry, physical environment, and observational beam size. Our study employs 3D magnetohydrodynamics (MHD) simulations of galactic disks with solar neighborhood conditions, where star formation and the three-phase interstellar medium (ISM) are self-consistently regulated by gravity and stellar feedback. Synthetic CO maps are obtained by post-processing the MHD simulations with chemistry and radiation transfer. We find that CO is only an approximate tracer of H_2. On parsec scales, W_CO is more fundamentally a measure of mass-weighted volume density, rather than H_2 column density. Nevertheless, $\langle X_\mathrm{CO} \rangle=0.7-1.0\times10^{20}~\mathrm{cm^{-2}K^{-1}km^{-1}s}$ consistent with observations, insensitive to the evolutionary ISM state or radiation field strength if steady-state chemistry is assumed. Due to non-steady-state chemistry, younger molecular clouds have slightly lower X_CO and flatter profiles of X_CO versus extinction than older ones. The CO-dark H_2 fraction is 26-79 %, anti-correlated with the average extinction. As the observational beam size increases from 1 pc to 100 pc, X_CO increases by a factor of ~ 2. Under solar neighborhood conditions, X_CO in molecular clouds is converged at a numerical resolution of 2 pc. However, the total CO abundance and luminosity are not converged even at the numerical resolution of 1 pc. Our simulations successfully reproduce the observed variations of X_CO on parsec scales, as well as the dependence of X_CO on extinction and the CO excitation temperature.

Cosmic ray pressure gradients transfer energy and momentum to extraplanar gas in disk galaxies, potentially driving significant mass loss as galactic winds. This may be particularly important for launching high-velocity outflows of "cool" (T < 10^4 K) gas. We study cosmic-ray driven disk winds using a simplified semi-analytic model assuming streamlines follow the large-scale gravitational potential gradient. We consider scaled Milky Way-like potentials including a disk, bulge, and halo with a range of halo velocities V_H = 50-300 km/s, and streamline footpoints with radii in the disk R_0=1-16 kpc at height 1 kpc. Our solutions cover a wide range of footpoint gas velocity u_0, magnetic-to-cosmic-ray pressure ratio, gas-to-cosmic-ray pressure ratio, and angular momentum. Cosmic ray streaming at the Alfv\'en speed enables the effective sound speed C_eff to increase from the footpoint to a critical point where C_eff,c = u_c ~ V_H; this differs from thermal winds in which C_eff decreases outward. The critical point is typically at a height of 1-6 kpc from the disk, increasing with V_H, and the asymptotic wind velocity exceeds the escape speed of the halo. Mass loss rates are insensitive to the footpoint values of the magnetic field and angular momentum. In addition to numerical parameter space exploration, we develop and compare to analytic scaling relations. We show that winds have mass loss rates per unit area up to ~ Pi_0 V_H^-5/3 u_0^2/3 where Pi_0 is the footpoint cosmic ray pressure and u_0 is set by the upwelling of galactic fountains. The predicted wind mass-loss rate exceeds the star formation rate for V_H < 200 km/s and u_0 = 50 km/s, a typical fountain velocity.

Gas blown away from galactic disks by supernova (SN) feedback plays a key role in galaxy evolution. We investigate outflows utilizing the solar neighborhood model of our high-resolution, local galactic disk simulation suite, TIGRESS. In our numerical implementation, star formation and SN feedback are self-consistently treated and well resolved in the multiphase, turbulent, magnetized interstellar medium. Bursts of star formation produce spatially and temporally correlated SNe that drive strong outflows, consisting of hot (T>5x10^5K) winds and warm (5050K < T < 2x10^4K) fountains. The hot gas at distance d>1kpc from the midplane has mass and energy fluxes nearly constant with d. The hot flow escapes our local Cartesian box barely affected by gravity and is expected to accelerate up to the terminal velocity of v_wind~350-500km/s. The mean mass and energy loading factors of the hot wind are 0.1 and 0.02, respectively. For warm gas, the mean outward mass flux through d=1kpc is comparable to the mean star formation rate, but only a small fraction of this gas is at velocity >50km/s. Thus, the warm outflows eventually fall back as inflows. The warm fountain flows are created by expanding hot superbubbles at d< 1kpc; at larger d neither ram pressure acceleration nor cooling transfers significant momentum or energy flux from the hot wind to the warm outflow. The velocity distribution at launching near d~1kpc better represents warm outflows than a single mass loading factor, potentially enabling development of subgrid models for warm galactic winds in arbitrary large-scale galactic potentials.

We present a kpc-scale analysis of the relationship between the molecular depletion time ($\tau_\mathrm{dep}^\mathrm{mol}$) and the orbital time ($\tau_\mathrm{orb}$) across the field of 39 face-on local galaxies, selected from the EDGE-CALIFA sample. We find that, on average, 5% of the available molecular gas is converted into stars per orbital time, or $\tau_\mathrm{dep}^\mathrm{mol}\sim20\tau_\mathrm{orb}$. The resolved relation shows a scatter of $\sim0.5$ dex. The scatter is ascribable to galaxies of different morphologies that follow different $\tau_\mathrm{dep}^\mathrm{mol}-\tau_\mathrm{orb}$ relations which decrease in steepness from early- to late-types. The morphologies appear to be linked with the star formation rate surface density, the molecular depletion time, and the orbital time, but they do not correlate with the molecular gas content of the galaxies in our sample. We speculate that in our molecular gas rich, early-type galaxies, the morphological quenching (in particular the disc stabilization via shear), rather than the absence of molecular gas, is the main factor responsible for their current inefficient star formation.

Momentum deposition by radiation pressure from young, massive stars may help to destroy molecular clouds and unbind stellar clusters by driving large-scale outflows. We extend our previous numerical radiation hydrodynamic study of turbulent, star-forming clouds to analyze the detailed interaction between non-ionizing UV radiation and the cloud material. Our simulations trace the evolution of gas and star particles through self-gravitating collapse, star formation, and cloud destruction via radiation-driven outflows. These models are idealized in that we include only radiation feedback and adopt an isothermal equation of state. Turbulence creates a structure of dense filaments and large holes through which radiation escapes, such that only ~50% of the radiation is (cumulatively) absorbed by the end of star formation. The surface density distribution of gas by mass as seen by the central cluster is roughly lognormal with sigma_ln(Sigma) = 1.3-1.7, similar to the externally-projected surface density distribution. This allows low surface density regions to be driven outwards to nearly 10 times their initial escape speed v_esc. Although the velocity distribution of outflows is broadened by the lognormal surface density distribution, the overall efficiency of momentum injection to the gas cloud is reduced because much of the radiation escapes. The mean outflow velocity is approximately twice the escape speed from the initial cloud radius. Our results are also informative for understanding galactic-scale wind driving by radiation, in particular the relationship between velocity and surface density for individual outflow structures, and the resulting velocity and mass distributions arising from turbulent sources.

We present an implementation of an adaptive ray tracing (ART) module in the Athena hydrodynamics code that accurately and efficiently handles the radiative transfer involving multiple point sources on a three-dimensional Cartesian grid. We adopt a recently proposed parallel algorithm that uses non-blocking, asynchronous MPI communications to accelerate transport of rays across the computational domain. We validate our implementation through several standard test problems including the propagation of radiation in vacuum and the expansions of various types of HII regions. Additionally, scaling tests show that the cost of a full ray trace per source remains comparable to that of the hydrodynamics update on up to $\sim 10^3$ processors. To demonstrate application of our ART implementation, we perform a simulation of star cluster formation in a marginally bound, turbulent cloud, finding that its star formation efficiency is $12\%$ when both radiation pressure forces and photoionization by UV radiation are treated. We directly compare the radiation forces computed from the ART scheme with that from the M1 closure relation. Although the ART and M1 schemes yield similar results on large scales, the latter is unable to resolve the radiation field accurately near individual point sources.

We present results from the EDGE survey, a spatially resolved CO(1-0) follow-up to CALIFA, an optical Integral Field Unit (IFU) survey of local galaxies. By combining the data products of EDGE and CALIFA, we study the variation in molecular gas depletion time ($\tau_{\rm dep}$) on kiloparsec scales in 52 galaxies. We divide each galaxy into two parts: the center, defined as the region within $0.1 \ R_{25}$, and the disk, defined as the region between $0.1$ and $0.7 \ R_{25}$. We find that 14 galaxies show a shorter $\tau_{\rm dep}$ ($\sim 1$ Gyr) in the center relative to that in the disk ($\tau_{\rm dep} \sim 2.4$ Gyrs), which means the central region in those galaxies is more efficient at forming stars per unit molecular gas mass. This finding implies that the centers with shorter $\tau_{\rm dep}$ resemble the intermediate regime between galactic disks and starburst galaxies. Furthermore, the central drop in $\tau_{\rm dep}$ is correlated with a central increase in the stellar surface density, suggesting that a shorter $\tau_{\rm dep}$ is associated with molecular gas compression by the stellar gravitational potential. We argue that varying the CO-to-H$_2$ conversion factor only exaggerates the central drop of $\tau_{\rm dep}$.

Winds arising from galaxies, star clusters, and active galactic nuclei are crucial players in star and galaxy formation, but it has proven remarkably difficult to use observations of them to determine physical properties of interest, particularly mass fluxes. Much of the difficulty stems from a lack of a theory that links a physically-realistic model for winds' density, velocity, and covering factors to calculations of light emission and absorption. In this paper we provide such a model. We consider a wind launched from a turbulent region with a range of column densities, derive the differential acceleration of gas as a function of column density, and use this result to compute winds' absorption profiles, emission profiles, and emission intensity maps in both optically thin and optically thick species. The model is sufficiently simple that all required computations can be done analytically up to straightforward numerical integrals, rendering it suitable for the problem of deriving physical parameters by fitting models to observed data. We show that our model produces realistic absorption and emission profiles for some example cases, and argue that the most promising methods of deducing mass fluxes are based on combinations of absorption lines of different optical depths, or on combining absorption with measurements of molecular line emission. In the second paper in this series, we expand on these ideas by introducing a set of observational diagnostics that are significantly more robust that those commonly in use, and that can be used to obtain improved estimates of wind properties.

We present interferometric CO observations made with the Combined Array for Millimeter-wave Astronomy (CARMA) of galaxies from the Extragalactic Database for Galaxy Evolution survey (EDGE). These galaxies are selected from the Calar Alto Legacy Integral Field Area (CALIFA) sample, mapped with optical integral field spectroscopy. EDGE provides good quality CO data (3$\sigma$ sensitivity $\Sigma_{\rm mol}$ $\sim$ 11 M$_\odot$ pc$^{-2}$ before inclination correction, resolution $\sim1.4$ kpc) for 126 galaxies, constituting the largest interferometric CO survey of galaxies in the nearby universe. We describe the survey, the data characteristics, the data products, and present initial science results. We find that the exponential scale-lengths of the molecular, stellar, and star-forming disks are approximately equal, and galaxies that are more compact in molecular gas than in stars tend to show signs of interaction. We characterize the molecular to stellar ratio as a function of Hubble type and stellar mass, present preliminary results on the resolved relations between the molecular gas, stars, and star formation rate, and discuss the dependence of the resolved molecular depletion time on stellar surface density, nebular extinction, and gas metallicity. EDGE provides a key dataset to address outstanding topics regarding gas and its role in star formation and galaxy evolution, which will be publicly available on completion of the quality assessment.

We introduce TIGRESS, a novel framework for multi-physics numerical simulations of the star-forming interstellar medium (ISM) implemented in the Athena MHD code. The algorithms of TIGRESS are designed to spatially and temporally resolve key physical features, including: (1) the gravitational collapse and ongoing accretion of gas that leads to star formation in clusters, (2) the explosions of supernovae (SNe) both near their progenitor birth sites and from runaway OB stars, with time delays relative to star formation determined by population synthesis, (3) explicit evolution of SN remnants prior to the onset of cooling, which leads to the creation of the hot ISM, (4) photoelectric heating of the warm and cold phases of the ISM that tracks the time-dependent ambient FUV field from the young cluster population, (5) large-scale galactic differential rotation, which leads to epicyclic motion and shears out overdense structures, limiting large-scale gravitational collapse, (6) accurate evolution of magnetic fields, which can be important for vertical support of the ISM disk as well as angular momentum transport. We present tests of the newly-implemented physics modules, and demonstrate application of TIGRESS in a fiducial model representing the Solar neighborhood environment. We use a resolution study to demonstrate convergence and evaluate the minimum resolution dx required to correctly recover several ISM properties, including the star formation rate, wind mass-loss rate, disk scale height, turbulent and Alfvenic velocity dispersions, and volume fractions of warm and hot phases. For the Solar neighborhood model, all these ISM properties are converged at dx <= 8pc.

We analyze synthetic neutral hydrogen (HI) absorption and emission spectral lines from a high- resolution, three-dimensional hydrodynamical simulation to quantify how well observational methods recover the physical properties of interstellar gas. We present a new method for uniformly decomposing HI spectral lines and estimating the properties of associated gas using the Autonomous Gaussian Decomposition (AGD) algorithm. We find that Hi spectral lines recover physical structures in the simulation with excellent completeness at high Galactic latitude, and this completeness declines with decreasing latitude due to strong velocity-blending of spectral lines. The temperature and column density inferred from our decomposition and radiative transfer method agree with the simulated values within a factor of < 2 for the majority of gas structures. We next compare synthetic spectra with observations from the 21-SPONGE survey at the Karl G. Jansky Very Large Array using AGD. We find more components per line of sight in 21-SPONGE than in synthetic spectra, which reflects insufficient simulated gas scale heights and the limitations of local box simulations. In addition, we find a significant population of low-optical depth, broad absorption components in the synthetic data which are not seen in 21-SPONGE. This population is not obvious in integrated or per-channel diagnostics, and reflects the benefit of studying velocity-resolved components. The discrepant components correspond to the highest spin temperatures (1000 < Ts < 4000 K), which are not seen in 21-SPONGE despite sufficient observational sensitivity. We demonstrate that our analysis method is a powerful tool for diagnosing neutral ISM conditions, and future work is needed to improve observational statistics and implementation of simulated physics.

We present a detailed study of a molecular outflow feature in the nearby starburst galaxy NGC 253 using ALMA. We find that this feature is clearly associated with the edge of NGC 253's prominent ionized outflow, has a projected length of ~300 pc, with a width of ~50 pc and a velocity dispersion of ~40 km s^-1, consistent with an ejection from the disk about 1 Myr ago. The kinematics of the molecular gas in this feature can be interpreted (albeit not uniquely) as accelerating at a rate of 1 km s^-1 pc^-1. In this scenario, the gas is approaching escape velocity at the last measured point. Strikingly, bright tracers of dense molecular gas (HCN, CN, HCO+, CS) are also detected in the molecular outflow: We measure an HCN(1-0)/CO(1-0) line ratio of ~1/10 in the outflow, similar to that in the central starburst region of NGC 253 and other starburst galaxies. By contrast, the HCN/CO line ratio in the NGC 253 disk is significantly lower (~1/30), similar to other nearby galaxy disks. This strongly suggests that the streamer gas originates from the starburst, and that its physical state does not change significantly over timescales of ~1 Myr during its entrainment in the outflow. Simple calculations indicate that radiation pressure is not the main mechanism for driving the outflow. The presence of such dense material in molecular outflows needs to be accounted for in simulations of galactic outflows.

Dynamic and thermal processes regulate the structure of the multi-phase interstellar medium (ISM), and ultimately establish how galaxies evolve through star formation. Thus, to constrain ISM models and better understand the interplay of these processes, it is of great interest to measure the thermal pressure ($P_{\rm th}$) of the diffuse, neutral gas. By combining [C II] 158 $\mu$m, HI, and CO data from 31 galaxies selected from the Herschel KINGFISH sample, we have measured thermal pressures in 534 predominantly atomic regions with typical sizes of $\sim$1 kiloparsec. We find a distribution of thermal pressures in the $P_{\rm th}/k\sim10^3-10^5$ K cm$^{-3}$ range. For a sub-sample of regions with conditions similar to those of the diffuse, neutral gas in the Galactic plane, we find thermal pressures that follow a log-normal distribution with a median value of $P_{\rm th}/k\approx3600$ K cm$^{-3}$. These results are consistent with thermal pressure measurements using other observational methods. We find that $P_{\rm th}$ increases with radiation field strength and star formation activity, as expected from the close link between the heating of the gas and the star formation rate. Our thermal pressure measurements fall in the regime where a two-phase ISM with cold and warm neutral medium could exist in pressure equilibrium. Finally, we find that the midplane thermal pressure of the diffuse gas is about $\sim30$% of the vertical weight of the overlying ISM, consistent with results from hydrodynamical simulations of self-regulated star formation in galactic disks.

We use numerical simulations to analyze the evolution and properties of superbubbles (SBs), driven by multiple supernovae (SNe), that propagate into the two-phase (warm/cold), cloudy interstellar medium (ISM). We consider a range of mean background densities n_avg=0.1-10 cm^{-3} and intervals between SNe dt_sn=0.01-1 Myr, and follow each SB until the radius reaches (1-2)H, where H is the characteristic ISM disk thickness. Except for embedded dense clouds, each SB is hot until a time t_sf,m when the shocked warm gas at the outer front cools and forms an overdense shell. Subsequently, diffuse gas in the SB interior remains at T_h 10^6-10^7K with expansion velocity v_h~10^2-10^3km/s (both highest for low dt_sn). At late times, the warm shell gas velocities are several 10's to ~100km/s. While shell velocities are too low to escape from a massive galaxy, they are high enough to remove substantial mass from dwarfs. Dense clouds are also accelerated, reaching a few to 10's of km/s. We measure the mass in hot gas per SN, M_h/N_SN, and the total radial momentum of the bubble per SN, p_b/N_SN. After t_sf,m, M_h/N_SN 10-100M_sun (highest for low n_avg), while p_b/N_SN 0.7-3x10^5M_sun km/s (highest for high dt_sn). If galactic winds in massive galaxies are loaded by the hot gas in SBs, we conclude that the mass-loss rates would generally be lower than star formation rates. Only if the SN cadence is much higher than typical in galactic disks, as may occur for nuclear starbursts, SBs can break out while hot and expel up to 10 times the mass locked up in stars. The momentum injection values, p_b/N_SN, are consistent with requirements to control star formation rates in galaxies at observed levels.

Radiation feedback from stellar clusters is expected to play a key role in setting the rate and efficiency of star formation in giant molecular clouds (GMCs). To investigate how radiation forces influence realistic turbulent systems, we have conducted a series of numerical simulations employing the {\it Hyperion} radiation hydrodynamics solver, considering the regime that is optically thick to ultraviolet (UV) and optically thin to infrared (IR) radiation. Our model clouds cover initial surface densities between $\Sigma_{\rm cl,0} \sim 10-300~M_{\odot}~{\rm pc^{-2}}$, with varying initial turbulence. We follow them through turbulent, self-gravitating collapse, formation of star clusters, and cloud dispersal by stellar radiation. All our models display a lognormal distribution of gas surface density $\Sigma$; for an initial virial parameter $\alpha_{\rm vir,0} = 2$, the lognormal standard deviation is $\sigma_{\rm ln \Sigma} = 1-1.5$ and the star formation rate coefficient $\varepsilon_{\rm ff,\bar\rho} = 0.3-0.5$, both of which are sensitive to turbulence but not radiation feedback. The net star formation efficiency $\varepsilon_\mathrm{final}$ increases with $\Sigma_{\rm cl,0}$ and decreases with $\alpha_{\rm vir,0}$. We interpret these results via a simple conceptual framework, whereby steady star formation increases the radiation force, such that local gas patches at successively higher $\Sigma$ become unbound. Based on this formalism (with fixed $\sigma_{\rm ln \Sigma}$), we provide an analytic upper bound on $\varepsilon_\mathrm{final}$, which is in good agreement with our numerical results. The final star formation efficiency depends on the distribution of Eddington ratios in the cloud and is strongly increased by turbulent compression of gas.

To understand the conditions under which dense, molecular gas is able to form within a galaxy, we post-process a series of three-dimensional galactic-disk-scale simulations with ray-tracing based radiative transfer and chemical network integration to compute the equilibrium chemical and thermal state of the gas. In performing these simulations we vary a number of parameters, such as the ISRF strength, vertical scale height of stellar sources, cosmic ray flux, to gauge the sensitivity of our results to these variations. Self-shielding permits significant molecular hydrogen (H2) abundances in dense filaments around the disk midplane, accounting for approximately ~10-15% of the total gas mass. Significant CO fractions only form in the densest, n>~10^3 cm^-3, gas where a combination of dust, H2, and self-shielding attenuate the FUV background. We additionally compare these ray-tracing based solutions to photochemistry with complementary models where photo-shielding is accounted for with locally computed prescriptions. With some exceptions, these local models for the radiative shielding length perform reasonably well at reproducing the distribution and amount of molecular gas as compared with a detailed, global ray tracing calculation. Specifically, an approach based on the Jeans Length with a T=40K temperature cap performs the best in regards to a number of different quantitative measures based on the H2 and CO abundances.

Dynamical expansion of H II regions around star clusters plays a key role in dispersing the surrounding dense gas and therefore in limiting the efficiency of star formation in molecular clouds. We use a semi-analytic method and numerical simulations to explore expansion of spherical dusty H II regions and surrounding neutral shells and the resulting cloud disruption. Our model for shell expansion adopts the static solutions of Draine (2011) for dusty H II regions and considers the contact outward forces on the shell due to radiation and thermal pressures as well as the inward gravity from the central star and the shell itself. We show that the internal structure we adopt and the shell evolution from the semi-analytic approach are in good agreement with the results of numerical simulations. Strong radiation pressure in the interior controls the shell expansion indirectly by enhancing the density and pressure at the ionization front. We calculate the minimum star formation efficiency $\epsilon_{min}$ required for cloud disruption as a function of the cloud's total mass and mean surface density. Within the adopted spherical geometry, we find that typical giant molecular clouds in normal disk galaxies have $\epsilon_{min} \lesssim 10$%, with comparable gas and radiation pressure effects on shell expansion. Massive cluster-forming clumps require a significantly higher efficiency of $\epsilon_{min} \gtrsim 50$% for disruption, produced mainly by radiation-driven expansion. The disruption time is typically of the order of a free-fall timescale, suggesting that the cloud disruption occurs rapidly once a sufficiently luminous H II region is formed. We also discuss limitations of the spherical idealization.

High surface density, rapidly star-forming galaxies are observed to have $\approx 50-100\,{\rm km\,s^{-1}}$ line-of-sight velocity dispersions, which are much higher than expected from supernova driving alone, but may arise from large-scale gravitational instabilities. Using three-dimensional simulations of local regions of the interstellar medium, we explore the impact of high velocity dispersions that arise from these disk instabilities. Parametrizing disks by their surface densities and epicyclic frequencies, we conduct a series of simulations that probe a broad range of conditions. Turbulence is driven purely horizontally and on large scales, neglecting any energy input from supernovae. We find that such motions lead to strong global outflows in the highly-compact disks that were common at high redshifts, but weak or negligible mass loss in the more diffuse disks that are prevalent today. Substantial outflows are generated if the one-dimensional horizontal velocity dispersion exceeds $\approx 35\,{\rm km\,s^{-1}},$ as occurs in the dense disks that have star formation rate densities above $\approx 0.1\,{\rm M}_\odot\,{\rm yr}^{-1}\,{\rm kpc}^{-2}.$ These outflows are triggered by a thermal runaway, arising from the inefficient cooling of hot material coupled with successive heating from turbulent driving. Thus, even in the absence of stellar feedback, a critical value of the star-formation rate density for outflow generation can arise due to a turbulent heating instability. This suggests that in strongly self-gravitating disks, outflows may be enhanced by, but need not caused by, energy input from supernovae.

We use three-dimensional magnetohydrodynamic (MHD) simulations to investigate the quasi-equilibrium states of galactic disks regulated by star formation feedback. We incorporate effects from massive-star feedback via time-varying heating rates and supernova (SN) explosions. We find that the disks in our simulations rapidly approach a quasi-steady state that satisfies vertical dynamical equilibrium. The star formation rate (SFR) surface density self-adjusts to provide the total momentum flux (pressure) in the vertical direction that matches the weight of the gas. We quantify feedback efficiency by measuring feedback yields, \eta_c\equiv P_c/\Sigma_SFR (in suitable units), for each pressure component. The turbulent and thermal feedback yields are the same for HD and MHD simulations, \eta_th~1 and \eta_ turb~4, consistent with the theoretical expectations. In MHD simulations, turbulent magnetic fields are rapidly generated by turbulence, and saturate at a level corresponding to \eta_mag,t~1. The presence of magnetic fields enhances the total feedback yield and therefore reduces the SFR, since the same vertical support can be supplied at a smaller SFR. We suggest further numerical calibrations and observational tests in terms of the feedback yields.

Recent hydrodynamic (HD) simulations have shown that galactic disks evolve to reach well-defined statistical equilibrium states. The star formation rate (SFR) self-regulates until energy injection by star formation feedback balances dissipation and cooling in the interstellar medium (ISM), and provides vertical pressure support to balance gravity. In this paper, we extend our previous models to allow for a range of initial magnetic field strengths and configurations, utilizing three-dimensional, magnetohydrodynamic (MHD) simulations. We show that a quasi-steady equilibrium state is established as rapidly for MHD as for HD models unless the initial magnetic field is very strong or very weak, which requires more time to reach saturation. Remarkably, models with initial magnetic energy varying by two orders of magnitude approach the same asymptotic state. In the fully saturated state of the fiducial model, the integrated energy proportions E_kin:E_th:E_mag,t:E_mag,o are 0.35:0.39:0.15:0.11, while the proportions of midplane support P_turb:P_th:\Pi_mag,t:\Pi_mag,o are 0.49:0.18:0.18:0.15. Vertical profiles of total effective pressure satisfy vertical dynamical equilibrium with the total gas weight at all heights. We measure the "feedback yields" \eta_c=P_c/\Sigma_SFR (in suitable units) for each pressure component, finding that \eta_turb~4 and \eta_th~1 are the same for MHD as in previous HD simulations, and \eta_mag,t~1. These yields can be used to predict the equilibrium SFR for a local region in a galaxy based on its observed gas and stellar surface densities and velocity dispersions. As the ISM weight (or dynamical equilibrium pressure) is fixed, an increase in $\eta$ from turbulent magnetic fields reduces the predicted \Sigma_SFR by ~25% relative to the HD case.

In giant molecular clouds (GMCs), shocks driven by converging turbulent flows create high-density, strongly-magnetized regions that are locally sheetlike. In previous work, we showed that within these layers, dense filaments and embedded self-gravitating cores form by gathering material along the magnetic field lines. Here, we extend the parameter space of our three-dimensional, turbulent MHD core formation simulations. We confirm the anisotropic core formation model we previously proposed, and quantify the dependence of median core properties on the pre-shock inflow velocity and upstream magnetic field strength. Our results suggest that bound core properties are set by the total dynamic pressure (dominated by large-scale turbulence) and thermal sound speed c_s in GMCs, independent of magnetic field strength. For models with Mach number between 5 and 20, the median core masses and radii are comparable to the critical Bonnor-Ebert mass and radius defined using the dynamic pressure for P_ext. Our results correspond to M_core = 1.2 c_s^4/sqrt(G^3 rho_0 v_0^2) and R_core = 0.34 c_s^2/sqrt(G rho_0 v_0^2) for rho_0 and v_0 the large-scale mean density and velocity. For our parameter range, the median M_core ~ 0.1-1 M_sun, but a very high pressure cloud could have lower characteristic core mass. We find cores and filaments form simultaneously, and filament column densities are a factor ~2 greater than the surrounding cloud when cores first collapse. We also show that cores identified in our simulations have physical properties comparable to those observed in the Perseus cloud. Superthermal cores in our models are generally also magnetically supercritical, suggesting that the same may be true in observed clouds.

Radiation feedback from young star clusters embedded in giant molecular clouds (GMCs) is believed to be important to the control of star formation. For the most massive and dense clouds, including those in which super star clusters (SSCs) are born, pressure from reprocessed radiation exerted on dust grains may disperse a significant portion of the cloud mass back into the interstellar medium (ISM). Using our radiaton hydrodynamics (RHD) code, Hyperion, we conduct a series of numerical simulations to test this idea. Our models follow the evolution of self-gravitating, strongly turbulent clouds in which collapsing regions are replaced by radiating sink particles representing stellar clusters. We evaluate the dependence of the star formation efficiency (SFE) on the size and mass of the cloud and $\kappa$, the opacity of the gas to infrared (IR) radiation. We find that the single most important parameter determining the evolutionary outcome is $\kappa$, with $\kappa \gtrsim 15 \text{ cm}^2 \text{ g}^{-1}$ needed to disrupt clouds. For $\kappa = 20-40 \text{ cm}^2 \text{ g}^{-1}$, the resulting SFE=50-70% is similar to empirical estimates for some SSC-forming clouds. The opacities required for GMC disruption likely apply only in dust-enriched environments. We find that the subgrid model approach of boosting the direct radiation force $L/c$ by a "trapping factor" equal to a cloud's mean IR optical depth can overestimate the true radiation force by factors of $\sim 4-5$. We conclude that feedback from reprocessed IR radiation alone is unlikely to significantly reduce star formation within GMCs unless their dust abundances or cluster light-to-mass ratios are enhanced.

We investigate prestellar core formation and accretion based on three-dimensional hydrodynamic simulations. Our simulations represent local $\sim 1$pc regions within giant molecular clouds where a supersonic turbulent flow converges, triggering star formation in the post-shock layer. We include turbulence and self-gravity, applying sink particle techniques, and explore a range of inflow Mach number ${\cal M}=2-16$. Two sets of cores are identified and compared: $t_1$-cores are identified of a time snapshot in each simulation, representing dense structures in a single cloud map; $t_\mathrm{coll}$-cores are identified at their individual time of collapse, representing the initial mass reservoir for accretion. We find that cores and filaments form and evolve at the same time. At the stage of core collapse, there is a well-defined, converged characteristic mass for isothermal fragmentation that is comparable to the critical Bonner-Ebert mass at the post-shock pressure. The core mass functions (CMFs) of $t_\mathrm{coll}$-cores show a deficit of high-mass cores ($\gtrsim 7M_\odot$) compared to the observed stellar initial mass function (IMF). However, the CMFs of $t_1$-cores are similar to the observed CMFs and include many low-mass cores that are gravitationally stable. The difference between $t_1$-cores and $t_\mathrm{coll}$-cores suggests that the full sample from observed CMFs may not evolve into protostars. Individual sink particles accrete at a roughly constant rate throughout the simulations, gaining one $t_\mathrm{coll}$-core mass per free-fall time even after the initial mass reservoir is accreted. High-mass sinks gain proportionally more mass at late times than low-mass sinks. There are outbursts in accretion rates, resulting from clumpy density structures falling into the sinks.

Supernova (SN) explosions deposit prodigious energy and momentum in their environments, with the former regulating multiphase thermal structure and the latter regulating turbulence and star formation rates in the interstellar medium (ISM). However, systematic studies quantifying the impact of SNe in realistic inhomogeneous ISM conditions have been lacking. Using three-dimensional hydrodynamic simulations, we investigate the dependence of radial momentum injection on both physical conditions (considering a range of mean density n=0.1-100) and numerical parameters. Our inhomogeneous simulations adopt two-phase background states that result from thermal instability in atomic gas. Although the SNR morphology becomes highly complex for inhomogeneous backgrounds, the radial momentum injection is remarkably insensitive to environmental details. For our two-phase simulations, the final momentum produced by a single SN is given by 2.8*10^5 M_sun*km/s n^{-0.17}. This is only 5% less than the momentum injection for a homogeneous environment with the same mean density, and only 30% greater than the momentum at the time of shell formation. The maximum mass in hot gas is quite insensitive to environmental inhomogeneity. Strong magnetic fields alter the hot gas mass at very late times, but the momentum injection remains the same. Initial experiments with multiple spatially-correlated SNe show a momentum per event nearly as large as single-SN cases. We also present a full numerical parameter study to assess convergence requirements. For convergence in the momentum and other quantities, we find that the numerical resolution dx and the initial size of the SNR r_init must satisfy dx, r_init<r_sf/3, where the shell formation radius is given by r_sf = 30 pc n^{-0.46} for two-phase models (or 30% smaller for a homogeneous medium).

We use ALMA to derive the mass, length, and time scales associated with the nuclear starburst in NGC 253. This region forms ~2 M_sun/yr of stars and resembles other starbursts in scaling relations, with star formation consuming the gas reservoir 10 times faster than in galaxy disks. We present observations of CO, the high effective density transitions HCN(1-0), HCO+(1-0), CS(2-1), and their isotopologues. We identify ten clouds that appear as peaks in line emission and enhancements in the HCN-to-CO ratio. These clouds are massive (~10^7 M_sun) structures with sizes (~30 pc) similar to GMCs in other systems. Compared to disk galaxy GMCs, they show high line widths (~20-40 km/s) given their size, with implied Mach numbers ~90. The clouds also show high surface (~6,000 M_sun/pc^2) and volume densities (n_H2~2,000 cm^-3). Given these, self-gravity can explain the line widths. This short free fall time (~0.7 Myr) helps explain the more efficient star formation in NGC 253. We also consider the starburst region as a whole. The geometry is confused by the high inclination, but simple models support a non-axisymmetric, bar-like geometry with a compact, clumpy region of high gas density embedded in an extended CO distribution. Even for the whole region, the surface density still exceeds that of a disk galaxy GMC. The orbital time (~10 Myr), disk free fall time (<~ 3 Myr), and disk crossing time (<~ 3 Myr) are each much shorter than in a normal spiral galaxy disk. Some but not all aspects of the structure correspond to predictions from assuming vertical dynamical equilibrium or a marginally stable rotating disk. Finally, the CO-to-H2 conversion factor implied by our cloud calculations is approximately Galactic, contrasting with results showing a low value for the whole starburst region. The contrast provides resolved support for the idea of mixed molecular ISM phases in starburst galaxies.

We present observations of N2H+(1-0), HCO+(1-0), and HCN(1-0) toward the Serpens Main molecular cloud from the CARMA Large Area Star Formation Survey (CLASSy). We mapped 150 square arcminutes of Serpens Main with an angular resolution of 7 arcsecs. The gas emission is concentrated in two subclusters (the NW and SE subclusters). The SE subcluster has more prominent filamentary structures and more complicated kinematics compared to the NW subcluster. The majority of gas in the two subclusters has subsonic to sonic velocity dispersions. We applied a dendrogram technique with N2H+(1-0) to study the gas structures; the SE subcluster has a higher degree of hierarchy than the NW subcluster. Combining the dendrogram and line fitting analyses reveals two distinct relations: a flat relation between nonthermal velocity dispersion and size, and a positive correlation between variation in velocity centroids and size. The two relations imply a characteristic depth of 0.15 pc for the cloud. Furthermore, we have identified six filaments in the SE subcluster. These filaments have lengths of 0.2 pc and widths of 0.03 pc, which is smaller than a characteristic width of 0.1 pc suggested by Herschel observations. The filaments can be classified into two types based on their properties. The first type, located in the northeast of the SE subcluster, has larger velocity gradients, smaller masses, and nearly critical mass-per-unit-length ratios. The other type, located in the southwest of the SE subcluster, has the opposite properties. Several YSOs are formed along two filaments which have supercritical mass per unit length ratios, while filaments with nearly critical mass-per-unit-length ratios are not associated with YSOs, suggesting that stars are formed on gravitationally unstable filaments.

We present details of the CARMA Large Area Star Formation Survey (CLASSy), while focusing on observations of Barnard 1. CLASSy is a CARMA Key Project that spectrally imaged N2H+, HCO+, and HCN (J=1-0 transitions) across over 800 square arcminutes of the Perseus and Serpens Molecular Clouds. The observations have angular resolution near 7" and spectral resolution near 0.16 km/s. We imaged ~150 square arcminutes of Barnard 1, focusing on the main core, and the B1 Ridge and clumps to its southwest. N2H+ shows the strongest emission, with morphology similar to cool dust in the region, while HCO+ and HCN trace several molecular outflows from a collection of protostars in the main core. We identify a range of kinematic complexity, with N2H+ velocity dispersions ranging from ~0.05-0.50 km/s across the field. Simultaneous continuum mapping at 3 mm reveals six compact object detections, three of which are new detections. A new non-binary dendrogram algorithm is used to analyze dense gas structures in the N2H+ position-position-velocity (PPV) cube. The projected sizes of dendrogram-identified structures range from about 0.01-0.34 pc. Size-linewidth relations using those structures show that non-thermal line-of-sight velocity dispersion varies weakly with projected size, while rms variation in the centroid velocity rises steeply with projected size. Comparing these relations, we propose that all dense gas structures in Barnard 1 have comparable depths into the sky, around 0.1-0.2 pc; this suggests that over-dense, parsec-scale regions within molecular clouds are better described as flattened structures rather than spherical collections of gas. Science-ready PPV cubes for Barnard 1 molecular emission are available for download.

We present the N2H+(J=1-0) map of the Serpens South molecular cloud obtained as part of the CARMA Large Area Star Formation Survey (CLASSy). The observations cover 250 square arcminutes and fully sample structures from 3000 AU to 3 pc with a velocity resolution of 0.16 km/s, and they can be used to constrain the origin and evolution of molecular cloud filaments. The spatial distribution of the N2H+ emission is characterized by long filaments that resemble those observed in the dust continuum emission by Herschel. However, the gas filaments are typically narrower such that, in some cases, two or three quasi-parallel N2H+ filaments comprise a single observed dust continuum filament. The difference between the dust and gas filament widths casts doubt on Herschel ability to resolve the Serpens South filaments. Some molecular filaments show velocity gradients along their major axis, and two are characterized by a steep velocity gradient in the direction perpendicular to the filament axis. The observed velocity gradient along one of these filaments was previously postulated as evidence for mass infall toward the central cluster, but these kind of gradients can be interpreted as projection of large-scale turbulence.

We use three-dimensional numerical hydrodynamic simulations of the turbulent, multiphase atomic interstellar medium (ISM) to construct and analyze synthetic HI 21 cm emission and absorption lines. Our analysis provides detailed tests of 21 cm observables as physical diagnostics of the atomic ISM. In particular, we construct (1) the "observed" spin temperature, $T_{s,obs}(v_{ch})\equiv T_B(v_{ch})/[1-e^{-{\tau}(v_{ch})}]$, and its optical-depth weighted mean T_s,obs; (2) the absorption-corrected "observed" column density, $N_{H,obs}\propto \int dv_{ch} T_B(v_{ch}){\tau}(v_{ch})/[[1-e^{-{\tau}(v_{ch})}]$; and (3) the "observed" fraction of cold neutral medium (CNM), $f_{c,obs}\equiv T_c/T_{s,obs}$ for T_c the CNM temperature; we compare each observed parameter with true values obtained from line-of-sight (LOS) averages in the simulation. Within individual velocity channels, T_s,obs(v_ch) is within a factor 1.5 of the true value up to ${\tau}(v_{ch})\approx10$. As a consequence, N_H,obs and T_s,obs are respectively within 5% and 12% of the true values for 90% and 99% of LOSs. The optically thin approximation significantly underestimates N_H for ${\tau}>1$. Provided that T_c is constrained, an accurate observational estimate of the CNM mass fraction can be obtained down to 20%. We show that T_s,obs cannot be used to distinguish the relative proportions of warm and thermally-unstable atomic gas, although the presence of thermally-unstable gas can be discerned from 21 cm lines with 200K<$T_{s,obs}(v_{ch})$<1000K. Our mock observations successfully reproduce and explain the observed distribution of the brightness temperature, optical depth, and spin temperature in Roy et al. (2013a). The threshold column density for CNM seen in observations is also reproduced by our mock observations. We explain this observed threshold behavior in terms of vertical equilibrium in the local Milky Way's ISM disk.

We investigate the roles of magnetic fields and ambipolar diffusion during prestellar core formation in turbulent giant molecular clouds (GMCs), using three-dimensional numerical simulations. Our simulations focus on the shocked layer produced by a converging flow within a GMC, and survey varying ionization and angle between the upstream flow and magnetic field. We also include ideal magnetohydrodynamic (MHD) and hydrodynamic models. From our simulations, we identify hundreds of self-gravitating cores that form within 1 Myr, with masses M ~ 0.04 - 2.5 solar-mass and sizes L ~ 0.015 - 0.07 pc, consistent with observations of the peak of the core mass function (CMF). Median values are M = 0.47 solar-mass and L = 0.03 pc. Core masses and sizes do not depend on either the ionization or upstream magnetic field direction. In contrast, the mass-to-magnetic flux ratio does increase with lower ionization, from twice to four times the critical value. The higher mass-to-flux ratio for low ionization is the result of enhanced transient ambipolar diffusion when the shocked layer first forms. However, ambipolar diffusion is not necessary to form low-mass supercritical cores. For ideal MHD, we find similar masses to other cases. These masses are 1 - 2 orders of magnitude lower than the value that defines a magnetically supercritical sphere under post-shock ambient conditions. This discrepancy is the result of anisotropic contraction along field lines, which is clearly evident in both ideal MHD and diffusive simulations. We interpret our numerical findings using a simple scaling argument which suggests that gravitationally critical core masses will depend on the sound speed and mean turbulent pressure in a cloud, regardless of magnetic effects.

Giant molecular clouds (GMCs) are the primary reservoirs of cold, star-forming molecular gas in the Milky Way and similar galaxies, and thus any understanding of star formation must encompass a model for GMC formation, evolution, and destruction. These models are necessarily constrained by measurements of interstellar molecular and atomic gas, and the emergent, newborn stars. Both observations and theory have undergone great advances in recent years, the latter driven largely by improved numerical simulations, and the former by the advent of large-scale surveys with new telescopes and instruments. This chapter offers a thorough review of the current state of the field.

The energy and momentum feedback from young stars has a profound impact on the interstellar medium (ISM), including heating and driving turbulence in the neutral gas that fuels future star formation. Recent theory has argued that this leads to a quasi-equilibrium self-regulated state, and for outer atomic-dominated disks results in the surface density of star formation $\Sigma_{SFR}$ varying approximately linearly with the weight of the ISM (or midplane turbulent + thermal pressure). We use three-dimensional numerical hydrodynamic simulations to test the theoretical predictions for thermal, turbulent, and vertical dynamical equilibrium, and the implied functional dependence of $\Sigma_{SFR}$ on local disk properties. Our models demonstrate that all equilibria are established rapidly, and that the expected proportionalities between mean thermal and turbulent pressures and $\Sigma_{SFR}$ apply. For outer disk regions, this results in $\Sigma_{SFR} \propto \Sigma \sqrt{\rho_{sd}}$, where $\Sigma$ is the total gas surface density and $\rho_{sd}$ is the midplane density of the stellar disk (plus dark matter). This scaling law arises because $\rho_{sd}$ sets the vertical dynamical time in our models (and outer disk regions generally). The coefficient in the star formation law varies inversely with the specific energy and momentum yield from massive stars. We find proportions of warm and cold atomic gas, turbulent-to-thermal pressure, and mean velocity dispersions that are consistent with Solar-neighborhood and other outer-disk observations. This study confirms the conclusions of a previous set of simulations, which incorporated the same physics treatment but was restricted to radial-vertical slices through the ISM.

The under-abundance of very massive galaxies in the universe is frequently attributed to the effect of galactic winds. Although ionized galactic winds are readily observable most of the expelled mass is likely in cooler atomic and molecular phases. Expanding molecular shells observed in starburst systems such as NGC 253 and M 82 may facilitate the entrainment of molecular gas in the wind. While shell properties are well constrained, determining the amount of outflowing gas emerging from such shells and the connection between this gas and the ionized wind requires spatial resolution <100 pc coupled with sensitivity to a wide range of spatial scales, hitherto not available. Here we report observations of NGC 253, a nearby starburst galaxy (D~3.4 Mpc) known to possess a wind, which trace the cool molecular wind at 50 pc resolution. At this resolution the extraplanar molecular gas closely tracks the H{\alpha} filaments, and it appears connected to molecular expanding shells located in the starburst region. These observations allow us to directly measure the molecular outflow rate to be > 3 Msun/yr and likely ~9 Msun/yr. This implies a ratio of mass-outflow rate to star formation rate of at least {\eta}~1-3, establishing the importance of the starburst-driven wind in limiting the star formation activity and the final stellar content.

We describe a module for the Athena code that solves the gray equations of radiation hydrodynamics (RHD), based on the first two moments of the radiative transfer equation. We use a combination of explicit Godunov methods to advance the gas and radiation variables including the non-stiff source terms, and a local implicit method to integrate the stiff source terms. We adopt the M1 closure relation and include all leading source terms. We employ the reduced speed of light approximation (RSLA) with subcycling of the radiation variables in order to reduce computational costs. Our code is dimensionally unsplit in one, two, and three space dimensions and is parallelized using MPI. The streaming and diffusion limits are well-described by the M1 closure model, and our implementation shows excellent behavior for a problem with a concentrated radiation source containing both regimes simultaneously. Our operator-split method is ideally suited for problems with a slowly varying radiation field and dynamical gas flows, in which the effect of the RSLA is minimal. We present an analysis of the dispersion relation of RHD linear waves highlighting the conditions of applicability for the RSLA. To demonstrate the accuracy of our method, we utilize a suite of radiation and RHD tests covering a broad range of regimes, including RHD waves, shocks, and equilibria, which show second-order convergence in most cases. As an application, we investigate radiation-driven ejection of a dusty, optically thick shell in the interstellar medium (ISM). Finally, we compare the timing of our method with other well-known iterative schemes for the RHD equations. Our code implementation, Hyperion, is suitable for a wide variety of astrophysical applications and will be made freely available on the Web.

We present the Evolution of molecular Gas in Normal Galaxies (EGNoG) survey, an observational study of molecular gas in 31 star-forming galaxies from z=0.05 to z=0.5, with stellar masses of (4-30)x10^10 M_Sun and star formation rates of 4-100 M_Sun yr^-1. This survey probes a relatively un-observed redshift range in which the molecular gas content of galaxies is expected to have evolved significantly. To trace the molecular gas in the EGNoG galaxies, we observe the CO(1-0) and CO(3-2) rotational lines using the Combined Array for Research in Millimeter-wave Astronomy (CARMA). We detect 24 of 31 galaxies and present resolved maps of 10 galaxies in the lower redshift portion of the survey. We use a bimodal prescription for the CO to molecular gas conversion factor, based on specific star formation rate, and compare the EGNoG galaxies to a large sample of galaxies assembled from the literature. We find an average molecular gas depletion time of 0.76 \pm 0.54 Gyr for normal galaxies and 0.06 \pm 0.04 Gyr for starburst galaxies. We calculate an average molecular gas fraction of 7-20% at the intermediate redshifts probed by the EGNoG survey. By expressing the molecular gas fraction in terms of the specific star formation rate and molecular gas depletion time (using typical values), we also calculate the expected evolution of the molecular gas fraction with redshift. The predicted behavior agrees well with the significant evolution observed from z~2.5 to today.

We describe implementation and tests of sink particle algorithms in the Eulerian grid-based code Athena. Introduction of sink particles enables long-term evolution of systems in which localized collapse occurs, and it is impractical (or unnecessary) to resolve the accretion shocks at the centers of collapsing regions. We discuss similarities and differences of our methods compared to other implementations of sink particles. Our criteria for sink creation are motivated by the properties of the Larson-Penston collapse solution. We use standard particle-mesh methods to compute particle and gas gravity together. Accretion of mass and momenta onto sinks is computed using fluxes returned by the Riemann solver. A series of tests based on previous analytic and numerical collapse solutions is used to validate our method and implementation. We demonstrate use of our code for applications with a simulation of planar converging supersonic turbulent flow, in which multiple cores form and collapse to create sinks; these sinks continue to interact and accrete from their surroundings over several Myr.

Star formation is self-regulated by its feedback that drives turbulence and heats the gas. In equilibrium, the star formation rate (SFR) should be directly related to the total (thermal plus turbulent) midplane pressure and hence the total weight of the diffuse gas if energy balance and vertical dynamical equilibrium hold simultaneously. To investigate this quantitatively, we utilize numerical hydrodynamic simulations focused on outer-disk regions where diffuse atomic gas dominates. By analyzing gas properties at saturation, we obtain relationships between the turbulence driving and dissipation rates, heating and cooling rates, the total midplane pressure and the total weight of gas, and the SFR and the total midplane pressure. We find a nearly linear relationship between the SFR and the midplane pressure consistent with the theoretical prediction.

We explore the self-regulation of star formation using a large suite of high resolution hydrodynamic simulations, focusing on molecule-dominated regions (galactic centers and [U]LIRGS) where feedback from star formation drives highly supersonic turbulence. In equilibrium the total midplane pressure, dominated by turbulence, must balance the vertical weight of the ISM. Under self-regulation, the momentum flux injected by feedback evolves until it matches the vertical weight. We test this flux balance in simulations spanning a range of parameters, including surface density $\Sigma$, momentum injected per stellar mass formed ($p_*/m_*$), and angular velocity. The simulations are 2D radial-vertical slices, including both self-gravity and an external potential that confines gas to the disk midplane. After the simulations reach a steady state in all relevant quantities, including the star formation rate $\Sigma_{SFR}$, there is remarkably good agreement between the vertical weight, the turbulent pressure, and the momentum injection rate from supernovae. Gas velocity dispersions and disk thicknesses increase with $p_*/m_*$. The efficiency of star formation per free-fall time at the mid-plane density is insensitive to the local conditions and to the star formation prescription in very dense gas. We measure efficiencies $\sim$0.004-0.01, consistent with low and approximately constant efficiencies inferred from observations. For $\Sigma\in$(100--1000) \msunpc, we find $\Sigma_{SFR}\in$(0.1--4) \sfrunits, generally following a $\Sigma_{SFR}\propto \Sigma^2$ relationship. The measured relationships agree very well with vertical equilibrium and with turbulent energy replenishment by feedback within a vertical crossing time. These results, along with the observed $\Sigma_{SFR}-\Sigma$ relation in high density environments, provide strong evidence for the self-regulation of star formation.

The most common means of converting an observed CO line intensity into a molecular gas mass requires the use of a conversion factor (Xco). While in the Milky Way this quantity does not appear to vary significantly, there is good reason to believe that Xco will depend on the larger-scale galactic environment. Utilising numerical models, we investigate how varying metallicities, gas temperatures and velocity dispersions in galaxies impact the way CO line emission traces the underlying H2 gas mass, and under what circumstances Xco may differ from the Galactic mean value. We find that, due to the combined effects of increased gas temperature and velocity dispersion, Xco is depressed below the Galactic mean in high surface density environments such as ULIRGs. In contrast, in low metallicity environments, Xco tends to be higher than in the Milky Way, due to photodissociation of CO in metal-poor clouds. At higher redshifts, gas-rich discs may have gravitationally unstable clumps which are warm (due to increased star formation) and have elevated velocity dispersions. These discs tend to have Xco values ranging between present-epoch gas-rich mergers and quiescent discs at low-z. This model shows that on average, mergers do have lower Xco values than disc galaxies, though there is significant overlap. Xco varies smoothly with the local conditions within a galaxy, and is not a function of global galaxy morphology. We combine our results to provide a general fitting formula for Xco as a function of CO line intensity and metallicity. We show that replacing the traditional approach of using one constant Xco for starbursts and another for discs with our best-fit function produces star formation laws that are continuous rather than bimodal, and that have significantly reduced scatter.

We analyze the properties of steady and time-dependent C shocks under conditions prevailing in giant molecular clouds. For steady C shocks, we show that ionization equilibrium holds and use numerical integrations to obtain the shock thickness mediated by ambipolar diffusion as a function of density, inflow velocity, magnetic field, and ionization fraction in the cloud. Our formula also agrees with an analytic estimate based on ion-neutral momentum exchange. Using time-dependent numerical simulations, we show that C shocks have a transient stage when the neutrals are compressed much more strongly than the magnetic field. The transient stage has a duration set by the neutral-ion collision time, t_AD ~ L_shock/v_drift ~ 0.1 - 1 Myr. This transient creates a strong enhancement in the mass-to-magnetic flux ratio. Under favorable conditions, supercritical prestellar cores may form and collapse promptly as a result of magnetic flux loss during the transient stage of C shocks.

We use vertically-resolved numerical hydrodynamic simulations to study star formation and the interstellar medium (ISM) in galactic disks. We focus on outer disk regions where diffuse HI dominates, with gas surface densities Sigma_SFR=3-20 Msun/kpc^2/yr and star-plus-dark matter volume densities rho_sd=0.003-0.5 Msun/pc^3. Star formation occurs in very dense, cold, self-gravitating clouds. Turbulence, driven by momentum feedback from supernova events, destroys bound clouds and puffs up the disk vertically. Time-dependent radiative heating (FUV) offsets gas cooling. We use our simulations to test a new theory for self-regulated star formation. Consistent with this theory, the disks evolve to a state of vertical dynamical equilibrium and thermal equilibrium with both warm and cold phases. The range of star formation surface densities and midplane thermal pressures is Sigma_SFR ~ 0.0001 - 0.01 Msun/kpc^2/yr and P_th/k_B ~ 100 -10000 cm^-3 K. In agreement with observations, turbulent velocity dispersions are ~7 km/s and the ratio of the total (effective) to thermal pressure is P_tot/P_th~4-5, across this whole range. We show that Sigma_SFR is not well correlated with Sigma alone, but rather with Sigma*(rho_sd)^1/2, because the vertical gravity from stars and dark matter dominates in outer disks. We also find that Sigma_SFR has a strong, nearly linear correlation with P_tot, which itself is within ~13% of the dynamical-equilibrium estimate P_tot,DE. The quantitative relationships we find between Sigma_SFR and the turbulent and thermal pressures show that star formation is highly efficient for energy and momentum production, in contrast to the low efficiency of mass consumption. Star formation rates adjust until the ISM's energy and momentum losses are replenished by feedback within a dynamical time.

Relating the observed CO emission from giant molecular clouds (GMCs) to the underlying H2 column density is a long-standing problem in astrophysics. While the Galactic CO-H2 conversion factor (Xco) appears to be reasonably constant, observations indicate that Xco may be depressed in high-surface density starburst environments. Using a multi-scale approach, we investigate the dependence of Xco on the galactic environment in numerical simulations of disc galaxies and galaxy mergers. Xco is proportional to the GMC surface density divided by the integrated CO intensity, Wco, and Wco is related to the kinetic temperature and velocity dispersion in the cloud. In disc galaxies (except within the central ~ kpc), the galactic environment is largely unimportant in setting the physical properties of GMCs provided they are gravitationally bound. The temperatures are roughly constant at ~10 K due to the balance of CO cooling and cosmic ray heating, giving a nearly constant CO-H2 conversion factor in discs. In mergers, the velocity dispersion of the gas rises dramatically during coalescence. The gas temperature also rises as it couples well to the warm (~50 K) dust at high densities (n > 10^4 cm^-3). The rise in velocity dispersion and temperature combine to offset the rise in surface density in mergers, causing Xco to drop by a factor of ~2-10 compared to the disc simulation. This model predicts that high-resolution ALMA observations of nearby ULIRGs should show velocity dispersions of ~10-100 km/s, and brightness temperatures comparable to the dust temperatures.

We compare atomic gas, molecular gas, and the recent star formation rate (SFR) inferred from H-alpha in the Small Magellanic Cloud (SMC). By using infrared dust emission and local dust-to-gas ratios, we construct a map of molecular gas that is independent of CO emission. This allows us to disentangle conversion factor effects from the impact of metallicity on the formation and star formation efficiency of molecular gas. On scales of 200 pc to 1 kpc we find a characteristic molecular gas depletion time of ~1.6 Gyr, similar to that observed in the molecule-rich parts of large spiral galaxies on similar spatial scales. This depletion time shortens on much larger scales to ~0.6 Gyr because of the presence of a diffuse H-alpha component, and lengthens on much smaller scales to ~7.5 Gyr because the H-alpha and H2 distributions differ in detail. We estimate the systematic uncertainties in our measurement to be a factor of 2-3. We suggest that the impact of metallicity on the physics of star formation in molecular gas has at most this magnitude. The relation between SFR and neutral (H2+HI) gas surface density is steep, with a power-law index ~2.2+/-0.1, similar to that observed in the outer disks of large spiral galaxies. At a fixed total gas surface density the SMC has a 5-10 times lower molecular gas fraction (and star formation rate) than large spiral galaxies. We explore the ability of the recent models by Krumholz et al. (2009) and Ostriker et al. (2010) to reproduce our observations. We find that to explain our data at all spatial scales requires a low fraction of cold, gravitationally-bound gas in the SMC. We explore a combined model that incorporates both large scale thermal and dynamical equilibrium and cloud-scale photodissociation region structure and find that it reproduces our data well, as well as predicting a fraction of cold atomic gas very similar to that observed in the SMC.

We investigate how the X factor, the ratio of H_2 column density (NH2) to velocity-integrated CO intensity (W), is determined by the physical properties of gas in model molecular clouds (MCs). We perform radiative transfer calculations on chemical-MHD models to compute X. Using integrated NH2 and W reproduces the limited range in X found in observations, resulting in a mean value X=2\times10^20 s/cm^2/K^1/km^1 from the Galactic MC model. However, in limited velocity intervals, X can take on a much larger range due to CO line saturation. Thus, X strongly depends on both the range in gas velocities and volume densities. The temperature (T) variations within individual MCs do not strongly affect X, as dense gas contributes most to setting X. For fixed velocity and density structure, gas with higher T has higher W, yielding X ~ T^-1/2 for T~20-100 K. We demonstrate that the linewidth-size scaling relation does not influence the X factor - only the range in velocities is important. Clouds with larger linewidths, regardless of the linewidth-size relation, have a higher W, corresponding to a lower value of X, scaling roughly as X ~ sigma^-1/2. The "mist" model, consisting of optically thick cloudlets with well-separated velocities, does not accurately reflect the conditions in a turbulent MC. We propose that the observed cloud-average values of X ~ XGal is simply a result of the limited range in NH2, temperatures, and velocities found in Galactic MCs - a ~constant value of X therefore does not require any linewidth-size relation, or that MCs are virialized objects. Since gas properties likely differ (slightly) between clouds, masses derived through a standard X should only be considered as a rough first estimate. For temperatures T~10-20 K, velocity dispersions ~1-6 km/s, and NH2~2-20\times10^21 cm^-2, we find cloud-averaged X ~ 2-4\times10^20 s/cm^2/K^1/km^1 for Solar-metallicity models.

Star formation rates in the centers of disk galaxies often vastly exceed those at larger radii. We investigate the idea that these central starbursts are self-regulated, with the momentum flux injected to the ISM by star formation balancing the gravitational force confining the gas. For most starbursts, supernovae are the largest contributor to the momentum flux, and turbulence provides the main pressure support for the predominantly-molecular ISM. If the momentum feedback per stellar mass formed is p_*/m_* ~ 3000 km/s, the predicted star formation rate is Sigma_SFR=2 pi G Sigma^2 m_*/p_* ~0.1(Sigma/100Msun/pc^2)^2 Msun/kpc^2/yr in regions where gas dominates the vertical gravity. We compare this prediction with numerical simulations of vertically-resolved disks that model star formation including feedback, finding good agreement for gas surface densities Sigma ~ 10^2-10^3 Msun/pc^2. We also compare to a compilation of star formation rates and gas contents from local and high-redshift galaxies (both mergers and normal galaxies), finding good agreement provided that X_CO decreases weakly as Sigma and Sigma_SFR increase. Star formation rates in dense, turbulent gas are also expected to depend on the gravitational free-fall time; if the efficiency per free-fall time is epsilon_ff ~ 0.01, the turbulent velocity dispersion driven by feedback is expected to be v_z = 0.4 epsilon_ff p_*/m_* ~ 10 km/s, relatively independent of Sigma or Sigma_SFR. Turbulence-regulated starbursts (controlled by kinetic momentum feedback) are part of the larger scheme of self-regulation; primarily-atomic low-Sigma outer disks may have star formation regulated by UV heating feedback, whereas regions at extremely high Sigma may be regulated by feedback of radiation that is reprocessed into trapped IR.

We use numerical hydrodynamic simulations to investigate prestellar core formation in the dynamic environment of giant molecular clouds, focusing on planar post-shock layers produced by colliding turbulent flows. A key goal is to test how core evolution and properties depend on the velocity dispersion in the parent cloud; our simulation suite consists of 180 models with inflow Mach numbers Ma=v/c_s=1.1-9. At all Mach numbers, our models show that turbulence and self-gravity collect gas within post-shock regions into filaments at the same time as overdense areas within these filaments condense into cores. This morphology, together with the subsonic velocities we find inside cores, is similar to observations. We extend previous results showing that core collapse develops in an ``outside-in'' manner, with density and velocity approaching the Larson-Penston asymptotic solution. The time for the first core to collapse varies as 1/sqrt(v), consistent with analytic estimates. Core building takes 10 times as long as core collapse, consistent with observed prestellar core lifetimes. Core shapes change from oblate to prolate as they evolve. To define cores, we use isosurfaces of the gravitational potential. We compare to cores defined using the potential computed from projected surface density, finding good agreement for core masses and sizes; this offers a new way to identify cores in observed maps. Cores with masses varying by three orders of magnitude (0.05 - 50 M_sun) are identified in our simulations. Stability analysis of post-shock layers predicts that the first core to collapse will have mass M \propto v^-1/2, and that the minimum mass for cores formed at late times will have M\propto v^-1. From our simulations, the median mass lies between these two relations.

Star formation depends on the available gaseous "fuel" as well as galactic environment, with higher specific star formation rates where gas is predominantly molecular and where stellar (and dark matter) densities are higher. The partition of gas into different thermal components must itself depend on the star formation rate, since a steady state distribution requires a balance between heating (largely from stellar UV for the atomic component) and cooling. In this presentation, I discuss a simple thermal and dynamical equilibrium model for the star formation rate in disk galaxies, where the basic inputs are the total surface density of gas and the volume density of stars and dark matter, averaged over ~kpc scales. Galactic environment is important because the vertical gravity of the stars and dark matter compress gas toward the midplane, helping to establish the pressure, and hence the cooling rate. In equilibrium, the star formation rate must evolve until the gas heating rate is high enough to balance this cooling rate and maintain the pressure imposed by the local gravitational field. In addition to discussing the formulation of this equilibrium model, I review the current status of numerical simulations of multiphase disks, focusing on measurements of quantities that characterize the mean properties of the diffuse ISM. Based on simulations, turbulence levels in the diffuse ISM appear relatively insensitive to local disk conditions and energetic driving rates, consistent with observations. It remains to be determined, both from observations and simulations, how mass exchange processes control the ratio of cold-to-warm gas in the atomic ISM.

We develop a model for regulation of galactic star formation rates Sigma_SFR in disk galaxies, in which ISM heating by stellar UV plays a key role. By requiring simultaneous thermal and (vertical) dynamical equilibrium in the diffuse gas, and star formation at a rate proportional to the mass of the self-gravitating component, we obtain a prediction for Sigma_SFR as a function of the total gaseous surface density Sigma and the density of stars + dark matter, rho_sd. The physical basis of this relationship is that thermal pressure in the diffuse ISM, which is proportional to the UV heating rate and therefore to Sigma_SFR, must adjust to match the midplane pressure set by the vertical gravitational field. Our model applies to regions where Sigma < 100 Msun/pc^2. In low-Sigma_SFR (outer-galaxy) regions where diffuse gas dominates, the theory predicts Sigma_SFR \propto Sigma (rho_sd)^1/2. The decrease of thermal equilibrium pressure when Sigma_SFR is low implies, consistent with observations, that star formation can extend (with declining efficiency) to large radii in galaxies, rather than having a sharp cutoff. The main parameters entering our model are the ratio of thermal pressure to total pressure in the diffuse ISM, the fraction of diffuse gas that is in the warm phase, and the star formation timescale in self-gravitating clouds; all of these are (in principle) direct observables. At low surface density, our model depends on the ratio of the mean midplane FUV intensity (or thermal pressure in the diffuse gas) to the star formation rate, which we set based on Solar neighborhood values. We compare our results to recent observations, showing good agreement overall for azimuthally-averaged data in a set of spiral galaxies. For the large flocculent spiral galaxies NGC 7331 and NGC 5055, the correspondence between theory and observation is remarkably close.

Galactic spiral shocks are dominant morphological features and believed to be responsible for substructure formation within spiral arms in disk galaxies. They can also contribute a substantial amount of kinetic energy to the interstellar gas by tapping the (differential) rotational motion. We use numerical hydrodynamic simulations to investigate dynamics and structure of spiral shocks with thermal instability in vertically stratified galactic disks, focusing on environmental conditions (of heating and the galactic potential) similar to the Solar neighborhood. We initially consider an isothermal disk in vertical hydrostatic equilibrium and let it evolve subject to interstellar cooling and heating as well as a stellar spiral potential. Due to thermal instability, a disk with surface density $\Sigma_0 \geq 6.7\Surf$ rapidly turns to a thin dense slab near the midplane sandwiched between layers of rarefied gas. The imposed spiral potential leads to a vertically curved shock that exhibits strong flapping motions in the plane perpendicular to the arm. The overall flow structure at saturation is comprised of arm, postshock expansion zone, and interarm regions that occupy typically 10\%, 20\%, and 70\% of the arm-to-arm distance, in which the gas resides for 15\%, 30\%, and 55\% of the arm-to-arm crossing time, respectively. The flows are characterized by transitions from rarefied to dense phases at the shock and from dense to rarefied phases in the postshock expansion zone, although gas with too-large postshock-density does not undergo this return phase transition, instead forming dense condensations. If self-gravity is omitted, the shock flapping drives random motions in the gas, but only up to $\sim 2-3 \kms$ in the in-plane direction and less than $2\kms$ in the vertical direction. Time-averaged shock profiles show... (abridged).

A method for implementing cylindrical coordinates in the Athena magnetohydrodynamics (MHD) code is described. The extension follows the approach of Athena's original developers and has been designed to alter the existing Cartesian-coordinates code as minimally and transparently as possible. The numerical equations in cylindrical coordinates are formulated to maintain consistency with constrained transport, a central feature of the Athena algorithm, while making use of previously implemented code modules such as the Riemann solvers. Angular-momentum transport, which is critical in astrophysical disk systems dominated by rotation, is treated carefully. We describe modifications for cylindrical coordinates of the higher-order spatial reconstruction and characteristic evolution steps as well as the finite-volume and constrained transport updates. Finally, we present a test suite of standard and novel problems in one-, two-, and three-dimensions designed to validate our algorithms and implementation and to be of use to other code developers. The code is suitable for use in a wide variety of astrophysical applications and is freely available for download on the web.

We present a unified model for molecular core formation and evolution, based on numerical simulations of converging, supersonic flows. Our model applies to star formation in GMCs dominated by large-scale turbulence, and contains four main stages: core building, core collapse, envelope infall, and late accretion. During the building stage, cores form out of dense, post-shock gas, and become increasingly centrally stratified as the mass grows over time. When the shock radius defining the core boundary exceeds $R\approx 4 a (4\pi G \rho_{mean})^{-1/2}$, where $a$ is the isothermal sound speed, a wave of collapse propagates from the edge to the center. During the building and collapse stages, density profiles can be fit by Bonnor-Ebert profiles with temperature 1.2 - 2.9 times the true value. As found previously for initially static equilibria, outside-in collapse leads to a Larson-Penston density profile $\rho \approx 8.86 a^2/(4 \pi G r^2)$. The third stage, consisting of an inside-out wave of gravitational rarefaction leading to $\rho\propto r^{-3/2}$, $v\propto r^{-1/2}$, is also similar to that for initially-static spheres, as originally described by Shu. We find that the collapse and infall stages have comparable duration, $\sim t_{ff}$, consistent with estimates for observed prestellar and protostellar (Class 0/I) cores. Core building takes longer, but does not produce high-contrast objects until shortly before collapse. The time to reach core collapse, and the core mass at collapse, decrease with increasing inflow Mach number. For all cases the accretion rate is $\gg a^3/G$ early on but sharply drops off; the final system mass depends on the duration of late-stage accretion, set by large-scale conditions in a cloud.

We use numerical simulations of turbulent, multiphase, self-gravitating gas orbiting in model disk galaxies to study the relationships among pressure, the vertical gas distribution, and the ratio of dense to diffuse gas. We show that the disk height and mean midplane pressure are consistent with effective hydrostatic equilibrium, provided that the turbulent vertical velocity dispersion and gas self-gravity are included. Mass-weighted pressures are an order of magnitude higher than the midplane pressure because self-gravity concentrates gas and increases the pressure in clouds. We also investigate the ratio Rmol=M(H2)/M(HI) for our simulations. Blitz and Rosolowsky (2006) showed that Rmol is proportional to the estimated midplane pressure. For model series in which the epicyclic frequency, kappa, and gas surface density, Sigma, are proportional, we recover the empirical relation. For other model series in which kappa and Sigma are independent, the midplane pressure and Rmol are not well correlated. We conclude that the molecular fraction -- and star formation rate -- of a galactic disk inherently depends on its rotational state, not just the local values of Sigma and the stellar density rho*. The empirical correlation between Rmol and midplane pressure implies that the "environmental parameters" kappa, Sigma, and rho* are interdependent in real galaxies, presumably as a consequence of evolution toward states with Toomre Q near unity. We note that Rmol in static models far exceeds both the values in our turbulent simulations and observed values, implying that turbulence is crucial to obtaining a realistic molecular fraction in the ISM.

Using numerical simulations of galactic disks resolving scales from ~1 to several hundred pc, we investigate dynamical properties of the multiphase ISM with turbulence driven by star formation feedback. We focus on HII region effects by applying intense heating in dense, self-gravitating regions. Our models are two-dimensional radial-vertical slices through the disk, and include sheared background rotation, vertical stratification, heating and cooling to yield temperatures T~10-10^4K, and thermal conduction. We separately vary the gas surface density Sigma, the stellar volume density rho_*, and the local angular rotation rate Omega to explore environmental dependencies, and analyze the steady-state properties of each model. Among other statistics, we evaluate turbulent amplitudes, virial ratios, Toomre Q parameters including turbulence, and the mass fractions at different densities. We find that the dense gas (n>100 cm^-3) has turbulence levels similar to observed GMCs and virial ratios ~1-2. The Toomre Q parameter in dense gas reaches near unity, demonstrating self-regulation via turbulent feedback. We also test how the star formation rate Sigma_SFR depends on Sigma, rho_*, and Omega. Under the assumption that the star formation rate is proportional to the mass at densities above n_th divided by the free-fall time at that threshold, we find that Sigma_SFR varies as Sigma^(1+p) with 1+p ~ 1.2-1.4 when n_th=10^2 or 10^3 cm^-3, consistent with observations. Estimated star formation rates based on large-scale properties (the orbital time, the Jeans time, or the free-fall time at the vertically-averaged density) however depart from rates computed using the dense gas mass, indicating that resolving the ISM structure in galactic disks at scales <<H is necessary for accurate predictions of the star formation rate.

Spiral arm spurs are prominent features that have been observed in extinction and 8$\mu$m emission in nearby galaxies. In order to understand their molecular gas properties, we used the Owens Valley Radio Observatory to map the CO(J=1--0) emission in three spurs emanating from the inner northwestern spiral arm of M51. We report CO detections from all three spurs. The molecular gas mass and surface density are M$_{H2} \sim3\times10^6$ M$_{\sun}$ and $\Sigma_{H2} \sim$50 M$_{\sun}$ pc$^{-2}$. Thus, relative to the spiral arms, the spurs are extremely weak features. However, since the spurs are extended perpendicular to the spiral arms for $\sim$500 pc and contain adequate fuel for star formation, they may be the birthplace for observed inter-arm HII regions. This reduces the requirement for the significant time delay that would be otherwise needed if the inter-arm star formation was initiated in the spiral arms. Larger maps of galaxies at similar depth are required to further understand the formation and evolution of these spurs and their role in star formation - such data should be forthcoming with the new CARMA and future ALMA telescopes and can be compared to several recent numerical simulations that have been examining the evolution of spiral arm spurs.

We include feedback in global hydrodynamic simulations in order to study the star formation properties, and gas structure and dynamics, in models of galactic disks. We extend previous models by implementing feedback in gravitationally bound clouds: momentum is injected at a rate proportional to the star formation rate. This mechanical energy disperses cloud gas back into the surrounding ISM, truncating star formation in a given cloud, and raising the overall level of ambient turbulence. Propagating star formation can however occur as expanding shells collide, enhancing the density and triggering new cloud and star formation. By controlling the momentum injection per massive star and the specific star formation rate in dense gas, we find that the negative effects of high turbulence outweigh the positive ones, and in net feedback reduces the fraction of dense gas and thus the overall star formation rate. The properties of the large clouds that form are not, however, very sensitive to feedback, with cutoff masses of a few million solar masses, similar to observations. We find a relationship between the star formation rate surface density and the gas surface density with a power law index ~2 for our models with the largest dynamic range, consistent with theoretical expectations for our model of disk flaring. We point out that the value of the "Kennicutt-Schmidt" index depends on the thickness of the disk. With our simple feedback prescription (a single combined star formation event per cloud), we find that global spiral patterns are not sustained; less correlated feedback and smaller scale turbulence appear to be necessary for spiral patterns to persist.

Using one-dimensional hydrodynamic simulations including interstellar heating, cooling, and thermal conduction, we investigate nonlinear evolution of gas flow across galactic spiral arms. We model the gas as a non-self-gravitating, unmagnetized fluid, and follow its interaction with a stellar spiral potential in a local frame comoving with the stellar pattern. Initially uniform gas rapidly separates into warm and cold phases as a result of thermal instability (TI), and also forms a quasi-steady shock that prompts phase transitions. After saturation, the flow follows a recurring cycle: warm and cold phases in the interarm region are shocked and immediately cool to become a denser cold medium in the arm; post-shock expansion reduces the mean density to the unstable regime in the transition zone and TI subsequently mediates evolution back into warm and cold interarm phases. For our standard model with n_0 = 2 cm^-3, the gas resides in the dense arm, thermally-unstable transition zone, and interarm region for 14%, 22%, 64% of the arm-to-arm crossing time. These regions occupy 1%, 16%, and 83% of the arm-to-arm distance, respectively. Gas at intermediate temperatures represents ~25-30% of the total mass, similar to the fractions estimated from HI observations. Despite transient features and multiphase structure, the time-averaged shock profiles can be matched to that of a diffusive isothermal medium with temperature 1,000 K and "particle" mean free path of l_0 = 100 pc. Finally, we quantify numerical conductivity associated with translational motion of phase-separated gas on the grid, and show that convergence of numerical results requires the numerical conductivity to be comparable to or smaller than the physical conductivity. (Abridged)

Velocity anisotropy induced by MHD turbulence is investigated using computational simulations and molecular line observations of the Taurus molecular cloud. A new analysis method is presented to evaluate the degree and angle of velocity anisotropy using spectroscopic imaging data of interstellar clouds. The efficacy of this method is demonstrated on model observations derived from three dimensional velocity and density fields from the set of numerical MHD simulations that span a range of magnetic field strengths. The analysis is applied to 12CO J=1-0 imaging of a sub-field within the Taurus molecular cloud. Velocity anisotropy is identified that is aligned within 10 degrees of the mean local magnetic field direction derived from optical polarization measurements. Estimated values of the field strength based on velocity anisotropy are consistent with results from other methods. When combined with new column density measurements for Taurus, our magnetic field strength estimate indicates that the envelope of the cloud is magnetically subcritical. These observations favor strong MHD turbulence within the low density, sub-critical, molecular gas substrate of the Taurus cloud.

We review current understanding of star formation, outlining an overall theoretical framework and the observations that motivate it. A conception of star formation has emerged in which turbulence plays a dual role, both creating overdensities to initiate gravitational contraction or collapse, and countering the effects of gravity in these overdense regions. The key dynamical processes involved in star formation -- turbulence, magnetic fields, and self-gravity -- are highly nonlinear and multidimensional. Physical arguments are used to identify and explain the features and scalings involved in star formation, and results from numerical simulations are used to quantify these effects. We divide star formation into large-scale and small-scale regimes and review each in turn. Large scales range from galaxies to giant molecular clouds (GMCs) and their substructures. Important problems include how GMCs form and evolve, what determines the star formation rate (SFR), and what determines the initial mass function (IMF). Small scales range from dense cores to the protostellar systems they beget. We discuss formation of both low- and high-mass stars, including ongoing accretion. The development of winds and outflows is increasingly well understood, as are the mechanisms governing angular momentum transport in disks. Although outstanding questions remain, the framework is now in place to build a comprehensive theory of star formation that will be tested by the next generation of telescopes.

We use CO and H alpha velocity fields to study the gas kinematics in the spiral arms and interarms of M51 (NGC 5194), and fit the 2D velocity field to estimate the radial and tangential velocity components as a function of spiral phase (arm distance). We find large radial and tangential streaming velocities, which are qualitatively consistent with the predictions of density wave theory and support the existence of shocks. The streaming motions are complex, varying significantly across the galaxy as well as along and between arms. Aberrations in the velocity field indicate that the disk is not coplanar, perhaps as far in as 20\arcsec\ (800 pc) from the center. Velocity profile fits from CO and H alpha are typically similar, suggesting that most of the H alpha emission originates from regions of recent star formation. We also explore vortensity and mass conservation conditions. Vortensity conservation, which does not require a steady state, is empirically verified. The velocity and density profiles show large and varying mass fluxes, which are inconsistent with a steady flow for a single dominant global spiral mode. We thus conclude that the spiral arms cannot be in a quasi-steady state in any rotating frame, and/or that out of plane motions may be significant.

We have performed time-dependent numerical simulations of the interstellar medium (ISM) which account for galactic shear and magnetic fields, vertical gravity, and a radiative cooling function for atomic gas. This allows us to study the magnetorotational instability (MRI) in cloudy, vertically-stratified disks. As in previous unstratified models, we find that thermal instability interacts with MRI-driven turbulence and galactic shear to produce a network of cold, dense, filamentary clouds embedded in a warm diffuse ambient medium. This structure strongly resembles the morphology of HI gas observed in the 21 cm line. There is significant thermally-unstable gas, but the density and temperature distributions retain the twin peaks of the classical two-phase ISM. We analyze the vertical distributions of density and various pressure terms, and address what supports the ISM vertically. Turbulent velocities of the cold gas increase as the cold mass fraction decreases, but are generally low ~1-3 km/s near the midplane; they increase to > 5 km/s at high z. Finally, we argue that in the outer parts of galactic disks, MRI is likely able to prevent the development of self-gravitating instabilities and hence suppress star formation, even if cold gas is present.

Galactic disks consist of both stars and gas. The gas is more dynamically responsive than the stars, and strongly nonlinear structures and velocities can develop in the ISM even while stellar surface density perturbations remain fractionally small. We use 2D numerical simulations to explore formation of bound clouds and turbulence generation in the gas of two-component galactic disks. We represent the stars with collisionless particles and follow their orbits using a PM method, and treat the gas as an isothermal, unmagnetized fluid. The two components interact through a combined gravity. Using stellar parameters typical of mid-disk conditions, we find that models with gaseous Toomre parameter Q_g < Q_c ~ 1.4 experience gravitational runaway and eventually form bound condensations. This Q_c value is nearly the same as previously found for razor-thin, gas-only models, indicating that the destabilizing effect of live stars is offsets the reduced self-gravity of thick disks. This result is also consistent with empirical studies showing that star formation is suppressed when Q_g > 1-2. The bound gaseous clouds that form have mass 6x10^7 Msun each; these represent superclouds that would subsequently fragment into GMCs. Self-gravity and sheared rotation also interact to drive turbulence in the gas when Q_g > Q_c. This turbulence is anisotropic, with more power in sheared than compressive motions. The gaseous velocity dispersion is ~ 0.6 times the thermal speed when Q_g ~ Q_c. This suggests that gravity is important in driving ISM turbulence in many spiral galaxies, since the low efficiency of star formation naturally leads to a state of marginal instability.

Spiral shocks are potentially a major source of turbulence in the interstellar medium. To address this problem quantitatively, we use numerical simulations to investigate gas flow across spiral arms in vertically stratified, self-gravitating, magnetized models of galactic disks. Our models are isothermal, quasi-axisymmetric, and local in the quasi-radial direction while global in the vertical direction. We find that a stellar spiral potential perturbation promptly induces a spiral shock in the gas flow. For vertically stratified gas disks, the shock front in the radial-vertical plane is in general curved, and never achieves a steady state. This behavior is in sharp contrast to spiral shocks in two-dimensional (thin) disks, which are generally stationary. The non-steady motions in our models include large-amplitude quasi-radial flapping of the shock front. This flapping feeds random gas motions on the scale of the vertical disk thickness, which then cascades to smaller scales. The induced gas velocity dispersion in quasi-steady state exceeds the sonic value for a range of shock strengths, suggesting that spiral shocks are indeed an important generator of turbulence in disk galaxies.

We present a survey of spiral arm extinction substructure referred to as feathers in 223 spiral galaxies using HST WFPC2 images. The sample includes all galaxies in the RC3 catalog with cz < 5000 km/s, B_T < 15, i < 60 degrees, and types Sa--Sd with well-exposed broadband WFPC2 images. The detection frequency of delineated, periodic feathers in this sample is 20% (45 of 223). This work is consistent with Lynds (1970), who concluded that feathers are common in prototypical Sc galaxies; we find that feathers are equally common in Sb galaxies. Sb--Sc galaxies without clear evidence for feathers either had poorer quality images, or flocculent or complex structure. We did not find clearly defined feathers in any Scd--Sd galaxy. The probability of detecting feathers was highest (83%) for spirals with well-defined primary dust lanes (PDLs; the lanes which line the inner edge of an arm); well-defined PDLs were only noted in Sab--Sc galaxies. Consistent with earlier work, we find that neighboring feathers tend to have similar shapes and pitch angles. OB associations are often found lining feathers, and many feathers transition to the stellar substructures known as spurs (Elmegreen 1980). We find that feathers are coincident with interarm filaments strikingly revealed in Spitzer 8 micron images. Comparison with CO 1-0 maps of NGC 0628 and NGC 5194 from BIMA SONG shows that feathers originate at the PDL coincident with gas surface density peaks. Contrary to the appearance at 8 microns, the CO maps show that gas surface density in feathers decreases rapidly with distance from the PDL. Also, we find that the spacing between feathers decreases with increasing gas surface density, consistent with formation via a gravitational instability.

We investigate the formation of substructure in spiral galaxies using global MHD simulations, including gas self-gravity. Our models extend previous local models by Kim and Ostriker (2002) by including the full effects of curvilinear coordinates, a realistic log-spiral perturbation, self-gravitational contribution from 5 radial wavelengths of the spiral shock, and variation of density and epicyclic frequency with radius. We show that with realistic Toomre Q values, self-gravity and galactic differential rotation produce filamentary gaseous structures with kpc-scale separations, regardless of the strength -- or even presence -- of a stellar spiral potential. However, the growth of sheared features distinctly associated with the spiral arms, described as spurs or feathers in optical and IR observations of many spiral galaxies, requires a sufficiently strong spiral potential in self gravitating models. Unlike independently-growing ''background'' filaments, the orientation of arm spurs depends on galactic location. Inside corotation, spurs emanate outward, on the convex side of the arm; outside corotation, spurs grow inward, on the concave side of the arm. Based on spacing, orientation, and the relation to arm clumps, it is possible to distinguish ''true spurs'' that originate as instabilities in the spiral arms from independently growing ''background'' filaments. Our models also suggest that magnetic fields are important in preserving grand design spiral structure when gas in the arms fragments via self-gravity into GMCs.

(Abridged) We investigate the growth of spiral-arm substructure in vertically stratified, self-gravitating, galactic gas disks, using local numerical MHD simulations. Our new models extend our previous 2D studies (Kim & Ostriker 2002), which showed that a magnetized spiral shock in a thin disk can undergo magneto-Jeans instability (MJI), resulting in interarm spur structures and massive fragments. Similar spur features have recently been seen in high-resolution observations of several galaxies. Here, we consider two sets of numerical models: 2D models that use a thick-disk gravitational kernel, and 3D runs with explicit vertical stratification. When disks are sufficiently magnetized and self-gravitating, the result in both sorts of models is the growth of spiral-arm substructure similar to that in our previous razor-thin models. Reduced self-gravity in thick disks increases the spur spacing to ~10 times the Jeans length at the arm peak. Bound clouds that form from spur fragmentation have masses ~(1-3)x10^7 Msun each, a factor ~3-8 times larger than in razor-thin models with the same gas surface density and stellar spiral arm strength. We find that unmagnetized or weakly magnetized 2D models are unstable to the wiggle instability (WI) previously identified by Wada & Koda (2004), and proposed as a potential spur- and clump-forming mechanism. However,our fully 3D models do not show this effect. Non-steady motions and strong vertical shear prevent coherent vortical structures from forming, evidently suppressing the WI that appears in 2D (isothermal) runs. We also find no clear traces of Parker instability in the nonlinear spiral arm substructures that emerge (in self-gravitating models), although conceivably Parker modes may help seed the MJI at early stages since azimuthal wavelengths are similar.

We study the stability of proto-planetary disks with vertical velocity gradients in their equilibrium rotation rates; such gradients are expected to develop when dust settles into the midplane. Using a linear stability analysis of a simple three-layer model, we show that the onset of instability occurs at a larger value of the Richardson number, and therefore for a thicker layer, when the effects of Coriolis forces are included. This analysis also shows that even-symmetry (midplane-crossing) modes develop faster than odd-symmetry ones. These conclusions are corroborated by a large number of nonlinear numerical simulations with two different parameterized prescriptions for the initial (continuous) dust distributions. Based on these numerical experiments, the Richardson number required for marginal stability is more than an order of magnitude larger than the traditional 1/4 value. The dominant modes that grow have horizontal wavelengths of several initial dust scale heights, and in nonlinear stages mix solids fairly homogeneously over a comparable vertical range. We conclude that gravitational instability may be more difficult to achieve than previously thought, and that the vertical distribution of matter within the dust layer is likely globally, rather than locally, determined.

This paper reports on 3D numerical simulations of dynamics and thermodynamics in the diffuse ISM. Our models are local, account for sheared galactic rotation, magnetic fields, and realistic cooling, and resolve scales ~1-200 pc. This combination permits the study of quasi-steady turbulence in a cloudy medium representing the warm/cold atomic ISM. Turbulence is driven by the magnetorotational instability (MRI); our models are the first to study the saturated state of MRI under strongly inhomogeneous conditions, with cloud/intercloud contrasts of about 100. For volume-averaged densities n=0.25 - 4 cm^3, the saturated-state velocity dispersion ranges from 8-1 km/s, scaling as n^-0.77. The MRI is therefore likely quite important in driving turbulence in low-density ISM regions, both at high latitudes in the inner Galaxy,and throughout the far outer Galaxy. The MRI may be key to suppressing star formation at large radii in spiral galaxies, where without MRI-driven turbulence, a thin, cold layer could form and become gravitationally instable. Our presentation also includes analysis of thermal structure (temperature and density PDFs are broadened and include extreme departures from equilibrium, but retain the classical two-phase bimodal character); results on the distribution of clump masses (the mass spectrum peaks at about 100 solar masses); comparisons of saturated-state MRI scalings with single-phase simulation results (we find B^2 is independent of n); and examples of synthetic HI line profile maps (showing that physical clumps are not easily distinguished in velocity components, and vice versa).

The structure and dynamics of diffuse gas in the Milky Way and other disk galaxies may be strongly influenced by thermal and magnetorotational instabilities (TI and MRI) on scales of about 1-100 pc. We initiate a study of these processes, using two-dimensional numerical hydrodynamic and magnetohydrodynamic (MHD) simulations with conditions appropriate for the atomic interstellar medium (ISM). We demonstrate, consistent with previous work, that nonlinear development of ``pure TI'' produces a network of filaments that condense into cold clouds at their intersections, yielding a distinct two-phase warm/cold medium within about 20 Myr. TI-driven turbulent motions of the clouds saturate at subsonic amplitudes for uniform initial P/k=2000 K cm^-3. MRI has previously been studied in near-uniform media; our simulations include both TI+MRI models, which begin from uniform-density conditions, and cloud+MRI models, which begin with a two-phase cloudy medium. Both the TI+MRI and cloud+MRI models show that MRI develops within a few galactic orbital times, just as for a uniform medium. The mean separation between clouds can affect which MRI mode dominates the evolution. Provided intercloud separations do not exceed half the MRI wavelength, we find the MRI growth rates are similar to those for the corresponding uniform medium. This opens the possibility, if low cloud volume filling factors increase MRI dissipation times compared to those in a uniform medium, that MRI-driven motions in the ISM could reach amplitudes comparable to observed HI turbulent linewidths.

Using local 3D MHD simulations, we investigate ways in which galactic turbulence associated with the magnetorotational instability (MRI) may influence the formation and properties of GMCs. Our disk models are vertically stratified and subject to uniform shear. Initial magnetic fields are weak and purely vertical. For simplicity, we adopt an isothermal equation of state. We find that MRI-driven turbulence develops rapidly, with the saturated-state Shakura & Sunyaev parameter alpha~(0.15-0.3) dominated by Maxwell stresses. Many of the dimensionless characteristics of the turbulence (e.g. the ratio of the Maxwell to Reynolds stresses) are similar to results from previous MRI studies of accretion disks, hence insensitive to the degree of vertical disk compression, shear rate, and the presence of self-gravity. The density-weighted velocity dispersions in non- or weakly self-gravitating disks are sigma_x ~ sigma_y ~ (0.4-0.6)c_s and sigma_z~(0.2-0.3)c_s, suggesting that MRI can contribute significantly to the observed level of galactic turbulence. The saturated-state magnetic field strength B ~ 2 \mu G is similar to typical galactic values. When self-gravity is strong enough, MRI-driven high-amplitude density perturbations are swing-amplified to form Jeans-mass (~10^7 Msun) bound clouds. Compared to previous unmagnetized or strongly-magnetized disk models, the threshold for nonlinear instability in the present models occurs for surface densities at least 50% lower, corresponding to the Toomre parameter Q_th~1.6. We present evidence that self-gravitating clouds like GMCs formed under conditions similar to our models can lose much of their original spin angular momenta by magnetic braking, preferentially via fields threading near-perpendicularly to their spin axes.

Observations reveal concentrations of molecular line emission on the sky, called ``clumps,'' in dense, star-forming molecular clouds. These clumps are believed to be the eventual sites of star formation. We study the three-dimensional analogs of clumps using a set of self-consistent, time-dependent numerical models of molecular clouds. The models follow the decay of initially supersonic turbulence in an isothermal, self-gravitating, magnetized fluid. We find the following. (1) Clumps are intrinsically triaxial. This explains the observed deficit of clumps with a projected axis ratio near unity, and the apparent prolateness of clumps. (2) Simulated clump axes are not strongly aligned with the mean magnetic field within clumps, nor with the large-scale mean fields. This is in agreement with observations. (3) The clump mass spectrum has a high-mass slope that is consistent with the Salpeter value. There is a low-mass break in the slope at $\sim 0.5 \msun$, although this may depend on model parameters including numerical resolution. (4) The typical specific spin angular momentum of clumps is $4 \times 10^{22} {\rm cm^2 s^{-1}}$. This is larger than the median specific angular momentum of binary stars. Scaling arguments suggest that higher resolution simulations may soon be able to resolve the scales at which the angular momentum of binary stars is determined.

We analyze the spectral properties of driven, supersonic compressible magnetohydrodynamic (MHD) turbulence obtained via high-resolution numerical experiments, for application to understanding the dynamics of giant molecular clouds. Via angle-averaged power spectra, we characterize the transfer of energy from the intermediate, driving scales down to smaller dissipative scales, and also present evidence for inverse cascades that achieve modal-equipartition levels on larger spatial scales. Investigating compressive versus shear modes separately, we evaluate their relative total power, and find that as the magnetic field strength decreases, (1) the shear fraction of the total kinetic power decreases, and (2) slopes of power-law fits over the inertial range steepen. To relate to previous work on incompressible MHD turbulence, we present qualitative and quantitative measures of the scale-dependent spectral anisotropy arising from the shear-Alfv\'{e}n cascade, and show how these vary with changing mean magnetic field strength. Finally, we propose a method for using anisotropy in velocity centroid maps as a diagnostic of the mean magnetic field strength in observed cloud cores.

We use 3D MHD simulations to investigate nonlinear development of the Parker, magneto-Jeans (MJI), and swing mechanisms in galactic disks. The model disks are local, isothermal, and begin from a vertically-stratified equilibrium. We first construct axisymmetric equilibria and examine their stability. Finite disk thickness reduces the critical Toomre Q parameter below unity; we find Q_c \~ 0.75, 0.72, and 0.57 for \beta=\infty, 10, and 1 cases, respectively. We then pursue fully 3D models. In non-self-gravitating cases, the peak density enhancement from the `pure' Parker instability is less a factor of two. The dominant growing modes have radial wavelengths comparable to the disk scale height H, much shorter than the azimuthal wavelength (~10-20 H). Shearing disks, being more favorable to midplane-symmetric modes, have somewhat different late-time magnetic field profiles from nonshearing disks, but otherwise saturated states are similar. Late-time velocity fluctuations at 10% of the sound speed persist, but no characteristic structural signatures of Parker modes remain in the new equilibria. In self-gravitating cases, the development of density structure is qualitatively similar to our previous results from thin-disk simulations. The Parker instability, although it may help seed structure or tip the balance under marginal conditions, appears to play a secondary role. In shearing disks with Q less than a threshold level ~ 1, swing amplification can produce bound clouds of a few times the local Jeans mass. The most powerful cloud-condensing mechanism appears to be the MJI. Our simulations show that condensations of a local Jeans mass (~3\times 10^7 M_sun) grow very rapidly, supporting the idea that MJI is at least partly responsible for the formation of bound cloud complexes in spiral galaxies.

An important aspect of astrophysical MHD turbulence research is developing diagnostics to connect simulations with the observable universe. Turbulent systems are by definition structurally complex in all fluid variables (density, velocity, and magnetic field), such that they must be described statistically. By developing and applying diagnostic tools to simulation data, it is possible to interpret empirical laws for the statistical properties of observed systems in terms of fundamental dynamical processes, and to identify and calibrate robust probes of physical parameters that cannot be measured directly. Using several different examples, I describe how structural diagnostic analyses have already yielded significant insights into the nature of turbulent molecular clouds. I review results from several different groups, and discuss directions for future diagnostics to enhance our understanding of cloud structure and constrain models of the evolutionary course that governs star formation.

Intermediate-scale spurs are common in spiral galaxies, but perhaps most distinctively evident in a recent HST image of M51 (Scoville & Rector 2001). We investigate, using time-dependent numerical MHD simulations, how such spurs could form (and subsequently fragment) from the interaction of a gaseous ISM with a stellar spiral arm. We model the gaseous medium as a self-gravitating, magnetized, differentially-rotating, razor-thin disk. The basic flow shocks and compresses as it passes through a local segment of a tightly-wound, trailing stellar spiral arm, modeled as a rigidly-rotating gravitational potential. We first construct 1D profiles for flows with spiral shocks. When the post-shock Toomre parameter Q_sp is sufficiently small, self-gravity is too large for one-dimensional steady solutions to exist. The critical values of Q_sp are 0.8, 0.5, and 0.4 for our models with zero, sub-equipartition, and equipartition magnetic fields, respectively. We then study the growth of self-gravitating perturbations in fully 2D flows, and find that spur-like structures rapidly emerge in our magnetized models. We associate this gravitational instability with the magneto-Jeans mechanism, in which magnetic tension forces oppose the Coriolis forces. The shearing and expanding velocity field shapes the condensed material into spurs as it flows downstream from the arms. Although we find swing amplification can help form spurs when the arm-interarm contrast is moderate, unmagnetized systems that are quasi-axisymmetrically stable are generally also stable to nonaxisymmetric perturbations, suggesting that magnetic effects are essential. In nonlinear stages, the spurs in our models undergo fragmentation to form 4\times 10^6 solar mass clumps, which we suggest could evolve into bright arm/interarm HII regions as seen in spiral galaxies.

As increasingly precise information about the spectrum of the cosmic microwave background fluctuations is gathered with balloon and satellite experiments, interest has grown in foreground sources of opacity affecting these observations. One potentially important source is electron scattering produced by a post-recombination luminosity source, which would significantly attenuate the higher harmonics in the spectrum. If such an ionization source exists, then it would also heat the universe, hence increasing the Jeans mass and suppressing early gravitational structure formation. Here we consider the effects of such heating. We concentrate on one type of ionization source: luminosity generated by accretion onto primordial compact objects. We show that if such objects generate enough luminosity to affect the CMB power spectrum, then they would produce enough heat to prevent the formation of 1 sigma collapsed objects until redshifts of about 5, significantly less than the redshift at which baryonic collapse could otherwise occur. Such processes would leave signatures detectable by upcoming instruments such as NGST, SIRTF, and SWIFT.

We investigate the susceptibility of gaseous, magnetized galactic disks to formation of self-gravitating condensations using two-dimensional, local models. We focus on two issues: (1) determining the threshold condition for gravitational runaway, taking into account nonlinear effects, and (2) distinguishing the magneto-Jeans instability (MJI) that arises under inner-galaxy rotation curves from the modified swing amplification (MSA) that arises under outer-galaxy rotation curves. For axisymmetric density fluctuations, instability is known to require a Toomre parameter Q<1. For nonaxisymmetric fluctuations, any nonzero shear $q \equiv -d\ln \Omega /d \ln R$ winds up wavefronts such that in linear theory amplification saturates. Any Q threshold for nonaxisymmetric gravitational runaway must originate from nonlinear effects. We use numerical magnetohydrodynamic simulations to demonstrate the anticipated threshold phenomenon, to analyze the nonlinear processes involved, and to evaluate the critical value $Q_c$ for stabilization. We find $Q_c \sim 1.2-1.4$ for a wide variety of conditions, with the largest values corresponding to nonzero but subthermal mean magnetic fields. Our findings for $Q_c$ are similar to those inferred from thresholds for active star formation in the outer regions of spiral galaxies. MJI is distinct from MSA in that opposition to Coriolis forces by magnetic tension, rather than cooperation of epicyclic motion with kinematic shear, enables nonaxisymmetric density perturbations to grow. We suggest that under low-shear inner-disk conditions, $Q_c$ will be larger than that in outer disks by a factor $\sim (v_A/q c_s)^{1/2}$, where $v_A$ and $c_s$ are the respective Alfven and sound speeds.

We analyze the dynamics of the shell produced when a bow shock from a collimated jet propagates into the surrounding medium. Under interstellar conditions, the shock is radiative, and a ballistic approximation for the shell flow is appropriate beyond the working surface where the jet impacts its surroundings. The solution is then determined by the ambient and jet densities and velocities and by the momentum impulse applied in the working surface. Using estimates for these impulses (confirmed by separate numerical simulations), we obtain solutions for the shell structure, and for the range of velocities in the shell at any point. We provide predictions for the position-velocity and mass-velocity relations expected for plane-of-sky bow-shock shells, and for the bulk shell properties. In a companion paper, we show that these analytic solutions are in excellent agreement with the results of direct numerical simulations. We argue that classical molecular (CO) outflows cannot be purely jet-driven, because the bow-shock shell solutions are much too elongated compared with observations. Finally, we suggest that the ``spur'' structures seen in position-velocity diagrams of observed molecular outflows are the manifestation of internal bow shocks which may be fit with our simple dynamical models.

We present two-dimensional hydrodynamical simulations of both jet- and wind-driven models for protostellar outflows in order to make detailed comparisons to the kinematics of observed molecular outflows. Comparing the different simulations with observations, we find that some outflows, e.g., HH 212, show features consistent with the jet-driven model, while others, e.g., VLA 05487, are consistent with the wind-driven model.

We analyze the temporal evolution of accretion onto rotating black holes subject to large-scale magnetic torques. Wind torques alone drive a disk towards collapse in a finite time $\sim t_{ff} E_k/E_B$, where $t_{ff}$ is the initial free-fall time and $E_k/E_B$ is the ratio of kinetic-to-poloidal magnetic energy. Additional spin-up torques from a rapidly rotating black hole can arrest the disk's inflow. We associate short/long gamma-ray bursts with hyperaccretion/suspended-accretion onto slowly/rapidly spinning black holes. This model predicts afterglow emission from short bursts, and may be tested by HETE-II.

We use 3D numerical MHD simulations to follow the evolution of cold, turbulent, gaseous systems with parameters representing GMC conditions. We study three cloud simulations with varying mean magnetic fields, but identical initial velocity fields. We show that turbulent energy is reduced by a factor two after 0.4-0.8 flow crossing times (2-4 Myr), and that the magnetically supercritical cloud models collapse after ~6 Myr, while the subcritical cloud does not collapse. We compare density, velocity, and magnetic field structure in three sets of snapshots with matched Mach numbers. The volume and column densities are both log-normally distributed, with mean volume density a factor 3-6 times the unperturbed value, but mean column density only a factor 1.1-1.4 times the unperturbed value. We use a binning algorithm to investigate the dependence of kinetic quantities on spatial scale for regions of column density contrast (ROCs). The average velocity dispersion for the ROCs is only weakly correlated with scale, similar to the mean size-linewidth relation for clumps within GMCs. ROCs are often superpositions of spatially unconnected regions that cannot easily be separated using velocity information; the same difficulty may affect observed GMC clumps. We analyze magnetic field structure, and show that in the high density regime, total magnetic field strengths increase with density with logarithmic slope 1/3 -2/3. Mean line-of-sight magnetic field strengths vary widely across a projected cloud, and do not correlate with column density. We compute simulated interstellar polarization maps at varying orientations, and determine that the Chandrasekhar-Fermi formula multiplied by a factor ~0.5 yields a good estimate of the plane-of sky magnetic field strength provided the dispersion in polarization angles is < 25 degrees.

We investigate shear and buoyancy instabilities in radially stratified, magnetized, cylindrical flows, for application to magnetocentrifugally driven winds and to magnetized accretion disks. We identify and study nine principal types of instabilities/overstabilities. When magnetic fields are predominantly toroidal, as in protostellar winds, the system exhibits axisymmetric fundamental (FM) and toroidal resonance modes, axisymmetric and non-axisymmetric toroidal buoyancy (TB) modes, and non-axisymmetric magnetorotational instability (MRI) modes. Winds with sufficiently steep field gradients are unstable to the FM, which promotes narrow dense jets in the centers of wider winds. The TB modes promote small-scale radial mixing. The MRIs have very low growth rates under low-temperature wind conditions. The stabilization of buoyancy modes by shear, and of MRIs by compressibility, may be important in allowing protostellar winds to propagate over vast distances in space. When magnetic fields are predominantly poloidal, as in winds close to their source or in astrophysical disks, the system exhibits axisymmetric Balbus-Hawley (BH), poloidal buoyancy (PB), non-axisymmetric geometric buoyancy (GPB), and poloidal resonance modes. The BH mode has the fastest growth rate. The PB mode promotes radial mixing on small scales. The GPB mode at high m is readily stabilized by shear. We extend previous MRI analyses to focus on compressibility. We introduce a ``coherent wavelet'' technique to derive closed-form expressions for instantaneous instability criteria, growth rates, and net amplification factors for non-axisymmetric MRIs in compressible flows with both poloidal and toroidal fields. We confirm that these are in excellent agreement with the results of shearing-sheet temporal integrations.

We review recent results from numerical simulations and related models of MHD turbulence in the interstellar medium (ISM) and in molecular clouds. We discuss the implications of turbulence for the processes of cloud formation and evolution, and the determination of clouds' physical properties. Numerical simulations of the turbulent ISM to date have included magnetic fields, self-gravity, parameterized heating and cooling, modeled star formation and other turbulent inputs. The structures which form reproduce well observed velocity-size scaling properties, while predicting the non-existence of a general density-size scaling law. Criteria for the formation of gravitationally-bound structures by turbulent compression are summarized. For flows with equations of state $P\propto \rho^\gamma$, the statistics of the density field depend on the exponent $\gamma$. Numerical simulations of both forced and decaying MHD compressible turbulence have shown that the decay rate is comparable to the non-magnetic case. For virialized clouds, the turbulent decay time is shorter than the gravitational free-fall time, so wholesale cloud collapse is only prevented by ongoing turbulent inputs and/or a strong mean magnetic field. Finally, perspectives for future work in this field are briefly discussed.

The molecular component of the Galaxy is comprised of turbulent, magnetized clouds, many of which are self-gravitating and form stars. To understand how these clouds' evolution may depend on their level of turbulence, mean magnetization, and degree of self-gravity, we perform a survey of direct numerical MHD simulations in 2.5 dimensions. Two of our independent survey parameters allow us to model clouds which either meet or fail conditions for magneto-Jeans stability and magnetic criticality. Our third survey parameter allows us to initiate turbulence of either sub- or super-Alfv\'enic amplitude. We evaluate the times for each cloud model to become gravitationally bound, and measure each model's kinetic energy loss over the fluid flow crossing time. We compare the evolution of density and magnetic field structural morphology, and quantify the differences in the density contrast generated by internal stresses, for models of differing mean magnetization. We find that the values of (c_s/v_A)^2 and L/L_Jeans, but not the initial Mach number v/c_s, determine the time for cloud gravitational binding and collapse. We find, contrary to some previous expectations, less than a factor of two difference between turbulent decay times for models with varying magnetic field strength. In all models, we find turbulent amplification in the magnetic field strength, with the turbulent magnetic energy between 25-60% of the turbulent kinetic energy after one flow crossing time. We find that for non-self-gravitating stages of evolution, the mass-averaged density contrast magnitudes <log(\rho/\bar\rho)> are in the range 0.2-0.5; we note that only the more strongly-magnetized models appear consistent with estimates of clump/interclump density contrasts inferred in Galactic GMCs.

Using time-dependent linear perturbation theory, we evaluate the dynamical friction force on a massive perturber M_p traveling at velocity V through a uniform gaseous medium of density rho_0 and sound speed c_s. This drag force acts in the direction -\hat V, and arises from the gravitational attraction between the perturber and its wake in the ambient medium. For supersonic motion (M=V/c_s>1), the enhanced-density wake is confined to the Mach cone trailing the perturber; for subsonic motion (M<1), the wake is confined to a sphere of radius c_s t centered a distance V t behind the perturber. Inside the wake, surfaces of constant density are hyperboloids or oblate spheroids for supersonic or subsonic perturbers, respectively, with the density maximal nearest the perturber. The dynamical drag force has the form F_df= - I 4\pi (G M_p)^2\rho_0/V^2. We evaluate I analytically; its limits are I\to M^3/3 for M<<1 and I\to ln(Vt/r_min) for M>>1. We compare our results to the Chandrasekhar formula for dynamical friction in a collisionless medium, noting that the gaseous drag is generally more efficient when M>1 but less efficient when M<1. To allow simple estimates of orbit evolution in a gaseous protogalaxy or proto-star cluster, we use our formulae to evaluate the decay times of a (supersonic) perturber on a near-circular orbit in an isothermal \rho\propto r^{-2} halo, and of a (subsonic) perturber on a near-circular orbit in a constant-density core. We also mention the relevance of our calculations to protoplanet migration in a circumstellar nebula.

We report results of a three dimensional, high resolution (up to 512^3) numerical investigation of supersonic compressible magnetohydrodynamic turbulence. We consider both forced and decaying turbulence. The model parameters are appropriate to conditions found in Galactic molecular clouds. We find that the dissipation time of turbulence is of order the flow crossing time or smaller, even in the presence of strong magnetic fields. About half the dissipation occurs in shocks. Weak magnetic fields are amplified and tangled by the turbulence, while strong fields remain well ordered.

This paper reviews current theoretical work on the various stages of accretion in protostars, and the relationship of these ideal stages to the spectral classes of observed protostellar systems. I discuss scaling relationships that have been obtained for models of pre-stellar cores as they evolve by ambipolar diffusion toward a central singularity, and expectations for the dynamical evolution as the cores collapse radially to form a rotating disk. I summarize work that suggests accretion in T Tauri systems may be limited by ionization to disk surface layers, and implications for the variation of the disk accretion rate with radius. Finally, I describe models for the asymptotic structure of free magnetocentrifugal winds, and show that the constant-density surfaces in these winds may be strongly collimated even when the streamlines are not.

We construct a two-parameter family of models for self-collimated, radially self-similar magnetized outflows from accretion disks. A flow at zero initial poloidal speed leaves the surface of a rotating disk and is accelerated and redirected toward the pole by helical magnetic fields threading the disk. At large distances from the disk, the flow streamlines asymptote to wrap around the surfaces of nested cylinders. In constrast to previous disk wind modeling, we have explicitly implemented the cylindrical asymptotic boundary condition to examine the consequences for flow dynamics. The solutions are characterized by the logarithmic gradient of the magnetic field strength and the ratios between the footpoint radius R_0 and asymptotic radius R_1 of streamlines; the Alfven radius must be found as an eigenvalue. Cylindrical solutions require the magnetic field to drop less steeply than 1/R. We find that the asymptotic poloidal speed on any streamline is typically just a few tenths of the Kepler speed at the corresponding disk footpoint. The asymptotic toroidal Alfven speed is, however, a few times the footpoint Kepler speed. We discuss the implications of the models for interpretations of observed optical jets and molecular outflows from young stellar systems. We suggest that the difficulty of achieving strong collimation in vector velocity simultaneously with a final speed comparable to the disk rotation rate argues against isolated jets and in favor of models with broader winds.

Using self-consistent magnetohydrodynamic (MHD) simulations, we explore the hypothesis that nonlinear MHD waves dominate the internal dynamics of galactic molecular clouds. We employ an isothermal equation of state and allow for self-gravity. We adopt ``slab-symmetry,'' which permits motions $\bf v_\perp$ and fields $\bf B_\perp$ perpendicular to the mean field, but permits gradients only parallel to the mean field. The Alfv\'en speed $v_A$ exceeds the sound speed $c_s$ by a factor $3-30$. We simulate the free decay of a spectrum of Alfv\'en waves, with and without self-gravity. We also perform simulations with and without self-gravity that include small-scale stochastic forcing. Our major results are as follows: (1) We confirm that fluctuating transverse fields inhibit the mean-field collapse of clouds when the energy in Alfv\'en- like disturbances remains comparable to the cloud's gravitational binding energy. (2) We characterize the turbulent energy spectrum and density structure in magnetically-dominated clouds. The spectra evolve to approximately $v_{\perp,\,k}^2\approx B_{\perp,\,k}^2/4\pi\rho\propto k^{-s}$ with $s\sim 2$, i.e. approximately consistent with a ``linewidth-size'' relation $\sigma_v(R) \propto R^{1/2}$. The simulations show large density contrasts, with high density regions confined in part by the fluctuating magnetic fields. (3) We evaluate the input power required to offset dissipation through shocks, as a function of $c_s/v_A$, the velocity dispersion $\sigma_v$, and the scale $\lambda$ of the forcing. In equilibrium, the volume dissipation rate is $5.5(c_s/v_a)^{1/2} (\lambda/L)^{-1/2}\times \rho \sigma_v^3/L$, for a cloud of linear size $L$ and density $\rho$. (4) Somewhat speculatively, we apply our results to a ``typical'' molecular cloud. The mechanical power input required

We describe hot, optically-thin solutions for one-temperature accretion disks around black holes. We include cooling by synchrotron, bremsstrahlung, and Comptonization. Our solutions are thermally and viscously stable, with gas temperatures on the order of $T \sim 10^9-10^{10.7}$K. The thermal stability is a direct result of the inclusion of synchrotron cooling. The new solution branch is related to the advection-dominated solution for a two-temperature gas described by Narayan \& Yi (1995b). It is present only for mass accretion rates less than some critical $\dot{M}_{crit}$ which depends on the radius $R$ and viscosity parameter $\alpha$. The solutions are advection-dominated for extremely low values of $\dot{M}$. However, for a range of intermediate accretion rates, the new solutions are both hot ($T \sim 10^{10}$K) and cooling-dominated. Because of this new feature, one-temperature solutions are significantly more luminous than the corresponding two temperature solutions. The radial profile of the new solutions is unusual. The inner parts of the flow are cooling-dominated and have a disk-like geometry, while the outer parts are fully advection-dominated and nearly quasi-spherical.