introduction to ipeps
play

Introduction to iPEPS Philippe Corboz, Institute for Theoretical - PowerPoint PPT Presentation

Sep 14, 2023 •553 likes •1.32k views

Introduction to iPEPS Philippe Corboz, Institute for Theoretical Physics, University of Amsterdam Overview: tensor networks in 1D and 2D 1D 1D MERA MPS Matrix-product state Multi-scale entanglement renormalization ansatz and more 1D tree

  1. Introduction to iPEPS Philippe Corboz, Institute for Theoretical Physics, University of Amsterdam

  2. Overview: tensor networks in 1D and 2D 1D 1D MERA MPS Matrix-product state Multi-scale entanglement renormalization ansatz and more ‣ 1D tree tensor 1 2 3 4 5 6 7 8 network Underlying ansatz of the ‣ correlator density-matrix renormalization product states group ( DMRG ) method i 1 i 2 i 3 i 4 i 5 i 6 i 7 i 8 i 9 i 10 i 11 i 12 i 13 i 14 i 15 i 16 i 17 i 18 ‣ ... 2D PEPS (TPS) 2D MERA and more projected entangled-pair state (tensor product state) ‣ Entangled- plaquette states ‣ 2D tree tensor network ‣ String-bond states ‣ ...

  3. Outline of the 2 lectures ‣ Part I: iPEPS ansatz ✦ Repetition: area law of the entanglement entropy ‣ Part II: Contraction of PEPS / iPEPS ✦ MPS-MPO approach, corner-transfer-matrix (CTM) method, Tensor Renormalization Group (TRG), Tensor network renormalization (TNR) ✦ Simple examples to get started: ➡ solving the 2D classical Ising model with the CTM method ➡ simple 2D quantum case (D=2, rotational symmetric) ‣ Interlude: Example application: the Shastry-Sutherland model ‣ Part III: Optimization ‣ Part IV: Computational cost & benchmarks ‣ Part V: iPEPS applications ‣ Outlook & summary

  4. PART I: iPEPS ansatz

  5. “Corner” of the Hilbert space Ground states (local H) ★ GS of local H’s are less entangled than a Hilbert random state in the Hilbert space ★ Area law of the entanglement entropy space

  6. Area law of the entanglement entropy . . . . E . . 1D 2D A E A E . . . . . . L L . . . . . . # relevant states � Entanglement entropy S ( A ) = − tr[ ρ A log ρ A ] = − λ i log λ i χ ∼ exp( S ) i General (random) state Ground state (local Hamiltonian) S ( L ) ∼ L d ( volume ) S ( L ) ∼ L d − 1 ( area law ) Critical ground states: 1D S ( L ) = const χ = const (all in 1D but not all in 2D) χ ∼ exp( α L ) 2D S ( L ) ∼ α L 1D S ( L ) ∼ log( L ) 2D S ( L ) ∼ L log( L )

  7. MPS & PEPS 1D MPS Matrix-product state Bond dimension D 1 2 3 4 5 6 7 8 Physical indices (lattices sites) S. R. White, PRL 69, 2863 (1992) Fannes et al., CMP 144, 443 (1992) Östlund, Rommer, PRL 75, 3537 (1995) ✓ Reproduces area-law in 1D S ( L ) = const

  8. MPS & PEPS 1D MPS Matrix-product state Bond dimension D ➡ One bond can contribute 1 2 3 4 5 6 7 8 at most log(D) to the entanglement entropy A E L S ( A ) ≤ log ( D ) = const rank ( ρ A ) ≤ D ✓ Reproduces area-law in 1D S ( L ) = const

  9. MPS & PEPS 1D 2D MPS can we use an MPS? Matrix-product state Bond dimension D 1 2 3 4 5 6 7 8 L Physical indices (lattices sites) S. R. White, PRL 69, 2863 (1992) Fannes et al., CMP 144, 443 (1992) Östlund, Rommer, PRL 75, 3537 (1995) !!! Area-law in 2D !!! ✓ Reproduces area-law in 1D S ( L ) ∼ L S ( L ) = const D ∼ exp ( L )

  10. MPS & PEPS 1D 2D PEPS (TPS) MPS projected entangled-pair state Matrix-product state (tensor product state) Bond dimension D Bond dimension D 1 2 3 4 5 6 7 8 Physical indices (lattices sites) S. R. White, PRL 69, 2863 (1992) Fannes et al., CMP 144, 443 (1992) Östlund, Rommer, PRL 75, 3537 (1995) F. Verstraete, J. I. Cirac, cond-mat/0407066 Nishino, Hieida, et al., Prog. Theor. Phys. 105 (2001). Nishio, Maeshima, Gendiar, Nishino, cond-mat/0401115 ✓ Reproduces area-law in 1D ✓ Reproduces area-law in 2D S ( L ) = const S ( L ) ∼ L

  11. PEPS: Area law D L one “thick” bond of dimension D L A B ... ... L S ( A ) ≤ L log D ∼ L each cut auxiliary bond can contribute (at most) log D to the entanglement entropy ✓ Reproduces area-law in 2D The number of cuts scales with the cut length S ( L ) ∼ L

  12. MPS & PEPS 1D 2D PEPS (TPS) MPS projected entangled-pair state Matrix-product state (tensor product state) Bond dimension D Bond dimension D 1 2 3 4 5 6 7 8 Physical indices (lattices sites) S. R. White, PRL 69, 2863 (1992) Fannes et al., CMP 144, 443 (1992) Östlund, Rommer, PRL 75, 3537 (1995) F. Verstraete, J. I. Cirac, cond-mat/0407066 Nishino, Hieida, et al., Prog. Theor. Phys. 105 (2001). Nishio, Maeshima, Gendiar, Nishino, cond-mat/0401115 ✓ Reproduces area-law in 1D ✓ Reproduces area-law in 2D S ( L ) = const S ( L ) ∼ L

  13. Infinite PEPS (iPEPS) 1D 2D iMPS iPEPS infinite matrix-product state infinite projected entangled-pair state A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A Jordan, Orus, Vidal, Verstraete, Cirac, PRL (2008) Nishio, Maeshima, Gendiar, Nishino, cond-mat/0401115 ★ Work directly in the thermodynamic limit: No finite size and boundary effects!

  14. Infinite PEPS (iPEPS) 1D 2D iMPS iPEPS infinite matrix-product state infinite projected entangled-pair state A B A B A B B A B A B A A B A B A B B A B A B A A B A B A B B A B A B A ★ Work directly in the thermodynamic limit: No finite size and boundary effects!

  15. iPEPS with arbitrary unit cells 1D 2D iMPS iPEPS infinite matrix-product state with arbitrary unit cell of tensors D A B C D A H E F G H E D B D A C A H E F G H E B D D A C A H E F G H E here: 4x2 unit cell PC, White, Vidal, Troyer, PRB 84 (2011) ★ Run simulations with different unit cell sizes and compare variational energies

  16. Overview: Tensor network algorithms (ground state) MPS PEPS 2D MERA TN ansatz (variational) 1D MERA Find the best Compute (ground) state observables | ˜ � ˜ Ψ | O | ˜ Ψ � Ψ ⇥ Contraction of the iterative optimization tensor network imaginary time of individual tensors exact / approximate evolution (energy minimization)

  17. PART II: Contraction

  18. Contracting a tensor network (repetition)

  19. Pairwise contractions...

  20. Pairwise contractions...

  21. Pairwise contractions...

  22. Pairwise contractions...

  23. Pairwise contractions...

  24. Pairwise contractions... done! the order of contraction matters for the computational cost!!!

  25. Contracting a tensor network ★ Reshape tensors into matrices and multiply them with optimized routines (BLAS) dimension D u u i cost D 5 w = A B T w j v v ( ij ) A ( uv ) B T ( uv ) = w w dimension D 2 ★ Computational cost: multiply the dimensions of all legs (connected legs only once)

  26. Contracting an MPS � Ψ | Ψ ⇥ BAD! = � Ψ | Ψ ⇥ Good! =

  27. Contracting the PEPS � Ψ | Ψ ⇥ reduced tensors D 2

  28. Contracting the PEPS dimension D 2 Problem: how do we contract this?? no matter how we contract, we will get intermediate tensors with O(L) legs number of coefficients D 2L Exponentially increasing with L! NOT EFFICIENT

  29. Contracting the PEPS ★ Exact contraction of an PEPS is exponentially hard! use controlled approximate contraction scheme TRG Corner transfer MPS-based Tensor Renormalization Group matrix method approach (+HOTRG, SRG, HOSRG) Murg,Verstraete,Cirac, PRA75 ’07 Levin, Nave, PRL99 (2007) Nishino, Okunishi, JPSJ65 (1996) Jordan,et al. PRL79 (2008) Xie et al. PRL 103, (2009) Orus, Vidal, PRB 80 (2009) ★ Accuracy of the approximate contraction is controlled by “boundary dimension” χ TNR Tensor Network Renormalization ★ Convergence in needs to be carefully checked χ Evenbly & Vidal, PRL 115 (2015) ★ Overall cost: with χ ∼ D 2 Loop-TNR: O ( D 10 ... 14 ) Yang, Gu & Wen, arXiv:1512.04938

  30. Contracting the PEPS Example: 2D Heisenberg model (CTM) ★ Fast convergence − 0.6685 ★ Effect of finite D is much larger! E s − 0.669 ★ Be careful with D=4 “variational” energy!!! D=5 D=6 − 0.6695 0 0.02 0.04 0.06 0.08 0.1 1/ χ

  31. Contracting the PEPS ★ Exact contraction of an PEPS is exponentially hard! use controlled approximate contraction scheme TRG Corner transfer MPS-based Tensor Renormalization Group matrix method approach (+HOTRG, SRG, HOSRG) Murg,Verstraete,Cirac, PRA75 ’07 Nishino, Okunishi, JPSJ65 (1996) Levin, Nave, PRL99 (2007) Jordan,et al. PRL79 (2008) Orus, Vidal, PRB 80 (2009) Xie et al. PRL 103, (2009) ★ Accuracy of the approximate contraction is controlled by “boundary dimension” χ TNR Tensor Network Renormalization ★ Convergence in needs to be carefully checked χ Evenbly & Vidal, PRL 115 (2015 ★ Overall cost: with χ ∼ D 2 Loop-TNR: O ( D 10 ... 14 ) Yang, Gu & Wen, arXiv:1512.04938

  32. Contracting the PEPS using an MPS Verstraete, Murg, Cirac, Adv. in Phys. 57, 143 (2008) dimension D 2 this is an MPS this is an MPO (matrix product operator)

  33. Contracting the PEPS using an MPS Verstraete, Murg, Cirac, Adv. in Phys. 57, 143 (2008) dimension D 2 xD 2 this is an MPS with bond dimension D 2 x D 2 truncate the bonds to χ there are different techniques for the efficient MPO-MPS multiplication (SVD, variational optimization, zip-up algorithm...) Schollwöck, Annals of Physics 326, 96 (2011) Stoudenmire, White, New J. of Phys. 12, 055026 (2010).

  34. Contracting the PEPS using an MPS Verstraete, Murg, Cirac, Adv. in Phys. 57, 143 (2008) dimension χ proceed... ★ We can do this from several directions ★ Similar procedure when computing an expectation value

Download Presentation
Download Policy: The content available on the website is offered to you 'AS IS' for your personal information and use only. It cannot be commercialized, licensed, or distributed on other websites without prior consent from the author. To download a presentation, simply click this link. If you encounter any difficulties during the download process, it's possible that the publisher has removed the file from their server.

Recommend

Introduction to iPEPS (second lecture) Philippe Corboz, Institute for Theoretical Physics,

Introduction to iPEPS (second lecture) Philippe Corboz, Institute for Theoretical Physics,

Introduction to iPEPS (second lecture) Philippe Corboz, Institute for Theoretical Physics, University of Amsterdam Outline Part I: iPEPS ansatz Part II: Contraction of PEPS / iPEPS Interlude: Example application: the

995 views • 74 slides
? iPEPS 0 0 . 189(2) 0 . 217(4) 0 . 25 2D tensor network ansatz I. Niesen and P. Corboz,

? iPEPS 0 0 . 189(2) 0 . 217(4) 0 . 25 2D tensor network ansatz I. Niesen and P. Corboz,

2D tensor network study of the S=1 bilinear-biquadratic Heisenberg model Philippe Corboz, Institute for Theoretical Physics, University of Amsterdam AF phase Haldane phase 3-SL 120 phase ? iPEPS 0 0 . 189(2) 0 . 217(4) 0 . 25 2D

758 views • 51 slides
INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION

INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION

TRAINING COURSE ON SOCIAL PROTECTION & FORMALIZATION TRINIDAD AND TOBAGO MARCH 15, 2017 INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION Design of the NIS Assistance from the

547 views • 23 slides
INTRODUCTION Introduction 2/42 INTRODUCTION Alternations I am giving bread to the

INTRODUCTION Introduction 2/42 INTRODUCTION Alternations I am giving bread to the

Patterns of variation in the expression of case and agreement Andrs Brny Leiden University Centre for Linguistics 9 November 2019, BaSIS Workshop a.barany@hum.leidenuniv.nl INTRODUCTION Introduction 2/42 INTRODUCTION

1.37k views • 42 slides
Introduction Introduction Introduction Introduction Outline Motivation Failures

Introduction Introduction Introduction Introduction Outline Motivation Failures

Introduction Introduction Introduction Introduction Outline Motivation Failures Definition and concepts Process and product properties Principles SE Approaches Motivation Software and the economy The economies of

1.03k views • 49 slides
Introduction Introduction Introduction Nationwide Cause for Concern 1

Introduction Introduction Introduction Nationwide Cause for Concern 1

Introduction Introduction Introduction Nationwide Cause for Concern 1 https://www.cdc.gov/std/stats17/infographic.htm Introduction Epidemiology in Western Colorado CPHE Surveillance Data Report, May 2018 Chlamydia Epidemiology in Western

977 views • 39 slides
DAQ introduction Purpose of this talk : (1) Introduction for those who have not been in every

DAQ introduction Purpose of this talk : (1) Introduction for those who have not been in every

DAQ introduction Purpose of this talk : (1) Introduction for those who have not been in every meeting (2) Some ideas of diagrams for section 7.1 (Introduction) of TP. Giles Barr 25 Jan 2018, Pembroke College On DUNE, the reality is there is a

474 views • 10 slides
INTRODUCTION cf. Schneider, Chapter 1 INTRODUCTION THEY ARE OUT TO GET YOU INTRODUCTION WHAT

INTRODUCTION cf. Schneider, Chapter 1 INTRODUCTION THEY ARE OUT TO GET YOU INTRODUCTION WHAT

ADVANCED COMPUTER SECURITY INTRODUCTION cf. Schneider, Chapter 1 INTRODUCTION THEY ARE OUT TO GET YOU INTRODUCTION WHAT IS SECURITY? A systems security policies describe What the system is supposed to do Store and provide access

479 views • 16 slides
Overview Introduction to SMIL Introduction to W3C and XML Introduction to SMIL

Overview Introduction to SMIL Introduction to W3C and XML Introduction to SMIL

Overview Introduction to SMIL Introduction to W3C and XML Introduction to SMIL Writing a SMIL file Adding Media Objects Martin Halvey Clipping media files January 2008 Timing and Synchronising Layout World Wide

387 views • 14 slides
14 Introduction Introduction Bad guys can put malware into Example: hosts

14 Introduction Introduction Bad guys can put malware into Example: hosts

Introduction Introduction source Encapsulation ISO/OSI reference model message M application segment H t H t M transport network datagram H n H n H t M presentation: allow applications to frame link H l H n H t M interpret

583 views • 4 slides
Previous Lecture Introduction to SMIL 2 Introduction to W3C and XML Introduction to SMIL

Previous Lecture Introduction to SMIL 2 Introduction to W3C and XML Introduction to SMIL

Previous Lecture Introduction to SMIL 2 Introduction to W3C and XML Introduction to SMIL Writing a SMIL file Adding Media Objects Martin Halvey Clipping media files January 2008 Timing and Synchronising Layout

322 views • 8 slides
Lecture 1 Andreas Habegger Introduction Zynq Introduction Zynq Introduction Zynq PS vs. PL

Lecture 1 Andreas Habegger Introduction Zynq Introduction Zynq Introduction Zynq PS vs. PL

Introduction Zynq Lecture 1 Andreas Habegger Introduction Zynq Introduction Zynq Introduction Zynq PS vs. PL Processing System Processor Peripherals Data Buses AXI Bus Conclusion BTE5380 - Embedded Systems Octobre 2014 Andreas Habegger

1.81k views • 30 slides
2018.06 01 SMILE5 Introduction S E 5 02 Alpha Cloud M I L 03 Company Introduction 04

2018.06 01 SMILE5 Introduction S E 5 02 Alpha Cloud M I L 03 Company Introduction 04

2018.06 01 SMILE5 Introduction S E 5 02 Alpha Cloud M I L 03 Company Introduction 04 Project References S E 5 PART.1 SMILE5 Introduction M I L Introduction STORION-SMILE5 Residential Energy Storage System Inverter 5kW Output Input

1.06k views • 25 slides
Introduction to CICS Course introduction Course introduction What is CICS? What is an

Introduction to CICS Course introduction Course introduction What is CICS? What is an

Introduction to CICS Course introduction Course introduction What is CICS? What is an application server? Why use an application server? Course introduction Services provided by CICS How CICS applications are defined Scaling CICS to meet

1.84k views • 146 slides
Network Visualization Introduction Presented by Shahed Introduction Introduction Basic

Network Visualization Introduction Presented by Shahed Introduction Introduction Basic

Network Visualization Introduction Presented by Shahed Introduction Introduction Basic building blocks Node Links (relationship between nodes) Spatial information Network data

584 views • 9 slides
Introduction Introduction (1) Main argument of the paper: universities are not only

Introduction Introduction (1) Main argument of the paper: universities are not only

Le jugement des pairs comme outil de management des universits Peer-review as a management tool Christine Musselin (CSO, Sciences Po et CNRS) Avril 2013, Lirrsistible ascension du capitalisme acadmique Introduction Introduction (1)

590 views • 28 slides
Human Action Recognition Using Semi-Latent Topic Models Yang Wang and Greg Mori , 2009 SE367

Human Action Recognition Using Semi-Latent Topic Models Yang Wang and Greg Mori , 2009 SE367

Human Action Recognition Using Semi-Latent Topic Models Yang Wang and Greg Mori , 2009 SE367 Paper Presentation - Deepak Pathak 10222 Introduction Human Action Recognition ( What ?) Still Images (eg: Poselets) v/s Video

378 views • 12 slides
Nuclear Level Density, Underlying Physics, and Constant Temperature Model Vladimir Zelevinsky

Nuclear Level Density, Underlying Physics, and Constant Temperature Model Vladimir Zelevinsky

Nuclear Level Density, Underlying Physics, and Constant Temperature Model Vladimir Zelevinsky NSCL/FRIB Michigan State University June 22, 2018 In collaboration with Mihai Horoi Sofia Karampagia Roman Senkov Antonio Renzaglia Alex

1.85k views • 83 slides
Higher-order tensor renormalization group with the corner transfer matrix Satoshi Morita (ISSP ,

Higher-order tensor renormalization group with the corner transfer matrix Satoshi Morita (ISSP ,

Higher-order tensor renormalization group with the corner transfer matrix Satoshi Morita (ISSP , UTokyo) TNSAA 2019-2020 @ National Cheng-Chi University, Taipei, Taiwan 1. Higher-order tensor renormalization group with the corner transfer

614 views • 32 slides
Computer Graphics (CS 543) Lecture 3 (Part 3): Implementing Transformations Prof Emmanuel Agu

Computer Graphics (CS 543) Lecture 3 (Part 3): Implementing Transformations Prof Emmanuel Agu

Computer Graphics (CS 543) Lecture 3 (Part 3): Implementing Transformations Prof Emmanuel Agu Computer Science Dept. Worcester Polytechnic Institute (WPI) Objectives Learn how to implement transformations in OpenGL Rotation

502 views • 38 slides
The Classical Turing Machine 2IT70 Finite Automata and Process Theory Technische Universiteit

The Classical Turing Machine 2IT70 Finite Automata and Process Theory Technische Universiteit

The Classical Turing Machine 2IT70 Finite Automata and Process Theory Technische Universiteit Eindhoven May 14, 2014 Classical Turing machine architecture Input yes/no Automaton Tape Tape aabbaa a b no input and no final states initial

237 views • 13 slides
Assessing the performance of RegCM4- chem in the simulation of ozone levels: Summer circulation

Assessing the performance of RegCM4- chem in the simulation of ozone levels: Summer circulation

Aristotle University Department of Meteorology of Thessaloniki and Climatology Assessing the performance of RegCM4- chem in the simulation of ozone levels: Summer circulation over Mediterranean and stratospheric intrusion events

612 views • 26 slides
Concurrency Theory Winter Semester 2019/20 Lecture 11: Trace Equivalence Joost-Pieter Katoen and

Concurrency Theory Winter Semester 2019/20 Lecture 11: Trace Equivalence Joost-Pieter Katoen and

Concurrency Theory Winter Semester 2019/20 Lecture 11: Trace Equivalence Joost-Pieter Katoen and Thomas Noll Software Modeling and Verification Group RWTH Aachen University https://moves.rwth-aachen.de/teaching/ws-19-20/ct/ Introduction

795 views • 66 slides
Open Government Initiative: Community Perspective Chip Rosenthal chip@unicom.com 512-573-5174

Open Government Initiative: Community Perspective Chip Rosenthal chip@unicom.com 512-573-5174

Open Government Initiative: Community Perspective Chip Rosenthal chip@unicom.com 512-573-5174 Open Austin www.open-austin.org After much excitement and anticipation, our Open Government efforts are flagging and may be on the brink of

216 views • 10 slides

More recommend