No Arabic abstract
Recent spatially resolved observations of galaxies at z=0.6-3 reveal that high-redshift galaxies show complex kinematics and a broad distribution of gas-phase metallicity gradients. To understand these results, we use a suite of high-resolution cosmological zoom-in simulations from the Feedback in Realistic Environments (FIRE) project, which include physically motivated models of the multi-phase ISM, star formation, and stellar feedback. Our simulations reproduce the observed diversity of kinematic properties and metallicity gradients, broadly consistent with observations at z=0-3. Strong negative metallicity gradients only appear in galaxies with a rotating disk, but not all rotationally supported galaxies have significant gradients. Strongly perturbed galaxies with little rotation always have flat gradients. The kinematic properties and metallicity gradient of a high-redshift galaxy can vary significantly on short time-scales, associated with starburst episodes. Feedback from a starburst can destroy the gas disk, drive strong outflows, and flatten a pre-existing negative metallicity gradient. The time variability of a single galaxy is statistically similar to the entire simulated sample, indicating that the observed metallicity gradients in high-redshift galaxies reflect the instantaneous state of the galaxy rather than the accretion and growth history on cosmological time-scales. We find weak dependence of metallicity gradient on stellar mass and specific star formation rate (sSFR). Low-mass galaxies and galaxies with high sSFR tend to have flat gradients, likely due to the fact that feedback is more efficient in these galaxies. We argue that it is important to resolve feedback on small scales in order to produce the diverse metallicity gradients observed.
We present a new model for the evolution of gas phase metallicity gradients in galaxies from first principles. We show that metallicity gradients depend on four ratios that collectively describe the metal equilibration timescale, production, transport, consumption, and loss. Our model finds that most galaxy metallicity gradients are in equilibrium at all redshifts. When normalized by metal diffusion, metallicity gradients are governed by the competition between radial advection, metal production, and accretion of metal-poor gas from the cosmic web. The model naturally explains the varying gradients measured in local spirals, local dwarfs, and high-redshift star-forming galaxies. We use the model to study the cosmic evolution of gradients across redshift, showing that the gradient in Milky Way-like galaxies has steepened over time, in good agreement with both observations and simulations. We also predict the evolution of metallicity gradients with redshift in galaxy samples constructed using both matched stellar masses and matched abundances. Our model shows that massive galaxies transition from the advection-dominated to the accretion-dominated regime from high to low redshifts, which mirrors the transition from gravity-driven to star formation feedback-driven turbulence. Lastly, we show that gradients in local ultraluminous infrared galaxies (major mergers) and inverted gradients seen both in the local and high-redshift galaxies may not be in equilibrium. In subsequent papers in this series, we show that the model also explains the observed relationship between galaxy mass and metallicity gradients, and between metallicity gradients and galaxy kinematics.
We examine radial and vertical metallicity gradients using a suite of disk galaxy simulations, supplemented with two classic chemical evolution approaches. We determine the rate of change of gradient and reconcile differences between extant models and observations within the `inside-out disk growth paradigm. A sample of 25 disks is used, consisting of 19 from our RaDES (Ramses Disk Environment Study) sample, realised with the adaptive mesh refinement code RAMSES. Four disks are selected from the MUGS (McMaster Unbiased Galaxy Simulations) sample, generated with the smoothed particle hydrodynamics (SPH) code GASOLINE, alongside disks from Rahimi et al. (GCD+) and Kobayashi & Nakasato (GRAPE-SPH). Two chemical evolution models of inside-out disk growth were employed to contrast the temporal evolution of their radial gradients with those of the simulations. We find that systematic differences exist between the predicted evolution of radial abundance gradients in the RaDES and chemical evolution models, compared with the MUGS sample; specifically, the MUGS simulations are systematically steeper at high-redshift, and present much more rapid evolution in their gradients. We find that the majority of the models predict radial gradients today which are consistent with those observed in late-type disks, but they evolve to this self-similarity in different fashions, despite each adhering to classical `inside-out growth. We find that radial dependence of the efficiency with which stars form as a function of time drives the differences seen in the gradients; systematic differences in the sub-grid physics between the various codes are responsible for setting these gradients. Recent, albeit limited, data at redshift z=1.5 are consistent with the steeper gradients seen in our SPH sample, suggesting a modest revision of the classical chemical evolution models may be required.
In this work, we explore the diversity of ionised gas kinematics (rotational velocity $v_{phi}$ and velocity dispersion $sigma_{mathrm{g}}$) and gas-phase metallicity gradients at $0.1 leq z leq 2.5$ using a compiled data set of 74 galaxies resolved with ground-based integral field spectroscopy. We find that galaxies with the highest and the lowest $sigma_{mathrm{g}}$ have preferentially flat metallicity gradients, whereas those with intermediate values of $sigma_{mathrm{g}}$ show a large scatter in the metallicity gradients. Additionally, steep negative gradients appear almost only in rotation-dominated galaxies ($v_{phi}/sigma_{mathrm{g}} > 1$), whereas most dispersion-dominated galaxies show flat gradients. We use our recently developed analytic model of metallicity gradients to provide a physical explanation for the shape and scatter of these observed trends. In the case of high $sigma_{mathrm{g}}$, the inward radial advection of gas dominates over metal production and causes efficient metal mixing, thus giving rise to flat gradients. For low $sigma_{mathrm{g}}$, it is the cosmic accretion of metal-poor gas diluting the metallicity that gives rise to flat gradients. Finally, the reason for intermediate $sigma_{mathrm{g}}$ showing the steepest negative gradients is that both inward radial advection and cosmic accretion are weak as compared to metal production, which leads to the creation of steeper gradients. The larger scatter at intermediate $sigma_{mathrm{g}}$ may be due in part to preferential ejection of metals in galactic winds, which can decrease the strength of the production term. Our analysis shows how gas kinematics play a critical role in setting metallicity gradients in high-redshift galaxies.
Within the standard model of hierarchical galaxy formation in a {Lambda}CDM Universe, the environment of galaxies is expected to play a key role in driving galaxy formation and evolution. In this paper we investigate whether and how the gas metallicity and the star formation surface density ({Sigma}_SFR) depend on galaxy environment. To this end we analyse a sample of 1162 local, star-forming galaxies from the galaxy survey Mapping Nearby Galaxies at APO (MaNGA). Generally, both parameters do not show any significant dependence on environment. However, in agreement with previous studies, we find that low-mass satellite galaxies are an exception to this rule. The gas metallicity in these objects increases while their {Sigma}SFR decreases slightly with environmental density. The present analysis of MaNGA data allows us to extend this to spatially resolved properties. Our study reveals that the gas metallicity gradients of low-mass satellites flatten and their {Sigma}SFR gradients steepen with increasing environmental density. By extensively exploring a chemical evolution model, we identify two scenarios that are able to explain this pattern: metal-enriched gas accretion or pristine gas inflow with varying accretion timescales. The latter scenario better matches the observed {Sigma}SFR gradients, and is therefore our preferred solution. In this model, a shorter gas accretion timescale at larger radii is required. This suggests that outside-in quenching governs the star formation processes of low-mass satellite galaxies in dense environments.
Galaxies at low-redshift typically possess negative gas-phase metallicity gradients (centres more metal-rich than their outskirts). Whereas, it is not uncommon to observe positive metallicity gradients in higher-redshift galaxies ($z gtrsim 0.6$). Bridging these epochs, we present gas-phase metallicity gradients of 84 star-forming galaxies between $0.08 < z < 0.84$. Using the galaxies with reliably determined metallicity gradients, we measure the median metallicity gradient to be negative ($-0.039^{+0.007}_{-0.009}$ dex/kpc). Underlying this, however, is significant scatter: $(8pm3)% [7]$ of galaxies have significantly positive metallicity gradients, $(38 pm 5)% [32]$ have significantly negative gradients, $(31pm5)% [26]$ have gradients consistent with being flat. (The remaining $(23pm5)% [19]$ have unreliable gradient estimates.) We notice a slight trend for a more negative metallicity gradient with both increasing stellar mass and increasing star formation rate (SFR). However, given the potential redshift and size selection effects, we do not consider these trends to be significant. Indeed, once we normalize the SFR relative to that of the main sequence, we do not observe any trend between the metallicity gradient and the normalized SFR. This is contrary to recent studies of galaxies at similar and higher redshifts. We do, however, identify a novel trend between the metallicity gradient of a galaxy and its size. Small galaxies ($r_d < 3$ kpc) present a large spread in observed metallicity gradients (both negative and positive gradients). In contrast, we find no large galaxies ($r_d > 3$ kpc) with positive metallicity gradients, and overall there is less scatter in the metallicity gradient amongst the large galaxies. These large (well-evolved) galaxies may be analogues of present-day galaxies, which also show a common negative metallicity gradient.