Electric dipole excitations near the neutron separation energies in $^{96}$Mo

Aaina Thapa; Emanuel Chimanski; Eun Jin In; Jutta Escher; Sophie P\'eru; Walid Younes

REVIEW 4 major objections 5 minor 13 references

Electric dipole excitations near the neutron separation energies in $^{96}$Mo

T0 review · 4 major / 5 minor · reviewed 2026-08-08 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read 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.

desk verdict 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. read the letter →

arxiv 2502.07113 v1 pith:FS6SKLOJ submitted 2025-02-10 nucl-th

Eun Jin In , Emanuel Chimanski , Jutta Escher , Sophie Péru , Aaina Thapa , Walid Younes This is my paper · ORCID

classification nucl-th MSC 81V35 PACS 21.60.Jz21.10.-k23.20.-g24.30.Cz

keywords pygmydipoleresonanceelectricstrengthtransitiondensitiesHartree-Fock-BogoliubovQuasiparticleRandomPhaseApproximationGognyD1Misovectorisoscalar

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

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.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the 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.

The 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.

This 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.

Who 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.

Referee Report

4 major / 5 minor

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.

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 (4)

[Sec. 3.1, Fig. 2] 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.
[Sec. 3.2, Fig. 3(b)] 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.
[Sec. 3.2] 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.
[Sec. 4, conclusion] 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.

minor comments (5)

[Fig. 2] 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.
[Fig. 3] 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.
[Sec. 1] 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.
[Sec. 2] 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.
[Sec. 3.1] The sentence 'Figures 3 displays...' contains a subject-verb agreement error; it should be 'Figure 3 displays...'.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the QRPA E1 response and transition densities are genuine predictions of the fixed Gogny D1M model, with no fitted inputs from 96Mo data.

full rationale

The derivation chain is self-contained. The paper computes the 96Mo dipole response by solving HFB+QRPA with the Gogny D1M finite-range interaction (Sec. 2), an interaction fixed in prior work and not adjusted to 96Mo observables. No experimental dipole strength, separation-energy, or transition-density data for 96Mo are used to set or tune any parameter; hence the near-separation-energy enhancement and the transition-density isoscalar/isovector classification are model predictions rather than fits. The Lorentzian folding width Gamma=1 MeV (Sec. 3.1) is a plotting and smoothing choice that does not enter the discrete QRPA spectrum; its lack of sensitivity analysis is a robustness limitation, not a circular step, because the continuous curves are not used to constrain any model parameter. The isospin characterization follows from the computed radial transition densities via the operators in Eqs. (1)-(2) and the standard PDR signature cited from independent references [13-15], not from assuming the conclusion. Reference [9] is authored in part by the present group but is cited only as background for the JLM transition-density application framework and does not carry the QRPA prediction. Therefore no load-bearing step reduces to its own input.

Assumptions & free parameters 1 free parameters · 4 assumptions · 0 invented entities

The central prediction rests on the standard nuclear-structure framework (HFB+QRPA) with an external effective interaction (Gogny D1M), plus one hand-chosen smoothing parameter. No new particles, forces, or conserved quantities are introduced. The main unpaid inputs are the reliability of the interaction and the QRPA approximation for the low-energy E1 response, and the representativeness of the states selected for the isospin classification.

free parameters (1)

Lorentzian folding width Gamma = 1 MeV
Chosen by hand in Sec. 3.1 to fold the discrete QRPA spectrum into continuous response functions in Fig. 2. No physical justification or sensitivity analysis is given, and the visibility of the claimed PDR enhancement depends on this value.

assumptions (4)

domain assumption The Gogny D1M effective interaction, whose parameters were fitted to nuclear ground-state observables in prior works, provides a reliable description of the HFB ground state and QRPA excitations of 96Mo.
Sec. 2 states the calculations use the same Gogny D1M finite-range effective interaction in a fully consistent manner; no re-fitting or validation against 96Mo dipole data is presented.
domain assumption The QRPA, built from two-quasiparticle configurations on the HFB ground state, captures the low-energy E1 response and its isospin character near the neutron separation energy.
Invoked in Sec. 2 to construct the excited spectrum; the paper does not assess beyond-QRPA correlations (phonon coupling, continuum) known to affect PDR fragmentation.
domain assumption The center-of-mass corrections in the E1 and isoscalar dipole operators (Eqs. 1 and 2) remove the spurious translational mode.
Sec. 2 states the corrections restore translational invariance; the paper provides no check that the spurious state is fully removed in the computed spectrum.
ad hoc to paper The states selected in Fig. 3 are representative of the PDR region and of the dominant component of the enhanced low-energy region.
Sec. 3.2 labels the displayed states as representative and generalizes from them to the isospin character of the whole enhancement region; the selection criteria are not specified.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Electric dipole excitations near the neutron separation energies in $^{96}$Mo." pith.science (2026). https://pith.science/paper/FS6SKLOJ

@misc{pith2026250207113,
  author       = {Pith},
  title        = {Pith review of: Electric dipole excitations near the neutron separation energies in $^96$Mo},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/FS6SKLOJ}},
  note         = {Machine review of arXiv:2502.07113}
}

abstract

Electric dipole strength near the neutron separation energy significantly impacts nuclear structure properties and astrophysical scenarios. These excitations are complex in nature and may involve the so-called pygmy dipole resonance (PDR). Transition densities play a crucial role in understanding the nature of nuclear excited states, including collective excitations, as well as in constructing transition potentials in DWBA or coupled-channels equations. In this work, we focus on electric dipole excitations in spherical molybdenum isotopes, particularly $^{96}$Mo, employing fully consistent Hartree-Fock-Bogoliubov (HFB) and Quasiparticle Random Phase Approximation (QRPA) methods. We analyze the dipole strength near the neutron separation energy, which represents the threshold for neutron capture processes, and examine the isospin characteristics of PDR states through transition density calculations. Examination of proton and neutron transition densities reveals distinctive features of each dipole state, indicating their isoscalar and isovector nature. We observe that the primary component in the enhanced low-energy region exhibits isovector character. The PDR displays a mixture of isoscalar and isovector nature, distinguishing it from the isovector giant dipole resonance (IVGDR). These findings lay the groundwork for future investigations into the role of transition densities in reaction models and for their application to inelastic scattering calculations.

Figures

Figures reproduced from arXiv: 2502.07113 by the authors.

**Figure 2.** The electric dipole (dB(E1)/dE) and isoscalar (dBIS (E1)/dE) response functions for 96Mo as a function of the QRPA excitation energy. The dashed vertical line represent the neutron separation energy for 96Mo. The PDR region is highlighted by a vertical band. observe significant excess in the neutron density in this space range, indicating the formation of a neutron skin. 84Mo, on the other hand, does not exhibit suc… view at source ↗

**Figure 3.** Radial transition densities for 96Mo isotope are shown for three cases: (a) the neutron PDR candidate, (b) the major peak of the enhancement region, and (c) the major GDR peak with the largest B(E1). the PDR [13–15]. On the other hand, the largest peak in the potential PDR region, at 13.5 MeV in 96Mo, exhibits an isovector nature similar to that of the IVGDR. Further investigation is needed to better understand thi… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

13 extracted references · 9 canonical work pages

[1]

N. Paar, D. Vretenar, E. Khan, G. Colò, Rep. Prog. Phys 70, R02 (2007). 10.1088/0034-4885/70/5/R02

work page doi:10.1088/0034-4885/70/5/r02 2007
[2]

Lanza, L

E. Lanza, L. Pellegri, A. Vitturi, M. Andrés, Prog. Part. Nucl. Phys. 129, 104006 (2023). https://doi.org/10.1016/j.ppnp.2022.104006

arXiv 2023
[3]

Tonchev, N

A. Tonchev, N. Tsoneva, C. Bhatia, C. Arnold, S. Goriely, S. Hammond, J. Kel- ley, E. Kwan, H. Lenske, J. Piekarewicz et al., Phys. Lett. B 773, 20 (2017). https://doi.org/10.1016/j.physletb.2017.07.062

work page doi:10.1016/j.physletb.2017.07.062 2017
[4]

Wieser, J.R.D

M.E. Wieser, J.R.D. Laeter, Phys. Rev. C 75, 055802 (2007). 10.1103/PhysRevC.75.055802

work page doi:10.1103/physrevc.75.055802 2007
[5]

Stephan, R

T. Stephan, R. Trappitsch, P . Hoppe, A.M. Davis, M.J. Pellin, O.S. Pardo, ApJ 877, 101 (2019). 10.3847/1538-4357/ab1c60

work page doi:10.3847/1538-4357/ab1c60 2019
[6]

Bauge, J.P

E. Bauge, J.P . Delaroche, M. Girod, Phys. Rev. C 58, 1118 (1998). 10.1103/PhysRevC.58.1118

work page doi:10.1103/physrevc.58.1118 1998
[7]

Bauge, J.P

E. Bauge, J.P . Delaroche, M. Girod, Phys. Rev. C 63, 024607 (2001). 10.1103/PhysRevC.63.024607

work page doi:10.1103/physrevc.63.024607 2001
[8]

Dupuis, G

M. Dupuis, G. Haouat, J.P . Delaroche, E. Bauge, J. Lachkar, Phys. Rev. C 100, 044607 (2019). 10.1103/PhysRevC.100.044607

work page doi:10.1103/physrevc.100.044607 2019

Show all 13 references

[9]

Thapa, J

A. Thapa, J. Escher, E. Chimanski, M. Dupuis, S. Péru, W . Y ounes, EPJ Web Conf. 292, 06003 (2024). 10.1051/epjconf/202429206003

2024
[10]

Escher, J.T

J.E. Escher, J.T. Harke, F.S. Dietrich, N.D. Scielzo, I.J. Thompson, W . Y ounes, Rev. Mod. Phys.84, 353 (2012). 10.1103/RevModPhys.84.353

2012 doi
[13]

Vretenar, N

D. Vretenar, N. Paar, P . Ring, G.A. Lalazissis, Phys. Rev. C 63, 047301 (2001). 10.1103/Phys- RevC.63.047301

2001 doi
[14]

Tsoneva, H

N. Tsoneva, H. Lenske, Phys. Rev. C 77, 024321 (2008). 10.1103/PhysRevC.77.024321

2008 doi
[15]

Co’, V .D

G. Co’, V .D. Donno, C. Maieron, M. Anguiano, A.M. Lallena, Phys. Rev. C 80, 014308 (2009). 10.1103/PhysRevC.80.014308

2009 doi

Pith tools

Reviewed August 8, 2026 · model on record in the stance chip above.

[1] [1]

N. Paar, D. Vretenar, E. Khan, G. Colò, Rep. Prog. Phys 70, R02 (2007). 10.1088/0034-4885/70/5/R02

work page doi:10.1088/0034-4885/70/5/r02 2007

[2] [2]

Lanza, L

E. Lanza, L. Pellegri, A. Vitturi, M. Andrés, Prog. Part. Nucl. Phys. 129, 104006 (2023). https://doi.org/10.1016/j.ppnp.2022.104006

arXiv 2023

[3] [3]

Tonchev, N

A. Tonchev, N. Tsoneva, C. Bhatia, C. Arnold, S. Goriely, S. Hammond, J. Kel- ley, E. Kwan, H. Lenske, J. Piekarewicz et al., Phys. Lett. B 773, 20 (2017). https://doi.org/10.1016/j.physletb.2017.07.062

work page doi:10.1016/j.physletb.2017.07.062 2017

[4] [4]

Wieser, J.R.D

M.E. Wieser, J.R.D. Laeter, Phys. Rev. C 75, 055802 (2007). 10.1103/PhysRevC.75.055802

work page doi:10.1103/physrevc.75.055802 2007

[5] [5]

Stephan, R

T. Stephan, R. Trappitsch, P . Hoppe, A.M. Davis, M.J. Pellin, O.S. Pardo, ApJ 877, 101 (2019). 10.3847/1538-4357/ab1c60

work page doi:10.3847/1538-4357/ab1c60 2019

[6] [6]

Bauge, J.P

E. Bauge, J.P . Delaroche, M. Girod, Phys. Rev. C 58, 1118 (1998). 10.1103/PhysRevC.58.1118

work page doi:10.1103/physrevc.58.1118 1998

[7] [7]

Bauge, J.P

E. Bauge, J.P . Delaroche, M. Girod, Phys. Rev. C 63, 024607 (2001). 10.1103/PhysRevC.63.024607

work page doi:10.1103/physrevc.63.024607 2001

[8] [8]

Dupuis, G

M. Dupuis, G. Haouat, J.P . Delaroche, E. Bauge, J. Lachkar, Phys. Rev. C 100, 044607 (2019). 10.1103/PhysRevC.100.044607

work page doi:10.1103/physrevc.100.044607 2019