Plaquette Ising Model Research Articles

Foliated order parameter in a fracton phase transition

Finding suitable indicators for characterizing quantum phase transitions plays an important role in understanding different phases of matter. It is especially important for fracton phases where a combination of topology and fractionalization leads to exotic features not seen in other known quantum phases. In this paper, we consider the above problem by studying phase transition in the X-cube model in the presence of a nonlinear perturbation. Using an analysis of the ground state fidelity and identifying a discontinuity in the global entanglement, we show there is a first-order quantum phase transition from a type-I fracton phase with a highly entangled nature to a magnetized phase. Accordingly, we conclude that the global entanglement, as a measure of the total quantum correlations in the ground state, can well capture certain features of fracton phase transitions. Then, we introduce a nonlocal order parameter in the form of a foliated operator which can characterize the above phase transition. We particularly show that such an order parameter has a geometric nature which captures specific differences of fracton phases with topological phases. Our study is specifically based on a well-known dual mapping to the classical plaquette Ising model where it shows the importance of such dualities in studying different quantum phases of matter.

Linked cluster expansions via hypergraph decompositions.

We propose a hypergraph expansion which facilitates the direct treatment of quantum spin models with many-site interactions via perturbative linked cluster expansions. The main idea is to generate all relevant subclusters and sort them into equivalence classes essentially governed by hypergraph isomorphism. Concretely, a reduced König representation of the hypergraphs is used to make the equivalence relation accessible by graph isomorphism. During this procedure we determine the embedding factor for each equivalence class, which is used in the final resummation in order to obtain the final result. As an instructive example we calculate the ground-state energy and a particular excitation gap of the plaquette Ising model in a transverse field on the three-dimensional cubic lattice.

(Four) Dual Plaquette 3D Ising Models.

A characteristic feature of the plaquette Ising model is its planar subsystem symmetry. The quantum version of this model has been shown to be related via a duality to the X-Cube model, which has been paradigmatic in the new and rapidly developing field of fractons. The relation between the plaquette Ising and the X-Cube model is similar to that between the quantum transverse spin Ising model and the Toric Code. Gauging the global symmetry in the case of the Ising model and considering the gauge invariant sector of the high temperature phase leads to the Toric Code, whereas gauging the subsystem symmetry of the quantum transverse spin plaquette Ising model leads to the X-Cube model. A non-standard dual formulation of the plaquette Ising model which utilises three flavours of spins has recently been discussed in the context of dualising the fracton-free sector of the X-Cube model. In this paper we investigate the classical spin version of this non-standard dual Hamiltonian and discuss its properties in relation to the more familiar Ashkin–Teller-like dual and further related dual formulations involving both link and vertex spins and non-Ising spins.

Correlation function diagnostics for type-I fracton phases

Fracton phases are recent entrants to the roster of topological phases in three dimensions. They are characterized by subextensively divergent topological degeneracy and excitations that are constrained to move along lower dimensional subspaces, including the eponymous fractons that are immobile in isolation. We develop correlation function diagnostics to characterize Type I fracton phases which build on their exhibiting {\it partial deconfinement}. These are inspired by similar diagnostics from standard gauge theories and utilize a generalized gauging procedure that links fracton phases to classical Ising models with subsystem symmetries. En route, we explicitly construct the spacetime partition function for the plaquette Ising model which, under such gauging, maps into the X-cube fracton topological phase. We numerically verify our results for this model via Monte Carlo calculations.

Exact solutions to plaquette Ising models with free and periodic boundaries

An anisotropic limit of the 3d plaquette Ising model, in which the plaquette couplings in one direction were set to zero, was solved for free boundary conditions by Suzuki (1972) [1], who later dubbed it the fuki-nuke, or “no-ceiling”, model. Defining new spin variables as the product of nearest-neighbour spins transforms the Hamiltonian into that of a stack of (standard) 2d Ising models and reveals the planar nature of the magnetic order, which is also present in the fully isotropic 3d plaquette model. More recently, the solution of the fuki-nuke model was discussed for periodic boundary conditions, which require a different approach to defining the product spin transformation, by Castelnovo et al. (2010) [2].We clarify the exact relation between partition functions with free and periodic boundary conditions expressed in terms of original and product spin variables for the 2d plaquette and 3d fuki-nuke models, noting that the differences are already present in the 1d Ising model. In addition, we solve the 2d plaquette Ising model with helical boundary conditions. The various exactly solved examples illustrate how correlations can be induced in finite systems as a consequence of the choice of boundary conditions.

Macroscopic degeneracy and order in the 3D plaquette Ising model

The purely plaquette 3D Ising Hamiltonian with the spins living at the vertices of a cubic lattice displays several interesting features. The symmetries of the model lead to a macroscopic degeneracy of the low-temperature phase and prevent the definition of a standard magnetic order parameter. Consideration of the strongly anisotropic limit of the model suggests that a layered, “fuki-nuke” order still exists and we confirm this with multi-canonical simulations. The macroscopic degeneracy of the low-temperature phase also changes the finite-size scaling corrections at the first-order transition in the model and we see this must be taken into account when analyzing our measurements.

Planar ordering in the plaquette-only gonihedric Ising model

In this paper we conduct a careful multicanonical simulation of the isotropic 3d plaquette (“gonihedric”) Ising model and confirm that a planar, fuki-nuke type order characterises the low-temperature phase of the model. From consideration of the anisotropic limit of the model we define a class of order parameters which can distinguish the low- and high-temperature phases in both the anisotropic and isotropic cases. We also verify the recently voiced suspicion that the order parameter like behaviour of the standard magnetic susceptibility χm seen in previous Metropolis simulations was an artefact of the algorithm failing to explore the phase space of the macroscopically degenerate low-temperature phase. χm is therefore not a suitable order parameter for the model.

Stretched exponentials and tensionless glass in the plaquette Ising model

Using Monte Carlo simulations, we show that the autocorrelation function C(t) in the d = 3 Ising model with a plaquette interaction has a stretched-exponential decay in a supercooled liquid phase. Such a decay characterizes also some ground-state probability distributions obtained from the numerically exact counting of up to 10(450) configurations. A related model with a strongly degenerate ground state but lacking glassy features does not exhibit such a decay. Although the stretched exponential decay of C(t) in the three-dimensional supercooled liquid is inconsistent with the droplet model, its modification that considers tensionless droplets might explain such a decay. An indication that tensionless droplets might play some role comes from the analysis of low-temperature domains that compose the glassy state. It shows that the energy of a domain of size l scales as l(1.15), hence these domains are indeed tensionless.

Nature of the glassy transition in simulations of the ferromagnetic plaquette Ising model

The homogeneous plaquette Ising model in two and three dimensions is investigated by means of Monte Carlo simulations. By introducing a suitable order parameter for the two-dimensional lattice, and the finite-size scaling of the corresponding fourth-order cumulant, it is found that, consistent with the previous theoretical indications, the model in two dimensions is disordered at finite temperature and exhibits a zero-temperature phase transition characteristic of the one-dimensional Ising model with an essential (exponential) singularity of the order-parameter susceptibility as opposed to a Curie-law (power-law) divergence. In three dimensions, however, the model is believed to have a first-order phase transition at Tc approximately 3.6 screened by strong metastability leading to a so-called "glassy transition" at T approximately 3.4 when subjected to slow cooling. By computing the configurational entropy Sc identical withS(liquid)-S(crystal) in the supercooled temperature range via thermodynamic integration of the internal energy results, the Kauzmann temperature defined as that temperature where the extrapolated configurational entropy Sc(T) vanishes, is estimated to be TK approximately 3.18 . By finding ways to estimate the equilibration time of the supercooled liquid and the nucleation time of the stable crystal droplets, it is shown that T approximately 3.4 is indeed the limit of stability or the effective spinodal temperature Tsp, at which the two time-scales associated with the quasiequilibration of the supercooled liquid, taueq, and the nucleation of the stable crystal droplets, taunuc, cross one another, with the former rising above the latter such that the supercooled liquid state becomes physically irrelevant below Tsp approximately 3.4 and the impending entropy crisis at TK approximately 3.18 (<Tsp) is thus avoided. Hence, what is sometimes called "glassy temperature," is really a kinetic spinodal temperature that may be regarded as the remnant of the mean-field spinodal.