pith. sign in

arxiv: 2604.11265 · v1 · submitted 2026-04-13 · ⚛️ nucl-th

Delta l =1 coupling of single-particle orbitals in octupole deformed nuclei

Xudong Wang , Bin Qi , Shouyu Wang , Chen Liu This is my paper

Pith reviewed 2026-05-10 15:42 UTC · model grok-4.3

classification ⚛️ nucl-th
keywords octupole deformationsingle-particle orbitalsΔl=1 couplingNilsson modelreflection asymmetryHellmann-Feynman relationnuclear structureparticle-rotor model
0
0 comments X

The pith

Δl=1 couplings between single-particle orbitals contribute substantially to octupole deformation alongside the conventional Δl=3 mode.

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

The paper challenges the long-standing emphasis on Δl=3 couplings as the main driver of octupole deformation by demonstrating that Δl=1 couplings are comparably significant near the octupole magic number N=134. Within the Nilsson model the authors compute mixing ratios of opposite-parity orbitals and introduce component-resolved energy contributions extracted via the Hellmann-Feynman relation to isolate the role of each (Δl, Δj) pair. Particle-rotor model calculations for 221Ra and 223Th then illustrate how the additional Δl=1 mixing alters rotational band structures. A sympathetic reader cares because the result implies that reflection asymmetry arises from the joint action of both couplings rather than from Δl=3 alone.

Core claim

In orbitals near N=134 the Δl=1 and Δl=3 octupole couplings between single-particle states act synergistically to produce reflection asymmetry, as shown by comparable mixing ratios in Nilsson wave functions, by matrix-element trends of the deformed potential, and by component-resolved single-particle energy shifts obtained from the Hellmann-Feynman theorem; the same couplings modify the rotational spectra of 221Ra and 223Th when treated in the particle-rotor model.

What carries the argument

Component-resolved single-particle octupole energy contributions derived from the Hellmann-Feynman relation applied to Nilsson-model wave functions, which separate the energetic role of each (Δl, Δj) coupling.

If this is right

  • Δl=1 and Δl=3 mixing ratios are of comparable magnitude for the relevant orbitals near N=134.
  • Both couplings together shape the rotational bands observed in 221Ra and 223Th.
  • A revised description of octupole correlations must incorporate the synergistic action of Δl=1 and Δl=3 terms.
  • Trends in the mixing ratios follow directly from the matrix elements of the deformed single-particle potential.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same synergy may operate in other octupole-deformed regions such as the actinides or rare-earth nuclei.
  • Collective models that currently emphasize only Δl=3 interactions could be extended to include explicit Δl=1 matrix elements.
  • Improved predictions for parity-mixing observables in heavy nuclei would follow from consistently treating both couplings.
  • Targeted experiments measuring specific inter-band transitions could quantify the relative strength of the two modes.

Load-bearing premise

That Nilsson-model wave functions and the Hellmann-Feynman energy components faithfully represent the actual mixing strengths in real nuclei near N=134 without large higher-order corrections.

What would settle it

A measurement of rotational energy levels or E3 transition strengths in 221Ra or 223Th that matches particle-rotor results obtained with Δl=1 coupling switched off but disagrees when the coupling is restored.

Figures

Figures reproduced from arXiv: 2604.11265 by Bin Qi, Chen Liu, Shouyu Wang, Xudong Wang.

Figure 1
Figure 1. Figure 1: FIG. 1. Spherical neutron single-particle levels and the dominant octupole-driving couplings. Blue [PITH_FULL_IMAGE:figures/full_fig_p020_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The neutron single-particle levels, obtained by diagonalizing the reflection-asymmetric axial [PITH_FULL_IMAGE:figures/full_fig_p021_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. The mixing ratio [PITH_FULL_IMAGE:figures/full_fig_p022_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Neutron single-particle wave function mixing ratio as functions of [PITH_FULL_IMAGE:figures/full_fig_p023_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Same plotting scheme as in Fig. 3 and Fig. 4, but with the vertical axis now showing the [PITH_FULL_IMAGE:figures/full_fig_p024_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Comparison between the PRM results (lines) and the experimental data (markers) [52] [PITH_FULL_IMAGE:figures/full_fig_p025_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Mixing ratios extracted from the PRM collective wave functions of [PITH_FULL_IMAGE:figures/full_fig_p026_7.png] view at source ↗
read the original abstract

Conventionally, octupole deformation in nuclei has been attributed to strong $\Delta l=3$ couplings between opposite-parity single-particle orbitals. In this work, we demonstrate that the often-overlooked $\Delta l=1$ mode also plays an important role. Taking orbitals near the octupole magic number $N = 134$ as a benchmark, we systematically evaluate the $\Delta l = 1$ and $\Delta l = 3$ mixing ratios of the wave functions within the Nilsson model, interpreting the trends through matrix elements of the deformed potential. We introduce component-resolved single-particle octupole energy contributions, based on the Hellmann--Feynman relation, to quantify the contributions of each $(\Delta l,\Delta j)$ coupling. Furthermore, the impact of $\Delta l = 1$ coupling on the rotational structure is demonstrated via particle-rotor model calculations for $^{221}$Ra and $^{223}$Th. Our work suggests that $\Delta l=1$ and $\Delta l=3$ octupole couplings act synergistically in driving reflection asymmetry, necessitating a revised paradigm for understanding octupole correlation.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 3 minor

Summary. The paper argues that octupole deformation in nuclei, conventionally attributed to strong Δl=3 couplings between opposite-parity single-particle orbitals, also involves important contributions from the often-overlooked Δl=1 mode. Using the Nilsson model for orbitals near the octupole magic number N=134, the authors evaluate Δl=1 and Δl=3 mixing ratios in the wave functions, interpret trends via matrix elements of the deformed potential, and introduce component-resolved single-particle octupole energy contributions derived from the Hellmann-Feynman relation. They further demonstrate the impact of Δl=1 coupling on rotational structure through particle-rotor model calculations for 221Ra and 223Th, concluding that the two couplings act synergistically to drive reflection asymmetry and that this necessitates a revised paradigm for octupole correlations.

Significance. If the model-based synergy between Δl=1 and Δl=3 couplings generalizes to realistic nuclei, the work would meaningfully revise the standard picture of octupole correlations by highlighting an additional mixing mechanism. The use of established Nilsson and particle-rotor frameworks with an interpretable Hellmann-Feynman decomposition provides clear insight into the contributions, which is a strength for pedagogical and interpretive value. However, the phenomenological nature of the approach limits the immediate impact on the broader field without further validation.

major comments (2)
  1. [Nilsson model analysis and Hellmann-Feynman section] The central claim of synergistic action between Δl=1 and Δl=3 couplings (abstract and concluding section) rests on Nilsson-model wave functions and Hellmann-Feynman energy decompositions near N=134. These quantities are derived within a phenomenological axially quadrupole-based potential extended by octupole terms, without self-consistent density rearrangements or pairing gaps. If pairing or configuration mixing alters the relative Δl=1 versus Δl=3 matrix elements, the synergistic picture would not hold beyond the model; the manuscript does not quantify this sensitivity or compare to self-consistent calculations.
  2. [Particle-rotor model calculations] In the particle-rotor model calculations for 221Ra and 223Th (rotational structure section), the demonstration of Δl=1 impact on the spectrum uses a single octupole deformation parameter without reported error estimates, sensitivity checks, or variation of other parameters such as pairing strength. This makes it difficult to assess whether the synergistic effect is robust or an artifact of the chosen parameter set, which is load-bearing for the revised-paradigm conclusion.
minor comments (3)
  1. [Abstract] The abstract states that the work 'systematically evaluate[s]' the mixing ratios but provides no numerical values, tables, or error bars; including at least one representative table of Δl=1/Δl=3 ratios for key orbitals would improve clarity.
  2. [Energy contributions section] Notation for the component-resolved energies (Hellmann-Feynman contributions) is introduced without an explicit equation defining the decomposition; adding the defining relation would aid readers.
  3. [Discussion] The manuscript would benefit from a brief discussion of how the results compare to existing experimental data on octupole moments or parity-doublet splittings in the cited nuclei.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major point below, providing clarifications on the scope and limitations of our phenomenological approach while strengthening the presentation where appropriate.

read point-by-point responses

  1. Referee: [Nilsson model analysis and Hellmann-Feynman section] The central claim of synergistic action between Δl=1 and Δl=3 couplings (abstract and concluding section) rests on Nilsson-model wave functions and Hellmann-Feynman energy decompositions near N=134. These quantities are derived within a phenomenological axially quadrupole-based potential extended by octupole terms, without self-consistent density rearrangements or pairing gaps. If pairing or configuration mixing alters the relative Δl=1 versus Δl=3 matrix elements, the synergistic picture would not hold beyond the model; the manuscript does not quantify this sensitivity or compare to self-consistent calculations.

    Authors: We agree that the Nilsson model employed is phenomenological and does not incorporate self-consistent density rearrangements or explicit pairing gaps. This choice was deliberate to isolate and interpret the single-particle mixing mechanisms driven by the octupole potential terms in a transparent manner, which is the central objective of the work. The Hellmann-Feynman decomposition provides a direct quantification of the Δl=1 and Δl=3 contributions within this framework. While pairing and configuration mixing could quantitatively modify the mixing ratios, the underlying synergy arises from the structure of the octupole operator itself and is expected to remain qualitatively robust. In the revised manuscript we will add an explicit discussion of these model limitations, including a note that self-consistent calculations (e.g., with density functional theory) would be a valuable next step to test the generality of the findings. We do not claim the results are model-independent but present them as a benchmark for such extensions. revision: partial

  2. Referee: [Particle-rotor model calculations] In the particle-rotor model calculations for 221Ra and 223Th (rotational structure section), the demonstration of Δl=1 impact on the spectrum uses a single octupole deformation parameter without reported error estimates, sensitivity checks, or variation of other parameters such as pairing strength. This makes it difficult to assess whether the synergistic effect is robust or an artifact of the chosen parameter set, which is load-bearing for the revised-paradigm conclusion.

    Authors: The particle-rotor calculations serve to demonstrate the qualitative impact of Δl=1 coupling on the rotational spectrum using established parameters for these nuclei. A single representative octupole deformation was selected to highlight the effect without introducing additional free parameters. We acknowledge that a more comprehensive sensitivity study would strengthen the presentation. In the revised manuscript we will include additional calculations varying the octupole deformation parameter around the chosen value and discuss the stability of the observed band structures. While formal error estimates are not standard in this class of models, we will clarify the parameter choices and their literature basis to allow readers to assess robustness. revision: partial

Circularity Check

0 steps flagged

No circularity: calculations derive from standard Nilsson matrix elements and Hellmann-Feynman theorem

full rationale

The derivation begins with the Nilsson deformed potential to evaluate Δl=1 and Δl=3 mixing ratios in wave functions near N=134, then applies the Hellmann-Feynman relation to obtain component-resolved energy contributions. These steps are direct evaluations of matrix elements and expectation values within the model; no parameter is fitted to the target synergy and then relabeled as a prediction, no self-definition equates inputs to outputs, and no load-bearing uniqueness theorem or ansatz is imported via self-citation. The synergistic role of the two couplings is an output of the computed trends and particle-rotor results, not presupposed by construction. The paper remains self-contained against external benchmarks for its stated scope.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The work rests on standard nuclear-physics models whose validity is assumed rather than re-derived; no new entities are postulated.

free parameters (1)
  • octupole deformation parameter
    Likely adjusted within the Nilsson model to reproduce observed trends near N=134.
axioms (2)
  • domain assumption Nilsson model provides accurate single-particle wave functions in axially deformed nuclei
    Basis for evaluating Δl=1 and Δl=3 mixing ratios.
  • standard math Hellmann-Feynman theorem can be applied component-wise to separate octupole energy contributions by (Δl,Δj)
    Used to quantify individual coupling strengths.

pith-pipeline@v0.9.0 · 5507 in / 1318 out tokens · 43401 ms · 2026-05-10T15:42:09.884901+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

52 extracted references · 52 canonical work pages

  1. [1]

    P. A. Butler and W. Nazarewicz, Rev. Mod. Phys.68, 349 (1996)

    work page 1996

  2. [2]

    P. A. Butler, J. Phys. G: Nucl. Part. Phys.43, 073002 (2016)

    work page 2016

  3. [3]

    P. A. Butler, Proc. R. Soc. A476,20200202 (2020)

    work page 2020

  4. [4]

    Nazarewicz, Nucl

    W. Nazarewicz, Nucl. Phys. A574, 27 (1994)

    work page 1994

  5. [5]

    X. D. Wang, B. Qi, C. Liu, A. Rohilla and Y. Zhang, Phys. Rev. C106, 064325 (2022)

    work page 2022

  6. [6]

    Engel, M

    J. Engel, M. J. Ramsey-Musolf, and U. van Kolck, Prog. Part. Nucl. Phys.71, 21 (2013)

    work page 2013

  7. [7]

    Bohr and B

    A. Bohr and B. R. Mottelson,Nuclear Structure(Benjamin, New York, 1975), Vol. II

    work page 1975

  8. [8]

    Frauendorf, Rev

    S. Frauendorf, Rev. Mod. Phys.73, 463 (2001)

    work page 2001

  9. [9]

    Ring and P

    P. Ring and P. Schuck,The Nuclear Many-body Problem(Springer-Verlag, Heidelberg, 1980)

    work page 1980

  10. [10]

    Bender, P.-H

    M. Bender, P.-H. Heenen, and P.-G. Reinhard, Rev. Mod. Phys.75, 121 (2003)

    work page 2003

  11. [11]

    Nazarewicz, P

    W. Nazarewicz, P. Olanders, I. Ragnarsson, J. Dudek, and G. A. Leander, Phys. Rev. Lett. 52, 1272 (1984)

    work page 1984

  12. [12]

    G. A. Leander, W. Nazarewicz, G. F. Bertsch and J. Dudek, Nucl. Phys. A453, 58 (1986)

    work page 1986

  13. [13]

    L. P. Gaffney, P. A. Butler, M. Scheck, A. B. Hayes, F. Wenander, M. Albers, B. Bastin, C. Bauer, A. Blazhev, and S. B¨oniget al., Nature (London)497, 199 (2013). 17

    work page 2013

  14. [14]

    Huang, S

    Y. Huang, S. J. Zhu, J. H. Hamilton, E. H. Wang, A. V. Ramayya, Z. G. Xiao, H. J. Li, Y. X. Luo, J. O. Rasmussen, G. M. Ter-Akopian and Y. T Oganessian, Phys. Rev. C93, 064321 (2016)

    work page 2016

  15. [15]

    V. E. Iacob, W. Urban, J. C. Bacelar, J. Jongman, J. Nyberg, G. Sletten, and L. Trache, Nucl. Phys. A596, 155 (1996)

    work page 1996

  16. [16]

    Y. J. Chen, S. J. Zhu, J. H. Hamilton, A. V. Ramayya, J. K. Hwang, M. Sakhaee, Y. X. Luo, J. O. Rasmussen, K. Li and I. Y. Leeet al., Phys. Rev. C73, 054316 (2006)

    work page 2006

  17. [17]

    P. A. Butleret al., Phys. Rev. Lett.124, 042503 (2020)

    work page 2020

  18. [18]

    J. F. C. Cockset al., Phys. Rev. Lett.78, 2920 (1997)

    work page 1997

  19. [19]

    Bucher, S

    B. Bucher, S. Zhu, C. Y. Wu, R. V. F. Janssens, D. Cline, A. B. Hayes, M. Albers, A. D. Ayangeakaa, P. A. Butler, C. M. Campbell, M. P. Carpenter, C. J. Chiara, J. A. Clark, H. L. Crawford, M. Cromaz, H. M. David, C. Dickerson, E. T. Gregor, J. Harker, C. R. Hoffman, B. P. Kay, F. G. Kondev, A. Korichi, T. Lauritsen, A. O. Macchiavelli, R. C. Pardo, A. Ri...

    work page 2016

  20. [20]

    Cline, A

    B.Bucher, S.Zhu, C.Y.Wu, R.V.F.Janssens, R.N.Bernard, L.M.Robledo, T.R.Rodríguez, D. Cline, A. B. Hayes, A. D. Ayangeakaa, M. Q. Buckner, C. M. Campbell, M. P. Carpenter, J. A. Clark, H. L. Crawford, H. M. David, C. Dickerson, J. Harker, C. R. Hoffman, B. P. Kay, F. G. Kondev, T. Lauritsen, A. O. Macchiavelli, R. C. Pardo, G. Savard, D. Seweryniak, and R....

    work page 2017

  21. [21]

    S. J. Zhu, E. H. Wang, J. H. Hamilton, A. V. Ramayya, Y. X. Liu, N. T. Brewer, Y. X. Luo, J. O. Rasmussen, Z. G. Xiao, Y. Huang, G. M. Ter-Akopian, and T. Oganessian, Phys. Rev. Lett.124, 032501 (2020)

    work page 2020

  22. [22]

    Urban, R

    W. Urban, R. M. Lieder, J. C. Bacelar, P. P. Singh, D. Alber, D. Balabanski, W. Gast, H. Grawe, G. Hebbinghaus and J. R. Jongman, Phys. Lett. B258, 293 (1991)

    work page 1991

  23. [23]

    R. K. Sheline, Phys. Rev. C37, 423 (1988)

    work page 1988

  24. [24]

    R. K. Sheline, Phys. Lett. B219, 222 (1989)

    work page 1989

  25. [25]

    C. Liu, S. Y. Wang, R. A. Bark, S. Q. Zhang, J. Meng, B. Qi, P. Jones, S. M. Wyngaardt, J. Zhao, and C. Xuet al., Phys. Rev. Lett.116, 112501 (2016)

    work page 2016

  26. [26]

    W. Z. Xuet al., Phys. Lett. B833, 137287 (2022)

    work page 2022

  27. [27]

    W. Z. Xuet al., Phys. Rev. C109, 044303 (2024). 18

    work page 2024

  28. [28]

    Xiaoet al., Phys

    X. Xiaoet al., Phys. Rev. C106, 064302 (2022)

    work page 2022

  29. [29]

    Fernández-Niello, C

    J. Fernández-Niello, C. Mittag, F. Riess, E. Ruchowska, and M. Stallknecht, Nucl. Phys. A 531, 164 (1991)

    work page 1991

  30. [30]

    R. K. Sheline, Phys. Lett. B166, 269 (1986)

    work page 1986

  31. [31]

    C. W. Reich, I. Ahmad, and G. A. Leander, Phys. Lett. B169, 148 (1986)

    work page 1986

  32. [32]

    Dahlinger, E

    M. Dahlinger, E. Kankeleit, D. Habs, D. Schwalm, B. Schwartz, R. S. Simon, J. D. Burrows, and P. A. Butler, Nucl. Phys. A484, 337 (1988)

    work page 1988

  33. [33]

    Morse, A

    C. Morse, A. O. Macchiavelli, H. L. Crawford et al., Phys. Rev. C102, 054328 (2020)

    work page 2020

  34. [34]

    S. J. Zhu, J. H. Hamilton, A. V. Ramayya et al., Phys. Rev. C65, 014307 (2001)

    work page 2001

  35. [35]

    Urban, W

    W. Urban, W. R. Phillips, J. L. Durell, M. A. Jones, M. Leddy, C. J. Pearson, A. G. Smith, B. J. Varley, I. Ahmad, L. R. Morss, M. Bentaleb, E. Lubkiewicz, and N. Schulz, Phys. Rev. C54, 945 (1996)

    work page 1996

  36. [36]

    S. J. Zhu, J. H. Hamilton, A. V. Ramayya, Y. X. Luo, J. O. Rasmussen, Z. G. Xiao, Y. Huang, G. M. Ter-Akopian, and T. Oganessian, Phys. Rev. C59, 1316 (1999)

    work page 1999

  37. [37]

    X. C. Han, C. Liu, and S. Y. Wang, Int. J. Mod. Phys. E32, 2340003 (2023)

    work page 2023

  38. [38]

    S. E. Agbemava, A. V. Afanasjev and P. Ring, Phys. Rev. C93, 044304 (2016)

    work page 2016

  39. [39]

    Rodríguez-Guzmán R, L. M. Robledo, K. Nomura and N. C. Hernandez, J. Phys. G: Nucl. Part. Phys.49, 015101 (2021)

    work page 2021

  40. [40]

    J. J. Sakurai and J. Napolitano,Modern Quantum Mechanics(Addison-Wesley, Boston, 2011)

    work page 2011

  41. [41]

    S. G. Nilsson, Mat. Fys. Medd. Dan. Vid. Selsk.29, No. 16 (1955)

    work page 1955

  42. [42]

    S. G. Nilsson, C. F. Tsang, A. Sobiczewski, Z. Szymański, S. Wycech, C. Gustafson, I.-L. Lamm, P. Møller, and B. Nilsson, Nucl. Phys. A131, 1 (1969)

    work page 1969

  43. [43]

    Y. Y. Wang, S. Q. Zhang, P. W. Zhao and J. Meng, Phys. Lett. B792, 454 (2019)

    work page 2019

  44. [44]

    G. A. Leander and R.K. Sheline, Nucl. Phys. A413, 375 (1984)

    work page 1984

  45. [45]

    G. A. Leander and Y.D. Chen, Phys. Rev. C35, 1145 (1987)

    work page 1987

  46. [46]

    G. A. Leander and Y. S. Chen, Phys. Rev. C37, 2744 (1988)

    work page 1988

  47. [47]

    Bengtsson and I

    T. Bengtsson and I. Ragnarsson, Nucl. Phys. A436, 14 (1985)

    work page 1985

  48. [48]

    Hamamoto, B

    I. Hamamoto, B. Mottelson, H. Xie, and X.Z. Zhang, Z. Phys. D21163 (1991)

    work page 1991

  49. [49]

    Zhang, B

    H. Zhang, B. Qi, X. D. Wang, H. Jia, and S. Y. Wang, Phys. Rev. C105, 034339 (2022)

    work page 2022

  50. [50]

    Strutinsky, At

    V. Strutinsky, At. Energ.4, 150 (1956); J. Nucl. Eng.4, 523 (1957)

    work page 1956

  51. [51]

    D. D. Zhang, Int. J. Mod. Phys. E32, 2340009 (2023). 19

    work page 2023

  52. [52]

    A. Jain, S. Singh, S. Kumar, and J. K. Tuli, Nucl. Data Sheets108, 883 (2007). Δ l=1d3 /2s1 /2g7 /2d5 /2j1 5/2i11/2g9 /2p1 /2f5 /2p3 /2i13/2h9 /2f7 /2d 3 /2h1 1/2s1 /2g7 /2d5 /24 g9 /23p1 /23f5 /23p3 /25 08 21 261 84Y 3 Y 3 Y 3 Y 3 1348 85 63 4Δ l=3 FIG. 1. Spherical neutron single-particle levels and the dominant octupole-driving couplings. Blue (red) li...

    work page 2007