A power-law coupled three-form dark energy model

We consider a field theory model of coupled dark energy which treats dark energy as a three-form field and dark matter as a spinor field. By assuming the effective mass of dark matter as a power-law function of the three-form field and neglecting the potential term of dark energy, we obtain three solutions of the autonomous system of evolution equations, including a de Sitter attractor, a tracking solution and an approximate solution. To understand the strength of the coupling, we confront the model with the latest Type Ia Supernova (SN Ia), Baryon Acoustic Oscillations (BAO) and Cosmic Microwave Backround (CMB) radiation observations, with the conclusion that the combination of these three databases marginalized over the present dark matter density parameter $\Omega_{m0}$ and the present three-form field $\kappa X_{0}$ gives stringent constraints on the coupling constant, $-0.017<\lambda<0.047$ ($2\sigma$ confidence level), by which we give out the model applicable parameter range.


I. INTRODUCTION
According to cosmological observations, the universe has entered a stage of an accelerated expansion with a redshift smaller than 1 [1,2]. Since all usual types of matter with positive pressure decelerate the expansion of the universe, a sector with negative pressure named as dark energy was suggested to account for the invisible fuel that accelerates the expansion rate of the current universe [3,4].
The simplest cosmological model of dark energy is the so called Lambda cold dark matter (ΛCDM) model, in which vacuum energy plays the role of dark energy. Although ΛCDM model provides an excellent fit to a wide range of astronomical data so far, such model in fact is theoretical problematic because of two cosmological constant problems, the fine tuning problem that why is the observational vacuum density so small compared to the theoretical one, and the coincidence problem that why is the observational vacuum density coincidentally comparable with the critical density at the present epoch in the long history of the universe. In order to alleviate the latter, various of evolving and spatially homogeneous scalar field, including quintessence [5], phantom [6], dilatonic [7], tachyon [8]and quintom [9] etc. were suggested to take vacuum energy's place. In these models, the resolution of the coincidence problem typically leads to a fine tuning of model parameters.
Since the experimental evidences of cosmology-specific scalars particles have not been discovered yet, there is no reason to exclude the possibility of some other high form field to be dark energy. Indeed, the three-form cosmology proposed in [10,11] could be a good alternative to scalar cosmology, because such high form field not only respects the Friedmann-Robertson-Walker(FRW) symmetry naturally but also can accelerates the expansion rate of the current universe without a slow-roll condition. Moreover, some interesting results, e.g.three-form with simple potentials lead to models of inflation with potentially large non-Gaussian signatures [12], etc., about three-form cosmology are obtained.
The coincidence problem mentioned above is just that the amount of dark matter is comparable to that of dark energy in the present universe, so it is natural to consider an interaction between these two components. As was pointed out in the paper [13], in comparison to coupled scalar dark energy model [14], some new features appear in coupled three-form dark energy model, including one that the stress tensor is modified by the interaction between two dark sectors, hence it is problematic to consider a coupled three-form dark energy model in a phenomenological way [15], and one needs to construct it in a Lagrangian formalism. Different from modeling dark matter as point particles [13], we follow the thread that describing the interaction between dark energy and dark matter from a fundamental field theory point of view [16] and consider dark matter as a Dirac spinor field.
The contents of this paper are as follows. In section II, we present a type of Lagrangians describing the interaction between a three-form field and a Dirac spinor field in curve spacetime and then derive the field equations from such Lagrangians. In section III, we consider these field equations in a FRW space-time by assuming the effective mass of dark matter as a power-law function of the three-form field and setting the potential of dark energy to be zero. In section IV, we carry out a simple likelihood analysis of the model with the use of 580 SN Ia data points from recently released Union2.1 compilation [17] and BAO data from the WiggleZ Survey [18], SDSS DR7 Galaxy sample [19]and 6dF Galaxy Survey datasets [20], together with CMB data from WMAP7 observations [21]. In the last section, we make a brief conclusion with this paper.

II. A TYPE OF FIELD THEORIES OF THREE-FORM AND DIRAC SPINOR IN CURVE SPACE-TIME
This section involves some concepts that are used to include fermionic sources in the Einstein theory of gravitation and for a more detailed analysis the reader is referred to [22][23][24][25].
A type of Lagrangians which describe the interaction between a canonical three-form field A αβγ with a potential V (A 2 ) and a Dirac spinor field ψ in a curve space-time can be constructed as where F = dA represents the field strength tensor and D µ is the covariant derivative of spinor which satisfies The Ω µ = 1 2 ω µab Σ ab appeared in (2),(3) denotes the spin connection which is constituted by the Ricci spin coefficients ω µab = e ν a ∇ µ e νb and the generators of the spinor representation of the Lorentz group Σ ab = 1 4 [γ a , γ b ] . γ a and Γ µ = e µ a γ a are the Dirac-Pauli matrices and their curve space-time counterparts respectively. Following the general covariance principle, the tetrad e ν a is related to the metric by g µν = e µ a e ν b η ab with η ab = diag(1, −1, −1, −1). In such Lagrangians, the coupling between two fields is demonstrated by the the function M(A 2 ) which is the effective mass of dark matter.
One now can obtains the field equations from the total action where L = L g + L m = R 2κ 2 + L m is the Lagrangian including gravity, R denotes the Ricci scalar and κ = √ 8πG is the inverse of the reduced Planck mass.
By varying the total action with respect to the three-form field and the Dirac field, we have the following equations of motion which are quite similar to that of electrodynamics The variation of the action with respect to the tetrad leads to Einstein field equation where the total energy-momentum tensor for two fields is given by In the end of this section, we show the dual description of the above three-form field theory to place this field theory in a context more familiar to most of the community, such dual description has the following action: where L = L g + L m and A represents the dual of the three-form. V ( A 2 ) and M ( A 2 ) represent the self-coupling of A and the coupling between two fields respectively. By varying the action with respect to the A, ψ andψ, we have the following equations of motion which are different from that of three-form model.

III. COSMOLOGICAL EVOLUTION OF THE POWER-LAW COUPLED THREE-FORM DARK ENERGY MODEL
We now consider the field equations in a homogeneous, isotropic, and spatially flat spacetime described by the metric where a(t) refers to the scale factor.
To be compatible with FRW symmetries, the three-form field is assumed as the time-like component of the dual vector field, i.e.
Since a three-form field without a potential can accelerates the expansion rate of the current universe 1 , for simplicity, we set the potential to be zero in the following discussions, together with choosing the coupling function as the following power-law form 2 we have the Friedmann equations 1 A three-form field without a potential is equivalent to a cosmological constant [11]. 2 For simplicity, we consider X ≥ 0 and neglect the absolute value sign in the following discussion. with λ is a dimensionless constant representing the strength of the coupling, this means that if λ = 0, such field theory becomes a free field theory, so the constant m with mass dimension is, in fact, the mass of dark matter in a free field theory. Since typically it is very difficult for dark energy to couple dark matter with mass bigger than milli-eV [26], we choose m to be smaller than milli-eV. Although the phenomenological bounds on dark matter mass coming from large scale structure require that most of dark matter is considerably heavier than 10 −3 eV [27], some dark matter's mass can be smaller than milli-eV. Indeed, Rajagopal, Turner and Wilczek considered axino in the keV range and they obtained the axino mass bound m a < 2 keV for axino to be warm dark matter [28].
In the FRW space-time, there is only one independent equation of motion of the threeform fieldẌ from the equations of motion of the spinor field and its Dirac adjoint, one can obtains the following equation with its simple solutionψ which shows that our model indeed returns to ΛCDM model when the coupling constant λ becomes 0.
Providing with the equations of motion, we have the continuity equations for both com- with the prime stands for derivative with respect to e-folding time N = ln a here and in the following.
In order to study cosmological dynamics in such coupled dark energy model, it is convenient to introduce the following dimensionless variable [29] x = κX, y = κ √ 6 (X ′ + 3X), By applying the Friedmann equations and equations of motion, one can obtains the autonomous system of evolution equations We note that ω 2 has been eliminated by the Friedmann constraint written in terms of these variables y 2 + ω 2 = 1.
There are two fixed points for such autonomous system. One of them is 2 3 , 1 , which is an attractor since its eigenvalues (−3, −3) are both negative. By rewriting density and pressure in term of the dimensionless variables, we have the total EOS indicating that such fixed point represents a three-form saturated de Sitter universe. The other one The trajectories with respect to x(N) and y(N) with a wide range of initial conditions and a assumption that λ = 0.01 (we will see that this is a good choice in the next section) can be visualized by Fig.1.
As is showed in the Fig.1, the trajectories run toward the de Sitter attractor, coasting along the saddle point. Also, one should note that the present value of x must be equal to or larger than a certain value x * that depends on both λ and y(0) to make sure in high redshift the autonomous system does not encounter singularity. After specifying λ and y(0), if x is i.e. x(N) > 0(−∞ < N < +∞), then x * is the lower limit of x. Now let us solve the the autonomous system of evolution equations by assuming a large x 0 . With the constraint y 2 < 1, it can be well approximated by two independent equations which have the following solutions we have replaced e-folding time by the redshift here.
Substituting the solutions into(27)-(29), we have noting that although δ is a constant, such solutions are different from the models proposed in [30][31][32] since is not a constant.
It can be inferred from (40) and (41) that the energy transfer between two dark sectors keeps the density of dark energy as a constant even when its EOS deviates from −1.
To compare the solutions with uncoupled dark energy model, let us rewrite Hubble parameter as is the normalized effective dark energy density and ω Xef f (z) = −1 + 1 3 is the effective EOS of dark energy.
By substituting the approximate solutions into (43) and (44), we have as one can see, depending on the sign of λ, there are two different high redshift approximate expressions for the effective energy density or the effective EOS. More specifically, effective energy density is approximate to Ω m0 1−Ω m0 (1 + z) 3(1+λ) (λ > 0) and − Ω m0 1−Ω m0 (1 + z) 3 (λ < 0) at high redshift, and the effective EOS of dark energy is approximate to λ(λ > 0) and At the end of this section, we place another restriction on the likelihood function of (x 0 , λ, Ω m0 ) with the aid of the approximate solution of the autonomous system of evolution equations, in fact we can see from (45) and (46) immediately that the Hubble parameter is independent of x 0 supposing x 0 take a large value, which means that the likelihood function becomes a none zero constant(in fact a large x 0 with proper values of λ and Ω m0 is quite favored by observations) with respect to x 0 if x 0 is large enough. Given this behavior of the likelihood function L, we have such formula which can be used to marginalize over x 0 without a prior.

IV. CONFRONT THE POWER-LAW COUPLED THREE-FORM DARK EN-ERGY MODEL WITH OBSERVATIONS
In this section, we perform a simple likelihood analysis on the free parameters of the model with the combination of data from Type Ia Supernova (SN Ia), Baryon Acoustic Oscillations (BAO) and Cosmic Microware Backround (CMB) radiation observations.
Firstly, we construct the following χ 2 function for SN Ia by using the recently released Union2.1 compilation with 580 data points where P , Q and R are defined as with µ th = 5 log 10 (1 + z) z 0 H 0 H(z ′ ) dz ′ +25 denotes the distance modulus predicted by theory and µ obs represents the observed one with a statistical uncertainty σ µ .
In the second step, we consider BAO data from the WiggleZ Survey, SDSS DR7 Galaxy sample and 6dF Galaxy Survey together with CMB data from WMAP 7 yeas observations to obtain the BAO/CMB constraints on the model parameters by defining χ 2 BAO/CM B as [33][34][35] 3 represent the co-moving angulardiameter distance and the dilation scale respectively, while z * ≈ 1091 is the decoupling time.
is the inverse of the correlation matrix.
Finally, the total χ 2 function for the combined observational datasets is given by χ 2 = χ 2 SN Ia + χ 2 BAO/CM A , from which we can construct the likelihood function as L = L 0 e − 1 2 χ 2 , here L 0 is a normalized constant which is independent of the free parameters.
Providing with the likelihood function, one can then obtain the best-fit values of the free parameters by maximizing it. However, as was mentioned in the previous section, the parameter x 0 can't be strictly restricted, so we leave it out of the discussions and consider the likelihood function that has been marginalized over x 0 without a prior. By the aid of (47), such function is proportional to wherex 0 is a large number, one can choose it as 10000, for example. We now present the fitting result in Fig.2 by analyzing such likelihood.
One may note that the marginalized likelihood in the right panel of Fig.2 appears a little non-Gaussian, this is mainly because of the non-Gaussian structure of the likelihood that haven't been marginalized. One also can see from the Fig.2 that observations favor a small positive coupling constant which, as we mentioned above, allows the existence of a tracking solution that can be used to alleviate the coincidence problem with a fine-turning of the model parameters.
Moreover, one thing here needs to be noticed, since x 0 can't be strictly restricted, the interact behavior between two dark sectors still remains uncertain, which can be inferred from Fig.3. However, one may decrease such uncertainty by taking into account observational constraints from future measurements.
From Fig.3, it can be inferred that the direction of energy transfer can be changed if x 0 is sufficiently close to x * , which is around 0.6 in such case, and the behavior of δ is almost the same between the choices of x 0 = 10 and x 0 = 100 if the redshift z > 0, which is consistent with the conclusion that the likelihood function becomes a constant with respect to x 0 if x 0 is adequately large that we have drawn in the section III. In fact we can prove this conclusion in a inductively way by plotting χ 2 (see Fig.4).

V. CONCLUSIONS
In this paper we have studied a power-law coupled dark energy model which considers dark energy as a three-form field and dark matter as a spinor field. By performing a dynamical analysis on the field equations with the introduction of three dimensionless variables, we obtained two fixed points of the autonomous system of evolution equations, among which one is a de Sitter attractor, and the other is a tracking solution, supposing λ > 0, that provides a possible solution of the coincidence problem.
By marginalizing over x 0 , we have also carried out a likelihood analysis on the free parameters λ and Ω m0 with the combination of SN Ia+BAO/CMB datasets, through which we have a best-fit value of the pair (λ, Ω m0 ) as (0.013, 0.269). In addition, the likelihood function marginalized over x 0 and Ω m0 showed that λ is restricted by −0.017 < λ < 0.047 (2σ confidence level, with a best-fit value 0.01), indicating that the measurements considered here are quite consistent between ΛCDM and our three-form model. However, future measurements might allow us to tell them apart.
Notwithstanding it can be told from the fitting result that λ and Ω m0 are strictly restricted, x 0 can be any value beyond x * . However, as mentioned above, future measurements might decrease the uncertainty on x 0 .