Source author record

Eve C. Ostriker

Eve C. Ostriker appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.CO astro-ph.SR astro-ph astro-ph.HE astro-ph.IM physics.plasm-ph

Catalog footprint

What is connected

51works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Cosmic-Ray Transport in Varying Galactic Environments

We study the propagation of mildly-relativistic cosmic rays (CRs) in multiphase interstellar medium environments with conditions typical of nearby disk galaxies. We employ the techniques developed in Armillotta+21 to post-process three high-resolution TIGRESS magnetohydrodynamic simulations modeling local patches of star-forming galactic disks. Together, the three simulations cover a wide range of gas surface density, gravitational potential, and star formation rate (SFR). Our prescription for CR propagation includes the effects of advection by the background gas, streaming along the magnetic field at the local ion Alfvén speed, and diffusion relative to the Alfvén waves, with the diffusion coefficient set by the balance between streaming-driven Alfvén wave excitation and damping mediated by local gas properties. We find that the combined transport processes are more effective in environments with higher SFR. These environments are characterized by higher-velocity hot outflows (created by clustered supernovae) that rapidly advect CRs away from the galactic plane. As a consequence, the ratio of midplane CR pressure to midplane gas pressures decreases with increasing SFR. We also use the post-processed simulations to make predictions regarding potential dynamical impacts of CRs. The relatively flat CR pressure profiles near the midplane argue that they would not provide significant support against gravity for most of the ISM mass. However, the CR pressure gradients are larger than the other pressure gradients in the extra-planar region (|z|>0.5 kpc), suggesting that CRs may affect the dynamics of galactic fountains and/or winds. The degree of this impact is expected to increase in environments with lower SFR.

preprint2022arXiv

Slow Star Formation in the Milky Way: Theory Meets Observations

The observed star formation rate of the Milky Way can be explained by applying a metallicity-dependent factor to convert CO luminosity to molecular gas mass and a star formation efficiency per free-fall time that depends on the virial parameter of a molecular cloud. These procedures also predict the trend of star formation rate surface density with Galactocentric radius. The efficiency per free-fall time variation with virial parameter plays the major role in bringing theory into agreement with observations for the total star formation rate, while the metallicity dependence of the CO luminosity to mass conversion is most notable in the variation with Galactocentric radius. Application of these changes resolves a factor of over 100 discrepancy between observed and theoretical star formation rates that has been known for nearly 50 years.

preprint2022arXiv

The Life and Times of Giant Molecular Clouds

Giant molecular clouds (GMCs) are the sites of star formation and stellar feedback in galaxies. Their properties set the initial conditions for star formation and their lifecycles determine how feedback regulates galaxy evolution. In recent years, the advent of high-resolution telescopes has enabled systematic GMC-scale studies of the molecular interstellar medium in nearby galaxies, now covering a wide range of physical conditions and allowing for the first studies of how GMC properties depend on galactic environment. These observational developments have been accompanied by numerical simulations of improving resolution that are increasingly accurately accounting for the effects of the galactic-scale environment on GMCs, while simultaneously improving the treatment of the small-scale processes of star-formation and stellar feedback within them. The combination of these recent developments has greatly improved our understanding of the formation, evolution, and destruction of GMCs. We review the current state of the field, highlight current open questions, and discuss promising avenues for future studies.

preprint2021arXiv

Cosmic Ray Transport in the Ionized and Neutral ISM: MHD-PIC Simulations and Effective Fluid Treatments

Cosmic rays (CRs) have critical impacts in the multiphase interstellar medium (ISM), driving dynamical motions in low-density plasma and modifying the ionization state, temperature, and chemical composition of higher-density atomic and molecular gas. We present a study of CR propagation between the ionized ISM and a neutral cloud. Using one-dimensional magnetohydrodynamic particle-in-cell simulations which include ion-neutral drag to damp Alfv$\acute{\text{e}}$n waves in the cloud, we self-consistently evolve the kinetic physics of CRs and fluid dynamics of the multiphase gas. By introducing the cloud in our periodic domain, our simulations break translational symmetry and allow the emergence of spatial structure in the CR distribution function. A negative spatial gradient forms across the fully-ionized ISM region while a positive gradient forms across the neutral cloud. We connect our results with CR hydrodynamics formulations by computing the wave-particle scattering rates as predicted by quasilinear, fluid, and Fokker-Planck theory. For momenta where the mean free path is short relative to the box size, we find excellent agreement among all scattering rates. By exploring different cloud sizes and ion-neutral collision rates, we show that our results are robust. Our work provides a first-principles verification of CR hydrodynamics when particles stream down their pressure gradient, and opens a pathway toward comprehensive calibrations of transport coefficients from self-generated Alfv$\acute{\text{e}}$n wave scattering with CRs.

preprint2021arXiv

Effects of Varying Mass Inflows on Star Formation in Nuclear Rings of Barred Galaxies

Observations indicate that the star formation rate (SFR) of nuclear rings varies considerably with time and is sometimes asymmetric rather than being uniform across a ring. To understand what controls temporal and spatial distributions of ring star formation, we run semi-global, hydrodynamic simulations of nuclear rings subject to time-varying and/or asymmetric mass inflow rates. These controlled variations in the inflow lead to variations in the star formation, while the ring orbital period ($18\,{\rm Myr}$) and radius ($600\,{\rm pc}$) remain approximately constant. We find that both the mass inflow rate and supernova feedback affect the ring SFR. An oscillating inflow rate with period $Δτ_\text{in}$ and amplitude 20 causes large-amplitude (a factor of $\gtrsim 5$), quasi-periodic variations of the SFR, when $Δτ_\text{in} \gtrsim 50\,{\rm Myr}$. We find that the time-varying ISM weight and midplane pressure track each other closely, establishing an instantaneous vertical equilibrium. The measured time-varying depletion time is consistent with the prediction from self-regulation theory provided the time delay between star formation and supernova feedback is taken into account. The supernova feedback is responsible only for small-amplitude (a factor of $\sim 2$) fluctuations of the SFR with a timescale $\lesssim 40\,{\rm Myr}$. Asymmetry in the inflow rate does not necessarily lead to asymmetric star formation in nuclear rings. Only when the inflow rate from one dust lane is suddenly increased by a large factor, the rings undergo a transient period of lopsided star formation.

preprint2020arXiv

Diffuse Ionized Gas in Simulations of Multiphase, Star-Forming Galactic Disks

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$α$ 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$α$ surface brightness profiles in edge-on galaxies.

preprint2020arXiv

Dynamical Equilibrium in the Molecular ISM in 28 Nearby Star-Forming 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, $Σ_\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 $Σ_\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.

preprint2020arXiv

First results from SMAUG: Characterization of Multiphase Galactic Outflows from a Suite of Local Star-Forming Galactic Disk Simulations

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 (Σ_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 Σ_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 Σ_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.

preprint2020arXiv

First results from SMAUG: Insights into star formation conditions from spatially-resolved ISM properties in TNG50

Physical and chemical 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 to explore the physical parameter space of 8 resolved ISM properties in star-forming regions to constrain the areas of this hyperspace over which most star-forming environments exist. We deconstruct our simulated galaxies spanning a wide range of mass (M$_\star = 10^{7-11}$ M$_\odot$) and redshift ($0 \leq z \leq 3$) into kpc-sized regions, and statistically analyze the gas/stellar surface densities, gas metallicity, vertical stellar velocity dispersion, epicyclic frequency and dark-matter volumetric density representative of each region in the context of their star formation activity and galactic environment (radial galactocentric location). By examining the star formation rate (SFR) weighted distributions of these properties, we show that stars primarily form in two spatially distinct environmental regimes, which are brought about by an underlying bi-component radial SFR surface density profile in galaxies. We examine how the relative prominence of these two regimes depends on host galaxy mass and cosmic time. We also compare our findings with those from integral field spectroscopy observations and achieve a good overall agreement. Further, using dimensionality reduction, we characterise the aforementioned hyperspace to reveal a high-degree of multicollinearity in relationships amongst ISM properties that drive the distribution of star formation at kpc-scales. Based on this, we show that a reduced 3D representation underpinned by a multi-variate radius relationship is sufficient to capture most of the variance in the original 8D space.

preprint2020arXiv

Local Simulations of Spiral Galaxies with the TIGRESS Framework: I. Star Formation and Arm Spurs/Feathers

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.

preprint2020arXiv

Multiphase Gas and the Fractal Nature of Radiative Turbulent Mixing Layers

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 $ρ_{\rm cold}/ρ_{\rm hot}$ of the system. Entrainment is expected to be enhanced in environments with short $t_{\rm cool}$, large $v_{\rm rel}$, and large $ρ_{\rm cold}/ρ_{\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.

preprint2020arXiv

Self-gravitating Filament Formation from Shocked Flows: Velocity Gradients across Filaments

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 {Δv_h}^2/(GM/L)$, where $Δ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).

preprint2020arXiv

The Impact of Type Ia Supernovae in Quiescent Galaxies: I. Formation of the Multiphase Interstellar medium

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.

preprint2020arXiv

The Impact of Type Ia Supernovae in Quiescent Galaxies: II. Energetics and Turbulence

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.

preprint2020arXiv

The Physical Nature of Starburst-Driven Galactic Outflows

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.

preprint2019arXiv

Cloud Properties and Correlations with Star Formation in Numerical Simulations of the Three-Phase ISM

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.

preprint2019arXiv

Kinematics and Dynamics of Multiphase Outflows in Simulations of the Star-Forming Galactic ISM

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.

preprint2016arXiv

Disruption of Molecular Clouds by Expansion of Dusty H II Regions

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 $ε_{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 $ε_{min} \lesssim 10$%, with comparable gas and radiation pressure effects on shell expansion. Massive cluster-forming clumps require a significantly higher efficiency of $ε_{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.

preprint2016arXiv

Galaxy Outflows Without Supernovae

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.

preprint2016arXiv

Superbubbles in the Multiphase ISM and the Loading of Galactic Winds

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.

preprint2015arXiv

Anisotropic Formation of Magnetized Cores in Turbulent Clouds

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.

preprint2015arXiv

Momentum Injection by Supernovae in the Interstellar Medium

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).

preprint2015arXiv

Numerical Simulations of Turbulent Molecular Clouds Regulated by Reprocessed Radiation Feedback from Nascent Super Star Clusters

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 $κ$, the opacity of the gas to infrared (IR) radiation. We find that the single most important parameter determining the evolutionary outcome is $κ$, with $κ\gtrsim 15 \text{ cm}^2 \text{ g}^{-1}$ needed to disrupt clouds. For $κ= 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.

preprint2015arXiv

Prestellar Core Formation, Evolution, and Accretion from Gravitational Fragmentation in Turbulent Converging Flows

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.

preprint2015arXiv

Vertical Equilibrium, Energetics, and Star Formation Rates in Magnetized Galactic Disks Regulated by Momentum Feedback from Supernovae

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:Π_mag,t:Π_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" η_c=P_c/Σ_SFR (in suitable units) for each pressure component, finding that η_turb~4 and η_th~1 are the same for MHD as in previous HD simulations, and η_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 $η$ from turbulent magnetic fields reduces the predicted Σ_SFR by ~25% relative to the HD case.

preprint2014arXiv

ALMA Reveals the Molecular Medium Fueling the Nearest Nuclear Starburst

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.

preprint2014arXiv

CARMA Large Area Star Formation Survey: Structure and Kinematics of Dense Gas in Serpens Main

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.

preprint2014arXiv

Formation of Magnetized Prestellar Cores with Ambipolar Diffusion and Turbulence

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.

preprint2014arXiv

Three Dimensional Hydrodynamic Simulations of Multiphase Galactic Disks with Star Formation Feedback: II. Synthetic HI 21 cm Line Observations

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^{-τ(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})τ(v_{ch})/[[1-e^{-τ(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 $τ(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 $τ>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.

preprint2013arXiv

A Two-moment Radiation Hydrodynamics Module in Athena Using a Time-explicit Godunov Method

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.

preprint2013arXiv

Formation of Molecular Clouds and Global Conditions for Star Formation

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.

preprint2013arXiv

The Starburst-Driven Molecular Wind in NGC 253 and the Suppression of Star Formation

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α 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 η~1-3, establishing the importance of the starburst-driven wind in limiting the star formation activity and the final stellar content.

preprint2013arXiv

Three Dimensional Hydrodynamic Simulations of Multiphase Galactic Disks with Star Formation Feedback: I. Regulation of Star Formation Rates

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 $Σ_{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 $Σ_{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 $Σ_{SFR}$ apply. For outer disk regions, this results in $Σ_{SFR} \propto Σ\sqrt{ρ_{sd}}$, where $Σ$ is the total gas surface density and $ρ_{sd}$ is the midplane density of the stellar disk (plus dark matter). This scaling law arises because $ρ_{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.

preprint2012arXiv

A General Model for the CO-H2 Conversion Factor in Galaxies with Applications to the Star Formation Law

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.

preprint2012arXiv

Implementation of Sink Particles in the Athena Code

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.

preprint2012arXiv

Maximally Star-Forming Galactic Disks II. Vertically-Resolved Hydrodynamic Simulations of Starburst Regulation

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 $Σ$, 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 $Σ_{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 $Σ\in$(100--1000) \msunpc, we find $Σ_{SFR}\in$(0.1--4) \sfrunits, generally following a $Σ_{SFR}\propto Σ^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 $Σ_{SFR}-Σ$ relation in high density environments, provide strong evidence for the self-regulation of star formation.

preprint2012arXiv

Numerical Modeling of Multiphase, Turbulent Galactic Disks with Star Formation Feedback

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.

preprint2011arXiv

Ambipolar Diffusion in Action: Transient C shock Structure and Prestellar Core Formation

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.

preprint2011arXiv

Dense core formation in supersonic turbulent converging flows

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.

preprint2011arXiv

Maximally Star-Forming Galactic Disks I. Starburst Regulation Via Feedback-Driven Turbulence

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.

preprint2011arXiv

Modeling CO Emission: II. The Physical Characteristics that Determine the X factor in Galactic Molecular Clouds

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.

preprint2011arXiv

Regulation of Star Formation Rates in Multiphase Galactic Disks: Numerical Tests of the Thermal/Dynamical Equilibrium Model

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.

preprint2011arXiv

The CO-H2 Conversion Factor in Disc Galaxies and Mergers

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.

preprint2010arXiv

Galactic Spiral Shocks with Thermal Instability in Vertically Stratified Galactic Disks

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 $Σ_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).

preprint2010arXiv

Regulation of Star Formation Rates in Multiphase Galactic Disks: a Thermal/Dynamical Equilibrium Model

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.

preprint2010arXiv

Star Formation and Gas Dynamics in Galactic Disks: Physical Processes and Numerical Models

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.

preprint2007arXiv

Gravitational Runaway and Turbulence Driving in Star-Gas Galactic Disks

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.

preprint2007arXiv

Models of Vertically Stratified Two-Phase ISM Disks with MRI-Driven Turbulence

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.

preprint2006arXiv

Global Modeling of Spur Formation in Spiral Galaxies

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.

preprint2006arXiv

Interstellar Turbulence Driving by Galactic Spiral Shocks

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.

preprint1996arXiv

Can Nonlinear Hydromagnetic Waves Support a Self-Gravitating Cloud?

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én speed $v_A$ exceeds the sound speed $c_s$ by a factor $3-30$. We simulate the free decay of a spectrum of Alfvén 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én- 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πρ\propto k^{-s}$ with $s\sim 2$, i.e. approximately consistent with a ``linewidth-size'' relation $σ_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 $σ_v$, and the scale $λ$ of the forcing. In equilibrium, the volume dissipation rate is $5.5(c_s/v_a)^{1/2} (λ/L)^{-1/2}\times ρσ_v^3/L$, for a cloud of linear size $L$ and density $ρ$. (4) Somewhat speculatively, we apply our results to a ``typical'' molecular cloud. The mechanical power input required

Eve C. Ostriker

What is connected

Connect this record

See the researcher in context

Building this map preview

51 published item(s)

Cosmic-Ray Transport in Varying Galactic Environments

Slow Star Formation in the Milky Way: Theory Meets Observations

The Life and Times of Giant Molecular Clouds

Cosmic Ray Transport in the Ionized and Neutral ISM: MHD-PIC Simulations and Effective Fluid Treatments

Effects of Varying Mass Inflows on Star Formation in Nuclear Rings of Barred Galaxies

Diffuse Ionized Gas in Simulations of Multiphase, Star-Forming Galactic Disks

Dynamical Equilibrium in the Molecular ISM in 28 Nearby Star-Forming Galaxies

First results from SMAUG: Characterization of Multiphase Galactic Outflows from a Suite of Local Star-Forming Galactic Disk Simulations

First results from SMAUG: Insights into star formation conditions from spatially-resolved ISM properties in TNG50

Local Simulations of Spiral Galaxies with the TIGRESS Framework: I. Star Formation and Arm Spurs/Feathers

Multiphase Gas and the Fractal Nature of Radiative Turbulent Mixing Layers

Self-gravitating Filament Formation from Shocked Flows: Velocity Gradients across Filaments

The Impact of Type Ia Supernovae in Quiescent Galaxies: I. Formation of the Multiphase Interstellar medium

The Impact of Type Ia Supernovae in Quiescent Galaxies: II. Energetics and Turbulence

The Physical Nature of Starburst-Driven Galactic Outflows

Cloud Properties and Correlations with Star Formation in Numerical Simulations of the Three-Phase ISM

Kinematics and Dynamics of Multiphase Outflows in Simulations of the Star-Forming Galactic ISM

Disruption of Molecular Clouds by Expansion of Dusty H II Regions

Galaxy Outflows Without Supernovae

Superbubbles in the Multiphase ISM and the Loading of Galactic Winds

Anisotropic Formation of Magnetized Cores in Turbulent Clouds

Momentum Injection by Supernovae in the Interstellar Medium

Numerical Simulations of Turbulent Molecular Clouds Regulated by Reprocessed Radiation Feedback from Nascent Super Star Clusters

Prestellar Core Formation, Evolution, and Accretion from Gravitational Fragmentation in Turbulent Converging Flows

Vertical Equilibrium, Energetics, and Star Formation Rates in Magnetized Galactic Disks Regulated by Momentum Feedback from Supernovae

ALMA Reveals the Molecular Medium Fueling the Nearest Nuclear Starburst

CARMA Large Area Star Formation Survey: Structure and Kinematics of Dense Gas in Serpens Main

Formation of Magnetized Prestellar Cores with Ambipolar Diffusion and Turbulence

Three Dimensional Hydrodynamic Simulations of Multiphase Galactic Disks with Star Formation Feedback: II. Synthetic HI 21 cm Line Observations

A Two-moment Radiation Hydrodynamics Module in Athena Using a Time-explicit Godunov Method

Formation of Molecular Clouds and Global Conditions for Star Formation

The Starburst-Driven Molecular Wind in NGC 253 and the Suppression of Star Formation

Three Dimensional Hydrodynamic Simulations of Multiphase Galactic Disks with Star Formation Feedback: I. Regulation of Star Formation Rates

A General Model for the CO-H2 Conversion Factor in Galaxies with Applications to the Star Formation Law

Implementation of Sink Particles in the Athena Code

Maximally Star-Forming Galactic Disks II. Vertically-Resolved Hydrodynamic Simulations of Starburst Regulation

Numerical Modeling of Multiphase, Turbulent Galactic Disks with Star Formation Feedback

Ambipolar Diffusion in Action: Transient C shock Structure and Prestellar Core Formation

Dense core formation in supersonic turbulent converging flows

Maximally Star-Forming Galactic Disks I. Starburst Regulation Via Feedback-Driven Turbulence

Modeling CO Emission: II. The Physical Characteristics that Determine the X factor in Galactic Molecular Clouds

Regulation of Star Formation Rates in Multiphase Galactic Disks: Numerical Tests of the Thermal/Dynamical Equilibrium Model

The CO-H2 Conversion Factor in Disc Galaxies and Mergers

Galactic Spiral Shocks with Thermal Instability in Vertically Stratified Galactic Disks

Regulation of Star Formation Rates in Multiphase Galactic Disks: a Thermal/Dynamical Equilibrium Model

Star Formation and Gas Dynamics in Galactic Disks: Physical Processes and Numerical Models

Gravitational Runaway and Turbulence Driving in Star-Gas Galactic Disks

Models of Vertically Stratified Two-Phase ISM Disks with MRI-Driven Turbulence

Global Modeling of Spur Formation in Spiral Galaxies

Interstellar Turbulence Driving by Galactic Spiral Shocks

Can Nonlinear Hydromagnetic Waves Support a Self-Gravitating Cloud?