{"id":"d88b0613-b167-421d-8d47-9a5904ae4f1a","arxiv_id":"2607.11017","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":5.5,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"Prompt D0 v2 is strongly linearly correlated with low-pT charged-particle v2 across q2-selected event classes, indicating initial-state geometry dominates charm-hadron elliptic flow in PbPb collisions.","lead":"CMS measures a strong linear correlation between D0 meson elliptic flow and charged-particle flow when events are selected by similar initial geometry via event-shape engineering. This shows that the initial collision shape largely drives charm-hadron collective motion in the quark-gluon plasma.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified","rationale":"The Reader correctly isolates the soft-particle proxy as the softest assumption and correctly judges that it does not invalidate the data-driven correlation. The manuscript supplies the necessary supporting measurements (q2–v2 linearity, large η-gap, normalized slopes/intercepts, systematic table) and the claim is phrased as a suggestion of substantial impact rather than a claim of exclusive hydrodynamic origin. High-pT deviations and the incomplete PHSD thermalization are acknowledged in the text and do not undercut the low-pT correlation that carries the central message. No adjustment to ACCEPT is warranted.","tokens_in":40859,"tokens_out":508,"duration_ms":4668,"concrete_test":"Recompute the Pearson r values of Fig. 4 after replacing the 1<pT<3 GeV charged-particle proxy with an independent mid-rapidity eccentricity estimator (e.g., two-particle cumulant v2{2} with a larger η-gap or participant-plane eccentricity from a Glauber+hydro hybrid on the same events); if r remains consistent with unity for 2<pT<10 GeV in 10–40% centrality, the proxy assumption is confirmed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that D0 v2 is strongly linearly correlated with low-pT charged-particle v2 when both are selected in common q2 intervals, implying that initial-state geometry substantially drives charm-hadron elliptic flow. The weakest link identified by the Reader (charged-particle v2 as eccentricity proxy) is already mitigated by the large rapidity gap (HF 3<|η|<5 vs. mid-rapidity POI), the observed linearity of charged-particle v2 vs. q2 itself (Fig. 2), the Pearson coefficients near unity for pT<10 GeV outside the most central bin (Fig. 4), and the normalized slopes consistent with 1 with intercepts consistent with 0 (Fig. 5). Residual non-flow or viscous non-linearity would have to conspire to preserve these relations across five centrality classes and four pT bins while still producing the reported systematics; that is not a load-bearing failure of the argument. The PHSD comparison underestimates thermalization but is used only for illustration, not as proof. No internal inconsistency or untested assumption that would reverse the claim is present.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript reports a CMS measurement of prompt D0 meson elliptic flow (v2) in PbPb collisions at √sNN = 5.02 TeV (0.607 nb−1) using event-shape engineering. Events are classified by a forward calorimeter asymmetry parameter q2 (HF, 3 < |η| < 5) within 1% centrality slices that are then recombined into 10% classes, suppressing centrality bias. Scalar-product v2 is extracted for D0 (2 < pT < 30 GeV, |y| < 1) and for low-pT charged particles (1 < pT < 3 GeV, |η| < 1) in common q2 intervals. The D0 v2 is found to be strongly linearly correlated with the charged-particle v2 (Pearson r near unity for pT < 10 GeV outside the most central bin; normalized slopes consistent with 1 and intercepts with 0). The authors interpret this as evidence that initial-state geometry substantially drives charm-hadron elliptic flow. PHSD model comparisons are shown and indicate underestimation of charm thermalization at low pT.","tokens_in":41102,"tokens_out":967,"duration_ms":7650,"significance":"The result supplies a clean, data-driven demonstration that charm-hadron v2 tracks the same event-by-event eccentricity variations that govern bulk flow, extending earlier ALICE ESE work with finer q2 binning, a larger pT reach, and explicit Pearson/slope/intercept quantification. The analysis is technically solid: scalar-product method with large rapidity gap, BDT selection, DCA-based 95% prompt purity, acceptance-efficiency corrections, and a documented systematic table. The normalized slopes ≈ 1 and intercepts ≈ 0 across five centralities and four pT bins constitute a falsifiable, nearly parameter-free statement of linear response. Residual non-flow or viscous non-linearity would have to conspire across many bins while remaining within the quoted systematics; that is not a load-bearing failure. The paper therefore strengthens the case that heavy quarks participate in the collective expansion of the QGP and provides a useful benchmark for transport models.","major_comments":[],"minor_comments":[{"comment":"Section 5.1: the modified scalar-product procedure that bins candidates in p(D0) and extracts a yield-weighted average is clear, but a short explicit statement that the same procedure is applied independently inside each q2 class would remove any residual ambiguity about order of operations.","section":null},{"comment":"Figure 3 and Appendix A: the linear fits are performed with statistical uncertainties only; a brief note that systematic uncertainties on the D0 points are shown as boxes but not included in the fit would help the reader assess the quoted slope/intercept errors.","section":null},{"comment":"Section 7 / Figure 4: the statement that r is consistent with unity for pT < 10 GeV outside 0–10% is supported by the data, but the larger uncertainties at 10–30 GeV make the claimed weakening of correlation only suggestive; a single clarifying sentence would avoid over-interpretation.","section":null},{"comment":"Table 1: the ranges of absolute differences are useful, but listing the dominant source for the slope and intercept in each centrality/pT bin (or at least noting that centrality and nonprompt dominate the slope) would improve transparency.","section":null},{"comment":"Typographical: “Crytsal Ball” appears twice in Section 4; correct to “Crystal Ball”.","section":null},{"comment":"References: the recent ALICE Λc and D-meson flow results (arXiv:2603.18966) and the CMS Ds paper (arXiv:2602.14221) are already cited; ensuring the final published versions are updated before production is advisable.","section":null}],"recommendation":"accept","confidential_remarks":"The central claim is robust and the technical execution meets the standard of the journal. The weakest link (charged-particle v2 as eccentricity proxy) is already mitigated by the large rapidity gap, the observed linearity of charged-particle v2 vs q2, and the near-unity Pearson coefficients; I do not regard it as a reason for revision. The paper is suitable for acceptance essentially as is."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"This is a clean experimental paper that does what it claims. CMS takes the ALICE event-shape engineering idea for D-meson flow and runs it with higher statistics, ten q2 classes per 10% centrality, four pT bins out to 30 GeV, Pearson coefficients, and normalized slopes/intercepts. The new quantitative pieces are the near-unity r values for pT < 10 GeV (outside the most central bin), the slopes consistent with 1 after q2-inclusive normalization, and the intercepts consistent with zero. That package is useful for anyone tuning heavy-quark transport or checking how much of the charm v2 is geometry-driven.\n\nWhat they do well is straightforward: large rapidity gap (HF 3 < |η| < 5 vs mid-rapidity POI), scalar-product method, BDT + DCA prompt purity at 95 %, acceptance-efficiency corrections, and a transparent systematics table. Fig. 2 already shows the soft-particle v2 tracking q2 linearly, so the proxy is not pure assumption. The PHSD comparison is honest about underestimating thermalization at low pT; they treat it as illustration, not proof.\n\nSoft spots are real but proportional. The charged-particle v2 (1–3 GeV) is still a proxy for eccentricity; residual non-flow or viscous non-linearity could in principle matter, especially in 0–10% and at pT > 10 GeV where the correlation visibly softens (Appendix figures). That is expected physics, not a hidden flaw, and the paper flags it. No circularity: the correlation is measured, not imposed. Citation pattern is normal for a CMS heavy-ion note.\n\nThis is for people who care about charm thermalization and initial-geometry response in the QGP. It will not rewrite the field, but it is a solid data point that deserves a serious referee. I would accept it for peer review without hesitation and would cite the slopes and Pearson numbers when I next need a quantitative handle on geometry-driven charm flow.","headline":"Solid CMS extension of ALICE event-shape engineering for D0 v2; data-driven linear correlation with soft-particle proxy is clean and well-documented, with only the expected high-pT/central softening.","tokens_in":41684,"tokens_out":517,"would_cite":true,"duration_ms":6320,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["25.75.Ld","25.75.Cj","12.38.Mh"],"model":"grok-4.5","headline":"D0 meson elliptic flow tracks charged-particle flow when events are sorted by the same initial geometry, showing that charm quarks inherit the collision shape.","keywords":["elliptic flow","D0 mesons","event-shape engineering","quark-gluon plasma","heavy-ion collisions","charm thermalization","initial-state geometry","PbPb 5.02 TeV"],"falsifier":"A high-statistics measurement in the same q2 classes that finds the D0–charged-particle v2 correlation coefficient dropping significantly below unity (or becoming non-linear) already at low pT (2–4 GeV) in mid-central collisions would undermine the claim that initial geometry is the dominant driver.","tokens_in":41758,"feed_emoji":"⚛️","tokens_out":1021,"duration_ms":11587,"temperature":0.7,"pith_summary":"In lead-lead collisions at 5.02 TeV, the CMS experiment uses event-shape engineering to test whether the elliptic flow of prompt D0 mesons is driven by the initial geometry of the quark-gluon plasma. An asymmetry parameter q2 is measured from forward transverse-energy deposits and used to bin events of similar centrality into classes that share similar elliptic eccentricity. Within those classes the D0 meson v2 (mid-rapidity, 2–30 GeV) is found to rise and fall in lock-step with the low-pT charged-particle v2 that serves as a proxy for the eccentricity. The resulting linear correlation, quantified by Pearson coefficients near unity for most centralities and pT below 10 GeV, implies that the initial spatial asymmetry is the dominant source of charm-hadron collective motion. A transport-model comparison underestimates the low-pT slope, suggesting incomplete thermalization of charm quarks in that calculation.","feed_headline":"Charm mesons inherit the shape of the quark-gluon plasma","feed_subtitle":"D0 elliptic flow tracks charged-particle flow when events share the same initial geometry","key_machinery":"Event-shape engineering with the reduced flow vector q2, constructed from transverse-energy deposits in the hadron-forward calorimeters (3 < |η| < 5). Events of fixed centrality are further sorted into ten q2 percentiles so that the mid-rapidity D0 and charged-particle v2 can be compared at matched initial eccentricity.","core_discovery":"When PbPb events are subdivided by the forward asymmetry parameter q2, the elliptic-flow coefficient v2 of prompt D0 mesons is strongly and linearly correlated with the v2 of low-pT charged particles measured in the same q2 intervals. The correlation holds across centrality classes and for D0 transverse momenta from 2 to 30 GeV, indicating that the initial-state geometry substantially controls the development of charm-hadron flow.","pith_inferences":["Because the correlation remains linear even after the D0 and charged-particle v2 values are normalized to their inclusive averages, the relative response of charm to geometry appears universal across the measured centrality and pT range.","The weaker correlation observed in the most central (0–10 %) and most peripheral (40–50 %) bins may already encode residual viscosity or non-flow contributions that become visible once the average eccentricity is small.","If the same linear relation is confirmed for bottom-flavor hadrons, the mass hierarchy of heavy-quark thermalization times would be tightly constrained by a single geometric observable."],"forward_implications":["Charm quarks participate in the collective expansion of the QGP sufficiently early and strongly that their final-state anisotropy mirrors the initial geometric eccentricity.","Transport models that under-predict the low-pT slope of the D0–charged-particle correlation must increase the degree of charm thermalization to match the data.","At high pT (10–30 GeV) any residual weakening of the correlation can be used as a quantitative handle on path-length-dependent energy loss of charm quarks.","The same q2-sorting technique can be applied to other heavy-flavor species (Ds, Λc) to test whether strangeness or baryon number alters the geometry-driven response."],"fun_headline_variants":["D0 elliptic flow tracks charged particles in shared geometry events","Initial geometry drives charm meson v2 in PbPb collisions","q2 selection links D0 flow tightly to charged-particle anisotropy","Event-shape engineering reveals geometry controls of D0 flow","Charm hadrons inherit elliptic flow from initial-state eccentricity"],"cache_read_input_tokens":32896,"weakest_assumption_plain":"The paper treats the measured low-pT charged-particle v2 as a faithful experimental stand-in for the true event-by-event initial eccentricity, relying on an approximately linear hydrodynamic response and a large rapidity gap that suppresses non-flow effects.","fun_headline_variants_meta":{"raw":{"variants":["D0 elliptic flow tracks charged particles in shared geometry events","Initial geometry drives charm meson v2 in PbPb collisions","q2 selection links D0 flow tightly to charged-particle anisotropy","Event-shape engineering reveals geometry controls of D0 flow","Charm hadrons inherit elliptic flow from initial-state eccentricity"]},"model":"grok-4.5","effort":"low","cost_usd":0.004218,"raw_usage":{"total_tokens":1331,"prompt_tokens":894,"num_sources_used":0,"completion_tokens":85,"cost_in_usd_ticks":42180000,"prompt_tokens_details":{"text_tokens":894,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":352,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":894,"tokens_out":85,"duration_ms":3734,"temperature":1.0,"reasoning_tokens":352,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-14T07:36:44.610725+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"A high-statistics measurement in the same q2 classes that finds the D0–charged-particle v2 correlation coefficient dropping significantly below unity (or becoming non-linear) already at low pT (2–4 GeV) in mid-central collisions would undermine the claim that initial geometry is the dominant driver.","supporting_citations":[],"review_version":1}