Stellar Mass-metallicity Relation Research Articles

Metallicities of young massive clusters in NGC 5236 (M83)

We present integrated-light spectra of 8 Young Massive Clusters (YMCs) in the metal-rich spiral galaxy NGC 5236 (M 83). The observations were taken with the X-Shooter spectrograph on the ESO Very Large Telescope. Through the use of theoretical isochrones and synthetic integrated-light (IL) spectra we derive metallicities and study the radial metallicity gradient observed through these young populations. For the inner regions of the galaxy we observe a relatively shallow metallicity gradient of $-$0.37 $\pm$0.29 dex R$_{25}^{-1}$, agreeing with chemical evolution models with an absence of infall material and a relatively low mass loss due to winds in the inner parts of the disk. We estimate a central metallicity of [$Z$] = $+$0.17 $\pm$ 0.12 dex, finding excellent agreement with that obtained via other methods (e.g. blue supergiants and J-band). We infer a metallicity of 12+log(O/H) = 8.75 $\pm$ 0.08 dex at R/R$_{25}$ = 0.4, which fits the stellar mass-metallicity relation (MZR) compilation of blue supergiants and IL studies.

Separate Ways: The Mass–Metallicity Relation Does Not Strongly Correlate with Star Formation Rate in SDSS-IV MaNGA Galaxies

We present the integrated stellar mass-metallicity relation (MZR) for more than 1700 galaxies included in the integral field area SDSS-IV MaNGA survey. The spatially resolved data allow us to determine the metallicity at the same physical scale (effective radius in arcsecs, $\mathrm{R_{eff}}$ ) using a heterogeneous set of ten abundance calibrators. Besides scale factors, the shape of the MZR is similar for all calibrators, consistent with those reported previously using single-fiber and integral field spectroscopy. We compare the residuals of this relation against the star formation rate (SFR) and specific SFR (sSFR). We do not find a strong secondary relation of the MZR with either SFR or the sSFR for any of the calibrators, in contrast with previous single-fiber spectroscopic studies. Our results agree with an scenario in which metal enrichment happens at local scales, with global outflows playing a secondary role in shaping the chemistry of galaxies and cold-gas inflows regulating the stellar formation.

Mufasa: the assembly of the red sequence

We examine the growth and evolution of quenched galaxies in the Mufasa cosmological hydrodynamic simulations that include an evolving halo mass-based quenching prescription, with galaxy colours computed accounting for line-of-sight extinction to individual star particles. Mufasa reproduces the observed present-day red sequence reasonably well, including its slope, amplitude, and scatter. In Mufasa, the red sequence slope is driven entirely by the steep stellar mass-stellar metallicity relation, which independently agrees with observations. High-mass star-forming galaxies blend smoothly onto the red sequence, indicating the lack of a well-defined green valley at M*>10^10.5 Mo. The most massive galaxies quench the earliest and then grow very little in mass via dry merging; they attain their high masses at earlier epochs when cold inflows more effectively penetrate hot halos. To higher redshifts, the red sequence becomes increasingly contaminated with massive dusty star-forming galaxies; UVJ selection subtly but effectively separates these populations. We then examine the evolution of the mass functions of central and satellite galaxies split into passive and star-forming via UVJ. Massive quenched systems show good agreement with observations out to z~2, despite not including a rapid early quenching mode associated with mergers. However, low-mass quenched galaxies are far too numerous at z<1 in Mufasa, indicating that Mufasa strongly over-quenches satellites. A challenge for hydrodynamic simulations is to devise a quenching model that produces enough early massive quenched galaxies and keeps them quenched to z=0, while not being so strong as to over-quench satellites; Mufasa's current scheme fails at the latter.

Radiation Backgrounds at Cosmic Dawn: X-Rays from Compact Binaries

We compute the expected X-ray diffuse background and radiative feedback on the intergalactic medium (IGM) from X-ray binaries prior and during the epoch of reionization. The cosmic evolution of compact binaries is followed using a population synthesis technique that treats separately neutron stars and black hole binaries in different spectral states and is calibrated to reproduce the observed X-ray properties of galaxies at z<4. Together with an updated empirical determination of the cosmic history of star formation, recent modeling of the stellar mass-metallicity relation, and a scheme for absorption by the IGM that accounts for the presence of ionized HII bubbles during the epoch of reionization, our detailed calculations provide refined predictions of the X-ray volume emissivity and filtered radiation background from "normal" galaxies at z>6. Radiative transfer effects modulate the background spectrum, which shows a characteristic peak between 1 and 2 keV. While the filtering of X-ray radiation through the IGM slightly increases the mean excess energy per photoionization, it also weakens the radiation intensity below 1 keV, lowering the mean photoionization and heating rates. Numerical integration of the rate and energy equations shows that the contribution of X-ray binaries to the ionization of the bulk IGM is negligible, with the electron fraction never exceeding 1%. Direct HeI photoionizations are the main source of IGM heating, and the temperature of the largely neutral medium in between HII cavities increases above the temperature of the cosmic microwave background (CMB) only at z<10, when the volume filling factor of HII bubbles is already >0.1. Therefore, in this scenario, it is only at relatively late epochs that the bulk of neutral intergalactic hydrogen may be observable in 21-cm emission against the CMB.

COLD-MODE ACCRETION: DRIVING THE FUNDAMENTAL MASS–METALLICITY RELATION AT z ∼ 2

ABSTRACT We investigate the star formation rate (SFR) dependence on the stellar mass and gas-phase metallicity relation at z = 2 with MOSFIRE/Keck as part of the ZFIRE survey. We have identified 117 galaxies (1.98 ≤ z ≤ 2.56), with 8.9 ≤ log(M/M ⊙) ≤ 11.0, for which we can measure gas-phase metallicities. For the first time, we show a discernible difference between the mass–metallicity relation, using individual galaxies, when dividing the sample by low (&lt;10 M ⊙ yr−1) and high (&gt;10 M ⊙ yr−1) SFRs. At fixed mass, low star-forming galaxies tend to have higher metallicity than high star-forming galaxies. Using a few basic assumptions, we further show that the gas masses and metallicities required to produce the fundamental mass–metallicity relation and its intrinsic scatter are consistent with cold-mode accretion predictions obtained from the OWLS hydrodynamical simulations. Our results from both simulations and observations are suggestive that cold-mode accretion is responsible for the fundamental mass–metallicity relation at z = 2 and it demonstrates the direct relationship between cosmological accretion and the fundamental properties of galaxies.

Galaxies in the EAGLE hydrodynamical simulation and in the Durham and Munich semi-analytical models

We compare global predictions from the EAGLE hydrodynamical simulation, and two semi-analytic (SA) models of galaxy formation, L-GALAXIES and GALFORM. All three models include the key physical processes for the formation and evolution of galaxies and their parameters are calibrated against a small number of observables at $z\approx 0$. The two SA models have been applied to merger trees constructed from the EAGLE dark matter only simulation. We find that at $z\leq 2$, both the galaxy stellar mass functions for stellar masses $M_*<10^{10.5} M_{\odot}$ and the median specific star formation rates (sSFRs) in the three models agree to better than $0.4$~dex. The evolution of the sSFR predicted by the three models closely follows the mass assembly history of dark matter haloes. In both EAGLE and L-GALAXIES there are more central passive galaxies with $M_*<10^{9.5} M_{\odot}$ than in L-GALAXIES. This difference is related to galaxies that have entered and then left a larger halo and which are treated as satellites in GALFORM. In the range $0<z<1$, the slope of the evolution of the star formation rate density in EAGLE is a factor of $\approx 1.5$ steeper than for the two SA models. The median sizes for galaxies with $M_*>10^{9.5} M_{\odot}$ differ in some instances by an order of magnitude, while the stellar mass-size relation in EAGLE is a factor of $\approx 2$ tighter than for the two SA models. Our results suggest the need for a revision of how SA models treat the effect of baryonic self-gravity on the underlying dark matter. The treatment of gas flows in the models needs to be revised based on detailed comparison with observations to understand in particular the evolution of the stellar mass-metallicity relation.

A MULTIPLE BURST ACCRETION MODEL TO DESCRIBE THE METALLICITY DISTRIBUTIONS AND STELLAR MASS–METALLICITY RELATION FOR LOCAL DWARF GALAXIES

A one parameter model to describe the individual metallicity distributions and stellar mass-metallicity relation for dwarf galaxies is presented. This multiple-burst model is based on an accretion scenario, accomodates the observational constraint between $\overline{z}$ and $\sigma_{z}^{2}$ recently established by Leaman (2012), and predicts a slope consistent with the stellar mass-metallicity relation of Kirby et al (2013) who showed that the local group dwarf spheroidal and dwarf irregular galaxies lie on the same relation. One interpretation of the model is that it describes star formation occuring either in gas rich mergers or at the intersection of colliding gas streams.

Linking dust emission to fundamental properties in galaxies: the low-metallicity picture

In this work, we aim at providing a consistent analysis of the dust properties from metal-poor to metal-rich environments by linking them to fundamental galactic parameters. We consider two samples of galaxies: the Dwarf Galaxy Survey (DGS) and KINGFISH, totalling 109 galaxies, spanning almost 2 dex in metallicity. We collect infrared (IR) to submillimetre (submm) data for both samples and present the complete data set for the DGS sample. We model the observed spectral energy distributions (SED) with a physically-motivated dust model to access the dust properties. Using a different SED model (modified blackbody), dust composition (amorphous carbon), or wavelength coverage at submm wavelengths results in differences in the dust mass estimate of a factor two to three, showing that this parameter is subject to non-negligible systematic modelling uncertainties. For eight galaxies in our sample, we find a rather small excess at 500 microns (< 1.5 sigma). We find that the dust SED of low-metallicity galaxies is broader and peaks at shorter wavelengths compared to more metal-rich systems, a sign of a clumpier medium in dwarf galaxies. The PAH mass fraction and the dust temperature distribution are found to be driven mostly by the specific star-formation rate, SSFR, with secondary effects from metallicity. The correlations between metallicity and dust mass or total-IR luminosity are direct consequences of the stellar mass-metallicity relation. The dust-to-stellar mass ratios of metal-rich sources follow the well-studied trend of decreasing ratio for decreasing SSFR. The relation is more complex for highly star-forming low-metallicity galaxies and depends on the chemical evolutionary stage of the source (i.e., gas-to-dust mass ratio). Dust growth processes in the ISM play a key role in the dust mass build-up with respect to the stellar content at high SSFR and low metallicity. (abridged)

Star formation in semi-analytic galaxy formation models with multiphase gas

We implement physically motivated recipes for partitioning cold gas into different phases (atomic, molecular, and ionized) in galaxies within semi-analytic models of galaxy formation based on cosmological merger trees. We then model the conversion of molecular gas into stars using empirical recipes motivated by recent observations. We explore the impact of these new recipes on the evolution of fundamental galaxy properties such as stellar mass, star formation rate (SFR), and gas and stellar phase metallicity. We present predictions for stellar mass functions, stellar mass versus SFR relations, and cold gas phase and stellar mass–metallicity relations for our fiducial models, from redshift z ∼ 6 to the present day. In addition we present predictions for the global SFR, mass assembly history, and cosmic enrichment history. We find that the predicted stellar properties of galaxies (stellar mass, SFR, metallicity) are remarkably insensitive to the details of the recipes used for partitioning gas into H i and H2. We see significant sensitivity to the recipes for H2 formation only in very low mass haloes (⁠|$M_{\rm h} \lesssim 10^{10.5}\, {{\rm M}_{{\odot }}}$|⁠), which host galaxies with stellar masses |$m_* \lesssim 10^8\, {{\rm M}_{{\odot }}}$|⁠. The properties of low-mass galaxies are also quite insensitive to the details of the recipe used for converting H2 into stars, while the formation epoch of massive galaxies does depend on this significantly. We argue that this behaviour can be interpreted within the framework of a simple equilibrium model for galaxy evolution, in which the conversion of cold gas into stars is balanced on average by inflows and outflows.

Metallicities in cosmological simulations with AGN feedback

AbstractIn our cosmological, chemodynamical simulations, (i) the black hole mass–velocity dispersion relation does not evolve, and black holes actually grow along the relation. (ii) the stellar mass–metallicity relation does not change its shape, while the gas-phase relation has a steeper slope at higher redshifts. (iii) While stellar metallicity gradients are made shallower by galaxy mergers, gas-phase gradients are affected more strongly by AGN feedback.

AN ANALYTICAL MODEL FOR GALAXY METALLICITY: WHAT DO METALLICITY RELATIONS TELL US ABOUT STAR FORMATION AND OUTFLOW?

We develop a simple analytical model that tracks galactic metallicities governed by star formation and feedback to gain insight from the observed galaxy stellar mass-metallicity relations over a large range of stellar masses and redshifts. The model reveals the following implications of star formation and feedback processes in galaxy formation. First, the observed metallicity relations provide a stringent upper limit for the averaged outflow mass-loading factors of local galaxies, which is ~20 for M_*~10^9Msun galaxies and monotonically decreases to ~1 for M_*~10^{11}Msun galaxies. Second, the inferred upper-limit for the outflow mass-loading factor sensitively depends on whether the outflow is metal-enriched with respect to the ISM metallicity. If half of the metals ejected from SNe leave the galaxy in metal-enriched winds, the outflow mass-loading factor for galaxies at any mass can barely be higher than ~10, which puts strong constraints on galaxy formation models. Third, the relatively lower stellar-phase to gas-phase metallicity ratio for lower-mass galaxies indicate that low-mass galaxies are still rapidly enriching their metallicities in recent times, while high-mass galaxies are more settled, which seems to show a downsizing effect in the metallicity evolution of galaxies. The analysis presented in the paper demonstrates the importance of accurate measurements of galaxy metallicities and the cold gas fraction of galaxies at different redshifts for constraining star formation and feedback processes, and demonstrates the power of these relations in constraining the physics of galaxy formation.

Chemical evolution on the scale of clusters of galaxies: a conundrum?

The metal content of clusters of galaxies and its relation to their stellar content is revisited making use of a cluster sample for which all four basic parameters are homogeneously measured within consistent radii, namely core-excised mass-weighted metallicity plus total, stellar and ICM masses. For clusters of total mass $M_{500} >$ $\sim 10^{14}$ $M_{\odot}$ nice agreement is found between their iron content and what expected from empirical supernova yields. For the same clusters, there also appears to be at least as much iron in the intracluster medium (ICM) as there is still locked into stars (i.e., the ICM/stars metal share is about unity). However, for more massive clusters the stellar mass fraction appears to drop substantially without being accompanied by a drop in the ICM metallicity, thus generating a major tension with the nucleosynthesis expectation and inflating the metal share to extremely high values (up to $\sim 6$). Various possible solutions of this conundrum are discussed, but are all considered either astrophysically implausible, or lacking an independent observational support. For this reason we still entertain the possibility that even some of the best cluster data may be faulty, though we are not able to identify any obvious bias. Finally, based on the stellar mass-metallicity relation for local galaxies we estimate the contribution of galaxies to the ICM enrichment as a function of their mass, concluding that even the most massive galaxies must have lost a major fraction of the metals they have produced.

INSIGHTS ON THE STELLAR MASS-METALLICITY RELATION FROM THE CALIFA SURVEY

We use spatially and temporally resolved maps of stellar population properties of 300 galaxies from the CALIFA integral field survey to investigate how the stellar metallicity (Z*) relates to the total stellar mass (M*) and the local mass surface density ($\mu$*) in both spheroidal and disk dominated galaxies. The galaxies are shown to follow a clear stellar mass-metallicity relation (MZR) over the whole 10$^9$ to 10$^{12}$ M$_{\odot}$ range. This relation is steeper than the one derived from nebular abundances, which is similar to the flatter stellar MZR derived when we consider only young stars. We also find a strong relation between the local values of $\mu$* and Z* (the $\mu$ZR), betraying the influence of local factors in determining Z*. This shows that both local ($\mu$*-driven) and global (M*-driven) processes are important in determining the metallicity in galaxies. We find that the overall balance between local and global effects varies with the location within a galaxy. In disks, $\mu$* regulates Z*, producing a strong $\mu$ZR whose amplitude is modulated by M*. In spheroids it is M* who dominates the physics of star formation and chemical enrichment, with $\mu$* playing a minor, secondary role. These findings agree with our previous analysis of the star formation histories of CALIFA galaxies, which showed that mean stellar ages are mainly governed by surface density in galaxy disks and by total mass in spheroids.

ACTIVE GALACTIC NUCLEI EMISSION LINE DIAGNOSTICS AND THE MASS-METALLICITY RELATION UP TO REDSHIFTz∼ 2: THE IMPACT OF SELECTION EFFECTS AND EVOLUTION

Emission line diagnostic diagrams probing the ionization sources in galaxies, such as the Baldwin-Phillips-Terlevich (BPT) diagram, have been used extensively to distinguish AGN from purely star-forming galaxies. Yet, they remain poorly understood at higher redshifts. We shed light on this issue with an empirical approach based on a z~0 reference sample built from ~300,000 SDSS galaxies, from which we mimic selection effects due to typical emission line detection limits at higher redshift. We combine this low-redshift reference sample with a simple prescription for luminosity evolution of the global galaxy population to predict the loci of high-redshift galaxies on the BPT and Mass-Excitation (MEx) diagnostic diagrams. The predicted bivariate distributions agree remarkably well with direct observations of galaxies out to z~1.5, including the observed stellar mass-metallicity (MZ) relation evolution. As a result, we infer that high-redshift star-forming galaxies are consistent with having "normal" ISM properties out to z~1.5, after accounting for selection effects and line luminosity evolution. Namely, their optical line ratios and gas-phase metallicities are comparable to that of low-redshift galaxies with equivalent emission-line luminosities. In contrast, AGN narrow-line regions may show a shift toward lower metallicities at higher redshift. While a physical evolution of the ISM conditions is not ruled out for purely star-forming galaxies, and may be more important starting at z>2, we find that reliably quantifying this evolution is hindered by selections effects. The recipes provided here may serve as a basis for future studies toward this goal. Code to predict the loci of galaxies on the BPT and MEx diagnostic diagrams, and the MZ relation as a function of emission line luminosity limits, is made publicly available.

The mass–metallicity relation at z ∼ 1.4 revealed with Subaru/FMOS★

We present a stellar mass–metallicity relation at z ∼ 1.4 with an unprecedentedly large sample of ∼340 star-forming galaxies obtained with Fibre Multi-Object Spectrograph (FMOS) on the Subaru Telescope. We observed K-band selected galaxies at 1.2 ≤ zph ≤ 1.6 in the Subaru XMM–Newton Deep Survey/Ultra Deep Survey fields with M* ≥ 109.5 M⊙, and expected F(Hα) ≥ 5 × 10−17 erg s−1 cm−2. Among the observed ∼1200 targets, 343 objects show significant Hα emission lines. The gas-phase metallicity is obtained from [N ii] λ6584/Hα line ratio, after excluding possible active galactic nuclei. Due to the faintness of the [N ii] λ6584 lines, we apply the stacking analysis and derive the mass–metallicity relation at z ∼ 1.4. Our results are compared to past results at different redshifts in the literature. The mass–metallicity relation at z ∼ 1.4 is located between those at z ∼ 0.8 and z ∼ 2.2; it is found that the metallicity increases with decreasing redshift from z ∼ 3 to z ∼ 0 at fixed stellar mass. Thanks to the large size of the sample, we can study the dependence of the mass–metallicity relation on various galaxy physical properties. The average metallicity from the stacked spectra is close to the local Fundamental Metallicity Relation (FMR) in the higher metallicity part but ≳ 0.1 dex higher in metallicity than the FMR in the lower metallicity part. We find that galaxies with larger E(B − V), B − R and R − H colours tend to show higher metallicity by ∼0.05 dex at fixed stellar mass. We also find relatively clearer size dependence that objects with smaller half-light radius tend to show higher metallicity by ∼0.1 dex at fixed stellar mass, especially in the low-mass part.

THE UNIVERSAL STELLAR MASS-STELLAR METALLICITY RELATION FOR DWARF GALAXIES

We present spectroscopic metallicities of individual stars in seven gas-rich dwarf irregular galaxies (dIrrs), and we show that dIrrs obey the same mass-metallicity relation as the dwarf spheroidal (dSph) satellites of both the Milky Way and M31: Z_* ~ M_*^(0.30 +/- 0.02). The uniformity of the relation is in contradiction to previous estimates of metallicity based on photometry. This relationship is roughly continuous with the stellar mass-stellar metallicity relation for galaxies as massive as M_* = 10^12 M_sun. Although the average metallicities of dwarf galaxies depend only on stellar mass, the shapes of their metallicity distributions depend on galaxy type. The metallicity distributions of dIrrs resemble simple, leaky box chemical evolution models, whereas dSphs require an additional parameter, such as gas accretion, to explain the shapes of their metallicity distributions. Furthermore, the metallicity distributions of the more luminous dSphs have sharp, metal-rich cut-offs that are consistent with the sudden truncation of star formation due to ram pressure stripping.

THE CHEMICAL EVOLUTION OF STAR-FORMING GALAXIES OVER THE LAST 11 BILLION YEARS

We calculate the stellar mass-metallicity relation at five epochs ranging to z~2.3. We quantify evolution in the shape of the mass-metallicity relation as a function of redshift; the mass-metallicity relation flattens at late times. There is an empirical upper limit to the gas-phase oxygen abundance in star-forming galaxies that is independent of redshift. From examination of the mass-metallicity relation and its observed scatter we show that the flattening at late times is a consequence of evolution in the stellar mass where galaxies enrich to this empirical upper metallicity limit; there is also evolution in the fraction of galaxies at a fixed stellar mass that enrich to this limit. The stellar mass where metallicities begin to saturate is ~0.7 dex smaller in the local universe than it is at z~0.8.

The role of cold gas and environment on the stellar mass-metallicity relation of nearby galaxies

We investigate the relationship between stellar mass, metallicity and gas content for a magnitude- and volume-limited sample of 260 nearby late-type galaxies in different environments, from isolated galaxies to Virgo cluster members. We derive oxygen abundance estimates using new integrated, drift-scan optical spectroscopy and the base metallicity calibrations of Kewley & Ellison (2008). Combining these measurements with ultraviolet to near-infrared photometry and HI 21 cm line observations, we examine the relations between stellar mass, metallicity, gas mass fraction and star formation rate. We find that, at fixed stellar mass, galaxies with lower gas fractions typically also possess higher oxygen abundances. We also observe a relationship between gas fraction and metal content, whereby gas-poor galaxies are typically more metal-rich, and demonstrate that the removal of gas from the outskirts of spirals increases the observed average metallicity by approximately 0.1 dex. Although some cluster galaxies are gas-deficient objects, statistically the stellar-mass metallicity relation is nearly invariant to the environment, in agreement with recent studies. These results indicate that internal evolutionary processes, rather than environmental effects, play a key role in shaping the stellar mass-metallicity relation. In addition, we present metallicity estimates based on observations of 478 nearby galaxies.

An empirical prediction for stellar metallicity distributions in nearby galaxies

We combine star-formation histories derived from observations of high redshift galaxies with measurements of the z ∼ 0 relation between gas-phase metallicity, stellar mass and star-formation rate to make an explicit and completely empirical connection between near-field and distant galaxy observations. Our approach relies on two basic assumptions: (1) galaxies’ average paths through time in stellar mass versus star-formation rate space are represented by a family of smooth functions that are determined by the galaxies’ final stellar mass and (2) galaxies grow and become enriched with heavy elements such that they always evolve along the mass–metallicity–star-formation rate relation. By integrating over these paths, we can track the chemical evolution of stars in galaxies in a model-independent way, without the need for explicit assumptions about gas inflow, outflow or star-formation efficiency. Using this approach, we present predictions of stellar metallicity (i.e. O/H) distribution functions for present day star-forming galaxies of different stellar masses and the evolution of the α-element stellar metallicity–mass relation since z ∼ 1. The metallicity distribution functions are fairly well described as Gaussians, truncated at high metallicity, with power-law tails to low metallicity. We find that the stellar metallicity distribution for Milky Way mass galaxies is in reasonable agreement with observations for our Galaxy, and that the predicted stellar mass versus mean stellar metallicity relation at z = 0 agrees quite well with results derived from galaxy surveys. This validates the assumptions that are implicit in our simple approach. Upcoming observations will further test these assumptions and their range of validity, by measuring the mean stellar mass–metallicity relation up to z ∼ 1, and by measuring the stellar metallicity distributions over a range of galaxy masses.

The role of cold gas on the stellar mass - metallicity relation of nearby galaxies

AbstractWe investigate the relationship between stellar mass, metallicity and gas content for a magnitude- and volume-limited sample of 260 nearby late-type galaxies in different environments. Combining new oxygen abundance measurements with ultraviolet to near-infrared photometry and Hi 21 cm line observations, we observe the relationship between stellar mass and metallicity. We also find that, at fixed stellar mass, galaxies with lower gas fractions typically possess higher oxygen abundances. Gas-poor galaxies are typically more metal-rich. Our results indicate that internal evolutionary processes, rather than environmental effects, play a key role in shaping the stellar mass-metallicity relation.