{"id":"94244d23-f4d5-4d39-8ec9-05a42ac60d13","arxiv_id":"2504.19208","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Imposing a minimum nucleon separation or sub-nucleon fluctuations in TRENTo initial conditions reduces initial eccentricities and lets CLVisc hydrodynamics with lower shear viscosity come closer to the measured v2/v3 ratio in ultra-central Pb+Pb collisions.","lead":"Hydrodynamic simulations still struggle to reproduce the near-equal elliptic and triangular flow seen in ultra-central lead-lead collisions. This paper tests two tweaks to the initial nuclear shape, a minimum nucleon spacing and sub-nucleon hotspots, and shows both reduce the discrepancy.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The dmin sampling is not shown to preserve the one-body Pb density, and eta/s is re-tuned per scenario; the nuclear-structure attribution of the v2/v3 improvement is therefore unproven.","rationale":"Good-faith read: the paper is a careful proof-of-principle hydrodynamic study. The TRENTo+CLVisc framework is standard, the calibration to pT spectra and centrality-dependent yields is reasonable, and the authors are candid in Section II.B that dmin does not imply a real lower bound on nucleon separation. However, the strongest claim—that initial-state nuclear structure is a critical factor in the ultra-central puzzle—rests on a knob whose physical content is not verified. The rejection method that enforces dmin does not preserve the one-body density, so the configurations may no longer represent 208Pb even at the level of the measured charge density. This directly impacts the eccentricities in Fig. 6 and the flow coefficients in Figs. 8-10; a flattened density profile would reduce epsilon_n regardless of any correlation physics. The per-scenario eta/s choices (0.16, 0.18, 0.22) complicate the interpretation further, since the reported v2/v3 improvement could be partly or wholly due to viscosity tuning rather than initial-state structure. The proposed tests are inexpensive: computing the one-body density from the existing sampled configurations is a post-processing step, and rerunning a few CLVisc events at fixed eta/s is small relative to the 1000-event ensembles already produced. If the one-body density is preserved and the scenario ranking persists at fixed eta/s, the conditional verdict can be upgraded; otherwise, the nuclear-structure interpretation is not justified. These checks are appropriate conditions for accepting the paper's central claim.","tokens_in":12311,"tokens_out":12825,"duration_ms":130797,"concrete_test":"Recompute the one-body density rho(r) and the rms radius of the sampled nucleon configurations for dmin = 0, 1.0, 1.4, and 1.7 fm using the same rejection algorithm described in Section II.B, and compare with the unmodified Woods-Saxon density of Eq. (2) (R = 6.62 fm, a = 0.546 fm) and with the measured 208Pb charge radius. If the central density or rms radius changes by more than a few percent, the dmin setup does not respect the constrained one-body distribution, and the eccentricity reduction in Fig. 6 cannot be attributed to two-body correlations. As a complementary check, run the 0-1% CLVisc comparisons at a common eta/s = 0.16 and at eta/s = 0.22 for all dmin and sub-nucleon cases: if the ranking of the scenarios by agreement with the CMS v2{2} and v3{2} data reverses with eta/s, then the viscosity choice, not the initial-state structure, is what drives the reported improvement.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central claim that nuclear and sub-nucleon structures are critical for the ultra-central v2/v3 puzzle rests on the minimum-separation knob in Section II.B. The sequential rejection sampling discards candidates closer than dmin to previously accepted nucleons; since the Woods-Saxon density of Eq. (2) peaks at the center, this preferentially rejects central candidates and flattens the radial profile. The authors state that the single-nucleon density and charge radius are experimentally constrained, but they never calculate, reweight, or otherwise verify that the dmin-modified configurations still match the measured 208Pb density or radius. If the dmin = 1.4 fm configurations are flatter and more extended than the real nucleus, then the reduced eccentricities in Fig. 6 and the improved v2/v3 agreement in Fig. 8 are an artifact of substituting a different mean density profile rather than evidence for realistic short-range correlations. The paper itself acknowledges that the dmin knob is a proof-of-principle toy, which is honest, but the abstract and discussion overreach by claiming implications for nuclear structure. A secondary confound is that the shear viscosity is re-tuned per scenario (eta/s = 0.22 in Figs. 8-9, 0.18 in Fig. 10, versus 0.16 in the calibration), so the improvement is not cleanly isolated from the viscosity choice.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the ultra-central v2/v3 puzzle in Pb+Pb collisions at 2.76 TeV using TRENTo initial conditions coupled to the CLVisc hydrodynamic model. Two modifications to the initial state are considered: imposing a minimum separation dmin between nucleons in the Woods-Saxon sampling, and adding sub-nucleon constituent-quark fluctuations. The authors report that increasing dmin reduces the root-mean-square initial eccentricities, and that combined with a re-tuned shear viscosity this brings vn{2} closer to CMS data in the 0-1% centrality bin. Sub-nucleon fluctuations are also shown to reduce eccentricities and to improve the v2/v3 trend. The manuscript concludes that initial-state nuclear and sub-nucleon structures are critical factors in the ultra-central puzzle.","tokens_in":12647,"tokens_out":4862,"duration_ms":51015,"significance":"If the central conclusion were fully established, the paper would be a valuable demonstration that short-range nuclear correlations and sub-nucleon granularity can affect ultra-central flow observables, with implications for both nuclear-structure studies and the extraction of QGP transport coefficients. The authors are transparent about the proof-of-principle nature of the minimum-distance prescription, and they compare against ALICE, ATLAS, and CMS data using an established hydrodynamic framework. The strength of the conclusion, however, is conditional on verifying that the dmin sampling preserves the known 208Pb density and on isolating the effect of the initial-state modification from the per-scenario re-tuning of eta/s. As written, the evidence supports a weaker statement: certain correlated initial geometries can improve the ultra-central vn description within this model setup.","major_comments":[{"comment":"The rejection sampling used to impose dmin is not shown to preserve the single-nucleon density of Eq. (2). Because the Woods-Saxon density peaks at r=0, sequentially discarding nucleons closer than dmin to previously accepted nucleons preferentially removes central candidates and flattens the radial profile. The statement in Section II.B that the single-nucleon density and charge radius are experimentally constrained is an assumption, not a verification; the paper never computes or reweights the one-body density or radius of the dmin=1.4 configurations. Without this check, the reduced eccentricities in Fig. 6 and the improved vn{2} in Fig. 8 could be caused by a modified mean radial profile rather than by short-range correlations.","section":"Section II.B, Fig. 1"},{"comment":"The shear viscosity is re-tuned for each initial-state scenario: eta/s=0.16 is used for the calibration and non-central comparisons in Figs. 3-5, eta/s=0.22 for the dmin scans in Figs. 8 and 9, and eta/s=0.18 for the sub-nucleon comparison in Fig. 10. Since both the initial geometry and the viscosity are changed simultaneously, the improvement in the v2/v3 ratio cannot be cleanly attributed to the initial-state structure alone. The abstract's claim that the initial-state modifications 'reduce required viscosity' is also not directly demonstrated by these plots, because the dmin=1.4 runs use a larger eta/s than the baseline calibration.","section":"Section III.C, Figs. 8-10"},{"comment":"The exclusion of sub-nucleon fluctuations for the dmin=1.4 case is justified post hoc: the text states that including them would make v3{2} and v4{2} significantly underestimate the data. Selecting the preferred configuration after inspecting the outcome removes the predictive content of the claim that sub-nucleon structure is a critical factor. The paper should show the dmin=1.4 plus sub-nucleon result and discuss the tension explicitly rather than omitting it from the main comparison.","section":"Section III.C, Fig. 10"},{"comment":"The abstract's conclusion that initial-state nuclear and sub-nucleon structures are 'critical factors' goes beyond the paper's own caveat in Section II.B that the minimum-separation prescription is a proof-of-principle study. The discussion should be reframed as a demonstration that a specific class of correlated initial geometries can improve the ultra-central v2/v3 description in this model setup, pending validation against measured nuclear densities and a consistent viscosity treatment.","section":"Abstract and Section IV"}],"minor_comments":[{"comment":"Please clarify how the centrality selection via the 'mult' parameter is affected by dmin; the statement in Section III.A that dmin does not influence centrality dependence is supported only for dmin=1.4 by Fig. 3, not for the full range of dmin values used in Fig. 6.","section":"Section III.B, Fig. 6"},{"comment":"The caption is confusing: it says the results consider 'nuclear-nuclear correlation (MUSIC + IP-glasma)' and 'without nuclear-nuclear correlation,' but the legend identifies blue triangles as MUSIC/IP-glasma results and red symbols as CLVisc results; please separate the model labels from the correlation labels.","section":"Fig. 8 caption"},{"comment":"The Gaussian distribution for constituent quarks appears to be missing parentheses around the coordinate differences; the expression '(x-x'2)+(y-y'2)' should be a squared Euclidean distance.","section":"Eq. (6)"},{"comment":"There are several typographical issues, including 'thev2-v3 gap' in the abstract, 'dij > 1/...' in the Fig. 1 caption, and reference 41 with 'Nature Communications' embedded in the title; a careful proofread is needed.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is best viewed as a proof-of-principle study rather than a definitive nuclear-structure claim. The main technical concerns, density preservation and viscosity re-tuning, are testable within the current framework, so a major revision with additional validation and a consistent comparison protocol is appropriate. If the authors instead prefer to keep the current framing, the claims in the abstract and discussion should be substantially weakened."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The one thing to know: this is a useful sensitivity study, not a demonstration that nuclear structure resolves the ultra-central v2/v3 puzzle. The authors scan a minimum nucleon separation dmin and add subnucleon fluctuations in TRENTo, feed into CLVisc, and show that both modifications lower the initial eccentricities and narrow the v2/v3 gap toward CMS data. That is a legitimate and well-executed numerical exercise, with proper calibration to pseudorapidity, pT spectra, and differential flow before the ultra-central comparisons.\n\nThe real soft spot is the dmin knob. The rejection sampling they use preferentially discards central nucleons because the Woods-Saxon density is highest there, so the sampled configurations are flatter and more extended than the physical 208Pb profile. They assert that single-nucleon density and charge radius are constrained, but they never compute the sampled rms radius or one-body density for dmin = 1.4 fm. If the profile has changed, the reduced eccentricity in Fig. 6 and the improved v2/v3 in Fig. 8 are partly an artifact of substituting a different mean shape, not evidence for short-range correlations. The paper's own disclaimer that dmin is proof-of-principle is honest, but the abstract and discussion overclaim by calling nuclear structure a 'critical factor.'\n\nSecond, the viscosity is re-tuned per scenario: eta/s = 0.16 for calibration, 0.22 for the dmin = 1.4 results, and 0.18 for the subnucleon case. That means the improvement is not isolated to the initial state; the authors are effectively scanning both dmin and eta/s simultaneously. Third, the exclusion of subnucleon fluctuations for dmin = 1.4 is post hoc — they say it would underestimate v3 and v4, which is cherry-picking rather than a systematic comparison.\n\nThat said, the paper is internally consistent, uses a well-known hydro code, and provides a clean demonstration that more uniform initial conditions lower epsilon_n and allow smaller viscosity. That is a useful constraint for the community, especially for Bayesian extractions of transport coefficients that assume a standard Woods-Saxon. The stress-test concern about the one-body density holds up and should be the central referee request.\n\nI would send this to peer review, because it addresses a real puzzle and the numerics deserve scrutiny. But the authors need to verify the sampled density profile, run a joint dmin–eta/s scan, and tone down the nuclear-structure claim. It is not yet a paper I would cite for a physical conclusion, only as a sensitivity benchmark.","headline":"A careful, well-calibrated simulation study that shows initial-state uniformity can reduce the ultra-central v2/v3 gap, but the dmin knob is unvalidated and the viscosity is re-tuned per scenario, so the nuclear-structure claim overreaches.","tokens_in":13168,"tokens_out":1891,"would_cite":false,"duration_ms":20137,"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":"The ultra-central $v_2\\{2\\}/v_3\\{2\\}$ puzzle in heavy-ion collisions is substantially an initial-state nuclear-structure effect.","keywords":["ultra-central heavy-ion collisions","anisotropic flow","v2/v3 puzzle","nucleon-nucleon correlations","sub-nucleon fluctuations","initial eccentricities","shear viscosity","TRENTo"],"falsifier":"Measure the two-nucleon separation distribution in $^{208}$Pb (from ab initio nuclear-structure calculations or from electron-scattering and short-range-correlation experiments). If it shows a substantial population of pairs closer than 1.4 fm, rerun the TRENTo+CLVisc calculation with that realistic distribution; loss of the $v_2\\{2\\}/v_3\\{2\\}$ agreement would falsify the claim that such uniformity is the physical cause.","tokens_in":12106,"feed_emoji":"🌀","tokens_out":6210,"duration_ms":55503,"temperature":0.7,"pith_summary":"This paper argues that the long-standing mismatch in ultra-central Pb+Pb collisions, where hydrodynamic models underpredict $v_3$ relative to $v_2$, is not just a transport-coefficient problem: it depends on how nucleons are distributed inside the colliding nuclei. Using TRENTo initial conditions and the CLVisc hydrodynamic code, the authors show that imposing a minimum separation distance between nucleons and adding sub-nucleon fluctuations both reduce the fluctuation-driven eccentricities of the initial fireball. With smaller eccentricities, a smaller shear viscosity can be used, which preferentially boosts the higher-order harmonic $v_3$ and narrows the gap between $v_2\\{2\\}$ and $v_3\\{2\\}$ toward the CMS data. If this is right, ultra-central flow becomes a sensitive probe of nuclear and sub-nucleon structure, and transport-coefficient extraction must account for initial-state correlations.","feed_headline":"Uniform nuclei close ultra-central flow puzzle gap","feed_subtitle":"Simulations with spaced nucleons and sub-nucleon hot spots cut required viscosity and match CMS data.","key_machinery":"The engine of the argument is a geometric one: the rejection-sampling constraint in TRENTo that enforces a minimum distance $d_{\\min}$ between nucleon centers, together with optional $N_c = 3$ Gaussian sub-nucleon constituents (\\'hot spots\\'). The former homogenizes the sampled nucleon distribution and suppresses the fluctuation-driven eccentricities $\\epsilon_n\\{2\\}$; the latter sharpens the entropy-density spots and further changes the eccentricities. These reduced eccentricities change the hydrodynamic response so that a lower shear-viscosity-to-entropy ratio $\\eta/s$ can reproduce the experimental $v_2\\{2\\}$ and $v_3\\{2\\}$, reversing the usual pattern where a large $\\eta/s$ needed to fit $v_2$ kills $v_3$.","core_discovery":"The paper's central claim is that the ultra-central $v_2\\{2\\}/v_3\\{2\\}$ ratio puzzle is substantially an initial-state nuclear-structure effect. The authors modify the TRENTo sampling of $^{208}$Pb nuclei by enforcing a minimum nucleon separation $d_{\\min} = 1.0$–$1.7$ fm to homogenize the nucleon distribution, and separately amplify sub-nucleon fluctuations by sampling $N_c = 3$ Gaussian constituent quarks inside each nucleon. Both modifications lower the event-averaged eccentricities $\\epsilon_n\\{2\\}$ for all orders $n = 2$–$6$. This lowering matters because viscous damping of higher harmonics is controlled by $\\eta/s$: with less eccentricity, the hydrodynamic run can use a smaller $\\eta/s$ without over-suppressing $v_3$, so the calculated $v_2\\{2\\}$ and $v_3\\{2\\}$ move closer to the CMS measurements in 0–1% centrality. The best agreement is obtained with $d_{\\min} = 1.4$ fm (or $d_{\\min} = 1.0$ fm plus sub-nucleon fluctuations) and a reduced shear viscosity, while $dN_{\\rm ch}/d\\eta$ and $p_T$ spectra remain consistent with ALICE data across centralities.","pith_inferences":["If the mechanism holds, the same correlated-nucleon treatment should affect ultra-central flow in other species, such as $^{129}$Xe or $^{238}$U, where nuclear deformation and correlations differ; comparing those data would test the universality of the $d_{\\min}$ knob.","A natural next step is to replace $d_{\\min}$ with a microscopically derived short-range repulsive core from nuclear forces; if a realistic nucleon-nucleon potential reproduces the $d_{\\min} \\approx 1.4$ fm effect, the proof-of-principle becomes a quantitative nuclear-structure statement.","Bayesian parameter estimations of $\\eta/s$ that use Woods-Saxon initial states may be systematically biased; reweighting the prior over nucleon configurations with correlations could shift the extracted transport coefficients and tighten their uncertainty."],"forward_implications":["Ultra-central $v_n\\{2\\}$ values can act as a probe of nucleon-nucleon correlations in heavy nuclei, since the minimum-separation parameter directly controls the flow pattern.","Extracted values of $\\eta/s$ depend on the assumed nuclear structure; using the conventional Woods-Saxon distribution may bias the shear-viscosity estimate.","Sub-nucleon structure is not a negligible correction: adding sub-nucleon fluctuations changes the eccentricities and flow harmonics enough to matter for the puzzle.","The modification does not spoil the centrality dependence of bulk observables ($dN_{\\rm ch}/d\\eta$, $p_T$ spectra), so the $d_{\\min}$ knob is consistent with global calibrations.","The higher harmonics $v_4$ and $v_5$ remain compatible with the data under the modified initial conditions, so the improvement is not achieved by distorting only the low-order harmonics."],"supporting_citations":[{"why":"Provides the earlier demonstration that nucleon-nucleon correlations affect collective flow in ultra-central collisions, motivating the $d_{\\min}$ hypothesis.","marker":"[16]"},{"why":"CMS measurement of dihadron correlations in ultra-central PbPb; supplies the $v_n\\{2\\}$ data the paper targets.","marker":"[17]"},{"why":"The TRENTo initial-condition model whose rejection sampling the paper modifies to enforce $d_{\\min}$.","marker":"[51]"},{"why":"The CLVisc (3+1)-dimensional viscous hydrodynamic code used for all final-state flow calculations.","marker":"[9]"},{"why":"Sub-nucleon Glauber treatment that supplies the $N_c=3$ constituent-quark sub-nucleon structure.","marker":"[27]"},{"why":"ALICE pseudorapidity distributions used to calibrate the hydrodynamic parameters.","marker":"[57]"},{"why":"ALICE transverse momentum spectra used to check the calibration.","marker":"[58]"},{"why":"ATLAS differential flow data used to validate $v_2$, $v_3$, $v_4$ in semi-central classes.","marker":"[59]"}],"fun_headline_variants":["Nucleon spacing and sub-nucleon spots fix ultra-central flow puzzle","Initial nuclear structure resolves QGP viscosity puzzle in ultra-central collisions","Tweaking nucleon positions closes gap in heavy-ion flow predictions","Simulations match CMS with modified nuclear structure in ultra-central Pb-Pb","Nuclear geometry key to solving ultra-central flow ratio puzzle"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The physical relevance of the minimum-separation parameter: the paper states it is a proof-of-principle knob and that real $^{208}$Pb contains short-range-correlated pairs closer than 1.4 fm; if the true nucleon distribution is not as uniform, the improved $v_2/v_3$ agreement would be an artifact of the toy sampling.","fun_headline_variants_meta":{"raw":{"variants":["Nucleon spacing and sub-nucleon spots fix ultra-central flow puzzle","Initial nuclear structure resolves QGP viscosity puzzle in ultra-central collisions","Tweaking nucleon positions closes gap in heavy-ion flow predictions","Simulations match CMS with modified nuclear structure in ultra-central Pb-Pb","Nuclear geometry key to solving ultra-central flow ratio puzzle"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000236,"raw_usage":{"total_tokens":1540,"prompt_tokens":1018,"completion_tokens":522,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":634,"completion_tokens_details":{"reasoning_tokens":432}},"tokens_in":634,"tokens_out":522,"duration_ms":5739,"temperature":1.0,"reasoning_tokens":432,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:58:20.453234+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the two-nucleon separation distribution in $^{208}$Pb (from ab initio nuclear-structure calculations or from electron-scattering and short-range-correlation experiments). If it shows a substantial population of pairs closer than 1.4 fm, rerun the TRENTo+CLVisc calculation with that realistic distribution; loss of the $v_2\\{2\\}/v_3\\{2\\}$ agreement would falsify the claim that such uniformity is the physical cause.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the earlier demonstration that nucleon-nucleon correlations affect collective flow in ultra-central collisions, motivating the $d_{\\min}$ hypothesis."},{"cited_title":"Chatrchyan, et al., Studies of Azimuthal Dihadron Correlations in Ultra-Central PbPb Collisions at√sNN = 2.76 TeV, JHEP 02 (2014) 088","cited_arxiv_id":null,"evidence_quote":"CMS measurement of dihadron correlations in ultra-central PbPb; supplies the $v_n\\{2\\}$ data the paper targets."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The TRENTo initial-condition model whose rejection sampling the paper modifies to enforce $d_{\\min}$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The CLVisc (3+1)-dimensional viscous hydrodynamic code used for all final-state flow calculations."},{"cited_title":"Loizides, Glauber modeling of high-energy nuclear col- lisions at the subnucleon level, Phys","cited_arxiv_id":null,"evidence_quote":"Sub-nucleon Glauber treatment that supplies the $N_c=3$ constituent-quark sub-nucleon structure."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"ALICE pseudorapidity distributions used to calibrate the hydrodynamic parameters."},{"cited_title":"Abelev, et al., Centrality Dependence of Charged Par- ticle Production at Large Transverse Momentum in Pb– Pb Collisions at √sNN = 2.76 TeV, Phys","cited_arxiv_id":null,"evidence_quote":"ALICE transverse momentum spectra used to check the calibration."},{"cited_title":"Aad, et al., Measurement of the azimuthal anisotropy for charged particle production in√sNN = 2.76 TeV lead- lead collisions with the ATLAS detector, Phys","cited_arxiv_id":null,"evidence_quote":"ATLAS differential flow data used to validate $v_2$, $v_3$, $v_4$ in semi-central classes."}],"review_version":1}