Tensor network renormalization

Transcription

1 Coogee'15 Sydney Quantum Information Theory Workshop Tensor network renormalization Guifre Vidal

2 In collaboration with GLEN EVENBLY IQIM Caltech UC Irvine

3 Quantum Mechanics Niels Bohr Albert Einstein Wolfgang Pauli Paul Dirac iħ t Ψ = HΨ Schrodinger equation Erwin Schrodinger Enrico Fermi Werner Heisenberg Richard Feynman

4 Exotic phases of quantum matter: Collective quantum many-body phenomena superfluids superconductors supersolids (?) spin liquids Bose-Einstein condensation fractional quantum Hall effect quantum criticality topological order

5 There is a problem Paul Dirac The fundamental laws necessary for the mathematical treatment of a large part of physics and the whole of chemistry are thus completely known, and the difficulty lies only in the fact that application of these laws leads to equations that are too complex to be solved. iħ t Ψ = HΨ Schrodinger equation? Emergence: long-distance physics is often radically different from short-distance physics

6 and there is a solution: the Renormalization Group Kenneth Wilson Renormalization group flow: λ=λ H c 3 λ<λ H c λ>λ H c 1 2 Leo Kadanoff m 0 m = 0 H order H crit H disorder

7 and there is a solution: the Renormalization Group Nice! But how do we do this in practice? Outline of this talk entanglement renormalization (old stuff) Real space RG transformation for ground state wave-functions Hamiltonians tensor network renormalization (new stuff) Real space RG of Euclidean path integral TNR MERA (+ thermal states!) Theory of minimal updates Conformal transformations arxiv: arxiv:1501.xxx arxiv:1501.yyy arxiv:1501.zzz

8 Outline entanglement renormalization (old stuff) Real space RG transformation for ground state wave-functions Hamiltonians tensor network renormalization (new stuff) Real space RG of Euclidean path integral TNR MERA (+ thermal states!) Theory of minimal updates Conformal transformations arxiv: arxiv:1501.xxx arxiv:1501.yyy arxiv:1501.zzz

9 Entanglement renormalization (old stuff) L L u disentangler w isometry L L

10 Entanglement renormalization (old stuff) removal of all short-range entanglement L L u disentangler w isometry some short-range entanglement is not removed L L

11 RG transformation on ground state wave-functions Ψ = Ψ

12 RG transformation on ground state wave-functions Ψ = L L L Multi-scale entanglement renormalization ansatz (MERA) topological order quantum criticality holography

13 RG flow in the space of wave-functions, as parameterized by tensors u, w u w (u, w) 3 λ=λ c (u, w)1 λ<λ c (u, w) 2 λ>λ c m 0 m = 0 (u, w) order (u, w) crit (u, w) disorder

14 RG transformation on Hamiltonians H = H

15 RG transformation on Hamiltonians Local Hamiltonian H = h i i h u u w w = =

16 RG transformation on Hamiltonians Local Hamiltonian H = h i i h u u w w = =

17 RG transformation on Hamiltonians Local Hamiltonian H = h i i h u u w = h = w =

18 RG flow in the space of local Hamiltonias, as parameterized by a strictly local Hamiltonian term h H = i h i h 1 λ<λ c h 3 λ=λ c h 2 λ>λ c m 0 m = 0 h order h crit h disorder

19 entanglement renormalization (old stuff) Real space RG transformation for ground state wave-functions Hamiltonians MANY OPEN QUESTIONS RG on Euclidean path integral (Hamiltonian Lagrangian) MERA for thermal states? Classical partition functions? When is MERA a good approximation? Theory of minimal updates (impurities, boundaries, etc) Better algorithms? ( branching, 2D)

20 Outline entanglement renormalization (old stuff) Real space RG transformation for ground state wave-functions Hamiltonians tensor network renormalization (new stuff) Real space RG of Euclidean path integral TNR MERA (+ thermal states!) Theory of minimal updates Conformal transformations arxiv: arxiv:1501.xxx arxiv:1501.yyy arxiv:1501.zzz

21 Tensor network renormalization (new stuff) local Hamiltonian H ground state Ψ Ψ = lim β e βh φ 0 e βh φ 0 Suzuki-Trotter decomposition e βh e εh Ae εh Be εh Ae εh B H A = H B = h i,i+1 i odd h i,i+1 i even e εh A e εh B e τh e τh e εh A e εh B e εh A

22 e εh A e εh B e εh A e τh e τh Ψ e εh B Ψ Ψ e εh A e εh A e εh B e τh e τh e εh A e εh B Ψ e εh A

23 e εh A e εh B e εh A e τh e τh Ψ e εh B Ψ o Ψ e εh A e εh A e τh o e εh B e τh e εh A e εh B Ψ e εh A

24 e τh A e τh u s v u s s v e τh = w d = w d w d w A =

25 Ψ o Ψ Ψ Ψ =

26 Tensor network Euclidean path integral / Euclidean time evolution operator A Classical partition function Overlap between two PEPS Monte Carlo sampling Plan: sample important amplitudes! e εh A e τh e εh B e τh sampling error e εh A (sign problem!) e εh B e εh A

27 Tensor network Euclidean path integral / Euclidean time evolution operator A Classical partition function Overlap between two PEPS Tensor network techniques (I) Plan: sum over all amplitudes boundary MPS (or CTM) poor approximation at/near criticality

28 Tensor network Euclidean path integral / Euclidean time evolution operator A Classical partition function Overlap between two PEPS Tensor network techniques (II) A A Plan: sum over all amplitudes by coarse-grainning the tensor network A A

29 Tensor network Euclidean path integral / Euclidean time evolution operator A Classical partition function Overlap between two PEPS Can we define an RG flow in the space of tensors A? A A A A 1 λ<λ c λ=λ A c 3 λ>λ A c 2 m 0 m = 0 A order A crit A disorder

30 Tensor Renormalization Group (TRG) (Levin, Nave, 2006) A A truncated SVD initial tensor network coarse-grained tensor network

31 Tensor Renormalization Group (TRG) (Levin, Nave, 2006) Accurate results for gapped 1D systems Many improvements and generalizations: Second Renormalization Group (SRG) Tensor Entanglement Filtering Renormalization (TEFR) Higher Order Tensor Renormalization Group (HOTRG) + many more (Xie, Jiang, Weng, Xiang, 2008) Improved accuracy, extensions to 2D systems (!?), etc (Gu, Wen, 2009) (Xie, Chen, Qin, Zhu, Yang, Xiang, 2012) However, TRG fails to remove some of the short-range correlations. Not a proper RG flow (wrong structure of fixed points) Computationally not sustainable at quantum criticality

32 Tensor Network Renormalization (Evenbly, Vidal, 2014) Introduce disentanglers to remove all short-range correlations. Proper RG flow (correct structure of fixed points) Works near/at quantum criticality Can be generalized to D > 1 dimensions

33

34 Proper RG flow: 2D classical Ising e.g. disordered phase T > T c TRG (Levin, Nave 2006) A 1 A 2 A 3 A 4 T = 1.1 T c T = 2.0 T c should converge to the same (trivial) fixed point, but they don t! TNR (new) T = 1.1 T c T = 2.0 T c converge to the same fixed point (containing only universal information of the phase)

35 Proper RG flow: 2D classical Ising Tensor Network Renormalization (TNR): Converges to one of three RG fixed points, consistent with a proper RG flow A 1 A 2 A 3 A 4 below critical T = 0.9 T c ordered (Z2) fixed point critical T = T c critical (scaleinvariant) fixed point above critical T = 1.1 T c disordered (trivial) fixed point

36 Sustainable RG flow: 2D classical Ising Does TRG give a sustainable RG flow? RG flow at criticality 10 0 s = 0 s = 2 s = 6 s = 1 Spectrum of A (s) χ Bond dimension required to maintain fixed truncation error (~10 3 ): Computational cost: TRG O(χ 6 ): TRG: ~10 ~20 ~40 > > Cost of TRG scales exponentially with RG iteration!

37 Sustainable RG flow: 2D classical Ising Spectrum of A (s) Does TRG give a sustainable RG flow? RG flow at criticality 10 0 s = 0 s = 2 s = 6 s = TRG TNR χ Bond dimension required to maintain fixed truncation error (~10 3 ): TRG: ~10 ~20 ~40 > 100 TNR: ~10 ~10 ~10 ~10 Sustainable Computational costs: O(χ 6 ) O(kχ 6 ) TRG : TNR : >

38 Proper RG flow: 2D classical Ising critical point: T = T c TNR bond dimension: χ = 4 A 1 A 2 A 3 A 4 A 8 A 7 A 6 A 5 A 9 A 10 A 11 A 12

39 Proper RG flow: 2D classical Ising more difficult! T = T c TNR bond dimension: χ = 4 A 1 A 2 A 3 A 4 A 8 A 7 A 6 A 5 A 9 A 10 A 11 A 12

40 Outline entanglement renormalization (old stuff) Real space RG transformation for ground state wave-functions Hamiltonians tensor network renormalization (new stuff) Real space RG of Euclidean path integral TNR MERA (+ thermal states!) Theory of minimal updates Conformal transformations arxiv: arxiv:1501.xxx arxiv:1501.yyy arxiv:1501.zzz

41 TNR yields MERA (Evenbly, Vidal, 2015, in prep)

42 TNR yields MERA (Evenbly, Vidal, 2015, in prep)

43 Traditional perspective MERA as a variational class of many-body states New perspective MERA as a byproduct of tensor network renormalization (TNR) reshape an existing representation of the ground state TNR truncation errors: accuracy certificate energy minimization trace optimization 1000 s of iterations on scale 1 iteration on scale non-local in scale local in scale

44 Extra bonus: thermal MERA (Evenbly, Vidal, 2015, in prep) e βh

45 Outline entanglement renormalization (old stuff) Real space RG transformation for ground state wave-functions Hamiltonians tensor network renormalization (new stuff) Real space RG of Euclidean path integral TNR MERA (+ thermal states!) Theory of minimal updates Conformal transformations arxiv: arxiv:1501.xxx arxiv:1501.yyy arxiv:1501.zzz

46 Extra bonus: theory of minimal updates (Evenbly, Vidal, 2015, in prep) homogeneous system: s w u x

47 Extra bonus: theory of minimal updates (Evenbly, Vidal, 2015, in prep) same system with impurity: C(R) s w L u w R w u x R

48 Example: critical Ising model with one modified bond H = ( X r X r+1 + Z r ) i H imp = 1 α X 0 X 1

49 Extra bonus: theory of minimal updates homogeneous system:

50 Extra bonus: theory of minimal updates same system with impurity:

51 Outline entanglement renormalization (old stuff) Real space RG transformation for ground state wave-functions Hamiltonians tensor network renormalization (new stuff) Real space RG of Euclidean path integral TNR MERA (+ thermal states!) Theory of minimal updates Conformal transformations arxiv: arxiv:1501.xxx arxiv:1501.yyy arxiv:1501.zzz

52 Extra bonus: conformal transformations (Evenbly, Vidal, 2015, in prep) Conformal transformation 1: Upper half-plane to AdS 2 dx 2 + dτ 2 = τ 2 dx 2 τ 2 + dτ2 τ 2 dx2 τ 2 + dτ2 τ 2 = dx2 + ds2 22s s = log 2 (τ) Holographic description

53 Extra bonus: conformal transformations (Evenbly, Vidal, 2015, in prep) plane Conformal transformation 2: Plane to cylinder τ x (radial quantization in CFT) z x + iτ z = 2 w w s + iθ s log 2 (x 2 + τ 2 ) s θ cylinder Extraction of scaling dimensions

54 Summary entanglement renormalization (old stuff) Real space RG transformation for ground state wave-functions Hamiltonians tensor network renormalization (new stuff) Real space RG of Euclidean path integral TNR MERA (+ thermal states!) Theory of minimal updates Conformal transformations arxiv: arxiv:1501.xxx arxiv:1501.yyy arxiv:1501.zzz

55 THANKS!