{"id":"f1382f0f-c98f-4648-b070-b5f3c39cf0fb","arxiv_id":"2607.08707","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":4,"one_line_summary":"Including a phenomenological heavy-light quark potential in the coalescence model enhances the heavy-quark coalescence probability to near unity at low momentum without ad hoc normalization.","lead":"This paper shows that including a heavy-quark interaction potential in the coalescence model for heavy-ion collisions raises the coalescence probability to near unity, resolving a known deficit. It matters because it removes the need for ad hoc normalization factors in heavy-quark hadronization models.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"The meson radius entering the Wigner function is treated as a free parameter despite being determined by the same potential used elsewhere in the calculation; self-consistent determination could move the coalescence probability away from unity.","rationale":"The reader's verdict of CONDITIONAL is appropriate, and the concern identified here sharpens rather than overturns it. The paper's central claim—that including the potential eliminates the need for ad hoc normalization—rests on the coalescence probability reaching near unity. This in turn depends on the meson radius being small. The paper treats this radius as a free parameter, but it is not free: it is determined by the same potential the paper already calibrates. This is a genuine internal inconsistency, not merely a sensitivity to an external unknown. If the self-consistent radius is small (≲0.6 fm), the claim strengthens considerably. If it is larger, the paper's main result weakens, and the need for a normalization factor may persist. The paper would be substantially strengthened by performing this self-consistent calculation, which requires no new physics input—only using the wave function already computed in Eq. (6) to determine the Wigner-function width. The in-medium screening analysis (Sec. IV) is secondary and does not change this core issue. The paper remains a solid and well-motivated contribution, but the central claim is not fully established without the self-consistent radius check.","tokens_in":10415,"tokens_out":2144,"duration_ms":101902,"concrete_test":"Solve the Schrödinger equation (Eq. 6) with the fitted vacuum potential (α_s ≈ 1.0, m_c = 1.14 GeV, σ = 0.445 GeV², µ = 0.35√σ) to obtain the 1S wave function ψ(r,0). Compute ⟨r²⟩ = ∫r²|ψ(r)|² d³r self-consistently. Insert this radius into Eq. (10) to fix σ, then recompute the coalescence probability from Eq. (15). If the self-consistent D-meson radius exceeds ~0.6 fm, the total charm coalescence probability at m_q ≈ 0.45 GeV falls below ~0.9, weakening the 'approaches unity' claim.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper uses the heavy-light potential in two places: (1) in the Schrödinger equation (Eq. 6) to fit meson masses, which yields a wave function ψ(r,T) with a definite spatial extent, and (2) in the thermal distribution of light antiquarks (Eq. 13) to enhance coalescence. The Wigner function (Eq. 9) requires a width σ related to the meson radius via Eq. (10), but this radius is treated as a free parameter (0.5–1.0 fm) rather than computed from the wave function already obtained by solving Eq. (6). This is an internal inconsistency: the same potential that reproduces meson masses also predicts ⟨r²⟩ = ∫r²|ψ(r)|² d³r, which should fix σ. The central claim that coalescence 'approaches unity' holds primarily at the lower end of the radius range (0.5 fm, upper band in Figs. 3–4). If the self-consistent radius from the fitted potential turns out to be 0.7–1.0 fm, the probability drops below unity even with the potential included, undermining the paper's main conclusion that no ad hoc normalization is needed. The reader correctly identified the radius sensitivity but did not flag that the radius is not independent of the potential—it is predicted by it.","agreement_with_reader":"partial"},"referee_report":{"model":"glm-5.2","summary":"This manuscript investigates the role of the heavy-light quark potential in heavy-quark coalescence at the QGP phase boundary. The authors construct a phenomenological potential (Coulomb + string + spin-spin) fitted to reproduce the vacuum masses of D, D*, B, and B* mesons. When this potential is included in the thermal distribution of light antiquarks (Eq. 13), the attractive interaction enhances the local density of light antiquarks around a static heavy quark, increasing the coalescence probability (Eq. 15) toward unity without an ad hoc normalization factor. The authors also study in-medium modifications via a temperature-dependent screening mass, finding that the coalescence probability remains near unity as long as the screening parameter satisfies mu/sqrt(sigma) < ~0.5, which corresponds to a ~60 MeV in-medium mass shift of D and B mesons.","tokens_in":10866,"tokens_out":1243,"duration_ms":176666,"significance":"The problem addressed is well-motivated: within the quasi-particle model framework constrained by the lattice-QCD equation of state, the total heavy-quark coalescence probability at T_c falls significantly below unity, which is physically problematic since a static heavy quark must hadronize via coalescence. The paper provides a concrete mechanism—potential-induced spatial clustering of light antiquarks—to resolve this without ad hoc rescaling. The potential is independently anchored to vacuum spectroscopy (Fig. 1), and the in-medium constraint (mu/sqrt(sigma) < 0.5 implying mass shifts < ~60 MeV) is a falsifiable prediction that can be compared against QCD sum-rule and effective-Lagrangian calculations. The finding that the coalescence probability is more sensitive to the screening mass than the meson mass itself is a useful phenomenological observation.","major_comments":[{"comment":"The heavy-meson radius entering the Wigner function width sigma (Eq. 10) is treated as a free parameter in the range 0.5-1.0 fm, but the same potential used in Eq. (13) is also solved via the Schrodinger equation (Eq. 6) to obtain the meson wave function psi(r, T). This wave function has a definite spatial extent, <r^2> = integral of r^2 |psi(r)|^2 d^3r, which should fix sigma self-consistently. The central claim that the coalescence probability 'approaches unity' holds primarily at the lower end of the radius range (0.5 fm, which corresponds to the upper boundary of the bands in Figs. 3-4 after the potential is included). If the self-consistent radius from the fitted potential is in the 0.7-1.0 fm range, the probability falls below unity even with the potential, which would undermine the conclusion that no ad hoc normalization is needed. The authors should either compute <r^2> from the波","section":"Sec. III, Eqs. (9), (10), and (6)"}],"minor_comments":[{"comment":"The caption states 'blue' and 'red' for the without- and with-potential curves, but the figure appears to use a different color scheme in the rendered version. Please verify consistency.","section":"Fig. 2 caption"},{"comment":"The notation D_3 is used for the spin degeneracy of the D meson, but D_3 was also used in Eq. (8) for the color-spin degeneracy factor of particle 3 (the produced hadron). Clarify whether these are the same quantity or distinct.","section":"Sec. III, Eq. (15)"},{"comment":"The sentence 'For pseudoscalar D meson, the spin degeneracy factor (D_3 in Eq. (15) equals 1' has an unmatched parenthesis. Please fix.","section":"Sec. III, below Eq. (18)"},{"comment":"The figure shows meson masses vs. mu/sqrt(sigma), but the range plotted (0.35-0.70) extends beyond the range discussed in the text (0.35-0.5). It would help to mark the mu/sqrt(sigma) = 0.5 value explicitly on the plot.","section":"Sec. IV, Fig. 5"},{"comment":"The sentence 'In this work, the heavy-light quark potential that reproduces the masses of both pseudoscalar and vector heavy mesons.' is grammatically incomplete; a verb is missing.","section":"Sec. V"},{"comment":"The phrase 'At the same time, resulting in a shallower potential...' is a sentence fragment. Please integrate with the preceding sentence.","section":"Sec. IV, last paragraph"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern about the self-consistency of the meson radius is, in my assessment, the single most important issue. The paper uses the potential in two places—for the thermal distribution (Eq. 13) and implicitly for the Wigner function width (Eq. 10 via the radius)—but does not close the loop by computing the radius from the same wave function. If the self-consistent radius turns out to be large, the main conclusion weakens considerably. This is fixable: the authors already solve Eq. (6), so computing <r^2> is straightforward. If the self-consistent radius is in the 0.5-0.7 fm range, the paper's conclusions are strengthened; if not, the claims need to be qualified. Either outcome is publishable, but the current treatment leaves a gap that a knowledgeable referee will flag."},"author_rebuttal":null,"desk_editor":{"model":"glm-5.2","letter":"The key result here is that including a heavy-light quark potential in the coalescence formalism boosts the total heavy-quark coalescence probability near $T_c$ to something close to unity, removing the need for the ad hoc normalization factor that previous models required. That's a real fix to a known problem in the field, and the authors deserve credit for it. The potential is well-motivated (lattice-QCD-inspired form, fitted to reproduce $D$, $D^*$, $B$, $B^*$ masses), and the medium modification via the screening mass is handled sensibly. The constraint that $P_{coal} ≈ 1$ implies $μ/√σ ≲ 0.5$ near $T_c$, corresponding to ~60 MeV in-medium mass shifts, is a useful, falsifiable prediction consistent with other approaches.","headline":"Solid idea, but the meson radius should be computed from the potential they already solved — not treated as a free parameter.","tokens_in":11143,"tokens_out":242,"would_cite":true,"duration_ms":101936,"reading_group":"no","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["25.75.-q","12.38.Mh","12.39.Jh"],"model":"glm-5.2","headline":"Attractive potential solves heavy-quark coalescence puzzle","keywords":["heavy quark","coalescence","hadronization","quark-gluon plasma","heavy-light potential","color screening","meson mass","charm"],"falsifier":"If lattice QCD or experimental determinations of the D- or B-meson radius near the phase transition yield values significantly larger than 0.5 fm, or if the in-medium screening mass exceeds the threshold identified here, the coalescence probability would fall substantially below unity, undermining the claim that the potential alone resolves the deficit.","tokens_in":10722,"feed_emoji":"🧩","tokens_out":1014,"duration_ms":154947,"temperature":0.7,"pith_summary":"Heavy quarks produced in heavy-ion collisions must hadronize through coalescence with light quarks as the plasma cools. Standard coalescence models, treating light quarks as uniformly distributed, predict a total coalescence probability far below the expected value of unity for a stationary heavy quark. This paper introduces a phenomenological heavy-light quark potential, fitted to reproduce the vacuum masses of D, D*, B, and B* mesons, into the coalescence formalism. The attractive potential concentrates light antiquarks around the heavy quark in coordinate space, naturally boosting the coalescence probability to near unity without any ad hoc normalization factor. The authors further show that this conclusion survives moderate in-medium color screening of the potential, but breaks down if screening is too strong, setting an upper bound on the in-medium mass shift of heavy mesons near the phase transition temperature.","feed_headline":"Attractive potential solves heavy-quark coalescence puzzle","feed_subtitle":"Including a realistic quark potential naturally raises the coalescence probability to unity, no ad hoc fix needed.","key_machinery":"The mechanism is a coordinate-space density enhancement: an attractive potential well, parameterized via a Coulomb-plus-string form with a spin-spin interaction, draws light antiquarks toward the heavy quark. The coalescence probability is computed using a Gaussian Wigner function for the meson wavefunction, integrated over the light-quark phase-space distribution modified by the potential. The potential itself is calibrated to reproduce vacuum meson masses, and medium effects enter through a temperature-dependent screening mass that controls the depth and range of the interaction.","core_discovery":"The central finding is that including a realistic attractive heavy-light quark potential in the coalescence model resolves the long-standing deficit in the total heavy-quark coalescence probability at low momentum. The potential enhances the local density of light antiquarks near the heavy quark, raising the probability from well below unity to approximately unity for both charm and bottom quarks. This result holds as long as the in-medium screening mass parameter does not exceed roughly half the square root of the string tension, which translates to a maximum mass reduction of about 60 MeV for D and B mesons near the critical temperature.","pith_inferences":["If the true meson radius is closer to the upper end of the considered range (1.0 fm), the coalescence probability falls below unity even with the potential, suggesting that tighter experimental or lattice constraints on meson radii would sharpen the conclusion.","The constraint on the screening parameter could be cross-checked against independent determinations from quarkonium dissociation patterns or heavy-quark diffusion coefficients in the QGP, potentially tightening or challenging the bound.","The reversal of the radius-dependence ordering when the potential is included (smaller radius gives higher probability) implies that the potential effect dominates over the geometric overlap effect, which could be tested by comparing coalescence yields for mesons of different sizes."],"forward_implications":["The coalescence probability constraint bounds the in-medium modification of heavy-light potentials near the QCD phase transition, complementing lattice QCD and QCD sum-rule estimates.","Heavy-flavor observables in heavy-ion collisions, such as the baryon-to-meson ratio and elliptic flow at intermediate transverse momentum, can be described without ad hoc normalization of the coalescence probability.","The framework can be extended to study coalescence of heavy quarks with strange quarks by adjusting the light-quark mass, with predictions for Ds and Bs production.","The sensitivity of the coalescence probability to the screening parameter provides a phenomenological handle on extracting the in-medium heavy-quark potential from experimental heavy-flavor data."],"fun_headline_variants":["Realistic quark potential fixes heavy-quark coalescence deficit","Attractive potential drives heavy-quark coalescence to unity","Heavy-quark coalescence probability reaches unity with model potential","Potential-based model resolves low-momentum coalescence shortfall","Constituent quark potential raises coalescence probability to unity"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The claim that the coalescence probability approaches unity depends on the heavy-meson radius being at the smaller end of the assumed range (0.5 fm). If the actual meson radius is larger, the probability drops below unity even with the potential included.","fun_headline_variants_meta":{"raw":{"variants":["Realistic quark potential fixes heavy-quark coalescence deficit","Attractive potential drives heavy-quark coalescence to unity","Heavy-quark coalescence probability reaches unity with model potential","Potential-based model resolves low-momentum coalescence shortfall","Constituent quark potential raises coalescence probability to unity"]},"model":"glm-5.2","effort":"low","cost_usd":0.0,"raw_usage":{"total_tokens":537,"prompt_tokens":451,"completion_tokens":86,"prompt_tokens_details":null},"tokens_in":451,"tokens_out":86,"duration_ms":54921,"temperature":1.0,"reasoning_tokens":null,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-10T02:35:37.093065+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"If lattice QCD or experimental determinations of the D- or B-meson radius near the phase transition yield values significantly larger than 0.5 fm, or if the in-medium screening mass exceeds the threshold identified here, the coalescence probability would fall substantially below unity, undermining the claim that the potential alone resolves the deficit.","supporting_citations":[],"review_version":1}