{"id":"de557f6d-c6d1-4acc-b43b-d61b2cebfa94","arxiv_id":"2603.04674","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A unified T-matrix transport framework for charm-quark diffusion and charmonium kinetics with off-shell spectral functions recovers the statistical equilibrium limit and gives fair LHC charmonium R_AA data.","lead":"This paper couples the diffusion of charm quarks to the dissociation and regeneration of charmonium states using one shared nonperturbative interaction, instead of treating open and hidden charm separately. It then applies the coupled framework to charmonium yields in lead-lead collisions at the LHC and finds fair agreement with data.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (21)'s off-shell regeneration rate inserts a thermal weight e^{-[ω'-ε_c(p')]/T} with no derivation; if that weight is wrong, the off-equilibrium regeneration and peripheral R_AA change.","rationale":"The reader's weakest assumption is exactly the ad hoc off-shell weighting in Eq. (21), and I agree. The paper's headline claims—'for the first time... same interactions... equilibrium limit recovered'—are supported for the thermal case: Eq. (20) follows from detailed balance with Eq. (16), and Fig. 10 is an internal consistency check. The only place where a new physical input enters without derivation is the off-equilibrium off-shell regeneration rate. The factor is chosen to satisfy the equilibrium limit, but that does not make it unique or correct for non-thermal charm quarks. Since the paper's phenomenological conclusions about off-equilibrium regeneration and peripheral collisions depend directly on this factor, the concern is load-bearing. It is, however, addressable: a first-principles derivation or a sensitivity study with a plausible alternative weighting would settle it. Therefore I would keep the reader's CONDITIONAL verdict rather than escalate to REJECT. No issue with the authors' conduct or the overall framework's novelty arises.","tokens_in":21815,"tokens_out":13908,"duration_ms":133514,"concrete_test":"Re-derive the 3→2 regeneration rate from the same half-off-shell T-matrix and spectral functions used for α (lower panel of Fig. 2) without imposing the factor e^{-[ω'-ε_c]/T}; e.g., sample the incoming charm energy ω' from ρ_c(ω',p') at each Langevin momentum p'_c and use f_c(ω',p',t)=f_c^{on}(p'_c,t) ρ_c(ω',p') (or, better, solve the off-shell Boltzmann equation for f_c(ω',p',t)). Compare the resulting β with Eq. (21) at T=250 MeV and a few P_Ψ. If the ratio deviates by >20%, recompute the 40–90% J/ψ R_AA(p_T) panel of Fig. 13 with the derived β; a shift of the low-p_T regenerated yield beyond the shadowing band would confirm the concern is load-bearing.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central phenomenological novelty is the off-equilibrium regeneration rate, Eq. (21), which replaces the thermal charm occupation n_c(ω') in Eq. (20) by f_c(p'_c,t) exp(-[ω'-ε_c(p'_c)]/T(t)). This exponential is not derived from the Kadanoff-Baym or Boltzmann dynamics used for the other rates; it is inserted so that in the thermal limit f_c -> exp(-ε_c/T) the product reduces to exp(-ω'/T), recovering Eq. (20). The text itself says only 'It turns out that...' (Sec. 3.2). For a genuinely non-thermal charm population, the off-shell energy dependence of f_c(ω',p',t) is an independent dynamical quantity that the on-shell Langevin simulation does not provide; the spectral function ρ_c(ω',p') already encodes the probability of off-shell fluctuations, so multiplying by an additional thermal weight can double-count medium effects. The suppression of the equilibrium limit in peripheral collisions (Fig. 11) is a direct consequence of this prescription. If the correct off-shell weighting is different—e.g., governed by the same relaxation kernel as A(p), or absent in the large-off-shellness regime—the magnitude and p_T dependence of regenerated charmonia would change. This is a missing derivation, not a numerical bug, and it is the most load-bearing assumption in the paper.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a coupled transport framework for open and hidden charm in the quark–gluon plasma formed in Pb–Pb collisions at the LHC. Charm-quark diffusion is simulated with relativistic Langevin dynamics, with the input relaxation rate A(p) computed from a thermodynamic T-matrix constrained by lattice-QCD Wilson-line correlators and including off-shell charm and light-parton spectral functions. Charmonium dissociation rates α for J/ψ, χ_c, and ψ(2S) are computed from the same T-matrix in a quasifree approximation, and the regeneration rate β is constructed from α through a detailed-balance relation (Eq. 17), which fixes the charmonium equilibrium limit to f_eq = γ_c² d_Ψ exp(−E_Ψ/T) (Eq. 26). The new element is the off-equilibrium extension of the regeneration rate (Eq. 21), in which the thermal charm occupation n_c(ω′) is replaced by the Langevin-evolved distribution f_c(p′) multiplied by an additional thermal off-shell weight exp[−(ω′−ε_c(p′))/T]; this is the paper's mechanism for controlling regeneration from non-thermalized charm. The framework is applied to Pb–Pb at √s_NN = 5.02 TeV using a schematic blast-wave fireball, producing R_AA as a function of centrality and p_T for J/ψ, χ_c, and ψ(2S), with fair agreement with ALICE data.","tokens_in":22302,"tokens_out":22567,"duration_ms":213129,"significance":"If the construction holds, this is a significant step: for the first time the same nonperturbative heavy-light interaction and the same in-medium spectral functions determine both open-charm diffusion and charmonium kinetics in a coupled Langevin–Boltzmann simulation, without K-factors. Strengths include the transparent derivation of the thermal regeneration rate (Eq. 20) from Eq. (16) via the identity n(E)=e^{−E/T}(1±n(E)); the internally consistent equilibrium limit, which is non-trivial to realize with broad spectral functions; the absence of tuned transport parameters; and falsifiable predictions for χ_c and ψ(2S) R_AA(p_T) and centrality dependence that can be confronted with LHC data. The improved peripheral low-p_T description relative to earlier blast-wave treatments is a concrete, credit-worthy result. However, the central novelty—the off-equilibrium regeneration rate in Eq. (21)—rests on an underexplained ansatz, the equilibrium-limit 'recovery' is enforced by construction, and the theoretical uncertainties on the rates are unquantified. These issues are addressable but currently limit the strength of the phenomenological claims.","major_comments":[{"comment":"The central new ingredient is introduced by replacing n_c(ω′) with f_c(p′_c) exp[−(ω′−ε_c(p′_c))/T], justified only by 'It turns out that…' and the requirement of the correct equilibrium limit. Since ρ_c(ω′,p′) already weights off-shell energies, the exponential is a further dynamical assumption: off-shell fluctuations are thermally populated at temperature T even when the momentum distribution f_c is far from equilibrium. The peripheral suppression of the equilibrium limit (dashed curves in Fig. 11) and the low-p_T regenerated R_AA in Figs. 13–15 inherit this prescription; no derivation (e.g., from a Kadanoff–Baym G^< ansatz) and no sensitivity estimate are given. Please derive the weight or test alternatives (no extra weight; on-shell projection; n(ω′)/n(ε(p′)) ratio) and report the effect on Fig. 11 and the peripheral R_AA. Note also that Eq. (23) yields exp(−ω′/T), whereas Eq. (20) c","section":"Sec. 3.2, Eq. (21)"},{"comment":"The equilibrium limit f^eq = β/α = γ_c² d_Ψ exp(−E_Ψ/T) is enforced by construction: β is defined in Eq. (17) as γ_c² d_Ψ exp(−E_Ψ/T) α, and Eq. (26) then follows from the stationary limit of Eq. (8). The comparison with the statistical-model curve in Fig. 10 therefore validates the internal consistency of the off-shell phase-space integrals, not the statistical-model equilibrium limit independently. Moreover, the two curves share the same charm fugacity from Eq. (27), the same masses, and the same binding energies, so the agreement is not an independent check. The abstract's wording ('the equilibrium limit of the statistical model is recovered') should be softened; the non-trivial content—consistency of the off-shell integral structure with the identity n=e^{−E/T}(1±n)—should be stated as such.","section":"Sec. 3.3, Eqs. (17),(26); Fig. 10"},{"comment":"The hadronic-phase regeneration is bridged by rescaling the QGP off-equilibrium suppression factor and holding it fixed throughout the hadronic evolution: 'we utilize the reduction of the equilibrium limit from the QGP phase at T_H ... and use the same rescaling factor throughout the hadronic medium evolution.' This is an ad hoc assumption; charm-quark distributions continue to evolve (or freeze) through hadronization, and the final R_AA values—especially for peripheral and semi-central collisions, where the hadronic phase occupies a sizeable fraction of the fireball lifetime—inherit this choice directly. A sensitivity estimate (e.g., varying the rescale factor between its central value and unity, or matching to a hadronic Langevin evolution) is needed to establish that the reported centrality trends do not depend on this bridge.","section":"Sec. 4.2 (hadronic-phase matching); Fig. 11"},{"comment":"The central transport inputs (A(p), α_Ψ, β, m*_c(T)) are presented as single curves with no propagated uncertainties from the WLC lattice constraints or the T-matrix fits. The shadowing bands in Figs. 12–15 are included, but the dominant theoretical uncertainty likely lies in the rates themselves. Since the paper's phenomenological statements are comparative ('capture the measured centrality and momentum dependence fairly well', 'slight underestimation', 'overshoot'), an estimate of the rate uncertainty, or a scan over, e.g., the WLC scenarios or the correlation-volume parameters (as done for κ in Sec. 3.3), would materially strengthen those conclusions.","section":"Secs. 2.2, 3.2, 4.3"}],"minor_comments":[{"comment":"The notation f_c(p′_c, T(t)) is misleading: the Langevin distribution has time, not temperature, as its second argument. Write f_c(p′_c, t), and state the normalization convention (per state vs total number) so that the γ_c² prefactor does not double-count charm number.","section":"Sec. 3.2, Eq. (21)"},{"comment":"Please clarify the volume factors: the canonical open-charm term carries a V_FB prefactor while the Bessel argument uses V_co. If this is the intended Zhao–Rapp convention, a parenthetical justification would prevent confusion.","section":"Eq. (27)"},{"comment":"The interference factor '1−e ik·r' should be typeset as 1−e^{ik·r}, and its sign/real-part convention relative to Ref. [19] should be stated explicitly.","section":"Eq. (16); Fig. 2"},{"comment":"The temperature list ('700 500 352 293 251 195 174') is printed below the legend boxes; a color bar or direct labeling of the curves would be clearer.","section":"Fig. 1 caption"},{"comment":"The phrase 'up to about p_T ≃ 5(2) GeV' is unclear; specify which states correspond to which values.","section":"Sec. 5"}],"recommendation":"major_revision","confidential_remarks":"This is a serious contribution from an established group, and the coupled framework is potentially important. The stress-test concern about Eq. (21) lands: the off-shell weight is load-bearing for the off-equilibrium regeneration and the peripheral equilibrium-limit suppression, and it is introduced without derivation. However, the ansatz is not unreasonable—it can be interpreted as a conditional thermal distribution of off-shellness at fixed momentum—so this is fixable with a derivation/appendix and sensitivity tests rather than a fatal error. I would not support rejection. The equilibrium-limit claim should also be reframed as internal consistency rather than independent validation. The phenomenological part is explicitly preliminary (schematic fireball), which is acceptable if the central claims are hedged accordingly."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper is a credible extension of the Rapp-group T-matrix program, not a breakthrough but a solid within-subfield step. The genuinely new piece is coupling Langevin charm diffusion and Boltzmann charmonium kinetics with the same nonperturbative interaction and spectral functions. That is the first time this is done in a single framework, and the detailed-balance derivation of the equilibrium regeneration rate (Eq. 20) is clean.\n\nWhat it does well: the paper is honest about its schematic fireball, the fits to lattice-QCD charm susceptibilities are neat, and the comparison to LHC data is fair, with no K-factors. The equilibrium limit is recovered when charm quarks thermalize, which is the right sanity check. The centrality and pT dependence of J/psi and excited states are reproduced reasonably well, and the authors flag their own deficits in the intermediate-pT region.\n\nSoft spots: Eq. (21) is the load-bearing piece. The thermal off-shell weighting factor exp(-[omega'-epsilon_c(p'_c)]/T) is inserted \"so that\" the equilibrium limit is recovered; it is not derived from the Kadanoff-Baym or Boltzmann dynamics used elsewhere. For a genuinely off-equilibrium charm spectrum, the off-shell energy dependence of f_c is a separate dynamical quantity, and the on-shell Langevin simulation does not provide it. The peripheral R_AA suppression depends on this factor. The stress-test note is right: this is a missing derivation, not a numerical bug. I would want that derivation, or at least a sensitivity study on the weighting, before accepting the off-equilibrium predictions.\n\nAlso: there are no quoted uncertainties on the transport rates, the medium is a schematic fireball, and the validation against the statistical model in Fig. 10 shares the same fugacity and density inputs, so it is not an independent check. These are addressable issues rather than fatal flaws. The central equilibrium-limit derivation holds up, and the framework is a genuine step forward.\n\nVerdict: I would send this to a serious referee. It deserves careful review, with special pressure on the derivation of Eq. (21). If that factor can be justified or shown to be phenomenologically robust, the paper is worth publishing; if not, the off-equilibrium conclusions would need to be walked back. Worth engaging with.","headline":"Same-T-matrix coupling of open and hidden charm is a real step; the off-shell off-equilibrium regeneration factor in Eq. (21) is imposed, not derived, and carries the peripheral R_AA conclusions.","tokens_in":22705,"tokens_out":1887,"would_cite":true,"duration_ms":19281,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Charm and charmonium transport can be described by one nonperturbative interaction, the thermodynamic T-matrix, coupled through Langevin and Boltzmann equations, with in-medium spectral functions providing a consistent link between charm qu","keywords":["charmonium transport","quark-gluon plasma","charm quark diffusion","T-matrix","Langevin dynamics","Boltzmann equation","regeneration","Wilson-line correlators"],"falsifier":"Compute the same off-shell 3-to-2 regeneration rate from a fully quantum nonequilibrium Green's function treatment and compare its thermalized limit to the statistical-model equilibrium limit: if the detailed-balance form is not recovered with the imposed weighting factor, the central phenomenological result fails. Alternatively, a high-precision measurement of peripheral J/psi R_AA at low pT, where off-equilibrium suppression is largest, would test the prediction.","tokens_in":21676,"feed_emoji":"⚛️","tokens_out":4644,"duration_ms":45935,"temperature":0.7,"pith_summary":"Charm quarks and charmonia (bound charm–anticharm states) have usually been treated with separate models in heavy-ion collisions: one for how charm quarks diffuse, another for how J/psi and its excited states dissociate and reform. This paper argues that a single nonperturbative interaction—the thermodynamic T-matrix with in-medium spectral functions constrained by lattice QCD—can govern both. It feeds the same interaction into a Langevin simulation of charm diffusion coupled to a Boltzmann equation for charmonium kinetics, and shows that the charmonium equilibrium limit is recovered when charm quarks thermalize. The framework then reproduces the measured centrality and transverse-momentum dependence of charmonium production in Pb-Pb collisions at the LHC without tuned rates. A sympathetic reader would care because this is a step toward one consistent description of all heavy-flavor probes in the strongly coupled quark-gluon plasma.","feed_headline":"One interaction now drives both charm flow and J/psi regeneration","feed_subtitle":"Coupled T-matrix transport ties charm diffusion to charmonium regeneration, reproducing LHC data without tuned rates.","key_machinery":"The central object is the thermodynamic T-matrix: a nonperturbative scattering amplitude that encodes how a heavy charm quark interacts with light quarks, antiquarks, and gluons, computed self-consistently with in-medium spectral functions and a heavy-quark potential constrained by lattice-QCD Wilson-line correlators. From this single input the paper computes the charm-quark relaxation rate, which feeds the Langevin equation, and the charmonium dissociation and regeneration rates, which feed the Boltzmann equation. The regeneration rate is the mechanism that connects the two sectors: it takes the time-dependent off-equilibrium charm distribution from the Langevin simulation and converts corr","core_discovery":"The paper constructs a closed transport description in which charm-quark and charmonium transport coefficients are computed from the same underlying heavy-light interaction: the nonperturbative T-matrix with broad in-medium spectral functions and a potential constrained by Wilson-line correlators from lattice QCD. It derives the charmonium regeneration rate using off-shell spectral functions, imposes detailed balance with the dissociation rate, and extends the rate to off-equilibrium charm-quark distributions from Langevin dynamics by inserting a thermal off-shell weighting factor. In the thermalized limit this factor reduces to the equilibrium detailed-balance relation and reproduces the st","pith_inferences":["A consequence the paper leaves implicit: if one interaction governs all heavy-flavor transport, then simultaneous measurements of D mesons, J/psi, chi_c, psi(2S), and bottomonia could be fit globally to pin down the in-medium potential, turning qualitative agreement into a quantitative constraint.","The off-shell weighting factor in the regeneration rate is an assumption, not a derivation; a fully quantum nonequilibrium Green's function calculation of the same process would either confirm it or replace it, and the same prescription would then carry over to bottomonium.","The centrality-dependent suppression of regeneration implies that the effective charmonium fugacity is not a constant; comparing psi(2S)/J/psi and chi_c/J/psi ratios across centrality could expose this off-equilibrium effect directly."],"forward_implications":["The same T-matrix interaction now links open-charm, e.g. D-meson, phenomenology to charmonium suppression, so both become tests of one in-medium QCD force without K-factors.","When charm quarks thermalize, the ratio of regeneration to dissociation rates reproduces the statistical-model charmonium equilibrium limit, validating the rate construction.","Incomplete charm thermalization lowers the effective equilibrium limit relative to the statistical model, most strongly in peripheral collisions, and flattens the regenerated pT spectra.","Charmonium states are predicted to survive as broad resonances up to higher temperatures (about 500, 350, and 250 MeV for J/psi, chi_c, and psi(2S)), extending the window in which regeneration can occur.","Primordial charmonia are essentially eliminated in central collisions up to pT around 5 GeV, with the final yield dominated by regeneration."],"fun_headline_variants":["One T-matrix interaction now drives charm flow and J/psi regeneration","Coupled charm-charmonium transport from same heavy-light force","Self-consistent transport ties charm diffusion to charmonium kinetics","Lattice-constrained T-matrix unifies charm and charmonium dynamics"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the off-shell regeneration rate can be modeled by taking the Langevin charm distribution and multiplying by an extra thermal factor exp[-[omega'-epsilon_c]/T] chosen to enforce detailed balance; this factor is not derived from the dynamics, and the predicted suppression of regeneration in peripheral collisions depends on it.","fun_headline_variants_meta":{"raw":{"variants":["One T-matrix interaction now drives charm flow and J/psi regeneration","Coupled charm-charmonium transport from same heavy-light force","Self-consistent transport ties charm diffusion to charmonium kinetics","Lattice-constrained T-matrix unifies charm and charmonium dynamics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000204,"raw_usage":{"total_tokens":1222,"prompt_tokens":739,"completion_tokens":483,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":483,"completion_tokens_details":{"reasoning_tokens":407}},"tokens_in":483,"tokens_out":483,"duration_ms":4870,"temperature":1.0,"reasoning_tokens":407,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T18:45:42.854423+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the same off-shell 3-to-2 regeneration rate from a fully quantum nonequilibrium Green's function treatment and compare its thermalized limit to the statistical-model equilibrium limit: if the detailed-balance form is not recovered with the imposed weighting factor, the central phenomenological result fails. Alternatively, a high-precision measurement of peripheral J/psi R_AA at low pT, where off-equilibrium suppression is largest, would test the prediction.","supporting_citations":[],"review_version":1}