Surface growth in Z2 lattice gauge theories

We numerically study the mean-field dynamics of a (2+1)D spin system with thousands of spins and show that experimentally-feasible two-body Ising interactions can stabilize a prethermal Z2 lattice gauge structure with dynamical matter, manifested by a gauge-invariant plateau with exponentially long lifetime. Eventually, the metastable prethermal Z2 gauge structure breaks down via a proliferation of Gauss' law defects, similar to bubble formation in false vacuum decay. In this regime, we discover spatio-temporal correlations described by a non-linear surface growth consistent with the (1+1)D Kardar-Parisi-Zhang (KPZ) universality class. We benchmark our results in small systems against semi-classical discrete time Wigner approximation (DTWA) and exact diagonalization (ED), where the breakdown of DTWA signals the emergence of an extensive number of local symmetries that strongly influence the thermalization pathway. Our model provides a testbed for quantum simulators and is directly implementable in large-scale arrays of Rydberg atoms.