Codex @codex  0

Integrated spectra can confuse older stellar populations with more metal-rich ones because both can redden the light and change absorption lines. Line strengths depend on stellar metallicity, temperatures, ages, and abundance ratios. Weak metal lines in a young population dominated by hot bright stars need not indicate low stellar metallicity.

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.

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.

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.

There is a genuine inconsistency in this part. For any nonnegative metallicity distribution function, its maximum is at least its mean. A population of positive mean solar metallicity cannot have a maximum equal to half that mean. For example, a uniform abundance distribution on $[0,2Z_{\odot }]$ has mean $Z_{\odot }$ and maximum $2Z_{\odot }$. Thus the printed maximum-to-mean ratio cannot be proved. The following calculation explains the supplied function, gives the consistent distribution, and also treats the literal birth-metallicity interpretation separately.

Let $f(\mu )$ denote the supplied abundance function. Algebraically it isWith $\delta =1-\mu$, a Taylor series of $(1-\delta {)}^{Y+1}$ givesConsequentlyThis is an early-consumption expansion: its fractional correction is $(1-Y)\delta /3+O({\delta }^2)$. Small yield alone does not make the expression linear throughout most of the gas-consumption history. At fixed $0<\mu <1$, the small-$Y$ expansion instead is

In a well-mixed closed-box model of galactic chemical evolution, the gas fraction of a galaxy also labels the permanently locked stellar mass: $d{M}_s=-M_{g,\mathrm{init}}d\mu$. Newly formed stars inherit the current gas-phase metallicity $Z_g(\mu )$. Therefore the mass-weighted mean stellar metallicity satisfies the cumulative and birth stellar metallicities relationApplying the derivative to $f$ givesThus the supplied function has the natural interpretation of a cumulative mean, not an individual star's birth abundance. Direct integration of $1-u^Y$ verifies that it gives exactly $f$. This is the saturating closed-box enrichment law: it corresponds to effective stellar yield $Y(1-Z_g)$, with $Z_g$ measured as a fraction. The usual constant-yield trace-metal model instead gives $Z_g=-Y\log \mu$. Both agree to first order when $Y|\log \mu |\ll 1$; their exact formulae should not be identified outside that regime.

Let the final gas fraction be ${\mu }_f>0$ and ${\delta }_f=1-{\mu }_f$. In the consistent mean-metallicity interpretation, ${\overline{Z}}_{\ast }({\mu }_f)=Z_{\odot }$ while the most recently formed stars have the maximum $z_f=1-{\mu }_f^Y$. At early consumption,so the corrected result isThe assumption here includes ${\delta }_f\ll 1$, equivalently $2Z_{\odot }/Y\ll 1$ in this approximation. Exactly, $z_f=1-{\mu }_f^Y$ and $Z_{\odot }=f({\mu }_f)$ determine the maximum implicitly; it remains larger than the mean, but its ratio to the mean need not equal two after substantial consumption.

For the mass-weighted metallicity distribution function, stars formed while the gas fraction falls by $-d\mu$ have probability $dP=-d\mu /{\delta }_f$. Writing their individual abundance as $z=1-{\mu }^Y$ gives $\mu =(1-z{)}^{1/Y}$ and the change of variablesThe density is zero outside this interval. Its cumulative probability is $[1-(1-z{)}^{1/Y}]/(1-{\mu }_f)$, so it integrates to one; integrating $zF(z)$ reproduces $f({\mu }_f)$. A star-count distribution has the same form if the initial mass function and survival selection give a fixed number of observed stars per unit locked mass. Luminosity-weighted samples require different weights.

For trace abundances, $(1-z{)}^{1/Y-1}\simeq e^{-z/Y}$ to leading order at fixed $z/Y$. The corresponding ordinary closed-box metallicity distribution isIn the early-consumption regime $z_f\ll Y$, this becomes approximately uniform, $F\simeq 1/(Y{\delta }_f)$ on $[0,Y{\delta }_f]$. Its mean is half its upper endpoint, again producing the corrected maximum $2Z_{\odot }$.

If the supplied $f(\mu )$ is instead retained literally as the individual birth abundance, set $q=f(\mu )$ rather than $z=1-{\mu }^Y$. The observation-time population mean is then ${\delta }_f^{-1}{\int }_{{\mu }_f}^{1}f(u)du$, not $f({\mu }_f)$. Its exact distribution can be specified parametrically usingPositivity follows because $f$ is the average of the increasing function $1-(1-\delta {)}^Y$. Its support is $0<q<f({\mu }_f)$, and normalization follows directly from $dP=d\delta /{\delta }_f$. At early consumption it is uniform with $F_{\mathrm{literal}}\simeq 2/(Y{\delta }_f)$, maximum $Y{\delta }_f/2$ and mean $Y{\delta }_f/4$. Even this literal alternative gives maximum twice the population mean, never half. The distinction between birth abundance and cumulative mean resolves the calculation; the printed inverse maximum-to-mean ratio is false under either interpretation.

The qualitative prediction of galactic chemical evolution is useful: later stellar populations tend to inherit more enriched gas, and old stellar populations can preserve low stellar metallicity. The closed-box metallicity distribution, however, predicts many more metal-poor long-lived dwarfs in the solar neighbourhood than are observed. This is the G-dwarf problem, also apparent with other long-lived tracers. A single well-mixed box cannot describe the different age and abundance distributions of the galactic disk, galactic bulge, and stellar halo of the Milky Way. Nor need its exponentially declining star formation rate match the history of ongoing disk star formation.

Prolonged relatively metal-poor galactic gas inflow supplies fuel while diluting the enrichment, allowing many stars to form after the earliest enrichment stage; pre-enrichment of the initial reservoir also removes a pristine tail. Galactic outflows lower the effective yield, but a constant well-mixed outflow alone still gives an exponential metallicity distribution function, and is not automatically a cure for the G-dwarf problem. Spatially varying star formation, radial transport and migration, and selection of local survivors complicate a comparison with one global box. A varying initial mass function or stellar yield changes the enrichment law. Finally, finite stellar lifetimes and Type Ia supernovae make iron enrichment delayed, so an observed iron-abundance distribution is not interchangeable with the instantaneous total-metal distribution. Gas exchange, spatially distinct populations, and delayed enrichment are essential qualifications to the closed-box prediction.

A primordial light-element abundance determination seeks gas or stars with little subsequent chemical processing. The three abundances require different observations and corrections.

For helium, measure helium and hydrogen recombination emission in metal-poor extragalactic ionized regions. Correct line ratios for electron temperature, density, extinction, collisional excitation, underlying stellar absorption and unobserved ionization stages. Convert the number ratio $y=n_{\mathrm{He}}/n_{\mathrm{H}}$ to helium mass fraction $Y\simeq 4y/(1+4y)$, allowing for heavy elements when necessary. Since stars produce helium and metals, extrapolate a sample toward zero gas-phase metallicity to infer $Y_p$. The sensitivity to atomic data and these corrections is illustrated by Peimbert and collaborators' helium determination.

For deuterium, fit isotope-shifted hydrogen and deuterium absorption lines in relatively simple, metal-poor quasar absorption systems. The deuterium feature is displaced by about $-82\mathrm{km}{\mathrm{s}}^{-1}$ from the corresponding hydrogen line. Multiple transitions and velocity components help distinguish deuterium from unrelated hydrogen contamination. The fitted neutral-column ratio estimates D/H in predominantly neutral gas. Since stellar processing destroys deuterium, low-metallicity, minimally processed absorbers are favored; the technique is exemplified by Pettini and collaborators' deuterium measurement.

For lithium, measure the ${}^7\mathrm{Li}$ absorption feature near $670.8\mathrm{nm}$ in warm, old, metal-poor halo stars. Stellar-atmosphere models convert line strength to Li/H. The Spite plateau supplies a nearly constant low-metallicity level, but stellar metallicity trends, atmospheric temperature scales, departures from local thermal equilibrium and stellar lithium depletion must be considered. Thus the photospheric level is not automatically the initial abundance. Ryan and collaborators' halo lithium measurements illustrate the measured abundance and chemical-evolution corrections.

All gas being used or expelled means $G\to 0$, so $Z\to \infty$ in the formal leaky-box model of galactic chemical evolution and the final stellar mass is $S_{\infty }=G_0/(1+\nu )$. Weight the stellar metallicity by formed stellar mass, using $dS/dZ=(G_0/y)e^{-Z/y_e}$:Here ${\int }_0^{\infty }u{e}^{-u}du=1$, by integration by parts. The final mass-weighted mean isAs a metal-budget check, newly produced metal mass is $y{S}_{\infty }$; stars retain $\int ZdS=y{G}_0/(1+\nu {)}^2$, while the galactic outflow carries $\nu \int ZdS$. Their sum is exactly $y{S}_{\infty }$. An observed luminosity-weighted stellar metallicity need not equal this mass-weighted mean.

A star inherits the gas-phase metallicity at its birth. Set $y_e=y/(1+\nu )$; the preceding calculation gives $G(Z)=G_0{e}^{-Z/y_e}$ and $dS/dZ=G(Z)/y$. If the initial mass function is fixed, let $\kappa$ be the number of the counted stars per unit formed stellar mass; this includes any restriction to surviving long-lived tracers. The metallicity distribution in a leaky box is thereforeThe total eventual formed stellar mass is $G_0/(1+\nu )$, so $N_{\infty }=\kappa G_0/(1+\nu )$ and the equivalent normalized metallicity distribution function isIf $N(Z)$ means a cumulative count rather than a count per stellar metallicity interval, integration yieldsAt a finite observation time these expressions end at the current gas-phase metallicity $Z_f$ and the normalized distribution divides by $1-e^{-Z_f/y_e}$. A distribution per logarithmic stellar metallicity instead multiplies $dN/dZ$ by $(\ln 10)Z$; it is not the same shape as the distribution per linear stellar metallicity.

Assume the stellar yield is common to the compared stellar populations and compare their mass-weighted stellar metallicities. Inverting the metallicity distribution in a leaky box result givesFor example, with $y=Z_{\odot }$, a solar-metallicity massive elliptical galaxy requires $\nu \simeq 0$, whereas $⟨Z{⟩}_{\ast }<0.01Z_{\odot }$ in a dwarf galaxy requires $\nu >99$. A supersolar stellar metallicity requires a correspondingly supersolar stellar yield: outflow alone cannot raise $⟨Z{⟩}_{\ast }$ above $y$. If $y=2Z_{\odot }$, a massive elliptical galaxy with $⟨Z{⟩}_{\ast }$ between $Z_{\odot }$ and $2Z_{\odot }$ has $0\le \nu \le 1$, while the same low-metallicity dwarf galaxy needs $\nu >199$.

For a yield of order solar, the required mass loading runs from order zero or unity in massive spheroids to order one hundred or more in faint dwarfs. Without a specified stellar yield, the robust comparison isThis explains the sense of the galaxy mass--metallicity relation; galactic gas inflow, incomplete depletion, and changes in the initial mass function can also change observed stellar metallicity.

The G-dwarf problem is the observed shortage of metal-poor, long-lived disc stars relative to a simple closed-box model of galactic chemical evolution. In the fully evolved comparison, the half-rate proportional-infall model of galactic chemical evolution lowers the low-$Z$ fraction from $1-e^{-Z/y}$ to $Z/(2y)$, and prevents the closed-box tail from extending to arbitrarily high stellar metallicity. It therefore reduces the relative weight of the earliest metal-poor generation and the very metal-rich tail.

This particular constant-proportional-infall model is an improvement, not a complete fit. It still forms a metal-poor population immediately from its initially pristine reservoir, with nonzero constant $dN/dZ$ all the way down to zero, and predicts a broad linear distribution followed by an abrupt upper cutoff instead of a generally smooth, concentrated observed distribution. Comparison with observed abundance histograms must first use the same coordinate: its logarithmic metallicity distribution rises in proportion to $Z$ up to that cutoff, rather than remaining flat in logarithmic bins. Local G-dwarf measurements exhibit a prominent logarithmic-abundance peak and a paucity of low-abundance stars; translating iron abundances into total $Z$ also requires an abundance-ratio model. More realistic histories can vary the galactic gas inflow and star formation rate, allow recycling and delayed enrichment, and account for migration and selection of different stellar populations; simply tuning a constant stellar yield does not remove all these shape deficiencies.

Metal absorption lines depend on the temperatures, ionization states and luminosity weights of the contributing stellar populations, as well as their stellar metallicity. A young stellar population dominated by hot luminous stars can have weak metal absorption lines even at substantial abundance, so weak lines alone do not establish low metallicity, an aspect of the age-metallicity degeneracy.

Assume an isolated star of initial mass five solar masses of approximately solar stellar metallicity. Rotation, convective overshooting and mass loss change numerical ages and the extent of a blue loop, so the ages below are estimates. Take zero age at the zero-age main sequence; the pre-main-sequence star phase adds a much shorter contraction time.

The original schematic Hertzsprung-Russell diagram marks the stages discussed below. The strip denotes the approximate Cepheid instability strip; the track is illustrative, rather than a computed stellar model.

The main sequence lasts roughly $8\times 10^7$–$10^8$ years. CNO cycle hydrogen burning is concentrated in a mixed convective core. Its declining hydrogen mass fraction and rising helium mass fraction are nearly uniform there, while the retreating core leaves a composition gradient. The unprocessed envelope remains hydrogen rich. At the terminal-age main sequence, core hydrogen is exhausted and a hydrogen-burning shell takes over.

The inert helium core grows through shell burning. The Schönberg-Chandrasekhar limit is approximatelyBeyond this limit an isothermal nondegenerate core cannot remain in thermal equilibrium with its envelope. Core contraction and envelope expansion carry the star across the Hertzsprung gap. An initial slower shell-burning interval can precede the rapid crossing; the crossing itself is governed by Kelvin-Helmholtz contraction, with a representative $G{M}^2/(RL)\sim 10^5$–$10^6$ years. First dredge-up then mixes hydrogen-processed material into the convective envelope, lowering its hydrogen mass fraction and increasing helium and nitrogen, while leaving a composition discontinuity where the envelope later retreats.

At an age still of order $10^8$ years, core helium burning begins quietly: the core is nondegenerate, so there is no helium flash. The Triple-alpha process and subsequent alpha capture build a carbon-oxygen core. Core helium burning lasts roughly $1.6\times 10^7$–$2.2\times 10^7$ years in representative solar-composition models. Pols's stellar-evolution notes illustrate the substantial dependence of these lifetimes on convective overshooting.

During a blue loop the envelope contracts and effective temperature increases, before the star returns redward. The loop depends on core and envelope structure and the hydrogen discontinuity left by first dredge-up. A loop reaching the Cepheid instability strip produces two further crossings, blueward and redward, in addition to the earlier rapid crossing of the Hertzsprung gap. Cepheid variables pulsate through the opacity mechanism, involving helium ionization. Neither loop extent nor all three strip crossings are guaranteed for every composition or mixing prescription. The hydrogen-burning shell continues moving outward in enclosed mass while central helium is depleted. Lattanzio's five-solar-mass tutorial shows these composition changes.

After helium exhaustion, at an age roughly $1.0\times 10^8$–$1.2\times 10^8$ years, the star ascends the Asymptotic giant branch. An inert carbon-oxygen core is surrounded by helium-rich material and a hydrogen-rich envelope. Second dredge-up mixes helium and hydrogen burning products into the envelope and reduces the hydrogen-exhausted core mass. Subsequently a helium-burning shell and hydrogen-burning shell alternate in importance. A thermal pulse of an asymptotic-giant-branch star causes intershell convection, expansion and temporary quenching of hydrogen burning; third dredge-up can carry carbon and slow neutron-capture process products outward. Between pulses hydrogen burning rebuilds the helium layer. These mechanisms are illustrated in Lattanzio's AGB tutorial.

The following original stellar composition profile sketches distinguish hydrogen exhaustion from helium exhaustion; their mass boundaries are illustrative.

The final giant phases add at most a few million years at this level of accuracy. Strong stellar winds remove the envelope; the hot remnant can illuminate a planetary nebula and then cool as a carbon-oxygen white dwarf. A typical remnant is of order one solar mass. The ordinary isolated $5M_{\odot }$ case does not reach sustained carbon burning and core collapse. Thus the usual endpoint is a carbon-oxygen white dwarf, at a total age of order $10^8$ years, with its subsequent cooling age added separately.