The linear–non-linear frontier for the Goldstone Higgs

The minimal SO(5) / SO(4) σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document}-model is used as a template for the ultraviolet completion of scenarios in which the Higgs particle is a low-energy remnant of some high-energy dynamics, enjoying a (pseudo) Nambu–Goldstone-boson ancestry. Varying the σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} mass allows one to sweep from the perturbative regime to the customary non-linear implementations. The low-energy benchmark effective non-linear Lagrangian for bosons and fermions is obtained, determining as well the operator coefficients including linear corrections. At first order in the latter, three effective bosonic operators emerge which are independent of the explicit soft breaking assumed. The Higgs couplings to vector bosons and fermions turn out to be quite universal: the linear corrections are proportional to the explicit symmetry-breaking parameters. Furthermore, we define an effective Yukawa operator which allows a simple parametrization and comparison of different heavy-fermion ultraviolet completions. In addition, one particular fermionic completion is explored in detail, obtaining the corresponding leading low-energy fermionic operators.


Introduction
The Higgs particle seems to be unnaturally light if there is new particle physics at higher scales to which the Higgs may couple. Barring a Copernican perspective on nature and the conclusion that our generation has completed the discovery of particle physics of the visible world, this puzzle -known as the "electroweak hierarchy problem" -constitutes a pressing question.
The persistent absence of evidence for new resonances in the vicinity of the electroweak scale calls for an in-depth exploration of beyond the Standard Model (BSM) theories which may separate and isolate the Higgs mass from the putative scale of exotic BSM resonances. Pseudo-Nambu-Goldstone bosons (PNGB) are naturally lighter than and decoupled from the rest of the spectrum of their mother theory. This suggested decades ago that the Higgs particle could be identified with the PNGB of some BSM high-energy theory [1][2][3].
In the initial proposal [1] a global SU (5) symmetry was considered for the high-energy strong dynamics. Recent attempts tend to start instead from a global SO(5) symmetry [4,5] at a high scale , spontaneously broken to SO (4) and producing at this stage an ancestor of the Higgs particle in the form of one of the resulting massless Goldstone bosons, with characteristic scale f and ≤ 4π f [6]. The coset SO(5)/SO (4) represents the minimal possibility to interpret the Higgs as a pseudo-Goldstone boson in the presence of a custodial symmetry. The explicit breaking of the global symmetry needed to generate the electroweak scale v = f and a mass for the Higgs usually stems from soft couplings of the high-energy dynamics to the Standard Model (SM) gauge bosons and fermions.
Most of the literature on composite Higgs models based on (or containing) SO (5) assumes from the start a strong dynamics and uses an effective non-linear formulation of the models [5,[7][8][9][10][11][12][13], often denominated "composite Higgs" scenario. The ratio encodes the degree of non-linearity of a given model and is a measure of the fine-tuning required to accommodate data. A complete renormalizable model was instead constructed in Refs. [14,15] (see also Refs. [9,16] for related results), which in its scalar part is a linear sigma model including a new singlet scalar σ . Furthermore, in Ref. [15] the procedure and first steps to obtain the non-linear effective Lagrangian were developed. Later work has impact on interesting phenomenological consequences [17,18] and other aspects [19,20]. That minimal sigma model allows one to gain intuition on the dependence on the ultraviolet (UV) completion scale: it can be considered either as an ultimate model made out of elementary fields, or as a renormalizable version of a deeper dynamics, much as the linear σ model [21] is to QCD. Upon spontaneous breaking of the electroweak symmetry, the Higgs and the σ field mix, and the resulting scalar sector is that of a system with two "Higgs-like" particles. 1 The implications on low-energy precision data and the expected signals at LHC have also been developed in Refs. [15,17,18], leading to a lower bound m σ 550 GeV and interesting σ decay channels into gauge bosons and tt pairs at LHC.
The linear-non-linear divide will be further explored here by varying the mass of the extra scalar 2 σ : a light σ particle corresponds to a weakly coupled regime, while in the high mass limit the theory should fall back onto a usual effective non-linear construction. The effective low-energy Lagrangian for non-linear realizations of electroweak symmetry breaking has been determined previously, but the number of couplings is very large in the most general case [22][23][24][25]. For the generic SO(5)/SO (4) construction, a very reduced subset of those couplings, constituting a complete basis of bosonic operators, was first established in Ref. [26]. This served to generically parametrize scenarios of electroweak symmetry breaking in which the Higgs particle is a low-energy remnant of some dynamics based on or containing SO(5) as a global symmetry. Here we will focus on the particular case of the minimal SO(5) sigma model, leading to an even more reduced subset of operators -the benchmark low-energy effective Lagrangian -which are expected to be common to the non-linear limit of any construction containing the SO(5)/SO(4) spontaneous breaking. We will con- 1 The heavier scalar is called "global Higgs" in Refs. [17,18]. 2 When using further below polar coordinates the extra scalar will be dubbed ρ as customary; see Sect. 2. sider here both bosonic and fermionic operators, though. The leading linear corrections and the leading dependence on the explicit SO(5)-breaking mechanism will be determined as well in this work.
While the bosonic couplings should be universal, the fermionic part may instead be quite model-dependent. Many different choices of exotic fermions have been explored in the literature, mainly within non-linear realizations of the global symmetry (see e.g. Ref. [8]). Nevertheless, a common characteristic of the so-called "partial compositeness" framework is that SM tree-level Yukawa couplings are forbidden by the global symmetry, while instead vertices coupling one or two heavy exotic fermions to the Higgs field are allowed. Effective Yukawa couplings between the SM and the Higgs field are thus induced at low energies, with a generic Seesawlike pattern for the mass generation of SM fermions, whose masses are then inversely proportional to those of the heavy fermions. We will analyze the problem in two approaches: • A rather model-independent one in which the field content of the SM is augmented exclusively by a singlet scalar within the minimal SO(5) setup mentioned, while the leading phenomenological impact of heavy fermions is encoded in an effective Yukawa coupling of the SM fields that we will define. This effective operator will serve to parametrize and disentangle among different choices of BSM fermion embeddings. • In a second step, a concrete choice for the heavy-fermion representations will be considered [15]. This sector will be integrated out explicitly.
Note that in Ref. [15] we had already integrated out the specific BSM heavy fermion mentioned in this second step, although leaving fully dynamical the scalar sector. That is, the effective Lagrangian made out of SM fields plus the σ particle was established. It was also proposed there to next integrate out the latter, which is a straightforward procedure starting from that result. This task will be completed here in order to compare with the case in which the order of integration of the heavy fields (bosons versus fermions) is inverted. The resulting benchmark couplings will also be compared with those stemming from the procedure indicated in the first bullet above. By furthermore keeping track of the linear and heavyfermion corrections, the analysis will provide candles to identify whether a renormalizable ultraviolet completion exists in nature or alternatively an underlying "composite" mechanism is at work at high energy, analogous to QCD for the chiral dynamics involving pions.
Note that the results may be relevant as well for other scenarios based on global groups larger than SO (5). Furthermore, a Goldstone-boson parenthood for the Higgs is not exclusive of strong interacting dynamical setups, but is also embedded in other constructions such as "little Higgs" models, extra-dimensional scenarios and others; our results will then apply as well to those constructions.
The structure of the paper can easily be inferred from the Table of Contents.

Model-independent analysis
Consider a Lagrangian comprising, in its scalar sector L s , a linear sigma model which exhibits a global SO(5) symmetry broken to SO (4) and includes a new scalar, σ , singlet under the SM gauge group [15] where φ = ( H , H, √ 2σ )/ √ 2 is a 5-plet of SO(5) encompassing the Higgs doublet degrees of freedom H in addition to σ . The terms which break softly 3 the SO(5) symmetry -proportional to α and β -endow the Higgs particle with a PNGB character, remaining naturally light as long as α, β λ. The embedding of the gauge group SU (2) L × U (1) Y inside SO(5) is purely conventional. L s contains as well the scalar interactions with gauge bosons, with the SU (2) L × U (1) Y covariant derivative given by where i L and i R denote, respectively, the generators of the SU (2) L and SU (2) R subgroups of the custodial SO(4) group contained in SO (5). Both h and σ acquire a vacuum expectation value (vev), leaving unbroken an SO(4) subgroup which is rotated with respect to the group SO(4) ≈ (2) encodes the kinetic terms for gauge bosons.
Consider now the fermion sector. A generic feature is that the phenomenological constraints on partial compositeness require additional vector-like fermions, which couple and act as mediators among the SM fields. The exact form of the effective coupling is model-dependent and varies according to how the SM fermions are embedded in SO (5). We will obviate until Sect. 3 the details of the heavy-fermion spectrum and use instead in this section a simplified -effectiveapproach to the dominant fermion-induced effects. The fermionic part of the Lagrangian in Eq.
(2) will be written as the sum of two terms, where L kin f,SM comprises the kinetic terms for only SM fermions. All what is needed here in addition is the fact that, in frameworks akin to "partial compositeness", the global symmetry is explicitly broken by couplings between the SM fermions and heavy exotic fermions, which are the source of: (1) non-zero values for the soft SO(5) breaking parameters α and β at one loop, inducing a potential and mass for the Higgs particle; (2) effective Yukawa couplings for the SM fermions and thus the generation of SM fermion masses.
The schematic effective Yukawa coupling in the presence of the σ particle is presented in Fig. 1. It follows that at low energies it is possible to write an effective Yukawa Lagrangian in terms of only the SM fermions, plus h and σ , which respects electroweak gauge invariance but not SO (5) invariance, where the constant y 0 f is a model-dependent coefficient 4 and we define the effective Yukawa operator for a given fermion f as with H ≡ iσ 2 H * . The ellipses in Eq. (6) refer to other SM fermion operators and possibly extra model-dependent terms coming from the heavy-fermion sector.
In the literature of composite Higgs models the notation MCHM A−B−C is often used to indicate their fermion composition, with A, B, C indicating the SO(5) representation in Table 1 Yukawa operators corresponding to particular embeddings (see e.g. Ref. [8]) of a SM quark doublet q L and right-handed q R fermion (either up-type or down-type right-handed quark) into SO (5). The coefficients y, y refer to distinct possible relative weights of SO(5) invariant operators allowed by the models which the SM doublet q L , up-type right-handed and downtype right-handed fermions are embedded, respectively; else, when only one subindex appears as in MCHM A it is understood to be of the type MCHM A−A−A . Table 1 summarizes the {n, m} parameter values for different models. 5 Equation (7) assumes that a given fermion mass corresponds to a single set of {n, m} values. This is often the case; for instance, the top and bottom Yukawa couplings in the Ref. [27]). Notice that, for the cases with a single Yukawa operator, the global coefficients and suppression scales in Eqs. (6), (7) are constrained by the fermion masses and therefore do not constitute any additional model dependence.
Nevertheless, in some scenarios a given fermion mass results instead from combining several operators of the type in Eq. (7) with different {n, m} values. The procedure derived can easily be extended to encompass it. A model-dependence remains then in the relative size of the y and y weights in Table 1. An example is the MCHM 14−14−10 scenario [8] in which different sets of {n, m} values are involved in generating the top mass, while the bottom mass only requires set {n, m} = {1, 0}. The cases of single and of multiple Yukawa operators contributing to a given mass will be further considered explicitly below. We focus in the following on the top Yukawa coupling unless otherwise explicitly stated, while the conclusions to be obtained are easily generalized to all light fermions.

Polar coordinates
Armed with the tools described, it is quite straightforward to derive the benchmark bosonic Lagrangian as well as the lead-ing couplings involving fermions. To this aim, it is convenient to rewrite the scalar degrees of freedom in polar coordinates, with where (x) denotes the Goldstone matrix corresponding to the longitudinal components of the electroweak gauge bosons. In this notation the scalar Lagrangian in Eq.
(2) reads where denotes the trace and V μ ≡ (D μ U )U † as is customary. The effective top Yukawa operator in Eqs. (6) and (7) is then given by where the right-handed SM fermions have been gathered in ) being a projector onto the up-type (downtype) right-handed SM fermions. The ρ and ϕ fields will develop vevs, where at the minimum of the potential the ϕ field corresponds to The connection between the vevs of the fields in the linear and polar parametrizations is The scalar resonance, which in the linear parametrization is customarily denoted σ , is traded by ρ in the polar parametrization, with m ρ = m σ exactly as expected for a physical observable, while the Higgs resonance h corresponds now to the excitation of the ϕ field; see Eq. 12. Finally, as the pure gauge Lagrangian L g and the weak coupling to fermions are not modified, the coefficient of the W μ mass term in Eq. (10) allows one to identify the electroweak scale v in terms of the Lagrangian parameters: Expansion in 1/λ The scalar quartic coupling λ can be conventionally traded by the ρ mass, given by m 2 ρ 8λ f 2 for negligible α and β, see Ref. [15] and further below; the non-linear model would be recovered in the limit m ρ f , that is, λ → ∞. Varying the ρ mass (that is, λ) allows one to sweep from the regime of perturbative ultraviolet completion to the non-linear one assumed in models in which the Higgs particle is a lowenergy remnant of some strong dynamics. We will explore this limit next.
The exact equation of motion for ρ reads where ≡ ∂ μ ∂ μ . In a 1/λ expansion, the ρ field can be expressed as where the leading terms are given by and subsequent ones can be written as polynomial functions of ρ 1 . Substituting these in Eq. (10) yields the 1/λ n Lagrangian corrections, where the different terms in this equation are given by L 0 coincides with the leading-order Lagrangian for the scalar sector of the minimal composite Higgs model [4], as expected. The expressions obtained for L 1 and L 2 are remarkably compact and a similar pattern holds for higher orders in 1/λ. The maximum number of derivatives of L n is 2 + 2n, although not all 2 + 2n derivative operators are generated at order n. This is as foreseen, as for large λ the non-linear regime is approached and ordering the operators by their 1/λ dependence does not coincide with the ordering given by mass dimensions. The ordering in which the operators appear is akin to the power counting of non-linear Higgs effective theory [6,[28][29][30].
Equations (18) and (21) suggest interesting correlations between operators involving the Higgs boson, gauge bosons and fermions. In particular, operators such as (∂ μ h) 2ψ ψ or V μ V μ ψ ψ, where ψ denotes a generic fermion, are weighted by the fermion mass and also bear a dependence on the SM embedding into SO (5), parametrized by the set {n, m} in Eq. (11). From those equations emerges the lowenergy effective Lagrangian in terms of SM fields at a given order in 1/λ, where the first two terms in the right-hand side contain, respectively, the kinetic terms for gauge bosons and fermions as in Eqs. (2) and (5), and the index i runs over all operator labels and coefficient functions F i (ϕ) in Table 2. The table collects all couplings corresponding to two and four "derivatives", where plain derivatives and gauge boson insertions are counted with equal weight, as they come together in the covariant derivative. The notation/basis for the purely bosonic operators was chosen according to Ref. [26] to facilitate the comparison with a model-independent approach. From this tree-level analysis we draw the following conclusions: • Up to first order in the linear corrections, the benchmark effective Lagrangian is determined to be composed of ten operators, five of them bosonic and the rest fermionic 6 including that responsible for Yukawa couplings. The coefficients of those operators are not free but intimately correlated by the coefficient functions explicitly determined in this work, and they are shown in the table.  are also shown. The right-hand column indicates the order in 1/λ at which a given couplings first appears. The Higgs field h is defined as the excitation of the field ϕ; see Eq. (12) Operator The two latter ones involve multiple Higgs insertions and are out of present experimental reach while the strength of P 6 , which involves vertices with four gauge bosons, is already tested directly by data, although the present sensitivity is very weak [31,32]. 7 • Operators involving SM fermions have an implicit dependence on the symmetry-breaking terms in the Lagrangian -they are weighted by the fermion masses in a pattern alike to that of the Minimal Flavor Violation setup [33][34][35]. Most interestingly, the corresponding F i (ϕ) coefficients, written as a function of the {n, m} parameters, 7 These bounds can be translated for instance in m ρ 70 GeV for f ≈ 600 GeV. allow one to differentiate the expected impact of different fermionic ultraviolet completions in the literature.  Table 2, including the only one appearing at O(1/λ 3 ). • The gauge field dependence is present only through powers of V μ V μ , consistent with its exclusively scalar covariant-derivative origin; see Eqs.
(3) and (10). Other Lorentz contractions such as V μ V ν 2 would be loopinduced, and thus expected to be subleading.
• The scalar functions F i (ϕ) obtained as operator weights of bosonic couplings are in agreement with those derived in Ref. [26] in the SO(5) invariant limit, for the subset of operators identified here as benchmarks; see their Eqs. (2.5)-(2.8). Table 2 provides in addition the leading deviations due to the presence of explicit SO(5)-breaking parameters α and β.

Bosonic sector
From Eqs. (20) and (21), the potential at order 1/λ reads with minimum at cos ϕ f − α 2β 1 − β 4λ , see Eq. (13). The kinetic energy of the physical Higgs excitation h (see Eq. (12)) gets then a correction given by which is reabsorbed by a field redefinition Renormalization The four independent parameters of the scalar Lagrangian Eq.
(3), f , λ, α and β, can be expressed in terms of the following observables [15]: with the Fermi constant G F as measured from muon decay, while κ V can be extracted from deviations of the Higgs couplings to two gauge bosons, for instance m h is determined from the Higgs pole mass and m ρ could in turn be determined from future measurements of the ρ mass, identifying them, respectively, with the light and heavy mass eigenvalues of the scalar sector [15] where the plus sign refers to the heavier eigenstate. Assuming the SO(5) explicit breaking to be small, |β|/4λ 1, the mass eigenvalues read with the measured value of m h implying β 0.13. In the non-linear limit λ → ∞ the ρ field decouples from the spectrum and the scalar sector would depend on just three renormalized parameters. It is now possible to foresee the impact of the linear corrections in terms of mass dependence. Precisely because the large λ and m ρ limits are in correspondence, dimensional arguments suggest the equivalence as expansion parameter; see Eqs. (29) and (30). In other words, the linear corrections are expected to be proportional to the two small parameters β and ξ and thus doubly suppressed.
Extending the renormalization scheme to the gauge sector, we choose the two extra observables needed to be the mass of the Z boson and the fine structure constant, with M Z and α em as determined from Z-pole mass measurements and from Thompson scattering, respectively [36]. In terms of the ensemble of renormalized parameters discussed above, predictions can now be made. For instance the relation between the gauge boson masses remains the same than that in the SM, where the weak angle is given at tree-level by As another example, the prediction for the Higgs→ Z Z width is modified with respect to the SM expectation by A generic expectation is the departure of κ V from 1. Indeed, the coupling between the Higgs and the gauge bosons which stems from the Lagrangian Eqs. (19)-(22) at order 1/λ is that encoded in the operator P C in Table 2 and reads (36) or in other words Assuming for illustrative purposes O(ξ ) ∼ O(1/λ) and expanding up to second order in these parameters, the result simplifies to The first term on the right-hand side of this equation is the well-known correction present in non-linear scenarios [4], while the second term encodes the linear correction linked to the scale of ultraviolet completion, which in terms of physical parameters we predict to be given by where Eq. (31) has been used. Higher-order corrections are expected to be very small, as they will stem from operators with at least four derivatives. For instance, the first extra treelevel contribution to κ V is the 1/λ 2 weight of the operator P 7 in Table 2,

Fermionic sector
Consider first the case in which the fermion mass is generated by a single Yukawa operator O (n,m) Yuk, f ; see Eqs. (6) and (7). From the Lagrangian in Eqs. (19)- (22), and more specifically from the Yukawa operator in the third line of Table 2, an expression for the fermion mass follows after applying Eqs. (12) and (26), Renormalization The renormalization scheme is now enlarged to the fermion sector choosing as observable precisely the fermion masses. The prediction that follows for the Higgs coupling to a given fermion f, then takes the form Encoding the deviations with respect to the SM expectations through the conventional κ f parameter, where g SM hff = m f /v, the exact and somewhat lengthy expression for κ f up to order 1/λ follows. The latter can be simply recast assuming again O(ξ ) ∼ O(1/λ), leading to where once again Eq. (31) has been used. It is straightforward to check that the first term on the right-hand side of this equation reproduces well-known κ f results for different models in the literature, which assume a non-linear realization. The second term gives instead the leading linear corrections. For instance, this equation leads to the following results for the MCHM 5−1−1 (corresponding to n = m = 0 in our parametrization) and MCHM 5 (corresponding to n = 1, m = 0): obtaining again at order 1/λ a correction doubly suppressed as proportional to both β and ξ ; see Eq. (31). Consider next the case in which a given fermion mass corresponds to the combination of several SO(5) invariant Yukawa operators, instead of just one as developed above, where c (n,m) are related to the generators of SO(5) and the fermion embedding in a given model. The procedure is still quite straightforward. The fermion mass will be a sum of contributions similar to that in Eq. (40) weighted by the coefficients c (n,m) , and a similar combination protocol will apply to the obtention of the fermion-Higgs coupling g hff and κ f . As an example, consider the MCHM 14−14−10 scenario [8], in which the third family quark doublet and the right-handed top are embedded each in a 14-plet of SO (5), denoted Q L and U R , respectively, while the right-handed bottom is included in a 10-plet representation denoted D R . Two SO(5) invariant operators [8] contribute in this case to the top quark mass, In contrast, in this same scenario only one effective Yukawa operator contributes to the bottom quark mass, and consequently All O(1/λ) corrections considered above show again the double suppression in ξ and β, which after Eq. (31) is tantamount to a m 2 h /m 2 ρ suppression factor, as expected.

Explicit fermion sector
In the previous section, the infinite mass limit for the heavyfermion sector was assumed from the start, while the corrections due to the heavy scalar singlets were explored. In this section we start instead of a complete (bosons plus fermions) renormalizable model, so as to estimate the impact of a fermionic ultraviolet completion beyond that related to the Yukawa couplings discussed earlier. The low-energy effective Lagrangian made out of SM fields will then be explicitly determined up to the leading corrections stemming from the heavy scalar and fermion sectors: respectively, up to where M i denotes generically the heavy-fermion masses.
The details of the fermion mass Lagrangian are quite model-dependent and derived from the specific SO(5) embedding of the light and heavy fermions. Many choices of fermion representations are possible. In addition to the use of heavy vectorial representations, a common trend is to avoid by construction direct Yukawa couplings of the Higgs field to the SM fermions, leading to a generalized see-saw pattern for light fermions with masses inversely proportional to those for the heavy fermions.
The fermionic Lagrangian L f in Eq. (2) needs to be redefined, where L kin f contains now kinetic terms for all fermions, light and heavy, and the fermion mass Lagrangian denoted by L Yuk f needs to be specified for a particular ultraviolet fermion completion. The model developed in Ref. [15] will be analyzed as an illustration, recalling first its main ingredients. In order to obtain the correct hypercharge assignments, the symmetry of the Lagrangian needs to be enlarged as customary to SO(5) × U (1) X which is broken down to R + X . The fermion fields that will generate the top mass are ψ (2/3) ∼ (X, Q, T (5) ) ∼ (2 +7/6 , 2 +1/6 , 1 +2/3 ), where ψ (x) and χ (x) belong, respectively, to the 5 and 1 representations of SO(5) with U (1) X charge x; their decomposition in terms of SU (2) L × U (1) Y charges is also shown. This choice of heavy-fermion representations corresponds to the MCHM 5−1−1 scenario, that is, to the entry 5-1 in the first row of Table 1, and thus to the effective Yukawa operator O (n,m) Yuk,f in Eq. (7) with {n, m} = {0, 0}. The fermionic Lagrangian for that field content reads where (x) n×m are spurion fields that break explicitly SO (5), while all other terms are SO (5) invariant. In terms of SU (2) L fields, the Yukawa and spurion terms in Eq. (52) provide now an explicit realization of the Lagrangian L Yuk f in Eq. (51): In Ref. [15] we had first integrated out the heavy fermions of this Lagrangian, determining then the effective Lagrangian made out of SM fields plus the singlet scalar present in the minimal SO(5) sigma model. Here we reverse the order of integration of the heavy fields, taking first the limit of heavy ρ and then that of heavy BSM fermions. We have explicitly checked that the final low-energy effective Lagrangian made out only of SM fields is independent of the order in which those limits are taken.
Using polar coordinates and integrating out the radial mode ρ does not bring about any novel complication with respect to the procedure carried out in the previous section, except for lengthier expressions. Nevertheless, L Yuk F can be compactly written prior to any integration procedure as (1) where O F s and O F c are heavy-fermion bilinears corresponding to the first four lines in Eq. (53):

L B
(1) where e + = (1, 0) and e − = (0, 1). Consider next the limit of very large scalar mass m ρ (that is, λ → ∞) and very heavy fermions. Implementing first the 1/λ corrections, the effective Lagrangian at this order takes exactly the form in Eq. (21), although ρ 1 shows now an explicit dependence on the heavy-fermion spectrum, instead of the effective dependence in Eq. (18). New operators beyond those previously considered appear, such as They are higher-order operators made out of both SM and heavy BSM fermions and related to the explicit fermionic ultraviolet completion. Furthermore, it is again easy to verify that the counting rule matches the NDA rule [6,30] by identifying λ f ∼ . Consider next the integration of the heavy-fermion sector in the results just obtained. This is an elaborated task, and the procedure and an explicit computation is described in Ref. [15]. To estimate the corrections, we adopt here a universal heavy-fermion mass scale M i associated with the mass generation mechanism of a given SM fermion, so that M 1 ∼ M 5 ∼ 1 ∼ · · · ∼ M t . Assuming this scale to be larger than f , f /M i is a good expansion parameter. The final set of five effective operators resulting up to first order in the 1/λ and f /M i expansions is shown in Table 3, where the a i σ 1 operator coefficients weighting the f /M i corrections are expected to be O(1) and their exact expressions can be found in Ref. [15]. 8 Noteworthy consequences include: • At tree level, the heavy fermions have no impact on the gauge-Higgs coupling and κ V is still given by Eq. (37). The coupling to top quarks, on the contrary, will receive fermionic contributions from the first operator in Table 3, Again, a double suppression acts on the leading heavyfermion corrections ∼ ξ f /M t , alike to the case for the bosonic ones in ∼ βξ/(2λ). It is important to note, though, that the tree-level fermionic contributions found may be larger than those induced by the scalar sector if f /M t > β/λ; this may occur specially for the top quark since the top partners, with characteristic mass scale M t , should be light enough in order not to generate a hierarchy problem. • On top of the above, higher-order effective operators involving SM fields are singled out at low scales: the dominant ones are the last three presented in Table 3. For these operators, the inclusion of an explicit heavyfermion sector does not change much the conclusions obtained previously by using an effective Yukawa coupling as defined in Eq. (7). • In the limit f /M i → 0, the operators in Table 3 coincide as expected with the fermion-Higgs and four fermion operators given previously in Table 2 using the effective Yukawa operator O (0,0) Yuk . 8 The a i σ 1 coefficients are a redefinition of the c i σ operator coefficients in Ref. [15] so as to extract explicitly the f /M i dependence: c i σ 1 → y 0 t a t σ 1 /M t ; the exact expressions for c i σ for the fermion model discussed here can be found in Table 3 of that reference. Table 3 Effective operators, up to order f /M i and 1/λ, after integrating out the radial mode ρ and the heavy fermions in a UV realization of partial compositeness. The coefficients a f σ 1 , with f = t, b can be found in Ref. [15]; see footnote 6. The Hermitian conjugate should be included for all operators here. The Higgs field h is defined as the excitation of the field ϕ; see Eq. (12) Operator

Conclusions
The linear sigma model for QCD allows one to monitor the transition from a completely renormalizable model in a weakly interacting regime to a non-linear regime in the high mass (λ → ∞) limit. In this work we have carried out an analogous exploration assuming that the Higgs particle may correspond to a pseudo-Goldstone boson of a spontaneously SO(5) global symmetry (or containing SO (5)) at high energies, completing the procedure first proposed and started in Ref. [15]. The minimal sigma model for SO (5) has been used as starting point. The results are independent of the relative order in which the high mass limit for the heavy boson ρ and for the heavy fermions are taken. In a first stage, the bosonic sector was left fully dynamical while we defined an effective Yukawa operator characterized by two parameters, which depend only on how the light SM fermion fields are embedded in representations of the SO(5) symmetry for any given model. Armed with this tool, the benchmark effective Lagrangian has been derived for the large m ρ limit, see Table 2. Up to first order in the linear corrections, it is shown to be composed of ten operators, five of them bosonic and the rest fermionic including that responsible for the usual fermion Yukawa coupling; for the fermionic operators the coefficients are given as an explicit function of the two parameters which define the effective Yukawa operator. Their simple form allows a direct prediction and comparison of the many models in the literature which differ by their fermionic embedding. It is straightforward to obtain from this result the expressions of the Higgs couplings to fermions, κ f , in general.
Among the O(1/λ) corrections, three bosonic operators have coefficients which are independent of the global symmetry-breaking mechanism and should thus be of special relevance; among them a four-gauge boson vertex is already being directed probed by present data, while the other two involve vertices with at least two Higgs fields. We have also proved that the leading phenomenological couplings of the Higgs particle to gauge bosons and SM fermions, κ V and κ f , are quite universal, the reason being that the linear corrections must be doubly suppressed as proportional to both 1/λ and to the explicit symmetry-breaking parameters, in a combination corresponding to a m 2 h /m 2 ρ suppression. The tower of higher-order operators obtained is shown to correspond to only a small subset of the most general nonlinear Lagrangian [26] for a generic non-linear realization of electroweak symmetry breaking. It is also consistent with the results for general SO(5)/SO(4) constructions in Ref. [26], singling out a fraction of operators found in the latter, with its expected coefficients. The minimal set identified here could serve to focus model-independent searches of a dynamical nature for the Higgs particle.
In a second stage, we have explored a complete renormalizable model with an explicit heavy-fermion ultraviolet completion and repeated the integration procedure; the leading corrections stemming both from the heavy-boson and from the heavy-fermion sector were then identified and their coefficients determined. New higher-order operators related to the specific fermion ultraviolet completion, made out of SM and heavy fermions and containing vertices with at least four fields, were identified as an intermediate step. Finally, the set of operators made out of SM fields and involving fermions was determined to consist of only five operators at the lead-ing order in both expansions: one coupling contains the usual SM Yukawa coupling, two are fermion-boson operators and the remaining two correspond to four fermion couplings. The results match those obtained in the first part using the effective Yukawa operator; interestingly, they show that deviations to the SM value for κ f due to tree-level exchange of heavy fermions may dominate over those stemming from the scalar (e.g. 1/λ) linear corrections.
The starting point of the analysis in this work is a minimal renormalizable sigma model for SO (5). Other renormalizable -more complicated -realizations are conceivable, in the same way that the linear sigma model for QCD could be extended. Although additional effective operators could be sourced in such constructions [26], the operators identified here are expected to be the tell-tale of a Goldstone-boson origin for the Higgs field and as such common to all realizations.