Constraints on Metastable Dark Energy Decaying into Dark Matter

. We revisit the proposal that an energy transfer from dark energy into dark matter can be described in field theory by a first order phase transition. We analyze the model proposed in Ref. [1], using updated constraints on the decay time of a metastable dark energy from the work of Ref. [2]. The results of our analysis show no prospects for potentially observable signals that could distinguish this scenario from the ΛCDM. We analyze, for the first time, the process of bubble nucleation in this model, showing that such model would not drive a complete transition to a dark matter dominated phase even in a distant future. Nevertheless, the model is not excluded by the latest data and we confirm that the mass of the dark matter particle that would result from such a process corresponds to an axion-like particle, which is currently one of the best motivated dark matter candidates. We argue that extensions to this model, possibly with additional couplings, still deserve further attention as it could provide an interesting and viable description for an interacting dark sector scenario based in a single scalar field.

We revisit the proposal that an energy transfer from dark energy into dark matter can be described in field theory by a first order phase transition.We analyze a metastable dark energy model proposed in the literature, using updated constraints on the decay time of a metastable dark energy from recent data.The results of our analysis show no prospects for potentially observable signals that could distinguish this scenario from the ΛCDM.We analyze, for the first time, the process of bubble nucleation in this model, showing that such model would not drive a complete transition to a dark matter dominated phase even in a distant future.Nevertheless, the model is not excluded by the latest data and we confirm that the mass of the dark matter particle that would result from such a process corresponds to the mass of an axion-like particle, which is currently one of the best motivated dark matter candidates.We argue that extensions to this model, possibly with additional couplings, still deserve further attention as it could provide an interesting and viable description for an interacting dark sector scenario based in a single scalar field.

I. Introduction 2 II. A Model for Dark Energy Decay 4
III. Dark matter from first order phase transition 5 IV. Bubble Nucleation 8

V. Conclusion and Prospects 10
A. Gravitational Effects 11

I. INTRODUCTION
One of the main goals of physical cosmology has been to understand the nature of the constituents of our Universe, especially dark matter (DM) and dark energy (DE) which, together, constitute nearly 96% of the total density of the Universe.Many theoretical models have been developed to explain the nature of the dark Universe [1].The DM component, due to its clustering properties, has been investigated in the context of several astrophysical and cosmological experiments.Moreover, unlike the case of DE, it has also been studied at particle physics level, being sought in direct detectors on Earth.In this context, one of the most interesting possibilities is that the DM can be an axion-like particle (a good review on axion-like dark matter can be found in [2], for instance), which is one of the main candidates for this component today 1 .On the other hand, the nature of the DE component remains yet very obscure, although there have been important advances in modeling its behavior beyond the simple cosmological constant scenario.
The explanation for a late-time accelerated phase in the Universe remains a topic of much debate, which is related to the well-known cosmological constant problem [12,13].In a different framework, there are many challenges when trying to embed such models in a more fundamental quantum gravity proposal.As an example, we can mention the Swampland conjecture, which describes a whole inhabitable landscape of field theories that are inconsistent with string theory, including the stable de Sitter vacua [14][15][16].In the process of trying to understand the DE properties, we are faced with the recurring discussion in the literature concerning the issue of the stability of a de Sitter phase.In the context of the late time Universe, a stable de Sitter phase has shown either to be hard to achieve from fundamental physics or even to be inconsistent in different theoretical contexts.For example, it has been subject of a long debate whether a de Sitter space is unstable due to infrared (IR) effects, as conjectured in [17][18][19][20][21], for instance.An instability of a de Sitter phase has also been obtained as a consequence of the backreaction effects of super-Hubble modes.As shown in several works [22][23][24][25][26][27][28], the backreaction of super-Hubble modes could give a negative contribution to the effective cosmological constant, causing the latter to relax.
From the observational point of view, one motivation for investigating possibilities beyond the cosmological constant solution is the current tension in H 0 measurements, which has shown to be alleviated in some quintessence models, as well as in models with dark sector interaction [29][30][31][32].Concerning the latter, in the framework of field theory, it can be argued that it is rather natural to consider an interaction between DM and DE, given that they are fundamental fields of the theory.Interacting models can also be claimed as a proposal for alleviating the coincidence problem [33], while also being able to provide a good fit to current data [19,[34][35][36][37][38][39][40].
In a related framework, when going beyond the cosmological constant scenario, it is natural to investigate the possibility of a metastable DE.As pointed out in [41], the remarkable qualitative similarity between the properties of the present DE and the component that supposedly drove inflation in the very early Universe makes it rather natural to put forward the hypothesis that the current DE can also be metastable [41][42][43][44][45][46][47][48][49][50].In this context, in [51], and more recently in [41], metastable DE phenomenological models were analyzed, in which the DE decay rate does not depend on external parameters, being assumed to be a constant depending only on DE intrinsic properties.In the latter analysis, they considered data from Pantheon compilation [52] in combination with BAO data from 6dFGS [53], MGS [54], BOSS DR12 [55], eBOSS DR14 [56], Lyα [57], and CMB [58,59], and found that the typical decay time in these scenarios must be many times larger than the age of the Universe.
One possible theoretical model that can provide a field theory description for the class of phenomenological scenarios considered in the analysis of [41] is the model proposed in [48], hereafter referred to as the MDE (Metastable Dark Energy) model.In this MDE model, a positive "cosmological constant" is modeled by a nonzero scalar vacuum energy with a potential of the expected order V ≈ 10 −47 GeV4 in the false vacuum.The potential of the scalar field in this model has a doubly degenerate energy minimum with small symmetry breaking terms that provide such a small energy difference2 .In this scenario, the field at the false vacuum represents DE.After the field passes the potential barrier, decaying from the false to the true vacuum, its equation of state is no longer that of dark energy, as it acquires non-negligible kinetic energy 3 .Analogously to what happens in the old inflationary scenario, the transition to the true vacuum occurs through the formation of bubbles of a new vacuum [65,66].Through this process, the energy released in the conversion of the false vacuum into the true can produce a new component, which has the properties of DM.In the MDE model, there is a single scalar field describing the dark sector of the Universe.Such a field is not responsible for inflation at early times and it has no significant coupling to the standard model sector (or to the inflation), except through the gravitational interaction.In addition to the absence of significant coupling to the standard model particles, despite there being coupling with gravity, we are in the small field regime.In this regime, where the field value is much smaller than the Planck scale, quantum corrections to the mass from gravity are expected to be small.
In the previous work of [48], it was shown that the mass of the DM in this scenario would correspond to the mass of an axion-like particle for a decay time of DE on the order of the age of the Universe, as assumed in the work of [48].The association of the equation of state of this real scalar field with the equation of state of a DM particle can be justified by the fact that after the phase transition, the oscillations about the quadratic minimum of the potential become the main important aspect [67,68].Apart from providing a unified description of the dark sector, the fact that scalar field in the MDE model has the mass of an axion-like particle, which is considered to be among the main candidates for DM today [2-11, 69, 70], motivates us to further explore this model in the light of the new constraints [41].
It is important to emphasize that the field after the phase transition is not the QCD axion.Following this transition, the oscillations around the quadratic minimum of the potential become the primary factor determining its effective equation of state.Consequently, the field behaves in a manner consistent with a dark matter equation of state, with a mass comparable to that of an axion-like particle.
However, despite this effective behavior 4 , the field may not possess all the characteristics of a standard axion dark matter candidate.In fact, it is debatable whether it can be classified within the broader category of axion-like particles, which encompasses a wide range of candidates [2].Nonetheless, due to similarities in mass and equation of state, leading to similarities in the cosmological behavior, we will refer to this "axion mass-like field" as an "axion-like field" for simplicity.
Here, we revisit the MDE model using the new constraints from [41] in order to test if the model is still viable and whether the newly constrained decay time still predicts the same mass for the DM particle produced in this scenario.As discussed above, due to the similarities of the current DE and the primordial inflation, one of our main goals is also to check if this model inherits the same problems of old inflation [71,72], i.e, potentially observable inhomogeneities from the process of bubble nucleation and evolution.
In particular, we aim to address the following questions: 1. Is MDE still a viable model considering the latest cosmological constraints?
2. What is the mass of the DM resulting from this process?
3. Would the bubble nucleation process in this model lead to inhomogeneities that could plague the model?
4. Could this model leave observational imprints that could be searched for in future experiments?
In order to address these questions, we organize this paper as follows: In Section II, we present the MDE model proposed in [48] to describe an interacting dark sector model based on field theory.In Section III, we compute the mass of the DM particle of this model that results from the DE decay with a characteristic time compatible with the recent constraints from [41].In Section IV, we analyze the bubble nucleation process and its subsequent dynamics in order to see if the model would predict observable inhomogeneities.In Section V, we conclude and discuss some future prospects.In Appendix A, we demonstrate the validity of the approximations considered.

II. A MODEL FOR DARK ENERGY DECAY
In this section, we present the MDE model proposed in [48] to describe an interacting dark energy model based on field theory.Since the model contemplates an energy exchange between dark energy and dark matter, none of these components is separately conserved.The equations describing the model are written as: where dot denotes the derivative with respect to cosmic time, ω i is the equation of state parameter (EoS), H is the Hubble parameter, Q i is the interaction term, and the subscript i indicates DE or DM.We consider that Q i is linearly proportional to the energy density of DE [48], as follows: In the model we are going to consider, the DE decay rate, Γ, does not depend on external parameters such as the curvature of the Universe or the scale factor.Instead, the DE decay rate is assumed to be a constant (similar to the case of the radioactive decay of unstable particles and nuclei).
Once the interaction term is defined, we can write an effective EoS for DE From Equation (2.4), we see that in a DE-DM interaction model, the effective EoS can have a phantom-like behavior (ω (eff) DE < −1) when the decay rate is negative (energy flow DM → DE), or a quintessence-like behavior (ω when the decay rate is positive (energy flow DE → DM), which is our focus.
In the MDE model, DE is represented by a scalar field with a potential energy endowed with doubly degenerate minima.When a small symmetry-breaking term is added, a small energy difference between the degenerate minima emerges, which we will denote by ϵ.The DE responsible for the current expansion of the Universe is represented by the scalar field at the metastable minima, where the potential energy is adjusted to ϵ ∼ 10 −47 GeV 4 .Such a potential can be written as where φ is a real scalar field of mass m and coupling constant λ and Q(φ) is the symmetry-breaking term, which is adjusted so that we have the value of cosmological constant at the metastable minimum 5 .Of course, this is a finetuned choice.This is the same fine tuning that exists in the standard ΛCDM model.Although the model analyzed here does not alleviate this fine tuning, it is an attempt to unify the dark sector by describing it through a single scalar field.In addition, this model can be viewed as a first step in obtaining a unified scenario that could describe a more dynamic dark energy.This possibility has become especially interesting after the recent BAO data released by the Dark Energy Spectroscopic Instrument (DESI) [73,74].Equation (2.5) above is inspired in the Wess-Zumino model [60].The potential is adjusted such that the stable minimum is the zero of the potential.The potential has minima at φ = 0 and φ = 2m 3λ .The general shape of the potential is shown in Figure 1.Following the previous work [48], we will consider the value λ = 10 −2 for the coupling constant.Although this choice can seem fine tuned at first, we will show later that the results are not very sensitive to the value of λ.In addition, values for the coupling term with a magnitude not much bigger or much smaller than 1 can be viewed as natural choices in the context of Quantum Field Theory.In the case we are interested in, with Γ > 0 and energy flowing from the cosmological constant into dark matter, it follows that the value of the dark matter density parameter extrapolated to high redshifts would be higher than the value predicted by ΛCDM, while the expansion rate would be lower than that in ΛCDM.If the energy transfer occurs at a non-negligible rate, this affects cosmological observables, which is why it is important to consider observational constraints in such a scenario.
In the work of [51], a generic model of metastable DE decaying into DM was analyzed and constraints on the value of the decay rate were obtained.Later in [41], a further analysis of this model was performed using more recent data.In the latter, they considered four combinations of datasets in their analysis: (1) Pantheon compilation [52] in combination with BAO data from 6dFGS [53], MGS [54], and BOSS DR12 [55].(2) They add BAO data from eBOSS DR14 [56] to the first dataset.(3) They add high redshift BAO measurement from Lyα [57] to the second dataset.(4) They included the CMB distance prior [58,59] to the full combination of datasets.These four analyses constrained, respectively, the following values: .More details on their analysis and methodology can be found in [41].We can use these constraints in order to estimate an upper limit on the decay rate and to analyze the scenario that results from this limiting value.Although the constraints differ depending on the dataset used, we can consider Γ/H 0 ≤ O(10 −1 ) as a reasonable estimate for the upper limit.The analysis we are going to perform in the following sections will provide us with the mass of the resulting dark matter particle (Section III), and the fraction of energy currently remaining in the false vacuum (Section IV).

III. DARK MATTER FROM FIRST ORDER PHASE TRANSITION
The energy transfer in the MDE model occurs due to the tunneling from the metastable (false) to the stable (true) vacuum of the potential, Equation (2.5).This process can be described by a first-order phase transition according to the semi-classical method developed in [65,66].Using this framework, we are going to compute the mass of the resulting DM particle as a function of Γ.
The decay rate (per unit volume) reads where φ + is the field amplitude at the false vacuum and φ(ρ) is, in analogy with the case of particles, the classical field amplitude in Euclidean space crossing the potential −V (φ) subject to boundary conditions φ initial = φ final = φ + .S E is the Euclidean action and S Λ is the Euclidean action of a particle at the false vacuum.Looking for an analytical solution, we will consider the so-called thin wall approximation [65,66], as conducted in the previous work of [48].
In this approximation, the energy difference between the two minima of the potential, given by the parameter ϵ, is considered to be small, then perturbative results can be obtained in terms of ϵ.
The classical equation of motion in the Euclidean space for the field φ subject to the potential V (φ) is obtained by the minimization of S E , which reads We also suppose that lim where τ = −it.
The solution is Euclidean invariant, which means that φ(⃗ x, τ , the equation of motion now reads which is analogous to the equation of motion of a particle at position φ as a function of time ρ, subject to a potential −V (φ) and a friction term.The decay process occurs by the formation of bubbles of a true vacuum surrounded by the false vacuum outside.The field is at rest both inside and outside, and the friction-like term, 3 ρ ∂ ∂ρ φ, is nonzero only at the bubble wall.Therefore, in the case the wall is thin, we can consider that ρ = R at the wall, being R the bubble radius (see [65,66] for further details).Finally, for a small energy difference ϵ, the quantity R = ρ is large and the friction coefficient 1/ρ → 0, then we can neglect the friction term also at the bubble wall.From these considerations, Equation (3.4) now reads Then, considering the potential in Equation (2.5), the value of the field is each region is which corresponds to the regions inside the bubble (0 < ρ ≪ R), at the thin wall (ρ ≈ R) and outside the bubble (ρ ≫ R).Now, we compute the action by adding the separated contribution of each region, which reads where ∆ represents the width of the wall, and we defined Minimizing the action S with respect to R, one obtains which gives the solution In the decay rate expression, Equation (3.1), the exponential term will dominate.We also know the pre-exponential term has the dimension of energy to the 4th power and that it will weight the overall value; then, we simply estimate it as 1 GeV 4 in order to provide correct units.Therefore, Equation (3.1) can be written as From the above equation, we can see that while λ is raised to the eighth power, the mass is raised to the twelfth power.This shows that the results we are going to obtain in this work are not very sensitive to the choice of λ.
Finally, by substituting the values of ϵ and λ, we obtain for the decay rate (per volume) One can verify that, within the validity of the semiclassical approach we are considering, maintaining the preexponential term in Equation (3.1) would have changed the result at most by a factor of 2. This factor of 2 would be multiplied by a factor of 10 157 in the equation above.Therefore, it is really sufficient for our purposes to consider the pre-exponential term to be of order 1 GeV 4 .
Due to the symmetry of our problem, we can invert Γ/V , take its fourth root, and interpret it as the decay time, t decay .If we equate the decay time to the age of universe, t 0 = 10 17 s, as was conducted in [48], the value m ∼ 10 −13 GeV is obtained.This result implies in a mass of an axion-like DM particle [48].The equation of state of this real scalar field can be associated with a DM particle since, after the phase transition, oscillations about the quadratic minimum of the potential are the important aspect to determine the equation of state [67,68].
We felt motivated to further explore this model in the light of the new cosmological data [41].In order to do so, let us now consider the case when the inverse decay rate is not equated to the age of the Universe, but to an arbitrary decay time.We parameterize this case by introducing the parameter α, assuming arbitrary values smaller than 1.If we consider the quantity (Γ/V ) −1/4 to be related to the age of the universe t 0 , the quantity (α Γ/V ) −1/4 will imply in a decay time of the form t 0 /(α) 1/4 .As a consequence, different choices of the parameter α (different choices of the decay rate/decay time) result in different values for the mass m of the DM particle.
In Figure 2, we illustrate the dependence of the mass of the resulting DM particle on α, with α ranging from 10 −10 to 10 −1 .From this figure, we notice that there is no change in the order of magnitude of the mass for this whole range of values of α.
The behavior shown in Figure 2 can be described by the analytical expression which shows the weak dependence of the mass on α.The order of magnitude of the mass only changes for α ∼ 10 −107598 , which corresponds to a decay time 10 26900 t 0 .Obviously, such a case is not interesting for us.Therefore, for any non-negligible (and observationally allowed) decay rate, the prediction for the DM mass, m ∼ 10 −13 GeV, remains consistent.In addition to that, one can verify that these results also do not change if one considers other similar values for the coupling constant, as λ = 10 −1 for instance.
As is well known, first-order phase transitions occur through the random nucleation of bubbles.In this process, firstly the false vacuum energy is transferred to the kinetic energy of the bubbles wall, which asymptotically expands at the speed of light [66,[75][76][77].However, in some cases, at some point, the walls of these bubbles may percolate, and the energy of the walls is converted into particles (energy density) that eventually thermalize [72,[78][79][80].Let us analyze this process further in the following section.1.650 × 10 -13 1.655 × 10 -13 1.660 × 10 -13 1.665 × 10 -13 1.670 × 10 -13 α m/GeV FIG. 2. Dependence of the mass on the parameter α, with the latter ranging from 10 −10 to 10 −1 .

IV. BUBBLE NUCLEATION
First-order phase transitions, as the one considered here, occur through the nucleation and growth of bubbles of the new phase 6 .This is similar to the process that happens in old inflationary models.In the case of old inflation, in order for the decay process to occur efficiently, it was necessary to have an inflation decay rate higher than a certain minimum value, in order to end inflation successfully.However, such a high decay rate would imply a bubble nucleation process, which leads to a highly inhomogeneous Universe.Given the similarities between the phase transition described here and the process that plagued the old inflationary model, we think it is important to discuss some implications of such a process in the context of the model considered here.
In our low-temperature late-time model for the dark sector, the bubble nucleation occurs through Coleman-Callan tunneling and the nucleation rate is essentially temperature-independent.Therefore, we consider the idealized model for the transition described in [71], in which the Universe is taken to be always at zero temperature, with negligible curvature.We can also approximately consider a de Sitter exponential expansion, a(t) ∝ e H t , in the late time Universe, H being the Hubble parameter assumed to be approximately constant.This is a sufficiently good approximation for our practical purposes.The de Sitter approximation can also be justified by the fact that the recent constraints that we are considering for the decay time implies a decay process that will only be effective in a distant future when the Universe is even more close to a de Sitter expansion.
We consider that bubble nucleation begins at a time t B , and afterward occurs at a constant rate per unit of physical volume Γ/V 0 .The decay rate per unit of coordinate volume is then given by Γ(t)/V = (Γ/V 0 )e 3Ht .We can consider the approximation in which a bubble starts expanding from a very small radius r 0 ≪ H −1 .In order to show that this is the case, let us estimate the initial radius of the bubble from which it starts expanding.Using Equations (3.10) and (3.11), we can estimate this radius to be where in the last equality we used the particle mass determined in Section III, m ≈ 10 −13 GeV, together with the values λ = 10 −2 and ϵ = 10 −47 GeV 4 .In a successful first-order phase transition, most of the bubbles are nucleated and collide in a time interval comparable to or shorter than the Hubble time.This is the so-called fast transition, in which the nucleation rate is comparable to the expansion rate of the Universe [71,72].After bubbles collide, in successful cases, the distribution of energy in the Universe becomes homogeneous as the relevant bubbles are sub-horizon sized when they collide.In the slow transitions, the nucleation rate is much smaller than the expansion rate of the Universe.We will show that in the case of the MDE model considered here, rare and very large bubbles are nucleated in a cosmological time.These bubbles can grow to astrophysical sizes and their dynamics cover a much longer period of time.After such a bubble is nucleated, it soon starts expanding at the speed of light [65,66].For a bubble that emerges at a point (t B , ⃗ x 0 ) with a very small radius, at a later time t its radius will be given by Evaluating the integral above gives us the result r(t, t B ) = H −1 (e −Ht B − e −Ht ).We can see that the bubble radius will asymptotically assume the finite value hence, the volume of the Universe occupied by this bubble will asymptotically be Due to the approximately de Sitter expansion, the bubble will grow in comoving size for only about a H −1 time, and after that, being super-Hubble, it simply stretches with the scale factor as the universe expands.This suggests that independently from the time past after the nucleation, two bubbles that emerge simultaneously at points (t B , ⃗ x ′ 0 ) and (t B , ⃗ x ′′ 0 ) at a proper distance D = a(t B ) × |⃗ x ′ 0 − ⃗ x ′′ 0 | > 2H −1 will never percolate.It is possible to analyze the evolution of the bubble distribution in the Universe by investigating the physical volume remaining in the false vacuum, given by V phys ∝ a 3 (t)p(t), where p(t) is the probability that a given point in space is in the false vacuum at a time t.The quantity p(t) can be written in the following form [71,72,[85][86][87] where I(t) is the expected volume of true-vacuum bubbles per unit volume of space at time t7 .For a phase transition beginning at time t 0 , we can expect this volume to be given, at time t, by the expression where a constant decay rate was considered Γ for the MDE model.Above, r(t, t ′ ) is the coordinate radius at time t of a bubble that was nucleated at a time t ′ , expressed in Equation (4.2).Note that we are integrating in the nucleation times.Above, we considered that the decay rate in the MDE model is a constant.The equation above can be better understood if we view the quantity I(t) as a function of the scale factor rather than t.In this case, we can rewrite the equation above as Above, we used Equation (4.2) for r(t, t ′ ).The integrals in this expression are dominated by the upper limit of the integration range.We can see that the quantity (Γ/V 0 )/H 4 sets the magnitude of I(t).Up to a numerical factor, this quantity measures the fraction of space occupied by large bubbles, which is given by 1 − p(t).Therefore, provided that this fraction is small, as it will be in all cases of interest, 1 − p(t) ≈ I(t).
Even expanding at the speed of light, a bubble which nucleates at a time t B can only grow to a finite comoving radius given by Equation (4.3) .Then, as shown above, if the separation between two bubbles at time t is greater than 2r A , the bubbles will never meet.Therefore, the bubbles nucleated in a time interval of duration of H −1 can never fill space by themselves, but instead only occupy a fraction of order (Γ/V 0 )/H 4 of the region which remained in the old phase at the time that they were nucleated.Although bubble nucleation continues indefinitely, and despite that the physical volume of the old phase region (proportional to a(t) 3 p(t)) is an increasing function of time, the bubbles produced have smaller and smaller comoving volume and so can fit in the remaining regions of old phase without overlapping.
Another simple way of understanding this is remembering that after a bubble is nucleated and grown to a size of about a Hubble radius, its size simply conformally stretches as the Universe expands.Then, from that point on, the volume fraction of the Universe that it occupies remains constant.The volume fraction occupied by the bubbles nucleated during the time interval ∆t = H −1 is roughly the volume of a bubble when it begins conformally stretching, which is around H −3 , times the number of such bubbles nucleated in this interval per unit volume given by Γ/(V 0 H).Therefore the quantity (Γ/V 0 )/H4 indicates the volume fraction of space occupied by bubbles nucleated over a Hubble time at a given epoch.We will denote by ξ this important quantity, A slow transition is considered to be the one in which the quantity above is much smaller than one, ξ ≪ 1, and the Coleman-Callan tunneling is the only significant mechanism of bubble nucleation.When this is not the case, the transition is considered to be fast.
The decay time in our model is given by the quantity t decay = (Γ/V ) −1/4 .As discussed in Section III, we can write this decay time as t decay ≡ α −1/4 t 0 .Therefore, we can write the quantity ξ as We can consider in the above equation the absolute mean values for Γ/H 0 described in Section II, which were obtained in the analysis of [41] (see Table 4 of this reference).For the values constrained with any of the four datasets used, Equation (4.9) gives us ξ ≪ 1.Therefore, one can safely consider the slow transition regime.
We can hence conclude that such a model will not drive a complete transition to a dark matter-dominated phase.Even in the future, there will be a dominant region of the universe in the metastable phase, where we would still have an approximate de Sitter expansion.However, unlike old inflation, there is no observational restriction imposing that a complete transition from the dark energy state into dark matter must take place.
As shown in this section, the initial radius of the bubble from which it starts expanding is R 0 ≈ 10 −2 cm.Although it expands very fast, the volume fraction of space occupied by bubbles nucleated over a Hubble time at a given epoch is much smaller than one, as shown above.Therefore, there are negligible cosmological consequences from the energy density associated with these bubbles.Since here we are only interested in the cosmological scales, we will not discuss further the issue of domain walls or gravitational waves.

V. CONCLUSION AND PROSPECTS
We considered the model proposed in [48], in which an energy transfer from dark energy into dark matter is described in field theory by a first-order phase transition.We further investigate this model in light of the recent cosmological data.We find that the model is not excluded by the data, although it is not distinguishable from ΛCDM.Since recent data constrain a decay time for metastable dark energy considerably larger than the age of the Universe, there is currently no prospect of observing the outcome of the DE decay process, which we show to be a dark matter particle with a mass corresponding to that of an axion-like particle, m ∼ 10 −13 GeV.We analyzed, for the first time, the process of bubble nucleation in this model, showing that such a model would not drive a complete transition to a dark matter-dominated phase even in a distant future.
In particular, in this work, we provided the following answers to some, until now, open questions: 1. Considering the recent cosmological data, the model proposed in [48] can still be considered, formally, a viable model for describing an unified dark sector.