Non-Parametric Spherical Jeans Mass Estimation with B-splines

Abstract

Spherical Jeans modeling is widely used to estimate mass profiles of systems from star clusters to galactic stellar haloes to clusters of galaxies. It derives the cumulative mass profile, M(<r), from kinematics of tracers of the potential under the assumptions of spherical symmetry and dynamical equilibrium. We consider the application of Jeans modeling to mapping the dark matter distribution in the outer reaches of the Milky Way using field halo stars. We present a novel non-parametric routine for solving the spherical Jeans equation by fitting B-splines to the velocity and density profiles of halo stars. While most implementations assume parametric forms for these profiles, B-splines provide non-parametric fitting curves with analytical derivatives. Our routine recovers the mass profiles of equilibrium systems with flattened haloes or a stellar disc and bulge excellently (<~ 10% error at most radii). Tests with non-equilibrium, Milky Way-like galaxies from the Latte suite of FIRE-2 simulations perform quite well (<~ 15% error for r <~ 100 kpc). We also create observationally motivated datasets for the Latte suite by imposing selection functions and errors on phase space coordinates characteristic of Gaia and the DESI Milky Way Survey. The resulting imprecise and incomplete data require us to introduce an MCMC-based subroutine to obtain deconvolved density and velocity dispersion profiles from the tracer population. With these observational effects taken into account, the accuracy of the Jeans mass estimate remains at the level 20% or better.

14 pages, 7 figures; accepted for publication in MNRAS; comments welcome

References in corpus (20)

• Array Programming with NumPy
• Gaia Data Release 1. Summary of the astrometric, photometric, and survey properties
• The DESI Experiment Part I: Science,Targeting, and Survey Design
• The mass distribution and gravitational potential of the Milky Way
• The Radial Velocity Experiment (RAVE): first data release
• The GALAH Survey: Scientific Motivation
• The Milky Way's Circular Velocity Curve to 60 kpc and an Estimate of the Dark Matter Halo Mass from Kinematics of ~2400 SDSS Blue Horizontal Branch Stars
• Gaia Data Release 1: Astrometry - one billion positions, two million proper motions and parallaxes
• Tracing Galaxy Formation with Stellar Halos II: Relating Substructure in Phase- and Abundance-Space to Accretion Histories
• APOGEE: The Apache Point Observatory Galactic Evolution Experiment
• The Mass Profile of the Galaxy to 80 kpc
• All-Sky Dynamical Response of the Galactic Halo to the Large Magellanic Cloud
• Detection of the Milky Way reflex motion due to the Large Magellanic Cloud infall
• Action-based distribution functions for spheroidal galaxy components
• Project overview and update on WEAVE: the next generation wide-field spectroscopy facility for the William Herschel Telescope
• When the Jeans don't fit: How stellar feedback drives stellar kinematics and complicates dynamical modeling in low-mass galaxies
• Bayesian Mass Estimates of the Milky Way: including measurement uncertainties with hierarchical Bayes
• A novel JEAnS analysis of the Fornax dwarf using evolutionary algorithms: mass follows light with signs of an off-centre merger
• Jeans that fit: weighing the mass of the Milky Way analogues in the $\Lambda \mathrm{CDM}$ universe
• Quantifying Departures from Equilibrium with the Spherical Jeans Equation