Abstract
We revise the dynamics of interacting vector-like dark energy, a theoretical framework proposed to explain the accelerated expansion of the universe. By investigating the interaction between vector-like dark energy and dark matter, we analyze its effects on the cosmic expansion history and the thermodynamics of the accelerating universe. Our results demonstrate that the presence of interaction significantly influences the evolution of vector-like dark energy, leading to distinct features in its equation of state and energy density. We compare our findings with observational data and highlight the importance of considering interactions in future cosmological studies.
Similar content being viewed by others
1 Introduction
As revealed by observational data, the accelerated expansion of the universe is one of the most perplexing phenomena in modern cosmology. The prevailing explanation for this cosmic acceleration lies in an enigmatic entity known as dark energy, which constitutes a significant portion of the universe’s energy content [1, 2]. Exploring this entity offers immense potential to unravel the mysteries surrounding the fundamental building blocks and evolution of our vast Universe. The enigma of dark energy lies in its ability to drive the universe’s accelerated expansion while eluding a comprehensive understanding of its physical origin and underlying mechanisms. This cosmic force permeates the fabric of space, counteracting the gravitational pull of matter and pushing galaxies farther apart. The precise nature of dark energy remains elusive, with candidates ranging from a cosmological constant, often associated with the energy of empty space [3, 4], to dynamical scalar fields evolving over cosmic time [5] and modified gravity theories [6, 7].
The simplest dark energy model, the \(\Lambda \)CDM model based on the cosmological constant and cold dark matter, faces theoretical issues like the fine-tuning problem tied to its energy scale, the cosmological constant problem, and the coincidence problem [8, 9]. Furthermore, notable discrepancies have emerged between the \(\Lambda \)CDM model predictions using cosmic microwave background (CMB) data and independent local measurements [10,11,12,13]. Specifically, the Hubble constant \(H_{0}\) from Planck data is \(4.0 \sigma \) to \(6.3 \sigma \) below local measurements [14], and the \(\sigma _8\) clustering amplitude is higher than values from low-redshift observations [13, 15,16,17]. Thus, these theoretical issues and recent observational findings have motivated the scientific community to explore and investigate alternative dark energy models. For instance, in the context of scalar field models and modified gravity theories, one can find in the literature non-minimally coupled scalar fields models [18,19,20,21,22,23,24,25], coupled dark energy [26,27,28,29,30,31], \(f(R,\phi )\) [32,33,34,35,36,37,38], and \(f(T,\phi )\) gravity [39,40,41,42,43,44,45,46], among others (see Refs. [2, 7, 47] and references therein).
Additionally, dark energy has not only been described using scalar fields but also through other types of fields, such as vector fields. In the context of generalized Proca theories, a massive vector field is introduced, which breaks the U(1) symmetry. These vector theories belong to the class of time-like vector models. The homogeneous version of this vector field has a non-zero temporal component that is a function of cosmic time. For the most general case, its dynamics exhibit an asymptotic de Sitter attractor [48,49,50,51,52,53]. However, when the field is in its canonical form and minimally coupled to gravity, the equation motion becomes trivial, lacking dynamics and rendering it incapable of acting as a source of dark energy [54]. Furthermore, space-like vector has also been explored as a means to model inflation and dark energy [54,55,56,57,58,59,60,61]. While space-like models are generally associated with generating a highly anisotropic universe, there are specific scenarios where this challenge can be overcome. For instance, one of these scenarios involves the assumption of many randomly oriented vector fields, which, on average, results in an isotropic cosmological background, as described in Refs. [62, 63]. Another method entails the consideration of three identical vector fields for each spatial direction, referred to as the ‘cosmic triad,’ a concept that also aligns with the background symmetry, as expounded upon [57].
Specifically, in Ref. [57], the author proposed a dark energy model based on three self-interacting vector fields minimally coupled to gravity, oriented in mutually orthogonal spatial directions, and sharing the same time-dependent length. This model, also known as the vector-like dark energy model, effectively drives the current accelerated expansion of the universe and introduces new tracking attractor solutions, making cosmic evolution insensitive to initial conditions [57] (see also Refs. [54, 58, 64, 65]). Furthermore, in Ref. [66], the authors examined the dark-like energy model using the Noether approach, while in Ref. [67], they investigated the unified first law within this model.
As an interesting extension of this latter model, the Interacting Vector-like Dark Energy (IVDE) model has been proposed. For instance, in Ref. [68], the authors constructed a vector-like dark energy model for which an interaction between dark energy and a background perfect fluid was assumed. They identified two types of cosmological coincidence problems: the first asks why we live in an epoch where dark energy and dust matter energy densities are comparable; the second ponders why we exist in an epoch with \(w_{de}<-1.\) In this way, they found that these cosmological coincidence problems can be alleviated in such IVDE models. Also, in Ref. [69], the predictions of the IVDE model were compared with observational H(z) data. In particular, they showed that the best IVDE models exhibit an oscillating feature for H(z) and the equation-of-state parameter (EoS), crossing \(-1\) around redshift \(z\sim 1.5.\)
The IVDE model offers a compelling avenue of exploration, providing a fresh approach to understanding the profound mysteries of dark energy. In this theoretical framework, dark energy is not an isolated entity but interacts with other cosmic components, such as dark matter, radiation, or neutrinos [70, 71]. These interactions, mediated by a vector-like field, can potentially influence the cosmic expansion history and the formation and evolution of large-scale structures.
Understanding the properties and dynamics of IVDE is essential for several reasons. To the present work, IVDE may offer an alternative explanation for the accelerated expansion of the universe that extends beyond the standard cosmological model. By examining specific interaction mechanisms and studying their effects on cosmic evolution, we can gain deeper insights into the nature of dark energy and its role in shaping the cosmos [72, 73]. Moreover, from a broader perspective in theoretical physics, IVDE has the potential to illuminate the complex interplay between dark energy and other fundamental forces and particles. For instance, exploring IVDE within the frameworks of quantum field theory and particle physics framework could reveal significant connections, thereby bridging the gaps between our understanding of the macroscopic universe and the microscopic realm of particles and fundamental interactions [74, 75].
Therefore, the investigation of IVDE represents a captivating and vital avenue to understanding the nature and properties of dark energy. Through an in-depth analysis of IVDE, we can strive towards a comprehensive understanding of the nature and properties of dark energy while also advancing our knowledge of fundamental physics and the evolution of the cosmos. This research promises to unlock new horizons in our quest to comprehend the underlying mechanisms that govern the vastness of our universe.
In this paper, we study the cosmological dynamics of IVDE, analyzing its effects on the cosmic expansion history and the thermodynamics of the accelerating universe. In particular, we have performed a detailed phase space analysis, assuming several different functional forms for the coupling function Q between vector-like dark energy and dark matter. We extend previous studies in the literature [68] by including not only linear functional forms but also nonlinear functions for Q in terms of the energy densities. In the case of nonlinear functions, we demonstrate that the Q function can exhibit sign changes throughout the cosmic evolution.
The paper is organized as follows: In Sect. 2, we introduce the IVDE models, present the total action of the model, and derive the field equations in a general space-time. We then obtain the cosmological equations and define the effective dark energy, along with its coupling to dark matter. In Sect. 3 we reformulate the complete set of cosmological equations using new phase-space variables to derive a closed autonomous system. Furthermore, a detailed phase-space analysis of the model is performed, where we identify the corresponding critical points and their stability conditions. In Sect. 4, we corroborate our analytical findings by numerically solving the field equations. In Sect. 5, we explore the thermodynamic aspects of the interacting scenario within the context of IVDE models, determining the evolution of the temperature of dark matter and dark energy as influenced by the coupling function Q. Finally, in Sect. 6, we summarize the results obtained.
2 Vector-like dark energy
The concept of “vector-like dark energy” refers to a “cosmic triad,” as described in reference [57]. This triad consists of three identical vectors oriented in mutually orthogonal directions, ensuring the preservation of isotropy. In line with the approach presented in reference [57], we examine the scenario where vector-like dark energy is minimally coupled to gravity. The corresponding action can be expressed as follows:
where \(\kappa ^2=8\pi G,\) with G being the gravitational constant, \({F^a}_{\mu \nu }=\partial _\mu A^a_\nu -\partial _\nu A^a_\mu \) and \({A^a}^2=g^{\mu \nu }A^a_\mu A^a_\nu .\) \(S_m\) is the action of matter \(S_r\) is the action of radiation. It is worth noting that the superscript a indicates each vector field constituting the cosmic triad.
Then, varying the action (1) with respect to the metric, we obtain the following equations:
And, varying with respect to the cosmic triad \(A^a_\mu \) we find the equation of motion
Below, we detail the basic equations for the cosmic triad in a cosmological background.
2.1 Cosmological dynamics
We consider a spatially flat Friedmann–Lemaître–Robertson–Walker (FLRW) universe with metric:
where a is the scale factor, a function of the cosmic time t. We also assume the following ansatz:
Thus, the modified Friedmann equations are given by
and the motion equation of the vector field A:
In the above, \(V\left( A^2\right) \) represents the scalar potential, \(H\equiv {\dot{a}}/a\) stands for the Hubble rate, where a dot denotes the derivative with respect to time, and a comma indicates derivative with respect to \(A^2.\) Furthermore, the functions \(\rho _i\) and \(p_i,\) with \(i=m,r\) denoting either non-relativistic matter (comprising cold dark matter and baryons) or radiation, respectively, serve as the energy and pressure densities. It is important to note that in the equations mentioned earlier, we have already incorporated the respective barotropic equations of state, namely, \(w_m=p_m/\rho _m=0\) and \(w_r=p_r/\rho _r=1/3.\)
Following Ref. [2] one can rewrite the Friedmann equations (6) and (7) in their standard form as
Therefore, we can define the effective energy and pressure densities of the model as:
On the other hand, from the conservation of the total energy-momentum tensor we have
each component (matter, radiation, and dark energy) satisfy (13) separately in a scenario without interaction. And, it is worth noticing that matter and radiation have standard behavior \(\rho _{m}\sim a^{-3}\) and \(\rho _{r}\sim a^{-4},\) respectively.
Additionally, we can define the effective dark energy EoS as
It is also convenient to introduce the total equation-of-state parameter as
which is related to the deceleration parameter q through
and, hence, it is easy to see that the acceleration of the Universe occurs for \(q<0\) or, equivalently, when \(w_{tot}<-1/3.\)
Finally, another set of cosmological parameters that can be introduced is the standard density parameters, which are defined as
and satisfy the constraint equation
In what follows, we extend the previous analysis to include the interaction between vector-like dark energy and dark matter [76].
2.2 Interacting dark energy
We assume that the vector-like dark energy interacts with the dark matter through a phenomenological interaction term Q, according to [76]
which preserves the total energy conservation law. By assuming this interaction between the dark components, the equation of motion (8) should be changed when \(Q\ne 0\) and a new term due to Q will appear in its r.h.s. [68].
In the following section, we conduct a thorough phase-space analysis for this dark energy model. Specifically, we derive the corresponding autonomous system based on the set of cosmological equations (9), (10), (19), (20) and (21).
3 Phase-space analysis
We introduce the following set of dimensionless variables [2]:
and the constraint equation:
Therefore, we obtain the dynamical system
where \(N\equiv \ln {a}\) is the number of e-folds, and it is used as the temporal parameter since it is a function of the cosmic scale factor. Using the above set of phase space variables, we can also write:
The equation of state (EoS) of dark energy \(w_{de}=p_{de}/\rho _{de}\) can be rewritten as
whereas the total EoS becomes
To obtain an autonomous system from the dynamical system (24), we need to define the potentials for the vector field. From now we concentrate on the exponential potential \(V(A^2)\sim e^{-\kappa ^2\lambda A^2}.\) Below, we study the critical points and their stability properties for various interaction cases.
In the literature, various models have been explored where the function Q, which represents the interaction term, is considered as a function of energy densities and the Hubble parameter [76, 77]. In this paper, we use three interacting scenarios: linear, non-linear, and a sign-changeable [78].
3.1 Critical points
In this section, we obtain the critical points from the conditions \(\text {d}{x}/\text {d}{N}=\text {d}{y}/\text {d}{N}=\text {d}{\varrho }/\text {d}{N}=\text {d}{u}/\text {d}{N}=0\) [2]. Where, if we consider the definition of each dynamical variable (22), the critical points that are physically possible are given by \(y_c\ge 0,\ \varrho _c\ge 0\) and \(u_c\ge 0.\)
Consequently, we find the critical points for each interaction term:
3.1.1 Case I: \(Q = 3 \alpha H \rho _m\)
Using the above set of phase space variables, we can rewrite the interaction term as:
Critical points of the system (24) with interaction case (30) are shown in Table 1 and the values of their cosmological parameters in Table 2. Also, in this subsection and so on, we define \(\Omega _{de}^{(r)}\) and \(\Omega _{de}^{(m)}\) as the amount of dark energy during a radiation and dark matter domination era, respectively.
The critical point \(a_R\) corresponds to a scaling radiation era, where \(\Omega ^{(r)}_{de}=u_{c}^2.\) When \(u_{c}=0,\) we obtain the radiation-dominated solution with \(\Omega _r=1\) and \(w_{de}=w_{tot}=1/3.\) For \(u_{c}\ne 0 ,\) this point represents a scaling radiation era. Thus, to satisfy the early constraint imposed by the physics of big bang nucleosynthesis (BBN) and ensure \(\Omega ^{(r)}_{de}<0.045\) [79, 80], we need to have \(u_{c}<0.122.\)
The critical point c represents a dark energy-dominated solution with a de Sitter EoS, where \(w_{de}=w_{tot}=-1.\) As a result, this critical point yields accelerated expansion for all parameter values.
On the other hand, for \(\alpha =0,\) the critical point labeled as \(b_M\) represents a matter-dominated era, with \(\Omega _m = 1\) and \(w_{de} = w_{tot} = 0.\) In this case, the energy density of matter dominates the universe, and the EoS parameter for dark energy \((w_{de}=1/3)\) and the total EoS parameter \((w_{tot})\) is zero. For \(\alpha \ne 0 ,\) we have a scaling matter era which is constrained to satisfy \(\Omega _{de}^{(m)}< 0.02\) \((95\%\) CL), at redshift de \(z\approx 50,\) according to CMB measurements [81]. Thus, this leads us to the constraint \(0<\alpha <0.01,\) which is also compatible with \(0<\Omega _m<1.\)
3.1.2 Case II: \(Q = 3\beta H\dfrac{{\rho _m}^2}{\rho _m+\rho _{de}}\)
With the phase space variables mentioned earlier, we express the interaction term in the following form
The critical points of the system described in Eq. (24) under the interaction case given by Eq. (31) are presented in Table 3, along with the corresponding values of their cosmological parameters in Table 4.
The critical point labeled as \(d_R\) corresponds to a scaling radiation era, characterized by \(\Omega ^{(r)}_{de}=u_{c}^2.\) When \(u_{c}=0,\) we obtain the radiation-dominated solution with \(\Omega _r=1\) and \(w_{de}=w_{tot}=1/3.\) For \(u_{c}\ne 0 ,\) this point describes a scaling radiation era. To satisfy the early constraint imposed by the physics of big bang nucleosynthesis (BBN) and ensure \(\Omega ^{(r)}_{de}<0.045\) [79, 80], the condition \(u_{c}<0.122\) must be met.
On another note, for \(\beta =0,\) the critical point denoted as \(e_M\) represents a matter-dominated era, where \(\Omega _m = 1\) and \(w_{de} = w_{tot} = 0.\) In this scenario, the energy density of matter dominates the universe, resulting in a dark energy EoS parameter of \(w_{de}=1/3\) and a total EoS parameter of \(w_{tot}=0.\) For \(\beta \ne 0 ,\) we have a scaling matter era. Applying the constraint \(\Omega _{de}^{(m)}< 0.02\) \((95\%\) CL), at redshift de \(z\approx 50,\) according to CMB measurements [81], we obtain \(0<\beta <0.01,\) which is also compatible with \(0<\Omega _m<1\).
Moreover, the critical point labeled as f signifies a dark energy-dominated solution characterized by a de Sitter EoS, with \(w_{de}=w_{tot}=-1.\) Consequently, this critical point leads to accelerated expansion regardless of the parameter values.
Finally, the points \(f_1\) and \(f_2\) are not physically viable, because the fractional density is constrained by \(0<\Omega _m<1,\) meaning that, \(0<-1/\beta <1\) or \(\beta <-1,\) but \(\beta \) must be positive and small to obtain a matter-dominated era as described by point \(e_{M}\). Therefore, aiming to reproduce the standard thermal history of the universe, points \(f_1\) and \(f_2\) are identified as unphysical.
3.1.3 Case III: \(Q = \dfrac{3}{2\kappa ^2}\eta H^3 (1+\Omega _r-3\Omega _{de})\Omega _m\)
Given the above set of phase space variables, we express the interaction term as follows:
Critical points of the system (24) with interaction case (32) are shown in Table 5 and the values of their cosmological parameters in Table 6
The critical point denoted as \(g_R\) corresponds to a scaling radiation era, where \(\Omega ^{(r)}_{de}=u_{c}^2.\) When \(u_{c}=0,\) it gives rise to the radiation-dominated solution with \(\Omega _r=1\) and \(w_{de}=w_{tot}=1/3.\) The case \(u_{c}\ne 0 ,\) corresponds to a scaling radiation era. Therefore, to meet the early constraint imposed by the physics of big bang nucleosynthesis (BBN) and ensure \(\Omega ^{(r)}_{de}<0.045\) [79, 80], the condition \(u_{c}<0.122\) must be satisfied.
Conversely, when \(\eta =0,\) the critical point labeled as \(h_{M}\) represents a matter-dominated era, characterized by \(\Omega _m = 1\) and \(w_{de} = w_{tot} = 0.\) In this scenario, the energy density of matter dominates the universe, leading to a dark energy EoS parameter of \(w_{de}=1/3.\) In contrast, the total equation of the state parameter remains at \(w_{tot}=0.\) When \(\eta \ne 0,\) this point represents a scaling matter era. Thus, from constraint \(\Omega _{de}^{(m)}< 0.02\) \((95\%\) CL), at redshift de \(z\approx 50,\) according to CMB measurements [81], we get \(-0.04<\eta <0,\) which is also compatible with \(0<\Omega _m<1\).
Additionally, the critical point i corresponds to a dark energy-dominated solution with a de Sitter EoS, where \(w_{de}=w_{tot}=-1.\) Consequently, this critical point results in accelerated expansion for all parameter values.
Finally, points \(i_{1}\) and \(i_{2}\) are not physically viable. For these points, the physical condition \(0<\Omega _{m}<0\) implies \(\eta <-3\) or \(\eta >6,\) but this result is not compatible with a dark matter-dominated era as described by point \(h_{M}\) \((-0.04<\eta <0)\). Hence, to reproduce the standard thermal history of the universe, we designate the points \(i_1\) and \(i_2\) as unphysical.
3.2 Stability of critical points
To study the stability of the critical points, we introduce time-dependent linear perturbations denoted as \(\delta x,\) \(\delta y,\) \(\delta \varrho ,\) and \(\delta u\) around each critical point. These perturbations take the form of \(x=x_c+\delta x,\) \(y=y_c+\delta y,\) \(\varrho =\varrho _c+\delta \varrho ,\) and \(u=u_c+\delta u.\) Substituting these expressions into the autonomous system (24) and linearizing the equations, we obtain the linear perturbation matrix \({\mathcal {M}},\) as outlined in [2]. The eigenvalues of \({\mathcal {M}},\) denoted as \(\mu _1,\) \(\mu _2,\) \(\mu _3,\) and \(\mu _4,\) evaluated at each fixed point, determine the stability conditions for those points. Typically, the classification of stability properties proceeds as follows: (i) A stable node exists when all the eigenvalues are negative, (ii) An unstable node emerges when all the eigenvalues are positive, (iii) A saddle point is characterized by having one, two, or three of the four eigenvalues as positive and the others as negative, (iv) A stable spiral is observed when the determinant of \({\mathcal {M}}\) is negative, and the real part of all the eigenvalues is negative. Points classified as stable nodes or stable spirals are referred to as attractor points, and these fixed points are reached during the cosmic evolution of the Universe, regardless of the initial conditions of the system, as long as they belong to the attraction basin of the critical point. In the following lines, we present the eigenvalues and stability conditions for each critical point in every interaction case.
3.2.1 Case I: \(Q = 3 \alpha H \rho _m\)
-
Point \(a_R\) has the eigenvalues
$$\begin{aligned} \mu _1=0, \ \mu _2=-1, \ \mu _3=2, \ \mu _4=1-3\alpha , \end{aligned}$$(33)then this point is always unstable for all the values of \(\alpha .\)
-
Point \(b_M\) has the eigenvalues
$$\begin{aligned}{} & {} \mu _1=-1, \ \mu _2=\frac{1}{2}(-1+3\alpha ), \ \mu _3=\frac{3}{2}(\alpha +1), \nonumber \\{} & {} \mu _4=-1+3\alpha , \end{aligned}$$(34)which is a saddle point when \(\alpha >\frac{1}{3}\) or \(-1<\alpha <\frac{1}{3}.\) On the other hand, it becomes a stable node when \(\alpha <-1.\) However, this point cannot account for the current accelerated expansion of the Universe.
-
Point c has the eigenvalues
$$\begin{aligned}{} & {} \mu _1=-2, \ \mu _2=-3(\alpha +1),\nonumber \\{} & {} \mu _{3,4}=\tfrac{-3\lambda \pm \sqrt{16\lambda -7\lambda ^2}}{2\lambda }. \end{aligned}$$(35)This is a de-Sitter solution, ensuring accelerated expansion for all parameter values. We observe that it is a stable node when \(\alpha >-1\) and \(1<\lambda \le \frac{16}{7}.\) Finally, this point never exhibits stable spiral behavior (Fig. 1).
3.2.2 Case II: \(Q = 3\beta H\dfrac{{\rho _m}^2}{\rho _m+\rho _{de}}\)
-
Point \(d_R\) has the eigenvalues
$$\begin{aligned} \mu _1=2, \ \mu _2=-1, \ \mu _3=1, \ \mu _4=0, \end{aligned}$$(36)which tells us that it is always an unstable node.
-
Point \(e_M\) has the eigenvalues
$$\begin{aligned}{} & {} \mu _1=-\tfrac{1}{2(1+3\beta )}, \ \mu _{2,3}=-1, \ \mu _4=\tfrac{3(1+4\beta )}{2(1+3\beta )}, \end{aligned}$$(37)which is a saddle point when \(\beta <-\frac{1}{3}\) or \(\beta >-\frac{1}{4}.\) On the other hand, it becomes a stable node when \(-\frac{1}{3}<\beta <-\frac{1}{4}.\) However, this point cannot account for the current accelerated expansion of the Universe.
-
Point f has the eigenvalues
$$\begin{aligned} \mu _1=-2, \ \mu _2=-3, \ \mu _{3,4}=\tfrac{-3\pm \sqrt{-7+\frac{16}{\lambda }}}{2}. \end{aligned}$$(38)This is a de-Sitter solution, hence ensuring accelerated expansion for all parameter values. We observe that it is a stable node when \(1<\lambda \le \frac{16}{7}.\) Finally, this point never exhibits stable spiral behavior.
-
Stability of points \(f_1\) and \(f_2\) is given by the solution of the characteristic polynomials (Fig. 2).
3.2.3 Case III: \(Q = \dfrac{3}{2\kappa ^2}\eta H^3 (1+\Omega _r-3\Omega _{de})\Omega _m\)
-
Point \(g_R\) has the eigenvalues
$$\begin{aligned} \mu _1=0, \ \mu _2=-1, \ \mu _3=2, \ \mu _4=1+\eta -6\eta {u_c}^2,\nonumber \\ \end{aligned}$$(39)which is always an unstable node.
-
Point \(h_M\) has the eigenvalues
$$\begin{aligned} \mu _1=-1, \ \mu _{2}=\tfrac{1-\eta }{3\eta -2}, \ \mu _{3}=\tfrac{3-5\eta }{2-3\eta }, \ \mu _4=\eta -1, \end{aligned}$$(40)which acts as a saddle point when \(\frac{2}{3}<\eta <1,\) \(\eta >1,\) or \(\eta <\frac{3}{5}.\) Conversely, it behaves as a stable node when \(\frac{3}{5}<\eta <\frac{2}{3}.\) Nonetheless, this point cannot explain the current accelerated expansion of the Universe.
-
Point i has the eigenvalues
$$\begin{aligned} \mu _1=-2, \ \mu _2=-3-\eta , \ \mu _{3,4}=\tfrac{-3\pm \sqrt{-7+\frac{16}{\lambda }}}{2}. \end{aligned}$$(41)This is a de-Sitter solution, ensuring accelerated expansion for all parameter values. We observe that it is a stable node when \(\eta >-3\) and \(1<\lambda \le \frac{16}{7}.\) Finally, this point never exhibits stable spiral behavior.
-
Stability of points \(i_1\) and \(i_2\) is given by the solution of the characteristic polynomials (Fig. 3).
A summary of this section is shown in Table 7, which contains the principal properties of critical points.
4 Numerical results
In this section, we aim to numerically solve the autonomous system represented by Eq. (24), associated with a set of cosmological equations (9)–(10). Our objective is to explore the characteristics of our model (for each interacting term Q) to explain the ongoing accelerated expansion of the Universe. Subsequently, we compare the predicted results with the most recent measurements of the Hubble parameter H(z), Appendix A. We consider various parameter values and initial conditions with the aim of studying the behavior of cosmological parameters throughout the thermal history of the universe at intermediate redshifts. By investigating these several different parameter values and initial conditions, we can identify the range of potential evolutions for these cosmological parameters, all converging towards a common final-state attractor point.
4.1 Case I: \(Q = 3 \alpha H \rho _m\)
We have presented the numerical results found for this interacting model in Figs. 1, 4, 5, 6 and 7.
In Fig. 1, we verify that the model provides a dark energy-dominated solution with attractor behavior. The attractor is a one-parameter solution, and thus, one can have different locations for this fixed point in the field space. This implies varying temporal durations for the transitions between the decelerating and accelerating phases. However, this is restricted by cosmological observations [82]. In Fig. 4, we depict the behavior of the energy density for radiation, matter, and dark energy. Remarkably, one can observe the scaling behavior of dark energy during the radiation and matter-dominated eras. This is an essential feature of the model because scaling solutions provide a natural mechanism to alleviate the energy scale problem of dark energy [41, 83, 84]. In Fig. 5, we show the evolution of both the EoS of dark energy and the total EoS.
Additionally, we added the corresponding curve associated with the \(\Lambda \)CDM model to contrast the predictions of our model. During the radiation and matter-dominated eras, the EoS of dark energy behaves like a radiation field, which is an expected result since the source of dark energy is a vector field [57]. Thus, in Fig. 6, we observe a slight discrepancy in the deceleration-acceleration transition redshift compared to the \(\Lambda \)CDM result. This latter result is validated in Fig. 7, where it shows the evolution of the Hubble rate H(z) as a function of the redshift, along with the corresponding results from \(\Lambda \)CDM. In the same plot, we have depicted the relative difference for the \(\Lambda \)CDM model, showing that our results are compatible with observations (for more detail, see Table 8).
4.2 Case II: \(Q = 3\beta H\dfrac{{\rho _m}^2}{\rho _m+\rho _{de}}\)
For this ansatz, we have shown the numerical results in Figs. 2, 8, 9, 10 and 11.
In Fig. 2 one can see that the system has an attractor point, which is a dark energy-dominated solution. Furthermore, this attractor point has a nature of cosmological constant with an EoS \(w_{de}=-1.\) However, the position of this fixed point in the phase space depends on the parameter \(\lambda .\) Thus, the transition time between the dark matter-dominated era and late-time acceleration depends on it, too. Let us remember that this transition time is constrained by current cosmological observations [82]. In Fig. 8, we can observe the scaling regimes during the dark-matter and radiation-dominated epochs. At early times, the dark energy component reaches higher energy scales, alleviating the energy scale problem of dark energy. In Fig. 9, we depict the behavior of the EoS of dark energy and total EoS for our model, comparing the corresponding results from \(\Lambda \)CDM.
Interestingly enough, the evolution of the total EoS of our model is very close to that of \(\Lambda \)CDM, with some slight differences during the dark matter-dark energy transition time. However, the EoS of dark energy behaves like a radiation field at early times and like a cosmological constant at late times, also, in Fig. 10, one can observe the decelerating-accelerating transition time as described by the evolution of the deceleration parameter. We found that this transition time is slightly displaced towards smaller redshift values than the result from \(\Lambda \)CDM. Nevertheless, this result is still compatible with observational data [82].
Finally, in Fig. 11, we numerically corroborated that our model is compatible with observational data. The evolution of the Hubble rate is consistent with Hubble rate data found in the literature [85, 86].
4.3 Case III: \(Q = \dfrac{3}{2\kappa ^2}\eta H^3 (1+\Omega _r-3\Omega _{de})\Omega _m\)
Finally, for our third model, the numerical results are shown in Figs. 3, 12, 13, 14 and 15.
In Fig. 3, we show the evolution curves in the phase space for this ansatz. As for previous models, we have an attractor fixed point, which is a cosmological constant solution. This attractor solution with accelerated expansion has a variable position in the phase space, which depends on the model parameter related to the slope of the vector’s potential. In Fig. 12, we show the scaling behavior for early times.
On the other hand, in Fig. 13 we depict the evolution of the EoS of dark energy and the total EoS. As before, the effective dark energy fluid behaves as radiation field density at early times, whereas it behaves as a quintessence field at late times, in Fig. 14, we observe that this model provides a deceleration-acceleration transition redshift closer to the value from the \(\Lambda \)CDM model, and therefore, it becomes compatible with observations [82]. This result is also corroborated in Fig. 15, where we plot the theoretical Hubble rate from our model along with the current Hubble rate data [85, 86].
5 Thermodynamics
The Universe, as a whole, can be seen as a classical thermodynamic system composed of interacting fluids. The temperature is defined through the Gibbs equation
where T is the temperature of the system and S is total entropy per comoving volume and \(N=n V=const\) with \(V=a^3,\) and n is the number density.
This equation can be written as
Thus, the integrability condition
can allow us to obtain
Since \(T=T(n,\rho ),\) and using the conservation laws
where \(\rho \) and p are the total energy and pressure densities, we obtain
Using the relation (45) we get
In the case of a barotropic perfect fluid with pressure \(p=w(a) \rho \) we obtain
For two interacting fluids (dark matter and dark energy) one can write
where it has been defined
and then Eq. (50) gives us
Thus, one can verify that for non-relativistic matter \(w_{m}=0,\) the temperature is constant in the absence of coupling. So, after specifying the coupling between dark matter and dark energy, we can evolve the temperature for both components as functions of redshifts. For the present interacting vector-like dark energy model, the time evolution of the effective EoS of dark energy \(w^{eff}_{de}\) depends on both the magnitude of the vector field and the coupling to matter. Therefore, temperature behavior is determined by the dynamics of the vector field. In Figs. 16, 17, and 18, we depict the temperature behavior of matter and dark energy, parameterized in terms of the ratios \(T_{de}(z)/T_{de}^{(0)}\) and \(T_{m}(z)/T_{m}^{(0)},\) respectively, as functions of redshift z for different initial conditions. Here, we define \(T_{m}^{(0)}\equiv T_{m}(z=0)\) and \(T_{de}^{(0)}\equiv T_{de}(z=0).\) Also, in Fig. 19, we depict the behavior of the coupling function Q(z) for the three models studied. From Figs. 16, 17, and 18, it is evident that the temperature of matter increases very slowly, while the temperature of dark energy rises more rapidly. Consequently, when we combine these findings with Fig. 19, we observe an energy transference from dark energy to dark matter, indicating that dark energy possesses a negative heat capacity. Furthermore, it is verified that the second law of thermodynamics is satisfied during the regime \(T_{m}<T_{de}\) for \(Q>0.\) From the Gibbs equation (42) and the continuity equations (19) and (20) one can demonstrate that \(T_{de} dS_{de}=-Q Vdt\) and \(T_{m} dS_{m}=Q Vdt,\) where \(S_{m}\) and \(S_{de}\) represent the entropy of matter and dark energy, respectively. Consequently, we can verify that \(d(S_{de}+S_{m})=dS_{tot}=QV(1/T_{m}-1/T_{de})dt,\) and the second law \(d(S_{de}+S_{m})>0,\) requires \(T_{m}<T_{de}\) if \(Q>0,\) or \(T_{m}>T_{de}\) if \(Q<0\) [87]. In Figs. 16, 17, and 18, we observe that \(T_{de}(z)/T_{de}^{(0)}<T_{m}(z)/T_{m}^{(0)}.\) Thus, for \(Q>0,\) the validity of the second law requires \(T_{m}<T_{de}\) [87], and therefore \(T_{m}^{(0)}<T_{de}^{(0)}\) (see Appendix B).
In the case of model III, we observe a sign change in Q, which suggests a change in the direction of the energy flux. This could imply a corresponding sign change in the heat capacity of dark energy, which must be in concordance with the temperature behavior in Fig. 18 (see Refs. [87, 88]).
6 Concluding remarks
In the present paper, we studied the cosmological dynamics of interacting vector-like dark energy. The vector-like dark energy refers to the “cosmic triad”, which consists of three identical vectors oriented in mutually orthogonal directions and thus preserving the space isotropy [57]. Furthermore, we assumed an interaction between dark energy and dark matter. This interaction is described by the function Q, which is a function of the energy densities of both dark energy and dark matter, as well as the Hubble rate. Thus, we have studied several different models for Q [76].
We have shown that the cosmic triad can explain the dynamics of dark energy, providing a dark energy-dominated solution with accelerated expansion. These dark energy solutions are attractors fixed points, and thus, the system will reach them over a wide range of initial conditions. Nevertheless, one can observe that the whole cosmological evolution is altered by the presence of interaction between the effective dark energy fluid and dark matter. In particular, we obtained new scaling solutions of radiation and dark matter, which originate from both the pure dynamics of the vector field and the matter-vector interaction function Q.
It is important to note that although the critical point describing the dark energy-dominated era is an attractor, reachable through a wide set of initial conditions, fine-tuning is still required. This is to ensure consistency with the current estimations of cosmological parameters at redshift \(z=0,\) as well as to accurately reproduce the thermal history of the universe in line with observational data. The need for fine-tuning arises because the attractor critical point is only asymptotically achieved in the future, not precisely at \(z=0.\) To address this issue of finite-tuning, also known as the cosmological coincidence problem [2], the existence of an attractor scaling solution featuring accelerated expansion may be necessary. However, obtaining such solutions is inherently difficult. Not all dark energy models offer this kind of solution [19, 89, 90], and in some cases, even when these solutions are present, the model may not successfully replicate the dark matter era [91].
As the scaling solutions contribute to the existence of small amounts of dark energy in the radiation and matter era, they can give rise to notable physical outcomes [83, 84]. This specifically changes the Hubble rate during that era, leading to adjustments in the theoretical forecasts for the abundances of primordial light elements [79, 80]. Furthermore, the shape of the Cosmic Microwave Background (CMB) anisotropies spectrum is profoundly affected by such a scaling field [80]. The influence of this scaling field also extends to the development of the universe’s large-scale structures, affecting the evolution of cosmic entities like galaxies and galaxy clusters [91].
We also have studied the thermodynamics of our model, considering the interaction function Q. Using the Gibbs equation applied to an expanding universe, one can obtain general expressions for the temperature of dark matter and dark energy as functions of redshift in the presence of interaction between them [87]. We have analyzed three types of interaction functions Q that depend on the energy density of dark energy and dark matter. Since the effective dark energy depends on the vector field and its dynamics, the interaction function Q does well. In particular, we found for each interaction model the expressions for Q in terms of the dynamical variables. By numerically solving the background equations and using the expressions for Q in terms of the phase-space variables, we have depicted the behavior of the temperatures of dark energy and dark matter. Our results showed us that the temperature of matter increases very slowly, whereas the temperature of dark energy rises more quickly. Thus, we corroborated that there is a transference of energy from dark energy to dark matter, as dark energy has a negative heat capacity. In this way, we verified that the second law of thermodynamics is satisfied during the regime \(T_{m}<T_{de}\) for \(Q>0\) [87], provided that \(T^{(0)}_{m}<T^{(0)}_{de}\) However, it is important to note that measuring the temperature of dark matter and dark energy is extremely challenging [87, 92, 93]. Currently, the temperatures of these components remain unknown. Additionally, developing a suitable thermodynamic framework for the dark energy sector continues to be an area of active research due to the existing gaps in our understanding of its intrinsic nature. It is conceivable that future astronomical and cosmological observations will enhance our understanding of these elusive components , thereby enriching our knowledge of their thermodynamic evolution [94, 95]. The analysis of the models studied through the interaction term (as shown in Fig. 19) determines that the models include a change of sign for the dark energy interaction function. Since the interaction terms change their signs during evolution, our results indicate that today, dark energy is transferred to dark matter, but in the past, the transfer was the opposite [78, 96, 97].
Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors’ comment: The present work is a theoretical study, and therefore no experimental data has been listed.]
References
P.A. Ade, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. Banday, R. Barreiro, J. Bartlett, N. Bartolo et al., Planck 2015 results-xiii. Cosmological parameters. Astron. Astrophys. 594, A13 (2016)
E.J. Copeland, M. Sami, S. Tsujikawa, Dynamics of dark energy. Int. J. Mod. Phys. D 15, 1753 (2006). arXiv:hep-th/0603057
S.M. Carroll, The cosmological constant. Living Rev. Relativ. 4, 1 (2001). arXiv:astro-ph/0004075
A. Padilla, Lectures on the Cosmological Constant Problem (2015). arXiv:1502.05296 [hep-th]
C. Wetterich, The cosmological constant and the time variation of the fundamental constants. Astron. Astrophys. 301, 321 (2015)
S. Tsujikawa, Modified gravity models of dark energy. Lect. Notes Phys. 800, 99 (2010). arXiv:1101.0191 [gr-qc]
S. Nojiri, S.D. Odintsov, Introduction to modified gravity and gravitational alternative for dark energy. eConf C0602061, 06 (2006). arXiv:hep-th/0601213
S. Weinberg, The cosmological constant problem. Rev. Mod. Phys. 61, 1 (1989)
J. Martin, Everything you always wanted to know about the cosmological constant problem (but were afraid to ask). Comptes Rendus Phys. 13, 566 (2012). arXiv:1205.3365 [astro-ph.CO]
E. Abdalla et al., Cosmology intertwined: a review of the particle physics, astrophysics, and cosmology associated with the cosmological tensions and anomalies, in 2022 Snowmass Summer Study JHEAp 34, 49–211 (2022). arXiv:2203.06142 [astro-ph.CO]. https://doi.org/10.1016/j.jheap.2022.04.002
E. Di Valentino et al., Snowmass 2021—letter of interest cosmology intertwined II: the Hubble constant tension. Astropart. Phys. 131, 102605 (2021). arXiv:2008.11284 [astro-ph.CO]
E. Di Valentino et al., Cosmology intertwined III: \(f\sigma _8\) and \(S_8\). Astropart. Phys. 131, 102604 (2021). arXiv:2008.11285 [astro-ph.CO]
L. Heisenberg, H. Villarrubia-Rojo, J. Zosso, Simultaneously solving the H0 and \({\sigma }8\) tensions with late dark energy. Phys. Dark Universe 39, 101163 (2023). arXiv:2201.11623 [astro-ph.CO]
A.G. Riess, The expansion of the universe is faster than expected. Nat. Rev. Phys. 2, 10 (2019). arXiv:2001.03624 [astro-ph.CO]
C. Heymans et al., KiDS-1000 cosmology: multi-probe weak gravitational lensing and spectroscopic galaxy clustering constraints. Astron. Astrophys. 646, A140 (2021). arXiv:2007.15632 [astro-ph.CO]
R.C. Nunes, S. Vagnozzi, Arbitrating the S8 discrepancy with growth rate measurements from redshift-space distortions. Mon. Not. R. Astron. Soc. 505, 5427 (2021). arXiv:2106.01208 [astro-ph.CO]
L. Heisenberg, H. Villarrubia-Rojo, J. Zosso, Can late-time extensions solve the H0 and \({\sigma }8\) tensions? Phys. Rev. D 106, 043503 (2022). arXiv:2202.01202 [astro-ph.CO]
U. Kasper, On the cosmological constant problem and inflation in the framework of nonminimal coupling. Nuovo Cim. B 103, 291 (1989)
L. Amendola, Scaling solutions in general nonminimal coupling theories. Phys. Rev. D 60, 043501 (1999). arXiv:astro-ph/9904120
J.-P. Uzan, Cosmological scaling solutions of nonminimally coupled scalar fields. Phys. Rev. D 59, 123510 (1999). arXiv:gr-qc/9903004
T. Chiba, Quintessence, the gravitational constant, and gravity. Phys. Rev. D 60, 083508 (1999). arXiv:gr-qc/9903094
N. Bartolo, M. Pietroni, Scalar tensor gravity and quintessence. Phys. Rev. D 61, 023518 (2000). arXiv:hep-ph/9908521
F. Perrotta, C. Baccigalupi, S. Matarrese, Extended quintessence. Phys. Rev. D 61, 023507 (1999). arXiv:astro-ph/9906066
R. Gannouji, D. Polarski, A. Ranquet, A.A. Starobinsky, Scalar–tensor models of normal and phantom dark energy. JCAP. 09, 016 (2006). arXiv:astro-ph/0606287
Z. Davari, V. Marra, M. Malekjani, Cosmological constrains on minimally and non-minimally coupled scalar field models. Mon. Not. R. Astron. Soc. 491, 1920 (2020). arXiv:1911.00209 [gr-qc]
J.R. Ellis, S. Kalara, K.A. Olive, C. Wetterich, Density dependent couplings and astrophysical bounds on light scalar particles. Phys. Lett. B 228, 264 (1989)
C. Wetterich, The Cosmon model for an asymptotically vanishing time dependent cosmological ‘constant’. Astron. Astrophys. 301, 321 (1995). arXiv:hep-th/9408025
T. Damour, K. Nordtvedt, Tensor–scalar cosmological models and their relaxation toward general relativity. Phys. Rev. D 48, 3436 (1993)
L. Amendola, Coupled quintessence. Phys. Rev. D 62, 043511 (2000)
E. Di Valentino, A. Melchiorri, O. Mena, S. Vagnozzi, Interacting dark energy in the early 2020s: a promising solution to the \(H_0\) and cosmic shear tensions. Phys. Dark Universe 30, 100666 (2020). arXiv:1908.04281 [astro-ph.CO]
E. Di Valentino, A. Melchiorri, O. Mena, S. Vagnozzi, Nonminimal dark sector physics and cosmological tensions. Phys. Rev. D 101, 063502 (2020). arXiv:1910.09853 [astro-ph.CO]
S. Tsujikawa, Matter density perturbations and effective gravitational constant in modified gravity models of dark energy. Phys. Rev. D 76, 023514 (2007). arXiv:0705.1032 [astro-ph]
A. De Felice, S. Mukohyama, S. Tsujikawa, Density perturbations in general modified gravitational theories. Phys. Rev. D 82, 023524 (2010). arXiv:1006.0281 [astro-ph.CO]
V. Faraoni, Cosmology in Scalar Tensor Gravity (Springer, Berlin, 2004)
H. Farajollahi, M. Setare, F. Milani, F. Tayebi, Cosmic dynamics in \(F(R,\phi )\) gravity. Gen. Relativ. Gravit. 43, 1657 (2011). arXiv:1005.2026 [physics.gen-ph]
F. Hammad, Density perturbations in f(R,\(\phi \)) gravity with an application to the varying-power-law model. Phys. Rev. D 96, 064006 (2017). arXiv:1709.02276 [gr-qc]
J. Solà Peracaula, A. Gomez-Valent, J. de Cruz Pérez, C. Moreno-Pulido, Brans–Dicke gravity with a cosmological constant smoothes out \(\Lambda \)CDM tensions. Astrophys. J. Lett. 886, L6 (2019). arXiv:1909.02554 [astro-ph.CO]
J. Solà Peracaula, A. Gómez-Valent, J. de Cruz Pérez, C. Moreno-Pulido, Brans–Dicke cosmology with a \(\Lambda \)-term: a possible solution to \(\Lambda \)CDM tensions. Class. Quantum Gravity 37, 245003 (2020). arXiv:2006.04273 [astro-ph.CO]
M. Hohmann, L. Järv, U. Ualikhanova, Covariant formulation of scalar-torsion gravity. Phys. Rev. D 97, 104011 (2018). arXiv:1801.05786 [gr-qc]
M. Gonzalez-Espinoza, G. Otalora, Generating primordial fluctuations from modified teleparallel gravity with local Lorentz-symmetry breaking. Phys. Lett. B 809, 135696 (2020). arXiv:2005.03753 [gr-qc]
M. Gonzalez-Espinoza, G. Otalora, Cosmological dynamics of dark energy in scalar-torsion \(f(T,\phi )\) gravity. Eur. Phys. J. C 81, 480 (2021). arXiv:2011.08377 [gr-qc]
M. Gonzalez-Espinoza, G. Otalora, J. Saavedra, Stability of scalar perturbations in scalar-torsion \(f(T,\phi )\) gravity theories in the presence of a matter fluid. JCAP. 10, 007 (2021). arXiv:2101.09123 [gr-qc]
M. Gonzalez-Espinoza, R. Herrera, G. Otalora, J. Saavedra, Reconstructing inflation in scalar-torsion \(f(T,\phi )\) gravity. Eur. Phys. J. C 81, 731 (2021). arXiv:2106.06145 [gr-qc]
L.K. Duchaniya, B. Mishra, J. Levi Said, Noether symmetry approach in scalar-torsion \(f(T,\phi )\) gravity. Eur. Phys. J. C 83, 613 (2023). arXiv:2210.11944 [gr-qc]
L.K. Duchaniya, S.A. Kadam, J. Levi Said, B. Mishra, Dynamical systems analysis in \(f(T,\phi )\) gravity. Eur. Phys. J. C 83, 27 (2023). arXiv:2209.03414 [gr-qc]
M. Gonzalez-Espinoza, G. Otalora, Y. Leyva, J. Saavedra, Phase-space analysis of torsion-coupled dilatonic ghost condensate (2023). arXiv preprint. arXiv:2306.03386 [gr-qc]
Y.-F. Cai, S. Capozziello, M. De Laurentis, E.N. Saridakis, f(T) teleparallel gravity and cosmology. Rep. Prog. Phys. 79, 106901 (2016). arXiv:1511.07586 [gr-qc]
A. De Felice, L. Heisenberg, R. Kase, S. Mukohyama, S. Tsujikawa, Y.-L. Zhang, Cosmology in generalized Proca theories. JCAP 06, 048 (2016). arXiv:1603.05806 [gr-qc]
A. De Felice, L. Heisenberg, R. Kase, S. Mukohyama, S. Tsujikawa, Y.-L. Zhang, Effective gravitational couplings for cosmological perturbations in generalized Proca theories. Phys. Rev. D 94, 044024 (2016). arXiv:1605.05066 [gr-qc]
S. Nakamura, R. Kase, S. Tsujikawa, Coupled vector dark energy. JCAP. 12, 032 (2019). arXiv:1907.12216 [gr-qc]
A. De Felice, S. Nakamura, S. Tsujikawa, Suppressed cosmic growth in coupled vector–tensor theories. Phys. Rev. D 102, 063531 (2020). arXiv:2004.09384 [gr-qc]
W. Cardona, J.L. Palacios-Córdoba, C.A. Valenzuela-Toledo, Scrutinizing coupled vector dark energy in light of data (2023). arXiv preprint. arXiv:2310.13877 [astro-ph.CO]
J.B. Orjuela-Quintana, C.A. Valenzuela-Toledo, Anisotropic k-essence. Phys. Dark Universe 33, 100857 (2021). arXiv:2106.06432 [gr-qc]
T. Koivisto, D.F. Mota, Vector field models of inflation and dark energy. JCAP 08, 021 (2008). arXiv:0805.4229 [astro-ph]
L.H. Ford, Inflation driven by a vector field. Phys. Rev. D 40, 967 (1989)
A.B. Burd, J.E. Lidsey, An analysis of inflationary models driven by vector fields. Nucl. Phys. B 351, 679 (1991)
C. Armendariz-Picon, Could dark energy be vector-like? JCAP 07, 007 (2004). arXiv:astro-ph/0405267
L.G. Gomez, Y. Rodriguez, Coupled multi-Proca vector dark energy. Phys. Dark Universe 31, 100759 (2021). arXiv:2004.06466 [gr-qc]
M. Gonzalez-Espinoza, G. Otalora, Y. Leyva, J. Saavedra, Dynamics of dark energy in a scalar–vector-torsion theory. Eur. Phys. J. Plus 138, 600 (2023). arXiv:2212.12071 [gr-qc]
Y. Rodríguez, A.A. Navarro, Non-abelian \(S\)-term dark energy and inflation. Phys. Dark Universe 19, 129 (2018). arXiv:1711.01935 [gr-qc]
J.C. Garnica, L.G. Gomez, A.A. Navarro, Y. Rodriguez, Constant-roll inflation in the generalized SU(2) Proca theory. Ann. Phys. 534, 2100453 (2022). arXiv:2109.10154 [gr-qc]
A. Golovnev, V. Mukhanov, V. Vanchurin, Vector inflation. JCAP 06, 009 (2008). arXiv:0802.2068 [astro-ph]
G. Otalora, A. Övgün, J. Saavedra, N. Videla, Inflation from a nonlinear magnetic monopole field nonminimally coupled to curvature. JCAP 06, 003 (2018). arXiv:1803.11358 [gr-qc]
T.S. Koivisto, N.J. Nunes, Inflation and dark energy from three-forms. Phys. Rev. D 80, 103509 (2009). arXiv:0908.0920 [astro-ph.CO]
R.C.G. Landim, Dynamical analysis for a vector-like dark energy. Eur. Phys. J. C 76, 480 (2016). arXiv:1605.03550 [hep-th]
Y. Zhang, Y.-G. Gong, Z.-H. Zhu, Noether symmetry approach in ‘cosmic triad’ vector field scenario. Class. Quantum Gravity 27, 135019 (2010). arXiv:0912.4766 [astro-ph.CO]
Y. Zhang, Y. Gong, Z.-H. Zhu, The unified first law in “cosmic triad’’ vector field scenario. Phys. Lett. B 700, 254 (2011). arXiv:1108.1046 [hep-th]
H. Wei, R.-G. Cai, Interacting vector-like dark energy, the first and second cosmological coincidence problems. Phys. Rev. D 73, 083002 (2006). arXiv:astro-ph/0603052
H. Wei, S.N. Zhang, Observational H(z) data and cosmological models. Phys. Lett. B 644, 7 (2007). arXiv:astro-ph/0609597
C. Wetterich, Phenomenological parameterization of quintessence. Phys. Lett. B 594, 17 (2004)
L.P. Chimento, A.S. Jakubi, D. Pavon, W. Zimdahl, Interacting quintessence solution to the coincidence problem. Phys. Rev. D 67, 083513 (2003)
C.G. Boehmer, G. Caldera-Cabral, R. Lazkoz, R. Maartens, Dynamics of dark energy with a coupling to dark matter. Phys. Rev. D 78, 023505 (2008)
A.A. Costa et al., Constraints on interacting vector-like dark energy from cosmological observations. JCAP 12, 050 (2021)
A. Morandi, M. Sun, Probing dark energy via galaxy cluster outskirts. Mon. Not. R. Astron. Soc. 457, 3266 (2016). arXiv:1601.03741 [astro-ph.CO]
A.A. Costa, X.-D. Xu, B. Wang, E. Abdalla, Constraints on interacting dark energy models from Planck 2015 and redshift-space distortion data. JCAP 01, 028 (2017). arXiv:1605.04138 [astro-ph.CO]
B. Wang, E. Abdalla, F. Atrio-Barandela, D. Pavon, Dark matter and dark energy interactions: theoretical challenges, cosmological implications and observational signatures. Rep. Prog. Phys. 79, 096901 (2016). arXiv:1603.08299 [astro-ph.CO]
S. Lepe, G. Otalora, J. Saavedra, Dynamics of viscous cosmologies in the full Israel–Stewart formalism. Phys. Rev. D 96, 023536 (2017). arXiv:1704.05625 [gr-qc]
F. Arévalo, A. Cid, Dynamics and statefinder analysis of a class of sign-changeable interacting dark energy scenarios. Eur. Phys. J. C 82, 946 (2022)
P.G. Ferreira, M. Joyce, Cosmology with a primordial scaling field. Phys. Rev. D 58, 023503 (1998). arXiv:astro-ph/9711102
R. Bean, S.H. Hansen, A. Melchiorri, Early universe constraints on a primordial scaling field. Phys. Rev. D 64, 103508 (2001). arXiv:astro-ph/0104162
P. Ade et al. (Planck), Planck 2015 results. XIV. Dark energy and modified gravity. Astron. Astrophys. 594, A14 (2016). arXiv:1502.01590 [astro-ph.CO]
N. Aghanim et al. (Planck), Planck 2018 results. VI. Cosmological parameters. Astron. Astrophys. 641, A6 (2020). arXiv:1807.06209 [astro-ph.CO]
I.S. Albuquerque, N. Frusciante, N.J. Nunes, S. Tsujikawa, New scaling solutions in cubic Horndeski theories. Phys. Rev. D 98, 064038 (2018). arXiv:1807.09800 [gr-qc]
J. Ohashi, S. Tsujikawa, Assisted dark energy. Phys. Rev. D 80, 103513 (2009). arXiv:0909.3924 [gr-qc]
S.-L. Cao, X.-W. Duan, X.-L. Meng, T.-J. Zhang, Cosmological model-independent test of \(\Lambda \) CDM with two-point diagnostic by the observational Hubble parameter data. Eur. Phys. J. C 78, 1 (2018). https://doi.org/10.1140/epjc/s10052-018-5796-y
O. Farooq, B. Ratra, Hubble parameter measurement constraints on the cosmological deceleration–acceleration transition redshift. Astrophys. J. Lett. 766, L7 (2013)
V.H. Cárdenas, D. Grandón, S. Lepe, Dark energy and dark matter interaction in light of the second law of thermodynamics. Eur. Phys. J. C 79, 357 (2019). arXiv:1812.03540 [astro-ph.CO]
S. Lepe, G. Otalora, Gauss–Bonnet dark energy on Hořava–Lifshitz cosmology. Eur. Phys. J. C 78, 331 (2018). arXiv:1801.04213 [gr-qc]
L. Amendola, M. Quartin, S. Tsujikawa, I. Waga, Challenges for scaling cosmologies. Phys. Rev. D 74, 023525 (2006). arXiv:astro-ph/0605488
G. Otalora, Scaling attractors in interacting teleparallel dark energy. JCAP 1307, 044 (2013). arXiv:1305.0474 [gr-qc]
L. Amendola, Coupled quintessence. Phys. Rev. D 62, 043511 (2000). arXiv:astro-ph/9908023
N. Mirabolfathi, Dark matter direct detection with low temperature detectors, in 48th Rencontres de Moriond on Electroweak Interactions and Unified Theories (2013), pp. 123–130. arXiv:1308.0044 [astro-ph.IM]
D. Obreschkow, C. Power, M. Bruderer, C. Bonvin, A robust measure of cosmic structure beyond the power-spectrum: cosmic filaments and the temperature of dark matter. Astrophys. J. 762, 115 (2013). arXiv:1211.5213 [astro-ph.CO]
P.-J. Wu, Y. Shao, S.-J. Jin, X. Zhang, A path to precision cosmology: synergy between four promising late-universe cosmological probes. JCAP 06, 052 (2023). arXiv:2202.09726 [astro-ph.CO]
Z.-W. Zhao, L.-F. Wang, J.-G. Zhang, J.-F. Zhang, X. Zhang, Probing the interaction between dark energy and dark matter with future fast radio burst observations. JCAP 04, 022 (2023). arXiv:2210.07162 [astro-ph.CO]
A. Cid, C. Rodriguez-Benites, M. Cataldo, G. Casanova, Bayesian comparison of interacting modified holographic Ricci dark energy scenarios. Eur. Phys. J. C 81, 31 (2021)
C. Rodriguez-Benites, M. Cataldo, M. Vásquez-Arteaga, Universe with holographic dark energy. Momento 1 (2020)
C. Zhang, H. Zhang, S. Yuan, S. Liu, T.-J. Zhang, Y.-C. Sun, Four new observational h (z) data from luminous red galaxies in the Sloan Digital Sky Survey data release seven. Res. Astron. Astrophys. 14, 1221 (2014)
J. Simon, L. Verde, R. Jimenez, Constraints on the redshift dependence of the dark energy potential. Phys. Rev. D 71, 123001 (2005)
M. Moresco, A. Cimatti, R. Jimenez, L. Pozzetti, G. Zamorani, M. Bolzonella, J. Dunlop, F. Lamareille, M. Mignoli, H. Pearce et al., Improved constraints on the expansion rate of the universe up to z 1.1 from the spectroscopic evolution of cosmic chronometers. J. Cosmol. Astropart. Phys. 2012(08), 006 (2012)
A.J. Cuesta, M. Vargas-Magaña, F. Beutler, A.S. Bolton, J.R. Brownstein, D.J. Eisenstein, H. Gil-Marín, S. Ho, C.K. McBride, C. Maraston et al., The clustering of galaxies in the SDSS-III baryon oscillation spectroscopic survey: baryon acoustic oscillations in the correlation function of LOWZ and CMASS galaxies in data release 12. Mon. Not. R. Astron. Soc. 457, 1770 (2016)
C. Blake, S. Brough, M. Colless, C. Contreras, W. Couch, S. Croom, D. Croton, T.M. Davis, M.J. Drinkwater, K. Forster et al., The WiggleZ Dark Energy Survey: joint measurements of the expansion and growth history at z \(<\) 1. Mon. Not. R. Astron. Soc. 425, 405 (2012)
A. Ratsimbazafy, S. Loubser, S. Crawford, C. Cress, B. Bassett, R. Nichol, P. Väisänen, Age-dating luminous red galaxies observed with the Southern African Large Telescope. Mon. Not. R. Astron. Soc. 467, 3239 (2017)
D. Stern, R. Jimenez, L. Verde, M. Kamionkowski, S.A. Stanford, Cosmic chronometers: constraining the equation of state of dark energy. I: H (z) measurements. J. Cosmol. Astropart. Phys. 2010(02), 008 (2010)
M. Moresco, Raising the bar: new constraints on the Hubble parameter with cosmic chronometers at \( z \sim 2\). Mon. Not. R. Astron. Soc. 450, L16 (2015)
T. Delubac, J.E. Bautista, J. Rich, D. Kirkby, S. Bailey, A. Font-Ribera, A. Slosar, K.-G. Lee, M.M. Pieri, J.-C. Hamilton et al., Baryon acoustic oscillations in the Ly\(\alpha \) forest of BOSS DR11 quasars. A &A 574, 59 (2015)
A. Font-Ribera, D. Kirkby, J. Miralda-Escudé, N.P. Ross, A. Slosar, J. Rich, É. Aubourg, S. Bailey, V. Bhardwaj, J. Bautista et al., Quasar-Lyman \(\alpha \) forest cross-correlation from BOSS DR11: baryon acoustic oscillations. J. Cosmol. Astropart. Phys. 05, 027 (2014)
Acknowledgements
M. Gonzalez-Espinoza acknowledges the financial support of FONDECYT de Postdoctorado, N\(^\circ \) 3230801. G. Otalora acknowledges Dirección de Investigación, Postgrado y Transferencia Tecnológica de la Universidad de Tarapacá for financial support through Proyecto UTA Mayor 4731-23. C. Rodriguez-Benites and M. Alva-Morales acknowledge the financial support of PE501082885-2023-PROCIENCIA.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Code availability
My manuscript has no associated code/software. [Author’s comment: The present work is a theoretical study, and no associated code/software is involved.]
Appendices
Appendix A: Hubble’s parameter data
In this appendix, we present Hubble’s parameter data for \(0.01< z < 2.360\):
Appendix B: Thermodynamics calculations
From Figs. 16, 17, and 18 we see that
where \(T_{de}^{(0)}=T_{de}(z=0)\) and \(T_{m}^{(0)}=T_{m}(z=0).\) Then, from (B1) we obtain
We know that for \(Q>0,\) the second law of thermodynamics requires \(T_{m}<T_{de}\) [87]. Therefore, we have:
Consequently, we get
Rights and permissions
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://creativecommons.org/licenses/by/4.0/.
Funded by SCOAP3.
About this article
Cite this article
Rodriguez-Benites, C., Gonzalez-Espinoza, M., Otalora, G. et al. Revisiting the dynamics of interacting vector-like dark energy. Eur. Phys. J. C 84, 276 (2024). https://doi.org/10.1140/epjc/s10052-024-12613-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-024-12613-z