Fact-checked by Grok 2 weeks ago

Galaxy formation and evolution

Galaxy formation and evolution encompasses the physical processes by which galaxies—immense assemblies of stars, gas, dust, and dark matter—emerged from the primordial fluctuations of the early universe and transformed over cosmic history. Following the Big Bang approximately 13.8 billion years ago, these structures began assembling within the first 500 million to 1 billion years, initially as small, irregular clumps of gas and stars that coalesced around denser pockets of matter influenced by cosmic inflation and gravitational instability.^[1]^[2]^[3] Driven by the hierarchical merging of dark matter halos and the accretion of intergalactic gas, galaxies evolved into diverse morphologies, including spirals, ellipticals, and irregulars, while undergoing cycles of star formation, feedback from stellar processes and supermassive black holes, and interactions that reshaped their structures and compositions.^[4]^[2] The initial formation phase was dominated by rapid gas cooling and collapse within dark matter-dominated halos, where ordinary baryonic matter followed the gravitational scaffolding provided by dark matter—estimated to be about five times more abundant than visible matter.^[4] This led to the birth of the first stars and protogalaxies around 300 million years after the Big Bang, as observed in deep-field images from the James Webb Space Telescope revealing compact, clumpy systems with high star formation rates.^[5]^[2] Over subsequent billions of years, mergers between galaxies played a crucial role in growth; for instance, many massive galaxies (with studies estimating ~40-70% having experienced at least one major merger) since the universe was about 6–7 billion years old, often triggering bursts of star formation by compressing interstellar gas.^[6]^[4]^[2] Star formation rates across the galaxy population peaked around 10–11 billion years ago, when the universe was roughly a third of its current age, with galaxies producing stars up to 10 times faster than modern rates due to abundant cold gas supplies from cosmic filaments and streams.^[3] This era, known as cosmic noon, saw intense activity in dusty, starburst galaxies, where feedback from supernovae, stellar winds, and accreting supermassive black holes—often millions to billions of solar masses at galactic centers—began regulating gas reservoirs and eventually quenching star formation in many systems.^[3]^[2] By contrast, present-day galaxies like the Milky Way form only a few stars per year, reflecting a steadier, more evolved state shaped by these regulatory mechanisms.^[3] Modern understanding of these processes relies on multi-wavelength observations from telescopes such as the Hubble Space Telescope, which has imaged over 10,000 galaxies in fields like the Hubble Ultra Deep Field dating back 13 billion years, and the James Webb Space Telescope (JWST), which probes the earliest epochs to clarify the interplay between galaxy assembly and black hole growth.^[2] These data, combined with cosmological simulations, highlight ongoing evolution: for example, the Milky Way is projected to collide with the Andromeda Galaxy in about 4.5 billion years, potentially forming a new elliptical galaxy.^[4]^[2] Such insights underscore the dynamic, interconnected nature of galaxy evolution within the expanding universe.^[4]

Cosmological Context

Lambda-CDM paradigm

The Lambda-CDM paradigm, often referred to as the concordance model of cosmology, posits that the universe consists of three primary components: cold dark matter (CDM), which dominates the non-baryonic matter content; a cosmological constant Λ, interpreted as dark energy driving the accelerated expansion; and ordinary baryonic matter, making up stars, gas, and other visible structures.^[7] This model successfully accounts for a wide range of cosmological observations, including the large-scale structure of the universe and the cosmic microwave background (CMB), though recent data as of 2025 reveal tensions, such as the Hubble constant (H_0) discrepancy between CMB and local measurements, and hints from the Dark Energy Spectroscopic Instrument (DESI) of possible time-varying dark energy.^[7]^[8] Within this framework, the present-day energy densities are characterized by key parameters derived from Planck satellite data: the total matter density parameter Ω_m ≈ 0.315, the dark energy density parameter Ω_Λ ≈ 0.685, and the Hubble constant H_0 ≈ 67.4 km/s/Mpc.^[7] Galaxy formation in the Lambda-CDM model is driven by the gravitational collapse of overdense regions into dark matter halos, which serve as the gravitational wells where baryonic matter cools and condenses to form stars and galaxies. These halos emerge from tiny initial density fluctuations, amplified by gravity over cosmic time, with dark matter providing the bulk of the mass and thus dictating the collapse dynamics. The process unfolds hierarchically through gravitational instability in an expanding universe, where small-scale structures form first and subsequently merge to build larger systems, such as galaxy clusters, reflecting the bottom-up assembly predicted by the cold dark matter power spectrum. In the matter-dominated era of cosmic history, the linear growth of density perturbations δ is described by the relation

\delta \propto a,

where a is the scale factor of the universe, indicating that overdensities grow proportionally with the expansion until nonlinear effects take over, leading to halo virialization. These primordial fluctuations trace back to quantum variations stretched by cosmic inflation in the early universe.

Initial conditions from CMB

The cosmic microwave background (CMB) serves as the relic radiation from the epoch of recombination, when the universe cooled sufficiently at redshift z \approx 1100 for electrons and protons to form neutral hydrogen, decoupling photons from matter and allowing the universe to become transparent. This snapshot of the early universe, approximately 380,000 years after the Big Bang, encodes the initial conditions for structure formation through tiny temperature fluctuations of order $10^{-5} K. Observations of these anisotropies have been pivotal, beginning with the Cosmic Background Explorer (COBE) satellite's detection in 1992, followed by the Wilkinson Microwave Anisotropy Probe (WMAP) from 2001 to 2010, and culminating in the high-precision measurements from the Planck satellite between 2009 and 2013.^[9]^[7] The angular power spectrum of CMB temperature anisotropies reveals a series of peaks arising from baryon acoustic oscillations (BAO) in the primordial plasma before recombination, where sound waves in the photon-baryon fluid left imprints on the distribution of matter and radiation. The first peak, located at a multipole moment \ell \approx 220, corresponds to the sound horizon at recombination and provides strong evidence for a spatially flat universe, consistent with the \LambdaCDM paradigm's curvature parameter \Omega_k \approx 0. Subsequent peaks reflect higher-order modes damped by diffusion processes, with the overall spectrum measured to exquisite precision by Planck, confirming the acoustic peak structure predicted by linear perturbation theory.^[7]^[10] These anisotropies trace scalar perturbations in the early universe, modeled as a nearly Gaussian random field with statistical properties inherited from cosmic inflation. The primordial power spectrum of these curvature perturbations is characterized by P(k) \propto k^{n_s - 4}, where k is the comoving wavenumber and the scalar spectral index n_s \approx 0.965 \pm 0.004 indicates near scale-invariance, deviating slightly from the Harrison-Zel'dovich limit of n_s = 1. The evolution of these perturbations from the inflationary epoch to recombination is described by the transfer function T(k), which modulates the power on different scales; small-scale modes (k \gtrsim 0.01 Mpc^{-1}) are suppressed due to Silk damping, a diffusive process where random photon scattering smears out fluctuations during the tight-coupling regime before recombination.^[7]^[10] The amplitude of matter fluctuations on scales of $8 h^{-1} Mpc, quantified by \sigma_8 \approx 0.811 \pm 0.006, is directly inferred from the CMB power spectrum normalization and sets the overall strength of the initial density seeds for gravitational collapse leading to galaxy formation.^[7] This parameter, combined with the power spectrum shape, provides the foundational input for \LambdaCDM simulations of cosmic structure growth.^[10]

Observational Foundations

Galaxy morphologies and properties

Galaxies exhibit a diverse range of morphologies, primarily classified using the Hubble sequence, which organizes them into a tuning-fork diagram based on visual appearance and structural features. This scheme divides galaxies into ellipticals (E), lenticulars (S0), spirals (S), barred spirals (SB), and irregulars (Irr). Elliptical galaxies, denoted E0 to E7 according to increasing ellipticity, feature smooth, featureless envelopes of old stars with no significant disk or spiral arms. Lenticular galaxies (S0) possess a prominent bulge and a thin disk but lack spiral structure, appearing as intermediate between ellipticals and spirals. Spiral galaxies range from Sa (tightly wound arms, large bulge) to Sd (loosely wound arms, small bulge), often with a central bar in SBa to SBd subtypes, characterized by prominent disks with ongoing star formation in arms. Irregular galaxies, including dwarf irregulars, show chaotic structures without clear symmetry, often rich in gas and young stars.^[11] Key observed properties of galaxies include their luminosity distribution and structural scaling. The galaxy luminosity function, which describes the number density of galaxies per unit luminosity interval, is well-fit by the Schechter function \phi(L) \propto L^{-\alpha} \exp(-L/L_*), where L_* is the characteristic luminosity and \alpha \approx 1.2-1.5 governs the faint-end slope, indicating a steeper decline toward lower luminosities.^[12] In terms of size and mass, the effective radius r_e of disk-dominated galaxies scales approximately as r_e \propto M^{0.5}, where M is the stellar mass, reflecting a roughly constant surface mass density in these systems.^[13] Several empirical scaling relations link luminosity to dynamical properties, providing insights into galaxy structure. For spiral galaxies, the Tully-Fisher relation correlates absolute luminosity L with the maximum rotation velocity v as L \propto v^4, enabling distance estimates and highlighting the role of rotational support in luminous disks. In contrast, elliptical galaxies follow the Faber-Jackson relation, where L \propto \sigma^4 and \sigma is the central stellar velocity dispersion, underscoring velocity dispersion as the primary support mechanism against gravity in these pressure-dominated systems. Galaxies also segregate in the color-magnitude diagram, a plot of rest-frame color (e.g., u-r) versus absolute magnitude, revealing bimodality in their stellar populations. The red sequence comprises quiescent, metal-rich, early-type galaxies with older stars and minimal ongoing star formation, occupying fainter magnitudes at redder colors due to lower luminosities from passive evolution. The blue cloud, conversely, hosts star-forming, late-type galaxies with younger, bluer stellar populations, extending to brighter magnitudes and indicating active dust-obscured star formation. Between these lies the green valley, a transitional region for galaxies quenching their star formation. A tight correlation exists between supermassive black hole masses and host galaxy properties, particularly the M_{\rm BH}-\sigma relation, where black hole mass M_{\rm BH} \propto \sigma^4 for the bulge velocity dispersion \sigma, suggesting co-evolution between black holes and their spheroidal hosts. At high redshifts, these morphologies evolve, with early galaxies showing more disturbed and clumpy structures compared to the smoother forms dominant locally.

High-redshift observations and JWST discoveries

High-redshift observations have provided critical insights into the reionization era, spanning redshifts z ≈ 6–10, when the intergalactic medium (IGM) transitioned from neutral to ionized due to the first luminous sources. The Lyman-alpha forest, consisting of absorption lines in quasar spectra from neutral hydrogen clouds, reveals a highly ionized IGM at z < 6, with the density of these absorbers decreasing toward lower redshifts, indicating progressive reionization.^[14] At z > 6, the Gunn-Peterson trough—a broad absorption feature due to resonant scattering of Lyman-alpha photons by neutral hydrogen—emerges in quasar spectra, signaling the presence of a significant neutral hydrogen fraction (x_HI > 10^{-3}) and the onset of reionization, as observed in multiple z > 6 quasars. These features collectively demonstrate that reionization was largely complete by z ≈ 6, driven by ultraviolet photons from early galaxies and quasars. The James Webb Space Telescope (JWST), operational since 2022, has revolutionized high-redshift galaxy studies through deep-field surveys like CEERS (Cosmic Evolution Early Release Science) and GLASS (Grism Lens-Amplified Survey from Space), uncovering a population of massive galaxies at z > 10, mere hundreds of millions of years after the Big Bang.^[15] For instance, in the CEERS field, JWST/NIRCam imaging has identified luminous galaxy candidates at z ≈ 10–12 with photometric redshifts confirmed spectroscopically, displaying rest-frame ultraviolet luminosities exceeding those expected in standard models.^[15] A prominent example is GN-z11, spectroscopically confirmed at z = 10.6 via JWST/NIRSpec observations in the JADES (JWST Advanced Deep Extragalactic Survey) program, revealing a stellar mass of approximately 10^{9.1} M_⊙ and active star formation with a specific star formation rate of ~10^{-7.7} yr^{-1}.^[16] These discoveries, spanning 2022–2025, indicate that massive galaxies (M_* > 10^9 M_⊙) were already assembled by z > 10, challenging predictions of gradual buildup. JWST has also revealed structured morphologies in these early galaxies, including clumpy, rotating disks at z ≈ 9–10. In 2025, observations identified a primordial rotating disk galaxy, approximately 900 million years post-Big Bang (z ≈ 10), composed of at least 15 star-forming clumps, with the clumpiness exceeding expectations from simulations and later-epoch galaxies.^[17] These features suggest rapid disk formation and dynamical maturity shortly after cosmic dawn, potentially driven by efficient gas accretion.^[17] These observations introduce tensions with hierarchical assembly models, as the prevalence of massive galaxies at z > 10 implies faster stellar mass growth than predicted, requiring either enhanced star formation efficiency or alternative formation pathways beyond standard cold dark matter scenarios.^[18] Additionally, galaxy size evolution, traced by effective radii r_e, follows a relation r_e ∝ (1 + z)^{-0.71 ± 0.19} at fixed stellar mass and rest-frame wavelength, which is flatter than the steeper (1 + z)^{-1} to -2 predicted for compact early systems, indicating less pronounced size contraction over cosmic time.^[19] JWST data further support the early emergence of galaxy clusters at z > 6, with 2024 studies identifying overdensities of massive galaxies in fields like Abell S1063, suggesting protoclusters assembling within the first billion years and accelerating large-scale structure formation.^[18]

Theoretical Models

Top-down formation

The top-down formation model, also known as the monolithic collapse scenario, posits that galaxies originate from the rapid gravitational collapse of a single, large-scale protogalactic gas cloud. This framework was first proposed by Eggen, Lynden-Bell, and Sandage in 1962, based on kinematic evidence from the motions of old stars in the Milky Way, suggesting an early collapse that funneled gas inward to form the galaxy's structure.^[20] In this model, the protogalactic cloud, initially part of the expanding universe post-Big Bang, decouples from the Hubble flow and undergoes contraction under its own gravity, leading to the simultaneous formation of the halo, disk, and bulge components.^[21] Key assumptions of the model include uniform density perturbations in the primordial gas cloud, which allow for a coherent, homogeneous collapse without significant fragmentation.^[21] Dissipation through radiative cooling of the gas is essential, enabling the material to shed angular momentum and settle into a thin, rotationally supported disk while the halo forms from stars generated during the infall.^[22] The collapse proceeds on the free-fall timescale, given by

t_{\rm ff} = \left( \frac{3\pi}{32 G \rho} \right)^{1/2},

where G is the gravitational constant and \rho is the mean density of the cloud; for a typical protogalactic density of \rho \approx 10^{-23} g cm^{-3}, this yields t_{\rm ff} \approx 10^8 yr, allowing rapid structure formation in the early universe.^[22] Angular momentum conservation during the collapse further shapes the disk, with initial cloud rotation preventing total central infall.^[20] The model predicts predominantly old stellar populations across galaxy components, as star formation occurs swiftly during the collapse phase, resulting in minimal subsequent evolution or mergers for isolated systems.^[21] It is particularly suited to explaining massive ellipticals or early-type disks in isolation, where a single dissipative event builds the bulk of the stellar mass.^[23] Despite its historical influence, the monolithic collapse model faces significant criticisms for incompatibility with the cold dark matter (CDM) paradigm, which favors hierarchical structure formation through the merging of smaller subunits rather than a single large cloud.^[24] The assumption of uniform density overlooks the small-scale power in CDM initial conditions, which promotes early fragmentation into substructures that the top-down scenario lacks mechanisms to incorporate.^[21]

Bottom-up hierarchical assembly

The bottom-up hierarchical assembly model describes galaxy formation as a progressive merging process where small dark matter halos and their embedded dwarf galaxies coalesce to build larger systems over cosmic history. In this paradigm, the earliest structures collapse from tiny density perturbations at high redshifts, with subsequent mergers driving the growth of massive galaxies at lower redshifts. This contrasts with monolithic collapse scenarios by emphasizing multi-scale accretion and merging as the dominant mechanism, naturally arising from gravitational instability in the cold dark matter (CDM) component of the Lambda-CDM cosmology. A foundational tool for quantifying halo abundances in this model is the Press-Schechter formalism, which statistically predicts the number density of collapsed objects as a function of mass. The mass function is given by

\frac{dn}{dM} = \sqrt{\frac{2}{\pi}} \frac{\bar{\rho}}{M^2} \frac{\delta_c}{\sigma(M)} \left| \frac{d \ln \sigma}{d \ln M} \right| \exp\left( -\frac{\delta_c^2}{2 \sigma^2(M)} \right),

where \bar{\rho} is the mean density, \delta_c \approx 1.686 is the linearly extrapolated critical overdensity threshold for spherical collapse in an Einstein-de Sitter universe, and \sigma(M) represents the root-mean-square fluctuation amplitude smoothed over the mass scale M. This expression, derived assuming Gaussian initial conditions, yields a characteristic mass scale below which small halos dominate at early times, enabling the hierarchical buildup.^[25] Hierarchical merging proceeds as dwarf galaxies in low-mass halos accrete onto larger hosts, with the merger rate in CDM cosmologies scaling roughly as (1+z)^{2.5}, reflecting faster coalescence at higher redshifts due to denser environments. The extended Press-Schechter formalism builds on this by providing the conditional probability distribution for progenitor masses in merger trees, allowing reconstruction of a halo's assembly history through random walks in the density field. Specifically, the probability that a halo of mass M at redshift z had progenitors of masses M_1, M_2, \dots at higher z' follows a multivariate Gaussian form tied to changes in \sigma(M), facilitating Monte Carlo simulations of merger sequences. Key predictions of this framework include the abundance of satellite galaxies, which emerge as surviving subhalos from minor mergers and accretion, matching observed dwarf populations around Milky Way-like systems. Additionally, major mergers between comparable-mass progenitors are predicted to disrupt disks and form spheroidal ellipticals, a process central to explaining the prevalence of early-type galaxies at intermediate redshifts. These features align closely with large-scale N-body simulations of CDM structure formation, which reproduce the merger-driven growth and substructure statistics. However, recent James Webb Space Telescope (JWST) observations of ultra-massive galaxies at z \gtrsim 7, with stellar masses exceeding $10^{10.5} M_\odot and high star formation rates, indicate an accelerated assembly phase in the first billion years that appears inconsistent with the standard hierarchical timeline, potentially requiring modifications to early-universe physics or feedback processes.^[26] As of 2025, continued JWST observations reinforce these challenges, with studies suggesting enhanced early star formation or adjustments to dark matter models to reconcile the data.^[27]

Disk Galaxy Formation

Angular momentum conservation

In the context of disk galaxy formation, protogalactic gas clouds acquire angular momentum primarily through tidal torques generated by gravitational interactions with neighboring density perturbations during the early stages of cosmic structure formation. This mechanism, first detailed by Peebles (1969), imparts a net rotation to collapsing gas and dark matter halos, with the total angular momentum L scaling as L \propto \lambda J_c, where \lambda \approx 0.05 is the dimensionless spin parameter characterizing the halo's rotational support relative to gravitational collapse, and J_c is the angular momentum for a circular orbit. Cosmological N-body simulations confirm that dark matter halos typically exhibit \lambda \approx 0.05, independent of mass and formation epoch in the \LambdaCDM paradigm. The specific angular momentum j = L/M of the halo, where M is the total mass, follows j \propto \lambda r_{\rm vir} v_{\rm vir}, with r_{\rm vir} and v_{\rm vir} denoting the virial radius and circular velocity, respectively. During the subsequent collapse of baryonic gas within the halo, this specific angular momentum is largely conserved in the absence of significant external torques, causing the material to flatten into a rotating disk rather than forming a spherical configuration. Fall & Efstathiou (1980) developed a seminal model assuming j_{\rm disk} \approx j_{\rm halo} (accounting for the cosmic baryon fraction), which predicts that the resulting disk surface density profile takes an exponential form \Sigma(r) \propto \exp(-r/R_d), where the scale length R_d is set by the inherited halo spin.^[28] Observational evidence for this conservation comes from the flat rotation curves of spiral galaxies, where orbital velocities remain roughly constant at v \approx 200 km/s out to large radii, consistent with the preservation of specific angular momentum from the protogalactic phase and implying a distribution of mass that supports extended rotation without rapid decline.^[29] However, dynamical processes such as non-axisymmetric structures can challenge pure conservation by transporting angular momentum outward: bars redistribute it on secular timescales, while transient spiral arms facilitate radial flows that alter the inner disk profile.

Secular processes and inside-out growth

Secular processes refer to the gradual, internal dynamical evolution of disk galaxies driven by gravitational instabilities and non-axisymmetric structures, such as bars and spiral arms, which redistribute gas and stars over gigayear timescales without requiring major mergers. These processes are particularly prominent in gas-rich disks, where instabilities lead to enhanced star formation and structural changes that shape the galaxy's morphology and radial profiles. In the context of disk galaxy formation, secular evolution facilitates the transformation of initially clumpy or unstable configurations into more stable, extended structures, with angular momentum acquired from the initial collapse playing a foundational role in setting the disk's scale. A key outcome of secular processes is inside-out growth, where star formation activity migrates outward as fresh gas accretes onto the outer disk, expanding the stellar disk radius over time. Theoretical models of gas infall predict that the disk scale length evolves as R_d \propto t^{0.5}, reflecting the buildup of mass with conserved specific angular momentum, which allows outer regions to form younger stars while inner regions age and quench earlier. This pattern is supported by simulations and observations of the Milky Way, where the low-α stellar disk exhibits a 43% increase in half-mass radius over the past 7 billion years, consistent with prolonged gas accretion fueling peripheral star formation. Recent James Webb Space Telescope (JWST) observations of a galaxy at redshift z ≈ 7.4 reveal a compact core surrounded by an extended star-forming disk, providing direct evidence of inside-out growth in the early universe, where star formation rates rise toward the outskirts.^[30]^[31]^[32] One mechanism driving inside-out growth involves the migration of giant star-forming clumps in high-redshift disks. In the clump migration model, massive, gravitationally bound clumps of gas and stars, with masses around $10^8–$10^9 M_\odot, form due to disk instabilities in gas-rich environments at z > 1 and migrate inward via dynamical torques, shedding angular momentum and depositing stars into the central bulge while leaving the outer disk to accrete fresh gas. Noguchi's 1999 simulations demonstrated that these clumps, arising from violent gravitational instabilities in young, clumpy disks, coalesce to form classical bulges, with survival rates enhanced by stellar feedback that resists disruption. Building on this, Dekel et al. (2009) incorporated cold gas streams feeding high-z galaxies, showing that clumps form compact spheroids centrally while the disk grows outward, with clump inward migration timescales of ~0.5–1 Gyr aligning with observed high-redshift morphologies.^[33] Bar instabilities further contribute to secular evolution by driving radial gas inflows that fuel central starbursts and bulge growth. When the Toomre stability parameter Q falls below 1 in the inner disk—indicating susceptibility to axisymmetric perturbations—a bar forms through gravitational collapse, rotating with a pattern speed that exerts torques on the gas, funneling it inward along dust lanes while transporting angular momentum outward. This process, observed in N-body simulations, leads to pseudobulge formation as gas accumulates centrally, with bar-driven inflows enhancing star formation rates by factors of 2–10 in the nuclear regions.^[34] Spiral arms, as quasi-stationary density waves, also play a role in secular processes by facilitating angular momentum transport and gas redistribution. In the Lin-Shu density wave theory, spiral patterns rotate at a constant pattern speed \Omega_p, slower than the orbital speed of stars in the inner disk but faster in the outer regions, causing material to pile up in the arms where star formation is triggered. This wave-like perturbation exchanges angular momentum between the disk and the halo, driving slow radial mixing of stars and gas over the disk lifetime, with arms acting as conduits for outward angular momentum flux that stabilizes the disk against further fragmentation. Observational evidence for inside-out growth manifests in radial color gradients, where disk outskirts appear bluer due to younger, more metal-poor stellar populations formed from recent gas accretion. Multi-wavelength surveys of nearby spirals show negative u-i and NUV-u gradients, with outer disks ~0.2–0.5 mag bluer than inner regions, ruling out dust reddening and instead attributing the trend to age differences of 2–4 Gyr, consistent with prolonged outer star formation. These gradients, measured across thousands of galaxies, strengthen support for secular inside-out assembly, as inner quenching precedes outer buildup without invoking external triggers.

Elliptical and Spheroidal Galaxies

Merger-induced formation

Major mergers between disk galaxies, typically involving mass ratios of 1:1 to 1:4, are a primary mechanism for the formation of elliptical and spheroidal systems. During such events, the gravitational interaction disrupts the orderly rotation of the progenitor disks, leading to a phase of violent relaxation where stars and dark matter redistribute into a hot, pressure-supported configuration. This process forms a stellar halo with a characteristic velocity dispersion of approximately 200-300 km/s, consistent with observed properties of elliptical galaxies.^[35] The inspiral of the satellite galaxy toward the primary is governed by dynamical friction, which decelerates the orbiting body through gravitational wake effects in the host's density field. The Chandrasekhar dynamical friction timescale is given by \tau_\mathrm{fric} \propto \frac{M_\mathrm{halo}}{M_\mathrm{sat}} \frac{v_c^2}{\sigma^2} \frac{1}{\rho}, where M_\mathrm{halo} and M_\mathrm{sat} are the masses of the host halo and satellite, v_c is the circular velocity, \sigma is the velocity dispersion, and \rho is the local density. This formula highlights how more massive satellites in denser environments experience shorter merger times, facilitating the coalescence into a single remnant.^[36] The resulting merger remnants exhibit surface brightness profiles well-fitted by the de Vaucouleurs r^{1/4} law, a hallmark of elliptical galaxies, with effective radii encompassing half the total light. N-body simulations demonstrate that these systems often display boxy isophotes, particularly in equal-mass mergers, arising from the orbital structure of the disrupted disks and the triaxiality induced by the interaction.^[35] In gas-rich (wet) mergers, the collision funnels interstellar gas to the center, triggering intense starbursts that contribute to the remnant's stellar mass. Conversely, gas-poor (dry) mergers primarily rearrange existing stars, preserving the kinematic signatures of the progenitors in the final velocity field.^[35] Within the bottom-up hierarchical assembly paradigm, the frequency of major mergers follows a rate dN/dt \propto (1+z)^3 in the local universe, indicating a decline toward lower redshifts as structures stabilize.^[37]

Wet and dry merger differences

Wet mergers, involving gas-rich progenitor galaxies, are characterized by significant gas dissipation during the interaction, which leads to the formation of compact remnants. This dissipation allows for efficient cooling and concentration of baryons toward the center, resulting in tightly bound stellar systems with high central densities. In contrast, dry mergers occur between gas-poor, spheroid-dominated galaxies, where the interaction is primarily collisionless among stars, leading to an adiabatic expansion and growth of the stellar envelope without substantial dynamical friction from gas.^[38] The presence of gas in wet mergers triggers intense bursts of star formation, with star formation rates enhanced by factors of 10 to 100 compared to isolated galaxies, driven by tidal torques that funnel gas into the central regions. This enhanced activity also promotes the formation of supermassive black hole binaries, as the dissipative environment facilitates the rapid sinking and pairing of central black holes. Dry mergers, lacking sufficient gas, produce minimal new stars and result in remnants with redder colors due to the dominance of older stellar populations, with little to no triggered star formation.^[39]^[38] Observationally, wet mergers manifest as post-starburst galaxies exhibiting spectral signatures of both young A-type and older K-type stars, indicative of a recent cessation of star formation following the burst. Dry mergers, on the other hand, contribute to the extended envelopes observed in massive ellipticals within clusters, where the remnants show diffuse, low-surface-brightness outskirts built from accreted stars.^[38] In the evolutionary sequence of massive galaxies, early wet mergers at high redshifts form compact cores through dissipative processes, while subsequent dry mergers at lower redshifts drive the growth of these cores' sizes, following the observed relation where the effective radius scales as r_e \propto M^{0.6}, with minor dry mergers contributing disproportionately to size increase. Numerical simulations indicate that wet mergers were more prevalent at z > 1, dominating the assembly of progenitors, whereas dry mergers become increasingly common at z < 0.5, shaping the present-day population of quiescent ellipticals.

Galaxy Evolution and Quenching

Star formation history

The cosmic star formation rate density (SFRD), denoted as \rho_{\rm SFR}, evolved dramatically over cosmic time, rising steeply with increasing redshift before peaking and then declining toward the present epoch. Observations indicate that \rho_{\rm SFR} \propto (1+z)^{3.5} up to the peak at z \approx 2, corresponding to approximately 3.3 Gyr after the Big Bang, after which it declines by nearly an order of magnitude to z=0. This behavior is encapsulated in the Madau-Lilly relation, an empirical fit derived from integrated ultraviolet (UV) and far-infrared (FIR) luminosity measurements across diverse galaxy populations, and then declining exponentially toward the present epoch with an e-folding timescale of approximately 4 Gyr.^[40] On individual galaxy scales, star formation histories (SFHs) similarly exhibit a peak between z=1 and z=3, reflecting the epoch of "cosmic noon" when galaxy assembly was most intense, followed by an exponential decline in the star formation rate (SFR). For Milky Way-like galaxies, this decline is characterized by a timescale \tau \approx 3 Gyr, consistent with chemical evolution models that trace the buildup of stellar populations from gas accretion and internal processes. A key empirical relation governing local star formation efficiency is the Kennicutt-Schmidt law, which states that the surface density of star formation \Sigma_{\rm SFR} \propto \Sigma_{\rm gas}^{1.4}, where \Sigma_{\rm gas} is the total gas surface density (atomic plus molecular), as derived from resolved observations of nearby star-forming regions.^[41] The cumulative outcome of this history is the present-day integrated stellar mass density, \Omega_* \approx 0.004, representing the fraction of the critical density contributed by stars in galaxies today, as measured from UV and IR surveys that account for dust-obscured star formation. Recent James Webb Space Telescope (JWST) observations at z > 10 reveal an unexpectedly higher SFRD in the early universe than previously extrapolated from lower-redshift data, suggesting a more rapid initial buildup of stellar mass during reionization.^[40]^[42]

Quenching mechanisms

Quenching mechanisms refer to the physical processes that suppress or halt star formation in galaxies, transitioning them from actively star-forming to quiescent states. These mechanisms are crucial for explaining the observed bimodality in galaxy properties, such as the separation between blue, star-forming disk galaxies and red, quiescent spheroids. Broadly, quenching can be driven by internal processes related to a galaxy's mass or morphology, or by external environmental effects in dense structures like clusters. The peak of cosmic star formation at redshift z ≈ 2 provides context for when many quenching processes become prominent, as galaxies begin to exhaust their gas reservoirs post-peak.^[43] Mass quenching primarily affects massive galaxies through feedback from active galactic nuclei (AGN), where energy and momentum from the central supermassive black hole expel or heat interstellar gas, preventing further star formation. In momentum-driven winds, the outflow momentum is proportional to the AGN luminosity divided by the speed of light, p \propto L_{\rm AGN}/c, allowing radiation pressure on dust to drive large-scale gas ejection without requiring energy conservation. This mechanism is particularly effective in galaxies above a stellar mass of about $10^{10.5} M_\odot, where AGN activity correlates with the shutdown of star formation, as evidenced by simulations and observations of outflow velocities matching predicted values. Environmental quenching dominates in group and cluster environments, where interactions with the surrounding medium remove or prevent the replenishment of cold gas. Ram-pressure stripping occurs when the intracluster medium (ICM) exerts a dynamic pressure on a infalling galaxy, given by \rho_{\rm IGM} v^2, exceeding the gravitational restoring force per unit area, \rho_{\rm gas} \sigma^2, leading to the truncation of the gas disk and suppression of star formation. This process, first analytically described for spiral galaxies in clusters, is observed in tails of stripped gas in systems like Virgo and Coma. Complementing this, strangulation halts the inflow of fresh gas from the cosmic web while allowing existing gas to be consumed, resulting in a gradual decline in star formation over gigayears, particularly for satellite galaxies with stellar masses around $10^9--$10^{10} M_\odot. Morphological quenching arises from the structural stability of galactic disks, where a massive stellar spheroid or bulge stabilizes the gas disk against gravitational collapse, inhibiting star formation even if gas is present. This is quantified by the Toomre stability parameter [Q](/page/Q) > 1, indicating that the disk's velocity dispersion and surface density prevent fragmentation into star-forming clouds. In simulated early-type galaxies, this mechanism naturally quenches star formation following bulge growth, linking morphology directly to quiescence without invoking gas removal. Observations confirm that quiescent galaxies with extended gas disks exhibit [Q](/page/Q) values above unity, supporting this as a key internal process.^[43] Observationally, quenching manifests as a transition through the "green valley" in color-magnitude diagrams, where galaxies evolve from blue to red over timescales of 0.5--2 Gyr, depending on morphology and environment. Spectroscopic studies of green valley galaxies reveal declining star formation rates consistent with rapid central quenching in early-types versus slower disk-wide suppression in late-types. Post-merger quenching often involves morphological transformation to early-type systems, where mergers drive gas inflows that fuel a final starburst before stabilizing the remnant into a quiescent spheroid. Recent analyses of post-merger samples show elevated quenching fractions within 500 Myr of coalescence, attributed to the combined effects of gas exhaustion and structural stabilization.^[44]^[45]

Environmental Influences

Group and cluster environments

Galaxies in dense group and cluster environments experience accelerated evolution compared to those in the field, primarily through interactions that suppress star formation and alter morphologies. Observations indicate that at redshifts z < 1, the fraction of red, quiescent galaxies (f_{\rm red}) in clusters reaches approximately 0.8 within the virial radius, significantly higher than the roughly 0.3 observed in comparable field populations of similar stellar mass. This environmental dependence highlights the role of cluster-specific processes in driving the buildup of passive galaxy populations over cosmic time. Galactic conformity manifests in these environments as a correlation where quenched satellite galaxies are more likely to surround passive central galaxies, extending beyond immediate pairwise interactions to scales influenced by the shared large-scale structure. This phenomenon suggests that the large-scale environment imprints quenching signatures on neighboring galaxies, with studies showing suppressed star formation in satellites around passive centrals in groups and clusters. In rich groups, particularly those with halo masses above $10^{13} \, M_\odot, galaxies exhibit reduced star formation rates (SFRs) by factors of up to 0.5 relative to poorer groups, alongside older stellar populations, as evidenced by recent analyses of local universe samples.^[46] These effects scale with group richness, emphasizing the cumulative impact of multiple interactions in denser halos. Harassment, involving repeated high-speed encounters between galaxies in clusters, heats stellar disks and disrupts morphologies without significantly removing gas reservoirs, leading to the transformation of spirals into lenticulars or dwarfs.^[47] This process is particularly effective for infalling galaxies on radial orbits that penetrate the cluster core, where cumulative dynamical heating truncates disks and scatters stars, contributing to the observed excess of early-type galaxies in cluster outskirts.^[48] Pre-processing occurs when galaxies evolve within infalling groups before reaching the cluster core, where group-scale interactions quench star formation and modify stellar populations prior to exposure to stronger cluster processes. Simulations and observations reveal that up to 50% of cluster members undergo this pre-processing, resulting in higher quenched fractions among group satellites compared to field galaxies at similar redshifts.^[49] This staged environmental influence underscores how hierarchical assembly amplifies quenching efficiency in bound structures.

Cosmic web effects

The cosmic web, comprising filaments, walls, and voids, exerts significant influence on galaxy formation by modulating gas accretion and angular momentum acquisition. Galaxies embedded within this large-scale structure experience anisotropic inflows that shape their morphological and dynamical evolution. Filaments, as the primary conduits of matter, channel cold gas streams toward massive halos, facilitating sustained star formation in early galaxies. These streams, penetrating shock-heated halos, deliver pristine, low-entropy gas that fuels disk growth and mergers. In filamentary environments, galaxy and halo spins tend to align preferentially with the filament axis, particularly for lower-mass systems. This alignment arises from tidal torques exerted by the surrounding web, where the direction of angular momentum acquisition correlates with the filament's orientation. Observational evidence supports this, with galaxy spin vectors showing coherence along cosmic filaments at redshifts up to z ≈ 1.^[50] Galaxies residing in cosmic voids, the underdense regions of the web, exhibit distinct properties due to their isolation from filamentary inflows. These void galaxies display higher specific star formation rates (SSFRs), bluer colors, and greater gas richness compared to those in denser environments at fixed stellar mass. The reduced environmental interactions in voids lead to lower quenching rates, allowing prolonged star formation activity.^[51] Across the cosmic web, an environmental quenching gradient manifests, with the quenched fraction of galaxies (f_\mathrm{quenched}) increasing progressively from voids to filaments, walls, and clusters. In voids, f_\mathrm{quenched} remains low (typically <20% for low-mass galaxies), reflecting minimal external suppression, while it rises sharply in cluster environments due to enhanced gas stripping and heating. This gradient highlights the web's role in modulating star formation cessation over cosmic scales. Numerical simulations reveal that the anisotropy of the cosmic web directly impacts halo spin properties. The spin parameter \lambda, a measure of rotational support, correlates with the angle \theta between the halo spin vector and the nearest filament, following \lambda \propto \cos \theta for low-mass halos. This relation stems from preferential accretion along filaments, which amplifies spin for aligned systems while suppressing it for perpendicular orientations. In hydrodynamical simulations like IllustrisTNG, this effect persists across redshifts, influencing disk stability and morphology.^[52] Recent observational studies from 2025 link a galaxy's position within the cosmic web to assembly bias, where halos of similar mass but different formation histories cluster differently based on their web environment. Galaxies near filaments show enhanced assembly bias, with earlier-forming systems exhibiting stronger clustering, as traced by spectroscopic surveys like DESI. This web-dependent bias underscores how large-scale structure imprints on galaxy evolution beyond local density effects.

Numerical Simulations

N-body and semi-analytic models

N-body simulations form the backbone of collisionless modeling in galaxy formation, focusing primarily on the gravitational dynamics of dark matter particles. These simulations employ algorithms such as particle-mesh (PM) or tree codes to solve the Poisson equation, \nabla^2 \Phi = 4\pi G \rho, which governs the gravitational potential \Phi due to the dark matter density \rho. Particle-mesh methods discretize the density on a grid to compute long-range forces efficiently via fast Fourier transforms, while tree codes, like the Barnes-Hut algorithm, hierarchically group particles into octrees for approximating short-range interactions, enabling scalable computations for billions of particles.^[53]^[54] A landmark example is the Millennium Simulation, conducted in 2005, which evolved $10^{10} dark matter particles within a comoving box of side length 500 h^{-1} Mpc from high redshift to the present day, assuming a \LambdaCDM cosmology. This simulation resolved structures down to halo masses of approximately $10^9 M_\odot and accurately predicted the halo mass function, providing merger trees that underpin subsequent galaxy formation studies.^[55] By tracing the hierarchical assembly of dark matter halos, N-body simulations like Millennium reveal the statistical properties of the cosmic web, such as the abundance and clustering of halos, which serve as the scaffolds for galaxy populations.^[56] Semi-analytic models (SAMs) build upon N-body outputs by applying simplified, analytic prescriptions to model baryonic processes on the resulting dark matter merger trees. These models track the flow of gas and the formation of stars across cosmic time, using recipes for processes like radiative cooling, where the cooling timescale is approximated as t_{\rm cool} = \frac{3/2 \, k T}{n \Lambda}, with n the gas density, T the temperature, k Boltzmann's constant, and \Lambda the cooling function. For instance, the L-Galaxies model, an evolution of the GALFORM framework, incorporates updated treatments of gas accretion, star formation, and feedback, with recent 2024 refinements improving matches to observed galaxy scaling relations from z \approx 0 to z \approx 10.^[57] Similarly, the GALFORM model predicts galaxy luminosity functions by parameterizing dynamical processes during mergers, achieving good agreement with observations up to z < 6. The GAEA (GAlaxy Evolution and Assembly) model extends this approach with detailed tracking of chemical enrichment and environmental effects, yielding predictions for galaxy clustering and quenching that align with surveys like SDSS.^[58] SAMs efficiently generate large mock catalogs of galaxy properties, facilitating comparisons with observations, but they inherently ignore full hydrodynamics by approximating baryonic effects through tunable recipes. This simplification limits their fidelity in capturing complex gas dynamics, though extensions can incorporate hydrodynamic insights from separate simulations.^[59]

Full hydrodynamical simulations

Full hydrodynamical simulations of galaxy formation and evolution solve the equations of ideal fluid dynamics coupled to self-gravity to model the behavior of baryonic gas alongside dark matter. These simulations treat gas as a compressible fluid governed by the Euler equations of hydrodynamics, which describe mass conservation,

\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0,

and momentum conservation,

\frac{\partial (\rho \mathbf{v})}{\partial t} + \nabla \cdot (\rho \mathbf{v} \mathbf{v} + P) = -\rho \nabla \Phi,

where \rho is the gas density, \mathbf{v} is the velocity field, P is the pressure, and \Phi is the gravitational potential. The gravitational potential \Phi is evolved using N-body methods for collisionless components like dark matter and stars, while gas dynamics are discretized using various numerical schemes. Common approaches include smoothed particle hydrodynamics (SPH), as implemented in the GADGET code, which represents gas as discrete particles with smoothed kernels; adaptive mesh refinement (AMR) methods, such as in the Enzo code, which solve the equations on a hierarchical grid; and moving-mesh techniques, like those in the AREPO code, which use a dynamic Voronoi tessellation to follow fluid elements with minimal numerical diffusion. Modern suites of full hydrodynamical simulations have achieved significant realism in reproducing observed galaxy properties by running large cosmological volumes over cosmic time. The IllustrisTNG project, initiated in 2018 with ongoing updates through 2025, employs the AREPO code to simulate boxes up to 300 Mpc on a side with initial gas masses of $1.0 \times 10^6 M_\odot h^{-1} for TNG100 and $7.4 \times 10^6 M_\odot h^{-1} for TNG300.^[60]^[61] Similarly, the EAGLE simulations use a modified SPH solver in GADGET to model a 100 Mpc box with initial gas particle masses of approximately $1.8 \times 10^6 M_\odot in their reference run, enabling detailed studies of galaxy assembly.^[62] The SIMBA suite, based on the meshless finite-mass method in the GIZMO code, simulates a 100 Mpc/h box with initial gas particle masses of $6.6 \times 10^6 M_\odot h^{-1}, focusing on black hole co-evolution with galaxies.^[63] These simulations incorporate subgrid models for unresolved processes like star formation and feedback to bridge scales between individual particles and cosmological volumes.^[63] Key outputs from these simulations include realistic galaxy morphologies, such as disks and bulges emerging from hierarchical merging, and accurate black hole growth histories driven by gas accretion. For instance, IllustrisTNG successfully predicts the observed correlation between supermassive black hole mass M_\mathrm{BH} and stellar velocity dispersion \sigma, with simulated galaxies matching empirical relations across a wide mass range. EAGLE and SIMBA similarly reproduce morphological diversity and the quenching of star formation in massive systems, providing insights into the baryonic Tully-Fisher relation and environmental effects on galaxy evolution.^[62]^[63] Despite advances, full hydrodynamical simulations face challenges like the overcooling catastrophe, where gas cools too rapidly in massive halos without sufficient heating, leading to unrealistically high star formation rates; this is mitigated by implementing stellar and active galactic nucleus (AGN) feedback to reheat or eject gas. Computational demands are substantial, with runs like IllustrisTNG-100 requiring tens of millions of CPU hours on supercomputers to evolve from high redshift to the present day.^[64] Recent developments include the 2025 P-Millennium simulation, which couples the GAEA semi-analytic model to high-resolution N-body merger trees for detailed galaxy assembly histories, complementing full hydro approaches by efficiently exploring parameter variations.^[65]

Baryonic Physics in Models

Gas thermodynamics and cooling

In galaxy formation simulations, the thermal evolution of gas is governed by heating and cooling processes that determine its ability to collapse and form structures. Radiative cooling functions, denoted as \Lambda(T, n), quantify the energy loss rate per unit volume as n^2 \Lambda(T, n), where T is temperature and n is density. For primordial hydrogen-helium gas, cooling primarily occurs through processes like Ly\alpha emission, collisional excitation of hydrogen and helium, and recombination, with \Lambda \approx 10^{-23} erg cm^3 s^{-1} at T = 10^4 K.^[66] These functions are typically tabulated or approximated using detailed atomic physics calculations for implementation in hydrodynamical codes.^[66] Metallicity significantly enhances cooling efficiency, particularly at temperatures around $10^4 K, where fine-structure lines of metals such as [O II] and [C II] dominate the energy loss. For metallicities Z > 0.01 Z_\odot, metal-line cooling boosts the total \Lambda by a factor of approximately 10 compared to primordial gas, enabling more rapid thermal relaxation and fragmentation in enriched environments.^[66] This enhancement arises from the excitation and de-excitation of low-lying atomic levels in ions like O^+ and C^+, which radiate efficiently in the interstellar medium. Heating of the gas counterbalances cooling and maintains thermal equilibrium in diffuse regions. Photoionization by the cosmic ultraviolet background, characterized by intensity J_{21} \approx 1 (in units of $10^{-21} erg cm^{-2} s^{-1} Hz^{-1} sr^{-1}) at the Lyman limit, ionizes hydrogen and heats the intergalactic medium to an equilibrium temperature of approximately $10^4 K. This background, primarily from quasars and early star-forming galaxies, prevents excessive cooling in low-density gas and influences the thermal state during hierarchical structure formation.^[67] The Jeans mass, M_J \propto T^{3/2} / \rho^{1/2}, sets the characteristic scale for gravitational instability and fragmentation of cooling gas clouds, where \rho is density. At high redshifts, where gas temperatures are elevated due to virial heating in collapsing halos, M_J determines the minimum mass for collapse, promoting fragmentation into substructures that seed galaxy formation. Simulations of galaxy formation often incorporate multiphase models for the interstellar medium (ISM) to capture its thermal structure, including cold molecular cores (T \lesssim 100 K), warm neutral hydrogen (HI) phases (T \sim 10^3 - 10^4 K), and hot halo gas (T > 10^5 K). These models treat the ISM as a composite of phases in pressure equilibrium, with cooling driving phase transitions and hot gas providing a reservoir for accretion.^[68] Such sub-resolution approaches enable realistic thermal evolution without resolving all microscopic scales.^[68]

Star formation and feedback

In subgrid models of galaxy formation simulations, star formation is prescribed as a local process occurring in dense, molecular gas phases, often following a recipe where the star formation rate density \rho_\mathrm{SFR} is proportional to the gas density \rho_\mathrm{gas} divided by the local free-fall time t_\mathrm{ff}, such that \rho_\mathrm{SFR} = \epsilon_* \rho_\mathrm{gas} / t_\mathrm{ff}. Here, \epsilon_* represents the star formation efficiency per free-fall time, calibrated to values between 0.01 and 0.1 to match observed molecular cloud efficiencies and global galaxy scaling relations. The free-fall time is defined as t_\mathrm{ff} = \left(3\pi / (32 G \rho_\mathrm{gas})\right)^{1/2}, reflecting the dynamical timescale for gravitational collapse in self-gravitating gas clouds. This formulation, rooted in theoretical expectations for turbulent, Jeans-limited fragmentation, is implemented in major simulation suites like EAGLE and FIRE after thresholding gas densities above a critical value, typically around 10-100 cm^{-3}, to mimic the shift to molecular conditions. Stellar feedback, which couples the energy and momentum from newly formed stars back into the interstellar medium (ISM), takes several forms to counteract rapid gas cooling and maintain realistic galaxy morphologies. Type II supernovae (SNe), exploding from massive stars with lifetimes of ~10 Myr, provide the dominant early feedback, with approximately one event per 100 M_\odot of stars formed according to the initial mass function (IMF).^[69] Each SN injects ~10^{51} erg of thermal energy, heating surrounding gas and driving shocks, while momentum injection reaches p_* \approx 3000 , \mathrm{km , s^{-1}} \times (M_* / 100 , M_\odot), where M_* is the stellar mass formed, arising from radiation pressure and SN shell momentum.^[70] Later phases include asymptotic giant branch (AGB) stellar mass loss, returning ~30-40% of stellar mass as enriched winds over 100-500 Myr, and radiation pressure from young stars, which imparts additional momentum (up to ~10^{36} dyn s per 100 M_\odot) via UV photons absorbed by dust.^[71] These mechanisms are stochastically distributed in simulations to approximate unresolved clustering. The combined effects of this feedback are crucial for regulating galaxy evolution by preventing overcooling, where dense gas would otherwise collapse unchecked into stars at efficiencies far exceeding observations. Thermal and momentum injections heat and pressurize the ISM, launching outflows with characteristic velocities of ~300 km s^{-1} that remove low-angular-momentum gas and suppress excessive central star formation.^[72] Metal enrichment from SN ejecta and AGB returns builds ISM abundances to near-solar levels, Z \approx 0.02 (Z_\odot), facilitating dust formation and further momentum transfer while tracing chemical evolution. Such models reproduce key observations, including bipolar galactic winds in local starburst galaxies like M82 and NGC 253, where outflow velocities and mass-loading rates align with detected blueshifted absorption lines and X-ray emitting plasma.^[73]

Supermassive black holes and AGN

In galaxy formation simulations, supermassive black holes (SMBHs) are typically initialized as seed black holes at high redshifts to model their early growth and influence on host galaxies. Common seeding mechanisms include direct collapse of pristine, metal-poor gas clouds in protogalaxies, which can produce seeds with masses around $10^5 M_\odot at z > 10, or remnants from the mergers of stellar-mass black holes formed by the first generations of Population III stars.^[74]^[75] These initial masses are chosen to reconcile the rapid appearance of billion-solar-mass SMBHs observed at z \approx 6-7 with theoretical growth limits, avoiding overly light seeds (\sim 10^2 M_\odot) that would require sustained super-Eddington accretion to match observations.^[76] SMBH growth in simulations is primarily driven by gas accretion, often limited by the Eddington rate to prevent unphysically rapid mass buildup. The Eddington accretion rate is given by \dot{M} = L_E / (\epsilon c^2), where L_E = 1.25 \times 10^{38} (M_\mathrm{BH}/M_\odot) erg s^{-1} is the Eddington luminosity, \epsilon \approx 0.1 is the radiative efficiency for a standard thin accretion disk, and c is the speed of light; this yields a characteristic e-folding (or doubling) timescale of about 45 Myr at the Eddington limit. In practice, accretion is modeled using the Bondi-Hoyle-Lyttleton formalism for spherically symmetric inflow from the surrounding hot circumgalactic medium, with \dot{M}_B \propto M_\mathrm{BH}^2 \rho / c_s^3, where \rho is the gas density and c_s the sound speed, though resolution limitations often require boosting factors to capture unresolved angular momentum transport.^[77] Mergers of SMBHs during galaxy interactions contribute additional growth, particularly at lower redshifts, but accretion dominates the mass assembly for most systems.^[76] AGN feedback from accreting SMBHs is implemented in two primary modes to regulate star formation and gas cooling in host galaxies. In the quasar (radiative) mode, active during high-accretion phases near the Eddington limit, energy is injected isotropically as thermal or radiative feedback with efficiency \epsilon_f \approx 0.05 L_\mathrm{bol}, where L_\mathrm{bol} is the bolometric luminosity, heating surrounding gas and potentially driving outflows that quench star formation in massive galaxies. The radio (kinetic or jet) mode activates at lower accretion rates (\dot{M} < 0.01 \dot{M}_\mathrm{Edd}), coupling momentum p \approx 20 L_\mathrm{bol} / c to the interstellar or circumgalactic medium via bipolar jets, which inflate bubbles and maintain hot halo atmospheres to suppress cooling flows in clusters.^[78] These modes are calibrated against observations of AGN outflows and cavity energies, ensuring realistic impact on galaxy evolution without over-ejecting baryons.^[79] The co-evolution of SMBHs and their host galaxies is captured through tight scaling relations in simulations, such as M_\mathrm{BH} \approx 0.001 M_\mathrm{bulge}, which emerges naturally from feedback-regulated growth where AGN activity limits bulge assembly by expelling or heating inflowing gas.^[80] This relation holds across cosmic time, with SMBHs growing in tandem with the stellar mass of spheroidal components, and AGN feedback playing a key role in quenching star formation once the relation is established, particularly in massive systems above M_* \sim 10^{10} M_\odot.^[81]

Radiative processes and other physics

In galaxy formation models, radiative processes are incorporated through radiation hydrodynamics (RHD), which solves the radiation transport equations to capture the propagation of photons from stars and other sources, interacting with gas via absorption, scattering, and emission.^[82] These simulations typically employ moment-based methods or Monte Carlo techniques to approximate the radiative transfer, accounting for opacities dominated by dust in the interstellar medium, where values around κ ≈ 1 cm²/g are common for ultraviolet and optical wavelengths. The inclusion of RHD significantly influences galaxy evolution by regulating gas heating and ionization, particularly during cosmic reionization, where radiative feedback from early stars and galaxies ionizes the intergalactic medium, altering the thermal state of baryons and suppressing low-mass galaxy formation.^[83] Magnetic fields play a crucial role in non-gravitational physics, amplified through turbulent dynamo processes in the collapsing gas of forming galaxies, following an exponential growth law B ∝ exp(γ t) with growth rates γ ≈ 10 Gyr⁻¹ driven by small-scale turbulence in the interstellar medium.^[84] In magnetohydrodynamical (MHD) simulations, such as those using the AREPO code, the Lorentz force from these fields resists gravitational collapse, shapes gas flows, and influences angular momentum transport, leading to more realistic disk morphologies and reduced fragmentation in protogalaxies. Primordial seed fields, often initialized at levels of 10⁻¹² G, are amplified to observed galactic strengths of μG over cosmic time, with dynamo action becoming efficient once sufficient turbulence is generated by hierarchical structure formation.^[85] Cosmic rays (CRs), high-energy particles accelerated by supernovae and other processes, contribute additional pressure and heating in galaxy models through their transport and interactions with the gas. CR pressure is given by P_CR ≈ e_CR / 3, where e_CR is the CR energy density, providing non-thermal support comparable to thermal gas pressure in dense regions and driving outflows via streaming instabilities that excite Alfvén waves, transferring energy to heat the plasma. In simulations, CRs modify the multiphase structure of the interstellar medium by heating diffuse gas and suppressing cooling, which reduces the amount of star-forming material and enhances galactic winds. Other physical processes, such as viscosity and thermal conduction, are often modeled at subgrid scales to account for unresolved microphysics in large-scale simulations. Viscosity, including artificial and physical forms, damps small-scale turbulence and prevents numerical instabilities in shocks, while Spitzer conduction transfers heat along magnetic field lines, smoothing temperature gradients in the circumgalactic medium.^[86] These are implemented with limiters to avoid over-dissipation, ensuring energy conservation in feedback-dominated environments. Recent simulations from 2023 onward have increasingly incorporated CR transport alongside these effects, yielding more realistic interstellar medium structures and reducing star formation rates by up to 20% in dwarf galaxies through enhanced gas heating and expulsion.

References

[1]
Galaxy Basics - NASA Science
Most galaxies are between 10 billion and 13.6 billion years old. Some are almost as old as the universe itself, which formed around 13.8 billion years ago. ...Galaxy Evolution · Galaxy Types · Galaxies Stories · Large Scale StructuresMissing: formation | Show results with:formation
[2]
Galaxy Evolution - NASA Science
Oct 22, 2024 · Astronomers know that some of the oldest galaxies observed formed when the universe was only about a billion years old. They think those ...Galaxy Evolution · Formation · Collisions
[3]
Herschel's chronicles of galaxy evolution - ESA Science & Technology
Sep 20, 2017 · Star formation in galaxies takes place within dense clouds of gas that, for most of cosmic history, also contain small amounts of dust. Newborn ...
[4]
Galaxies Over Time - NASA Science
Sep 28, 2023 · Computer models that scientists have made to understand galaxy formation indicate that galaxies are created when dark matter merges and clumps ...
[5]
[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.
[6]
[astro-ph/9902027] Summary of Results from COBE - arXiv
Feb 1, 1999 · The results include a precise measurement of the Cosmic Microwave Background (CMB) radiation intensity, discovery and maps of the CMB anisotropy ...
[7]
Planck 2018 results - VI. Cosmological parameters
We present cosmological parameter results from the final full-mission Planck measurements of the cosmic microwave background (CMB).Missing: n_s sigma_8
[8]
Spectral Classification of Galaxies Along the Hubble Sequence - arXiv
Jul 12, 1995 · We develop a straightforward and quantitative two-step method for spectroscopically classifying galaxies from the low signal-to-noise (S/N) optical spectra.Missing: seminal | Show results with:seminal
[9]
The Luminosity Function of Galaxies in SDSS Commissioning Data
The result is well-described by a Schechter function with parameters ϕ* = (1.46 ± 0.12) × 10-2 h3 Mpc-3, M* = -20.83 ± 0.03, and α = -1.20 ± 0.03. The ...<|separator|>
[10]
size distribution of galaxies in the Sloan Digital Sky Survey
Galaxies more massive than M0 have and σln R∼ 0.3, while less massive galaxies have and σln R∼ 0.5. For early-type galaxies, the relation is significantly ...
[11]
Determining the escape fraction of ionizing photons during ...
Observations of the Gunn–Peterson optical depth in the high-redshift Lyα forest imply a very highly ionized IGM at z≲ 6, suggesting that reionization was ...Missing: trough | Show results with:trough
[12]
Two Remarkably Luminous Galaxy Candidates at z ≈ 10–12 ...
Nov 17, 2022 · Here we present a search for luminous z > 10 galaxies with JWST/NIRCam photometry spanning ≈1–5 μm and covering 49 arcmin 2 from the public JWST Early Release ...
[13]
JADES Imaging of GN-z11: Revealing the Morphology and ...
We find that GN-z11 has a stellar mass of 109.1±0.3 M⊙, is actively forming stars (sSFR = 10−7.7±0.3 yr−1) with a young stellar age (half-mass time ...Missing: Msun | Show results with:Msun
[14]
(PDF) Primordial Rotating Disk Composed of ≥15 Star Forming ...
Aug 9, 2025 · The clumpiness of the Cosmic Grapes significantly exceeds that of galaxies at later epochs and the predictions from current simulations for ...
[15]
The Early Emergence of Massive Galaxies and Clusters of Galaxies
Nov 12, 2024 · Galaxies in the early Universe appear to have grown too big too fast, assembling into massive, monolithic objects more rapidly than ...
[16]
EPOCHS VI: the size and shape evolution of galaxies since z ∼ 8 ...
We present the results of a size and structural analysis of 1395 galaxies at 0.5 ≤ z ≲ 8 with stellar masses log (M*/M⊙)> 9.5 within the James Webb Space ...
[17]
https://www.researchgate.net/publication/388636794_Primordial_Rotating_Disk_Composed_of_15_Star_Forming_Clumps_at_Cosmic_Dawn
[18]
https://iopscience.iop.org/article/10.3847/1538-4357/ad834d
[19]
The New Galaxy: Signatures of Its Formation - Ken Freeman & Joss ...
ELS are widely viewed as advocating a smooth monolithic collapse of the protocloud with a timescale of order 108 years. But Sandage (1990) stresses that this is ...
[20]
Monolithic View of Galaxy Formation and Evolution - MDPI
Since in Cold Dark Matter (CDM) cosmogony small scale structures typically collapse early and form low-mass haloes that subsequently can merge to assembly ...Missing: incompatible | Show results with:incompatible
[21]
[astro-ph/0111097] Star Streams in the Milky Way Halo - arXiv
Nov 5, 2001 · The last 10-20 years has seen a profound shift in views of how the Galaxy's halo formed. The idea of a monolithic early collapse of a single ...
[22]
https://ned.ipac.caltech.edu/level5/Sept16/Freeman/Freeman1.html
[23]
Accelerated formation of ultra-massive galaxies in the first billion years
Nov 13, 2024 · Here we report a systematic study of 36 massive dust-obscured galaxies with spectroscopic redshifts between 5 and 9 from the JWST FRESCO survey.Missing: red monsters paper
[24]
Formation and rotation of disc galaxies with haloes - Oxford Academic
(2) The halo rotation curves, computed on the assumption of detailed conservation of angular momentum in the discs, are fairly fiat. (3) Most discs are globally ...
[25]
rotation curves of spiral galaxies
Three dashed lines are averaged rotation curves of galaxies with maximum velocities greater than 200 km s−1 , between 200 and 250 km s−1 , and below 200 km s−1 ...
[26]
[1909.07118] The Inside-out Growth of the Galactic Disk - arXiv
Sep 16, 2019 · Our model predicts the half-mass and half-light sizes of the Galactic disk at earlier epochs, which can be compared to the observed redshift - ...Missing: R_d proportional t^{ 0.5} gas infall
[27]
The Inside-out Growth of the Galactic Disk - IOPscience
Oct 15, 2019 · We quantify the inside-out growth of the Milky Way's low-α stellar disk, modeling the ages, metallicities, and Galactocentric radii of APOGEE red clump stars.
[28]
A core in a star-forming disc as evidence of inside-out growth in the ...
Oct 11, 2024 · This evidence suggests a detection of the inside-out growth of a galaxy as a proto-bulge and a star-forming disc in the epoch of reionization.
[29]
https://ned.ipac.caltech.edu/level5/Sept16/Sofue/Sofue3.html
[30]
Evolution of violent gravitational disc instability in galaxies
The gas and stars flow inward through the disc to a central bulge due to torques that drive angular momentum outwards. The gravitational energy released by the ...<|separator|>
[31]
[astro-ph/0301385] Major Mergers and the Origin of Elliptical Galaxies
Jan 20, 2003 · The formation of elliptical galaxies as a result of the merging of spiral galaxies is discussed. We analyse a large set of numerical N-Body ...
[32]
Dynamical friction and galaxy merging time-scales - Oxford Academic
In this paper, we study the merging time-scales of extended dark matter haloes using N-body simulations. We compare these results to standard estimates.INTRODUCTION · SATELLITE DECAY FROM... · COMPUTING MERGER TIME...
[33]
MAJOR-MERGER GALAXY PAIRS IN THE COSMOS FIELD—MASS ...
Feb 16, 2012 · The best-fitting function for the merger rate is Rmg (Gyr−1) = 0.053 × (Mstar/1010.7 M☉)0.3(1 + z)2.2/(1 + z/8). This rate implies that galaxies ...
[34]
THE ROLE OF DRY MERGERS FOR THE FORMATION AND ...
We demonstrate that the brightest cluster galaxy (BCG) evolves away from the Kormendy relation as defined by the smaller mass galaxies.
[35]
star formation triggering in interactions - from mergers to starbursts
The triggering of star formation in merger largely results from tidal inflows driven by the interaction, which fuels a relatively concentrated or nuclear ...
[36]
[1403.0007] Cosmic Star Formation History - arXiv
Feb 28, 2014 · Here we review the range of complementary techniques and theoretical tools that allow astronomers to map the cosmic history of star formation.
[37]
[astro-ph/9712213] The Global Schmidt Law in Star Forming Galaxies
Dec 16, 1997 · The Schmidt law provides a surprisingly tight parametrization of the global star formation law, extending over several orders of magnitude in SFR and gas ...
[38]
Probing the Origin of the Star Formation Excess Discovered by ...
Nov 18, 2024 · Remarkably, as shown in Figure 2, the JWST revealed that the SFR density beyond z ≳ 10 is higher than the extrapolation of the pre-JWST results ...
[39]
MORPHOLOGICAL QUENCHING OF STAR FORMATION
ABSTRACT. We point out a natural mechanism for quenching of star formation in early-type galaxies (ETGs). It automatically links the color of a ...
[40]
green valley is a red herring: Galaxy Zoo reveals two evolutionary ...
We have shown that the star formation histories of late-type galaxies are consistent with a very gradual quenching, corresponding to an exponential time-scale ...
[41]
Post-merger quenching peaks within 500 Myr of coalescence - arXiv
Oct 8, 2024 · Our results demonstrate that 1) galaxy-galaxy interactions can lead to rapid post-merger quenching within 0.5 Gyr of coalescence, 2) the majority of (but not ...
[42]
Galaxy Evolution in the Local Universe: Group Richness Effects on ...
Sep 16, 2025 · (2024). Fernandez et al. (2024) observed a reduced star formation activity and older stellar population in rich groups when compared to poorer ...
[43]
Galaxy harassment and the evolution of clusters of galaxies - Nature
Feb 15, 1996 · Here we propose that multiple highspeed encounters between galaxies—'galaxy harassment'— drives the morphological evolution in clusters.
[44]
sensitivity of harassment to orbit: mass loss from early-type dwarfs in ...
We find harassment is only effective at stripping stars or truncating their stellar discs for orbits that enter deep into the cluster core. Comparing to the ...
[45]
The Preprocessing of Galaxies in the Early Stages of Cluster ...
Jan 31, 2022 · The suppression of star formation (“quenching”) occurs in the infall regions of these structures even before the galaxies enter the denser group ...
[46]
Spin alignments within the cosmic web: a theory of constrained tidal ...
This misalignment implies in particular an azimuthal orientation for the spins of more massive galaxies and a spin alignment with the filament for less massive ...<|separator|>
[47]
ON THE STAR FORMATION PROPERTIES OF VOID GALAXIES
The number of red galaxies within the walls is reduced, but there is still a higher ratio of red to blue galaxies in walls than voids. For clarity, the ...
[48]
Connecting the cosmic web to the spin of dark haloes
We investigate the alignment of the spin of dark matter haloes relative (i) to the surrounding large-scale filamentary structure, and (ii) to the tidal tensor ...
[49]
TreePM: A code for Cosmological N-Body Simulations - arXiv
Nov 3, 1999 · We describe the TreePM method for carrying out large N-Body simulations to study formation and evolution of the large scale structure in the Universe.
[50]
tree particle-mesh: an adaptive, efficient, and parallel code for ...
An improved implementation of an N-body code for simulating collisionless cosmological dynamics is presented. TPM (tree particle-mesh) combines the PM ...
[51]
Simulations of the formation, evolution and clustering of galaxies ...
Here we present a simulation of the growth of dark matter structure using 2,1603 particles, following them from redshift z = 127 to the present in a cube-shaped ...
[52]
Millennium Simulation Project - MPA Garching
The Millennium Run used more than 10 billion particles to trace the evolution of the matter distribution in a cubic region of the Universe over 2 billion light ...Movies of the simulation · Pictures of the galaxy distribution
[53]
Probing galaxy evolution from z = 0 to z ≃ 10 through galaxy scaling ...
We present a comprehensive examination of the three latest versions of the L-galaxies semi-analytic galaxy formation model, focusing on the evolution of galaxy ...
[54]
The rise of active galactic nuclei in the galaxy evolution ... - NASA ADS
We present a new implementation of the GAlaxy Evolution and Assembly (GAEA) semi-analytic model, that features an improved modelling of the process of cold gas ...
[55]
GAEA - Google Sites
The GAEA (GAlaxy Evolution and Assembly) semi-analytic model is a state-of-the-art theoretical tool that simulates galaxy formation and evolution in ...
[56]
IllustrisTNG - Main
The IllustrisTNG project is an ongoing series of large, cosmological magnetohydrodynamical simulations of galaxy formation.Images + Videos · Results · Project Description · Data Access
[57]
First results from the IllustrisTNG simulations: the stellar mass ...
The IllustrisTNG project is a new suite of cosmological magnetohydrodynamical simulations of galaxy formation performed with the arepo code and updated ...
[58]
EAGLE project: simulating the evolution and assembly of galaxies ...
Table 2 lists box sizes and resolutions of the main EAGLE simulations. All simulations were run to redshift z = 0. Note that contrary to convention, box sizes, ...
[59]
[1901.10203] Simba: Cosmological Simulations with Black Hole ...
Jan 29, 2019 · We introduce the Simba simulations, the next generation of the Mufasa cosmological galaxy formation simulations run with Gizmo's meshless finite mass ...
[60]
About - Illustris simulation
The largest was run on 8,192 compute cores, and took 19 million CPU hours ... Hydrodynamical cosmological simulations, due to their high computational cost ...
[61]
Galaxy Assembly and Evolution in the P-Millennium simulation - arXiv
Sep 3, 2024 · We present results from the latest version of the GAEA theoretical model of galaxy formation coupled with merger trees extracted from the Planck Millennium ...
[62]
https://academic.oup.com/mnras/article/446/1/521/1316115
[63]
Cosmological smoothed particle hydrodynamics simulations
Springel (2000) tried to overcome this problem by developing a feedback method based on an effective model for the star-forming interstellar medium (ISM), where ...
[64]
Supernova feedback in numerical simulations of galaxy formation
We present a novel implementation of six sub-grid SN feedback schemes in the moving-mesh code Arepo, including injections of thermal and/or kinetic energy.
[65]
[1102.1446] Maximally Star-Forming Galactic Disks I. Starburst ...
Feb 7, 2011 · If the momentum feedback per stellar mass formed is p_*/m_* ~ 3000 km/s, the predicted star formation rate is Sigma_SFR=2 pi G Sigma^2 m_*/p_* ...Missing: M_sun | Show results with:M_sun
[66]
Modeling for Stellar Feedback in Galaxy Formation Simulations
Various heuristic approaches to model unresolved supernova (SN) feedback in galaxy formation simulations exist to reproduce the formation of spiral galaxies ...Missing: seminal | Show results with:seminal
[67]
How to model supernovae in simulations of star and galaxy formation
Abstract. We study the implementation of mechanical feedback from supernovae (SNe) and stellar mass loss in galaxy simulations, within the Feedback In Real.
[68]
Starburst-driven galactic winds — I. Energetics and intrinsic X-ray ...
Galactic winds are unambiguously detected in many local edge-on starburst galaxies (Lehnert & Heckman 1996), and their presence can even be inferred in ...
[69]
Direct Collapse Black Hole Formation with Realistic Lyman–Werner ...
We study the birth of supermassive black holes from the direct collapse process and characterize the sites where these black hole seeds form.Missing: seminal | Show results with:seminal
[70]
Massive black hole seeds born via direct gas collapse in galaxy ...
Abstract. We study the statistics and cosmic evolution of massive black hole seeds formed during major mergers of gas-rich late-type galaxies. Generalizing.Missing: seminal | Show results with:seminal
[71]
Supermassive black holes in cosmological simulations I: M BH
ABSTRACT. The past decade has seen significant progress in understanding galaxy formation and evolution using large-scale cosmological simulations. While t.
[72]
Bondi or not Bondi: the impact of resolution on accretion and drag ...
In this paper, we investigate the impact of resolution on two commonly used SMBH sub-grid algorithms; the Bondi–Hoyle–Lyttleton (BHL) formula for accretion on ...
[73]
case for large-scale AGN feedback in galaxy formation simulations
In this paper, we therefore explore a range of empirical active galactic nucleus (AGN) feedback models, within the FABLE simulation suite, with different ...
[74]
The co-evolution of supermassive black holes and galaxies in ...
A manifestation of this co-evolution is the correlation that has been found between the SMBH mass, MBH, and the galaxy bulge or stellar mass, M*. The cosmic ...
[75]
The evolution of supermassive black hole mass–bulge mass relation ...
May 13, 2024 · In this study, we focus on the evolution with redshift of the relationship between the mass of an SMBH and the mass of its host bulge, in order ...
[76]
[2309.06475] The thesan project: public data release of radiation ...
Sep 12, 2023 · This project focuses on the z \geq 5.5 Universe, combining a radiation-hydrodynamics solver (AREPO-RT), a well-tested galaxy formation model (IllustrisTNG) and ...
[77]
[1811.11192] Cosmic Dawn II (CoDa II) - arXiv
Nov 27, 2018 · CoDa II is a new, fully-coupled radiation-hydrodynamics simulation of cosmic reionization and galaxy formation and their mutual impact, to redshift z < 6.
[78]
[2503.19131] Turbulent dynamos in a collapsing cloud - arXiv
Mar 24, 2025 · The amplification of magnetic fields is crucial for understanding the observed magnetization of stars and galaxies. Turbulent dynamo is the ...
[79]
[2007.12837] Time-resolved fast turbulent dynamo in a laser plasma
Jul 25, 2020 · Understanding magnetic-field generation and amplification in turbulent plasma is essential to account for observations of magnetic fields in the ...<|control11|><|separator|>
[80]
Sphenix: Smoothed Particle Hydrodynamics for the next generation ...
Dec 7, 2020 · Sphenix uses a Density-Energy equation of motion, along with variable artificial viscosity and conduction, including limiters designed to work ...