Galaxy rotation curves via conformal factors

We propose a new formula to explain circular velocity profiles of spiral galaxies obtained from the Starobinsky model in the Palatini formalism. It is based on the assumption that the gravity can be described by two conformally related metrics: one of them is responsible for the measurement of distances, while the other, the so-called dark metric, is responsible for a geodesic equation and therefore can be used for the description of the velocity profile. The formula is tested against a subset of galaxies taken from the HI Nearby Galaxy Survey (THINGS).


Introduction and motivation
Since there exist many issues that recently appeared in fundamental physics, astrophysics and cosmology which cannot be explained by General Relativity (GR) [1][2][3], one looks for other approaches which allow one to understand their mechanism. Classical GR is a well-posed theory. Many astronomical observations tested GR and have confirmed that it is the best matching theory that we have had so far for explaining gravitational phenomena. Unfortunately, GR is not enough to describe many unsolved problems such as late-time cosmic acceleration [4,5] (which one explains by existing an exotic fluid called Dark Energy introduced to the standard Einstein's field equations as cosmological constant), the Dark Matter puzzle [6][7][8][9][10][11][12][13]16], inflation [17,18], and the renormalization problem [19].
There are two main ideas competing for an explanation of the Dark Matter problem: geometric modification of the gravitational field equations (see e.g. [13][14][15]) or going beyond the Standard Model of elementary particles and introducing weakly interacting particles, which have failed to be detected [20]. In fact, these two ideas do not contradict each other and can be combined in some future successful theory. The existence of Dark Matter is mainly indicated by anomalies in a e-mail: ciprian.sporea@e-uvt.ro the observed galactic rotation curves. It interacts only gravitationally with visible matter and radiation, and also it has effects on the large-scale structure of the universe [21,22].
Many interesting and promising models have faced the Dark Matter problem. The most famous one is Modified Newtonian Dynamics (MOND) [23][24][25][26][27][28][29][30]. It has predicted many galactic phenomena and hence it is widely used by astrophysicists. Closely related is the so-called Tensor/Vector/Scalar (TeVeS) theory of gravity [31,32] which is, roughly speaking, the relativistic version of MOND. Another approach is to consider Extended Theories of Gravity (ETGs)-one modifies the geometric part of the field equations [33,34,36]. There are also attempts to obtain a MOND result from ETGs; see for example [37][38][39][40][41][42]. Another interesting proposal for explaining rotation curves is by using Weyl conformal gravity [43][44][45]. It should be noted that there is a model based on a quantum effective action and large-scale renormalization group effects [46,47], later on constrained by Solar System tests [48]. It provides, up to the first nontrivial order, a conformal transformation of the spacetime metric and a logarithmic term in the modified Newtonian potential.
In the following paper we would like to show how the socalled "Dark Metric" h μν [49,50], that is, the metric which is conformally related to the physical metric g μν appearing in an action of a theory under consideration, may explain the flatness problem of the galaxies' rotation curves. To this aim we employ the Ehlers-Pirani-Schild approach (EPS) [51] in the way it is presented and considered in Refs. [49,50]. The formalism assumes that the geometry of spacetime can be described by two structures, that is, conformal and projective ones. The first one is a class of Lorentzian metrics related to each other by the conformal transformation where Ω is a positive defined function (such a transformation is often interpreted as a change of frame, e.g. in scalar-tensor theories). The projective structure instead is a class of connections such that with A μ being a 1-form. Because of the positivity of the conformal function Ω, the conformal structure defines light cones as well as timelike, lightlike and spacelike directions in the spacetime considered. It should be noticed that it does not determine lengths of timelike and spacelike curves unless one chooses a representative of the conformal class. Geodesics in spacetime are defined by a connection. The different connections belonging to a considered projective structure define the same geodesics which are parameterized in two different ways [49]. The choice of a parametrization is related to a choice of clock, that is, a metric. We will say that the two discussed structures are EPS-compatible if the following holds (cf. (2)): A triple consisting of the spacetime manifold M and EPScompatible structures is called an EPS geometry [49].
We would like to emphasize that GR is a very special case of the formalism described. One assumes at the very beginning that the connectionΓ is a Levi-Civita connection of the metric g (the 1-form A is zero) which as a result means to treat the action of the theory as just metric-dependent. One may also treat the Einstein-Hilbert action as the one depending on two independent objects, that is, the metric g and the connectionΓ . This approach is called the Palatini formalism. Considering the simplest gravitational Lagrangian, linear in the scalar curvature R, the Palatini approach leads to the dynamical result thatΓ is a Levi-Civita connection of the metric g. It is not so in the case of more complicated Lagrangians appearing in ETGs. Moreover, as shown in [52], all Palatini connections of the form (2) are singled out by the variational principle.
Our aim is to show how important the EPS interpretation can be for an explanation of the galaxy rotation curves. It turns out, as expected, that at the end we deal with an expression which consists of the Newtonian part and some modification which depends on the theory one wants to study in the EPS approach. As the simplest example which we want to examine is the Starobinsky quadratic Lagrangian [17] in the Palatini formalism, which currently reaches very good results in the cosmological applications [53][54][55]. The starting point will be the standard geodesic equation from which we will derive the rotational velocity. It will be shown that the velocity can be written as GR plus extra terms coming from the conformal factor. Its usefulness is tested on a sample of six HSB galaxies. The conclusions and future ideas will be presented in the last part of the paper.

Velocities via conformal factors
Let us now derive a formula for the velocity of a star moving on a periodical trajectory in a given galaxy. For simplicity (and in good agreement with astronomical observations [56]) we will assume the orbit to be circular. In this case the centripetal acceleration and the velocity are related by On the other hand the Einstein equivalence principle remains valid for a theory of gravity that is conformally related with standard GR. This implies that a test particle (a star in our considerations) will satisfy the geodesic equation Stars can move around the galactic center at very high velocities. However, compared with the speed of light, the velocities are still very small; the condition v << c is always satisfied. Using the coordinate parametrization ( Under these conditions, together with the weak field limit of the geodesic equation (5) for a static spacetime (Γ 0 00 = 0), one obtains for the radial component Inserting now Eq. (6) into (4) we simply get v 2 (r ) = rc 2 Γ r 00 .
We have already discussed the idea of projective structures, that is, the class of connections related to each other by Eq.
(2). As already mentioned, connections belonging to the same projective structure describe the same geodesics but are differently parameterized. One needs to choose which metric from the conformal structure is connected to the geodesic motion. Let us consider the case when one deals with the Weyl geometry: the connection appearing above is a Levi-Civita connection of the conformal metricg ≡ h: One sees that Γ r 00 entering into Eq. (7) is which for a spherically symmetric metric takes the form and needs to be computed for a chosen model of gravity. If we consider any modified Einstein field equations (using the convention of κ = −8π Gc −4 from [57]) of the form where σ represents a coupling to gravity (for example a scalar field), one can write [58] Γ r 00 = We assume that the functionsB(r ) = B(r )/Ω andÃ(r ) = A(r )/Ω will take the Schwarzschild form in the weak field limit, that is, when we consider distances much smaller than the core size of a galaxy. In that case, the term g (1) 00 ≡ B (1) in the approximation h μν = Ω(η μν + g (1) μν ) represents the first order correction coming from the weak field limit of GR (see for example [57]) In the next section we will derive the exact form of (13) for a particular model of gravity which admits EPS interpretation.

An example: Starobinsky model
In principle there exists an entire class of gravity theories [49,51] that are conformally related with Einstein general relativity via Eq. (8). Our aim is to explain the observed galaxy rotation curves using Eqs. (13) and (7) without assuming the existence of Dark Matter.
In what follows we would like to propose a model that fits well the astronomical observed data on galaxy rotation curves. Our analysis is performed on a subset of galaxies obtained from THINGS: the HI Nearby Galaxy Survey catalogue [61,62], which is a high spectral and spatial resolution survey of HI emission lines from 34 nearby galaxies.
Any new model must take into account and reproduce the observed flatness of the galaxy rotation curves. At short distances (at least of the order of the size of the solar system) the velocity should have as a limit the Newtonian result v 2 (r ) = GM/r . This imposes some constraints on the functions A(r ) and B(r ).
Before we move on to the Starobinsky model, let us briefly recall the Palatini formalism. The action is where f (R) is a function of the Ricci scalarR = g μνR μν (Γ ), while S m is a matter action independent of the connection. The Ricci scalar is constructed by the metric-independent torsion-free connectionΓ . Varying the action with respect to the metric gives where the prime means differentiation with respect toR and as usually T μν is the standard (symmetric) energy-momentum tensor given by the variation of the matter action S m with respect to g μν . The g-trace of (17) arises as the structural equation of the spacetime controlling (17), Assuming that we are able to solve (18) asR(T ) we see that f (R) is a function of the energy-momentum tensor trace T , that is, T = g μν T μν . The variation of (16) with respect to the connection iŝ From the above it immediately follows that the connection ∇ is the Levi-Civita connection for the conformally related metric f (R)g μν . For a more detailed discussion we suggest to refer for example to [63][64][65]. The matter part of the modified Einstein equations will be considered as the perfect fluid energy-momentum tensor with the trace In the following work we will consider the pressureless case, that is, for dust p = 0. Moreover, we will assume only the radial dependence of the energy density, that is, ρ = ρ(r ). Let us notice [66] that Eq. (17) can be transformed into a so-called Einstein frame, where the Einstein tensor is constructed with the conformal metric h μν , whileT μν = T μν / f (R) and the effective poten- ). Thus, we see that (22) can be treated in a similar manner to Eq. (12) with W μν = − 1 2 h μν U . Now, we are equipped with all the tools needed in order to examine the model due to our assumption as regards the dynamics. As was already mentioned, the easiest modification of the Einstein-Hilbert action is adding the Starobinsky quadratic term [17]: where γ is a very small parameter having significance in the case of a strong gravitational field. From the cosmological consideration of the model in Ref. [55], with cosmological constant added additionally to the matter part, one finds that γ is of order 10 −11 . Since in the considered case the conformal factor is simply Ω = f (R) = 1 + 2γR and since the structural equation (18) for the dust matter gives rise tô we get particularly for the Starobinsky model Using the results from [58] one finds that where the modified mass distribution has the following form: Let us now assume a simple galaxy model obtained from the following matter distribution: (28) where r c can be interpreted as the "core radius", M 0 is the total mass of the galaxy and R 0 is the scale length of the galaxy to be matched with the observed one. Like in Ref. [59], the values of the parameter β are 1 for high surface brightness galaxies (HSBs) and β = 2 in the case of low surface brightness galaxies (LSBs). This matter distribution is a slightly modified version of the model used in Ref. [59] in order to be closer to the actual mass profile as inferred from the observed photometric profile. Although Eq. (28) is still a very rough approximation it is good enough for the purposes of this paper.
It can be noticed that the mass distributions (27) and (28) can be identified for a suitable choice of energy density ρ(r ) in (27); more exactly, comparing the derivativesM (r ) = M (r ) one gets an algebraic equation for ρ(r ).
Moreover, with the solution (26) we are able to see that Eq. (13) is now and therefore we write the quadratic velocity as Since the exact form of B(r ) is very complex, let us take the approximated value, that is, It follows immediately that the circular velocity of a star around the galactic center can be approximated by Table 1 Best-fit results according to Eq. (32) using the parametric mass distribution (28). These numerical values correspond to rotation curves presented in Fig  Equation (32) is the main result of this paper: the circular velocity obtained from the Starobinsky Lagrangian in the Palatini formalism. It is now used to obtain plots for a sample of six HSB galaxies after determining the parameters M 0 and r c using the NonlinearModelFit function in Wolfram Mathematica. They are presented in Fig.  2.
The corresponding values for M 0 and r c are presented in Table 1. We observe very good agreement between the data points and the fitted continuous black curve. Unfortunately, after plotting the Newtonian curves using the same values of M 0 and r c from Table 1 we have found that there is almost no difference (see Fig. 1) between the Newtonian curve and the one derived using Eq. (32). Furthermore, the mass-to-light M/L values inferred from the fit in Table 1 are also too large compared to what is expected based on stellar population synthesis models [60].

Conclusions
In this paper we have considered the possible explanation of the observed galaxy rotation curves by the assumption that the metric and the connection are independent objects in the spirit of EPS formalism. We studied the case when the connection is a Levi-Civita connection of a metric conformally related to the metric which is responsible for the measurement of distances and angles. Due to that interpretation, masses moving in a gravitational field should follow geodesics determined by the connection, providing different equations of motion. It turns out that the rotational velocity formula obtained under this formalism differs from the Newtonian one by the presence of extra terms coming indirectly from the conformal factor of the metrics. This term is treated as a deviation from the Newtonian limit of General Relativity.
In Sect. 2.1 we used Palatini gravity, which is a representation of the EPS formalism, and as a working example we took the Starobinsky Lagrangian f (R) =R + γR 2 in order to derive a rotational velocity formula given by the expression in Eq. (32) for a star moving in a circular trajectory around the galactic center. Our results are presented in Table 1 together with Figs. 1 and 2. Although the galaxy masses resulting from the fitting of the data sub-sample proved to be too high, giving thus rise to unsatisfactory values for the mass-to-light ratio, nonetheless we have showed that the approach of obtaining galaxy rotation curves via conformal factors can be valid. It should also be noticed that we have used a very simple mat- Fig. 2 (color online) Rotational velocities (in km/s) as a function of distance (in kpc). The black curve represents the best parametric fit of galaxy rotation curves using Eq. (32) for the sub-sample of six THINGS galaxies. The values of the best-fit parameters can be found in Table 1 ter distribution (28) in order to be able to obtain an expression for the energy density ρ(r ). More complex distributions could possibly give corrections which would provide different mass-to-light M/L values. This is one of the tasks of our future work.
Thus, we would like to briefly conclude by saying that the approach of obtaining galaxy rotation curves using two conformally related metrics can be valid and deserves further investigations. Furthermore, by trying other Lagrangians and other gravity models in the future, it is definitely possible to improve the findings reported here.