⚛️ nuclear theory

Heavy quark coalescence probability in the presence of a potential

This study develops a phenomenological heavy quark potential based on the constituent quark model to demonstrate that the heavy quark coalescence probability remains close to unity at low momentum, provided that medium-induced modifications in the quark-gluon plasma are sufficiently moderate.

Original authors: Taesoo Song, Jiaxing Zhao

Published 2026-07-10

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Taesoo Song, Jiaxing Zhao

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 a giant, super-hot soup called the Quark-Gluon Plasma (QGP). It's the stuff the universe was made of just after the Big Bang, and scientists recreate it by smashing heavy atoms together at nearly the speed of light. Inside this soup, there are "heavy" particles called heavy quarks (like charm and bottom quarks) and a sea of lighter particles.

As the soup cools down, these heavy quarks need to grab a partner from the light particles to form a new "family" called a heavy meson. This process is called coalescence. Think of it like a heavy dancer in a crowded room trying to find a light partner to dance with.

The Big Problem: The "Empty Room" Theory
For a long time, scientists tried to calculate how likely it is for a heavy quark to find a partner. They assumed the light particles were spread out evenly, like people randomly scattered across a huge dance floor. When they did the math using this "even spread" idea, the results were disappointing. The calculation suggested that even at the very end, when the heavy quark is standing still, it would only find a partner about 60% to 70% of the time.

This didn't make sense. If a heavy quark is just sitting there, it should be able to grab a partner 100% of the time. It's like saying a dancer standing still in a room full of people can't find anyone to dance with. To fix this, previous models had to cheat by adding a "magic number" to force the answer to be 100%. But that felt like a cheat code, not real physics.

The New Idea: The Invisible Magnet
In this study, Taesoo Song and Jiaxing Zhao asked a simple question: What if the heavy quark isn't just sitting in an empty room, but is actually pulling the light particles toward it?

They introduced a heavy-quark potential, which you can think of as an invisible magnet or a gravity well. Just like a planet pulls in nearby asteroids, a heavy quark pulls in the light anti-quarks around it. This means the light particles aren't spread out evenly; they are crowded tightly around the heavy quark, waiting to be picked up.

The Results: Finding the Partner
When the scientists ran their simulations with this "invisible magnet" included, the story changed completely.

Without the magnet: The chance of finding a partner was low (far below 100%).
With the magnet: The chance of finding a partner jumped up to nearly 100% (unity).

This happened without needing any "magic numbers" or cheat codes. The attraction of the potential naturally packed the light particles close enough to the heavy quark that coalescence became almost guaranteed.

What Happens in the Hot Soup? (The Screening Effect)
The scientists also wondered what happens if the "soup" gets really hot and dense. In the QGP, the heat can weaken that invisible magnet. They call this color screening. Imagine the magnet getting rusty or the soup getting so thick it blocks the magnetic pull.

They simulated what happens if the magnet gets weaker (by increasing a parameter called the screening mass, $\mu$ ).

They found that if the magnet gets too weak, the light particles scatter away again, and the chance of finding a partner drops below 100%.
However, as long as the magnet stays strong enough (specifically, if the screening parameter $\mu/\sqrt{\sigma}$ stays below about 0.5), the coalescence probability stays close to 100%.

The Mass Shift
There's one more cool side effect. When the magnet gets weaker due to the heat, the heavy mesons they form get slightly lighter. The study suggests that near the temperature where the soup turns back into normal matter (the crossover temperature, $T_c$ ), the mass of these heavy mesons might drop by about 60 MeV compared to their weight in a vacuum.

The Bottom Line
The paper suggests that the reason heavy quarks successfully form particles is because they act like magnets, gathering their partners close by. If the heat of the universe gets too intense and breaks this magnetic pull too much, the process would fail. But based on their simulations, the pull is likely strong enough to keep the coalescence probability near 100%, meaning the heavy quarks almost always find a dance partner before the music stops.

Technical Summary: Heavy Quark Coalescence Probability in the Presence of a Potential

Problem Statement
In relativistic heavy-ion collisions, heavy quarks (charm and bottom) produced in initial hard scatterings hadronize as the Quark-Gluon Plasma (QGP) cools. While fragmentation dominates at high transverse momentum ( $p_T$ ), coalescence (recombination) with thermal light (anti)quarks is the primary mechanism at low $p_T$ . In the limit $p_T \to 0$ , a static heavy quark cannot hadronize via fragmentation, implying the coalescence probability should approach unity.

However, previous phenomenological studies utilizing a quasi-particle model framework found that the total heavy-quark coalescence probability at the critical temperature ( $T_c$ ) remained significantly below unity. This discrepancy arises because the density of light (anti)quarks, constrained by the lattice-QCD equation of state, is insufficient to yield a probability of unity, even when treating light quarks as massless. Consequently, existing models often rely on ad hoc normalization factors to force the probability to unity at $p_T=0$ , lacking a rigorous dynamical justification.

Methodology
To address this insufficiency, the authors develop a phenomenological heavy-quark potential based on the constituent quark model, motivated by lattice-QCD calculations of the heavy-quark free energy. The methodology proceeds as follows:

Potential Construction: The interaction between a heavy quark and a light antiquark is modeled using a potential $V(r, T)$ comprising a Coulombic term, a string term, and a spin-spin interaction. The potential includes a temperature-dependent screening mass ( $\mu$ ) to account for medium effects. The parameters (strong coupling $\alpha_s$ , heavy quark mass $m_Q$ ) are tuned to reproduce the vacuum masses of pseudoscalar ( $D, B$ ) and vector ( $D^*, B^*$ ) mesons.
Wave Function and Wigner Function: The potential is employed in the Schrödinger equation to determine the heavy meson wave function $\psi(r, T)$ . This wave function defines the Wigner function $\Phi(r, k)$ , which characterizes the spatial and momentum distribution of the bound state.
Coalescence Calculation: The coalescence probability for a static heavy quark is calculated by integrating the product of the light antiquark distribution function (modified by the potential) and the Wigner function over phase space. The light antiquark distribution follows a Fermi-Dirac statistics modified by the potential energy $V^*(r)$ .
Total Probability Estimation: The total coalescence probability is estimated by calculating the probability for the ground-state pseudoscalar meson and scaling it by the ratio of total heavy-hadron density to the ground-state density, as predicted by the statistical hadronization model.
Medium Effects: The study investigates the impact of color screening in the QGP by varying the screening mass parameter $\mu$ . The results are constrained to maintain consistency with the lattice-QCD energy density at $T_c$ ( $\approx 0.4 \text{ GeV/fm}^3$ ).

Key Contributions and Results

Enhancement via Potential: The study demonstrates that incorporating the attractive heavy-quark potential significantly enhances the coalescence probability. The potential concentrates light antiquarks in the vicinity of the heavy quark, increasing the local density available for recombination.
Resolution of the Unity Problem: Unlike previous models that yielded probabilities well below unity, the inclusion of the potential allows the total coalescence probability to approach unity across a broad range of parameters (specifically for light quark masses consistent with lattice-QCD quasiparticle descriptions, $m_q \approx 0.45 \text{ GeV}$ ). This removes the need for ad hoc normalization factors.
Dependence on Meson Radius: The results show a reversal in the uncertainty band ordering compared to models without a potential. Smaller meson radii (0.5 fm) yield higher coalescence probabilities because the attractive potential more effectively concentrates light quarks within the smaller spatial volume.
Sensitivity to Screening: The coalescence probability is highly sensitive to the screening mass $\mu$ . As $\mu$ increases (representing stronger color screening in the medium), the potential well becomes shallower, reducing the spatial concentration of light quarks and lowering the coalescence probability.
Constraints on In-Medium Mass Shifts: The authors find that for the coalescence probability to remain close to unity, the screening parameter $\mu/\sqrt{\sigma}$ must not exceed approximately 0.5. At this limit, the masses of $D$ and $B$ mesons are reduced by roughly 60 MeV relative to their vacuum values. Larger screening masses lead to a substantial suppression of the probability, suggesting that in-medium mass shifts for open heavy-flavor mesons near $T_c$ are unlikely to exceed this magnitude.

Significance
The paper claims that the inclusion of the heavy-quark interaction potential is essential for a consistent description of heavy-quark hadronization via coalescence. By naturally reproducing the vacuum meson spectrum and accounting for the spatial correlation between heavy and light quarks, the model resolves the long-standing issue of insufficient coalescence probability without arbitrary adjustments. Furthermore, the study provides a quantitative constraint on the medium modification of the heavy-quark potential, linking the survival of the coalescence mechanism to specific limits on color screening and in-medium mass shifts near the phase transition.

More like this