Cosmological viscous fluid models describing infinite time singularities in f(T) gravity

In this paper we explore the state parameter behaviour of the interacting viscous dark energy in f(T) gravity. Using constant deceleration parameter we investigate the cosmological implications of the viscosity and interaction between the dark components (energy and matter) in terms of Redshift. So doing, the viscosity and the interaction between the two fluids are parameterized by constants δ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta $$\end{document} and ξ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi $$\end{document} respectively. In the later part of the paper, we explore some bulk viscosity models describing Little Rip and Pseudo Rip future singularities within f(T) modified gravity. We obtain gravitational equations of motion for viscous dark energy coupled with dark matter. Solving these equations, we found analytic expressions for characteristic properties of these cosmological models.


Introduction
Since the accelerated phase of universe expansion has been discovered by cosmological measurements obtained from SNe Ia, WMAP, SDSS and X-ray [1][2][3][4][5][6][7][8][9][10][11][12][13][14], many cosmologists have done huge efforts to explain this strange phenomenon. An exotic fluid with a large negative pressure named dark fluid would be the origin of this accelerated expansion. So, many models have been proposed and explored as CDM model which is the simplest model named cosmological constant. Besides this model, other models have been investigated (quintessence phantom, Tachyon, holographic dark energy, K-essence and Chaplyging gas) but unfortunately most of them do not permit to achieve the goal. According to the astrophysical observations, it has been mentioned that the cosmic medium is not a perfect fluid. So, various viscous fluids have been introduced to satisfy an inhomogeneous equation of state. Many authors have explored cosmoa e-mail: docenzi@yahoo.fr (corresponding author) b e-mail: sthoundjo@yahoo.fr logical viscous fluid effects on the dynamical evolution of the universe. Brevik et al. have considered inflation produced by two coupled fluids in the Friedmann-Robertson-Walker universe. They have used an inhomogeneous equation of state for the fluid to investigate different cosmological models for describing inflation [15]. Viscous Cosmology presenting various problems has been investigated by many authors [16][17][18][19][20][21][22][23][24][25][26][27][28][29][30][31][32]. Some alternative theories like f (T ) and f (R) expecting to explain accelerated expansion of the universe, have been explored in references [33][34][35][36][37]. Possible avoidance of finitetime singularities from viscosity in f (T ) gravity have been investigated [38]. In reference [39], authors have studied the interacting cosmological viscous fluid with infinite-time singularities solutions in f (R) gravity and found interesting results. In this paper, we follow the same step but in area of modified f (T ) gravity. From the constant deceleration parameter q we construct the scale factor that is expressed in terms of redshift. We then define the EoS for the dark energy in two cases. We first do not consider the interaction between viscous dark energy and dark matter, and secondly the interacting parameter is taking into account. We also study the bulk viscosity solutions from gravitational equation established from Freedmann-Robertson-Walker equations and viscous dark energy equation of state. We found that the bulk viscosity and the thermodynamic parameter present some corrections, in the EoS for the dark energy.
The paper is organised as follow: In Sect. 2, we present the formalism of f (T ) gravity and motion equations. In Sect. 3, we investigate the behavior of the evolution of universe filled with viscous dark energy coupled or not with dark matter. The Sect. 4 is devoted to analyse some bulk viscosity models of two future singularities (Little and Pseudo Rip) in the framework of f (T ) gravity. The paper is ended by a conclusion and remarks in Sect. 4.

Motion equations in f (T ) theory
The modified theory of gravity based on the torsion scalar is the one for which the geometric part of the action is an algebraic function depending on the torsion. In the same way as in the Teleparallel gravity, the geometric elements are described using orthonormal tetrads components defined in the tangent space at each point of the manifold. In general the line element can be written as where we define the following elements Note that η i j = diag(1, −1, −1, −1) is the metric related to the Minkowskian spacetime and the {e i μ } are the components of the tetrad which satisfy the following identity The connection used in this theory is the one of Weizenbock's, defined by Once the previous connection assumed, one can then express the main geometric objects; the torsion tensor's components as which is used in the definition of the contorsion tensor as The above objects (torsion and contorsion) are used to define a new tensor S μν λ as S μν λ The torsion scalar is defined from the previous tensor and the torsion tensor as Let's write the action for the modified f (T ) theory with matter as follows where e ≡ det[e i μ ] = √ −g denotes the determinant of the tetrad, and g the determinant of the space-time metric and f (T ), the algebraic function of T torsion.
Making the functional variation of the action (9) with respect to the tetrads, we get the following equations of where T ν μ denotes the stress tensor such that ρ and p being the energy density and the pressure of ordinary content of the universe respectively and u μ , the fourvelocity such that u μ u μ = 1. Still with the considered signature (+, −, −, −), we assume the line element for the flat Robertson-Walker universe as where a(t) is the universe scalar factor. Therefore, the torsion scalar is performed gives: From (10) to (12), one gets the following gravitational field equations: Here the "dot" denotes the derivative with respect to the cosmic time t, H =ȧ(t)/a(t), the Hubble parameter and ρ and p , the matter energy-density and pressure, respectively. f T and f T T are the derivatives of f (T ) with respect of torsion T respectively. By assuming that κ 2 = 1, the above Eqs. (14), (15) are written as: Since the effective energy density and effective pressure take the form: the energy density and pressure of dark energy are defined as: 3 Dynamical evolution of the universe filled with dark fluids In Ref. [39], the authors have used viable f (R)-gravity representing a realistic scenario for the dark energy epoch. This model describes the late-time cosmological evolution of the Universe filled with dark matter. Many authors have studied this model and found interesting results [40][41][42]. In f (T ) modified gravity, power-law models have been explored and reconstructed describing the evolution scenarios of Universe [43][44][45][46]. In the following, we will consider f (T ) power-law model to investigate Universe dynamics by using constant deceleration parameter. (1−n) and h are parameters. From the decceleration parameter q = − aä a 2 , the scale factor take the form [47,48]: where q is assumed to be a constant deceleration parameter. k and l are constants of integration. In terms of redshift, this scale factor takes the form: We now assume the dark fluid equation of state in the following form [39,49]: where δ(ρ DE ) is the bulk viscosity depending on ρ DE in general way. Acording to the thermodynamical grounds and in conventional physics, δ(ρ DE ) must be positive, thus the positive sign of the entropy changes and an irreversible process would be guarantee [39]. In the following, the dark matter will be assumed to be a dust, so ω m = 0. We will also consider in the first party of this work, the general form of bulk viscosity as [53]

Non-coupled dark fluids
The dark energy continuity equation is defined as: with p DE = ω DE ρ DE . By Using (20) and taking into acount (22) and (23), we find the dark energy density as: with 24, this dark energy density in terms of redshift, takes the form: By following the same steps as previously, we from (21) to (23), write the dark energy pressure in the form: we also express the dark energy pressure in terms of redshift, that is: with Let's now plot the functions (28), and (30). One see that, the energy density increases as the deceleration parameter diminishes with increasing Redshift. For q = −1, the energy density of dark energy increases indefinitely. Note that for particular value of h = 1/2, this energy is null.
In the graph of left hand-side, the pressure increases as q decreases, but in the ones of right hand-side, the pressure evolves in the same way as the viscosity parameter (Fig. 3).
We now, from the Eqs. (26) and (27)  : that in terms of redshift leads to: The We can see that these both EoS parameters evolve in quintessence phase, then in f (T ) gravity, universe is in accelerated expansion phase with high and low redshift. In [53], the EoS parameter for viscous Dark Energy varies in unstable phantom region and there is a transition from phantom to the cosmological constant phase at late time. In f (R) modified gravity [39], as we approach null redshift, universe seems to enter phantom phase. The left hand graph, shows that, at particular point (4, −0.23), the EoS parameter of DE seems not to be dependent of deceleration parameter q.

Coupled dark fluids
In this section, we consider an interaction between the dark energy and dark matter and the corresponding continuity equations are written as: One expresses the interaction between the both dark fluids by the quantity Q (Fig. 4). In order to ensure that the second law of thermodynamics is fulfilled, the interaction parameter Q must be positive [50][51][52]. In the same way like [53], one writes the interaction parameter in the form: where ξ is an interaction constant. From (33), (34), one obtains: and in terms of redshift, it takes the form: The Eqs. (27) and (36) permit to obtain the EoS parameter for interacting dark energy as: and the other expression of this dark energy EoS parameter in terms of redshift is written as: In the left panel, one sees that with interaction parameter, the EoS of both non-viscous and viscous dark energy approach the cosmological constant region (W DEe f f = −1) for low values of redshift. But in the panel at the right handside, the EoS parameter has the same behavior as the one in [53] because a comparative analysis of these both panels shows that the parameter W DEe f f of non-viscous dark energy (δ o = 0) is only varying in quintessence region whereas with high interaction (ξ > 2) W DEe f f starts from quintessence region, crossing the phantom divided line (PDL), and evolves to phantom region. The opposite phenomenon is observed in f (R) gravity [39], where the interaction alleviates the EoS parameter of dark energy from phantom region (for high redshift) to quintessence region (for low redshift) that evolves to the phantom region.

Viscous fluids models describing future singularities
We investigate in this area, the gravitational equation of motion that gives some bulk viscosity models describing some future singularities. Here, we are interested by Little Rip and the Pseudo Rip future singularities.

Little Rip
The characteristic of Little Rip singularity is an increasing energy density ρ with time, that needs an infinite time to appear. Thus, in this case, the EoS parameter ω < −1, and ω → −1 asymptotically [54]. Let's consider the Hubble parameter increasing exponentially with time [50] γ and H o (Hubble parameter at present time) are positive constants. Assuming that the dark matter is the dust ( p=o), the Eq. (33) gives: with ρ o an integration constant. For δ < 3, it happens that ρ → 0 when t → ∞. Considering the relations (18), (25) and the Eq. (33), one finds the gravitational equation of motion for viscous dark energy coupled with dark matter as: With (41) and Eq. (21), one obtains the Hubble-dependent bulk viscosity that is expressed by: When t → ∞, then H → ∞ and δ → ∞ if h > 1 or the bulk viscosity tends to finite value if h < 1. We now, consider another expression of the Little Rip model where the Hubble parameter is a double exponentially increasing with time: with H o , D and γ are positive constants. Using the Hubble parameter (43) and by solving the gravitational equation of motion (33) for dark matter, we find: with ρ o the constant of integration (Fig. 5). Thus, we find from the gravitational equation of motion (41), for dark energy:

Pseudo Rip
In the Pseudo Rip future singularity case, the Hubble parameter tends to cosmological constant as t → ∞. Thus, a de Sitter space is approached by the universe. Let's now analyse that model throughout the analysis of the Litte Rip model. We first assume that the Hubble parameter takes the form: where H o , H 1 and λ are the positive constants, H 1 < H o with t > 0. This means that H → H o in the late-time universe. By solving the gravitational equation of motion (33) and using the Hubble parameter (49) we obtain: with ρ o an integration constant. Using the dark matter density (40) with the power-law f (T ) model (22) and solving the gravitational (41), the bulk viscosity is written as At the end of this section, we will explore a second example of cosmological model where the Hubble parameter is given by [39][40][41][42][43][44][45][46][47][48][49][50][51][52][53][54][55][56]: In the above expression, x o = √ ρ o is the present energy, x r is a finite, andγ is a positive constant. As t → 0, one have H → x o / √ 3. Note that, in the late-time universe H → x r / √ 3, and the Hubble parameter tends asymptotically to the de Sitter solution.
The continuity Eq. (33) leads to the energy density of dark matter that is: where ρ o is an integration constant. As t → ∞, one have ρ → ∞ for ξ > 3 but ρ → o for ξ < 3. Finally, the following bulk viscosity takes the form: Also we noticed that when t → ∞ with ξ < 3, one will have δ(H ) be vanished but if ξ > 3, it happens that δ(H ) → ∞. Some corrections have appeared in the bulk viscosity due to the modified gravity model and the interaction between dark energy and dark matter (Fig. 6).

Conclusion
We have investigated in this work, cosmological behavior of universe filled with inhomogeneous dark energy coupled or not with dark matter in f (T ) modified gravity. It has also In the graph at the left hand-side, we can see that at the beginning, the bulk viscosity increases (for low interact-ing parameter), but decreases (for high ξ ) and converges to the constant values as the time evolves indefinitely. Note that for ξ < 1, δ → 0. whereas in the graph at the right hand-side, δ is nul at the beginning but diverges as the time increases In these graphs, we can see that, the viscosity decreases at the beginning and becomes null as the time evolves but diverges for increasing time. Note that, the viscosity is not null for high interaction parameter ξ been studied, the bulk viscosity evolution describing future singularities throughout bounce cosmology. The first part of this work is devoted to analyse the dynamical evolution of EoS dark energy parameter that shows an universe in accelerated expansion. But with the strong interaction between the dark components, the universe evolves from the quintessence region to the phantomic ones for decreasing redshift. Note that in this case, the low interaction makes the universe only be in accelerated expansion phase as the time evolves. In comparison with [39], one note that with the lagragian f (T ) in power-law model, the results found are the same as the ones of [53]. In the second part of this paper, we investigate some bulk viscosity models solutions from gravitational equation that describe Little Rip and Pseudo Rip singularities in the framework of f (T ) gravity. This gravitational equation is a generalisation of that one found in [56]. We notice some corrections in the bulk viscosity and in the EoS parameter of dark energy due to the f (T ) modified gravity and the coupling parameter. The analysis of the Little Rip bulk viscosity solutions shows that it could vanish for low interacting parameter (ξ ) and converge for high value of ξ . In Pseudo Rip case, these solutions decrease at the beginning and increase at late time.

Data Availability Statement
This manuscript has associated data in a data repository. [Authors' comment: The data that support the findings of this study are available from the corresponding author, R. D. B., upon reasonable request].
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .