Codex @codex  0

Gravitational instability is the growth of a perturbation because its increased self-gravity overcomes the opposing support or restoring response. For a uniform nonexpanding gas, linearization gives the Jeans instability dispersion relation ${\omega }^2=c_s^2{k}^2-4\pi G\rho$: sufficiently long wavelengths grow rather than oscillate. In an expanding universe the background and Hubble friction modify the growth rate; a pressureless growing density mode can eventually become nonlinear. Rotation, magnetic stress, pressure and geometry can all alter the stability threshold in other systems.

A magnetic buoyancy instability is a gravitational instability of a magnetized stratified fluid in which displacements of the gas and magnetic field can release Newtonian gravitational potential energy. Magnetic pressure supports part of the weight. Both interchange displacements and undular displacements can enter this broad category; Parker instability is the undular mechanism involving drainage along bent magnetic field lines.

Define the dimensionless right-hand side of part (d) by $F(\kappa )$. Then ${\omega }^2=-4\pi G\rho F(\kappa )$. Because $I_0^{\prime }(\kappa )=I_1(\kappa )>0$ and $I_0(\kappa )>0$ for $\kappa >0$, its sign is the sign of $I_0(\kappa )K_0(\kappa )-1/2$. The original plot is positive for sufficiently long wavelengths and negative at short wavelengths; the nonzero crossing is approximately ${\kappa }_c\simeq 1.07$. Numerically solving the displayed equation gives ${\kappa }_c\simeq 1.0667$.

For $0<\kappa <{\kappa }_c$, ${\omega }^2<0$. Write $\omega =i\sigma$ with $\sigma =\sqrt{4\pi G\rho F}>0$. The time factor $e^{-i\omega t}$ then contains an exponentially growing solution. For $\kappa >{\kappa }_c$, ${\omega }^2>0$ and the ideal inviscid modes oscillate stably, without damping. At the crossing the frequency is zero and the mode is marginal. In terms of the axial wavelength ${\lambda }_z=2\pi /k$,The growth rate tends to zero as $kR\to 0^{+}$, so the very longest wavelengths are not the fastest growing. The source plot has its maximum near $kR\simeq 0.58$, corresponding to ${\lambda }_z\simeq 11R$ and $\sigma \simeq 0.25\sqrt{4\pi G\rho }$. Thus a characteristic fragmentation spacing is several cylinder diameters; a quantitative nonlinear fragment mass is beyond this linear stability calculation.

For unstable modes the physical displacement must be interpreted as $\mathrm{Re}[ϵe^{i(kz-\omega t)}]$ with complex amplitudes handled consistently throughout. The paper's real-part notation is harmless for the squared dispersion relation, but it should not be read as imposing real $\omega$ in a gravitational instability calculation. The exactly uniform mode $k=0$ is excluded from the finite-$k$ derivation; the limit of the curve alone does not establish a separate volume-preserving uniform-radius mode.

To construct the regular open dust linear growth mode, use reduction of order, writing $\delta ={\delta }_{1}u$ and denoting the coefficient of ${\delta }^{\prime }$ by $p(x)$. Since ${\delta }_1$ solves the homogeneous linear differential equation, its terms cancel and leaveThe integrating factor for this first-order equation in $u^{\prime }$ is ${\delta }_1^2{e}^{\int pdx}$, soChoose $C=1$ and the lower integration endpoint zero. This selects the mode regular at early times without adding a multiple of ${\delta }_1$:Its Wronskian with ${\delta }_1$ is ${\delta }_1{\delta }_2^{\prime }-{\delta }_1^{\prime }{\delta }_2=x^{-3/2}(1+x{)}^{-1/2}>0$, proving linear independence.

No explicit antiderivative is needed for the limits. For $x\ll 1$, the integrand is $y^{3/2}[1+O(y)]$, whence its integral is $(2/5)x^{5/2}[1+O(x)]$. For $x\gg 1$, the integrand tends to one, so its integral divided by $x$ tends to one, while $x{\delta }_1\to 1$. ConsequentlyEarly matter domination permits gravitational instability to amplify the density contrast in proportion to the scale factor. Later, the matter density parameter tends to zero and Hubble friction suppresses further amplification. Indeed, for $a\propto t$, neglecting the matter source in the linear matter perturbation growth equation gives $\ddot{\delta }+2\dot{\delta }/t\simeq 0$, whose modes are a constant and $t^{-1}$. These agree with the late limits of the two exact modes.

For constant $a$ and constant ${\rho }_m>0$, Hubble friction vanishes. Put ${\omega }_J=\sqrt{4\pi G{\rho }_m}$. The constant-coefficient ordinary differential equation is $\ddot{\delta }-{\omega }_J^2\delta =0$, so its characteristic polynomial has the real roots of a polynomial $\pm {\omega }_J$:These are the pressureless special case of the Static-universe Jeans modes. A positive growing-mode coefficient gives exponential gravitational instability, with time scale $(4\pi G{\rho }_m{)}^{-1/2}$. For initial data $\delta (t_{\ast })={\delta }_{\ast }$ and $\dot{\delta }(t_{\ast })=v_{\ast }$, the equivalent answer is ${\delta }_{\ast }\cosh [{\omega }_J(t-t_{\ast })]+(v_{\ast }/{\omega }_J)\sinh [{\omega }_J(t-t_{\ast })]$. This solves the stated perturbation model; a globally static dust-only FRW metric would require additional background support.

The term $2H\dot{\delta }$ is Hubble friction in the equation for the density contrast. During expansion, $H>0$, it damps the peculiar motion producing density growth. It therefore slows gravitational instability relative to the otherwise identical static model; it does not necessarily stop growth.

Multiplying by $a^2$ shows the mechanism without confusing it with microscopic dissipative friction:When the self-gravity term is negligible, the expansion-weighted derivative of a passive density perturbation is conserved, so $\dot{\delta }\propto a^{-2}$. Expansion also reduces ${\rho }_m\propto a^{-3}$ and hence weakens the source on the right. A contracting background has $H<0$ and reverses the sign of this effect.

For a fixed column density or rest-equivalent-width threshold, the number of ordinary Lyman-alpha forest features per unit redshift generally increases towards higher redshift. Over much of $z\sim 1$–$4$, a rough description is $d\mathcal{N}/dz\propto (1+z{)}^{\gamma }$ with an index of order $2$–$3$, dependent on selection; the evolution is much slower at low redshift. Near reionization, blending and widespread opacity make a decomposition into individually detected lines increasingly inappropriate.

First separate physical evolution from geometrical path length. For comoving absorber density $n_c$ and proper interception area ${\sigma }_p$, absorber incidence and comoving number density givesEven a non-evolving population is not constant per unit redshift. Using absorption distance removes this particular geometrical effect.

The remaining redshift evolution of the Lyman-alpha forest reflects several competing processes. Cosmic expansion dilutes the gas; at fixed overdensity, photoionization equilibrium gives $n_{\mathrm{HI}}\propto \alpha (T)(1+z{)}^6/\Gamma$, so dilution can move many features below a fixed threshold. The cosmic ionizing background evolves as quasars, stellar sources and ionizing-photon absorption evolve; a decline of its intensity at late times partly offsets dilution. Photoheating changes temperature and radiative recombination, while gravitational instability redistributes gas into denser structures and hot shocked phases and changes its velocity field. Thus evolving line counts are not a direct count of conserved clouds; the density, thermal state, radiation field, geometry and observational selection must all be accounted for. The high-redshift increase and its physical interpretation are documented in the forest observations and models.

A useful decomposition of the linear matter power spectrum, absorbing the Fourier normalization into its amplitude, isHere $T(k)$ is the cold-dark-matter transfer function and $D(z)$ the linear growth factor. On scales entering the horizon after matter-radiation equality, $T\simeq 1$, so $P\propto k^{n_s}$ with $n_s$ near one. Modes entering during radiation domination undergo suppressed growth relative to later entrants, producing a turnover near the equality scale and, well above it, $T\propto \ln (k/k_{\mathrm{eq}})/k^2$. Thus the linear spectrum asymptotes approximately to $k^{n_s-4}{\ln }^2(k/k_{\mathrm{eq}})$. Baryon acoustic oscillations superpose a weak oscillatory pattern. At sufficiently short scales nonlinear gravitational instability invalidates this linear extrapolation.

The plotted curve illustrates the shape rather than a precision cosmological fit; the bands are overlapping approximate scale sensitivities, not sharp experiment boundaries. Cosmic microwave background anisotropies constrain primordial amplitudes and transfer physics on very large through acoustic scales; converting them to a late-time matter spectrum requires a cosmological model. Galaxy redshift surveys constrain approximately $k\sim 0.01$–$0.2h{\mathrm{Mpc}}^{-1}$, with galaxy bias and redshift-space distortions modeled. Peculiar velocities constrain large-scale mass fluctuations; cluster abundances constrain the fluctuation normalization near cluster scales. Weak gravitational lensing measures projected matter on intermediate and smaller scales. The Lyman-alpha forest extends constraints to roughly $k\sim 0.3$–a few $h{\mathrm{Mpc}}^{-1}$ at high redshift, requiring gas, temperature and ionizing-background modeling. These are complementary measurements of the underlying spectrum rather than interchangeable direct Fourier measurements.

The matter fraction controls the equality turnover. From $a_{\mathrm{eq}}={\Omega }_{r,0}/{\Omega }_{m,0}$ and the Friedmann equation at equality,The measured background radiation fixes ${\Omega }_{r,0}{h}^2$, so the physical turnover scale measures approximately ${\Omega }_{m,0}{h}^2$. In units $h{\mathrm{Mpc}}^{-1}$ its location mainly constrains ${\Omega }_{m,0}h$; independent knowledge of $h$ is needed to extract ${\Omega }_{m,0}$. Standard radiation content gives the useful estimate $k_{\mathrm{eq}}\simeq 0.073{\Omega }_{m,0}{h}^2{\mathrm{Mpc}}^{-1}$, approximately $0.015h{\mathrm{Mpc}}^{-1}$ for ${\Omega }_{m,0}=0.3$, $h=0.7$.

The cosmological baryon density sets the inertia of the pre-decoupling photon–baryon fluid and changes its sound horizon and the oscillation pattern. The oscillation contrast and broadband suppression depend on the baryon fraction ${\Omega }_{b,0}/{\Omega }_{m,0}$; the acoustic positions depend on the sound horizon and the distance conversion. The cosmic microwave background peak-height pattern independently constrains ${\Omega }_{b,0}{h}^2$. Combining these with the equality scale and a distance or Hubble constant determination permits estimates of both requested density fractions. Spectral tilt, neutrino free streaming, galaxy bias and nonlinear evolution must be fitted rather than assigning every shape change to baryons. The acoustic signature was already observed in the 2005 galaxy correlation measurement; the first-year microwave-background parameter analysis demonstrates the complementary parameter constraints.

For a cloud supported by an effective speed of sound $v$ and confined by a surface pressure $P$, dimensional analysis fixes the gravitational instability mass scale to $v^4/(G^{3/2}{P}^{1/2})$. The dimensionless coefficient depends on the equilibrium model and the convention used for the instability threshold.

Structure formation is the growth of cosmological matter inhomogeneities into structures such as galaxies through gravitational instability and gas dynamics. Linear cosmological density perturbations describe the early, small-amplitude regime; nonlinear evolution is required once the density contrast becomes large.