Fully charm tetraquark production at hadronic collisions with gluon radiation effects

This paper presents the first complete next-to-leading order QCD calculation for fully charm tetraquark production at hadronic collisions, incorporating gluon radiation effects through resummation and utilizing LHCb and CMS data to extract universal long-distance matrix elements and predict kinematic distributions for the $X(6900)$ and its spin-zero partners.

The Gist

The Big Picture: Building a "Super-Heavy" Lego Castle

Imagine the universe is a giant construction site. Usually, the builders (particles) stick to simple rules: they build small houses (protons) or tiny cars (mesons) using just two or three bricks (quarks).

For a long time, physicists thought it was impossible to build a house out of four heavy bricks (charm quarks) stuck together. But in 2020, the LHCb experiment at CERN found a strange new structure called X(6900). It's a "tetraquark"—a four-brick castle made entirely of heavy charm quarks.

This paper is like a team of master architects (Wang and Zhu) who just finished the most detailed blueprint ever created for how to build this specific castle. They didn't just guess; they used the most advanced math available to predict exactly how often these castles appear when you smash particles together, and what they look like when they break apart.

---

1. The Blueprint: How the Bricks Fit Together

The authors had to figure out the internal structure of this four-brick castle. They asked: How are these four heavy bricks holding hands?

• The "Good" vs. "Bad" Handshake: Imagine two pairs of twins.
  • Scenario A: The twins hold hands tightly in a specific, symmetrical way (like a dance). This is the "color-symmetry-antisymmetry" basis. It's like a very stable, compact knot.
  • Scenario B: The twins form two separate couples that just happen to be standing near each other, maybe holding a balloon between them. This is the "color-singlet-octet" basis. It's like two separate cars parked next to each other.
• The Discovery: The authors calculated that the universe prefers the "tight knot" version for this specific particle. They proved that the math describing how these bricks stick together is surprisingly simple (the "renormalization constant is exactly unity"), meaning nature isn't trying to trick us with hidden complexity here.

2. The Construction Site: Smashing Particles

To build these castles, you need a massive construction site. In this case, it's the Large Hadron Collider (LHC), where they smash protons together at near light speed.

• The Glue (Gluons): When protons smash, they don't just hit head-on; they spray out a lot of "glue" (gluons).
• The Problem: If you try to calculate the result of the crash using simple math, the "glue" creates a mess. It's like trying to count the raindrops in a hurricane; the numbers get huge and nonsensical because of the "soft" and "collinear" (gliding) gluons.
• The Solution (Resummation): The authors used a special mathematical technique called Resummation. Think of this as a noise-canceling headphone for the math. It filters out the chaotic, infinite noise of the extra gluons and sums them up into a clean, predictable signal. This allowed them to predict the results with extreme precision, even when the particles are moving slowly or very fast.

3. The "Fingerprint" of the Castle

The authors didn't just calculate the math; they compared it to real-world photos taken by the LHCb and CMS experiments.

• The Match: They looked at how often the X(6900) appears compared to a standard pair of J/ψ particles (which are like two-brick cars).
• The Result: By matching their complex calculations to the real data, they were able to extract a "fingerprint" called the Long-Distance Matrix Element (LDME).
  • Analogy: Imagine you see a car driving down the street. You can't see the engine, but by watching how fast it accelerates and how it handles turns, you can deduce exactly how powerful the engine is.
  • The Paper's Achievement: They deduced the "engine power" (the internal binding strength) of the X(6900) tetraquark. This is a universal number that can now be used to predict how this particle behaves in any experiment, not just the one they studied.

4. Predicting the Future: What to Look For

Now that they have the blueprint and the fingerprint, they made predictions for what other experiments should see:

• Speed and Direction: They predicted exactly how many of these castles will be made if you look at them moving at different speeds (transverse momentum) or at different angles (rapidity).
• The Spin: They confirmed that the X(6900) is likely a "spin-2" object (like a spinning top that wobbles in a specific way), not a "spin-0" ball.
• The Partners: They also predicted that there should be a "sibling" particle (a spin-0 version) that is harder to find but should exist. They provided a map for where to look for it.

Why This Matters

Before this paper, we knew the X(6900) existed, but we didn't fully understand how it was built or how it was made.

• The "First" Achievement: This is the first time anyone has done a "Next-to-Leading Order" (NLO) calculation for this specific particle. In the world of physics, this is like moving from a sketch on a napkin to a full 3D architectural model with stress tests.
• The Impact: It confirms that our understanding of the "strong force" (the glue holding atoms together) is correct, even for these exotic, heavy four-quark states. It gives experimentalists a precise target: "Look here, at this speed, and you will find the particle."

In short: These physicists built the ultimate math model for a mysterious four-quark particle, proved it matches real-world photos, and gave the rest of the scientific community a treasure map to find more of them.

Technical Summary

1. Problem Statement

The discovery of the $X(6900)$ state by the LHCb experiment in 2020 provided the first evidence of a fully charm tetraquark ($T_{4c}$, composed of $cc\bar{c}\bar{c}$). However, several critical theoretical challenges remained:

• Theoretical Precision: Previous calculations were often limited to Leading Order (LO) or lacked a complete treatment of QCD corrections, which are crucial for hadronic collisions where gluon radiation is significant.
• Internal Structure: The internal quark configuration (compact diquark-antidiquark vs. molecular charmonium pairs) and the specific spin-parity ($J^{PC}$) quantum numbers of the $X(6900)$ and its potential partners were not fully constrained by theory.
• Kinematic Discrepancies: Experimental data showed signal enhancements at medium-to-high transverse momentum ($P_\perp$) that standard perturbative QCD (pQCD) struggles to describe accurately due to large logarithmic terms ($\ln(P_\perp^2/M^2)$) arising from soft and collinear gluon emissions.
• Non-perturbative Parameters: The Long-Distance Matrix Elements (LDMEs) required for Non-Relativistic QCD (NRQCD) factorization for fully heavy tetraquarks were unknown and needed extraction from experimental data.

2. Methodology

The authors employed a rigorous theoretical framework combining Non-Relativistic QCD (NRQCD) factorization with Transverse Momentum Dependent (TMD) factorization and resummation techniques.

• Basis Expansion: The fully charm tetraquark states were expanded in the color-symmetry-antisymmetry basis ($6 \otimes \bar{6}$ and $\bar{3} \otimes 3$). This allows for the description of both "good" (axial-vector diquark) and "bad" (scalar diquark) configurations, as well as color-singlet molecular-like states.
• Perturbative Calculation (NLO):
  • Calculated the differential and total cross-sections for $pp \to T_{4c} + X$ at Next-to-Leading Order (NLO) in $\alpha_s$.
  • Included all possible spin configurations ($0^{++}$ and $2^{++}$) and mixing effects.
  • Generated and computed 2,008 one-loop diagrams and 618 real emission diagrams for the gluon-gluon fusion channel, and 170/98 diagrams for quark-antiquark annihilation/scattering.
  • Used dimensional regularization ($D=4-2\epsilon$) to handle UV and IR divergences, utilizing packages like FeynArts, FeynCalc, and Kira.
• Resummation (NLL):
  • Applied Transverse Momentum Resummation to Next-to-Leading Logarithmic (NLL) accuracy.
  • This technique resums large logarithms induced by soft and collinear gluon radiation to all orders in $\alpha_s$, specifically addressing the divergence of perturbative results at low $P_\perp$.
  • Utilized the $b$-space formalism with a non-perturbative Sudakov factor ($W_{NP}$) fitted to global data (BLNY and SIYY forms) to handle the impact parameter region $b \geq 1/\Lambda_{QCD}$.
• Renormalization: Investigated the renormalization constants of the color-singlet four-charm-quark operators.

3. Key Contributions

• First Complete NLO Calculation: This is the first complete NLO QCD calculation for fully charm tetraquark production at hadronic colliders.
• Renormalization Constant Discovery: The authors discovered that the renormalization constant of the color-singlet four-charm-quark operator is exactly unity at NLO. This is a non-trivial result not previously explored in literature, simplifying the theoretical framework.
• LDME Extraction: By combining the NLO+NLL theoretical predictions with LHCb data on the total cross-section of $X(6900)$ and CMS measurements of its spin-parity ($2^{++}$), the authors successfully extracted the non-perturbative but universal Long-Distance Matrix Element (LDME) for the $2^{++}$ state in a model-independent way.
• Prediction of Partners: The study provides predictions for the production of the spin-zero ($0^{++}$) partner of the $X(6900)$ and other quark configurations, including their rapidity and transverse momentum distributions.

4. Key Results

• LDME Value: The extracted LDME for the $2^{++}$ state is:
  $$\langle \mathcal{O}_{3\bar{3}}^{T_{4c}} (5S_2, 2^{++}) \rangle \mathcal{B}(T_{4c}) = (2.22 \pm 0.80^{+1.29+0.62}_{-0.57-0.00}) \times 10^{-4} \, \text{GeV}^9$$
  This value was determined using LHCb data ($R \approx 1.1\%$) and validated against CMS spin-parity constraints.
• Cross-Section Ratios: The theoretical prediction for the ratio $R = \sigma(T_{4c})/\sigma(2J/\psi) \times \mathcal{B}$ with a transverse momentum cut $P_\perp > 5.2$ GeV is $R_{theo} = [3.2 \pm 2.0]\%$, which shows good agreement with LHCb data ($R_{LHCb} = [2.6 \pm 0.6 \pm 0.8]\%$).
• Kinematic Distributions:
  • $P_\perp$ Distribution: The NLO+NLL resummation successfully removes the singular structure at $P_\perp \to 0$ seen in fixed-order NLO calculations, providing a smooth, physical distribution.
  • Enhancement: The $2^{++}$ tetraquark cross-section is significantly enhanced compared to the $0^{++}$ state, driven by both differences in LDMEs and short-distance coefficients.
  • Rapidity: The differential cross-section exhibits a relatively flat profile across low and intermediate rapidity ranges.
• Angular Distribution: The polar angular distribution of the $J/\psi$ pairs in the decay of the $2^{++}$ state was fitted to CMS data, yielding coefficients ($L_1 \approx 170, L_3 \approx -59$, etc.) that support the $J^{PC}=2^{++}$ assignment.

5. Significance

• Validation of QCD Factorization: The work demonstrates that the combination of NRQCD and TMD resummation is a robust framework for describing heavy tetraquark production, even for states with four heavy quarks.
• Guidance for Future Experiments: The predictions for the $0^{++}$ partner and the differential distributions (rapidity and $P_\perp$) provide specific targets for the LHCb, ATLAS, and CMS collaborations to verify the existence and nature of the tetraquark family.
• Understanding Color Confinement: By distinguishing between compact diquark-antidiquark structures and molecular states via the color basis expansion and comparing with data, the paper offers insights into the color confinement mechanism for multi-heavy quark systems.
• Foundation for Spectroscopy: The methodology established here serves as a foundation for systematically studying the production and inner structure of the entire family of fully-charmed tetraquark states (including P-wave and higher excitations).

In summary, this paper bridges the gap between experimental observations of the $X(6900)$ and fundamental QCD theory by providing the first high-precision (NLO+NLL) calculation, extracting crucial non-perturbative parameters, and predicting the properties of yet-to-be-observed tetraquark partners.