Cold Dark Matter Substructure and Galactic Disks II: Dynamical Effects of Hierarchical Satellite Accretion

(Abridged) We perform dissipationless N-body simulations to elucidate the dynamical response of thin disks to bombardment by cold dark matter (CDM) substructure. Our method combines (1) cosmological simulations of the formation of Milky Way (MW)-sized CDM halos to derive the properties of substructure and (2) controlled numerical experiments of consecutive subhalo impacts onto an initially-thin, fully-formed MW type disk galaxy. The present study is the first to account for the evolution of satellite populations over cosmic time in such an investigation of disk structure. We find that accretions of massive subhalos onto the central regions of host halos, where the galactic disks reside, since z~1 should be common. One host halo accretion history is used to initialize the controlled simulations of satellite-disk encounters. We show that these accretion events severely perturb the thin galactic disk and produce a wealth of distinctive dynamical signatures on its structure and kinematics. These include (1) considerable thickening and heating at all radii, with the disk thickness and velocity ellipsoid nearly doubling at the solar radius; (2) prominent flaring associated with an increase in disk thickness greater than a factor of 4 in the disk outskirts; (3) surface density excesses at large radii, beyond ~5 disk scale lengths, resembling those of observed antitruncated disks; (4) lopsidedness at levels similar to those measured in observational samples of disk galaxies; and (5) substantial tilting. The interaction with the most massive subhalo drives the disk response while subsequent bombardment is much less efficient at disturbing the disk. We conclude that substructure-disk encounters of the kind expected in the LCDM paradigm play a significant role in setting the structure of disk galaxies and driving galaxy evolution.

💡 Research Summary

This paper investigates how cold‑dark‑matter (CDM) substructure perturbs the dynamical state of thin galactic disks by means of high‑resolution, dissipationless N‑body simulations. The authors adopt a two‑step methodology. First, they run cosmological simulations of Milky Way‑sized dark‑matter halos to extract realistic subhalo populations from redshift z ≈ 1 to the present. The resulting subhalos span masses of 10⁹–10¹⁰ M⊙, have a broad distribution of orbital eccentricities, and frequently pass within the inner ~30 kpc of the host, in agreement with ΛCDM predictions. Second, they construct a fully formed, equilibrium stellar disk (scale length ≈ 3 kpc, vertical scale height ≈ 300 pc) and subject it to a sequence of controlled encounters using the exact masses, positions, and velocities of the cosmologically derived subhalos. Five impacts are modeled, each separated by ~1 Gyr, allowing the disk to relax between events. No gas dynamics or star formation is included; the study isolates pure gravitational effects.

The simulations reveal a suite of robust, observable signatures. (1) Thickening and heating: The vertical scale height roughly doubles across the disk, and the vertical velocity dispersion at the solar radius (≈ 8 kpc) rises from ~15 km s⁻¹ to ~30 km s⁻¹. (2) Flaring: In the outer disk (> 3 scale lengths) the thickness increases by a factor of four or more, reproducing the flared profiles seen in many edge‑on galaxies. (3) Antitruncated surface‑density profiles: Beyond ~5 scale lengths the originally exponential decline flattens, creating an excess of stars at large radii that resembles observed Type III (antitruncated) disks. (4) Lopsidedness: The m = 1 Fourier mode of the surface density reaches amplitudes of ~0.1, comparable to the asymmetries measured in large samples of late‑type spirals. (5) Disk tilting: The entire stellar disk is reoriented by 8–12°, reflecting angular‑momentum transfer from the infalling subhalos.

A key insight is that the most massive subhalo (≈ 1 × 10¹⁰ M⊙) dominates the dynamical response; its first pericentric passage drives the bulk of the thickening, flaring, and tilting. Subsequent impacts by lower‑mass subhalos produce comparatively modest additional heating, indicating that once the disk has been pre‑heated and puffed up, it becomes more resilient to further perturbations. The authors argue that these results naturally explain several longstanding observational puzzles: the prevalence of thick disks, the existence of flared outer disks, the frequency of antitruncated profiles, and the commonality of mild lopsidedness in isolated spirals.

By explicitly incorporating the time‑dependent subhalo accretion history, the study improves upon earlier work that assumed a single, ad‑hoc satellite encounter. It demonstrates that the hierarchical assembly of CDM substructure, as expected in the standard cosmological model, can by itself generate the bulk of the structural and kinematic features observed in present‑day disk galaxies, even in the absence of gas physics. The authors conclude that satellite‑disk interactions are a fundamental driver of disk galaxy evolution within ΛCDM, and that future work should explore how gas dissipation, star formation, and feedback may modulate or amplify these gravitational effects.