We study the antiferromagnetic kagome Heisenberg model with additional scalar-chiral interaction by using the infinite projected entangled-pair state (iPEPS) ansatz. We discuss in detail the implementation of optimization algorithm in the framework of the single-layer tensor network based on the corner-transfer matrix technique. Our benchmark based on the full-update algorithm shows that the single-layer algorithm is stable, which leads to the same level of accuracy as the double-layer ansatz but with much less computation time. We further apply this algorithm to study the nature of the kagome Heisenberg model with a scalar-chiral interaction by computing the bond dimension scaling of magnetization, bond energy difference, chiral order parameter and correlation length. In particular, we find that for strong chiral coupling the correlation length, which is extracted from the transfer matrix, saturates to a finite value for large bond dimension, representing a gapped spin-liquid state. Further comparison with density matrix renormalization group results supports that our iPEPS faithfully represents the time-reversal symmetry breaking chiral state. Our iPEPS simulation results shed new light on constructing PEPS for describing gapped chiral topological states.
We study the antiferromagnetic kagome Heisenberg model with additional scalar-chiral interaction by using the infinite projected entangled-pair state (iPEPS) ansatz. We discuss in detail the implementation of optimization algorithm in the framework of the single-layer tensor network based on the corner-transfer matrix technique. Our benchmark based on the full-update algorithm shows that the single-layer algorithm is stable, which leads to the same level of accuracy as the double-layer ansatz but with much less computation time. We further apply this algorithm to study the nature of the kagome Heisenberg model with a scalar-chiral interaction by computing the bond dimension scaling of magnetization, bond energy difference, chiral order parameter and correlation length. In particular, we find that for strong chiral coupling the correlation length, which is extracted from the transfer matrix, saturates to a finite value for large bond dimension, representing a gapped spin-liquid state. Further comparison with density matrix renormalization group results supports that our iPEPS faithfully represents the time-reversal symmetry breaking chiral state. Our iPEPS simulation results shed new light on constructing PEPS for describing gapped chiral topological states.
We study the antiferromagnetic kagome Heisenberg model with additional scalar-chiral interaction by using the infinite projected entangled-pair state (iPEPS) ansatz. We discuss in detail the implementation of optimization algorithm in the framework of the single-layer tensor network based on the corner-transfer matrix technique. Our benchmark based on the full-update algorithm shows that the single-layer algorithm is stable, which leads to the same level of accuracy as the double-layer ansatz but with much less computation time. We further apply this algorithm to study the nature of the kagome Heisenberg model with a scalar-chiral interaction by computing the bond dimension scaling of magnetization, bond energy difference, chiral order parameter and correlation length. In particular, we find that for strong chiral coupling the correlation length, which is extracted from the transfer matrix, saturates to a finite value for large bond dimension, representing a gapped spin-liquid state. Further comparison with density matrix renormalization group results supports that our iPEPS faithfully represents the time-reversal symmetry breaking chiral state. Our iPEPS simulation results shed new light on constructing PEPS for describing gapped chiral topological states.
We study the antiferromagnetic kagome Heisenberg model with an additional scalar-chiral interaction by using the infinite projected entangled-pair state (iPEPS) ansatz. We discuss in detail the implementation of the optimization algorithm in the framework of the single-layer tensor network based on the corner transfer matrix technique. Our benchmark based on the full-update algorithm shows that the single-layer algorithm is stable, which leads to the same level of accuracy as the double-layer ansatz but with much less computation time. We further apply this algorithm to study the nature of the kagome Heisenberg model with a scalar-chiral interaction by computing the bond dimension scaling of magnetization, bond energy difference, chiral order parameter, and correlation length. In particular, we find that for strong chiral coupling the correlation length, which is extracted from the transfer matrix, saturates to a finite value for large bond dimension, representing a gapped spin-liquid state. Further comparison with density matrix renormalization group results supports that our iPEPS faithfully represents the time-reversal symmetry breaking chiral state. Our iPEPS simulation results shed new light on constructing PEPS for describing gapped chiral topological states.
The infinite projected entangled simplex state (iPESS), a type of tensor network (TN) state, has been used successfully for simulating and characterizing {\it non-chiral} spin liquids on the kagome lattice. Here, we demonstrate that iPESS also provides a faithful representation of a {\it chiral} spin liquid (CSL) on the same lattice, namely the ground state of the spin-1/2 kagome Heisenberg antiferromagnet with a scalar chirality. By classifying local tensors according to SU(2) and point group symmetries, we construct a chiral ansatz breaking reflection P and time reversal T symmetries while preserving PT. The variational TN states are shown to host, for bond dimension D≥8, a chiral gapless entanglement spectrum following SU(2)1 conformal field theory. The correlation function shows a small weight long-range tail complying with the prediction of the TN bulk-edge correspondence. %{\color{red}With more constraints included, the chiral ansatz is reduced to a non-chiral one which preserves full point group symmetries due to an emergent {\it tensor conservation law} and is of relevance to the ground state at the Heisenberg point. Lastly, by simulations in the complete ansatz family we discuss the transition from the non-chiral spin liquid to the CSL induced by the scalar chirality term.} We identify a non-chiral manifold spanned by only a subset of symmetric tensors where a new emergent {\it tensor conservation law} is realized. This allows us to both probe the stability of the non-chiral spin liquid and discuss its transition to CSL induced by a scalar chirality term.
The infinite projected entangled simplex state (iPESS), a type of tensor network (TN) state, has been used successfully for simulating and characterizing nonchiral spin liquids on the kagome lattice. Here, we demonstrate that iPESS also provides a faithful representation of a chiral spin liquid (CSL) on the same lattice, namely the ground state of the spin-1/2 kagome Heisenberg antiferromagnet with a scalar chirality. By classifying local tensors according to SU(2) and point group symmetries, we construct a chiral ansatz breaking reflection P and time reversal T symmetries while preserving PT. The variational TN states are shown to host, for bond dimension $D\ensuremath{\ge}8$, a chiral gapless entanglement spectrum following SU(2)1 conformal field theory. The correlation function shows a small weight long-range tail complying with the prediction of the TN bulk-edge correspondence. We identify a nonchiral manifold spanned by only a subset of symmetric tensors where a new emergent tensor conservation law is realized. This allows us to both probe the stability of the nonchiral spin liquid and discuss its transition to CSL induced by a scalar chirality term.
We develop an improved variant of U(1)-symmetric infinite projected entangled-pair states (iPEPS) ansatz to investigate the ground-state phase diagram of the spin-1/2 square ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model. In order to improve the accuracy of the ansatz, we discuss a simple strategy to select automatically relevant symmetric sectors and also introduce an optimization method to treat second-neighbor interactions more efficiently. We show that variational ground-state energies of the model obtained by the U(1)-symmetric iPEPS ansatz (for a fixed bond dimension D) set a better upper bound, improving previous tensor-network-based results. By studying the finite-D scaling of the magnetically order parameter, we find a N\'eel phase for J2/J1<0.53. For 0.53<J2/J1<0.61, a nonmagnetic columnar valence bond solid (VBS) state is established as observed by the pattern of local bond energy. The divergent behavior of correlation length $\ensuremath{\xi}\ensuremath{\sim}{D}^{1.2}$ and vanishing order parameters are consistent with a deconfined N\'eel-to-VBS transition at J2c1/J1=0.530(5), where estimated critical anomalous exponents are ${\ensuremath{\eta}}_{s}\ensuremath{\sim}0.6$ and ${\ensuremath{\eta}}_{d}\ensuremath{\sim}1.9$ for spin and dimer correlations, respectively. We show that the associated VBS order parameter monotonically increases with J2/J1 and finally a first-order quantum phase transition takes place at J2c2/J1=0.610(2) to the conventional Stripe phase. We compare our results with earlier DMRG and PEPS studies and suggest future directions for resolving remaining issues.
A typical quantum state obeying the area law for entanglement on an infinite two-dimensional (2D) lattice can be represented by a tensor network Ansatz, known as an infinite projected entangled pair state (iPEPS), with a finite bond dimension D. Its real-imaginary time evolution can be split into small time steps. An application of a time step generates a new iPEPS with a bond dimension k times the original one. The new iPEPS does not make optimal use of its enlarged bond dimension kD; hence, in principle, it can be represented accurately by a more compact Ansatz, preferably with the original D. In this work we show how the more compact iPEPS can be optimized variationally to maximize its overlap with the new iPEPS. To compute the overlap we use the corner-transfer-matrix renormalization group. By simulating sudden quench of the transverse field in the 2D quantum Ising model with the proposed algorithm, we provide a proof of principle that real-time evolution can be simulated with iPEPS. A similar proof is provided with the same model for imaginary-time evolution of purification of its thermal states.