GravSphere2: A higher-order Jeans method for mass-modeling spherical stellar systems
arXiv:2509.24103 · doi:10.1051/0004-6361/202557195
Abstract
Mass-modeling methods are used to infer the gravitational field of stellar systems, from globular clusters to giant elliptical galaxies. While many methods exist, most require assumptions about the form of the underlying distribution function or data binning that leads to loss of information. With only line-of-sight (LOS) data, many methods suffer from the well-known mass-anisotropy degeneracy. To overcome these limitations, we develop a new, publicly available mass-modeling method, GravSphere2. This combines individual stellar velocities from LOS and proper motion (PM) measurements to solve the Jeans equations up to fourth order, without any data binning. Using flexible functional forms for the anisotropy profiles at second and fourth order, we show how including additional constraints from a new observable - fourth-order PMs - fully closes the system of equations, breaking the mass-anisotropy degeneracy at all orders. We test our method on mock data for dwarf galaxies, showing how GravSphere2 improves on previous methods. GravSphere2 recovers the mass density, stellar velocity anisotropy, and logarithmic slope of the mass density profile within its quoted 95% confidence intervals across almost all mocks over a wide radial range (0.1 < R/Rhalf < 10). We find GravSphere2 outperforms simple mass estimators, suggesting that it is worth using even when only a few LOS velocities are available. With 1,000 tracers without PMs, GravSphere2 recovers the logarithmic density slope at Rhalf with 12% (25%) statistical errors for cuspy (cored) mock data, enabling a distinction between the two. Including PMs, this improves to 8% (12%). With just 100 tracers and no PMs, we recover slopes with ~ 30% (20%) errors. GravSphere2 will be a valuable new tool to hunt for black holes and dark matter in spherical stellar systems, from globular clusters and dwarf galaxies to giant ellipticals and galaxy clusters.
Accepted for publication in A&A, corrected typos found in eqs. 17, 18, 20, 23
References in corpus (74)
- SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
- Array Programming with NumPy
- emcee: The MCMC Hammer
- MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics
- dynesty: A Dynamic Nested Sampling Package for Estimating Bayesian Posteriors and Evidences
- Small-Scale Challenges to the CDM Paradigm
- Intermediate-Mass Black Holes
- A Universal Mass Profile for Dwarf Spheroidal Galaxies?
- A Method for Measuring (Slopes of) the Mass Profiles of Dwarf Spheroidal Galaxies
- AGAMA: Action-based galaxy modelling architecture
- Measuring the inclination and mass-to-light ratio of axisymmetric galaxies via anisotropic Jeans models of stellar kinematics
- PolyChord: nested sampling for cosmology
- PolyChord: next-generation nested sampling
- The unexpected diversity of dwarf galaxy rotation curves
- Nuclear Star Clusters
- Dark matter cores all the way down
- Triaxial orbit based galaxy models with an application to the (apparent) decoupled core galaxy NGC 4365
- Dark matter in elliptical galaxies: II. Estimating the mass within the virial radius
- Dark matter distribution in the Coma cluster from galaxy kinematics: breaking the mass-anisotropy degeneracy
- Observational evidence for intermediate-mass black holes
- Dark matter heats up in dwarf galaxies
- New Limits on an Intermediate Mass Black Hole in Omega Centauri: I. Hubble Space Telescope Photometry and Proper Motions
- New Limits on an Intermediate Mass Black Hole in Omega Centauri: II. Dynamical Models
- MAMPOSSt: Modelling Anisotropy and Mass Profiles of Observed Spherical Systems. I. Gaussian 3D velocities
- Dark matter distribution in dwarf spheroidal galaxies
- NMAGIC: Fast Parallel Implementation of a Chi-Squared-Made-to-Measure Algorithm for Modelling Observational Data
- The case for a cold dark matter cusp in Draco
- Dearth of dark matter or massive dark halo? Mass-shape-anisotropy degeneracies revealed by NMAGIC dynamical models of the elliptical galaxy NGC 3379
- Systematics in virial mass estimators for pressure-supported systems
- Dark Matter in Dwarf Spheroidals I: Models
- A family of lowered isothermal models
- Baryonic clues to the puzzling diversity of dwarf galaxy rotation curves
- Internal kinematics and dynamical models of dwarf spheroidal galaxies around the Milky Way
- Dark matter distribution in the Draco dwarf from velocity moments
- Fast-moving stars around an intermediate-mass black hole in Omega Centauri
- Diversity of dark matter density profiles in the Galactic dwarf spheroidal satellites
- The Nuclear Cluster of the Milky Way: Total Mass and Luminosity (long version)
- How to break the density-anisotropy degeneracy in spherical stellar systems
- Orbit-based dynamical models of the Sculptor dSph galaxy
- On the dark matter profile in Sculptor: Breaking the degeneracy with Virial shape parameters
- Dark Matter Cores and Cusps in Spiral Galaxies and their Explanations
- Discrete dynamical modelling of omega Centauri
- A Cusp Slope -- Central Anisotropy Theorem
- Tailoring triaxial N-body models via a novel made-to-measure method
- Dynamical models with a general anisotropy profile
- Reconciling mass-estimates of ultra-diffuse galaxies
- Knowing the unknowns: uncertainties in simple estimators of galactic dynamical masses
- New Constraints on the Mass of Fermionic Dark Matter from Dwarf Spheroidal Galaxies
- Andromeda XXI -- a dwarf galaxy in a low density dark matter halo
- To beta or not to beta: can higher-order Jeans analysis break the mass-anisotropy degeneracy in simulated dwarfs?
- Analytical Solutions to the Mass-Anisotropy Degeneracy with Higher Order Jeans Analysis: A General Method
- Breaking Beta: A comparison of mass modelling methods for spherical systems
- Dwarf galaxies imply dark matter is heavier than
- Constraining the Mass Profiles of Stellar Systems: Schwarzschild Modeling of Discrete Velocity Datasets
- Stellar mass segregation as separating classifier between globular clusters and ultra-faint dwarf galaxies
- Confronting fuzzy dark matter with the rotation curves of nearby dwarf irregular galaxies
- Disc galaxy modelling with a particle-by-particle M2M method
- Kinematics of the Tucana Dwarf Galaxy: An Unusually Dense Dwarf in the Local Group
- Dark and luminous mass components of Omega Centauri with stellar kinematics
- Measuring galaxy cluster mass profiles into the low acceleration regime with galaxy kinematics
- Multimass modelling of Milky Way globular clusters -- I. Implications on their stellar initial mass function above 1 M
- An elusive dark central mass in the globular cluster M4
- HSTPROMO Internal Proper Motion Kinematics of Dwarf Spheroidal Galaxies: I. Velocity Anisotropy and Dark Matter Cusp Slope of Draco
- New constraints on the central mass contents of Omega Centauri from combined stellar kinematics and pulsar timing
- oMEGACat. VI. Analysis of the overall kinematics of Omega Centauri in 3D: velocity dispersion, kinematic distance, anisotropy, and energy equipartition
- Near-Gaussian distributions for modelling discrete stellar velocity data with heteroskedastic uncertainties
- The MUSE-Faint survey. III. Constraining scalar field dark matter with Antlia B
- Mass Modeling the Andromeda Dwarf Galaxies: Andromeda VI and Andromeda XXIII
- Trial by FIRE: Probing the dark matter density profile of dwarf galaxies with GraphNPE
- Leo I: the classical dwarf spheroidal galaxy with the highest dark-matter density
- The radial acceleration relation at the EDGE of galaxy formation: testing its universality in low-mass dwarf galaxies
- Galaxy mass profiles with convolutional neural networks
- Quantifying Departures from Equilibrium with the Spherical Jeans Equation
- Probing populations of dark stellar remnants in the globular clusters 47 Tuc and Terzan 5 using pulsar timing
Cited by in corpus (3)
- Multi-component, axisymmetric dynamical models of dSphs based on distribution functions: inferences on dark matter and intermediate-mass black holes in Draco and Ursa Minor
- Beyond single tracers: CNN-based inference of galaxy mass profiles from combined gas and stellar kinematics
- Broken Expectations: The Effects of Modelling Assumptions on the Inferred Dark Matter Distribution in the Milky Way's Satellites