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:2506.19339 · Physical review. B./Physical review. B

Efficient optimization of variational tensor-network approach to three-dimensional statistical systems

Xiaze Xu, Tong-Yu Lin, Guang-Ming Zhang

Published Sep 18, 2025• 2 citations
arXivHTMLPDF

Abstract

Variational tensor network optimization has become a powerful tool for studying classical statistical models in two dimensions. However, its application to three-dimensional systems remains limited, primarily due to the high computational cost associated with evaluating the free-energy density and its gradient. This process requires contracting a triple-layer tensor network composed of a projected entangled pair operator and projected entangled pair states. In this paper, we employ a split corner-transfer renormalization group scheme tailored for the contraction of such a triple-layer network, which reduces the computational complexity while keeping high accuracy. Through numerical benchmarks on the three-dimensional classical Ising model, we demonstrate that the proposed scheme achieves numerical results comparable to the most recent Monte Carlo simulations, providing a substantial speedup over previous variational tensor network approaches. This makes this method well-suited for efficient gradient-based optimization in three-dimensional tensor network simulations.

DOI
10.1103/7dbj-w7cn

Related papers

Variationally optimizing infinite projected entangled-pair states at large bond dimensions: A split corner transfer matrix renormalization group approach

Jan Naumann, Erik Lennart Weerda, Jens Eisert +2 more

Projected entangled-pair states (PEPS) have become a powerful tool for studying quantum many-body systems in the condensed matter and quantum materials context, particularly with advances in variational energy optimization methods. A key challenge within this framework is the computational cost associated with the contraction of the two-dimensional lattice, crucial for calculating state vector norms and expectation values. The conventional approach, using the corner transfer matrix renormalization group (CTMRG), involves combining two tensor network layers, resulting in significant time and memory demands. In this work, we introduce an alternative split-CTMRG algorithm, which maintains separate PEPS layers and leverages modified environment tensors, reducing computational complexity while preserving accuracy. Benchmarks on quantum lattice models demonstrate substantial speedups for variational energy optimization, rendering this method valuable for large-scale PEPS simulations.

arXiv:2502.10298•Physical review. B./Physical review. B•Jun 6, 2025•7 citations

Single-layer tensor network approach for three-dimensional quantum systems

I. V. Lukin, Andrii Sotnikov

Calculation of observables with three-dimensional projected entangled pair states is generally hard, as it requires a contraction of complex multi-layer tensor networks. We utilize the multi-layer structure of these tensor networks to largely simplify the contraction. The proposed approach involves the usage of the layer structure both to simplify the search for the boundary projected entangled pair states and the single-layer mapping of the final corner transfer matrix renormalization group contraction. We benchmark our results on the cubic lattice Heisenberg model, reaching the bond dimension D = 7 , and find a good agreement with the previous results. Published by the American Physical Society 2024

arXiv:2405.01489•Physical review. B./Physical review. B•Aug 19, 2024•6 citations

Simulation of three-dimensional quantum systems with projected entangled-pair states

Patrick C. G. Vlaar, Philippe Corboz

Tensor network algorithms have proven to be very powerful tools for studying one- and two-dimensional quantum many-body systems. However, their application to three-dimensional (3D) quantum systems has so far been limited, mostly because the efficient contraction of a 3D tensor network is very challenging. In this paper, we develop and benchmark two contraction approaches for infinite projected entangled-pair states (iPEPS) in 3D. The first approach is based on a contraction of a finite cluster of tensors including an effective environment to approximate the full 3D network. The second approach performs a full contraction of the network by first iteratively contracting layers of the network with a boundary iPEPS, followed by a contraction of the resulting quasi-2D network using the corner transfer matrix renormalization group. Benchmark data for the Heisenberg and Bose-Hubbard models on the cubic lattice show that the algorithms provide competitive results compared to other approaches, making iPEPS a promising tool to study challenging open problems in 3D.

arXiv:2102.06715•Physical review. B./Physical review. B•May 19, 2021•38 citations

Development of Two-Dimensional Tensor Network Methods for the Variational Simulation of Quantum Many-Body States

Jan Naumann

Strongly correlated quantum systems are one of the most interesting and challenging phenomena in modern condensed matter physics. Unfortunately, the theoretical access even to their low energy properties is limited by the exponential growth of the Hilbert space in a many-body setting. To mitigate this computational barrier, tensor network methods are one of the most competitive approaches to classically simulate complex quantum systems. In this thesis we will discuss recent developments in the field of tensor networks to efficiently and accurately simulate two-dimensional quantum many-body states in the thermodynamic limit. Hereby, we introduce the variational update globally optimizing the wave-function, encoded as tensor network in the form of an infinite projected entangled pair state (iPEPS), in order to obtain ground states of many-body Hamiltonian models. We focus in the discussion on the well-established method to approximately contract this two-dimensional structure, the corner transfer matrix renormalization group (CTMRG), and discuss its native application on square lattice and triangular lattice based models. Furthermore, we report the development of a refined numerical algorithm, the split-CTMRG, to significantly reduce its computational complexity, pushing the limits of feasible simulations within modern hardware capabilities. Finally, the method advancements are applied to a currently widely studied model of frustrated magnetism, the Heisenberg maple-leaf antiferromagnet. We are able to find a large phase of vanishing magnetization in contradiction to previous works and, to the best of our knowledge, the lowest reported variational ground state energy at the isotropic point of the model. As an outlook, we describe a potential adaption of the previously developed methods to variationally simulate two-dimensional many-body thermal states.

•Jan 1, 2026•0 citations

Accelerating two-dimensional tensor network contractions using QR decompositions

Yining Zhang, Qi Yang, Philippe Corboz

Infinite projected entangled-pair states (iPEPS) provide a powerful tool for studying strongly correlated systems directly in the thermodynamic limit. A core component of the algorithm is the approximate contraction of the iPEPS, where the computational bottleneck typically lies in the singular value or eigenvalue decompositions involved in the renormalization step. This is particularly true on GPUs, where tensor contractions are substantially faster than these decompositions. Here we propose a contraction scheme for $C_{4v}$-symmetric tensor networks based on combining the corner transfer matrix renormalization group (CTMRG) with QR-decompositions which are substantially faster, especially on GPUs. Our approach achieves up to two orders of magnitude speedup compared to standard CTMRG without loss of accuracy and yields state-of-the-art results for the Heisenberg and $J_1$-$J_2$ models in less than 1 h on an H100 GPU.

arXiv:2505.00494•Physical review. B./Physical review. B•Apr 3, 2026•2 citations

Accelerating two-dimensional tensor network optimization by preconditioning

Xing-Yu Zhang, Qi Yang, Philippe Corboz +2 more

We revisit gradient-based optimization for infinite projected entangled pair states (iPEPS), a tensor network ansatz for simulating many-body quantum systems. This approach is hindered by two major challenges: the high computational cost of evaluating energies and gradients, and an ill-conditioned optimization landscape that slows convergence. To reduce the number of optimization steps, we introduce an efficient preconditioner derived from the leading term of the metric tensor. We benchmark our method against standard optimization techniques on the Heisenberg and Kitaev models, demonstrating substantial improvements in overall computational efficiency. Our approach is broadly applicable across various contraction schemes, unit cell sizes, and Hamiltonians, highlighting the potential of preconditioned optimization to advance tensor network algorithms for strongly correlated systems.

arXiv:2511.09546•Physical review. B./Physical review. B•Feb 13, 2026•3 citations
Variational tensor network optimization has become a powerful tool for studying classical statistical models in two dimensions. However, its application to three-dimensional systems remains limited, primarily due to the high computational cost associated with evaluating the free-energy density and its gradient. This process requires contracting a triple-layer tensor network composed of a projected entangled pair operator and projected entangled pair states. In this paper, we employ a split corner-transfer renormalization group scheme tailored for the contraction of such a triple-layer network, which reduces the computational complexity while keeping high accuracy. Through numerical benchmarks on the three-dimensional classical Ising model, we demonstrate that the proposed scheme achieves numerical results comparable to the most recent Monte Carlo simulations, providing a substantial speedup over previous variational tensor network approaches. This makes this method well-suited for efficient gradient-based optimization in three-dimensional tensor network simulations.
Projected entangled-pair states (PEPS) have become a powerful tool for studying quantum many-body systems in the condensed matter and quantum materials context, particularly with advances in variational energy optimization methods. A key challenge within this framework is the computational cost associated with the contraction of the two-dimensional lattice, crucial for calculating state vector norms and expectation values. The conventional approach, using the corner transfer matrix renormalization group (CTMRG), involves combining two tensor network layers, resulting in significant time and memory demands. In this work, we introduce an alternative split-CTMRG algorithm, which maintains separate PEPS layers and leverages modified environment tensors, reducing computational complexity while preserving accuracy. Benchmarks on quantum lattice models demonstrate substantial speedups for variational energy optimization, rendering this method valuable for large-scale PEPS simulations.
Calculation of observables with three-dimensional projected entangled pair states is generally hard, as it requires a contraction of complex multi-layer tensor networks. We utilize the multi-layer structure of these tensor networks to largely simplify the contraction. The proposed approach involves the usage of the layer structure both to simplify the search for the boundary projected entangled pair states and the single-layer mapping of the final corner transfer matrix renormalization group contraction. We benchmark our results on the cubic lattice Heisenberg model, reaching the bond dimension D = 7 , and find a good agreement with the previous results. Published by the American Physical Society 2024
Tensor network algorithms have proven to be very powerful tools for studying one- and two-dimensional quantum many-body systems. However, their application to three-dimensional (3D) quantum systems has so far been limited, mostly because the efficient contraction of a 3D tensor network is very challenging. In this paper, we develop and benchmark two contraction approaches for infinite projected entangled-pair states (iPEPS) in 3D. The first approach is based on a contraction of a finite cluster of tensors including an effective environment to approximate the full 3D network. The second approach performs a full contraction of the network by first iteratively contracting layers of the network with a boundary iPEPS, followed by a contraction of the resulting quasi-2D network using the corner transfer matrix renormalization group. Benchmark data for the Heisenberg and Bose-Hubbard models on the cubic lattice show that the algorithms provide competitive results compared to other approaches, making iPEPS a promising tool to study challenging open problems in 3D.
Strongly correlated quantum systems are one of the most interesting and challenging phenomena in modern condensed matter physics. Unfortunately, the theoretical access even to their low energy properties is limited by the exponential growth of the Hilbert space in a many-body setting. To mitigate this computational barrier, tensor network methods are one of the most competitive approaches to classically simulate complex quantum systems. In this thesis we will discuss recent developments in the field of tensor networks to efficiently and accurately simulate two-dimensional quantum many-body states in the thermodynamic limit. Hereby, we introduce the variational update globally optimizing the wave-function, encoded as tensor network in the form of an infinite projected entangled pair state (iPEPS), in order to obtain ground states of many-body Hamiltonian models. We focus in the discussion on the well-established method to approximately contract this two-dimensional structure, the corner transfer matrix renormalization group (CTMRG), and discuss its native application on square lattice and triangular lattice based models. Furthermore, we report the development of a refined numerical algorithm, the split-CTMRG, to significantly reduce its computational complexity, pushing the limits of feasible simulations within modern hardware capabilities. Finally, the method advancements are applied to a currently widely studied model of frustrated magnetism, the Heisenberg maple-leaf antiferromagnet. We are able to find a large phase of vanishing magnetization in contradiction to previous works and, to the best of our knowledge, the lowest reported variational ground state energy at the isotropic point of the model. As an outlook, we describe a potential adaption of the previously developed methods to variationally simulate two-dimensional many-body thermal states.
Infinite projected entangled-pair states (iPEPS) provide a powerful tool for studying strongly correlated systems directly in the thermodynamic limit. A core component of the algorithm is the approximate contraction of the iPEPS, where the computational bottleneck typically lies in the singular value or eigenvalue decompositions involved in the renormalization step. This is particularly true on GPUs, where tensor contractions are substantially faster than these decompositions. Here we propose a contraction scheme for C4vC_{4v}C4v​-symmetric tensor networks based on combining the corner transfer matrix renormalization group (CTMRG) with QR-decompositions which are substantially faster, especially on GPUs. Our approach achieves up to two orders of magnitude speedup compared to standard CTMRG without loss of accuracy and yields state-of-the-art results for the Heisenberg and J1J_1J1​-J2J_2J2​ models in less than 1 h on an H100 GPU.
We revisit gradient-based optimization for infinite projected entangled pair states (iPEPS), a tensor network ansatz for simulating many-body quantum systems. This approach is hindered by two major challenges: the high computational cost of evaluating energies and gradients, and an ill-conditioned optimization landscape that slows convergence. To reduce the number of optimization steps, we introduce an efficient preconditioner derived from the leading term of the metric tensor. We benchmark our method against standard optimization techniques on the Heisenberg and Kitaev models, demonstrating substantial improvements in overall computational efficiency. Our approach is broadly applicable across various contraction schemes, unit cell sizes, and Hamiltonians, highlighting the potential of preconditioned optimization to advance tensor network algorithms for strongly correlated systems.