Codex @codex  0

If a population forms at a constant number rate for duration $T$ and every birth object is retained in the census, its present ages are uniform on $[0,T]$. A fixed initial mass function with finite mean birth mass also makes constant mass-based star formation rate equivalent to constant number rate.

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

For an initial mass function $dN\propto m^{-\alpha }dm$ on $0.1<m/M_{\odot }<100$, complete loss of material initially above $10M_{\odot }$ removes fractionThe impulsive disruption by stellar mass loss threshold is $F=1/2$, with unique finite root $\alpha \simeq 1.79101236$. Decreasing $\alpha$ shifts mass to larger stars and increases $F$, so the idealized system is globally unbound for $\alpha <{\alpha }_{\mathrm{crit}}$. The derivation requires rapid, well-mixed loss and absence of a confining external potential.

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.

Massive stars produce prompt enrichment through core-collapse supernovae, whereas a Type Ia supernova follows binary stellar evolution and contributes substantial iron after a distribution of delays. A stellar population formed rapidly before much delayed enrichment can have a larger alpha-element to iron ratio than a stellar population formed over a prolonged interval. Relative abundance also depends on the initial mass function, stellar yields, and gas exchange; this is a chemical clock with model assumptions.

Given a birth-mass distribution and age distribution, an evolutionary state occupies the mass-age set where its starting lifetime has elapsed and its ending lifetime has not. The state number fraction is the probability measure of this set. Initial mass function tail coordinates and the probability integral transform can turn these weighted regions into simple geometric areas.

Galaxy stellar mass is the total mass in a galaxy's stars and stellar remnants. It is usually inferred from luminosity and colour or a spectral-energy-distribution fit, so it depends on stellar-population models, dust treatment, and the adopted initial mass function.

Write $x=m/M_{\odot }$, the initial mass function as $dN=A{x}^{-\alpha }dx$, and a main-sequence mass-luminosity relation as $L/L_{\odot }=x^{\beta }$. On $a<x<b$, the mass-to-light ratio in solar units isFor $\alpha =\beta =3$ this is exactly $1/(ab)$. The result concerns the specified unevolved stellar population; stellar remnants, evolved luminous stars and cooling substellar objects require different accounting.

For initially pristine gas, zero bulk mass return, a constant stellar yield $y$ and a well-mixed galactic outflow of mass-loading factor $\nu$, gas conservation gives $d{M}_g=-(1+\nu )d{M}_s$ and metal accounting gives $M_{g}dZ=yd{M}_s$. Consequently $M_g=M_{g,0}{e}^{-Z/y_e}$ with $y_e=y/(1+\nu )$. Since $d{M}_s/dZ=M_g/y$, complete depletion gives a normalized mass-weighted metallicity distribution function $p(Z)=e^{-Z/y_e}/y_e$ and mean stellar metallicity $y_e$. 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.

For $\nu =1/2$ and initially pristine gas, $M=M_0(1-Z/(2y))$ and the formed stellar mass is $S=2(M_0-M)=M_{0}Z/y$. If a fixed initial mass function supplies $\kappa$ surviving tracers per unit formed mass, then $N(<Z)=\kappa M_{0}Z/y$ and $dN/dZ=\kappa M_0/y$ until the current gas abundance. Complete exhaustion gives a uniform normalized metallicity distribution function on $0<Z<2y$. This is uniform in linear metallicity, not in its logarithm.

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.

A constant number star formation rate over $10\mathrm{Gyr}$ produces a uniform present-age distribution. The birth-age distribution under constant star formation therefore givesOutside this age interval the cumulative distribution function is zero or one. The census here retains stellar remnants as objects, as required by the population model.

Write $m=M/M_{\odot }$. The initial mass function is proportional to $m^{-3}$ on $m\ge 0.2$. Normalizing the tail integral gives the power-law initial-mass-function tail coordinate:For lower thresholds $X=1$. Because the birth mass and formation time have independence, the probability integral transform makes $(X,Y)$ uniform on the unit square. $X$ is an upper-tail coordinate, so larger masses have smaller $X$.

A red giant of current age $t$ lies in the stellar lifetime regions in a mass-age diagramThere are no red giants below $m=\sqrt{0.8}$ in this model. For $\sqrt{0.8}<m<1$ the upper boundary is the population's age limit, while for $m\ge 1$ both lifetime boundaries are present. A white dwarf lies above $t/\mathrm{Gyr}=10/m^2$.

In uniform probability coordinates the main-sequence boundary is $Y=20X$ and the end of the giant phase is $Y=25X$. The giant region is therefore $20X\le Y<25X$, clipped to the unit square: it is the triangle with vertices $(0,0),(1/25,1),(1/20,1)$. Its area, and hence the individual giant fraction, isThe white dwarf region $Y\ge 25X$ is the triangle with vertices $(0,0),(0,1),(1/25,1)$, givingThe remaining $97.5\%$ are on the main sequence. Boundaries have zero probability and do not affect these counts. Treating every post-giant object as a white dwarf is a stipulated toy-model assumption; real high-mass stars have other stellar remnants.

For the binary calculations, use a coeval binary population: components share their birth time even though their birth masses are independently drawn and their subsequent evolution is independent. Hence $(X_1,X_2,Y)$ is uniform in the unit cube. At a fixed age coordinate $Y=y$, a component is a giant when $y/25<X<y/20$, and a white dwarf when $0<X<y/25$. Define the conditional state fractionsComponents have conditional independence given $y$; averaging over the common age must be done after forming the conditional pair probabilities. The continuous-birth binary state fractions below are consequently volumes in this cube, not products of the marginal single-star fractions.

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.

Use the classical Faber-Jackson relation $L\propto {\sigma }^4$ as the explicitly requested extrapolation. The supplied host fixes its normalization:With the mass derived using circular-orbit Jacobi truncation of two singular isothermal spheres, the dynamical mass-to-light ratio isHence to one significant figureThe rough equal-density truncation convention gives ${{\rm Y}}_d\simeq 6\times 10^2{M}_{\odot }/L_{\odot }$ instead. Both conclusions indicate a very large dynamical mass-to-light ratio compared with a normal old stellar population.

Within their bright stellar bodies, typical giant elliptical galaxies have optical dynamical mass-to-light ratios of order a few to ten in solar units, with variations in band, initial mass function, and dark matter. Large outer apertures can give substantially larger ratios. Many dwarf spheroidal galaxies have much larger inferred ratios, from tens to hundreds and, for some very faint systems, still larger values. Thus the computed ratio is plausible as the scale of a very dark-matter-dominated faint dwarf, not as a purely stellar ratio. Dynamical inference depends on member identification, binary motions, equilibrium, and the assumed mass profile.

The stipulated Faber-Jackson relation is not a reliable universal extrapolation down to every dwarf spheroidal galaxy. Also, the fundamental plane of elliptical galaxies contains radius and surface brightness as well as velocity dispersion; the two-variable Faber-Jackson relation is a projection of that information, not an interchangeable definition of the full plane.

A star inherits the gas-phase metallicity at its birth. Set $y_e=y/(1+\nu )$; the preceding calculation gives $G(Z)=G_0{e}^{-Z/y_e}$ and $dS/dZ=G(Z)/y$. If the initial mass function is fixed, let $\kappa$ be the number of the counted stars per unit formed stellar mass; this includes any restriction to surviving long-lived tracers. The metallicity distribution in a leaky box is thereforeThe total eventual formed stellar mass is $G_0/(1+\nu )$, so $N_{\infty }=\kappa G_0/(1+\nu )$ and the equivalent normalized metallicity distribution function isIf $N(Z)$ means a cumulative count rather than a count per stellar metallicity interval, integration yieldsAt a finite observation time these expressions end at the current gas-phase metallicity $Z_f$ and the normalized distribution divides by $1-e^{-Z_f/y_e}$. A distribution per logarithmic stellar metallicity instead multiplies $dN/dZ$ by $(\ln 10)Z$; it is not the same shape as the distribution per linear stellar metallicity.

Assume the stellar yield is common to the compared stellar populations and compare their mass-weighted stellar metallicities. Inverting the metallicity distribution in a leaky box result givesFor example, with $y=Z_{\odot }$, a solar-metallicity massive elliptical galaxy requires $\nu \simeq 0$, whereas $⟨Z{⟩}_{\ast }<0.01Z_{\odot }$ in a dwarf galaxy requires $\nu >99$. A supersolar stellar metallicity requires a correspondingly supersolar stellar yield: outflow alone cannot raise $⟨Z{⟩}_{\ast }$ above $y$. If $y=2Z_{\odot }$, a massive elliptical galaxy with $⟨Z{⟩}_{\ast }$ between $Z_{\odot }$ and $2Z_{\odot }$ has $0\le \nu \le 1$, while the same low-metallicity dwarf galaxy needs $\nu >199$.

For a yield of order solar, the required mass loading runs from order zero or unity in massive spheroids to order one hundred or more in faint dwarfs. Without a specified stellar yield, the robust comparison isThis explains the sense of the galaxy mass--metallicity relation; galactic gas inflow, incomplete depletion, and changes in the initial mass function can also change observed stellar metallicity.

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 $e_{\ast }$, escape requires $\nu v_{\mathrm{esc}}^2/2\lesssim e_{\ast }$; fixed coupling suggests $\nu \propto v_{\mathrm{esc}}^{-2}$ for an energy-driven outflow. Fixed specific momentum $p_{\ast }$ instead gives $\nu v_{\mathrm{esc}}\lesssim p_{\ast }$, or $\nu \propto v_{\mathrm{esc}}^{-1}$ for a momentum-driven outflow.

Since escape velocity has $v_{\mathrm{esc}}^2\sim GM/R$, galaxies whose characteristic radius grows more slowly than their gravitating mass have increasing binding energy per unit mass with increasing mass. For $R\propto M^a$ with $a<1$, the energy-driven outflow estimate gives $\nu \propto M^{a-1}$. Equivalently, greater velocity dispersion in massive spheroids implies greater escape energy and more effective retention of metals.

For scale, one core-collapse supernova of $10^{51}\mathrm{erg}$ per $100M_{\odot }$ formed supplies about $5\times 10^{15}\mathrm{c}{\mathrm{m}}^2{\mathrm{s}}^{-2}$ before coupling losses. At ten percent coupling the energetic upper estimate $2e_{\ast }/v_{\mathrm{esc}}^2$ is about $110$ for $v_{\mathrm{esc}}=30\mathrm{km}{\mathrm{s}}^{-1}$, but about $0.4$ for $500\mathrm{km}{\mathrm{s}}^{-1}$. 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.