{"id":"31d04d62-92db-458f-8efe-ecab8678797f","arxiv_id":"2502.07113","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A fully consistent HFB+QRPA calculation with the Gogny D1M interaction predicts a pygmy dipole resonance in 96Mo near the neutron separation energy, with the dominant low-energy component being isovector while a candidate PDR state has mixed isoscalar/isovector character.","lead":"This paper uses a standard microscopic nuclear model to calculate how the molybdenum-96 nucleus responds to electric dipole radiation, finding a small enhancement of dipole strength near the neutron separation energy. The result matters because this low-energy dipole strength, called the pygmy dipole resonance, affects neutron capture rates important for making heavy elements in stars and for nuclear data applications.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed near-S_n enhancement and its isovector/isoscalar classification rest on an untested 1 MeV folding width and on selected states, one at 13.5 MeV which appears inconsistent with the empirical S_n of 96Mo.","rationale":"The paper presents a legitimate fully consistent HFB+QRPA calculation with the Gogny D1M interaction, and the transition-density analysis is a useful diagnostic. However, the specific claims about a pygmy enhancement near S_n and its isovector/isoscalar nature depend on two unverified choices: the Lorentzian folding width and the selection of representative states. The apparent placement of the largest 'PDR region' peak at 13.5 MeV, well above the empirical S_n of 96Mo, strengthens the concern that the analysis may not actually pertain to the near-threshold region. These are not internal contradictions in the QRPA machinery, but they are load-bearing for the interpretation. The reader's weakest-assumption assessment correctly identified the smoothing and state-selection dependence, and the recommended CONDITIONAL verdict remains appropriate. A revision that shows the discrete spectrum, a Gamma-sensitivity study, the numerical S_n used, and integrated strengths in a well-defined S_n window would materially strengthen the claims.","tokens_in":4528,"tokens_out":10462,"duration_ms":92994,"concrete_test":"Report the numerical S_n used for Fig. 2 and overlay the discrete QRPA B(E1) spectrum with S_n marked; repeat the Lorentzian folding with Gamma=0.1, 0.5, 1, and 2 MeV. Then integrate B(E1) and isoscalar strength over an S_n +/- 2 MeV window and compare with the single 13.5 MeV state. If the 13.5 MeV state lies outside that window, or if the window-integrated isovector fraction differs from the single-state classification, the near-S_n enhancement claim is not supported by the presented analysis.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that the small bump near the neutron separation energy in Fig. 2 is a robust concentration of E1 strength and that the states analyzed in Fig. 3 are representative of that bump. Two features undermine this. First, the continuous dB(E1)/dE and dB_IS(E1)/dE curves are obtained by folding the discrete QRPA spectrum with Lorentzians of width Gamma=1 MeV (Sec. 3.1), but the zero-width spectrum is not shown and no sensitivity to Gamma is reported. If the QRPA eigenvalues are sparse near S_n, a 1 MeV width can merge isolated transitions into an apparent bump whose height and centroid are smoothing artifacts. Second, Sec. 3.2 calls both panels (a) and (b) 'states in the potential PDR region' and identifies panel (b) as the largest peak of that region at 13.5 MeV, yet the empirical neutron separation energy of 96Mo is about 9.15 MeV. A 13.5 MeV state lies more than 4 MeV above the threshold and cannot support a claim about strength near S_n unless the dashed vertical line in Fig. 2 uses a different definition, which is not stated. Thus the isovector classification of the 13.5 MeV peak does not by itself establish the nature of the near-threshold enhancement, and the conclusion that the dominant low-energy component is isovector depends on an unverified representativeness assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports fully consistent Hartree-Fock-Bogoliubov (HFB) and Quasiparticle Random Phase Approximation (QRPA) calculations with the Gogny D1M interaction for 96Mo, focusing on electric dipole (E1) and isoscalar dipole excitations near the neutron separation energy. The authors compute ground-state densities, fold the discrete QRPA dipole spectrum with Lorentzian functions of width 1 MeV, and identify a small enhancement near the separation energy, which they associate with the pygmy dipole resonance (PDR). They then analyze the radial transition densities of two selected states in the low-energy enhancement region and one giant dipole resonance state, concluding that the dominant low-energy components are isovector while a PDR candidate exhibits mixed isoscalar-isovector character. The introduction and conclusions also claim a correlation between the enhancement and neutron excess, based on the comparison of ground-state densities of 84Mo and 96Mo.","tokens_in":4746,"tokens_out":5403,"duration_ms":45437,"significance":"If substantiated, the calculation provides a parameter-free prediction (no parameters fitted to 96Mo data) of the dipole response and transition densities for a spherical nucleus of astrophysical interest. The transition densities are potentially valuable as inputs for inelastic-scattering and neutron-capture reaction models, which is a stated goal of the authors. However, the central claims about the PDR enhancement and its isospin classification currently rest on an arbitrary smoothing width and on visual inspection of two states, which limits the paper's scientific impact until these points are strengthened.","major_comments":[{"comment":"The continuous dipole strength curves are obtained by folding the discrete QRPA spectrum with Lorentzian functions of width Γ = 1 MeV, but the discrete spectrum is not shown and no sensitivity study of Γ is reported. If the QRPA eigenvalues are sparse in the region near the separation energy, a 1 MeV width can merge isolated transitions into an apparent enhancement whose height and centroid are smoothing artifacts. To support the central claim of a small enhancement near S_n, the authors should show the zero-width spectrum and repeat the folding with at least one additional width (e.g., Γ = 0.5 and 2 MeV) to demonstrate that the enhancement persists.","section":"Sec. 3.1, Fig. 2"},{"comment":"The text identifies the state at 13.5 MeV as the largest peak in the potential PDR region and uses its transition density to characterize the enhancement near the neutron separation energy. However, the empirical neutron separation energy of 96Mo is about 9.15 MeV, so a 13.5 MeV state lies more than 4 MeV above S_n. The manuscript does not state the numerical value of the dashed vertical line in Fig. 2 or the calculated S_n from the HFB model. If the enhancement is meant to be concentrated near S_n, the analysis must either focus on states within a narrow window around S_n or explain why a state at 13.5 MeV is relevant to the near-threshold region.","section":"Sec. 3.2, Fig. 3(b)"},{"comment":"The conclusion that the dominant low-energy component is isovector while the PDR state is mixed isospin rests on visual inspection of the transition densities of two selected states (panels (a) and (b) of Fig. 3). No quantitative criterion is given, such as the ratio of integrated neutron to proton transition density in the surface region or the overlap of each QRPA state with the isovector and isoscalar operators of Eqs. (1) and (2). Because the manuscript asserts a general property of the enhancement region, a quantitative measure should be applied to all QRPA states in that energy interval, not just to the two displayed states.","section":"Sec. 3.2"},{"comment":"The conclusion states that the observed enhancement 'correlates with neutron excess,' but no dipole response for 84Mo is shown; only ground-state densities are presented in Fig. 1. The correlation with neutron excess is therefore not demonstrated by the displayed results. The paper should either show the E1 response for 84Mo (or another neutron-deficient isotope) or temper the claim to apply only to 96Mo.","section":"Sec. 4, conclusion"}],"minor_comments":[{"comment":"The caption does not define the dashed vertical line or the shaded band; please state the numerical value of S_n used and the energy range of the PDR region.","section":"Fig. 2"},{"comment":"The caption does not identify the energies or B(E1) values of the states shown in panels (a), (b), and (c); please add this information so that the reader can connect the panels to the text.","section":"Fig. 3"},{"comment":"The phrase 'serve as a potential singular piece of experimental evidence' is awkward; consider rewording to 'serve as a potential piece of experimental evidence' or similar.","section":"Sec. 1"},{"comment":"A brief description of the QRPA configuration space (two-quasiparticle basis, energy cutoffs, and the treatment of spurious modes beyond the center-of-mass corrections in Eqs. (1)-(2)) would improve reproducibility.","section":"Sec. 2"},{"comment":"The sentence 'Figures 3 displays...' contains a subject-verb agreement error; it should be 'Figure 3 displays...'.","section":"Sec. 3.1"}],"recommendation":"major_revision","confidential_remarks":"This manuscript is a short proceedings-style contribution. The HFB+QRPA calculation with the Gogny D1M interaction is a standard and credible approach, and the transition densities are potentially useful for reaction modeling. The main weaknesses are the lack of sensitivity analysis for the Lorentzian folding width, the apparent inconsistency between the 13.5 MeV state and the empirical S_n of 96Mo, and the purely qualitative isospin classification. These points are load-bearing for the central claims, but they are addressable within the scope of a revised manuscript. If the authors provide the requested quantitative analyses, the paper could be suitable for publication in a proceedings volume or a theory journal focused on nuclear structure and reactions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely useful part here is the fully consistent HFB+QRPA calculation with Gogny D1M for 96Mo, a nucleus that seems not to have been treated this way before. The paper computes both electromagnetic and isoscalar dipole responses and gives radial transition densities for representative states, with no parameters fitted to 96Mo data. That makes the results a real prediction rather than a fit, and the transition densities could be directly useful for future DWBA or coupled-channels work. Credit is due for that.\n\nThe soft spots are real but not fatal. The \"small enhancement\" called a PDR is obtained by folding the discrete QRPA spectrum with a Lorentzian width of 1 MeV, and no zero-width spectrum or sensitivity test is shown. If the QRPA levels are sparse near threshold, the bump could be partly an artifact of the smoothing. The isospin classification is drawn from two states, and the paper calls both \"PDR region\" even though one of them, the largest peak, is at 13.5 MeV—more than 4 MeV above the empirical neutron separation energy of 96Mo (~9.15 MeV). That is a genuine inconsistency in how the enhancement region is defined and described. The authors themselves say further investigation is needed, which is honest, but the current text does not establish that the near-S_n strength is dominated by isovector states or that the 13.5 MeV peak is part of the same structure. The neutron-excess correlation is also based on just two isotopes (84Mo and 96Mo), so that claim is suggestive at best.\n\nThis is a proceedings paper, so the depth is limited, but the method is standard and the calculation is defensible as a first step. The main fixes are cheap: show the unfolded spectrum, test the folding width, and clarify what \"near the neutron separation energy\" and \"potential PDR region\" mean. A comparison with any available experimental dipole data for 96Mo would also strengthen the interpretation.\n\nWho is it for? Specialists in nuclear structure and reaction theory who want transition densities for 96Mo or who are mapping PDR systematics across isotopes. It is not a landmark, but it is a legitimate data point. I would send it to peer review with a request for revision rather than desk reject, because the underlying calculation is reproducible in principle and the transition densities could be re-used.","headline":"A solid but thin proceedings-style QRPA study of the E1 response in 96Mo; the PDR identification rests on an untested smoothing width and a confusingly labeled 13.5 MeV state, so the central claim needs sharper support.","tokens_in":5367,"tokens_out":1509,"would_cite":false,"duration_ms":15107,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81V35"],"pacs":["21.60.Jz","21.10.-k","23.20.-g","24.30.Cz"],"model":"deepseek-v4-flash","headline":"Fully consistent HFB+QRPA with Gogny D1M predicts a pygmy dipole enhancement near the neutron separation energy in 96Mo, with a dominant isovector low-energy component and a candidate state of mixed isoscalar-isovector character.","keywords":["pygmy dipole resonance","electric dipole strength","transition densities","Hartree-Fock-Bogoliubov","Quasiparticle Random Phase Approximation","Gogny D1M","isovector","isoscalar"],"falsifier":"Measure the electric dipole strength of 96Mo near the neutron separation energy with photon-scattering or inelastic-proton-scattering experiments and compare the energy distribution and angular distributions to the predicted isovector and isoscalar responses; if the observed enhancement is absent, or the angular distributions require a different isospin mixture than the one extracted from the transition densities, the claim is falsified.","tokens_in":4241,"feed_emoji":"⚛️","tokens_out":7776,"duration_ms":61769,"temperature":0.7,"pith_summary":"The paper predicts the electric dipole response of 96Mo near the neutron separation energy using fully consistent Hartree-Fock-Bogoliubov plus Quasiparticle Random Phase Approximation calculations with the Gogny D1M interaction. It finds a small enhancement in the dipole strength just below the neutron separation energy—a signature of the pygmy dipole resonance—and uses proton and neutron transition densities to characterize the excited states. The dominant component of the enhanced low-energy region is isovector, while a candidate PDR state displays mixed isoscalar and isovector character, distinguishing it from the isovector giant dipole resonance. These transition densities are the ingredients needed for folding-model predictions of inelastic scattering, which is why the isospin classification matters.","feed_headline":"96Mo shows predicted pygmy dipole strength near neutron separation","feed_subtitle":"Fully consistent HFB+QRPA calculations reveal an isovector-dominated enhancement and a mixed-mode PDR candidate.","key_machinery":"The central machinery is the fully consistent Hartree-Fock-Bogoliubov (HFB) plus Quasiparticle Random Phase Approximation (QRPA) framework, in which the same Gogny D1M finite-range effective interaction generates both the ground state and the correlated two-quasiparticle excited states. The electric dipole and isoscalar dipole operators (with center-of-mass corrections) define the response functions, and the radial transition densities extracted from the QRPA transition matrix elements are the tools that reveal the isoscalar versus isovector character of each state. The paper uses these transition densities to argue that the low-energy enhancement is mainly isovector while the candidate PDR state is mixed.","core_discovery":"Using the same Gogny D1M finite-range effective interaction for both the Hartree-Fock-Bogoliubov ground state and the Quasiparticle Random Phase Approximation excited states, the paper's calculations yield a small enhancement in the electric dipole strength of 96Mo near the neutron separation energy, a signature of a pygmy dipole resonance. The radial transition densities of the states in this region reveal that the dominant low-energy component is isovector, while a candidate PDR state—whose proton and neutron transition densities oscillate in phase in the interior and are neutron-dominated at the surface—has mixed isoscalar and isovector character. The paper presents this mixture as the distinguishing feature that sets the PDR apart from the isovector giant dipole resonance.","pith_inferences":["The isospin classification rests on only two representative states; a systematic decomposition of the full energy interval near the neutron separation energy would be needed to confirm that the mixed character is a property of the mode rather than of the selected states.","Because the discrete QRPA spectrum is folded with an arbitrary Lorentzian width of 1 MeV, checking the sensitivity of the enhancement to widths in the 0.5–2 MeV range would show whether the PDR signature is robust or a smoothing artifact.","The same transition densities could be used to generate isoscalar and isovector inelastic-scattering observables, which would allow a direct experimental discrimination of the predicted mixed character through angular-distribution comparisons.","If the mixed isoscalar-isovector character is confirmed, it would complicate the simple picture of the PDR as a pure neutron-skin oscillation and instead place 96Mo's pygmy mode in a transitional regime."],"forward_implications":["The calculated transition densities can be folded with a microscopic interaction to generate transition potentials for DWBA or coupled-channels calculations of inelastic scattering on 96Mo.","The prediction that the low-energy enhanced region is dominantly isovector means that electromagnetic probes and hadronic probes sensitive to isoscalar components should see different relative strengths across the PDR region.","The mixed isoscalar-isovector character of the candidate PDR state provides a specific signature that future (p,p') or (α,α') experiments can test.","The enhancement near the neutron separation energy implies an increase in neutron-capture cross sections relevant to the astrophysical s-process, as the PDR boosts the dipole strength available for capture.","The correlation of the neutron skin with the enhancement suggests that more neutron-rich molybdenum isotopes would exhibit an even more pronounced PDR."],"supporting_citations":[{"why":"Surveys microscopic models of the pygmy dipole resonance, establishing the context the paper's HFB+QRPA prediction extends.","marker":"[1]"},{"why":"Reviews transition-density analyses of the PDR and provides the isoscalar/isovector classification tools the paper applies.","marker":"[2]"},{"why":"Reports evidence that PDR strength enhances neutron-capture cross sections, motivating the astrophysical relevance of the predicted enhancement.","marker":"[3]"},{"why":"Provides the expected transition-density pattern of a PDR (in-phase oscillations in the interior, neutron-dominated surface) used to identify the candidate state.","marker":"[13]"},{"why":"Shows how transition densities encode isovector versus isoscalar character in PDR calculations, supporting the interpretation of the 96Mo states.","marker":"[14]"},{"why":"Analyzes transition densities of dipole states and serves as a reference for the radial behavior used to classify the 96Mo results.","marker":"[15]"}],"fun_headline_variants":["96Mo pygmy dipole: isovector boost, mixed nature","Near-threshold E1 in 96Mo: isovector dominant","Pygmy dipole in 96Mo shows mixed isoscalar-isovector","96Mo: isovector low-energy dipole, mixed PDR"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central claim assumes that the arbitrary 1 MeV Lorentzian folding width and the two selected states used for the transition-density analysis faithfully represent the physical pygmy dipole mode rather than an artifact of smoothing.","fun_headline_variants_meta":{"raw":{"variants":["96Mo pygmy dipole: isovector boost, mixed nature","Near-threshold E1 in 96Mo: isovector dominant","Pygmy dipole in 96Mo shows mixed isoscalar-isovector","96Mo: isovector low-energy dipole, mixed PDR"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00068,"raw_usage":{"total_tokens":3087,"prompt_tokens":942,"completion_tokens":2145,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":558,"completion_tokens_details":{"reasoning_tokens":2064}},"tokens_in":558,"tokens_out":2145,"duration_ms":14651,"temperature":1.0,"reasoning_tokens":2064,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T13:47:31.986045+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the electric dipole strength of 96Mo near the neutron separation energy with photon-scattering or inelastic-proton-scattering experiments and compare the energy distribution and angular distributions to the predicted isovector and isoscalar responses; if the observed enhancement is absent, or the angular distributions require a different isospin mixture than the one extracted from the transition densities, the claim is falsified.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Surveys microscopic models of the pygmy dipole resonance, establishing the context the paper's HFB+QRPA prediction extends."},{"cited_title":"Tonchev, N","cited_arxiv_id":null,"evidence_quote":"Reports evidence that PDR strength enhances neutron-capture cross sections, motivating the astrophysical relevance of the predicted enhancement."},{"cited_title":"Tsoneva, H","cited_arxiv_id":null,"evidence_quote":"Shows how transition densities encode isovector versus isoscalar character in PDR calculations, supporting the interpretation of the 96Mo states."},{"cited_title":"Co’, V .D","cited_arxiv_id":null,"evidence_quote":"Analyzes transition densities of dipole states and serves as a reference for the radial behavior used to classify the 96Mo results."}],"review_version":1}