{"id":"13fb1014-2a8d-4a58-bf42-fff82ed37579","arxiv_id":"2607.21875","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A dispersive-optical-model analysis with a pairing extension reproduces single-particle structure and (e,e'p) cross sections in 90Zr and predicts a neutron skin of 0.078 fm.","lead":"This paper extends the dispersive optical model to the open-shell nucleus 90Zr by adding a pairing gap that splits proton orbits near the Fermi energy, and uses it to describe elastic scattering, charge density, and proton-knockout cross sections. A nuclear-physics specialist would read it to see whether the DOM framework can be pushed from closed-shell nuclei to semi-magic open-shell systems and to weigh its predicted neutron skin.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The (e,e'p) validation is partly circular: the 1p3/2 and 0f5/2 spectroscopic factors are renormalized to the measured strengths via Eq. (24), leaving the 1p1/2 transition as the only independent test of the pairing extension.","rationale":"The central claim is that the DOM+pairing yields a good description of the (e,e'p) cross sections. This claim is supported by the agreement in Figs. 10-13. However, the comparison is not a clean prediction for the 1p3/2 and 0f5/2 transitions. The paper explicitly renormalizes the DOM spectroscopic factors using the experimental fragmentation (Eq. 24) and an additional theoretical ratio. The final strengths 0.577 and 0.522 are therefore not outputs of the model alone. This does not make the shape comparison meaningless—the distorted-wave and overlap shapes are still predicted—but it removes the magnitude as a test. The only unrenormalized tests are the 1p1/2 ground-state transition (0.396 vs 0.360, within the ±0.05 error) and the small 0g9/2 fragment (0.047 vs 0.054). The 0g9/2 cross sections have large relative uncertainty. Thus the pairing extension's predictive power for spectroscopic factors rests on essentially one transition. The reader's concern about the BCS ansatz is related but secondary: if the single-gap ansatz were wrong, the 1p1/2 split would likely fail, but it does not; the more immediate logical flaw is that the model's (e,e'p) validation is partly circular. The proposed test—removing the renormalizations and comparing the unrenormalized cross sections to the data—would directly show whether the model can predict these magnitudes. If they overshoot, the abstract's phrasing 'yield a good description' is too strong.","tokens_in":20420,"tokens_out":8933,"duration_ms":80250,"concrete_test":"Recompute the 1p3/2 and 0f5/2 momentum distributions using only the first-step DOM spectroscopic factors split by the pairing formulas (Eqs. 14-15), i.e., omit the Eq. (24) experimental renormalization and the spectral-strength ratio correction. Compare the resulting cross sections to the Nikhef data at both Ep=70 and 100 MeV (Figs. 10b, 11b, 12a, 13a). If the integrated cross sections exceed the data by more than the combined statistical and systematic uncertainties (e.g., >20%, given Z≈0.711 vs 0.465), the model does not by itself describe these transitions and the abstract's claim should be qualified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the DOM+pairing framework yields a good description of the 90Zr(e,e'p) cross sections is weakened because for two of the four transitions (1p3/2 and 0f5/2) the spectroscopic factors are not predicted by the model but renormalized to the experimental data. In Sec. IV, Eq. (24) enforces that the integrated DOM strength in the fragment region equals the experimental strength, and an additional correction is applied using the ratio of the DOM spectral strength integrated from -20 MeV to the Fermi energy to the normalization of the Eq. (6) quasihole solution. The quoted DOM spectroscopic factors 0.577 (1p3/2) and 0.522 (0f5/2) therefore incorporate these ad hoc renormalizations plus a small pairing correction. Consequently, the agreement in Figs. 10(b), 11(b), 12(a), and 13(a) tests only the shape of the momentum distributions, not the model's ability to predict the cross-section magnitudes. The only unrenormalized tests of the pairing extension are the 1p1/2 ground-state transition (0.396 vs 0.360, within the ±0.05 error) and the weak 0g9/2 fragment (0.047 vs 0.054). The latter cross section has large relative uncertainty. If the unrenormalized DOM+pairing strengths for 1p3/2 and 0f5/2 (≈0.711 and ≈0.670 from Eq. (6), before corrections) significantly exceed the Nikhef values, the central claim overstates the predictive power of the model.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends the nonlocal dispersive optical model (DOM) to the semi-magic nucleus 90Zr by adding a BCS-like pairing treatment for the open proton shell. The DOM self-energy is fitted to proton and neutron elastic scattering, reaction and total cross sections, single-particle energies, particle numbers, charge density, and binding energy. The pairing extension splits the strength of the 1p1/2, 1p3/2, 0f5/2, and 0g9/2 orbits across the Fermi energy using one global gap parameter Δ = 1 MeV. The resulting overlaps and spectroscopic factors are used in DWIA calculations of 90Zr(e,e'p)89Y momentum distributions and compared to Nikhef data. A neutron skin R_skin = 0.078 ± 0.039 fm is predicted.","tokens_in":20927,"tokens_out":5758,"duration_ms":58604,"significance":"If the central claim holds, this is a valuable step: the DOM framework, previously applied to doubly magic nuclei, is extended to a semi-magic open-shell nucleus while maintaining a simultaneous description of scattering, charge density, single-particle strength, and knockout cross sections. The paper is transparent about its parametrization (Appendix A), provides extensive comparisons to scattering and density data, and gives a falsifiable neutron-skin prediction. The most important independent evidence is the 1p1/2 ground-state transition (0.396 vs. 0.360, within the claimed uncertainty) and the weak 0g9/2 fragment (0.047 vs. 0.054). However, as detailed below, the predictive content of the (e,e'p) comparisons is substantially weakened because two of the four transitions are renormalized to the experimental integrated strength via Eq. (24). The framework is promising, but the paper's central claim overstates what is actually predicted.","major_comments":[{"comment":"Eq. (24) explicitly enforces equality between the DOM integrated spectral strength and the experimental integrated strength in the fragment region for the 1p3/2 and 0f5/2 orbits. The quoted DOM spectroscopic factors in Table II (0.577 and 0.522) therefore inherit this experimental normalization, together with an additional ad hoc 'theoretical correction' described in the text. Consequently, Figs. 10(b), 11(b), 12(a), and 13(a) test only the shape of the momentum distributions, not the model's ability to predict cross-section magnitudes for those transitions. The abstract's claim that the model yields 'a good description of the (e,e'p) cross sections' should be tempered, or the unrenormalized DOM values (≈0.711 and ≈0.670 from Eq. (6)) should be reported so the reader can judge the predictive content. As written, the central claim is partly circular for two of the four transitions.","section":"Sec. IV, Eq. (24)"},{"comment":"The pairing extension relies on a single global pairing gap Δ = 1 MeV chosen by hand ('close to the empirical value') and not determined by the fit. The uncertainty bands in Figs. 10–13 are obtained by varying Δ = 1 ± 0.1169 MeV, which is a sensitivity test, not a χ²-constrained uncertainty. Because this parameter directly controls the splitting in Eqs. (11)–(17) and hence the spectroscopic factors used in the (e,e'p) calculation, the paper should state explicitly that the pairing treatment introduces one free parameter and should quantify how much of the observed agreement for the independent 1p1/2 and 0g9/2 transitions depends on this choice. A comparison with the empirical estimate from Eq. (13) would also help.","section":"Sec. II B and Sec. III"},{"comment":"The 'theoretical correction' applied to the 1p3/2 and 0f5/2 strengths is an ad hoc scaling: the DOM spectral strength integrated from -20 MeV to the Fermi energy (0.797 and 0.757) is divided by the quasihole normalization from Eq. (6) (0.711 and 0.670). The paper does not derive this correction from a controlled many-body approximation or from an independent observable. If both Eq. (24) and this correction were omitted, the unrenormalized DOM strengths are about 0.711 and 0.670, roughly 50% larger than the Nikhef values 0.465 and 0.462. The manuscript should report these unrenormalized values and discuss whether the discrepancy indicates a genuine deficiency of the DOM imaginary part or of the simple BCS pairing ansatz, rather than absorbing it into normalization.","section":"Sec. IV, paragraph after Eq. (24)"}],"minor_comments":[{"comment":"The fragment region FR is not defined explicitly. Please specify the integration limits used in Eq. (24), e.g., the energy interval around each quasihole peak corresponding to the experimentally identified fragments.","section":"Sec. IV, Eq. (24)"},{"comment":"The row for ρwb(p,n) lists the unit as [MeV]; this should presumably be [fm], since it is a radial parameter in Eq. (A4). Please correct.","section":"Table IV"},{"comment":"There appears to be a typo: 'aHF aym' should read 'aHF asy'.","section":"Appendix A, Eq. (A3)"},{"comment":"The column header 'Znlj DOM' does not indicate that the 1p3/2 and 0f5/2 entries include the Eq. (24) renormalization and the additional theoretical correction. A more precise header or a footnote would prevent misinterpretation.","section":"Sec. IV, Table II"},{"comment":"The four columns are labeled 'Theory', 'Exp', 'Exp Split', and 'Theory Split' in the caption, but the text refers to 'the second column', 'the third column', and 'the fourth column'. Please make the column numbering explicit and consistent.","section":"Sec. III, Fig. 6"},{"comment":"The notation Sn−lj(E) in Eq. (18) and Slj(E) in Eq. (5) is visually similar; consider a clearer distinction between the spectral strength of a quasihole state and the angle-integrated spectral function.","section":"Sec. II A, Eq. (6)"}],"recommendation":"major_revision","confidential_remarks":"The circularity concern raised by Eq. (24) is real and load-bearing: two of the four (e,e'p) transitions are renormalized to experiment, so the central claim is stronger than the evidence. However, the paper itself discloses the procedure, and there is genuine independent content in the 1p1/2 and 0g9/2 comparisons plus a substantial body of scattering and density fits. I do not recommend rejection; rather, the authors should rewrite the abstract/conclusion to state which transitions are predictive and which are shape-only, report the unrenormalized strengths, and quantify the role of the hand-picked pairing gap. With those revisions, the paper could meet the journal's standard."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a solid extension of the nonlocal DOM framework to a semi-magic, open-shell nucleus (90Zr) by tacking a BCS-like pairing step onto the usual DOM self-energy. The new pieces are the open-shell proton treatment, the two-step pairing implementation, and the first DOM-based neutron-skin prediction for 90Zr. The fits to elastic scattering, total/reaction cross sections, charge density, and elastic electron scattering are extensive and convincing. The paper also does something rare: it clearly states its own limitations (the ad hoc gap, the fragmentation correction, the 1.5% binding-energy miss) and gives enough detail to see exactly what was done.\n\nThe soft spots are real but not fatal. The stress-test note holds up: for the 1p3/2 and 0f5/2 fragments, Eq. (24) renormalizes the DOM strength to the experimental integrated strength, so the good agreement in those momentum distributions tests shapes, not absolute magnitudes. The paper is transparent about this, but the abstract's claim to \"a good description\" overstates the predictive content for those two transitions. The independent tests are the 1p1/2 ground-state transition (0.396 vs 0.360) and the weak 0g9/2 fragment (0.047 vs 0.054); those are encouraging, and the 1p1/2 comparison is the strongest single check of the pairing split. The pairing gap itself is chosen by hand at 1 MeV, with a sensitivity band, rather than fitted or derived from Eq. (13); that is a minor weakness because the gap mostly shifts where the splitting sits and the data constrain it only loosely.\n\nThe citation pattern is fair and the work builds honestly on the authors' prior DOM results and on Migdal's pairing formalism. The binding energy being off by 1.5% is a known DOM trait and is not load-bearing for the conclusions, though it would be nice to see it discussed more.\n\nOverall: this is a useful capability demonstration, not a paradigm shift. The framework is credible, the methodology is reproducible in principle, and the neutron-skin prediction is a falsifiable target. The circularity concern weakens the central claim but does not sink it. I would take this seriously in peer review and would ask the authors to either compute unrenormalized strengths for the 1p3/2 and 0f5/2 transitions or label them explicitly as constrained rather than predicted.\n\nRecommendation: send to peer review. The paper deserves referee time, and the authors have the tools and track record to address the normalization question directly.","headline":"Extends DOM to an open-shell nucleus with a simple pairing ansatz, but the two strongest (e,e'p) tests are renormalized to experiment; the 1p1/2 fragment is the real independent check.","tokens_in":21425,"tokens_out":1108,"would_cite":true,"duration_ms":13743,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["24.10.Ht","25.30.Fj","21.10.Jx"],"model":"deepseek-v4-flash","headline":"Adding a single pairing gap to the dispersive optical model reproduces the 90Zr(e,e'p) cross sections and yields a neutron-skin prediction of 0.078 ± 0.039 fm.","keywords":["dispersive optical model","pairing correlations","semi-magic nucleus","90Zr","spectroscopic factors","(e,e'p) knockout","single-particle structure","neutron skin"],"falsifier":"Measure the 90Zr(e,e'p) cross sections with higher missing-energy resolution and extract spectroscopic factors for the 1p1/2, 1p3/2, 0f5/2, and 0g9/2 fragments; if the extracted factors deviate from the predicted values (0.396, 0.577, 0.522, and 0.047, respectively) by more than the stated uncertainties, the single-gap pairing ansatz is ruled out. Alternatively, a parity-violating electron-scattering measurement of the 90Zr neutron skin would test the predicted R_skin = 0.078 ± 0.039 fm; a skin outside this range would falsify the model's ground-state densities.","tokens_in":20302,"feed_emoji":"⚛️","tokens_out":5774,"duration_ms":53125,"temperature":0.7,"pith_summary":"This paper extends the dispersive optical model—a framework that fits elastic scattering, reaction cross sections, and ground-state properties with one energy-dependent self-energy—to the semi-magic nucleus 90Zr, whose 40 protons form an open shell. The authors show that adding a BCS-like pairing step with a single gap parameter of about 1 MeV splits the strength of the 1p1/2, 1p3/2, 0f5/2, and 0g9/2 orbits across the Fermi energy, producing the smooth mostly-full-to-mostly-empty occupation pattern expected in a paired system. With this extension, the model reproduces the measured (e,e'p) momentum distributions for four discrete transitions, improves the description of the charge density, and predicts a neutron skin of 0.078 ± 0.039 fm. If correct, this gives a unified description of scattering, structure, and knockout for an open-shell nucleus and opens a path to studying single-particle strength in semi-magic systems.","feed_headline":"One pairing gap reproduces 90Zr knockout data","feed_subtitle":"Extended dispersive optical model matches (e,e'p) cross sections and predicts a 0.078 fm neutron skin.","key_machinery":"The central object is the irreducible self-energy of the nonlocal dispersive optical model, a complex one-body potential whose energy dependence is generated by a subtracted dispersion relation from a fitted imaginary part. The new mechanism is a two-step pairing treatment: first solve the Dyson equation with the DOM self-energy to get quasihole states and spectroscopic factors; then apply a BCS-like splitting using a single gap parameter of about 1 MeV and chemical potentials from experimental binding energies, redistributing each orbit's strength into upper and lower fragments. These split spectroscopic factors, combined with DOM overlap functions and distorted waves in a distorted-wave im","core_discovery":"The central claim is that the nonlocal dispersive optical model, previously validated for doubly closed-shell nuclei like 40Ca, 48Ca, and 208Pb, can be extended to the semi-magic nucleus 90Zr by treating the open proton shell with pairing correlations in a two-step procedure. In the first step the DOM self-energy is fitted to elastic scattering, total and reaction cross sections, single-particle energies, particle numbers, charge density, and binding energy. In the second step, a BCS-like pairing ansatz with a single gap parameter of about 1 MeV splits the spectroscopic strength of the four orbits nearest the Fermi energy (1p1/2, 1p3/2, 0f5/2, and 0g9/2) into fragments above and below the Fe","pith_inferences":["If the BCS-like pairing treatment holds for heavier semi-magic nuclei, the same two-step DOM could be applied to nuclei like 92Mo or 96Ru, where proton pairing is expected to be stronger, providing a systematic map of single-particle strength.","The neutron-skin prediction of 0.078 fm for 90Zr could be tested with parity-violating electron scattering, offering an asymmetry benchmark between the 48Ca and 208Pb measurements already reported.","The smooth occupation transition predicted here suggests that higher-resolution (e,e'p) experiments could resolve the small predicted fragments (for example 0.047 for 0g9/2) and discriminate the BCS ansatz from phonon-coupling models.","A natural extension is to treat the gap parameter as orbit-dependent or compute it self-consistently instead of fixing it at 1 MeV; testing the sensitivity of the (e,e'p) cross sections to the gap would reveal how much physics the single parameter is carrying."],"forward_implications":["The DOM with pairing reproduces the 90Zr(e,e'p)89Y momentum distributions for the 1p1/2 ground-state transition and the 1p3/2, 0f5/2, and 0g9/2 excited transitions, with spectroscopic factors close to those extracted from experiment.","Including pairing improves the calculated charge density and elastic electron-scattering cross sections, extending the DOM's density description from closed-shell to open-shell nuclei.","The model yields a neutron-skin prediction of 0.078 ± 0.039 fm for 90Zr, a benchmark between 48Ca and 208Pb.","A single gap parameter near the empirical pairing gap suffices to describe the fragmentation of four proton orbits, suggesting a simple route to other semi-magic nuclei.","The occupation numbers of proton orbits show a smooth mostly-full-to-mostly-empty transition characteristic of pairing, distinct from the abrupt transition in closed-shell nuclei like 40Ca."],"fun_headline_variants":["Single gap extends optical model to 90Zr","90Zr's open shell tamed by one pairing gap","Pairing gap unlocks 90Zr knockout data","Nonlocal DOM with one gap fits 90Zr","One pairing gap explains 90Zr's shell structure"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The analysis assumes that a single energy-gap parameter, the same for all four proton orbits near the Fermi energy, is enough to describe how the open proton shell splits into occupied and empty fragments; if the real nucleus distributes this strength through orbit-dependent or more complicated correlations, the predicted spectroscopic factors and knockout cross sections would not hold.","fun_headline_variants_meta":{"raw":{"variants":["Single gap extends optical model to 90Zr","90Zr's open shell tamed by one pairing gap","Pairing gap unlocks 90Zr knockout data","Nonlocal DOM with one gap fits 90Zr","One pairing gap explains 90Zr's shell structure"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000752,"raw_usage":{"total_tokens":3151,"prompt_tokens":683,"completion_tokens":2468,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":427,"completion_tokens_details":{"reasoning_tokens":2393}},"tokens_in":427,"tokens_out":2468,"duration_ms":16614,"temperature":1.0,"reasoning_tokens":2393,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T06:25:09.906949+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the 90Zr(e,e'p) cross sections with higher missing-energy resolution and extract spectroscopic factors for the 1p1/2, 1p3/2, 0f5/2, and 0g9/2 fragments; if the extracted factors deviate from the predicted values (0.396, 0.577, 0.522, and 0.047, respectively) by more than the stated uncertainties, the single-gap pairing ansatz is ruled out. Alternatively, a parity-violating electron-scattering measurement of the 90Zr neutron skin would test the predicted R_skin = 0.078 ± 0.039 fm; a skin outside this range would falsify the model's ground-state densities.","supporting_citations":[],"review_version":1}