Cosmological structure formation

Cosmological structure formation is the growth of primordial density perturbations into sheets, filaments, halos, galaxies, and clusters under gravity.

In hierarchical galaxy formation, low-mass dark-matter structures typically collapse before larger structures because the smoothed linear density variance is larger on smaller mass scales. Accretion and mergers assemble increasingly massive haloes. Cooling, feedback and angular momentum determine how baryons form galaxies within that halo hierarchy; the ordering is statistical rather than an absolute rule for every object.

The galactic angular-momentum catastrophe is the excessive loss or redistribution of baryonic specific angular momentum in early simulated galaxy formation, producing overly compact disks or large central concentrations. Premature cooling, gravitational torques, interactions, and resolution contributed. Improved stellar feedback and gas accretion histories can mitigate it; it is a historical simulation problem rather than a theorem against cold dark matter.

Downsizing denotes the tendency for the characteristic mass of actively star-forming galaxies to decrease with time, together with related evidence that more massive quiescent systems formed their stars earlier or over shorter intervals. It does not assert that their dark-matter haloes assembled earlier: stellar formation and halo or stellar assembly are separate histories.

Monolithic galaxy formation idealizes a galaxy as forming in one early rapid dissipative collapse, followed mainly by passive evolution. Rapid formation of the stars is distinct from early assembly of all those stars into their final object; hierarchical galaxy formation can assemble old stellar populations later.

The integrated halo star-formation efficiency is the fraction of a dark matter halo's reference baryonic mass converted into stars: ${ϵ}_{\ast }=M_{\ast }/(f_b{M}_h)$, with cosmic baryon fraction $f_b$. It is distinct from star-formation efficiency per free-fall time. Cooling, gas supply and retention, stellar feedback and active-galactic-nucleus feedback make it strongly dependent on dark matter halo mass. With stellar mass-to-light ratio ${{\rm Y}}_{\ast }$, a simple central-galaxy mapping is $L={ϵ}_{\ast }{f}_b{M}_h/{{\rm Y}}_{\ast }$. Changing this efficiency changes both the mapping and the logarithmic Jacobian between the dark matter halo mass function and the luminosity function.

Selecting a baryonic mass $M_b$ from a normalized cumulative distribution $M_b(<j)\propto j$ fixes a cutoff specific angular momentum $j_d=(M_b/M_{b,\mathrm{tot}})j_{\mathrm{vir}}$. With circular support and dominant enclosed baryonic gravity, the outer radius is $r_d\simeq j_d^2/(G{M}_b)$ and speed $v_d\simeq G{M}_b/j_d$. The mean specific angular momentum is $j_d/2$, which defines a different one-zone radius. The estimate assumes energy dissipation but negligible redistribution of specific angular momentum.

The baryon conversion efficiency of a halo is the fraction of its baryonic mass incorporated into stars. If the halo retains the cosmic baryon fraction $f_b$, its stellar-to-total-mass ratio is $f_b{f}_{\ast }$. This efficiency depends on halo mass, epoch, gas supply and feedback and must not be confused with a star-formation efficiency per free-fall time.

A monotonic luminosity assignment to a counted halo population converts its differential abundance per mass into a luminosity function by a change of variable. If $n_h\propto M_h^{-a}$ and $L\propto M_h^p$, then $\Phi \propto L^{-1-(a-1)/p}$. This simple mapping assumes one counted galaxy per halo, no scatter and a fixed selection. A mass-dependent baryon conversion efficiency of a halo, satellites and stellar-population differences modify it.

The stellar-to-halo mass ratio compares a galaxy's galactic stellar mass with its host's total halo mass. In a simple baryon-conversion model it is ${ϵ}_{\ast }=f_b{f}_{\ast }$. Observational comparisons must distinguish this total-mass ratio from the fraction of baryons converted, and must specify whether central, satellite or all stellar populations are included.

If every counted halo hosts one counted galaxy and the baryon fraction and conversion efficiency are mass-independent without scatter, galactic stellar mass is a constant multiple of halo mass. The galaxy mass function is obtained by changing variable with $dM/d{M}_{\ast }=1/{ϵ}_{\ast }$. Satellites, efficiency variation and scatter prevent this mapping from being generally exact.

Gas in a collapsing halo is heated toward a virial temperature. If its cooling time is short compared with the dynamical time, thermal support is removed and rapid condensation can form a galaxy. If cooling is slower, a hot atmosphere can persist or condense only gradually. This criterion links baryonic galaxy formation to the dark-matter assembly hierarchy without making every halo into one luminous object.

Primordial atomic line cooling becomes effective near the temperature where electronic excitations are accessible, of order ten thousand kelvin. Molecular hydrogen has lower-energy transitions and can cool colder gas if it forms and survives. The corresponding virial temperature thresholds translate into epoch-dependent halo mass thresholds. Radiation, molecular chemistry, metal-line cooling and gas density qualify a single universal minimum galaxy mass.

In simple hierarchical models, increasing system mass raises the virial temperature and can carry primordial gas beyond its efficient atomic line-cooling interval. At assembly densities, the cooling-time constraint then distinguishes galaxy-sized condensations from larger hot groups or clusters. The resulting characteristic mass is approximate, depends on composition and structure, and does not prohibit subsequent mergers of already formed stellar systems.

The Zeldovich approximation maps a fluid element from its Lagrangian coordinate $\mathbf{q}$ to the comoving position $\mathbf{x}=\mathbf{q}+D(t)\mathbf{s}(\mathbf{q})$, where $D$ is the linear growth factor. It follows the initial gravitational displacement along a straight comoving trajectory until trajectories cross.

For $\mathbf{r}=a[\mathbf{q}+b(t)(p(q_1),0,0)]$, put $J=1+b{p}^{\prime }$. Before shell crossing, mass conservation gives $\rho =\bar{\rho }/J$. The dust background equation $3\ddot{a}/a=-4\pi G\bar{\rho }$ and the growing-mode equation $\ddot{b}+2H\dot{b}=4\pi G\bar{\rho }b$ imply ${\nabla }_r\cdot \ddot{\mathbf{r}}=-4\pi G\bar{\rho }+4\pi G\bar{\rho }b{p}^{\prime }/J=-4\pi G\rho$. Thus planar Zeldovich approximation trajectories satisfy the full nonlinear dust force-divergence equation until the map loses invertibility. All three physical coordinates are included in the divergence; the two unperturbed directions still contribute background acceleration.

The cosmological Lagrangian displacement sends an initially labelled matter element from its Lagrangian coordinate $\mathbf{q}$ to $\mathbf{x}(\mathbf{q},t)=\mathbf{q}+\mathbf{\psi }(\mathbf{q},t)$. Its time derivative determines the peculiar velocity through $\mathbf{v}=a\dot{\mathbf{\psi }}$.

For an irrotational displacement $\mathbf{s}=-{\nabla }_q\Phi$, the deformation tensor is the symmetric Hessian $d_{ij}=-{\partial }_i{\partial }_j\Phi$. Its real eigenvalues determine the principal collapse directions.

Shell crossing occurs when the Lagrangian-to-Eulerian map ceases to be one-to-one. Its Jacobian determinant vanishes, the pressureless single-stream density formally diverges, and several streams subsequently occupy the same position.

For the Zeldovich approximation $\mathbf{x}=\mathbf{q}-b{\nabla }_q\psi$, let ${\lambda }_i$ be the eigenvalues of the potential Hessian, with this sign convention. Before shell crossing, mass conservation gives the displayed density. The first positive eigenvalue reaches $b{\lambda }_i=1$ and collapses one direction, creating a sheetlike cosmological caustic. Extrapolating straight trajectories after crossing lets particles pass through without self-consistent gravitational trapping, so the predicted multistream sheet spreads excessively. Collisionless gravitational evolution keeps matter more tightly bound; the linear displacement extrapolation cannot determine its late thickness reliably.

A cosmological caustic is the image of a singular point of the Lagrangian map. The ideal collisionless density diverges there, although velocity dispersion and other small-scale physics regularize the singularity.

Tidal torque theory explains the early angular momentum of a protohalo through the misalignment of its mass second-moment tensor $I_{ij}$ and the Hessian matrix $T_{ij}={\partial }_i{\partial }_j\Phi$ of the surrounding peculiar gravitational potential. At leading order,

The cosmic web is the network of matter sheets, filaments, nodes, and intervening voids produced by anisotropic gravitational collapse.

For comoving position $\mathbf{x}$, canonical momentum $\mathbf{p}=am{\mathbf{x}}^{\prime }$, and Newtonian potential $\varphi$, the dark-matter phase-space density obeys

This is the gravitational, expanding-universe specialization of the Vlasov equation.

The velocity-dispersion tensor is the second central moment of a phase-space distribution,It contributes the stress-divergence term ${\rho }^{-1}{\partial }_j(\rho {\sigma }^{ij})$ to the Euler equation and vanishes in the ideal single-stream limit.