A Golden Probe of Nonlinear Higgs Dynamics

TL;DR

This paper performs the first fit of Wilson coefficients in nonlinear Higgs dynamics using LHC data from the golden H to 4L decay channel, constraining the decay constant and exploring the structure of the Higgs sector.

Contribution

It introduces a novel analysis combining rate and shape information to constrain Wilson coefficients and the decay constant in nonlinear Higgs models using LHC Run 2 data.

Findings

Current LHC data constrains the decay constant $f$ and Wilson coefficients.

Including shape information broadens the allowed parameter space.

Future high-luminosity LHC could improve constraints and enable multi-coefficient fits.

Abstract

The most salient generic feature of a composite Higgs boson resides in the nonlinearity of its dynamics, which arises from degenerate vacua associated with the pseudo-Nambu-Goldstone (PNGB) nature of the Higgs boson. It has been shown that the nonlinear Higgs dynamics is universal in the IR and controlled only by a single parameter $f$, the decay constant of the PNGB Higgs. In this work we perform a fit, for the first time, to Wilson coefficients of ${\cal O}(p^4)$ operators in the nonlinear Lagrangian using the golden H $\to$ 4L decay channel. By utilizing both the "rate" information in the signal strength and the "shape" information in the fully differential spectra, we provide limits on the Goldstone decay constant $f$, as well as ${\cal O}(p^4)$ Wilson coefficients, using Run 2 data at the LHC. In rate measurements alone, the golden channel prefers a negative $\xi =v^2/f^2$ corresponding to a non-compact coset structure. Including the shape information, we identify regions of parameter space where current LHC constraint on $f$ is still weak, allowing for $\xi \lesssim 0.5$ or $\xi \gtrsim -0.5$. We also comment on future sensitivity at the high-luminosity upgrade of the LHC which could allow for simultaneous fits to multiple Wilson coefficients.

Tables

Table 1. Table 1: Summary of CMS rate measurements of H → → \to 4L from ggh+ttH production channels Khachatryan  et al. ( 2015 ); Sirunyan  et al. ( 2019b ) . Each measurement is optimized with respect to a particular anomalous HVV coupling as indicated. However, the fit is performed assuming A i V = 0 superscript subscript 𝐴 𝑖 𝑉 0 A_{i}^{V}=0 .

${\mu }_F^{\text{4L}}\left(A_i\right)$	Best Fit at $A_i=0$	$95\%$ C.L. Sirunyan et al. (2019b)
${\mu }_F^{\text{4L}}\left(A_2^{ZZ}\right)$	1.19	$[0.857,1.602]$
${\mu }_F^{\text{4L}}\left(A_4^{ZZ}\right)$	1.26	$[0.907,1.652]$
${\mu }_F^{\text{4L}}\left(A_4^{Z\gamma }\right)$	1.24	$[0.907,1.612]$

Table 2. Table 2: Summary of CMS shape measurements on ratios of anomalous HVV couplings Khachatryan  et al. ( 2015 ); Sirunyan  et al. ( 2019b ) . In the middle column we indicate the assumptions on the Higgs total width.

${ℛ}_i=A_i/A_1^{ZZ}$	Constraint	$95\%$ CL
${ℛ}_2^{ZZ}$	${\Gamma }_H={\Gamma }_H^{\mathrm{SM}}$	$[-0.381,0.180]$ Sirunyan et al. (2019b)
${ℛ}_4^{ZZ}$	${\Gamma }_H={\Gamma }_H^{\mathrm{SM}}$	$[-0.272,0.129]$ Sirunyan et al. (2019b)
${ℛ}_4^{Z\gamma }$	On-shell Events	$[-0.792,0.287]$ Sirunyan et al. (2019b)

Equations

$$\Lambda \sim 4\pi f,$$

$$S_{nl\sigma m}=\int d^{4}x{\Lambda }^2{f}^2\mathcal{L}\left(\frac{\pi }{f},\frac{\partial }{\Lambda }\right)=\int d^{4}x\left[{\mathcal{L}}^{(2)}+{\mathcal{L}}^{(4)}+\cdots \right],$$

$${\mathcal{L}}^{(2)}$$

$$m_W=\frac{m_Z}{\cos {\theta }_W}=\frac{1}{2}gv=\frac{1}{2}gf\sin \theta .$$

$$\frac{1}{2}{\partial }_{\mu }h{\partial }^{\mu }h+\left[2\sqrt{1-\xi }\frac{h}{v}+(1-2\xi )\frac{h^2}{v^2}+\cdots \right]\left(m_W^2{W}_{\mu }^{+}{W}^{-\mu }+\frac{1}{2}{m}_Z^2{Z}_{\mu }{Z}^{\mu }\right).$$

$$S_{SILH}=\int d^{4}x{M}_{\rho }^2{f}^2\mathcal{L}\left(\frac{\pi }{f},\frac{\partial }{M_{\rho }}\right),$$

$$\frac{v}{f}\lesssim 1\text{and}\frac{E}{\Lambda }\ll 1,$$

$$S_{SILH}^{(4)}=\int d^{4}x{M}_{\rho }^2{f}^2{\mathcal{L}}^{(4)}\left(\frac{\pi }{f},\frac{D}{M_{\rho }}\right)=\int d^{4}x\sum _i\frac{c_i}{g_{\rho }^2}{O}_i,$$

$${\mathcal{L}}^{(1h)}$$

$$C_1^h$$

$$C_2^h$$

$$C_3^h$$

$$C_4^h$$

$$C_5^h$$

$$C_6^h$$

$${\Gamma }_V^{\mu \nu }=\frac{1}{v}\left[A_1^V{m}_Z^2{g}^{\mu \nu }+A_2^V(k_1^{\nu }{k}_2^{\mu }-k_1\cdot k_2{g}^{\mu \nu })+A_3^V{ϵ}^{\mu \nu \alpha \beta }{k}_{1\alpha }{k}_{2\beta }+(A_4^V{k}_1^2+{\bar{A}}_4^V{k}_2^2)g^{\mu \nu }\right],$$

$$A_1^{ZZ}$$

$$A_2^{ZZ}$$

$$A_4^{ZZ}$$

$$A_2^{Z\gamma }$$

$$A_4^{Z\gamma }$$

$$\frac{d{\Gamma }_{\text{H→4L}}}{d\mathcal{O}}=\sum _{ij}{A}_i{A}_j^{\ast }\times \frac{d{\hat{\Gamma }}_{ij}}{d\mathcal{O}},$$

$${\hat{\Gamma }}_{ij}=\int \frac{d{\hat{\Gamma }}_{ij}}{d\mathcal{O}}d\mathcal{O},$$

$${\Gamma }_{\text{H→4L}}=\sum _{ij}{A}_i{A}_j^{\ast }\times {\hat{\Gamma }}_{ij},$$

$${\mathcal{R}}_{\text{4L}}\equiv \frac{{\Gamma }_{\text{H→4L}}}{{\Gamma }_{\mathrm{H}\to 4L}^{\text{SM}}}=\sum _{ij}{A}_i{A}_j^{\ast }\times \frac{{\hat{\Gamma }}_{ij}}{{\Gamma }_{\mathrm{H}\to 4L}^{\text{SM}}}.$$

$${\mathcal{R}}_{2e2\mu }$$

$${\mathcal{R}}_{4e/4\mu }$$

$${\mu }_i^{\text{4L}}=\frac{\sigma (i\to \text{H}\to \text{4L})}{\sigma (i\to \text{H}\to \text{4L}{)}_{\mathrm{SM}}}=\frac{\sigma (i\to \text{H})}{\sigma (i\to \text{H}{)}_{\mathrm{SM}}}\times \frac{{\Gamma }_{\text{H}}^{\mathrm{SM}}}{{\Gamma }_{\text{H}}}\times \frac{{\Gamma }_{\text{H→4L}}}{{\Gamma }_{\mathrm{H}\to 4L}^{\text{SM}}},$$

$${\mu }_F^{\text{4L}}=\frac{{\Gamma }_{\text{H→4L}}}{{\Gamma }_{\mathrm{H}\to 4L}^{\text{SM}}}={\mathcal{R}}_{4L}.$$

$${\mathcal{R}}_i^V\equiv A_i^V/A_1^{ZZ},$$

$${\mu }_F^{\text{4L}}({\mathcal{R}}_i)$$

$${\mu }_F^{4L}(C_i)$$

$${\mu }_{Z\gamma }=\frac{\Gamma (\text{H}\to \text{Z}\gamma )}{\Gamma (\text{H}\to \text{Z}\gamma {)}_{SM}}\simeq {\left|1-1.1\frac{(1\text{TeV}{)}^2}{M_{\rho }^2}{C}_4^h\right|}^2,$$

$$\frac{M_{\rho }}{\sqrt{|C_4^h|}}>0.56(0.84)\text{TeV},$$

$$\frac{M_{\rho }}{\sqrt{|C_4^h|}}>2.2(2.5)\text{TeV}$$

Full text

1. I Introduction

2. II Universal Higgs Nonlinearity

3. II.1 ${\cal O}(p^{4})$ Operators

4. II.2 H$\to$ 4L Tensor Structure

5. III Experimental Constraints

6. III.1 The Golden Channel: H $\to$ 4L

7. III.1.1 Constraints From Rate Measurements

8. III.1.2 Projections with Multi-dimensional Parameter Likelihood Fits

9. III.2 Direct H$\,\to\,$Z$\gamma$ Decay

10. III.3 Precision Electroweak Constraints

11. IV Conclusions

I Introduction

The Higgs boson plays a central role in our understanding of physics at the electroweak scale, the energy scale being probed by the Large Hadron Collider (LHC). Precise measurement of its property is among the top priorities of experimental programs at the LHC as well as any possible future electron-positron and hadron colliders. A major goal of these efforts is to understand the microscopic nature of the 125 GeV Higgs boson:

• •

Is the Higgs boson a fundamental particle like the electron or a composite particle such as the pion?

Such a question could carry far reaching implications in our understanding of the Universe at the most fundamental scale, especially given the fact that no other fundamental scalar particles have been observed in nature. While current data at the LHC is consistent with the expectation of a Standard Model (SM) Higgs boson, the experimental uncertainty is still sizeable, on the order of 10 - 20% or more Aad et al. (2016, 2019); Sirunyan et al. (2019a). Such a large uncertainty can hardly be characterized as “precise” and the question of whether the 125 GeV Higgs boson is indeed the SM Higgs remains open.

There is a long history in the idea of a composite Higgs boson, dating back to the classic papers Kaplan and Georgi (1984); Kaplan et al. (1984) several decades ago, where the Higgs boson arises as a pseudo-Nambu-Goldstone boson (PNGB) like the pion in low-energy QCD. Models with a PNGB Higgs were revived and refined much later, via the little Higgs theories Arkani-Hamed et al. (2001, 2002a, 2002b) and the holographic Higgs models Contino et al. (2003); Agashe et al. (2005). By now they are collectively referred to as the composite Higgs models.

There are numerous composite Higgs models Bellazzini et al. (2014) which differ in several aspects. For example, the choice of symmetry breaking pattern $G/H$ where $G$ is the broken group in the UV and $H$ is the unbroken group containing the electroweak $SU(2)_{L}\times U(1)_{Y}$ gauge groups. The Higgs arises as a PNGB when $G$ is spontaneously broken to $H$ at a scale $\Lambda\sim 10$ TeV. One other aspect that varies greatly is the implementation of additional fermions associated with the third generation quarks in the SM, which are introduced to reduce the UV sensitivity in the Higgs mass originating from the SM top quark. In some cases the new fermions do not even have to carry SM color charge Chacko et al. (2006) or electroweak quantum numbers Cheng et al. (2018); Cohen et al. (2018).

In spite of all these variations in model-building, there is one salient prediction that is generic to the entire class of composite Higgs models:

• •

Nonlinear dynamics in Higgs interactions with the electroweak gauge boson, as well as self-interactions carrying derivatives.

The origin of nonlinear Higgs interactions goes to the essence of a composite Higgs: the PNGB nature of the Higgs boson. It is the same nonlinear dynamics appearing in interactions of pions in chiral symmetry breaking, or any other Nambu-Goldstone bosons observed in nature. The theoretical tool employed to construct effective Lagrangians of Nambu-Goldstone bosons was developed half a century ago, in the seminal papers by Callan, Coleman, Wess and Zumino (CCWZ) Coleman et al. (1969); Callan et al. (1969). In the CCWZ approach each symmetry breaking pattern $G/H$ results in a seemingly different effective Lagrangian, each with its own set of nonlinear interactions and experimental predictions.

Only in recent years was it realized that the nonlinear interaction of a PNGB has very little to do with the details of the UV group $G$ that is being spontaneously broken. This can be seen either by imposing shift symmetries on the Nambu-Goldstone bosons in the IR Low (2015a, b), or by soft bootstrapping tree-level amplitudes of nonlinear sigma model (nl$\sigma$m) Cheung et al. (2015, 2017). More specifically, the nonlinear interaction of Nambu-Goldstone bosons owes its presence to the existence of degenerate vacua in the deep IR and knows very little about the details of the broken Group $G$ Low and Yin (2018). The only information on $G$ resides in the overall normalization of the Goldstone decay constant $f$, which can be taken as an input parameter in the low-energy.

Following this progress, the complete list of modifications to the couplings of one Higgs boson to two electroweak gauge bosons (HVV), two Higgs bosons with two electroweak gauge bosons (HHVV), one Higgs coupled to three electroweak gauge bosons (HVVV), as well as triple gauge boson couplings (TGC), were studied in Refs. Liu et al. (2018a, b) up to four-derivative order and all orders in $1/f$. These corrections are the universal predictions of a composite Higgs boson. In particular, a set of “universal relations” among the couplings were proposed which depend on only one input parameter $f$. Experimental confirmation of these universal relations would be a striking signal of the PNGB nature of the 125 GeV Higgs.

In this work we continue with the exploration of universal Higgs nonlinearity and its implications in Higgs coupling measurements. In particular, we will study the possibility of measuring and constraining Wilson coefficients of four-derivative operators in the effective Lagrangian of a composite Higgs boson using the H $\to$ 4L decay channel, the so-called “Golden channel” because it is the clearest and cleanest among all Higgs decay channels. The small background contamination combined with the availability of full kinematic distributions of decay products offers not only a powerful probe of the spin and CP property of the Higgs boson Choi et al. (2003); Buszello et al. (2004), but also a unique opportunity to employ advanced multivariate techniques Gainer et al. (2011); Stolarski and Vega-Morales (2012); Chen et al. (2013); Bolognesi et al. (2012); Chen and Vega-Morales (2014); Chen et al. (2015a, 2014a); Gainer et al. (2015); Chen et al. (2015b, c, 2016); Chen and Vega-Morales (2017) to enhance the experimental sensitivity 111Analyses of Higgs couplings based on information geometry have also been shown to be powerful Brehmer et al. (2017, 2018).. We will use the publicly available LHC analyses in the 4L channel to probe and constrain, for the first time, Higgs nonlinear dynamics and the associated Wilson coefficients. Furthermore, we demonstrate that the conventional approach of bounding the decay constant $f$ using the signal strength (total rate) in HVV measurements is incomplete, especially when effects of ${\cal O}(p^{4})$ operators are included.

This work is organized as the follows. In Section II we briefly review the universal Higgs nonlinearity, listing all operators modifying the HVV couplings up to ${\cal O}(p^{4})$. We also map out the correspondence to the tensor structure basis used in the LHC analyses. Then in Section III we study the experimental constraints from both the rate and shape measurements in the H$\to$ 4L channel, as well as in H$\to$ Z$\gamma$ two body decays. Future sensitivity projections at the high-luminosity (HL) LHC are also provided in this section as well as a brief study of precision electroweak constraints. Finally we conclude in Section IV.

II Universal Higgs Nonlinearity

The effective Lagrangian of a PNGB Higgs boson contains two expansion parameters,

• •

Nonlinear expansion characterized by ${\pi}/{f}$, where $\pi$ denotes a generic Nambu-Goldstone field like a composite Higgs boson. The expansion in $\pi/f$ is highly nonlinear for a PNGB Higgs boson.

• •

Derivative expansion chararcterized by ${\partial_{\mu}}/{\Lambda}$. If some of the unbroken symmetry is gauged, one replaces the ordinary derivative by the gauge covariant derivative, $\partial_{\mu}\to D_{\mu}=\partial_{\mu}-iA_{\mu}$. In this power counting the gauge field $A_{\mu}$ then counts as one derivative.

This is similar to the chiral Lagrangian in low-energy QCD, which describes interactions of pions as PNGB’s arising from spontaneously broken $SU(2)_{L}\times SU(2)_{R}$ chiral symmetry Georgi (1984). In the chiral Lagrangian $f=f_{\pi}\sim 93$ MeV is the pion decay constant and $\Lambda=\Lambda_{\chi SB}\sim 1$ GeV is the scale where QCD becomes strongly coupled and chiral symmetry is spontaneously broken. In naive dimensional analysis (NDA) Manohar and Georgi (1984), requiring quantum corrections from higher loops to be comparable in size to the derivative expansion in $\partial/\Lambda$ leads to the relation

[TABLE]

which is saturated in a strongly interacting theories with no relevant small adjustable parameters. Then the effective Lagrangian of Nambu-Goldstone bosons based on a nonlinear sigma model (nl$\sigma$m) is organized as follows

[TABLE]

where ${\cal L}^{(n)}$ represent operators containing $n$ derivatives. In this work we focus on $n\leq 4$. While ${\cal L}^{(n)}$ incorporates the derivative expansion at a well-defined order, the nonlinear expansion is resummed to all orders in $\pi/f$ in ${\cal L}^{(n)}$. Conventional wisdom from the CCWZ approach Coleman et al. (1969); Callan et al. (1969) has it that the nonlinear expansion realizes the spontaneously broken symmetries in the UV. The summation of operators to all orders in $\pi/f$ and only at a certain order in $\partial/\Lambda$ stems from requiring the resummed operators to have well-defined transformation properties under the broken symmetry. The nonlinear Higgs dynamics then follows. However, the CCWZ perspective obscures the universality in the Nambu-Goldstone interactions. It turns out that these nonlinear interactions arise entirely from the presence of degenerate vacua in the IR and are insensitive to the coset structure $G/H$ in the UV Low (2015a, b); Low and Yin (2018); Liu et al. (2018a). The only parameter dependent on the coset structure is the normalization of $f$. In the end, by probing and measuring the Higgs nonlinear dynamics, one could potentially gain insight on a broad range of composite Higgs models, regardless of the $G/H$ coset.

Invoking the custodial invariance and focusing on the scenario where the 125 GeV Higgs transforms as a fundamental representation of an unbroken $SO(4)$ group, the leading two-derivative Lagrangian is simple, upon gauging an $SU(2)_{L}\times U(1)_{Y}$ subgroup of $SO(4)$ Liu et al. (2018a, b),

[TABLE]

where $h$ is the 125 GeV Higgs boson. In the above $\sin\theta\equiv v/f$, where $v=246$ GeV, is the vacuum misalignment angle and $W^{\pm}$ and $Z$ are the electroweak massive gauge bosons, whose masses can be read off from Eq. (3):

[TABLE]

As emphasized previously, ${\cal L}^{(2)}$ contains two-derivative operators that are to all orders in $h/f$. In particular, coefficients of operators at different orders in $h/f$ are all predetermined by a single parameter $f$ and resummed to the form $\sin(\theta+h/f)$. Although not transparent at all, this particular form is enforced by four ”shift symmetries” acting on the four components of the scalar $SU(2)_{L}$ doublet containing the neutral scalar $h$ Low (2015b).

Expanding Eq. (3) in $h/v$, one obtains corrections to the SM HVV and HHVV couplings at the two-derivative level,

[TABLE]

where $\xi\equiv v^{2}/f^{2}$. Given that $f$ is an input parameter whose normalization is UV dependent, Eq. (5) is the basis for extracting $f$ using the signal strength in HVV coupling measurements.

In a strongly coupled theory there are typically composite resonances below the cutoff scale $\Lambda$, as is evident in low-energy QCD. In composite Higgs models these resonances are represented by the scale $M_{\rho}=g_{\rho}f$, where $1\lesssim g_{\rho}\lesssim 4\pi$. Assuming there is a mass gap between $M_{\rho}$ and $f$, one could further integrate out the composite resonances and the resulting effective Lagrangian, which is valid at energy $E\sim f$, inherits the nonlinear interactions from above the scale $M_{\rho}$ Giudice et al. (2007). This effective theory is often referred to as the SILH Lagrangian,

[TABLE]

which has the same nonlinear structure as in Eq. (2), but with $\Lambda\to M_{\rho}$. In particular, at the two-derivative level, Eq. (3) remains unchanged even after integrating out $M_{\rho}$, and corrections to HVV and HHVV couplings are also unchanged. There are effects that could potentially spoil the nonlinear structure of the effective Lagrangian. These are related to “explicit” symmetry breaking effects and examples include the Higgs potential and the Higgs coupling to fermions. However, as argued in Ref. Liu et al. (2018a), they would modify the nonlinearity HVV and HHVV couplings only at the loop-level.222These loop effects are non-universal and beyond the scope of current work. The typical spectrum in a composite Higgs model is displayed in Fig. 1.

II.1 ${\cal O}(p^{4})$ Operators

The two-derivative nonlinear Lagrangian in Eq. (3) is valid when

[TABLE]

where $E$ is the typical energy scale probed by the experiment. In particular, resumming contributions to all orders in $v/f$ allows for an $f$ not too far above the weak scale set by $v\approx 246$ GeV. We can further extend the range of validity of the effective action by including higher order terms in the derivative expansion, such as ${\cal O}(p^{4})$ operators. Chiral Lagrangian operators that are ${\cal O}(p^{4})$ and to all orders in $\pi/f$ have been enumerated in Ref. Gasser and Leutwyler (1984). In the context of composite Higgs bosons the relevant four-derivative operators for Higgs couplings were given in Refs. Liu et al. (2018a, b) in the unitary gauge.

More specifically, the SILH Lagrangian at ${\cal O}(p^{4})$, including the full nonlinearity structure, can be written as

[TABLE]

where $c_{i}$ are expected to be order unity constants parameterizing the incalculable UV physics at the scale $\Lambda\sim 4\pi f$. In some cases operators contributing to couplings of neutral particles and an on-shell photon are further suppressed by additional loop factors. In total there are 7 operators labelled by $O_{1},O_{2},O_{3},O_{4}^{\pm}$ and $O_{5}^{\pm}$, each with the corresponding unknown Wilson coefficients $c_{i}$. These 7 Wilson coefficients contribute to a dozen different observables that can be measured in HVV and HHVV couplings Liu et al. (2018a, b). Focusing on those relevant for HVV couplings we have, in the unitary gauge,

[TABLE]

where ${\mathcal{D}}=\partial^{\mu}\partial^{\nu}-\eta^{\mu\nu}\partial^{2}$. Although there are six coefficients $C_{i}^{h},i=1,\cdots,6$ in the unitary gauge, they are secretly related by only five Wilson coefficients residing in Eq. (8), as well as the input parameter defined by $\sin\theta=\sqrt{\xi}=v/f$, in a composite Higgs model:

[TABLE]

In the above $c_{w},c_{2w},t_{w}$ denote $\cos\theta_{W},\cos 2\theta_{W},\tan\theta_{W}$ respectively, where $\theta_{W}$ is the weak mixing angle.

II.2 H$\to$ 4L Tensor Structure

The operators listed in Eq. (9) enter into H$\to$ 4L decays and manifest themselves through kinematic distributions of the decay product. For our analysis we utilize the “Golden Channel” analysis framework developed in Refs. Gainer et al. (2011); Chen et al. (2013); Chen and Vega-Morales (2014); Chen et al. (2015a, 2014a, b, c) which parametrizes the effective Higgs boson couplings to pairs of neutral vector boson pairs in terms of the Lorentz tensor structures,

[TABLE]

where $V=(ZZ,Z\gamma,\gamma Z,\gamma\gamma)$. We assume massless leptons, using the notation and conventions defined in Chen et al. (2014a), but have also included the $A_{4}^{V}$ tensor structures which were not included previously. Note that electromagnetic gauge invariance requires $A_{1}^{\gamma\gamma}=A_{1}^{Z\gamma}=A_{4}^{\gamma\gamma}=0$. In general the $A_{i}^{V}$ can be momentum dependent form factors which are functions of Lorentz invariant products of the external momenta. For our purpose it is sufficient to take them to be real and constant coefficients in which case we have $A_{4}=\bar{A}_{4}$. Furthermore, in the massless lepton case we have $Z^{\mu}\partial^{\nu}V_{\mu\nu}=Z_{\mu}\mathcal{D}^{\mu\nu}V_{\nu}$. Then the relation between $A_{i}^{V}$ and the Wilson coefficients defined in Eq. (9) is quite simple:

[TABLE]

For simplicity we have kept only CP-even terms and neglected effective couplings to pairs of photons in Eq. (16), which do not appear in universal Higgs nonlinearity, though it would be straightforward to include them in the fit. Note that as $\xi\to 0$ we recover the SM value for the tree level HZZ coupling, $A_{1}^{ZZ}=2$. The important observation is:

• •

Only $A_{1}^{ZZ}$ is directly related to $\xi$; the other tensor structures are dependent on both $M_{\rho}$ and $C_{i}^{h}$.

At the leading order in derivative expansion, the signal strength in HVV coupling measurements is attributed entirely to $A_{1}^{ZZ}$, which is then used to constrain $\xi$ de Blas et al. (2018). We see this is not the case anymore when ${\cal O}(p^{4})$ effects are included. It turns out that all couplings in Eq. (II.2) can be extracted from the fully differential spectra of H$\to$ 4L decays. In what follows we use CMS H$\to$ 4L data to constrain the $A_{i}^{V}$ coefficients, which then translate into limits on $\xi$ and the Wilson coefficients $C_{i}^{h}$.

III Experimental Constraints

III.1 The Golden Channel: H $\to$ 4L

In this subsection we will use both the “rate information,” which pertains to the signal strength measured in H$\to$ 4L channel, and the “shape information,” as contained in the differential spectra of final state leptons, to place constraints on the Wilson coefficients in the nonlinear Higgs Lagrangian in Eq. (9). In particular, we utilize the fully differential decay width analytically computed in Chen et al. (2013); Chen and Vega-Morales (2014) for H $\to 2e2\mu,4e$, and $4\mu$ assuming on-shell decay of the Higgs boson. Interference effects between the different tensor structures in Eq. (16), as well as among identical final states in the case of $4e$ and $4\mu$, have been fully accounted for. Before examining current CMS constraints, we briefly review the H$\to$ 4L partial decay width and fully differential spectra.

For our purpose it is convenient to single out the $A_{i}^{V}$ dependence in the differential decay width, which can be written schematically as

[TABLE]

where $\mathcal{O}$ represents all observables available in the H$\to$ 4L decay channel Chen et al. (2013); Chen and Vega-Morales (2014) and $i,j$ sum over all of the possible tensor structures in Eq. (16). The advantage of doing so is then the remaining quantities ${d\hat{\Gamma}_{ij}}/{d\mathcal{O}}$ can be calculated and integrated over a particular phase space. We can then define the ‘sub-widths’ for each combination of $A_{i}A_{j}^{\ast}$ as

[TABLE]

which is just a numerical constant once a selection cut over the phase space is chosen. It is worth emphasizing that the sub-widths could be negative for certain combinations of tensor structures when they interfere destructively. However, the partial width written as the sum

[TABLE]

must be positive and is now a function of the effective couplings and the phase space cuts. Eq. (18) and Eq. (20) allow for a full reconstruction of the differential decay spectra and the partial width of H$\to$ 4L, respectively.

It will be convenient to normalise $\Gamma_{\rm H\to 4L}$ to the (tree level) SM expectation, which corresponds to $A_{1}^{ZZ}=2$ and all other couplings set to zero,

[TABLE]

Performing the integration over phase space we obtain the normalized partial width in Eq. (21) for the $2e2\mu$ and $4e/4\mu$ channels respectively,

[TABLE]

with cuts and reconstruction corresponding to ‘CMS-like’ phase space selections Khachatryan et al. (2015); Sirunyan et al. (2019b, 2017). In this work we have included the $A_{4}^{ZZ/Z\gamma}$ couplings, which were not included in previous studies Chen et al. (2014a, 2015c). For completeness we have included all of the operators in Eq. (16), though below we will focus on those relevant for the nonlinear Higgs dynamics we are interested in. Note terms linear in the CP-odd couplings (${A}_{3}^{VV}$) do not appear because they integrate to (nearly) zero when choosing CMS-like selection cuts, which reflects the fact that “rate” measurements are not sensitive to CP violation. In this regard, one could employ the shape information in the differential spectra Chen et al. (2014a, 2015b, 2015c) or the construction of forward-backward asymmetries Chen et al. (2014b) to probe CP violation, but we do not explore this possibility here.

In the CMS analyses in Refs. Khachatryan et al. (2015); Sirunyan et al. (2019b), the signal strength in H$\to$ 4L is measured in two categories, depending on whether the production channels involve Higgs couplings to fermions (ggH and ttH channels) or Higgs couplings to electroweak bosons (VBF and VH channels),

[TABLE]

where $i=F,V$ represents the ggH+ttH channels and VBF+VH channels, respectively. In the following we will focus on the fermionic production channels, ggH+ttH, for two reasons: 1) ggH is by far the dominant production channel of the 125 GeV Higgs at the LHC and the uncertainty is smaller and 2) any modification in HVV couplings will enter into both the production and decay amplitudes in the VBF+VH channels. For simplicity we assume the production cross-sections in ggH+ttH channels are equal to their SM values as well as the total Higgs width. Under these assumptions the deviation in “rate measurements” in $\mu_{F}^{\text{4L}}$ arises entirely from the decay amplitudes:

[TABLE]

For completeness we briefly summarize the procedures taken by the CMS collaboration in Refs. Khachatryan et al. (2015); Sirunyan et al. (2019b), which we refer the reader to for details. For CP-even couplings, which we focus on in this work, CMS built three multi-variate likelihood functions, each optimized for a particular anomalous Higgs coupling in Eq. (16): $A_{2}^{ZZ},A_{4}^{ZZ}$ and $A_{4}^{Z\gamma}$. In each likelihood function all anomalous couplings, other than the one the likelihood function is specifically optimized for, are set to zero Gritsan (ions). In other words, in these analyses the anomalous HVV couplings are turned on only one at a time. (See Refs. Chen and Vega-Morales (2014); Chen et al. (2015a) for a framework that allows for simultaneous measurements of several anomalous Higgs couplings at the same time.) The likelihood function allows one to constrain and fit: 1) the rate (signal strength) in H$\to$ 4L decays and 2) the “fraction” of the observed 4L events originated from the anomalous HVV coupling. It turns out that the best fit value for the three CP-even anomalous HVV couplings analyzed in Refs. Khachatryan et al. (2015); Sirunyan et al. (2019b) all have central values extremely close to zero, albeit with varying degrees of uncertainties. As a result, CMS provided four different fits to the signal strength in 4L channel, under the assumption of vanishing anomalous HVV couplings. The relevant CMS results are summarized in Table 1, where we have computed the 95% C.L. from the 67% C.L. by assuming a Gaussian distribution. As for the shape measurements, CMS provided limits on the ratio 333Note we use a slightly different normalization than in the CMS analysis Khachatryan et al. (2015); Sirunyan et al. (2019b).

[TABLE]

from the differential spectra in the 4L decays. This is possible because different tensor structures in Eq. (16) result in different shapes in the differential distributions. We quote the results in Table 2 and give a visual representation of the bounds in Fig. 2.

In their Run 2 analyses CMS did not provide a fit to $A_{2}^{\gamma\gamma/Z\gamma}$ Khachatryan et al. (2015); Sirunyan et al. (2019b), because these are better constrained from direct H$\to\gamma\gamma$ and H$\to$Z$\gamma$ two body decays. It turns out that only $A_{2}^{Z\gamma}$ is universal in nonlinear Higgs dynamics, as $A_{2}^{\gamma\gamma}$ requires shift-symmetry breaking effects. In Section III.2 we perform a separate fit to $A_{2}^{Z\gamma}$ using the direct Z$\gamma$ channel.

It is worth emphasising that, in order to properly constrain nonlinear Higgs dynamics, one would need to adapt the present CMS analyses. In particular, we would like to

fit the magnitude of the $A_{i}^{V}$ coefficient in Eq. (16), instead of just the fraction. 2. 2.

allow all anomalous HVV couplings to be present at the same time.

It would be of interest to study the impact on the signal strength ($\mu_{F}$) measurement by relaxing the assumption of vanishing anomalous HVV couplings. In what follows we will perform a limited analysis within the framework of the CMS analyses, and await a more comprehensive experimental study in the future.

III.1.1 Constraints From Rate Measurements

Using the CMS bounds given in Table 1 and Table 2 we can now study limits derived from rate measurements. We restrict ourselves to the CP-even HVV couplings $A_{1}^{ZZ},A_{2}^{ZZ},A_{4}^{ZZ}$ and $A_{4}^{Z\gamma}$, and trade $A_{1}^{ZZ}$ for the nonlinear parameter $\xi$. Then Eq. (24) can be written as

[TABLE]

where we have used only the $2e2\mu$ normalized partial width in Eq. (22) with similar results obtained using the $4e$ channel for the operators of interest in this study. With Eq. (26) we can perform two kinds of fits to $\xi$:

In the right-panel of Fig. 3, we present bounds on $\xi$ by setting all ${\cal R}_{i}=0$ and applying the measured $\mu_{F}^{\text{4L}}$ in Table 1, where the constraints are obtained assuming that the anomalous HVV coupling vanishes. The values of $\xi$ for the different analysis vary because different categorization and observables are utilized. 2. 2.

In the left-panel of Fig. 3, we use the central value of $\mu_{F}^{\text{4L}}$ from Table 1 in the left-hand side of Eq. (26) and plug in the experimental limits on the corresponding ${\cal R}_{i}$ from Table 2 in the right-hand side. We then invert Eq. (26) to derive a limit on $\xi$.

The moral of these two fits on $\xi$ is very different, and neither is perfect.

In the first fit, we assume all the observed 4L events came from the $g^{\mu\nu}$ tensor structure in Eq. (16), which is similar in spirit to the conventional approach of using the signal strength to constrain $\xi$ without including ${\cal O}(p^{4})$ effects. We have seen in Section II that this holds only at leading order in the derivative expansion. One could justify this approach somewhat by pointing out that the central values of the anomalous HVV couplings from the shape information are all extremely close to zero Sirunyan et al. (2019b). However the uncertainties remain significant, as can be seen in Table 2, and it is not clear how to interpret $\mu_{F}^{\text{4L}}$ when the anomalous couplings are turned on.

In the second fit, we attempted to subtract out the 4L events originating from the anomalous HVV coupling, by using the experimental limit on ${\cal R}_{i}^{V}$ in Table 2. The remaining events then can be interpreted as arising entirely from the $g^{\mu\nu}$ tensor structure, whose coefficient is given by $2\sqrt{1-\xi}$ in Eq. (II.2). However, it is not clear what “signal strength” one should use to subtract out the anomalous 4L events, as $\mu_{F}^{\text{4L}}$ in Table 1 is extracted assuming the anomalous coupling vanishes.

Having made these qualifications on our fitting procedures, there are interesting features in Fig. 3 which are likely to survive even after a more rigorous fitting is adopted:

• •

The most prominent feature in the right panel of Fig. 3 is that $\xi$ is preferred to be negative, due to the fact that $\mu_{F}^{\text{4L}}$ in Table 1 is larger than 1 in all three different likelihood fits.

• •

In the left panel, the uncertainty in $\xi$ is still rather large, allowing for $\xi\lesssim 0.5$ when ${\cal R}_{4}^{Z\gamma}$ is turned on or $\xi\gtrsim-0.5$ when ${\cal R}_{4}^{ZZ}$ is allowed. It turns out that the reason behind the rather loose limits on $\xi$ is different between these two scenarios.

In the case of ${\cal R}_{4}^{Z\gamma}$, the loose constraint is due to a combination to two factors: the rather large experimental uncertainty in Table 2 and a large numerical coefficient, from phase space integration, in front of the $|{\cal R}_{4}^{Z\gamma}|^{2}$ term in Eq. (26). The large experimental uncertainty implies the differential spectra from $A_{4}^{Z\gamma}$ are quite similar to those from the leading order $A_{1}^{ZZ}$. As a consequence, the likelihood fit is unable to separate the two tensor structures.

On the other hand, while the limit on ${\cal R}_{4}^{ZZ}$ is the strongest in Table 2, implying the likelihood fit is capable of distinguishing this particular tensor structure efficiently, any small presence of $A_{4}^{ZZ}$ is amplified by the large numerical coefficient from phase space integration defined in Eq. (19). This can be seen explicitly either in Eq. (26), where the interference term linear in ${\cal R}_{4}^{ZZ}$ has a numerical coefficient as large as the leading order coefficient.

It is possible to further obtain constraints on the Wilson coefficients defined in Eq. (9), by using the mapping in Eq. (II.2). Again we can perform two different fits using the rate and the shape information, respectively.

For rate measurements, we re-write Eq. (26) in terms of $\xi,g_{\rho}$ and $C_{i}^{h}$,

[TABLE]

where we have plugged in $M_{\rho}=g_{\rho}f$. With this we can use the rate measurements in 4L channel to examine bounds in the $(C_{i},\xi)$ and $(g_{\rho},\xi)$ two dimensional planes, by fixing the third variable. In the top row of Fig. 4, we fix $g_{\rho}=1.5$ and examine the $(C_{i},\xi)$ plane. (Recall that in QCD $\alpha_{s}(m_{b})\approx 0.22$, which corresponds to $g_{s}\approx 1.7$.) We see that turning on $C_{i}^{h}$ could have a non-negligible impact on the extraction of $\xi$. However, the impact gets diminished as $g_{\rho}$ becomes larger and larger because of the $1/g_{\rho}^{2}$ and $1/g_{\rho}^{4}$ dependence in Eq. (27). This feature is demonstrated explicitly in the second row of Fig. 4 where we fix $C_{i}^{h}=1$ and plot the constraints in the $(g_{\rho},\xi)$ plane.

As for the shape measurements, since the CMS constraints are presented in terms of the percentage of 4L events originating from the anomalous HVV couplings Khachatryan et al. (2015); Sirunyan et al. (2019b), we find it convenient to obtain bounds on the ratios $C_{i}^{h}/A_{1}^{ZZ}$, which are independent of the 4L signal strength. The outcome is presented in  Fig. 5. We see that current bounds are still week and, depending on which coupling is varied in the fit, still allow for large positive or negative values of the $C_{i}^{h}$.

III.1.2 Projections with Multi-dimensional Parameter Likelihood Fits

Currently at the LHC only one parameter fits are performed in the 4L channel Khachatryan et al. (2015); Sirunyan et al. (2019b). In the future this channel offers a golden opportunity to conduct multi-parameter fits, as demonstrated in the matrix element method (MEM) framework developed in Refs. Chen et al. (2013); Chen and Vega-Morales (2014); Chen et al. (2015a, 2014a, c). In this framework all decay observables are utilized and a combined likelihood for the $2e2\mu,4e,4\mu$ final states is constructed from the normalized fully differential decay width. The dominant $q\bar{q}\to$ 4L background, computed analytically in Refs. Gainer et al. (2011); Chen et al. (2013); Chen and Vega-Morales (2014), is also included in the likelihood. This likelihood function is a function of all anomalous HVV couplings in Eq. (16), and allowing all them to vary simultaneously we can obtain projection curves for future sensitivity to the Wilson coefficients in the nonlinear Higgs Lagrangian by the mapping in Eq. (II.2). Details of the statistical analysis and likelihood maximization procedure can be found in Chen et al. (2013); Chen and Vega-Morales (2014); Chen et al. (2015a, 2014a, c, 2016).

In the left panel of Fig. 6 we show projections of the $95\%$ C.L. contours for $C_{i}^{h}/A_{1}^{ZZ}$ from the shape measurement as a function of luminosity $(\mathcal{L})$, or number of signal events $(N_{S})$, at the LHC, assuming $M_{\rho}=1$ TeV. When calculating the necessary luminosity we include the ggH+VBF production channels at $\sqrt{S}=14$ TeV. The solid lines are projections for an 8-dimensional parameter fit, allowing all anomalous HVV couplings in Eq. (16) to vary simultaneously while the dashed lines are obtained allowing only a single coupling to vary, as done in current CMS analyses. In the projection we have included $C_{4}^{h}$, although it is not currently included in the CMS shape measurements in the 4L channel using Run 2 data.

On the right-panel we present projections for $M_{\rho}$ from the shape measurements using $C_{i}^{h}/A_{1}^{ZZ}=1/2$. We can see that, in both plots, at large statistics single parameter and multi-parameter fits converge to similar values. However, at low statistics we see that single parameter fits could lead to an overly optimistic sensitivity to the Wilson coefficients and in particular for $C_{4}^{h}$. Eventually at the HL LHC Cepeda et al. (2019), the sensitivity to each coupling could improve by around an order of magnitude from current limits. As can also be seen, the strongest sensitivity corresponds to $C_{4}^{h}$ where values $|C_{4}^{h}/A_{1}^{ZZ}|\sim 4$ can eventually be probed or correspondingly, scales $M_{\rho}\sim 750$ GeV. This stronger sensitivity is due to the fact the differential spectra corresponding to the $A_{2}^{Z\gamma}$ can be efficiently distinguished from the dominant tree level $A_{1}^{ZZ}$ operator Chen et al. (2014a). This is primarily due to the photon propagator which introduces a pole at low di-lepton invariant mass. However, in the case of the $A_{4}^{Z\gamma}$ (corresponding to $C_{3}^{h}$), the $k_{1,2}^{2}$ dependence in the tensor structure (see Eq. (16)) cancels this pole thus eliminating the distinguishing feature in the differential spectra.

III.2 Direct H$\,\to\,$Z$\gamma$ Decay

The Wilson coefficient $C_{4}^{h}$ contributes to on-shell H $\rightarrow$ Z$\gamma$ decays, which can be parametrized as Giudice et al. (2007):

[TABLE]

where the $\mathcal{O}(1)$ coefficient in front of $C_{4}^{h}$ is due to the fact that the SM H $\rightarrow$ Z$\gamma$ process arises at the one-loop level. Although present measurements from the LHC on the signal strength $\mu_{Z\gamma}$ only gives $6.6$ at $95\%$ CL upper limit Aaboud et al. (2017); Sirunyan et al. (2018), it already puts strong constraint on the coefficient $C_{4}^{h}$:

[TABLE]

with $C_{4}^{h}>0(C_{4}^{h}<0)$ respectively. The HL-LHC will be able to measure this channel with $20\%$ uncertainty Dawson et al. (2013), which will in turn give the constraint on $C_{4}^{h}$:

[TABLE]

In Fig. 7, we show the $95\%$ C.L. constraint on $C_{4}^{h}$ with $M_{\rho}=$ 1 TeV from the present measurement with integrated luminosity $L=36.1$ fb*-1* Aaboud et al. (2017) and from the HL-LHC prospective Dawson et al. (2013). Note also that when combining the $\mu_{Z\gamma}$ signal strength with the $H\to V^{\ast}V$ signal strength, their ratio is directly sensitive to the nonlinearity of a composite Higgs boson Cao et al. (2019).

III.3 Precision Electroweak Constraints

While we have demonstrated that precision Higgs measurements, in particular the golden H $\to$ 4L channel, can directly constrain nonlinear Higgs dynamics, it is well-known that there are other low-energy, indirect constraints on anomalous HVV couplings from precision electroweak test (EWPT), which we consider in this section. The $S,T$ parameters are related to $\xi$ as follows Barbieri et al. (2007):

[TABLE]

where $\Lambda$ is the UV cutoff that regularizes the logarithmic divergence. In the following, we will choose $\Lambda=4\pi f$. The fact that the IR contributions to the $S,T$ parameters have different sign put a strong constraint on the value of $\xi$, as the EWPT gives strong positive correlation $92\%$ among the $S$ and $T$ parameter Tanabashi et al. (2018). The present fit to the EWPT gives Tanabashi et al. (2018):

[TABLE]

which will be used in the following analysis.

In addition to the IR contribution, there are potential UV contributions arising from the presence of $\mathcal{O}(p^{4})$ operators in the effective Lagrangian. In particular, $O_{5}^{+}$ will contribute to the $S$ parameter at the tree-level.

[TABLE]

There are also contributions from the one-loop threshold corrections involving the vector or fermionic resonances, which are model-dependent Contino and Salvarezza (2015a, b); Ghosh et al. (2017). Instead of going into the detail of a particular UV construction, we simply choose some benchmark values for the these contributions. Note that the one loop effects can give negative and positive contributions to the $S$ and $T$ parameters, respectively, which are in the opposite direction from the IR contribution in Eq. (31). This observation can potentially relax the bound on $\xi$ from EWPT.

In addition, since $\mathcal{O}_{5}^{+}$ contribution to the $S$ parameter arises at the tree-level, as can been seen from the large coefficient in Eq. (33), only a small negative $c_{5}^{+}$ is needed to relax the bound. They can be achieved by noticing that there are potential cancellations between different contributions to $c_{5}^{+}$ in the Lagrangians of the composite spin-1 resonances Contino et al. (2011):

[TABLE]

where $\alpha_{2}$ is the coefficient of the operator $Q_{2}$ defined in Ref. Contino et al. (2011). There are also potential positive contribution to the $T$ parameter coming from the loop of the fermionic resonances, for example, the electroweak singlet top partner Lavoura and Silva (1993); Panico et al. (2011). We are not going to discuss in detail this possibility, but only take a benchmark value of $T_{f}=0.08$ as an illustration. In the left panel of Fig. 8, we show the $\Delta\chi^{2}$ as function of $\xi$ under different three different assumptions:

${c_{5}^{+}}/{g_{\rho}^{2}}=0,\ T_{f}=0$  , 2. 2.

${c_{5}^{+}}/{g_{\rho}^{2}}=-{1}/{16^{2}},\ T_{f}=0$ , 3. 3.

${c_{5}^{+}}/{g_{\rho}^{2}}=-{1}/{16^{2}},\ T_{f}=0.08$.

We can see that the constraint on $\xi$ is very strong in the absence of other contributions: $\xi\in[-0.11,0.002]$ at $95\%$ C.L., which again prefers a negative value of $\xi$. A small negative value of $c_{5}^{+}$ relax the bound on the negative $\xi$ to $-0.19$ and additional positive contribution from $T_{f}$ can relax the bound on positive value of $\xi\in[-0.015,0.16]$. In the right panel we present the $95\%$ C.L. allowed region in the $\xi-c_{5}^{+}$ plane with $g_{\rho}=1.5$. We can see that the bound on $\xi$ is very strong with only IR contribution and can be significantly relaxed in the presence of small negative value of $c_{5}^{+}$ and positive contribution $T_{f}$, which is consistent with the results discussed above.

We emphasize that the precision electroweak constraints are orthogonal to direct measurements of Higgs couplings, as they involve different sets of assumptions. Therefore it is important to pursue both of them independently.

IV Conclusions

Nonlinear dynamics of the 125 GeV Higgs boson is the most salient feature of a composite Higgs boson. The nonlinear interaction realize the PNGB nature of the Higgs, in the same way pions in low-energy QCD manifest their PNGB nature through the nonlinear interactions. In fact, the nonlinear dynamics can serve as the defining property of a composite Higgs boson, independent of how the fermionic sector is implemented, even when the so-called top partner is neutral under all SM charges. In particular, recent theoretical advances showed the nonlinear interaction is an IR property of the composite Higgs boson that is insensitive to the symmetry-breaking pattern invoked in the UV. At the two-derivative level, the nonlinear Lagrangian is determined by a single parameter $f$, the Goldstone decay constant. In the low-energy regime, $f$ needs to be taken as an input from the experimental data. Traditionally this is done by relating $f$ to the signal strength of HVV couplings assuming the observed events arise entirely from ${\cal O}(p^{2})$ operators in the nonlinear Higgs Lagrangian. This is obviously an over-simplification as ${\cal O}(p^{4})$ operators introduce several new anomalous tensor structures to HVV couplings.

In order to disentangle effects of ${\cal O}(p^{2})$ operators from those of ${\cal O}(p^{4})$ operators, it is crucial to include shape information, as different tensor structures lead to different shapes in the fully differential spectra. In this regard, the H$\to$ 4L “golden channel” decay is the ideal probe of these nonlinear Higgs dynamics. In this work we have performed, for the first time, experimental fits to the Wilson coefficients of ${\cal O}(p^{4})$ operators. In addition, we demonstrated the limitation of using only the signal strength to constrain the nonlinear parameter $\xi$ at leading order in derivative expansions. We showed that it is important to include both the rate information in the signal strength measurements and the shape information in the fully differential spectra, in order to constrain the nonlinear parameter $\xi$ and the Wilson coefficients $C_{i}^{h}$. The fitting procedures we adopted are less than ideal and limited by the methods employed in current LHC analyses.

We found that in rate measurements $\xi$ is preferred to be negative, pointing to a non-compact coset structure in the UV Low (2015a, b). Moreover, using the shape measurements we identified scenarios where $\xi$ could be as large as +0.5 or -0.5, corresponding to turning on $A_{4}^{Z\gamma}$ and $A_{4}^{ZZ}$, respectively. Such a large $\xi$ could be further constrained in a specific UV model, either by precision electroweak measurements or direct searches in the fermionic sector. It would be an interesting model-building challenge to devise a concrete UV model realizing such a large $\xi$ while avoiding other experimental constraints.

Our work points to a new experimental frontier in using precision Higgs measurements to probe nonlinear dynamics of a composite Higgs boson, beyond the conventional signal strength measurements. Given that the Higgs boson remains a top priority in current and future experimental programs in high-energy colliders, it is desirable to introduce new tools and techniques to further the experimental analysis in measurements of Higgs properties. We leave this direction for future works.

Acknowledgements.

We are grateful to Yi Chen for his help in obtaining some of the results presented in this work and Daniel Stolarski for useful comments. This work is supported in part by the U.S. Department of Energy under contracts No. DE-AC02-06CH11357 and No. DE-SC0010143. The work of R.V.M. is supported by MINECO grants FPA 2016-78220-C3-1-P, FPA 2013-47836-C3-2/3-P (including ERDF), and the Juan de la Cierva program, as well as by the Junta de Andalucia Project FQM-101.