More on the holographic Ricci dark energy model: smoothing Rips through interaction effects?

The background cosmological dynamics of the late Universe is analysed on the framework of a dark energy model described by an holographic Ricci dark energy component. Several kind of interactions between the dark energy and the dark matter components are considered herein. We solve the background cosmological dynamics for the different choices of interactions with the aim to analyse not only the current evolution of the universe but also its asymptotic behaviour and, in particular, possible future singularities removal. We show that in most of the cases, the Big Rip singularity, a finger print of this model in absence of an interaction between the dark sectors, is substituted by a de Sitter or a Minkowski state. Most importantly, we found two new future bouncing solutions leading to two possible asymptotic behaviours, we named Little Bang and Little Sibling of the Big Bang. At a Little Bang, as the size of the universe shrinks to zero in an infinite cosmic time, the Hubble rate and its cosmic time derivative blow up. In addition, at a Little sibling of the Big Bang, as the size of the universe shrinks to zero in an infinite cosmic time, the Hubble rate blows up but its cosmic time derivative is finite. These two abrupt events can happen as well in the past.

Introduction

Several astrophysical observations (cf. for example supernovae type Ia [1, 2], cosmic microwave background [3], large scale structure [4], etc.) have confirmed that the universe is undergoing a state of accelerated expansion if homogeneity and isotropy are assumed on large scales. In addition, those experiments indicate that the matter content (i.e., the total mass-energy) of the universe, leading this accelerated expansion, must contain an exotic energy which is characterised with a sufficiently negative pressure. Dark energy (DE) is the most accepted hypothesis to explain the current observations, and constitutes roughly of the total matter content of the universe. However, so far there is no clear understanding of the true fundamental nature of DE. Indeed, the mysterious nature of dark energy is still among the long-standing problems in theoretical physics.

There are many dynamical models trying to explain the nature of DE [5, 6]. Among them, there is an attractive model which is inspired on the holographic principle rooted on quantum gravity [7-9]; it is the so called holographic dark energy [10, 11]. We next summarise the ideas behind this model. A well known fact is that the entropy of a given system with finite volume, , has an upper bound which is not proportional to its volume, but rather to its surface area, [12, 13]. In addition, for an effective quantum field theory with a given ultra-violet (UV) cutoff, , the entropy of that system scales as . Consequently, there is always a scale or a length where the quantum field theory with UV cutoff is expected to fail. This is expected to happen for large volumes or lengths. To solve this problem a link between UV and infrared (IR) cutoffs was proposed in [9]: By this mean the validity of the quantum field theory within this regime is assured. When the inequality is saturated, we can define an energy density which is inversely proportional to the square of the characteristic length of the system. These ideas have been applied to the universe giving rise to what is known as the holographic dark energy scenario [10]. The appealing holographic Ricci dark energy (HRDE) model consists in taking the square of the length characterising the Universe, L, as the inverse of the Ricci scalar curvature [14] (see also [15-21]).

It has been proven that the HRDE is suitable to describe the current acceleration of the universe as shown in [14, 22–25] . It has been equally shown that this model might induce a big rip (BR) singularity [14]; i.e., the scale factor, the Hubble parameter and its first cosmic time derivative reach very large values in a finite future cosmic time [26-28]. This model has been also constrained observationally [29]. More recent observational constraints on the HRDE can be found in Refs. [20, 21, 24, 30, 31] (cf. Ref. [32] for an extended list of references on the HRDE scenario).

On this paper, we intend to see how the BR1 present on the HRDE can be removed or appeased by the inclusion of interactions between cold dark matter (CDM) and the HRDE. An interaction on the dark sector and within the HRDE model has been previously analysed in [19–21, 30, 34–36] where the main goal of these papers was to study the adequacy (from an observational point of view) of these models to describe the late time acceleration rather than analysing the asymptotic behaviour of the universe. On this work, we will carry a thorough analytical analysis of the HRDE when a CDM and DE, given through the HRDE, are interacting. The goal of this work is to identify those interactions that are able to remove or smooth the BR. It is worth mentioning that an interaction on the dark sector has been favoured observationally [37, 38] and has been shown to be extremely helpful to mitigate the coincidence problem (cf. the recent review [39] on this topic).

The paper is organised as follows. In Sect. 2, we review the general setup of the HRDE model in a Friedmann-Lemaître-Robertson-Walker (FLRW) background in presence of an interaction term between DM and DE. In Sect.  3, we carry out a careful analysis of the asymptotic behaviour of the universe in this framework. Finally, in Sect. 4, we present our conclusions.

General setup

We consider a spatially flat FLRW universe, filled with matter with the energy density , whose evolution is described by the Friedmann equation:where is the reduced Planck mass. We assume that the total energy density, , of the cosmic fluid is described through a CDM component with the energy density and a HRDE component . The HRDE density is defined as [14]where and is a positive dimensionless parameter that measures the strength of the holographic component. From now on a zero subindex stands for quantities evaluated at present.

It is convenient to rewrite the Friedmann equation (2.1) in terms of the dimensionless energy densities:where . Therefore, Eq. (2.1) becomesThis equation constrains the cosmological parameters of the model at present time, , asin which the present value of the dimensionless Hubble rate E and its derivative with respect to x are governed by the following equations:In addition, at present time, the deceleration parameter reads,which must be negative (i.e., ) because the universe is accelerating currently. So that, the holographic parameter is constrained through the inequalityA characteristic of the HRDE model is that the energy density of any matter component is self-conserved. However, we will further assume in this model an interaction between the energy densities of CDM and the HRDE components. Therefore, the corresponding conservation equations read where the function Q denotes the interaction between the energy density of CDM and the holographic dark energy density components. Furthermore, positive Q represents energy transfer from CDM to DE, and vice versa for negative Q.

General equations in the presence of interaction

It is expected physically and also from Eqs. (2.9) and (2.10) that the interaction is defined through the energy densities involved in the system, in particular Q should be a function of the energy densities of CDM and HRDE components multiplied by a quantity with the unit of inverse of time. For convenience, we choose the Hubble rate H as the characteristic magnitude with units of inverse of time, and hence where is the critical energy density. Since the value of the interaction parameter Q is small (cf. Ref. [39]), a power law expansion of Q in terms of the energy densities of the system is doable. The first order terms of this interaction; i.e. linear interaction, corresponds towhere , and are constants.

Substituting Eq. (2.11) in Eq. (2.10), and replacing the energy densities in terms of dimensionless parameters, the conservation equation ( 2.10) can be written asOn the other hand, the Friedmann equation (2.4) can be rewritten as a differential equation of the dimensionless Hubble rate E with respect to x asNotice that the equations (2.12) and (2.13) form a system of coupled equations. However, by differentiating both sides of Eq. (2.13) with respect to x, and using Eq. ( 2.12), we obtain a second order differential equation for :The total conservation law for the system is given by where the total energy density is given by . The equation of state (EoS) for the HRDE isThis EoS can be rewritten by using Eqs. (2.3)-(2.4) asNotice that Eq. (2.16) is valid for any HRDE model where the total energy density of the universe is conserved.

In order to study the behaviour of the universe within the context of the HRDE, we will analyse the solutions of the differential equation (2.14) which depend on , , and parameters. We remind that we are interested in analysing the phase space of the parameters of the model for which the universe has a smooth future behaviour.

Expected behaviour: general solutions

The general solution of the equation (2.14) can be written as:with being constants,andThe parameter readsOn the other hand, for the case in which and , the solution for readswhere and are constants. The Hubble rate at the present time must satisfy the conditions of Eq. (2.6), i.e., the constants , and can be expressed in terms of , Eq. (2.7), as Furthermore, by substituting the solutions (2.17) and (2.21) in the second equation of (2.3), the dimensionless HRDE densities read, respectivelyandIn the rest of this paper, we will analyse the solutions of Eqs. (2.17) and (2.21), and we will discuss their asymptotic behaviours, which depend on the choices of the holographic parameter and the interaction constants and . We will show that, in the far future, the fate of the universe may end up in one of the following states: a BR singularity [26-28]; a Minkowskian or a de Sitter behaviour; a little sibling of the BR (LSBR) [40]; we will as well show for the first time the presence of two other possible asymptotic behaviours corresponding to what we named the little bang (LB) and the little sibling of the big bang (LSBB) (cf. Sects. 3.2 and 3.3).

Asymptotic behaviour and interaction effects

In this section, we will analyse the asymptotic behaviour of a FLRW filled with an HRDE fluid interacting with CDM within the model introduced on the previous section. In particular, we will analyse potential future singularities that might appear on the model by considering different interaction functions . First of all, we will start with a brief review of the standard HRDE model in the absence of interactions.

Standard HRDE model:

For vanishing parameters and , in Eq. (2.11), there is no interaction (i.e., ) and the standard HRDE model is recovered [14]. Then, by setting in solution (2.17) we obtain the dimensionless Hubble rate E(x) (for the case ) as (cf. Ref [14])The second term in Eq. (3.1) vanishes as , and the Hubble rate is governed only by the first term. This indicates that, in the far future, the HRDE energy density mimicking matter (second term in Eq. (3.1)) is practically zero, and the universe converges asymptotically to a universe filled with the dominant HRDE component. Therefore, the properties of the solution (3.1) depend on the different ranges of the holographic parameter . Let us summarise those behaviours as follows:Finally, for the case (which is not included in Eq. (3.1)), since , the universe tends to a Minkowski state in the far future.

For (), the HRDE density tends to zero and the universe heads to a Minkowski state in the far future.

For , the universe tends to a de Sitter state in the far future.

If , the HRDE is dominant at late time. Using Eq. (2.16), the equation of state for HRDE energy can be written as which is always smaller than . Therefore, the HRDE component behaves as phantom-like matter at late time. By integrating Eq. (3.1), we obtain the evolution of the scale factor a(t): where Notice that, we have set as the value of the scale factor at present. On the other hand, from Eq. (2.25), we find that the dimensionless energy density of the HRDE increases negatively with time (since ). Then, at a finite time, namely , the scale factor (Eq. (3.3)), the Hubble parameter and its cosmic time derivative blow up at . Therefore, the universe hits a BR singularity at . In fact, it can be seen that is a finite time and depends on the holographic parameter .

On the one hand, the latest observational data (e.g., Planck results [3]), implies that and , thus, (for a CDM universe). Then, by considering the condition (2.7), we can estimate the holographic parameter to be of the order . By setting Gyrs, we find that the BR would take place when Gyrs. On the other hand, we have seen that for there would be no abrupt events or singularities at late time. Given that we are interested in analysing DE singularities in this model, which are observationally favoured, from now on we will disregard values of such that .

The solution with

Let us now consider two particular classes of solutions for Eq. (2.14). For one of these solutions, the right hand side of Eq. (2.14) vanishes for any value of the scale factor. This corresponds to the solution (2.21) for , and . Consequently, we haveThen, the corresponding solution for E(x) in this case can be written aswhere A and B are constants.

By applying the condition (2.6), and using Eq. (2.7), we obtainSince we expect (this can be proven from a rough estimation based on the CDM model), from the right hand side of Eq. (3.7), we expect A to be always negative. Positiveness of in Eq. (3.6) implies that, x always lies in the range .

By using Eq. (3.6) in the relation and integrating both sides, we obtain, for negative values of A, a relation for the time dependence of x, aswhere denotes the present time for which . Notice that, we have set in the above relation. By taking the time derivative of x , and replacing it on the left hand side of equation , we obtain the dimensionless Hubble parameterConsequently, the time derivative of the Hubble rate (3.9) readswhich is constant during the evolution of the universe.

Since in this case, at wherex(t) reaches its upper limit , where the Hubble rate vanishes, and the first time derivative of the Hubble rate, Eq. (3.10), remains constant. Thus, at this point the universe hits a bounce. Hereafter, the universe starts to collapse during an infinite time, and as , for which , the Hubble rate (3.9) diverges while its time derivative (3.10) remains finite. This represents a new abrupt event for the future fate of the universe which happens at , and is different from the big bang singularity (for which the Hubble rate and its time derivative diverge). This abrupt event is smoother than a big bang: we named it the “little sibling of the big bang” (LSBB).

By using Eq. (2.1) and the total conservation equation , the total EoS readsIn the limit , the total EoS tends to . We would like to stress that the evolution of the universe is symmetric with regards to the bounce. Therefore, the universe evolves from a LSBB to a bounce and recollapse heading back to a LSBB.

The second class of solution is the case in which parameters and satisfy the conditions (3.5) but . In this case Eq. (2.14) reduces toWith respect to the cosmic time t, Eq. (3.12) becomeswhere, . By imposing the conditions (2.6) and (2.7), the dimensionless Hubble rate readswhere, from Eq. (2.6) at the present time , we getMoreover, the scale factor is given by2 Here, is a constant of integration. Notice that, as the energy density (2.2) of the universe diverges.

Likewise, the Hubble rate blows up in this case.

In order to better understand the behaviour of the Hubble rate (3.12), it is convenient to write its solution in the following form:where and are defined as For , the only physically interesting situations takes place when x belongs to the range : this describes a universe that starts its evolution at an infinite time in the past (where and ), at which the Hubble rate, E(x), and its time derivative, , diverge. Therefore, the universe hits a new abrupt event in this case which is similar to the big bang singularity, but which occurs at an infinite cosmic time. We thus name this new abrupt event as the “little bang” (LB).After this point, the universe expands till and then it bounces back to a LB. The other situation, , corresponds to an expansion that starts at a bounce when and heads to a little rip, when and [43-50]. Finally, when , there is a unique Lorentzian solution that interpolate between two bounces located at .

Interacting HRDE model

In this section, and in order to illustrate our purpose, we consider an interacting HRDE model with only the arbitrary parameters and . By introducing a new quantity r, the ratio between the energy densities of CDM and HRDE [41, 42]:and using the conservation Eqs. (2.9) and (2.10), we getBy substituting the interaction in Eq. (3.21), we obtainThe solutions for the differential equation (3.22) depend on the sign of the discriminant , Eq. (2.20), which can be rewritten aswhere , and are: To get the physical solutions of Eq. (3.22), it is helpful to distinguish three cases and :We will henceforth, study the late time behaviour of the universe predicted by the two physical solutions corresponding to and .

When , the parameter fulfils or (where ). In this case, the solution of the differential equation (3.22) reads, where , and is a constant of integration. As , the ratio r(x) in Eq. (3.27) tends to a constant: A physical solution implies that r(x) must be always a positive valued function, that is, ; this corresponds to . Consequently, this relation imposes a constraint on the holographic parameter . Since is always larger than , the positivity of the function implies that must be on the range for the solution (3.27) to be meaningful. In the far future (as ), by setting and for the range of the parameters satisfying , the ratio of the energy densities r(x) vanishes () and the budget content of the universe becomes dominated by the HRDE, as would be expected. On the other hand, if , the ratio r(x) becomes negative which is not physically possible. While by setting , in the far future, r(x) approaches the value Since r(x) is always positive, the parameter is constrained to fulfil .

For the case , the possible values of are or . In this case, the solution of Eq. (3.22) for r(x) is given by where is a constant of integration. On the limit , the solution (3.31) tends to Again r(x) must be positive, then, the holographic parameter must be on the range . Since is larger than , the only possible range for the holographic parameter in this case is

Finally, if then . The solution of Eq. (3.22) in this case is given as where is an integration constant. Equation (3.34) shows that when the argument of the ‘tangent’ term reaches the values where , the ratio r diverges with positive and negative signs, respectively. The former corresponds to a universe filled with CDM at late time, which does not match with the cosmological observations. In addition the latter limiting case corresponds to a negative ratio r(x) for the energy densities of the universe. Therefore, the solution given by Eq. (3.34) (that is, the solution provided by the case ) is not physically relevant.

The case

The dimensionless Hubble rate E(x) can be obtained by using Eq. (2.17) in which and are given bywith given by Eq. (3.23). Moreover, readsWhen , the Hubble rate (2.17) at late time becomesFor the two ranges of the parameter in which or , the argument of the exponential term in Eq. (3.37) is positive. In this case, the universe will undergo a BR singularity in the far future. By integrating Eq. (3.37), we can find the behaviour of the scale factor a(t) at late times asThis indicates that the BR occurs at whereand the scale factor diverges as . On the other hand, for the range of the holographic parameter such that , the argument of the exponential term is negative, the Hubble rate vanishes and the universe converges to a Minkowski state in the far future. In addition, for and the Hubble rate (3.37) reduces to a constant , indicating that the universe approaches a de Sitter state at late time.

By setting , Eq. (2.17) reduces toNow, we summarise the properties of the solution (3.40), for the case :The late time behaviour of the Hubble rate E(x) for the case , is the same as the general one given by Eq. (3.37); see Table 1 for more details.

Table 1

Summary of the behaviours of the universe at late times, for the physical range of holographic parameters , for different DM and DE interactions

Sec.	Interacting model	Δ	β	λ	σ_0,σ_±	Late time behaviour
3.1	λ_m=λ_H=0	Δ>0	β=12	–	σ_-=-3,  σ_+=0	de Sitter
		Δ>0	β<12	–	σ_-=-3,  σ_+>0	BR
3.2	λ_m=8β-1,	Δ=0	β<12	λ_c=0	σ_0=0	LSBB
	λ_H=2β(2β-1)^2			λ_c≠0		LR
3.3	λ_m≠0, λ_H=0	Δ=0	β<12	–	σ_0>0	Minkowski
		Δ>0	β<21+λ_m≤12	λ_m≥3	σ_±>0,	Minkowski
		Δ>0	β<12<21+λ_m	λ_m<1-2q_0<3	σ_-<0,  σ_+>0	BR
		Δ>0	β<12<21+λ_m	1-2q_0<λ_m<3	σ_-<0,  σ_+>0	Minkowski
3.3	λ_m=0, λ_H≠0	Δ=0	β<21+λ_H	λ_H<9	σ_0<-2(1+q_0)<0	Minkowski
		Δ=0	β<21+λ_H	λ_H<9	-2(1+q_0)<σ_0<0	LB
		Δ=0	β<21+λ_H	λ_H>9	-2(1+q_0)<0<σ_0	Minkowski
		Δ=0	β=18	λ_H=9	σ_0=0	LSBB^a
		Δ>0	β<27+λ_H	λ_H>3	σ_+>0,	BR
					σ_-<0
		Δ>0	27+λ_H<β<32(3+λ_H)	λ_H>0	σ_+>0,	BR
					σ_-<0
		Δ>0	32(3+λ_H)<β<12	all λ_H	σ_±<0	Minkowski
		Δ>0	β=32(3+λ_H)<12	all λ_H	σ_+=0,	de Sitter
					σ_-<0
3.3	λ_m≠0, λ_H≠0	Δ=0	β=β_2<27+λ_H-λ_m	–	σ_0>0	Minkowski
		Δ=0	27+λ_H-λ_m<β=β_2<21+λ_m+λ_H	–	σ_0<-2(1+q_0)<0	Minkowski
		Δ=0	27+λ_H-λ_m<β=β_2<21+λ_m+λ_H	–	-2(1+q_0)<σ_0<0	LB
		Δ>0	β<27+λ_H-λ_m or	–	σ_+>0	BR
			27+λ_H-λ_m<β<3-λ_m6+2(λ_H-λ_m)		σ_-<0
		Δ>0	β>3-λ_m6+2(λ_H-λ_m)	–	σ_±<0	Minkowski
		Δ>0	β=3-λ_m6+2(λ_H-λ_m)	–	σ_+=0, σ_-<0	de Sitter

This is a specific case of the solution of Eq. (3.5) given in the Sect. 3.2 (cf. see the third row in the table above)

For the case , the holographic parameter is within the range  . In this range, the first term in Eq. (3.40) becomes negative; , and the second one becomes positive; . Thus, the Hubble rate decreases and vanishes at which means that the universe hits a bounce at In the far future, after the bounce (), the dimensionless Hubble rate vanishes and the universe tends to a Minkowski state.

For the case , the holographic parameter satisfies the condition and we have always . Two possibles late time behaviour of the universe can be found. The case , i.e. , where the universe hits a BR singularity in the far future. The case , i.e. , where the dimensionless Hubble rate decreases and vanishes at (given in Eq. (3.41)) and the universe bounces at this point. In the far future, after the bounce (), the Hubble rate vanishes and the universe tends to a Minkowski state.

In this case, where or , the Hubble rate is given by Eq. (2.21):where and are defined in Eq. (2.24) as , where is the argument of the exponential term in the equation above. The total EoS of the universe in this case readswhich converges to as .

Following the discussion in item 3.3 above, we consider only the case . For , the argument of the exponential term is always positive (), while . This indicates that the universe undergoes a bounce at some in the future. Hereafter, the universe starts to collapse and as it converges into a Minkowski state.

For (cf. Eq. (3.33)), the argument of the exponential term is always negative ():

if , then ; the exponential term decays faster than the linear term , thus, the universe tends to a Minkowski state in the far future.

However, when (so that ), the universe will bounce at in the future, afterwards, it will start collapsing and in an infinite time, where (i.e., ), the Hubble rate and its time derivative diverge. Therefore, the universe meets a LB at late time. Notice that, in this case (), the matter content of the universe at late time behaves as phantom matter.

In the case of interaction parameter , if and , the term in the bracket in Eq. (3.42) is negative, while the argument of the exponential on the same equation is positive (for ). Therefore, the universe bounces in the future at some . After the bounce, the universe starts collapsing towards . Hereafter, the term in the bracket will be always positive, while the argument of the exponential term will evolve negatively. Whence, in the far future (), as the Hubble rate tends to zero and the universe tends to a flat Minkowskian state.

Finally, for the choice of , we get , and . In this case, Eq. (3.42 ) reduces to Eq. (3.6), we haveThis implies a bouncing scenario for the universe at . Hereafter, the universe starts to collapse during an infinite time, and as x tends to the scale factor tends to zero, the Hubble rate diverges while its time derivative remains finite; . Therefore, in the far future (), the universe hits a LSBB abrupt event. In summary, what happens is that the universe evolves from a LSBB in the past, expands until it bounces and heads back to a LSBB.

Whenever , the dimensionless parameter satisfies . Consequently, by substituting in the Eq. (2.16) we find that the EoS lies in the rangeThis relation implies that, the universe starts from a LSBB in the past with EoS for the HRDE , and reaches a bounce when at some in the future. After the bounce, the universe will collapse and will hit a LSBB abrupt event in the far future () with the same EoS, , of the initial state. Since the EoS of the dark energy evolves in the region , the corresponding interacting HRDE model represents a phantom-like behaviour, although its EoS is too negative to describe nowadays universe. Notice that, when the universe hits a LSBB in the far future, the total EoS, , of the universe tends to which has the same value for the total EoS at the initial LSBB abrupt event.

If , for the particular case , the interaction parameters are equal and the condition (3.5) is satisfied for . Then, the Hubble rate is obtained from Eq. (3.6) asOn the case where the holographic dimensionless parameter satisfies , the universe undergoes a bounce followed by a LB in the far future. It can be shown that at the LB , whereas at the bounce it reads . Therefore, and the corresponding interacting HRDE model has a phantom-like behaviour. Moreover, when the universe hits a LB in the far future, the EoS of HRDE tends to again. In fact, there is a symmetric evolution with respect to the source; i.e. the universe heads from a LB to a bounce and then back to a LB.

All the solutions we have analysed in Sect. 3 are summarised in Table 1.

Summary and conclusions

In this paper, we have considered a HRDE [14], as the dark energy component of the universe, coupled to the CDM component of the universe. As it is well known, in the absence of interaction () and depending on the physically relevant range of the holographic parameter , the late time universe will end up in a big rip (BR) singularity if [14]. We remind as well that those values of are consistent with the latest observations.

On this work, we are interested in investigating whether or not different interactions between HRDE and CDM components of the universe could resolve or smoothen the BR singularity. We considered different interaction functions such as , and to study the late time behaviour of the universe (r is defined in Eq. (3.20)).

In the presence of the general function , the corresponding differential equation (2.14) governing the dynamical evolution of the universe, possesses two types of generic solutions (cf. Eqs. (2.17) and (2.21)) depending on the discriminant function , given by Eq. (2.20).

For different choices of the interaction constants , , , and the HRDE parameter , the physically relevant solutions corresponding to (cf. table 1) lead to an asymptotic behaviour of the universe with: (i) a BR singularity, (ii) a Minkowski or (iii) a de Sitter state in the far future (see table 1 for a summary of the late time behaviour of the different solutions).

In particular, we have shown that for proper combinations of , and , the BR singularity can be removed.

On the other hand, when (and ), the general Hubble rate is given by Eq. (2.21). When and (or and ), the universe approaches a Minkowski state, asymptotically (cf. Table 1).

However, the case with where , represents two new abrupt events in the far future of the universe. We next summarise carefully these two cases:Finally, a similar analysis for the case and , for a non-vanishing interaction parameters , represents a universe which starts its evolution from an infinite past with a LB, and afterwards it expands until it bounces at a finite time in the future, and then it recollapses again to a LB; or it can start from a bounce and heads to a little rip at an infinite cosmic time. For , there is a unique Lorentzian solution that interpolate between two bounces.

For , the Hubble rate reads , while . This solution bounces at . After this point, the universe starts to collapse and as , the Hubble rate diverges while its time derivative remains finite. This corresponds to a new abrupt event, which is smoother than the big bang singularity and happens in an infinite cosmic time. We have called it the “little sibling of the big bang” (LSBB).

For , the Hubble rate reads , therefore . This solution indicates that, the universe bounces first at , then will collapse, and as , the Hubble rate and its time derivative diverge at an infinite time. Therefore, in the far future (), the universe tends to another type of abrupt event which we have called the “little bang” (LB) (cf. table 1).

These new types of abrupt events arise only when an interaction between HRDE and CDM is present. We have further shown that, under these conditions the HRDE may have a phantom like behaviour.