- Cliniques Universitaires Saint-Luc, Université Catholique de Louvain, Brussels, Belgium
A novel axially symmetric metric is proposed to solve the Einstein field equations. This provides an analytical solution within the matter in the equatorial plane for any galaxy density profile. The solution predicts the observed increase in rotational velocity up to the edge of the galaxy’s bulge. However, beyond the bulge, the rotational velocity remains constant, which contradicts the observed peak curves. The existence of the Universe is then considered by approximating the gravitational fields within the galaxy as the sum of those generated by the galaxy and the Universe. The resulting solution explicitly includes a Universe frame-dragging term, aligning with the sixth version of Mach’s principle proposed by Bondi and Samuel: “inertial mass is affected by the global distribution of matter”. Neglecting the presence of the Universe, stars only have a relative rotation to the bulge, and their rotational velocities monotonically increase with the radial distance r to balance the increasing mass contained in distances < r. At larger distances, the bulge’s attraction and its frame-dragging effect decrease, resulting in a constant rotational velocity. When the Universe is considered, stars also have a relative rotation to the non-rotating Universe and experience an additional centrifugal force at any distance from the bulge. This component induces a decrease in rotational velocity as the gravitational influence of the bulge diminishes with r. This model predicts the observed rotational velocity curves for the galaxies M31, M101, and M81 without requiring any dark matter halo or adjustable parameters. This success substantiates Mach’s idea as an alternative to the dark matter halo theory.
1 Introduction
Most galaxies are dominated by a central mass according to their observed luminous material. However, the observed rotational velocity in these galaxies increases up to a radial distance significantly greater than what Newtonian dynamics predicts. Solving this discrepancy has been the subject of numerous studies.
Slight modifications to the law of gravitation, such as modified Newtonian dynamics (MOND), were proposed to address this issue (Milgrom, 2002). However, the most accepted theory is the presence of a dark matter halo, with several proposed profiles (Navarro, 1996; Merritt et al., 2005; Merritt et al., 2006).
Numerous axially symmetric solutions to the Einstein field equations have been derived and studied (Bronnikov et al., 2020). Recently, it has been demonstrated that the full resolution of the Einstein field equations using the Weyl–Lewis–Papapetrou coordinates could explain the flatness of the rotational velocity curve (Cooperstock and Tieu, 2007, Beordo et al., 2024). A simplified Mach’s principle added to Newtonian dynamics has also been proposed as an explanation for the flatness of the velocity curve (Darabi, 2013).
In a lecture on the origin of the general theory of relativity, Albert Einstein mentioned that he was already familiar with the intriguing idea of Ernst Mach that the inertia of anybody is the result of the interaction of that body with the rest of the bodies in the Universe (though this formulation by Einstein is not found in Ernst Mach’s work (Mach 1901)). Einstein further said: “This idea fascinated me; but it did not provide a basis for a new theory” (Einstein 1933). Although general relativity is not directly based on this idea, Einstein was enthusiastic about the existence of the Lense–Thirring effect, widely considered as a Machian manifestation of the general relativity (Lense and Thirring, 1918; Bondi and Samuel, 1997). However, some exact solutions, such as the Gödel rotating Universe (Gödel, 1949) and the completely empty hyperboloid Universe found by De Sitter in 1917 which was of particular concern to Einstein (Janssen, 2016), are in complete contradiction to Mach’s principle.
The relationship between Mach’s principle and general relativity has been a controversial topic, intensely discussed throughout the 20th century (Sciama, 1953; Rindler, 1994; Bondi and Samuel, 1997; Narlikar, 2003). At the Tübingen conference on Mach’s principle in 1995, an exit poll showed that three participants believed general relativity was Machian, while 21 did not; 14 participants believed general relativity using appropriate boundary conditions was Machian, while seven did not (Barbour et al., 1995). Several attempts have been made to integrate Mach’s principle with general relativity (Brans, 1962; Licata et al., 2016), Einstein himself attempted to reformulate it to fully agree with this principle.
In recent decades, researchers have focused more on searching for a perturbing body, i.e., dark matter, rather than on improving or developing theories. A winning strategy to explain the anomaly of Uranus’ orbit, but which failed to explain the perihelion anomaly of Mercury via the hypothetical planet Vulcan (Baum and Sheehan 2013).
This study proposes an axially symmetric metric that enables an analytical solution to the Einstein field equations and investigates whether a solution explaining the observed galactic rotational velocity profiles can emerge by adding Machian properties.
2 Method
2.1 General solution derivation
We consider the following stationary axially symmetric diagonal metric in the coordinate system [t, r, ϕ, z]:
where λ, β, and γ depend only on r and z.
All the differentials and tensors were computed using wxMaxima 22.04.0 (https://wxMaxima-developers.github.io/wxmaxima/) and analytical integrations with https://www.integral-calculator.com/.
The non-null components of the mixed Einstein tensor for this metric are as follows (see the cyl_metrics.wxmx file, which in this paper ranges from zero to three rather than from one to four in Maxima):
With regard to the length of the equations, we maintained the compact Maxima notation in Eqs 2–6, i.e., the coordinates in the subscript refer to the differential against these coordinates.
The r-collapse is prevented by imposing:
where
Notably, by using the metric in Eq. 1, we neglected the handling of the components
The introduction of an integrating factor is a well-known method that involves multiplying nonlinear differential equations by a factor to facilitate their solution (Hermann et al., 2016). To facilitate the solutions of Eqs. 7, 8, we use a similar strategy, except for the introduction here of an additive integrating factor
To overcome the nonlinearity of Eq. 4, we focus on systems modeled by:
As explained at the end of this section, in the solution,
Eq. 11 can be directly integrated into:
where A is an arbitrary function of r. By introducing Eqs. 10, 12 into Eq. 9, we derive:
By inserting Eqs. 9, 10, 13 into Eq. 7 and multiplying by
The general solution is as follows:
where
After substituting Eqs. 15, 16 into eq. 14, a straightforward calculation proved that these equations are the solutions.
As mentioned earlier, the resolution is based on the introduction of an additive integrating factor
Eqs. 15, 16 can be analytically computed only for some simple power functions of r, such as:
Using Eq. 17 results in the following galaxy density and rotation velocity equations (see cyl_metric_1_r_z.wxmx):
Eq. 19 shows that the rotational velocity monotonically increases with r, which contradicts the observed profiles.
2.2 Rotational velocity on the equatorial plane
Eq. 18 shows that the transverse density profile is directly given by
We consider a galaxy for which the transverse
where ε is any strictly positive dimensionless small constant. For simplicity, and owing to the dimensionless nature of coefficients λ, β, and γ, we did not choose the normalized distribution form. Notably, even if the transverse density profile
With this distribution, the first- and second-order differentials of the metric components versus z vanish in the z = 0 plane, facilitating to an analytical solution to the Einstein field equations on the equatorial plane of the galaxy.
To handle the components
where N denotes a function of r. We could have kept the term
All the following equations in this study are formulated for the equatorial plane, with only the lowest order in each equation explicitly maintained (see ax_gal_z0. wxmx for the derivation of the following equations). On this plane, the cancellation of
We now begin with a domino cascade. Using Eq. 22, the resulting density equation at the lowest order is:
Eq. 23 shows that
The right-hand side of Eq. 24 shows that
The rotational velocity equation is:
As the classical centrifugal force is
The cross terms equation
The right-hand side member shows that in reality
Eqs. 26, 29 are identical and define the function N, the explicit knowledge of which is not required to predict the rotational velocity.
For typical galaxy density
Using Eqs. 20, 21 and directly integrating Eq. 27 to derive
Eq. 30 shows that the rotation velocity monotonically increases for any density distribution
Consequently, the universe-frame-dragging was missing. Indeed, a uniformly dense Universe is modeled using
2.3 Machian solution to the rotating Universe problem
The Einstein formulation of Mach’s principle presented in Section 1 implies that the complete treatment of a specific body cannot be performed without considering the impact of all the other bodies. Unfortunately, computing the metric of a rotating disk surrounded by uniformly dense matter at rest from the Einstein field equations is unfeasible. An easy way to address this issue in the study of finite-size bodies (stars, galaxies, and galaxy clusters) is to consider the following modified equation:
where
where
With Eq. 32, the null-Einstein tensor, Minkowski metric now corresponds to a uniformly dense Universe as it involves
Choosing a reference frame co-moving with the Universe, we have
Finally, Eqs. 28, 33 provide the rotational velocity of the galaxy:
The first square root of Eq. 34 represents the gravitational force, which is a monotonically increasing term. Interestingly, this term contains a virtual negative density
Figure 1 illustrates the prediction results of rotational velocity for the three galaxies M31, M81, and M101 (see M81_M101_M31_velocity.xlsx file) using Eq. 34 and the following density profile:
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1429235/fspas-11-1429235-HTML/image_m/fspas-11-1429235-g001.jpg)
Figure 1. Left vertical axis: circles are the rotational velocities observed by Roberts and Rots in 1973; solid lines are velocities predicted by Eq. 34 in the absence of dark matter halo; blue dotted line is predicted velocity of M31 without universe matter, i.e.,
The rotational velocities of the three galaxies were predicted using the same value of the current estimate of the mean Universe density present in Eq. 35, i.e.,
The measurement of the rotational velocity in 1973 of these three galaxies by Roberts and Rots (Roberts and Rots, 1973) initiated the dark matte halo genesis.
3 Discussion
The proposed solutions (Eqs 28, 33) accurately predicted the rotational velocity shapes of the galaxies. After reaching a peak, the velocity declined and finally vanished when the galaxy density
Although the model provides the gravitational field and motion equations (Eqs 27, 28), it does not provide any information regarding the actual metrics within the galaxy. The metrics in Eq. 1 or 21 correspond to the metrics for a lone galaxy in empty space. Therefore, no physical implications could be derived directly from these metrics.
As galaxies have low densities and stars move at low velocities, galactic rotation is widely considered as a non-relativistic problem, for which the Einstein field equations should not provide results very different from the Newtonian limit. However, special relativity tells us that making such a shortcut can be hazardous: the low-speed limit of kinetic energy in special relativity does not vanish as it does in the Newtonian mechanics, rather, it becomes
Mach and Einstein were concerned that, similar to constant motion, acceleration should make sense only relative to the other mass bodies. A later study showed that, according to the Einstein field equations, a body at rest but surrounded by masses in rotation undergoes inertial forces induced by the rotating masses (Pfister, 2007). Therefore, even if the stars are moving at low velocities, the centrifugal force experienced by a star can strongly depend on its relative rotation with the bulge of the galaxy and with the other celestial bodies too.
The proposed solution based on Eq. 34 clearly integrates this process. If the galaxy alone was modeled, i.e.,
When the Universe density
Eq. 36 differs from the Newtonian theory in three aspects:
- The term
This partial Machian property of the metric resulted from the fact that in the domino cascade Eqs. 23–26, the priority has been made to consider that χ should determine the galaxy density. This choice at the end results in a marginal impact of the term N in the final solution. In other axial metrics (Cooperstock and Tieu, 2007; Bronnikov et al., 2020), the cross term directly represented galactic rotation, which resulted in different solutions. Both solutions are mathematically valid but have different physical implications. Based on these considerations, Eq. 31 does not work with the other axially symmetric metrics.
- At the center of the galaxy, the inertial mass is primarily generated by the galaxy itself. When r increases, the inertial mass is gradually more impacted by the Universe and becomes equal to the classical value when:
The divergence of this term when
- The fact that it is not the whole galaxy mass included in the sphere radius r which appears in the gravitational force (
The two first properties showed that the present model satisfies the sixth version of Mach’s principle proposed in (Bondi and Samuel, 1997), i.e., “inertial mass is affected by the global distribution of matter,” with the scaling factor being
Definitely, Eqs 37, 38 shows that the inertial mass is affected by both, the galactic distribution of matter
Mach’s principle was intuitively proposed to explain the curve flatness with Newtonian dynamics (Darabi, 2013), where
An analytical solution of the first G order for a free-pressure rotating fluid has already been obtained to explain the flatness of the rotation velocity curve (Cooperstock and Tieu, 2007). An explanation for the rotation curve of the Milky Way as a consequence of the frame-dragging effect was recently proposed (Beordo et al., 2024). Both studies used a non-diagonal metric solution based on the Weyl–Lewis–Papapetrou coordinates.
The proposed metric expressed by Eq. 21 on the equatorial plane was obtained through an intuitive deduction process. First, a diagonal axial metric (Eq. 1) enabled the resolution of the Einstein field equations within a free-pressure stationary rotating fluid, i.e., without radial collapse. In this resolution, the cross-time-angular components of the stress tensor were neglected. Subsequently, the metric was slightly modified by adding a cross-time angular term. By keeping the fact that the galaxy density is determined by χ, the resolution of the field equations reveals that this cross-time term does not appear in the equatorial rotation velocity in the first density order (Eqs. 30, 34). We do not see any physical reason why this cross term could become significant when moving away from the equatorial plane. However, we could not prove this mathematically.
By transposing the density term of the Universe to the right-hand side of Eq. 32, it could be considered as an undetectable virtual negative density matter filling the entire space. However, it has some benefits compared to the dark matter halo thesis: it does not contain any adjustable parameter, while the observer is free to choose the dark matter halo distribution to explain their observations; it has a straightforward physical interpretation as the impact of the gravitational field of the Universe, which results in the implementation of Mach’s principle.
Eq. 34 includes the frame-dragging, which reduces the rotational velocity and cancels it when the galaxy density
Our success in predicting the rotational velocity profiles substantiates Mach’s ideas about the nature of inertial forces as an alternative to the dark matter halo thesis. Further investigations should be performed to determine whether the obtained solution does not violate any physical constraints. Finally, the impact of Eq. 32 in galaxy genesis and formation should be evaluated to support or discard the model. In this evaluation, statistical fluctuations could be added to
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://data.mendeley.com/datasets/x32hn4kxsh/1.
Author contributions
SW: Writing–original draft, Writing–review and editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Acknowledgments
The author thanks the reviewers for their valuable comments which helped in improving the paper.
Conflict of interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fspas.2024.1429235/full#supplementary-material
References
Barbour, J. B., Pfister, H., Beordo, W., Crosta, M., Lattanzi, M. G., Re Fiorentin, P., et al. (1995). Geometry-driven and dark-matter-sustained Milky way rotation curves with gaia DR3. Mon. Notices R. Astronomical Soc. 529, 4681–4698. stae855. doi:10.1093/mnras/stae855
Baum, R. P., and Sheehan, W. (2013). In search of planet Vulcan: the ghost in Newton’s clockwork universeSpringer.
Beordo, W., Crosta, M., Lattanzi, M. G., Re Fiorentin, P., and Spagna, A. (2014). Geometry-driven and dark-matter-sustained Milky Way rotation curves with Gaia DR3Monthly Notices of the Royal Astronomical Society529 (4), 4681–4698.
Beordo, W, Crosta, M, Lattanzi, MG, Re Fiorentin, P, and Spagna, A. Geometry-driven and dark-matter-sustained Milky Way rotation curves with Gaia DR3. Monthly Notices of the Royal Astronomical Society. 529(4), 4681–98.
Bondi, H., and Samuel, J. (1997). The Lense-Thirring effect and Mach's principle. Phys. Lett. A 228, 121–126. doi:10.1016/s0375-9601(97)00117-5
Brans, C. H. (1962). Mach's principle and the locally measured gravitational constant in general relativity. Phys. Rev. 125, 388–396. doi:10.1103/physrev.125.388
Cooperstock, F. I., and Tieu, S. (2007). Int. J. Mod. Phys. A 22, 2293–2325. doi:10.1142/s0217751x0703666x
Einstein, A., Feng, C.-C., Lin, L.-H., Wang, H.-H., and Taam, R. E. (1933). Glasgow: jackson wylie and Co. Astrophysical J. 785, 103. doi:10.1088/0004-637x/785/2/103
Feng, C. C., Lin, L. H., Wang, H. H., and Taam, R. E. (2014). Spiral density waves in M81. I. Stellar spiral density wavesThe Astrophysical Journal785 (2), 103.
Gödel, K. (1949). An example of a new type of cosmological solutions of einstein's field equations of gravitation. Rev. Mod. Phys. 21, 447–450. doi:10.1103/revmodphys.21.447
Hermann, M., Saravi, M., Hermann, M., and Saravi, M. (2016). Nonlinear ordinary differential equations: analytical approximation and numerical methods, 1.
Mach, E., Brockhaus, M. D., Graham, A. W., Moore, B., Diemand, J., and Terzić, B. (1901). Empirical models for dark matter halos. I. Nonparametric construction of density profiles and comparison with parametric models. Astronomical J. 132, 2685–2700. doi:10.1086/508988
Merritt, D., Navarro, J. F., Ludlow, A., and Jenkins, A. (2005). Astrophysical J. 624, L85–L88. doi:10.1086/430636
Milgrom, M. (2002). MOND—theoretical aspects. New Astron. Rev. 46, 741–753. doi:10.1016/s1387-6473(02)00243-9
Narlikar, J. (2003). Action at a distance and cosmology: a historical perspective. Annu. Rev. Astronomy Astrophysics 41, 169–189. doi:10.1146/annurev.astro.41.112202.151716
Navarro, J. F. (1996). The structure of cold dark matter halos. Symposium-international Astron. union. Camb. Univ. Press 171, 255–258. doi:10.1017/s0074180900232452
Keywords: gravitation, galaxy dynamics, dark matter, general relativity, Mach’s principle
Citation: Walrand S (2024) A machian model as potential alternative to dark matter halo thesis in galactic rotational velocity prediction. Front. Astron. Space Sci. 11:1429235. doi: 10.3389/fspas.2024.1429235
Received: 07 May 2024; Accepted: 03 July 2024;
Published: 01 August 2024.
Edited by:
Tonatiuh Matos, National Polytechnic Institute of Mexico (CINVESTAV), MexicoReviewed by:
José Domingo Arbañil Vela, Universidad Privada del Norte, PeruElmo Benedetto, University of Salerno, Italy
Copyright © 2024 Walrand. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Stephan Walrand, stephan.walrand@uclouvain.be