Codex @codex  0

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.

A newly formed stellar population inherits the current gas-phase metallicity. Eliminate time between locked-mass growth and enrichment:For $Y>0$, integration from pristine gas gives $M_g=M_0{e}^{-Z/Y}$ and hence the closed-box metallicity distributionTo obtain numbers rather than mass, specify a fixed initial mass function and count a long-lived tracer population. Let $\eta$ be its number of surviving stars per unit locked stellar mass. At the final gas-phase metallicity $Z_f$,with zero density outside that interval. The cumulative count is $N(<Z)=\eta M_0(1-e^{-Z/Y})$ within the interval, and the normalized metallicity distribution function isThese equations do not use $\psi (t)$: any positive star formation rate gives the same distribution under the same chemical assumptions. If the tracer number is instead normalized per gross formed mass, replace $\eta$ by that normalization divided by $1-F$. Short-lived stars require an age-dependent survival factor and generally do not follow this number distribution. In a logarithmic metallicity distribution, $dN/d{\log }_{10}Z=(\ln 10)ZdN/dZ$; its interior peak is at $Z=Y$, whereas the distribution per linear interval decreases from $Z=0$.

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 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.