Codex @codex  0

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

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.

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.

Let $Z$ be the well-mixed gas-phase metallicity, measured as a metal mass fraction. Removing gas into stars takes metals at rate $Z\psi$; prompt return restores existing metals at rate $ZF\psi$, together with newly synthesized metals at rate $Y(1-F)\psi$. Thus the closed-box model of galactic chemical evolution has metal budgetExpand the left side and use ${\dot{M}}_g=-(1-F)\psi$. The terms proportional to existing gas-phase metallicity cancel, leavingSince $\psi =k{M}_g$, initially pristine gas obeysThe stellar yield here is per permanently locked mass, so inserting another factor $1-F$ into $Y$ would double-count recycling. The approximation treats metals as a trace constituent and ignores finite stellar-return delays. Apply it while $Y\lambda t\ll 1$; the formal unlimited linear growth is not a physical prediction for almost entirely metallic gas.

Divide metal enrichment by gas depletion in the closed-box model of galactic chemical evolution:The logarithm appears on integrating from $M_0$ and $Z=0$:With the gas fraction of a galaxy $\mu =M_g/M_0$, this isIt gives $Z=0$ at $\mu =1$ and increasing gas-phase metallicity as gas is consumed, independently of the time dependence of the star formation rate. For this closed initially star-free system, $M_0=M_g+M_s$, so this initial-reservoir fraction equals the present gas-to-total baryonic fraction. That equality fails in an open system.

The five printed subparts are separated below. All use one well-mixed closed-box model of galactic chemical evolution with a stellar yield defined per permanently locked mass; the gross star formation rate must be distinguished from that locked-mass rate.

Let $G$ be the remaining gas mass, $S$ the formed stellar mass, and $G_0$ the initial gas mass; initially $S=0$ and the gas-phase metallicity is zero. Take the stellar yield $y$ and the mass-loading factor $\nu$ to be constant, with homogeneous mixing and no galactic gas inflow. With the specified zero bulk return, mass conservation givesNew metals add $ydS$, while star formation and the galactic outflow remove metals at the existing gas-phase metallicity $Z$. ThusThe leaky-box model of galactic chemical evolution therefore givesFor the usual gas fraction of a galaxy, $\mu =G/(G+S)$, the denominator is the baryonic mass still in the system, which decreases under the galactic outflow. Since $G_0=G+(1+\nu )S$, we have $G_0/G=1+(1+\nu )(1-\mu )/\mu$. In the remaining-mass convention,If instead $\mu$ denotes $G/G_0$ relative to the fixed initial reservoir, the same calculation gives $\boxed{Z=y\ln (1/\mu )/(1+\nu )}$. Stating the denominator matters: these two gas fraction of a galaxy conventions coincide for the closed-box model of galactic chemical evolution, recovered at $\nu =0$, but differ when gas escapes. Neglecting bulk return while retaining a nonzero stellar yield is the intended trace-metal approximation.

The closed-box model of galactic chemical evolution has $M=M_0{e}^{-Z/y}$ and $S=M_0-M$, giving $dN/dZ=(\kappa M_0/y)e^{-Z/y}$. To compare shapes without confusing different total stellar counts, the figure normalizes both metallicity distribution functions to unit integral after complete gas exhaustion and uses $x=Z/y$:The right panel also shows the corresponding cumulative fractions, $1-e^{-x}$ and $\min (x/2,1)$.

Pristine galactic gas inflow replenishes the fuel and dilutes freshly produced metals, so enrichment per unit newly formed stellar mass is slower and asymptotes to $2y$. More stars are formed at later abundances; the proportional-infall model of galactic chemical evolution makes twice as many stars as the closed box for the same initial gas mass in this limit. Consequently the normalized low-abundance fraction is smaller: for $Z\ll y$, it is $Z/(2y)$ rather than approximately $Z/y$. Their unnormalized $dN/dZ$ actually agree at $Z=0$ for the same $M_0$, which is why stating the normalization matters. At a finite stopping abundance $Z_f$, use $p_Z=1/Z_f$ for infall and $p_Z=e^{-Z/y}/[y(1-e^{-Z_f/y})]$ for the closed box, each truncated at its own chosen stopping point.

Let $\psi (t)$ be the star formation rate, $M$ the remaining gas mass, and $S$ the mass put into stars. With zero returned bulk gas, pristine galactic gas inflow, and no galactic outflow, mass conservation givesA well-mixed gas reservoir loses existing metals at rate $Z\psi$ and gains newly synthesized metals at rate $y\psi$; the incoming gas brings none. ThusFor $0<\nu <1$, eliminate time and integrate:The proportional-infall model of galactic chemical evolution therefore gives, for initially pristine gas,The displayed form requires $Z_0=0$; more generally $Z=y/\nu +(Z_0-y/\nu )(M/M_0{)}^{\nu /(1-\nu )}$. Also $S=(M_0-M)/(1-\nu )$ for an initially star-free reservoir, so $M/M_0$ is not the gas fraction of a galaxy $M/(M+S)$. As a check, the $\nu \to 0$ limit is the closed-box model of galactic chemical evolution relation $Z=y\log (M_0/M)$; at $\nu =1$ the gas mass is constant and one must instead solve $\dot{Z}=(y-Z)\psi /M_0$. Retaining a nonzero stellar yield while setting bulk mass return to zero is the prescribed trace-metal approximation, not a literal claim that stellar ejecta have no mass.

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.

Use the instantaneous recycling approximation with a well-mixed gas reservoir, constant stellar yield $p$, no gas inflow, and an irreversible galactic outflow. Let $M_s$ be the mass locked into long-lived stars and remnants after prompt recycling, and put $\psi (t)={\dot{M}}_s(t)\ge 0$. The yield $p$ is the newly synthesized mass of heavy elements returned to the gas per unit increase of $M_s$. Thus $M_h$ in $Z=M_h/M_g$ denotes the heavy-element mass in the gas, excluding metals already locked into stars.

Assume $M_s(0)=0$, $M_g(0)=M_0$, and $Z(0)=Z_0$; initially pristine gas means $Z_0=0$. Take the outflow rate to be $\alpha \psi$ with constant $\alpha \ge 0$, and assume its metallicity equals the current gas-phase metallicity $Z$. This is the well-mixed leaky-box model of galactic chemical evolution. The parameter $\alpha$ is a mass-loading factor relative to the net rate of locking mass into stars. If a gross star formation rate $\Psi$ and a prompt returned fraction $\mathcal{R}$ are used instead, $\psi =(1-\mathcal{R})\Psi$ and a wind rate $\eta \Psi$ corresponds to $\alpha =\eta /(1-\mathcal{R})$.

Conservation of mass givesand the baryonic mass remaining in the box is $M_g+M_s=M_0-\alpha M_s$. Newly made metals enter the gas at rate $p\psi$, while pre-existing metals are locked into stars at rate $Z\psi$ and leave in the wind at rate $\alpha Z\psi$. HenceUsing $M_h=Z{M}_g$ and the product rule,The cancellation expresses the fact that a well-mixed wind removes gas and metals in the same proportion, while fresh stellar production raises the abundance in the remaining gas.

Divide this equation by ${\dot{M}}_g=-(1+\alpha )\psi$ along the evolving reservoir, or integrate in time through intervals with no star formation. The resulting gas-phase metallicity isEquivalently, its fully explicit dependence on a prescribed net star formation rate isvalid while $M_g(t)>0$. With pre-existing stars, replace $M_s(t)$ in the gas-mass relation by $M_s(t)-M_s(0)$. Without specifying $\psi (t)$ or a gas-consumption law, the model fixes $Z$ as a function of gas mass but does not determine a unique function of time.

For the closed-box model of galactic chemical evolution, set $\alpha =0$. Then $M_g=M_0-M_s$ andAt the same remaining fraction $g=M_g/M_0$ of the initial gas reservoir, the leaky box has the logarithmic coefficient $p/(1+\alpha )$, lower than the true stellar yield by the wind factor. This coefficient is often called an effective yield, but the definition of the gas fraction must be specified.

In particular, the observable gas fraction of a galaxy is often $\mu =M_g/(M_g+M_s)$, using the mass still present, rather than $g$. In our initially star-free model,ThereforeFor $0<\mu <1$ and $\alpha >0$, the leaky-box enrichment is smaller at the same $\mu$: $\ln [1+\alpha (1-\mu )]\le \alpha (1-\mu )<\alpha \ln (1/\mu )$. If effective yield is defined observationally as $(Z-Z_0)/\ln (1/\mu )$, it is consequently not simply the constant $p/(1+\alpha )$ at all $\mu$.

As an explicit time example, choose a linear gas-consumption law $\psi =ϵM_g$, with constant $ϵ>0$. The mass and metallicity equations giveThe closed box with the same $ϵ$ has the same metallicity at a given time, while the leaky box exhausts its gas faster. Thus a comparison at fixed time depends on the chosen star formation rate, even though the comparison at fixed gas fraction above is unambiguous.

Finally, the assumption about wind metallicity is essential. If the expelled gas has abundance $Z_w$ instead, its metal-loss rate is $\alpha Z_w\psi$, and the same bookkeeping yieldsA wind enriched relative to the ambient gas lowers the enrichment rate; the logarithmic solution derived above applies to $Z_w=Z$.

The closed-box model of galactic chemical evolution assumes a fixed total baryonic mass, no gas inflow or outflow, instantaneous and homogeneous mixing, a constant stellar yield, a fixed initial mass function, and usually the instantaneous recycling approximation. It predicts $Z=y\log (1/\mu )$ for gas fraction $\mu$, but real galaxies violate these assumptions.

Evidence for galactic outflows includes blueshifted absorption, broad or split emission lines, extraplanar ionized and molecular gas, X-ray bubbles, and metal-enriched circumgalactic material. The low baryon fractions and low effective yields of dwarf galaxies, together with the galaxy mass--metallicity relation, provide indirect evidence for preferential gas and metal loss. Sustained star formation for longer than a closed reservoir's depletion time, low-metallicity high-velocity clouds, metallicity dilution during starbursts, the G-dwarf problem, and circumgalactic or intergalactic accretion signatures imply continuing galactic gas inflow.

Let $f(t)=A{e}^{-t/\tau }$ be the inflow rate, let the outflow rate be $\lambda \psi$, and let $R$ be the promptly returned fraction. The total baryonic, gas, and metal masses obeyThe first metal term is newly synthesized material; the second locks existing metals into long-lived stars and removes them in an ambient-composition wind. Expanding the final derivative and substituting ${\dot{M}}_g$ cancels those common terms, leaving

The observed Kennicutt–Schmidt law relates star-formation and gas surface densities approximately by ${\Sigma }_{\mathrm{SFR}}\propto {\Sigma }_g^{1.4}$, with a nearly linear relation to molecular gas over many resolved regimes. For the simple integrated model requested here, take a constant depletion coefficient $S$ so thatThen, with $\alpha =S(1-R+\lambda )$,The integrating factor $e^{\alpha t}$ gives, for $\alpha \ne 1/\tau$,At resonance, $\alpha =1/\tau$, the continuous limit isThe first term is depletion of the initial reservoir; the second is gas supplied by the exponentially declining inflow and subsequently consumed or expelled.

The Kennicutt–Schmidt law is the empirical relationfor disk-averaged total gas, with a nearly linear molecular-gas relation in many resolved observations. Atomic gas is mapped through the H I 21-cm line, molecular gas mainly through carbon-monoxide line emission and a CO-to-${\mathrm{H}}_2$ conversion factor, and star formation through combinations of ultraviolet continuum, H-alpha recombination emission, and infrared dust emission. Inclination, dust attenuation, the initial mass function, tracer lifetimes, and conversion factors must be treated consistently.

The gas-depletion time $M_g/{\dot{M}}_{\ast }$ is typically of order a gigayear, whereas a giant molecular cloud has a dynamical or free-fall time of order a few megayears. Star formation is therefore inefficient per collapse time, commonly at the percent level, rather than converting an entire cloud in one free fall.

For the first closed-box model of galactic chemical evolution, neglect returned mass or absorb it into the definitions. Thenand henceThus a region reaching $Z_{\odot }=0.014$ after an enrichment time $t_{\odot }$ with $y_Z=0.006$ hasabout $4.3\mathrm{Gyr}$ for a fiducial $t_{\odot }\simeq 10\mathrm{Gyr}$ enrichment age of the Galactic disk.

The second prescription implies ${\dot{M}}_{\ast }=M_g^2/(M_{g0}{\tilde{\tau }}_{\ast })$. Solving the gas-consumption equation givesTherefore

For long-lived stars, $dN$ is proportional to $d{M}_{\ast }=-d{M}_g$. In the exponential model,In the second model, $1+t/{\tilde{\tau }}_{\ast }=e^{Z/y_Z}$; multiplying its star-formation rate by $dt/dZ$ gives exactly the same result:The metallicity distribution is fixed by the closed-box relation $M_g/M_{g0}=e^{-Z/y_Z}$ and is independent of the star-formation history. Merely changing the time law therefore does not cure the G-dwarf problem; gas inflow, outflow, variable yields, or selection effects must alter the closed-box assumptions.

The galaxy mass--metallicity relation is the observed tendency for more massive galaxies to have larger gas-phase metallicity and stellar metallicity. Gas metallicity is commonly inferred from nebular emission-line ratios in star-forming H II regions, often quoted as $12+{\log }_{10}(\mathrm{O}/\mathrm{H})$ and measured within a finite spectroscopic aperture. Stellar metallicity comes from stellar absorption features or population-synthesis fits and is luminosity weighted unless the analysis explicitly reconstructs a mass-weighted distribution. Galaxy stellar mass is inferred from photometry or a spectral-energy-distribution fit and depends on the adopted initial mass function. Radial metallicity gradients, dust, line calibration, and aperture selection must consequently be matched before samples are compared.

The usual physical explanation is that a shallow potential well lets a low-mass galaxy lose a larger fraction of newly synthesized metals in galactic outflows. The closed-box model of galactic chemical evolution is therefore replaced by a leaky-box model of galactic chemical evolution withwhere $\lambda$ is the mass-loading factor. Under the instantaneous recycling approximation, let $y_z$ be the stellar yield, absorb the returned mass fraction into the definitions, and suppose the escaping gas has the current gas metallicity $Z_g$. Then mass conservation and metal conservation areSubstitution of the first equation into the second cancels the terms that merely transfer pre-existing metals and leavesBecause the mass of metals locked into stars obeys $d{M}_{z,\ast }=Z_{g}d{M}_{\ast }$, integration gives $M_{z,\ast }=M_{\ast }{Z}_{\ast }=\int Z_{g}d{M}_{\ast }$. The total newly made metal mass is partitioned between present gas, stars, and the outflow:Writing the gas-to-stellar mass ratio as $r_g=M_g/M_{\ast }$ therefore producesThus simultaneous gas and stellar metallicities, together with the gas fraction and an assumed nucleosynthetic yield, estimate the integrated mass loading. The corresponding effective yield is $y_{\mathrm{eff}}=y_z/(1+\lambda )$, and the leaky box has

The G-dwarf problem is that the local Milky Way disk contains far fewer low-metallicity long-lived G dwarfs than the constant-yield closed-box metallicity distribution predicts. A yield that rises with metallicity may initially sound promising because enrichment would accelerate after the first generations. In fact it worsens the problem. In a closed box, $d{M}_{\ast }=-d{M}_g$ andFor $y_z=kZ$ and a nonzero initial metallicity $Z_0$ at gas mass $M_{g,0}$,Hence the cumulative mass of stars born below metallicity $Z$ isIf the system begins at $Z_0=0$, the assumed yield also vanishes and enrichment never starts. For $Z_0>0$, the metallicity distribution function haswhich puts still more stellar mass near the low-metallicity floor. Metal-poor gas inflow, pre-enrichment, and selective outflow are therefore more plausible ingredients in resolving the G-dwarf problem.

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.