Issue |
A&A
Volume 698, May 2025
|
|
---|---|---|
Article Number | A167 | |
Number of page(s) | 12 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361/202554793 | |
Published online | 13 June 2025 |
The p-Laplacian as a framework for generalizing Newtonian gravity and Milgromian gravitation
1
Universität Bonn, 53115 Bonn, Germany
2
Helmholtz-Institut für Strahlen- und Kernphysik, Universität Bonn, Nussallee 14-16, 53115 Bonn, Germany
3
Astronomical Institute, Charles University, V Holesovickach 2, 18000 Praha, Czech Republic
4
School of Astronomy and Space Science, Nanjing University, Nanjing 210093, People's Republic of China
5
Key Laboratory of Modern Astronomy and Astrophysics (Nanjing University), Ministry of Education, Nanjing 210093, People's Republic of China
⋆ Corresponding author: s6dvsche@uni-bonn.de
Received:
27
March
2025
Accepted:
21
April
2025
Context. The radial acceleration relation (RAR) follows from Milgromian gravitation (MoND) and velocity dispersion data of many dwarf spheroidal galaxies (dSphs) and galaxy clusters have been reported to be in tension with it.
Aims. We consider the generalized Poisson equation (GPE), expressed in terms of the p-Laplacian, which has been applied in electrodynamics, and investigate whether it can address these tensions.
Methods. From the GPE we derive a generalized RAR characterized by the p parameter from the p-Laplacian and a velocity dispersion formula for a Plummer model. We apply these models to Milky Way and Andromeda dSphs and HIFLUGS galaxy clusters and derive a p parameter for each dSph and galaxy cluster. We explore a relation of p to the mass density of the bound system, and alternatively a relation of p to the external field predicted from Newtonian gravity
Results. This ansatz allows the deviations of dSphs and galaxy clusters from the RAR without the need of introducing dark matter. Data points deviate from the Milgromian case, p = 3, with up to 5σ-confidence. Also, we find the model predicts velocity dispersions, each of which lies in the 1σ-range of their corresponding data point allowing the velocity dispersion to be predicted for early-type dwarf satellite galaxies from their baryonic density. The functional relation between the mass density of the bound system and p suggests p to increase with decreasing density. We find for the critical cosmological density p(ρcrit) = 12.27±0.39. This implies significantly different behaviour of gravitation on cosmological scales. Alternatively, the functional relation between p and the external Newtonian gravitational field suggests p to decrease with increasing field strength.
Conclusions. The GPE fits the RAR data of dSphs and galaxy clusters, reproduces the velocity dispersions of the dSphs, gives a prediction for the velocity dispersion of galaxy clusters from their baryonic density and may explain the non-linear behaviour of galaxies in regions beyond the Newtonian regime.
Key words: gravitation / galaxies: clusters: general / galaxies: dwarf / galaxies: groups: general / galaxies: kinematics and dynamics / dark matter
© The Authors 2025
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1. Introduction
Modified Newtonian dynamics (MoND) was introduced by Milgrom (1983). It is a non-relativistic modification of the Newtonian law of gravity, which can largely explain the dynamics in galaxies without the need to introduce dark matter. This constitutes an important research program (Merritt 2020) given the tensions between observations and dark-matter-based models (Kroupa et al. 2023). Still, the acceleration, a0, known as Milgrom's constant needs to be introduced in the theory and must be obtained from data. The value a0=(1.20±0.02)×10−10 m s−2 (Lelli et al. 2017) is used for data analysis in Chapter 3. MoND replicates Newtonian behaviour for large gravitational accelerations of a≫a0, while postulating a logarithmic potential for a≪a0, recently verified through weak lensing by Mistele et al. (2024). If the gravitational acceleration predicted by MoND is plotted against the gravitational acceleration predicted by Newton's law of gravity, gN, a deviation from the line of identity for a<a0 is obtained, so that a>gN. The functional relation, a(gN), is called the radial acceleration relation (RAR). The velocity dispersions of many dwarf spheroidal satellite galaxies (dSphs) and galaxy clusters have been reported to be in tension with the RAR (Safarzadeh & Loeb 2021; Lelli et al. 2017; Sánchez-Salcedo & Hernandez 2007; Angus 2008; McGaugh & Milgrom 2013a, b; Alexander et al. 2017). This tension can be interpreted as implying that MoND might not be a generally valid effective description of non-relativistic gravitation. However, this tension needs to be seen in view of the multiple flavours of dSphs among the data of Lelli et al. (2017) and Safarzadeh & Loeb (2021): those that lie along the RAR, those that should not be on the RAR in MoND due to the external field effect (EFE, predicted by MoND), those that are subject to tidal disruption (a relevant process also in Newtonian dynamics, as was shown by Kroupa 1997 and applied to dSph data by McGaugh & Wolf 2010), and those that might be genuinely offset. The dSphs that fall on the RAR generally (though not always) have the best data and those that are offset have consistently less reliable data. In MoND, one must distinguish between isolated dSphs and those subject to the EFE. The former should follow the RAR but the latter should not. Both of these situations assume dynamical equilibrium, but there is a narrow window in which that holds before the EFE becomes tidally disruptive. Judging from the simulations by Brada & Milgrom (2000), some dSphs should be out of equilibrium and off the RAR (Kroupa 1997). The complexity of the gravitational stellar-dynamics allows for new ansatzes (Sanders 2021; McGaugh et al. 2021). Here, we show that generalizing gravitation in the regime of a<a0 with the help of the p-Laplacian permits a comprehensive accommodation of dSphs and galaxy clusters.
The generalized Poisson equation (GPE) reads (Evans 1982):
Φ(x) is the gravitational potential, such that the acceleration is a(x) = −∇Φ with a=|a|. G is the gravitational constant, ρ(x) the mass density distribution, and x the position vector. p is a real dimensionless parameter with p>1 (Lindqvist 2019a). For p = 2, it follows the Newtonian potential. For p = 3, the logarithmic Milgromian potential is obtained outside of ρ(x). Aside from these two p values, there is no physical meaning of p known in the context of gravitational theory. The data analysed here allows for the possibility that the p-Laplacian approach combines the non-equilibrium dynamics (EFE and tidal disruption) through a varying p value. In electrodynamics, the p-Laplacian with arbitrary p is used to describe electrorheological fluids, which are a special sort of non-linear dielectrics (Růžička 2000). There, p is a material constant and can be a function of the electric field. Interpreting p as the dimension of three-dimensional or higher-dimensional Euclidean space, Milgrom (1997) showed the GPE to be conformally invariant. The analysis presented here does not follow his interpretation, but investigates a functional relation between p and ρ(x), assuming ρ(x) to be isolated, and a possible functional relation between p and an external Newtonian gravitational field, gtot. This analysis considers simple solutions for arbitrary p in spherical symmetry in Sect. 2, applies them to data on dSphs and galaxy clusters in Sect. 3, elaborates on MoND being possibly not fundamental dynamics in Sect. 4, and discusses possible physical implications in Sect. 5.
2. Derivation of the RAR, gravitational potential, and velocity dispersion assuming Plummer models for arbitrary p
2.1. RAR
In this subsection, the RAR is derived for arbitrary p from the GPE, Eq. (1), as a static spherically symmetric problem. This generalized RAR enables to obtain the p parameters in Sect. 3. The GPE simplifies in spherical symmetry with to the expression
The Newtonian acceleration, gN(r), at a distance, r, in a spherical volume, , is given by
Solving the GPE, Eq. (2) by radial integration leads to the p-dependent relation between the fundamental and Newtonian accelerations:
and is shown in Fig. 1. In the Newtonian case p = 2, the RAR becomes the identity. For large p:
Solving Eq. (5) for a(r) yields::
In the cases p = 2 and p = 3, Eq. (7) simplifies to the Newtonian and Milgromian gravitational accelerations:
2.2. The velocity dispersion of the generalized Plummer mass distribution
In this subsection, the velocity dispersion for arbitrary p and a Plummer model is derived. The Plummer model is a simple and proper density distribution to describe dSphs. The velocity dispersion cannot be calculated independently of the density distribution, if p≠3 (for a proof, see Appendix F). We considered the gravitational acceleration for arbitrary p (Eq. (7)) and the Plummer model (p. 81 in Heggie & Hut 2003):
with the Plummer radius, rpl, and stellar mass, Mst. The radial component of the isotropic velocity dispersion, , obtained from the Jeans equation, reads (Chap. 4.8 in Binney & Tremaine 2008):
Inserting Eqs. (7), (10) and (11) into Eq. (12) yields:
after introducing the following definitions:
where 2F1(a,b;c;z) is the hypergeometric function.
The characteristic 1D velocity dispersion, , was obtained.
is the kinetic energy (Chap. 4.8 in Binney & Tremaine 2008) for the Plummer sphere.
Inserting Eq. (15) for σr(r) and Eq. (11) for ρ(r) yields:
with
Converting the integration variable to q(r) and integrating by parts results in:
with
In the last step, the lower limit q(r) of the integral (Eq. (20)) does not contribute: for p>1, limq→0(Λ(q,p)) is finite, limq→0(q3) = 0, , and
. Using l’H
spital's rule for
leads to the final expression for :
In Sect. 3, velocity dispersions for dSphs are calculated with Eq. (30). The Plummer model requires spherical symmetry and the velocity distribution is assumed to be isotropic. In general, this is not likely to be true for dSphs (Kroupa 1997).
2.3. Newtonian (p = 2) and Milgromian (p = 3) limit for the velocity dispersion
In the Newtonian (p = 2) and Milgromian (p = 3) limit, Eq. (30) reduces to the known expressions for the 1D velocity dispersion (cf. Binney & Tremaine 2008 and Milgrom 1994):
This result for p = 3 is only valid in the field limit for a one-dimensional component of isotropic velocity dispersion. Also, it can be derived from other principles, such as scale invariance in the modified-gravity framework (Milgrom 2014).
3. Data analysis
The data of Milky Way and Andromeda satellite dSphs from Lelli et al. (2017) and a selection from the HIghest X-ray FLUx Galaxy Cluster Sample (HIFLUGS) from Li et al. (2023) were used to check whether they hint at a relation between p and the density of the satellite galaxy or the galaxy cluster. To accomplish this, what we interpret to be gN and a from observations was obtained for each dSph and galaxy cluster. Since there is no interest in the radial dependence of gN(r) and a(r), from now on gN and a is written. With the RAR for arbitrary p derived in Chapter 2, p was calculated. From these p parameters, velocity dispersions were obtained for the Milky Way dSphs analysed. All relevant data are shown in Table A.1 for the dSphs and in Table B.1 for the galaxy clusters. Figure 2 shows the data compared to the RAR for arbitrary p derived in Chapter 2.
![]() |
Fig. 2. RAR for selected Milky Way (blue) and Andromeda (orange) satellite dSphs and galaxy clusters (green). Uncertainties are 1σ. The data are taken from Lelli et al. (2017) for the dSphs (Milky Way and Andromeda satellites) and Li et al. (2023) for the galaxy clusters, as is tabulated in Tables A.1 and B.1, respectively. The bold grey lines for different p are the same as in Fig. 1. |
For error propagation, normally distributed and uncorrelated values were assumed. The error Δα of an arbitrary quantity α(β1,…,βn) was calculated from arbitrary quantities, β1,…,βn, , with uncertainties of Δβ1,…,Δβn:
Uncertainties were calculated to two significant digits. If the third significant digit was smaller than 5, uncertainties were rounded down; otherwise, they were rounded up.
3.1. RAR for the data at hand
What we interpret to be gN and a of the dSphs has been inferred by Lelli et al. (2017) from observational quantities. To accomplish this, they collected distances, V-band magnitudes, half-light-radii, and velocity dispersions from various sources. Velocity dispersions were obtained from high-resolution individual-star spectroscopy. We adopted their values, interpreting their data in the following way:
They calculated gbar with Newton's law of gravity and stellar mass, only contributing and gobs from velocity dispersion using an estimator by Wolf et al. (2010). Here, the analysis was limited to their high-quality sample, reducing the impact of tidal effects and neglecting galaxies with ellipticity larger than 0.45 (for details, see Lelli et al. 2017).
From the HIFLUGS clusters, the gas mass, Mgas,500, stellar mass, M*,500, and dynamical mass, Mdyn,500, within the virial radius, r500, were used. Li et al. (2023) obtained the gas mass from integrating surface brightness fits and the gas mass by assuming a correlation to the stellar mass. They inferred the dynamical mass from X-ray and optical data, taking into account various morphological properties, dynamical properties, and the temperature, assuming hydrostatic equilibrium. They took their values for the virial radius from Zhang et al. (2011), who obtained it from the X-ray measured mass distribution, assuming hydrostatic equilibrium.
The same derivation as in Li et al. (2023) was followed here. From the gas mass, Mgas,500, stellar mass, M*,500, and dynamical mass, Mdyn,500, enclosed in the virial radius, r500, for the galaxy clusters (gc), a and gN were computed according to Newton's law of gravity:
with the gravitational constant G=(4.300 917 270 6±0.0000000003) pc/M⊙(km s)2. Uncertainties were not given explicitly by Li et al. (2023). Uncertainties were taken for Mgas,500 from Zhang et al. (2011) and calculated for M*,500 according to Eq. (23) in Li et al. (2023):
Uncertainties were not reconstructed for r500 and Mdyn,500, and so underestimated the uncertainties for gN and a. The uncertainty from the gravitational constant, G, was neglected, since it is so small that it would not contribute significantly to the uncertainties for gN and a.
Figure 2 shows the RAR for the data at hand. Deviations from the p = 3 line can clearly be seen. For Canes Venatici I, the deviation reaches a 5σ significance. The galaxy cluster data show a 5σ deviation at least from the p = 3 limit, subject to the caveat that their uncertainties are underestimated.
3.2. Calculating p from the RAR and density for the data at hand
To calculate p for the data at hand, Eq. (5) was used. Solving Eq. (5) for p yields:
Error propagation according to Eq. (33) demands symmetrical uncertainties, which is not fulfilled the values for a and gN. Uncertainties are therefore symmetrized by taking the larger uncertainty, overestimating the uncertainties. The error of a0 does not significantly influence the uncertainties of p for the dSphs, but for galaxy clusters. This is because not all relevant uncertainties for the galaxy clusters are considered. The baryonic mass of the dSphs was computed from their V-band luminosity, LV, and the stellar mass-to-light-ratio, , given by Lelli et al. (2017):
Lelli et al. (2017) inferred the mass-to-light-ratio as a statistical mean from studies of Fornax and Sculptor dSphs by de Boer et al. (2012). The mean density, , of the dSphs was estimated from their half-light-radius, r1/2, from Lelli et al. (2017):
The mean density, , of the galaxy clusters was estimated as:
The exact shape was considered for neither the dSphs nor the galaxy clusters, which leads to a systematic uncertainty.
Figure 3 shows the density of the analysed dSphs and galaxy clusters against p. The p parameters deviate from p = 3, in many cases with a 1σ significance, in some with a 2σ significance, and for Canes Venatici I with a 5σ significance. The uncertainties on LV and r1/2, and therefore on p of Canes Venatici I, are relatively small. For the galaxy clusters, we refrained from making a claim, since the uncertainties were underestimated. Furthermore, the Python package kafe2 was used to fit an exponential function of the following form and with the parameters ρα and α to the data, with ndf being the number of degrees of freedom:
![]() |
Fig. 3. Relation between density of dSphs (blue and orange) and galaxy clusters (green) and p parameter from GPE, Eq. (1). Uncertainties are 1σ. p was obtained from the RAR, as is shown in Fig. 2. See Tables A.1 and B.1 for values. The red line is an exponential function fit to the data with χ2/ndf = 8.001. |
8 M⊙ pc−3 was used as a starting value for ρα. The large χ2/ndf was due to the small uncertainties for the galaxy clusters that were underestimated. Without the galaxy clusters, χ2/ndf = 3.115 is obtained1. This behaviour is extrapolated to the critical density, ρcrit, of the Universe at redshift z = 0 from the WMAP 7-year cosmology implemented in the astropy.cosmology Python package:
3.3. Infinity Laplacian for weak fields
When p→∞, the p-Laplace operator transitions to a limiting form. The non-linear operator, , tends to the infinity Laplacian (Lindqvist 2015):
where ∇(∇ϕ) represents the Hessian matrix of second derivatives, and the shorthand ∇ϕ·∇(∇ϕ) is equivalent to ∇ϕ⊤(∇2ϕ)∇ϕ. The indices i and j range over the three spatial dimensions, and ∂ denotes partial differentiation. The infinity Laplacian governs the minimization of the maximum gradient magnitude of the potential, ϕ, a problem known as the absolutely minimal Lipschitz extension problem (Camilli et al. 2017). Lu & Wang (2008) prove that when the infinity Laplacian is applied to a positive source function, Δ∞ϕ=f(x), and when this function tends to zero, the infinity Laplace operator suppresses variations in the magnitude of ∇ϕ (Lu & Wang 2008, see their Theorem 3.1). Consequently, the inhomogeneous infinity Laplacian converges to the homogeneous case, i.e., Δ∞ϕ = 0. In the case of matter, the source function is always f(x)≥0. This limiting equation implies that the infinity-Laplace operator smooths out the steepest gradients (Lindqvist 2019b), yielding a gravitational acceleration field with nearly constant magnitude, as anticipated in Eq. (9). This result independently shows that in the weakest gravitational potentials, the GPE predicts that gravity smooths out any sharp variations of the gradient magnitude across the non-linear domain.
3.4. Limitations of this analysis
The present work serves as an invitation for future observational studies, which are within reach instrumentally but which have yet to undergo systematic and thorough investigation. To reduce the uncertainty for the mass-to-light-ratios of the dSphs, we suggest an analysis that considers individual mass-to-light-ratios for the dSphs. This could be achieved within a homogeneous study by one single team. In such an analysis, uncertainties on a and gN might also be reduced. To estimate the density more accurately, we suggest considering the shape of the individual objects analysed. That the half-light-radius is used for the dSphs and the virial radius for the galaxy clusters might complicate comparability. We also suggest in the future an analysis with all uncertainties for the galaxy clusters. These improvements have not been possible here as the corresponding data are not available.
Finally, the present work constitutes a static analysis of the dynamics generated by the GPE, Eq. (1). At present, live simulations of the studied objects are only possible for the p = 2 and p = 3 cases (e.g. using the code Phantom of Ramses, Lüghausen et al. 2015). Simulations of the satellite galaxies orbiting their host galaxy can thus be performed for these values of p. For example, a detailed study for p = 2 by Kroupa (1997) already suggested that even Newtonian solutions without dark matter exist for most of the galaxy's satellites. It will be important to extend such work to the p = 3 case to study if apparent p>3 solutions emerge.
3.5. Calculating the velocity dispersion for the Plummer model from an obtained p
To infer velocity dispersions of the dSphs from the data by Lelli et al. (2017), Eq. (30) was used, identifying (Kroupa 2008) and Mst=MdSphs. Now, the full sample of Milky Way satellites was used and compared to the Lelli et al. (2017) and Safarzadeh & Loeb (2021) data neglecting errors in Fig. 4. The sample of Milky Way satellites chosen by Safarzadeh & Loeb (2021) was used. The latter adopted velocity dispersion data from various astrophysical literature sources. Any calculated velocity dispersion lies within the 1σ regime of the corresponding data point by Lelli et al. (2017) and often also of the corresponding Safarzadeh & Loeb (2021) point. The outliers are Ursa Minor with a 2σ tension and Bootes I as well as Fornax with a 3σ tension between the Lelli et al. (2017) and Safarzadeh & Loeb (2021) points. The ellipticity of Ursa Minor is 0.54, which might not be negligible. Still, for Ursa Major I with an ellipticity of 0.80, there is no tension at all. Bootes I has no peculiar properties; there might be a systematic uncertainty that has not been considered. Fornax is the most massive Milky Way dSph in the sample at hand.
![]() |
Fig. 4. 1D velocity dispersion for selected Milky Way satellite dSphs. Blue are data by Lelli et al. (2017) and orange data by Safarzadeh & Loeb (2021). The black triangles were calculated from the Lelli et al. (2017) data with Eq. (30). |
4. Interpreting GPE as a modification of gravitational acceleration under the influence of external Newtonian gravity
In this section, the following hypothesis is examined: the GPE and MoND are not a modification of fundamental dynamics, but just a modification of acceleration. Drawing an analogy to electrorheological fluids suggests that gravity remains Newtonian on a fundamental level, but the p parameter describes a non-linear response of the ‘induced’ field, a(r), inside the satellite galaxy to an external Newtonian field, ghost(r), by the host galaxy and gN(r) by the satellite itself. This would mean that the Poisson equation still holds true while the gravitational field is modified:
In this picture, gtot takes the role of electric field strength, E, and a the role of electric displacement, D. The p parameter is used to describe viscosity behaviour of electrorheological fluids. An analogue physical process in galaxies is unknown. The theory of electrorheological fluids postulates a functional relation between the p parameter and the external field. Here, we likewise inquire whether an analogue relation might exist for the analysed dSphs. For this task, the p parameter is re-evaluated:
Equation (54) holds thanks to superposition and spherical symmetry. Mhost is the mass of the host galaxy (Mhost = 1×1012 M⊙ for the Milky Way and Mhost = 2×1012 M⊙ for Andromeda). Dhost is the distance of the satellite to its host. All values were again taken from Lelli et al. (2017). Since no uncertainties for Mhost and Dhost were given, the uncertainties for gtot were the same as for gN. The uncertainty on G was still neglected. For this task, all 60 satellites analysed by Lelli et al. (2017) were investigated. A logarithmic fit leads to the following result in Fig. 5:
![]() |
Fig. 5. Relation between p parameter and external field gtot for Milky Way (blue) and Andromeda (orange) satellite dSphs. The dashed line describes Milgrom's constant a0. Uncertainties are not shown to keep the plot legible (see Tables C.1 and D.1 for uncertainties on p). The external field is Newtonian and consists of the field by the host galaxy and the internal field of the satellite. A logarithmic fit function (solid green curve) described by Eq. (56) suggests a functional relation, p(gtot). |
The fit model (Eq. (56)) was chosen to be identical to Eq. (53), if a0=g0 and . Due to large uncertainties, these identities cannot be confidently defended. This would imply that a∝a0. Extrapolating to large gtot is unphysical, since Eq. (56) does not reproduce Newtonian behaviour and allows p≤1, which is mathematically forbidden (Lindqvist 2019a). Still, extrapolating to low accelerations yields a testable prediction. The small χ2/ndf is due to large uncertainties. There seem to be five outliers, including two Milky Way dSphs with p≤1. Still, the uncertainties are so large that they cannot be identified as such.
5. Discussion
This analysis leads to the following conclusions:
-
The gravitational acceleration, a(r), generated by the source distribution can be given by Eq. (7), derived from the GPE, Eq. (1), with p depending on
according to Eq. (43)2.
-
For large p, the Milgromian gravitational acceleration, a, becomes independent of the central mass and distance and equals a0 (Eq. (6)).
-
The analysed data hint at a relation between p and the density of the bound system,
(Eq. (43)). The lower the density, the higher p. For the critical density, p(ρcrit) = 12.27±0.39 is found.
-
Combining the GPE with a Plummer model predicts the observed velocity dispersions in Milky Way dSphs with a less than 1σ deviation without exception (Fig. 4). That is, once the density of the satellite galaxy is known, we can determine the p parameter using Eq. (43), and from Eq. (30) we obtain an estimate for the velocity dispersion.
-
Interpreting the GPE as fundamentally Newtonian, but the gravitational field, a, as a non-linear response to an external field, gtot, requires a relation, p(gtot). The analysed data hints at such a relation (Fig. 5).
Items 2 and 3 imply that for densities of the bound system on cosmological scales, the Milgromian gravitational acceleration, a, becomes non-directional and independent of central mass and distance and equals a0. This increases the gravitational acceleration compared to Newtonian gravity, reducing the need for dark matter. Free particles are not possible for p>3, since does not exist. The exponential fit (Eq. (43)) cannot be extrapolated to densities related to p<2, because it is known that there is Newtonian gravitation (p = 2) for far larger densities than ρ(p = 2) = (2.14±0.78) M⊙ pc−3.
Items 1 and 2 were derived from the model presented in Sect. 2. This model makes many simplifying assumptions (e.g. spherical symmetry) discussed in Sect. 2. Item 3 is suggested by the data presented in Sect. 3. From the data, claims of up to 5σ significance are stated.
Item 4 might solve the problem of MoND not being able to predict all velocity dispersions in Milky Way satellites raised by Safarzadeh & Loeb (2021). It must be considered that the p parameters have been inferred from the observed accelerations, which have been inferred from velocity dispersions by (Lelli et al. 2017). An independent test might be beneficial to check whether the data correctly represent the Milgromian acceleration, a, which is fundamentally independent of the velocity dispersion.
Apart from items 3 and 5, the physical meaning of p is to be explored. p might be explained by drawing the parallel to electrorheological fluids. This has been done for standard MoND by Blanchet (2007). The data suggests a relation p(gtot), which makes a prediction for low accelerations. Still, an analysis with more precise data is desirable. If the analogy to electrorheological fluids were justified, any data analysis would always have to explore the effect of external fields. Furthermore, this raises the question of whether there is a physical process leading to a non-linear response of dSphs to external gravitational fields. The analysis presented here may suggest that the GPE describes an effective coupling between the dynamics of a low-acceleration system to the gravitational field in the non-relativistic limit.
This analysis is yet another hint that there might be a deeper underlying gravitational theory generalizing MoNDian and Newtonian gravity. This justifies further reviewing current cosmological models and gravitational theories.
Acknowledgments
We thank the anonymous referee and Mordehai Milgrom for their constructive comments. P.K. acknowledges support from the DAAD-Eastern European exchange programme at the University of Bonn and Charles University in Prague. E.G. acknowledges the support of the National Natural Science Foundation of China (NSFC) under grants NOs. 12173016, 12041305. E.G. acknowledges the science research grants from the China Manned Space Project with NOs. CMS-CSST-2021-A08, CMS-CSST-2021-A07.
References
- Alexander, S. G., Walentosky, M. J., Messinger, J., et al. 2017, ApJ, 835, 233 [Google Scholar]
- Angus, G. W. 2008, MNRAS, 387, 1481 [Google Scholar]
- Binney, J., & Tremaine, S. 2008, Galactic Dynamics: Second Edition (Princeton: Princeton University Press) [Google Scholar]
- Blanchet, L. 2007, Classical Quantum Gravity, 24, 3529 [Google Scholar]
- Brada, R., & Milgrom, M. 2000, ApJ, 541, 556 [NASA ADS] [CrossRef] [Google Scholar]
- Camilli, F., Capitanelli, R., & Vivaldi, M. A. 2017, Nonlinear Anal., 163, 71 [Google Scholar]
- de Boer, T. J. L., Tolstoy, E., Hill, V., et al. 2012, A&A, 544, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Evans, L. C. 1982, J. Differ. Equations, 45, 356 [Google Scholar]
- Heggie, D., & Hut, P. 2003, The Gravitational Million-Body Problem: A Multidisciplinary Approach to Star Cluster Dynamics (Cambridge, UK: Cambridge University Press) [Google Scholar]
- Kroupa, P. 1997, NA, 2, 139 [Google Scholar]
- Kroupa, P. 2008, Initial Conditions for Star Clusters (Netherlands: Springer), 181 [Google Scholar]
- Kroupa, P., Gjergo, E., Asencio, E., et al. 2023, The Many Tensions with Dark-matter Based Models and Implications on the Nature of the Universe [Google Scholar]
- Lelli, F., McGaugh, S. S., Schombert, J. M., & Pawlowski, M. S. 2017, ApJ, 836, 152 [Google Scholar]
- Li, P., Tian, Y., Júlio, M. P., et al. 2023, A&A, 677, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lindqvist, P. 2015, Notes on the Infinity-Laplace Equation [Google Scholar]
- Lindqvist, P. 2019a, Notes on the Stationary p-Laplace Equation (Cham: Springer) [CrossRef] [Google Scholar]
- Lindqvist, P. 2019b, The Infinity Laplacian (Cham: Springer International Publishing), 69 [Google Scholar]
- Lu, G., & Wang, P. 2008, Commun. Partial Differ. Equations, 33, 1788 [Google Scholar]
- Lüghausen, F., Famaey, B., & Kroupa, P. 2015, CJP, 93, 232 [Google Scholar]
- McGaugh, S., & Milgrom, M. 2013a, ApJ, 766, 22 [NASA ADS] [CrossRef] [Google Scholar]
- McGaugh, S., & Milgrom, M. 2013b, ApJ, 775, 139 [NASA ADS] [CrossRef] [Google Scholar]
- McGaugh, S. S., & Wolf, J. 2010, ApJ, 722, 248 [NASA ADS] [CrossRef] [Google Scholar]
- McGaugh, S. S., Lelli, F., Schombert, J. M., et al. 2021, AJ, 162, 202 [NASA ADS] [CrossRef] [Google Scholar]
- Merritt, D. 2020, A Philosophical Approach to MOND: Assessing theMilgromian Research Program in Cosmology (Cambridge University Press) [CrossRef] [Google Scholar]
- Milgrom, M. 1983, ApJ, 270, 365 [Google Scholar]
- Milgrom, M. 1994, ApJ, 429, 540 [NASA ADS] [CrossRef] [Google Scholar]
- Milgrom, M. 1997, Phys. Rev. E, 56, 1148 [NASA ADS] [CrossRef] [Google Scholar]
- Milgrom, M. 2014, Phys. Rev. D, 89, 024016 [NASA ADS] [CrossRef] [Google Scholar]
- Mistele, T., McGaugh, S., Lelli, F., Schombert, J., & Li, P. 2024, ApJ, 969, L3 [NASA ADS] [CrossRef] [Google Scholar]
- Růžička, M. 2000, Electrorheological Fluids Modeling and Mathematical Theory [Google Scholar]
- Safarzadeh, M., & Loeb, A. 2021, ApJ, 914, L37 [Google Scholar]
- Sánchez-Salcedo, F. J., & Hernandez, X. 2007, ApJ, 667, 878 [CrossRef] [Google Scholar]
- Sanders, R. H. 2021, MNRAS, 507, 803 [NASA ADS] [CrossRef] [Google Scholar]
- Wolf, J., Martinez, G. D., Bullock, J. S., et al. 2010, MNRAS, 406, 1220 [NASA ADS] [Google Scholar]
- Zhang, Y. Y., Andernach, H., Caretta, C. A., et al. 2011, A&A, 526, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A: RAR data
All relevant quantities for the dSphs and galaxy clusters.
Cluster | ![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
p | ![]() |
---|---|---|---|---|---|---|---|---|
A0085 | 1217 | 6.67(32) | 4.39(13) | 8.66 | -11.175(20) | -10.089 | 8.46(30) | 1.652(74) |
A0262 | 755 | 1.08(12) | 1.47(10) | 1.35 | -11.523(43) | -10.481 | 3.858(80) | 1.19(12) |
A0496 | 967 | 2.79(12) | 2.60(7) | 4.65 | -11.342(17) | -10.159 | 6.97(17) | 1.413(56) |
A0576 | 869 | 2.00(71) | 2.13(45) | 4.00 | -11.39(14) | -10.132 | 7.96(69) | 1.413(45) |
A1795 | 1085 | 4.95(14) | 3.68(6) | 4.54 | -11.201(11) | -10.270 | 4.670(64) | 1.744(46) |
A2029 | 1247 | 8.24(29) | 4.99(11) | 8.38 | -11.106(14) | -10.124 | 6.83(19) | 1.887(63) |
A2142 | 1371 | 13.40(80) | 6.68(24) | 8.17 | -10.982(25) | -10.218 | 4.57(10) | 2.29(13) |
A2589 | 837 | 1.77(12) | 1.98(8) | 3.11 | -11.407(27) | -10.208 | 6.17(13) | 1.406(86) |
A3158 | 1013 | 3.75(33) | 3.11(16) | 4.56 | -11.258(35) | -10.208 | 5.66(15) | 1.64(13) |
A3571 | 1133 | 5.16(28) | 3.77(12) | 7.47 | -11.221(22) | -10.091 | 8.64(31) | 1.595(81) |
Notes. The columns are, from left to right: the name of the cluster, the virial radius, r500, the gas mass, Mgas,500, the stellar mass, M*,500, the dynamical mass, Mdyn,500, Newtonian acceleration, gN, observationally deduced acceleration, a, the p parameter from p-Laplacian, p (Eq. (1)), and the mean density, . We adopted the first four quantities from Li et al. (2023) and calculated from them gN, a, p, and ρ. Uncertainties for Mgas,500 were taken from Zhang et al. (2011). Uncertainties for M*,500 have been reconstructed via Eq. (38).
Appendix B: Velocity dispersion
Velocity dispersion for Milky Way dSphs.
Appendix C: External field Milky Way
Data on how Milky Way affects its dSphs via gravity.
Appendix D: External field Andromeda
Relevant data on testing for the effect of the Andromeda gravitational field on Andromeda dSphs.
Appendix E: External field additional plots
![]() |
Fig. E.1. Relation between p parameter and external field for Milky Way (blue) and Andromeda (orange) satellite dSphs. Uncertainties are not shown to keep the plot legible. The dashed line describes Milgrom's constant, a0. The external field is Newtonian and consists of either the field by the host galaxy E.1a or the internal field of the satellite E.1b. The dependence of p on the host field is more obvious. |
Appendix F: Derivation of velocity dispersion independently of density profile
In the following, the velocity dispersion for arbitrary p neglecting the EFE is derived. Furthermore, it is shown that, for arbitrary p, numerical values for the velocity dispersion can only be calculated by specifying a density profile. Milgrom (1994) derived the velocity dispersion for p = 3. This derivation is analogous to his. Milgrom considered the quantity, Q, integrating over phase-space:
mk are the masses of arbitrary particles (e.g. stars). fk(r,v,t) is a phase-space-distribution function. Milgrom showed that its time derivative, , can be written as:
Here, the finite total mass, M, of the system and the momentary, mass-weighted, 3-D, mean-square velocity, 〈v2〉(t), is introduced. The GPE is used:
As opposed to Milgrom, μ(x) is here defined in terms of a fixed p>1. Milgrom inserted the density distribution, ρ, described by Eq. (F.3) into Eq. (F.2). Furthermore, F(y) was introduced, such that
Since in the case for arbitrary p, μ(x) is defined by Eq. (F.4):
Milgrom showed that can be written as
Milgrom solved Eq. (F.8) for 〈v2〉(t) and took the long-time average, . Doing so for arbitrary p:
We assume the density distribution to have a sphere with a radius, R, as compact support and ∇Φ at R for arbitrary p to behave like the gravitational acceleration produced by a point mass:
Equation F.9 was evaluated at R:
Equation F.11 was evaluated at R:
Equation F.12 was evaluated at R:
I4 depends on the density distribution. It vanishes only for p = 3 (Milgromian case). This explains why the velocity dispersion is independent of the density distribution only for Milgromian gravity.
All Tables
Relevant data on testing for the effect of the Andromeda gravitational field on Andromeda dSphs.
All Figures
![]() |
Fig. 1. RAR for arbitrary p. The bold grey lines indicate the RARs for different p (Eq. (5)). |
In the text |
![]() |
Fig. 2. RAR for selected Milky Way (blue) and Andromeda (orange) satellite dSphs and galaxy clusters (green). Uncertainties are 1σ. The data are taken from Lelli et al. (2017) for the dSphs (Milky Way and Andromeda satellites) and Li et al. (2023) for the galaxy clusters, as is tabulated in Tables A.1 and B.1, respectively. The bold grey lines for different p are the same as in Fig. 1. |
In the text |
![]() |
Fig. 3. Relation between density of dSphs (blue and orange) and galaxy clusters (green) and p parameter from GPE, Eq. (1). Uncertainties are 1σ. p was obtained from the RAR, as is shown in Fig. 2. See Tables A.1 and B.1 for values. The red line is an exponential function fit to the data with χ2/ndf = 8.001. |
In the text |
![]() |
Fig. 4. 1D velocity dispersion for selected Milky Way satellite dSphs. Blue are data by Lelli et al. (2017) and orange data by Safarzadeh & Loeb (2021). The black triangles were calculated from the Lelli et al. (2017) data with Eq. (30). |
In the text |
![]() |
Fig. 5. Relation between p parameter and external field gtot for Milky Way (blue) and Andromeda (orange) satellite dSphs. The dashed line describes Milgrom's constant a0. Uncertainties are not shown to keep the plot legible (see Tables C.1 and D.1 for uncertainties on p). The external field is Newtonian and consists of the field by the host galaxy and the internal field of the satellite. A logarithmic fit function (solid green curve) described by Eq. (56) suggests a functional relation, p(gtot). |
In the text |
![]() |
Fig. E.1. Relation between p parameter and external field for Milky Way (blue) and Andromeda (orange) satellite dSphs. Uncertainties are not shown to keep the plot legible. The dashed line describes Milgrom's constant, a0. The external field is Newtonian and consists of either the field by the host galaxy E.1a or the internal field of the satellite E.1b. The dependence of p on the host field is more obvious. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.