An official website of the United States government
ABSTRACT Understanding the evolution of satellite galaxies of the Milky Way (MW) and M31 requires modelling their orbital histories across cosmic time. Many works that model satellite orbits incorrectly assume or approximate that the host halo gravitational potential is fixed in time and is spherically symmetric or axisymmetric. We rigorously benchmark the accuracy of such models against the FIRE-2 cosmological baryonic simulations of MW/M31-mass haloes. When a typical surviving satellite fell in ($$3.4\!-\!9.7\, \rm {Gyr}$$ ago), the host halo mass and radius were typically 26–86 per cent of their values today, respectively. Most of this mass growth of the host occurred at small distances, $$r\lesssim 50\, \rm {kpc}$$, opposite to dark matter only simulations, which experience almost no growth at small radii. We fit a near-exact axisymmetric gravitational potential to each host at z = 0 and backward integrate the orbits of satellites in this static potential, comparing against the true orbit histories in the simulations. Orbital energy and angular momentum are not well conserved throughout an orbital history, varying by 25 per cent from their current values already $$1.6\!-\!4.7\, \rm {Gyr}$$ ago. Most orbital properties are minimally biased, ≲10 per cent, when averaged across the satellite population as a whole. However, for a single satellite, the uncertainties are large: recent orbital properties, like the most recent pericentre distance, typically are ≈20 per cent uncertain, while earlier events, like the minimum pericentre or the infall time, are ≈40–80 per cent uncertain. Furthermore, these biases and uncertainties are lower limits, given that we use near-exact host mass profiles at z = 0.
more »
« less
This content will become publicly available on November 29, 2025
Efficient and Accurate Force Replay in Cosmological-baryonic Simulations
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 (four 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 ≤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 ≤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 ∼10% even after multiple orbital periods. These techniques have diverse applications, including studying satellite disruption in cosmological contexts.
more »
« less
- Award ID(s):
- 2045928
- PAR ID:
- 10627118
- Publisher / Repository:
- The Astrophysical Journal
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 977
- Issue:
- 1
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 23
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT The orbits of satellite galaxies encode rich information about their histories. We investigate the orbital dynamics and histories of satellite galaxies around Milky Way (MW)-mass host galaxies using the FIRE-2 cosmological simulations, which, as previous works have shown, produce satellite mass functions and spatial distributions that broadly agree with observations. We first examine trends in orbital dynamics at z = 0, including total velocity, specific angular momentum, and specific total energy: the time of infall into the MW-mass halo primarily determines these orbital properties. We then examine orbital histories, focusing on the lookback time of first infall into a host halo and pericentre distances, times, and counts. Roughly 37 per cent of galaxies with $$M_{\rm star}\lesssim 10^7\, {\rm M}_{\odot }$$ were ‘pre-processed’ as a satellite in a lower-mass group, typically $$\approx 2.7\, {\rm Gyr}$$ before falling into the MW-mass halo. Half of all satellites at z = 0 experienced multiple pericentres about their MW-mass host. Remarkably, for most (67 per cent) of these satellites, their most recent pericentre was not their minimum pericentre: the minimum typically was ∼40 per cent smaller and occurred $$\sim 6\, {\rm Gyr}$$ earlier. These satellites with growing pericentres appear to have multiple origins: for about half, their specific angular momentum gradually increased over time, while for the other half, most rapidly increased near their first apocentre, suggesting that a combination of a time-dependent MW-mass halo potential and dynamical perturbations in the outer halo caused these satellites’ pericentres to grow. Our results highlight the limitations of idealized, static orbit modelling, especially for pericentre histories.more » « less
-
Abstract Recent observations of the stellar halo have uncovered the debris of an ancient merger, Gaia–Sausage–Enceladus (GSE), estimated to have occurred ≳8 Gyr ago. Follow-up studies have associated GSE with a large-scale tilt in the stellar halo that links two well-known stellar overdensities in diagonally opposing octants of the Galaxy (the Hercules–Aquila Cloud and Virgo Overdensity; HAC and VOD). In this paper, we study the plausibility of such unmixed merger debris persisting over several gigayears in the Galactic halo. We employ the simulated stellar halo from Naidu et al., which reproduces several key properties of the merger remnant, including the large-scale tilt. By integrating the orbits of these simulated stellar halo particles, we show that adoption of a spherical halo potential results in rapid phase mixing of the asymmetry. However, adopting a tilted halo potential preserves the initial asymmetry in the stellar halo for many gigayears. The asymmetry is preserved even when a realistic growing disk is added to the potential. These results suggest that HAC and VOD are long-lived structures that are associated with GSE and that the dark matter halo of the Galaxy is tilted with respect to the disk and aligned in the direction of HAC–VOD. Such halo–disk misalignment is common in modern cosmological simulations. Lastly, we study the relationship between the local and global stellar halo in light of a tilted global halo comprised of highly radial orbits. We find that the local halo offers a dynamically biased view of the global halo due to its displacement from the Galactic center.more » « less
-
Abstract A significant fraction of Milky Way (MW) satellites exhibit phase-space properties consistent with a coherent orbital plane. Using tailored N -body simulations of a spherical MW halo that recently captured a massive (1.8 × 10 11 M ⊙ ) LMC-like satellite, we identify the physical mechanisms that may enhance the clustering of orbital poles of objects orbiting the MW. The LMC deviates the orbital poles of MW dark matter particles from the present-day random distribution. Instead, the orbital poles of particles beyond R ≈ 50 kpc cluster near the present-day orbital pole of the LMC along a sinusoidal pattern across the sky. The density of orbital poles is enhanced near the LMC by a factor δ ρ max = 30% (50%) with respect to underdense regions and δ ρ iso = 15% (30%) relative to the isolated MW simulation (no LMC) between 50 and 150 kpc (150–300 kpc). The clustering appears after the LMC’s pericenter (≈50 Myr ago, 49 kpc) and lasts for at least 1 Gyr. Clustering occurs because of three effects: (1) the LMC shifts the velocity and position of the central density of the MW’s halo and disk; (2) the dark matter dynamical friction wake and collective response induced by the LMC change the kinematics of particles; (3) observations of particles selected within spatial planes suffer from a bias, such that measuring orbital poles in a great circle in the sky enhances the probability of their orbital poles being clustered. This scenario should be ubiquitous in hosts that recently captured a massive satellite (at least ≈1:10 mass ratio), causing the clustering of orbital poles of halo tracers.more » « less
-
Abstract There is an intriguing and growing population of Neptune-sized planets with stellar obliquities near ∼90°. One previously proposed formation pathway is a disk-driven resonance, which can take place at the end stages of planet formation in a system containing an inner Neptune, outer cold Jupiter, and protoplanetary disk. This mechanism occurs within the first ∼10 Myr, but most of the polar Neptunes we see today are ∼Gyr old. Up until now, there has not been an extensive analysis of whether the polar orbits are stable over ∼Gyr timescales. Tidal realignment mechanisms are known to operate in other systems, and if they are active here, this would cause theoretical tension with a primordial misalignment story. In this paper, we explore the effects of tidal evolution on the disk-driven resonance theory. We use bothN-body and secular simulations to study tidal effects on both the initial resonant encounter and long-term evolution. We find that the polar orbits are remarkably stable on ∼Gyr timescales. Inclination damping does not occur for the polar cases, although we do identify subpolar cases where it is important. We consider two case study polar Neptunes, WASP-107 b and HAT-P-11 b, and study them in the context of this theory, finding consistency with present-day properties if their tidal quality factors areQ≳ 104andQ≳ 105, respectively.more » « less
