Codex @codex  0

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 }}$.

The cosmological redshift is the stretching of a freely propagating photon's wavelength between emission at cosmic time $t_e$ and reception at $t_0$ in an FLRW metric. For a comoving observer at each endpoint,Here ${\lambda }_e,{\lambda }_0$ are emitted and received wavelengths, and ${\nu }_e,{\nu }_0$ are the corresponding frequencies. A positive cosmological redshift means that the scale factor grew while the photon travelled. The name describes a shift toward longer wavelengths; it applies to light beyond the visible range as well.

For a computed example, let the scale factor double from $a(t_e)=0.5$ to $a(t_0)=1$. An emitted wavelength of $500\mathrm{nm}$ becomes $1000\mathrm{nm}$, giving $z=1000/500-1=1$. The received frequency is half the emitted frequency. This determines a ratio of spatial rulers, not a recession speed: using $cz$ as an actual velocity at $z=1$ is unjustified. The wavelength panel in the illustration under Hubble's law shows this doubling.

A short derivation follows successive wave crests along a radial null geodesic in FLRW spacetime. Their source and receiver have fixed comoving coordinates, so each crest traverses the same comoving distance:Shift both endpoint times to the next crest while keeping $\chi$ fixed. To first order in the very short periods $\delta t_e,\delta t_0$, this gives $\delta t_0/a(t_0)=\delta t_e/a(t_e)$. Locally, wavelength equals the speed of light times the crest period, yielding the displayed cosmological redshift relation. These periods are much shorter than the timescale on which the scale factor changes.

The relation isolates the smooth expansion contribution. An observed redshift can also contain a Doppler effect from source or observer peculiar velocity and a gravitational redshift. See combined kinematic and cosmological redshift for the multiplicative combination of endpoint motion and expansion. Treating every observed redshift as pure cosmological redshift can therefore misassign a nearby source's distance.

The useful insight is that a received wavelength ratio records how much the scale factor changed during the photon's journey. Turning that ratio into a distance requires an expansion history, through the Hubble parameter, and a specified convention such as luminosity distance or proper distance in cosmology; the ratio alone supplies neither.

A galaxy survey measures positions and other properties of a selected galaxy sample. A spectroscopic survey uses cosmological redshift to estimate radial distance, with corrections for peculiar velocity. Its volume, selection, galaxy bias, and sampling determine which matter power spectrum modes can be inferred.

Hubble's law states that the instantaneous rate of change of proper distance in cosmology between comoving observers is proportional to that distance:This is exact for fixed comoving distance $\chi$ in an FLRW metric: differentiating $d=a\chi$ gives $\dot{d}=\dot{a}\chi =(\dot{a}/a)d$. The dot uses cosmic time, and both distances are compared on spatial slices at the same cosmic time. The rate $\dot{d}$ is not a relative velocity measured in a single inertial frame of special relativity. Thus the instantaneous relation is not an exact observed redshift-distance formula. The associated Hubble flow describes the ideal expansion; peculiar velocities describe departures from it.

For nearby sources, the observed redshift approximately satisfieswhen source and observer peculiar velocities are small compared with $H_{0}d$, additional gravitational redshift is negligible, and the changes in the scale factor and Hubble parameter during light travel are small. Here $d$ can be taken as the present proper distance in cosmology; familiar distance conventions agree at leading order in small redshift. Briefly, $a(t_e)\simeq a(t_0)[1-H_0(t_0-t_e)]$, while $d\simeq c(t_0-t_e)$, so cosmological redshift gives $z\simeq H_{0}d/c$. At greater redshift, use an expansion history and an explicit distance convention instead.

Stipulate the illustrative Hubble constant $H_0=70\mathrm{km}{\mathrm{s}}^{-1}\mathrm{Mp}{\mathrm{c}}^{-1}$ and speed of light $c=299792.458\mathrm{km}{\mathrm{s}}^{-1}$. A $\mathrm{Mpc}$ is a million parsecs. The following are computed examples, not a catalogue of observed galaxies:

The displayed $H_{0}d$ values are exact arithmetic under the stipulated Hubble constant; the cosmological redshift estimates use the nearby approximation and become less accurate as distance grows.

Hubble's law turns a nearby observed redshift into an approximate distance indicator and makes the uniformity of cosmic expansion tangible: twice the proper distance in cosmology means twice the instantaneous separation rate. Recentring the ideal Hubble flow on another comoving observer preserves the same proportionality, as illustrated under scale factor, so recession in every direction does not select a central observer.

In his 1929 PNAS paper, Edwin Hubble presented a relation between the distances and radial velocities of extragalactic objects. His fits gave coefficients of $465$ and $513\mathrm{km}{\mathrm{s}}^{-1}\mathrm{Mp}{\mathrm{c}}^{-1}$ for 24 individual objects and nine groups respectively, and he used about $500\mathrm{km}{\mathrm{s}}^{-1}\mathrm{Mp}{\mathrm{c}}^{-1}$ as a rounded representative value. Those historical fits used the paper's distance estimates; they are distinct from the stipulated modern-scale illustration above. ^[ref].

A kinematic redshift is the relativistic Doppler effect from relative motion of emitter and receiver. With positive radial recession speed $v$, its small-velocity limit is $z_{\mathrm{kin}}=v/c+O(v^2/c^2)$. A peculiar velocity gives such a factor in addition to cosmological redshift.

The Lyman-alpha forest is the multitude of foreground Lyman-alpha absorption features in a distant quasar's spectrum, at wavelengths shorter than its own redshifted Lyman-alpha emission. The quasar supplies a bright background; intervening neutral hydrogen atoms remove photons near their local Lyman-alpha line resonance. Many dips can therefore come from one atomic transition encountered at different redshifts, rather than from many different chemical elements. Some features are instead transitions of other ions, so a dip alone does not identify Lyman-alpha absorption. This makes faint foreground gas detectable through the light it removes, even when its own emission would be difficult to see. ^[ref].

Ignoring peculiar velocity for the moment, an absorber at cosmological redshift $z_a$ puts its Lyman-alpha line atThe rest wavelength is in the ultraviolet; at $z_a=2$ it appears at $3647.01\text{A˚}$ in the optical spectrum. In a toy sightline to a quasar at $z_Q=3$, absorbers at $z_a=2,2.35,2.7$ appear at $3647.01,4072.4945,4497.979\text{A˚}$, all shortward of the quasar's Lyman-alpha emission position $4862.68\text{A˚}$. The absorbing atoms see the selected photons at their own resonance; the observed dip's redshift belongs to the foreground gas. Peculiar velocity adds a Doppler shift, so the inferred $z_a$ is not a pure distance coordinate. The rest wavelength and resonance interpretation are given in Rauch (1998), section 1.

The diagram uses invented optical depths and 26 Doppler cores; three selected foreground regions are marked. Their drawn positions and sizes are schematic. The dashed comparison imposes a uniform optical depth $\tau =2$, giving transmission $e^{-2}\simeq 0.135$ throughout the chosen foreground band. It is an idealized smooth-absorption model, not a measured spectrum or a prediction of the actual intergalactic medium.

The smooth comparison explains the Gunn-Peterson effect: neutral hydrogen distributed throughout an expanding intergalactic medium can absorb continuously over a wavelength interval, because different parts of the path resonate with different observed wavelengths. The resolved Lyman-alpha forest instead has strong transmission fluctuations and distinguishable absorption lines. These need not represent isolated compact clouds: diffuse sheets and filaments in the cosmic web can generate peaks in optical depth along one line of sight. A continuously varying gas distribution can look like a collection of lines. ^[ref]. ^[ref].

An early observation makes the identification concrete. Lynds's 1971 paper analyzed six spectrograms of the quasar 4C 05.34, with emission redshift $2.877$, taken in spring 1970 using an image-tube spectrograph on Kitt Peak's 84-inch telescope. Five covered $3200$–$7500\text{A˚}$ with approximately $10\text{A˚}$ resolution; the most informative covered about $3500$–$6000\text{A˚}$ at $5\text{A˚}$ resolution. Its reproduced spectrogram shows the strong Lyman-alpha emission flanked mostly on the short-wavelength side by numerous absorption lines. Lynds tabulated their wavelengths and estimated strengths, identified five absorption redshift systems, and proposed that many remaining unidentified lines were the same Lyman-alpha line at additional redshifts. He explicitly warned about blending, an incomplete census of weak lines, and wavelength errors around $0.5$–$1\text{A˚}$. The result established a striking spectral pattern; the paper's proposed interpretation did not measure a separate physical cloud for every dip. ^[ref].

A contrasting observation demonstrates smooth, very strong absorption. Becker and collaborators (2001) used the Echelle Spectrograph and Imager on Keck II to observe the quasar SDSSp J103027.10+052455.0 at $z_Q=6.28$. Sky lines gave an instrumental full width at half maximum of about $1.8\text{A˚}$ near $8500\text{A˚}$, or spectral resolving power $R\simeq 4700$. After sky subtraction and flux calibration, the mean continuum-normalized transmission over $8450$–$8710\text{A˚}$ was $0.0038\pm 0.0026$, consistent with zero; this corresponds to absorption redshifts approximately $5.95$–$6.16$. They identified a Gunn-Peterson trough. Crucially, resonant Lyman-alpha absorption is so sensitive that even a small neutral fraction can make transmission undetectable: the trough alone does not establish a mostly neutral intergalactic medium or uniquely date reionization. ^[ref].

What a spectrum constrains is a radiative transfer problem. For an ideal foreground screen with negligible emission into the observed beam, the intrinsic continuum-normalized transmission isHere $F$ is spectral radiative flux, $F_{\mathrm{cont}}$ is its unabsorbed continuum, $\tau$ is optical depth, $n_{\mathrm{HI}}$ is neutral hydrogen number density, and $N_{\mathrm{HI}}$ is neutral hydrogen column density along proper path length $\ell$. A Voigt profile model constrains redshift, column density and a velocity width together. A nearly black line has lost much of its sensitivity to further increases in optical depth; blending and an uncertain continuum can also admit different fitted components. Finite spectral resolving power and noise limit which widths and weak features can be measured. ^[ref]. ^[ref].

For a purely thermal Doppler core, use the velocity convention $\tau (v)\propto \exp [-(v-v_0{)}^2/b^2]$. The Maxwell-Boltzmann velocity distribution then giveswhere $k_B$ is the Boltzmann constant, $T$ is temperature and $m_{\mathrm{H}}$ is a hydrogen atom's mass. This $b$ is $\sqrt{2}$ times the Gaussian distribution's standard deviation. Higher temperature spreads the same column density over more velocities. An instrument also broadens a line: locally, for a nearly flat continuum, the recorded transmission is approximately the convolution of $\mathcal{T}$ with the normalized instrumental response, not the exponential of a convolved optical depth. A Gaussian response with full width at half maximum $c/R$ has standard deviation $c/[R\sqrt{8\log 2}]$, where $c$ is the speed of light and $R$ the spectral resolving power.

The upper model fixes the integral $\int \tau (v)dv$, hence the neutral hydrogen column density for this transition, while changing temperature. A cooler line blurred by an instrument resembles a warmer intrinsic line. The lower model has two equal Doppler cores separated by $24\mathrm{km}{\mathrm{s}}^{-1}$: reducing spectral resolving power merges them into one apparent feature. Both panels use invented profiles, Gaussian instrumental responses, and no noise or damping wings. They illustrate an inverse problem, rather than measurements of temperature or a universal relation between a dip and a cloud.

Better observations provide the inputs to that inverse problem. Kim, Cristiani and D'Odorico (2001) analyzed three quasars with the VLT's UVES spectrograph. For HE2217–2818, Table 1 records $z_Q=2.413$ and $16200\mathrm{s}$ of exposure in the $3050$–$3860\text{A˚}$ setup in October 1999. Isolated thorium-argon lamp lines measured $R\simeq 45000$; the data were resampled into $0.05\text{A˚}$ bins, with signal-to-noise ratios about $45$–$50$ per pixel over $3550$–$4050\text{A˚}$. They fitted a local polynomial continuum, fitted Voigt profiles, and compared statistics of the transmitted flux. Their analysis excluded the region near the quasar to avoid the quasar proximity effect and used the interval between its Lyman-beta and Lyman-alpha emission to reduce confusion with other Lyman-series absorption. Features below $N_{\mathrm{HI}}=10^{12.5}\mathrm{c}{\mathrm{m}}^{-2}$ became confused with noise, setting their adopted lower column density limit. Figure 1 shows the measured spectrum, fitted profiles and residuals. This was a calibrated absorption measurement, with explicit selection limits. ^[ref].

Even after instrumental correction, a width need not be purely thermal Doppler broadening. Hubble flow across extended gas, peculiar velocity structure and unresolved blends can contribute. If an isolated line's intrinsic width obeys $b^2=b_{\mathrm{th}}^2+b_{\mathrm{nt}}^2$ with an independent Gaussian nonthermal contribution $b_{\mathrm{nt}}^2\ge 0$, then $T\le m_{\mathrm{H}}{b}^2/(2k_B)$; $b=25\mathrm{km}{\mathrm{s}}^{-1}$ gives a computed upper bound about $3.79\times 10^4\mathrm{K}$. A complex blend need not satisfy this single-component model. Noise can also bias the fitted widths of weak lines. This is why the thermal cutoff of Lyman-alpha line widths, calibrated against gas models, is more informative than assigning a temperature to every dip. ^[ref].

Finally, neutral hydrogen column density is not total gas mass density. In a pure hydrogen approximation, with gas optically thin to ionizing radiation and nearly fully ionized in photoionization equilibrium, the Lyman-alpha forest photoionization scaling giveswhere $n_{\mathrm{H}}$ counts all hydrogen nuclei, $\Gamma$ is the photoionization rate, $\alpha (T)$ is the radiative recombination coefficient and $x_{\mathrm{HI}}=n_{\mathrm{HI}}/n_{\mathrm{H}}$ is the neutral fraction. At fixed temperature and $\Gamma$, doubling $n_{\mathrm{H}}$ quadruples $n_{\mathrm{HI}}$; doubling $\Gamma$ halves it. These computed examples explain both the sensitivity to density contrast and the ambiguity: stronger absorption can mean denser gas, weaker ionizing radiation, a different temperature, or more path length. The Jeans-length ionization correction for the Lyman-alpha forest adds a size model rather than an independently measured depth. Optically thick Lyman limit systems and damped Lyman-alpha systems require different treatment. ^[ref]. ^[ref].

The beauty of the Lyman-alpha forest is that a spectrum converts otherwise faint foreground gas into a sequence of markers through cosmic time. With a model of radiative transfer, velocities, ionizing radiation and the temperature-density relation of photoionized intergalactic gas, its transmission statistics probe the cosmic web and gas thermal history. A single line of sight remains one narrow sample, not a unique three-dimensional density map. ^[ref].

The projected quasar-to-cloud separation isFor the usual transverse-separation estimate, take the cloud to lie at the same line-of-sight proper distance as the quasar. Since the cloud and galaxy are gravitationally bound, the local source-to-cloud path does not acquire the redshift factor $1+z=3$ that separates that system from us. Neglect internal relative Doppler shifts and attenuation. The observer in the cloud seesThis uses the quasar's isotropic intrinsic emission at the same rest frequency; anisotropic radiation or obscuration would modify it.

An angular separation only fixes $R_{\perp }$. If $s$ is a source-to-cloud displacement toward us, the true local distance is $r=\sqrt{R_{\perp }^2+s^2}$ and the general flux is $L_{\nu }/(4\pi r^2)$. The quoted transverse value is therefore an upper estimate at fixed projected separation. Equality of cosmological redshift and gravitational binding do not by themselves measure the small internal line-of-sight displacement; a peculiar velocity can also mask it.

The Helmholtz decomposition separates the peculiar velocity into a gradient field and a transverse vector field:With decaying or periodic spatial boundary conditions this is equivalently the Fourier transform projection, for $k\ne 0$,The transverse projection operator annihilates a gradient, because $({\delta }_{ij}-k_i{k}_j/k^2)k_j=0$. Both forces in the linear cosmological Euler equation are gradients, with spatially homogeneous coefficients. Applying the transverse projection operator therefore givesSince $H=\dot{a}/a$, this is $d(a{\mathbf{v}}_{\perp })/dt=0$, and henceAt this order the scalar pressure and gravitational forces cannot generate vorticity; an initially vanishing transverse vector field remains zero. Its nonzero initial value is damped by cosmic expansion.

Define the momentum density $q_i=(\bar{\rho }+\bar{P})v_{\perp i}$, lowering its spatial label with ${\delta }_{ij}$, so the mixed stress-energy tensor obeys $T^i{}_0=q_i$ and $T^0{}_i=-q_i$. All products of $J_{ij}$ or its derivatives with $q_i$ are second order and must be dropped in linear cosmological perturbation theory.

For $\nu =0$, the three terms in stress-energy conservation areHere the trace $J_{ii}^{\prime }$ vanishes, ${\Gamma }_{00}^i={\Gamma }_{i0}^0=0$, and ${\partial }_i{q}_i=(\bar{\rho }+\bar{P}){\partial }_i{v}_{\perp }^i=0$. ThusFor $\nu =i$, spatial pressure gradients vanish because this is a purely vector cosmological perturbation. The first two terms areIn the remaining contraction, the two background-connection times momentum density terms cancel:The other terms are ${\Gamma }_{0i}^0\bar{\rho }=0$ and $-\bar{P}{\Gamma }_{ji}^j=0$. ConsequentlyExpanding this conservation law and using the background cosmological continuity equation gives the equivalent peculiar velocity equation, when $\bar{\rho }+\bar{P}\ne 0$,

The $\nu =0$ result is the background cosmological perfect-fluid continuity equation. For a comoving region of physical volume proportional to $a^3$, it is preciselyexpressing local conservation of energy including the pressure-volume work done during cosmic expansion. A transverse vector field has no velocity divergence, and the vector spatial metric tensor has no trace perturbation, so neither produces a first-order compression or density change.

The spatial result is vector-fluid momentum dilution and integrates immediately:Thus, for any specified background perfect fluid,A statement about $\mathbf{q}$ alone does not fix the velocity scaling independently of the equation of state. For constant equation-of-state parameter $w\ne -1$, the cosmological continuity equation gives $\bar{\rho }\propto a^{-3(1+w)}$ and henceFor pressureless matter, $w=0$ and ${\mathbf{v}}_{\perp }\propto a^{-1}$, exactly as in the Newtonian calculation in part (a). For radiation in cosmology, $w=1/3$ and the vector peculiar velocity remains constant, although its momentum density still falls as $a^{-4}$. More generally, for a barotropic perfect fluid with cosmological adiabatic sound speed $c_a^2={\bar{P}}^{\prime }/{\bar{\rho }}^{\prime }$, the equivalent equation is ${\mathbf{v}}_{\perp }^{\prime }+\mathcal{H}(1-3c_a^2){\mathbf{v}}_{\perp }=0$. The Newtonian damping law describes nonrelativistic inertia; it need not apply to a relativistic perfect fluid. A pure cosmological constant has $\bar{\rho }+\bar{P}=0$ and supplies no such fluid velocity degree of freedom.

The four contributions have distinct physical origins. The intrinsic photon temperature perturbation is ${\delta }_{\gamma }/4$: a local blackbody radiation energy density satisfies ${\rho }_{\gamma }\propto T^4$, so its first-order fractional perturbation is four times the fractional temperature perturbation. The endpoint cosmological Newtonian potential term ${\psi }_{\ast }$ is the gravitational redshift between the emitting plasma and the observer; escaping a negative potential well lowers the measured photon energy. Together with the intrinsic contribution it produces the ordinary Sachs-Wolfe effect. The term $-\hat{\mathbf{n}}\cdot {\mathbf{v}}_{\ast }$ is the Doppler CMB anisotropy, caused by the plasma's line-of-sight peculiar velocity at last scattering, with the source's direction convention.

Finally, $2\int d\chi {\psi }^{\prime }$ is the Integrated Sachs-Wolfe effect. A photon passing through a static potential well gains energy on entry and loses it on exit, with no net shift. When the potential evolves while the photon crosses, the two shifts fail to cancel. The factor two is the sum of the temporal and spatial scalar-potential contributions when anisotropic stress is negligible, so the two potentials are equal. The prime is a partial derivative at fixed comoving coordinate; it is not the total derivative along the light ray.

For a low-cosmological redshift galaxy sample, the intrinsic, endpoint-potential and last-scattering Doppler CMB anisotropy terms are localized near $z\simeq 1100$, whereas the galaxy density contrast is measured much nearer us. In the Limber approximation their radial kernels have no overlap, so their leading cross correlations vanish. The late Integrated Sachs-Wolfe effect samples the same nearby evolving potentials as the galaxies, giving the desired late-time CMB-matter cross correlation. The early Integrated Sachs-Wolfe effect near matter-radiation equality is likewise radially separated.

This is the leading linear, separated-kernel description intended here. The primary Cosmic microwave background anisotropy and late density contrast share primordial initial conditions, so their exact large-angle correlation need not be identically zero outside the Limber approximation. Additional secondary effects, such as the Sunyaev-Zeldovich effect, can also correlate with galaxies in real maps; they are absent from the stipulated linear expression.

Type Ia supernovae are standardizable candles. Repeated multiband photometry measures peak apparent magnitude, light-curve width and colour; spectra identify the supernova class and determine redshift, often from the host galaxy. Nearby objects calibrate the relation between light-curve shape, colour and absolute luminosity. Applying those relations produces a corrected apparent magnitude $m_{\mathrm{corr}}$, not an assumption that every event has identical luminosity. The distance modulus then givesA common uncertainty in $M$ is degenerate with the absolute distance scale, but does not erase the redshift-dependent shape of the Hubble diagram.

luminosity distance probes the expansion history. In a spatially flat model,For fixed low-redshift calibration, a negative $q_0$ increases the second-order distance relative to a decelerating model. Objects at a given redshift are therefore fainter. A Hubble diagram spanning enough redshift can distinguish this curvature from the common magnitude intercept. In matter-plus-cosmological-constant models, $q_0={\Omega }_m/2-{\Omega }_{\Lambda }$; fitting the distances can test whether $q_0<0$. The two original supernova teams found distances supporting an accelerating expansion in that framework.^[ref][ref]

The comparison requires matched rest-frame passbands and K corrections, correction for redshift time dilation in the light curves, a common photometric calibration, and explicit uncertainty propagation. Dust extinction can imitate dimming; multicolour and infrared observations test reddening and the adopted extinction law. Evolution in the progenitor population or host-galaxy dependence can alter standardization, so one compares matched populations and fits host-dependent terms where warranted. Flux-limited discovery introduces Malmquist bias and other selection effects, which must be modelled through the survey detection efficiency. peculiar velocities matter at low redshift, while gravitational lensing and intrinsic scatter broaden the distant sample. Classification and outlier tests reduce contamination by other transient types.

Independent data reductions, nearby anchors, spectroscopic checks and a fit over the full redshift range distinguish an expansion-history signal from one instrumental offset. Complementary cosmological observations help separate matter density, curvature and dark-energy parameters. The acceleration inference comes from the relative distance-redshift relation, conditional on standardization and an expansion model. Distant supernovae integrate the past expansion history; they do not directly measure the instantaneous second time derivative of today's scale factor without such modelling.

Write $\tau =t-t_0$, let $R$ be the scale factor, $H=\dot{R}/R$ the Hubble parameter, and $\mu =\sqrt{\Lambda /3}$ the asymptotic expansion rate for the positive cosmological constant. The continuity equation for a homogeneous Hubble flow gives $\rho R^3={\rho }_0$. Isotropy and the Poisson equation for Newtonian gravity give a relative acceleration proportional to displacement from any chosen origin:This uses the local relative gravitational acceleration, rather than trying to sum an absolutely convergent potential over an infinite universe. A fluid element with constant comoving coordinate $\mathbf{x}$ has $\mathbf{r}=R\mathbf{x}$. The spatially uniform pressure has no gradient, so the Euler equations for an inviscid fluid giveMultiply by $2\dot{R}$ and integrate. With $A=8\pi G{\rho }_0/3$, the resulting first integral isFor an expanding branch that reaches arbitrarily large $R$, $\dot{R}/R\to \mu$; more precisely, $H=\mu +O(R^{-2})$ and the correction to $\ln R-\mu t$ is integrable. Therefore $R\sim C{e}^{\mu t}$. Positive $\Lambda$ alone does not make every initially expanding branch escape: sufficiently large positive $\kappa$ can give a turning point or a limiting static radius. Under the stated vacuum-dominance inequalities at $R=1$, however, $A/R$ and $|\kappa |$ remain small compared with ${\mu }^2{R}^2$ for the future expanding branch, and $R\simeq e^{\mu \tau }$ is the required approximation.

To derive the linear cosmological density perturbation equation, write $\delta ={\rho }^{\prime }/\rho$ and $\mathbf{v}=H\mathbf{r}+\mathbf{u}$, where $\mathbf{u}$ is the peculiar velocity. At fixed $\mathbf{x}=\mathbf{r}/R$, the linearized continuity equation, Euler equations for an inviscid fluid and perturbation Poisson equation areTaking a divergence of the second equation and differentiating the first eliminates the peculiar velocity. For a Fourier mode of comoving wavenumber $k$, this givesThe adiabatic sound speed for the four-thirds polytropic equation of state is $c_s^2=dp/d\rho \propto {\rho }^{1/3}={\rho }_0^{1/3}/R$, hence $c_s^2=c_0^2/R$. Both the restoring pressure term and the gravitational term are proportional to $R^{-3}$. Thus, in the exponential approximation,The sign of ${\eta }^2$ is the Jeans instability distinction, but exponential expansion makes the total subsequent amplification finite.

Remove Hubble friction by setting $\epsilon =e^{-\mu \tau }y$. The standard-form equation isWhere $Q>0$ and $|Q^{\prime }|\ll Q^{3/2}$, the WKB approximation givesWhere $Q<0$, put $q=\sqrt{-Q}$; away from a turning point the corresponding WKB approximation isFor ${\eta }^2=-s^2<0$, $q=({\mu }^2+s^2{e}^{-3\mu \tau }{)}^{1/2}$ never vanishes. The exponential integral can be evaluated explicitly:Differentiate this expression, using $q^{\prime }=-3\mu (q^2-{\mu }^2)/(2q)$, to verify it. Early on, if coefficients change little over the interval considered, the two local exponents for $\epsilon$ are $-\mu \pm \sqrt{{\mu }^2+s^2}$: one increases and one decreases. At late times $q\to \mu$ and the displayed integral is $\mu \tau +O(1)$. The nominally growing WKB approximation branch of $y$ therefore becomes a constant branch of $\epsilon$, while the other branch decays like $e^{-2\mu \tau }$. The WKB approximation is not automatically accurate at all intermediate times; its stated derivative criterion must be checked.

For ${\eta }^2>0$, the early frozen-coefficient exponents are $-\mu \pm \sqrt{{\mu }^2-{\eta }^2}$. There are damped oscillations if ${\eta }^2>{\mu }^2$, a repeated-root solution $(A+B\tau )e^{-\mu \tau }$ at equality, and two real decaying local branches if $0<{\eta }^2<{\mu }^2$. In the oscillatory case $Q$ vanishes at ${\tau }_t=\ln ({\eta }^2/{\mu }^2)/(3\mu )$: the ordinary WKB approximation breaks down there, and a turning-point connection is needed. After that transition the solutions again approach a constant or decay. A positive pressure-minus-gravity restoring term does not imply indefinitely continuing oscillation; its physical frequency is redshifted away.

The limiting conclusion does not depend on the accuracy of WKB approximation. With $x=2|\eta |e^{-3\mu \tau /2}/(3\mu )$, the equation for $y$ isFor ${\eta }^2>0$ its basis is $J_{2/3}(x),Y_{2/3}(x)$, the Bessel functions of the first kind and Bessel functions of the second kind. For ${\eta }^2<0$ use $I_{2/3}(x),K_{2/3}(x)$, the Modified Bessel functions of the first kind and modified Bessel functions of the second kind. Their small-$x$ branches are proportional to $x^{2/3}$ and $x^{-2/3}$. Multiplication by $e^{-\mu \tau }$ gives the decaying and constant branches respectively. At ${\eta }^2=0$ the exact solution is $A+B{e}^{-2\mu \tau }$. Thus Bessel density modes during exponential expansion proveThe same conclusion holds in the full four-thirds perturbation equation on an exponentially escaping background: $(R^2\dot{\epsilon }{)}^{\prime }=-{\eta }^2{R}^{-1}\epsilon$. Its integrated equation has an absolutely integrable Volterra kernel, bounded by a constant times $e^{-3\mu \tau }$. The Gronwall inequality first bounds $\epsilon$; then $R^{-1}\epsilon$ is integrable and $\dot{\epsilon }=O(R^{-2})$ is integrable. A finite limiting amplitude follows without setting $H$ exactly constant. This argument requires positive $\Lambda$ and the escaping branch; it is not a claim about $\Lambda =0$ or recollapse.

The requested strict increase of that limit needs an initial-velocity qualification. For ${\eta }^2=-s^2<0$,If $\epsilon (0)>0$ and $\dot{\epsilon }(0)\ge 0$, this identity prevents a first loss of positivity and makes the amplitude strictly increase for $\tau >0$. Integrating it twice givesThis is monotone Jeans growth with positive cosmological constant. Without the extra initial condition the printed assertion is false: $\epsilon =e^{-\mu \tau }{I}_{2/3}(2s{e}^{-3\mu \tau /2}/(3\mu ))$ is positive at $\tau =0$ and tends to zero. An arbitrarily small multiple remains an infinitesimal linear cosmological density perturbation, so linearity does not remove the counterexample.

For a radial null geodesic in the FRW metric, define $\chi (r)={\int }_0^r(1-k{r}^{\prime 2}{)}^{-1/2}d{r}^{\prime }$. Two successive wave crests emitted and received at the same comoving positions obeySubtracting and keeping the first order in the wave periods gives $\Delta t_0/a(t_0)=\Delta t_e/a(t_e)$. Thus the cosmological time dilation and the cosmological redshift have the same factor:Here the emitter and observer are comoving; an additional peculiar velocity would supply a separate kinematic redshift. The derivation uses the null geodesic propagation of the crests rather than interpreting a finite cosmological distance as an ordinary Doppler effect.

The cosmological principle asserts that, after averaging over sufficiently large scales, the universe is spatially homogeneous and isotropic on constant-cosmic time slices. Homogeneity means no distinguished location; isotropy means no distinguished spatial direction. These are statistical large-scale statements, not assertions that individual galaxies, clusters or the cosmic web are uniform, and they do not require the universe to be unchanging in time.

The nearly isotropic cosmic microwave background is strong evidence for large-scale isotropy about us. Its dominant dipole is explained by our peculiar velocity; the much smaller residual fluctuations provide the seeds and record of structure formation, rather than invalidating the principle. Large-area galaxy and radio-source distributions and deep galaxy surveys support statistical isotropy and a tendency towards statistical homogeneity after smoothing beyond the strongly clustered scales. An approximately isotropic Hubble flow and broadly uniform background radiation supply complementary checks.

There is an important logical distinction: isotropy seen by one observer alone does not prove homogeneity. A spherically symmetric inhomogeneous universe centered on that observer could look isotropic. Combining the observations with the assumption that our location is not privileged, so comparable observers would see comparable large-scale isotropy, motivates the homogeneous FLRW metric. Selection effects, foregrounds, finite survey volumes and the fact that observations lie on our past light cone must be accounted for. The evidence supports a large-scale statistical principle, not exact uniformity on every scale.

Let $\delta =(\rho -\bar{\rho })/\bar{\rho }$ be the density contrast, let $\mathbf{x}$ be a comoving coordinate, and use $a_0=1$. Define $C=3H_0^2{\Omega }_0/2$, so $4\pi G\bar{\rho }=C/a^3$ for the pressureless background. The perturbation Poisson equation isFor the comoving velocity field $\mathbf{w}=\dot{\mathbf{x}}$, local mass conservation givesOne obtains this by subtracting $\dot{\bar{\rho }}+3H\bar{\rho }=0$ from ${\partial }_t\rho +3H\rho +{\nabla }_x\cdot (\rho \mathbf{w})=0$. In the linear regime it reduces to ${\partial }_t\delta =-{\nabla }_x\cdot \mathbf{w}$, the linearized cosmological continuity equation. The spatial derivatives are with respect to $\mathbf{x}$, and the time derivative in the conservative continuity equation is at fixed $\mathbf{x}$. In contrast the acceleration of an individual fluid element is a material derivative. The physical peculiar velocity is $a\mathbf{w}$.

The condition is that the velocity be the pure irrotational growing-mode velocity, with the homogeneous component in part (c) absent. It is enough to impose ${\mathbf{u}}_{\mathrm{in}}=-\nabla {\psi }_{\mathrm{in}}$ and exclude additional harmonic or transverse motions through the boundary conditions. The linear equation for $\mathbf{u}+\nabla \psi$ then keeps that difference zero. A time-independent velocity also forces the same equality in the linear motion equation when its coefficient is nonzero.

At first order the gradient can be evaluated at the initial position: shifting its argument by the first-order displacement changes the velocity only at second order. Integrate $d\mathbf{x}/db=-{\nabla }_{{\mathbf{x}}_{\mathrm{in}}}{\psi }_{\mathrm{in}}$ from $b_{\mathrm{in}}$ to $b$, givingThe trajectory is straight in comoving position as a function of growth time. The physical trajectory $a(t)\mathbf{x}(t)$ need not be straight or have constant physical peculiar velocity. This linear displacement relation is the basis of the Zeldovich approximation.

The strongest evidence for large-scale cosmological isotropy is the nearly direction-independent blackbody temperature of the cosmic microwave background. After removing the dipole associated mainly with our peculiar velocity, its residual temperature fluctuations are of order $10^{-5}$. Nearly isotropic counts of distant radio sources and galaxies, after correcting for sky coverage and selection effects, provide independent evidence. Local galaxies, clusters and voids clearly are not isotropic structures; the claim concerns a statistical description on sufficiently large scales.

Large redshift surveys show that the mean galaxy density stabilizes as increasingly large volumes are sampled, and that the excess pair correlations weaken at large separations. This supports statistical homogeneity beyond the characteristic clustering scales. The consistency of large-scale observations in different directions and the success of a common FRW metric description are additional checks.

Isotropy around one observer does not alone prove spatial homogeneity: a radially varying distribution centered on that observer could look isotropic. The Copernican principle supplies the assumption that we occupy no preferred position; comparable isotropy about typical locations then supports the cosmological principle. Surveys observe a finite past light cone, with source evolution and selection, so they cannot establish homogeneity everywhere on a simultaneous spatial slice. The evidence supports large-scale statistical isotropy and homogeneity, with substantial small-scale structure and an explicit typical-observer assumption.

Use a radial null geodesic of the FRW metric in the form $d{s}^2=-d{t}^2+a(t{)}^{2}d{\chi }^2$ along a fixed spatial direction, with $c=1$. For an emitter and observer at fixed comoving coordinates, the comoving separation isApply this to two adjacent wave crests. Since their endpoint separation is the same, subtraction to first order gives $\delta t_o/a(t_o)=\delta t_e/a(t_e)$. These intervals are the local periods measured by the comoving observers, so the received wavelength is stretched by $a(t_o)/a(t_e)$. Consequently the cosmological redshift obeyswhen $a(t_o)=1$. This relation assumes comoving emitter and observer; additional peculiar velocities supply kinematic redshift factors.

A genuine nondegenerate bounce has ${\dot{a}}_{\ast }=0$ and ${\ddot{a}}_{\ast }>0$; using ${\ddot{a}}_{\ast }\ge 0$ also includes the limiting stationary case. Put $s=a_{\ast }^{-1}=1+z_{\ast }$, and define the present matter density parameter ${\Omega }_m=8\pi G{\rho }_0/(3H_0^2)$, ${\Omega }_{\Lambda }=\Lambda /(3H_0^2)$ and ${\Omega }_k=-k/H_0^2$. The dust-plus-cosmological-constant Friedmann equation and its acceleration equation becomeAt the bounce these giveFor a past bounce with $s>1$, eliminate ${\Omega }_k$ from the first equation:Subtracting ${\Omega }_m{s}^3/2$ and multiplying by the positive denominator givesThus the matter bound from a Friedmann bounce isStrict positive bounce acceleration makes this inequality strict. The zero of the Friedmann equation also requires positive curvature for positive matter and positive $\Lambda$.

An observed source at $z=6$ has $a_e=1/7$. Since the minimum scale factor cannot exceed the scale factor at emission, $z_{\ast }\ge 6$. The denominator of the bound increases with positive $z_{\ast }$, soThe observational lower bound is much larger. Cluster mass-to-light ratios, inferred independently from motions, gas temperatures and gravitational lensing, combined with the observed luminosity density imply substantial cosmic mass density. Measurements of the cosmic baryon fraction from representative clusters, combined with the baryon cosmological density parameter inferred from Big Bang nucleosynthesis, provide another estimate: ${\Omega }_m\simeq {\Omega }_b/f_b$. Large-scale clustering and peculiar velocities provide independent dynamical checks, although their interpretation must account for galaxy bias. These methods support matter densities of order a few tenths, rather than below one percent. The dust-plus-cosmological-constant bouncing model is incompatible with both a quasar at redshift six and the stated observational evidence for ${\Omega }_m\gtrsim 0.2$.