Cosmic ray transport on galactic scales depends on the detailed properties of the magnetized, multiphase interstellar medium (ISM). In this work, we post-process a high-resolution TIGRESS magnetohydrodynamic simulation modeling a local galactic disk
patch with a two-moment fluid algorithm for cosmic ray transport. We consider a variety of prescriptions for the cosmic rays, from a simple purely diffusive formalism with constant scattering coefficient, to a physically-motivated model in which the scattering coefficient is set by critical balance between streaming-driven Alfven wave excitation and damping mediated by local gas properties. We separately focus on cosmic rays with kinetic energies of $sim 1$ GeV (high-energy) and $sim 30$~MeV (low-energy), respectively important for ISM dynamics and chemistry. We find that simultaneously accounting for advection, streaming, and diffusion of cosmic rays is crucial for properly modeling their transport. Advection dominates in the high-velocity, low-density, hot phase, while diffusion and streaming are more important in higher density, cooler phases. Our physically-motivated model shows that there is no single diffusivity for cosmic-ray transport: the scattering coefficient varies by four or more orders of magnitude, maximal at density $n_mathrm{H} sim 0.01, mathrm{cm}^{-3}$. Ion-neutral damping of Alfven waves results in strong diffusion and nearly uniform cosmic ray pressure within most of the mass of the ISM. However, cosmic rays are trapped near the disk midplane by the higher scattering rate in the surrounding lower-density, higher-ionization gas. The transport of high-energy cosmic rays differs from that of low-energy cosmic rays, with less effective diffusion and greater energy losses for the latter.
In a companion paper, we develop a theory for the evolution of stellar wind driven bubbles in dense, turbulent clouds. This theory proposes that turbulent mixing at a fractal bubble-shell interface leads to highly efficient cooling, in which the vast
majority of the input wind energy is radiated away. This energy loss renders the majority of the bubble evolution momentum-driven rather than energy-driven, with expansion velocities and pressures orders of magnitude lower than in the classical Weaver77 solution. In this paper, we validate our theory with three-dimensional, hydrodynamic simulations. We show that extreme cooling is not only possible, but is generic to star formation in turbulent clouds over more than three orders of magnitude in density. We quantify the few free parameters in our theory, and show that the momentum exceeds the wind input rate by only a factor ~ 1.2-4. We verify that the bubble/cloud interface is a fractal with dimension ~ 2.5-2.7. The measured turbulent amplitude (v_t ~ 200-400 km/s) in the hot gas near the interface is shown to be consistent with theoretical requirements for turbulent diffusion to efficiently mix and radiate away most of the wind energy. The fraction of energy remaining after cooling is only 1-Theta ~ 0.1-0.01, decreasing with time, explaining observations that indicate low hot-gas content and weak dynamical effects of stellar winds.
Winds from massive stars have velocities of 1000 km/s or more, and produce hot, high pressure gas when they shock. We develop a theory for the evolution of bubbles driven by the collective winds from star clusters early in their lifetimes, which invo
lves interaction with the turbulent, dense interstellar medium of the surrounding natal molecular cloud. A key feature is the fractal nature of the hot bubbles surface. The large area of this interface with surrounding denser gas strongly enhances energy losses from the hot interior, enabled by turbulent mixing and subsequent cooling at temperatures T = 10^4-10^5 K where radiation is maximally efficient. Due to the extreme cooling, the bubble radius scales differently (R ~ t^1/2) from the classical Weaver77 solution, and has expansion velocity and momentum lower by factors of 10-10^2 at given R, with pressure lower by factors of 10^2 - 10^3. Our theory explains the weak X-ray emission and low shell expansion velocities of observed sources. We discuss further implications of our theory for observations of the hot bubbles and cooled expanding shells created by stellar winds, and for predictions of feedback-regulated star formation in a range of environments. In a companion paper, we validate our theory with a suite of hydrodynamic simulations.
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.
We explore the physics of the gyro-resonant cosmic ray streaming instability (CRSI) including the effects of ion-neutral (IN) damping. This is the main damping mechanism in (partially-ionized) atomic and molecular gas, which are the primary component
s of the interstellar medium (ISM) by mass. Limitation of CRSI by IN damping is important in setting the amplitude of Alfven waves that scatter cosmic rays and control galactic-scale transport. Our study employs the MHD-PIC hybrid fluid-kinetic numerical technique to follow linear growth as well as post-linear and saturation phases. During the linear phase of the instability -- where simulations and analytical theory are in good agreement -- IN damping prevents wave growth at small and large wavelengths, with the unstable bandwidth lower for higher ion-neutral collision rate $ u_{rm in}$. Purely MHD effects during the post-linear phase extend the wave spectrum towards larger $k$. In the saturated state, the cosmic ray distribution evolves toward greater isotropy (lower streaming velocity) by scattering off of Alven waves excited by the instability. In the absence of low-$k$ waves, CRs with sufficiently high momentum are not isotropized. The maximum wave amplitude and rate of isotropization of the distribution function decreases at higher $ u_{rm in}$. When the IN damping rate approaches the maximum growth rate of CSRI, wave growth and isotropization is suppressed. Implications of our results for CR transport in partially ionized ISM phases are discussed.
Molecular clouds are supported by turbulence and magnetic fields, but quantifying their influence on cloud lifecycle and star formation efficiency (SFE) remains an open question. We perform radiation MHD simulations of star-forming giant molecular cl
ouds (GMCs) with UV radiation feedback, in which the propagation of UV radiation via ray-tracing is coupled to hydrogen photochemistry. We consider 10 GMC models that vary in either initial virial parameter ($1lealpha_{v,0}le 5$) or dimensionless mass-to-magnetic flux ratio (0.5-8 and $infty$); the initial mass $10^5M_{odot}$ and radius 20pc are fixed. Each model is run with five different initial turbulence realizations. In most models, the duration of star formation and the timescale for molecular gas removal (primarily by photoevaporation) are 4-8Myr. Both the final SFE ($epsilon_*$) and time-averaged SFE per freefall time ($epsilon_{ff}$) are reduced by strong turbulence and magnetic fields. The median $epsilon_*$ ranges between 2.1% and 9.5%. The median $epsilon_{ff}$ ranges between 1.0% and 8.0% and anticorrelates with $alpha_{v,0}$, in qualitative agreement with previous analytic theory and simulations. However, the time-dependent $alpha_{v}(t)$ and $epsilon_{ff,obs}(t)$ based on instantaneous gas properties and cluster luminosity are positively correlated due to rapid evolution, making observational validation of star formation theory difficult. Our median $epsilon_{ff,obs}(t)approx$ 2% is similar to observed values. We show that the traditional virial parameter estimates the true gravitational boundedness within a factor of 2 on average, but neglect of magnetic support and velocity anisotropy can sometimes produce large departures. Magnetically subcritical GMCs are unlikely to represent sites of massive star formation given their unrealistic columnar outflows, prolonged lifetime, and low escape fraction of radiation.
Galactic outflows have density, temperature, and velocity variations at least as large as that of the multiphase, turbulent interstellar medium (ISM) from which they originate. We have conducted a suite of parsec-resolution numerical simulations usin
g the TIGRESS framework, in which outflows emerge as a consequence of interaction between supernovae (SNe) and the star-forming ISM. The outflowing gas is characterized by two distinct thermal phases, cool (T<10^4 K) and hot (T>10^6 K), with most mass carried by the cool phase and most energy and newly-injected metals carried by the hot phase. Both components have a broad distribution of outflow velocity, and especially for cool gas this implies a varying fraction of escaping material depending on the halo potential. Informed by the TIGRESS results, we develop straightforward analytic formulae for the joint probability density functions (PDFs) of mass, momentum, energy, and metal loading as distributions in outflow velocity and sound speed. The model PDFs have only two parameters, SFR surface density Sigma_SFR and the metallicity of the ISM, and fully capture the behavior of the original TIGRESS simulation PDFs over Sigma_SFR~(10^{-4},1)M_sun/kpc^2/yr. Employing PDFs from resolved simulations will enable galaxy formation subgrid model implementations with wind velocity and temperature (as well as total loading factors) that are based on theoretical predictions rather than empirical tuning. This is a critical step to incorporate advances from TIGRESS and other high-resolution simulations in future cosmological hydrodynamics and semi-analytic galaxy formation models. We release a python package to prototype our model and to ease its implementation.
CO is the most widely used observational tracer of molecular gas. The observable CO luminosity is translated to H_2 mass via a conversion factor, X_CO, which is a source of uncertainty and bias. Despite variations in X_CO, the empirically-determined
solar neighborhood value is often applied across different galactic environments. To improve understanding of X_CO, we employ 3D magnetohydrodynamics simulations of the interstellar medium (ISM) in galactic disks with a large range of gas surface densities, allowing for varying metallicity, far-ultraviolet (FUV) radiation, and cosmic ray ionization rate (CRIR). With the TIGRESS simulation framework we model the three-phase ISM with self-consistent star formation and feedback, and post-process outputs with chemistry and radiation transfer to generate synthetic CO(1--0) and (2--1) maps. Our models reproduce the observed CO excitation temperatures, line-widths, and line ratios in nearby disk galaxies. X_CO decreases with increasing metallicity, with a power-law slope of -0.8 for the (1--0) line and -0.5 for the (2--1) line. X_CO also decreases at higher CRIR, and is insensitive to the FUV radiation. As density increases, X_CO first decreases due to increasing excitation temperature, and then increases when the emission is fully saturated. We provide fits between X_CO and observable quantities such as the line ratio, peak antenna temperature, and line brightness, which probe local gas conditions. These fits, which allow for varying beam size, may be used in observations to calibrate out systematic biases. We also provide estimates of the CO-dark H_2 fraction at different gas surface densities, observational sensitivities, and beam sizes.
It has been hypothesized that photons from young, massive star clusters are responsible for maintaining the ionization of diffuse warm ionized gas seen in both the Milky Way and other disk galaxies. For a theoretical investigation of the warm ionized
medium (WIM), it is crucial to solve radiation transfer equations where the ISM and clusters are modeled self-consistently. To this end, we employ a Solar neighborhood model of TIGRESS, a magnetohydrodynamic simulation of the multiphase, star-forming ISM, and post-process the simulation with an adaptive ray tracing method to transfer UV radiation from star clusters. We find that the WIM volume filling factor is highly variable, and sensitive to the rate of ionizing photon production and ISM structure. The mean WIM volume filling factor rises to ~0.15 at |z|~1 kpc. Approximately half of ionizing photons are absorbed by gas and half by dust; the cumulative ionizing photon escape fraction is 1.1%. Our time-averaged synthetic H$alpha$ line profile matches WHAM observations on the redshifted (outflowing) side, but has insufficient intensity on the blueshifted side. Our simulation matches the Dickey-Lockman neutral density profile well, but only a small fraction of snapshots have high-altitude WIM density consistent with Reynolds Layer estimates. We compute a clumping correction factor C = <n_e>/sqrt<n_e^2>~0.2 that is remarkably constant with distance from the midplane and time; this can be used to improve estimates of ionized gas mass and mean electron density from observed H$alpha$ surface brightness profiles in edge-on galaxies.
Spiral arms greatly affect gas flows and star formation in disk galaxies. We use local three-dimensional simulations of the vertically-stratified, self-gravitating, differentially-rotating, interstellar medium (ISM) subject to a stellar spiral potent
ial 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.
