An Effective Sign Switching Dark Energy: Lotka-Volterra Model of Two Interacting Fluids

One of the recent attempts to address the Hubble and $S_8$ tensions is to consider the Universe started out not as a de Sitter-like spacetime, but rather anti-de Sitter-like. That is, the Universe underwent an"AdS-to-dS"transition at some point. We study the possibility that there are two dark energy fluids, one of which gave rise to the anti-de Sitter-like early Universe. The interaction is modeled by the Lotka-Volterra equations, commonly used in population biology. We consider"competition"models that are further classified as"unfair competition"and"fair competition". The former involves a quintessence in competition with a phantom, and the second involves two phantom fluids. Surprisingly, even in the latter scenario it is possible for the overall dark energy to cross the phantom divide. The latter model also allows a constant $w$"AdS-to-dS"transition, thus evading the theorem that such a dark energy must possess a singular equation of state. We also consider a"conversion"model in which a phantom fluid still manages to achieve"AdS-to-dS"transition even if it is being converted into a negative energy density quintessence. In these models, the energy density of the late time effective dark energy is related to the coefficient of the quadratic self-interaction term of the fluids, which is analogous to the resource capacity in population biology.

One of the recent attempts to address the Hubble and S8 tensions is to consider the Universe started out not as a de Sitter-like spacetime, but rather anti-de Sitter-like.That is, the Universe underwent an "AdS-to-dS" transition at some point.We study the possibility that there are two dark energy fluids, one of which gave rise to the anti-de Sitter-like early Universe.The interaction is modeled by the Lotka-Volterra equations, commonly used in population biology.We consider "competition" models that are further classified as "unfair competition" and "fair competition".The former involves a quintessence in competition with a phantom, and the second involves two phantom fluids.Surprisingly, even in the latter scenario it is possible for the overall dark energy to cross the phantom divide.The latter model also allows a constant w "AdS-to-dS" transition, thus evading the theorem that such a dark energy must possess a singular equation of state.We also consider a "conversion" model in which a phantom fluid still manages to achieve "AdS-to-dS" transition even if it is being converted into a negative energy density quintessence.In these models, the energy density of the late time effective dark energy is related to the coefficient of the quadratic self-interaction term of the fluids, which is analogous to the resource capacity in population biology.

I. INTRODUCTION: COSMOLOGY WITH SIGN SWITCHING DARK ENERGY
The Hubble tension [1][2][3][4] and the S 8 tension [5][6][7] in cosmology continue to be highly debated [8][9][10][11][12][13][14][15][16][17].The former is the mismatch between the locally measured expansion rate and the inferred rate via the cosmic microwave background (CMB), while the latter concerns the measurement of galaxy clusters on a scale of 8 h −1 Mpc, which revealed that matter has not clumped as much as expected assuming the concordance ΛCDM cosmology and its parameters constraints given by the CMB data.(For a review on these issues, as well as other challenges facing ΛCDM cosmology, see [18].See [19] for an introduction to various dark energy scenarios).If these effects are real, they could be due to modified gravity or other new physics [20,21].
Note that in the ΛCDM model, the Hubble parameter as a function of the redshift H(z), is specified by two constant fitting parameters [22]: (H 0 , Ω m ) or equivalently (A, B) as follows: = A + B(1 + z) 3 , where A is the term associated with dark energy.The aforementioned tensions could mean that ΛCDM is not correct and thus the "constant" fitting parameters could evolve with redshift (or equivalently, with time).In [22] it was noted that using constraints from H(z), observed data exhibit an increasing Ω m (decreasing H 0 ) trend with * Electronic address: ycong@yzu.edu.cnincreasing bin redshift, and yields a 'pile up' around A = 0 or Ω = 1.If Ω can increase beyond unity, this amounts to the dark energy density switching sign.(See also [23,24].)Indeed, one of the possible ways to ameliorate these tensions is to consider the possibility that the Universe was originally more anti-de Sitter-(AdS)-like than de Sitter-(dS)-like [25][26][27][28][29][30][31][32][33][34][35][36][37][38][39][40][41].That is, the physics of dark energy (DE) might be more complicated than we initially expected 1 .The main idea is to reduce the tension between the higher value of H 0 obtained from CMB with the lower value obtained by local measurements, by changing cosmological models either at the recombination epoch or at late time [11].For example, a phantom energy at late time would accelerate cosmic expansion faster than a cosmological constant would.In addition, Baryon Oscillation Spectroscopic Survey (BOSS) of SDSS-III probing the Lyα forest of quasars also indicates preference for a positive dark energy density at late time but a negative one at early time [48].Future observations such as SKA [49], BINGO [50,51], and Euclid [52] could potentially further constrain this possibility.Another observational support for negative energy density comes from Pantheon+ data of high redshift supernovae [53].
Such a possibility can be realized by simply promoting the cosmological constant Λ to 2 where Λ s0 is the present value of the cosmological constant, and z † is the value of redshift at which the sign switching abruptly happened [34].See also [54].Another model considers a "graduated dark energy" with energy density and pressure satisfying ρ + p ∝ ρ λ , which provides a continuous transition (controlled by the parameter λ) from AdS-like to dS-like Universe [28].Other options include the possibility that the dark energy sector could consist of a negative cosmological constant and a phantom dark energy [27] or a quintessence [35].(See, however, [55].)A different approach based on fractal modification to the entropy via a running Barrow index [56] could also give rise to an effective sign-changing dark energy [38].
In this work, we shall consider what happens if instead of a scalar field on top of a negative cosmological constant, we have either (1) a quintessence with negative energy density, which competes with a phantom dark energy with a positive energy density 3 , or (2) a phantom with a negative energy density that competes with another phantom with a positive energy density, or (3) a phantom with positive energy density being converted into a negative energy density quintessence.We shall refer to these scenarios, respectively, as "unfair competition", "fair competition", and "conversion" models.This is inspired by the interacting models between dark matter (DM) and dark energy [59][60][61][62][63][64], as well as from biological species interactions.In fact, the connection between these two subjects has been noticed in the literature.For example, in [65] Perez et al., as well as Aydiner in [66], pointed out that the DM-DE interaction can be re-written as the Lotka-Volterra equation, which is commonly used in population biology to model the interactions between various species.In cosmological contexts, Lotka-Volterra equation was also studied in [67,68].In the population model, it is of course required that the numbers of the species involved are non-negative, whereas in our model, the corresponding quantities are the energy densities of the DE fluids, which can be negative by assumption.It should furthermore be mentioned that if one considers non-minimal interactions between DE and DM to reduce the Hubble tension, the S 8 tension would in turn be exacerbated, hence we consider the alternative of non-minimal interaction between two DE fields to get one effective DE that exhibits sign-switching energy density, which in principle can address both tensions simultaneously.Our work is only meant to be an illustration of concept with the simplest models.Some assumptions would need to be relaxed or improved before a more realistic model can be used for data fitting the actual universe.
We will work out the conditions on the interactions between two dark energy fluids ("DE-DE interaction") in order to obtain a late time accelerated expansion with a very small but positive energy density.In our models, unlike the single fluid models in the literature, neither fluid exhibits any singular behavior in their equation 3 The idea that different dark energy components might interact with each other is not new.For example, models in which a quintessence interacts with a Chaplygin gas was considered in [57], in an attempt to explain the coincidence problem [58].
of state, although if the phantom divide is crossed, the combined effective dark energy typically does exhibit such a singular behavior during the AdS-to-dS transition.Remarkably, we found that in the fair competition model, it is possible for the effective dark energy to cross the phantom divide despite both component fluids satisfy w < −1.In addition, in this scenario there are evolutions that allow AdS-to-dS transition without crossing the phantom divide, which therefore is free of singular behavior in its equation of state.Finally, it is often said that the fact that the dark energy density is extremely small is "unnatural".We shall see that in the two-fluid model, this value is related to the coefficient of the quadratic self-interaction term of the fluids, which mathematically plays the same role as the resource capacity of a biological population.

II. THE UNFAIR COMPETITION MODEL
In [65] and [66], the authors considered a model in which dark energy is being converted into dark matter via where γ > 0. Likewise, dark matter can be converted into dark energy with γ < 0. See [69] for generalizations.
The equations of state considered in [66] are w DM 0 and w DE < −1, i.e., dark matter is a "normal matter" while dark energy is a phantom fluid.Thus, one can define two positive quantities: Furthermore, assuming that the Hubble parameter is slowly varying (so that R 1 and R 2 are both approximately constant) and upon introducing4 x 1 := γρ DE /r 2 and x 2 := γρ DM /r 1 , Eq.( 2) can be re-written as which is explicitly a Lotka-Volterra equation that describes a predator-prey dynamic with x 1 being the "prey" and x 2 the "predator".The system thus displays an oscillatory behavior.Such an interacting model could potentially help to resolve the coincidence problem.
In our case, we have two interacting dark energy fluids, which we will denote by Λ 1 < 0 and Λ 2 > 0 (despite the notation, they are not constant; the notation is meant to remind us that they are mimicking cosmological constants).Their energy densities would be denoted by ρ Λ1 and ρ Λ2 , respectively.The transition from an early time AdS-like Universe to a late time dS-like spacetime thus amounts to Λ 1 becoming subdominant to Λ 2 .Let us first consider an unrealistic model (to be improved upon later) that is analagous to Eq.( 2), so we can point out the differences: where γ > 0. As before, we assume that the late time dark energy is a phantom, thus w Λ2 < −1.On the other hand 5 , w Λ1 > −1 but with ρ Λ1 < 0. This is justified by [71,72], in which it was argued that solving both the Hubble and S 8 tensions require the overall effective dark energy to cross the phantom divide (if we assume that Newton's gravitational constant is not varying).See also [73].This is also similar to the model in [28].
Thus we define, analogous to Eq.( 3), two parameters We assume that both w Λ1 and w Λ2 are constant.Upon introducing the dimensionless quantities x := γρ Λ1 /r 2 and y := γρ Λ2 /r 1 , we obtain the Lotka-Volterra equation of the form whose some signs are different from that of interacting DM-DE model in Eq.( 4), and with x < 0, y > 0. This system does not oscillate, but rather there is an attractor x → 0 − and y → ∞.Note that in contrast to the DM-DE interaction model, it is not quite right to say that Λ 1 is being "converted" into Λ 2 here (a true conversion model will be studied in Sec.(IV)), since x < 0 implies that the interaction term is negative for both x and y.In other words, the two fluids are competing, but unlike two competing species whose birth rates are both positive, x is itself diminishing exponentially due to the "death rate" term −r 1 x.Hence the name "unfair competition".It is clear that the phantom fluid thus dominates over the quintessence.That is, ρ Λ1 is asymptotically vanishing while ρ Λ2 becomes large (and eventually diverges) at late time.This is not a desired property since we know from observation that dark energy density is only of the order of 10 −30 g/cm 3 .We can improve upon this model by modifying the dy/ dt term so that where K > 0 is a constant.The attractor is then x → 0 − and y → K, so we can prescribe to K the observed value.At this point this seems rather ad hoc, but later on we will give it a physical interpretation.
FIG. 1: The phase diagram of the two dark energy fluid model in Eq.( 8).
The analogy to population biology can also be made here: with x small at late time, in the absence of y/K term, what we have is analogous to an exponential growth population dy/dt = r 2 y, whereas Eq.( 8) corresponds to the logistic model with a resource capacity (or "carrying capacity") K, so that the actual population size cannot diverge but rather asymptotes to a constant value.Note that since x → 0 − anyway, there is no good reason to add the resource term to dx/ dt in this model.A typical phase diagram is given in Fig. (1).
Note that the effective dark energy density is the sum ρ Λ1 + ρ Λ2 .Thus, in order that this quantity starts out negative, we need the initial condition to satisfy y(0) < |x(0)|.Then, since the attractor is (x, y) = (0, K) with K > 0, it follows that by continuity x(t) + y(t) must cross-over to y > 0 at some point.The exact profile of ρ Λ1 (t) + ρ Λ2 (t) or the re-scaled equation x(t) + y(t) would depend on the initial conditions, but the transition from FIG. 2: Top: The evolution of x(t) (bottom curve in blue) and y(t) (top curve in red) with the initial condition set to be x(0) = −3, y(0) = 0.1, and with r1 = 1 = r2, K = 2. Bottom: The evolution of x(t) + y(t), which is essentially the re-scaled version of ρΛ 1 (t) + ρΛ 2 (t) (in this example, they are in fact equal).
an overall negative energy density to a positive one can be smoother than the model in [28].One example is given in Fig. (2).
If we know what kind of profile is desired from observational constraints, this would in turn provides us a mean to choose the coefficients r 1 and r 2 , as well as the initial conditions of the Lotka-Volterra equation.We also note that the equation of states of both dark energy components are never singular, though the combined effective dark energy density has to pass through ρ Λ1 (t) + ρ Λ2 (t) = 0, and the effective equation of state is singular at that point [74].This can be seen as follows: if we consider the combined fluid to still satisfy the equation of state of the form p = wρ, where ρ := ρ Λ1 + ρ Λ2 and p := p Λ1 + p Λ2 , then which implies that the effective varying w is the weighted average: Therefore, evidently w → ±∞ when the denominator is zero.Again, the situation is similar to the "graduated dark energy" model of [28], though in that case there is only one dark energy fluid, whose equation of state becomes singular.An example of the evolution of w is shown in Fig. (3).Note that the values of w Λ1 and w Λ1 are not freely prescribed since they are constrained by Eq.( 6), in the sense that once we fix r 1 and r 2 , the relations between w Λ1 and w Λ2 are also determined.In our example, choosing We remark that the sign of x(t) + y(t) is not necessarily the same as the sign of ρ Λ1 (t) + ρ Λ2 (t).In fact, Thus the sign of x(t) + y(t) is the same as the sign of r 1 ρ Λ1 (t) + r 2 ρ Λ2 (t).For simplicity our example deals with r 1 = r 2 , and the two expressions do have the same sign, and furthermore w is singular when x + y = 0 or equivalently at r 1 ρ Λ1 + r 2 ρ Λ2 = 0.

III. THE FAIR COMPETITION MODEL
In population biology, two species that are very similar (i.e., fulfilling the same ecological niche) will compete for the same resources.To model such a situation, we consider both species to have a positive "birth rate", so that in this sense the competition is fair.In place of Eq.( 8), we have: where now r 1 := −3H(1 + w Λ2 ) > 0. Note the coefficient pre-multiplying r 1 is now 1 instead of −1.Since x < 0, we consider the resource term to be 1 + x/K 1 instead of 1 − x/K 1 , keeping K 1 > 0. Note again that "competition" means that the interaction is harmful for both species, so the interaction term is negative for both fluids (r 2 xy < 0 because x < 0).For a fair competition we also include the carrying capacities K 1 and K 2 for both fluids, with O(K 1 ) = O(K 2 ) and both being positive.A typical phase diagram is shown in Fig. (4).Two families of flows are observed: those that flow towards K 1 and those that flow towards K 2 .We are concerned with the latter.
Inserting a few terms that cancel each other yields: FIG. 5: Top: The evolution of x(t) (bottom curve in blue) and y(t) (top curve in red) with the initial condition set to be x(0) = −1, y(0) = 0.9, and with r1 = 1, r2 = 1.5;K1 = 2, K2 = 3. Bottom: The evolution of x(t) + y(t), which is essentially the re-scaled version of ρΛ 1 (t) + ρΛ 2 (t). Therefore, If r 2 > r 1 , we can write the last equation as Likewise, if r 1 > r 2 , we can write Thus, for example, if r 2 > r 1 (and hence r 2 2 > r 2 1 -recall that r 1 , r 2 are positive), and we observe that x + y and x are both monotonicaly increasing, then so must ρ Λ1 + ρ Λ2 : The phantom divide can be crossed as shown in Fig. (6), where w is given by Eq. (10).The values of w Λ1 and w Λ2 are constrained by the defining equations r 1 = −3H(1 + w Λ1 ) and r 2 = −3H(1 + w Λ2 ).With r 1 = 1 and r 2 = 1.5, if we take w Λ1 = −1.01,say, then w Λ1 = −1.015.Note that there are two phantom crossings 6 here: the first one FIG.6: The evolution of the overall effective w of the example in Fig. (5) with wΛ 1 = −1.01 and wΛ 1 = −2.01.There are two phantom crossings.The first crossing is smooth and corresponds to the time when x + y = 0.The second crossing is singular and corresponds to ρΛ 1 + ρΛ 2 = 0.These two conditions are not the same since r1 = r2.
occurs without any singularity.From Eq.( 10) it can shown that if the denominator is not zero, then such a smooth crossing occurs precisely when x + y = 0.This cannot happen for the unfair competition model as the condition would be x − y = 0 instead (which cannot occur since x < 0 and y > 0).Note that despite the fact that there are two phantom crossings, there is only one AdS-to-dS transition in this example.
Even more surprising is the fact that the overall w < −1 can stay constant, yet there is still a AdS-to-dS transition.
To achieve this we simply need to choose r 1 = r 2 .The plots of x(t), y(t) and their sum is qualitatively the same as Fig. (5) and are thus not shown.However, from the defining relations of r 1 and r 2 we would have w 1 = w 2 , and so in Eq.( 10), we obtain Strictly speaking during the transition point ρ Λ1 +ρ Λ2 = 0, which otherwise would give rise to a singular behavior, we get an indeterminate form 0/0, but both the left and right limit is well-defined and equal to w 1 , so physically it makes sense to say that w ≡ w 1 = w 2 for all t.As such this evades the recent theorem that a sign-changing dark energy must have a singular equation of state [74].The reason this does not really violate the theorem therein is because the proof in [74] is strictly for DE fluids that obey the usual continuity equation ρ + 3H(ρ + p) = 0, whereas in our model it can be checked that the combined DE does not satisfy the continuity equation; the "carrying capacity" term K breaks the continuity equation.This is equivalent to saying that ∇ µ T µν DE = 0.This is not surprising -as we will see in the Discussion section, our models have nontrivial nonlinear self-interaction term that acts as a source term for the continuity equation.

IV. THE CONVERSION MODEL
Given the results above, one might wonder whether the AdS-to-dS transition can still happen if we restrict the growth of Λ 2 by it into Λ 1 , or equivalently, by giving Λ 1 an advantage.This is achieved by the following model involving a quintessence Λ 1 and a phantom Λ 2 : in which we note that the second term of dx/ dt is now +r 1 xy instead of −r 1 xy.The interaction is therefore beneficial to x but harmful to y.This sounds like the complete opposite of what we wish to achieve (to have Λ 2 being the dominant term at late time).Surprisingly even in such a scenario it is possible to have a phantom crossing.The only difference being the attractor is now a stable spiral centered at as can be seen in the example depicted in Fig. (7).What happens is that, despite the conversion term, Λ 2 can still dominate at late time.After all, we do not need Λ 2 to grow too big.The evolution of x(t) and y(t) as well as their sum (which in this example is the same as ρ Λ1 + ρ Λ2 ) are shown in Fig. (8).The phantom crossing is shown in Fig. (9).
The reason we consider a resource term 1−x/K 1 instead of 1 + x/K 1 as in the previous section (incidentally, this puts the fixed point x = K 1 > 0 outside the physical phase space) is that otherwise, with 1 + x/K 1 and the coefficient pre-multiplying r 1 being negative (quintessence) instead curve in red) with the initial condition set to be x(0) = −3, y(0) = 0.1, and with r1 = 1 = r2; K1 = 1, K2 = 2. Bottom: The evolution of x(t) + y(t), which in this case is the same as the overall ρ. of positive (phantom), we observe that for "most" initial conditions, and thus the magnitude of x is increasing and x → −∞ instead of 0, which is not the behavior that we want.This can be seen in Fig. (10).However, even in this scenario there is one interesting feature worth mentioning.In the neighborhood of the origin, there exists a center around which the flows are cyclic.This implies both x and y, as well as their sum, are oscillatory.As a result, there are multiple (infinitely many) phantom crossings, and infinitely many transitions between AdS-like to dS-like cosmology.These are shown in Fig. (11) and Fig. (12).Indeed, multiple transition scenarios have been considered in the literature [39].

V. DISCUSSION: SIGN SWITCHING DARK ENERGY AND NATURALNESS
One of the longstanding questions about dark energy density is why its value is so small, which is some 10 −120 times smaller than the natural scale for a quantum vacuum energy if it is indeed a cosmological constant (for a dynamical field, the problem translates into an extremely light mass of the field).Of course, it is debatable whether this is indeed a problem [75].In any case, it would be interesting to see what this value corresponds to in the Lotka-Volterra equations in these models.
Take for example, the unfair competition model.We note that the evolution equation for y(t), namely is equivalent to the following fluid equation: In other words, the "resource term" in the Lotka-Volterra equation corresponds to a quadratic self-interaction term.How might one interpret this term?Such a term was also considered in [66] and [76].As commented therein, pressure and density may not be FIG.11: Top: The evolution of x(t) (bottom curve in blue) and y(t) (top curve in red) with the initial condition set to be x(0) = −1, y(0) = 0.5, for the model that corresponds to Fig. (10).Bottom: The evolution of x(t) + y(t), which in this case is the same as the overall ρ. linearly related in more complicated and more realistic systems.If we assume p = p(ρ) for any barotropic fluid to be an analytic function, we can consider equation of state of the form p = p 0 + A 1 ρ + A 2 ρ 2 + O(ρ 3 ).This is a Taylor expansion of p = p(ρ) about ρ = 0, or upon regrouping of terms, the expansion about the present energy density [76,77].If this is the correct interpretation, then the self-interaction term in Eq.( 25) can be interpreted as the result of the first order non-linear term in the expansion.However, in a series expansion, typically the coefficients of the subsequent terms are roughly of the same order of magnitude, so the "natural" expectation is that γ K 1+wΛ 2 1+wΛ 1 ∼ O(γ).Even if w Λ1 and w Λ2 can be very close to −1, generically we would have the ratio (1 + w Λ2 )/(1 + w Λ1 ) to be of order 1.This means that K being small O(ε) 1 would typically result in the quadratic coefficient being unnaturally large and dominate over the linear term, which in turn suggests that we should not, in fact, interpret this term as a term in a Taylor series expansion of p Λ2 = p Λ2 (ρ Λ2 ).Note that if we do not interpret ρ Λ2 and ρ 2 Λ2 term as part of a Taylor series, we can still absorb the ρ 2 Λ2 term as part of the pressure so that p Λ2 = w 2 ρ Λ2 + const.ρ 2 Λ2 .In which case K is related to the mass scale M of Λ 2 via K ∼ M 4 , see [78,79]; so this is just the aforementioned fact that in the case of dark energy being dynamical, the naturalness problem is its small mass scale.In the conversion model, the situation is similar.The attractor of the spiral is given in Eq. (22).In which we see that y and x are both small if K 1 and K 2 are small.Obviously, our models do not solve the naturalness problem, unless one could explain dynamically why the attractor has such a small value.Perhaps a fundamental understanding of the nature of the phantom fluid or an entropic argument could provide such a mechanism.(In the context of cosmological constant, it was argued in [80] that gravitational entropy is maximized by Λ → 0 + .) To conclude, in this work, motivated by the idea that a sign switching dark energy from an early time AdS-like Universe to a late time dS-like Universe can help to ameliorate the Hubble tension and the S 8 tension, we consider a scenario in which the dark energy sector consists of two interacting fluids.We found that AdS-to-dS transition can happen under various models, even if both fluids are phantom, and even if the combined effective dark energy has a constant w.Of course, these are only toy models serve to illustrate the qualitative features.The profiles of these fluids need to be constrained by observations.Still, the possibility that the dark energy sector may contain various interacting components deserves a closer look (see [79,81] for other aspects of self-interacting dark energy) as it can realize many different interesting features.
For generalizations, one could also consider interactions between the two dark energy components with dark matter and/or dark radiation in a more complicated model (a quintom model was considered in [82], with the phantom component interacting with dark matter, but not with the quintessence sector; see also [83]).Another possibility is to consider other forms of interaction terms in place of ±γρ Λ1 ρ Λ2 ; see for example [84] and the references therein for the case of DM-DE interaction.Most importantly, more realistic models need to go beyond the assumption that the Hubble parameter is slowly varying when setting up the Lotka-Volterra equations.

FIG. 10 :
FIG. 10: The phase diagram of the two dark energy fluid model in Eq.(21) but with the resource term changed to 1 − x/K1 instead of 1 + x/K1.Here we use r1 = 1 = r2, K1 = 2 = K2.Most trajectories flow towards x → −∞ and y → 0. Note, however, the presence of a center surrounded by cyclic flows.