Two-loop analysis of classically scale-invariant models with extended Higgs sectors

We present the first explicit calculation of leading two-loop corrections to the Higgs trilinear coupling λhhh in models with classical scale invariance (CSI), using the effective-potential approximation. Furthermore, we also study — for the first time at two loops — the relation that appears between the masses of all states in CSI theories, due to the requirement of reproducing correctly the 125-GeV Higgs-boson mass. In addition to obtaining analytic results for general CSI models, we consider two particular examples of Beyond-the-Standard-Model theories with extended Higgs sectors, namely an N-scalar model (endowed with a global O(N) symmetry) and a CSI version of the Two-Higgs-Doublet Model, and we perform detailed numerical studies of these scenarios. While at one loop the value of the Higgs trilinear coupling is identical in all CSI models, and deviates by approximately 82% from the (one-loop) SM prediction, we find that the inclusion of two- loop corrections lifts this universality and allows distinguishing different BSM scenarios with CSI. Taking into account constraints from perturbative unitarity and the relation among masses, we find for both types of scenarios we consider that at two loops λhhh deviates from its SM prediction by 100 ± 10% — i.e. a quite significant further deviation with respect to the one-loop result of ∼ 82%.


Introduction
The discovery of a Higgs boson of mass 125 GeV at the CERN Large Hadron Collider [1,2] has completed the particle spectrum of the Standard Model (SM) of particle physics, and has established the Higgs sector as the origin of electroweak symmetry breaking (EWSB).
Nevertheless, numerous open questions remain unsolved by the SM, and point to the need for Beyond-the-Standard-Model (BSM) physics. Among these issues, one can mention, on the theory side, for instance the hierarchy problem, the lack of a description of quantum gravity, or inflation. Additionally, a number of experimental results have contradicted prior assumptions -upon which the SM was built -such as for example (to cite only a few) the need to realise baryogenesis, the evidence for dark matter and dark energy, or the discovery that neutrinos have tiny masses.
There is no doubt that the Higgs boson plays a central role in new physics, as well as in searches for it. Indeed, many (if not most) problems of the SM can be related to Higgs physics. In order to address these, BSM theories commonly feature a Higgs sector with extended particle content or modified dynamics -as there is no compelling argument for it to be minimal. Therefore, shedding light on the structure of the Higgs sector is one of the most pressing tasks for particle physicists. Unfortunately, for the time being, all measurements of Higgs properties seem to be in excellent agreement with SM predictions [3,4], and this forces us to investigate how new physics can appear to solve problems of the SM while simultaneously evading detection. A first, but unexciting, possibility is decoupling -i.e. that all new states are well beyond our experimental reach and their contributions to SM(-like) observables are vanishingly small due to the decoupling theorem [5]. Another more interesting option is alignment without decoupling [6], which occurs 1 when the couplings of the 125-GeV Higgs boson to fermions and gauge bosons are JHEP03(2021)297 SM-like at tree level, possibly because of some mechanism or symmetry of the Lagrangian -see for instance refs. [7][8][9][10][11][12][13][14]. We will return to this latter case in the following. 2 Among the issues of the SM, one that has received considerable attention is the (gauge) hierarchy problem [16], or in other words the need for a mechanism to stabilise the mass of the Higgs boson M h at the electroweak (EW) scale. Because the Higgs boson is a scalar, its mass is in the SM not protected by any symmetry, meaning that it typically receives huge -quadratically divergent -radiative corrections. In turn, this implies that a cancellation of tremendous -and unrealistic -accuracy would have to occur between tree-level and loop-level contributions to M h . Solutions to the hierarchy problem usually involve some new symmetry or mechanism that protects the Higgs mass. For example, supersymmetry [17] introduces a new symmetry between bosons and fermions, so that scalar masses become protected by the chiral symmetry. Meanwhile, in composite Higgs models [18,19], there are no fundamental scalars and the Higgs boson appears as a composite state of elementary fermions. But because new physics has so far eluded discovery, experimental constraints on most popular BSM scenarios have become more stringent (see for instance refs. [20,21] and references therein) and the need to explore new avenues of research is increasingly evident.
Another interesting possibility to address the hierarchy problem, which we will consider in this paper, comes from the concept of classical scale invariance (CSI), first proposed by Bardeen in ref. [22]. The idea of CSI is to forbid all mass terms in the Lagrangian, so that the theory is scale invariant at the tree (i.e. classical) level. Of course, as this is only imposed at tree level, scale invariance is explicitly broken by radiative corrections, and EWSB can occur radiatively, via the Coleman-Weinberg mechanism [23] -which was extended to multi-scalar cases by Gildener and Weinberg in ref. [24]. The Higgs vacuum expectation value (VEV) v and all particle masses (proportional to v) are then generated via dimensional transmutation. If the CSI model is taken to be UV complete, or in other words if scale invariance is assumed at the Planck scale, the absence of a tree-level Higgs mass parameter together with the generation of the Higgs mass entirely at loop level can ensure that there is no hierarchy problem -see for instance refs. [22,[25][26][27][28][29][30][31]. On the other hand, models in which scale invariance is imposed at the EW scale can often not be extended up to the Planck scale, because of Landau poles appearing in the running of dimensionless parameters. However, even if a Landau pole appears at some intermediate scale Λ and thus requires the existence of some further new physics beyond Λ, tree-level scalar masses are prohibited up to this scale thanks to the CSI. Therefore there is no hierarchy problem up to the scale Λ, and a problem only remains between Λ and the Planck scale -arguably a less severe issue than in the SM. The hierarchy problem is thus not entirely solved in such theories, although it can be somewhat alleviated.
Nonetheless, motivated by a number of phenomenological considerations, we restrict our attention in this paper on the latter type of CSI models, defined with extended scalar sectors at the EW scale. A first interesting feature of such models is that, in the absence of any (BSM) mass term, new states cannot be decoupled, nor very heavy, and hence should JHEP03(2021)297 be found in a foreseeable future -if they exist. Additionally, in CSI models, the 125-GeV Higgs boson, which corresponds in the language of ref. [24] to the "scalon" direction along which EWSB is radiatively realised, is automatically aligned at tree level -as pointed out in ref. [32] for a CSI variant of the Two-Higgs-Doublet Model (2HDM) -thereby explaining current experimental results.
A further discriminative characteristic of CSI models comes from the Higgs trilinear coupling λ hhh . This coupling is of great importance as it controls the shape of the Higgs potential, and in turn, the nature and strength of the EW phase transition (EWPT). In particular, one of the requirements to make the scenario of electroweak baryogenesis work [33][34][35] is for the EWPT to be of strong first order. In turn, a large deviation of the Higgs trilinear coupling from its SM prediction is one of the main signs of a strong firstorder EWPT -for massive models, a deviation of at least 20% was found to be necessary in refs. [36,37] -together with the special spectrum of gravitational waves that would be produced by bubble collisions during the EWPT [38].
Couplings of the Higgs boson -and especially its trilinear coupling -can potentially differ significantly from their SM values in BSM theories with extended scalar sectors, because of non-decoupling effects from loop corrections involving the BSM states -this was found at one loop, first in the 2HDM in refs. [39,40], and later in various BSM models with extra singlets [41][42][43][44], doublets [44][45][46][47], or triplets [48]. The limit when these BSM corrections are maximal is when the BSM scalars acquire large masses exclusively from the Higgs VEV (rather from some BSM mass parameter). Incidently, this is exactly what happens in CSI models, and as a matter of fact, the Higgs trilinear coupling has been found [49][50][51] to be at one loop (λ CSI hhh ) (1) = 5/3(λ SM hhh ) (0) in any CSI model (without mixing), where (λ SM hhh ) (0) is the tree-level SM prediction for λ hhh -in other words (λ CSI hhh ) (1) deviates by ∼ 67% from (λ SM hhh ) (0) . The universality of this prediction is a particularly strong and unique property of CSI models, and follows (as we will see in section 2) from the especially simple form that the effective potential takes at one-loop order along the "scalon" direction. This value of (λ CSI hhh ) (1) also means that the EWPT is always of strong first order in CSI theories -this, together with the synergy between the measurements of gravitational waves and of the Higgs trilinear coupling was investigated in ref. [52] (for the example of O(N )-symmetric CSI models).
One may then naturally wonder what is the current experimental status for the determination of λ hhh as well as future perspectives. It turns out that large deviations of λ hhh from its SM prediction are still allowed by experiments: indeed, at the present time the most stringent limits on the trilinear coupling come from the ATLAS collaboration using single-Higgs production searches with LHC Run 2 data and give −3.2 < λ hhh /λ SM hhh < 11.9 at 95% confidence level (CL) [53] (see also refs. [54,55] and [56,57] respectively for ATLAS and CMS results using double-Higgs production). However, these limits will be greatly improved at future colliders, thus making λ hhh an ideal target to search for BSM effects. We refer the interested reader to ref. [58] for a detailed discussion of the current prospects for the determination of λ hhh at upcoming experiments, and we only summarise here some of the main results. First of all, the high-luminosity upgrade of the LHC (HL-LHC) is expected to reach an accuracy of about 50 − 60% with 3 ab −1 of data [59] (see also ref. [60])
The universality of the value of λ hhh in CSI models at one loop also raises the question of the impact of two-loop corrections. Indeed, as the form of the effective potential is changed at two loops -with new squared-logarithmic terms appearing, see e.g. refs. [68][69][70] -one may expect the effective Higgs trilinear coupling derived from it to be modified as well. Moreover, two-loop corrections to λ hhh have already been investigated in a number of extension of the SM with enlarged scalar sectors, 4 namely in the Inert Doublet Model [73][74][75], in the 2HDM [74,75], and in a real-singlet extension of the SM [75] (this last model actually corresponds to the N -scalar theory we will consider in this work for N = 1). It is therefore also important to compute two-loop corrections to λ hhh in the CSI theories, to match the accuracy achieved in non-CSI scenarios and allow precise comparisons, as the computational method differs between the two types of models.
In this paper, we have performed the first explicit 5 calculation of dominant two-loop corrections to λ hhh in theories with CSI: we extend the known one-loop (effective-potential) calculation to two loops, so as to include the same level of corrections as in our previous works on non-CSI models [74,75]. We investigate in particular two types of BSM scenarios: first, O(N )-symmetric CSI models, and second, a CSI variant of the 2HDM [76]. Moreover, we include the renormalisation scheme conversion, from the MS scheme (in which the effective potential is obtained) to the OS scheme, in order to express our results in terms of physical quantities -note that we will be referring to the MS and OS results respectively as λ hhh andλ hhh to avoid confusions. We also take into account the requirement of generating the 125-GeV mass of the Higgs boson correctly, at the two-loop order. This yields a relation between the masses of all the states in the CSI models, which we use to constrain the BSM parameters (masses and couplings) appearing in our calculations. We also ensure that perturbative unitarity is maintained and that the EW vacuum remains the true vacuum of the effective potential for the parameter points we consider. As mentioned already, the Higgs trilinear coupling has a universal prediction at one loop in all CSI models -(λ CSI hhh ) (1) = 5/3(λ SM hhh ) (0) -and our main result in this paper is finding that the inclusion of two-loop corrections breaks this universality. Through our detailed numerical studies,

JHEP03(2021)297
we observe that the different theoretical constraints we impose on the considered scenarios limit the allowed range of values for λ hhh at two loops -the strongest constraint coming from the requirement of correctly reproducing the 125-GeV Higgs mass. Nevertheless, the two-loop Higgs trilinear coupling retains dependence on the parameters of the BSM scalar sectors, and may provide an opportunity to distinguish different BSM scenarios with CSI. This paper is organised as follows: in section 2 we present general results, applicable to a wide range of CSI models, and we illustrate our calculational scheme with the example of the CSI Standard Model. Next, in sections 3 and 4, we study the Higgs trilinear coupling and the relation between masses in scenarios of N -scalar models with O(N ) global symmetries and also in a CSI variant of the 2HDM. We discuss implications of our present in section 5, before concluding in section 6. A number of appendices provide details about our computations and some results we deemed too long for the main text. Definitions of one-and two-loop functions are given in appendix A, complementary expressions for the O(N )-symmetric models are provided in appendix B, and results in the CSI-2HDM for general masses are presented in appendix C. Finally, we provide in appendix D generic results for the coefficients in the effective potential that can be applied for any CSI model without mixing.

Generalities about classically scale-invariant theories
We begin this paper by recalling known features of CSI models, before deriving some new results applicable for all such models.

Classically scale-invariant theories
Classically scale-invariant models arise from the assumption that at some energy scalewhich we will in this paper consider to be at or near the EW scale -the scalar potential exhibits scale invariance, by which we mean that all mass-dimensionful quantities vanish at the classical (i.e. tree) level. Radiative corrections explicitly violate the CSI and allow the EW gauge symmetry to be broken dynamically, giving rise to a mass scale via dimensional transmutation. This mechanism to realise spontaneous EWSB radiatively was first studied for a model of a single scalar by S. Coleman and E. Weinberg [23]. However, it is known that it is impossible to explain the Higgs boson mass in this simplest single-scalar case: neglecting the top-quark effects, the original paper of Coleman and Weinberg [23] predicted a scalar of mass less than 10 GeV, and including the top quark renders the (one-loop) potential unstable. In order to reproduce the 125-GeV Higgs mass in models with CSI, it is therefore necessary to consider extended scalar sectors. The treatment of this general case of multi-scalar theories requires more care as one must also identify the direction in field space along which the EW gauge symmetry is broken. This is the purpose of the method devised by E. Gildener and S. Weinberg [24], which we will review briefly now.

The Gildener-Weinberg method
The tree-level scalar potential of a generic CSI theory with scalars {ϕ i } takes the form

JHEP03(2021)297
where the Λ ijkl are (dimensionless) quartic couplings -tadpole, mass, and trilinear terms all have a non-zero mass dimension and must therefore vanish. Clearly, this potential does not admit any other minimum than the origin {ϕ i } = 0 and in order for the EW gauge symmetry to be broken spontaneously using the Coleman-Weinberg mechanism, it is necessary that the tree-level potential vanish along a ray in field space. This particular flat direction corresponds to the SM-like Higgs boson ϕ h (called scalon in ref. [24]) which carries the EW vacuum expectation value (VEV) v and acquires a mass of 125 GeV at loop level. Along the flat direction, the scalar potential is generated at one-loop order, and it can be written as a function of the order parameter (i.e. the Higgs field) h ≡ ϕ h − v using the supertrace formula [77] where M S,f,V (h) denote respectively the scalar, fermion, and gauge-boson field-dependent mass matrices (at tree level), and Q is the renormalisation scale. Note that, this expression is given in the MS renormalisation scheme. Additionally, compared to the paper of Gildener and Weinberg [24], we have kept here all the terms of the supertrace formula, rather than just the logarithmic terms. As we will see in the next section, this has no impact on any physical quantity or discussion. Due to the absence of any mass term in CSI models, the field-dependent masses of all particles in theory -scalars, fermions, gauge bosons -can be related to the fieldindependent masses as We can exploit this particularity of CSI models to write the one-loop effective potential under the simple form where κ is the loop factor defined in eq. (A.1) and (2.5) We will discuss in the following how two-loop corrections modify this form of V eff , however, first we will review the one-loop result for the effective Higgs trilinear coupling that is computed with this potential.

The Higgs trilinear coupling at one-loop order
We study in this paper corrections to an effective Higgs trilinear coupling, calculated as the third derivative of the effective potential with respect to the Higgs field, and evaluated at its minimum In such a computation based on the effective potential, lower derivatives of the potential also play an important role. First of all, from the first derivative of V eff we obtain the tadpole equation which allows us -as v = 0 -to eliminate 6 the quantity A (1) Next, the second derivative of V eff defines the effective-potential, or curvature, mass of the Higgs boson, which we denote This equation has a very important consequence, unique to models with classical scale invariance: it relates on the left-hand side the Higgs curvature mass, close to its known pole mass of 125 GeV, to the tree-level masses of all particles in the theory through the quantity B (1) on the right-hand side. If the model we consider contains no extra fermions and gauge bosons but only additional scalars, we obtain [24,50,76,78] using eq. (2.5) having taken into account the N c = 3 colour factor for the top quark and with m t , m W , and m Z the tree-level masses of the top quark and of the W and Z bosons respectively. Once one considers specific scenarios for the particle contents of BSM CSI theories, it becomes possible to obtain constraints on the masses of the new states -as studied for example in refs. [32,50,79]. We will return to this type of relations among masses in the following, when we consider particular BSM models.
Returning now to the one-loop potential (cf. eq. (2.4)) we find that the only two parameters on which it depends -A (1) and B (1) -are fixed by the tadpole equation and the knowledge of the Higgs mass. This implies that if we take further derivatives of the JHEP03(2021)297 one-loop potential, the computed quantities are also fixed in terms of the Higgs curvature mass. Indeed, the effective Higgs trilinear coupling is at one-loop order 12) and using equations (2.7) and (2.9), we can eliminate A (1) and B (1) , and rewrite λ hhh in terms of the tree-level result for λ hhh in the Standard Model as [50] Similarly, we can compute a one-loop (effective) quartic coupling λ hhhh (2.14) Note that for these two last equations, we made use of the known tree-level results for the Standard Model The predictions in equations (2.13) and (2.14) deviate significantly from the usual SM values -for the case of λ hhh this effect is sufficiently large to be observed already at the HL-LHC. A further crucial aspect of these results is that, in CSI models, the Higgs trilinear and quartic coupling are, at one loop, completely independent of the particle content of the theory. However, as we will see now, this universality is lost when two-loop corrections are included in the effective potential.

Two-loop corrections to the Higgs trilinear coupling
At two loops, corrections to the effective potential are obtained by computing one-particle irreducible vacuum bubble diagrams. In addition to logarithmic and non-logarithmic terms -as are present at one loop -such two-loop diagrams also contain squared logarithms, meaning that the effective potential takes the form where the coefficient C appears only at two-loop order. Expanding this expression in terms of one-and two-loop quantities, it becomes

JHEP03(2021)297
While A (1) and B (1) are exactly the quantities given in eq. (2.5), A (2) , B (2) , and C (2) must be derived from the two-loop diagrams contributing to the effective potential -and will depend on the model one is considering. Before investigating specific models, it is enlighting to repeat the procedure of the previous section, and derive general expressions for the derivatives of V eff , using now the two-loop form in eq. (2.16). First, the tadpole equation reads from which we can isolate A as Next, the Higgs curvature mass is, after eliminating A Equation (2.20) is also important as it constitutes the two-loop version of the relation between masses in CSI models -which we first encountered at one loop in the previous section. Once we consider particular BSM models in sections 3 and 4, we will see that this relation places stringent constraints on the allowed BSM mass and parameter ranges. Now, if we compute the effective Higgs trilinear coupling, we obtain at two loops once A and B are replaced using eqs. (2.19) and (2.21). We see here that because of the new term C, appearing from two loops in V eff , non-universal, model-dependent corrections arise in the expression of λ hhh . In practice we will be computing C from the two-loop diagrams contributing to V eff and use this result to obtain λ hhh at two loops. However, we note that we can formally use the last equation to determine C in terms of [M 2 h ] V eff and λ hhh , as We can finally use this expression of C to compute λ hhhh at two loops in terms of [M 2 h ] V eff and λ hhh only Figure 1. Dominant two-loop diagrams contributing to the SM effective potential.
This equation implies that once λ hhh is measured, the Higgs quartic coupling can be predicted in a universal way for all CSI theories, up to two loops. This being said, once three-loop corrections are included, we can expect this universality to be lost because of log 3 terms appearing in V eff -similarly to what happens with λ hhh once two-loop effects are taken into account.

Results for the classically scale-invariant Standard Model
We end this section by calculating the leading two-loop contributions to V eff and λ hhh for the classically scale-invariant version of the SM, which will also provide an occasion to present the renormalisation scheme conversion that we perform to obtain expressions in terms of physical quantities. We must emphasise that the CSI-SM is already known not to be a valid theory of Nature as it cannot reproduce the Higgs mass correctly [23], however, its simplicity makes it the ideal setting to explain our calculational setup. The tree-level scalar potential of the CSI-SM does not contain a mass term, and reads simply where the Higgs doublet Φ is defined as . (2.26) G and G ± are respectively the neutral and charged would-be Goldstone bosons, while h is the 125-GeV Higgs boson. At tree level, the field-dependent masses of the Higgs and Goldstone bosons read while the top-quark mass is (like for the usual SM) The need for a flat direction in the tree-level potential imposes λ = 0, which is also the result of solving the tree-level tadpole equation.

JHEP03(2021)297
The dominant contributions to the CSI-SM two-loop effective potential arise from the diagrams shown in figure 1 and involving the SU(3) c gauge coupling g 3 and the top Yukawa coupling y t -the two EW gauge couplings g 2 and g Y , and the light fermion Yukawa couplings are negligible in comparison. These contributions read (in the MS scheme) where J and I are respectively one-and two-loop functions, for which definitions 7 and useful limits are given in appendix A.
In the previous section, we found that to obtain the two-loop corrections to λ hhh , we need to compute C, i.e. the coefficient of the log 2 term in V eff . For the CSI-SM, after expanding the expression in equation (2.29), we find Using then equation (2.22), we obtain for the effective Higgs trilinear coupling in the CSI-SM, to leading two-loop order, At this point, we should also discuss the choice of renormalisation scheme for the parameters appearing in the results presented in this work. As the effective potential, from which we derive λ hhh , is calculated in the MS scheme, the expression in equation (2.31) should be understood as being written in terms of MS parameters. In the following, we will however choose to convert MS expressions into on-shell (OS) ones -i.e. written in terms of physical quantities and including also finite corrections from wave-function renormalisation (WFR). An OS version of the Higgs trilinear coupling can be obtained aŝ are respectively the on-shell and MS scheme Higgs WFR constants, Π (1) hh is the finite part of the one-loop Higgs boson self-energy, and with λ hhh expressed in terms of physical quantities. Note that it suffices to compute the WFR constants to one-loop order, because λ hhh is generated at one loop.
In the two-loop result for λ hhh in eq. (2.31), m t and g 3 only enter at two loops, and therefore they require no scheme conversion here -this would only produce corrections at three-loop order. However, it is necessary to include a one-loop scheme conversion for

JHEP03(2021)297
the Higgs MS VEV v and curvature mass [M 2 h ] V eff , both appearing from one loop. First, the MS VEV can be related to the Fermi constant at (leading) one loop as For notational convenience, we also define an OS Higgs VEV as v 2 Next, we recall that the Higgs pole (i.e. physical) and curvature masses are related by the following formula In the CSI-SM, the Higgs mass is generated at loop level, and thus this relation simplifies into Combining all these intermediate results together, we obtain finallŷ It is important to note that, because M 2 h is generated at loop level, the second term, which is proportional to κM 2 h M 2 t /v 3 OS , is formally of two-loop order -only the first term in this expression is of one-loop order.
Taking from the PDG [80], the following values for the physical inputs We find numerically 38) or in other words, the two-loop corrections to λ hhh yield a positive shift of the order of 2%.

Model definitions
The first class of BSM CSI theories that we consider are extensions of the SM with N additional real singlet scalar states. We will furthermore assume that these models are endowed with a global O(N ) symmetry, under which only the singlet scalars are charged while the SM particles do not carry any charge. Although this assumption greatly simplifies

JHEP03(2021)297
the expressions we obtain -by reducing the number of different masses and coupling constants -it is important to note that it does not significantly affect our discussion about corrections to the Higgs trilinear coupling. Such O(N )-symmetric CSI models have been discussed in many references -see refs. [81][82][83][84][85] -from various points of view.
In particular, the exact O(N ) symmetry provides a symmetry to stabilise dark matter candidates in these models. It was found in ref. [82] that there is an upper bound on the value of N if one considers direct dark matter search results. However, in this paper, we choose to study O(N ) models exclusively from the point of view of Higgs physics and corrections to λ hhh , and we do not consider potential dark matter aspects of these models.
In other words, we consider here O(N ) as a phenomenological symmetry -which greatly simplifies the study of Higgs physics in N -scalar theories -rather than a symmetry to stabilise dark matter. We can write the scalar potential of O(N )-symmetric models with CSI as where Φ is the SM-like Higgs doublet -same as that shown in eq. (2.26) -and S = (S 1 , S 2 , · · · S N ) is a real scalar, singlet under the SM gauge groups and belonging to an Ndimensional representation of the global O(N ) symmetry group -another way to see this is to think of N real scalar singlets. Only the neutral component of the Higgs doublet acquires a non-zero VEV -the singlets do not as this would break the global O(N ) symmetry.
The scalar masses at tree-level are and the tree-level tadpole condition simply gives λ = 0. Finally, it follows immediately from this equation that we can replace λ ΦS by the tree-level BSM scalar mass m S , as

Dominant two-loop corrections to the Higgs trilinear coupling
We can obtain for the leading BSM contributions to V (2) by considering the new diagrams involving λ ΦS and λ S , shown in figure 2. Their expressions are The squared-logarithm terms in these two-loop contributions are, in turn Using once again equation (2.22), we can obtain the leading two-loop order MS result for λ hhh Among the terms in the brackets, the first two (those involving m t ) are the SM-like contributions, while the last two are the new BSM effects. Note that, like in the SM case, m S and λ S only enter the expression of λ hhh at two loops, and hence we do not need to specify their renormalisation scheme. If we want to express the one-loop result in terms of M 2 h and v OS we must however include scheme conversions for both [M 2 h ] V eff and v and take into account WFR effects -as discussed already in section 2.5. The only difference we must consider arises from the new momentum-dependent term in the Higgs self-energy in the O(N )-symmetric model, which reads where the one-loop Passarino-Veltman function [88] B 0 is defined in appendix A. Adding this additional term, we find the OS expression forλ hhh at leading two loops to bê Once again, we emphasise that the terms in the first pair of brackets -proportional to κM 2 h /v 3 OS -are formally of two-loop order. In principle the above result forλ hhh involves 3 BSM parameters, namely M S , λ S , and N . However, the classical scale invariance produces a further constraint on the BSM parameters, by relating the masses of the particles in the theory, as we will see in the following section.

Mass relation
When investigating derivatives of the effective potential in general CSI models in section 2.4, we had found that the curvature mass of the Higgs boson is related to the B and C coefficients in V eff . In terms of one-and two-loop quantities, this relation reads As we have computed the effective potential in O(N )-symmetric CSI models, it is a simple matter to derive all the necessary ingredients to apply this relation to this class of models.
In addition to C (2) , which was given in the previous section, we also need the logarithmic terms at one and two loops -B (1) and B (2) . Including both SM-like and BSM terms, these are Inserting these into equation (3.8), we obtain in terms of MS-renormalised parameters (3.10) The first line is simply the one-loop result, while the second and third lines are respectively the SM-like and BSM two-loop corrections to the relation between masses.
In order to convert this result in terms of physical masses for the SM and BSM states, we require the Higgs and BSM scalar self-energies -provided in appendix B.1 -as well as equation (2.33) for the Higgs VEV. We obtain finally If one considers this relation from the point of view of computing M S in terms of λ S , N , and the known SM inputs, the two-loop corrections decrease the extracted value of M S . However, in our numerical investigations, we will instead use this relation to compute λ S as a function of M S and N .

Numerical study
At this point, we have derived all the analytical results necessary to study the possible magnitude of the leading two-loop corrections to the Higgs trilinear coupling in the O(N )symmetric CSI model. Rather than examining the absolute value ofλ hhh , we prefer to investigate by how much it deviates from the prediction in the (non-CSI) Standard Model, at the same leading two-loop order. Therefore, we present in the following results for the BSM deviation δR, defined as whereλ hhh is the leading two-loop result for the Higgs trilinear coupling in the OS scheme and the SM -its analytical expression can be found for instance in ref. [75] (see also refs. [73,74]). Taking the same input values as in eq. (2.37), we obtain numericallŷ λ SM hhh 177.4 GeV. We present in figure 3 our results for δR using only the expressions from section 3.2 -in other words, we do not impose the relation between masses from eq. (3.11) yet. The red dot-dashed line shows the one-loop result, comparing the one-loop CSI result (λ hhh ) (1) = 5M 2 h /v OS to the (leading) one-loop expression in the SM -see for instance ref. [40]. We should point out here that this result differs from the ∼ 67% mentioned in the introduction, and shown with the black dotted curve: we have (δR) (1) 82%. This is because the value of ∼ 67% is obtained by comparing the one-loop (λ CSI hhh ) (1) to the tree-level

JHEP03(2021)297
SM result for (λ SM hhh ) (0) , but once we include the one-loop contribution from the top quark in the SM expression we obtain a smaller value for (λ SM hhh ) (1) < (λ SM hhh ) (0) , and hence a larger deviation. We must emphasise that with our definition of (δR) (1) we are comparing oneloop effective-potential results, which do not include any dependence on external momenta. It should however be noted that, as studied in ref. [40], the momentum dependence in the SM result can be quite significant, and could modify significantly the value of the BSM deviation. In ref. [85], momentum-dependent corrections were also included in the study of the Higgs trilinear coupling in O(N )-symmetric CSI models, and were found to be nonnegligible. Nevertheless, we still choose to define (δR) (1) in terms of effective Higgs trilinear couplings, as this allows simpler comparisons with our new two-loop results, calculated with the effective potential. The result we find for (δR) (1) = (λ CSI hhh ) (1) /(λ SM hhh ) (1) − 1 remains totally independent of the particle content of the CSI model.
Next, the blue and green curves in figure 3 are the values of δR at two loops, for different choices of N and λ S : blue and green curves correspond respectively to N = 1 and N = 4; solid curves are made for λ S = 0 while the dashed ones are for λ S = 2. For the lowest values of M S -i.e. where the BSM effects are the smallest -the two-loop deviation (δR) (2) is only slightly larger than at one loop, however, as one increases M S , the BSM contributions grow rapidly. This growth is stronger for larger values of N and λ S , and, crucially, the universality of the BSM shift found at one loop is entirely lost at two loops. Looking at concrete numerical values of (δR) (2) , for example for M S = 500 GeV, we have on the one hand for N = 1 and λ S = 0 a deviation of ∼ 94%, while on the other hand for N = 4 and λ S = 2 we find (δR) (2) ∼ 168%. Anticipating a little on the discussion in the next pages, we note that all parameter points in this figure 3 fulfill the criterion of tree-level perturbative unitarity [52,89].
Another interesting comparison we can perform is to investigate the difference in the BSM deviation inλ hhh between CSI and non-CSI versions of O(N )-symmetric models. For this purpose, we derive in appendix B.2 the dominant two-loop corrections toλ hhh in a non-CSI variant of the O(N )-symmetric model, and for generic values of N -this calculation follows the procedure developed in refs. [74,75] (note also that the case of N = 1 was already treated in ref. [75], under the name of "Higgs singlet model"). We present in figure 4 a comparison of our results for the BSM deviation δR computed with respect to the SM result in both CSI and non-CSI O(N )-symmetric models, for N = 1 (left side) and for N = 4 (right side). Note that, as there is no mass term in the CSI O(N )-symmetric model, we must setμ S = 0 in the result for the non-CSI version of the model in order to be able to compare both variants consistently. The solid curves in figure 4 are obtained for the versions of the models with CSI, while the results for the non-CSI variants are given by the dashed curves. Additionally, in both plots, the red curves show one-loop values, while the blue and black (green and brown) curves are the two-loop results for N = 1 (N = 4) for λ S = 0 and λ S = 2.
One clearly notices that CSI and non-CSI variants of this BSM model behave quite differently. First, for low BSM masses, the BSM deviation in the non-CSI case is minute, while in the CSI case the deviation is already large (80% or so) -this is because of the difference between the normal SM and CSI SM. Then, when M S increases, the non-CSI corrections grow much faster than their CSI counterparts, as in the non-CSI models the one-loop corrections are proportional to M 4 S -cf. equation (B.9). Until this point, we have not taken into account the mass relation -as derived in section 3.3 -nor the theoretical constraint from unitarity. For the latter, we choose to use as our criterion tree-level 8 perturbative unitarity [89], and we can for the O(N )-symmetric model employ the results given in ref. [52], in particular equation (22) therein. We show in figure 5 the regions of parameter space -in the M S -λ S plane -allowed (light green) and excluded (light red) under this criterion of tree-level perturbative unitarity, for both N = 1 (left) and N = 4 (right). Moreover, we superimpose on top of this plane the curves obtained when computing λ S as a function of M S (and N ) using the mass relation we derived in eq. (3.11) -the blue line is for N = 1 while the red one is for N = 4. We recall also that negative values of λ S are forbidden, otherwise the potential would not be bounded from below. Although a wide range of BSM masses seems to be allowed by unitarity, once one takes into account the relation between masses in the CSI model -or in other words once one imposes that M h = 125 GeV -only a narrow interval remains. Concretely, we find for N = 1 that 509 GeV M S 524 GeV, and for N = 4 that 365 GeV M S 375 GeV. 8 We want to constrain the values of λS and MS, both of which appearing inλ hhh only at two-loop order. Therefore, the difference from including higher-order (instead of tree-level) perturbative unitarity constraints should be formally a three-loop effect, and thus within the expected theoretical uncertainty of our result. Another important criterion that we use to verify the validity of the parameter points we consider is the true-vacuum condition [90] -i.e. the requirement that the EW vacuum is the true minimum of the (two-loop) potential and corresponds to a lower value of V eff than the origin of the potential (at ϕ h = v + h = 0). In CSI theories, as all states acquire their masses entirely from their coupling to the Higgs boson, all masses vanish at the origin v + h = 0, and hence V eff (ϕ h = 0) = 0 and we only need to make sure that V eff (ϕ h = v) is negative. The effective potential depends not only on the field ϕ h , but also on the renormalisation scale Q and therefore we examine the value of V eff (ϕ h = v) for Q = v and Q = M S -the two choices that we believe to be the most natural. We find for N = 1 (where M S ∼ 515 GeV) that the true-vacuum condition is fulfilled for Q = M S but not for Q = v, while for N = 4 (where M S ∼ 370 GeV) the condition is verified for both values of Q. As it is unclear which value of Q would be optimal (both values being around the EW scale anyway), the case of N = 1 is not conclusively excluded by this criterion.
Finally, we present in figure 6 the BSM deviations δR we obtain once we take into account the mass relation and the constraint from perturbative unitarity. As for previous figures, the left and right panes are for N = 1 and N = 4 respectively, and the red dot- We emphasise that there is now only a single result for δR at two loops because once M S and N are fixed, the only remaining parameterλ S -is fixed via equation (3.11). The ranges of M S for these plots correspond to λ S varying between 0 and 4π, but the light-red shaded regions are excluded by the criterion of tree-level perturbative unitarity, as discussed above.
Taking into consideration all these theoretical constraints, we find that the two-loop BSM deviation toλ hhh must be within the intervals from ∼ 99% to ∼ 114% and from ∼ 93% to ∼ 109% respectively for N = 1 and N = 4. This means that while the universality found at one loop for (δR) (1) is lost at two loops, only a limited range of values of (δR) (2) is possible. Moreover, our computation finds the Higgs trilinear coupling in CSI O(N )symmetric model to be about twice as large as in the SM -(δR) (2) ∼ 100% ± 10%and is larger than what is obtained at one loop, which makes the BSM deviation easier to access in experimental searches.
Before ending our discussion of the CSI O(N )-symmetric models, it is also important to comment on the experimental status of the parameter points we have considered. We have verified using HiggsBounds [91][92][93][94][95]

Model definition
We now turn to the case of a CSI version of the Two-Higgs-Doublet Model, as devised in ref. [76] (see also ref. [103]). We will be assuming here that CP is conserved in the Higgs sector, and thus we can write the tree-level potential of the model in terms of two SU(2) L doublets Φ 1 and Φ 2 of hypercharge 1/2 as where all the quartic couplings λ i are real. To prevent the occurrence of tree-level flavourchanging neutral currents, we impose a discrete Z 2 symmetry [104,105] under which the scalar doublets transform as while the fermions remain unaffected. We expand the two doublets as where v i (i = 1, 2) are the vacuum expectation values of the neutral components of the doublets. These can be taken to be real -owing to the assumption of CP conservationand furthermore verify the relation v 2 1 + v 2 2 = v 2 (246 GeV) 2 . From the ratio of v 1 and v 2 , we can also define the usual quantity tan β ≡ v 2 /v 1 .
Requiring both VEVs to be non-vanishing, we find for the tree-level tadpole equations with the shorthand notation λ 345 ≡ λ 3 + λ 4 + λ 5 . These equations can be solved as

JHEP03(2021)297
Next, mass eigenstates are obtained by rotating the gauge eigenstate fields with the angle β having defined Among these fields, h, H are CP-even Higgs boson, A is a CP-odd Higgs boson, and H ± is the charged Higgs boson. As in the SM, G and G ± are respectively the neutral and charged would-be Goldstone bosons.
In this new basis, the tree-level scalar mass eigenvalues can be found to be Interestingly, and to the difference of the usual (non-CSI) 2HDM, the CP-even mass matrix is already diagonal after the rotation of angle β. In other words, no additional rotation of angle α − β is necessary, and the CSI-2HDM is naturally aligned (at tree level). Moreover, finding m h = 0 also tells us that h is indeed the flat direction along which we want to work in what follows. We also follow the common choice of trading the free Lagrangian quartic couplings λ 3 , λ 4 , λ 5 for the tree-level masses of the scalars. From eq. (4.9), one has We should note also that in the following, we neglect loop-induced deviations from the alignment found at tree level. Such effects have been studied for instance in ref. [32] in the CSI 2HDM (see also ref. [106] for the usual 2HDM) and were found to be relatively small for most of the parameter space, so that we consider them subleading before the two-loop corrections from the BSM scalar states. This enables us to write the effective potential under the form of equation (2.16) and thus to follow the calculational procedure described in section 2. As in the usual 2HDM, the two scalar doublets in the CSI-2HDM can couple to quarks and leptons. Under the Z 2 symmetry, four types of coupling assignments are possible, corresponding to the four types of 2HDMs [107][108][109]. We note however that, for the major part of the discussion in this section, we do not need to specify a type as we only consider effects from the top quark -only when verifying experimental constraints do we assume interactions of type I in order to have less severe limits from flavour physics than for instance for types II and Y [110].

Dominant two-loop corrections to the Higgs trilinear coupling
At one loop, the result for λ hhh is model-independent and was discussed in section 2. We can now consider the dominant two-loop contributions to the effective potential, and in turn to the Higgs trilinear coupling. We choose in this section to present results for the case of degenerate BSM scalar masses M H = M A = M H ± ≡ M Φ , because the expressions for general values of the BSM scalar masses are quite long and cumbersome, and moreover because this will not affect the physics we discuss here. We provide the complete expressions for V (2) , λ hhh , andλ hhh in appendix C.1. We should also note that as the BSM scalar masses will only appear in the Higgs trilinear coupling starting from the two-loop order, it does not make a difference whether we use pole or tree-level masses for the BSM scalars -the difference in λ hhh from using one or the other is formally of three-loop order.
The BSM diagrams contributing to leading order to V (2) are shown in figure 7. In the limit of degenerate masses, their expressions read From these equations we can extract the BSM logarithm-squared terms, and including also JHEP03(2021)297 the SM-like contributions, we find in total for C (2) (4.12) Applying equation (2.22), we obtain the dominant two-loop corrections to λ hhh , involving g 3 , m t , m Φ , and tan β as Lastly, we can -like in the previous sections -express the one-loop result 5[M 2 h ] V eff /v in terms of the Higgs pole mass and the physical Higgs VEV, and take into account finite WFR effects, to obtain an OS result. The only 2HDM-specific corrections come from the momentum-dependent BSM scalar contributions to the Higgs self-energy. These can be found to be where n H,A = 1 and n H ± = 2. From this, we have dΠ (1) hh dp 2 Finally, we obtain in terms of OS-renormalised parameterŝ (4.16)

Mass relation
Next, we also consider the mass relation for the case of the CSI-2HDM, which will serve to relate the BSM scalar masses M H , M A , M H ± and tan β. As in the previous section, we present here only expressions in the limit of degenerate scalar masses, and leave the complete results for appendix C.2. To derive corrections to λ hhh , we have already used the expression of C (2) given in eq. (4.12). We must now also derive results for B at one-and two-loop orders. First, at one loop it is straightforward to calculate

JHEP03(2021)297
At two loops, we have FF S , (4.18) where the first two terms are SM-like contributions and the latter three -B (2) SSS , B (2) FF S -are the BSM ones. These can be found to be in the limit of degenerate BSM masses Applying equation (3.8), the mass relation at two loops is, in terms of MS-renormalised quantities, (4.20) The first line in this equation is the one-loop result, while the additional lines are the twoloop corrections. Among these, the first (square) brackets are the SM-like contributions, and the remains terms (in the curly brackets) are the BSM terms, proper to the CSI-2HDM.
To convert equation (4.20) from the MS to the OS scheme, we require as previously the correction to the Higgs VEVc.f equation (2.33) and the top-quark and Higgs-boson one-loop self-energies, as well as the self-energies for the BSM scalars H, A, H ± (also at one loop). All necessary expressions are provided in appendix C.2. With these, we are able to obtain the following result, in terms only of physical parameters where the function f B is defined in equation (A.6), and Re[x] denotes the real part of a given quantity x.

Numerical study
We have now derived all the results we require to study numerically the leading two-loop corrections to the Higgs trilinear coupling in the CSI-2HDM. We follow here the same outline as we did for the O(N )-symmetric model, and we provide results for the same type of BSM deviation δR -as defined in eq. consider the three BSM scalars H, A, H ± to be degenerate in mass -this ensures that the ρ parameter remains close to one, and furthermore this conveniently reduces the number of BSM mass scales to one.
We present first in figure 8 the deviation δR at one loop (red line) and two loops (grey, teal, blue, and green curves) for different values of tan β. We employ here the expression in equation (4.16), without for now taking into account the relation between masses. As we had seen for the N -scalar model, while at one loop the BSM deviation inλ hhh has a universal value of δR ∼ 82%, the two-loop corrections introduce new dependences on both M Φ and tan β, thereby spoiling the universality of δR. We consider in this plot different values of tan β, ranging from tan β = 1 to tan β = 1.4, which we find to be the largest value of tan β for which tree-level perturbative unitarity [89,111] is maintained all the way to M Φ = 500 GeV -see figure 10 and the discussion in the following. Considering for instance M Φ = 500 GeV, we find (δR) (2) = 134% for tan β = 1, and (δR) (2) = 162% for tan β = 1.4, in stark contrast to the one-loop result.
At this point, we should comment briefly on our interpretation of the quantity tan β. While tan β is originally defined as the ratio of the VEVs of the two scalar doublets in the model, it appears in two-loop corrections to the Higgs trilinear coupling via the dependence of four-point scalar interactions at tree level and of the coupling between BSM scalars and the top quark. It should be emphasised however that because tan β appears only from the two-loop order, any scheme difference in how its value could be extracted from some experimental data (in the event a 2HDM Higgs sector is discovered) only has an impact from three-loop order and beyond -i.e. below the level of accuracy to which we are working. Next, it is also interesting to compare predictions for the Higgs trilinear coupling in variants of the 2HDM with or without CSI. For this purpose, we present in figure 9 the BSM deviations inλ hhh computed with respect to the usual-SM prediction in both variants of the 2HDM. For the usual 2HDM, we use the analytic expression forλ hhh derived in refs. [74,75] -see specifically equation (5.15) in ref. [75] -in which we set the mass term M (representing the scale of the soft breaking of the Z 2 symmetry acting on the Higgs doublets) to zero to allow the comparison with the CSI scenario.
We can observe that the Higgs trilinear coupling behaves very differently in CSI and non-CSI versions of the 2HDM. The main discrepancy occurs at one-loop level: indeed, while in the CSI-2HDM (δR) (1) is constant, in the usual (non-CSI) 2HDM, one-loop corrections involving the BSM scalars are minute in the low-mass range, but grow rapidly (as M 4 Φ ) and reach 240% for M Φ = 500 GeV. At two loops, the radiative corrections toλ hhh behave quite similarly in the CSI and non-CSI scenarios -the main difference arising from effects due to the scheme conversion of the one-loop BSM corrections in the non-CSI 2HDM case (cf. ref. [75]).
At this point, we should discuss the theoretical constraints on the BSM parameters -M Φ and tan β -coming from unitarity and from the mass relation in eq. (4.21). For the former, we choose to take tree-level perturbative unitarity [89] as our criterion, and we employ results from ref. [111] (see also ref. [112]). by the dashed red (at one loop) and solid blue (at two loops) curves. On the one hand, at one loop tan β does not appear in the mass relation and hence at a constant value is found for M Φ . On the other hand, at two loops, the result for M Φ obtained from eq. (4.21) becomes a function of tan β. Along the two-loop (blue) curve, we find that the maximal possible value of tan β is 1.855, which corresponds to a maximal value of the degenerate mass of the BSM scalars M Φ ≤ 382.2 GeV. As for the O(N )-symmetric models, we have also verified that the true-vacuum condition is respected in the CSI-2HDM: we have found that the value of the two-loop potential is, at the EW minimum, indeed lower than the value of the potential at the origin, and this is true both for Q = v and Q = M Φ ∼ 378 GeV.
Having considered the constraints on M Φ and tan β, we can finally investigate the theoretically-allowed values of the Higgs trilinear coupling. figure 11 therefore presents our results for the deviation δR at one loop (red dot-dashed line) and at two loops (blue solid curve), once tan β is extracted as a function of M Φ from the mass relation in eq. (4.21). As can be seen in figure 10, the possible values of M Φ are all within the narrow interval from 373.9 GeV to 382.1 GeV. We find that δR varies between 89.9% at the lowest, and 112.5% at most. In other words, as for the O(N )-symmetric model earlier, we obtain a larger deviation ofλ hhh from its deviation at two loops, and the total BSM deviation is (approximately) in the range 100 ± 10%.

JHEP03(2021)297
We have verified that the corresponding parameter points are not excluded by experimental searches for BSM scalars. Assuming Yukawa interactions of type I (to benefit from less severe constraints from flavour physics than for types II or Y, see e.g. ref. [110]), we create a SPheno-based spectrum generator for the CSI-2HDM using SARAH, which we in turn used to provide inputs for HiggsBounds to verify the experimental status of the parameter choices we took. We ensured that the Lagrangian inputs provided to SPheno yielded the corrected spectrum of scalar masses -a 125-GeV SM-like Higgs boson and BSM scalars at 374 GeV, with an allowed error of 1 GeV -using the Mathematica package SSP [113].

Discussion
The computations presented in this paper allow studying the Higgs trilinear coupling in models with classical scale invariance to the same level of accuracy as for non-CSI extensions of the SM -as done for instance in refs. [73][74][75] (and appendix B.2 of this work). Classical scale invariance is an attractive concept for model building, and may relate to the solution to the hierarchy problem. It is therefore of paramount importance to find new ways to investigate CSI models, in complement of collider and dark matter searches -see e.g. refs. [32,79,[81][82][83][84][85][86][87].
Due to the simple form of the potential at one loop, containing only two free quantities -recall equation (2.4) -λ hhh is found to be predicted universally at the one-loop level, in all CSI models. This universality of (λ hhh ) (1) is lifted by the inclusion of two-loop effects, which thus allows distinguishing different scenarios of CSI models. We found that the new corrections at two loops are quite significant, leading to a further 5-30% deviation ofλ hhh from the (non-CSI) SM prediction, compared to the one-loop CSI result. We should emphasise that this is by no means a problem from the point of view of the validity of the perturbative loop expansion: indeed the two-loop corrections toλ hhh involve new parameters that are not present at one loop -namely M S , λ S , N for the O(N )-symmetric models and M H , M A , M H ± , tan β in the CSI-2HDM. Furthermore, as there are numerous arguments in favour of extended scalar sectors, it is interesting to note that the two-loop contributions to the Higgs trilinear coupling receives from scalars are always positive, which makes it deviate further from the SM value so that it is more easily accessible in collider searches than what is expected from the one-loop result.
We have also paid particular attention to the relation among masses arising in theories with CSI, which is often over-looked in the literature. More than being simply bounded from above, the masses of BSM states in CSI models are actually strongly constrained by the requirement of reproducing the correct 125-GeV mass of the Higgs boson -which is entirely generated at loop level. In the scenarios considered in the main text of this paper, we imposed additional constraints -a global O(N ) symmetry for the N -scalar model, and mass degeneracy of the BSM scalars in the CSI-2HDM -which resulted in severe bounds on the allowed BSM masses. 9 In turn, this allows us to find predictions forλ hhh JHEP03(2021)297 within a narrow range: in both types of BSM models, we obtain 90% (δR) (2) 115%. Looking back at figures 4 and 9, it is worth noting that the values obtained for (δR) (2) for the allowed BSM mass ranges are very close to those found for the same masses in the non-CSI counterparts of the scenarios. Taking the example of the 2HDM, while in the CSI version of the model we have 90% (δR) (2) 110%, in the non-CSI variant and taking the non-decoupling limit (i.e.M = 0) the value of (δR) (2) computed at M Φ 374 GeV varies from ∼ 90% (for tan β = 1) to 105% (for tan β = 1.9, close to the upper limit from perturbative unitarity). This implies that the magnitude of a potential BSM deviation in the Higgs trilinear coupling may not provide a clear indication of whether the associated new physics exhibits classical scale invariance -instead, a better criterion might come from whether masses and couplings fulfill the mass relations imposed by CSI.
Rather than considering only a single type of observable, a powerful strategy to investigate CSI theories, and to discriminate theories with or without CSI, would be to study the correlation between several observables. In particular, synergies between measurements of the Higgs trilinear coupling and direct searches of new particles at colliders may prove crucial to ascertain whether Nature exhibits classical scale invariance or not. Indeed, we should recall that due to the relation among masses in CSI models, the new BSM states cannot be arbitrarily heavy (even in non-degenerate scenarios) and will be within reach of experiments in a foreseeable future (see for instance the discussions for the CSI-2HDM in refs. [32,79]). In the eventuality that a new scalar state is found in direct searches, determining the value of the Higgs trilinear coupling and verifying whether the relation among masses is fulfilled or not can potentially provide clear evidence of the CSI nature of the underlying theory. Furthermore, if one considers a model with charged scalars (e.g. H ± in the CSI-2HDM), sizeable one-loop corrections to the Higgs-to-two-photon (hγγ) coupling can also be expected, as studied in ref. [50]. Such effects could also help confirm the existence (or absence) of CSI in BSM models. In this context, depending on the achieved accuracy of the determination of the hγγ coupling at future colliders -possibly down to percent level, see e.g. refs. [63,114] for prospects at the ILC -it may become important (if not unavoidable) to compute theory predictions for this coupling at two-loop level in order to match the precision of experimental results and allow consistent comparisons between theory and experiment. We leave this endeavour for future work.
A caveat that should be mentioned is that the scenarios considered in this work require relatively large scalar quartic couplings -once again, this comes from the need to generate the correct mass of the Higgs boson, in the absence of any BSM mass term in the Lagrangian. Consequently, under renormalisation-group running, we can expect these couplings to grow rapidly and encounter Landau poles well before the Planck scale -in a similar way as what occurs in the non-decoupling limit of (non-CSI) extensions of the SM. This is consistent with our expectation that the assumption of CSI at the electroweak scale can a priori not be extended up to the Planck scale, and therefore does not solve the hierarchy problem by itself.

JHEP03(2021)297
Finally, we should also point out the importance of computing the Higgs trilinear coupling to high precision when studying the behaviour of a number of BSM phenomena occuring in models with extended scalar sectors -in particular the possibility of a strongfirst order electroweak phase transition and the spectrum of gravitational waves that may be produced during the former. In complementarity with collider searches, the future measurement of primordial gravitational waves (GW) at the LISA and DECIGO spacebased interferometers will provide a new way to probe the value of the Higgs trilinear coupling and the shape of the Higgs potential. The synergy between the two -collider and GW measurements -was studied at one loop for O(N )-symmetric CSI models in ref. [52] (see also refs. [115][116][117] for similar studies in non-CSI theories). Given the large deviation from the SM in the one-loop predictions for λ hhh , one expects in CSI models the EWPT to be of strong first order (cf. ref. [37]). However, two-loop corrections to the Higgs trilinear coupling certainly affect the strength of the EWPT to some extent, although for the time being this type of effects have only been studied for a non-CSI IDM in ref. [73] -the EWPT was then found to be slightly weaker once two-loop corrections toλ hhh where included. In turn, this also modifies the spectrum of gravitational waves produced during the EWPT if it is of strong first-order nature. There is therefore strong motivation to consider two-loop corrections to the Higgs trilinear coupling at finite temperatures in CSI theories.

Summary
In this paper, we have performed the first explicit calculation of leading two-loop corrections to the Higgs trilinear coupling in models with classical scale invariance, using the effectivepotential approximation. For the wide range of CSI models in which the scalon/Higgs direction does not mix with other states -so that we can compute V eff as in eq. (2.16)we have derived general two-loop results relating the coefficients of the effective potential to the corrections to the Higgs trilinear coupling. Importantly, we find that while at one loop the prediction forλ hhh is the same in all CSI theories, this is not any more the case at two loops.
We focused our investigations on two particular types of CSI theories: first an Nscalar model with a global O(N ) symmetry (where we took either N = 1 or N = 4 for numerical applications), and next a CSI-2HDM. For both models, we computed the leading corrections at two loops involving the BSM scalars -and in the CSI-2HDM, the top quark as well. Because these expressions are derived from V (2) , which is computed in the MS scheme, results are obtained at first in the MS scheme. Therefore, we have included the necessary scheme conversion to rewrite the corrections to the Higgs trilinear coupling in terms of physical parameters (pole masses and physical Higgs VEV) and to take into account the effects from finite WFR. Once we turn to the numerical study of these models, we observe that the two-loop corrections clearly lift the degeneracy inλ hhh among different BSM masses and parameters. Furthermore, the results forλ hhh in CSI and non-CSI variants of the same model strongly differ in their theoretical behaviours, mainly JHEP03(2021)297 because of the difference at one loop -constant (mass-independent) value in the CSI case, and corrections proportional to the fourth power of the BSM mass(es) in the non-CSI case.
However, when we also consider the constraints on the allowed ranges of the BSM parameters coming from the criterion of tree-level perturbative unitarity and from the mass relation in CSI theories, the possible values ofλ hhh become quite restricted, namely 90% (δR) (2) 110 − 115%. In other words, while the Higgs trilinear coupling at two loops is not universally predicted in all CSI models (without mixing), its allowed values are still severely limited by the classical scale invariance. Nevertheless, the additional positive deviation of order 5 − 30% of the Higgs trilinear coupling from its SM prediction at two loops can make the BSM discrepancy easier to find at future experiments, and is of the same order of magnitude as the expected accuracy of the determination of λ hhh at future colliders (as discussed in the introduction). Therefore, two-loop calculations will in the future prove necessary for CSI theories in order to consistently compare theoretical predictions with experimental results.
While it may be difficult to distinguish between CSI scenarios and non-CSI scenarios (in the non-decoupling limit) simply from the size of the Higgs trilinear coupling, our work demonstrates the important role of the mass relation that strongly constrains masses and parameters in CSI theories, and may prove to be one of the most powerful tools to probe the CSI nature of a potential BSM discovery. Additionally, as we have mentioned in the previous section, the synergy of the measurement of the spectrum of primordial gravitational waves produced during the EWPT together with that of the Higgs trilinear coupling can also help distinguish scenarios with or without CSI [52].
Finally, in appendix D we have provided generic MS results for the coefficients B (2) and C (2) (see eq. (2.16)) applicable to all CSI models without mixing, and which served to cross-check some of our calculations. We hope these can be of use for the community to study further CSI scenarios.

Acknowledgments
This work is supported by JSPS, Grant-in-Aid for Scientific Research, No. 16H06492, 18F18022, 18F18321 and 20H00160. This work is also partly supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany's Excellence Strategy -EXC 2121 "Quantum Universe" -390833306.

A Loop functions
This appendix summarises notations and definitions of loop functions used throughout this paper.
First of all, the loop factor is defined as (A.1)

JHEP03(2021)297
Next, we denote the regularisation and renormalisation scales respectively µ and Q. These two scales are related as where γ E 0.577 is the Euler constant. For convenience, we also make extensive use of the notation

A.1 One-loop functions
At one-loop order, we employ the following two integrals Their finite (and -independent) part yield the usual Passarino-Veltmann functions [88] J Some simple limits of particular interest include Finally, when computing contributions from finite WFR, it is useful to expand the B 0 function for low momentum p 2 with the relation

A.2 Two-loop functions
When considering two-loop corrections to the effective potential, we also come to encounter the sunrise integral I (see e.g. ref. [69]) defined as .
Its finite part is with J the one-loop integral from eq. (A.4). Expressions for I, equivalent up to relations among dilogarithms can be found in many references (see for instance refs. [69,70,118,119]). Ref. [119] provides the result in a convenient way, which we reproduce here (A.13) Li 2 (z) ≡ − z 0 log(1 − t)/tdt is the dilogarithm function, and Note that the definition of the function Φ in eq. (A.13) is only valid for x, y < z, and the other branches of the function can be found with the relations Furthermore, a number of useful limits are included in refs. [118,120], among which Next, the one-loop self-energy of the top quark is the same as in the SM (the BSM singlet scalars do not couple to it) and it reads For the conversion from the Higgs curvature to pole mass and the Higgs WFR, we only require the momentum-dependent contributions to the one-loop Higgs self-energymomentum-independent terms cancel out when taking the difference Π hh (p 2 = m 2 h ) − Π hh (p 2 = 0) and moreover their derivative with respect to p 2 is also zero of course. The momentum-dependent corrections are where the first term is the SM-like correction from the top quark and the second term is the BSM contribution from the scalars S i . Finally, we find for the one-loop self-energy of the BSM scalars

B.2 Without classical scale invariance
In order to compare predictions for the Higgs trilinear coupling in O(N )-symmetric models with and without CSI, we derive in this appendix the leading two-loop contributions tô λ hhh in the non-CSI version of these models. Note that the case for N = 1 was already considered in ref. [75] -where it was referred to as the "Higgs-Singlet Model" (HSM). If we do not impose classical scale invariance, the tree-level scalar potential in a theory invariant under a global O(N ) symmetry reads

JHEP03(2021)297
The dominant two-loop contributions to the effective potential in this model are the same as shown in figure 2, and the expressions are (like for the CSI version of the model): Following the calculation method explained in detail in ref. [75] and using the differential operator D 3 defined therein, we find for the dominant two-loop BSM corrections to the Higgs trilinear coupling λ hhh To convert this result to the on-shell scheme, we need the one-loop self-energy of the scalars S i , given in eq. (B.3), as well as the finite counter-term for µ 2 S (to ensure proper decoupling of the OS result when takingμ 2 S → ∞ as discussed at length in ref. [75]). The latter can be obtained straightforwardly by adapting the result in eq. (5.29) of ref. [75] (originally for the HSM) to our case. We find δ OS µ 2 S = 1 16π 2 (N + 2)λ Sμ . (C.9)

C.2.1 Expression in the MS scheme
Finally, we also consider the mass relation for the CSI-2HDM with general BSM scalar masses, and therefore we start by deriving expressions for the quantities C (2) , B (1) , B (2) . First of all, we have that SSS , and C (2) FF S given respectively in equations (C.4), (C.5), and (C.6). Next, B (1) is simply The expression for B (2) is significantly longer, and can be written as

JHEP03(2021)297
potential in general renormalisable theories (note that we employ the MS expressions). The expressions presented here can in principle be applied in order to study corrections to the Higgs trilinear coupling as well as the relation between masses in any CSI theory where the "scalon" direction does not mix with other states.
The couplings are defined among real scalars φ i , Dirac fermions ψ I , and gauge bosons A a µ (and also ghosts ω a ) as follows Note that we work in the Landau gauge, and hence we do not need to care about couplings between scalars and ghosts.
Our notations for the different contributions to the effective potential correspond to those in ref. [70], and we refer the reader to equations (4.2) to (4.21) in this reference for all definitions. We simply list in the following the results we obtain (in the MS scheme) for the coefficients B (2) and C (2) for each type of contribution. Note that we do not give results for A (2) as we do not require at any point in our analysis, and also because they can be obtained straightforwardly by the replacement Q → v in the expressions of ref. [70].