pith. sign in

arxiv: 2606.03552 · v1 · pith:KMO54PZNnew · submitted 2026-06-02 · ✦ hep-ph · nucl-th

D⁰-D_s^+ Elliptic-Flow Splitting under Event-Shape Engineering: A Probe of Sequential Charm Hadronization

Pith reviewed 2026-06-28 09:25 UTC · model grok-4.3

classification ✦ hep-ph nucl-th
keywords sequential charm hadronizationelliptic flow splittingevent-shape engineeringD0 and Ds+ mesonscharm flow responseQCD transition temperaturePb-Pb collisions at 5.02 TeV
0
0 comments X

The pith

Event-shape engineering shows sequential charm hadronization by making the D0-Ds+ elliptic flow splitting grow positive with initial geometry.

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

The paper tests whether open-charm mesons form at different times in the cooling quark-gluon plasma. It applies event-shape engineering to select Pb-Pb collisions with varying initial eccentricity via the q2 parameter. In the sequential picture, Ds+ forms earlier than D0, which produces a positive flow difference that strengthens with larger q2 and a steeper response slope for D0. The simultaneous formation baseline instead yields near-zero or negative splitting without the same geometry scaling. The 30-50 percent centrality window maximizes the signal because of the interplay between QGP lifetime and eccentricity.

Core claim

By selecting events on q2 in 0-10 percent and 30-50 percent centrality classes at 5.02 TeV, the positive Delta v2(D0-Ds+) grows systematically with q2 in the sequential scenario, and the response slope chi obeys chi(D0) greater than chi(Ds+), a hierarchy absent when both species form simultaneously; the q2 ratios of Ds+/D0 yields stay near unity, confirming the splitting is a dynamical flow effect rather than a yield change.

What carries the argument

The response slope chi of elliptic flow to the event-shape parameter q2, which isolates the geometry-driven flow conversion from the timing of hadronization.

If this is right

  • The positive splitting and species-dependent chi hierarchy appear only under sequential formation and scale with q2.
  • Semi-central collisions provide the clearest window due to the non-monotonic lifetime-eccentricity interplay.
  • Yield ratios remaining near unity isolates the effect as flow conversion rather than chemical modification.
  • Delta v2 and chi under ESE serve as differential probes of formation timing near the QCD transition.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same ESE method could be applied to other charm species such as D+ or Lambda_c to map a fuller formation sequence.
  • If confirmed, the hierarchy would constrain the temperature window between 1.2 Tc and Tc for different bound states.
  • This differential probe might resolve tensions in inclusive heavy-flavor flow data by separating geometry from timing.
  • Extension to lower beam energies could test whether the sequential window shrinks or widens with changing QGP lifetime.

Load-bearing premise

The hydrodynamic or transport model used to generate the v2 and q2 distributions correctly encodes the space-time evolution and species-dependent hadronization times without extra biases that could create the reported splitting and slope hierarchy.

What would settle it

Data showing that Delta v2(D0-Ds+) stays near zero or decreases with rising q2, or that the response slopes lack the D0-greater-than-Ds+ ordering in 30-50 percent collisions, would falsify the sequential prediction.

Figures

Figures reproduced from arXiv: 2606.03552 by Ben-Wei Zhang, Enke Wang, Jiaxing Zhao, Tan Luo, Wei Dai, Yu-Jie Huang.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: ). The q2 dependence is examined in the 3.25 < pT < 3.75 GeV/c interval, where the positive sequential ∆v2(D0−D+ s ) is already prominent and rising toward its maximum (see [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6 [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: shows the D+ s /D0 ratio for large- and small-q2 event classes in two centrality classes in the sequential scenario. The ratio rises with decreasing pT at low pT , reaching a pronounced peak in the 0–10% centrality class, and then falls toward high pT . The curves for the 0–10% centrality lie systematically above those for the 30–50% centrality, indicating stronger strange-meson production in the denser ce… view at source ↗
Figure 8
Figure 8. Figure 8: makes this control test more explicit by form￾ing the q2 double ratio of D+ s /D0 , defined as Rq2 ≡ (D+ s /D0 )large-q2 /(D+ s /D0 )small-q2 . If the event-shape se￾lection strongly modified charm hadron chemistry, this observable would move substantially away from unity. Instead, the double ratio remains close to one over the full pT range for both centrality classes. The sequential and simultaneous calc… view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9 [PITH_FULL_IMAGE:figures/full_fig_p009_9.png] view at source ↗
read the original abstract

Recent work has proposed sequential hadronization of open-charm hadrons in the quark-gluon plasma, wherein more tightly bound species such as $D_s^+$ form earlier near $1.2 T_c$ and $D^0$ forms later at $T_c$. That work showed that this mechanism naturally reverses the sign of the $D^0-D_s^+$ elliptic-flow splitting relative to the conventional simultaneous baseline. In this work, we demonstrate that event-shape engineering (ESE) provides a sharper discrimination between the two pictures than inclusive measurements alone. By selecting large-$q_2$ and small-$q_2$ events in 0--10\% and 30--50\% centrality classes in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV, we show that the geometry-driven enhancement of charm-meson $v_2$ can be separated from the hadronization-time response: the positive $\Delta v_2(D^0-D_s^+)$ in the sequential scenario grows systematically with $q_2$, while the corresponding response slope $\chi$ reveals a species-dependent hierarchy $\chi(D^0) > \chi(D_s^+)$ that is robust against the overall flow normalization and absent in the simultaneous baseline. In the simultaneous case, the splitting is near zero or negative and does not follow the same geometry scaling. Notably, the semi-central 30--50\% class emerges as the optimal window, because the non-monotonic interplay between QGP lifetime and initial eccentricity maximizes the late-stage flow conversion. The $q_2$ ratios of the $D_s^+/D^0$ yield ratio remain close to unity, confirming that the splitting is a dynamical flow effect rather than a chemical yield modification. These results establish $\Delta v_2(D^0-D_s^+)$ and the response slope $\chi$ under ESE as complementary differential probes of the space-time structure of charm hadronization near the QCD transition temperature.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript claims that event-shape engineering (ESE) via q2 selection in 0-10% and 30-50% centrality classes of Pb-Pb collisions at 5.02 TeV provides a sharper probe of sequential charm hadronization than inclusive v2 measurements. In the sequential scenario (Ds+ near 1.2 Tc, D0 at Tc), positive Δv2(D0-Ds+) grows systematically with q2 while the response slope χ exhibits the hierarchy χ(D0) > χ(Ds+), robust against overall flow normalization; the simultaneous baseline instead yields near-zero or negative splitting without the same geometry scaling. The 30-50% class is identified as optimal due to QGP lifetime and eccentricity interplay, and Ds+/D0 yield ratios near unity confirm the effect is dynamical rather than chemical.

Significance. If the central discrimination holds, the work supplies a useful differential observable for mapping the space-time structure of charm hadronization near the QCD transition. The separation of geometry-driven flow enhancement from hadronization-time response via ESE, together with the identification of an optimal semi-central window, adds practical guidance for future measurements.

major comments (2)
  1. [Model and simulation description] The discrimination between sequential and simultaneous scenarios rests on a single hydrodynamic/transport realization without reported variations of charm-medium coupling (drag/diffusion coefficients) or alternative space-time mappings of formation hypersurfaces. This leaves open whether the reported χ(D0) > χ(Ds+) hierarchy and positive q2 scaling of Δv2 arise specifically from hadronization timing or from species-dependent transport biases already present in the model.
  2. [Results on response slopes and robustness tests] The assertion that the χ hierarchy is robust against overall flow normalization is stated but not demonstrated by explicit rescaling or parameter variation; without such checks, it is unclear whether the hierarchy survives changes to the flow response parameters that are independent of the hadronization temperatures.
minor comments (1)
  1. [Abstract] The abstract references specific hadronization temperatures but supplies no model equations, parameter values, or statistical uncertainties, making it difficult to assess the quantitative strength of the claimed q2 dependence and χ hierarchy.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below and indicate where revisions will be made to improve the presentation and robustness of the results.

read point-by-point responses
  1. Referee: [Model and simulation description] The discrimination between sequential and simultaneous scenarios rests on a single hydrodynamic/transport realization without reported variations of charm-medium coupling (drag/diffusion coefficients) or alternative space-time mappings of formation hypersurfaces. This leaves open whether the reported χ(D0) > χ(Ds+) hierarchy and positive q2 scaling of Δv2 arise specifically from hadronization timing or from species-dependent transport biases already present in the model.

    Authors: We agree that the analysis uses one hydrodynamic plus transport realization. However, both scenarios employ identical drag and diffusion coefficients; the sole difference is the formation temperature (and thus the hypersurface) for each species. The ESE q2 selection isolates the timing effect because the simultaneous baseline produces no positive Δv2 growth with q2, while the sequential case does. Species-dependent transport biases would be present equally in both scenarios and cannot explain the differential q2 scaling. We will add a paragraph in Sec. II clarifying this point and referencing earlier sensitivity studies on coupling variations, while noting that a full parameter scan lies beyond the present scope. revision: partial

  2. Referee: [Results on response slopes and robustness tests] The assertion that the χ hierarchy is robust against overall flow normalization is stated but not demonstrated by explicit rescaling or parameter variation; without such checks, it is unclear whether the hierarchy survives changes to the flow response parameters that are independent of the hadronization temperatures.

    Authors: We thank the referee for this observation. Although the differential response is expected to be independent of overall normalization, we did not show explicit rescaling. In the revised manuscript we will add a short subsection (or appendix) performing a constant-factor rescaling of the v2 values (e.g., ×0.8 and ×1.2) and demonstrate that the χ(D0) > χ(Ds+) ordering remains unchanged. This will make the robustness explicit. revision: yes

Circularity Check

0 steps flagged

No significant circularity; model predictions remain independent of fitted inputs

full rationale

The provided abstract and context describe hydrodynamic/transport simulations that implement sequential (Ds+ at ~1.2 Tc, D0 at Tc) versus simultaneous hadronization and then compute ESE-selected v2 splitting and response slopes χ. No equations, parameter fits, or self-citations are quoted that reduce the reported Δv2(q2) growth or χ(D0) > χ(Ds+) hierarchy to a tautology or to the same fitted values by construction. The discrimination is presented as an output of the space-time evolution under the two scenarios, with the q2 ratios of yields stated to remain near unity as an internal consistency check rather than a redefinition of the observable. This constitutes a standard model-based prediction chain without the enumerated circular patterns.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on hydrodynamic modeling of QGP evolution and a temperature-dependent sequential hadronization prescription whose specific temperatures are introduced to differentiate the two scenarios.

free parameters (2)
  • Ds+ hadronization temperature = 1.2 Tc
    Set near 1.2 Tc to realize earlier formation in the sequential picture
  • D0 hadronization temperature = Tc
    Set at Tc to realize later formation in the sequential picture
axioms (2)
  • domain assumption Hydrodynamic evolution accurately captures the collective flow response to initial eccentricity
    Invoked to translate q2 selection into geometry-driven v2 enhancement
  • domain assumption Charm quarks thermalize sufficiently to inherit the medium flow before hadronization
    Required for v2 to develop differently according to formation time

pith-pipeline@v0.9.1-grok · 5932 in / 1712 out tokens · 36756 ms · 2026-06-28T09:25:45.171905+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

57 extracted references · 26 linked inside Pith

  1. [1]

    Beraudo et al

    A. Beraudo et al. , Extraction of Heavy-Flavor Trans- port Coefficients in QCD Matter, Nucl. Phys. A 979, 21 (2018) , arXiv:1803.03824 [nucl-th]

  2. [2]

    Dong, Y.-J

    X. Dong, Y.-J. Lee, and R. Rapp, Open Heavy-Flavor Production in Heavy-Ion Collisions, Annu. Rev. Nucl. Part. Sci. 69, 417 (2019) . 10

  3. [3]

    Greco, C

    V. Greco, C. M. Ko, and P. Levai, Parton coalescence at RHIC, Phys. Rev. C 68, 034904 (2003) , arXiv:nucl- th/0305024

  4. [4]

    R. J. Fries, B. Muller, C. Nonaka, and S. A. Bass, Hadron production in heavy ion collisions: Fragmentation and recombination from a dense parton phase, Phys. Rev. C 68, 044902 (2003) , arXiv:nucl-th/0306027

  5. [5]

    R. J. Fries, V. Greco, and P. Sorensen, Coalescence Models For Hadron Formation From Quark Gluon Plasma, Ann. Rev. Nucl. Part. Sci. 58, 177 (2008) , arXiv:0807.4939 [nucl-th]

  6. [6]

    T. Song, H. Berrehrah, D. Cabrera, J. M. Torres-Rincon, L. Tolos, W. Cassing, and E. Bratkovskaya, Tomogra- phy of the Quark-Gluon-Plasma by Charm Quarks, Phys. Rev. C 92, 014910 (2015) , arXiv:1503.03039 [nucl-th]

  7. [7]

    Plumari, V

    S. Plumari, V. Minissale, S. K. Das, G. Coci, and V. Greco, Charmed Hadrons from Coalescence plus Fragmentation in relativistic nucleus-nucleus collisions at RHIC and LHC, Eur. Phys. J. C 78, 348 (2018) , arXiv:1712.00730 [hep-ph]

  8. [8]

    Beraudo, A

    A. Beraudo, A. De Pace, M. Monteno, M. Nardi, and F. Prino, Development of heavy-flavour flow-harmonics in high-energy nuclear collisions, JHEP 02 (02), 043, arXiv:1712.00588 [hep-ph]

  9. [9]

    He and R

    M. He and R. Rapp, Charm-Baryon Production in Proton-Proton Collisions, Phys. Lett. B 795, 117 (2019) , arXiv:1902.08889 [nucl-th]

  10. [10]

    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, Charmed hadron chemistry in rel- ativistic heavy-ion collisions, Phys. Lett. B 807, 135561 (2020), arXiv:1911.00456 [nucl-th]

  11. [11]

    Andronic, P

    A. Andronic, P. Braun-Munzinger, M. K. Köh- ler, A. Mazeliauskas, K. Redlich, J. Stachel, and V. Vislavicius, The multiple-charm hierarchy in the statistical hadronization model, JHEP 07 (07), 035, arXiv:2104.12754 [hep-ph]

  12. [12]

    J. Zhao, J. Aichelin, P. B. Gossiaux, V. Ozvenchuk, and K. Werner, Heavy-flavor hadron production in relativistic heavy ion collisions at energies available at BNL RHIC and at the CERN LHC in the EPOS4HQ framework, Phys. Rev. C 110, 024909 (2024) , arXiv:2401.17096 [hep- ph]

  13. [13]

    Nahrgang, J

    M. Nahrgang, J. Aichelin, P. B. Gossiaux, and K. Werner, Influence of hadronic bound states above Tc on heavy- quark observables in Pb + Pb collisions at at the CERN Large Hadron Collider, Phys. Rev. C 89, 014905 (2014) , arXiv:1305.6544 [hep-ph]

  14. [14]

    Bellwied, S

    R. Bellwied, S. Borsanyi, Z. Fodor, S. D. Katz, and C. Ratti, Flavor dependence of the QCD transition on physical quark masses, Phys. Rev. Lett. 111, 202302 (2013), arXiv:1305.6297 [hep-lat]

  15. [15]

    Cao, G.-Y

    S. Cao, G.-Y. Qin, and S. A. Bass, Energy loss, hadronization and hadronic interactions of heavy fla- vors in relativistic heavy-ion collisions, Phys. Rev. C 92, 024907 (2015) , arXiv:1505.01413 [nucl-th]

  16. [16]

    Beraudo, A

    A. Beraudo, A. De Pace, M. Monteno, M. Nardi, and F. Prino, Heavy flavors in heavy-ion collisions: quench- ing, flow and correlations, Eur. Phys. J. C 75, 121 (2015) , arXiv:1410.6082 [hep-ph]

  17. [17]

    He and R

    M. He and R. Rapp, Hadronization and Charm-Hadron Ratios in Heavy-Ion Collisions, Phys. Rev. Lett. 124, 042301 (2020) , arXiv:1905.09216 [nucl-th]

  18. [18]

    Zhao et al., Hadronization of heavy quarks, Phys

    J. Zhao et al., Hadronization of heavy quarks, Phys. Rev. C 109, 054912 (2024) , arXiv:2311.10621 [hep-ph]

  19. [19]

    Z.-X. Xu, W. Dai, B.-W. Zhang, J. Zhao, and P. Zhuang, Sequential hadronization of open-charm mesons in Pb-Pb collisions, arXiv:2510.16299 [nucl-th] (2025), arXiv:2510.16299

  20. [20]

    Acharya et al

    S. Acharya et al. (ALICE), Event-shape engineering for the D-meson elliptic flow in mid-central Pb-Pb col- lisions at √sNN = 5 .02 TeV, JHEP 02 (02), 150, arXiv:1809.09371 [nucl-ex]

  21. [21]

    Acharya et al

    S. Acharya et al. (ALICE), Transverse-momentum and event-shape dependence of D-meson flow harmonics in Pb–Pb collisions at √sN N = 5.02 TeV, Phys. Lett. B 813, 136054 (2021) , arXiv:2005.11131 [nucl-ex]

  22. [22]

    M. L. Sambataro, Y. Sun, V. Minissale, S. Plumari, and V. Greco, Event-shape engineering analysis of D meson in ultrarelativistic heavy-ion collisions, Eur. Phys. J. C 82, 833 (2022) , arXiv:2206.03160 [hep-ph]

  23. [23]

    L.-G. Pang, H. Petersen, and X.-N. Wang, Pseudo- rapidity distribution and decorrelation of anisotropic flow within the open-computing-language implementa- tion CL Visc hydrodynamics, Phys. Rev. C 97, 064918 (2018), arXiv:1802.04449 [nucl-th]

  24. [24]

    Wu, G.-Y

    X.-Y. Wu, G.-Y. Qin, L.-G. Pang, and X.-N. Wang, (3+1)-D viscous hydrodynamics at finite net baryon den- sity: Identified particle spectra, anisotropic flows, and flow fluctuations across energies relevant to the beam- energy scan at RHIC, Phys. Rev. C 105, 034909 (2022) , arXiv:2107.04949 [hep-ph]

  25. [25]

    Cacciari, S

    M. Cacciari, S. Frixione, N. Houdeau, M. L. Mangano, P. Nason, and G. Ridolfi, Theoretical predictions for charm and bottom production at the LHC, JHEP 10 (10), 137, arXiv:1205.6344 [hep-ph]

  26. [26]

    M. L. Miller, K. Reygers, S. J. Sanders, and P. Stein- berg, Glauber modeling in high energy nuclear collisions, Ann. Rev. Nucl. Part. Sci. 57, 205 (2007) , arXiv:nucl- ex/0701025

  27. [27]

    S. Wang, W. Dai, B.-W. Zhang, and E. Wang, Diffusion of charm quarks in jets in high-energy heavy-ion colli- sions, Eur. Phys. J. C 79, 789 (2019) , arXiv:1906.01499 [nucl-th]

  28. [28]

    W. Dai, S. Wang, S.-L. Zhang, B.-W. Zhang, and E. Wang, Transverse Momentum Balance and Angular Distribution of b¯b Dijets in Pb+Pb collisions, Chin. Phys. C 44, 104105 (2020) , arXiv:1806.06332 [nucl-th]

  29. [29]

    S. Wang, W. Dai, B.-W. Zhang, and E. Wang, Ra- dial profile of bottom quarks in jets in high-energy nuclear collisions, Chin. Phys. C 45, 064105 (2021) , arXiv:2012.13935 [nucl-th]

  30. [30]

    S. Wang, W. Dai, B.-W. Zhang, and E. Wang, Z 0 boson associated b-jet production in high-energy nu- clear collisions*, Chin. Phys. C 47, 054102 (2023) , arXiv:2005.07018 [hep-ph]

  31. [31]

    Cao, G.-Y

    S. Cao, G.-Y. Qin, and S. A. Bass, Heavy-quark dynam- ics and hadronization in ultrarelativistic heavy-ion colli- sions: Collisional versus radiative energy loss, Phys. Rev. C 88, 044907 (2013) , arXiv:1308.0617 [nucl-th]

  32. [32]

    P. B. Gossiaux, R. Bierkandt, and J. Aichelin, Tomogra- phy of a quark gluon plasma at RHIC and LHC energies, Phys. Rev. C 79, 044906 (2009) , arXiv:0901.0946 [hep- ph]

  33. [33]

    Kubo, The fluctuation-dissipation theorem, Rept

    R. Kubo, The fluctuation-dissipation theorem, Rept. Prog. Phys. 29, 255 (1966)

  34. [34]

    Altenkort, O

    L. Altenkort, O. Kaczmarek, R. Larsen, S. Mukherjee, P. Petreczky, H.-T. Shu, and S. Stendebach (HotQCD), Heavy Quark Diffusion from 2+1 Flavor Lattice QCD 11 with 320 MeV Pion Mass, Phys. Rev. Lett. 130, 231902 (2023), arXiv:2302.08501 [hep-lat]

  35. [35]

    Xue, Z.-X

    X.-F. Xue, Z.-X. Xu, W. Dai, J. Zhao, and B.-W. Zhang, Bayesian inference of heavy-quark dissipation and jet transport parameters from D-meson observables in heavy-ion collisions at energies available at the CERN Large Hadron Collider, Phys. Rev. C 113, 054902 (2026) , arXiv:2512.07169 [nucl-th]

  36. [36]

    Guo and X.-N

    X.-f. Guo and X.-N. Wang, Multiple scattering, par- ton energy loss and modified fragmentation functions in deeply inelastic e A scattering, Phys. Rev. Lett. 85, 3591 (2000), arXiv:hep-ph/0005044

  37. [37]

    Zhang, E

    B.-W. Zhang, E. Wang, and X.-N. Wang, Heavy quark energy loss in nuclear medium, Phys. Rev. Lett. 93, 072301 (2004) , arXiv:nucl-th/0309040

  38. [38]

    Zhang, E.-k

    B.-W. Zhang, E.-k. Wang, and X.-N. Wang, Multiple par- ton scattering in nuclei: Heavy quark energy loss and modified fragmentation functions, Nucl. Phys. A 757, 493 (2005) , arXiv:hep-ph/0412060

  39. [39]

    Majumder, Hard collinear gluon radiation and mul- tiple scattering in a medium, Phys

    A. Majumder, Hard collinear gluon radiation and mul- tiple scattering in a medium, Phys. Rev. D 85, 014023 (2012), arXiv:0912.2987 [nucl-th]

  40. [40]

    Deng and X.-N

    W.-t. Deng and X.-N. Wang, Multiple Parton Scatter- ing in Nuclei: Modified DGLAP Evolution for Frag- mentation Functions, Phys. Rev. C 81, 024902 (2010) , arXiv:0910.3403 [hep-ph]

  41. [41]

    Xie, S.-Y

    M. Xie, S.-Y. Wei, G.-Y. Qin, and H.-Z. Zhang, Extract- ing jet transport coefficient via single hadron and di- hadron productions in high-energy heavy-ion collisions, Eur. Phys. J. C 79, 589 (2019) , arXiv:1901.04155 [hep- ph]

  42. [42]

    Ravagli and R

    L. Ravagli and R. Rapp, Quark Coalescence based on a Transport Equation, Phys. Lett. B 655, 126 (2007) , arXiv:0705.0021 [hep-ph]

  43. [43]

    J. Zhao, S. Shi, N. Xu, and P. Zhuang, Sequential Coales- cence with Charm Conservation in High Energy Nuclear Collisions (2018), arXiv:1805.10858 [hep-ph]

  44. [44]

    J. Zhao, S. Shi, N. Xu, and P. Zhuang, Sequential hadronization in heavy ion collisions, Nucl. Phys. A 1005, 121898 (2021) , arXiv:2004.12305 [hep-ph]

  45. [45]

    J. Zhao, J. Aichelin, P. B. Gossiaux, and K. Werner, Cor- relations between heavy mesons and the creation of the charmonia, bottomonia, and Bc mesons in high energy pp collisions (2025), arXiv:2511.08796 [hep-ph]

  46. [46]

    S. Shi, J. Zhao, and P. Zhuang, Heavy flavor dissociation in framework of multi-body Dirac equations, Chin. Phys. C 44, 084101 (2020) , arXiv:1905.10627 [nucl-th]

  47. [47]

    H. w. Crater and P. Van Alstine, TWO-BODY DIRAC EQUATIONS, Annals Phys. 148, 57 (1983)

  48. [48]

    H. W. Crater and P. Van Alstine, Two-body Dirac Equa- tions for Particles Interacting Through World Scalar and Vector Potentials, Phys. Rev. D 36, 3007 (1987)

  49. [49]

    Sazdjian, The Connection of Two Particle Relativistic Quantum Mechanics With the Bethe-Salpeter Equation, J

    H. Sazdjian, The Connection of Two Particle Relativistic Quantum Mechanics With the Bethe-Salpeter Equation, J. Math. Phys. 28, 2618 (1987)

  50. [50]

    Sazdjian, N BODY BOUND STATE RELATIVISTIC W A VE EQUATIONS,Annals Phys

    H. Sazdjian, N BODY BOUND STATE RELATIVISTIC W A VE EQUATIONS,Annals Phys. 191, 52 (1989)

  51. [51]

    S. P. Klevansky, The Nambu–Jona-Lasinio model of quantum chromodynamics, Rev. Mod. Phys. 64, 649 (1992)

  52. [52]

    Bazavov et al

    A. Bazavov et al. (HotQCD), Chiral crossover in QCD at zero and non-zero chemical potentials, Phys. Lett. B 795, 15 (2019) , arXiv:1803.02706 [hep-lat]

  53. [53]

    Peterson, D

    C. Peterson, D. Schlatter, I. Schmitt, and P. M. Zerwas, Scaling Violations in Inclusive e+ e- Annihilation Spec- tra, Phys. Rev. D 27, 105 (1983)

  54. [54]

    S. K. Das, J. M. Torres-Rincon, L. Tolos, V. Minissale, F. Scardina, and V. Greco, Propagation of heavy baryons in heavy-ion collisions, Phys. Rev. D 94, 114039 (2016) , arXiv:1604.05666 [nucl-th]

  55. [55]

    Lisovyi, A

    M. Lisovyi, A. Verbytskyi, and O. Zenaiev, Combined analysis of charm-quark fragmentation-fraction measure- ments, Eur. Phys. J. C 76, 397 (2016) , arXiv:1509.01061 [hep-ex]

  56. [56]

    M. He, R. J. Fries, and R. Rapp, Ds-Meson as Quan- titative Probe of Diffusion and Hadronization in Nu- clear Collisions, Phys. Rev. Lett. 110, 112301 (2013) , arXiv:1204.4442 [nucl-th]

  57. [57]

    J. M. Torres-Rincon, G. Montaña, À. Ramos, and L. To- los, In-medium kinetic theory of D mesons and heavy- flavor transport coefficients, Phys. Rev. C 105, 025203 (2022), arXiv:2106.01156 [hep-ph]