On the Catalysis of the Electroweak Vacuum Decay by Black Holes at High Temperature

We study the effect of primordial black holes on the classical rate of nucleation of AdS regions within the standard electroweak vacuum at high temperature. We find that the energy barrier for transitions to the new vacuum, which determines the exponential suppression of the nucleation rate, can be reduced significantly, or even eliminated completely, in the black-hole background if the Standard Model Higgs is coupled to gravity through the renormalizable term $\xi {\cal R} h^2$.


Introduction
It is well-known that for the current measured values of the Higgs and top quark masses the Standard Model (SM) effective potential develops an instability. Due to the running of the quartic coupling, the effective potential of the Higgs reaches a maximum and then becomes unbounded from below at values of the Higgs field of about 5 × 10 10 GeV. Our electroweak vacuum has a lifetime which is many orders of magnitude larger than the age of the Universe and is therefore stable against decay through quantum tunnelling.
A pertinent question concerns the fate of the electroweak vacuum during the evolution of the Universe in situations in which gravity plays a pronounced role. The gravitational background can have a significant effect on the rate of vacuum decay, leading to its enhancement or suppression [1]. This issue is crucial for the decay of the electroweak vacuum, because of the extreme sensitivity of the Higgs potential to the Higgs and top masses [2,3,4,5,6,7]. The fact that the gravitational effects are significant during inflation has led to an extensive investigation of the stability of the electroweak vacuum during this era, leading to bounds on the inflationary scale [4,8,9,10,11] in order to avoid a catastrophic singularity characterized as the AdS crunch [1,10,12].
Strong gravitational fields are also sourced by black holes, which raises the question of electroweak vacuum stability in their vicinity. It has been argued that black holes can act as impurities enhancing the quantum decay rate to a level incompatible with the age of the Universe [13]. It is natural to expect a similar enhancement for classical transitions to the true Higgs vacuum, induced by the high-temperature environment of the early Universe [14]. This second mechanism can be explored through a more intuitive approach, with fewer technical uncertainties than the calculation of the quantum tunnelling rate. In any case, the question whether the existence of primordial black holes [15,16] is consistent with the stability of the electroweak vacuum has important implications for the compatibility of the Standard Model of particle physics and the cosmological model [17].
In this paper we explore further the scenario of temperature-induced vacuum decay in the black-hole background [14]. Firstly, we repeat the analysis using the temperature-corrected effective potential of the Higgs field, instead of the zero-temperature one employed in [14]. Secondly, we allow for a nonminimal coupling of the Higgs field to gravity. This additional coupling can have strong effects on the action of the bounce in the case of quantum tunnelling, or the energy of the critical bubble in the case of thermal tunnelling, because it modifies the effective Higgs mass in regions of large curvature. Depending on the sign of the coupling, the vacuum decay rate may be reduced if the mass grows and fluctuations are suppressed, or an instability may appear if the mass is made to vanish. Such behavior is known to occur in the context of inflation [10] and similar features are expected within the strong gravitational field of a black hole.
The paper is organized as follows. In section 2 we derive the equations of motion and the relevant expressions for the energy of the bubble configuration. In section 3 we study the finite-temperature effects on the Higgs potential and determine the size of the barrier to be overcome for vacuum decay in the presence of black holes. In section 4 we discuss the stability of the electroweak vacuum and derive constraints for the allowed range of the nonminimal coupling of the Higgs field to gravity. In the final section we give our conclusions.

Equations of motion
Our starting point is the action of the SM Higgs coupled to gravity, where we also allow the renormalizable term whose strength is parametrized by the coupling ξ The equations of motion are The use of the equations of motion allows to simplify the action further. In particular, taking the trace of Eq. (2) expresses the curvature scalar R in terms the scalar field. The derivative of the potential can be eliminated through Eq. (3). One eventually finds Substitution in Eq. (1) gives The second term in the above expression is reduced to a boundary term upon integration.
We consider now the SM Higgs field in the presence of black holes that can trigger the decay of the electroweak vacuum by providing nucleation sites [13]. The decay takes place through the formation of bubbles of the new vacuum around the black holes. The critical bubble is a static configuration, whose profile can be obtained by solving the equation of motion of the Higgs field with appropriate boundary conditions. The field interpolates between values on either side of the maximum of the potential. The ADM mass for the critical bubble is a measure of the energy barrier that must be overcome for the field to fluctuate beyond the potential maximum. In a thermal environment, the ratio of this mass and the temperature is expected to determine the exponential suppression of the vacuum decay rate.
The appropriate metric for our analysis has the form With this ansatz, one has to solve the Einstein and Higgs-field equations for the static critical bubbles. The presence of the horizon at r = R h requires appropriate boundary conditions at this point. 1 The function N (r) vanishes on the horizon, which gives 2GM (R h ) = R h . The quantity R h /(2G) can be identified with the black-hole mass M bh . In order to obtain a finite ADM mass for the bubble configuration we must avoid a singularity in Eq. (9) at r = R h . This can be achieved by choosing h (R h ) appropriately. We also require that h(r) → 0 for r → ∞. We are interested in the quantity δM (r) = M (r) − R h /(2G). The asymptotic value δM tot ≡ δM (∞) gives an estimate of the barrier associated with the production of the bubble around a central black hole with horizon radius R h . The characteristic scale of the solutions is set by the largest of h max , the value of the SM Higgs at the maximum of its potential 2 and the temperature T , which are expected to be much smaller than the Planck scale M Pl . This leads to a significant simplification of the setup, resulting from the fact that δM (r) R h /(2G). Since the black-hole mass is the leading factor determining the gravitational background, we write M (r) = R h /(2G) + δM (r) and keep the leading contributions in G. We 1 When the spacetime has a boundary, the consistency of the variational principle requires the presence of a boundary term [18] with (g ind ) αβ the induced metric on the boundary surface, = ±1, depending on whether the surface is timelike or spacelike, respectively, and K the trace of its extrinsic curvature. We solve the equations of motion, starting slightly outside the horizon, assuming the presence of such a term with = −1.
2 While this definition is not gauge-invariant, one could define h max as V 1/4 max which, thanks to the Nielsen identity, is gauge-invaviant [19].
find that the equations of motion are reduced to with the prime denoting a derivative with respect to r. The equation for δ is reduced to δ = 0, and we set δ = 1. Avoiding the singularity at r = R h in Eq. (8) requires that we impose This boundary condition correctly reproduces the standard condition h (0) = 0 in the absence of a black hole. Our definition of δM (r) imposes δM (R h ) = 0. The value of h(R h ) must be tuned so that the condition h(r) → 0 for r → ∞ is satisfied. The form of the above equations demonstrates that the Higgs configuration is not modified by the nonminimal coupling to gravity, and, therefore, is independent of ξ. The bubble has a negligible effect on the gravitational background, which is determined by the black-hole mass to a very good approximation. Allowing for a nonminimal Higgs coupling to gravity through a term ∼ h 2 R in the action does not modify this conclusion, because the curvature R vanishes in a black-hole background. As a result, all modifications to the background due to the Higgs are suppressed by M Pl .
However, the mass of the bubble configuration receives a correction proportional to ξ. It can be expressed as with It is possible to put the above expressions in an alternative form. The equation of motion (8) can be written as Substitution into the expression (12) gives Figure 1: The quantities F 1 and F 2 , that determine the bubble mass δM tot , as a function of the black-hole radius. All quantities in units of the value of the Higgs field h max at the maximum of the potential.
Also, the use of Eq. (14) in (13) gives In Ref. [14] the zero-temperature potential was used in order to compute the mass of the bubble configuration as a function of the black-hole mass for ξ = 0. The bubble mass was then used in order to estimate the vacuum decay rate at non-zero temperature. In the following section we shall improve on this estimate by considering the temperaturecorrected Higgs potential, which corresponds to the free-energy density of the system. Before doing so, we present here the generalization of the results of Ref. [14] for nonzero ξ. For the Higgs potential we use the approximate form V ∼ λ(h) h 4 /4, which is a good approximation for values of the Higgs field around its maximum at h max ∼ 5 × 10 10 GeV. The quartic coupling λ varies from 0.02 to −0.02 for Higgs values between 10 6 GeV and 10 20 GeV, respectively. Near the maximum the potential can be approximated as [10] In Fig. 1 we plot the functions F 1 and F 2 that determine the bubble mass δM tot , as a function of the black-hole radius. All quantities are measured in units of h max . We have normalized F 1 (R h ) with respect to its value F 1 (0) in the absence of a black hole and we depict it in the left panel of Fig. 1. It is apparent that for ξ = 0 the ADM mass of the bubble configuration is reduced by an approximate factor of 2 within a certain range of black-hole masses around R h = 10 h max . In the right panel of Fig. 1 we depict the ratio It vanishes for R h → 0, but quickly grows and remains almost constant for R h > ∼ 20 h max , with a value 6.6.
One may wonder if the problem for ξ = 0 can be analyzed in the Einstein frame as well. There are subtleties in the transition to this frame, related to the geometry that we consider. In our analysis we treat the horizon as a boundary, assuming the presence of a Gibbons-Hawking term there in order to make the variational problem well-defined. The reformulation of the problem in the Einstein frame has to account for both the form of the background and the boundary term. We find that our approach is the most transparent.

Finite-temperature effects
The next step is to consider the immersion of the black hole-Higgs system in the thermal environment of the early Universe. This framework raises certain conceptual issues that need to be clarified. In a thermal environment, the quantity relevant for transitions between different states is the free energy, which accounts for the effect of entropy. In this respect, it is natural to employ the temperature-dependent effective potential, which can be identified with the free-energy density. On the other hand, the gravitational field is sourced by the energy-momentum tensor, which includes the energy density and is identified with the zero-temperature potential. The formal derivation of the appropriate expressions on a strong gravitational background is difficult. For example, the effects of temperature are usually taken into account by compactifying the Euclidean time direction. On a black-hole background, the compactification scale is set by the Hawking temperature, which does not coincide necessarily with the ambient temperature.
As we saw in the previous section, the backreaction of the Higgs field is negligible, so that the background is described by the Schwarzschild metric to a very good approximation. As a result, the contributions of the thermal environment to the Einstein equations need not be considered. The equations relevant for our problem, Eqs. (8), (9), determine the shape and energy of the Higgs configuration that the fluctuating system has to go through for the transition to occur. At nonzero temperature, energy is expected to be replaced by free energy. It is then justifiable intuitively to replace the zero-temperature potential with the high-temperature one. The resulting equations have the correct limits for vanishing temperature or black-hole mass. We emphasize, however, that a formal derivation is lacking.
The characteristic scale of the solutions is set by the largest of h max and the temperature T . For T h max , the form of the zero-temperature potential near its maximum at h max is not relevant for our calculation. The temperature effects shift the maximum of the potential to a value proportional to T . The only dimensionful scale in the problem is the temperature, which determines the scale at which the running couplings must be evaluated. The temperature corrections to the Higgs potential are summarized in Ref. Figure 2: The rescaled Higgs field outside a black hole with rescaled radiusR h = 0.01, 0.5, 2 (lines from left to right). [6]. The Higgs field develops a thermal mass m 2 T = κ 2 T 2 , with κ 2 = 1 12 The running gauge, τ -Yukawa and quartic running couplings must be evaluated at a scale set by the temperature.
A semi-analytical treatment of the effects on thermal tunnelling arising from the blackhole background and the nonminimal Higgs-gravity coupling is possible if we approximate the full potential as The quartic λ is taken to be constant with a value corresponding to the running coupling at a scale set by the temperature. In the temperature range between 10 14 and 10 18 GeV, we have κ 0.3, with a 10% decrease for increasing T . Also, λ −0.015 with a 15% decrease for increasing T [6]. Neglecting the logarithmic running of the quartic coupling is a reasonable approximation. Our discussion can be repeated easily for lower reheating temperatures, but a fully numerical analysis is necessary if quantitative precision is required.
Through the rescalings the total free energy of the bubble configuration can be expressed as κT /|λ| times a numerical factor [20]. Implementing the above, along with results in Eq. (8) being replaced by a similar equation for the rescaled variables and V (h) =h 2 /2−h 4 /4. Its solution is presented in fig. 2 for three different black-hole masses. It is apparent that the characteristic scale of the Higgs configuration is comparable to the horizon radius forR h = O(1). For such black holes, the Hawking temperature satisfies . There is no thermal equilibrium between the black hole and the environment, while the effect of the Hawking radiation on the ambient thermal bath can be neglected. We shall see below that black holes withR h 0.5 have a strong effect on the vacuum decay rate [14]. This effect is purely gravitational and is localized in the vicinity of the black hole. As noted in Ref. [17], the analysis of the quantum decay rate [13] results in an expression that indicates that the compactification of the Euclidean time direction at a radius ∼ 1/T H , linked to the Hawking temperature, plays a role at arbitrarily large distances from the black hole. The effect we are discussing is not of similar nature, and has a more robust intuitive interpretation. The free energy of the bubble configuration can be put in the form with κ/|λ| ∼ 20 and We find F 1 (0) 18.9, in agreement with Ref. [20]. Also F 2 (0) = 0, as can be checked through a partial integration, use of the equation of motion, and remembering that h (0) = 0 for R h = 0. Our approximations result in a value for the ratio δM tot /T = 18.9κ/|λ| in the absence of a black hole that is independent of the temperature: δM tot /T 380. The complete analysis of the renormalized action of the SM Higgs field is consistent with this value, but also indicates a 20% decrease of the action for temperatures between 10 14 and 10 18 GeV [6,7]. At lower temperatures, the logarithmic running of the quartic term in the potential becomes important and numerical determination of the bubble profile is needed for every value of T . However, we expect that the effect of the black hole on the ratio δM tot /T can be described again by a multiplicative factor (F 1 (R h ) + ξ F 2 (R h ))/F 1 (0), as in eq. (22). This expectation is supported by the analysis of Ref. [14], in which the zero-temperature potential was used. Figure 3: The dimensionless quantities F 1 and F 2 , that determine the ratio δM tot /T , as a function of the rescaled black-hole radius.
In Fig. 3 we depict the functions F 1 (R h ) and F 2 (R h ). We have normalized F 1 (R h ) with respect to its value F 1 (0) in the absence of a black hole. The form of F 1 (R h ) is similar to that observed in Ref. [14]: the free energy of the bubble configuration for ξ = 0 is reduced by an approximate factor of 2 within a certain range of black-hole masses around R h = 0.5/(κT ) 1.7/T . As a result, the nucleation rate for ξ = 0 can be reduced significantly in the black-hole background. In the right panel of Fig. 3 we depict the ratio F 2 (R h )/F 1 (R h ). It vanishes forR h → 0, but quickly grows and remains almost constant forR h > ∼ 0.5, with a value 5.3.
A very interesting feature is the effect of the nonminimal coupling on the bubble free energy. The dependence on ξ is linear, as can be seen from Eq. (22). Large positive values of ξ result in the growth of δM tot /T and the suppression of the nucleation rate. As we mentioned in the introduction, the nonminimal coupling can have a strong effect because it modifies the effective Higgs mass in regions of large curvature. Depending on the sign of the coupling, the rate may be reduced if the mass grows and fluctuations are suppressed, or an instability may appear if the mass is made to vanish. Such behavior has been discussed in the context of inflation [10].
From Eq. (22) it is apparent that, for every black-hole mass at a given temperature, there is a critical negative value ξ cr = −F 1 (R h )/F 2 (R h ), for which the exponential suppression of the nucleation rate is eliminated. Clearly, the saddle-point approximation breaks down before this point is reached. On the other hand, the probability of vacuum decay in the early Universe depends also on the total number of primordial black holes that can be generated. In the following section we discuss this point in detail.

Vacuum (in)stability in the presence of black holes
The nucleation probability per unit time in the background of a black hole during the radiation-dominated era of the early Universe is dP/dt ∼ T exp(−δM tot /T ). There is no volume factor because of the absence of translational invariance. The characteristic time interval that can be associated with the temperature scale T is the Hubble time 1/H ∼ M Pl /T 2 . The nucleation can be efficient over longer times, but we are interested in a lower bound for the probability. Neglecting the evaporation of the black hole, we have a bubble-nucleation probability P ∼ M Pl /T exp(−δM tot /T ). The most uncertain factor in our discussion is the number of primordial black holes within our past light cone. The number N of causally independent regions at a time during the radiation-dominated era with temperature T , which are currently within our horizon, is N ∼ 10 34 (T /GeV) 3 [14]. One needs to estimate the probability p for a black hole to exist within one of these regions. This is a strongly model-dependent quantity and little progress can be achieved without assuming a specific mechanism for black-hole creation. In any case, putting everything together, we find that the logarithm of the probability of electroweak vacuum decay in the presence of black holes of typical radius R h can be written as where M Pl = (8πG) −1/2 2.43 × 10 18 , (δM tot /T ) 0 is computed in the absence of black holes, and Fig. 3.
In Ref. [14] it was pointed out that A(R h T ) can lead to the reduction of the effect of (δM tot /T ) 0 by a factor of roughly 2. Our present analysis demonstrates that the nonminimal coupling ξ to gravity can have a more dramatic effect: For ξ = −1/B(R h T ), the exponential suppression of the nucleation rate can be eliminated completely. Clearly, the saddle-point expansion around the bubble configuration breaks down before this point is reached. On the other hand, the presence of a large prefactor, coming from the huge number of causally independent regions in the early Universe, indicates that the electroweak vacuum is likely to become totally unstable.
As we have mentioned above, the largest uncertainty in the calculation is connected with the probability p to find a black hole within each causally independent region at the time when the ambient temperature is T . A primordial black hole can form when the density fluctuations are sufficiently large for an overdense region to collapse [15]. It is usually assumed that its typical mass is of the order of its maximal possible mass. The latter is given by the total mass within the particle horizon ∼ M 2 Pl /H, while the maximal radius is ∼ 1/H. It is not clear, however, if these assumptions are consistent with the typical size of density fluctuations in the early Universe, or the constraints imposed by the observations of the microwave background, for example.
As a concrete application, we consider the scenario of Ref. [17], which is consistent with the current bounds on the size of primordial density fluctuations. It is assumed that the inflationary era is followed by a period during which the inflaton oscillates and decays into particles. The effective equation of state is similar to that of a matter-dominated Universe. It is possible, therefore, for perturbations that were generated during inflation to reenter the horizon, grow gravitationally and collapse into black holes. When the thermalization of the decay products takes place, the equation of state changes and the growth of perturbations is suppressed. In such a scenario, the black holes most relevant for our discussion are those produced just before thermalization, because they are the most massive ones. This must be contrasted with the scenario of the quantum decay of the electroweak vacuum in which the most relevant black holes are the light ones, because of their large Hawking temperature [13,17].
The modes that result in black-hole creation start from a primordial value δρ/ρ = δ i ∼ 10 −4 at horizon crossing and grow until δρ/ρ ∼ 1, at which point they decouple from the Hubble flow, turn around and collapse. The probability for them to form black holes has been estimated as P BH ∼ 2 × 10 −2 (R h /R) 13/2 [21], with R h the Schwarzschild and R the turnaround radius. The black-hole mass is equal to that contained within the turnaround radius: M = 4πρR 3 /3. The Friedmann equation for a matter-dominated era then gives: i /H. As one expects ∼ (HR) −3 regions of size R within a volume ∼ 1/H 3 , the probability to find a black hole within the horizon is [17] p ∼ 2 × 10 −22 δ i 10 −4

5
. (25) The bubble-nucleation rate is enhanced if, at the onset of the radiation-dominated era, there are black holes for whichR h = κR h T > ∼ 0.5. At that time ρ = (g * π 2 /30)T 4 , where we take g * = 106.75, assuming only the particle content of the Standard Model. We obtainR If the reheating temperature is larger than ∼ 10 −7 M Pl ,R h is too small at the onset of the radiation-dominated era for the black holes to have a siginificant effect on the rate. At later times, black-hole creation is suppressed, the temperature drops, while the black-hole mass is reduced through evaporation. As a resultR h becomes even smaller. It seems, therefore, that a high reheating temperature eliminates the danger posed by the presence of black holes. On the other hand, for T > ∼ 10 −7 M Pl , the typical mass of a black hole in this scenario is M < ∼ 10 27 GeV 10 3 gr. Such small black holes evaporate very quickly and are of little phenomenological interest. In this respect, the possibility of a reheating temperature smaller than 10 −7 M Pl seems more exciting. As an example, let us consider the possibility of a reheating temperature T 5 × 10 11 GeV, for which Eq. (26) with δ i 10 −4 givesR h 0.5, so that A(R h T ) 0.5 and B(R h T ) 5.3 in Eq. (24). Using Eq. (25), we obtain ln(N pP ) 124 − 0.5(1 + 5.3 ξ)(δM tot /T ) 0 , with (δM tot /T ) 0 computed in the absence of black holes. The probability becomes of order one for ξ < ∼ −0.19 + 47/(δM tot /T ) 0 . The stability of the electroweak vacuum in the presence of primordial black holes imposes a strong constraint on the nonminimal coupling of the Higgs field to gravity, by forbidding this range.

Conclusions
Given the absence of new physics in the LHC data, the stability of the SM electroweak vacuum has become a more pressing issue. Special conditions in the early Universe may not be left out when studying the rate of vacuum decay. As a particular example, black holes may trigger the decay, as first suggested in Ref. [13]. In this paper we have analyzed how the presence of primordial black holes may induce the SM electroweak vacuum decay in the early Universe at finite temperature. In particular, our findings indicate that a nonminimal, but renormalizable coupling between the SM Higgs field and gravity may alter considerably the decay rate.
Ultimately, the final fate of the electroweak vacuum is a model-dependent issue, which suffers from our ignorance of the precise early Universe dynamics and the exact black hole mass function as a function of time. Nevertheless, our results indicate that even moderate values of the coupling ξ between the Higgs and gravity can render more stable or unstable the electroweak vacuum, depending on the sign of the coupling. While this calls for a refinement of our analysis, e.g. by determining more precisely the probability p for a black hole to exist in one of the many causally independent regions which are currently within our visible Universe, it demonstrates once more the importance of gravity for the issue of vacuum decay.