Cosmology of the Nambu-Jona-Lasinio Model

We review the Nambu and Jona-Lasinio model (NJL), proposed long time ago, in the sixties, as a fermion interaction theory with chiral symmetry. The theory is not renormalizable and presents a symmetry breaking due to quantum effects which depends on the strength of the coupling constant. We may associate a phase transition with this symmetry breaking, leading from fermion states to a fermion condensate which can be described effectively by a scalar field. Our purpose in this paper is to exploit the interesting properties of NJL in a different context other than particle physics by studying its cosmological dynamics. We are interested in finding whether possibly the NJL model could be used to describe the still unknown dark energy and/or dark matter, from up to 95% of the energy content of the universe at present time.


Introduction
In the last years the study of our universe has received a great deal of attention since, on the one hand fundamental theoretical cosmological questions remain unanswered and, on the other hand we have now the opportunity to measure the cosmological parameters with an extraordinary precision.In the last decades, research in cosmology has revealed the presence of unexplained forms of matter and energy called Dark Energy "DE" and Dark Matter "DM" making up to 95% of the energy content of the universe at present time.The study of supernovas SNIa shows that the universe is not only expanding, but besides it is accelerating [1]- [6].Such behaviour can be explained by the existence of a new form of energy, Dark Energy with an anti-gravitational property, which would be explained by a fluid with negative pressure.Independent evidence for Dark Matter (DM) and Dark Energy (DE), is provided through the analysis of the Cosmic Microwave Background radiation (CMB) [7]- [10], which has been measured by satellite WMAP [11], and more recently by Planck mission [12], and the dynamics of galaxies, clusters and super clusters, and the study of the formation of Large Scale Structure [13]- [16] in the universe and weak lensing (the gravitational deviation of light), which point out the existence of matter that do interacts with ordinary standard model matter only weakly, as due to gravity.Other important measurements are the Baryon Acoustic Oscillations "BAO" [17]- [19].⋅ ⋅ [12].However, the nature and dynamics of Dark Energy and Dark Matter are topics of major interest in the field [20].The equation of state "EOS" of DE is at present time 0.93 0.13 o w − ±  but we still do not have a precise measurement of ( ) w z as a function of redshift z [12] [16].
Since the properties of Dark Energy are still under investigation, different DE parametrizations have been proposed to help discern on the dynamics of DE [20]- [23].Some of these DE parametrizations have the advantage of having a reduced number of parameters, but they may lack a physical motivation and may also be too restrictive.Perhaps the best physically motivated candidates for Dark Energy are scalar fields which can be minimally coupled, only via gravity, to other fluids [20]- [23] or can interact weakly in interacting Dark Energy "IDE" [24]- [27].Scalar fields have been widely studied in the literature [20]- [23] and special interest was devoted to tracker fields [22] [23] since in this case the behavior of the scalar field φ is very weakly dependent on the in- itial conditions at a very early epoch and well before matter-radiation equality.In this class of models the fundamental question of why DE is relevant now, also called the coincidence problem, can be ameliorated by the insensitivity of the late time dynamics on the initial conditions of φ .Nowadays there are a huge number of ideas aimed to explain these unknown cosmological fluids DE and DM, from the theoretical point of view, none of them being still conclusive.This situation supports and motivates our research.Given that our most successful theory of matter, the Standard Model of particle physics (SM), which is settled within the theoretical frame of Quantum Field Theory (QFT), it would be reasonable to ask a theory attempting to describe dark fluids to be based on QFT as well.
In this paper we study a fermion interaction theory with a chiral symmetry, the Nambu-Jona-Lasinio (NJL) model.Though this is an old and well known model in the context of hadron physics, it has interesting properties and it is worth to consider it with a different perspective, by studying its possible relevance for Cosmological Physics.Other examples of QFT models of DE and DM have been proposed using gauge groups, similar to QCD in particle physics, and have been studied to understand the nature of Dark Energy [28] [29] and also Dark Matter [30] [31].
We organized the present work as follows: In Section 2 we present the NJL model.In Section 3 we review the pertinent cosmological theory.Sections 4 and 5 present a study of the cosmological dynamics of a NJL fluid with a weak and strong coupling, respec-tively.In Section 6 we consider the addition of a cosmological constant to our NJL fluid, and analyze the different possible behaviours.In Section 7 we comment an interesting possible way to modify the original NJL model, obtaining an additional term in the effective potential which could be related with a Cosmological Constant.Finally, in Section 8 we summarize our results and present the conclusions.

The Nambu-Jona-Lasinio Model
Inspired by a, by then recently explained phenomenon in Superconductivity research, professors Y. Nambu and Jona-Lasinio, suggested that the mass of fermion particles (described by a Dirac equation) could be generated from a primary four-fermion self interaction, leading to a chiral symmetry breaking.The proposed Lagrangian, invariant under chiral transformations, has the form ( ) ( ) , 2 where ψ is a four-component spinor, and g is a coupling constant.From Equation (1) the four fermion interaction term is given by ( ) ( ) with no original mass term for the fermions.Since the coupling has dimension-2 in mass units, the theory is non-renormalizable.However, we are interested in considering the NJL model as an effective theory, useful below certain energy scale.The theory (1) describes a four-fermion interaction which can be expanded following conventional perturbation theory, and represented by Feynman diagrams (Figure 1).The infinite number of fermion loops can be resumed giving a non-perturbative potential.This can be easily done by introducing an auxiliary scalar field φ and an equi- valent Lagrangian for Equation (2) in the form The field φ plays the role of a Lagrange multiplier which can be eliminated using the Euler-Lagrange equations, ( ) . For the Lagrangian above we find where φ has mass dimensions, and by substituting Equation (4) in Equation ( 3) one can recover the original Lagrangian Equation (2).Note that we introduced the parame- ter m with a physical dimension of mass, so that mg is a dimensionless coefficient, and we have dimensional consistency for all the physical quantities. 1 The term 5 ψγ ψ in Equation ( 1) represents a pseudo-scalar quantity, and we have allowed ourselves to ignore the field contribution associated with it in the new Lagrangian in Equation ( 3), as we would like to start to study the simplest possible model. 2 From the equivalent Lagrangian one may read the fermion mass and the tree level scalar potential 0 V .We have respectively: ( ) The effect of quantum processes (represented by loop diagrams) may be taken into account through the well known Coleman-Weinberg potential ( ) the minus sign in 1 V is because it corresponds to the fermionic contribution to the Coleman-Weinberg potential.As we will see, for the strong coupling case, this will enable the effective potential to adopt a negative value when the field stabilizes at the minimum.The integral grows up indefinitely as the upper limit goes to infinity, i.e. it has an ultraviolet divergence.Because of the non-renormalizability of the theory, we cannot avoid this divergence, so we regularize by introducing a cut-off Λ .This para- meter defines the energy scale below of which the theory is valid.We define the x variable as and the potential becomes ( ) Notice that the one-loop potential 1 V is negative since it corresponds to the contribution of the original fermion field ψ , and we choose to parameterize it in terms of the effective scalar field φ c.f. Equation (4).
For the sake of concision we also define ( ) ( ) In this way, taking quantum corrections into account we obtain an effective potential given by ( ) 1 Remember that the dimension of a scalar field equals that of mass.

2
It will become clear that by doing so does not affect qualitatively the implied physical processes.
with the complete potential ( ) As a function of φ it can be written explicitly as ( ) Equation (13) gives the complete NJL scalar potential, and we are interested in studying its cosmological implications.Let us determine the asymptotic behaviour of the scalar potential V in Equation (13).To analyze the potential we seek for extremum points.For the function ( ) f x in Equation (10) we have the derivative ( ) and for the derivative of V we have The condition implies the following equations: The first one says that the origin 0 φ = is an extremum, and if we take the second derivative we see that if then the extremum at 0 φ = corresponds to a minimum, while for we have a maximum at the origin.The equation above suggest to define a critical value of the coupling c g as so that we see that for a weak coupling c g g < we have a minimum at the origin, while at strong coupling c g g > we have a maximum.The type of extrema at the origin of the potential corresponds to the value of the coupling.Now let us determine the second (possible) extreme of the potential.Since the r.h.s of the second equation in Equation ( 17) is negative (i.e.log 0 1 ) this equation has a solution only for a strong coupling c g g > .A value for x (or that of the scalar field φ ), at the minimum cannot be solved analytically, since the second equation in Equation ( 17) is a transcendental equation.One way to determine a solution is to seek for the intersection between the curve of the function log 1 r.h.s. in the second Equation (17), and the constant in the l.h.s.In this case do exist an intersection (only one, as the r.h.s. is a monotonic function), giving a solution for the x variable, leading in its turn to a non-trivial solution in min φ φ = which is a minimum. 3The extremum in this case corresponds to a minimum.Notice that in all cases we have at large x the limit V → ∞ for x → ∞ regardless of the value of the coupling g.
Therefore, we have: if  17) should be solved.However, since it is a transcendental equation in the variable x, an algebraic expression cannot be written, and we need to use numerical procedures.Let us introduce a parameter s to write g in the form 2π .
In this way we make sure to have a strong coupling by taking 1 s > .Now, for a given value of the coupling with c g g > , there exists a definite value of x, say 0 x , satisfying the condition ii) in Equation (17), which is the solution for the minimum.Therefore, the potential valuated at this point yields the minimized potential, i.e.

( )
Then, by substituting ii) Equation (17), and using Equation (20) in the expression for the potential Equation ( 12), we can write min V in the suitable form ( ) ( ) which provides a good idea of how min V is related to the energy scale.
From Equation (4) the field φ ψψ , is a Lorentz invariant quantity, so φ is scalar field.When the field φ is stabilized, a non trivial expectation value reflects the pres- ence of a fermion condensate.Now, if the field has an expectation value 0 φ = , it means that the state of paired fermions ψψ is not present, so we have a system consisting in the original massless fermion particles with a 4-Fermi interaction, and a condensate is not energetically favoured.This happens for a "weak" coupling c g g < .On the other hand, if the expectation value 0 φ ≠ , then we have a fermion condensate represented effectively by the scalar field.This happens for a "strong" coupling c g g > , and a fermion condensate is dynamically formed since it reduces the energy of the system.Thus, we see that two different fluid phases (massless fermions or fermion condensate) are obtained depending on the strength of the coupling.Next, we investigate the cosmological dynamics of each of these fluids.

Standard Cosmology
The widely accepted current standard cosmological model (the Big Bang theory) is based in Einstein's theory of General Relativity.If conditions of spatial homogeneity and isotropy are assumed, the space-time metric adopt the well-known simple form ( ) where the variables , , r θ φ are comoving coordinates parameterizing the spatial sec- tion of space-time, and k takes the values +1, 0, −1 for spaces of constant positive (spherical), zero (flat), or negative (hyperbolic) curvature.When this metric is used in the Einstein's equations, the so called FRWL equations (Friedmann-Robertson-Walker-Lemaitre) can be obtained.As these assumptions agree with observations 4 to a very high precision, we will use this same theoretical framework.Because the necessary equations are well known and their deduction can be found in standard text books, in the 100 Mpc  .following we limit ourselves to write them and to give only a brief explanation.
The equation relates the expansion rate (in time) of the scale factor a, and the curvature k of the universe, to the total energy density ρ .Along this paper we will always take a flat geome- try 0 k = , as suggested on the one hand from the theory of early cosmological inflation, and on the other hand (and most important) from observation of the CMBR.
Introducing the usual definition relating the Hubble parameter H with the rate of change in time of the scale factor a , a a H = ⋅  Equation ( 23) (with The continuity equation for a fluid with energy density ρ and pressure P is ( ) For a perfect fluid "α" satisfying a barotropic equation of state P w α α α ρ = , with w α a constant, Equation ( 26) can be solved analytically.We sometimes will refer to such a fluid with the name of "barotropic fluid".From the cosmological point of view, the substances contained in the universe can be described as radiation, which has , .
A scalar field φ , with a self-interaction potential ( ) V φ , has energy density φ ρ and pressure P φ given by where we have also defined the kinetic energy k E in the third equation.Considering an universe containing radiation, matter and a scalar field, the total energy density is written .
For a given component fluid "α", it is useful to know its relative density, defined as the ratio of its energy density to the total energy density: where we have used Equation (25) in the second equality.In a flat universe one has the condition 1.
It is interesting to note that while Equation ( 31) remains valid even when we have a negative α ρ , the quantity α Ω is no longer constrained to the values 0 1 α ≤ Ω ≤ .In the work presented here, the fluids can have a negative energy density, giving 0 α Ω < , or a total energy density ρ that vanish at finite values of the scale factor ( ) a t , in which case we would have α Ω → ±∞ .
Taking the time derivative in Equation ( 25) and using Equation ( 26), it can be found Note that the r.h.s. in Equation ( 32) is always negative.The equation of motion for a spatially homogeneous scalar field, (a modified Klein-Gordon equation) is given by It is also useful an equation for the acceleration of the scale factor: ( ) Differential Equations ( 24), ( 32), ( 33), together with (27) constitute a complete set which can be solved numerically (since we cannot always write an analytical solution).Nevertheless, it is convenient to attempt to outline the general behaviour of the dynamical system.Thus, before going to solve for our NJL potential, let us point out the following generic facts: The evolution of the scalar field is such that it will minimize the scalar potential ( ) V φ , so for an arbitrary initial value i φ , the field will roll to lower values of the poten- tial, in such a way that eventually it will adopt a constant value ( min φ φ = being the minimum).Given than the scale factor is a positive defined quantity, the energy densities for matter and radiation Equation ( 27) are always positive quantities and never equal to zero for finite values of the scale factor ( ) a t .So, the total energy density Equation (29) remains always positive as long as the condition is satisfied.Thus, Equation (25) says that 0 H = , that is 0 a =  , never happens (Equation ( 24)) as long as 0 ρ ≠ .This implicates that 0 a >  always.This means that the scale factor ( ) a t never reaches an extremum value along its time evolution (taking an initial condition 0 i H > , since we know that the universe is expanding at present time).Nevertheless, it is interesting o note that there is no known physical principle forbidding the existence of a fluid with a negative potential ( ) 0 V φ < , at least for some val- ues of the field φ .In this case, it could well happen that Equation (35) become an equality, meaning 0 ρ = for finite values of ( ) a t , which in turn implies 0 H = , and 0 a =  ; i.e., the scale factor reaches an extremum value (indeed a maximum, since as seen before, it was initial growing).Now, Equation ( 32) imposes an always decreasing Hubble parameter H (because the right hand side is always negative), so that after being 0 H = it must be 0 H < , and therefore 0 a <  , i.e. the scale factor decrease.In other words, the universe must be contracting after reaching its maximum size.Observe that this result is a consequence only of the negativity of the potential, and it is independent of its specific form.This collapsing universe is valid even for a flat universe 0 k = .To conclude, if a fermion condensate is energetically favored then the minimum of potential ( ) V φ is negative and the universe will recolapse.

Dynamics of Massless Fermions Phase (Weak Coupling g < g c )
As we have seen in Section 2, for a weak coupling c g g < the minimum of the potential ( ) is located at the origin with min 0 φ = , and V does not take negative values.Therefore, the total energy density and H never vanish for finite values of the scale factor a, and we have 0 a >  due to Equation (24).So the scale factor ( ) a t is always growing, going to an infinite size in an infinite time.Now, from Equation (34), it can be seen that, in order to have 0 a <  , i.e. the universe to slow down its expansion rate, then ( ) is a condition to be satisfied.This, of course, in not always the case: we could take an initial field amplitude i φ as big to make the initial value of the potential ( ) big enough so that inequality (36) does not hold, and we would have instead In this case we could have an acceleration of the scale factor, i.e.
an accelerating universe, though it would be an "early" acceleration, as it would be present an initial times, i.e. before letting the fluid densities to dilute and field to evolve.
As time passes, the field rolls down minimizing the potential, and eventually acquires some value i φ φ < such that condition (36) becomes fulfilled.Given that the densities of matter and radiation never reach a null value in a finite time, and that the field amplitude tends to be stabilized around the minimum (i.e.0 φ → ), for a big enough amount of time, we expect a vanishing potential and velocity, ~0 V , ~0 φ  to be a good approximation to a final situation, in which ( 36) is still satisfied.
We show an example of numerical solution in the figures.In Figure 3 we see that the field has a damped oscillation around 0 φ = , and in consistency with this, its kinetic energy (velocity) diminish in time and we show in Figure 4 the evolution of the relative densities , , Ω Ω Ω for radiation, matter and φ .Simultaneously, the potential valuated at ~0 φ goes to lower values (according to ( ) We can see that although the universe is expanding, it always ends up in a non-accelerating regime (Figure 5).A Taylor expansion for the potential about where the whole coefficient multiplying on 2 φ , is a positive quantity, as 2 2 2 4π g < Λ .
The coefficient of state φ ω defined below Equation ( 26), for the field φ , written ex- plicitly is   .
Since at late times, when the field oscillates around its minimum with a quadratic potential, the average value is 0 Within the context of Early Cosmic Inflation theory, the so called Slow Roll parameters are defined as follows: which have to satisfy the conditions η < in order to the potential may cause a positive acceleration.Even though they are valid for a single field, without additional fluids (matter and/or radiation), we show them in Figure 6 the Slow Roll parameters, for the seek of completeness.

Fermions Condensate Dynamics (Strong Coupling, g > g c )
The strong coupling case leads to a fermion condensate and therefore to a negative potential V at its minimum.The potential has at the origin ( ) when minimized.If the initial velocity , then the kinetic energy of the field has a null value as well, so we have for the initial energy density i i V φ ρ = .The initial amplitude for the field i φ cannot be such that makes ( ) 0 i V φ < , because it would lead to an im- aginary value for H, according to Equation (25).Thus, we must take always i φ such that 0 i V > .As before we begin with 1 0 i H = > , therefore Equation ( 24) says that ( ) a t initially is increasing in time.The equation ( 32) is written , so that H always diminish in time.As the potential is minimized, it goes from positive to negative values, and from Equation ( 25) eventually it will be 0 H = , and after this 0 H < , corresponding respectively to 0 a =  and 0 a <  .In words this means that after an initial period of expansion (increasing scale factor), a maximum value is reached, followed by a period of contraction.Since H  remains always negative, then ( ) a t will con- tinue decreasing, so that it necessarily will collapse.In other words, it will be 0 a = in a finite time in the future (because the evolution is forward in time: the field minimizes, not otherwise).Now, while the expanding phase is taking place, the field is rolling down, eventually entering in a damped oscillatory regime nearly the minimum, where the potential has become negative, min 0 V < .Because of the damping, the kinetic energy tends to a zero value, 0  a t is decreasing in the contracting phase, the energy densities behave the opposite way, i.e. they all increase in time.Therefore, we expect 0 φ ρ = to happen twice.In its turn, this implicate that the coefficient of state φ ω , Equation become a divergent quantity also twice, this two points, and near them, φ ω is not anymore a useful parameter to characterize the fluid represented by the field φ .Below we show a numerical solution example (Figures 7-11).
As we mention before, in Section 3, a similar circumstance arises in dealing with the relative densities α Ω : it is considered that in order to this parameter to make sense, a relative density should adopt values 0 1 α ≤ Ω ≤ .However, as can be seen in Equation (30), if at some time is 0 H = , then nearly this value, each α Ω turns into a divergent variable.The situation is even weirdest for the field, because near the minimum it is min ~0 V φ ρ < , the energy density of the field is similar to the potential, which is nega- tive.This would make α Ω → −∞ (a divergent and negative relative density!).
Consider now a universe containing matter and radiation in addition to our NJL fluid.An interesting question is, may the presence of these fluids prevent the universe to collapse?Remember that the condition for an increasing scale factor can be reduced to the inequality (35).If the scale factor is supposed to grow forever, this condition must be hold always.Now, according to the explanations given above, initially the scale factor is growing indeed.Thus, from Equation (27) we see that the densities of both barotropic fluids (matter and radiation) must be decreasing.At the same time, because the field is stabilizing in the minimum of the potential, the kinetic energy of the field  is diminishing to zero, whereas the potential is going to a constant value , in such a way that necessarily, condition (35) ceases to hold.Therefore, even in presence of additional barotropic fluids (does not matter the relative amount with respect to that of the fluid associated with the field), the collapsing universe situation cannot be avoided.
The previous qualitative generic analysis is verified by the numerical solution for our NJL potential in particular (Figures 7-11).By observing the graphics, we found an unpredicted interesting non-trivial behaviour of the field amplitude: while the scale factor undergoes the expanding, and contracting phases successively, an damped oscillating phase around min φ is taking place, as expected.But then, at some point in the contracting phase, the field amplitude goes to bigger values, and as the scale factor approach to 0 a = , the field is taken out from the minimum and it begins to increase monotonically! 5 Is this an acceptable result?Intuitively, as a is decreasing, it is reasonable to expect all densities to be growing.In particular, if the field density is getting bigger, it should be due to an increase in the field velocity (so k E gets bigger), or in the field amplitude (so V gets bigger); or both.This behaviour can indeed be      Figure 8 show both kinetic and potential energies, and we can see that even though the initial kinetic energy is zero, it overtakes the potential energy and remains so until the collapsing moment final t when ( ) 0 a t = at late times.Nevertheless, the potential energy also grows as the time is approaching final t , so the field amplitude is eventually expelled from oscillating about the minimum.

NJL Fluid with a Cosmological Constant
Due to its theoretical properties and observational requirements, a Cosmological Constant is a very usual and useful ingredient included in cosmological models, and it is worth to consider such contribution in our model.Its defining property is an energy density ρ Λ which does not vary in time, and a coefficient of state gives a pressure P ρ Λ Λ = − .In a universe containing only a Cosmological Constant, the equation ( 34) is written which, as 0 ρ Λ > , implicates ( ) 0 a t >  always.Therefore, such an universe is always accelerating its expansion.In fact, in this case the Equation ( 24) may be solved analytically, after substituting Equation ( 25), giving the well known solution ( ) ( ) . How do the presence of a Cosmological Constant affect our previous considerations of a universe including our NJL fluid, besides matter and radiation components?Will the universe accelerate or collapse, even in the presence of a scalar field with a negative potential 0 V < ?Because the density ρ Λ is constant, we have that the differential equations are not modified, other than just adding a term in the expression for H, equation (25).In particular, the equation of motion Equation ( 33) remains unchanged, so the field dynamics is not affected.As before, we have to deal with two cases.a) Free Fermions ( As studied before, the potential is 0 V ≥ , and its min- imum value is min 0 V = .Also, with the pass of time, both matter and radiation densities dilute, going to vanish.From Equation ( 34), it can be deduced the condition for universe to decelerate: ( ) ( ) Given that the left hand side in this inequality is diminishing in time, whereas the right hand side remains constant, we have that eventually this inequality cannot hold anymore, and becomes an equality, meaning 0 a =  .This points the beginning of the ac- celeration period, i.e. 0 a >  , where the inequality (42) gets inverted.Had the initial conditions been such that inequality (42) were the opposite, then there would be always an acceleration holding always, because the LHS would never go back to grow.
Thus, we see that for a free fermions NJL fluid with a Cosmological Constant, the universe necessarily accelerate, the precise moment depending on the amount of energy densities m ρ , r ρ , with respect to that of ρ Λ .This can be specified in the initial con- ditions, which in their turn can be chosen to solve for a realistic model fitting the observations.b) Fermion Condensate ( c g g > ).We found before that, for a strong coupling, the potential is negative when minimized, min 0 V < .Do the universe necessarily accelerate also in this case?In order for this to happen, condition (42) eventually must turn into an equality, meaning 0 a =  .This is a minimal condition to be satisfied, because it points at least the beginning of an acceleration; it remains to be sure that acceleration will be sustained.Let us label all quantities with a subindex "ac" at time ac t , when 0 a =  (vgr.( ) ).From Equation ( 36), we have6 ( ) Remember that the potential take positive values as well as negative ones, so both possibilities must be taken into account.Certainly one can find such set of values of V for a given ρ Λ to satisfy the inequality.However, if we rather want to consider realis- tic models, we should consider plausible values from observations (besides, we would not like to complicate our lives by considering unrealistic generic situations).
From definitions (30) it can be found that ( ) , where z is the redshift, and 0 0 r m r = Ω Ω says the amount of radiation with respect to that of matter.The subindex "0" refers to current values, i.e. quantities measured "today".Now, the estimate for z (the time when acceleration begins) is around ~1 z ; and it been measured


. Now, remember that a decelerating period dominated by matter is supposed to have taken place before 0 a =  .In order for this to happen, condition (42) should have to be true before condi- tion (43).For ~1 z (it could be even as big as, let's say ~10 z , as this would not change the essence of the argument) and using condition (42) we would have in order to be 0 .
If a positive acceleration eventually come up, the above expression is expected to become an equality.Now, suppose 0 V > .Then, unless k E decrease even fast, the RHS in the inequality should be decreasing as time passes, because the potential is minimizing.But k E cannot behave like that indeed, as the field is under a damped rolling, not to mention that k E is never a negative quantity, so the sum of terms 2 k V E − will end up decreasing (would the values of these terms been such that the equality somehow would be accomplished at some time, in this case the acceleration could not be attached to ρ Λ anyway).On the other hand, for 0 V < , the inequality would become even more strong in time, because again, the potential is minimizing: , and . Therefore, if initially the inequality (44) begins being satisfied, it will remain being so always; in other words, the universe will never accelerate.
What about a collapse in the future?May the presence of a cosmological constant prevent a decreasing scale factor (time going forward)?For a growing scale factor we have 0 a >  , which is true indeed because we take 0 i H > is the initial value of H. 7 As we explained before, if the scale factor is to reach a maximum max a a = , it must be 0 a =  .Let us name am t the time when this is accomplished (if so), and label with a subindex "am" the variables valuated at this time.We have for the total energy density The only way in which this could happen is for 0 , so the equation, as a condition to be satisfied by ρ Λ , can be written in the more intelligible form ( ) If we want to keep our analysis as simple as possible, we may ignore the contribution from radiation, 0 ram ρ = (observe that, had an acceleration would be possible, then we should assume ac am t t < , i.e. acceleration before receding, otherwise the model would not be useful.So, if Now, nothing forbids to exist a potential sufficiently deep min 0 V < , so that the equality (45) can be accomplished.The exact time at which this is achieved will depend on the relative amounts kam E , mam ρ , with respect to ρ Λ , i.e. on the initial conditions.
However, we can estimate a limit value by making 0 mam ρ → , 0 kam E → , and a stabi- . Then we have

( )
Maximum allowed value for the universe to collapse .
After 0 a =  , i.e.H = 0 (Equation ( 24)), the universe must enter into a contraction phase because H is always decreasing (Equation ( 32)), meaning . So, eventually the universe collapse in the future in a finite lapse of time.For , the scale factor would never go to contract, as in this case the total energy density ρ would never vanish.
It is interesting to observe that a Cosmological Constant may be seen as a particular case of a scalar field evolving under a potential stabilized with a positive minimum.As we have seen, the NJL model has two different behaviours depending on the value of the coupling constant g.For weak coupling c g g < the potential ( ) V φ has a mini- mum at the origin with ( ) and ( ) 0 V φ ≥ otherwise.On the other hand, at strong coupling , 2 with o V a constant value (it would be 0 o V = at weak coupling and 0 o V < at strong coupling) and o φ a constant.We can now ask ourself if we can have an accelerating universe.The evolution of the scalar field is just corresponds to a massive scalar field with energy density   47)) a barotropic fluid, which we now take for simplicity as matter (without lose of generality), and a cosmological constant ρ Λ , we have ( ) . Since the potential o V vanishes at weak coupling and is negative at strong coupling, there is a cancelation between the two cosmological constants ρ Λ and o V , and the NJL model plays therefore against an accelerating phase around the minimum of the potential, since o V is negative.

Dark Energy from NJL and SUSY Gauge Theory
As we have seen until now, the original NJL model has interesting cosmological conse-quences.However, the model by itself does not reproduce the observed feature of an accelerated expansion of the universe, and it is not desirable to introduce a cosmological constant by hand, without a good explanation.We rather ask for any model to be motivated from a deeper fundamental theory.Nowadays, a paradigm for such fundamental theory is played by Super Symmetric Field Theories, and a lot of work has been done in attempting to explain Dark Matter as well as Dark Energy as some super symmetric particle (references are given in the introduction, sec.1).Nevertheless, any conclusive theory has been established yet to present date.We would like now to generalize the NJL potential to include a physically motivated potential from supersymmetric gauge theories.These class of models have been previously studied in Dark Energy models derived from gauge theory [28] [29] (and references therein), [21] [28]- [35] and are based on ADS (Affleck-Dine-Seiberg) superpotential [35]- [38].The derived potential is of the form which is obtained from a non-perturbative super potential in a gauge theory, e.g.

( ) ( )
with f N flavours, and φ represent a fermion condensate, i.e. φ ψψ ≈ .The condensation energy Λ  is the scale of breaking of the gauge symmetry. 8We now add the potential in Equation (48) to our NJL model.Since the effective NJL potential in Equation ( 13) has a quadratic term .Then, at some lower scale Λ , the self interaction of the field φ becomes more involved and the dynamics of the field is also governed by the effective NJL potential.By adding Equations ( 13) and (48), we would have the total potential ( ) ( ) shown in Figure 12.Of course, this is an effective theory which is plausible to the ex-Figure 12. Graph of the total potential V = VNJL + U, as a function of the variable x.
tent that NJL and Equation (48) are valid or useful theories.This is on the same footing than using the NJL model to study the dynamics of hadrons, without having obtained the model directly from the QCD Lagrangian.Since we are simply adding a term to the already studied NJL potential, we use the results of the previous section 2. Using eq.ii), (17), the condition to be satisfied by the minimum x is now written (remember that wrote before where we have defined the function ( ) x α , which is seen to be parameterized through . This function has 1 α −∞ < < + and it is a monotonous growing function, regardless of the values of the parameters (they all are positive definite).Now, we know that the NJL potential NJL V is minimized in a non-trivial minimum min 0 φ ≠ when c g g > , or 1 s > .In this case, the LHS in Equation ( 50) is a positive quantity, corres- ponding to 0 1 α < < in the RHS, determining a solution 0 x x = .This means that the total potential (49) still is minimized for some 0 x , giving in its turn a non-trivial 0 φ .
Given that the minimum 0 x satisfies (50), the minimized potential can be written ( ) From this equation we see that, it is possible to obtain min 0 V > (which would behave like a cosmological constant), if the parameters satisfy ( ) ( ) 8π m Λ = Λ  .In this way the Equation (50) is written which has the solution 0 1.83 x x =  .Then, RHS Equation (52) gives the number ( ) This means that, in order for the potential to be positive at the minimum, the parameters must satisfy .We can use Equation (51) to obtain min V ; in order of magnitude we have ).Given that this theory allows min 0 V > as a result of the dynamics of the field, we have a possible explanation for the presence of a cosmological constant, and an accelerating universe.

Summary of Results and Discussion
The fermion model of Nambu and Jona-Lasinio (NJL) includes two different fermion states resulting from quantum effects, each one being associated with two different physical phases.For a weak coupling c g g < we have massless fermion fluid, whereas for a strong coupling c g g > a massive fermion condensate fluid is obtained.In this later case we can determine the mass of fermions and it is due to non-perturbative effects due to the strong coupling.A very convenient way to describe the system is to consider an equivalent scalar field φ moving under an effective potential + , which has a different form depending on the coupling strength.
Notice that in the strong coupling case c g g > , the potential has a non-trivial negative minimum due to the negative contribution one-loop potential 1 V in Equation ( 6).
The negative sign of this potential is due to the fermionic origin of ψ field, and we have chosen to parameterize the fermion condensate in terms of an effective scalar field φ ψψ , as in Equation ( 4).
Here we studied the potential and solved the cosmological evolution for each fluid in presence of additional barotropic fluids (e.g.matter-dust or radiation).
For a weak coupling, we found a coefficient of state φ ω with oscillating values around zero, in such a way that the average value 0 φ ω = .Also, because the potential goes as 2 Ṽ φ near the minimum, we have that the NJL fluid in the form of free fer- mions dilutes as a matter.A universe containing such a fluid (with or without matter and/or radiation) will expand forever without accelerating.On the other hand, a universe containing this NJL fluid besides a cosmological constant (with or without matter and/or radiation), will eventually accelerate necessarily, expanding forever.
On the other hand, the strong coupling case (without a cosmological constant) always causes an eventually vanishing energy density.This is due to the fact that the potential is negative when minimized, and even the additional presence of matter and/or radiation does not prevent this to happen.Since the vanishing energy (which is associated with the scale factor reaching a maximum), is followed by a contracting period, this means that a fermion condensate always makes the universe collapsed.The energy density of the field φ ρ vanishes a couple of times (one in the expanding phase, and the another one in the contracting phase).Because of this, some quantities ( φ Ω , φ ω ) be- come inadequate to describe the fluid.It is important to point out the following interesting fact: Equation ( 23) has been known and well studied since long time ago.If the curvature parameter is 1 k = + , the universe is said to have a spherical geometry; the scale factor is expected to get a null value eventually, so we have a collapsing universe.Because a spherical universe is also finite or closed, a collapsing universe was always associated with a closed universe.On the other hand, if 0 k = , the universe has a flat geometry.
For ordinary matter the total energy density could be diminishing, but it could never vanish effectively in a finite time, so the scale factor in this case is expected to be always increasing.Because divergent geodesic lines in a plane never meet again, a flat universe is said to be open.So, an open universe was thought to be infinite in size (although not necessarily, but in any case, always growing).Now, remember that from the beginning, in our present study, we have taken the curvature parameter to be 0 k = , so we have been treating with a flat universe all the time.Nevertheless, we found that, if the universe contains a scalar field with a negative potential, then a future collapse cannot be avoided, giving a collapsing flat universe!In particular, because a negative potential arises naturally for the NJL model, a collapsing flat universe is also a natural consequence.
We also studied a variant of the strong coupling model, consisting in the addition of a cosmological constant.We found that, if the energy density ρ Λ is not big enough to overtake at least the minimized potential min V , the eventual receding of the scale factor cannot be avoided, and the universe will collapse inevitably.But if ρ Λ exceeds min V , then the scale factor will accelerate eventually, and the collapse will be absent.Perhaps it is worth to emphasize that, in both cases of weak and strong coupling and without considering a cosmological constant, one may induce an acceleration of the scale factor by manipulating the initial condition for the field amplitude i φ , but we do not interest in it because 1) it has to be fine-tuned, and 2) it does not allow to include realistic models in which a previous deceleration period of matter dominance took place.
It is important to keep in mind that, once we settle a coupling strength (weak or strong), there is nothing in the theory to allow to switch between them, so actually a phase transition cannot be considered.
A very appealing feature of the NJL model is, in our opinion, the fact that 1) it is based on a "fundamental" symmetry (chiral symmetry), 2) the model leads to a potential which, due to quantum corrections, can adopt negative values in a natural way, and 3) it includes only one parameter: the coupling constant g (two parameters if we count the cut-off Λ ).In return we obtain interesting consequences, as allowing more than one physical phase (each having different cosmological implications), and the possibility of a collapsing universe.This is to be compared with other models involving a symmetry breaking 10 or introducing new kinds of fluids aimed to be relevant to cosmological problems, but at the expense of introducing several fields or parameters. 11 Finally, we saw that by considering an additional term besides the NJL potential, in 10 For instance in Higgs-like models are required two parameters "m" and "λ" in order to get a potential  For instance, to "justify" the existence of scalar fields with useful potentials, frequently one has to invoke more sophisticated theories, like String, Kaluza-Klein, GUT's, etc. which demand a bigger effort to derive relevant results, and often implicate new exotic physics.the form of an inverse power (which is motivated from some supersymmetric theories), then it is possible to obtain a total potential with a positive minimum, thus allowing to explain a cosmological constant as a consequence of a field dynamics, which is a fermion particle (instead of a scalar field) governed by simple basic symmetries.
It has been established that our universe is flat and dominated at present time by Dark Energy "DE" and Dark Matter "DM" with 0
value separating between both behaviours of the potential (in Figure2we show all the three cases , , f. Equation (4).To estimate the value of the potential at the minimum for c g g > , the equation ii) in Equation (

Figure 2 .
Figure 2. Effective potential (13) as a function of φ .The critical value of the coupling c g , separates two kinds of behaviours.

Figure 3 .
Figure 3. Left: Scalar field amplitude φ .Right: State equation coefficient φ ω .Both variables are shown as functions of time.

Figure 4 .
Figure 4. Left: Relative densities α Ω for radiation, matter and φ .Right: Total relative density for barotropic fluids (matter and radiation) and for the field φ .The horizontal axis in both graphics represents time.Note that we show a different scale of time in each plot for the same solution.

Figure 5 .
Figure 5. Left: Scale factor ( ) a t .Right: Acceleration of the scale factor ( ) a t  .Both variables are , eventually passing from negative to positive values.Let us consider at first the simpler approach of a universe containing only a scalar field ( 0 r m ρ ρ = = , i.e. no additional fluids), evolving under a generic potential possessing a negative value ( )

kEFigure 6 .
Figure 6.Slow roll parameters  (dashed-red curve), and η (continuous-blue).Left: From 0 φ = to The total energy density, as well as the individual densities for each fluid (if there were additional fluids), would go to diminish in time (as can be seen for radiation and matter in Equation (27) with ~n a α ρ − , and ( ) a t increasing).By a similar rea- soning, because ( )

Figure 8 .
Figure 8. Left: Although the kinetic energy (red-upper curve) is zero initially, it overtakes the potential energy (blue-lower curve) and remains dominant all the way even to the collapsing time when ( ) 0 a t = .Right: The field oscillates around min φ and is becoming divergent as getting

Figure 9 .
Figure 9. Left: Hubble parameter.It is a null quantity about 200 t  approx.Right: Relative density of the field.As

Figure 11 .
Figure 11.Left: Scale factor.Right: Acceleration ( ) a t  .Both plots are to be interpreted as de- scribing a universe which expands without acceleration (note that ( ) a t  is never greater than zero), reaching a maximum value about 200 t  approx., thereafter falling in a contracting phase all the way long to collapse.
seek of simplicity, here we are interested only in orders of magnitude).Then we have negative minimum ( ) min 0 V φ < .So let us approximate the potential V around the minimum and take the ansatz massive scalar field may accelerate the universe only at large values of φ′ (larger than the Planck mass) when the Slow Roll parameters  and η are smaller than one, while at a late time when the scalar field oscillates around the minimum the energy density φ ρ ′ redshifts as matter, i.e.
field (with potential given in Equation ( us take 2 n = in Equation (48) so that we have the symmetry under 1 φ φ → , we need to say something about the parameters, so let us take 4 2 3 Let us now estimate some real physical values.The total energy density today is about , in the limit of stabilized fields about the minimum9 the energy density of our NJL fluid is value of coefficient  depending on the value min V , as shown above (for our example we have 2 ~10 −  have to have a "correct" relation between them in order to break the symmetry.

11
Since we have seen that H  is negative, this implies that H is always a decreasing function of time.If H can vanish