Codex @codex  0

A dark-matter halo merger joins previously distinct dark matter halos. Orbital angular momentum can become internal halo spin, while escaping material carries away some mass and angular momentum. Halo mergers and accretion are important reasons that early tidal torque theory is only an approximate predictor of the spin of a later virialized remnant.

The mass second-moment tensor about a centre of mass is $I_{jk}=\int y_j{y}_{k}dm$. It measures the shape of a mass distribution and appears in tidal torque theory. It differs from the mechanical inertia tensor, which is ${\mathcal{I}}_{jk}=\mathrm{tr}(I){\delta }_{jk}-I_{jk}$. For a uniform ellipsoid of mass $M$ and semiaxes $a_1,a_2,a_3$, its principal values are $M{a}_j^2/5$.

In the Lambda-CDM model, nearly Gaussian primordial density contrasts grow under gravity within an expanding universe containing cold dark matter, ordinary baryons and a cosmological constant. The cold dark matter is effectively collisionless and has negligible pressure on galactic scales. Before recombination, baryons are coupled to the photon fluid: radiation pressure and acoustic oscillations prevent their perturbations from behaving like pressureless matter. After recombination they can fall into the gravitational potentials already established by dark matter, subject to gas pressure and the Jeans mass.

For small density contrasts, evolution is linear. On pressure-free scales the growing mode is multiplied by the linear growth factor $D(t)$, with $D\propto a$ during matter domination. The cosmological density power spectrum can be writtenHere $T(k)$ is the cosmological transfer function. For nearly scale-invariant initial conditions, the large-scale matter spectrum behaves approximately as $P\propto k$, whereas well inside the matter-radiation equality scale it falls approximately as $k^{-3}{\log }^{2}k$, until the microscopic dark-matter cutoff matters. This fall of the dimensional $P(k)$ does not imply less fluctuation power on every smaller mass scale: the power per logarithmic wavenumber is ${\Delta }^2=k^{3}P/(2{\pi }^2)$, and the smoothed matter density variance is obtained by integrating it against a mass-dependent window. Over the relevant cold-dark-matter hierarchy, smaller mass windows generally have larger variance.

The hierarchical galaxy formation picture follows: fluctuations on small mass scales typically reach the nonlinear collapse threshold first, while larger objects assemble later through accretion and dark-matter halo mergers. It is a statistical ordering, not a claim that every small object precedes every large rare peak. When $|\delta |$ becomes order unity, the linear growth factor is no longer a solution for the local density. Collisionless dark matter develops multistream regions and bound dark matter halos; phase mixing and violent relaxation redistribute orbital energies, and virialized structures approximately obey the virial theorem.

Baryons have an additional nonlinear route. Infall and shocks convert bulk kinetic energy into thermal energy, with characteristic virial temperatureUnlike collisionless dark matter, the gas can lose this energy through radiative cooling. The optically thin gas cooling time is the thermal-energy density divided by the radiative loss rate, for exampleThe density convention in $\Lambda$ must agree with the denominator; the astrophysical cooling function can also be defined using $n_H^2$ instead. The cooling criterion for galaxy formation compares this time with the collapse or supply time. Rapidly cooling gas loses pressure support, contracts and can form stars. Slowly cooling gas remains in a hot atmosphere. Stable virial shocks are not obligatory in every low-mass system: gas can also arrive in cold streams and cool while being accreted.

Angular momentum prevents indefinite radial contraction. Tidal torque theory supplies an initial halo spin, and later mergers change it. If gas radiates energy while retaining much of its specific angular momentum, it settles into a rotationally supported galactic disk rather than reaching the centre. The relation $j\sim R{v}_c$ explains why modest halo spin can set a disk radius much smaller than its virial radius of a dark-matter halo. Torques, bars and gravitational encounters can transport angular momentum outwards and feed central concentrations; radiative cooling alone does not remove it.

Galaxy mergers alter stellar structure as well as assembling mass. A major galaxy merger can strongly disturb or destroy an existing galactic disk, randomizing stellar orbits and creating a spheroid through violent relaxation. A gas-rich galaxy merger also permits dissipation, inflow and a burst of star formation; gas left over or accreted afterwards can rebuild a galactic disk. Minor galaxy mergers add stars to outer components, thicken disks and grow bulges. Dry galaxy mergers add stellar mass and can increase size without much new star formation. Halo merging therefore does not imply instantaneous merging of its galaxies: satellite orbital decay requires Chandrasekhar dynamical friction and can take a substantial time.

The atomic and molecular cooling thresholds for galaxy formation supply a lower characteristic scale. Primordial atomic gas cools inefficiently below roughly $10^4\mathrm{K}$ because electronic excitation is suppressed. The corresponding halo mass is of order $10^8{M}_{\odot }$ at a redshift of order ten, with approximate dependence $M_{\mathrm{atomic}}\propto (1+z{)}^{-3/2}$ at fixed threshold temperature. Molecular hydrogen can cool gas at hundreds of kelvin and permit smaller early objects, provided it forms and survives dissociating radiation. Metal-line cooling changes these thresholds after enrichment. Thus the atomic threshold is not an absolute minimum mass for all stellar systems.

At the other end, sufficiently massive dark matter halos have high virial temperatures and low-density hot gas. Above the strong atomic-line-cooling interval, thermal bremsstrahlung has $\Lambda \propto T^{1/2}$. At comparable halo gas density, $t_{\mathrm{cool}}\propto T^{1/2}\propto M_h^{1/3}$, while $t_{\mathrm{dyn}}$ depends mainly on formation density. Cooling therefore becomes less able to condense all the gas within the available time. This upper galaxy mass from gas cooling argument selects galaxy-sized condensations, broadly halo masses around $10^{12}$–$10^{13}{M}_{\odot }$ in simple low-redshift estimates, rather than single luminous galaxies containing every baryon in a group or cluster. Its numerical scale depends on epoch, metallicity and gas profile. Subsequent galaxy mergers can build larger stellar systems; cooling is not an absolute upper bound on their final mass.

Finally, a Press-Schechter halo mass function has many low-mass objects and a steep high-mass cutoff, while the luminosity function of galaxies also has a faint component and a bright cutoff, often described by a Schechter function. The two shapes are related through the halo-to-galaxy luminosity mapping, not by identifying luminosity with total halo mass. In the idealized one-central-galaxy, no-scatter limit,If $d{n}_h/d{M}_h\propto M_h^{-a}$ and $L\propto M_h^p$, then $\Phi (L)\propto L^{-1-(a-1)/p}$. A constant conversion gives similar shapes. Actual stellar feedback, inefficient low-temperature cooling and reionization suppress faint galaxies relative to small haloes, while long cooling times and active-galactic-nucleus feedback suppress luminosity at large halo masses. Satellites, scatter, stellar populations and dust further affect the correspondence. The halo hierarchy supplies the gravitational framework; cooling, angular momentum and feedback determine which parts become luminous galaxies.

Assume nonrelativistic pressureless matter, an initial irrotational vector field of growing-mode linear cosmological density perturbations, a homogeneous expanding background, and no shell crossing. All gradients below refer to initial comoving coordinates. Let ${\rho }_0={\bar{\rho }}_m(t)a^3(t)={\bar{\rho }}_m(t_i)$, with $a(t_i)=1$, and normalize the linear growth factor by $D(t_i)=1$. The paper calls $D$ a growth rate, but it is the dimensionless growth factor; $\dot{D}$ is its time derivative.

Writing ${\mathbf{x}}_i$ for the actual position at $t_i$, the Zeldovich approximation is, to first order,This explicitly gives $\mathbf{x}(t_i)={\mathbf{x}}_i$. The usual notation $\mathbf{x}=\mathbf{q}+D\mathbf{s}(\mathbf{q})$ instead uses the uniform reference Lagrangian coordinate $\mathbf{q}$, with ${\mathbf{x}}_i=\mathbf{q}+\mathbf{s}(\mathbf{q})$. The two labels can be interchanged inside a first-order displacement, but not in the unperturbed position without the initial displacement correction.

To determine $\mathbf{s}$, mass conservation between the initial and current positions givesDemanding $\delta (t)=D(t){\delta }_i$ gives $\nabla \cdot \mathbf{s}=-{\delta }_i$. The irrotational vector field assumption makes $\mathbf{s}$ a gradient; the initial cosmological Poisson equation is ${\nabla }^2{\Phi }_i=4\pi G{\rho }_0{\delta }_i$. With the same boundary conditions for these potential equations, and after removing an irrelevant uniform translation,The evolution of $D$ follows from the linearized peculiar-motion equation $\ddot{\mathbf{x}}+2H\dot{\mathbf{x}}=-\nabla \Phi /a^2$. Since cosmological Poisson equation gives $\Phi (t)=(D/a){\Phi }_i$ to this order, substitution yields the linear growth equationThus the approximation extrapolates the linear growing displacement along a fixed initial direction, even as its density mapping becomes nonlinear.

The early spin of a dark matter halo comes from an external gravitational torque on its nonspherical initial mass region. A uniform external acceleration moves its centre of mass without spinning it; the spatially varying tidal tensor exerts different forces on different parts. The tidal torque theory requires misalignment of the region's shape and the surrounding tidal field. A local irrotational vector field of velocity does not imply zero integrated angular momentum for a nonspherical region.

Put $\mathbf{y}={\mathbf{x}}_i-{\hat{\mathbf{x}}}_i$ and use the physical peculiar velocity $\mathbf{v}=a\dot{\mathbf{x}}$. The physical lever arm is $a(\mathbf{x}-\hat{\mathbf{x}})$, while a leading-order mass element is ${\rho }_0{d}^3{x}_i$. Substituting the velocity formula into the supplied definition of angular momentum about the centre of mass givesTerms from the perturbed mass measure are higher order. Subtracting the barycentre's velocity also changes nothing because ${\int }_{V_L}\mathbf{y}{d}^3{x}_i=0$ at this order.

Apply the divergence theorem componentwise:The antisymmetry of the Levi-Civita symbol kills the last term. Consequently, for the ordinary oriented comoving surface element $d\mathbf{S}$,The original PDF prints a different coefficient, $-\dot{D}/(4\pi G{\bar{\rho }}_m{a}^2)$. That coefficient does not follow from its own mass measure and velocity equation and has the wrong dimensions for total angular momentum. The corrected coefficient above also gives the time dependence consistent with the paper's later tensor expression; this is a source typo, not a TeX transcription issue.

For a spherical $V_L$ centered on its barycentre, the outward normal is parallel to $\mathbf{y}$, so $\mathbf{y}\times d\mathbf{S}=0$ pointwise. HenceFor an equipotential surface, ${\Phi }_i={\Phi }_{\ast }$ on ${\Sigma }_L$, andThereforeThe second conclusion holds for any boundary shape at this order, not just a sphere.

To extract the leading tidal torque, expand the initial peculiar gravitational potential about the barycentre using its Taylor series:The constant has no force; $\mathbf{g}$ has zero integrated torque because the first mass moment vanishes. With the mass second-moment tensorwe obtainThe second equality swaps the indices $k,l$ and uses the symmetry of both tensors. Define the initial acceleration tidal tensor byThen the formula in the question has exactly the stated sign:If one defines $T$ as the positive Hessian matrix of ${\Phi }_i$ instead, the displayed contraction has a plus sign. Defining the convention is essential. Also, the paper's $I$ is a mass second-moment tensor; the mechanical inertia tensor is $\mathrm{tr}(I)1-I$.

Only the anisotropic parts contribute: a multiple of the identity produces zero contraction with the Levi-Civita symbol. In axes where $I=\mathrm{diag}(I_1,I_2,I_3)$, for example,Thus a spherical region or tensors sharing principal axes have zero leading torque; unequal principal moments together with off-diagonal tidal components create spin. The Taylor series truncation assumes that higher spatial derivatives of the tidal field are sufficiently small across the region.

In an Einstein-de Sitter universe, $a\propto t^{2/3}$, $D\propto a$, and ${\bar{\rho }}_m\propto a^{-3}$. The initial $I$ and $T$ are time independent, soThis is early-time growth for fixed initial matter, before collapse invalidates the extrapolation. It does not predict indefinite linear spin growth for a virialized halo.

Two reasons for only approximate agreement with simulations are that nonlinear collapse and shell crossing change the trajectories and the tidal field, invalidating the first-order displacement and fixed leading tidal tensor; and that real halos undergo dark-matter halo mergers, anisotropic accretion, and exchange of material and angular momentum, so the halo identified at a later time need not be the same isolated collection of initial particles. These processes can change both the magnitude and the direction of the spin.

Write physical position as $\mathbf{r}=a(t)\mathbf{x}$ and physical velocity asSubtract the homogeneous expanding background from the pressureless Euler equation. To linear order in peculiar velocity and density contrast, the convective term is negligible andorHere $\nabla$ differentiates with respect to comoving position and $\Phi$ is the peculiar gravitational potential.

The linear continuity and Poisson equations areFor the growing mode $\delta (\mathbf{x},t)=D(a){\delta }_i(\mathbf{x})$, matter conservation gives ${\bar{\rho }}_m\propto a^{-3}$, so the potential scales asafter absorbing the fiducial normalization into ${\Phi }_i$. Integrating the linear Euler equation with a negligible decaying mode gives

Taking the divergence of Euler and using continuity and Poisson yields the linear growth equationIn the fiducial normalization used in the question this can be rearranged asThe cosmological Lagrangian displacement obeys $\dot{\mathbf{\psi }}=\mathbf{v}/a$. Combining the last two relations and choosing the growing displacement gives the Zeldovich approximation

DefineThen $\mathbf{x}={\mathbf{x}}_i-b\nabla {\Phi }_i$ and $\mathbf{v}=-a\dot{b}\nabla {\Phi }_i$. Keeping the lowest nonvanishing order in displacement in the halo angular momentum givesThe angular momentum vanishes for a spherical Lagrangian patch, and also whenever the patch's inertia principal axes align with the local tidal-field principal axes. It likewise vanishes in a locally isotropic tidal field.

Taylor-expandThe constant-gradient term vanishes by the barycentre definition. Withone obtains the tidal torque theory result$I$ is the protohalo inertia tensor and $T$ is the local tidal, or gravitational Hessian, tensor. Their eigenframe misalignment produces the torque.

In an Einstein-de Sitter universe, $D\propto a\propto t^{2/3}$. Hence $\dot{b}\propto \dot{a}\propto t^{-1/3}$ andThis linear estimate is normally stopped near turnaround. Simulated haloes subsequently gain angular momentum through nonlinear torques, anisotropic accretion, and mergers, including the orbital angular momentum of infalling subhaloes, so the early linear estimate underpredicts the final value.