Simulations of the Milky Way’s central molecular zone – I. Gas dynamics
ABSTRACT We use hydrodynamical simulations to study the Milky Way’s central molecular zone (CMZ). The simulations include a non-equilibrium chemical network, the gas self-gravity, star formation, and supernova feedback. We resolve the structure of the interstellar medium at sub-parsec resolution while also capturing the interaction between the CMZ and the bar-driven large-scale flow out to $R\sim 5\, {\rm kpc}$. Our main findings are as follows: (1) The distinction between inner (R ≲ 120 pc) and outer (120 ≲ R ≲ 450 pc) CMZ that is sometimes proposed in the literature is unnecessary. Instead, the CMZ is best described as single structure, namely a star-forming ring with outer radius R ≃ 200 pc which includes the 1.3° complex and which is directly interacting with the dust lanes that mediate the bar-driven inflow. (2) This accretion can induce a significant tilt of the CMZ out of the plane. A tilted CMZ might provide an alternative explanation to the ∞-shaped structure identified in Herschel data by Molinari et al. (3) The bar in our simulation efficiently drives an inflow from the Galactic disc (R ≃ 3 kpc) down to the CMZ (R ≃ 200 pc) of the order of $1\rm \, M_\odot \, yr^{-1}$, consistent with observational determinations. (4) Supernova feedback can more »
Authors:
; ; ; ; ; ; ;
Award ID(s):
Publication Date:
NSF-PAR ID:
10283766
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
499
Issue:
3
Page Range or eLocation-ID:
4455 to 4478
ISSN:
0035-8711
Previous studies of fueling black holes in galactic nuclei have argued (on scales ${\sim}0.01{-}1000\,$pc) accretion is dynamical with inflow rates $\dot{M}\sim \eta \, M_{\rm gas}/t_{\rm dyn}$ in terms of gas mass Mgas, dynamical time tdyn, and some η. But these models generally neglected expulsion of gas by stellar feedback, or considered extremely high densities where expulsion is inefficient. Studies of star formation, however, have shown on sub-kpc scales the expulsion efficiency fwind = Mejected/Mtotal scales with the gravitational acceleration as $(1-f_{\rm wind})/f_{\rm wind}\sim \bar{a}_{\rm grav}/\langle \dot{p}/m_{\ast }\rangle \sim \Sigma _{\rm eff}/\Sigma _{\rm crit}$ where $\bar{a}_{\rm grav}\equiv G\, M_{\rm tot}(\lt r)/r^{2}$ and $\langle \dot{p}/m_{\ast }\rangle$ is the momentum injection rate from young stars. Adopting this as the simplest correction for stellar feedback, $\eta \rightarrow \eta \, (1-f_{\rm wind})$, we show this provides a more accurate description of simulations with stellar feedback at low densities. This has immediate consequences, predicting the slope and normalization of the MBH − σ and MBH − Mbulge relation, LAGN −SFR relations, and explanations for outliers in compact Es. Most strikingly, because star formation simulations show expulsion is efficient (fwind ∼ 1) below total-mass surface density $M_{\rm tot}/\pi \, r^{2}\lt \Sigma _{\rm crit}\sim 3\times 10^{9}\, \mathrm{M}_{\odotmore » 2. ABSTRACT We report on our combined analysis of HST, VLT/MUSE, VLT/SINFONI, and ALMA observations of the local Seyfert 2 galaxy, NGC 5728 to investigate in detail the feeding and feedback of the active galactic nucleus (AGN). The data sets simultaneously probe the morphology, excitation, and kinematics of the stars, ionized gas, and molecular gas over a large range of spatial scales (10 pc to 10 kpc). NGC 5728 contains a large stellar bar that is driving gas along prominent dust lanes to the inner 1 kpc where the gas settles into a circumnuclear ring. The ring is strongly star forming and contains a substantial population of young stars as indicated by the lowered stellar velocity dispersion and gas excitation consistent with H ii regions. We model the kinematics of the ring using the velocity field of the CO (2–1) emission and stars and find it is consistent with a rotating disc. The outer regions of the disc, where the dust lanes meet the ring, show signatures of inflow at a rate of 1 M$\odot$yr−1. Inside the ring, we observe three molecular gas components corresponding to the circular rotation of the outer ring, a warped disc, and the nuclear stellar bar. The AGN is driving an ionized gasmore » 3. Abstract Large-scale bars can fuel galaxy centers with molecular gas, often leading to the development of dense ringlike structures where intense star formation occurs, forming a very different environment compared to galactic disks. We pair ∼0.″3 (30 pc) resolution new JWST/MIRI imaging with archival ALMA CO(2–1) mapping of the central ∼5 kpc of the nearby barred spiral galaxy NGC 1365 to investigate the physical mechanisms responsible for this extreme star formation. The molecular gas morphology is resolved into two well-known bright bar lanes that surround a smooth dynamically cold gas disk (Rgal∼ 475 pc) reminiscent of non-star-forming disks in early-type galaxies and likely fed by gas inflow triggered by stellar feedback in the lanes. The lanes host a large number of JWST-identified massive young star clusters. We find some evidence for temporal star formation evolution along the ring. The complex kinematics in the gas lanes reveal strong streaming motions and may be consistent with convergence of gas streamlines expected there. Indeed, the extreme line widths are found to be the result of inter-“cloud” motion between gas peaks;ScousePydecomposition reveals multiple components with line widths of 〈σCO,scouse〉 ≈ 19 km s−1and surface densities of, similar to the properties observed throughout the rest of the central molecular gas structure. Tailored hydrodynamical simulations exhibit many of the observed properties and imply that the observed structures are transient and highly time-variable. From our study of NGC 1365, we conclude that it is predominantly the high gas inflow triggered by the bar that is setting the star formation in its CMZ. 4. Abstract The Milky Way’s central molecular zone (CMZ) has emerged in recent years as a unique laboratory for the study of star formation. Here we use the simulations presented in Tress et al. 2020 to investigate star formation in the CMZ. These simulations resolve the structure of the interstellar medium at sub-parsec resolution while also including the large-scale flow in which the CMZ is embedded. Our main findings are as follows. (1) While most of the star formation happens in the CMZ ring at R ≳ 100 pc, a significant amount also occurs closer to SgrA* at R ≲ 10 pc. (2) Most of the star formation in the CMZ happens downstream of the apocentres, consistent with the “pearls-on-a-string” scenario, and in contrast to the notion that an absolute evolutionary timeline of star formation is triggered by pericentre passage. (3) Within the timescale of our simulations (∼100 Myr), the depletion time of the CMZ is constant within a factor of ∼2. This suggests that variations in the star formation rate are primarily driven by variations in the mass of the CMZ, caused for example by AGN feedback or externally-induced changes in the bar-driven inflow rate, and not by variations in the depletion time. (4)more » 5. ABSTRACT It remains a major challenge to derive a theory of cloud-scale ($\lesssim100$pc) star formation and feedback, describing how galaxies convert gas into stars as a function of the galactic environment. Progress has been hampered by a lack of robust empirical constraints on the giant molecular cloud (GMC) lifecycle. We address this problem by systematically applying a new statistical method for measuring the evolutionary timeline of the GMC lifecycle, star formation, and feedback to a sample of nine nearby disc galaxies, observed as part of the PHANGS-ALMA survey. We measure the spatially resolved (∼100 pc) CO-to-H α flux ratio and find a universal de-correlation between molecular gas and young stars on GMC scales, allowing us to quantify the underlying evolutionary timeline. GMC lifetimes are short, typically$10\!-\!30\,{\rm Myr}$, and exhibit environmental variation, between and within galaxies. At kpc-scale molecular gas surface densities$\Sigma _{\rm H_2}\ge 8\,\rm {M_\odot}\,{{\rm pc}}^{-2}$, the GMC lifetime correlates with time-scales for galactic dynamical processes, whereas at$\Sigma _{\rm H_2}\le 8\,\rm {M_\odot}\,{{\rm pc}}^{-2}$GMCs decouple from galactic dynamics and live for an internal dynamical time-scale. After a long inert phase without massive star formation traced by H α (75–90 per cent of the cloud lifetime), GMCs disperse within just$1\!-\!5\,{\rm Myr}\$ oncemore »