401. Exact results for the Hubbard model on bipartite lattices in spatial dimensions d>1: Seven theorems from the full [SU(2)×SU(2)×U(1)]/Z2^2 symmetry

Abstract

The Hubbard model on bipartite lattices in spatial dimension d>1 remains one of the central paradigms for strongly correlated electrons, despite the scarcity of exact results beyond one dimension and special limits. On lattices such as the square, honeycomb, cubic, body-centered cubic, face-centered cubic, and diamond lattices, the model with nearest-neighbor hopping t and on-site repulsion U captures the competition between itinerancy and interaction-driven localization in a form simple enough to define precisely yet rich enough to encode magnetism, Mott physics, and pairing tendencies. A powerful route to exact statements is provided by the full global symmetry group [SU(2)×SU(2)×U(1)]/Z2^2, where the two SU(2) factors are associated with physical spin and the η-spin (pseudospin) sector, and the U(1) factor fixes the hidden charge degree of freedom. This symmetry structure allows one to derive seven exact theorems that hold on any bipartite lattice in d>1 under standard Hubbard assumptions. Together, these theorems establish rigorous constraints on eigenstates, quantum numbers, degeneracies, and selection rules, and they sharpen the physical interpretation of the model in terms of spin and η-spin algebra. The result is an exact framework that is broadly relevant to correlated materials, including cuprates, graphene, and graphene-derived systems.

1. Introduction

The Hubbard model is the canonical lattice model of correlated electrons. Its Hamiltonian reads
H = -t Σ_{⟨ij⟩,σ}(c†_{iσ}c_{jσ}+h.c.) + U Σ_i n_{i↑}n_{i↓} - μN - hS^z,
with fermionic creation and annihilation operators c†_{iσ}, c_{iσ}, number operators n_{iσ}, and total particle number N. On a bipartite lattice, the sites can be divided into two sublattices A and B such that nearest-neighbor hopping connects only opposite sublattices. This geometric property underlies the existence of the η-spin SU(2) symmetry and makes bipartite lattices qualitatively distinct from non-bipartite ones.

In dimensions d>1, exact solutions are rare. Unlike the one-dimensional Hubbard model, which is solvable by Bethe ansatz, higher-dimensional bipartite lattices generally resist exact diagonalization except for small clusters. Nonetheless, the symmetry structure of the model is exact and powerful. The full global symmetry is not merely SU(2)spin×SU(2)η×U(1), but rather the quotient [SU(2)×SU(2)×U(1)]/Z2^2, reflecting the fact that certain discrete transformations act trivially on the physical Hilbert space. This group organizes the spectrum into irreducible representations labeled by spin S, spin projection S^z, η-spin Sη, η-spin projection Sη^z, and an additional U(1) quantum number related to the number of singly occupied sites.

The seven theorems discussed here are exact consequences of this symmetry. They do not solve the model in the sense of providing all eigenvalues, but they provide rigorous structural information: the decomposition of the Hilbert space, the classification of states, the existence of highest-weight states, the degeneracy structure, and the selection rules governing operators and excitations. These results are especially valuable because they apply to physically relevant lattices in d>1, including square and honeycomb lattices relevant to layered oxides and graphene, and cubic and diamond lattices relevant to three-dimensional correlated materials.

2. The Hubbard model on bipartite lattices

Consider a bipartite lattice with Na sites. The Hamiltonian in its canonical form is
H = -t Σ_{⟨ij⟩,σ}(c†_{iσ}c_{jσ}+h.c.) + U Σ_i (n_{i↑}-1/2)(n_{i↓}-1/2).
This particle-hole symmetric form is especially convenient because it makes the spin and η-spin structures manifest. The operators
S^α = 1/2 Σ_i Σ_{σσ'} c†_{iσ} τ^α_{σσ'} c_{iσ'}
generate the physical spin SU(2), where τ^α are Pauli matrices. The η-spin generators are
η^+ = Σ_i ε_i c†_{i↑}c†_{i↓}, η^- = (η^+)†, η^z = 1/2(N - Na),
where ε_i = +1 on sublattice A and ε_i = -1 on sublattice B. The alternating phase ε_i is essential: it ensures that η^± commute with the kinetic term on a bipartite lattice.

The full symmetry is larger than the obvious spin symmetry because the Hamiltonian is invariant under independent rotations in spin and η-spin spaces, while the remaining U(1) symmetry counts the number of singly occupied sites, often associated with the c fermion number in rotated-electron formulations. The quotient by Z2^2 encodes the identification of simultaneous 2π rotations in the two SU(2) sectors with a U(1) phase.

3. Theorem 1: Exact global symmetry of the bipartite Hubbard Hamiltonian

The first theorem states that for any bipartite lattice with nearest-neighbor hopping and on-site interaction, the Hubbard Hamiltonian commutes with all generators of spin SU(2), η-spin SU(2), and the additional U(1) generator. Therefore the full global symmetry is [SU(2)×SU(2)×U(1)]/Z2^2.

This theorem is exact and independent of dimension as long as the lattice is bipartite and hopping connects only opposite sublattices. Its importance lies in the fact that it fixes the algebraic structure of the model before any approximation is made. The spin symmetry implies conservation of total spin and its projection. The η-spin symmetry implies conservation of charge-sector pseudospin and organizes states by particle-hole structure. The U(1) factor adds a third conserved label, enabling a refined classification of states beyond the usual spin and charge quantum numbers.

Physically, this theorem explains why the Hubbard model can support both magnetism and pairing-related charge fluctuations within a single exact algebraic framework. It also clarifies why the model is especially well suited to describing low-energy physics in cuprates and related materials, where spin and charge degrees of freedom are strongly intertwined.

4. Theorem 2: Decomposition of the Hilbert space into irreducible multiplets

The second theorem establishes that the full Hilbert space decomposes into irreducible representations of the symmetry group labeled by (S, S^z, Sη, Sη^z, SU(2)×SU(2)×U(1) quantum number). Every eigenstate belongs to one such multiplet, and all members of a multiplet are exactly degenerate when the Hamiltonian is SU(2)-invariant in both sectors and does not include symmetry-breaking fields.

This decomposition is not merely formal. It determines the allowed multiplicities of states at fixed particle number and fixed lattice size. In particular, the total number of many-body states can be reconstructed by summing over all allowed multiplets with their representation dimensions (2S+1)(2Sη+1) times the dimension associated with the U(1) label. The theorem thus provides a complete group-theoretical bookkeeping of the many-body spectrum.

A useful consequence is that one can classify states by highest-weight representatives only. Once a highest-weight state is known, the remaining states in the multiplet are generated by repeated application of lowering operators S^- and η^-.

5. Theorem 3: Existence and uniqueness of highest-weight states

The third theorem states that for each irreducible representation of the symmetry group, there exists a unique highest-weight state characterized by maximal spin projection S^z = S and maximal η-spin projection η^z = η. All other states in the multiplet are obtained by symmetry lowering.

This theorem is a direct consequence of the SU(2) algebra but has strong physical significance in the Hubbard context. It means that the entire spectrum can, in principle, be organized by a reduced set of canonical representatives. In particular, many-body eigenstates with complicated occupation patterns can be traced back to a highest-weight configuration with definite numbers of doubly occupied, singly occupied, and empty sites.

The uniqueness of the highest-weight state is especially helpful in constructing exact eigenstates in special sectors. For example, the fully polarized ferromagnetic state at fixed filling is a highest-weight spin state. Similarly, η-paired states at special fillings are highest-weight η-spin states. This theorem therefore connects abstract representation theory to concrete many-body wavefunctions.

6. Theorem 4: Exact degeneracy and multiplet structure under symmetry operations

The fourth theorem asserts that all states within the same irreducible multiplet are exactly degenerate under the Hubbard Hamiltonian, provided the symmetry is unbroken. In addition, the degeneracy pattern is determined entirely by the dimensions of the spin and η-spin multiplets and by the U(1) label.

This theorem implies that the energy spectrum is highly structured. If an eigenstate has quantum numbers (S, η), then there are exactly (2S+1)(2η+1) states related by the action of the SU(2) ladder operators, all with the same energy. This degeneracy survives on any bipartite lattice, including those with nontrivial coordination and geometry, such as the honeycomb and diamond lattices.

The theorem also yields selection rules for perturbations. A weak magnetic field breaks spin degeneracy but preserves η-spin symmetry; a chemical potential shift affects the η-spin sector. Thus the pattern of splittings under external fields can be predicted exactly from the symmetry content. This is a powerful diagnostic for interpreting numerical spectra and experimental response functions.

7. Theorem 5: Constraints on ground-state quantum numbers

The fifth theorem gives rigorous restrictions on the possible quantum numbers of the ground state. Depending on filling, interaction strength, and lattice bipartiteness, the ground state must lie in a symmetry sector compatible with the Lieb-type constraints on spin and η-spin. In particular, at half filling on a balanced bipartite lattice, the ground state is constrained to have η = 0, while the spin sector is determined by sublattice imbalance and interaction sign.

For repulsive U>0 on a bipartite lattice with equal sublattice sizes, the half-filled ground state typically has total spin S = 0 in the absence of frustration or explicit symmetry breaking. More generally, the theorem establishes that the ground state cannot carry arbitrary spin and charge pseudospin; its quantum numbers are restricted by the exact symmetry and by lattice topology.

This result is crucial because it provides a rigorous baseline for discussing antiferromagnetism and Mott insulating behavior. It also shows that the Hubbard model does not generically support spontaneous symmetry breaking in finite systems; rather, order emerges only in the thermodynamic limit, while the finite-size spectrum remains symmetry constrained.

8. Theorem 6: Selection rules for operators and excitations

The sixth theorem identifies exact selection rules for local and nonlocal operators. Any operator transforming as a tensor under SU(2)spin×SU(2)η×U(1) can connect only states whose quantum numbers differ according to the tensor product decomposition of the corresponding representations. As a result, spectral functions, response functions, and transition amplitudes are strongly constrained.

For example, single-electron creation and annihilation operators carry spin-1/2 and η-spin-1/2 character in a specific combined representation. Therefore, they can only connect multiplets satisfying precise changes in S, S^z, η, η^z, and the U(1) label. Similarly, spin operators connect states within the same charge sector but can change spin by one unit, while η operators connect states differing in particle number by two.

These selection rules are essential for understanding measurable quantities such as optical conductivity, neutron-scattering spectra, and photoemission. They explain why certain transitions are forbidden even when energetically allowed, and why the low-energy response of correlated systems often exhibits sharp symmetry-imposed structures.

9. Theorem 7: Exact framework for rotated electrons and hidden U(1) structure

The seventh theorem concerns the existence of a unitary transformation to rotated electrons, in which the number of singly occupied sites becomes a good quantum number associated with the hidden U(1) symmetry. This theorem provides the algebraic basis for separating charge fluctuations into doubly occupied, singly occupied, and empty-site sectors in an exact manner.

Although the explicit unitary transformation is highly nontrivial, its existence implies that the Hubbard model can be reformulated in a basis where the conserved U(1) quantum number is manifest. In that basis, the spin and η-spin degrees of freedom are carried by distinct elementary objects, allowing a clearer physical interpretation of low-energy excitations. The theorem therefore bridges the abstract symmetry group and a more intuitive quasiparticle picture.

This is especially valuable in the strong-coupling regime U/t ≫ 1, where the Hubbard model approaches effective spin models at half filling, but retains charge fluctuations away from half filling. The rotated-electron framework provides a controlled exact language for connecting these regimes without losing the underlying symmetry.

10. Physical implications for correlated materials

The seven theorems together define a rigorous skeleton for the physics of the Hubbard model in d>1. They show that much of the model’s complexity is organized by exact symmetry rather than by approximate dynamics. This has direct implications for cuprate superconductors, where square-lattice Hubbard physics is often invoked to describe antiferromagnetism, pseudogap phenomena, and d-wave pairing tendencies. It also applies to graphene and graphene-derived systems, where the honeycomb lattice and bipartite structure make η-spin considerations relevant to charge ordering and particle-hole symmetry.

In three-dimensional systems, such as cubic, body-centered cubic, face-centered cubic, and diamond lattices, the same algebraic framework applies, though the detailed ordering tendencies differ because of coordination and band structure. The exact symmetry results do not determine the phase diagram by themselves, but they provide nonnegotiable constraints that any approximate theory or numerical method must respect.

11. Conclusion

The Hubbard model on bipartite lattices in spatial dimensions d>1 is one of the most important unsolved problems in condensed matter physics. While exact solutions remain limited, the full [SU(2)×SU(2)×U(1)]/Z2^2 symmetry yields seven exact theorems that substantially clarify the structure of the model. These theorems establish the exact symmetry, irreducible multiplet decomposition, highest-weight construction, degeneracy pattern, ground-state constraints, operator selection rules, and the hidden U(1) structure associated with rotated electrons.

Taken together, they provide a robust and universal framework for analyzing the Hubbard model on square, honeycomb, cubic, body-centered cubic, face-centered cubic, diamond, and related bipartite lattices. They do not replace dynamical calculations, but they sharply delimit what those calculations can and cannot produce. In this sense, the seven theorems constitute a foundational exact result set for the study of correlated electrons in dimensions higher than one, and they offer a lasting theoretical basis for future work on magnetism, Mott transitions, pairing, and quantum materials described by the Hubbard model.