{"id":"9d4b9dbe-54ae-46ff-9465-362dfb3fd4f7","arxiv_id":"2411.16014","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"IMSRG(3)-N7 calculations for calcium-44, -48, and -52 show that three-body flow corrections substantially lower the 48Ca 2+ energy, improving the description of the N=28 shell closure, while leaving the 52Ca charge-radius puzzle unresolved.","lead":"This paper applies a more complete many-body method, IMSRG(3)-N7, to neutron-rich calcium isotopes and finds that it fixes part of a long-standing mismatch in the excited states of calcium-48. It also shows that the method does not explain the experimentally large charge radius of calcium-52, so that puzzle must come from somewhere else.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 48Ca 2+ correction is not converged in the three-body model-space truncation, so the claimed ~0.68 MeV improvement and its closeness to EOM-CCSD(T) are not yet quantitatively established.","rationale":"The reader's weakest assumption identified the same two coupled issues: the three-body model-space truncation and the dropping of three-body valence operators. I focus on the model-space truncation as the most load-bearing because the paper itself concedes the 2+ result is not converged, and the specific -0.677 MeV correction is central to the claimed improvement. The concrete test of extending to emax,3b=7 directly determines whether the correction is robust. The cluster-hierarchy approximation is a real secondary concern, but the ground-state consistency between IMSRG(3)-N7 and VS-IMSRG(3)-N7 provides some indirect evidence for its validity, and it is less directly tied to the headline number. The paper is transparent, provides honest caveats, and includes valuable cross-checks (e.g., the IMSRG/VS-IMSRG difference reduction, the spin-orbit radius correction). No internal inconsistency or unsupported factual claim was found. The reader's CONDITIONAL verdict is appropriate: the qualitative improvement is likely, but the quantitative magnitude requires confirmation. Therefore I recommend no change to the verdict, while emphasizing that the convergence check is necessary to upgrade to ACCEPT.","tokens_in":22402,"tokens_out":4751,"duration_ms":44311,"concrete_test":"Compute VS-IMSRG(3)-N7 for 48Ca at emax,3b=7 with E3max up to 21 (or at minimum E3max=20 with emax,3b=6) and compare the resulting 2+ excitation energy to the emax,3b=6, E3max=18 value of 4.253 MeV. If the energy shifts by more than ~0.15 MeV (i.e., >20% of the claimed 0.68 MeV correction) or deviates from the established downward trend, then the headline numerical correction is not converged and the quantitative claim must be revised. Additionally, if feasible, compare against a valence-space diagonalization that retains the three-body valence-space operators for a lighter isotope (e.g., 44Ca) to quantify the cluster-hierarchy truncation error for excited states.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central positive claim rests on the VS-IMSRG(3)-N7 correction lowering the 48Ca 2+ energy by 0.677 MeV at emax,3b=6, E3max=18 (Table I), bringing it from 4.930 MeV to about 4.253 MeV, closer to the experimental 3.832 MeV. However, the paper explicitly states (Sec. III A, Fig. 4) that this result is 'far from fully converged' at the largest three-body truncation. The trend in Fig. 4 is monotonically downward as both emax,3b and E3max increase, with no sign of saturation at emax,3b=6, E3max=18. Therefore the 0.68 MeV correction could be substantially larger (or, less plausibly, smaller) at full convergence, which would change the quantitative comparison with EOM-CCSD(T) and with experiment. The cluster-hierarchy approximation of dropping three-body valence-space operators before the final shell-model diagonalization (Sec. II B, W(s→infinity)=0) is a second, unquantified source of error for excitation energies; the paper's consistency argument between IMSRG(3)-N7 and VS-IMSRG(3)-N7 applies only to ground-state energy and charge radius, not to the 2+ energy. Thus the headline improvement is not yet pinned down: the sign and qualitative direction appear robust, but the magnitude—and hence the strength of the claim—depends on untested convergence and on the dropped three-body valence operators.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper applies the recently developed IMSRG(3)-N7 and VS-IMSRG(3)-N7 methods to the calcium isotopes 44Ca, 48Ca, and 52Ca using the 1.8/2.0 (EM) chiral Hamiltonian. It reports that the VS-IMSRG(3)-N7 corrections to the first 2+ excitation energy of 48Ca are large and bring the prediction closer to experiment and to coupled-cluster results with triples, improving the description of the N=28 shell closure. For ground-state energies, charge radii, and neutron skins, the IMSRG(3)-N7 corrections are small, and the charge-radius differences between 52Ca, 48Ca, and 44Ca remain underpredicted. The paper also presents size-extensivity-based estimates of IMSRG(2) many-body uncertainties that are intended to be applicable beyond calcium.","tokens_in":22683,"tokens_out":4562,"duration_ms":44935,"significance":"If the main 2+ result holds, the paper is significant because it identifies a many-body truncation effect as the origin of a long-standing VS-IMSRG(2) discrepancy in a key closed-shell nucleus and demonstrates that the IMSRG(3)-N7 framework can quantify and reduce such uncertainties. The calculations are parameter-free with respect to the calcium observables, since the Hamiltonian is fixed from prior chiral EFT fits, and the manuscript is unusually transparent about basis and truncation convergence, including explicit statements of what is not converged. The correction to the spin-orbit charge-radius operator in Appendix B is a useful technical contribution. However, the headline 2+ correction is not converged at the largest truncation studied, so the quantitative significance of the improvement, and the closeness to EOM-CCSD(T), are not yet fully established.","major_comments":[{"comment":"The central positive claim rests on a 0.677 MeV lowering of the 48Ca 2+ energy at emax,3b=6, E3max=18 (Table I), but the paper itself states in Sec. III A that the 2+ energy is \"far from fully converged\" and that \"a quantitative assignment of the VS-IMSRG(3)-N7 prediction is not possible.\" The trend in Fig. 4 is monotonically downward with no sign of saturation at the largest truncation, so the converged correction could be substantially larger (or, less plausibly, smaller) than the reported value. As written, the abstract's claim of a \"significantly better description\" and the quantitative comparison with EOM-CCSD(T) and experiment go beyond what the convergence evidence supports. The authors should either provide additional truncation points or a conservative extrapolation of the remaining correction, or clearly reframe the 2+ result as a qualitative/directional improvement and adjust the abstract and conclusion accordingly.","section":"Sec. III A, Fig. 4, Table I"},{"comment":"The VS-IMSRG(3)-N7 calculation drops three-body valence-space operators before the final shell-model diagonalization, motivated by the cluster hierarchy. The authors justify this approximation by the improved consistency between IMSRG(3)-N7 and VS-IMSRG(3)-N7 results for the ground-state energy and charge radius (Figs. 2 and 3), but this consistency check is not presented for excitation energies. Since the central claim of the paper concerns the 2+ excitation energy, the effect of the discarded three-body valence-space operators on Eex(2+) is an unquantified source of error at potentially the same order as the claimed correction. A test of this approximation in a tractable system, or at least an explicit statement of its expected size for excitation energies, is needed before the magnitude of the 2+ correction can be considered established.","section":"Sec. II B, W(s -> infinity) = 0"}],"minor_comments":[{"comment":"The phrase \"significantly better description of the first 2+ excitation energy of 48Ca\" should be hedged to reflect the non-convergence stated in Sec. III A, e.g., \"indicates a significantly better description\" or \"is consistent with a substantially improved description.\"","section":"Abstract and Conclusion"},{"comment":"Table I should include a footnote that the IMSRG(3)-N7 corrections are computed at emax,3b=6, E3max=18 and that the 48Ca 2+ correction is not fully converged in the three-body model-space truncation.","section":"Table I"},{"comment":"The general uncertainty estimates for IMSRG(2) charge radii and neutron skins are derived from corrections that are not fully converged for all observables; the text should state more explicitly that these are estimates based on truncated IMSRG(3)-N7 corrections and may be lower bounds.","section":"Sec. III D"},{"comment":"It would improve clarity to add a visual or textual marker on the emax,3b=6, E3max=18 point noting that this result is not converged, so that the figure cannot be misread as a final prediction.","section":"Fig. 4"},{"comment":"The correction to the spin-orbit radius operator is clearly described, but the authors should consider adding a sentence on the numerical impact of this correction on the reported charge radii, since the manuscript otherwise focuses on small differences.","section":"Appendix B"}],"recommendation":"major_revision","confidential_remarks":"The paper is a good fit for the journal and the authors are admirably transparent about the limitations of their calculations. The main issue is that the abstract and the framing of the 48Ca 2+ result are somewhat stronger than the convergence evidence allows, and the dropped three-body valence-space operators are not directly tested for excitation energies. I do not see grounds for rejection, as the qualitative direction appears robust and the methodological contribution is valuable. If the authors can either add one more truncation point or a conservative extrapolation for the 2+ correction, and soften the quantitative claims in the abstract and Table I, the paper would be acceptable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a solid, honest paper. It delivers the first systematic IMSRG(3)-N7 study of the calcium chain, and its main positive claim—that the 2+ energy of 48Ca comes down substantially when three-body operators are included—is real in direction, though the magnitude is not yet pinned down. The negative result on charge radii is the more robust finding: the IMSRG(3)-N7 corrections to radii are strongly correlated across 44,48,52Ca, so they barely change the underpredicted 52Ca–48Ca difference.\n\nWhat is actually new: the application of the recently developed IMSRG(3)-N7 method to 44,48,52Ca, including a systematic map of the three-body model-space truncation, a demonstration that 3N flow corrections lower excitation energies by a system-dependent factor, and practical uncertainty estimates for IMSRG(2) based on size extensivity. The paper also corrects a published error in the spin-orbit radius operator. No parameters are fitted to calcium observables; the Hamiltonian is fixed from elsewhere, so all reported numbers are genuine predictions. The transparency is a real strength: convergence in emax, emax,3b, and E3max is shown, and the authors explicitly state where results are not converged.\n\nThe soft spots are real but mostly acknowledged. The 48Ca 2+ correction is not converged at the largest truncation; the paper itself says a quantitative assignment is not possible. That means the apparent near-agreement with EOM-CCSD(T) and experiment is suggestive, not established. The second soft spot is the cluster-hierarchy approximation of dropping three-body valence-space operators before the final diagonalization. The consistency argument between IMSRG and VS-IMSRG is indirect and applies to ground-state energy and radius, not directly to excitation energies. The 'common factor' reduction of the spectrum is an empirical observation, not yet explained. The uncertainty estimates are rules of thumb based on a limited set of systems, but they are presented as such.\n\nNone of these are load-bearing flaws. The qualitative conclusions—that IMSRG(3)-N7 moves 48Ca's 2+ down substantially and does not resolve the radius puzzle—are well supported. The paper's claims are appropriately hedged, and the stress-test concern about non-convergence is already visible in the text; it is a limitation, not a hidden defect.\n\nThis paper deserves a serious referee. The ideal referee would ask for an independent check of the dropped three-body valence operators, or a larger three-body space in a tractable case, but the work is technically careful and the community will cite it. Send it to peer review.","headline":"A careful, transparent IMSRG(3)-N7 study of calcium isotopes with a genuine but not yet converged correction to the 48Ca 2+ energy, and a robust negative result on the charge-radius puzzle.","tokens_in":23287,"tokens_out":1866,"would_cite":true,"duration_ms":19391,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["21.60.Cs","27.40.+z","21.10.Ft"],"model":"deepseek-v4-flash","headline":"Including three-body operators in IMSRG calculations lowers the predicted $^{48}$Ca first $2^+$ excitation energy by about 0.68 MeV, bringing it close to experiment and strengthening the description of the $N=28$ shell closure.","keywords":["calcium isotopes","in-medium similarity renormalization group","three-body operators","charge radii","2+ excitation energy","shell closure N=28","many-body uncertainty","ab initio nuclear structure"],"falsifier":"Recompute the $^{48}$Ca first $2^+$ energy retaining the three-body valence-space operators in the final diagonalization, or with $e_{\\max,3b}>6$ (e.g., $e_{\\max,3b}=7$) on a tractable model space; if the converged value moves back above roughly 4.5 MeV, the claimed improvement would be a truncation artifact rather than a many-body correction.","tokens_in":22177,"feed_emoji":"⚛️","tokens_out":13468,"duration_ms":111272,"temperature":0.7,"pith_summary":"This paper tests whether a known weakness of the in-medium similarity renormalization group (IMSRG) at two-body truncation, the overprediction of the first $2^+$ excitation energy of $^{48}$Ca, is cured by including three-body operators in the flow. In the IMSRG(3)-$N^7$ calculation the $2^+$ energy drops by about 0.68 MeV at the largest three-body model space used, moving from 4.930 MeV toward the measured 3.832 MeV and toward coupled-cluster results with triples. The authors then ask whether the same upgrade fixes the underpredicted charge radius of $^{52}$Ca relative to $^{48}$Ca, and find that it does not: the corrections are nearly identical in the two isotopes and largely cancel in the difference. The same calculations provide size-extensivity-based estimates that IMSRG(2) carries roughly 2-3% uncertainty on correlation energies, 1-1.5% on charge radii, and 5-7.5% on neutron skins.","feed_headline":"Three-body operators pull 48Ca's 2+ energy close to experiment","feed_subtitle":"IMSRG(3)-N7 cuts the overprediction by ~0.7 MeV; the 52Ca-48Ca radius gap stays unexplained.","key_machinery":"The central object is the IMSRG(3)-$N^7$: an in-medium similarity renormalization group evolution that includes normal-ordered three-body operators $W(s)$ throughout the flow, truncated to all terms that scale no worse than $N^7$ in basis size. After decoupling the $^{40}$Ca core and the $0\\hbar\\omega$ neutron valence space, the residual three-body valence-space operators are set to zero following the cluster hierarchy (one-body effects dominate two-body, which dominate three-body), and a standard shell-model diagonalization is performed. The mechanism carrying the argument is that induced three-body operators feed back into the effective one- and two-body valence-space interactions, lowering the computed excitation spectrum; because the IMSRG is size extensive, meaning errors scale with the system rather than growing disproportionately, the observed correction sizes can be quoted as transferable uncertainty estimates for IMSRG(2) in other medium-mass systems.","core_discovery":"The paper's central finding is that the valence-space IMSRG(3)-$N^7$, which keeps normal-ordered three-body operators while solving the renormalization-group flow and then omits them in the final shell-model diagonalization, changes calcium spectra substantially while leaving isotope shifts of charge radii almost untouched. For $^{48}$Ca the first $2^+$ state falls from 4.930 MeV in VS-IMSRG(2) by $-0.68$ MeV at $e_{\\max,3b}=6$, $E_{3\\max}=18$, reducing the gap to experiment from about 1.1 MeV to about 0.42 MeV and matching the direction and size of triple-excitation corrections in coupled-cluster theory. All predicted levels in $^{44}$Ca, $^{48}$Ca, and $^{52}$Ca are shifted downward by a roughly common factor that is largest for $^{48}$Ca. For charge radii, the IMSRG(3)-$N^7$ corrections to $^{44}$Ca, $^{48}$Ca, and $^{52}$Ca are strongly correlated, so the differences $R_{\\rm ch}(^{52}{\\rm Ca})-R_{\\rm ch}(^{48}{\\rm Ca})$ and $R_{\\rm ch}(^{44}{\\rm Ca})-R_{\\rm ch}(^{48}{\\rm Ca})$ barely move and the large measured $^{52}$Ca-$^{48}$Ca radius difference remains underpredicted. The paper reads these results as showing that the $2^+$ discrepancy was largely a many-body truncation artifact, while the radius puzzle is not resolved at this order.","pith_inferences":["A testable extension: apply the same VS-IMSRG(3)-$N^7$ treatment to other closed-shell nuclei with overpredicted $2^+$ states, such as $^{78}$Ni, to see whether the calcium-sized downward shift generalizes.","Because the radius corrections cancel almost completely in isotope differences, the persistent $^{52}$Ca-$^{48}$Ca puzzle likely sits in the Hamiltonian or in physics outside this truncation, such as multishell valence-space excitations, rather than in the IMSRG(2) approximation.","One could directly test the cluster-hierarchy assumption by keeping the three-body valence-space operators in the final diagonalization with a three-body-capable shell-model solver; the paper predicts their effect is small, so a noticeable shift would mark a boundary of the approximation."],"forward_implications":["VS-IMSRG(2) $2^+$ energies in closed-shell nuclei such as $^{48}$Ca can be off by roughly 0.7 MeV from the higher-order result, so comparing such states to experiment requires the three-body truncation or a comparable uncertainty estimate.","The $^{48}$Ca $2^+$ value is not fully converged at $e_{\\max,3b}=6$, $E_{3\\max}=18$, so the final prediction will move somewhat as the three-body model space grows, but the correction has a definite sign and approximate size.","Charge-radius isotope shifts in the calcium chain change by less than 10% when going from IMSRG(2) to IMSRG(3)-$N^7$, so the underprediction of the $^{52}$Ca-$^{48}$Ca difference is not a many-body truncation effect at this order.","For soft chiral Hamiltonians, IMSRG(2) results should be assigned uncertainties of roughly 2-3% on correlation energies, 1-1.5% on charge radii, and 5-7.5% on neutron skins."],"supporting_citations":[{"why":"Introduced the IMSRG formalism with three-body operators that this work pushes to the calcium isotopes.","marker":"[36]"},{"why":"Defined the IMSRG(3)-N7 truncation with flowing three-body operators and its relation to other truncations.","marker":"[37]"},{"why":"Coupled-cluster study of 48Ca and 78Ni showing that triples lower the 2+ energy; it is the benchmark this result matches.","marker":"[25]"},{"why":"Provided the valence-space IMSRG framework and ensemble normal ordering used for the shell-model diagonalization.","marker":"[27]"},{"why":"Supplied the cluster-hierarchy argument for dropping three-body valence-space operators before the final diagonalization.","marker":"[35]"},{"why":"The 1.8/2.0 (EM) chiral Hamiltonian used for the main calculations.","marker":"[56]"},{"why":"Experimental measurement establishing the large 52Ca-48Ca charge radius difference that the paper seeks to understand.","marker":"[5]"},{"why":"Coupled-cluster benchmark for 48Ca charge radius and neutron skin used in the comparison.","marker":"[26]"},{"why":"Experimental excitation energies, including the 48Ca 2+ value of 3.832 MeV.","marker":"[69]"}],"fun_headline_variants":["New three-body terms fix 48Ca's 2+ shift","Three-body operators close 48Ca excitation gap","IMSRG(3) explains 48Ca's 2+ but not radius","Calcium twist: spin gap fixed, radius puzzle remains","Three-body effects fix calcium 48, miss radius gap"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The headline correction assumes that the three-body model-space truncation ($e_{\\max,3b}\\le 6$, $E_{3\\max}\\le 18$) and the dropping of three-body valence-space operators before diagonalization are sufficient to fix the sign and size of the $2^+$ shift, even though that shift is not fully converged.","fun_headline_variants_meta":{"raw":{"variants":["New three-body terms fix 48Ca's 2+ shift","Three-body operators close 48Ca excitation gap","IMSRG(3) explains 48Ca's 2+ but not radius","Calcium twist: spin gap fixed, radius puzzle remains","Three-body effects fix calcium 48, miss radius gap"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000681,"raw_usage":{"total_tokens":3190,"prompt_tokens":1140,"completion_tokens":2050,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":756,"completion_tokens_details":{"reasoning_tokens":1963}},"tokens_in":756,"tokens_out":2050,"duration_ms":13959,"temperature":1.0,"reasoning_tokens":1963,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:37:35.820277+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the $^{48}$Ca first $2^+$ energy retaining the three-body valence-space operators in the final diagonalization, or with $e_{\\max,3b}>6$ (e.g., $e_{\\max,3b}=7$) on a tractable model space; if the converged value moves back above roughly 4.5 MeV, the claimed improvement would be a truncation artifact rather than a many-body correction.","supporting_citations":[{"cited_title":"A nucleus-dependent valence-space approach to nuclear structure","cited_arxiv_id":"1607.03229","evidence_quote":"Provided the valence-space IMSRG framework and ensemble normal ordering used for the shell-model diagonalization."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Experimental excitation energies, including the 48Ca 2+ value of 3.832 MeV."}],"review_version":1}