The Hilbert Series and the Flavor Invariants of the 3HDM

This paper presents a systematic study of invariant operators in the three-Higgs-doublet model (3HDM) by computing the associated Hilbert series and constructing explicit expressions for flavor invariants up to cubic order in the couplings.

The Gist

Imagine you are a master chef trying to create a new, complex dish. You have three types of special ingredients (let's call them "Higgs Doublets") and a massive pantry of spices (the "couplings" and "masses"). Your goal is to figure out every possible unique flavor combination you can make that tastes the same no matter how you rotate your plate or rearrange the ingredients on the table.

This paper is essentially a recipe book and a counting machine for a theoretical physics model called the Three-Higgs-Doublet Model (3HDM).

Here is a breakdown of what the authors did, using simple analogies:

1. The Problem: A Kitchen Too Big to Cook In

In the standard model of physics, we usually deal with one or two of these "ingredients." But in this model, there are three.

• The Challenge: When you have three ingredients, the number of ways to mix them explodes. It's like trying to count every possible sentence you can write with a dictionary of 10,000 words.
• The Goal: Physicists need to know which combinations of these ingredients represent "real" physical laws. These combinations must be invariant, meaning they don't change even if you look at the system from a different angle (a concept called "symmetry"). If a combination changes when you rotate the view, it's just an illusion, not a fundamental law.

2. The Counting Machine: The "Hilbert Series"

The authors first built a mathematical tool called a Hilbert Series.

• The Analogy: Imagine you have a giant vending machine. You put in a coin (representing a specific type of ingredient), and the machine tells you exactly how many unique, valid recipes you can make with that coin.
• The Difficulty: Because the 3HDM is so complex, this "vending machine" was jammed. The math involved was like trying to solve a puzzle where the pieces keep changing shape and multiplying into thousands of new pieces. The authors had to invent a new way to "unclog" the machine, using clever computer tricks to handle the massive numbers without crashing their software.
• The Result: They successfully counted every single unique recipe possible. They found that the "numerator" (the list of recipes) is so huge it has 31 million terms. It's too big to print in the paper, so they put it on a website, like a massive digital library.

3. The Recipe Book: Building the Invariants

Counting the recipes is one thing; actually writing them down is another. The authors then set out to write the actual "recipes" (the mathematical formulas) for the most important combinations.

• The Strategy: Instead of trying to write down every single recipe at once (which would take forever), they used a technique called the "Background Field" method.
• The Analogy: Imagine you want to find all the symmetrical patterns on a spinning globe. It's hard to see them while it's spinning. So, you freeze the globe in one specific position (the "background"). Now, the patterns are easier to spot. Once you find the patterns in this frozen state, you use math to "unfreeze" the globe and prove that those patterns hold true even when it spins.
• The Outcome: They wrote down the explicit formulas for the most important combinations (up to a certain complexity level). They organized them like a menu:
  • Zero "Special" Ingredients: Simple mixes.
  • One "Special" Ingredient: Slightly more complex.
  • Two or Three "Special" Ingredients: The most complex, exotic dishes.

4. Why Does This Matter?

You might ask, "Why do we need a list of 31 million recipes?"

• Finding New Physics: Just like a chef might discover a new flavor by mixing ingredients in a way no one has before, physicists use these "invariant recipes" to look for new particles or forces.
• Dark Matter: Some of these recipes might explain what Dark Matter is (the invisible stuff holding galaxies together). If a recipe is stable and doesn't decay, it could be a dark matter candidate.
• Time Travel (Sort of): Some recipes might explain why time moves forward and not backward (CP violation), which is a mystery in our universe.
• The "Universal Language": By writing these recipes in a way that doesn't depend on how you label your ingredients, the authors created a "universal language" for this model. Any scientist, anywhere, can use this list to check if their theory is valid without getting confused by different naming conventions.

Summary

Think of this paper as the ultimate index card system for a cosmic kitchen.

1. The Problem: The kitchen is too big to navigate.
2. The Solution: They built a super-computer to count every possible dish (Hilbert Series).
3. The Application: They then wrote down the instructions for the most delicious and stable dishes (Invariant Operators).
4. The Payoff: This gives physicists a map to navigate the complex world of three Higgs fields, helping them hunt for dark matter, understand the universe's symmetry, and potentially discover new physics beyond what we currently know.

The authors didn't just count the stars; they drew a map of the entire galaxy so others can explore it.

Technical Summary

1. Problem Statement

The Three-Higgs-Doublet Model (3HDM) is a significant extension of the Standard Model (SM) that introduces three scalar doublets. While the Two-Higgs-Doublet Model (2HDM) has been extensively studied, the 3HDM presents a vastly more complex parameter space and symmetry structure.

• The Core Challenge: Identifying physical observables in the 3HDM is difficult because the theory possesses a large global $SU(3)$ flavor symmetry. Physical quantities (such as CP-violating observables or renormalization group invariants) must be expressible as operators invariant under this symmetry and independent of the basis chosen for the Higgs fields.
• The Gap: Unlike the 2HDM, where a minimal generating set of invariants was established using bilinear and "birdtrack" techniques, no systematic study existed for the 3HDM. The sheer number of parameters and the complexity of the $SU(3)$ representations (specifically the presence of a 27-dimensional representation) made a full classification of invariant operators computationally prohibitive using standard methods.

2. Methodology

The authors employed a two-pronged approach combining algebraic geometry (Hilbert Series) and background field techniques to solve the problem.

A. Hilbert Series Calculation

The Hilbert series is a generating function that counts the number of independent invariant operators at a given order.

• Representation Content: The 3HDM scalar parameters decompose into three adjoint representations ($\mathbf{8}$) and one 27-dimensional representation ($\mathbf{27}$) of the global $SU(3)$ flavor group.
• Technical Obstacles: Direct computation of the Hilbert series via contour integration (using the Molien-Weyl formula) faced severe hurdles:
  1. High-Order Poles: The integrand contained quadratic and quartic factors in the denominator, leading to high-order poles requiring complex derivative evaluations.
  2. Branch Cuts: The order of integration introduced branch points and cuts, complicating the residue theorem application.
  3. Combinatorial Explosion: The number of poles grew into the thousands, making symbolic summation in standard software (Mathematica, Python/SymPy) memory-intensive and unstable.
• Novel Computational Strategy: To overcome these, the authors developed a hybrid numerical-symbolic approach:
  1. Derivative Trick: They replaced high-order poles with derivatives with respect to auxiliary variables, converting the problem into one involving only simple poles.
  2. Partial Fraction Decomposition: They decomposed the integrand into a sum of simpler terms, reducing the problem to calculating eight distinct types of contour integrals.
  3. Coefficient Array Manipulation: Recognizing that the resulting polynomials had coefficients as large as $O(10^{27})$ and thousands of terms, they abandoned standard symbolic polynomial manipulation. Instead, they represented the numerator and denominator as multi-dimensional arrays of arbitrary-precision Python integers. They performed polynomial multiplication and division (synthetic division) directly on these arrays to handle the massive scale of the data.

B. Construction of Explicit Invariants

To find the explicit form of the operators, the authors utilized a background field method:

1. Symmetry Breaking: They moved to the mass basis where the mass matrix is diagonal. This breaks the global $SU(3)$ symmetry down to a residual $U(1) \times U(1)$ symmetry.
2. Charge Analysis: They determined the $U(1) \times U(1)$ charges of the field components (specifically the quartic couplings $\lambda$).
3. Matching: They identified all polynomials invariant under the reduced $U(1) \times U(1)$ symmetry.
4. Reconstruction: By treating the mass terms as background fields, they systematically matched these reduced invariants to full $SU(3)$-invariant structures. They constructed a basis of independent invariants by ensuring the generated set could reproduce all $U(1) \times U(1)$ invariants, using linear algebra tools (NullSpace and SubspaceBasis) to eliminate linear dependencies.

3. Key Contributions

1. First Full Hilbert Series for 3HDM: The paper presents the first complete calculation of the multigraded Hilbert series for the 3HDM, encoding the number and structure of all invariant operators associated with the $\mathbf{8} \oplus \mathbf{8} \oplus \mathbf{8} \oplus \mathbf{27}$ representation content.
2. Computational Breakthrough: The authors developed and demonstrated a robust algorithm for computing Hilbert series for large representations involving high-order poles and massive polynomial expansions, a technique applicable to other complex models in particle physics.
3. Explicit Basis of Invariants: They constructed an explicit basis of independent invariant operators up to cubic order in the quartic couplings ($O(\lambda^3)$).
   • For invariants involving only the adjoint representations ($\mathbf{8}$), they provided results to arbitrary order.
   • For invariants involving the $\mathbf{27}$ representation, they provided the full set up to $O(\lambda^3)$.
4. Categorization: The invariants are systematically categorized by the number of $\mathbf{27}$ fields ($W$) and $\mathbf{8}$ fields ($V$) they contain, providing a clear organizational structure for future phenomenological work.

4. Results

• Hilbert Series: The authors derived the rational function $H(s, t, u, q) = N/D$, where $s, t, u$ correspond to the three $\mathbf{8}$ representations and $q$ to the $\mathbf{27}$. The denominator $D$ is a product of factors reflecting the relations among the invariants. The numerator $N$ contains over 31 million terms (fully expanded), which is too large to print but is available in an online repository. The expansion of the series matches previous ungraded results, serving as a consistency check.
• Invariant Operators:
  • Zero 27s: A complete set of invariants built solely from the adjoint fields ($V$) was found to arbitrary order, including traces of products like $\text{Tr}(V^I V^J)$, $\text{Tr}(V^I V^J V^K)$, and commutator-based structures.
  • With 27s: The paper lists the minimal generating set for operators involving the $\mathbf{27}$ field ($W$) coupled to $V$ fields up to cubic order. Examples include structures like $W^{kl}_{ij} (V^I)^i_k (V^J)^j_l$ and higher-order contractions.
  • Independence: The authors verified that their constructed basis is linearly independent and sufficient to generate all possible invariants up to the specified order.

5. Significance

• Phenomenological Utility: The identified invariants provide a basis-independent framework for analyzing the 3HDM. This is crucial for:
  • CP Violation: Constructing CP-odd invariants to distinguish between explicit and spontaneous CP violation without relying on specific field bases.
  • Dark Matter: Identifying stable scalar candidates stabilized by residual symmetries (e.g., $S_3$ or unbroken $U(1)$).
  • Flavor Physics: Understanding the connection between the 3HDM symmetry structure and the flavor structure of quarks and leptons.
• Model Building: The results offer the necessary building blocks for constructing effective Lagrangians and studying the renormalization group evolution of the 3HDM.
• Methodological Impact: The computational techniques developed (handling massive polynomial expansions via array manipulation and resolving high-order pole integrals) are transferable to other areas of theoretical physics involving large symmetry groups and complex invariant rings, such as Supersymmetric QCD or higher-dimensional Higgs sectors ($N > 3$).

In conclusion, this work resolves a long-standing technical bottleneck in the study of multi-Higgs models, providing both a theoretical map (via the Hilbert series) and practical tools (explicit invariants) for exploring the rich phenomenology of the Three-Higgs-Doublet Model.