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

arXiv:2311.13433 · SciPost Physics Core

State diagrams to determine tree tensor network operators

Richard M. Milbradt, Qunsheng Huang, Christian B. Mendl

Published Jun 19, 2024• 6 citations
arXivHTMLPDF

Abstract

This work is concerned with tree tensor network operators (TTNOs) for representing quantum Hamiltonians. We first establish a mathematical framework connecting tree topologies with state diagrams. Based on these, we devise an algorithm for constructing a TTNO given a Hamiltonian. The algorithm exploits the tensor product structure of the Hamiltonian to add paths to a state diagram, while combining local operators if possible. We test the capabilities of our algorithm on random Hamiltonians for a given tree structure. Additionally, we construct explicit TTNOs for nearest neighbour interactions on a tree topology. Furthermore, we derive a bound on the bond dimension of tensor operators representing arbitrary interactions on trees. Finally, we consider an open quantum system in the form of a Heisenberg spin chain coupled to bosonic bath sites as a concrete example. We find that tree structures allow for lower bond dimensions of the Hamiltonian tensor network representation compared to a matrix product operator structure. This reduction is large enough to reduce the number of total tensor elements required as soon as the number of baths per spin reaches 3.

DOI
10.21468/scipostphyscore.7.2.036

Related papers

Optimal symbolic construction of matrix product operators and tree tensor network operators

Hazar Çakır, Richard M. Milbradt, Christian B. Mendl

This research introduces an improved framework for constructing matrix product operators and tree tensor network operators (TTNOs), crucial tools in quantum simulations. A given (Hamiltonian) operator typically has a known symbolic “sum of operator strings” form that can be translated into a tensor network structure. Combining the existing bipartite-graph-based approach and a newly introduced symbolic Gaussian elimination preprocessing step, our proposed method improves upon earlier algorithms in cases when Hamiltonian terms share the same prefactors. We test the performance of our method against established ones for benchmarking purposes. Finally, we apply our methodology to the model of a cavity filled with molecules in a solvent. This open quantum system is cast into the hierarchical equation of motion setting to obtain an effective Hamiltonian. Construction of the corresponding TTNO demonstrates a sublinear increase of the maximum bond dimension.

arXiv:2502.18630•Physical review. B./Physical review. B•Jun 10, 2025•9 citations

Low-Rank Tree Tensor Network Operators for Long-Range Pairwise Interactions

Gianluca Ceruti, Daniel Kreßner, Dominik Sulz

Abstract. Compactly representing and efficently applying linear operators are fundamental ingredients in tensor network methods for simulating quantum many-body problems and solving high-dimensional problems in scientific computing. In this work, we study such representations for tree tensor networks, the so-called tree tensor network operators (TTNOs), paying particular attention to Hamiltonian operators that involve long-range pairwise interactions between particles. Generalizing the work by Lin, Tong, and others on matrix product operators, we establish a direct connection between the hierarchical low-rank structure of the interaction matrix and the TTNO property. This connection allows us to arrive at very compact TTNO representations by compressing the interaction matrix into a hierarchically semiseparable matrix. Numerical experiments for different quantum spin systems validate our results and highlight the potential advantages of TTNOs over matrix product operators.

arXiv:2405.09952•SIAM Journal on Scientific Computing•Aug 1, 2025•4 citations

Optimal tree tensor network operators for tensor network simulations: Applications to open quantum systems

Weitang Li, Jiajun Ren, Hengrui Yang +2 more

Tree tensor network states (TTNS) decompose the system wavefunction to the product of low-rank tensors based on the tree topology, serving as the foundation of the multi-layer multi-configuration time-dependent Hartree method. In this work, we present an algorithm that automatically constructs the optimal and exact tree tensor network operators (TTNO) for any sum-of-product symbolic quantum operator. The construction is based on the minimum vertex cover of a bipartite graph. With the optimal TTNO, we simulate open quantum systems, such as spin relaxation dynamics in the spin-boson model and charge transport in molecular junctions. In these simulations, the environment is treated as discrete modes and its wavefunction is evolved on equal footing with the system. We employ the Cole-Davidson spectral density to model the glassy phonon environment and incorporate temperature effects via thermo-field dynamics. Our results show that the computational cost scales linearly with the number of discretized modes, demonstrating the efficiency of our approach.

arXiv:2407.13098•The Journal of Chemical Physics•Aug 6, 2024•20 citations

TTNOpt: Tree tensor network package for high-rank tensor compression

Ryo Watanabe, Hidetaka Manabe, Toshiya Hikihara +1 more

We have developed TTNOpt, a software package that utilizes tree tensor networks (TTNs) for quantum spin systems and high-dimensional data analysis. TTNOpt provides efficient and powerful TTN computations by locally optimizing the network structure, guided by the entanglement pattern of the target tensors. For quantum spin systems, TTNOpt searches for the ground state of Hamiltonians with bilinear spin interactions and magnetic fields, and computes physical properties of these states, including the variational energy, bipartite entanglement entropy (EE), single-site expectation values, and two-site correlation functions. Additionally, TTNOpt can target the lowest-energy state within a specified subspace, provided that the Hamiltonian conserves total magnetization. For high-dimensional data analysis, TTNOpt factorizes complex tensors into TTN states that maximize fidelity to the original tensors by optimizing the tensors and the network. When a TTN is provided as input, TTNOpt reconstructs the network based on the EE without referencing the fidelity of the original state. We present three demonstrations of TTNOpt:(1) Ground-state search for the hierarchical chain model with a system size of $256$. The entanglement patterns of the ground state manifest themselves in a tree structure, and TTNOpt successfully identifies the tree.(2) Factorization of a quantic tensor of the $2^{24}$ dimensions representing a three-variable function where each variant has a weak bit-wise correlation. The optimized TTN shows that its structure isolates the variables from each other.(3) Reconstruction of the matrix product network representing a $16$-variable normal distribution characterized by a tree-like correlation structure. TTNOpt can reveal hidden correlation structures of the covariance matrix.

arXiv:2505.05908•Computer Physics Communications•Oct 30, 2025•2 citations

Tensor operators: Constructions and applications for long-range interaction systems

Florian Fröwis, Volckmar Nebendahl, Wolfgang Dür

We consider the representation of operators in terms of tensor networks and their application to the ground-state approximation and time evolution of systems with long-range interactions. We provide an explicit construction to represent an arbitrary many-body Hamilton operator in terms of a one-dimensional tensor network (i.e., as a matrix product operator). For pairwise interactions, we show that such a representation is always efficient and requires a tensor dimension growing only linearly with the number of particles. For systems obeying certain symmetries or restrictions we find optimal representations with minimal tensor dimension. We discuss the analytic and numerical approximation of operators in terms of low-dimensional tensor operators. We demonstrate applications for time evolution and the ground-state approximation, in particular for long-range interaction with inhomogeneous couplings. The operator representations are also generalized to other geometries such as trees and two-dimensional lattices, where we show how to obtain and use efficient tensor network representations respecting a given geometry.

arXiv:1003.1047•Physical Review A•Jun 25, 2010•85 citations

Tree tensor network state approach for solving hierarchical equations of motion

Yaling Ke

The hierarchical equations of motion (HEOM) method is a numerically exact open quantum system dynamics approach. The method is rooted in an exponential expansion of the bath correlation function, which in essence strategically reshapes a continuous environment into a set of effective bath modes that allow for more efficient cutoff at finite temperatures. Based on this understanding, one can map the HEOM method into a Schrödinger-like equation, with a non-Hermitian super-Hamiltonian for an extended wave function being the tensor product of the central system wave function and the Fock state of these effective bath modes. In this work, we explore the possibility of representing the extended wave function as a tree tensor network state (TTNS) and the super-Hamiltonian as a tree tensor network operator of the same structure as the TTNS, as well as the application of a time propagation algorithm using the time-dependent variational principle. Our benchmark calculations based on the spin-boson model with a slow-relaxing bath show that the proposed HEOM+TTNS approach yields consistent results with those of the conventional HEOM method, while the computation is considerably sped up. In addition, the simulation with a genuine TTNS is four times faster than a one-dimensional matrix product state decomposition scheme.

arXiv:2304.05151•The Journal of Chemical Physics•Jun 1, 2023•46 citations
This work is concerned with tree tensor network operators (TTNOs) for representing quantum Hamiltonians. We first establish a mathematical framework connecting tree topologies with state diagrams. Based on these, we devise an algorithm for constructing a TTNO given a Hamiltonian. The algorithm exploits the tensor product structure of the Hamiltonian to add paths to a state diagram, while combining local operators if possible. We test the capabilities of our algorithm on random Hamiltonians for a given tree structure. Additionally, we construct explicit TTNOs for nearest neighbour interactions on a tree topology. Furthermore, we derive a bound on the bond dimension of tensor operators representing arbitrary interactions on trees. Finally, we consider an open quantum system in the form of a Heisenberg spin chain coupled to bosonic bath sites as a concrete example. We find that tree structures allow for lower bond dimensions of the Hamiltonian tensor network representation compared to a matrix product operator structure. This reduction is large enough to reduce the number of total tensor elements required as soon as the number of baths per spin reaches 3.
This research introduces an improved framework for constructing matrix product operators and tree tensor network operators (TTNOs), crucial tools in quantum simulations. A given (Hamiltonian) operator typically has a known symbolic “sum of operator strings” form that can be translated into a tensor network structure. Combining the existing bipartite-graph-based approach and a newly introduced symbolic Gaussian elimination preprocessing step, our proposed method improves upon earlier algorithms in cases when Hamiltonian terms share the same prefactors. We test the performance of our method against established ones for benchmarking purposes. Finally, we apply our methodology to the model of a cavity filled with molecules in a solvent. This open quantum system is cast into the hierarchical equation of motion setting to obtain an effective Hamiltonian. Construction of the corresponding TTNO demonstrates a sublinear increase of the maximum bond dimension.
Abstract. Compactly representing and efficently applying linear operators are fundamental ingredients in tensor network methods for simulating quantum many-body problems and solving high-dimensional problems in scientific computing. In this work, we study such representations for tree tensor networks, the so-called tree tensor network operators (TTNOs), paying particular attention to Hamiltonian operators that involve long-range pairwise interactions between particles. Generalizing the work by Lin, Tong, and others on matrix product operators, we establish a direct connection between the hierarchical low-rank structure of the interaction matrix and the TTNO property. This connection allows us to arrive at very compact TTNO representations by compressing the interaction matrix into a hierarchically semiseparable matrix. Numerical experiments for different quantum spin systems validate our results and highlight the potential advantages of TTNOs over matrix product operators.
Tree tensor network states (TTNS) decompose the system wavefunction to the product of low-rank tensors based on the tree topology, serving as the foundation of the multi-layer multi-configuration time-dependent Hartree method. In this work, we present an algorithm that automatically constructs the optimal and exact tree tensor network operators (TTNO) for any sum-of-product symbolic quantum operator. The construction is based on the minimum vertex cover of a bipartite graph. With the optimal TTNO, we simulate open quantum systems, such as spin relaxation dynamics in the spin-boson model and charge transport in molecular junctions. In these simulations, the environment is treated as discrete modes and its wavefunction is evolved on equal footing with the system. We employ the Cole-Davidson spectral density to model the glassy phonon environment and incorporate temperature effects via thermo-field dynamics. Our results show that the computational cost scales linearly with the number of discretized modes, demonstrating the efficiency of our approach.
We have developed TTNOpt, a software package that utilizes tree tensor networks (TTNs) for quantum spin systems and high-dimensional data analysis. TTNOpt provides efficient and powerful TTN computations by locally optimizing the network structure, guided by the entanglement pattern of the target tensors. For quantum spin systems, TTNOpt searches for the ground state of Hamiltonians with bilinear spin interactions and magnetic fields, and computes physical properties of these states, including the variational energy, bipartite entanglement entropy (EE), single-site expectation values, and two-site correlation functions. Additionally, TTNOpt can target the lowest-energy state within a specified subspace, provided that the Hamiltonian conserves total magnetization. For high-dimensional data analysis, TTNOpt factorizes complex tensors into TTN states that maximize fidelity to the original tensors by optimizing the tensors and the network. When a TTN is provided as input, TTNOpt reconstructs the network based on the EE without referencing the fidelity of the original state. We present three demonstrations of TTNOpt:(1) Ground-state search for the hierarchical chain model with a system size of 256256256. The entanglement patterns of the ground state manifest themselves in a tree structure, and TTNOpt successfully identifies the tree.(2) Factorization of a quantic tensor of the 2242^{24}224 dimensions representing a three-variable function where each variant has a weak bit-wise correlation. The optimized TTN shows that its structure isolates the variables from each other.(3) Reconstruction of the matrix product network representing a 161616-variable normal distribution characterized by a tree-like correlation structure. TTNOpt can reveal hidden correlation structures of the covariance matrix.
We consider the representation of operators in terms of tensor networks and their application to the ground-state approximation and time evolution of systems with long-range interactions. We provide an explicit construction to represent an arbitrary many-body Hamilton operator in terms of a one-dimensional tensor network (i.e., as a matrix product operator). For pairwise interactions, we show that such a representation is always efficient and requires a tensor dimension growing only linearly with the number of particles. For systems obeying certain symmetries or restrictions we find optimal representations with minimal tensor dimension. We discuss the analytic and numerical approximation of operators in terms of low-dimensional tensor operators. We demonstrate applications for time evolution and the ground-state approximation, in particular for long-range interaction with inhomogeneous couplings. The operator representations are also generalized to other geometries such as trees and two-dimensional lattices, where we show how to obtain and use efficient tensor network representations respecting a given geometry.
The hierarchical equations of motion (HEOM) method is a numerically exact open quantum system dynamics approach. The method is rooted in an exponential expansion of the bath correlation function, which in essence strategically reshapes a continuous environment into a set of effective bath modes that allow for more efficient cutoff at finite temperatures. Based on this understanding, one can map the HEOM method into a Schrödinger-like equation, with a non-Hermitian super-Hamiltonian for an extended wave function being the tensor product of the central system wave function and the Fock state of these effective bath modes. In this work, we explore the possibility of representing the extended wave function as a tree tensor network state (TTNS) and the super-Hamiltonian as a tree tensor network operator of the same structure as the TTNS, as well as the application of a time propagation algorithm using the time-dependent variational principle. Our benchmark calculations based on the spin-boson model with a slow-relaxing bath show that the proposed HEOM+TTNS approach yields consistent results with those of the conventional HEOM method, while the computation is considerably sped up. In addition, the simulation with a genuine TTNS is four times faster than a one-dimensional matrix product state decomposition scheme.