⚛️ phenomenology

Interpreting the top quark mass using simulations of top-flavoured mesons

This paper utilizes Pythia 8.3 simulations of fictitious top-flavoured mesons within Heavy Quark Effective Theory to determine that the top quark mass parameter in the generator is approximately 240–330 MeV lower than the top quark pole mass, a difference consistent with non-perturbative QCD effects.

Original authors: Alexander Lind

Published 2026-10-23

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Alexander Lind

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 subatomic world, the top quark holds a unique place. It is the heaviest known elementary particle, weighing as much as a gold atom, yet it exists for only a fleeting instant before decaying into other particles. Because it is so heavy, its mass is a critical number for physicists. It helps determine whether the universe is stable or if it might one day collapse, and it influences how the Higgs field, which gives particles their mass, behaves. However, measuring this mass is tricky. Since the top quark vanishes before it can be directly weighed on a scale, scientists must infer its mass by watching the debris it leaves behind in particle collisions. They compare the patterns of this debris against complex computer simulations to find the best match. The number that emerges from this process is called the "Monte Carlo mass." For years, a quiet debate has simmered among theorists: does this computer-generated number represent the true, fundamental mass of the top quark, or is it slightly skewed by the messy, unpredictable forces that govern how particles stick together?

A recent study by Alexander Lind at the INFN in Rome tackles this question by taking a bold, imaginative step. Since the top quark decays too quickly to form stable particles in nature, the researcher created a fictional scenario within a computer program called Pythia. In this simulation, he forced the top quarks to behave as if they were stable enough to bind with lighter quarks, forming heavy, short-lived particles known as top-flavored mesons. By treating these made-up particles as if they were real, Lind could apply a well-established set of rules called Heavy Quark Effective Theory. These rules describe how the mass of a heavy particle relates to the mass of the larger particle it helps form. In the real world, this relationship is used to understand particles like the B-meson, but here it was applied to the hypothetical top meson to create a bridge between the computer's input settings and the theoretical "pole mass," which is the standard definition of a particle's true mass.

The simulation allowed the researcher to run thousands of virtual collisions, mimicking the conditions of both the Large Hadron Collider at CERN and a proposed future lepton collider. In these virtual events, the top quarks were made to decay into the fictional mesons, which then broke apart in a predictable way. The team focused on a specific measurement: the combined weight of the resulting B-meson and a lepton, a type of particle similar to an electron. By adjusting the input mass of the top quark in the computer program and watching how this specific weight changed, they could trace a direct line between the simulation's settings and the theoretical mass. The results revealed a consistent gap. The mass value fed into the computer program was systematically lower than the theoretical pole mass. Specifically, the simulations showed that the computer mass is about 240 to 330 million electronvolts lower than the pole mass. This difference is small in the grand scheme of particle physics, but it is significant enough to be measured, and it aligns with the scale of the forces that bind quarks together.

This finding challenges the idea that the mass extracted from current experiments is identical to the fundamental theoretical mass. Instead, the study suggests that the standard way of measuring the top quark mass in particle colliders carries a built-in offset caused by the complex dynamics of how particles form and decay. The research does not claim to have solved the mystery of the top quark's mass, but it provides a concrete estimate of the shift introduced by the modeling of these non-perturbative effects. The work remains a simulation, a controlled experiment within a computer, but it offers a clear path forward. By understanding this offset, physicists can better interpret the data from real-world collisions, ensuring that the fundamental numbers they use to describe the universe are as precise as possible. The study concludes that the difference is real and measurable, consistent with the expected scale of quantum chromodynamics, the theory describing the strong nuclear force.

Technical Summary: Interpreting the Top Quark Mass via Top-Flavoured Meson Simulations

Problem Statement
The precise determination of the top quark mass ( $m_t$ ) is critical for Standard Model (SM) consistency checks, including electroweak precision tests, vacuum stability, and Higgs inflation models. Currently, the top mass is extracted from collider data (e.g., $t\bar{t}$ kinematics) by comparing experimental observables to Monte Carlo (MC) event generators like Pythia or Herwig. This yields the "Monte Carlo mass" ( $m_t^{\text{MC}}$ ). A central theoretical challenge is relating $m_t^{\text{MC}}$ to well-defined theoretical mass schemes, such as the pole mass ( $m_t^{\text{pole}}$ ) or the $\overline{\text{MS}}$ mass.

The relationship is generally expressed as $m_t^{\text{MC}} = m_t^{\text{pole}} + \delta m \pm \Delta m$ , where $\delta m$ represents a potential theoretical shift and $\Delta m$ represents uncertainty from non-perturbative QCD modeling (hadronisation). Literature presents two competing views:

$\delta m \simeq 0$ , implying $m_t^{\text{MC}}$ is effectively the pole mass within an uncertainty of $\mathcal{O}(\Lambda_{\text{QCD}})$ .
$\delta m \neq 0$ , suggesting a fundamental difference with shifts estimated between 200 MeV and 1 GeV.

Methodology
To investigate the size of non-perturbative effects on mass extraction, the author proposes a simulation framework where top quarks are forced to hadronise into fictitious heavy-light top-flavoured mesons ( $T$ -mesons, e.g., $T^+ = t\bar{d}$ ), despite the physical top quark decaying before hadronisation.

Simulation Setup: Using a modified version of Pythia 8.3, the author simulates $t\bar{t}$ production at the LHC ( $\sqrt{s} = 13.6$ TeV) and a future lepton collider (FLC, $\sqrt{s} = 1$ TeV). The top width is set to zero (narrow-width approximation) to allow the top quark to form a $T$ -meson.
Decay Model: The $T$ -mesons decay via a spectator model analogous to $B$ -mesons. The top quark inside the meson decays ( $t \to Wb \to \ell\nu b$ ), while the light spectator quark carries a fraction of the meson's momentum. The resulting $b$ -quark and spectator undergo standard hadronisation.
Theoretical Connection: Heavy Quark Effective Theory (HQET) is used to relate the mass of the fictitious $T$ -meson ( $m_T$ ) to the top pole mass ( $m_t^{\text{pole}}$ ). In the limit of large top mass, this relation is approximated as a constant shift: $m_T \approx m_t^{\text{pole}} + \bar{\Lambda}$ , where $\bar{\Lambda} \approx 0.473$ GeV is derived from $B$ -meson parameters.
Observable: The analysis focuses on the invariant mass of the $B$ -hadron and the lepton ( $m_{B\ell}$ ). The author varies the input mass parameters in Pythia ( $m_t^{\text{Pythia}}$ for standard samples and $m_T$ for $T$ -meson samples) and fits the resulting $\langle m_{B\ell} \rangle$ distributions.

Key Results
The study establishes a linear dependence of the average observable $\langle m_{B\ell} \rangle$ on the input mass parameters:

For standard samples: $\langle m_{B\ell} \rangle \propto m_t^{\text{Pythia}}$
For $T$ -meson samples: $\langle m_{B\ell} \rangle \propto m_T$

By equating the slopes of these linear fits (which agree within errors) and substituting the HQET relation ( $m_T \approx m_t^{\text{pole}} + \bar{\Lambda}$ ), the author extracts the difference between the input Pythia mass and the pole mass:

LHC ($pp$ collisions): $m_t^{\text{Pythia}} - m_t^{\text{pole}} = -0.24 \pm 0.07$ GeV.
FLC ( $e^+e^-$ collisions): $m_t^{\text{Pythia}} - m_t^{\text{pole}} = -0.33 \pm 0.07$ GeV.

These results indicate a negative mass difference of approximately 240–330 MeV, consistent with an $\mathcal{O}(\Lambda_{\text{QCD}})$ shift. The study also examined the sensitivity to the "recoil treatment" in top decays (whether the $W$ boson or the $b$ -quark absorbs the recoil). While this choice affects the extracted mass difference (by $\sim 40$ MeV for $pp$ and $\sim 120$ MeV for $e^+e^-$ ), the effect remains within the statistical uncertainties of the simulation.

Significance and Claims
The paper claims to provide a quantitative estimate of the systematic uncertainty introduced by non-perturbative modelling in the interpretation of the top quark mass. By simulating a scenario where top quarks hadronise, the author isolates the impact of hadronisation dynamics on mass extraction.

The primary conclusion is that the Monte Carlo mass extracted from current and future collider data is likely lower than the pole mass by a value in the range of 240–330 MeV. This finding supports the view that the Monte Carlo mass is not identical to the pole mass but differs by a shift of order $\Lambda_{\text{QCD}}$ , rather than being a distinct parameter with a large, uncontrolled shift. The work serves as a proof-of-concept for using fictitious heavy-light meson simulations to constrain theoretical mass definitions, with ongoing work planned to extend these simulations to Herwig 7.3.0.

Technical Summary: Interpreting the Top Quark Mass via Top-Flavoured Meson Simulations

More like this