⚛️ phenomenology

Perturbative unitarity for models with singlet and doublet scalars

This paper establishes a comprehensive framework for perturbative unitarity bounds in models featuring extra $SU(2)$ doublet, neutral singlet, and charged singlet scalars by proposing a minimal set of scattering matrices and providing a public Mathematica tool, BounDS, to implement and verify these constraints.

Original authors: Carolina T. Lopes, André Milagre, João P. Silva

Published 2026-03-03

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Carolina T. Lopes, André Milagre, João P. Silva

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant, incredibly complex video game. The "Standard Model" is the current rulebook that explains how the characters (particles) interact. It works great for most things, but it has a glaring plot hole: it can't explain Dark Matter, the invisible stuff that holds galaxies together.

To fix this plot hole, physicists invent new "expansion packs" for the game. These packs add new characters, often in the form of extra scalar particles (think of them as invisible, weightless ghosts that can pop in and out of existence). The paper you're asking about is a quality control manual for these expansion packs.

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

1. The Problem: "Too Much Power, Too Fast"

In physics, there's a rule called Unitarity. Think of it like a "probability budget." If you roll a die, the chances of getting a 1, 2, 3, 4, 5, or 6 must add up to exactly 100%. You can't have a 120% chance of something happening; that breaks the laws of logic.

When physicists add new particles to their theories, they have to make sure that if these new particles smash into each other at super-high speeds (like in the Large Hadron Collider), the math doesn't break. If the math says there's a 200% chance of a collision happening, the theory is broken. It's like a video game character having a "jump" power that is so high they fly off the map and crash the game engine.

2. The Solution: The "BounDS" Tool

The authors, Carolina, André, and João, realized that checking these rules for every possible new particle combination is a nightmare. It's like trying to manually check every single line of code in a massive video game to find bugs.

So, they built a digital mechanic called BounDS (a Mathematica notebook).

What it does: You tell the tool, "I want to add 2 new doublet particles and 1 singlet particle."
What it outputs: The tool automatically writes down all the possible ways these particles can crash into each other, calculates the "probability budget" for every crash, and tells you exactly which rules (coupling constants) you are allowed to break and which ones you must keep safe.

3. The "Traffic Light" System (Scattering Matrices)

To check if the game is stable, the authors had to organize the chaos. Imagine a busy intersection with cars (particles) coming from all directions.

The Old Way: You tried to watch every single car and every possible crash at once. It was confusing and messy.
The New Way (This Paper): They invented a smart traffic system. They grouped cars by their License Plate (Charge), Color (Hypercharge), and Size (Isospin).
- Instead of watching a red car crash into a blue car, they realized that a "Red-Size-1" car crashing into a "Red-Size-1" car behaves exactly the same as a "Blue-Size-1" car crashing into a "Blue-Size-1" car.
- By grouping them this way, they reduced a massive, impossible-to-solve puzzle into a few small, manageable checklists (called scattering matrices).

4. The "Speed Limit" Signs

Once they organized the traffic, they calculated the Speed Limits.

In their paper, they derived specific mathematical "speed limits" for the strength of the interactions between these new particles.
If a model suggests the particles interact too strongly (like a car going 500 mph in a school zone), the theory is invalid. The authors provide the exact speed limit signs (equations) for dozens of different "expansion packs" (models).
They even checked their work against old maps (previous scientific papers) and found that their new GPS (BounDS) agrees with the old maps but is much faster and covers more territory.

5. Why Should You Care?

You might think, "I don't care about scalar singlets." But here is the metaphor:

Dark Matter is the biggest mystery in physics. We know it's there, but we don't know what it is.
Most theories about Dark Matter involve adding these extra "ghost" particles.
Before we build a giant machine to hunt for Dark Matter, we need to know which theories are even possible.
This paper is the filter. It tells theorists: "Don't waste your time building a model with these specific settings; it breaks the laws of physics. But if you tweak the settings this way, it might work."

Summary

The authors created a universal rulebook and a software tool to ensure that any new theory trying to explain Dark Matter (or other mysteries) doesn't break the fundamental laws of probability. They organized the chaos of particle collisions into neat categories, calculated the safety limits, and gave the scientific community a free, automated tool to check their work.

In short: They built a "Safety Inspector" for the universe's expansion packs, ensuring that the new characters we add to our cosmic story don't crash the simulation.

1. Problem Statement

The Standard Model (SM) of particle physics, while successful, fails to explain Dark Matter (DM). Consequently, many Beyond the Standard Model (BSM) scenarios introduce new scalar fields, typically SU(2) doublets (like in the Two-Higgs-Doublet Model, 2HDM) and SU(2) singlets (neutral or charged).

For any consistent quantum field theory, perturbative unitarity must be preserved. This ensures that scattering probabilities do not exceed unity at high energies. In scalar extensions, the presence of additional degrees of freedom can cause scattering amplitudes to grow with energy, violating unitarity unless the quartic couplings of the scalar potential are constrained.

The Challenge:

Deriving these constraints manually for models with arbitrary numbers of doublets and singlets is tedious and error-prone.
Existing literature often treats specific cases (e.g., 2HDM, specific singlet extensions) individually, lacking a unified framework for general combinations of $n_D$ doublets, $n_c$ charged singlets, and $n_n$ neutral singlets.
Identifying the minimal set of independent scattering channels is non-trivial due to redundancies arising from gauge symmetries and identical particle statistics.

2. Methodology

The authors develop a systematic framework to derive perturbative unitarity bounds for the most general renormalizable scalar potential involving an arbitrary number of doublets and singlets.

A. Theoretical Framework

Model Definition: The study considers an $SU(2) \times U(1)$ $S U (2) \times U (1)$ gauge theory with:
- $n_D$ complex scalar doublets ( $\Phi_i$ ) with hypercharge $Y=1/2$ .
- $n_c$ complex scalar singlets ( $\phi^+_i$ ) with $Y=1$ .
- $n_n$ real scalar singlets ( $\chi_i$ ) with $Y=0$ .
High-Energy Limit: Calculations are performed in the high-energy limit ( $s, t, u \to \infty$ ). In this regime, contributions from propagators ( $s, t, u$ -channels) vanish, and the scattering amplitude is dominated solely by the quartic interactions of the scalar potential.
Partial Wave Analysis: The $2 \to 2$ scattering amplitudes are decomposed into partial waves. The strongest constraint comes from the zeroth partial wave ( $a_0$ ). The unitarity condition requires $|a_0| \leq 1/2$ , which translates to $|M| \leq 8\pi$ for the scattering matrix eigenvalues.

B. Classification and Basis Construction

A key methodological innovation is the classification of two-particle states.

Quantum Numbers: States are labeled by $|Q, Y, T\rangle$ , where $Q$ is electric charge, $Y$ is hypercharge, and $T$ is total isospin.
Redundancy Reduction: The authors demonstrate that labeling by $Q$ and $Y$ alone (as done in some previous works) leads to redundant scattering matrices. By including total isospin $T$ , they identify a minimal basis of independent channels.
Scattering Matrices: They derive explicit formulas for the scattering matrices ( $M_{|Q,Y,T\rangle}$ ) by taking the fourth derivative of the quartic potential with respect to the fields. They account for symmetry factors ( $1/\sqrt{2}$ ) for identical particles.

C. Automation Tool (BounDS)

To handle the combinatorial complexity of models with many fields and symmetries, the authors developed a Mathematica notebook called BounDS.

Input: The user specifies the number of fields ( $n_D, n_c, n_n$ ) and any additional discrete or continuous symmetries (e.g., $Z_2, Z_3, CP$ ).
Process: The tool automatically:
1. Constructs the most general renormalizable quartic potential consistent with the symmetries.
2. Identifies the minimal set of independent quartic couplings.
3. Constructs the 7 independent scattering matrices corresponding to the minimal basis states.
4. Block-diagonalizes these matrices.
5. Computes eigenvalues and outputs the unitarity bounds ( $|\lambda_{eigen}| \leq 8\pi$ ).

3. Key Contributions

Unified Formalism: Provided a complete, general description of perturbative unitarity bounds for models with arbitrary numbers of $SU(2)$ doublets and singlets (neutral and charged).
Minimal Basis Identification: Proved that using the $|Q, Y, T\rangle$ classification yields a minimal set of scattering matrices, eliminating redundant channels and simplifying calculations compared to previous $|Q, Y\rangle$ approaches.
Generalized Coupling Relations: Derived the full set of relations among quartic couplings required by Hermiticity and gauge invariance for this general class of models.
Open-Source Tool (BounDS): Released a public Mathematica package that automates the derivation of unitarity bounds, significantly lowering the barrier for phenomenologists to test new BSM models.
Comprehensive Validation: Applied the framework to numerous specific cases (SM, 2HDM, 3HDM, singlet extensions) and verified that the results match or generalize existing literature.

4. Results

The paper presents explicit scattering matrices and unitarity bounds for several benchmark models:

Standard Model (SM): Recovers the classic bound $\lambda \leq 4\pi/3$ .
SM + 1 Neutral Singlet: Derives bounds involving the mixing coupling $\gamma$ and singlet self-coupling $\beta$ .
2HDM (with and without $Z_2$ symmetry): Recovers the well-known bounds on $\lambda_{1-7}$ , confirming agreement with Refs. [17, 18].
General 2HDM: Provides the full set of matrices for the unconstrained 2HDM (including complex $\lambda_{6,7}$ ).
Multi-Singlet Extensions: Analyzed models with 2 neutral singlets and 2HDM + 1 singlet, deriving complex eigenvalue conditions.
Charged Singlets: Treated the case of 2 doublets + 1 neutral singlet + 1 charged singlet, explicitly including the $\kappa$ couplings (which mix doublets and singlets via the $\sigma_2$ matrix).
3HDM with $Z_3$ symmetry: Derived bounds for a 3-doublet model, showing how symmetries reduce the dimensionality of the scattering matrices.

In all cases, the authors provided the explicit eigenvalues (where analytically solvable) or the matrix forms required for numerical diagonalization.

5. Significance

Theoretical Consistency: The work provides a necessary "sanity check" for any BSM model involving scalars. Before simulating a model's phenomenology (e.g., DM relic density or collider signatures), one must ensure the parameter space satisfies unitarity.
Efficiency: The BounDS tool automates a process that previously required manual derivation of Lagrangians and scattering matrices, which is prone to human error and extremely time-consuming for complex models.
Dark Matter Phenomenology: Since many Dark Matter models rely on scalar singlets or inert doublets, this framework allows for a rigorous exploration of the viable parameter space for these theories, ensuring they remain within the perturbative regime.
Future-Proofing: The methodology is scalable. As new models with higher multiplets or complex symmetry structures are proposed, the BounDS framework can be immediately applied without re-deriving the fundamental scattering formalism.

In summary, this paper establishes a robust, automated, and generalizable standard for enforcing perturbative unitarity in scalar extensions of the Standard Model, bridging the gap between abstract theoretical constraints and practical model-building tools.