Spiral Structure Properties, Dynamics, and Evolution in MW-mass Galaxy Simulations

The structure of spiral galaxies is essential to understanding the dynamics and evolution of disk galaxies; however, the precise nature of spiral arms remains uncertain. Two challenges in understanding the mechanisms driving spirals are how galactic environment impacts spiral morphology and how they evolve over time. We present a catalog characterizing the properties, dynamics, and evolution of $m=2$ spiral structure in 10 Milky Way-mass galaxies from the FIRE-2 cosmological zoom-in simulations. Consistent with previous literature, we find that FIRE-2 spirals are transient, recurring features simultaneously present in the disk at varying pattern speeds ($Ω_p$) that broadly decrease with radius. These spirals persist on Gyr timescales (mean duration 1.90 Gyr), but fluctuate in amplitude on timescales of hundreds of Myr. Tidal interactions and bar episodes impact the resulting $m=2$ spiral structure; strong satellite interactions generally produce shorter-lived, stronger spirals with larger radial extent, and bars can increase $Ω_p$. Galactic environment influences spiral structure; kinematically colder disks can support longer-lived, stronger spirals. The properties of identified spirals in FIRE-2 vary widely in radial extent (0.3-10.8 kpc), duration (1.00-6.00 Gyr), and amplitudes ($a_{2,\text{max}}$=0.018-0.192). We find the presence of spirals in all age populations, suggesting these are density wave-driven features. This work represents the first time that spiral structure has been cataloged in this manner in cosmological simulations; the catalog can be leveraged with current and forthcoming observational surveys, enabling systematic comparisons to further our understanding of galaxy evolution.

💡 Research Summary

This paper presents a comprehensive catalog of m = 2 spiral structure in ten Milky‑Way‑mass galaxies simulated with the FIRE‑2 cosmological zoom‑in framework. The authors aim to quantify the properties, dynamics, and temporal evolution of spiral arms and to explore how galactic environment and evolutionary history shape these features.
Simulation suite and methodology – The study uses five isolated galaxies from the Latte suite and five Local‑Group‑like pairs from the Elvis on FIRE suite, each resolved with baryonic particle masses of ~5 × 10³ M⊙ and adaptive gas softening down to ~1 pc. Snapshots are saved every 25 Myr (≈2 Myr in the last 100 Myr) from 6.8 Gyr ago to the present. At each snapshot the stellar component within 10 kpc is re‑oriented so that the angular momentum vector aligns with the z‑axis, ensuring a consistent disk frame. The authors then apply a windowed discrete Fourier transform (WDFT) to the 2‑D stellar mass distribution, binning stars in 0.3 kpc annuli without any vertical cut. This yields time‑resolved amplitudes a₂(R,t) and pattern speeds Ωₚ(R,t) for the m = 2 mode.
Key findings –

1. Transient, recurring spirals: All ten galaxies satisfy a stellar Toomre Q > 1 beyond the central few kiloparsecs, allowing self‑gravitating spiral modes to appear repeatedly. The spirals are transient, with individual episodes lasting 1–6 Gyr (mean ≈ 1.9 Gyr) but with amplitude fluctuations on ∼100–300 Myr timescales.
2. Pattern speed gradients: Ωₚ generally declines with radius, consistent with a winding, density‑wave picture. In galaxies that develop a central bar, the overall Ωₚ is higher, indicating bar‑driven reinforcement of the spiral pattern.
3. Environmental impact: Strong satellite encounters (e.g., in m12f) generate spirals that extend to larger radii (up to 10.8 kpc) and have higher peak amplitudes (a₂,max ≈ 0.19) but are shorter‑lived (≈1 Gyr). Conversely, kinematically colder disks—characterized by low radial velocity dispersion and high stellar surface density—host longer‑lived, stronger spirals.
4. Age‑independent presence: The catalog detects the same m = 2 pattern in young (<1 Gyr), intermediate (∼3 Gyr), and old (>5 Gyr) stellar populations, supporting the interpretation that the arms are genuine density waves rather than transient star‑formation clumps.
5. Amplitude distribution: The maximum amplitudes span 0.018–0.192, covering the range observed in external galaxies and indicating a broad diversity of spiral strengths even within a narrow mass range.
   Implications and future work – By providing a publicly available, time‑resolved spiral catalog, the authors enable direct comparisons with upcoming large‑scale surveys such as Gaia DR3/DR4, SDSS‑V, the Roman Space Telescope, and JWST. Matching observed pattern speeds, radial extents, and lifetimes to the simulated catalog will allow the community to test competing theories (classical density‑wave, swing amplification, bar‑driven, tidal triggering, and self‑excited modes) in a fully cosmological context. Moreover, the study demonstrates that high‑resolution, physics‑rich simulations like FIRE‑2 can capture the nuanced interplay between internal dynamics (bars, Q‑stability) and external perturbations (satellites) that governs spiral morphology. This work thus establishes a benchmark for future investigations into the physical drivers of spiral structure and their role in redistributing angular momentum, fueling star formation, and shaping the secular evolution of disk galaxies.