Starobinsky-like inflation and running vacuum in the context of Supergravity

We describe the primeval inflationary phase of the early Universe within a quantum field theoretical (QFT) framework that can be viewed as the effective action of vacuum decay in the early times. Interestingly enough, the model accounts for the"graceful exit"of the inflationary phase into the standard radiation regime. The underlying QFT framework considered here is Supergravity (SUGRA), more specifically an existing formulation in which the Starobinsky-type inflation (de-Sitter background) emerges from the quantum corrections to the effective action after integrating out the gravitino fields in their (dynamically induced) massive phase. We also demonstrate that the structure of the effective action in this model is consistent with the generic idea of renormalization group (RG) running of the cosmological parameters, specifically it follows from the corresponding RG equation for the vacuum energy density as a function of the Hubble rate, $\rho_{\Lambda}(H)$. Overall our combined approach amounts to a concrete-model realization of inflation triggered by vacuum decay in a fundamental physics context which, as it turns out, can also be extended for the remaining epochs of the cosmological evolution until the current dark energy era.


I. INTRODUCTION
In the last two years we have witnessed extraordinary developments on experimental tests of inflationary models [1], based on studies of photons in the cosmic microwave background radiation. In particular, the results of Planck collaboration [2] and the associated non-observation of B-mode polarizations of primordial light fluctuations, have imposed very stringent restrictions on single scalar-field models of slow-roll inflation, allowing basically models with very low tensor-to-scalar fluctuation ratio r = n T /n s ≪ 1, with a scalar spectral index n s ≃ 0.96 and no appreciable running. In fact the upper bound set by Planck Collaboration [2] on this ratio, as a consequence of the non-observation of B-modes, is r < 0.11, but their favored regions point towards r ≤ 10 −3 . This is a feature that characterizes the so-called Starobinsky-type (or R 2 -inflation, with R denoting the scalar space-time curvature) inflationary models [3]. The estimated energy scale E I of inflation, which in inflaton-type models is related to the -approximately constant-scalar potential during inflation through E I = V (2) The above can be rephrased as H I /m P 6.39 × 10 −6 , where m P = 1/G 1/2 = √ 8πM Pl ≃ 1.22 × 10 19 GeV is the Planck mass in natural units. This result is consistent with the well known CMB bound H I /m P 10 −5 on the temperature fluctuations induced by the tensor modes. As we will see, the actual value of H during inflation for the class of models under study satisfies H H I , and hence the CMB bound is preserved by them.
The recent joint BICEP2-Planck analysis [5] confirmed the early Planck result, namely the likelihood curve for r yields an upper limit r < 0.12 at 95%. Moreover, the present BICEP2-Planck data are consistent with a scalar spectral index n s ≃ 0.96 and no appreciable running, in agreement with the previous Planck data [2]. Using the aforementioned new upper limit r max = 0.12, the Hubble parameter during slow-roll inflation H I is estimated to be below H Bicep2+Planck and hence H I /m P 6.64 × 10 −6 . Because of the low significance of the new limit on r, the possibility that r is actually much smaller than the current upper limit r max remains as natural as it was before. In fact, nothing actually prevents at present that the typical value of the tensor to scalar ratio can be, for example, r = O(10 −3 ), and in this sense the Starobinsky-type scenarios can still be considered as a serious possibility to describe the inflationary universe. Following this point of view, we continue in this paper with the investigation of Starobinskylike models as potential candidates for realistic implementation of inflation compatible with the data.
In previous publications one of us (N.E.M.) with collaborators [6] discussed the dynamical breaking of supergravity (SUGRA) theories via gravitino condensation, and demonstrated [7] the compatibility of this scenario with Starobinsky-like [3] inflationary scenarios. As we discussed, this phase is characterized by the dynamical emergence of a de-Sitter background. As argued in [6], the Starobinsky-type inflation appears much more natural (from the point of view of the order of the parameters involved) than a hill-top inflation scenario [8] in which the gravitino condensate itself is the inflaton field. In the latter, very large values of the wave function renormalization of the condensate field are required to ensure slow-roll inflation if one insists on (phenomenologically realistic) sub-Planckian supersymmetry breaking scales. It is important to notice at this point that in the original Starobinsky model [3] the R 2 terms crucial for inflation arise from the conformal anomaly in the path integral of massless (conformal) matter in a de Sitter background, and thus their coefficient is arbitrary, and can only be fixed phenomenologically. A similar, although not identical, situation occurs in the context of anomaly-induced inflation [9,10], where the term R 2 is absent at the classical level but is generated from the conformal anomaly. In this case, however, the coefficient of R (entering the β-functions and controlling the stability of inflation) presents also some arbitrariness which can only be fixed by a special renormalization condition. Par contrast, in the considered SUGRA scenario, such terms arise in the one-loop effective action of the gravitino condensate field, evaluated in a de Sitter background, after integrating out massive gravitino fields, whose mass was generated dynamically. The order of the de Sitter cosmological constant, Λ > 0 that breaks supersymmetry, and the gravitino mass are all evaluated dynamically (selfconsistently) in our approach from the minimization of the effective potential. Thus, the resulting R 2 coefficient, which determines the phenomenology of the inflationary phase, is calculable [7].
In regard to the early universe we emphasize that the RVM defines a class of non-singular inflationary scenarios with graceful exit into the standard radiation regime. These models are related to Starobinsky inflation models, although they are not equivalent. We will discuss in this paper the correspondence between them, and most particularly with the dynamically broken SUGRA model with gravitino condensation that we have mentioned above. It is especially remarkable that such specific implementation of the SUGRA model leads, as we will show in this paper, to the effective behavior of the RVM with calculable coefficients. In this way the former automatically benefits from the successful consequences of the latter. Let us mention that the RVM also provides some important clues for alleviating the cosmological constant problem [14].
Finally, we should like to mention that the RVM's have been tested against the wealth of accurate SNIa+BAO+H(z)+LSS+BBN+CMB data -see [30] for a recent summary review -and they turn out to provide a quality fit that is significantly better than the ΛCDM. This fact has become especially prominent in the light of the most recent works [31,32]. Therefore, there is every motivation for further investigating these dynamical vacuum models from different perspectives, with the hope of finding possible connections with fundamental aspects of the cosmic evolution. In point of fact, this is the main aim of this work.
The structure of the article is as follows. The general framework of the RVM is introduced in II. The basic theoretical elements of the Starobinsky inflation are presented in section III. The main properties of the dynamical breaking of local SUGRA theory and its connection to Starobinsky-type inflation are reviewed in section IV. In section V we demonstrate how the RVM describes the effective framework of the Starobinsky [3] and the dynamically broken SUGRA [8]  It is the purpose of this work to go one step further from demonstrating compatibility of the dynamically broken SUGRA scenario [6,7] with inflation and discuss the possibility of a dynamical evolutionof the inflationary phase ground state to the standard radiation regime within the context of the running vacuum model (RVM) of the cosmic evolution utilizing an effective "renormalization group (RG) approach", see [10,11] -and [14,15] for comprehensive expositions. Specifically, we wish to show that the behavior of the aforementioned SUGRA scenario effectively mimics the RVM. Once this link is elucidated, the general "decaying" vacuum description inherent to the RVM formulation allows to smoothly connect inflation to the standard Fridman-Lemaître-Robertson-Walker (FLRW) radiation era, which subsequently proceeds into a matter and dark energy domination in the present era, in which it still carries a mild dynamical behavior compatible with the current cosmological data [21,22]. Such an expansion history of the Universe has been put forward in previous works by the authors in various collaborations and contexts, e.g. non-equilibrium string-inspired cosmologies [33] or conventional field-theoretic cosmologies in which the above mentioned RVM is extensively applied for the study of the early cosmic history [16][17][18]20].
In the effective RG approach underlying the RVM one can write down an evolution equation for the effective vacuum energy density ρ Λ (t) = ρ Λ (µ c (t)), treated as a dynamical quantity whose cosmic time evolution is inherited from its dependence on a characteristic cosmic scale variable µ c = µ c (t). This variable plays the role of running (mass) scale of the renormalization group approach, and a natural candidate for such scale in FLRW cosmology is the Hubble parameter H(t). Therefore the proposed RG equation is [14]: In general µ 2 c can be associated to a linear combination of H 2 andḢ and the variety of terms appearing on the r.h.s. of (4) can be richer [15], but the canonical possibility is the previous one and hereafter we restrict to it. The coefficients a i , b i , c i . . . appearing in (4) are dimensionless and receive contributions from loop corrections of boson and fermion matter fields with different masses M i . It must be stressed that the general covariance of the action [10,12,13] necessitates the appearance of only even powers of the (cosmic-time t dependent) Hubble parameter H(t) on the right-hand-side of (4). For a specific framework where the above RG is concretely realized and the β-function coefficients can be computed, see [10].
We note at this stage that, if the evolution of the Universe is restricted to eras below the Grand Unified Theory (GUT) scale, then for all practical purposes it is at most the H 4 terms (those with dimensionless coefficients b i ) that can contribute significantly. The H 2 term is of course negligible at this point, and the higher powers of H n for n = 6, 8, .. are suppressed by the corresponding inverse powers of the heavy masses M i , which go to the denominator, as required by the decoupling theorem. In the scenarios of dynamical breaking of local supergravity discussed in ref. [6][7][8], the breaking and the associated inflationary scenarios could occur around the GUT scale, in agreement with the inflationary phenomenology suggested by the Planck satellite data [2], provided Jordan-frame supergravity models (with broken conformal symmetry) are used, in which the conformal frame function acquired, via appropriate dynamics, some non trivial vacuum expectation value. For these situations, therefore, corrections in (4) involving higher powers than H 4 will be ignored.
In the next sections, after revising the general framework of Starobinsky inflation, we shall compute ρ Λ in such supergravity models and study their evolution from the exit from the Starobinsky inflationary phase, that occurs in the massive gravitino phase, until today. The computation of ρ Λ will be made via the corresponding calculation of the one-loop effective action after massive gravitinos are integrated out in a path integral. Then, an identification of the effective equation of state can be derived by integrating (4), following the approach of [16][17][18]20]. Before doing so, it is instructive to review first the emergence of Starobinsky-type inflation.

III. GENERIC STAROBINSKY INFLATION
Starobinsky inflation is the oldest model of inflation [3], prior to the traditional, scalar-field-based, inflaton models. It is characterized for being able to realize the de Sitter (inflationary) phase from the gravitational field equations derived from a four-dimensional action that includes higher curvature terms, specifically of the type involving the quadratic curvature correction ∼ R 2 [3] 1 : Here κ 2 = 8πG = 1/M 2 Pl (in units of = c = 1 we are working on), G = 1/m 2 P is Newton's (gravitational) constant in four space-time dimensions, with m P the Planck mass, and M is a constant of mass dimension one, characteristic of the model. Notice that the curvature terms in the action are just the dimension-4 combination m 2 P R/16π + R 2 m 2 P /(6M 2 ). With this normalization, M gives the value of the so-called scalaron mass. The smaller is M in Planck mass units (i.e. the larger is the dimensionless parameter m 2 P /M 2 in front of R 2 ) the longer is the inflationary time (cf. Fig. 2 of [15]). Of course M cannot be much below the natural scale of inflation, and in fact it should be of the same order, i.e. M ∼ M X , where M X is some GUT scale below the Planck mass. Typically M X ∼ 10 16 GeV ∼ 10 −3 m P .
The most relevant feature of this model is that inflationary dynamics is driven by the purely gravitational sector, through the R 2 terms. From a microscopic point of view, these terms can be viewed as the result of quantum fluctuations (at one-loop level) of conformal (massless or high energy) matter fields of various spins, which have been integrated out in the relevant path integral in a curved background space-time [35]. The model in fact is to be understood in the context of QFT in curved spacetime. The quantum mechanics of this model, by means of tunneling of the Universe from a state of "nothing" to the inflationary phase of ref. [3] has been discussed in detail in [36]. The above considerations necessitate truncation to one-loop quantum order and to curvaturesquare (four-derivative) terms, which implies that there 1 Our metric signature is (−, +, +, +) and the definitions of the Ricci and Riemann curvature tensors are Rµν = R λ µλν and R λ µνρ = ∂ν Γ λ µρ − . . . , respectively, i.e. we follow the exact threesign conventions (+, +, +) of Misner-Thorn-Wheeler [34]. must be a region of validity for curvature invariants such that O R 2 /m 4 P ≪ 1. Recalling that R ∼ 12 H 2 I in the inflationary phase (where H I is the nearly constant Hubble rate in that phase), we observe that this is indeed a condition satisfied in phenomenologically realistic scenarios of inflation [1,2], for which the inflationary Hubble scale H I is typically constrained to obey (2) (Planck data [2]) or (3) (BICEP2 data [4]), which are at present essentially the same.
Although the inflation in this model is not driven by fundamental rolling scalar fields, nevertheless the model (5) (and for that matter, any other model where the Einstein-Hilbert space-time Lagrangian density is replaced by an arbitrary function f (R) of the scalar curvature) is conformally equivalent to that of an ordinary Einstein-gravity coupled to a scalar field with a potential that drives inflation [37]. To see this, one firstly linearises the R 2 terms in (5) by means of an auxiliary (Lagrangemultiplier) fieldα(x), before rescaling the metric by a conformal transformation and redefining the scalar field (so that the final theory acquires canonically-normalised Einstein and scalar-field terms): where again κ = √ 8πG. These steps may be understood schematically via where the arrows have the meaning that the corresponding actions appear in the appropriate path integrals. The ensuing effective potential V eff (ϕ) is given by: One can check that the mass of the scalaron, which can be seenn as the new gravitational degree of freedom that the conformal transformation was able to elucidate from the Starobinsky action, is indeed given by the parameter M = 8π/3β: Note that for ϕ = 0 one hasα = 0 and the two conformally equivalent metrics coincide at this point. The effective potential for the scalar d.o.f. that conformally replaces the effect of the R 2 term is plotted in Fig. 1. We observe that V (ϕ) is sufficiently flat for κϕ = ϕ/M P l ≫ 1 FIG. 1. The effective potential (8) of the collective scalar field ϕ that describes the one-loop quantum fluctuations of matter fields, leading to the higher-order scalar curvature corrections in the Starobinsky model for inflation (5). Notice that according to (8)  (i.e. for sufficiently large values of ϕ as compared to the reduced Planck scale) to produce phenomenologically acceptable inflation. Obviously the scalaron field ϕ is effectively playing the role of the inflaton in this context. The difference with the usual inflaton is that ϕ is not a new scalar d.o.f. imported from outside the gravitational action, but just an integral part of it, namely, it is just a gravitational d.o.f. that describes in an effective (and very convenient) way the ∼ R 2 term of (5). The Starobinsky model based on the action (5) indeed fits excellently the Planck data on inflation [2], and also the corresponding data from the joint BICEP2-Planck analysis [5]. Quantum-gravity corrections in the original Starobinsky model (5) have been considered recently in [38] from the point of view of an exact renormalisation-group analysis [39]. It was shown that the non-perturbative betafunctions for the 'running' of Newton's 'constant' G and the dimensionless inverse R 2 coupling κ 2 β −1 ∼ M 2 /M 2 P l in (5) imply an asymptotically safe Ultraviolet (UV) fixed point for the former (that is, G(k → ∞) → constant, for some 4-momentum cutoff scale k), in the spirit of Weinberg [40], and an attractive asymptotically-free (κ 2 β −1 (k → ∞) → 0) point for the latter. In this sense, the smallness of the (inverse) R 2 coupling, required for agreement with inflationary observables [2], is naturally ensured by the presence of the asymptotically free UV fixed point.
The agreement of the model of [3] with the Planck data triggered an enormous interest in the current literature, and indeed Starobinsky inflation has been revisited from various points of view, such as its connection with no-scale supergravity [41] and (super)conformal versions of supergravity and related areas [42]. In the latter works however the Starobinsky scalaron field is fundamental, arising from the appropriate scalar component of some chiral superfield that appears in the superpotentials of the model.
Although of great value, illuminating a strong connection between supergravity models and inflationary physics, and especially for explaining the low-scale of inflation compared to the Planck scale, these works contradict the original spirit of the Starobinsky model (5) where, as mentioned previously, the higher curvature corrections are viewed as arising from quantum fluctuations of matter fields in a curved space-time background such that inflation is driven by the pure gravity sector in the absence of fundamental scalars. On the other hand, the scenario of [7], in which a Starobinsky-type inflation arises in the massive gravitino phase of SUGRA models, after integrating out the massive degrees of freedom is in the same spirit of Starobinsky, and even better, in the sense that the model does not have to assume dominance of conformal matter during inflation.
We next proceed to summarize the construction of the one-loop effective action of the massless degrees of freedom after massive gravitino integration in this dynamically-broken SUGRA model with spontaneous breaking of global supersymmetry (SUSY) [6].

IV. STAROBINSKY-TYPE INFLATION IN DYNAMICALLY BROKEN SUGRA
Dynamical breaking of SUGRA, in the sense of the generation of a mass for the gravitino field ψ µ , whilst the gravitons remain massless, occurs in the model as a result of the four-gravitino interactions characterizing the SUGRA action, arising from the torsionful contributions of the spin connection, characteristic of local supersymmetric theories.
Our starting point is the N = 1 D = 4 (on-shell) action for 'minimal' Poincaré supergravity in the second order formalism [43]: where R(e) and D ν ψ ρ ≡ ∂ ν ψ ρ + 1 4 ω νab (e) γ ab ψ ρ are defined via the torsion-free connection; and, given the gauge condition γ · ψ = 0, arising from the fermionic torsion parts of the spin connection. Extending the action off-shell requires the addi-tion of auxiliary fields to balance the graviton and gravitino degrees of freedom. These fields however are nonpropagating and may only contribute through the development of scalar vacuum expectation values, which would ultimately be resummed into the cosmological constant. Making further use of the above gauge condition together with the Fierz identities (as detailed in [6]), we may write where the couplings λ S , λ PS and λ PV express the freedom we have to rewrite each quadrilinear in terms of the others via Fierz transformation. This freedom in turn leads to a known ambiguity in the context of (perturbative) mean field theory [44] and can only be resolved by a non-perturbative treatment. Specifically, we wish to linearise these four-fermion interactions via suitable auxiliary fields, e.g.
where the equivalence (at the level of the action) follows as a consequence of the subsequent Euler-Lagrange equation for the auxiliary scalar σ. Our task is then to look for a non-zero vacuum expectation value σ which would induce as an effective mass m 3/2 ∼ σ κ √ λ S for the gravitino. This is however complicated by the fact that our coupling λ S into this particular channel is, by virtue of Fierz transformations, ambiguous at a perturbative level and, as mentioned, in order to fix them a fully nonperturbative treatment of SUGRA-like models would be required, which are not currently at hand. Nevertheless, there is another way out [6,8] whereby the Fierz ambiguities may be absorbed by dilaton-expectation-value shifts in an extension of N = 1 SUGRA which incorporates local supersymmetry in the Jordan frame, enabled by an associated dilaton superfield [45]. The (logarithm of the) scalar component ϕ of the latter can be either a fundamental space-time scalar mode of the gravitational multiplet, i.e. the trace of the graviton (as happens, for instance, in supergravity models that appear in the lowenergy limit of string theories), or a composite scalar field constructed out of matter multiplets. In the latter case these could include the standard model fields and their superpartners that characterise the Next-to-Minimal Supersymmetric Standard Model [46], which can be consistently incorporated in such Jordan frame extensions of SUGRA.
Upon appropriate breaking of conformal symmetry, induced by specific dilaton potentials (which we do not discuss here), one may then assume that the dilaton field acquires a non-trivial vacuum expectation value φ = 0, thus absorbing any ambiguities in the value of the appropriate coefficient λ S induced by Fierz (12). One consequence of this is then that in the broken conformal symmetry phase, the resulting supergravity sector, upon passing (via appropriate field redefinitions) to the Einstein frame is described by an action of the form (10), but with the coupling of the gravitino four-fermion interaction terms being replaced bỹ while the Einstein term in the action carries the standard gravitational coupling 1/2κ 2 . For phenomenological reasons, associated with gravitino masses in the ballpark of Grand Unified Theory (GUT) scales, one must haveκ ≫ κ. This is assumed to be guaranteed by appropriate microscopic dilaton potentials that break the (super)conformal symmetry of the Jordan-frame SUGRA appropriately.
To induce the super-Higgs effect [47] we couple to the action (10) the Goldstino associated to global supersymmetry breaking via the addition of where λ is the Goldstino, √ f expresses the scale of global supersymmetry breaking, and . . . represents higher order terms which may be neglected in our weak-field expansion of the determinant. It is worth emphasising at this point the universality of (15); any model containing a Goldstino may be related to L λ via a non-linear transformation [48], and thus the generality of our approach is preserved.
Upon the aforementioned gauge choice for the gravitino field γ µ ψ µ = 0 , and an appropriate redefinition, one may eliminate any presence of the Goldstino field from the final effective action describing the dynamical breaking of local supersymmetry, except the cosmological constant term f 2 in (15), which serves as a reminder of the pertinent scale of supersymmetry breaking. The non-trivial energy scale this introduces, along with the disappearance (through field redefinitions) of the Goldstino field from the physical spectrum and the concomitant development of a gravitino mass, characterises the super-Higgs effect.
The linearisation of the four-gravitino terms (13), when combined with the f 2 term of the super-Higgs effect implies a tree-level cosmological constant which must be negative due to the incompatibility of supergravity with de Sitter vacua (notice that in our conventions both σ and f have dimension +2 in natural units).
The one-loop effective potential for the scalar gravitino condensate field σ(x) (with vacuum expectation value σ c ∝ ψ µ ψ µ ) has a double-well shape as a function of σ(x) which is symmetric about the origin (cf. fig. 2), as dictated by the fact that the sign of a fermion mass does not have physical significance. Dynamical generation of the gravitino mass occurs at the non-trivial minima corresponding to σ c = 0. The potential of the σ c field is also flat near the origin, and this has been identified in [8] with a first inflationary phase.
In [6] the one-loop effective potential was derived by first formulating the theory on a curved de Sitter background [50,51], with cosmological constant (one-loop induced) Λ > 0, not to be confused with the (negative) tree-level one Λ 0 (16), and then integrating out spin-2 (graviton) and spin 3/2 (gravitino) quantum fluctuations in a given class of gauges (physical ), before considering the flat limit Λ → 0 in a self-consistent way. The detailed analysis in [6], performed in the physical gauge, has demonstrated that the dynamically broken phase is then stable (in the sense of the effective action not being characterized by imaginary parts) provided the scale of the gravitino condensate is equal or below the scale of spontaneous breaking of global SUSY: which guarantees the aforementioned result on the necessity of the negative nature of the tree-level cosmological constant (16). The former result demonstrates the importance of the existence of global SUSY breaking scale for the stability of the phase where dynamical generation of gravitino masses occurs, which was not considered in the previous literature [49]. In super-conformal versions of SUGRA, e.g. those in ref. [45,46], phenomenologically realistic scales for f 2 and gravitino mass of order of the GUT scale, appear for appropriate values of the expectation value of the conformal factor. These imply inflationary scenarios in perfect agreement with the Planck data [2,8], on equal footing to the original Starobinsky model.
In [7] we considered an extension of the analysis of [6] to the case where the de Sitter parameter Λ is perturbatively small compared to m 2 P , but non zero, so that truncation of the series to order Λ 2 suffices. This is in the spirit of the original Starobinsky model [3], with the rôle of matter fulfilled by the now-massive gravitino field. Specifically, we were interested in the behavior of the effective potential near the non-trivial minimum, where σ ≃ σ c is a non-zero constant (cf Fig. 2). The one-loop effective potential, obtained by integrating out [50,51] gravitons and (massive) gravitino fields in the scalar channel (after appropriate euclideanisation), may be expressed as a power series in Λ: where S cl denotes the classical action with tree-level cosmological constant Λ 0 (to be contrasted with the one- Generic shape of the one-loop effective potential (in dimensionless unitsκ 4 V eff , whereκ is the conformally rescaled gravitational coupling, see the text) for the gravitino condensate field σ in dynamically broken (conformal) Supergravity models in the presence of a non-trivial de Sitter background with cosmological constant Λ > 0 [7]. The Starobinsky inflationary phase is associated with fluctuations of the condensate and gravitational field modes near the non-trivial minimum of the potential, where the condensate σc = 0, and the potential assumes the value Λ > 0, consistent with supersymmetry breaking. The dashed green lines denote "forbidden" areas of the condensate field values, violating the condition (17), for which imaginary parts appear in the effective potential, thereby destabilizing the broken symmetry phase.
loop cosmological constant Λ 2 ): with R denoting the fixed S 4 background we expand around ( R = 4Λ, and the 4-dimensional Euclidean Volume is 24π 2 /Λ 2 ), and the α's indicate the bosonic (graviton) and fermionic (gravitino) quantum corrections at each order in Λ. The leading order term in Λ is then the effective action found in [6] in the limit Λ → 0, with 2 The reader should notice that, upon the restriction (17) guaranteeing the absence of imaginary parts in the one-loop effective action, the tree-level cosmological constant (16) Λ 0 < 0, while the one-loop one Λ > 0, as appropriate for a de Sitter background. Thus, Λ 0 should not be confused with the current-epoch positive cosmological constantΛ 0 , which we introduce later on, in section V.2, when we discuss RVM (cf. (57)). where and indicate the leading (as Λ → 0) contributions to the effective potential from bosonic (graviton) and fermionic (gravitino) quantum fluctuations respectively, to oneloop order. Above, µ is a RG scale, associated with a short-distance proper time cutoff [6], not to be confused with the RG scale of the RVM µ c (t) (cf. Sect. II), which is such that the flow from Ultraviolet (UV) to Infrared (IR) corresponds to the direction of increasing µ, σ c denotes the gravitino scalar condensate σ c ∝κ ψ µ ψ µ at the non-trivial minimum of the one-loop effective potential (cf. Fig. 2),κ is the conformally-rescaled gravitational constant in the Jordan-frame SUGRA model of [45], defined in (14), corresponding to a non-trivial v.e.v. of the conformal ('dilaton') factor, φ = 0, assumed to be stabilized by means of an appropriate potential, leading to the breaking of the conformal symmetry. In the case of standard N = 1 SUGRA, φ = 0. The remaining (higher order in Λ) one-loop quantum corrections then, proportional to Λ and Λ 2 may be identified respectively with Einstein-Hilbert R-type and Starobinsky R 2 -type terms in an effective action of the form 3 where we have combined terms of order Λ 2 into curvature scalar square terms. For general backgrounds such terms would correspond to invariants of the form R µνρσ R µνρσ , R µν R µν and R 2 , which for a de Sitter background all combine to yield R 2 terms 4 . The coefficients α 1 and 3 The reader should recall at this stage that the sign of g in √ g and the overall minus sign in front of the right-hand-side of (24) is due to the Euclidean-signature formulation of the path integral and disappears upon analytic continuation back to the Minkowski space-time at the end of the computations, which is necessary in order to make contact with phenomenology/cosmology (see, e.g. ). This should be understood in what follows, and especially in the context of linking the SUGRA model with the RVM in section V. α 2 in (24) absorb the non-polynomial (logarithmic) in Λ contributions, so that we may then identify (24) with (18) via where we note that α 1 is dimensionless whereas α 2 has dimension of inverse mass squared. The coefficients α F,B i , i = 1, 2 in Eq. (25) can be computed using the results of [6], derived via an asymptotic expansion: and To identify the conditions for phenomenologically acceptable Starobinsky inflation around the non-trivial minima of the broken SUGRA phase of our model, we impose first the cancellation of the "classical" Einstein-Hilbert space term R by the "cosmological constant" term Λ 1 , i.e. that This condition should be understood as a necessary one characterizing our background in order to produce phenomenologically-acceptable Starobinsky inflation in the broken SUGRA phase following the first inflationary stage, as discussed in [8]. This may naturally be understood as a generalization of the relation R = 2Λ 1 = 0, imposed in [6] as a self-consistency condition for the dynamical generation of a gravitino mass in the flat (zero Λ) limit. From Eq. (28) it follows that the (positive) cosmological constant Λ > 0 satisfies the four-dimensional Einstein equations in the non-trivial minimum, and in fact coincides with the value of the one-loop effective potential of the gravitino condensate at this minimum. As we discussed in [6], this non-vanishing positive value of the effective potential is consistent with the generic features of dynamical breaking of supersymmetry [52]. In terms of the Starobinsky inflationary potential (8), the value Λ > 0 corresponds to the approximately constant value of this potential in the high ϕ-field regime (κϕ ≫ 1) of Fig. 1, in the flat region where Starobinsky-type inflation takes place. Thus we may set where H I the (approximately) constant Hubble scale during inflation, which is constrained by the current data to satisfy (2) or (3). In the SUGRA context under discussion H I is linked to the scale of global SUSY breaking through H I ≃ f /M P l .
The effective Newton's constant in (24), after the imposition of (28), is then defined as and from this, we can express the effective Starobinsky parameter (5) in terms of κ eff as This condition thus makes a direct link between the action (18) with a Starobinsky type action (5). Comparing with (5) we can determine the effective scalaron mass in this case: As we know, this mass parameter also sets the order of magnitude of the inflationary scale in the Starobinsky model. We may then determine the coefficients α 1 and α 2 in order to evaluate the scale 1/ √ β ∼ M of the effective Starobinsky potential given in Fig. 1 in this case, and thus the scale of the second inflationary phase.
In [7] we searched numerically for points in the parameter space such that: • The effective equations are satisfied, together with the condition (28), • the cosmological constant Λ is small and positive, satisfying (29), and for phenomenological reasons it should be of order to ensure the validity of our expansion in Λ, consistent with the phenomenology of Planck-satellite data [2].
• The scalaron mass should also be of order M ∼ 10 −5 M Pl , hence allowing us to achieve phenomenologically acceptable Starobinsky inflation in the massive gravitino phase, consistent with the Planck-satellite data [2]. Forκ = κ (i.e. for non-conformal supergravity), we were unable to find any solutions satisfying these constraints. This of course may not be surprising, given the previously demonstrated non-phenomenological suitability of this simple model [6]. If we considerκ ≫ κ however, we find that we are able to satisfy the above constraints for a range of values 5 .
In general, typical values obtained in phenomenological realistic conformal SUGRA models satisfyκ ≫ κ (e.g. of orderκ/κ = O(10 3 − 10 4 )), under the constraints (29) and (34), in such a way that: Since the scale of SUSY breaking must be in the ballpark of the typical GUT scale associated to the inflation, namely √ f ∼ 10 16 GeV∼ 10 −2 M P l , from the above we have Λ ∼ κ 2 f 2 = f 2 /M 2 P l ∼ 10 27 GeV 2 . As a result the scale of the gravitino is some two to three orders of magnitude below the GUT scale, that is to say, m 3/2 ∼ √ Λ 10 13 GeV∼ 10 −5 M P l . These values are compatible with both the combined Planck and Bicep2 bound (3) and the typical mass of the gravitino in this framework [7].
Exit from the inflationary phase is, of course, a complicated issue which we shall not discuss here at the level of the SUGRA model itself, aside from the observation 5 A comment concerning SUGRA models in the Jordan frame with such large values for their frame functions is in order here. In our approach, the dilaton 2φ could be a genuine (dimensionless) dilation scalar field arising in the gravitational multiplet of string theory, whose low-energy limit may be identified with some form of SUGRA action. In our normalization the string coupling would be gs ≡ e φ = κ/κ , which may not be far from values attained in realistic phenomenological string models. On the other hand, in the Jordan-frame SUGRA models of [45], the frame function reads Φ ≡ e −2φ = c. , in the notation of [46] for the various matter super fields of the next-to-minimal supersymmetric standard model that can be embedded in such supergravities. The quantity χ is a constant parameter. At energy scales much lower than GUT, it is expected that the various fields take on subplanckian values, in which case the frame function is almost one, and henceκ ≃ κ for such models today. To ensureκ ≫ κ, and thus large values of the frame function, Φ ≫ 1, as required in our analysis, one needs to invoke trasnplanckian values for some of the fields, H 0 u,d , and large values of χ, which may indeed characterize the inflationary phase of such theories. A similar situation occurs for the values of the Higgs field (playing the role of the inflaton) in the non-supersymmetric Higgs inflation models [53]. that it can be achieved by coherent oscillations of the gravitino condensate field around its minima and subsequent decays to radiation and matter fields (thus requiring detailed knowledge of the matter content of the SUGRA models in order to arrive at quantitative predictions for the exit phase), or tunnelling processesà la Vilenkin [36]. However, in the next section we will show that the SUGRA model can be represented by an effective running vacuum model along the lines indicated in Sect.II, and from this point of view the exiting from the inflationary phase into the standard radiation phase can be guaranteed on very general grounds.
Before doing so, though, we should make some important remarks concerning the presence of logarithms of the de Sitter scale Λ in the coefficients α i of the curvature terms of the effective action (24). When one computes the effective action in a fixed de Sitter background, it is tempting to identify a Λ term with the Ricci scalar, which eventually will be allowed to depend on time. Thus, naively, the presence of logarithms would imply non-polynomial terms of the form R lnR which would be problematic for any RVM interpretation of the exit from the inflationary phase, as it would contradict the spirit of the approach where only integer powers of the curvature terms would be allowed in the respective flow equations [10,11,[13][14][15]. Fortunately this is not the case. To understand this, we first remark that any effective action obtained by integrating out massive degrees of freedom such as gravitino fields, that we restrict ourselves here, must consists for reasons of covariance and consistency of the weak gravitational fluctuations about the de Sitter background only of polynomial structures of the curvature tensors, for instance to fourth order in derivatives terms involving the squares of the Ricci scalar and Ricci tensors and covariant derivatives thereof. Any R lnR term would be incompatible with the weak gravity perturbative expansion about a background, say of constant non-zero curvature.
Thus, the coefficients α 1 and α 2 in the action (24) are kept fixed, not undergoing temporal evolution, which is guaranteed by the fixing of the two free scales in the problem µ (35) and Λ (29). Notice that the scale µ should not be confused with the subsequent RG scale µ c (t) that describes the cosmological evolution of the RVM vacuum (cf. Sect. II). Indeed, the scale µ first of all is a high energy cut-off. As already mentioned, it plays the role of a proper-time cutoff scale [6], appearing in the integral representations of some ζ-functions that are part of the determinants arising in the path integral of the SUGRA action arising from integrating out massive spin 3/2 (gravitino) and spin 2 (graviton) fluctuations about the de Sitter background. The scale µ is therefore, in contrast to µ c (t), an inverse renormalization group scale. Its value has to be fixed so as to guarantee SUGRA breaking and to generate a fixed gravitino mass which should not depend on time. This implies that the spontaneous breaking of SUGRA and the inflationary phase are characterised by such fixed scales, which implies the time independence of Λ (or, equivalently, the Hubble parameter) during inflation, the gravitino mass, related to the gravitino condensate vacuum expectation value σ c , and thus the coefficients α 1 and α 2 . On the other hand, integer positive powers of Λ, appearing in the effective action may be replaced by higher order tensorial structures involving the square of the curvature tensors, which are allowed to vary with the cosmic time during the RVM phase after exit from inflation. Notice that microscopically the exit phase is characterised by an unknown sort of phase transition, either through decays of the gravitino condensates to matter parts and reheating of the inflated universe, or tunnelling, as mentioned previously, and thus using different RG running to relate various eras of the Universe after inflation is to be expected.

V. "DECAY" OF EFFECTIVE VACUUM ENERGY: RUNNING VACUUM MODEL
The main aim of this section is to demonstrate that there exists a family of time-dependent effective vacuum energy decaying models of running type, i.e. the class of the running vacuum models (RVM's) introduced in Sect. II, which characterize the evolution of the Universe from the exit of the Starobinsky inflationary phase till the present era. In fact, the RVM's are able to interpolate on very general grounds the primeval de Sitter epoch with the late time de Sitter era, i.e. the dark energy one, where a much smaller cosmological constant essentially dominates. We shall follow the approach of the RVM outlined in Sect. II, in which the vacuum energy density ρ Λ (H) varies with time through its dependence on H = H(t). The Hubble parameter, having dimension of energy in natural units, acts as the natural running scale via the RG equation Eq.(4). As mentioned in Sect. II, only the even powers of H can be involved in that equation, owing to the general covariance of the effective action. This is an important point to make possible a general QFT description of this RG approach and is essential for the connection with the SUGRA model under discussion.
It is evident that the expansion (4) quickly converges at low energies, where H is rather small -certainly much smaller than any particle mass. No other H 2n -term beyond H 2 (not even H 4 ) can contribute significantly on the r.h.s. of equation (4) at any stage of the cosmological history below the GUT scale M X ∼ 10 16 GeV, where presumably inflation occurs.
On the other hand, if we want to deal with the physics of inflation and in general to the very early states of the cosmic evolution, we have to keep at least the term H 4 , which in fact is the dominant term in the series (4) during the high energy regime. In contrast, the terms H 6 /M 2 i and above are less and less important because these higher and higher powers of H are suppressed by the inverse powers of the heavy fermion and boson masses in the GUT, as required by the Appelquist-Carazzone decoupling theorem. Therefore, the dominant part of the series (4) is expected to be naturally truncated at the H 4 term. Higher order terms should contain the bulk of the high energy contributions within Quantum Field Theory in curved spacetime, namely within a semi-classical description of gravity near but (possibly a few orders) below the Planck scale. Models of inflation based on higher order terms inspired by the RG framework exist since long in the literature (see [11]) as well as the unified inflationdark energy framework of [10]. For a more phenomenological treatment unrelated to the RG, see [54][55][56][57].

V.1. A distinct class of running vacuum models
Based on the above arguments, it is natural to consider the case in which the highest power of the Hubble rate in the RG Eq. (4) is H 4 . Integrating the RG equation provides the simplest realization of RVM that can describe inflation and the various stages of the FLRW regime: Here c 0 is an integration constant (with dimension +2 in natural units, i.e. energy squared) which can be fixed from the low energy data of the current universe [21,24]. On the other hand the dimensionless coefficients are given as follows: and α = 1 96π 2 At this point, we would like to make some comments which will hopefully make the reader appreciate the physical interpretation of the running vacuum scenario. The coefficient ν behaves as the reduced (dimensionless) beta-function for the RG running of ρ Λ at low energies, whereas α plays a similar role at high energies. Notice that the index i depends on whether bosons (B) or fermions (F ) dominate in the loop contributions. Of course, since the coefficients (ν, α) play the role of oneloop beta-functions (at the respective low and high energy scales) they are expected to be naturally small because M 2 i ≪ M 2 Pl for all the particles, even for the heavy fields of a typical GUT. Indeed, an estimate of ν within a generic GUT is found in the range |ν| = 10 −6 − 10 −3 [10]. The dimensionless coefficient α is also small, |α| ≪ 1, because the inflationary scale H I is certainly below the Planck scale, see Eq. (3). From the observational viewpoint, utilizing a joint likelihood analysis of the recent supernovae type Ia data, the CMB shift parameter, and the Baryonic Acoustic Oscillations, it has been found |ν| = O(10 −3 ) [21,22,24,25], which is nicely in accordance with the aforementioned theoretical expectations as well as it insures a mild dynamical behavior of the vacuum energy at low energies. As we have already stated in section II, the Quantum-gravity corrections in the Starobinsky model have been found in the context of the RG analysis [38,39] through the appropriate betafunctions. The fact that the nature of the main coefficients of both theories (running vacuum and Starobinsky) is based on the RG approach is a hint that perhaps there is a possible connection between the two models.
Indeed as we will confirm below, this is the case. In particular, let us start with the effective action (24), with coefficients (25), (26), (27) and the constraints (28), (35). As we have seen, this action was obtained after integrating out both quantum-gravity (metric) fluctuations and massive gravitino fields. The action admits dynamical solutions of de-Sitter vacua with small cosmological constant Λ of order less than the GUT scale at early epochs of the Universe. At the non-trivial minimum, the effective potential takes on a value of order Λ (cf. Fig. 2).
Around the minimum, we should replace Λ by an effective Hubble parameter during inflation, H I , Eq. (29).
At the end of the day the effective (dynamical) vacuum energy density, ρ Λ (H), during the inflationary phase of our SUGRA model can be extracted from the SUGRA effective action Γ (24), upon applying the constraint (28) and analytically continuing the results back to Minkowski space-time signature. In particular, the effective potential is defined as . Doing so, we observe that the so obtained ρ Λ (H), remarkably, adopts precisely the generic RVM structure (36) around that phase, in which the Ricci scalar -see Eq. (50) below -boils down to R ≃ 12H 2 since H remains (approximately) constant in this phase. Some important remarks are in order here. The imposition of the constraint (28) during the Starobinsky inflationary phase implies, as already mentioned, that the correct phenomenology is attained as a result of the effective gravitational coupling (30) that characterises that phase. So, if the constraint was an exact result, the effective vacuum energy density of the SUGRA model would then correspond to the R 2 → 144 H 4 terms in (24) with (30) playing the rôle of the effective gravitational constant, where we used (25), (26), (27). The form (39) constitutes an admissible class of RVM (cf. Eq. (36)). Notice that in Eq. (39) there is no ν term. This is important, in the sense that in such a model, as a result of the effective gravitational constant (30) entering the game, which in this scenario [7] is viewed as the 'physical' reduced Planck mass of order 10 18 GeV, the gravitino mass and global SUSY breaking scales, (35), when expressed in terms of κ eff are of order one, that is one encounters a Planckscale gravitino. Despite this, the vanishing of ν makes the renormalization-group equation (39) a consistent one within the perturbative class of (36).
However the above construction leads to the absence of a present-era (small, positive) cosmological constant c 0 . This arises from the fact that we imposed the constraint (28) exactly. It may well be that such a condition leaves (non-perturbatively, when all the higher than one-loop contributions are taken into account) a very small (constant in cosmic time) contribution c 0 > 0 which is preserved until the present day. Unfortunately, our one-loop construction does not allow us to explain the magnitude and the sign of this constant term, but this is equivalent to offering a solution to the cosmological constant problem, which of course our approximate one-loop analysis cannot provide. While we do not have a quantitative calculation at this point, the above argument provides at least an interesting qualitative explanation, to wit: the origin of the current cosmological term ρ 0 Λ in the model might well be attributed to quantum (non-perturbative) effects in the SUGRA effective action, which prevent the complete cancellation (28) from being realised. The constant residue c 0 is then transferred throughout the cosmic history and pops up in our days in the form of the tiny vacuum energy ρ 0 Λ = (3/κ 2 )(c 0 + νH 2 0 ) ≃ 3c 0 /κ 2 . Under this assumption, then, the initial gravitational coupling κ, and thus the Einstein term 1 2κ 2 d 4 x √ −g R would enter the game during the exit phase from inflation 6 . In such a case, in the exit phase, the effective vacuum energy of the SUGRA model at the inflationary phase should correspond to both α 1 and α 2 terms of (24), with the constraint (28) failing by a tiny amountc 0 > 0 corresponding to the present-era cosmological constant.
where the explicit form of the coefficients α i , i = 1, 2, given by Eqs. (25), (26) and (27), in which the scales µ and Λ are fixed through (29) and (35) respectively. As mentioned already, fixing of the scale µ and Λ, implies fixed values for the gravitino mass. It is worth stressing that the aforementioned ambiguity concerning the failure of the exact constraint (28) can be avoided altogether by observing that the corresponding ρ Λ (H) ultimately derives from the RG equation (4) discussed in Sect. II. It is therefore more appropriate, and elegant, if one performs the matching between the running vacuum energies ρ Λ (H) in the SUGRA model and RVM by equating the corresponding "RG beta" func- 6 The reader should bear in mind that, since during the inflationary phase the scalar degree of freedom of the Starobinsky action is slowly rolling, if there is inflation in the conformally rescaled metric (6), there is also inflation in the initial metric. The Starobinsky inflation arguments are also not affected if a small contribution to the cosmological constant, of order of the present-era one, enters the effective action (7), as this is negligible compared to the Hubble scale of inflation.
It is remarkable that the effective dynamical vacuum energy density ρ Λ (H) associated to the SUGRA model under consideration turns out to follow the general RG of the RVM, see Eq. (4), in which the coefficients of the H 2 and H 4 terms can be computed precisely from the underlying SUGRA framework. From Eq. (41) it follows that ρ Λ = ρ Λ (H) evolves ("runs") with the (time-evolving) value of H. Such evolution will be studied in more detail in the next section, but is relatively small in the beginning, namely the varying H is only slightly below the initial value given in Eq. (29). As a result the Universe can start with an initial inflationary phase, which is dominated by the ∼ H 4 term of (41). However, well after the inflationary period the ∼ H 2 -term takes over and remains in force until the present time, thereby providing a mild evolution of the current cosmological "constant".
Eventually one has to add to the effective action loop contributions from other matter fields, including particle multiplicities, but at the moment we take a mass of order of the gravitino mass and shall comment on the possible additional effects below. The value of m 3/2 is of order of the GUT scale, since it is proportional to the gravitino condensate through σ c ∼κ −1 m 3/2 , the latter being bounded from above by the GUT scale: (17)]. We must also keep in mind that for phenomenologically acceptable solutions of the broken SUGRA model the ratio r ≡κ/κ is forced to stay in the range r = O(10 3 − 10 4 ).
From the generic values (35) we adopt here we find: The integrated form of Eq.(42), yields of course the effective vacuum energy density at the scale H, ρ Λ (H), Eq. (40), within the current SUGRA scenario. We thus have: wherec 0 is the integration constant, which will play the rôle of the current-era cosmological constant, as already mentioned. The result naturally adopts the generic form of the canonical RVM, Eq. (36) with c 0 = κ 2c 0 /3 and the effective values for the coefficients ν and α given by Thus we observe that, within the context of the pure SUGRA model, where only the gravitino plays the rôle of "matter", both coefficients are small, of typical order 10 −9 , in accordance with their interpretation as βfunction coefficients of the running vacuum energy density. Let us also note that the above estimate for ν eff nicely fits with the formal expression obtained in the different context of anomaly-induced inflation, where it also takes the structure (37), namely a quantity proportional to the (squared) ratio of a heavy particle scale (in general a collection of them) to the Planck masssee [10] for details. In the full SUGRA case, after its coupling to ordinary matter and radiation fields, other SUSY heavy fermions with masses near the SUSY breaking scale should also contribute to ν eff and α, and this should enhance the value of this coefficient very significantly, thus bringing the obtained result even closer to the situation studied in Ref. [10]. Overall, such considerations in phenomenologically realistic SUGRA situations could bring these parameters to a range ∼ 10 −4 accessible to current observations [21,22]. We next remark that, in the general case where the parameters of the SUGRA model are varied from the generic values considered in (35), but within the allowed range, the values of ν eff and α eff can also undergo some variation and the sign of ν eff could change. However we stress that the sign of α eff remains always positive, which is essential for a correct description of inflation. This can be seen explicitly by comparing the various logarithms involved in the structure of the coefficient α 2 of the H 4 -term in Eq. (41), together with the size and sign of their respective numerical coefficients (cf. Eqs. (25), (26), (27)). The positivity of α 2 is maintained throughout the physically allowed parameter space, and derives essentially from the fact that H 2 I ∼κ 2 σ 2 c ∼ κ 2 m 2 3/2 , in agreement with the generic result (35).
Finally, let us note that the circumstance that ν eff could have either sign can only affects the dynamics of the vacuum energy in the late universe. The phenomenological implications for both signs have actually been explored recently in [21,22], see also [15].
The upshot of the above considerations points to the existence of a remarkable relation between the running vacuum model Eq. (36) with that of SUGRA Eq. (43). In the next section we discuss the predicted inflationary scenario [16,17] in the context of the general RVM [14,15], and provide some interesting phenomenology that can be tested for the low energy regime, namely for the current Universe.

V.2. Running Vacuum Evolution: from current to inflationary era
In this section we investigate the conditions under which the running vacuum model can provide an inflationary era. The point of this session is first to demonstrate that is, if one starts from an inflationary era, at an early epoch, obtained in the context of a microscopic model, such as the Starobinsky inflation induced in the SUGRA model, then the RVM can smoothly connect it with the current era, characterised by a very small value of the vacuum energy, with a cosmology of ΛCDM type. We shall follow a "bottom-up" approach, in which, by starting from a late epoch FLRW Universe and applying RVM evolution ("backwards" in cosmic time, or, in a RG sense, an IR to UV flow), one arrives at an inflationary era in the early Universe. As we shall see, however, in this bottom-up approach, there is no unique way to identify the underlying microscopic model during the de Sitter era, which was to be expected in view of the rather generic features encapsulated in the RVM evolution.
To this end, let us first reproduce the Friedmann equations in the framework of a running ρ Λ . The resulting equations are expected to be formally equivalent to the ΛCDM case, inasmuch as the Cosmological Principle, which is embedded in the FLRW metric, perfectly allows the possibility of a time-evolving cosmological term. In general, the Einstein-Hilbert action is given by (here and in what follows we are back in Minkowski-signature space-time, described by a metric g µν ): where in our case ρ Λ (t) = Λ(t)/κ 2 represents the effective vacuum energy, which is allowed to vary with the cosmic time (more specifically as a function of a dynamical cosmological variable that evolves with time), and L m is the Lagrangian of matter. Varying the action (46) with respect to the metric we arrive at where the totalT µν is given byT µν ≡ T µν − g µν ρ Λ , with T µν = −2∂L m /∂g µν +g µν L m the energy-momentum tensor corresponding to the matter Lagrangian. The extra piece is ρ Λ = Λ/κ 2 , that is to say, the vacuum energy density associated to the presence of Λ(t) (with pressure p Λ = −ρ Λ ). Let us remark that this equation of state (EoS) does not depend on whether the vacuum is dynamical or not. In contrast to other forms of dark energy, the vacuum is defined as that for which the EoS parameter is precisely ω = −1 in any circumstance. Modeling the expanding universe as a perfect fluid with velocity 4-vector field U µ , we obtain T µν = p m g µν +(ρ m + p m ) U µ U ν , where ρ m is the density of matter-radiation and p m = ω m ρ m is the corresponding pressure, in which ω m is the EoS of matter. Obviously,T µν takes the same form as T µν with ρ tot = ρ m + ρ Λ and In the context of a spatially flat FLRW metric, we derive the Friedmann equations in the presence of a dynamical Λ-term: and the Ricci scalar where the overdot denotes derivative with respect to cosmic time t. Note, that the Bianchi identities ▽ µT µν = 0, insure the covariance of the theory and, if the Newtonian coupling is strictly G =const., entail an energy exchange between vacuum and matter.
Combining equations (48), (49) and (51), we infer the basic differential equation that governs the dynamics of the Universe, namely the equation for the Hubble rate: Inserting in it the expression (36) for the dynamical vacuum energy we arrive at the following equation: The dynamics of this model has been thoroughly discussed in [16,17], see also [15]. We can summarize it as follows. First of all we identify the presence of an inflationary epoch (de Sitter phase) associated to the constant value solution H 2 = (1 − ν)H 2 I /α of Eq. (53), which is valid for the very early epoch of the universe (in which we can neglect c 0 /H 2 ≪ 1). In this regime, solving Eq.(53) we find where D is a positive constant of integration. For the early universe we assume that matter is essentially relativistic, thus we take ω m = 1/3 at this point. Overall, one can see from (54) that for Da 4(1−ν) ≪ 1 the universe starts from an unstable inflationary phase [early de Sitter era, H 2 = (1 − ν)H 2 I /α] powered by the huge value H I presumably connected to the scale of a Grand Unified Theory (GUT). Well after the primeval inflationary era, specifically for Da 4(1−ν) ≫ 1, the Universe enters the standard radiation phase. Subsequently the radiation component becomes subdominant and the matter dominated era appears. This is confirmed from the evolution of the vacuum energy and radiation energy densities. If we neglect ν and c 0 /H 2 in this early epoch, which is justified, we can insert (54) into (36) and we find: Then solving (51) we obtain: Here ρ I = 3H 2 I /κ 2 is the critical density in the inflationary epoch. As it is obvious from the above expressions, there is no singularity in the initial state: the Universe starts at a = 0 with a huge vacuum energy density ρ I /α (and zero radiation) which is progressively converted into relativistic matter. In the asymptotic radiation regime we indeed retrieve the standard behavior ρ r ∼ a −4 with essentially negligible vacuum energy density: ρ Λ ∼ a −8 ≪ ρ r . Graceful exit is, therefore, implemented.
Subsequently the radiation component becomes subdominant and the matter dominated era appears. This is the point when the c 0 /H 2 term in Eq. (53) surfaces and starts to dominate over αH 2 /H 2 I because the early de Sitter era is left well behind (H ≪ H I ). In this case Eq.(36) boils down to Λ(H) =Λ 0 + 3ν(H 2 − H 2 0 ) which corrects the concordance ΛCDM modelà posteriori. Notice, that is the vacuum (cosmological constant) energy density at the present time, which is positive, and should not be confused with the negative tree-level Λ 0 of the SUGRA model (19). This can be understood by studying the evolution of the universe at a time after recombination, therefore consisting of dust (ω m = 0) plus the running vacuum fluid with H ≪ H I . In this case, using d/dt = aH d/da, we can rewrite Eq.(53) as The solution satisfying the boundary condition H = H 0 at present (a = 1) is: Note that the aforementioned boundary condition fixes the value of the parameter c 0 as follows: c 0 = H 2 0 (Ω 0 Λ −ν). For ν = 0 we correctly recover the behavior of the ΛCDM. However, for small ν the Universe possesses a mildly evolving vacuum energy that could appear as dynamical dark energy without invoking spurious scalar fields. Furthermore, the above vacuum model is in agreement with the latest cosmological data and it predicts a growth rate of clustering which is in agreement with the observations (for more details concerning the late dynamics, see [21,22,24,25]).
Let us note that the main stage of the cosmic evolution where we can match the SUGRA model of Sect. IV with the RVM is the early period comprising inflation and the incipient radiation epoch, to which it leads after graceful exit, as described in this section. Later on the microscopic description of the SUGRA model is more difficult to analyze and we adopt here the point of view that the subsequent effective behavior of the Universe still follows the RVM flow dictated by the general RG equation (4), which, to order H 4 , entails the dynamical vacuum energy density (36). As previously mentioned, at low energies this implies that only the dynamical part ∼ H 2 is active and may lead to interesting phenomenological implications for the dynamical DE of the current universe [21,22,25])

V.3. Geometrical description: RVM versus Starobinsky
Finally, let us focus now on some aspects of the inflationary era that are especially relevant for the present study. As we have already mentioned, in this epoch we have a de Sitter solution H 2 ≃ (1−ν)H 2 I /α =const. Now, as previously indicated,Ḣ ≃ 0 in this period and hence R ≃ 12H 2 . Finally, neglecting the matter component from the action (46), which is justified in the inflationary period, and using Λ(H) ≃ 3αH 4 /H 2 I [see Eq.(36)] we schematically find This demonstrates our point that an inflationary vacuum can be connected smoothly, under the RVM, with a late epoch ΛCDM Universe. However, there is no unique way by means of which we can associate the inflationary era RVM effective action (58) to a microscopic model, which, as already mentioned, is to be expected due to the generic features of the RVM that describe classes of models and therefore may correspond to more than one microscopic theories, as far as the exit from inflationary phase is concerned. An interesting point concerns Eq. (58) if one replaces H 4 by the square of the Ricci scalar. In this case one may . Notice that, since α > 0 in our case, the RVM model is not formally and directly equivalent to a Starobinski-type model, for which the effective Lagrangian has the form (5) corresponding to a negative α coefficient in (58). This point has also been discussed in [15]. The root of the problem lies in the fact that the metric tensors of the two models, (5) and (58) are different, related by a non-trivial conformal transformation (6) involving the linearising Hubbard-Stratonovich field ϕ, which plays the rôle of the "physical" inflaton. The RVM metric is identified with the Einstein-frame metric g E µν in (6), while the original one-loop effective SUGRA action is described in terms of the g µν metric. Nevertheless, contact with Starobinskytype models, like the one induced within the context of SUGRA model examined here, can be achieved by observing that it is precisely the passage from the Einstein to Jordan-frame actions, via (6), which guarantees the opposite sign, relative to the Ricci scalar term, of the effective potential (8) of the Hubbard-Strstonovich inflaton field in (7). Upon making the identification for large In Fig. 3 we present the RVM effective potential ακ 2 U/H 2 I (solid curve) as a function of κϕ. In the same figure we plot the effective Starobinsky (dashed curve) potential which is showed in Fig. 1. From the comparison it becomes clear that, although the RVM and Starobinsky models live in different geometrical backgrounds, namely GR and R 2 , the two models are similar from the point of view of those features of inflation that can be described by an effective scalar-field dynamics. However, in other important aspects they are different. We should mention that the RVM model provides a simple description of the graceful exit and reheating problem, see [16,17,19,58] for details. As for the Starobinsky model, the reheating of the universe after the exit of the inflationary phase it has been discussed for example in [59]. In contrast to the Starobinsky model, a general effective action from where the RVM can be derived is not known [12], and currently this has been achieved only in some cases [10].

VI. CONCLUSIONS
In the light of the latest Planck+Bicep2 results [5] it has been proposed that the Starobinsky inflation plays a key role because it fits quite well the Cosmic Microwave Background (CMB) data on inflation. In the present paper we have further investigated the class of the running vacuum models (RVM) [14] (based on renormalizationgroup approach in curved spacetimes) and their implications on the inflationary universe [16,17]. In particular, we have addressed the possibility that they can mimic both the original Starobinsky model and the spontaneously broken SUGRA models based on dynamically induced gravitino condensates [8].
We have shown that the vacuum energy density ρ Λ (H) of these SUGRA models can be expressed as an even power series (4) of the Hubble parameter, which can be naturally truncated at the H 4 term. This is exactly the generic form expected in the simplest class of running vacuum models and therefore we can apply the known implications of these models for inflation [16,17]. Namely, after computing the modified form of the Friedmann equation, we find that the physics of inflation (which in our case occurs for H ≃ H I , a value associated to the spontaneously broken SUGRA model) is mainly described by the H 4 -term. Furthermore, being H 4 of order R 2 , we can trace some relationship of this model with Starobinsky inflation, although of course there is not a full identification or equivalence. Most noticeably we point out the distinguishing feature that within the entire class of running vacuum models -and hence, in particular, the SUGRA model that we have studied (which adapts to the same pattern) -the RVM performs successful graceful exit from the inflationary phase into the standard radiation regime [16,17]. This feature is characteristic of the running type of vacuum models, in contrast to the original Starobinsky model. Nevertheless, we have also shown here that the RVM model admits a scalar field description as well, via the vacuumon field, and its potential can be made equivalent to the Starobinsky potential upon appropriate scalar field redefinitions, despite the fact that the geometric backgrounds of the two models are very different. This dynamical equivalence implies that the two models should provide the same inflationary features, at least in all the aspects that can be described through an effective scalar field potential. Not so in other aspects which may differ from one model to the other. In particular let us emphasize that the Starobinsky model derives from a local effective action whereas the structure of the effective action in the general case RVM is not presently known, except in particular cases in which it is found to be non-local.
The low energy physics, on the other hand, and in particular the evolution of the Universe in the current epoch, is determined by the constant additive term of ρ Λ (H) and the power H 2 , which provides a remnant dynamical evolution still in our days, which is of the form ρ Λ (H) = ρ 0 Λ + (3 ν/κ 2 )(H 2 − H 2 0 ). Such evolution is mild because the coefficient of H 2 is small (it is interpreted as the β-function coefficient of the running vacuum energy at present). The signature of the RVM at present is precisely that mild quadratic dynamical behavior of the vacuum energy density around the current value ρ 0 Λ which is parameterized by the small parameter ν. The model has been thoroughly put to the test recently and it allows values of |ν| = O(10 −3 ) [21,22,24,25]). On the other hand, its successful performance in describing the physics of the early Universe (in particular the graceful exit of the inflationary phase into the standard radiation one) is also quite encouraging, especially after realizing that specific QFT models lead to this kind of behavior. In this paper we have shown that SUGRA models with a dynamically induced massive gravitino phase lead to the RVM behavior and therefore provide a strong support for a fundamental description of the cosmic history.
Finally, we would like to stress that, in the context of the running vacuum model, the universe evolution, and especially its accelerated phase either during inflation or at late times, is not attributed to an ad hoc scalar field, or to a modification of the gravitational interaction, but rather arises from the modification of the vacuum itself, which is endowed with a dynamical nature. Remarkably, the SUGRA framework studied here provides a concrete realization of this possibility within the fundamental context of quantum field theory in curved spacetime.