New features in the Z2xZ2 3HDM two component DM model

This paper investigates the vacuum structure and phenomenology of a $Z_2 \times Z_2$ symmetric three Higgs doublet model with two inert scalar doublets, revealing new viable parameter space where two-component dark matter candidates contribute comparably to the observed relic density and offer distinctive experimental signatures.

The Gist

Imagine the universe as a giant, complex machine. For a long time, scientists have used a "Standard Model" blueprint to explain how the visible parts of this machine work (like stars, planets, and us). But we know there's a huge missing piece: Dark Matter. We can't see it, but we know it's there because of its gravity.

This paper is like a team of architects (Jorge, Rafael, Pedro, and Joao) proposing a new, more detailed blueprint to explain how Dark Matter works. Instead of just adding one new part to the machine, they are adding three new layers (called "Higgs doublets") and using a specific set of rules (a "Z2 x Z2 symmetry") to keep everything stable.

Here is the breakdown of their work in simple terms:

1. The Setup: Two "Invisible" Twins

In their model, two of these new layers are "inert." Think of them like ghosts that don't interact with normal light or matter, but they do have mass. Because they are "inert," they are stable and can't just disappear. This makes them perfect candidates for Dark Matter.

Usually, scientists look for a single type of Dark Matter particle. But this model suggests a two-component scenario: there are two different types of these "ghost" particles (let's call them Ghost A and Ghost B) living in the same universe.

2. The Challenge: Finding the "Lowest Point"

Imagine a hilly landscape. The universe wants to settle in the deepest valley (the lowest energy state). If it settles in the wrong valley, the whole model falls apart.

The authors spent a lot of time mapping this landscape. They found:

• Old Map: Scientists previously knew about a few valleys where the universe could settle.
• New Discovery: The authors found two new valleys (which they named F0DM0' and F0CB) that no one had noticed before.
• The Goal: They had to prove that the "Two Ghost" valley is the deepest one (the global minimum). If the universe falls into a different valley, our model of reality breaks. They used complex math to ensure the "Two Ghost" valley is indeed the winner.

3. The Rules of the Game

Before they could say their model works, they had to check if it followed the "laws of physics" (like not having infinite energy or breaking the speed of light). They ran a massive simulation (a "scan") checking:

• The Big Bang: Does it match the amount of Dark Matter we see today? (The answer is yes).
• The Particle Collider (LHC): Have we already smashed these particles at the Large Hadron Collider? (They checked and found regions where we haven't seen them yet, so the model is still safe).
• Direct Detection: If we try to catch these ghosts in a lab, will they bounce off our detectors? They checked this against current experiments (like LZ) and future ones (like DARWIN).

4. The Big Surprise: A Perfect Team

The most exciting finding is about how the two ghosts share the job.

• The Old Idea: Usually, one ghost does all the work, and the other is just a bystander.
• The New Finding: The authors found specific scenarios where Ghost A and Ghost B share the workload equally. They both contribute 50% to the total amount of Dark Matter in the universe.

This is like a relay race where, instead of one runner doing the whole lap, two runners split the distance perfectly. This creates a unique "signature" that future experiments might be able to spot.

5. The Mass Range

They found that these Dark Matter particles can have almost any weight, from very light (half the weight of the Higgs boson) to very heavy (1,000 times heavier).

• If they are light, they might be hiding in a specific mass range.
• If they are heavy, they might be hiding in another.
• Crucially, they found that even if one particle is light and the other is heavy, they can still work together to create the right amount of Dark Matter.

Summary

The authors built a more complex, robust version of the Dark Matter model. They mapped out all the possible ways the universe could settle down, found two new possibilities, and proved that a model with two different types of Dark Matter particles working together is a very viable option.

They didn't just say "it's possible"; they showed exactly where to look for it in the data and highlighted that the future might reveal a universe where Dark Matter isn't a solo act, but a duet.

Technical Summary

Technical Summary: New features in the Z2×Z2 3HDM two component DM model

Problem Statement
The paper investigates the phenomenology and theoretical consistency of a Three Higgs Doublet Model (3HDM) extended by a $Z_2 \times Z_2$ symmetry. This specific configuration introduces two inert scalar doublets, creating a two-component Dark Matter (DM) scenario. The primary challenge addressed is the rigorous determination of the model's vacuum structure to ensure the "double inert" vacuum (where only the SM-like Higgs doublet acquires a vacuum expectation value, while the other two remain inert) is the global minimum. Furthermore, the authors aim to map the viable parameter space under comprehensive theoretical constraints (boundedness, unitarity) and experimental limits (collider data, relic density, direct and indirect detection), specifically looking for regions where two distinct DM candidates contribute comparably to the observed relic abundance.

Methodology
The authors employ a multi-stage analytical and numerical approach:

1. Potential Analysis and Boundedness:
   The scalar potential is defined with 15 parameters. To ensure the potential is Bounded From Below (BFB), the authors utilize a parameterization of the scalar fields and decompose the quartic part of the potential into neutral, charge-breaking, and $Z_2 \times Z_2$ specific terms. They derive sufficient conditions for BFB by constructing a lower-bound potential and analyzing the copositivity of a specific $3 \times 3$ matrix ($A$) formed from the quartic couplings.

2. Vacuum Structure Identification:
   The study systematically identifies all possible minima, categorizing them into neutral and charge-breaking configurations.

   • Neutral Vacua: The authors review known configurations (e.g., $2$-Inert, $2$HDM + 1 DM) and identify a new neutral vacuum configuration, denoted as F0DM0' ($(v_1, iv_2, 0)$).
   • Charge-Breaking Vacua: Similarly, they list all possible charge-breaking vacua and identify a new configuration, F0CB.
   • Global Minimum Verification: To guarantee the 2-Inert vacuum is the global minimum, the authors compare the potential values of all identified minima (including the new ones) against the 2-Inert case. They employ numerical minimization using the Minuit algorithm with random initial conditions to verify that no lower vacuum exists.

3. Constraint Implementation:
   The model is subjected to a rigorous set of constraints:

   • Theoretical: Perturbative unitarity (S-matrix), BFB conditions, and electroweak precision parameters ($S, T, U$).
   • Collider: Higgs signal strengths ($h_{125}$) from ATLAS, limits on additional scalars from LEP, Tevatron, and LHC (via HiggsTools), and limits on the total Higgs decay width.
   • Dark Matter:
     • Relic Density: The sum of contributions from both DM candidates ($H_1$ and $H_2$) must match the Planck value ($\Omega_{DM}h^2 \approx 0.12$).
     • Direct Detection: Spin-independent scattering cross-sections are compared against current (LZ) and projected (DARWIN/XLZD, PandaX-xT) limits, rescaled by the relative relic density of each component.
     • Indirect Detection: Constraints from Fermi-LAT (dwarf galaxies), AMS-02 (antiprotons), and H.E.S.S. (Galactic Center) are applied to the thermally averaged annihilation cross-section $\langle \sigma v \rangle$.

4. Parameter Scan:
   A FORTRAN program generates random parameter sets within defined ranges (masses from 50–1000 GeV, couplings within $\pm [10^{-3}, 10]$). The model is implemented in micrOMEGAs 6.0.5 via FeynRules and CalcHEP to calculate observables and test constraints.

Key Contributions

• New Vacuum Configurations: The authors identify and characterize two previously unexplored vacuum configurations: F0DM0' (a neutral vacuum with no DM candidate and massless fermions) and F0CB (a charge-breaking vacuum). They provide explicit expressions for these cases to facilitate direct comparison with the global minimum.
• Comprehensive Vacuum Stability: The work establishes a robust procedure to ensure the 2-Inert vacuum is the absolute global minimum by including these new configurations in the stability analysis.
• Two-Component DM Phenomenology: The study demonstrates that the $Z_2 \times Z_2$ 3HDM allows for a broader range of viable scenarios than previously explored, specifically where two DM components contribute significantly to the relic density.

Results

• Mass Spectrum Viability: The analysis reveals that viable DM candidate masses can exist across the entire range $[\frac{1}{2}m_h, 1000 \text{ GeV}]$.
• Relic Density Scenarios:
  • In the intermediate mass range, it is possible for one component (specifically $H_2$) to dominate the relic density while the other plays a secondary role.
  • Equal Contribution: The authors identify specific regions where both DM candidates contribute equally to the observed relic density. This occurs either when $m_{H_1}$ is in the range $\frac{1}{2}m_h < m_{H_1} < 80 \text{ GeV}$ or when $m_{H_1} \gtrsim 500 \text{ GeV}$.
• Detection Prospects:
  • Direct Detection: For low $m_{H_1}$, direct detection experiments can probe $H_1$ without necessarily constraining $H_2$. The authors note that future experiments (DARWIN) may reach the "neutrino floor" for high-mass regions, necessitating complementary probes.
  • Indirect Detection: The study finds that the constraints from the Planck relic density measurement are stringent enough that indirect detection limits generally do not further exclude the viable parameter space.
• Constraints: The model remains consistent with all current collider constraints, including Higgs signal strengths and searches for additional scalars.

Significance and Claims
The paper claims that simplistic conclusions drawn from narrow regions of the parameter space in multi-component DM models are insufficient. By rigorously analyzing the full vacuum structure (including new minima) and scanning the parameter space, the authors reveal a "much broader and more diverse range of possibilities."

The primary significance lies in demonstrating that the $Z_2 \times Z_2$ 3HDM can accommodate scenarios where two distinct DM particles contribute comparably to the relic density, offering distinctive experimental signatures. The work highlights that while direct detection experiments will probe significant portions of the parameter space, the model's viability in high-mass regions (near the neutrino floor) requires complementary detection strategies. The identification of new vacuum configurations ensures that previous stability analyses were incomplete, and the updated constraints provide a more accurate map of the model's phenomenological viability.