Papers
Topics
Authors
Recent
Search
2000 character limit reached

General Two Higgs Doublet Model (G2HDM)

Updated 7 July 2026
  • The General Two Higgs Doublet Model (G2HDM) is the most general renormalizable extension of the Standard Model with two complex SU(2) doublets, featuring a complete scalar potential and flexible Yukawa sector.
  • It incorporates both tree-level flavor-changing neutral currents and multiple CP violation sources, enabling rich phenomenology from collider signals to electroweak baryogenesis.
  • The model’s basis independence, unitarity constraints, and effective field theory formulations offer a comprehensive framework for exploring UV completions and experimental tests.

Searching arXiv for the cited G2HDM literature and related terminology. The General Two Higgs Doublet Model (G2HDM) is the most general renormalizable extension of the Standard Model scalar sector based on two complex SU(2)LSU(2)_L doublets, without imposing a discrete Z2\mathbb{Z}_2 symmetry and without requiring CP conservation. In that sense it is the unconstrained 2HDM: the scalar potential contains the full set of renormalizable operators, the Yukawa sector contains two independent Yukawa structures for each fermion type, and tree-level flavor-changing neutral scalar interactions and new CP-violating phases are generic unless controlled by alignment, textures, or symmetries. After electroweak symmetry breaking, the spectrum contains five physical Higgs states, conventionally denoted H±H^\pm, AA, hh, and HH, with basis freedom playing a central role in how parameters are represented (Bahl et al., 2022, Sartore et al., 2022, Branco et al., 2011). A recurrent source of ambiguity is nomenclature: in the arXiv literature the acronym G2HDM is also used for the Gauged Two Higgs Doublet Model, which is a distinct framework with an extra SU(2)H×U(1)XSU(2)_H\times U(1)_X gauge sector (Huang et al., 2015, Yang et al., 2024).

1. Scalar-sector definition and basis structure

The scalar potential of the general 2HDM is conventionally written as

V=m112Φ1Φ1+m222Φ2Φ2(m122Φ1Φ2+h.c.) +λ12(Φ1Φ1)2+λ22(Φ2Φ2)2+λ3(Φ1Φ1)(Φ2Φ2)+λ4(Φ1Φ2)(Φ2Φ1) +[λ52(Φ1Φ2)2+λ6(Φ1Φ1)(Φ1Φ2)+λ7(Φ2Φ2)(Φ1Φ2)+h.c.].\begin{aligned} V &= m_{11}^2 \Phi_1^\dagger \Phi_1 + m_{22}^2 \Phi_2^\dagger \Phi_2 - (m_{12}^2 \Phi_1^\dagger \Phi_2 + h.c.) \ &\quad + \frac{\lambda_1}{2} (\Phi_1^\dagger \Phi_1)^2 + \frac{\lambda_2}{2} (\Phi_2^\dagger \Phi_2)^2 + \lambda_3 (\Phi_1^\dagger \Phi_1)(\Phi_2^\dagger \Phi_2) + \lambda_4 (\Phi_1^\dagger \Phi_2)(\Phi_2^\dagger \Phi_1) \ &\quad + \left[\frac{\lambda_5}{2} (\Phi_1^\dagger \Phi_2)^2 + \lambda_6 (\Phi_1^\dagger \Phi_1)(\Phi_1^\dagger \Phi_2) + \lambda_7 (\Phi_2^\dagger \Phi_2)(\Phi_1^\dagger \Phi_2) + h.c.\right]. \end{aligned}

Here m112m_{11}^2, m222m_{22}^2, and Z2\mathbb{Z}_20 are real, while Z2\mathbb{Z}_21, Z2\mathbb{Z}_22, Z2\mathbb{Z}_23, and Z2\mathbb{Z}_24 may be complex, giving 14 real parameters in the scalar potential before exploiting basis freedom (Kanemura et al., 2015, Bahl et al., 2022).

A central structural fact is that Z2\mathbb{Z}_25 and Z2\mathbb{Z}_26 are not physical observables. They may be redefined by Z2\mathbb{Z}_27 transformations in Higgs-family space, and many apparently distinct parameter sets therefore describe the same physics. The Higgs basis is particularly useful: only one linear combination acquires a vacuum expectation value, while the orthogonal combination contains the non-SM scalar degrees of freedom. In one standard parametrization,

Z2\mathbb{Z}_28

with Z2\mathbb{Z}_29 and H±H^\pm0 (Kanemura et al., 2015). In the alignment limit, the observed 125 GeV state is dominantly aligned with the vacuum direction, and its couplings approach the SM pattern.

The gauge-invariant bilinear formulation sharpens this basis structure. Defining the hermitian matrix

H±H^\pm1

one obtains the basis-independent constraints

H±H^\pm2

In this language, stability, electroweak symmetry breaking, CP properties, scalar masses, trilinear and quartic couplings, gauge interactions, and Yukawa couplings can all be reformulated in terms of gauge-invariant quantities. An especially compact result is the all-order relation

H±H^\pm3

for the charged Higgs mass at a charge-conserving minimum (Sartore et al., 2022).

2. Yukawa structure, flavor violation, and restricted limits

The defining phenomenological feature of the G2HDM is its Yukawa generality. Unlike type-I, type-II, type-X, or type-Y models, both scalar doublets may couple to all fermions. In the Higgs-basis language used in flavor analyses, this introduces, in addition to the SM-like Yukawas H±H^\pm4, a second set of Yukawa matrices H±H^\pm5 for each fermion sector H±H^\pm6. In the alignment regime, the observed Higgs H±H^\pm7 remains approximately SM-like, while the extra neutral and charged scalars H±H^\pm8, H±H^\pm9, and AA0 retain unsuppressed AA1-driven interactions (Hou et al., 2023, Teunissen, 2023).

Tree-level flavor-changing neutral currents are then generic, because the two Yukawa structures for a given fermion type are not generically simultaneously diagonalizable. This is the key difference between the general model and the natural-flavor-conserving constructions enforced by discrete symmetries. The standard restricted models are recovered as symmetry limits: type I, type II, lepton-specific, and flipped models impose specific Yukawa assignments so that all fermions of a given charge and chirality couple to only one doublet, thereby eliminating tree-level FCNCs (Branco et al., 2011). In a broader flavor taxonomy, the fully general setup has also been denoted 2HDM-X or, in flavor-centered discussions, associated with type-III 2HDM language (Diaz-Cruz et al., 2010, Herrero-Garcia et al., 2019).

Several control mechanisms have been developed. Yukawa alignment assumes that the two Yukawa matrices for each fermion type are proportional at a high scale, so that tree-level FCNCs vanish there. However, one-loop running spoils exact alignment. In leading-log form, the induced misalignment generates off-diagonal neutral-Higgs couplings proportional to structures such as AA2, yielding a calculable lower bound on exotic FCNCs in generic aligned 2HDMs barring cancellations and exact discrete symmetries (Braeuninger et al., 2010). Other suppression strategies discussed in the literature include Cheng–Sher-type four-texture ansätze and approximate family symmetries (Diaz-Cruz et al., 2010, Branco et al., 2011).

In recent g2HDM phenomenology, a heuristic flavor pattern has often been adopted in which the extra Yukawas echo the SM hierarchy,

AA3

This permits AA4, allows AA5, and suppresses first-generation couplings. The intent is not to impose natural flavor conservation, but to encode a hierarchical flavor organization that can simultaneously soften FCNC and EDM constraints (Hou et al., 2023, Teunissen, 2023).

3. CP violation and electroweak baryogenesis

The G2HDM admits several qualitatively different sources of CP violation. CP may be violated explicitly or spontaneously in the scalar potential, or generated in the Yukawa sector through complex extra Yukawa matrices. In the language of the Higgs basis, the CP-odd invariants

AA6

vanish if and only if CP is conserved. A notable feature of the full general model is that exact alignment can coexist with CP violation; this does not occur in the softly broken AA7-symmetric case, where exact alignment forces CP conservation in the scalar sector (Bahl et al., 2022).

This enlarged CP structure is one reason the G2HDM has been studied as a framework for electroweak baryogenesis. In one widely studied realization, the Higgs potential is taken CP-conserving and all new CP violation is placed in the Yukawa sector. The baryogenesis source is then dominated by the top sector, especially

AA8

while a strong first-order electroweak phase transition is supplied by AA9 Higgs quartics. The favored exotic-scalar spectrum is sub-TeV, with benchmarks such as

hh0

and the observation that masses around hh1 GeV would likely make baryogenesis even more efficient (Hou et al., 2023).

Electric dipole moments provide the sharpest low-energy probe of these CP phases. In the g2HDM, the dominant electron EDM contribution is typically the two-loop Barr–Zee diagram, top dominated and sensitive to hh2 and hh3. A central result is the flavor-structured cancellation condition

hh4

with hh5. This “natural flavor cancellation” can keep the eEDM below present bounds while retaining baryogenesis-scale CP violation. The corresponding phenomenological expectation is that the eEDM may first emerge around hh6, possibly even around hh7, while the neutron EDM may appear around hh8, with hh9-class experiments probing even cancellation regions (Hou et al., 2023, Teunissen, 2023).

Spontaneous CP violation also arises in UV completions that reduce to the G2HDM at low energies. In the two Higgs bi-doublet left-right symmetric model, the decoupling regime HH0 yields a general 2HDM with spontaneous CP violation, nontrivial vacuum phases, three light neutral Higgs bosons, and one light charged Higgs pair at the electroweak scale (Liu et al., 2012). This illustrates that the general 2HDM CP structure is not only a bottom-up parametrization but also a natural low-energy limit of more complete theories.

4. Theoretical consistency: unitarity, boundedness from below, and vacuum structure

The G2HDM scalar sector is strongly constrained by perturbative unitarity. At high energy, the relevant amplitudes are two-body scalar scatterings, equivalently including longitudinal gauge bosons through the equivalence theorem. In the most general model, the full problem can be organized by hypercharge and isospin into block-diagonal scattering matrices. One explicit construction yields 14 neutral channels, 8 singly charged channels, and 3 doubly charged channels, reducing the unitarity problem to the eigenvalues of the blocks HH1, HH2, HH3, and the invariant HH4, with perturbative unitarity imposed on all eigenvalues (Kanemura et al., 2015).

Recent work has supplemented the exact numerical treatment with analytic necessary and sufficient criteria tailored to the fully complex potential. New tools include principal minors, Gershgorin disks, and Frobenius-norm bounds. One practical conclusion is that the crude estimate HH5 is too loose: unitarity may fail for significantly smaller values, while HH6 is much safer in numerical scans (Bahl et al., 2022).

Boundedness from below is similarly nontrivial once HH7 and HH8 are nonzero and complex. In the restricted HH9 limit, the familiar necessary and sufficient tree-level conditions are

SU(2)H×U(1)XSU(2)_H\times U(1)_X0

but these become only necessary in the general case (Jurciukonis et al., 2018). Exact bounded-from-below conditions for the fully general model have been derived by rewriting the quartic potential as a positivity problem for a quartic polynomial in field directions, and practical fast-scan tests have also been obtained (Bahl et al., 2022).

Vacuum consistency is most transparently described in the bilinear formalism. Writing the potential as

SU(2)H×U(1)XSU(2)_H\times U(1)_X1

one obtains the discriminant

SU(2)H×U(1)XSU(2)_H\times U(1)_X2

with the criterion

SU(2)H×U(1)XSU(2)_H\times U(1)_X3

A physically useful relation is

SU(2)H×U(1)XSU(2)_H\times U(1)_X4

which links the neutrality-enforcing Lagrange multiplier directly to the charged Higgs mass. This renders vacuum stability tests interpretable in terms of physical masses (Bahl et al., 2022). The same bilinear machinery also clarifies the distinction between charge-conserving, charge-breaking, and CP-breaking vacua in a basis-independent way (Sartore et al., 2022).

5. Flavor, EDM, and collider phenomenology

Low-energy flavor observables probe the extra Yukawa sector with exceptional sensitivity. In kaon physics, charged-Higgs loops driven by SU(2)H×U(1)XSU(2)_H\times U(1)_X5 couplings generate contributions to SU(2)H×U(1)XSU(2)_H\times U(1)_X6, SU(2)H×U(1)XSU(2)_H\times U(1)_X7, and rare decays such as SU(2)H×U(1)XSU(2)_H\times U(1)_X8 and SU(2)H×U(1)XSU(2)_H\times U(1)_X9. Among these, V=m112Φ1Φ1+m222Φ2Φ2(m122Φ1Φ2+h.c.) +λ12(Φ1Φ1)2+λ22(Φ2Φ2)2+λ3(Φ1Φ1)(Φ2Φ2)+λ4(Φ1Φ2)(Φ2Φ1) +[λ52(Φ1Φ2)2+λ6(Φ1Φ1)(Φ1Φ2)+λ7(Φ2Φ2)(Φ1Φ2)+h.c.].\begin{aligned} V &= m_{11}^2 \Phi_1^\dagger \Phi_1 + m_{22}^2 \Phi_2^\dagger \Phi_2 - (m_{12}^2 \Phi_1^\dagger \Phi_2 + h.c.) \ &\quad + \frac{\lambda_1}{2} (\Phi_1^\dagger \Phi_1)^2 + \frac{\lambda_2}{2} (\Phi_2^\dagger \Phi_2)^2 + \lambda_3 (\Phi_1^\dagger \Phi_1)(\Phi_2^\dagger \Phi_2) + \lambda_4 (\Phi_1^\dagger \Phi_2)(\Phi_2^\dagger \Phi_1) \ &\quad + \left[\frac{\lambda_5}{2} (\Phi_1^\dagger \Phi_2)^2 + \lambda_6 (\Phi_1^\dagger \Phi_1)(\Phi_1^\dagger \Phi_2) + \lambda_7 (\Phi_2^\dagger \Phi_2)(\Phi_1^\dagger \Phi_2) + h.c.\right]. \end{aligned}0 is one of the strongest constraints on V=m112Φ1Φ1+m222Φ2Φ2(m122Φ1Φ2+h.c.) +λ12(Φ1Φ1)2+λ22(Φ2Φ2)2+λ3(Φ1Φ1)(Φ2Φ2)+λ4(Φ1Φ2)(Φ2Φ1) +[λ52(Φ1Φ2)2+λ6(Φ1Φ1)(Φ1Φ2)+λ7(Φ2Φ2)(Φ1Φ2)+h.c.].\begin{aligned} V &= m_{11}^2 \Phi_1^\dagger \Phi_1 + m_{22}^2 \Phi_2^\dagger \Phi_2 - (m_{12}^2 \Phi_1^\dagger \Phi_2 + h.c.) \ &\quad + \frac{\lambda_1}{2} (\Phi_1^\dagger \Phi_1)^2 + \frac{\lambda_2}{2} (\Phi_2^\dagger \Phi_2)^2 + \lambda_3 (\Phi_1^\dagger \Phi_1)(\Phi_2^\dagger \Phi_2) + \lambda_4 (\Phi_1^\dagger \Phi_2)(\Phi_2^\dagger \Phi_1) \ &\quad + \left[\frac{\lambda_5}{2} (\Phi_1^\dagger \Phi_2)^2 + \lambda_6 (\Phi_1^\dagger \Phi_1)(\Phi_1^\dagger \Phi_2) + \lambda_7 (\Phi_2^\dagger \Phi_2)(\Phi_1^\dagger \Phi_2) + h.c.\right]. \end{aligned}1, often stronger than the corresponding V=m112Φ1Φ1+m222Φ2Φ2(m122Φ1Φ2+h.c.) +λ12(Φ1Φ1)2+λ22(Φ2Φ2)2+λ3(Φ1Φ1)(Φ2Φ2)+λ4(Φ1Φ2)(Φ2Φ1) +[λ52(Φ1Φ2)2+λ6(Φ1Φ1)(Φ1Φ2)+λ7(Φ2Φ2)(Φ1Φ2)+h.c.].\begin{aligned} V &= m_{11}^2 \Phi_1^\dagger \Phi_1 + m_{22}^2 \Phi_2^\dagger \Phi_2 - (m_{12}^2 \Phi_1^\dagger \Phi_2 + h.c.) \ &\quad + \frac{\lambda_1}{2} (\Phi_1^\dagger \Phi_1)^2 + \frac{\lambda_2}{2} (\Phi_2^\dagger \Phi_2)^2 + \lambda_3 (\Phi_1^\dagger \Phi_1)(\Phi_2^\dagger \Phi_2) + \lambda_4 (\Phi_1^\dagger \Phi_2)(\Phi_2^\dagger \Phi_1) \ &\quad + \left[\frac{\lambda_5}{2} (\Phi_1^\dagger \Phi_2)^2 + \lambda_6 (\Phi_1^\dagger \Phi_1)(\Phi_1^\dagger \Phi_2) + \lambda_7 (\Phi_2^\dagger \Phi_2)(\Phi_1^\dagger \Phi_2) + h.c.\right]. \end{aligned}2-physics limits, especially when the charged Higgs is heavy. A distinctive result is the special sensitivity of V=m112Φ1Φ1+m222Φ2Φ2(m122Φ1Φ2+h.c.) +λ12(Φ1Φ1)2+λ22(Φ2Φ2)2+λ3(Φ1Φ1)(Φ2Φ2)+λ4(Φ1Φ2)(Φ2Φ1) +[λ52(Φ1Φ2)2+λ6(Φ1Φ1)(Φ1Φ2)+λ7(Φ2Φ2)(Φ1Φ2)+h.c.].\begin{aligned} V &= m_{11}^2 \Phi_1^\dagger \Phi_1 + m_{22}^2 \Phi_2^\dagger \Phi_2 - (m_{12}^2 \Phi_1^\dagger \Phi_2 + h.c.) \ &\quad + \frac{\lambda_1}{2} (\Phi_1^\dagger \Phi_1)^2 + \frac{\lambda_2}{2} (\Phi_2^\dagger \Phi_2)^2 + \lambda_3 (\Phi_1^\dagger \Phi_1)(\Phi_2^\dagger \Phi_2) + \lambda_4 (\Phi_1^\dagger \Phi_2)(\Phi_2^\dagger \Phi_1) \ &\quad + \left[\frac{\lambda_5}{2} (\Phi_1^\dagger \Phi_2)^2 + \lambda_6 (\Phi_1^\dagger \Phi_1)(\Phi_1^\dagger \Phi_2) + \lambda_7 (\Phi_2^\dagger \Phi_2)(\Phi_1^\dagger \Phi_2) + h.c.\right]. \end{aligned}3 to a TeV-scale charged Higgs: the V=m112Φ1Φ1+m222Φ2Φ2(m122Φ1Φ2+h.c.) +λ12(Φ1Φ1)2+λ22(Φ2Φ2)2+λ3(Φ1Φ1)(Φ2Φ2)+λ4(Φ1Φ2)(Φ2Φ1) +[λ52(Φ1Φ2)2+λ6(Φ1Φ1)(Φ1Φ2)+λ7(Φ2Φ2)(Φ1Φ2)+h.c.].\begin{aligned} V &= m_{11}^2 \Phi_1^\dagger \Phi_1 + m_{22}^2 \Phi_2^\dagger \Phi_2 - (m_{12}^2 \Phi_1^\dagger \Phi_2 + h.c.) \ &\quad + \frac{\lambda_1}{2} (\Phi_1^\dagger \Phi_1)^2 + \frac{\lambda_2}{2} (\Phi_2^\dagger \Phi_2)^2 + \lambda_3 (\Phi_1^\dagger \Phi_1)(\Phi_2^\dagger \Phi_2) + \lambda_4 (\Phi_1^\dagger \Phi_2)(\Phi_2^\dagger \Phi_1) \ &\quad + \left[\frac{\lambda_5}{2} (\Phi_1^\dagger \Phi_2)^2 + \lambda_6 (\Phi_1^\dagger \Phi_1)(\Phi_1^\dagger \Phi_2) + \lambda_7 (\Phi_2^\dagger \Phi_2)(\Phi_1^\dagger \Phi_2) + h.c.\right]. \end{aligned}4 amplitude contains a double CKM enhancement of the V=m112Φ1Φ1+m222Φ2Φ2(m122Φ1Φ2+h.c.) +λ12(Φ1Φ1)2+λ22(Φ2Φ2)2+λ3(Φ1Φ1)(Φ2Φ2)+λ4(Φ1Φ2)(Φ2Φ1) +[λ52(Φ1Φ2)2+λ6(Φ1Φ1)(Φ1Φ2)+λ7(Φ2Φ2)(Φ1Φ2)+h.c.].\begin{aligned} V &= m_{11}^2 \Phi_1^\dagger \Phi_1 + m_{22}^2 \Phi_2^\dagger \Phi_2 - (m_{12}^2 \Phi_1^\dagger \Phi_2 + h.c.) \ &\quad + \frac{\lambda_1}{2} (\Phi_1^\dagger \Phi_1)^2 + \frac{\lambda_2}{2} (\Phi_2^\dagger \Phi_2)^2 + \lambda_3 (\Phi_1^\dagger \Phi_1)(\Phi_2^\dagger \Phi_2) + \lambda_4 (\Phi_1^\dagger \Phi_2)(\Phi_2^\dagger \Phi_1) \ &\quad + \left[\frac{\lambda_5}{2} (\Phi_1^\dagger \Phi_2)^2 + \lambda_6 (\Phi_1^\dagger \Phi_1)(\Phi_1^\dagger \Phi_2) + \lambda_7 (\Phi_2^\dagger \Phi_2)(\Phi_1^\dagger \Phi_2) + h.c.\right]. \end{aligned}5 term, making this decay far more powerful than naive heavy-mass suppression would suggest (Hou et al., 2022).

The same scan studies show that V=m112Φ1Φ1+m222Φ2Φ2(m122Φ1Φ2+h.c.) +λ12(Φ1Φ1)2+λ22(Φ2Φ2)2+λ3(Φ1Φ1)(Φ2Φ2)+λ4(Φ1Φ2)(Φ2Φ1) +[λ52(Φ1Φ2)2+λ6(Φ1Φ1)(Φ1Φ2)+λ7(Φ2Φ2)(Φ1Φ2)+h.c.].\begin{aligned} V &= m_{11}^2 \Phi_1^\dagger \Phi_1 + m_{22}^2 \Phi_2^\dagger \Phi_2 - (m_{12}^2 \Phi_1^\dagger \Phi_2 + h.c.) \ &\quad + \frac{\lambda_1}{2} (\Phi_1^\dagger \Phi_1)^2 + \frac{\lambda_2}{2} (\Phi_2^\dagger \Phi_2)^2 + \lambda_3 (\Phi_1^\dagger \Phi_1)(\Phi_2^\dagger \Phi_2) + \lambda_4 (\Phi_1^\dagger \Phi_2)(\Phi_2^\dagger \Phi_1) \ &\quad + \left[\frac{\lambda_5}{2} (\Phi_1^\dagger \Phi_2)^2 + \lambda_6 (\Phi_1^\dagger \Phi_1)(\Phi_1^\dagger \Phi_2) + \lambda_7 (\Phi_2^\dagger \Phi_2)(\Phi_1^\dagger \Phi_2) + h.c.\right]. \end{aligned}6 is achievable after imposing V=m112Φ1Φ1+m222Φ2Φ2(m122Φ1Φ2+h.c.) +λ12(Φ1Φ1)2+λ22(Φ2Φ2)2+λ3(Φ1Φ1)(Φ2Φ2)+λ4(Φ1Φ2)(Φ2Φ1) +[λ52(Φ1Φ2)2+λ6(Φ1Φ1)(Φ1Φ2)+λ7(Φ2Φ2)(Φ1Φ2)+h.c.].\begin{aligned} V &= m_{11}^2 \Phi_1^\dagger \Phi_1 + m_{22}^2 \Phi_2^\dagger \Phi_2 - (m_{12}^2 \Phi_1^\dagger \Phi_2 + h.c.) \ &\quad + \frac{\lambda_1}{2} (\Phi_1^\dagger \Phi_1)^2 + \frac{\lambda_2}{2} (\Phi_2^\dagger \Phi_2)^2 + \lambda_3 (\Phi_1^\dagger \Phi_1)(\Phi_2^\dagger \Phi_2) + \lambda_4 (\Phi_1^\dagger \Phi_2)(\Phi_2^\dagger \Phi_1) \ &\quad + \left[\frac{\lambda_5}{2} (\Phi_1^\dagger \Phi_2)^2 + \lambda_6 (\Phi_1^\dagger \Phi_1)(\Phi_1^\dagger \Phi_2) + \lambda_7 (\Phi_2^\dagger \Phi_2)(\Phi_1^\dagger \Phi_2) + h.c.\right]. \end{aligned}7-physics and kaon-mixing bounds, while light-V=m112Φ1Φ1+m222Φ2Φ2(m122Φ1Φ2+h.c.) +λ12(Φ1Φ1)2+λ22(Φ2Φ2)2+λ3(Φ1Φ1)(Φ2Φ2)+λ4(Φ1Φ2)(Φ2Φ1) +[λ52(Φ1Φ2)2+λ6(Φ1Φ1)(Φ1Φ2)+λ7(Φ2Φ2)(Φ1Φ2)+h.c.].\begin{aligned} V &= m_{11}^2 \Phi_1^\dagger \Phi_1 + m_{22}^2 \Phi_2^\dagger \Phi_2 - (m_{12}^2 \Phi_1^\dagger \Phi_2 + h.c.) \ &\quad + \frac{\lambda_1}{2} (\Phi_1^\dagger \Phi_1)^2 + \frac{\lambda_2}{2} (\Phi_2^\dagger \Phi_2)^2 + \lambda_3 (\Phi_1^\dagger \Phi_1)(\Phi_2^\dagger \Phi_2) + \lambda_4 (\Phi_1^\dagger \Phi_2)(\Phi_2^\dagger \Phi_1) \ &\quad + \left[\frac{\lambda_5}{2} (\Phi_1^\dagger \Phi_2)^2 + \lambda_6 (\Phi_1^\dagger \Phi_1)(\Phi_1^\dagger \Phi_2) + \lambda_7 (\Phi_2^\dagger \Phi_2)(\Phi_1^\dagger \Phi_2) + h.c.\right]. \end{aligned}8 benchmarks can also produce negative values down to about V=m112Φ1Φ1+m222Φ2Φ2(m122Φ1Φ2+h.c.) +λ12(Φ1Φ1)2+λ22(Φ2Φ2)2+λ3(Φ1Φ1)(Φ2Φ2)+λ4(Φ1Φ2)(Φ2Φ1) +[λ52(Φ1Φ2)2+λ6(Φ1Φ1)(Φ1Φ2)+λ7(Φ2Φ2)(Φ1Φ2)+h.c.].\begin{aligned} V &= m_{11}^2 \Phi_1^\dagger \Phi_1 + m_{22}^2 \Phi_2^\dagger \Phi_2 - (m_{12}^2 \Phi_1^\dagger \Phi_2 + h.c.) \ &\quad + \frac{\lambda_1}{2} (\Phi_1^\dagger \Phi_1)^2 + \frac{\lambda_2}{2} (\Phi_2^\dagger \Phi_2)^2 + \lambda_3 (\Phi_1^\dagger \Phi_1)(\Phi_2^\dagger \Phi_2) + \lambda_4 (\Phi_1^\dagger \Phi_2)(\Phi_2^\dagger \Phi_1) \ &\quad + \left[\frac{\lambda_5}{2} (\Phi_1^\dagger \Phi_2)^2 + \lambda_6 (\Phi_1^\dagger \Phi_1)(\Phi_1^\dagger \Phi_2) + \lambda_7 (\Phi_2^\dagger \Phi_2)(\Phi_1^\dagger \Phi_2) + h.c.\right]. \end{aligned}9. In contrast, m112m_{11}^20 are much less informative in the studied parameter space: m112m_{11}^21 is dominated by long-distance two-photon physics, and m112m_{11}^22 changes by at most about m112m_{11}^23 (Hou et al., 2022).

Higgs-mediated quark flavor violation is another characteristic signature. In full scans of the general 2HDM consistent with perturbativity, vacuum stability, unitarity, oblique parameters, Higgs signal strengths, m112m_{11}^24 mixing, m112m_{11}^25, and direct m112m_{11}^26 limits, the allowed rates remain substantial. At m112m_{11}^27,

m112m_{11}^28

while at m112m_{11}^29,

m222m_{22}^20

If the mild m222m_{22}^21-mass-splitting discrepancy were confirmed and dominated by tree-level scalar exchange, the model would predict

m222m_{22}^22

This phenomenology is difficult to reproduce in vector-like-quark-only completions because correlated m222m_{22}^23-flavor violation suppresses the corresponding Higgs flavor-violating rates much more strongly (Herrero-Garcia et al., 2019).

The Higgs self-interactions are also unusually flexible. In a Higgs-basis analysis imposing tree-level unitarity, boundedness from below, the m222m_{22}^24-parameter constraint, and m222m_{22}^25, the cubic self-coupling m222m_{22}^26 of the observed Higgs may even flip sign, while the quartic coupling satisfies

m222m_{22}^27

The same analysis finds that the extra scalar masses are typically below about m222m_{22}^28 GeV for m222m_{22}^29 and below about Z2\mathbb{Z}_200 GeV for Z2\mathbb{Z}_201, unless alignment is approached very closely (Jurciukonis et al., 2018). This suggests a phenomenology in which alignment does not imply decoupling, but often coexists with a relatively light extended scalar spectrum.

6. Effective-field-theory form, UV origins, and conceptual boundaries

The general 2HDM is frequently treated as an effective field theory. The 2HDM EFT supplements the renormalizable Lagrangian with higher-dimensional operators and retains both doublets explicitly in the low-energy theory: Z2\mathbb{Z}_202 A complete Warsaw-like construction yields 228 dimension-six operators, together with three dimension-five neutrino-mass operators. Rewriting the EFT in the Higgs basis separates operators that modify masses and standard-model-like couplings from operators that contribute only to scattering amplitudes, and makes correlations among multi-Higgs processes more transparent. The same framework also provides specific operator subsets for the type-I, II, X, and Y symmetry-restricted models, with 76 operators common to all four types (Dermisek et al., 2024).

Several UV completions illuminate how the G2HDM can emerge nontrivially. The two Higgs bi-doublet left-right symmetric model has a genuine decoupling regime in which, for Z2\mathbb{Z}_203, the theory reduces to a general 2HDM with spontaneous CP violation and a light electroweak-scale two-doublet sector (Liu et al., 2012). A composite near-conformal construction based on dilaton EFT yields instead a highly constrained inert 2HDM corner, with

Z2\mathbb{Z}_204

and a specific hierarchy among quartic couplings, showing how strong dynamics can populate only a narrow subspace of the general parameter space (Appelquist et al., 2022). In 6D gauge-Higgs unification, two Higgs doublets can arise from extra-dimensional gauge fields, but the quartic potential is then a restricted subset of the full 2HDM potential, often of the form

Z2\mathbb{Z}_205

rather than the fully general scalar interaction (Chang et al., 2012).

These examples clarify an important conceptual point. The G2HDM is best viewed not as one model but as the maximal renormalizable 2HDM parameter space consistent with gauge symmetry. Symmetry-restricted 2HDMs, aligned models, BGL constructions, EFT deformations, and UV completions all occupy identifiable submanifolds inside that space. Equally important, the general 2HDM should not be conflated with the Gauged Two Higgs Doublet Model, which embeds the two doublets into an Z2\mathbb{Z}_206 gauge multiplet and uses the same acronym in a different sense (Huang et al., 2015, Yang et al., 2024). The former is a general low-energy field-theoretic framework; the latter is a distinct gauge extension.

Topic to Video (Beta)

No one has generated a video about this topic yet.

Whiteboard

No one has generated a whiteboard explanation for this topic yet.

Follow Topic

Get notified by email when new papers are published related to General Two Higgs Doublet Model (G2HDM).