{"id":"9e9f6683-ba62-46ec-88b2-9c2d378f3b62","arxiv_id":"1908.03561","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The paper implements a separable pairing force in the cranking CDFT shell-model-like approach and reproduces the rotational bands of 60Fe.","lead":"A new nuclear modeling code combines a finite-range separable pairing force with a cranking relativistic density functional approach that conserves particle number exactly. Applied to three rotational bands in iron-60, it reproduces the measured energies and shows better convergence with model space size than the older monopole-pairing version.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"MPC-convergence check is quoted at a single low frequency; spectra and crossing at hbar-omega ≈ 0.85 MeV may be unconverged at dimension 1000.","rationale":"The reader's weakest assumption -- transferability of G and a from nuclear-matter Gogny fits to finite 60Fe -- is a fair physical question, and a refit or a comparison with Gogny-HFB in 60Fe would address it. But it does not directly threaten the internal consistency of the method: even if the parameters were refitted, the algorithm would still be valid and its convergence property would likely survive. The more load-bearing check for the paper's own headline is numerical convergence of the observables that are actually compared with data. The current convergence statement is limited to one frequency where the band structure is still simple; the crossing region is exactly where MPC truncation is most dangerous. The paper itself argues that monopole pairing is unreliable because of strong dimension dependence, so the analogous question for the separable force near the crossing must be settled before 'better convergence' and 'excellent agreement' can be accepted. This is a concrete and finite computational task. The CONDITIONAL verdict is therefore unchanged, with the condition sharpened to require a full spectrum-level convergence test.","tokens_in":14165,"tokens_out":14207,"duration_ms":156540,"concrete_test":"Recompute the three bands with MPC dimensions 2000 (and 3000 if feasible) for both protons and neutrons, keeping G = -728 MeV fm^3 and a = 0.644 fm; compare the yrast energies, I(omega) relations, and band-crossing frequencies with the published 1000-configuration results. If the crossing frequency in band A shifts by more than about 50 keV, or if level energies shift by more than about 100 keV, the 1000-configuration results are not sufficiently converged and the 'excellent agreement' claim cannot be evaluated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of better convergence and excellent agreement rests on the MPC-space truncation at 1000 configurations for neutrons and protons. The convergence test in Sec. III quotes changes in total energy and alignment at hbar-omega = 0.2 MeV only (0.151% and 0.667% between 1000 and 2000). The physics of interest, however, lies near the band crossing in band A at hbar-omega ≈ 0.85 MeV (Sec. IV.B), where configuration mixing is maximal and truncation errors are expected to be largest. No convergence of the actual spectra, yrast energies, or crossing frequencies with MPC dimension is reported. If the 1000-configuration results are not converged in that region, the 'delayed' band crossing and the claimed agreement with the 60Fe data would be truncation artifacts rather than properties of the separable-pairing SLAP.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops the shell-model-like approach (SLAP) based on cranking covariant density functional theory (CDFT) with a separable pairing force of the Tian-Ma-Ring form, replacing the monopole pairing force used in the earlier cranking CDFT-SLAP of Ref. [42]. The method is applied to 60Fe to study the positive-parity yrast band A and the two negative-parity signature partner bands B and C. The calculated total energies, angular momenta, pairing energies, and deformation parameters are compared with data from Ref. [68] and with the monopole-pairing and no-pairing results of Ref. [42]. The paper reports better MPC-dimension convergence than the monopole variant and claims excellent agreement with the data.","tokens_in":14355,"tokens_out":3856,"duration_ms":39363,"significance":"The implementation is a genuine technical step: the separable pairing force is finite-range, avoids ultraviolet divergences, and is computationally cheaper than other finite-range forces, and the Appendix provides a detailed derivation of the pairing matrix elements in the 3D HO basis. The calculation is predictive in that the pairing parameters G and a are taken from a nuclear-matter fit rather than fitted to 60Fe data, and the comparison with the monopole-pairing variant is a useful benchmark. However, the quantitative evidence for the central claims is thin: convergence is tested only at one low frequency, the 'excellent agreement' is asserted without error measures, and the band-A crossing is admitted to be delayed. If these points are addressed, the paper would be a valid contribution.","major_comments":[{"comment":"The MPC-dimension convergence test is reported only for ℏω = 0.2 MeV, where the changes between dimension 1000 and 2000 are 0.151% for the total energy and 0.667% for the alignment. The physics of interest, however, lies near the band crossing in band A at ℏω ≈ 0.85 MeV (Sec. IV.B), where configuration mixing is maximal and truncation errors are expected to be largest. Please report the same convergence data at frequencies across the crossing region, including total energies, alignments, and the crossing frequency itself as functions of MPC dimension. Without this, the delayed crossing and the claimed agreement could be truncation artifacts.","section":"§III, Fig. 1"},{"comment":"The statement 'excellent agreement with the available data' is not supported by any quantitative error measure. No rms deviations or average differences for excitation energies, transition energies, or alignments are given. Please include a table comparing calculated and experimental level energies for bands A, B, and C with deviations. The text also concedes that the band crossing in band A occurs 'a bit later' than the data; quantify this delay and discuss whether it lies within the typical accuracy of the approach.","section":"§IV.A, Fig. 2"},{"comment":"For band A, the crossing frequency from the separable-pairing calculation is ℏω ≈ 0.85 MeV, versus ≈ 0.75 MeV from both the monopole-pairing and no-pairing calculations of Ref. [42]. Since the data indicate a structural change around I = 8ℏ, the separable force produces a larger crossing delay than the simpler monopole force. Please state explicitly whether this is an improvement or a known deficiency, and test whether the delay is sensitive to the MPC truncation (see the first major comment).","section":"§IV.B, Fig. 3"},{"comment":"The separable-pairing parameters G = −728 MeV fm^3 and a = 0.644 fm are taken from a nuclear-matter fit to Gogny pairing gaps in Ref. [65]. The paper provides no evidence that these parameters are transferable to the finite nucleus 60Fe when combined with the PC-PK1 functional and the SLAP diagonalization. Please add a sensitivity study (for example, vary G and a by a few percent and show the effect on yrast energies and crossing frequencies) or cite previous finite-nucleus applications that validate this parameter set in the same framework.","section":"§II.B and §III"}],"minor_comments":[{"comment":"The affiliation line contains visible spacing artifacts ('den sity', 'Chin a') that should be corrected.","section":"Title page"},{"comment":"The experimental data are shown as solid dots in all panels, which makes it hard to distinguish bands A, B, and C at a glance; consider using distinct point styles for the data sets.","section":"Figs. 2 and 3"},{"comment":"The comparison of neutron pairing energies from the separable and monopole forces is informative, but the text could clarify that these are model-dependent quantities rather than directly observable pairing condensates.","section":"§IV.C"},{"comment":"Ref. [42] is cited for both the monopole-pairing SLAP results and the no-pairing results; it would be helpful to state explicitly in the text which curves are reproduced from that paper.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is a modest but sound extension of Ref. [42]. The main risk is that the claims of 'excellent agreement' and 'better convergence' are stronger than the presented evidence; the requested convergence checks and quantitative comparisons are straightforward and should be within the scope of a revision. The single-author format and citation pattern are typical for this research group."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does what it says: it implements a separable pairing force in the cranking CDFT-SLAP framework, derives the matrix elements explicitly, and shows that the 60Fe spectra can be described without fitting any parameter to this nucleus. The genuinely new piece is the combination — prior SLAP work used monopole pairing, and separable pairing had only been used in TAC-CDFT, not in SLAP. That is a useful step for the subfield.\n\nWhat the paper does well: the derivation in the appendix is explicit and the separable form makes the MPC-space convergence much better than the monopole case, where the effective strength had to be re-tuned with dimension. That is a real practical improvement. The comparison with data for the three bands is honest: the band-A crossing is noted as a bit late, and the negative-parity bands look reasonable. No code or data are shipped, but for this kind of method paper that is common.\n\nSoft spots, in proportion: the convergence check is quoted at ħω = 0.2 MeV, where changes between MPC dimensions 1000 and 2000 are 0.151% and 0.667%. But the physics they care about — the band crossing in band A — sits at ħω ≈ 0.85 MeV, where configuration mixing is maximal. There is no convergence statement for energies, alignments, or the crossing frequency itself in that region. That is a genuine gap. If the 1000-dimension results shift at higher frequency, the delayed crossing could be partly a truncation artifact. The paper should show the crossing frequency as a function of MPC dimension, or at least quote the dimension dependence at the frequencies where the bands actually cross.\n\nThe “excellent agreement” claim is also stronger than what the comparison supports. The fits are eyeballed, not quantified. A few rms deviations in level spacings or I(ω) would let the reader judge whether separable pairing actually beats monopole pairing. Right now the evidence is suggestive, not conclusive.\n\nOn the transferability of G and a from nuclear matter to 60Fe: that is a real assumption, but it is not a fitted parameter, so the circularity burden is low. Still, the paper could acknowledge explicitly that the strength was fitted to Gogny pairing gaps in nuclear matter, not to this nucleus.\n\nBottom line: a workmanlike extension with a clean derivation and a fair application. The central argument holds — separable pairing in SLAP is feasible and gives better convergence than the monopole version. The missing convergence check near the crossing and the lack of quantitative error measures are fixable in revision. I would send it to a referee; a specialist will get real value from it, especially the matrix-element derivation and the convergence comparison. I would not cite it in my own work, but that is because I do not work in this subfield.\n\nRecommendation: accept for peer review, with the expectation that the referee asks for the high-frequency convergence check and quantitative comparisons with data.","headline":"A solid, honest method paper that adds a realistic separable pairing force to cranking CDFT-SLAP and tests it on one nucleus; the main soft spot is that convergence is only demonstrated at low rotational frequency, not where the band crossing happens.","tokens_in":14840,"tokens_out":2089,"would_cite":false,"duration_ms":22260,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["21.10.-k","21.60.Cs","21.60.Jz","27.50.+e"],"model":"deepseek-v4-flash","headline":"The cranking shell-model-like approach (SLAP) is extended to a finite-range separable pairing force and reproduces all three observed rotational bands of 60Fe with cleaner convergence than the monopole-pairing version.","keywords":["shell-model-like approach","cranking covariant density functional theory","separable pairing force","pairing correlations","rotational bands","60Fe","many-particle configuration space","band crossing"],"falsifier":"Apply the same formalism and the same G and a to a neighbouring even-mass Fe isotope such as 58Fe or 62Fe and compare the calculated yrast energies and band-crossing frequencies with measured data; a systematic growing offset, or disagreement with measured odd-even mass differences, would show that the nuclear-matter-fitted parameters do not transfer to finite nuclei.","tokens_in":13986,"feed_emoji":"⚛️","tokens_out":16826,"duration_ms":143860,"temperature":0.7,"pith_summary":"Rotating nuclei are usually modelled with pairing treated by BCS-like approximations that do not conserve particle number and can artificially quench pairing at high spin. This paper extends the cranking shell-model-like approach (SLAP), which diagonalizes the rotating many-body Hamiltonian in a truncated many-particle configuration space and conserves particle number exactly, by replacing the monopole pairing force with a finite-range separable pairing force from Ref. [65]. Applied to 60Fe, the method reproduces the positive-parity yrast (lowest-energy) band and the two negative-parity signature partner bands, matching available data at least as well as the monopole version and improving the negative-parity bandheads. The separable force also converges more cleanly as the configuration space grows, so the effective pairing strength no longer has to be re-fitted when the truncation is enlarged; that is what makes the calculation more predictive rather than phenomenological.","feed_headline":"Separable pairing force reproduces all three 60Fe bands","feed_subtitle":"Exact particle-number diagonalization now uses a finite-range force and matches all observed 60Fe bands.","key_machinery":"SLAP is the shell-model-like approach: instead of a BCS quasiparticle vacuum, the cranking many-body Hamiltonian $\\hat{H}=\\hat{H}'+\\hat{H}_{\\rm pair}$ is diagonalized exactly in a truncated many-particle configuration (MPC) space built from cranked single-particle Routhians (single-particle energies in the rotating frame), so particle number is conserved and blocking effects are handled automatically. The pair interaction is the separable finite-range force of Ref. [65], $\\hat{V}_{\\rm pair}=G\\,\\delta(\\mathbf{R}-\\mathbf{R}')\\,P(\\mathbf{r})P(\\mathbf{r}')\\,\\frac{1}{2}(1-P^{\\sigma})$, with $P(\\mathbf{r})=e^{-r^{2}/4a^{2}}/(4\\pi a^{2})^{3/2}$. In the three-dimensional harmonic-oscillator basis its antisymmetrized matrix elements factorize into a sum of separable terms built from one-dimensional Talmi-Moshinsky brackets, which makes the many-body diagonalization practical and gives the cleaner convergence with MPC dimension. The finite Gaussian range avoids the ultraviolet divergence of zero-range forces and, because the force is not restricted to $J=0$ pairs, the pairing energy stays finite after band crossing; occupation probabilities from the diagonalization are fed back into the CDFT densities and currents for self-consistency.","core_discovery":"Using the PC-PK1 point-coupling density functional in the particle-hole channel and the separable pairing parameters $G = -728$ MeV fm$^3$ and $a = 0.644$ fm from Ref. [65] in the particle-particle channel, the cranking CDFT-SLAP reproduces the three observed rotational bands of 60Fe: band A (positive-parity yrast) and bands B and C (negative-parity signature partners). For band A the agreement is comparable to the monopole-pairing calculation and clearly better than ignoring pairing, with a band crossing that appears slightly later, about $\\hbar\\omega = 0.85$ MeV instead of $0.75$ MeV. For bands B and C the separable force gives a better bandhead description. The calculations also show smoother convergence with the dimension of the many-particle configuration space, and pairing energies that fall gradually with rotation without collapsing to zero, because the separable force carries correlations in pairs beyond $J = 0$.","pith_inferences":["Editorial inference: if the clean convergence with MPC dimension holds in other nuclei, the method turns the pairing strength into a fixed input and enables systematic large-scale surveys of rotating nuclei without truncation-dependent re-normalization.","Editorial inference: a sharper test of the separable force would be to compute odd-even mass differences or two-particle transfer strengths in odd and even Fe isotopes; the force's $J>0$ components should leave fingerprints beyond level energies.","Editorial inference: the slightly later band-crossing frequency in band A may be a systematic tendency of this nuclear-matter-fitted force; comparing 58Fe and 62Fe with the same parameters would show whether the offset grows in a predictable way with mass."],"forward_implications":["For the positive-parity band A of 60Fe, the separable-pairing calculation places the band crossing at about $\\hbar\\omega=0.85$ MeV, slightly above the monopole and no-pairing values near $0.75$ MeV, and still reproduces the observed sudden structure change at spin $I=8\\hbar$.","For the two negative-parity signature partner bands, the separable force yields a better description of the bandheads, and band B's crossing at $\\hbar\\omega\\approx1.1$ MeV is reproduced well.","Total energy and angular momentum converge as the many-particle configuration space grows; enlarging the MPC dimension from 1000 to 2000 changes the total energy by only about 0.15 percent at $\\hbar\\omega=0.2$ MeV, whereas the monopole version needs its effective pairing strength re-tuned with dimension.","Neutron and proton pairing energies decrease smoothly with rotational frequency rather than collapsing, and for band A the separable-force neutron pairing energy stays near 1 MeV after crossing while the monopole result drops to almost zero.","Because the separable force includes pair correlations with angular momentum beyond $J=0$, the calculated quadrupole deformation is slightly larger than in the monopole calculation even where the pairing energy is larger."],"supporting_citations":[{"why":"Introduces the separable finite-range pairing force used here and determines its parameters G and a from nuclear-matter pairing gaps.","marker":"[65]"},{"why":"Develops the cranking CDFT-SLAP with monopole pairing and supplies the 60Fe results, shape evolutions, and band-crossing picture that this paper extends and compares against.","marker":"[42]"},{"why":"Provides the experimental data for the three 60Fe rotational bands used as the comparison benchmark.","marker":"[68]"},{"why":"Supplies the PC-PK1 point-coupling density functional used in the particle-hole channel.","marker":"[71]"},{"why":"Introduces the SLAP diagonalization within CDFT and the feedback of occupation probabilities into densities and currents.","marker":"[41]"},{"why":"Supplies the TAC-CDFT calculation used to validate the cranking code when pairing is switched off.","marker":"[23]"},{"why":"Gives the harmonic-oscillator formalism for separable pairing matrix elements that the appendix follows.","marker":"[80]"}],"fun_headline_variants":["Separable pairing fits all three 60Fe bands","Finite-range pairing matches 60Fe spectra","Cranking CDFT-SLAP with separable force fits 60Fe","Three 60Fe bands reproduced by separable pairing"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the separable pairing parameters fitted to nuclear-matter pairing gaps, $G = -728$ MeV fm$^3$ and $a = 0.644$ fm, retain their validity inside the finite nucleus 60Fe when used with the PC-PK1 functional and the SLAP diagonalization.","fun_headline_variants_meta":{"raw":{"variants":["Separable pairing fits all three 60Fe bands","Finite-range pairing matches 60Fe spectra","Cranking CDFT-SLAP with separable force fits 60Fe","Three 60Fe bands reproduced by separable pairing"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001261,"raw_usage":{"total_tokens":5104,"prompt_tokens":826,"completion_tokens":4278,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":442,"completion_tokens_details":{"reasoning_tokens":4212}},"tokens_in":442,"tokens_out":4278,"duration_ms":27920,"temperature":1.0,"reasoning_tokens":4212,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:27:47.632943+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Apply the same formalism and the same G and a to a neighbouring even-mass Fe isotope such as 58Fe or 62Fe and compare the calculated yrast energies and band-crossing frequencies with measured data; a systematic growing offset, or disagreement with measured odd-even mass differences, would show that the nuclear-matter-fitted parameters do not transfer to finite nuclei.","supporting_citations":[{"cited_title":"Molique and J","cited_arxiv_id":null,"evidence_quote":"Introduces the separable finite-range pairing force used here and determines its parameters G and a from nuclear-matter pairing gaps."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the experimental data for the three 60Fe rotational bands used as the comparison benchmark."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the PC-PK1 point-coupling density functional used in the particle-hole channel."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the SLAP diagonalization within CDFT and the feedback of occupation probabilities into densities and currents."},{"cited_title":"Meng and S.-G","cited_arxiv_id":null,"evidence_quote":"Supplies the TAC-CDFT calculation used to validate the cranking code when pairing is switched off."},{"cited_title":"Koepf and P","cited_arxiv_id":null,"evidence_quote":"Gives the harmonic-oscillator formalism for separable pairing matrix elements that the appendix follows."}],"review_version":1}