No Arabic abstract
We construct models hosting classical fractal spin liquids on two realistic three-dimensional (3D) lattices of corner-sharing triangles: trillium and hyperhyperkagome (HHK). Both models involve the same form of three-spin Ising interactions on triangular plaquettes as the Newman-Moore (NM) model on the 2D triangular lattice. However, in contrast to the NM model and its 3D generalizations, their degenerate ground states and low-lying excitations cannot be described in terms of scalar cellular automata (CA), because the corresponding fractal structures lack a simplifying algebraic property, often termed the Freshmans dream. By identifying a link to matrix CAs -- that makes essential use of the crystallographic structure -- we show that both models exhibit fractal symmetries of a distinct class to the NM-type models. We devise a procedure to explicitly construct low-energy excitations consisting of finite sets of immobile defects or fractons, by flipping arbitrarily large self-similar subsets of spins, whose fractal dimensions we compute analytically. We show that these excitations are associated with energetic barriers which increase logarithmically with system size, leading to fragile glassy dynamics, whose existence we confirm via classical Monte Carlo simulations. We also discuss consequences for spontaneous fractal symmetry breaking when quantum fluctuations are introduced by a transverse magnetic field, and propose multi-spin correlation function diagnostics for such transitions. Our findings suggest that matrix CAs may provide a fruitful route to identifying fractal symmetries and fracton-like behaviour in lattice models, with possible implications for the study of fracton topological order.
We introduce fractal liquids by generalizing classical liquids of integer dimensions $d = 1, 2, 3$ to a fractal dimension $d_f$. The particles composing the liquid are fractal objects and their configuration space is also fractal, with the same non-integer dimension. Realizations of our generic model system include microphase separated binary liquids in porous media, and highly branched liquid droplets confined to a fractal polymer backbone in a gel. Here we study the thermodynamics and pair correlations of fractal liquids by computer simulation and semi-analytical statistical mechanics. Our results are based on a model where fractal hard spheres move on a near-critical percolating lattice cluster. The predictions of the fractal Percus-Yevick liquid integral equation compare well with our simulation results.
Close-packed, classical dimer models on three-dimensional, bipartite lattices harbor a Coulomb phase with power-law correlations at infinite temperature. Here, we discuss the nature of the thermal phase transition out of this Coulomb phase for a variety of dimer models which energetically favor crystalline dimer states with columnar ordering. For a family of these models we find a direct thermal transition from the Coulomb phase to the dimer crystal. While some systems exhibit (strong) first-order transitions in correspondence with the Landau-Ginzburg-Wilson paradigm, we also find clear numerical evidence for continuous transitions. A second family of models undergoes two consecutive thermal transitions with an intermediate paramagnetic phase separating the Coulomb phase from the dimer crystal. We can describe all of these phase transitions in one unifying framework of candidate field theories with two complex Ginzburg-Landau fields coupled to a U(1) gauge field. We derive the symmetry-mandated Ginzburg-Landau actions in these field variables for the various dimer models and discuss implications for their respective phase transitions.
A generalized set of Clifford cellular automata, which includes all Clifford cellular automata, result from the quantization of a lattice system where on each site of the lattice one has a $2k$-dimensional torus phase space. The dynamics is a linear map in the torus variables and it is also local: the evolution depends only on variables in some region around the original lattice site. Moreover it preserves the symplectic structure. These are classified by $2ktimes 2k$ matrices with entries in Laurent polynomials with integer coefficients in a set of additional formal variables. These can lead to fractal behavior in the evolution of the generators of the quantum algebra. Fractal behavior leads to non-trivial Lyapunov exponents of the original linear dynamical system. The proof uses Fourier analysis on the characteristic polynomial of these matrices.
We study matrix product unitary operators (MPUs) for fermionic one-dimensional (1D) chains. In stark contrast with the case of 1D qudit systems, we show that (i) fermionic MPUs do not necessarily feature a strict causal cone and (ii) not all fermionic Quantum Cellular Automata (QCA) can be represented as fermionic MPUs. We then introduce a natural generalization of the latter, obtained by allowing for an additional operator acting on their auxiliary space. We characterize a family of such generalized MPUs that are locality-preserving, and show that, up to appending inert ancillary fermionic degrees of freedom, any representative of this family is a fermionic QCA and viceversa. Finally, we prove an index theorem for generalized MPUs, recovering the recently derived classification of fermionic QCA in one dimension. As a technical tool for our analysis, we also introduce a graded canonical form for fermionic matrix product states, proving its uniqueness up to similarity transformations.
The present paper considers some classical ferromagnetic lattice--gas models, consisting of particles that carry $n$--component spins ($n=2,3$) and associated with a $D$--dimensional lattice ($D=2,3$); each site can host one particle at most, thus implicitly allowing for hard--core repulsion; the pair interaction, restricted to nearest neighbors, is ferromagnetic, and site occupation is also controlled by the chemical potential $mu$. The models had previously been investigated by Mean Field and Two--Site Cluster treatments (when D=3), as well as Grand--Canonical Monte Carlo simulation in the case $mu=0$, for both D=2 and D=3; the obtained results showed the same kind of critical behaviour as the one known for their saturated lattice counterparts, corresponding to one particle per site. Here we addressed by Grand--Canonical Monte Carlo simulation the case where the chemical potential is negative and sufficiently large in magnitude; the value $mu=-D/2$ was chosen for each of the four previously investigated counterparts, together with $mu=-3D/4$ in an additional instance. We mostly found evidence of first order transitions, both for D=2 and D=3, and quantitatively characterized their behaviour. Comparisons are also made with recent experimental results.