Fact-checked by Grok 2 weeks ago

Structure formation

Structure formation in cosmology describes the evolution of the universe from a nearly uniform, hot plasma in the aftermath of the Big Bang into the complex web of galaxies, galaxy clusters, filaments, and voids observed today, primarily through the gravitational amplification of primordial density perturbations.^[1] This process is governed by general relativity and the dynamics of dark matter, baryonic matter, and radiation in an expanding universe, leading to hierarchical assembly where smaller structures merge into larger ones over cosmic time.^[1] In the prevailing Lambda cold dark matter (ΛCDM) model, initial quantum fluctuations—stretched to macroscopic scales by cosmic inflation—are the seeds for these perturbations, which grow via Jeans instability once the universe cools sufficiently for gravity to dominate over pressure.^[1] Key parameters shaping this formation include the total density parameter Ω_total ≈ 1, matter density Ω_m ≈ 0.309, and Hubble constant H_0 ≈ 67.6 km/s/Mpc, consistent with observations from the cosmic microwave background (CMB).^[2] The timeline of structure formation spans from recombination at z ≈ 1100 (about 380,000 years after the Big Bang), when photons decouple and free-streaming begins, allowing perturbations to collapse unhindered, to the present epoch at z = 0, roughly 13.8 billion years later.^[1] Dark matter, comprising about 85% of the universe's matter content, plays a dominant role by clustering first and providing gravitational wells for baryonic gas to cool, fragment, and form stars within galaxies.^[1] Baryonic processes, including hydrodynamics, radiative cooling, and feedback from supernovae and active galactic nuclei, introduce complexities on smaller scales, influencing galaxy morphologies and the suppression of structure in low-mass halos.^[3] Large-scale surveys, such as those from the Sloan Digital Sky Survey and Euclid, alongside CMB data from Planck, validate the ΛCDM framework while highlighting tensions, like the σ_8 discrepancy between early-universe and late-universe probes of clustering amplitude.^[2] Numerical N-body simulations, evolving millions to billions of particles under gravity, have been instrumental since the 1980s in predicting structure growth, power spectra, and the abundance of dark matter halos via algorithms like those in the GADGET or Enzo codes.^[3] These models incorporate cold dark matter's collisionless nature, enabling efficient collapse into cuspy profiles like the Navarro-Frenk-White (NFW) halo, though variants such as warm dark matter address potential overproduction of small-scale structures.^[1] Observational signatures include the CMB's acoustic peaks, which encode the initial perturbation spectrum, and the galaxy correlation function, tracing nonlinear evolution.^[2] Ongoing challenges encompass reconciling the Hubble tension (discrepant H_0 measurements) and exploring extensions to ΛCDM, such as modified gravity or dynamical dark energy, to better explain late-time acceleration and structure on cosmic scales.^[2]

Early Universe Foundations

Inflation and Primordial Perturbations

Cosmic inflation refers to a brief period of accelerated exponential expansion in the early universe, driven by the potential energy of a scalar field known as the inflaton. This phase resolves key issues in standard Big Bang cosmology, such as the horizon and flatness problems, by rapidly stretching quantum-scale fluctuations to cosmological sizes. The expansion occurs approximately $10^{-36} to $10^{-32} seconds after the Big Bang, corresponding to about 60 e-folds of growth at energy scales around $10^{15} GeV.^[4]^[5] During inflation, the universe approximates a quasi-de Sitter spacetime, where the Hubble parameter H is nearly constant, and the inflaton field \phi slowly rolls down its potential V(\phi). Quantum fluctuations in the inflaton field, initially sub-Hubble scale, are amplified as they exit the Hubble horizon during this expansion. These fluctuations transition from quantum to classical density perturbations through the squeezing of their state, establishing the seeds for all large-scale structure. The initial vacuum state for these modes is the Bunch-Davies vacuum, which matches the Minkowski vacuum at short wavelengths and ensures a unique, invariant choice for de Sitter-like spacetimes.^[6]^[7] The primordial perturbations generated during inflation consist primarily of scalar modes, which correspond to density contrasts \delta \rho / \rho, and tensor modes, which represent gravitational waves. Scalar perturbations are the dominant contributors to structure formation, as they source the gravitational potentials that influence matter clustering, while tensor modes provide a stochastic gravitational wave background but have negligible direct impact on density evolution. In single-field slow-roll inflation, the amplitude of scalar perturbations is given by \Delta_\zeta^2 \approx (H^2 / (8\pi^2 \epsilon M_{\rm Pl}^2)), where \epsilon = -\dot{H}/H^2 is the slow-roll parameter and M_{\rm Pl} is the Planck mass; tensor amplitudes follow \Delta_t^2 \approx (2 H^2 / (\pi^2 M_{\rm Pl}^2)).^[6]^[4] The power spectrum of these primordial curvature perturbations \zeta is nearly scale-invariant, following the Harrison-Zel'dovich form P(k) \propto k^{n_s - 4}, where n_s is the scalar spectral index. Observations from the Planck satellite confirm n_s = 0.9649 \pm 0.0042 (68% CL), indicating a slight red tilt (n_s < 1) that favors slower expansion on larger scales. This spectrum arises naturally from the quasi-de Sitter dynamics, with deviations parameterized by slow-roll parameters.^[8] The curvature perturbation \zeta is conserved on superhorizon scales for adiabatic modes, making it a robust tracer of primordial conditions. In the comoving gauge, it is defined as \zeta = -\psi - \frac{1}{3} \frac{\delta \rho / \rho}{1 + w}, where \psi is the metric perturbation, \delta \rho / \rho is the density contrast, and w = p / \rho is the equation of state. This quantity remains constant outside the horizon after inflation, providing the initial conditions for subsequent cosmological evolution.^[9]

Post-Inflation Expansion and Initial Conditions

Following the end of inflation, the universe undergoes reheating, during which the oscillating inflaton field decays into Standard Model particles, thereby populating the cosmos with a hot plasma of relativistic matter and marking the onset of the radiation-dominated era. This process typically occurs over a brief period around $10^{-32} seconds after the big bang, converting the vacuum energy of the inflaton into the thermal energy density of radiation and particles. The efficiency and details of reheating depend on the inflaton potential and coupling strengths, but it establishes the initial thermal bath necessary for subsequent cosmological evolution.^[10] The initial conditions for density perturbations post-reheating are dominated by adiabatic modes generated during inflation, where the curvature perturbation \zeta remains conserved on superhorizon scales and evolves into metric perturbations in gauges suitable for post-inflation dynamics, such as the Newtonian gauge. In this framework, \zeta sets the primordial amplitude for scalar perturbations, with nearly scale-invariant power ensuring uniformity on large scales while allowing small-scale variations that seed structure formation. These adiabatic perturbations imply that relative entropy fluctuations between components like radiation and matter are negligible, providing a coherent starting point for the coupled evolution of metric and fluid perturbations during the early expansion. During the radiation-dominated era, the primordial power spectrum P(k) is modified by the transfer function T(k), which accounts for the suppression of perturbations due to free-streaming of relativistic particles (such as neutrinos) and Silk damping from photon diffusion in the baryon-photon plasma. For large scales (k \to 0), T(k) \approx 1, preserving the primordial spectrum, whereas on small scales, free-streaming and damping lead to T(k) \propto k^{-2} or stronger suppression, limiting power transfer to subhorizon modes. This scale-dependent filtering shapes the initial conditions for matter clustering once non-relativistic components dominate. The dynamics of metric perturbations in the Newtonian gauge during this early phase are governed by the constraint equation

k^2 \Phi + 3\mathcal{H} (\Phi' + \mathcal{H} \Phi) = -4\pi G a^2 \delta \rho,

which relates the gravitational potential \Phi to the density contrast \delta \rho under the influence of Hubble expansion \mathcal{H}, highlighting the competition between gravitational clustering and cosmic dilution. As perturbation modes evolve, those larger than the horizon (k < aH) remain frozen, with \Phi constant, while subhorizon modes (k > aH) enter the horizon and begin oscillating as acoustic waves in the relativistic plasma, preventing significant growth until the matter era.^[11]

Linear Growth Regime

Density Perturbations in an Expanding Universe

In the standard model of cosmology, small initial density perturbations in the early universe serve as the seeds for large-scale structure formation, evolving under gravitational instability within the framework of general relativity. These perturbations, arising from quantum fluctuations during cosmic inflation, are characterized by the density contrast \delta = \frac{\delta \rho}{\rho}, where \rho is the background density and \delta \rho is the perturbation. In a homogeneous and isotropic expanding universe, their linear growth is governed by gravitational attraction competing with expansion and, for collisional fluids, pressure support. Gravitational instability drives the amplification of overdensities on scales larger than the Jeans length, while collisionless components like dark matter experience unhindered growth due to the absence of pressure.^[12] The background geometry of the universe is described by the Friedmann-Lemaître-Robertson-Walker (FLRW) metric, which assumes spatial homogeneity and isotropy:

ds^2 = -dt^2 + a(t)^2 \left[ \frac{dr^2}{1 - kr^2} + r^2 d\Omega^2 \right],

where a(t) is the scale factor, k is the curvature parameter (k = 0, +1, -1), and d\Omega^2 = d\theta^2 + \sin^2\theta \, d\phi^2. This metric encapsulates the expansion history, with a(t) evolving according to the Friedmann equations derived from Einstein's field equations. For small perturbations around this background, the dynamics on sub-horizon scales (k \gg H a, where H = \dot{a}/a is the physical Hubble parameter) can be approximated using the Newtonian limit, valid for non-relativistic matter velocities much less than the speed of light.^[12] In the Newtonian gauge, the evolution of scalar perturbations is captured by the linearized fluid equations in comoving coordinates. The continuity equation relates the density perturbation to the peculiar velocity \mathbf{v}:

\dot{\delta} + \frac{1}{a} \nabla \cdot \mathbf{v} = 0.

The Euler equation describes the acceleration due to the gravitational potential \Phi and Hubble drag:

\dot{\mathbf{v}} + H \mathbf{v} + \frac{1}{a} \nabla \Phi = 0,

where primes denote derivatives with respect to conformal time \eta (with dt = a d\eta), but here dots indicate physical time derivatives for clarity. The Poisson equation relates the potential to the density perturbation:

\nabla^2 \Phi = 4\pi G \rho \delta a^2.

For collisionless fluids, such as dark matter, these equations reveal pure gravitational instability without pressure terms, leading to growth on all sub-horizon scales. In contrast, for collisional baryonic gas, an additional pressure term in the Euler equation introduces the Jeans instability, suppressing growth below the Jeans wavenumber k_J = \sqrt{4\pi G \rho / c_s^2}, where c_s is the sound speed; perturbations with k > k_J oscillate rather than grow.^[12] Combining these equations yields a second-order differential equation for \delta:

\ddot{\delta} + 2H \dot{\delta} - 4\pi G \rho \delta = 0,

valid in the sub-horizon limit for pressureless matter. During matter domination, where \rho \propto a^{-3} and H^2 = 8\pi G \rho / 3, the solutions consist of a growing mode \delta_+ \propto a(t) and a decaying mode \delta_- \propto a^{-3/2}. The growing mode dominates at late times, amplifying initial perturbations proportionally to the expansion, setting the stage for structure formation once \delta \sim 1. Initial conditions from inflation provide \delta \sim 10^{-5} at horizon entry.^[13]^[12] In the late universe, dark energy modifies this growth through accelerated expansion. In the \LambdaCDM model, the growing mode is suppressed relative to pure matter domination, with the growth rate f = d \ln \delta / d \ln a approximated as f \approx \Omega_m^{0.55}, where \Omega_m is the present-day matter density parameter. This parametrization, accurate to better than 0.5% for \Omega_m \in [0.2, 1], reflects the reduced effectiveness of gravity in pulling matter together amid faster expansion, leading to slower structure growth at low redshifts.^[14]

Growth Factor and Power Spectrum Evolution

The linear growth factor D(a) quantifies the evolution of density perturbations \delta(\mathbf{k}, a) = D(a) \delta(\mathbf{k}, a_{\rm ini}) in the linear regime, where a is the cosmic scale factor. It is conventionally normalized such that D(a=1) = 1 at the present epoch. The growth factor satisfies the second-order differential equation \ddot{D} + 2 H \dot{D} - \frac{3}{2} \Omega_m \frac{H_0^2}{a^3} D = 0, derived from the perturbed Friedmann equations in an expanding universe, where dots denote derivatives with respect to cosmic time t, H = \dot{a}/a is the Hubble parameter, H_0 is its present-day value, and \Omega_m is the present-day matter density parameter.^[15]^[16] In the matter-dominated era, where \Omega_m \approx 1 and the cosmological constant is negligible, the growing mode solution approximates D(a) \propto a, reflecting proportional growth with the expanding scale factor. This approximation holds from matter-radiation equality until the onset of dark energy domination, providing a simple benchmark for perturbation amplification.^[15] The linear matter power spectrum P(k, a), which encodes the statistical distribution of density fluctuations as a function of comoving wavenumber k, evolves as P(k, a) = D^2(a) T^2(k) P_{\rm prim}(k), where T(k) is the transfer function that accounts for sub-horizon evolution after horizon entry, and P_{\rm prim}(k) is the primordial power spectrum set during inflation. This scaling ensures that the power spectrum amplitude grows quadratically with the growth factor in the linear regime. The transfer function, computed numerically from Boltzmann codes, briefly modulates the spectrum but is referenced here only for its role in shaping P(k, a).^[17] Across cosmic eras, the growth factor exhibits distinct behaviors: during radiation domination, dark matter perturbations experience logarithmic growth D(a) \propto \ln a for sub-horizon modes due to the Mészáros effect, where radiation pressure suppresses clustering. In the matter-dominated era, growth accelerates linearly with a, with baryonic contributions enhancing the total matter clustering after recombination as baryons settle into dark matter potential wells, increasing the effective \Omega_m.^[17] The dimensionless power spectrum \Delta^2(k) = (k^3 / 2\pi^2) P(k) highlights the variance per logarithmic interval in k, revealing a characteristic turnover at the equality scale k_{\rm eq} \approx 0.01 \, h \, \rm Mpc^{-1}, corresponding to modes entering the horizon at matter-radiation equality. This turnover marks the transition from scale-invariant primordial fluctuations to damped small-scale power due to free-streaming and Silk damping.^[17] The variance \sigma_8, defined as the root-mean-square density fluctuation smoothed over $8 \, h^{-1} \, \rm Mpc spheres, serves as a key normalization for the overall growth, with \sigma_8 = \sqrt{(2\pi^2 / k^3) \int_0^\infty \Delta^2(k) W^2(k R) \, dk} where W is the top-hat window function and R = 8 \, h^{-1} \, \rm Mpc. Observations from the cosmic microwave background constrain \sigma_8 \approx 0.81, anchoring the amplitude of structure formation to the present day.^[17]

Recombination Era

Physics of Recombination

Recombination in the early universe marks the epoch when the plasma of free electrons and protons transitioned to neutral hydrogen atoms, decoupling photons from baryons around redshift z \approx 1100. This process is governed by the Saha ionization equation, which balances the rates of ionization and recombination under thermal equilibrium assumptions:

\frac{x_e^2}{1 - x_e} = \frac{1}{n_b} \left( \frac{2\pi m_e k T}{h^2} \right)^{3/2} \exp\left( -\frac{I}{k T} \right),

where x_e is the electron ionization fraction, n_b is the baryon number density, m_e is the electron mass, k is Boltzmann's constant, T is the temperature, h is Planck's constant, and I = 13.6 eV is the hydrogen ionization energy. This equation predicts the equilibrium ionization state, but in the expanding universe, deviations arise due to non-equilibrium effects.^[18] The timeline of recombination begins around z \approx 1500 (corresponding to T \approx 4000 K), when the universe cools sufficiently for electrons to bind with protons, and proceeds rapidly but incompletely, completing the bulk by z \approx 800. The expansion of the universe causes a freeze-out, where the recombination rate lags behind the decreasing density and temperature, preventing full equilibrium and leaving a residual ionized fraction. This freeze-out is captured in detailed rate equations that account for the cosmological redshift, shifting the effective recombination dynamics away from simple Saha predictions.^[18] Direct recombination to the hydrogen ground state (1s) is inefficient because the emitted Lyman-alpha photons (n=2 to n=1) are resonantly scattered, maintaining coupling between matter and radiation. Instead, the two-photon decay from the 2s excited state to the 1s ground state provides a key channel for decoupling, allowing photons to escape without immediate reabsorption. High-lying excited states (n \gg 2) contribute to residual ionization, as thermal photons can ionize these loosely bound electrons, sustaining a small electron fraction post-recombination of x_e \approx 10^{-4}.^[18] The Thomson scattering cross-section, \sigma_T = \frac{8\pi}{3} \left( \frac{e^2}{m_e c^2} \right)^2, quantifies the interaction probability between photons and free electrons, determining the thickness of the last scattering surface during recombination. This surface, where most photons last interact with the plasma, spans a width set by the scattering rate and expansion, imprinting key scales in the density field. Notably, the sound horizon at recombination, r_s \approx 150 Mpc, arises from acoustic oscillations in the baryon-photon fluid prior to decoupling, where pressure support drives waves that freeze in as a characteristic scale in the matter power spectrum upon recombination. These baryon acoustic oscillations (BAO) originate from the tight coupling regime, with the horizon size reflecting the integrated sound speed up to z \approx 1100.^[18]^[19]

Transition to Matter-Dominated Era

The epoch of radiation-matter equality, occurring at a redshift z_{\rm eq} \approx 3400, marks the transition where the energy density of matter becomes comparable to that of radiation in the expanding universe.^[20] At this point, the scale factor a_{\rm eq} \propto \Omega_r / \Omega_m \approx 1/3600, beyond which matter begins to dominate the expansion dynamics.^[21] Prior to equality, density perturbations in matter grow logarithmically due to the radiation-dominated background, but afterward, gravitational instability allows for linear growth, setting the stage for the amplification of initial fluctuations into the large-scale structure observed today. Following recombination at z \approx 1090, when the optical depth to Thomson scattering drops to \tau = 1 and the universe becomes transparent, the plasma of ionized hydrogen and helium transitions to a neutral gas primarily composed of atomic hydrogen and helium.^[22] Photons, now decoupled from baryons, begin to free-stream without significant interactions, removing the pressure support that previously suppressed gravitational clustering on sub-horizon scales.^[23] This shift enables matter perturbations, which were largely frozen during the tight-coupling phase of recombination, to resume growth under gravity alone. A key feature of this transition is the Silk damping effect, arising from the random walk of photons during the final stages of decoupling, which erases density perturbations on small scales through diffusion. The damping scale is characterized by k_d \approx 0.2 \, h \, \rm Mpc^{-1}, suppressing power in the matter spectrum for comoving scales smaller than approximately 30 Mpc. On larger scales, unaffected by this damping, the density contrast \delta evolves as \delta \propto a in the matter-dominated era, with cold dark matter providing the primary gravitational potential wells that drive the clustering of both dark matter and baryons.^[24]^[25] This linear growth regime persists until nonlinear effects become prominent at lower redshifts, laying the foundation for hierarchical structure formation.

Dark Matter Structure Development

Linear Dark Matter Clustering

Dark matter in cosmological structure formation is modeled as collisionless, non-relativistic particles that dominate the matter content of the universe, with leading candidates including weakly interacting massive particles (WIMPs) and axions.^[26] These particles decouple early from the thermal bath, behaving as pressureless dust with negligible self-interaction, which permits unrestricted gravitational infall and efficient clustering on all relevant scales.^[27] In the cold dark matter (CDM) paradigm, the particles exhibit very low velocity dispersion, with typical speeds much less than the speed of light (v \ll c), contrasting with hot dark matter candidates like massive neutrinos that have relativistic velocities at early times and suppress small-scale structure formation. This cold nature ensures that primordial density perturbations in dark matter can grow freely via gravity without thermal pressure opposing collapse, seeding the hierarchical buildup of cosmic structures from large to small scales. The statistical distribution of dark matter density fluctuations is characterized by its power spectrum, which in the linear regime evolves as P_{\rm dm}(k, a) \approx D^2(a) P_{\rm lin}(k), where D(a) is the linear growth factor normalized to unity at the present scale factor a = 1, k is the comoving wavenumber, and P_{\rm lin}(k) is the initial linear power spectrum derived from primordial perturbations. For dark matter itself, the bias parameter b = 1, reflecting its direct role as the underlying matter tracer without additional biasing effects seen in galaxies. This scaling captures how fluctuations amplify proportionally to the square of the growth factor during the matter-dominated era, preserving the shape of the initial spectrum on large scales until nonlinear effects emerge. To extend beyond the strict Eulerian linear perturbation theory, the Zeldovich approximation provides a first-order Lagrangian description of dark matter particle trajectories, where the displacement field \psi maps initial Lagrangian positions \mathbf{q} to Eulerian positions \mathbf{x} = \mathbf{q} + \psi(\mathbf{q}, a), with \psi = -\nabla \Phi / (a H f) for the growing mode.^[28] Here, \Phi is the gravitational potential, H is the Hubble parameter, and f = d \ln D / d \ln a is the logarithmic growth rate. This approximation analytically predicts the formation of sheet-like structures (pancakes) by following ballistic motion of particles along initial velocity fields, offering insights into the transition to mildly nonlinear regimes on scales beyond strict linearity.^[28] On large scales, dark matter clustering induces coherent peculiar velocities that drive bulk flows, given in linear theory by \mathbf{v} = -f H a \nabla^{-1} \delta_{\rm dm}, where \delta_{\rm dm} is the dark matter density contrast and the inverse Laplacian extracts the velocity potential.^[29] The growth rate f is well-approximated by f \approx \Omega_m^{0.55} in flat \LambdaCDM cosmologies, where \Omega_m is the present matter density parameter, enabling quantitative predictions for velocity statistics. These velocities manifest observationally as redshift-space distortions, elongating structures along the line of sight in galaxy surveys and providing a probe of the underlying linear dark matter distribution.^[29]

Nonlinear Halo Formation and Collapse

As density perturbations in the dark matter component grow beyond the linear regime, where the linear growth factor no longer accurately describes their evolution, gravitational instability drives the formation of bound structures known as dark matter halos. This nonlinear phase begins when the density contrast δ exceeds unity on small scales, leading to gravitational collapse and the emergence of virialized halos that serve as the gravitational wells for galaxy formation. The transition is characterized by the breakdown of perturbation theory, requiring semi-analytic models and simulations to predict halo properties.^[30] The spherical collapse model provides a foundational semi-analytic description of this nonlinear evolution, assuming a spherically symmetric top-hat overdensity in an otherwise homogeneous Einstein-de Sitter universe. In this model, the overdense region expands with the Hubble flow but decelerates due to self-gravity, reaching a maximum expansion (turnaround) when the linearly extrapolated density contrast δ_ta ≈ 1.06, after which it collapses to form a singular density at the center. The full nonlinear collapse occurs when the linear extrapolation of the density contrast reaches δ_coll ≈ 1.686, marking the threshold for halo formation; this critical value is independent of scale in the Einstein-de Sitter approximation but adjusts slightly in more general cosmologies.^[31]^[30] Building on the spherical collapse threshold, the Press-Schechter formalism offers a statistical framework for predicting the abundance of halos of different masses at a given epoch. It assumes that the fraction of mass in regions exceeding the critical density contrast δ_c = 1.686 in the linear regime corresponds to the collapsed fraction, leading to a halo mass function of the form

\frac{dn}{dM} \propto \frac{\delta_c}{\sigma(M)} \exp\left(-\frac{\delta_c^2}{2\sigma^2(M)}\right) \frac{1}{M},

where σ(M) is the root-mean-square density fluctuation smoothed over a mass scale M, derived from the linear power spectrum. This exponential cutoff ensures rarity of massive halos, while the prefactor captures the abundance of low-mass objects, providing a first-order prediction for the halo distribution that matches early simulation results reasonably well despite its simplicity.^[30] The extended Press-Schechter formalism and its excursion set theory generalization address limitations of the original approach by incorporating the hierarchical merging history of halos. In this framework, halo formation is modeled as a random walk in the space of density variance, where trajectories cross a constant barrier at δ_c to form halos; merger trees are constructed by considering upcrossings of moving barriers that account for progenitor distributions. This enables predictions of halo assembly rates and conditional mass functions, revealing that typical halos form through successive mergers of smaller progenitors rather than monolithic collapse. Once collapsed, dark matter halos achieve virial equilibrium, where the virial theorem dictates that twice the kinetic energy equals the absolute value of the potential energy: 2K + W = 0. For a halo of mass M, this balance defines a virial radius r_vir enclosing an average overdensity Δ relative to the critical density ρ_c, approximated as

r_\mathrm{vir} \approx \left( \frac{3M}{4\pi \Delta \rho_c} \right)^{1/3},

with Δ ≈ 200 commonly adopted for Navarro-Frenk-White (NFW) profiles in low-redshift clusters, reflecting the point where infall transitions to virialization.^[32] High-resolution N-body simulations reveal that the equilibrium density profiles of dark matter halos are well-fitted by the Navarro-Frenk-White (NFW) form:

\rho(r) = \frac{\rho_s}{(r/r_s)(1 + r/r_s)^2},

characterized by a scale density ρ_s and scale radius r_s, with a cuspy inner slope of -1 and an outer slope of -3. This universal profile emerges from hierarchical clustering across a wide range of masses and redshifts, underscoring the self-similar nature of nonlinear gravitational collapse in cold dark matter cosmologies.^[32]

Baryonic Gas Dynamics

Post-Recombination Gas Evolution

Following recombination at z ≈ 1100, the baryonic gas decouples from the photons and becomes predominantly neutral, allowing it to respond freely to gravitational potentials dominated by dark matter. The gas falls into these dark matter wells, tracing the underlying density perturbations on large scales with a linear bias b_b \approx 1 + \frac{3}{5} \frac{\Omega_b}{\Omega_m}, where \Omega_b and \Omega_m are the present-day baryon and total matter density parameters, respectively. This infall amplifies the initial primordial fluctuations, initiating the hierarchical buildup of cosmic structure as the universe transitions to matter domination.^[33] The thermal evolution of the gas is initially governed by residual Compton scattering with cosmic microwave background (CMB) photons, maintaining the gas temperature close to the CMB temperature, T_\mathrm{gas} \approx T_\mathrm{CMB}, until approximately z ≈ 200 when the coupling becomes negligible. Beyond this redshift, the gas undergoes adiabatic cooling due to cosmic expansion, with its temperature scaling as T \propto 1/a^2, where a is the scale factor; this results in a cooler intergalactic medium (IGM) compared to the CMB, reaching temperatures around 10 K by z = 10.^[33] The dynamics of the baryonic gas are described by the standard set of hydrodynamic equations in an expanding universe: the continuity equation for mass conservation, the Euler equation for momentum, which includes a pressure support term -\frac{1}{\rho} \nabla P (where ρ is the gas density and P is the pressure), and the Poisson equation for the gravitational potential sourced primarily by dark matter. These equations couple the gas motion to dark matter gravity while accounting for the gas's thermal pressure, which suppresses collapse on small scales until cooling processes become efficient.^[33] The earliest collapsed objects, known as minihalos, form through this gas infall at redshifts z ≈ 100–1000, with characteristic masses of $10^5–$10^6 \, M_\odot, where molecular hydrogen (H₂) line cooling enables fragmentation and collapse within dark matter halos. These structures represent the precursors to the first stars but are limited in number and impact due to their small size.^[33] Later, at z ≈ 6–10, the onset of reionization is triggered by ultraviolet photons from these early stars (and subsequent generations), ionizing the neutral IGM and raising its temperature to 10⁴ K, which modifies the gas dynamics by reducing neutral hydrogen density and enhancing photoheating effects. Recent James Webb Space Telescope (JWST) observations suggest that faint active galactic nuclei may contribute substantially to the ionizing photon budget, potentially starting reionization earlier than previously thought and highlighting tensions in standard models.^[33]^[34]

Cooling Processes and Feedback Mechanisms

In the early universe, baryonic gas within dark matter halos cools primarily through atomic processes once the virial temperature exceeds approximately 8000 K, enabling gravitational collapse in halos with masses greater than about $10^8 \, M_\odot at redshifts z > 10. The dominant cooling channel is collisional excitation of neutral hydrogen atoms followed by radiative decay via the Lyman-alpha line at 1216 Å, supplemented by recombination radiation and bremsstrahlung from ionized gas. For smaller halos below this threshold, molecular hydrogen (H_2) lines provide the primary cooling mechanism, though atomic cooling becomes efficient in more massive structures where the gas can reach higher densities and temperatures. The overall cooling rate follows \Lambda \approx n^2 \alpha(T), where n is the hydrogen number density and \alpha(T) is the temperature-dependent cooling function; for atomic hydrogen, \alpha(T) \propto T^{-0.7} in the regime around $10^4 K, reflecting the balance between excitation and radiative de-excitation. Subsequent metal enrichment from the first supernovae dramatically enhances cooling efficiency, particularly at lower temperatures below 8000 K, by introducing elements like carbon and oxygen that enable fine-structure line emission. These Population III supernovae eject metals into the surrounding intergalactic medium (IGM) and nearby halos, with the enriched gas exhibiting a cooling rate \Lambda_\mathrm{metal} \propto Z n^2 T^{0.5} at temperatures T \lesssim 10^4 K, where Z is the metallicity.^[35] This enhancement lowers the temperature floor for collapse, promoting fragmentation and the formation of lower-mass Population II stars compared to the metal-free case. Feedback mechanisms play a crucial role in regulating these cooling processes and modulating structure growth. Radiative feedback, primarily from ultraviolet photons during reionization, photoheats the IGM to temperatures around $10^4 K, suppressing cooling and collapse in low-mass halos by increasing the thermal pressure. Mechanical feedback arises from supernova outflows, which inject kinetic energy and heat the surrounding gas, driving outflows that enrich and heat the IGM while preventing excessive gas accretion onto galaxies.^[35] In more massive systems, active galactic nuclei (AGN) provide powerful mechanical feedback through relativistic jets and winds, expelling gas on kiloparsec scales and quenching star formation in host galaxies.^[36] The Jeans mass, defined as M_J \propto T^{3/2} / \rho^{1/2} where \rho is the gas density, evolves significantly under these influences, setting a dynamic threshold for gravitational instability. During reionization, photoheating raises the gas temperature, increasing M_J by factors of 10–100 and delaying star formation in halos below $10^8 \, M_\odot until larger structures can form. Radiative transfer effects further shape ionized regions around early sources, as described by the Strömgren sphere radius R_s = \left( \frac{3 \dot{N}_\gamma}{4\pi \alpha_B n_H^2} \right)^{1/3}, where \dot{N}_\gamma is the ionizing photon emission rate and \alpha_B is the case-B recombination coefficient. For the first massive stars or miniquasars, this yields R_s \approx 10 kpc, carving out discrete H II regions that contribute to patchy reionization and influence nearby gas dynamics.^[37]

Observational Evidence

Large-Scale Structure Surveys

Large-scale structure surveys have been instrumental in mapping the distribution of galaxies and galaxy clusters, providing empirical evidence for the hierarchical formation of cosmic structures predicted by the ΛCDM model. These surveys measure the clustering of matter through redshift-space distortions, baryon acoustic oscillations (BAO), and weak gravitational lensing, tracing the evolution of density perturbations from redshifts z ≈ 0 to z ≈ 3.5. By compiling spectroscopic redshifts for millions of objects, they enable precise determinations of the matter power spectrum and the growth rate of structure, offering constraints on cosmological parameters such as the matter density Ω_m and the amplitude of fluctuations σ_8.^[38] The Sloan Digital Sky Survey (SDSS) stands as a foundational effort in this domain, having mapped the three-dimensional positions of approximately 10^6 galaxies across a significant portion of the sky to probe large-scale structure on scales up to 100 Mpc. Through its Baryon Oscillation Spectroscopic Survey (BOSS) phase, SDSS measured the galaxy power spectrum, confirming the acoustic peak feature corresponding to the BAO scale imprinted at recombination, which aligns with ΛCDM predictions for the expansion history and matter clustering. These measurements, derived from luminous red galaxies and quasars at redshifts z < 1, validated the standard model's success in describing the observed filamentary web of the cosmic large-scale structure.^[38]^[39] Building on SDSS, the extended BOSS (eBOSS) component of SDSS-IV extended these measurements to higher redshifts, achieving precision BAO detections using emission-line galaxies, quasars, and Lyman-α forest absorbers over 0.5 < z < 3.5. The survey's analysis yielded robust constraints on the dimensionless combinations H(z) r_d and D_M(z)/r_d at intermediate redshifts, enabling tight bounds on H_0 ≈ 67-70 km/s/Mpc and Ω_m ≈ 0.3 when combined with other probes. These results refined the distance-redshift relation and tested the consistency of dark energy dominance in the late universe.^[40] The Dark Energy Spectroscopic Instrument (DESI) collaboration released Data Release 2 (DR2) results in March 2025, presenting baryon acoustic oscillation measurements from over 14 million galaxies and quasars across 0.1 < z < 2.3. These data achieve 0.5% precision on the BAO scale, constraining the expansion history and revealing hints of evolving dark energy, while contributing to the σ_8 tension debate with lower clustering amplitudes than CMB predictions.^[41] The Dark Energy Survey (DES) complemented spectroscopic efforts with wide-field imaging, employing weak lensing tomography to dissect the matter distribution across redshift bins up to z ≈ 1. DES Year 3 results integrated galaxy clustering and cosmic shear measurements over 5000 deg², quantifying the growth rate parameter fσ_8(z) — the product of the linear growth factor f and the smoothed density contrast σ_8 — which probes the rate of structure amplification under gravity. This approach also revealed the cosmic web's topology, including the detection of large voids (diameters > 100 Mpc) and prominent filaments, highlighting deviations from uniformity on scales of 10-100 Mpc/h.^[42] Galaxy cluster abundance provides a complementary probe of nonlinear structure growth, with surveys targeting massive halos (M > 10^{14} M_⊙) via X-ray emission and the Sunyaev-Zel'dovich (SZ) effect. The eROSITA telescope, launched in 2019, has detected thousands of clusters in its all-sky survey, enabling measurements of the cluster number counts dN/dz that evolve with redshift according to the halo mass function, approximately following dN/dz ∝ exp(-E(z)^2 / 2σ_8^2) in the high-mass limit, where E(z) is the normalized expansion rate. These data test nonlinear evolution and constrain σ_8 ≈ 0.8 with percent-level precision when calibrated against weak lensing masses. Similarly, Planck's SZ cluster catalog, comprising over 1600 clusters out to z ≈ 1, used thermal SZ decrement to derive abundance trends sensitive to σ_8, confirming growth suppression consistent with ΛCDM but highlighting mild tensions with primary CMB inferences. As of 2025, the Euclid space mission's early data releases have advanced 3D clustering analyses, leveraging spectroscopic and photometric redshifts for millions of galaxies over 15,000 deg² to map baryon acoustic features and redshift-space distortions with unprecedented volume coverage. The Quick Data Release 1 (Q1) in March 2025 provided initial 3D maps from 63 deg², enabling early constraints on σ_8 with sub-1% precision when extrapolated to full survey forecasts, significantly tightening bounds on matter clustering and dark energy equation-of-state parameters. These advances, combined with theoretical power spectrum models, underscore the ongoing refinement of structure formation paradigms.^[43]^[44]

Cosmic Microwave Background Anisotropies

The cosmic microwave background (CMB) anisotropies serve as a snapshot of the universe at the epoch of recombination, approximately 380,000 years after the Big Bang, when photons decoupled from the baryon-photon plasma and the universe became transparent. These tiny temperature fluctuations, on the order of $10^{-5} relative to the mean CMB temperature of 2.725 K, encode the initial density perturbations that seeded large-scale structure formation. The power spectrum of these anisotropies, characterized by angular scales via multipole moments \ell, reveals the evolution of these perturbations through gravitational instability, acoustic oscillations, and diffusion damping up to the last scattering surface. Observations from satellites like COBE, WMAP, and Planck have mapped these anisotropies with increasing precision, providing stringent constraints on cosmological parameters and confirming the adiabatic, nearly scale-invariant primordial spectrum predicted by inflation. On large angular scales (\ell \lesssim 10), corresponding to superhorizon modes at recombination, the dominant contribution to temperature anisotropies arises from the Sachs-Wolfe effect. This effect describes how photons climbing out of gravitational potential wells lose energy, leading to a temperature perturbation \Delta T / T = \frac{1}{3} \Phi, where \Phi is the primordial gravitational potential at last scattering, plus smaller integrated terms along the line of sight. For modes larger than the horizon at recombination, the ordinary Sachs-Wolfe term dominates, directly linking observed CMB fluctuations to the primordial potential depth, which is tied to the amplitude of scalar perturbations A_s \approx 2.1 \times 10^{-9}. This effect is crucial for probing the initial conditions of structure formation, as it imprints the gravitational redshift from early density contrasts without significant post-recombination evolution.^[45] At smaller scales (\ell \gtrsim 100), the CMB power spectrum exhibits a series of acoustic peaks resulting from baryon-photon acoustic oscillations before recombination. In the tightly coupled baryon-photon fluid, primordial overdensities drive compressions and rarefactions, with odd-numbered peaks (first, third, etc.) arising from enhanced compression phases in potential wells, while even peaks stem from rarefaction phases. The position of the first acoustic peak at \ell \approx 220 corresponds to the angular scale of the sound horizon at recombination, \theta_s \approx 0.6^\circ, which measures the distance sound waves traveled in the plasma, roughly 150 Mpc comoving. The relative heights of these peaks are sensitive to the baryon density \Omega_b h^2, with higher baryon loading suppressing even peaks relative to odd ones due to increased inertia during rarefactions. These oscillations provide a standard ruler for the universe's expansion history and confirm the tight coupling of baryons and photons until recombination. CMB polarization anisotropies, at about 10% of the temperature signal amplitude, offer complementary insights into the same primordial perturbations through Thomson scattering at the recombination surface. The dominant E-mode polarization, characterized by curl-free patterns aligned with scalar perturbations, arises from the quadrupole moment of the photon distribution scattering off free electrons, producing a power spectrum C_\ell^{EE} that peaks at \ell \approx 200 near the first acoustic peak scale. E-modes trace the velocity gradients in the baryon-photon fluid, enhancing our understanding of acoustic damping and diffusion. In contrast, B-mode polarization, with curl patterns, is primarily sourced by primordial tensor (gravitational wave) modes from inflation or secondary effects like gravitational lensing of E-modes by large-scale structure; detecting primordial B-modes would directly probe the energy scale of inflation via the tensor-to-scalar ratio r. Current limits from ground-based experiments like BICEP/Keck, combined with CMB data, place r < 0.036 at 95% confidence.^[46] Late-time evolution introduces the integrated Sachs-Wolfe (ISW) effect, a secondary anisotropy where photons experience a net energy shift as they traverse time-varying gravitational potentials during the transition to dark energy domination. The temperature perturbation is given by \Delta T / T \approx -2 \int \dot{\Phi} \, dl / c, where \dot{\Phi} is the time derivative of the potential along the photon path, reflecting the decay of potentials as matter density dilutes relative to dark energy. This effect contributes on large scales (\ell \lesssim [20](/page/2point0)), boosting the low-\ell tail of the CMB power spectrum and providing evidence for cosmic acceleration. Cross-correlations between the ISW signal and large-scale structure tracers, such as galaxy surveys, confirm this decay and constrain dark energy properties.^[45] High-precision measurements from the Planck satellite have solidified these features as evidence for inflationary initial conditions. The 2018 Planck analysis yields a scalar spectral index n_s = 0.9649 \pm 0.0043 (68% CL), deviating from scale invariance (n_s = 1) at over 8\sigma significance and favoring a tilted spectrum consistent with slow-roll inflation. Combined with BICEP/Keck data, it tightens the upper limit to r < 0.036 (95% CL), ruling out many large-field inflation models while allowing small-field ones. Planck also extracts the CMB lensing potential using quadratic estimators, which reconstruct the deflection field from higher-order correlations in the observed anisotropies, confirming the lensing power spectrum and validating the integrated effects of structure growth on the CMB. These results underscore the CMB as a cornerstone for structure formation, linking primordial quantum fluctuations to the observed cosmic web.^[17]^[47]

Theoretical and Numerical Modeling

Analytical Perturbation Theory

Analytical perturbation theory extends the linear regime of cosmological density perturbations into mildly nonlinear scales by expanding the density field and velocity potential in powers of the initial fluctuations. This approach, known as standard perturbation theory (SPT), provides semi-analytic predictions for statistics like the power spectrum and higher-order correlations, capturing mode coupling effects that linear theory neglects. In SPT, the density contrast \delta(\mathbf{k}, \eta) is expanded as \delta = \delta^{(1)} + \delta^{(2)} + \cdots, where \delta^{(1)}(\mathbf{k}, \eta) = D(\eta) \delta_0(\mathbf{k}) is the linear growing mode with growth factor D(\eta), and higher orders involve convolutions over lower-order terms. The second-order contribution is given by

\delta^{(2)}(\mathbf{k}) = \int \frac{d^3\mathbf{k_1} d^3\mathbf{k_2}}{(2\pi)^3} \, \delta^{(1)}(\mathbf{k_1}) \delta^{(1)}(\mathbf{k_2}) \, F_2(\mathbf{k_1}, \mathbf{k_2}),

where the symmetric kernel F_2(\mathbf{k_1}, \mathbf{k_2}) encodes nonlinear interactions and is approximately F_2(\mathbf{k_1}, \mathbf{k_2}) = \frac{5}{7} + \frac{1}{2} \frac{\mathbf{k_1} \cdot \mathbf{k_2}}{k_1 k_2} (\frac{k_1}{k_2} + \frac{k_2}{k_1}) + \frac{2}{7} \left( \frac{\mathbf{k_1} \cdot \mathbf{k_2}}{k_1 k_2} \right)^2 in the Einstein-de Sitter limit. This kernel arises from solving the perturbed continuity, Euler, and Poisson equations to second order, capturing the generation of power on small scales from large-scale modes. The resulting one-loop power spectrum P(k) = P_L(k) + P_{13}(k) + P_{22}(k) includes terms like P_{22}(k) \propto \int d^3q \, [F_2(\mathbf{k}, -\mathbf{q})]^2 P_L(q) P_L(|\mathbf{k - q}|), which quantifies the transfer of power across wavenumbers. These expressions are accurate up to k \sim 0.1 \, h/\mathrm{Mpc} at low redshifts, as validated against simulations. Redshift-space distortions (RSD) arise because galaxy redshifts mix true distances with peculiar velocities, distorting the observed clustering along the line of sight. In linear theory, the redshift-space power spectrum is P^s(k, \mu) = P(k) (1 + \beta \mu^2)^2, where \mu = \hat{\mathbf{k}} \cdot \hat{\mathbf{n}} is the angle to the line of sight, \beta = f/b with growth rate f \approx \Omega_m^{0.55} and bias b. This Kaiser formula predicts enhanced power along the line of sight (\mu \approx 1) due to coherent infall, forming "pancakes," while small-scale random motions cause elongation into "fingers-of-God." Higher-order PT extends this to include terms like \mu^4 and \mu^6, improving fits to data on scales k \lesssim 0.2 \, h/\mathrm{Mpc}. Observations from surveys like 2dFGRS confirm \beta \approx 0.4-0.5 at z=0, consistent with \LambdaCDM. The bispectrum B(\mathbf{k_1}, \mathbf{k_2}, \mathbf{k_3}), the Fourier transform of the three-point correlation function, measures deviations from Gaussianity in the density field. At tree-level in PT, it is B(\mathbf{k_1}, \mathbf{k_2}, \mathbf{k_3}) = 2 F_2(\mathbf{k_1}, \mathbf{k_2}) P(k_1) P(k_2) + \mathrm{cycl.}, where the cyclic permutations account for all mode couplings. This gravitational contribution generates hierarchical non-Gaussianity, with amplitude scaling as the square of the linear power. Primordial non-Gaussianity, parameterized by f_\mathrm{NL}, adds a term \propto f_\mathrm{NL} [P(k_1) P(k_2) + \mathrm{cycl.}]; current constraints from CMB and large-scale structure yield f_\mathrm{NL} = -0.9 \pm 5.1 (68% CL), ruling out significant local-type deviations at the percent level. The bispectrum's squeezed limit (k_3 \ll k_1 \approx k_2) is particularly sensitive to both gravitational evolution and primordial signals, aiding parameter inference. Time-dependent or renormalized perturbation theory (RPT) addresses the breakdown of SPT on intermediate scales by resumming infinite series of diagrams, focusing on propagators that describe how initial perturbations evolve. The propagator G(k, \eta_1, \eta_2) = \langle \delta(\mathbf{k}, \eta_2) \delta_0(-\mathbf{k}, \eta_1) \rangle / P_L(k) decays exponentially for large k, as G(k, \eta) \propto \exp\left( -k^2 \sigma_d^2(\eta)/2 \right), where \sigma_d^2 is the displacement variance from the Zeldovich approximation. This resummation damps small-scale power due to large-scale displacements, improving predictions for the power spectrum by up to 20% at k \sim 0.5 \, h/\mathrm{Mpc} compared to SPT, and matches N-body results over a wider range of scales and redshifts. RPT thus bridges linear and nonlinear regimes analytically. The halo model provides a complementary framework by decomposing the power spectrum into contributions from within halos and between them: P(k) = P^{1h}(k) + P^{2h}(k). The one-halo term P^{1h}(k) = \int dm \, n(m) |u(k|m)|^2, with halo mass function n(m) and profile Fourier transform u(k|m), dominates on small scales (k \gtrsim 1 \, h/\mathrm{Mpc}) and captures intra-halo clustering via NFW profiles. The two-halo term P^{2h}(k) = [b_1 P_L(k)]^2 + \cdots uses linear bias b_1(m) for large scales, with higher-order biases from peak-background split. This ansatz bridges nonlinear halo profiles with linear large-scale clustering, accurately reproducing the power spectrum across scales when calibrated to simulations, and extends to higher-point functions via similar decompositions.^[48]

N-Body and Hydrodynamic Simulations

N-body simulations provide a computational framework for modeling the nonlinear gravitational clustering of collisionless dark matter particles, which dominate structure formation in the universe. These methods discretize dark matter into a large number of particles and evolve their trajectories under self-gravity by solving the Poisson equation for the gravitational potential, typically using particle-mesh algorithms for long-range forces or tree-based methods like the Barnes-Hut approximation for hierarchical force computations. The GADGET code exemplifies this approach, employing a parallel tree algorithm to efficiently calculate interactions in simulations containing up to $10^{10} particles, allowing the tracking of structure growth from initial conditions derived from linear perturbation theory over cosmic timescales.^[49]^[50] To incorporate baryonic gas dynamics alongside dark matter, hydrodynamic extensions couple N-body gravity solvers with fluid descriptions of gas evolution. Smoothed particle hydrodynamics (SPH), as implemented in GADGET, represents gas as particles with smoothed kernel estimates for density and pressure, enabling the simulation of shocks, cooling, and heating processes. Alternatively, the moving-mesh code AREPO uses a dynamic Voronoi tessellation that adapts to the flow, combining Lagrangian flexibility with Eulerian accuracy to better resolve gas instabilities, turbulence, and angular momentum transport in galaxy formation scenarios; it includes modules for radiative cooling, star formation based on local density thresholds, and feedback from supernovae and active galactic nuclei.^[51]^[52] The IllustrisTNG simulation suite demonstrates the power of these combined methods, running magnetohydrodynamical simulations in periodic boxes of 100 Mpc side length from redshift z \approx 20 to z = 0 using AREPO, with dark matter particle masses around $5 \times 10^7 M_\odot and gas cell masses of order $10^6 M_\odot. These runs reproduce key aspects of galaxy evolution, including realistic morphologies from disk settling and bulge formation, supermassive black hole growth via accretion and mergers, and the amplification of magnetic fields through dynamo effects in turbulent interstellar media.^[53]^[54] Numerical resolution imposes fundamental limits on these simulations, with gravitational softening lengths \epsilon \approx 1 kpc preventing artificial collapse in dense regions while preserving two-body relaxation times longer than the Hubble time. This achieves a dynamic mass resolution range of approximately $10^6, from individual stars to galaxy clusters, facilitating subgrid models for unresolved processes like supernova feedback, where each event injects thermal energy E_\mathrm{SN} = 10^{51} erg to drive galactic winds and regulate star formation rates.^[53] Recent advancements leverage exascale computing for unprecedented scale and fidelity, as seen in the AbacusSummit suite of over 150 N-body simulations with up to $3.3 \times 10^{11} particles across varied cosmologies, enabling precise modeling of massive neutrino contributions to clustering suppression and tests of modified gravity theories for Bayesian parameter inference from upcoming surveys. By 2025, these simulations continue to integrate such extensions, supporting high-precision predictions for dark energy constraints while pushing resolution limits through adaptive refinement techniques.^[55]^[56]

References

[1]
[1505.02821] Cosmological Structure Formation - arXiv
May 11, 2015 · Most star-forming galaxies are very clumpy; that galaxies often undergo compaction which reduces their radius and increases their central density.
[2]
Cosmological parameters derived from the final (PR4) Planck data ...
We present constraints on cosmological parameters using maps from the last Planck data release (PR4). In particular, we detail an upgraded version of the ...
[3]
Simulations of Structure Formation in the Universe - E. Bertschinger
This article reviews the techniques and results of cosmological structure formation simulations since the early 1980s.
[4]
[PDF] 23. Inflation - Particle Data Group
In the following we will review scalar field cosmology in general relativity and the spectra of primordial fluctuations produced during inflation, before ...
[5]
None
Nothing is retrieved...<|separator|>
[6]
None
### Summary of Key Concepts from arXiv:0904.4584
[7]
Quantum field theory in de Sitter space: renormalization by point ...
We examine the modes of a scalar field in de Sitter space and construct quantum two-point functions. These are then used to compute a finite stress tensor.Missing: URL | Show results with:URL
[8]
None
Below is a merged summary of the Planck 2018 data for the primordial power spectrum spectral index $ n_s $, combining all the provided segments into a single, comprehensive response. To retain all details efficiently, I will use a table in CSV format for the key values and details, followed by a narrative summary and a consolidated list of useful URLs. This ensures all information is preserved while maintaining clarity and density.
[9]
Reheating in Inflationary Cosmology: Theory and Applications - arXiv
Jan 15, 2010 · Reheating is an important part of inflationary cosmology. It describes the production of Standard Matter particles after the phase of accelerated expansion.Missing: seminal | Show results with:seminal
[10]
[astro-ph/0402129] Cosmological perturbation theory - arXiv
Feb 5, 2004 · This is a review on cosmological perturbation theory. After an introduction, it presents the problem of gauge transformation.
[11]
The growth of density perturbations in Friedmann model universes ...
We examine the linear growth of density perturbations in homogeneous isotropic (Friedmann) model universes, including the effect of a decoupled radiation ...
[12]
[astro-ph/0507263] Cosmic Growth History and Expansion ... - arXiv
Jul 11, 2005 · The cosmic expansion history tests the dynamics of the global evolution of the universe and its energy density contents.
[13]
https://link.springer.com/article/10.1007/BF00639395
[14]
[1807.06209] Planck 2018 results. VI. Cosmological parameters - arXiv
Jul 17, 2018 · Abstract:We present cosmological parameter results from the final full-mission Planck measurements of the CMB anisotropies.
[15]
How Exactly Did the Universe Become Neutral? - IOPscience
We present a refined treatment of H, He I, and He II recombination in the early universe. The difference from previous calculations is that we use multilevel ...
[16]
Baryonic Features in the Matter Transfer Function - IOPscience
Baryonic Features in the Matter Transfer Function. Daniel J. Eisenstein and Wayne Hu. © 1998. The American Astronomical Society. All rights reserved.
[17]
[PDF] arXiv:1410.2509v1 [astro-ph.CO] 9 Oct 2014
Oct 9, 2014 · the radiation density drops faster than the dark-matter density. At a redshift zeq ∼ 3400 the equality of radia- tion and matter densities ...
[18]
[PDF] Model-Independent Measurement of the Matter-Radiation Equality ...
The equality redshift zeq is given by the point where the blue (total relativistic energy density) and green (total non- relativistic energy density) lines ...
[19]
[PDF] Cosmological parameters after WMAP5: forecasts for Planck ... - arXiv
Nov 26, 2008 · In most cosmologies, the size of the sound horizon at last scatter- ing provides a standard ruler at a redshift of zLS ≃ 1090. Indeed, most ...
[20]
[PDF] arXiv:2304.06072v1 [hep-ph] 12 Apr 2023
Apr 12, 2023 · In ΛCDM, the IGM is almost completely neutral after recombination and the photons free-stream to good approximation; ionization of the IGM ...
[21]
[PDF] Cosmological Perturbations - arXiv
Thus, after recombination large-scale perturbations in the non-relativistic component grow as δb ∝ a ∝ t2/3. 3.4.2 The radiation dominated era. Before ...
[22]
[PDF] arXiv:1504.03587v1 [astro-ph.CO] 14 Apr 2015
Apr 14, 2015 · After recombination, baryons and photons evolved independently. The main point of this part of lectures is that by analysing the density ...
[23]
[PDF] 27. Dark Matter - Particle Data Group
Dec 1, 2023 · The presence of a dominant non-relativistic (“cold”) pressure-less matter component decoupled from the thermal bath well before ...
[24]
Axion dark matter: What is it and why now? | Science Advances
Feb 23, 2022 · The axion has emerged in recent years as a leading particle candidate to provide the mysterious dark matter in the cosmos.
[25]
An approximate theory for large density perturbations. - ADS
Zel'dovich, Ya. B. Abstract. An approximate solution is given for the problem of the growth of perturbations during the expansion of matter without pressure.Missing: URL | Show results with:URL
[26]
Clustering in real space and in redshift space - Oxford Academic
Nick Kaiser, Clustering in real space and in redshift space, Monthly Notices of the Royal Astronomical Society, Volume 227, Issue 1, July 1987, Pages 1–21 ...
[27]
https://www.science.org/doi/10.1126/sciadv.abj3618
[28]
http://ui.adsabs.harvard.edu/abs/1970A&A.....5...84Z/abstract
[29]
https://academic.oup.com/mnras/article/227/1/1/1065830
[30]
https://ui.adsabs.harvard.edu/abs/1974ApJ...187..425P/abstract
[31]
II. Clustered star formation and the influence of metal line cooling
The introduction of metals by the first supernovae is generally thought to have driven the transition, given the enhanced cooling and thus fragmentation in ...Star Formation In The First... · Results · Star Formation Efficiency
[32]
Mechanical feedback from active galactic nuclei in galaxies, groups ...
May 23, 2012 · We show that cooling x-ray atmospheres and the ensuing star formation and nuclear activity are probably coupled to a self-regulated feedback ...
[33]
The observational signature of the first H ii regions - Oxford Academic
We here carry out radiation hydrodynamics simulations of the evolution of H ii regions around massive Pop III stars in minihaloes, giving us a detailed ...
[34]
[1607.03155] The clustering of galaxies in the completed SDSS-III ...
Jul 11, 2016 · This paper presents cosmological results from galaxy clustering data, measuring DM and H, and using BAO, AP, and RSD methods, and finding ...
[35]
BOSS: Dark Energy and the Geometry of Space - SDSS
The SDSS's Baryon Oscillation Spectroscopic Survey (BOSS) mapped the spatial distribution of luminous red galaxies (LRGs) and quasars to detect the ...
[36]
[2007.08991] The Completed SDSS-IV extended Baryon Oscillation ...
Jul 17, 2020 · The survey uses data from SDSS, SDSS-II, BOSS, and eBOSS to measure baryon acoustic oscillation (BAO) and growth rate parameter (f\sigma_8) for ...
[37]
[2105.13549] Dark Energy Survey Year 3 Results: Cosmological ...
Abstract:We present the first cosmology results from large-scale structure in the Dark Energy Survey (DES) spanning 5000 deg^2.
[38]
[2503.19196] Euclid Quick Data Release (Q1). First detections from ...
Mar 24, 2025 · This paper showcases, for the first time, the performance of cluster finders on Euclid data and presents examples of validated clusters in the Quick Release 1 ...Missing: 2024 | Show results with:2024
[39]
https://www.sdss.org/surveys/boss/
[40]
https://arxiv.org/abs/2007.08991
[41]
[astro-ph/9706147] A CMB Polarization Primer - arXiv
We present a pedagogical and phenomenological introduction to the study of cosmic microwave background (CMB) polarization to build intuition about the ...Missing: EB | Show results with:EB
[42]
[PDF] 25. Cosmological Parameters - Particle Data Group
Dec 1, 2023 · ... ns) gives a 95% confidence upper limit of r < 0.036 at the scale 0.05 Mpc. −1. [63]. The density perturbation spectral index is clearly ...Missing: n_s | Show results with:n_s
[43]
[astro-ph/0206508] Halo Models of Large Scale Structure - arXiv
Jun 28, 2002 · Halo models associate mass with dark matter halos, using models of halo distribution to estimate large-scale density and velocity fields.
[44]
Simulating cosmic structure formation with the gadget-4 code
Here, we discuss recent methodological progress in the gadget code, which has been widely applied in cosmic structure formation over the past two decades.
[45]
GADGET-4 - MPA Garching
This site provides basic information about the GADGET-4 simulation code, which is a parallel cosmological N-body and SPH code meant for simulations of cosmic ...
[46]
Moving mesh cosmology: the hydrodynamics of galaxy formation
Sep 15, 2011 · We present a detailed comparison between the well-known SPH code GADGET and the new moving-mesh code AREPO on a number of hydrodynamical test problems.
[47]
AREPO - Code - HITS - Heidelberg Institute for Theoretical Studies
Oct 28, 2014 · The AREPO code is a cosmological hydrodynamical simulation code on a fully dynamic unstructured mesh. At present, hydrodynamic cosmological ...
[48]
First results from the IllustrisTNG simulations: the stellar mass ...
IllustrisTNG (Fig. 1) is an ongoing suite of MHD cosmological simulations that model the formation and evolution of galaxies within the ΛCDM paradigm. It builds ...
[49]
Data Access - Specifications - IllustrisTNG
This page contains a comprehensive description of the simulation snapshots, group catalogs, merger trees, and supplementary data sets.
[50]
AbacusSummit: a massive set of high-accuracy, high-resolution N ...
This simulation suite, named AbacusSummit, consists of 150 simulation boxes, spanning 97 cosmological models, most of which have 330 billion particles with ...
[51]
[PDF] JCAP09(2025)007 - OSTI.GOV
Using the AbacusSummit simulations described in section 3, these models are tuned to approximately reproduce the clustering of the DESI One-Percent Survey [8] ...