Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

SciPost Physics Codebases

Codebase release 4.0 for QSpace

Andreas Weichselbaum

Published Nov 14, 2024• 6 citations
Publisher

Abstract

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 (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.

DOI
10.21468/scipostphyscodeb.40-r4.0

Related papers

QSpace - An open-source tensor library for Abelian and non-Abelian symmetries

Andreas Weichselbaum

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.

arXiv:2405.06632•SciPost Physics Codebases•Nov 14, 2024•23 citations

Tensor network states and algorithms in the presence of a global SU(2) symmetry

Sukhwinder Singh, Guifré Vidal

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 $\mathcal{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.

arXiv:1208.3919•Physical Review B•Nov 9, 2012•114 citations

SeQuant Framework for Symbolic and Numerical Tensor Algebra. I. Core Capabilities

Bimal Gaudel, Robert G. Adam, Ajay Melekamburath +5 more

SeQuant is an open-source library for symbolic algebra of tensors over commutative (scalar) and non-commutative (operator) rings. The key innovation supporting most of its functionality is a graph-theoretic tensor network (TN) canonicalizer that can handle tensor networks with symmetries faster than their standard group-theoretic counterparts. The TN canonicalizer is used for routine simplification of conventional tensor expressions, for optimizing application of Wick's theorem (used to canonicalize products of tensors over operator fields), and for manipulation of the intermediate representation leading to the numerical evaluation. Notable features of SeQuant include support for noncovariant tensor networks (which often arise from tensor decompositions) and for tensors with modes that depend parametrically on indices of other tensor modes (such dependencies between degrees of freedom are naturally viewed as nesting of tensors, "tensors of tensors" arising in block-wise data compressions in data science and modern quantum simulation). SeQuant blurs the line between pure symbolic manipulation/code generation and numerical evaluation by including compiler-like components to optimize and directly interpret tensor expressions using external numerical tensor algebra frameworks. The SeQuant source code is available at https://github.com/ValeevGroup/SeQuant.

cs.MSCross-listedarXiv:2511.09943•The Journal of Chemical Physics•Nov 14, 2025•1 citation

Tensor network states and algorithms in the presence of a global U(1) symmetry

Sukhwinder Singh, Robert N. C. Pfeifer, Guifré Vidal

Tensor network decompositions offer an efficient description of certain many-body states of a lattice system and are the basis of a wealth of numerical simulation algorithms. In a recent paper [Phys. Rev. A 82, 050301 (2010)] we discussed how to incorporate a global internal symmetry, given by a compact, completely reducible group $\mathcal{G}$, into tensor network decompositions and algorithms. Here we specialize to the case of Abelian groups and, for concreteness, to a U(1) symmetry, associated, e.g., with particle number conservation. We consider tensor networks made of tensors that are invariant (or covariant) under the symmetry, and explain how to decompose and manipulate such tensors in order to exploit their symmetry. In numerical calculations, the use of U(1)-symmetric tensors allows selection of a specific number of particles, ensures the exact preservation of particle number, and significantly reduces computational costs. We illustrate all these points in the context of the multiscale entanglement renormalization Ansatz.

arXiv:1008.4774•Physical Review B•Mar 15, 2011•287 citations

A general non-Abelian density matrix renormalization group algorithm with application to the C2 dimer

Sandeep Sharma

We extend our previous work [S. Sharma and G. K.-L. Chan, J. Chem. Phys. 136, 124121 (2012)], which described a spin-adapted (SU(2) symmetry) density matrix renormalization group algorithm, to additionally utilize general non-Abelian point group symmetries. A key strength of the present formulation is that the requisite tensor operators are not hard-coded for each symmetry group, but are instead generated on the fly using the appropriate Clebsch-Gordan coefficients. This allows our single implementation to easily enable (or disable) any non-Abelian point group symmetry (including SU(2) spin symmetry). We use our implementation to compute the ground state potential energy curve of the C2 dimer in the cc-pVQZ basis set (with a frozen-core), corresponding to a Hilbert space dimension of 10(12) many-body states. While our calculated energy lies within the 0.3 mEh error bound of previous initiator full configuration interaction quantum Monte Carlo and correlation energy extrapolation by intrinsic scaling calculations, our estimated residual error is only 0.01 mEh, much more accurate than these previous estimates. Due to the additional efficiency afforded by the algorithm, the excitation energies (Te) of eight lowest lying excited states: a(3)Πu, b(3)Σg (-), A(1)Πu, c(3)Σu (+), B(1)Δg, B(') (1)Σg (+), d(3)Πg, and C(1)Πg are calculated, which agree with experimentally derived values to better than 0.06 eV. In addition, we also compute the potential energy curves of twelve states: the three lowest levels for each of the irreducible representations (1)Σg (+), (1)Σu (+), (1)Σg (-), and (1)Σu (-), to an estimated accuracy of 0.1 mEh of the exact result in this basis.

The Journal of Chemical Physics•Jan 9, 2015•36 citations

Simulation of Anyons Using Symmetric Tensor Network Algorithms

Robert N. C. Pfeifer

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.

arXiv:1202.1522•Feb 7, 2012•6 citations
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 (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.
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\mathcal{G}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.
SeQuant is an open-source library for symbolic algebra of tensors over commutative (scalar) and non-commutative (operator) rings. The key innovation supporting most of its functionality is a graph-theoretic tensor network (TN) canonicalizer that can handle tensor networks with symmetries faster than their standard group-theoretic counterparts. The TN canonicalizer is used for routine simplification of conventional tensor expressions, for optimizing application of Wick's theorem (used to canonicalize products of tensors over operator fields), and for manipulation of the intermediate representation leading to the numerical evaluation. Notable features of SeQuant include support for noncovariant tensor networks (which often arise from tensor decompositions) and for tensors with modes that depend parametrically on indices of other tensor modes (such dependencies between degrees of freedom are naturally viewed as nesting of tensors, "tensors of tensors" arising in block-wise data compressions in data science and modern quantum simulation). SeQuant blurs the line between pure symbolic manipulation/code generation and numerical evaluation by including compiler-like components to optimize and directly interpret tensor expressions using external numerical tensor algebra frameworks. The SeQuant source code is available at https://github.com/ValeevGroup/SeQuant.
Tensor network decompositions offer an efficient description of certain many-body states of a lattice system and are the basis of a wealth of numerical simulation algorithms. In a recent paper [Phys. Rev. A 82, 050301 (2010)] we discussed how to incorporate a global internal symmetry, given by a compact, completely reducible group G\mathcal{G}G, into tensor network decompositions and algorithms. Here we specialize to the case of Abelian groups and, for concreteness, to a U(1) symmetry, associated, e.g., with particle number conservation. We consider tensor networks made of tensors that are invariant (or covariant) under the symmetry, and explain how to decompose and manipulate such tensors in order to exploit their symmetry. In numerical calculations, the use of U(1)-symmetric tensors allows selection of a specific number of particles, ensures the exact preservation of particle number, and significantly reduces computational costs. We illustrate all these points in the context of the multiscale entanglement renormalization Ansatz.
We extend our previous work [S. Sharma and G. K.-L. Chan, J. Chem. Phys. 136, 124121 (2012)], which described a spin-adapted (SU(2) symmetry) density matrix renormalization group algorithm, to additionally utilize general non-Abelian point group symmetries. A key strength of the present formulation is that the requisite tensor operators are not hard-coded for each symmetry group, but are instead generated on the fly using the appropriate Clebsch-Gordan coefficients. This allows our single implementation to easily enable (or disable) any non-Abelian point group symmetry (including SU(2) spin symmetry). We use our implementation to compute the ground state potential energy curve of the C2 dimer in the cc-pVQZ basis set (with a frozen-core), corresponding to a Hilbert space dimension of 10(12) many-body states. While our calculated energy lies within the 0.3 mEh error bound of previous initiator full configuration interaction quantum Monte Carlo and correlation energy extrapolation by intrinsic scaling calculations, our estimated residual error is only 0.01 mEh, much more accurate than these previous estimates. Due to the additional efficiency afforded by the algorithm, the excitation energies (Te) of eight lowest lying excited states: a(3)Πu, b(3)Σg (-), A(1)Πu, c(3)Σu (+), B(1)Δg, B(') (1)Σg (+), d(3)Πg, and C(1)Πg are calculated, which agree with experimentally derived values to better than 0.06 eV. In addition, we also compute the potential energy curves of twelve states: the three lowest levels for each of the irreducible representations (1)Σg (+), (1)Σu (+), (1)Σg (-), and (1)Σu (-), to an estimated accuracy of 0.1 mEh of the exact result in this basis.
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.