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

Woong-Tae Kim

Woong-Tae 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-phastro-ph.EPastro-ph.SRastro-ph.CO

Catalog footprint

What is connected

35works
5topics
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

35 published item(s)

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
preprint2021arXiv

OGLE-2019-BLG-0468Lb,c: two microlensing giant planets around a G-type star

With the aim of interpreting anomalous lensing events with no suggested models, we conducted a project of reinvestigating microlensing data in and before the 2019 season. In this work, we report a multi-planet system OGLE-2019-BLG-0468L found from the project. The light curve of the lensing event OGLE-2019-BLG-0468, which consists of three distinctive anomaly features, could not be explained by the usual binary-lens or binary-source interpretation. We find a solution explaining all anomaly features with a triple-lens interpretation, in which the lens is composed of two planets and their host, making the lens the fourth multi-planet system securely found by microlensing. The two planets have masses $\sim 3.4~M_{\rm J}$ and $\sim 10.2~M_{\rm J}$, and they are orbiting around a G-type star with a mass $\sim 0.9~M_\odot$ and a distance $\sim 4.4$ kpc. The host of the planets is most likely responsible for the light of the baseline object, although the possibility for the host to be a companion to the baseline object cannot be ruled out.

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

One Planet or Two Planets? The Ultra-sensitive Extreme-magnification Microlensing Event KMT-2019-BLG-1953

We present the analysis of a very high-magnification ($A\sim 900$) microlensing event KMT-2019-BLG-1953. A single-lens single-source (1L1S) model appears to approximately delineate the observed light curve, but the residuals from the model exhibit small but obvious deviations in the peak region. A binary lens (2L1S) model with a mass ratio $q\sim 2\times 10^{-3}$ improves the fits by $Δχ^2=181.8$, indicating that the lens possesses a planetary companion. From additional modeling by introducing an extra planetary lens component (3L1S model) and an extra source companion (2L2S model), it is found that the residuals from the 2L1S model further diminish, but claiming these interpretations is difficult due to the weak signals with $Δχ^2=16.0$ and $13.5$ for the 3L1S and 2L2L models, respectively. From a Bayesian analysis, we estimate that the host of the planets has a mass of $M_{\rm host}=0.31^{+0.37}_{-0.17}~M_\odot$ and that the planetary system is located at a distance of $D_{\rm L}=7.04^{+1.10}_{-1.33}~{\rm kpc}$ toward the Galactic center. The mass of the securely detected planet is $M_{\rm p}=0.64^{+0.76}_{-0.35}~M_{\rm J}$. The signal of the potential second planet could have been confirmed if the peak of the light curve had been more densely observed by followup observations, and thus the event illustrates the need for intensive followup observations for very high-magnification events even in the current generation of high-cadence surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

OGLE-2016-BLG-1227L: A Wide-separation Planet from a Very Short-timescale Microlensing Event

We present the analysis of the microlensing event OGLE-2016-BLG-1227. The light curve of this short-duration event appears to be a single-lens event affected by severe finite-source effects. Analysis of the light curve based on single-lens single-source (1L1S) modeling yields very small values of the event timescale, $t_{\rm E}\sim 3.5$ days, and the angular Einstein radius, $θ_{\rm E}\sim 0.009$ mas, making the lens a candidate of a free-floating planet. Close inspection reveals that the 1L1S solution leaves small residuals with amplitude $ΔI\lesssim 0.03$ mag. We find that the residuals are explained by the existence of an additional widely-separated heavier lens component, indicating that the lens is a wide-separation planetary system rather than a free-floating planet. From Bayesian analysis, it is estimated that the planet has a mass of $M_{\rm p} = 0.79^{+1.30}_{-0.39} M_{\rm J}$ and it is orbiting a low-mass host star with a mass of $M_{\rm host}=0.10^{+0.17}_{-0.05} M_\odot$ located with a projected separation of $a_\perp=3.4^{+2.1}_{-1.0}$ au. The planetary system is located in the Galactic bulge with a line-of-sight separation from the source star of $D_{\rm LS}=1.21^{+0.96}_{-0.63}$ kpc. The event shows that there are a range of deviations in the signatures of host stars for apparently isolated planetary lensing events and that it is possible to identify a host even when a deviation is subtle.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Origin of Non-axisymmetric Features of Virgo Cluster Early-type Dwarf Galaxies. II. Tidal Effects on Disk Features and Stability

A fraction of dwarf galaxies in the Virgo cluster contain disk features like bars and spiral arms. Using $N$-body simulations, we investigate the effects of tidal forces on the formation of such disk features in disk dwarf galaxies resembling VCC856. We consider 8 Cluster-Galaxy models in which disk dwarf galaxies with differing pericenter distance and spin orientation experience the tidal gravitational force of a Virgo-like NFW halo, and additional 8 Galaxy-Galaxy models in which two dwarf galaxies undergo tidal interactions with different strength. We find that the cluster tidal effect is moderate due to the small galaxy size, making the bars form earlier by $\sim1$--$1.5\Gyr$ compared to the cases in isolation. While the galactic halos significantly lose their mass within the virial radius due to the cluster tidal force, the mass of the stellar disks is nearly unchanged, suggesting that the inner regions of a disk-halo system is secured from the tidal force. The tidal forcing from either the cluster potential or a companion galaxy triggers the formation of two-armed spirals at early time before a bar develops. The tidally-driven arms decay and wind with time, suggesting that they are kinematic density waves. In terms of the strength and pitch angle, the faint arms in VCC856 are best matched with the arms in a marginally unstable galaxy produced by a distant tidal encounter with its neighbor $\sim0.85\Gyr$ ago.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Nuclear Filaments inside the Circumnuclear Disk in the Central 0.5 pc of the Galactic center

We present CS(7-6) line maps toward the central parsec of the Galactic Center (GC), conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The primary goal is to find and characterize the gas structure in the inner cavity of the circumnuclear disk (CND) in high resolution (1.3"=0.05 pc). Our large field-of-view mosaic maps -- combining interferometric and single-dish data that recover extended emission - provide a first homogeneous look to resolve and link the molecular streamers in the CND with the neutral nuclear filaments newly detected within the central cavity of the CND. We find that the nuclear filaments are rotating with Keplerian velocities in a nearly face-on orbit with an inclination angle of ~10-20 degree (radius <= 0.5 pc). This is in contrast to the CND which is highly inclined at ~65-80 degree (radius ~2-5 pc). Our analysis suggests a highly warped structure from the CND to the nuclear filaments. This result may hint that the nuclear filaments and the CND were created by different external clouds passing by Sgr A*.

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Equilibrium Sequences and Gravitational Instability of Rotating Isothermal Rings

Nuclear rings at centers of barred galaxies exhibit strong star formation activities. They are thought to undergo gravitational instability when sufficiently massive. We approximate them as rigidly-rotating isothermal objects and investigate their gravitational instability. Using a self-consistent field method, we first construct their equilibrium sequences specified by two parameters: alpha corresponding to the thermal energy relative to gravitational potential energy, and R_B measuring the ellipticity or ring thickness. Unlike in the incompressible case, not all values of R_B yield an isothermal equilibrium, and the range of R_B for such equilibria shrinks with decreasing alpha. The density distributions in the meridional plane are steeper for smaller alpha, and well approximated by those of infinite cylinders for slender rings. We also calculate the dispersion relations of nonaxisymmetric modes in rigidly-rotating slender rings with angular frequency Omega_0 and central density rho_max. Rings with smaller alpha are found more unstable with a larger unstable range of the azimuthal mode number. The instability is completely suppressed by rotation when Omega_0 exceeds the critical value. The critical angular frequency is found to be almost constant at ~ 0.7 sqrt(G*rho_c) for alpha > 0.01 and increases rapidly for smaller alpha. We apply our results to a sample of observed star-forming rings and confirm that rings without a noticeable azimuthal age gradient of young star clusters are indeed gravitationally unstable.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hydrodynamical Simulations of Nuclear Rings in Barred Galaxies

Dust lanes, nuclear rings, and nuclear spirals are typical gas structures in the inner region of barred galaxies. Their shapes and properties are linked to the physical parameters of the host galaxy. We use high-resolution hydrodynamical simulations to study 2D gas flows in simple barred galaxy models. The nuclear rings formed in our simulations can be divided into two groups: one group is nearly round and the other is highly elongated. We find that roundish rings may not form when the bar pattern speed is too high or the bulge central density is too low. We also study the periodic orbits in our galaxy models, and find that the concept of inner Lindblad resonance (ILR) may be generalized by the extent of $x_2$ orbits. All roundish nuclear rings in our simulations settle in the range of $x_2$ orbits (or ILRs). However, knowing the resonances is insufficient to pin down the exact location of these nuclear rings. We suggest that the backbone of round nuclear rings is the $x_2$ orbital family, i.e. round nuclear rings are allowed only in the radial range of $x_2$ orbits. A round nuclear ring forms exactly at the radius where the residual angular momentum of infalling gas balances the centrifugal force, which can be described by a parameter $f_{\rm ring}$ measured from the rotation curve. The gravitational torque on gas in high pattern speed models is larger, leading to a smaller ring size than in the low pattern speed models. Our result may have important implications for using nuclear rings to measure the parameters of real barred galaxies with 2D gas kinematics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Physical Properties of Tidal Features of Interacting Disk Galaxies: Three-dimensional Self-consistent Models

Using self-consistent three-dimensional (3D) N-body simulations, we investigate the physical properties of non-axisymmetric features in a disk galaxy created by a tidal interaction with its companion. The primary galaxy consists of a stellar disk, a bugle, and a live halo, corresponding to Milky-Way type galaxies, while the companion is represented by a halo alone. We vary the companion mass and the pericenter distance to explore situations with differing tidal strength parameterized by either the relative tidal force P or the relative imparted momentum S. We find that the formation of a tidal tail in the outer parts requires P > 0.05 or S > 0.07. A stronger interaction results in a stronger, less wound tail that forms earlier. Similarly, a stronger tidal forcing produces stronger, more loosely wound spiral arms in the inner parts. The arms are approximately logarithmic in shape, with both amplitude and pitch angle decaying with time. The derived pattern speed decreases with radius and is close to the Omega-kappa/2 curve at late time, with Omega and kappa denoting the angular and epicycle frequencies, respectively. This suggests that the tidally-induced spiral arms are most likely kinematic density waves weakly modified by self-gravity. Compared to the razor-thin counterparts, arms in the 3D models are weaker, have a smaller pitch angle, and wind and decay more rapidly. The 3D density structure of the arms is well described by the concentrated and sinusoidal models when the arms are in the nonlinear and linear regimes, respectively. We demonstrate that dynamical friction between interacting galaxies transfers the orbital angular momentum of one galaxy to the spin angular momentum of the companion halo.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Wiggle Instability of Galactic Spiral Shocks: Effects of Magnetic Fields

It has been suggested that the wiggle instability (WI) of spiral shocks in a galactic disk is responsible for the formation of gaseous feathers observed in grand-design spiral galaxies. We perform both a linear stability analysis and numerical simulations to investigate the effect of magnetic fields on the WI. The disk is assumed to be infinitesimally-thin, isothermal, and non-self-gravitating. We control the strengths of magnetic fields and spiral-arm forcing using the dimensionless parameters $β$ and $\mathcal{F}$, respectively. By solving the perturbation equations as a boundary-eigenvalue problem, we obtain dispersion relations of the WI for various values of $β=1-\infty$ and $\mathcal{F}=5\%$ and $10\%$. We find that the WI arising from the accumulation of potential vorticity at disturbed shocks is suppressed, albeit not completely, by magnetic fields. The stabilizing effect of magnetic fields is not from the perturbed fields but from the unperturbed fields that reduce the density compression factor in the background shocks. When $\mathcal{F}=5\%$ and $β\lesssim 10$ or $\mathcal{F}=10\%$ and $β\sim5-10$, the most unstable mode has a wavelength of $\sim0.1-0.2$ times the arm-to-arm separation, which appears consistent with a mean spacing of observed feathers.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Effects of Spiral Arms on Star Formation in Nuclear Rings of Barred-spiral Galaxies

We use hydrodynamic simulations to study the effect of spiral arms on the star formation rate (SFR) occurring in nuclear rings of barred-spiral galaxies. We find that spiral arms can be an efficient means of gas transport from the outskirts to the central parts, provided that the arms are rotating slower than the bar. While the ring star formation in models with no-arm or corotating arms is active only during about the bar growth phase, arm-driven gas accretion makes the ring star formation both enhanced and prolonged significantly in models with slow-rotating arms. The arm-enhanced SFR is larger by a factor of ~ 3-20 than in the no-arm model, with larger values corresponding to stronger and slower arms. Arm-induced mass inflows also make dust lanes stronger. Nuclear rings in slow-arm models are ~ 45% larger than in the no-arm counterparts. Star clusters that form in a nuclear ring exhibit an age gradient in the azimuthal direction only when the SFR is small, whereas no noticeable age gradient is found in the radial direction for models with arm-induced star formation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gaseous Spiral Structure and Mass Drift in Spiral Galaxies

We use hydrodynamic simulations to investigate nonlinear gas responses to an imposed stellar spiral potential in disk galaxies. The gaseous medium is assumed to be infinitesimally thin, isothermal, and unmagnetized. We consider various spiral-arm models with differing strength and pattern speed. We find that the extent and shapes of gaseous arms as well as the related mass drift rate depend rather sensitively on the arm pattern speed. In models where the arm pattern is rotating slow, the gaseous arms extend across the corotation resonance (CR) all the way to the outer boundary, with a pitch angle slightly smaller than that of the stellar counterpart. In models with a fast rotating pattern, on the other hand, spiral shocks are much more tightly wound than the stellar arms, and cease to exist in the regions near and outside the CR where $\mathcal{M}_\perp/{\rm sin} p_* \ge 25-40$, with $\mathcal{M}_\perp$ denoting the perpendicular Mach number of a rotating gas relative to the arms with pitch angle $p_*$. Inside the CR, the arms drive mass inflows at a rate of $\sim 0.05-3.0 {\rm M}_\odot {\rm yr}^{-1}$ to the central region, with larger values corresponding to stronger and slower arms. The contribution of the shock dissipation, external torque, and self-gravitational torque to the mass inflow is roughly 50%, 40%, and 10%, respectively. We demonstrate that the distributions of line-of-sight velocities and spiral-arm densities can be a useful diagnostic tool to distinguish if the spiral pattern is rotating fast or slow.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Instability of Magnetized Ionization Fronts Surrounding H II Regions

An ionization front (IF) surrounding an H II region is a sharp interface where a cold neutral gas makes transition to a warm ionized phase by absorbing UV photons from central stars. We investigate the instability of a plane-parallel D-type IF threaded by parallel magnetic fields, by neglecting the effects of recombination within the ionized gas. We find that weak D-type IFs always have the post-IF magnetosonic Mach number $\mathcal{M}_{\rm M2} \leq 1$. For such fronts, magnetic fields increase the maximum propagation speed of the IFs, while reducing the expansion factor $α$ by a factor of $1+1/(2β_1)$ compared to the unmagnetized case, with $β_1$ denoting the plasma beta in the pre-IF region. IFs become unstable to distortional perturbations due to gas expansion across the fronts, exactly analogous to the Darrieus-Landau instability of ablation fronts in terrestrial flames. The growth rate of the IF instability is proportional linearly to the perturbation wavenumber as well as the upstream flow speed, and approximately to $α^{1/2}$. The IF instability is stabilized by gas compressibility and becomes completely quenched when the front is D-critical. The instability is also stabilized by magnetic pressure when the perturbations propagate in the direction perpendicular to the fields. When the perturbations propagate in the direction parallel to the fields, on the other hand, it is magnetic tension that reduces the growth rate, completely suppressing the instability when $\mathcal{M}_{\rm M2}^2 < 2/(β_1 - 1)$. When the front experiences an acceleration, the IF instability cooperates with the Rayleigh-Taylor instability to make the front more unstable.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nature of the Wiggle Instability of Galactic Spiral Shocks

Gas in disk galaxies interacts nonlinearly with an underlying stellar spiral potential to form galactic spiral shocks. While numerical simulations typically show that spiral shocks are unstable to wiggle instability (WI) even in the absence of magnetic fields and self-gravity, its physical nature has remained uncertain. To clarify the mechanism behind the WI, we conduct a normal-mode linear stability analysis as well as nonlinear simulations assuming that the disk is isothermal and infinitesimally thin. We find that the WI is physical, originating from the generation of potential vorticity at a deformed shock front, rather than Kelvin-Helmholtz instabilities as previously thought. Since gas in galaxy rotation periodically passes through the shocks multiple times, the potential vorticity can accumulate successively, setting up a normal mode that grows exponentially with time. Eigenfunctions of the WI decay exponentially downstream from the shock front. Both shock compression of acoustic waves and a discontinuity of shear across the shock stabilize the WI. The wavelength and growth time of the WI depend on the arm strength quite sensitively. When the stellar-arm forcing is moderate at 5%, the wavelength of the most unstable mode is about 0.07 times the arm-to-arm spacing, with the growth rate comparable to the orbital angular frequency, which is found to be in good agreement with the results of numerical simulations.

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

Instability of Evaporation Fronts in the Interstellar Medium

The neutral component of the interstellar medium is segregated into the cold neutral medium (CNM) and warm neutral medium (WNM) as a result of thermal instability. It was found that a plane-parallel CNM-WNM evaporation interface, across which the CNM undergoes thermal expansion, is linearly unstable to corrugational disturbances, in complete analogy with the Darrieus-Landau instability (DLI) of terrestrial flames. We perform a full linear stability analysis as well as nonlinear hydrodynamic simulations of the DLI of such evaporation fronts in the presence of thermal conduction. We find that the DLI is suppressed at short length scales by conduction. The length and time scales of the fastest growing mode are inversely proportional to the evaporation flow speed of the CNM and its square, respectively. In the nonlinear stage, the DLI saturates to a steady state where the front deforms to a finger-like shape protruding toward the WNM, without generating turbulence. The evaporation rate at nonlinear saturation is larger than the initial plane-parallel value by a factor of 2.4 when the equilibrium thermal pressure is 1800 k_B cm^-3 K. The degrees of front deformation and evaporation-rate enhancement at nonlinear saturation are determined primarily by the density ratio between the CNM and WNM. We demonstrate that the Field length in the thermally unstable medium should be resolved by at least four grid points to obtain reliable numerical outcomes involving thermal instability.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Star Formation in Nuclear Rings of Barred Galaxies

Nuclear rings in barred galaxies are sites of active star formation. We use hydrodynamic simulations to study temporal and spatial behavior of star formation occurring in nuclear rings of barred galaxies where radial gas inflows are triggered solely by a bar potential. The star formation recipes include a density threshold, an efficiency, conversion of gas to star particles, and delayed momentum feedback via supernova explosions. We find that star formation rate (SFR) in a nuclear ring is roughly equal to the mass inflow rate to the ring, while it has a weak dependence on the total gas mass in the ring. The SFR typically exhibits a strong primary burst followed by weak secondary bursts before declining to very small values. The primary burst is associated with the rapid gas infall to the ring due to the bar growth, while the secondary bursts are caused by re-infall of the ejected gas from the primary burst. While star formation in observed rings persists episodically over a few Gyr, the duration of active star formation in our models lasts for only about a half of the bar growth time, suggesting that the bar potential alone is unlikely responsible for gas supply to the rings. When the SFR is low, most star formation occurs at the contact points between the ring and the dust lanes, leading to an azimuthal age gradient of young star clusters. When the SFR is large, on the other hand, star formation is randomly distributed over the whole circumference of the ring, resulting in no apparent azimuthal age gradient. Since the ring shrinks in size with time, star clusters also exhibit a radial age gradient, with younger clusters found closer to the ring. The cluster mass function is well described by a power law, with a slope depending on the SFR. Giant gas clouds in the rings have supersonic internal velocity dispersions and are gravitationally bound.

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

Gaseous Structures in Barred Galaxies: Effects of the Bar Strength

Using hydrodynamic simulations, we investigate the physical properties of gaseous substructures in barred galaxies and their relationships with the bar strength. The gaseous medium is assumed to be isothermal and unmagnetized. The bar potential is modeled as a Ferrers prolate with index n. To explore situations with differing bar strength, we vary the bar mass fbar relative to the spheroidal component as well as its aspect ratio. We derive expressions as functions of fbar and the aspect ratio for the bar strength Qb and the radius r(Qb) where the maximum bar torque occurs. When applied to observations, these expressions suggest that bars in real galaxies are most likely to have fbar=0.25-0.5 and n<1. Dust lanes approximately follow one of x1-orbits and tend to be more straight under a stronger and more elongated bar, but are insensitive to the presence of self-gravity. A nuclear ring of a conventional x2 type forms only when the bar is not so massive or elongated. The radius of an x2-type ring is generally smaller than the inner Lindblad resonance, decreases systematically with increasing Qb, and slightly larger when self-gravity is included. This evidences that the ring position is not determined by the resonance but by the amount of angular momentum loss at dust-lane shocks. Nuclear spirals exist only when the ring is of the x2-type and sufficiently large in size. Unlike the other features, nuclear spirals are transient in that they start out as being tightly-wound and weak, and then due to the nonlinear effect unwind and become stronger until turning into shocks, with an unwinding rate higher for larger Qb. The mass inflow rate to the galaxy center is found to be less than 0.01 Msun/yr for models with Qb<0.2, while becoming larger than 0.1 Msun/yr when Qb>0.2 and self-gravity is included.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gravitational Instability of Rotating, Pressure-Confined, Polytropic Gas Disks With Vertical Stratification

We investigate gravitational instability (GI) of rotating, vertically-stratified, pressure-confined, polytropic gas disks using linear stability analysis as well as analytic approximations. The disks are initially in vertical hydrostatic equilibrium and bounded by a constant external pressure. We find that GI of a pressure-confined disk is in general a mixed mode of the conventional Jeans and distortional instabilities, and is thus an unstable version of acoustic-surface-gravity waves. The Jeans mode dominates in weakly confined disks or disks with rigid boundaries. When the disk has free boundaries and is strongly pressure-confined, on the other hand, the mixed GI is dominated by the distortional mode that is surface-gravity waves driven unstable under own gravity and thus incompressible. We demonstrate that the Jeans mode is gravity-modified acoustic waves rather than inertial waves and that inertial waves are almost unaffected by self-gravity. We derive an analytic expression for the effective sound speed c_eff of acoustic-surface-gravity waves. We also find expressions for the gravity reduction factors relative to a razor-thin counterpart, appropriate for the Jeans and distortional modes. The usual razor-thin dispersion relation after correcting for c_eff and the reduction factors closely matches the numerical results obtained by solving a full set of linearized equations. The effective sound speed generalizes the Toomre stability parameter of the Jeans mode to allow for the mixed GI of vertically-stratified, pressure-confined disks.

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
preprint2012arXiv

Two-Dimensional Magnetohydrodynamic Simulations of Barred Galaxies

Barred galaxies are known to possess magnetic fields that may affect the properties of bar substructures such as dust lanes and nuclear rings. We use two-dimensional high-resolution magnetohydrodynamic (MHD) simulations to investigate the effects of magnetic fields on the formation and evolution of such substructures as well as on the mass inflow rates to the galaxy center. The gaseous medium is assumed to be infinitesimally-thin, isothermal, non-self-gravitating, and threaded by initially uniform, azimuthal magnetic fields. We find that there exists an outermost x1-orbit relative to which gaseous responses to an imposed stellar bar potential are completely different between inside and outside. Inside this orbit, gas is shocked into dust lanes and infalls to form a nuclear ring. Magnetic fields are compressed in dust lanes, reducing their peak density. Magnetic stress removes further angular momentum of the gas at the shocks, temporarily causing the dust lanes to bend into an 'L' shape and eventually leading to a smaller and more centrally distributed ring than in unmagnetized models. The mass inflow rates in magnetized models correspondingly become larger, by more than two orders of magnitude when the initial fields have an equipartition value with thermal energy, than in the unmagnetized counterparts. Outside the outermost x1-orbit, on the other hand, an MHD dynamo due to the combined action of the bar potential and background shear operates near the corotation and bar-end regions, efficiently amplifying magnetic fields. The amplified fields shape into trailing magnetic arms with strong fields and low density. The base of the magnetic arms has a thin layer in which magnetic fields with opposite polarity reconnect via a tearing-mode instability. This produces numerous magnetic islands with large density which propagate along the arms to turn the outer disk into a highly chaotic state.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Central Regions of Barred Galaxies: Two-Dimensional Non-self-gravitating Hydrodynamic Simulations

The inner regions of barred galaxies contain substructures such as off-axis shocks, nuclear rings, and nuclear spirals. These substructure may affect star formation, and control the activity of a central black hole (BH) by determining the mass inflow rate. We investigate the formation and properties of such substructures using high-resolution, grid-based hydrodynamic simulations. The gaseous medium is assumed to be infinitesimally-thin, isothermal, and non-self-gravitating. The stars and dark matter are represented by a static gravitational potential with four components: a stellar disk, the bulge, a central BH, and the bar. To investigate various galactic environments, we vary the gas sound speed c_s as well as the mass of the central BH M_BH. Once the flow has reached a quasi-steady state, off-axis shocks tend to move closer to the bar major axis as c_s increases. Nuclear rings shrink in size with increasing c_s, but are independent of M_BH, suggesting that ring position is not determined by the Lindblad resonances. Rings in low-c_s models are narrow since they are occupied largely by gas on x2-orbits and well decoupled from nuclear spirals, while they become broad because of large thermal perturbations in high-c_s models. Nuclear spirals persist only when either c_s is small or M_BH is large; they would otherwise be destroyed completely by the ring material on eccentric orbits. The shape and strength of nuclear spirals depend sensitively on c_s and M_BH such that they are leading if both c_s and M_BH are small, weak trailing if c_s is small and M_BH is large, and strong trailing if both c_s and M_BH are large. While the mass inflow rate toward the nucleus is quite small in low-c_s models because of the presence of a narrow nuclear ring, it becomes larger than 0.01 Msun/yr when c_s is large, providing a potential explanation of nuclear activity in Seyfert galaxies.

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

A Semi-empirical Mass-loss Rate in Short-period Cataclysmic Variables

The mass-loss rate of donor stars in cataclysmic variables (CVs) is of paramount importance in the evolution of short-period CVs. Observed donors are oversized in comparison with those of isolated single stars of the same mass, which is thought to be a consequence of the mass loss. Using the empirical mass-radius relation of CVs and the homologous approximation for changes in effective temperature T_2, orbital period P, and luminosity of the donor with the stellar radius, we find the semi-empirical mass-loss rate M2_dot of CVs as a function of P. The derived M2_dot is at ~10^(-9.5)-10^(-10) Msun/yr and depends weakly on P when P > 90 min, while it declines very rapidly towards the minimum period when P < 90 min, emulating the P-T_2 relation. Due to strong deviation from thermal equilibrium caused by the mass loss, the semi-empirical M2_dot is significantly different from, and has a less-pronounced turnaround behavior with P than suggested by previous numerical models. The semi-empirical P-M2_dot relation is consistent with the angular momentum loss due to gravitational wave emission, and strongly suggests that CV secondaries with 0.075 Msun < M_2 < 0.2 Msun are less than 2 Gyrs old. When applied to selected eclipsing CVs, our semi-empirical mass-loss rates are in good agreement with the accretion rates derived from the effective temperatures T_1 of white dwarfs, suggesting that M2_dot can be used to reliably infer T_2 from T_1. Based on the semi-empirical M2_dot, SDSS 1501 and 1433 systems that were previously identified as post-bounce CVs have yet to reach the minimal period.

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
preprint2010arXiv

Nonlinear Dynamical Friction of a Circular-Orbit Perturber in a Gaseous Medium

We use three-dimensional hydrodynamic simulations to investigate the nonlinear gravitational responses of gas to, and the resulting drag forces on, very massive perturbers moving on circular orbits. This work extends our previous studies that explored the cases of low-mass perturbers on circular orbits and massive perturbers on straight-line trajectories. The background medium is assumed to be non-rotating, adiabatic with index 5/3, and uniform with density rho0 and sound speed a0. We model the gravitating perturber using a Plummer sphere with mass Mp and softening radius rs in a uniform circular motion at speed Vp and orbital radius Rp, and run various models with differing R=rs/Rp, Mach=Vp/a0, and B=G*Mp/(a0^2*Rp). A quasi-steady density wake of a supersonic model consists of a hydrostatic envelope surrounding the perturber, an upstream bow shock, and a trailing low-density region. The continuous change in the direction of the perturber motion makes the detached shock distance reduced compared to the linear-trajectory cases, while the orbit-averaged gravity of the perturber gathers the gas toward the center of the orbit, modifying the background preshock density to rho1=(1 + 0.46B)*rho0 depending weakly on Mach. For sufficiently massive perturbers, the presence of a hydrostatic envelope makes the drag force smaller than the prediction of the linear perturbation theory, resulting in F = 4*pi*rho1*(G*Mp/Vp)^2 * (0.7/etaB) for etaB = B/(Mach^2 -1) > 0.1; the drag force for low-mass perturbers with etaB < 0.1 agrees well with the linear prediction. The nonlinear drag force becomes independent of R as long as R < etaB/2, which places an upper limit on the perturber size for accurate evaluation of the drag force in numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Internal Structure and Apsidal Motions of Polytropic Stars in Close Binaries

We consider a synchronized, circular-orbit binary consisting of a polytrope with index n and a point-mass object, and use a self-consistent field method to construct the equilibrium structure of the polytrope under rotational and tidal perturbations. Our self-consistent field method is distinct from others in that the equilibrium orbital angular velocity is calculated automatically rather than being prescribed, which is crucial for obtaining apsidal motion rates accurately. We find that the centrifugal and tidal forces make perturbed stars more centrally condensed and larger in size. For n=1.5 polytopes with fixed entropy, the enhancement factor in stellar radii is about 23% and 4-8% for mu=1 and sim0.1-0.9, respectively, where mu is the fractional mass of the polytrope relative to the total. The centrifugal force dominates the tidal force in determining the equilibrium structure provided mu > 0.13-0.14 for n > 1.5. The shape and size of rotationally- and tidally-perturbed polytropes are well described by the corresponding Roche models as long as n > 2. The apsidal motion rates calculated for circular-orbit binaries under the equilibrium tide condition agree well with the predictions of the classical formula only when the rotational and tidal perturbations are weak. When the perturbations are strong as in critical configurations, the classical theory underestimates the real apsidal motion rates by as much as 50% for n=1.5 polytropes, although the discrepancy becomes smaller as n increases. For practical uses, we provide fitting formulae for various quantities including the density concentration, volume radius, and effective internal structure constant, as functions of mu and the perturbation parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Nonlinear Dynamical Friction in a Gaseous Medium

Using high-resolution, two-dimensional hydrodynamic simulations, we investigate nonlinear gravitational responses of gas to, and the resulting drag force on, a very massive perturber M_p moving at velocity V_p through a uniform gaseous medium of adiabatic sound speed a_0. We model the perturber as a Plummer potential with softening radius r_s, and run various models with differing A=GM_p/(a_0^2 r_s) and M=V_p/a_0 by imposing cylindrical symmetry with respect to the line of perturber motion. For supersonic cases, a massive perturber quickly develops nonlinear flows that produce a detached bow shock and a vortex ring, which is unlike in the linear cases where Mach cones are bounded by low-amplitude Mach waves. The flows behind the shock are initially non-steady, displaying quasi-periodic, overstable oscillations of the vortex ring and the shock. The vortex ring is eventually shed downstream and the flows evolve toward a quasi-steady state where the density wake near the perturber is in near hydrostatic equilibrium. We find that the detached shock distance $δ$ and the nonlinear drag force F depend solely on η=A/(M^2-1) such that δ/r_s=ηand F/F_{lin}=(η/2)^{-0.45} for η>2, where F_{lin} is the linear drag force of Ostriker (1999). The reduction of F compared with F_{lin} is caused by front-back symmetry in the nonlinear density wakes. In subsonic cases, the flows without involving a shock do not readily reach a steady state. Nevertheless, the subsonic density wake near a perturber is close to being hydrostatic, resulting in the drag force similar to the linear case. Our results suggest that dynamical friction of a very massive object as in a merger of black holes near a galaxy center will take considerably longer than the linear prediction.

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

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
preprint2005arXiv

Dynamical Friction and Cooling Flows in Galaxy Clusters

We investigate a model of galaxy clusters in which the hot intracluster gas is efficiently heated by dynamical friction (DF) of galaxies. We allow for both subsonic and supersonic motions of galaxies and use the gravitational drag formula in a gaseous medium presented by Ostriker (1999). The energy lost by the galaxies is either redistributed locally or into a Gaussian centered on the galaxy. We find that the condition of hydrostatic equilibrium and strict energy balance yields a trivial isothermal solution T_iso, independent of radius, or rising temperature distributions provided T_iso/gamma < T < T_iso, where gamma is the adiabatic index of the gas. The isothermal temperature corresponds to the usual scaling relation between the gas temperatures and the velocity dispersions of galaxies. However the minimal temperature associated with the rising solutions is ~ 0.5 T_vir, larger than that inferred from observations, the radial distribution of galaxy masses notwithstanding. Heating by supersonically moving galaxies cannot suppress thermal instability, although it can lengthen the growth time up to the level comparable to the ages of clusters when Mach number of galaxies is less than about two. We show using numerical hydrodynamic simulations that DF-induced heating is generally unable to produce stable equilibrium cores by evolving arbitrary non-equilibrium clusters, although it can lengthen the cooling time. We conclude that DF-induced heating alone is an unlikely solution to the cooling flow problem, although it can still be an important heat supplier, considerably delaying cooling catastrophe. We discuss other potential consequences of DF of galaxies in galaxy clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Dynamical-Friction Galaxy-Gas Coupling and Cluster Cooling Flows

We revisit the notion that galaxy motions can efficiently heat intergalactic gas in the central regions of clusters through dynamical friction. For plausible values of the galaxy mass-to-light ratio, the heating rate is comparable to the cooling rate due to X-ray emission. Heating occurs only for supersonic galaxy motions, so the mechanism is self-regulating: it becomes efficient only when the gas sound speed is smaller than the galaxy velocity dispersion. We illustrate with the Perseus cluster, assuming a stellar mass-to-light ratio for galaxies in the very central region with the dark-matter contribution becoming comparable to this at some radius $r_s$. For $r_s \la 400 {\rm kpc} \sim 3 r_{\rm cool}$--corresponding to an average mass-to-light ratio of $\sim10$ inside that radius--the dynamical-friction coupling is strong enough to provide the required rate of gas heating. The measured sound speed is smaller than the galaxy velocity dispersion, as required by this mechanism. With this smaller gas temperature and the observed distribution of galaxies and gas, the energy reservoir in galactic motions is sufficient to sustain the required heating rate for the lifetime of the cluster. The galaxies also lose a smaller amount of energy through dynamical friction to the dark matter implying that non--cooling-flow clusters should have flat-cored dark-matter density distributions.

0 citations0 reviewsNo code linkedOpen access