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
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,
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
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
, ,..., ,..., ,..., ,...,
1 1 1 1
j i N N j j i i
N N
, 1 1 ,..., ,...,
1 1
i N i i
N
1 ,...,
1
N
i
1
i
N
i
1
i
N
j
1
j
1,..., N
i i
1 1
,..., ,...,
N N
j j i i
Why? Efficient description:
a b c
N
i
1
i
1
i
N
i
Tensor network states Tensor product states
N
i
1
i
11
i
NN
i
1N
i
MERA: G.Vidal PEPS: F.Verstraete, I. Cirac
Thermal equilibrium Local interactions Arbitrary dimensions (Hastings) Numerical algorithms
N
physically relevant
Thermal equilibrium Local interactions Arbitrary dimensions (Hastings) Many-body physics
n
i α β
n
i
αβγδ δ γ
N
physically relevant
Physical spins:
N ⊗
Physical spins: Auxiliary spins:
N ⊗
Physical spins: Auxiliary spins:
i
αβγδ
P‘s act locally Similar to AKLT construction Contain the information about the state
N ⊗
Physical spins: Auxiliary spins:
i
αβγδ
n
i α β
n
i
αβγδ δ γ
A A
∂
Only the auxiliary particles at the boundary contribute Linear maps P cannot increase entanglement
Only the auxiliary particles at the boundary contribute Linear maps P cannot increase entanglement
A A
∂
# degrees of freedom in the bulk scale with the size of boundary
A
Bulk Boundary
A A
⊗ ⊗∂
isommetry
† A A
∂ = A
† A A
∂ =
†
†
A
H A
−
Bulk Boundary
† A A
∂ = A
† A A
∂ =
Expectation values:
† †
A A A A A A A A
∂ ∂ ∂
Boundary Hamiltonian:
A
H A
∂
− ∂ = † A A
∂ = A
Entanglement spectrum:: (
A A
∂
The standard ES is exactly the spectrum of the boundary Hamiltonian The boundary Hamltonian has a physical meaning
The boundary operators can be determined using PEPS algorithms Symmetries:
A
H A
∂
− ∂ =
The boundary Hamiltonian reflects properties of the original state | Ψ〉
g
i g
θ
† g A g A
∂ ∂
Criticality: If is the ground state of a GAPPED LOCAL Hamiltonian, then the boundary Hamiltonian is LOCAL
Topology: Non-local projector
AKLT model in 2D
A
is the 1D Heisenberg Hamiltonian ES corresponds to c=1 CFT Kitaev‘s toric code
2
A
is a non-local projector ES is flat
2
RVB on a Kagome lattice
A
contains a non-local projector 2
A
is a 1D t-J model
Kitaev‘s toric code (square lattice)
2
(Kagome lattice) RVB (Kagome lattice)
Kitaev‘s toric code (square lattice) They correspond to the same phase
2
(Kagome lattice) RVB (Kagome lattice) local unitary local invertible RVB is ground state of local (FF) Hamiltonian (4-fold degeneracy)
Combine the tensors of regions A and B
boundary Polar decomposition:
A
B
Reduced state:
A B
† A B A
A A B A
∂ =
1.- Contract the tensor A with ist complex conjugate
2.- Determine
A B
Exact calculations N ×∞ Reflection symmetry:
A B
2 A A
∂ =
with MPS algorithms ∞×∞
Different tensors give rise to the same state:
n
i
αβγδ
n
i
1
1
Under general conditions, the above is the only possibility
(Perez-Garcia, Sanz, Gonzalez, Wolf, IC, 2009)
Global symmetry:
where the v‘s and w‘s are (projective) representations of the same group
g
i g
θ
must be related by a local Gauge trafo
N g g
g
† g
g
† g
g
Boundary operator has the same symmetry:
† A A A g A g
⊗∂ ⊗∂ ∂ ∂
1,
α
2,
α
3,
α
=
A
For pure states (MPS): full classification
(Verstraete, IC, Latorre, Rico, Wolf, 2005)
Boundary theory: entangled states Trace-preserving CPM
, 1 n n
h A
+
− ∂
Local Hamiltonian Degeneracy and topology
h
v
v
spin Deformed AKLT Hamiltonian
, ,
n m n m m n n m
< >
projector onto S=4 subspace nematic deformation
2
8
z
S
Boundary Hamiltonian:
A r r
∂ =∑ ∑ all possible terms with range-r interacctions
Ground state: PEPS with D=2 Symmetry:
Similar results with other models For AKLT the boundary Hamiltonian is s=1/2 Heisenberg
A even
Non-local operator
Boundary Hamiltonian: trivial (up to the projector) Boundary state Ground state: PEPS with D=2 Degenerate ground state. Gapped. Symmetry: 2
Boundary Hamiltonian: t-J model
Parent Hamiltonian acting on two stars PEPS with D=3
1 2 ⊕
representation
correlation function fidelity RVB and toric code seem to be in the same phase
Fernandez, Schuch, Wolf, IC, Perez-Garcia, arXiv:1111.5817 Haegeman, Perez-Garcia, IC, Schuch, arXiv:1201.4174
3 1 2 1 2 1
13
N N N N N N g g h h
⊗ ⊗ ⊗ ⊗ ⊗ ⊗
Anne Nielse, IC, German Sierra, arXiv:1201.3096
N
physically relevant
Thermal equilibrium and local interaction spins can be efficiently described by PEPS Numerical algorithms New perspective
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
i
αβγδ
,
n m
M h H N N M
β β β
− − ⊗ ⊗ →∞
, , n m n m
< >
g
g
† g
g
† g
g
† g
g
† g
g
g
† g
g
g
g
g
they can be moved to any column
g
g
† g
g
g
g
g
g
g
g
g
g
Are locally indistinguishable. Any Hamiltonian for which one is the ground state is degenerate.
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
α α β β
i
2
α
1
β
2
β