skip to main content

Title: A Comprehensive Perturbative Formalism for Phase Mixing in Perturbed Disks. I. Phase Spirals in an Infinite, Isothermal Slab

Galactic disks are highly responsive systems that often undergo external perturbations and subsequent collisionless equilibration, predominantly via phase mixing. We use linear perturbation theory to study the response of infinite isothermal slab analogs of disks to perturbations with diverse spatiotemporal characteristics. Without self-gravity of the response, the dominant Fourier modes that get excited in a disk are the bending and breathing modes, which, due to vertical phase mixing, trigger local phase-space spirals that are one- and two-armed, respectively. We demonstrate how the lateral streaming motion of slab stars causes phase spirals to damp out over time. The ratio of the perturbation timescale (τP) to the local, vertical oscillation time (τz) ultimately decides which of the two modes is excited. Faster, more impulsive (τP<τz) and slower, more adiabatic (τP>τz) perturbations excite stronger breathing and bending modes, respectively, although the response to very slow perturbations is exponentially suppressed. For encounters with satellite galaxies, this translates to more distant and more perpendicular encounters triggering stronger bending modes. We compute the direct response of the Milky Way disk to several of its satellite galaxies and find that recent encounters with all of them excite bending modes in the solar neighborhood. The encounter with Sagittarius triggers a response that is at least 1–2 orders of magnitude larger than that due to any other satellite, including the Large Magellanic Cloud. We briefly discuss how ignoring the presence of a dark matter halo and the self-gravity of the response might impact our conclusions.

more » « less
Author(s) / Creator(s):
; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Medium: X Size: Article No. 135
Article No. 135
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We develop a linear perturbative formalism to compute the response of an inhomogeneous stellar disk embedded in a nonresponsive dark matter (DM) halo to various perturbations like bars, spiral arms, and encounters with satellite galaxies. Without self-gravity to reinforce it, the response of a Fourier mode phase mixes away due to an intrinsic spread in the vertical (Ωz), radial (Ωr), and azimuthal (Ωϕ) frequencies, triggering local phase-space spirals. The detailed galactic potential dictates the shape of phase spirals: phase mixing occurs more slowly and thus phase spirals are more loosely wound in the outer disk and in the presence of an ambient DM halo. Collisional diffusion due to scattering of stars by structures like giant molecular clouds causes superexponential damping of the phase spiral amplitude. Thezvzphase spiral is one-armed (two-armed) for vertically antisymmetric (symmetric) bending (breathing) modes. Only transient perturbations with timescales (τP) comparable to the vertical oscillation period (τz∼ 1/Ωz) can trigger vertical phase spirals. Each (n,l,m) mode of the response to impulsive (τP<τ= 1/(nΩz+lΩr+mΩϕ)) perturbations is power-law (∼τP/τ) suppressed, but that to adiabatic (τP>τ) perturbations is exponentially weak (expτP/τα) except for resonant (τ→ ∞ ) modes. Slower (τP>τz) perturbations, e.g., distant encounters with satellite galaxies, induce stronger bending modes. Sagittarius (Sgr) dominates the solar neighborhood response of the Milky Way (MW) disk to satellite encounters. Thus, if the Gaia phase spiral was triggered by a MW satellite, Sgr is the leading contender. However, the survival of the phase spiral against collisional damping necessitates an impact ∼0.6–0.7 Gyr ago.

    more » « less
  2. Abstract

    Signatures of vertical disequilibrium have been observed across the Milky Way’s (MW’s) disk. These signatures manifest locally as unmixed phase spirals inzvzspace (“snails-in-phase”), and globally as nonzero meanzandvz, wrapping around the disk into physical spirals in thexyplane (“snails-in-space”). We explore the connection between these local and global spirals through the example of a satellite perturbing a test-particle MW-like disk. We anticipate our results to broadly apply to any vertical perturbation. Using azvzasymmetry metric, we demonstrate that in test-particle simulations: (a) multiple local phase-spiral morphologies appear when stars are binned by azimuthal actionJϕ, excited by a single event (in our case, a satellite disk crossing); (b) these distinct phase spirals are traced back to distinct disk locations; and (c) they are excited at distinct times. Thus, local phase spirals offer a global view of the MW’s perturbation history from multiple perspectives. Using a toy model for a Sagittarius (Sgr)–like satellite crossing the disk, we show that the full interaction takes place on timescales comparable to orbital periods of disk stars withinR≲ 10 kpc. Hence such perturbations have widespread influence, which peaks in distinct regions of the disk at different times. This leads us to examine the ongoing MW–Sgr interaction. While Sgr has not yet crossed the disk (currently,zSgr≈ −6 kpc,vz,Sgr≈ 210 km s−1), we demonstrate that the peak of the impact has already passed. Sgr’s pull over the past 150 Myr creates a globalvzsignature with amplitude ∝MSgr, which might be detectable in future spectroscopic surveys.

    more » « less

    The second data release of ESA’s Gaia mission revealed numerous signatures of disequilibrium in the Milky Way’s disc. These signatures are seen in the planar kinematics of stars, which manifest as ridges and ripples in R–vϕ, and in vertical kinematics, where a prominent spiral is seen in the z–vz phase space. In this work, we show an equivalent ΔR–vR phase spiral forms following a perturbation to the disc. We demonstrate the behaviour of the ΔR–vR phase spirals in both a toy model and a high-resolution N-body simulation of a satellite interaction. We then confront these models with the data, where we find partial ΔR–vR phase spirals in the Solar neighbourhood using the most recent data from Gaia DR3. This structure indicates ongoing radial phase mixing in the Galactic disc, suggesting a history of recent perturbations, either through internal or external (e.g. satellite) processes. Future work modelling the z–vz and ΔR–vR phase spirals in tandem may help break degeneracy’s between possible origins of the perturbation.

    more » « less

    Gaia Data Release 2 revealed that the Milky Way contains significant indications of departures from equilibrium in the form of asymmetric features in the phase space density of stars in the Solar neighbourhood. One such feature is the z–vz phase spiral, interpreted as the response of the disc to the influence of a perturbation perpendicular to the disc plane, which could be external (e.g. a satellite) or internal (e.g. the bar or spiral arms). In this work, we use Gaia Data Release 3 to dissect the phase spiral by dividing the local data set into groups with similar azimuthal actions, Jϕ, and conjugate angles, θϕ, which selects stars on similar orbits and at similar orbital phases, thus having experienced similar perturbations in the past. These divisions allow us to explore areas of the Galactic disc larger than the surveyed region. The separation improves the clarity of the z–vz phase spiral and exposes changes to its morphology across the different action-angle groups. In particular, we discover a transition to two armed ‘breathing spirals’ in the inner Milky Way. We conclude that the local data contain signatures of not one, but multiple perturbations with the prospect to use their distinct properties to infer the properties of the interactions that caused them.

    more » « less
  5. Abstract

    We analyze quiet‐time data from the Gravity Field and Ocean Circulation Explorer satellite as it overpassed the Southern Andes atz≃275 km on 5 July 2010 at 23 UT. We extract the 20 largest traveling atmospheric disturbances from the density perturbations and cross‐track winds using Fourier analysis. Using gravity wave (GW) dissipative theory that includes realistic molecular viscosity, we search parameter space to determine which hot spot traveling atmospheric disturbances are GWs. This results in the identification of 17 GWs having horizontal wavelengthsλH = 170–1,850 km, intrinsic periodsτIr = 11–54 min, intrinsic horizontal phase speedscIH = 245–630 m/s, and density perturbations 0.03–7%. We unambiguously determine the propagation direction for 11 of these GWs and find that most had large meridional components to their propagation directions. Using reverse ray tracing, we find that 10 of these GWs must have been created in the mesosphere or thermosphere. We show that mountain waves (MWs) were observed in the stratosphere earlier that day and that these MWs saturated atz∼ 70–75 km from convective instability. We suggest that these 10 Gravity Field and Ocean Circulation Explorer hot spot GWs are likely tertiary (or higher‐order) GWs created from the dissipation of secondary GWs excited by the local body forces created from MW breaking. We suggest that the other GW is likely a secondary or tertiary (or higher‐order) GW. This study strongly suggests that the hot spot GWs over the Southern Andes in the quiet‐time middle winter thermosphere cannot be successfully modeled by conventional global circulation models where GWs are parameterized and launched in the troposphere or stratosphere.

    more » « less