The Galactic bar plays a critical role in the evolution of the Milky Ways Central Molecular Zone (CMZ), as its potential drives mass inward toward the Galactic Center via gas flows known as dust lanes. To explore the interaction between the CMZ and t
he dust lanes, we run hydrodynamic simulations in Arepo, modeling the potential of the Milky Ways bar in the absence of gas self-gravity and star formation physics, and we study the flows of mass using Monte Carlo tracer particles. We estimate the efficiency of the inflow via the dust lanes, finding that only about a third (30 +/- 12%) of the dust lanes mass initially accretes onto the CMZ, while the rest overshoots and accretes later. Given observational estimates of the amount of gas within the Milky Ways dust lanes, this suggests that the true total inflow rate onto the CMZ is 0.8 +/- 0.6 Msun yr$^{-1}$. Clouds in this simulated CMZ have sudden peaks in their average density near apocenter, where they undergo violent collisions with inflowing material. While these clouds tend to counter-rotate due to shear, co-rotating clouds occasionally occur (~52% are strongly counter-rotating, and ~7% are strongly co-rotating of the 44 cloud sample), likely due to the injection of momentum from collisions with inflowing material. We investigate the formation and evolution of these clouds, finding that they are fed by many discrete inflow events, providing a consistent source of gas to CMZ clouds even as they collapse and form stars.
To investigate how molecular clouds react to different environmental conditions at a galactic scale, we present a catalogue of giant molecular clouds resolved down to masses of $sim 10$~M$_{odot}$ from a simulation of the entire disc of an interactin
g M51-like galaxy and a comparable isolated galaxy. Our model includes time-dependent gas chemistry, sink particles for star formation and supernova feedback, meaning we are not reliant on star formation recipes based on threshold densities and can follow the physics of the cold molecular phase. We extract giant molecular clouds at a given timestep of the simulations and analyse their properties. In the disc of our simulated galaxies, spiral arms seem to act merely as snowplows, gathering gas and clouds without dramatically affecting their properties. In the centre of the galaxy, on the other hand, environmental conditions lead to larger, more massive clouds. While the galaxy interaction has little effect on cloud masses and sizes, it does promote the formation of counter-rotating clouds. We find that the identified clouds seem to be largely gravitationally unbound at first glance, but a closer analysis of the hierarchical structure of the molecular interstellar medium shows that there is a large range of virial parameters with a smooth transition from unbound to mostly bound for the densest structures. The common observation that clouds appear to be virialised entities may therefore be due to CO bright emission highlighting a specific level in this hierarchical binding sequence. The small fraction of gravitationally bound structures found suggests that low galactic star formation efficiencies may be set by the process of cloud formation and initial collapse.
Sgr E is a massive star formation complex found toward the Galactic center that consists of numerous discrete, compact HII regions. It is located at the intersection between the Central Molecular Zone (CMZ) and the far dust lane of the Galactic bar,
similar to hot spots seen in external galaxies. Compared with other Galactic star formation complexes, the Sgr E complex is unusual because its HII regions all have similar radio luminosities and angular extents, and they are deficient in ~10micron emission from their photodissociation regions (PDRs). Our Green Bank Telescope (GBT) radio recombination line observations increase the known membership of Sgr E to 19 HII regions. There are 43 additional HII region candidates in the direction of Sgr E, 26 of which are detected for the first time here using MeerKAT 1.28GHz data. Therefore, the true HII region population of Sgr E may number >60. Using APEX SEDIGISM 13CO 2-1 data we discover a 3.0x10^5 Solar mass molecular cloud associated with Sgr E, but find few molecular or far-infrared concentrations at the locations of the Sgr E HII regions. Comparison with simulations and an analysis of its radio continuum properties indicate that Sgr E formed upstream in the far dust lane of the Galactic bar a few Myr ago and will overshoot the CMZ, crashing into the near dust lane. We propose that the unusual infrared properties of the Sgr E HII regions are caused by their orbit about the Galactic center, which have possibly stripped their PDRs.
The nuclear stellar disc (NSD) is a flattened stellar structure that dominates the gravitational potential of the Milky Way at Galactocentric radii $30 lesssim R lesssim 300{, rm pc}$. In this paper, we construct axisymmetric Jeans dynamical models o
f the NSD based on previous photometric studies and we fit them to line-of-sight kinematic data of APOGEE and SiO maser stars. We find that (i) the NSD mass is lower but consistent with the mass independently determined from photometry by Launhardt et al. (2002). Our fiducial model has a mass contained within spherical radius $r=100{, rm pc}$ of $M(r<100{, rm pc}) = 3.9 pm 1 times 10^8 {, rm M_odot}$ and a total mass of $M_{rm NSD} = 6.9 pm 2 times 10^8 {, rm M_odot}$. (ii) The NSD might be the first example of a vertically biased disc, i.e. with ratio between the vertical and radial velocity dispersion $sigma_z/sigma_R>1$. Observations and theoretical models of the star-forming molecular gas in the central molecular zone suggest that large vertical oscillations may be already imprinted at stellar birth. However, the finding $sigma_z/sigma_R > 1$ depends on a drop in the velocity dispersion in the innermost few tens of parsecs, on our assumption that the NSD is axisymmetric, and that the available (extinction corrected) stellar samples broadly trace the underlying light and mass distributions, all of which need to be established by future observations and/or modelling. (iii) We provide the most accurate rotation curve to date for the innermost $500 {, rm pc}$ of our Galaxy.
The Milky Ways 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 res
olve 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 $Rgtrsim100 {, rm pc}$, a significant amount also occurs closer to SgrA* at $R lesssim 10{, rm 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 ($sim100$ Myr), the depletion time of the CMZ is constant within a factor of $sim2$. 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) We study the trajectories of newly born stars in our simulations. We find several examples that have age and 3D velocity compatible with those of the Arches and Quintuplet clusters. Our simulations suggest that these prominent clusters originated near the collision sites where the bar-driven inflow accretes onto the CMZ, at symmetrical locations with respect to the Galactic centre, and that they have already decoupled from the gas in which they were born.
We use hydrodynamical simulations to study the Milky Ways 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 interstel
lar medium at sub-parsec resolution while also capturing the interaction between the CMZ and the bar-driven large-scale flow out to $Rsim 5kpc$. Our main findings are as follows: (1) The distinction between inner ($Rlesssim120$~pc) and outer ($120lesssim Rlesssim450$~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 $Rsimeq 200$~pc which includes the 1.3$^circ$ 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 $infty$-shaped structure identified in Herschel data by Molinari et al. 2011. (3) The bar in our simulation efficiently drives an inflow from the Galactic disc ($Rsimeq 3$~kpc) down to the CMZ ($Rsimeq200$~pc) of the order of $1rm,M_odot,yr^{-1}$, consistent with observational determinations. (4) Supernova feedback can drive an inflow from the CMZ inwards towards the circumnuclear disc of the order of $sim0.03,rm M_odot,yr^{-1}$. (5) We give a new interpretation for the 3D placement of the 20 and 50 km s$^{-1}$ clouds, according to which they are close ($Rlesssim30$~pc) to the Galactic centre, but are also connected to the larger-scale streams at $Rgtrsim100$~pc.
We present here the first of a series of papers aimed at better understanding the evolution and properties of giant molecular clouds (GMCs) in a galactic context. We perform high resolution, three-dimensional {sc arepo} simulations of an interacting
galaxy inspired by the well-observed M51 galaxy. Our fiducial simulations include a non-equilibrium, time-dependent, chemical network that follows the evolution of atomic and molecular hydrogen as well as carbon and oxygen self-consistently. Our calculations also treat gas self-gravity and subsequent star formation (described by sink particles), and coupled supernova feedback. In the densest parts of the simulated interstellar medium (ISM) we reach sub-parsec resolution, granting us the ability to resolve individual GMCs and their formation and destruction self-consistently throughout the galaxy. In this initial work we focus on the general properties of the ISM with a particular focus on the cold star-forming gas. We discuss the role of the interaction with the companion galaxy in generating cold molecular gas and controlling stellar birth. We find that while the interaction drives large-scale gas flows and induces spiral arms in the galaxy, it is of secondary importance in determining gas fractions in the different ISM phases and the overall star-formation rate. The behaviour of the gas on small GMC scales instead is mostly controlled by the self-regulating property of the ISM driven by coupled feedback.
Observations of molecular gas near the Galactic centre ($| l | < 10^circ$, $| b | < 1^circ$) reveal the presence of a distinct population of enigmatic compact clouds which are characterised by extreme velocity dispersions ($Delta v > 100, rm km/s$).
These Extended Velocity Features (EVFs) are very prominent in the datacubes and dominate the kinematics of molecular gas just outside the Central Molecular Zone (CMZ). The prototypical example of such a cloud is Bania Clump 2. We show that similar features are naturally produced in simulations of gas flow in a realistic barred potential. We analyse the structure of the features obtained in the simulations and use this to interpret the observations. We find that the features arise from collisions between material that has been infalling rapidly along the dust lanes of the Milky Way bar and material that belongs to one of the following two categories: (i) material that has `overshot after falling down the dust lanes on the opposite side; (ii) material which is part of the CMZ. Both types of collisions involve gas with large differences in the line-of-sight velocities, which is what produces the observed extreme velocity dispersions. Examples of both categories can be identified in the observations. If our interpretation is correct, we are directly witnessing (a) collisions of clouds with relative speeds of $sim 200, rm km/s$ and (b) the process of accretion of fresh gas onto the CMZ.
Fitting equilibrium dynamical models to observational data is an essential step in understanding the structure of the gaseous hot haloes that surround our own and other galaxies. However, the two main categories of models that are used in the literat
ure are poorly suited for this task: (i) simple barotropic models are analytic and can therefore be adjusted to match the observations, but are clearly unrealistic because the rotational velocity $v_phi(R,z)$ does not depend on the distance $z$ from the galactic plane, while (ii) models obtained as a result of cosmological galaxy formation simulations are more realistic, but are impractical to fit to observations due to high computational cost. Here we bridge this gap by presenting a general method to construct axisymmetric baroclinic equilibrium models of rotating galactic coronae in arbitrary external potentials. We consider in particular a family of models whose equipressure surfaces in the $(R,z)$ plane are ellipses of varying axis ratio. These models are defined by two one-dimensional functions, the axial ratio of pressure $q_{rm axis}(z)$ and the value of the pressure $P_{rm axis}(z)$ along the galaxys symmetry axis. These models can have a rotation speed $v_phi(R,z)$ that realistically decreases as one moves away from the galactic plane, and can reproduce the angular momentum distribution found in cosmological simulations. The models are computationally cheap to construct and can thus be used in fitting algorithms. We provide a python code that given $q_{rm axis}(z)$, $P_{rm axis}(z)$ and $Phi(R,z)$ returns $rho(R,z)$, $T(R,z)$, $P(R,z)$, $v_phi(R,z)$. We show a few examples of these models using the Milky Way as a case study.
We develop a dynamical theory for the origin of nuclear rings in barred galaxies. In analogy with the standard theory of accretion discs, our theory is based on shear viscous forces among nested annuli of gas. However, the fact that gas follows non c
ircular orbits in an external barred potential has profound consequences: it creates a region of reverse shear in which it is energetically favourable to form a stable ring which does not spread despite dissipation. Our theory allows us to approximately predict the size of the ring given the underlying gravitational potential. The size of the ring is loosely related to the location of the Inner Lindblad Resonance in the epicyclic approximation, but the predicted location is more accurate and is also valid for strongly barred potentials. By comparing analytical predictions with the results of hydrodynamical simulations, we find that our theory provides a viable mechanism for ring formation if the effective sound speed of the gas is low ($cslesssim1kms$), but that nuclear spirals/shocks created by pressure destroy the ring when the sound speed is high ($cssimeq10kms$). We conclude that whether this mechanism for ring formation is relevant for real galaxies ultimately depends on the effective equation of state of the ISM. Promising confirmation comes from simulations in which the ISM is modelled using state-of-the-art cooling functions coupled to live chemical networks, but more tests are needed regarding the role of turbulence driven by stellar feedback. If the mechanism is relevant in real galaxies, it could provide a powerful tool to constrain the gravitational potential, in particular the bar pattern speed.
