REVIEW 4 major objections 4 minor 34 references
Significance of Gravitational Nonlinearities on the Dynamics of Disk Galaxies
T0 review · 4 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read The paper claims that general relativity's nonlinear self-interaction makes baryonic matter alone sufficient to explain the tight relation between galaxies' baryonic mass and observed acceleration, and that the MOND acceleration scale…
desk verdict A bold, honest attempt to explain the radial acceleration relation with GR self-interaction, but the load-bearing disk force law is assumed rather than derived, and the claimed a0 emergence inherits that assumption. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is gravitational field self-interaction, the $n>0$ terms in the expanded Einstein–Hilbert Lagrangian $\mathcal{L}=\sum_n (16\pi G M)^{n/2}[\phi^n(\partial\phi\partial\phi-(16\pi G M)^{1/2}\phi T)]$. By analogy with QCD flux tubes, self-interaction makes the gravitational field lines collapse into the disk plane for a flattened mass distribution, giving a logarithmic potential $\Phi_d(r)=G' M_d^{\mathrm{enc}}(r)\ln r$ and a $1/r$ force in the disk, while a spherical bulge restores the Newtonian $1/r^2$ force inside the transition radius $r_t$. The two-dimensional coupling $G'$ is fixed to $G/r_t$ by continuity at $r_t$. Two complementary galaxy models sample observed parameter correlations among bulge effective radius, Sersic index, disk scale length, and bulge/disk masses, so the resulting $g_{\rm SI}$–$g_{\rm N}$ correlation is parameter-free.
What would settle it
Take a sample of disk galaxies with independent baryonic mass maps and compute $g_{\rm SI}(r)=G M_b^{\mathrm{enc}}(r)/r^2+(G/r_t)M_d^{\mathrm{enc}}(r)/r$ with $r_t$ determined from the bulge–disk force balance, then compare with observed accelerations from rotation curves; if the residuals exceed the observational uncertainties or correlate with galaxy parameters, the central claim fails.
Extended reading notes
Core claim
The central claim is that the nonlinear terms in the Einstein–Hilbert Lagrangian, which make the gravitational field self-interact, grow important at galactic scales and change the effective force law in a disk galaxy from $1/r^2$ to $1/r$ in the disk-dominated region. The transition happens at a radius $r_t$ where the bulge and disk forces balance, and matching the effective two-dimensional coupling $G'=G/r_t$ leaves the models without adjustable parameters. Computed with this rule, the acceleration including self-interaction, $g_{\rm SI}$, plotted against the Newtonian baryonic acceleration $g_{\rm N}$, reproduces the empirical relation between observed and baryonic accelerations; the direct lattice calculation already does so for bulge-less galaxies, and the two sampling models extend this across the observed range of disk morphologies. The acceleration at the transition radius peaks at $1.25\pm 0.06\times 10^{-10}\,\mathrm{m/s^2}$, consistent with the MOND scale $a_0\approx 1.2\times 10^{-10}\,\mathrm{m/s^2}$. Thus, in the paper's own terms, baryonic matter alone suffices to explain disk galaxy dynamics once general relativity's nonlinearity is included.
Load-bearing premise
The load-bearing premise is that gravitational field lines, once distorted by the dense central mass, remain confined to the disk plane in real disk galaxies, so the disk force truly follows $1/r$ with $G'=G/r_t$; the paper explicitly flags this as the assumption that extending the one-dimensional trapping result to a two-dimensional disk does not compromise the trapping of the field.
Editorial extensions
If this is right
- The observed tight relation between baryonic mass and acceleration in disk galaxies stops being evidence for dark matter; it becomes a direct consequence of general relativity's nonlinearity.
- The MOND acceleration scale $a_0$ is not a new constant of nature but the dynamically selected acceleration at the bulge–disk transition, so its value can vary with galaxy morphology and mass distribution.
- Rotation curves should flatten naturally in disk-dominated regions because the GR self-interaction force $1/r$ exceeds the Newtonian $1/r^2$ expectation increasingly with radius.
- Dark-matter-free baryonic models with the $1/r$ force law and $G'=G/r_t$ predict the small intrinsic scatter of the empirical relation, since the residual depends only weakly on galaxy parameters.
Reading between the lines
- If the mechanism is right, gravitational lensing reconstructions of disk galaxies should show the extra apparent mass following the $1/r$ potential rather than a spherical dark halo, a signature testable with current lensing data.
- Because $G'=G/r_t$ depends on galaxy structure, the transition acceleration should shift systematically with bulge-to-disk ratio; MOND predicts a strictly universal $a_0$, so this shift is a distinguishing observation.
- The same field self-interaction should act in galaxy clusters, but their less flattened geometry predicts a weaker or differently shaped correction; cluster dynamics are therefore a natural place to look for the effect's limits.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes that nonlinear self-interaction of the gravitational field in General Relativity can explain the baryonic mass--acceleration relation (RAR) of disk galaxies reported by McGaugh et al. 2016, without dark matter or modified gravity. The central mechanism is that, for an axisymmetric disk, gravitational field self-interaction traps the field lines in two dimensions, producing a logarithmic potential and a 1/r force outside a transition radius rt, while a spherical bulge retains the Newtonian 1/r^2 force. The paper uses direct lattice calculations for a few galaxies, then constructs two phenomenological disk models (Model 1 with uniform sampling, Model 2 with random sampling) using a Sersic bulge plus exponential disk, with G' = G/rt fixed by continuity at rt. It reports that the resulting g_SI versus g_N relation reproduces the MLS2016 correlation and that the acceleration at rt peaks at about 1.2e-10 m/s^2, which it identifies with the MOND scale a0.
Significance. If the central premise were established, the paper would be highly significant: it would offer a baryon-only, dark-matter-free and MOND-free explanation of the RAR, and would provide a dynamical origin for the MOND acceleration scale. The manuscript includes useful cross-checks: lattice calculations that recover known free-field potentials and the Cornell potential, a mean-field disk calculation, and quantitative fits to the MLS2016 data with comparison of residuals. However, the significance is entirely contingent on the assumed 1/r gravitational force law in disks and on the transition prescription. Since that force law is not derived from the Einstein-Hilbert Lagrangian for a realistic galaxy and is explicitly acknowledged in the text to be an assumption, the paper does not currently demonstrate its central claim.
major comments (4)
- [Sec. 4; Sec. 5; Appendix C] The load-bearing input of the models is the assumed force law: Newtonian 1/r^2 inside rt and 1/r outside rt, with effective coupling G' = G/rt (Eqs. 9, 12--14). This force law is not derived from Eq. (1) for a disk mass distribution. The manuscript itself concedes in Sec. 4 that "Extending the one-dimensional result to the two-dimensional disk case assumes that the spread of the mass within the disk area does not compromise the trapping of the field in two dimensions," and Appendix C states that rt (equivalently G') "cannot presently be analytically calculated from first principles" and is "assessed phenomenologically." Since every downstream result, including the reproduction of the MLS2016 relation and the claimed emergent a0, is an output of this imposed force law, the paper does not independently support its central claim that GR nonlinearities make baryonic matter alone sufficient.
- [Sec. 6.2; Eq. (12)] The claimed "dynamically emerging" acceleration scale is largely built into the model definition. In Model 2, rt is defined as the radius where the bulge and disk accelerations are equal, and with G' = G/rt this reduces to M_b^enc(rt) = M_d^enc(rt) (Eq. 13). The acceleration at rt, whose distribution is shown in Fig. 4, is therefore determined by the chosen transition condition, the mass profiles, and the empirical scaling relations Eqs. (3)--(5). It is not an independent prediction from the GR calculation. The statement that a0 can be explained as the acceleration where "the disk mass overtakes the bulge mass" overstates what the model actually derives.
- [Sec. 5; Eq. (10)] Despite the abstract's claim of "parameter-free galactic models," the construction in Eqs. (6)--(14) contains several free or adjustable choices: the transition radius rt (2Re in Model 1; force equality in Model 2), the effective coupling G', and the width of the Fermi-Dirac smoothing function D(r). G' is fixed by continuity, G' = G/rt, rather than predicted from the Lagrangian. The statement in Sec. 5 that "there are no adjustable parameters" is therefore misleading: the parameters are not fitted to the MLS2016 data, but they are also not derived from first principles, as Appendix C acknowledges.
- [Sec. 4; Fig. 1 (top)] The only calculation based directly on Eq. (1) agrees with the observed RAR for bulge-less Hubble types 5 and 6, while types 3 and 4 overestimate g_SI and "lie on the edge of the observed distribution." Because many observed disk galaxies contain bulges, the direct GR-based evidence covers only a subset of the relevant population. The agreement for the full population therefore rests entirely on Models 1 and 2, which impose the aforementioned force law and transition prescription rather than deriving them.
minor comments (4)
- [Eq. (1)] The notation in Eq. (1) is compressed: the coefficients an, the sign conventions, and the meaning of the bracket shorthand are not fully specified, which makes the central equation difficult to audit.
- [Fig. 5] The figure caption lists the Kendall coefficient as ck = 0.207, while the text of Sec. 6.3 reports ck = -0.207. The sign inconsistency should be corrected.
- [Sec. 2] The statement that the ratio GM/L "becomes 10^-2 for galactic systems" is given without explicit unit conventions or a sample calculation. In natural units the numerical value is not transparent, and a reader cannot verify the claimed threshold without additional definitions.
- [Sec. 6.1] The fit of the Model 2 simulation to the MLS2016 functional form, Eq. (2), shows a compatible g†, but the paper should clarify whether the fit is treated as a model comparison or as a phenomenological diagnostic, given that the same functional form is used for both the data and the simulation.
Circularity Check
The claimed 'predictions' are built into the model: the 1/r disk force is assumed (Sec. 4), G' = G/rt is imposed by continuity (Sec. 5.1), and the 'emergent' a0 is the acceleration at the defined transition radius (Sec. 6.2).
-
ansatz smuggled in via citation
[Section 4, direct calculations (paragraph after two-body lattice result)]
"Extending the one-dimensional result to the two-dimensional disk case of galaxies assumes that the spread of the mass within the disk area does not compromise the trapping of the field in two dimensions."
The 1/r disk force is the load-bearing input: it produces the outer-region gSI proportional to sqrt(gN) behavior that matches the MLS2016 relation. The paper does not derive this force law from the Einstein-Hilbert action for a disk; it imports it from prior same-author work (Deur 2009, 2017; Deur 2020) built on a QCD flux-tube analogy and point-source lattice calculations. The agreement with the observed radial acceleration relation is therefore a consistency check of the assumed force law, not a derivation of it from first principles.
-
self definitional
[Sec. 5.1 (before Eq. 10); Sec. 6.2]
"G′ is determined by requiring the accelerations to match at rt: gSI,r<rt(rt) = gSI,r>rt(rt). Thus, G′ = G/rt, by construction."
G' sets the normalization of the 1/r force; fixing it by continuity at the modeler-chosen transition radius means the transition scale is an input, not an output. In Model 2, rt is itself defined by the force-equality condition (Eq. 12), and Sec. 6.2 identifies the 'dynamically emerging' acceleration a(rt) with MOND's a0. Thus a0 is the acceleration at a radius whose definition and continuity condition were chosen to produce a single characteristic scale; its peak value inherits the empirical scaling relations (3)-(5), so it is not an independent prediction.
1 more flagged steps
-
fitted input called prediction
[Appendix C]
"Therefore, even for infinitely thin disks, rt—or equivalently G′—is non-universal and cannot presently be analytically calculated from first principles. It can be obtained from numerical calculations such as those in Refs. (Deur 2009, 2017), or assessed phenomenologically as done in this article."
The key scale rt/G'—which controls the 1/r-force normalization and hence both the RAR reproduction and the value of a0—is 'assessed phenomenologically' rather than derived. The numerical references are same-author works and are for pointlike/two-body sources. The abstract's 'parameter-free galactic models' and 'direct calculations based on General Relativity's Lagrangian' are therefore overstated: the central scale is supplied from phenomenology, and the claimed emergence of a0 is downstream of that phenomenological input.
full rationale
The paper contains genuine calculations: lattice evaluation of the GR Lagrangian for point sources, and a comparison of the resulting gSI-gN curves for bulge-less galaxies to MLS2016. Those direct results for Hubble types 5 and 6 agree with the data. However, the generalization to disk galaxies—the step that produces the full radial acceleration relation—is an explicit assumption (Sec. 4), and the transition scale/coupling is imposed by continuity or 'assessed phenomenologically' (Sec. 5.1, App. C). The emergent a0 is identified with the acceleration at the defined transition radius (Sec. 6.2), so that prediction reduces by construction to the model's chosen transition condition and to the empirical scaling relations used to populate the galaxy sample. Because these load-bearing elements are assumed or imported from prior self-citations rather than independently derived, the central claims are partially circular. The paper's own text candidly states the key limitation, which supports this reading. The score is 7 rather than higher because there is independent content in the direct lattice calculations and the galaxy scaling relations are external to MLS2016, but the central 'parameter-free' predictions are substantially shaped by the model's own definitions.
Assumptions & free parameters
free parameters (3)
- G' (effective 2D gravitational coupling) =
G/rt per galaxy
- Transition radius rt =
Model 1: rt = 2Re; Model 2: solution of M_b(rt) = M_d(rt)
- Transition width of the Fermi-Dirac smoothing function =
rt/2 (Model 1)
assumptions (5)
- standard math Einstein-Hilbert Lagrangian can be expanded as in Eq. (1), with the Fierz-Pauli term as the n=0 part.
- domain assumption GR self-interaction becomes important when GM reaches a fraction of the system size L, around 10^-3 to 10^-2.
- ad hoc to paper Field-line trapping confines gravity to two dimensions in a disk, yielding a logarithmic potential and a 1/r force.
- domain assumption Spherical symmetry cancels self-interaction, so the bulge remains Newtonian with a 1/r^2 force.
- domain assumption The empirical scaling relations in Eqs. (3)-(5) define a representative galaxy parameter space.
Cite this review
Pith. "Pith review of Significance of Gravitational Nonlinearities on the Dynamics of Disk Galaxies." pith.science (2026). https://pith.science/paper/QDUMAVB6
@misc{pith2026190900095,
author = {Pith},
title = {Pith review of: Significance of Gravitational Nonlinearities on the Dynamics of Disk Galaxies},
year = {2026},
howpublished = {\url{https://pith.science/paper/QDUMAVB6}},
note = {Machine review of arXiv:1909.00095}
}
read the original abstract
The discrepancy between the visible mass in galaxies or galaxy clusters, and that inferred from their dynamics is well known. The prevailing solution to this problem is dark matter. Here we show that a different approach, one that conforms to both the current Standard Model of Particle Physics and General Relativity, explains the recently observed tight correlation between the galactic baryonic mass and its observed acceleration. Using direct calculations based on General Relativity's Lagrangian, and parameter-free galactic models, we show that the nonlinear effects of General Relativity make baryonic matter alone sufficient to explain this observation.
Figures
Figures from the paper (6 more)
Reference graph
Works this paper leans on
-
[1]
, " * write output.state after.block = add.period write newline
ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...
-
[2]
write newline
" write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...
-
[3]
- [1] #1 = = ^ ^ ^ .\!\!^ d .\!\!^ h .\!\!^ m .\!\!^ s .\!\!^ @mss
thebibliography [1] 20pt to REFERENCES 6pt =0pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key on reference command Each re...
2019
-
[4]
Brodsky, S. J., de Teramond, G. F., & Deur, A. 2010, Phys.\ Rev.\ D, 81, 096010. arXiv:1002.3948
arXiv 2010
-
[5]
Brodsky, S. J., de Teramond, G. F., Dosch, H. G., & Erlich, J. 2015, Phys.\ Rept., 584, 1, doi:10.1016/j.physrep.2015.05.001
-
[6]
1998, Nucl.\ Phys.\ B, 521, 219, doi:10.1016/S0550-3213(98)00215-6
Buchmuller, W., & Jakovac, A. 1998, Nucl.\ Phys.\ B, 521, 219, doi:10.1016/S0550-3213(98)00215-6
-
[7]
1993, , 265, 1013, doi:10.1093/mnras/265.4.1013
Caon , N., Capaccioli , M., & D'Onofrio , M. 1993, , 265, 1013, doi:10.1093/mnras/265.4.1013
-
[8]
2006, The Millenium Prize Problems (American Mathematical Society)
Carlson, J., Jaffe, A., & Wiles, A. 2006, The Millenium Prize Problems (American Mathematical Society)
work page 2006
Show all 34 references
-
[9]
2009, Physics Letters B, 676, 21, 10.1016/j.physletb.2009.04.060
Deur, A. 2009, Physics Letters B, 676, 21, 10.1016/j.physletb.2009.04.060
2009 doi
-
[10]
2014, Monthly Notices of the Royal Astronomical Society, 438, 1535, doi:10.1093/mnras/stt2293
---. 2014, Monthly Notices of the Royal Astronomical Society, 438, 1535, doi:10.1093/mnras/stt2293
2014 doi
- [11]
-
[12]
2019, European Physics Journal C, in press, doi:10.1140/epjc/s10052-019-7393-0
---. 2019, European Physics Journal C, in press, doi:10.1140/epjc/s10052-019-7393-0
2019 doi
-
[13]
---. 2020. arXiv:2004.05905
2020 arXiv
-
[14]
J., & de Teramond, G
Deur, A., Brodsky, S. J., & de Teramond, G. F. 2016, Prog. Part. Nucl. Phys., 90, 1, doi:10.1016/j.ppnp.2016.04.003
2016 doi
-
[15]
D., Hoyer, P., & Jrvinen, M
Dietrich, D. D., Hoyer, P., & Jrvinen, M. 2013, Phys.\ Rev.\ D, 87, 065021
2013
-
[16]
1975, Phys.\ Rev.\ Lett., 34, 369, doi:10.1103/PhysRevLett.34.369
Eichten, E., Gottfried, K., Kinoshita, T., et al. 1975, Phys.\ Rev.\ Lett., 34, 369, doi:10.1103/PhysRevLett.34.369
1975 doi
-
[17]
1938, Annals of Mathematics, 39, 65
Einstein, A., Infeld, L., & Hoffmann, B. 1938, Annals of Mathematics, 39, 65. http://www.jstor.org/stable/1968714
1938
-
[18]
1939, Proc
Fierz, M., & Pauli, W. 1939, Proc. Roy. Soc. Lond. A, A173, 211, doi/10.1098/rspa.1939.0140
1939
-
[19]
2011, J.\ Nonlin.\ Math.\ Phys., 18, 291, doi:10.1142/S1402925111001441
Frasca, M. 2011, J.\ Nonlin.\ Math.\ Phys., 18, 291, doi:10.1142/S1402925111001441
2011 doi
-
[20]
W., & Worley , C
Graham , A. W., & Worley , C. C. 2008, , 388, 1708, doi:10.1111/j.1365-2966.2008.13506.x
2008
-
[21]
G., Wadadekar, Y., & Kembhavi, A
Khosroshahi, H. G., Wadadekar, Y., & Kembhavi, A. 2000, The Astrophysical Journal, 533, 162, 10.1086/308654
2000 doi
-
[22]
S., Schombert, J
Lelli, F., McGaugh, S. S., Schombert, J. M., & Pawlowski, M. S. 2017, The Astrophysical Journal, 836, 152, 10.3847/1538-4357/836/2/152
2017 doi
-
[23]
D., Ben\' tez-Llambay, A., Schaller, M., et al
Ludlow, A. D., Ben\' tez-Llambay, A., Schaller, M., et al. 2017, Phys. Rev. Lett., 118, 161103, doi:10.1103/PhysRevLett.118.161103
2017 doi
-
[24]
S., Lelli, F., & Schombert, J
McGaugh, S. S., Lelli, F., & Schombert, J. M. 2016, Phys. Rev. Lett., 117, 201101, doi:10.1103/PhysRevLett.117.201101
2016 doi
-
[25]
M \'e ndez-Abreu , J., Aguerri , J. A. L., Corsini , E. M., & Simonneau , E. 2008, , 487, 555, 10.1051/0004-6361:20078089e
2008 doi
-
[26]
1983, Astrophysical Journal, 270, 365, doi:10.1086/161130
Milgrom, M. 1983, Astrophysical Journal, 270, 365, doi:10.1086/161130
1983 doi
-
[27]
2008, Int.\ J.\ Mod.\ Phys.\ D, 17, 367, doi:10.1142/S0218271808012085
Padmanabhan, T. 2008, Int.\ J.\ Mod.\ Phys.\ D, 17, 367, doi:10.1142/S0218271808012085
2008 doi
-
[28]
1959, Nuovo Cim., 14, 951, doi:10.1007/BF02728177
Regge, T. 1959, Nuovo Cim., 14, 951, doi:10.1007/BF02728177
1959 doi
-
[29]
1974, in IC/74/55
Salam, A. 1974, in IC/74/55. http://inspirehep.net/record/90219
1974
-
[30]
S \'e rsic , J. L. 1963, Boletin de la Asociacion Argentina de Astronomia La Plata Argentina, 6, 41
1963
-
[31]
2015, Publications of the Astronomical Society of Japan, 68, doi:10.1093/pasj/psv103
Sofue, Y. 2015, Publications of the Astronomical Society of Japan, 68, doi:10.1093/pasj/psv103
2015 doi
-
[32]
B., & Fisher , J
Tully , R. B., & Fisher , J. R. 1977, , 54, 661
1977
-
[33]
2013, Einstein Gravity in a Nutshell (Princeton University Press)
Zee, A. 2013, Einstein Gravity in a Nutshell (Princeton University Press)
2013
-
[34]
WEŎXQ (E Aѹ wr[ЅV 6 ͝L䗐ϙ3g(
eqnarray [ - 16 GM T ]= 1 2 ^ _ _ ^ - 1 2 ^ _ ^ _ _ ^ - ^ _ _ ^ + ^ _ ^ ^ _ - (16 GM) ^ T_ . eq:Fierz-Pauli Lagrangian eqnarray While Eq. ( eq:EH ) is often used to study quantum gravity---with questions raised regarding its applicability in that context, see e.g. Padmanabhan:...
2020
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.