In this thesis we extend the formalism of tensor network algorithms to incorporate global internal symmetries. We describe how to both numerically protect the symmetry and exploit it for computational gain in tensor network simulations. Our formalism is independent of the details of a specific tensor network decomposition since the symmetry constraints are imposed at the level of individual tensors. Moreover, the formalism can be applied to a wide spectrum of physical symmetries described by any discrete or continuous group that is compact and reducible. We describe in detail the implementation of the conservation of total particle number (U(1) symmetry) and of total angular momentum (SU(2) symmetry). Our formalism can also be readily generalized to incorporate more exotic symmetries such as conservation of total charge in anyonic systems.
This is the documentation for the tensor library QSpace (v4.0) that provides a toolbox to exploit `quantum symmetry spaces’ in tensor network states in the quantum many-body context. QSpace permits arbitrary combinations of symmetries including the Abelian symmetries \mathbb{Z}_n ℤ n and U(1) U ( 1 ) , as well as all non-Abelian symmetries represented by the simple classical Lie algebras: A_n A n , B_n B n , C_n C n , and D_n D n , or respectively, the special unitary group SU( n n ), the odd orthogonal group SO( 2n+1 2 n + 1 ), the symplectic group Sp( 2n 2 n ), and the even orthogonal group SO( 2n 2 n ). The code (C++ embedded via the MEX interface into Matlab) is available open-source as of QSpace v4.0 on bitbucket (https://bitbucket.org/qspace4u) under the Apache 2.0 license. QSpace is designed as a bottom-up approach for non-Abelian symmetries. It starts from the defining representation and the respective Lie algebra. By explicitly computing and tabulating generalized Clebsch-Gordan coefficient tensors, QSpace is versatile in its operations across all symmetries. At the level of an application, much of the symmetry-related details are hidden within the QSpace C++ core libraries. Hence when developing tensor network algorithms with QSpace, these can be coded (nearly) as if there are no symmetries at all, despite being able to fully exploit general non-Abelian symmetries.
The benefits of exploiting the presence of symmetries in tensor network algorithms have been extensively demonstrated in the context of matrix product states (MPSs). These include the ability to select a specific symmetry sector (e.g., with a given particle number or spin), to ensure the exact preservation of total charge, and to significantly reduce computational costs. Compared to the case of a generic tensor network, the practical implementation of symmetries in the MPS is simplified by the fact that tensors only have three indices (they are trivalent, just as the Clebsch-Gordan coefficients of the symmetry group) and are organized as a one-dimensional array of tensors, without closed loops. Instead, a more complex tensor network, one where tensors have a larger number of indices and/or a more elaborate network structure, requires a more general treatment. In two recent papers, namely, (i) [Singh, Pfeifer, and Vidal, Phys. Rev. A 82, 050301 (2010)] and (ii) [Singh, Pfeifer, and Vidal, Phys. Rev. B 83, 115125 (2011)], we described how to incorporate a global internal symmetry into a generic tensor network algorithm based on decomposing and manipulating tensors that are invariant under the symmetry. In (i) we considered a generic symmetry group
G that is compact, completely reducible, and multiplicity free, acting as a global internal symmetry. Then, in (ii) we described the implementation of Abelian group symmetries in much more detail, considering a U(1) symmetry (e.g., conservation of global particle number) as a concrete example. In this paper, we describe the implementation of non-Abelian group symmetries in great detail. For concreteness, we consider an SU(2) symmetry (e.g., conservation of global quantum spin). Our formalism can be readily extended to more exotic symmetries associated with conservation of total fermionic or anyonic charge. As a practical demonstration, we describe the SU(2)-invariant version of the multiscale entanglement renormalization ansatz and apply it to study the low-energy spectrum of a quantum spin chain with a global SU(2) symmetry.
In this thesis I develop a formalism whereby a tensor network may be understood in terms of a unitary braided tensor category, and represented in a particularly efficient manner corresponding to the exploitation of this mathematical structure. This approach permits the exploitation of the global Abelian and non-Abelian internal group symmetries of a lattice model, both to facilitate the study of particular symmetry sectors and to provide a considerable reduction in computational cost. It also enables the study of models with non-trivial exchange statistics (fermions, Abelian and non-Abelian anyons) using tensor network algorithms which scale polynomially in the system size.