Boundary theories for spins in lattices J. IGNACIO CIRAC GGI, - - PowerPoint PPT Presentation

boundary theories for spins in lattices
SMART_READER_LITE
LIVE PREVIEW

Boundary theories for spins in lattices J. IGNACIO CIRAC GGI, - - PowerPoint PPT Presentation

Boundary theories for spins in lattices J. IGNACIO CIRAC GGI, Florence, May 24, 2012 IC, Poilblanc, Schuch, and Verstraete, Phys. Rev. B 83, 245134 (2011) Poilblanc, Schuch, Perez-Garcia, IC, arXiv:1202.0947 Schuch, Poilblanc, IC, Perez-Garcia,


slide-1
SLIDE 1

Boundary theories for spins in lattices

  • J. IGNACIO CIRAC

IC, Poilblanc, Schuch, and Verstraete, Phys. Rev. B 83, 245134 (2011) Poilblanc, Schuch, Perez-Garcia, IC, arXiv:1202.0947 Schuch, Poilblanc, IC, Perez-Garcia, arXiv:1203.4816

GGI, Florence, May 24, 2012

slide-2
SLIDE 2

N spins:

States and observables can be written in terms of tensors Expectation values are tensor contractions:

= 〉 Ψ Ψ 〈

j i i i j j i i j j

N N N N

c X c X

, ,..., ,..., ,..., ,...,

1 1 1 1

| |

〉〈 =

j i N N j j i i

i i j j X X

N N

, 1 1 ,..., ,...,

| ,..., ,..., |

1 1

〉 = 〉 Ψ

i N i i

i i c

N

,..., | |

1 ,...,

1

N

i

1

i

N

i

1

i

N

j

1

j

1,..., N

i i

c

1 1

,..., ,...,

N N

j j i i

X

TENSOR NETWORKS

slide-3
SLIDE 3

=

Why? Efficient description:

a b c

N d D

N

i

1

i

1

i

N

i

Rewrite tensors in terms of smaller tensors: STATES: OBSERVABLES: similarly

Tensor network states Tensor product states

  • r

TENSOR NETWORKS

Guiding principle: entanglement

slide-4
SLIDE 4

Multi-scale ENTANGLEMENT renormalization ansatz

TENSOR NETWORKS EXAMPLES

Projected ENTANGLED-pair states

N

i

1

i

11

i

NN

i

1N

i

MERA: G.Vidal PEPS: F.Verstraete, I. Cirac

slide-5
SLIDE 5

TENSOR NETWORKS PEPS

Projected ENTANGLED-pair states

Thermal equilibrium Local interactions Arbitrary dimensions (Hastings) Numerical algorithms

2

N

physically relevant

slide-6
SLIDE 6

TENSOR NETWORKS PEPS

Projected ENTANGLED-pair states

Thermal equilibrium Local interactions Arbitrary dimensions (Hastings) Many-body physics

n

i α β

n

i

A

αβγδ δ γ

2

N

physically relevant

slide-7
SLIDE 7

Projected entangled-pairs:

PEPS

| Ψ〉

Physical spins:

slide-8
SLIDE 8

Projected entangled-pairs:

PEPS

| Ψ〉 |

N ⊗

Φ〉

Physical spins: Auxiliary spins:

slide-9
SLIDE 9

Projected entangled-pairs:

PEPS

| Ψ〉 |

N ⊗

Φ〉

Physical spins: Auxiliary spins:

| , , , |

i

P A i = 〉〈

αβγδ

α β γ δ

P‘s act locally Similar to AKLT construction Contain the information about the state

slide-10
SLIDE 10

Projected entangled-pairs:

PEPS

| Ψ〉 |

N ⊗

Φ〉

Physical spins: Auxiliary spins:

| , , , |

i

P A i = 〉〈

αβγδ

α β γ δ

n

i α β

n

i

A

αβγδ δ γ

slide-11
SLIDE 11

Area law:

PEPS

A

( )

A A

S N ρ

slide-12
SLIDE 12

Area law:

PEPS

slide-13
SLIDE 13

Area law:

PEPS

A

Only the auxiliary particles at the boundary contribute Linear maps P cannot increase entanglement

slide-14
SLIDE 14

Area law:

PEPS

Only the auxiliary particles at the boundary contribute Linear maps P cannot increase entanglement

( )

A A

S N ρ

# degrees of freedom in the bulk scale with the size of boundary

slide-15
SLIDE 15

PEPS

Bulk-boundary correspondence:

A

ρ

Bulk Boundary

U :

A A

U H h

⊗ ⊗∂

isommetry

† A A

U U σ ρ

∂ = A

X

† A A

x UX U

∂ =

A

1h UU =

1H U U =

slide-16
SLIDE 16

Bulk-boundary correspondence:

PEPS

A

H A

e ρ

=

Bulk Boundary

† A A

U U σ ρ

∂ = A

X

† A A

x UX U

∂ =

Expectation values:

† †

tr( ) tr( ) tr( )

A A A A A A A A

X X UX U U U x x ρ ρ σ

∂ ∂ ∂

〈 〉 = = = = 〈 〉

Boundary Hamiltonian:

A

H A

e σ

− ∂ = † A A

H UH U

∂ = A

ρ

Entanglement spectrum:: (

) ( )

A A

H H σ σ

=

The standard ES is exactly the spectrum of the boundary Hamiltonian The boundary Hamltonian has a physical meaning

slide-17
SLIDE 17

PEPS

The boundary operators can be determined using PEPS algorithms Symmetries:

Boundary theory:

A

H A

e σ

− ∂ =

The boundary Hamiltonian reflects properties of the original state | Ψ〉

| |

g

i g

u e

θ

Ψ〉 = Ψ〉 ⇒

† g A g A

U H U H

∂ ∂

=

Criticality: If is the ground state of a GAPPED LOCAL Hamiltonian, then the boundary Hamiltonian is LOCAL

Ψ

Topology: Non-local projector

slide-18
SLIDE 18

PEPS

AKLT model in 2D

Examples:

  • Auxiliary particles s=1/2
  • Symmetry:
  • Finite correlation length

(2) su

A

H∂

is the 1D Heisenberg Hamiltonian ES corresponds to c=1 CFT Kitaev‘s toric code

  • Auxiliary particles s=1/2
  • Symmetry:
  • Finite correlation length
  • Topological

2

Z

A

σ ∂

is a non-local projector ES is flat

2

( ) Z

RVB on a Kagome lattice

  • Auxiliary particles s=1
  • Symmetry:
  • Finite correlation length
  • Topological

A

σ ∂

contains a non-local projector 2

( ) Z (2) su

A

H∂

is a 1D t-J model

slide-19
SLIDE 19

PEPS

spin liquids

Kitaev‘s toric code (square lattice)

2

Z

  • -RVB (RK)

(Kagome lattice) RVB (Kagome lattice)

slide-20
SLIDE 20

PEPS

spin liquids

Kitaev‘s toric code (square lattice) They correspond to the same phase

2

Z

  • -RVB (RK)

(Kagome lattice) RVB (Kagome lattice) local unitary local invertible RVB is ground state of local (FF) Hamiltonian (4-fold degeneracy)

slide-21
SLIDE 21

OUTLINE

How to determine the boundary theory for a PEPS Symmetries Finite correlation length Examples

slide-22
SLIDE 22

Reduced state:

PEPS BOUNDARY THEORY

A B

Combine the tensors of regions A and B

A B

boundary Polar decomposition:

U

A

σ

B

σ V | Ψ〉

Reduced state:

tr (| |)

A B

ρ = Ψ〉〈Ψ

=

† A B A

U U σ σ σ

A A B A

σ σ σ σ

∂ =

slide-23
SLIDE 23

In practice:

PEPS BOUNDARY THEORY

1.- Contract the tensor A with ist complex conjugate

=

2.- Determine

,

A B

σ σ

slide-24
SLIDE 24

Cylinder: A

Exact calculations N ×∞ Reflection symmetry:

A B

σ σ =

2 A A

σ σ

∂ =

with MPS algorithms ∞×∞

PEPS BOUNDARY THEORY

slide-25
SLIDE 25

Gauge :

Different tensors give rise to the same state:

PEPS SYMMETRIES

n

i

A

αβγδ

n

i

Bαβγδ

=

1

X − X

1

Y − Y

Under general conditions, the above is the only possibility

(Perez-Garcia, Sanz, Gonzalez, Wolf, IC, 2009)

slide-26
SLIDE 26

Symmetry :

Global symmetry:

PEPS SYMMETRIES

where the v‘s and w‘s are (projective) representations of the same group

| |

g

i g

U e

θ

Ψ〉 = Ψ〉

must be related by a local Gauge trafo

N g g

U u⊗ =

g

u

=

† g

w

g

v

† g

v

g

w

Boundary operator has the same symmetry:

† A A A g A g

v v σ σ

⊗∂ ⊗∂ ∂ ∂

=

slide-27
SLIDE 27

Cylinder :

PEPS FINITE CORRELATION LENGTH

1,

α 

2,

α

 3,

α

=

A

σ

RG ?

For pure states (MPS): full classification

(Verstraete, IC, Latorre, Rico, Wolf, 2005)

slide-28
SLIDE 28

RG for mixed states MPDO :

PEPS FINITE CORRELATION LENGTH

Boundary theory: entangled states Trace-preserving CPM

, 1 n n

h A

e σ

+

− ∂

∑ = ⊕

Local Hamiltonian Degeneracy and topology

slide-29
SLIDE 29

2D AKLT in a 2-leg ladder/square lattice:

PEPS EXAMPLES

h

N

v

N

v

N 2 S =

spin Deformed AKLT Hamiltonian

, ,

( ) ( ) ( ) ( )

n m n m m n n m

H Q Q P Q Q

< >

= ∆ ∆ ∆ ∆

projector onto S=4 subspace nematic deformation

2

8

( )

z

S

Q e− ∆ ∆ =

Boundary Hamiltonian:

A r r

H d

∂ =∑ ∑ all possible terms with range-r interacctions

Ground state: PEPS with D=2 Symmetry:

(2) / (1) su u

slide-30
SLIDE 30

2D AKLT in a 2-leg ladder/square lattice:

PEPS EXAMPLES

Similar results with other models For AKLT the boundary Hamiltonian is s=1/2 Heisenberg

slide-31
SLIDE 31

Kitaev toric code:

PEPS EXAMPLES

A even

  • dd

P P σ ∂ = ⊕

Non-local operator

Boundary Hamiltonian: trivial (up to the projector) Boundary state Ground state: PEPS with D=2 Degenerate ground state. Gapped. Symmetry: 2

Z

slide-32
SLIDE 32

2D RVB on a Kagome lattice:

PEPS EXAMPLES

Boundary Hamiltonian: t-J model

  • 1. single spin ½ at each edge

Parent Hamiltonian acting on two stars PEPS with D=3

(2) su

1 2 ⊕

representation

  • 2. Three spins ½ at each edge: dimers are orthogonal (related to KR model)
slide-33
SLIDE 33

2D RVB on a Kagome lattice:

PEPS EXAMPLES

  • 1. single spin ½ at each edge
slide-34
SLIDE 34

2D: interpolation RVB-oRVB:

PEPS EXAMPLES

| ( ) θ Ψ 〉 | (0) | RVB Ψ 〉 = 〉 | (1) | oRVB Ψ 〉 = 〉

correlation function fidelity RVB and toric code seem to be in the same phase

slide-35
SLIDE 35

‚Uncle‘ Hamiltonians

PEPS EXAMPLES

An order parameter for gapped phases in 1D Parent Hamiltonian for Laughlin spin state in a lattice

Fernandez, Schuch, Wolf, IC, Perez-Garcia, arXiv:1111.5817 Haegeman, Perez-Garcia, IC, Schuch, arXiv:1201.4174

3 1 2 1 2 1

13

| ( 1 ) (1 ) |

N N N N N N g g h h

  • u

u u u

⊗ ⊗ ⊗ ⊗ ⊗ ⊗

= 〈Ψ ⊗ ⊗ ⊗ ⊗ Ψ〉 F

Anne Nielse, IC, German Sierra, arXiv:1201.3096

poster

slide-36
SLIDE 36

2

N

physically relevant

Thermal equilibrium and local interaction spins can be efficiently described by PEPS Numerical algorithms New perspective

HERE: CONCLUSION:

SUMMARY and CONCLUSIONS

Area law: bulk-boundary correspondence Boundary reflects properties of the bulk: criticality, topology, etc Finite correlation length implies locality of boundary Hamiltonian Locality + symmetries dictate entanglement spectrum

Applicaton: contraction of PEPS is efficient

slide-37
SLIDE 37
slide-38
SLIDE 38

Physical interpretation:

TENSOR NETWORKS PROJECTED ENTANGLED-PAIR STATES

| , , , |

i

P A i = 〉〈

αβγδ

α β γ δ

Why do they provide efficient descriptions?

,

| lim | |

n m

M h H N N M

e e

β β β

ϕ ϕ

− − ⊗ ⊗ →∞

  Ψ 〉 〉 = 〉    

, , n m n m

H h

< >

= ∑

slide-39
SLIDE 39

g

u =

g

v

† g

v

g

w

† g

w Symmetries: Gauge symmetries: =

g

v

† g

v

g

w

† g

w ⇔ =

g

v

g

v

† g

w

g

w

TOPOLOGICAL PHASES

slide-40
SLIDE 40

g

v

g

v

g

v

TOPOLOGICAL PHASES

Wilson strings: =

they can be moved to any column

=

g

v

g

v

† g

w

g

w

g

v

g

v

g

v

slide-41
SLIDE 41

g

v

g

v

g

v

g

v

g

v Ψ ' Ψ

Are locally indistinguishable. Any Hamiltonian for which one is the ground state is degenerate.

TOPOLOGICAL PHASES

Ground state degeneracy:

slide-42
SLIDE 42

1,

α 

2,

α

 3,

α

 1,1

i

1,2

i

1,

i 

2,1

i

1, 1

i

+  1,

h

N

i

1,1

i

1,2

i

2,1

i

1,

i 

1,

h

N

i

/2,1

v

N

i

1

α

1 2 1 2

, ; , i

A

α α β β

i

2

α

1

β

2

β