BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Guifre Vidal

Guifre Vidal appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.str-elquant-phhep-thcond-mat.stat-mechcond-mat.dis-nncond-mat.mes-hallcond-mat.quant-gashep-latMachine Learningphysics.comp-phDistributed, Parallel, and Cluster Computing

Catalog footprint

What is connected

53works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

53 published item(s)

preprint2022arXiv

Simulation of quantum many-body dynamics with Tensor Processing Units: Floquet prethermalization

Tensor Processing Units (TPUs) are specialized hardware accelerators developed by Google to support large-scale machine-learning tasks, but they can also be leveraged to accelerate and scale other linear-algebra-intensive computations. In this paper we demonstrate the usage of TPUs for massively parallel, classical simulations of quantum many-body dynamics on long timescales. We apply our methods to study the phenomenon of Floquet prethermalization, i.e., exponentially slow heating in quantum spin chains subject to high-frequency periodic driving. We simulate the dynamics of L=34 qubits for over $10^5$ Floquet periods, corresponding to circuits with millions of two-qubit gates. The circuits simulated have no additional symmetries and represent a pure-state evolution in the full $2^L$-dimensional Hilbert space. This is achieved by distributing the computation over 128 TPU cores. On that size TPU cluster, we find speedups in wall-clock runtime of 230x and 15x when compared to reference CPU and single-GPU simulations, respectively, for shorter 30-qubit simulations that can be handled by all three platforms. We study the computational cost of the simulations, as a function of both the number of qubits and the number of TPU cores used, up to our maximum capacity of L=40 qubits, which requires a ``full pod" of 2048 TPU cores with tens of terabytes of memory in total. For these simulations, an 8-TPU-core machine is comparable to a single A100 GPU, and thus the full TPU pod is comparable to a machine with hundreds of GPUs. However, the TPU pod is more energy and cost efficient, and readily accessible (via Google Cloud), unlike such large many-GPU configurations. We also study the accumulation of numerical error as a function of circuit depth in very deep circuits. Our work demonstrates that TPUs can offer significant advantages for state-of-the-art simulations of quantum many-body dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Continuous entanglement renormalization on the circle

The continuous multi-scale entanglement renormalization ansatz (cMERA) is a variational class of states for quantum fields. As originally formulated, the cMERA applies to infinite systems only. In this paper we generalize the cMERA formalism to a finite circle, which we achieve by wrapping the action of the so-called entangler around the circle. This allows us to transform a cMERA on the line into a cMERA on the circle. In addition, in the case of a Gaussian cMERA for non-interacting quantum fields, the method of images allow us to prove the following result: if on the line a cMERA state is a good approximation to a ground state of a local QFT Hamiltonian, then (under mild assumptions on their correlation functions) the resulting cMERA on a circle is also a good approximation to the ground state of the same QFT Hamiltonian on the circle.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Entanglement renormalization for gauge invariant quantum fields

The continuous multiscale entaglement renormalization ansatz (cMERA) [Haegeman et al., Phys. Rev. Lett., 110, 100402 (2013)] is a variational wavefunctional for ground states of quantum field theories. So far, only scalar bosons and fermions have been considered. In this paper we explain how to generalize the cMERA framework to gauge invariant quantum fields. The fundamental difficulty to be addressed is how to make the gauge constraints (local linear constraints in the Hilbert space) compatible with the UV structure of the cMERA wavefunctional (which is generated by a quasi-local entangler). For simplicity, we consider $U(1)$ gauge theory in $d+1$ spacetime dimensions, a non-interacting theory with massless Hamiltonian $H_{U(1)}$ and Gaussian scale invariant ground state $|Ψ_{U(1)}\rangle$. We propose a gauge invariant cMERA wavefunctional $|Ψ^Λ_{U(1)}\rangle$ that, by construction, accurately reproduces the long distance properties of $|Ψ_{U(1)}\rangle$ while remaining somewhat unentangled at short distances. Moreover, $|Ψ^Λ_{U(1)}\rangle$ is the exact ground state of a gauge invariant, local Hamiltonian $H^Λ_{U(1)}$ whose low energy properties coincide with those of $H_{U(1)}$. Our construction also extends the cMERA formalism to massive (non-gauge invariant) vector boson quantum fields.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Large Scale Distributed Linear Algebra With Tensor Processing Units

We have repurposed Google Tensor Processing Units (TPUs), application-specific chips developed for machine learning, into large-scale dense linear algebra supercomputers. The TPUs' fast inter-core interconnects (ICI)s, physically two-dimensional network topology, and high-bandwidth memory (HBM) permit distributed matrix multiplication algorithms to rapidly become computationally bound. In this regime, the matrix-multiply units (MXU)s dominate the runtime, yielding impressive scaling, performance, and raw size: operating in float32 precision, a full 2048-core pod of third generation TPUs can multiply two matrices with linear size $N= 220= 1 048 576$ in about 2 minutes. Via curated algorithms emphasizing large, single-core matrix multiplications, other tasks in dense linear algebra can similarly scale. As examples, we present (i) QR decomposition; (ii) resolution of linear systems; and (iii) the computation of matrix functions by polynomial iteration, demonstrated by the matrix polar factorization.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Multi-boundary generalization of thermofield double states and their realization in critical quantum spin chains

We propose a multi-boundary generalization of thermofield double states (TFD) of a two-dimensional conformal field theory (CFT) and show, through a conformal map to the complex plane, that they are closely related to multi-point correlation functions. We then also describe how to approximately realize these multi-boundary TFD states numerically on the lattice, starting from a critical quantum spin chain Hamiltonian. In addition, finite size corrections on the lattice are seen to be significantly reduced by the use of \textit{smoothers} -- numerically optimized unitary transformations that locally re-arrange the quantum spin degrees of freedom. One merit of the spin chain realization is that it allows us to probe the properties of the proposed multi-boundary TFD states through numerical experiments, including the characterization of their entanglement structure. As an illustration, we explicitly construct generalized TFD states with three and four boundaries for the Ising CFT and compute entanglement quantities using novel free fermion techniques. We find ranges of parameters where their multipartite entanglement is significant or negligible.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Anomaly Detection with Tensor Networks

Originating from condensed matter physics, tensor networks are compact representations of high-dimensional tensors. In this paper, the prowess of tensor networks is demonstrated on the particular task of one-class anomaly detection. We exploit the memory and computational efficiency of tensor networks to learn a linear transformation over a space with dimension exponential in the number of original features. The linearity of our model enables us to ensure a tight fit around training instances by penalizing the model's global tendency to a predict normality via its Frobenius norm---a task that is infeasible for most deep learning models. Our method outperforms deep and classical algorithms on tabular datasets and produces competitive results on image datasets, despite not exploiting the locality of images.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Determining topological order from infinite projected entangled pair states

We present a method of extracting information about the topological order from the ground state of a strongly correlated two-dimensional system computed with the infinite projected entangled pair state (iPEPS). For topologically ordered systems, the iPEPS wrapped on a torus becomes a superposition of degenerate, locally indistinguishable ground states. Projectors in the form of infinite matrix product operators (iMPO) onto states with well-defined anyon flux are used to compute topological $S$ and $T$ matrices (encoding mutual- and self-statistics of emergent anyons). The algorithm is shown to be robust against a perturbation driving string-net toric code across a phase transition to a ferromagnetic phase. Our approach provides accurate results near quantum phase transition, where the correlation length is prohibitively large for other numerical methods. Moreover, we used numerically optimized iPEPS describing the ground state of the Kitaev honeycomb model in the toric code phase and obtained topological data in excellent agreement with theoretical prediction.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Entanglement and Tensor Networks for Supervised Image Classification

Tensor networks, originally designed to address computational problems in quantum many-body physics, have recently been applied to machine learning tasks. However, compared to quantum physics, where the reasons for the success of tensor network approaches over the last 30 years is well understood, very little is yet known about why these techniques work for machine learning. The goal of this paper is to investigate entanglement properties of tensor network models in a current machine learning application, in order to uncover general principles that may guide future developments. We revisit the use of tensor networks for supervised image classification using the MNIST data set of handwritten digits, as pioneered by Stoudenmire and Schwab [Adv. in Neur. Inform. Proc. Sys. 29, 4799 (2016)]. Firstly we hypothesize about which state the tensor network might be learning during training. For that purpose, we propose a plausible candidate state $|Σ_{\ell}\rangle$ (built as a superposition of product states corresponding to images in the training set) and investigate its entanglement properties. We conclude that $|Σ_{\ell}\rangle$ is so robustly entangled that it cannot be approximated by the tensor network used in that work, which must therefore be representing a very different state. Secondly, we use tensor networks with a block product structure, in which entanglement is restricted within small blocks of $n \times n$ pixels/qubits. We find that these states are extremely expressive (e.g. training accuracy of $99.97 \%$ already for $n=2$), suggesting that long-range entanglement may not be essential for image classification. However, in our current implementation, optimization leads to over-fitting, resulting in test accuracies that are not competitive with other current approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Conformal fields and operator product expansion in critical quantum spin chains

We propose a variational method for identifying lattice operators in a critical quantum spin chain with scaling operators in the underlying conformal field theory (CFT). In particular, this allows us to build a lattice version of the primary operators of the CFT, from which we can numerically estimate the operator product expansion coefficients $C_{αβγ}^{\textrm{ CFT}}$. We demonstrate the approach with the critical Ising quantum spin chain.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Emergence of conformal symmetry in quantum spin chains: anti-periodic boundary conditions and supersymmetry

Universal properties of a critical quantum spin chain are encoded in the underlying conformal field theory (CFT). This underlying CFT is fully characterized by its conformal data. We propose a method to extract the conformal data from a critical quantum spin chain with both periodic and anti-periodic boundary conditions (PBC and APBC) based on low-energy eigenstates, generalizing previous work on spin chains with only PBC. First, scaling dimensions and conformal spins are extracted from the energies and momenta of the eigenstates. Second, the Koo-Saleur formula of lattice Virasoro generators is generalized to APBC and used to identify conformal towers. Third, local operators and string operators on the lattice are identified with CFT operators with PBC and APBC, respectively. Finally, operator product expansion coefficients are extracted by computing matrix elements of lattice primary operators in the low-energy subspaces with PBC and APBC. To go beyond exact diagonalization, tensor network methods based on periodic uniform matrix product states are used. We illustrate our approach with critical and tricritical Ising quantum spin chains. In the latter case, we propose lattice operators that correspond to supervirasoro generators and verify their action on low-energy eigenstates. In this way we explore the emergence of superconformal symmetry in the quantum spin chain.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Explicit construction of local conserved operators in disordered many-body systems

The presence and character of local integrals of motion -- quasi-local operators that commute with the Hamiltonian -- encode valuable information about the dynamics of a quantum system. In particular, strongly disordered many-body systems can generically avoid thermalisation when there are extensively many such operators. In this work, we explicitly construct local conserved operators in $1$D spin chains by directly minimising their commutator with the Hamiltonian. We demonstrate the existence of an extensively large set of local integrals of motion in the many-body localised phase of the disordered XXZ spin chain. These operators are shown to have exponentially decaying tails, in contrast to the ergodic phase where the decay is (at best) polynomial in the size of the subsystem. We study the algebraic properties of localised operators, and confirm that in the many-body localised phase they are well-described by "dressed" spin operators.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Topological conformal defects with tensor networks

The critical 2d classical Ising model on the square lattice has two topological conformal defects: the $\mathbb{Z}_2$ symmetry defect $D_ε$ and the Kramers-Wannier duality defect $D_σ$. These two defects implement antiperiodic boundary conditions and a more exotic form of twisted boundary conditions, respectively. On the torus, the partition function $Z_{D}$ of the critical Ising model in the presence of a topological conformal defect $D$ is expressed in terms of the scaling dimensions $Δ_α$ and conformal spins $s_α$ of a distinct set of primary fields (and their descendants, or conformal towers) of the Ising CFT. This characteristic conformal data $\{Δ_α, s_α\}_{D}$ can be extracted from the eigenvalue spectrum of a transfer matrix $M_{D}$ for the partition function $Z_D$. In this paper we investigate the use of tensor network techniques to both represent and coarse-grain the partition functions $Z_{D_ε}$ and $Z_{D_σ}$ of the critical Ising model with either a symmetry defect $D_ε$ or a duality defect $D_σ$. We also explain how to coarse-grain the corresponding transfer matrices $M_{D_ε}$ and $M_{D_σ}$, from which we can extract accurate numerical estimates of $\{Δ_α, s_α\}_{D_ε}$ and $\{Δ_α, s_α\}_{D_σ}$. Two key new ingredients of our approach are (i) coarse-graining of the defect $D$, which applies to any (i.e. not just topological) conformal defect and yields a set of associated scaling dimensions $Δ_α$, and (ii) construction and coarse-graining of a generalized translation operator using a local unitary transformation that moves the defect, which only exist for topological conformal defects and yields the corresponding conformal spins $s_α$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A tensor network quotient takes the vacuum to the thermal state

In 1+1-dimensional conformal field theory, the thermal state on a circle is related to a certain quotient of the vacuum on a line. We explain how to take this quotient in the MERA tensor network representation of the vacuum and confirm the validity of the construction in the critical Ising model. This result suggests that the tensors comprising MERA can be interpreted as performing local scale transformations, so that adding or removing them emulates conformal maps. In this sense, the optimized MERA recovers local conformal invariance, which is explicitly broken by the choice of lattice. Our discussion also informs the dialogue between tensor networks and holographic duality.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A theory of minimal updates in holography

Consider two quantum critical Hamiltonians $H$ and $\tilde{H}$ on a $d$-dimensional lattice that only differ in some region $\mathcal{R}$. We study the relation between holographic representations, obtained through real-space renormalization, of their corresponding ground states $\left.| ψ\right\rangle$ and $\left.| \tildeψ \right\rangle$. We observe that, even though $\left.| ψ\right\rangle$ and $\left.| \tildeψ \right\rangle$ disagree significantly both inside and outside region $\mathcal{R}$, they still admit holographic descriptions that only differ inside the past causal cone $\mathcal{C}(\mathcal{R})$ of region $\mathcal{R}$, where $\mathcal{C}(\mathcal{R})$ is obtained by coarse-graining region $\mathcal{R}$. We argue that this result follows from a notion of directed influence in the renormalization group flow that is closely connected to the success of Wilson's numerical renormalization group for impurity problems. At a practical level, directed influence allows us to exploit translation invariance when describing a homogeneous system with e.g. an impurity, in spite of the fact that the Hamiltonian is no longer invariant under translations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Edge-entanglement spectrum correspondence in a nonchiral topological phase and Kramers-Wannier duality

In a system with chiral topological order, there is a remarkable correspondence between the edge and entanglement spectra: the low-energy spectrum of the system in the presence of a physical edge coincides with the lowest part of the entanglement spectrum (ES) across a virtual cut of the system, up to rescaling and shifting. In this paper, we explore whether the edge-ES correspondence extends to nonchiral topological phases. Specifically, we consider the Wen-plaquette model which has Z_2 topological order. The unperturbed model displays an exact correspondence: both the edge and entanglement spectra within each topological sector a (a = 1,...,4) are flat and equally degenerate. Here, we show, through a detailed microscopic calculation, that in the presence of generic local perturbations: (i) the effective degrees of freedom for both the physical edge and the entanglement cut consist of a spin-1/2 chain, with effective Hamiltonians H_edge^a and H_ent.^a, respectively, both of which have a Z_2 symmetry enforced by the bulk topological order; (ii) there is in general no match between their low energy spectra, that is, there is no edge-ES correspondence. However, if supplement the Z_2 topological order with a global symmetry (translational invariance along the edge/cut), i.e. by considering the Wen-plaquette model as a symmetry enriched topological phase (SET), then there is a finite domain in Hamiltonian space in which both H_edge^a and H_ent.^a realize the critical Ising model, whose low-energy effective theory is the c = 1/2 Ising CFT. This is achieved because the presence of the global symmetry implies that both Hamiltonians, in addition to being Z_2 symmetric, are Kramers-Wannier self-dual. Thus, the bulk topological order and the global translational symmetry of the Wen-plaquette model as a SET imply an edge-ES correspondence at least in some finite domain in Hamiltonian space.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Lieb-Liniger model with exponentially-decaying interactions: a continuous matrix product state study

The Lieb-Liniger model describes one-dimensional bosons interacting through a repulsive contact potential. In this work, we introduce an extended version of this model by replacing the contact potential with a decaying exponential. Using the recently developed continuous matrix product states techniques, we explore the ground state phase diagram of this model by examining the superfluid and density correlation functions. At weak coupling superfluidity governs the ground state, in a similar way as in the Lieb-Liniger model. However, at strong coupling quasi-crystal and super-Tonks-Girardeau regimes are also found, which are not present in the original Lieb-Liniger case. Therefore the presence of the exponentially-decaying potential leads to a superfluid/super-Tonks-Girardeau/quasi-crystal crossover, when tuning the coupling strength from weak to strong interactions. This corresponds to a Luttinger liquid parameter in the range $K \in (0, \infty)$; in contrast with the Lieb-Liniger model, where $K \in [1, \infty)$, and the screened long-range potential, where $K \in (0, 1]$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Local scale transformations on the lattice with tensor network renormalization

Consider the partition function of a classical system in two spatial dimensions, or the Euclidean path integral of a quantum system in two space-time dimensions, both on a lattice. We show that the tensor network renormalization (TNR) algorithm [\emph{G. Evenbly and G. Vidal, Phys. Rev. Lett. 115, 180405}] can be used to implement local scale transformations on these objects, namely a lattice version of conformal maps. Specifically, we explain how to implement the lattice equivalent of the logarithmic conformal map that transforms the Euclidean plane into a cylinder. As an application, and with the 2D critical Ising model as a concrete example, we use this map to build a lattice version of the scaling operators of the underlying conformal field theory, from which one can extract their scaling dimensions and operator product expansion coefficients.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

NCON: A tensor network contractor for MATLAB

A fundamental process in the implementation of any numerical tensor network algorithm is that of contracting a tensor network. In this process, a network made up of multiple tensors connected by summed indices is reduced to a single tensor or a number by evaluating the index sums. This article presents a MATLAB function ncon(), or "Network CONtractor", which accepts as its input a tensor network and a contraction sequence describing how this network may be reduced to a single tensor or number. As its output it returns that single tensor or number. The function ncon() may be obtained by downloading the source of this preprint.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tensor Network Renormalization

We introduce a coarse-graining transformation for tensor networks that can be applied to study both the partition function of a classical statistical system and the Euclidean path integral of a quantum many-body system. The scheme is based upon the insertion of optimized unitary and isometric tensors (disentanglers and isometries) into the tensor network and has, as its key feature, the ability to remove short-range entanglement/correlations at each coarse-graining step. Removal of short-range entanglement results in scale invariance being explicitly recovered at criticality. In this way we obtain a proper renormalization group flow (in the space of tensors), one that in particular (i) is computationally sustainable, even for critical systems, and (ii) has the correct structure of fixed points, both at criticality and away from it. We demonstrate the proposed approach in the context of the 2D classical Ising model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tensor network renormalization yields the multi-scale entanglement renormalization ansatz

We show how to build a multi-scale entanglement renormalization ansatz (MERA) representation of the ground state of a many-body Hamiltonian $H$ by applying the recently proposed \textit{tensor network renormalization} (TNR) [G. Evenbly and G. Vidal, arXiv:1412.0732] to the Euclidean time evolution operator $e^{-βH}$ for infinite $β$. This approach bypasses the costly energy minimization of previous MERA algorithms and, when applied to finite inverse temperature $β$, produces a MERA representation of a thermal Gibbs state. Our construction endows TNR with a renormalization group flow in the space of wave-functions and Hamiltonians (and not just in the more abstract space of tensors) and extends the MERA formalism to classical statistical systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Constructing local integrals of motion in the many-body localized phase

Many-body localization provides a generic mechanism of ergodicity breaking in quantum systems. In contrast to conventional ergodic systems, many-body localized (MBL) systems are characterized by extensively many local integrals of motion (LIOM), which underlie the absence of transport and thermalization in these systems. Here we report a physically motivated construction of local integrals of motion in the MBL phase. We show that any local operator (e.g., a local particle number or a spin flip operator), evolved with the system's Hamiltonian and averaged over time, becomes a LIOM in the MBL phase. Such operators have a clear physical meaning, describing the response of the MBL system to a local perturbation. In particular, when a local operator represents a density of some globally conserved quantity, the corresponding LIOM describes how this conserved quantity propagates through the MBL phase. Being uniquely defined and experimentally measurable, these LIOMs provide a natural tool for characterizing the properties of the MBL phase, both in experiments and numerical simulations. We demonstrate the latter by numerically constructing an extensive set of LIOMs in the MBL phase of a disordered spin chain model. We show that the resulting LIOMs are quasi-local, and use their decay to extract the localization length and establish the location of the transition between the MBL and ergodic phases.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Disentangling Theorem and Monogamy for Entanglement Negativity

Entanglement negativity is a measure of mixed-state entanglement increasingly used to investigate and characterize emerging quantum many-body phenomena, including quantum criticality and topological order. We present two results for the entanglement negativity: a disentangling theorem, which allows the use of this entanglement measure as a means to detect whether a wave-function of three subsystems $A$, $B$, and $C$ factorizes into a product state for parts $AB_1$ and $B_2C$; and a monogamy relation, which states that if $A$ is very entangled with $B$, then $A$ cannot be simultaneaously very entangled also with $C$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Entanglement contour

In the context of characterizing the structure of quantum entanglement in many-body systems, we introduce the entanglement contour, a tool to identify which real-space degrees of freedom contribute, and how much, to the entanglement of a region A with the rest of the system B. The entanglement contour provides a complementary, more refined approach to characterizing entanglement than just considering the entanglement entropy between A and B, with several concrete advantages. We illustrate this in the context of ground states and quantum quenches in fermionic quadratic systems. For instance, in a quantum critical system in $D = 1$ spatial dimensions, the entanglement contour allows us to determine the central charge of the underlying conformal field theory from just a single partition of the system into regions A and B, (using the entanglement entropy for the same task requires considering several partitions). In $D \geq 2$ dimensions, the entanglement contour can distinguish between gapped and gapless phases that obey a same boundary law for entanglement entropy. During a local or global quantum quench, the time-dependent contour provides a detailed account of the dynamics of entanglement, including propagating entanglement waves, which offers a microscopic explanation of the behavior of the entanglement entropy as a function of time.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Algorithms for entanglement renormalization: boundaries, impurities and interfaces

We propose algorithms, based on the multi-scale entanglement renormalization ansatz, to obtain the ground state of quantum critical systems in the presence of boundaries, impurities, or interfaces. By exploiting the theory of minimal updates [G. Evenbly and G. Vidal, arXiv:1307.0831], the ground state is completely characterized in terms of a number of variational parameters that is independent of the system size, even though the presence of a boundary, an impurity, or an interface explicitly breaks the translation invariance of the host system. Similarly, computational costs do not scale with the system size, allowing the thermodynamic limit to be studied directly and thus avoiding finite size effects e.g. when extracting the universal properties of the critical system.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Entanglement negativity and topological order

We use the entanglement negativity, a measure of entanglement for mixed states, to probe the structure of entanglement in the ground state of a topologically ordered system. Through analytical calculations of the negativity in the ground state(s) of the toric code model, we explicitly show that the entanglement of a region $A$ and its complement $B$ is the sum of two types of contributions. The first type of contributions consists of \textit{boundary entanglement}, which we see to be insensitive to tracing out the interior of $A$ and $B$. It therefore entangles only degrees of freedom in $A$ and $B$ that are close to their common boundary. As it is well-known, each boundary contribution is proportional to the size of the relevant boundary separating $A$ and $B$ and it includes an additive, universal correction. The second contribution appears only when $A$ and $B$ are non-contractible regions (e.g. on a torus) and it consists of long-range entanglement, which we see to be destroyed when tracing out a non-contractible region in the interior of $A$ or $B$. Only the long-range contribution to the entanglement may depend on the specific ground state under consideration.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Global symmetries in tensor network states: symmetric tensors versus minimal bond dimension

Tensor networks offer a variational formalism to efficiently represent wave-functions of extended quantum many-body systems on a lattice. In a tensor network N, the dimension χof the bond indices that connect its tensors controls the number of variational parameters and associated computational costs. In the absence of any symmetry, the minimal bond dimension χ^{min} required to represent a given many-body wave-function |Ψ> leads to the most compact, computationally efficient tensor network description of |Ψ>. In the presence of a global, on-site symmetry, one can use a tensor network N_{sym} made of symmetric tensors. Symmetric tensors allow to exactly preserve the symmetry and to target specific quantum numbers, while their sparse structure leads to a compact description and lowers computational costs. In this paper we explore the trade-off between using a tensor network N with minimal bond dimension χ^{min} and a tensor network N_{sym} made of symmetric tensors, where the minimal bond dimension χ^{min}_{sym} might be larger than χ^{min}. We present two technical results. First, we show that in a tree tensor network, which is the most general tensor network without loops, the minimal bond dimension can always be achieved with symmetric tensors, so that χ^{min}_{sym} = χ^{min}. Second, we provide explicit examples of tensor networks with loops where replacing tensors with symmetric ones necessarily increases the bond dimension, so that χ_{sym}^{min} > χ^{min}. We further argue, however, that in some situations there are important conceptual reasons to prefer a tensor network representation with symmetric tensors (and possibly larger bond dimension) over one with minimal bond dimension.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Matrix product states for anyonic systems and efficient simulation of dynamics

Matrix product states (MPS) have proven to be a very successful tool to study lattice systems with local degrees of freedom such as spins or bosons. Topologically ordered systems can support anyonic particles which are labeled by conserved topological charges and collectively carry non-local degrees of freedom. In this paper we extend the formalism of MPS to lattice systems of anyons. The anyonic MPS is constructed from tensors that explicitly conserve topological charge. We describe how to adapt the time-evolving block decimation (TEBD) algorithm to the anyonic MPS in order to simulate dynamics under a local and charge-conserving Hamiltonian. To demonstrate the effectiveness of anyonic TEBD algorithm, we used it to simulate (i) the ground state (using imaginary time evolution) of an infinite 1D critical system of (a) Ising anyons and (b) Fibonacci anyons both of which are well studied, and (ii) the real time dynamics of an anyonic Hubbard-like model of a single Ising anyon hopping on a ladder geometry with an anyonic flux threading each island of the ladder. Our results pertaining to (ii) give insight into the transport properties of anyons. The anyonic MPS formalism can be readily adapted to study systems with conserved symmetry charges, as this is equivalent to a specialization of the more general anyonic case.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quantum Criticality with the Multi-scale Entanglement Renormalization Ansatz

The goal of this manuscript is to provide an introduction to the multi-scale entanglement renormalization ansatz (MERA) and its application to the study of quantum critical systems. Only systems in one spatial dimension are considered. The MERA, in its scale-invariant form, is seen to offer direct numerical access to the scale-invariant operators of a critical theory. As a result, given a critical Hamiltonian on the lattice, the scale-invariant MERA can be used to characterize the underlying conformal field theory. The performance of the MERA is benchmarked for several critical quantum spin chains, namely Ising, Potts, XX and (modified) Heisenberg models, and an insightful comparison with results obtained using a matrix product state is made. The extraction of accurate conformal data, such as scaling dimensions and operator product expansion coefficients of both local and non-local primary fields, is also illustrated.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Scaling of entanglement entropy in the (branching) multi-scale entanglement renormalization ansatz

We investigate the scaling of entanglement entropy in both the multi-scale entanglement renormalization ansatz (MERA) and in its generalization, the branching MERA. We provide analytical upper bounds for this scaling, which take the general form of a boundary law with various types of multiplicative corrections, including power-law corrections all the way to a bulk law. For several cases of interest, we also provide numerical results that indicate that these upper bounds are saturated to leading order. In particular we establish that, by a suitable choice of holographic tree, the branching MERA can reproduce the logarithmic multiplicative correction of the boundary law observed in Fermi liquids and spin-Bose metals in $D\geq 2$ dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Symmetry protected entanglement renormalization

Entanglement renormalization is a real-space renormalization group (RG) transformation for quantum many-body systems. It generates the multi-scale entanglement renormalization ansatz (MERA), a tensor network capable of efficiently describing a large class of many-body ground states, including those of systems at a quantum critical point or with topological order. The MERA has also been proposed to be a discrete realization of the holographic principle of string theory. In this paper we propose the use of symmetric tensors as a mechanism to build a symmetry protected RG flow, and discuss two important applications of this construction. First, we argue that symmetry protected entanglement renormalization produces the proper structure of RG fixed-points, namely a fixed-point for each symmetry protected phase. Second, in the context of holography, we show that by using symmetric tensors, a global symmetry at the boundary becomes a local symmetry in the bulk, thus explicitly realizing in the MERA a characteristic feature of the AdS/CFT correspondence.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A class of highly entangled many-body states that can be efficiently simulated

We describe a quantum circuit that produces a highly entangled state of N qubits from which one can efficiently compute expectation values of local observables. This construction yields a variational ansatz for quantum many-body states that can be regarded as a generalization of the multi-scale entanglement renormalization ansatz (MERA), and to which we refer as the branching MERA. In a lattice system in D dimensions, the scaling of entanglement of a region of size L^D in the branching MERA is not subject to restrictions such as a boundary law L^{D-1}, but can be proportional to the size of the region, as we demonstrate numerically for D=1,2 dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A real space decoupling transformation for quantum many-body systems

We propose a real space renormalization group method to explicitly decouple into independent components a many-body system that, as in the phenomenon of spin-charge separation, exhibits separation of degrees of freedom at low energies. Our approach produces a branching holographic description of such systems that opens the path to the efficient simulation of the most entangled phases of quantum matter, such as those whose ground state violates a boundary law for entanglement entropy. As in the coarse-graining transformation of [Phys. Rev. Lett. 99, 220405 (2007)], the key ingredient of this decoupling transformation is the concept of entanglement renormalization, or removal of short-range entanglement. We demonstrate the feasibility of the approach, both analytically and numerically, by decoupling in real space the ground state of a critical quantum spin chain into two. Generalized notions of RG flow and of scale invariance are also put forward.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Characterizing topological order by studying the ground states of an infinite cylinder

Given a microscopic lattice Hamiltonian for a topologically ordered phase, we describe a tensor network approach to characterize its emergent anyon model and, in a chiral phase, also its gapless edge theory. First, a tensor network representation of a complete, orthonormal set of ground states on a cylinder of infinite length and finite width is obtained through numerical optimization. Each of these ground states is argued to have a different anyonic flux threading through the cylinder. In a chiral phase, the entanglement spectrum of each ground state is seen to reveal a different sector of the corresponding gapless edge theory. A quasi-orthogonal basis on the torus is then produced by chopping off and reconnecting the tensor network representation on the cylinder. Elaborating on the recent proposal of [Y. Zhang et al. Phys. Rev. B 85, 235151 (2012)], a rotation on the torus yields an alternative basis of ground states and, through the computation of overlaps between bases, the modular matrices S and U (containing the mutual and self statistics of the different anyon species) are extracted. As an application, we study the hard-core boson Haldane model by using the two-dimensional density matrix renormalization group. A thorough characterization of the universal properties of this lattice model, both in the bulk and at the edge, unambiguously shows that its ground space realizes the ν=1/2 bosonic Laughlin state.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Infinite boundary conditions for matrix product state calculations

We propose a formalism to study dynamical properties of a quantum many-body system in the thermodynamic limit by studying a finite system with infinite boundary conditions (IBC) where both finite size effects and boundary effects have been eliminated. For one-dimensional systems, infinite boundary conditions are obtained by attaching two boundary sites to a finite system, where each of these two sites effectively represents a semi-infinite extension of the system. One can then use standard finite-size matrix product state techniques to study a region of the system while avoiding many of the complications normally associated with finite-size calculations such as boundary Friedel oscillations. We illustrate the technique with an example of time evolution of a local perturbation applied to an infinite (translationally invariant) ground state, and use this to calculate the spectral function of the S=1 Heisenberg spin chain. This approach is more efficient and more accurate than conventional simulations based on finite-size matrix product state and density-matrix renormalization-group approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Perfect Sampling with Unitary Tensor Networks

Tensor network states are powerful variational ansätze for many-body ground states of quantum lattice models. The use of Monte Carlo sampling techniques in tensor network approaches significantly reduces the cost of tensor contractions, potentially leading to a substantial increase in computational efficiency. Previous proposals are based on a Markov chain Monte Carlo scheme generated by locally updating configurations and, as such, must deal with equilibration and autocorrelation times, which result in a reduction of efficiency. Here we propose a perfect sampling scheme, with vanishing equilibration and autocorrelation times, for unitary tensor networks -- namely tensor networks based on efficiently contractible, unitary quantum circuits, such as unitary versions of the matrix product state (MPS) and tree tensor network (TTN), and the multi-scale entanglement renormalization ansatz (MERA). Configurations are directly sampled according to their probabilities in the wavefunction, without resorting to a Markov chain process. We also describe a partial sampling scheme that can result in a dramatic (basis-dependent) reduction of sampling error.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Tensor network states and algorithms in the presence of a global SU(2) symmetry

The benefits of exploiting the presence of symmetries in tensor network algorithms have been extensively demonstrated in the context of matrix product states (MPSs). These include the ability to select a specific symmetry sector (e.g. with a given particle number or spin), to ensure the exact preservation of total charge, and to significantly reduce computational costs. Compared to the case of a generic tensor network, the practical implementation of symmetries in the MPS is simplified by the fact that tensors only have three indices (they are trivalent, just as the Clebsch-Gordan coefficients of the symmetry group) and are organized as a one-dimensional array of tensors, without closed loops. Instead, a more complex tensor network, one where tensors have a larger number of indices and/or a more elaborate network structure, requires a more general treatment. In two recent papers, namely (i) [Phys. Rev. A 82, 050301 (2010)] and (ii) [Phys. Rev. B 83, 115125 (2011)], we described how to incorporate a global internal symmetry into a generic tensor network algorithm based on decomposing and manipulating tensors that are invariant under the symmetry. In (i) we considered a generic symmetry group G that is compact, completely reducible and multiplicity free, acting as a global internal symmetry. Then in (ii) we described the practical implementation of Abelian group symmetries. In this paper we describe the implementation of non-Abelian group symmetries in great detail and for concreteness consider an SU(2) symmetry. Our formalism can be readily extended to more exotic symmetries associated with conservation of total fermionic or anyonic charge. As a practical demonstration, we describe the SU(2)-invariant version of the multi-scale entanglement renormalization ansatz and apply it to study the low energy spectrum of a quantum spin chain with a global SU(2) symmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Translation invariance, topology, and protection of criticality in chains of interacting anyons

Using finite size scaling arguments, the critical properties of a chain of interacting anyons can be extracted from the low energy spectrum of a finite system. In Phys. Rev. Lett. 98, 160409 (2007), Feiguin et al. showed that an antiferromagnetic (AFM) chain of Fibonacci anyons on a torus is in the same universality class as the tricritical Ising model, and that criticality is protected by a topological symmetry. In the present paper we first review the graphical formalism for the study of anyons on the disc and demonstrate how this formalism may be consistently extended to the study of systems on surfaces of higher genus. We then employ this graphical formalism to study finite rings of interacting anyons on both the disc and the torus, and show that analysis on the disc necessarily yields an energy spectrum which is a subset of that which is obtained on the torus. For a critical Hamiltonian, one may extract from this subset the scaling dimensions of the local scaling operators which respect the topological symmetry of the system. Related considerations are also shown to apply for open chains.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Variational Monte Carlo with the Multi-Scale Entanglement Renormalization Ansatz

Monte Carlo sampling techniques have been proposed as a strategy to reduce the computational cost of contractions in tensor network approaches to solving many-body systems. Here we put forward a variational Monte Carlo approach for the multi-scale entanglement renormalization ansatz (MERA), which is a unitary tensor network. Two major adjustments are required compared to previous proposals with non-unitary tensor networks. First, instead of sampling over configurations of the original lattice, made of L sites, we sample over configurations of an effective lattice, which is made of just log(L) sites. Second, the optimization of unitary tensors must account for their unitary character while being robust to statistical noise, which we accomplish with a modified steepest descent method within the set of unitary tensors. We demonstrate the performance of the variational Monte Carlo MERA approach in the relatively simple context of a finite quantum spin chain at criticality, and discuss future, more challenging applications, including two dimensional systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Entanglement renormalization and gauge symmetry

A lattice gauge theory is described by a redundantly large vector space that is subject to local constraints, and can be regarded as the low energy limit of an extended lattice model with a local symmetry. We propose a numerical coarse-graining scheme to produce low energy, effective descriptions of lattice models with a local symmetry, such that the local symmetry is exactly preserved during coarse-graining. Our approach results in a variational ansatz for the ground state(s) and low energy excitations of such models and, by extension, of lattice gauge theories. This ansatz incorporates the local symmetry in its structure, and exploits it to obtain a significant reduction of computational costs. We test the approach in the context of the toric code with a magnetic field, equivalent to Z2 lattice gauge theory, for lattices with up to 16 x 16 sites (16^2 x 2 = 512 spins) on a torus. We reproduce the well-known ground state phase diagram of the model, consisting of a deconfined and spin polarized phases separated by a continuous quantum phase transition, and obtain accurate estimates of energy gaps, ground state fidelities, Wilson loops, and several other quantities.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Simulation of fermionic lattice models in two dimensions with Projected Entangled-Pair States: Next-nearest neighbor Hamiltonians

In a recent contribution [Phys. Rev. B 81, 165104 (2010)] fermionic Projected Entangled-Pair States (PEPS) were used to approximate the ground state of free and interacting spinless fermion models, as well as the $t$-$J$ model. This paper revisits these three models in the presence of an additional next-nearest hopping amplitude in the Hamiltonian. First we explain how to account for next-nearest neighbor Hamiltonian terms in the context of fermionic PEPS algorithms based on simulating time evolution. Then we present benchmark calculations for the three models of fermions, and compare our results against analytical, mean-field, and variational Monte Carlo results, respectively. Consistent with previous computations restricted to nearest-neighbor Hamiltonians, we systematically obtain more accurate (or better converged) results for gapped phases than for gapless ones.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Simulation of anyons with tensor network algorithms

Interacting systems of anyons pose a unique challenge to condensed matter simulations due to their non-trivial exchange statistics. These systems are of great interest as they have the potential for robust universal quantum computation, but numerical tools for studying them are as yet limited. We show how existing tensor network algorithms may be adapted for use with systems of anyons, and demonstrate this process for the 1-D Multi-scale Entanglement Renormalisation Ansatz (MERA). We apply the MERA to infinite chains of interacting Fibonacci anyons, computing their scaling dimensions and local scaling operators. The scaling dimensions obtained are seen to be in agreement with conformal field theory. The techniques developed are applicable to any tensor network algorithm, and the ability to adapt these ansaetze for use on anyonic systems opens the door for numerical simulation of large systems of free and interacting anyons in one and two dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Simulation of interacting fermions with entanglement renormalization

We propose and test an algorithm to simulate a lattice system of interacting fermions in two spatial dimensions. The approach is an extension of the entanglement renormalization technique [Phys. Rev. Lett. 99, 220405 (2007)] and the related multi-scale entanglement renormalization ansatz. Benchmark calculations for free and interacting fermions on a lattice of $6\times 6$ sites and periodic boundary conditions confirm the validity of this proposal.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Simulation of strongly correlated fermions in two spatial dimensions with fermionic Projected Entangled-Pair States

We explain how to implement, in the context of projected entangled-pair states (PEPS), the general procedure of fermionization of a tensor network introduced in [P. Corboz, G. Vidal, Phys. Rev. B 80, 165129 (2009)]. The resulting fermionic PEPS, similar to previous proposals, can be used to study the ground state of interacting fermions on a two-dimensional lattice. As in the bosonic case, the cost of simulations depends on the amount of entanglement in the ground state and not directly on the strength of interactions. The present formulation of fermionic PEPS leads to a straightforward numerical implementation that allowed us to recycle much of the code for bosonic PEPS. We demonstrate that fermionic PEPS are a useful variational ansatz for interacting fermion systems by computing approximations to the ground state of several models on an infinite lattice. For a model of interacting spinless fermions, ground state energies lower than Hartree-Fock results are obtained, shifting the boundary between the metal and charge-density wave phases. For the t-J model, energies comparable with those of a specialized Gutzwiller-projected ansatz are also obtained.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Tensor network decompositions in the presence of a global symmetry

Tensor network decompositions offer an efficient description of certain many-body states of a lattice system and are the basis of a wealth of numerical simulation algorithms. We discuss how to incorporate a global symmetry, given by a compact, completely reducible group G, in tensor network decompositions and algorithms. This is achieved by considering tensors that are invariant under the action of the group G. Each symmetric tensor decomposes into two types of tensors: degeneracy tensors, containing all the degrees of freedom, and structural tensors, which only depend on the symmetry group. In numerical calculations, the use of symmetric tensors ensures the preservation of the symmetry, allows selection of a specific symmetry sector, and significantly reduces computational costs. On the other hand, the resulting tensor network can be interpreted as a superposition of exponentially many spin networks. Spin networks are used extensively in loop quantum gravity, where they represent states of quantum geometry. Our work highlights their importance also in the context of tensor network algorithms, thus setting the stage for cross-fertilization between these two areas of research.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Tensor network states and algorithms in the presence of a global U(1) symmetry

Tensor network decompositions offer an efficient description of certain many-body states of a lattice system and are the basis of a wealth of numerical simulation algorithms. In a recent paper [arXiv:0907.2994v1] we discussed how to incorporate a global internal symmetry, given by a compact, completely reducible group G, into tensor network decompositions and algorithms. Here we specialize to the case of Abelian groups and, for concreteness, to a U(1) symmetry, often associated with particle number conservation. We consider tensor networks made of tensors that are invariant (or covariant) under the symmetry, and explain how to decompose and manipulate such tensors in order to exploit their symmetry. In numerical calculations, the use of U(1) symmetric tensors allows selection of a specific number of particles, ensures the exact preservation of particle number, and significantly reduces computational costs. We illustrate all these points in the context of the multi-scale entanglement renormalization ansatz.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Entanglement renormalization in two spatial dimensions

We propose and test a scheme for entanglement renormalization capable of addressing large two-dimensional quantum lattice systems. In a translationally invariant system, the cost of simulations grows only as the logarithm of the lattice size; at a quantum critical point, the simulation cost becomes independent of the lattice size and infinite systems can be analysed. We demonstrate the performance of the scheme by investigating the low energy properties of the 2D quantum Ising model on a square lattice of linear size L={6,9,18,54,inf} with periodic boundary conditions. We compute the ground state and evaluate local observables and two-point correlators. We also produce accurate estimates of the critical magnetic field and critical exponent beta. A calculation of the energy gap shows that it scales as 1/L at the critical point.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Fermionic multi-scale entanglement renormalization ansatz

In a recent contribution [arXiv:0904:4151] entanglement renormalization was generalized to fermionic lattice systems in two spatial dimensions. Entanglement renormalization is a real-space coarse-graining transformation for lattice systems that produces a variational ansatz, the multi-scale entanglement renormalization ansatz (MERA), for the ground states of local Hamiltonians. In this paper we describe in detail the fermionic version of the MERA formalism and algorithm. Starting from the bosonic MERA, which can be regarded both as a quantum circuit or in relation to a coarse-graining transformation, we indicate how the scheme needs to be modified to simulate fermions. To confirm the validity of the approach, we present benchmark results for free and interacting fermions on a square lattice with sizes between $6 \times 6$ and $162\times 162$ and with periodic boundary conditions. The present formulation of the approach applies to generic tensor network algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Numerical study of the hard-core Bose-Hubbard Model on an Infinite Square Lattice

We present a study of the hard-core Bose-Hubbard model at zero temperature on an infinite square lattice using the infinite Projected Entangled Pair State algorithm [Jordan et al., Phys. Rev. Lett. 101, 250602 (2008)]. Throughout the whole phase diagram our values for the ground state energy, particle density and condensate fraction accurately reproduce those previously obtained by other methods. We also explore ground state entanglement, compute two-point correlators and conduct a fidelity-based analysis of the phase diagram. Furthermore, for illustrative purposes we simulate the response of the system when a perturbation is suddenly added to the Hamiltonian.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Simulation of two dimensional quantum systems on an infinite lattice revisited: corner transfer matrix for tensor contraction

An extension of the projected entangled-pair states (PEPS) algorithm to infinite systems, known as the iPEPS algorithm, was recently proposed to compute the ground state of quantum systems on an infinite two-dimensional lattice. Here we investigate a modification of the iPEPS algorithm, where the environment is computed using the corner transfer matrix renormalization group (CTMRG) method, instead of using one-dimensional transfer matrix methods as in the original proposal. We describe a variant of the CTMRG that addresses different directions of the lattice independently, and use it combined with imaginary time evolution to compute the ground state of the two-dimensional quantum Ising model. Near criticality, the modified iPEPS algorithm is seen to provide a better estimation of the order parameter and correlators.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Entanglement renormalization

In the context of real-space renormalization group methods, we propose a novel scheme for quantum systems defined on a D-dimensional lattice. It is based on a coarse-graining transformation that attempts to reduce the amount of entanglement of a block of lattice sites before truncating its Hilbert space. Numerical simulations involving the ground state of a 1D system at criticality show that the resulting coarse-grained site requires a Hilbert space dimension that does not grow with successive rescaling transformations. As a result we can address, in a quasi-exact way, tens of thousands of quantum spins with a computational effort that scales logarithmically in the system's size. The calculations unveil that ground state entanglement in extended quantum systems is organized in layers corresponding to different length scales. At a quantum critical point, each rellevant length scale makes an equivalent contribution to the entanglement of a block with the rest of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Efficient classical simulation of slightly entangled quantum computations

We present a scheme to efficiently simulate, with a classical computer, the dynamics of multipartite quantum systems on which the amount of entanglement (or of correlations in the case of mixed-state dynamics) is conveniently restricted. The evolution of a pure state of n qubits can be simulated by using computational resources that grow linearly in n and exponentially in the entanglement. We show that a pure-state quantum computation can only yield an exponential speed-up with respect to classical computations if the entanglement increases with the size n of the computation, and gives a lower bound on the required growth.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Entanglement monotones

In the context of quantifying entanglement we study those functions of a multipartite state which do not increase under the set of local transformations. A mathematical characterization of these monotone magnitudes is presented. They are then related to optimal strategies of conversion of shared states. More detailed results are presented for pure states of bipartite systems. It is show that more than one measure are required simultaneously in order to quantify completely the non-local resources contained in a bipartite pure state, while examining how this fact does not hold in the so-called asymptotic limit. Finally, monotonicity under local transformations is proposed as the only natural requirement for measures of entanglement.

0 citations0 reviewsNo code linkedOpen access