Galactic chemical evolution

Galactic chemical evolution follows the production, transport, dilution, locking, and loss of chemical elements as stars form and die and gas enters or leaves a galaxy.

An alpha element is an element such as oxygen, neon or magnesium associated with nucleosynthesis involving helium nuclei and massive-star burning. Core-collapse supernovae promptly return substantial abundances of these elements. Their abundance relative to iron can distinguish a rapid enrichment history from one with substantial delayed Type Ia supernova enrichment.

Massive stars produce prompt enrichment through core-collapse supernovae, whereas a Type Ia supernova follows binary stellar evolution and contributes substantial iron after a distribution of delays. A stellar population formed rapidly before much delayed enrichment can have a larger alpha-element to iron ratio than a stellar population formed over a prolonged interval. Relative abundance also depends on the initial mass function, stellar yields, and gas exchange; this is a chemical clock with model assumptions.

With well-mixed gas-phase metallicity $Z$, zero mass return, star formation rate $\psi$, pristine galactic gas inflow $\nu \psi$, and constant stellar yield $y$, mass conservation and metal accounting give $\dot{M}=-(1-\nu )\psi$ and $M\dot{Z}=(y-\nu Z)\psi$. Thus, for $0<\nu <1$,The initial gas-phase metallicity $Z_0=0$ gives the usual pristine-infall expression. Neglecting returned bulk mass while retaining a nonzero metal yield is a trace-metal idealization.

For $\nu =1/2$ and initially pristine gas, $M=M_0(1-Z/(2y))$ and the formed stellar mass is $S=2(M_0-M)=M_{0}Z/y$. If a fixed initial mass function supplies $\kappa$ surviving tracers per unit formed mass, then $N(<Z)=\kappa M_{0}Z/y$ and $dN/dZ=\kappa M_0/y$ until the current gas abundance. Complete exhaustion gives a uniform normalized metallicity distribution function on $0<Z<2y$. This is uniform in linear metallicity, not in its logarithm.

The gas fraction of a galaxy is commonly $\mu =M_g/(M_g+M_s)$, using the gas and stellar masses remaining in the system. In a leaky-box model of galactic chemical evolution, it differs from the fraction $g=M_g/M_{g,0}$ of the initial gas reservoir. For an initially star-free box with outflow rate $\alpha {\dot{M}}_s$, $\mu =(1+\alpha )g/(1+\alpha g)$. These different fractions give different-looking metallicity formulae and effective yields.

Galactic metallicity measures the abundance of elements heavier than helium in a galaxy. Gas and stars retain different time-weighted records of enrichment.

Gas-phase metallicity is the heavy-element abundance of interstellar gas. In star-forming galaxies it is commonly inferred from nebular emission lines in H II regions, often through oxygen abundance $12+{\log }_{10}(\mathrm{O}/\mathrm{H})$.

Stellar metallicity is the heavy-element abundance locked in stars. Integrated-light measurements are luminosity weighted unless a stellar-population reconstruction converts them to a mass-weighted quantity.

In a well-mixed closed-box model of galactic chemical evolution, gas consumed into locked stellar mass supplies a natural cumulative mass coordinate $\delta$. The mean stellar metallicity is ${\overline{Z}}_{\ast }(\delta )={\delta }^{-1}{\int }_0^{\delta }{Z}_{\mathrm{birth}}(s)ds$, giving the displayed derivative relation. If the birth abundance increases linearly, $Z_{\mathrm{birth}}=b\delta$, its cumulative mean is half the current maximum. More generally the maximum of any nonnegative abundance distribution cannot lie below its mean. A fixed initial mass function and the same survival selection are required to identify this mass-weighted distribution with a star-count distribution.

The galaxy mass--metallicity relation is the empirical tendency for galaxies with greater galaxy stellar mass to have higher gas-phase metallicity and stellar metallicity. Mass-dependent metal loss, gas inflow, star-formation efficiency, and recycling all contribute.

The instantaneous recycling approximation treats short-lived stars as returning their ejecta immediately after formation while long-lived stars retain their mass indefinitely. It turns delayed chemical enrichment into local differential equations.

Under the instantaneous recycling approximation, forming mass $dM$ locks $ds=(1-\beta )dM$ into stars. If a fraction $f$ of returned material remains in the interstellar medium, the retained original metal mass is $\beta fZdM$, and retained newly made metal mass is $fyds$, where $y$ is the stellar yield per locked mass. Subtracting the metal mass $ZdM$ initially consumed gives the displayed equation. Pristine galactic gas inflow adds gas but no metals; enriched inflow adds its metal mass separately. The retained fraction applies to both old and newly made metals in this idealization.

Let $S=s/s_{\infty }$, let the interstellar medium have mass $g=s_{\infty }S(1-S)$, and retain fraction $f=S$ of returned material. For $0\le \beta <1$ put $q=(2-\beta )/(1-\beta )$. The integrating factor for the metal mass is $S^{1/(1-\beta )}/(1-S)$. A finite initial gas-phase metallicity removes the homogeneous solution and gives the displayed expression, equivalently $Z=y{\sum }_{n\ge 1}{S}^n/(q+n)$ for $0\le S<1$. The increasing initial gas mass requires pristine galactic gas inflow: gas accounting gives $d{M}_{\mathrm{in}}/ds=q(1-S)$. A fixed initial mass function makes the birth metallicity distribution function proportional to $ds/dZ$ when enrichment is monotone.

The stellar yield is the newly synthesized metal mass returned to gas per unit mass locked into long-lived stars, under a specified recycling convention and initial mass function.

The closed-box model evolves a fixed total baryonic mass with no inflow or outflow. Under instantaneous recycling and constant yield it predicts $Z=y\log (1/\mu )$ for gas fraction $\mu$.

A finite-metal-fraction variant of the closed-box model of galactic chemical evolution uses effective stellar yield $Y(1-Z_g)$, rather than constant yield. With pristine initial gas and $Y>0$, metal accounting gives $d{Z}_g/d\log (1/\mu )=Y(1-Z_g)$ and hence $Z_g=1-{\mu }^Y$. The mass-weighted mean stellar metallicity is $[Y-(Y+1)\mu +{\mu }^{Y+1}]/[(Y+1)(1-\mu )]$. At final gas fraction ${\mu }_f$, the metallicity distribution function is $(1-z{)}^{1/Y-1}/[Y(1-{\mu }_f)]$ on $0<z<1-{\mu }_f^Y$. For trace abundances, $Z_g\simeq -Y\log \mu$ and the usual closed-box metallicity distribution is recovered to leading order.

For a pristine closed-box model of galactic chemical evolution, $M_{g}dZ=Yd{M}_s$ and $d{M}_s=-d{M}_g$ imply $M_g=M_0{e}^{-Z/Y}$ independently of the star formation rate. The resulting metallicity distribution function is exponential on $0\le Z\le Z_f$. Its normalized density is $e^{-Z/Y}/[Y(1-e^{-Z_f/Y})]$. A fixed initial mass function and survival selection are needed to turn locked stellar mass into a stellar number distribution. Logarithmic abundance intervals have the additional factor $(\ln 10)Z$.

With gross star formation rate $\psi =k{M}_g$, constant returned mass fraction $F$, and the instantaneous recycling approximation, the net locked mass rate is $(1-F)\psi$. Mass conservation gives ${\dot{M}}_g=-(1-F)k{M}_g$. For constant stellar yield $Y$ per unit locked mass and initially pristine gas, metal conservation gives $M_g\dot{Z}=Y(1-F)\psi$, so gas-phase metallicity grows as $Z=Y(1-F)kt$. Gross star formation and permanent locking have different normalizations.

The leaky-box model adds a galactic outflow to the closed box. If the outflow has the ambient gas metallicity and constant mass-loading factor $\lambda$, its effective yield is reduced to $y/(1+\lambda )$ when returned fractions are absorbed into the definitions.

For initially pristine gas, zero bulk mass return, a constant stellar yield $y$ and a well-mixed galactic outflow of mass-loading factor $\nu$, gas conservation gives $d{M}_g=-(1+\nu )d{M}_s$ and metal accounting gives $M_{g}dZ=yd{M}_s$. Consequently $M_g=M_{g,0}{e}^{-Z/y_e}$ with $y_e=y/(1+\nu )$. Since $d{M}_s/dZ=M_g/y$, complete depletion gives a normalized mass-weighted metallicity distribution function $p(Z)=e^{-Z/y_e}/y_e$ and mean stellar metallicity $y_e$. A time-independent initial mass function makes the distribution of surviving tracers proportional to this distribution; a sample selected by light need not have the same weights.

The effective yield is the yield that a simple closed-box relation would need to reproduce an observed gas fraction and metallicity. Outflow and pristine inflow generally lower it below the nucleosynthetic yield.

The G-dwarf problem is the deficit of low-metallicity long-lived G dwarfs near the Sun relative to the prediction of a constant-yield closed-box model. Prolonged accretion of metal-poor gas is a standard ingredient in its resolution.

A metallicity distribution function gives the number or mass of stars per metallicity interval. Its interpretation depends on sample selection, stellar lifetimes, spatial migration, and whether metallicity means total heavy-element abundance or a particular abundance ratio.

Changing from linear stellar metallicity to its logarithm multiplies the metallicity distribution function by its Jacobian determinant. A flat distribution in linear metallicity rises exponentially in logarithmic metallicity. The commonly used iron abundance coordinate is equivalent to total-metallicity scaling only if relative element abundances are fixed.