⚛️ phenomenology

Hadron production through Higgs decay at next-to-leading order in the general-mass variable-flavor-number scheme

This paper presents the first next-to-leading order study of Higgs decay into B-mesons within the general-mass variable-flavor-number scheme, demonstrating that incorporating the masses of both b-quarks and B-mesons significantly enhances the partial decay width in the low and peak scaled-energy regions, respectively, unlike previous massless approximations.

Original authors: S. Mohammad Moosavi Nejad

Published 2026-03-20

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: S. Mohammad Moosavi Nejad

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 Higgs boson as a tiny, unstable firework that explodes the moment it's created. When it explodes, it doesn't just vanish; it shoots out a pair of heavy "bottom" quarks (let's call them Bottoms). But these Bottoms are shy and don't stay alone for long. They immediately grab onto other particles to form a "B-meson" (a type of heavy particle), kind of like a Bottom putting on a heavy winter coat before it can even run away.

Scientists at the Large Hadron Collider (LHC) want to study these fireworks to understand the rules of the universe. To do this, they measure how much energy the resulting B-mesons carry. This is like measuring how fast the shrapnel flies after an explosion.

The Old Way: The "Weightless" Guess

For years, scientists calculated these energy patterns using a simplified model. They assumed the Bottom quarks and the B-mesons had zero weight (like ghosts).

The Analogy: Imagine trying to predict the flight path of a bowling ball by pretending it's a feather. It's easy to do the math, but the result is wrong because a feather floats differently than a bowling ball.
The Problem: Because Bottom quarks are actually quite heavy, ignoring their weight meant the old calculations missed important details, especially at the very low and very high ends of the energy spectrum.

The New Way: The "Real-World" Calculation

This paper introduces a new, more sophisticated method called the GM-VFN scheme. Think of this as upgrading from a sketch to a high-definition 3D simulation.

The Analogy: Instead of pretending the bowling ball is a feather, this new method weighs the ball, measures its size, and calculates exactly how the air resistance affects it. It acknowledges that the "coat" (the B-meson) has weight and that the "person inside" (the Bottom quark) also has weight.

What Did They Find?

By doing the math with the actual weights included, the authors discovered two major surprises:

The "Heavy Coat" Effect (B-meson mass):
- What happened: When they included the weight of the B-meson, the number of particles with low energy shot up significantly.
- The Metaphor: Imagine a runner wearing a heavy backpack. If they try to sprint, they slow down. Similarly, because the B-meson is heavy, it "drags" the energy down, creating a much larger crowd of slow-moving particles than the old "feather" model predicted. This creates a "threshold"—a minimum speed limit below which these particles simply cannot exist.
The "Heavy Core" Effect (Bottom quark mass):
- What happened: When they included the weight of the Bottom quark itself, the number of particles with medium-to-high energy (the peak of the curve) increased.
- The Metaphor: Think of a heavy engine in a car. A heavier engine changes how the car accelerates and where the peak speed occurs. The heavy quark pushes the "sweet spot" of the energy distribution to a higher level.

Why Does This Matter?

The LHC is about to get much more powerful, producing billions of these Higgs explosions.

The Goal: Scientists want to know if the Higgs boson behaves exactly as the Standard Model (our current rulebook) predicts, or if there are "glitches" that hint at new physics.
The Impact: If we use the old "feather" math, we might think we see a glitch when it's actually just a calculation error. By using this new "heavy" math, scientists can strip away the confusion. They can now say with much higher confidence: "Yes, this is exactly how nature works," or "Wait, this doesn't match our heavy-weight calculation either—maybe we found something new!"

In a Nutshell

This paper is like updating the instruction manual for a complex machine. The old manual said, "Ignore the weight of the parts to make the math easier." This new manual says, "Actually, the weight matters a lot. If you include it, you'll see that the machine behaves differently at the slow end and the fast end than we thought." This update is crucial for the next generation of experiments to ensure they don't miss any hidden secrets of the universe.

1. Problem Statement

The Higgs boson, discovered at the CERN-LHC, decays into a pair of bottom quarks ( $H \to b\bar{b}$ ) with a branching ratio of approximately 60%. These bottom quarks subsequently hadronize into bottom-flavored mesons ( $B$ -mesons) before decaying. Studying the scaled-energy distribution ( $x_B$ ) of these $B$ -mesons provides a crucial channel for probing Higgs characteristics and measuring Yukawa couplings.

The Core Issue: Previous theoretical studies of this process ( $H \to B + \text{Jets}$ ) have predominantly utilized the Zero-Mass Variable-Flavor-Number Scheme (ZM-VFNS). In this scheme, the masses of both the bottom quark ( $m_b$ ) and the produced $B$ -meson ( $m_B$ ) are set to zero ( $m_b = m_B = 0$ ) to simplify calculations.

Limitation: This approximation ignores significant mass effects, particularly near the kinematic threshold and in the peak region of the energy spectrum.
Goal: To provide a more precise theoretical prediction by retaining the finite mass effects of both the $b$ -quark and the $B$ -meson, which is essential for high-precision analyses at the High-Luminosity LHC (HL-LHC) and future colliders (FCC-ee, ILC).

2. Methodology

The author employs the General-Mass Variable-Flavor-Number Scheme (GM-VFNS), which combines the advantages of the Fixed-Flavor-Number Scheme (FFNS) and the ZM-VFNS while avoiding their respective flaws.

A. Theoretical Framework

Factorization: The calculation relies on the QCD-improved parton model factorization theorem (Collins' theorem). The differential decay rate is expressed as a convolution of hard scattering partonic rates and non-perturbative fragmentation functions (FFs).
Mass Treatment: Unlike ZM-VFNS, GM-VFNS retains the full dependence on $m_b$ and $m_B$ . It resums large logarithms ( $\ln(m_H^2/m_b^2)$ ) while preserving finite mass terms.
Process: The study focuses on the inclusive process $H \to b\bar{b} \to B + X$ at Next-to-Leading Order (NLO) in perturbative QCD (pQCD).

B. Calculation Steps

Parton-Level Calculation (FFNS):
- Calculated the differential decay widths $d\tilde{\Gamma}/dx_i$ for the process $H \to b\bar{b} (+g)$ at NLO.
- Included Born level, Virtual radiative corrections (vertex corrections, mass/wave-function renormalization), and Real radiative corrections (gluon bremsstrahlung).
- Used dimensional regularization ( $D=4-2\epsilon$ ) to handle ultraviolet (UV) and infrared (IR) divergences.
- Derived analytical expressions for the differential rates normalized to the Born rate.
Transition to GM-VFNS:
- Constructed subtraction terms by taking the limit $m_b \to 0$ of the FFNS results and subtracting the known ZM-VFNS results.
- The final GM-VFNS rate is obtained by:
  $\left(\frac{d\Gamma}{dx}\right)_{\text{GM-VFNS}} = \left(\frac{d\tilde{\Gamma}}{dx}\right)_{\text{FFNS}} - \left(\frac{d\Gamma}{dx}\right)_{\text{Subtraction}}$
- This ensures the recovery of ZM-VFNS results in the massless limit while retaining finite mass effects.
Hadron-Level Convolution:
- Applied an improved factorization theorem to account for the finite mass of the $B$ -meson ( $m_B$ ).
- The relation between the partonic scaling variable ( $x_a$ ) and the hadronic scaling variable ( $x_B$ ) was modified to include mass terms:
  $z = \frac{x_B + \sqrt{x_B^2 - \rho_B^2}}{x_a + \sqrt{x_a^2 - \rho_a^2}}$
  where $\rho_i = 2m_i/m_H$ .
- Convolved the partonic rates with non-perturbative fragmentation functions ( $D_B^a(z, \mu)$ ) extracted from $e^+e^-$ annihilation data (ALEPH, OPAL, SLD).

C. Numerical Inputs

Parameters: $m_H = 125.3$ GeV, $m_b = 4.78$ GeV, $m_B = 5.279$ GeV.
Coupling: $\alpha_s$ calculated at NLO with $\Lambda_{\overline{MS}}^{(5)} = 231.0$ MeV.
Fragmentation Functions: Used a power-law model ( $N z^\alpha (1-z)^\beta$ ) fitted to experimental data, evolved via DGLAP equations.

3. Key Contributions

First NLO GM-VFNS Calculation: This is the first study to compute the $B$ -meson energy spectrum from Higgs decay at NLO while simultaneously retaining the finite masses of both the $b$ -quark and the $B$ -meson.
Analytical Derivation: Provided new analytical expressions for the NLO differential decay rates in the massive scheme, including the necessary subtraction terms to match ZM-VFNS.
Improved Factorization: Demonstrated the application of an improved factorization theorem that correctly handles the kinematic thresholds introduced by non-zero hadron masses.

4. Key Results

The numerical analysis reveals significant deviations from the massless (ZM-VFNS) predictions:

Threshold Effect: The finite mass of the $B$ -meson ( $m_B$ ) introduces a kinematic threshold at $x_B \approx 0.093$ . The spectrum is strictly zero below this value, a feature absent in massless schemes.
Low- $x_B$ Enhancement: The inclusion of the finite $B$ -meson mass leads to a significant enhancement (approx. 27% at the threshold) of the partial decay width in the low- $x_B$ region ( $0.093 \le x_B \le 0.25$ ).
Peak and High- $x_B$ Enhancement: The finite $b$ -quark mass ( $m_b$ ) is responsible for enhancements in the peak region ( $x_B \approx 0.8$ ) and the high- $x_B$ tail. Specifically, NLO corrections increase the decay width in the peak region by up to 46% compared to Leading Order (LO).
Gluon Fragmentation: The contribution from gluon fragmentation ( $g \to B$ ) is negative and only appreciable in the low- $x_B$ region ( $x_B < 0.3$ ). At higher $x_B$ , the spectrum is dominated by $b \to B$ fragmentation.
Scale Dependence: The results show considerable sensitivity to the factorization scale $\mu$ in the peak region and above, with variations observed for $\mu$ in the range $[m_H/2, 2m_H]$ .

5. Significance

Precision Physics: The findings are critical for the upcoming High-Luminosity LHC (HL-LHC) and future lepton colliders (FCC-ee, ILC), where statistical uncertainties will be drastically reduced, making theoretical precision the limiting factor.
Higgs Characterization: Accurate modeling of the $B$ -meson energy spectrum is essential for extracting precise Yukawa couplings and testing the Standard Model. Ignoring mass effects could lead to misinterpretation of experimental data.
Theoretical Robustness: The work validates the GM-VFNS as a superior framework for heavy hadron production, bridging the gap between fixed-flavor calculations (valid at low scales) and massless approximations (valid at very high scales).
General Applicability: The formalism developed is not limited to $B$ -mesons but can be extended to other heavy hadrons (e.g., $D$ -mesons, protons) produced in Higgs decays.

In conclusion, this paper establishes that mass effects are not negligible corrections but dominant features in specific kinematic regions of Higgs decay. The proposed GM-VFNS framework provides the necessary theoretical tool to exploit the high-statistics data expected from future collider runs.