{"id":"41456793-ca66-4bd3-abd9-ec6ce826c1a4","arxiv_id":"2502.04999","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"IBM-CM with configuration mixing reproduces the charge radius dip at 48Ca and connects the isotopic shift to E0 transition strengths in Ca, Ar, and Ti isotopes.","lead":"This paper uses the interacting boson model with two mixed configurations to explain why calcium nuclei shrink at 48Ca and grow sharply at 50Ca. It connects this size change to electric monopole transitions between 0+ states, linking nuclear radii to the location of excited 0+ states.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The E0–radius correlation stands or falls with the κ_int = nb/(nb+2) κ_reg scaling: without it, κ terms enter the E0 matrix element and the η values extracted from radii no longer determine ρ²(E0).","rationale":"The central claim is that a direct correlation between nuclear size and E0 transitions is established. In the IBM-CM, this claim is operationalized by using the same η_reg and η_int in the radius (Eq. 5) and E0 (Eq. 6) operators. The only formal link that lets the E0 expression omit the κ terms is the scaling relation before Eq. (5). The reader identified this as the weakest assumption; I agree and sharpen it: without the scaling, the κ operator is not proportional to the identity over the mixed space, so it contributes to the off-diagonal matrix element. The paper provides no microscopic justification beyond citing ref. [33] and the requirement of linearity, which is a phenomenological convenience rather than a derivation. This does not invalidate the interesting reproduction of the 48Ca dip and 50Ca kink, which appears robust to the choice of radius parametrization because it is driven by the wave-function mixing, but it does mean the quantitative E0 predictions are contingent on the scaling. A targeted refit with κ_int free would directly test whether the central correlation survives. The paper's own caveats (38,42Ca failures) and the unpublished IBM-1 code reinforce the need for a conditional rather than unconditional acceptance. No issue with the authors' integrity; this is a standard phenomenological assumption that needs a sensitivity test.","tokens_in":14091,"tokens_out":7903,"duration_ms":83624,"concrete_test":"Re-fit the Ca radius data without imposing Eq. (5)'s scaling: treat κ_int as a free parameter (or set κ_int = κ_reg or to a value estimated microscopically for the 2p–2h intruder space), refit κ_reg, η_reg, η_int to the 42–52Ca isotopic shifts, and recompute ρ²(E0;0+2→0+1) using the full E0 operator including κ terms, for 44,48Ca. If the predicted ρ²(E0) moves outside the experimental error bars, or if the refit changes the η values by more than the spread needed to maintain the 44,48Ca agreement, the scaling assumption is load-bearing and the correlation claim is conditional on it.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section II, Eq. (5), imposes κ_int = nb/(nb+2) κ_reg [33] so that the κ-part of the radius operator is a configuration-independent constant (κ_reg nb). This makes the ground-state radius mixing-dependent only through the η n_d terms, and it is also what removes the κ contribution from the E0 matrix element in Eq. (6): under the scaling, κ_reg nb P_reg + κ_int(nb+2) P_int = κ_reg nb (P_reg+P_int), which is proportional to the identity and therefore has zero matrix element between orthogonal 0+ states. Conversely, if the intruder configuration does not obey this exact scaling, the κ part contributes to ρ(E0;0+2→0+1) and the fitted η_reg, η_int values from isotopic shifts do not, by themselves, determine the E0 strength. The paper cites ref. [33] for the relation but gives no microscopic derivation or empirical test that the 2p–2h intruder bosons follow the same linear scaling as the regular ones. Since the claimed agreement for 44,48Ca depends on cancellations between the η terms, an unjustified deviation from the scaling would change both the fitted η values and the predicted ρ²(E0). This is the load-bearing assumption for the paper's central correlation statement.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses the interacting boson model with configuration mixing (IBM-CM) to study the even-even Ca, Ar, and Ti isotopic chains. Hamiltonian parameters are fitted to low-energy spectra, and the resulting wave functions are combined with the radius operator of Eq. (4) to compute charge-radius isotopic shifts, including the scaling relation κ_int = nb/(nb+2) κ_reg between regular and intruder configurations. The same η parameters are then used in Eq. (6) to compute ρ²(E0; 0+2 → 0+1) for the Ca isotopes. The calculation reproduces the isotopic-shift dip at 48Ca and the rise at 50Ca, and gives E0 strengths in agreement with data for 44Ca and 48Ca. The paper also reports NSM energy comparisons and discusses the role of binding energy and shell gaps in the Ar and Ti chains.","tokens_in":14585,"tokens_out":9968,"duration_ms":112633,"significance":"If the proposed correlation could be rigorously established, it would provide a simple phenomenological bridge between ground-state charge radii and 0+2 → 0+1 electric monopole strengths in shape-coexisting nuclei, which is of practical interest for interpreting isomerism and configuration mixing. The paper is transparent about its caveats, and the E0 comparison for 44,48Ca is genuinely nontrivial because those E0 data are not used in the fits. The main weaknesses are that the radius comparison is partly circular, since the same isotopic-shift data determine the fitted parameters, and that the E0 formula depends on an unverified scaling relation for the intruder bosons. The Ar and Ti results are presented as exploratory and are not backed by E0 data. Overall, the manuscript is a reasonable phenomenological study, but the central claim needs to be supported by additional justification and quantitative assessment.","major_comments":[{"comment":"The cancellation of the κ contribution to ρ(E0; 0+2 → 0+1) rests on the exact relation κ_int = nb/(nb+2) κ_reg introduced before Eq. (5). Under this relation, κ_reg nb P_reg + κ_int(nb+2) P_int is proportional to the identity, so the κ term has zero matrix element between orthogonal 0+ states. The paper cites ref. [33] for this scaling but gives no microscopic derivation or empirical test for the 2p-2h intruder bosons. If the scaling is violated, a term proportional to (κ_reg nb − κ_int(nb+2)) times the regular-space overlap enters Eq. (6), and the η values fitted to isotopic shifts no longer determine the E0 strength by themselves. Because the claimed radius–E0 correlation is the central result, the authors should either justify the scaling specifically for the intruder configuration or present a sensitivity analysis with respect to deviations from it.","section":"Sec. II, Eqs. (4)-(6)"},{"comment":"The agreement in the radius plots is not an independent test of the model. For each chain, κ_reg, η_reg, and η_int are fitted to the isotopic-shift data displayed in those figures, so the comparison largely reflects the fitting procedure. The paper is honest about the fitting, but the abstract and conclusions should not present the radius agreement as independent evidence for the correlation without distinguishing it from the non-fitted E0 results in Fig. 5. To make the assessment quantitative, the authors should add a residual or rms deviation and state explicitly that only the E0 strengths in 44,48Ca, which were not used in any fit, provide independent support.","section":"Sec. III A, Fig. 4(a); Sec. III B, Fig. 7"},{"comment":"The successful E0 points are only two (44Ca and 48Ca), and the model fails for 42Ca in both radius and E0, for 38Ca and 42Ca in radius, and for 44Ti in radius. These failures are not marginal but are attributed to structural reasons, such as the 'limited number of active bosons' in Sec. III A. The claim that a 'direct correlation' is established is therefore too strong. The authors should either give a systematic condition for when the correlation is expected to hold or reduce the claim to one of a correlation suggested by the 44,48Ca data and the IBM-CM mechanism.","section":"Sec. III A, Fig. 5; Sec. III B"}],"minor_comments":[{"comment":"Equation (5) has a factor 2 multiplying κ_reg nb, whereas Eq. (4) has κ_reg nb. The origin of this factor, presumably the number of nucleons per boson pair, should be stated explicitly so the two equations do not appear inconsistent.","section":"Eqs. (4) and (5)"},{"comment":"The sentence 'The calculated results are in an overall agreement with the data for the long chain of Ca isotopes' is difficult to reconcile with the immediately following admission that the isotopic shifts of 38Ca and 42Ca are not explained; please qualify the claim.","section":"Sec. III A"},{"comment":"The text refers to the experimental 2+1 level of 52Ca at 2.563 MeV while also saying that experimental data for 52Ca are quite limited; please provide the source or clarify that this value is taken from ENSDF.","section":"Sec. III A"},{"comment":"Several levels in Table II are marked with asterisks for unconfirmed spin-parity assignments, but the text does not discuss how this affects the 50Ca comparison; a brief note would help the reader.","section":"Table II"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope for a nuclear-structure journal and contains useful IBM-CM results, but the central radius–E0 correlation rests on the κ_int scaling assumption, which is neither derived nor tested. If the authors can justify that scaling from the cited literature or quantify its sensitivity, I would be willing to accept after a further round. The paper's transparency about its failures is a positive feature."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version: this is the first IBM-CM calculation of the Ca charge-radius puzzle, and it gets the qualitative structure right—the 48Ca dip, the 50Ca kink—by letting intruder mixing drive the size change. It also uses the same η parameters fitted to radii to compute ρ²(E0), and the E0 data for 44,48Ca come out in reasonable agreement. That partial independence is the strongest point: the E0 values are not fitted, they follow from the radius fit and the wave functions.\n\nWhat's actually new is the application, not the framework. The IBM-CM machinery comes from the authors' earlier spectral study of Ca isotopes, and the radius operator is standard. But connecting the size puzzle to E0 strengths in this way is a genuine step, and the paper is honest about where it fails: 38,42Ca and 44Ti are not explained, and the 42Ca E0 deviates from data.\n\nThe soft spot that matters is the κ_int = nb/(nb+2) κ_reg scaling in Eq. (5), cited to ref [33]. This is what makes the κ part of the radius operator configuration-independent and what removes it from the E0 matrix element in Eq. (6). The η's fitted to isotopic shifts then determine ρ²(E0). If the intruder 2p–2h configuration doesn't obey that exact linear scaling, the κ terms re-enter the E0 matrix element and the fitted η values no longer pin down the E0 strength. The paper gives no microscopic derivation or empirical test for this scaling. That doesn't kill the paper—the 44,48Ca E0 agreement is real—but it does mean the central 'direct correlation' claim is conditional on an ad hoc assumption, not established. A referee should ask for a sensitivity test: vary κ_int around the scaling and show how much ρ²(E0) moves.\n\nMinor issues: the agreement with isotopic shifts is visual, with no error bars or χ² metric; the IBM-1 code is unpublished (ref [29]); and the per-chain fitting of κ and η means the radius prediction is partly a fit, as the reader's report notes. None of these are disqualifying for a phenomenological paper.\n\nBottom line: this deserves a serious referee. It is a solid, useful step for the Ca radius puzzle, and the E0 link is worth testing further. The main thing I'd want from a revision is a robustness analysis of the κ_int scaling and a quantitative goodness-of-fit. I'd bring it to a reading group, and I'd cite it for the IBM-CM radii application.","headline":"First IBM-CM treatment of the Ca charge-radius puzzle reproduces the 48Ca dip and links radii to E0 strengths, but the central correlation rests on an untested κ_int scaling assumption that a referee should probe.","tokens_in":15015,"tokens_out":2743,"would_cite":true,"duration_ms":28073,"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":"The paper claims that one configuration-mixing mechanism, fixed by energy levels, reproduces both the calcium charge-radius isotopic shift and the electric monopole strengths $\\rho^2(E0)$, establishing a direct size–E0 correlation.","keywords":["interacting boson model","configuration mixing","charge radii","isotopic shift","electric monopole transition","0+ states","calcium isotopes","shape coexistence"],"falsifier":"Measure the electric monopole strength $\\rho^2(E0;0_2^+\\to0_1^+)$ in $^{46}$Ca: the model's parameters are already fixed by the Ca isotopic shifts, so Eq. (6) gives a definite value, and a significant deviation would show the claimed size–E0 correlation is not quantitative.","tokens_in":13830,"feed_emoji":"⚛️","tokens_out":11487,"duration_ms":96959,"temperature":0.7,"pith_summary":"Charge radii of the calcium isotopes behave oddly: $^{40}$Ca and $^{48}$Ca have nearly equal radii, $^{44}$Ca is a local maximum, and the radius rises steeply beyond $^{48}$Ca. This paper argues that the interacting boson model with configuration mixing (IBM-CM), with wave functions fixed by fitting energy levels, reproduces this isotopic shift, in particular the $^{48}$Ca dip and the $^{50}$Ca rise. It then uses the same deformation parameters that fix the radii to compute electric monopole strengths $\\rho^2(E0;0_2^+\\to0_1^+)$; the results agree with data for $^{44}$Ca and $^{48}$Ca. If correct, the argument establishes a direct correlation between nuclear size and monopole transitions and makes unmeasured $E0$ strengths predictable from measured radii.","feed_headline":"Same mixing explains calcium radius dip and E0 strengths","feed_subtitle":"Wave functions fitted to energy levels reproduce the 48Ca radius dip and the E0 strengths of 44Ca and 48Ca.","key_machinery":"The central machinery is the IBM-CM Hamiltonian $\\hat H' = \\hat H_{\\rm reg} + (\\hat H_{\\rm int}+\\Delta)+\\hat V_{\\rm mix}$ acting in the direct sum of the regular $[n_b]$ and intruder $[n_b+2]$ boson spaces, together with the charge-radius operator $\\hat T(r^2) = \\langle r^2\\rangle_c + \\kappa\\,\\hat n_b + \\eta\\,\\hat n_d/\\hat n_b$. The load-bearing identity is the linear-scaling constraint $\\kappa_{\\rm int} = \\frac{n_b}{n_b+2}\\kappa_{\\rm reg}$, which keeps the radius linear when moving from the regular to the intruder space. Once the wave functions are fixed from energy levels, the isotopic-shift formula and the E0 formula express both observables in terms of the same $\\eta_{\\rm reg}$ and $\\eta_{\\rm int}$, so the radius and the monopole strength are locked together by construction. The mixing probabilities of the ground and $0_2^+$ states carry the physical content: the $^{48}$Ca ground state is 78% intruder and $^{50}$Ca is 98% intruder.","core_discovery":"The central claim is that the location of the $0_2^+$ state and the size of the nucleus are not independent: both are governed by the mixing between the regular $[n_b]$ (0p-0h) and intruder $[n_b+2]$ (2p-2h) boson spaces. In IBM-CM, wave functions fitted to level energies give ground-state mixing probabilities in which $^{48}$Ca is 78% intruder and $^{50}$Ca is 98% intruder, and the radius operator with the imposed linear-scaling constraint then predicts the Ca isotopic shift, including the $^{48}$Ca dip. Taking the same $\\eta_{\\rm reg}$ and $\\eta_{\\rm int}$ from the isotopic-shift fit, the E0 formula yields $\\rho^2(E0)$ values for $^{44}$Ca and $^{48}$Ca in reasonable agreement with experiment; thus the paper establishes a direct correlation between nuclear size and electric monopole transitions. The mechanism is that the $^{48}$Ca ground state is strongly intruder-dominated but blocked, giving a small radius, while $^{50}$Ca is almost pure intruder, giving the steep rise; the same wave functions explain why the $0_2^+$ state is isomeric in $^{40}$Ca but less so in $^{48}$Ca. Calculations for the Ar and Ti isotopes show a similar, though weaker, connection.","pith_inferences":["A testable extension is to apply the same parameter-sharing to other mass regions with known charge radii and E0 data, such as the tin or lead isotopes; if the correlation persists there, it is a structural feature of configuration mixing rather than a Ca-specific fit.","The linear-scaling constraint is a modeling assumption; a shell-model calculation that computed the intruder 2p-2h space's radius directly could verify or falsify the $\\kappa_{\\rm int}$ relation without new experiment.","The paper's failure to describe $^{42}$Ca suggests the two-configuration truncation may be too restrictive there; adding 4p-4h admixtures or a different intruder Hamiltonian could be the next step.","A sharp, implied prediction is the sign and magnitude of $\\rho^2(E0)$ in $^{46}$Ca, which is measurable with current electron-conversion techniques and would provide a clean test of the size–E0 link."],"forward_implications":["If the correlation is correct, the $^{48}$Ca radius dip and the $^{50}$Ca rise require no new physics beyond two-configuration mixing; the same wave functions that fit energies already produce them.","Because the same $\\eta_{\\rm reg}$ and $\\eta_{\\rm int}$ enter the radius and E0 formulas, measured charge radii become predictions for unmeasured $\\rho^2(E0)$ values, for example in $^{46}$Ca.","The $0_2^+$ energy becomes a size-relevant observable: interactions that misplace the $0_2^+$ state, as the paper notes for the $^{36,38}$Ca shell-model results, will also fail on radii.","The binding-energy argument yields a qualitative rule for other chains: isotopes with similar binding energy per nucleon, such as $^{40}$Ca and $^{48}$Ca, should have similar radii, while unequal pairs such as $^{38}$Ar and $^{46}$Ar should differ."],"supporting_citations":[{"why":"It introduces the configuration-mixing version of the interacting boson model (IBM-CM) that defines the regular-plus-intruder framework used throughout the paper.","marker":"[26]"},{"why":"It provides the interacting boson model and the standard charge-radius operator form on which Eq. (4) is based.","marker":"[9]"},{"why":"It supplies the linear-scaling constraint $\\kappa_{\\rm int} = n_b/(n_b+2)\\kappa_{\\rm reg}$ used to keep the calculated radius linear across configurations.","marker":"[33]"},{"why":"It gives the fitting procedure and the IBM-CM Hamiltonian parameters for the $^{42-46}$Ca isotopes from which the present calculations start.","marker":"[19]"},{"why":"It provides the experimental $\\rho^2(E0;0_2^+\\to0_1^+)$ values for $^{40,42,44,48}$Ca against which the computed E0 strengths are compared.","marker":"[25]"},{"why":"It supplies the measured charge radius of the proton-rich isotope $^{36}$Ca used in the isotopic-shift comparison.","marker":"[7]"},{"why":"It supplies the measured charge radii of the neutron-rich calcium isotopes, including the $^{48}$Ca dip and the steep rise beyond, which the model aims to reproduce.","marker":"[8]"},{"why":"It provides the shell-model isotopic-shift results used as a comparison and the interaction whose $0_2^+$ location is tied to the radius behavior.","marker":"[46]"},{"why":"It supplies the Hartree-Fock mass-table isotopic shifts shown for comparison, which miss the $^{48}$Ca dip.","marker":"[51]"}],"fun_headline_variants":["One mixing mechanism links 48Ca size dip and E0 strengths","Mixing predicts calcium radii and monopole strengths","Calcium radius dip and monopole strengths share a single origin","Nuclear size and E0 strengths: one mixing controls both"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the radius contribution of the intruder $[n_b+2]$ configuration scales with the regular $[n_b]$ configuration through the imposed linear relation $\\kappa_{\\rm int} = \\frac{n_b}{n_b+2}\\kappa_{\\rm reg}$; if a microscopic calculation showed otherwise, the fitted $\\eta$ values and the derived $\\rho^2(E0)$ strengths would lose their meaning.","fun_headline_variants_meta":{"raw":{"variants":["One mixing mechanism links 48Ca size dip and E0 strengths","Mixing predicts calcium radii and monopole strengths","Calcium radius dip and monopole strengths share a single origin","Nuclear size and E0 strengths: one mixing controls both"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000567,"raw_usage":{"total_tokens":2685,"prompt_tokens":941,"completion_tokens":1744,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":557,"completion_tokens_details":{"reasoning_tokens":1676}},"tokens_in":557,"tokens_out":1744,"duration_ms":11470,"temperature":1.0,"reasoning_tokens":1676,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T20:37:53.438978+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the electric monopole strength $\\rho^2(E0;0_2^+\\to0_1^+)$ in $^{46}$Ca: the model's parameters are already fixed by the Ca isotopic shifts, so Eq. (6) gives a definite value, and a significant deviation would show the claimed size–E0 correlation is not quantitative.","supporting_citations":[{"cited_title":"nndc.bnl.gov/ensdf/","cited_arxiv_id":null,"evidence_quote":"It introduces the configuration-mixing version of the interacting boson model (IBM-CM) that defines the regular-plus-intruder framework used throughout the paper."},{"cited_title":"Iachello and A","cited_arxiv_id":null,"evidence_quote":"It provides the interacting boson model and the standard charge-radius operator form on which Eq. (4) is based."},{"cited_title":"Elliott and J.A","cited_arxiv_id":null,"evidence_quote":"It supplies the linear-scaling constraint $\\kappa_{\\rm int} = n_b/(n_b+2)\\kappa_{\\rm reg}$ used to keep the calculated radius linear across configurations."},{"cited_title":"Maheshwari and K","cited_arxiv_id":null,"evidence_quote":"It gives the fitting procedure and the IBM-CM Hamiltonian parameters for the $^{42-46}$Ca isotopes from which the present calculations start."},{"cited_title":"Schulz, Ann","cited_arxiv_id":null,"evidence_quote":"It provides the experimental $\\rho^2(E0;0_2^+\\to0_1^+)$ values for $^{40,42,44,48}$Ca against which the computed E0 strengths are compared."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the measured charge radius of the proton-rich isotope $^{36}$Ca used in the isotopic-shift comparison."},{"cited_title":"Otsuka, T","cited_arxiv_id":null,"evidence_quote":"It supplies the measured charge radii of the neutron-rich calcium isotopes, including the $^{48}$Ca dip and the steep rise beyond, which the model aims to reproduce."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It provides the shell-model isotopic-shift results used as a comparison and the interaction whose $0_2^+$ location is tied to the radius behavior."},{"cited_title":"Retamosa, E","cited_arxiv_id":null,"evidence_quote":"It supplies the Hartree-Fock mass-table isotopic shifts shown for comparison, which miss the $^{48}$Ca dip."}],"review_version":1}