GravSphere2: A higher-order Jeans method for mass-modeling spherical stellar systems
Pith reviewed 2026-05-18 12:34 UTC · model grok-4.3
The pith
GravSphere2 solves the Jeans equations to fourth order using unbinned stellar velocities and proper motions to recover mass density profiles while breaking the mass-anisotropy degeneracy.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
GravSphere2 combines unbinned line-of-sight and proper-motion velocities to solve the Jeans equations to fourth order. Flexible functional forms for the second- and fourth-order anisotropy profiles are adopted, and fourth-order proper-motion constraints close the system, eliminating the mass-anisotropy degeneracy at all orders. On mock data the method recovers the mass density, stellar velocity anisotropy, and logarithmic slope of the mass density profile within its quoted 95 percent confidence intervals across almost all mocks over 0.1 < R/Rhalf < 10.
What carries the argument
Flexible functional forms for the second- and fourth-order anisotropy profiles together with fourth-order proper-motion constraints that close the Jeans system without data binning.
If this is right
- With 1,000 tracers and no proper motions the logarithmic density slope at the half-light radius is recovered to 12 percent (25 percent) statistical error for cuspy (cored) mocks.
- Adding proper motions improves the slope errors to 8 percent (12 percent) for the same number of tracers.
- Even with only 100 tracers and no proper motions the slope is recovered to roughly 30 percent (20 percent) error.
- The method outperforms simple mass estimators and is therefore worth applying even when proper-motion data are unavailable.
- The recovered quantities remain reliable over nearly the full radial range 0.1 < R/Rhalf < 10 in almost all mocks.
Where Pith is reading between the lines
- The same fourth-order closure could be tested on real globular-cluster data to place tighter limits on intermediate-mass black holes than are possible with second-order Jeans modeling alone.
- Because the method works with unbinned data it can be applied directly to sparse samples from future wide-field surveys without first grouping stars into radial bins.
- Extending the same logic to non-spherical systems would require replacing the spherical Jeans equations with the axisymmetric or triaxial versions while retaining the higher-order proper-motion constraints.
- The recovered density slopes at the half-light radius could be compared directly with independent estimates from strong lensing or X-ray gas in massive ellipticals to test consistency across methods.
Load-bearing premise
That the chosen flexible functional forms for the anisotropy profiles introduce no significant bias when they are used with fourth-order proper-motion data to close the Jeans equations.
What would settle it
Running GravSphere2 on mock data sets whose true anisotropy profiles lie outside the family of flexible forms adopted in the method and checking whether the recovered mass densities fall outside the quoted 95 percent confidence intervals over a substantial fraction of the radial range.
Figures
read the original 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces GravSphere2, a publicly available method for mass-modeling spherical stellar systems that solves the Jeans equations to fourth order using unbinned individual stellar line-of-sight and proper-motion velocities. It employs flexible functional forms for the second- and fourth-order anisotropy profiles together with fourth-order proper-motion moments to close the system and break the mass-anisotropy degeneracy. Mock tests on dwarf-galaxy data demonstrate recovery of the mass density, stellar velocity anisotropy, and logarithmic slope of the mass-density profile within the quoted 95% confidence intervals across almost all mocks over the radial range 0.1 < R/R_half < 10, with improved performance relative to simple mass estimators even for modest numbers of tracers (100–1000).
Significance. If the recovery statistics hold under the stated assumptions, GravSphere2 provides a practical advance for dynamical modeling by extending Jeans analysis to fourth order without binning and by incorporating proper-motion constraints. The public code release and the quantitative demonstration that the logarithmic slope at R_half can be recovered to 8–12% (with PMs) or 12–25% (without) for cuspy versus cored profiles constitute clear strengths that would make the tool useful for constraining dark matter and black holes in globular clusters, dwarfs, and larger spherical systems.
major comments (2)
- [Section 3] Section 3 (method development): the claim that the chosen flexible functional forms for the second- and fourth-order anisotropy profiles, together with fourth-order PM moments, fully close the Jeans system and break the degeneracy without significant bias from the parametrization choice is load-bearing for the central recovery result. The manuscript should include explicit tests with alternative parametrizations or a non-parametric anisotropy model to quantify any residual bias in the recovered mass density and logarithmic slope.
- [Mock-data tests section] Mock-data tests section: the 95% confidence intervals are reported to contain the true values for most mocks, yet the details of error propagation for the fourth-order moments and any post-hoc adjustments to the fitting procedure are not fully specified. These elements directly affect the reliability of the quoted recovery fractions and should be documented with the relevant equations or pseudocode.
minor comments (2)
- [Abstract] Abstract: the phrase 'across almost all mocks' should be replaced by the exact fraction or number of mocks in which recovery falls outside the 95% intervals to give readers a precise performance metric.
- [Figures] Figure captions and text: ensure that all result panels explicitly label the radial range, mock type (cuspy/cored), and whether PM data are included so that the wide-range recovery claim can be verified at a glance.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for the constructive comments. We address each major point below and have revised the manuscript accordingly to improve clarity and robustness.
read point-by-point responses
-
Referee: [Section 3] Section 3 (method development): the claim that the chosen flexible functional forms for the second- and fourth-order anisotropy profiles, together with fourth-order PM moments, fully close the Jeans system and break the degeneracy without significant bias from the parametrization choice is load-bearing for the central recovery result. The manuscript should include explicit tests with alternative parametrizations or a non-parametric anisotropy model to quantify any residual bias in the recovered mass density and logarithmic slope.
Authors: We appreciate the referee's emphasis on this central claim. The functional forms were chosen for their flexibility and physical motivation, and the inclusion of fourth-order PM moments is intended to close the system at all orders. To directly quantify any residual bias arising from the specific parametrization, we will add explicit tests in the revised manuscript using an alternative parametrization for the fourth-order anisotropy profile. These tests will report the impact on recovered mass density and logarithmic slope, allowing readers to assess the robustness of the results. revision: yes
-
Referee: [Mock-data tests section] Mock-data tests section: the 95% confidence intervals are reported to contain the true values for most mocks, yet the details of error propagation for the fourth-order moments and any post-hoc adjustments to the fitting procedure are not fully specified. These elements directly affect the reliability of the quoted recovery fractions and should be documented with the relevant equations or pseudocode.
Authors: We agree that greater transparency on these technical details is needed. In the revised manuscript we will add the explicit equations for computing the fourth-order moments and propagating their uncertainties into the likelihood function. We will also include pseudocode for the overall fitting procedure in an appendix so that the error treatment and any adjustments are fully reproducible. revision: yes
Circularity Check
No significant circularity in derivation chain
full rationale
The paper introduces GravSphere2 as a new method extending the Jeans equations to fourth order, employing flexible parametrizations of second- and fourth-order anisotropy profiles together with fourth-order proper-motion moments to close the system and break the mass-anisotropy degeneracy. All central claims concern recovery performance on independently generated mock data sets whose true mass profiles, anisotropies, and density slopes are known a priori and are not derived from the same functional forms or fitted parameters used in the modeling. No self-definitional steps, fitted inputs relabeled as predictions, or load-bearing self-citations appear in the provided material; the validation statistics are reported directly from the mock tests rather than being forced by construction. The derivation therefore remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (1)
- parameters in flexible anisotropy functional forms
axioms (1)
- domain assumption The stellar system is spherical and in steady-state dynamical equilibrium
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
solve the Jeans equations up to fourth order... flexible functional forms for the anisotropy profiles at second and fourth order... fourth-order PMs
-
IndisputableMonolith/Foundation/DimensionForcing.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
coreNFWtides profile... generalized αβγ tracer density
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Alvarez, A., Calore, F., Genina, A., et al. 2020, JCAP, 2020, 004
work page 2020
-
[2]
Alvey, J., Sabti, N., Tiki, V ., et al. 2021, MNRAS, 501, 1188
work page 2021
-
[3]
An, J. H. & Evans, N. W. 2006, ApJ, 642, 752
work page 2006
-
[4]
Anderson, J. & Van der Marel, R. P. 2010, The Astrophysical Journal, 710, 1032 Article number, page 15 of 19 A&A proofs:manuscript no. main
work page 2010
-
[5]
M., Pascale, R., Battaglia, G., et al
Arroyo-Polonio, J. M., Pascale, R., Battaglia, G., et al. 2025, A&A, 699, A347
work page 2025
-
[6]
Askar, A., Baldassare, V . F., & Mezcua, M. 2023, arXiv:2311.12118 Bañares-Hernández, A., Calore, F., Martin Camalich, J., & Read, J. I. 2025, A&A, 693, A104 Bañares-Hernández, A., Castillo, A., Martin Camalich, J., & Iorio, G. 2023, A&A, 676, A63
- [7]
-
[8]
2013, New Astronomy Reviews, 57, 52
Battaglia, G., Helmi, A., & Breddels, M. 2013, New Astronomy Reviews, 57, 52
work page 2013
- [9]
-
[10]
Baumgardt, H., Faller, J., Meinhold, N., McGovern-Greco, C., & Hilker, M. 2022, MNRAS, 510, 3531
work page 2022
- [11]
-
[12]
Theia: Faint objects in motion or the new astrometry frontier
Boehm, C., Krone-Martins, A., Amorim, A., et al. 2017, arXiv:1707.01348
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[13]
1978, PhD thesis, University of Groningen, Netherlands
Bosma, A. 1978, PhD thesis, University of Groningen, Netherlands
work page 1978
-
[14]
A., Helmi, A., van den Bosch, R
Breddels, M. A., Helmi, A., van den Bosch, R. C. E., van de Ven, G., & Battaglia, G. 2013, MNRAS, 433, 3173
work page 2013
-
[15]
Bullock, J. S. & Boylan-Kolchin, M. 2017, ARA&A, 55, 343
work page 2017
-
[16]
Campbell, D. J. R., Frenk, C. S., Jenkins, A., et al. 2017, MNRAS, 469, 2335
work page 2017
- [17]
-
[18]
Cappellari, M. 2015, arXiv:1504.05533 Chanamé, J., Kleyna, J., & van der Marel, R. 2008, ApJ, 682, 841
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[19]
Collins, M. L. M., Read, J. I., Ibata, R. A., et al. 2021, MNRAS, 505, 5686 De Leo, M., Read, J. I., Noël, N. E. D., et al. 2024, MNRAS, 535, 1015 de Lorenzi, F., Debattista, V . P., Gerhard, O., & Sambhus, N. 2007, MNRAS, 376, 71 de Lorenzi, F., Gerhard, O., Coccato, L., et al. 2009, MNRAS, 395, 76
work page 2021
- [20]
- [21]
- [22]
-
[23]
Dickson, N., Hénault-Brunet, V ., Baumgardt, H., Gieles, M., & Smith, P. J. 2023, MNRAS, 522, 5320
work page 2023
-
[24]
J., Hénault-Brunet, V ., Gieles, M., & Baumgardt, H
Dickson, N., Smith, P. J., Hénault-Brunet, V ., Gieles, M., & Baumgardt, H. 2024, MNRAS, 529, 331
work page 2024
-
[25]
Errani, R., Peñarrubia, J., & Walker, M. G. 2018, MNRAS, 481, 5073
work page 2018
- [26]
-
[27]
Evslin, J. & Del Popolo, A. 2017, ApJ, 841, 90 Expósito-Márquez, J., Brook, C. B., Huertas-Company, M., et al. 2023, MN- RAS, 519, 4384
work page 2017
- [28]
-
[29]
2016, The Journal of Open Source Software, 1, 24
Foreman-Mackey, D. 2016, The Journal of Open Source Software, 1, 24
work page 2016
-
[30]
Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, Publications of the Astronomical Society of the Pacific, 125, 306
work page 2013
-
[31]
K., Chatzopoulos, S., Gerhard, O., et al
Fritz, T. K., Chatzopoulos, S., Gerhard, O., et al. 2016, ApJ, 821, 44
work page 2016
-
[32]
Genina, A., Read, J. I., Fattahi, A., & Frenk, C. S. 2022, MNRAS, 510, 2186
work page 2022
- [33]
- [34]
-
[35]
D., Das, P., Heber, D., & Izzard, R
Gration, A., Hendriks, D. D., Das, P., Heber, D., & Izzard, R. G. 2025, arXiv:2509.14316
-
[36]
Greene, J. E., Strader, J., & Ho, L. C. 2020, ARA&A, 58, 257
work page 2020
-
[37]
Gregory, A. L., Collins, M. L. M., Read, J. I., et al. 2019, MNRAS, 485, 2010 Häberle, M., Neumayer, N., Clontz, C., et al. 2025, ApJ, 983, 95 Häberle, M., Neumayer, N., Seth, A., et al. 2024, Nature, 631, 285
work page 2019
-
[38]
Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357–362
work page 2020
- [39]
- [40]
-
[41]
Higson, E., Handley, W., Hobson, M., & Lasenby, A. 2019, SC, 29, 891
work page 2019
-
[42]
Hunt, J. A. S. & Kawata, D. 2013, MNRAS, 430, 1928
work page 2013
-
[43]
Hunter, J. D. 2007, Computing in Science & Engineering, 9, 90
work page 2007
- [44]
- [45]
-
[46]
Kluyver, T., Ragan-Kelley, B., Pérez, F., et al. 2016, in IOS Press, 87–90
work page 2016
-
[47]
2022, joshspeagle/dynesty: v1.2.3
Koposov, S., Speagle, J., Barbary, K., et al. 2022, joshspeagle/dynesty: v1.2.3
work page 2022
-
[48]
Laporte, C. F. P., Agnello, A., & Navarro, J. F. 2019, MNRAS, 484, 245
work page 2019
- [49]
- [50]
-
[51]
H., Hayashi, K., Horigome, S., Matsumoto, S., & Nojiri, M
Lim, S. H., Hayashi, K., Horigome, S., Matsumoto, S., & Nojiri, M. M. 2025, arXiv:2505.00763 Łokas, E. L. 2002, MNRAS, 333, 697 Łokas, E. L. & Mamon, G. A. 2003, MNRAS, 343, 401 Łokas, E. L., Mamon, G. A., & Prada, F. 2005, MNRAS, 363, 918
-
[52]
Mamon, G. A., Biviano, A., & Boué, G. 2013, MNRAS, 429, 3079
work page 2013
-
[53]
Mamon, G. A. & Łokas, E. L. 2005, MNRAS, 363, 705
work page 2005
- [54]
-
[55]
Merrifield, M. R. & Kent, S. M. 1990, AJ, 99, 1548
work page 1990
- [56]
- [57]
-
[58]
2017, International Journal of Modern Physics D, 26, 1730021
Mezcua, M. 2017, International Journal of Modern Physics D, 26, 1730021
work page 2017
- [59]
-
[60]
Neal, R. M. 2003, The Annals of Statistics, 31, 705
work page 2003
- [61]
-
[62]
Nguyen, T., Read, J., Necib, L., et al. 2025, arXiv:2503.03812
-
[63]
Oman, K. A., Navarro, J. F., Fattahi, A., et al. 2015, MNRAS, 452, 3650
work page 2015
-
[64]
Osipkov, L. P. 1979, Pisma v Astronomicheskii Zhurnal, 5, 77
work page 1979
-
[65]
Pascale, R., Nipoti, C., Calura, F., & Della Croce, A. 2025, A&A, 700, A77
work page 2025
-
[66]
Pickett, C. S., Collins, M. L. M., Rich, R. M., et al. 2025, MNRAS, 540, 1701
work page 2025
-
[67]
Read, J. I., Agertz, O., & Collins, M. L. M. 2016, MNRAS, 459, 2573
work page 2016
- [68]
-
[69]
Read, J. I. & Steger, P. 2017, MNRAS, 471, 4541
work page 2017
- [70]
- [71]
- [72]
- [73]
- [74]
-
[75]
Sanders, J. L. & Evans, N. W. 2020, MNRAS, 499, 5806
work page 2020
- [76]
- [77]
-
[78]
2004, AIP Conference Series, 735, 395
Skilling, J. 2004, AIP Conference Series, 735, 395
work page 2004
- [79]
-
[80]
Speagle, J. S. 2020, MNRAS, 493, 3132
work page 2020
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.