Modeling the Rising Tails of Galaxy Rotation Curves

: It is well known, but under-appreciated in astrophysical applications, that it is possible for gravity to take on a life of its own in the form of Weyl-curvature-only metrics (note that we are referring to the Weyl-only solutions of ordinary General Relativity; we are not considering Weyl conformal gravity or any other modiﬁed gravity theories), as numerous examples demonstrate the existence of gravitational ﬁelds not being sourced by any matter. In the weak ﬁeld limit, such autonomous gravitational contents of our universe manifest as solutions to the homogeneous Poisson’s equation. In this note, we tentatively explore the possibility that they may perhaps account for some phenomenologies commonly attributed to dark matter. Speciﬁcally, we show that a very simple solution of this kind exists that can be utilized to describe the rising tails seen in many galaxy rotation curves, which had been difﬁcult to reconcile within the cold dark matter or modiﬁed Newtonian dynamics frameworks. This solution may also help explain the universal ∼ 1 Gyr rotation periods of galaxies in the local universe.


Introduction
A confirmation of the nature of dark matter (DM) remains elusive. Weakly-interacting massive particles have recently been running up against strong constraints set by the null results of direct [1][2][3] and indirect [4] detection experiments. Astrophysical observations have also ruled out the original version of axions proposed to solve the strong CP problem [5,6], and searches for its "invisible" reincarnations have also become severely constrained [7,8]. Regardless of particle specifics, the cold DM (CDM) implied hierarchical structure formation paradigm faces a number of challenges such as the angular momentum problem [9], the missing clusters [10] and satellites [11] problems, as well as the issue of missing stochastic gravitational wave background [12]. On the other hand, the modified gravity theories have to overcome well-posedness issues [13,14] and problems such as unstable stars [15] (see also Famaey and McGaugh [16] for other complications with which specific theories have to contend).
Since new physics are not readily forthcoming, it is perhaps worthwhile to revisit the question of whether the proven General Relativity (and its weak field limit), together with the known Standard Model particles, can possibly already account for some of the apparent DM (aDM) phenomenologies. This prospect may be viable due to the presence of autonomous gravitational fields not sourced by any matter, whose existence has in fact been known since the early days [17,18] (see also, e.g., Brill and Hartle [19] or Misner and Taub [20] for a curvature-singularity-free warped universe that is completely devoid of matter; for complex manifolds, we also have the well-studied Calabi-Yau manifolds [21,22]). However, when trying to explain the observation of the mysterious gravitational content of the universe, which does not appear to interact through other forces, this most natural of possibilities significances, so additional topological constraints arise to fill the gap (spacetime cannot just bend arbitrarily if it is to close up into the correct topology). Imagine an idealized compact boundary-less universe (similar considerations can be applied to individual spatial slices), whose Chern-Pontryagin density [30] (R a bcd is the Riemann tensor, the Levi-Civita pseudotensor, and C a bcd the Weyl curvature tensor): and Gauss-Bonnet invariant (R is the Ricci scalar): integrate into the instanton number [31] and 32π 2 times the Euler characteristic [32], respectively, both of which are small integers for simple topologies. In other words, the average amplitudes of ρ CP and ρ GB should be on par with the inverse of the spacetime's volume. However, a cosmological constant Λ can create a surfeit of ρ GB (via R ab = Λg ab ) if the size of the universe is greater than 1/ √ Λ (∼10 Gly for the real thing). In order to achieve a recalibration of ρ GB without inadvertently bloating ρ CP , two new gravitational components should be introduced, at commensurate abundances to each other and to Λ, so near-cancellations can occur. The EWM and ordinary Standard Model matter are the obvious choices.
Incidentally, this desirability for two separate new ingredients is not apparent when examining the specific Friedman-Lemaitre-Robertson-Walker metric, whose oversimplifying assumption of exact isotropy artificially takes the Weyl half of GR and ρ CP out of action (there is no way to arrange the principal null directions of the Weyl tensor [33] without breaking isotropy, or equivalently, one can invoke the previous sectional curvature argument). In reality, the anisotropy on small scales is boosted by the two derivatives (multiplications by wavenumbers in momentum space) taking us from metric to curvature, so Weyl is far from negligible even with near-isotropy on large scales. In any case, regardless of whether this balancing act required by the compact example describes our actual universe, there is no reasoning that prevents the EWM from being present in it anyway.

Overview
To execute a preliminary assessment of the admissibility of the EWM as an aDM candidate, we examine the galaxy rotation curves that helped launch the field of DM research in the first place. As the assumptions of a weak stationary field and slow motion (as compared to the speed of light c) are reasonable, we can simply use the effective potential formalism (recall this is more general than Newton's law) and superimpose the potentials or accelerations from different origins.
Before going into any details, it is worthwhile re-emphasizing that we are not proposing a new modified theory of gravity. We are working completely within GR (and its weak-field slow-motion Newtonian limit). The governing equation that we solve is simply the Poisson equation ∇Φ = 4πρ M (whose solution is Φ = Φ H + Φ N , with Φ H being the solution to the homogeneous version of the equation and Φ N the particular solution to the inhomogeneous equation) satisfied by the Newtonian gravitational potential, which has dominated gravitational physics for centuries. What we are doing differently is that we also consider those solutions to this equation that have previously been thrown away due to oversimplified boundary conditions (i.e., we are not changing any theories, just trying to make sure that the relevant solutions are not inadvertently discarded).
The potential issue with boundary conditions we explore is that traditionally in astrophysics, it had been second nature to assume that the object in focus can be studied in isolation. For example, when trying to solve the gravitational field inside a star, we often assume that it is sitting all by itself in an otherwise empty universe, so that the gravitational field comes only from the star itself and drops to zero when we are far away from the stellar surface (such an asymptotically-vanishing boundary condition ends up setting Φ H = 0). This is not always valid of course, since the surroundings of the star may not really be negligible (e.g., the star may be in a binary with a black hole) and may exert strong influences on the interior of the star (e.g., tidally deform it). Inside the star then, such external influences manifest through a non-vanishing Φ H , since they are not sourced by matter inside of that star. Similarly with galaxies, their rotation curves had traditionally been examined by assuming that the galaxies exist in isolation (i.e., the gravitational field is implicitly assumed to vanish at large distances, thus setting Φ H = 0), even though we now know that they must be subsumed into larger-scale gravitational structures like clusters and filaments, so we must have Φ H = 0. In fact, the aDM seen in these larger structures could also be of Φ H in nature (i.e., not sourced by the matter content within the clusters, etc.), which are themselves determined by the boundary and initial (since at the largest scales, we cannot assume stationarity) conditions of the universe (please refer to the discussions in Section 1). When zooming in on an individual galaxy then, these clusters Φ H serve to produce non-trivial boundary conditions on the region around that galaxy, which we carve out as the computational domain to solve the Poisson equation.
Our new ingredient helping to explain the rising tails of the rotation curves is simply this non-vanishing Φ H , which has always been allowed by the Newtonian limit of GR (although Newton's law is narrower in scope than Poisson's equation and only gives Φ N ), but had unfortunately been thrown away by the implicit assumption that galaxies exist in isolation. When elevated back to GR language, this Φ H must be a part of the Weyl curvature tensor, because it is not sourced by in situ matter (the other half of the Riemann tensor, namely the Ricci tensor, must equate to the matter stress-energy tensor and thus cannot correspond to Φ H ). Therefore, in the bigger picture, the Φ H we consider below in the more specialized context of individual galaxies (explaining galactic-scale aDM phenomenology) is also an integral part of the autonomous gravitational content of the universe (explaining aDM phenomenology on all scales).
With autonomy, Φ H brings with it genuine additional freedom, which means that the present framework is more flexible than modified gravity theories such as modified Newtonian dynamics (MOND), which keeps Φ H = 0 and instead adjusts the dependence of Φ N on matter, and therefore does not really add any freedom unless this dependence is allowed to vary from situation to situation. As we have discussed in Section 1, this freedom allowed a decoupling between Φ H and Φ N , and thus consistency with the Bullet cluster observations. Another place where such flexibility may prove useful is when aDM exhibits variations across galaxies, e.g., if the claim of a DM-less galaxy by van Dokkum et al. [34] turns out to be valid, then the present framework can easily accommodate it if the local environment within the galaxy cluster is such that the galaxy in question experiences boundary conditions that lead to a vanishing Φ H (i.e., the galaxy is more isolated from its siblings than usual). We caution though that an alternative explanation is that since the EWM only becomes important at large galactocentric radii (see below) and observations of distant galaxies may simply be seeing only the brighter central core, so the lack of aDM may just be that the EWM-dominated regions are not seen.

A Galactic Effective Weyl Matter Solution
Our autonomous gravity component mimicking a galactic aDM halo is the very simple axisymmetric solution (this form provides a rising tail without introducing any singularities or discontinuities in the potential; in particular, the lack of a singularity at the origin fixes the powers): to the homogeneous Poisson's equation, which in the cylindrical coordinates {ρ, φ, z} reads: The gravitational acceleration associated with this solution is given by: and it is simple to double-check that the flux of this acceleration across the surface of any pillbox surrounding the galaxy (the galactic plane is assumed to be at z = 0) vanishes, so there is no mass sourcing it. This gravitational field is supported instead by the boundary conditions where the galactic aDM halo subsumes into the aDM distribution for the entire cluster (note that studies revealing a collusion between the halos of different galaxies in the lensing data [35], in their orientations [36], and on missing aDM in the local universe [37] suggest that these individual halos should not exist in isolation, but instead form an interconnected whole).
It is interesting to note that such boundary conditions (can be recovered by taking ρ and z to large values in Equation (3)) would entail larger acceleration values in the galactic corona, providing a possible mechanism to produce the highly energetic (10 5 -10 6 K) gases found in those remote regions. Concentrating on the galactic plane where we have more tracers and data, we note that the component g ρ (g z will be discussed in Section 2.4) along the plane increases linearly with ρ, translating into an "extended DM core" in the CDM terminology (i.e., a constant ρ aDM spherical density profile; there is thus no cuspy halo problem [38] in the galactic center), which has occasionally been noted in the literature (see, e.g., [39,40]) to be able to provide better fits to data than CDM-based profiles (e.g., [41,42]). To see the direct effects of g ρ , we need to go to very large ρ and study the tails of the galaxy rotation curves, which we turn to examine now.

Anatomy of Rotation Curves
Rising tails at large galactocentric radii ρ are seen for a large number of rotation curves (see the many references in Ruiz-Granados et al. [43]), yet such a trait cannot be explained with the theoretically-modeled CDM profiles, nor MOND [44], which predict curves that are flat or declining at very large ρ [39,45]. This is a cause for concern given that the tail is particularly informative, because while uncertainties in ordinary matter distribution (e.g., those associated with stellar mass-to-light ratios and the molecular gas densities) complicate interpretation of the inner and intermediate segments of the rotation curves, their impacts usually do not extend to the far outer regions, which consequently reflect the properties of the aDM in a more transparent manner.
Unfortunately, perhaps precisely because rising tails do not support either CDM or MOND, neither camp has paid much attention to it, especially in the written literature, sometimes chalking it up instead as due to low accuracy of measurements carried out at the extreme galactocentric distances (this is the reason why we included the dwarf galaxy in our Figure 1c below, which shows that the linearly-rising EWM profile is required not just for the tails, but also for the aDM contribution to not overshoot observed curves on the inside, where the accuracy argument does not apply; this issue is most salient for dwarf galaxies since the aDM contribution is appreciable even on the inside for them, so there is not a large baryonic contribution that masks the problem). In other cases, continuously-rising curves are interpreted as "have not yet reached the flat part by the last measured point" (see, e.g., the first paragraph of Section 7.1 in de Blok et al. [46], and also their Figure 57, which in fact displays many curves with rising tails). In other words, it appears that the extreme incompatibility of rising tails with CDM or MOND had prompted many authors to disregard them even as a possibility, and instead unconsciously interpreted them out of existence, even in cases where data for sufficiently large galactocentric radii are available and rising tails are unambiguously seen.
We have to caution nevertheless though that whether a rising tail is universal across all galaxies has not been firmly established observationally. Exploratory examinations of the extreme outer edges of the galaxies are not receiving sufficient observation time allocations, so it is possible that the absence of a rising tail from some curves may simply be due to data points not extending far enough out (by definition, the rising tails are regions where baryonic matter becomes subdominant, so there are few baryonic tracer objects to measure the rotation speeds with; further in from the tails, an ∼1/ √ r drop-off from baryonic contribution combines naturally with an ∼r rising aDM contribution into a rather flat curve in the transition zone where the two are of comparable strengths). We note in addition that our study can be used to infer an estimate on when the rising tails would appear (since κ appears common across galaxies; see below for details), and so, it is easy to verify or falsify with large-scale observation campaigns of galaxy outer regions.
More confusingly, with the existing surveys, people had been more interested in declines in the curves associated with mass distribution cut-offs, so even when data already exhibit early signs of rising tails, they are instead interpreted as "asymptotically approaching a flat curve after a temporary dip" (see many examples of this morphology in the agnostic sample of THINGS curves in Wojnar et al. [47]) and used to argue that earlier results showing declining curves do not provide the full picture, and making observations further out instead produces once again "flat" profiles. We do not quite agree with this "dip" interpretation (which assumes that the curves will flatten off rather than keep rising), because to produce this "rising back up" feature on top of the CDM or MOND profiles (which are smoothly declining or flat at large distances), baryonic matter must form a high density shell at extremely large galactocentric distances, and we have not seen any discussions as to why and how such a shell should form (especially since this morphology occurs quite frequently, the shell must be a ubiquitous robust feature of galaxies), or indeed seen any observational evidence for the existence of such shells at all. We believe that if observational efforts are carried further out to even larger distances, undeniable tails rising above any sensible definition of "flat" would become clear. The blue crosses are the HI -implied rotation speeds taken from Chemin et al. [45] and Oh et al. [48] and the blue curves their polynomial interpolations. The black curves in panels (b,c) are the measured gas contributions from Oh et al. [48]. The green curve in panel (c) is a CDM profile fit (see Equation (6)) demonstrating that such profiles tend to overshoot the observed rotation curve. Note that the rotation speeds add in quadrature, so accelerations add linearly.
In fact, the EWM model quite naturally produces such apparent "dip" features, as well as the declining curves that constitute another large population of the curves in, e.g., the SPARC sample. Because the EWM's |g ρ | profile rises linearly in a fashion more slowly than the CDM and MOND profiles, the inside regions of the rotation curves are more dominated by baryonic matter in this model. Therefore, when the baryonic matter is fairly concentrated (e.g., into the bulge), the rotation curve would naturally contain a declining mid-section where the baryonic contribution quickly declines (nearly Keplerian), but where the EWM contribution is still relatively small and so has not taken over yet. Only after going further out would the EWM contribution begin to show up. Therefore, depending on the galactocentric radii for which we have good data (as exemplified by truncating the dashed red line in Figure 1a at different places), one may see a (truncated) curve ending up in either a decline (in some extreme cases, the curves may appear to contain no dark matter contribution, which is difficult to explain with CDM or MOND, but quite natural in the present model; the EWM simply has not risen to an appreciable level yet at the observed distance), a "rising back up" convex shape, or a continuously-rising tail that rises above any sense of "flatness". In other words, we think that the EWM model is not only suitable for explaining the rising tails, it is also compatible with the other common morphological features of the rotation curves seen further in.
Coming back to the tail segment that is most reflective of the aDM properties, we note that a rising tail is quite natural with an EWM-generated linearly-rising |g ρ |. For concreteness, we carry out fits for a few example galaxies where high quality data for the rotation curve are available. Before starting though, we caution that there is significant difficulty with fitting to rotation curves, since besides the aDM contributions, there are also other unknowns like the mass-to-light ratio or molecular gas and compact object abundances, as well as more esoteric features such as whether a galaxy has recently gone through a merger so there are baryonic density waves. Even if measurements of "all" relevant quantities are available (which as far as we know is not the case for any galaxy, since at the very least, the molecular gases are difficult to see), they would contain large uncertainties, so it is not quite optimal to produce least squares fitting wellness parameters for entire curves, since the goodness of fits will have to be interpreted conditional on the aforementioned uncertainties associated with all the baryonic parameters. Therefore, the statistical significances of the fits are substantially degraded, rendering the fits much less informative. In short, due to the large uncertainties, there are too many handles in the baryonic models that one can crank to yield essentially any conclusion one wants about the aDM profile.
However, it is possible to make more robust quantitative measurements based on particular segments of the rotation curves. In particular, because baryonic matter density drops with increasing galactocentric radius ρ, their contribution to the overall rotation curve must rapidly decline, so for the tail section of the curves that we are interested in, uncertainties related to baryonic matter fortunately become subdominant. Therefore, by concentrating on the tail section, it becomes possible to make robust assessments, for which the impact of baryonic matter uncertainty is minimized (if for specific galaxies, high quality data on certain baryonic components are available, it also does not hurt to take them into account explicitly, so as to further reduce the uncertainty, a strategy that we will adopt with our examples below). This lack of contamination by baryonic matter uncertainties is fortunate, since our aim in this paper is precisely to explain, with the EWM, this rising tail end, which is the feature that is in strong conflict with CDM and MOND.
Specifically, to gauge quantitatively how well the EWM profile performs, we note that Equation (5) predicts a linearly-rising EWM contribution to the rotation curve that, in addition, must pass through the origin (i.e., v = 0 when ρ = 0). In other words, while even the simplest linear fit to the tail will produce two parameters, the slope and the intercept, the EWM profile is so rigid that it only has one parameter κ that can be adjusted. This "over-determinacy" translates into a rather stringent test for the EWM profile; if not correct, there is no way it will match both the slope and the intercept of the linear fit with a single parameter. We can also produce a quantitative measurement of this statement, namely by how much the slope of the EWM fit (varying only κ) differs from the slope of the best linear fit (varying both slope and intercept).
Although not the main focus of this paper, definitive assessments can also be made for the inner regions, even without prescribing specific models (thereby bringing in unwarranted assumptions) for the uncertain baryonic matter. After already fixing the definitive baryonic matter and EWM contributions, we place the remaining baryonic contribution that would produce a perfect fit (so the blue curve threading the observational points is the overall fitting curve if one prefers) into the "remaining baryonic" component in Figure 1. Definitive physical constrains exist for this left-over (it is not the fitting error and is supposed to be large since all the disc stars, etc., are in it), namely that it should have a shape resembling that arising from a disk occupying the small galactocentric radii regions. More stringently, this contribution has to be a real, as opposed to imaginary, number. In other words, the EWM plus definitive baryonic contribution should not have already overshot the observed curve, leaving negative headroom for the uncertain baryonic matter to contribute.
We begin with the nearby (thus good data quality) M31 (Andromeda). Our detailed fitting procedure is as follows 1.
Because M31 exhibits a strong bulge influence, it becomes possible to infer a bulge mass by fitting to the very inner part (ρ < 5 kpc) of the rotation curve. This is done because as discussed in the last paragraph, any reduction in the uncertainty about baryonic contributions will further reduce the errors when we fit to the tail section. The bulge contribution is modeled as a mass monopole Φ Nb = −M/ρ, and by inspection, we obtain M ≈ 4 × 10 10 M . See Figure 1a for the fitting results.

2.
We then turn to the tail section (beyond the dotted vertical line in Figure 1a at around 33 kpc). We begin by taking off the monopolar bulge contribution, so that the remainder contains the contribution from the EWM (dominate) and the other non-bulge baryonic components (subdominant).

3.
A fit for κ is then carried out. Equation (5) predicts an EWM contribution to the tail of the rotation curve at v = c √ κρ (where the speed of light c is in units of km/s), and we fit this functional form to the monopole-removed tail section of the rotation curve. The fitting is done using the standard computation package Mathematica, yielding the best fit κ = 5.035 × 10 −10 kpc −2 . For comparison with a full linear fit, it is more convenient to present the result in terms of the overall slope c √ κ, whose best fit value is 6.732, with standard error 0.0498. 4.
The highly rigid nature of the EWM profile (containing only one free parameter) as compared to CDM profiles (containing many tunable parameters) allows us to carry out a test on whether its rigid functional form matches data. Namely, as discussed previously, we can carry out a full linear fit to the tail with v = aρ + b and see whether a is close to the fit for c √ κ we obtained during the last step, as well as whether b is close to zero. Recall though that the tail still contains non-bulge baryonic contributions, and they would contribute to b; nevertheless, since they are subdominant in the tail section, b should not be large. Our fitting yields the best fit slope at a = 5.96, with standard error 1.29, which matches rather well with the c √ κ obtained earlier (well within one standard error). The parameter b has a best fit value of 27.65 and a standard error of 46.36, which is also consistent with a nearly-vanishing intercept.

5.
Finally, a few words regarding the mid-section of the rotation curve are in order. The remaining miscellaneous contributions, from the stars, dust, and gas residing outside of the central core, are collected into the "remaining baryonic" in Figure 1a. We caution that the study by Bournaud et al. [49] on recycled dwarf galaxies suggests the presence of large amounts of difficult-to-see cold molecular gas in the discs of their parent spiral galaxies (see also, e.g., Li et al. [50] for more direct observational evidences); thus, there is likely still an invisible matter component in the aDM (but of a mundane variety; being ordinary matter, their amounts and distributions could also differ substantially between galaxies of different types and ages, just as stellar matter would) that contributes to the inner to intermediate regions of the rotation curves.
Note further that all these aforementioned components are spread out on a disc, and not distributed in a spherically-symmetric manner, so the "remaining baryonic" curve is not bounded from below by a Keplerian profile. Instead, a rapid ∼1/∆ρ decline near a high matter density strip, ∆ρ being the distance to the sharp edge of the strip, provides a better approximation (and a slower decline in the orbiting speed is to be expected if the drop-off in density is more gradual), and is consistent with Figure 1a.
Beyond bulgy disc galaxies like M31, it is also interesting to examine dwarf galaxies. Because they are aDM-dominated, their entire curves should behave much like the tail, continuously rising throughout the available data range, a trend that is indeed seen in, e.g., Figure 57 of de Blok et al. [46]. Here, we examine in more detail two galaxies investigated by Oh et al. [48], which had been careful to scrub the contaminations due to non-circular motion (a significant problem with dwarfs) from their rotation curves. Beginning with NGC 2366 (Figure 1b), we note that it exhibits a pronounced tail section, which is once again well accounted for by the EWM profile. In contrast, even empirical CDM profiles proposed to resolve the cuspy halo problem cannot fit to it (see Figure 22 in Oh et al. [48]). Our detailed fitting procedure follows closely that of M31 and is as follows: 1.
The rotation curve does not contain a clear bulge-dominated inner segment, so we cannot infer a bulge monopole size, but this also means the bulge should not be significant (recall that we are dealing with a dwarf galaxy). Instead, the carefully-analyzed gas contribution derived from the integrated HI map is available for this galaxy, so we explicitly account for this contribution.

2.
Turning to the tail section, we once again remove the contribution from the known baryonic component, gas this time, so as to scrub down the tail maximally to a cleaner EWM domination. We then fit for κ with the same procedure as was done with M31, yielding best fit κ at 6.215 × 10 −10 kpc −2 . Note that despite the drastically different nature of the galaxies (massive spiral versus dwarf) and different rotation speeds (hundreds of km/s versus dozens), the κ values are remarkably similar between NGC 2366 and M31, differing by only 25% rather than orders of magnitude.

3.
The κ value translates into a slope c √ κ at 7.47916, with standard error 0.0518. The full linear fit, on the other hand, yields a slope of a = 6.907 with standard error 0.600, as well as an intercept b = 4.169 with standard error 4.356. Once again, the EWM's rigid form is consistent with the morphology of the tail section of the rotation curve (a and c √ κ agree within one standard error, and b is consistent with being nearly zero).

4.
The uncertain stellar contributions for this faint dwarf galaxy are collected into the "remaining baryonic" contribution.
For the fainter IC 2574 (Figure 1c), we need a somewhat different fitting procedure: 1.
Once again, the gas contribution from the integrated HI map is available. However, this time, the gas contribution rises very rapidly on the outer regions, so we end up with a large and varying gas contribution in the tail. This unfortunately obscures the linear EWM tail there (the curve bends downwards following the trend in the gas contribution). Nevertheless, the gas contribution declines at the extreme large ρ end, while EWM keeps rising to greater dominance there, so we can adopt a simpler procedure by letting the EWM match the last point in the rotation curve. This yields κ = 4.5 × 10 −10 kpc −2 , differing from the value for M31 by only around 10%, despite the very different galaxy types and rotation curve morphology.

2.
The very different (as compared to Figure 1a for massive spiral galaxies), continuously-rising rotation curve morphology (also more or less shared by NGC 2366) seen for the dwarfs is because for these aDM-dominated dwarf galaxies, the EWM contribution is significant throughout the entire curve, even on the inside. The almost linearly-rising curves thus provide a rather direct support for the rigid linear EWM profile. In contrast, they pose a serious challenge to common theoretical CDM profiles that flatten off at large ρ, as they cannot provide sufficient speed on the outside without overshooting the inside. Since it is difficult to produce outward acceleration, overshooting is a more troublesome problem. In Figure 1c, we provide a demonstration of this difficulty by making a Navarro-Frenk-White [41] profile fit as a green curve. The fitting is done by inspection of the tail of the curve beyond ρ ∼ 7 kpc, yielding parameter values v 200 = 160 km/s, C = 1, and R 200 = 100 kpc, feeding into: where x ≡ ρ/R 200 . Note that even though we have tried to suppress the overshooting by reducing C to a perhaps unrealistically small value (it is around 10-15 for the Milky Way) and by making R 200 very large (beyond this value, the CDM profile turns downwards), the inner regions of the CDM predicted curve still rest significantly above the observed rotation curve.
As already alluded to, the κ values for all three galaxies are quite similar, which could be understood in the EWM context as due to κ being determined by the fairly common general conditions in a local neighborhood of the cluster that set the boundary conditions for the individual galactic halos (both dwarf galaxies reside in the nearby M81 group, and M31 is in our local group; furthermore, the supra-galactic structures within the local group-e.g., planes of dwarf galaxies-have been shown to possess curious alignments [51], so the individual galaxies must all subsume into a common larger cluster-scale structure). We obviously cannot claim a quasi-universality of κ based on data from only three galaxies (the scrubbing technique of Oh et al. [48] should be applied to more galaxies to boost the statistics), but nevertheless note that it would be quite a coincidence if the commonality shared by these three rather arbitrarily-chosen cases is purely accidental, given that the galaxies belong to different classes, as well as exhibit different rotation curve morphologies and average speeds. Even more fortunately, establishing the universality of κ does not require us to obtain the full rotation curves. At very large galactocentric radii, EWM would dominate, leading to a linear v = √ κρ and thus a constant angular velocity √ κ (using our average κ ≈ 5.3 × 10 −10 kpc −2 from the three galaxies examined above, this translates into a rotation period for the outer rims of galaxies at about 0.89 Gyr). Therefore, with HI and other measurements of the outer regions of the galaxies, one should record a universal rotation period, shared between all galaxies regardless of their masses and the radii at which we happen to take readings observationally (such universality would not be present with truly asymptotically flat rotation curves), which matches and provides a robust (without needing additional assumptions such as sharp truncations of the discs) explanation for the observational results reported in Meurer et al. [52] for large sample sizes. Alternatively, recall that the EWM profile resembles that of a spherical CDM halo with a constant density proportional to √ κ (i.e., an extended core); the universality of κ across different galaxies then implies that the characteristic volume density (the overall scaling factor) in best CDM fits should be roughly constant, which is indeed observed to be the case for samples spanning over five decades of galaxy luminosity [53].
Furthermore, this universality in κ also provides a simple explanation for the well-known relationship between specific angular momentum j = J/M of a galaxy and its total stellar mass M, namely that j ∝ ∼ M 0.6 [54]. In the present consideration, we have that, for measurements where baryonic influence does not overwhelm the EWM contribution to the rotation curve (failure of this simplifying condition feeds into the rather significant scattering in the observed j vs. M relation), we can approximate the galaxy as a rigidly-rotating disk with angular velocity √ κ and some distribution of stellar mass density, say exponential ρ * = ρ core exp(ρ/H), where H is some scale distance. Then, simple integration gives: where the central stellar mass density ρ core is unsurprisingly dependent on M, as well, with ρ core ∝ M 0.42−0.74 observed for early-type galaxies [55]. Therefore, despite extreme crudeness, our estimate can already produce a M 0.6 power law, while it should be clear as well that, if the angular velocities of galaxies are completely arbitrary, we would not be able to obtain a clean dependence on M alone.
Note also that the functional form of Equation (7) is not sensitive to the detailed morphology of the galaxy due to dimensionality, but the precise numerical coefficient is, which implies that disc and bulge-dominated galaxies should naturally lie on two different parallel lines in a log j vs. log M plot, as is indeed observed. Another simple estimate shows that the universality of κ also leads to the baryonic Tully-Fisher relation [56]. Note first that the flat part of the rotation curve corresponds to where the decline and increase in v arising from the baryonic and EWM contributions balance out, leading to (at sufficiently large ρ, the bulge and disc contribution to v can be approximated by a 1/ √ ρ profile): where ρ f signifies the flat part of the rotation curve. Solving for ρ f and substituting into the rotation speed for that flat region, we have then: or in other words, M bar ∝ v 3 f when κ is common across galaxies. The observational value for the power index in the relationship between M bar and v f is currently under active debate, with the claimed results landing between three (e.g., Gurovich et al. [57] gives a range of 3.0 ± 0.2, and the recent study by Ponomareva et al. [58] yields 2.99 ± 0.2) and four (see, e.g., McGaugh [59]). The difference is largely attributable to how the mass-to-light ratio effects are handled, and the ongoing effort to extract a more accurate value is driven by the fact that these two values correspond to those anticipated under CDM and MOND, respectively. Our prediction is more in line with CDM, but if a definitive consensus value of four is eventually reached, it would serve to falsify our proposed model.

Off the Galactic Plane
A very distinguishing feature of the EWM as described (to lowest order approximation) by Equation (5) is the fact that the EWM is "repulsive" off the galactic plane. More specifically, g z is pointed away from the galactic plane, with an increasing strength at larger |z|, so it tends to (very slowly; the typical value for EWM acceleration is κ × 1 kpc ∼ 10 −12 ms −2 in SI units) remove matter that is not close to this plane or become dislocated from it (e.g., during galactic mergers); thus, it may possibly help to chisel out and maintain super thin pure-disc galaxies through the many mergers that they likely would have experienced (in order to grow to their present sizes), which the CDM framework struggles to account for [9,60].
Note though that the gravitational pull from regular disc matter points in the opposite direction to g z , so the stable disc plane is thickened. As some of that disc matter is dim and possibly lumped into aDM contribution, one may even on occasions infer that there is an overall dark attractive force towards the disc. The relative weakness of the "repulsive" g z is even more pronounced in the core region, where the high concentration of ordinary matter dominates (see Figure 2 below). EWM thus contributes little to bulge dynamics and certainly would not send the bulge matter flying away (large elliptical galaxies share similar dynamical characteristics to the bulges of disc galaxies, and their older age likely implies greater proportions of dim ordinary matter-as molecular gases, as well as compact objects-whose distribution and density are unfortunately uncertain). In short, absent major events such as mergers knocking (mostly gaseous) matter into high |z| regions, one should not expect to see significant outflows away from the galaxy when it is resting in a quiescent state.
Even during extreme events such as mergers, the dislocated material does not simply launch vertically into deep space. The combined acceleration directions from the EWM and the massive core of ordinary matter is plotted in Figure 2. We see that there are two regimes 4 : (1) matter starting off closer to the disc gets pushed towards the galactic center, but does so by first being carried to higher latitudes and then compressed nearly radially towards the galactic core; (2) matter starting off/knocked out sufficiently far from the disc would eventually get launched into deep space (which possibly explains how metals arrive at large distances of over 150 kpc from the galactic centers [61]), but in a direction more closely hugging the rotation axis of the galaxy. Spiral galaxies such as the Milky Way would frequently experience minor mergers with dwarf galaxies, during which stars (and gas) from the spiral galaxy would inevitably become dislocated off the disc and into region (1), and subsequently become lifted to moderately high galactic latitudes by the EWM field. The region (1) halo stars thus produced should, and have in fact been observed to, retain some features indicative of their disc ancestry, e.g., in their chemical composition [62], and prograde rotation [63]. Much of the stellar debris originating in the tidally-disrupted dwarf galaxy will however fall into region (2) and form a distinct outer population of halo stars with lower metallicity. The existence of this dichotomy of halo star populations is supported by observations (see e.g., Carollo et al. [63], Hartwick [64] and references therein), and the boundary between their spatial distributions is indeed located at around 10-15 kpc from the galactic center, matching that dividing EWM regions (1) and (2) shown in Figure 2 (i.e., the two EWM regions correspond to the observed inner and outer halos, respectively).
In the radial direction, the stars and gas knocked into Region (1) are guided towards the galactic center, helping to feed the growth of the supermassive black hole, leading to the consequence of an under-abundance of luminous matter in the disc as compared to the sizes of the supermassive black holes, which would then appear more massive than expected from a traditional CDM-based 4 Note that for simplicity, we had not included the matter contribution from the extended disc. Including it will not change the basic picture apart from pushing the separatrix streamline further up, enhancing the size of Regime (1). The plot also does not include the centrifugal forces arising from any circular motion, which would be particle specific. co-evolution scenario (see, e.g., Wu et al. [65]). Because merger events are episodic, the galactic core is fed material in discrete bursts. This could provide clues to the missing pieces in the formation mechanism of the mysterious Fermi bubbles [66,67]. These are giant structures seen above and below the galactic center (schematically marked out by the red dashed line in Figure 2) in γ-rays, microwave, and polarized radio signals [68][69][70], as well as in X-rays [71]. Digging through the archaeological evidence found within the bubbles, it has been suggested that they are outflows driven by a single past episode of active galactic nucleus (AGN) outburst and/or starburst activity occurring in the galactic core a few million years ago [66,72,73]. However, there is currently no contemplation as to where the extra material falling into Sgr A * and/or driving the star formation activities came from in the first place and why it is supplied in a short discrete burst. We note though that there is evidence for several recent minor mergers where the Milky Way had absorbed smaller galaxies, and when driven by the EWM, dislocated gases will take at least (centrifugal forces will likely prolong the journey) a quarter of a billion years to arrive at the galactic core, largely regardless of where they started off (solution to d 2 ρ/dt 2 = −κρ is ρ = ρ 0 cos( √ κt), so ρ = 0 when t = π/(2 √ κ), independent of initial ρ 0 ). Therefore, an initially-diffuse plume will automatically compactify into a small extent in ρ and feed either or both of the aforementioned AGN and starburst activities in a bursty fashion. Furthermore, this timetable is broadly in line with, e.g., the previous crossing of the Milky Way's galactic disc by the Sagittarius dwarf galaxy (estimated to have happened around 0.85 Gyr ago; see Purcell et al. [74]).
In addition, gravitational lensing may also provide clues to the structure of the aDM off of the galactic plane, if the lens galaxy is viewed edge-on. Figure 2 shows that the EWM should produce a quadrupolar contribution to the overall lensing potential (the acceleration field is the gradient to the potential; also recall that Figure 2 includes a monopolar contribution from the galactic core, which should be removed for the present consideration), essentially acting as an enhancement to the disc contribution (a monopole plus a quadrupole at the lowest orders), which may be powerful enough to resolve the radio flux-ratio anomalies. Specifically, the brightnesses of the multiple images of a lensed object often do not match traditional ellipsoidal CDM predictions based on the locations of the images (see, e.g., Goobar et al. [75] for an arbitrary example). This discrepancy is difficult to explain even when additional dark substructures are introduced into the modeling of the lens galaxies [76], as the CDM predicted substructure abundance is insufficient to account for the high incidence rate of such anomalies. It has been noted though by Hsueh et al. [77][78][79] that an edge-on disc can possibly help make up for the residual, provided the discs are massive enough (such masses are currently observationally unavailable). The present consideration offers another alternative, namely that the smooth galactic EWM, behaving much like a very massive disc in terms of producing higher multipoles in the lens potential, accounts for most of the anomalies (other factors such as free-free absorption make up the remainder), leaving little need for injecting additional low mass dark substructures.

Conclusions
Galaxy rotation curves exhibit features that require gravitational contributions beyond those arising from Newton's law of universal gravitation. Through concentrating on places most directly reflective of said new ingredients, we demonstrate in this paper that they can possibly be introduced within our verified understanding of gravity and particle physics, by noting that the gravitational field, like in any other field theories, does not always require source terms to exist; it can also be sustained by non-trivial boundary and/or initial conditions. On even larger length scales, the loss of symmetries poses a technical challenge to the analytical analysis of the EWM, but constructs such as the integral curves of the Weyl tensor's characteristic directions [80] may provide useful tools for extracting general attributes, e.g., along the lines of [81]. On a broader note, as any numerical relativist working on binary black hole mergers would attest, gravity is not always a surrendering slave to matter (black holes can be excised from the computational domain since no information comes out, and the violently twitching exterior vacuum is what is being simulated), and if some of the aDM turns out to indeed be the EWM, the reverse may in fact be more true of our universe. Substantial efforts likely need to be invested to ascertain whether this is indeed the case, with the first steps perhaps being the establishment of whether a rising tail is indeed a universal feature among all galaxies, by carefully examining their large galactocentric radius regions with the next generation of more powerful radio telescopes. The present paper aims to solicit interest from experts of diverse backgrounds to investigate the crevices where GR may be hiding further surprises for us, resulting in the more nuanced features of this more mundane gravity theory having to be considered in place or in addition to other DM considerations.
We must however avoid being overly optimistic and offer up a large dose of caution. The alternative, more established, CDM and MOND proposals have each gone through substantial amounts of scrutiny, and tremendous amounts of efforts have gone into resolving any difficulties they may face. In contrast, the present paper only handles a very specific phenomenology, with no attempt to explain things such as structure formation and other cosmological observations. These questions will have to be answered with extensive studies in the future, by experts in various fields who are more capable than the author, with results that may instantly disprove the proposed model. The present paper can only hope to generate some interest in carrying out those studies, and many more studies are required for the proposed model to be deemed viable at all.
Particularly challenging for the EWM model would be structure formation. This is because the autonomous gravity would behave rather differently from CDM (if one adopts a particle view of gravity, gravitons would be massless, but these massless particles are not like neutrinos or other hot dark matter, since they directly impress long-range forces-gravity-on ordinary matter). It is vitally important that the aDM and ordinary matter couple closely, so that the aDM halo and galaxies are at least spatially colocated. We believe that the autonomous gravitational field would still couple quite directly and closely with ordinary matter, because even though the autonomous gravitational field is not sourced by any matter, it would exert force and accelerate matter just as any gravitational field would. Therefore, just as CDM would attract ordinary matter into the potential well it creates, matter would also sink into any potential well in the autonomous gravitational field and arrange itself to align with the structures of the latter (so the EWM fields would still be centered on galaxies).
On the other hand, the feedback of ordinary matter to the autonomous gravitational field would be much weaker than for CDM and certainly MOND. In principal, GR is nonlinear, so ordinary matter can influence the EWM via the nonlinear coupling between the autonomous field and the gravitational field that the matter sources. However, at the galactic scales, both fields are weak in strength, and the coupling may need a very long accumulation timescale to produce anything significant at all. In contrast, the matter-sourced gravitational field would directly accelerate CDM particles (and pick up MOND corrections even more directly) and thus produce a much stronger feedback. We do not know if strong feedback is still required in the EWM context, since, e.g., we do not have the cuspy halo problem (the EWM profile is like a very extended core; there is not the large amounts of pressure-less CDM-like particles sinking into the center), so the viability of the autonomous gravity proposal in the structure formation context is essentially a very open question.