Pith. sign in

REVIEW 1 major objections 6 minor 42 references

Heavy quark coalescence probability in the presence of a potential

T0 review · 1 major / 6 minor · reviewed 2026-07-10 · glm-5.2

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Attractive potential solves heavy-quark coalescence puzzle

desk verdict Solid idea, but the meson radius should be computed from the potential they already solved — not treated as a free parameter. read the letter →

arxiv 2607.08707 v1 pith:4L3BNNHB submitted 2026-07-09 hep-ph nucl-th

Taesoo Song , Jiaxing Zhao This is my paper · ORCID
classification hep-phnucl-th PACS 25.75.-q12.38.Mh12.39.Jh
keywords heavyquarkcoalescencehadronizationquark-gluonplasmaheavy-lightpotentialcolorscreeningmesonmasscharm
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Heavy quarks produced in heavy-ion collisions must hadronize through coalescence with light quarks as the plasma cools. Standard coalescence models, treating light quarks as uniformly distributed, predict a total coalescence probability far below the expected value of unity for a stationary heavy quark. This paper introduces a phenomenological heavy-light quark potential, fitted to reproduce the vacuum masses of D, D*, B, and B* mesons, into the coalescence formalism. The attractive potential concentrates light antiquarks around the heavy quark in coordinate space, naturally boosting the coalescence probability to near unity without any ad hoc normalization factor. The authors further show that this conclusion survives moderate in-medium color screening of the potential, but breaks down if screening is too strong, setting an upper bound on the in-medium mass shift of heavy mesons near the phase transition temperature.

What carries the argument

The mechanism is a coordinate-space density enhancement: an attractive potential well, parameterized via a Coulomb-plus-string form with a spin-spin interaction, draws light antiquarks toward the heavy quark. The coalescence probability is computed using a Gaussian Wigner function for the meson wavefunction, integrated over the light-quark phase-space distribution modified by the potential. The potential itself is calibrated to reproduce vacuum meson masses, and medium effects enter through a temperature-dependent screening mass that controls the depth and range of the interaction.

What would settle it

If lattice QCD or experimental determinations of the D- or B-meson radius near the phase transition yield values significantly larger than 0.5 fm, or if the in-medium screening mass exceeds the threshold identified here, the coalescence probability would fall substantially below unity, undermining the claim that the potential alone resolves the deficit.

Watch

Extended reading notes

Core claim

The central finding is that including a realistic attractive heavy-light quark potential in the coalescence model resolves the long-standing deficit in the total heavy-quark coalescence probability at low momentum. The potential enhances the local density of light antiquarks near the heavy quark, raising the probability from well below unity to approximately unity for both charm and bottom quarks. This result holds as long as the in-medium screening mass parameter does not exceed roughly half the square root of the string tension, which translates to a maximum mass reduction of about 60 MeV for D and B mesons near the critical temperature.

Load-bearing premise

The claim that the coalescence probability approaches unity depends on the heavy-meson radius being at the smaller end of the assumed range (0.5 fm). If the actual meson radius is larger, the probability drops below unity even with the potential included.

Editorial extensions

If this is right

  • The coalescence probability constraint bounds the in-medium modification of heavy-light potentials near the QCD phase transition, complementing lattice QCD and QCD sum-rule estimates.
  • Heavy-flavor observables in heavy-ion collisions, such as the baryon-to-meson ratio and elliptic flow at intermediate transverse momentum, can be described without ad hoc normalization of the coalescence probability.
  • The framework can be extended to study coalescence of heavy quarks with strange quarks by adjusting the light-quark mass, with predictions for Ds and Bs production.
  • The sensitivity of the coalescence probability to the screening parameter provides a phenomenological handle on extracting the in-medium heavy-quark potential from experimental heavy-flavor data.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the true meson radius is closer to the upper end of the considered range (1.0 fm), the coalescence probability falls below unity even with the potential, suggesting that tighter experimental or lattice constraints on meson radii would sharpen the conclusion.
  • The constraint on the screening parameter could be cross-checked against independent determinations from quarkonium dissociation patterns or heavy-quark diffusion coefficients in the QGP, potentially tightening or challenging the bound.
  • The reversal of the radius-dependence ordering when the potential is included (smaller radius gives higher probability) implies that the potential effect dominates over the geometric overlap effect, which could be tested by comparing coalescence yields for mesons of different sizes.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

1 major / 6 minor

Summary. This manuscript investigates the role of the heavy-light quark potential in heavy-quark coalescence at the QGP phase boundary. The authors construct a phenomenological potential (Coulomb + string + spin-spin) fitted to reproduce the vacuum masses of D, D*, B, and B* mesons. When this potential is included in the thermal distribution of light antiquarks (Eq. 13), the attractive interaction enhances the local density of light antiquarks around a static heavy quark, increasing the coalescence probability (Eq. 15) toward unity without an ad hoc normalization factor. The authors also study in-medium modifications via a temperature-dependent screening mass, finding that the coalescence probability remains near unity as long as the screening parameter satisfies mu/sqrt(sigma) < ~0.5, which corresponds to a ~60 MeV in-medium mass shift of D and B mesons.

Significance. The problem addressed is well-motivated: within the quasi-particle model framework constrained by the lattice-QCD equation of state, the total heavy-quark coalescence probability at T_c falls significantly below unity, which is physically problematic since a static heavy quark must hadronize via coalescence. The paper provides a concrete mechanism—potential-induced spatial clustering of light antiquarks—to resolve this without ad hoc rescaling. The potential is independently anchored to vacuum spectroscopy (Fig. 1), and the in-medium constraint (mu/sqrt(sigma) < 0.5 implying mass shifts < ~60 MeV) is a falsifiable prediction that can be compared against QCD sum-rule and effective-Lagrangian calculations. The finding that the coalescence probability is more sensitive to the screening mass than the meson mass itself is a useful phenomenological observation.

major comments (1)
  1. [Sec. III, Eqs. (9), (10), and (6)] The heavy-meson radius entering the Wigner function width sigma (Eq. 10) is treated as a free parameter in the range 0.5-1.0 fm, but the same potential used in Eq. (13) is also solved via the Schrodinger equation (Eq. 6) to obtain the meson wave function psi(r, T). This wave function has a definite spatial extent, <r^2> = integral of r^2 |psi(r)|^2 d^3r, which should fix sigma self-consistently. The central claim that the coalescence probability 'approaches unity' holds primarily at the lower end of the radius range (0.5 fm, which corresponds to the upper boundary of the bands in Figs. 3-4 after the potential is included). If the self-consistent radius from the fitted potential is in the 0.7-1.0 fm range, the probability falls below unity even with the potential, which would undermine the conclusion that no ad hoc normalization is needed. The authors should either compute <r^2> from the波
minor comments (6)
  1. [Fig. 2 caption] The caption states 'blue' and 'red' for the without- and with-potential curves, but the figure appears to use a different color scheme in the rendered version. Please verify consistency.
  2. [Sec. III, Eq. (15)] The notation D_3 is used for the spin degeneracy of the D meson, but D_3 was also used in Eq. (8) for the color-spin degeneracy factor of particle 3 (the produced hadron). Clarify whether these are the same quantity or distinct.
  3. [Sec. III, below Eq. (18)] The sentence 'For pseudoscalar D meson, the spin degeneracy factor (D_3 in Eq. (15) equals 1' has an unmatched parenthesis. Please fix.
  4. [Sec. IV, Fig. 5] The figure shows meson masses vs. mu/sqrt(sigma), but the range plotted (0.35-0.70) extends beyond the range discussed in the text (0.35-0.5). It would help to mark the mu/sqrt(sigma) = 0.5 value explicitly on the plot.
  5. [Sec. V] The sentence 'In this work, the heavy-light quark potential that reproduces the masses of both pseudoscalar and vector heavy mesons.' is grammatically incomplete; a verb is missing.
  6. [Sec. IV, last paragraph] The phrase 'At the same time, resulting in a shallower potential...' is a sentence fragment. Please integrate with the preceding sentence.

Circularity Check

0 steps flagged · score 2.0 of 10

No strict circularity; internal inconsistency in radius treatment noted but coalescence probability is genuinely derived

full rationale

The paper's derivation chain is not circular in the strict sense. The heavy-light potential (Eq. 3) is fitted to external data (vacuum D, D*, B, B* meson masses), and the coalescence probability (Eq. 15) is a genuinely derived quantity that depends on the interplay of the potential, the thermal Fermi-Dirac distribution (Eq. 13), and the Wigner function (Eq. 9). The result is not forced to equal the input by construction. The self-citation [24] (Song & Zhao) is used for motivation ('Recently we demonstrated...'), not as a load-bearing mathematical step in the current derivation. The energy-density constraint comes from external lattice QCD calculations [37, 38]. The statistical hadronization ratio (Eq. 19, from [23]) is an external model input, not a self-derived quantity. However, there is a notable internal inconsistency: the Schrödinger equation (Eq. 6) with the fitted potential yields a wave function ψ(r,T) that determines ⟨r²⟩ and hence the Wigner function width σ via Eq. (10), but the paper treats the meson radius as a free parameter (0.5–1.0 fm) rather than computing it self-consistently. The central claim that coalescence probability 'approaches unity' holds primarily at the lower end of this range (0.5 fm, upper band in Figs. 3–4). If the self-consistent radius from the fitted potential were larger (0.7–1.0 fm), the probability would fall below unity even with the potential included, potentially undermining the conclusion that no ad hoc normalization is needed. This is a parameter-selection concern and an internal inconsistency, but it does not constitute circularity as defined: the prediction (coalescence probability) does not reduce to the fitted input (meson masses) by construction, and the paper transparently shows the radius sensitivity rather than hiding it. Score 2 reflects the minor concern that the paper avoids computing a quantity (the radius) that its own framework predicts, which could have constrained or invalidated its main conclusion.

Assumptions & free parameters 4 free parameters · 3 assumptions · 0 invented entities

The model relies on standard constituent quark model parameters and lattice-QCD motivated potential forms. No new particles or forces are invented. The main model dependence is on the meson radius and the statistical hadronization ratios.

free parameters (4)
  • alpha_s = ~1.0 (0.97 for charm, 1.03 for bottom)
    Strong coupling constant fitted to reproduce D/D* and B/B* mass splitting.
  • m_Q = 1.14 GeV (charm), 4.59 GeV (bottom)
    Constituent heavy quark mass fitted to reproduce spin-averaged meson masses.
  • sigma (meson radius) = 0.5-1.0 fm
    Gaussian Wigner function width parameter, treated as a range. The central result depends on this being at the lower end.
  • alpha, beta (in r_0) = 2.2 fm^-1, 0.277
    Parameters in the spin-spin interaction form factor, taken from Ref. [29].
assumptions (3)
  • domain assumption Coalescence probability should approach unity at low p_T.
    Stated in Introduction: 'a static heavy quark cannot hadronize through fragmentation.' This is the physical motivation for the entire study.
  • domain assumption Statistical hadronization model accurately describes relative abundances of charm hadrons.
    Used in Eq. (19) to extrapolate total coalescence probability from D-meson probability alone.
  • domain assumption Gaussian wavefunction is a valid approximation for 1-S state mesons.
    Used in Eq. (9) for the Wigner function. The result's sensitivity to the radius parameter stems from this choice.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Heavy quark coalescence probability in the presence of a potential." pith.science (2026). https://pith.science/paper/4L3BNNHB

@misc{pith2026260708707,
  author       = {Pith},
  title        = {Pith review of: Heavy quark coalescence probability in the presence of a potential},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/4L3BNNHB}},
  note         = {Machine review of arXiv:2607.08707}
}
read the original abstract

In this study, we explore the role of the heavy quark potential in heavy quark coalescence, whose probability is expected to be unity at low momentum. To this end, we develop a phenomenological heavy quark potential based on the constituent quark model that reproduces the vacuum masses of pseudoscalar and vector heavy mesons. Using this potential, we demonstrate its enhancement effect on the coalescence probability. We also investigate how medium-induced modifications of the heavy quark potential in the quark gluon plasma affect the coalescence process. Our results indicate that the coalescence probability remains close to unity as long as the modification of the potential is sufficiently moderate.

Figures

Figures reproduced from arXiv: 2607.08707 by the authors.

Figure 1
Figure 1. FIG. 1: Masses of [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: Charm-quark coalescence probabilities at [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3: (Upper) the energy density correction and (lower) [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figures from the paper (3 more)
Figure 5
Figure 5. Figure 5: shows the masses of the D, D∗ , B and B∗ mesons as functions of the screening mass µ. As ex￾pected, the meson masses decrease with increasing µ. For example, when µ/√ σ increases from 0.35 to 0.5, the D￾meson mass decreases from 1.87 GeV to 1.81 GeV, while the B-meson …
Figure 4
Figure 4. Figure 4: FIG. 4: Bottom-quark coalescence probabilities at [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]
Figure 6
Figure 6. Figure 6: FIG. 6: Charm-quark (upper panel) and bottom-quark (lower [PITH_FULL_IMAGE:figures/full_fig_p006_6.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

42 extracted references · 42 canonical work pages

  1. [1]

    as functions of the strong couplingα s form c = 1.14 GeV andm b = 4.59 GeV. III. HEA VY QUARK COALESCENCE PROBABILITY Heavy quarks hadronize in heavy-ion collisions through either fragmentation or coalescence [21]. Fragmentation dominates at highp T , whereas coalescence becomes the primary hadronization mechanism at lowp T . Conse- quently, a static heav...

    work page

  2. [2]

    Heavy Quarks in the Quark-Gluon Plasma

    R. Rapp and H. van Hees (2010), pp. 111–206, 0903.1096

    work page Pith review arXiv 2010

  3. [3]

    Heavy-flavour and quarkonium production in the LHC era: from proton-proton to heavy-ion collisions

    A. Andronic et al., Eur. Phys. J. C76, 107 (2016), 1506.03981

    work page Pith review arXiv 2016

  4. [4]

    Dong and V

    X. Dong and V. Greco, Prog. Part. Nucl. Phys.104, 97 7 (2019)

    work page 2019

  5. [5]

    Open Heavy-Flavor Production in Heavy-Ion Collisions

    X. Dong, Y.-J. Lee, and R. Rapp, Ann. Rev. Nucl. Part. Sci.69, 417 (2019), 1903.07709

    work page Pith review arXiv 2019

  6. [6]

    J. Zhao, K. Zhou, S. Chen, and P. Zhuang, Prog. Part. Nucl. Phys.114, 103801 (2020), 2005.08277

    work page Pith review arXiv 2020

  7. [7]

    T. Song, I. Grishmanovskii, O. Soloveva, and E. Bratkovskaya, Phys. Rev. C110, 034906 (2024), 2404.00425

    work page Pith review arXiv 2024

  8. [8]

    Transport coefficients of heavy quarks by elastic and radiative scatterings in the strongly interacting quark-gluon plasma

    I. Grishmanovskii, T. Song, C. Greiner, and E. Bratkovskaya, Phys. Rev. D112, 014042 (2025), 2503.22311

    work page Pith review arXiv 2025

Show all 42 references
  1. [9]

    R. J. Fries, B. Muller, C. Nonaka, and S. A. Bass, Phys. Rev. C68, 044902 (2003), nucl-th/0306027

    2003 arXiv

  2. [10]

    Greco, C

    V. Greco, C. M. Ko, and P. Levai, Phys. Rev. C68, 034904 (2003), nucl-th/0305024

    2003 arXiv

  3. [11]

    Greco, C

    V. Greco, C. M. Ko, and P. Levai, Phys. Rev. Lett.90, 202302 (2003), nucl-th/0301093

    2003 arXiv

  4. [12]

    Plumari, V

    S. Plumari, V. Minissale, S. K. Das, G. Coci, and V. Greco, Eur. Phys. J. C78, 348 (2018), 1712.00730

    2018 arXiv

  5. [13]

    Cao, K.-J

    S. Cao, K.-J. Sun, S.-Q. Li, S. Y. F. Liu, W.-J. Xing, G.- Y. Qin, and C. M. Ko, Phys. Lett. B807, 135561 (2020), 1911.00456

    2020 arXiv

  6. [14]

    Greco, C

    V. Greco, C. M. Ko, and R. Rapp, Phys. Lett. B595, 202 (2004), nucl-th/0312100

    2004 arXiv

  7. [15]

    He and R

    M. He and R. Rapp, Phys. Rev. Lett.124, 042301 (2020), 1905.09216

    2020 arXiv

  8. [16]

    T. Song, H. Berrehrah, D. Cabrera, J. M. Torres-Rincon, L. Tolos, W. Cassing, and E. Bratkovskaya, Phys. Rev. C92, 014910 (2015), 1503.03039

    2015 arXiv

  9. [17]

    T. Song, H. Berrehrah, D. Cabrera, W. Cassing, and E. Bratkovskaya, Phys. Rev. C93, 034906 (2016), 1512.00891

    2016 arXiv

  10. [18]

    J. Zhao, J. Aichelin, P. B. Gossiaux, V. Ozvenchuk, and K. Werner, Phys. Rev. C110, 024909 (2024), 2401.17096

    2024 arXiv

  11. [19]

    Adams et al

    J. Adams et al. (STAR), Phys. Rev. Lett.92, 052302 (2004), nucl-ex/0306007

    2004 arXiv

  12. [20]

    S. S. Adler et al. (PHENIX), Phys. Rev. Lett.91, 182301 (2003), nucl-ex/0305013

    2003 arXiv

  13. [21]

    Beraudo et al., Nucl

    A. Beraudo et al., Nucl. Phys. A979, 21 (2018), 1803.03824

    2018 arXiv

  14. [22]

    Zhao et al., Phys

    J. Zhao et al., Phys. Rev. C109, 054912 (2024), 2311.10621

    2024 arXiv

  15. [23]

    Altmann, A

    J. Altmann, A. Dubla, V. Greco, A. Rossi, and P. Skands, Eur. Phys. J. C85, 16 (2025), 2405.19137

    2025 arXiv

  16. [24]

    Song and G

    T. Song and G. Coci, Nucl. Phys. A1028, 122539 (2022), 2104.10987

    2022 arXiv

  17. [25]

    Song and J

    T. Song and J. Zhao, Phys. Rev. C113, 024916 (2026), 2512.06339

    2026

  18. [26]

    H. Satz, J. Phys. G32, R25 (2006), hep-ph/0512217

    2006 arXiv

  19. [27]

    Digal, O

    S. Digal, O. Kaczmarek, F. Karsch, and H. Satz, Eur. Phys. J. C43, 71 (2005), hep-ph/0505193

    2005 arXiv

  20. [28]

    Kaczmarek, Eur

    O. Kaczmarek, Eur. Phys. J. C61, 811 (2009)

    2009

  21. [29]

    R. K. Bhaduri, L. E. Cohler, and Y. Nogami, Nuovo Cim. A65, 376 (1981)

    1981

  22. [30]

    A. Park, P. Gubler, M. Harada, S. H. Lee, C. Nonaka, and W. Park, Phys. Rev. D93, 054035 (2016), 1601.01250

    2016 arXiv

  23. [31]

    Gubler, T

    P. Gubler, T. Song, and S. H. Lee, Phys. Rev. D101, 114029 (2020), 2003.09073

    2020 arXiv

  24. [32]

    Song and H

    T. Song and H. Berrehrah, Phys. Rev. C94, 034901 (2016), 1601.04449

    2016 arXiv

  25. [33]

    Song and Q

    T. Song and Q. Zhou, Phys. Scripta99, 125304 (2024), 2405.08535

    2024 arXiv

  26. [34]

    Cho, Phys

    S. Cho, Phys. Rev. C91, 054914 (2015), 1408.4756

    2015 arXiv

  27. [35]

    Kordell, II, R

    M. Kordell, II, R. J. Fries, and C. M. Ko, Annals Phys. 443, 168960 (2022), 2112.12269

    2022 arXiv

  28. [36]

    J. Zhao, P. B. Gossiaux, T. Song, E. Bratkovskaya, and J. Aichelin, EPJ Web Conf.296, 09014 (2024), 2312.11349

    2024 arXiv

  29. [37]

    Andronic, P

    A. Andronic, P. Braun-Munzinger, M. K. K¨ ohler, A. Mazeliauskas, K. Redlich, J. Stachel, and V. Vislavi- cius, JHEP07, 035 (2021), 2104.12754

    2021 arXiv

  30. [38]

    Borsanyi, G

    S. Borsanyi, G. Endrodi, Z. Fodor, A. Jakovac, S. D. Katz, S. Krieg, C. Ratti, and K. K. Szabo, JHEP11, 077 (2010), 1007.2580

    2010 arXiv

  31. [39]

    Borsanyi, G

    S. Borsanyi, G. Endrodi, Z. Fodor, S. D. Katz, S. Krieg, C. Ratti, and K. K. Szabo, JHEP08, 053 (2012), 1204.6710

    2012 arXiv

  32. [40]

    Plumari, W

    S. Plumari, W. M. Alberico, V. Greco, and C. Ratti, Phys. Rev. D84, 094004 (2011), 1103.5611

    2011 arXiv

  33. [41]

    Moreau, O

    P. Moreau, O. Soloveva, L. Oliva, T. Song, W. Cassing, and E. Bratkovskaya, Phys. Rev. C100, 014911 (2019), 1903.10257

    2019 arXiv

  34. [42]

    Monta˜ na, `A

    G. Monta˜ na, `A. Ramos, L. Tolos, and J. M. Torres- Rincon, Phys. Lett. B806, 135464 (2020), 2001.11877

    2020 arXiv

Pith tools

Reviewed July 10, 2026 · model on record in the stance chip above.