Angular Dipole Deformation in Nuclei & Atoms

• Angular dipole deformation is the anisotropic response of dipole modes due to broken spherical symmetry, resulting in K-resolved splitting.
• It leads to strength fragmentation and centroid shifts in L=1 transitions, affecting both Gamow-Teller and isovector electric dipole responses.
• Advanced microscopic frameworks like the axially deformed pnRQRPA capture these effects, yielding improved agreement with experimental observations and astrophysical predictions.

Angular dipole deformation denotes deformation- or orientation-induced anisotropy of dipole response or of dipole-like collective variables. In the nuclear setting treated most explicitly in the axially deformed proton-neutron relativistic quasiparticle random-phase approximation, it is the response of $L=1$ spin-isospin modes to axial quadrupole deformation: the spherical $2J+1$ degeneracy is lifted, only the projection $K$ of angular momentum on the symmetry axis remains a good quantum number, and the $J=1$ sector splits into $K=0$ and $K=\pm1$, with consequent centroid shifts, $K$-dependent collectivity, and fragmentation of the strength function (Ravlić et al., 2024). Closely related usages appear in isovector electric dipole response, magnetic-rotational dipole bands, anisotropic dipole-dipole interactions, beyond-dipole high-harmonic generation, and gauge-invariant $l=1$ perturbation theory, where reduced symmetry or directional structure resolves a formerly isotropic dipole degree of freedom into distinct angular components (Wang et al., 2016, Trivedi et al., 2012, Ravets et al., 2015, Albar et al., 19 Mar 2025, Nakamura, 2021).

1. Symmetry reduction and the $K$-resolved dipole

The central mechanism is symmetry breaking. In a spherical nucleus, dipole excitations are classified by total angular momentum $J$ and magnetic quantum number $2J+1$0, and the familiar $2J+1$1 degeneracy holds. In an axially deformed nucleus, spherical symmetry is broken and only the projection $2J+1$2 on the symmetry axis remains conserved. For charge-exchange modes this yields the selection rules

$2J+1$3

so the response separates into $2J+1$4 blocks rather than $2J+1$5 blocks. For $2J+1$6 modes, the spherical triplet becomes $2J+1$7 and $2J+1$8, with the $2J+1$9 branches degenerate under time reversal. In the deformed Gamow-Teller case,

$K$0

which makes the $K$1-decomposition explicit (Ravlić et al., 2024).

The same kinematic idea governs the isovector electric dipole response in an axially symmetric weakly bound nucleus. There the good quantum number is again $K$2, now interpreted as the projection of the dipole multipolarity $K$3. The $K$4 branch is longitudinal, polarized along the symmetry axis, whereas $K$5 denotes the two degenerate transverse polarizations. Angular dipole deformation therefore appears as different centroid energies, line shapes, and transition-current topologies for longitudinal and transverse dipole motion (Wang et al., 2016).

This suggests a broad but coherent classification. In deformed nuclear spectroscopy, angular dipole deformation usually means $K$6-splitting and redistribution of strength. In magnetic rotation it refers to the angular orientation of a magnetic dipole vector relative to the total spin. In AMO settings it often denotes anisotropic dipole-dipole coupling or the deflection of an induced dipole away from a principal polarization axis. In gravitational perturbation theory the same language is extended to $K$7 axial or mass-dipole sectors, where angular structure encodes spin or center-of-mass dipole charges (Trivedi et al., 2012, Ravets et al., 2015, Albar et al., 19 Mar 2025, Nakamura, 2021).

2. Axially deformed pnRQRPA and the microscopic treatment of $K$8 modes

A fully covariant microscopic realization is provided by the axially deformed relativistic Hartree-Bogoliubov plus proton-neutron relativistic quasiparticle random-phase approximation framework. The ground state is obtained with the axially deformed relativistic Hartree-Bogoliubov model using relativistic density-dependent point-coupling energy density functionals, specifically DD-PC1 and DD-PCX. The excitations are then constructed in an axially deformed pnRQRPA formulated in linear response, so that axial symmetry reduces the problem to separate $K$9 blocks and enables explicit strength decomposition by $J=1$0 (Ravlić et al., 2024).

The external charge-exchange operators are the Fermi operator,

$J=1$1

the Gamow-Teller operator,

$J=1$2

and the spin-dipole family,

$J=1$3

The associated strength function is

$J=1$4

and in axial symmetry the operator and response are resolved into $J=1$5-components evaluated with axially deformed Dirac spinors (Ravlić et al., 2024).

Self-consistency is enforced at the level of both particle-hole and particle-particle channels. The particle-hole interaction contains the isovector-vector term inherited from the ground-state functional and an isovector-pseudovector pion-like term added only in the residual interaction; its Landau-Migdal parameter is adjusted to reproduce Gamow-Teller centroids, with $J=1$6 for DD-PC1 and $J=1$7 for DD-PCX. The pairing sector uses a separable finite-range interaction in both isovector $J=1$8 and isoscalar $J=1$9 channels, and the application to unnatural-parity transitions sets the isoscalar pairing scale to $K=0$0 (Ravlić et al., 2024).

Computationally, the formal pnRQRPA amplitudes satisfy the standard QRPA matrix equation in each $K=0$1 block, but the practical solver avoids explicit diagonalization of the full two-quasiparticle space. Instead it solves the Bethe-Salpeter equation for the response matrix in a reduced space of separable interaction channels. This substantially lowers the cost, and singular value decomposition further reduces the interaction-channel dimension without loss of accuracy. The axial harmonic-oscillator basis converges around $K=0$2 for $K=0$3 nuclei, while production calculations use $K=0$4; strength distributions are smeared with $K=0$5 MeV, and anti-particle contributions in the charge-exchange channel are neglected as negligible (Ravlić et al., 2024).

Benchmarking in the spherical limit is an essential part of the construction. For doubly magic $K=0$6O, the deformed solver constrained to spherical shape reproduces spherical pnRQRPA results for the isobaric analog state and GT$K=0$7 response, including the expected $K=0$8 and $K=0$9 degeneracy and the factor-of-three relation between total and $K=\pm1$0-resolved Gamow-Teller strength (Ravlić et al., 2024).

3. Spin-isospin angular dipole deformation in deformed nuclei

The resulting phenomenology is sharply mode dependent. Fermi transitions are almost shape independent because only $K=\pm1$1 contributes and the operator carries isospin but no spin or coordinate dependence. In $K=\pm1$2Fe the isobaric analog state energy differs by at most $K=\pm1$3 MeV between oblate, spherical, and prolate shapes. For $K=\pm1$4Fe, the experimental IAS centroid is around $K=\pm1$5 MeV, about $K=\pm1$6 MeV above calculated values, and DD-PCX shifts the theory upward by about $K=\pm1$7 MeV (Ravlić et al., 2024).

Gamow-Teller transitions are markedly more sensitive. In $K=\pm1$8Ni GT$K=\pm1$9, the spherical strength is a single peak at $K$0 MeV. Constraining the nucleus to the oblate minimum at $K$1 produces broad fragmentation, whereas the prolate shape at $K$2 generates two peaks near $K$3 MeV and $K$4 MeV; with $K$5 MeV these merge into a broadened resonance, and with $K$6 MeV they separate clearly. Across the deformed ground states of $K$7Ni and $K$8Zn, which are oblate, and $K$9Fe and $l=1$0Ti, which are prolate, the deformed pnRQRPA strength functions display fragmented main peaks of reduced height and agree better with measured GT$l=1$1 distributions than spherical calculations. DD-PCX generally shifts centroids a few hundred keV higher than DD-PC1 (Ravlić et al., 2024).

The $l=1$2-resolved GT$l=1$3 pattern in the Fe chain shows a robust correlation with the sign of $l=1$4. For prolate shapes, low-lying GT$l=1$5 strength is dominated by $l=1$6, whereas the giant GT resonance region is dominated by $l=1$7. For oblate shapes the trend is reversed: $l=1$8 is pushed to lower energies and $l=1$9 to higher energies, so the high-energy giant-resonance region carries more $K$0 content. The magnitude of the $K$1-splitting grows with $K$2. The mechanism is explicitly traced to Nilsson-level splitting and mixing, the increased two-quasiparticle level density in a deformed mean field, and pairing-induced changes in the coherence of low-lying configurations; the TV and TPV residual interactions set collectivity and centroid energies, while isoscalar pairing modifies low-lying GT strength in the unnatural-parity channel (Ravlić et al., 2024).

Spin-dipole transitions exhibit a related but more structured behavior because they combine $K$3, $K$4, and explicit radial dependence. In $K$5Fe, the centroid energies $K$6 satisfy the ordering

$K$7

with the oblate values $K$8, $K$9, and $J$0 MeV, respectively. This ordering persists across shapes and isotopes. With increasing neutron number from $J$1Fe to $J$2Fe to $J$3Fe, all three spin-dipole centroids move to lower energy and the total strengths increase. The $J$4-decomposition is again shape dependent: $J$5 carries only $J$6; $J$7 contains $J$8 and $J$9, with the relative ordering of their centroids reversing between prolate and oblate $2J+1$00Fe; $2J+1$01 contains $2J+1$02, with $2J+1$03 dominating near $2J+1$04 MeV in prolate $2J+1$05Fe but controlling low-lying strength near $2J+1$06 MeV in the oblate shape. Within the deformed pnRQRPA, the non-energy-weighted SD sum rule $2J+1$07 is preserved, consistent with the self-consistent residual interaction (Ravlić et al., 2024).

These redistributions have direct astrophysical implications. Deformation-induced fragmentation modifies the low-energy GT$2J+1$08 and SD$2J+1$09 strength available within the $2J+1$10 window and therefore can change beta-decay half-lives. The improved match to GT$2J+1$11 strength in deformed $2J+1$12-shell nuclei indicates that deformation is needed for reliable electron-capture rates. The deformation dependence of SD modes also affects first-forbidden contributions relevant to neutrino-nucleus reactions and hot astrophysical environments (Ravlić et al., 2024).

4. Electric dipole anisotropy, giant resonances, and shape diagnostics

The same angular logic governs the isovector electric dipole response, but here weak binding and surface diffuseness can substantially alter the deformation pattern. In $2J+1$13Mg, a fully self-consistent continuum finite-amplitude QRPA on a large deformed coordinate-space mesh was used to study both prolate and oblate shapes. The prolate ground state has $2J+1$14, while an oblate local minimum at $2J+1$15 lies about $2J+1$16 MeV higher. The $2J+1$17-resolved strengths show a conventional giant dipole splitting of about $2J+1$18 MeV for both shapes, consistent with the near-linear hydrodynamic expectation. The pygmy dipole resonance behaves differently: in the prolate case the lowest $2J+1$19 branch appears at $2J+1$20 MeV and the $2J+1$21 branch is centered around $2J+1$22 MeV, giving $2J+1$23 MeV, larger than the simple hydrodynamic estimate of about $2J+1$24 MeV; in the oblate case the $2J+1$25 component dominates the low-energy strength and the splitting is compressed to about $2J+1$26 MeV, much smaller than the estimate of about $2J+1$27 MeV. Transition-current analysis identifies the lowest prolate $2J+1$28 pygmy mode as a compressional core-halo oscillation with a single nodal boundary around $2J+1$29 fm, whereas higher-energy PDR and GDR branches display progressively more complicated nodal and circulation-like patterns. Accurate angular decomposition requires a very large box, with $2J+1$30 fm found necessary; smaller boxes generate false peaks and distort the surface flow (Wang et al., 2016).

A common misconception is that deformation should always appear as a clean double-hump IVGDR in well-deformed nuclei. Proton inelastic scattering at $2J+1$31 MeV and $2J+1$32 on $2J+1$33Nd and $2J+1$34Sm shows a more nuanced behavior. Along the chain the adopted axial quadrupole deformation increases from $2J+1$35 in $2J+1$36Nd to $2J+1$37 in $2J+1$38Sm. In the well-deformed nuclei $2J+1$39Nd and $2J+1$40Sm, the extracted photoabsorption cross sections show a pronounced asymmetry rather than the distinct double-hump structure expected from simple $2J+1$41-splitting. The interpretation advanced in the self-consistent SLy6 Skyrme separable RPA analysis is that these nuclei lie near the critical point of the phase transition from vibrators to rotors and possess a soft quadrupole deformation potential; shape fluctuations smear the $2J+1$42 and $2J+1$43 branches and reduce the distinctness of the low-energy component (Donaldson et al., 2016).

Polarized $2J+1$44-decay provides an even more differential probe of angular dipole deformation. In $2J+1$45Sm, linearly polarized quasimonochromatic photon beams were used to disentangle elastic scattering and $2J+1$46-Smekal-Raman scattering from the isovector giant dipole resonance through their distinct angular distributions,

$2J+1$47

$2J+1$48

A simultaneous fit to $2J+1$49-decay branching ratios and photoabsorption data yielded three IVGDR components and extracted the ground-state deformation parameters

$2J+1$50

Photoabsorption alone gave a flat $2J+1$51 posterior with a $2J+1$52 credible upper limit $2J+1$53, whereas inclusion of the $2J+1$54-decay observables reduced that upper limit to $2J+1$55. Here angular dipole deformation is not only the splitting of dipole eigenmodes along principal axes, but also the anisotropic emission pattern of their decay amplitudes (Kleemann et al., 2024).

At finite temperature and angular momentum, deformation can also be inferred indirectly from the giant dipole width. Pandit and collaborators established a non-linear, mass-dependent relation for $2J+1$56,

$2J+1$57

with

$2J+1$58

$2J+1$59

The subtraction of $2J+1$60 encodes the GDR-induced quadrupole moment and explains why the width-deformation relation is not simply linear. The correlation holds for $2J+1$61 MeV and below the approximate Jacobi-transition scale $2J+1$62 (Pandit et al., 2013).

5. Intrinsic dipoles, rotation, and parity mixing

In magnetic-rotational spectroscopy, angular dipole deformation refers less to collective $2J+1$63-splitting than to the directional geometry of a magnetic dipole built from high-$2J+1$64 quasiparticles. In $2J+1$65In, three dipole bands were observed up to $2J+1$66 MeV and $2J+1$67. The positive-parity band A shows a monotonic decrease of $2J+1$68 with spin, from $2J+1$69 at $2J+1$70 to $2J+1$71 at $2J+1$72, $2J+1$73 at $2J+1$74, and $2J+1$75 at $2J+1$76. Tilted-axis cranking assigns this band the configuration $2J+1$77, with a favored static solution near $2J+1$78, $2J+1$79, corresponding to $2J+1$80, and an almost constant tilt angle of the rotational axis near $2J+1$81. Band C is assigned $2J+1$82, with $2J+1$83 and $2J+1$84. The shears interpretation uses

$2J+1$85

so the decrease of $2J+1$86 tracks the closing of the proton and neutron blades and the concomitant reduction of the transverse magnetic moment. The inferred intrinsic quadrupole moments remain small, $2J+1$87 eb for band A and $2J+1$88 eb for band C, which is why the angular, rather than strongly collective quadrupole, character dominates the dipole bands (Trivedi et al., 2012).

A different but related body-fixed dipole appears in octupole-deformed nuclei. For reflection-asymmetric, axially symmetric shapes, the surface is expanded as

$2J+1$89

and a nonzero octupole parameter $2J+1$90 induces a collective intrinsic electric dipole moment aligned with the symmetry axis. To leading order in $2J+1$91, the intrinsic electric dipole and Schiff moments are

$2J+1$92

$2J+1$93

In a laboratory eigenstate of good parity and angular momentum, rotational averaging makes these expectation values vanish. They reappear only when $2J+1$94-violating interactions mix opposite-parity rotational doublets, with the laboratory moments

$2J+1$95

For $2J+1$96Ra the tabulated values are $2J+1$97, $2J+1$98, $2J+1$99, $K$00, and $K$01. The quoted uncertainty of the estimates may exceed a factor of $K$02, but the formal point is clear: angular dipole deformation here is the alignment problem of a collective intrinsic dipole with respect to the spin axis under parity mixing (Flambaum et al., 2024).

6. Atomic, optical, and relativistic extensions

In Rydberg physics, angular dipole deformation becomes the anisotropic reshaping of the interaction graph by the angle between the internuclear axis and the quantization axis. At a Stark-tuned Förster resonance between two individual $K$03Rb atoms, the isolated exchange channel obeys

$K$04

with a magic angle

$K$05

where the resonant coupling vanishes. For $K$06, the measured oscillation frequency at $K$07 is $K$08 MHz, minimal oscillations are seen near $K$09, and the $K$10 frequency is reduced consistently with the $K$11 form. In many-atom simulations of resonant $K$12 transport, the same anisotropy enters through

$K$13

leading to faster propagation for transverse fields than for axial fields and to a calculated localization length of $K$14 in one representative geometry (Ravets et al., 2015, Bigelow et al., 2015).

In strong-field and structured-light physics, the phrase denotes a dynamical rotation of the driven electronic dipole away from the principal polarization axis. Beyond the electric-dipole approximation, a linearly polarized OAM Bessel beam is treated with the full minimal-coupling Hamiltonian,

$K$15

so that spatial gradients and magnetic fields induce nonzero transverse dipole components. A practical deformation angle is defined by

$K$16

Using the reported off-axis dipole amplitudes at the hot spot, one infers $K$17 rad $K$18. Despite this small time-domain deflection, the spectral consequences are large: the second-harmonic emission along $K$19 is approximately doubled relative to a matched plane wave, and the beyond-dipole correction reaches about $K$20 in amplitude, while the third-harmonic correction remains about $K$21. Angular dipole deformation is therefore observed through symmetry breaking, transverse emission, and the appearance of even harmonics (Albar et al., 19 Mar 2025).

In fast charged-particle scattering from atoms, the angular distribution of secondary electrons relative to the momentum-transfer direction $K$22 departs from the familiar photoionization dipole pattern because the transition operator contains a full multipole expansion,

$K$23

The differential generalized oscillator strength is expanded in Legendre polynomials,

$K$24

For $K$25-subshells in the optical limit, the scattering pattern is $K$26, not the photoionization form $K$27, and the non-dipole terms scale with $K$28 rather than $K$29. Quadrupole contributions are therefore much less suppressed than in photoionization. Calculations for He, Ar, and Xe, including the collectivized Xe $K$30 subshell, show prominent deviations even at very small $K$31 (Amusia et al., 2010).

Relativistic perturbation theory extends the terminology further. In Nakamura’s gauge-invariant any-order treatment of perturbations on Schwarzschild spacetime, angular dipole deformation is the $K$32 odd-parity sector. In vacuum and the stationary limit it yields

$K$33

equivalently $K$34 with $K$35, reproducing the slow-rotation limit of Kerr. In the eikonal framework for classical gravitational scattering, the same language is applied to changes in the system’s spatial angular momentum and center-of-mass dipole charges. The ambiguity-free center-of-mass dipole is

$K$36

which removes recoil drift and Coulombic ambiguities from the boost component $K$37. Here angular dipole deformation no longer refers to a material dipole oscillation, but to the $K$38 charge content of the gravitational field and its radiative loss (Nakamura, 2021, Heissenberg et al., 2024).

Across these domains, the unifying structure is the same. A system that is isotropic or effectively dipolar in a high-symmetry limit acquires direction-dependent branches, couplings, or observables once symmetry is reduced by deformation, orientation, gradients, or rotation. In nuclear structure this appears most canonically as $K$39-splitting and strength fragmentation of $K$40 modes; in neighboring fields it persists as anisotropic polarization, angularly selective coupling, or $K$41 charge transport.