Pith. sign in

REVIEW 3 major objections 4 minor 16 references

Dissociation and regeneration of charmonia within microscopic Langevin simulations

T0 review · 3 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash

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

Pith's one-line read This paper claims that a classical Langevin model, in which charm and anticharm quarks move as Brownian particles and bind through a screened Coulomb potential, reproduces the measured elliptic flow of regenerated charmonia in heavy-ion…

desk verdict The charmonium v2 is the only new result, but a bound-state criterion that counts arbitrarily separated pairs as charmonia undermines the central claim. read the letter →

arxiv 2411.17417 v1 pith:XWHBFNTV submitted 2024-11-26 hep-ph nucl-th

Naomi Oei , Juan M. Torres-Rincon , Hendrik van Hees , Carsten Greiner This is my paper · ORCID
classification hep-phnucl-th
keywords quark-gluonplasmacharmoniumellipticflowLangevinsimulationheavyquarkdiffusiondissociationandregenerationstatisticalhadronizationmodelrelativisticheavy-ioncollisions
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

This paper sets out to show that the full life cycle of charmonia in the quark-gluon plasma---binding, dissociation, and recombination---can be generated microscopically from a classical Langevin description, with no primordial charmonium input. Charm and anticharm quarks are Brownian particles interacting through a screened Coulomb potential, and a pair counts as a bound charmonium when its instantaneous relative energy is negative. In fixed-temperature box simulations the model reaches the expected thermal distribution of relative energies and produces charmonium yields that match the statistical hadronization model once the strong coupling is set to $\alpha_s=0.7$. In an expanding fireball the model predicts that all charmonia are regenerated, and the resulting elliptic flow $v_2$ of inclusive $J/\psi$ rises with transverse momentum and reaches about $11.6\%$ for 20-40% centrality at $\sqrt{s_{NN}}=2.76$ TeV, matching the measured LHC values within uncertainties. If this is right, quarkonium flow in heavy-ion collisions can be understood as the flow of regenerated heavy-quark pairs, and the same classical machinery can be extended to bottomonium and to the nuclear modification factor.

What carries the argument

The load-bearing object is the complex heavy-quark potential $V(r)=-\frac{g^2}{4\pi m_D}-\frac{g^2}{4\pi}\frac{e^{-m_D r}}{r}-i\frac{g^2 T}{4\pi}\phi(m_D r)$: its real part supplies the screened Coulomb attraction that can bind a charm-anticharm pair, and its imaginary part feeds the drag coefficient $\gamma$ that governs Brownian motion in the Langevin equation $d\mathbf{p}=-\gamma\mathbf{p}\,dt+\mathbf{F}\,dt+\sqrt{2\gamma E T\,dt}\,\boldsymbol{\rho}$. The pair is considered a charmonium whenever $E_{c\bar c}=E_c+E_{\bar c}+\operatorname{Re}V(r)-E_{\mathrm{tot}}<0$, a criterion checked at every time step. Box simulations at fixed temperature and volume verify that the relative-energy distribution converges to the classical density of states, while an expanding fireball (boost-invariant longitudinal flow with elliptical transverse expansion) turns the same dynamics into a prediction for $v_2$.

What would settle it

Measure the nuclear modification factor $R_{AA}$ of $J/\psi$ at LHC energies: because the model contains no primordial charmonia, its $R_{AA}$ is determined solely by recombination and would show a characteristic centrality and $p_T$ pattern; observation of strong high-$p_T$ suppression or a different centrality trend than the recombination-only prediction would contradict the claim. A cheaper check is to replace the step-function bound-state criterion with a smooth binding probability and see whether the reproduced $v_2$ leaves the quoted uncertainties.

Watch

Extended reading notes

Core claim

The central discovery is that a single complex potential, with a Debye-screened Coulomb real part and an imaginary part encoding collisional dissociation, is sufficient to drive charmonium dissociation and regeneration in relativistic Langevin simulations. Bound states are identified by the sign of $E_{c\bar c}=E_c+E_{\bar c}+\operatorname{Re}V(r)-E_{\mathrm{tot}}$, and this classical criterion produces the equilibrium density of states in a box and charmonium multiplicities consistent with the statistical hadronization model. When the same dynamics are embedded in a boost-invariant elliptic fireball, the elliptic flow of charm quarks and charmonia is obtained from the transverse momentum anisotropy $v_2=\langle(p_x^2-p_y^2)/(p_x^2+p_y^2)\rangle$. Since the simulation contains no primordial charmonia, the $J/\psi$ $v_2$ is generated entirely by recombination of heavy quarks that have developed flow in the plasma, and for semi-central Pb+Pb collisions at $\sqrt{s_{NN}}=2.76$ TeV it reaches $11.6\%$, consistent with the measured inclusive $J/\psi$ elliptic flow at forward rapidity.

Load-bearing premise

The entire dissociation/regeneration balance rests on the classical step-function criterion that a charm-anticharm pair is bound exactly when $E_c+E_{\bar c}+\operatorname{Re}V(r)-E_{\mathrm{tot}}<0$, with only the real part of the potential acting as the force; if that criterion does not faithfully represent the quantum bound-state spectrum, the rates and the resulting $v_2$ all inherit the error.

Editorial extensions

If this is right

  • If the central claim is correct, the $J/\psi$ elliptic flow measured at the LHC can be produced entirely by regenerated charmonia, so the observed flow does not require a large primordial component.
  • Charmonium $v_2$ then inherits the heavy-quark $v_2$ built up by Langevin diffusion, making quarkonium flow a probe of heavy-quark transport coefficients in the quark-gluon plasma.
  • The quadratic scaling of charmonium yield with the number of charm pairs, validated against the statistical hadronization model, indicates that recombination is the rate-setting mechanism in this picture rather than an instantaneous chemical freeze-out.
  • The same framework can be applied to bottomonium by changing the quark mass, yielding a predicted $v_2$ for $\Upsilon$ states that can be checked against existing and future LHC data.
  • Adding primordial charmonia and computing the nuclear modification factor $R_{AA}$, listed as future work in the paper, would turn the model into a complete prediction for quarkonium suppression and regeneration.

Reading between the lines

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

  • Inference: a model with zero primordial charmonia that still matches the measured $J/\psi$ $v_2$ implies that the elliptic-flow observable alone cannot separate recombination from primordial survival; distinguishing them requires $p_T$-dependent $R_{AA}$ data, which the paper does not yet provide.
  • Inference: the bound-state criterion is a hard sign flip on the relative energy. Replacing it with a smooth probability weight derived from the imaginary part of the potential would test whether the $v_2$ match is robust or an artifact of the classical step function.
  • Inference: because the drag coefficient and the binding force are both derived from the same complex potential, the data comparison is a consistency check of the potential itself; a different screening-mass or coupling parametrization would change both friction and attraction, so the reported agreement is not two independent tunings.
  • Inference: extending the same Langevin machinery to smaller systems (e.g., p-Pb collisions) would give a testable prediction for how the recombination contribution to quarkonium $v_2$ scales with system size and multiplicity.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 4 minor

Summary. The paper presents a classical Langevin simulation in which charm and anticharm quarks propagate through a quark-gluon plasma and interact via a screened, complex Coulomb-like potential. Bound states are identified by a negative relative-energy criterion, and box simulations at fixed temperature are compared with a classical density-of-states formula. The strong coupling is set to alpha_s = 0.7 to match charmonium yields from the statistical hadronization model. The model is then embedded in a boost-invariant fireball, and the elliptic flow of charm quarks and charmonia is computed at RHIC and LHC energies, with the charmonium v2 compared to ALICE inclusive J/psi data. The central claim is that this microscopic Langevin framework can describe charmonium dissociation, regeneration, and elliptic flow in heavy-ion collisions.

Significance. If the bound-state criterion and validation strategy were sound, the model would offer a simple, falsifiable kinetic description of quarkonium regeneration, and the equilibrium check in Fig. 1 would be a useful internal consistency test. The paper also makes a concrete prediction for the centrality and energy dependence of regenerated-charmonium v2. However, the current validation chain rests on a potentially misdefined bound-state criterion, a tuned strong coupling, and a comparison to inclusive data without statistical uncertainties, so the significance as stated is not yet established.

major comments (3)
  1. [Sec. 3, Eq. (2) and Eq. (5)] The real part of the potential in Eq. (2) contains the constant term -g^2/(4 pi m_D), so ReV(r) tends to a negative constant as r goes to infinity. With the stated criterion E_c_bar_c = E_c + E_bar_c + ReV(r) - E_tot < 0, a pair at arbitrarily large separation can be counted as bound unless E_tot is defined to cancel this constant, and the text does not state such a definition. Correspondingly, in Eq. (5) the integrand tends to sqrt(E_rel + g^2/(4 pi m_D)) at large r, so the "bound" yield receives a volume-proportional contribution from unlocalized pairs. Because alpha_s = 0.7 is adjusted to match the SHM benchmark, the comparison in Fig. 2b can absorb this artifact. Please redefine the binding energy relative to the asymptotic potential, for example using V(r) - V(infinity) or an explicit binding-energy threshold, and show that the equilibrium bound yield is independent of the box volume.
  2. [Sec. 3, paragraph 'To be consistent with the SHM results...'] The agreement with the SHM shown in Fig. 2b is not an independent validation because the strong coupling alpha_s(T_c) is tuned to that benchmark. This makes the yield comparison circular as a test of the model. To support the claim that bound-state formation and regeneration occur in the expected manner, the authors should fix alpha_s from an independent observable, or demonstrate that the SHM comparison is insensitive to the tuning within a physically motivated range.
  3. [Sec. 3, Fig. 4 and accompanying text] The comparison to ALICE inclusive J/psi v2 is not quantitative as presented: no statistical uncertainties are shown for the simulation, only 5 initial heavy-quark pairs are used, and the simulation contains only regenerated charmonia while the data are inclusive J/psi. The sentence claiming that the results "match well within uncertainties" therefore lacks a quantitative basis. The authors should show simulation error bars, overlay the ALICE data points directly, and discuss how the absence of primordial charmonia affects the comparison with inclusive data.
minor comments (4)
  1. [Sec. 2] The phrase '𝒑− →0' appears to be a typo for '𝒑→0' in the description of the static limit.
  2. [Sec. 3, Fig. 3 caption] The text says that v2 is positive in all three panels, but then describes a negative dip at low pT for RHIC; please clarify whether the statement refers to the overall trend or to individual pT bins.
  3. [Sec. 3, fireball parametrization] The boundary accelerations a_a and a_b are introduced as free parameters, but their numerical values for the different centrality classes are not given; please provide them or cite a table in Ref. [9] or [10].
  4. [Sec. 3, Fig. 1] The equilibrium comparison in Fig. 1 would be more informative with error bars on the simulation histogram and a statement of the number of events or statistics used.

Circularity Check

2 steps flagged · score 6.0 of 10

Charmonium-yield validation is calibrated to SHM after fitting α_s, and the bound-state criterion counts arbitrarily distant pairs due to the constant in ReV; the ALICE v2 comparison retains partial independence.

  1. fitted input called prediction [Section 3, paragraph after Eq. (6) and Fig. 2b caption]
    "To validate our results, we compare them to the SHM [6–8] in Fig. 2b... To be consistent with the SHM results, we adjust the strong coupling constant, α_s(T_c), selecting a value of α_s = 0.7, to align our results with the SHM predictions."

    The charmonium yield is the output being validated: the section first says 'To validate our results, we compare them to the SHM', then fixes the model's free strong coupling α_s specifically 'to align our results with the SHM predictions'. The agreement in Fig. 2b is therefore produced by the fit, not tested by it. The yield prediction is not independent of the SHM input; it is an interpolation to that benchmark. The v2 comparison with ALICE is not fitted in this way and retains some independent content, but the SHM-yield validation itself reduces to the calibration.

  2. self definitional [Section 3, bound-state criterion paragraph and Eq. (2)]
    "we adopt the criterion that a charm-anticharm pair is considered to be bound if its relative energy is smaller than zero... V(r) = -g^2/(4π m_D) - g^2/(4π) exp(-m_D r)/r - i g^2 T/(4π) φ(m_D r)"

    In Eq. (2) the real part of the potential contains the constant term -g^2/(4π m_D), so ReV(r) tends to -g^2/(4π m_D) as r→∞. With E_c\bar c defined as E_c+E_\bar c+V(r)-E_tot and V=ReV, any pair at arbitrarily large separation has E_c\bar c ≈ -g^2/(4π m_D) < 0 (unless E_tot is defined to cancel the constant, which the text does not specify) and is counted as a bound charmonium. Eq. (5)'s left-hand 'bound' integral then receives a volume-proportional contribution from unlocalized pairs. The resulting 'charmonium' yield and v2 are thus fixed by the operational definition plus the potential's constant shift, not by genuine bound-state dynamics; the tuned α_s then compensates for this threshold artifact. This is a self-definitional construction of the claimed observable.

full rationale

The paper contains two genuine reductions to its own inputs. First, the box charmonium yield is not a prediction: after saying 'To validate our results, we compare them to the SHM', the authors adjust the model's strong coupling α_s specifically 'to align our results with the SHM predictions'; the agreement in Fig. 2b is then produced by the calibration, so the yield comparison is circular. Second, the definition of a charmonium is self-referential in a harmful way: with V(r)=ReV(r) from Eq. (2), the constant -g^2/(4π m_D) makes ReV(∞)<0, so even infinitely separated c-cbar pairs satisfy E_c\bar c<0 and are labeled bound; Eq. (5)'s 'bound' branch includes these unlocalized configurations, and the SHM-tuned α_s absorbs the resulting threshold ambiguity. The central v2 result is partly independent: the fireball is calibrated to light-hadron data and the ALICE charmonium v2 is not used as a fit target, so the v2 comparison provides genuine external contact. That is why the circularity score is intermediate (6) rather than complete (10). Self-citations to Refs. [4] and [10] are used for details and calibration, but they do not themselves smuggle in the claimed v2 result; the load-bearing circularity is the α_s-to-SHM calibration and the bound-state definition.

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

The model's central claims rely on the classical Brownian treatment and the complex potential from Ref. [5], on a bound-state criterion defined by the real part of that potential, and on a fireball tuned to light-hadron data. The only explicit fit to the benchmark used for validation is alpha_s = 0.7, tuned to SHM yields. No new entities are invented.

free parameters (4)
  • strong coupling alpha_s(Tc) = 0.7
    Adjusted in Sec. 3 to align the simulated charmonium yield with the statistical hadronization model; used in the potential (Eq. 2) and drag coefficient (Eq. 4), so it controls both binding and transport.
  • cut-off Lambda in drag coefficient = 4 GeV
    Introduced in Eq. (4) following Ref. [5] to regulate the drag coefficient; no independent constraint is given in this paper.
  • fireball boundary accelerations a_a, a_b = not reported
    Chosen 'to achieve a good agreement with experimental data of pT-spectra and elliptic flow of light hadrons' (Sec. 3), i.e., tuned to external data; values are not reported in this contribution.
  • initial heavy-quark pair multiplicity = 5
    The fireball simulations are carried out with 5 initial HQ pairs; the yield scales quadratically with pair number, and the v2 results depend on this modest number.
assumptions (6)
  • domain assumption Heavy quarks in the quark-gluon plasma obey the classical relativistic Langevin equation (Eq. 1) with drag gamma and Gaussian noise.
    Sec. 2 formulates the HQ dynamics as Brownian motion; no quantum master equation or wave-function evolution is used.
  • domain assumption The c-cbar interaction is described by the complex screened Coulomb potential of Eq. (2) from Ref. [5], with the real part giving the force and the imaginary part giving dissociation and drag.
    Sec. 2 adopts this potential and derives gamma from its imaginary part; the paper does not derive the potential from QCD.
  • ad hoc to paper Bound states are defined by E_c_bar_c < 0 with V = Re V(r) (Sec. 3).
    This is a classical energy criterion adopted for this simulation; it is not justified from the quarkonium spectrum.
  • domain assumption The QGP background is a boost-invariant fireball with volume V(tau) = pi a(tau) b(tau) (z0 + c tau), whose boundary accelerations are tuned to light-hadron data.
    Sec. 3 describes the fireball parametrization; the tuned accelerations are not reported here.
  • domain assumption The statistical hadronization model (SHM) provides the correct equilibrium charmonium yields used to benchmark the simulation.
    Sec. 3 compares simulated yields to SHM (Eq. 6) and uses this comparison to fix alpha_s.
  • standard math The classical density of states with the Boltzmann factor (Eq. 5) describes the equilibrium distribution of pair energies.
    Used to validate the box simulation in Fig. 1.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Dissociation and regeneration of charmonia within microscopic Langevin simulations." pith.science (2026). https://pith.science/paper/XWHBFNTV

@misc{pith2026241117417,
  author       = {Pith},
  title        = {Pith review of: Dissociation and regeneration of charmonia within microscopic Langevin simulations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/XWHBFNTV}},
  note         = {Machine review of arXiv:2411.17417}
}
read the original abstract

We present a classical model to study the formation of charmonia, as well as dissociation and regeneration processes of heavy-quark bound states in the quark gluon plasma using Langevin simulations. The charm and anticharm quarks are described as Brownian particles in the background medium of light quarks and gluons and interact among them over a Coulomb-like screened potential to form bound states, which can dissociate again due to interactions with the medium. Box simulations at fixed temperature and volume are used to verify that the system reaches the expected thermal distribution in the equilibrium limit and to test bound state properties. The medium evolution is then parametrized by a boost-invariant fireball. In this configuration, the elliptic flow of charm and anticharm quarks as well as of charmonia is studied at RHIC and LHC energies.

Figures

Figures reproduced from arXiv: 2411.17417 by the authors.

Figure 1
Figure 1. Simulation of the distribution of the relative energy of the heavy-quark pair in a cubic box with volume 𝑉 = (8 fm) 3 and temperature 𝑇 = 160 MeV, compared to the analytic expectation, Eq. (5). Furthermore, we study the time evolution of the bound-state formation for several numbers of initial pairs, as shown in the left part of Fig. 2a. (a) (b) [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. (a): Time evolution of the fraction of bound states for different numbers of initial HQ pairs in a box of volume 𝑉 = (10 fm) 3 and temperature 𝑇 = 160 MeV. (b): Final charmonium yields for different numbers of HQ pairs from the simulation compared to the results from the SHM. The HQs are initially placed randomly in the box, which allows for a small chance of forming bound states immediately. Over time, recombinatio… view at source ↗
Figure 3
Figure 3. displays the elliptic flow of individual charm and anticharm quarks as a function of the transverse momentum for different energies and centralities. The simulations are carried out with 5 initial HQ pairs and we consider fireball parametrizations for RHIC energies at 20 − 40% centrality and for central (0 − 20%) as well as semi-central (20 − 40%) collisions at LHC energies. (a) HQ 𝑣2 for semi-central collisions (20… view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Charmonium elliptic flow as a function of 𝑝𝑇 for different collision energies and centralities. 4. Conclusions and Outlook We developed a microscopic model to describe charmonia based on relativistic Langevin simulations, which we tested in box simulations at fixed tem…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

16 extracted references · 8 canonical work pages

  1. [4]

    Oei et al.,Formation, dissociation and regeneration of charmonia within microscopic Langevin simulations, 2410.19619

    N. Oei et al.,Formation, dissociation and regeneration of charmonia within microscopic Langevin simulations, 2410.19619

    arXiv

  2. [1]

    Matsui and H

    T. Matsui and H. Satz,𝐽/𝜓 Suppression by Quark-Gluon Plasma Formation, Phys. Lett. B 178 (1986)

    work page 1986

  3. [2]

    M.C.Abreuetal.(NA50Collaboration), Anomalous𝐽/𝜓 suppressioninPb-Pbinteractions at 158 GeV/c per nucleon,Phys. Lett. B410 (1997) 337

    work page 1997

  4. [3]

    PHENIX Collaboration,𝐽/𝜓 Production vs Centrality, Transverse Momentum, and Rapidity in Au+Au Collisions at√𝑠𝑁𝑁 = 200 GeV, Phys. Rev. Lett.98(2007) 232301 [nucl-ex/0611020]

    arXiv 2007

  5. [5]

    Blaizot et al.,Heavy quark bound states in a quark–gluon plasma: Dissociation and recombination,Nuclear Physics A946 (2016) 49

    J.-P. Blaizot et al.,Heavy quark bound states in a quark–gluon plasma: Dissociation and recombination,Nuclear Physics A946 (2016) 49

    work page 2016

  6. [6]

    Andronic et al.,Statistical hadronization of charm in heavy ion collisions at SPS, RHIC and LHC, Phys

    A. Andronic et al.,Statistical hadronization of charm in heavy ion collisions at SPS, RHIC and LHC, Phys. Lett. B571 (2003) 36 [nucl-th/0303036]

    arXiv 2003

  7. [7]

    Andronic et al.,Decoding the phase structure of QCD via particle production at high energy,Nature 561 (2018) 321 [1710.09425]

    A. Andronic et al.,Decoding the phase structure of QCD via particle production at high energy,Nature 561 (2018) 321 [1710.09425]

    arXiv 2018

  8. [8]

    Andronic et al.,The multiple-charm hierarchy in the statistical hadronization model, JHEP 07 (2021) 035 [2104.12754]

    A. Andronic et al.,The multiple-charm hierarchy in the statistical hadronization model, JHEP 07 (2021) 035 [2104.12754]

    arXiv 2021

Show all 16 references
  1. [9]

    van Hees, M

    H. van Hees, M. He and R. Rapp,Pseudo-critical enhancement of thermal photons in relativistic heavy-ion collisions?,Nuclear Physics A933 (2015) 256

    2015

  2. [10]

    Oei,𝐽/𝜓 formation within microscopic Langevin simulations, Master’s thesis, Goethe Universität Frankfurt am Main, 2022

    N. Oei,𝐽/𝜓 formation within microscopic Langevin simulations, Master’s thesis, Goethe Universität Frankfurt am Main, 2022

    2022

  3. [11]

    Cacciari, M

    M. Cacciari, M. Greco and P. Nason,The𝑝𝑇 spectrum in heavy-flavour hadroproduction., JHEP 05 (1998) 007 [hep-ph/9803400]

    1998 arXiv

  4. [12]

    Cacciari, S

    M. Cacciari, S. Frixione and P. Nason,The p(T) spectrum in heavy flavor photoproduction, JHEP 03 (2001) 006 [hep-ph/0102134]

    2001 arXiv

  5. [13]

    T. Song, J. Aichelin and E. Bratkovskaya,Production of primordial𝐽/𝜓 in relativistic𝑝+𝑝 and heavy-ion collisions, Phys. Rev. C96(2017) 014907 [1705.00046]

    2017 arXiv

  6. [14]

    ALICE Collaboration,𝐷 Meson Elliptic Flow in Noncentral Pb-Pb Collisions at√𝑠𝑁𝑁 = 2.76 TeV, Phys. Rev. Lett.111 (2013) 102301

    2013

  7. [15]

    ALICE Collaboration,𝐽/𝜓 Elliptic Flow in Pb-Pb Collisions at√𝑠𝑁𝑁 = 2.76 TeV, Phys. Rev. Lett.111 (2013) 162301

    2013

  8. [16]

    𝐽/𝜓 Elliptic Flow in Pb-Pb Collisions at√𝑠NN= 2.76 TeV

    ALICE Collaboration, “𝐽/𝜓 Elliptic Flow in Pb-Pb Collisions at√𝑠NN= 2.76 TeV.” HEPData (collection), 2013. 7

    2013

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.