We study the Rényi entanglement entropy (EE) of the two-dimensional J-Q model, the emblematic quantum spin model of deconfined criticality at the phase transition between antiferromagnetic and valence-bond-solid ground states. State-of-the-art quantum Monte Carlo calculations of the EE reveal critical corner contributions that scale logarithmically with the system size, with a coefficient in remarkable agreement with the form expected from a large-N conformal field theory with SO(N=5) symmetry. However, details of the bipartition of the lattice are crucial in order to observe this behavior. If the subsystem for the reduced density matrix does not properly accommodate valence-bond fluctuations, logarithmic contributions appear even for cornerless bipartitions. We here use a 45° tilted cut on the square lattice. Beyond supporting an SO(5) deconfined quantum critical point, our results for both the regular and tilted cuts demonstrate important microscopic aspects of the EE that are not captured by conformal field theory.
At a quantum critical point, the universal scaling behavior of Rényi entanglement entropy is controlled by the universality class of the codimension-two Rényi (or conical) defects in the infrared theory. In this work we perform a systematic study of critical correlations along Rényi defect lines in (2+1)
D quantum spin models realizing quantum phase transitions described by the O(3) Wilson-Fisher universality class, using large-scale quantum Monte Carlo simulations. We present numerical evidence that, for a fixed Rényi index
n, there exist multiple Rényi defect universality classes, with distinct critical exponents for the O(3) order parameter on the defect. These universality classes are realized by choosing microscopically different entanglement cuts in lattice models, which we classify as ordinary, special and extraordinary according to their relation to surface criticality. For the extraordinary entanglement cut, we further find evidence for a phase transition on the defect as a function of the Rényi index. Our results highlight the key role of defect universality classes in determining the universal scaling of Rényi entropy, and provide a framework for understanding the previously observed dependence of Rényi entropy scaling on microscopic lattice details.
The deconfined quantum critical point (DQCP) in square lattice S=1/2 quantum antiferromagnets has been extensively studied with a large body of evidence pointing to a weakly first-order transition scenario. Recent studies, which focused on entanglement at this nearly continuous DQCP in square lattice J-Q models, have observed conflicting bipartite entanglement entropy (EE) scaling behavior. One bipartition choice gave scaling coefficients in remarkable agreement with predictions from the unitary CFT corresponding to the putative DQCP. While another equally natural choice gave scaling coefficients in complete violation of unitary CFT that may be attributed to lack of scale invariance at the known weakly first-order behavior of the model. This motivates the exploration of DQCP behavior via entanglement measures in lattice models with distinct crystalline symmetries. Here we study a S=1/2 honeycomb model that hosts a nearly continuous transition between Néel and valence-bond-solid ground states relevant to probing DQCP. Using large-scale quantum Monte Carlo simulations, we compute the Rényi EE for a variety of bipartitions and test the CFT based description of the DQCP on the honeycomb lattice. For smooth bipartitions, we find no evidence of logarithmic corrections, in accordance with CFT, thereby essentially ruling out contributions from Goldstone modes. For subsystems with corners, CFT predicts universal logarithmic contributions, which we extract for corners with 60 and 120 degree angles and find close agreement with an emergent SO(5) CFT. While we observe scaling consistent with a critical system in the majority of cases, we also demonstrate an intriguing counterexample of the hexagon subsystem that exhibits a subtle period three oscillation. This results in three separate finite-size series, where the sign of the logarithmic term apparently changes depending on the series.
We study the entanglement properties of a class of ground states defined by matrix product states, which are generalizations of the valence bond solid (VBS) state in one dimension. It is shown that the transfer matrix of these states can be related to representations of the Temperley–Lieb algebra, allowing an exact computation of Renyi entropy. For an alternating bipartition, we find that the Renyi entropy can be mapped to an eight vertex model partition function on a rotated lattice. We also show that for the VBS state, the Renyi entropy of the alternating partition is described by a critical field theory with central charge c = 1. The generalization to SU ( n ) VBS and its connection with a dimerization transition in the entanglement Hamiltonian is discussed.
We perform a quantum Monte Carlo simulation of the resonating valence bond wavefunction on a two-dimensional square lattice with periodic boundary conditions. Using two replicas of the system, we calculate the second Renyi entropy on a spatial bipartition with a square geometry. Through a finite-size scaling analysis, we extract the logarithmic correction to the area law due to the presence of the sharp corners in the entangling surface. We find that the coefficient of this logarithm is positive with a value of 0.073 for a single
90∘ corner.
We use quantum Monte Carlo methods to study the ground-state phase diagram of a
S=1/2 honeycomb lattice magnet in which a nearest-neighbor antiferromagnetic exchange
J (favoring N\'eel order) competes with two different multispin interaction terms: a six-spin interaction
Q3 that favors columnar valence-bond solid (VBS) order, and a four-spin interaction
Q2 that favors staggered VBS order. For ${Q}_{3}\ensuremath{\sim}{Q}_{2}\ensuremath{\gg}J$, we establish that the competition between the two different VBS orders stabilizes N\'eel order in a large swath of the phase diagram even when
J is the smallest energy scale in the Hamiltonian. When ${Q}_{3}\ensuremath{\gg}({Q}_{2},J)$ $\mathrm{[}{Q}_{2}\ensuremath{\gg}({Q}_{3},J)\text{]}$, this model exhibits at zero temperature phase transition from the N\'eel state to a columnar (staggered) VBS state. We establish that the N\'eel-columnar VBS transition is continuous for all values of
Q2, and that critical properties along the entire phase boundary are well characterized by critical exponents and amplitudes of the noncompact
CP1 (NCCP1) theory of deconfined criticality, similar to what is observed on a square lattice. However, a surprising threefold anisotropy of the phase of the VBS order parameter at criticality, whose presence was recently noted at the
Q2=0 deconfined critical point, is seen to persist all along this phase boundary. We use a classical analogy to explore this by studying the critical point of a three-dimensional
XY model with a fourfold anisotropy field which is known to be weakly irrelevant at the three-dimensional
XY critical point. In this case, we again find that the critical anisotropy appears to saturate to a nonzero value over the range of sizes accessible to our simulations.