Mass and Decay-Constant Evolution of Heavy Quarkonia and $B_c$ States from Thermal QCD Sum Rules

This paper utilizes finite-temperature QCD sum rules with updated quark masses and lattice-informed condensates to predict the thermal evolution of masses and decay constants for $J/\psi$, $\Upsilon$, and $B_c$ mesons up to near-critical temperatures, revealing a suppression hierarchy consistent with binding energies and confirming the $B_c$ $1P$--$1S$ splitting observed by LHCb.

The Gist

Imagine the universe as a giant, invisible soup. Under normal conditions, this soup is cool and clear. But inside a particle accelerator like the Large Hadron Collider (LHC), scientists smash atoms together so hard that they create a tiny, fleeting drop of this soup that is hotter than the center of the sun. This is called the Quark-Gluon Plasma.

In this super-hot soup, the rules of how particles stick together change. The paper you asked about is a detailed recipe for predicting how specific "heavy" particles behave as they get cooked in this soup.

Here is the breakdown using simple analogies:

1. The Characters: Heavy Quark "Balls"

Inside atoms, there are tiny particles called quarks. Usually, they are glued together in pairs to form particles called mesons.

• $J/\psi$ and $\Upsilon$ (Upsilon): These are like two identical twins holding hands. One is a "charm" twin, the other is a "bottom" twin. They are very heavy and hold on tight.
• $B_c$: This is a mixed pair: one "charm" and one "bottom." Because they are different sizes, they hold on a bit differently.

The paper looks at three specific types of these pairs:

1. $J/\psi$ (Charm-Charm)
2. $\Upsilon$ (Bottom-Bottom)
3. $B_c$ (Charm-Bottom)

2. The Method: The "Thermal QCD Sum Rule"

How do scientists study these particles when they are melting in a super-hot soup? They can't just look at them with a microscope; the soup is too chaotic.

Instead, the authors use a mathematical tool called QCD Sum Rules. Think of this like predicting the weather.

• You can't see the future, but you can look at current data (pressure, humidity, wind) and use a complex formula to predict if it will rain.
• In this paper, the "formula" is the Sum Rule. It takes known facts about the particles (their mass, how they are built) and adds in the "heat" of the soup to predict what happens to them as the temperature rises.

The authors updated their "formula" with the latest data (like a weather app updated with 2024 satellite data) to make their predictions more accurate.

3. The Experiment: Heating the Soup

The scientists simulated heating this soup from room temperature up to the "boiling point" (called the Critical Temperature, or $T_c$). At this boiling point, the soup changes phase, and the particles that were glued together might fall apart.

They tracked two things as the heat rose:

• Mass ($m$): How heavy the particle feels.
• Decay Constant ($f$): How "strongly" the particle holds together (its glue strength).

4. The Results: Who Melts First?

The most exciting part of the paper is the hierarchy of melting. Imagine three different types of ice cubes in a hot pan:

• A thick block of granite ice ($\Upsilon$): It barely changes. Even when the pan is hot, it stays solid.
  • Result: The $\Upsilon$ particle is the most stable. It barely loses its "glue" even near the boiling point.
• A standard ice cube ($J/\psi$): It starts to sweat and shrink a bit, but it holds together for a while.
  • Result: The $J/\psi$ gets weaker as it heats up, but it survives longer than the mixed pair.
• A delicate snowflake ($B_c$): This one melts the fastest. Because the two parts are different sizes, they don't hold on as tightly in the heat.
  • Result: The $B_c$ particle loses its strength the quickest and is the first to "melt" or dissolve into the soup.

The Order of Melting:
$\Upsilon$ (Survives best) > $J/\psi$ (Middle) > $B_c$ (Melts first)

5. Why This Matters

• Checking the Theory: The authors used this method to predict the mass of a recently discovered excited version of the $B_c$ particle. Their prediction matched the real-world data from the LHCb experiment almost perfectly. This proves their "recipe" works.
• Understanding the Early Universe: Just after the Big Bang, the entire universe was this hot soup. Understanding which particles melt first helps us understand how the universe cooled down and formed the atoms we have today.
• Heavy Ion Collisions: When scientists smash heavy atoms together to recreate the Big Bang, they see these particles disappearing. This paper gives them a baseline to understand why they are disappearing and when.

Summary

This paper is like a survival guide for heavy particles in a heatwave. It tells us that if you throw a heavy particle into the hottest soup imaginable:

1. The heavy, identical twins ($\Upsilon$) will survive the longest.
2. The mixed pair ($B_c$) will fall apart first.
3. The authors updated their math to match the latest experimental data, confirming that their predictions are reliable up to about 90% of the boiling point.

It's a bridge between the abstract math of the universe's fundamental forces and the real-world data we see in particle accelerators today.

Technical Summary

1. Problem Statement and Motivation

The paper addresses the need for a precise, calibrated theoretical framework to describe the behavior of heavy quarkonia ($J/\psi$, $\Upsilon$) and the mixed-flavor $B_c$ meson in a thermal QCD medium (finite temperature $T$). Previous analyses (e.g., from 2016) faced limitations due to:

• Outdated input parameters (quark masses, condensates).
• Ad-hoc modeling of the temperature-dependent gluon condensate and continuum threshold.
• Lack of direct comparison with recent experimental data, specifically the 2025 LHCb observation of orbitally excited $B_c$ states.
• Uncertainties regarding the validity of the Operator Product Expansion (OPE) near the critical temperature ($T_c$).

The primary goal is to re-evaluate the thermal evolution of masses $m(T)$ and decay constants $f(T)$ up to $T/T_c \lesssim 0.9$, using modern inputs and a rigorous uncertainty analysis to establish a reliable baseline for future higher-order corrections.

2. Methodology

The study utilizes the Finite-Temperature QCD Sum Rule (TQCDSR) formalism. The methodology involves several key technical steps:

• Correlator and Currents: The analysis focuses on the thermal two-point correlator for vector ($J/\psi, \Upsilon, B_c^{(1S)}$) and axial-vector ($B_c^{(1P)}$) currents. The ground-state pole contribution is parameterized by mass and decay constant.
• OPE and Thermal Inputs:
  • Perturbative Side: Uses Leading Order (LO) spectral density.
  • Non-Perturbative Side: Includes dimension-4 ($D=4$) gluon condensate terms. Higher-dimensional operators ($D=6$) are neglected but estimated to be subleading ($<10\%$) within the working Borel windows.
  • Gluon Condensate: The temperature dependence of the gluon condensate $\langle \alpha_s G^2 \rangle_T$ is anchored to modern HotQCD lattice QCD determinations of the trace anomaly and equation of state, replacing previous phenomenological fits.
  • Continuum Threshold ($s_0$): Modeled as $s_0(T) = s_0(0)[1 - (T/T_c)^{n_C}] + (m_1+m_2)^2(T/T_c)^{n_C}$. The exponent $n_C$ is channel-dependent to reflect different binding strengths.
• Calibration at $T=0$:
  • Uses PDG 2024 quark masses.
  • Calibrates the vacuum continuum threshold $s_0(0)$ and decay constants $f(0)$ to reproduce experimental masses (PDG/LHCb) and lattice reference values.
  • This ensures that subsequent finite-temperature predictions are model-based extrapolations from a verified vacuum baseline, rather than parameter-free ab initio calculations.
• Borel Stability: The analysis is restricted to Borel windows where pole dominance ($P_{pole} \gtrsim 0.6$) and OPE convergence ($F_{D=4} \lesssim 0.3$) are satisfied. The validity range is strictly limited to $T < 0.9 T_c$ to avoid deconfinement effects that invalidate the single-condensate approximation.

3. Key Contributions

• Updated Inputs: First application of PDG 2024 masses and lattice-informed gluon condensate evolution to heavy quarkonia sum rules.
• $B_c$ Spectroscopy: Direct comparison with the LHCb 2025 observation of excited $B_c$ states. The model successfully predicts the $1P-1S$ mass splitting ($\Delta m \approx 0.477$ GeV), which aligns with the experimental range of $0.430–0.478$ GeV.
• Systematic Uncertainty Quantification: A comprehensive error budget is provided, identifying five main sources: Borel window variation, continuum threshold modeling, gluon condensate normalization, heavy-quark mass errors, and perturbative truncation. The total systematic uncertainty is estimated at $\sim 10-12\%$.
• Refined Thermal Hierarchy: The study clarifies the sequential melting pattern, correcting previous "uniform shift" criticisms by demonstrating distinct suppression rates based on binding energy.

4. Key Results

A. Zero-Temperature Calibration

The model reproduces experimental benchmarks with high precision (within $\sim 0.6\%$ for masses and $\sim 2\%$ for decay constants), validating the calibration procedure:

• $m_{J/\psi}(0) = 3.103$ GeV (Exp: 3.097 GeV)
• $m_{\Upsilon}(0) = 9.517$ GeV (Exp: 9.460 GeV)
• $m_{B_c(1S)}(0) = 6.239$ GeV (Exp: 6.275 GeV)
• $m_{B_c(1P)}(0) = 6.716$ GeV (Consistent with LHCb peaks at 6.705 and 6.752 GeV).

B. Thermal Evolution ($T \to 0.9 T_c$)

The study reveals a clear hierarchy of thermal suppression driven by binding energies ($\epsilon_{bind}$):

1. $\Upsilon$ (Bottomonium): Most stable.
   • Mass suppression: $\Delta m/m \approx -0.5\%$.
   • Decay constant survival: $f(0.9T_c)/f(0) \approx 0.79$.
   • Dissociation temperature: $T_{diss} > 0.9 T_c$.
2. $J/\psi$ (Charmonium): Moderate suppression.
   • Mass suppression: $\Delta m/m \approx -6.4\%$.
   • Decay constant survival: $f(0.9T_c)/f(0) \approx 0.75$.
   • Dissociation temperature: $T_{diss} \approx 0.87 T_c$.
3. $B_c$ (Mixed Charm-Bottom): Strongest suppression.
   • Vector ground state survival: $f(0.9T_c)/f(0) \approx 0.64$.
   • Axial-vector ($1P$) survival: $f(0.9T_c)/f(0) \approx 0.54$.
   • Dissociation temperature: $T_{diss} \approx 0.80 T_c$ (earlier than $J/\psi$).

C. Physical Interpretation

The sequential melting is attributed to color screening in the thermal medium. The $B_c$ system, having unequal quark masses, possesses a larger spatial extent (Bohr radius) and reduced reduced mass compared to symmetric quarkonia, making it more susceptible to the drop in the non-perturbative gluon condensate $\langle \alpha_s G^2 \rangle_T$.

5. Significance and Future Directions

• Theoretical Consistency: The paper resolves previous methodological concerns by strictly limiting the analysis to the confined phase ($T < T_c$) and anchoring inputs to lattice QCD.
• Experimental Relevance: The results provide a coherent baseline for interpreting heavy-ion collision data (e.g., ALICE and LHCb nuclear modification factors $R_{AA}$), where the early melting of $B_c$ states is consistent with observed suppression patterns.
• Foundation for Extensions: The calibrated LO+$D=4$ framework serves as a reference for future work, which the author plans to extend by including:
  • Radiative $O(\alpha_s)$ corrections.
  • Finite-width effects (Breit-Wigner smearing).
  • Higher-dimensional ($D=6$) condensates.
  • Running coupling $\alpha_s(T)$.

In conclusion, this work establishes a robust, calibrated phenomenological model for heavy quarkonia in thermal media, successfully bridging lattice thermodynamics, QCD sum rules, and recent experimental spectroscopy.