-
Inflow and outflow properties, not total gas fractions, drive the evolution of the mass-metallicity relation
Authors:
Luigi Bassini,
Robert Feldmann,
Jindra Gensior,
Claude-André Faucher-Giguère,
Elia Cenci,
Jorge Moreno,
Mauro Bernardini,
Lichen Liang
Abstract:
Observations show a tight correlation between the stellar mass of galaxies and their gas-phase metallicity (MZR). This relation evolves with redshift, with higher-redshift galaxies being characterized by lower metallicities. Understanding the physical origin of the slope and redshift evolution of the MZR may provide important insight into the physical processes underpinning it: star formation, fee…
▽ More
Observations show a tight correlation between the stellar mass of galaxies and their gas-phase metallicity (MZR). This relation evolves with redshift, with higher-redshift galaxies being characterized by lower metallicities. Understanding the physical origin of the slope and redshift evolution of the MZR may provide important insight into the physical processes underpinning it: star formation, feedback, and cosmological inflows. While theoretical models ascribe the shape of the MZR to the lower efficiency of galactic outflows in more massive galaxies, what drives its evolution remains an open question. In this letter, we analyze how the MZR evolves over $z=0-3$, combining results from the FIREbox cosmological volume simulation with analytical models. Contrary to a frequent assertion in the literature, we find that the evolution of the gas fraction does not contribute significantly to the redshift evolution of the MZR. Instead, we show that the latter is driven by the redshift-dependence of the inflow metallicity, outflow metallicity, and mass loading factor, whose relative importance depends on stellar mass. These findings also suggest that the evolution of the MZR is not explained by galaxies moving along a fixed surface in the space spanned by stellar mass, gas phase metallicity, and star formation rate.
△ Less
Submitted 24 January, 2024;
originally announced January 2024.
-
Starburst-induced gas-stars kinematic misalignment
Authors:
Elia Cenci,
Robert Feldmann,
Jindra Gensior,
James S. Bullock,
Jorge Moreno,
Luigi Bassini,
Mauro Bernardini
Abstract:
A kinematic misalignment of the stellar and gas components is a phenomenon observed in a significant fraction of galaxies. However, the underlying physical mechanisms are not well understood. A commonly proposed scenario for the formation of a misaligned component requires any pre-existing gas disc to be removed, via fly-bys or ejective feedback from an active galactic nucleus. In this Letter, we…
▽ More
A kinematic misalignment of the stellar and gas components is a phenomenon observed in a significant fraction of galaxies. However, the underlying physical mechanisms are not well understood. A commonly proposed scenario for the formation of a misaligned component requires any pre-existing gas disc to be removed, via fly-bys or ejective feedback from an active galactic nucleus. In this Letter, we study the evolution of a Milky Way mass galaxy in the FIREbox cosmological volume that displays a thin, counter-rotating gas disc with respect to its stellar component at low redshift. In contrast to scenarios involving gas ejection, we find that pre-existing gas is mainly removed via the conversion into stars in a central starburst, triggered by a merging satellite galaxy. The newly-accreted, counter-rotating gas eventually settles into a kinematically misaligned disc. About 4.4 (8 out of 182) of FIREbox galaxies with stellar masses larger than 5e9 Msun at z=0 exhibit gas-star kinematic misalignment. In all cases, we identify central starburst-driven depletion as the main reason for the removal of the pre-existing co-rotating gas component, with no need for feedback from, e.g., a central active black hole. However, during the starburst, the gas is funneled towards the central regions, likely enhancing black hole activity. By comparing the fraction of misaligned discs between FIREbox and other simulations and observations, we conclude that this channel might have a non-negligible role in inducing kinematic misalignment in galaxies.
△ Less
Submitted 12 December, 2023;
originally announced December 2023.
-
Starbursts driven by central gas compaction
Authors:
Elia Cenci,
Robert Feldmann,
Jindra Gensior,
Jorge Moreno,
Luigi Bassini,
Mauro Bernardini
Abstract:
Starburst (SB) galaxies are a rare population of galaxies with star formation rates (SFRs) greatly exceeding those of the majority of star-forming galaxies with similar stellar mass. It is unclear whether these bursts are the result of either especially large gas reservoirs or enhanced efficiencies in converting gas into stars. Tidal torques resulting from gas-rich galaxy mergers are known to enha…
▽ More
Starburst (SB) galaxies are a rare population of galaxies with star formation rates (SFRs) greatly exceeding those of the majority of star-forming galaxies with similar stellar mass. It is unclear whether these bursts are the result of either especially large gas reservoirs or enhanced efficiencies in converting gas into stars. Tidal torques resulting from gas-rich galaxy mergers are known to enhance the SFR by funneling gas towards the centre. However, recent theoretical works show that mergers do not always trigger a SB and not all SB galaxies are interacting systems, raising the question of what drives a SB. We analyse a large sample of SB galaxies and a mass- and redshift-matched sample of control galaxies, drawn from the FIREbox cosmological volume at z=0-1. We find that SB galaxies have both larger molecular gas fractions and shorter molecular depletion times than control galaxies, but similar total gas masses. Control galaxies evolve towards the SB regime by gas compaction in their central regions, over timescales of about 70 Myr, accompanied by an increase in the fraction of ultra-dense and molecular gas. The driving mechanism behind the SB varies depending on the mass of the galaxy. Massive (Mstar > 1e10 Msun) galaxies undergoing intense, long-lasting SBs are mostly driven by galaxy interactions. Conversely, SBs in non-interacting galaxies are often triggered by a global gravitational instability, that can result in a breathing mode in low-mass galaxies.
△ Less
Submitted 16 September, 2023;
originally announced September 2023.
-
The inefficiency of stellar feedback in driving galactic outflows in massive galaxies at high redshift
Authors:
L. Bassini,
R. Feldmann,
J. Gensior,
C. C. Hayward,
C. -A. Faucher-Giguère,
E. Cenci,
L. Liang,
M. Bernardini
Abstract:
Recent observations indicate that galactic outflows are ubiquitous in high redshift galaxies, including normal star forming galaxies, quasar hosts, and dusty star forming galaxies (DSFGs). However, the impact of outflows on the evolution of their hosts is still an open question. Here, we analyse the star formation histories (SFH) and galactic outflow properties of galaxies in massive haloes (…
▽ More
Recent observations indicate that galactic outflows are ubiquitous in high redshift galaxies, including normal star forming galaxies, quasar hosts, and dusty star forming galaxies (DSFGs). However, the impact of outflows on the evolution of their hosts is still an open question. Here, we analyse the star formation histories (SFH) and galactic outflow properties of galaxies in massive haloes ($10^{12}M_{\odot}<M_{\rm vir} <5\times 10^{12}M_{\odot}$) at $z\gtrsim5.5$ in three zoom-in cosmological simulations from the MassiveFIRE suite, as part of the Feedback In Realistic Environments (FIRE) project. The simulations were run with the FIRE-2 model, which does not include feedback from active galactic nuclei (AGN). The simulated galaxies resemble $z>4$ DSFGs, with SFRs of $\sim 1000\ M_{\odot}\rm yr^{-1}$ and molecular gas masses of $M_{\rm mol}\sim 10^{10}\ M_{\odot}$. However, the simulated galaxies are characterised by higher circular velocities than those observed in high-z DSFGs. The mass loading factors from stellar feedback are of the order of $\sim 0.1$, implying that stellar feedback is inefficient in driving galactic outflows and gas is consumed by star formation on much shorter time-scales than it is expelled from the interstellar medium (ISM). We also find that stellar feedback is highly inefficient in self-regulating star formation in this regime, with an average integrated star formation efficiency (SFE) per dynamical time of $30\%$. Finally, compared to FIRE-2 galaxies hosted in similarly massive haloes at lower redshift, we find lower mass loading factors and higher SFEs in the high redshift sample. We argue that both effects originate from the higher total and gas surface densities that characterise high$-z$ massive systems.
△ Less
Submitted 15 November, 2022;
originally announced November 2022.
-
FIREbox: Simulating galaxies at high dynamic range in a cosmological volume
Authors:
Robert Feldmann,
Eliot Quataert,
Claude-André Faucher-Giguère,
Philip F. Hopkins,
Onur Çatmabacak,
Dušan Kereš,
Luigi Bassini,
Mauro Bernardini,
James S. Bullock,
Elia Cenci,
Jindra Gensior,
Lichen Liang,
Jorge Moreno,
Andrew Wetzel
Abstract:
We introduce a suite of cosmological volume simulations to study the evolution of galaxies as part of the Feedback in Realistic Environments project. FIREbox, the principal simulation of the present suite, provides a representative sample of galaxies (~1000 galaxies with Mstar > 10^8 Msun at z=0) at a resolution (~20 pc, m_b ~ 6x10^4 Msun) comparable to state-of-the-art galaxy zoom-in simulations.…
▽ More
We introduce a suite of cosmological volume simulations to study the evolution of galaxies as part of the Feedback in Realistic Environments project. FIREbox, the principal simulation of the present suite, provides a representative sample of galaxies (~1000 galaxies with Mstar > 10^8 Msun at z=0) at a resolution (~20 pc, m_b ~ 6x10^4 Msun) comparable to state-of-the-art galaxy zoom-in simulations. FIREbox captures the multiphase nature of the interstellar medium in a fully cosmological setting (L=22.1 Mpc) thanks to its exceptionally high dynamic range (~10^6) and the inclusion of multi-channel stellar feedback. Here, we focus on validating the simulation predictions by comparing to observational data. We find that simulated galaxies with Mstar < 10^{10.5-11} Msun have star formation rates, gas masses, and metallicities in broad agreement with observations. These galaxy scaling relations extend to low masses (Mstar ~ 10^7 Msun) and follow a (broken) power-law relationship. Also reproduced are the evolution of the cosmic HI density and the HI column density distribution at z~0-5. At low z, FIREbox predicts a peak in the stellar-mass--halo-mass relation, but also a higher abundance of massive galaxies and a higher cosmic star formation rate density than observed, showing that stellar feedback alone is insufficient to reproduce the properties of massive galaxies at late times. Given its high resolution and sample size, FIREbox offers a baseline prediction of galaxy formation theory in a $Λ$CDM Universe while also highlighting modeling challenges to be addressed in next-generation galaxy simulations.
△ Less
Submitted 21 April, 2023; v1 submitted 30 May, 2022;
originally announced May 2022.
-
Non-isotropic feedback from accreting spinning black holes
Authors:
Luca Sala,
Elia Cenci,
Pedro R. Capelo,
Alessandro Lupi,
Massimo Dotti
Abstract:
Active galactic nuclei (AGN) are massive black holes (BHs) caught in the act of accreting gas at the centre of their host galaxies. Part of the accreting mass is converted to energy and released into the surrounding medium, in a process loosely referred to as AGN feedback. Most numerical simulations include AGN feedback as a sub-grid model, wherein energy or momentum (or both) is coupled to the ne…
▽ More
Active galactic nuclei (AGN) are massive black holes (BHs) caught in the act of accreting gas at the centre of their host galaxies. Part of the accreting mass is converted to energy and released into the surrounding medium, in a process loosely referred to as AGN feedback. Most numerical simulations include AGN feedback as a sub-grid model, wherein energy or momentum (or both) is coupled to the nearby gas. In this work, we implement a new momentum-driven model in the hydrodynamics code GIZMO, in which accretion from large scales is mediated by a sub-grid accretion disc model, and gas particles are stochastically kicked over a bi-conical region, to mimic observed kinetic winds. The feedback cone's axis can be set either parallel to the angular momentum of the gas surrounding the BH or to the BH spin direction, which is self-consistently evolved within the accretion-disc model. Using a circumnuclear disc (CND) as a test bed, we find that (i) the conical shape of the outflow is always visible and is weakly dependent on the launching orientation and aperture, resulting in comparable mass inflows and outflows; (ii) the cone's orientation is also similar amongst our tests, and it is not always the same as the initial value, due to the interaction with the CND playing a crucial role in shaping the outflow; and (iii) the velocity of the outflow, instead, differs and strongly depends on the interplay with the CND.
△ Less
Submitted 13 April, 2022; v1 submitted 12 November, 2020;
originally announced November 2020.
-
Black hole spin evolution in warped accretion discs
Authors:
Elia Cenci,
Luca Sala,
Alessandro Lupi,
Pedro R. Capelo,
Massimo Dotti
Abstract:
Massive black holes (BHs) inhabiting galactic nuclei can be described by two parameters only, i.e. mass and spin, that change through cosmic time in response to accretion and merger events. While most numerical simulations accurately track the BH mass, spin evolution is rarely taken into account. In this work, we implement and validate a self-consistent sub-grid model for the evolution of the BH m…
▽ More
Massive black holes (BHs) inhabiting galactic nuclei can be described by two parameters only, i.e. mass and spin, that change through cosmic time in response to accretion and merger events. While most numerical simulations accurately track the BH mass, spin evolution is rarely taken into account. In this work, we implement and validate a self-consistent sub-grid model for the evolution of the BH mass and spin via gas accretion in the hydrodynamics code GIZMO. The model assumes that accretion from resolved scales does not occur instantaneously, but is mediated by a sub-grid geometrically thin $α$-disc. After validating our model semi-analytically, we test it in an idealized environment consisting of a circumnuclear disc, where gas accretion onto the accretion disc is consistently determined by GIZMO. In the absence of any accretion-related feedback, the spin evolution closely traces that observed in the semi-analytical models and depends on the free parameters of our implementation, such as the initial BH spin, angular momentum of the accretion disc, and the radius at which the gas inflow circularises. In GIZMO, we also couple our model with the biconical-outflow model presented in a companion paper, wherein the feedback axis is always aligned with the BH spin. In this last case, the evolution of the central BH differs significantly from the previous cases, since the feedback process modifies the gas dynamics and its inflow rates from resolved scales. Such an interaction cannot be modeled by simple semi-analytical models and should be treated using full N-body hydrodynamical simulations.
△ Less
Submitted 13 April, 2022; v1 submitted 12 November, 2020;
originally announced November 2020.