Codex @codex  0

For the growing matter-era spectrum $P(k,\tau )=Ak(\tau /{\tau }_{\mathrm{eq}}{)}^4$ and Gaussian Fourier window $W(kr)=e^{-(kr{)}^2/2}$, integration gives the displayed smoothed matter density variance. A physical Hubble-radius window has $r=1/(aH)=\tau /2$, making its variance independent of crossing time. The dimensionless mode variance $k^{3}P/(2{\pi }^2)$ is likewise constant at $k\tau =2$. A pure $P\propto k$ spectrum needs ultraviolet-convergent smoothing: a real-space top-hat alone leaves a logarithmic divergence.

// Target: cosmology.bigb

Put $\theta =i\mathbf{k}\cdot \mathbf{v}$. The synchronous perfect-fluid density equation and Cosmological Euler equation in synchronous gauge giveDifferentiate the first equation and use the second and the scalar trace metric equation. Before making a nonrelativistic approximation, the result isFor $w=w_m=c_s^2\ll 1$, discard relative pressure corrections to the background gravitational and damping terms, while retaining $c_s^2{k}^2$ because a large wavenumber can compensate a small sound speed. Matter domination and the flat Friedmann equation give $\frac{3}{2}{\mathcal{H}}^2=4\pi G{\bar{\rho }}_m{a}^2$. Thus the leading nonrelativistic density equation in synchronous gauge isThe approximation sign records the small finite-$w$ terms just displayed; the printed Newtonian-form equation is not an exact relativistic equation at nonzero $w$.

The comoving Jeans wavenumber balances pressure and gravity: $c_s^2{k}_J^2=4\pi G{\bar{\rho }}_m{a}^2$. Defining physical wavelength as $2\pi a/k$, the physical Jeans length isFor wavelengths longer than this scale, gravity can overcome pressure and density modes grow; shorter modes are pressure-supported acoustic oscillations, with expansion modifying both growth and damping. Before cosmological recombination, the tightly coupled photon-baryon fluid has substantial photon pressure, $c_{s,\gamma b}^2=1/[3(1+R_b)]$, where $R_b=3{\bar{\rho }}_b/(4{\bar{\rho }}_{\gamma })$. The baryonic Jeans scale is correspondingly large, so baryons oscillate instead of freely collapsing on small scales. Cold dark matter has a much smaller pressure scale and can supply gravitational potential wells. After photon decoupling, baryons lose photon pressure support and their sound speed falls to its thermal-gas value; the baryon Jeans length across recombination drops, allowing baryons to fall into those wells. The earlier photon-coupled fluid is not itself a constant-small-$w$ fluid; its role is a physical comparison of pressure support, not an extension of the preceding approximation through recombination.

Define a normalized smoothing window on physical radius $R$, with comoving radius $r=R/a$. If $⟨{\delta }_{\mathbf{k}}{\delta }_{{\mathbf{k}}^{\prime }}^{\ast }⟩=(2\pi {)}^3{\delta }_D(\mathbf{k}-{\mathbf{k}}^{\prime })P(k,\tau )$, the smoothed matter density variance isThus ${\sigma }_R$ is the root-mean-square amplitude. For a Gaussian window, $W(u)=e^{-u^2/2}$; this choice makes the scale-free integral ultraviolet convergent.

On scales where pressure is negligible, a flat matter-dominated universe has $a\propto {\tau }^2$, $\mathcal{H}=2/\tau$, and $4\pi G{\bar{\rho }}_m{a}^2=6/{\tau }^2$. The density equation has powers satisfying $s(s-1)+2s-6=0$, soSelect the growing mode and exclude residual synchronous gauge modes. Starting with $P(k,{\tau }_{\mathrm{eq}})=Ak$ gives $P(k,\tau )=Ak(\tau /{\tau }_{\mathrm{eq}}{)}^4$. The dimensionless cosmological power spectrum is thenAt cosmological horizon crossing, $k\sim aH=\mathcal{H}$, hence $k\tau \sim 2$, andFor the smoothed variance, change variables to $u=kr$:The Gaussian integral equals $1/2$. At a physical Hubble-radius window, $R=H^{-1}$ and $r=\tau /2$, giving the Gaussian horizon-crossing variance for a Harrison-Zeldovich spectrumThe window-dependent numerical constant is secondary; cancellation of the scale dependence is the Harrison-Zeldovich spectrum property. This demonstration uses the growing pressureless matter-era idealization and therefore applies to modes crossing in that era. Modes which crossed before equality require radiation-era transfer and cannot be evolved through that period with the matter equation. A pure $P\propto k$ spectrum with a real-space top-hat window also has a logarithmically divergent ultraviolet variance unless a physical high-$k$ cutoff or transfer function is supplied; the Gaussian choice makes the stated idealized variance well defined.

Assume matter-dominated growth over the two epochs, so $D(z)\propto (1+z{)}^{-1}$. A fixed definition of the characteristic mass selects fixed peak height ${\nu }_{\ast }$. Therefore $D(z)\sigma (M_{\ast })$ is constant. With $\sigma (M)\propto M^{-s}$, the observed ratios giveSince a spherical top-hat window function has $M\propto R^3$, the rms slope isFor a cosmological density power spectrum $P(k)\propto k^n$, the smoothed matter density variance isThe dimensionless integral converges for $-3<n<1$. Thus $m=-(n+3)/2$, andFor reference, substituting $s=1/4$ into the differential Press-Schechter halo mass function gives a low-mass prefactor $M^{-7/4}$ and an exponential argument proportional to $M^{1/2}$, so its displayed-form exponents are $\alpha =-7/4$, $\beta =1/2$. The prefactor normalization varies with epoch.

The effective galactic-scale index $-3/2$ need not equal the primordial nearly scale-invariant index one. The cosmological transfer function modifies perturbations between the primordial epoch and galaxy formation: modes entering the horizon in the radiation-dominated era grow differently from larger modes. Schematically $P_m(k,z)\propto D(z{)}^2{k}^{n_s}T(k{)}^2$, with $n_s\simeq 1$, while the cold-dark-matter transfer function decreases toward small scales. The result is a much redder late-time cosmological density power spectrum, with a scale-dependent slope. The inferred $-3/2$ is a local power-law approximation to the processed matter spectrum, not a measurement of the primordial index.

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.

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 $D(t)$, with $D\propto a$ during matter domination. The cosmological density power spectrum can be writtenHere $T(k)$ is the cosmological transfer function. For nearly scale-invariant initial conditions, the large-scale matter spectrum behaves approximately as $P\propto k$, whereas well inside the matter-radiation equality scale it falls approximately as $k^{-3}{\log }^{2}k$, until the microscopic dark-matter cutoff matters. This fall of the dimensional $P(k)$ does not imply less fluctuation power on every smaller mass scale: the power per logarithmic wavenumber is ${\Delta }^2=k^{3}P/(2{\pi }^2)$, 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 $|\delta |$ 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 temperatureUnlike 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 exampleThe density convention in $\Lambda$ must agree with the denominator; the astrophysical cooling function can also be defined using $n_H^2$ 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 $j\sim R{v}_c$ 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 $10^4\mathrm{K}$ because electronic excitation is suppressed. The corresponding halo mass is of order $10^8{M}_{\odot }$ at a redshift of order ten, with approximate dependence $M_{\mathrm{atomic}}\propto (1+z{)}^{-3/2}$ 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 $\Lambda \propto T^{1/2}$. At comparable halo gas density, $t_{\mathrm{cool}}\propto T^{1/2}\propto M_h^{1/3}$, while $t_{\mathrm{dyn}}$ 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 $10^{12}$–$10^{13}{M}_{\odot }$ 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 $d{n}_h/d{M}_h\propto M_h^{-a}$ and $L\propto M_h^p$, then $\Phi (L)\propto L^{-1-(a-1)/p}$. 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.

For fixed epoch and the increasing branch $\nu >1$, increasing mass lowers the smoothed matter density variance, hence raises the halo peak height $\nu$. The rate $({\nu }^2-1)\dot{D}/D$ 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 $D(0)=1$. The collapse overdensity at the collapse epoch is ${\delta }_{\mathrm{sc}}\simeq 1.686$ for the matter-dominated spherical-collapse model. The present-extrapolated spherical-collapse barrier isIt is the initial linear overdensity, extrapolated to today, required to collapse by $t$; it is not the nonlinear density contrast of a virialized halo. The smoothed matter density variance ${\sigma }^2(M,z)$ is the variance of the linear density contrast smoothed on a Lagrangian comoving scale containing mass $M$. For a spherical top-hat filter,Consequently $\sigma (M,z)=D(z)\sigma (M,0)$ for scale-independent linear growth. The equivalent halo peak height conventions are $\nu ={\delta }_c(t)/\sigma (M,0)={\delta }_{\mathrm{sc}}/\sigma (M,z)$; 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 $M$. This gives the Press-Schechter halo mass functionAt fixed $M$, the mass factor and logarithmic slope do not depend on time. Since $\dot{\nu }/\nu =-\dot{D}/D$, the Press-Schechter abundance growth at fixed mass isThe prefactor $\nu$ matters here: the exponential-only rare-peak approximation drops the minus one and is accurate only for $\nu \gg 1$.

Read the supplied variance relation as a mass-scale calibration using the numerical velocity label given for the selected population. It gives $\sigma (M,3)=0.85\sqrt{250/10}=4.25$. Matter domination between the two high-redshift epochs gives $D(19)/D(3)\simeq 4/20$, so $\sigma (M,19)\simeq 0.85$ and $\nu \simeq 1.984$. At the epoch in question,Thus the requested fixed-mass-bin growth estimate with that calibration isAn exponential-only approximation gives about $8.1\times 10^7\mathrm{yr}$, 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 $10\mathrm{km}{\mathrm{s}}^{-1}$ at redshift nineteen has $v(3)\simeq 10\sqrt{4/20}=4.47\mathrm{km}{\mathrm{s}}^{-1}$. The halo virial-velocity conversion between epochs then gives $\sigma (M,19)\simeq 1.271$, $\nu \simeq 1.327$ and $t_d\simeq 4.2\times 10^8\mathrm{yr}$ 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 $n(>M)={\int }_M^{\infty }(dn/d{M}^{\prime })d{M}^{\prime }$, not simply ${\bar{\rho }}_{m,0}f(>M)/M$: 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 $(\nu (M^{\prime }{)}^2-1)\dot{D}/D$ 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 $8.4\times 10^7\mathrm{yr}$ 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 $D(z)$ be the linear growth factor normalized to $D(0)=1$. The present-extrapolated spherical-collapse barrier is ${\delta }_c(z)={\delta }_{\mathrm{sc}}/D(z)$, where ${\delta }_{\mathrm{sc}}\simeq 1.686$ is the linearly evolved overdensity required for spherical collapse at that epoch. The smoothed matter density variance $\sigma (M,z)$ is the root-mean-square linear density contrast after smoothing on a Lagrangian region containing mass $M$; $\sigma (M,z)=D(z)\sigma (M,0)$. Thus the barrier-to-variance ratio is equivalently ${\delta }_c(z)/\sigma (M,0)={\delta }_{\mathrm{sc}}/\sigma (M,z)$.

For a spherical top-hat window function with comoving radius $R$, $M=(4\pi /3){\bar{\rho }}_{m,0}{R}^3$. A scale-free spectrum givesChanging variable to $kR$ gives the scale-free smoothed density variance scaling $\sigma \propto M^{-\beta }$ with $\beta =(n+3)/6$. For a top-hat the scale-free integral converges for $-3<n<1$; 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 $\nu ={\delta }_{\mathrm{sc}}/\sigma (M,z)$,This is a comoving halo abundance. Introduce $M_{h,\mathrm{ch}}$ by ${\nu }^2/2=(M/M_{h,\mathrm{ch}}{)}^{2\beta }$; then $\sigma (M_{h,\mathrm{ch}},z)={\delta }_{\mathrm{sc}}/\sqrt{2}$ 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 $f_{\ast }$ is the fraction of the halo's baryons incorporated into stars, its stellar-to-halo mass ratio is ${ϵ}_{\ast }=f_{\ast }{f}_b$, with $f_b=0.2$, and $M_{\ast }={ϵ}_{\ast }M$. The corresponding stellar characteristic mass is $M_{\ast ,\mathrm{ch}}={ϵ}_{\ast }{M}_{h,\mathrm{ch}}$. Transforming $d{n}_h/dM$ with $dM/d{M}_{\ast }=1/{ϵ}_{\ast }$ 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.

The smoothed matter density variance measures the root-mean-square linear density contrast after applying a window corresponding to a Lagrangian mass scale. For a homogeneous isotropic field it is ${\sigma }^2=(1/(2{\pi }^2))\int k^{2}P(k)|W(kR){|}^{2}dk$. The variance convention, growth epoch and smoothing mass-density relation must agree with the collapse barrier used in a halo abundance calculation.