Search-E1: Nuclear PDR and Bayesian Modeling

• Search-E1 is a framework that decomposes low-energy E1 excitations in neutron-rich nuclei into skin-mode and pn-mode, clarifying the physical nature of these oscillations.
• The method employs the Random Phase Approximation and transition density analysis to quantify the crossover from pure neutron-skin to pn-mode dominance as energy increases.
• Its integration of Bayesian search techniques and nuclear modeling provides actionable insights into symmetry energy correlations and predictive failure patterns.

The term "Search-E1" encompasses several highly technical concepts across nuclear, atomic, condensed matter, computational, and machine learning domains. Prominently, it refers to both advanced nuclear structure phenomena—most notably the decomposition of low-energy electric dipole ($E1$) excitations in neutron-rich nuclei—as well as cutting-edge innovations in search-augmented reasoning and Bayesian modeling for the analysis of component failures. The following sections provide a comprehensive, technically rigorous exploration of all primary usages, with a focus on the quantifiable decomposition of $E1$ nuclear modes as formulated by Nakada, Inakura, and Sawai (Nakada et al., 2012), and their methodological and physical implications.

1. $E1$ Excitations in Neutron-Rich Nuclei: Mode Decomposition and Physical Picture

Electric dipole ($E1$) excitations in neutron-rich nuclei exhibit a low-energy component below the giant dipole resonance (GDR), often referred to as the "pygmy dipole resonance" (PDR). These excitations manifest as collective oscillations involving both the neutron-rich surface ("skin") and the isovector proton–neutron core oscillation. The Random Phase Approximation (RPA), built atop a Hartree–Fock (HF) ground state parametrized by effective interactions (Skyrme, Gogny, M3Y, etc.), provides the microscopic framework for their analysis.

From the HF ground state, small-amplitude 1p–1h oscillations are treated in RPA, which yields a spectrum of $1^-$ eigenmodes $|\alpha\rangle$ with discrete excitation energies $\omega_\alpha$ and transition amplitudes $\langle \alpha|\hat O^{(E1)}|0\rangle$, where the operator includes the center-of-mass correction: $$S^{(E1)}(\omega) = \frac{\gamma}{\pi} \sum_\alpha \left[\frac{1}{(\omega-\omega_\alpha)^2 + \gamma^2} - \frac{1}{(\omega+\omega_\alpha)^2 + \gamma^2}\right] |\langle\alpha|\hat O^{(E1)}|0\rangle|^2$$ with typical smearing width $2\gamma = 1$ MeV.

2. Transition Density Formulation and Decomposition Criterion

For each $E1$0 eigenstate $E1$1, the proton and neutron transition densities,

$E1$2

$E1$3

are combined into an isovector E1 transition density with exact center-of-mass correction: $E1$4 By construction, $E1$5 for each state, ensuring absence of spurious c.m. admixture.

The decomposition into "pn-mode" (core proton–neutron oscillation) and "skin-mode" (neutron-skin against core) components proceeds as follows:

• At each $E1$6, if $E1$7 ($E1$8, typically 0.05), classify as "skin-like"; otherwise, as pn-like.
• The E1 transition density is then split:
  • $E1$9 in "skin-like" regions, zero elsewhere.
  • $E1$0 in pn-like regions, zero elsewhere.
• The matrix element decomposes:

  $E1$1

  where $E1$2 for $E1$3.

This methodology enables a transparent, quantitative partitioning of each RPA state into neutron-skin and pn-mode oscillation content.

3. Strength Function Decomposition and Mixing Ratios

Each excitation contributes partial strengths to the pn- and skin-modes, as well as an interference term: $E1$4 with partial strengths defined by replacing $E1$5 in the summation by $E1$6, $E1$7, and $E1$8 respectively.

The normalized fractions are: $E1$9 with $E1$0.

4. Universal Crossover Phenomenon: Skin-PN Mixing and Energy Dependence

A key empirical result is the universal crossover behavior in neutron-rich, (near-)doubly-magic nuclei—from $E1$1O and $E1$2O up to $E1$3Ca, $E1$4Zr, and $E1$5Sn—across a variety of effective interactions (SkM*, SkI2, M3Y-P7, …):

• At the lowest $E1$6 energies, $E1$7: almost pure neutron-skin oscillations.
• As $E1$8 increases, $E1$9 decreases monotonically; $1^-$0 increases, crossing at $1^-$1 at $1^-$2 MeV.
• Near the GDR peak, $1^-$3, and the E1 strength is dominated by pn-mode motion.

This crossover energy is remarkably insensitive to mass number (between $1^-$4Ca and $1^-$5Sn) and the choice of effective interaction. For $1^-$6Zr, even the observable pygmy bump at $1^-$7 MeV is not of pure skin-mode character—the decomposition reveals substantial pn-mode admixture, and the corresponding RPA state occurs near 13 MeV in SkI2-based calculation (Nakada et al., 2012).

5. Transition Density Structure and Experimental Discrimination

The identification of the dominant oscillation type at each energy is informed directly by the computed transition densities:

• In the pn-mode, $1^-$8 and $1^-$9 are everywhere out of phase ($|\alpha\rangle$0), representing a collective core proton–neutron oscillation.
• In the skin-mode, proton and neutron transition densities are in phase in the core (interior $|\alpha\rangle$1), but the neutron density protrudes at the nuclear surface, producing a characteristic "neutron-skin vs. core" pattern.

Experimentally, the separation of pn- and skin-mode contributions in the PDR region (low-energy E1) is only possible through an analysis of transition densities—e.g., via comparative studies of $|\alpha\rangle$2 and $|\alpha\rangle$3 cross sections, which have different sensitivity to isoscalar (skin-like) and isovector (pn-mode) oscillations.

6. Implications for Nuclear Structure and Symmetry Energy

The realization that the low-energy E1 strength (PDR) is always a mixture of neutron-skin and pn-mode oscillations, with the degree of mixing controlled mainly by excitation energy:

• Validates but also qualifies the macroscopic picture of the PDR as a pure neutron-skin vibration.
• Demonstrates that below $|\alpha\rangle$4 MeV, skin-mode content dominates; above, the classical proton–neutron GDR emerges.
• Provides a systematic, model-independent way to relate low-energy E1 strength to neutron-skin thickness and the symmetry-energy sector of the nuclear equation of state.

Consequently, measurement and interpretation of the PDR region must account for this continuous, energy-dependent crossover. The observed "pygmy" strengths in stable and moderately neutron-rich nuclei are not exclusively skin-mode in character; only direct analysis of transition densities in both experiment and theory can disentangle the true mode content (Nakada et al., 2012).

7. Broader Methodological Significance and Further Research Directions

The Search-E1 decomposition propounded by Nakada et al. is highly general and robust:

• Universality of the crossover energy and mode-fraction curves has been confirmed for a range of doubly-magic nuclei, effective forces, and nuclear masses.
• The technique is independent of the specific ground-state interaction, provided self-consistent HF+RPA is used.
• The framework enables re-examination of the physical nature of the PDR and its correlation with neutron-skin thickness and the density dependence of the symmetry energy.
• Further extensions could combine this decomposition method with experimental probes explicitly sensitive to spatial current patterns, as well as Energy Density Functional (EDF) calculations in open-shell or deformed nuclei.

This paradigm fundamentally constrains how low-energy dipole collectivity and its observables are connected to isovector nuclear matter properties, and supplies an essential reference point for both ab initio and phenomenological modeling of exotic nuclear systems (Nakada et al., 2012).