Minimally Entangled Typical Thermal States (METTS)
|
|
- Arnold Flynn
- 6 years ago
- Views:
Transcription
1 Minimally Entangled Typical Thermal States (METTS) Vijay B. Shenoy Centre for Condensed Matter Theory, IISc Bangalore Quantum Condensed Matter Journal Club April 17, / 33
2 Acknowledgements 2 / 33
3 Acknowledgements Funding: DST, DAE 2 / 33
4 Acknowledgements Funding: DST, DAE 2 / 33
5 Acknowledgements Funding: DST, DAE Amal Medhi 2 / 33
6 Papers in Focus 3 / 33
7 Papers in Focus 3 / 33
8 Papers in Focus 3 / 33
9 Papers in Focus 3 / 33
10 Papers in Focus Main source: Talk by M. Stoudenmire at the KSK Discussion Meeting, March 19-21, 2012, MatScience, Chennai 3 / 33
11 Overview Part I: Difficulties with quantum many body simulations 4 / 33
12 Overview Part I: Difficulties with quantum many body simulations Part II: Ideas of entanglement and entanglement entropy 4 / 33
13 Overview Part I: Difficulties with quantum many body simulations Part II: Ideas of entanglement and entanglement entropy Part III: METTS Algorithm 4 / 33
14 Overview Part I: Difficulties with quantum many body simulations Part II: Ideas of entanglement and entanglement entropy Part III: METTS Algorithm Part IV: (Brief) Technical remarks and our plan 4 / 33
15 Overview Part I: Difficulties with quantum many body simulations Part II: Ideas of entanglement and entanglement entropy Part III: METTS Algorithm Part IV: (Brief) Technical remarks and our plan 4 / 33
16 I. Numerical Simulations of Quantum Many Body Systems Difficulties 5 / 33
17 A Simple Model - TFIM 6 / 33
18 A Simple Model - TFIM The transverse field Ising model (TFIM) 6 / 33
19 A Simple Model - TFIM The transverse field Ising model (TFIM) A lattice of N sites with a spin- 1 2 at each site 6 / 33
20 A Simple Model - TFIM The transverse field Ising model (TFIM) A lattice of N sites with a spin- 1 2 at each site Ferro-magnetic Ising interaction between neighbouring spins 6 / 33
21 A Simple Model - TFIM The transverse field Ising model (TFIM) A lattice of N sites with a spin- 1 2 at each site Ferro-magnetic Ising interaction between neighbouring spins A homogeneous transverse magnetic field h 6 / 33
22 A Simple Model - TFIM The transverse field Ising model (TFIM) A lattice of N sites with a spin- 1 2 at each site Ferro-magnetic Ising interaction between neighbouring spins A homogeneous transverse magnetic field h H = S z i Sz j + h ij i S x i 6 / 33
23 A Simple Model - TFIM The transverse field Ising model (TFIM) A lattice of N sites with a spin- 1 2 at each site Ferro-magnetic Ising interaction between neighbouring spins A homogeneous transverse magnetic field h H = S z i Sz j + h ij i S x i Z 2 symmetry (spin rotation by π about the x-axis) 6 / 33
24 A Simple Model - TFIM The transverse field Ising model (TFIM) A lattice of N sites with a spin- 1 2 at each site Ferro-magnetic Ising interaction between neighbouring spins A homogeneous transverse magnetic field h H = S z i Sz j + h ij i S x i Z 2 symmetry (spin rotation by π about the x-axis) Goal: Obtain the h T phase diagram of TFIM numerically (T is temperature) 6 / 33
25 A Simple Model - TFIM The transverse field Ising model (TFIM) A lattice of N sites with a spin- 1 2 at each site Ferro-magnetic Ising interaction between neighbouring spins A homogeneous transverse magnetic field h H = S z i Sz j + h ij i S x i Z 2 symmetry (spin rotation by π about the x-axis) Goal: Obtain the h T phase diagram of TFIM numerically (T is temperature) Size of Hilbert space M = 2 N = e N ln / 33
26 A Simple Model - TFIM The transverse field Ising model (TFIM) A lattice of N sites with a spin- 1 2 at each site Ferro-magnetic Ising interaction between neighbouring spins A homogeneous transverse magnetic field h H = S z i Sz j + h ij i S x i Z 2 symmetry (spin rotation by π about the x-axis) Goal: Obtain the h T phase diagram of TFIM numerically (T is temperature) Size of Hilbert space M = 2 N = e N ln 2...Hilbert space grows exponentially! 6 / 33
27 A Simple Model - TFIM The transverse field Ising model (TFIM) A lattice of N sites with a spin- 1 2 at each site Ferro-magnetic Ising interaction between neighbouring spins A homogeneous transverse magnetic field h H = S z i Sz j + h ij i S x i Z 2 symmetry (spin rotation by π about the x-axis) Goal: Obtain the h T phase diagram of TFIM numerically (T is temperature) Size of Hilbert space M = 2 N = e N ln 2...Hilbert space grows exponentially! Matrix diagonalization difficulty goes as M 3 e3 ln 2N 6 / 33
28 A Simple Model - TFIM The transverse field Ising model (TFIM) A lattice of N sites with a spin- 1 2 at each site Ferro-magnetic Ising interaction between neighbouring spins A homogeneous transverse magnetic field h H = S z i Sz j + h ij i S x i Z 2 symmetry (spin rotation by π about the x-axis) Goal: Obtain the h T phase diagram of TFIM numerically (T is temperature) Size of Hilbert space M = 2 N = e N ln 2...Hilbert space grows exponentially! Matrix diagonalization difficulty goes as M 3 e3 ln 2N Difficulty generic to all quantum problems! 6 / 33
29 TFIM in 1D 7 / 33
30 TFIM in 1D TFIM on a 1D chain is exactly soluble 7 / 33
31 TFIM in 1D TFIM on a 1D chain is exactly soluble At T = 0, there are two phases 7 / 33
32 TFIM in 1D TFIM on a 1D chain is exactly soluble At T = 0, there are two phases For h 1: Z 2 -broken ferromagnetic phase S z 0 7 / 33
33 TFIM in 1D TFIM on a 1D chain is exactly soluble At T = 0, there are two phases For h 1: Z 2 -broken ferromagnetic phase S z 0 For h 1: Z2 symmetric (unique) state S x 0 7 / 33
34 TFIM in 1D TFIM on a 1D chain is exactly soluble At T = 0, there are two phases For h 1: Z 2 -broken ferromagnetic phase S z 0 For h 1: Z2 symmetric (unique) state S x 0 Quantum phase transition at h = / 33
35 TFIM in 1D TFIM on a 1D chain is exactly soluble At T = 0, there are two phases For h 1: Z 2 -broken ferromagnetic phase S z 0 For h 1: Z2 symmetric (unique) state S x 0 Quantum phase transition at h = 1 2 Critical point described by a gapless free fermion theory (CFT with central charge 1/2) 7 / 33
36 TFIM in 1D TFIM on a 1D chain is exactly soluble At T = 0, there are two phases For h 1: Z 2 -broken ferromagnetic phase S z 0 For h 1: Z2 symmetric (unique) state S x 0 Quantum phase transition at h = 1 2 Critical point described by a gapless free fermion theory (CFT with central charge 1/2)... 7 / 33
37 TFIM in 1D TFIM on a 1D chain is exactly soluble At T = 0, there are two phases For h 1: Z 2 -broken ferromagnetic phase S z 0 For h 1: Z2 symmetric (unique) state S x 0 Quantum phase transition at h = 1 2 Critical point described by a gapless free fermion theory (CFT with central charge 1/2)... Strategy: Use TFIM as a tool to understand METTS 7 / 33
38 TFIM in 1D TFIM on a 1D chain is exactly soluble At T = 0, there are two phases For h 1: Z 2 -broken ferromagnetic phase S z 0 For h 1: Z2 symmetric (unique) state S x 0 Quantum phase transition at h = 1 2 Critical point described by a gapless free fermion theory (CFT with central charge 1/2)... Strategy: Use TFIM as a tool to understand METTS All calculations shown here with N = 10 (M = 1024) done on this laptop!... Essentially impossible if N = 14 7 / 33
39 II. All is not lost -Entanglement and Entanglement Entropy 8 / 33
40 Difficulty: A Second Look 9 / 33
41 Difficulty: A Second Look Use the states s 1, s 2,..., s n the product states as the basis (s i = ± 1 2 ) 9 / 33
42 Difficulty: A Second Look Use the states s 1, s 2,..., s n the product states as the basis (s i = ± 1 2 ) A generic state ψ = ψ(s 1, s 2,..., s n ) s 1, s 2,..., s n 9 / 33
43 Difficulty: A Second Look Use the states s 1, s 2,..., s n the product states as the basis (s i = ± 1 2 ) A generic state ψ = ψ(s 1, s 2,..., s n ) s 1, s 2,..., s n...we need 2 N pieces of information to describe a typical state...looks complicated! 9 / 33
44 Difficulty: A Second Look Use the states s 1, s 2,..., s n the product states as the basis (s i = ± 1 2 ) A generic state ψ = ψ(s 1, s 2,..., s n ) s 1, s 2,..., s n...we need 2 N pieces of information to describe a typical state...looks complicated! Question: How complicated is a state ψ? Is there a way to quantify this? 9 / 33
45 Difficulty: A Second Look Use the states s 1, s 2,..., s n the product states as the basis (s i = ± 1 2 ) A generic state ψ = ψ(s 1, s 2,..., s n ) s 1, s 2,..., s n...we need 2 N pieces of information to describe a typical state...looks complicated! Question: How complicated is a state ψ? Is there a way to quantify this? Idea: If the state ψ is made up by gluing states of different parts of the full system then the wave function is not complicated 9 / 33
46 Difficulty: A Second Look Use the states s 1, s 2,..., s n the product states as the basis (s i = ± 1 2 ) A generic state ψ = ψ(s 1, s 2,..., s n ) s 1, s 2,..., s n...we need 2 N pieces of information to describe a typical state...looks complicated! Question: How complicated is a state ψ? Is there a way to quantify this? Idea: If the state ψ is made up by gluing states of different parts of the full system then the wave function is not complicated Need to find a way to see if the state can be split up into pieces that reside on different (spatial) parts of the system 9 / 33
47 Splitting up a state 10 / 33
48 Splitting up a state 10 / 33
49 Splitting up a state Imagine a system broken up in to two parts A and B 10 / 33
50 Splitting up a state Imagine a system broken up in to two parts A and B The product states in A(B) are i A ( j B ) 10 / 33
51 Splitting up a state Imagine a system broken up in to two parts A and B The product states in A(B) are i A ( j B ) The state ψ of the system can written as ψ = ψ(i A, j B ) i A, j B, i A, j B = i A j B i A,j B since i A, j B span the Hilbert space of the full system 10 / 33
52 Splitting up a state Imagine a system broken up in to two parts A and B The product states in A(B) are i A ( j B ) The state ψ of the system can written as ψ = ψ(i A, j B ) i A, j B, i A, j B = i A j B i A,j B since i A, j B span the Hilbert space of the full system Question: Suppose I want to keep information about ψ using i A, what should I keep so that it most faithfully represents ψ? 10 / 33
53 The Reduced Density Matrix 11 / 33
54 The Reduced Density Matrix Most faithfully means: Let O A be an observable that is defined only on sub-system A..ask that ψ O A ψ to be reproduced exactly for every O A with the information that we keep! 11 / 33
55 The Reduced Density Matrix Most faithfully means: Let O A be an observable that is defined only on sub-system A..ask that ψ O A ψ to be reproduced exactly for every O A with the information that we keep! It is easy to show that the information about ψ that we need to keep is the reduced density matrix 11 / 33
56 The Reduced Density Matrix Most faithfully means: Let O A be an observable that is defined only on sub-system A..ask that ψ O A ψ to be reproduced exactly for every O A with the information that we keep! It is easy to show that the information about ψ that we need to keep is the reduced density matrix ρ A = i A,i A ρ A (i A, i A) i A i A, ρ A (i A, i A) = j B ψ (i A, j B)ψ(i A, j B ) which is an operator defined on A by tracing out the states of B in ψ 11 / 33
57 The Reduced Density Matrix Most faithfully means: Let O A be an observable that is defined only on sub-system A..ask that ψ O A ψ to be reproduced exactly for every O A with the information that we keep! It is easy to show that the information about ψ that we need to keep is the reduced density matrix ρ A = i A,i A ρ A (i A, i A) i A i A, ρ A (i A, i A) = j B ψ (i A, j B)ψ(i A, j B ) which is an operator defined on A by tracing out the states of B in ψ For every O A, ψ O A ψ = tr(ρ A O A ) 11 / 33
58 The Reduced Density Matrix Most faithfully means: Let O A be an observable that is defined only on sub-system A..ask that ψ O A ψ to be reproduced exactly for every O A with the information that we keep! It is easy to show that the information about ψ that we need to keep is the reduced density matrix ρ A = i A,i A ρ A (i A, i A) i A i A, ρ A (i A, i A) = j B ψ (i A, j B)ψ(i A, j B ) which is an operator defined on A by tracing out the states of B in ψ For every O A, ψ O A ψ = tr(ρ A O A ) ρ A is non-negative Hermitian 11 / 33
59 The Reduced Density Matrix Most faithfully means: Let O A be an observable that is defined only on sub-system A..ask that ψ O A ψ to be reproduced exactly for every O A with the information that we keep! It is easy to show that the information about ψ that we need to keep is the reduced density matrix ρ A = i A,i A ρ A (i A, i A) i A i A, ρ A (i A, i A) = j B ψ (i A, j B)ψ(i A, j B ) which is an operator defined on A by tracing out the states of B in ψ For every O A, ψ O A ψ = tr(ρ A O A ) ρ A is non-negative Hermitian trρ A = 1 11 / 33
60 Entanglement Spectrum and Entanglement Entropy 12 / 33
61 Entanglement Spectrum and Entanglement Entropy We can diagonalize ρ A = α wα A α A α A (note 0 w α A 1) 12 / 33
62 Entanglement Spectrum and Entanglement Entropy We can diagonalize ρ A = α wα A α A α A (note 0 w α A 1) The spectrum w α A is called the entanglement spectrum 12 / 33
63 Entanglement Spectrum and Entanglement Entropy We can diagonalize ρ A = α wα A α A α A (note 0 w α A 1) The spectrum w α A is called the entanglement spectrum Back to the question: How complicated is ψ? / 33
64 Entanglement Spectrum and Entanglement Entropy We can diagonalize ρ A = α wα A α A α A (note 0 w α A 1) The spectrum w α A is called the entanglement spectrum Back to the question: How complicated is ψ?...if w α A is such that it is non-zero only for, say, α = 0, then ψ can be broken up...in fact, ψ = 0 A φ B...need much less 2 N numbers to describe ψ! 12 / 33
65 Entanglement Spectrum and Entanglement Entropy We can diagonalize ρ A = α wα A α A α A (note 0 w α A 1) The spectrum w α A is called the entanglement spectrum Back to the question: How complicated is ψ?...if w α A is such that it is non-zero only for, say, α = 0, then ψ can be broken up...in fact, ψ = 0 A φ B...need much less 2 N numbers to describe ψ! If there are many w α A 0, there is non trivial entanglement 12 / 33
66 Entanglement Spectrum and Entanglement Entropy We can diagonalize ρ A = α wα A α A α A (note 0 w α A 1) The spectrum w α A is called the entanglement spectrum Back to the question: How complicated is ψ?...if w α A is such that it is non-zero only for, say, α = 0, then ψ can be broken up...in fact, ψ = 0 A φ B...need much less 2 N numbers to describe ψ! If there are many w α A 0, there is non trivial entanglement We can quantify the amount of entanglement by the entanglement entropy S A = tr(ρ A log ρ A ) = α w α A log wα A 12 / 33
67 Entanglement Spectrum and Entanglement Entropy We can diagonalize ρ A = α wα A α A α A (note 0 w α A 1) The spectrum w α A is called the entanglement spectrum Back to the question: How complicated is ψ?...if w α A is such that it is non-zero only for, say, α = 0, then ψ can be broken up...in fact, ψ = 0 A φ B...need much less 2 N numbers to describe ψ! If there are many w α A 0, there is non trivial entanglement We can quantify the amount of entanglement by the entanglement entropy S A = tr(ρ A log ρ A ) = α w α A log wα A Example for the 1D chain shown: 12 / 33
68 Entanglement Spectrum and Entanglement Entropy We can diagonalize ρ A = α wα A α A α A (note 0 w α A 1) The spectrum w α A is called the entanglement spectrum Back to the question: How complicated is ψ?...if w α A is such that it is non-zero only for, say, α = 0, then ψ can be broken up...in fact, ψ = 0 A φ B...need much less 2 N numbers to describe ψ! If there are many w α A 0, there is non trivial entanglement We can quantify the amount of entanglement by the entanglement entropy S A = tr(ρ A log ρ A ) = α w α A log wα A Example for the 1D chain shown: Only one w α A 0 mean S A = 0...unentangled 12 / 33
69 Entanglement Spectrum and Entanglement Entropy We can diagonalize ρ A = α wα A α A α A (note 0 w α A 1) The spectrum w α A is called the entanglement spectrum Back to the question: How complicated is ψ?...if w α A is such that it is non-zero only for, say, α = 0, then ψ can be broken up...in fact, ψ = 0 A φ B...need much less 2 N numbers to describe ψ! If there are many w α A 0, there is non trivial entanglement We can quantify the amount of entanglement by the entanglement entropy S A = tr(ρ A log ρ A ) = α w α A log wα A Example for the 1D chain shown: Only one w α A 0 mean S A = 0...unentangled If all w α A are nearly equal...s A L...badly entangled 12 / 33
70 Entanglement Spectrum and Entanglement Entropy We can diagonalize ρ A = α wα A α A α A (note 0 w α A 1) The spectrum w α A is called the entanglement spectrum Back to the question: How complicated is ψ?...if w α A is such that it is non-zero only for, say, α = 0, then ψ can be broken up...in fact, ψ = 0 A φ B...need much less 2 N numbers to describe ψ! If there are many w α A 0, there is non trivial entanglement We can quantify the amount of entanglement by the entanglement entropy S A = tr(ρ A log ρ A ) = α w α A log wα A Example for the 1D chain shown: Only one w α A 0 mean S A = 0...unentangled If all w α A are nearly equal...s A L...badly entangled 12 / 33
71 How complicated is a state?...and Area Law 13 / 33
72 How complicated is a state?...and Area Law To study this vary the subsystem size L and obtain S A (L) 13 / 33
73 How complicated is a state?...and Area Law To study this vary the subsystem size L and obtain S A (L) If S A (L) grows as L the state is highly entangled! 13 / 33
74 How complicated is a state?...and Area Law To study this vary the subsystem size L and obtain S A (L) If S A (L) grows as L the state is highly entangled! Bad news: 13 / 33
75 How complicated is a state?...and Area Law To study this vary the subsystem size L and obtain S A (L) If S A (L) grows as L the state is highly entangled! Bad news: A typical state in the Hilbert space is highly entangled! 13 / 33
76 How complicated is a state?...and Area Law To study this vary the subsystem size L and obtain S A (L) If S A (L) grows as L the state is highly entangled! Bad news: A typical state in the Hilbert space is highly entangled! Give up? 13 / 33
77 How complicated is a state?...and Area Law To study this vary the subsystem size L and obtain S A (L) If S A (L) grows as L the state is highly entangled! Bad news: A typical state in the Hilbert space is highly entangled! Give up? No: 13 / 33
78 How complicated is a state?...and Area Law To study this vary the subsystem size L and obtain S A (L) If S A (L) grows as L the state is highly entangled! Bad news: A typical state in the Hilbert space is highly entangled! Give up? No: Ground states are special! 13 / 33
79 How complicated is a state?...and Area Law To study this vary the subsystem size L and obtain S A (L) If S A (L) grows as L the state is highly entangled! Bad news: A typical state in the Hilbert space is highly entangled! Give up? No: Ground states are special! Area Law The entanglement entropy of the ground states of nice Hamiltonians in d-dimensions satisfies the area law S A (L) L d 1 d where L is the volume of the subsystem A! Ground states are such that information shared between subsystems is proportional to the interfacial area! 13 / 33
80 How complicated is a state?...and Area Law To study this vary the subsystem size L and obtain S A (L) If S A (L) grows as L the state is highly entangled! Bad news: A typical state in the Hilbert space is highly entangled! Give up? No: Ground states are special! Area Law The entanglement entropy of the ground states of nice Hamiltonians in d-dimensions satisfies the area law S A (L) L d 1 d where L is the volume of the subsystem A! Ground states are such that information shared between subsystems is proportional to the interfacial area! In 1D this means that S A (L) = Const! This is because the interface is a set of two points! 13 / 33
81 How complicated is a state?...and Area Law To study this vary the subsystem size L and obtain S A (L) If S A (L) grows as L the state is highly entangled! Bad news: A typical state in the Hilbert space is highly entangled! Give up? No: Ground states are special! Area Law The entanglement entropy of the ground states of nice Hamiltonians in d-dimensions satisfies the area law S A (L) L d 1 d where L is the volume of the subsystem A! Ground states are such that information shared between subsystems is proportional to the interfacial area! In 1D this means that S A (L) = Const! This is because the interface is a set of two points! Crux: It is much easier to describe ground states! 13 / 33
82 Area Law in Action 14 / 33
83 Area Law in Action Check the area law in TFIM 14 / 33
84 Area Law in Action Check the area law in TFIM Exact diagonalization with N = / 33
85 Area Law in Action Check the area law in TFIM Exact diagonalization with N = S (Entanglement Entropy) L h(transverse magnetic field) 14 / 33
86 Area Law in Action Check the area law in TFIM Exact diagonalization with N = S (Entanglement Entropy) L h(transverse magnetic field) Away from the critical h, area law holds exceedingly well! 14 / 33
87 Area Law in Action Check the area law in TFIM Exact diagonalization with N = S (Entanglement Entropy) L h(transverse magnetic field) Away from the critical h, area law holds exceedingly well! Critical ground states seem to be more complicated! 14 / 33
88 Critical Ground States 15 / 33
89 Critical Ground States Critical ground states are more complicated 15 / 33
90 Critical Ground States Critical ground states are more complicated Good news: they are not highly entangled 15 / 33
91 Critical Ground States Critical ground states are more complicated Good news: they are not highly entangled In a 1D critical ground state (roughly) S A (L) c log(l) 15 / 33
92 Critical Ground States Critical ground states are more complicated Good news: they are not highly entangled In a 1D critical ground state (roughly) S A (L) c log(l) c is the central charge of the CFT describing the critical point! 15 / 33
93 Excited States? 16 / 33
94 Excited States? Does area law hold for excited states? 16 / 33
95 Excited States? Does area law hold for excited states? For h = 0.25, vary L = 2, 3, 4, 5 (N = 10) 16 / 33
96 Excited States? Does area law hold for excited states? For h = 0.25, vary L = 2, 3, 4, 5 (N = 10) S(Ent. Entropy) EigenState Number 16 / 33
97 Excited States? Does area law hold for excited states? For h = 0.25, vary L = 2, 3, 4, 5 (N = 10) S(Ent. Entropy) EigenState Number 16 / 33
98 Excited States? Does area law hold for excited states? For h = 0.25, vary L = 2, 3, 4, 5 (N = 10) S(Ent. Entropy) EigenState Number 16 / 33
99 Excited States? Does area law hold for excited states? For h = 0.25, vary L = 2, 3, 4, 5 (N = 10) S(Ent. Entropy) EigenState Number S ground and highest state (and their neighbours) are unaltered...consistent with area law (if H is nice, so is H) 16 / 33
100 Excited States? Does area law hold for excited states? For h = 0.25, vary L = 2, 3, 4, 5 (N = 10) S(Ent. Entropy) EigenState Number S ground and highest state (and their neighbours) are unaltered...consistent with area law (if H is nice, so is H) But S of a typical excited state grows with L...very bad news for finite temperature simulations! 16 / 33
101 Excited States? Does area law hold for excited states? For h = 0.25, vary L = 2, 3, 4, 5 (N = 10) S(Ent. Entropy) EigenState Number S ground and highest state (and their neighbours) are unaltered...consistent with area law (if H is nice, so is H) But S of a typical excited state grows with L...very bad news for finite temperature simulations! 16 / 33
102 A Question 17 / 33
103 A Question Ground state obeys area law Excited states are highly entangled... Question: Is there a critical excited state? Is there an Anderson Transition somewhere in the spectrum? My prejudice: YES! 17 / 33
104 18 / 33
105 Summary of Part II 18 / 33
106 Summary of Part II Ideas of entanglement, entanglement entropy 18 / 33
107 Summary of Part II Ideas of entanglement, entanglement entropy Area Law: ground states are easy(easier) to describe...the amount of information required to describe a ground state scales as N (linear scaling) 18 / 33
108 Summary of Part II Ideas of entanglement, entanglement entropy Area Law: ground states are easy(easier) to describe...the amount of information required to describe a ground state scales as N (linear scaling) Critical ground states are harder 18 / 33
109 Summary of Part II Ideas of entanglement, entanglement entropy Area Law: ground states are easy(easier) to describe...the amount of information required to describe a ground state scales as N (linear scaling) Critical ground states are harder Excited states...tough! 18 / 33
110 Summary of Part II Ideas of entanglement, entanglement entropy Area Law: ground states are easy(easier) to describe...the amount of information required to describe a ground state scales as N (linear scaling) Critical ground states are harder Excited states...tough! Not a good omen for finite temperature calculations! We may be stuck in the cold confines of the ground state! 18 / 33
111 III. The way to warmth: METTS 19 / 33
112 Finite T Party Line 20 / 33
113 Finite T Party Line Hamiltonian H = n E n Ψ n Ψ n 20 / 33
114 Finite T Party Line Hamiltonian H = n E n Ψ n Ψ n Any observable O (β = 1/T) O = tr(e βh O) tre βh = n e βen Z Ψ n O ψ n Not possible in practise since it is very hard to describe (and find) Ψ n 20 / 33
115 METTS Idea 21 / 33
116 METTS Idea 21 / 33
117 METTS Idea Idea of METTS is to find a (complete) set of states with minimal entanglement that describe the Boltzmann density matrix 21 / 33
118 METTS Idea Idea of METTS is to find a (complete) set of states with minimal entanglement that describe the Boltzmann density matrix Minimum entanglement means that the entanglement entropy of any such state should follow the area law as closely as possible 21 / 33
119 METTS Idea Idea of METTS is to find a (complete) set of states with minimal entanglement that describe the Boltzmann density matrix Minimum entanglement means that the entanglement entropy of any such state should follow the area law as closely as possible White provides a prescription to do this...he calls the state METTS 21 / 33
120 Schrödinger vs. Feynmann...battle in Steven s mind...schrödinger s seems to win! 22 / 33
121 METTS formulation 23 / 33
122 METTS formulation Say we have a complete set of states c 23 / 33
123 METTS formulation Say we have a complete set of states c O = 1 c e βh O c = 1 c e βh 2 Oe βh 2 c = φ(c) O φ(c) Z Z c c c where φ(c) = e βh 2 Z c is an unnormalized state 23 / 33
124 METTS formulation Say we have a complete set of states c O = 1 c e βh O c = 1 c e βh 2 Oe βh 2 c = φ(c) O φ(c) Z Z βh where φ(c) = e 2 Z Normalize φ(c) c c c is an unnormalized state φ(c) φ(c) = P(c) = c e βh c, φ(c) = 1 φ(c) Z P(c) c 23 / 33
125 METTS formulation Say we have a complete set of states c O = 1 c e βh O c = 1 c e βh 2 Oe βh 2 c = φ(c) O φ(c) Z Z βh where φ(c) = e 2 Z Normalize φ(c) c c c is an unnormalized state φ(c) φ(c) = P(c) = c e βh c, φ(c) = 1 φ(c) Z P(c) We can rewrite this as c O = c P(c) φ(c) O φ(c) 23 / 33
126 METTS formulation Say we have a complete set of states c O = 1 c e βh O c = 1 c e βh 2 Oe βh 2 c = φ(c) O φ(c) Z Z βh where φ(c) = e 2 Z Normalize φ(c) c c c is an unnormalized state φ(c) φ(c) = P(c) = c e βh c, φ(c) = 1 φ(c) Z P(c) We can rewrite this as c O = c P(c) φ(c) O φ(c) Key requirement φ(c) s should have as little entanglement as possible! 23 / 33
127 METTS Formulation 24 / 33
128 METTS Formulation METTS idea I: Choose c to be the product states s 1, s 2,..., s N 24 / 33
129 METTS Formulation METTS idea I: Choose c to be the product states s 1, s 2,..., s N These c have zero entanglement...and represent the typical high energy states of the system 24 / 33
130 METTS Formulation METTS idea I: Choose c to be the product states s 1, s 2,..., s N These c have zero entanglement...and represent the typical high energy states of the system The state φ(c) will also be expected to have minimal amount of entanglement...(this can certainly be checked a posteriori) 24 / 33
131 METTS Formulation METTS idea I: Choose c to be the product states s 1, s 2,..., s N These c have zero entanglement...and represent the typical high energy states of the system The state φ(c) will also be expected to have minimal amount of entanglement...(this can certainly be checked a posteriori) Since O = c P(c) φ(c) O φ(c) we now need to device a strategy to sample the state c with the probability distribution P(c) 24 / 33
132 METTS Formulation METTS idea I: Choose c to be the product states s 1, s 2,..., s N These c have zero entanglement...and represent the typical high energy states of the system The state φ(c) will also be expected to have minimal amount of entanglement...(this can certainly be checked a posteriori) Since O = c P(c) φ(c) O φ(c) we now need to device a strategy to sample the state c with the probability distribution P(c) Develop a transition matrix W(c c) such that (Markov stationary condition) W(c c)p(c) = P(c ) c 24 / 33
133 METTS Formulation 25 / 33
134 METTS Formulation METTS idea II: Transition matrix W(c c) = c φ(c) 2 25 / 33
135 METTS Formulation METTS idea II: Transition matrix W(c c) = c φ(c) 2 Stated in words: The next classical state is selected by collapsing the state φ(c) via an S z measurement at every site! 25 / 33
136 METTS Formulation METTS idea II: Transition matrix W(c c) = c φ(c) 2 Stated in words: The next classical state is selected by collapsing the state φ(c) via an S z measurement at every site! This automatically satisfies the detailed balance condition W(c c)p(c) = W(c c )P(c ) and the Markov stationary condition is automatically satisfied! 25 / 33
137 METTS Formulation METTS idea II: Transition matrix W(c c) = c φ(c) 2 Stated in words: The next classical state is selected by collapsing the state φ(c) via an S z measurement at every site! This automatically satisfies the detailed balance condition W(c c)p(c) = W(c c )P(c ) and the Markov stationary condition is automatically satisfied! Note that there is no rejection! 25 / 33
138 METTS Formulation METTS idea II: Transition matrix W(c c) = c φ(c) 2 Stated in words: The next classical state is selected by collapsing the state φ(c) via an S z measurement at every site! This automatically satisfies the detailed balance condition W(c c)p(c) = W(c c )P(c ) and the Markov stationary condition is automatically satisfied! Note that there is no rejection!?scheme fails for the classical Ising model! Product states should not diagonalize the Hamiltonian 25 / 33
139 METTS Algorithm 26 / 33
140 METTS Algorithm In a picture (Schollwock 2009) 26 / 33
141 METTS Algorithm In a picture In Fortran (Schollwock 2009) 26 / 33
142 TIFM with METTS 27 / 33
143 TIFM with METTS METTS simulations with h = / 33
144 TIFM with METTS METTS simulations with h = 0.6 e βh 2 c is obtained exactly (cheating) 27 / 33
145 TIFM with METTS METTS simulations with h = 0.6 e βh 2 c is obtained exactly (cheating) Exact METTS Energy T (Temperature) 27 / 33
146 TIFM with METTS METTS simulations with h = 0.6 e βh 2 c is obtained exactly (cheating) Exact METTS Energy T (Temperature) Algorithm works remarkably well (even with a few hundred steps)! 27 / 33
147 M of METTS How entangled are the METTS? 0.24 S (Entanglement Entropy of METTS) L h= MTTES Step Remarkably little entanglement (lower than the ground state even!) 28 / 33
148 Part III Summary: METTS 29 / 33
149 Part III Summary: METTS Need to use minimally entangled states to obtain partition function trace 29 / 33
150 Part III Summary: METTS Need to use minimally entangled states to obtain partition function trace Minimally entangled states are φ(c) e βh 2 c where c s are the product state 29 / 33
151 Part III Summary: METTS Need to use minimally entangled states to obtain partition function trace Minimally entangled states are φ(c) e βh 2 c where c s are the product state The induced probability distribution on c is sampled by a collapsing scheme 29 / 33
152 Part III Summary: METTS Need to use minimally entangled states to obtain partition function trace Minimally entangled states are φ(c) e βh 2 c where c s are the product state The induced probability distribution on c is sampled by a collapsing scheme The most difficult step: Calculation of e βh 2 c 29 / 33
153 IV. Brief Remarks for Experts and Plan for Future Work 30 / 33
154 MPS Implementation of METTS Matrix product states (MPS) play a crucial role in practical implementation of METTS Operators can also be represented a matrix products Calculation of e βh 2 c involves Breaking up of H = n H n...h n acts locally say at most on two sites Trotter decompose e τh = n e τhn (approximation) Each e τh n has a simple MPO representation Apply e τh sufficient number of times to generate e βh 2 Scales as Nm 3 β where m is the size of matrices in the matrix products 31 / 33
155 Plan for Future Work Strongly correlated fermions Work in grand canonical ensemble Product states are site basis with four states per site Develop efficient imaginary time evolution strategies Do high temperature calculations and check with virial expansion results 32 / 33
156 Summary Part II Summary: Entanglement ideas 33 / 33
157 Summary Part II Summary: Entanglement ideas Ideas of entanglement, entanglement entropy 33 / 33
158 Summary Part II Summary: Entanglement ideas Ideas of entanglement, entanglement entropy Area Law: ground states are easy(easier) to describe...the amount of information required to describe a ground state scales as N 33 / 33
159 Summary Part II Summary: Entanglement ideas Ideas of entanglement, entanglement entropy Area Law: ground states are easy(easier) to describe...the amount of information required to describe a ground state scales as N Critical ground states are harder 33 / 33
160 Summary Part II Summary: Entanglement ideas Ideas of entanglement, entanglement entropy Area Law: ground states are easy(easier) to describe...the amount of information required to describe a ground state scales as N Critical ground states are harder Excited states...tough! 33 / 33
161 Summary Part II Summary: Entanglement ideas Ideas of entanglement, entanglement entropy Area Law: ground states are easy(easier) to describe...the amount of information required to describe a ground state scales as N Critical ground states are harder Excited states...tough! Part III Summary: METTS Need to use minimally entangled states to obtain partition function trace Minimally entangled states are φ(c) e βh 2 c where c s are the product state The induced probability distribution on c is sampled by a collapsing scheme Lots of interesting things to do! 33 / 33
Classical Monte Carlo Simulations
Classical Monte Carlo Simulations Hyejin Ju April 17, 2012 1 Introduction Why do we need numerics? One of the main goals of condensed matter is to compute expectation values O = 1 Z Tr{O e βĥ} (1) and
More informationTensor network simulations of strongly correlated quantum systems
CENTRE FOR QUANTUM TECHNOLOGIES NATIONAL UNIVERSITY OF SINGAPORE AND CLARENDON LABORATORY UNIVERSITY OF OXFORD Tensor network simulations of strongly correlated quantum systems Stephen Clark LXXT[[[GSQPEFS\EGYOEGXMZMXMIWUYERXYQGSYVWI
More information3 Symmetry Protected Topological Phase
Physics 3b Lecture 16 Caltech, 05/30/18 3 Symmetry Protected Topological Phase 3.1 Breakdown of noninteracting SPT phases with interaction Building on our previous discussion of the Majorana chain and
More informationPHY305: Notes on Entanglement and the Density Matrix
PHY305: Notes on Entanglement and the Density Matrix Here follows a short summary of the definitions of qubits, EPR states, entanglement, the density matrix, pure states, mixed states, measurement, and
More information4 Matrix product states
Physics 3b Lecture 5 Caltech, 05//7 4 Matrix product states Matrix product state (MPS) is a highly useful tool in the study of interacting quantum systems in one dimension, both analytically and numerically.
More informationThe 1+1-dimensional Ising model
Chapter 4 The 1+1-dimensional Ising model The 1+1-dimensional Ising model is one of the most important models in statistical mechanics. It is an interacting system, and behaves accordingly. Yet for a variety
More informationEnsembles and incomplete information
p. 1/32 Ensembles and incomplete information So far in this course, we have described quantum systems by states that are normalized vectors in a complex Hilbert space. This works so long as (a) the system
More informationScaling analysis of snapshot spectra in the world-line quantum Monte Carlo for the transverse-field Ising chain
TNSAA 2018-2019 Dec. 3-6, 2018, Kobe, Japan Scaling analysis of snapshot spectra in the world-line quantum Monte Carlo for the transverse-field Ising chain Kouichi Seki, Kouichi Okunishi Niigata University,
More informationTypical quantum states at finite temperature
Typical quantum states at finite temperature How should one think about typical quantum states at finite temperature? Density Matrices versus pure states Why eigenstates are not typical Measuring the heat
More informationEntanglement in Many-Body Fermion Systems
Entanglement in Many-Body Fermion Systems Michelle Storms 1, 2 1 Department of Physics, University of California Davis, CA 95616, USA 2 Department of Physics and Astronomy, Ohio Wesleyan University, Delaware,
More informationEfficient time evolution of one-dimensional quantum systems
Efficient time evolution of one-dimensional quantum systems Frank Pollmann Max-Planck-Institut für komplexer Systeme, Dresden, Germany Sep. 5, 2012 Hsinchu Problems we will address... Finding ground states
More informationSPACETIME FROM ENTANGLEMENT - journal club notes -
SPACETIME FROM ENTANGLEMENT - journal club notes - Chris Heinrich 1 Outline 1. Introduction Big picture: Want a quantum theory of gravity Best understanding of quantum gravity so far arises through AdS/CFT
More informationMatrix Product States
Matrix Product States Ian McCulloch University of Queensland Centre for Engineered Quantum Systems 28 August 2017 Hilbert space (Hilbert) space is big. Really big. You just won t believe how vastly, hugely,
More informationStudying quantum systems using a quantum computer Ph/CS 219, 23 February 2009
Studying quantum systems using a quantum computer Ph/CS 219, 23 February 2009 Estimating energy eigenvalues and preparing energy eigenstates Ph/CS 219, 2 March 2009 We have argued that a quantum
More informationLecture 6: Fluctuation-Dissipation theorem and introduction to systems of interest
Lecture 6: Fluctuation-Dissipation theorem and introduction to systems of interest In the last lecture, we have discussed how one can describe the response of a well-equilibriated macroscopic system to
More informationH ψ = E ψ. Introduction to Exact Diagonalization. Andreas Läuchli, New states of quantum matter MPI für Physik komplexer Systeme - Dresden
H ψ = E ψ Introduction to Exact Diagonalization Andreas Läuchli, New states of quantum matter MPI für Physik komplexer Systeme - Dresden http://www.pks.mpg.de/~aml laeuchli@comp-phys.org Simulations of
More informationQuantum Many Body Systems and Tensor Networks
Quantum Many Body Systems and Tensor Networks Aditya Jain International Institute of Information Technology, Hyderabad aditya.jain@research.iiit.ac.in July 30, 2015 Aditya Jain (IIIT-H) Quantum Hamiltonian
More informationSimulating Quantum Systems through Matrix Product States. Laura Foini SISSA Journal Club
Simulating Quantum Systems through Matrix Product States Laura Foini SISSA Journal Club 15-04-2010 Motivations Theoretical interest in Matrix Product States Wide spectrum of their numerical applications
More informationUniversal phase transitions in Topological lattice models
Universal phase transitions in Topological lattice models F. J. Burnell Collaborators: J. Slingerland S. H. Simon September 2, 2010 Overview Matter: classified by orders Symmetry Breaking (Ferromagnet)
More informationMagnets, 1D quantum system, and quantum Phase transitions
134 Phys620.nb 10 Magnets, 1D quantum system, and quantum Phase transitions In 1D, fermions can be mapped into bosons, and vice versa. 10.1. magnetization and frustrated magnets (in any dimensions) Consider
More informationMany-Body Localization. Geoffrey Ji
Many-Body Localization Geoffrey Ji Outline Aside: Quantum thermalization; ETH Single-particle (Anderson) localization Many-body localization Some phenomenology (l-bit model) Numerics & Experiments Thermalization
More informationS.K. Saikin May 22, Lecture 13
S.K. Saikin May, 007 13 Decoherence I Lecture 13 A physical qubit is never isolated from its environment completely. As a trivial example, as in the case of a solid state qubit implementation, the physical
More informationFirst Problem Set for Physics 847 (Statistical Physics II)
First Problem Set for Physics 847 (Statistical Physics II) Important dates: Feb 0 0:30am-:8pm midterm exam, Mar 6 9:30am-:8am final exam Due date: Tuesday, Jan 3. Review 0 points Let us start by reviewing
More informationPhysics 239/139 Spring 2018 Assignment 2 Solutions
University of California at San Diego Department of Physics Prof. John McGreevy Physics 39/139 Spring 018 Assignment Solutions Due 1:30pm Monday, April 16, 018 1. Classical circuits brain-warmer. (a) Show
More informationHolographic entanglement entropy
Holographic entanglement entropy Mohsen Alishahiha School of physics, Institute for Research in Fundamental Sciences (IPM) 21th Spring Physics Conference, 1393 1 Plan of the talk Entanglement entropy Holography
More informationAn introduction to tensornetwork
An introduction to tensornetwork states and MERA Sissa Journal Club Andrea De Luca 29/01/2010 A typical problem We are given: A lattice with N sites On each site a C d hilbert space A quantum hamiltonian
More informationLecture 11: Long-wavelength expansion in the Neel state Energetic terms
Lecture 11: Long-wavelength expansion in the Neel state Energetic terms In the last class we derived the low energy effective Hamiltonian for a Mott insulator. This derivation is an example of the kind
More informationTemperature Dependence of Entanglement Negativity in Lattice Models: Area Laws and Sudden Death
Temperature Dependence of Entanglement egativity in attice Models: Area aws and Sudden Death icholas E. Sherman Department of Physics, University of California Davis, CA 9566, USA (Dated: ovember 4, 5)
More information1 Quantum field theory and Green s function
1 Quantum field theory and Green s function Condensed matter physics studies systems with large numbers of identical particles (e.g. electrons, phonons, photons) at finite temperature. Quantum field theory
More informationJournal Club: Brief Introduction to Tensor Network
Journal Club: Brief Introduction to Tensor Network Wei-Han Hsiao a a The University of Chicago E-mail: weihanhsiao@uchicago.edu Abstract: This note summarizes the talk given on March 8th 2016 which was
More informationarxiv:quant-ph/ v2 24 Dec 2003
Quantum Entanglement in Heisenberg Antiferromagnets V. Subrahmanyam Department of Physics, Indian Institute of Technology, Kanpur, India. arxiv:quant-ph/0309004 v2 24 Dec 2003 Entanglement sharing among
More informationFabian Essler (Oxford)
Quantum Quenches in Integrable Models Fabian Essler (Oxford) Collaborators: J.S. Caux (Amsterdam), M. Fagotti (Oxford) earlier work: P. Calabrese (Pisa), S. Evangelisti (Oxford), D. Schuricht (Aachen)
More informationQuantum Annealing in spin glasses and quantum computing Anders W Sandvik, Boston University
PY502, Computational Physics, December 12, 2017 Quantum Annealing in spin glasses and quantum computing Anders W Sandvik, Boston University Advancing Research in Basic Science and Mathematics Example:
More informationDecoherence and Thermalization of Quantum Spin Systems
Copyright 2011 American Scientific Publishers All rights reserved Printed in the United States of America Journal of Computational and Theoretical Nanoscience Vol. 8, 1 23, 2011 Decoherence and Thermalization
More informationarxiv:quant-ph/ v2 23 Mar 2001
A PHYSICAL EXPLANATION FOR THE TILDE SYSTEM IN THERMO FIELD DYNAMICS DONG MI, HE-SHAN SONG Department of Physics, Dalian University of Technology, Dalian 116024, P.R.China arxiv:quant-ph/0102127v2 23 Mar
More informationQuantum mechanics in one hour
Chapter 2 Quantum mechanics in one hour 2.1 Introduction The purpose of this chapter is to refresh your knowledge of quantum mechanics and to establish notation. Depending on your background you might
More information8.334: Statistical Mechanics II Problem Set # 4 Due: 4/9/14 Transfer Matrices & Position space renormalization
8.334: Statistical Mechanics II Problem Set # 4 Due: 4/9/14 Transfer Matrices & Position space renormalization This problem set is partly intended to introduce the transfer matrix method, which is used
More informationTensor network methods in condensed matter physics. ISSP, University of Tokyo, Tsuyoshi Okubo
Tensor network methods in condensed matter physics ISSP, University of Tokyo, Tsuyoshi Okubo Contents Possible target of tensor network methods! Tensor network methods! Tensor network states as ground
More informationStatistical Thermodynamics Solution Exercise 8 HS Solution Exercise 8
Statistical Thermodynamics Solution Exercise 8 HS 05 Solution Exercise 8 Problem : Paramagnetism - Brillouin function a According to the equation for the energy of a magnetic dipole in an external magnetic
More informationDensity Matrices. Chapter Introduction
Chapter 15 Density Matrices 15.1 Introduction Density matrices are employed in quantum mechanics to give a partial description of a quantum system, one from which certain details have been omitted. For
More informationThermalization and Revivals after a Quantum Quench in a Finite System
after a Quantum Quench in a Finite System John Cardy University of Oxford EPSRC Workshop Nottingham May 2014 arxiv:1403.3040, to appear in PRL (Global) Quantum Quench prepare an extended system (in thermodynamic
More informationQuantum Information Types
qitd181 Quantum Information Types Robert B. Griffiths Version of 6 February 2012 References: R. B. Griffiths, Types of Quantum Information, Phys. Rev. A 76 (2007) 062320; arxiv:0707.3752 Contents 1 Introduction
More informationAdvanced Computation for Complex Materials
Advanced Computation for Complex Materials Computational Progress is brainpower limited, not machine limited Algorithms Physics Major progress in algorithms Quantum Monte Carlo Density Matrix Renormalization
More informationDiffusion Monte Carlo
Diffusion Monte Carlo Notes for Boulder Summer School 2010 Bryan Clark July 22, 2010 Diffusion Monte Carlo The big idea: VMC is a useful technique, but often we want to sample observables of the true ground
More informationEntanglement Entropy in Extended Quantum Systems
Entanglement Entropy in Extended Quantum Systems John Cardy University of Oxford STATPHYS 23 Genoa Outline A. Universal properties of entanglement entropy near quantum critical points B. Behaviour of entanglement
More informationQuantum Quench in Conformal Field Theory from a General Short-Ranged State
from a General Short-Ranged State John Cardy University of Oxford GGI, Florence, May 2012 (Global) Quantum Quench prepare an extended system at time t = 0 in a (translationally invariant) pure state ψ
More informationExam TFY4230 Statistical Physics kl Wednesday 01. June 2016
TFY423 1. June 216 Side 1 av 5 Exam TFY423 Statistical Physics l 9. - 13. Wednesday 1. June 216 Problem 1. Ising ring (Points: 1+1+1 = 3) A system of Ising spins σ i = ±1 on a ring with periodic boundary
More informationPhysics 127b: Statistical Mechanics. Lecture 2: Dense Gas and the Liquid State. Mayer Cluster Expansion
Physics 27b: Statistical Mechanics Lecture 2: Dense Gas and the Liquid State Mayer Cluster Expansion This is a method to calculate the higher order terms in the virial expansion. It introduces some general
More informationIn-class exercises. Day 1
Physics 4488/6562: Statistical Mechanics http://www.physics.cornell.edu/sethna/teaching/562/ Material for Week 8 Exercises due Mon March 19 Last correction at March 5, 2018, 8:48 am c 2017, James Sethna,
More informationScale invariance on the lattice
Coogee'16 Sydney Quantum Information Theory Workshop Feb 2 nd - 5 th, 2016 Scale invariance on the lattice Guifre Vidal Coogee'16 Sydney Quantum Information Theory Workshop Feb 2 nd - 5 th, 2016 Scale
More informationThe density matrix renormalization group and tensor network methods
The density matrix renormalization group and tensor network methods Outline Steve White Exploiting the low entanglement of ground states Matrix product states and DMRG 1D 2D Tensor network states Some
More informationMarkov Chain Monte Carlo The Metropolis-Hastings Algorithm
Markov Chain Monte Carlo The Metropolis-Hastings Algorithm Anthony Trubiano April 11th, 2018 1 Introduction Markov Chain Monte Carlo (MCMC) methods are a class of algorithms for sampling from a probability
More informationStatistical Mechanics
Franz Schwabl Statistical Mechanics Translated by William Brewer Second Edition With 202 Figures, 26 Tables, and 195 Problems 4u Springer Table of Contents 1. Basic Principles 1 1.1 Introduction 1 1.2
More informationTime-dependent DMRG:
The time-dependent DMRG and its applications Adrian Feiguin Time-dependent DMRG: ^ ^ ih Ψ( t) = 0 t t [ H ( t) E ] Ψ( )... In a truncated basis: t=3 τ t=4 τ t=5τ t=2 τ t= τ t=0 Hilbert space S.R.White
More informationQuantum quenches in 2D with chain array matrix product states
Quantum quenches in 2D with chain array matrix product states Andrew J. A. James University College London Robert M. Konik Brookhaven National Laboratory arxiv:1504.00237 Outline MPS for many body systems
More informationCompression and entanglement, entanglement transformations
PHYSICS 491: Symmetry and Quantum Information April 27, 2017 Compression and entanglement, entanglement transformations Lecture 8 Michael Walter, Stanford University These lecture notes are not proof-read
More informationQuantum field theory and Green s function
1 Quantum field theory and Green s function Condensed matter physics studies systems with large numbers of identical particles (e.g. electrons, phonons, photons) at finite temperature. Quantum field theory
More informationEntanglement in Topological Phases
Entanglement in Topological Phases Dylan Liu August 31, 2012 Abstract In this report, the research conducted on entanglement in topological phases is detailed and summarized. This includes background developed
More informationAn Inverse Mass Expansion for Entanglement Entropy. Free Massive Scalar Field Theory
in Free Massive Scalar Field Theory NCSR Demokritos National Technical University of Athens based on arxiv:1711.02618 [hep-th] in collaboration with Dimitris Katsinis March 28 2018 Entanglement and Entanglement
More informationQuantum and classical annealing in spin glasses and quantum computing. Anders W Sandvik, Boston University
NATIONAL TAIWAN UNIVERSITY, COLLOQUIUM, MARCH 10, 2015 Quantum and classical annealing in spin glasses and quantum computing Anders W Sandvik, Boston University Cheng-Wei Liu (BU) Anatoli Polkovnikov (BU)
More informationPhase Transitions and Critical Behavior:
II Phase Transitions and Critical Behavior: A. Phenomenology (ibid., Chapter 10) B. mean field theory (ibid., Chapter 11) C. Failure of MFT D. Phenomenology Again (ibid., Chapter 12) // Windsor Lectures
More informationNumerical diagonalization studies of quantum spin chains
PY 502, Computational Physics, Fall 2016 Anders W. Sandvik, Boston University Numerical diagonalization studies of quantum spin chains Introduction to computational studies of spin chains Using basis states
More informationSpin liquids in frustrated magnets
May 20, 2010 Contents 1 Frustration 2 3 4 Exotic excitations 5 Frustration The presence of competing forces that cannot be simultaneously satisfied. Heisenberg-Hamiltonian H = 1 J ij S i S j 2 ij The ground
More informationMany-Body Fermion Density Matrix: Operator-Based Truncation Scheme
Many-Body Fermion Density Matrix: Operator-Based Truncation Scheme SIEW-ANN CHEONG and C. L. HENLEY, LASSP, Cornell U March 25, 2004 Support: NSF grants DMR-9981744, DMR-0079992 The Big Picture GOAL Ground
More informationQuantum Phase Transition
Quantum Phase Transition Guojun Zhu Department of Physics, University of Illinois at Urbana-Champaign, Urbana IL 61801, U.S.A. (Dated: May 5, 2002) A quantum system can undergo a continuous phase transition
More informationRenormalization of Tensor- Network States Tao Xiang
Renormalization of Tensor- Network States Tao Xiang Institute of Physics/Institute of Theoretical Physics Chinese Academy of Sciences txiang@iphy.ac.cn Physical Background: characteristic energy scales
More informationSupersymmetry breaking and Nambu-Goldstone fermions in lattice models
YKIS2016@YITP (2016/6/15) Supersymmetry breaking and Nambu-Goldstone fermions in lattice models Hosho Katsura (Department of Physics, UTokyo) Collaborators: Yu Nakayama (IPMU Rikkyo) Noriaki Sannomiya
More informationRenormalization Group: non perturbative aspects and applications in statistical and solid state physics.
Renormalization Group: non perturbative aspects and applications in statistical and solid state physics. Bertrand Delamotte Saclay, march 3, 2009 Introduction Field theory: - infinitely many degrees of
More informationEntanglement Dynamics for the Quantum Disordered XY Chain
Entanglement Dynamics for the Quantum Disordered XY Chain Houssam Abdul-Rahman Joint with: B. Nachtergaele, R. Sims, G. Stolz AMS Southeastern Sectional Meeting University of Georgia March 6, 2016 Houssam
More informationin-medium pair wave functions the Cooper pair wave function the superconducting order parameter anomalous averages of the field operators
(by A. A. Shanenko) in-medium wave functions in-medium pair-wave functions and spatial pair particle correlations momentum condensation and ODLRO (off-diagonal long range order) U(1) symmetry breaking
More informationFloquet Topological Insulators and Majorana Modes
Floquet Topological Insulators and Majorana Modes Manisha Thakurathi Journal Club Centre for High Energy Physics IISc Bangalore January 17, 2013 References Floquet Topological Insulators by J. Cayssol
More informationExcursion: MPS & DMRG
Excursion: MPS & DMRG Johannes.Schachenmayer@gmail.com Acronyms for: - Matrix product states - Density matrix renormalization group Numerical methods for simulations of time dynamics of large 1D quantum
More informationQuantum spin systems - models and computational methods
Summer School on Computational Statistical Physics August 4-11, 2010, NCCU, Taipei, Taiwan Quantum spin systems - models and computational methods Anders W. Sandvik, Boston University Lecture outline Introduction
More informationMBL THROUGH MPS. Bryan Clark and David Pekker. arxiv: Perimeter - Urbana. Friday, October 31, 14
MBL THROUGH MPS Bryan Clark and David Pekker arxiv:1410.2224 Perimeter - Urbana Many Body Localization... It has been suggested that interactions + (strong) disorder LX H = [h i Si z X + J Ŝ i Ŝi+1] i=1
More informationObservation of topological phenomena in a programmable lattice of 1800 superconducting qubits
Observation of topological phenomena in a programmable lattice of 18 superconducting qubits Andrew D. King Qubits North America 218 Nature 56 456 46, 218 Interdisciplinary teamwork Theory Simulation QA
More informationarxiv: v3 [quant-ph] 23 Jun 2011
Feasibility of self-correcting quantum memory and thermal stability of topological order Beni Yoshida Center for Theoretical Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139,
More informationQuantum wires, orthogonal polynomials and Diophantine approximation
Quantum wires, orthogonal polynomials and Diophantine approximation Introduction Quantum Mechanics (QM) is a linear theory Main idea behind Quantum Information (QI): use the superposition principle of
More informationThe general reason for approach to thermal equilibrium of macroscopic quantu
The general reason for approach to thermal equilibrium of macroscopic quantum systems 10 May 2011 Joint work with S. Goldstein, J. L. Lebowitz, C. Mastrodonato, and N. Zanghì. Claim macroscopic quantum
More informationUnusual ordered phases of magnetized frustrated antiferromagnets
Unusual ordered phases of magnetized frustrated antiferromagnets Credit: Francis Pratt / ISIS / STFC Oleg Starykh University of Utah Leon Balents and Andrey Chubukov Novel states in correlated condensed
More informationLanguage of Quantum Mechanics
Language of Quantum Mechanics System s Hilbert space. Pure and mixed states A quantum-mechanical system 1 is characterized by a complete inner-product space that is, Hilbert space with either finite or
More informationTypicality paradigm in Quantum Statistical Thermodynamics Barbara Fresch, Giorgio Moro Dipartimento Scienze Chimiche Università di Padova
Typicality paradigm in Quantum Statistical Thermodynamics Barbara Fresch, Giorgio Moro Dipartimento Scienze Chimiche Università di Padova Outline 1) The framework: microcanonical statistics versus the
More informationAshvin Vishwanath UC Berkeley
TOPOLOGY + LOCALIZATION: QUANTUM COHERENCE IN HOT MATTER Ashvin Vishwanath UC Berkeley arxiv:1307.4092 (to appear in Nature Comm.) Thanks to David Huse for inspiring discussions Yasaman Bahri (Berkeley)
More informationPolynomial-time classical simulation of quantum ferromagnets
Polynomial-time classical simulation of quantum ferromagnets Sergey Bravyi David Gosset IBM PRL 119, 100503 (2017) arxiv:1612.05602 Quantum Monte Carlo: a powerful suite of probabilistic classical simulation
More informationQuantum annealing by ferromagnetic interaction with the mean-field scheme
Quantum annealing by ferromagnetic interaction with the mean-field scheme Sei Suzuki and Hidetoshi Nishimori Department of Physics, Tokyo Institute of Technology, Oh-okayama, Meguro, Tokyo 152-8551, Japan
More information19 Entanglement Entropy in Quantum Field Theory
19 Entanglement Entropy in Quantum Field Theory So far we have discussed entanglement in ordinary quantum mechanics, where the Hilbert space of a finite region is finite dimensional. Now we will discuss
More informationIntroduction to Heisenberg model. Javier Junquera
Introduction to Heisenberg model Javier Junquera Most important reference followed in this lecture Magnetism in Condensed Matter Physics Stephen Blundell Oxford Master Series in Condensed Matter Physics
More informationComplexity of the quantum adiabatic algorithm
Complexity of the quantum adiabatic algorithm Peter Young e-mail:peter@physics.ucsc.edu Collaborators: S. Knysh and V. N. Smelyanskiy Colloquium at Princeton, September 24, 2009 p.1 Introduction What is
More informationMP 472 Quantum Information and Computation
MP 472 Quantum Information and Computation http://www.thphys.may.ie/staff/jvala/mp472.htm Outline Open quantum systems The density operator ensemble of quantum states general properties the reduced density
More informationQuantum entanglement, it s entropy, and why we calculate it
Quantum entanglement, it s entropy, and why we calculate it Piotr Witkowski Max Planck Institute for Physics 14.7 2016 Munich 1 What is entanglement? 2 Quantifying entanglement the entropy 3 The (very)
More informationarxiv:quant-ph/ v5 10 Feb 2003
Quantum entanglement of identical particles Yu Shi Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, United Kingdom and Theory of
More informationDensity matrix. (3) where in the last relation I have used the fact that Ĥ is Hermitian. In matrix notation we have
Density matrix. It seems very strange and uncomfortable that our fundamental equation is a linear equation for the wavefunction while all observed quantities involve bi-linear combinations of the wavefunction.
More informationPhysics 7240: Advanced Statistical Mechanics Lecture 1: Introduction and Overview
Physics 7240: Advanced Statistical Mechanics Lecture 1: Introduction and Overview Leo Radzihovsky Department of Physics, University of Colorado, Boulder, CO 80309 (Dated: 30 May, 2017) Electronic address:
More informationEntanglement in Quantum Field Theory
Entanglement in Quantum Field Theory John Cardy University of Oxford Landau Institute, June 2008 in collaboration with P. Calabrese; O. Castro-Alvaredo and B. Doyon Outline entanglement entropy as a measure
More informationPhysics 342 Lecture 17. Midterm I Recap. Lecture 17. Physics 342 Quantum Mechanics I
Physics 342 Lecture 17 Midterm I Recap Lecture 17 Physics 342 Quantum Mechanics I Monday, March 1th, 28 17.1 Introduction In the context of the first midterm, there are a few points I d like to make about
More informationTensor network simulation of QED on infinite lattices: learning from (1 + 1)d, and prospects for (2 + 1)d
Tensor network simulation of QED on infinite lattices: learning from (1 + 1)d, and prospects for (2 + 1)d Román Orús University of Mainz (Germany) K. Zapp, RO, Phys. Rev. D 95, 114508 (2017) Goal of this
More information1 Mathematical preliminaries
1 Mathematical preliminaries The mathematical language of quantum mechanics is that of vector spaces and linear algebra. In this preliminary section, we will collect the various definitions and mathematical
More information6. Auxiliary field continuous time quantum Monte Carlo
6. Auxiliary field continuous time quantum Monte Carlo The purpose of the auxiliary field continuous time quantum Monte Carlo method 1 is to calculate the full Greens function of the Anderson impurity
More informationEntanglement in Quantum Field Theory
Entanglement in Quantum Field Theory John Cardy University of Oxford DAMTP, December 2013 Outline Quantum entanglement in general and its quantification Path integral approach Entanglement entropy in 1+1-dimensional
More informationA. Time dependence from separation of timescales
Lecture 4 Adiabatic Theorem So far we have considered time independent semiclassical problems. What makes these semiclassical is that the gradient term (which is multiplied by 2 ) was small. In other words,
More informationKitaev honeycomb lattice model: from A to B and beyond
Kitaev honeycomb lattice model: from A to B and beyond Jiri Vala Department of Mathematical Physics National University of Ireland at Maynooth Postdoc: PhD students: Collaborators: Graham Kells Ahmet Bolukbasi
More information