Next Article in Journal
Quantum Field Theory of Black Hole Perturbations with Backreaction: I General Framework
Previous Article in Journal
Bayesian Knowledge Infusion for Studying Historical Sunspot Numbers
Previous Article in Special Issue
The de Sitter Swampland Conjectures in the Context of Chaplygin-Inspired Inflation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Constraints on Metastable Dark Energy Decaying into Dark Matter

by
Jônathas S. T. de Souza
1,
Gustavo S. Vicente
2 and
Leila L. Graef
1,*
1
Physics Institute, Fluminense Federal University, Niterói 24210-346, RJ, Brazil
2
Faculty of Technology, Rio de Janeiro State University, Resende 27537-000, RJ, Brazil
*
Author to whom correspondence should be addressed.
Universe 2024, 10(9), 371; https://doi.org/10.3390/universe10090371
Submission received: 8 August 2024 / Revised: 11 September 2024 / Accepted: 11 September 2024 / Published: 18 September 2024

Abstract

:
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.

1. 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 today1. 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 GeV 4 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 energy3. 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,3,4,5,6,7,8,9,10,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 behavior4, 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:
  • Is MDE still a viable model considering the latest cosmological constraints?
  • What is the mass of the DM resulting from this process?
  • Would the bubble nucleation process in this model lead to inhomogeneities that could plague the model?
  • 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 2, we present the MDE model proposed in [48] to describe an interacting dark sector model based on field theory. In Section 3, 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 4, we analyze the bubble nucleation process and its subsequent dynamics in order to see if the model would predict observable inhomogeneities. In Section 5, we conclude and discuss some future prospects. In Appendix A, we demonstrate the validity of the approximations considered.

2. 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:
ρ ˙ DE + 3 H ρ DE ( 1 + ω DE ) = Q DE ,
ρ ˙ DM + 3 H ρ DM ( 1 + ω DM ) = Q DM ,
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:
Q DM = Q DE = Γ ρ DE .
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
ω DE ( eff ) = ω DE + Γ 3 H .
From Equation (4), we see that in a DE-DM interaction model, the effective EoS can have a phantom-like behavior ( ω DE ( eff ) < 1 ) when the decay rate is negative (energy flow DM → DE), or a quintessence-like behavior ( ω DE ( eff ) > 1 ) 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
V ( φ ) = 2 m φ 3 λ φ 2 2 + Q ( φ ) ,
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 minimum5. Of course, this is a fine-tuned 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 (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 φ = 2 m 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: (1) Γ / H 0 = 0 . 47 0.5 + 0.5 (mean with 1 σ deviation); (2) Γ / H 0 = 0 . 07 0.3 + 0.4 ; (3) Γ / H 0 = 0 . 51 0.25 + 0.27 ; (4) Γ / H 0 = 0 . 02 0.01 + 0.01 . 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 3), and the fraction of energy currently remaining in the false vacuum (Section 4).

3. 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 (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
Γ V = S E 2 ( φ ˜ ( ρ ) ) ( 2 π ) 2 × d e t ( μ μ + V ( φ ˜ ( ρ ) ) d e t ( μ μ + V ( φ + ) ) 1 2 × e S E S Λ ,
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
δ S E ( φ ( x ) ) δ φ = μ μ φ ( x ) + V ( φ ) = 0 .
We also suppose that
lim τ ± φ ( x , τ ) = φ + ,
where τ = i t .
The solution is Euclidean invariant, which means that φ ( x , τ ) = φ ( ( | x | 2 + τ 2 ) 1 2 ) . Defining ρ ( | x | 2 + τ 2 ) 1 2 , the equation of motion now reads
2 φ ρ 2 + 3 ρ ρ φ V ( φ ) = 0 ,
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 (9) now reads
2 φ ρ 2 = V ( φ ) .
Then, considering the potential in Equation (5), the value of the field is each region is
φ = 0 i f   0 < ρ R , φ ˜ i f   R Δ < ρ < R + Δ , 2 m / 3 λ i f   ρ R ,
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
S S E S Λ 2 π 2 0 R Δ d ρ ρ 3 ( ϵ ) R Δ R + Δ d ρ ρ 3 1 2 d φ ˜ d ρ 2 + U + R + Δ d ρ ρ 3 0 ( GeV 4 ) = 2 π 2 ϵ R 4 4 + 2 π 2 R 3 R Δ R + Δ d ρ 1 2 d φ ˜ d ρ 2 + U = π 2 R 4 2 ϵ + 2 π 2 R 3 S 1 ,
where Δ represents the width of the wall, and we defined
S 1 R Δ R + Δ d ρ 1 2 d φ ˜ d ρ 2 + U .
Minimizing the action S with respect to R, one obtains
d S d R = 2 π 2 R 3 ϵ + 6 π 2 R 2 S 1 = 0 ,
which gives the solution
R = 3 S 1 ϵ .
For small ϵ , integrating Equation (10), we obtain ρ φ = 2 U . Using this expression and considering the potential from Equation (5), one obtains from Equation (13)
S 1 = 2 4 m 3 27 λ 2 .
Finally, using Equation (16) in Equation (12), the action will result in
S m 12 λ 8 ϵ 3 .
In the decay rate expression, Equation (6), the exponential term will dominate. We also know the pre-exponential term has the dimension of energy to the 4 t h 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 (6) can be written as
Γ V = exp m 12 λ 8 ϵ 3 GeV 4 .
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)
Γ V = exp 10 157 m GeV 12 GeV 4 .
One can verify that, within the validity of the semiclassical approach we are considering, maintaining the pre-exponential term in Equation (6) 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 GeV4.
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
m ( α ) = 1.12246 × 10 13 94.8245 ln α 12 ,
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.

4. Bubble Nucleation

First-order phase transitions, as the one considered here, occur through the nucleation and growth of bubbles of the new phase6. 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 3 H t . 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 (15) and (16), we can estimate this radius to be
R 0 = 2 12 m 3 27 λ 2 ϵ 10 2 cm ,
where in the last equality we used the particle mass determined in Section 3, 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
r ( t , t B ) = x x 0 = t B t d t a ( t ) 1 .
Evaluating the integral above gives us the result r ( t , t B ) = H 1 ( e H t B e H t ) . We can see that the bubble radius will asymptotically assume the finite value
r A lim t r ( t , t B ) = H 1 e H t B ,
hence, the volume of the Universe occupied by this bubble will asymptotically be
V asym = 4 π 3 H 3 e 3 H t B .
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 > 2 H 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 p h y s 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]
p ( t ) = e I ( t ) ,
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
I ( t ) = 4 π 3 t 0 t d t ( Γ / V 0 ) a 3 ( t ) r 3 ( t , t ) .
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 (22). 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
I ( t ) = 4 π 3 R 0 R ( t ) d a ( Γ / V 0 ) H a 2 a R ( t ) d a 1 a 2 H 3 .
Above, we used Equation (22) 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 (23). Then, as shown above, if the separation between two bubbles at time t is greater than 2 r 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 ) / H 4 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,
ξ Γ / V 0 H 4 .
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 3, we can write this decay time as t d e c a y α 1 / 4 t 0 . Therefore, we can write the quantity ξ as
ξ = t decay 4 t 0 4 = α .
We can consider in the above equation the absolute mean values for Γ / H 0 described in Section 2, 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 (29) 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.

5. Conclusions 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:
  • Considering the recent cosmological data, the model proposed in [48] can still be considered, formally, a viable model for describing an unified dark sector.
  • The recent constraints in the decay time of the metastable dark energy imply in a resulting DM with a mass of an axion-like particle, although this resulting DM would only appear in the far future.
  • We do not expect this model to lead to observational imprints that could be searched for in future experiments, unless extra couplings are added to the Lagrangian of the model.
  • The bubble nucleation process was analyzed and we showed that the model considered, apart from not leading to current observable inhomogeneities, would not drive a complete transition to a dark matter-dominated phase, even in the far future.
We can conclude that in order for this model to successfully describe a transition to the true vacuum, the required value of the DE decay rate would need to be bigger than the range allowed by the current observations. Although this model does not currently inherit the characteristics of a typical interacting dark sector model due to the large decay time, it still presents some qualitative advantages, as it is able to describe a dark sector in a unified manner through a single scalar field. In addition, the fact that around the true vacuum, this field behaves as a DM with a mass consistent with the one of an axion-like particle, gives us an indication of the potential of the model. For these reasons, we believe that possible extensions of this model deserve further investigation, as they could lead to potentially observable signatures in case additional couplings are included in the model Lagrangian (in this context see for example the works of [44,45]). The model analyzed here can be viewed as a first step in obtaining a unified scenario that could describe a more dynamic dark energy. Furthermore, recently the new BAO data released by the Dark Energy Spectroscopic Instrument (DESI) [73,74] have encouraged the analysis of more dynamical DE models. In addition to that, there is evidence from different theoretical contexts that exact de Sitter solutions with a positive cosmological constant may not be suitable to describe the late-time universe. DE models based on scalar fields evolving in time are more promising in this regard, although in the context of the Swampland conjecture, for example, they still have to satisfy certain criteria [14,15,16]. Verifying whether extensions of the model here could satisfy these and other theoretical conjectures, and leave observable traces, is an issue left for a future work.

Author Contributions

Conceptualization, J.S.T.d.S., G.S.V. and L.L.G.; Methodology, J.S.T.d.S., G.S.V. and L.L.G.; Formal analysis, J.S.T.d.S., G.S.V. and L.L.G.; Writing—original draft, J.S.T.d.S., G.S.V. and L.L.G. All authors have read and agreed to the published version of the manuscript.

Funding

J.S.T.S. is supported by the Fundação Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES). L.L.G. is supported by research grants from Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Grant No. 307636/2023-2 and from the Fundacao Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro (FAPERJ), Grant No. E-26/201.297/2021. L.L.G. also thank Elisa Ferreira for the important discussions.

Data Availability Statement

No data was generated.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A. Gravitational Effects

In this appendix, we will discuss whether it is necessary to include corrections in the action due to gravitational effects. In order to investigate these effects, let us work with the following action:
S ¯ = d 4 x g 1 2 g μ ν μ φ ν φ V ( φ ) R 16 π G ,
where R is the curvature scalar. Using the thin wall approximation, which is still valid for the action (A1) in the cases of interest, it is straightforward to show the following relation between the action S ¯ and the action S 0 without gravitational effects [88]
S ¯ = S 0 1 + R 0 2 Δ 2 2 .
Above R 0 is the initial radius of the bubble calculated in Equation (21) and Δ is the Schwarzschild radius associated with the bubble of the new vacuum, which is given by Δ = 2 G ϵ ( 4 π R 0 3 / 3 ) . The expression for Δ can be understood from the fact that in the decay process from the metastable vacuum with energy ϵ to the stable vacuum with zero energy, there exists liberation of energy proportional to the energy density of the metastable vacuum and the volume of the nucleated new vacuum bubble. From Equation (A2), we can see that if R 0 / Δ 1 , we need to consider gravitational effects in our calculations. The gravitational effects will be important when the radius of the new vacuum bubble is of the order of the Schwarzschild radius. We can obtain the radius of a nucleated bubble that would be equal to its Schwarzschild radius by equating R 0 = Δ , which gives us
R 0 = 3 8 π G ϵ .
Using Equations (15) and (16), we can see that such an initial radius would correspond to the decay of a particle with a mass of order m 10 2 , 8 GeV , which is much heavier than the particle we obtained. Therefore, the addition of gravitational effects would not significantly alter the results we obtained.

Notes

1
For constraints on DM mass in related contexts see for instance [3,4,5,6,7,8,9,10,11].
2
There are examples in which this configuration can naturally appear as, for example, in the case of the Wess–Zumino model [60], which has a double degenerate bosonic vacuum due to super-symmetry, presumably broken only non perturbatively.
3
Other interesting models which consider a unified dark sector through a single field can be found for instance in [61,62,63,64] and references therein.
4
We can think of any field theory in the context of cosmology as an effective field theory valid until some energy scale. The same is true for our model. We can think of our model as an effective model at low energies. Since this field in the metastable vacuum accounts for today’s dark energy, having an energy density of the order of 10 47 GeV4, it has a negligible cosmological contribution at earlier times. Any quantum correction at early times would be associated with a field with totally negligible contribution to the total energy, having no impact on cosmology, which is the reason why we do not explore these issues in the present work.
5
Another interesting possibility is to consider the scalar field with the even self-interactions up to sixth order, as analyzed in the work of [44].
6
Evolution of bubbles of new vacuum in de Sitter backgrounds has been studied in different context (see, for instance, [75,76,77,80,81,82]). In addition, first-order phase transitions have recently been considered as one of the possible explanations for the positive evidence of a low-frequency stochastic gravitational-wave background found in PTA experiments, see for example [83]. Another recent interesting application of first-order transitions is the New Early dark Energy models, see for instance [84].
7
As discussed in [72], the exponentiation of I ( t ) corrects for some effects, like the fact that when calculating I ( t ) , regions in which bubbles overlap are counted twice. Furthermore, the virtual bubbles which would have nucleated and had their point of nucleation not already be in a true-vacuum region are also included.

References

  1. Sahni, V. Dark matter and dark energy. Lect. Notes Phys. 2004, 653, 141–180. [Google Scholar] [CrossRef]
  2. Ferreira, E.G.M. Ultra-light dark matter. Astron. Astrophys. Rev. 2021, 29, 7. [Google Scholar] [CrossRef]
  3. Amin, M.A.; Mirbabayi, M. A lower bound on dark matter mass. arXiv 2022, arXiv:2211.09775. [Google Scholar] [CrossRef] [PubMed]
  4. Nadler, E.O.; Birrer, S.; Gilman, D.; Wechsler, R.H.; Du, X.; Benson, A.; Nierenberg, A.M.; Treu, T. Dark Matter Constraints from a Unified Analysis of Strong Gravitational Lenses and Milky Way Satellite Galaxies. Astrophys. J. 2021, 917, 7. [Google Scholar] [CrossRef]
  5. Iršič, V.; Viel, M.; Haehnelt, M.G.; Bolton, J.S.; Becker, G.D. First constraints on fuzzy dark matter from Lyman-α forest data and hydrodynamical simulations. Phys. Rev. Lett. 2017, 119, 031302. [Google Scholar] [CrossRef] [PubMed]
  6. Dalal, N.; Kravtsov, A. Excluding fuzzy dark matter with sizes and stellar kinematics of ultrafaint dwarf galaxies. Phys. Rev. D 2022, 106, 063517. [Google Scholar] [CrossRef]
  7. Powell, D.M.; Vegetti, S.; McKean, J.P.; White, S.D.M.; Ferreira, E.G.M.; May, S.; Spingola, C. A lensed radio jet at milli-arcsecond resolution II: Constraints on fuzzy dark matter from an extended gravitational arc. arXiv 2023, arXiv:2302.10941. [Google Scholar] [CrossRef]
  8. Semertzidis, Y.K.; Youn, S. Axion dark matter: How to see it? Sci. Adv. 2022, 8, abm9928. [Google Scholar] [CrossRef]
  9. Nakai, Y.; Namba, R.; Obata, I. Peaky Production of Light Dark Photon Dark Matter. arXiv 2022, arXiv:2212.11516. [Google Scholar] [CrossRef]
  10. Nakatsuka, H.; Morisaki, S.; Fujita, T.; Kume, J.; Michimura, Y.; Nagano, K.; Obata, I. Stochastic effects on observation of ultralight bosonic dark matter. arXiv 2022, arXiv:2205.02960. [Google Scholar] [CrossRef]
  11. Alesini, D.; Braggio, C.; Carugno, G.; Crescini, N.; D’Agostino, D.; Di Gioacchino, D.; Di Vora, R.; Falferi, P.; Gambardella, U.; Gatti, C.; et al. Realization of a high quality factor resonator with hollow dielectric cylinders for axion searches. Nucl. Instrum. Meth. A 2021, 985, 164641. [Google Scholar] [CrossRef]
  12. Carroll, S.; Press, W.; Turner, E. The cosmological constant. Annu. Rev. Astron. Astrophys. 1992, 30, 499–542. [Google Scholar] [CrossRef]
  13. Martin, J. Everything You Always Wanted To Know About The Cosmological Constant Problem (However, Were Afraid To Ask). Comptes Rendus Phys. 2012, 13, 566–665. [Google Scholar] [CrossRef]
  14. Palti, E. The Swampland: Introduction and Review. Fortsch. Phys. 2019, 67, 1900037. [Google Scholar] [CrossRef]
  15. Heisenberg, L.; Bartelmann, M.; Brandenberger, R.; Refregier, A. Dark Energy in the Swampland II. Sci. China Phys. Mech. Astron. 2019, 62, 990421. [Google Scholar] [CrossRef]
  16. Heisenberg, L.; Bartelmann, M.; Brandenberger, R.; Refregier, A. Dark Energy in the Swampland. Phys. Rev. D 2018, 98, 123502. [Google Scholar] [CrossRef]
  17. Polyakov, A.M. De Sitter space and eternity. Nucl. Phys. B 2008, 797, 199–217. [Google Scholar] [CrossRef]
  18. Polyakov, A.M. Infrared instability of the de Sitter space. arXiv 2012, arXiv:1209.4135. [Google Scholar]
  19. Valiviita, J.; Majerotto, E.; Maartens, R. Instability in interacting dark energy and dark matter fluids. JCAP 2008, 07, 020. [Google Scholar] [CrossRef]
  20. Mazur, P.; Mottola, E. Spontaneous Breaking of De Sitter Symmetry by Radiative Effects. Nucl. Phys. B 1986, 278, 694–720. [Google Scholar] [CrossRef]
  21. Mottola, E. A Quantum Fluctuation Dissipation Theorem for General Relativity. Phys. Rev. D 1986, 33, 2136. [Google Scholar] [CrossRef] [PubMed]
  22. Brandenberger, R.; Graef, L.L.; Marozzi, G.; Vacca, G.P. Backreaction of super-Hubble cosmological perturbations beyond perturbation theory. Phys. Rev. D 2018, 98, 103523. [Google Scholar] [CrossRef]
  23. Abramo, L.R.W.; Brandenberger, R.H.; Mukhanov, V.F. The Energy-momentum tensor for cosmological perturbations. Phys. Rev. D 1997, 56, 3248–3257. [Google Scholar] [CrossRef]
  24. Mukhanov, V.F.; Abramo, L.R.W.; Brandenberger, R.H. On the Back reaction problem for gravitational perturbations. Phys. Rev. Lett. 1997, 78, 1624–1627. [Google Scholar] [CrossRef]
  25. Finelli, F.; Marozzi, G.; Vacca, G.P.; Venturi, G. Energy momentum tensor of field fluctuations in massive chaotic inflation. Phys. Rev. D 2002, 65, 103521. [Google Scholar] [CrossRef]
  26. Finelli, F.; Marozzi, G.; Vacca, G.P.; Venturi, G. Energy momentum tensor of cosmological fluctuations during inflation. Phys. Rev. D 2004, 69, 123508. [Google Scholar] [CrossRef]
  27. Marozzi, G. Back-reaction of Cosmological Fluctuations during Power-Law Inflation. Phys. Rev. D 2007, 76, 043504. [Google Scholar] [CrossRef]
  28. Brandenberger, R.H. Back reaction of cosmological perturbations. In Proceedings of the 3rd International Conference on Particle Physics and the Early Universe; World Scientific: Singapore, 2000; pp. 198–206. [Google Scholar] [CrossRef]
  29. de Sá, R.; Benetti, M.; Graef, L.L. An empirical investigation into cosmological tensions. Eur. Phys. J. Plus 2022, 137, 1129. [Google Scholar] [CrossRef]
  30. Di Valentino, E.; Ferreira, R.Z.; Visinelli, L.; Danielsson, U. Late time transitions in the quintessence field and the H0 tension. Phys. Dark Univ. 2019, 26, 100385. [Google Scholar] [CrossRef]
  31. Di Valentino, E.; Melchiorri, A.; Mena, O.; Vagnozzi, S. Interacting dark energy in the early 2020s: A promising solution to the H0 and cosmic shear tensions. Phys. Dark Univ. 2020, 30, 100666. [Google Scholar] [CrossRef]
  32. Zhao, G.B.; Raveri, M.; Pogosian, L.; Wang, Y.; Crittenden, R.G.; Handley, W.J.; Percival, W.J.; Beutler, F.; Brinkmann, J.; Chuang, C.H.; et al. Dynamical dark energy in light of the latest observations. Nat. Astron. 2017, 1, 627–632. [Google Scholar] [CrossRef]
  33. Weinberg, S. The Cosmological constant problems. In Proceedings of the 4th International Symposium on Sources and Detection of Dark Matter in the Universe (DM 2000), Marina del Rey, CA, USA, 23–25 February 2000; pp. 18–26. [Google Scholar]
  34. Ferreira, E.G.M.; Quintin, J.; Costa, A.A.; Abdalla, E.; Wang, B. Evidence for interacting dark energy from BOSS. Phys. Rev. D 2017, 95, 043520. [Google Scholar] [CrossRef]
  35. Amendola, L. Coupled quintessence. Phys. Rev. D 2000, 62, 043511. [Google Scholar] [CrossRef]
  36. Abdalla, E.; Abramo, L.R.; de Souza, J.C.C. Signature of the interaction between dark energy and dark matter in observations. Phys. Rev. D 2010, 82, 023508. [Google Scholar] [CrossRef]
  37. Faraoni, V.; Dent, J.B.; Saridakis, E.N. Covariantizing the interaction between dark energy and dark matter. Phys. Rev. D 2014, 90, 063510. [Google Scholar] [CrossRef]
  38. He, J.H.; Wang, B. Effects of the interaction between dark energy and dark matter on cosmological parameters. JCAP 2008, 06, 010. [Google Scholar] [CrossRef]
  39. Costa, A.A.; Xu, X.D.; Wang, B.; Ferreira, E.G.M.; Abdalla, E. Testing the Interaction between Dark Energy and Dark Matter with Planck Data. Phys. Rev. D 2014, 89, 103531. [Google Scholar] [CrossRef]
  40. Benetti, M.; Borges, H.; Pigozzo, C.; Carneiro, S.; Alcaniz, J. Dark sector interactions and the curvature of the universe in light of Planck’s 2018 data. JCAP 2021, 8, 14. [Google Scholar] [CrossRef]
  41. Li, X.; Shafieloo, A.; Sahni, V.; Starobinsky, A.A. Revisiting Metastable Dark Energy and Tensions in the Estimation of Cosmological Parameters. Astrophys. J. 2019, 887, 153. [Google Scholar] [CrossRef]
  42. Urbanowski, K. Cosmological “constant” in a universe born in the metastable false vacuum state. Eur. Phys. J. C 2022, 82, 242. [Google Scholar] [CrossRef]
  43. Urbanowski, K. A universe born in a metastable false vacuum state needs not die. Eur. Phys. J. C 2023, 83, 55. [Google Scholar] [CrossRef]
  44. Landim, R.G.; Abdalla, E. Metastable dark energy. Phys. Lett. B 2017, 764, 271–276. [Google Scholar] [CrossRef]
  45. Landim, R.G.; Marcondes, R.J.F.; Bernardi, F.F.; Abdalla, E. Interacting Dark Energy in the Dark SU(2)R Model. Braz. J. Phys. 2018, 48, 364–369. [Google Scholar] [CrossRef]
  46. Stojkovic, D.; Starkman, G.D.; Matsuo, R. Dark energy, the colored anti-de Sitter vacuum, and LHC phenomenology. Phys. Rev. D 2008, 77, 063006. [Google Scholar] [CrossRef]
  47. Greenwood, E.; Halstead, E.; Poltis, R.; Stojkovic, D. Dark energy, the electroweak vacua and collider phenomenology. Phys. Rev. D 2009, 79, 103003. [Google Scholar] [CrossRef]
  48. Abdalla, E.; Graef, L.L.; Wang, B. A Model for Dark Energy decay. Phys. Lett. B 2013, 726, 786–790. [Google Scholar] [CrossRef]
  49. Casey, R.; Ilie, C. Dark Sector Tunneling Field Potentials for a Dark Big Bang. arXiv 2024, arXiv:2407.05752. [Google Scholar]
  50. Freese, K.; Winkler, M.W. Dark matter and gravitational waves from a dark big bang. Phys. Rev. D 2023, 107, 083522. [Google Scholar] [CrossRef]
  51. Shafieloo, A.; Hazra, D.K.; Sahni, V.; Starobinsky, A.A. Metastable Dark Energy with Radioactive-like Decay. Mon. Not. Roy. Astron. Soc. 2018, 473, 2760–2770. [Google Scholar] [CrossRef]
  52. Scolnic, D.M.; Jones, D.O.; Rest, A.; Pan, Y.C.; Chornock, R.; Foley, R.J.; Huber, M.E.; Kessler, R.; Narayan, G.; Riess, A.G.; et al. The Complete Light-curve Sample of Spectroscopically Confirmed SNe Ia from Pan-STARRS1 and Cosmological Constraints from the Combined Pantheon Sample. Astrophys. J. 2018, 859, 101. [Google Scholar] [CrossRef]
  53. Beutler, F.; Blake, C.; Colless, M.; Jones, D.H.; Staveley-Smith, L.; Campbell, L.; Parker, Q.; Saunders, W.; Watson, F. The 6dF Galaxy Survey: Baryon Acoustic Oscillations and the Local Hubble Constant. Mon. Not. Roy. Astron. Soc. 2011, 416, 3017–3032. [Google Scholar] [CrossRef]
  54. Ross, A.J.; Samushia, L.; Howlett, C.; Percival, W.J.; Burden, A.; Manera, M. The clustering of the SDSS DR7 main Galaxy sample – I. A 4 per cent distance measure at z = 0.15. Mon. Not. R. Astron. Soc. 2015, 449, 835–847. [Google Scholar] [CrossRef]
  55. Alam, S.; Ata, M.; Bailey, S.; Beutler, F.; Bizyaev, D.; Blazek, J.A.; Bolton, A.S.; Brownstein, J.R.; Burden, A.; Chuang, C.-H.; et al. The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: Cosmological analysis of the DR12 galaxy sample. Mon. Not. Roy. Astron. Soc. 2017, 470, 2617–2652. [Google Scholar] [CrossRef]
  56. Zhao, G.B.; Wang, Y.; Saito, S.; Gil-Marín, H.; Percival, W.J.; Wang, D.; Chuang, C.-H.; Ruggeri, R.; Mueller, E.-M.; Zhu, F.; et al. The clustering of the SDSS-IV extended Baryon Oscillation Spectroscopic Survey DR14 quasar sample: A tomographic measurement of cosmic structure growth and expansion rate based on optimal redshift weights. Mon. Not. Roy. Astron. Soc. 2019, 482, 3497–3513. [Google Scholar] [CrossRef]
  57. Des Bourboux, H.D.M.; Le Goff, J.M.; Blomqvist, M.; Guy, J.; Rich, J.; Yèche, C.; Bautista, J.E.; Bertin, É.; Dawson, K.S.; Eisenstein, D.J.; et al. Baryon acoustic oscillations from the complete SDSS-III Lyα-quasar cross-correlation function at z= 2.4. Astron. Astrophys. 2017, 608, A130. [Google Scholar] [CrossRef]
  58. Chen, L.; Huang, Q.G.; Wang, K. Distance Priors from Planck Final Release. JCAP 2019, 2, 028. [Google Scholar] [CrossRef]
  59. Aghanim, N.; Akrami, Y.; Ashdown, M.; Aumont, J.; Baccigalupi, C.; Ballardini, M.; Banday, A.J.; Barreiro, R.B.; Bartolo, N.; Basak, S.; et al. Planck 2018 results. VI. Cosmological parameters. Astron. Astrophys. 2020, 641, A6, Erratum in Astron. Astrophys. 2021, 652, C4. [Google Scholar] [CrossRef]
  60. Wess, J.; Zumino, B. Supergauge transformations in four dimensions. Nucl. Phys. B 1974, 70, 39. [Google Scholar] [CrossRef]
  61. Brandenberger, R.; Fröhlich, J.; Namba, R. Unified Dark Matter, Dark Energy and baryogenesis via a “cosmological wetting transition”. JCAP 2019, 9, 69. [Google Scholar] [CrossRef]
  62. Brandenberger, R.; Cuzinatto, R.R.; Fröhlich, J.; Namba, R. New Scalar Field Quartessence. JCAP 2019, 2, 43. [Google Scholar] [CrossRef]
  63. Bertacca, D.; Bartolo, N.; Matarrese, S. Unified Dark Matter Scalar Field Models. Adv. Astron. 2010, 2010, 904379. [Google Scholar] [CrossRef]
  64. Frion, E.; Camarena, D.; Giani, L.; Miranda, T.; Bertacca, D.; Marra, V.; Piattella, O. Bayesian analysis of Unified Dark Matter models with fast transition: Can they alleviate the H0 tension? arXiv 2023, arXiv:2307.06320. [Google Scholar]
  65. Callan, C.G., Jr.; Coleman, S.R. The Fate of the False Vacuum. 2. First Quantum Corrections. Phys. Rev. D 1977, 16, 1762–1768. [Google Scholar] [CrossRef]
  66. Coleman, S.R. The Fate of the False Vacuum. 1. Semiclassical Theory. Phys. Rev. D 1977, 15, 2929–2936, Erratum in Phys. Rev. D 1977, 16, 1248. [Google Scholar] [CrossRef]
  67. Marsh, D.J.E. Axion Cosmology. Phys. Rept. 2016, 643, 1–79. [Google Scholar] [CrossRef]
  68. Magana, J.; Matos, T. A brief Review of the Scalar Field Dark Matter model. J. Phys. Conf. Ser. 2012, 378, 012012. [Google Scholar] [CrossRef]
  69. Cicoli, M.; Guidetti, V.; Righi, N.; Westphal, A. Fuzzy Dark Matter candidates from string theory. JHEP 2022, 5, 107. [Google Scholar] [CrossRef]
  70. Harigaya, K.; Leedom, J.M. QCD Axion Dark Matter from a Late Time Phase Transition. JHEP 2020, 06, 034. [Google Scholar] [CrossRef]
  71. Guth, A.H.; Weinberg, E.J. Could the universe have recovered from a slow first-order phase transition? Nucl. Phys. B 1983, 212, 321. [Google Scholar] [CrossRef]
  72. Turner, M.S.; Weinberg, E.J.; Widrow, L.M. Bubble nucleation in first-order inflation and other cosmological phase transitions. Phys. Rev. D 1992, 46, 2384–2403. [Google Scholar] [CrossRef]
  73. Adame, A.G.; Aguilar, J.; Ahlen, S.; Alam, S.; Alexander, D.M.; Alvarez, M.; Alves, O.; Anand, A.; Andrade, U.; Armengaud, E.; et al. DESI 2024 III: Baryon Acoustic Oscillations from Galaxies and Quasars. arXiv 2024, arXiv:2404.03000. [Google Scholar]
  74. Lodha, K.; Shafieloo, A.; Calderon, R.; Linder, E.; Sohn, W.; Cervantes-Cota, J.L.; de Mattia, A.; García-Bellido, J.; Ishak, M.; Matthewson, W.; et al. DESI 2024: Constraints on Physics-Focused Aspects of Dark Energy using DESI DR1 BAO Data. arXiv 2024, arXiv:2405.13588. [Google Scholar]
  75. Teppa Pannia, F.A.; Perez Bergliaffa, S.E. Evolution of Vacuum Bubbles Embeded in Inhomogeneous Spacetimes. JCAP 2017, 3, 26. [Google Scholar] [CrossRef]
  76. Simon, D.; Adamek, J.; Rakic, A.; Niemeyer, J.C. Tunneling and propagation of vacuum bubbles on dynamical backgrounds. JCAP 2009, 11, 008. [Google Scholar] [CrossRef]
  77. Fischler, W.; Paban, S.; Zanic, M.; Krishnan, C. Vacuum bubble in an inhomogeneous cosmology: A Toy model. JHEP 2008, 5, 41. [Google Scholar] [CrossRef]
  78. Hawking, S.W.; Moss, I.; Stewart, J. Bubble collisions in the very early universe. Phys. Rev. D 1982, 26, 2681. [Google Scholar] [CrossRef]
  79. Kosowsky, A.; Turner, M.S.; Watkins, R. Gravitational radiation from colliding vacuum bubbles. Phys. Rev. D 1992, 45, 4514. [Google Scholar] [CrossRef]
  80. Casadio, R.; Orlandi, A. Bubble dynamics: (Nucleating) radiation inside dust. Phys. Rev. D 2011, 84, 024006. [Google Scholar] [CrossRef]
  81. Pannia, F.A.T.; Bergliaffa, S.E.P.; Pinto-Neto, N. Particle production in accelerated thin bubbles. JCAP 2022, 4, 15. [Google Scholar] [CrossRef]
  82. Aguirre, A.; Johnson, M.C. A Status report on the observability of cosmic bubble collisions. Rept. Prog. Phys. 2011, 74, 074901. [Google Scholar] [CrossRef]
  83. Afzal, A.; Agazie, G.; Anumarlapudi, A.; Archibald, A.M.; Arzoumanian, Z.; Baker, P.T.; Bécsy, B.; Blanco-Pillado, J.J.; Blecha, L.; Boddy, K.K.; et al. The NANOGrav 15 yr Data Set: Search for Signals from New Physics. Astrophys. J. Lett. 2023, 951, L11. [Google Scholar] [CrossRef]
  84. Niedermann, F.; Sloth, M.S. New early dark energy. Phys. Rev. D 2021, 103, L041303. [Google Scholar] [CrossRef]
  85. Guth, A.H.; Tye, S.H.H. Phase Transitions and Magnetic Monopole Production in the Very Early Universe. Phys. Rev. Lett. 1980, 44, 631, Erratum in Phys. Rev. Lett. 1980, 44, 963. [Google Scholar] [CrossRef]
  86. Guth, A.H.; Weinberg, E.J. Cosmological Consequences of a First Order Phase Transition in the SU(5) Grand Unified Model. Phys. Rev. D 1981, 23, 876. [Google Scholar] [CrossRef]
  87. Guth, A.H. The Inflationary Universe: A Possible Solution to the Horizon and Flatness Problems. Phys. Rev. D 1981, 23, 347–356. [Google Scholar] [CrossRef]
  88. Coleman, S.R.; De Luccia, F. Gravitational Effects on and of Vacuum Decay. Phys. Rev. D 1980, 21, 3305. [Google Scholar] [CrossRef]
Figure 1. Shape of the potential V ( φ ) .
Figure 1. Shape of the potential V ( φ ) .
Universe 10 00371 g001
Figure 2. Dependence of the mass on the parameter α , with the latter ranging from 10 10 to 10 1 .
Figure 2. Dependence of the mass on the parameter α , with the latter ranging from 10 10 to 10 1 .
Universe 10 00371 g002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

de Souza, J.S.T.; Vicente, G.S.; Graef, L.L. Constraints on Metastable Dark Energy Decaying into Dark Matter. Universe 2024, 10, 371. https://doi.org/10.3390/universe10090371

AMA Style

de Souza JST, Vicente GS, Graef LL. Constraints on Metastable Dark Energy Decaying into Dark Matter. Universe. 2024; 10(9):371. https://doi.org/10.3390/universe10090371

Chicago/Turabian Style

de Souza, Jônathas S. T., Gustavo S. Vicente, and Leila L. Graef. 2024. "Constraints on Metastable Dark Energy Decaying into Dark Matter" Universe 10, no. 9: 371. https://doi.org/10.3390/universe10090371

APA Style

de Souza, J. S. T., Vicente, G. S., & Graef, L. L. (2024). Constraints on Metastable Dark Energy Decaying into Dark Matter. Universe, 10(9), 371. https://doi.org/10.3390/universe10090371

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop