Higgs Effective Field Theory (HEFT)

• Higgs Effective Field Theory (HEFT) is a general low-energy framework that models non-linear interactions of the Higgs boson, gauge bosons, and fermions.
• It systematically organizes chiral expansions using derivative power counting and accommodates both decoupling and non-decoupling dynamics.
• HEFT underpins precise phenomenology, from vector boson scattering to resonance searches, with significant implications for UV model matching and electroweak constraints.

Higgs Effective Field Theory (HEFT) is the most general low-energy effective field theory capturing the interactions of the Standard Model (SM) gauge bosons, fermions, and a light scalar singlet (the observed Higgs boson), under a non-linear realization of electroweak symmetry breaking. Unlike the Standard Model Effective Field Theory (SMEFT), which assumes a linear $SU(2)_L$ doublet structure for the scalar sector, HEFT allows arbitrary couplings, including non-decoupling and strongly coupled scenarios, by treating the Higgs as an $SU(2)_L \times U(1)_Y$ singlet and parameterizing Goldstone bosons independently. This framework enables a systematic expansion in chiral (derivative) order and powers of $h/v$, and provides a robust platform for analyzing phenomenology ranging from vector boson scattering, Higgs pair production, electroweak precision measurements, to nonperturbative soliton states.

1. Structural Foundations and Lagrangian Formulation

In HEFT, the scalar sector after electroweak symmetry breaking is described by a set of three Goldstone bosons (denoted $\omega^i$ or $\pi^a$) embedded nonlinearly via a unitary matrix $U(x) = \exp(i \omega^i(x) \sigma^i / v)$ and an independent real singlet scalar $h$ corresponding to the observed 125 GeV Higgs boson. The most general custodial-symmetric bosonic Lagrangian to next-to-leading chiral order takes the form (Dobado et al., 2019, Alonso et al., 2015, Brivio et al., 2017):

$$\begin{aligned} \mathcal{L}_{\rm HEFT} &= \frac{1}{2}\Big[1 + 2a \frac{h}{v} + b \Big(\frac{h}{v}\Big)^2\Big] \partial_\mu\omega^i\,\partial^\mu\omega^j\,\Big[\delta_{ij} + \frac{\omega^i\omega^j}{v^2}\Big] \ & + \frac{1}{2} \partial_\mu h \, \partial^\mu h - V(h) \ &+ \sum_i \alpha_i \mathcal{O}_i\,[{\rm dim}=4] + \cdots \end{aligned}$$

with $v\simeq246$ GeV the electroweak vacuum expectation value, $a$, $b$ encoding $hVV$ and $hhVV$ anomalous couplings, and $\alpha_i$ Wilson coefficients for all independent higher-order chiral operators $\mathcal{O}_i$ built from $\{U,h\}$, gauge fields, and covariant derivatives. The minimal SM is recovered for $a=b=1$ and all higher $\alpha_i=0$.

HEFT power-counting employs the chiral dimension, counting $\partial_\mu$ and masses as $O(p)$, and organizes operators accordingly. The theory does not require an underlying linear multiplet for the Higgs, so it can encompass the full range of possible scalar manifold geometries, including cases with curved or flat field spaces (Alonso et al., 2015, Nagai et al., 2019).

2. Geometry, Nonlinear Symmetry, and Physical Significance

A central insight in modern HEFT formulations is the geometric interpretation of the scalar field space, $\mathcal{M}$. The general HEFT kinetic term is a nonlinear sigma model on $\mathcal{M}$ with an $SU(2)_L \times U(1)_Y$-invariant metric $g_{ij}(\phi)$, where $\phi^i$ denotes coordinates spanning both Goldstone and Higgs directions (Alonso et al., 2015, Nagai et al., 2019):

$$\mathcal{L}_{\text{scalar}} = \frac{1}{2} g_{ij}(\phi) D_\mu \phi^i D^\mu \phi^j$$

Deviations from the SM are characterized by nonzero Riemann curvature tensors $R^i_{\,\,jkl}$ on $\mathcal{M}$. Important physical observables—the $hVV$ coupling, double-Higgs production amplitudes, and the electroweak $S$ parameter—directly probe these geometrical invariants:

• The $hVV$ and $hhVV$ couplings are tied to the expansion coefficients of the radial form-factor $F(h)$, with $1-F'(0)^2v^2$ proportional to the curvature in the Goldstone sector.
• The growth of high-energy vector boson scattering amplitudes (e.g., $W_LW_L \to W_LW_L$) is fixed by components of $R_{ijkl}$, with unitarity violation scale $\Lambda \sim 4\pi v/\sqrt{r_i}$.
• The one-loop contribution to the $S$ parameter is proportional to the curvature $r_4$, controlling log-enhanced high-scale sensitivity (Alonso et al., 2015, Nagai et al., 2019).

A requirement of tree-level unitarity in longitudinal scattering fixes $\mathcal{M}$ to be flat, enforcing relations among couplings (the so-called unitarity sum rules). Flatness in this context guarantees the absence of UV divergences in the one-loop $S$ and $U$ parameters, closing the link between symmetry, geometry, and renormalizability (Nagai et al., 2019). These geometric considerations underpin the theoretical consistency of UV completions and EFT convergence.

3. Operator Bases and Classification: On-Shell and Hilbert Series Techniques

The systematic construction and enumeration of independent HEFT operators have advanced through both amplitude-based and algebraic techniques:

• On-shell amplitude construction: Utilizing the massive spinor-helicity and Young-Tableau methods, complete $n$-point operator bases up to any operator dimension are systematically built, with each operator mapped directly to a physical S-matrix element (Dong et al., 2022). This approach streamlines the identification of redundancies (IBP/EOM) and provides physical meaning through direct connection with scattering amplitudes.
• Hilbert series enumeration: The operator basis can be counted and constructed using Hilbert series techniques extended to massive, nonlinearly realized settings (Gráf et al., 2022). This method incorporates the Goldstone equivalence theorem at the partition-function level and allows for consistent inclusion of spurion fields, automatically capturing custodial symmetry breaking.

At NNLO, the HEFT operator basis grows rapidly, with $>10^4$ structures for a single fermion family (Sun et al., 2022). Matching at dimension-8 and higher between SMEFT and HEFT involves nontrivial field redefinitions and operator projection, with the NNLO HEFT basis providing the target for SMEFT dimension-8 amplitude decompositions.

4. Matching to UV Models and Interpretation of Wilson Coefficients

Matching HEFT to specific ultraviolet scenarios (such as singlet scalar extensions, real Higgs triplet models, or the Two-Higgs-Doublet Model) reveals both decoupling and nondecoupling phenomena (Arco et al., 2023, Buchalla et al., 2023, Song et al., 2 Mar 2025, Dawson et al., 2023):

• In the decoupling regime, HEFT Wilson coefficients scale as powers of $v^2/\Lambda^2$, and the SMEFT expansion is reliable.
• In nondecoupling or strongly coupled scenarios (e.g., large scalar quartics), leading deviations in $a,b,\kappa_3,\kappa_4$ and loop-induced $h\to\gamma\gamma$ coefficients can remain even as new scalar masses are sent to infinity.
• The choice of expansion scheme or "power-counting" (scaling assumptions for UV parameters) is not unique, and different parameter scalings yield inequivalent truncations of the low-energy HEFT. Consequently, the interpretation of HEFT measurements in terms of UV parameters carries a systematic ambiguity (Dawson et al., 2023).

Functional matching algorithms, both at tree and one loop, permit the resummation of all powers of $h/v$ in the coefficient functions for the Goldstone kinetic term, scalar potential, and Yukawa interactions (Buchalla et al., 2023). Explicit formulas for the dependence of HEFT coefficients on the UV parameters are tabulated in the referenced works.

5. Phenomenological Applications: Vector Boson Scattering, Resonances, and Beyond

HEFT provides a robust, unitarized framework for computing both resonant and nonresonant $VV$ and $Vh$ cross sections at LHC energies, notably in high-energy vector boson scattering (VBS) (Dobado et al., 2019, Asiáin et al., 2021). The workflow involves:

• Computing elastic longitudinal $VV$ scattering amplitudes (using the Goldstone Equivalence Theorem), expanded at one loop with renormalized partial waves.
• Employing dispersion relations and solving unitarity constraints via either the Inverse Amplitude Method (IAM) or the $N/D$ technique to yield analytic, unitary amplitudes.
• The unitarized amplitudes can support dynamical resonances, manifest as poles in the second Riemann sheet. Their masses and widths, $M_R$ and $\Gamma_R$, are directly calculable as functions of HEFT parameters ($a,b,a_4,a_5,...$).

Phenomenologically, HEFT allows for realistic predictions of VBS and $Wh$ invariant mass distributions at the LHC, with the position and existence of resonances highly sensitive to parameter combinations such as $(e-2d)$. As this approaches zero, resonances decouple to infinity (The SM limit); higher values enhance new-physics signals in LHC processes (Dobado et al., 2019). The framework can be extended to multiple Higgs production, Higgs-plus-jet final states, non-perturbative objects such as electroweak skyrmions (Criado et al., 2021), and global fits to EW precision, single, and double Higgs data.

6. UV Resonance Saturation, Low-Energy Constants, and Experimental Constraints

Matching to explicit resonance models (spin-0 and spin-1 states) at high energies, integrating out the heavy multiplets under short-distance constraints, saturates the bosonic HEFT low-energy constants (LECs): $F_1, F_3, F_4, F_5$ (Rosell et al., 2021, Rosell et al., 2020). These are analytically expressible as functions of resonance couplings and masses:

$$\begin{aligned} F_1 &= -\frac{F_V^2}{4M_V^2} + \frac{F_A^2}{4M_A^2}, \ F_3 &= -\frac{F_V G_V}{2M_V^2}, \ F_4 &= \frac{G_V^2}{4M_V^2} - \frac{G_A^2}{4M_A^2}, \ F_5 &= -\frac{G_V^2}{4M_V^2} + \frac{c_d^2}{4M_S^2} \end{aligned}$$

Imposing vanishing of high-energy form-factors and Weinberg sum rules (stemming from asymptotic UV behavior), the parameter space collapses to functions of resonance masses. Current experimental constraints from oblique parameters ($S$), triple and quartic gauge couplings, and LHC direct measurements require these resonance masses to lie generically above $2$ TeV (Rosell et al., 2021, Rosell et al., 2020). This illustrates the power of precision measurements to directly probe or exclude patterns of strong coupling in the EWSB sector.

7. Theoretical Flexibility and Phenomenological Outlook

HEFT subsumes both weakly coupled, decoupling BSM scenarios (retrieving SMEFT in the appropriate limit) and strongly coupled, composite, or non-decoupling paradigms. It allows arbitrary functional dependence in $h/v$ and organizing independent multi-Higgs vertices at each order in chiral and $h/v$ expansion. The operator basis, constructed via modern amplitude or algebraic approaches, enables thorough phenomenological studies at single-, double-, and multi-Higgs level, as well as in high-energy VBS observables.

Careful matching to UV models must specify scaling/power-counting assumptions, as ambiguities in parameter assignment propagate directly into the inferred patterns of HEFT coefficients (Dawson et al., 2023). While the flexibility of HEFT is phenomenologically essential, it necessitates comprehensive and systematic treatments in both theory and global data analyses.

HEFT thus stands as the modern paradigm for model-independent analyses of Higgs and electroweak data, capable of encompassing the full landscape of symmetry-breaking dynamics in the post-Higgs discovery era (Alonso et al., 2015, Dobado et al., 2019, Rosell et al., 2021, Brivio et al., 2017).