Revisiting the dynamics of interacting vector-like dark energy

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.


I. 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].
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 Ref. [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 timedependent length.This model, also known as the vectorlike 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 ∼ 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 illu-minate 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 vectorlike 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 section II, we introduce the IVDE models, present the total action of the model, and derive the field equations in a general spacetime.We then obtain the cosmological equations and define the effective dark energy, along with its coupling to dark matter.In Section III 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 Section IV, we corroborate our analytical findings by numerically solving the field equations.In Section V, 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 Section VI, we summarize the results obtained.

II. 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 κ 2 = 8πG, with G being the gravitational constant, F a µν = ∂ µ A a ν − ∂ ν A a µ and A a2 = g µν A a µ A a ν .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 µ we find the equation of motion Below, we detail the basic equations for the cosmic triad in a cosmological background.

A. 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 A 2 represents the scalar potential, H ≡ ȧ/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 ρ 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 /ρ m = 0 and w r = p r /ρ 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 ρ + 3H(ρ + p) = 0, (13) 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 ρ m ∼ a −3 and ρ r ∼ a −4 , respectively.
Additionally, we can define the effective dark energy EoS as It is also convenient to introduce the total equation-ofstate parameter as which is related to the deceleration parameter q through q = 1 2 (1 + 3w tot ), (16) 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].

B. 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] ρde + 3H(ρ de + p de ) = −Q, 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 ̸ = 0 and a new term due to Q will appear in its r.h.s [68].
In the following section, we conduct a thorough phasespace 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).

III. PHASE-SPACE ANALYSIS
We introduce the following set of dimensionless variables [2]: and the constraint equation: Therefore, we obtain the dynamical system where N ≡ 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 /ρ 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 ) ∼ e −κ 2 λ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 signchangeable [78].

A. Critical points
In this section, we obtain the critical points from the conditions dx/dN = dy/dN = dϱ/dN = du/dN = 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 ≥ 0, ϱ c ≥ 0 and u c ≥ 0.
Consequently, we find the critical points for each interaction term: 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 I and the values of their cosmological parameters in Table II.Also, in this subsection and so on, we define Ω (r) de and Ω (m) de 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 Ω (r) de = u 2 c .When u c = 0, we obtain the radiation-dominated solution with Ω r = 1 and w de = w tot = 1/3.For u c ̸ = 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 Ω (r) de < 0.045 [79,80], we need to have u c < 0.122.
The critical point c represents a dark energydominated 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 α = 0, the critical point labeled as b M represents a matter-dominated era, with Ω 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 α ̸ = 0, we have a scaling matter era which is constrained to satisfy Ω (m) de < 0.02 (95% C.L.), at redshift de z ≈ 50, according to CMB measurements [81].Thus, this leads us to the constraint 0 < α < 0.01, which is also compatible with 0 < Ω m < 1.
With the phase space variables mentioned earlier, we express the interaction term in the following form The critical points of the system described in Equation ( 24) under the interaction case given by Equation ( 31) are presented in Table III, along with the corresponding values of their cosmological parameters in Table IV.
The critical point labeled as d R corresponds to a scaling radiation era, characterized by Ω (r) de = u 2 c .When u c = 0, we obtain the radiation-dominated solution with Ω r = 1 and w de = w tot = 1/3.For u c ̸ = 0, this point describes a scaling radiation era.To satisfy the early constraint imposed by the physics of big bang nucleosynthesis (BBN) and ensure Ω (r) de < 0.045 [79,80], the condition u c < 0.122 must be met.
On another note, for β = 0, the critical point denoted as e M represents a matter-dominated era, where Ω 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 β ̸ = 0, we have a scaling matter era.Applying the constraint Ω (m) de < 0.02 (95% C.L.), at redshift de z ≈ 50, according to CMB measurements [81], we obtain 0 < β < 0.01, which is also compatible with 0 < Ω 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 < Ω m < 1, meaning that, 0 < −1/β < 1 or β < −1, but β must be positive and small to obtain a matterdominated 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.Where Cosmological parameters for the critical points in Table III.
Name Ω de Ωm Ωr 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 V and the values of their cosmological parameters in Table VI The critical point denoted as g R corresponds to a scaling radiation era, where Ω (r) de = u 2 c .When u c = 0, it gives rise to the radiation-dominated solution with Ω r = 1 and w de = w tot = 1/3.The case u c ̸ = 0, corresponds to a scaling radiation era.Therefore, to meet the early constraint imposed by the physics of big bang nucleosynthesis (BBN) and ensure Ω (r) de < 0.045 [79,80], the condition u c < 0.122 must be satisfied.
Conversely, when η = 0, the critical point labeled as h M represents a matter-dominated era, characterized by Ω 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 η ̸ = 0, this point represents a scaling matter era.Thus, from constraint Ω (m) de < 0.02 (95% C.L.), at redshift de z ≈ 50, according to CMB measurements [81], we get −0.04 < η < 0, which is also compatible with 0 < Ω 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 < Ω m < 0 implies η < −3 or η > 6, but this result is not compatible with a dark matter-dominated era as described by point h M (−0.04 < η < 0).Hence, to reproduce the standard thermal history of the universe, we designate the points i 1 and i 2 as unphysical.

B. Stability of critical points
To study the stability of the critical points, we introduce time-dependent linear perturbations denoted as δx, δy, δϱ, and δu around each critical point.These perturbations take the form of x = x c + δx, y = y c + δy, ϱ = ϱ c + δϱ, and u = u c + δu.Substituting these expressions into the autonomous system (24) and linearizing the equations, we obtain the linear perturbation matrix M, as outlined in [2].The eigenvalues of M, denoted as µ 1 , µ 2 , µ 3 , and µ 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 character- unstable for η < 3 5 and η > 2 ized 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 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.
1. Case I: Q = 3αHρm • Point a R has the eigenvalues then this point is always unstable for all the values of α.
• Point b M has the eigenvalues which is a saddle point when α > 1 3 or −1 < α < 1 3 .On the other hand, it becomes a stable node when α < −1.However, this point cannot account for the current accelerated expansion of the Universe.
• Point c has the eigenvalues This is a de-Sitter solution, ensuring accelerated expansion for all parameter values.We observe that It is shown that both trajectories converge to the attractor c, which is a stable node that describes the dark-energy-dominated universe.Also, we have matched the current values for the fractional energy densities of dark energy Ω (0) de ≈ 0.68 and dark matter Ω (0) dm ≈ 0.32 , at redshift z = 0, according to Plancks results [82] it is a stable node when α > −1 and 1 < λ ≤ 16  7 .Finally, this point never exhibits stable spiral behavior.

Case II
which tells us that it is always an unstable node.
• Point e M has the eigenvalues which is a saddle point when β < − 1 3 or β > − 1 4 .On the other hand, it becomes a stable node when − 1 3 < β < − 1 4 .However, this point cannot account for the current accelerated expansion of the Universe.
• Point f has the eigenvalues This is a de-Sitter solution, hence ensuring accelerated expansion for all parameter values.We observe that it is a stable node when 1 < λ ≤ 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.3. Case III: • Point g R has the eigenvalues which is always an unstable node.
• Point h M has the eigenvalues which acts as a saddle point when 2 3 < η < 1, η > 1, or η < 3  5 .Conversely, it behaves as a stable node when 3  5 < η < 2 3 .Nonetheless, this point cannot explain the current accelerated expansion of the Universe.
• Point i has the eigenvalues This is a de-Sitter solution, ensuring accelerated expansion for all parameter values.We observe that it is a stable node when η > −3 and 1 < λ ≤ 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.
A summary of this section is shown in Table VII, which contains the principal properties of critical points.

IV. NUMERICAL RESULTS
In this section, we aim to numerically solve the autonomous system represented by equations ( 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.

A. Case I: Q = 3αHρ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 Λ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 Λ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 ΛCDM.In the same plot, we have depicted the relative difference for the ΛCDM model, showing that our results are compatible with observations.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 λ.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 ΛCDM.
Interestingly enough, the evolution of the total EoS of our model is very close to that of Λ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 We depict the total EoS parameter wtot (orange line), the EOS parameter of dark energy w de (blue line), and the total EOS of the ΛCDM model (green line) as redshift functions.We also used the same initial conditions as shown in Figure 8 to obtain both the solid and dashed blue lines.For the dashed line, a value of w de ≈ −0.968 is observed at the current time z = 0.Meanwhile, for the solid line, a value of w de ≈ −0.999 is observed at the current time z = 0, which is consistent with the observational constraint w (0) de = −1.028± 0.032.
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 Λ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].
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 Λ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].

V. 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 = nV = 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, ρ), and using the conservation laws where ρ 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)ρ we obtain For two interacting fluids (dark matter and dark energy) one can write ρde + 3H(1 + w ef f de )ρ de = 0, (51) 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 ef f 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 (0) de and T m (z)/T (0) m , respectively, as functions of redshift z for different initial conditions.Here, we define T (0) m ≡ T m (z = 0) and T (0) de ≡ 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 = −QV dt and T m dS m = QV dt, 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 (0) de < T m (z)/T (0) m .Thus, for Q > 0, the validity of the second law requires T m < T de [87], and therefore 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]).

VI. 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", FIG.18.We depict the evolution of the temperature of dark energy and dark matter as a function of the redshift for interaction "Case III" using the same set of initial conditions used in Figure 12. 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 energydominated 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 phasespace 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) 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].

FIG. 6 .FIG. 7 .
FIG.6.We show the evolution of the deceleration parameter q(z), for the same initial conditions used in figure4.

FIG. 10 .FIG. 11 .
FIG.10.We show the evolution of the deceleration parameter q(z), for the same initial conditions used in figure8.

FIG. 15 .
FIG.15.We present the evolution of the Hubble rate H(z) and its relative difference ∆rH(z) = 100 × |H − HΛCDM|/HΛCDM with respect to the ΛCDM model, as functions of redshift, using the same initial conditions as shown in Figure12.This is complemented by the evolution of the Hubble rate HΛCDM in the ΛCDM model and the Hubble data from references[85,86].We have employed the current value of the Hubble rate, H0 = 67.4km/(Mpc•s), from Planck 2018[82].

FIG. 16 .FIG. 17 .
FIG.16.We depict the evolution of the temperature of dark energy and dark matter as a function of the redshift for interaction "Case I" using the same set of initial conditions used in Figure4.

FIG. 19 .
FIG.19.We depict the evolution of the interaction term Q for each case considered.

TABLE I .
Critical points for the autonomous system.

TABLE II .
Cosmological parameters for the critical points in TableI.

TABLE III .
Critical points for the autonomous system.

TABLE V .
Critical points for the autonomous system.Where, Rc =

TABLE VI .
Cosmological parameters for the critical points in TableV.

TABLE VII .
Properties of the critical points