Exact cosmological solutions for MOG

We find some new exact cosmological solutions for the covariant scalar-tensor-vector gravity theory, the so-called MOdified Gravity (MOG). The exact solution of the vacuum field equations has been derived. Also, for non vacuum cases we have found some exact solutions with the aid of the Noether symmetry approach. More specifically, the symmetry vector and also the Noether conserved quantity associated to the point-like Lagrangian of the theory have been found. Also we find the exact form of the generic vector field potential of this theory by considering the behavior of the relevant point-like Lagrangian under the infinitesimal generator of the Noether symmetry. Finally, we discuss the cosmological implications of the solutions.


Introduction
We study a Scalar-Tensor-Vector theory which has been introduced as an alternative to particle dark matter [1]. This theory is referred as MOG (MOdified Gravity) in the current literature. Mathematically, MOG is complicated than GR in the sense that it postulates more gravitational filed than GR. In fact, MOG is a Scalar-Tensor-Vector theory of gravity, while GR is a tensor theory. In other words, in MOG in addition to the metric tensor, there are two scalar fields (µ(x α ) and G(x α )) and also a massive Proca vector field φ α . In fact the gravitational constant is variable as is in scalar-tensor theories of gravity. Therefore, the existence of these fields will provide some degrees of freedom which may help to handle the dark matter problem without invoking dark matter particles. It is worth mentioning that adding non-minimal scalar fields to the effective Lagrangian of gravitational filed dates back to Brans-Dicke theory [2]. In this theory the nonminimal scalar filed was introduced in order to incorporate MOdified Gravity theory a e-mail: mroshan@um.ac.ir the Mach's principle into GR. However, there are also other motivations for adding into play the scalar fields. For example, non-minimally coupled terms between scalar fields and geometric quantities, such as the Ricci scalar R, have to be added to the gravitational action when quantum corrections are taken into account. Also scalar-tensor theories show up as the dimensionally reduced effective theories of higher dimensional theories, such as Kaluza-Klein theory. For a comprehensive review of the subject, we refer the reader to [3].
On the other hand, the existence of the vector fields in the gravitational action has its own long story. In vectortensor theories, in the form of Einstein-ather theories [4], there is a Lorentz-violating vector field. The existence of this vector field can dramatically affect the cosmology. For example, it may leave an imprint on perturbations in the early universe [5], and it can even affect the growth rate of structure in the universe. Also the vector fields are a simple and natural candidates to explain some certain anomalies in the Cosmic Microwave Background [6]. For another example of theories containing vector fields, one may refer to Tensor-Vector-Scalar theory (TeVeS) [7]. TeVeS is a relativistic covariant theory for Modified Newtonian Dynamics (MOND) [8]. TeVeS has both types of above mentioned fields, i.e. scalar and vector fields.
As we mentioned before, MOG is also another example of theories where scalar and vector fields participate in the gravitational sector. This theory has been applied to explain the rotation curves of spiral galaxies and the mass discrepancy in the galaxy clusters [9]- [14]. Recently, it has been claimed that this theory is in excellent agreement with the rotation curve data for the Milky Way, while MOND does not fit the data [15]. In fact, the observed rotation curve data for the Milky Way extends as far as 200kpc from the center of the galaxy. Data for such a large distances from the core of the galaxy may provide a critical test of modified theories of gravity [15].
There are some papers in the relevant literature which have studied the consequences of MOG at the cosmological scale, for example see [16]- [20]. However, to the best of our knowledge, the astrophysical aspects of this theory has been investigated more extensively than its cosmological consequences. For example, there is no exact cosmological solution for this theory in the literature. This point motivated us to study the cosmological behavior of this theory by looking for exact solutions for the modified Friedmann equations. To do so, we use the Noether symmetry approach [21]. In fact, as every conservation theorem, the Noether conservation theorem allows to reduce dynamics of the cosmological model and gives insight into conserved quantities. This approach has been extensively used in various cosmological models and modified theories of gravity, for example see [22]- [33]. In [33] the Noether symmetry approach is used in non-minimal scalar-tensor theories and in higher order theories of gravity. Also a general review of the Noether symmetry approach can be found in this paper. In the current paper, i.e. in the context of MOG, there is also a non-minimal scalar field. We shall see that there are some cosmological exact solutions in MOG which correspond to some specific solutions presented in [33], for example see subsection 6.4. It is worth mentioning at this point that the Hojman's conservation theorem [34] has also been used to study the dynamics of cosmological models and to find exact cosmological solutions for them [35].
The layout of the paper is the following. In Section (2), we briefly review MOG's field equations and derive the corresponding modified Friedmann equations. In section (4), we derive the point-like Lagrangian which is needed for applying the Noether symmetry approach. Section (5) is devoted to the application of the Noether conservation theorem to MOG. In this section we find the explicit form of the generic vector field potential by symmetry issues. In Section (6), for the first time in the literature, we find some new exact solutions for MOG and discuss their cosmological implications. Conclusions are drawn in Section (7).

Modified gravity (MOG)
We start with the generic action of MOG where R is the Ricci scalar, S M is the matter action, Λ is the cosmological constant, ω 0 is a positive coupling constant, B µν = ∇ µ φ ν − ∇ ν φ µ , and the scalar fields χ and ψ are related to those of the original paper [1] as χ 2 = 2/G and ψ = ln µ. Also V φ is the self-interaction potential for the vector field φ µ and can be a function of φ µ φ µ and ψ. Variation of (1) with respect to g µν , φ α , χ, ψ leads to the following field equations respectively [36] G µν + Λ g µν = 1 Where G µν is the Einstein tensor and J α is a "fifth force" matter current defined as nonzero J α means that there is a coupling between matter and the Proca vector field φ µ . Consequently, one can verify that the matter energy-momentum tensor is not conserved and the Einstein's equivalence principle is violated [37]. Also, the total energy-momentum tensor is defined as where T µν is the energy-momentum tensor for the ordinary matter, and 32π ∇ µ ψ∇ ν ψ − 1 2 g µν ∇ α ψ∇ α ψ Furthermore, we assume a perfect fluid energy-momentum tensor for the ordinary matter. Also we consider the cosmological behavior of MOG in a flat Friedmann-Robertson-Walker (FRW) universe for which the line element is given by where a(t) is the cosmic scale factor. Now, by inserting this metric to the field equations (2)-(5), and also noting the relation between χ, ψ and G, µ, we find the modified version of the Friedmann equations. the 0-0 component of (2) giveṡ is the zeroth covariant component of the vector field and ρ is the total energy density of the matter/radiation. Also dot denotes the derivative with respect to time. The terms inside the bracket can be considered as corrections to the corresponding standard Friedmann equation.
Using the equation of motion of the test particles in MOG (the generalized version of the geodesic equation in metric theories), one can find an action for a pressure-less dust. Therefore, it can be shown that J 0 = κρ (see [37] for more detail). Where κ is another coupling constant. Although it seems that there are two independent coupling constants ω 0 and κ, they almost always appear as the combination κ 2 /ω 0 through the calculations of the physical quantities. Also it should be mentioned that because of the cosmological principle, spatial components of J α and all components of antisymmetric tensor B µν are zero. Furthermore the i-i components of (2) takes the following form where p is the pressure of the cosmic fluid, and we have used equation (9) for terms includingȧ 2 . Equation (3), i.e. the field equation of the vector field, takes the following form Therefore the potential V φ should satisfy the following general criterion for all values of µ(t) and φ 0 (t) Furthermore, equation (4) can be written as and finally equation (5) reads If the generic potential V φ and the equation of state p = p(ρ) are known, then equations (9),(10), (11), (13) and (14) are five non-linear partial differential equations for five unknown functions a(t), φ 0 (t), µ(t), G(t) and ρ(t). It is interesting to remark that the vector field is not "dynamical" in the sense that its field equation does not contain any time derivative. In other words, the vector field is algebraically related to the energy density ρ.

Vacuum exact solution
Before starting the Noether symmetry approach, we consider the vacuum exact solution (ρ = 0 and Λ = 0). In this case the modified Friedman equations can be analytically integrated for a wide range of generic potential V φ . Let us rewrite Equations (9),(10), (11), (13) and (14), respectively, as where H =ȧ/a and prime denotes derivative with respect to ln a. Also it should be mentioned that we have assumed that the Proca potential has a general form where f 1 and f 2 are arbitrary functions. f 2 satisfies f 2 (0) = 0 and d f 2 (0) dφ 0 = 0. Note that in the absence of matter/radiation the vector field freezes out at a finite value. We have assumed models for which the vector field freezes at φ 0 = 0. This is the case for original potential of MOG, i.e.
The exact solutions to equations (15)-(18) are where G 0 , µ 0 , ε, ε 1 and ε 2 are constants of integration. One may consider this solution as a late time solution where the energy density of matter and radiation is zero. If 2 < ε < 3 then we have an accelerated expansion, i.e.ä > 0. In this case G(t) is an increasing function of time and µ(t) can be decreasing or increasing. It is important mentioning that this late time acceleration is not a de Sitter universe where the cosmic scale factor grows as a(t) ∼ e √ Λt .

Canonical point-like Lagrangian for MOG
Equations (9),(10), (11), (13) and (14) can also be deduced from a canonical point-like Lagrangian L (a,ȧ, χ,χ, φ 0 ,φ 0 , ψ,ψ). For notational simplicity, in this section we use ψ and χ instead of µ and G. The point-like Lagrangian L is derived from action (1). We denote the configuration space by where q can be a, χ or ψ. Also the energy equation is therefore combining equations (22) and (23), we have five equations which make a complete set of partial differential equations for five unknown functions.
In order to find L one must insert the flat FRW metric into action (1) and simplify the result using an integration by parts. Note that the Ricci scalar curvature takes the form R = 6(ȧ 2 /a 2 +ä/a). As a result, one can easily verify that, in a FRW manifold, the Lagrangian related to the action (1) takes the point-like form where ρ = ρ m + ρ r is the proper energy density. Where ρ m is the matter energy density, and ρ r is the radiation energy density. It is worth mentioning that, in general, the mater energy-momentum tensor in MOG is not conserved (see equation A31 in [37]). However in the case of FRW universe we have ∇ µ T µν = 0 and consequently one can write ρ m = ρ m0 a −3 and ρ r = ρ r0 a −4 , where ρ m0 and ρ r0 are the current values of the corresponding energy densities. The current magnitude of the cosmic scale factor has been normalized to unity. It is interesting that although there is a coupling between matter and the vector field, the energy-momentum tensor of the ordinary matter is conserved. This is an enormous simplification made by the cosmological principle. Now by writing the Euler-Lagrange equations (22) for the Lagrangian (24), one can straightforwardly derive the equations (10)- (14). On the other hand, the energy condition (23) will give the equation (9).

Noether symmetry approach
Now let us introduce the lift vector field X [38]. X is an infinitesimal generator of the Noether symmetry in the tangent space T Q, and is defined as where the functions α, β , γ and ξ depend on the configuration space variables a, χ, ψ and φ 0 . One can straightforwardly show that where L X L is the Lie derivative of the Lagrangian along the lift vector X, and Σ is defined as Therefor one may conclude that Σ is a conserved quantity if L X L = −16πξ a 3 κρ (this is the Noether's theorem). It is needed to recall that in theories which there is no coupling between matter and the gravitational fields, the right hand side of this criterion is zero. By applying this condition to the Lagrangian, the analytic form of the functions α, β , γ and ξ will be specified. In fact, one may simply equate to zero the coefficients of terms likeψ 2 ,χ 2 ,ψχ and so on in order to find some differential equations for the lift vector's components. This is a common procedure in the relevant literature of the Noether symmetry approach, for example see [21]- [33]. Applying this procedure to the Lagrangian (24), we find the following differential equations Furthermore, there are some terms in L X L which do not contain any time derivatives of configuration space variables.
The summation of these terms should be separately zero: It is obvious from equations (28) that the function ξ can be an arbitrary function of a, χ and ψ. In other words, the existence of the Noether symmetry does not restrict the functional form of ξ . Therefore, one can construct infinite number of symmetry generators for this cosmological model. It is also worth mentioning that, there are nine differential equations in (28) which are enough for finding the exact form of α(a, χ, ψ), β (a, χ, ψ), γ(a, χ, ψ). After finding this functions and inserting them into (29), we shall find the form of the generic potentials V φ for which the Noether symmetry exists. After some straightforward algebraic calculations, one can verify that the general solutions to equation (28) are Where α 1 , α 2 and α 3 are integration constants. Although the Lagrangian (24) seems complicated, it possesses simple Noether symmetry generators. Furthermore, using the second equation in (30), the term containing the cosmological constant Λ vanishes in (29). Therefore, equation (29) can be rewritten as In order to find the exact form of the potential V φ it is needed to find the scale factor a(t) and the scalar field χ with respect to ψ and φ 0 . To do so, let us first derive the conserved quantity Σ associated to the Noether symmetry generator X. We recall that Noether's theorem states that every symmetry of the Lagrangian implies the existence of a conservation law, and therefore a conserved quantity. Using the Lagrangian (24) and equation (27), we find For the sake of simplicity we assume that Σ = 0. Thus one can easily verify that when α 1 = 0, we have where α 4 is a new integration constant. When α 1 is zero (i.e. when γ is constant), it is easy to show that Substituting these results into equation (31), we find where A is a constant and is equal to 1 2α 1 α 4 when α 1 = 0 and A = (3α 2 )/α 3 when α 1 is zero. It is important to note that, in principle, one may find the cosmic scale factor as a function of ∂V φ ∂ φ 0 from the field equation of the vector field (i.e. equation (11)). Therefore, equation (35) is a partial differential equation for potential V φ (φ 0 , ψ). This equation will determine the general form of V φ for which the Noether symmetry exists.

5.1
The case ρ r = 0 and ρ m , Λ = 0 In this case the energy budget consists of the cosmological constant Λ , non-relativistic matter, the scalar fields G and µ, and the vector field φ µ . Equation (35) takes a simple form and can be easily integrated. The result is where f (φ 0 ) is an arbitrary function of φ 0 , and we have used ψ(t) = ln(µ(t)). This potential covers a wide range of potentials. For example the original potential of MOG, i.e. V φ = − 1 2 µ 2 φ ν φ ν [1], lies in this subclass. Another, straightforward example is which in an isotropic and homogeneous space-time takes the form V φ = V 0 µ n φ m 0 . Also, it is obvious that if α 2 =0 (or equivalently A = 0) then V φ is a function of the vector field and does not depend on the scalar field µ. More importantly, every function of φ µ φ µ is a solution.

5.2
The case ρ m = 0 and ρ r , Λ = 0 In this case, we rewrite equation (35) as follows and the general form of the solution is in which y = exp(φ 0 )µ −A/3κ , and f is an arbitrary function of y. Simple examples of this class of potentials,

Cosmological exact solutions for MOG
Now we present some new exact cosmological solutions to MOG's modified Friedmann equations. Also we restrict ourselves to solutions with zero radiation density. We recall that by setting the Noether constant of motion Σ to zero, we found equations (33) and (34). Let us rewrite these equations as where c is a constant. This is a key equation to derive exact solutions. Using this equation we find some solutions for specific potentials introduced in the previous sections as potentials for which the Noether symmetry exists.

Model
Where λ 0 is a constant. We have used a new parameter q instead of A as A = −2q. This model corresponds to a scalartensor theory with two scalar fields µ(t) and G(t). In fact, since the matter and radiation are absent in the model, there is no source for producing the vector field. With the aid of equation (40) and assuming a power-law type solutions for the scalar fields as µ(t) = µ 0 t n and G(t) = G 0 t m , we find the following exact solution where G 0 is related to other constants as One can easily show that this solution corresponds to an accelerated expansion (i.e.ä > 0 andȧ > 0) if 3/20 < q < 1/ √ 5 or −1/ √ 5 < q < − 3/20. In this case, the magnitude of the scalar G(t) decreases with time and for µ(t) increases.
This solution has an interesting feature: the model can be contracting (i.e.ȧ < 0) whileä > 0. This situation happens when q > 1/ √ 5 or q < −1/ √ 5. In this case µ(t) is a decreasing function of time and G(t) is increasing. In fact the existence of the scalar fields µ(t) and G(t) can lead to repulsive gravitational force in the early universe.
6.2 Model V φ = − 1 2 µ 2 φ α φ α , ρ r = Λ = 0 and ρ m = 0 As we have already mentioned, this potential is the original potential postulated in MOG [1]. For this model we use equation (40) and also assume power-law type solution for G(t) and µ(t). One can easily verify the following exact solution a(t) = 3 544 9 π α 0 G N ρ m0 µ 0 2 where α 0 = κ 2 /ω 0 G N is one of the free parameters of MOG, and G N is the gravitational constant. Its observational values is α 0 = 8.89 ± 0.34 [12]. One can easily verify that the condition (12) is satisfied. It is obvious that the "velocity" of the vector field is zero and field's value remains constant in this matter dominated phase. Also this solution corresponds a decelerating universe. If we compare this solution with the matter dominated universe in GR (flat FRW space time with Λ = 0) where a(t) ∼ t 2/3 , then we see that the scale factor increases slower than the standard Friedman model. Thus the matter density also decreases with slower rate.
Where V 0 is a constant. With the aid of equation (40), one may easily verify the following solution where m can has two values: m = −18/7 and m = −3, and G 0 is related to the other constants as it is worth mentioning that the condition (12) holds if G 0 > 0. As we see the vector field is again constant. Also, the existence of two different values for m shows that our power-law type solutions are not the complete and unique solutions to the cosmological field equations. For m = −3 the cosmic scale factor increases linearly with time and so without deceleration, i.e.ä = 0. therefore scale factor increases faster than the standard matter dominated Friedmann model. Also we recall that in GR, for a matter dominated universe with flat spatial curvature and zero cosmological constant,ä can not vanishes. For m = −18/7 the scale factor increases as a(t) ∼ t 6/7 . therefore the scale factor increases faster than the standard case. However, unlike the case m = −3, we have a negative acceleration (ä < 0). As a final remark for this solution, we mention that the repulsive nature of the theory is again evident. In other words, the combination of the scalar fields and the vector field MOG can yield a repulsive gravitational force at the cosmological scale.

Model
Where m and n are real constant, and m = 0, 1, 2. Therefore this model is different from those presented in the previous subsections. Again, by using equation (40) which was obtained from the conserved quantity of the Noether symmetry, we find the following solution where K = ±2 √ 15 corresponding to n = ∓ √ 15/5 and the relations between constants are The criterion (12) leads to (2 − m)/G 0 > 0. If we assume that G 0 > 0 (note that G(t) is related to the gravitational constant), then everything is well defined provided that m < 2. It is interesting that this solution represent a constant acceleration for the cosmic expansion. We mention again that such a behavior is not achievable in the corresponding standard Friedmann model. As it is clear from equations (46), the scalar field G(t) is an increasing function of time, and the scalar field µ(t) depending on the value of parameter K can be an increasing function of time or a decreasing one. However, the rate of expansion is independent of the behavior of µ(t).
It is worth mentioning that G is a non-minimal scalar field. As we have already mentioned, non-minimal scalartensor theories have been already investigated using the Noether symmetry approach, for example see [33]. Although there is an extra vector field in MOG which is coupled to matter distribution, there are some similarities between the exact solutions of these theories. For example the exact solution presented in this subsection is reminiscent of the exact solution presented in equation (6.42) in [33]. The condition (12) is satisfied if µ 0 > 0 and G 0 > 0. For m = 5/4 the scale factor decreases as a(t) ∼ t −9/8 , and sȯ a < 0 andä > 0. However, for m = −20/11 the scale factor grows as a(t) ∼ t.

Conclusion
In this paper we have studied the cosmology of a specific scalar-vector-tensor theory of gravity known as MOG using a well-known approach. We have found the Noether symmetry generators and also the conserved quantity associated to the symmetry. More specifically, this approach helped us to find the explicit form of the vector field potential V φ (see equations (36) and (39)). We showed that the original potential postulated in this theory, is one of those potentials for which the Noether symmetry exists. Also, using the conserved quantity associated to the Noether symmetry, we have found some exact solutions for the modified Friedmann equations of MOG. We have found five set of exact solutions and except in one of them, we have not neglected the ordinary matter contribution. More specifically, our solutions correspond to matter dominated universes. As we mentioned before, power-law type exact solutions are not necessarily unique. However, these solutions are very often used in cosmological theories in order to insight into the physical content of them. In the context of MOG, there are some interesting features in the exact solutions which may help to understand the cosmological behavior of this theory. For example, the existence of the scalar and vector fields in MOG can lead to contracting universe with positive acceleration. In other words, these fields effect can behave like repulsive gravitational force at the cosmological scales. Also, there exist some solutions which the "gravitational constant" (the scalar field G(t)) increases with time as well as solutions where G(t) is decreasing.