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-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–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-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 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-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-42]. Another interesting proposal for explaining rotation curves is by using Weyl conformal gravity [43-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 so-called “Dark Metric” [49, 50], that is, the metric which is conformally related to the physical metric 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 transformationwhere 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 thatwith 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 EPS-compatible 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-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.

The metric signature convention is .

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 byOn 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 equationStars 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 is always satisfied. Using the coordinate parametrization it follows immediately that if then also . Under these conditions, together with the weak field limit of the geodesic equation (5) for a static spacetime (), one obtains for the radial componentInserting now Eq. (6) into (4) we simply getWe 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 metric :One sees that entering into Eq. (7) iswhich for a spherically symmetric metrictakes the formand needs to be computed for a chosen model of gravity. If we consider any modified Einstein field equations (using the convention of from [57]) of the formwhere represents a coupling to gravity (for example a scalar field), one can write [58]We assume that the functions and 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 in the approximationrepresents 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 . 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 iswhere is a function of the Ricci scalar , while 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 giveswhere the prime means differentiation with respect to and as usually is the standard (symmetric) energy- momentum tensor given by the variation of the matter action with respect to . The g-trace of (17) arises as the structural equation of the spacetime controlling (17),Assuming that we are able to solve (18) as we see that is a function of the energy-momentum tensor trace T, that is, . The variation of (16) with respect to the connection isFrom the above it immediately follows that the connection is the Levi-Civita connection for the conformally related metric For a more detailed discussion we suggest to refer for example to [63-65].

The matter part of the modified Einstein equations will be considered as the perfect fluid energy-momentum tensor with the traceIn the following work we will consider the pressureless case, that is, for dust . Moreover, we will assume only the radial dependence of the energy density, that is, .

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 , while and the effective potential . Thus, we see thatcan be treated in a similar manner to Eq. (12) with .

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 . Since in the considered case the conformal factor is simply and since the structural equation (18) for the dust matter gives rise towe get particularly for the Starobinsky modelUsing the results from [58] one finds thatwhere the modified mass distribution has the following form:Let us now assume a simple galaxy model obtained from the following matter distribution:where can be interpreted as the “core radius”, is the total mass of the galaxy and 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 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 in (27); more exactly, comparing the derivatives one gets an algebraic equation for .

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

Moreover, with the solution (26) we are able to see that Eq. (13) is nowand therefore we write the quadratic velocity asSince 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 byEquation (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 and using the NonlinearModelFit function in Wolfram Mathematica. They are presented in Fig. 2.

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. 2. Col. (1) galaxy name; Col. (2) total gas mass, in units of , given by , with the data taken from [61]; Col. (3) measured scale length of the galaxy in kpc; Col. (4) galaxy luminosity in the B-band, in units of , calculated from [61]; Col. (5) presents the best-fit results for the predicted total mass of the galaxy (in units); col. (6) gives the predicted core radius in kpc; Col. (7) reduced ; and Col. (8) the stellar mass-to-light ratio in units of . Note: all six galaxies are of type HSB

Galaxy	M_gas	R_0	L_B	M_0	r_c	χr2	M / L
(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)
NGC 3031	0.48	2.6	3.049	14.86	2.10	4.88	4.71
NGC 3521	1.07	3.3	3.698	38.45	3.69	1.84	10.10
NGC 3627	0.11	3.1	3.076	8.68	2.25	0.45	2.78
NGC 4736	0.05	2.1	1.294	0.53	0.59	2.41	0.37
NGC 6946	0.55	2.9	2.729	78.19	5.09	2.18	28.44
NGC 7793	0.12	1.7	0.511	18.24	3.36	4.82	35.45

The corresponding values for and 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 and 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].

Fig. 1

The difference (vertical axis in Km/s) as a function of distance (horizontal axis in kpc), where v is given by Eq. (32) and is the Newtonian velocity . The galaxies from top to bottom are as follows: NGC3031, NGC3627, NGC3521, NGC6946 and NGC7793. The galaxy NGC4736 is not shown, but for it always

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 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 matter distribution (28) in order to be able to obtain an expression for the energy density . 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.