symtenet¶
A Python library for symmetric tensor networks: NumPy-ish API on the surface, category theory under the hood.
tenet gives you block-sparse tensors that carry a symmetry exactly — SU(N), SU(2),
U(1), Z2, fermion parity and their products — through exact recoupling coefficients.
A tensor carries legs, contracts with tenet.einsum and factorizes with tenet.linalg,
the way an ndarray does. Blocks live in NumPy, JAX or PyTorch arrays through
autoray, so the same tensor runs under
jax.jit and jax.grad with its symmetry structure intact. On top of the tensor layer,
tenet.network ships finite DMRG and CTMRG.
Install¶
The core install pulls numpy, autoray, opt-einsum and racah-py, and every
symmetry works on it. Two optional extras:
uv add "symtenet[jax]" # jax>=0.10 — pytrees, jit, grad
uv add "symtenet[torch]" # torch>=2.0 — eager blocks
Quickstart¶
A 20-site spin-1/2 Heisenberg chain, U(1)-graded by 2 S^z, to its ground state:
from tenet.models import spin_half
from tenet.network import MPO, MPS, dmrg_
from tenet.symmetry import U1Sector
site, n = spin_half(), 20
terms = []
for i in range(n - 1):
terms.append((1.0, [(site.ops["Sz"], i), (site.ops["Sz"], i + 1)]))
terms.append((0.5, [(site.ops["S+"], i), (site.ops["S-"], i + 1)]))
terms.append((0.5, [(site.ops["S-"], i), (site.ops["S+"], i + 1)]))
h = MPO.from_terms(n, terms)
psi = MPS.product(site.phys, [U1Sector(1 if i % 2 else -1) for i in range(n)])
out = dmrg_(psi, h, chi=64)
print(out.sweeps, out.energy) # 6 -8.682473334398...
print(out.psi.entanglement_entropy()) # {bond: S}, in nats
The Néel product state's own charges put the run in the S^z_tot = 0 sector, and the
site tensors' invariance keeps it there — no projector, no penalty term.
Getting started reads this
example line by line.
With the jax extra, the same tensors are JAX pytrees: jit, grad and vmap reach
through to the blocks while the symmetry structure stays static.
import jax
import tenet
tenet.enable_jax()
t = out.psi[0].to_backend("jax")
g = jax.jit(jax.grad(lambda x: tenet.norm(x) ** 2))(t)
assert g.legs == t.legs
What it supports¶
- Symmetries. SU(N), SU(2), U(1), Z2, fermion parity (
fZ2) and Deligne products of any of them. Non-Abelian sectors are multiplets with exact Clebsch-Gordan, F- and R-symbols; fermionic wires carry their Koszul signs. - Tensors.
SymmetricTensorover a flat tuple ofLegs, one reduced block per allowed fusion channel.einsum,tensordot,compose,transpose,fuse,repartition,trace, andtenet.linalg'ssvd,qr,lq,eigh,eig,polar,expm,left_null, plus the truncatingsvd_truncated/eigh_truncated. - Algorithms. Finite two-site DMRG (
dmrg_, schedules, noise, excited states) with MPS/MPO containers, environment caches and measurement; CTMRG on the directionalEnvCTMand its C4v specializationEnvCTMc4v, with a differentiable fixed-bond move. - Backends. NumPy, JAX and PyTorch blocks through
autoray.tenet.enable_jax()registersSymmetricTensoras a JAX pytree, sojit,gradandvmapreach the blocks while the structure stays static metadata.
Docs¶
- Getting started — install and the first example.
- User guide — tensors, symmetries, contraction, Hamiltonians, DMRG, truncation, JAX, files.
- Tutorials — DMRG, fermions, SU(2), quantum chemistry, CTMRG, VMC.
- Examples — runnable files with their committed output.
- API reference — every public name.
docs/design.md— the categorical model underneath.REPOSITORY_RULES.md— process rules for contributing.
Citation¶
If you use symtenet in your research, please cite it:
@software{symtenet,
author = {Watanabe, Ryo},
title = {{symtenet}: a {Python} library for symmetric tensor networks
--- a {NumPy}-style {API} on the surface, category theory under
the hood},
url = {https://github.com/Ryo-wtnb11/symtenet},
license = {Apache-2.0},
year = {2026}
}
The same metadata lives in
CITATION.cff —
GitHub's "Cite this repository" button renders it as BibTeX or APA.
License¶
Apache License 2.0 — see LICENSE.