Supersymmetry, Cosmological Constant and Inflation: Towards a fundamental cosmic picture via"running vacuum"

On the occasion of a century from the proposal of General relativity by Einstein, I attempt to tackle some open issues in modern cosmology, via a toy but non-trivial model. Specifically, I would like to link together: (i) the smallness of the cosmological constant today, (ii) the evolution of the universe from an inflationary era after the big-bang till now, and (iii) local supersymmetry in the gravitational sector (supergravity) with a broken spectrum at early eras, by making use of the concept of the"running vacuum"in the context of a simple toy model of four-dimensional N=1 supergravity. The model is characterised by dynamically broken local supersymmetry, induced by the formation of gravitino condensates in the early universe. As I will argue, there is a Starobinsky-type inflationary era characterising the broken supersymmetry phase in this model, which is compatible with the current cosmological data, provided a given constraint is satisfied among some tree-level parameters of the model and the renormalised cosmological constant of the de Sitter background used in the analysis. Applying the"running vacuum"concept, then, to the effective field theory at the exit of inflation, makes a smooth connection (in cosmic time) with the radiation dominance epoch and subsequently with the current era of the Universe, characterised by a small (but dominant) cosmological-constant contribution to the cosmic energy density. In this approach, the smallness of the cosmological constant today is attributed to the failure (due to quantum gravity non-perturbative effects) of the aforementioned constraint.


Introduction and Summary
The year 2015 celebrates a century since the proposal by Einstein for the theory of General Relativity (GR), an elegant and quite successful (from a phenomenological point of view) classical field theory of Gravitation, generalising Newtonian gravity. GR pioneered the idea of linking the geometry of the cosmos with its dynamics, by attributing the curvature of space time to mass/energy of matter/radiation. The theory has two versions, local -describing the dynamics of space-time in the neighbourhood of bodies, including celestial objects, and global, dealing with the entire Universe as a whole, the latter being described today by the celebrated Friedmann-Lemaitre-Robertson-Walker (FLRW) model for Cosmology. The FLRW model describes quite successfully the plethora of the currently available data on Cosmology in the framework of the so-called ΛC(old)D(ark)M(atter) a e-mail: Nikolaos.Mavromatos@kcl.ac.uk which is five orders of magnitudes smaller than the reduced Planck mass. This hierarchy between inflation scale and Planck scale, already prompted researchers to seek for a natural reason for it, and already in the early days after the inflation paradigm was at play, links of the low scale of inflation with supersymmetry have been made [8].
If supersymmetry is realised in nature however, it is certainly broken, with partners masses above TeV scale, given that they have not been discovered at the large hadron collider at present. Simple realisations of global supersymmetry (SUSY) breaking, such as in the Wess-Zumino model [9], can provide, when embedded in gravitational environments, slow-roll models for inflation consistent with the Planck data [6]. Rigorous embeddings of global SUSY to local supersymmetry (SUGRA) have also been considered and explored in the literature over the years in connection with various scenarios for inflation [10], such as hybrid [11], chaotic [12], no-scale SUGRA [13], Jordan-frame-SUGRA [14] and others [15]. In the case of some no-scale SUGRA models, the inflationary potential can be made similar to that of the so-called Starobinsky model of inflation [16], by a specific choice of parameters. Unlike the SUGRA case, however, in the original Starobinsky model, the inflaton scalar mode represents collectively the effects of the quadratic-order scalar curvature terms (R 2 ) in the effective gravitational action, obtained after integrating out (conformal) matter in the early hot Universe [16]. Starobinsky-type models seem to be preferred by the Planck data [6].
In this talk I shall review a rather minimal inflationary scenario which is associated indirectly with a Starobinsky type inflation in its original sense [16]. Specifically, I will discuss the appearance of scalar inflationary modes associated with higher curvature terms of the (one-loop) effective action of simple (3+1)-dimensional N = 1 SUGRA models without matter [17,18], obtained by integrating out massive gravitino degrees of freedom in the broken supergavity phase. The supersymmetry breaking occurs dynamically via the formation of gravitino condensates, as a consequence of the four-gravitino interactions that characterise (any) supergravity action, via the fermionic torsion parts of the spin connection. The approach is documented in a series of previous publications [19][20][21] and will be reviewed briefly here. We should stress that our approach does not supersymmetrise higher curvature terms in the action to obtain Starobinksy-like effective actions, as e.g. in the case of [22]. For us, the higher curvature terms emerge, as we mentioned, as a result of quantum loop corrections in the broken-local-supersymmetry phase of the N = 1 simple SUGRA action, after path-integrating out massive gravitino fields and graviton fluctuations in a de Sitter background (where such integration becomes meaningful [23]). The dynamical breaking process may be concretely realised by means of a phase transition from the supersymmetric phase where the bilinear ψ μ ψ μ representing the effective scalar degree of freedom has zero vacuum expectation value, to one where σ ≡ ψ μ ψ μ 0. The quantum excitations about this condensate vacuum are then identified with a gravitino condensate scalar field. Since this must be an energetically favourable process to occur, it then follows that the effective potential experienced by the gravitino condensate must be locally concave about the origin. The corresponding one-loop effective potential of the gravitino condensate scalar field, therefore, has the characteristic form of a Coleman-Weinberg double well potential, offering the possibility of hilltop-type inflation, with the condensate field playing the role of the inflaton [19,20], which would be the simplest scenario. However, in order to guarantee a slow-roll inflation one needs unnaturally large values of the gravitino-condensate wave function renormalisation, unless transplanckian scales for supersymmetry breaking are invoked. This prompted us to discuss a second [21], rather indirect way, by means of which the gravitino condensate field is associated with inflation. This is realised via the higher (in particular quadratic-order) curvature corrections of the (quantum) effective action of the gravitino condensate field, obtained after path-integration of graviton and gravitino degrees of freedom in the massive gravitino phase. These corrections induce a Starobinsky-type inflation [16], which occurs for quite natural values of the parameters of the N = 1 SUGRA model (actually, its Jordan-frame variants [14]), with the inflation scale and gravitino mass in the ball park of the GUT scale.
The structure of the talk is as follows: In Section 2 we review the formalism of the dynamical breaking of SUGRA in four dimensional N = 1 models, including superconformal extensions thereof (with broken conformal symmetry) that are necessitated for phenomenological reasons, as we shall see. In Section 3 we discuss scenarios for inflation of Starobinsky type that may occur in the massive gravitino phase, near the non-trivial minimum of the effective potential. In such scenarios, which are compatible with the Planck satellite data [6], the rôle of the inflaton field is played by the scalar mode that describes collectively the effects of scalar-curvature-squared terms that characterise the gravitational sector of the one-loop effective action in the broken SUGRA phase, after integrating out the massive gravitinos. In section 4 we apply [24] the running vacuum concept [25] in the effective field theory at the exit of the Starobinsky inflation, in order to ensure a smooth (in cosmic time) connection of that era with radiation-and matter-dominance eras of the Universe, and subsequently with the current era characterised by a small but dominant cosmological-constant type contribution to the cosmic energy density. Finally, conclusions and outlook are presented in section 5.

Dynamical breaking of N = 1 D = 4 SUGRA models
This section is based on material presented in refs. [20]. We consider the (on-shell) action for 'minimal' N = 1 D = 4 supergravity in the second order formalism [17,18]: where κ 2 = 8πG = 1/M 2 Pl , in units = c = 1, M Pl is the reduced Planck mass in four space-time dimensions, e is the vierbein determinant and the scalar space-time curvature R(e) and gravitational covariant derivative D ν ψ ρ ≡ ∂ ν ψ ρ + 1 4 ω νab (e) γ ab ψ ρ are defined via the torsion-free connection 1 . In the gauge condition γ μ ψ μ = 0 , we shall work with in this article, the torsion term L torsion , arising from the fermionic torsion parts of the spin connection, reads: We notice at this stage that the Fierz identities among the gravitinos, make the actual coefficients of these four-fermion terms ambiguous. This is a known ambiguity in the context of mean field theory, and only the non-perturbative physics can settle the value of the four-fermion terms [26]. Thus we may extend the definition of the torsion four-fermion interaction terms to include such ambiguities and 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. They satisfy the constraint: λ S − λ PS + 4λ PV = −3κ 2 /4. This freedom was addressed in the second reference in [20], where we refer the reader for details. 1 Our metric signature is (−, +, +, +) and the definitions of the Ricci and Riemann curvature tensors are R μν = R λ μλν and R λ μνρ = ∂ ν Γ λ μρ −. . . , respectively, with Γ λ μν = Γ λ νμ the torsion-free Christofel symbol, and the Ricci scalar is given by R = g μν R μν . Following the original ideas of dynamical symmetry breaking by Nambu and Jona-Lasinio [27], 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 serve as an effective mass for the gravitino. To induce the super-Higgs effect [28] we also couple in the Goldstino, associated to global supersymmetry breaking, via [3] L 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 (6); any model containing a Goldstino may be related to L λ via a non-linear transformation [29], and thus the generality of our approach is preserved. Upon the specific gauge choice (3) for the gravitino field 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 (6), 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 [28]. We may then identify in the broken phase an effective action where Λ is the renormalised cosmological constant, to be contrasted with the (negative) tree level cosmological constant and m dyn ∝ σ is a dynamically generated gravitino mass, the origin of which will be explained presently. It is worth stressing at this point that Λ 0 must be negative due to the incompatibility of supergravity with de Sitter vacua; if SUGRA is broken at tree level, then of course no further dynamical breaking may take place. For phenomenological reasons which have been outlined in detail in refs. [19,20], and we shall discuss below, we adopt an extension of N = 1 SUGRA which incorporates local supersymmetry in the Jordan frame, enabled by an associated dilaton superfield [14]. 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 low-energy 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, which can be consistently incorporated in such Jordan frame extensions of SUGRA [14,30], leading to supersymmetric generalisations of the Higgs inflation [31]. Upon appropriate breaking of conformal symmetry, induced by specific dilaton potentials (which we do not discuss here), one may assume that the dilaton field acquires a non-trivial vacuum expectation value ϕ 0. 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 DOI: 10.1051/ , 126 12602020 EPJ Web of Conferences epjconf/2016 02020 (2016) described by an action of the form (2), but with the coupling of the gravitino four-fermion interaction terms being replaced byκ while the Einstein term in the action carries the standard gravitational coupling 1/2κ 2 . In this way the v.e.v. of the dilaton may be combined with the Fierz-induced ambiguities (5), so that the nonperturbative physics that determines the latter also determines the scale of the conformal symmetry breaking. As in our approach we do not discuss details of the dilaton potential, we shall consider from now on such a v.e.v., as well the coefficient λ S of the four fravitino terms in (5) leading to the formation of scalar gravitino condensates, of interest to us here, as phenomenological. Expanding the graviton field about a de Sitter background [23] (under the assumption that it is a solution of the one-loop effective equations) with renormalised cosmological constant Λ > 0, and integrating out both bosonic and fermionic quantum fluctuations to one loop yields the following effective potential for the gravitino condensate field σ in the flat space-time limit Λ → 0, as detailed in ref. [20], where Λ 0 is given by (8) and indicate the contributions to the effective potential from bosonic and fermionic fields respectively, and μ is an inverse renormalisation group (RG) scale, associated with a proper time short distance cut-off [20], such that the flow from the ultraviolet to the infrared limit corresponds to the direction of an increasing μ. The one-loop effective potential (10) is depicted in fig. 1 and it has the characteristic double-well shape of dynamical breaking of a symmetry. We note that as we flow from UV to IR (i.e. in the direction of increasing μ), we obtain the correct double-well shape required for the super-Higgs effect, and secondly that tuning f allows us to shift V eff and thus attain the correct vacuum structure (i.e. non-trivial minima σ c such that V eff (σ c ) = 0). Moreover, the shape of the effective potential changes, as one varies the (renormalisation) scale μ from ultraviolet to infrared values (i.e. flowing in the direction of increasing μ), in such a way that the broken symmetry phase (double-well shaped potential) is reached in the IR. This indicates that the dynamical generation of a gravitino mass is actually an IR phenomenon, in accordance with the rather general features of dynamical mass in field theory. In the broken phase, the mass of the gravitino condensate is then given by m 2 σ ≡ V eff (σ c ), where σ c is the minimum of V eff and a prime denotes a functional derivative with respect to the gravitino-condensate field. As observed from (11), the bosonic contributions to the effective potential contain logarithmic terms which would contribute imaginary terms, leading to instabilities, unless From (10) it is straightforward to see that this condition is equivalent to the negativity of the treelevel cosmological constant Λ 0 , which is entirely sensible; if Λ 0 > 0 then SUGRA is broken at tree level (given the incompatibility of supersymmetry with de Sitter vacua) and there can be no dynamical breaking. As such, we must then tune f for a given value of μ to find self consistent minima σ c satisfying (12), thereby ensuring a real V eff . In fact, here lies the importance of the super-Higgs effect, and thus of a non-zero positive f 2 > σ 2 c > 0, in allowing dynamical breaking of local supersymmetry 2 .
As discussed in refs. [19,20], phenomenologically realistic situations, where one avoids transplanckian gravitino masses, for supersymmetry breaking scales f at most of order of the Grand Unification (GUT) scale 10 15−16 GeV, as expected from arguments related to the stability of the electroweak vacuum, can occur only for largeκ couplings, typically of orderκ ∼ 10 3 − 10 4 κ. Given the relation (9) this corresponds to dilaton vev of O (−10), where the negative sign may be familiar in the context of dilaton-influenced cosmological scenarios [33].
If we consider for concreteness the caseκ = 10 3 κ, which is a value dictated by the inflationary phenomenology of the model [19], we may find solutions with a vanishing one-loop effective potential at the non-trivial minima corresponding to:κ 2 σ c 3.5 ,κ 2 f 3.7 ,κ μ 4.0, which leads to a global supersymmetry breaking scale 2 In refs. [32], the importance of the super-Higgs effect was ignored, which led to the incorrect conclusion that imaginary parts exist necessarily in the one-loop effective potential, in the same class of gauges as in ref. [20] and here, and hence dynamical breaking of SUGRA was not possible. As we have seen above, such imaginary parts are absent when the condition (12) is satisfied, and thus dynamical breaking of SUGRA occurs. and dynamical gravitino mass At the non-trivial minima we findκ 4 V (1) 5.9×10 −13 , with tree-level cosmological constantκ 2 Λ 0 −1.4. We thus observe that fermion contributions to the effective potential are much stronger than the corresponding bosonic contributions for the cases of large couplingsκ κ. These values are phenomenologically realistic, thereby pointing towards the viability (from the point of view of producing realistic results of relevance to phenomenology) of the scenarios of dynamical breaking of local supersymmetry in conformal supergravity models 3 . On the other hand, in standard SUGRA scenarios, whereκ = κ, one finds, as already mentioned, transplanckian values for the dynamically generated gravitino mass [20]: m dyn 2.0 × 10 19 GeV, and a global supersymmetry breaking scale f 4.7 × 10 18 GeV, far too high to make phenomenological sense.

Starobinsky-type inflation in the broken SUGRA phase
Starobinsky inflation is a model for obtaining a de Sitter (inflationary) cosmological solution to gravitational equations arising from a (four space-time-dimensional) action that includes higher curvature terms. Specifically, an action of the type in which the quadratic curvature corrections consist only of scalar curvature terms [16] where M is a constant of mass dimension one, characteristic of the model, which cannot be determined by theoretical considerations, and hence should be considered phenomenological.
The important feature of this model is that inflationary dynamics is driven purely by the gravitational sector, through the R 2 terms, and that the scale of inflation is linked to M. From a microscopic point of view, the scalar curvature-squared terms in (15) are 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 [34].
The above considerations necessitate truncation to one-loop quantum order and to curvaturesquare (four-derivative) terms, which implies that there must be a region of validity for curvature invariants such that O R 2 /M 4 Pl 1. This is of course a condition satisfied in phenomenologically realistic scenarios of inflation [6,7], for which the inflationary Hubble scale H I satisfies (1) (the reader should recall that R = 12H 2 I in the inflationary phase). Although the inflation in this model is not driven by rolling scalar fields, nevertheless the model (15) (and for that matter, any other model where the Einstein-Hilbert space-time Lagrangian density 3 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 ϕ could be a genuine (dimensionless) dilation scalar field ϕ = 2φ arising in the gravitational multiplet of strings, whose low-energy limit may be identified with some form of SUGRA action. In our normalization the string coupling would be g s ≡ e φ =κ −1/2 . In such a case, a value ofκ = e − ϕ κ = O(10 3−4 ) would imply a large negative v.e.v. of the (four-dimensional) dilaton field of order φ = −O(5) < 0, and thus a weak string coupling squared g s = O(10 −2 ), which may not be far from values attained in phenomenologically realistic string cosmologies [33]. On the other hand, in the Jordanframe SUGRA models of [14], the frame function reads Φ ≡ e −ϕ = 1 − 1 c. , in the notation of [30] 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 [31]. 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 [35]. This can be seen by first linearising the R 2 terms in (15) by means of an auxiliary (Lagrange-multiplier) fieldα(x), and then 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: These steps may be understood schematically via where the arrows have the meaning that the corresponding actions appear in the appropriate path integrals, with the potential V(ρ) given by: This potential is sufficiently flat for large values of ρ (compared to the Planck scale) to produce phenomenologically acceptable inflation, with the scalar field ρ playing the role of the inflaton. In fact, the Starobinsky model, with its characteristically low value of the tensor-to-scalar ratio r, provides an excellent fit to the recent Planck data on inflation [6]. Below we consider an extension [21] of the dynamically broken SUGRA analysis of the previous section to the case where the de Sitter parameter Λ is perturbatively small compared to m 2 P , but not zero, so that truncation of the series to order Λ 2 suffices. This is in the spirit of the original Starobinsky model [16], with the rôle of matter fulfilled by the now-massive gravitino field. Specifically, we are interested in the behaviour of the effective potential near the non-trivial minimum, where σ c is a non-zero constant. It is important to notice at this point that, in contrast to the original Starobinsky model [16], where the crucial for inflation R 2 terms have been argued to arise from the conformal anomaly in the path integral of massless (conformal) matter in a de Sitter background, and thus their coefficient was arbitrary, in our 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 (self-consistently) 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 [21]. Moreover, in our analysis, unlike Starobinsky's original work, we will keep the contributions from both graviton (spin-two) and gravitino quantum fluctuations. Specifically, we are interested in the behaviour 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 gravitons and (massive) gravitino fields in the scalar channel (after appropriate euclideanisation), may be expressed as a power series in the renormalised cosmological constant Λ [21]: where S cl = − 1 2κ 2 Generic shape of the one-loop effective potential for the gravitino condensate field σ c in dynamically broken (conformal) Supergravity models in the presence of a non-trivial de Sitter background with cosmological constant Λ > 0 [21]. 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 (12), for which imaginary parts appear in the effective potential, thereby destabilising the broken symmetry phase. (19) indicate the bosonic (B) and fermionic (F) quantum corrections at a given order (i = 1, 2) in Λ, respectively. The leading order term in Λ is then the effective action found in [20] in the limit Λ → 0, The remaining quantum corrections , proportional to Λ and Λ 2 , may be identified respectively with Einstein-Hilbert R-type and Starobinsky R 2 -type terms in an effective action (21) of the form 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 α 2 absorb the non-polynomial (logarithmic) in Λ contributions, so that we may then identify (21) with (19) via 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 characterising our background in order to produce phenomenologically-acceptable Starobinsky inflation in the broken SUGRA phase following the first inflationary stage, as discussed in ref. [19]. This may naturally be understood as a generalisation of the relation R = 2Λ 1 = 0, imposed in ref. [20] as a self-consistency condition for the dynamical generation of a gravitino mass in a Minkowski space-time, discussed in the previous section. From (23) it follows that the cosmological constant Λ 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 [20], this non-vanishing positive value of the effective potential is consistent with the generic features of dynamical breaking of supersymmetry [36]. In terms of the Starobinsky inflationary potential (18), the value Λ > 0 corresponds to the approximately constant value of this potential in the high ϕ-field regime, where Starobinsky-type inflation takes place. Thus we may set Λ ∼ 3 H 2 I , where H I the (approximately) constant Hubble scale during inflation, which is constrained by the current data to satisfy (1).
The effective Newton's constant in (21) is then and from this, we can express the effective Starobinsky scale (15) in terms of κ eff as β eff ≡ α 2 /α 1 . This condition thus makes a direct link between the action (19) with a Starobinsky type action (15). Comparing with (15), we can then identify the Starobinsky inflationary scale in this case as The coefficients α 1 and α 2 (22) have been determined in ref. [20]: more independent invariant, R μν R μν , in addition to R 2 , since the Gauss-Bonnet combination, R μναβ R μναβ − 4 R μν R μν + R 2 , is a total space-time derivative. This is not the case for Jordan-frame SUGRA, where non-trivial-dilaton factors accompany the higher-order curvature terms. Nevertheless, in the minimum of the dilaton potential, where the conformal symmetry is broken, the dilaton equals its v.e.v. which is a space-time constant and, hence, this case is equivalent to the constant-dilaton one. where σ c denotes the gravitino scalar condensate at the non-trivial minimum of the one-loop effective potential (cf. fig. 2), andκ = e − ϕ κ is the conformally-rescaled gravitational constant in the model of ref. [14], defined previously via (9).
In ref. [21] we searched numerically for points in the parameter space such that the effective equations ∂Γ ∂Λ = 0 , ∂Γ ∂σ = 0, are satisfied; Λ is small and positive (0 < Λ < 10 −5 M 2 Pl , to ensure the validity of our expansion in Λ) and 10 −6 < M/M Pl < 10 −4 , to match with known phenomenology of Starobinsky inflation [6]. Forκ = κ (i.e. for non-conformal supergravity), we were unable to find any solutions satisfying these constraints, which was to be expected, given the previously demonstrated non-phenomenological suitability of this simple model [20]. If we considerκ >> κ, however, for instance of orderκ/κ = O(10 3 − 10 4 )), we find that we are able to satisfy the above constraints for a range of values [21], with gravitino mass and global SUSY breaking scale in the ball-park of GUT scale, and the Starobinsky scale of order M ∼ 10 −5 M Pl , which leads to phenomenologically acceptable inflation in the massive gravitino phase, consistent with the Planck-satellite data [6]. In particular, typical values obtained for the parameters of our conformal SUGRA models satisfy [20,24] Λ ∼ 3H 2 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 Pl , from (27) we obtain Λ ∼ κ 2 f 2 = f 2 /M 2 Pl ∼ 10 27 GeV 2 . As a result, the scale of the gravitino is some two to three orders of magnitude below the GUT scale, m 3/2 ∼ √ Λ 10 13 GeV∼ 10 −5 M Pl , which is compatible with the bound (1).

"Running Vacuum" and Broken SUGRA effective field theory
Exit from the inflationary phase is a complicated issue which we shall not discuss here, aside from the observation that it can be achieved by coherent oscillations of the gravitino condensate field around its minima. This is still an open issue, which may be addressed via construction of more detailed supersymmetric models, including coupling of the matter sector to gravity, which will determine the pattern of the inflaton decays.

The "Running Vacuum" Scenario and post inflationary evolution of the Universe
Nevertheless, non-trivial conclusions about the evolution of the Universe after the phase transition that describes graceful exit from inflation and a smooth connection (in the sense of an evolution in cosmic time) of this phase with the radiation-dominated and subsequent eras of the Universe till the current epoch, can be made by applying the concept of the "running vacuum" [25] (RV) immediately after the exit from inflation [24]. According to the RV hypothesis [25], the dynamical cosmic evolution ("decay") of the inflationary phase ground state ("vacuum") to the standard radiation regime can be described using a "renormalization-group-like" approach, whereby the time evolution of the vacuum energy density ρ Λ (t) 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 (RG) approach, and a natural candidate for such scale in FLRW cosmology is the Hubble parameter H(t). Therefore the proposed RG equation is [25]: The coefficients a i , b i , c i . . . appearing in (28) 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 where D is a positive constant of integration and a = a(t) is the scale factor of the universe. The universe will enter the standard radiation phase, with w = 1/3, in the case when Da 4(1−ν) 1, which is confirmed by substituting this solution (33) into (29) to obtain: where ρ I = 3H 2 I /κ 2 is the critical density in the inflationary epoch. Then solving (31) we obtain for the radiation energy density ρ m : From the above expressions, it is apparent that 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 from inflation is, therefore, naturally implemented in this formalism.
Subsequently, the universe will enter a matter dominated phase, and then the current era, in which the term c 0 in (29) becomes dominant. In this case Eq.(29) amounts tõ which plays the rôle of the cosmological constant at the present time [24], which is positive. The reader should notice the RV corrections to the cosmic-concordance ΛCDM model.

Applying the "Running Vacuum" scenario to the Dynamically Broken SUGRA model
Let us now apply the above reasoning to the effective potential of the broken N = 1 SUGRA model, at the exit of its inflationary phase [24]. The effective (dynamical) vacuum energy density, ρ Λ (H), during the inflationary phase of our SUGRA model can be extracted from the SUGRA effective action Γ (21), upon applying the constraint (23) and analytically continuing the results back to Minkowski spacetime signature. In particular, the effective potential is defined as . Moreover, at the exit from inflation, any logarithmic dependence of the coefficients α i on Λ, μ is kept fixed, and only positive integer powers of Λ are allowed to vary with cosmic time t and be set equal to 3 H(t) 2 . This is for covariance reasons explained in [24,25]. Doing so, we observe that the so obtained ρ Λ (H), remarkably, adopts precisely the generic RV structure (29) around that phase, in which the Ricci scalar is approximated by R 12H 2 , since the Hubble parameter H remains (approximately) constant in this phase.
The imposition of the constraint (23) during the Starobinsky inflationary phase implies, as discussed previously, that the correct phenomenology is attained. So, if the constraint was an exact result, the effective vacuum energy density of the SUGRA model would correspond to the R 2 → 144 H 2 terms in (21) with (24) playing the rôle of the effective gravitational constant, where we used (22), (26). The form (37) constitutes an admissible class of RV models (cf. Eq. (29)). Notice that in Eq. (37) there is no ν term. This is important, in the sense that in such a model, as a result of the effective gravitational constant (24) entering the game, which in this scenario [21] is viewed as the 'physical' reduced Planck mass of order 10 18 GeV, the gravitino mass and global SUSY breaking scales, (27), when expressed in terms of κ eff are of order one, that is one encounters a Planck-scale gravitino. Despite this, the vanishing of ν makes the renormalization-group equation (37) a consistent one within the perturbative class of (29).
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 (23) 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, nevertheless, the above argument provides at least an interesting qualitative explanation for it: 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 (23) 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 contribution (36), as discussed above.
Under this assumption, 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 5 . 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 (21), with the constraint (23) failing by a tiny amountc 0 > 0 corresponding to the present-era cosmological constant: where we have used the explicit form of the coefficients α i , i = 1, 2 (Eqs. (22), (26)), with the scales μ and Λ fixed through (27). Comparing (38) with (29) with c 0 =c 0 /3, one obtains the effective values for the coefficients ν and α in the matched-SUGRA/RV model [24]: 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. In the general case where the parameters of the SUGRA model are varied from the generic values considered in (27), 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 [24]. 5 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 (16), 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 (17), as this is negligible compared to the Hubble scale of inflation. Having matched the evolution of the SUGRA model with the RV flow at the exit of the inflationary Starobinsky phase, one can apply the analysis outlined in the previous subsection 4.1, in order to smoothly connect the inflationary and the current (cosmological-constant dominated) eras of this Universe, with the interpolation of radiation and matter dominated epochs. In this approach, the smallness of the cosmological constant today is attributed to a failure of the constraint (23) due to non-perturbative quantum (supergravity) effects. This offers an interesting new insight into the cosmological constant problem, which, needless to say however, remains unsolved.

Conclusions
In this talk we considered a minimal inflationary scenario, by means of which a gravitino condensate in supergravity models is held responsible for breaking local supersymmetry dynamically and inducing inflation in an indirect way by means of a Starobinsky-type inflation in the massive gravitino phase. The inflaton field in this approach is associated with the scalar mode that collectively parametrizes the effects of the quadratic-curvature contributions to the one-loop quantum effective action of the gravitino condensate, after integrating out graviton and massive gravitino degrees of freedom, in a de Sitter background. The model involves parameters that assume values of a natural and phenomenologically relevant order of magnitude, specifically global supersymmetry scale and gravitino masses of the order of GUT mass scales or less. Such a scenario is a truly minimal scenario for natural inflation, in the sense that it involves two scalar primordial composite modes to achieve dynamical breaking of a gauge symmetry (supergravity) and inflation.
Moreover, the effective potential of this model has been cast in a form that allows a running vacuum scenario to be in operation, thereby implying a smooth connection of the inflationary phase to the current era, characterised by a small cosmological constant. The smallness of the latter has been attributed to non perturbative quantum (supergravity) effects that lead to the failure of the constraint (23) characterising the Starobinsky inflationary phase of the model. From our analysis it becomes clear that, in order to ensure phenomenologically relevant supersymmetry breaking scales and gravitino masses, one needs to apply the above ideas to Jordan-frame extensions of the minimal N = 1 D = 4 SUGRA model, which involve a third scalar field -the dilaton. The latter is responsible, through its appropriate potential, for the breaking of the local scale symmetry of the theory. In the context of the next to minimal supersymmetric standard model, which Jordan-frame SUGRA models can incorporate, such dilatons may be composite of appropriate matter superfields, involving Higgs (supermultiplets).
Details of the microscopic matter model are important in order to ensure the correct cosmological evolution, in particular satisfaction of the Big-Bang-Nucleosynthesis constraints. A GUT scale gravitino can be made to decay fast enough so as not to disturb the BBN, but this depends on the details of the matter sector of the theory. Such topics are left for future investigations. Nevertheless, we believe that the above simple models for dynamical breaking of supergravity and their links with Starobinsky inflation and Running Vacuum models, deserve further study, and stand a serious chance of leading to realistic phenomenological scenarios compatible with the cosmological and particle physics data.