Let denote the dual edge variable. The PEPO implements the domain-wall map
Because edge domain walls arise from vertex spins, their product around every plaquette is constrained:
Thus the transverse-field Ising operators map to the lattice-gauge operators, while the image is projected into the positive eigenspace of all . At the commuting-projector fixed point, adjoining these automatic projectors gives
the toric code Hamiltonian. Conversely, solving the zero-flux constraint writes locally and recovers the Ising variables, establishing the duality on the supported symmetry sectors.
Solved by gpt-5.6-sol high.
The original Ising symmetry is the global spin flip
The dual system has a one-form symmetry generated by products of edge Pauli operators along closed noncontractible loops; equivalently, closed Wilson loops label its electric and magnetic topological sectors.
The PEPO annihilates the nonsymmetric Ising sector because each virtual index is summed with the global parity constraint. On the dual side it has support only on configurations obeying all contractible zero-flux constraints and, for a fixed untwisted PEPO, one choice of noncontractible loop eigenvalues. Therefore the duality is invertible only after restricting both Hilbert spaces to corresponding symmetry sectors.
On a torus the toric code has four ground states distinguished by two independent noncontractible loop eigenvalues. A single untwisted duality operator reaches only one of them; inserting the two possible Ising twists supplies the other three. This explains why the global-symmetry Ising description and topologically degenerate toric-code description do not contradict each other.
Solved by gpt-5.6-sol high.