---
title: 'Search-E1: Nuclear PDR and Bayesian Modeling'
url: https://www.emergentmind.com/topics/search-e1
type: topic
---

# Search-E1: Nuclear PDR and Bayesian Modeling

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 [1211.0057], 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 $1^-$ eigenstate $|\alpha\rangle$, the proton and neutron transition densities,
\[
r^2\,\delta\rho_p^{(1)}(r;\alpha) = \langle\alpha|\sum_{i\in p}\delta(r-r_i) r_i Y^{(1)}(\hat r_i)|0\rangle,
\]
\[
r^2\,\delta\rho_n^{(1)}(r;\alpha) = \langle\alpha|\sum_{i\in n}\delta(r-r_i) r_i Y^{(1)}(\hat r_i)|0\rangle,
\]
are combined into an isovector E1 transition density with exact center-of-mass correction:
\[
\delta\rho^{(E1)}(r;\alpha) = \frac{N}{A}\delta\rho_p^{(1)}(r;\alpha) - \frac{Z}{A}\delta\rho_n^{(1)}(r;\alpha).
\]
By construction, $\int r^2 dr\,[\delta\rho_p^{(1)}(r;\alpha)+\delta\rho_n^{(1)}(r;\alpha)] = 0$ 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 $r$, if $\delta\rho_p^{(1)}(r;\alpha)/\delta\rho_n^{(1)}(r;\alpha) > -\lambda_s$ ($\lambda_s \ll 1$, typically 0.05), classify as "skin-like"; otherwise, as pn-like.
- The E1 transition density is then split:
    - $\delta\rho_{\rm skin}^{(E1)}(r;\alpha) = \delta\rho^{(E1)}(r;\alpha)$ in "skin-like" regions, zero elsewhere.
    - $\delta\rho_{\rm pn}^{(E1)}(r;\alpha) = \delta\rho^{(E1)}(r;\alpha)$ in pn-like regions, zero elsewhere.
- The matrix element decomposes:
    \[
    \langle\alpha|\hat O^{(E1)}|0\rangle = M_{\rm pn}(\alpha) + M_{\rm skin}(\alpha)
    \]
    where $M_X(\alpha) = \int r^2 dr\, \delta\rho_X^{(E1)}(r;\alpha)$ for $X = \rm pn, skin$.

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:
\[
S^{(E1)}(\omega) = S_{\rm pn}(\omega) + S_{\rm skin}(\omega) + S_{\rm intf}(\omega)
\]
with partial strengths defined by replacing $|\langle\alpha|\hat O|0\rangle|^2$ in the summation by $|M_{\rm pn}(\alpha)|^2$, $|M_{\rm skin}(\alpha)|^2$, and $2\Re[M_{\rm pn}^* M_{\rm skin}]$ respectively.

The normalized fractions are:
\[
\mathcal R_{\rm pn}(\omega) = \frac{S_{\rm pn}(\omega)}{S^{(E1)}(\omega)},\quad
\mathcal R_{\rm skin}(\omega) = \frac{S_{\rm skin}(\omega)}{S^{(E1)}(\omega)},\quad
\mathcal R_{\rm intf}(\omega) = \frac{S_{\rm intf}(\omega)}{S^{(E1)}(\omega)},
\]
with $\mathcal R_{\rm pn}+\mathcal R_{\rm skin}+\mathcal R_{\rm intf}=1$.

## 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 $^{22}$O and $^{24}$O up to $^{52}$Ca, $^{90}$Zr, and $^{132}$Sn—across a variety of effective interactions (SkM*, SkI2, M3Y-P7, …):

- At the lowest $E1$ energies, $\mathcal R_{\rm skin}(\omega)\approx1$: almost pure neutron-skin oscillations.
- As $\omega$ increases, $\mathcal R_{\rm skin}(\omega)$ decreases monotonically; $\mathcal R_{\rm pn}(\omega)$ increases, crossing at $\mathcal R_{\rm skin}=\mathcal R_{\rm pn}=0.5$ at $\omega\approx10$ MeV.
- Near the GDR peak, $\mathcal R_{\rm skin}(\omega)\rightarrow0$, and the E1 strength is dominated by pn-mode motion.

This crossover energy is remarkably insensitive to mass number (between $^{52}$Ca and $^{132}$Sn) and the choice of effective interaction. For $^{90}$Zr, even the observable pygmy bump at $\sim9$ 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 [1211.0057].

## 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, $\delta\rho_p^{(1)}(r;\alpha)$ and $\delta\rho_n^{(1)}(r;\alpha)$ are everywhere out of phase ($\delta\rho_p/\delta\rho_n<0$), representing a collective core proton–neutron oscillation.
- In the skin-mode, proton and neutron transition densities are in phase in the core (interior $r$), 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,\alpha'\gamma)$ and $(\gamma,\gamma')$ 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 $\sim10$ 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 [1211.0057].

## 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 [1211.0057].

Source: https://www.emergentmind.com/topics/search-e1