Several observations such as Super Novae type Ia  -  , Cosmic Microwave Background Radiation  , Baryon Acoustic Oscillations    and PLANK collaborations  indicate that our universe is in the phase of an accelerating expansion. It is believed that the reason of the late time acceleration is Dark Energy (DE). However, it is still not completely known what the DE is. According to the abovementioned observations, there is an important known fact that ~70% rate of the energy content of the universe is DE. Scientists generally try to explain dark energy in two ways. One is to adopt some exotic matter sources. The other is to modify the original field equations of Einstein’s general theory of relativity. The first and the simple modification of Einstein’s field equations (EFEs) as a DE candidate is the cosmological constant Λ. Nevertheless, the cosmological constant is suffered by fine tuning problem and cosmic coincidence problem. Except the cosmological constant Λ, there are several modifications of EFEs. Besides others, the modified gravity theory has become very popular for last ten years  .
Many researchers have considered many problems of the cosmology in this theory. Adhav  has obtained the solutions of locally rotationally symmetric (LRS) Bianchi type-I universe in theory. Houndjo et al.  have constructed a cosmological scenario in gravity and discoursed a matter dominated era transition. Taking a variable deceleration parameter that depends on the Hubble parameter, Tiwari et al. considered the LRS Bianchi type-I model in theory  . Sofuoğlu showed that the Gödel universe is a solution of the field equations in theory  . In theory of gravity, Tiwari et al. examined Bianchi type-I model in which the cosmological term is decaying  . Recently, Tiwari et al.  have examined the gravity using a time dependent cosmological term. Tiwari et al.  have explained phase transition in expansion of LRS Bianchi type-I universe in theory. Tiwari et al.  have considered a quadratically varying deceleration parameter in this modified theory.
In this work, adopting a periodically varying deceleration parameter, we investigate LRS Bianchi-type I universe in gravity theory by adopting a particular form of function as , where is a constant. For finding an exact solution of the field equations of the model, the model presented is based on a unique condition of periodically time varying deceleration parameter.
The outline of this paper as follows: Basic formalism of theory is given in Section 2; the solutions of the field equations for LRS Bianchi-type I universe are obtained and the physical and geometrical characterization of the model is represented in Section 3; and the conclusions are given in Section 4.
Over the study, we use the natural units as .
2. Basic Formalism of Theory
The action of theory is described by 
Here R is the Ricci scalar and T is the trace of the matter energy-momentum tensor , is an arbitrary function of these two scalars, is the Lagrange density of matter. is defined as
where is the covariant components of the metric tensor and g is the determinant of the metric tensor. The variation of the gravitational action (1) with respect to , yields the field equations of this modified theory:
where , and is the four-dimensional covariant derivative, is defined in Equation (2) and .
Here we note that if one takes then Equation (3) give the field equations of theory.
The contraction of Equation (3) gives
Now, being any function of T, we take as
By putting Equation (5) in Equation (3), we obtain the following field equations of gravity
where is the Einstein tensor and .
For a perfect fluid source, may be written as
Here is the matter-energy density of the fluid, p is isotropic pressure, is the four-velocity vector of the observer. In this case the field Equation (6) become  
3. The Model and Solutions
The LRS Bianchi-Type-I universe is defined by the following ansatz
where A and B are time-dependent scale factors.
Now, being a constant, choosing the function as
Equation (8) gives the following field equations for Equation (9)
Here a dot (⋅) represents a derivative with respect to time t.
Equations (11)-(13) define a system of equations which has four unknowns (A, B, ρ, p) with three independent equations. In order to find the solution such a system, one more relation is required. Hence, we carry out a law of variation of deceleration parameter (DP). The time varying DP is important in evolution of the universe. Its phase transition in expansion may be well explained by the time varying DP. Now, we adopt the following periodic time varying DP 
Here m and k are positive real numbers. Figure 1 displays the behavior of DP with time for different values of the constants. One can observe from Figure 1 that the model starts with a decelerating phase of expansion ( ) and later transits to an accelerating phase of expansion ( ) and then goes to super-exponential phase of expansion ( ) in a cyclic history.
Using the definition of DP as , the integration of Equation (14) gives the Hubble parameter H as
Here is a constant of integration. Here we may choose and then Hubble parameter becomes
Figure 1. DP vs time in the units of Gyrs.
Using the definition of Hubble parameter as in Equation (16), the average scale factors a is obtained as
where a0 is a constant of integration.
Figures 1-3 show that evolution of the deceleration parameter, the Hubble
Figure 2. Hubble parameter vs time.
Figure 3. Scale factor vs time.
parameter and the scale factor in time with the unşts of gigayears.
The spatial volume is given by
where a is the average scale factor.
The average Hubble parameter H of LRS Bianchi Type-I universe can be written as
The directional Hubble parameters in x, y and z directions are defined as, respectively,
and , (20)
For our model, the directional Hubble parameters are obtained as follows:
The anisotropisation in expansion of the model is given by the parameter Δ which is defined and found as
where i runs from 1 to 3.
The expansion scalar is
The shear scalar is defined as and found as
Using Equations (11)-(13), we obtain the pressure p as
and we obtain the energy density as
We see that the average scale factor is zero at initially. It increases in cosmic time and changes periodically. The metric potentials are vanish initially it means our model has point type singularity. Also, the singularities occur periodically at
( ). All the cosmological parameters ρ, p, θ, σ and Δ are
infinite initially and they preserve their periodic behavior in time. It is interesting that the density and pressure have large values initially and decrease to a minimum value and then again increase with the evolution of time. So all the quantities are infinite initially and they preserve their periodic behavior against the cosmic time.
In the present study, the spatially homogeneous and anisotropic LRS Bianchi type-I universe in modified theory of gravity has been investigated. The gravitational field equations are derived for the function , for a periodically variable DP as a function of cosine trigonometric form. We have held the exact solution of the model. We have shown that the obtained model has a point-type singularity initially and the similar singularities occur periodically. Also, we have explained and discussed the kinematical and dynamical character of the model that all the quantities are infinite initially and they preserve their periodic behavior against the cosmic time.
 Riess, A.G., Filippenko, A.V., Challis, P., Clocchiatti, A., Diercks, A., Garnavich, P.M., et al. (1998) Observational Evidence from Supernovae for an Accelerating Universe and a Cosmological Constant. The Astronomical Journal, 116, 1009-1038.
 Perlmutter, S., Aldering, G., Goldhaber, G., Knop, R.A., Nugent, P., Castro, P.G., et al. (1999) Measurements of Ω and Λ from 42 High-Redshift Supernovae. The Astronomical Journal, 517, 565-586.
 Clocchiatti, A., Schmidt, B.P., Filippenko, A.V., Challis, P., Coi, A.L., Covarrubias, R., et al. (1999) Hubble Space Telescope and Ground-Based Observations of Type Ia Supernovae at Redshift 0.5: Cosmological Implications. The Astronomical Journal, 642, 1-21.
 Spergel, D.N., Bean, R., Doré, O., Nolta, M.R., Bennett, C.L., Dunkley, J., et al. (2007) Three-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Implications for Cosmology. The Astrophysical Journal Supplement Series, 170, 377-408.
 Riess, A.G., Kirshner, R.P., Schmidt, B.P., Jha, S., Challis, P., Garnavich, P.M., et al. (1999) BVRI Light Curves for 22 Type Ia Supernovae. The Astronomical Journal, 117, 707-724.
 Knop, R.A., Aldering, G., Amanullah, R., Astier, P., Blanc, G., Burns, M.S., et al. (2003) New Constraints onΩM, ΩΛ, and w from an Independent Set of 11 High-Redshift Supernovae Observed with the Hubble Space Telescope. The Astrophysical Journal, 598, 102-137.
 Spergel, D.N., Verde, L., Peiris, H.V., Komatsu, E., Nolta, M.R., Bennett, C.L., et al. (2003) First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Determination of Cosmological Parameters. The Astrophysical Journal Supplement Series, 148, 175-194.
 Nojiri, S., Odintsov, S.D. and Sami, M. (2006) Dark Energy Cosmology from Higher-Order, String-Inspired Gravity, and Its Reconstruction. Physical Review D, 74, Article ID: 046004.
 Houndjo, M.J.S. (2012) Reconstruction of f(R, T) Gravity Describing Matter Dominated and Accelerated Phases. International Journal of Modern Physics D, 21, Article ID: 1250003.
 Tiwari, R.K., Beesham, A. and Shukla, B.K. (2018) Cosmological Model with Variable Deceleration Parameter in f(R, T) Modified Gravity. International Journal of Geometric Methods in Modern Physics, 15, Article ID: 1850115.
 Tiwari, R.K., Sofuoğlu, D. and Mishra S.K. (2020) Accelerating Universe with Varying Λ in f(R, T) Theory of Gravity. New Astronomy, 83, Article No. 101476.
 Tiwari, R.K., Sofuoğlu, D. and Dubey, V.K. (2020) Phase Transition of LRS Bianchi Type-I Cosmological Model in f(R, T) Gravity. International Journal of Geometric Methods in Modern Physics, 17, Article ID: 2050187.
 Tiwari, R.K. and Sofuoğlu, D. (2020) Quadratically Varying Deceleration Parameter in f(R, T) Gravity. International Journal of Geometric Methods in Modern Physics, 17, Article ID: 2030003.
 Shen, M. and Zhao, L. (2014) Oscillating Quintom Model with Time Periodic Varying Deceleration Parameter. Chinese Physics Letters, 31, Article ID: 010401.