{"id":"188a42e4-af9f-4612-8d3d-889682161aaa","arxiv_id":"2505.01072","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A shell model driven by Gogny-D1S density-functional matrix elements reproduces low-lying spectra of sd-shell and Ca isotopes with accuracy comparable to empirical interactions, with known systematic deviations.","lead":"This paper tests a hybrid method that builds nuclear shell-model interactions from a Gogny density functional without fitting each nucleus, and applies it to oxygen, neon, magnesium, and calcium isotopes. The calculated spectra and ground-state energies match experiment reasonably well in most cases, while revealing a systematic underestimate near neutron-rich and closed-shell nuclei.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The broad 'comparable accuracy' claim is undercut by the paper's own O and Ca spectra, where 2+ excitation energies are systematically low; this is attributed to 0ℏω truncation but not quantitatively demonstrated.","rationale":"The reader's weakest-assumption analysis identifies the 0ℏω valence-space truncation and the corresponding omission of core excitations as the key weak point. My reading converges on the same structural issue from the direction of the paper's own displayed results: Figures 5 and 8 show systematically low excitation energies in O and Ca isotopes, and Figure 3 shows growing ground-state deviations in neutron-rich Ca. The paper explicitly acknowledges this in Section III and the Summary, so the weakness is not hidden, but the Summary's 'accuracy comparable ... in O, Ne, Mg, and Ca isotopes' is stated without the quantitative caveat that the figures demand. The strongest claim as phrased in the reader's report includes a hedge, 'in the regions where the 0 hbar-omega valence space is adequate,' but the paper does not define that region. This makes the central claim hard to test and potentially unfalsifiable as written. I do not see a more damaging technical flaw: the TBME comparison with USDB and SDPF-MU is reasonable, the iterative density treatment follows the established approach of Ref. [12], and the B(E2) effective-charge concern is secondary because the paper itself limits the B(E2) discussion to Mg isotopes and shows sensitivity to the effective-charge choice. The reproducibility issues and the unstated core-energy input are real but addressable. A conditional verdict is appropriate, with the condition that the claim be restricted to the mass regions where the 0ℏω truncation is adequate, or that the paper provide a quantitative estimate of the core-polarization correction.","tokens_in":10051,"tokens_out":8376,"duration_ms":99206,"concrete_test":"Compute the root-mean-square deviation between calculated and experimental low-lying excitation energies (for example, 2+_1, 4+_1, and yrast levels up to about 5 MeV) separately for the O, Ne, Mg, and Ca chains, using both the D1S hybrid interaction and the appropriate empirical interaction (USDB for sd-shell nuclei, SDPF-MU for Ca isotopes). If the D1S rms exceeds the empirical-interaction rms by more than roughly 50% for the O or Ca chains while Ne and Mg are comparable, the 'comparable accuracy' claim should be restricted to the latter or made conditional on including core-excitation effects. As a second, more expensive check, repeat the 18-22O and 42-50Ca calculations in a space that includes 2ℏω cross-shell configurations; if the compressed 2+ energies move upward toward experiment, the inert-core truncation is confirmed as the immediate cause.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires the D1S-based hybrid interaction to match the empirical USDB and SDPF-MU interactions in reproducing experimental spectra. The paper's own Figures 5 and 8 show this is not the case in two of the four isotope chains: in 18O, 20O, and 22O the calculated 2+ excitation energies are sizably below both experiment and USDB, and the Ca isotopes are systematically compressed, including 48Ca whose large 2+ energy is underestimated. The text attributes these deficits to the inert 16O and 40Ca cores and the 0ℏω valence-space truncation. That attribution is plausible, but it means the summary statement that the model reproduces spectra 'with an accuracy comparable to the existing empirical interactions in O, Ne, Mg, and Ca isotopes' is either overbroad or circularly restricted to nuclei where the truncation happens to be acceptable. No quantitative criterion is given for where the 0ℏω space is adequate, so the central claim cannot be evaluated as stated. The neutron-rich Ca ground-state energies also deviate increasingly with neutron number (Figure 3), and the paper itself lists core excitations as a possible source. The load-bearing weakness is therefore not the TBME construction, which appears internally consistent, but the unsupported breadth of the headline claim relative to the known model-space truncation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports shell-model calculations with a Hamiltonian derived from the Gogny-D1S density functional. The single-particle energies and two-body matrix elements are obtained from the Gogny interaction in a 0ℏω harmonic-oscillator basis, with the density in the density-dependent term determined iteratively from shell-model ground states. Ground-state energies, low-lying spectra, and B(E2) values for O, Ne, Mg, and Ca isotopes are compared with experiment and with the empirical USDB and SDPF-MU interactions. The authors find good agreement for Ne and Mg, systematic underestimation of excitation energies in O and Ca, and reasonable but underestimated B(E2) values in Mg. The paper concludes that the hybrid model reproduces spectra with accuracy comparable to empirical interactions across all four isotope chains, and suggests future extensions to heavier nuclei.","tokens_in":10324,"tokens_out":5942,"duration_ms":56827,"significance":"The paper proposes a hybrid shell-model density-functional approach in which the two-body matrix elements of a Gogny-D1S EDF are evaluated in a 0ℏω valence space with a self-consistently determined density. The method is a step toward non-empirical shell-model interactions applicable across a broad mass range. Strengths include the use of a well-established EDF, the reproduction of the N=16 subshell closure in 24O, the overall good agreement for Ne and Mg spectra, and the careful treatment of center-of-mass corrections. However, the headline claim of accuracy comparable to empirical interactions is not supported for O and Ca isotopes, and the B(E2) comparison relies on effective charges fitted to a different interaction. With a revised claim and clarified methodology, the paper would be a useful contribution.","major_comments":[{"comment":"The summary claim that the model reproduces energy spectra \"with an accuracy comparable to the existing empirical interactions in O, Ne, Mg, and Ca isotopes\" is not supported by the results shown. In 18O, 20O, and 22O the calculated 2+ energies are sizably below both experiment and USDB (Fig. 5), and in the Ca isotopes the 2+ energies are systematically compressed, including 48Ca where a large 2+ energy is a key feature (Fig. 8). The paper attributes these deficiencies to the inert-core/0ℏω truncation, but does not quantify where the 0ℏω space is adequate. Please either narrow the claim to the nuclei where the agreement holds (notably Ne and Mg), or provide a quantitative measure (e.g., RMS deviations relative to experiment and to USDB/SDPF-MU) and a demonstration that core excitations are the cause.","section":"§IV and Figs. 5, 8"},{"comment":"The B(E2) comparison uses effective charges (ep,en)=(1.36,0.45)e that were fitted to the USDB interaction, not to the D1S-based Hamiltonian. Judging the D1S wavefunctions with USDB-fitted charges is not a fair test of the hybrid model and is likely responsible for part of the systematic underestimation. Please fit effective charges to the D1S interaction or present the results as a qualitative sensitivity study; in either case, the statement that the B(E2) values show \"reasonable agreement\" needs to be justified quantitatively.","section":"§III, Fig. 9"},{"comment":"The iterative construction of ρ(r) is incompletely specified. It is not stated whether the density ρ used in the density-dependent term includes the core density plus the valence-shell density, or only the valence density. Since the initial Woods-Saxon density presumably describes the full nucleus, while the \"ground-state density given by the shell model wave function\" contains only valence nucleons, the procedure is ambiguous. Please specify the total density construction and state the convergence criterion for the iteration.","section":"§II, Eq. (3) and following paragraph"},{"comment":"The core energy E_Core is introduced but never defined. It should be clarified how E_Core is computed (e.g., from a spherical Hartree-Fock calculation with the same interaction) and whether it is consistent with the single-particle energies TSPE in Eq. (2). The empirical Coulomb formula of Ref. [17] should be written out explicitly, since the ground-state-energy comparison in Figs. 2 and 3 depends on these choices.","section":"§II, Eq. (4)"}],"minor_comments":[{"comment":"In Eq. (2) and the following line, \"Valance\" should be \"Valence.\"","section":"§II"},{"comment":"The claim that the USDB mass-dependence factor (A/18)^{-0.3} \"is well described by the Gogny ones\" is not supported by any quantitative analysis in Fig. 1; consider fitting the A-dependence or removing the claim.","section":"§III"},{"comment":"The caption of Fig. 1 states that the Ca calculations adopt only the pf shell; this important model-space statement should appear in the main text as well.","section":"Fig. 1 caption"},{"comment":"The sentence \"The points are somewhat scattered but close to the diagonal line, and the Gogny-D1S TBMEs agree with the empirical ones roughly within 2 MeV\" would be more informative with an RMS deviation or correlation coefficient for the TBME comparison.","section":"§III"},{"comment":"The discussion of the center-of-mass correction states that the two-body correction vanishes in the 0ℏω model space; a brief justification (e.g., the momentum operator changes the oscillator quantum number by one) would aid the reader.","section":"§III, Eq. (5)"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the journal's scope and the hybrid approach is worth publishing after revision. The main issue is the overbroad summary claim, which the authors' own figures contradict; the revision should rephrase the claim and ideally add a quantitative accuracy metric. The density-dependence specification in Sec. II also needs clarification, as it is central to the method. I see no fundamental correctness error in the TBME construction itself."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe one thing to know: this is a solid but incremental application of the hybrid scheme from Jiang et al. (Ref [12]) to sd-shell nuclei and Ca isotopes. The method is unchanged; what's new is the numerical territory. The paper is worth a look because it does the comparisons honestly, but the summary sentence claiming accuracy comparable to empirical interactions in all four isotope chains overstates what Figures 5 and 8 actually show.\n\nWhat's good: the TBME scatter plots against USDB/SDPF-MU give a quick read on how close the D1S-derived matrix elements are. The ground-state energies are compared to HFB and AME masses, and the center-of-mass correction discussion is clear. The dripline prediction for 24O is a nice qualitative check. For Ne and Mg, the spectra do land close to USDB and experiment. The paper states its own limitations explicitly—core excitations, empirical ℏω, Coulomb formula—so you never feel misled while reading.\n\nThe soft spots: the load-bearing claim is too broad. In 18-22O and in the Ca chains, the 2+ states come out sizably below both experiment and USDB/SDPF-MU, and the paper says this is likely due to the inert 16O/40Ca cores. That may be right, but no quantitative estimate of the truncation error is given, and no criterion is offered for where 0ℏω is adequate. As written, the phrase 'comparable accuracy in O, Ne, Mg, and Ca isotopes' is only true if you restrict it to the well-behaved cases. Also, the B(E2) analysis uses effective charges fitted to USDB, so the D1S results are not a clean test; the paper acknowledges this, but it means the E2 section is weaker than the spectra section. Reproducibility is so-so: no input TBMEs or convergence thresholds are provided, though KSHELL is open.\n\nNothing here is fatal. The TBME construction is internally consistent, the empirical comparisons are transparent, and the limitations are stated. It's a legitimate extension of an existing program, not a breakthrough. Who for: people working on non-empirical shell-model interactions and Gogny-EDF practitioners. It deserves peer review because the method is promising and the application to Ca is useful, but the referee should require the authors to either narrow the conclusion or quantify the model-space truncation effect, and to add enough input details for reproduction. My recommendation: send it to review, with a request for modest revision.","headline":"A careful but incremental extension of the Gogny-to-shell-model hybrid to sd-shell and Ca nuclei, with an accuracy claim that outruns the paper's own O and Ca spectra.","tokens_in":10854,"tokens_out":3718,"would_cite":true,"duration_ms":34587,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81V35"],"pacs":[],"model":"deepseek-v4-flash","headline":"A shell-model Hamiltonian built from the Gogny D1S density functional, with no fitting to shell-model spectra, reproduces experimental ground-state energies and low-lying spectra of O, Ne, Mg, and Ca isotopes about as well as the…","keywords":["nuclear density functional theory","shell-model calculation","Gogny interaction","sd-shell nuclei","calcium isotopes","effective interaction","beyond-mean-field correlations","electromagnetic transition"],"falsifier":"Repeat the calculation for 18O, 20O, 22O, and the Ca isotopes with a model space enlarged to include core-excited configurations, allowing nucleons to leave the 16O or 40Ca cores. If the 2+ excitation energies remain systematically below experiment even with core excitations included, then the frozen-core assumption is not the cause, and the claim of accuracy comparable to empirical interactions would be contradicted in those nuclei.","tokens_in":9811,"feed_emoji":"⚛️","tokens_out":7165,"duration_ms":69690,"temperature":0.7,"pith_summary":"This paper attempts to show that the shell model can be powered by an energy density functional instead of an interaction fitted to nuclear spectra. The authors build the shell-model Hamiltonian from the Gogny D1S functional, computing the two-body matrix elements with harmonic-oscillator wave functions and iterating the density self-consistently from the shell-model ground state. They compare with the empirical USDB interaction in the sd shell and SDPF-MU in the pf shell, finding comparable accuracy for ground-state energies and low-lying spectra of even-even O, Ne, Mg, and Ca isotopes. If this stands, it would make shell-model calculations available for mass regions where empirical interactions have not been fit. The paper also identifies a clear limitation: excitation energies near the closed 16O and 40Ca cores are systematically underestimated, which it attributes to neglected core excitations.","feed_headline":"Shell model from Gogny functional matches empirical interactions","feed_subtitle":"Interactions built from a density functional, not fitted to data, reproduce spectra of four isotope chains.","key_machinery":"The load-bearing object is the shell-model Hamiltonian built from the Gogny D1S finite-range density-dependent interaction: single-particle energies come from the one-body terms normal-ordered with respect to the inert core, and the two-body matrix elements are computed with harmonic-oscillator basis states at a frequency set by an empirical mass-dependent formula. The density entering the density-dependent term is updated iteratively from the shell-model ground-state wave function until convergence, mimicking self-consistent mean-field practice. This is what turns a mean-field energy functional into a configuration-interaction Hamiltonian without empirical shell-model fitting.","core_discovery":"The central claim is that the Gogny D1S two-body matrix elements, evaluated with harmonic-oscillator orbitals and a self-consistently iterated ground-state density, produce shell-model results for even-even O, Ne, Mg, and Ca isotopes whose low-lying spectra and ground-state energies match experiment with accuracy comparable to the empirically fitted USDB interaction in the sd shell and SDPF-MU in the pf shell. The hybrid calculation correctly reproduces the N=16 subshell gap in 24O, the drop of the 2_2+ state in 26Mg, and the doubly magic character of 48Ca, while systematically underestimating excitation energies in 18-22O and the Ca isotopes, an effect the authors attribute to omitted core excitations.","pith_inferences":["If the density-functional-derived two-body matrix elements remain stable when the valence space is extended, the hybrid approach could be pushed into medium-mass and deformed nuclei, where empirical fits are harder to obtain; that is a natural next test not performed here.","Because the B(E2) values in Mg are systematically low with the USDB-tuned effective charges, a fairer comparison might renormalize charges for the Gogny interaction or test whether transition strengths improve once core-excited configurations are added.","The mass dependence of the USDB interaction is already mimicked by the Gogny functional, so the hybrid scheme may carry an implicit mass dependence; comparing isotopes far from stability would test whether that dependence is quantitatively right."],"forward_implications":["A single density functional can supply useful shell-model Hamiltonians for more than one valence space, here the sd shell and the pf shell, removing the need to fit single-particle energies and two-body matrix elements separately for each space.","The model reproduces the neutron dripline at 24O, a case where the Gogny mean-field calculation fails, showing that the beyond-mean-field correlations of the shell model add essential physics.","The N=16 subshell closure in 24O and the doubly magic character of 48Ca emerge from the interaction without being fitted, supporting the claim that the functional contains the relevant shell structure.","The systematic underestimation of excitation energies in 18-22O and the Ca isotopes indicates that enlarging the valence space to include core excitations would be the next step, as the paper states."],"supporting_citations":[{"why":"Establishes the hybrid scheme of deriving shell-model two-body matrix elements from a Gogny interaction with the density self-consistently obtained from the shell-model wave function.","marker":"[12]"},{"why":"Provides the D1S parameter set used to construct the interaction in this work.","marker":"[14]"},{"why":"Supplies the empirical sd-shell Hamiltonian whose spectra and B(E2) values serve as the accuracy benchmark for O, Ne, and Mg isotopes.","marker":"[5]"},{"why":"Supplies the empirical pf-shell Hamiltonian used as the benchmark for Ca isotopes.","marker":"[6]"},{"why":"Performs the large-scale shell-model diagonalizations needed for the calculations.","marker":"[15]"},{"why":"Provides the experimental ground-state energies against which the model is compared.","marker":"[16]"},{"why":"Supplies the measured B(E2) values used to test transition strengths in the Mg isotopes.","marker":"[32]"}],"fun_headline_variants":["Hybrid shell model from Gogny DFT matches experiment","Gogny-functional shell model rivals empirical fits","Density functional shell model reproduces isotope spectra","Shell model from Gogny DFT matches data without empirical fits","Unfitted Gogny shell model captures nuclear spectra"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that the inert 16O and 40Ca cores stay completely frozen, so only the valence orbits participate; if those cores need to be excited, the claimed accuracy near closed-shell and neutron-rich nuclei would break down.","fun_headline_variants_meta":{"raw":{"variants":["Hybrid shell model from Gogny DFT matches experiment","Gogny-functional shell model rivals empirical fits","Density functional shell model reproduces isotope spectra","Shell model from Gogny DFT matches data without empirical fits","Unfitted Gogny shell model captures nuclear spectra"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000542,"raw_usage":{"total_tokens":2536,"prompt_tokens":821,"completion_tokens":1715,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":437,"completion_tokens_details":{"reasoning_tokens":1640}},"tokens_in":437,"tokens_out":1715,"duration_ms":14752,"temperature":1.0,"reasoning_tokens":1640,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:27:23.628550+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the calculation for 18O, 20O, 22O, and the Ca isotopes with a model space enlarged to include core-excited configurations, allowing nucleons to leave the 16O or 40Ca cores. If the 2+ excitation energies remain systematically below experiment even with core excitations included, then the frozen-core assumption is not the cause, and the claim of accuracy comparable to empirical interactions would be contradicted in those nuclei.","supporting_citations":[{"cited_title":"Gogny-force-derived effective shell-model Hamiltonian","cited_arxiv_id":null,"evidence_quote":"Establishes the hybrid scheme of deriving shell-model two-body matrix elements from a Gogny interaction with the density self-consistently obtained from the shell-model wave function."},{"cited_title":"Tables of E2 transition probabilities from the first 2+ states in even–even nuclei","cited_arxiv_id":null,"evidence_quote":"Supplies the measured B(E2) values used to test transition strengths in the Mg isotopes."}],"review_version":1}