Resonating Valence Bond State Research Articles

Twisted graphene superlattices: resonating valence bond states and magnetic properties

Resonating valence bond (RVB) states are fundamental for understanding quantum spin liquids in two-dimensional (2D) systems. The RVB state is a collective phenomenon in which spins are uncoupled. 2D lattices such as triangular, honeycomb, and dice lattices were investigated using the Hubbard model and exact diagonalization method. We analyzed the total spin, spin–spin correlation functions, local magnetic moments, and spin and charge gaps as a function of on-site Coulomb repulsion, electron concentration, and electronic hopping parameters. Phase diagrams showed that RVB states can live in half-filled and hole-doped anisotropic triangular lattices. We found two types of RVB states: one in the honeycomb sublattice and the other in the centered hexagons in the triangular lattices. Owing to the novel discovery of exotic magnetic ordering in triangular moiré patterns in twisted bilayer graphene and transition metal dichalcogenide systems, our results provide physical insights into the onset of magnetism and possible spin liquid states in these layered materials.

Anyon condensation and confinement transition in a Kitaev spin liquid bilayer

Transitions between quantum spin liquids (QSLs) are fundamental problems lying beyond the Landau paradigm and requiring a deep understanding of the entanglement structures of QSLs called topological orders. The novel concept of anyon condensation has been proposed as a theoretical mechanism, predicting various possible transitions between topological orders, but it has long been elusive to confirm the mechanism in quantum spin systems. Here, we introduce a concrete spin model that incarnates the mechanism of the anyon condensation transition. Our model harbors two topological QSLs in different parameter regions, a non-Abelian Kitaev spin liquid (KSL) bilayer state and a resonating valence bond (RVB) state. The bilayer-KSL-to-RVB transition indeed occurs by the mechanism of anyon condensation, which we identify by using parton theories and exact diagonalization studies. Moreover, we observe “anyon confinement” phenomena in our numerical results, akin to the quark confinement in high-energy physics. Namely, non-Abelian Ising anyons of the bilayer KSL are confined in the transition to the RVB state. Implications and extensions of this study are discussed in various aspects such as (i) anyon-condensed multilayer construction of Kitaev's 16-fold way of anyon theories, (ii) an additional vison condensation transition from the RVB to a valence bond solid in the Kitaev bilayer system, (iii) dynamical anyon condensation in a non-Hermitian Kitaev bilayer, (iv) generalizations of our model to other lattice geometries, and (v) experimental realizations. This work puts together the two fascinating QSLs that are extensively studied in modern condensed matter and quantum physics into a concrete spin model, offering a comprehensive picture that unifies the anyon physics of the Kitaev spin liquids and the resonating valence bonds. Published by the American Physical Society 2024

Thermodynamics of resonating-valence-bond states toward the understanding of quantum spin liquid phenomena.

A quantum spin liquid (QSL) is a state of matter in which spins do not exhibit magnetic order. In contrast to paramagnets, the spins in a QSL interact strongly, similar to conventional ordered magnets; however the thermodynamic stability of QSLs is rarely studied. Here, the thermodynamic properties of centered hexagon nanoclusters were investigated using the Hubbard model, which was solved using exact numerical diagonalization. The total spin, spin-spin correlation functions, local magnetic moments, charge and spin gaps, and magnetocaloric effect were analyzed for a half-filled band as a function of the ratio between the on-site Coulomb repulsion and electronic hopping (U/t). The centered hexagon nanocluster exhibited an antiferromagnetic (AFM) behavior with exotic magnetic ordering. Resonating-valence-bond (RVB) states were observed for intermediate values of U/t, in which short-range spin-spin correlation functions were suppressed to minimize spin frustration. The AFM order was examined in terms of the Néel-like temperature derived from the temperature dependence of the magnetic susceptibility. An interesting result is that the systems under external magnetic fields exhibited an inverse magnetocaloric effect, which was remarkable for intermediate values of U/t, where the RVB state was observed. Owing to the novel discovery of exotic magnetic ordering in triangular moiré patterns in twisted bilayer graphene (TBLG) systems, these results provide insights into the onset of magnetism and the possible spin liquid states in these graphene moiré materials.

Physics-Inspired Quantum Simulation of Resonating Valence Bond States─A Prototypical Template for a Spin-Liquid Ground State.

Spin liquids─an emergent, exotic collective phase of matter─have garnered enormous attention in recent years. While experimentally many prospective candidates have been proposed and realized, theoretically modeling real materials that display such behavior may pose serious challenges due to the inherently high correlation content of such phases. Over the last few decades, the second-quantum revolution has been the harbinger of a novel computational paradigm capable of initiating a foundational evolution in computational physics. In this report, we strive to use the power of the latter to study a prototypical model, a spin-1/2-unit cell of a Kagome antiferromagnet. Extended lattices of such unit cells are known to possess a magnetically disordered spin-liquid ground state. We employ robust classical numerical techniques such as the density-matrix renormalization group (DMRG) to identify the nature of the ground state through a matrix-product state (MPS) formulation. We subsequently use the gained insight to construct an auxiliary Hamiltonian with reduced measurables and also design an ansatz that is modular and gate-efficient. With robust error-mitigation strategies, we are able to demonstrate that the said ansatz is capable of accurately representing the target ground state even on a real IBMQ backend within 1% accuracy in energy. Since the protocol is linearly scaling O(n) in the number of unit cells, gate requirements, and the number of measurements, it is straightforwardly extendable to larger Kagome lattices that can pave the way for efficient construction of spin-liquid ground states on a quantum device.

Separability and entanglement of resonating valence-bond states

We investigate separability and entanglement of Rokhsar-Kivelson (RK) states and resonating valence-bond (RVB) states. These states play a prominent role in condensed matter physics, as they can describe quantum spin liquids and quantum critical states of matter, depending on their underlying lattices. For dimer RK states on arbitrary tileable graphs, we prove the exact separability of the reduced density matrix of kk disconnected subsystems, implying the absence of bipartite and multipartite entanglement between the subsystems. For more general RK states with local constraints, we argue separability in the thermodynamic limit, and show that any local RK state has zero logarithmic negativity, even if the density matrix is not exactly separable. In the case of adjacent subsystems, we find an exact expression for the logarithmic negativity in terms of partition functions of the underlying statistical model. For RVB states, we show separability for disconnected subsystems up to exponentially small terms in the distance dd between the subsystems, and that the logarithmic negativity is exponentially suppressed with dd. We argue that separability does hold in the scaling limit, even for arbitrarily small ratio d/Ld/L, where LL is the characteristic size of the subsystems. Our results hold for arbitrary lattices, and encompass a large class of RK and RVB states, which include certain gapped quantum spin liquids and gapless quantum critical systems.

50 years of quantum spin liquids

In 1973, Philip Anderson published a paper introducing the resonating valence bond state, which can be recognized in retrospect as a topologically ordered phase of matter — one that cannot be classified in the conventional way according to its patterns of spontaneously broken symmetry. Steven Kivelson and Shivaji Sondhi reflect on the impact of this paper over the past 50 years.

Resonating Valence Bond States in an Electron-Phonon System.

We study a simple electron-phonon model on square and triangular versions of the Lieb lattice using an asymptotically exact strong coupling analysis. At zero temperature and electron density n=1 (one electron per unit cell), for various ranges of parameters in the model, we exploit a mapping to the quantum dimer model to establish the existence of a spin-liquid phase with Z_{2} topological order (on the triangular lattice) and a multicritical line corresponding to a quantum critical spin liquid (on the square lattice). In the remaining part of the phase diagram, we find a host of charge-density-wave phases (valence-bond solids), a conventional s-wave superconducting phase, and with the addition of a small Hubbard U to tip the balance, a phonon-induced d-wave superconducting phase. Under a special condition, we find a hidden pseudospin SU(2) symmetry that implies an exact constraint on the superconducting order parameters.

Dynamical Preparation of Quantum Spin Liquids in Rydberg Atom Arrays.

We theoretically analyze recent experiments [Semeghini etal., Science 374, 1242 (2021)SCIEAS0036-807510.1126/science.abi8794] demonstrating the onset of a topological spin liquid using a programmable quantum simulator based on Rydberg atom arrays. In the experiment, robust signatures of topological order emerge in out-of-equilibrium states that are prepared using a quasiadiabatic state preparation protocol. We show theoretically that the state preparation protocol can be optimized to target the fixed point of the topological phase-the resonating valence bond state of hard dimers-in a time that scales linearly with the number of atoms. Moreover, we provide a two-parameter variational manifold of tensor network states that accurately describe the many-body dynamics of the preparation process. Using this approach we analyze the nature of the nonequilibrium state, establishing the emergence of topological order.

Gutzwiller approximation approach to the SU(4) t−J model

We develop the Gutzwiller approximation method to obtain the renormalized Hamiltonian of the SU(4) $t\text{\ensuremath{-}}J$ model with the corresponding renormalization factors. Subsequently, a mean-field theory is employed on the renormalized Hamiltonian of the model on the honeycomb lattice under the scenario of a cooperative condensation of carriers moving in the resonating valence bond state of flavors. In particular, we find that the extended $s$-wave superconducting state is more favorable than the $d\ifmmode\pm\else\textpm\fi{}id$-wave superconducting state in the doping range close to quarter filling. The pairing states of the SU(4) case reveal the property that the spin-singlet pairing and the spin-triplet pairing can coexist simultaneously. Our results might provide new insights into the twisted bilayer graphene system.

Probing resonating valence bond states in artificial quantum magnets

Designing and characterizing the many-body behaviors of quantum materials represents a prominent challenge for understanding strongly correlated physics and quantum information processing. We constructed artificial quantum magnets on a surface by using spin-1/2 atoms in a scanning tunneling microscope (STM). These coupled spins feature strong quantum fluctuations due to antiferromagnetic exchange interactions between neighboring atoms. To characterize the resulting collective magnetic states and their energy levels, we performed electron spin resonance on individual atoms within each quantum magnet. This gives atomic-scale access to properties of the exotic quantum many-body states, such as a finite-size realization of a resonating valence bond state. The tunable atomic-scale magnetic field from the STM tip allows us to further characterize and engineer the quantum states. These results open a new avenue to designing and exploring quantum magnets at the atomic scale for applications in spintronics and quantum simulations.

Efficient variational contraction of two-dimensional tensor networks with a non-trivial unit cell

Tensor network states provide an efficient class of states that faithfully capture strongly correlated quantum models and systems in classical statistical mechanics. While tensor networks can now be seen as becoming standard tools in the description of such complex many-body systems, close to optimal variational principles based on such states are less obvious to come by. In this work, we generalize a recently proposed variational uniform matrix product state algorithm for capturing one-dimensional quantum lattices in the thermodynamic limit, to the study of regular two-dimensional tensor networks with a non-trivial unit cell. A key property of the algorithm is a computational effort that scales linearly rather than exponentially in the size of the unit cell. We demonstrate the performance of our approach on the computation of the classical partition functions of the antiferromagnetic Ising model and interacting dimers on the square lattice, as well as of a quantum doped resonating valence bond state. Tensor network states provide an efficient class of states that faithfully capture strongly correlated quantum models and systems in classical statistical mechanics. While tensor networks can now be seen as becoming standard tools in the description of such complex many-body systems, close to optimal variational principles based on such states are less obvious to come by. In this work, we generalize a recently proposed variational uniform matrix product state algorithm for capturing one-dimensional quantum lattices in the thermodynamic limit, to the study of regular two-dimensional tensor networks with a non-trivial unit cell. A key property of the algorithm is a computational effort that scales linearly rather than exponentially in the size of the unit cell. We demonstrate the performance of our approach on the computation of the classical partition functions of the antiferromagnetic Ising model and interacting dimers on the square lattice, as well as of a quantum doped resonating valence bond state.

Triplet Resonating Valence Bond State and Superconductivity in Hund's Metals.

A central idea in strongly correlated systems is that doping a Mott insulator leads to a superconductor by transforming the resonating valence bonds (RVBs) into spin-singlet Cooper pairs. Here, we argue that a spin-triplet RVB (tRVB) state, driven by spatially, or orbitally anisotropic ferromagnetic interactions can provide the parent state for triplet superconductivity. We apply this idea to the iron-based superconductors, arguing that strong on site Hund's interactions develop intra-atomic tRVBs between the t_{2g} orbitals. On doping, the presence of two iron atoms per unit cell allows these interorbital triplets to coherently delocalize onto the Fermi surface, forming a fully gapped triplet superconductor. This mechanism gives rise to a unique staggered structure of on site pair correlations, detectable as an alternating π phase shift in a scanning Josephson tunneling microscope.

Magnetism and topological phases in an interacting decorated honeycomb lattice with spin-orbit coupling

We study the interplay between spin-orbit coupling (SOC) and Coulomb repulsion in a Hubbard model on a decorated honeycomb lattice which leads to a plethora of phases. While a quantum spin hall insulator is stable at weak Coulomb repulsion and moderate SOC, a semimetallic phase emerges at large SOC in a broad range of Coulomb repulsion. This semimetallic phase has topological properties not observed in conventional metals such as a finite, non-quantized spin Hall conductivity. At large Coulomb repulsion and negligible spin-orbit coupling, electronic correlations stabilize a resonance valence bond (RVB) spin liquid state in contrast to the classical antiferromagnetic state predicted by mean-field theory. Under sufficiently strong SOC, such RVB state is transformed into a magnetic insulator consisting on S~3/2 localized moments on a honeycomb lattice with antiferromagnetic order and topological features.

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.

Gapped Z2 spin liquid in the breathing kagome Heisenberg antiferromagnet

We investigate the spin-1/2 Heisenberg antiferromagnet on the kagome lattice with breathing anisotropy (i.e. with weak and strong triangular units), constructing an improved simplex Resonating Valence Bond (RVB) ansatz by successive applications (up to three times) of local quantum gates which implement a filtering operation on the bare nearest-neighbor RVB state. The resulting Projected Entangled Pair State involves a small number of variational parameters (only one at each level of application) and preserves full lattice and spin-rotation symmetries. Despite its simple analytic form, the simplex RVB provides very good variational energies at strong and even intermediate breathing anisotropy. We show that it carries $Z_2$ topological order which does not fade away under the first few applications of the quantum gates, suggesting that the RVB topological spin liquid becomes a competing ground state candidate for the kagome antiferromagnet at large breathing anisotropy.

SU(4) topological resonating valence bond spin liquid on the square lattice

We generalize the construction of the spin-1/2 SU(2) Resonating Valence Bond (RVB) state to the case of the self-conjugate $\bf 6$-representation of SU(4). As for the case of SU(2) [J-Y. Chen and D. Poilblanc, Phys. Rev. B $\textbf{97}$, 161107(R) (2018)], we use the Projected Entangled Pair State (PEPS) formalism to derive a simple (two-dimensional) family of generalized SU(4) RVB states on the square lattice. We show that, when longer-range SU(4)-singlet bonds are included, a local gauge symmetry is broken down from U(1) to $\mathbb{Z}_2$, leading to the emergence of a short-range spin liquid. Evidence for the topological nature of this spin liquid is provided by the investigation of the Renyi entanglement entropy of infinitely-long cylinders and of the modular matrices. Relevance to microscopic models and experiments of ultracold atoms is discussed.

Solving the quantum dimer and six-vertex models one electric field line at a time

The nature and the very existence of the resonant plaquette valence bond state that separates the classical columnar phase and the Rokhsar and Kivelson point in the quantum dimer model remains unsettled. Here we take a different line of attack on this model, and on the closely related six vertex model, by exploiting the global conservation law of the number of electric field lines. This allows us to study a single fluctuating electric field line which we show maps exactly onto a one dimensional spin chain. In the case of the six vertex model, the electric field line maps onto the celebrated spin 1/2 XXZ model which can be solved exactly. In the quantum dimer model, the electric field line is mapped onto a two-leg spin 1/2 ladder, which we study using numerical exact diagonalization. Our findings are consistent with the existence of three distinct phases including a Luttinger liquid phase, the one-dimensional precursor to the two-dimensional plaquette valence bond solid. The uncanny resemblance of our quasi-one-dimensional electric field line problem to the full two-dimensional problem suggests that much of the behavior of the latter might be understood by thinking of it as a closely packed array of field lines which themselves are undergoing non-trivial phase transitions.

Noise-tolerant signature of ZN topological order in quantum many-body states

Topologically ordered states are fundamentally important in theoretical physics, which are also suggested as promising candidates to build fault-tolerant quantum devices. However, it is still elusive how topological orders can be affected or detected under noises. In this work, we find a quantity, termed as the ring degeneracy $\mathcal{D}$, which is robust under pure noise to detect both trivial and intrinsic topological orders. The ring degeneracy is defined as the degeneracy of the solutions of the self-consistent equations that encode the contraction of the corresponding tensor network(TN). For the $Z_N$ orders, we find that the ring degeneracy satisfies a simple relation $\mathcal{D} = (N + 1)/2 + d$, with $d = 0$ for odd $N$ and $d = 1/2$ for even $N$. Simulations on several non-trivial states (two-dimensional Ising model, $Z_N$ topological states, and resonating valence bond states) show that the ring degeneracy can tolerate noises up to a strength associated to the gap of the TN boundary theory.

Symmetry-Projected Jastrow Mean-Field Wave Function in Variational Monte Carlo.

We extend our low-scaling variational Monte Carlo (VMC) algorithm to optimize the symmetry-projected Jastrow mean-field (SJMF) wave functions. These wave functions consist of a symmetry-projected product of a Jastrow and a general broken-symmetry mean-field reference. Examples include Jastrow antisymmetrized geminal power, Jastrow-Pfaffians, and resonating valence bond states, among others, all of which can be treated with our algorithm. We will demonstrate using benchmark systems including the nitrogen molecule, a chain of hydrogen atoms, and 2-D Hubbard model that a significant amount of correlation can be obtained by optimizing the energy of the SJMF wave function. This can be achieved at a relatively small cost when multiple symmetries including spin, particle number, and complex conjugation are simultaneously broken and projected. We also show that reduced density matrices can be calculated using the optimized wave functions, which allows us to calculate other observables such as correlation functions and will enable us to embed the VMC algorithm in a complete active-space self-consistent field calculation.

Gapped SU(3) spin liquid with Z3 topological order

We construct a topological spin liquid (TSL) model on the kagome lattice, with SU(3) symmetry with the fundamental representation at each lattice site, based on Projected Entangled Pair States (PEPS). Using the PEPS framework, we can adiabatically connect the model to a fixed point model (analogous to the dimer model for Resonating Valence Bond states) which we prove to be locally equivalent to a $Z_3$ quantum double model. Numerical study of the interpolation reveals no sign of a phase transition or long-range order, characterizing the model conclusively as a gapped TSL. We further study the entanglement spectrum of the model and find that while it is gapped, it exhibits branches with vastly different velocities, with the slow branch matching the counting of a chiral $SU(3)_1$ CFT, suggesting that it can be deformed to a model with chiral $SU(3)_1$ entanglement spectrum.