Dynamically Driven Inflow onto the Galactic Center and its Effect upon Molecular Clouds

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 the 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.

Simulations of the star-forming molecular gas in an interacting M51-like galaxy: cloud population statistics

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 interacting 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.

Simulations of the Milky Ways central molecular zone -- II. Star formation

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 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 $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.

Simulations of the Milky Ways central molecular zone -- I. Gas dynamics

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 interstellar 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.

Simulations of the star-forming molecular gas in an interacting M51-like galaxy

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.

The geometry of the gas surrounding the Central Molecular Zone: on the origin of localised molecular clouds with extreme velocity dispersions

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.

A theoretical explanation for the Central Molecular Zone asymmetry

It has been known for more than thirty years that the distribution of molecular gas in the innermost 300 parsecs of the Milky Way, the Central Molecular Zone, is strongly asymmetric. Indeed, approximately three quarters of molecular emission comes from positive longitudes, and only one quarter from negative longitudes. However, despite much theoretical effort, the origin of this asymmetry has remained a mystery. Here we show that the asymmetry can be neatly explained by unsteady flow of gas in a barred potential. We use high-resolution 3D hydrodynamical simulations coupled to a state-of-the-art chemical network. Despite the initial conditions and the bar potential being point-symmetric with respect to the Galactic Centre, asymmetries develop spontaneously due to the combination of a hydrodynamical instability known as the wiggle instability and the thermal instability. The observed asymmetry must be transient: observations made tens of megayears in the past or in the future would often show an asymmetry in the opposite sense. Fluctuations of amplitude comparable to the observed asymmetry occur for a large fraction of the time in our simulations, and suggest that the present is not an exceptional moment in the life of our Galaxy.