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:2007.03638 · SciPost Physics

Riemannian optimization of isometric tensor networks

Markus Hauru, Maarten Van Damme, Jutho Haegeman

Published Feb 17, 2021• 68 citations
arXivHTMLPDF

Abstract

Several tensor networks are built of isometric tensors, i.e. tensors satisfying W\dagger W = \mathbb{1} W † W = 1 . Prominent examples include matrix product states (MPS) in canonical form, the multiscale entanglement renormalization ansatz (MERA), and quantum circuits in general, such as those needed in state preparation and quantum variational eigensolvers. We show how gradient-based optimization methods on Riemannian manifolds can be used to optimize tensor networks of isometries to represent e.g. ground states of 1D quantum Hamiltonians. We discuss the geometry of Grassmann and Stiefel manifolds, the Riemannian manifolds of isometric tensors, and review how state-of-the-art optimization methods like nonlinear conjugate gradient and quasi-Newton algorithms can be implemented in this context. We apply these methods in the context of infinite MPS and MERA, and show benchmark results in which they outperform the best previously-known optimization methods, which are tailor-made for those specific variational classes. We also provide open-source implementations of our algorithms.

DOI
10.21468/scipostphys.10.2.040

Related papers

Riemannian geometry and automatic differentiation for optimization problems of quantum physics and quantum technologies

Ilia A Luchnikov, Mikhail E Krechetov, Sergey N. Filippov

Abstract Optimization with constraints is a typical problem in quantum physics and quantum information science that becomes especially challenging for high-dimensional systems and complex architectures like tensor networks. Here we use ideas of Riemannian geometry to perform optimization on the manifolds of unitary and isometric matrices as well as the cone of positive-definite matrices. Combining this approach with the up-to-date computational methods of automatic differentiation, we demonstrate the efficacy of the Riemannian optimization in the study of the low-energy spectrum and eigenstates of multipartite Hamiltonians, variational search of a tensor network in the form of the multiscale entanglement-renormalization ansatz, preparation of arbitrary states (including highly entangled ones) in the circuit implementation of quantum computation, decomposition of quantum gates, and tomography of quantum states. Universality of the developed approach together with the provided open source software enable one to apply the Riemannian optimization to complex quantum architectures well beyond the listed problems, for instance, to the optimal control of noisy quantum systems.

arXiv:2007.01287•New Journal of Physics•Jun 14, 2021•47 citations

Optimizing two-dimensional isometric tensor networks with quantum computers

Sebastian Leontica, Alberto Baiardi, Julian Schuhmacher +2 more

We propose a hybrid quantum-classical algorithm for approximating the ground state of two-dimensional quantum systems using an isometric tensor network ansatz, which maps naturally to quantum circuits. Inspired by the density matrix renormalization group, we optimize tensors sequentially by diagonalizing a series of effective Hamiltonians. These are constructed using a tomography-inspired method on a qubit subset whose size depends only on the bond dimension. Our approach leverages quantum computers to enable accurate solutions without relying on approximate contractions, circumventing the exponential complexity faced by classical techniques. We demonstrate our method on the two-dimensional (2D) transverse-field Ising model, achieving ground-state optimization on up to 25 qubits with modest quantum overhead -- significantly less than standard solutions based on variational quantum eigensolvers. Overall, our results offer a path towards scalable variational quantum algorithms in both noisy and fault-tolerant regimes.

quant-pharXiv:2511.13827•Nov 19, 2025•0 citations

Differentiable programming of isometric tensor networks

Chenhua Geng, Hong-Ye Hu, Yijian Zou

Abstract Differentiable programming is a new programming paradigm which enables large scale optimization through automatic calculation of gradients also known as auto-differentiation. This concept emerges from deep learning, and has also been generalized to tensor network optimizations. Here, we extend the differentiable programming to tensor networks with isometric constraints with applications to multiscale entanglement renormalization ansatz (MERA) and tensor network renormalization (TNR). By introducing several gradient-based optimization methods for the isometric tensor network and comparing with Evenbly–Vidal method, we show that auto-differentiation has a better performance for both stability and accuracy. We numerically tested our methods on 1D critical quantum Ising spin chain and 2D classical Ising model. We calculate the ground state energy for the 1D quantum model and internal energy for the classical model, and scaling dimensions of scaling operators and find they all agree with the theory well.

arXiv:2110.03898•Machine Learning Science and Technology•Jan 6, 2022•12 citations

Optimization schemes for unitary tensor-network circuit

Reza Haghshenas

An efficient representation of a quantum circuit is of great importance in achieving a quantum advantage on current noisy intermediate scale quantum (NISQ) devices and the classical simulation of quantum manybody systems.The quantum circuits are playing the key ingredient in the performance of variational quantum algorithms and quantum dynamics in problems of physics and chemistry.In this paper, we study the role of the network structure of a quantum circuit in its performance.We discuss the variational optimization of quantum circuit (a unitary tensor-network circuit) with different network structures.The ansatz is performed based on a generalization of well-developed multiscale entanglement renormalization algorithm and also the conjugategradient method with an effective line search.We present the benchmarking calculations for different network structures by studying the Heisenberg model in a strongly disordered magnetic field and a tensor-network QR decomposition.Our work can contribute to achieve the most out of NISQ hardware and to classically develop isometric tensor network states.

arXiv:2009.02606•Physical Review Research•May 26, 2021•15 citations

Variational Power of Quantum Circuit Tensor Networks

Reza Haghshenas, Johnnie Gray, Andrew C. Potter +1 more

We characterize the variational power of quantum circuit tensor networks in the representation of physical many-body ground states.Such tensor networks are formed by replacing the dense block unitaries and isometries in standard tensor networks by local quantum circuits.We explore both quantum circuit matrix product states and the quantum circuit multiscale entanglement renormalization Ansatz, and introduce an adaptive method to optimize the resulting circuits to high fidelity with more than 10 4 parameters.We benchmark their expressiveness against standard tensor networks, as well as other common circuit architectures, for the 1D and 2D Heisenberg and 1D Fermi-Hubbard models.We find quantum circuit tensor networks to be substantially more expressive than other quantum circuits for these problems, and that they can even be more compact than standard tensor networks.Extrapolating to circuit depths which can no longer be emulated classically, this suggests a region of advantage in quantum expressiveness in the representation of physical ground states.

arXiv:2107.01307•Physical Review X•Mar 11, 2022•80 citations

Combining non-parametric quantum states and MERA tensor networks for ground-state optimization

Julian Schuhmacher, Alberto Baiardi, Francesco Tacchino +1 more

Hybrid tensor networks offer a promising route to enhance the expressivity of classical tensor network methods by incorporating quantum states prepared on a quantum computer. Existing approaches are limited by the variational optimization of the quantum component of the tensor network. In this work, we introduce an alternative strategy that combines a non-parametric quantum state prepared through quantum annealing and a classical isometric tensor network. The latter is variationally optimized while the former is used as a fixed, boundary tensor resource in the form of classical shadows. We demonstrate the feasibility of this approach through extensive numerical simulations on the transverse-field Ising model, showing that the optimization procedure remains robust under statistical and hardware noise. Moreover, our results indicate that our newly proposed setup improves the accuracy of the obtained ground state approximation compared to the original quantum simulation, without increasing the depth of the applied quantum circuits. Therefore, this setup offers a practical route to scale variational quantum algorithms towards the quantum utility scale.

quant-pharXiv:2605.21447•May 21, 2026•0 citations
Several tensor networks are built of isometric tensors, i.e. tensors satisfying W\dagger W = \mathbb{1} W † W = 1 . Prominent examples include matrix product states (MPS) in canonical form, the multiscale entanglement renormalization ansatz (MERA), and quantum circuits in general, such as those needed in state preparation and quantum variational eigensolvers. We show how gradient-based optimization methods on Riemannian manifolds can be used to optimize tensor networks of isometries to represent e.g. ground states of 1D quantum Hamiltonians. We discuss the geometry of Grassmann and Stiefel manifolds, the Riemannian manifolds of isometric tensors, and review how state-of-the-art optimization methods like nonlinear conjugate gradient and quasi-Newton algorithms can be implemented in this context. We apply these methods in the context of infinite MPS and MERA, and show benchmark results in which they outperform the best previously-known optimization methods, which are tailor-made for those specific variational classes. We also provide open-source implementations of our algorithms.
Abstract Optimization with constraints is a typical problem in quantum physics and quantum information science that becomes especially challenging for high-dimensional systems and complex architectures like tensor networks. Here we use ideas of Riemannian geometry to perform optimization on the manifolds of unitary and isometric matrices as well as the cone of positive-definite matrices. Combining this approach with the up-to-date computational methods of automatic differentiation, we demonstrate the efficacy of the Riemannian optimization in the study of the low-energy spectrum and eigenstates of multipartite Hamiltonians, variational search of a tensor network in the form of the multiscale entanglement-renormalization ansatz, preparation of arbitrary states (including highly entangled ones) in the circuit implementation of quantum computation, decomposition of quantum gates, and tomography of quantum states. Universality of the developed approach together with the provided open source software enable one to apply the Riemannian optimization to complex quantum architectures well beyond the listed problems, for instance, to the optimal control of noisy quantum systems.
We propose a hybrid quantum-classical algorithm for approximating the ground state of two-dimensional quantum systems using an isometric tensor network ansatz, which maps naturally to quantum circuits. Inspired by the density matrix renormalization group, we optimize tensors sequentially by diagonalizing a series of effective Hamiltonians. These are constructed using a tomography-inspired method on a qubit subset whose size depends only on the bond dimension. Our approach leverages quantum computers to enable accurate solutions without relying on approximate contractions, circumventing the exponential complexity faced by classical techniques. We demonstrate our method on the two-dimensional (2D) transverse-field Ising model, achieving ground-state optimization on up to 25 qubits with modest quantum overhead -- significantly less than standard solutions based on variational quantum eigensolvers. Overall, our results offer a path towards scalable variational quantum algorithms in both noisy and fault-tolerant regimes.
Abstract Differentiable programming is a new programming paradigm which enables large scale optimization through automatic calculation of gradients also known as auto-differentiation. This concept emerges from deep learning, and has also been generalized to tensor network optimizations. Here, we extend the differentiable programming to tensor networks with isometric constraints with applications to multiscale entanglement renormalization ansatz (MERA) and tensor network renormalization (TNR). By introducing several gradient-based optimization methods for the isometric tensor network and comparing with Evenbly–Vidal method, we show that auto-differentiation has a better performance for both stability and accuracy. We numerically tested our methods on 1D critical quantum Ising spin chain and 2D classical Ising model. We calculate the ground state energy for the 1D quantum model and internal energy for the classical model, and scaling dimensions of scaling operators and find they all agree with the theory well.
An efficient representation of a quantum circuit is of great importance in achieving a quantum advantage on current noisy intermediate scale quantum (NISQ) devices and the classical simulation of quantum manybody systems.The quantum circuits are playing the key ingredient in the performance of variational quantum algorithms and quantum dynamics in problems of physics and chemistry.In this paper, we study the role of the network structure of a quantum circuit in its performance.We discuss the variational optimization of quantum circuit (a unitary tensor-network circuit) with different network structures.The ansatz is performed based on a generalization of well-developed multiscale entanglement renormalization algorithm and also the conjugategradient method with an effective line search.We present the benchmarking calculations for different network structures by studying the Heisenberg model in a strongly disordered magnetic field and a tensor-network QR decomposition.Our work can contribute to achieve the most out of NISQ hardware and to classically develop isometric tensor network states.
We characterize the variational power of quantum circuit tensor networks in the representation of physical many-body ground states.Such tensor networks are formed by replacing the dense block unitaries and isometries in standard tensor networks by local quantum circuits.We explore both quantum circuit matrix product states and the quantum circuit multiscale entanglement renormalization Ansatz, and introduce an adaptive method to optimize the resulting circuits to high fidelity with more than 10 4 parameters.We benchmark their expressiveness against standard tensor networks, as well as other common circuit architectures, for the 1D and 2D Heisenberg and 1D Fermi-Hubbard models.We find quantum circuit tensor networks to be substantially more expressive than other quantum circuits for these problems, and that they can even be more compact than standard tensor networks.Extrapolating to circuit depths which can no longer be emulated classically, this suggests a region of advantage in quantum expressiveness in the representation of physical ground states.
Hybrid tensor networks offer a promising route to enhance the expressivity of classical tensor network methods by incorporating quantum states prepared on a quantum computer. Existing approaches are limited by the variational optimization of the quantum component of the tensor network. In this work, we introduce an alternative strategy that combines a non-parametric quantum state prepared through quantum annealing and a classical isometric tensor network. The latter is variationally optimized while the former is used as a fixed, boundary tensor resource in the form of classical shadows. We demonstrate the feasibility of this approach through extensive numerical simulations on the transverse-field Ising model, showing that the optimization procedure remains robust under statistical and hardware noise. Moreover, our results indicate that our newly proposed setup improves the accuracy of the obtained ground state approximation compared to the original quantum simulation, without increasing the depth of the applied quantum circuits. Therefore, this setup offers a practical route to scale variational quantum algorithms towards the quantum utility scale.