Generalized Real-Space Quadrupole Moment

• Generalized real-space quadrupole moment is a second-rank observable capturing anisotropy in both quantum surfaces and crystalline lattice systems.
• It employs coordinate-independent constructions—from embedded surfaces and lattice operators to discrete quantum geometries—ensuring robustness under symmetry and disorder.
• The observable reveals practical insights in topological phases and loop quantum gravity by isolating shape information beyond traditional area or dipole metrics.

Searching arXiv for the cited topic and papers. arxiv_search(query="generalized real-space quadrupole moment", max_results=10) arxiv_search(query="quadrupole moment real-space operator multipole insulator", max_results=10) arxiv_search({"query":"quadrupole moment real-space operator multipole insulator", "max_results": 10}) Searching arXiv for related papers: (Goeller et al., 2018, Wheeler et al., 2018, Yang et al., 2022, Onaya et al., 21 Mar 2025, Deng et al., 1 Aug 2025, Li et al., 2019). The generalized real-space quadrupole moment denotes a family of coordinate-space constructions that extend ordinary dipole observables to second-rank moments without relying exclusively on momentum-space topology. In one line of work, it is a geometric tensor built from the unit normal field of a two-dimensional surface embedded in $\mathbb R^3$, designed to probe the anisotropy of quantum surfaces in loop quantum gravity (LQG). In another, it is a many-body or projector-based observable for crystalline and disordered lattice systems, formulated directly in real space so that quadrupole moments remain computable under periodic, open, infinite-lattice, or finite-temperature settings. Across these settings, the common purpose is to isolate shape information beyond total area or total charge, while preserving the appropriate invariances of the problem (Goeller et al., 2018, Wheeler et al., 2018, Onaya et al., 21 Mar 2025, Deng et al., 1 Aug 2025).

1. Geometric definition on embedded surfaces

For a smooth two-dimensional surface $S\subset\mathbb R^3$ with local coordinates $\sigma=(u,v)\mapsto x(u,v)$, the induced metric is

$g_{AB}(\sigma)=\partial_A x\cdot \partial_B x,$

and the oriented surface element is determined by

$$\vec N(\sigma)=\partial_u x\times \partial_v x,\qquad N(\sigma)=|\vec N(\sigma)|,\qquad \vec n(\sigma)=\vec N(\sigma)/N(\sigma).$$

The area element is $dA=N\,du\,dv$. The dual quadrupole tensor is then defined by

$$Q_{ij}=\int_S dA\, n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N(\sigma)\,n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N^{-1}(\sigma)\,N_i(\sigma)\,N_j(\sigma).$$

Its trace is

$\mathrm{Tr}\,Q=Q_{ii}=\int_S dA=\mathrm{Area}(S).$

In this formulation, the monopole is the area, the dipole is the closure vector, and the quadrupole captures the leading anisotropy of the surface, distinguishing, for example, ellipsoidal deformations from the round sphere (Goeller et al., 2018).

The construction is explicitly reparametrization-invariant. Under a change of coordinates $\sigma\mapsto \tilde\sigma(\sigma)$, the measure $dA=N\,du\,dv$ is invariant, while the unit normal $S\subset\mathbb R^3$0 behaves as a scalar under reparametrizations. Consequently, the integrand $S\subset\mathbb R^3$1 is a scalar times the invariant measure, and $S\subset\mathbb R^3$2 is coordinate-independent. This invariance is not a secondary technicality: it is the condition that makes the quadrupole a genuine observable of shape rather than of parametrization (Goeller et al., 2018).

2. Real-space quadrupole operators in lattice systems

In crystalline insulators, the direct position operator and its products are not well-defined on periodic Hilbert spaces, so real-space quadrupole observables are typically introduced in exponentiated form. A normalized lattice quadrupole operator may be written as

$S\subset\mathbb R^3$3

with corresponding many-body unitary

$S\subset\mathbb R^3$4

and quadrupole moment

$S\subset\mathbb R^3$5

A more general multipole construction on periodic lattices uses

$S\subset\mathbb R^3$6

from which the bulk quadrupole density is extracted through the phase of $S\subset\mathbb R^3$7 in the thermodynamic limit (Wheeler et al., 2018, Li et al., 2019).

A central consistency condition is that lower multipoles vanish modulo their own quanta. The many-body quadrupole operator transforms trivially under lattice translations only if the filling per $S\subset\mathbb R^3$8-plane is an integer and the dipole moment vanishes mod $S\subset\mathbb R^3$9; if the lower moments are nontrivial, the quadrupole expectation value must vanish. This constraint is fundamental to the interpretation of a bulk quadrupole as an invariant rather than an artifact of unresolved dipolar structure (Wheeler et al., 2018).

A distinct real-space route avoids momentum-space discretization altogether by using the spectral projector $\sigma=(u,v)\mapsto x(u,v)$0 and a square-root projector $\sigma=(u,v)\mapsto x(u,v)$1. On an infinite lattice, the quadrupole is represented by

$\sigma=(u,v)\mapsto x(u,v)$2

where $\sigma=(u,v)\mapsto x(u,v)$3 and $\sigma=(u,v)\mapsto x(u,v)$4. For a finite system under open boundary conditions,

$\sigma=(u,v)\mapsto x(u,v)$5

The infinite-lattice expressions assume translational invariance of $\sigma=(u,v)\mapsto x(u,v)$6 or $\sigma=(u,v)\mapsto x(u,v)$7, while the open-boundary formula is numerically useful, explicitly gauge-invariant, and emphasizes the role of edge states in polarization and quadrupole moments (Onaya et al., 21 Mar 2025).

3. Discrete and quantum-surface formulations in LQG

In LQG, a continuous surface is replaced by a polyhedral patch with $\sigma=(u,v)\mapsto x(u,v)$8 faces, each carrying a normal $\sigma=(u,v)\mapsto x(u,v)$9 and area $g_{AB}(\sigma)=\partial_A x\cdot \partial_B x,$0. The discrete quadrupole is

$g_{AB}(\sigma)=\partial_A x\cdot \partial_B x,$1

The inverse-norm weighting is the discrete analogue of the continuum factor $g_{AB}(\sigma)=\partial_A x\cdot \partial_B x,$2 and is what preserves the shape-sensitive character of the construction (Goeller et al., 2018).

Upon quantization, the face normals become SU(2) generators $g_{AB}(\sigma)=\partial_A x\cdot \partial_B x,$3 acting on the intertwiner Hilbert space

$g_{AB}(\sigma)=\partial_A x\cdot \partial_B x,$4

Because $g_{AB}(\sigma)=\partial_A x\cdot \partial_B x,$5 must be regularized, a natural choice is the shift $g_{AB}(\sigma)=\partial_A x\cdot \partial_B x,$6. Two operator orderings are then considered:

$g_{AB}(\sigma)=\partial_A x\cdot \partial_B x,$7

and

$g_{AB}(\sigma)=\partial_A x\cdot \partial_B x,$8

Both satisfy

$g_{AB}(\sigma)=\partial_A x\cdot \partial_B x,$9

so the trace reproduces the area operator. Higher invariants such as $$\vec N(\sigma)=\partial_u x\times \partial_v x,\qquad N(\sigma)=|\vec N(\sigma)|,\qquad \vec n(\sigma)=\vec N(\sigma)/N(\sigma).$$0 and $$\vec N(\sigma)=\partial_u x\times \partial_v x,\qquad N(\sigma)=|\vec N(\sigma)|,\qquad \vec n(\sigma)=\vec N(\sigma)/N(\sigma).$$1 are likewise well-defined SU(2)-invariant operators on the intertwiner space (Goeller et al., 2018).

This quantum-surface quadrupole is not redundant with area or volume observables. The standard LQG toolkit already contains area, volume, and holonomy operators, and the motivation for introducing dual multipoles is precisely that those operators do not carry information on the global shape of intertwiners. In this sense, the generalized real-space quadrupole moment functions as a shape observable for fundamental quanta of geometry (Goeller et al., 2018).

4. Relation to coherent intertwiners, Wannier topology, and bulk invariants

For Livine–Speziale coherent intertwiners,

$$\vec N(\sigma)=\partial_u x\times \partial_v x,\qquad N(\sigma)=|\vec N(\sigma)|,\qquad \vec n(\sigma)=\vec N(\sigma)/N(\sigma).$$2

the norm takes the form

$$\vec N(\sigma)=\partial_u x\times \partial_v x,\qquad N(\sigma)=|\vec N(\sigma)|,\qquad \vec n(\sigma)=\vec N(\sigma)/N(\sigma).$$3

with

$$\vec N(\sigma)=\partial_u x\times \partial_v x,\qquad N(\sigma)=|\vec N(\sigma)|,\qquad \vec n(\sigma)=\vec N(\sigma)/N(\sigma).$$4

Expanding about the stationary point $$\vec N(\sigma)=\partial_u x\times \partial_v x,\qquad N(\sigma)=|\vec N(\sigma)|,\qquad \vec n(\sigma)=\vec N(\sigma)/N(\sigma).$$5, the first derivative enforces closure,

$$\vec N(\sigma)=\partial_u x\times \partial_v x,\qquad N(\sigma)=|\vec N(\sigma)|,\qquad \vec n(\sigma)=\vec N(\sigma)/N(\sigma).$$6

and the second derivative yields the Hessian

$$\vec N(\sigma)=\partial_u x\times \partial_v x,\qquad N(\sigma)=|\vec N(\sigma)|,\qquad \vec n(\sigma)=\vec N(\sigma)/N(\sigma).$$7

where $$\vec N(\sigma)=\partial_u x\times \partial_v x,\qquad N(\sigma)=|\vec N(\sigma)|,\qquad \vec n(\sigma)=\vec N(\sigma)/N(\sigma).$$8 is the continuum-inspired quadrupole. The Gaussian approximation becomes

$$\vec N(\sigma)=\partial_u x\times \partial_v x,\qquad N(\sigma)=|\vec N(\sigma)|,\qquad \vec n(\sigma)=\vec N(\sigma)/N(\sigma).$$9

The quadrupole therefore controls the Gaussian spread of coherent intertwiners and enters directly into the semiclassical regime of spinfoam amplitudes (Goeller et al., 2018).

In higher-order topological insulators, an analogous structural role is played by nested Wilson loops and Wannier-sector polarizations. In the generalized separable BBH model,

$dA=N\,du\,dv$0

where $dA=N\,du\,dv$1 are Wannier-sector polarizations and $dA=N\,du\,dv$2 are the winding numbers of the constitutive one-dimensional chains. The same quantity appears in a Bloch representation,

$dA=N\,du\,dv$3

and its fractional part agrees with the nested-Wilson-loop quadrupole. In that model, translation symmetry, chiral symmetry, and separability are sufficient; no mirror or $dA=N\,du\,dv$4 rotation is required (Yang et al., 2022).

The many-body operator and the nested-Wilson-loop invariant are related but not interchangeable in all regimes. In generalized electric quadrupole insulators, the nested-loop construction can fail to diagnose certain real transitions, particularly when edge loops close and reopen while the bulk remains gapped or when an indirect-gap phase hides corner modes in bulk bands. By contrast, the real-space many-body formula tracks the presence of fractional corner charge and vanishes when the many-body gap closes, even indirectly (Li et al., 2019).

5. Symmetry, finite temperature, and disorder

The quantization mechanism for a real-space quadrupole moment depends on the framework. Mirror symmetries $dA=N\,du\,dv$5 or $dA=N\,du\,dv$6 enforce $dA=N\,du\,dv$7 and thus quantize $dA=N\,du\,dv$8 to $dA=N\,du\,dv$9 or $$Q_{ij}=\int_S dA\, n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N(\sigma)\,n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N^{-1}(\sigma)\,N_i(\sigma)\,N_j(\sigma).$$0 modulo $$Q_{ij}=\int_S dA\, n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N(\sigma)\,n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N^{-1}(\sigma)\,N_i(\sigma)\,N_j(\sigma).$$1. Chiral symmetry can also quantize the nested-loop result to half-integers in certain models, and generalized BBH constructions show that quadrupole moments can remain quantized even when mirror symmetries are absent. In a separable chiral-symmetric model, the quantization into $$Q_{ij}=\int_S dA\, n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N(\sigma)\,n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N^{-1}(\sigma)\,N_i(\sigma)\,N_j(\sigma).$$2 follows from the winding numbers of the one-dimensional building blocks, while in extended BBH models one numerically observes $$Q_{ij}=\int_S dA\, n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N(\sigma)\,n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N^{-1}(\sigma)\,N_i(\sigma)\,N_j(\sigma).$$3 and half-quantized edge polarizations even after breaking both $$Q_{ij}=\int_S dA\, n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N(\sigma)\,n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N^{-1}(\sigma)\,N_i(\sigma)\,N_j(\sigma).$$4 and $$Q_{ij}=\int_S dA\, n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N(\sigma)\,n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N^{-1}(\sigma)\,N_i(\sigma)\,N_j(\sigma).$$5 but preserving inversion $$Q_{ij}=\int_S dA\, n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N(\sigma)\,n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N^{-1}(\sigma)\,N_i(\sigma)\,N_j(\sigma).$$6 (Yang et al., 2022, Li et al., 2019).

A finite-temperature generalization replaces ground-state expectation values by ensemble averages. For spinless fermions on an $$Q_{ij}=\int_S dA\, n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N(\sigma)\,n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N^{-1}(\sigma)\,N_i(\sigma)\,N_j(\sigma).$$7 lattice with $$Q_{ij}=\int_S dA\, n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N(\sigma)\,n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N^{-1}(\sigma)\,N_i(\sigma)\,N_j(\sigma).$$8 orbitals,

$$Q_{ij}=\int_S dA\, n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N(\sigma)\,n_i(\sigma)\,n_j(\sigma) =\int du\,dv\,N^{-1}(\sigma)\,N_i(\sigma)\,N_j(\sigma).$$9

and with

$\mathrm{Tr}\,Q=Q_{ii}=\int_S dA=\mathrm{Area}(S).$0

the thermal density matrix

$\mathrm{Tr}\,Q=Q_{ii}=\int_S dA=\mathrm{Area}(S).$1

leads to

$\mathrm{Tr}\,Q=Q_{ii}=\int_S dA=\mathrm{Area}(S).$2

If the single-particle Hamiltonian obeys chiral symmetry $\mathrm{Tr}\,Q=Q_{ii}=\int_S dA=\mathrm{Area}(S).$3 with $\mathrm{Tr}\,Q=Q_{ii}=\int_S dA=\mathrm{Area}(S).$4, then $\mathrm{Tr}\,Q=Q_{ii}=\int_S dA=\mathrm{Area}(S).$5, every factor inside the argument is real, and

$\mathrm{Tr}\,Q=Q_{ii}=\int_S dA=\mathrm{Area}(S).$6

At $\mathrm{Tr}\,Q=Q_{ii}=\int_S dA=\mathrm{Area}(S).$7 this reproduces the zero-temperature formula in terms of occupied states, while at $\mathrm{Tr}\,Q=Q_{ii}=\int_S dA=\mathrm{Area}(S).$8 one finds $\mathrm{Tr}\,Q=Q_{ii}=\int_S dA=\mathrm{Area}(S).$9 (Deng et al., 1 Aug 2025).

Temperature and disorder alter the phase structure in ways not visible at zero temperature alone. In the isotropic BBH case $\sigma\mapsto \tilde\sigma(\sigma)$0, there is a single transition from $\sigma\mapsto \tilde\sigma(\sigma)$1 at low $\sigma\mapsto \tilde\sigma(\sigma)$2 to $\sigma\mapsto \tilde\sigma(\sigma)$3 at high $\sigma\mapsto \tilde\sigma(\sigma)$4, with critical temperature decreasing to zero as $\sigma\mapsto \tilde\sigma(\sigma)$5. In the anisotropic case $\sigma\mapsto \tilde\sigma(\sigma)$6, there can be a reentrant nontrivial window at intermediate-high temperature. With quasi-periodic disorder $\sigma\mapsto \tilde\sigma(\sigma)$7 and $\sigma\mapsto \tilde\sigma(\sigma)$8, an initially trivial system can be driven into a topological phase for moderate $\sigma\mapsto \tilde\sigma(\sigma)$9, then returned to $dA=N\,du\,dv$0 at larger $dA=N\,du\,dv$1; at intermediate temperatures the phase can exhibit multiple re-entries as $dA=N\,du\,dv$2 is swept (Deng et al., 1 Aug 2025). In the infinite-lattice projector formalism, disorder may also be included in hopping or onsite energies without destroying the validity of the finite-system formulas, provided a bulk gap remains (Onaya et al., 21 Mar 2025).

6. Physical meaning, bulk–boundary correspondence, and common pitfalls

The physical interpretation of the generalized real-space quadrupole moment is consistently shape-theoretic, but the object whose shape is being probed differs by field. In LQG, the quadrupole tensor measures the anisotropy of a quantum surface, its eigenvalues becoming quantized shape observables, and its appearance as a Hessian means it governs Gaussian shape fluctuations around a classical geometry. The proposed extension to spin-network states is to correlate quadrupole operators at different vertices so as to describe propagating shape fluctuations, viewed as discrete, gauge-invariant “shape-waves” and as a possible avenue toward modeling quantum gravitational waves in the kinematical framework (Goeller et al., 2018).

In higher-order topological matter, the quadrupole moment is tied to edge polarizations, corner charges, and adiabatic response to electric-field gradients. Changes in the phase of the many-body quadrupole operator are tied to flows of dipole current, and the operator encodes the adiabatic evolution of the system in the presence of an $dA=N\,du\,dv$3st gradient of the electric field. In the BBH family, the quadrupole pump evolves the ground-state quadrupole from $dA=N\,du\,dv$4, and the many-body operator tracks the corner charge and edge polarization throughout the cycle (Wheeler et al., 2018).

Two recurrent misconceptions are explicitly ruled out by the literature. First, fractional corner charge is not by itself sufficient to establish a nontrivial bulk quadrupole: there are mirror-symmetric models with fractional corner charge but vanishing bulk quadrupole, for which the operator $dA=N\,du\,dv$5 remains trivial (Wheeler et al., 2018). Second, symmetry requirements are model-dependent rather than universal. Early quadrupole-insulator constructions emphasized crystalline symmetries, but later results show that translation symmetry plus chiral symmetry and separability can already yield quantized $dA=N\,du\,dv$6, and finite-temperature quantization can be enforced by chiral symmetry alone (Yang et al., 2022, Deng et al., 1 Aug 2025).

Taken together, these formulations establish the generalized real-space quadrupole moment as a next-to-leading observable beyond area or dipole order. In geometry, it resolves global shape information not captured by area and volume. In condensed matter, it provides a bulk invariant and response diagnostic that survives real-space disorder, open boundaries, and thermal ensemble averaging. The common structure is the extraction of anisotropy from a coordinate-space second moment while preserving the invariances—reparametrization, lattice translation, gauge covariance, or symmetry quantization—required by the underlying theory.