Gapless quantum spin liquid and global phase diagram of the spin-1/2 J1-J2 square antiferromagnetic Heisenberg model.

The nature of the zero-temperature phase diagram of the spin-1/2J1-J2 Heisenberg model on a square lattice has been debated in the past three decades, and it remains one of the fundamental problems unsettled in the study of quantum many-body theory. By using the state-of-the-art tensor network method, specifically, the finite projected entangled pair state (PEPS) algorithm, to simulate the global phase diagram of the J1-J2 Heisenberg model up to 24×24 sites, we provide very solid evidences to show that the nature of the intermediate nonmagnetic phase is a gapless quantum spin liquid (QSL), whose spin-spin and dimer-dimer correlations both decay with a power law behavior. There also exists a valence-bond solid (VBS) phase in a very narrow region 0.56≲J2/J1≤0.61 before the system enters the well known collinear antiferromagnetic phase. We stress that we make the first detailed comparison between the results of PEPS and the well-established density matrix renormalization group (DMRG) method through one-to-one direct benchmark for small system sizes, and thus give rise to a very solid PEPS calculation beyond DMRG. Our numerical evidences explicitly demonstrate the huge power of PEPS for highly frustrated spin systems. Finally, an effective field theory is also proposed to understand the physical nature of the discovered gapless QSL and its relation to deconfined quantum critical point (DQCP).

Simulating Chiral Spin Liquids with Projected Entangled-Pair States.

Doubts have been raised on the representation of chiral spin liquids exhibiting topological order in terms of projected entangled pair states (PEPSs). Here, starting from a simple spin-1/2 chiral frustrated Heisenberg model, we show that a faithful representation of the chiral spin liquid phase is in fact possible in terms of a generic PEPS upon variational optimization. We find a perfectly chiral gapless edge mode and a rapid decay of correlation functions at short distances consistent with a bulk gap, concomitant with a gossamer long-range tail originating from a PEPS bulk-edge correspondence. For increasing bond dimension, (i) the rapid decrease of spurious features-SU(2) symmetry breaking and long-range tails in correlations-together with (ii) a faster convergence of the ground state energy as compared to state-of-the-art cylinder matrix-product state simulations involving far more variational parameters, prove the fundamental relevance of the PEPS ansatz for simulating systems with chiral topological order.

SU(3)_{1} Chiral Spin Liquid on the Square Lattice: A View from Symmetric Projected Entangled Pair States.

Quantum spin liquids can be faithfully represented and efficiently characterized within the framework of projected entangled pair states (PEPS). Guided by extensive exact diagonalization and density matrix renormalization group calculations, we construct an optimized symmetric PEPS for a SU(3)_{1} chiral spin liquid on the square lattice. Characteristic features are revealed by the entanglement spectrum (ES) on an infinitely long cylinder. In all three Z_{3} sectors, the level counting of the linear dispersing modes is in full agreement with SU(3)_{1} Wess-Zumino-Witten conformal field theory prediction. Special features in the ES are shown to be in correspondence with bulk anyonic correlations, indicating a fine structure in the holographic bulk-edge correspondence. Possible universal properties of topological SU(N)_{k} chiral PEPS are discussed.

Topological order in the projected entangled-pair states formalism: transfer operator and boundary Hamiltonians.

We study the structure of topological phases and their boundaries in the projected entangled-pair states (PEPS) formalism. We show how topological order in a system can be identified from the structure of the PEPS transfer operator and subsequently use these findings to analyze the structure of the boundary Hamiltonian, acting on the bond variables, which reflects the entanglement properties of the system. We find that in a topological phase, the boundary Hamiltonian consists of two parts: A universal nonlocal part which encodes the nature of the topological phase and a nonuniversal part which is local and inherits the symmetries of the topological model, which helps to infer the structure of the boundary Hamiltonian and thus possibly of the physical edge modes.

Constructing a gapless spin-liquid state for the spin-1/2 J(1)-J(2) Heisenberg model on a square lattice.

We construct a class of projected entangled pair states which is exactly the resonating valence bond wave functions endowed with both short range and long range valence bonds. With an energetically preferred resonating valence bond pattern, the wave function is simplified to live in a one-parameter variational space. We tune this variational parameter to minimize the energy for the frustrated spin-1/2 J(1)-J(2) antiferromagnetic Heisenberg model on the square lattice. Taking a cylindrical geometry, we are able to construct four topological sectors with an even or odd number of fluxes penetrating the cylinder and an even or odd number of spinons on the boundary. The energy splitting in different topological sectors is exponentially small with the cylinder perimeter. We find a power law decay of the dimer correlation function on a torus, and a lnL correction to the entanglement entropy, indicating a gapless spin-liquid phase at the optimum parameter.

Fractionalization of itinerant anyons in one-dimensional chains.

We construct models of interacting itinerant non-Abelian anyons moving along one-dimensional chains, focusing, in particular, on itinerant Ising anyon chains, and derive effective anyonic t-J models for the low-energy sectors. Solving these models by exact diagonalization, we find a fractionalization of the anyons into charge and (non-Abelian) anyonic degrees of freedom--a generalization of spin-charge separation of electrons which occurs in Luttinger liquids. A detailed description of the excitation spectrum by combining spectra for charge and anyonic sectors requires a subtle coupling between charge and anyonic excitations at the microscopic level (which we also find to be present in Luttinger liquids), despite the macroscopic fractionalization.

Entanglement spectra of quantum Heisenberg ladders.

Bipartite entanglement measures are surprisingly useful tools to investigate quantum phases of correlated electrons. Here, I analyze the entanglement spectrum of gapped two-leg quantum Heisenberg ladders on a periodic ribbon partitioned into two identical periodic chains. The entanglement spectrum closely reflects the low-energy gapless spectrum of each individual edge. This extends the conjecture initially drawn for fractional quantum Hall systems to the field of quantum magnetism, stating a direct correspondence between the low-energy entanglement spectrum of a partitioned system and the true spectrum of the virtual edges. A mapping of the reduced density matrix to a thermodynamic density matrix is also proposed via the introduction of an effective temperature.

Quantum melting of valence-bond crystal insulators and novel supersolid phase at commensurate density.

Bosonic and fermionic Hubbard models on the checkerboard lattice are studied numerically for infinite on-site repulsion. At particle density n=1/4 and strong nearest-neighbor repulsion, insulating Valence-Bond crystals (VBC) of resonating particle pairs are stabilized. Their melting into superfluid or metallic phases under increasing hopping is investigated at T=0 K. We identify a novel and unconventional commensurate VBC supersolid region, precursor to the melting of the bosonic crystal. Hardcore bosons (spins) are compared to fermions (electrons), as well as positive to negative (frustrating) hoppings.

Supersolid phases of a doped valence-bond quantum antiferromagnet: evidence for a coexisting superconducting order parameter.

Motivated by numerical evidence of the valence-bond ground state of the two-dimensional Heisenberg pyrochlore lattice, we argue using a t-J model that it evolves under doping into novel phases characterized by superconductivity coexisting with the underlying valence-bond solid order. A fermionic mean-field theory supplemented by exact diagonalization results provide strong arguments in favor of the stability of such supersolid phases. The resemblance with modulated superconducting patterns in high-Tc cuprates as well as possible relevance to frustrated noncuprate superconductors such as spinels and pyrochlores is discussed.

Magnetic field induced transition in a quantum magnet described by the quantum dimer model.

The effect of a magnetic field on a gapped quantum magnet is described within the framework of the quantum dimer model. A minimal model describing the proliferation of itinerant spinons above a critical field is proposed and investigated by Lanczos exact diagonalizations and quantum Monte Carlo simulations. For both square and triangular lattices, it is shown that spinons are fully polarized and Bose condense. This offers a novel scenario of a quantum critical point in the dimer-liquid phase (triangular lattice) characterized by the continuous appearance of a spinon superfluid density, contrasting with the usual triplet condensation picture. The possible role of other spinon kinetic terms neglected in the model are discussed.

Properties of holons in the quantum dimer model.

I introduce a doped two-dimensional quantum dimer model describing a doped Mott insulator and retaining the original Fermi statistics of the electrons. This model shows a rich phase diagram including a d-wave hole-pair unconventional superconductor at small enough doping and a bosonic superfluid at large doping. The hole kinetic energy is shown to favor binding of topological defects to the bare fermionic holons turning them into bosons, in agreement with arguments based on resonating valence bond wave function. Results are discussed in the context of cuprate superconductors.

Phase separation and flux quantization in the doped quantum dimer model on square and triangular lattices.

The doped two-dimensional quantum dimer model is investigated by numerical techniques on the square and triangular lattices, with significantly different results. On the square lattice, at small enough doping, there is always a phase separation between an insulating valence-bond solid and a uniform superfluid phase, whereas on the triangular lattice, doping leads directly to a uniform superfluid in a large portion of the resonating-valence-bond (RVB) phase. Under an applied Aharonov-Bohm flux, the superfluid exhibits quantization in terms of half-flux quanta, consistent with Q=2e elementary charge quanta in transport properties.

Enhanced pairing in doped quantum magnets with frustrated hole motion.

Evidence for strong pairing at arbitrarily small J/t is provided in a t-J model on the checkerboard lattice for a specific sign of the hopping amplitude. Destructive quantum interferences suppress Nagaoka ferromagnetism when J/t-->0 and drastically reduce coherent hole motion in the fluctuating singlet background. It is shown that, by pairing in various orbital symmetry channels, holes can benefit from a large gain of kinetic energy.

Spin-charge separation in two-dimensional frustrated quantum magnets.

The dynamics of a mobile hole in two-dimensional frustrated quantum magnets is investigated by exact diagonalization techniques. Our results provide evidence for spin-charge separation upon doping the kagome lattice, a prototype of a spin liquid. In contrast, in the checkerboard lattice, a symmetry broken valence bond crystal, a small quasiparticle peak is seen for some crystal momenta, a finding interpreted as a restoration of weak holon-spinon confinement.

Superconducting gap for a two-leg t-J ladder.

Single-particle diagonal and off-diagonal Green's functions of a two-leg t-J ladder at 1/8 doping are investigated by exact diagonalizations techniques. A numerically tractable expression for the superconducting gap is proposed and the frequency dependence of the real and imaginary parts of the gap are determined. The role of the low-energy gapped spin modes, whose energies are computed by a (one-step) contractor renormalization procedure, is discussed.

Tetramerization of a frustrated spin-1/2 chain.

We investigate a model of a frustrated spin-1/2 Heisenberg chain coupled to adiabatic phonons with a general form of magnetoelastic coupling. For large enough frustration and lattice coupling, a new tetramerized phase with three different bond lengths is found. We argue that the zigzag spin-1/2 chain LiV2O5 might be a good candidate to observe such a phase.

Doped coupled frustrated spin-1/2 chains with four-spin exchange.

The role of various magnetic interchain couplings is investigated by numerical methods in doped frustrated quantum spin chains. A nonmagnetic dopant introduced in a gapped spin chain releases a free spin-1/2 soliton. The formation of a local magnetic moment is analyzed in terms of soliton confinement. A four-spin coupling which might originate from cyclic exchange is shown to produce such a confinement. Dopants on different chains experience an effective space-extended nonfrustrating pairwise spin interaction.

Li-induced spin and charge excitations in a spin ladder.

A lithium dopant in a cuprate spin ladder acts as a vacant (nonmagnetic) site accompanied by an extra hole bound by a Coulomb potential. We find that, although the undoped ladder spin gap is not essentially altered by Li doping, a dopant-magnon bound state appears just below the gap. Contrary to Zn-doped ladders, the spin liquid character of a Li-doped ladder should be robust against antiferromagnetism. We give a possible explanation for the similarity of Zn and Li dopants in underdoped cuprates and propose experimental tests.