BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Chang-Goo Kim

Chang-Goo Kim appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph

Catalog footprint

What is connected

20works
2topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

20 published item(s)

preprint2026arXiv

The impact of ram pressure on the radio spectral index and magnetic field of NGC 4522: A high-resolution VLA continuum study

We present high-resolution Very Large Array (VLA) continuum observations at S-band ($3$ GHz, $560$ pc scale) and X-band ($10$ GHz, $200$ pc scale) of the ram-pressure-stripped Virgo galaxy NGC 4522, to investigate the characteristics of its radio continuum, spectral index, and magnetic field under the influence of the intracluster medium (ICM). The total radio continuum shows an asymmetry that extends northwest, mirroring the HI gas distribution, but showing distinct features in the extraplanar regions. The spectral index steepens systematically from $α\sim-0.6$ in the main disk to $α\sim-1.1$ in the outer disk. We find that the spectral index behavior of the outer disk is mainly due to an ICM shock that can re-accelerate electrons and a significant reduction of thermal emission. Intriguingly, extraplanar clouds exhibit exceptionally flat spectral indices ($α\sim-0.2$ to $0$), resulting from a combination of significantly enhanced thermal emission and pronounced spectral aging of the non-thermal component. Although some of these regions correlate with H$α$, others do not. We propose that the mixing between the ICM and interstellar medium (ISM) is an alternative mechanism that enhances thermal emission independently of star formation. Polarized continuum emissions are highly asymmetric, preferentially distributed along the ICM wind side, and the polarization fraction increases radially outward from the galactic midplane, indicating that the polarized emission is strongly influenced by the ICM wind. Our results show how and where the ICM substantially affects the ISM, and also demonstrate that high-frequency observations are crucial for analyzing the radio continuum of ram pressure stripping galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Ram pressure stripping of the multiphase ISM: a detailed view from TIGRESS simulations

Ram pressure stripping (RPS) is a process that removes the interstellar medium (ISM) quickly, playing a vital role in galaxy evolution. Previous RPS studies have treated the ISM as single-phase or lack the resolution and physical processes to properly capture the full multiphase ISM. To improve this simplification, we introduce an inflowing, hot intracluster medium (ICM) into a self-consistently modeled ISM in a local patch of star-forming galactic disks using the TIGRESS framework. Our simulations reveal that the workings of RPS are not only direct acceleration of the ISM by ICM ram pressure but also mixing-driven momentum transfer involving significant phase transition and radiative cooling. The hot ICM passes through the low-density channels of the porous, multiphase ISM, shreds the cool ISM, and creates mixing layers. The ICM momentum is transferred through the mixing layers while populating the intermediate temperature gas and radiating thermal energy away. The mixed gas extends beyond galactic disks and forms stripped tails that cool back unless the ICM fluxes are large enough to prevent cooling until they escape the simulation domain. The mixing-driven momentum transfer predicts that the more ICM mixes in, the faster the ISM moves, resulting in the anti-correlation of outflow velocity and gas metallicity of the stripped ISM. The compression of the ISM disks due to the ICM ram pressure enhances star formation rates up to 50% compared to the model without ICM. With the ICM ram pressure higher than the disk anchoring pressure, star formation is quenched within ~100 Myr.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Extracting the cold neutral medium from HI emission with deep learning: Implications for Galactic foregrounds at high latitude

Resolving the phase structure of neutral hydrogen (HI) is crucial for understanding the life cycle of the interstellar medium (ISM). However, accurate measurements of HI temperature and density are limited by the availability of background continuum sources for measuring HI absorption. Here we test the use of deep learning for extracting HI properties over large areas without optical depth information. We train a 1D convolutional neural network using synthetic observations of 3D numerical simulations of the ISM to predict the fraction of cold neutral medium (f_CNM) and the correction to the optically-thin HI column density for optical depth (R_HI) from $21\rm\,cm$ emission alone. We restrict our analysis to high Galactic latitudes ($|b|>30 deg), where the complexity of spectral line profiles is minimized. We verify that the network accurately predicts f_CNM and R_HI by comparing the results with direct constraints from 21 cm absorption. By applying the network to the GALFA-HI survey, we generate large-area maps of f_CNM and R_HI. Although the overall contribution to the total HI column of cold neutral medium (CNM)-rich structures is small (~5%), we find that these structures are ubiquitous. Our results are consistent with the picture that small-scale structures observed in 21 cm emission aligned with the magnetic field are dominated by CNM. Finally, we demonstrate that the observed correlation between HI column density and dust reddening (E(B-V)) declines with increasing R_HI, indicating that future efforts to quantify foreground Galactic E(B-V) using HI, even at high latitudes, should increase fidelity by accounting for HI phase structure.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Ly-alpha Radiative Transfer: Monte-Carlo Simulation of the Wouthuysen-Field Effect

A three-dimensional Monte Carlo Ly-alpha radiative transfer (RT) code, named LaRT, is developed to study the Ly-alpha RT and the Wouthuysen-Field (WF) effect. Using the code, we calculate the line profile of Ly-alpha radiation within the multiphase interstellar medium (ISM), with a particular emphasis on gas at low densities. We show that the WF effect is in action: the central portion of the line profile tends to approach a small slice of the Planck function with a color temperature equal to the kinetic temperature of the gas, even in a system with an optical thickness as low as tau_0~100-500. We also investigate the effects of the turbulent motion of the ISM on the emergent Ly-alpha spectrum and color temperature. The turbulent motion broadens, as generally expected, the emergent spectrum, but the color temperature is not affected by the turbulent motion in typical astrophysical environments. We utilize two multiphase ISM models, appropriate for the vicinity of the Sun, to calculate the 21-cm spin temperature of neutral hydrogen, including excitation via the Ly-alpha resonant scattering. The first ISM model is a simple clumpy model, while the second is a self-consistent magnetohydrodynamics simulation model using the TIGRESS framework. Ly-alpha photons originating from both H II regions and the collisionally cooling gas are taken into account. We find that the Ly-alpha radiation field is, in general, likely to be strong enough to bring the 21-cm spin temperature of the warm neutral medium close to the kinetic temperature. The escape fraction of Ly-alpha in our ISM models is estimated to be ~7-20%.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access