{"id":"6a43d327-2086-4d18-91ca-26d50deada01","arxiv_id":"2507.23244","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"HFB+QRPA calculations for 82-98Mo show low-lying dipole strength correlates with skin thickness, with skin-oscillation states exhibiting moderate two-quasiparticle fragmentation but limited coherence.","lead":"This paper studies low-energy dipole excitations in molybdenum isotopes using a self-consistent HFB+QRPA model, finding that the strength of these excitations correlates with the development of neutron or proton skins. It offers a microscopic explanation of the so-called pygmy dipole resonance in terms of skin oscillations and moderate collectivity, distinct from the giant dipole resonance.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The collectivity quantifier N* (Eq. 13) depends explicitly on basis size through 1/N2qp, yet QRPA convergence is only checked for HFB energies; the moderate-collectivity claim could be an artifact of the n=10 cutoff.","rationale":"The reader identifies the 2qp truncation and the 11-shell basis as the weakest premise; my read sharpens this to a specific, testable defect: the quantitative collectivity metric is defined relative to N2qp but never subjected to QRPA convergence tests. This is load-bearing because the abstract's central claim is exactly the 'moderate collectivity' statement, and the only quantitative support is N*, R, and delta E, all computed at n=10. The 2qp truncation (no phonon coupling) is a separate, well-known limitation, but the paper's conclusions are framed within QRPA; the basis-size issue is a correctable internal check that should precede any comparison with beyond-QRPA models. I do not see internal inconsistency: the QRPA implementation appears standard and the HFB+QRPA framework is fully consistent. However, calibrating n=10 solely on HFB energies is insufficient, because the central observable is a fragmentation statistic that depends explicitly on the number of 2qp configurations. The proposed n=12/14 rerun would settle whether the collectivity result is stable; if it is, the main claim withstands, and if not, the verdict should be conditional on reanalysis. Thus I agree with the reader's CONDITIONAL verdict and recommend no change to it.","tokens_in":15390,"tokens_out":6069,"duration_ms":66493,"concrete_test":"Repeat the HFB+QRPA calculation for representative isotopes 82Mo, 94Mo, and 98Mo with maximum oscillator shell n=12 and n=14 (same Gogny D1M, same convergence criteria and 8-16 MeV window). Recompute the low-energy B(E1) distribution, N* for skin-oscillation and major PDR states, R=N*_P/N*_G, and delta E/<E2qp>. If R changes by more than roughly 10%, or if the relative ordering/classification of skin-oscillation vs major-peak states changes, the moderate-collectivity conclusion is not robust to basis truncation. A secondary check: compare the radial transition densities for the major PDR peaks at n=10 and n=12 to ensure the surface/IS-IV pattern used to classify states is stable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central load-bearing premise is not the 2qp truncation per se, but the basis-size dependence of the metric used to justify the headline. Section II.D selects n=10 by HFB energy convergence (Table I, Fig. 1), where the n=10 HFB energy differs from n=16 by less than 0.2%. However, the collectivity conclusion in Section III.D rests on the fragmentation number N* of Eq. (13), which counts 2qp configurations with weight above 1/N2qp. At n=10, N2qp=12152 for K^pi=0^-; increasing the basis to n=12 or n=14 would enlarge N2qp, lower the threshold 1/N2qp, and generally change N* for skin-oscillation states, major PDR peaks, and the GDR reference. The ratio R=N*_P/N*_G and the energy shifts delta E (Eq. 14), which are the paper's quantitative evidence for 'moderate collectivity ... substantial configuration mixing, but limited coherence,' are therefore not shown to be converged. No QRPA-level convergence test for B(E1), transition densities, N*, R, or delta E is reported. The QRPA response, especially the high-lying GDR used as reference, can be more sensitive to the basis than the HFB ground-state energy. If R or the relative energy shifts move significantly with basis size, the headline distinction between skin-oscillation states and GDR states would be an artifact of the n=10 cutoff rather than a physical finding.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents fully consistent HFB+QRPA calculations with the Gogny D1M interaction for the electric and isoscalar dipole responses of the spherical even-even molybdenum isotopes 82Mo to 98Mo. It reports a low-energy dipole enhancement that correlates with the development of neutron or proton skins, and decomposes the enhancement into skin-oscillation states and a major PDR peak on the basis of radial transition densities. The collectivity of these states is quantified through a fragmentation number N* (Eq. 13), a relative fragmentation ratio R = N*_P/N*_G, and energy shifts delta E (Eq. 14), leading to the central claim that skin-oscillation states show moderate collectivity with substantial configuration mixing but limited coherence, whereas GDR states show strong coherence and large energy shifts.","tokens_in":15626,"tokens_out":4947,"duration_ms":61978,"significance":"If the central claim is robust, the paper provides a useful systematic survey of dipole excitations across a long isotopic chain, including proton-rich and neutron-rich sides, using a forward calculation in which the Gogny D1M interaction is taken from earlier work and not adjusted to the Mo dipole response. The fully consistent HFB+QRPA setup, the inclusion of all 2qp configurations without an energy cutoff, the explicit formulas for transition densities, and the study of HFB basis convergence are strengths. The decomposition into skin-oscillation and major-PDR states, together with quantitative collectivity indicators, is of interest for ongoing debates about the nature of the PDR. However, the central collectivity conclusion rests on quantities whose basis-size convergence is not demonstrated; this is the main gap that needs to be addressed.","major_comments":[{"comment":"The basis-size convergence of the collectivity metrics is not established. The convergence study in Section II.D (Table I, Fig. 1) is performed only at the HFB level, showing that the HFB energy for n=10 differs from n=16 by less than 0.2%. However, Eq. (13) defines N* by counting 2qp configurations with a_w_2qp above the threshold 1/N2qp, and N2qp depends explicitly on the number of major oscillator shells (N2qp = 12152 for K^pi = 0^- at n=10). Increasing the basis to n=12 or n=14 enlarges N2qp, lowers the threshold, and can change N* for skin-oscillation states, major PDR peaks, and the GDR reference in different ways. Since the ratio R = N*_P/N*_G and the energy shifts delta E are the quantitative basis for the claim of moderate collectivity in Section III.D, the absence of QRPA-level convergence tests for B(E1), transition densities, N*, R, and delta E leaves open the possibility that the central distinction is an artifact of the n=10 cutoff. I request such tests for representative isotopes (at least 82Mo and 94Mo) or a demonstration that R and delta E are stable under basis enlargement.","section":"Section II.D, Eq. (13), and Section III.D"},{"comment":"The identification of 'skin oscillation states' is not defined by a quantitative criterion. The paper selects representative states by visual inspection of their transition densities, but Figs. 4(c), 6, and 7 aggregate over all skin-oscillation states. It is not stated how many states are classified in this way, what threshold for surface neutron or proton dominance is used, or whether the classification is stable under small changes in the energy window (8-16 MeV) or in the smoothing procedure. Since the decomposition of the low-energy enhancement into skin-oscillation and major-PDR contributions underlies Fig. 4(c) and the subsequent collectivity analysis, a precise, reproducible selection criterion is needed, together with a check of robustness against reasonable variations of that criterion.","section":"Section III.C and Fig. 5"},{"comment":"The paper explicitly restricts the excitation space to 2qp configurations ('Excitations beyond 2qp configurations are not considered'), and the conclusion about moderate collectivity is therefore a statement about the QRPA 2qp space. This is acknowledged in the methods, but the abstract and conclusion present the finding without this qualifier. Given the literature cited in the introduction showing that phonon coupling can change the fragmentation and collectivity of PDR states, I recommend that the conclusion be framed as a QRPA-level result and that the possible impact of beyond-QRPA correlations be stated explicitly as a limitation.","section":"Section II.A and Section IV"}],"minor_comments":[{"comment":"The energy shift is defined as delta E_w = E_w - sum_mu E_mu |a_mu|, but the weights |a_mu| are not normalized to unity because Eq. (12) normalizes the signed sum of a_mu to 1. Using |a_mu| / sum_mu |a_mu| in Eq. (14) would make delta E independent of the arbitrary overall normalization of the amplitudes.","section":"Eq. (14) and surrounding text"},{"comment":"The choice of 16 MeV as the upper boundary of the low-energy dipole region is motivated by the EWSR fractions, but the paper does not discuss how the main conclusions would change if a different cutoff, for example one tied to the experimental PDR region around 6-10 MeV in 94Mo, were adopted. A short robustness statement would be helpful.","section":"Section III.B"},{"comment":"There are several language issues, for example 'This results correlates' in Section III.D and 'the feature discussed here indicate' in Section IV. The text would benefit from a careful proofread.","section":"General"},{"comment":"The caption mentions 'gray lines' for the discrete spectra, but the gray lines are not clearly visible in the figure. Please check the visibility and clarify the description of the discrete spectra.","section":"Fig. 3"},{"comment":"The paper does not report quantitative comparisons with experimental data for the dipole strength in Mo isotopes, despite citing experimental work on 94Mo and other isotopes. A comparison of the calculated B(E1) distributions or gamma-strength functions with at least one measured case would strengthen the credibility of the predictions.","section":"Section III.D"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern about the basis-size dependence of N* is real and load-bearing. The paper is honest about the 2qp truncation and does not fit parameters to the dipole response, which is good. The main issue is that the central collectivity claim is based on a metric that is explicitly basis dependent, and the manuscript only checks basis convergence of HFB ground-state energies. Adding QRPA-level convergence tests for the collectivity indicators should resolve the concern. The classification criterion for skin-oscillation states also needs to be made quantitative. I believe the paper is suitable for the journal after a major revision addressing these points."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a clean, internally consistent QRPA calculation. The authors take Gogny D1M, build HFB ground states, run fully consistent QRPA, and analyze transition densities across a whole chain of spherical Mo isotopes from 82 to 98. That systematic span, including the proton-rich side where proton skins appear, is genuinely new and useful. They do not tune anything to the dipole response, which keeps the calculation honest. The EWSR fractions are given, the qualitative correlation between skin thickness and low-energy strength is clearly shown, and the transition-density decomposition into skin-oscillation states versus a major PDR peak is a reasonable way to organize the results. Credit is due for stating the 2qp limitation up front and for making the formulas for the collectivity metrics explicit.\n\nThe soft spots are real but not fatal. The most important one is precisely the stress-test concern: the fragmentation number N* in Eq. (13) depends explicitly on N2qp, which grows with the number of oscillator shells. The paper checks HFB energy convergence with n=10 versus n=16, but never checks whether the QRPA spectrum, B(E1) distributions, N*, the ratio R, or the energy shifts are converged with basis size. Since the headline distinction between \"moderate collectivity\" in skin states and \"strong coherence\" in the GDR rests on those metrics, an actual basis-size test is needed. The 2qp truncation is the other obvious caveat; the authors admit it, and the literature suggests phonon coupling can change fragmentation patterns. The 8–16 MeV window is somewhat arbitrary, though the EWSR argument gives it some support. I also note there is no quantitative comparison to experimental Mo data, which would help anchor the interpretation.\n\nNone of this sinks the paper. The central claim is plausible and conditionally stated. A serious referee should ask for QRPA-level convergence checks and, ideally, a comparison with available (gamma,xn) or NRF data. The paper is worth engaging with, especially for people working on PDR structure or needing transition densities for reaction calculations. I would send it to review, and I would cite it as a reference calculation for Mo isotopes, though I would not rely on the collectivity classification until the convergence question is settled.","headline":"A solid forward HFB+QRPA study of low-energy dipole states in Mo isotopes with a plausible skin-oscillation interpretation, but the collectivity metric's basis-size dependence is not checked, so the central claim stays conditional.","tokens_in":16268,"tokens_out":1725,"would_cite":true,"duration_ms":21645,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["21.60.Jz"],"model":"deepseek-v4-flash","headline":"This paper argues that the low-energy dipole enhancement in molybdenum isotopes is not a single collective pygmy resonance but a mix of two distinct types of states: moderately collective skin oscillations and one more GDR-like peak.","keywords":["pygmy dipole resonance","quasiparticle random phase approximation","Hartree-Fock-Bogoliubov","Gogny D1M","neutron skin","proton skin","transition densities","nuclear collectivity"],"falsifier":"Measure the E1 strength function in 82–98Mo with resolution sufficient to resolve individual low-energy states, for example through high-resolution (γ,γ′) or (p,p′) experiments, and compare the number and B(E1) values of transitions in the 8–16 MeV window: if the observed fragmentation is far denser than the few dominant QRPA states, or if the energy centroids do not follow the predicted skin-correlation trend, the two-quasiparticle-only picture is falsified.","tokens_in":15100,"feed_emoji":"⚛️","tokens_out":6145,"duration_ms":70224,"temperature":0.7,"pith_summary":"This paper argues that the long-debated pygmy dipole resonance in molybdenum isotopes is not a single collective vibration. Using a fully consistent Hartree-Fock-Bogoliubov plus quasiparticle random phase approximation with the Gogny D1M interaction, the authors separate the low-energy electric dipole enhancement into two classes: skin oscillation states, in which the neutron or proton skin moves against an in-phase core, and a single major low-energy peak whose transition density resembles the giant dipole resonance. They conclude that the skin oscillation states are only moderately collective, with substantial configuration mixing but limited coherence, while the GDR states are strongly coherent. The work matters because it gives a microscopic picture of where the enhancement comes from and ties its size to whether a neutron or proton skin develops.","feed_headline":"Pygmy dipole in Mo is really two modes","feed_subtitle":"QRPA shows skin oscillations are only moderately collective; the main low-energy peak mimics the GDR.","key_machinery":"The load-bearing object is the QRPA excited-state wave function built on HFB ground states in an 11-oscillator-shell basis, with the same Gogny D1M force used in both steps. Two diagnostics carry the argument: the radial neutron and proton transition densities $\\delta\\rho_n(r)$ and $\\delta\\rho_p(r)$, which classify each state as in-phase/isoscalar, surface-dominated skin oscillation, or out-of-phase/isovector GDR-like; and the pair of coherence measures consisting of the fragmentation number $N^* = \\sum_{2qp}\\Theta(a^{w}_{2qp}-1/N_{2qp})$ and the relative energy shift $\\delta E/\\langle E_{2qp}\\rangle$, which together separate having many configurations from having coherent motion.","core_discovery":"In the spherical even-even molybdenum isotopes 82Mo to 98Mo, the enhancement of E1 strength near the neutron separation energy is correlated with the development of a neutron skin (92–98Mo) or a proton skin (82–88Mo). Radial transition densities show that inside the nucleus proton and neutron densities oscillate in phase, while beyond the surface one nucleon species dominates, matching the skin type. The states carrying this surface oscillation have fragmentation numbers comparable to the GDR, but small relative energy shifts and suppressed B(E1) values because their isoscalar character causes partial cancellations. The single strongest low-energy peak, by contrast, shows a GDR-like out-of-phase pattern. The central claim is that the low-energy dipole response in Mo is a mixture of moderately collective skin oscillations and one more GDR-like peak, not a single coherent pygmy mode.","pith_inferences":["If this decomposition survives comparison with data, it suggests that energy-density-functional estimates of the symmetry energy from low-energy E1 strength should focus on the skin-oscillation component rather than the total strength, since the major peak tracks a different mode.","The same transition-density classification could be applied to deformed nuclei by following the K-quantum number; shape deformation is expected to mix the isoscalar and isovector patterns and may split the single major peak into a fragmented multiplet, a testable prediction.","Because transition densities feed reaction calculations, the densities produced here could be used to predict (p,p′) or (α,α′) cross sections for 82–98Mo, providing an independent test of the skin-oscillation versus GDR-like assignment that does not rely on B(E1) alone."],"forward_implications":["In 92–98Mo the skin-oscillation contribution to the 8–16 MeV EWSR grows as the neutron skin thickens; in 84–90Mo, where a proton skin develops, it is nearly absent, and the enhancement is dominated by the single major PDR peak.","Skin oscillation states will not behave like textbook collective modes in reactions: their isoscalar-dominated transition densities suppress B(E1), so experimental signatures should appear mainly in isoscalar probes such as (α,α′γ) rather than in photoabsorption alone.","The average energy of the low-energy states lies below the neutron separation energy for proton-skin isotopes and above it for neutron-skin isotopes, so the location of the enhancement relative to threshold is skin-driven rather than universal.","If the major PDR peak is removed from the E1 response, the remaining low-energy enhancement almost disappears for 84–90Mo, meaning that for those isotopes the PDR label applies to a single state rather than to a resonance-like accumulation.","The GDR-like major peak and the skin oscillation states respond differently to the dipole operator, so total B(E1) alone cannot be used as a clean measure of pygmy collectivity."],"supporting_citations":[{"why":"Supplies the HFB+QRPA framework and the fragmentation-number diagnostics used to quantify collectivity of low-energy dipole states.","marker":"[16]"},{"why":"Defines the Gogny effective interaction on which both the HFB and QRPA calculations are built.","marker":"[36]"},{"why":"Provides the D1M parametrization chosen for its improved mass predictions and performance in neutron-rich nuclei.","marker":"[37]"},{"why":"Gives the experimental 94Mo gamma-ray strength function that motivates comparison with QRPA E1 results.","marker":"[33]"},{"why":"Presents quasiparticle-phonon model results that challenge a collective interpretation of the PDR, the debate this paper engages.","marker":"[17]"},{"why":"Tracks the PDR peak energy relative to the neutron separation energy across isotope chains, a trend the Mo results extend.","marker":"[19]"},{"why":"Shows that destructive interference among configurations can keep PDR transition strengths below strongly collective GDR values, a mechanism central to the coherence analysis.","marker":"[14]"}],"fun_headline_variants":["Mo dipole: skin mode plus GDR-like peak","QRPA splits Mo's low-energy dipole into two","Molybdenum's dipole response is two-faced","Not one pygmy mode: Mo shows two dipole types","Skin oscillations and a collective peak in Mo"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The results depend on the assumption that the restriction to two-quasiparticle excitations, with no coupling to more complicated configurations, is enough to capture the fragmentation and coherence of these low-energy dipole states.","fun_headline_variants_meta":{"raw":{"variants":["Mo dipole: skin mode plus GDR-like peak","QRPA splits Mo's low-energy dipole into two","Molybdenum's dipole response is two-faced","Not one pygmy mode: Mo shows two dipole types","Skin oscillations and a collective peak in Mo"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000306,"raw_usage":{"total_tokens":1782,"prompt_tokens":1001,"completion_tokens":781,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":617,"completion_tokens_details":{"reasoning_tokens":707}},"tokens_in":617,"tokens_out":781,"duration_ms":8331,"temperature":1.0,"reasoning_tokens":707,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T10:53:31.624768+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the E1 strength function in 82–98Mo with resolution sufficient to resolve individual low-energy states, for example through high-resolution (γ,γ′) or (p,p′) experiments, and compare the number and B(E1) values of transitions in the 8–16 MeV window: if the observed fragmentation is far denser than the few dominant QRPA states, or if the energy centroids do not follow the predicted skin-correlation trend, the two-quasiparticle-only picture is falsified.","supporting_citations":[{"cited_title":"Martini, S","cited_arxiv_id":null,"evidence_quote":"Supplies the HFB+QRPA framework and the fragmentation-number diagnostics used to quantify collectivity of low-energy dipole states."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the experimental 94Mo gamma-ray strength function that motivates comparison with QRPA E1 results."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Tracks the PDR peak energy relative to the neutron separation energy across isotope chains, a trend the Mo results extend."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows that destructive interference among configurations can keep PDR transition strengths below strongly collective GDR values, a mechanism central to the coherence analysis."}],"review_version":1}