An exotic k-essence interpretation of interactive cosmological models

We define a generalization of scalar fields with non-canonical kinetic term which we call exotic k-essence or briefly, exotik. These fields are generated by the global description of cosmological models with two interactive fluids in the dark sector and under certain conditions, they correspond to usual k-essences. The formalism is applied to the cases of constant potential and of inverse square potential and also we develop the purely exotik version for the modified holographic Ricci type of dark energy (MHR), where the equations of state are not constant. With the kinetic function $F=1+mx$ and the inverse square potential we recover, through the interaction term, the identification between k-essences and quintessences of exponential potential, already known for Friedmann-Robertson-Walker and Bianchi type I geometries. Worked examples are shown that include the self-interacting MHR and also models with crossing of the phantom divide line (PDL).


Introduction
The cosmological acceleration suggested by astrophysical data [1] can be explained by very different models, among which, the simplest one is the Λ CDM [2]. A significant number of them reject the possibility of interaction between the modeler components of dark matter (DM) and dark energy (DE), and perhaps for that reason, they fail to justify the coincidence problem. On the other hand, there is no fundamental reason to assume an underlying symmetry which would suppress the coupling. Whereas interactions between DE and normal matter particles are heavily constrained by observations (e.g. in the solar system and gravitational experiments on Earth), this is not the case for DM particles. Moreover, the possibility of DE-DM interaction is a phenomenon consistent not only with recent Planck cosmological data [3] but also it looks like theoretically possible when coupled scalar fields are considered [4], with many different types of interactions [5]. With respect to DE, among the several candidates to play its role, such as vacuum polarization, vector models, tachyons, Chaplygin gas, k-essences, Cardassian expansion, quasi-steady state cosmology and scalar-tensor models [6][7][8][9][10][11][12][13][14], the k-essence cosmology has received a lot of effort [15][16][17] also playing both, the dark matter and dark energy roles [18]. The k-essence would explain the coincidence problem because the transition between the tracker behavior during the radiation-matter doma e-mail:forte.monica@gmail.com arXiv:1507.03658v3 [gr-qc] 27 Dec 2015 ination and a cosmological constant-like behavior, appears to arise for purely dynamical reasons without fine-tuning necessity.
In this work we try to connect both issues, cosmological systems of two interactive fluids on one side and unified models handled by a single field on the other side, with the intention to use the results of both fields of study, regarding them as different approaches to the same model. The idea of describing a universe filled with two interactive fluids with constant equations of state ω 1 and ω 2 , through a unified model driven by a single scalar field was first implemented in [19] and that produced a new type of field, dubbed exotic quintessence. This field, affected by an exponential potential, differed from the usual quintessence because of the inclusion of the parameters ω 1 and ω 2 in the expressions of its evolution equation, its density of energy and its pressure. For the values ω 1 = 1 and ω 2 = −1, that is, when fluids are assumed to be stiff matter and vacuum energy, the usual quintessence was recovered. A special case of this field, when cold dark matter ω 1 = 0 and vacuum energy ω 2 = −1 are considered, was used by Liddle et al. in [20]. In this paper, we implemented the same idea through a new scalar field φ with non-canonical kinetic term, coming from a Lagrangian of [17]. We nicknamed exotic k-essence or exotik to this new class of non-canonical scalar, because of the new equations of motion that they fulfill, are reduced to a common k-essence under certain conditions. The connection with the work in [19] holds in the special case in which the quintessence approach is reduced to a common quintessence with exponential potential and the k-essence approach is reduced to a common k-essence defined by the kinetic function F(x) = 1 + mx and is affected by a potential of inverse square. In that case both models are described by the same scale factor and so they are geometrically equivalent. [21], [22]. The mechanism of generation of mathematical expressions for the exotic field, favors the role of one of the two interacting fluids. However, either of the two fluids is also able to perform the mathematical description. The defining choice will depend on the particular interest of each model because equations will include terms such as (1 + ω 1 ) −1 if the key field is the second or (1 + ω 2 ) −1 if the key field is the first. In this paper we have in mind cosmological systems in which the fluid 1 is never considered as vacuum energy, and therefore, the logical choice is to take the fluid 2 as a key field generating factors (1 + ω 1 ) −1 that do not falsely divergent the physical magnitudes into consideration. The next step is giving to the key fluid (from now on fluid 2, playing the role of dark energy) a functional form that describes the specific description is intended to work. At this stage is where any anzats supported on specific physical models or described as combinations of physical quantities, appears. For example, the anzats taking linear combinations of density of energy and pressure is formally equivalent to the expression used for the dark energy fluid in modified holographic Ricci (MHR) type models (if [23]. The last step corresponds to express the above process in terms of known elements, that is, expressions of the energy density and the pressure of a common k essence, ρ = V (F − 2xF x ) and p = −V F. The various possible exotik descriptions (each described by a particular ansatz), are obtained by solving the equation of evolution for the key fluid affected by an interaction Q. Conversely, we can discover new types of interactions associated with known k-essences solving the same equation of evolution written in terms of known F and V.
Our paper is organized as follows. In Section II we consider the general interacting twofluid cosmological model and introduce the definition of exotic k-essence. In Section III we gain deeper insight into the subject analyzing the representation with constant potential and with square inverse potential and also giving first integrals in both cases. In Sections IV we establish the equation which connects the exotik with an arbitrary interaction. In sections V and VI we show worked examples for constant potential (including the modified holographic Ricci DE [24]), and for the square inverse potential (including the ratification of the equivalence between the linear k-essence F = 1 + mx and the quintessence with exponential potential) respectively.

Exotik
We consider a model consisting of two perfect fluids with an energy-momentum tensor T ik = T (1) ik + T (2) ik . Here T (n) ik = (ρ n + p n )u i u k + p n g ik , where ρ n and p n are the energy density and the equilibrium pressure of fluid n and u i is the four-velocity. Assuming that the two fluids interact between them in a spatially flat homogeneous and isotropic Friedmann-Robertson-Walker (FRW) cosmological model, the Einstein's equations reduce to: where H =ȧ/a and a stand for the Hubble expansion rate and the scale factor respectively and where we consider constant equations of state ω i = (p i /ρ i ) for i = 1, 2. The whole equation of conservation (2) shows the interaction between both fluid components allowing the mutual exchange of energy and momentum. Then, we assume an overall perfect fluid description with an effective equation of state, ω = p/ρ = −2Ḣ/3H 2 −1, where p = p 1 + p 2 and ρ = ρ 1 + ρ 2 . The dot means derivative with respect to the cosmological time and from Eqs. (1), (2) we get As it was done with an exotic canonical scalar field in [19], we propose that the interactive system as a whole be represented by an exotic field φ (labeled by the potential function V (φ ) and the kinetic function F(x), x = −φ 2 , through the relationship The global density of energy ρ and the global pressure p = ωρ are where ω 1 = −1, ∆ = ω 1 − ω 2 and the field φ satisfies the conservation equation for The equations (5)- (7) define the exotik field, a generalization of the k-essence field. So many classes of exotiks exist as functional forms we adopted for ρ 2 . The anzats with constants α and β , can be thought as the general linear combination of the density of energy and pressure of a common k-essence 1 On the other hand, a similar expression arises naturally in modified holographic Ricci (MHR) dark energy models [23,25] if we think that the global model is driving by a k-essence 2 This proposal allows us to understand the reason of the name Exotik for this k-field as follow. The application of (8) to the equations (5)-(7) transforms them into and In this context, our exotik representation of the interacting system results in a habitual k-essence when α = (1 + ω 1 )/∆ and β = ω 1 /∆ for arbitrary pairs of fluids. This motivates the name exotik or exotic kessence, for the general case where the two independent parameters α and β can not be considered superfluous or included into either the function F, or the potential V. Note that these very particular identifications for α and β imply that ω 1 and ω 2 are constant and so it is not a recommended option when we are using a MHR fluid as DE. 3 Because of in these holographic models always it is verified that ω 2 = (A−ω 1 −1)ρ 1 /ρ 2 +B−1 , a constant ratio of dark densities of energy arises, condition not necessarily true for general interactions. In the general statement and from the equations (9), the global equation of state ω is read as

Potentials and asymptotic behaviors
The evolution equation (9c) allows us to find the functional form of the exotik field once the potential and the kinetic function are given. The choice of the potential is subject to the type of comparison in which we are interested. For example, the constant potential is needed to contrast with purely k-essences and the inverse square potential is used to collate with quintessences with exponential potentials.
In the constant case, the equation (9c) has the first integral with m 0 being a constant of integration. In turn, this expression leads to see a novel feature of the equation of state (10) of these models, with the definitions C = α∆ /2m for graphical economy. It is known that the usual k-essences whose kinetic functions have a root [18], drive models with a dust behavior at epoch around t 0 where F(x(t 0 )) = 0. These unified dark A represents the maximum possible value of the overall equation of state (ω MHR ≤ A − 1) and therefore A ≤ 2. 3 Instead of (5)- (7) , the correct equations for the case of MHR fluid as DE, are energy models are included here, independently of the pair of fluids in interaction, when β = ω 1 /∆ . Actually, even if the kinetic functions have no roots, the same behavior is observed if a ν FF σ −1 x 1 when ν > 0. But our purely exotiks let go further and include the stages of radiative dominance at early times, properly choosing the representation with β ∆ = (3ω 1 − 1)/4. For ν > 0 all these models show an accelerated behavior (ω → −1) when a → ∞, providing appropriate behavior for our actual expanding universe.
For the inverse square case V = V 0 /φ 2 , equation (9c) can be rewritten in terms of the global equation of state ω asφ where u = α∆ /(1 + β ∆ ) and λ = (1 + ω 1 )/(1 + β ∆ ). This difficult equation has a first integral when we are dealing with common k-essence and more generally, asking u = 1, which is equivalent to having a single free parameter in (8). This first integral can be written in three different ways with m 1 a positive constant of integration. For the case u = 1 is λ = (1 + ω 1 )/(α∆ ) and the equation (9a) gives us the simple relation V F = (λ − 1 − ω)ρ which allows to see that ω = λ − 1 whenever F vanishes. Using (14b), the global EoS (10) is written as From (15) it follows that when F has a root, the global EoS is ω = λ − 1 whatever the sign of λ be. For times outside those cases, when F(x) is monotone and non-zero, and considering that in cases of interest is λ > 0, the asymptotic behaviors are: ω = λ − 1 at early times (a → 0) and 1 + ω = 2/(3F x ) at late times (a 1). The latter behavior generates a necessary condition to be satisfied by possible kinetic functions when the phantom regime is excluded. In that cases, from (10) and λ > 0, the only admissible kinetic functions are the monotonous increasing F.
The intention behind the attitude of not fixing α and β , is that the representation may still be used in the modified holographic case where these constants shape the upper limit of the value of the global barotropic index. Also, as we will show in the next section, some of them do not allow for adequate representation while others facilitate the resolution of the problem.

The interactions
The above results are quite general and apply to any kinetic function F(x), but the particular choice of the function will be determined by the interaction Q amending the evolution of both fluids. We define the interaction Q, through the partition of the global conservation equation (2), aṡ Then, the expressions (8), (9c) and (16b) let us write the equation that must be fulfilled by the kinetic function F(x) once the interaction Q(V, F) and the potential V are fixed.
with M = F x + xF xx and N = (α − 2β )F x − 2xF xx β . The expression Q(V, F) means that the interaction, often expressed as a function of ρ and its derivatives, should be given using equations (1), (2), (8) and (9). 4 The extremely complex equation (17) is made simple in several important cases as those where V = V 0 , constant.

Purely exotiks
The equation (17) for V = V 0 , is reduced to which is a highly nonlinear equation for F. However, the change of variables ζ = ρ us to obtain the more simple differential equation for ρ, with ρ = dρ/dζ and ρ = d 2 ρ/dζ 2 . This is the already known source equation for the energy density described in [27].
-Examples Q → F -Λ CDM The trivial case Q = 0, ω 1 = 0 and ω 2 = −1, which corresponds to the Λ CDM model can be represented by the kinetic function F(x) = F 0 + F 1 (−x) α 2β with F 0 and F 1 , two constant of integration. Using the first integral (−x) 1+β −α/2 F 1+β x = m 0 /a 3 we recover the known expression ρ = ρ 01 /a 3 + ρ 02 with ρ 01 = 2V 0 m 0 (−2β ) β (F 1 α) −β and ρ 02 = αV 0 F 0 . This case shows that in general, a purely k-essence will not be a correct interpretation of that interaction because it leads to a constant density of energy. Therefore, in principle, we must consider the two parameters α and β not fixed.

-CDM and MHRDE
This is a very interesting case because, although the dark energy EoS ω 2 is not a constant, the option for the density of dark energy ρ MHR relates ω 1 and ω 2 with A and B through the expression ω 1 ρ 1 + ω 2 ρ 2 = Aρ 1 + Bρ 2 . Taking ω 1 = 0 we define the modified interaction Q M througḣ and so, the equation (18) changes to The expression of ρ MHR 2 shows that there is a perfect agreement with the anzats (8) for α = A/(A − B) and β = (A − 1)/(A − B) and so the exotik function is a common k-essence. We apply this to the 4 The following expressions are useful to express the interactions null interaction Q = 0, which is equivalent to replacing Q M = Q + (1 − A)ρ 1 in (21). The solution is 1) and lets writing the densities of energy ρ MHR = b 1 a −3 + b 2 a −3B and ρ MHR 2 = ((A − 1)/(A − B))b 1 a −3 + b 2 a −3B . So, the purely exotik includes the modified holographic Ricci DE model with cold dark matter, where it can be seen that, even when Q is null, the dark energy component is far from remaining independent of the CDM. This is a consequence of the "holography" of the model introduced into the expression of ρ MHR 2 [25]. In other words, the MHR fluid, is always a self-interacting component.
Halfway between choosing none or both parameters, sometimes, set just one makes it easier to reach the goal. Select which and how to fix it, it is evident through the process of solving the equations. For example, in the case of interaction proportional to ρ , Q = τρ /∆ it is highly advisable to set β = −1/∆ . Then, we obtain the first integral (−x) −α∆ /2 = m 0 a −3(1+ω 1 ) and the appropriate exotik with The corresponding global density of energy ρ = ρ 01 a −3n+ + ρ 02 a −3n− and global EoS show that the global behavior is determined by the relationship between the strength of the coupling τ and the EoS of the considered fluids. As (n + − n − ) is always positive, the asymptotic values are ω → n + − 1 at early times and ω → n − − 1 at late times. The demeanor of the model is like a quintessence or phantom according to whether τ < ω 1 + ω 2 + 2 or τ > ω 1 + ω 2 + 2, but there is no crossing of the phantom divide line (PDL). The option τ = ω 1 + ω 2 + 2 is only admissible in the case of a cosmological constant (ω 2 = −1) and there, the effect of this interaction is to freeze the overall density of energy. Could this interaction be a mechanism to freeze the densities of elementary particles in the primordial times, assuming the existence of a cosmological constant? It is important to note that α and β determine the representation and we must work within it. However, the results of the corresponding models do not depend on these parameters. This should be clear in the next section where we go from F to Q, which is expressed through the magnitudes given by the particular representation. The effects are studied within the representation but not depend on it.
-Examples F → Q -An interesting example in the representation of common k-essence corresponds to the Chimento function F Ch (x) = 1 V 0 (2n−1) 2nn 0 √ −x − (−x) n , n = 0 and n = 1/2. Using (18) we obtain the associated interaction Q Ch Q Ch ∆ = (1 + ω 1 )(1 + ω 2 )ρ + (ω 1 + ω 2 + 1)ρ + 1 2n F Ch have a root at x r = −(2nn 0 ) 2/(2n−1) and an extreme at x ex = −n 2/(2n−1) 0 . The equation (12), now with σ = 1, lets us know that the cosmological model driven by Q Ch has an era with like dust behavior around a 3 root = m 1/(1+ω 1 ) 0 /nn 0 and a cosmological constant behavior at late times, regardless of the duet of fluid considered. The corresponding total density of energy is easily expressible from (9a) y (11) as Moreover, writing (12) as with v = (a/a root ) 3 it can be seen that these models present a dust like behavior at early times. Also, they look as a Λ CDM models with no interaction at all for n >> 1 in which case, the dust-like behavior at the root of F coincides with the dust-like behavior at the epoch a << 1.
-Another interesting example arises, again in the common k-essence representation, α = (1 + ω 1 )/∆ and β = ω 1 /∆ , when we are really dealing with a Chaplygin gas, through the kinetic function (18) we can express the interaction as and the corresponding density of energy ρ is obtained directly from (9a) and (11), as 6 Exotiks with inverse square potential For cases with inverse square potential V = V 0 /φ 2 the expression (17) can be written as -F(x) = 1 + mx As a first example we get the interaction term in the general representation (with α and β free) for the simplest kinetic function that supports a non-constant equation of state, F(x) = 1 + mx, with m constant. In this case (26) is written as The usual k-essence version (α = (1 + ω 1 )/∆ , β = ω 1 /∆ ) of (27) allows us to show, through the interaction term, that the models driven by k-essence with inverse square potential and kinetic function F(x) = 1 + mx, and those dominated by a quintessence ϕ with exponential potential U(ϕ) = U 0 exp((− 2/mV 0 ϕ) share the same geometry. This result was obtained in FRW [21] and also in Bianchi I cosmologies [22]. To test this feature, we must consider that the exotic quintessence used in [19] becomes a usual one when ω 1 = 1 and ω 2 = −1 and so, it is the scalar representation of two fluids with the interaction Q = −ρ 2 2(1 + ω)/(3mV 0 ). Replacing the corresponding values of the equations of state at the usual k-essence version of (27) we find the same equivalence because the interaction terms in both models agree. Moreover, note that the (27) vanishes identically the source equation [27] so that it is not possible to solve the system using that method. But here, using F = 1 + mx in (14c) we obtain the equation mV 0Ḣ + H 2 + 2m 1 Ha −3 + m 2 1 a −6 = 0 whose solutions for the factor of scale and the k-essence field have already been given in [21].
formerly used for constant potential and now for inverse square potential, shows the power of this interpretation in usual k-essence representation to resolve the evolution of systems with two arbitrary fluids controlled by a nonlinear interaction. From (26) this function leads to the highly nonlinear interaction with ν = (3V 0 ) −1/2 a positive constant. The general solution for the total density of energy in interactive systems with this interaction is ρ = (ρ 01 a −3 + ρ 02 a −3ν ) 2 . Since this kinetic function transforms the equation (9a) in 3H 2 = ρ = −xV = V 0φ 2 /φ 2 , the exotik field is φ = φ 0 a 3ν , φ 0 = φ (a = 1). The corresponding global EoS, ω = −1 − ρ /ρ can be written as and shows that the global system interpolates between a stiff fluid and a fluid with ω asym = 1 − 1/νφ 0 ρ 02 in either direction, depending on whether ν is less than or greater than 1. The interesting case is ν < 1, with stiff behavior at early times, that even admits the crossing of the phantom divide line, at late times.

Conclusions
We have described cosmological systems composed by two interactive fluids in the dark sector, with constant equations of state ω 1 and ω 2 , by means of unified models that are controlled by non-canonical versions of a scalar field. The description produces a generalization of k-essence field we have called exotic k-essence or exotik. There are many kinds of exotiks according to the expression proposed for the fluid representing the dark energy ρ 2 . The anzats used here is a linear combination of the density of energy and the pressure of usual k-essences, parameterized by the constant values α and β . Then, the different representations or exotiks are labeled by the particular choice of the potential V and by those parameters, noting that the exotik description is reduced to a common k-essence when α = (1 + ω 1 )/(ω 1 − ω 2 ) and β = ω 1 /(ω 1 − ω 2 ). The derivation of the equations for the exotik field requires that the equations of state are constant, which is not the case in the modified holographic Ricci model. However, there is a simple relationship between the two EoS and the constants A and B in the DE expression ρ MHR 2 and so the overall description can be carried out according to the process indicated at the footnotes in the end of bibliography. The choice constant V = V 0 simplifies the evolution equation for exotik field and its first integral is useful to observe novel features compared to the usual k-essence. In the latter case it is proven that the kinetic functions that have a root F(x(t 0 )) = 0 have a dust behavior in the epoch around t 0 and also the same happens when a ν FF σ −1 x << 1, that is, at early times, for α∆ > 0. The exotik includes these cases but also allows us to include the time of radiation domination. Choosing the appropriate representation β = (3ω 1 − 1)/4(ω 1 − ω 2 ) and 2α(ω 1 − ω 2 ) < 3(1 + ω 1 ) the model interpolates between the radiation like era at early times and ω → −1 at late times. The choice V = V 0 φ −2 is also an interesting representation because, even discarding solutions with stationary global EoS, there are first integrals for the evolution equation of common k-essence and more generally for representations with α∆ = 1 + β ∆ , that is only one free parameter. The ability to model the cosmological behavior at early times still exists in this more complicated pose although here the realistic models should be conducted with strictly increasing kinetic functions if we want to avoid the phantom regime. Having chosen the representation, α and β , the particular kinetic function F or the interaction Q are defined as the solution of equation (17) for a given potential. For V = V 0 we have shown the "self-interaction" of the modified holographic Ricci dark energy models and have found the exotik kinetic function in the case Q = τρ /∆ where we have found that the phantom or non-phantom regime depends on the strength of the coupling constant τ not allowing the crossing of the phantom divide line. Also we have found the interaction associated with the Chimento function and the associated with a purely Chaplying gas. Interestingly, the two complex interactions have an easy resolution from the point of view of exotik representation. For V = V 0 φ −2 , we have found that the observed equivalence (in FRW and in Bianchi I backgrounds) between usual k-essence with F(x) = 1 + mx and usual quintessences with exponential potential is still observed at the level of associated interactions. Finally, we found the interaction Q ν corresponding to a simplification of the Chimento function, for which we express the exotik field and the global EoS as functions of the factor of scale showing that for ν < 1 the models admit the crossing of the PDL, at late times. These exotik representations allow leverage all the work done with k-essence unified models giving them an interactive justification inside the dark sector.