Origin of the electric hexadecapole isomer in $^{93}$Mo

Basic Information

• Paper ID: 2510.22352
• Title: Origin of the electric hexadecapole isomer in $^{93}$Mo
• Authors: B. Maheshwari (GANIL), P. Van Isacker (GANIL), P. M. Walker (University of Surrey)
• Classification: nucl-th (nuclear theory), nucl-ex (nuclear experiment)
• Submission Date: October 25, 2025
• Paper Link: https://arxiv.org/abs/2510.22352

Abstract

This paper investigates the microscopic origin of the anomalous $21/2^+$ isomeric state in $^{93}$Mo through shell model analysis. This state is an important candidate for nuclear electron capture excitation (NEEC) processes. Its unique characteristic lies in the fact that decay is dominated by slow electric hexadecapole (E4) transitions, while the typically faster electric quadrupole (E2) decay pathway is energetically forbidden. The study reveals the microscopic mechanism through detailed analysis of initial and final state wave function structures and E4 transition matrix elements, and compares with its particle-hole conjugate nucleus $^{99}$Cd (which lacks E4 transitions).

Research Background and Motivation

Research Questions

1. Core Physics Problem: Why does the $21/2^+$ isomeric state in $^{93}$Mo decay through the rare E4 transition rather than the common E2 transition? What is the microscopic nuclear structure origin of this phenomenon?
2. Practical Application Context: This isomeric state is the only nuclear state reported to exhibit the NEEC phenomenon. The NEEC process involves using a 5 keV trigger to induce transitions from the long-lived $21/2^+$ state (half-life 6.85 hours) to the short-lived $17/2^+$ state (half-life 3.53 nanoseconds), potentially achieving at least 50-fold energy gain.

Importance of the Problem

1. Fundamental Physics Significance: High-multipolarity transitions are extremely sensitive to radial and angular overlaps of wave functions, making them excellent probes for nuclear structure
2. Energy Application Prospects: Understanding the NEEC mechanism is crucial for developing potential nuclear energy storage and release technologies
3. Isomeric State Physics: The formation mechanism of long-lived excited states is a fundamental problem in nuclear physics

Limitations of Previous Research

1. Large-scale shell model calculations, while reproducing experimental features, lack transparent quantitative analysis
2. Lack of clear explanation of the microscopic mechanism of E4 transitions themselves (wave function components, coupling modes)
3. Previous estimates of the critical $B(E2; 17/2^+ → 21/2^+)$ value may contain deviations

Research Motivation

Through detailed analysis in a simplified shell model space ($\nu1d_{5/2} \otimes \pi0g^2_{9/2}$), to explicitly quantify the role of neutron-proton interactions in level inversion, and verify the universality of the mechanism through comparison with $^{99}$Cd.

Core Contributions

1. Revealing the Level Inversion Mechanism: Quantitatively demonstrates that the level inversion of $21/2^+$ and $17/2^+$ states originates from strong attractive neutron-proton (νπ) interactions in the maximum spin alignment configuration
2. Microscopic Decomposition of E4 Transitions: For the first time, explicitly decomposes the various contributing components of the E4 transition matrix element, finding that it is primarily contributed by the proton two-body matrix element $\pi0g^2_{9/2}(J_\pi=8) → \pi0g^2_{9/2}(J_\pi=4)$, with neutron and proton parts exhibiting constructive interference
3. Correction of Critical NEEC Parameters: Reduces the theoretical value of $B(E2)$ for $17/2^+ → 21/2^+$ by approximately 40% compared to previous estimates, which is significant for NEEC process assessment
4. Particle-Hole Comparison Verification: Through systematic comparison with $^{99}$Cd, demonstrates that under Pandya transformation, the νπ interaction in particle-hole configurations becomes repulsive, leading to no level inversion and no E4 isomeric state, thereby verifying the correctness of the proposed mechanism
5. Detailed Wave Function Analysis: Accurately calculates the wave function composition of each energy level, with the lowest $13/2^+$ state approximately 80% composed of $J_\pi=4$ proton coupling, and the $17/2^+$ state primarily $J_\pi=6$ coupling

Detailed Methodology

Task Definition

Input: Nuclear shell model Hamiltonian parameters (single-particle energies, two-body interaction matrix elements)
Output: Energy level structure, wave function composition, electromagnetic transition probabilities B(EL)
Constraints: Limited to the $\nu1d_{5/2} \otimes \pi0g^2_{9/2}$ configuration space

Theoretical Framework

1. Shell Model Hamiltonian

Standard form is adopted: $\hat{H} = \sum_i \epsilon_i \hat{n}_i + \frac{1}{4}\sum_{ijkl} V_{ijkl} a^\dagger_i a^\dagger_j a_l a_k + \cdots$

Parameter Determination Strategy:

• Like-particle interactions $V_{\pi\pi}$: Empirically determined from the experimental energy spectrum of $^{92}$Mo (two valence protons in the $0g_{9/2}$ orbit)
  • $V_{\pi\pi}(J=0)=0.000$, $V_{\pi\pi}(J=2)=1.510$, $V_{\pi\pi}(J=4)=2.283$, $V_{\pi\pi}(J=6)=2.612$, $V_{\pi\pi}(J=8)=2.761$ MeV
• Neutron-proton interaction $V_{\nu\pi}$: Fixed from the energy spectrum of $^{92}$Nb ($\nu1d_{5/2}$ coupled with $\pi0g_{9/2}$)
  • $V_{\nu\pi}(J=2)=0.135$, $V_{\nu\pi}(J=3)=0.286$, $V_{\nu\pi}(J=4)=0.480$, $V_{\nu\pi}(J=5)=0.357$, $V_{\nu\pi}(J=6)=0.501$, $V_{\nu\pi}(J=7)=0.000$ MeV
  • Key feature: The maximum alignment state $J=7$ exhibits the strongest attractive interaction (set as reference point)

2. Electromagnetic Transition Operators

Reduced transition probability: $B(EL;\alpha_i J_i \to \alpha_f J_f) = \frac{|\langle \alpha_f J_f ||\hat{O}(EL)||\alpha_i J_i\rangle|^2}{2J_i+1}$

Electric transition operator: $\hat{O}(ELM) = \sum_k e_k r^L_k Y^M_L(\theta_k, \phi_k)$

Many-body state expansion: $|\alpha_i J_i M_i\rangle = \sum_k a_k(\alpha_i J_i)|kJ_i M_i\rangle$

Leading to generalized transition probability: $B(EL;\alpha_i J_i \to \alpha_f J_f) = \frac{1}{2J_i+1}\left|\sum_{kl} a_k(\alpha_i J_i)b_l(\alpha_f J_f)\langle lJ_f||\hat{O}(EL)||kJ_i\rangle\right|^2$

This expression explains the interference effects of neutron and proton contributions.

3. Single-Particle Reduced Matrix Elements

For E4 transitions: $\langle \nu1d_{5/2}||\hat{O}_\nu(E4)||\nu1d_{5/2}\rangle = \frac{513}{4}\sqrt{\frac{1}{7\pi}}e_\nu b^4$$\langle \pi0g_{9/2}||\hat{O}_\pi(E4)||\pi0g_{9/2}\rangle = \frac{9}{4}\sqrt{\frac{715}{\pi}}e_\pi b^4$

where $b$ is the harmonic oscillator length parameter: $b^2 = 41.46/(45A^{-1/3} - 25A^{-2/3})$ fm²

Technical Innovations

1. Validity of Simplified Model Space: Although considering only one neutron orbit and one proton orbit, by carefully selecting two-body interaction parameters, the model accurately reproduces experimental energy spectra and transition strengths
2. Systematic Parameter Determination: Rather than fitting $^{93}$Mo itself, parameters are determined from experimental data of neighboring nuclei, ensuring reliability of predictions
3. Particle-Hole Symmetry Analysis: Utilizes Pandya transformation to establish theoretical connections between $^{93}$Mo and $^{99}$Cd, verifying the mechanism through comparison
4. Complete Decomposition of Transition Matrix Elements: Explicitly provides contributions of various wave function components to transition strength and their phase relationships

Experimental Setup

Computational Tools

• Mathematica code shell.m used for calculations in J-scheme
• Coefficients of fractional parentage method employed
• Numerical diagonalization required to solve eigenvalue problems

Effective Charge Selection

E2 Transitions:

• Proton effective charge: $e_\pi = 1.32$ (standard value)
• Neutron effective charge: Three values adopted
  • $e_\nu = 1.18$ (determined from $Q(5/2^+)$ quadrupole moment of $^{91}$Zr)
  • $e_\nu = 1.73$ (determined from $B(E2; 2^+_1 \to 0^+_1)$ of $^{92}$Zr)
  • $e_\nu = 1.48$ (average value)

E4 Transitions:

• Proton effective charge: $e_\pi = 1.20$ (determined from measured B(E4) value of $^{92}$Mo)
• Neutron effective charge: Same values as E2 operator

Comparison Method

Primarily compared with large-scale shell model calculation results from Hasegawa et al. (2011)

Experimental Results

Main Results

1. Energy Level Structure (Table 1)

Key Finding: Successfully reproduces the level inversion where the $21/2^+$ state (2.437 MeV) lies below the $17/2^+$ state (2.454 MeV), with an error of only 13 keV, quite close to the experimental energy difference of 5 keV.

2. Electromagnetic Transition Probabilities (Table 2, units: W.u.)

Important Results:

• E4 transition strength in excellent agreement with experiment (1.37 W.u. vs 1.449 W.u., using average effective charge)
• $B(E2)$ value for $17/2^+ \to 21/2^+$ approximately 40% lower than Hasegawa's estimate (2.21 vs 3.5 W.u.)
• E2 transition $17/2^+ \to 13/2^+$ as critical test of wave function accuracy, theoretical value 3.98 W.u. agrees well with experimental 4.48 W.u.

Wave Function Analysis (Figure 2)

$13/2^+_1$ State Composition:

• $J_\pi=4$: ~80% (dominant component)
• $J_\pi=6$: ~15%
• $J_\pi=8$: ~5%

$17/2^+_1$ State Composition:

• $J_\pi=6$: ~85% (dominant component, lower alignment)
• $J_\pi=4$ and $J_\pi=8$: ~7-8% each

$21/2^+$ State:

• 100% pure $J_\pi=8$ configuration (maximum alignment)

E4 Transition Matrix Element Decomposition (Figure 5)

Proton contributions (units: $b^4$):

• $\pi0g^2_{9/2}(J_\pi=8) \to \pi0g^2_{9/2}(J_\pi=4)$: Maximum contribution (~3.0)
• $\pi0g^2_{9/2}(J_\pi=8) \to \pi0g^2_{9/2}(J_\pi=6)$: Moderate contribution (~1.5)
• $\pi0g^2_{9/2}(J_\pi=8) \to \pi0g^2_{9/2}(J_\pi=8)$: Small contribution with opposite phase (~-0.5)

Neutron contribution: $1d_{5/2} \to 1d_{5/2}$, matrix element 23.32 $b^4$

Total Transition Probability: $B(E4; 21/2^+ \to 13/2^+) = (0.83e_\nu + 3.53e_\pi)^2 b^8$

Neutron and proton contributions have the same sign, exhibiting constructive interference, which enhances the transition strength.

Comparison with $^{99}$Cd (Figure 1)

$^{93}$Mo:

• $21/2^+$ state (2.424 MeV) < $17/2^+$ state (2.429 MeV) → E2 forbidden → E4 isomeric state
• νπ interaction: Particle-particle, attractive at maximum alignment

$^{99}$Cd:

• $21/2^+$ state far above $17/2^+$ state → E2 allowed → No E4 isomeric state
• νπ interaction: Particle-hole (relative to $^{100}$Sn core), repulsive at maximum alignment (Pandya transformation)

This comparison perfectly verifies the mechanism that the nature of νπ interaction determines the level sequence.

Physical Picture Analysis

Angular Momentum Characteristics (Figure 3)

Displays the three-dimensional form of the fourth-order spherical harmonic $Y^M_4(\theta,\phi)$ (M=-4 to +4):

• Complex nodal structure and angular symmetry
• Nucleons must undergo significant orbital rearrangement to match these angular distributions
• Red/blue represents positive/negative phases

Radial Overlap (Figure 4)

Compares normalized radial integrated functions for E2 and E4 transitions:

• The $r^4$ dependence of E4 makes the transition highly sensitive to nuclear surface wave function behavior
• E4 contribution peaks at larger radii (~4-5 fm), deviating from the most probable position of protons inside the nucleus (~2-3 fm)
• The $0g_{9/2}$ orbit is more compact than $1d_{5/2}$, but both belong to the same harmonic oscillator shell

Radial Integral Values:

• $\langle \nu1d_{5/2}|r^4|\nu1d_{5/2}\rangle = 42.75 b^4$
• $\langle \pi0g_{9/2}|r^4|\pi0g_{9/2}\rangle = 35.75 b^4$

The realization of E4 transitions requires both complex angular rearrangement and radial overlap to be simultaneously satisfied.

Related Work

Nuclear Isomeric State Research

1. Classical Theory: Auerbach and Talmi (1964) first proposed the qualitative explanation that νπ interaction in $^{93}$Mo leads to level inversion
2. Large-Scale Calculations: Hasegawa et al. (2011) reproduced spectral features using shell model with larger configuration space
3. This Work's Improvement: Provides transparent quantitative analysis and corrected transition strengths within simplified space

NEEC Research

1. Experimental Report: Chiara et al. (2018, Nature) first reported NEEC phenomenon in $^{93}$Mo
2. Controversial Verification: Guo et al. (2022, PRL) questioned this result, still under active investigation
3. Application Prospects: Carroll and Chiara (2024) discussed energy applications of NEEC

High-Multipolarity Transitions

• High-multipolarity transitions are rare due to strong suppression by Weisskopf estimates
• E4 transitions are extremely sensitive to nuclear structure, serving as stringent tests of theoretical models

Particle-Hole Symmetry

• Pandya transformation (1956) establishes the relationship between particle-particle and particle-hole interactions
• This work represents a successful application of this symmetry in practical nuclear structure problems

Conclusions and Discussion

Main Conclusions

1. Level Inversion Mechanism Clarified: The inversion of $21/2^+$ and $17/2^+$ states in $^{93}$Mo originates entirely from the strong attractive νπ interaction at maximum spin alignment (J=7) in the $\nu1d_{5/2} \otimes \pi0g^2_{9/2}$ configuration
2. Origin of E4 Isomeric State: Level inversion makes the fast E2 decay pathway energetically forbidden, forcing the $21/2^+$ state to decay through slow E4 transition to the $13/2^+$ state
3. Quantification of Transition Mechanism: E4 transition is dominated by proton two-body matrix element ($J_\pi=8 \to 4$ contribution largest), with neutron and proton parts exhibiting constructive interference, total strength $(0.83e_\nu + 3.53e_\pi)^2 b^8$
4. NEEC Parameter Correction: The critical $B(E2; 17/2^+ \to 21/2^+)$ value is approximately 40% lower than previous estimates, significantly impacting NEEC process assessment
5. Mechanism Universality Verification: In $^{99}$Cd, particle-hole configuration makes νπ interaction repulsive, with no level inversion and no E4 isomeric state, perfectly verifying the proposed mechanism

Limitations

1. Model Space Constraints: Only considers $\nu1d_{5/2} \otimes \pi0g^2_{9/2}$ configuration, neglecting possible contributions from other orbits (such as $\pi0g_{7/2}$, $\nu0g_{7/2}$, etc.)
2. Effective Charge Uncertainty: Neutron effective charge varies in the range 1.18-1.78, and while having smaller impact on E4, still introduces systematic error
3. Energy Difference Precision: Calculated $21/2^+-17/2^+$ energy difference is 13 keV, while experiment shows only 5 keV, with relatively large relative error (though absolute error is small)
4. Three-Body and Higher Interactions: Shell model truncation at two-body level may miss many-body correlation effects
5. Core Excitations Neglected: Assumes inert $^{88}$Sr core, while actual core polarization and core excitations may contribute

Future Directions

1. M4 Isomeric States: Proposes that M4 isomeric state in $^{84}$Rb requires similar investigation (mentioned as ongoing in the paper)
2. More Candidate Nuclei: Search for similar high-multipolarity isomeric states in other regions of the nuclear chart
3. NEEC Experimental Verification: Provide more precise nuclear structure input for ongoing NEEC experimental controversies
4. Isomeric State Depletion: Study possible pathways to induce isomeric state depletion
5. Extended Model Space: Verify conclusions of simplified model in larger configuration space

In-Depth Evaluation

Strengths

1. Clear Physical Picture: Through simplified model space and systematic parameter determination, the complex nuclear structure problem is reduced to a clear physical mechanism—strong attractive νπ interaction in maximum alignment configuration
2. Thorough Quantitative Analysis:
   • Explicitly provides percentages of various wave function components
   • Complete decomposition of E4 transition matrix element contributions
   • Quantifies interference effects between neutron and proton parts
3. Excellent Experimental Agreement:
   • Energy level structure error <100 keV
   • E4 transition strength agrees with experiment within 10%
   • Critical E2 transition also shows good agreement
4. Powerful Comparative Analysis: Systematic comparison between $^{93}$Mo and $^{99}$Cd through particle-hole symmetry verifies the mechanism, enhancing credibility of conclusions
5. Excellent Visualization:
   • Figure 3 displays three-dimensional form of spherical harmonics, intuitively understanding angular momentum characteristics
   • Figure 4's radial distribution comparison reveals surface sensitivity of E4 transitions
   • Figure 5's matrix element decomposition clearly shows contributions of each component
6. High Practical Value: Provides corrected critical parameters for NEEC research, offering direct guidance for energy application studies

Shortcomings

1. Insufficient Model Space Verification: While simplified model gives good results, lacks systematic investigation of conclusion stability when including more orbits
2. Lacking Uncertainty Analysis:
   • No theoretical error estimates provided for calculated results
   • While effective charge selection effects are discussed, analysis is not systematic enough
   • Error propagation from parameter fitting not quantified
3. Simplified $^{99}$Cd Analysis: While mentioning ability to reproduce its energy level structure, lacks detailed data and figures; comparative analysis could be more thorough
4. Missing Core Effect Discussion: Completely neglects core polarization and core excitations, without assessing their possible impact
5. E2 Forbiddenness Degree: While level inversion causes E2 "energy forbiddenness," actual energy difference is only 5 keV, so this transition, though slow, is not completely forbidden; discussion of this detail insufficient
6. Limited Experimental Predictions: Beyond correcting B(E2) value, lacks new predictions for experimental verification

Impact

1. Academic Contribution:
   • Provides in-depth analysis of typical case for nuclear isomeric state physics
   • Demonstrates effectiveness of simplified shell model in understanding complex nuclear structure
   • Quantifies crucial role of νπ interaction in level structure
2. Practical Value:
   • Directly serves resolution of NEEC experimental controversy
   • Provides theoretical support for new concepts in nuclear energy storage and release
   • Can guide search for other high-multipolarity isomeric states
3. Reproducibility:
   • Model space and parameters completely explicit
   • Computational methods standard and mature
   • Results easily verifiable through repetition
4. Potential Applications:
   • Nuclear clock development (utilizing long-lived isomeric states)
   • Nuclear battery technology
   • Refinement of fundamental nuclear data

Applicable Scenarios

1. Theoretical Research:
   • Studies requiring understanding of microscopic origin of nuclear isomeric states
   • Quantitative studies of νπ interaction effects
   • Exploration of high-multipolarity transition mechanisms
2. Experimental Guidance:
   • Experimental design and data analysis for NEEC processes
   • Selection of candidate nuclei for isomeric state search
   • Theoretical support for precision nuclear spectroscopy experiments
3. Methodological Reference:
   • Strategy for effective use of simplified model space
   • Method for determining interactions from neighboring nuclei data
   • Application of particle-hole symmetry in nuclear structure analysis
4. Teaching Case:
   • Typical example of shell model application
   • Instance of nuclear level structure and transition theory
   • Example of theory-experiment integration

References

Key references cited in this paper include:

1. Isomeric State Review: Walker, Jain, Maheshwari (2024) - Nuclear isomeric states special topic
2. First NEEC Report: Chiara et al., Nature 554, 216 (2018)
3. NEEC Controversy: Guo et al., PRL 128, 242502 (2022)
4. Classical Theory: Auerbach & Talmi, Phys. Lett. 9, 153 (1964)
5. Previous Calculations: Hasegawa et al., Phys. Lett. B 696, 197 (2011)
6. Authors' Previous Work: Maheshwari et al., Phys. Rev. C 112, 024308 (2025)
7. Nuclear Structure Theory: Talmi, "Simple Models of Complex Nuclei" (1993)
8. Pandya Transformation: Pandya, Phys. Rev. 103, 956 (1956)

Overall Assessment: This is a high-quality nuclear theory research paper that, through simplified yet elegant shell model analysis, deeply reveals the microscopic origin of the rare E4 isomeric state in $^{93}$Mo. The physical picture is clear, quantitative analysis thorough, experimental agreement excellent, and the mechanism verified through particle-hole comparison. The paper not only solves fundamental nuclear structure problems but also provides important corrected parameters for the cutting-edge NEEC research. Main shortcomings lie in systematic verification of model space and uncertainty analysis, which could be further strengthened. This work has significant academic value and practical significance for nuclear isomeric state physics and NEEC research.