{"id":"62dfe48d-88de-4663-bfdb-18f9f439fcfb","arxiv_id":"2501.19037","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Closed-form perturbative expressions for spin and pseudo-spin energy splittings in relativistic mean-field models are derived, decomposing pseudo-spin breaking into a mass-dependent central-potential part and a nearly constant spin-orbit part.","lead":"Spin and pseudo-spin symmetries shape the energy levels of atomic nuclei, and this paper derives simple perturbative formulas for how those symmetries break. The formulas separate the roles of the central and spin-orbit potentials and explain how the energy gaps scale with nuclear mass number.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (12)'s ΔEg term appears sign-flipped: integration by parts of Eq. (4) gives −∫(Ga−Gb)Σ′dr, not +∫(Ga−Gb)Σ′dr, so the central PSS closed form and its A-scaling interpretation need verification.","rationale":"The reader's weakest assumption concerns unquantified Woods-Saxon fits and the single N=50 chain; that is a legitimate robustness issue but secondary. The sign of the leading closed-form PSS term is more fundamental: the paper's main deliverable is Eq. (12), and the entire pseudo-spin mechanism section interprets the sign and overlap of the two integrand terms. If the sign is wrong, the A-scaling conclusions and the identification ΔEHO=ΔEg, ΔESO=−ΔEf in Eq. (15) inherit the error. I do not claim the numerical calculations are necessarily wrong; a typographical sign slip in Eq. (12) is plausible. But because no code or data are provided, the reader cannot tell from the manuscript whether Fig. 5 was produced with the printed formula or with the corrected sign. The proposed test settles this and is inexpensive. Therefore the verdict remains conditional, but the condition is now a definite algebraic check on the central formula rather than only a request for fit uncertainty quantification.","tokens_in":1172,"tokens_out":809,"duration_ms":189859,"concrete_test":"For the ν(2d5/2−1g7/2) PSS pair in the N=50 chain, take the SS-RHO reference states and RHB Σ,Δ potentials, and compute the first-order gap three ways: (1) direct evaluation of Eq. (4); (2) Eq. (12) as printed; (3) Eq. (12) with the sign of the ΔE_g term reversed. If the values of (1) and (3) match and (2) does not, the sign error is confirmed and the Fig. 5 decomposition and the 'positive lobe' overlap argument must be re-derived. Also evaluate ∫r^2(g_a^2−g_b^2)Σ_HO to confirm it vanishes; if not, Eq. (12) misses a state-dependent reference contribution.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central PSS formula in Eq. (12) is stated as ΔE_PSS = ∫(G_a−G_b) dΣ/dr dr + ∫(F_a−F_b) dΔ/dr dr. The lower-component term is consistent with Eq. (4): with B = d0−Δ, integration by parts gives +∫(F_a−F_b)Δ′ dr because B′ = −Δ′. The upper-component term does not follow the same sign. Starting from Eq. (4), the g/Σ contribution is M = ∫ r^2(g_a^2−g_b^2)(Σ−Σ_HO) dr. If the Σ_HO integral vanishes, M = ∫ r^2(g_a^2−g_b^2)Σ dr. Defining A = G_a−G_b with A′ = r^2(g_a^2−g_b^2), and using Σ(∞)=0, integration by parts yields M = −∫ A Σ′ dr, the negative of the first term printed in Eq. (12). This sign matters because Section IV.C.1 interprets the growth of ΔE_g with A as Σ′ increasingly probing the positive lobe of G_a−G_b; with the corrected sign that overlap would decrease ΔE_g. Thus either the closed form or the mechanism argument is reversed, and the reader's quoted central formula is not supported as written. The issue is directly checkable numerically.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a first- and second-order perturbative treatment of spin-symmetry (SS) and pseudo-spin-symmetry (PSS) breaking in relativistic mean-field models, using the spin-symmetric relativistic harmonic oscillator as the reference. It derives closed-form, first-order expressions for the SS and PSS energy splittings, analyzes their mass-number dependence along the N=50 isotonic chain, and concludes that the SS splitting scales as A^(-2/3) through two separate A^(-1/3) factors, while the PSS splitting is governed by a g/Σ contribution that grows with A and an f/Δ contribution that is approximately constant. The results are benchmarked against RHB calculations and experimental data.","tokens_in":12364,"tokens_out":19827,"duration_ms":193800,"significance":"If the sign issue in Eq. (12) is resolved, the unified perturbative framework is a valuable qualitative tool: it is not fitted to the energy gaps, it cleanly separates the g/Σ and f/Δ contributions to PSS breaking, and it provides a transparent covariant account of the known A^(-2/3) spin-orbit trend. The N=50 benchmark and the explicit checks of the perturbation parameter are strengths. The generality of the PSS scaling claims, however, currently depends on a single orbital pair and on Woods-Saxon fits whose quality is not reported, so the quantitative conclusions should be treated with caution until those points are addressed.","major_comments":[{"comment":"The upper-component term in Eq. (12) has the wrong sign. From Eq. (4), the g-dependent contribution is I_g = ∫ r²(g_a²−g_b²)(Σ−Σ_HO) dr. The Σ_HO part vanishes for degenerate PSS partners, as stated in Section IV.B. Defining A(r)=∫_0^r r′²(g_a²−g_b²) dr′ (so A′=r²(g_a²−g_b²)), integration by parts gives I_g = ∫ A′Σ dr = [AΣ]_0^∞ − ∫ A Σ′ dr = −∫(G_a−G_b)Σ′ dr, because Σ(∞)=0 and A(0)=0. The printed Eq. (12) has +∫(G_a−G_b)Σ′ dr. This is not a cosmetic sign: with the shapes displayed in Fig. 6, where G_a−G_b is negative in the core and positive near the surface while Σ′ is surface-peaked, the corrected formula produces an overlap trend with A that is opposite to the mechanism described in Section IV.C.1. Please correct the sign and re-evaluate the ΔE_g curve in Fig. 5, the interpretation in Section IV.C.1, and the decomposition in Eq. (15).","section":"IV.B, Eq. (12)"},{"comment":"The A-scaling conclusions for ΔE_g and ΔE_f rely on fitting the self-consistent RHB potentials Σ and Δ to Woods-Saxon forms, but the paper reports no fit parameters, residuals, or uncertainties, and no sensitivity test of the surface derivatives Σ′ and Δ′. Because the closed-form overlaps in Eq. (12) weight precisely the surface region, a small misfit in the Woods-Saxon form could change the sign or magnitude of the integrated contributions. The generalization is also drawn from a single PSS doublet, ν(2d5/2−1g7/2), along the N=50 chain. Please provide fit-quality information (e.g., parameter errors or χ²), a direct comparison of the fitted derivatives with the RHB derivatives, and at least one additional isotonic/isotopic chain or orbital pair before claiming generic scaling.","section":"IV.C"}],"minor_comments":[{"comment":"The abstract contains a typo: \"theses spin\" should read \"these spin\".","section":"Abstract"},{"comment":"The definition of Δ is inconsistent: the text defines Δ ≡ V − S, while the caption of Fig. 1 states Δ = S − V. Please unify the convention, since the sign of Δ′ enters Eqs. (8) and (12).","section":"Introduction and Fig. 1"},{"comment":"Please clarify that Δ0 denotes the integrated jump of Δ across the surface, not necessarily Δ(r=0), to make the sign in Eq. (9) unambiguous.","section":"Eq. (9)"},{"comment":"The in-figure title labels the plotted gap as ν(1g7/2 − 1g9/2), whereas the text defines the benchmark PSS splitting as ν(2d5/2 − 1g7/2); please correct the label and verify that the plotted quantity is the intended PSS gap.","section":"Fig. 2"},{"comment":"The condition W_m < 1 is not quantified for the PSS partners used in the figures; please report the numerical values of W_m for the orbitals shown in Figs. 5–7.","section":"II.C"},{"comment":"The numerical rescaling used to isolate effect ii) and to infer its A^(-1/3) scaling is not described; please specify the procedure and the fit used to extract the exponent.","section":"III.B.1"}],"recommendation":"major_revision","confidential_remarks":"The sign error in Eq. (12) is the decisive issue. I would ask the authors to evaluate ΔE_g both directly from Eq. (4) and from the corrected integration-by-parts formula, and to regenerate Fig. 5 and the related discussion. If the corrected ΔE_g(A) trend reverses, the paper's main PSS scaling claim will need substantial reworking. Please also double-check the Fig. 2 label before resubmission."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Honest take: this paper has a useful framework but one of its central formulas, Eq. (12), is printed with the wrong sign in the ΔEg term. The total first-order expression, Eq. (4), is correct, so the N=50 benchmark in Fig. 2 is not compromised. The sign error sits specifically in the decomposition that drives the paper's interpretation of why ΔEg grows with A.\n\nWhat is genuinely new here is the closed-form first-order decomposition of pseudo-spin splitting into a g/Σ part and an f/Δ part, and the identification of two distinct A-behaviors. The spin-symmetry part is well executed: the recovery of the A−2/3 scaling as a product of a norm decrease and an overlap decrease is clean and gives a mechanistic picture that goes beyond quoting the usual non-relativistic derivation. The optimization of the HO reference against the perturbation parameter follows Ref. [24] and is done properly.\n\nThe soft spot is the sign. Starting from Eq. (4), with A = Ga−Gb, the g contribution is ∫ r²(ga²−gb²)Σ dr = ∫ A' Σ dr = −∫ A Σ' dr, after dropping the ΣHO piece that vanishes by degeneracy. So the first term in Eq. (12) should be −∫(Ga−Gb)Σ' dr, not +. The f term is correctly +∫(Fa−Fb)Δ' dr because the boundary term cancels. The consequence is material: the text around Fig. 6 argues that as A increases, Σ' increasingly probes the positive lobe of Ga−Gb and therefore ΔEg increases. With the minus sign, that same overlap would make ΔEg decrease. The story inverts. Fig. 5, if computed with the printed plus sign, would not actually represent the separate g and f contributions.\n\nTwo other things are softer. The scaling conclusions for ΔEg and ΔEf rest on Woods-Saxon fits to the self-consistent potentials along a single isotonic chain, with no fit uncertainties or sensitivity test reported. That makes the A-scaling claims suggestive rather than established. Also, no code or data are provided, which limits reproducibility.\n\nThe paper is aimed at nuclear structure theorists working on shell evolution, bubble nuclei, and relativistic mean-field. It deserves a serious referee: the framework is relevant, the error is local and fixable, and the SS part alone is worth publishing. I would send it to review with a request that the referee verify the signs and ask for a robustness check on the fits.","headline":"Eq. (12) has a sign error in the ΔEg term that inverts the paper's main A-scaling interpretation; the total first-order expression is fine and the paper is salvageable.","tokens_in":12851,"tokens_out":10171,"would_cite":false,"duration_ms":93169,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["21.10.Pc","21.60.Jz"],"model":"deepseek-v4-flash","headline":"The paper derives closed-form first-order perturbative formulas for spin and pseudo-spin energy gaps in relativistic mean-field nuclear models, and shows how each gap scales with mass number.","keywords":["spin symmetry","pseudo-spin symmetry","relativistic mean-field","perturbation theory","spin-orbit splitting","nuclear shell structure","mass dependence","harmonic oscillator"],"falsifier":"Compute the first-order closed-form gaps using the actual self-consistent $\\Sigma$ and $\\Delta$ potentials without Woods-Saxon fitting for several isotonic chains and several pseudo-spin pairs; if the upper-component term no longer grows roughly linearly with $A$ or the lower-component term no longer stays roughly constant, the claimed $A$-dependence pattern is an artifact of the fits.","tokens_in":11791,"feed_emoji":"⚛️","tokens_out":13220,"duration_ms":103830,"temperature":0.7,"pith_summary":"This paper proposes a common perturbative framework for two related symmetries in atomic nuclei: spin symmetry and pseudo-spin symmetry. Working at first order around a spin-symmetric relativistic harmonic oscillator reference state, it obtains explicit closed expressions for the energy gaps between spin doublets and pseudo-spin doublets. For spin doublets, the gap depends only on the lower component of the Dirac wavefunction and on the potential difference $\\Delta = S - V$, and it reproduces the well-known $A^{-2/3}$ mass scaling as the product of two separate $A^{-1/3}$ effects. For pseudo-spin doublets, both upper and lower components and both potentials $\\Sigma = S + V$ and $\\Delta$ contribute: the upper-component term grows roughly linearly with $A$ while the lower-component term stays nearly constant. The result is a unified, parameter-light explanation of how these symmetries are broken across the nuclear chart.","feed_headline":"One formula explains nuclear spin and pseudo-spin gaps","feed_subtitle":"Around a spin-symmetric oscillator, spin gaps shrink as A^-2/3 while pseudo-spin terms separate.","key_machinery":"The central object is the spin-symmetric relativistic harmonic oscillator (SS-RHO) reference state: a Dirac Hamiltonian with harmonic-oscillator central potential $\\Sigma_{HO} = c_0 + c_2 r^2$ and constant $\\Delta_{HO} = d_0$, whose eigenfunctions $g$ (upper) and $f$ (lower) are known analytically. The perturbation operator $W = \\mathrm{diag}(\\Sigma - \\Sigma_{HO}, d_0 - \\Delta)$ is applied with the three reference parameters fixed by minimizing a perturbation parameter; in practice only the oscillator frequency is optimized. The derivation then expresses the energy gap at first order as integrals over the integrated density differences $F_i(r) = \\int_0^r dr'\\, r'^2 f_i(r')^2$ and $G_i(r) = \\int_0^r dr'\\, r'^2 g_i(r')^2$, weighted by the surface derivatives $\\Sigma'$ and $\\Delta'$. These integrated densities and their overlaps with the potential derivatives carry the whole argument about $A$-dependence.","core_discovery":"At first order in perturbation theory around the spin-symmetric relativistic harmonic oscillator, the spin-orbit splitting of a doublet collapses to $\\Delta E_{SS} \\approx (F_\\downarrow(R) - F_\\uparrow(R)) \\Delta_0$, where $F_i(r)$ is the integrated lower-component density and $\\Delta_0$ is the depth of $\\Delta = S - V$. The pseudo-spin splitting takes the analogous but two-term form $\\Delta E_{PSS} = \\int dr\\, (G_a - G_b) \\Sigma' + \\int dr\\, (F_a - F_b) \\Delta'$, where $G_i$ is the integrated upper-component density. The paper argues that these closed forms capture the dominant mechanisms: the spin gap inherits its $A^{-2/3}$ scaling from the relativistic weakening of the lower component (norm $\\sim \\omega \\sim A^{-1/3}$) and from the decreasing overlap between the surface-peaked $\\Delta'$ and the integrated density difference (another $A^{-1/3}$); the pseudo-spin gap is governed by the competition between a $\\Sigma'$-driven upper-component term that grows roughly linearly with $A$ and a $\\Delta'$-driven lower-component term that remains almost constant. This gives a common footing for spin and pseudo-spin symmetry breaking and explains the decomposition $\\Delta E_{PSS} = \\Delta E_{HO} - \\Delta E_{SO}$, where the central-potential term and the spin-orbit term have opposite signs.","pith_inferences":["Since the paper's A-dependence analysis rests on a single orbital pair in the N=50 chain, an immediate test is whether the linear-in-A growth of $\\Delta E_g$ and the constancy of $\\Delta E_f$ survive for other pseudo-spin pairs (e.g., in N=82 or N=126 isotones) and for proton doublets, where the Coulomb potential modifies $\\Sigma$ and $\\Delta$.","The framework suggests a practical way to predict shell evolution far from stability: if the surface diffusivity of $\\Sigma$ changes (e.g., in neutron-rich nuclei with thick neutron skins), the linear term $\\Delta E_g$ will grow, pushing pseudo-spin partners apart; the paper's mechanism would then predict where magic gaps weaken.","The same perturbative machinery could be carried to second order systematically (the paper shows second order shifts values but preserves trends) to build an analytic map of where pseudo-spin symmetry is accidentally restored across the nuclear chart."],"forward_implications":["If the first-order closed forms are correct, the spin-orbit splitting of any doublet can be computed directly from the integrated lower-component density and the surface slope of $\\Delta$, without solving the full Dirac equation.","The $A^{-2/3}$ scaling of spin-orbit splittings is explained as the product of two kinematic effects: the relativistic weakening of the lower component and the geometric mismatch between $F_\\downarrow - F_\\uparrow$ and $\\Delta'$, making the scaling a consequence of relativistic dynamics rather than a phenomenological input.","For pseudo-spin doublets, the near-constancy of the lower-component term and the roughly linear growth of the upper-component term imply that the total pseudo-spin splitting can change sign along a chain, which would show up as accidental (re)appearance of pseudo-spin degeneracy.","Because the perturbation is controlled ($W_m \\ll 1$) for both spin and pseudo-spin partners, the same reference state supports both calculations, so spin and pseudo-spin phenomena can be studied with one consistent expansion."],"supporting_citations":[{"why":"Introduces the spin-symmetric relativistic harmonic oscillator as a reference state and the perturbation parameter whose minimization fixes the reference parameters.","marker":"[24]"},{"why":"Provides the analytic eigenfunctions and eigenenergies of the spin-symmetric relativistic harmonic oscillator, including the norm relations used to derive the A-dependence.","marker":"[25]"},{"why":"Establishes the connection between the spin-orbit potential and the derivative of $\\Delta$, and the $A^{-2/3}$ scaling that the present first-order formula recovers.","marker":"[19]"},{"why":"Gives the relativistic derivation of the spin-orbit interaction and the condition that constant $\\Delta$ ensures spin symmetry.","marker":"[23]"},{"why":"Supplies the experimental spin-orbit splitting trend that the first-order formula reproduces and benchmarks against.","marker":"[30]"},{"why":"Previous non-perturbative analysis of pseudo-spin symmetry breaking whose divergent integrals the present closed forms avoid.","marker":"[33]"},{"why":"The relativistic Hartree-Bogoliubov code that produces the self-consistent $\\Sigma$ and $\\Delta$ potentials used in the perturbative calculation.","marker":"[27]"},{"why":"The covariant energy density functional parametrization (DD-MEV) used in the calculations.","marker":"[28]"},{"why":"Shows that pseudo-spin symmetry originates in the smallness of the $\\Sigma = S + V$ potential, the premise behind the perturbative treatment.","marker":"[9]"}],"fun_headline_variants":["Nuclear spin gaps scale as A^-2/3; pseudo-spin has two terms","Single formula captures spin and pseudo-spin symmetry breaking","Spin doublet splitting shrinks as A^-2/3 in relativistic mean-field","Perturbative expansion explains spin and pseudo-spin gaps","Spin and pseudo-spin: common footing from covariant approach"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusions about how pseudo-spin splitting changes with mass number depend on the assumption that Woods-Saxon fits to the self-consistent potentials faithfully reproduce the surface slopes of $\\Sigma$ and $\\Delta$, and that the behavior seen in one neutron orbital pair along the $N=50$ chain is generic.","fun_headline_variants_meta":{"raw":{"variants":["Nuclear spin gaps scale as A^-2/3; pseudo-spin has two terms","Single formula captures spin and pseudo-spin symmetry breaking","Spin doublet splitting shrinks as A^-2/3 in relativistic mean-field","Perturbative expansion explains spin and pseudo-spin gaps","Spin and pseudo-spin: common footing from covariant approach"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000246,"raw_usage":{"total_tokens":1575,"prompt_tokens":1016,"completion_tokens":559,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":632,"completion_tokens_details":{"reasoning_tokens":466}},"tokens_in":632,"tokens_out":559,"duration_ms":5092,"temperature":1.0,"reasoning_tokens":466,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T21:32:39.757116+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the first-order closed-form gaps using the actual self-consistent $\\Sigma$ and $\\Delta$ potentials without Woods-Saxon fitting for several isotonic chains and several pseudo-spin pairs; if the upper-component term no longer grows roughly linearly with $A$ or the lower-component term no longer stays roughly constant, the claimed $A$-dependence pattern is an artifact of the fits.","supporting_citations":[{"cited_title":"Delafosse, D","cited_arxiv_id":null,"evidence_quote":"Introduces the spin-symmetric relativistic harmonic oscillator as a reference state and the perturbation parameter whose minimization fixes the reference parameters."},{"cited_title":"Ebran, E","cited_arxiv_id":null,"evidence_quote":"Provides the analytic eigenfunctions and eigenenergies of the spin-symmetric relativistic harmonic oscillator, including the norm relations used to derive the A-dependence."},{"cited_title":"Nazarewicz, P","cited_arxiv_id":null,"evidence_quote":"Establishes the connection between the spin-orbit potential and the derivative of $\\Delta$, and the $A^{-2/3}$ scaling that the present first-order formula recovers."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the relativistic derivation of the spin-orbit interaction and the condition that constant $\\Delta$ ensures spin symmetry."},{"cited_title":"Liang, P","cited_arxiv_id":null,"evidence_quote":"Supplies the experimental spin-orbit splitting trend that the first-order formula reproduces and benchmarks against."},{"cited_title":"Liang, J","cited_arxiv_id":null,"evidence_quote":"The relativistic Hartree-Bogoliubov code that produces the self-consistent $\\Sigma$ and $\\Delta$ potentials used in the perturbative calculation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The covariant energy density functional parametrization (DD-MEV) used in the calculations."},{"cited_title":"Duguet, H","cited_arxiv_id":null,"evidence_quote":"Shows that pseudo-spin symmetry originates in the smallness of the $\\Sigma = S + V$ potential, the premise behind the perturbative treatment."}],"review_version":1}