Impact of New Physics on the EW vacuum stability in a curved spacetime background

It has been recently shown that, contrary to an intuitive decoupling argument, the presence of new physics at very large energy scales (say around the Planck scale) can have a strong impact on the electroweak vacuum lifetime. In particular, the vacuum could be totally destabilized. This study was performed in a flat spacetime background, and it is important to extend the analysis to curved spacetime since these are Planckian-physics effects. It is generally expected that under these extreme conditions gravity should totally quench the formation of true vacuum bubbles, thus washing out the destabilizing effect of new physics. In this work we extend the analysis to curved spacetime and show that, although gravity pushes toward stabilization, the destabilizing effect of new physics is still (by far) the dominating one. In order to get model independent results, high energy new physics is parametrized in two different independent ways: as higher order operators in the Higgs field, or introducing new particles with very large masses. The destabilizing effect is observed in both cases, hinting at a general mechanism that does not depend on the parametrization details for new physics, thus maintaining the results obtained from the analysis performed in flat spacetime.


I. INTRODUCTION
One of the most important goals of present theoretical and experimental particle physics is the search for New Physics (NP) beyond the Standard Model (BSM), even though direct experimental searches up to now have not revealed any sign of it. When looking for patterns towards BSM theories, the stability analysis [1][2][3][4][5][6][7][8][9][10][11] of the electroweak (EW) vacuum plays a crucial role. Earlier studies were mainly focused on establishing bounds for the Higgs boson mass, based either on the requirement that the Higgs effective potential V (φ) could not take values lower than the EW minimum v (so that the latter is the stable vacuum of the theory) or on the possibility that our Universe sits on a metastable (false) vacuum state (i.e. V (v) is not the absolute minimum of V (φ)), with a lifetime larger than its own age [12], [13], [2], [10].
On the theoretical side, the stability analysis has its roots in a pioneering work of Bender and collaborators [40], where the tunneling for a quantum mechanical system with several degrees of freedom was studied with the help of saddle point techniques. This was later extended to quantum field theory by Coleman and Callan, who studied the decay of the false vacuum in a flat spacetime background [41], and then by Coleman and De Luccia [42], who included gravity in their analysis.
Physically the false vacuum decay is triggered by quantum fluctuations, that induce a finite probability for a bubble of true vacuum to materialize in a false vacuum sea. Both in flat and curved spacetime backgrounds, Coleman and collaborators considered a scalar theory where the potential V (φ) has a relative and an absolute minimum, at φ false and φ true respectively, such that the energy density difference V (φ false ) − V (φ true ) is much smaller than the height of the "potential barrier" V (φ top ) − V (φ false ), where V (φ top ) is the maximum of the potential between the two minima. Under these conditions the true vacuum bubble is separated from the false vacuum sea by a "thin wall", and this allows to treat the problem analytically, within the so called "thin wall" approximation.
Going back to the SM, it is known that due to the top loop corrections the Higgs potential V (φ) bends down for values of φ > v, where v ∼ 246 GeV is the location of the EW minimum, and for the present experimental values of M H and M t , namely M H ∼ 125.09 GeV and M t ∼ 173.34 GeV [43,44], it develops a second minimum, much deeper than the EW one and at a much larger value of the field, φ true v = φ false . The instability scale of the Higgs potential is then identified as the value φ inst of the field such that V (φ inst ) = V (v) and V (φ) < V (v) for φ > φ inst . For the values of the Higgs and top masses reported above, it turns out that φ inst ∼ 10 11 GeV. Clearly the conditions under which the thin wall approximation can be applied are not fulfilled in the SM case, so the results of [41] and [42] cannot be directly applied.
The EW vacuum stability condition was first studied in a flat spacetime background, and the interesting possibility that the SM is valid all the way up to the Planck scale M P , or more generally to some very large scale M large , meaning that NP shows up only at this scale, was investigated. In such a scenario, naturally prompted by the lack of direct observation of hints of new physics, the analysis was performed under the assumption that the presence of NP at M large could be neglected for the computation of the tunneling time τ from the false to the true vacuum of the SM, so that τ was calculated by considering SM interactions only [10,11,14,[30][31][32]. In fact it was argued that the relevant scale for tunneling is the instability scale φ inst ∼ 10 11 GeV, and that the contribution to the tunneling rate coming from NP that lives at the scale M large should be suppressed. In other words, as φ inst ∼ 10 11 GeV M large , a decoupling effect was expected [11]. It was later realized that the assumption that NP lives at M large ( φ inst ) does not imply that it cannot affect the stability condition of the EW vacuum. On the contrary, the latter turns out to be very sensitive to unknown NP even if it lives at scales far away from φ inst , and the expected decoupling phenomenon does not take place [33][34][35]38].
The reason why the decoupling theorem does not hold in this case is that tunneling is a non-perturbative phenomenon [38], while the former applies when calculating scattering amplitudes in perturbation theory at energies E much lower than M large . In this case the contributions to scattering amplitudes from physics that lives at M large is suppressed by factors of E/M large to the appropriate power, and this is how physics at the scale M large is decoupled from physics at the scale E.
For our tunneling phenomenon however, the bulk of the contribution to τ comes from the exponential of the (Euclidean) action calculated at the saddle point of the path integral for the tunneling rate, the so called bounce solution to the (Euclidean) Euler-Lagrange equation [41], and for this tree level contribution no suppression factors of the kind (E/M large ) n can ever appear. If the Higgs potential is modified by the presence of NP at M large , the new bounce is certainly different from the one obtained when these terms are neglected. The action calculated for this new bounce solution is also modified and (once exponentiated) it can give rise to a value of τ enormously different from the one obtained when the NP terms are neglected.
The inclusion of gravity in the vacuum stability analysis was pioneered in [42], where the case of the thin wall regime was studied. For the transition from a false Minkowski vacuum to a true Anti-de Sitter (AdS) vacuum, it was shown that, when the size of the Schwarzschild radius of the true vacuum bubble is much smaller than its size, i.e. when gravitational effects are weak, the probability of materialization of such a bubble is close to the flat spacetime result, while when the Schwarzschild radius becomes comparable to the bubble size, i.e. in a strong gravitational regime, the presence of gravity stabilizes the false vacuum, preventing the materialization of a true vacuum bubble. In other words, gravity tends to stabilize the false vacuum, and in a strong gravity regime the materialization of bubbles of true vacuum is quenched.
The above results are obtained in the thin wall regime, but they are commonly considered of more general validity. In particular it was recently claimed [45][46][47] that in the presence of Planckian physics (strong gravity regime) the new bounce solutions [33][34][35]38] to the (Euclidean) equation of motion modified by the presence of new physics at M large should disappear, and that as a consequence these NP terms could not induce any modification to the tunneling rate.
As noted above however, the SM case is very far from the thin wall regime analyzed in [42], and before jumping to any conclusion, and also prior to studies that consider the inclusion of new physics, the stability analysis for the SM in the presence of gravity has to be performed. An early attempt to study the impact of gravity on the EW vacuum decay rate was done in [48], where a perturbative expansion of the bounce around the flat spacetime solution was considered. As later noted in [49] and then in [50] however, the boundary conditions for the bounce solution, that are essential in the calculation of the decay rate, are not respected already at the first order of the expansion, and the results of this work are at least questionable.
In order to get close to the SM case, but still keeping a simple model as in [42], a scalar theory with a potential whose parameters can be adjusted to explore cases far from the thin wall regime was considered in [49], and a numerical analysis of the false vacuum stability condition was performed. The main result is that for the potential that well approximates the SM case, the stabilizing effect of gravity is hardly seen even in very strong gravity regimes. As suggested in [51][52][53], the total quenching of the vacuum decay rate can eventually be reached at some very high scale. As shown in [49] however, for the SM case such an effect takes place in a far transplanckian regime where the theory has already lost its validity (similarly to what happens for the Landau pole in QED, where the latter occurs at such a high energy scale that the theory has lost its significance several orders of magnitudes below that scale). The results obtained with the simple model considered in [49] were later confirmed in [50], where a bona fide SM Higgs effective potential was used.
The issue raised in [45][46][47] however, namely the possibility that strong gravity effects should make the presence of Planckian NP harmless in the calculation of the tunneling time, is a crucial open question. In order to complete the stability analysis of the EW vacuum it is of the greatest importance to understand if gravity really cancels the destabilizing effect induced by NP at high energies, as claimed in [45][46][47].
In the present work we address these issues, that are very important for current studies and for model building of BSM physics, where we are often confronted with the possibility of considering NP at Planckian and/or trans-Planckian scales. Anticipating on the results of the following sections, we will see that the tunneling time from the false to the true vacuum is still strongly dependent on NP even if it lives at very high ( φ inst ) scales, thus confirming the results of the analysis performed in the flat spacetime background.
The rest of the paper is organized as follows. In Section II the general theoretical set-up for our work, mainly consisting of the equations that will be used in the subsequent numerical analysis, is presented. Moreover, in order to keep the present paper as self-contained as possible, and also to check our tools against known results, the EW vacuum stability analysis in the absence of new physics in both flat and curved spacetime backgrounds is briefly sketched, and the known results are recovered. In section III we study the impact of NP at the Planck scale on the stability condition of the EW vacuum when the presence of gravity is taken into account, parametrizing NP in terms of higher order operators. In Section IV a different parametrization for NP at high energy scales is used, namely we introduce a new boson and a new fermion (with very large masses) coupled to the Higgs boson. Section V is for our conclusions.

II. THEORETICAL BACKGROUND
In the present section we briefly review the theoretical background for the computation of the tunneling decay rate from a Minkowski false vacuum (minimum of the potential with vanishing energy density) to an Anti-De Sitter (AdS) true vacuum (minimum of the potential with negative energy density), considering both the flat and the curved spacetime background cases.
Flat spacetime. Let us begin by considering the flat spacetime Euclidean action for a single component real scalar field φ: where V (φ) is a potential with a local minimum (false vacuum) at φ = φ fv , and an absolute minimum (true vacuum) at φ = φ tv . In order to calculate the false vacuum lifetime we have to look for the so called bounce solution to the Euclidean Euler-Lagrange equation that have O(4) symmetry and satisfy certain boundary conditions [41]. If r is the radial coordinate, the equation takes the form: where the dot indicates derivative with respect to r, and the boundary condition are: Denoting with φ b (r) the bounce solution, the action at φ b is: and the decay rate Γ of the false vacuum is given by: where is the so called tunneling exponent and the exponential of −B gives the "tree-level" contribution to the decay rate, while D is the quantum fluctuation determinant. If V (φ fv ) = 0, the action S[φ fv ] vanishes, and the tunneling exponent is simply In order to determine the false vacuum decay rate in the flat spacetime case, in the following we integrate numerically Eq. (2) with boundary conditions (3), and use (4) to get the tunneling exponent B of (5). Curved spacetime. The next step is to study the impact of gravity on the vacuum decay rate, and to this end we consider the previous theory in a curved spacetime background. Including the Einstein-Hilbert term, the Euclidean action becomes: where R is the Ricci scalar and G is the Newton constant. Requiring again O(4) symmetry, the (Euclidean) metric takes the form: where dΩ 2 3 is the unit 3-sphere line element and ρ(r) is the volume radius of the 3-sphere at fixed r coordinate [42]. The bounce configuration needed to calculate the false vacuum transition rate is now given by the field and the metric solution, φ b (r) and ρ b (r) respectively, of the coupled equations (κ ≡ 8πG): where the first equation replaces (2), while the second is the only Einstein equation left by the symmetry. For the decay of a Minkowski false vacuum to a true AdS vacuum, the case of interest to us, the boundary conditions are: Asymptotically (r → ∞) the bounce φ b (r) approaches the constant false vacuum solution φ fv , where V (φ fv ) = 0 (Minkowski vacuum). In the same limit, from (8) 2 we see that the bounce solution metric ρ b (r) approaches the flat spacetime metric: In the thin wall regime the constant c is obtained analytically [42], while in general it is determined from the numerical integration of (8). Differentiating (8) 2 with respect to r we get: that is a useful equation that we will use in our analysis as it is more robust than (8) 2 for numerical integration [54]. Finally the Ricci scalar R in terms of ρ is given by: Inserting the above results in (6), the action for the bounce (φ b , ρ b ) becomes: As the false vacuum action S[φ fv , ρ fv ] contains the same boundary term (the last one) of (13), in the tunneling exponent B this term does not appear. Neglecting it, and using (8) in (13) we finally get: From (14), the false vacuum action S fv ≡ S[φ fv , ρ fv ] can be easily derived. As already said, we are interested in the decay of a Minkowski false vacuum (the state that corresponds to the minimum of the potential with vanishing energy density, V (φ fv ) = 0) to an AdS true vacuum. For the false vacuum case, φ(r) = φ fv , andρ 2 = 1, i.e. ρ fv (r) = r + const. This latter constant is fixed and is c of Eq. (10), as the bounce metric asymptotically tends to the flat spacetime metric of the Minkowski false vacuum. In calculating the action S fv we must then integrate over the Minkowski spacetime from a real radius equal to zero up to infinity. This corresponds to the integration range [−c, ∞[ in the r coordinate, so that for the tunneling exponent B = S b − S fv we have: In the following, the decay rate of the Minkowski false vacuum into the AdS vacuum will be calculated by integrating numerically Eqs. (8) 1 and (11) with boundary conditions (9). Then with the help of (16) the tunneling exponent will be obtained.
As we are interested in the stability analysis of the EW vacuum, in our case the scalar field φ is the Higgs field, and the potential V (φ) is the Higgs effective potential. More specifically it is the renormalization group improved potential, that can be written as: where λ SM (φ) is the quartic running coupling λ SM (µ) (µ is the running scale) with µ = φ [2,55].
In order to get λ SM (µ), the system of RG equation of the SM couplings has to be run. In [50], where a stability analysis of the SM in the presence of gravity was performed, an improved Higgs potential obtained with the help of three-loops beta functions was used, while in previous analyses two-loop beta functions were considered [14]. In this respect, we note that the consistent counting of loops in the beta functions of different SM particles poses a problem in itself, as it was found that taking the same order for all the different components of the SM (Yukawa, gauge and quartic couplings) actually leads to inconsistencies [56].
The purpose of the present work however is to study the impact that NP at high energies can have on the stability condition of the EW vacuum when the SM coupling to gravity is taken into account. We are then not interested in precision measurements and/or refinements of previous analyses. For our illustrative scopes, the differences between these cases (the two loops counting, the three loops counting, and the "consistent" counting proposed in [56]) are minimal, and they have no impact on the results and conclusions of our analysis.
We can then leave aside these questions and work in a simplified yet very robust framework, by using a good approximation of the SM effective potential that was obtained in [57] by fitting the two-loops improved Higgs potential with the three parameter function [57]: where M P = 1/ √ G is the Planck mass. The fit gives: In the following we work with the Higgs potential (17) with λ SM (φ) given by (18) and (19). Both in the flat and curved spacetime cases, an important parameter is the size R of the bounce, defined as the value of r such that Going back to (5) for the vacuum decay rate, we note that a good approximation to the prefactor for the case that we are considering is given in terms of the bounce size R and of T U , the age of the Universe, and the EW vacuum tunneling time τ = Γ −1 turns out to be [35]: In the following we use (21) to calculate the false vacuum lifetime. Before ending this section and moving to the study of the impact of NP on the EW vacuum stability, we would like to test our tools starting with the known cases of the flat and curved spacetime backgrounds in the absence of NP (i.e. considering the SM alone), and briefly sketch the analysis for these cases.
Flat spacetime. In order to proceed with the numerical solution of the bounce equation (2), we begin by scaling the dimensionful field φ and the radial coordinate r to dimensionless quantities, x and ϕ(x) respectively, by using Planck units: Eq. (2), the boundary conditions (3) and the potential (17) then become: where the prime indicates the derivative respect to x. After the rescaling (22), the tunneling exponent (4) becomes: Solving numerically the bounce equation (23), with the Higgs potential given by (25) and (19), and inserting the result for ϕ b (x) in (26), after using the values found for B and R, namely B = 2025.27 and R = 10.7597, we finally get for the lifetime τ of the EW vacuum: in very good agreement with the results known in the literature. This is the first test of our numerical method, and also shows that we are considering a good approximation for the Higgs potential.
Curved spacetime. As in the case of flat spacetime, we move to dimensionless quantities. Defining the dimensionless curvature a(x) = M P ρ(r), Eqs. (8) 1 and (11) (that will be used in the following numerical integration) become: where the potential U (ϕ) is the same as in (25). The corresponding boundary conditions are: As we have already said, ρ(r) ∼ r for r → ∞, and the asymptotic (x → ∞) behavior of the bounce solution in the presence of gravity is the same as in the flat spacetime case.  (25) with the parameters λ * , α, β given in (19). The center of the bounce is at ϕ(0) = 0.0712, its size is R = 350.2996 and the tunneling exponent is B = 2062.5836. Right panel: Difference between the curvature radius and its asymptotic value, a(x) − x, for the same parameters as in the left panel.
In terms of dimensionless quantities, from (16) we find for tunneling exponent: where c (as already said above) is the constant that determines the asymptotic behavior of the metric a(x) for x → ∞, while (ϕ b , a b ) is the bounce solution to the system (28).
In the left panel of Fig. 1 the bounce profile ϕ b (x) is plotted. The right panel shows the difference a b (x) − x: we clearly see how asymptotically a b (x) reaches the Minkowskian behavior a(x) ∼ x + c, and we can read the value of the constant c. Finally, with the help of (21), we obtain the tunneling time in the presence of gravity: Once again we observe that the above result is in good agreement with known results [50]. Moreover, comparing (31) with the corresponding flat spacetime tunneling time (27), we see that gravity (as expected) tends to stabilize the EW vacuum.

III. NEW PHYSICS: HIGHER ORDER OPERATORS
The results briefly presented in Section II are known and concern the stability analysis under the assumption that, in the event that unknown NP lives at a scale M large much greater than the instability scale φ inst , it has no impact on the stability condition of the EW vacuum. In other words it is assumed that new physics at M large is decoupled from the physics that triggers the EW vacuum decay, and that it should be possible to calculate the tunneling rate ignoring these terms.
The analysis of the previous section is essential to set the proper framework where the effects of the presence of NP at M large can be properly investigated. To be specific, we consider the case where NP lives at the Planck scale M P (i.e. we choose M large = M P ) and parametrize it as in [33][34][35] with the help of higher powers of φ added to the Higgs potential: It was shown in [33][34][35] for the flat spacetime case that when λ 6 < 0 and λ 8 > 0 the potential (32) destabilizes the EW vacuum. In other words, these NP terms favor the nucleation of true vacuum bubbles and, depending on the specific values of λ 6 and λ 8 , this destabilization effect could dramatically reduce the EW vacuum lifetime τ in (27) and make it even shorter than the age of the Universe T U . We now consider the same kind of analysis in the presence of gravity.
Adding the NP terms (32) to the SM Higgs potential (17), and moving again to dimensionless quantities, the new dimensionless potential U (ϕ) becomes: We are now ready to study the impact of high energy NP on the EW vacuum stability condition in the presence of gravity.
A first important result of our analysis is that for each value of the couple (λ 6 , λ 8 ) there is a different bounce solution to Eqs. (28), all of them being different from the solution obtained for the SM alone, i.e. the case λ 6 = 0, λ 8 = 0.
Contrary to the expectations of [45][46][47] then, gravity does not induce the disappearance of the "new" bounce solutions related to the presence of new physics, here parametrized in   With the help of (21) we now calculate the EW vacuum lifetime for different values of the NP couplings λ 6 and λ 8 . The fourth column of Table I contains different values of the tunneling time obtained for different couples (λ 6 , λ 8 ). For comparison, the third column contains the corresponding values of τ for the flat spacetime analysis. First of all we note that the effect already seen in the previous section (also reported in the first line of the table, the case λ 6 = 0, λ 8 = 0), namely that the presence of gravity tends to stabilize the EW vacuum, is maintained even in the presence of new physics.
However, a simple inspection of this table shows that even though the presence of gravity tends to stabilize the EW vacuum as compared to the corresponding flat spacetime case, still for O(1) values of the new physics couplings λ 6 and λ 8 the tunneling time can be made smaller than the age of the Universe T U . Let us consider just a couple of examples. For λ 6 = −0.3 and λ 8 = 0.3 for instance, the EW vacuum in the flat spacetime background is unstable, being τ ∼ 10 −52 T U , but for the corresponding case with gravity included we observe a stabilization of the EW vacuum: τ ∼ 10 287 T U . There is a competition between the destabilizing effect of NP and the stabilizing effect of gravity. In this example, gravity takes over new physics and as a result the EW vacuum turns out to be stable. However for larger (absolute) values of the NP couplings, the destabilizing effect of NP takes over the stabilizing effect of gravity. For instance, for λ 6 = −1.2 and λ 8 = 1.0, despite the stabilizing effect of gravity (τ grav τ flat ), the EW vacuum turns out to be unstable: τ grav ∼ 10 −58 T U . The results discussed above with the help of Table I are better summarized in Fig. 3 where stability diagrams in the (λ 6 , λ 8 ) plane are presented for the range of values −1.5 < λ 6 < 0.4 Let us summarize the results of the present section. On the one hand we have shown that the inclusion of gravity in the stability analysis of the EW vacuum does not cause the disappearance of the "new" bounce solutions found in the flat spacetime case [33][34][35] that in turn cause the destabilization of the EW vacuum to the point that its lifetime can become smaller than the age of the Universe. On the other hand we have seen that the stability condition of the EW vacuum is the result of a competition between the destabilizing effect of NP and the stabilizing tendency of gravity. In particular we have also seen that outside a certain range of values of the NP couplings λ 6 and λ 8 , the destabilizing effect of Planckian NP always takes over, thus making the EW vacuum unstable.

IV. NEW PHYSICS: FERMIONS AND BOSONS WITH LARGE MASSES
In the present section the stability analysis of the EW vacuum will be performed by considering a different parametrization for NP at high energy scales. Actually in [33][34][35] and in the previous section the analysis was performed by parametrizing NP at the Planck scale in terms of few higher order (non-renormalizable) operators. This is just a convenient and efficient way of mimicking the presence of new physics, clearly not an (illegitimate) truncation of the UV completion of the SM. Some authors however expressed a certain skepticism on these results, suggesting that this effect should disappear when the infinite tower higher dimensional operators of the renormalizable UV completion of the SM is taken into account, so that the expected decoupling of very high energy physics from the mechanism that triggers the decay of the false vacuum should be recovered. It was actually suspected that this effect takes place above the physical cutoff, where the control of the theory is lost [58].
Although it is understandable that the parametrization of NP in terms of higher order operators could be the source of a certain confusion, the destabilizing effect has nothing to do with this parametrization. For the case of a flat spacetime background in [38] the stability analysis was performed by parametrizing NP in terms of renormalizable additional terms, with a fermion and a boson with very high masses that interact with the Higgs field, and it was shown that the destabilizing effect found in [33][34][35] is still present.
In this section we present the same kind of analysis of [38] taking into account the presence of gravity (i.e. considering the case of a curved spacetime background), and show that as for the case of the parametrization used in the previous section, gravity does not produce any washing out of the destabilizing effect of new physics, although it slightly mitigates it.
In order to illustrate the destabilization phenomenon we consider as in [38] a renormalizable model that is not a realistic high energy UV-completion of the SM but is very appropriate to the purposes of the present work. NP lives at very high energy scales and is parametrized by adding to the SM a scalar field S and a fermion field ψ that interact in a simple way with the Higgs field φ, with masses M S and M f of the scalar and fermion respectively well above the instability scale: M S , M f φ inst . Apart from the kinetic terms, the additional terms in the Lagrangian are: To understand how a NP Lagrangian of this kind can arise in a physical setup, we note that the large mass term M f can be thought as a sort of heavy right handed "neutrino" in the framework of a see-saw mechanism. While the corresponding light "neutrino" is totally harmless for the stability of the EW vacuum, the heavy "neutrino" can play an important role in destabilizing the vacuum. The scalar field S counterbalances the destabilizing effect of ψ. Note that models with new scalar fields coupled to the Higgs (although admittedly unrealistic) have already been used to provide a stabilization mechanism for the Higgs effective potential [59,60].
Before proceeding with our work, it is worth to mention a problem concerning the physical observables involved in the stability analysis, namely the gauge dependence of the effective potential away from the extrema [61][62][63]. In particular, absolute stability bounds on the Higgs mass (formally gauge independent) turn out to be gauge dependent at any order of perturbation theory. Only when a consistent resummation is considered the result can be made gauge independent, providing a slight improvement in mass bounds [63]. The gauge dependence of the instability scale has also been investigated and the range of uncertainty identified [62], but it is not known if it is possible to calculate this quantity in a gauge independent manner. Moreover, for the main quantity of interest to us, namely the false vacuum decay rate, again we know that it is a formally gauge-invariant quantity. In a truncated perturbative expansion, however, order-by-order gauge independence can possibly be achieved only after resumming the appropriate terms, as it was done for the energies at the minima of the effective potential [61]. Having these warnings in mind, and waiting for improvements on these gauge dependence issues, we proceed now with the analysis of our model following the usual pattern.
For the purposes of the present work, it is sufficient to consider the impact of these additional terms on the Higgs effective potential V (φ) at the one-loop level only. In fact we do not need a better level of precision as we are not interested in extracting numbers but we only want to illustrate the destabilization effect that arises from very high energy physics (see also the considerations developed below Eq. (17)). The one-loop contribution to V (φ) from these terms is: where the renormalization scale µ is taken as µ = M S . In this respect we note that at very high values of the running scale the SM quartic coupling reaches a plateau: λ SM (µ) has practically the same value in the whole range [M f , M S ], and this is why we can use λ SM (M S ) as the threshold value for the coupling (even though strictly speaking we should use λ SM (M f ) (see below)), and can choose µ = M S as the renormalization (threshold) scale. The presence of the high energy NP of (34) is then taken into account by adding to the SM potential V SM (φ) in (17) and (18) the contribution coming from V 1 (φ). To this end we have to implement the matching conditions described below. First of all we expand the potential V 1 (φ) in powers of φ and isolate the constant, the φ 2 and the φ 4 terms. Then at the threshold scale M f we require that: (i) the renormalized cosmological constant Λ, given by the sum of all the constant terms (those coming from the SM potential and those coming from V 1 (φ)) vanishes, Λ(µ = M f ) ∼ 0; (ii) the renormalized mass term, given by the sum of all the coefficients of φ 2 , and identified with the SM mass parameter m 2 SM (µ = M f ) at the scale M f , vanishes: m 2 SM (µ = M f ) ∼ 0 (more precisely we neglect this term to a very high degree of accuracy for the large values of φ considered); (iii) the renormalized quartic coupling, given by the sum of all the coefficients of φ 4 , is identified with the SM quartic coupling at the scale M f , λ SM (µ = M f ). In other words, at the scale M f this coefficient is matched with the value of the quartic coupling obtained by considering the running of the renormalization group equations for the SM couplings alone.
The above requirements for the renormalized cosmological constant and mass are well known features. For the renormalized Λ (apart from the fine tuning problem) we can practically consider that Λ(µ = 0) ∼ Λ(µ = M f ) ∼ 0. The same is true for the renormalized mass, for which we take m 2 (µ = 0) ∼ m 2 (µ = M f ) ∼ 0, meaning that we neglect the φ 2 term as compared to the φ 4 and other terms for these large values of φ, and that the running of the renormalized mass is totally harmless in this respect. For the quartic coupling we have a true matching condition. In fact we require that at the threshold scale µ = M f the quartic coupling coincides with λ SM (µ = M f ), that is obtained by running the renormalization group equations for the SM couplings only. Practically starting from the scale M f , the potential is given by the SM contribution V SM (φ) plus the contribution of V 1 (φ) subtracted of its constant, quadratic and quartic powers of φ, that we call V 1 (φ) from now on: We  (28) for ϕ(x) and a(x). In Fig. 5 the profiles of the bounce solutions ϕ b (x) for the four different cases (i), (ii), (iii) and (iv) and the corresponding plots of a b (x) − x are presented, with colors yellow, blue, green and red respectively.
These are the first relevant results of the present section. They show that, contrary to the claims in [45][46][47], the presence of gravity does not cause the disappearance of the new bounce solutions due to the presence of high energy NP, that were already found in the case of the analysis carried in the flat spacetime background [38]. Therefore they confirm the results of the previous section, where high energy NP was parametrized in terms of higher order operators, and once again show that the appearance of new bounce solutions is not an artifact of the specific parametrization used in Section III.
Using (21) to calculate the vacuum lifetime, for the examples considered above we find in units of T U (going from (i) to (iv)): Eqs. (37) and (38) together with Fig. 5 contain the main lesson of the present section. They definitely show that, even when gravity is included in the analysis, the presence of NP at high energy scales (even though much higher than the instability scale) can have an enormous impact on the vacuum lifetime. It is worth to remind here that when the calculation is performed in the curved spacetime background and the presence of high energy new physics is not considered, the tunneling time is given by (31) (τ ∼ 10 661 T U ), while from (37) we see that τ strongly depends on the parameters of new physics, and can turn out to be even shorter than the age of the Universe.
Moreover, by comparing (37) and (38) we see that gravity slightly pushes toward stabilization (that is also what we observe for the case when new physics is not taken into account by comparing (27) with (31)), but this is a "tiny" effect, that (as we have just seen) does not generate the claimed [45][46][47] disappearance of the new bounce solutions.
Despite the fact that in our toy model NP lives at very high energy scales, much higher than the instability scale φ inst , the expectation that the tunneling time should be insensitive to it, in other words that the result shown in (31) should not be modified by the presence of NP at high energies, is not fulfilled. These results confirm the analysis of the previous section. Here, with the help of a fully renormalizable toy UV completion of the SM, we have shown that the EW vacuum lifetime strongly depends on NP even if the latter lives at very high energy scales.
These findings are at odds with a widely diffused expectation, based on a naive application of the decoupling argument, and show that the fact that the vacuum stability condition depends on physics that lives at very high energy scales is not due to an illegitimate extrapolation of the theory beyond its validity, as it was previously thought [64]. On the contrary, it is an illegitimate application of the decoupling argument to a phenomenon to which it cannot be applied, namely the (non-perturbative) tunneling phenomenon, that leads to the expectation that physics at scales much higher than the instability scale φ inst should have no impact on the stability condition.
Before ending this section, we would like to stress once again that with respect to the previous section, where NP interactions were parametrized with the help of higher order non-renormalizable operators, here NP is given in terms of a fully renormalizable theory, thus showing that the effect that we present is a genuine physical effect and has nothing to do with the specific parametrization of NP.

V. SUMMARY, CONCLUSIONS, AND OUTLOOK
We studied the impact of very high energy NP (around the Planck scale M P ) on the stability condition of the EW vacuum by carrying the analysis in a curved spacetime background, i.e. by taking into account the presence of gravity. Despite the expectation [45][46][47] that strong gravity effects should act against the formation of new true vacuum bubbles, thus invalidating the results of previous analyses [33][34][35]38] carried in a flat spacetime background, we found that these new solutions to the bounce equations persist even in the presence of gravity, and as in the case of the flat spacetime background analysis, they can have an enormous impact on the EW vacuum lifetime.
As in [33][34][35] we first performed the analysis by adding to the SM potential higher powers of the the Higgs field, more precisely terms as φ 6 /M 2 P and φ 8 /M 4 P (Section III) that are certainly generated in a quantum gravity context [65]. Following [38] we then parametrized high energy new physics in a different manner, namely by adding to the SM potential a boson S and a fermion ψ, with very large masses M S and M f , coupled to the Higgs boson. As for the analysis carried in flat spacetime, in both cases we find that the presence of new physics can have an enormous impact on the EW vacuum stability condition.
These results definitely show that, irrespectively of the parametrization used to describe high energy new physics, it is not possible to ignore its presence when the stability analysis is performed. They are of the greatest importance for current studies and for model building of Beyond Standard Model physics, where we often have to take into account new physics at very high (Planckian and/or trans-Planckian) scales.
A question that is left open by the present analysis is the role that could be played by a non-minimal coupling of gravity to the Higgs boson. Work in this direction is in progress [66]. 2. Ifφ(x max ) −φ(x ref ) < − , the scalar field is overshooting the top of the hill. The initial guess for ϕ(x min ) was therefore too high, and B is correspondingly decreased by a quantity δ.
Furthermore, the algorithm stores a state consisting of both B and its value for the two preceding iterations. It is therefore possible to detect oscillations in B and bisect (or otherwise decrease) δ. We found that decreasing δ to δ/10 each time B changes trend leads to a particularly efficient search, that converges exponentially to the solution, as will be illustrated below.
The inclusion of gravity does not modify our algorithm, as we found that, for x → ∞, the areal radius goes as a ∼ x + c (see text), i.e. the curvature tends asymptotically to zero (the value of the constant c is given by a(x) − x when the plateau is reached).
For this reason, the criteria for tuning the initial value of ϕ, introduced in the flat case, can also be used on a curved spacetime. Again, we implemented the boundary condition in x min : ϕ(x min ) = B ϕ (x min ) = 0 a(x min ) = a (x min ) = 1 , with 1 (clearly we cannot use a(x min ) = 0 due to the factor a −1 in (28) 1 ). After solving the equations, the size R of the bounce solution is obtained and we can compute the integral in Eq. (30).
We go now back to the flat spacetime case and include the NP terms of (32), thus obtaining the (dimensionless) potential (33). The bonus for our analysis is that this potential does not modify the asymptotic behavior of the bounce solution ϕ b (x) for x → ∞, as in this limit we still have U (ϕ(x)) → 0. Therefore the inclusion of NP does not lead to any substantial change in our numerical method. The only modification with respect to the flat spacetime case concerns the expansion of the bounce solution around the origin x = 0. In the flat case, by considering the integration range [x min , x max ], we found that the initial values for solving Eq. (2) with the shooting method are obtained once we take the expansion (A2) of ϕ(x) and its first derivative at x min . Repeating the same analysis for the potential (33), we find that NP simply leads to additional terms in these expressions. Thus, the new expansion for ϕ(x) is given by (again up to O(x 2 )): ϕ(x) = B + B 3 8 λ * + λ 6 B 2 + λ 8 B 4 + α 2 ln B + α ln 2 B + β ln 3 B + β ln 4 B x 2 + · · · (A5) which we use to set initial conditions for ϕ(x min ) and ϕ (x min ). An analogous approach is followed when we consider the alternative parametrization of NP given in Eq. (36). Just like in the case without new physics, when we include gravity we use the initial values (A4), and the numerical integration of the equations of motion is reduced to the tuning of the parameter B. Fig. 6 illustrates the exponential convergence of our algorithm in the four cases presented in this appendix.