{"id":"101a84e5-2d24-42d3-991b-b347f333f81d","arxiv_id":"2501.14872","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Sub-nucleonic hotspots in the initial state increase longitudinal flow decorrelation and reduce baryon stopping in simulated Pb+Pb collisions, but the model still underestimates decorrelation in mid-central events.","lead":"Using a hydrodynamics-coupled initial state model with sub-nucleonic hotspots, the authors compute rapidity-dependent flow and decorrelation in Pb+Pb collisions at 2.76 TeV. The work suggests that missing event-by-event fluctuations in the quark sector may explain why simulations underestimate longitudinal decorrelation in mid-central collisions.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Per-scenario K_g re-tuning shifts the quark/gluon energy partition, so the claimed effects of sub-nucleonic fluctuations are not isolated from a change in normalization.","rationale":"The reader's verdict correctly identifies the quark-fluctuation attribution as weakly supported, but the more foundational threat is the scenario-dependent K_g. The claimed new physics—that sub-nucleonic fluctuations alter longitudinal structure—rests on comparing five variants that differ simultaneously in the fluctuation prescription and in the gluon normalization. Since K_g enters only the gluon term in Eq. (1), a 43% increase from nucleon to hotspots+fluct sigma=1.2 (2.23 to 3.20) substantially reduces the quark-energy fraction at fixed midrapidity multiplicity. Baryon stopping is generated by the quark term (Eq. (2)), so the observed reduction in stopping may be a direct consequence of this re-weighting rather than of the spatial hotspot geometry. Likewise, the larger K_g amplifies the fluctuating gluon component relative to the smoother quark contribution, which could directly increase v2, v3 and r2. Without a control at fixed K_g, the paper's central attribution is underdetermined. The paper openly discloses the K_g values and the energy shifts, and the code is public, so this control is straightforward. If the trends survive fixed-K_g reruns, the paper's conclusions are strengthened; if not, the primary claim about sub-nucleonic fluctuation effects must be weakened to a statement about tuned-model variants. This is a concrete, settled-by-computation concern, and it does not require new physics or external data.","tokens_in":18124,"tokens_out":9231,"duration_ms":78763,"concrete_test":"Fix K_g to a common value, e.g., K_g=2.23 (the nucleon-case value) or 2.25 (the value used in Fig. 1), and rerun the five scenario sets of Sec. II D through the same CLVisc pipeline. Repeat the centrality selection and recompute the baryon-stopping ratios (Sec. III A), v2(eta), v3(eta), and r2(eta) (Secs. III C-III E). If the qualitative trends—reduced baryon stopping, enhanced flow, enhanced decorrelation with hotspots—persist at fixed K_g, the confound is not responsible. If they disappear or reverse, the central claim about sub-nucleonic fluctuations is not supported and the paper should be revised to separate normalization effects from fluctuation-structure effects.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central comparison across the five scenarios in Sec. II D is not at fixed initial-state physics: each scenario uses a separately tuned K_g (nucleon 2.23, nucleon+fluct 2.34, hotspots 2.61, hotspots+fluct sigma=0.637 2.71, hotspots+fluct sigma=1.2 3.20). Because K_g multiplies only the gluon term in Eq. (1), this re-tuning changes the relative weight of gluon vs. quark energy deposition. The paper itself notes in Sec. III A that 'the energy carried by gluons increases (due to the increase of Kg) and tends to concentrate in the mid-rapidity region while the energy carried by quarks decreases.' Consequently, the reported reduction of baryon stopping with hotspots may be caused by the decreased weight of quark-stopping energy, not by the geometric non-overlap of hotspots. Similarly, the enhanced v2(eta) and r2(eta) in hotspot runs could stem from the increased weight of the fluctuating gluon component, rather than from the sub-nucleonic spatial fluctuations per se. The comparisons are between tuned models with different physics content, so the causal attribution to sub-nucleonic fluctuations is not established. A control at fixed K_g is needed.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper couples the 3D resolved McDipper initial-state model, extended with sub-nucleonic fluctuations (hotspots and thickness fluctuations), to the 3+1D viscous hydrodynamics code CLVisc, and studies Pb+Pb collisions at sqrt(s_NN)=2.76 TeV. The authors compare five initial-state scenarios (smooth nucleons, hotspots, nucleon thickness fluctuations, and hotspot plus thickness fluctuations with two fluctuation strengths) and present results for charged-hadron and net-proton rapidity distributions, pT spectra, directed flow v1, elliptic and triangular flow v2/v3, and the longitudinal decorrelation observable r2 across centrality classes. Their central finding is that sub-nucleonic fluctuations reduce baryon stopping and enhance anisotropic flow and longitudinal decorrelation, especially in central collisions, while the model still underestimates decorrelation in mid-central and peripheral collisions, which they attribute to missing quark-sector fluctuations.","tokens_in":18417,"tokens_out":4393,"duration_ms":41052,"significance":"If the central claim holds, the paper provides a valuable step toward a more complete 3D description of heavy-ion collisions, showing that sub-nucleonic degrees of freedom leave an imprint on rapidity-dependent observables and should be included in initial-state models. The main strengths are the first coupling of the resolved McDipper initial state to 3+1D viscous hydrodynamics, the breadth of observables compared with ALICE and CMS data, and the public availability of the McDipper code (Ref. [55]). The paper is also honest about its limitations, explicitly noting the absence of a hadronic afterburner, the lack of a 3+1D pre-equilibrium stage, and the need for further tuning. However, the causal interpretation is weakened by the scenario-dependent re-tuning of K_g, which changes the quark/gluon energy partition and is therefore entangled with the geometric effects that the paper aims to isolate.","major_comments":[{"comment":"The central scenario comparisons are not at fixed physics content because K_g is re-tuned for each scenario (Sec. II D: K_g = 2.23, 2.34, 2.61, 2.71, 3.20). Since K_g multiplies only the gluon term in Eq. (1), increasing K_g does more than change the overall normalization: it increases the gluon energy share and decreases the quark energy share. The authors themselves note in Sec. III A that the energy carried by gluons increases due to the increase of K_g and tends to concentrate at mid-rapidity while the quark energy decreases. Consequently, the reported reduction of baryon stopping and the enhanced v2(eta), v3(eta), and r2 in the hotspot runs could be driven by the altered quark/gluon partition rather than by the geometric non-overlap of hotspots. A control calculation at fixed K_g, or a decomposition separating geometry effects from energy-sharing effects, is needed before the causal statements in Sec. IV can be supported.","section":"Sec. II D and Sec. III A"},{"comment":"The attribution of the remaining longitudinal-decorrelation deficit to missing quark-sector fluctuations is an interpretation, not a demonstrated result. The same deficit could receive significant contributions from the neglected 3+1D pre-equilibrium evolution (which the authors flag as desirable in Secs. III E and IV), from the omitted hadronic afterburner (whose absence they already connect to the too-hard pT spectra and overestimated flow in Secs. III B and III D), or from hydrodynamic thermal fluctuations, which they cite in Sec. III E as essential (Refs. [88, 89]). Since quark fluctuations are nowhere implemented in the model (Eq. (4) uses averaged collinear quark PDFs), the statement that the disagreement 'should be attributed to missing sources of fluctuations in the initial condition, most likely in the quark sector' goes beyond the evidence presented in this work. Please rephrase this as one candidate explanation among several, or perform a dedicated test that isolates the quark sector (e.g., sampling quark positions according to the local PDF density).","section":"Sec. III E and Sec. IV"},{"comment":"The model overestimates the magnitudes of v2 and v3 in several centrality classes, and the authors attribute this to the too-hard pT spectra. Because the longitudinal decorrelation r2 is defined through ratios of flow vectors, the quantitative agreement with the CMS r2 data is affected by this same model deficiency. The claim that the inclusion of sub-nucleonic fluctuations 'brings the shape of the rapidity dependence of elliptic flow closer to experimental data' is therefore only qualitative; the paper would benefit from a quantitative goodness-of-fit comparison, or from an explicit statement that the r2 comparison is illustrative rather than a quantitative validation of the missing-fluctuation hypothesis.","section":"Sec. III D and Figs. 9-12"}],"minor_comments":[{"comment":"In the abstract and Sec. II A, 'resolvedMcDipper' should be written as 'resolved McDipper' with a space.","section":"Abstract and Sec. II A"},{"comment":"The p+p multiplicity comparison in Fig. 2 uses K_g values tuned to Pb+Pb collisions and parameters in Eq. (12) from 5.02 TeV pp data, while the ALICE data shown are at 7 TeV; the authors acknowledge this, but the caption should state it explicitly to avoid misleading readers.","section":"Fig. 2 and Sec. II B"},{"comment":"The function a_J in Eq. (12) is not defined; please define it or provide a reference.","section":"Eq. (12)"},{"comment":"In the text below Fig. 9, 'the integrated McDipper+CLVisc model generally works well at low pT' is too strong given that the v2{2} and v2{4} curves visibly overshoot the data; consider using 'gives a reasonable description' or similar.","section":"Sec. III D"},{"comment":"The GitHub reference for McDipper is cited without a version tag or commit identifier; since the paper describes 'McDipper v1.2', please provide a versioned DOI or a specific release tag to ensure reproducibility.","section":"Ref. [55]"}],"recommendation":"major_revision","confidential_remarks":"The K_g re-tuning is the main technical obstacle to accepting the paper's causal claim. A control at fixed K_g or a clear re-framing of the conclusions as scenario comparisons under separate tuning would resolve the core issue. The quark-fluctuation attribution in Sec. III E should be softened, as the paper itself lists several other missing ingredients. The framework and the breadth of comparisons are solid and suitable for the journal; the revision should focus on the interpretation and on adding the suggested control or caveats."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is a genuinely useful piece of machinery. The authors couple the McDipper initial state, including sub-nucleonic hotspots and thickness fluctuations, to the CLVisc 3+1D viscous hydro code, and use it to study baryon stopping, v_n(eta), and r_2 decorrelation in 2.76 TeV Pb-Pb. The code is public, the methods are described in detail, and the paper is conspicuously honest about what is and is not included. That alone earns it a serious read.\n\nWhat it does well: the scan over fluctuation scenarios is systematic, and the qualitative comparisons to ALICE and CMS are reasonable. The observation that thickness fluctuations are needed to reproduce the pp multiplicity tail is a nice consistency check. The central numerical result—sub-nucleonic fluctuations increase longitudinal decorrelation at central rapidities—is interesting and worth knowing.\n\nThe soft spot is the one the stress-test flags, and it is real. The five scenarios are compared at different values of K_g, which multiplies only the gluon contribution in Eq. (1). The values run from 2.23 (nucleon) to 3.20 (hotspots+fluct sigma=1.2). The paper itself notes in Sec. III A that as K_g increases, gluon energy grows and quark energy drops. So the reported reductions in baryon stopping and the changes in v2(eta) and r2(eta) are not cleanly attributable to the spatial granularity of hotspots; they may partly reflect the altered quark/gluon partition. A control at fixed K_g, or at least a demonstration that the trends persist when the normalization is held fixed, is needed before the causal claim is solid. The speculation that missing quark-sector fluctuations explain the remaining underprediction of r_2 is plausible but untested; with no pre-equilibrium dynamics and no hadronic afterburner, the deficit has other candidate explanations. These limitations are stated openly, but they do weaken the headline interpretation.\n\nWho is this for: heavy-ion phenomenologists working on initial-state modeling and longitudinal dynamics. The paper deserves a serious referee. My recommendation: accept for review, but the referee should ask for a fixed-K_g control or a physical justification for the re-tuning, and for a softer statement of the quark-fluctuation conclusion. The machinery is valuable; the causal story needs another pass.","headline":"Solid new McDipper+CLVisc machinery with an honest but unproven central claim, confounded by per-scenario K_g re-tuning.","tokens_in":18965,"tokens_out":3029,"would_cite":true,"duration_ms":26516,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["25.75.-q","24.10.Nz"],"model":"deepseek-v4-flash","headline":"Sub-nucleonic fluctuations reshape longitudinal flow and baryon stopping in heavy-ion collisions.","keywords":["sub-nucleonic fluctuations","heavy-ion collisions","longitudinal decorrelation","baryon stopping","anisotropic flow","Color Glass Condensate","3+1D viscous hydrodynamics","Pb-Pb collisions at 2.76 TeV"],"falsifier":"Run the same McDipper+CLVisc setup with event-by-event fluctuating valence and sea quark distributions, as in the prescription of Ref. [87], and recompute r2 in 30-50% central Pb-Pb collisions; if the longitudinal decorrelation still falls short of the CMS data, the missing-fluctuations-in-the-quark-sector explanation is ruled out.","tokens_in":17884,"feed_emoji":"⚛️","tokens_out":4630,"duration_ms":36481,"temperature":0.7,"pith_summary":"This paper asks whether the lumpy, sub-nucleonic structure of colliding nuclei changes how energy and baryon number are deposited along the collision axis, and whether that structure leaves a measurable imprint on final flow patterns. It builds a complete simulation chain in which the Color Glass Condensate based McDipper initial state, extended with hotspot and thickness fluctuations inside each nucleon, feeds the 3+1 dimensional viscous hydrodynamic code CLVisc. The authors find that sub-nucleonic fluctuations reduce baryon stopping and enhance anisotropic flow and longitudinal decorrelation, with the strongest effects in central collisions. Comparing against CMS and ALICE data, the model describes central-collision decorrelation but underpredicts it in mid-central and peripheral collisions, which the authors interpret as evidence that quark-sector fluctuations are still missing.","feed_headline":"Sub-nucleon hotspots boost flow and cut baryon stopping","feed_subtitle":"First 3D McDipper+CLVisc run links initial-state granularity to rapidity-dependent flow and decorrelation data.","key_machinery":"The working machinery is the 3D resolved McDipper initial-state model, a kT-factorized Color Glass Condensate formulation in which gluon production comes from dipole amplitudes and quark production from valence-quark stopping, extended with sub-nucleonic fluctuations: each nucleon is a weighted superposition of Gaussian hotspots whose weights carry log-normal thickness fluctuations. These energy and net-baryon densities are matched onto the 3+1D viscous hydrodynamics code CLVisc, which evolves the fireball including baryon diffusion, and the longitudinal decorrelation observable rn is computed from the forward/backward ratio of flow vectors. The hotspot granularity is the ingredient that breaks the fireball into steeper density gradients, while thickness fluctuations add event-by-event weight variations; together they shape the rapidity dependence of vn and the decorrelation rate.","core_discovery":"The central claim, stated in Sec. IV, is that sub-nucleonic fluctuations in the initial state measurably change the longitudinal structure of heavy-ion collisions: including three hotspots and log-normal thickness fluctuations per nucleon suppresses baryon stopping, increases anisotropic flow, and strengthens longitudinal decorrelation, especially in central collisions. The paper further claims that the residual underprediction of the longitudinal decorrelation ratio r2 in mid-central and peripheral collisions points to a missing source of fluctuations in the quark sector, because at forward rapidity a growing share of the deposited energy comes from quark stopping and the current implementation uses only averaged collinear quark parton distributions.","pith_inferences":["If quark-sector fluctuations do close the r2 gap, the same mechanism should also affect net-proton cumulants and charge-balance observables at forward rapidity, giving independent tests beyond flow decorrelations.","The hotspot prescription could be tested by varying the number of hotspots Nq and their width Bq separately to see whether r2 constrains them individually or only in combination.","A quantitative match to the CMS data will likely require adding pre-equilibrium dynamics and a hadronic afterburner, so the attribution of the deficit to quark fluctuations should be revisited once those stages are included."],"forward_implications":["Sub-nucleonic fluctuations must be included in initial-state models if rapidity-dependent observables such as v2(η), v3(η), and r2 are to be described quantitatively.","Baryon stopping is sensitive to the granularity of the initial state, so baryon-number transport measurements can constrain the number and width of hotspots.","Because sub-nucleonic fluctuations mostly enhance flow in central collisions and suppress it in peripheral ones, the centrality dependence of vn serves as a clean diagnostic of initial-state granularity.","The underprediction of r2 in mid-central and peripheral collisions motivates a statistical treatment of quark parton distributions in saturation-based initial conditions.","Longitudinal decorrelation at large rapidity receives a growing contribution from quark stopping, so forward measurements and mid-rapidity references are needed to disentangle gluon and quark deposition mechanisms."],"supporting_citations":[{"why":"Provides the 3D resolved McDipper initial-state model that is the foundation of the simulation chain.","marker":"[34]"},{"why":"Introduces the hotspot and log-normal thickness fluctuation prescription used to add sub-nucleonic structure.","marker":"[35, 36]"},{"why":"Supplies the CLVisc (3+1)D viscous hydrodynamics code that evolves the initial state to final hadrons.","marker":"[7, 8]"},{"why":"Provides the CMS longitudinal decorrelation data r2 that the paper compares against and underpredicts in mid-central collisions.","marker":"[18]"},{"why":"Provides the ALICE rapidity-dependent anisotropic flow data used to assess v2(η) and v3(η).","marker":"[17]"},{"why":"Fixes the nucleon width and hotspot geometry conventions taken from the TRENTo model.","marker":"[56]"},{"why":"Is the statistical quark PDF prescription cited as the likely missing source of quark-sector fluctuations.","marker":"[87]"},{"why":"Documents the growing relative contribution of quark stopping at forward rapidity that motivates the quark-fluctuation deficit.","marker":"[86]"}],"fun_headline_variants":["Sub-nucleon hotspots boost flow and cut stopping","Hotspots inside nucleons alter heavy-ion longitudinal flow","Rapidity decorrelation tied to sub-nucleon hotspots","Sub-nucleon fluctuations change flow and stopping","3D model links sub-nucleon granules to flow decorrelation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's interpretation rests on the assumption that the longitudinal structure of the initial state is faithfully described by the kT-factorized CGC formulas with smooth, averaged collinear quark parton distributions, and that the omitted pre-equilibrium evolution and hadronic afterburner do not account for the remaining decorrelation deficit.","fun_headline_variants_meta":{"raw":{"variants":["Sub-nucleon hotspots boost flow and cut stopping","Hotspots inside nucleons alter heavy-ion longitudinal flow","Rapidity decorrelation tied to sub-nucleon hotspots","Sub-nucleon fluctuations change flow and stopping","3D model links sub-nucleon granules to flow decorrelation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000837,"raw_usage":{"total_tokens":3573,"prompt_tokens":792,"completion_tokens":2781,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":408,"completion_tokens_details":{"reasoning_tokens":2701}},"tokens_in":408,"tokens_out":2781,"duration_ms":18775,"temperature":1.0,"reasoning_tokens":2701,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T14:50:42.946872+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same McDipper+CLVisc setup with event-by-event fluctuating valence and sea quark distributions, as in the prescription of Ref. [87], and recompute r2 in 30-50% central Pb-Pb collisions; if the longitudinal decorrelation still falls short of the CMS data, the missing-fluctuations-in-the-quark-sector explanation is ruled out.","supporting_citations":[{"cited_title":"Pseudorapidity dependence of the anisotropic flow of charged particles in Pb-Pb collisions at $\\sqrt{s_{\\rm NN}}=2.76$ TeV","cited_arxiv_id":"1605.02035","evidence_quote":"Provides the ALICE rapidity-dependent anisotropic flow data used to assess v2(η) and v3(η)."}],"review_version":1}