Pygmy Dipole Resonance (PDR) Overview

• Pygmy Dipole Resonance (PDR) is defined as low-energy E1 strength in neutron-rich nuclei, representing oscillations of excess neutrons against an isospin-saturated core.
• Microscopic studies using transition-density patterns and γ-decay reveal that PDR exhibits a heterogeneous, mixed isoscalar/isovector character distinct from the giant dipole resonance.
• Experimental probes and QRPA models show that PDR strength evolves with neutron number and shell structure, offering insights into neutron-skin thickness and symmetry energy.

Searching arXiv for recent and foundational PDR papers to ground the article. arxiv_search query: "pygmy dipole resonance neutron skin transition densities QRPA isoscalar isovector" The pygmy dipole resonance (PDR) is a concentration of low-lying electric dipole ($E1$) strength that appears on the low-energy tail of the isovector giant dipole resonance (IVGDR), typically around the neutron emission threshold and most prominently in nuclei with neutron excess. It is conventionally associated with oscillation of excess neutrons, or of a neutron skin, against a more isospin-saturated core, but modern microscopic analyses show that the low-energy dipole region is structurally heterogeneous: genuinely pygmy-like states coexist with more ordinary isovector dipole excitations, and the PDR itself often exhibits mixed isoscalar/isovector character rather than a purely isovector or purely isoscalar identity (Lanza et al., 2022, In et al., 10 Feb 2025).

1. Definition, scope, and phenomenological status

In current usage, the PDR denotes the low-lying $E1$ strength located below the main IVGDR and carrying only a small fraction of the total dipole strength. In neutron-rich nuclei this strength is often found around the neutron separation energy, below or around particle threshold, and is commonly discussed as a mode of the neutron-rich surface rather than of the nuclear bulk (Lanza et al., 2022). The associated strength is typically a few percent of the dipole energy-weighted sum rule (EWSR), not the dominant share characteristic of the GDR (Baran et al., 2013).

This low-energy strength is not uniform. One recurring theme in both experiment and theory is a splitting of the low-lying $E1$ distribution into a lower-energy component that responds to both isoscalar and isovector probes and a higher-energy component that remains visible mainly to electromagnetic probes and is structurally closer to the GDR tail (Lanza et al., 2014). A related systematic observation is that the PDR strength tends to grow with neutron number in isotopic chains such as Pd and Cd, while remaining modest in absolute EWSR fraction (Eriksen et al., 2014).

A macroscopic skin-based description remains useful for systematics but is not sufficient as a full microscopic definition. A recent macroscopic treatment based on a modified Isacker–Nagarajan–Warner approach and a neutron-skin prescription can reproduce the main trends of PDR energies and EWSR fractions in Ni, Sn, and Pb chains, yet explicitly concludes that this agreement does not justify identifying the PDR as a pure collective state (Plujko et al., 8 Apr 2026).

2. Microscopic signatures and isospin structure

The most widely used microscopic diagnostic of the PDR is the transition-density pattern. In the canonical picture, proton and neutron transition densities oscillate in phase in the nuclear interior, while the neutron component dominates at the surface. This distinguishes pygmy-like states from the IVGDR, where proton and neutron motion is globally out of phase (In et al., 10 Feb 2025). The review literature treats this transition-density pattern as the defining structural hallmark of the PDR rather than the low excitation energy alone (Lanza et al., 2022).

A fully self-consistent HFB+QRPA study of $^{96}$Mo with the Gogny D1M interaction makes this point explicitly. The translationally corrected operators used there are

$$\hat{O}_{E1} = \frac{eZ}{A} \sum_{n=1}^{N} r_n Y_{1m}(\hat{r}_n) - \frac{eN}{A} \sum_{p=1}^{Z} r_p Y_{1m}(\hat{r}_p),$$

$$\hat{O}_{IS} = \sum_{i=1}^{A} \left(r_i^3 - \frac{5}{3} \left<r^2\right> r_i \right) Y_{1m}(\hat{r}_i).$$

For the identified PDR-like state in $^{96}$Mo, proton and neutron transition densities are in phase in the interior, while for $r \gtrsim 7\ \mathrm{fm}$ the proton transition density vanishes and the neutron transition density remains significant; in that outer region the isoscalar and isovector transition densities have nearly equal magnitude. The same calculation also shows that the largest peak in the low-energy enhancement, at $E_x=13.5\ \mathrm{MeV}$, is instead IVGDR-like, so the low-energy enhancement as a whole cannot be equated with the PDR (In et al., 10 Feb 2025).

The mixed isospin character of the PDR has also been probed through $\gamma$-decay. In $E1$0Pb, a Skyrme particle-vibration coupling study found that $E1$1 decay from PDR states to the low-lying $E1$2 state is strongly suppressed relative to analogous decay from the IVGDR. This suppression was interpreted as evidence that the selected near-threshold PDR states are predominantly isoscalar, even though they reside in the $E1$3 channel. The same work showed that the surviving decay strength is sensitive to $E1$4-$E1$5 admixtures, placing the PDR between a simple $E1$6-$E1$7 excitation and a strongly collectivized giant mode (Lv et al., 22 Feb 2026).

3. Theoretical descriptions and the question of collectivity

Mean-field and small-amplitude approaches remain the principal microscopic framework for PDR studies. Self-consistent HF+RPA, HFB+QRPA, continuum RPA, relativistic QRPA, and their extensions all describe the PDR as a low-energy dipole mode distinct from the IVGDR, but they differ in the predicted degree of collectivity, fragmentation, and sensitivity to shell structure (Lanza et al., 2022).

A Gogny HF+RPA study of medium-heavy nuclei concluded that the PDR is less collective than the GDR. In that treatment, PDR states show protons and neutrons vibrating in phase and are dominated by particle-hole excitations involving neutrons in excess, whereas the GDR involves out-of-phase motion of all nucleons and strong sensitivity to the isovector residual interaction. The same study argued that the PDR is more strongly linked to shell structure than to fully developed collective motion (Co' et al., 2013).

A more recent QRPA analysis of spherical Mo isotopes sharpened this distinction by introducing explicit collectivity measures based on two-quasiparticle fragmentation and energy shifts. There, skin-oscillation states were found to have moderate collectivity, with substantial configuration mixing but limited coherence, while GDR states exhibited strong coherence and large energy shifts characteristic of fully developed collective motion. The same work showed that the major low-energy peak can dominate the EWSR fraction in the enhancement region without being the purest realization of a skin mode (In et al., 31 Jul 2025).

Other frameworks emphasize different aspects. The canonical-basis time-dependent HFB approach identified abrupt changes in low-energy $E1$8 strength with shell filling and showed that pairing smooths those shell-driven onsets without erasing them (Ebata et al., 2011). A symmetry-unrestricted Skyrme RPA survey for $E1$9 found strong correlations between PDR emergence, occupation of low-$E1$0 neutron orbitals, and neutron-skin development, with pronounced kinks at shell closures and a local linear correlation between the PDR energy-weighted fraction and $E1$1 within isotopic chains (Inakura et al., 2011). By contrast, Vlasov dynamics, which omits shell effects, still produces a low-energy dipole mode with the empirical characteristics of the PDR and identifies it with a low-lying isoscalar toroidal mode, suggesting that some PDR features are generic collective consequences of neutron excess rather than exclusively shell-driven (Urban, 2011).

4. Neutron skin, shell structure, and symmetry energy

The relation between the PDR and neutron-skin thickness is one of the most persistent themes in the literature. A transport-model study based on the Landau–Vlasov equation found that the PDR centroid is well described by

$E1$2

while the strength exhausted in the pygmy region is correlated with neutron-skin thickness. In that analysis, an increase of $E1$3 in the pygmy EWSR corresponds to a $E1$4 increase of the neutron skin, and the centroid is far less sensitive to the symmetry-energy parametrization than the strength itself (Baran et al., 2013).

A broad Skyrme RPA survey reached a related but more local conclusion. Within isotopic chains, especially in the shell-filling intervals $E1$5 and $E1$6, the fraction $E1$7 and the neutron-skin thickness $E1$8 show an approximately linear correlation with slope $E1$9. The same survey tied the onset of stronger PDR strength to occupation of weakly bound low-$^{96}$0 orbitals, with notable increases at $^{96}$1, $^{96}$2, and $^{96}$3 (Inakura et al., 2011).

A complementary interpretation emphasizes the density dependence of the symmetry energy. In a schematic quantum many-body treatment with extended separable interactions, the PDR in $^{96}$4Sn emerges around $^{96}$5–$^{96}$6 MeV and exhausts about $^{96}$7 of the EWSR. The proposed mechanism is a combined dynamics of the neutron skin and of core isovector polarization: the residual dipole interaction is weaker in the low-density surface than in the bulk, producing a low-energy branch with pronounced isoscalar features in the transition densities (Baran et al., 2020).

Recent Mo-isotope studies generalize the skin argument beyond neutron skins alone. Along the spherical $^{96}$8Mo chain, enhanced low-energy dipole strength correlates with development of either proton or neutron skins, with $^{96}$9Mo near the crossover and the clearest systematic growth on the neutron-rich side. This suggests that surface asymmetry itself is a key control parameter, while the detailed microscopic realization of the enhancement remains isotope dependent (In et al., 31 Jul 2025).

5. Experimental probes and structural selectivity

The PDR is unusual in that no single probe exhausts its structural information. Photon scattering, relativistic Coulomb excitation, forward-angle proton scattering, $$\hat{O}_{E1} = \frac{eZ}{A} \sum_{n=1}^{N} r_n Y_{1m}(\hat{r}_n) - \frac{eN}{A} \sum_{p=1}^{Z} r_p Y_{1m}(\hat{r}_p),$$0 scattering, heavy-ion inelastic scattering, and transfer reactions each emphasize different components of the low-lying dipole response (Lanza et al., 2022).

Forward-angle $$\hat{O}_{E1} = \frac{eZ}{A} \sum_{n=1}^{N} r_n Y_{1m}(\hat{r}_n) - \frac{eN}{A} \sum_{p=1}^{Z} r_p Y_{1m}(\hat{r}_p),$$1 MeV proton scattering on $$\hat{O}_{E1} = \frac{eZ}{A} \sum_{n=1}^{N} r_n Y_{1m}(\hat{r}_n) - \frac{eN}{A} \sum_{p=1}^{Z} r_p Y_{1m}(\hat{r}_p),$$2Pb established a benchmark case. In the angular interval $$\hat{O}_{E1} = \frac{eZ}{A} \sum_{n=1}^{N} r_n Y_{1m}(\hat{r}_n) - \frac{eN}{A} \sum_{p=1}^{Z} r_p Y_{1m}(\hat{r}_p),$$3–$$\hat{O}_{E1} = \frac{eZ}{A} \sum_{n=1}^{N} r_n Y_{1m}(\hat{r}_n) - \frac{eN}{A} \sum_{p=1}^{Z} r_p Y_{1m}(\hat{r}_p),$$4, the $$\hat{O}_{E1} = \frac{eZ}{A} \sum_{n=1}^{N} r_n Y_{1m}(\hat{r}_n) - \frac{eN}{A} \sum_{p=1}^{Z} r_p Y_{1m}(\hat{r}_p),$$5 cross section is essentially Coulomb excitation, allowing direct extraction of $$\hat{O}_{E1} = \frac{eZ}{A} \sum_{n=1}^{N} r_n Y_{1m}(\hat{r}_n) - \frac{eN}{A} \sum_{p=1}^{Z} r_p Y_{1m}(\hat{r}_p),$$6 strengths. A multipole decomposition analysis based on QPM angular distributions identified PDR-like states below $$\hat{O}_{E1} = \frac{eZ}{A} \sum_{n=1}^{N} r_n Y_{1m}(\hat{r}_n) - \frac{eN}{A} \sum_{p=1}^{Z} r_p Y_{1m}(\hat{r}_p),$$7 and yielded for the PDR in $$\hat{O}_{E1} = \frac{eZ}{A} \sum_{n=1}^{N} r_n Y_{1m}(\hat{r}_n) - \frac{eN}{A} \sum_{p=1}^{Z} r_p Y_{1m}(\hat{r}_p),$$8Pb a summed strength

$$\hat{O}_{E1} = \frac{eZ}{A} \sum_{n=1}^{N} r_n Y_{1m}(\hat{r}_n) - \frac{eN}{A} \sum_{p=1}^{Z} r_p Y_{1m}(\hat{r}_p),$$9

and centroid

$$\hat{O}_{IS} = \sum_{i=1}^{A} \left(r_i^3 - \frac{5}{3} \left<r^2\right> r_i \right) Y_{1m}(\hat{r}_i).$$0

The same study showed that Coulomb–nuclear interference in the angular distributions is sensitive to whether a given $$\hat{O}_{IS} = \sum_{i=1}^{A} \left(r_i^3 - \frac{5}{3} \left<r^2\right> r_i \right) Y_{1m}(\hat{r}_i).$$1 state has PDR-like or GDR-like transition densities (Poltoratska et al., 2012).

Hadronic isoscalar probes reveal the complementary side of the mode. In $$\hat{O}_{IS} = \sum_{i=1}^{A} \left(r_i^3 - \frac{5}{3} \left<r^2\right> r_i \right) Y_{1m}(\hat{r}_i).$$2Sn, semiclassical $$\hat{O}_{IS} = \sum_{i=1}^{A} \left(r_i^3 - \frac{5}{3} \left<r^2\right> r_i \right) Y_{1m}(\hat{r}_i).$$3 calculations using RQTBA transition densities reproduced the experimentally observed splitting of the low-lying $$\hat{O}_{IS} = \sum_{i=1}^{A} \left(r_i^3 - \frac{5}{3} \left<r^2\right> r_i \right) Y_{1m}(\hat{r}_i).$$4 strength: strong $$\hat{O}_{IS} = \sum_{i=1}^{A} \left(r_i^3 - \frac{5}{3} \left<r^2\right> r_i \right) Y_{1m}(\hat{r}_i).$$5-excitation for lower-energy states, followed by a sharp suppression above about $$\hat{O}_{IS} = \sum_{i=1}^{A} \left(r_i^3 - \frac{5}{3} \left<r^2\right> r_i \right) Y_{1m}(\hat{r}_i).$$6 MeV. The integrated cross section up to $$\hat{O}_{IS} = \sum_{i=1}^{A} \left(r_i^3 - \frac{5}{3} \left<r^2\right> r_i \right) Y_{1m}(\hat{r}_i).$$7 MeV was calculated as $$\hat{O}_{IS} = \sum_{i=1}^{A} \left(r_i^3 - \frac{5}{3} \left<r^2\right> r_i \right) Y_{1m}(\hat{r}_i).$$8, compared with the experimental $$\hat{O}_{IS} = \sum_{i=1}^{A} \left(r_i^3 - \frac{5}{3} \left<r^2\right> r_i \right) Y_{1m}(\hat{r}_i).$$9, showing that the lower-energy group carries a stronger isoscalar interior component while the higher group is structurally closer to the GDR tail (Lanza et al., 2014).

The Oslo-method $^{96}$0-strength-function studies add a different perspective. In $^{96}$1Pd, all four isotopes show a PDR-like enhancement above $^{96}$2 MeV, modeled with a Gaussian centered at

$^{96}$3

while only $^{96}$4Pd shows a low-energy enhancement below $^{96}$5 MeV. Under the assumption that the excess strength is entirely $^{96}$6, the integrated PDR exhausts about $^{96}$7 to $^{96}$8 of the Thomas–Reiche–Kuhn sum rule and increases with neutron number (Eriksen et al., 2014).

Transfer reactions probe another layer of structure. In $^{96}$9Pb, the $r \gtrsim 7\ \mathrm{fm}$0 reaction, combined with $r \gtrsim 7\ \mathrm{fm}$1, showed that many low-lying $r \gtrsim 7\ \mathrm{fm}$2 states carry substantial neutron single-particle strength, especially the 5292-keV and 5947-keV states, but that this single-particle character does not map one-to-one onto $r \gtrsim 7\ \mathrm{fm}$3. The low-lying states can therefore have dominant neutron $r \gtrsim 7\ \mathrm{fm}$4-$r \gtrsim 7\ \mathrm{fm}$5 building blocks and still display transition densities associated with neutron-skin oscillation (Spieker et al., 2020).

6. Extensions beyond the standard ground-state $r \gtrsim 7\ \mathrm{fm}$6 picture and unresolved issues

One active extension concerns PDR-like modes built on excited states. In $r \gtrsim 7\ \mathrm{fm}$7Ge, a phonon-damping-model study with exact pairing predicted two enhanced $r \gtrsim 7\ \mathrm{fm}$8-transitions at $r \gtrsim 7\ \mathrm{fm}$9 and $E_x=13.5\ \mathrm{MeV}$0 MeV at $E_x=13.5\ \mathrm{MeV}$1 MeV, interpreted as a hot PDR built mainly on the first $E_x=13.5\ \mathrm{MeV}$2 state. The same work argued that the appearance of this hot PDR in the narrow interval $E_x=13.5\ \mathrm{MeV}$3–$E_x=13.5\ \mathrm{MeV}$4 MeV implies a violation of the Brink–Axel hypothesis in the PDR region (Phuc et al., 8 Apr 2025).

That theme was strengthened experimentally by the first evidence for $E_x=13.5\ \mathrm{MeV}$5 components of the PDR. In the $E_x=13.5\ \mathrm{MeV}$6-delayed $E_x=13.5\ \mathrm{MeV}$7 decay of the $E_x=13.5\ \mathrm{MeV}$8 isomer $E_x=13.5\ \mathrm{MeV}$9Ga, subthreshold resonant structures in $\gamma$0Ge at $\gamma$1 and $\gamma$2 keV were identified as $\gamma$3 components of a PDR built on the low-lying $\gamma$4 state. This extends the PDR concept beyond the standard $\gamma$5 ground-state picture and implies that, in thermally excited astrophysical environments, low-energy dipole strength can reside in multiplets built on excited configurations (Li et al., 31 Oct 2025).

Despite the breadth of existing work, several issues remain open. The degree of collectivity is still debated; low-energy enhancement near threshold is not automatically pygmy in nature; deformation effects remain incompletely resolved; and standard QRPA-based approaches do not include the full fragmentation induced by multiphonon admixtures and particle-vibration coupling. Recent Mo studies state this explicitly: the low-energy enhancement region contains both skin-oscillation states and stronger peaks with more intricate structure, and further investigation is needed to distinguish genuine pygmy modes from ordinary GDR-tail excitations on a state-by-state basis (In et al., 31 Jul 2025). A plausible implication is that the PDR is best regarded not as a single universal mode, but as a class of low-energy dipole excitations generated by neutron-rich surface dynamics, mixed isospin character, and incomplete coherence among many microscopic configurations.