It has long been believed that doped quantum spin liquids (QSLs) can give rise to fascinating quantum phases, including the possibility of high-temperature superconductivity (SC) as proposed by P. W. Anderson's resonating valence bond (RVB) scenario. The Kagome lattice $t$-$J$ model is known to exhibit spin liquid behavior at half-filling, making it an ideal system for studying the properties of doped QSL. In this study, we employ the fermionic projected entangled simplex state (PESS) method to investigate the ground state properties of the Kagome lattice $t$-$J$ model with $t/J = 3.0$. Our results reveal a phase transition from charge density wave (CDW) states to uniform states around a critical doping level $\delta_c \approx 0.27$. Within the CDW phase, we observe different types of Wigner crystal (WC) formulated by doped holes that are energetically favored. As we enter the uniform phase, a non-Fermi liquid (NFL) state emerges within the doping range $0.27 < \delta < 0.32$, characterized by an exponential decay of all correlation functions. With further hole doping, we discover the appearance of a pair density wave (PDW) state within a narrow doping region $0.32 < \delta < 1/3$. We also discuss the potential experimental implications of our findings.
In recent years, tensor network renormalization (TNR) has emerged as an efficient and accurate method for studying (1+1)D quantum systems or 2D classical systems using real-space renormalization group (RG) techniques. One notable application of TNR is its ability to extract central charge and conformal scaling dimensions for critical systems. In this paper, we present the implementation of the Loop-TNR algorithm, which allows for the computation of dynamical correlation functions. Our algorithm goes beyond traditional approaches by not only calculating correlations in the spatial direction, where the separation is an integer, but also in the temporal direction, where the time difference can contain decimal values. Our algorithm is designed to handle both imaginary-time and real-time correlations, utilizing a tensor network representation constructed from a path-integral formalism. Additionally, we highlight that the Loop-TNR algorithm can also be applied to investigate critical properties of non-Hermitian systems, an area that was previously inaccessible using density matrix renormalization group(DMRG) and matrix product state(MPS) based algorithms.
We investigate the optical behavior of a single Laguerre-Gaussian cavity optomechanical system consisting of two mechanically rotating mirrors. We explore the effects of various physical parameters on the double optomechanically induced transparency (OMIT) of the system and provide a detailed explanation of the underlying physical mechanism. We show that the momentum is not the cause of the current double-OMIT phenomena; rather, it results from the orbital angular momentum between the optical field and the rotating mirrors. Additionally, the double-OMIT is simply produced using a single Laguerre-Gaussian cavity optomechanical system rather than by integrating many subsystems or adding the atomic medium as in earlier studies. We also investigate the impact of fast and slow light in this system. Finally, we show that the switching between ultrafast and ultraslow light can be realized by adjusting the angular momentum, which is a new source of regulating fast-slow light.
We propose a minimal interacting lattice model for two-dimensional class-$D$ higher-order topological superconductors with no free-fermion counterpart. A Lieb-Schultz-Mattis-type constraint is proposed and applied to guide our lattice model construction. Our model exhibits a trivial product ground state in the weakly interacting regime, whereas, increasing electron interactions provoke a novel topological quantum phase transition to a $D_4$-symmetric higher-order topological superconducting state. The symmetry-protected Majorana corner modes are numerically confirmed with the matrix-product-state technique. Our theory paves the way for studying interacting higher-order topology with explicit lattice model constructions.
In recent years, fermionic topological phases of quantum matter has attracted a lot of attention. In a pioneer work by Gu, Wang and Wen, the concept of equivalence classes of fermionic local unitary(FLU) transformations was proposed to systematically understand non-chiral topological phases in 2D fermion systems and an incomplete classification was obtained. On the other hand, the physical picture of fermion condensation and its corresponding super pivotal categories give rise to a generic mathematical framework to describe fermionic topological phases of quantum matter. In particular, it has been pointed out that in certain fermionic topological phases, there exists the so-called q-type anyon excitations, which have no analogues in bosonic theories. In this paper, we generalize the Gu, Wang and Wen construction to include those fermionic topological phases with q-type anyon excitations. We argue that all non-chiral fermionic topological phases in 2+1D are characterized by a set of tensors $(N^{ij}_{k},F^{ij}_{k},F^{ijm,\alpha\beta}_{kln,\chi\delta},n_{i},d_{i})$, which satisfy a set of nonlinear algebraic equations parameterized by phase factors $\Xi^{ijm,\alpha\beta}_{kl}$, $\Xi^{ij}_{kln,\chi\delta}$, $\Omega^{kim,\alpha\beta}_{jl}$ and $\Omega^{ki}_{jln,\chi\delta}$. Moreover, consistency conditions among algebraic equations give rise to additional constraints on these phase factors which allow us to construct a topological invariant partition for an arbitrary triangulation of 3D spin manifold. Finally, several examples with q-type anyon excitations are discussed, including the Fermionic topological phase from Tambara-Yamagami category for $\mathbb{Z}_{2N}$, which can be regarded as the $\mathbb{Z}_{2N}$ parafermion generalization of Ising fermionic topological phase.
The nature of the zero-temperature phase diagram of the spin-$1/2$ $J_1$-$J_2$ Heisenberg model on a square lattice has been debated in the past three decades, which may hold the key to understand high temperature superconductivity. 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 the $J_1$-$J_2$ Heisenberg model up to $24\times 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\lesssim J_2/J_1\leq0.61$ before the system enters the well known collinear antiferromagnetic phase. The physical nature of the discovered gapless QSL and potential experimental implications are also addressed. 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 precisely capturing long-range physcis for highly frustrated systems, and also demonstrate the finite PEPS method is a very powerful approach to study strongly corrleated quantum many-body problems.
Periodically driven quantum systems manifest various non-equilibrium features which are absent at equilibrium. For example, discrete time-translation symmetry can be broken in periodically driven quantum systems leading to an exotic phase of matter, called discrete time crystal(DTC). For open quantum systems, previous studies showed that DTC can be found only when there exists a meta-stable state in the undriven system. However, by investigating the simplest Bose-Hubbard model with dissipation and time periodically tunneling, we find in this paper that a $2T$ DTC can appear even when the meta-stable state is absent in the undriven system. This observation extends the understanding of DTC and shed more light on the physics behind the DTC. Besides, by the detailed analysis of simplest two-sites model, we show further that the two-sites model can be used as basic building blocks to construct large rings in which a $nT$ DTC might appear. These results might find applications into engineering exotic phases in driven open quantum systems.
Based on the scheme of variational Monte Carlo sampling, we develop an accurate and efficient two-dimensional tensor-network algorithm to simulate quantum lattice models. We find that Monte Carlo sampling shows huge advantages in dealing with finite projected entangled pair states, which allows significantly enlarged system size and improves the accuracy of tensor network simulation. We demonstrate our method on the square-lattice antiferromagnetic Heisenberg model up to $32 \times 32$ sites, as well as a highly frustrated $J_1-J_2$ model up to $24\times 24$ sites. The results, including ground state energy and spin correlations, are in excellent agreement with those of the available quantum Monte Carlo or density matrix renormalization group methods. Therefore, our method substantially advances the calculation of 2D tensor networks for finite systems, and potentially opens a new door towards resolving many challenging strongly correlated quantum many-body problems.
In this paper, we construct $2n-1$ locally indistinguishable orthogonal product states in $\mathbb{C}^n\otimes\mathbb{C}^{4}~(n>4)$ and $\mathbb{C}^n\otimes\mathbb{C}^{5}~(n\geq 5)$ respectively. Moreover, a set of locally indistinguishable orthogonal product states with $2(n+2l)-8$ elements in $\mathbb{C}^n\otimes\mathbb{C}^{2l}~(n\geq 2l>4)$ and a class of locally indistinguishable orthogonal product states with $2(n+2k+1)-7$ elements in $\mathbb{C}^n\otimes\mathbb{C}^{2k+1}~(n\geq 2k+1>5)$ are also constructed respectively. These classes of quantum states are then shown to be distinguishable by local operation and classical communication (LOCC) using a suitable $\mathbb{C}^2\otimes\mathbb{C}^2$ maximally entangled state respectively.
Subnatural-linewidth single-photon source is a potential candidate for exploring the time degree of freedom in photonic quantum information science. This type of single-photon source has been demonstrated to be generated and reshaped in atomic ensembles without any external cavity or filter, and is typically characterized through photon-counting technology. However, the full complex temporal mode function(TMF) of the photon source is not able to be revealed from direct photon counting measurement. Here, for the first time, we demonstrate the complete temporal mode of the subnatural-linewidth single photons generated from a cold atomic cloud. Through heterodyne detection between the single photon and a local oscillator with various central frequencies, we recover the temporal density matrix of the single photons at resolvable time bins. Further we demonstrate that the reduced autocorrelation function measured through homodyne detection perfectly reveals the pure temporal-spectral state of the subnatural-linewidth single photons.
We propose to utilize the sub-system fidelity (SSF), defined by comparing a pair of reduced density matrices derived from the degenerate ground states, to identify and/or characterize symmetry protected topological (SPT) states in one-dimensional interacting many-body systems. The SSF tells whether two states are locally indistinguishable (LI) by measurements within a given sub-system. Starting from two polar states (states that could be distinguished on either edge), the other combinations of these states can be mapped onto a Bloch sphere. We prove that a pair of orthogonal states on the equator of the Bloch sphere are LI, independently of whether they are SPT states or cat states (symmetry-preserving states by linear combinations of states that break discrete symmetries). Armed with this theorem, we provide a scheme to construct zero-energy exitations that swap the LI states. We show that the zero mode can be located anywhere for cat states, but is localized near the edge for SPT states. We also show that the SPT states are LI in a finite fraction of the bulk (excluding the two edges), whereas the symmetry-breaking states are distinguishable. This can be used to pinpoint the transition from SPT states to the symmetry-breaking states.
In this paper, we mainly characterize the structure of product bases of the complex vector space $\mathbb{C}^{2}\bigotimes\mathbb{C}^{n}$. It gives an answer to the conjecture in case of $d=2n$ proposed by McNulty et al in 2016. As the application of the result, we obtain all the product bases of a bipartite system $\mathbb{C}^{2}\bigotimes\mathbb{C}^{n}$. It is helpful to review the structure of all the product bases of $\mathbb{C}^{2}\bigotimes\mathbb{C}^{2}$ and $\mathbb{C}^{2}\bigotimes\mathbb{C}^{3}$, which given by McNulty et al.
We introduce a tensor renormalization group scheme for coarse-graining a two-dimensional tensor network that can be successfully applied to both classical and quantum systems on and off criticality. The key innovation in our scheme is to deform a 2D tensor network into small loops and then optimize the tensors on each loop. In this way, we remove short-range entanglement at each iteration step and significantly improve the accuracy and stability of the renormalization flow. We demonstrate our algorithm in the classical Ising model and a frustrated 2D quantum model.
Quantum-disordering a discrete-symmetry breaking state by condensing domain-walls can lead to a trivial symmetric insulator state. In this work, we show that if we bind a 1D representation of the symmetry (such as a charge) to the intersection point of several domain walls, condensing such modified domain-walls can lead to a non-trivial symmetry-protected topological (SPT) state. This result is obtained by showing that the modified domain-wall condensed state has a non-trivial SPT invariant -- the symmetry-twist dependent partition function. We propose two different kinds of field theories that can describe the above mentioned SPT states. The first one is a Ginzburg-Landau-type non-linear sigma model theory, but with an additional multi-kink domain-wall topological term. Such theory has an anomalous $U^k(1)$ symmetry but an anomaly-free $Z_N^k$ symmetry. The second one is a gauge theory, which is beyond Abelian Chern-Simons/BF gauge theories. We argue that the two field theories are equivalent at low energies. After coupling to the symmetry twists, both theories produce the desired SPT invariant.
Bosonic topological insulators (BTI) in three dimensions are symmetry-protected topological phases (SPT) protected by time-reversal and boson number conservation symmetries. BTI in three dimensions were first proposed and classified by the group cohomology theory which suggests two distinct root states, each carrying a $\mathbb{Z}_2$ index. Soon after, surface anomalous topological orders were proposed to identify different root states of BTI, which even leads to a new BTI root state beyond the group cohomology classification. In this paper, we propose a universal physical mechanism via \textitvortex-line condensation from a 3d superfluid to achieve all three root states. It naturally produces bulk topological quantum field theory (TQFT) description for each root state. Topologically ordered states on the surface are \textitrigorously derived by placing TQFT on an open manifold, which allows us to explicitly demonstrate the bulk-boundary correspondence. Finally, we generalize the mechanism to $Z_N$ symmetries and discuss potential SPT phases beyond the group cohomology classification.
Searching for $p+ip$ superconducting(SC) state has become a fascinating subject in condensed matter physics, as a dream application awaiting in topological quantum computation. In this paper, we report a theoretical discovery of a $p+ip$ SC ground state (coexisting with ferromagnetic order) in honeycomb lattice Hubbard model with infinite repulsive interaction at low doping($\delta< 0.2$), by using both the state-of-art Grassmann tensor product state(GTPS) approach and a quantum field theory approach. Our discovery suggests a new mechanism for $p+ip$ SC state in generic strongly correlated systems and opens a new door towards experimental realization. The $p+ip$ SC state has an instability towards a potential non-Fermi liquid with a large but finite $U$. However, a small Zeeman field term stabilizes the $p+ip$ SC state. Relevant realistic materials are also proposed.
The Z_2 topological order in Z_2 spin liquid and in exactly soluble Kitaev toric code model is the simplest topological order for 2+1D bosonic systems. More general 2+1D bosonic topologically ordered states can be constructed via exact soluble string-net models. However, the most important topologically ordered phases of matter are arguably the fermionic fractional quantum Hall states. Topological phases of matter for fermion systems are strictly richer than their bosonic counterparts because locality has different meanings for the two kinds of systems. In this paper, we describe a simple fermionic version of the toric code model to illustrate many salient features of fermionic exactly soluble models and fermionic topologically ordered states.
We study the effect of interactions on 2D fermionic symmetry-protected topological (SPT) phases using the recently proposed braiding statistics approach. We focus on a simple class of examples: superconductors with a Z2 Ising symmetry. Although these systems are classified by Z in the noninteracting limit, our results suggest that the classification collapses to Z8 in the presence of interactions -- consistent with previous work that analyzed the stability of the edge. Specifically, we show that there are at least 8 different types of Ising superconductors that cannot be adiabatically connected to one another, even in the presence of strong interactions. In addition, we prove that each of the 7 nontrivial superconductors have protected edge modes.
We report the theoretical discovery of a systematic scheme to produce topological flat bands (TFBs) with arbitrary Chern numbers. We find that generically a multi-orbital high Chern number TFB model can be constructed by considering multi-layer Chern number C=1 TFB models with enhanced translational symmetry. A series of models are presented as examples, including a two-band model on a triangular lattice with a Chern number C=3 and an $N$-band square lattice model with $C=N$ for an arbitrary integer $N$. In all these models, the flatness ratio for the TFBs is larger than 30 and increases with increasing Chern number. In the presence of appropriate inter-particle interactions, these models are likely to lead to the formation of novel Abelian and Non-Abelian fractional Chern insulators. As a simple example, we test the C=2 model with hardcore bosons at 1/3 filling and an intriguing fractional quantum Hall state is observed.
Symmetry-protected topological (SPT) phases are gapped short-range-entangled quantum phases with a symmetry $G$, which can all be smoothly connected to the trivial product states if we break the symmetry. It has been shown that a large class of interacting bosonic SPT phases can be systematically described by group cohomology theory. In this paper, we introduce a (special) group supercohomology theory which is a generalization of the standard group cohomology theory. We show that a large class of short-range interacting fermionic SPT phases can be described by the group supercohomology theory. Using the data of super cocycles, we can obtain the ideal ground state wave function for the corresponding fermionic SPT phase. We can also obtain the bulk Hamiltonian that realizes the SPT phase, as well as the anomalous (ie, non-on-site) symmetry for the boundary effective Hamiltonian. The anomalous symmetry on the boundary implies that the symmetric boundary must be gapless for 1+1D boundary, and must be gapless or topologically ordered beyond 1+1D. As an application of this general result, we construct a new SPT phase in 3D, for interacting fermionic superconductors with coplanar spin order (which have $T^2=1$ time-reversal $Z_2^T$ and fermion-number parity $Z_2^f$ symmetries described by a full symmetry group $Z_2^T\times Z_2^f$). Such a fermionic SPT state can neither be realized by free fermions nor by interacting bosons (formed by fermion-pairs), and thus are not included in the K-theory classification for free fermions or group cohomology description for interacting bosons. We also construct three interacting fermionic SPT phases in 2D with a full symmetry group $Z_2\times Z_2^f$. Those 2D fermionic SPT phases all have central-charge $c=1$ gapless edge excitations, if the symmetry is not broken.
The ground state phase of spin-1/2 $J_1$-$J_2$ antiferromagnetic Heisenberg model on square lattice around the maximally frustrated regime ($J_2\sim 0.5J_1$) has been debated for decades. Here we study this model using the cluster update algorithm for tensor product states (TPSs). The ground state energies at finite sizes and in the thermodynamic limit (with finite size scaling) are in good agreement with exact diagonalization study. Through finite size scaling of the spin correlation function, we find the critical point $J_2^{c_1}=0.572(5)J_1$ and critical exponents $\nu=0.50(8)$, $\eta_s=0.28(6)$. In the range of $0.572 < J_2/J_1 \leqslant 0.6 $ we find a paramagnetic ground state with exponentially decaying spin-spin correlation. Up to $24\times 24$ system size, we observe power law decaying dimer-dimer and plaquette-plaquette correlations with an anomalous plaquette scaling exponent $\eta_p=0.24(1)$ and an anomalous columnar scaling exponent $\eta_c=0.28(1)$ at $J_2/J_1=0.6$. These results are consistent with a potential gapless $U(1)$ spin liquid phase. However, since the $U(1)$ spin liquid is unstable due to the instanton effect, a VBS order with very small amplitude might develop in the thermodynamic limit. Thus, our numerical results strongly indicate a deconfined quantum critical point (DQCP) at $J_2^{c_1}$. Remarkably, all the observed critical exponents are consistent with the $J-Q$ model.
Inspired by recent theoretical discovery of robust fractional topological phases without a magnetic field, we search for the non-Abelian quantum Hall effect (NA-QHE) in lattice models with topological flat bands (TFBs). Through extensive numerical studies on the Haldane model with three-body hard-core bosons loaded into a TFB, we find convincing numerical evidence of a stable $\nu=1$ bosonic NA-QHE, with the characteristic three-fold quasi-degeneracy of ground states on a torus, a quantized Chern number, and a robust spectrum gap. Moreover, the spectrum for two-quasihole states also shows a finite energy gap, with the number of states in the lower energy sector satisfying the same counting rule as the Moore-Read Pfaffian state.
Recently, the Grassmann-tensor-entanglement renormalization group(GTERG) approach was proposed as a generic variational approach to study strongly correlated boson/fermion systems. However, the weakness of such a simple variational approach is that generic Grassmann tensor product states(GTPS) with large inner dimension $D$ will contain a large number of variational parameters and be hard to be determined through usual minimization procedures. In this paper, we first introduce a standard form of GTPS which significantly simplifies the representations. Then we describe a simple imaginary-time-evolution algorithm to efficiently update the GTPS based on the fermion coherent state representation and show all the algorithm developed for usual tensor product states(TPS) can be implemented for GTPS in a similar way. Finally, we study the environment effect for the GTERG approach and propose a simple method to further improve its accuracy. We demonstrate our algorithms by studying a simple 2D free fermion system on honeycomb lattice, including both off-critical and critical cases.
Symmetry protected topological (SPT) phases are gapped short-range-entangled quantum phases with a symmetry G. They can all be smoothly connected to the same trivial product state if we break the symmetry. The Haldane phase of spin-1 chain is the first example of SPT phase which is protected by SO(3) spin rotation symmetry. The topological insulator is another exam- ple of SPT phase which is protected by U(1) and time reversal symmetries. It has been shown that free fermion SPT phases can be systematically described by the K-theory. In this paper, we show that interacting bosonic SPT phases can be systematically described by group cohomology theory: distinct d-dimensional bosonic SPT phases with on-site symmetry G (which may contain anti-unitary time reversal symmetry) can be labeled by the elements in H^1+d[G, U_T(1)] - the Borel (1 + d)-group-cohomology classes of G over the G-module U_T(1). The boundary excitations of the non-trivial SPT phases are gapless or degenerate. Even more generally, we find that the different bosonic symmetry breaking short-range-entangled phases are labeled by the following three mathematical objects: (G_H, G_\Psi, H^1+d[G_\Psi, U_T(1)], where G_H is the symmetry group of the Hamiltonian and G_\Psi the symmetry group of the ground states.
Quantum phases with different orders exist with or without breaking the symmetry of the system. Recently, a classification of gapped quantum phases which do not break time reversal, parity or on-site unitary symmetry has been given for 1D spin systems in [X. Chen, Z.-C. Gu, and X.-G. Wen, Phys. Rev. B \textbf83, 035107 (2011); arXiv:1008.3745]. It was found that, such symmetry protected topological (SPT) phases are labeled by the projective representations of the symmetry group which can be viewed as a symmetry fractionalization. In this paper, we extend the classification of 1D gapped phases by considering SPT phases with combined time reversal, parity, and/or on-site unitary symmetries and also the possibility of symmetry breaking. We clarify how symmetry fractionalizes with combined symmetries and also how symmetry fractionalization coexists with symmetry breaking. In this way, we obtain a complete classification of gapped quantum phases in 1D spin systems. We find that in general, symmetry fractionalization, symmetry breaking and long range entanglement(present in 2 or higher dimensions) represent three main mechanisms to generate a very rich set of gapped quantum phases. As an application of our classification, we study the possible SPT phases in 1D fermionic systems, which can be mapped to spin systems by Jordan-Wigner transformation.
Recent proposals of topological flat band (TFB) models have provided a new route to realize the fractional quantum Hall effect (FQHE) without Landau levels. We study hard-core bosons with short-range interactions in two representative TFB models, one of which is the well known Haldane model (but with different parameters). We demonstrate that FQHE states emerge with signatures of even number of quasi-degenerate ground states on a torus and a robust spectrum gap separating these states from higher energy spectrum. We also establish quantum phase diagrams for the filling factor 1/2 and illustrate quantum phase transitions to other competing symmetry-breaking phases.
The string-net approach by Levin and Wen, and the local unitary transformation approach by Chen, Gu, and Wen, provide ways to classify topological orders with gappable edge in 2D bosonic systems. The two approaches reveal that the mathematical framework for 2+1D bosonic topological order with gappable edge is closely related to unitary fusion category theory. In this paper, we generalize these systematic descriptions of topological orders to 2D fermion systems. We find a classification of 2+1D fermionic topological orders with gappable edge in terms of the following set of data $(N^{ij}_k, F^{ij}_k, F^{ijm,\alpha\beta}_{jkn,\chi\delta},d_i)$, that satisfy a set of non-linear algebraic equations. The exactly soluble Hamiltonians can be constructed from the above data on any lattices to realize the corresponding topological orders. When $F^{ij}_k=0$, our result recovers the previous classification of 2+1D bosonic topological orders with gappable edge.
Quantum many-body systems divide into a variety of phases with very different physical properties. The question of what kind of phases exist and how to identify them seems hard especially for strongly interacting systems. Here we make an attempt to answer this question for gapped interacting quantum spin systems whose ground states are short-range correlated. Based on the local unitary equivalence relation between short-range correlated states in the same phase, we classify possible quantum phases for 1D matrix product states, which represent well the class of 1D gapped ground states. We find that in the absence of any symmetry all states are equivalent to trivial product states, which means that there is no topological order in 1D. However, if certain symmetry is required, many phases exist with different symmetry protected topological orders. The symmetric local unitary equivalence relation also allows us to obtain some simple results for quantum phases in higher dimensions when some symmetries are present.
Two gapped quantum ground states in the same phase are connected by an adiabatic evolution which gives rise to a local unitary transformation that maps between the states. On the other hand, gapped ground states remain within the same phase under local unitary transformations. Therefore, local unitary transformations define an equivalence relation and the equivalence classes are the universality classes that define the different phases for gapped quantum systems. Since local unitary transformations can remove local entanglement, the above equivalence/universality classes correspond to pattern of long range entanglement, which is the essence of topological order. The local unitary transformation also allows us to define a wave function renormalization scheme, under which a wave function can flow to a simpler one within the same equivalence/universality class. Using such a setup, we find conditions on the possible fixed-point wave functions where the local unitary transformations have \emphfinite dimensions. The solutions of the conditions allow us to classify this type of topological orders, which generalize the string-net classification of topological orders. We also describe an algorithm of wave function renormalization induced by local unitary transformations. The algorithm allows us to calculate the flow of tensor-product wave functions which are not at the fixed points. This will allow us to calculate topological orders as well as symmetry breaking orders in a generic tensor-product state.
The projective construction (the slave-particle approach) has played an very important role in understanding strongly correlated systems, such as the emergence of fermions, anyons, and gauge theory in quantum spin liquids and quantum Hall states. Recently, fermionic Projected Entangled Pair States (fPEPS) have been introduced to effciently represent many-body fermionic states. In this paper, we show that the strongly correlated bosonic/fermionic states obtained both from the projective construction and the fPEPS approach can be represented systematically as Grassmann tensor product states. This construction can also be applied to all other tensor network states approaches. The Grassmann tensor product states allow us to encode many-body bosonic/fermionic states effciently with a polynomial number of parameters. We also generalize the tensor-entanglement renormalization group (TERG) method for complex tensor networks to Grassmann tensor networks. This allows us to approximate the norm and average local operators of Grassmann tensor product states in polynomial time, and hence leads to a variational approach for describing strongly correlated bosonic/fermionic systems in higher dimensions.
The tensor product representation of quantum states leads to a promising variational approach to study quantum phase and quantum phase transitions, especially topological ordered phases which are impossible to handle with conventional methods due to their long range entanglement. However, an important issue arises when we use tensor product states (TPS) as variational states to find the ground state of a Hamiltonian: can arbitrary variations in the tensors that represent ground state of a Hamiltonian be induced by local perturbations to the Hamiltonian? Starting from a tensor product state which is the exact ground state of a Hamiltonian with $\mathbb{Z}_2$ topological order, we show that, surprisingly, not all variations of the tensors correspond to the variation of the ground state caused by local perturbations of the Hamiltonian. Even in the absence of any symmetry requirement of the perturbed Hamiltonian, one necessary condition for the variations of the tensors to be physical is that they respect certain $\mathbb{Z}_2$ symmetry. We support this claim by calculating explicitly the change in topological entanglement entropy with different variations in the tensors. This finding will provide important guidance to numerical variational study of topological phase and phase transitions. It is also a crucial step in using TPS to study universal properties of a quantum phase and its topological order.
Many-body entangled quantum states studied in condensed matter physics can be primary resources for quantum information, allowing any quantum computation to be realized using measurements alone, on the state. Such a universal state would be remarkably valuable, if only it were thermodynamically stable and experimentally accessible, by virtue of being the unique ground state of a physically reasonable Hamiltonian made of two-body, nearest neighbor interactions. We introduce such a state, composed of six-state particles on a hexagonal lattice, and describe a general method for analyzing its properties based on its projected entangled pair state representation.
Traditional mean-field theory is a simple generic approach for understanding various phases. But that approach only applies to symmetry breaking states with short-range entanglement. In this paper, we describe a generic approach for studying 2D quantum phases with long-range entanglement (such as topological phases). Our approach is a variational method that uses tensor product states (also known as projected entangled pair states) as trial wave functions. We use a 2D real space RG algorithm to evaluate expectation values in these wave functions. We demonstrate our algorithm by studying several simple 2D quantum spin models.