Researcher profile

Guifre Vidal

Guifre Vidal contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
11topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

14 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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.