A study of positive energy condition in Bianchi V spacetimes via Noether symmetries

In this paper we use Noether symmetries of the geodesic Lagrangian in Bianchi V spacetimes to study various cosmological solutions of Einstein's field equations. Our first result is the identification of the subalgebras of Noether symmetries of the equations of motions in such spacetimes with dimension 4, 5, 6, 7, 9 or 10 of the maximal algebra of Lie point symmetries of dimension 13. Secondly we give physical interpretation of new cosmological solutions which satisfy positive energy condition and yield critical bounds on the expansion coefficient $\alpha$, in which the underlying non-flat spacetimes carry interesting physical properties. Specifically the energy density behaves in one of the following ways. (i) It is positive and constant for all time. (ii) It varies with time and attains a global maximum after some time and then asymptotically converges to zero. (iii) It increases for all time and attains a maximum value at the asymptotic limit $t\rightarrow \infty$. In particular a non-flat spacetime is obtained that mimics the expansion in a flat FRW universe dominated by vacuum energy such that the expansion factor has the same form in both. However, the energy density is dynamical in the former.


Introduction
Finding exact solutions of the Einstein's field equations (EFEs) is one of the important and an old problem in General Relativity (GR). The search for new solutions opens new avenues to our understanding of the universe. These equations are the result of Einstein's revolutionary idea that the existence of matter induces curvature in a spacetime. This is responsible for an inward pull of all surrounding objects, thereby completely changing the standard perspective of gravitational force [1,2]. Recent developments in astrophysics and cosmology reveal that our universe is mainly composed of a mysterious form of matter (dark energy and dark matter) which is the main cause of cosmic acceleration at a faster pace today [3,4]. This astonishing observation laid down the basis of numerous phenomenological models and led to the departure from the standard solutions of EFEs in the thrust of new exact solutions possessing revived geometrical and dynamical properties. Our attempt in this direction is to exploit the geometrical features in homogeneous cosmological solutions of these equations which mainly arise from Noether symmetries.
Suppose M is a four-dimensional smooth manifold equipped with a non-degenerate metric g of Lorentzian signature (+, −, −, −). The EFEs of the universe in standard gravitational 1 units c = 1 = G, are a system of partial differential equations G µν = T µν (κ = 1), where G µν is the Einstein tensor which contains all the basic geometric properties of spacetimes and T µν is the stress-energy tensor that describes the density and flux of energy and momentum in spacetimes [2]. The exact solutions of EFEs either arise from the geometrical consideration of spacetimes or by studying physical characteristics of matter or a combination of both. In cosmology, homogeneous spacetimes of dimension 1 + 3 are characterized using Bianchi classification of real Lie algebras of dimension three which resulted into nine classes I, II, ...,IX [5]. The famous cosmological solutions of EFEs known as Friedmann-Lematre-Robertson-Walker (FRW) metrics are both isotropic and homogeneous and are particular cases of Bianchi types I, V, VII and IX [1]. Therefore it is important to investigate other Bianchi homogeneous cosmologies with regard to their invariant geometrical properties. Indeed the authors [6], have done the symmetry analysis of all Bianchi spacetimes in the presence of a dynamical field. In [7], authors have classified Bianchi V spacetimes using symmetries of the curvature and Weyl tensors known as curvature collineations (CC) and Weyl collineations (WC), respectively.
The concept of symmetries is closely related to the invariance of geometrical quantities under certain diffeomorphism on the manifold that could change the manifold however it keeps the underlying structure intact. This is a well known problem of equivalence of the two geometric objects under diffeomorphisms [8,9]. The invariance of geometrical quantities which are described by tensor fields on the manifold under certain diffeomorphisms is ensured if the Lie derivative vanishes along the vector field representing a flow. Suppose T is a tensor field of any type on the manifold M. The key to check the invariance of certain geometrical quantities is to use the definition of Lie derivatives and require that £ X T = 0, where £ X is the Lie derivative along the vector field X [5]. In case of the metric tensor the symmetries are known as conformal Killing vectors (CKV) provided £ X g µν = 2Ψg µν , where Ψ(x µ ) is known as a conformal factor. If Ψ ,µν = 0, then a CKV is said to be proper otherwise it reduces to a special CKV if Ψ ,µν = 0 and Ψ ,µ = 0. The homotheties (HV) and Killing vectors (KV) arise if Ψ ,µ = 0 and Ψ = 0, respectively. Furthermore the symmetries of the curvature tensor are obtained by replacing the tensor field with curvature tensor £ X R µ νρσ = 0. In order to investigate pure gravitational field of spacetimes Weyl tensor plays a significant role as it is conformally invariant and its symmetries are obtained by the same condition that the Lie derivative of the Weyl tensor C µ νρσ vanishes.
Noether symmetries play an essential role in finding conservation laws of the equations of motion with the use of Noether theorem [8,10]. Sometimes it is difficult to obtain exact solutions of the equations of motion, however the dynamics can be reduced by investigating the invariant properties of the system provided the problem under consideration is variational and there exist a Lagrangian. Noether's theorem provides an explicit formula of a conserved quantity for each continuous symmetry transformation that leaves the action invariant [11]. In addition Noether symmetry analysis reduces the problem by specifying the unknown functions that appear in the Lagrangian. This line of approach has been followed by several authors, notably in the pioneering works of S. Capozziello [12,13], G. E. Prince [14,15] and M. Tsamparlis [16,17]. The connection between symmetries of the underlying manifold and that of the differential equations was discussed in [18]. The classification of spherically symmetric static spacetimes via Noether symmetries is done in [19]. Besides approximate Lie symmetries was used to resolve the problem of energy in general relativity in [20,21]. In this paper our main interest is the classification of Bianchi V spacetimes using Noether symmetries. These symmetries provide crucial physical information about the conserved quantities of a physical system.
We start our investigation by considering the geodesic motion in Bianchi type V spacetimes.
In such spacetimes there are three arbitrary functions whose specification arise from the presence of Noether symmetries. We give a complete classification of Noether algebras along with first integrals for each case. It is found that Noether algebra of the symmetries of the geodesic Lagrangian in such spacetimes can have dimension, 4 − 7, 9 or 10. The algbera of Lie point symmetries of the equations of motion is of dimension, 5 − 8, 10, 12 or 13. The connection between Lie and Noether symmetries with symmetries of the spacetimes like HVs, projective collineations (PCs) and CCs has already been discussed in [18,16]. We highlight the important features of new cosmological solutions in the light of above results.
It is well known that there are three cosmological models of our universe which are filled with vacuum energy, radiation and matter [1]. The large astrophysical data suggest a widely accepted view that our universe is nearly flat. It would be interesting to investigate as to what degree we can relax the condition of flatness of the underlying manifold such that the non-flat spacetime still carry all the important features of a realistic and viable cosmological model. It turns out that using Noether symmetries we can specify certain non-flat spacetimes which possess very nice physical properties. In particular we determine a non-flat spacetime that mimics the behavior of a flat spacetime filled with vacuum energy. Interestingly the expansion factor (so is the Hubble parameter) in the non-flat and flat spacetimes is same. The energy density is fixed in the flat model (filed with vacuum energy) and dynamical in the other case. Based on our finding we further investigate non-flat spacetimes and compare these with flat spacetimes dominated by matter (Einstein-de Sitter universe) or radiation.
The paper is subsequently divided into three main parts. The Noether symmetry analysis is presented in the next Section where we give preliminaries to the symmetry approach. In the third Section we investigate the physical characteristics of such spacetimes and examine the implications of positive energy condition on the solutions. To confront our results with other extended theories of gravity we consider a particular Bianchi V spacetime which meets the positive energy condition and modify our results in f (R)−gravity. It is found that f (R) ∝ R 3/2 , for a dust cloud in an anisotropic Bianchi V spacetime. Lastly we summarize our results in the last Section. 3

Noether symmetry Analysis
The Bianchi type V spacetimes have the form where A, B, C, are three non-zero arbitrary functions of the cosmic time t. The expansion factor α, is a non-zero constant which has the units of inverse length. The above spacetimes represent a specially homogeneous and anisotropic cosmologies as the coefficients depend on the time variable t. Since these are non-static spacetimes and thus do not admit a time like Killing vector field. All Binachi V spacetimes admit a three-dimensional Lie algebra of spacelike KVs (corresponding to the isometry group G 3 , that acts transitively on the spacelike hypersurfaces) [5] Subsequently we identify this algebra with K 3 = {X 1 , X 2 , X 3 }, corresponding to the basic isometry algebra of the underlying Bianchi V spacetime.
An action A = L ds, of the geodesic motion possesses the Lagrangian L = L(s, x µ ,ẋ µ ), in the background of Bianchi V spacetimes, equivalently where an over dot represents the derivative with respect to the geodetic parameter s. A Noether symmetry is a vector field 1 which leaves the action invariant such that it satisfies the condition where ξ = ξ(s, x µ ), η µ = η µ (s, x µ ), (µ = 0, .., 3) and G(s, x µ ) is an unknown function of the corresponding symmetry. The first integrals of the equations of motion are determined by famous Noether theorem which gives an explicit formula [11] where x µ = (t, x, y, z) and η µ = (η 0 , η 1 , η 2 , η 3 ) denote the coefficients of the Noether symmetry (3). Note that D s I = 0, upon using the equations of motion, where 1 Note that the action of the underlying Lie group is extended to the product manifold M × R, which is (4 + 1)−dimensional so is to bring the dynamical symmetries (Lie and Noether) and geometrical symmetries (CKV,CC,PC etc) on equal footing. The action naturally includes the parametrization with respect to the geodetic parameter s, where this line of approach was followed in [14].
The equations of motion (geodesic equations) for the geodesic Lagrangian in Bianchi V spacetimes comprise of a system of four second-order ordinary differential equations in the field variablesẗ + e αx (BB ′ẏ2 + CC ′ż2 ) + AA ′ẋ2 = 0, where ′ represents the derivative with respect to t. It is known that Noether algebra is a subalgebra of the Lie algebra of point symmetries of the differential equations [8], therefore our study yields classification of Noether algebras of the equations of motion for the geodesic Lagrangian in Bianchi V spacetimes. In order to find Lie point symmetries of the differential equations we extend the vector field (3) to the jet space of second-order where we require that where E, is a given system of differential equations and X [2] is the second-order prolongation which is defined below where [1] η µ and [2] η µ are determined from the formulas [1] [2] η µ = d [1] η µ ds −ẍ µ dξ ds .
The Lie symmetries are obtained by applying the operator (9) on the geodesic equations while replacing the second order derivative terms in the last equation from the equations of motion (7).
Since all Bianchi V spacetimes admit K 3 therefore the minimum dimension of the Noether algebra of the geodesic Lagrangian can be obtained easily. Indeed it is easy to prove that the minimal Noether algebra is N 4 = K 3 ⊕ {∂ s }, which is of dimension four. Proof. The Lie algebra of isometries, K 3 , leaves the action invariant therefore all KVs are Noether symmetries of the Lagrangian (2). It is easy to see that the global action of full isometry group G 3 on the Lagrangian (2) is L(s, t, x, y, z) −→ L(s, t, x + 2ǫ 3 , e −αǫ 3 y + ǫ 1 , e −αǫ 3 z + ǫ 2 ), which clearly leaves the action invariant. Therefore all KVs are Noether symmetries of the Lagrangian (2). Finally Lagrangian (2) does not explicitly depend on s therefore invariance under s−translation is trivial, i.e., L(s, t, x, y, z) −→ L(s + ǫ 4 , t, x, y, z).
Theorem 2. The minimum dimension of algebra of Lie point symmetries of the geodesic equations (7) in Bianchi V spacetimes is 5.
Proof. Since Noether algebra is a subalgebra of Lie symmetry algebra therefore N 4 is a Lie subalgebra. Besides the equations of motion (7) are invariant under scaling symmetry s∂ s , therefore for all Bianchi V spacetimes the minimum dimension of Lie algebra is five which is identified as N 4 ⊕ {s∂ s }.
In order to identify Bianchi V spacetimes that possess minimal Noether algebra N 4 of dimension four and other algebras we solve the determining equations for Noether symmetries which are obtained by splitting equation (4) in monomials and obtain a set of nineteen linear partial differential equations The above system can be integrated for various forms of the arbitrary functions. We use CAS (Maple) to categorize all the cases using the command 'rifsimp'. It is interesting to see that in most cases the function A(t) turns out to be an affine function. Therefore the above system of PDEs uniquely characterize the affine form of A(t) which specify the forms of other arbitrary functions. We denote Noether algebra with N i , where i refers to its dimensions and the distinction of subcases is made through another subscript. We observe that the following cases arise and in each case we also obtain the first integrals of the equations of motion (7) using Noether theorem (5).

Case 1. (4-dimensional algebra)
It is identified as N 4 , which arise for the geodesic Lagrangian in Bianchi V spacetimes where evolutionary functions are arbitrary and not among those that arise in the subsequent cases. The four independent invariants of the equations of motion corresponding to N 4 , include three invariants for algebra K 3 and fourth invariant is for ∂ s , which is the Lagrangian itself, or equivalently The Lie algebra of point symmetries of the geodesic equations includes all Noether symmetries and an additional scaling symmetry L 1 = s∂ s , which is not a Noether symmetry. Here we use L, to denote a Lie symmetry which is not a Noether symmetry with respect to a geodesic Lagrangian. Therefore, the Lie algebra is N 4 ⊕ {L 1 } and is of dimensions five.
There are three cases of five dimensional Noether algebras which are possessed by the Lagrangian in Bianchi V spacetimes with where c i , ∀ i are non-zero constants. In the first case both powers m and n can not be equal to zero or one, simultaneously. Similarly the powers β and γ can not be zero, simultaneously in case (2c). In the first case we obtain two Noether algebras Besides invariants I 1 − I 4 , an additional invariant for each case is however for other two cases the invariants include The algebra of Lie point symmetries is six in all cases where the identification is N 5, Case 3. (6-dimensional algebra) We obtain two cases in which the geodesic Lagrangian admits a 6−dimensional Noether 7 algebra. These are identified as up to a redefinition of constant in case 2, which are possessed by the Lagrangians in Bianchi V spacetimes respectively. The invariants I 1 − I 5 are same as before with an extra invariant We find that there are seven Lie symmetries for the geodesic equations where we get a decomposition of the extra vector fields obtained in each case. In particular, the vector field 2αs∂ s + αt∂ t + 2(1 − n)∂ x , divides into two independent Lie symmetries such that there linear combination is a Noether symmetry but independently these are not Noether symmetries. It is easy to check that the algebra N 5,b ⊕ {L 1,a , L 2,a }, is a sevendimensional Lie algebra of point symmetries for equations of motion in case (3a). For the other case (3b) the Lie algebra of point symmetries is Case 4. (7-dimensional algebra) The geodesic Lagrangian with all evolutionary functions being equal admits N 5,b algebra along with two additional symmetries The two new invariants are The Lie algebra of point symmetries is N 7 ⊕ {L 1 = s∂ s }, which has dimension eight.
Case 5. (9-dimensional algebra) In this case there appear three subcases in which all evolutionary functions are specified. In the first two subcases, N 7 is a subalgebra of N 9 , however in the last subcase only N 5,b is a subalgebra.
5a. It is identified as N 9,a = N 7 ⊕ {X 8 , X 9 }, where evolutionary functions are specified by Additional symmetries include where the gauge term is constant in all cases except for X 9 , which is obtained with nonconstant gauge G = 2t. The corresponding conserved quantities are The Lie symmetry algebra of the corresponding geodesic equations is of twelve dimension which is identified as Additional symmetries include with invariants Now in this case the algebra of Lie point symmetries is only ten dimensional unlike the previous case (5a), where the extra symmetry generator arise from the decomposition of X 9 , i.e.
In this case we again obtain twelve dimensional Lie algebra of the geodesic equations which is Case 6. (10-dimensional algebra) There are two subcases of ten dimensional algebras of Noether symmetries corresponding to α = 2c 1 β and α = −2c 1 β, respectively.
6a. It is identified as N 9,c ⊕ A N 1 , with where N 9 is the same as in case (5c) by substituting the value α = 2c 1 β. An additional symmetry A N 1 that arises is with corresponding invariant 2 ) − c 2 2 yẏ − c 2 3 zż e 2β(t+c 1 x) . (47) In this case we obtain thirteen dimensional Lie algebra of the geodesic equations which is 6b. It is identified as where N 9 is the same as in case (5c) by substituting the value α = −2c 1 β. An additional symmetry A N 1 that arise is with corresponding invariant In this case we again obtain thirteen dimensional Lie algebra of the geodesic equations which is This completes the classification of Noether symmetries of the geodesic Lagrangian in Bianchi V spacetimes. The equations of motion (7)    The summary of our results is given in Table 1, where d(N ) and d(L) refers to the dimensions of Noether and Lie algebras, respectively.

Physical Interpretation of New Solutions
The study of inhomogeneous and anisotropic cosmologies started soon after the birth of GR [22] and a detailed analysis of exact solutions in terms of the asymptotic of singularities appeared in the works [23,24] followed by numerous attempts [25,26,27,28]. For a more detailed survey on such cosmologies the reader is referred to [29]. The class of Bianchi V spacetimes contains anisotropic and homogeneous cosmologies which are crucial to investigate for several reasons. The universe is homogeneous at very large scales and the question of whether it had started with a little bit irregularity require us to analyze small perturbations away from the high symmetry of Friedmann models. This could help us investigate on the present-day anisotropy of the microwave background radiation and irregularities in density and temperature at the early epochs when the radiation was emitted. Besides adiabatic cooling, viscous dissipation and particle creation of an anisotropic universe can be studied with new cosmological models.
The information of matter content in Bianchi V spacetimes is contained in the Riemann and Weyl tensors which can be used to examine the effect of tidal forces due to curvature in the manifold along the geodesics. The former helps us to monitor the change in volume while later provides change in shape of the observer along geodesics. The basic requirement for a spacetime to be physical is that the positive energy condition is met. The positive energy condition requires that T 00 which corresponds to the energy density is non-negative therefore T 00 ≥ 0. It is worth pointing out that the existence of such Noether symmetries completely specify the cosmological models which is in agreement with the results in [30], where pointlike Noether symmetries were employed to determine feasible models in extended gravity quantum cosmology. Our prime interest here is to investigate the positive energy condition in the models specified by Noether symmetries.
The components of Einstein tensor, Weyl and curvature tensors involving arbitrary functions are already given in [7]. It was also shown that the rank of a 6 × 6 curvature matrix is 3, 4, 5 or 6 in Bianchi V spacetimes where the case of rank 3 give rise to infinite dimensional Lie algebra of proper CCs. On the other hand the rank of 6 × 6 Weyl matrix is 0, 4 or 6 and such spacetimes do not admit proper WCs except for the trivial rank zero case. The Lie algebra of proper conformal KVs in such spacetimes is four-dimensional. We now discuss the physical interpretation of the cosmological solutions obtained in the last section in the light of positive energy condition.

Case 1. (4−dimensional algebra)
The class of Bianchi V spacetimes admitting 4-dimensional algebra of Noether symmetries contains all three arbitrary scale factors and the corresponding invariants are given in (12). We investigate general physical characteristics of this spacetime by considering the energy density where the positive energy condition requires that ρ(t) ≥ 0. It is convenient to bring ρ(t) into a more useful form by introducing a simple change of quantities a = ln A, b = ln B, c = ln C. The first three terms in the above equation does have the qualitative behavior of kinetic energies in Newtonian mechanics [1]. Since a ′ , encodes the information how fast or slow the expansion takes place, therefore the mixed term a ′ b ′ , can be regarded as the kinetic energy of the composite system of both a and b. Similar argument holds for the composite systems (b,c) 13 and (a,c). However it is important to note that in the above relation we do not have the contribution of individual kinetic energies due to a, b and c, respectively. In fact the product terms can be set equivalently to where the first term is positive and corresponds to the kinetic energy T c , of the composite system a, b and c. The second term is the total kinetic energy T e , of individual systems which is negative. Therefore in all Bianchi V spacetimes the total kinetic energy of individual and composite systems is irrelevant and the quantity which is crucial is the difference T c − T e . For a realistic model we require that T c > T e . The energy density is negative in the case when T c = T e , where the spacetime is unrealistic. The last term could be regarded as the potential term which is an exponential function of the scale factor a. The significance of the first scale factor A(t), upon the others is apparent and notably Noether symmetries also characterized it.
Note that the contribution of the last term is small compared to other terms as long as |A(t)| grows with time. However the evolution of A(t) is critical in the interval A(t) ∈ [0, 1] , in which case a(t) is negative and the last term plays a significant contribution in decreasing the overall energy density of the spacetime. On the other hand the slopes of b(t) and c(t) play a significant role on the energy density of spacetimes. In particular from the equation (51), it is clear that the contribution of the product b ′ c ′ , is larger than the sum b ′ + c ′ . Therefore if both b and c, simultaneously accelerate or decelerate then both slopes are positive or negative and we expect the spacetime to be realistic. Moreover if the slopes are opposite then the only possibility for a physical spacetime is that when A(t) ∈ [0, 1] , in which case the last term contribute in a positive energy density so is to balance the effect of other terms. Now we consider those cases in which the evolutionary functions are specified by the existence of Noether symmetries and identify the critical bounds on the expansion parameters for realistic Bianchi V spacetimes.

Case 2. (5-dimensional algebra)
In this case we obtained three subcases in which one of the subcase contains arbitrary functions. We consider the cases in which all functions are completely specified which has non-zero components of the Weyl tensor therefore it is not conformally flat spacetime. There are eleven non-zero curvature invariants where we mention only Ricci scalar and denote it by I R , given by 14 which is singular at t = 0. The Einstein tensor is given by where (56) For a realistic cosmological model, we impose positive energy condition on the dynamical energy density which is positive if 4c 2 1 (m + mn + n) − 3α 2 ≥ 0, therefore we obtain a critical bound on the evolution factor α Since the quantity in the square-root must be positive therefore we obtain an extra condition on the powers m and n, i.e., m > −n/(n + 1), n = −1. If n is a positive number then there are two possibilities that either m is positive or negative. If it is negative then the spacetime becomes singular at t = 0, while it is non-singular in the other case. For −1 < n < 0, the evolutionary function B(t) = c 1 t m , has a positive power. However in the other case n < −1, the evolutionary function can have a positive or negative power as before. Therefore the evolution of above spacetime is such that it started at an initial time t = c 4 = 0, then the energy density continue to decrease and vanishes as t → ∞. The flux across x and y surfaces is zero therefore the density of first two components of linear momentum is zero however the z− component of linear momentum density is which could be positive or negative depending on the choice of α, m and n and it asymptotically decays as t → ∞.
Similarly in the other case the spacetime has the form which also has non-zero components of the Weyl tensor therefore it is not conformally flat.
Again there are eleven non-zero curvature invariants and the Ricci scalar is which is non-singular and non-dynamical unlike the previous case. The matter tensor is given by In order to obtain a bound on the coefficient α, we consider the density which is positive if Now there are two cases that either β and γ are both positive or negative. Unlike the previous case the above spacetime has a fixed energy density for all time. On the other hand the momentum density is also fixed which could be positive if both β and γ are positive and negative otherwise.
The spacetime is where the Weyl tensor vanishes therefore it is conformally flat and Petrov type O. There are four curvature invariants where the Ricci scalar is The matter tensor becomes therefore the energy density is positive G 00 = ρ(t) > 0, for the following critical bound on α where the above spacetime is realistic. The quantity n(n + 2) must be non-negative which is true if n > 0 or n < −2. The energy density vanishes as t → ∞, so is the momentum density for the above spacetime.
In the other subcase the spacetime is where the Weyl tensor vanishes therefore it is conformally flat and Petrov type O. The Einstein tensor becomes therefore the energy density is negative G 00 = ρ < 0, for all time and the above spacetime is unrealistic.

Case 4. (7-dimensional algebra)
In this case the spacetime involve one arbitrary function which turns out to be the case where Weyl tensor is zero thus the above spacetime is conformally flat and is of the Petrov type O. This is an interesting case as it can be regarded as an inhomeogeneuos extension of the FRW spactime. There are four curvature invariants including the basic Ricci scalar The Einstein tensor assumes the form where there is no non-zero component of the linear momentum density. The energy density which could be positive or negative, in general. We now consider two dynamical evolutions of this non-flat spacetime in terms of a power-law and an exponential law. In case of a power-law A(t) = c 1 t m , the density ρ(t) → ∞ for m < 0 and ρ(t) → 0 for m > 0, as t → ∞, where the former is not an interesting case while the later contains some physical information. For example the density takes the form and its behavior could be examined by choosing an initial point where it is zero and it was negative before. Therefore we obtain three different cases (i) 0 < m < 1 (ii) m = 1 and (iii) m > 1.
In the first case we choose m = 1/2, without loss of generality and obtain which is zero at t = c 2 1 /α 2 > 0. The critical point of ρ(t) is t c = 2c 2 1 /α 2 , where the density attains a global minimum because ρ ′′ (t c ) > 0. Since ρ(t c ) < 0, is negative therefore it is not a realistic spacetime. For the second case m = 1, the behavior of energy density is such that it starts from an initial value and decreases continuously till it vanishes and the global maximum is the initial point where the energy density has started at some non-zero time. The graph of energy density is depicted in Figure 1. In the last case m > 1, it attains a global maximum value at after which it continue to decrease till it vanishes as t → ∞, as is shown in Figure 2. By introducing a Hubble type parameter responsible for the possible expansion or contraction of this spacetime we find that therefore H(t) ∝ t −1 , and this model mimics the behavior of a flat matter dominated universe (Einstein-de Sitter universe). Note that in which the expansion function carry the same form for m = 2/3, in the above case. It is interesting that our model is not flat and carries an energy density which is different from the energy density in a flat Einstein-de Sitter universe which varies as ρ m (t) ∝ t −2 . In our non-flat spacetime the energy density varies ρ(t) ∝ (4/3 − 3α 2 t 2/3 /4c 2 1 )t −2 in the case m = 2/3, where the convergence of ρ(t) is faster in comparison to ρ m . In the case of radiation dominated universe the energy density also varies in the same way. However it is constant for a flat spacetime dominated only by vacuum energy. Note that in the case of radiation dominated flat universe the expansion factor is proportional to t 1/2 , that lies in our first case. The energy density varies in the same way ρ γ ∝ t −2 , which in the above non-flat spacetime is ρ(t) ∝ (1 − α 2 t/c 2 1 )t −2 .
In the case of an exponential form A(t) = c 1 e βt , the energy density becomes which has the positive initial value ρ(0) ≥ 0, if α ≤ √ 3c 1 β. The case in which β < 0, represents a contracting non-flat spacetime that carries a positive energy for a short time and continue to decrease such that ρ → −∞ as t → ∞, because the contribution of the dynamical term will dominate the first term at some point. On the other hand an expanding non-flat spacetime requires that β > 0, in which case we observe that the energy density increases rapidly for a short time where after it remains constant such that at an asymptotic limit t → ∞, it attains a finite positive value 3β 2 /4. The graph of energy density for this case is given in Figure 3. It is unexpected because in an expanding flat spacetime we expect that the energy density decreases indefinitely, however the curvature in this spacetime confine it to attain a positive definite value at an asymptotic limit. The Hubble parameter for this spacetime is non-dynamical H(t) = β > 0, therefore this non-flat spacetime mimics a flat universe dominated by vacuum energy which has fixed density. However this non-flat spacetime has a varying energy density such that it attains a fixed positive value after a short interval of time.
19 Figure 3. The graph of energy density for α = 1, β = 1, c 1 = 1 The curvature rank for the spacetimes in case (5b) is 3 and the Lie algebra of both curvature and Weyl collineations is infinite dimensional [7]. In the first case we obtain following nonzero components of the Einstein tensor (81) Since the energy density is negative therefore the above spacetime is unrealistic although it is conformally flat. The components of the matter tensor in case (5b) can be obtained directly by placing the value of A(t) = c 1 t + c 2 , in case 4, equivalently The positive energy condition requires that |α| < 2|c 1 |, which indicates that the expansion coefficient α must be smaller than the twice the slope of graph of A(t) = c 1 t+ c 2 . The energy density vanishes as t → ∞.
Lastly the case (5c) has a distinct difference from the previous two cases that there arise non-zero off-diagonal terms in the matter tensor. Following are the non-zero components of 20 the Einstein tensor where we must have |α| < 2/ √ 3 |c 1 β|, for a realistic model. Note that the energy density in this Bianchi V spacetime is fixed and it attains non-zero components of the momentum density unlike the previous cases where these were zero.
Case 6. (10-dimensional algebra) There are two subcases of 10-dimensional Noether algebra which are Although the above models are conformally flat but these are unrealistic spacetimes as the Einstein tensor is and the energy density in both sub cases is negative, G 00 = −2β 2 < 0.

Some Considerations into f (R)−Gravity Theory
We now develop a few considerations for the anisotropic cosmological solutions that are obtained as a consequence of the existence of Noether symmetries in the context of f (R)−gravity.
The idea is to probe into earlier epochs in the evolution of the Universe using the results obtained in the previous section. As there arises several cases therefore we confine ourself to a few most relevant cases for our purpose. For example let us consider the case (5b) which meets the positive energy condition for (|α| < 2|c 1 |) given by in which case the Ricci scalar is which remains negative for all time due to positive energy condition therefore the underlying manifold is hyperbolic. We now obtain the field equations using the action of f (R)−gravity involving a matter term in standard gravitational units [31,32] where κ = 1, g denotes the background metric and f (R) is a general function of first curvature invariant which in our case is given by (86). Here ψ denotes all matter fields coupled to gravity. The field equations now assume the form where f ′ = df (R)/dR, ∇ µ is the covariant derivative with respect to the Levi-Civita connection of the metric and = ∇ µ ∇ µ . An immediate consequence of the above field equations is that their trace yield an important equation which indicates how matter part T , in f (R)-gravity is differentially connected to the curvature in the given spacetime unlike in standard GR, where we have an algebraic relationship R = −κT . We now employ Bianchi V spacetime (85) and solve the fields equations in a relatively simple case of a dust cloud. As in our case the Ricci scalar R is a function of t, therefore the field equations yield a system of ordinary differential equations f ′ (R(t)) = F (t), where we subsequently represent the time derivative with an overdot (not to be confused with an overdot in Section 1). In this case the non-zero components of Ricci tensor are which upon using the field equations (87) result into We now assume a perfect fluid whose energy-momentum tensor is given by where ρ(t) and p(t) are energy and pressure densities of the fluid which satisfy the equation while u µ = √ g 00 (1, 0, 0, 0), is the four velocity in comoving coordinates. The conservation equation yieldsρ which for the underlying spacetime (85) and in case of a pressure-less fluid ω = 0, assumes the formρ + 3c 1 ρ c 1 t + c 2 = 0, 22 whose solution involving a constant of integration C 1 is By substituting the above value in (91) and (92), we obtain an ordinary differential equation F (t), given byF which has a solution where C 2 and C 3 are constants of integration. The exponents λ + and λ − , are given by which are both real for positive energy condition. In order to find the explicit form of f (R), we employ the equation (89) to obtain which in terms of R, using (86) is given by This is quite unexpected because the above relationship is found for f (R)−gravity in the background of FRW universe and it is observed that all other models other than this remain obscure [31]. Therefore, an analysis of a dust cloud in the background of an anisotropic and homogeneous Bianchi V spacetime (85) which does not correspond to an FRW universe, reveal the same form of f (R). Here we have demonstrated a successful implementation of our results obtained in the previous section in f (R)−gravity. It would be interesting to investigate more physical models with different equations of state in the background of other Bianchi V spacetimes following the same lines. A detailed analysis of rest of the cases of such spacetimes and a brief comparison with other extended theories of gravity will be discussed elsewhere.

Summary
We used Noether symmetries to study Bianchi V cosmologies. It is identified that the algebra of Lie point symmetries of the geodesic equations in such spacetimes can have dimensions 5 − 8, 10, 12 or 13. On the other hand the dimensions of Noether algebras is 4 − 7, 9 or 10. The presence of Noether symmetries helps us to obtain first integrals or constant of motion using Noether theorem corresponding to the geodesic Lagrangian. A key advantage of this approach is the determination of unknown functions in the Lagrangian which are specified while solving the determining equations for Noether symmetries that reduces the dynamics significantly. In all cases 2 − 6, we explicitly specified the unknown evolutionary functions in these spacetimes. Besides one can use first integrals to obtain closed form of exact solutions of the equations of motion, i.e., geodesics in non-flat spacetimes.
In order to interpret our results physically we used positive energy condition which holds for all realistic cosmological models. Interestingly it imposes constraints on the solutions and provide us critical bounds on the constants for which a Bianchi V spacetime is physical. The expansion coefficient α, in all Bianchi V spacetimes is crucial. By imposing the positive energy condition on the obtained solutions it turns out that we can investigate different possible critical bounds on α, in which the non-flat spacetimes carry interesting physical characteristics. The positive energy condition further imposes constraints on the spacetimes and we obtain cases in which the energy density behaves in one of the following ways. It is positive and constant for all time. It varies with time and attains a global maximum after sometime where after it asymptotically converges to a relatively smaller but positive value. It increases for all time and attains a maximum value at the asymptotic limit t → ∞. A brief comparison of flat models (vacuum, radiation, matter) with the non-flat spacetimes is given in detail.
Another important consequence of this study is the identification of f (R) ∝ R 3/2 , for a dust cloud in an anisotropic and homogeneous spacetime (85) which is proven to exist in FRW universe using standard f (R)−gravity approach. It indicates us to briefly probe into earlier epochs of our Universe when it was least isotropic. Therefore those Bianchi V spacetimes where the positive energy condition successfully holds are the best candidate for further investigation with more physical models with different equations of state.