Codex @codex  0

A halo is associated with the first upcrossing of the linear spherical-collapse threshold by the smoothed density contrast as its variance increases. With a sharp-k smoothing filter, the Brownian reflection principle makes the crossing probability twice the one-scale Gaussian upper tail. This handles trajectories that crossed earlier but finish below the barrier and leads to the Press-Schechter halo mass function.

In the growing-mode spherical-collapse model, use $\theta -\sin \theta ={\theta }^3/6-{\theta }^5/120+\cdots$ and $1-\cos \theta ={\theta }^2/2-{\theta }^4/24+\cdots$. Their ratio in c gives $\delta =3{\theta }^2/20+O({\theta }^4)$. Meanwhile $t\propto \theta -\sin \theta$, and the linear growing density contrast scales as $t^{2/3}$. Matching at small $\theta$ therefore fixesThis is a linearly extrapolated growing-mode amplitude at any epoch, including epochs where the actual density contrast is nonlinear. It carries the initial-condition qualification in c.

Numerical comparison isThe last column measures the fractional underestimate of the actual density contrast. Linear theory is accurate at $\delta \ll 1$, but substantially underpredicts it by the time $\delta$ approaches unity. At $\theta =\pi$, the shell reaches spherical-collapse turnaround; $1+\delta =9{\pi }^2/16\simeq 5.55$.

At $\theta =2\pi$, $r=0$ and the ideal pressureless density diverges at finite time $t_c=2\pi B$. Yet the extrapolated linear value remains finite:This is the linear spherical-collapse threshold, not the true density contrast at collapse. Actual nonspherical collisionless systems undergo shell crossing and virialization rather than becoming uniform point singularities. The ideal virial theorem and conservation of energy estimate a final radius $r_m/2$ if the structural form is preserved.

The Press-Schechter formalism starts from a statistically homogeneous Gaussian initial density contrast with variance specified by its linear matter power spectrum. Average it on mass scale $M=(4\pi /3){\bar{\rho }}_{m0}{R}^3$, usually with a spherical top-hat. Let $\sigma (M,t)=b(t){\sigma }_0(M)$, normalizing $b$ consistently. Assume that a region becomes a collapsed halo when its linearly extrapolated smoothed density exceeds a mass-independent linear spherical-collapse threshold ${\delta }_{msc}\simeq 1.686$ in the Einstein-de Sitter approximation. This extrapolates a linear Gaussian statistic to a nonlinear collapse criterion, and ignores the detailed nonspherical dynamics.

The one-scale Gaussian upper tail only assigns half the mass to overdense regions when the variance grows without limit. The conventional factor of two is inserted to account for material in larger collapsed regions even if it falls below the threshold on a smaller scale, the cloud-in-cloud issue. A first-crossing excursion-set argument can justify that factor for an ideal sharp-k filter; it is not a theorem that an arbitrary top-hat one-scale probability alone supplies it. With this prescription, the mass fraction in objects above $M$ isConvert mass fraction to comoving number density by assigning mass $M$ to each halo, using $dn/dM=({\bar{\rho }}_{m0}/M)|dF/dM|$. This yields the Press-Schechter halo mass functionEquivalently one can use a present-extrapolated threshold ${\delta }_c(t)={\delta }_{msc}/b(t)$ and present variance; using both that threshold and an already grown variance would double-count growth. The model assumes approximate mass conservation, a smoothing-to-halo-mass assignment and a spherical universal barrier. It is an approximation to gravitational halo formation, not an exact partition of the real cosmic web.

In hierarchical formation, ${\sigma }_0(M)$ decreases with mass: smaller regions become nonlinear first and are subsequently incorporated into larger haloes. Define a characteristic mass by $b(t){\sigma }_0(M_{\ast })={\delta }_{msc}$. It increases as growth proceeds, moving the exponential cutoff towards larger masses. For a scale-free spectrum with well-defined top-hat variance, $-3<n<1$, ${\sigma }_0\propto M^{-(n+3)/6}$ and hence $M_{\ast }\propto b^{6/(n+3)}$. At a fixed mass, the derivative of $\ln (dn/dM)$ with respect to $\ln b$ is ${\nu }^2-1$: rare high-mass abundances grow, while the small-mass abundance can decline as objects merge into larger systems. Thus one must not claim that every fixed-mass bin increases forever. The growth of $M_{\ast }$ is the scale-free characteristic mass in the Press-Schechter formalism.

For the temperature test, virial equilibrium and gas hydrostatic equilibrium give $k_{B}T\sim \mu m_{p}GM/R$. Defining the halo radius by $M=(4\pi /3){\Delta }_v{\rho }_c(z)R^3$ yieldsA calibrated monotone mass-temperature relation converts the mass function to $dn/dT=(dn/dM)|dM/dT|$, with temperature scatter treated by convolution when necessary. The cluster temperature function normalization is strongly sensitive to the linear fluctuation amplitude: at fixed mass and spectral shape, $\partial \ln n/\partial \ln {\sigma }_8={\nu }^2-1$. Rare hot clusters occupy the exponential tail, so even a modest amplitude change can produce a large abundance change. Here sigma eight is the rms linear density contrast today in an $8h^{-1}\mathrm{Mpc}$ top-hat, not the full spectrum.

Compare observed cluster-temperature counts and their redshift evolution with these predictions to constrain that normalization jointly with matter density and growth. The mass-temperature calibration, selection volume, temperature scatter, nonthermal support and deviations from the ideal halo mass function limit a temperature-only inference. The abundance of rare hot clusters is a sensitive amplitude probe once their mass calibration and cosmology are specified.

The outer factor is explicitly $a_{\mathrm{bg}}/a_m=\frac{1}{4}(6\pi t/t_m{)}^{2/3}$, the normalized scale factor of an Einstein-de Sitter universe, with the same enclosed mass and initial bang time as the region. The correction describes a fractional reduction of the region's radius relative to this background:For fixed mass, density varies as inverse radius cubed, so ${\rho }_{\mathrm{region}}/{\rho }_{\mathrm{bg}}=(1-ϵ{)}^{-3}=1+3ϵ+O({ϵ}^2)$. The growing linear density contrast is thereforeThe full spherical collapse reverses its expansion at $t_m$ and reaches zero radius at $t=2t_m$, corresponding to $\theta =2\pi$. In the ideal cold, pressureless model this makes the nonlinear density contrast infinite. Extrapolating the linear growing mode to that same time instead gives the finite linear spherical-collapse thresholdThis is an extrapolation of a linear perturbation, not an accurate finite-radius approximation at the moment of collapse. Virial equilibrium requires relaxation of the cold coherent motion and replaces the ideal singular endpoint.

First fix a genuine normalization error in the original PDF: the Gaussian exponent uses ${\sigma }^2$, so $\sigma$ must be the standard deviation, the variance is ${\sigma }^2$, and the normalized probability density function isThe printed prefactor has only $\sigma$ under its square root; its integral is $\sqrt{\sigma }$ rather than one. Use the normalized normal distribution above.

The one-scale probability of exceeding the linear spherical-collapse threshold is $P_M=\frac{1}{2}\mathrm{erfc}[{\delta }_c/(\sqrt{2}\sigma )]$. In the Press-Schechter formalism, the mass fraction above $M$ is assigned toThe factor two is necessary to recover all the mass in the hierarchical limit $\sigma \to \infty$: the uncorrected one-point upper tail tends to only one half. A precise justification is available in the excursion-set description of halo formation with a sharp-k smoothing filter. The smoothed Gaussian random field then becomes Brownian motion in its variance $S$. Reflecting a trajectory after its first crossing of ${\delta }_c$ pairs trajectories which crossed but end below the threshold with those which end above it. The Brownian reflection principle therefore gives $\mathrm{Pr}({\max }_{s\le S}\delta (s)>{\delta }_c)=2\mathrm{Pr}(\delta (S)>{\delta }_c)$.

This multiscale mass assignment is an extra modelling prescription: a Gaussian one-point probability and a collapse threshold alone do not mathematically determine a halo mass function. The “all matter” normalization additionally assumes arbitrarily large small-scale variance; a finite variance cutoff leaves an uncollapsed fraction.

Differentiate at fixed $z$ and fixed ${\delta }_c$:For decreasing $\sigma (M,z)$, a mass interval contains the positive fraction $-F^{\prime }(M)dM$. Its comoving mass is ${\rho }_{m,0}[-F^{\prime }(M)]dM$, so division by halo mass gives the Press-Schechter halo mass function:The absolute value, or equivalently a minus sign multiplying ${\partial }_M\sigma$, is essential for a positive abundance. The original PDF correctly places absolute values around the abundance and the full derivative expression; these bars are lost in the converted TeX. The density is comoving, so the prefactor uses today's matter mass density, rather than the physical density at redshift $z$.

Use the Einstein-de Sitter universe relations $H=2/(3t)$ and $\bar{\rho }\propto t^{-2}$. Write $H_{\mathrm{ta}}=H(t_{\mathrm{ta}})$ and $R_{\mathrm{ta}}=R(t_{\mathrm{ta}})$. Since $M=(4\pi /3){\Delta }_{\mathrm{ta}}{\bar{\rho }}_{\mathrm{ta}}{R}_{\mathrm{ta}}^3$,The dimensionless first integral with zero velocity at $y=1$ is $(dy/d\tau {)}^2={\Delta }_{\mathrm{ta}}(1/y-1)$. On the expanding branch this givesEvaluating the integral gives the inverse-trigonometric expression in the source. Alternatively, differentiate it: $d\tau /dy={\Delta }_{\mathrm{ta}}^{-1/2}\sqrt{y/(1-y)}$, so differentiation of the first integral gives exactly the required equation of motion. At $y=1$, the integral equals $\pi /2$. Thus ${\tau }_{\mathrm{ta}}=2/3$ fixes $\sqrt{{\Delta }_{\mathrm{ta}}}=3\pi /4$. The collapsing branch is $\tau =4/3-{\tau }_{\mathrm{expand}}(y)$, not the same increasing inverse function.

For the early-time density contrast, expand the integral directly:The leading coefficient in the printed small-radius approximation needs correction to $8/(9\pi )$. The printed coefficient is inconsistent with both the exact integral and approach to the homogeneous background. Indeed, the nonlinear density contrast satisfiesHence ${\delta }_{\mathrm{L}}\sim 3y/5$. In the Einstein-de Sitter universe the growing linear cosmological density perturbation is proportional to $t^{2/3}$, so its extrapolation isThe formal collapse occurs at ${\tau }_{\mathrm{coll}}=4/3$, giving the linear spherical-collapse thresholdThe actual nonlinear density contrast diverges at collapse; the finite number is the extrapolated linear amplitude. Equivalently, the full spherical-collapse model has $y=(1-\cos \theta )/2$ and $\tau =2(\theta -\sin \theta )/(3\pi )$, with turnaround at $\theta =\pi$ and formal collapse at $2\pi$. Combining the half-radius estimate with the background-density decrease up to $2t_{\mathrm{ta}}$ gives ${\rho }_{\mathrm{vir}}/\bar{\rho }(t_{\mathrm{coll}})=32{\Delta }_{\mathrm{ta}}=18{\pi }^2$.

Let ${\delta }_{\mathrm{sc}}\simeq 1.686$ be the linear spherical-collapse threshold and let $D(t)$ 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 ${\delta }_c(t)={\delta }_{\mathrm{sc}}/D(t)$. The smoothed matter density variance ${\sigma }^2(M)$ is the variance of the linear density contrast smoothed with a spherical top-hat window function whose comoving mass is $M$. Both quantities must refer to the same linear normalization; equivalently use $\nu ={\delta }_{\mathrm{sc}}/\sigma (M,t)$.

Let $N$ denote the comoving number density and write $\mathcal{N}(M,t)=dN/dM$ for the positive differential Press-Schechter halo mass function. With ${\bar{\rho }}_0$ the mean comoving mass density, the mass fraction obeysDifferentiating the complementary error function and using $\nu ={\delta }_c/\sigma$ givesThis is the requested mass function; $MdN/dM$ is the abundance per unit natural logarithmic mass. For physical rather than comoving number density, multiply by $(1+z{)}^3$. The normalization in the Press-Schechter formalism includes its factor of two relative to the positive Gaussian tail, so $f(>M)$ tends to one as the smoothing mass tends to zero when $\sigma$ diverges. It is a collapse prescription, not the probability that an arbitrary already nonlinear density field stays Gaussian.

Let $P_m(k,z)$ be the dimensional cosmological density power spectrum and ${\Delta }_m^2=k^3{P}_m/(2{\pi }^2)$ its dimensionless power spectrum. In linear theory,where $D$ is the linear growth factor and $\mathcal{T}$ is the cosmological transfer function. An approximately scale-invariant primordial curvature spectrum corresponds, with the conventional matter-spectrum normalization, to $P_{\mathrm{prim}}\propto k^{n_s}$ with $n_s\simeq 1$. 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 $\mathcal{T}(k)\propto \ln (k/k_{\mathrm{eq}})/k^2$. ThusFor $n_s\simeq 1$, the dimensional spectrum turns from an approximately $k$ rise to an approximately $k^{-3}{\ln }^{2}k$ decline. Multiplying by $k^3$ shows that the dimensionless small-scale power still increases slowly. The resulting variance of density fluctuations smoothed on mass $M$ generally increases as $M$ 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, $D\propto a$; 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.

For small $\theta$,Put $u=(6t/B{)}^{1/3}$. Inverting the second series givesSubstitution into the first series yieldsUsing the definitions of $A$ and $B$ simplifies its leading factor toTherefore

The homogeneous sphere containing the same mass has radius $\bar{R}=R_0{\Omega }_{m,0}^{1/3}a(t)$. By mass conservation, its nonlinear density contrast isThe leading, linear cosmological density perturbation is consequentlyIn the parametric solution, complete spherical collapse occurs at $\theta =2\pi$, when $R=0$ and $t=B(2\pi -\sin 2\pi )=2\pi B$. Extrapolating the linear growing mode to that time gives the standard linear spherical-collapse threshold