Efficient and accurate force replay in cosmological-baryonic simulations

Abstract

We construct time-evolving gravitational potential models for a Milky Way-mass galaxy from the FIRE-2 suite of cosmological-baryonic simulations using basis function expansions. These models capture the angular variation with spherical harmonics for the halo and azimuthal harmonics for the disk, and the radial or meridional plane variation with splines. We fit low-order expansions (4 angular/harmonic terms) to the galaxy's potential for each snapshot, spaced roughly 25 Myr apart, over the last 4 Gyr of its evolution, then extract the forces at discrete times and interpolate them between adjacent snapshots for forward orbit integration. Our method reconstructs the forces felt by simulation particles with high fidelity, with 95% of both stars and dark matter, outside of self-gravitating subhalos, exhibiting errors $\le$4% in both the disk and the halo. Imposing symmetry on the model systematically increases these errors, particularly for disk particles, which show greater sensitivity to imposed symmetries. The majority of orbits recovered using the models exhibit positional errors $\le$10% for 2-3 orbital periods, with higher errors for orbits that spend more time near the galactic center. Approximate integrals of motion are retrieved with high accuracy even with a larger potential sampling interval of 200 Myr. After 4 Gyr of integration, 43% and 70% of orbits have total energy and angular momentum errors within 10%, respectively. Consequently, there is higher reliability in orbital shape parameters such as pericenters and apocenters, with errors $\sim$10% even after multiple orbital periods. These techniques have diverse applications, including studying satellite disruption in cosmological contexts.

29 pages, 13 figures, and 4 tables. Published in APJ

References in corpus (68)

• Array Programming with NumPy
• Planck 2015 results. XIII. Cosmological parameters
• galpy: A Python Library for Galactic Dynamics
• The Rockstar Phase-Space Temporal Halo Finder and the Velocity Offsets of Cluster Cores
• Concentrations of Dark Halos from their Assembly Histories
• FIRE-2 Simulations: Physics versus Numerics in Galaxy Formation
• The mass distribution and gravitational potential of the Milky Way
• The Inner Structure of LambdaCDM Halos I: A Numerical Convergence Study
• Reconciling dwarf galaxies with LCDM cosmology: Simulating a realistic population of satellites around a Milky Way-mass galaxy
• Formation and evolution of galaxy dark matter halos and their substructure
• Are the Magellanic Clouds on their First Passage about the Milky Way?
• AGAMA: Action-based galaxy modelling architecture
• The Sagittarius Dwarf Galaxy: a Model for Evolution in a Triaxial Milky Way Halo
• A direct dynamical measurement of the Milky Way's disk surface density profile, disk scale length, and dark matter profile at 4 kpc < R < 9 kpc
• Origin of the system of globular clusters in the Milky Way
• Constraining the Milky Way potential with a 6-D phase-space map of the GD-1 stellar stream
• The total mass of the Large Magellanic Cloud from its perturbation on the Orphan stream
• ELVIS: Exploring the Local Volume in Simulations
• CMasher: Scientific colormaps for making accessible, informative and 'cmashing' plots
• Tango for three: Sagittarius, LMC, and the Milky Way
• Proper motions and dynamics of the Milky Way globular cluster system from Gaia DR2
• Gaia DR2 Proper Motions of Dwarf Galaxies within 420 kpc: Orbits, Milky Way Mass, Tidal Influences, Planar Alignments, and Group Infall
• The Influence of Sagittarius and the Large Magellanic Cloud on the Milky Way Galaxy
• And yet it moves: The dangers of artificially fixing the Milky Way center of mass in the presence of a massive Large Magellanic Cloud
• A dark matter profile to model diverse feedback-induced core sizes of $\Lambda$CDM haloes
• A class of symplectic integrators with adaptive timestep for separable Hamiltonian systems
• Gaia Proper Motions and Orbits of the Ultra-Faint Milky Way Satellites
• Constraints on the shape of the Milky Way dark matter halo from the Sagittarius stream
• Public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation
• The Planes of Satellite Galaxies Problem, Suggested Solutions, and Open Questions
• A profile in FIRE: resolving the radial distributions of satellite galaxies in the Local Group with simulations
• Constraining the Milky Way Halo Potential with the GD-1 stellar stream
• Orbits of Massive Satellite Galaxies: I. A Close Look at the Large Magellanic Cloud and a New Orbital History for M33
• Quantifying the impact of the Large Magellanic Cloud on the structure of the Milky Way's dark matter halo using Basis Function Expansions
• Proper Motions, Orbits, and Tidal Influences of Milky Way Dwarf Spheroidal Galaxies
• Measuring the Mass of the Large Magellanic Cloud with Stellar Streams Observed by $S^5$
• $S^5$: Probing the Milky Way and Magellanic Clouds potentials with the 6-D map of the Orphan-Chenab stream
• The Orbit and Mass of the Sagittarius Dwarf Galaxy
• Milky Way Mass and Potential Recovery Using Tidal Streams in a Realistic Halo
• Signatures of resonance and phase mixing in the Galactic disc
• Determining orbits for the Milky Way's dwarfs
• Reflex motion in the Milky Way stellar halo resulting from the Large Magellanic Cloud infall
• A Milky Way with a massive, centrally concentrated thick disc: new Galactic mass models for orbit computations
• Modelling the Tucana III stream - a close passage with the LMC
• The Orbital Ellipticity of Satellite Galaxies and the Mass of the Milky Way
• The effect of the deforming dark matter haloes of the Milky Way and the Large Magellanic Cloud on the Orphan-Chenab stream
• Quantifying the Stellar Halo's Response to the LMC's Infall with Spherical Harmonics
• The physical nature of the cosmic accretion of baryons and dark matter into halos and their galaxies
• Dark matter halo shapes in the Auriga simulations
• A new code for orbit analysis and Schwarzschild modelling of triaxial stellar systems
• Estimating the Mass of the Milky Way Using the Ensemble of Classical Satellite Galaxies
• A Halo Expansion (HEX) Technique for Approximating Simulated Dark Matter Haloes
• Streams on FIRE: Populations of Detectable Stellar Streams in the Milky Way and FIRE
• Uncertainties associated with the backward integration of dwarf satellites using simple parametric potentials
• The infall of dwarf satellite galaxies are influenced by their host's massive accretions
• The Milky Way's stellar streams and globular clusters do not align in a Vast Polar Structure
• Simulating tidal streams in a high resolution dark matter halo
• On the stability of tidal streams in action space
• EXP: N-body integration using basis function expansions
• DM halo morphological types of MW-like galaxies in the TNG50 simulation: Simple, Twisted, or Stretched
• Orbital dynamics and histories of satellite galaxies around Milky Way-mass galaxies in the FIRE simulations
• The observable properties of galaxy accretion events in Milky Way-like galaxies in the FIRE-2 cosmological simulations
• Milky Way globular clusters on cosmological timescales. I. Evolution of the orbital parameters in time-varying potentials
• Models of Distorted and Evolving Dark Matter Halos
• LMC-driven anisotropic boosts in stream--subhalo interactions
• Basis function expansions for galactic dynamics: Spherical versus cylindrical coordinates
• Bar formation and destruction in the FIRE-2 simulations
• Planes of satellites, no longer in tension with LCDM

Cited by in corpus (4)

• Shaping the Milky Way: The interplay of mergers and cosmic filaments
• Dynamics of tidal spiral arms: Machine learning-assisted identification of equations and application to the Milky Way
• Variable gravitational potential of Milky Way analogues in HESTIA suite
• No Stream Left Unscathed: The imprint of a host galaxy