⚛️ phenomenology

Two-loop QCD corrections to $CP$ -even Higgs boson decays into $V\gamma$ in the Type-I Two-Higgs-Doublet Model

This paper presents a comprehensive analysis of next-to-leading order QCD corrections to the loop-induced decays of the heavy CP-even Higgs boson into $Z\gamma$ and $\gamma\gamma$ within the Type-I Two-Higgs-Doublet Model, revealing that these corrections are dominated by top-quark loops, typically range from -20% to +20% (with larger effects in non-alignment scenarios), and may become crucial for interpreting future HL-LHC exclusion bounds despite current data not yet constraining the parameter space.

Original authors: Lin-Qiao Qu, Zhong-Yuan Liu, Peng-Fei Duan, Yu Zhang, Wei-Long Duan

Published 2026-09-01

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Lin-Qiao Qu, Zhong-Yuan Liu, Peng-Fei Duan, Yu Zhang, Wei-Long Duan

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

In the vast, invisible landscape of particle physics, scientists are constantly searching for the rules that govern how matter and energy interact. At the heart of this exploration sits a particle known as the Higgs boson, a fundamental piece of the universe that gives other particles their mass. While the standard version of this particle was discovered over a decade ago, theorists have long suspected that nature might be hiding a heavier, more elusive cousin. This hypothetical heavier particle would behave similarly to its lighter counterpart but could hold the key to understanding mysteries that the current laws of physics cannot explain, such as why the universe is made of matter rather than antimatter. To find this hidden partner, researchers look for the rare moments when these particles decay, or break apart, into flashes of light or other particles. Because these decays happen through complex, invisible loops of virtual particles, even tiny changes in the underlying laws of nature can leave a distinct fingerprint on the outcome, making these rare events powerful windows into the unknown.

A team of researchers has recently taken a deep dive into the theoretical behavior of this heavier Higgs boson, specifically focusing on how it might decay into a photon, a particle of light, and a Z boson, a heavy carrier of the weak force, or into two photons. They worked within a specific theoretical framework called the Type-I Two-Higgs-Doublet Model, which proposes that the universe contains two sets of fields that generate mass, rather than just one. By running sophisticated computer simulations, the team calculated the probability of these decays occurring with a high degree of precision, accounting for the strong nuclear force that binds quarks together. Their work reveals that the heavier Higgs boson is far more likely to decay into two photons than into a photon and a Z boson, with the two-photon channel occurring roughly four times more often. This difference is driven by the mass of the top quark, the heaviest known elementary particle, which dominates the process, while the contributions from other particles, such as the W boson, effectively vanish under the specific conditions the team analyzed.

The researchers discovered that the likelihood of these decays is not a fixed number but changes depending on the mass of the heavy Higgs boson and a parameter known as tan beta, which describes the relationship between the two mass-generating fields. As the value of tan beta increases, the probability of the decay drops steadily. Conversely, as the mass of the heavy Higgs boson increases, the decay rate generally grows because there is more room for the resulting particles to move. A particularly interesting finding emerged when the mass of the heavy Higgs boson approached the combined mass of a top quark and an anti-top quark. In this specific energy range, the researchers observed a noticeable peak in the correction factors, a phenomenon where the calculated probability shifts by about fifteen percent. This peak occurs because the particles involved are momentarily created in a specific quantum state that allows them to interact strongly, yet the team confirmed that the mathematical description of this interaction remains stable and does not break down, ensuring the predictions are reliable.

To ensure their findings were robust, the team tested thousands of different possible scenarios within their model, checking them against known experimental limits and theoretical constraints. They found that in the vast majority of viable scenarios, the corrections they calculated were relatively small, typically ranging between negative twenty percent and positive twenty percent. However, in certain specific configurations where the model deviates slightly from the standard alignment, these corrections can become much larger, reaching up to thirty-eight percent in the two-photon channel. This level of precision is crucial because current experiments at the Large Hadron Collider have not yet found evidence of this heavy Higgs boson, but they are getting closer. The team compared their new, more accurate predictions with the latest data from the ATLAS experiment and found that while the current limits do not yet rule out the existence of this particle, the improved calculations will be essential for interpreting future data. As the collider prepares to collect even more data in the coming years, these refined theoretical predictions will help scientists determine whether the subtle signals they are looking for are truly signs of new physics or just a reflection of the complex, well-understood forces of the standard model.

Technical Summary: Two-loop QCD corrections to CP-even Higgs boson decays into $V\gamma$ in the Type-I Two-Higgs-Doublet Model

Problem Statement
The discovery of the 125 GeV Higgs boson ( $h$ ) at the LHC has established the Standard Model (SM) as a robust framework, yet loop-induced decay channels such as $h \to V\gamma$ ( $V=Z, \gamma$ ) remain sensitive probes for Beyond the Standard Model (BSM) physics. While the SM predictions for these decays are known to high precision (up to N $^3$ LO/N $^4$ LO for $h \to \gamma\gamma$ ), theoretical predictions for heavy neutral Higgs bosons ( $H$ ) in extended scalar sectors, specifically the Two-Higgs-Doublet Model (THDM), often rely on approximations. Previous studies of heavy Higgs decays in the THDM have frequently utilized the large-top-quark-mass approximation, restricting their validity to light Higgs masses, or have omitted specific channels like $H \to Z\gamma$ . Furthermore, the behavior of these decays near the top-quark pair ( $t\bar{t}$ ) production threshold requires careful scrutiny regarding the reliability of fixed-order perturbation theory, particularly concerning Coulomb corrections. This work addresses the need for a complete, exact mass-dependent calculation of the Next-to-Leading Order (NLO) QCD corrections to the decays $H \to \gamma\gamma$ and $H \to Z\gamma$ within the Type-I THDM, retaining full quark mass dependence and analyzing the perturbative stability near the $t\bar{t}$ threshold.

Methodology
The authors perform a complete two-loop QCD calculation for the decay amplitudes of a heavy CP-even Higgs boson ( $H$ ) into a photon and a $Z$ boson or two photons.

Framework: The analysis is conducted within the Type-I Two-Higgs-Doublet Model, where all fermions couple to a single Higgs doublet. The study focuses on the alignment limit ( $s_{\beta-\alpha} = 1$ ), where the light Higgs couplings mimic the SM, and extends to scenarios with slight deviations from this limit.
Calculation Setup: The calculation is performed in the 't Hooft–Feynman gauge using the on-shell (OS) renormalization scheme. The amplitude is decomposed into tensor structures, with the Ward identity ensuring the cancellation of ultraviolet divergences.
Tools: The Feynman diagrams (31 for LO $H \to \gamma\gamma$ and 50 for LO $H \to Z\gamma$ ; 24 two-loop diagrams for NLO) are generated using FeynArts 3.11 and simplified with FeynCalc 10.0.0. One-loop integrals are evaluated using Package-X. For the two-loop NLO QCD corrections, Integration-by-Parts (IBP) reduction is performed with Kira 2.3, and the resulting master integrals are computed numerically using the Auxiliary Mass Flow (AMF) method via the AMFlow package.
Parameter Space: A scan of 5,000 viable parameter points is conducted, subject to theoretical constraints (vacuum stability, perturbativity, unitarity) and experimental constraints (electroweak precision observables, Higgs coupling measurements, direct searches for charged Higgs). The scan covers $m_H \in [126, 1000]$ GeV, $\tan\beta \in [1.5, 50]$ , and assumes mass degeneracy among non-SM Higgs bosons ( $m_H = m_A = m_{H^\pm}$ ).

Key Contributions and Results

Dominance of Top-Quark Loops: In the alignment limit, the $W$ -boson contribution to the $H \to V\gamma$ amplitude vanishes. The decay widths are dominated by the top-quark loop, with charged-Higgs contributions being suppressed by the effective trilinear coupling and the factor $\cot^2\beta$ . Consequently, the decay widths scale as $\cot^2\beta$ , decreasing monotonically as $\tan\beta$ increases.
Relative Decay Rates: The decay width for $H \to \gamma\gamma$ is found to exceed that of $H \to Z\gamma$ by a factor of approximately four across the viable parameter space.
NLO QCD Corrections: The QCD corrections span a range of approximately $-20\%$ $- 20%$ to $+20\%$ $+ 20%$ .
- Mass Dependence: The corrections exhibit a strong dependence on the heavy Higgs mass $m_H$ . They attain a maximum of roughly $+15\%$ near the $t\bar{t}$ production threshold ( $m_H \approx 2m_t$ ) and become negative for heavier masses, reaching $\approx -20\%$ at $m_H = 1000$ GeV.
- $\tan\beta$ Dependence: In the strict alignment limit, the relative corrections are largely independent of $\tan\beta$ . However, in non-alignment scenarios ( $s_{\beta-\alpha} \neq 1$ ), the corrections can reach up to $+38\%$ for $H \to \gamma\gamma$ at large $\tan\beta$ , and the $H \to Z\gamma$ corrections can cross zero to become positive.
Threshold Behavior and Perturbative Stability: A critical finding is the stability of the fixed-order prediction near the $t\bar{t}$ threshold. Due to CP conservation, the intermediate $t\bar{t}$ pair in a CP-even Higgs decay must be in a $P$ -wave ( $^3P_0$ ) configuration. This suppresses the Coulomb singularity (Sommerfeld enhancement) that typically plagues $S$ -wave processes. The authors demonstrate via threshold expansion that the two-loop QCD correction remains finite, with the imaginary part proportional to $\beta_t^2$ and the real part remaining finite. This confirms that the observed enhancement near threshold is a genuine perturbative effect rather than a breakdown of fixed-order perturbation theory.
Phenomenological Implications: Comparing the NLO predictions for $\sigma(gg \to H) \times \text{Br}(H \to V\gamma)$ with current ATLAS limits (139 fb $^{-1}$ for $\gamma\gamma$ and 140 fb $^{-1}$ for $Z\gamma$ ), the authors find that current data do not yet constrain the Type-I THDM parameter space. The $H \to \gamma\gamma$ predictions reach about 83% of the observed limit near the threshold, while $H \to Z\gamma$ remains roughly 30 times below the sensitivity.

Significance
The paper provides a rigorous, complete numerical evaluation of NLO QCD corrections for heavy Higgs decays in the Type-I THDM, removing previous reliance on the large-top-mass approximation. By establishing the perturbative stability of the prediction near the $t\bar{t}$ threshold through the suppression of Coulomb effects in the $P$ -wave configuration, the work validates the use of fixed-order calculations for interpreting future data. The authors conclude that while current LHC data does not exclude the model, the $\mathcal{O}(15\%)$ NLO corrections will become relevant for interpreting exclusion bounds at the High-Luminosity LHC (HL-LHC), where projected experimental precisions for $H \to \gamma\gamma$ and $H \to Z\gamma$ are expected to reach $\sim 1.8\%$ and $\sim 9.8\%$ , respectively. These results serve as essential input for future experimental searches and precision measurements of heavy neutral Higgs bosons.

More like this