Avoiding the Goldstone Boson Catastrophe in general renormalisable field theories at two loops

We show how the infra-red divergences associated to Goldstone bosons in the minimum condition of the two-loop Landau-gauge effective potential can be avoided in general field theories. This extends the resummation formalism recently developed for the Standard Model and the MSSM, and we give compact, infra-red finite expressions in closed form for the tadpole equations. We also show that the results at this loop order are equivalent to (and are most easily obtained by) imposing an “on-shell” condition for the Goldstone bosons. Moreover, we extend the approach to show how the infra-red divergences in the calculation of the masses of neutral scalars (such as the Higgs boson) can be eliminated. For the mass computation, we specialise to the gaugeless limit and extend the effective potential computation to allow the masses to be determined without needing to solve differential equations for the loop functions — opening the door to fast, infra-red safe determinations of the Higgs mass in general theories.


JHEP12(2016)056
The discovery of the Higgs boson has added a wealth of electroweak precision observables, chief among them being its mass, which is remarkably known to within a few hundred MeV. The combined data can be used to determine the fundamental lagrangian parameters of the theory, such as the Higgs mass-squared parameter and quartic coupling, and then make predictions for future measurements, such as the Higgs self-coupling; or to provide the starting point for extrapolations of the potential to high energy scales to study unification or vacuum stability.
In the context of theories beyond the Standard Model, the Higgs mass-squared parameter and also its mass are very sensitive to new particles, and can thus be used to constrain new physics. However, precisely because of this sensitivity, the accuracy of the theoretical calculation is typically much poorer than the experimental measurement, in particular for supersymmetric field theories (which remain renormalisable up to high scales) and there has therefore been a significant effort to improve these calculations.
Typically, a subset of the scalar mass-squared parameters in the lagrangian are determined from the tadpole equations, which are the (first) derivatives of the effective potential. In the Standard Model the full effective potential was computed to two-loop order in [1], with the 3-loop leading contributions involving the strong and Yukawa couplings found in [2], and the 4-loop part at leading order in QCD in [3]. However, for general renormalisable theories the potential is only known to two-loop order 1 via the expressions given in [6] which were used in [7] to derive the tadpole equations (while the diagrams for the masses were already given in [8]).
For reasons of calculational simplicity, the effective potential beyond one loop has been calculated only in the Landau gauge, which means that the would-be Goldstone boson is treated as an actual massless Goldstone boson. Unfortunately, this leads to a technical problem known as the Goldstone Boson Catastrophe: the mass-squared lagrangian parameter of the Goldstone boson determined from the tadpole equations 2 is small and can even be negative (as opposed to the pole mass, which is always zero) and this causes the loop integrals for the tadpoles to diverge or be complex. While this problem can in principle be circumvented by dropping the complex parts and changing the renormalisation scale to attempt to find non-negligible positive mass-squareds, this is not easy to implement consistently.
A solution for the tadpoles was proposed in [9,10] for the Standard Model and applied to the MSSM in [11] (see also [12][13][14] for recent related work): (a subset of) the terms involving the Goldstone boson should be resummed to all orders, roughly speaking replacing its mass-squared parameter (which appears in the loop functions) with the equivalent parameter derived from the total effective potential (i.e. zero, since it is a Goldstone boson). In section 3 of this work we show how this can be extended to general 1 Although note that results in the zero momentum approximation are available for the Higgs boson mass in the Minimal Supersymmetric Standard Model (MSSM) up to partial three-loop order [4,5]. 2 Note that we take the expectation values to be fixed and loop-correct the mass-squared terms rather than vice versa.

JHEP12(2016)056
renormalisable theories. Note that this approach is related to the (symmetry-improved) two-particle-irreducible potential approach pioneered in [15][16][17], where essentially all particle propagators are resummed -which is somewhat more difficult to automate. In [9,10] it was noted that the Goldstone resummation would not regulate divergences in the second derivatives of the effective potential, and so to have a divergence-free calculation of the neutral scalar (i.e. Higgs) masses it would be necessary to include the external momentum in the self energies rather than using an effective potential approximation. This is particularly important because the zero momentum approximation is widely used to calculate the Higgs mass [4,5,7, -indeed there are few publicly available implementations of diagrammatic calculations of the Higgs mass beyond one loop in theories beyond the Standard Model which do not use it (some momentum-dependent diagrammatic calculations are available for the MSSM [44][45][46]). While the Goldsone Boson Catastrophe can be avoided in the MSSM in the gaugeless limit (where the Goldstone boson does not couple to the Higgs, and so generates no divergences) it is of pressing concern for more general theories, since the two-loop computation has recently become publicly available through SARAH [7,34,[47][48][49][50][51][52]; the Goldstone Boson Catastrophe as it affects that implementation was discussed in [42,53], and recently manifested itself in [54][55][56]. Indeed, while the numerical impact of the problem in the Standard Model seems to be small (at least away from the divergent points, simply neglecting the imaginary part of the potential seems to give results close to those of the full solutions), in more complicated theories it can cause divergent contributions to the masses for many regions of the parameter space; in [54,55] it was necessary to restrict to only the two-loop corrections proportional to the strong gauge coupling for those regions in performing parameter scans.
In section 4, we shall show that the inclusion of external momentum in the scalar self-energies does not by itself avoid all divergences. In fact, it is necessary to resum the Goldstone boson contributions in the mass diagrams too -to cancel the divergences in a class of diagrams which do not depend on momentum. We will also show that the resummation can be implemented most easily to two loop order by using an "on-shell" scheme for the Goldstone bosons. With these modifications, to cure the remaining divergences the diagrammatic implementation in [7] could in principle be extended to include the external momentum by changing the loop functions to those implemented in TSIL [57,58]. However, analytic expressions for general loop functions with momenta are not known: they are in general obtained by solving differential equations, which is numerically expensive. Therefore, in appendix B we give a complete set of analytic expressions for expansions of the necessary functions including all divergent and constant terms in an expansion of the four-momentum-squared s around zero (neglecting those of O(s)). This allows fast evaluation of a generalised effective potential approximation for the neutral scalar masses -although for this part we shall be restricted to the gaugeless limit (setting the couplings of all broken gauge groups to zero) since the mass diagrams are known only up to second order in the gauge couplings.
Once the Goldstone Boson Catastrophe has been solved, using similar techniques it was shown in [9][10][11] that it is also possible to improve the solution of the tadpole equations for the other mass-squared parameters (not just the one corresponding to the tree-level

JHEP12(2016)056
Goldstone boson mass). In general, the same mass-squared parameters m 2 appear both as solutions of the tadpole equations, and in the loop functions, in the schematic form where m 2 0 is the tree-level solution of the tadpole equation, v is some expectation value and ∆V are the loop corrections to the effective potential. Although resummation is not required for them (except perhaps for the Higgs boson, where the quantum corrections are so large that they force its tree-level mass to become negative -we shall not discuss such a case here), these other mass-squared parameters can be expanded perturbatively in the loop functions so that the equations can be solved directly rather than iteratively. In other words, we find only the tree-level values of the parameters on the right-hand side of the equation, and the loop-corrected solution on the left (as opposed to the loop-corrected value on both sides): we shall refer to these throughout as "self-consistent equations". In section 5, we will show how to carry out this procedure in general, showing that the formulae can be given in simpler form than in, e.g., [11] for the MSSM case. We shall also go further and show how this shifts the mass diagrams.
Finally, we have endeavoured to keep the paper as self-contained as possible, and for that purpose we provide in appendix A a set of all of the loop functions used throughout.

Abelian Goldstone model
Let us begin by recalling the problem of the Goldstone Boson Catastrophe. For simplicity we shall take the simplest abelian Goldstone model defined by a complex scalar field Φ (and no gauge group) with potential and expand around an expectation value v as Φ = 1

JHEP12(2016)056
Defining m 2 G ≡ µ 2 + λv 2 , M 2 h ≡ µ 2 + 3v 2 λ, we can then compute the effective potential up to two loops: where the one-loop functions f (x), A(x) and the two-loop function I(x, y, z) are defined in the appendix, equations (A.3), (A.6), and (A.24). The potential is regular as m G → 0 but does contain terms of order m 2 G log m 2 G (where log x is also defined in the appendix -eq. (A.2)) so that when we derive the tadpole equation and expand the derivative of The log m 2 G terms on the second line are the manifestation of the Goldstone Boson Catastrophe: we cannot insert the tree-level solution m 2 G = 0 into them, and will have a complex potential if we find m 2 G < 0. The solution proposed in [9,10] is to resum the Goldstone boson propagators -in the one-loop effective potential we make the substitution where C is a constant defined in equation (A.1), and Π GG (k 2 ) is the Goldstone boson self energy, given here at one loop by With zero external momentum, this becomes The term involving only the Goldstone mass-squared will not have a well-defined derivative, and this also leads to divergences when we resum the effective potential at three loops and above. The prescription of [10] is to drop it in favour of Π g = 1

JHEP12(2016)056
Note that this does not correspond to dropping one particular class of diagrams (at one loop it is a combination of the one-and two-propagator diagrams) but instead must be defined in terms of dropping contributions from "soft" Goldstone bosons. Defining we then should use instead the resummed potential where l is the loop order to which V eff has been calculated; the terms in square brackets simply ensure that the potentials are identical up to l loops and only differ at higher orders. Performing this procedure for the potential above we find With the above procedure, we have resummed the leading divergences at two loops, i.e. the terms of order m 2 G log m 2 G for small m 2 G (we expect m 2 G to be of order a one-loop quantity at the minimum). If we are interested in the first derivative of the potential then this is sufficient; to find the minimum to two-loop order we can expand the potential to order m 2 G with the help of eq. (A.54): (2.12) making the regularity apparent, although note that the higher-order terms still contain a m 4 G log m 2 G term. The tadpole equation, neglecting terms of three-loop order, is then Noting that the solution to the one-loop equation is m 2 G + 1 16π 2 ∆ 1 = 0, we see that we can neglect the A(m 2 G + ∆) term as it gives a correction of order three loops. We ought then to find that we can identify the term in curly brackets with ∆ 2 : for a Goldstone boson we should find m 2 G + Π GG (0) = 0, so we expect that we should be able in general to identify

JHEP12(2016)056
This leads to the prescription in [9,11], which is somewhat simpler: we expand the potential V eff as a series in m 2 G : We can then use this as the definition of ∆ 1 instead of equation (2.9). We then resum the effective potential aŝ By doing this, we immediately find the expression in (2.12), with Ω ≡ 4λ 2 v 2 R SS (0, M 2 h ). When we take the derivative and expand up to two-loop order then the minimum is at m 2 G + ∆ = 0 with We shall follow this second procedure to find the minimum condition in general renormalisable field theories at two loop order. We shall also consider a hybrid approach, which is to adopt an on-shell condition for the Goldstone boson: we define This is particularly effective at two loops, where we only need Π GG ; furthermore, since (m 2 G ) OS = 0, at this loop order there is no difference between Π GG and Π g . Making the above substitution in the potential we find exactly the same result as our resummed version in equation (2.12). However, we also have the advantage that we can make this substitution directly in the tadpole equation: which gives exactly the expression that we found above in (2.13). We shall find in the following that this simple approach is also exactly what we need for the mass diagrams. However, we must first introduce some notation and formalism to handle the general case when (potentially several) Goldstone bosons and neutral scalars can mix.

Notation for general field theories
In the previous subsection we considered the simplest possible model where there were only two real scalars which cannot mix. Once we consider more general theories, there can be JHEP12(2016)056 more Goldstone bosons and, even when they have been identified at tree level, they can in general mix with other scalars (only pseudoscalars in the case of CP conservation) once loop corrections are included. This problem does not arise in the Standard Model as treated in references [9,10], because all of the pseudoscalars are would-be Goldstone bosons and the neutral and charged Goldstones cannot mix, so can be treated as two separate sectors. In the MSSM, there are additional scalars and pseudoscalars, but in the CP-conserving case considered in [11] the mixing is at most among pairs of fields, and could be written in each case in terms of mixing angles and 2 × 2 matrices. Furthermore, the same applies for all of the scalars treated so far in those references: at most pairs of fields could mix. These complications are particularly important because in the previous cases the simplest way to derive the tadpole equations was to write down the potential and take the derivatives; once we consider more complicated cases this is no longer true and we will want to be able to directly write down expressions for the derivatives as in [7]. Starting with the scalar sector, since we will need to take the derivatives of the potential with respect to scalar fields, we shall follow the procedure in [7]. We introduce first the unrotated scalar potential in terms of real scalar fields ϕ 0 and their fluctuations around Herem 2 0,ij satisfies the tree-level tadpole equations. From this we can define the fielddependent masses and couplings, We then introduce a new basis {φ i } and an orthogonal matrixR to diagonalise the tree-level mass matrix as and obtain the new masses and couplings Next we need to define what happens when we introduce the loop corrections to the effective potential ∆V and modify the tadpole equations. We shall take the expectation values v i to be fixed (i.e. they are the true values at the minimum of the full quantumcorrected potential) and instead correct the mass-squared parameters in the lagrangian,

JHEP12(2016)056
passing fromm 2 0,ij (which satisfy the tree-level tadpole equations) to new quantities m 2 ij . Using the minimisation conditions, the relationship between them is Diagonalising these requires the introduction of a new basis via φ 0 i = R ij φ j , having masses m i and couplings λ ijk , λ ijkl .
For the couplings involving fermions and scalars, we shall use the notation for a general renormalisable field theory used in [6,8]; we repeat here the scalar, scalar-fermion and scalar-gauge-boson interactions: The fermions here are in Weyl notation and are supposed to be defined in a basis where the masses of all fields are diagonal.

Goldstone bosons in general field theories
To deal with Goldstone boson mixing in general theories, we will need some notation and simple results. We start from a theory with a global symmetry such that the scalars transform under a set of infinitesimal shifts as and differentiate the relation once: When we sit at the minimum of the potential ∂V ∂φ 0 i = 0 but for a spontaneously broken symmetry α G i is not zero for all i, and thus we have a null eigenvector of the scalar mass matrix -i.e. the Goldstone boson. For more than one symmetry broken then there will be multiple null eigenvectors and these should be formed into an orthonormal set. Let us write the symmetry shifts as linear coefficients α G i = a G ij φ 0 j after this has been performed so that i α G i α G i = δ GG and then We use the index "G" now to refer to the Goldstone boson(s) in the diagonal basis. The first identity that we need arises from taking a further derivative of the above equations to give i.e. there are no three-Goldstone couplings.

JHEP12(2016)056
If we were able to work at the true minimum of the potential and with self-consistent values of all the parameters then this would be sufficient. However, we must use the minimum conditions to determine the parameters -a subset of the mass-squared parameters, in our case -and this means that the above equations will be violated by loop corrections. In particular, the mass-squared parameter -in the diagonal basis -for the would-be Goldstone boson is no longer zero. To see this, let us define the loop tadpoles so that we can solve (2.28) with the commonly-made choice of Note that this is the value at the minimum of the potential -so δ i is not regarded as a function of {φ 0 i } when we take derivatives below. Now i.e. we can use the tree-level rotation matrices to obtain the Goldstone mass from the loop tadpoles up to corrections of two-loop order, which is all we shall require in the following. This generalises, for example, equations (2.26) and (2.27) of [11]. Following equation (2.32) above, we then see that in general. This is a crucial result in the following, even if in theories that preserve CP both couplings are zero to all orders. For theories breaking CP that could generate such a term at one or two loops, when we expand the potential as a series in m 2 G as in section 2.1 (justified by it being a one-loop quantity) we shall also implicitly expand the Goldstone self-coupling λ GG G for the same reason; implicitly because we shall not need the higherorder terms and this just corresponds to setting λ GG G = 0 everywhere. Note that this is automatic once we also employ re-expansion of the tadpoles and masses in terms of tree-level parameters to obtain consistent tadpole equations in section 5.
In practice when we are considering the broken gauge groups to be SU(2) × U(1) Y the unbroken U(1) QED allows the Goldstones to be separated into one neutral and one (complex) charged Goldstone that cannot mix. Hence in the following to simplify the notation we will restrict to a single neutral Goldstone boson and drop the lower index G, but the treatment of the charged Goldstone is identical. In this case we can also write where we now allow the normalisation of a ij v j to be arbitrary) for the linearly realised symmetries considered here.

Small m 2
G expansion of the effective potential for general theories To close this section we can now apply the notation and machinery from the previous subsections to resum the general effective potential at two loops, generalising the procedure of [9,10].

JHEP12(2016)056
The total potential up to two loops expands as (2.37) For use in the elimination of the infrared divergences in the derivatives of the effective potential, we expand V eff for small m 2 G . More precisely, we want to write the two-loop part of V eff as where the quantities ∆ 1 and Ω are to be determined. The two-loop potential splits into contributions [6]: where the subscripts denote the propagators in the loops as scalar, fermion or vector (gauge sector). The terms in the brackets will not be resummed (since they contain no scalars) and so can be taken to be unchanged from the expressions in [6]. The loop functions appearing in the other terms are recalled in the MS and DR schemes and Landau gauge in appendix A.1.2. First, the scalar contributions to the effective potential at two-loop order V and these functions can be expanded using formulae (3.7), (3.8) of [11]. Separating terms with one or more Goldstone bosons from the terms without any, and using the fact that λ GGG vanishes at leading order -see the discussion around equation (2.36) -we find the expansion of V S : from which we can identify the scalar part of ∆ 1 and Ω ,

JHEP12(2016)056
Next, the terms in V (2) involving fermions and scalars are Here, there are only two cases to consider, either k = G or k = G, and for the latter case, we can use eqs. (3.9) and (3.10) from [11] to expand the loop functions for small m 2 G . We then obtain for (∆ 1 ) F S and Ω F S Finally, the terms with scalars and gauge bosons read As previously, we can expand these terms and separate the contributions of the Goldstone boson, and we find The expansion (2.38) of V (2) enables us to rewrite the two-loop effective potential after resummation of the leading Goldstone boson contributions aŝ The minimum of this potential will be found at m 2 G + ∆ G = 0 (along with the minimisation conditions for the additional scalars) and clearly contains no logarithmic divergences for small m 2 G . The above expression could now be used for studies of general theories: the simplest would be for numerical studies where the potential is evaluated as a function of the expectation values and the derivatives taken numerically, as performed for the MSSM in [28,29] and implemented generally in [34]. However, there are potential numerical instabilities JHEP12(2016)056 when the expectation values of additional scalars are small, and for complicated models many evaluations of the potential are required which can be slow: it is therefore useful to have explicit expressions for the tadpoles, as were derived at two loops in [7]. In the next section we shall compute these for the resummed potential.

Removing infra-red divergences in the minimum condition
In the previous section we derived the resummed two-loop effective potential expanded in m 2 G that explicitly contains no infra-red divergences in its derivatives. In this section we shall present these derivatives. However, we shall also present a new approach to the problem which allows us to calculate the derivatives simply, and so we shall also give our derivations. For the scalar-only diagrams we do this by three methods: (i) The first method is to generalise the approach of [9,11], and simply take the derivatives of the resummed potential (2.53). However, this has the disadvantage of requiring us to compute the derivative of the rotation matrix elements (ii) To avoid the derivatives of rotation matrix elements, we instead take the derivatives ofV eff before diagonalising the mass matrix and singling out the Goldstone boson and expanding the potential in m 2 G . This leads to a simpler derivation of the results.
(iii) For our third method, we introduce a new approach: we set the Goldstone boson mass "on-shell" in the (non-resummed) effective potential. We shall show that this gives the same result as the other methods but (much) more simply, and does not suffer from the problem of needing to exclude Goldstone self interactions by hand. Furthermore, in the next section we shall employ this approach to compute the mass digrams, which would be more complicated using the alternative methods.

Elimination of the divergences by method (i)
Generalising the approach of [11] to extract the tadpoles we take the derivatives of equation (2.42). Starting with the one-loop potential, we note that, since m 2 G + ∆ G = 0 at the minimum, the derivative of f (m 2 G + ∆ G ) will vanish. Hence we only require Note that throughout we shall adopt the Einstein convention for summing repeated indices when all indices are to be summed over; when there is an index that is summed over only a subset (i.e. excluding the Goldstone boson indices) we shall write an explicit sum symbol.

JHEP12(2016)056
For the two-loop terms, recall the scalar part Treating each of these pieces in turn we find: and similarly we see Putting this all together we see that they combine to give the compact expression .
Next we turn to the SS terms: where the two terms again combine into a single compact expression. The final piece is

Elimination of the divergences by method (ii)
From inspection it is clear that the one-loop tadpole is not divergent when we send m 2 G → 0. However, at two loops we found that the process of isolating the divergences in the potential, expanding it in the Goldstone mass, and then taking the derivatives was rather cumbersome due to the derivatives of the mixing matrix elements R ij . Instead we could consider taking the derivatives before having cancelled out the divergent parts, and then ensure the cancellations later. Hence we rewrite the resummed effective potential aŝ using formulae (2.38) and (2.53). We expect the terms from the derivative of − 1 2 A(m 2 G )∆ G to cancel off the IR divergences in the derivatives of V eff . To show this, we use the expression of ∆ G derived in eq. (2.43). The relevant contribution to the minimum condition at twoloop order is The purely scalar contribution to the non-resummed tadpoles is, at one-loop order and at two loops where [7] T p SS = 14) with the notation f

JHEP12(2016)056
In these formulae, we can then consider separately the Goldstone contributions and investigate the divergent terms. We find two types of divergent terms in eq. (3.12): • The first type of divergent term comes from T SS , for j = k = G, and 3 l = G, and reads • The other divergent terms, coming from T p SSSS with j = j = G, are • A potentially more dangerous element of those terms, for the particular case k = l = G is not present as λ GGG = 0 (at least up to terms of one-loop order).
All the other terms in After relabelling of the indices in the sums, we observe that the log m 2 G divergences from the terms in eqs. (3.16) and (3.17) cancel out perfectly with the ones from eq. (3.10). We can then take the limit m 2 G → 0 in the one-loop and two-loop parts of the minimum condition: this limit is regular in the one-loop tadpole (3.11) so we recover eq. (3.1), while we find at two-loop order. It is important to notice that all three functions f SSS , U 0 and R SS are regular when one of their arguments goes to zero, hence the result we find is indeed free of infrared divergences.

Elimination of the divergences by setting the Goldstone boson on-shell
Here we shall introduce a new approach to the Golstone Boson Catastrophe: we shall treat the Goldstone boson mass as an on-shell parameter and enforce that it is identically zero. This means replacing the dimensionally regularised (DR or MS) Goldstone mass by the on-shell (or pole) mass in the following way where the pole mass is (m 2 G ) OS = 0. Note that we only need the one-loop relation here, so any mixing in the mass terms between the Goldstone boson and other (pseudo-)scalars is irrelevant -it would be proportional to (Π (1) iG ) 2 and thus a two-loop effect. When we write the effective potential in terms of the on-shell Goldstone boson mass we should find that it is free of divergences. To do this, we shall start from the dimensionally regularised potential and substitute the Goldstone boson mass in equation (3.19), expanding out to the appropriate loop order; this gives the result that we would obtain by performing the calculation using the on-shell mass with the appropriate counterterms. For our case, we only need to use the one-loop self-energy in the one-loop tadpole; the scalar contribution to the Goldstone boson self-energy at one-loop order is where we again require the result λ GGG = 0 to leading order -although in this case we could (if desired) make it an on-shell condition. Applying the above relation to the tadpole in eq. (3.11) we obtain the following shift to the two-loop tadpole: Since B(0, m 2 j , m 2 k ) = −P SS (m 2 j , m 2 k ), these shifts correspond exactly to the divergent terms we saw in equations (3.16) and (3.17) and so when we formally take the limit (m 2 G ) OS → 0 we find exactly the same tadpole given explicitly in (3.18) that we found by the two other methods. This derivation is certainly much faster than the first method, but note that the principle is different to the previous calculations: there is no ad-hoc resummation, nor are we required to expand the potential as a series in m 2 G . However, perhaps remarkably, we find exactly the same result for the tadpole that remains, implying that, at least at two loops, the two approaches are equivalent. This new approach will prove to be simpler than both previous methods when we turn our attention to mass diagrams; for now we shall simply complete the set of tadpole equations.
Before moving on to diagrams with fermions, we shall comment on the prescription to follow when there is more than one Goldstone boson. In that case, since the Goldstone

JHEP12(2016)056
bosons are all degenerate the mutual mixing between them becomes a leading-order effect and we must diagonalise the self-energies Π GG on the subspace of indices G, G which run over all Goldstones. However, we can also easily write this in the non-diagonalised basis as a generalisation of (3.19): where formally all Goldstone bosons have the same mass m 2 G which we set to zero. Then we can rewrite the tadpole as If the gauge group of the model of interest is just that of the Standard Model, then clearly the charged and neutral Goldstone bosons cannot mix, so this becomes trivial -hence in the following we shall restrict for clarity to the one-Goldstone case. However, we shall later write the full result in the general case.

Diagrams with scalars and fermions
The one-loop tadpoles involving fermions are and these do not present any divergence in the limit of vanishing Goldstone boson mass. The two-loop contributions are [7], (3.27) with the loop functions from eq. (II.38) of [7]. The second term T F F F S is regular when m 2 G → 0, because the loop functions, B 0 , I, U 0 , that appear in its expression are all regular when only one of their argument goes to JHEP12(2016)056 zero. However, the k = l = G terms in T SSF F are divergent: After either resummation or setting the Goldstone boson on-shell we find the total, finite, two-loop contribution T p SSF F in equation (3.41) and note that T p F F F S is not modified from eq. (3.27).

Diagrams with scalars and gauge bosons
The one-loop tadpole involving (massive) gauge bosons is (3.29) which contains no scalar propagators so has no divergences in the Goldstone boson mass. However, the gauge boson contribution to the one-loop scalar self-energy in Landau gauge is [8]: where the loop functions are given in [8,57] but simplify for zero momentum in Landau gauge to Recall that there are six scalar-gauge boson contributions to the two-loop tadpole [7]: Of these only three are potentially singular -f

JHEP12(2016)056
i.e. they correspond exactly to the potentially singular terms. However, note that B SV term is zero -and indeed we find that f We give the final finite tadpoles in equation (3.43).

Total tadpole
Here we gather the results of the previous subsections and rewrite them for the most general case, that of multiple Goldstone bosons. The total tadpole, after curing the Goldstone boson catastrophe and taking m 2 G → 0 everywhere, is The all-scalar diagrams are λ pGG λ Gkl λ G kl R SS (m 2 k , m 2 l ), (3.40) where by (j, j ) = (G, G ) we mean that j, j are not both Goldstone indices. The fermionscalar diagrams are

JHEP12(2016)056
while the gauge boson-scalar tadpoles are Finally the gauge boson-fermion and gauge diagrams are not affected by the Goldstone boson catastrophe, as scalar masses do not appear in them, and can be found in the appendix C.2 of [7] T

Mass diagrams in the gaugeless limit
As discussed in the introduction, the scalar masses are among the most interesting electroweak precision observables, and their calculation also suffers from the Goldstone Boson Catastrophe. Earlier literature pointed out that the calculation in the effective potential approximation contains more severe divergences that cannot be solved by resummation, and thus the inclusion of the external momentum is necessary. However, we shall find that there are also divergences that are not regulated by external momentum -and thus both setting the Goldstone boson on-shell and external momentum are required to obtain finite, accurate results. On the other hand, the effective potential approximation is still useful and has advantages over a full momentum-dependent result, chief among these being simplicity and speed of calculation. In particular, the evaluation of the loop functions at arbitrary external momentum requires the numerical solution of differential equations [57] which, although implemented in the fast package TSIL [58], is still much slower than the zero-momentum functions, and when the functions must be repeatedly called can lead to times orders of magnitude longer for complicated models. Hence we shall consider expanding the two-loop JHEP12(2016)056 self-energies as a series in s ≡ −p 2 (for metric signature (−, +, +, +)) as s k k! (4.1) and we shall neglect terms of O(s), giving a "generalised effective potential" approximation: for loop functions where the singular terms Π l,ij vanish the result is identical to the second derivative of the effective potential. This approximation is particularly good when the mass of the scalars considered is smaller than the scale of other particles that they couple to; but even when they are similar we find that typically the difference is only a few percent. This should then be within other uncertainties in the calculation for most purposes.
We shall perform our calculations using our procedure of taking the Goldstone boson mass(es) on-shell as before, working in the general case now of allowing multiple Goldstone bosons throughout. We shall make heavy use of the existing expressions for two-loop scalar self energies from [8]; however, these are only available up to second order in the gauge coupling. Hence we shall be restricted to work in the very popular "gaugeless limit" where we neglect the gauge couplings of broken gauge groups (including electromagnetism, since hypercharge and weak SU(2) are both broken so their gauge couplings are neglected). The two-loop self-energy in this limit can be decomposed as follows: This consists of scalar-only propagators, diagrams with scalar and fermion propagators, diagrams with scalar and vector propagators, and fermions and vectors. We find that Π

All-scalar terms
The two-loop scalar self-energy contribution with only scalar propagators is given by [8]: The loop integral functions are recalled in (A.52).

JHEP12(2016)056
When at most one of the propagators is a Goldstone boson, we can set m 2 G → 0, s → 0 and use the simplified expressions below (B.2) of [7]. However, for cases including more Goldstone bosons we must look for singularities since, in general, only the S SSS term is regular. Furthermore, we can divide the functions into those regulated by the momentum and those that are not. In particular, by inspection we see that for two or more Goldstone bosons W, X, Y, V can be divergent as m 2 G → 0, even for finite momentum; this means those terms must be regulated by resummation -or, in our case, by shifts from the oneloop self energy by putting the Goldstone bosons on shell. On the other hand, the terms U, M and Z must be regulated by including finite momentum.
It should be noted that the divergences that are not regulated by momentum all involve a Goldstone boson self-energy as a subdiagram. It is then logical to consider how they relate to the divergent terms in the tadpole graphs. If we consider the effective potential approximation and take the derivatives of the tadpoles as in [7], then we see that the topologies X, Y, Z descend from the T SS graphs; S, U arise from T SSS ; and M, V, W from T SSSS . Then it is clear that, since the T SSS graphs contain no divergences, resummation is irrelevant for S and U , while T SS and T SSSS are both divergent when there is part of a Goldstone boson self-energy as a subdiagram. We also see that W and X topologies arise from T SSSS and T SS respectively by replacing a three-point vertex with a four-point one, and likewise V and Y arise by adding a leg connected directly by a propagator to the other leg; we illustrate this whole discussion in figure 1. Hence we expect that these special divergences should follow the same pattern as the tadpoles, and be cured in the same way. However, we shall also find below some subtleties remain in the V topology.

Goldstone shifts
To determine the effect on the mass diagrams, let us make the shifts using the method of an on-shell Goldstone boson. Recall that the contribution to the one-loop self-energy is where B is defined in eq. (A.13), and are all regular, the diagrams V SSSSS (0, 0, x, y, z) and V SSSSS (0, 0, x, 0, y) are divergent, and their expression may be found simply by using those for U (0, 0, x, y) and U (0, 0, 0, x) given in appendix B: (4.14) All other V SSSSS diagrams are either regular or vanish due to the prefactor λ GGG . The remining functions U SSSS , M SSSSS and Z SSSS require regulation by momentum: we give expressions for the expansion of these in appendix B.

Fermion-scalar diagrams
The potentially singular mass diagrams are Π , but among these there are only a subset once more that are regulated by the Goldstone boson shifts; indeed, as in the purely scalar case we find that the topology M is purely regulated by momentum for which all of the limits of the loop functions are provided in appendix B. For the other two, there are exactly four diagrams to regulate, which will match exactly. They have the form [8]: As in the scalar case, we look at the shift in the one-loop scalar mass contribution involving Goldstone bosons: and we can expand as usual in perturbation theory using ∆ ij = 1 ij + . . . to find that we should shift the tadpoles according to By (i, i ) = (G, G ) we mean that the sum over (i, i ) excludes the cases where both i and i are Goldstone boson indices. This allows us to express the δ i entirely in terms of the treelevelm 2 parameters and obtain a perturbative expansion for m 2 -note that we should also replace all of the couplings λ ijk , λ ijkl etc and rotation matrices R ij with their tree-level valuesλ ijk ,λ ijkl ,R ij (we already implicitly used this to disregard the λ GGG terms). The only subtlety occurs whenm 2 i =m 2 j for some i, j which is not ensured by a symmetry so that ∆ ii = ∆ jj ; in that case as usual theR matrix must be modified to diagonalise ∆ (1) ij on those indices. However the expression above is still valid in that case. Note that the shift only occurs for scalar propagators in the one-loop diagrams, which is why there is no modification of the fermionic or vector tadpole diagrams.
We can apply the same procedure to use the tree-level masses in the mass diagrams: after some algebra we find (in the gaugeless limit -otherwise we will have some additional shifts from scalar-vector diagrams) that Π (2) ij (s, m 2 ) = Π where we used the usual C function defined in eq. (A.17). These together then allow us to determine the scalar masses to be the values of s that give solutions to: Typically in spectrum generators the two-loop corrections are computed at fixed momentum and then the eigenvalues of the above matrix computed iteratively. Since we have given the expansion of all the loop functions relevant for the two-loop corrections up to terms of order

JHEP12(2016)056
O(s), this could be generalised to include our simple momentum dependence for the twoloop part as in equation (4.1) without significant loss of speed since the computationally expensive parts of the two-loop functions would only need to be evaluated once. However, since all of the expansions are strictly valid only up to two-loop order, the equation above could be solved perturbatively itself with no significant loss of accuracy.

Conclusions
We have presented a solution to the Goldstone Boson Catastrophe in general renormalisable theories to two-loop order. We showed that the approach of Goldstone boson resummation is equivalent (at least at two-loop order) to an on-shell scheme for the Goldstone boson(s), the latter being much more convenient calculationally. We then showed how there are a set of self-energy diagrams that also exhibit the Goldstone Boson Catastrophe even when external momentum is included -but that our solution naturally avoids those singularities. We were then able to give expressions for a "generalised effective potential approximation" for neutral scalar masses in the gaugeless limit, that are free of infra-red divergences and give a good approximation to the full momentum-dependent result. This also included the re-expansion of the masses in terms of the values obtained from the tree-level tadpole equations, allowing a self-consistent solution of the tadpole equations (i.e. equations where no terms to be solved for appear on both left and right hand sides). The expressions contained in this paper should now allow simple infra-red safe calculations in a wide variety of theories. Most practically, it would be simple to implement them in a package such as SARAH, to enable automated calculations for any model and avoid the problems seen, for example, in [42,53,56], with the existing implementation. This should also enable more reliable and accurate explorations of the parameter space of many models; in particular for non-supersymmetric models (such as the two-Higgs-doublet model), where the existing "solution" 4 to the Goldstone Boson Catastrophe is not particularly successful, relying as it does on there being a gauge-coupling dependent part of the scalar potential (as in supersymmetric theories).
However, it would also be interesting to explore further many aspects of the problem more generally: the two-loop mass-diagram calculation to quartic order in the gauge couplings; the link between resummation and our on-shell scheme; and also the extension to higher orders. Indeed, these three topics are linked: in [10], it was shown that the momentum dependence of the self-energy in the resummation was necessary for the resummation of certain subleading divergences. By reorganising the expansion in terms ofM 2 G ≡ m 2 G + Π g (0), they showed that the one-loop resummed potential (2.5) -which contributes the most divergent parts -can be rewritten as

JHEP12(2016)056
This shows that the momentum dependence cannot contribute a more divergent term than O(M 4 G logM 2 G ). Hence if we rewrite the diagrams in terms ofM 2 G -similar to our onshell scheme -then the momentum dependence of the self-energy will disappear from the tadpole condition; indeed, [9] did not require momentum dependence. This also shows that for higher-order contributions it may be most efficient to perform the calculations directly in such a scheme, rather than work in a pure minimal subtraction scheme and then apply the shifts. On the other hand, it could be relevant for the mass diagrams: [10] showed that a term of orderM 4 G logM 2 G in the two-loop potential arising from a diagram with a W-boson, charged Goldstone boson and photon could be resummed by including the momentum dependence in the self-energy; the masslessness of the photon giving rise to additional infra-red divergences. This issue did not arise here because we worked up to two-loop order, and in the gaugeless limit for the mass diagrams (so that the photon cannot contribute, and its contribution to the tadpoles is benign). It would certainly be interesting to explore this in the future. Acknowledgments M.D.G. thanks Florian Staub for collaboration and discussion on related topics. We both thank Pietro Slavich for many interesting discussions, and for very detailed comments on the first version of this paper. This work was supported in part by French state funds managed by the Agence Nationale de la Recherche (ANR), in the context of the LABEX ILP (ANR-11-IDEX-0004-02, ANR-10-LABX-63), and M.D.G. acknowledges support from the ANR grant "HiggsAutomator" (ANR-15-CE31-0002). J.B. was supported by a scholarship from the Fondation CFM.

A Loop functions
Throughout our work, we have followed closely the notations of [57], however we present in this appendix the loop functions and the notations that were used. These definitions of loop functions use Euclidean momentum integrals in dimensional reduction to d = 4 − 2 dimensions, and involve the loop factor We also recall the following shorthand notations where Q 2 = 4πe −γ E µ 2 is the renormalisation scale squared.

A.1.1 One-loop functions
In the expression of the one-loop effective potential, we make use of the function f defined as

JHEP12(2016)056
Two important one-loop functions that will appear in the expression of the effective potential, of its derivatives and in the self-energies are the finite parts of ) and x ± = p 2 + x + y ± (p 2 + x + y) 2 − 4p 2 x 2p 2 . (A.9) In two-loop order expressions, the function J is sometimes used, although it is equal to A J(x) = A(x). (A.10) A limit of particular interest of B is the limit of vanishing external momentum, that we denote B 0 , and is related to the P SS function we have used The derivative of the B function with respect to one of the mass arguments is also used, with the notation B (p 2 , x, y) = ∂ ∂x B(p 2 , x, y). (A.13) For the fermion and gauge boson contributions to the scalar self-energy we also use the functions P F F , P F F and P V V related to A and P SS as In the context of the reexpansion of the mass diagrams, we also make use of the oneloop three-point function C(p 2 1 , p 2 2 , (p 1 +p 2 ) 2 , x, y, z), which is the finite part of the following integral C(p 2 1 , p 2 2 , (p 1 + p 2 ) 2 , x, y, z) ≡ −C d d k (k 2 +x)((k−p 1 ) 2 +y)((k−p 1 −p 2 ) 2 + z) . R V V (x, y) = 1 4xy(x − y) 2 3A(x)A(y) x 2 + y 2 + 6xy − 24xy xA(x) + yA(y) + 14xy(x 2 + y 2 ) + 20x 2 y 2 + 3(x + y) 3 I(0, x, y) − 3(x − y) 2 xI(0, 0, x) + yI(0, 0, y) . (A.66)

JHEP12(2016)056
taken from expanding expressions in [57,59] or found by newly solving or expanding the integral equations in [57]. Hence we stress that (most of) this section contains new results not found elsewhere. First, for Z, we only need the fact that Then, for the U function, taking one argument to zero does not cause any divergence, and we find, looking at the integral definition (A.19) of U , that U (x, y, 0, 0), U (0, x, y, 0), U (x, 0, y, 0), U (x, 0, 0, 0) and U (0, y, 0, 0) are all regular so we can substitute them for U 0 + O(s). The only divergent function is U (0, 0, x, y) that has the form is not required, as it always appears with λ GGG as a factor, which is zero up to higher order corrections. Turning now to the M function, there are more cases to consider. In the case of only one argument vanishing, we see from the integral expression (A.20) that the function is regular. From eqs. (6.28) and (6.31) in [57], we also find that M (x, y, z, 0, 0), M (x, y, 0, 0, v) and M (x, y, 0, 0, 0) are finite.
The expression of M (x, 0, 0, 0, 0) with full momentum dependence can be found in equation (6.31) of [57], and becomes when expanding to leading order for small s M (x, 0, 0, 0, 0) = 1 6x 18 + π 2 − 12 log (−s/x) + 3 log 2 (−s/x) . The approximate formulae for U (0, 0, x, y), M (0, y, 0, u, v) and M (0, 0, 0, 0, v) have been checked against the numerical results from TSIL [58] and show excellent agreement until s becomes of the order of the arguments in the functions -even when s is of the order of the mass parameters, the difference between the approximate result and the numerical from TSIL is about 10%.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.