We introduce quantum circuits in two and three spatial dimensions which are classically simulable, despite producing a high degree of operator entanglement. We provide a partial characterization of these "automaton" quantum circuits, and use them to study operator growth, information spreading, and local charge relaxation in quantum dynamics with subsystem symmetries, which we define as overlapping symmetries that act on lower-dimensional submanifolds. With these symmetries, we discover the anomalous subdiffusion of conserved charges; that is, the charges spread slower than diffusion in the dimension of the subsystem symmetry. By studying an effective operator hydrodynamics in the presence of these symmetries, we predict the charge autocorrelator to decay (i) as log(t)/ √ t in two dimensions with a conserved U (1) charge along intersecting lines, and (ii) as 1/t 3/4 in three spatial dimensions with intersecting planar U (1) symmetries. Through large-scale studies of automaton dynamics with these symmetries, we numerically observe charge relaxation that is consistent with these predictions. In both cases, the spatial charge distribution is distinctly non-Gaussian, and reminiscent of the diffusion of charges along a fractal surface. We numerically study the onset of quantum chaos in the spreading of local operators under these automaton dynamics, and observe power-law broadening of the ballistically-propagating fronts of evolving operators in two and three dimensions, and the saturation of out-of-time-ordered correlations to values consistent with quantum chaotic behavior.

In this paper, we explore the relationship between strong spin-orbit coupling and spin liquid physics. We study a very general model on the triangular lattice where spin-orbit coupling leads to the presence of highly anisotropic interactions. We use variational Monte Carlo to study both U(1) quantum spin liquid states and ordered ones, via the Gutzwiller projected fermion construction. We thereby obtain the ground state phase diagram in this phase space. We furthermore consider effects beyond the Gutzwiller wavefunctions for the spinon Fermi surface quantum spin liquid, which are of particular importance when spin-orbit coupling is present.

By developing a method to represent the Renyi entropies via a replica trick on classical statistical mechanical systems, we introduce a procedure to calculate the Renyi mutual information (RMI) in any Monte Carlo simulation. Through simulations on several classical models, we demonstrate that the RMI can detect finite-temperature critical points, and even identify their universality class, without knowledge of an order parameter or other thermodynamic estimators. Remarkably, in addition to critical points mediated by symmetry breaking, the RMI is able to detect topological vortex-unbinding transitions, as we explicitly demonstrate on simulations of the XY model.

No abstract