Entanglement entropy:
From Field Theory to Condensed Matter
Based on collaboration with:
- J. Cardy, V. Alba, M. Fagotti, E. Tonni....
Entanglement entropy: From Field Theory to Condensed Matter - - PowerPoint PPT Presentation
Entanglement entropy: From Field Theory to Condensed Matter Pasquale Calabrese SISSA-Trieste ERG 2016, Trieste, 23/9/2016 Based on collaboration with: J. Cardy, V. Alba, M. Fagotti, E. Tonni.... Why entanglement entropy? hep-th arXiv
0" 2" 0" 0" 4" 0" 0" 6" 3" 3" 3" 0" 11" 13" 22" 40" 39" 37" 35" 56" 62" 76" 119" 170" 220" 0" 50" 100" 150" 200" 250" 91" 92" 93" 94" 95" 96" 97" 98" 99" 00" 01" 02" 03" 04" 05" 06" 07" 08" 09" 10" 11" 12" 13" 14" 15" Year%
hep-th arXiv preprints with “entanglement” in the title
0" 2" 0" 0" 4" 0" 0" 6" 3" 3" 3" 0" 11" 13" 22" 40" 39" 37" 35" 56" 62" 76" 119" 170" 220" 0" 50" 100" 150" 200" 250" 91" 92" 93" 94" 95" 96" 97" 98" 99" 00" 01" 02" 03" 04" 05" 06" 07" 08" 09" 10" 11" 12" 13" 14" 15" Year%
hep-th arXiv preprints with “entanglement” in the title
PW Anderson 1972
Interacting particles Mott insulators Topological states Unconventional Superconductors Spin-Charge separation
Non-abelian Statistics
Confined phases (QCD)
2N coefficients: too many for a classical PC
si=±
2N coefficients: too many for a classical PC
si=±
Alice can measure only in A, while Bob in the remainder B Alice measures are entangled with Bob’s ones: Schmidt deco |Ψ =
cn|ΨnA|ΨnB cn ⌅ 0,
c2
n = 1
Alice can measure only in A, while Bob in the remainder B Alice measures are entangled with Bob’s ones: Schmidt deco |Ψ =
cn|ΨnA|ΨnB cn ⌅ 0,
c2
n = 1
= -∑ cn ln cn
2 2
[Srednicki ’93 +many more]
[Srednicki ’93 +many more]
Full Hilbert space Area law states
Subir Sachdev
s1,...,sN=±
s1 . . . A[N] sN
±
⇤Φ1(x)|ρA|Φ2(x)⌅ =
PC, J Cardy 2004
The density matrix at temperature β-1
= Z [dφ(x, τ)] Z Y
x
δ(φ(x, 0)φ2(x)) Y
x
δ(φ(x, β)φ1(x)) e−SE
R
The trace sews together the edges along τ = 0 and τ = β to form a cylinder
A = (u, v): ρA sews together only those points x which are not in A, leaving an
PC, J Cardy 2004
A =
n
n→1
A
= cn|u − v|− c
6 (n−1/n)
w → = w−u
w−v ; → z = 1/n⇒ w → z =
w−v
⇥1/n
Trn
A = ⌅Tn(u) ¯
Tn(v)⇧
⇤T T ⌅ ∆Tn = c 12 ⇧ n 1 n ⌃
SA = − lim
n→1
⇤ ⇤nTr⇥n
A = c
3 log ⌅
n
n
PC, J Cardy 2004
doi:10.1038/nature15750
Rajibul Islam1, Ruichao Ma1, Philipp M. Preiss1, M. Eric Tai1, Alexander Lukin1, Matthew Rispoli1 & Markus Greiner1 Entanglement is one of the most intriguing features of quantum mechanics. It describes non-local correlations between quantum objects, and is at the heart of quantum information sciences. Entanglement is now being studied in diverse fields ranging from condensed matter to quantum gravity. However, measuring entanglement remains a challenge. This is especially so in systems of interacting delocalized particles, for which a direct experimental measurement of
spatial entanglement has been elusive. Here, we measure entanglement in such a system of itinerant particles using quantum interference of many-body twins. Making use of our single-site-resolved control of ultracold bosonic atoms in optical lattices, we prepare two identical copies of a many-body state and interfere them. This enables us to directly measure quantum purity, Rényi entanglement entropy, and mutual information. These experiments pave the way for using entanglement to characterize quantum phases and dynamics of strongly correlated many-body systems.
Figure 1 | Bipartite entanglement and partial measurements. A generic pure quantum many-body state has quantum correlations (shown as arrows) between different parts. If the system is divided into two subsystems A and B, the subsystems will be bipartite entangled with each other when there are quantum correlations between them (right column). Only when there is no bipartite entanglement present, the partitioned system | ψAB〉 can be described as a product of subsystem states | ψA〉 and | ψB〉 (left column). A path for measuring the bipartite entanglement emerges from the concept of partial measurements: ignoring all information about subsystem B (indicated as ‘Trace’) will put subsystem A into a statistical mixture, to a degree given by the amount of bipartite entanglement present. Finding ways of measuring the many-body quantum state purity of the system and comparing that of its subsystems would then enable measurements of entanglement. For an entangled state, the subsystems will have less purity than the full system.
Entangled state Product state
B A B Trace Pure Trace Mixed |= |A ⊗ | B |≠ |A ⊗ | B
Figure 2 | Measurement of quantum purity with many-body bosonic interference of quantum twins. a, When two N-particle bosonic systems that are in identical pure quantum states are interfered on a 50%–50% beam splitter, they always produce output states with an even number
a b
Even particle number in Output 2 Even particle number in Output 1 Odd particle number Odd particle number Two identical N-particle bosonic states |N |N Pi Average parity Tr (1 2) Quantum state
Tr (2) Purity 1 2 P1 P2 1 = 2 pi Pi = =
(k) k∏
II 0.5 1 1.5 2 2.5 3 3.5 0.2 0.4 0.6 1 2 3 4 0.5 1 1.5 2 System A size 1 2 3 1 2 Boundary size 1 2 3 0.5 1 Distance, d S2(A) S2(B) S2(A) S2(B) II I U/J 1.2 2.7 39 12 7.2 A B A B d Mutual Information, IAB IAB
a b
Mott insulator Superfluid Adiabatic melt A B A B A B
c
A B A A B A A B B + + + A B I Superfluid Mott Renyi entropy, S2 Renyi entropy, S2 Renyi entropy, S2 Superfluid Mott insulator II Superfluid 100 101 102 U/Jx IAB IAB IAB
A = Disconnected regions: More complex Riemann surface: Rn,2 of genus (n − 1)
[Rn,N has genus (n − 1)(N − 1)] PC, J Cardy, E Tonni 2009/10
Trρn
A = c2 n
„ |u1 − u2||v1 − v2| |u1 − v1||u2 − v2||u1 − v2||u2 − v1| « c
6 (n−1/n)
Fn(x)
„ | − || − x = (u1−v1)(u2−v2)
(u1−u2)(v1−v2) = 4 − point ratio
PC, J Cardy, E Tonni 2009/10
Tr ⇥n
A = c2 n(⌃1⌃2)− c
6 (n− 1 n ) ⇧
{kj}
⇤ ⌃1⌃2 n2r2 ⌅P
j(∆j+∆j)
⇥
n
⌃
j=1
⇤kj
⇤2
C
Γ is an (n − 1) × (n − 1) matrix Γrs = 2i n
n−1
X
k = 1
sin „ π k n « β k
n cos
» 2π k n (r − s) – with βy = Hy(1 − x) Hy(x) , Hy(x) =
2F1(y, 1 − y; 1; x)
Riemann theta function Θ(z|Γ) ≡ X
m ∈ Zn−1
exp ˆ iπ m · Γ · m + 2πim · z ˜
PC, J Cardy, E Tonni 2009/2010
M Fagotti, PC, 2010
The RDM of two intervals is not trivial because of JW string Igloi-Peschel
−
Refael and Moore, Laflorencie, Santachiara, Jacobsen, Saleur ...
Kitaev and Preskill, Levin and Wen, Fradkin and Moore, Schoutens et al....
PC, Essler, Cardy, Ravanini, Franchini,Ercolessi, Alcaraz.... Haldane, Regnault, Read, Ludwig, Bernevig, Poilblanc, Rezayi, Haque.........
1 1 2 3 5
Casini and Huerta, Myers,
Ryu and Takayanagi. Headrick, Maldacena, Myers PC and Cardy, Vidal, Schollwoeck, Kollath, Eisert, Cirac....
Fazio, Amico, Vidal....
Casini and Huerta, Myers,
Ryu and Takayanagi. Headrick, Maldacena, Myers PC and Cardy, Vidal, Schollwoeck, Kollath, Eisert, Cirac....
Fazio, Amico, Vidal....
von Neumann in 1929 posed the question [1003.2133]
t→∞
t→∞
SA(l,t)
vt SA(l,t)=
πct 6 πc ℓ 12
v
PC, Cardy 2005
SA(l,t)
vt SA(l,t)=
πct 6 πc ℓ 12
v
PC, Cardy 2005
SA(l,t)
Fagotti, PC 2008
dk
vmax exists
PC, Cardy 2005
dk
vmax exists
PC, Cardy 2005
Slower particles change entanglement and correlations after t = l/2vmax: large t is driven by slowest particles
~ ~ position time b d = v t a
quench
Mott insulator. a, A 1d ultracold gas of bosonic atoms (black balls) in an optical lattice is initially prepared deep in the Mott-insulating phase with unity filling. The lattice depth is then abruptly lowered, bringing the system out of
pairs emerge at all sites. Each of these pairs consists of a doublon (red ball) and a holon (blue ball) on top of the unity- filling background, which propagate ballistically in opposite
site occupancy builds up at time t between any pair of sites separated by a distance d = vt, where v is the relative velocity
PC, J Cardy 2005
contribution to the entanglement entropy, this implies
SA(t) ⇡ Z
x02A
dx0 Z
x002B
dx00 Z 1
1
dx Z f(p)dp
Z 1 dpf(p)2vp✓(` 2vpt) + ` Z 1 dpf(p)✓(2vpt `) (
term is vanishing for 2 vmax t<l and the entanglement entropy grows linearly with time up to a value linear in l !
PC, J Cardy 2005
S(t) = t Z
2|✏0|t<`
d' 2⇡ 2|✏0|H(cos ∆')+ ` Z
2|✏0|t>`
d' 2⇡ H(cos ∆') (2)
M Fagotti, PC 2008
H(x) = 1 + x 2 log 1 + x 2 1 x 2 log 1 x 2
= 1 cos '(h + h0) + hh0 ✏'✏0
'
.
cos ∆' = ( contains al
SA(l,t)
vt
t→∞
ge s s ,
s e c s s al
T0
global unitary dynamics local thermalization
1 2 3 4 5 6 0.0 0.2 0.4 0.6 0.8 1.0 Observable A P(A)
pure state pure state
quantum quench
1 2 3 4 5 6 0.0 0.2 0.4 0.6 0.8 1.0 Observable A P(A)
T>0
Schematic of thermalization dynamics in closed systems. An isolated quantum system at zero tem- perature can be described by a single pure wavefunction |Ψi. Subsystems of the full quantum state appear pure, as long as the entanglement (indicated by grey lines) between sub- systems is negligible. If suddenly perturbed, the full system evolves unitarily, developing significant entanglement between all parts of the system. While the full system remains in a pure, zero-entropy state, the entropy of entanglement causes the subsystems to equilibrate, and local, thermal mixed states appear to emerge within a globally pure quantum state.
Quench
1
1
Expand and Measure Local and Global Purity
1
1
1 1
Expand and Measure Local Occupation Number
1 2 2 1 2 1 1 2
~ 50 Sites ~ 50 Sites Mott insulator Even Odd
680 nm
Initialize Many-body interference 45 Er 6 Er Global thermal state purity Locally thermal Locally pure Globally pure On-site Statistics Particle number time after quench (ms) 0.2 0.4 0.6 0.8 1 1 2 3 4 5 6 Entropy: -Log(Tr[2]) 100 10-2 4.6 2.3 10 20
A B C
y x Initial state
quench
Purity: Tr[2] 10-1 P(n) 0.2 0.4 0.6 0.8 1 P(n) Particle number 1 2 3 4 5 6 On-site Statistics Many-body purity t=0 ms t=16 ms
Adam M. Kaufman, M. Eric Tai, Alexander Lukin, Matthew Rispoli, Robert Schittko, Philipp M. Preiss, Markus Greiner*
http://science.sciencemag.org/ Downloaded from
Science 353, 794 (2016)
1 2 3 4
Renyi entropy S2
5 10 15 20 1 2 3 4
time (ms) Renyi entropy S2
5 10 15 20
time (ms)
1 2 3 0.5 1 1.5 Subsystem size slope (ms-1)
Starting from a low-entanglement ground state, a global quantum quench leads to the development of large-scale entanglement between all subsystems. We quench a six-site system from the Mott insulating product state (J/U ⌧ 1) with one atom per site to the weakly interacting regime of J/U = 0.64 and measure the dynamics of the entanglement entropy. As it equilibrates, the system acquires local entropy while the full system entropy remains constant and at a value given by measurement imperfections. The dynamics agree with exact numerical simulations with no free parameters (solid lines). Error bars are the standard error of the mean (S.E.M.). For the largest entropies encountered in the three-site system, the large number of populated microstates leads to a significant statistical uncertainty in the entropy, which is reflected in the upper error bar extending to large entropies or being unbounded. Inset: slope of the early time dynamics, extracted with a piecewise linear fit (see Supplementary Material). The dashed line is the mean of these measurements.
Alba & PC, 2016
SY Y = L
∞
X
n=1
Z dλ[ρn,t(λ) ln ρn,t(λ)−ρn,p(λ) ln ρn,p(λ)−ρn,h(λ) ln ρn,h(λ)]
S(t) = X
n
h 2t Z
2|vn|t<`
dvn()sn() + ` Z
2|vn|t>`
dsn() i ,
S(t) = X
n
h 2t Z
2|vn|t<`
dvn()sn() + ` Z
2|vn|t>`
dsn() i ,
2 4 6 8
time t
0.25 0.5 0.75
S’
∆=2 ∆=4 2 4 6 8 10 12 0.2 0.4
S’
ϑ=0 ∆=1 ϑ=0 ∆=2 ϑ=0 ∆=4 ϑ=0 ∆=10 ϑ=π/6 ∆=2 ϑ=π/9 ∆=8 4 8 12 16
0.2 0.4 0.6
S’
ϑ=π/3 ∆=4 ϑ=π/6 ∆=4 ϑ=π/3 ∆=2
(a) Neel (b) Ferromagnet (c) Dimer (MG)
0.5 1 1.5 2
vMt/ℓ
0.1 0.2 0.3 0.4
S/ℓ
tDMRG Conjecture Extrapolations 8
(c)
ℓ=5 - 20
[Vidal-Werner 02]
case is y |e(1)
i ⌥
d H of
d |e(2)
j ⌥
ach part,
i e(2) j |⇥T2|e(1) k e(2) l
i e(2) l
k e(2) j ⌥
⇤T T ⌅ T T Tr|⇥T2| =
i
|i| =
λi>0
i
λi<0
i
Tr(⇥T2)ne =
ne
i
=
|i|ne +
|i|ne , Tr(⇥T2)no =
no
i
=
|i|no
|i|no
ne→1 ln Tr(⇥T2)ne,
no→1 Tr(T2)no = TrT2 = 1
2
2
PC, J Cardy, E Tonni 2012
A = ⌦Tn(u1) ¯
A )n =
A )n = ⌦T 2 n (u2) ¯
n (v2)↵
2
2
2
Tr(T2
A )ne = (⌦Tne/2(u2) ¯
Tne/2(v2)↵)2 = (Trne/2
A2
)2 Tr(T2
A )no = ⌦Tno(u2) ¯
Tno(v2)↵ = Trno
A2 ,
PC, J Cardy, E Tonni 2012
Tr(T2
A )ne = (⌦Tne/2(u2) ¯
Tne/2(v2)↵)2 = (Trne/2
A2
)2 Tr(T2
A )no = ⌦Tno(u2) ¯
Tno(v2)↵ = Trno
A2 ,
A
no has dimension ∆T 2
no = c
ne has dimension ∆T 2
ne = c
A )n = ⌦T 2 n (u2) ¯
n (v2)↵
2
2
1 4 ln sin( ⇥⇤1
L ) sin(⇥⇤2 L )
sin ⇥(⇤1+⇤2)
L
+ cnst
A )n = ↵Tn(1) ¯
n (0)Tn(2)
⌅ ||⇤T2
A || ⌃
✓ 12 1 + 2 ◆ c
4
⇧ E = c 4 ln 12 1 + 2 + cnst
PC, J Cardy, E Tonni 2012