X-Cube Model: Fracton Topological Order

• X-Cube Model is a 3D lattice system displaying type-I fracton order, defined on a cubic lattice with qubits on links and commuting stabilizers.
• It exhibits unique mobility constraints where fractons are immobile, lineons move along one axis, and planons are mobile within planes.
• The model features subextensive ground-state degeneracy that intertwines lattice geometry, foliation, and boundary conditions, influencing phase transitions.

The X-cube model is a canonical exactly solvable lattice model of type-I fracton topological order. In its standard form, it is defined on a three-dimensional cubic lattice with qubits on links and commuting stabilizer terms attached to cubes and planar vertex “crosses.” Its defining features are subextensive ground-state degeneracy, immobile fractons, lineons constrained to one-dimensional motion, and planons mobile only within planes. Unlike conventional topological phases, its long-distance structure depends not only on topology but also on geometry, lattice foliation, and boundary conditions (Zhou et al., 2022, Slagle et al., 2017).

1. Lattice definition and stabilizer structure

A common lattice presentation places qubits on the edges of a cubic lattice and uses the Hamiltonian

$$H=-\sum_{v\in V}\left(A_v^x+A_v^y+A_v^z\right)-\sum_{c\in C} B_c.$$

Here $A_v^i$ applies Pauli $Z$ on the four edges perpendicular to the $i$-axis and incident to $v$, while $B_c$ applies Pauli $X$ on the twelve edges of cube $c$. These operators satisfy

$(A_v^i)^2=B_c^2=1,\qquad [A_v^i,B_c]=0,$

so the model is an exactly solvable commuting-projector stabilizer code. The ground state may be written as

$$|GS\rangle=\prod_c \frac{1+B_c}{2}\,|\phi_0\rangle, \qquad |\phi_0\rangle=|00\cdots 0\rangle,$$

with the vertex projectors omitted because $A_v^i$0 acts trivially on $A_v^i$1 in this convention (Chen et al., 2022).

A $A_v^i$2 formulation uses generalized clock operators $A_v^i$3 and $A_v^i$4 obeying

$A_v^i$5

and writes the lattice Hamiltonian as

$A_v^i$6

In that notation, $A_v^i$7 is a cube operator and $A_v^i$8 is a vertex operator orthogonal to direction $A_v^i$9. The local relation

$Z$0

encodes the redundancy among the three vertex terms at a site and is central to the excitation and constraint structure of the model (Slagle et al., 2017, Lake et al., 2021).

The model is often described as the representative exactly solvable realization of type-I fracton order because the commuting stabilizers encode restricted-mobility excitations without invoking conventional symmetry breaking. This suggests that the stabilizer algebra, rather than a continuum order parameter, is the primary organizing principle of the phase (Zhou et al., 2022).

2. Excitations and mobility constraints

The elementary excitations fall into three standard classes. Fractons are cube-term violations and are completely immobile as isolated excitations. Lineons arise from violations of vertex-cross terms and come in three species, each mobile only along one Cartesian axis. Planons are bound states, typically of two fractons separated along a fixed axis, and can move within the plane perpendicular to that axis (Zhou et al., 2022).

The mobility constraints are tied to the geometry of the creation operators. A membrane operator

$Z$1

creates four fractons at the corners of a square membrane, while a straight string operator

$Z$2

creates lineons at its endpoints. If the string bends, extra lineons appear at the turning points. The absence of spatially deformable string-like operators is therefore the direct origin of subdimensional motion: moving a fracton or lineon in a forbidden direction necessarily creates additional excitations and costs finite energy (Zhou et al., 2022).

In the $Z$3 description, lineons may be denoted $Z$4, with motion constrained respectively to the $Z$5, $Z$6, and $Z$7 directions. Their fusion rule is

$Z$8

Lineon dipoles $Z$9 and fracton dipoles $i$0 are planons: they move freely in planes but not in the direction of their dipole moment. The continuum field-theory analysis encodes the same kinematics through higher-derivative charge constraints, with fracton density $i$1 and lineon density $i$2 obeying generalized continuity equations rather than ordinary particle-number conservation (Lake et al., 2021, Slagle et al., 2017).

These mobility constraints are robust to arbitrary local perturbations in the sense emphasized by the continuum description. A plausible implication is that the defining content of the X-cube phase is not merely the existence of unusual excitations, but the higher-rank kinematic structure that forbids generic transport processes (Slagle et al., 2017).

3. Ground-state degeneracy, foliation, and geometric order

On an $i$3 periodic lattice, the $i$4 X-cube model has ground-state degeneracy

$i$5

For $i$6, this becomes

$i$7

The $i$8 reflects three global redundancies among the nonlocal operators. This subextensive scaling is one of the model’s defining features and distinguishes it sharply from conventional three-dimensional topological orders such as the 3D toric code, whose degeneracy is finite and topology-only (Slagle et al., 2017, Mühlhauser et al., 2021).

A central conceptual result is that this degeneracy is not purely topological. The continuum analysis shows that even on a manifold with trivial topology, spatial curvature can induce a robust degeneracy. In a curved cubic lattice with angular defects and “straight loops,” the simplest $i$9 construction yields a $v$0 degeneracy, and with $v$1 layers this becomes $v$2. The splitting under local perturbations is exponentially small in the size of the curved region because connecting different sectors requires a large nonlocal loop operator. Geometry, not topology alone, therefore protects part of the degeneracy (Slagle et al., 2017).

The generic-lattice construction makes this dependence explicit. Starting from intersecting two-dimensional “i-surfaces” in three dimensions, one defines vertices where three i-surfaces intersect, links where two intersect, and 3-cells as enclosed volumes. The resulting generalized X-cube models inherit their mobility constraints from the geometry of these surfaces: dimension-1 particles move along intersection lines, and dimension-2 particles move on the i-surfaces themselves. The degeneracy scales as

$v$3

where $v$4 is the genus of an orientable i-surface $v$5. On this basis, the authors propose that fracton orders should be regarded as a geometric order (Slagle et al., 2017).

The foliation perspective refines this picture. The X-cube model contains hidden layers of 2+1D topological states and can be viewed through a coupled-layer construction of stacked toric codes. On the 3-torus,

$v$6

with the $v$7 arising because Wilson loops of the minimal intersecting toric-code layers are tied together along their intersections. Screw dislocations provide a direct probe of this foliation structure: they change $v$8 by a finite amount and also induce tunneling of subdimensional excitations along the defect line, revealing structure that is not reducible to a simple stack of 2D layers (Manoj et al., 2020).

4. Continuum field theory and boundary structure

The X-cube model admits a continuum field-theory description, but this theory is not and cannot be a topological quantum field theory. The reason is that X-cube order is not invariant under arbitrary smooth spacetime deformations; instead, the continuum action is invariant only under the restricted “subconformal” transformations

$v$9

This restricted invariance reflects the anisotropic locality of the lattice model and the fact that bending coordinate axes changes the allowed mobility of excitations (Slagle et al., 2017).

In one continuum formulation, the Lagrangian takes the form

$B_c$0

with the constraint

$B_c$1

The gauge transformations are

$B_c$2

together with temporal-component shifts. These relations reproduce the lattice commutation algebra, braiding-like statistical phases, and the system-size-dependent degeneracy (Slagle et al., 2017).

Boundary physics is comparably rich. For the $B_c$3 X-cube model on a $B_c$4 constant surface, the low-energy boundary theory reduces to a generalized two-component $B_c$5-matrix theory,

$B_c$6

with a momentum subsystem symmetry

$B_c$7

The double spatial derivative in the kinetic term is the boundary signature of subsystem structure. The theory reproduces a subset of the bulk exchange statistics, particularly those involving mobile composites such as fracton dipoles and lineons (Luo et al., 2022).

The boundary admits multiple gapped terminations. On $B_c$8, smooth boundaries on both ends give

$B_c$9

rough boundaries on both ends give

$X$0

and one smooth plus one rough boundary gives

$X$1

Anisotropic and dyonic boundaries also exist, and the anomaly inflow analysis shows that the boundary ’t Hooft anomaly is not canceled uniquely by the X-cube bulk. This suggests that the boundary theory probes only part of the three-dimensional fractonic data (Luo et al., 2022).

5. Dynamics, finite-temperature behavior, and phase transitions

Perturbing the model by Zeeman fields exposes distinctive dynamical signatures of subdimensional motion. For the Hamiltonian

$X$2

large-scale quantum Monte Carlo and stochastic analytic continuation show that the fracton phase exhibits strong anisotropy in both real-space correlations and dynamical structure factors. For example, in the fracton phase at $X$3, the nearest-neighbor lineon correlation satisfies

$X$4

The lineon spectral function disperses only along the allowed direction, whereas the fracton spectral function is essentially flat in momentum and remains gapped, with a peak energy around $X$5. A first-order transition into a trivial paramagnetic phase occurs near

$X$6

depending on the perturbation channel (Zhou et al., 2022).

At finite temperature, the equilibrium thermodynamics are unexpectedly simple. Exact partition-function calculations and bond-algebraic dualities show that the X-Cube model and its $X$7 generalization have no finite-temperature phase transitions. For open boundaries, the partition function is

$X$8

with $X$9 and $c$0, and the thermodynamic free-energy density agrees with the cylindrical and periodic cases. Under a solvable Glauber-type bath, autocorrelations decay exponentially with finite rates at all nonzero temperatures, so the model is thermally fragile rather than glassy in the sense of system-size-dependent positive-temperature memory (Weinstein et al., 2018).

Quantum phase transitions out of the X-cube phase depend strongly on which excitations condense. Condensing fracton dipoles or lineon dipoles produces stacked deconfined gauge theories; condensing one lineon species yields a stack of 2D $c$1 gauge theories; condensing two or more lineon species yields a trivial paramagnet. For $c$2, the dipole-condensation transitions can be governed by stable critical points described by decoupled stacks of $c$3 conformal field theories, while lineon condensation can produce a gapless intermediate phase described as an array of $c$4 conformal field theories (Lake et al., 2021).

A complementary route through the phase diagram arises when the X-cube model competes with the 3D toric code. In the combined Hamiltonian, the exact commuting decomposition

$c$5

allows a detailed analysis. Besides the toric-code and X-cube phases, two further phases appear, both adiabatically connected to classical limits with nontrivial sub-extensive degeneracies. All phase transitions are found to be first order (Mühlhauser et al., 2021).

6. Realizations, algorithms, defects, and extensions

The X-cube ground state can be prepared exactly by a Clifford-only quantum circuit. On an $c$6 cubic lattice on a 3D torus, the stated preparation depth is

$c$7

The construction reformulates the projector application problem geometrically, using Hadamard and CNOT gates in a layered pattern. A related gluing method extends the preparation strategy via measurements and membrane-based correction operators, reflecting the fact that X-cube excitations are fractonic rather than freely mobile (Chen et al., 2022).

A dynamical measurement-based realization is provided by the X-Cube Floquet code. Built from intersecting 2D 4.8.8 Floquet-code layers in the $c$8, $c$9, and $(A_v^i)^2=B_c^2=1,\qquad [A_v^i,B_c]=0,$0 directions, it uses a period-six sequence of local weight-two measurements:

1. yellow checks and on-site checks,
2. blue checks,
3. green checks,
4. yellow checks,
5. blue checks,
6. green checks. Within one Floquet period, the codespace alternates between X-cube fracton order and layers of entangled 2D toric codes. The encoded logical-qubit count is

$(A_v^i)^2=B_c^2=1,\qquad [A_v^i,B_c]=0,$1

and the model is argued to have a non-zero error threshold. The same work also gives a two-body Hamiltonian realization of the coupled-layer X-cube limit (Zhang et al., 2022).

Defects and non-Euclidean geometries further enlarge the family of X-cube-like phenomena. Screw dislocations change $(A_v^i)^2=B_c^2=1,\qquad [A_v^i,B_c]=0,$2 by a finite amount and can enable fractons or lineons to tunnel along the defect line, exposing both the foliated structure and genuinely fractonic mobility effects (Manoj et al., 2020). On hyperbolic lattices embedded in $(A_v^i)^2=B_c^2=1,\qquad [A_v^i,B_c]=0,$3, the Y-cube generalization replaces some X-shaped vertex terms by Y-shaped ones and supports treeons, excitations confined to a fractal tree rather than a line. In flat limits, these treeons reduce to either lineons or planeons, showing explicitly that subdimensional mobility can itself be geometry-dependent (Yan et al., 2022).

Taken together, these constructions show that the X-cube model functions both as a fixed exactly solvable phase and as a template for a broader class of geometry-sensitive fracton systems. This suggests that its most durable significance lies in unifying stabilizer-code realizations, continuum higher-rank gauge structure, and lattice-geometric mechanisms of restricted mobility within a single paradigm (Slagle et al., 2017).