
Unidirectional subsystem symmetry in a hole-doped honeycomb-lattice Ising magnet
Sambuddha Sanyal,1, ∗Alexander Wietek,2, 3, †and John Sous4, ‡
1Department of Physics, Indian Institute of Science Education and Research (IISER) Tirupati, Tirupati 517507, India
2Center for Computational Quantum Physics, Flatiron Institute, 162 Fifth Avenue, New York, NY 10010, USA
3Max Planck Institute for the Physics of Complex Systems,
N¨othnitzer Strasse 38, Dresden 01187, Germany
4Department of Physics, Stanford University, Stanford, CA 93405, USA
(Dated: October 4, 2022)
We study a model of a hole-doped collinear Ising antiferromagnet on the honeycomb lattice as a
route toward realization of subsystem symmetry. We find nearly exact conservation of dipole sym-
metry verified both numerically with exact diagonalization (ED) on finite clusters and analytically
with perturbation theory. The emergent symmetry forbids the motion of single holes – or fractons
– but allows hole pairs – or dipoles – to move freely along a one-dimensional line, the antiferromag-
netic direction, of the system; in the transverse direction both fractons and dipoles are completely
localized. This presents a realization of a ‘unidirectional’ subsystem symmetry. By studying inter-
actions between dipoles, we argue that the subsystem symmetry is likely to continue to persist up
to finite (but probably small) hole concentrations.
Introduction. Fractons are the newest addition in the
line of exotic quasiparticles in condensed matter physics.
Recent years have witnessed tremendous progress in un-
derstanding fractons [1–40], which unmasked connec-
tions with other areas of physics including topological
order [2,4,8,9], gauge theory [12,13], quantum comput-
ing [4], glasses and soft matter [21,22], These exotic prop-
erties result from unusual mobility constraints whereby
a fracton is a charge-like excitation which has restricted
mobility when in isolation, but which, nonetheless, can
easily move in a subdimension of space when bound to an
oppositely charged fracton in a dipolar bound state [41].
The realization of these aberrant mobility constraints
and subsystem symmetry in physical systems is, however,
not naturally available, stimulating many interesting the-
oretical proposals [42–44]. One particularly appealing
approach is based on the idea that defect motion-induced
frustration in an ordered background leads to emergent
immobility constraints on the motion of defects them-
selves which in turn gives rise to fractonic quasiparticles,
with hole-doped Ising antiferromagnets [45–54] being a
prime example of this mechanism [27,28]. Unlike other
constructions proposed to realize fracton topological or-
der (e.g. [55]), this approach does not require an exten-
sive number of locally conserved quantities (which results
in an extensively degenerate ground state), making the
physical realization of fracton conservation laws (with-
out topological order) in quantum magnets significantly
more accessible. However, the emergence of fracton-like
quasiparticles in fully two-dimensional (2D) hole-doped
Ising antiferromagnets has been demonstrated only in the
asymptotic limit of tJ(where tis the hole hopping
and Jis the Ising coupling) [27], where already at the
leading order in t2/J, the direction of the dipole moment
is not conserved even though its magnitude is because
hole pairs can rotate as they move in the 2D plane.
In this letter we overcome this limitation and pro-
FIG. 1. Ground state of the Hamiltonian HIsing (Eq. (2))
with one electron per site. The three bonds of the honey-
comb lattice are denoted as x0, y0and z0respectively (xand
yare the Cartesian directions). The Hamiltonian describes
an Ising magnet on a honeycomb lattice with antiferromag-
netic exchange along the x0and y0bonds and ferromagnetic
exchange along the z0bonds. The vectors δx0,δy0and δz0
connect near-neighbor sites in the x0,y0and z0directions, re-
spectively. The unit cell of the honeycomb lattice with A and
B sublattices is shown in the dotted region, and a1and a2are
the primitive lattice vectors.
pose an essentially exact physical realization of fractonic
quasiparticles with subsystem dipolar symmetry in a 2D
non-degenerate, ordered spin system with local two-spin
interactions. The key ingredient in our proposal is the
collinear antiferromagnetic order in which defect motion-
induced frustration completely prevents hole pairs or
dipoles from moving in the perpendicular direction, re-
sulting in exact conservation of both the magnitude and
direction of dipole moment in the asymptotic limit tJ.
This further endows the system with a subsystem symme-
try since dipole motion is restricted to a one-dimensional
(1D) submanifold of the system. This symmetry mani-
fests only along the antiferromagnetic direction, thus we
arXiv:2210.00012v1 [cond-mat.str-el] 30 Sep 2022