Thermodynamical aspects of running vacuum models

The thermal history of a large class of running vacuum models in which the effective cosmological term is described by a truncated power series of the Hubble rate, whose dominant term is $\Lambda (H) \propto H^{n+2}$, is discussed in detail. Specifically, by assuming that the ultra-relativistic particles produced by the vacuum decay emerge into space-time in such a way that its energy density $\rho_r \propto T^{4}$, the temperature evolution law and the increasing entropy function are analytically calculated. For the whole class of vacuum models explored here we findthat the primeval value of the comoving radiation entropy density (associated to effectively massless particles) starts from zero and evolves extremely fast until reaching a maximum near the end of the vacuum decay phase, where it saturates. The late time conservation of the radiation entropy during the adiabatic FRW phase also guarantees that the whole class of running vacuum models predicts thesame correct value of the present day entropy, $S_{0} \sim 10^{87-88}$ (in natural units), independently of the initial conditions. In addition, by assuming Gibbons-Hawking temperature as an initial condition, we find that the ratio between the late time and primordial vacuum energy densities is in agreement with naive estimates from quantum field theory, namely, $\rho_{\Lambda 0}/\rho_{\Lambda I} \sim10^{-123}$. Such results are independent on the power $n$ and suggests that the observed Universe may evolve smoothly between two extreme, unstable, nonsingular de Sitter phases.


Introduction
A non-singular early de Sitter phase driven by a decaying vacuum energy density was phenomenologically proposed long ago to solve some problems of the Big-Bang cosmology [1,2]. The basic idea is closely related to early attempts aimed at solving (or at least alleviating) some cosmic mysteries, such as the "graceful exit" problem, which plague many inflationary scenarios (for recent reviews see [3,4]), and also the CCP or cosmological constant problem [5,6,7,8,9,10,11], probably, the deepest conundrum of all inflationary theories describing the very early Universe.
Nevertheless, new theoretical developments are suggesting a possible way to circumvent such problems. Results based on the renormalization group (RG) theoretical techniques of quantum field theory (QFT) in curved spacetimes combined with some phenomenological insights provided a set of dynamical Λ(H)-models (or running vacuum cosmologies) described by an even power series of the Hubble rate [12,13,14] (cf. [11] for a review). In this line, we have discussed in a series of recent works a unified class of models accounting for a complete cosmological history evolving between two extreme (primeval and late time) de Sitter phases whose spacetime dynamics is supported by a dynamical decaying (or running) vacuum energy density [15,16,17,18,19,20,21]. In such models, the effective vacuum energy density is a truncated power series of the Hubble rate, whose dominant term is Λ(H) ∝ H n+2 , where for the sake of generality the power n > 0. Unlike several inflationary variants endowed with a preadiabatic phase, this decaying vacuum is responsible for an increasing entropy evolution since the very early Universe, described by the primeval nonsingular de Sitter spacetime.
Several theoretical and observational properties of this large class of nonsingular running vaccum scenarios have been discussed in the literature. In particular, Mimoso and Pavón shown their thermodynamic consistency based on the generalized second law of thermodynamics (GSLT) by taking into account quantum corrections to the Bekenstein-Hawking entropy [22]. Many details regarding the late-time dynamics can be found in Refs. [23,24,25,26,27], and especially in the recent, updated and very comprehensive works [28] and [29]. More recently, even the solution of the coincidence problem has been detailed discussed in the present framework [30], as well as in generic decaying vacuum cosmologies [31].
Finally, let us mention that the running vacuum models under study have recently been tested against the wealth of SNIa+BAO+H(z)+LSS+BBN+CMB data -see [32] for a short review -and they turn out to provide a quality fit that is significantly better than the ΛCDM. This fact has become most evident in the recent works [33,34], where it is shown that the significance of the fit improvement is at ∼ 4σ c.l. Therefore, there is plenty of motivation for further investigating these running vacuum models from different perspectives, with the hope of finding possible connections with fundamental aspects of the cosmic evolution In the present work, we focus our attention on the entropy of the cosmic microwave background radiation (CMBR) generated by this large class of non-singular decaying vacuum cosmology. By considering that the decaying vacuum process occurs under adiabatic conditions [in the sense that the specific entropy (per particle) is preserved] this means that the radiation produced satisfies the standard scaling laws, namely, n r ∝ T 3 and ρ r ∝ T 4 [36]. Under these conditions, the final value of the entropy produced by the decaying vacuum supporting the unstable primeval de Sitter phase is exactly the present radiation entropy existing within the current Hubble radius. Within this framework, the early decaying vacuum process is not plagued with "graceful exit" problem of most inflationary variants and generates the correct number, S 0 ≃ 10 88 [35], regardless of the power n present in the phenomenological decaying Λ (H)-term. In addition, the ratio between the primeval and the present day vacuum energy densities is ρ Λ0 /ρ ΛI ≃ 10 −123 , as required by some naive estimates from quantum field theory.
The article is structured in the following manner. In section 2 we justify the phenomenological decaying vacuum law adopted in the paper, whereas in section 3 we set up the basic set of equations and the transition from the early de Sitter to the radiation phase is addressed. How inflation ends and the temperature evolution law are presented in section 4, while in section 5 we discuss the entropy production generated by the decaying vacuum medium. Finally, the main conclusions are summarized in section 6.

General model for a complete cosmic history
The general Λ(H)-scenario accounting for a complete chronology of the universe (from de Sitter to de Sitter) is based on the following expression for the dynamical cosmological term defining the relevant class of running vacuum models under consideration [15,16,17]: Here H ≡ȧ/a is the Hubble rate, a = a(t) is the scale factor, and the over-dot denotes derivative with respect to the cosmic time t. By definition ρ Λ (H) = Λ(H)/8π G is the corresponding vacuum energy density (G being Newton's constant). The even powers of H (therefore n = 2, 4, ...) are thought to be of more fundamental origin due the general covariance of the effective action, as required by the QFT treatment in curved spacetime [12,13,14,11]. In the numerical analysis, however, we will explore also the cases n = 1, 3 and 4 for comparison.
The dimensionless free parameters α and ν have distinct status. The first can be absorbed (for each value of n in the arbitrary scale H I so that it can be fixed to unity without loss of generality (if the scale of inflation is not precisely known) [20,30], whereas ν has been determined from observations based on a joint likelihood analysis involving SNe Ia, Baryonic Acoustic Oscillations (BAO) and Cosmic Background Radiation (CMB) data, with the result |ν| ≡ O(10 −3 ) [23,24,25,26,27,28,29] -see especially the most recent analyses in which the ν = 0 result (associated to the ΛCDM in the post inflationary time) is excluded at ∼ 4σ c.l. [33,34]. The small value of ν is natural since at late times, the dynamical model of the vacuum energy cannot depart too much from ΛCDM. In this connection, by using a generic Grand Unified Theory (GUT), it has been shown that |ν| ∼ 10 −6 − 10 −3 [13]. Finally, the constant c 0 with the same dimension of Λ yields the dominant term at very low energies, when H approaches the measured value H 0 (from now on the index "0" denotes the present day values of the quantities).

Basic equations of the Λ(H) model
It is well known that the Einstein field equations (EFE) are valid either for a strictly constant Λ or a dynamical one [15,17]. Therefore, using the vacuum energy density ρ Λ = Λ/(8πG) and the nominal pressure law p Λ = −ρ Λ one can write the EFE in the framework of a spatially flat FLRW (Friedmann-Lemaître-Robertson-Walker) space-time where R is the Ricci scalar andT µν = (p m + ρ m )u µ u ν − (p m − ρ Λ )g µν is the total energy-momentum tensor and the index m refers to the dominant fluid component (nonrelativistic matter or radiation).
Obviously, in our case, the only difference with respect to the more conventional field equations is the fact that Λ = Λ(H). In this framework, the local energy-conservation law ∇ µT µν = 0 which insures the covariance of the theory reads: where we have used p m = ωρ m for the ordinary cosmic fluid, namely ω = 0 for dust and ω = 1/3 for radiation. In this enlarged framework, the Friedmann equations are given by By combining the above equations with the class of vacuum models (2.1) one obtains the equation driving the evolution of the Hubble parameter: A solution of this equation in the high energy regime [where the term c 0 /H 2 ≪ 1 of (3.5) can be neglected] is given by:  The quantityρ I in the above equations is the critical energy density defining the primeval de Sitter stageρ (3.10)

From initial de Sitter stage to radiation phase
For ω = 1/3 the energy density for the vacuum and radiation read: We can see from (3.1.1) that the value (3.10) just provides the vacuum energy density for a → 0, namely ρ Λ (0) =ρ I . As |ν| ≪ 1 we also see thatρ I /ρ I ∼ α −2/n , thereby effectively showing that the constant α can be absorbed in the scale H I , as remarked before. Let us also emphasize from the previous formula that for a → 0 we have ρ r /ρ Λ ∝ a 2n(1−ν) → 0, i.e. the very early universe is indeed vacuum dominated with a negligible amount of radiation. In the rest of the paper, we neglect the effects proportional to ν (which, since it is the coefficient of H 2 , is not essential for the study of the early universe, the epoch where the H 4 -term is fully dominant). Thus, without loss of generality, H I will be hereafter rescaled so that α = 1 and we set ν = 0in all the formulae. Within this framework, we haveH I = H I and Eq.(3.10) becomes In addition, from Eq. (3.6) it follows that the scale factor of the universe takes on an exponential form a(t) ∼ a i e H I t as long as the condition Da 3(1+ω) ≪ 1 is fulfilled. Obviously, this means that the universe is initially driven by a pure nonsingular de Sitter vacuum state, and therefore is inflating. However, the above mentioned de-Sitter inflationary phase is only ephemeral. Indeed it is easy to check that in the post-inflationary regime, i.e. for a ≫ a i (with Da 3n(1+ω) ≫ 1), we are led to H ∝ a −2 (or a ∝ t 1/2 ) for ω = 1/3 [see Eq. (3.6)]. Therefore, the present model evolves smoothly from inflation towards the conventional radiation stage, thereby insuring that the initial very large amount of vacuum energy density does not preclude the standard picture of the primordial big-bang nucleossynthesis.
On the other hand, using equation (3.5) one may check that the decelerating parameter, q = −ä/aH 2 = −1 −Ḣ/H 2 , for arbitrary values of the power index n, reads: Naturally, the existence of the radiation stage is not enough to identify precisely the end of inflation since the deceleration parameter varied from q = −1 (de Sitter) to q = 1 (radiation) and the inflationary period finished when q = 0, or equivalently,ä = 0. We will discuss this point in the next section.

When exactly Inflation Ends?
In order to answer this question we first combine Eqs. (3.1.1) and (3.1.2) so as to obtain the ratio of the radiation energy density (ω = 1/3) to the vacuum energy density: In principle, inflation must end when both components -the decreasing vacuum energy density and the created radiation energy density -contribute alike. Assuming that the scale factor at the point of vacuum-radiation "equality" is a = a eq the above expression implies that D a 2n eq = 1. This relation enables us to rewrite the Hubble parameter (3.6) in the following way: H(a, n) = H I 1 + (a/a eq ) 2n 1/n . It follows that the value of the Hubble function at the vacuum-radiation equality depends on the value of the parameter n and the initial scale H I : Now, by inserting this value of H eq into Eq. (3.1.4) we obtain effectively that q = 0 as should be expected. Hence, once the arbitrary scale H I is fixed, the energy scale or the moment for which the inflation ends is readily defined.
As remarked before, the start of the radiation phase in this model is not characterized by the canonical radiation value q = 1, as one might have naively expected. Due to the continuous energy exchange between vacuum and radiation there is a short period of time in which q goes from q = 0 to the standard result q = 1. In the begin of the standard adiabatic radiation regime the deceleration parameter is written as In the approach to the the radiation phase one may safely assume that it started when q(H rad ) ∼ 0.9999 with the Hubble parameter given by At this point it is appropriate to make the following comments concerning the cosmic history. First, in the full radiation era the value of the scale factor of the universe a rad becomes at least one order of magnitude larger than the corresponding value at the vacuum-radiation equality. Actually, the total entropy generated by the vacuum decaying process does not depend on the exact value of the ratio a rad /a eq . Second, when the radiation epoch is well left behind the Universe goes into the cold dark matter dominated era [Einstein-de Sitter, a(t) ∝ t 2/3 ]; and, after some billion years (∼ 7Gyrs) it enters the present vacuum dominated phase, in which Λ ≃Λ =const -confer [24,28] and [33,34].
In the upper panel of Fig.1 we provide the evolution of the normalized Hubble parameter E(a, n) = H(a, n)/H I for the following vacuum models n = 1 (dashed line), n = 2 (solid line), n = 3 (long dashed line) and n = 4 (dotted line) [see Eq.(2.1)]. We observe that for a ≪ a eq the cosmic evolution begins from an unstable inflationary phase [early de Sitter era, H ≃ H I ] powered by the huge value H I presumably connected to the scale of a Grand Unified Theory (GUT). Obviously, when the primeval inflationary era is left behind, particularly for a ≫ a eq , the cosmic evolution enters smoothly in the standard radiation period H ∝ a −2 . Overall, we would like to emphasize that the above natural mechanism for graceful exit is universal for the whole class of vacuum models which obey the restriction n > 0. Subsequently when the c 0 /H 2 quantity in Eq. (3.5) [23,24,26,27], and especially in the recent, updated and very comprehensive works [28] and [29].
Finally, by using n = 2 (corresponding to Λ(H) ∝ H 4 /H 2 I in the early universe) as a fiducial model in the large class (2.1), we can appreciate in the bottom panel of Fig.1 the relative deviation of the normalized Hubble parameters E(a, n) with respect to the n = 2 solution E(a, 2). Obviously, when we are far away from the epoch of the vacuum-radiation equality the deviations from the fiducial E(a, 2) case are extremely small. On the other hand there is a visible deviation from the latter around the transition region a → a eq . This deviation becomes at the level of ∼ −30% and ∼ +20% for n = 1 and n > 2 respectively.

Radiation and its Temperature Evolution Law
Assuming an "adiabatic" decaying vacuum it has been found that the specific entropy of the produced massless particles remains constant, despite the fact that the total entropy can be increasing 1 . This implies that the energy and the number density as a function of the temperature are given by the the standard expressions, namely: ρ r ∝ T 4 r and n r ∝ T 3 r , but the temperature does not obey the scaling relation T r (t) ∼ a(t) −1 . Such results were firstly derived based on a covariant nonequilibrium thermodynamic description [36,37]), and, more recently, using a kinetic theoretic approach [38]. In what follows, we discuss the temperature evolution law in the present framework along these lines.
Let us start with Eq.(3.1.2) which is rewritten in terms of a eq as follows ρ r = ρ I (a/a eq ) 2n where in the last equality we included the degrees of freedom (d.o.f.) of the created massless modes through the g * -factor (see [35]). Now, by solving for the temperature we find: 4n T eq (a/a eq ) n/2 where T eq = T r (a eq ). Obviously, T eq is the maximum value of the radiation temperature (4.2) which is defined by As it is expected T eq is given in terms of the arbitrary initial scale of the primeval de Sitter phase (ρ I , or equivalently H I ). Unlike the value of H eq (see Eq.(3.2.3)), this value of T eq is valid for the whole class of models since is does not depend on the power n. In particular, this unique maximum temperature suggests that the total entropy generated within the horizon and presently observed is basically the same for all models (see next section).
In the top panel of Fig. 2 we present the temperature evolution for several values of n. Note that in the very beginning of the evolution (when a → 0), the temperature of the created photons is also zero in accordance with the fact that ρ r = 0 (see Eq. (4.1)). However, for finite values of a = 0, we observe the existence of two regimes. In the first, a ≪ a eq , the radiation temperature (4.2) increases as T r ∝ a n/2 (it is linear for n = 2 [19]), reaching of course its maximum value at a = a eq . In this non-adiabatic regime, the vacuum instability guides constantly the model to the standard radiation epoch from a process that started in the non-singular de Sitter phase. Note that the evolution is different from inflationary models where a highly non-adiabatic "reheating" process happens immediately after the adiabatic evolution of the inflaton field (see [35] and references therein). In the opposite regime, a ≫ a eq , we are well within the radiation epoch when all running vacuum models decrease in the classical way, that is, following an adiabatic scaling law, T r ∝ a −1 . As we shall discuss bellow, the power-law increasing of the temperature up to the vacuum-radiation equality, beyond of which the universe enters in the standard temperature regime, is the reason of the large radiation entropy observed in the present epoch.
For the purpose of the present study it is also important to find a way to calculate the primeval parameters (ρ I , H I ) from first principles and indeed one may use the following approaches (from now on we consider natural units, = k B = c = 1). A first possibility is to consider that the initial de Sitter energy density ρ I is related with the Planck scale M P through ρ I = M 4 P . An alternative approach is to connect ρ I to the GUT energy scale M X ∼ 10 16 GeV and thus ρ I = M 4 X (see Refs. [20,21] It should be stressed that all the above characteristic scales are slightly below the corresponding Planck scale, which means that we are inside in the semi-classical QFT regime. For instance, by taking g * = 106.75 which corresponds to the particle content of the standard model of particle physics we find from Eq.(4.5) that T eq ≃ 0.26 M P which as expected does not depend on the power n. In this respect if we take into account the number of light d.o.f in the GUT then we find that the characteristic temperature is still smaller (nearly 10% of Planck mass). At this point it should be stressed that the primordial Gibbons-Hawking thermal bath was used only as a peculiar initial condition to fix the arbitrary scale H I . As we shall see in the next section, the ratio between the very early and late time vacuum energy densities depends only on the pair (H I , H F ) characterizing the extreme de Sitter phases. It has the expected magnitude thereby also contributing to alleviate the so-called cosmological constant problem in the context of the such models (see Eq. (5.7)).

The Radiation Entropy
The total entropy of the radiation included in the present Hubble radius (d H ≃ H −1 0 ) reads: In order to check the viability of our model we need to take into account that the total entropy should be measured from the initial entropy generated by the decaying vacuum. Since in our vacuum model the BBN proceeds fully standard, the equilibrium entropy formula remains valid because only the temperature law is modified [36,37,38]. Therefore, the radiation entropy per comoving volume is given by the well known expressions: where g s is the entropy factor at temperature T r (at very high temperature g s is essentially equal to the effective number of massless species, g * . However, for lower values there is a correction related to the freeze out of neutrinos and electron-positron annihilation).
With the help of equation (4.2) the comoving entropy (5.2) can be expressed as a function of the scale factor as follows: where r ≡ a/a eq and the function f n (r) depends on the parameter n as given below: The obtained result for the comoving entropy boils down to the one derived in Ref. [19] for n = 2, as it should. Note also that lim a→0 S r = 0, as it ought to be expected from the fact that the initial de Sitter state is supported by a pure vacuum (no radiation fluid). We also see that during the inflationary phase (a ≪ a eq ), that is, at the early stages of the evolution, the total comoving radiation entropy of the Universe increases very fast; in fact, proportional to a 3(1+n/2) . For instance, for n = 2 (corresponding to H 4 -driven inflation) the initial entropy raises as ∼ a 6 . Note also that for a = a eq , f n (r = 1) = 1/8 (n+2)/4n and the associated value S(a eq ) is still not the total comoving entropy that the decaying vacuum is able to generate (see discussion below Eq. (3.1.4)). This occurs only when a = a rad so that r 2n = (a rad /a eq ) 2n ≫ 1 and f n (r) ≃ 1 for all practical purposes. The evolution of the radiation temperature (normalized with respect to its maximum value) during the inflationary period for several values of the free parameter n (see caption of Fig.1 for definitions) and its transition into the FLRW radiation era. As in Fig.1, the lines correspond to the following scenarios n = 1 (dashed), n = 2 (solid), n = 3 (dot-dashed) and n = 4 (dotted). Bottom Panel: The evolution of the normalized comoving entropy from the inflationary period (where it increases) until reaching the saturation plateau for a/a rad ≥ 1. The asymptotic value corresponds the total entropy at present.
At this point the generated comoving entropy S r reaches the final value, S f rad = S r (a rad ), when the standard radiation phase begins.
It thus follows that the asymptotic (adiabatic) value of (5.3), defined by f (r) ≃ 1 for r ≫ 1, is given by: a saturated value that must be compared with the present day entropy since the subsequent evolution of the model is isentropic In the bottom panel of Fig. 2 we show the entropy as a function of the ratio r = a/a eq for several values of n. Notice, that the entropy is scaled to its value at the vacuum-radiation equality. Initially, for a ≪ a eq the amount of entropy differs among the vacuum models but when a → a eq the corresponding entropies start to converge and subsequently they reach a plateau, namely S r (a)/S r (a eq ) → 8 (n+2)/4n for a ≫ a eq that characterizes the standard adiabatic phase, which is sustained until the present days because the bulk of the vacuum energy Λ(H) ∝ H n+2 /H n I already decayed. Armed with the above expressions we now compute the prediction of the total entropy inside the current horizon volume ∼ H −3 0 . Using the temperature evolution law [see Eq. (4.2)] one may see that the expression T 3 eq a 3 eq which appears in the final entropy as given by (5.5) is equal to 2 −3(n+2)/4n T 3 rad a 3 rad , where T rad = T r (a rad ) and a rad was defined in Sect. 3.1 [see discussion below Eq. (3.2.4)]. Note that the n-dependence cancels out and we arrive at to the same final result (the power n is important only in the inflationary phase since it determines the time scale in which the radiation equilibrium phase is attained): where S 0 is given by (5.1). In the last step we used the entropy conservation law of the standard adiabatic radiation phase, which implies that g s T 3 rad a 3 rad = g s,0 T 3 r0 a 3 0 .
It should be stressed that in the very early de Sitter era, the radiation entropy is zero. However, it increases steadily as S r ∼ a 3(1+n/2) and, finally, as shown in the bottom panel Fig. 2, deep inside the radiation stage becomes constant and approaching its asymptotic present day observed value, S 0 .
In other words the running vacuum model provides an overall past evolution to the present ΛCDM cosmology by connecting smoothly between two extreme cosmic eras (early inflation and dark energy) driven by the vacuum medium [see Eq.(3.5)]. Specifically, using Eq. (4.4) we find that the ratio between the associated vacuum energy densities becomes in agreement with traditional estimates based on quantum field theory (see Refs. [5] and [7]). Note that the late time de Sitter scale, H F = H 2 0 Ω Λ 0 , was used in the above expression [30]. To conclude, the decaying vacuum model explains the amount of the radiation entropy and simultaneously it also alleviates the so-called cosmological constant problem.
Needless to say, for a final resolution of the problem one needs to understand the ultimate origin of the current value of the cosmological constant. Remarkably, the obtained description of the cosmic history is based on a unified dynamical Λ(H)-term accounting for both the vacuum energy of the early and of the current universe. An alternative approach to such description, based on the Grand Unified Theory framework, can be shown to provide similar results, see [20,21]. Somehow this shows that the obtained results are truly robust and independent of the initial conditions reigning in the primeval Universe.

Discussion and Conclusions
In this paper we have addressed a fundamental issue concerning the thermodynamics of the early Universe. It is well-known that within the context of the concordance cosmological model, or ΛCDM model, the thermal history is incomplete and it leads to inconsistencies with the present observations. One of the main problems is the well-known horizon problem, which can be rephrased thermodynamically as the entropy problem. The concordance model, in fact, cannot offer an explanation for the large entropy enclosed in our Hubble sphere, which is tantamount to say that it cannot explain the large amount of causally disconnected regions contained in it.
The well-known solution to these problems is inflation, a patch that has to be added to the incomplete ΛCDM description. While in the more traditional approach inflation is accomplished by postulating the existence of a new fundamental scalar field called inflaton, in the present work we have proposed an entirely different (but no less efficient) framework. It is based on the properties of a large class of non-singular decaying vacuum models whose structure is that of a truncated power series of the Hubble rate, Λ(H). The involved powers to describe inflation must be higher than H 2 since the latter can be relevant only for the current Universe.
In this work we have explored an inflationary dynamics where the decaying vacuum is triggered by an arbitrary higher power of the Hubble rate, H n+2 (n > 0) recently proposed [15,16,17].
In such unified model of the vacuum energy, the effective cosmological term is a dynamical quantity, which evolves extremely fast in the early Universe and goes through an approximate de Sitter phase in our time -the dark energy epoch. In principle, the late time Λ(H)-Universe as described by Eq. (2.1) remains very close to the concordance model, but it stills carries a mild vacuum evolution (hence a mildly evolving cosmological term) compatible with observations. In principle, such a term may act as a smoking gun of its lengthy and energetic history; indeed, a much richer history than that associated to the idle Λ-term inherent to the ΛCDM model. However, since the late time entropy production is very small, for the sake of simplicity, we have taken the ν parameter equal to zero. Therefore, the model discussed here can be seen as a primeval nonsingular phase of the standard ΛCDM model. Once we leave well back the early times, the νparameter recovers an important role which, in point of fact, makes the running vacuum models not only compatible with the current cosmological data but fully competitive with the ΛCDM description [33,34].
We have studied in detail some important thermodynamical aspects of that class of dynamical vacuum models. Most noticeably we have focused on the issue of the radiation entropy, its origin and generation in the first stages of the primeval Universe, and then its final impact on the current epoch. Our calculations were based on the assumption that the produced radiation from vacuum decay satisfy the standard relations, n r ∝ T 3 and ρ r ∝ T 4 [36], a hypothesis related with the idea of an "adiabatic" decaying vacuum and the fact that the specific entropy is preserved during the process [36]. The basic result is that at early times the temperature of the radiation increases (T r ∝ a n/2 ) until its maximum value determined by the equilibrium temperature T eq of the vacuum-radiation transition (see Fig. 2) and the same happens with radiation energy density. As a consequence, the entropy also increases at very early stages (S r ∝ a 3(n+2)/2 ) being later on conserved during the radiation epoch (neglecting the photon entropy produced in the electron-positron annihilation). In this context we have found that the large amount of radiation entropy now (S 0 ∼ 10 87−88 in natural units) can be fully accountable in our dynamical vacuum context. We have first shown that the entropy was produced during the inflationary process itself at the expense of the continuous vacuum decay. Subsequently, its production stagnated and this occurred shortly after the vacuum had lost its energetic power and the Universe entered the standard radiation phase. From this point onwards the adiabatic evolution of the cosmos carried the comoving entropy unscathed until our days. Overall, the wide class of running Λ(H)-vacuum models provides not only an alternative scenario for inflation (beyond the traditional inflationary scalar field models, some of them in serious trouble after the analysis of Planck 2015 data [40]), but also a new clue for graceful exit, which is indeed fully guaranteed within the Λ(H)-cosmology context and does not depend on the power n of the Hubble rate. The remarkable independence of both the graceful exit and the entropy prediction from the power n singles out such class of dynamical vacuum models from the rest.
Finally, as originally pointed out, the thermodynamical solution insures that no horizon problem exists because all points of the current Hubble sphere remained causally connected as of the early times when the huge entropy was generated by the decaying dynamics of the primeval vacuum. Interestingly enough, the above features are not only universal for the entire class of Λ(H)-models but also independent of the initial conditions of the early Universe. They provide a rather robust basis for the dynamical Λ(H)-cosmology, which is currently being tested and will continue being tested against the increasingly accurate observations.

Acknowledgments:
JASL is partially supported by the Brazilian Research Agencies, CAPES, CNPq and FAPESP (PROCAD2013, INCT-A and LLAMA, projects). JS has been supported in part by FPA2013-46570 (MICINN), CSD2007-00042 (CPAN), 2014-SGR-104 (Generalitat de Catalunya) and MDM-2014-0369 (ICCUB). SB also acknowledges support by the Research Center for Astronomy of the Academy of Athens in the context of the program "Tracing the Cosmic Acceleration". JASL and SB are also grateful to the Department ECM (Universitat de Barcelona) for the hospitality and support when this work was being finished.
A FRW radiation phase: Wien's law and the relation T rad /T I In this appendix a simple argument based on the equilibrium Wien's law is adopted to show that the decaying vacuum drives the model progressively to the radiation phase without reheating period (no exit problem).
To begin with we recall that the initial conditions in our picture are T r = 0 and T GH = T I = H I /2π. However, due to the evolution of the Universe, the Hubble parameter and the horizon temperature diminish while the temperature of the created radiation increases due to the continuous vacuum decay.
The model evolves out of equilibrium because the entropy generation is concomitant with the inflationary process. In principle, the standard radiation FRW phase is reached when the Wien law becomes strictly valid. In what follows we show that it provides an useful constraint on the value of the temperature in the begin of the FRW phase thereby suggesting the physical consistency of the model.
In natural units, the standard Wien's law (λ m T = 0.290 cm. K) reads: In order to equalize the horizon temperature, the wavelength at the maximum black body intensity (in the begin of the FRW phase, H = H rad , T = T rad ) is expected to be λ m < H −1 rad . Hence, the Wien's law takes the form: Now, in our model we know that H rad ∼ H I /(2.10 4 ) 1/n (see section 3.2) and H I = 2πT I . Thus, it follows that 1.27 < T rad /H rad < (2.10 4 ) 1/n T rad 2πT I , (A.3) which can be rewritten as: T rad T I > 7.98 (2.10 4 ) 1/n (A.4) The above relation shows two interesting aspects of the model: (i) T rad is much smaller than T eq = T I /2 1/n , but it can assume values not dramatically too low in comparison with the both characteristics scales (T I and T eq ); a result in agreement with the demonstrated progressive approach to the FRW phase (see Eq. (4.2) and the associated comments), and (ii) it also suggests that there is no exit problem in our model (the reheating process is not needed for this class of decaying vacuum models). For instance, by taking n= 1, 2 we see that T rad > 4.10 −4 T I , and T rad > 0.056T I , respectively. Note also that the inequality is still more safely satisfied for values of n ≤ 2. This is comprehensible because inflation ends faster for higher values of n.