Tensor-network methods have opened a powerful route for the study of dynamical spectral functions in two-dimensional quantum systems. However, existing approaches within the framework of infinite projected entangled-pair states construct the required renormalization tensors solely from the ground-state environment and can suffer from severe numerical instability. We identify the origin of this instability and introduce an excitation-tailored corner-transfer-matrix renormalization-group (ET-CTMRG) method to resolve it. By incorporating excitation tensors into the renormalization procedure, the method constructs a substantially more accurate effective Hamiltonian matrix and thereby yields reliable and well-converged excitation spectra. For Heisenberg antiferromagnets, it reduces truncation errors by orders of magnitude and for the particularly complex case of the supersolid phase in the triangular-lattice XXZ magnet $\mathrm{K_2Co(SeO_3)_2}$, it achieves excellent quantitative agreement with inelastic neutron-scattering measurements. ET-CTMRG therefore provides a robust framework for investigating the dynamical properties of strongly correlated quantum systems.
Optimized Tensor-Network Renormalization for Quantum Dynamics: Resolving the Spectral Function of K2Co(SeO3)2
Tensor-network methods have opened a powerful route for the study of dynamical spectral functions in two-dimensional quantum systems. However, existing approaches within the framework of infinite projected entangled-pair states construct the required renormalization tensors solely from the ground-state environment and can suffer from severe numerical instability. We identify the origin of this instability and introduce an excitation-tailored corner-transfer-matrix renormalization-group (ET-CTMRG) method to resolve it. By incorporating excitation tensors into the renormalization procedure, the method constructs a substantially more accurate effective Hamiltonian matrix and thereby yields reliable and well-converged excitation spectra. For Heisenberg antiferromagnets, it reduces truncation errors by orders of magnitude and for the particularly complex case of the supersolid phase in the triangular-lattice XXZ magnet K2Co(SeO3)2, it achieves excellent quantitative agreement with inelastic neutron-scattering measurements. ET-CTMRG therefore provides a robust framework for investigating the dynamical properties of strongly correlated quantum systems.
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.
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.
Within the Projected Entangled Pair State (PEPS) tensor network formalism, a simple update (SU) method has been used to investigate the time evolution of a two-dimensional U(1) critical spin-1/2 spin liquid under Hamiltonian quench [Phys. Rev. B 106, 195132 (2022)]. Here we introduce two different variational frameworks to describe the time dynamics of SU(2)-symmetric translationally-invariant PEPS, aiming to improve the accuracy. In one approach, after using a Trotter-Suzuki decomposition of the time evolution operator in term of two-site elementary gates, one considers a single bond embedded in an environment approximated by a Corner Transfer Matrix Renormalization Group (CTMRG). A variational update of the two tensors on the bond is performed under the application of the elementary gate and then, after symmetrization of the site tensors, the environment is updated. In the second approach, a cluster optimization is performed on a finite (periodic) cluster, maximizing the overlap of the exact time-evolved state with a symmetric finite-size PEPS ansatz. Observables are then computed on the infinite lattice contracting the infinite-PEPS (iPEPS) by CTMRG. We show that the variational schemes outperform the SU method and remain accurate over a significant time interval before hitting the entanglement barrier. Studying the spectrum of the transfer matrix, we find that the asymptotic correlations are very well preserved under time evolution, including the critical nature of the singlet correlations, as expected from the Lieb-Robinson (LR) bound theorem. Consistently, the system (asymptotic) boundary is found to be described by the same Conformal Field Theory of central charge c = 1 c = 1 during time evolution. We also compute the time-evolution of the short distance spin-spin correlations and estimate the LR velocity.
In the context of tensor network states, we for the first time reformulate the corner transfer matrix renormalization group (CTMRG) method into a variational bilevel optimization algorithm. The solution of the optimization problem corresponds to the fixed-point environment pursued in the conventional CTMRG method, from which the partition function of a classical statistical model, represented by an infinite tensor network, can be efficiently evaluated. The validity of this variational idea is demonstrated by the high-precision calculation of the residual entropy of the dimer model, and is further verified by investigating several typical phase transitions in classical spin models, where the obtained critical points and critical exponents all agree with the best known results in literature. Its extension to three-dimensional tensor networks or quantum lattice models is straightforward, as also discussed briefly.
We revisit the corner transfer matrix renormalization group (CTMRG) method of Nishino and Okunishi for contracting two-dimensional (2D) tensor networks and demonstrate that its performance can be substantially improved by determining the tensors using an eigenvalue solver as opposed to the power method used in CTMRG. We also generalize the variational uniform matrix product state (VUMPS) ansatz for diagonalizing 1D quantum Hamiltonians to the case of 2D transfer matrices and discuss similarities with the corner methods. These two new algorithms will be crucial to improving the performance of variational infinite projected entangled pair state (PEPS) methods.
The study of interacting quantum many-body systems is one of the central challenges in modern physics. Interestingly, by understanding the physical properties of such systems from a microscopic description, one can explain many complex quantum phenomena such as the high-temperature superconductivity, fractional quantum Hall effect, etc. However, due to high complexity of these systems, theoretical investigations such as analytical and perturbation theories are limited to studying only a small subset of problems. In the past few decades, numerical simulation has emerged as an indispensable method for studying many phenomena, and has led to immense success and progress in recent years. The main obstacle in addressing quantum many-body systems is that the number of parameters required to represent a state of the system scales exponentially with the system size. However, one hardly needs the full Hilbert space but only a small subspace to describe physical properties of the system. From the point of view in quantum information theory, quantum entanglement is one of the main factors that is responsible for the complexity of studying quantum many-body systems. In particular, the amount of quantum entanglement in such system quantifies the correlation between constituents and therefore it determines how big the subspace is required to represent the system. In this thesis, we are interested in using tensor network states, a new emerging numerical technique to study dynamical and ground state properties of the system. The key ingredient is that a tensor network state which consists of a fairly small number of parameters can represent very well the complete wave function of a system. In one spatial dimension, Matrix Product State (MPS) has become a standard tool, combined with related algorithms such as Time Evolving Block Decimation algorithm (TEBD) to study dynamical properties of a locally-perturbed quantum lattice system. In this thesis, we introduce a technique called ``infinite boundary conditions'' that allows us to use a finite-size MPS to represent the state of an infinite-size system. As a result, one can simulate dynamics of a system for a longer time with a significant computational cost reduction compared to traditional methods. In two spatial dimensions, the tensor network known as Projected Entangled Pair States (PEPS) is emerging as a useful approach. However algorithms for optimizing the matrix elements of these tensors are complex and not very efficient. In particular, the efficiency of the optimization depends critically on the quality of the so-called ``environment'', which is expensive to compute. We have introduced a new method to accelerate the convergence of calculating the PEPS ground state by recycling the environment tensors for many optimization steps, vastly reducing the number of times that the full environment needs to be calculated.