{"id":"dac3c1d8-faa6-4a9f-9ab8-c38be8802c5e","arxiv_id":"2412.08431","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Using Gogny-HFB calculations, the authors show that collective inertias computed with the one-body part of the particle number fluctuation operator closely match those from the pairing gap operator, supporting the use of ΔN^2 as a fission collective coordinate.","lead":"This paper compares three ways of including pairing correlations in nuclear fission calculations and shows that particle number fluctuations and the pairing gap can be used interchangeably as collective variables. It also shows that artificially strengthening the pairing interaction lowers fission barriers and half-lives, in contrast to the other two methods.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The ΔN^2 one-body inertia benchmark is partly circular: 'Δ-constrained' states are sign-flipped ΔN^2 states (Appendix A), so the Fig. 8 agreement does not independently validate the one-body reduction.","rationale":"The paper's strongest claim relies entirely on the comparison in Fig. 8. The authors are honest about the Appendix A recipe, but that recipe means the 'Δ' coordinate is not an independent constraint. In the BCS/HFB canonical basis, flipping v_k signs changes the pairing tensor but not the density, so the set of states explored is essentially the ΔN^2-constrained family (up to discrete sign choices). The observed agreement between the ΔN^2 one-body inertia and the Δ inertia is then a near-necessary consequence of the quasi-linear relation between the two order parameters (Fig. 9), rather than a test of the central approximation. This does not refute the claim—the one-body reduction may well be accurate—but it removes the stated evidence. A direct Δ-constrained calculation or a full two-body cranking evaluation would settle it. The reader's conditional verdict already anticipates a missing independent benchmark; our analysis identifies the specific mechanism (the sign-flip recipe) that makes the benchmark circular, so the condition should be explicitly tied to repeating the calculation with a genuine Δ constraint. The QF analysis and half-life comparisons are not the load-bearing part of the central claim; they are consistent with prior literature and are not affected by this concern. Hence we do not ask for rejection, but for an added validation step before the conclusion is stated as 'solid ground.'","tokens_in":12208,"tokens_out":6534,"duration_ms":69171,"concrete_test":"Perform self-consistent HFB calculations with a genuine Lagrange-multiplier constraint on Δ = G Σ u_k v_k, maintaining the sign of v_k consistently per K block (e.g., by monitoring the K-block contributions to the pairing tensor and applying a small symmetry-restoring field or a positive-definite constraint on Σ |u_k v_k| with the appropriate sign assignment). Recompute the (Q20, Δ) PES, least-energy path, and B_eff. If the resulting energies and inertias differ from the Appendix A sign-flip recipe beyond a few percent, the Fig. 8 agreement is a reparameterization artifact and the one-body ΔN^2 inertia claim is unsubstantiated. Alternatively, compute the ΔN^2 inertia without the one-body reduction by finite differences of the constrained HFB states in the ΔN^2 direction, using the full two-body operator in the moments; compare with the one-body-part result.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (Sec. III.E) is that effective inertias from the one-body part of ΔN^2 agree with those from Δ, justifying the one-body cranking formula for a two-body operator. But the validation is not independent. As Appendix A explains, direct Δ-constrained HFB is avoided because v_k signs within K blocks can flip; instead, the authors solve with a ΔN^2 constraint and then flip the signs of v_k to make their contribution to Δ positive. Thus the states labelled by Δ are the ΔN^2-constrained states (with sign-adjusted occupancies), not states generated by constraining the one-body operator Δ. The PES and inertias in the (Q20, Δ) plane are therefore a reparameterization of the (Q20, ΔN^2) family, not an independent benchmark. Since Fig. 9 shows an almost linear Δ–ΔN^2 relation, the near-unity overlap and the close agreement of the two effective inertias can be understood as a Jacobian-scaled consistency between two coordinates on the same states; it does not test whether the one-body reduction of ΔN^2 reproduces the inertia one would obtain from a truly two-body cranking response. A different HFB state with the same Δ but different density could give different inertias. Consequently, the 'solid ground' claimed in the conclusions is weaker than stated and requires an independent check.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies how three different pairing-related quantities—the pairing gap Δ, the particle-number fluctuation ΔN^2, and a pairing-strength quenching factor QF—affect the static mean-field description of fission. Using HFB calculations with the Gogny D1S interaction, the authors compute potential energy surfaces, pairing energies, collective inertias, and WKB spontaneous-fission half-lives for 262Rf, 260Fm, and chains of U and Pu isotopes. The main methodological conclusion is that Δ and ΔN^2 can be used interchangeably as collective variables, and specifically that the one-body reduction of the two-body operator ΔN^2 in the perturbative cranking inertia is validated by comparison with the Δ-constrained inertias. The QF study shows that increasing the pairing strength lowers barriers and inertias, shortening half-lives, with an optimal QF around 1.05–1.10 for the D1S parametrization.","tokens_in":12533,"tokens_out":6441,"duration_ms":72375,"significance":"If the central validation claim held, the paper would provide useful justification for using ΔN^2 as a collective coordinate in fission studies, an approach used in several recent publications. The work has clear strengths: state-of-the-art HFB calculations with Gogny D1S, a consistent comparative framework for three pairing prescriptions, standard WKB half-life methodology, and a concrete empirical trend in half-lives versus QF. The paper is also honest about the numerical difficulty of constraining Δ directly, devoting Appendix A to the sign-of-v_k issue. However, as argued below, the central benchmark is not independent, so the significance of the main claim is currently weaker than the conclusions state.","major_comments":[{"comment":"The benchmark underlying the central claim is not independent. As Appendix A explains, the “Δ-constrained” HFB states are not obtained by constraining the one-body operator Δ; they are ΔN^2-constrained states with the sign of v_k adjusted block by block so that the contribution to Δ is positive. Because Fig. 9 shows Δ and ΔN^2 are almost linearly related over the entire (Q20, ΔN^2) grid, the agreement of the effective inertias in Fig. 8 is largely a consistency check between two coordinate labels on the same state family. It does not test whether the one-body part of ΔN^2 reproduces the inertia of the full two-body operator, since no independent states with the same Δ but different densities are involved. The conclusion in Sec. IV that this “gives a solid ground to the use of ΔN^2 as collective variable” is therefore stronger than the evidence. A direct test would be to compute the cranking inertia of ΔN^2 retaining the two-body terms (or using the full ATDHFB response) and compare it with the one-body-part result.","section":"Sec. III.E and Appendix A"},{"comment":"The central quantitative claim of “very good agreement” is supported only by visual inspection. No deviation metric is provided, such as the relative RMS difference between solid and dashed curves at matched Δ values, or the average absolute difference as a function of Q20 and Δ. Because this figure is the main evidence for the paper’s central conclusion, a quantitative measure is necessary to judge whether the agreement is actually satisfactory and whether it improves with pairing strength as the text asserts.","section":"Fig. 8"},{"comment":"The definition of the path used in the effective inertia is ambiguous. If Δ or ΔN^2 is held fixed along the Q20 path, then dq_pair/ds = 0 and Beff reduces to B_Q20Q20; if instead the path is a least-energy or least-action path in the two-dimensional space, this should be stated explicitly and the derivative dΔ/ds should be provided. The comparison in Fig. 8 is difficult to interpret without this specification, because the magnitude of the pairing-coordinate contribution depends on which path is actually followed.","section":"Sec. III.E, Eq. (5)"}],"minor_comments":[{"comment":"The text contains several typos and grammatical slips, including “Univeristy” in the affiliations, “Bohadilla” for Boadilla, “paring strength” in the Fig. 5 caption, and “spontaneous fission half-live” in Sec. IV.","section":"Throughout"},{"comment":"The phrase “the constrain on given ⟨ΔN^2⟩” should read “the constraint on a given ⟨ΔN^2⟩.”","section":"Fig. 9 caption"},{"comment":"The sentence “both following a quadratic trend as can be observed in Fig 1 in Ref [15]” is unclear; please specify which quantity follows the quadratic trend and which panel of Ref. [15] is meant.","section":"Sec. III.B"},{"comment":"The notation in Eq. (5) uses s for the path coordinate but the text identifies s with Q20; please make the path parametrization explicit so that the derivatives dq_i/ds are unambiguous.","section":"Eq. (5)"},{"comment":"The statement that ΔN^2 is “proportional to the pairing gap” is stronger than what Fig. 9 supports; the relation is approximately linear over the studied range, so the wording should be softened accordingly.","section":"Sec. III.B"}],"recommendation":"major_revision","confidential_remarks":"The paper is competently executed and the systematic comparison of Δ, ΔN^2, and QF is a useful contribution. However, the paper’s central methodological claim in Sec. III.E rests on a comparison that Appendix A itself shows is not independent: the Δ-constrained states are sign-adjusted ΔN^2-constrained states. The agreement in Fig. 8 is therefore less probative than the conclusions imply. I would ask the authors to provide an independent test of the one-body reduction, or at least to substantially soften the claim and supply quantitative agreement metrics. With those changes, the paper would be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nTwo things you should know about this paper. First, it is a competent and genuinely useful systematic comparison of the three ways pairing is usually handled in static fission calculations: constraining the pairing gap Δ, constraining particle-number fluctuation ΔN^2, and scaling the pairing strength with a quenching factor QF. Second, its central claim—that Δ and ΔN^2 are interchangeable as collective coordinates, validated by comparing their cranking inertias—is weaker than the text admits, because the comparison is not independent.\n\nWhat is new and good: the QF analysis reconciles an apparent contradiction in the literature, where increasing pairing strength lowers fission barriers while imposing Δ or ΔN^2 constraints raises them. The paper shows both effects are real and distinct. The half-life systematics across U, Pu, Cm, Fm, and Rf isotopes are useful, and the authors are explicit that the optimal QF is an empirical fit, not a prediction. Appendix A's discussion of the sign ambiguity of the v_k occupancies under a Δ constraint is a solid, nontrivial observation about HFB/BCS that deserves to be widely known.\n\nThe soft spot is the validation of the one-body cranking approximation for ΔN^2. As Appendix A makes clear, the states labelled by Δ are actually ΔN^2-constrained states whose v_k signs are flipped to give positive Δ. So Figure 8 compares two different one-body operators evaluated on the same set of wavefunctions, not the one-body approximation against an independent two-body cranking result. The near-linear Δ–ΔN^2 relation and the near-unity overlaps then do most of the work. The conclusion that this 'gives a solid ground' for using ΔN^2 as a collective variable is overstated. That said, the practical recommendation—that ΔN^2 is a reasonable pairing coordinate for fission—is probably correct for the high-pairing regime that matters for least-action paths, and the authors do note that the agreement improves as pairing grows.\n\nWho this is for: fission theorists and code developers using collective inertial models. It deserves a serious referee; the calculations are mature, the dataset is useful, and the sign-ambiguity result alone is publishable. The referee should require an explicit statement that the inertia comparison is a consistency check on the same states, not an independent benchmark, and ideally a test against a full two-body cranking response.\n\nMy recommendation: send it to peer review with a request for revision that directly addresses the circularity concern. The paper will be better for it.","headline":"Solid systematic study of pairing in fission; the headline Δ–ΔN^2 inertia validation is partly circular and the interchangeability claim is stronger than the evidence.","tokens_in":13010,"tokens_out":2956,"would_cite":true,"duration_ms":30852,"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":"Particle-number fluctuations match the pairing gap as a fission coordinate, so ΔN² can replace Δ in static fission calculations.","keywords":["spontaneous fission","pairing gap","particle-number fluctuations","collective inertia","cranking approximation","HFB","Gogny D1S","fission barrier"],"falsifier":"Compute the collective inertia for ΔN² without the one-body reduction—for example, with a generator-coordinate method using the full two-body fluctuation operator or with a full two-body cranking formula—along the same fission path in 262Rf; if the resulting inertia differs substantially from the one-body cranking value, the interchangeability of Δ and ΔN² would fail exactly where the paper claims it holds.","tokens_in":12048,"feed_emoji":"⚛️","tokens_out":5521,"duration_ms":52824,"temperature":0.7,"pith_summary":"This paper studies how three different ways of representing pairing correlations—the pairing gap Δ, the particle-number fluctuation ΔN², and the quenching factor QF that scales pairing strength—affect static fission predictions. Its central result is that Δ and ΔN² behave nearly identically as collective coordinates: both raise and broaden the fission barrier while reducing the collective inertia, and the effective cranking inertia computed from the one-body part of ΔN² agrees well with that computed from the one-body operator Δ. This agreement justifies the common practice of using ΔN² as a dynamical variable in fission calculations, even though it is formally a two-body operator. The quenching factor acts differently: it lowers the barrier and inertia together, monotonically shortening half-lives, with a D1S optimum near QF = 1.05–1.10.","feed_headline":"Particle-number fluctuations can replace the pairing gap in fission","feed_subtitle":"Matching collective inertias make ΔN² a reliable fission coordinate, and a 5–10% pairing boost best fits experiment.","key_machinery":"The central object is the effective collective inertia B_eff(s) = Σ_{ij} B_ij (dq_i/ds)(dq_j/ds) along the fission path s = Q20, where B_ij is the inverse of the perturbative cranking mass tensor. The cranking formula B(Q20) = M_{-3}/M_{-1}^2 assumes a one-body collective operator, so for the two-body operator ΔN² the paper uses its quasiparticle one-body part and compares the resulting inertia with that of the genuine one-body operator Δ. This comparison validates the one-body approximation and makes ΔN² a usable dynamical coordinate. A supporting piece of machinery is the Appendix-A sign-fixing recipe, which ensures a consistent sign for the occupation factors when constraining Δ, because otherwise individual K-blocks can contribute negative terms to the gap.","core_discovery":"Using the HFB method with the Gogny D1S interaction for 262Rf, 260Fm, and several U and Pu isotopes, the paper shows that constraining on Δ or on ΔN² produces nearly identical potential-energy surfaces, pairing energies, and collective inertia profiles along the least-energy fission path. The load-bearing comparison is the effective inertia B_eff(s) along Q20: for fixed values of Δ and the corresponding ΔN², the inertia obtained from ΔN² with the one-body part of the operator in the cranking formula agrees with the inertia obtained from the genuine one-body operator Δ, with the agreement improving for larger pairing correlations. The authors conclude that ΔN² can be used interchangeably with Δ as a collective variable in fission studies, and that the one-body assumption for its inertia is solid. They also quantify how the quenching factor differs: increasing QF strengthens the attractive pairing interaction, lowering the barrier and inertia and reducing half-lives, with the best reproduction of experimental half-lives at QF between 1.05 and 1.10.","pith_inferences":["Because ΔN² is sign-independent while Δ requires an ad hoc sign-fixing recipe, ΔN² may be the more robust order parameter for constrained HFB fission calculations; the paper's agreement suggests this robustness carries over to inertias.","The one-body reduction of ΔN² could be tested non-perturbatively with a generator-coordinate inertia computed from the full two-body operator; if the cranking agreement persists, the result would validate a much cheaper dynamical treatment.","The QF optimum near 1.05–1.10 for D1S implies that the fitted pairing strength of this interaction is slightly too weak for fission observables; re-fitting the interaction itself could remove the need for a manual quenching factor.","The same Δ-versus-ΔN² inertial comparison could be extended to asymmetric fission paths and to octupole degrees of freedom, where the paper already notes the QF effect is weaker."],"forward_implications":["Δ and ΔN² can be used interchangeably as collective coordinates in static fission calculations without changing the predicted barrier, inertia, or half-life.","The one-body cranking formula can be applied to ΔN², making dynamic least-action fission calculations with pairing fluctuations tractable.","Increasing pairing via Δ or ΔN² broadens and raises the barrier but lowers the inertia; the net effect is shorter spontaneous-fission half-lives.","Scaling pairing strength via QF lowers both barrier and inertia, monotonically reducing half-lives; the D1S force needs QF ≈ 1.05–1.10 to match measured half-lives, with the effect concentrated on the symmetric first barrier.","Since the agreement improves at larger Δ values, the effective inertia approximation is most reliable in the high-pairing region that dominates the least-action path."],"supporting_citations":[{"why":"Supplies the perturbative cranking formula for the collective inertia that the paper applies to both Δ and the one-body part of ΔN².","marker":"[40]"},{"why":"Introduced ΔN² as a collective degree of freedom in fission studies, the approach whose validity the present comparison tests.","marker":"[15]"},{"why":"Established the 1/Δ² dependence of collective inertia on the pairing gap, the baseline behavior that the ΔN² results reproduce.","marker":"[8]"},{"why":"Used ΔN² as a dynamic pairing coordinate in fission, providing the context for validating one-body inertias.","marker":"[16]"},{"why":"Extended ΔN² to beyond-mean-field fission studies, motivating the need to justify its one-body cranking inertia.","marker":"[18]"}],"fun_headline_variants":["Particle-number fluctuations can replace pairing gap in fission","ΔN² matches Δ as fission coordinate, inertias agree","Fission: Use ΔN² instead of pairing gap for collective inertia","Quenching factor 5–10% boost best reproduces fission half-lives","Collective inertias confirm ΔN² as viable fission variable"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper assumes that the perturbative cranking inertia computed from the one-body part of ΔN² is a faithful proxy for the inertia of the true two-body operator, with no independent non-perturbative benchmark.","fun_headline_variants_meta":{"raw":{"variants":["Particle-number fluctuations can replace pairing gap in fission","ΔN² matches Δ as fission coordinate, inertias agree","Fission: Use ΔN² instead of pairing gap for collective inertia","Quenching factor 5–10% boost best reproduces fission half-lives","Collective inertias confirm ΔN² as viable fission variable"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001379,"raw_usage":{"total_tokens":5522,"prompt_tokens":816,"completion_tokens":4706,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":432,"completion_tokens_details":{"reasoning_tokens":4616}},"tokens_in":432,"tokens_out":4706,"duration_ms":35873,"temperature":1.0,"reasoning_tokens":4616,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T17:50:23.346213+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the collective inertia for ΔN² without the one-body reduction—for example, with a generator-coordinate method using the full two-body fluctuation operator or with a full two-body cranking formula—along the same fission path in 262Rf; if the resulting inertia differs substantially from the one-body cranking value, the interchangeability of Δ and ΔN² would fail exactly where the paper claims it holds.","supporting_citations":[{"cited_title":"Baran, J","cited_arxiv_id":null,"evidence_quote":"Supplies the perturbative cranking formula for the collective inertia that the paper applies to both Δ and the one-body part of ΔN²."},{"cited_title":"Sadhukhan, J","cited_arxiv_id":null,"evidence_quote":"Used ΔN² as a dynamic pairing coordinate in fission, providing the context for validating one-body inertias."},{"cited_title":"Rodríguez-Guzmán, L","cited_arxiv_id":null,"evidence_quote":"Extended ΔN² to beyond-mean-field fission studies, motivating the need to justify its one-body cranking inertia."}],"review_version":1}