Minimum-entropy constraints on galactic potentials

Abstract

A tracer sample in a gravitational potential, starting from a generic initial condition, phase-mixes towards a stationary state. This evolution is accompanied by an entropy increase, and the final state is characterized by a distribution function (DF) that depends only on integrals of motion (Jeans' theorem). We present a method to constrain a gravitational potential assuming a stationary (phase mixed) sample by minimizing the entropy the sample would have if it were allowed to phase-mix in trial potentials. This method avoids modeling the DF, and is applicable to any sets of integrals. We provide expressions for the entropy of DFs depending on energy, $f(E)$, energy and angular momentum, $f(E,L)$, or three actions, $f(\vec{J})$, and investigate the bias and statistical uncertainties in their estimates. We show that the method correctly recovers the parameters for spherical and axisymmetric potentials. We also present a methodology to characterize the posterior probability distribution of the parameters with an Approximate Bayesian Computation, indicating a pathway for application to observational data. Using $10^4$ tracers with $10\%(20\%)$-uncertainties in the 6D coordinates, we recover the flattening parameter $q$ of an axisymmetric potential with ${\sigma }_q/q\sim 5\%(10\%)$. The python module for the entropy estimators, \texttt{tropygal}, is made publicly available.

Accepted for publication by ApJ. In comparison to previous version, some change in the text and a few changes in the analysis, but the results and conclusions are the same

References in corpus (48)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• The Gaia mission
• The Apache Point Observatory Galactic Evolution Experiment (APOGEE)
• The mass distribution and gravitational potential of the Milky Way
• The GALAH Survey: Scientific Motivation
• Are the Magellanic Clouds on their First Passage about the Milky Way?
• AGAMA: Action-based galaxy modelling architecture
• The Masses of the Milky Way and Andromeda galaxies
• Actions for axisymmetric potentials
• A class of Rényi information estimators for multidimensional densities
• Nonequilibrium Statistical Mechanics of Systems with Long-Range Interactions: Ubiquity of Core-Halo Distributions
• Overview of the DESI Milky Way Survey
• Quantifying the impact of the Large Magellanic Cloud on the structure of the Milky Way's dark matter halo using Basis Function Expansions
• A review of action estimation methods for galactic dynamics
• Detection of the Milky Way reflex motion due to the Large Magellanic Cloud infall
• Detection of the LMC-induced sloshing of the Galactic halo
• Phase-space mixing and the merging of cusps
• cosmoabc: Likelihood-free inference via Population Monte Carlo Approximate Bayesian Computation
• The mass of the Milky Way out to 100 kpc using halo stars
• Action-space clustering of tidal streams to infer the Galactic potential
• A non-parametric k-nearest neighbour entropy estimator
• Approximate Bayesian Computation in Large Scale Structure: constraining the galaxy-halo connection
• Analysing surveys of our Galaxy -- II. Determining the potential
• Galactic potential constraints from clustering in action space of combined stellar stream data
• Hypervirial Models of Stellar Systems
• Action-based distribution function modelling for constraining the shape of the Galactic dark matter halo
• A statistical method for measuring the Galactic potential and testing gravity with cold tidal streams
• Analysing surveys of our Galaxy I: basic astrometric data
• Bayes versus the virial theorem: inferring the potential of a galaxy from a kinematical snapshot
• A Comprehensive Perturbative Formalism for Phase-Mixing in Perturbed Disks. I. Phase spirals in an Infinite, Isothermal Slab
• Constraining the Galactic potential via action-based distribution functions for mono-abundance stellar populations
• Orbital Torus Imaging: Using Element Abundances to Map Orbits and Mass in the Milky Way
• The orbital PDF: general inference of the gravitational potential from steady-state tracers
• Action-based Dynamical Modelling for the Milky Way Disk
• Orbital roulette: a new method of gravity estimation from observed motions
• Generalized entropy production in collisionless plasma flows and turbulence
• Deep Potential: Recovering the gravitational potential from a snapshot of phase space
• Eccentricity dynamics of wide binaries -- I. The effect of Galactic tides
• The discreteness-driven relaxation of collisionless gravitating systems: entropy evolution and the Nyquist-Shannon theorem
• Phase-space entropy cascade and irreversibility of stochastic heating in nearly collisionless plasma turbulence
• The discreteness-driven relaxation of collisionless gravitating systems: entropy evolution in external potentials, N-dependence and the role of chaos
• The Arrow of Time in the collapse of collisionless self-gravitating systems: non-validity of the Vlasov-Poisson equation during violent relaxation
• Non-Parametric Spherical Jeans Mass Estimation with B-splines
• Symplectic coarse graining approach to the dynamics of spherical self-gravitating systems
• Entropy Estimation via Uniformization
• Dimensional measures of generalized entropy
• Measuring the mass distribution in stellar systems

Cited by in corpus (1)

• Evolution of action-space coherence in a Milky Way-like simulation