{"id":"3ff2c69f-934f-4c65-991c-ce6d787ac40b","arxiv_id":"2508.11897","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":6.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":3,"one_line_summary":"Heavy quarkonium equilibrates with the quark-gluon plasma in coupled Boltzmann and Langevin simulations, with timescales close to LHC fireball lifetimes.","lead":"Using a Boltzmann-plus-Langevin simulation in a static QGP box, the authors show heavy quarkonium can reach kinetic and chemical equilibrium with single heavy quarks, with equilibration times comparable to LHC QGP lifetimes. The work offers a dynamical route to explaining why statistical hadronization models work for charmonium, and a path to more realistic quarkonium transport phenomenology.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equilibration-time claim relies on unvalidated perturbative dipole rates near T_c; abstract-only prevents any check of the central numbers.","rationale":"The reader's weakest_assumption — the perturbative color-electric dipole amplitudes — is indeed the single most load-bearing concern. The quantitative claim of equilibration within QGP lifetime scales essentially with these reaction rates. Since the manuscript is abstract-only, this concern cannot be resolved now, and the UNVERDICTED verdict is appropriate. My concrete test would settle the matter by comparing the rates to a non-perturbative benchmark and rerunning with modified rates to gauge sensitivity.","tokens_in":679,"tokens_out":4461,"duration_ms":52492,"concrete_test":"Obtain the full manuscript and extract the thermally averaged dissociation and recombination cross sections used in the Boltzmann equation. Replace them with rates from a non-perturbative calculation (e.g., lattice QCD or pNRQCD with resummed couplings) and rerun the box simulation. If the resulting equilibration time changes by more than a factor of 2, the claim about LHC QGP lifetime is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim—that quarkonia reach full equilibrium with single heavy quarks within the LHC QGP lifetime—depends entirely on the dissociation and regeneration rates used in the Boltzmann–Langevin simulation. The abstract states these rates come from perturbative computations based on an effective color-electric dipole coupling. Near T_c, where quarkonium binding is marginal, such perturbative rates become unreliable. A moderate error of a factor of 2 in the rates would proportionally shift the equilibration time, potentially moving it outside the ~10 fm/c lifetime of the most central LHC fireballs. Because the full text is unavailable, there is no way to check the rate formulas, the temperature at which the box is run, or the sensitivity of the conclusion. This is the load-bearing point: the headline claim is only as strong as the accuracy of these perturbative amplitudes.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The abstract reports a simulation of heavy quarkonium equilibration in a static, homogeneous QGP box using a semi-classical Boltzmann equation for quarkonia coupled in real time to a Langevin simulation of single heavy-quark diffusion. Dissociation and regeneration amplitudes are said to come from perturbative computations based on effective color-electric dipole coupling of the quarkonium to thermal gluons. The central claim is that kinetic and chemical equilibrium between quarkonia and single heavy quarks is achieved in both yields and momentum distributions, with an equilibration time comparable to the QGP lifetime in central LHC heavy-ion collisions. This is offered as a dynamical explanation for the success of the statistical hadronization model. The abstract contains no equations, no numerical results, no error bars, and no comparison to data; the full text was not available for this review.","tokens_in":959,"tokens_out":3585,"duration_ms":43643,"significance":"If substantiated, the result would be significant for heavy-ion phenomenology: a microscopic transport calculation that dynamically drives quarkonia toward equilibrium with open heavy quarks within the fireball lifetime would provide a mechanistic underpinning for statistical hadronization and would distinguish between kinetic freeze-out and chemical equilibrium scenarios. The coupled Boltzmann-Langevin approach, combining open-heavy-quark diffusion with bound-state transport, is a promising methodology and is worth pursuing. However, on the basis of the abstract alone, the central claim is not checkable. No numerical timescale is quoted, no temperature or coupling is specified, and no sensitivity or validation exercise is described. The claim is falsifiable in principle—for instance, by comparing computed yields and momentum distributions to LHC data—so the potential value is clear, but the current submission does not allow a referee to verify any of it.","major_comments":[{"comment":"The central quantitative claim—that equilibration time is comparable to the QGP lifetime in central LHC collisions—is unsupported by any numerical value, temperature, box volume, or uncertainty. The word 'demonstrate' is not accompanied by a figure, table, or quantitative statement. Without these, a referee cannot assess whether the claim is correct or whether the timescale is within a factor of a few of the fireball lifetime.","section":"Abstract"},{"comment":"The dissociation and regeneration amplitudes are stated to be taken from perturbative computations based on effective color-electric dipole coupling, but no form of the amplitudes, no cutoff scale, and no domain of validity are given. Near the quarkonium dissociation temperature, where binding is marginal, perturbative dipole rates are known to be sensitive to the potential and to non-perturbative effects. A factor-of-two error in these rates would proportionally shift the equilibration time, potentially moving it outside the ~10 fm/c LHC QGP lifetime. The manuscript must quantify this sensitivity before the headline claim can be accepted.","section":"Abstract (scattering amplitudes)"},{"comment":"The real-time coupling of the Boltzmann transport for quarkonia and Langevin simulation for single heavy quarks is the core of the method, but no implementation details are given. In particular, the abstract does not state how the 'real-time fashion' coupling is realized (e.g., operator splitting, timestep, back-reaction of quarkonia on the open-heavy-quark distribution), how detailed balance is enforced, or how the box is initialized and thermalized. Without such information, the equilibration claim cannot be reproduced or checked.","section":"Abstract (simulation setup)"},{"comment":"No comparison is made to measured quarkonium yields or momentum spectra in heavy-ion collisions, nor to existing statistical hadronization or transport calculations. The abstract argues that the simulation 'provides a dynamical way of understanding' the phenomenological success of the statistical hadronization model, but without showing that the computed equilibrium yields and spectra match the data (or known model results), this conclusion is not established. A validation panel, even for a box calculation, is necessary.","section":"Abstract (validation)"}],"minor_comments":[{"comment":"'Comparable to the lifetime' is vague. The manuscript should state the numerical equilibration time and its criterion (e.g., 95% approach to equilibrium yields) and the QGP lifetime value used.","section":"Abstract"},{"comment":"'Heavy quarkonium with single heavy quarks' is ambiguous. Specify whether equilibrium is reached in the combined phase space of open and hidden heavy flavor, and define the chemical equilibrium condition.","section":"Abstract"},{"comment":"The phrase 'paves the way for realistic phenomenological applications' is promotional rather than informative; either remove it or indicate a specific next step.","section":"Abstract"}],"recommendation":"uncertain","confidential_remarks":"The review is necessarily provisional because the full manuscript was not available. The abstract alone is not sufficient to either accept or reject the work. If the full text is available, I would be glad to re-review with attention to the perturbative rate formulas, the coupled Boltzmann-Langevin implementation, and the numerical equilibration timescale. The editor may wish to verify that the submission includes the promised simulation details and results before sending it to a full review."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe abstract promises a genuinely new technical combination: real-time coupling of Langevin transport for single heavy quarks with Boltzmann transport for quarkonium, all in a static QGP box, with the aim of showing that the open and hidden sectors equilibrate together. If the simulation actually does that, it is a useful step toward a dynamical explanation of the statistical hadronization picture—a real contribution, not a routine extension. The paper also deserves credit for naming the key issue: the interplay between open and hidden heavy flavor in the equilibration process.\n\nThe soft spots are exactly where the stress-test note lands. The headline numerical claim—equilibration time comparable to the LHC QGP lifetime—rests entirely on the dissociation and regeneration rates, which come from perturbative color-electric dipole amplitudes. Near T_c those amplitudes are suspect; a factor-of-2 error in rates can shift the equilibration time by a comparable factor and push it outside the fireball lifetime. The abstract gives no equations, no numbers, no sensitivity analysis, and no comparison with known charmonium data, so we cannot check any of this. That is not a flaw in the paper's argument on its face—it is an absence of evidence in the abstract. The paper may well include the full rate formulas and tests; we just cannot see them.\n\nMy only other reservation is the list of free parameters (drag/diffusion coefficients, reaction rates, box temperature). If those are tuned to get the claimed equilibration, the result is less impressive; if they come from independent computations, fine. We cannot tell from the abstract.\n\nBottom line: the idea is worth referee time. Send it to someone who can verify the rate formulas and the coupling scheme, and ask the authors to report sensitivity of the equilibration time to rate uncertainties. I would not cite it without seeing the full text, but I would bring it to a reading group if someone gets a preprint with the numbers.\n\nRecommendation: accept for peer review, with serious scrutiny of the reaction rates and the equilibration-time sensitivity.","headline":"A new coupled Langevin–Boltzmann simulation of quarkonium in a QGP box, with a plausible but currently unverifiable claim about equilibration timescales that hinges on perturbative reaction rates near T_c.","tokens_in":1316,"tokens_out":2714,"would_cite":false,"duration_ms":32161,"reading_group":"maybe","serious_thinker":"unclear","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Heavy quarkonia reach full kinetic and chemical equilibrium with the quark-gluon plasma within the fireball lifetime, according to a coupled Boltzmann–Langevin simulation.","keywords":["heavy quarkonium","quark-gluon plasma","chemical equilibration","kinetic equilibration","Boltzmann transport","Langevin diffusion","statistical hadronization","charmonium"],"falsifier":"Compute the same dissociation and regeneration rates at temperatures around 1.5–2.0 T_c directly from the color-electric dipole amplitudes and compare with the rates built into the simulation; a sizable mismatch would invalidate the equilibration time. Experimentally, if the J/ψ yield in the most central Pb-Pb collisions at the LHC lay significantly below the statistical-hadronization expectation based on measured open-charm production, or if its momentum spectrum remained strongly anisotropic after full thermalization, the claimed equilibration would be ruled out.","tokens_in":665,"feed_emoji":"⚛️","tokens_out":3765,"duration_ms":43560,"temperature":0.7,"pith_summary":"This paper tries to establish that heavy quarkonia—bound states such as charmonium—reach both kinetic and chemical equilibrium with the surrounding quark-gluon plasma within the plasma's lifetime in central heavy-ion collisions at LHC energies. The authors simulate the process with a semi-classical Boltzmann transport for the quarkonium coupled in real time to a Langevin simulation of single heavy-quark diffusion. They find that yields and momentum distributions of the bound states converge to the values expected in statistical hadronization, offering a dynamical explanation for why that phenomenological model has worked. If true, the result would mean quarkonium production in such collisions is not a simple superposition of initial production and survival but is actively equilibrated by the medium.","feed_headline":"Quarkonia equilibrate within the quark-gluon plasma lifetime","feed_subtitle":"Simulation couples quark diffusion with bound-state transport to explain charmonium yields at the LHC.","key_machinery":"The coupled real-time transport: a Boltzmann equation for heavy quarkonium dissociation and regeneration (with LO and NLO reaction rates from effective color-electric dipole coupling to thermal gluons) and a Langevin equation for single heavy-quark diffusion, iterated together so that the heavy-quark phase space feeds the regeneration rates while quarkonium decay feeds back into the open-heavy-quark distribution.","core_discovery":"The central claim is that, in a static and homogeneous box of quark-gluon plasma, heavy quarkonium equilibrates with the single heavy quarks through repeated dissociation and regeneration reactions, reaching chemical equilibrium in both the bound-state yields and their momentum spectra. The equilibration time extracted from the simulation is comparable to the lifetime of the QGP fireball in the most central Pb-Pb collisions at the LHC, which is why the statistical hadronization model reproduces charmonium production so well. The dynamics hinge on the coupling between the open heavy sector (single heavy quarks diffusing through the medium) and the hidden sector (the bound states), so that reg","pith_inferences":["The same coupled transport could be extended to bottomonium; because of the larger b-quark mass and smaller dipole coupling, the equilibration time should be longer, predicting incomplete bottomonium equilibration in the same collisions.","The result suggests that statistical hadronization is not an independent principle but an emergent attractor: any plasma hot enough and long-lived enough will drive quarkonium yields to equilibrium regardless of the initial production ratio.","A testable scaling follows: the equilibration time should be inversely related to the heavy-quark diffusion coefficient and the in-medium dissociation rate, so measurements of quarkonium momentum broadening could be used to pin down these transport coefficients."],"forward_implications":["Central Pb-Pb collisions at LHC should yield charmonium abundances and momentum spectra close to the statistical-hadronization equilibrium values set by the open-charm density.","The degree of quarkonium equilibration should vary with collision centrality and system size, since the QGP lifetime in smaller systems may be shorter than the calculated equilibration time.","Quarkonium momentum distributions should be thermalized with the plasma, so reduced anisotropies relative to a purely primordial production scenario would be expected.","The open and hidden heavy flavor sectors cannot be modeled independently: regeneration and dissociation couple the two, and any realistic transport description must include both."],"supporting_citations":[],"fun_headline_variants":["Heavy quarkonia reach equilibrium inside QGP fireball","Simulation: quarkonia equilibrate with heavy quarks in QGP","Quarkonium equilibration time matches QGP lifetime at LHC","Coupling open and hidden heavy sector drives quarkonia equilibration","Quarkonium equilibration explains charmonium yields at LHC"],"cache_read_input_tokens":2816,"weakest_assumption_plain":"The predicted equilibration time rests on perturbative scattering amplitudes for quarkonium dissociation and regeneration based on an effective color-electric dipole coupling to thermal gluons; if those amplitudes are inaccurate near the dissociation temperature, the computed rates and timescales would shift.","fun_headline_variants_meta":{"raw":{"variants":["Heavy quarkonia reach equilibrium inside QGP fireball","Simulation: quarkonia equilibrate with heavy quarks in QGP","Quarkonium equilibration time matches QGP lifetime at LHC","Coupling open and hidden heavy sector drives quarkonia equilibration","Quarkonium equilibration explains charmonium yields at LHC"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000614,"raw_usage":{"total_tokens":2684,"prompt_tokens":728,"completion_tokens":1956,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":472,"completion_tokens_details":{"reasoning_tokens":1866}},"tokens_in":472,"tokens_out":1956,"duration_ms":13913,"temperature":1.0,"reasoning_tokens":1866,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T19:40:47.214524+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the same dissociation and regeneration rates at temperatures around 1.5–2.0 T_c directly from the color-electric dipole amplitudes and compare with the rates built into the simulation; a sizable mismatch would invalidate the equilibration time. Experimentally, if the J/ψ yield in the most central Pb-Pb collisions at the LHC lay significantly below the statistical-hadronization expectation based on measured open-charm production, or if its momentum spectrum remained strongly anisotropic after full thermalization, the claimed equilibration would be ruled out.","supporting_citations":[],"review_version":1}