GLOBAL KNOWLEDGE ARK

Cosmology · From density fluctuations to growing structure

An acoustic imprint in the statistics of structure

Use an explicit finite teaching spectrum to distinguish a preferred pair separation from a physical shell, and amplitude smoothing from a complete evolution model.

Preparation: ratios, powers of ten and ordinary arithmetic. Physical quantities and notation are introduced as they are needed.

Back to the chapter

After reading, you can

  • Explain how an early acoustic process can leave a statistical length in later structure.
  • Relate a specified isotropic power spectrum to a two-point correlation with explicit units.
  • Calculate the different effects of a density-amplitude filter on amplitude and power.
01

A preferred separation is a statistical imprint

Before photons and baryons cease to exchange momentum efficiently, pressure and gravity support acoustic motion in their coupled medium. The resulting perturbations can leave a characteristic comoving length in later matter statistics.

The relevant sound-travel integral has the form ∫c_s dt/a over a specified early interval. Photon last scattering and baryon drag describe related but distinct endpoints; a sound horizon must state which endpoint and physical history it uses.

Later structure contains a superposition of responses to many initial perturbations. An acoustic feature therefore concerns a small statistical preference in pair separations. It does not require every galaxy to sit on one spherical wall surrounding a central object.

Bulk motions, nonlinear growth and tracer properties affect the measured feature. Here a supplied analytic spectrum and a supplied smoothing operation isolate some statistical ideas without fitting an observed cosmology or calculating its early sound horizon.

02

Define one hypothetical isotropic field

Let L be an arbitrary positive comoving length. Choose P₀ = 0.1L³, α = 0.5 and r_s = 5L. These are teaching parameters; r_s is a chosen length label, not an evaluated physical sound horizon or a value in megaparsecs.

Define j₀(z) = sin z/z with j₀(0) = 1, and choose P(k) = P₀ exp(−σ²k²/2)[1 + αj₀(kr_s)]. The wavenumber k has units of inverse length, P has units of length cubed and σ is a comoving smoothing length.

Use the Fourier convention ξ(r) = ∫₀∞ k²P(k)j₀(kr) dk/(2π²). The two-point correlation ξ is dimensionless. All analytic integrals run over the full positive wavenumber range, even though a drawing may show only a finite window.

Because |j₀| ≤ 1 and α = 0.5, this spectrum is nonnegative. It specifies a hypothetical statistically homogeneous and isotropic perturbation field. Its small variance alone does not guarantee positive physical density in every possible random realisation.

03

The transform separates a central term from an acoustic excess

For this spectrum, ξ(r) = P₀[G_σ(r) + αK_σ(r,r_s)], where G_σ(r) = (2π)⁻³ᐟ²σ⁻³ exp[−r²/(2σ²)]. The Gaussian central term and the acoustic excess are explicitly separate.

For positive r, K_σ = (2π)⁻³ᐟ²{exp[−(r−r_s)²/(2σ²)] − exp[−(r+r_s)²/(2σ²)]}/(2σrr_s). At r = 0 use its regular limit K_σ(0,r_s) = G_σ(r_s), rather than dividing by zero.

Equivalently, K_σ is the uniform angular average of G_σ(|r − r_s n|) over unit vectors n. This makes the three-dimensional geometry explicit. The maximum of the acoustic excess need not fall exactly at r_s.

The table shows E(r) = P₀αK_σ, rather than the full ξ. For σ = L, the full variance is ξ(0) ≈ 0.0063493754; for σ = 2L it is about 0.0008111062. A finite radial plot from 0 to 10L does not redefine those full-transform results.

Acoustic excess E = P₀αK_σ in the supplied correlation model; values are dimensionless and exclude the central Gaussian term
r/LE with σ = LE with σ = 2L
00.00000001183090.0000174357
30.00001432150.0000320746
40.00004813850.0000350190
50.00006349360.0000317467
60.00003209240.0000233470
80.0000004408440.00000644168
100.0000000001183090.000000697429
04

Filtering density amplitude squares the filter in power

Now prescribe a density-mode amplitude filter W(k) = exp(−Σ²k²/2), with Σ = √(3/2)L. Multiplying each density Fourier amplitude by W multiplies its power by W² = exp(−Σ²k²).

Starting from σ = L, the filtered spectrum has the same analytic form with σ_eff² = σ² + 2Σ² = 4L². Thus σ_eff = 2L. Using W only once in the power would give the wrong smoothing length.

At r = 5L, the acoustic excess decreases from about 0.0000634936 to 0.0000317467, while excess at some more distant separations increases as the feature broadens. The central Gaussian contribution also changes, so these excess values are not the full correlation.

For example, at kL = 1 the initial power is about 0.0548369L³ and the filtered power about 0.0122358L³. Their ratio is exp(−1.5), the square of the chosen amplitude filter at that wavenumber.

A chosen power spectrum and its three-dimensional acoustic excessLet L be an arbitrary positive comoving length, P₀ = 0.1 L³, α = 0.5 and r_s = 5L. These are chosen parameters, not a measured sound horizon or megaparsec calibration. With j₀(z) = sin(z)/z and j₀(0) = 1, prescribe the full nonnegative isotropic spectrum P(k) = P₀ exp(−σ²k²/2)[1+αj₀(kr_s)] for widths σ = L and 2L. k has inverse-length units, P has length-cubed units. Panel A displays only kL from zero to four with a linear P/L³ ordinate starting at zero. Solid blue circles identify σ = L and dashed orange squares σ = 2L. At kL = one the powers are approximately 0.0548369 L³ and 0.0122358 L³. The stated inverse Fourier convention is ξ(r) = integral from zero to infinity of k²P(k)j₀(kr) dk/(2π²), dimensionless. All analytic integrals are over the full positive wavenumber range; finite drawings do not truncate them. The full correlation is P₀[G_σ(r)+αK_σ(r, r_s)], where G_σ(r) = (2π)^(−3/2)σ^(−3)exp[−r²/(2σ²)]. The three-dimensional K is the uniform angular average of G_σ(|r−r_s n|) over unit vectors, not a one-dimensional Gaussian shell. For positive r it equals (2π)^(−3/2){exp[−(r−r_s)²/(2σ²)]−exp[−(r+r_s)²/(2σ²)]}/(2σ r r_s), with regular limit G_σ(r_s) at r = zero. Panel B displays only the acoustic excess E = P₀αK, multiplied by 10^5 on a linear zero-baseline ordinate, for r/L zero to ten. E remains dimensionless; multiply displayed heights by 10^(−5) to recover E. Circles and squares sample r/L = 0, 3, 4, 5, 6, 8, 10. At r = 5L the excesses are approximately 0.0000634936 and 0.0000317467. The violet dotted r_s = 5L guide marks a supplied input, not the exact peak; three-dimensional geometry shifts the maximum. The excluded central term changes too: full variances are approximately 0.0063493754 and 0.0008111062. A supplied density-amplitude filter W = exp(−Σ²k²/2), Σ² = 1.5L², multiplies power by W². Starting at σ = L, this gives σ_eff² = σ²+2Σ² = 4L², hence σ_eff = 2L and power ratio exp(−1.5) at kL = one. Small spectrum variance does not guarantee positive physical density in every realization. No measured correlation, shell around every galaxy, complete transfer function, current cosmology, calibrated nonlinear displacement law or physical sound-horizon calculation is supplied. Drag and photon last scattering require distinct endpoints in any real propagation integral.A statistical acoustic excess changes under a filterChosen P₀ = 0.1 L³, α = 0.5, r_s = 5L; L is arbitrary.P(k) = P₀ exp(−σ²k²/2) [1 + α j₀(kr_s)].j₀(z) = sin(z)/z with j₀(0) = 1.σ = Lσ = 2LA · P / L³ · supplied nonnegative spectrum, linear scale00.050.10.1501234Comoving wavenumber kL · dimensionlessB · Acoustic excess E × 10⁵ · dimensionless; not full ξ024680246810Pair separation r / L · supplied r_s guide at 5ξ(r) = ∫₀∞ k² P(k) j₀(kr) dk / (2π²); full integral, not plot cut.E = P₀ α K_σ; the central Gaussian term is excluded here.The 3D peak need not coincide with the chosen r_s = 5L.Amplitude W = exp(−Σ²k²/2) enters power as W².Σ² = 1.5 L² gives σ_eff² = L² + 2Σ² = 4L².At kL = 1: power ratio = exp(−1.5), not exp(−0.75).No measured sound horizon, physical shell or calibrated evolution.
Let L be an arbitrary positive comoving length, P₀ = 0.1 L³, α = 0.5 and r_s = 5L. These are chosen parameters, not a measured sound horizon or megaparsec calibration. With j₀(z) = sin(z)/z and j₀(0) = 1, prescribe the full nonnegative isotropic spectrum P(k) = P₀ exp(−σ²k²/2)[1+αj₀(kr_s)] for widths σ = L and 2L. k has inverse-length units, P has length-cubed units. Panel A displays only kL from zero to four with a linear P/L³ ordinate starting at zero. Solid blue circles identify σ = L and dashed orange squares σ = 2L. At kL = one the powers are approximately 0.0548369 L³ and 0.0122358 L³. The stated inverse Fourier convention is ξ(r) = integral from zero to infinity of k²P(k)j₀(kr) dk/(2π²), dimensionless. All analytic integrals are over the full positive wavenumber range; finite drawings do not truncate them. The full correlation is P₀[G_σ(r)+αK_σ(r, r_s)], where G_σ(r) = (2π)^(−3/2)σ^(−3)exp[−r²/(2σ²)]. The three-dimensional K is the uniform angular average of G_σ(|r−r_s n|) over unit vectors, not a one-dimensional Gaussian shell. For positive r it equals (2π)^(−3/2){exp[−(r−r_s)²/(2σ²)]−exp[−(r+r_s)²/(2σ²)]}/(2σ r r_s), with regular limit G_σ(r_s) at r = zero. Panel B displays only the acoustic excess E = P₀αK, multiplied by 10^5 on a linear zero-baseline ordinate, for r/L zero to ten. E remains dimensionless; multiply displayed heights by 10^(−5) to recover E. Circles and squares sample r/L = 0, 3, 4, 5, 6, 8, 10. At r = 5L the excesses are approximately 0.0000634936 and 0.0000317467. The violet dotted r_s = 5L guide marks a supplied input, not the exact peak; three-dimensional geometry shifts the maximum. The excluded central term changes too: full variances are approximately 0.0063493754 and 0.0008111062. A supplied density-amplitude filter W = exp(−Σ²k²/2), Σ² = 1.5L², multiplies power by W². Starting at σ = L, this gives σ_eff² = σ²+2Σ² = 4L², hence σ_eff = 2L and power ratio exp(−1.5) at kL = one. Small spectrum variance does not guarantee positive physical density in every realization. No measured correlation, shell around every galaxy, complete transfer function, current cosmology, calibrated nonlinear displacement law or physical sound-horizon calculation is supplied. Drag and photon last scattering require distinct endpoints in any real propagation integral. Original Ark teaching model.
05

A clean filter is only one part of real evolution

The supplied filter illustrates suppression of small-scale power and broadening of a chosen correlation feature. It is not a derivation of gravitational displacements, a calibrated nonlinear evolution law or a source of all new power produced by mode coupling.

A measured galaxy correlation additionally depends on sample selection, tracer bias, redshift-space effects, distance-coordinate choices and measurement covariance. Extracting a distance requires a physical ruler model and a likelihood, rather than reading r_s directly from this toy.

The preferred length in the early transfer process is also distinct from the exact maximum of a later, smoothed three-dimensional correlation function. Geometry, the broad spectrum and subsequent evolution matter when translating one into the other.

An acoustic imprint can link early plasma physics with later statistical structure. Establishing that link in data requires the full transfer calculation, observational controls and model comparison; the analytic example supplies no current detection, precision claim or uniquely confirmed origin of perturbations.

Worked reasoning

Recover the effective smoothing length

For initial σ = L and amplitude-filter width Σ = √(3/2)L, what is the new σ_eff?

  1. Power is quadratic in amplitude, so P_new = W²P and W² = exp(−Σ²k²).
  2. Combine exponentials: exp[−(σ² + 2Σ²)k²/2].
  3. Insert σ² = L² and 2Σ² = 3L² to obtain σ_eff² = 4L².

The effective width is 2L. Applying the amplitude filter only once to power would miss the factor of two in the added variance.

Try the reasoning

  1. Does a preferred pair separation require every galaxy to lie on one physical shell?

    Show the answer and reasoning

    No.

    The correlation describes statistical pair relationships in superposed perturbations, rather than one universal material wall.

  2. At r = 5L, is the table entry E the complete correlation ξ?

    Show the answer and reasoning

    No.

    The central contribution P₀G_σ must also be added; the table deliberately isolates the acoustic excess.

  3. Does this Gaussian amplitude filter predict all nonlinear gravitational power or an observed acoustic distance?

    Show the answer and reasoning

    No.

    The operation is supplied for teaching. Actual evolution, tracer effects, physical ruler calibration and a data likelihood require additional work.

What to carry forward

  • An early acoustic process can leave a preferred statistical separation without forming one universal galaxy shell.
  • A Fourier convention, units and the full transform distinguish power from the correlation feature.
  • A density-amplitude filter enters power squared; a prescribed smoothing model does not replace full nonlinear evolution or observation.

References

Sources supporting the explanations, physical distinctions and examples in this lesson.

  1. Large-Scale Structure of the Universe and Cosmological Perturbation TheoryF. Bernardeau; S. Colombi; E. Gaztañaga; R. ScoccimarroTeaching resourcearxiv.orgRead source
  2. On the Robustness of the Acoustic Scale in the Low-Redshift Clustering of MatterDaniel J. Eisenstein; Hee-Jong Seo; Martin WhiteResearch paperarxiv.orgRead source
  3. Baryonic Features in the Matter Transfer FunctionDaniel J. Eisenstein; Wayne HuResearch paperarxiv.orgRead source