{"id":"feb8792f-80b2-462e-b65c-ab41cf3ede75","arxiv_id":"2608.03663","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A relativistic density functional based statistical model computes nuclear level densities and neutron resonance spacings for 67 even-even nuclei, with accuracy comparable to existing microscopic methods.","lead":"This paper develops a way to compute the density of excited nuclear states using relativistic mean-field theory with pairing and collective effects. It tests the method on 67 nuclei and finds it matches measured level densities and neutron resonance spacings about as well as other microscopic models.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Reported f_rms values may be in-sample: dU is fit to resonance spacings at Sn, then the same type of data is used to claim reproduction.","rationale":"The reader's verdict is CONDITIONAL, and I agree with that overall assessment. However, I identify the in-sample dU fit as the most load-bearing concern rather than the parity equipartition assumption. The dU fit is explicitly stated in Section II, and the validation in Section IIID uses the same kind of data at the same energy. Without a clear statement that the 67 nuclei were excluded from the fit, the f_rms values cannot be interpreted as independent tests of the model. The parity assumption is also a genuine limitation for the p-wave spacings, and both concerns warrant qualification of the central claim. Because the model still provides a systematic RDFT-based framework with reasonable performance for many nuclei, and because both concerns can be addressed with straightforward additional analyses (cross-validation for dU; parity-projected level densities for the parity assumption), a CONDITIONAL verdict remains appropriate. Thus I recommend no change to the reader's verdict.","tokens_in":17050,"tokens_out":11390,"duration_ms":130038,"concrete_test":"Obtain the code or ask the authors whether the 67 s-wave resonance spacings used to compute f_rms were also used to fit dU. If they were (or if this cannot be ruled out), perform a cross-validation: randomly split the 67 nuclei into a fitting set and a test set, refit dU on the fitting set, and compute f_rms on the test set; repeat at least 100 times. If the mean held-out f_rms is substantially larger than the reported 2.29 (e.g., >20% relative increase), the reported validation is in-sample and the predictive claim is not supported. If dU was fit on independent data, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section II (Eq. 18) states that dU = 6.27 MeV is determined by fitting experimental data at the neutron separation energy. Section IIID then reports f_rms = 2.29 (DD-ME2) and 2.56 (DD-PC1) for the 67 even-even nuclei with experimental s-wave neutron resonance spacings, and similar p-wave and combined values. The manuscript nowhere states that the dU fit excluded these 67 nuclei or used an independent dataset. If dU was adjusted to match these same resonance spacings, the f_rms values are in-sample measures of fit quality, not out-of-sample assessments of predictive power. The central claim in Section IV that the approach 'can reasonably reproduce' experimental spacings then rests on a potentially circular validation. This directly affects all reported f_rms values and the comparison to non-relativistic methods (f_rms = 2.14) and phenomenological models (f_rms = 1.78). The parity equipartition assumption is a separate concern for the p-wave extraction, but the in-sample dU fit is more load-bearing because it affects the primary quantitative evidence for the central claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a microscopic statistical model for nuclear level density (NLD) based on relativistic density functional theory (RDFT). It uses self-consistent single-particle levels from RDFT, finite-temperature BCS pairing, and rotational/vibrational collective enhancement factors, with a microscopically computed spin cut-off parameter. The model is applied to 67 even–even nuclei with experimental s-wave neutron resonance spacings, comparisons are made with phenomenological and non-relativistic microscopic models, and a detailed study of shape coexistence in 98Sr is presented. The central claim is that this RDFT-based statistical approach can reasonably reproduce experimental s- and p-wave neutron resonance spacings and NLDs, with f_rms = 2.29 (DD-ME2) and 2.56 (DD-PC1) for s-wave spacings.","tokens_in":17344,"tokens_out":2226,"duration_ms":26671,"significance":"If the reported validation is genuinely predictive, the paper would establish an RDFT-based statistical NLD method with a microscopically computed spin cut-off as a viable tool for practical applications, comparable to existing non-relativistic microscopic methods and useful for nuclei where phenomenological fits are unreliable. The systematic comparison of two covariant interactions, the treatment of deformation effects, and the concrete study of 98Sr are valuable. However, the paper's central quantitative evidence is weakened by a potential in-sample fitting of the damping parameter dU to the same type of neutron resonance data used in the validation, and the assumption of parity equipartition for p-wave spacings. These issues need to be addressed before the central claim is fully supported.","major_comments":[{"comment":"The damping parameter dU = 6.27 MeV is stated to be 'determined by fitting experimental data at the neutron separation energy.' The same paper then reports f_rms values for s-wave neutron resonance spacings (Sec. III.D) for 67 even–even nuclei. The manuscript does not state whether these 67 nuclei were excluded from the dU fit. If they were not excluded, the quoted f_rms values (2.29 and 2.56) are in-sample measures of fit quality, not out-of-sample tests of predictive power. This would directly affect the central claim in Sec. IV that the approach 'can reasonably reproduce experimental s- and p-wave neutron resonance spacings' and the comparison to non-relativistic (f_rms=2.14) and phenomenological (f_rms=1.78) methods. The authors must clarify the fitting protocol and, if necessary, perform an out-of-sample validation (e.g., leave-one-out or fit on an independent subset).","section":"Sec. III.D"},{"comment":"The p-wave resonance spacing analysis assumes that positive and negative parities are equally distributed. The paper does not compute parity-dependent level densities, despite the introduction's claim to treat both spin and parity dependence. Since the p-wave f_rms and the combined s+p f_rms depend on this assumption, the authors should either justify the equipartition assumption with concrete estimates for the nuclei considered or compute parity-dependent NLDs. If parity equipartition is inaccurate for deformed or mid-shell nuclei, the reported p-wave f_rms (2.84/2.93) and combined f_rms (2.41/2.64) could shift significantly.","section":"Sec. III.D"},{"comment":"The perpendicular spin cut-off parameter σ²⊥ in Eq. (14) uses a simplified moment of inertia I⊥,k = (2/5) m R², with R presumably the nuclear radius. This rigid-body estimate is not derived from the RDFT density or the actual mass distribution. Since the spin cut-off is a claimed improvement in this work, the sensitivity of σ²⊥ to the choice of R (e.g., from RDFT vs. empirical formula) should be discussed. This is not fatal, but it is a load-bearing input for the resonance-spacing predictions and should be validated.","section":"Sec. II, Eq. (14)"}],"minor_comments":[{"comment":"Typo: 'thes- andp-wave' should be 'the s- and p-wave'.","section":"Abstract"},{"comment":"The f_rms comparison involves different sample sizes (67 vs. 278 vs. 295 nuclei). A statistical measure that accounts for sample size or a match to the same nuclei would make the comparison more meaningful.","section":"Sec. III.D"},{"comment":"The symbol Pa′ in Eq. (19) is not defined in the text; indicate that it is the deuteron pairing energy and give its source.","section":"Sec. III.A, Eq. (19)"},{"comment":"Figure 3 legend is difficult to read; please increase font size and use distinct line styles for the six deformations.","section":"Sec. III.B, Fig. 3"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the journal's scope and presents a useful methodological contribution. The main reason for major revision is the unresolved in-sample fitting of dU, which undermines the quantitative validation. If the authors can show that dU was fixed independently of the 67 nuclei (or present a leave-one-out analysis), a minor revision might suffice. The parity assumption is also a concern but is secondary. I recommend the editor seek a clear statement on the fitting protocol."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a competent, useful paper—a systematic RDFT-based statistical NLD calculation over 67 even-even nuclei, with microscopically computed spin cut-off and a genuinely instructive 98Sr deformation study. But the central validation numbers are less clean than they look. The damping parameter dU = 6.27 MeV is fitted to experimental data at the neutron separation energy, and then Section IIID reports f_rms values against s- and p-wave resonance spacings at the same separation energy, with no statement that the fit excluded these nuclei. That makes the f_rms figures a measure of fit quality, not out-of-sample prediction. It isn't fatal—dU is one global parameter—but the comparison to other models' f_rms values is not apples-to-apples unless those models also fitted damping parameters to the same data. At minimum the paper should say explicitly which data determined dU and re-label the validation claims accordingly.\n\nThe introduction advertises explicit treatment of both spin and parity dependence, but Eqs. (11) and (12) are J-dependent only, and Section IIID simply assumes positive and negative parities are equally distributed. The p-wave spacings hinge on that assumption. If they want to claim parity treatment, they need to actually compute parity-dependent densities or soften the claim.\n\nWhat deserves credit: the method follows the well-established Demetriou-Goriely statistical framework, the RDFT input is cleanly described, two density functionals are compared, and the 98Sr shape-coexistence analysis gives concrete microscopic insight into how single-particle level density and pairing gaps drive deformation-dependent state density and collective enhancement. The comparison to Oslo data and to MSk7/BSk14 models is informative. The authors also flag the neglect of temperature-dependent single-particle levels in their own summary, which is honest.\n\nOverall this is a solid new application, not a methodological breakthrough, and it doesn't resolve a long-open question. It does provide a useful systematic dataset and a credible reference point for RDFT-based NLD work. With revision addressing the dU fit and the parity overclaim, it earns a peer-reviewed slot. I would send it to review rather than desk reject.","headline":"Solid systematic RDFT-based NLD application, but the f_rms validation is partly in-sample because dU is fit to neutron-resonance data, and the claimed parity dependence is really just parity equipartition.","tokens_in":17802,"tokens_out":2526,"would_cite":true,"duration_ms":31599,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81V35"],"pacs":["21.10.Ma"],"model":"deepseek-v4-flash","headline":"A relativistic-density-functional statistical model reproduces nuclear level densities and neutron resonance spacings from self-consistent single-particle levels alone.","keywords":["nuclear level density","relativistic density functional theory","finite-temperature BCS","neutron resonance spacing","spin cut-off parameter","collective enhancement","shape coexistence","even-even nuclei"],"falsifier":"Measure or calculate parity-resolved level densities near the neutron separation energy for a subset of the 67 nuclei, for example from neutron resonances with assigned parities or from parity-projected shell-model calculations; if the positive-to-negative parity ratio at the relevant spin and energy departs substantially from one, the equipartition assumption fails and the p-wave comparisons become an unreliable test of the model.","tokens_in":16997,"feed_emoji":"⚛️","tokens_out":8648,"duration_ms":92768,"temperature":0.7,"pith_summary":"The paper tries to show that a microscopic statistical model built on relativistic density functional theory can predict nuclear level densities — how many quantum states a nucleus has at a given excitation energy — without fitting to level-density data. It takes self-consistent single-particle levels from RDFT, adds pairing in finite-temperature BCS theory, and multiplies in collective rotational and vibrational enhancement. The spin cut-off parameter that spreads level density over angular momenta is computed microscopically from the same single-particle levels, so shell structure and deformation are preserved. Tested against neutron resonance spacings for 67 even-even nuclei, the model agrees within a factor of about 2.3 (DD-ME2), comparable to established non-relativistic microscopic methods and better than a previous relativistic combinatorial calculation without pairing. If right, this gives a parameter-free-from-fitting route to level densities for nuclei where resonance data do not exist.","feed_headline":"Model matches neutron resonance spacings to factor 2.3","feed_subtitle":"Relativistic mean-field single-particle levels plus pairing and collective boosts reproduce data for 67 even-even nuclei.","key_machinery":"The load-bearing object is the spin cut-off parameter σ², computed by a microscopic sum over the magnetic-substate projections m² of the RDFT single-particle levels, weighted by the thermal occupation factor sech²(E/2T). This replaces fitted empirical spin-cutoff formulas and lets shell closures and deformation imprint themselves on the angular-momentum distribution. Around it sits the finite-temperature BCS saddle-point state density: quasiparticle energies from RDFT levels plus pairing gaps fix entropy and excitation energy, and collective effects enter as a rotational enhancement factor built from the perpendicular moment of inertia, suppressed by pairing, plus a vibrational enhancement f","core_discovery":"Central claim: a statistical level-density model built from relativistic density functional theory can reasonably reproduce measured nuclear level densities and s- and p-wave neutron resonance spacings. The recipe uses self-consistent RDFT single-particle levels, finite-temperature BCS pairing, saddle-point state density, and rotational plus vibrational enhancement; the spin cut-off parameter is computed microscopically from the same levels. For 67 even-even nuclei with s-wave data, f_rms = 2.29 (DD-ME2) or 2.56 (DD-PC1); including 18 p-wave sets gives 2.41/2.64. The paper reads these as comparable to non-relativistic microscopic models, better than a relativistic combinatorial calculation w","pith_inferences":["The parity-equipartition assumption blurs the distinction between spin and parity dependence: the p-wave comparison mostly tests the spin distribution, so a parity-projected version could sharpen the test.","Because the model computes σ² from thermalized single-particle levels, it could be inserted directly into Hauser–Feshbach reaction codes, replacing empirical spin-cutoff systematics and improving extrapolations to short-lived nuclei.","The single-particle levels are temperature-independent in this paper; at excitation energies well above the neutron separation energy, thermal reshaping of the mean field may change the high-energy slope, so the method's reach beyond measured energies is the main open question.","A sharper per-nucleus test than the global f_rms would compare the predicted spin cut-off parameter with values extracted from angular distributions or from shell-model Monte Carlo, since σ² is the channel through which all structure dependence flows."],"forward_implications":["For even-even nuclei without measured resonances, the approach provides level densities with accuracy comparable to existing non-relativistic microscopic models, so it can stand in for them in reaction and astrophysics calculations.","Because the spin cut-off parameter comes out of the structure calculation, the model automatically carries shell-closure dips and deformation effects that empirical spin-cutoff formulas miss or misplace.","The DD-ME2 and DD-PC1 interactions give consistent results, with DD-ME2 slightly closer to resonance data; remaining differences quantify functional uncertainty.","The 98Sr analysis shows that near shape coexistence, level density and rotational enhancement vary with deformation through single-particle level density and pairing gaps, not through collective-phase-space arguments alone.","The quoted agreement assumes equal positive- and negative-parity level densities, so the parity-blind part of the model is what is tested; parity-dependent extensions are the natural next step."],"supporting_citations":[{"why":"supplies the experimental s-wave and p-wave neutron resonance spacings for the 67 nuclei used as the quantitative benchmark.","marker":"[3]"},{"why":"provides the vibrational enhancement expressions and the damping function interpolating between deformed and spherical level densities.","marker":"[35]"},{"why":"sets the non-relativistic HFB-combinatorial comparison with pairing, at f_rms = 2.3, against which the present accuracy is judged.","marker":"[36]"},{"why":"the prior relativistic combinatorial calculation without pairing whose f_rms = 3.62 the present method improves on.","marker":"[39]"},{"why":"supplies the microscopic statistical state-density formalism and the microscopic spin cut-off formula adapted here.","marker":"[46]"},{"why":"earlier RDFT-based microscopic level-density treatment that this work extends to systematic spin- and parity-dependent comparison.","marker":"[50]"},{"why":"QRPA-plus-boson-expansion calculation whose f_rms = 1.65 marks the best microscopic benchmark cited.","marker":"[54]"},{"why":"supplies the finite-range separable pairing interaction used in the RDFT single-particle calculations.","marker":"[86]"},{"why":"provides the RDFT mean-field framework and the DD-ME2 interaction that generate the single-particle input.","marker":"[92]"},{"why":"defines the DD-PC1 density-dependent interaction used for the alternative calculations.","marker":"[94]"}],"fun_headline_variants":["RDFT level density model matches resonance spacings to 2.3","Microscopic level density from relativistic theory reproduces data","Relativistic model for nuclear level density: test on 67 nuclei","Relativistic DFT level density: matches neutron spacings (f=2.3)","67 even-even nuclei: relativistic level density model fits resonance data"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The calculation assumes positive- and negative-parity states are equally numerous at every excitation energy when comparing with neutron resonance data; if real parity distributions are skewed, the quoted p-wave accuracy and the mixed f_rms values would shift.","fun_headline_variants_meta":{"raw":{"variants":["RDFT level density model matches resonance spacings to 2.3","Microscopic level density from relativistic theory reproduces data","Relativistic model for nuclear level density: test on 67 nuclei","Relativistic DFT level density: matches neutron spacings (f=2.3)","67 even-even nuclei: relativistic level density model fits resonance data"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000818,"raw_usage":{"total_tokens":3408,"prompt_tokens":725,"completion_tokens":2683,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":469,"completion_tokens_details":{"reasoning_tokens":2588}},"tokens_in":469,"tokens_out":2683,"duration_ms":21366,"temperature":1.0,"reasoning_tokens":2588,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T15:04:54.240032+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure or calculate parity-resolved level densities near the neutron separation energy for a subset of the 67 nuclei, for example from neutron resonances with assigned parities or from parity-projected shell-model calculations; if the positive-to-negative parity ratio at the relevant spin and energy departs substantially from one, the equipartition assumption fails and the p-wave comparisons become an unreliable test of the model.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the vibrational enhancement expressions and the damping function interpolating between deformed and spherical level densities."},{"cited_title":"Hilaire and S","cited_arxiv_id":null,"evidence_quote":"sets the non-relativistic HFB-combinatorial comparison with pairing, at f_rms = 2.3, against which the present accuracy is judged."},{"cited_title":"Uhrenholt, S","cited_arxiv_id":null,"evidence_quote":"the prior relativistic combinatorial calculation without pairing whose f_rms = 3.62 the present method improves on."},{"cited_title":"Decowski, W","cited_arxiv_id":null,"evidence_quote":"supplies the microscopic statistical state-density formalism and the microscopic spin cut-off formula adapted here."},{"cited_title":"Rahmatinejad, T","cited_arxiv_id":null,"evidence_quote":"earlier RDFT-based microscopic level-density treatment that this work extends to systematic spin- and parity-dependent comparison."},{"cited_title":"Rossignoli, A","cited_arxiv_id":null,"evidence_quote":"QRPA-plus-boson-expansion calculation whose f_rms = 1.65 marks the best microscopic benchmark cited."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the RDFT mean-field framework and the DD-ME2 interaction that generate the single-particle input."},{"cited_title":"Døssing and S","cited_arxiv_id":null,"evidence_quote":"defines the DD-PC1 density-dependent interaction used for the alternative calculations."}],"review_version":1}