For pressureless linear cosmological density perturbations in an Einstein-de Sitter universe, and . The growth equation is an Euler differential equation whose indicial equation is . Hence the independent modes are and . The first is the growing linear growth factor, while the second is a decaying transient.
Use a monotone linear growth factor as time, , and with and . The cosmological Poisson equation becomes ; the continuity equation becomes . The material derivative is . Substitution into the pressureless cosmological Euler equation, followed by the growth equation and , proves the displayed motion equation, with . The density derivative in the conservative continuity equation is Eulerian, not material.
Linear growth equation Created 2026-09-29 Updated 2026-10-07
For pressureless matter on subhorizon scales, the linear growth factor obeys
The second term is dilution by cosmic expansion, while the final term drives Jeans instability.
Matter-spectrum response to extra dark energy Created 2026-10-07 Updated 2026-10-14
With physical matter and radiation densities, primordial amplitudes and spatial flatness fixed, extra smooth late dark energy increases but keeps fixed. Provided the dark energy is negligible at equality, the early cosmological transfer function and equality wavenumber in physical units are unchanged. Earlier accelerated expansion suppresses the linear growth factor, lowering present dimensional matter power by its squared growth ratio. Units involving can additionally relabel the horizontal and vertical axes.
In an Einstein-de Sitter universe, and . The Friedmann equation gives
This is an Euler differential equation. Substitution of gives , namely . The Einstein-de Sitter density-growth modes are therefore
The linear growth factor for the growing mode is proportional to , rather than an exponential in cosmic time. These expressions apply while the density contrast is small enough for linear cosmological density perturbations.
Use and , where the growing linear growth factor obeys . The Poisson relation in part (a) immediately becomes
Along a fluid trajectory,
Substitute this into the pressureless equation of motion to obtain
The growth equation cancels the second coefficient. By the definition of the logarithmic growth rate, , while . Therefore
This is the printed when interpreted along the particle. Dividing the Eulerian continuity equation by instead gives
The density derivative in this conservative equation must be Eulerian. If both derivatives were read as material derivatives, the density equation would double-count advection; its material form is . These distinctions make the growth-time form of cosmological dust equations valid beyond the first-order terms, up to the single-stream pressureless-fluid limit. The reparameterization assumes .
A useful decomposition of the linear matter power spectrum, absorbing the Fourier normalization into its amplitude, is
Here is the cold-dark-matter transfer function and the linear growth factor. On scales entering the horizon after matter-radiation equality, , so with near one. Modes entering during radiation domination undergo suppressed growth relative to later entrants, producing a turnover near the equality scale and, well above it, . Thus the linear spectrum asymptotes approximately to . Baryon acoustic oscillations superpose a weak oscillatory pattern. At sufficiently short scales nonlinear gravitational instability invalidates this linear extrapolation.
Figure 1.
Schematic matter power spectrum with equality turnover, baryon oscillations and approximate observational scale ranges
. Details.
The plotted curve illustrates the shape rather than a precision cosmological fit; the bands are overlapping approximate scale sensitivities, not sharp experiment boundaries. Cosmic microwave background anisotropies constrain primordial amplitudes and transfer physics on very large through acoustic scales; converting them to a late-time matter spectrum requires a cosmological model. Galaxy redshift surveys constrain approximately –, with galaxy bias and redshift-space distortions modeled. Peculiar velocities constrain large-scale mass fluctuations; cluster abundances constrain the fluctuation normalization near cluster scales. Weak gravitational lensing measures projected matter on intermediate and smaller scales. The Lyman-alpha forest extends constraints to roughly –a few at high redshift, requiring gas, temperature and ionizing-background modeling. These are complementary measurements of the underlying spectrum rather than interchangeable direct Fourier measurements.
The matter fraction controls the equality turnover. From and the Friedmann equation at equality,
The measured background radiation fixes , so the physical turnover scale measures approximately . In units its location mainly constrains ; independent knowledge of is needed to extract . Standard radiation content gives the useful estimate , approximately for , .
The cosmological baryon density sets the inertia of the pre-decoupling photon–baryon fluid and changes its sound horizon and the oscillation pattern. The oscillation contrast and broadband suppression depend on the baryon fraction ; the acoustic positions depend on the sound horizon and the distance conversion. The cosmic microwave background peak-height pattern independently constrains . Combining these with the equality scale and a distance or Hubble constant determination permits estimates of both requested density fractions. Spectral tilt, neutrino free streaming, galaxy bias and nonlinear evolution must be fitted rather than assigning every shape change to baryons. The acoustic signature was already observed in the 2005 galaxy correlation measurement; the first-year microwave-background parameter analysis demonstrates the complementary parameter constraints.
The broad turnover near marks the matter-radiation equality scale. Modes entering the Hubble radius after matter-radiation equality grow efficiently during matter domination. Modes entering earlier have their cold-matter growth delayed by radiation domination: their subhorizon growth is logarithmic rather than proportional to . This is the Mészáros effect. The smaller undulations on the declining side are baryon acoustic oscillations inherited from the photon-baryon fluid, with an oscillatory wavenumber spacing set approximately by divided by the sound horizon. They are damped at sufficiently short scales by photon diffusion and the transfer of the baryon perturbation into the total matter field. The plotted spectrum is a linear prediction; nonlinear evolution changes the actual sufficiently small-scale spectrum at redshift zero.
For an adiabatic growing mode, the comoving curvature perturbation transfers to the matter overdensity as
up to an irrelevant overall sign convention for . Normalize the cosmological transfer function by at small and the linear growth factor by in matter domination. The relation distinguishes a dimensional matter power spectrum from a dimensionless cosmological power spectrum: the dimensional curvature spectrum is proportional to . Therefore
For locally constant scalar spectral index , the asymptotic cold-dark-matter power spectrum is
because the cold-dark-matter transfer function tends to one at small and scales as at large . The asymptotes refer to the smooth envelope, not the individual baryonic wiggles. The large- local logarithmic slope is . For the quadratic potential at , derived below, , so the limiting powers are approximately and . Multiplying by to form dimensionless matter power would change both exponents by three; that is not the plotted quantity.
Let be the linear spherical-collapse threshold and let be the linear growth factor, normalized to one at the epoch to which the random density field has been extrapolated. The time-dependent barrier is . The smoothed matter density variance is the variance of the linear density contrast smoothed with a spherical top-hat window function whose comoving mass is . Both quantities must refer to the same linear normalization; equivalently use .
Let denote the comoving number density and write for the positive differential Press-Schechter halo mass function. With the mean comoving mass density, the mass fraction obeys
Differentiating the complementary error function and using gives
This is the requested mass function; is the abundance per unit natural logarithmic mass. For physical rather than comoving number density, multiply by . The normalization in the Press-Schechter formalism includes its factor of two relative to the positive Gaussian tail, so tends to one as the smoothing mass tends to zero when diverges. It is a collapse prescription, not the probability that an arbitrary already nonlinear density field stays Gaussian.
Let and define the cosmological density power spectrum by . The contribution per logarithmic wave-number interval is the dimensionless cosmological power spectrum . The linear cosmological mass variance on mass scale is
with a comoving smoothing radius and the top-hat filter in Fourier space. On halo scales the observed spectrum is consistent with greater variance at smaller masses. Growth therefore brings smaller objects to the collapse threshold earlier, statistically; larger structures assemble by accretion and dark-matter halo mergers. This statistical growth from smaller bound systems to larger ones is hierarchical galaxy formation. It is an ordering of dark-matter assembly, not a rule that every small visible galaxy must precede every large visible galaxy.
For linear modes the linear growth factor gives . During matter domination , so and the characteristic nonlinear mass increases. Once modes become nonlinear, coupling between scales and halo formation change the shape, so the same multiplication cannot be used for the entire late-time spectrum. At low redshift accelerated expansion suppresses linear growth. Galaxy clustering traces the matter spectrum with galaxy bias; it should not be equated directly with an unbiased matter measurement.
The broad linear shape was set by the cosmological transfer function before and around matter-radiation equality, with baryonic acoustic structure also imprinted before recombination. A nearly scale-invariant primordial curvature spectrum has , with close to one. After converting curvature perturbations to matter-density perturbations,
Thus nearly scale-invariant primordial curvature does not mean a constant density . Modes with enter the horizon after equality and have . Modes with enter during radiation domination, when radiation controls the expansion and cold-matter perturbations grow only slowly. The cold-dark-matter transfer function behaves approximately as at large . Hence the density spectrum turns over near :
The turnover records the equality horizon, while the late nonlinear excess records gravitational clustering. On galactic scales the effective slope is greater than , giving the growing small-scale variance needed for hierarchical galaxy formation.
Cold dark matter has negligible primordial thermal velocities and a very short collisionless free streaming length on galactic scales. Warm dark matter has appreciable residual velocities while structure is being seeded; particles stream across small fluctuations and reduce their contrast. Its cosmological transfer function is consequently cut off below a characteristic length, suppressing low-mass halos and delaying their formation. Above that cutoff its assembly can still be hierarchical. The important distinction is the free-streaming scale, rather than the present temperature or an arbitrary particle-mass label.
Linear evolution assumes . When becomes of order unity, overdense regions depart strongly from the Hubble flow and can turn around and collapse. Collisionless dark matter develops multistream motion after trajectories cross; gravitational mixing redistributes energy and produces a bound dark matter halo. The formal infinite-density collapse of an ideal spherical pressureless solution is not the physical endpoint. A roughly virialized halo has and a characteristic virial velocity . Its gas virial temperature is conventionally
It measures the thermal energy associated with the gravitational potential; the numerical factor depends on the velocity-dispersion convention. It does not imply that the collisionless dark matter has a thermodynamic gas temperature.
The unheaded request about baryon conversion efficiency of a halo is also answered here. Define using the cosmic baryon fraction . In small halos, shallow potentials let stellar feedback drive outflows or repeatedly heat star-forming gas; supernova energy per stellar mass is roughly fixed while binding energy per gas mass scales as . Photoheating during reionization also prevents very small halos from retaining or accreting cool gas. Molecular/atomic cooling thresholds further reduce star formation in the smallest systems. These effects make fall toward low mass.
Near , gas can cool efficiently and the potential is deep enough to retain more of it, while a long-lived hot atmosphere and maintenance heating are less effective than in larger systems. At high mass, higher virial temperature and lower cooling efficiency let a substantial hot atmosphere persist. Active-galactic-nucleus feedback can prevent that atmosphere from supplying cold gas and can expel some gas; the cooling-time bottleneck alone is not an adequate explanation for the low stellar fractions of massive groups and clusters. The peak reflects a competition between gas supply/cooling and feedback, rather than complete conversion of all baryons at a sharply universal mass. Its exact location and height depend on epoch, metallicity, gas history and the stellar population included.
The prescribed photon budget requires a collapsed mass fraction , assuming baryons trace the collapsed mass fraction and using the supplied effective photon yield. No extra photon-escape or recombination correction should be counted on top of this stipulated budget.
For spectral slope , . The mass-scale measurement gives
Both relevant epochs are matter-dominated, so the linear growth factor ratio is , independent of any present-day normalization convention. Thus . The Press-Schechter formalism gives
Using the supplied inverse value,
This is the photon-budget reionization threshold in the specified toy model, not a measurement of the full astrophysical reionization history.
To estimate the characteristic mass, use the non-dissipative spherical-collapse model result . At this high redshift the background matter density is practically the critical density. Circular virial velocity then obeys
Consequently
For , this gives and , with physical virial radius about .
At the threshold mass the halo peak height is , so
The present mean matter density is . Hence the differential abundance per logarithmic mass interval is about and
Rounded inputs and the supplied approximate factor justify quoting about . The corresponding physical separation is about at this epoch. The requested halo spacing from a differential mass function uses a logarithmic mass bin of order unity; it is not obtained by identifying with that differential abundance. A cumulative number density above the threshold requires a separate mass integral.
In the Lambda-CDM model, nearly Gaussian primordial density contrasts grow under gravity within an expanding universe containing cold dark matter, ordinary baryons and a cosmological constant. The cold dark matter is effectively collisionless and has negligible pressure on galactic scales. Before recombination, baryons are coupled to the photon fluid: radiation pressure and acoustic oscillations prevent their perturbations from behaving like pressureless matter. After recombination they can fall into the gravitational potentials already established by dark matter, subject to gas pressure and the Jeans mass.
For small density contrasts, evolution is linear. On pressure-free scales the growing mode is multiplied by the linear growth factor , with during matter domination. The cosmological density power spectrum can be written
Here is the cosmological transfer function. For nearly scale-invariant initial conditions, the large-scale matter spectrum behaves approximately as , whereas well inside the matter-radiation equality scale it falls approximately as , until the microscopic dark-matter cutoff matters. This fall of the dimensional does not imply less fluctuation power on every smaller mass scale: the power per logarithmic wavenumber is , and the smoothed matter density variance is obtained by integrating it against a mass-dependent window. Over the relevant cold-dark-matter hierarchy, smaller mass windows generally have larger variance.
The hierarchical galaxy formation picture follows: fluctuations on small mass scales typically reach the nonlinear collapse threshold first, while larger objects assemble later through accretion and dark-matter halo mergers. It is a statistical ordering, not a claim that every small object precedes every large rare peak. When becomes order unity, the linear growth factor is no longer a solution for the local density. Collisionless dark matter develops multistream regions and bound dark matter halos; phase mixing and violent relaxation redistribute orbital energies, and virialized structures approximately obey the virial theorem.
Baryons have an additional nonlinear route. Infall and shocks convert bulk kinetic energy into thermal energy, with characteristic virial temperature
Unlike collisionless dark matter, the gas can lose this energy through radiative cooling. The optically thin gas cooling time is the thermal-energy density divided by the radiative loss rate, for example
The density convention in must agree with the denominator; the astrophysical cooling function can also be defined using instead. The cooling criterion for galaxy formation compares this time with the collapse or supply time. Rapidly cooling gas loses pressure support, contracts and can form stars. Slowly cooling gas remains in a hot atmosphere. Stable virial shocks are not obligatory in every low-mass system: gas can also arrive in cold streams and cool while being accreted.
Angular momentum prevents indefinite radial contraction. Tidal torque theory supplies an initial halo spin, and later mergers change it. If gas radiates energy while retaining much of its specific angular momentum, it settles into a rotationally supported galactic disk rather than reaching the centre. The relation explains why modest halo spin can set a disk radius much smaller than its virial radius of a dark-matter halo. Torques, bars and gravitational encounters can transport angular momentum outwards and feed central concentrations; radiative cooling alone does not remove it.
Galaxy mergers alter stellar structure as well as assembling mass. A major galaxy merger can strongly disturb or destroy an existing galactic disk, randomizing stellar orbits and creating a spheroid through violent relaxation. A gas-rich galaxy merger also permits dissipation, inflow and a burst of star formation; gas left over or accreted afterwards can rebuild a galactic disk. Minor galaxy mergers add stars to outer components, thicken disks and grow bulges. Dry galaxy mergers add stellar mass and can increase size without much new star formation. Halo merging therefore does not imply instantaneous merging of its galaxies: satellite orbital decay requires Chandrasekhar dynamical friction and can take a substantial time.
The atomic and molecular cooling thresholds for galaxy formation supply a lower characteristic scale. Primordial atomic gas cools inefficiently below roughly because electronic excitation is suppressed. The corresponding halo mass is of order at a redshift of order ten, with approximate dependence at fixed threshold temperature. Molecular hydrogen can cool gas at hundreds of kelvin and permit smaller early objects, provided it forms and survives dissociating radiation. Metal-line cooling changes these thresholds after enrichment. Thus the atomic threshold is not an absolute minimum mass for all stellar systems.
At the other end, sufficiently massive dark matter halos have high virial temperatures and low-density hot gas. Above the strong atomic-line-cooling interval, thermal bremsstrahlung has . At comparable halo gas density, , while depends mainly on formation density. Cooling therefore becomes less able to condense all the gas within the available time. This upper galaxy mass from gas cooling argument selects galaxy-sized condensations, broadly halo masses around – in simple low-redshift estimates, rather than single luminous galaxies containing every baryon in a group or cluster. Its numerical scale depends on epoch, metallicity and gas profile. Subsequent galaxy mergers can build larger stellar systems; cooling is not an absolute upper bound on their final mass.
Finally, a Press-Schechter halo mass function has many low-mass objects and a steep high-mass cutoff, while the luminosity function of galaxies also has a faint component and a bright cutoff, often described by a Schechter function. The two shapes are related through the halo-to-galaxy luminosity mapping, not by identifying luminosity with total halo mass. In the idealized one-central-galaxy, no-scatter limit,
If and , then . A constant conversion gives similar shapes. Actual stellar feedback, inefficient low-temperature cooling and reionization suppress faint galaxies relative to small haloes, while long cooling times and active-galactic-nucleus feedback suppress luminosity at large halo masses. Satellites, scatter, stellar populations and dust further affect the correspondence. The halo hierarchy supplies the gravitational framework; cooling, angular momentum and feedback determine which parts become luminous galaxies.
At fixed mass, decreasing redshift increases the linear growth factor and reduces the halo peak height . During matter domination also falls as the universe expands. As long as , both effects reduce the positive abundance-growth rate, so grows towards lower redshift. In the very rare-tail limit provides the approximate trend.
There is a limit to calling this an increase time. At , the fixed-mass Press-Schechter halo mass function reaches its maximum as a function of time; its logarithmic growth rate vanishes and diverges. For , the derivative becomes negative, corresponding to net transfer of small objects into more massive systems, so the signed inverse is no longer a positive doubling or increase time. At still later epochs the cosmological constant further suppresses the linear growth factor. These qualifications prevent extrapolating the positive high-redshift growth trend indefinitely.
For fixed epoch and the increasing branch , increasing mass lowers the smoothed matter density variance, hence raises the halo peak height . The rate therefore increases and its inverse decreases. Rarer, more massive haloes have a shorter fractional abundance-growth time, even though their actual number density is much smaller. They lie farther into the exponential tail of the Press-Schechter halo mass function, so a small change in the linear growth factor causes a large fractional change. This statement compares fixed-mass bins in the same cosmology and concerns relative growth, not the time for one halo to assemble all its mass.
Normalize the linear growth factor to . The collapse overdensity at the collapse epoch is for the matter-dominated spherical-collapse model. The present-extrapolated spherical-collapse barrier is
It is the initial linear overdensity, extrapolated to today, required to collapse by ; it is not the nonlinear density contrast of a virialized halo. The smoothed matter density variance is the variance of the linear density contrast smoothed on a Lagrangian comoving scale containing mass . For a spherical top-hat filter,
Consequently for scale-independent linear growth. The equivalent halo peak height conventions are ; using both an evolved barrier and an evolved variance would count growth twice.
To calculate a number-density growth time, use a narrow fixed-mass bin, not the collapsed mass fraction itself. Differentiate the Press-Schechter formalism mass fraction and divide the mass density in the resulting interval by . This gives the Press-Schechter halo mass function
At fixed , the mass factor and logarithmic slope do not depend on time. Since , the Press-Schechter abundance growth at fixed mass is
The prefactor matters here: the exponential-only rare-peak approximation drops the minus one and is accurate only for .
Read the supplied variance relation as a mass-scale calibration using the numerical velocity label given for the selected population. It gives . Matter domination between the two high-redshift epochs gives , so and . At the epoch in question,
Thus the requested fixed-mass-bin growth estimate with that calibration is
An exponential-only approximation gives about , which is somewhat shorter because this is only a roughly two-sigma population.
The wording leaves two sample conventions worth distinguishing. First, an actual virial velocity of a spherical-overdensity halo is epoch-dependent at fixed mass. If the variance fit is instead calibrated using physical virial velocities at redshift three, the same mass whose velocity is at redshift nineteen has . The halo virial-velocity conversion between epochs then gives , and for a fixed-mass bin. The two numerical answers reflect the velocity-label convention in the supplied fit, not two ways of differentiating one fixed fit.
Second, a cumulative number density is , not simply : the latter is a mass fraction divided by a threshold mass, not the number of objects. The cumulative derivative is an abundance-weighted average of over that integral. If the supplied power-law variance fit is extended over all larger masses, with the same velocity-label convention, direct integration gives a cumulative-number growth time of about instead. A sample maintained at fixed physical velocity at successive epochs additionally moves its mass boundary and needs a selection convention. The numerical estimate above explicitly uses the ordinary fixed-mass differential interpretation. These distinctions are important when an exact growth time rather than a rare-tail estimate is intended.
Let be the dimensional cosmological density power spectrum and its dimensionless power spectrum. In linear theory,
where is the linear growth factor and is the cosmological transfer function. An approximately scale-invariant primordial curvature spectrum corresponds, with the conventional matter-spectrum normalization, to with . This does not say that the dimensional density spectrum is white noise or that the early curvature and density perturbations are the same variable.
The crucial scale dependence develops around horizon entry and matter-radiation equality. Modes entering during radiation domination undergo only logarithmic cold-matter growth, whereas modes entering later avoid that suppression. For cold dark matter, the large-scale transfer tends to a constant and the small-scale transfer behaves schematically as . Thus
For , the dimensional spectrum turns from an approximately rise to an approximately decline. Multiplying by shows that the dimensionless small-scale power still increases slowly. The resulting variance of density fluctuations smoothed on mass generally increases as decreases, over the scales relevant to galaxy assembly.
After equality, and over scales where linear growth is scale-independent, the transfer shape is approximately preserved while the amplitude grows. During matter domination, ; late accelerated expansion slows that growth. A region collapses when its linearly evolved density contrast reaches the linear spherical-collapse threshold. Smaller mass scales, having larger variance, typically reach this threshold earlier. Small haloes form first, then accrete matter and merge into larger haloes; galaxies form from cooling baryons in those wells and themselves merge. This small-to-large assembly is hierarchical galaxy formation. It is a statistical trend: rare high peaks can produce unusually massive early objects, and gas cooling and feedback prevent a one-to-one identification of halo collapse with star formation.
For cold dark matter, random particle velocities are small enough that collisionless free streaming erases little power on galactic scales. Warm dark matter retains a larger early velocity dispersion and a larger free-streaming length, suppressing power below a finite scale. There are consequently fewer low-mass haloes, delayed formation of the smallest galaxies, and a lower limit to the hierarchy that can develop. Above the cutoff, warm-matter structure can still assemble hierarchically. The distinction concerns the particles' velocity history and transfer function, not the thermal temperature of gas in a galaxy today. The broad late-time shape was largely established by early horizon-entry physics, with subsequent nonlinear collapse and merging altering the spectrum at high wavenumber.
Let be the linear growth factor normalized to . The present-extrapolated spherical-collapse barrier is , where is the linearly evolved overdensity required for spherical collapse at that epoch. The smoothed matter density variance is the root-mean-square linear density contrast after smoothing on a Lagrangian region containing mass ; . Thus the barrier-to-variance ratio is equivalently .
For a spherical top-hat window function with comoving radius , . A scale-free spectrum gives
Changing variable to gives the scale-free smoothed density variance scaling with . For a top-hat the scale-free integral converges for ; the inferred index below lies in this interval.
Differentiate the Press-Schechter formalism mass fraction and divide the resulting mass density by halo mass. Writing ,
This is a comoving halo abundance. Introduce by ; then and
For the constant baryon conversion mapping of halo and stellar mass, identify one counted galaxy with each halo and neglect scatter and subhalo multiplicity. If is the fraction of the halo's baryons incorporated into stars, its stellar-to-halo mass ratio is , with , and . The corresponding stellar characteristic mass is . Transforming with gives the power-law Press-Schechter stellar mass function:
It has the required power-law low-mass behavior and stretched exponential cutoff. Matching both exponents and the normalization determines the effective spectral index and stellar conversion efficiency.
A larger cosmological constant changes both the expansion history and the growth history: it alters cosmological distances, the baryon acoustic oscillation scale as observed in angle and redshift, and the time at which growth begins to freeze. The new interaction in part (d) leaves the background Friedmann equation unchanged but changes clustering abruptly from onward.
Redshift-resolved measurements of the linear growth factor, weak gravitational lensing, galaxy clustering, and cosmological distance measures can therefore compare geometry with growth. Equal present amplitudes do not imply equal , and only the higher-dark-energy model changes the distance-redshift relation. This breaks the degeneracy.