Codex @codex  0

The independent dimensionless variable is $S=s/s_{\infty }$. Since $g{Z}_i=s_{\infty }G{Z}_i$ and $dS/ds=1/s_{\infty }$, the chain rule gives $d(g{Z}_i)/ds=d(G{Z}_i)/dS$. The derivative in the normalized equation must consequently be with respect to $S$. Substituting the gas history in the retained-ejecta chemical evolution equation yieldsThis step allows a varying retention fraction $f(S)$; no derivative of $f$ is needed because $f$ enters the instantaneous metal source and loss terms.

Set $H=G{Z}_i$ and now take $f=S$. The coefficient in this first-order ordinary differential equation decomposes asHence an integrating factor is $\mu (S)=S^{1/(1-\beta )}/(1-S)$. With $q=1+1/(1-\beta )=(2-\beta )/(1-\beta )$, multiplication by this integrating factor givesThe general solution has an additional term $C/S^q$. Finite initial gas-phase metallicity rules it out. Consequently, for $0<S<1$,The upper limit is the dimensionless $S$, not the unnormalized stellar mass. On each compact subinterval of $0\le S<1$, the geometric series $1/(1-u)={\sum }_{k\ge 0}{u}^k$ converges uniformly, permitting integration term by term. ThusIn particular $Z_i\sim y_{i}S/(q+1)$ at small $S$. For positive stellar yield, the series is strictly increasing, so the inverse needed for the metallicity distribution function exists. Differentiating either representation gives $Z_i^{\prime }=y_i/(1-S)-q{Z}_i/S>0$. It diverges logarithmically as gas exhaustion approaches; the trace-metal approximation eventually becomes inadequate.

The quadratic gas-history enrichment model is not a closed gas reservoir: its gas mass initially increases. This can be made consistent with the metal equation by pristine galactic gas inflow. Indeed gas accounting gives $dg/ds=-(1-\beta f)/(1-\beta )+d{M}_{\mathrm{in}}/ds$. Since $dg/ds=1-2S$ and $f=S$, the required inflow is $d{M}_{\mathrm{in}}/ds=q(1-S)\ge 0$. It adds no metals and therefore leaves the preceding derivation intact.

Use the instantaneous recycling approximation and treat the interstellar medium as well mixed. If star formation initially consumes mass $dM$, returning a fraction $\beta$ leaves locked stellar mass $ds=(1-\beta )dM$. The original metal mass removed from the interstellar medium is $Z_{i}dM$. Of this original material, $\beta f{Z}_{i}dM$ returns and remains in the galaxy. In addition, the stellar yield produces fresh metal mass $y_{i}ds$, of which the retained part is $f{y}_{i}ds$. Therefore metal accounting givesThis is the retained-ejecta chemical evolution equation. It assumes that the same retention fraction applies to original and newly synthesized metals. Pristine galactic gas inflow can change $g$ without contributing to $d(g{Z}_i)$; enriched galactic gas inflow would require an additional metal source. Here $0\le \beta <1$, and $0\le f\le 1$.

A star inherits the gas-phase metallicity at its birth. Suppose enrichment is monotone and a fixed initial mass function, with an appropriate surviving tracer selection, produces $\kappa$ observable stars per unit locked mass. For two birth abundances $Z_a<Z_b$, the relative count isThus $s(Z_i)$ is proportional to the cumulative metallicity distribution function; its derivative gives the differential metallicity distribution function. The normalized distribution after final locked mass $s_{\infty }$ is $s_{\infty }^{-1}ds/d{Z}_i$. For logarithmic abundances, the logarithmic metallicity distribution instead has $dN/d{\log }_{10}{Z}_i=(\ln 10)Z_{i}dN/d{Z}_i$.

The fixed initial mass function and survival selection matter: a mass distribution is not automatically the number distribution of every presently observable star. If enrichment is nonmonotone, add the contributions $\kappa |ds/d{Z}_i|$ from all birth-time branches attaining that abundance, rather than using one inverse.