An unified cosmological evolution driven by a mass dimension one fermionic field

An unified cosmological model for an Universe filled with a mass dimension one (MDO) fermionic field plus the standard matter fields is considered. After a primordial quantum fluctuation the field slowly rolls down to the bottom of a symmetry breaking potential, driving the Universe to an inflationary regime that increases the scale factor for about 71 e-folds. After the end of inflation, the field starts to oscillate and can transfer its energy to the standard model particles through a reheating mechanism. Such a process is briefly discussed in terms of the admissible couplings of the MDO field with the electromagnetic and Higgs fields. We show that even if the field loses all its kinetic energy during reheating, it can evolve as dark matter due a gravitational coupling (of spinorial origin) with baryonic matter. Since the field acquires a constant value at the bottom of the potential, a non-null, although tiny, mass term acts as a dark energy component nowadays. Therefore, we conclude that MDO fermionic field is a good candidate to drive the whole evolution of the Universe, in such a way that the inflationary field, dark matter and dark energy are described by different manifestations of a single field.


Introduction
A new class of mass dimension one fermionic fields has been proposed by Ahluwalia and Grumiller [1][2][3], which are constructed by means of charge conjugation spinors. In its first formulation the spin sums of the quantized fields showed to be Lorentz violating due to the intrinsic definition of their dual fields, which rose several doubts on the applicaa e-mail: s.pereira@unesp.br b e-mail: marcio.alves@unesp.br c e-mail: thiago.mogui@gmail.com bility of this new field. Such an old construction were named Elko 1 fields. However, recently [4,5] a subtle deformation in the dual structure solved the problems of Lorentz violation, putting the theory in solid bases from the point of view of the quantum field theory, and the new fields are called just mass dimension one (MDO) fermionic fields. In this work we are interested just in the classical analogue of MDO fermionic fields, thus we can use Elko or MDO fields without loss of generality.
The fermionic fields constructed from charge conjugation spinors are natural candidates to describe dark matter particles in the universe, since they are neutral and have canonical mass dimension one, which make them to couple very weakly to other particles of the standard model. Indeed, the only admissible couplings are with scalars and Higgs fields [6][7][8] and with electromagnetic stress tensor [9], in addition to quartic self-couplings. Different from Dirac fermion, which has mass dimension 3 2 , is parity conjugation invariant and admits several couplings with standard model particles, the MDO field born with an intrinsic dark character. Consequently, detection and constraints in its physical properties is difficult, once it does not couple directly with standard model particles.
Several cosmological applications of this new field have been done, first in torsion free frameworks [10][11][12][13][14][15][16][17][18][19] and more recently considering its coupling to torsion [20][21][22][23] in an Einstein-Cartan framework, which turns the system of equations richer than the previous case, especially for inflationary applications. In particular, in Refs. [22,23] the numerical results describing inflation and dark matter evolution were presented together with the energy density evolution, in good agreement with the required inflationary number of e-folds and also with the expected energy density for the dark matter component. Notwithstanding, we refer the reader to other important recent results regarding the MDO fermionic field in the context of the quantum field theory [24][25][26][27], and also in thermal field theory [28]. The study of scalar and tensor perturbations for Elko field has been done in [13,14] in the torsion free case and the study of first order vector perturbations has been done in [15].
In the present article we consider the fermionic MDO field as a candidate to drive the complete evolution of the Universe. A symmetry breaking potential is responsible for the inflationary phase, while a reheating like mechanism drives the field to the radiation phase. The effective mass of the field plays an essential role in the energy transfer process. After the field rolls down to the true vacuum of the potential, the system follows a dark matter evolution due to a natural (gravitational) coupling between the MDO field and the baryonic energy density. Finally, a dark energy accelerated phase can be obtained from the constant quadratic mass term, which acts exactly like a cosmological constant term. All different phases are connected and the free parameters can be constrained by observational data and by some reasonably physical assumptions. The cosmic coincidence problem can also be naturally understood in such a scenario.
It is worth to stress that the search for an unified cosmological model is an old task and most of the models are based on a single scalar field [29][30][31]. The advantage to use scalar fields is that the inflationary phase and the generation of primordial perturbations are in good agreement with observations [32][33][34]. After reheating, if the decay of the scalar field is incomplete, it may act as dark matter, while its zero point energy acts as dark energy. Although this is a standard model framework, scalar fields does not furnish any physical interpretation for inflation and also no scalar field has been observed yet, except the Higgs boson.
On the other hand, the general idea of considering the usual Dirac fermions to explain inflation, dark matter and dark energy have already been considered in the literature (see, e.g., [35][36][37][38]). In this spirit, the aim of the present article is to show that the MDO fermionic field is a more reasonable field to drive inflation and the subsequent stages of evolution of the Universe in an unified fashion. Being neutral and not interacting with other particles of the standard model, it is a good candidate to describe dark matter without the need of supposing an incomplete decay during reheating, an it happens in the case of scalar dark matter particles.
Furthermore, as pointed out by Pereira et al. [22,23], the "MDO inflation" can be interpreted in light of the Pauli exclusion principle. When the fermionic field rolls down to the minimum energy state of the potential, the Pauli exclusion principle starts to act, not allowing all particles to occupy the lowest energy state. If the potential is stronger than the degeneracy pressure, the whole system responds with an abrupt expansion in order to accommodate all particles in the lowest energy state, since the spacing between energy levels in a bound system is inversely proportional to the size of the system. This effect allows all the particles to accommodate very close to the lower energy state after inflation.
Moreover, since the MDO particles satisfy a Klein-Gordon like equation, the set of equations describing the evolution of the Universe is very similar to the scalar field case. Thus, several features of the scalar field cosmology are recovered, but for a fermionic field instead.
The paper is organized as follows. Section 2 presents the basic Friedmann-Lemaître-Robertson-Walker (FLRW) equations already derived in [21][22][23] in the presence of torsion terms. The dark matter and the dark energy behavior of the MDO field are presented in the Sect. 3, where some parameters are constrained with observational data. In the Sect. 4 the MDO inflation is studied by means of a numerical solution of the evolution of the field subject to a symmetry breaking potential. Two possible ways of implement the reheating phase in the MDO cosmological scenario is briefly discussed in the Sect. 5, and in the Sect. 6 we finish with our conclusions.

Dynamic equations for the MDO fermionic field
The action for the model is [20,21] where κ 2 ≡ 8π G = 8π/m 2 pl with c =h = 1 and S m is the usual action for other matter fields, as baryonic matter or radiation. The tilde represents torsion terms into the Ricci scalar R and covariant derivatives, namely, ΛΓμ, whereΓ μ are the spin connections. In the Einstein-Cartan framework the contorsion terms generalizes the affine connectionΓ ρ μν = Γ ρ μν + K ρ μν and are given by where the only non-vanishing torsion terms obeying cosmological principle are [39] where the torsion functions h(t) and f (t) must be determined.
Considering the flat FLRW metric ds 2 = N (t) 2 dt 2 − a(t) 2 [dx 2 + dy 2 + dz 2 ], where N (t) is the lapse function, the two Friedmann equations, the dynamic field equation for φ(t) and the torsion functions h(t) and f (t) can be obtained by taking the variation of the Lagrangian of the model with respect to N (t), a(t), φ(t), h(t) and f (t) respectively (see [20][21][22][23] for further details). Thus we obtain (setting N → 1 at the end) h(t) = − 1 8 where H =ȧ/a, as usual, and ρ m and p m are the energy density and the pressure of other matter components, which satisfies a conservation equation for a perfect fluiḋ On the other hand, the energy density and the pressure of the MDO field is given by The set of Eqs. (5-11) looks-like a generalization of the standard scalar field model, but it is important to keep in mind that here φ(t) is just the temporal part of the fermionic field Λ. Now, by substituting the Eqs. (5) and (6) in the Eqs. (10) and (11), it is possible to write the above expressions of ρ φ and p φ in a different useful form which can be straightforwardly combined to show thaṫ Notice also the presence of a coupling between the MDO field with the energy density and pressure of the standard matter into the first terms of Eqs. (12) and (13). Such terms are a manifestation of the coupling between the spin components of the MDO field with gravity, such that if one takes the Minkowski spacetime (H = 0) then ρ φ and p φ for the scalar field are exactly recovered, as can be directly checked by means of the Eqs. (10) and (11).
From the above formulation, we see that the present model is equivalent in form to a model with zero torsion containing only perfect fluid components, such that one of them has the energy density and pressure defined by the Eqs. (12) and (13). Moreover, the Eq. (14) is obviously equivalent to the Eq. (7), and the Eqs. (5) and (6) can now be written as: Regarding the potential V (φ), two forms were previously considered in the literature. The first one is a symmetry breaking potential [22] where A and φ c are positive constants. It was showed that as the field rolls down to the true vacuum of the potential at φ = φ c the inflation occurs with the correct number of e-folds, depending on the initial value of φ and also on the constants A and φ c . After inflation, a dark matter evolution follows naturally, leading to correct energy densities for different phases.
On the other hand, in the Refs. [21,23] the potential was chosen to be of the form: where m is the physical mass of the field and α is a dimensionless coupling constant. Considering such a potential, the inflationary and dark matter evolution were also obtained [23], with the correct numerical estimates to the energy density at different epochs. Furthermore, a scenario where φ(t) is a slowly varying function at late times was interpreted as a time varying cosmological term proportional to H 2 [21]. Despite the successful to describe some individual phases of evolution with correct numerical estimates, the radiation phase and a smooth transition to late cosmic acceleration were not completely addressed by the model. In order to describe all the phases of evolution of the Universe in a con-sistent unified fashion, here we will consider a potential of the form with V 1 and V 2 given by (17) and (18), and also the presence of the standard matter fields. As it will be seen in the following sections, all the stages of evolution of the Universe can be recovered in a natural way. In order to better understand the role of the above potential in the dynamic of the field let us write it in the expanded form as where c is an effective self coupling constant and C = A 4 a constant. The potential is attractive whatever the φ value ifᾱ > 0, which leads to α > −4 A 4 /φ 4 c , and the minimum occurs at For m 2 << 4 A 4 /φ 2 c and α << 4 A 4 /φ 4 c the minimum occurs at φ min ≈ φ c , which will be the case when we fix the parameters with observational constraints.

Dark matter and dark energy evolution
From now on we will consider that the only additional standard matter present is of pressureless baryonic type, namely ρ m = ρ b and p m = p b = 0. Let us start with the late time evolution of the Universe in order to constrain the parameters φ c , m and α with observational data. The value of φ c is of particular interest in the inflationary epoch for which we give a detailed description in the next section. For now, it is enough to suppose that, at early times, the field φ is initially at rest around the false vacuum of V 1 and that after a quantum fluctuation, it rolls down to the bottom of the potential acquiring, after a period long enough, the constant value φ c . Hence, when the kinetic energy of the MDO field is finally negligible, with the field satisfyingφ 0, we have V 2 >> V 1 . Therefore, at late times when only the baryonic matter and the MDO field are relevant for the cosmic dynamics, we can use the Eqs. (12) and (13) to obtain and where ρ b is the baryonic energy density and the potential has the fixed value V 2 (φ c ) = 1 2 m 2 φ 2 c + 1 4 αφ 4 c . From the above expressions, it is clear that at late times the energy density of the MDO field is an addition of two distinct contributions. The first one behaves as a pressureless fluid following the evolution of the baryonic matter, and the second one is an effective cosmological constant given by Inserting ρ φ and ρ b into the Friedmann equation (15), the Hubble function can be obtained as follows where 2 Ω b = ρ b,0 /ρ crit,0 and now it is clear that the second term in the right-hand-side of the above equation can be interpreted as a dark matter component with the present density parameter defined as which comes from the gravitational coupling of the MDO field with standard baryonic matter. The last term act as a cosmological constant term, with Ω Λ,φ ≡ Λ eff /κ 2 ρ crit,0 . Notice that Eq. (25) has exactly the same form as the ΛCDM model, with baryonic and dark matter components following the same evolution as a −3 . At the end of the inflationary and reheating stages, gravity is the only non-negligible interaction of the fermionic MDO field (besides its self-interaction) which can cluster, enabling the posterior formation of structures in the Universe. The above definition of the dark matter density parameter of MDO shows a remarkable feature of this field, namely, its natural coupling with other matter fields which comes directly from the structure of the equations of motion in curved spaces. It is worth to stress that we have not included any direct coupling in the action (1). Due to its gravitational origin, such a coupling disappears in the Minkowski space.
The usual standard scalar field can also behaves as dust matter, and unified inflation-dark matter scenarios have been explored in the literature (see, e.g., [29,31] and references therein). But in this case, the zero pressure is obtained by a time average of the field while it is coherently oscillating after the end of the slow-roll regime. Nevertheless, a dark matter contribution can only survive if the decay of the inflaton is incomplete. The MDO field also oscillates at the end of inflation as we will see in the Sect. 5, however it can evolve as dark matter at late times even if its kinetic energy decays completely during reheating. This is explained by the presence of the term evolving as a −3 in the Eq. (25) coming from the MDO energy density. Therefore, the scenario described in this Section would not be possible if the field goes to zero after inflation as, e.g., in the chaotic MDO inflation studied in the Ref. [23], since the coupling of the MDO field with the baryonic energy density is cancelled [see Eq. (22)]. The dark matter behaviour in this case is obtained by averaging the coherent oscillations of the MDO field in the present epoch leading to a zero pressure, in a similar manner to what happens to the scalar field inflationary model.
Notwithstanding, assuming that the symmetry breaking potential V 1 (φ) was dominant in the inflationary regime, the MDO potential acquires its true vacuum value as the amplitude of oscillations becomes virtually zero. Although small, such zero point energy leads to a non-negligible contribution today, which results in a non-null cosmological constant like parameter Ω Λ,φ .
Finally, in order to estimate the values for the parameters of the theory, let us now consider, for instance, the last results from the Planck collaboration on the cosmological parameters [40]. Considering a flat spatial section and H 0 = 67.4 km s −1 Mpc −1 we have Ω b 0.05, Ω DM,φ 0. 27 and With these values in the relation (26) we obtain φ c 1.3 m pl , and the value of the potential at its minimum can be obtained from (24) and (26) in terms of the cosmological parameters as From (24) we obtain directly the value of the effective cosmological constant, namely Λ e f f = 4.37 × 10 −66 eV 2 , almost exactly the value obtained from Planck collaboration 3 [40]. From (27), (18) and the values of α andᾱ from (20) and (21) we can also constraint the value for the physical mass of the field. With A ∼ 10 14 GeV (see next section) and φ c 1.3 m pl , the α parameter must be in the interval −10 −18 α 10 −125 . Thus, the range for the mass is 0 < m 2.4 × 10 10 GeV. It is important to remark that, contrary to scalar field dark matter models [29], here the contribution to dark matter comes just from the value of the incomplete field decay, namely φ c , through (26), and not from the mass m of the field. The mass and α parameter are related to the cosmological constant like term, Ω Λ,φ , which represents the dark energy sector. This opens the possibility to mass be in a very large range, which must be constrained by other methods, as the growth of primordial perturbations. 3 From [40], Λ = (4.24 ± 0.11) × 10 −66 eV 2 .

MDO inflation with a symmetry breaking potential
Now, we can ask ourselves if the energy density of the MDO fermionic field can be dominant at the very early stages of evolution of the Universe. If we remember that one of the admissible couplings of the MDO field Λ is with the Higgs field Φ, we may give an affirmative answer as follows. A scenario analogous to the hybrid inflationary model [41] can be pictured out for the present case. An interaction between Λ and Φ makes the latter goes to zero rapidly, while the Λ energy density could remain large for a much longer time since Λ does not interact appreciably with other fields. In the present article, nevertheless, we do not consider the details of such a coupling, but we just consider the stage of inflation at large energy density for the MDO field while Φ = 0, as it has already been discussed in the literature by one of us and other authors [22,23]. The coupling with the Higgs field could be particularly important in the last stage of inflation, and we intend to take this into account in forthcoming investigations.
Therefore, let us analyze the inflationary epoch by assuming that the energy density of the Universe is initially dominated by the MDO field. At the high energy regime of inflation the only relevant potential is the symmetry breaking potential V 1 (φ). At the time that inflation begins, t i , the initial value of the field φ i φ c , and the constant A can be chosen as A (3H 2 i m 2 pl /8π) 1/4 ∼ 5.3 × 10 14 GeV, such that the Hubble parameter assumes the following value H (t i ) ≡ H i ∼ 10 35 s −1 [34] at the GUT scale. Conversely, it is remarkable that the constant φ c is not fixed by inflation, but by the ratio between the amount of dark matter and baryonic matter in the Universe as it was shown in the last section [Eq. (26)].
The MDO field is initially in a nearly equilibrium point of the false vacuum of V 1 (φ). Since the field has mass dimension one (or energy dimension), we can use a quantum uncertainty relation in order to establish its initial value. A quantum fluctuation Δφ in a time interval Δt must satisfy ΔφΔt ≥ 1. Since inflation starts at t i ∼ 10 −35 s, we can choose Δt of this order to arrive at Δφ ∼ 1/t i ∼ 5.4 × 10 −9 m pl , which will be assumed as the initial condition for the field in what follows.
On the other hand, the initial value of the first derivative of the fieldφ(t i ) ≡φ i can be obtained by taking the equations of motion (7) at the time t i . Let us consider that the condition φ Hφ is satisfied when inflation starts, then we are lead tȯ If one takes φ c = 1.3 m pl , for instance, the initial value iṡ φ i = 1.86 × 10 26 m pl s −1 .
Having settled the initial conditions, one can numerically solve the system of differential equations (5) to study the evolution of the dynamical parameters of the system. In this Section we consider just the MDO field as a source of energy, not including baryonic matter or radiation, which are assumed to be not dominant in that phase. In the Fig. 1, the evolution of the field φ is shown as a function of the cosmic time t. From this figure it is clear that the field evolves from φ i to φ c while rolling down to the bottom of the potential, with small oscillations at the end of its evolution, which will be discussed in the next Section.
In order to study the duration and the kinematics of the inflation, it is useful to analyze the evolution of the slow roll parameters governing the inflationary epoch. They are defined in the usual way as Moreover, it is convenient to analyze the evolution of the relevant quantities as functions of the number of e-folds achieved at the time t, defined as follows In the Fig. 2a, the evolution of and η are shown. It has been also included the evolution of the Hubble parameter along with the torsion function and the field φ itself. Notice that assuming the above specified (well motivated) initial conditions, at the end of inflation ( = 1) we have the total e-folds N 71. But, from the uncertainty relation, the initial value of the field can be higher yielding a smaller number of e-folds during inflation. In this case, the field spends less time in the slow-roll regime which leads, therefore, to a decrease of N .
Although we have imposed |η| 1 initially, this parameter grows very rapidly in the beginning of inflation and assumes a constant value (η 0.3) during almost the entire inflationary period. In the remaining 10 e-folds until the end of inflation |η| is nearly zero, and finally it becomes large in the last e-fold of inflation, corresponding to a substantial variation of , which indicates the end of inflation.
Notice that the change in the parameter η is associated with the increase of the MDO field in the last 15 e-folds and also with a "bump" in the Hubble parameter in the same period. Such a bump, corresponding to a moderated increase in the energy density of the MDO field with respect to its initial value, is notably correlated with an increase of the torsion function h(t). From the Eq. (8) we see that initially, |h(t i )| is much smaller than the Hubble parameter, but it becomes of the same order of H as φ → φ c . In the Fig. 2b it is also possible to note that the height of the bump and the corresponding transient increasing of |h(t)| are strongly dependent on the value of φ c . The bump can even disappears for a value of φ c small enough, but in this case changes in the initial conditions would be required in order to achieve the desired total e-folds of inflation. An interesting question about the bump at the end of inflation is whether it could leave an imprint in the cosmic microwave background (CMB) radiation, since it is characterized by a sudden increase in the energy density at that time. A better quantitative treatment must be done in The field remains oscillating during all the reheating stage. The horizontal dotted line is the one period average of the deceleration parameter q(t) = 1/2 indicating that the MDO field has a matter behavior in the reheating epoch as it happens for a typical massive scalar field for instance order to study the constraints in CMB spectrum due to such bump.
We have also studied the evolution of the deceleration parameter at the inflationary period. In the left panel of the Fig. 3 one can notice that q(t) < −1 for a short period preceding the maximum of the bump, for which we have again q max = −1 since (Ḣ ) max = 0. Therefore, the present inflationary scenario for the MDO field exhibits a phase with a phantom like behavior. Again, it is explained by the increase of the function |h(t)| which contributes significantly for the increase of the effective energy density at that epoch. After the end of inflation, the torsion remains important until the present time and its contribution scales with the Hubble function H (t), since the field acquires a constant value of the order of the Planck mass. Such a behavior is expected, since we are assuming the whole Universe homogeneously filled with the MDO field, and, therefore, each point of space containing a fermionic field interacts with torsion, contributing to the energy-momentum tensor and providing a contribution to the evolution. On the other hand, in the case of the chaotic inflation with a potential described by a sum of a quadratic and a quartic self-interacting term, the field goes to zero after inflation, thus the torsion which was initially important vanishes after inflation [23]. In this case, there is no bump at the final stage of inflation.

Reheating
After the end of inflation ( = 1 and q = 0) the amplitude of the MDO fermionic field φ(t) starts to oscillate coher-ently around the minimum of the potential with a frequency ω φ = √ V (φ) as can be seen in the Fig. 1 and a reheating mechanism takes place. During this period, the deceleration parameter oscillates around its average value q(t) = 1/2 as shown in the right panel of the Fig. 3. Therefore, during the regime of coherent oscillations, the field behaves in average as non-relativistic particles, i.e. p φ = 0, in the same way it happens for a typical scalar field subject to a quadratic potential. Such a result has already been verified in the Ref [23] for the MDO inflation with a quadratic mass term plus a self-interacting potential, and also for a symmetry breaking potential as can be seen in the Ref. [22]. Now, let us briefly discuss two possible mechanisms for a successful reheating process after the MDO inflation, namely: (i) the electromagnetic coupling of the MDO field at high energies, and (ii) the coupling of the MDO field with the Higgs boson in a preheating phase preceding the reheating. Now, if the reheating process occurs with a MDO decay rate Γ we need to add a damping term Γφ in the left-handside of the Eq. (7) which is negligible in the inflationary epoch. Therefore, the energy density of the MDO field is not conserved, but there is a term −(1+κ 2 φ 2 /8)Γφ 2 in the righthand-side of the Eq. (14). From the Eqs. (12) and (13) one can show that averaging in one period of oscillation (1 + κ 2 φ 2 /8)φ 2 = ρ φ . Hence, by considering the mechanism (i), the field decays in photons and the evolution of the MDO field and the decay product's energy density are described by the coupled equationṡ where ρ r is the energy density of radiation. Notice that the above equations are identical in form to the case of reheating with a scalar field. As a result while the MDO field dominates the energy density, the deceleration parameter is 1/2 in average and the Universe is decelerating during all the reheating stage. Then, there is a transition to a radiation dominated epoch until the deceleration parameter reaches q rad = 1, which indicates the end of the reheating stage.
In the Ref. [9], the interaction of the MDO field with photons at the tree level was studied. It was shown that electromagnetic interactions are excluded at the 95 % CL by the LHC data for MDO masses up to 1 TeV. However, such a coupling might be relevant at the high energy levels achieved at the reheating phase if we use the superior limit of mass discussed earlier, m 2.4 × 10 10 GeV. Let us suppose that the MDO decay rate is approximately given by Γ g e m, where g e is a dimensionless coupling constant. In order to the damping term does not modify the evolution of the Universe in the inflationary epoch one has Γ H i , and than g e 1.8. If one considers the LHC upper bound g e 10 −5 , for example, we have Γ 2.2 × 10 5 GeV. Therefore, in this reheating scenario, the Universe entropy can be generated by transferring the energy of the MDO field to photons. Then, the high energetic photons can generate other particles via pair production process.
The second possible mechanism for energy transfer after inflation is justified by the coupling of the MDO field with the Higgs field through the addition of the interaction term in the Lagrangian density of the action (1). Such an interaction has been considered to investigate the possibility of discovering the MDO field in the LHC [6][7][8]. Now, the full Lagrangian is given by where the last term is the Higgs Lagrangian. At the time of reheating, the MDO field oscillates as with a slowly varying amplitude φ 0 (t) which is initially of the order of the Planck mass. Then, if we vary the above Lagrangian with respect to Φ, one can obtain the equations of motion of the Higgs field coupled with φ(t), which will be similar to a Mathieu equation. The solutions of this equation are characterized by a parametric resonance which results in a fast amplification of Φ. After such a period (called preheating) we have the reheating since the coupling of the Higgs field Φ with the standard-model particles leads to thermalization. Obviously, the temperature at the end of reheating depends on the details of the model. In practice, both of the above described processes can happen simultaneously, although one might expect that one of them is dominant. Moreover, the success of the reheating process in the MDO cosmological scenario depends essentially in how strong are the couplings with the electromagnetic field and with the boson field. In this respect, a further investigation is necessary in order to constrain the parameter space of the models in a cosmological perspective.

Conclusions
In the present article we have studied an unified cosmological evolution driven by a MDO fermionic field.
We have shown that the MDO inflation with a symmetry breaking potential can successfully achieve a number of efolds as large as 71, with an initial condition for the field given by the smallest possible quantum fluctuation. The slow-roll regime is maintained during the whole inflation with a constant parameter η 0.3 during almost the entire period. The torsion function plays and important role in the remaining 15 e-folds until the end of inflation generating a bump in the Hubble function at that epoch.
Such a sudden increase of the energy density of the Universe to a value greater than the initial value rises the question if it could generate any imprint in the cosmic microwave background anisotropies spectra. Similarly, one may expect a modification in the form of the primordial spectrum of gravitational waves. These issues will be answered by forthcoming investigations.
It is worth to stress that the dark matter behavior of the MDO fermionic field in the way it is understood in the present article is quite different from the preceding works on the topic [22,23]. In those works, the coherent oscillations of the field after inflation is interpreted as dark matter. This is because the oscillating field is pressureless in average in the same manner it happens with an oscillating scalar field subject to a quadratic potential. However, this is possible only in the case the MDO field interaction with any other field is negligibly small. Conversely, if the field couples with radiation or with the Higgs boson at the still high energetic regime after the inflationary epoch, than its energy density will decay exponentially in a reheating phase as discussed in general in the Sect. 5. In this sense, the MDO dark matter behavior only survives if the field acquires a non-null constant value after reheating which is of the order of the Planck mass or, more precisely, φ c = m pl Ω DM,φ /π Ω b . Therefore, if this condition is satisfied, it does not matter what kind of potential the field is subject (neither the inflationary mechanism), the gravitational coupling of the field with baryonic matter leads to a contribution in the form of a pressureless matter. It is worth to emphasize that such a coupling has its origin in the definition of the spin connections in curved spacetimes. No explicit coupling between the MDO field with other fields were assumed.
Moreover, if the potential is not null in φ = φ c , than the MDO field has an additional constant contribution to the energy density that works exactly as a cosmological constant. Such term is also responsible to constrain the mass of the field in the range 0 < m 2.4 × 10 10 GeV.
Finally, the scenario described in this article provide us with a natural explanation for the cosmological coincidence problem. This is because the evolution of the energy density of the MDO field scales with the energy density evolution of other matter fields, resulting in a present value of the same order of magnitude of the usual matter.
Therefore, we conclude that MDO fermionic field is a good candidate to drive the whole evolution of the Universe in an unified fashion, in such a way that the inflationary field, dark matter and dark energy are described by different manifestations of a single fermionic field. Additionally, as already pointed out at the Introduction, a possible interpretation of the inflationary phase as a consequence of the Pauli exclusion principle put this model on very interesting physical grounds.