---
title: Band-Projected Electric Quadrupole
url: https://www.emergentmind.com/topics/band-projected-electric-quadrupole
type: topic
---

# Band-Projected Electric Quadrupole

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 [1508.04055, 1502.04405, 1912.12487, 1812.06990].

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

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
\[
P(k)\equiv \sum_{n=1}^{N_{\rm occ}} |u_n(k)\rangle\langle u_n(k)|.
\]
This projector enters the Wilson-loop and nested-Wilson-loop constructions that define the quadrupole moment \(Q_{xy}\), and it is also the practical noninteracting reduction of the many-body quadrupole operator formalism [1912.12487, 1812.06990].

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 \(B(E2)\) expressions

Within collective Bohr-type descriptions, the electric quadrupole transition operator is written as
\[
\hat T_\mu(E2)=t\,\beta\left[D^2_{\mu0}(\Omega)\cos\gamma+\frac{1}{\sqrt2}\bigl(D^2_{\mu2}(\Omega)+D^2_{\mu,-2}(\Omega)\bigr)\sin\gamma\right].
\]
After Wigner–Eckart reduction and separation of \(\beta\)-, \(\gamma\)-, and Euler-angle integrals, the general transition probability becomes
\[
B(E2;nLK\to n'L'K')=
\frac{1}{2L+1}\,\big|\langle n'L'K'\Vert T(E2)\Vert nLK\rangle\big|^2,
\]
and, in the factorized band-projected form,
\[
B(E2)=\frac{t^2}{16\pi}\,(L\,2\,L'|K\,q\,K')^2\,I_{nL,n'L'}^2\,C_{n_\gamma K,n'_\gamma K'}^2.
\]
Here the Clebsch–Gordan coefficient carries the rotational selection rules, \(I_{nL,n'L'}\) is the radial overlap, and \(C_{n_\gamma K,n'_\gamma K'}\) is the \(\gamma\)-overlap [1108.0994].

For the ground-state band of an axially symmetric rotor, \(n=n'=0\), \(K=K'=0\), and \(L\to L-2\), the result reduces to
\[
B(E2;L\to L-2)=\frac{t^2}{16\pi}\;\frac{3\,(L-1)\,L}{(2L-1)\,(2L+1)}\;[I_{0L,0\,L-2}]^2,
\]
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 \(10\%\) up to \(L=10\), while some interband transitions are systematically underpredicted in \(\gamma\)-unstable nuclei and overpredicted in rotational nuclei [1108.0994].

The Manning–Rosen realization yields a closely related structure. In the axial prolate \((\gamma\approx 0)\) limit,
\[
B(E2;\,n,L,n_\gamma,K\to n',L',n'_\gamma,K')
=\frac{5}{16\pi}\,t^2\,
\bigl\langle L\,K;\,2\,K'-K\bigm\|L'\,K'\bigr\rangle^2
\bigl[I^{(\beta)}_{nL\to n',L'}\bigr]^2
\bigl[C_{n_\gamma,K\to n'_\gamma,K'}\bigr]^2.
\]
For intra-band \((\Delta K=0)\) transitions, \(C=\delta_{n_\gamma,n'_\gamma}\delta_{K,K'}\); for \(\Delta K=2\), the extra \(\gamma\)-overlap introduces a \(1/\sqrt{3c}\) suppression. The paper emphasizes that the main difference between the \(\gamma\)-unstable and axial-rotor limits is precisely this extra \(\Delta K=2\) coupling [1605.05270].

## 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 \(|\Phi\rangle\) obtained by a Nilsson+BCS calculation at deformation \(\epsilon_2\), with \(\epsilon_2\approx 0.48\) for \(^{40}\)Ar. Good angular momentum is restored with
\[
P^I_{MK}=\frac{2I+1}{8\pi^2}\int d\Omega\,D^I_{MK}(\Omega)\,R(\Omega),
\]
and the final yrast wave function is
\[
|\Psi^I_M\rangle=\sum_{K,\kappa} f^I_{K\kappa}\,P^I_{MK}|\Phi_\kappa\rangle,
\]
where \(|\Phi_\kappa\rangle\) spans 0-, 2-, and 4-quasiparticle configurations [1508.04055].

The laboratory-frame quadrupole operator is
\[
\hat Q_{2\mu}
= e_\pi \sum_p r_p^2 Y_{2\mu} + e_\nu \sum_n r_n^2 Y_{2\mu},
\]
with the standard effective charges \(e_\pi=1.5\,e\) and \(e_\nu=0.5\,e\). After evaluating projected matrix elements and applying the Wigner–Eckart theorem, the transition probability is
\[
B(E2;I\to I-2)=\frac{1}{2I+1}\,|\langle I-2\Vert \hat Q_2\Vert I\rangle|^2.
\]
Numerically, the reduced matrix element is obtained from the Euler-angle integral on a discrete Gauss mesh, typically \(16\times16\times16\) or finer, with rotated-overlap kernels evaluated by standard BCS-vacuum techniques and summed over the \(K\)-mixing coefficients \(f\) [1508.04055].

The intrinsic moments are parameterized by
\[
Q_0\equiv \langle\Phi|\hat Q_{20}|\Phi\rangle
\simeq \frac{3ZR_0^2}{4\pi}\beta_2\cos\gamma,\qquad
Q_2\equiv \langle\Phi|\hat Q_{22}|\Phi\rangle
\simeq \frac{3ZR_0^2}{4\pi}\beta_2\sin\gamma/\sqrt2,
\]
with \(R_0\approx 1.2\,A^{1/3}\,\mathrm{fm}\). In the PSM, \(\beta_2\) is fixed in the Nilsson+BCS basis with \(\gamma=0\), while triaxial components appear dynamically through \(K\neq 0\) admixtures in the final wave function. The active single-particle space usually includes three major shells, \(N=1,2,3\) for both protons and neutrons in the \(40\) mass region, and the BCS vacuum is generated with \(G_M\approx 18.65/A\) MeV and \(G_Q\approx 0.2\,G_M\) [1508.04055].

The PSM analysis of the superdeformed yrast band in \(^{40}\)Ar 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 \(B(E2)\) values of order \(O(100)\,\mathrm{W.u.}\) confirm strongly deformed collective rotation with \(\beta_2\approx 0.5\), while successive band crossings and pair breaking lead to a smooth reduction of \(B(E2)\) by \(\lesssim 30\%\) up to \(I\approx 18\hbar\). The deduced triaxial deformation parameters are small throughout the entire band, suggesting that triaxiality is not very important for this superdeformed band [1508.04055].

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

At leading and subleading order, the quadrupole operator has the schematic form
\[
\hat Q_{2\mu}
= e\,Q_0\,D^2_{\mu0}(\alpha,\beta,\gamma)
+ e\,\frac{Q_2}{\Lambda}\,\{\mathbf I^2,D^2_{\mu0}\}+\ldots,
\]
where \(Q_0\) is the leading quadrupole moment, \(Q_2/Q_0\sim O((\xi/\omega)^2)\), and \(\Lambda\sim\omega\) is the breakdown scale. For intra-band transitions \(J_i\to J_f=J_i-2\), the NNLO result is
\[
B(E2;J_i\to J_f)
=\frac{(eQ_0)^2}{60}\,\bigl|C^{J_f0}_{J_i\,0,\;2\,0}\bigr|^2
\Bigl[1+\frac{b}{a}\,J_i(J_i-1)\Bigr],
\]
with \(b/a\sim(\xi/\omega)^2\). In this language the reduced quadrupole moment acquires a quadratic spin dependence,
\[
Q_{if}^2 \equiv \frac{B(E2)}{|C^{J_f0}_{J_i0,20}|^2}
=Q_0^2\,[1+(b/a)\,J_i(J_i-1)].
\]
For inter-band \(\beta\to g\) transitions at leading order,
\[
B(E2;\beta,J\to g,J-2)
=\frac{C_\beta^2}{2\,C_0^2\,\omega_0}\,(eQ_0)^2\,
\bigl|C^{J-2,0}_{J\,0,\,2\,0}\bigr|^2,
\]
with analogous expressions for \(\gamma\to g\) transitions [1502.04405].

A central feature of this band-projected formulation is its explicit uncertainty model. LO truncation errors scale as \((\xi/\omega)^2 f(J)\), NLO errors as \((\xi/\omega)^4 g(J)\), and the errors grow with spin \(J\), signaling breakdown near \(J\sim \omega/\xi\). 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 \(\beta\to g\) transitions in Er and Sm that collective models overestimate by factors \(2\)–\(10\), once one measured \(\beta\to g\) matrix element is used to fix \(C_\beta\) [1502.04405].

## 5. Occupied-band quadrupole moments in crystalline insulators

In extended systems with periodic boundary conditions, a naive quadrupole \(Q_{xy}=\int \rho(r)\,x\,y\,dr\) is not gauge-invariant. The many-body formulation replaces it by an exponentiated operator. One real-space definition introduces
\[
\hat q_{xy}\equiv \frac1{L_xL_y}\sum_{R_x,R_y}R_xR_y\,\hat n_{R_x,R_y},
\]
and then
\[
Q_{xy}=\frac1{2\pi}\,\mathrm{Im}\ln\langle \Psi_G|\,e^{2\pi i \hat q_{xy}}|\Psi_G\rangle,
\]
where \(|\Psi_G\rangle\) is the many-body ground state. A related many-body multipole construction defines
\[
\hat{\mathcal U}^Q_{ab}
=\exp\!\Bigl[\frac{2\pi i}{N_aN_b}\,b_a^i\,\hat q^{\,ij}\,b_b^j\Bigr],
\]
with \(\hat q^{\,ij}=\sum_{n=1}^{N_e}\hat x_n^i\hat x_n^j\), and in \(d=2\)
\[
Q^{ab}=\frac{e}{2\pi}\,\mathrm{Im}\bigl[\ln z^{(Q)}_{xy}\bigr]\quad (\bmod\, e).
\]
These constructions require vanishing lower multipoles, specifically zero net charge and zero polarization modulo their quanta, together with integer filling per \((a,b)\) plane [1912.12487, 1812.06990].

For noninteracting insulators, the practical formulation is band-projected. With
\[
P(k)=\sum_{n=1}^{N_{\rm occ}} |u_n(k)\rangle\langle u_n(k)|,
\]
one builds the Wilson loop along \(y\) at fixed \(k_x\),
\[
W_y(k_x)\equiv P(k_x,k_{y,N_y-1})\cdots P(k_x,k_{y,0}),
\]
whose eigenvalues are written as
\[
\mathrm{Spec}[W_y(k_x)] = e^{i2\pi \nu_y^{(j)}(k_x)}.
\]
If the Wannier bands \(\nu_y^{(j)}\) split into separated sectors, one defines the sector projector
\[
P^\nu(k_x)=\sum_{j:\nu_y^{(j)}<1/2} |\nu^{(j)}(k_x)\rangle\langle \nu^{(j)}(k_x)|,
\]
and then the nested Wilson loop along \(x\),
\[
W_x^\nu(k_y)\equiv P^\nu(k_{x,N_x-1},k_y)\cdots P^\nu(k_{x,0},k_y).
\]
From this nested loop one extracts
\[
p_x^\nu = -\frac1{2\pi}\sum_{m=0}^{N_y-1}\mathrm{Arg}\bigl[\det W_x^\nu(k_{y,m})\bigr],
\]
and, in a \(C_2\)-symmetric insulator with vanishing total dipole,
\[
Q_{xy}\equiv p_x^\nu \;(\mathrm{mod}\,1)=p_y^\nu\;(\mathrm{mod}\,1).
\]
This is the canonical occupied-band version of the electric quadrupole moment [1912.12487].

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 \(W_a^{(1)}(k_b)\) to obtain a polarization band \(P^a(k_b)\), and then a second, nested loop
\[
W_{ab}^{(2)}=\prod_{k_b}\exp\!\bigl[i\,2\pi P^a(k_b)\bigr]
\]
to obtain
\[
Q^{ab}=\frac{e}{2\pi}\,\mathrm{Im}\ln W_{ab}^{(2)}.
\]
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 [1812.06990].

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

The numerical procedure for the occupied-band formulation is explicit: diagonalize \(H(k)\), form \(P(k)\), compute \(W_y(k_x)\), verify that the Wannier bands split into gapped sectors, construct \(P^\nu\), build \(W_x^\nu(k_y)\), and then evaluate \(p_x^\nu\) and \(Q_{xy}\). The nested-loop construction fails when the Wannier-sector gap closes, because the nested projector \(P^\nu\) 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 \(U_2\)-operator expectation value [1912.12487].

Symmetry controls the quantization of the band-projected quadrupole. Mirror symmetries \(M_x\) and \(M_y\) force
\[
Q_{xy}\in\{0,\tfrac12\}\quad (\mathrm{mod}\,1).
\]
Chiral symmetry does not directly quantize \(Q_{xy}\), but it pins corner states to zero energy when \(Q_{xy}=1/2\). Inversion \(C_2=M_xM_y\) enforces vanishing total polarization and is a necessary precondition for a well-defined \(Q_{xy}\). 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 \(M_x\) and \(M_y\) but preserved inversion, one still finds \(Q_{xy}=1/2\) together with sharp corner modes. The same analysis reports that \(Q_{xy}=1/2\) survives weak on-site disorder and remains quantized even for strong \(\Gamma_4\)-type disorder until the spectral gap closes [1912.12487].

A related, but technically distinct, usage of band projection appears in molecular spectroscopy. For CO\(_2\), a purely vibrational band strength is defined by summing over all allowed \(J\to J'\) branches with Boltzmann and nuclear-spin weights,
\[
S_{\mathrm{band}}(T)
= \sum_{J,J'} (2J+1)\,e^{-E_{vJ}/(kT)}\,S_Q(v,J\to v',J')/Z(T),
\]
with the sum restricted to \(J\le 40\) in the reported calculations [2106.00585]. In the Duo implementation for diatomic molecules, the band-projected line strength \(S_{fi}(J',J'')\) is defined by summing over \(M_i\), \(M_f\), and tensor components \(m\), including the nuclear-spin statistical weight \(g_{\rm ns}\) [2111.05840]. 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.

Source: https://www.emergentmind.com/topics/band-projected-electric-quadrupole