Measure the number-density profile and vertical velocity dispersion of a tracer stellar population. The vertical Jeans equation estimates the restoring field from
The second term accounts for the tilt of the velocity ellipsoid; an exponential tracer with constant velocity dispersion and negligible tilt gives above the plane. The tracer need not supply all the gravitating mass. To infer the dynamical surface density from vertical stellar motions, integrate the axisymmetric Poisson equation for Newtonian gravity:
For a nearly flat galaxy rotation curve with a small radial correction, this reduces to . Estimate stellar surface density independently from counts, the initial mass function and faint objects, and add the gas surface density and stellar remnants. The residual
limits the allowed unseen material within the same height interval. A fitted dark matter halo contributes too: measurements at several heights help distinguish its extended density from a thin dark component in the galactic disk. Under a specified decomposition, the disc dark-matter fraction is the inferred dark-disc surface density divided by the total disc surface density, with uncertainty or an upper limit if the residual is consistent with zero. Selection effects, disequilibrium, the tilt term and uncertain baryonic counts must enter the error budget; an excess alone does not establish that all of it is in the galactic disk.
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- fraction from to , 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 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 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 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.
Metal absorption lines depend on the temperatures, ionization states and luminosity weights of the contributing stellar populations, as well as their stellar metallicity. A young stellar population dominated by hot luminous stars can have weak metal absorption lines even at substantial abundance, so weak lines alone do not establish low metallicity, an aspect of the age-metallicity degeneracy.
Low-mass main sequence stars live for much of the Milky Way's age, retaining representatives of many generations, while phase mixing disperses their birth associations around their orbits. In an approximately axisymmetric galactic disk, a sufficiently complete local census can therefore sample the mixed annular stellar population at the Sun's radius despite its small volume. This is a representative local-radius sample, not automatically an unbiased whole-disc census: different vertical scale heights, selection effects and radial migration require corrections when inferring integrated population fractions.
Older disc stellar populations generally have larger velocity dispersions, partly through accumulated galactic disk heating, so random-motion stresses supply more of their support against gravity. The radial Jeans equation then requires a smaller mean rotation than the circular speed, producing larger stellar asymmetric drift than in dynamically colder young stars.
Typical present elliptical galaxies are dominated by old stellar populations: their brightest short-lived stars have disappeared, leaving faint long-lived stars and stellar remnants, so luminosity per unit stellar mass is low. Ongoing star formation in many spiral galaxies and irregular galaxies supplies luminous young stars, giving a smaller stellar mass-to-light ratio in a fixed band; dark matter, the initial mass function and aperture must also be specified if the quoted ratio is dynamical.
Many open clusters are weakly bound and are dispersed by two-body relaxation, stellar evaporation and tidal heating, including encounters with molecular clouds in the galactic disk, so their present age histogram preferentially retains young survivors. More massive compact globular clusters can survive much longer, while the released low-mass stars remain as old field stellar populations: disruption removes clusters from the catalogue without removing their stars from the disc.
For this stellar population, let denote initial mass in units of the solar mass, and let denote present stellar age. Constant formation of equal numbers of stars per unit time makes the stellar age a uniform distribution on Gyr, so is uniform on . The normalized initial mass function has probability density function for , since . Consequently
These formulae have the stated age and mass domains; outside them the relevant cumulative fractions saturate at zero or one. The probability integral transform makes uniform on : gives . The time-independent initial mass function and constant number formation rate give independence of and . All fractions here count objects, including white dwarfs, using the stipulated stellar evolution law.
A red giant has
Thus the mass boundaries are and for positive . Equivalently, for the red giant region runs from to the age cap , and for it runs from to . There are no red giants with . Boundaries have zero probability and their endpoint convention does not affect the fractions. In the uniform square the red giant region is , with triangle vertices , and . The white dwarf region is the triangle .
Figure 1.
Mass-age regions and their uniform-coordinate images
. The right panel magnifies the evolved part of the unit square. All of the remaining region at larger X is main sequence. Details.
At fixed , define the conditional probabilities of a red giant, white dwarf and main sequence star by , and . Their interval widths are
Integration over the uniform distribution of age gives the individual-star fractions
The systems form a coeval binary population: the two binary star components have the same age. Their masses are independent, so has uniform probability density function on and their states have conditional independence given . Unconditional independence of their states would be incorrect: older systems make both evolved states more likely. The law of total probability now gives
The subtraction removes the double counting of systems containing two red giants.
Assume virial equilibrium, structurally homologous early-type galaxies, a fixed relation between the measured central velocity dispersion and the global random stellar motions, and a constant total mass-to-light ratio . The virial theorem gives , where the dimensionless structure factor is common to the sample. With denoting mean linear surface brightness inside the effective radius, . Thus
The usual observed fundamental plane of elliptical galaxies has a tilt relative to these exponents because mass-to-light ratio varies with stellar population, chemical composition, and dark matter fraction; structural and dynamical nonhomology can also change .
Printed convention. The source places the logarithm of a value expressed in in the empirical fit. This differs from the usual logarithm of linear surface brightness. If is the quoted magnitude surface brightness,
A coefficient on therefore becomes on , with an adjusted intercept; it does not become . For a concrete distinction, changing from 20 to 25 decreases linear surface brightness by a factor 100, so the linear-intensity fit changes by , whereas the literal printed logarithm changes it by .
The following part evaluates the printed numerical fit literally and labels that convention. If the fit was intended for linear surface brightness, its calibration must be specified; a conditional result for using the same coefficients with linear units is also given. No equality between these two conventions is assumed.
Evaluate the source's printed logarithm using the numerical magnitude value 21.0, with base-ten logarithms:
This is the literal printed-fit convention identified above. Convert the mean visual surface brightness to a linear luminosity per projected area. At 10 parsecs, one square arcsecond subtends . The definition of absolute magnitude consequently gives
The constant is the conversion between the angular and physical projected area. Since the effective radius encloses half the total luminosity,
For the stipulated stellar population of Solar-type stars, approximate the stellar mass-to-light ratio by . Summing 100 such galaxies then gives
This assumes negligible interstellar extinction and that the supplied surface brightness is an intrinsic or already corrected value; an observed cosmological surface brightness requires additional dimming and passband corrections before this conversion.
For comparison only, if the empirical coefficients and intercept were instead meant for , substituting in that fit would give , , and . Those are conditional alternative results, not an unannounced correction of the printed units.
Use an explicit Tully-Fisher-type scaling for comparable spiral populations and photometric band. Relative to the preceding galaxy,
The exponent and normalization are empirical and have scatter; one galaxy alone does not establish the relation. Different stellar populations, inclination errors and extinction change the photometric inference.
For an isolated final stellar population with isotropic one-dimensional velocity dispersion , kinetic energy is . The virial theorem gives and total energy . Collapse from a diffuse nearly unbound gas reservoir therefore requires removal of approximately of energy. The radiated fraction can differ if there are outflows, appreciable initial binding energy, or boundary stresses.
A plane-symmetric, self-gravitating stellar population with constant vertical velocity dispersion satisfies the vertical Jeans equation and the Poisson equation for Newtonian gravity. With , their solution is
Its surface density is , so . Its outer density tail is exponential with scale height ; its central profile is not exponential.
Spiral galaxy 2026-10-14
A spiral galaxy has a galactic disk with spiral structure, commonly cold interstellar gas and ongoing star formation as well as an older stellar population.
Star formation history 2026-10-25
The time-dependent rate at which a stellar population formed. Convolving it with an event delay-time distribution predicts the event rate; the time since star formation ended does not alone specify all stellar ages.
Stellar asymmetric drift Created 2026-10-14 Updated 2026-10-25
Stellar asymmetric drift is the lag of a pressure-supported stellar population's mean rotation behind the local circular speed. Neglecting the tilted velocity ellipsoid, radial Jeans equation balance gives
A declining density and larger random velocities generally increase the lag; gravitational balance is supplied partly by random-motion stresses.
Stellar halo 2026-10-14
A stellar halo is the extended, approximately spheroidal stellar population surrounding a galaxy. It differs from the total dark matter halo. If its outer luminosity density is , its projected light follows the projection of a spherical power-law emissivity and scales as away from breaks and outer truncation.
Stellar remnant 2026-10-14
A stellar remnant is the object left after a star finishes its principal nuclear-burning evolution, such as a white dwarf, neutron star or black hole. Stellar remnants retain mass while often contributing little optical light, affecting an old stellar population's mass-to-light ratio.
The High-Z Supernova Search Team analysis fitted multicolour light curves and interstellar extinction, and compared a second decline-rate method. Similar local and distant spectroscopy supported, but did not establish, population stability. Nearly colourless attenuation remained a qualified alternative.[ref]
The Supernova Cosmology Project compared low- and high-redshift colours, removed red outliers in its primary fit, and tested Malmquist bias with artificial-source recovery.[ref]
The inference can be organized as a chain of independently testable failure modes. A magnitude zero point offset between nearby and distant observations imitates a distance modulus offset; a common offset instead changes the intercept. An erroneous photometric passband response changes a K correction in a redshift-dependent way because different emitted wavelengths cross the detector band. Residual host-galaxy light changes the recovered radiative flux. Color excess constrains wavelength-dependent interstellar extinction only after intrinsic colour is modeled; grey attenuation is not bounded by colour alone. Changes in stellar population or metallicity could shift standardized absolute magnitude without large changes in the chosen shape parameter. These are logical degeneracies, so increasing the sample size alone cannot remove them.
A useful control strategy therefore compares subsets with different host environments, fits consistent rest-frame photometric passbands, and checks artificial-source recovery through the same image processing and selection. A Malmquist bias tends to select brighter detectable sources, but the relative bias between the nearby and distant samples is what matters. Peculiar velocity mainly perturbs the nearby redshift anchor; weak gravitational lensing changes individual brightnesses and the selection distribution. Neither can be replaced by an arbitrary extra scatter term when its mean or selection dependence is uncertain.
Vertical Jeans equation 2026-10-14
For a steady axisymmetric stellar population with no mean vertical motion, the vertical Jeans equation balances the gradient of the vertical random-motion stress and the radial gradient of the tilted velocity ellipsoid against the vertical Newtonian gravitational field. Neglecting the mixed moment gives . A tracer density need not equal the density sourcing the gravitational potential.