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.

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

References in corpus (61)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• Astropy: A Community Python Package for Astronomy
• emcee: The MCMC Hammer
• Gaia Data Release 2. Summary of the contents and survey properties
• The Astropy Project: Building an inclusive, open-science project and status of the v2.0 core package
• Gaia Data Release 1. Summary of the astrometric, photometric, and survey properties
• The Galaxy in Context: Structural, Kinematic and Integrated Properties
• Galaxies on FIRE (Feedback In Realistic Environments): Stellar Feedback Explains Cosmologically Inefficient Star Formation
• The DESI Experiment Part I: Science,Targeting, and Survey Design
• GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
• GAIA: Composition, Formation and Evolution of the Galaxy
• The mass distribution and gravitational potential of the Milky Way
• Tracing Galaxy Formation with Stellar Halos I: Methods
• The Radial Velocity Experiment (RAVE): first data release
• Resolving Cosmic Structure Formation with the Millennium-II Simulation
• 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
• Reconciling dwarf galaxies with LCDM cosmology: Simulating a realistic population of satellites around a Milky Way-mass galaxy
• Gaia Data Release 1: Astrometry - one billion positions, two million proper motions and parallaxes
• The APOSTLE simulations: solutions to the Local Group's cosmic puzzles
• AGAMA: Action-based galaxy modelling architecture
• Tracing Galaxy Formation with Stellar Halos II: Relating Substructure in Phase- and Abundance-Space to Accretion Histories
• The Masses of the Milky Way and Andromeda galaxies
• The radial velocity dispersion profile of the Galactic halo: Constraining the density profile of the dark halo of the Milky Way
• The APOSTLE project: Local Group kinematic mass constraints and simulation candidate selection
• 4MOST - 4-metre Multi-Object Spectroscopic Telescope
• APOGEE: The Apache Point Observatory Galactic Evolution Experiment
• The Mass Profile of the Galaxy to 80 kpc
• The shape of the Galactic halo with $Gaia$ DR2 RR Lyrae. Anatomy of an ancient major merger
• The origin of the diverse morphologies and kinematics of Milky Way-mass galaxies in the FIRE-2 simulations
• The mass of our Milky Way
• On determining the shape of matter distributions
• Synthetic Gaia surveys from the FIRE cosmological simulations of Milky Way-mass galaxies
• All-Sky Dynamical Response of the Galactic Halo to the Large Magellanic Cloud
• Estimating the dark matter halo mass of our Milky Way using dynamical tracers
• Detection of the Milky Way reflex motion due to the Large Magellanic Cloud infall
• The Gravitational Force Field of the Galaxy Measured From the Kinematics of RR Lyrae in Gaia
• Detection of the LMC-induced sloshing of the Galactic halo
• How to break the density-anisotropy degeneracy in spherical stellar systems
• 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
• The cumulative mass profile of the Milky Way as determined by globular cluster kinematics from Gaia DR2
• Discrete dynamical modelling of omega Centauri
• The mass of the Milky Way out to 100 kpc using halo stars
• Equilibrium models of the Milky Way mass are biased high by the LMC
• The formation of shell galaxies similar to NGC 7600 in the cold dark matter cosmogony
• When the Jeans don't fit: How stellar feedback drives stellar kinematics and complicates dynamical modeling in low-mass galaxies
• Self-consistent flattened isochrone models
• Preliminary Target Selection for the DESI Milky Way Survey (MWS)
• Constraining the Milky Way Mass Profile with Phase-Space Distribution of Satellite Galaxies
• To beta or not to beta: can higher-order Jeans analysis break the mass-anisotropy degeneracy in simulated dwarfs?
• First Gaia dynamical model of the Milky Way disc with six phase space coordinates: a test for galaxy dynamics
• Action-based distribution function modelling for constraining the shape of the Galactic dark matter halo
• Bayesian Mass Estimates of the Milky Way: including measurement uncertainties with hierarchical Bayes
• What to expect from dynamical modelling of galactic haloes II: the spherical Jeans equation
• Estimating the Milky Way's Mass via Hierarchical Bayes: A Blind Test on MUGS2 Simulated Galaxies
• A novel JEAnS analysis of the Fornax dwarf using evolutionary algorithms: mass follows light with signs of an off-centre merger
• Reliability of the measured velocity anisotropy of the Milky Way stellar halo
• 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

Cited by in corpus (7)

• Overview of the DESI Milky Way Survey
• Milky Way Mass with K Giants and BHB Stars Using LAMOST, SDSS/SEGUE, and Gaia: 3D Spherical Jeans Equation and Tracer Mass Estimator
• What to expect from dynamical modelling of cluster haloes II. Investigating dynamical state indicators with Random Forest
• Measuring Galactic Dark Matter through Unsupervised Machine Learning
• Is the core-cusp problem a matter of perspective: Jeans Anisotropic Modeling against numerical simulations
• Local dark matter density from Gaia DR3 K-dwarfs using Gaussian processes
• Minimum-entropy constraints on galactic potentials