{"id":"6c62e8f4-8a4b-4d3c-8026-dddc2318fad2","arxiv_id":"2505.22412","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Applying a neutron-proton correction around the Fermi surface in a relativistic Hartree-Bogoliubov model sharpens charge-radius kinks at N=8, 20, and 28 and improves radii for meson-exchange interactions, but not for density-dependent interactions.","lead":"This paper tests whether an ad hoc neutron-proton pairing correction to the charge radius formula, with constants fitted to heavier nuclei, improves agreement with measured radii in light oxygen-to-argon isotopes and reveals shell closures at N=8, 20, and 28. A generalist reader might care because it probes how well simple mean-field models can capture nuclear magic numbers, a basic organizing feature of nuclei.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The charge-radius kinks at N=8, 20, and 28 are largely built into Eq. (6): adding a positive term proportional to |D_n - D_p| automatically raises open-shell radii more than closed-shell ones, and D_n, D_p are not mixed neutron-proton pair amplitudes, so the 'np-pairing enhancement' claim is…","rationale":"The reader's weakest_assumption was transferability of a0 and delta and the spherical-32Mg problem. My main concern is more fundamental and partially overlaps: the correction term's functional form, not just its fitted constants, generates the kinks. This is the load-bearing point because the abstract's central sentence attributes the kinks to neutron-proton pairing. If the term were refit to this dataset alone, it might still fit, but the physical statement requires that DeltaD specifically encodes neutron-proton correlations. Eq. (7) shows it does not: it is a difference of two like-particle pairing traces. The BCS-vs-Bogoliubov comparison for 26Mg is a genuine piece of support, because it shows the treatment of pairing changes D_p and hence the radius in the right direction. However, it does not validate the 'np pairing' label, nor does it establish that the N=20 or N=28 kinks are physical. The deformation problem for 32Mg is a separate symptom of the same underdetermination: a spherical mean field plus an ad hoc radius shift can mimic the experimental radius of a deformed nucleus. My proposed refit and null-proxy comparison would settle whether the constants are transferable and whether the kinks are specific to DeltaD; if they are not, the paper should be revised to present the correction as a phenomenological radius adjustment rather than evidence for enhanced shell closures. This keeps the verdict CONDITIONAL: the calculation is reproducible in structure and the BCS/Bogoliubov comparison is informative, but the central interpretation needs an out-of-sample test.","tokens_in":15988,"tokens_out":7357,"duration_ms":96156,"concrete_test":"Refit a0 (and delta if odd-odd data are included) using only the 24 even-even O-Ar charge radii in Table 1, with leave-one-out cross-validation, and compare the predictive error against the same procedure with a null correction based on a smooth proxy such as |N-Z|/A or |lambda_n - lambda_p|. If the refitted a0 is within 1 sigma of zero, or if the proxy yields the same kinks at N=8, 20, and 28, the claimed enhancement is not specific to the np-pairing ansatz.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on the correction term in Eq. (6), r_ch^2 = <r_p^2> + 0.7056 + a0 * DeltaD / sqrt(A) + delta / sqrt(A), with DeltaD = |D_n - D_p| defined in Eq. (7). Two properties of this term carry most of the weight. First, D_n and D_p are separate like-particle pairing traces, sum u_k v_k; they contain no mixed neutron-proton anomalous density. The label 'neutron-proton pairing correlation' is therefore an interpretation, not a consequence of the calculation. Second, the correction is always positive and is small when both pairing traces are small, which is precisely the closed-shell situation. Adding such a term will deepen existing minima at magic neutron numbers even if no physical np correlation is present. The observed 'sudden strengthening' at N=8, 20, and 28 is thus to a large extent guaranteed by the ansatz's functional form. What is missing is a demonstration that the effect survives when the constants are not globally pre-fit, or when a genuine np-pairing observable is used. The paper's own Table 1 shows the correction helps meson-exchange interactions but worsens DD-PC1, so the effect is interaction-dependent; without an out-of-sample test, the shell-closure enhancement claim is underdetermined rather than demonstrated.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes ground-state charge radii for even-even O, Ne, Mg, Si, and Ar isotopes using the multidimensionally-constrained relativistic Hartree-Bogoliubov model with four effective interactions (DD-ME2, DD-PC1, NL3, PK1). It augments the standard charge-radius formula with a correction term, Eq. (6), built from the absolute difference of neutron and proton pairing traces defined in Eq. (7), with constants taken from a previous global fit. The authors claim that this 'neutron-proton correction around the Fermi surface' produces sudden enhancements of charge radii at N=8, 20, and 28, thereby reinforcing shell closures at these neutron numbers, and that for Mg isotopes the Bogoliubov treatment of pairing is more consistent with the measured N=14 radius than BCS. They quantify the overall agreement with 24 measured radii through chi-square and rms deviations in Table 1.","tokens_in":16270,"tokens_out":4222,"duration_ms":53435,"significance":"If the central claim were established, the paper would offer a simple phenomenological route to identify shell closures in the sd shell from charge radii, and it would pinpoint a specific deficiency of BCS pairing in neutron-rich Mg isotopes. The work has concrete strengths: four effective interactions are compared, the BCS-versus-Bogoliubov comparison for the Mg chain in Fig. 2 is instructive, and the tabulated deviations in Table 1 give a transparent measure of global performance. However, the central physical interpretation is not supported by the computed quantity: Eq. (7) contains separate like-particle pairing traces, not a neutron-proton mixed pairing amplitude, and the correction term in Eq. (6) is positive-definite with constants calibrated on the same type of data, so the appearance of kinks at magic numbers is partially built into the ansatz rather than emerging from a physically independent neutron-proton pairing effect.","major_comments":[{"comment":"The quantities D_n and D_p are defined as sums of u_k v_k over neutron and proton quasiparticle states separately. These are like-particle pairing traces; they contain no neutron-proton anomalous density or mixed isospin amplitude. The paper repeatedly refers to 'neutron-proton pairing correlations around the Fermi surface' and claims that these correlations 'lead to a sudden strengthening' of charge radii. That causal statement is not supported by the calculation. Either a genuine neutron-proton pairing observable must be computed, or the terminology and all associated claims must be revised to describe a difference of like-particle pairing strengths.","section":"Sec. 2, Eq. (7)"},{"comment":"The constants a0 = 0.561 and delta = 0.355 are taken unchanged from the global fit of Ref. [62], which was itself fit to charge-radii data. The correction term is therefore not an independent prediction; its magnitude is calibrated on the same class of observables that the paper then compares with experiment. The manuscript provides no out-of-sample test, no refit restricted to the sd shell, and no sensitivity analysis for these constants. This is load-bearing for the claim that the correction 'enhances' shell closures, because the size of the kinks depends directly on the pre-fitted values.","section":"Sec. 2, Eq. (6)"},{"comment":"The central assertion of sudden strengthening at N=8, 20, and 28 is based on visual inspection of kinks in Fig. 1, without any quantitative measure such as a second difference of r_ch, a slope change, or an uncertainty estimate. Since the correction term in Eq. (6) is positive and depends on |D_n - D_p|, it will generically tend to raise radii more in open-shell nuclei than in doubly magic or closed-neutron-shell nuclei, so visual kinks are not by themselves evidence of a physical shell-closure enhancement. A quantitative definition of 'kink' and a demonstration that the pattern is not an artifact of the functional form are needed.","section":"Sec. 3, Fig. 1"},{"comment":"Table 1 shows that the RHB* correction improves the average deviation for the meson-exchange interactions NL3 and PK1, but the rms deviation Δ increases for both density-dependent interactions DD-ME2 and DD-PC1, and the chi-square for DD-PC1 worsens substantially (952.21 to 1083.08). This interaction dependence is acknowledged in the text but not explained. Because the claimed universal enhancement of shell closures is expected to be interaction-independent, the paper should either identify why density-dependent interactions behave differently or restrict the claim to meson-exchange functionals.","section":"Sec. 3, Table 1"},{"comment":"The manuscript notes that the RHB model yields spherical ground states for 32Mg for all four forces, despite the well-established deformation in the island of inversion, and that only beyond-mean-field effects restore deformation. Nevertheless, the paper uses the RHB* charge radius of 32Mg as evidence for an N=20 shell-closure effect. If the mean-field ground state is qualitatively wrong in this region, the inferred enhancement at N=20 cannot be interpreted as physical without a deformation or projection treatment. This should be addressed before any conclusion about N=20 magicity is drawn.","section":"Sec. 3, Mg isotopes and N=20"}],"minor_comments":[{"comment":"The phrase 'A ansatz' should read 'An ansatz', and 'Acnowledgements' in the back matter is a typo for 'Acknowledgements'.","section":"Abstract and Sec. 1"},{"comment":"The symbol N is used both for the neutron number and for the number of data points in the definition of the average deviation; this is confusing and should be relabeled, for example as N_data.","section":"Sec. 3, Eq. (9)"},{"comment":"The caption states that the gray band represents N=14 and 20, but the text also discusses N=16 as a relevant magic number; the figure shading and caption should be clarified so that all discussed shell closures are identified consistently.","section":"Fig. 1 caption"},{"comment":"The notation 'r_ch^2 = <r_p^2> + ...' is dimensionally consistent, but the paper should state explicitly that the second term 0.7056 fm^2 is the standard proton finite-size correction and cite the source consistently with Eq. (6) of Ref. [62].","section":"Sec. 2, Eq. (6)"}],"recommendation":"major_revision","confidential_remarks":"The paper is an incremental phenomenological study with a significant terminological problem: the computed D_n and D_p are not neutron-proton pairing amplitudes, so the central physical claim is overstated. A major revision that reframes the correction as an empirical difference of like-particle pairing strengths, adds quantitative kink diagnostics and an out-of-sample test, and restricts the conclusions to the interactions for which the correction actually helps could make the manuscript publishable. I would also flag that the self-citation pattern to Refs. [54, 62] is heavy but not inappropriate given that the ansatz originates there."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper is a straightforward extension of the authors' earlier charge-radius correction ansatz to sd-shell nuclei. What is new: they apply the same neutron-proton 'correction' term to O, Ne, Mg, Si, Ar chains within an RHB framework, compare BCS and Bogoliubov pairing treatments for Mg, and test four interactions. The BCS-versus-Bogoliubov comparison for the N=14 radius in Mg is the cleanest part, and they are admirably honest that the correction helps meson-exchange interactions but not DD-PC1.\n\nThe soft spot is the central claim. The correction term in Eq. (6) is |D_n - D_p|, where D_n and D_p are separate like-particle pairing traces. There is no mixed neutron-proton anomalous density in the calculation, so calling it 'neutron-proton pairing correction' is an interpretation, not something the calculation demonstrates. And because the term is positive and small when both pairing traces are small, it will automatically deepen minima at closed shells. The kinks at N=8, 20, 28 are therefore largely built into the functional form, not evidence for a new physical mechanism. The constants a0 and delta come from a global fit in Ref. [62], so the magnitude is not independent of the data being compared. That said, the kinks are driven by the model-computed Delta D, so it is not fully circular; the paper just does not establish that the ansatz captures genuine np pairing.\n\nOther soft spots: the model restricts to beta20 symmetry and gives a spherical ground state for 32Mg, which the authors themselves note. That weakens the N=20 conclusion. There is no uncertainty quantification on the kinks. These are real but not fatal if the paper is read as a phenomenological study of a specific ansatz rather than a derivation.\n\nWho is this for? Practitioners in nuclear charge radii and mean-field models who want to see how this correction behaves across the sd shell. It deserves a serious referee; the interaction dependence and the BCS-RHB comparison are worth publishing, but the shell-closure-enhancement language should be tempered and the np-pairing interpretation justified or softened.\n\nRecommendation: send to peer review, but the referee should ask for a clearer statement that the correction is an empirical ansatz, and ideally an out-of-sample test or a demonstration with a genuine mixed-pairing observable.","headline":"A useful systematic extension of the authors' charge-radius ansatz to the sd shell, but the headline claim that the correction enhances N=8, 20, 28 shell closures is largely built into the functional form of the correction term.","tokens_in":16844,"tokens_out":2059,"would_cite":false,"duration_ms":22671,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A neutron-proton correction near the Fermi surface strengthens the charge-radius signature of shell closures at N=8, 20, and 28 in light nuclei.","keywords":["charge radius","magic numbers","neutron-proton pairing","relativistic Hartree-Bogoliubov","sd shell","shell closure","N=14 subshell","N=20 island of inversion"],"falsifier":"Refit the two constants $a_0$ and $\\delta$ using the 24 measured radii of the O–Ar chains (or a subset excluding the magic nuclei) and recompute the radii: if the sharp kinks at $N=8$, 20, and 28 disappear or shift when the constants are determined from non-magic nuclei only, the claimed shell-closure enhancement is an artifact of the global fit rather than a physical signal. Alternatively, measure charge radii for neutron-rich Ne or Mg isotopes beyond $N=20$ (e.g., $^{34}$Mg or $^{34-36}$Ne): the RHB* model predicts a specific kink pattern there, and disagreement would falsify the ansatz's predictive power.","tokens_in":15670,"feed_emoji":"⚛️","tokens_out":13929,"duration_ms":125248,"temperature":0.7,"pith_summary":"This paper proposes that a simple correction term, representing neutron-proton correlations among quasiparticle states near the Fermi surface, makes nuclear charge radii into sharper fingerprints of shell closures. Applied to even-even O, Ne, Mg, Si, and Ar isotopes within a relativistic Hartree-Bogoliubov model, the term produces sudden upward steps in the radius chains at neutron numbers $N=8$, 20, and 28, which the authors interpret as enhancement of the corresponding shell closures. The same correction resolves a known failure at $N=14$ in the Mg chain: BCS pairing overestimates the shell effect, while the Bogoliubov treatment, which creates fractional proton occupations near the Fermi surface, matches the measured $^{26}$Mg radius. The correction improves agreement with 24 measured radii when meson-exchange effective interactions (NL3, PK1) are used, but not for density-dependent interactions (DD-ME2, DD-PC1), indicating the effect is not universal across functionals.","feed_headline":"Sharpened charge-radius kinks flag magic numbers at N=8, 20, 28","feed_subtitle":"A Fermi-surface pairing term amplifies radius jumps at magic neutron numbers and fixes the N=14 puzzle in Mg.","key_machinery":"The load-bearing object is the modified rms charge-radius formula Eq. (6), $r_{\\rm ch}^2 = \\langle r_p^2 \\rangle + 0.7056\\ {\\rm fm}^2 + \\frac{a_0}{\\sqrt{A}}\\Delta D\\ {\\rm fm}^2 + \\frac{\\delta}{\\sqrt{A}}\\ {\\rm fm}^2$, where $\\Delta D = |D_n - D_p|$ measures neutron-proton correlations around the Fermi surface through $D_{n,p} = \\sum_{k>0} u_k v_k$ summed over quasiparticle levels with $|E_k - \\lambda| < 20$ MeV. The term converts the non-integer occupation of proton and neutron quasiparticle orbitals into a radius shift, and it is this $\\Delta D$ that injects a kink into the charge-radius chain whenever a shell closure rearranges the Fermi-surface occupations. The correction is applied on top of the multidimensionally-constrained relativistic Hartree-Bogoliubov model restricted to axial/reflection symmetry (only $\\beta_{20}$), using a separable finite-range pairing force with strength $G = 728$ MeV fm$^3$ and range $a = 0.644$ fm, and tested with four effective interactions (PK1, NL3, DD-ME2, DD-PC1) to gauge parameter dependence.","core_discovery":"The central claim is that the modified charge-radius formula Eq. (6), with the neutron-proton correlation term $\\Delta D = |D_n - D_p|$ constructed from quasiparticle occupations around the Fermi surface, turns charge-radii evolution in the $sd$ shell into a clear signal of shell structure. For the five isotopic chains with proton numbers $Z = 8, 10, 12, 14, 18$, the correction suddenly strengthens the charge radii at $N = 8$, 20, and 28, and the authors read this as an enhancement of the shell closure at those neutron numbers. In the Mg isotopes, the way pairing is treated determines the $N=14$ signal: BCS overestimates the shell effect, whereas the Bogoliubov transformation yields fractional proton occupations near the Fermi surface and reproduces the measured $^{26}$Mg radius. The authors further report that the correction improves the description of 24 measured charge radii when meson-exchange effective interactions (NL3, PK1) are used, but does not significantly help the density-dependent interactions (DD-ME2, DD-PC1).","pith_inferences":["A direct refit of the constants $a_0$ and $\\delta$ from the 24 $sd$-shell radii alone would test whether the $N=8$, 20, and 28 kinks are driven by the chosen global fit or by shell structure itself; large changes in the constants would weaken the physical interpretation.","Because the unprojected model gives a spherical ground state for $^{32}$Mg while the measured radius is better reproduced with the correction, angular-momentum projected RHB* calculations could show whether the $N=20$ kink survives when deformation is restored.","The failure of the correction for density-dependent interactions suggests the microscopic content may be an effective isovector surface term rather than literal neutron-proton pairing; comparing with ab initio charge radii for $^{24-32}$Mg would decide which interpretation is more natural.","A systematic scan of the same $\\Delta D$ term across heavier mass regions (for example around $N=50$ or $N=82$) would show whether the ansatz predicts spurious kinks at non-magic neutron numbers, constraining how widely the formula can be used."],"forward_implications":["If the correction is physical, charge-radius kinks at $N=8$, 20, and 28 in the O–Ar chains can serve as a relatively inexpensive experimental signature of shell closures in neutron-rich nuclei.","The $N=14$ charge radius in the Mg isotopes becomes a diagnostic for how pairing is treated: only the Bogoliubov treatment with fractional proton occupations near the Fermi surface matches the measured value, indicating BCS-type approximations are inadequate for radii in this region.","Because the correction improves meson-exchange effective interactions but not density-dependent ones, the work implies that meson-exchange covariant functionals are missing an isovector surface term that density-dependent functionals already effectively contain.","Applying the same correction to heavier isotopic chains, which the authors state is in progress, should reveal whether additional or shifted shell closures appear in regions with new magic numbers.","The results support the general statement that neutron-proton correlation near the Fermi surface leaves a measurable imprint on the nuclear charge radius, not just like-particle pairing."],"supporting_citations":[{"why":"Supplies the modified charge-radius formula and the fitted constants $a_0=0.561$ and $\\delta=0.355$ used in Eq. (6).","marker":"[62]"},{"why":"Introduces the neutron-proton correlation extracted from quasiparticle states around the Fermi surface in the MDC-RHB framework, the method this paper applies to the $sd$ shell.","marker":"[58]"},{"why":"Provides the earlier RMF(BCS)* ansatz that overestimates the $N=14$ charge radius in Mg, motivating the switch to the Bogoliubov treatment.","marker":"[54]"},{"why":"Defines the separable finite-range pairing force with strength $G=728$ MeV fm$^3$ and range $a=0.644$ fm used to build the pairing field.","marker":"[61]"},{"why":"Gives the experimental ground-state charge radii of the O, Ne, Si, and Ar isotopic chains against which the RHB and RHB* results are compared.","marker":"[3]"},{"why":"Provides the experimental Mg charge radii, including the key $N=14$ datum at $^{26}$Mg that distinguishes the BCS and Bogoliubov pairing treatments.","marker":"[22]"}],"fun_headline_variants":["Fermi-surface pairing sharpens charge radii at magic numbers","Neutron-proton correction flags shell closures in sd-shell radii","Charge-radius jumps at N=8,20,28 traced to Fermi-surface pairing","Pairing term in charge radii solves Mg N=14 puzzle","Shell quenching visible in charge radii via neutron-proton term"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The correction term in Eq. (6) is applied with the constants $a_0 = 0.561$ and $\\delta = 0.355$ taken from a global fit in Ref. [62], so the whole analysis assumes these constants are transferable to $sd$-shell nuclei and that $\\Delta D$ really captures physical neutron-proton pairing; if the constants are not universal, or if deformation and beyond-mean-field effects are essential (the model returns spherical ground states even for $^{32}$Mg), the inferred shell-closure enhancement may be an artifact of the ansatz.","fun_headline_variants_meta":{"raw":{"variants":["Fermi-surface pairing sharpens charge radii at magic numbers","Neutron-proton correction flags shell closures in sd-shell radii","Charge-radius jumps at N=8,20,28 traced to Fermi-surface pairing","Pairing term in charge radii solves Mg N=14 puzzle","Shell quenching visible in charge radii via neutron-proton term"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000234,"raw_usage":{"total_tokens":1569,"prompt_tokens":1088,"completion_tokens":481,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":704,"completion_tokens_details":{"reasoning_tokens":391}},"tokens_in":704,"tokens_out":481,"duration_ms":5255,"temperature":1.0,"reasoning_tokens":391,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:07:54.652911+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Refit the two constants $a_0$ and $\\delta$ using the 24 measured radii of the O–Ar chains (or a subset excluding the magic nuclei) and recompute the radii: if the sharp kinks at $N=8$, 20, and 28 disappear or shift when the constants are determined from non-magic nuclei only, the claimed shell-closure enhancement is an artifact of the global fit rather than a physical signal. Alternatively, measure charge radii for neutron-rich Ne or Mg isotopes beyond $N=20$ (e.g., $^{34}$Mg or $^{34-36}$Ne): the RHB* model predicts a specific kink pattern there, and disagreement would falsify the ansatz's predictive power.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the modified charge-radius formula and the fitted constants $a_0=0.561$ and $\\delta=0.355$ used in Eq. (6)."},{"cited_title":"Yang, Y .-T","cited_arxiv_id":null,"evidence_quote":"Introduces the neutron-proton correlation extracted from quasiparticle states around the Fermi surface in the MDC-RHB framework, the method this paper applies to the $sd$ shell."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the earlier RMF(BCS)* ansatz that overestimates the $N=14$ charge radius in Mg, motivating the switch to the Bogoliubov treatment."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the separable finite-range pairing force with strength $G=728$ MeV fm$^3$ and range $a=0.644$ fm used to build the pairing field."},{"cited_title":"Angeli, K","cited_arxiv_id":null,"evidence_quote":"Gives the experimental ground-state charge radii of the O, Ne, Si, and Ar isotopic chains against which the RHB and RHB* results are compared."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the experimental Mg charge radii, including the key $N=14$ datum at $^{26}$Mg that distinguishes the BCS and Bogoliubov pairing treatments."}],"review_version":1}