0000000000003171
AUTHOR
Roman Orus
Simulation of many-qubit quantum computation with matrix product states
Matrix product states provide a natural entanglement basis to represent a quantum register and operate quantum gates on it. This scheme can be materialized to simulate a quantum adiabatic algorithm solving hard instances of a NP-Complete problem. Errors inherent to truncations of the exact action of interacting gates are controlled by the size of the matrices in the representation. The property of finding the right solution for an instance and the expected value of the energy are found to be remarkably robust against these errors. As a symbolic example, we simulate the algorithm solving a 100-qubit hard instance, that is, finding the correct product state out of ~ 10^30 possibilities. Accum…
Entanglement continuous unitary transformations
Continuous unitary transformations are a powerful tool to extract valuable information out of quantum many-body Hamiltonians, in which the so-called flow equation transforms the Hamiltonian to a diagonal or block-diagonal form in second quantization. Yet, one of their main challenges is how to approximate the infinitely-many coupled differential equations that are produced throughout this flow. Here we show that tensor networks offer a natural and non-perturbative truncation scheme in terms of entanglement. The corresponding scheme is called "entanglement-CUT" or eCUT. It can be used to extract the low-energy physics of quantum many-body Hamiltonians, including quasiparticle energy gaps. We…
Topological Minimally Entangled States via Geometric Measure
Here we show how the Minimally Entangled States (MES) of a 2d system with topological order can be identified using the geometric measure of entanglement. We show this by minimizing this measure for the doubled semion, doubled Fibonacci and toric code models on a torus with non-trivial topological partitions. Our calculations are done either quasi-exactly for small system sizes, or using the tensor network approach in [R. Orus, T.-C. Wei, O. Buerschaper, A. Garcia-Saez, arXiv:1406.0585] for large sizes. As a byproduct of our methods, we see that the minimisation of the geometric entanglement can also determine the number of Abelian quasiparticle excitations in a given model. The results in …
Benchmarking global SU(2) symmetry in two-dimensional tensor network algorithms
We implement and benchmark tensor network algorithms with $SU(2)$ symmetry for systems in two spatial dimensions and in the thermodynamic limit. Specifically, we implement $SU(2)$-invariant versions of the infinite projected entangled pair states and infinite projected entangled simplex states methods. Our implementation of $SU(2)$ symmetry follows the formalism based on fusion trees from Schmoll et al. [Ann. Phys. 419, 168232 (2020)]. In order to assess the utility of implementing $SU(2)$ symmetry, the algorithms are benchmarked for three models with different local spin: the spin-1 bilinear-biquadratic model on the square lattice, and the kagome Heisenberg antiferromagnets (KHAFs) for spi…
The classical two-dimensional Heisenberg model revisited: An $SU(2)$-symmetric tensor network study
The classical Heisenberg model in two spatial dimensions constitutes one of the most paradigmatic spin models, taking an important role in statistical and condensed matter physics to understand magnetism. Still, despite its paradigmatic character and the widely accepted ban of a (continuous) spontaneous symmetry breaking, controversies remain whether the model exhibits a phase transition at finite temperature. Importantly, the model can be interpreted as a lattice discretization of the $O(3)$ non-linear sigma model in $1+1$ dimensions, one of the simplest quantum field theories encompassing crucial features of celebrated higher-dimensional ones (like quantum chromodynamics in $3+1$ dimensio…
Topological transitions from multipartite entanglement with tensor networks: a procedure for sharper and faster characterization
Topological order in a 2d quantum matter can be determined by the topological contribution to the entanglement R\'enyi entropies. However, when close to a quantum phase transition, its calculation becomes cumbersome. Here we show how topological phase transitions in 2d systems can be much better assessed by multipartite entanglement, as measured by the topological geometric entanglement of blocks. Specifically, we present an efficient tensor network algorithm based on Projected Entangled Pair States to compute this quantity for a torus partitioned into cylinders, and then use this method to find sharp evidence of topological phase transitions in 2d systems with a string-tension perturbation…
Tensor Network Annealing Algorithm for Two-Dimensional Thermal States
Tensor network methods have become a powerful class of tools to capture strongly correlated matter, but methods to capture the experimentally ubiquitous family of models at finite temperature beyond one spatial dimension are largely lacking. We introduce a tensor network algorithm able to simulate thermal states of two-dimensional quantum lattice systems in the thermodynamic limit. The method develops instances of projected entangled pair states and projected entangled pair operators for this purpose. It is the key feature of this algorithm to resemble the cooling down of the system from an infinite temperature state until it reaches the desired finite-temperature regime. As a benchmark we …
Symmetry-protected intermediate trivial phases in quantum spin chains
Symmetry-protected trivial (SPt) phases of matter are the product-state analogue of symmetry-protected topological (SPT) phases. This means, SPt phases can be adiabatically connected to a product state by some path that preserves the protecting symmetry. Moreover, SPt and SPT phases can be adiabatically connected to each other when interaction terms that break the symmetries protecting the SPT order are added in the Hamiltonian. It is also known that spin-1 SPT phases in quantum spin chains can emerge as effective intermediate phases of spin-2 Hamiltonians. In this paper we show that a similar scenario is also valid for SPt phases. More precisely, we show that for a given spin-2 quantum cha…
Infinite projected entangled pair states algorithm improved: Fast full update and gauge fixing
© 2015 American Physical Society. ©2015 American Physical Society. The infinite projected entangled pair states (iPEPS) algorithm [J. Jordan, Phys. Rev. Lett. 101, 250602 (2008)PRLTAO0031-900710.1103/PhysRevLett.101.250602] has become a useful tool in the calculation of ground-state properties of two-dimensional quantum lattice systems in the thermodynamic limit. Despite its many successful implementations, the method has some limitations in its present formulation which hinder its application to some highly entangled systems. The purpose of this paper is to unravel some of these issues, in turn enhancing the stability and efficiency of iPEPS methods. For this, we first introduce the fast f…
Toward Pricing Financial Derivatives with an IBM Quantum Computer
Pricing interest-rate financial derivatives is a major problem in finance, in which it is crucial to accurately reproduce the time evolution of interest rates. Several stochastic dynamics have been proposed in the literature to model either the instantaneous interest rate or the instantaneous forward rate. A successful approach to model the latter is the celebrated Heath-Jarrow-Morton framework, in which its dynamics is entirely specified by volatility factors. In its multifactor version, this model considers several noisy components to capture at best the dynamics of several time-maturing forward rates. However, as no general analytical solution is available, there is a trade-off between t…
Spin-$\frac{1}{2}$ Heisenberg antiferromagnet on the star lattice: Competing valence-bond-solid phases studied by means of tensor networks
Using the infinite Projected Entangled Pair States (iPEPS) algorithm, we study the ground-state properties of the spin-$1/2$ quantum Heisenberg antiferromagnet on the star lattice in the thermodynamic limit. By analyzing the ground-state energy of the two inequivalent bonds of the lattice in different unit-cell structures, we identify two competing Valence-Bond-Solid (VBS) phases for different antiferromagnetic Heisenberg exchange couplings. More precisely, we observe (i) a VBS state which respects the full symmetries of the Hamiltonian, and (ii) a resonating VBS state which, in contrast to previous predictions, has a six-site unit-cell order and breaks $C_3$ symmetry. We also studied the g…
Fine Grained Tensor Network Methods.
We develop a strategy for tensor network algorithms that allows to deal very efficiently with lattices of high connectivity. The basic idea is to fine-grain the physical degrees of freedom, i.e., decompose them into more fundamental units which, after a suitable coarse-graining, provide the original ones. Thanks to this procedure, the original lattice with high connectivity is transformed by an isometry into a simpler structure, which is easier to simulate via usual tensor network methods. In particular this enables the use of standard schemes to contract infinite 2d tensor networks - such as Corner Transfer Matrix Renormalization schemes - which are more involved on complex lattice structu…
A Practical Introduction to Tensor Networks: Matrix Product States and Projected Entangled Pair States
This is a partly non-technical introduction to selected topics on tensor network methods, based on several lectures and introductory seminars given on the subject. It should be a good place for newcomers to get familiarized with some of the key ideas in the field, specially regarding the numerics. After a very general introduction we motivate the concept of tensor network and provide several examples. We then move on to explain some basics about Matrix Product States (MPS) and Projected Entangled Pair States (PEPS). Selected details on some of the associated numerical methods for 1d and 2d quantum lattice systems are also discussed.
Systematic construction of spin liquids on the square lattice from tensor networks with SU(2) symmetry
We elaborate a simple classification scheme of all rank-5 SU(2)-spin rotational symmetric tensors according to i) the on-site physical spin-$S$, (ii) the local Hilbert space $V^{\otimes 4}$ of the four virtual (composite) spins attached to each site and (iii) the irreducible representations of the $C_{4v}$ point group of the square lattice. We apply our scheme to draw a complete list of all SU(2)-symmetric translationally and rotationally-invariant Projected Entangled Pair States (PEPS) with bond dimension $D\leqslant 6$. All known SU(2)-symmetric PEPS on the square lattice are recovered and simple generalizations are provided in some cases. More generally, to each of our symmetry class can…
The spin-1/2 Kagome XXZ model in a field: competition between lattice nematic and solid orders
We study numerically the spin-1/2 XXZ model in a field on an infinite Kagome lattice. We use different algorithms based on infinite Projected Entangled Pair States (iPEPS) for this, namely: (i) with simplex tensors and 9-site unit cell, and (ii) coarse-graining three spins in the Kagome lattice and mapping it to a square-lattice model with nearest-neighbor interactions, with usual PEPS tensors, 6- and 12-site unit cells. Similarly to our previous calculation at the SU(2)-symmetric point (Heisenberg Hamiltonian), for any anisotropy from the Ising limit to the XY limit, we also observe the emergence of magnetization plateaus as a function of the magnetic field, at $m_z = \frac{1}{3}$ using 6-…
Holographic encoding of universality in corner spectra
In numerical simulations of classical and quantum lattice systems, 2d corner transfer matrices (CTMs) and 3d corner tensors (CTs) are a useful tool to compute approximate contractions of infinite-size tensor networks. In this paper we show how the numerical CTMs and CTs can be used, {\it additionally\/}, to extract universal information from their spectra. We provide examples of this for classical and quantum systems, in 1d, 2d and 3d. Our results provide, in particular, practical evidence for a wide variety of models of the correspondence between $d$-dimensional quantum and $(d+1)$-dimensional classical spin systems. We show also how corner properties can be used to pinpoint quantum phase …
Quantum criticality on a chiral ladder: An SU(2) infinite density matrix renormalization group study
In this paper we study the ground-state properties of a ladder Hamiltonian with chiral $\text{SU}(2)$-invariant spin interactions, a possible first step toward the construction of truly two-dimensional nontrivial systems with chiral properties starting from quasi-one-dimensional ones. Our analysis uses a recent implementation by us of $\text{SU}(2)$ symmetry in tensor network algorithms, specifically for infinite density matrix renormalization group. After a preliminary analysis with Kadanoff coarse graining and exact diagonalization for a small-size system, we discuss its bosonization and recap the continuum limit of the model to show that it corresponds to a conformal field theory, in agr…
Infinite projected entangled-pair state algorithm for ruby and triangle-honeycomb lattices
The infinite Projected Entangled-Pair State (iPEPS) algorithm is one of the most efficient techniques for studying the ground-state properties of two-dimensional quantum lattice Hamiltonians in the thermodynamic limit. Here, we show how the algorithm can be adapted to explore nearest-neighbor local Hamiltonians on the ruby and triangle-honeycomb lattices, using the Corner Transfer Matrix (CTM) renormalization group for 2D tensor network contraction. Additionally, we show how the CTM method can be used to calculate the ground state fidelity per lattice site and the boundary density operator and entanglement entropy (EE) on an infinite cylinder. As a benchmark, we apply the iPEPS method to th…
A universal tensor network algorithm for any infinite lattice
We present a general graph-based Projected Entangled-Pair State (gPEPS) algorithm to approximate ground states of nearest-neighbor local Hamiltonians on any lattice or graph of infinite size. By introducing the structural-matrix which codifies the details of tensor networks on any graphs in any dimension $d$, we are able to produce a code that can be essentially launched to simulate any lattice. We further introduce an optimized algorithm to compute simple tensor updates as well as expectation values and correlators with a mean-field-like effective environments. Though not being variational, this strategy allows to cope with PEPS of very large bond dimension (e.g., $D=100$), and produces re…
Spin-S Kagome quantum antiferromagnets in a field with tensor networks
Spin-$S$ Heisenberg quantum antiferromagnets on the Kagome lattice offer, when placed in a magnetic field, a fantastic playground to observe exotic phases of matter with (magnetic analogs of) superfluid, charge, bond or nematic orders, or a coexistence of several of the latter. In this context, we have obtained the (zero temperature) phase diagrams up to $S=2$ directly in the thermodynamic limit thanks to infinite Projected Entangled Pair States (iPEPS), a tensor network numerical tool. We find incompressible phases characterized by a magnetization plateau vs field and stabilized by spontaneous breaking of point group or lattice translation symmetry(ies). The nature of such phases may be se…
A programming guide for tensor networks with global SU(2) symmetry
Abstract This paper is a manual with tips and tricks for programming tensor network algorithms with global S U ( 2 ) symmetry. We focus on practical details that are many times overlooked when it comes to implementing the basic building blocks of codes, such as useful data structures to store the tensors, practical ways of manipulating them, and adapting typical functions for symmetric tensors. Here we do not restrict ourselves to any specific tensor network method, but keep always in mind that the implementation should scale well for simulations of higher-dimensional systems using, e.g., Projected Entangled Pair States, where tensors with many indices may show up. To this end, the structur…
All spin-1 topological phases in a single spin-2 chain
Here we study the emergence of different Symmetry-Protected Topological (SPT) phases in a spin-2 quantum chain. We consider a Heisenberg-like model with bilinear, biquadratic, bicubic, and biquartic nearest-neighbor interactions, as well as uniaxial anisotropy. We show that this model contains four different effective spin-1 SPT phases, corresponding to different representations of the $(\mathbb{Z}_2 \times \mathbb{Z}_2) + T$ symmetry group, where $\mathbb{Z}_2$ is some $\pi$-rotation in the spin internal space and $T$ is time-reversal. One of these phases is equivalent to the usual spin-1 Haldane phase, while the other three are different but also typical of spin-1 systems. The model also …