Codex @codex  0

Baryon acoustic oscillations are the relic spatial scale of sound waves in the pre-recombination photon-baryon fluid. The scale appears both as a preferred separation in galaxy clustering and as oscillatory structure in the cosmological density power spectrum.

Modes with $k\gg k_{\mathrm{eq}}$ enter the horizon during radiation domination. Their matter perturbations grow only logarithmically until matter-radiation equality, rather than as $a$. Matching the radiation-era logarithmic growth onto the matter-era growing mode gives the large-$k$ cosmological transfer function scaling shown, up to constants and the matching scale inside the logarithm. For a primordial scalar spectral index $n_s$, the late-time linear cosmological density power spectrum consequently changes from $P(k)\propto k^{n_s}$ at small $k$ to $P(k)\propto k^{n_s-4}{\log }^2(k/k_{\mathrm{eq}})$ at large $k$. These scalings omit baryonic acoustic structure and nonlinear evolution.

Here the cold-dark-matter transfer function is normalized so that ${\delta }_m(k,a)=A(a)k^{2}T(k)\mathcal{R}(k)$, where $\mathcal{R}$ is the primordial comoving curvature perturbation and $A(a)$ is independent of $k$. During radiation domination, horizon entry occurs at $a_{\mathrm{ent}}\propto k^{-1}$. An adiabatic entry amplitude of order $\mathcal{R}$ then acquires a factor $\log (a_{\mathrm{eq}}/a_{\mathrm{ent}})\sim \log (k/k_{\mathrm{eq}})$ by matter-radiation equality. Subsequent matter domination supplies a common growth factor; division by the conventional $k^2$ prefactor therefore gives $T(k)\propto k^{-2}\log (k/k_{\mathrm{eq}})$.

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 =(\rho -\bar{\rho })/\bar{\rho }$ and define the cosmological density power spectrum by $⟨{\delta }_{\mathbf{k}}{\delta }_{{\mathbf{k}}^{\prime }}^{\ast }⟩=(2\pi {)}^3{\delta }_{\mathrm{D}}(\mathbf{k}-{\mathbf{k}}^{\prime })P(k)$. The contribution per logarithmic wave-number interval is the dimensionless cosmological power spectrum ${\Delta }^2(k)=k^{3}P(k)/(2{\pi }^2)$. The linear cosmological mass variance on mass scale $M$ iswith $R$ a comoving smoothing radius and $W$ 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 $P(k,z)=D^2(z)P(k,0)$. During matter domination $D\propto a$, so $P\propto a^2$ 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 ${\mathcal{P}}_{\mathcal{R}}(k)\propto k^{n_s-1}$, with $n_s$ close to one. After converting curvature perturbations to matter-density perturbations,Thus nearly scale-invariant primordial curvature does not mean a constant density $P(k)$. Modes with $k\ll k_{\mathrm{eq}}$ enter the horizon after equality and have $T\simeq 1$. Modes with $k\gg k_{\mathrm{eq}}$ 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 $T\propto \ln (k/k_{\mathrm{eq}})/(k/k_{\mathrm{eq}}{)}^2$ at large $k$. Hence the density spectrum turns over near $k_{\mathrm{eq}}=a_{\mathrm{eq}}{H}_{\mathrm{eq}}/c$:The turnover records the equality horizon, while the late nonlinear excess records gravitational clustering. On galactic scales the effective slope is greater than $-3$, 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 $|\delta |\ll 1$. When $\delta$ 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 $2K+W\simeq 0$ and a characteristic virial velocity $V_{\mathrm{vir}}^2=G{M}_h/r_{\mathrm{vir}}$. Its gas virial temperature is conventionallyIt 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 $f_{\ast }=M_{\mathrm{stars}}/(f_b{M}_h)$ using the cosmic baryon fraction $f_b$. 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 $V_{\mathrm{vir}}^2$. 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 $f_{\ast }$ fall toward low mass.

Near $M_h\sim 10^{12}{M}_{\odot }$, 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.

Take regular adiabatic initial conditions and a common scale-invariant primordial amplitude, $k^{3/2}{\Phi }_{\mathrm{prim}}(k)=\text{constant}$. Let $k_{\mathrm{eq}}={\mathcal{H}}_{\mathrm{eq}}$. The three histories are distinguished by their time of cosmological horizon crossing:

The $9/10$ factor follows from conserving the large-scale adiabatic curvature: the constant potential is proportional to $3(1+w)/(5+3w)$, whose matter-to-radiation ratio is $(3/5)/(2/3)=9/10$. Define the cold-dark-matter transfer function by ${\Phi }_{\mathrm{late}}=(9/10)T(k){\Phi }_{\mathrm{prim}}$. Then $T\to 1$ on large scales, while $T\sim (k_{\mathrm{eq}}/k{)}^2\ln (k/k_{\mathrm{eq}})$ far below the equality length, as explained by the density growth below.

The upper panel sketches $k^{3/2}\Phi$ for the three modes. The lower panel provides the corresponding density contrast histories used in the next part. A negative common primordial potential was chosen so the growing density is positive; this arbitrary phase has no effect on a cosmological density power spectrum.

The curves integrate the ideal coupled radiation-fluid and pressureless-matter equations, rather than patching a pure-radiation solution onto a matter solution. They neglect baryons, free-streaming anisotropic stress, dark energy and nonlinear evolution, consistently with the stated mixture. Dots mark $k=\mathcal{H}$; equality is the dashed vertical line.

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.

Write ${\mathcal{P}}_c(k,\tau )=k^3|{\Delta }_c(k,\tau ){|}^2/(2{\pi }^2)$ for the dimensionless power spectrum. A radiation-era mode enters the horizon at ${\tau }_h\sim k^{-1}$. The early superhorizon law gives ${\mathcal{P}}_c(k,{\tau }_h)\sim A$: the ${\tau }_h^4$ and $k^4$ factors cancel.

For $k{\tau }_{\mathrm{eq}}\ll 1$, the mode remains outside the horizon throughout the radiation era. Its amplitude grows as ${\tau }^2$ both before and after equality, and continues with the same scale-independent growing factor even if it later enters during matter domination. Hence the low-wavenumber shape is preserved:For $k{\tau }_{\mathrm{eq}}\gg 1$, horizon entry occurs during radiation domination. Its subsequent logarithmic amplitude growth supplies a factor $\ln ({\tau }_{\mathrm{eq}}/{\tau }_h)\sim \ln (k{\tau }_{\mathrm{eq}})$ by equality. During the matter era the amplitude grows by $(\tau /{\tau }_{\mathrm{eq}}{)}^2$. Squaring gives the logarithmic high-wavenumber density-spectrum transfer:The horizon-entry and equality matching fixes order-one coefficients and additive constants inside the logarithm. The displayed branches capture the requested leading spectrum and share its schematic normalization; their asymptotic forms must not be joined literally at $k{\tau }_{\mathrm{eq}}=1$, where the high-$k$ logarithm alone vanishes.

Thus the present dimensionless power spectrum rises as $k^4$ on large scales and only as ${\ln }^{2}k$ on small scales, with a smooth bend near the matter-radiation equality scale, $k_{\mathrm{eq}}\sim {\tau }_{\mathrm{eq}}^{-1}$. The ordinary dimensional cosmological density power spectrum instead satisfiesIt has a turnover near equality; the dimensionless spectrum does not. The following comparison of dimensional and dimensionless matter spectra across equality uses a smooth guide with the correct asymptotes, not an exact solution through the transition. Overall vertical normalization and the later scale-independent growth factor are suppressed.

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.

The Mészáros effect suppresses the growth of modes which enter the horizon during radiation domination. Modes with comoving wavenumber $k\gg k_{\mathrm{eq}}$ enter earlier, but gain only a logarithmic growth factor until matter-radiation equality; modes entering during matter domination grow as $a$ after entry. Consequently small-scale amplitudes are reduced relative to an extrapolation of the unsuppressed large-scale cosmological density power spectrum, creating a turnover near the matter-radiation equality scale $k_{\mathrm{eq}}={\mathcal{H}}_{\mathrm{eq}}$ in units $c=1$.

More quantitatively, for adiabatic primordial fluctuations with scalar spectral index $n_s$, the late-time linear spectrum can be written, up to a time-dependent normalization,The cold-dark-matter transfer function approaches a constant at $k\ll k_{\mathrm{eq}}$, whereas logarithmic small-scale matter transfer gives $T(k)\propto (k_{\mathrm{eq}}/k{)}^2\log (k/k_{\mathrm{eq}})$ at large $k$. ThusFor approximately scale-invariant initial curvature fluctuations, $n_s\simeq 1$, the large-$k$ scaling is $k^{-3}{\log }^{2}k$, rather than $k$. This discussion ignores baryonic acoustic structure and nonlinear clustering. The Mészáros equation solution basis includes both a constant and a logarithmic branch during radiation domination; the physical logarithmic growth is selected by matching the earlier radiation forcing, so not every mathematical solution grows logarithmically.

Interpret $P(k)$ in the variance integral as the unsmoothed cosmological density power spectrum of $\delta$. A linear smoothing window gives $P_R(k)=P(k)|\tilde{W}(kR){|}^2$ for the smoothed field; the PDF's description of $P(k)$ as already belonging to ${\delta }_R$ would otherwise count the window twice.

For a scale-free matter power spectrum $P(k)=A{k}^{n_{\mathrm{eff}}}$, substitute $q=kR$ to obtainWhenever the dimensionless integral is finite and nonzero, this proves the scale-free smoothed density varianceThe printed condition $n_{\mathrm{eff}}<3$ is insufficient by itself. For a normalized window with $\tilde{W}(0)=1$, infrared convergence requires $n_{\mathrm{eff}}>-3$. Ultraviolet convergence depends on the filter: a spherical top-hat window function requires $n_{\mathrm{eff}}<1$, a Gaussian filter suppresses every ultraviolet power, and a sharp-k smoothing filter has finite Fourier support. In particular, the monotone relation $R\to \infty \Rightarrow S={\sigma }_R^2\to 0$ requires $n_{\mathrm{eff}}>-3$ in the scale-free model.

At a single variance $S_{\ast }$, the Gaussian random field has ${\delta }_{S_{\ast }}\sim N(0,S_{\ast })$. Thus its endpoint tail isHere $\mathrm{erfc}$ is the complementary error function.

For the crossing probability, the key assumption is independent, symmetric increments as the smoothing scale changes. A sharp-k smoothing filter, $\tilde{W}(kR)=1_{k<1/R}$, supplies this property: decreasing $R$ adds independent Fourier shells of the underlying Gaussian field. Parameterizing their accumulated variance by $S$ gives a Brownian motion with covariance $\min (S,S^{\prime })$. An arbitrary real-space window gives correlated increments and does not justify this argument.

Reflect a trajectory after its first hit of ${\delta }_c$. The Strong Markov property and symmetric future increments preserve its probability, and reflection maps an endpoint below the barrier to one above it. Every endpoint above the barrier has already crossed; the reflected paths give an equally probable set that crossed but ended below. The Brownian reflection principle therefore yields the excursion-set description of halo formation:Its derivative is the Brownian first-passage time densityThe factor of two resolves trajectories that crossed a larger-scale barrier but finish below it on the chosen smaller scale.

Define the halo peak height $\nu ={\delta }_c/{\sigma }_R$. The cumulative probability $\mathcal{P}(>M)$ is a dimensionless mass fraction. With $V_M=M/\bar{\rho }$, its conversion to halo number density givesConsequently the Press-Schechter halo mass function isIt is positive when ${\sigma }_R$ decreases with mass. A chosen mass-radius convention $M\propto R^3$ gives $d\log {\sigma }_R/d\log M=-(n_{\mathrm{eff}}+3)/6$. For a sharp Fourier filter that mass assignment needs a convention, since its real-space kernel is not a localized top-hat volume.

For warm dark matter, the cutoff introduces a new length and invalidates the pure power-law variance scaling near that length. With the same sharp-k filter and $n_{\mathrm{eff}}>-3$, the variance is exactlyThe Brownian walk runs only up to $S_{\max }$ and then stops. The cutoff does not destroy independence of the Fourier shells that are still present. Thus the first-crossing density and derivative mass-function formula remain valid on the invertible part of $S(M)$, using the actual cutoff variance. For $R<k_{\mathrm{WDM}}^{-1}$, $d{\sigma }_R/dM=0$ and this idealized sharp-filter model has no new crossings or new low-mass haloes. Its total collapsed fraction isThis is halo first crossing with a finite variance cutoff: some trajectories never cross, so one must not force the collapsed fraction to one. The uncut power-law abundance cannot be extrapolated to small masses. For other windows, a variance plateau still occurs but the Markov/reflection derivation is not exact; the same mass function is then a modelling approximation, rather than a result established by this calculation.

Evaluate all derivatives at zero long perturbation, where $\nu ={\delta }_c/{\sigma }_{G,s}$, and write $\bar{n}(\nu )$ for the unperturbed abundance. The peak-background split changes the halo peak height toThe chain rule in the specified bias expansion therefore givesNo explicit abundance law is specified, so these derivatives are the general answer. For example, if $\bar{n}$ has the Press-Schechter halo mass function dependence $\nu e^{-{\nu }^2/2}$, then $b_{10}=({\nu }^2-1)/{\delta }_c$ and $b_{01}=2f_{\mathrm{NL}}({\nu }^2-1)$.

For the linear long mode, the cosmological Poisson equation gives ${\phi }_{G,L}(\mathbf{k})={\delta }_{G,L}(\mathbf{k})/\alpha (k)$. Consequently the deterministic galaxy bias and power spectrum areHere $P_m$ is the linear cosmological density power spectrum; stochastic tracer noise has been omitted. This scale-dependent halo bias from local non-Gaussianity is proportional to $k^{-2}$ in the specified $\alpha (k)=2k^2/(3H^2{\Omega }_m)$ convention.

For $b_{10}\ne 0$, plot the ratio $P_g/(b_{10}^2{P}_m)$ to separate the effect from the shape of $P_m$. Positive $f_{\mathrm{NL}}$ increases the ratio as $k$ decreases. Negative $f_{\mathrm{NL}}$ first suppresses it, gives a zero when $\alpha (k)=2|f_{\mathrm{NL}}|{\delta }_c$, and then makes it rise again after the effective galaxy bias changes sign. The auto-power spectrum never becomes negative. A truncation $P_g\simeq b_{10}^2[1+4f_{\mathrm{NL}}{\delta }_c/\alpha (k)]P_m$ is valid only when $|2f_{\mathrm{NL}}{\delta }_c/\alpha (k)|\ll 1$ and must not be extrapolated through the zero. Within the deterministic expression, both signs give $P_g\propto k^{-4}{P}_m$ at sufficiently small $k$. For a nearly scale-invariant primordial spectrum, $P_m\propto k^{n_s}$ in the large-scale transfer limit, so this formal asymptote is $k^{n_s-4}$.

There is a convention issue in calling the specified abundance response an observed galaxy overdensity. The threshold response is naturally a linear Lagrangian halo bias. If $b_{10}$ instead denotes the physical linear Eulerian halo bias, the Lagrangian-to-Eulerian linear bias relation adds one to the density coefficient: $b_{01}=2f_{\mathrm{NL}}{\delta }_c(b_{10}-1)$, and $P_g=[b_{10}+b_{01}/\alpha {]}^2{P}_m$. The boxed formulas follow the abundance expansion explicitly supplied in the paper. The Newtonian treatment also cannot be extrapolated beyond its physical large-scale domain without relativistic projection effects.

Use signed one-dimensional external momenta $k_1+k_2+k_3=0$, with nonzero $k_i$. Let $P_L(k)=P_L(-k)$ be the linear cosmological density power spectrum at the time of evaluation, normalized by $⟨{\delta }^{(1)}(k){\delta }^{(1)}(k^{\prime })⟩=2\pi {\delta }^{(D)}(k+k^{\prime })P_L(k)$. Using $P_L$ absorbs the factor $a^6$ multiplying a product of three initial power spectra. All kernels below are symmetrized one-dimensional cosmological density kernels.

The connected B222 contribution to the one-loop matter bispectrum has a triangular Feynman diagram: each second-order external field has two Gaussian inputs, joined pairwise to the other two vertices. The connected B411 contribution to the one-loop matter bispectrum has a fourth-order vertex joined to both linear external fields, with its remaining two inputs contracted into a loop at that vertex. External stubs indicate the measured momenta; blue internal lines represent linear power spectra.

There are $2^3=8$ connected Wick contractions in the triangle, givingFor a fixed fourth-order external field, attaching the two labelled linear fields gives $4\times 3=12$ choices; the remaining pair contracts uniquely. ThereforeThese are contributions to the reduced connected one-loop matter bispectrum, with its external momentum-conserving Dirac delta function removed. Disconnected contractions have a zero external momentum and are excluded.

For the ultraviolet behaviour, the three $F_2$ denominators multiply to a negative square. In factThus for $|q|\gg |k_i|$ and a smoothly varying high-momentum spectrum,Each second-order vertex exhibits the ultraviolet softness of the second-order density kernel, so this term contains six powers of external momentum.

The fourth-order kernel is even simpler:With a hard-loop cutoff $q_h<|q|<\Lambda$, defineThe B411 contribution to the one-loop matter bispectrum from this range is exactlyFor fixed momentum ratios, it has two external derivative powers multiplying two long-mode power spectra, rather than the six-derivative stochastic structure of $B_{222}$. It is the leading deterministic UV-sensitive contribution for external modes below the nonlinear wavenumber. This is a comparison of derivative orders and generic cutoff sensitivity; an arbitrary specially chosen spectrum or vanishing shape can change their numerical ranking. A convergent integral can still depend on nonperturbative short scales. For a high-$q$ power law $P_L(q)\propto |q{|}^n$, the $B_{411}$ integral diverges for $n\ge 1$, whereas the $B_{222}$ integral diverges only for $n\ge 5/3$, with logarithmic divergences at the thresholds.

At the field level, there are six ways to contract a hard pair within ${\delta }^{(4)}$. Comparing $6F_4$ with $F_2$ givesThe appropriate second-order density Laplacian counterterm is thereforeIndeed its contribution with two linear fields is $-2C(\Lambda ){\sum }_{\mathrm{cyclic}}{k}_3^2{F}_2{P}_L(k_1)P_L(k_2)$, which cancels the displayed hard-loop term. The counterterm sign follows from ${\partial }_x^2\leftrightarrow -k^2$. Its finite coefficient must be matched to short-scale dynamics, as in the effective field theory of large-scale structure. A term proportional to ${\partial }_x^2{\delta }^{(1)}$ similarly renormalizes a linear response, but a purely linear insertion cannot by itself cancel this $⟨{\delta }^{(4)}{\delta }^{(1)}{\delta }^{(1)}⟩$ shape, since an odd Gaussian three-point function vanishes.

For an approximately scale-invariant primordial spectrum, the late cosmological density power spectrum behaves schematically asModes with $k\ll k_{\mathrm{eq}}$ enter the horizon after matter-radiation equality and retain $T(k)\simeq 1$, so $P(k)\sim k^{n_s}$. Modes with $k\gg k_{\mathrm{eq}}$ enter during radiation domination and grow only logarithmically until equality. Their transfer function behaves as $T(k)\sim \log (k/k_{\mathrm{eq}})/(k/k_{\mathrm{eq}}{)}^2$, giving a strongly falling small-scale spectrum. The result is a turnover near the matter-radiation equality scale $k_{\mathrm{eq}}$.

Yes. The baryons carry the oscillatory baryon transfer function inherited from the pre-recombination photon-baryon fluid. Their later infall toward cold-dark-matter overdensities suppresses the relative mode and drives ${\delta }_b/{\delta }_c\to 1$, but gravitational evolution does not remove every scale-dependent acoustic phase. The total cosmological density power spectrum therefore contains the small wiggles known as a baryon acoustic oscillation in the matter power spectrum.

Relative to massless-neutrino growth after $z_{\nu }$,The total contrast has the additional factor $1-f_{\nu }$. Squaring its amplitude to obtain the cosmological density power spectrum givesThe two terms respectively encode the unclustered neutrino fraction and the accumulated slowing of cold-matter growth.