Loop optimization for tensor network renormalization

Similar documents
Tensor network renormalization

Introduction to tensor network state -- concept and algorithm. Z. Y. Xie ( 谢志远 ) ITP, Beijing

Renormalization of Tensor- Network States Tao Xiang

Scale invariance on the lattice

Renormalization of Tensor Network States

Tensor network renormalization

Tensor network methods in condensed matter physics. ISSP, University of Tokyo, Tsuyoshi Okubo

c 2017 Society for Industrial and Applied Mathematics

Renormalization of Tensor Network States. Partial order and finite temperature phase transition in the Potts model on irregular lattice

arxiv: v3 [cond-mat.str-el] 15 Sep 2015

Efficient time evolution of one-dimensional quantum systems

Tensor Networks, Renormalization and Holography (overview)

arxiv: v1 [quant-ph] 18 Jul 2017

arxiv: v1 [cond-mat.str-el] 24 Sep 2015

From Path Integral to Tensor Networks for AdS/CFT

An introduction to tensornetwork

Holographic Branching and Entanglement Renormalization

Simulating Quantum Systems through Matrix Product States. Laura Foini SISSA Journal Club

Fermionic tensor networks

The density matrix renormalization group and tensor network methods

Machine Learning with Quantum-Inspired Tensor Networks

Disentangling Topological Insulators by Tensor Networks

Matrix product states for the fractional quantum Hall effect

It from Qubit Summer School

Quantum Hamiltonian Complexity. Itai Arad

The uses of Instantons for classifying Topological Phases

Lecture 3: Tensor Product Ansatz

Frustration-free Ground States of Quantum Spin Systems 1

arxiv: v1 [quant-ph] 6 Jun 2011

Time Evolving Block Decimation Algorithm

Quantum Fields, Gravity, and Complexity. Brian Swingle UMD WIP w/ Isaac Kim, IBM

arxiv: v2 [quant-ph] 31 Oct 2013

Manifestly diffeomorphism invariant classical Exact Renormalization Group

4 Matrix product states

Fermionic topological quantum states as tensor networks

Many-body entanglement witness

IPAM/UCLA, Sat 24 th Jan Numerical Approaches to Quantum Many-Body Systems. QS2009 tutorials. lecture: Tensor Networks.

Rigorous free fermion entanglement renormalization from wavelets

MERA for Spin Chains with Critical Lines

Machine Learning with Tensor Networks

3 Symmetry Protected Topological Phase

News on tensor network algorithms

Matrix-Product states: Properties and Extensions

Fermionic tensor networks

Quantum simulation with string-bond states: Joining PEPS and Monte Carlo

Entanglement signatures of QED3 in the kagome spin liquid. William Witczak-Krempa

Classification of Symmetry Protected Topological Phases in Interacting Systems

Introduction to Tensor Networks: PEPS, Fermions, and More

Regularization Physics 230A, Spring 2007, Hitoshi Murayama

Quantum quenches in 2D with chain array matrix product states

Holographic Geometries from Tensor Network States

Zhong-Zhi Xianyu (CMSA Harvard) Tsinghua June 30, 2016

Positive Tensor Network approach for simulating open quantum many-body systems

Towards a manifestly diffeomorphism invariant Exact Renormalization Group

NUMERICAL METHODS FOR QUANTUM IMPURITY MODELS

Loop corrections in Yukawa theory based on S-51

Shunsuke Furukawa Condensed Matter Theory Lab., RIKEN. Gregoire Misguich Vincent Pasquier Service de Physique Theorique, CEA Saclay, France

The Density Matrix Renormalization Group: Introduction and Overview

VI.D Self Duality in the Two Dimensional Ising Model

arxiv: v4 [cond-mat.stat-mech] 13 Mar 2009

Renormalization Group for the Two-Dimensional Ising Model

Matrix Product Operators: Algebras and Applications

Gravity Actions from Tensor Networks

One Loop Tests of Higher Spin AdS/CFT

VI.D Self Duality in the Two Dimensional Ising Model

Momentum-space and Hybrid Real- Momentum Space DMRG applied to the Hubbard Model

Entanglement in quantum phase transition

Richard Williams C. S. Fischer, W. Heupel, H. Sanchis-Alepuz

Quantum many-body systems and tensor networks: simulation methods and applications

Manuscript code: BH11160 RECVD: Mon Feb 8 16:14: Resubmission to: Physical Review B Resubmission type: resubmit

Entanglement in Topological Phases

Plaquette Renormalized Tensor Network States: Application to Frustrated Systems

Matrix product approximations to multipoint functions in two-dimensional conformal field theory

Spectral action, scale anomaly. and the Higgs-Dilaton potential

RG Limit Cycles (Part I)

Z2 topological phase in quantum antiferromagnets. Masaki Oshikawa. ISSP, University of Tokyo

Universal Dynamics from the Conformal Bootstrap

Exact Functional Renormalization Group for Wetting Transitions in Dimensions.

arxiv: v2 [cond-mat.str-el] 6 Nov 2013

Boundary Degeneracy of Topological Order

(r) 2.0 E N 1.0

Local Algorithms for 1D Quantum Systems

arxiv: v1 [cond-mat.str-el] 7 Aug 2011

Gapped spin liquid with Z[subscript 2] topological order for the kagome Heisenberg model

Isotropic Entanglement

4 4 and perturbation theory

Critical Region of the QCD Phase Transition

Implications of Poincaré symmetry for thermal field theories

Dephasing, relaxation and thermalization in one-dimensional quantum systems

Excursion: MPS & DMRG

Tensor network simulations of strongly correlated quantum systems

Approaching conformality in lattice models

REVIEW REVIEW. Quantum Field Theory II

Quantum Field Theory II

Tensor network simulation of QED on infinite lattices: learning from (1 + 1)d, and prospects for (2 + 1)d

Phase Transitions and Renormalization:

Inverse square potential, scale anomaly, and complex extension

Valence Bonds in Random Quantum Magnets

The Ising model Summary of L12

Dynamical mean field approach to correlated lattice systems in and out of equilibrium

Transcription:

Yukawa Institute for Theoretical Physics, Kyoto University, Japan June, 6 Loop optimization for tensor network renormalization Shuo Yang! Perimeter Institute for Theoretical Physics, Waterloo, Canada Zheng-Cheng Gu (CUHK, PI) Xiao-Gang Wen (MIT, PI) arxiv:5.98

Overview tensor network + renormalization group = tensor network renormalization! Partition function of classical statistical system!!! Euclidean path integral of D quantum system!!! Physical properties of D quantum system e h Z = X { } exp( X hiji i j)

Overview tensor network + renormalization group = tensor network renormalization boundary MPS fully isotropic coarse-graining

Overview tensor network + renormalization group = tensor network renormalization How the physics change with scale scale transformation Non-perturbative approach, suitable for strongly correlated systems Reproduce long-range physics

Overview tensor network + renormalization group = tensor network renormalization How the physics change with scale scale transformation Non-perturbative approach, suitable for strongly correlated systems Reproduce long-range physics RG flow Aim! Remove short-range entanglement / correlations Generate proper RG flow & correct fixed points Recover scale invariance at criticality A C B

Real space renormalization group D quantum: real space NRG entanglement DMRG S. White (99) Real space analogue of NRG Wilson (975) breakdown for particle in a box problem D classical: block spin entanglement TRG Levin & Nave (7) L.P. Kadanoff (966)

Tensor renormalization group Levin & Nave (7) LN-TNR Three steps. Deform tensors, make a truncation. Coarse graining. Renormalize tensors (multiply tensor by a constant factor)

Tensor renormalization group SVD singular value decomposition (SVD) keep only keep the largest keep SVD original keep = S S T U singular values cost function keep keep p = p V S keep = S keep = /

Drawbacks of LN-TNR Off criticality cannot remove conner-double-line (CDL) tensors cannot give the correct structures of fixed points D D D D = D B @ D D... D C A A C B

Drawbacks of LN-TNR At criticality cannot explicitly recover scale invariance cannot completely remove short-range entanglement

Drawbacks of LN-TNR At criticality cannot explicitly recover scale invariance cannot completely remove short-range entanglement Example classical Ising model partition function Z = X { } exp( X hiji i j)

Drawbacks of LN-TNR At criticality cannot explicitly recover scale invariance cannot completely remove short-range entanglement Example classical Ising model partition function Z = X { } exp( X i j) hiji local tensor T = T Ising u,l,d,r T Ising,,, = e, T Ising,,, = e, T Ising,,, = e, T Ising,,, = e, others =.

Drawbacks of LN-TNR Ising CFT central charge c =/ scaling dimensions I

Drawbacks of LN-TNR Ising CFT central charge c =/ scaling dimensions I Calculate scaling dimensions transfer matrix eigenvalues of the transfer matrix c, Zheng-Cheng Gu and Xiao-Gang Wen, Phys. Rev. B 8, 55 (9).

Drawbacks of LN-TNR Ising CFT central charge c =/ scaling dimensions I Calculate scaling dimensions transfer matrix eigenvalues of the transfer matrix c, Zheng-Cheng Gu and Xiao-Gang Wen, Phys. Rev. B 8, 55 (9). LN-TNR c, even, odd Loop-I, D= scaling dimensions change with RG step cannot explicitly recover scale invariance 5 5 RG step

Drawbacks of LN-TNR Ising CFT central charge c =/ scaling dimensions I Calculate scaling dimensions transfer matrix eigenvalues of the transfer matrix c, Zheng-Cheng Gu and Xiao-Gang Wen, Phys. Rev. B 8, 55 (9). LN-TNR c, even, odd Loop-I, D= scaling dimensions change with RG step cannot explicitly recover scale invariance high-index parts will destroy low-index parts accuracy is fine, stability is bad 5 5 RG step

Improvements of TRG off-critical critical TRG (Levin & Nave, 6) SRG (Xie, Jiang, Weng, Xiang, 8) TEFR (Gu & Wen, 9) TNR (Evenbly & Vidal, ) Loop-TRG! (current)

Improvements of TRG off-critical critical TRG (Levin & Nave, 6) SRG (Xie, Jiang, Weng, Xiang, 8) TEFR (Gu & Wen, 9) TNR (Evenbly & Vidal, ) Loop-TRG (current)

Renormalization Group Momentum space Tree level approximation no loop Real space Tree level approximation no loop Beyond tree level one loop Beyond tree level one loop

Renormalization Group Momentum space Tree level approximation no loop Real space Tree level approximation no loop momentum b ignored fast slow simply ignore the fast modes Beyond tree level one loop Beyond tree level one loop

Renormalization Group Momentum space Tree level approximation no loop Real space Tree level approximation no loop momentum b ignored fast slow simply ignore the fast modes Beyond tree level one loop Beyond tree level one loop momentum b ignored fast slow integrate over the fast modes

Renormalization Group Momentum space Tree level approximation no loop Real space Tree level approximation no loop momentum b ignored fast slow simply ignore the fast modes LN-TNR remove short-range entanglement by a local SVD Beyond tree level one loop Beyond tree level one loop momentum b ignored fast slow integrate over the fast modes

Renormalization Group Momentum space Tree level approximation no loop Real space Tree level approximation no loop momentum b ignored fast slow simply ignore the fast modes LN-TNR remove short-range entanglement by a local SVD Beyond tree level one loop Beyond tree level one loop momentum b ignored fast slow integrate over the fast modes LN-TNR +loop optimization further remove short-range entanglement inside a loop! This work!

Algorithms of Loop-TNR

Algorithms of Loop-TNR Part One: Entanglement filtering

Algorithms of Loop-TNR Part One: Entanglement filtering Part Two: Optimizing tensors on a loop

Algorithms of Loop-TNR Part One: Entanglement filtering Together: Complete remove short-range entanglement Part Two: Optimizing tensors on a loop

Part One Entanglement filtering How. Find & insert projectors. Define new tensors 8 7 6 5 5 6 7 8 P L P L P L P L T T T T P R P R P R P R

Part One Entanglement filtering How. Find & insert projectors. Define new tensors 8 7 6 5 5 6 7 8 P L P L P L P L T T T T P R P R P R P R Aim Remove conner double line (CDL) tensors Generate local canonical gauge Zheng-Cheng Gu and Xiao-Gang Wen, Phys. Rev. B 8, 55 (9).

Part Two Optimizing tensors on a loop

Part Two Optimizing tensors on a loop cost function LN-TNR Loop-TNR

Part Two Optimizing tensors on a loop solve the linear equation

Part Two Optimizing tensors on a loop

Part Two Optimizing tensors on a loop Loop-TNR Results L =

Part Two Optimizing tensors on a loop Loop-TNR Results scaling dimensions does not change with scale ~ scale invariance a clear gap between high-level parts and low-level parts L =

Stability c, even, odd (a) LN-TNR =6, L= (b) Loop-TNR, =6, L= 5 5 c, even, odd (c) LN-TNR =, L= (d) Loop-TNR, =, L= 5 5 6 (e) LN-TNR, =6, L= 6 (f) Loop-TNR, =6, L= c, even, odd 5 5 5 5 Iteration step Iteration step

Stability c, even, odd (a) LN-TNR =6, L= (b) Loop-TNR, =6, L= = 6 remain accurate up to iteration steps 5 5 c, even, odd (c) LN-TNR =, L= 5 5 (d) Loop-TNR, =, L= even longer for = the proper RG flow last longer for larger 6 (e) LN-TNR, =6, L= 6 (f) Loop-TNR, =6, L= c, even, odd 5 5 5 5 Iteration step Iteration step

Stability c, even, odd (a) LN-TNR =6, L= (b) Loop-TNR, =6, L= 5 5 c, even, odd (c) LN-TNR =, L= (d) Loop-TNR, =, L= increasing, more scaling dimensions can be resolved 5 5 6 (e) LN-TNR, =6, L= 6 (f) Loop-TNR, =6, L= c, even, odd 5 5 5 5 Iteration step Iteration step

Stability c, even, odd (a) LN-TNR =6, L= (b) Loop-TNR, =6, L= effectively, = 6 5 5 c, even, odd (c) LN-TNR =, L= (d) Loop-TNR, =, L= c, even, odd 6 5 5 5 (e) LN-TNR, =6, L= 6 5 (f) Loop-TNR, =6, L= effectively, = 6 = 56 much more scaling dimensions can be read off accuracy is higher 5 5 Iteration step Iteration step

Stability c, even, odd (a) LN-TNR =6, L= (b) Loop-TNR, =6, L= effectively, = 6 5 5 c, even, odd (c) LN-TNR =, L= (d) Loop-TNR, =, L= c, even, odd 6 5 5 5 (e) LN-TNR, =6, L= 5 5 Iteration step 6 5 (f) Loop-TNR, =6, L= Iteration step effectively, = 6 = 56 much more scaling dimensions can be read off accuracy is higher infinite ~ infinite dimensional fixed point tensor described by Ising CFT

Accuracy

Relative error of the free energy per site f (Tc) - -5-6 -7-8 -9 - Loop-TNR (a) Critical LN-TNR 5 6 f (T) - -5-7 (b) Off-critical, =8 LN-TNR Loop-TNR -9.9.95..5. T/T c

Relative error of the free energy per site At critical point, the error of Loop-TNR decays much faster than the error of LN-TNR f (Tc) - -5-6 -7-8 -9 - Loop-TNR (a) Critical LN-TNR 5 6 f (T) - -5-7 (b) Off-critical, =8 LN-TNR Loop-TNR -9.9.95..5. T/T c almost decays exponentially with bond dimension

Relative error of the free energy per site At critical point, the error of Loop-TNR decays much faster than the error of LN-TNR f (Tc) - -5-6 -7-8 -9 - Loop-TNR (a) Critical LN-TNR 5 6 f (T) - -5-7 (b) Off-critical, =8 LN-TNR Loop-TNR -9.9.95..5. T/T c almost decays exponentially with bond dimension At off-critical points, the errors almost remains a constant for all temperatures near the critical point

Loop TRG with reflection symmetry i i l j j k k l k j j k i i l l (a) (b) axis of symmetry global C symmetry (f) (e) (d) (c) fix the gauge (g) (h) j i k l T (i) S ideal case: explicitly invariant fixed point tensor S

Fixed point tensors LN-TNR absolute value LN-TNR absolute difference low-index parts are approximately invariant Loop-TNR absolute value Loop-TNR absolute difference all tensor elements are nearly invariant ~ scale invariance

Summary

Summary Real space TNR! Tree level: LN-TNR (Levin & Nave, 7) One loop: GW-TNR, EV-TNR, Loop-TNR

Summary Real space TNR! Tree level: LN-TNR (Levin & Nave, 7) One loop: GW-TNR, EV-TNR, Loop-TNR D algorithm LN-TNR + D algorithm itebd GW-TNR (Gu & Wen, 9) Loop-TNR Part One DMRG Varational MPS Loop-TNR Part Two MERA EV-TNR (Evenbly & Vidal, )

Thank you!