{"id":"63a7ea61-c58c-400e-8258-f7ba7e8dc1f2","arxiv_id":"2411.17417","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Charm quarks treated as Brownian particles in a hot medium form, dissociate, and reform charmonia, and the resulting elliptic flow roughly matches ALICE measurements once the strong coupling is tuned to the statistical hadronization model.","lead":"A simulation treats charm quarks as particles randomly jostled by the hot nuclear medium, lets them bind into charmonia through a screened attraction, and computes how these charmonia flow in heavy-ion collisions. The results are compared with ALICE measurements to see if this simple microscopic picture can describe charmonium regeneration.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Bound-state criterion counts far-separated pairs as charmonia because ReV(r) tends to a negative constant; the claimed yield and elliptic flow may reflect mislabeled pairs.","rationale":"The reader's weakest assumption is that the classical E_rel<0 criterion may not faithfully represent the quantum bound-state spectrum. My concern is more specific and more internal: even within the classical treatment, the energy threshold is inconsistent with the potential's asymptotic value. The constant self-energy term in Eq. (2) shifts the continuum threshold to -g^2/(4 pi m_D), so the stated criterion labels spatially unlocalized pairs as charmonia. This directly affects the SHM validation (where alpha_s is tuned) and the ALICE v2 comparison, because the simulated charmonium population is partly a numerical artifact of the threshold choice. The proposed radial-distribution check would settle whether this is actually happening. The paper's central claim is still plausible in principle, but the current evidence is conditional on correcting the bound-state definition and showing that the v2 result survives that correction. Therefore the verdict should remain CONDITIONAL, with a sharper and more concrete condition than the reader's. I do not see grounds for outright rejection without the test, because the underlying Langevin framework and the independent v2 comparison could still be valid if the threshold issue is resolved.","tokens_in":5638,"tokens_out":9376,"duration_ms":88960,"concrete_test":"Re-analyze the box simulation of Fig. 1 (or rerun it) and plot the radial distribution of pairs satisfying E_rel<0. For a genuine charmonium bound state, the pair separation should cluster at r of order a few times 1/m_D; if a large fraction (e.g., >10%) have r >> 1/m_D, the criterion misclassifies unbound pairs. As a second check, repeat the alpha_s=0.7 tuning and the fireball v2 calculation with the potential shifted so that ReV(infinity)=0 and require E_rel<V(infinity); if the SHM match requires a different alpha_s or the ALICE v2 comparison degrades, the present comparison does not support the claim.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim rests on the operational definition of a bound pair in Sec. 3: a pair is bound if E_rel = E_c + E_anticharm + V(r) - E_tot < 0, with V(r)=ReV(r). The real part of Eq. (2) contains the constant term -g^2/(4 pi m_D), so ReV(r) -> -g^2/(4 pi m_D) as r -> infinity. A pair at rest at arbitrarily large separation therefore has E_rel ~ -g^2/(4 pi m_D) < 0 and is counted as a charmonium state. In Eq. (5), the density of states with E_rel<0 receives a volume-proportional contribution from unlocalized configurations (the integrand tends to sqrt(E_rel + g^2/(4 pi m_D)) at large r), so the 'bound' yield is not a genuine quarkonium yield. Adjusting alpha_s = 0.7 to match the SHM then absorbs the arbitrary threshold, and the v2 from these mislabeled pairs does not test the claimed regeneration mechanism. Unless E_tot is defined to cancel the constant (not stated in the text), this is an internal inconsistency of the present criterion, independent of the quantum-versus-classical question.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":5900,"tokens_out":4687,"duration_ms":48336,"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":[{"comment":"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.","section":"Sec. 3, Eq. (2) and Eq. (5)"},{"comment":"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.","section":"Sec. 3, paragraph 'To be consistent with the SHM results...'"},{"comment":"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.","section":"Sec. 3, Fig. 4 and accompanying text"}],"minor_comments":[{"comment":"The phrase '𝒑− →0' appears to be a typo for '𝒑→0' in the description of the static limit.","section":"Sec. 2"},{"comment":"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.","section":"Sec. 3, Fig. 3 caption"},{"comment":"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].","section":"Sec. 3, fireball parametrization"},{"comment":"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.","section":"Sec. 3, Fig. 1"}],"recommendation":"major_revision","confidential_remarks":"This is a proceedings contribution, and some of the missing statistics and error bars could be considered acceptable for a short proceedings paper if framed as preliminary. However, the bound-state criterion issue in Sec. 3 is not merely a presentation problem; it directly affects the meaning of the charmonium yield and v2. The paper should not be published in its present form until that criterion is clarified or corrected and the tuning of alpha_s is justified as more than a fit to the SHM benchmark."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the only genuinely new result here is the charmonium v2 computed in a fireball, and it is compromised by a broken bound-state criterion. The stress-test note holds up. The real part of the potential in Eq. (2) contains the constant -g^2/(4π m_D), so ReV(r) tends to that constant at large r. A stationary pair at large separation then has E_rel near V(∞) < 0 and is counted as a charmonium state. The density of states in Eq. (5) picks up a volume-proportional contribution from these mislabeled pairs. So the yield comparison to SHM and the v2 curves are not actually testing regeneration physics; they reflect an arbitrary threshold plus the tuned α_s = 0.7.\n\nWhat the paper does well: it implements the Langevin framework transparently, shows the box equilibrium distribution is internally consistent with the adopted density of states, and honestly flags the limited statistics. The charm-quark v2 as a proxy for open charm flow is a sensible application of the established machinery. The extension to charmonia, though, inherits the flawed criterion.\n\nOther soft spots are less severe. Five heavy-quark pairs per event is very low; the ALICE comparison is without error bars and for inclusive J/ψ while the simulation has no primordial charmonia. The α_s tuning to SHM means the agreement in Fig. 2b is not independent validation. The companion paper [4] is the right reference, but this proceedings alone does not give the reader enough detail to judge the transport coefficients.\n\nFor a proceedings note, the presentation is clear and the goal is worthwhile. But the central observable is undermined by the criterion. I would not cite this in its current form, and I would not send it to a referee as is. The right fix is to define the binding energy relative to the asymptotic real potential, rerun, and see whether the v2 result survives. If the authors do that, the framework may be worth revisiting.","headline":"The charmonium v2 is the only new result, but a bound-state criterion that counts arbitrarily separated pairs as charmonia undermines the central claim.","tokens_in":6450,"tokens_out":3996,"would_cite":false,"duration_ms":38588,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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…","keywords":["quark-gluon plasma","charmonium","elliptic flow","Langevin simulation","heavy quark diffusion","dissociation and regeneration","statistical hadronization model","relativistic heavy-ion collisions"],"falsifier":"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.","tokens_in":30,"feed_emoji":"⚛️","tokens_out":11622,"duration_ms":161516,"temperature":0.7,"pith_summary":"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.","feed_headline":"Langevin quark pairs reproduce J/psi flow at the LHC","feed_subtitle":"A screened Coulomb force plus Brownian motion reproduces LHC v2 data without primordial charmonia.","key_machinery":"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$.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"establishes the color-screening mechanism for quarkonium suppression that motivates the whole simulation strategy.","marker":"[1]"},{"why":"companion paper supplying the Langevin formalism, drag coefficient, and box-simulation method used here.","marker":"[4]"},{"why":"derives the complex heavy-quark potential from which the screened Coulomb force and drag are taken.","marker":"[5]"},{"why":"provides the statistical hadronization model yields used to calibrate the strong coupling and validate charmonium multiplicities.","marker":"[6]"},{"why":"provides the boost-invariant fireball parametrization that turns the box dynamics into a heavy-ion collision environment.","marker":"[9]"},{"why":"supplies the perturbative-QCD-tuned initial momentum distribution of heavy quarks used in the fireball simulations.","marker":"[13]"},{"why":"supplies the measured inclusive J/psi elliptic flow at sqrt(s_NN)=2.76 TeV used for the final comparison.","marker":"[15]"},{"why":"data collection version of the same J/psi elliptic flow measurement used in the comparison.","marker":"[16]"}],"fun_headline_variants":["J/psi v2 without primordial charmonia","Langevin recombines charm into flowing J/psi","Charm quarks flow, then recombine into J/psi","Recombination-driven J/psi v2 at LHC","No primordials? Langevin still gets J/psi flow"],"cache_read_input_tokens":8576,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["J/psi v2 without primordial charmonia","Langevin recombines charm into flowing J/psi","Charm quarks flow, then recombine into J/psi","Recombination-driven J/psi v2 at LHC","No primordials? Langevin still gets J/psi flow"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000309,"raw_usage":{"total_tokens":1760,"prompt_tokens":933,"completion_tokens":827,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":549,"completion_tokens_details":{"reasoning_tokens":741}},"tokens_in":549,"tokens_out":827,"duration_ms":8584,"temperature":1.0,"reasoning_tokens":741,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:07:51.072103+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Matsui and H","cited_arxiv_id":null,"evidence_quote":"establishes the color-screening mechanism for quarkonium suppression that motivates the whole simulation strategy."},{"cited_title":"Blaizot et al.,Heavy quark bound states in a quark–gluon plasma: Dissociation and recombination,Nuclear Physics A946 (2016) 49","cited_arxiv_id":null,"evidence_quote":"derives the complex heavy-quark potential from which the screened Coulomb force and drag are taken."},{"cited_title":"van Hees, M","cited_arxiv_id":null,"evidence_quote":"provides the boost-invariant fireball parametrization that turns the box dynamics into a heavy-ion collision environment."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the measured inclusive J/psi elliptic flow at sqrt(s_NN)=2.76 TeV used for the final comparison."},{"cited_title":"𝐽/𝜓 Elliptic Flow in Pb-Pb Collisions at√𝑠NN= 2.76 TeV","cited_arxiv_id":null,"evidence_quote":"data collection version of the same J/psi elliptic flow measurement used in the comparison."}],"review_version":1}