Rokhsar-Kivelson Point Research Articles

Fully packed quantum loop model on the square lattice: Phase diagram and application for Rydberg atoms

The quantum dimer and loop models attract great attentions, partially because the fundamental importance in the phases and phase transitions emerging in these prototypical constrained systems, and partially due to their intimate relevance toward the on-going experiments on Rydberg atom arrays in which the blockade mechanism naturally enforces the local constraint. Here we show, by means of the sweeping cluster quantum Monte Carlo method, the complete ground state phase diagram of the fully packed quantum loop model on the square lattice. We find between the lattice nematic (LN) phase with strong dimer attraction and the staggered phase (SP) with strong dimer repulsion, there emerges a resonating plaquette (RP) phase with off-diagonal translational symmetry breaking. Such a quantum phase is separated from the LN via a first order transition and from the SP by the famous Rokhsar-Kivelson point. Our renormalization group analysis reveals the different flow directions, fully consistent with the order parameter histogram in Monte Carlo simulations. The realization and implication of our phase diagram in Rydberg experiments are proposed.

Classical vertex model dualities in a family of two-dimensional frustrated quantum antiferromagnets

We study a general class of easy-axis spin models on a lattice of corner sharing even-sided polygons with all-to-all interactions within a plaquette. The low energy description corresponds to a quantum dimer model on a dual lattice of even coordination number with a multi dimer constraint. At an appropriately constructed frustration-free Rokhsar-Kivelson (RK) point, the ground state wavefunction can be exactly mapped onto a classical vertex model on the dual lattice. When the dual lattice is bipartite, the vertex models are bonded and are self dual under Wegner's duality, with the self dual point corresponding to the RK point of the original multi-dimer model. We argue that the self dual point is a critical point based on known exact solutions to some of the vertex models. When the dual lattice is non-bipartite, the vertex model is arrowed, and we use numerical methods to argue that there is no phase transition as a function of the vertex weights. Motivated by these wavefunction dualities, we construct two other distinct families of frustration-free Hamiltonians whose ground states can be mapped onto these vertex models. Many of these RK Hamiltonians provably host $\mathbb{Z}_2$ topologically ordered phases.

Reentrance effect in the high-temperature critical phase of the quantum dimer model on the square lattice

We present a quantum Monte Carlo investigation of the finite-temperature phase diagram of the quantum dimer model on the square lattice. We use the sweeping cluster algorithm, which allows to implement exactly the dimer constraint, supplemented with a equal-time directed loop move that allows to sample winding sectors. We find a high-temperature critical phase with power-law correlations that extend down to the Rokshar-Kivelson point, in the vicinity of which a re-entrance effect in the lines of constant exponent is found. For small values of the kinetic energy strength, we find finite-temperature transitions to ordered states (columnar and staggered) which match those of interacting classical dimer models.

Multicritical point with infinite fractal symmetries

Recently, ``Pascal's triangle model'' constructed with $\mathrm{U}(1)$ rotor degrees of freedom was introduced and it was shown that ($\mathit{i}$) this model possesses an infinite series of fractal symmetries and ($\mathit{ii}$) it is the parent model of a series of ${Z}_{p}$ fractal models each with its own distinct fractal symmetry. In this work we discuss a multicritical point of Pascal's triangle model that is analogous to the Rokhsar-Kivelson point of the better known quantum dimer model. We demonstrate that the expectation value of the characteristic operator of each fractal symmetry at this multicritical point decays as a power law of space and this multicritical point is shared by the family of descendent ${Z}_{p}$ fractal models. Afterwards, we generalize our discussion to a $(3+1)d$ model termed ``Pascal's tetrahedron model'' that has both planar and fractal subsystem symmetries. We also establish a connection between Pascal's tetrahedron model and the $\mathrm{U}(1)$ Haah's code.

Quantum-critical electrodynamics of Luttinger fermions

We study the quantum electrodynamics of Luttinger fermions with quadratic band-crossing dispersion in three dimensions. The model can be viewed as the low-energy effective theory of a putative $U(1)$ quantum spin liquid with fermionic Luttinger spinons, or as an extension of the Luttinger-Abrikosov-Beneslavskii (LAB) model that accounts for transverse gauge fluctuations with finite photon velocity. Aided by a renormalization group analysis below four dimensions, we elucidate the presence and stability of quantum critical phenomena in this model. We find that the non-Fermi liquid LAB phase is stable against gauge fluctuations, and can thus also be viewed as a $U(1)$ spin liquid with gapless Luttinger spinons. We discover a multicritical point with Lifshitz scaling that corresponds to a time-reversal symmetry-breaking quantum phase transition from the LAB state to a chiral spin liquid with spinon Landau levels and birefringent emergent photons. This multicritical point is characterized by a finite fermion-photon coupling in the infrared and can be viewed as a fermionic analog of the Rokhsar-Kivelson point in three-dimensional quantum dimer models.

Quantum dimer model containing Rokhsar-Kivelson point expressed by spin- 12 Heisenberg antiferromagnets

We obtain a quantum dimer model (QDM) containing a Rokhsar-Kivelson (RK) point expressed by spin-$\frac{1}{2}$ Heisenberg antiferromagnets on a diamond-like decorated square lattice. This lattice has macroscopically degenerated nonmagnetic ground states, which are equivalent to the Hilbert space of a square-lattice QDM. Then, a square-lattice QDM containing the RK point as a second-order effective Hamiltonian is obtained by introducing further neighbor couplings as perturbation. Our model can provide a new method for the experimental realization of a resonating valence bond state in the QDM.

Deconfining Disordered Phase in Two-Dimensional Quantum Link Models.

We explore the ground-state physics of two-dimensional spin-1/2 U(1) quantum link models, one of the simplest nontrivial lattice gauge theories with fermionic matter within experimental reach for quantum simulations. Whereas in the large mass limit we observe Neél-like vortex-antivortex and striped crystalline phases, for small masses there is a transition from the striped phases into a disordered phase whose properties resemble those at the Rokhsar-Kivelson point of the quantum dimer model. This phase is characterized on ladders by boundary Haldane-like properties, such as vanishing parity and finite string ordering. Moreover, from studies of the string tension between gauge charges, we find that, whereas the striped phases are confined, the novel disordered phase present clear indications of being deconfined. Our results open exciting perspectives of studying highly nontrivial physics in quantum simulators, such as spin-liquid behavior and confinement-deconfinement transitions, without the need of explicitly engineering plaquette terms.

Critical colored-RVB states in the frustrated quantum Heisenberg model on the square lattice

We consider a family of SU(2)-symmetric Projected Entangled Paired States (PEPS) on the square lattice, defining colored-Resonating Valence Bond (RVB) states, to describe the quantum disordered phase of the J_1-J_2J1−J2 frustrated Heisenberg model. For J_2/J_1\sim 0.55J2/J1∼0.55 we show the emergence of critical (algebraic) dimer-dimer correlations – typical of Rokhsar-Kivelson (RK) points of quantum dimer models on bipartite lattices – while, simultaneously, the spin-spin correlation length remains short. Our findings are consistent with a spin liquid or a weak Valence Bond Crystal in the neighborhood of an RK point.

Lieb-Schultz-Mattis type theorem with higher-form symmetry and the quantum dimer models

The Lieb-Schultz-Mattis theorem dictates that a trivial symmetric insulator in lattice models is prohibited if lattice translation symmetry and $U(1)$ charge conservation are both preserved. In this paper, we generalize the Lieb-Schultz-Mattis theorem to systems with higher-form symmetries, which act on extended objects of dimension $n > 0$. The prototypical lattice system with higher-form symmetry is the pure abelian lattice gauge theory whose action consists only of the field strength. We first construct the higher-form generalization of the Lieb-Schultz-Mattis theorem with a proof. We then apply it to the $U(1)$ lattice gauge theory description of the quantum dimer model on bipartite lattices. Finally, using the continuum field theory description in the vicinity of the Rokhsar-Kivelson point of the quantum dimer model, we diagnose and compute the mixed 't Hooft anomaly corresponding to the higher-form Lieb-Schultz-Mattis theorem.

Novel constructive method for the quantum dimer model in spin-1/2 Heisenberg antiferromagnets with frustration on a diamond-like-decorated square lattice

We study spin-1/2 Heisenberg antiferromagnets on a diamond-like-decorated square lattice. The diamond-like-decorated square lattice is a lattice in which the bonds in a square lattice are replaced with diamond units. The diamond unit has two types of antiferromagnetic exchange interactions, and the ratio λ of the diagonal bond strength to that of the other four edges controls the frustration strength. For 0.974 < λ < 2, the present system has a nontrivial macroscopic degeneracy, which is called the macroscopically degenerated tetramer-dimer (MDTD) states. The MDTD states are identical to the Hilbert space of the Rokhsar–Kivelson (RK) quantum dimer model (QDM). By introducing further neighbor couplings in the MDTD states, we calculate the second-order effective Hamiltonian, which is exactly the same as the square-lattice QDM with a finite hopping amplitude t and dimer-dimer interaction v. Furthermore, we calculate v/|t| as a function of the ratio λ in the Heisenberg model and examine which phases of the square-lattice QDM appear in our obtained states. Our obtained QDM has a region where λ exhibits a finite hopping amplitude (|t| > 0) and repulsive interaction between dimers (v > 0). This suggests the possibility of realizing the resonating valence bond (RVB) state because the RVB state is obtained at v = |t|, which is known as the RK point.

Phases of quantum dimers from ensembles of classical stochastic trajectories

We study the connection between the phase behaviour of quantum dimers and the dynamics of classical stochastic dimers. At the so-called Rokhsar-Kivelson (RK) point a quantum dimer Hamiltonian is equivalent to the Markov generator of the dynamics of classical dimers. A less well understood fact is that away from the RK point the quantum-classical connection persists: in this case the Hamiltonian corresponds to a non-stochastic "tilted" operator that encodes the statistics of time-integrated observables of the classical stochastic problem. This implies a direct relation between the phase behaviour of quantum dimers and properties of ensembles of stochastic trajectories of classical dimers. We make these ideas concrete by studying fully packed dimers on the square lattice. Using transition path sampling - supplemented by trajectory umbrella sampling - we obtain the large deviation statistics of dynamical activity in the classical problem, and show the correspondence between the phase behaviour of the classical and quantum systems. The transition at the RK point between quantum phases of distinct order corresponds, in the classical case, to a trajectory phase transition between active and inactive dynamical phases. Furthermore, from the structure of stochastic trajectories in the active dynamical phase we infer that the ground state of quantum dimers has columnar order to one side of the RK point. We discuss how these results relate to those from quantum Monte Carlo, and how our approach may generalise to other problems.

Quantum Slow Relaxation and Metastability due to Dynamical Constraints.

One of the general mechanisms that give rise to the slow cooperative relaxation characteristic of classical glasses is the presence of kinetic constraints in the dynamics. Here we show that dynamical constraints can similarly lead to slow thermalization and metastability in translationally invariant quantum many-body systems. We illustrate this general idea by considering two simple models: (i)a one-dimensional quantum analogue to classical constrained lattice gases where excitation hopping is constrained by the state of neighboring sites, mimicking excluded-volume interactions of dense fluids; and (ii)fully packed quantum dimers on the square lattice. Both models have a Rokhsar-Kivelson (RK) point at which kinetic and potential energy constants are equal. To one side of the RK point, where kinetic energy dominates, thermalization is fast. To the other, where potential energy dominates, thermalization is slow, memory of initial conditions persists for long times, and separation of timescales leads to pronounced metastability before eventual thermalization. Furthermore, in analogy with what occurs in the relaxation of classical glasses, the slow-thermalization regime displays dynamical heterogeneity as manifested by spatially segregated growth of entanglement.

Finite-volume energy spectrum, fractionalized strings, and low-energy effective field theory for the quantum dimer model on the square lattice

We present detailed analytic calculations of finite-volume energy spectra, mean field theory, as well as a systematic low-energy effective field theory for the square lattice quantum dimer model. The analytic considerations explain why a string connecting two external static charges in the confining columnar phase fractionalizes into eight distinct strands with electric flux $\frac{1}{4}$. An emergent approximate spontaneously broken $SO(2)$ symmetry gives rise to a pseudo-Goldstone boson. Remarkably, this soft phonon-like excitation, which is massless at the Rokhsar-Kivelson (RK) point, exists far beyond this point. The Goldstone physics is captured by a systematic low-energy effective field theory. We determine its low-energy parameters by matching the analytic effective field theory with exact diagonalization results and Monte Carlo data. This confirms that the model exists in the columnar (and not in a plaquette or mixed) phase all the way to the RK point.

Phase Diagram of an Extended Quantum Dimer Model on the Hexagonal Lattice.

We introduce a quantum dimer model on the hexagonal lattice that, in addition to the standard three-dimer kinetic and potential terms, includes a competing potential part counting dimer-free hexagons. The zero-temperature phase diagram is studied by means of quantum Monte Carlo simulations, supplemented by variational arguments. It reveals some new crystalline phases and a cascade of transitions with rapidly changing flux (tilt in the height language). We analyze perturbatively the vicinity of the Rokhsar-Kivelson point, showing that this model has the microscopic ingredients needed for the "devil's staircase" scenario [Eduardo Fradkin et al. Phys. Rev. B 69, 224415 (2004)], and is therefore expected to produce fractal variations of the ground-state flux.

Double-semion topological order from exactly solvable quantum dimer models

We construct a generalized quantum dimer model on two-dimensional nonbipartite lattices including the triangular lattice, the star lattice and the kagome lattice. At the Rokhsar-Kivelson (RK) point, we obtain its exact ground states that are shown to be a fully gapped quantum spin liquid with the double-semion topological order. The ground-state wave function of such a model at the RK point is a superposition of dimer configurations with a nonlocal sign structure determined by counting the number of loops in the transition graph. We explicitly demonstrate the double-semion topological order in the ground states by showing the semionic statistics of monomer excitations. We also discuss possible implications of such double-semion resonating valence bond states to candidate quantum spin-liquid systems discovered experimentally and numerically in the past few years.

Transition of aZ3topologically ordered phase to trivial and critical phases

Topologically ordered quantum systems have robust physical properties, such as quasiparticle statistics and ground-state degeneracy, which do not depend on the microscopic details of the Hamiltonian. We consider topological phase transitions under a deformation such as an effective string tension on a $\mathbb{Z}_3$ topological state. This is studied in terms of the gauge-symmetry preserved quantum state renormalization group, first proposed by He, Moradi and Wen [Phys. Rev. B {\bf 90}, 205114 (2014)]. In this approach modular matrices $S$ and $T$ can be obtained and used as order parameters to characterize the topological properties of the phase and determine phase transitions. From a mapping to a classical 2D Potts model on the square lattice, the critical string tension, at which the transition to a topologically trivial phase takes place, can be obtained analytically and agrees with the numerically determined value. Such a transition can be generalized to a $\mathbb{Z}_N$ topological model under a string tension and determined in the same way. With different deformations, the $\mathbb{Z}_3$ topological phase can also be driven to a critical phase which contains, in the large deformation limit, a point analogous to Rokhsar-Kivelson point in the quantum dimer model.

Interfaces, strings, and a soft mode in the square lattice quantum dimer model

The quantum dimer model on the square lattice is equivalent to a $U(1)$ gauge theory. Quantum Monte Carlo calculations reveal that, for values of the Rokhsar-Kivelson (RK) coupling $\lambda < 1$, the theory exists in a confining columnar phase. The interfaces separating distinct columnar phases display plaquette order, which, however, is not realized as a bulk phase. Static "electric" charges are confined by flux tubes that consist of multiple strands, each carrying a fractionalized flux $\frac{1}{4}$. A soft pseudo-Goldstone mode emerges around $\lambda \approx 0$, long before one reaches the RK point at $\lambda = 1$.

The Rényi entanglement entropy of a general quantum dimer model at the RK point: a highly efficient algorithm

A highly efficient and simple-to-implement Monte Carlo algorithm is proposed for the evaluation of the Rényi entanglement entropy (REE) of the quantum dimer model (QDM) at the Rokhsar–Kivelson (RK) point. It makes possible the evaluation of REE at the RK point to the thermodynamic limit for a general QDM. We apply the algorithm to a QDM defined on the triangular and the square lattice in two dimensions and the simple and the face centered cubic (fcc) lattice in three dimensions. We find the REE on all these lattices follows perfect linear scaling in the thermodynamic limit, apart from an even–odd oscillation in the case of the square lattice. We also evaluate the topological entanglement entropy (TEE) with both a subtraction and an extrapolation procedure. We find the QDMs on both the triangular and the fcc lattice exhibit robust Z2 topological order. The expected TEE of ln2 is clearly demonstrated in both cases. Our large scale simulation also proves the recently proposed extrapolation procedure in cylindrical geometry to be a highly reliable way to extract the TEE of a topologically ordered system.

Geometric properties of loop condensed phases on the square lattice

Loop condensed phases are scale-invariant quantum liquid phases of matter. These phases include topologically ordered liquid phases such as the toric code as well as critical liquids such as the Rokhsar-Kivelson point of the quantum dimer model on the square lattice. To investigate the extent to which nonlocal geometric observables capture a signature of the nonlocal quantum order present in these phases, we compute geometric properties of such loop condensed states using directed loop Monte Carlo calculations. In particular, we investigate the loop condensed nature of ground state of the square lattice quantum dimer model at the Rokhsar-Kivelson point and compare with other loop condensed states on the square lattice, including those of the toric code and fully packed loop model. The common features of such liquids are a scale invariant distribution of loops and a fractal dimensionality of spanning loops. We find that the fractal dimension of the loop condensate of the square lattice quantum dimer model at the Rokhsar-Kivelson point is $3/2$, which provides quantitative confirmation of the effective height model that is commonly used to describe this critical dimer liquid.

Fermionic quantum dimer and fully-packed loop models on the square lattice

We consider fermionic fully-packed loop and quantum dimer models which serve as effective low-energy models for strongly correlated fermions on a checkerboard lattice at half and quarter filling, respectively. We identify a large number of fluctuationless states specific to each case, due to the fermionic statistics. We discuss the symmetries and conserved quantities of the system and show that for a class of fluctuating states in the half-filling case, the fermionic sign problem can be gauged away. This claim is supported by numerical evaluation of the low-lying states and can be understood by means of an algebraic construction. The elimination of the sign problem then allows us to analyze excitations at the Rokhsar-Kivelson point of the models using the relation to the height model and its excitations, within the single-mode approximation. We then discuss a mapping to a U(1) lattice gauge theory which relates the considered low-energy model to the compact quantum electrodynamics in 2+1 dimensions. Furthermore, we point out consequences and open questions in the light of these results.