Band-Projected Electric Quadrupole

• Band-projected electric quadrupole is defined by restricting the quadrupole operator to a chosen band subspace in both nuclear and condensed-matter systems.
• It enables computation of key observables such as B(E2) in nuclear transitions and the gauge-invariant Q_xy in crystalline insulators through projection techniques.
• The approach unifies different projection procedures, facilitating effective theory constructions, numerical evaluations, and experimental data matching across various domains.

Searching arXiv for recent and directly relevant papers on band-projected electric quadrupole. Band-projected electric quadrupole denotes a quadrupole operator, quadrupole moment, or quadrupole transition observable restricted to a selected band subspace. In nuclear-structure theory, the relevant subspace is a rotational band, a projected quasiparticle configuration space, or a band-to-band sector defined by angular-momentum projection, and the main observable is usually $B(E2)$. In crystalline insulators, the relevant subspace is the occupied Bloch-band manifold, and band projection is the essential step that makes the bulk quadrupole moment $Q_{xy}$ gauge-invariant under periodic boundary conditions. The term is therefore domain-dependent, but in each usage the central operation is the same: projection of the electric quadrupole structure onto a physically distinguished sector before computing observables (Yang et al., 2015, Pérez et al., 2015, Li et al., 2019, Wheeler et al., 2018).

1. Conceptual scope and formal meaning

In nuclear applications, band projection is tied to reduced matrix elements between collective states. In the effective-theory formulation for axially symmetric systems, the projector onto band $b$ is defined as

$$P_b=\sum_{M=-I_b}^{I_b}|b,I_bM\rangle\langle b,I_bM|,$$

and the band-to-band quadrupole operator is

$$Q_{2\mu}^{(b\to b')} \equiv P_{b'}\,\hat Q_{2\mu}\,P_b,$$

with $\hat Q_{2\mu}\equiv \hat{\mathscr M}(E2)$ in Coulomb gauge. For intra-band transitions, $P_b$ enforces the band restriction and one recovers the usual $D$-function structure; for inter-band transitions, the same construction isolates the relevant collective coupling (Pérez et al., 2015).

In crystalline insulators, band projection is performed in momentum space. If $|u_n(k)\rangle$, $n=1,\dots,N_{\rm occ}$, are cell-periodic occupied-band eigenstates, the occupied-band projector is

$Q_{xy}$0

This projector enters the Wilson-loop and nested-Wilson-loop constructions that define the quadrupole moment $Q_{xy}$1, and it is also the practical noninteracting reduction of the many-body quadrupole operator formalism (Li et al., 2019, Wheeler et al., 2018).

A plausible implication is that “band-projected electric quadrupole” is best understood not as a single formalism, but as a family of projection procedures that render quadrupole observables well-defined within a chosen collective, quasiparticle, or occupied-band subspace.

2. Nuclear rotational bands and factorized $Q_{xy}$2 expressions

Within collective Bohr-type descriptions, the electric quadrupole transition operator is written as

$Q_{xy}$3

After Wigner–Eckart reduction and separation of $Q_{xy}$4-, $Q_{xy}$5-, and Euler-angle integrals, the general transition probability becomes

$Q_{xy}$6

and, in the factorized band-projected form,

$Q_{xy}$7

Here the Clebsch–Gordan coefficient carries the rotational selection rules, $Q_{xy}$8 is the radial overlap, and $Q_{xy}$9 is the $b$0-overlap (Inci et al., 2011).

For the ground-state band of an axially symmetric rotor, $b$1, $b$2, and $b$3, the result reduces to

$b$4

so all model dependence is pushed into the radial overlap. In the Morse-potential realization, ground-state-band transition rates are generally reproduced to better than $b$5 up to $b$6, while some interband transitions are systematically underpredicted in $b$7-unstable nuclei and overpredicted in rotational nuclei (Inci et al., 2011).

The Manning–Rosen realization yields a closely related structure. In the axial prolate $b$8 limit,

$b$9

For intra-band $$P_b=\sum_{M=-I_b}^{I_b}|b,I_bM\rangle\langle b,I_bM|,$$0 transitions, $$P_b=\sum_{M=-I_b}^{I_b}|b,I_bM\rangle\langle b,I_bM|,$$1; for $$P_b=\sum_{M=-I_b}^{I_b}|b,I_bM\rangle\langle b,I_bM|,$$2, the extra $$P_b=\sum_{M=-I_b}^{I_b}|b,I_bM\rangle\langle b,I_bM|,$$3-overlap introduces a $$P_b=\sum_{M=-I_b}^{I_b}|b,I_bM\rangle\langle b,I_bM|,$$4 suppression. The paper emphasizes that the main difference between the $$P_b=\sum_{M=-I_b}^{I_b}|b,I_bM\rangle\langle b,I_bM|,$$5-unstable and axial-rotor limits is precisely this extra $$P_b=\sum_{M=-I_b}^{I_b}|b,I_bM\rangle\langle b,I_bM|,$$6 coupling (Chabab et al., 2016).

3. Projected shell model formulation and superdeformed bands

In the Projected Shell Model (PSM), band-projected electric quadrupole observables are computed from an axially deformed quasiparticle vacuum $$P_b=\sum_{M=-I_b}^{I_b}|b,I_bM\rangle\langle b,I_bM|,$$7 obtained by a Nilsson+BCS calculation at deformation $$P_b=\sum_{M=-I_b}^{I_b}|b,I_bM\rangle\langle b,I_bM|,$$8, with $$P_b=\sum_{M=-I_b}^{I_b}|b,I_bM\rangle\langle b,I_bM|,$$9 for $$Q_{2\mu}^{(b\to b')} \equiv P_{b'}\,\hat Q_{2\mu}\,P_b,$$0Ar. Good angular momentum is restored with

$$Q_{2\mu}^{(b\to b')} \equiv P_{b'}\,\hat Q_{2\mu}\,P_b,$$1

and the final yrast wave function is

$$Q_{2\mu}^{(b\to b')} \equiv P_{b'}\,\hat Q_{2\mu}\,P_b,$$2

where $$Q_{2\mu}^{(b\to b')} \equiv P_{b'}\,\hat Q_{2\mu}\,P_b,$$3 spans 0-, 2-, and 4-quasiparticle configurations (Yang et al., 2015).

The laboratory-frame quadrupole operator is

$$Q_{2\mu}^{(b\to b')} \equiv P_{b'}\,\hat Q_{2\mu}\,P_b,$$4

with the standard effective charges $$Q_{2\mu}^{(b\to b')} \equiv P_{b'}\,\hat Q_{2\mu}\,P_b,$$5 and $$Q_{2\mu}^{(b\to b')} \equiv P_{b'}\,\hat Q_{2\mu}\,P_b,$$6. After evaluating projected matrix elements and applying the Wigner–Eckart theorem, the transition probability is

$$Q_{2\mu}^{(b\to b')} \equiv P_{b'}\,\hat Q_{2\mu}\,P_b,$$7

Numerically, the reduced matrix element is obtained from the Euler-angle integral on a discrete Gauss mesh, typically $$Q_{2\mu}^{(b\to b')} \equiv P_{b'}\,\hat Q_{2\mu}\,P_b,$$8 or finer, with rotated-overlap kernels evaluated by standard BCS-vacuum techniques and summed over the $$Q_{2\mu}^{(b\to b')} \equiv P_{b'}\,\hat Q_{2\mu}\,P_b,$$9-mixing coefficients $\hat Q_{2\mu}\equiv \hat{\mathscr M}(E2)$0 (Yang et al., 2015).

The intrinsic moments are parameterized by

$\hat Q_{2\mu}\equiv \hat{\mathscr M}(E2)$1

with $\hat Q_{2\mu}\equiv \hat{\mathscr M}(E2)$2. In the PSM, $\hat Q_{2\mu}\equiv \hat{\mathscr M}(E2)$3 is fixed in the Nilsson+BCS basis with $\hat Q_{2\mu}\equiv \hat{\mathscr M}(E2)$4, while triaxial components appear dynamically through $\hat Q_{2\mu}\equiv \hat{\mathscr M}(E2)$5 admixtures in the final wave function. The active single-particle space usually includes three major shells, $\hat Q_{2\mu}\equiv \hat{\mathscr M}(E2)$6 for both protons and neutrons in the $\hat Q_{2\mu}\equiv \hat{\mathscr M}(E2)$7 mass region, and the BCS vacuum is generated with $\hat Q_{2\mu}\equiv \hat{\mathscr M}(E2)$8 MeV and $\hat Q_{2\mu}\equiv \hat{\mathscr M}(E2)$9 (Yang et al., 2015).

The PSM analysis of the superdeformed yrast band in $P_b$0Ar found a perfect collective-rotor behavior for the energy levels, but the detailed wave functions at high spin are dominated by mixed 0-, 2-, and 4-quasiparticle configurations. The calculated electric quadrupole transition probabilities reproduce the known experimental data and suggest a reduced, but still significant, collectivity in the high-spin region. Large low-spin $P_b$1 values of order $P_b$2 confirm strongly deformed collective rotation with $P_b$3, while successive band crossings and pair breaking lead to a smooth reduction of $P_b$4 by $P_b$5 up to $P_b$6. The deduced triaxial deformation parameters are small throughout the entire band, suggesting that triaxiality is not very important for this superdeformed band (Yang et al., 2015).

4. Effective-theory construction of band-to-band quadrupole operators

The effective theory of axially symmetric systems reformulates electric quadrupole transitions in a model-independent way. The leading couplings to electromagnetic fields arise from minimal coupling, while subleading corrections use gauge-invariant non-minimal couplings. This construction yields transition operators that are consistent with the Hamiltonian, and its power counting provides theoretical uncertainty estimates (Pérez et al., 2015).

At leading and subleading order, the quadrupole operator has the schematic form

$P_b$7

where $P_b$8 is the leading quadrupole moment, $P_b$9, and $D$0 is the breakdown scale. For intra-band transitions $D$1, the NNLO result is

$D$2

with $D$3. In this language the reduced quadrupole moment acquires a quadratic spin dependence,

$D$4

For inter-band $D$5 transitions at leading order,

$D$6

with analogous expressions for $D$7 transitions (Pérez et al., 2015).

A central feature of this band-projected formulation is its explicit uncertainty model. LO truncation errors scale as $D$8, NLO errors as $D$9, and the errors grow with spin $|u_n(k)\rangle$0, signaling breakdown near $|u_n(k)\rangle$1. The effective theory describes ground-state-band transitions in well-deformed nuclei within LO uncertainty, captures visible quadratic deviations in transitional nuclei at NLO, and reproduces faint $|u_n(k)\rangle$2 transitions in Er and Sm that collective models overestimate by factors $|u_n(k)\rangle$3–$|u_n(k)\rangle$4, once one measured $|u_n(k)\rangle$5 matrix element is used to fix $|u_n(k)\rangle$6 (Pérez et al., 2015).

5. Occupied-band quadrupole moments in crystalline insulators

In extended systems with periodic boundary conditions, a naive quadrupole $|u_n(k)\rangle$7 is not gauge-invariant. The many-body formulation replaces it by an exponentiated operator. One real-space definition introduces

$|u_n(k)\rangle$8

and then

$|u_n(k)\rangle$9

where $n=1,\dots,N_{\rm occ}$0 is the many-body ground state. A related many-body multipole construction defines

$n=1,\dots,N_{\rm occ}$1

with $n=1,\dots,N_{\rm occ}$2, and in $n=1,\dots,N_{\rm occ}$3

$n=1,\dots,N_{\rm occ}$4

These constructions require vanishing lower multipoles, specifically zero net charge and zero polarization modulo their quanta, together with integer filling per $n=1,\dots,N_{\rm occ}$5 plane (Li et al., 2019, Wheeler et al., 2018).

For noninteracting insulators, the practical formulation is band-projected. With

$n=1,\dots,N_{\rm occ}$6

one builds the Wilson loop along $n=1,\dots,N_{\rm occ}$7 at fixed $n=1,\dots,N_{\rm occ}$8,

$n=1,\dots,N_{\rm occ}$9

whose eigenvalues are written as

$Q_{xy}$00

If the Wannier bands $Q_{xy}$01 split into separated sectors, one defines the sector projector

$Q_{xy}$02

and then the nested Wilson loop along $Q_{xy}$03,

$Q_{xy}$04

From this nested loop one extracts

$Q_{xy}$05

and, in a $Q_{xy}$06-symmetric insulator with vanishing total dipole,

$Q_{xy}$07

This is the canonical occupied-band version of the electric quadrupole moment (Li et al., 2019).

The many-body paper shows how this band-projected result emerges from the noninteracting reduction of the many-body quadrupole operator. In that derivation, one first computes a Wilson loop $Q_{xy}$08 to obtain a polarization band $Q_{xy}$09, and then a second, nested loop

$Q_{xy}$10

to obtain

$Q_{xy}$11

The formulation reproduces quantized values in topological quadrupole insulators and nonquantized continuous values in trivial insulators, provided the crystal is gapped, the many-body ground state is nondegenerate, lower multipoles vanish modulo their quanta, and Wannier functions are localized (Wheeler et al., 2018).

6. Symmetry constraints, numerical failure modes, and related usages

The numerical procedure for the occupied-band formulation is explicit: diagonalize $Q_{xy}$12, form $Q_{xy}$13, compute $Q_{xy}$14, verify that the Wannier bands split into gapped sectors, construct $Q_{xy}$15, build $Q_{xy}$16, and then evaluate $Q_{xy}$17 and $Q_{xy}$18. The nested-loop construction fails when the Wannier-sector gap closes, because the nested projector $Q_{xy}$19 becomes ill-defined. The paper therefore warns that, in indirect-gap or symmetry-broken cases, one must revert to direct real-space evaluation through the $Q_{xy}$20-operator expectation value (Li et al., 2019).

Symmetry controls the quantization of the band-projected quadrupole. Mirror symmetries $Q_{xy}$21 and $Q_{xy}$22 force

$Q_{xy}$23

Chiral symmetry does not directly quantize $Q_{xy}$24, but it pins corner states to zero energy when $Q_{xy}$25. Inversion $Q_{xy}$26 enforces vanishing total polarization and is a necessary precondition for a well-defined $Q_{xy}$27. The generalized BBH models illustrate several distinct failure or persistence modes: chiral-symmetry breaking can produce an indirect-gap phase in which corner modes are hidden in the bulk bands; mirror-symmetry breaking can make the Wannier bands gapless and render nested Wilson loops arbitrary; and, in a model with broken $Q_{xy}$28 and $Q_{xy}$29 but preserved inversion, one still finds $Q_{xy}$30 together with sharp corner modes. The same analysis reports that $Q_{xy}$31 survives weak on-site disorder and remains quantized even for strong $Q_{xy}$32-type disorder until the spectral gap closes (Li et al., 2019).

A related, but technically distinct, usage of band projection appears in molecular spectroscopy. For CO$Q_{xy}$33, a purely vibrational band strength is defined by summing over all allowed $Q_{xy}$34 branches with Boltzmann and nuclear-spin weights,

$Q_{xy}$35

with the sum restricted to $Q_{xy}$36 in the reported calculations (Yachmenev et al., 2021). In the Duo implementation for diatomic molecules, the band-projected line strength $Q_{xy}$37 is defined by summing over $Q_{xy}$38, $Q_{xy}$39, and tensor components $Q_{xy}$40, including the nuclear-spin statistical weight $Q_{xy}$41 (Somogyi et al., 2021). This usage concerns band-integrated or band-resolved intensities rather than projector-defined quadrupole operators, and it should not be conflated with the occupied-band quadrupole moment of crystalline insulators or with band-to-band quadrupole operators in nuclear collective theory.