XY-Plaquette Model in Quantum Lattices

• The XY-plaquette model is a quantum and classical lattice system featuring multi-site plaquette interactions that generate fractonic excitations and diverse phase diagrams.
• Its Hamiltonian combines ring-exchange terms in spatial planes with standard XY-type interactions along the temporal direction, revealing rich duality and unconventional criticality.
• Numerical and analytical studies uncover vortex-wall phases and emergent symmetries, informing research on exotic quantum matter and fractonic field theories.

The XY-plaquette model encompasses a broad family of quantum and classical lattice systems whose Hamiltonians feature multi-spin or multi-boson interactions localized on plaquettes. In particular, the "XY-plaquette model" typically refers to a compact bosonic field, with dynamics governed by ring-exchange (plaquette) interactions on spatial planes, often supplemented by standard XY-type terms along an auxiliary (imaginary time) direction. The model has become pivotal in the study of fractonic field theories, lattice models of Bose liquids, exotic magnets, and quantum dimer models. XY-plaquette models have rich phase diagrams, exhibit unconventional long-range order, and manifest dualities distinct from standard XY models.

1. Lattice Formulations and Hamiltonians

Several concrete realizations exist under the XY-plaquette model umbrella. The minimal 2+1D formulation is defined on a cubic Euclidean lattice with coordinates $(x, y, \tau)$ and a compact phase $\theta_{x,y,\tau} \in (-\pi, \pi]$ at each site. The Hamiltonian is

$$H = -K\sum_{x,y,\tau} \cos\bigl(\Delta_x\Delta_y\,\theta_{x,y,\tau}\bigr) - J\sum_{x,y,\tau} \cos\bigl(\theta_{x,y,\tau+1} - \theta_{x,y,\tau}\bigr),$$

where the mixed difference operator $\Delta_x\Delta_y$ is defined by

$$\Delta_x\Delta_y\,\theta_{x,y,\tau} = [\theta_{x+1,y+1,\tau} - \theta_{x+1,y,\tau}] - [\theta_{x,y+1,\tau} - \theta_{x,y,\tau}].$$

This Hamiltonian encodes a ring-exchange (plaquette) term $K$ in the xy-planes and a standard nearest-neighbor XY-term $J$ along the $\tau$-direction. In the quantum context, related models involve O(2) or SU(2) spins with similar multi-site plaquette couplings (Begun et al., 2024, Biskup et al., 2010).

Alternate formulations arise in studies of valence-bond-ordered phases, such as the spin-1/2 J₁–J₂ XY model on the square lattice: $$H = J_1\sum_{\langle i,j \rangle} (S^x_i S^x_j + S^y_i S^y_j) + J_2\sum_{\langle\langle i,k\rangle\rangle} (S^x_i S^x_k + S^y_i S^y_k),$$ where $S_i^\alpha$ are spin-1/2 operators and the sums run over nearest and next-nearest neighbors (Chan et al., 2022).

2. Physical Motivations and Connections

XY-plaquette interactions appear naturally in several physical settings:

• Fractonic field theories: The continuum limit yields action densities of the form $$\int d^3x\left[(\partial_\tau\phi)^2 + (\partial_x\partial_y\phi)^2\right]$$, producing multipolar conservation laws and fundamentally immobile ('fracton') charges. Compactification ($\phi \to e^{i\theta}$) on the lattice induces vortex excitations, directly leading to the XY-plaquette model structure (Begun et al., 2024, Yoneda, 2022).
• Exciton Bose liquids and cold atomic gases: Four-body ring-exchange arises as a low-energy effective interaction in bosonic systems, realizing the XY-plaquette Hamiltonian in optical lattices (Begun et al., 2024).
• Quantum dimer models: The resonant dynamics of dimers require ring-exchange terms, yielding models mathematically equivalent to certain variants of the XY-plaquette system (Begun et al., 2024).
• Orbital and orbital-compass models: Checkerboard and alternating-plaquette models on the square lattice generalize the XY-plaquette interactions with directionally dependent couplings (Biskup et al., 2010).

3. Vortex Structure, Fractonicity, and Exotic Ordering

The core feature of the XY-plaquette model is the emergence of vortex-like topological excitations governed by unique mobility and conservation constraints:

• Vorticity on plaquettes is defined by

$$v_{P} = \frac{1}{2\pi}\left[\bar{\theta}_{x,y+1,\tau} - \bar{\theta}_{x,y,\tau} - \bar{\theta}_{x+1,y+1,\tau} + \bar{\theta}_{x+1,y,\tau}\right] \in \mathbb{Z},$$

imposing a conservation law (closed vortex loops).

• Phases: Monte Carlo studies identify (i) a disordered, vortex-dominated phase with percolating 3D vortex loops ($D_f \approx 3$), and (ii) a partially disordered "vortex-wall" phase where vortices confine to flat $x\tau$ or $y\tau$ planes, spontaneously breaking internal dipole-like symmetry, with fractal dimension $D_f \approx 2$ (Begun et al., 2024).
• Orientational Long-Range Order (OLRO): Some classical models show spontaneous selection of nematic axes (e.g., spins aligning along x or z directions), breaking rotational symmetry but not yielding net magnetization, due to local $Z_2$ symmetries—an "orientational" rather than magnetic order (Biskup et al., 2010).

4. Duality and Villain Formulation

The XY-plaquette model exhibits a rich web of dualities manifest in multiple path-integral and Hamiltonian formalisms:

• Modified Villain formulation (MVF): Incorporates a dual phase field $\tilde\theta$ and a topological BF-type coupling, making duality manifest at the cost of introducing a second field and an explicit topological term $i\tilde{\theta}\omega$.
• Standard Villain formulation (SVF): Achieves self-duality through Poisson summation and introduces divergence constraints on "magnetic currents," hiding the topological structure in delta-function constraints.
• Dual Hamiltonian method (DHM): Bypasses explicit path integrals and instead matches "pseudo-Josephson" current-voltage relations to derive the same self-dual correspondence (Yoneda, 2022).

These formulations are mathematically equivalent and yield the same mapping of couplings, winding, and momentum currents, rigorously demonstrating the XY-plaquette model's self-duality structure (Yoneda, 2022).

5. Phase Diagrams and Criticality

Extensive Monte Carlo and density-matrix renormalization group (DMRG) studies provide detailed phase diagrams and order parameter characterizations:

• In the minimal 2+1D XY-plaquette model, the diagram features a transition line $J_c(K)$ separating vortex-dominated and vortex-wall phases, with both first-order ($\gamma/\nu=3$) and continuous ($\gamma/\nu<3$) transitions observed, e.g., at representative points $$(J,K)=(1.5,0.757), (3,0.322), (0.757,1.5), (0.359,3)$$ (Begun et al., 2024).
• Along the diagonal $J=K=\beta$, the critical point is $\beta_c=1.00(7)$. For $\beta>\beta_c$, susceptibility plateaus at large values, signifying static vortex-wall order (Begun et al., 2024).
• In the spin-1/2 J₁–J₂ XY model, three phases are realized: Néel AFM ($g<0.50$), a narrow PVB window ($0.50 \leq g \leq 0.54$), and stripy AFM ($g>0.54$), with sharply first-order transitions and absence of quantum spin liquid regimes (Chan et al., 2022).

6. Quantum Extensions and Numerical Methods

Quantum generalizations with SU(2) or O(2) spins admit similar order phenomena at large spin $S$, with bounds on orientational order established using coherent-state and reflection-positivity methods (Biskup et al., 2010). For spin-1/2 models, rigorous results remain inaccessible due to insufficient spin-wave approximations at low $S$.

Advanced numerical algorithms address critical slowing down in simulations, particularly at low $T$:

• Plaquette-flip updates: Augment standard Metropolis or heat-bath algorithms with random “plaquette-flips” (i.e., simultaneous flipping of all components on a plaquette), leveraging exact symmetries to increase sampling efficiency, especially in symmetry-broken phases (Biskup et al., 2010).
• Cluster methods analogous to Swendsen-Wang algorithms for Potts models are effective in deep-ordered regimes by targeting non-local degrees of freedom.

7. Comparison, Broader Impact, and Open Problems

The XY-plaquette model is distinguished from standard XY and Heisenberg models by the nature of its ordering and excitations:

• Absence of magnetic or quantum spin liquid order: No region is found with all conventional (magnetic or dimer) order parameters vanishing; crystal-like symmetry breaking (plaquette valence bond or OLRO) prevails where intermediate spin-liquid phases might appear in Heisenberg analogues (Chan et al., 2022).
• Highly degenerate ground states with local Z₂ symmetries: Leads to fundamentally different low-temperature physics, including unconventional criticality and fractonic behavior.
• Emergent symmetry breaking: Internal polynomial shift symmetry or 90° rotation symmetry can be spontaneously broken, resulting in domain wall formation in the vortex sector or nematic-type orientational transitions.

Ongoing research is focused on:

• Extending rigorous results to spin-1/2 quantum systems;
• Fully classifying universality classes of observed transitions;
• Connecting lattice XY-plaquette models to continuum fractonic and higher-rank gauge theories;
• Investigating the interplay of restricted mobility, topology, and quantum entanglement in engineered cold-atom systems and quantum materials (Begun et al., 2024, Chan et al., 2022, Yoneda, 2022, Biskup et al., 2010).