Search

Nuclear rings at the centers of barred galaxies are known to be strongly magnetized. To explore the effects of magnetic fields on star formation in these rings and nuclear gas flows, we run magnetohydrodynamic simulations in which there is a temporally-constant magnetized inflow to the ring, representing a bar-driven inflow. The mass inflow rate is $1\,M_\odot\,\mathrm{yr}^{-1}$, and we explore models with a range of field strength in the inflow. We adopt the TIGRESS framework developed by Kim & Ostriker to handle radiative heating and cooling, star formation, and resulting supernova (SN) feedback. We find that magnetic fields are efficiently amplified in the ring due to rotational shear and SN feedback. Within a few $100\,\mathrm{Myr}$, the turbulent component $B_\mathrm{trb}$ in the ring saturates at $\sim 35\,\mu\mathrm{G}$ (in rough equipartition with the turbulent kinetic energy density), while the regular component $B_\mathrm{reg}$ exceeds $50\,\mu\mathrm{G}$. Expanding superbubbles created by clustered SN explosions vertically drag predominantly-toroidal fields from near the midplane to produce poloidal fields in high-altitude regions. The growth of magnetic fields greatly suppresses star formation at late times. Simultaneously, strong magnetic tension in the ring drives radially inward accretion flows from the ring to form a circumnuclear disk in the central region; this feature is absent in the unmagnetized model., Comment: 29 pages, 17 figures, accepted for publication in ApJ

While bars are common in disk galaxies, their formation conditions are not well understood. We use $N$-body simulations to study bar formation and evolution in isolated galaxies consisting of a stellar disk, a classical bulge, and a dark halo. We consider 24 galaxy models that are similar to the Milky Way but differ in the mass and compactness of the classical bulge and halo concentration. We find that the bar formation requires $(Q_{T,\text{min}}/1.2)^2 + (\text{CMC}/0.05)^2\lesssim 1$, where $Q_{T,\text{min}}$ and CMC refers to the minimum value of the Toomre stability parameter and the central mass concentration, respectively. Bars tend to be stronger, longer, and rotate slower in galaxies with a less massive and less compact bulge and halo. All bars formed in our models correspond to slow bars. A model with the bulge mass of $\sim10$--$20$\% of the disk under a concentrated halo produces a bar similar to the Milky-Way bar. We discuss our findings in relation to other bar formation criteria suggested by previous studies., Accepted for publication in ApJ

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 $\Delta \tau_\text{in}$ and amplitude 20 causes large-amplitude (a factor of $\gtrsim 5$), quasi-periodic variations of the SFR, when $\Delta \tau_\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., Comment: 15 pages, 6 figures, accepted for publication in ApJ. Replaced with minor text corrections

We investigate the gravitational microlensing event KMT-2019-BLG-1715, of which light curve shows two short-term anomalies from a caustic-crossing binary-lensing light curve: one with a large deviation and the other with a small deviation. We identify five pairs of solutions, in which the anomalies are explained by adding an extra lens or source component in addition to the base binary-lens model. We resolve the degeneracies by applying a method, in which the measured flux ratio between the first and second source stars is compared with the flux ratio deduced from the ratio of the source radii. Applying this method leaves a single pair of viable solutions, in both of which the major anomaly is generated by a planetary-mass third body of the lens, and the minor anomaly is generated by a faint second source. A Bayesian analysis indicates that the lens comprises three masses: a planet-mass object with $\sim 2.6~M_{\rm J}$ and binary stars of K and M dwarfs lying in the galactic disk. We point out the possibility that the lens is the blend, and this can be verified by conducting high-resolution followup imaging for the resolution of the lens from the source., Comment: 11 pages, 13 figures, 7 tables

Nuclear rings are sites of intense star formation at the centers of barred galaxies. To understand what determines the structure and star formation rate (SFR; $\dot{M}_{\rm SF}$) of nuclear rings, we run semi-global, hydrodynamic simulations of nuclear rings subject to constant mass inflow rates $\dot{M}_{\rm in}$. We adopt the TIGRESS framework of Kim \& Ostriker to handle radiative heating and cooling, star formation, and related supernova (SN) feedback. We find that the SN feedback is never strong enough to destroy the ring or quench star formation everywhere in the ring. Under the constant $\dot{M}_{\rm in}$, the ring star formation is very steady and persistent, with the SFR exhibiting only mild temporal fluctuations. The ring SFR is tightly correlated with the inflow rate as $\dot{M}_{\rm SF}\approx 0.8\dot{M}_{\rm in}$, for a range of $\dot{M}_{\rm in}=0.125-8\,M_\odot\,{\rm yr}^{-1}$. Within the ring, vertical dynamical equilibrium is maintained, with the midplane pressure (powered by SN feedback) balancing the weight of the overlying gas. The SFR surface density is correlated nearly linearly with the midplane pressure, as predicted by the pressure-regulated, feedback-modulated star formation theory. Based on our results, we argue that the ring SFR is causally controlled by $\dot{M}_\text{in}$, while the ring gas mass adapts to the SFR to maintain the vertical dynamical equilibrium under the gravitational field arising from both gas and stars., 31 pages, 22 figures, accepted for publication in ApJ

Utilizing the Atacama Large Millimeter/submillimeter Array (ALMA), we present CS line maps in five rotational lines ($J_{\rm u}=7, 5, 4, 3, 2$) toward the circumnuclear disk (CND) and streamers of the Galactic Center. Our primary goal is to resolve the compact structures within the CND and the streamers, in order to understand the stability conditions of molecular cores in the vicinity of the supermassive black hole (SMBH) Sgr A*. Our data provide the first homogeneous high-resolution ($1.3'' = 0.05$ pc) observations aiming at resolving density and temperature structures. The CS clouds have sizes of $0.05-0.2$ pc with a broad range of velocity dispersion ($\sigma_{\rm FWHM}=5-40$ km s$^{-1}$). The CS clouds are a mixture of warm ($T_{\rm k}\ge 50-500$ K, n$_{\rm H_2}$=$10^{3-5}$ cm$^{-3}$) and cold gas ($T_{\rm k}\le 50$ K, n$_{\rm H_2}$=$10^{6-8}$ cm$^{-3}$). A stability analysis based on the unmagnetized virial theorem including tidal force shows that $84^{+16}_{-37}$ % of the total gas mass is tidally stable, which accounts for the majority of gas mass. Turbulence dominates the internal energy and thereby sets the threshold densities $10-100$ times higher than the tidal limit at distance $\ge 1.5$ pc to Sgr A*, and therefore, inhibits the clouds from collapsing to form stars near the SMBH. However, within the central $1.5$ pc, the tidal force overrides turbulence and the threshold densities for a gravitational collapse quickly grow to $\ge 10^{8}$ cm$^{-3}$., Comment: accepted for publication in ApJ

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., Comment: 10 pages, 3 tables, 11 figures

Microlensing planets occurring on faint source stars can escape detection due to their weak signals. Occasionally, detections of such planets are not reported due to the difficulty of extracting high-profile scientific issues on the detected planets. For the solid demographic census of microlensing planetary systems based on a complete sample, we investigate the microlensing data obtained in the 2016 and 2017 seasons to search for planetary signals in faint-source lensing events. From this investigation, we find four unpublished microlensing planets including KMT-2016-BLG-2364Lb, KMT-2016-BLG-2397Lb, OGLE-2017-BLG-0604Lb, and OGLE-2017-BLG-1375Lb. We analyze the observed lensing light curves and determine their lensing parameters. From Bayesian analyses conducted with the constraints from the measured parameters, it is found that the masses of the hosts and planets are in the ranges $0.50\lesssim M_{\rm host}/M_\odot\lesssim 0.85$ and $0.5 \lesssim M_{\rm p}/M_{\rm J}\lesssim 13.2$, respectively, indicating that all planets are giant planets around host stars with subsolar masses. The lenses are located in the distance range of $3.8 \lesssim \dl/{\rm kpc}\lesssim 6.4$. It is found that the lenses of OGLE-2017-BLG-0604 and OGLE-2017-BLG-1375 are likely to be in the Galactic disk., 11 pages, 8 tables, 12 figures

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., Accepted for publication in the ApJ

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., Accepted for publication in the ApJ

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, 6 pages, 4 figures, accepted in ApJL

We run self-consistent simulations of Milky Way-sized, isolated disk galaxies to study formation and evolution of a stellar bar as well as a nuclear ring in the presence of gas. We consider two sets of models with cold or warm disks that differ in the radial velocity dispersions, and vary the gas fraction $f_{\rm gas}$ by fixing the total disk mass. A bar forms earlier and more strongly in the cold disks with larger $f_{\rm gas}$, while gas progressively delays the bar formation in the warm disks . The bar formation enhances a central mass concentration which in turn makes the bar decay temporarily, after which it regrows in size and strength, eventually becoming stronger in models with smaller $f_{\rm gas}$. Although all bars rotate fast in the beginning, they rapidly turn to slow rotators. In our models, only the gas-free, warm disk undergoes rapid buckling instability, while other disks thicken more gradually via vertical heating. The gas driven inward by the bar potential readily forms a star-forming nuclear ring. The ring is very small when it first forms and grows in size over time. The ring star formation rate is episodic and bursty due to feedback, and well correlated with the mass inflow rate to the ring. Some expanding shells produced by star formation feedback are sheared out in the bar regions and collide with dust lanes to appear as filamentary interbar spurs. The bars and nuclear rings formed in our simulations have properties similar to those in the Milky Way.

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

We aim to find missing microlensing planets hidden in the unanalyzed lensing events of previous survey data. For this purpose, we conduct a systematic inspection of high-magnification microlensing events, with peak magnifications $A_{\rm peak}\gtrsim 30$, in the data collected from high-cadence surveys in and before the 2018 season. From this investigation, we identify an anomaly in the lensing light curve of the event KMT-2018-BLG-1025. The analysis of the light curve indicates that the anomaly is caused by a very low mass-ratio companion to the lens. We identify three degenerate solutions, in which the ambiguity between a pair of solutions (solutions B) is caused by the previously known close--wide degeneracy, and the degeneracy between these and the other solution (solution A) is a new type that has not been reported before. The estimated mass ratio between the planet and host is $q\sim 0.8\times 10^{-4}$ for the solution A and $q\sim 1.6\times 10^{-4}$ for the solutions B. From the Bayesian analysis conducted with measured observables, we estimate that the masses of the planet and host and the distance to the lens are $(M_{\rm p}, M_{\rm h}, D_{\rm L})\sim (6.1~M_\oplus, 0.22~M_\odot, 6.7~{\rm kpc})$ for the solution A and $\sim (4.4~M_\oplus, 0.08~M_\odot, 7.5~{\rm kpc})$ for the solutions B. The planet mass is in the category of a super-Earth regardless of the solutions, making the planet the eleventh super-Earth planet, with masses lying between those of Earth and the Solar system's ice giants, discovered by microlensing., 10 pages, 10 figures, 5 tables

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 $\Delta\chi^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 $\Delta\chi^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., Comment: 10 pages, 9 figures

Utilizing James Clark Maxwell Telescope (JCMT) 850 $\mu$m SCUPOL dust polarization data, we investigate the configuration of the magnetic ($B$) field in the circumnuclear disk (CND) of the Galactic Center (GC). The SCUPOL data show a highly improved polarization coverage and resolution compared to earlier 100 $\mu$m observations. The 850 $\mu$m data have a resolution and coverage similar to previous 350 $\mu$m polarimetry data. However, with a proper sampling on a 10$"$ grid, we find the 850 $\mu$m data trace the morphological structures of the CND substantially better. Furthermore, because the 850 $\mu$m trace the field deeper into the material near Sgr A*, they represent the highest resolution submillimeter probe to date of the CND magnetic field. The observed $B$-field morphology is well described by a self-similar axisymmetric disk model where the radial infall velocity is one quarter of the rotational velocity. A detailed comparison with higher-resolution interferometric maps from the Submillimeter Array further reveals that the $B$-field aligns with the neutral gas streamers connecting to the CND. Moreover, the innermost observed $B$-field structure also appears to trace and align with the mini-spiral located inside the CND. This suggests that there is one underlying $B$-field structure that is connecting the CND with its streamers and the inner mini-spiral. An estimate of $\beta_{\rm Plasma} \lesssim 1-$based on the global B-field morphology that constrains the azimuthal-to-vertical field strength ratio of around 40 combined with a measurement of the azimuthal velocity--indicates that the B-field appears dynamically significant towards the CND and also onwards to the inner mini-spiral., Comment: 24 pages, 17 figure, accepted for publication in ApJ

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

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

We present a new face-on map of dense neutral atomic hydrogen (HI) gas in the outer Galaxy. Our map has been produced from the Leiden/Argentine/Bonn (LAB) HI 21-cm line all-sky survey by finding intensity maxima along every line of sight and then by projecting them on the Galactic plane. The resulting face-on map strikingly reveals the complex spiral structure beyond the solar circle, which is characterized by a mixture of distinct long arcs of HI concentrations and numerous 'interarm' features. The comparison with more conventional spiral tracers confirms the nature of those long arc structures as spiral arms. Our map shows that the HI spiral structure in the outer Galaxy is well described by a four-arm spiral model (pitch angle of 12 deg) with some deviations, and gives a new insight into identifying HI features associated with individual arms., 24 pages, 10 figures, 3 tables; accepted for publication in PASP

A fraction of early-type dwarf galaxies in the Virgo cluster have a disk component and even possess disk features such as bar, lens, and spiral arms. In this study, we construct 15 galaxy models that resemble VCC856, considered to be an infalling progenitor of disk dwarf galaxies, within observational error ranges, and use $N$-body simulations to study their long-term dynamical evolution in isolation as well as the formation of bar in them. We find that dwarf disk galaxies readily form bars unless they have an excessively concentrated halo or a hot disk. This suggests that infalling dwarf disk galaxies are intrinsically unstable to bar formation, even without any external perturbation, accounting for a population of barred dwarf galaxies in the outskirts of the Virgo cluster. The bars form earlier and stronger in galaxies with a lower fraction of counter-streaming motions, lower halo concentration, lower velocity anisotropy, and thinner disk. Similarly to normal disk galaxies, dwarf disk galaxies also undergo recurrent buckling instabilities. The first buckling instability tends to shorten the bar and to thicken the disk, and drives a dynamical transition in the bar pattern speed as well as mass inflow rate. In nine models, the bars regrow after the mild first buckling instability due to the efficient transfer of disk angular momentum to the halo, and are subject to recurrent buckling instabilities to turn into X-shaped bulges.

The interaction between a supermassive black hole (SMBH) and the surrounding material is of primary importance in modern astrophysics. The detection of the molecular 2-pc circumnuclear disk (CND) immediately around the Milky Way SMBH, SgrA*, provides an unique opportunity to study SMBH accretion at sub-parsec scales. Our new wide-field CS(2-1) map toward the Galactic center (GC) reveals multiple dense molecular streamers originated from the ambient clouds 20-pc further out, and connecting to the central 2 parsecs of the CND. These dense gas streamers appear to carry gas directly toward the nuclear region and might be captured by the central potential. Our phase-plot analysis indicates that these streamers show a signature of rotation and inward radial motion with progressively higher velocities as the gas approaches the CND and finally ends up co-rotating with the CND. Our results might suggest a possible mechanism of gas feeding the CND from 20 pc around 2 pc in the GC. In this paper, we discuss the morphology and the kinematics of these streamers. As the nearest observable Galactic nucleus, this feeding process may have implications for understanding the processes in extragalactic nuclei., Comment: 32 pages, 13 figures. Accepted for publication in ApJ

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., 17 figures and 2 tables; Accepted for publication in the ApJ

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., Accepted for publication in the ApJ

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., 22 pages, 17 figures; accepted for publication in the ApJ

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 $\beta$ 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 $\beta=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 $\beta\lesssim 10$ or $\mathcal{F}=10\%$ and $\beta\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., Comment: 41 pages, 12 figures, 3 tables, Accepted for publication in ApJ

Nuclear spirals are ubiquitous in galaxy centers. They exist not only in strong barred galaxies but also in galaxies without noticeable bars. We use high-resolution hydrodynamic simulations to study the properties of nuclear gas spirals driven by weak bar-like and oval potentials. The amplitude of the spirals increases toward the center by a geometric effect, readily developing into shocks at small radii even for very weak potentials. The shape of the spirals and shocks depends rather sensitively on the background shear. When shear is low, the nuclear spirals are loosely wound and the shocks are almost straight, resulting in large mass inflows toward the center. When shear is high, on the other hand, the spirals are tightly wound and the shocks are oblique, forming a circumnuclear disk through which gas flows inward at a relatively lower rate. The induced mass inflow rates are enough to power black hole accretion in various types of Seyfert galaxies as well as to drive supersonic turbulence at small radii., 6 pages, 5 figures; accepted for the publication in the ApJL

We use three-dimensional numerical hydrodynamic simulations of the turbulent, multiphase atomic interstellar medium (ISM) to construct and analyze synthetic HI 21 cm emission and absorption lines. Our analysis provides detailed tests of 21 cm observables as physical diagnostics of the atomic ISM. In particular, we construct (1) the "observed" spin temperature, $T_{s,obs}(v_{ch})\equiv T_B(v_{ch})/[1-e^{-{\tau}(v_{ch})}]$, and its optical-depth weighted mean T_s,obs; (2) the absorption-corrected "observed" column density, $N_{H,obs}\propto \int dv_{ch} T_B(v_{ch}){\tau}(v_{ch})/[[1-e^{-{\tau}(v_{ch})}]$; and (3) the "observed" fraction of cold neutral medium (CNM), $f_{c,obs}\equiv T_c/T_{s,obs}$ for T_c the CNM temperature; we compare each observed parameter with true values obtained from line-of-sight (LOS) averages in the simulation. Within individual velocity channels, T_s,obs(v_ch) is within a factor 1.5 of the true value up to ${\tau}(v_{ch})\approx10$. As a consequence, N_H,obs and T_s,obs are respectively within 5% and 12% of the true values for 90% and 99% of LOSs. The optically thin approximation significantly underestimates N_H for ${\tau}>1$. Provided that T_c is constrained, an accurate observational estimate of the CNM mass fraction can be obtained down to 20%. We show that T_s,obs cannot be used to distinguish the relative proportions of warm and thermally-unstable atomic gas, although the presence of thermally-unstable gas can be discerned from 21 cm lines with 200K, Comment: 34 pages, 12 figures. Accepted for publication in ApJ. For Paper I, see http://arxiv.org/abs/1308.3231

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., 18 pages, 15 figures, 2 tables; Accepted for publication in the MNRAS

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., Comment: 17 pages, 11 figures, 2 tables; Accepted for publication in the ApJ

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., 17 pages, 13 figures; Accepted for publication in the ApJ

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

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., Comment: 22 pages, 19 figures, 3 tables; Accepted for publication in the ApJ; Version with full-resolution figures available at http://mirzam.snu.ac.kr/~wkim/Bar/mhdbar.pdf

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., 64 pages, 15 figures, accepted by the ApJ

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 $\delta$ and the nonlinear drag force F depend solely on \eta=A/(M^2-1) such that \delta/r_s=\eta and F/F_{lin}=(\eta/2)^{-0.45} for \eta>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., Comment: 40 pages, 18 figures, accepted for publication in Astrophysical Journal

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., Comment: ApJ accepted version; we expanded the discussion of the limitations of this work in Section 4.7, and included a new subsection (Section 4.8) to demonstrate the convergence test for the resolution effects; 15 pages; emulateapj format. A movie showing the gas evolution in the canonical model is available on the ApJ website and at http://hubble.shao.ac.cn/~shen/nuclear_rings/canonicalmodel2.gif

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 $\alpha$ by a factor of $1+1/(2\beta_1)$ compared to the unmagnetized case, with $\beta_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 $\alpha^{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/(\beta_1 - 1)$. When the front experiences an acceleration, the IF instability cooperates with the Rayleigh-Taylor instability to make the front more unstable., Comment: 35 pages, 6 figures; Accepted for publication in the ApJ

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., 2 pages, 1 figure. To appear in the proceeding of the IAU GA XXVIII, Special Session 12: Modern Views of the Interstellar Medium

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., Comment: 10 pages, 1 table, 6 figures, Accepted for publication in the ApJ

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

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., 18 pages, 13 figures; Accepted for publication in the ApJ

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 n0.2 and self-gravity is included., 24 pages, 17 figures, 5 tables; Accepted for publication in the ApJ; Version with full-resolution figures available at http://mirzam.snu.ac.kr/~wkim/Bar/barHDn.pdf

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., Comment: Accepted for publication in the ApJ; Version with full-resolution figures available at http://mirzam.snu.ac.kr/~wkim/Bar/bar.pdf