Charting Galactic Accelerations: When and How to Extract a Unique Potential from the Distribution Function
arXiv:2106.05981 · doi:10.1093/mnras/stab2049
Abstract
The advent of datasets of stars in the Milky Way with six-dimensional phase-space information makes it possible to construct empirically the distribution function (DF). Here, we show that the accelerations can be uniquely determined from the DF using the collisionless Boltzmann equation, providing the Hessian determinant of the DF with respect to the velocities is non-vanishing. We illustrate this procedure and requirement with some analytic examples. Methods to extract the potential from datasets of discrete positions and velocities of stars are then discussed. Following Green & Ting (arXiv:2011.04673), we advocate the use of normalizing flows on a sample of observed phase-space positions to obtain a differentiable approximation of the DF. To then derive gravitational accelerations, we outline a semi-analytic method involving direct solutions of the over-constrained linear equations provided by the collisionless Boltzmann equation. Testing our algorithm on mock datasets derived from isotropic and anisotropic Hernquist models, we obtain excellent accuracies even with added noise. Our method represents a new, flexible and robust means of extracting the underlying gravitational accelerations from snapshots of six-dimensional stellar kinematics of an equilibrium system.
10 pages, 4 figures, MNRAS, in press
References in corpus (9)
- The Gaia mission
- Measuring the inclination and mass-to-light ratio of axisymmetric galaxies via anisotropic Jeans models of stellar kinematics
- The local high velocity tail and the Galactic escape speed
- Dipping Our Toes in the Water: First Models of GD-1 as a Stream
- Galactoseismology and the Local Density of Dark Matter
- The Stellar Distribution Function and Local Vertical Potential from Gaia DR2
- Temporal Normalizing Flows
- Deep Potential: Recovering the gravitational potential from a snapshot of phase space
- Reflection symmetries of Isolated Self-consistent Stellar Systems