Galaxy rotation curves and preferred reference frame effects

TL;DR

This paper explores an alternative gravity model involving scalar and vector fields to explain galaxy rotation curves without dark matter, showing good agreement with observed data and the baryonic Tully-Fisher relation.

Contribution

It introduces a generalized Jordan-Brans-Dicke scalar-vector-tensor gravity model with a preferred reference frame to account for galaxy rotation curves without dark matter.

Findings

Modified Schwarzschild-de Sitter metric in weak fields.

The model predicts galaxy rotation velocities consistent with observations.

Good correspondence with the baryonic Tully-Fisher relation for 12 galaxies.

Abstract

As an alternative to dark matter models we use generalized Jordan-Brans-Dicke scalar-vector-tensor (JBD-SVT) gravity model to study the behavior of the rotational velocities of test particles moving around galaxies. To do so we consider an interaction potential $U(\phi,N_{\mu})$ between the Brans-Dicke scalar field $\phi$ and time like dynamical four-vector field $N_\mu$ which plays as four velocity of a preferred reference frame. We show that at in weak field limits metric solution of the galaxy under consideration reaches to a modified Schwarzschild-de Sitter space in which mass of the vector field plays as an effective cosmological constant. In fact the present work proposes modification on the formulation of Newton's gravitational acceleration. This is used to explain circular velocity of galaxies without postulating dark matter. We also check our theoretical results with empirical baryonic Tully Fisher relation which states a linear relations between the rotational speed of galaxies and their mass. Mathematical calculation predict a good correspondence between our theoretical results and experimental observations for a set of 12 spiral galaxies.

Tables

Table 1. Table 1: Fitting parameters ( σ , r c , r q ​ f ) 𝜎 subscript 𝑟 𝑐 subscript 𝑟 𝑞 𝑓 (\sigma,r_{c},r_{qf}) for 12 spiral galaxies. Their observed data which are used from [ 5 , 45 ] are the relative stellar baryonic mass μ s ​ t ​ e ​ l ​ l ​ a ​ r = M s ​ t ​ e ​ l ​ l ​ a ​ r M s ​ u ​ n × 10 10 , subscript 𝜇 𝑠 𝑡 𝑒 𝑙 𝑙 𝑎 𝑟 subscript 𝑀 𝑠 𝑡 𝑒 𝑙 𝑙 𝑎 𝑟 subscript 𝑀 𝑠 𝑢 𝑛 superscript 10 10 \mu_{stellar}=\frac{M_{stellar}}{M_{sun}}\times 10^{10}, the relative mass of neutral Hydrogen μ H ​ I = M H ​ I M s ​ u ​ n × 10 10 subscript 𝜇 𝐻 𝐼 subscript 𝑀 𝐻 𝐼 subscript 𝑀 𝑠 𝑢 𝑛 superscript 10 10 \mu_{HI}=\frac{M_{HI}}{M_{sun}}\times 10^{10} in which galactic total mass is M H ​ I + M s ​ t ​ e ​ l ​ l ​ a ​ r subscript 𝑀 𝐻 𝐼 subscript 𝑀 𝑠 𝑡 𝑒 𝑙 𝑙 𝑎 𝑟 M_{HI}+M_{stellar}

$Galaxy$	$type$	${\mu }_{star}$	${\mu }_{HI}$	$\sigma$	$r_c$	$r_{qf}$	$\varsigma$
NGC 247	$LSB$	$0.4$	$0.13$	$-1.71138\times {10}^7$	$0.63315$	$0.193396$	$2$
NGC 1003	$LSB$	$0.30$	$0.82$	$48381.3$	$0.520752$	$0.431666$	$2$
NGC 2403	$LSB$	$1.1$	$0.47$	$1.70255$	$0.470889$	$0.416284$	$2$
UGC 6930	$LSB$	$0.42$	$0.31$	$-4652.3$	$0.393393$	$0.269871$	$2$
UGC 6983	$LSB$	$0.57$	$0.29$	$7.75637\times {10}^{10}$	$0.279081$	$0.321622$	$2$
NGC 300	$HSB$	$0.22$	$0.13$	$15468.$	$0.979994$	$0.155275$	$1$
NGC 3992	$HSB$	$15.3$	$0.92$	$9790.58$	$-0.299233$	$1.33564$	$1$
NGC 4157	$HSB$	$4.83$	$0.79$	$1.35773\times {10}^7$	$0.149683$	$0.793509$	$1$
NGC 5907	$HSB$	$9.7$	$1.1$	$3.89838\times {10}^7$	$-0.873617$	$1.26357$	$1$
NGC 6503	$HSB$	$0.83$	$0.24$	$20052.9$	$0.0169275$	$0.404241$	$1$
NGC 6946	$HSB$	$2.7$	$2.7$	$3.2588$	$0.196318$	$0.931982$	$1$
NGC 2903	$HSB$	$5.5$	$0.31$	$-11.4416$	$0.198349$	$0.726036$	$1$

Equations

$$I_T=I_{BD}+I_N$$

$$I_{BD}=\frac{1}{16\pi }\int d^{4}x\sqrt{-g}\left\{\varphi R-\frac{\omega }{\varphi }{g}^{\mu \nu }{\nabla }_{\mu }\varphi {\nabla }_{\nu }\varphi \right\}$$

$$I_N=\frac{1}{16\pi }\int d^{4}x\sqrt{-g}\{\zeta (x^{\nu })(g^{\mu \nu }{N}_{\mu }{N}_{\nu }+1)+2\varphi F_{\mu \nu }{F}^{\mu \nu }+U(\varphi ,N_{\mu })$$

$$-\varphi N_{\mu }{N}^{\nu }(2F^{\mu \lambda }{\Omega }_{\nu \lambda }+F^{\mu \lambda }{F}_{\nu \lambda }+{\Omega }^{\mu \lambda }{\Omega }_{\nu \lambda }-2R_{\nu }^{\mu }+\frac{2\omega }{{\varphi }^2}{\nabla }^{\mu }\varphi {\nabla }_{\nu }\varphi )\},$$

$$F_{\mu \nu }=2({\nabla }_{\mu }{N}_{\nu }-{\nabla }_{\nu }{N}_{\mu }),{\Omega }_{\mu \nu }=2({\nabla }_{\mu }{N}_{\nu }+{\nabla }_{\nu }{N}_{\mu })$$

$$\underset{\omega \to \infty }{\lim }\varphi \approx \frac{1}{G}$$

$$d{s}^2=-e^{2ϵ\alpha (r)}d{t}^2+e^{2ϵ\beta (r)}[d{r}^2+r^{2}d{\theta }^2+r^2{\sin }^2\theta d{\phi }^2]$$

$$e^{2ϵ\alpha (r)}\approx 1+2ϵ\alpha (r)+O({ϵ}^2)$$

$$e^{2ϵ\beta (r)}\approx 1+2ϵ\beta (r)+O({ϵ}^2),$$

$$\phi(r)\approx\frac{1}{G}\bigg{[}1+\epsilon\psi(r)+O(\epsilon^{2})\bigg{]}.$$

$$N_{\mu }(r)=(b(r),q(r),0,0)$$

$$-b^2(r)+q^2(r)\approx -1.$$

$$F_{tr}(r)=-2b^{\prime }(r),{\Omega }_{tr}(r)=2b^{\prime }(r)-4ϵb(r){\alpha }^{\prime }(r),$$

$${\Omega }_{tt}\approx -4ϵq{\alpha }^{\prime },{\Omega }_{rr}\approx 4q^{\prime }-4ϵq{\beta }^{\prime }$$

$$R_{tt}\approx -\frac{ϵ}{r}[r{\alpha }^{\prime \prime }+2{\alpha }^{\prime }]+O({ϵ}^2),R_{rr}\approx \frac{ϵ}{r}[2{\beta }^{\prime }+r{\alpha }^{\prime \prime }+2r{\beta }^{\prime \prime }]+O({ϵ}^2)$$

$$R\approx \frac{2ϵ}{r}[2{\alpha }^{\prime }+4{\beta }^{\prime }+r{\alpha }^{\prime \prime }+2r{\beta }^{\prime \prime }]+O({ϵ}^2).$$

$${\mathcal{L}}_{BD}=\frac{ϵ}{G}\left[r{\alpha }^{\prime }+2r{\beta }^{\prime }+\frac{r^2{\alpha }^{\prime \prime }}{2}+r^2{\beta }^{\prime \prime }\right]$$

$$+\frac{{ϵ}^2}{G}[(\alpha +3\beta +\psi )[r{\alpha }^{\prime }+2r{\beta }^{\prime }+\frac{r^2{\alpha }^{\prime \prime }}{2}+r^2{\beta }^{\prime \prime }]-\frac{\omega }{4}{r}^2{\psi }^{\prime 2}]$$

$$\mathcal{L}_{BD}^{eff}=-\frac{\epsilon^{2}r^{2}}{G}\bigg{[}\frac{\omega\psi^{\prime 2}}{4}+\frac{\alpha^{\prime 2}}{2}+3\beta^{\prime 2}+\frac{5}{2}\alpha^{\prime}\beta^{\prime}+\frac{\alpha^{\prime}\psi^{\prime}}{2}+\beta^{\prime}\psi^{\prime}\bigg{]}.$$

$${\mathcal{L}}_N={\mathcal{L}}_N^{(0)}+ϵ{\mathcal{L}}_N^{(1)}+{ϵ}^2{\mathcal{L}}_N^{(2)}+\cdots$$

$$\mathcal{L}_{N}^{(0)}=\frac{r^{2}}{G}\bigg{(}-10b^{\prime 2}+\frac{V_{0}(b)}{4}\bigg{)},$$

$$\mathcal{L}_{N}^{(1)}=\frac{r^{2}}{G}\bigg{\{}\frac{V_{1}}{4}+\frac{(\alpha+3\beta)V_{0}}{4}+[14\beta+2\alpha-10\psi+12b^{2}(\alpha-\beta)]b^{\prime 2}$$

$$-\frac{\alpha^{\prime\prime}}{2}+(b^{2}-1)\beta^{\prime\prime}+\frac{(b^{2}-1)\beta^{\prime}}{r}-\frac{b^{2}\alpha^{\prime}}{r}+8b(1-b^{2})b^{\prime}\alpha^{\prime}+4b(3-2b^{2})b^{\prime}\beta^{\prime}\bigg{\}}$$

$$\mathcal{L}_{N}^{(2)}=\frac{r^{2}}{G}\bigg{\{}\frac{V_{2}}{4}+\frac{(\alpha+3\beta)V_{1}}{4}+16b^{2}(b^{2}-1)\alpha^{\prime}\beta^{\prime}-8(b^{2}-1)^{2}\beta^{\prime 2}-8b^{4}\alpha^{\prime 2}$$

$$+\frac{(1-b^{2})\omega\psi^{\prime 2}}{2}+(b^{2}-1)(\alpha-\beta+\psi)\beta^{\prime\prime}+\bigg{[}\frac{(\beta-\alpha-\psi)}{2}+2b^{2}(\alpha-\beta)\bigg{]}\alpha^{\prime\prime}$$

$$+\frac{(\alpha -\beta )(b^2-1){\beta }^{\prime }}{r}+\frac{[b^2(\alpha -3\beta )-1]{\alpha }^{\prime }}{r}+[4(\alpha +3\beta )(3\alpha +4\beta )+2\psi (5\alpha +\beta )$$

$$+12b^2(\alpha -\beta )(\alpha +3\beta -\psi )]b^{\prime 2}+4[8\beta +3\psi +2b^2(5\beta -3\alpha -\psi )]b{b}^{\prime }{\alpha }^{\prime }$$

$$+4[\alpha+29\beta-2\psi+b^{2}(2\psi-22\beta+2\alpha)]bb^{\prime}\beta^{\prime}\bigg{\}}.$$

$$\mathcal{L}_{N}^{(1)eff}=\frac{r^{2}}{G}\bigg{\{}\frac{V_{1}}{4}+\frac{(\alpha+3\beta)V_{0}}{4}+[14\beta+2\alpha-10\psi+12b^{2}(\alpha-\beta)]b^{\prime 2}$$

$$-\frac{(b^{2}-1)\beta^{\prime}}{r}+\frac{(1-b^{2})\alpha^{\prime}}{r}+8b(1-b^{2})b^{\prime}\alpha^{\prime}+2b(5-4b^{2})b^{\prime}\beta^{\prime}\bigg{\}}$$

Full text

1 Introduction

Observations of the dynamics of galaxies as well as the dynamics of the whole Universe indicates that a main part of the Universe’s mass must be missing. According to the works done by authors in Ref. [1], this missing mass is possibly made by (the unknown) dark matter. So galactic scale dynamics is one of the important systems which is subjected to the dark matter studies. The observations of galaxies show that there is a discrepancy between the observed dynamics and the mass inferred from luminous matter [2, 3]. Since no dark matter has been detected so far, replacing dark matter by a modified gravity theory, is an alternative approach to the problem of missing mass [4]. There are different approaches to solve this problem, such as Modified Newtonian Dynamics (MOND) [5] and it’s relativistic extension [6], Modified Gravity (MOG) or its generally covariant version that called scalar-vector-tensor gravity theory(STVG) [4], conformal gravity [7], nonsymmetric gravity theory (NGT) [8] and nonlocal gravity [9, 10]. In STVG theory the dynamics of a test particle is given by a modified equation of motion. Since the metric field is coupled to scalar fields and a massive vector field, the solution of the field equations for a point mass is different from the point mass Schwarzschild solution of general relativity [11]. The predictions of STVG theory for the rotation curves of galaxies have been compared to observational data [12], by using a static spherically symmetric point mass metric solution for galaxies under consideration. The same approach has also been applied to the dynamics of globular clusters and clusters of galaxies [13, 14].   According to the works done by authors in refs. [15, 16, 17, 18, 19, 20] a new scalar-vector-tensor gravity theory is introduced where a time-like dynamical vector field can change a Lorentzian signature of the background to an Euclidian form. The dynamical vector field is usually considered as four velocity of preferred reference frames. Another approach was used to extend the Jordan-Brans-Dicke scalar-tensor gravity theory [21] and called as Jordan-Brans-Dicke scalar-vector-tensor (JBD-SVT) model by one of the authors in refs. [22, 23]. Several applications of this model are studied in refs.[24, 25, 26, 27]. The main motivations to present this JBD-SVT gravity model are as follows: The exact Lorentz invariance is impossible to test uniformly, as the boost parameter of this group is unrestrained. It also leads to the problem of divergences in quantum field theory associated with states of arbitrarily high energy and momentum [22]. This problem can be solved by a short distance high energy cutoff length which violates the Lorentz invariance [28]. Lorentz invariance violation causes to change the metric signature of space-time. For these reasons, we scrutinize the possibility that there is a preferred rest frame at each space-time point. If this frame was to be a fixed structure, it would violate general covariance, then the matter energy-momentum tensor is not divergence-less and the Einstein field equations is inconsistent. In order to preserve general covariance, the preferred frame should be considered dynamical. In this case dynamical frame could be defined by a vector field [29, 30, 31, 32] or by the gradient of a scalar field [33, 34]. One of the remarkable results of a suitable dynamical preferred frame, is to have two metrics with inequivalent causal structure named as Lorentzian and Euclidean metrics, respectively. The quantum cosmology viewpoint of the very early universe models describe its origin as a quantum tunneling from Euclidian to the Lorentzian space-time [22]. In both of them the space-time coordinates are still real coordinates. Whereas an ordinary way to obtain an Euclidean metric solution of the Einstein field equations, is to introduce a complex time variable $\tau=it$, so the Lorentzian signature of the metric is changed to its Euclidean version [22].

In accordance with the works has been done previously [4, 7, 8, 9, 10], in the weak field approximation the potential for a matter distribution of an extended object behaves like the Newtonian potential with an enhanced gravitational constant and an additional Yukawa potential. Applying a general potential $U(\phi,N_{\mu})$ and weak field approximation of the background metric obtained from dynamical field equations of the JBD-SVT gravity model [22] we study rotation curves of galaxies without postulating exotic dark matter which usually is described by dynamical scalar fields. This paper is organized as follows:

In section 2 we defined the JBD-SVT gravity [22] and calculate weak field limit of the Lagrangian density of the model for spherically symmetric state metric equation. Then we obtained linear order solutions of the Euler-Lagrange equations of the fields. In section 3 the acceleration equation and rotational velocity of a test particle is calculated in the weak field limit. Section 4 is devoted to the observational tests of the model. It is shown that the model is consistent with observational rotation curves of a sample of 12 spiral galaxies and the empirical Tully-Fisher relation. Last section denotes to concluding remarks and outlook of the work.

2 The model

Let us start with the following JBD-SVT gravity model [22]

[TABLE]

in which

[TABLE]

is JBD scalar-tensor gravity [21]. $g$ is absolute value of determinant of the metric tensor $g_{\mu\nu}$ where we use the metric signature convention as (-,+,+,+). $\phi$ is the Brans-Dicke scalar field and $\omega$ is the Brans-Dicke adjustable coupling constant. The second term of the action (2.1) is

[TABLE]

[TABLE]

for which

[TABLE]

describes action of a unit time-like dynamical four vector field $N_{\mu}.$ In other words it is the action of a preferred reference frame with four velocity $N_{\mu}(x^{\nu})$. Up to the terms of $\zeta(x^{\nu})$ and the scalar potential $U(\phi,N_{\mu})$, the action (2) is obtained by transforming the JBD action (2.2) under the metric transformation $g_{\mu\nu}\to g_{\mu\nu}+2N_{\mu}N_{\nu}.$ One can follow ref. [22] for more details. Action (2) shows that the vector field $N_{\mu}$ is coupled non-minimally to the JBD scalar field $\phi$ and the metric field $g_{\mu\nu}.$ Action (2.1) is written in units $c=G=\hbar=1$ . The undetermined Lagrange multiplier $\zeta(x^{\nu})$ controls $N_{\mu}$ to be an unit time-like vector field. $\phi$ describes inverse of variable Newton’s gravitational coupling parameter and its dimension is $(lenght)^{-2}$ in units $c=G=\hbar=1$. Present limits of dimensionless BD parameter $\omega$ based on time-delay experiments [35, 36, 37, 38] requires $\omega\geq 4\times 10^{4}.$ General relativistic approach of the BD gravity action (2.2) is obtained by setting $\omega\to\infty$ [40] for which one can infer

[TABLE]

where $G$ is the Newtonian coupling constant.

Let us consider a general spherically symmetric static metric equation which can be expressed in the following isotropic form [21]

[TABLE]

where $\epsilon$ is a constant. In order to study the behavior of JBD-SVT in astrophysical scales we should apply weak field approximation for the dynamics of the fields by perturbing the fields around Minkowski spacetime. In this case $\epsilon$ should have small values and will be order parameter of the perturbation. Furthermore, $\epsilon$ should be defined versus the parameter of the model which comes from non-linear counterpart of the action. In the BD action (2.2), $\omega$ determines nonlinear counterpart of the action. Hence it is suitable to choose $\epsilon=\epsilon(\omega).$ One can write perturbation series form of the fields which up to the second order terms are

[TABLE]

[TABLE]

and

[TABLE]

For static spherically symmetric metric (2.6) the vector field $N_{\mu}$ and the tensor fields $F_{\mu\nu}$ and $\Omega_{\mu\nu}$ should depend on the radial coordinate $r$ for which one can choose [39]

[TABLE]

where the time like condition $g^{\mu\nu}N_{\mu}N_{\nu}=-1$ in weak field limits reads

[TABLE]

Substituting (2.10) into (2.4) we obtain

[TABLE]

and

[TABLE]

where prime denotes to derivative with respect to $r$ coordinate and we used linear order terms of non-vanishing Christoffel symbols $\Gamma_{tt}^{r}(r)=\Gamma^{t}_{tr}(r)\approx\epsilon\alpha^{\prime}(r)$ and $\Gamma_{rr}^{r}\approx\epsilon\beta^{\prime}(r)$ for small $\epsilon$. Substituting (2.7) and (2.8) we obtain $tt$ and $rr$ components of the Ricci tensor for small $\epsilon$ as

[TABLE]

where the 4D Ricci scalar $R^{\mu}_{\mu}$ reads

[TABLE]

Substituting the above perturbative functions into the BD action (2.2) we obtain $I_{BD}=\int dtdr\mathcal{L}_{BD}$ in which $\mathcal{L}_{BD}$ is the Brans Dicke Lagrangian density such that

[TABLE]

[TABLE]

Integrating by part, we remove $\alpha^{\prime\prime}$ and $\beta^{\prime\prime}$ terms of the Lagrangian density (2) to obtain its effective counterpart which up to third order term $O(\epsilon^{3})$ is as follows

[TABLE]

If we use similar calculations for the action functional (2) then we will have $I_{N}=\int dtdr\mathcal{L}_{N}$ where $\mathcal{L}_{N}$ is

the Lagrangian density of the non-minimal interacting vector field which is defined by

[TABLE]

where

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

Integrating by part and removing $\alpha^{\prime\prime}$ and $\beta^{\prime\prime}$, the effective counterpart of equations (2) and (2) will have the form

[TABLE]

[TABLE]

and

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

Without reducing the generality of the issue, we assume that the interaction potential has the following form

[TABLE]

Adding (2.17), (2.19), (2) and (2) we can obtain total effective Lagrangian density of the system as follows

[TABLE]

where

[TABLE]

[TABLE]

[TABLE]

and

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

for

[TABLE]

The effective total lagrangian (2.25) shows that in the weak field limit and slowly varying Brans-Dicke scalar field where $\epsilon(\omega)\to 0$ and $\omega\to\infty$ the lagrangian density of the vector field is dominated instead of the Brans Dicke scalar field lagrangian density. In fact for $\omega\to\infty$ the Brans-Dicke action (2.2) reaches to the Einstein-Hilbert counterpart which can not support the galactic rotation curves alone. In the weak field limits $(\epsilon\to 0)$ one can see that the Brans Dicke Lagrangian counterpart (2.17) vanishes and so the vector field lagrangian density (2.26) is dominated to determine the galactic rotation curves at $\epsilon=0.$ Now we obtain Euler lagrange equations of the fields $b(r),~{}\alpha(r),~{}\beta(r)$ and $\psi(r)$ by varying the total Lagrangian density (2.25) and solve them order by order as follows. Setting $\epsilon=0$ the zero order term of the Euler-Lagrange equation reads

[TABLE]

which is independent of other fields and so can be solved alone. To solve the above equation we should choose some suitable potentials for which solutions of the above equation satisfy the time-like condition $q=\sqrt{b^{2}-1}$ (i.e. $|b|>1)$ with boundary condition $q(\infty)=0.$ Thus we choose the following ansatz for the potential $V_{0}(b).$

[TABLE]

in which $\mu$ can be described as mass of the vector field $N_{\mu}.$ We will see that $\mu$ is related with a suitable effective cosmological constant $\Lambda$ which means the Hubble constant of the expansion of the universe. Also the above potential makes the equation (2.30) as linear differential equation. We should point that, nonlinear differential equations in general form, have usually unstable solutions which may be reach to chaos. Substituting the above potential into the equation (2.30) one can obtain a particular convergent solution for $b(r)$ as follows.

[TABLE]

where $A$ is integral constant which should be considered as fitting parameter when we study its effect on the galactic circular velocity. The equation (2.30) with potential $(\ref{qeq})$ has other solution as $\frac{e^{\mu r}}{r}$ for $\mu>0$ which we do not consider here because it diverges to infinity at $r\to\infty.$ It is clear that only the decaying exponentials should be a part of the physical solution where the energy does not diverge to some infinite values and as we should retrieve the standard Newtonian gravitational potential at $r\to\infty$. Substituting (2.32), the equation $q(r)=\pm\sqrt{b^{2}-1}$ reads

[TABLE]

which satisfies the boundary condition $q(\infty)=0.$ At large distances $r\to\infty$ we see the vector field components approach to $(b,q)\approx(1,0)$ and $(b^{\prime},q^{\prime})\approx(0,0).$ Substituting these asymptotic solutions into the effective Lagrangian densities (2.26), (2) and (2) we obtain asymptotic behavior of the total effective Lagrangian density (2.25) as follows

[TABLE]

[TABLE]

Varying the Lagrangian density (2) with respect to the fields $\psi(r),\beta(r)$ and $\alpha(r)$ we obtain their Euler-Lagrange equations which up to second order terms $O(\epsilon^{2})$ become respectively as

[TABLE]

[TABLE]

and

[TABLE]

To solve the above equations we first substitute $2r\alpha^{\prime}$ from (2.35) into the equation (2.37) and integrate it to obtain

[TABLE]

in which $C_{1}$ is integral constant. The equation (2.36) can be integrated alone to obtain

[TABLE]

Now we substitute $\beta$ from the above equation and $r\alpha^{\prime}$ from the equation (2.35) into the equation (2.38) to obtain

[TABLE]

in which $C_{2,3}$ are integral constants. Substituting $\alpha^{\prime}$ from the equation (2.37) into the derivative of the above equation we obtain a linear differential equation for the Brans Dicke scalar field $\psi(r)$ as follows.

[TABLE]

One can solve explicitly the equations (2.41), (2.40) and (2.39) synchronously to obtain exact solutions for the fields $\psi(r),$ $\alpha(r)$ and $\beta(r)$ respectively as follows.

[TABLE]

where $C_{4,5,6}$ are integral constants and we defined

[TABLE]

and

[TABLE]

where we defined

[TABLE]

with

[TABLE]

[TABLE]

in which we defined

[TABLE]

[TABLE]

and

[TABLE]

Substituting the above solutions into the relations (2.7), (2.8) and (2.9) we obtain respectively

[TABLE]

[TABLE]

and

[TABLE]

Now we should fix the integral constants of the above solutions by regarding physical boundary conditions.

Appearance of the square term in the above metric potentials as $\mu^{2}r^{2}$ remember us the de Sitter universe at large scale structure of the metric solution where $\mu^{2}$ behaves as an effective cosmological constant $\Lambda>0$. This means that our solutions treat as de Sitter metric at large distances $\mu r>>1$ while at small distances $\mu r<<1$ the inverse distance factor $\frac{1}{r}$ is dominated instead of the term $\mu^{2}r^{2}$ which can be related to the Schwarzschild counterpart of the metric in absence of the terms $C_{4,5}$. We know that the Schwarzschild de Sitter black hole metric with $\Lambda>0$ is $ds^{2}=-(1-2GM/r-\Lambda r^{2}/3)dt^{2}+dr^{2}/(1-2GM/r-\Lambda r^{2}/3)+r^{2}d\theta^{2}+r^{2}\sin^{2}\theta d\varphi^{2}$ which is obtained from $G_{\mu\nu}+\Lambda g_{\mu\nu}=0$ and asymptotically reduces to the vacuum de Sitter space at large distances $r>>2GM.$ In weak field limit it reads $ds^{2}=-(1-2GM/r-\Lambda r^{2}/3)dt^{2}+(1+2GM/r+\Lambda r^{2}/3)dr^{2}+r^{2}d\theta^{2}+r^{2}\sin^{2}\theta d\varphi^{2}$ which can be compared with solutions (2.51) and (2.52) by setting

[TABLE]

where $M$ is total mass of the central black hole of the Schwarzschild de Sitter space time. In the standard $\Lambda CDM$ cosmological model the unknown exotic dark matter/energy is dominated to support the cosmic inflation and it behaves as an effective cosmological constant. Substituting (2.45), (2.48), (2.49) into the condition (2.54) the equation $K_{3}=K_{1}$ reads

[TABLE]

and the equation $K_{2}+K_{4}=0$ reduces to the following condition.

[TABLE]

in which

[TABLE]

Substituting (2.54) into the solutions (2.51), (2.52) and (2.53) we will have

[TABLE]

[TABLE]

and

[TABLE]

in which we defined

[TABLE]

Now we should say about physical situations of the integral constants $C_{4,5}.$ We can show that they are related to a radial length $r_{qf}$ scale which determines quasi flat regions of the space time where

[TABLE]

for which the equations (2.58) and (2.59) read

[TABLE]

and

[TABLE]

We know from the equation (2.5) which the solution (2.60) should not diverge to infinity at $\omega\to\infty$. In fact for quasi flat regions the solution (2.60) should reach to the boundary condition (2.5) for which we should set

[TABLE]

It is easy to check for large $\omega$ we will have asymptotically

[TABLE]

[TABLE]

and

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

and

[TABLE]

The above calculations show that for $\omega\to\infty$ we will have $J(\infty)=0=C_{4,5}(\infty)$ and so the equation (2.65) leads to small $\epsilon_{+}(\infty)\approx+0.5$ instead a zero value. This means the choice $\epsilon_{+}$ is physical but not $\epsilon_{-}.$ In fact, if we substitute $\epsilon_{-}(\infty)\approx-3$ into the equation (2.65) then we can result as $\lim_{\omega\to\infty}G\phi\to-2$ which is not a physical boundary condition ($\phi$ is a positive valued field). However for a large $\omega$ we substitute the above approximated equations together with $\epsilon_{+}\approx 0.5293+\frac{0.0020788}{\omega}$ into the equation (2.65) to obtain an equation for $r_{qf}$ as follows.

[TABLE]

Regarding $\vartheta$ as a free parameter one can obtain the following identities from (2.75).

[TABLE]

It is remarkable that the left side of the above equation shows that for large $\omega$ the quasi flat region $r_{qf}$ should reach to the Schwarzschild radius $2GM$ of the central black hole of the galaxy under consideration while the right side of them shows that by raising $\omega$ the quasi flat region $r_{qf}$ reach to the cosmological event horizon of the de Sitter vacuum space $\sqrt{\frac{3}{\Lambda}}\approx 1.833\times 10^{9}(light~{}year)=0.5623\times 10^{9}(Parsec).$ In other words for a fixed $\omega$ the left side equation shows that by rasing $r_{qf}$ then $\vartheta\to\frac{\pi}{2}.$ These helps us to obtain an approximated solution for $r_{qf}$ by setting $\vartheta=\frac{\pi}{2}$ in which we have

[TABLE]

At last we are in a position to write metric solutions for large $\omega$ as follows.

[TABLE]

and

[TABLE]

In the following section we obtain equation of motion of a test particle orbiting around galaxy in weak limit of the gravitational field.

3 Modified equations of motion

For galactic scales studies one can consider a test particle (a star in our consideration) as a tracer of gravitational field. At first approximation we often assume the galaxies are spherically symmetric objects. When we investigate rotation curves at significantly larger radii than the central region of the galaxy, it is justifiable to use a point-like source approximation for any physical model that does not involve an extended dark matter halo. These approximations are surprisingly useful because of the often very large uncertainty in the baryonic mass-to-light ratio(M/L) of many galaxies. We assume a test star with mass parameter $m$ orbiting around a central body and interact with the 4-vector field $N^{\mu}$ for which one can write its action functional as follows

[TABLE]

where $N_{\mu}=(b(r),q(r),0,0)$ given by (2.10) is 4-velocity of the preferred reference frame, $\tau$ is the proper time along the world line of the test particle with the dimension of length. $\lambda$ denotes to the interaction coupling constant between the vector field $N^{\mu}$ and the test particle with four velocity $V^{\mu}=\frac{dx^{\mu}}{d\tau}$. Varying (3.1) with respect to the coordinates $x^{\mu}$ one can obtain Euler Lagrange equations of the test particle such that

[TABLE]

For the line element (2.6), the equation of motion (3.2) in weak field approximation for small $\epsilon$ read

[TABLE]

[TABLE]

[TABLE]

[TABLE]

and

[TABLE]

Because of spherically symmetric property of the metric equation, we choose $\theta=\frac{\pi}{2}$ for planar orbit of the test particle which satisfies trivially equation (3.5). Therefore equation (3.6) leads to a conserved angular momentum of the test particle as follows

[TABLE]

in which $L$ is constant angular momentum of the test particle. In the weak field limit, slow motions of the test particle and for circular orbits we can substitute the approximations $\frac{dt}{d\tau}\approx 1$ and $\frac{dr}{d\tau}\approx 0$ into the geodesic equations (3.3) and (3). Equation (3.3) vanishes trivially while equation (3) leads to the acceleration equation of the test particle, which up to second order term $O(\epsilon^{2})$ is

[TABLE]

To obtain circular stable orbits of test particle (the test star) we should set $\frac{dr^{2}}{dt^{2}}=0$ and $L=mrv(r)$ in the above equation where $v(r)$ is circular velocity of the test particle moving on a circular orbit with constant radius $r.$ Regarding the latter conditions the equation (3.8) reads

[TABLE]

Substituting (2.32), (2.78) and (2.79) into the above formula we obtain galactic circular velocity as follows.

[TABLE]

where we defined

[TABLE]

[TABLE]

[TABLE]

and

[TABLE]

where we used (2.76) and

[TABLE]

It is easy to check that the above formula for the circular velocity reaches to the Newtonian regime $v_{N}=\sqrt{\frac{GM}{r}}$ at the quasi flat regions $r\approx r_{qf},$ if we set

[TABLE]

which reads

[TABLE]

for large $\omega.$ The equation (3.15) has two different solution as

[TABLE]

and

[TABLE]

$\vartheta_{1}$ takes $\frac{\pi}{2}$ for $\omega\to\infty$ and so it is not a suitable physical case because $\tan\vartheta$ given in the equation (3) diverges to a infinite value. Thus we will continue the work by using (3.17) to calculate the galactic circular velocity (3.10) which for large $\omega$ we choose $\sin^{2}\vartheta_{2}\sim\frac{1}{3}$ where (3) and (3.12) read respectively

[TABLE]

[TABLE]

[TABLE]

and

[TABLE]

It is easy to check

[TABLE]

and

[TABLE]

The above galactic circular velocity is for a point like source with total mass $M.$ Let us now extend the equation (3.10) for a mass distribution of spherically symmetric galaxies where the mass function

[TABLE]

is the total amount of ordinary visible matter within a sphere of radius $r$ and $\rho(\eta)$ is the density of visible matter of the spherical galaxy contains an inner core at radius $r=r_{c}$. There are different models which are proposed for mass distribution of visible galaxies. According to the model presented in ref. [12] we consider here a simple power-law mass distribution function as

[TABLE]

where $\varsigma$ is a constant parameter. The values for $\varsigma$ are $1$ and $2$ for HSB (high surface brightness) and LSB (low surface brightness) $\&$ Dwarf galaxies respectively. This difference is due to the fact that rotation curves of LSB and dwarf galaxies rise more slowly than those of HSB galaxies [43]. $r_{c}$ is radius of the galaxy core. Inside the core $r<r_{c}$ there is a constant mass density for HSB galaxies. For LSB galaxies it is a raising function as $\rho(r)\approx(r/r_{c})^{3}$ inside the core $r<r_{c}$ [12]. Well outside the core radius, where $r>>r_{c}$, equation (3.23) implies that

[TABLE]

in which $M$ is total mass of the galaxy under consideration. Numerical values for $r_{c}$ are given in the table 1 for a set of 12 observed galaxies. They are obtained by fitting the observational data and our solutions for circular velocity of the galaxies called in the table 1. At last we substitute (3.23) into the circular velocity given by (3.10) to obtain circular velocity for a mass distributed galaxy as follows.

[TABLE]

where $P(r)$ and $Q(r)$ should be inserted from (3) and (3.19). One can show that the Newton‘s gravity coupling constant $G\simeq 6.67\times 10^{-11}m^{3}/kg.s^{2}$ can be rewritten as G=4.3\times 10^{-6}\big{(}\frac{kpc}{M_{sun}}\big{)}\big{(}\frac{km}{s}\big{)}^{2}. In the latter case the galactic circular velocity can be described in units $km/s$ if we choose radial distance $r$ in units $kpc$ instead of the $metter$ and the galactic mass $M$ is given as relative mass with respect to the mass of sun as $(M/M_{sun})\times 10^{10}.$ Because v=\sqrt{\frac{GM}{r}}\to\big{(}\frac{km}{s}\big{)}\sqrt{4.3\times 10^{-6}\big{(}\frac{kpc}{r}\big{)}\big{(}\frac{M}{M_{sun}}\big{)}}. To plot numerical diagram of the circular velocity (3.25) we substitute radial distances in units $kpc$, and $G=4.3\times 10^{-6}$ and $M$ as $(M/M_{sun}\times 10^{10})$ (see table 1). In the latter case the velocity is described in $(km/s)$ units and $r$ in units $kpc.$ We plot numerical diagrams for the Brans Dicke parameter as $\omega=40000$ (see refs. [35, 36, 37, 38]). Other parameters given in the above formula called as $r_{qf},\sigma,r_{c}$ should be considered as fitting parameters when we set result of our theoretical model with observational data given in the table 1. We obtained numerical values for the fitting parameters $(r_{qf},\sigma,r_{c})$ and collect them into the table 1 for different galaxies which we used to study. circular velocity diagrams for 12 observed galaxies are given in the figure 1. There are three charts for each galaxy. That is, the experimental data (black dots) and the theoretical results (red dots) of our model and the Newtonian limit (blue dots) of velocity are plotted in terms of distance from the center of the galaxy. The galaxy mass is given from [5, 45]. There is two components for each galaxy mass which one of them is related to the neutral Hydrogen gases $\mu_{HI}$ surrounds the galaxy and the other is the central core star-like mass $\mu_{star}$ and so we should use the galactic relative total mass as $M=\mu_{HI}+\mu_{star}$ (see the table 1) for each galaxy rotation curve in the equation (3.25).

4 Spiral galaxy rotation curves

and observational data

Sample selection

We employ new database of SPARC (Spitzer Photometry and Accurate Rotation Curves) [44] to select observed data of circular velocity of a sample of 12 observed galaxies. They can be seen with $black~{}dotes$ in figure 1. In fact SPARC is sample of 175 nearby galaxies with new surface photometry at 3.6 $\mu$m and high-quality rotation curves from previous studies about the atomic hydrogen HI which is one of the best kinematical tracer of the gravitational potential in nearby galaxies [50]. In fact it is representing all rotationally supported morphological types of galaxies. To minimize the star-halo degeneracy, the best approach is to use near-infrared surface photometry (K-band or 3.6 $\mu$m), which provides the closest proxy to the stellar mass (see [50] and references therein). In short the SPARC is the largest galaxy sample to date with spatially resolved data on the distribution of both stars and gas. In order to test our model, we have therefore considered a sample of 12 galaxies with well-measured rotation curves extracted from Ref.[5] and presented in table 1.

Model fit to rotation curves of galaxies

To investigate the rotation curves of galaxies within the framework of JBD-SVT gravity, we suppose that there is no actual dark matter, therefore, such a galaxy essentially consists of baryonic matter containing stars and interstellar gas. Hence the observed circular speeds and the Newtonian ones are derived only from the observed mass of the stellar objects and HI (Neutral Hydrogen) components of galaxies. Moreover we ignore dust in our analysis, since the mass of the dust is at most a few percent of the mass of the interstellar matter.

Armed with equation (3.25) we are in a position to plots rotation curves for a sample of 12 spiral galaxies after determining the fitting parameters $\sigma,r_{c},r_{qf}$ (see table 1). To do so we use the ‘Nonlinear Model Fit‘ (NMF) function in Wolfram Mathematica software . We observe very good agreement between the observed data points and the fitted curve which are shown in figure 1. However, we must use a more sophisticated analysis for the general case where we do not have spherical symmetry, which is left for our next work. In fact the latter considerations bring some higher precision on the results. This is done for some spiral galaxies in ref. [46] but for MOG gravity model. It is worth to mention that for a broad range of even well researched galaxies there is no consensus on galaxy mass among various sources. Since the galaxy mass is not directly measured hence their estimations are based on various galactic models. So it is not surprising such a broad range of estimations on the mass and core radius of galaxies. However references of mass estimations which we have used in this work, are mentioned in table 1 and extracted from refs. [5, 45].

Baryonic Tully-Fisher relation

The baryonic Tully-Fisher relation (BTF) is an empirical relation between baryonic mass $M_{B}$ and maximum rotation velocity $v_{max}$ which may be roughly expressed as

[TABLE]

where $3.5\leq a\leq 4$ is obtained from data analysis by applying the different gravity models such as MOG, MOND and etc. [47]. Since we are interested in fitting rotation curves without any dark matter halo, the baryonic mass of a galaxy contains stellar and gaseous components only such as $M_{B}=M_{star}+M_{gas}$. For galaxies whose rotation curves can be resolved, $v_{max}$ can be chosen as different forms for instance: the maximum observed velocity, the average velocity, or the velocity at a fixed radius where the rotation curves seams flat. In ref. [48] one can see different choices for $v_{max}$ which are applicable to study the empirical Tully-Fisher equation. Since the rotation curves are close to the flat region so all of these choices are effectively equivalent. In order to construct a baryonic Tully-Fisher relation, one needs to estimate the maximum rotation velocity of a galaxy. There are several ways in which this can be done. We choose $v_{max}$ to be asymptotically value of the equation (3.25) at infinity $v(r>>r_{qf}).$ To do so we substitute (3.20) into the relation (3.25) for which we obtain

[TABLE]

in which we defined

[TABLE]

In the following we intend to investigate whether the results of JBD-SVT theory is in general agreement with the BTF?. A log-log diagram of the maximum speed $v_{max}$ given by (4.3) is plotted versus the baryonic mass of 12 observed galaxies in figure 2. The least-square fit to the data points is given by the following straight-line equation

[TABLE]

in which

[TABLE]

[TABLE]

and

[TABLE]

One can see that the equations (4.5) and (4.6) are in a good agreement with BTF given by (4.1) alone and slope of them are $4$ approximately, while a combined LSB $\&$ HSB given by (4.7) dose not satisfy the BTF and its slope reaches to some small values with respect to $4.$ Difference of slopes with respect to the prediction in ref. [47] as $3.5<a<4$, maybe resolved by regarding possibility observational errors, or by using more samples of the observational data, regarding cylindrical symmetry of the used spiral galaxies or by choosing other samples for the maximum circular velocity defined in ref. [48]. Authors in most articles are used more than 100 sample of observed galaxies data to obtain slope of the empirical Tully Fisher relation. Since BTF is an empirical relation so the existence of the Baryonic Tully-Fisher relation may implies that the mass observed in baryons is the total mass [49] and it challenges the dark matter hypothesis.

5 Concluding remarks

The rotation curves of galaxies still remain as one of the profound challenges in the present day physics. In this paper we have considered an alternative view to the dark matter problem. In fact we used generalized Jordan-Brans-Dicke scalar-vector-tensor (JBD-SVT) gravity model to generate galaxy rotation curves where the time-like vector field $N_{\mu}$ is coupled non-minimally to the BD scalar field $\phi$ and the background metric. In weak field limits, we considered a power-law self-interaction potential for the vector field, which plays an important role to produce a repulsive Yukawa like metric potential. In fact we obtained that the vector field mass parameter is related to an effective cosmological constant which supports the acceleration of the universe in the large scalae structure. In other word we obtained that in weak field limits metric around a galaxy under consideration behaves as modified Schwarzschild de Sitter space time. However this corrects galactic rotation curves by regarding the experimental observational data. The results were surprisingly the same as ones which are given in [12] but provided their physical basis are totally different. To test the observational consequences of the circular velocities in our model, we used the well-measured SPARC database [44] to fit the theoretical rotation curves predicted by JBD-SVT to the observational data. Our results are appropriately consistent with the observations. We then demonstrate that our results are in a good agreement with the baryonic Tully-Fisher relation for spiral galaxies. In short, the outlook of this work can be as follows: A vacuum sector of the Brans Dicke scalar tensor gravity in weak field limit which reduces to a Newtonian approach of the general relativity can not give a good fit to the galactic rotations curves (see blue dotes in figure 1). While a time-like dynamical self interacting vector field moving on the curved background metric of a galaxy can produce a good correspondence between theoretical results (see red-dotes in figure 1) of galactic circular velocities and observational data (see black dotes in figure 1). On the other hand, baryonic Tully Fisher relation confirms that one can describe galactic rotation curves correctly without using unknown cold dark matter and just with visible baryonic matter, using JBD-SVT theory. As further investigations one can study the solar system tests of the GBD-SVT theory, cosmological implications of the theory such as domain walls and study the form of modified virial theorem in the theory.

**Acknowledgments

**The authors are grateful to the editor and anonymous referees for their valuable comments and suggestions which cause to improve the article for readers.