{"id":"19a3922f-1aad-4543-97e8-a5f1e67a5373","arxiv_id":"2501.00682","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Short-range S-wave contact terms in the chiral nucleon-nucleon force substantially change computed quadrupole collectivity in 6Li and 12C by shifting surface oscillations within one dominant nuclear shape.","lead":"Nuclear theorists show that two short-range forces between protons and neutrons can noticeably alter the shape oscillations that drive collective rotation-like motion in light nuclei. The result suggests that measurements of nuclear shapes could help pin down the detailed short-distance physics of the strong nuclear force.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 12C leg of the central claim is not validated: the 14-shape, 8ℏΩ model space captures only ~80% of the 2+1 state, and the shape-selection independence shown for 6Li (Supplemental Figs. 7–8) is never demonstrated for 12C.","rationale":"The reader's weakest assumption correctly identifies the pivotal gap. The paper's headline is explicitly two-nucleus ('6Li and 12C'), and the quantitative 'up to ~50%' claim for the 12C 2+1 state is computed in a truncated symmetry-adapted space that captures only about 80% of the state. The authors do provide convergence evidence for 6Li in the supplement (Figs. 7–9), and that evidence is genuine and relevant: it shows that for 6Li the sensitivity pattern is stable across shape selection and model-space size. But the same validation is absent for 12C, where the complete 8ℏΩ space is far larger and where the shape composition is less dominated by a single irrep (60–70%, not 82–86%). Since the central claim depends on the 12C response being large and mechanistically clean, the absence of a complete-space check is the most load-bearing soft spot. The secondary issue—the abstract's 'without changing that shape's overall contribution' being an overstatement for 12C—would not by itself change the verdict, but it reinforces the need for caution. I therefore keep the reader's CONDITIONAL verdict: the 6Li result appears well supported, but the 12C leg requires the proposed complete-space (or at least shape-converged) check before the claim can be accepted as stated. This is not a rejection: the paper is internally consistent, uses a credible method with real convergence checks for 6Li, and the sensitivity pattern is reproducible in principle. The concern is a missing validation, not a demonstrated error.","tokens_in":14556,"tokens_out":10466,"duration_ms":100981,"concrete_test":"Compute the 12C 2+1 quadrupole moment in the complete 8-ℏΩ no-core shell-model space (all Sp(3,R) shapes, no truncation) for the NNLOopt LEC set and for C_3S1^LO and C_1S0^NLO each shifted by ±10%, holding all other LECs fixed (five diagonalizations). Extract dQ2/dC for each contact and compare with the 14-shape slopes in Fig. 4(c). If either complete-space slope differs from the truncated slope by more than the ~15% many-body uncertainty quoted in the paper, or changes sign, the 12C part of the central claim is not established; if the slopes agree within that band, the truncation concern is resolved. If a full 8ℏΩ complete run is impractical, the same comparison at 6ℏΩ with a 14-shape-to-all-shape sequence, plus a 14-shape-to-all-shape check at 8ℏΩ for the two endpoints, would provide a quantitative bound on the omitted-shape contribution.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing assumption is that the 14-shape, 8-ℏΩ symmetry-adapted space used for the 12C 2+1 state (main text: 'fourteen shapes that account for about 80% of the state') preserves the magnitude and sign of the quadrupole response to the S-wave contacts C_3S1^LO and C_1S0^NLO. The authors explicitly validate shape-selection independence only for 6Li: Supplemental Fig. 7 compares 1-shape, 14-shape, and all-shape results at 8ℏΩ, and Fig. 8 extends the 14-shape calculation to 10 and 12 shells for 6Li. No equivalent complete-space or larger-space check is given for 12C. Because Q2 is a sum over Sp(3,R) irreps, the omitted ~20% of the 12C wave function contributes additively, and nothing bounds the size or LEC-dependence of its Q2 matrix elements; the 80% overlap does not imply 80% of the response. The resemblance of the 12C sensitivity pattern to 6Li (Fig. 3) is suggestive but not a validation. A secondary, related overstatement is in the abstract: for 12C the dominant oblate shape varies from ~60% to ~70% and the second shape from ~5% to ~15% across the 32 samples (Supplemental Fig. 6c), so 'without changing that shape's overall contribution' is not literally true for carbon, and some of the Q2 response may be shape mixing rather than purely internal surface oscillations. The first issue is the one that would falsify the quantitative claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript reports symmetry-adapted no-core shell model (SA-NCSM) calculations for the 6Li 1+ ground state, 6Li 3+1, and 12C 2+1 states, in which the 14 low-energy constants (LECs) of the NNLOopt chiral NN potential are varied within ±10% and analyzed with Sobol' global sensitivity analysis using 300,000 samples. The central claim is that quadrupole collectivity, as measured by the quadrupole moment Q2, is strongly influenced by two S-wave contact couplings, C_3S1^LO and C_1S0^NLO, not by changing the overall nuclear shape composition but by redistributing polar versus equatorial surface oscillations within the dominant Sp(3,R) shape. The authors support this with first- and total-order sensitivity indices, single-LEC response curves, shape-probability plots, and an energy decomposition that links the Q2 response to kinetic- and potential-energy changes.","tokens_in":14938,"tokens_out":7465,"duration_ms":72938,"significance":"If correct, this is a significant result. It identifies a direct, sizable footprint of short-range NN contact physics on a long-range collective observable, with implications for chiral-potential fitting, uncertainty quantification, and the microscopic interpretation of nuclear deformation. The paper's strengths include a very large 300,000-evaluation GSA, explicit convergence checks for 6Li against complete model spaces and larger harmonic-oscillator spaces, and a consistent energy decomposition (Fig. 5) that supports the proposed polar/equatorial mechanism. The study is computational and contains no fitted observables; the LEC variation range is an input assumption. The result is falsifiable: the predicted Q2 response to the two S-wave contacts should persist in complete-space and larger-space calculations for 12C.","major_comments":[{"comment":"The 12C leg of the central claim is not validated. The fourteen-shape, 8-ℏΩ model space used for the 12C 2+1 state accounts for about 80% of the state, and the authors explicitly validate shape-selection independence for 6Li only (Supplemental Figs. 7 and 8). Because Q2 is additive over Sp(3,R) irreps and the omitted ~20% of the wave function contributes additively to the expectation value, an 80% overlap does not bound the magnitude or LEC dependence of the omitted contribution; the resemblance between the 12C and 6Li sensitivity patterns in Fig. 3 is suggestive but not a substitute. Please provide a complete-space or larger-Nmax check for the 12C response, or explicitly restrict the 'up to ~50%' claim to 6Li.","section":"Sensitivity of collective observables; Supplemental Figs. 7-8"},{"comment":"The abstract's claim that the LEC variations alter Q2 'without changing that shape's overall contribution within the nucleus' is not literally true for 12C. Supplemental Fig. 6c shows the dominant oblate shape varying from about 60% to 70% and the second shape from about 5% to 15% across the samples. Some of the Q2 response in 12C may therefore be due to shape mixing rather than purely internal surface oscillations of a fixed dominant shape. Please rephrase the abstract and Sec. 'Short-range footprints on long-range physics' to state that shape probabilities change modestly, and quantify the split between shape-probability changes and intra-shape oscillation changes.","section":"Abstract; Supplemental Fig. 6c"},{"comment":"The Sobol' indices are prior-dependent. The uniform, independent ±10% sampling of all 14 LECs is an input assumption, not a physical uncertainty estimate, and the first-order indices are normalized to the total variance under this prior. The 'up to ~50%' magnitude is therefore conditional on the chosen range. Please either justify the ±10% range (e.g., from LEC covariance or naturalness) or state explicitly that the quantitative claim is a response to this defined range; the response slopes in Fig. 4, which are range-independent, should be presented as the primary quantitative result.","section":"Sensitivity of collective observables; Fig. 3"}],"minor_comments":[{"comment":"The gray symbols for '300 simultaneously varied LEC samples' are not described in the text; please clarify whether these are the same samples as in the GSA and how the reader should compare them with the one-at-a-time curves.","section":"Fig. 4 caption and text"},{"comment":"The caption uses 'probability amplitudes' for quantities that appear to be probabilities; please use consistent terminology.","section":"Fig. 2 caption"},{"comment":"The phrase 'without loss of generality, our calculations use a single ℏΩ = 15 MeV' is not strictly lossless for finite-space observables; please add a brief caveat or reference showing that the LEC-response conclusions are ℏΩ-insensitive.","section":"Resilience of nuclear shapes"},{"comment":"The absence of three-nucleon forces is noted to minimally affect binding energies, but the manuscript does not discuss their possible effect on the Q2 response; a sentence noting the caveat (cf. Ref. [43]) would be helpful.","section":"Resilience of nuclear shapes"},{"comment":"Supplemental Fig. 7(d) shows that for the 6Li 1+ ground state the sensitivity indices differ between the complete space and the 14-shape space for some LECs (e.g., C3P2, c3); the main text's statement that the sensitivity patterns are 'practically unchanged' should acknowledge that this holds for the collective states and is less clean for the weakly-collective ground-state Q2.","section":"Supplemental Fig. 7(d)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a strong candidate for publication if the 12C truncation issue is resolved. The missing validation is in principle a modest additional computation, not a conceptual flaw. The main concern is the abstract's overstatement for 12C and the lack of a complete-space check for that nucleus. I would not reject the manuscript on the current evidence."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What you should know: this paper does something new and does it mostly well. It presents the first global sensitivity analysis of quadrupole moments in ab initio nuclear structure, using 300,000 SA-NCSM evaluations with a fixed chiral potential and no fitted parameters. The result is that two S-wave contact couplings, C_3S1^LO and C_1S0^NLO, can change Q2 of the 6Li 3+ and 12C 2+ states by up to ~50%, and the mechanism is not shape mixing but a redistribution between polar and equatorial surface oscillations within the dominant Sp(3,R) shape. That is a genuinely new connection between short-range NN physics and collective observables, and it is backed by real numerical evidence: convergence checks against complete model spaces and larger shell spaces for 6Li, a consistent energy decomposition in Fig. 5, and total-order indices that match first-order indices, suggesting the LEC effects are decoupled. I believe the 6Li leg of the claim is solid.\n\nThe soft spots are in proportion. The biggest one is the 12C leg. The 14-shape, 8ℏΩ space captures only about 80% of the 2+ state, and the shape-selection independence is validated only for 6Li. The stress-test note is right: an 80% wave-function overlap does not imply 80% of the Q2 response, since the omitted shapes contribute additively and nothing bounds their LEC dependence. A complete-space check for 12C, or at least a larger-space check, would close this. Second, the abstract says the shape's overall contribution is unchanged, but the supplement shows the dominant oblate shape in 12C varies from ~60% to ~70% and the second shape from ~5% to ~15% across samples. So for carbon some of the response may be shape mixing after all. That is an overstatement, not a fatal flaw. Third, the ±10% uniform sampling is artificial—real LEC uncertainties are correlated and nonuniform—but this is a sensitivity analysis, not a prediction, so I do not weigh it heavily. Minor circularity: the polar/equatorial decomposition uses the same Sp(3,R) quantum numbers that define Q2, so the mechanistic story is partly built into the basis, but the direction and size of the effect are still a computational observation.\n\nWho is this for? Nuclear structure people working on chiral EFT and collectivity, and anyone doing uncertainty quantification in ab initio methods. It deserves a serious referee. I would send it to review, and ask for a 12C validation or a bound on the omitted-shape contribution, plus a correction of the abstract's wording. The central claim is credible enough that the paper should be in the literature, with those caveats addressed.","headline":"A serious ab initio sensitivity study that credibly links short-range S-wave contacts to quadrupole collectivity in 6Li and 12C, with a genuinely new mechanism and one real soft spot: the 12C leg lacks full-space validation.","tokens_in":15460,"tokens_out":1511,"would_cite":true,"duration_ms":15889,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Short-range contact couplings of the nucleon-nucleon force can change nuclear quadrupole collectivity by up to roughly 50 percent by rebalancing surface oscillations within the dominant nuclear shape, without altering that shape's…","keywords":["nuclear collectivity","quadrupole moments","short-range nucleon-nucleon interactions","chiral effective field theory","symmetry-adapted no-core shell model","Sobol sensitivity analysis","nuclear shapes","surface oscillations"],"falsifier":"Run the same ±10% LEC variation analysis for the 12C 2+1 state in a complete (all-shape) model space spanning 8 or more harmonic-oscillator shells and compare the resulting quadrupole response; if the combined effect of C(LO)3S1 and C(NLO)1S0 drops well below 50 percent or changes sign relative to the 14-shape result, the paper's central claim for carbon fails.","tokens_in":14368,"feed_emoji":"⚛️","tokens_out":5747,"duration_ms":48135,"temperature":0.7,"pith_summary":"The paper claims that the short-range contact part of the nucleon-nucleon force, not just the long-range correlations traditionally credited with deforming nuclei, can substantially change nuclear collectivity. Using ab initio calculations for low-lying states of 6Li and 12C, it shows that varying two S-wave contact couplings changes electric quadrupole moments by up to roughly 50 percent. The effect arises because these couplings tip the balance between polar and equatorial surface oscillations within a single dominant nuclear shape, while leaving the shape composition nearly unchanged. This matters because it identifies a new, short-distance route to collective nuclear behavior and suggests quadrupole moments can help pin down the contact parameters of chiral potentials.","feed_headline":"Short-range nuclear forces shift collectivity by up to 50%","feed_subtitle":"Two S-wave couplings rebalance surface oscillations and shift quadrupole moments, linking short-range physics to collective motion.","key_machinery":"The load-bearing object is the Sp(3,R)-adapted 'shape' basis of the symmetry-adapted no-core shell model, which groups the many-body Hilbert space into subspaces that each represent a nuclear shape with a static deformation plus its dynamical surface oscillations (2ℏΩ particle-hole excitations). Decomposing a state this way lets the paper separate changes in shape composition from changes inside the dominant shape. The analysis machinery is Sobol' global sensitivity analysis applied to 300,000 samples of the fourteen low-energy constants, which yields first- and total-order indices showing that Q2 responds mainly to C(LO)3S1 and C(NLO)1S0, and this is confirmed by one-at-a-time response curves.","core_discovery":"On its own terms, the paper establishes that the LO 3S1 and NLO 1S0 S-wave contact couplings of the NNLOopt chiral NN potential act in opposition on the quadrupole moments Q2 of the 6Li 1+ ground state, 6Li 3+1, and 12C 2+1 states. Strengthening the attractive LO 3S1 coupling decreases Q2, while strengthening the repulsive NLO 1S0 coupling increases it; together they can change Q2 by up to about 50 percent for the collective states. The mechanism is not a rearrangement among nuclear shapes—the dominant shape's probability stays essentially fixed—but a redistribution within that shape's dynamical surface oscillations, with polar oscillations (along the symmetry axis) enhancing collectivity and equatorial oscillations suppressing it. The paper further finds that first-order and total-order Sobol' sensitivity indices nearly coincide, so the two couplings act largely independently, and that the quadrupole response is carried almost entirely by the dominant shape.","pith_inferences":["If the mechanism persists in larger spaces, one would expect the polar/equatorial balance to leave a signature in B(E2) transition strengths, not just static quadrupole moments; this is a natural next calculation the paper does not perform.","The same rebalancing argument suggests that other short-range operators—such as tensor or spin-orbit contacts—could have comparably large collective effects in heavier nuclei where shape coexistence is richer.","A practical extension would be to use the dominant-shape response as a cheap emulator for uncertainty quantification of Q2 across a chiral-potential family, since the paper shows the response lives almost entirely in one shape.","Experimental measurements of Q2 or B(E2) with percent-level precision in light nuclei could directly discriminate between chiral potentials that differ mainly in these two S-wave contacts."],"forward_implications":["Quadrupole moments of collective states become observables that can constrain the contact LECs of chiral potentials, since short-range coupling changes show up as sizable Q2 shifts.","The two S-wave contacts affect binding energies and quadrupole moments through similar physics, so fits of these LECs can in principle be informed by both types of data simultaneously.","For states whose Q2 is dominated by a single shape, the shape composition is not the only route to collectivity; surface-oscillation balance within that shape is an equally important lever.","Because the sensitivity pattern for 12C matches 6Li, the mechanism may be common across light deformed nuclei, and not specific to one isotope."],"supporting_citations":[{"why":"Supplies the NNLOopt chiral NN potential whose low-energy constants are varied in the sensitivity analysis.","marker":"[41]"},{"why":"Establishes the Sp(3,R) shape decomposition and the dominant-shape picture used to interpret the response.","marker":"[1]"},{"why":"Describes the symmetry-adapted no-core shell model framework that makes the 300,000 calculations feasible.","marker":"[25]"},{"why":"Provides the symmetry-guided shell-model methodology and shape basis underlying the calculations.","marker":"[26]"},{"why":"Gives the Sobol' variance-based sensitivity indices used for the global sensitivity analysis.","marker":"[23]"},{"why":"Provides the Saltelli estimation procedure for total-order sensitivity indices.","marker":"[24]"},{"why":"Earlier global sensitivity analysis showing binding-energy sensitivity to the same S-wave contacts in 16O, which the paper links to its Q2 finding.","marker":"[52]"},{"why":"Prior ab initio demonstration of collective modes in light nuclei that motivates the quadrupole-moment focus.","marker":"[68]"}],"fun_headline_variants":["Short-range nuclear forces stir surface oscillations, shift collectivity by 50%","Opposing S-wave couplings swing nuclear collectivity by 50%","Short-range contacts rebalance nuclear collectivity without shape change","Short-range nuclear forces redistribute surface oscillations, boosting collectivity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The nearly 50 percent quadrupole response is computed in a 14-shape model space for 12C that captures about 80 percent of the state, and the size and sign of that response is assumed to survive in the complete space—a convergence check performed for 6Li, not for 12C.","fun_headline_variants_meta":{"raw":{"variants":["Short-range nuclear forces stir surface oscillations, shift collectivity by 50%","Opposing S-wave couplings swing nuclear collectivity by 50%","Short-range contacts rebalance nuclear collectivity without shape change","Short-range nuclear forces redistribute surface oscillations, boosting collectivity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000786,"raw_usage":{"total_tokens":3427,"prompt_tokens":862,"completion_tokens":2565,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":478,"completion_tokens_details":{"reasoning_tokens":2492}},"tokens_in":478,"tokens_out":2565,"duration_ms":18848,"temperature":1.0,"reasoning_tokens":2492,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:45:49.448947+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same ±10% LEC variation analysis for the 12C 2+1 state in a complete (all-shape) model space spanning 8 or more harmonic-oscillator shells and compare the resulting quadrupole response; if the combined effect of C(LO)3S1 and C(NLO)1S0 drops well below 50 percent or changes sign relative to the 14-shape result, the paper's central claim for carbon fails.","supporting_citations":[{"cited_title":"Dytrych, K","cited_arxiv_id":null,"evidence_quote":"Establishes the Sp(3,R) shape decomposition and the dominant-shape picture used to interpret the response."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the Sobol' variance-based sensitivity indices used for the global sensitivity analysis."},{"cited_title":"Saltelli, P","cited_arxiv_id":null,"evidence_quote":"Provides the Saltelli estimation procedure for total-order sensitivity indices."},{"cited_title":"(#$)!!!\"(&#$)!!'#","cited_arxiv_id":null,"evidence_quote":"Prior ab initio demonstration of collective modes in light nuclei that motivates the quadrupole-moment focus."}],"review_version":1}