Retained-ejecta chemical evolution ($d(gZ)/ds=-(1-\beta f)Z/(1-\beta )+yf$)

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.