Supernova Feedback Research Articles

The Single-cloud Star Formation Relation

Abstract One of the most important and well-established empirical results in astronomy is the Kennicutt–Schmidt relation between the density of interstellar gas and the rate at which that gas forms stars. A tight correlation between these quantities has long been measured at galactic scales. More recently, using surveys of YSOs, a KS relationship has been found within molecular clouds relating the surface density of star formation to the surface density of gas; however, the scaling of these laws varies significantly from cloud to cloud. In this Letter, we use a recently developed, high-accuracy catalog of young stellar objects from Spitzer combined with high-dynamic-range gas column density maps of 12 nearby (<1.5 kpc) molecular clouds from Herschel to re-examine the KS relation within individual molecular clouds. We find a tight, linear correlation between clouds’ star formation rate per unit area and their gas surface density normalized by the gas freefall time. The measured intracloud KS relation, which relates star formation rate to the volume density, extends over more than two orders of magnitude within each cloud and is nearly identical in each of the 12 clouds, implying a constant star formation efficiency per freefall time ϵ ff ≈ 0.026. The finding of a universal correlation within individual molecular clouds, including clouds that contain no massive stars or massive stellar feedback, favors models in which star formation is regulated by local processes such as turbulence or stellar feedback such as protostellar outflows, and disfavors models in which star formation is regulated only by galaxy properties or supernova feedback on galactic scales.

A little FABLE: exploring AGN feedback in dwarf galaxies with cosmological simulations

ABSTRACT Contrary to the standard lore, there is mounting observational evidence that feedback from active galactic nuclei (AGN) may also play a role at the low-mass end of the galaxy population. We investigate this using the cosmological simulation suite fable, with a particular focus on the dwarf regime (Mstellar < 109.5 M⊙). We find that overmassive black holes (BHs), with respect to the mean scaling relations with their host galaxies, drive hotter and faster outflows and lead to significantly reduced gas mass fractions. They are also more likely to display a kinematically misaligned ionized gas component in our mock MaNGA velocity maps, although we caution that cosmic inflows and mergers contribute to misalignments as well. While in the local Universe the majority of AGN in dwarfs are much dimmer than the stellar component, for z ≥ 2 there is a significant population that outshines their hosts. These high-redshift overmassive BHs contribute to the quenching of dwarfs, whereas at late cosmic times supernova (SN) feedback is more efficient. While our results are overall in good agreement with X-ray observations of AGN in dwarfs, the lack of high-luminosity X-ray AGN in fable at low redshifts highlights an interesting possibility that SN feedback could be too strong in fable’s dwarfs, curtailing AGN growth and feedback. We predict that future observations may uncover many more AGN in dwarfs with lower luminosities and at higher redshifts.

Astraeus I: the interplay between galaxy formation and reionization

ABSTRACT We introduce a new self-consistent model of galaxy evolution and reionization, astraeus (seminumerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dArk mattEr simUlationS), which couples a state-of-the-art N-body simulation with the semi-analytical galaxy evolution delphi and the seminumerical reionization scheme cifog. astraeus includes all the key processes of galaxy formation and evolution (including accretion, mergers, supernova, and radiative feedback) and follows the time and spatial evolution of the ionized regions in the intergalactic medium (IGM). Importantly, it explores different radiative feedback models that cover the physically plausible parameter space, ranging from a weak and delayed to a strong and immediate reduction of gas mass available for star formation. From our simulation suite that covers the different radiative feedback prescriptions and ionization topologies, we find that radiative feedback continuously reduces star formation in galaxies with $M_\mathrm{ h}\lesssim 10^{9.5}\, {\rm \rm M_\odot }$ upon local reionization; larger mass haloes are unaffected even for the strongest and immediate radiative feedback cases during reionization. For this reason, the ionization topologies of different radiative feedback scenarios differ only on scales smaller than 1–2 comoving Mpc, and significant deviations are found only when physical parameters (e.g. the escape fraction of ionizing photons) are altered based on galactic properties. Finally, we find that observables (the ultraviolet luminosity function, stellar mass function, reionization histories and ionization topologies) are hardly affected by the choice of the used stellar population synthesis models that model either single stars or binaries.

The Star Formation History of Eridanus II: On the Role of Supernova Feedback in the Quenching of Ultrafaint Dwarf Galaxies*

Abstract Eridanus II (Eri II) is an ultrafaint dwarf (UFD) galaxy (M V = −7.1) located at a distance close to the Milky Way virial radius. Early shallow color–magnitude diagrams (CMDs) indicated that it possibly hosted an intermediate-age or even young stellar population, which is unusual for a galaxy of this mass. In this paper, we present new Hubble Space Telescope/Advanced Camera for Surveys CMDs reaching the oldest main-sequence turnoff with excellent photometric precision and derive a precise star formation history (SFH) for this galaxy through CMD fitting. This SFH shows that the bulk of the stellar mass in Eri II formed in an extremely short star formation burst at the earliest possible time. The derived star formation rate profile has a width at half maximum of 500 Myr and reaches a value compatible with null star formation 13 Gyr ago. However, tests with mock stellar populations and with the CMD of the globular cluster M92 indicate that the star formation period could be shorter than 100 Myr. From the quantitative determination of the amount of mass turned into stars in this early star formation burst ( ∼2 × 105 M ⊙) we infer the number of supernova (SN) events and the corresponding energy injected into the interstellar medium. For reasonable estimates of the Eri II virial mass and values of the coupling efficiency of the SN energy, we conclude that Eri II could be quenched by SN feedback alone, thus casting doubts on the need to invoke cosmic reionization as the preferred explanation for the early quenching of old UFD galaxies.

Radiation Hydrodynamics of Turbulent H ii Regions in Molecular Clouds: A Physical Origin of LyC Leakage and the Associated Lyα Spectra

Abstract We examine Lyman continuum (LyC) leakage through H ii regions regulated by turbulence and radiative feedback in a giant molecular cloud in the context of fully coupled radiation hydrodynamics (RHD). The physical relations of the LyC escape with H i covering fraction, kinematics, ionizing photon production efficiency, and emergent Lyα line profiles are studied using a series of RHD turbulence simulations performed with ramses-rt. The turbulence-regulated mechanism allows ionizing photons to leak out at early times before the onset of supernova feedback. The LyC photons escape through turbulence-generated low column density channels that are evacuated efficiently by radiative feedback via photoheating-induced shocks across the D-type ionization fronts. The Lyα photons funnel through the photoionized channels along the paths of LyC escape, resulting in a diverse Lyα spectral morphology including narrow double-peaked profiles. The Lyα peak separation is controlled by the residual H i column density of the channels, and the line asymmetry correlates with the porosity and multiphase structure of the H ii region. This mechanism through the turbulent H ii regions can naturally reproduce the observed Lyα spectral characteristics of some of the LyC-leaking galaxies. This RHD turbulence origin provides an appealing hypothesis to explain high LyC leakage from very young (∼3 Myr) star-forming galaxies found in the local universe without need of extreme galactic outflows or supernova feedback. We discuss the implications of the turbulent H ii regions on other nebular emission lines and a possible observational test with the Magellanic System and local blue compact dwarf galaxies as analogs of reionization-era systems.

An extended halo around an ancient dwarf galaxy

The Milky Way is surrounded by dozens of ultra-faint (< $10^5$ solar luminosities) dwarf satellite galaxies. They are the surviving remnants of the earliest galaxies, as confirmed by their ancient (~13 billion years old) and chemically primitive stars. Simulations suggest that these systems formed within extended dark matter halos and experienced early galaxy mergers and supernova feedback. However, the signatures of these events would lie outside their core regions (>2 half-light radii), which are spectroscopically unstudied due to the sparseness of their distant stars. Here we identify members of the Tucana II ultra-faint dwarf galaxy in its outer region (up to 9 half-light radii), demonstrating the system to be dramatically more spatially extended and chemically primitive than previously found. These distant stars are extremely metal-poor (<[Fe/H]>=-3.02; less than ~1/1000th of the solar iron abundance), affirming Tucana II as the most metal-poor known galaxy. We observationally establish, for the first time, an extended dark matter halo surrounding an ultra-faint dwarf galaxy out to one kiloparsec, with a total mass of >$10^7$ solar masses. This measurement is consistent with the expected ~2x$10^7$ solar masses using a generalized NFW density profile. The extended nature of Tucana II suggests that it may have undergone strong bursty feedback or been the product of an early galactic merger. We demonstrate that spatially extended stellar populations, which other ultra-faint dwarfs hint at hosting as well, are observable in principle and open the possibility for detailed studies of the stellar halos of relic galaxies.

The Horizon Run 5 Cosmological Hydrodynamical Simulation: Probing Galaxy Formation from Kilo- to Gigaparsec Scales

Abstract Horizon Run 5 (HR5) is a cosmological hydrodynamical simulation that captures the properties of the universe on a Gpc scale while achieving a resolution of 1 kpc. Inside the simulation box, we zoom in on a high-resolution cuboid region with a volume of 1049 × 119 × 127 cMpc3. The subgrid physics chosen to model galaxy formation includes radiative heating/cooling, UV background, star formation, supernova feedback, chemical evolution tracking the enrichment of oxygen and iron, the growth of supermassive black holes, and feedback from active galactic nuclei in the form of a dual jet-heating mode. For this simulation, we implemented a hybrid MPI-OpenMP version of RAMSES, specifically targeted for modern many-core many-thread parallel architectures. In addition to the traditional simulation snapshots, lightcone data were generated on the fly. For the post-processing, we extended the friends-of-friend algorithm and developed a new galaxy finder PGalF to analyze the outputs of HR5. The simulation successfully reproduces observations, such as the cosmic star formation history and connectivity of galaxy distribution, We identify cosmological structures at a wide range of scales, from filaments with a length of several cMpc, to voids with a radius of ∼ 100 cMpc. The simulation also indicates that hydrodynamical effects on small scales impact galaxy clustering up to very large scales near and beyond the baryonic acoustic oscillation scale. Hence, caution should be taken when using that scale as a cosmic standard ruler: one needs to carefully understand the corresponding biases. The simulation is expected to be an invaluable asset for the interpretation of upcoming deep surveys of the universe.

Hydrogen–helium chemical and nuclear galaxy collision: Hydrodynamic simulations on AVX-512 supercomputers

A computational hydrodynamic model of interacting galaxies is presented. The interstellar medium (ISM) is described by a model of gravitational multicomponent single-velocity hydrodynamics. A model based on first moments of the collisionless Boltzmann equation is used to describe the stellar component and dark matter. Subgrid processes of star formation and supernovae feedback, as well as cooling and heating functions are added to the hydrodynamic model. The computational model includes chemical and nuclear reactions of the basic forms of hydrogen and helium before the formation of the helium hydride ion. The hydrodynamic equations have a suitable form for solving by a unified computational method based on the “Harten – Lax – van Leer” (HLL) method. The equations and the numerical method are written in vector form to use advanced vector extensions AVX-512 to speed up the calculations. A parallel implementation of the model and some computational simulation experiments on the chemical dynamics of interacting galaxies are presented.

The origin of low-surface-brightness galaxies in the dwarf regime

ABSTRACT Low-surface-brightness galaxies (LSBGs) – defined as systems that are fainter than the surface-brightness limits of past wide-area surveys – form the overwhelming majority of galaxies in the dwarf regime (M⋆ &amp;lt; 109 M⊙). Using NewHorizon, a high-resolution cosmological simulation, we study the origin of LSBGs and explain why LSBGs at similar stellar mass show the large observed spread in surface brightness. NewHorizon galaxies populate a well-defined locus in the surface brightness–stellar mass plane, with a spread of ∼3 mag arcsec−2, in agreement with deep Sloan Digital Sky Survey (SDSS) Stripe 82 data. Galaxies with fainter surface brightnesses today are born in regions of higher dark matter density. This results in faster gas accretion and more intense star formation at early epochs. The stronger resultant supernova feedback flattens gas profiles at a faster rate, which, in turn, creates shallower stellar profiles (i.e. more diffuse systems) more rapidly. As star formation declines towards late epochs ( z &amp;lt; 1), the larger tidal perturbations and ram pressure experienced by these systems (due to their denser local environments) accelerate the divergence in surface brightness, by increasing their effective radii and reducing star formation, respectively. A small minority of dwarfs depart from the main locus towards high surface brightnesses, making them detectable in past wide surveys (e.g. standard-depth SDSS images). These systems have anomalously high star formation rates, triggered by recent fly-by or merger-driven starbursts. We note that objects considered extreme or anomalous at the depth of current data sets, e.g. ‘ultra-diffuse galaxies’, actually dominate the predicted dwarf population and will be routinely visible in future surveys like the Legacy Survey of Space and Time (LSST).

Cosmological direct-collapse black hole formation sites hostile for their growth

ABSTRACT The direct collapse (DC) is a promising mechanism that provides massive seed black holes (BHs) with ∼105 M⊙ in the early universe. To study a long-term accretion growth of a direct-collapse black hole (DCBH), we perform cosmological radiation-hydrodynamics simulations, extending our previous work where we investigated its formation stage. With a high spatial resolution down below the Bondi radius, we show that the accretion rate on to the BH is far below the Eddington value. Such slow mass growth is partly because of the strong radiative feedback from the accreting BH to the surrounding dense gas. Even after it falls into the first galaxy, the accretion rate is substantially suppressed due to the supernova feedback associated with the intense star formation. Moreover, the BH has a large velocity of ∼100 km s−1 relative to the gas, which further reduces the accretion rate. This large relative velocity stems from the fact that the DCBHs form in metal-free environments typically at ∼1 kpc from the galaxy. The BH accelerates as it approaches the galactic centre due to the gravity. The relative velocity never damps and the BH wanders around the outer galactic region. An analytic estimate predicts that the DCBH formation within ∼100 pc around the galactic centre is necessary to decelerate the BH with dynamical friction before z = 7. Since metal enrichment with Z ∼ 10−5−10−3 Z⊙ is expected there, the formation of DCBHs in the metal-enriched environments is preferable for the subsequent rapid growth.

The nature of giant clumps in high-z discs: a deep-learning comparison of simulations and observations

ABSTRACT We use deep learning to explore the nature of observed giant clumps in high-redshift disc galaxies, based on their identification and classification in cosmological simulations. Simulated clumps are detected using the 3D gas and stellar densities in the VELA zoom-in cosmological simulation suite, with ${\sim}25\ \rm {pc}$ maximum resolution, targeting main-sequence galaxies at 1 &amp;lt; z &amp;lt; 3. The clumps are classified as long-lived clumps (LLCs) or short-lived clumps (SLCs) based on their longevity in the simulations. We then train neural networks to detect and classify the simulated clumps in mock, multicolour, dusty, and noisy HST-like images. The clumps are detected using an encoder–decoder convolutional neural network (CNN), and are classified according to their longevity using a vanilla CNN. Tests using the simulations show our detector and classifier to be ${\sim}80{{\ \rm per\ cent}}$ complete and ${\sim}80{{\ \rm per\ cent}}$ pure for clumps more massive than ∼107.5 M⊙. When applied to observed galaxies in the CANDELS/GOODS S+N fields, we find both types of clumps to appear in similar abundances in the simulations and the observations. LLCs are, on average, more massive than SLCs by ∼0.5 dex, and they dominate the clump population above Mc ≳ 107.6 M⊙. LLCs tend to be found closer to the galactic centre, indicating clump migration to the centre or preferential formation at smaller radii. The LLCs are found to reside in high-mass galaxies, indicating better clump survivability under supernova feedback there, due to clumps being more massive in these galaxies. We find the clump masses and radial positions in the simulations and the observations to agree within a factor of 2.

Most of the cool CGM of star-forming galaxies is not produced by supernova feedback

Abstract The characterization of the large amount of gas residing in the galaxy halos, the so called circumgalactic medium (CGM), is crucial to understand galaxy evolution across cosmic time. We focus here on the the cool (T ∼ 104 K) phase of this medium around star-forming galaxies in the local universe, whose properties and dynamics are poorly understood. We developed semi-analytical parametric models to describe the cool CGM as an outflow of gas clouds from the central galaxy, as a result of supernova explosions in the disc (galactic wind). The cloud motion is driven by the galaxy gravitational pull and by the interactions with the hot (T ∼ 106 K) coronal gas. Through a bayesian analysis, we compare the predictions of our models with the data of the COS-Halos and COS-GASS surveys, which provide accurate kinematic information of the cool CGM around more than 40 low-redshift star-forming galaxies, probing distances up to the galaxy virial radii. Our findings clearly show that a supernova-driven outflow model is not suitable to describe the dynamics of the cool circumgalactic gas. Indeed, to reproduce the data, we need extreme scenarios, with initial outflow velocities and mass loading factors that would lead to unphysically high energy coupling from the supernovae to the gas and with supernova efficiencies largely exceeding unity. This strongly suggests that, since the outflows cannot reproduce most of the cool gas absorbers, the latter are likely the result of cosmological inflow in the outer galaxy halos, in analogy to what we have previously found for early-type galaxies.

Radiative Supernova Remnants and Supernova Feedback

Abstract Supernova (SN) explosions are a major feedback mechanism, regulating star formation in galaxies through their momentum input. We review the observations of SNRs in radiative stages in the Milky Way, to validate theoretical results regarding the momentum/energy injection from a single SN explosion. For seven supernova remnants (SNRs) where we can observe fast-expanding, atomic radiative shells, we show that the shell momentum inferred from H i 21 cm line observations is in the range of (0.5–4.5) × 105 M ⊙ km s−1. In two SNRs (W44 and IC 443), shocked molecular gas with momentum comparable to that of the atomic SNR shells has also been observed. We compare the momentum and kinetic/thermal energy of these seven SNRs with the results from 1D and 3D numerical simulations. The observation-based momentum and kinetic energy agree well with the expected momentum/energy input from an SN explosion of ∼1051 erg. It is much more difficult to use data/model comparisons of thermal energy to constrain the initial explosion energy, however, due to rapid cooling and complex physics at the hot/cool interface in radiative SNRs. We discuss the observational and theoretical uncertainties of these global parameters and explosion energy estimates for SNRs in complex environments.

Galactic 26Al traces metal loss through hot chimneys

ABSTRACT Radioactive 26Al is an excellent tracer for metal ejection in the Milky Way, and can provide a direct constraint on the modelling of supernova feedback in galaxy evolution. Gamma-ray observations of the 26Al decay line have found high velocities and hence require a significant fraction of the Galactic 26Al in the hot component. At the same time, meteoritic data combined with simulation results suggest that a significant amount of 26Al makes its way into stars before decay. We investigated the distribution into hot and cold channels with a simulation of a Milky-Way-like galaxy with massive-star feedback in superbubbles and with ejecta traced by 26Al. About 30–40 per cent of the ejecta remain hot, with typical cooling times of the order Gyr. 26Al traces the footpoints of a chimney-fed outflow that mixes metals turbulently into the halo of the model galaxy on a scale of at least 50 kpc. The rest diffuses into cold gas ≲ 104 K, and may therefore be quickly available for star formation. We discuss the robustness of the result by comparison to a simulation with a different global flow pattern. The branching ratio into hot and cold components is comparable to that of longer term average results from chemical evolution modelling of galaxies, clusters, and the intracluster medium.

The Cloud Factory II: gravoturbulent kinematics of resolved molecular clouds in a galactic potential

ABSTRACTWe present a statistical analysis of the gravoturbulent velocity fluctuations in molecular cloud complexes extracted from our ‘Cloud Factory’ Galactic-scale interstellar medium (ISM) simulation suite. For this purpose, we produce non-local thermodynamic equilibrium 12CO J = 1 − 0 synthetic observations and apply the principal component analysis (PCA) reduction technique on a representative sample of cloud complexes. The velocity fluctuations are self-consistently generated by different physical mechanisms at play in our simulations, which include Galactic-scale forces, gas self-gravity, and supernova feedback. The statistical analysis suggests that, even though purely gravitational effects are necessary to reproduce standard observational laws, they are not sufficient in most cases. We show that the extra injection of energy from supernova explosions plays a key role in establishing the global turbulent field and the local dynamics and morphology of molecular clouds. Additionally, we characterize structure function scaling parameters as a result of cloud environmental conditions: some of the complexes are immersed in diffuse (interarm) or dense (spiral-arm) environments, and others are influenced by embedded or external supernovae. In quiescent regions, we obtain time-evolving trajectories of scaling parameters driven by gravitational collapse and supersonic turbulent flows. Our findings suggest that a PCA-based statistical study is a robust method to diagnose the physical mechanisms that drive the gravoturbulent properties of molecular clouds. Also, we present a new open source module, the pcafactory, which smartly performs PCA to extract velocity structure functions from simulated or real data of the ISM in a user-friendly way.

Massive Black Hole Merger Rates: The Effect of Kiloparsec Separation Wandering and Supernova Feedback

Abstract We revisit the predictions for the merger rate of massive black hole binaries detectable by the Laser Interferometer Space Antenna (LISA) and their background signal for pulsar-timing arrays. We focus on the effect of the delays between the merger of galaxies and the final coalescence of black hole binaries, and on supernova feedback on the black hole growth. By utilizing a semianalytic galaxy formation model, not only do we account for the driving of the evolution of binaries at separations ≲1 pc (gas-driven migration, stellar hardening, and triple/quadruple massive black hole systems), but we also improve on previous studies by accounting for the time spent by black hole pairs from kiloparsec down to parsec separation. We also include the effect of supernova feedback, which may eject gas from the nuclear region of low-mass galaxies, thus hampering the growth of black holes via accretion and suppressing their orbital migration in circumbinary disks. Despite including these novel physical effects, we predict that the LISA detection rate should still be , irrespective of the model for the black hole seeds at high redshifts. Scenarios where black holes form from seeds are more significantly impacted by supernova feedback. We also find that for detectable events, the merging black holes typically have mass ratios between ∼0.1 and 1. Predictions for the stochastic background in the band of pulsar-timing array experiments are instead rather robust and show only a mild dependence on the model.

The Effect of Supernovae on the Turbulence and Dispersal of Molecular Clouds

Abstract We study the impact of supernovae on individual molecular clouds, using a high-resolution magnetohydrodynamic simulation of a 250 pc region where we resolve the formation of individual massive stars. The supernova feedback is implemented with real supernovae, meaning supernovae that are the natural evolution of the resolved massive stars, so their position and timing are self-consistent. We select a large sample of molecular clouds from the simulation to investigate the supernova energy injection and the resulting properties of molecular clouds. We find that molecular clouds have a lifetime of a few dynamical times, less than half of them contract to the point of becoming gravitationally bound, and the dispersal time of bound clouds of order one dynamical time is a factor of 2 shorter than that of unbound clouds. We stress the importance of internal supernovae, that is, massive stars that explode inside their parent cloud, in setting the cloud dispersal time, and their huge overdensity compared to models where the supernovae are randomly distributed. We also quantify the energy injection efficiency of supernovae as a function of supernova distance to the clouds. We conclude that intermittent driving by supernovae can maintain molecular cloud turbulence and may be the main process for cloud dispersal and that the full role of supernovae in the evolution of molecular clouds cannot be fully accounted for without a self-consistent implementation of the supernova feedback.

Cosmological simulations of the same spiral galaxy: the impact of baryonic physics

ABSTRACT The interplay of star formation (SF) and supernova (SN) feedback in galaxy formation is a key element for understanding galaxy evolution. Since these processes occur at small scales, it is necessary to have sub-grid models that recover their evolution and environmental effects at the scales reached by cosmological simulations. In this work, we present the results of the Mochima simulation, where we simulate the same spiral galaxy inhabiting a Milky Way (MW) size halo in a cosmological environment changing the sub-grid models for SN feedback and SF. We test combinations of the Schmidt law and a multifreefall based SF with delayed cooling feedback or mechanical feedback. We reach a resolution of 35 pc in a zoom-in box of 36 Mpc. For this, we use the code $\rm{\small RAMSES}$ with the implementation of gas turbulence in time and trace the local hydrodynamical features of the star-forming gas. Finally, we compare the galaxies at redshift 0 with global and interstellar medium observations in the MW and local spiral galaxies. The simulations show successful comparisons with observations. Nevertheless, diverse galactic morphologies are obtained from different numerical implementations. We highlight the importance of detailed modelling of the SF and feedback processes, especially for simulations with a resolution that start to reach scales relevant for molecular cloud physics. Future improvements could alleviate the degeneracies exhibited in our simulated galaxies under different sub-grid models.

The effects of subgrid models on the properties of giant molecular clouds in galaxy formation simulations

ABSTRACT Recent cosmological hydrodynamical simulations are able to reproduce numerous statistical properties of galaxies that are consistent with observational data. Yet, the adopted subgrid models strongly affect the simulation outcomes, limiting the predictive power of these simulations. In this work, we perform a suite of isolated galactic disc simulations under the SMUGGLE framework and investigate how different subgrid models affect the properties of giant molecular clouds (GMCs). We employ astrodendro, a hierarchical clump-finding algorithm, to identify GMCs in the simulations. We find that different choices of subgrid star formation efficiency, ϵff, and stellar feedback channels, yield dramatically different mass and spatial distributions for the GMC populations. Without feedback, the mass function of GMCs has a shallower power-law slope and extends to higher mass ranges compared to runs with feedback. Moreover, higher ϵff results in faster molecular gas consumption and steeper mass function slopes. Feedback also suppresses power in the two-point correlation function (TPCF) of the spatial distribution of GMCs. Specifically, radiative feedback strongly reduces the TPCF on scales below 0.2 kpc, while supernova feedback reduces power on scales above 0.2 kpc. Finally, runs with higher ϵff exhibit a higher TPCF than runs with lower ϵff, because the dense gas is depleted more efficiently, thereby facilitating the formation of well-structured supernova bubbles. We argue that comparing simulated and observed GMC populations can help better constrain subgrid models in the next generation of galaxy formation simulations.

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 drive an inflow from the CMZ inwards towards the circumnuclear disc of the order of ${\sim}0.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 (R ≲ 30 pc) to the Galactic Centre, but are also connected to the larger scale streams at R ≳ 100 pc.