The mass-metallicity relation of tidal dwarf galaxies

Dwarf galaxies generally follow a mass-metallicity (MZ) relation, where more massive objects retain a larger fraction of heavy elements. Young tidal dwarf galaxies (TDGs), born in the tidal tails produced by interacting gas-rich galaxies, have been thought to not follow the MZ relation, because they inherit the metallicity of the more massive parent galaxies. We present chemical evolution models to investigate if TDGs that formed at very high redshifts, where the metallicity of their parent galaxy was very low, can produce the observed MZ relation. Assuming that galaxy interactions were more frequent in the denser high-redshift universe, TDGs could constitute an important contribution to the dwarf galaxy population. The survey of chemical evolution models of TDGs presented here captures for the first time an initial mass function (IMF) of stars that is dependent on both the star formation rate and the gas metallicity via the integrated galactic IMF (IGIMF) theory. As TDGs form in the tidal debris of interacting galaxies, the pre-enrichment of the gas, an underlying pre-existing stellar population, infall, and mass dependent outflows are considered. The models of young TDGs that are created in strongly pre-enriched tidal arms with a pre-existing stellar population can explain the measured abundance ratios of observed TDGs. The same chemical evolution models for TDGs, that form out of gas with initially very low metallicity, naturally build up the observed MZ relation. The modelled chemical composition of ancient TDGs is therefore consistent with the observed MZ relation of satellite galaxies.Comment: 7 pages, 3 figures, MNRAS accepte

Resolution criteria to avoid artificial clumping in Lagrangian hydrodynamic simulations with a multi-phase interstellar medium

Large-scale cosmological galaxy formation simulations typically prevent gas in the interstellar medium (ISM) from cooling below $\approx 10^4$ K. This has been motivated by the inability to resolve the Jeans mass in molecular gas (>>$10^5\,\mathrm{M}_{\odot}$) which would result in undesired artificial clumping. We show that the classical Jeans criteria derived for Newtonian gravity are not applicable in the simulated ISM if the spacing of resolution elements representing the dense ISM is below the gravitational force softening length and gravity is therefore softened and not Newtonian. We re-derive the Jeans criteria for softened gravity in Lagrangian codes and use them to analyse gravitational instabilities at and below the hydrodynamical resolution limit for simulations with adaptive and constant gravitational softening lengths. In addition, we define criteria for which a numerical runaway collapse of dense gas clumps can occur caused by over-smoothing of the hydrodynamical properties relative to the gravitational force resolution. This effect is illustrated using simulations of isolated disk galaxies with the smoothed particle hydrodynamics code Swift. We also demonstrate how to avoid the formation of artificial clumps in gas and stars by adjusting the gravitational and hydrodynamical force resolutions.Comment: 24 pages, 15 figures, accepted for publication in MNRAS, smaller updates to match published versio

A thermal-kinetic subgrid model for supernova feedback in simulations of galaxy formation

We present a subgrid model for supernova feedback designed for simulations of galaxy formation. The model uses thermal and kinetic channels of energy injection, which are built upon the stochastic kinetic and thermal models for stellar feedback used in the OWLS and EAGLE simulations, respectively. In the thermal channel, the energy is distributed statistically isotropically and injected stochastically in large amounts per event, which minimizes spurious radiative energy losses. In the kinetic channel, we inject the energy in small portions by kicking gas particles in pairs in opposite directions. The implementation of kinetic feedback is designed to conserve energy, linear momentum and angular momentum, and is statistically isotropic. To test and validate the model, we run simulations of isolated Milky Way-mass and dwarf galaxies, in which the gas is allowed to cool down to 10 K. Using the thermal and kinetic channels together, we obtain smooth star formation histories and powerful galactic winds with realistic mass loading factors. Furthermore, the model produces spatially resolved star formation rates and velocity dispersions that are in agreement with observations. We vary the numerical resolution by several orders of magnitude and find excellent convergence of the global star formation rates and the mass loading of galactic winds. We show that large thermal-energy injections generate a hot phase of the interstellar medium (ISM) and modulate the star formation by ejecting gas from the disc, while the low-energy kicks increase the turbulent velocity dispersion in the neutral ISM, which in turn helps suppress star formation.Comment: 22 pages, 17 figures (including appendix); submitted to MNRA

Tests of subgrid models for star formation using simulations of isolated disk galaxies

We use smoothed-particle hydrodynamics simulations of isolated Milky Way-mass disk galaxies that include cold, interstellar gas to test subgrid prescriptions for star formation (SF). Our fiducial model combines a Schmidt law with a gravitational instability criterion, but we also test density thresholds and temperature ceilings. While SF histories are insensitive to the prescription for SF, the Kennicutt-Schmidt (KS) relations between SF rate and gas surface density can discriminate between models. We show that our fiducial model, with an SF efficiency per free-fall time of 1 per cent, agrees with spatially-resolved and azimuthally-averaged observed KS relations for neutral, atomic and molecular gas. Density thresholds do not perform as well. While temperature ceilings selecting cold, molecular gas can match the data for galaxies with solar metallicity, they are unsuitable for very low-metallicity gas and hence for cosmological simulations. We argue that SF criteria should be applied at the resolution limit rather than at a fixed physical scale, which means that we should aim for numerical convergence of observables rather than of the properties of gas labelled as star-forming. Our fiducial model yields good convergence when the mass resolution is varied by nearly 4 orders of magnitude, with the exception of the spatially-resolved molecular KS relation at low surface densities. For the gravitational instability criterion, we quantify the impact on the KS relations of gravitational softening, the SF efficiency, and the strength of supernova feedback, as well as of observable parameters such as the inclusion of ionized gas, the averaging scale, and the metallicity.Comment: Submitted to MNRAS, 23 pages, 20 figure

Diagnosing the interstellar medium of galaxies with far-infrared emission lines I. The [C II] 158 microns line at z~0

Atomic fine structure lines have been detected in the local Universe and at high redshifts over the past decades. The [C II] emission line at 158 $\mu$m is an important observable as it provides constraints on the interstellar medium (ISM) cooling processes. We develop a physically motivated framework to simulate the production of far-infrared line emission from galaxies in a cosmological context. This first paper sets out our methodology and describes its first application, simulating the [C II] 158 $\mu$m line emission in the local Universe. We combine the output from EAGLE cosmological hydrodynamical simulations with a multi-phase model of the ISM. Gas particles are divided into three phases: dense molecular gas, neutral atomic gas and diffuse ionised gas (DIG). We estimate the [C II] line emission from the three phases using a set of Cloudy cooling tables. Our results agree with previous findings regarding the contribution of these three ISM phases to the [C II] emission. Our model shows good agreement with the observed ${\rm L_{[C II]}}$-star formation rate (SFR) relation in the local Universe within 0.4 dex scatter. The fractional contribution to the [C II] line from different ISM phases depends on the total SFR and metallicity. The neutral gas phase dominates the [C II] emission in galaxies with $\rm{SFR}\sim0.01$-$1\,\rm{M_{\odot}}\,\rm{yr^{-1}}$, but the ionised phase dominates at lower SFRs. Galaxies above solar metallicity exhibit lower ${\rm L_{[C II]}}$/SFR ratios for the neutral phase. In comparison, the ${\rm L_{[C II]}}$/SFR ratio in the DIG is stable when metallicity varies. We suggest that the reduced size of the neutral clouds, caused by increased SFRs, is the likely cause for the ${\rm L_{[C II]}}$ deficit at high infrared luminosities, although EAGLE simulations do not reach these luminosities at $z=0$.Comment: Accepted for publication in A&A; 16 pages, 11 figures, 4 tables (plus appendix

The impact of stochastic modeling on the predictive power of galaxy formation simulations

All modern galaxy formation models employ stochastic elements in their sub-grid prescriptions to discretise continuous equations across the time domain. In this paper, we investigate how the stochastic nature of these models, notably star formation, black hole accretion, and their associated feedback, that act on small ($<$ kpc) scales, can back-react on macroscopic galaxy properties (e.g. stellar mass and size) across long ($>$ Gyr) timescales. We find that the scatter in scaling relations predicted by the EAGLE model implemented in the SWIFT code can be significantly impacted by random variability between re-simulations of the same object, even when galaxies are resolved by tens of thousands of particles. We then illustrate how re-simulations of the same object can be used to better understand the underlying model, by showing how correlations between galaxy stellar mass and black hole mass disappear at the highest black hole masses ($M_{\rm BH} > 10^8$ M$_\odot$), indicating that the feedback cycle may be interrupted by external processes. We find that although properties that are collected cumulatively over many objects are relatively robust against random variability (e.g. the median of a scaling relation), the properties of individual galaxies (such as galaxy stellar mass) can vary by up to 25\%, even far into the well-resolved regime, driven by bursty physics (black hole feedback) and mergers between galaxies. We suggest that studies of individual objects within cosmological simulations be treated with caution, and that any studies aiming to closely investigate such objects must account for random variability within their results.Comment: Accepted for publication in MNRA

The importance of the way in which supernova energy is distributed around young stellar populations in simulations of galaxies

Horizon 2020(H2020)860744Galaxie

TangoSIDM: tantalizing models of self-interacting dark matter

Large scale structure and cosmologyGalaxie

Tests of subgrid models for star formation using simulations of isolated disk galaxies

We use smoothed-particle hydrodynamics simulations of isolated Milky Way-mass disk galaxies that include cold, interstellar gas to test subgrid prescriptions for star formation (SF). Our fiducial model combines a Schmidt law with a gravitational instability criterion, but we also test density thresholds and temperature ceilings. While SF histories are insensitive to the prescription for SF, the Kennicutt-Schmidt (KS) relations between SF rate and gas surface density can discriminate between models. We show that our fiducial model, with an SF efficiency per free-fall time of 1 per cent, agrees with spatially-resolved and azimuthally-averaged observed KS relations for neutral, atomic and molecular gas. Density thresholds do not perform as well. While temperature ceilings selecting cold, molecular gas can match the data for galaxies with solar metallicity, they are unsuitable for very low-metallicity gas and hence for cosmological simulations. We argue that SF criteria should be applied at the resolution limit rather than at a fixed physical scale, which means that we should aim for numerical convergence of observables rather than of the properties of gas labelled as star-forming. Our fiducial model yields good convergence when the mass resolution is varied by nearly 4 orders of magnitude, with the exception of the spatially-resolved molecular KS relation at low surface densities. For the gravitational instability criterion, we quantify the impact on the KS relations of gravitational softening, the SF efficiency, and the strength of supernova feedback, as well as of observable parameters such as the inclusion of ionized gas, the averaging scale, and the metallicity