Recognition: unknown
The Pulsar Radial Acceleration Relation
Pith reviewed 2026-05-08 07:30 UTC · model grok-4.3
The pith
Pulsar timing accelerations fit the three-dimensional radial acceleration relation from galaxies better than Newtonian baryonic gravity, but the comparison is dominated by the solar acceleration value.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By comparing the measured SPARC RAR generalized to 3D with literature accelerations for 26 binary-system pulsars, the data give a reduced chi-squared of 3.58, versus 10.86 for Newtonian baryonic gravity alone. Setting all accelerations equal to the solar value instead produces a reduced chi-squared of 3.75, demonstrating that the vector RAR test with current pulsar timing data is dominated by the solar acceleration.
What carries the argument
The vector generalization of the radial acceleration relation, g_obs = nu(|g_bar|) g_bar, where nu is taken from SPARC galaxy data and applied to the line-of-sight differential accelerations between the Sun and pulsars.
If this is right
- The RAR holds for accelerations measured along lines of sight outside the galactic disk.
- Pulsar timing supplies an independent test of the acceleration relation in the Milky Way environment.
- The vector form of the relation can be probed in three dimensions using differential accelerations.
- If additional pulsars yield smaller chi-squared values under the RAR than under the solar reference, the test will become constraining rather than solar-dominated.
Where Pith is reading between the lines
- Future catalogs of binary pulsars at larger average distances from the Sun could reduce the dominance of the solar acceleration and allow the shape of the interpolating function to be tested locally.
- The method could be cross-checked against accelerations inferred from stellar streams or globular clusters in the same volume to test consistency across different tracers.
- If the relation continues to hold, it would imply that the acceleration law is insensitive to the specific galactic environment sampled by pulsars.
- The near-equivalence of the RAR and solar-only fits suggests that expanding the sample size or precision will be necessary before the test can distinguish between competing explanations for the RAR.
Load-bearing premise
The line-of-sight differential accelerations from pulsar timing can be directly compared to the three-dimensional vector form of the RAR without large projection effects, unaccounted systematics, or deviations from the SPARC relation in the Milky Way.
What would settle it
Measuring independent accelerations for pulsars at distances much larger than the Sun's galactocentric radius and finding that their values deviate systematically from the RAR prediction while remaining inconsistent with the solar acceleration.
Figures
read the original abstract
The radial acceleration relation (RAR) links observed and baryonic accelerations, and is best established in rotation curves of late-type galaxies. Pulsar timing, which measures line-of-sight (LOS) differential accelerations between the Sun and pulsars, provides a novel probe of this relation, including along directions outside the Galactic disc. By combining these pulsar differential accelerations with the acceleration at the Sun, we test whether current pulsar timing data carry information on a vector generalisation of the RAR, ${g}_{\rm obs}=\nu(|{g}_{\rm bar}|){g}_{\rm bar}$. Comparing the measured SPARC RAR (generalised to 3D) to 26 binary-system pulsars with literature accelerations, we find a reduced $\chi^2$ of 3.58, compared with 10.86 for Newtonian baryonic gravity alone. However, setting all accelerations to that of the Sun gives a reduced $\chi^2$ of 3.75, showing that this vector RAR test is dominated by the Solar acceleration with current data.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes using line-of-sight differential accelerations measured via pulsar timing for 26 binary pulsars as a test of a 3D vector generalization of the radial acceleration relation (RAR) established by SPARC galaxy data. Combining these with the Solar acceleration, the authors report a reduced χ² of 3.58 for the RAR model versus 10.86 for Newtonian baryonic gravity alone; however, a constant-acceleration model fixed to the Solar value yields χ² = 3.75, leading to the conclusion that current pulsar data are dominated by the Solar term and do not yet provide a discriminating test.
Significance. If the statistical comparison holds after full data release, the work supplies an independent, out-of-plane probe of the RAR that is orthogonal to rotation-curve analyses. It correctly imports the SPARC relation without circular fitting and demonstrates that pulsar timing can in principle constrain vector generalizations of the relation. The present sample size and precision, however, render the test inconclusive, so the immediate scientific impact is modest but the methodological direction is useful for future larger datasets.
major comments (2)
- [Abstract / Results] The abstract (and presumably the results section) presents reduced χ² values of 3.58, 10.86 and 3.75 without tabulating the 26 individual pulsar accelerations, their uncertainties, the precise selection cuts, or the covariance matrix used in the fit. Without these, the claim that the RAR model is only marginally better than the constant-Solar model cannot be verified and the statistical robustness remains unassessed.
- [Methods / Discussion] The direct mapping of observed line-of-sight differential accelerations onto the vector RAR g_obs = ν(|g_bar|) g_bar assumes negligible projection effects, negligible deviations of the Milky Way from the SPARC relation, and that the Solar acceleration vector is known to sufficient precision. These assumptions are load-bearing for the conclusion that the test is “dominated by the Solar acceleration”; a quantitative assessment of projection biases and environment-dependent scatter is required.
minor comments (1)
- [Introduction] The notation for the vector generalization (g_obs = ν(|g_bar|) g_bar) is introduced without an explicit equation number or reference to the precise functional form of ν taken from SPARC; adding an equation label would improve clarity.
Simulated Author's Rebuttal
We thank the referee for the careful review and constructive comments on our manuscript. We address each major comment point by point below.
read point-by-point responses
-
Referee: [Abstract / Results] The abstract (and presumably the results section) presents reduced χ² values of 3.58, 10.86 and 3.75 without tabulating the 26 individual pulsar accelerations, their uncertainties, the precise selection cuts, or the covariance matrix used in the fit. Without these, the claim that the RAR model is only marginally better than the constant-Solar model cannot be verified and the statistical robustness remains unassessed.
Authors: We agree that the underlying data should be tabulated to permit independent verification. In the revised manuscript we will add a table in the results section listing the 26 pulsars together with their measured line-of-sight differential accelerations, uncertainties, the exact selection cuts applied, and a description of the covariance matrix used for the χ² fits. revision: yes
-
Referee: [Methods / Discussion] The direct mapping of observed line-of-sight differential accelerations onto the vector RAR g_obs = ν(|g_bar|) g_bar assumes negligible projection effects, negligible deviations of the Milky Way from the SPARC relation, and that the Solar acceleration vector is known to sufficient precision. These assumptions are load-bearing for the conclusion that the test is “dominated by the Solar acceleration”; a quantitative assessment of projection biases and environment-dependent scatter is required.
Authors: The referee correctly notes that the mapping rests on several assumptions. The near-equality of reduced χ² for the RAR model (3.58) and the constant-Solar model (3.75) already indicates that the Solar term dominates with the present sample. In revision we will add a short quantitative discussion of projection effects, estimating the typical misalignment angle for the selected pulsars and the resulting bias on the inferred acceleration. We will also state that the Solar acceleration is taken from the most precise literature values and that Milky-Way deviations from the SPARC relation are expected to be modest. A full Monte-Carlo treatment of environment-dependent scatter lies beyond the scope of this work but will be pursued with future larger datasets. revision: partial
Circularity Check
No circularity: external SPARC RAR tested on independent pulsar data
full rationale
The paper imports the SPARC RAR (and its 3D generalization) from external literature and performs a direct statistical comparison (reduced χ²) against 26 independent pulsar timing accelerations. No equation defines the RAR using pulsar data, fits parameters to the pulsars then renames them as predictions, or relies on self-citations for load-bearing uniqueness or ansatz. The internal check that Solar acceleration alone yields similar χ² is an explicit caveat, not a circular reduction. The derivation chain is a standard external test and remains self-contained against benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The SPARC radial acceleration relation generalizes to a vector form g_obs = nu(|g_bar|) g_bar that can be applied to line-of-sight pulsar accelerations in the Milky Way.
Reference graph
Works this paper leans on
-
[1]
Vertical waves in the solar neighbourhood in Gaia DR2
Bennett, M., & Bovy, J. 2019, MNRAS, 482, 1417, doi: 10.1093/mnras/sty2813
-
[2]
2016, ARA&A, 54, 529, doi: 10.1146/annurev-astro-081915-023441
Bland-Hawthorn, J., & Gerhard, O. 2016, ARA&A, 54, 529, doi: 10.1146/annurev-astro-081915-023441
-
[3]
Brouwer, M. M., Oman, K. A., Valentijn, E. A., et al. 2021, A&A, 650, A113, doi: 10.1051/0004-6361/202040108
-
[4]
Quillen, A. C. 2021, ApJL, 907, L26, doi: 10.3847/2041-8213/abd635
-
[5]
Desmond, H. 2023a, MNRAS, 526, 3342, doi: 10.1093/mnras/stad2762
-
[6]
Desmond, H. 2023b, MNRAS, 521, 1817, doi: 10.1093/mnras/stad597
-
[7]
2025, PhRvD, 111, 103036, doi: 10.1103/PhysRevD.111.103036
Donlon, T., Chakrabarti, S., Vanderwaal, S., et al. 2025, PhRvD, 111, 103036, doi: 10.1103/PhysRevD.111.103036
-
[8]
I., Chakrabarti, S., Widrow, L
Donlon, T. I., Chakrabarti, S., Widrow, L. M., et al. 2024, PhRvD, 110, 023026, doi: 10.1103/PhysRevD.110.023026
-
[9]
Famaey, B., & McGaugh, S. S. 2012, Living Reviews in Relativity, 15, 10, doi: 10.12942/lrr-2012-10
-
[10]
Freundlich, J., Famaey, B., Oria, P.-A., et al. 2022, A&A, 658, A26, doi: 10.1051/0004-6361/202142060 J´ ulio, M. P., Read, J. I., Pawlowski, M. S., et al. 2025, A&A, 704, A330, doi: 10.1051/0004-6361/202557106
-
[11]
Lelli, F., McGaugh, S. S., Schombert, J. M., & Pawlowski, M. S. 2017, ApJ, 836, 152, doi: 10.3847/1538-4357/836/2/152
-
[12]
McGaugh, S. S., Lelli, F., & Schombert, J. M. 2016, Physical Review Letters, 117, 201101, doi: 10.1103/PhysRevLett.117.201101
-
[13]
McMillan, P. J. 2017, MNRAS, 465, 76, doi: 10.1093/mnras/stw2759
-
[14]
Mistele, T., McGaugh, S., Lelli, F., Schombert, J., & Li, P. 2024, Journal of Cosmology and Astroparticle Physics, 2024, 020, doi: 10.1088/1475-7516/2024/04/020
-
[15]
Moran, A., Mingarelli, C. M. F., Van Tilburg, K., & Good, D. 2024, PhRvD, 109, 123015, doi: 10.1103/PhysRevD.109.123015
-
[16]
Paranjape, A., & Sheth, R. K. 2021, MNRAS, 507, 632, doi: 10.1093/mnras/stab2141
-
[17]
F., Ravi, A., Ebadi, R., & Walsworth, R
Phillips, D. F., Ravi, A., Ebadi, R., & Walsworth, R. L. 2021, PhRvL, 126, 141103, doi: 10.1103/PhysRevLett.126.141103
-
[18]
Stiskalek, R., & Desmond, H. 2023, MNRAS, 525, 6130, doi: 10.1093/mnras/stad2675 V˘ ar˘ a¸ steanu, A. A., Jarvis, M. J., Ponomareva, A. A., et al. 2025, MNRAS, 541, 2366, doi: 10.1093/mnras/staf1079
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.