If stellar feedback supplies energy per unit formed stellar mass, an energy-driven outflow requires . Fixed energy coupling therefore suggests mass-loading factor . A fixed momentum supply instead gives and momentum-driven outflow scaling . These are launch and escape estimates, not universal wind laws; subsequent entrainment, cooling and fallback change the escaping mass.
Leaky-box model of galactic chemical evolution Created 2026-09-28 Updated 2026-10-14
The leaky-box model adds a galactic outflow to the closed box. If the outflow has the ambient gas metallicity and constant mass-loading factor , its effective yield is reduced to when returned fractions are absorbed into the definitions.
For initially pristine gas, zero bulk mass return, a constant stellar yield and a well-mixed galactic outflow of mass-loading factor , gas conservation gives and metal accounting gives . Consequently with . Since , complete depletion gives a normalized mass-weighted metallicity distribution function and mean stellar metallicity . 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.
Let be the remaining gas mass, the formed stellar mass, and the initial gas mass; initially and the gas-phase metallicity is zero. Take the stellar yield and the mass-loading factor to be constant, with homogeneous mixing and no galactic gas inflow. With the specified zero bulk return, mass conservation gives
New metals add , while star formation and the galactic outflow remove metals at the existing gas-phase metallicity . Thus
The leaky-box model of galactic chemical evolution therefore gives
For the usual gas fraction of a galaxy, , the denominator is the baryonic mass still in the system, which decreases under the galactic outflow. Since , we have . In the remaining-mass convention,
If instead denotes relative to the fixed initial reservoir, the same calculation gives . Stating the denominator matters: these two gas fraction of a galaxy conventions coincide for the closed-box model of galactic chemical evolution, recovered at , but differ when gas escapes. Neglecting bulk return while retaining a nonzero stellar yield is the intended trace-metal approximation.
The escape-speed scaling of wind mass loading makes the proposed trend dynamically plausible: the shallower Newtonian gravitational potential of a dwarf galaxy lets stellar feedback expel more gas per unit stellar mass. If coupled feedback supplies specific energy , escape requires ; fixed coupling suggests for an energy-driven outflow. Fixed specific momentum instead gives , or for a momentum-driven outflow.
Since escape velocity has , galaxies whose characteristic radius grows more slowly than their gravitating mass have increasing binding energy per unit mass with increasing mass. For with , the energy-driven outflow estimate gives . Equivalently, greater velocity dispersion in massive spheroids implies greater escape energy and more effective retention of metals.
For scale, one core-collapse supernova of per formed supplies about before coupling losses. At ten percent coupling the energetic upper estimate is about for , but about for . Efficient loss from shallow potentials and retention in deep potentials can span the required range. These are illustrative energetics rather than a unique prediction: cooling, entrainment, fallback and the actual initial mass function determine the escaping mass-loading factor.
Use the instantaneous recycling approximation with a well-mixed gas reservoir, constant stellar yield , no gas inflow, and an irreversible galactic outflow. Let be the mass locked into long-lived stars and remnants after prompt recycling, and put . The yield is the newly synthesized mass of heavy elements returned to the gas per unit increase of . Thus in denotes the heavy-element mass in the gas, excluding metals already locked into stars.
Assume , , and ; initially pristine gas means . Take the outflow rate to be with constant , and assume its metallicity equals the current gas-phase metallicity . This is the well-mixed leaky-box model of galactic chemical evolution. The parameter is a mass-loading factor relative to the net rate of locking mass into stars. If a gross star formation rate and a prompt returned fraction are used instead, and a wind rate corresponds to .
Conservation of mass gives
and the baryonic mass remaining in the box is . Newly made metals enter the gas at rate , while pre-existing metals are locked into stars at rate and leave in the wind at rate . Hence
Using 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 along the evolving reservoir, or integrate in time through intervals with no star formation. The resulting gas-phase metallicity is
Equivalently, its fully explicit dependence on a prescribed net star formation rate is
valid while . With pre-existing stars, replace in the gas-mass relation by . Without specifying or a gas-consumption law, the model fixes 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 . Then and
At the same remaining fraction of the initial gas reservoir, the leaky box has the logarithmic coefficient , 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 , using the mass still present, rather than . In our initially star-free model,
Therefore
For and , the leaky-box enrichment is smaller at the same : . If effective yield is defined observationally as , it is consequently not simply the constant at all .
As an explicit time example, choose a linear gas-consumption law , with constant . The mass and metallicity equations give
The 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 instead, its metal-loss rate is , and the same bookkeeping yields
A wind enriched relative to the ambient gas lowers the enrichment rate; the logarithmic solution derived above applies to .
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 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 with
where is the mass-loading factor. Under the instantaneous recycling approximation, let be the stellar yield, absorb the returned mass fraction into the definitions, and suppose the escaping gas has the current gas metallicity . Then mass conservation and metal conservation are
Substitution of the first equation into the second cancels the terms that merely transfer pre-existing metals and leaves
Because the mass of metals locked into stars obeys , integration gives . The total newly made metal mass is partitioned between present gas, stars, and the outflow:
Writing the gas-to-stellar mass ratio as therefore produces
Thus 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 , 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, and
For and a nonzero initial metallicity at gas mass ,
Hence the cumulative mass of stars born below metallicity is
If the system begins at , the assumed yield also vanishes and enrichment never starts. For , the metallicity distribution function has
which 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.