Thanks to Pfaffian techniques, we study the Rényi entanglement entropies and the entanglement spectrum of large subsystems for two-dimensional Rokhsar–Kivelson wavefunctions constructed from a dimer model on the triangular lattice. By including a fugacity t on some suitable bonds, one interpolates between the triangular lattice ( t = 1) and the square lattice ( t = 0). The wavefunction is known to be a massive topological liquid for t > 0 whereas it is a gapless critical state at t = 0. We mainly consider two geometries for the subsystem: that of a semi-infinite cylinder and the disc-like set-up proposed by Kitaev and Preskill (2006 Phys. Rev. Lett. 96 110404). In the cylinder case, the entropies contain an extensive term—proportional to the length of the boundary—and a universal subleading constant s n ( t ). Fitting these cylinder data (up to a perimeter of L = 32 sites) provides s n with a very high numerical accuracy (10 − 9 at t = 1 and 10 − 6 at t = 0.5). In the topological liquid phase we find s n ( t > 0) = − ln2, independent of the fugacity t and the Rényi parameter n . At t = 0 we recover a previously known result, for n < 1 and s n ( t = 0) = − ln(2)/( n − 1) for n > 1. In the disc-like geometry—designed to get rid of the boundary contributions—we find an entropy s n KP ( t 0) = − ln2 in the whole massive phase whatever n > 0, in agreement with the result of Flammia et al (2009 Phys. Rev. Lett. 103 261601). Some results for the gapless limit are discussed.
We study the
n=2 R\'enyi entanglement entropy of the triangular quantum dimer model via Monte Carlo sampling of Rokhsar-Kivelson- (RK-) like ground-state wave functions. Using the construction proposed by Kitaev and Preskill [Phys. Rev. Lett. 96, 110404 (2006)] and an adaptation of the Monte Carlo algorithm described by Hastings et al. [Phys. Rev. Lett. 104, 157201 (2010)], we compute the topological entanglement entropy (TEE) at the RK point $\ensuremath{\gamma}=(1.001\ifmmode\pm\else\textpm\fi{}0.003)\mathrm{ln}2$, confirming earlier results. Additionally, we compute the TEE of the ground state of a generalized RK-like Hamiltonian and demonstrate the universality of TEE over a wide range of parameter values within a topologically ordered phase approaching a quantum phase transition. For system sizes that are accessible numerically, we find that the quantization of TEE depends sensitively on correlations. We characterize corner contributions to the entanglement entropy and show that these are well described by shifts proportional to the number and types of corners in the bipartition.
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.
An outstanding problem in the study of topologically ordered phases is to find methods to distinguish topologically ordered from unordered phases. The topological entanglement entropy, a negative shift to the entanglement entropy scaling, has emerged as one possible distinguishing probe. The n=2 Rényi topological entanglement entropy of the triangular quantum dimer model is studied by classical Monte Carlo sampling of Rokhsar-Kivelson-like ground state wave functions, and using quantum Monte Carlo sampling for Hamiltonian parameter values away from the Rokhsar-Kivelson point. Both results demonstrate the universality of topological entanglement entropy over a wide range of parameter values within a topologically ordered phase. For systems sizes that are accessible numerically, it is found that the quantization of topological entanglement entropy depends sensitively on correlations. At the Rokhsar-Kivelson point, competing corner contributions to the entanglement entropy are also characterized.Another class of topologically ordered materials existing in three dimensions are topological insulators, which exhibit Dirac modes on its surface. The anomalous currents arising from these modes are investigated. The Dirac energy eigenstates are computed on a flat torus (genus one topology) and closed cap cylinder (genus zero topology) for various mass-term geometries. From the resulting properties of the surface spectra, a potential application for a flux-charge qubit is presented.
We consider a family of generalized Rokhsar-Kivelson (RK) Hamiltonians, which are reverse-engineered to have an arbitrary edge-weighted superposition of dimer coverings as their exact ground state at the RK point. We focus on a quantum dimer model on the triangular lattice, with doubly periodic edge weights. For simplicity we consider a
2×1 periodic model in which all weights are set to one except for a tunable horizontal edge weight labeled
α. We analytically show that the model exhibits a continuous quantum phase transition at
α=3, changing from a topological
Z2 quantum spin liquid (
α<3) to a columnar ordered state (
α>3). The dimer-dimer correlator decays exponentially on both sides of the transition with the correlation length
ξ∝1/∣α−3∣ and as a power-law at criticality. The vison correlator exhibits an exponential decay in the spin liquid phase, but becomes a constant in the ordered phase, which we explain in terms of loop statistics of the double-dimer model. Using finite-size scaling of the vison correlator, we extract critical exponents consistent with the 2D Ising universality class. Additionally, we analytically show that the topological Rényi entropy of order
∞ (topological min-entropy) changes from
log2 for the quantum spin liquid phase
α<3, to
0 for the ordered phase
α>3, thereby analytically confirming the topological nature of the phase transition.
Entanglement is one of the most intriguing features of quantum theory and a main resource in quantum information science. Ground states of quantum many-body systems with local interactions typically obey an “area law” which means that the entanglement entropy is proportional to the boundary length. It is exceptional when the system is gapless, and the area law had been believed to be violated by at most a logarithm over two decades. Recent discovery of Motzkin and Fredkin spin chain models is striking, since these models provide significant violation of the entanglement beyond the belief, growing as a square root of the volume in spite of local interactions. In this paper, we first analytically compute the Rényi entropy of the Motzkin and Fredkin models by careful treatment of asymptotic analysis. The Rényi entropy is an important quantity, since the whole spectrum of an entangled subsystem is reconstructed once the Rényi entropy is known as a function of its parameter. We find nonanalytic behavior of the Rényi entropy with respect to the parameter, which is a novel phase transition never seen in any other spin chain studied so far. Interestingly, similar behavior is seen in the Rényi entropy of Rokhsar–Kivelson states in two dimensions.
We consider the two-dimensional Ising model on an infinitely long cylinder and study the probabilities
pi to observe a given spin configuration
i along a circular section of the cylinder. These probabilities also occur as eigenvalues of reduced density matrices in some Rokhsar-Kivelson wave functions. We analyze the subleading constant to the R\'enyi entropy ${R}_{n}=1/(1\ensuremath{-}n)\text{ln}({\ensuremath{\sum}}_{i}{p}_{i}^{n})$ and discuss its scaling properties at the critical point. Studying three different microscopic realizations, we provide numerical evidence that it is universal and behaves in a steplike fashion as a function of
n with a discontinuity at the Shannon point
n=1. As a consequence, a field theoretical argument based on the replica trick would fail to give the correct value at this point. We nevertheless compute it numerically with high precision. Two other values of the R\'enyi parameter are of special interest:
n=1/2 and $n=\ensuremath{\infty}$ are related in a simple way to the Affleck-Ludwig boundary entropies associated to free and fixed boundary conditions, respectively.