Entanglement in Spin Chains Vladimir Korepin Yang Institute for Theoretical Physics, State University of New York at Stony Brook

Total Page:16

File Type:pdf, Size:1020Kb

Entanglement in Spin Chains Vladimir Korepin Yang Institute for Theoretical Physics, State University of New York at Stony Brook Entanglement in Spin Chains Vladimir Korepin Yang Institute for Theoretical Physics, State University of New York at Stony Brook, Stony Brook -10 Entanglement is a resource for quantum computation. How much quantum effects we can use to control one quantum sys- tem by another. OUTLINE 1. Entanglement of two quantum subsystems A and B. Sim- plest case when the whole system fA; Bg is in a pure state: a unique wave function jΨA;Bi. Entropy of a subsystem is a measure of entanglement. 2. The wave function will be a unique ground state of a dynam- ical model: interacting spins, Bose gas or strongly correlated electrons. Universal properties of entropy of a large subsystem. -9 Easy to explain when entanglement is absent: jΨA;Bi = jΨAi ⊗ jΨBi () no entanglement But if the wave function is a sum of several such terms d jΨA;Bi = X jΨAi ⊗ jΨBi entangled. j=1 j i A B Here d > 1 and jΨj i are linear independent; as well as jΨj i. Measure of entanglement? C.H. Bennett, H.J. Bernstein, S.Popescu, and B.Schumacher 1996 entropy of a subsystem: 0 A;B A;B 1 S = −trA (ρA ln ρA) ; ρA = trB @jΨ ihΨ jA -8 Contra-intuitive features: Two spins 1=2 forming spin 1. j "i ⊗ j "i no entanglement but middle component with Sz = 0 is maximally entangled: j "i ⊗ j #i + j "i ⊗ j #i maximum entanglement Entanglement depends on Sz. Ferromagnetic XXX: L X x x y y z z Hf = − fσ σ + σj σj+1 + σ σ g j=1 j j+1 j j+1 0 1 PL z σ are Pauli matrices. [Hf;Sz] = 0, Sz = @ j=1 σj A We can fix Sz to fix ground state uniquely: 0 1 PL z Maximum of Sz = @ j=1 σj A j "i ⊗ j "i ::: ⊗ j "i no entanglement. -7 0 1 − PL − We can lower Sz by applying the operator S = @ j=1 σj A: − M jΨfmi = (S ) (j "i ⊗ j "i ::: ⊗ j "i) entangled Non-trivial magnetization Sz=L if M ! 1 proportionally to L ! 1. We shall artificially represent the ground state as by partite system: Block of spins on an space interval [1; x] is the subsystem A , the rest of the ground state is the subsystem B. Entanglement of a block of spins on a space interval [1; x] with the rest of the ground state jΨfmi. Density matrix of the block ρ(x) the entropy of the block S(x). 1 S(x) ! ln(x); x ! 1; F 2 Double scaling limit 1 << x << L. Salerno and Popkov 2004 -6 Anti - ferromagnetic XXX: X x x y y z z Haf = fσ σ + σj σj+1 + σ σ g j j j+1 j j+1 The ground state is unique and has spin 0. Also a gapless case, but for low lying excitations energy is proportional to momentum ∼ p. So one can use conformal field theory for evaluation of asymptotic of entropy of a large block of spins, belonging to the ground state: 1 S(x) ! ln(x); x ! 1; AF 3 Holzhey, Larsen and Wilczek 1994. Finite size corrections 1 X j S(x) = ln(x)=3 + x0 + aj=x j=1 -5 Zamolodchikov, Fatteev, Takhatajan and Babujian: H = X fX − X2g spin s = 1 1 n n n ~ ~ x x y y z z Xn = SnSn+1 = SnSn+1 + SnSn+1 + SnSn+1 Solvable by Bethe Ansatz. Higher spin s Faddeev 1983: H = X F (X );F (X)is a polynomial of degree 2s s n n 2s 2s 1 2s X − yj F (X) = 2 X X Y ; y = l(l+1)=2−s(s+1) j = 0 l l=0 k=l+1 k j 6= l yl − yj The entropy of a block of x spins V.K. 2003: s S (x) = ln x; as x ! 1 s + 1 -4 Other models: Bose gas with delta interaction: 2 3 Z y y y H = dx 4@ x@ x + g 5 : Here is a canonical Bose field and g > 0 is a coupling constant. @2 NX X HN = − + 2g δ(xk − xl) j=1 @xj N≥k>l≥1 Bethe Ansaz: E. Lieb, W.Liniger in 1963. The entropy of gas on a space interval [0; x] also scales as 1 S(x) ! ln(x); x ! 1 3 V.K, 2004 -3 The Hubbard model H: X y y X H = − (cj,σcj+1,σ + cj+1,σcj,σ) + u nj;"nj;# j=1 j=1 σ=";# y Here cj,σ is a canonical Fermi operator on the lattice [creates y of an electron] and nj,σ = cj,σcj,σ is an operator on number of electrons in cite number j with spin σ, u > 0. Lieb and Wu in 1968. Below half filling [less then one electron per lattice cite] both charge and spin degrees of freedom can be described by Virasoro algebra with central charge equal to 1, Korepin, Frahm 1990 : cc = 1 and cs = 1. 2 S (x) = ln x 3 V.K. 2004. At half filed band S (x) = (1=3) ln x -2 R´enyi entropy also scales logarithmically 0 −1 1 1 α B1+α C Sα = 1−α ln T r(ρA) = @ 6 A ln x Bai Qi Jin , Korepin 2004 -1 The entropy of a block of spins for a model with a gap ap- proaches a constant S1 as the size of the block increases: S(x) ! S1 as x ! 1 The Hamiltonian of XY model : 1 X x x y y z H = − (1 + γ)σ σ + (1 − γ)σ σn+1 + hσ n=−∞ n n+1 n n Here γ is anisotropy parameter: 0 < γ < 1 and 0 < h is a magnetic field and σn are Pauli matrices. Solution, phase transitions: Lieb, Schultz, Mattis, Barouch and McCoy. Toeplitz determinants and integrable Fredholm operators were used for evaluation of correlation functions, 0 The model has three different cases: s 1. Case 1a: 2 1 − γ2 < h < 2. medium magnetic field s 2. Case 1b: 0 ≤ h < 2 1 − γ2 , small magnetic field 3. Case 2: h > 2 , strong magnetic field A.Its, B.-Q. Jin, V. Korepin 2004: Case 1b 2 0 2 1 0 2 1 0 3 6 B C B C 7 1 6 B k C B k C 4I(k)I(k )7 6 B C B C 7 6 B C B C 7 S1 = 6 ln B C + B1 − C 7 + ln 2 6 4 @16k0 A @ 2 A π 5 p with k0 = 1 − k2, and v u u u 2 2 u u1 − (h=2) − γ Z π=2 dα u u s k = u I(k) = 0 t 1 − (h=2)2 1 − k2 sin2 α I(k) is the complete elliptic integral of the first kind: 1 I. Peschel 2004: Case 1a: medium magnetic field 2 0 2 1 0 2 1 0 3 6 B C B C 7 1 6 B k C B k C 4I(k)I(k )7 6 B C B C 7 6 B C B C 7 S1 = 6 ln B C + B1 − C 7 + ln 2; 6 4 @16k0 A @ 2 A π 5 and in Case 2 : strong magnetic field 2 0 3 6 7 1 6 16 4I(k)I(k )7 6 2 02 7 6 7 S1 = 6 ln + (k − k ) 7 ; (1) 12 4 (k2k02) π 5 8 > s > > 2 2 > > (h=2) + γ − 1 = γ; Phase 1a <> k = > (2) > s > > 2 2 > :> γ = (h=2) + γ − 1; Phase 2 Consider range of variation: S1 reaches minimum as ordered states and it reaches maximum at phase transitions [disordered states]. 2 Local minimum S1 = ln 2 at boundary between cases 1a and s 1b: h = 2 1 − γ2 The ground state is doubly degenerated: Y n jG1i = [cos(θ)j "ni + (−1) sin(θ)j #ni] n2lattice Y n jG2i = [cos(θ)j "ni − (−1) sin(θ)j #ni] n2lattice Each state is factorized and has no entropy. cos2(2θ) = (1 − γ)=(1 + γ) Absolute minimum of is reached at h ! 1 : S1 ! 0 as the state becomes ferromagnetic " ::: ". These are most ordered states. These is a gap. 3 Phase transitions: as the gap closes S1 ! +1 1. Critical magnetic field: h ! 2 and γ 6= 0: 1 1 2 S1 ! − ln j2 − hj + ln 4γ + O(j2 − hj ln j2 − hj) 6 3 2. An approach to XX model: γ ! 0 and 0 < h < 2: 1 1 2 1 2 S1 ! − ln γ + ln(4 − h ) + ln 2 + O(γ ln γ) 3 6 3 This agrees with conformal approach of P. Calabrese, J. Cardy and Toeplitz determinant approach of B.-Q. Jin, V.Korepin. 4 Let us consider another gapped spin chain introduced by Affleck, Kennedy, Lieb, and Tasaki: AKLT model also known as VBS. It consists of a chain of N spin-1's in the bulk, and two spin-1/2 ~ on the boundary. We shall denote by Sk the vector of spin-1 operators and by ~sb spin-1/2 operators at boundaries. 0 1 N−1 B 1 C X B 2C B ~ ~ ~ ~ C H = @SkSk+1 + (SkSk+1) A + π0;1 + πN;N+1: k=1 3 1 ~ ~ 1 ~ ~ 2 1 Notice that 2SkSk+1 + 6(SkSk+1) + 3 is a projector on a state of spin 2. The terms π describe interaction of boundary spin 1/2 and next spin 1; π is a projector on a state with spin 3/2: 2 0 1 2 0 1 ~ ~ π0;1 = @1 + ~s0S1A ; π = @1 + ~s S A : 3 N;N+1 3 N+1 N 5 Ground state is unique and there is a gap: # # # # # # v vv vv vv vv vv vv v 1"!"!"!"!"!"! a dot is spin-2; circle means symmetrisation [makes spin 1].
Recommended publications
  • (Oxon) on the THEORY of MODERN QUANTUM ALGORITHMS
    SKOLKOVO INSTITUTE OF SCIENCE AND TECHNOLOGY draft version Jacob Daniel Biamonte, BS, DPhil (Oxon) ON THE THEORY OF MODERN QUANTUM ALGORITHMS Specialization 01.01.03 Mathematical Physics Dissertation submitted for the degree of Doctor of Physical and Mathematical Sciences arXiv:2009.10088v1 [quant-ph] 21 Sep 2020 Moscow | 2020 2 Contents Page Chapter 1. Introduction: Programming Ground States .......5 1.1 Computation and the Ising model . .9 1.2 Low-energy subspace embedding . 15 1.3 Two-body reductions . 19 1.4 P- vs. NP problems and physics . 23 1.5 Computational phase transitions . 24 1.5.1 Thermal states at the SAT phase transition . 29 Chapter 2. Quantum vs Probabilistic Computation ......... 32 2.1 Defining Mechanics . 32 2.1.1 Stochastic time evolution . 33 2.1.2 Quantum time evolution . 33 2.1.3 Properties in stochastic mechanics . 34 2.1.4 Hamiltonian properties in quantum mechanics . 35 2.1.5 Observables . 36 2.2 Walks on Graphs: quantum vs stochastic . 44 2.2.1 Normalized Laplacians . 47 2.2.2 Stochastic walk . 50 2.2.3 Quantum walk . 51 2.2.4 Perron's theorem . 54 2.3 Google page rank|a ground eigenvector problem . 59 2.4 Kitaev's quantum phase estimation algorithm . 61 2.5 Finding the ground state . 66 Chapter 3. Tensor Networks and Quantum Circuits ......... 69 3.1 Clifford gates . 69 3.2 Tensor network building blocks . 72 3.2.1 Reversible logic . 75 3.2.2 Heisenberg picture . 77 3.2.3 Stabilizer tensor theory . 79 3 Chapter 4. Variational Search and Optimization ..........
    [Show full text]
  • Algebraic Bethe Ansatz and Tensor Networks
    Algebraic Bethe Ansatz and Tensor Networks A Dissertation presented by You Quan Chong to The Graduate School in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy in Physics Stony Brook University May 2015 Stony Brook University The Graduate School You Quan Chong We, the dissertation committee for the above candidate for the Doctor of Philosophy degree, hereby recommend acceptance of this dissertation Vladimir E. Korepin - Dissertation Advisor Professor, C. N. Yang Institute for Theoretical Physics, Department of Physics and Astronomy Tzu-Chieh Wei - Chairperson of Defense Assistant Professor, C. N. Yang Institute for Theoretical Physics, Department of Physics and Astronomy Dominik A. Schneble - Committee Member Associate Professor, Department of Physics and Astronomy Scott Sutherland - Outside Member Associate Professor, Department of Mathematics This dissertation is accepted by the Graduate School Charles Taber Dean of the Graduate School ii Abstract of the Dissertation Algebraic Bethe Ansatz and Tensor Networks by You Quan Chong Doctor of Philosophy in Physics Stony Brook University 2015 We consider several exactly solvable models of strongly correlated elec- trons in one dimension, such as the Heisenberg XXX model, the supersym- metric t-J model and the Hubbard model. These models can be solved by using the method of graded algebraic Bethe ansatz. We use it to design graded tensor networks which can be contracted approximately to obtain a Matrix Product State. This overcomes a major shortcoming of current density matrix renormal- ization group (DMRG) methods which work well on the ground states, but have difficulty working with the excited states of such models. In addition, observables such as correlation functions are important as they are experimentally measurable, but have been analytically described in the double scaling limit only.
    [Show full text]
  • Entanglement and Quantum Spin Chains
    ABitofHistoryofQuantumEntanglement Quantum Information and Dynamical Systems ENTANGLEMENT AND QUANTUM SPIN CHAINS V. Korepin YITP of Stony Brook University New exactly solvable spin chain came out of quantum information V. Ko r e p in F-H Formula ABitofHistoryofQuantumEntanglement Quantum Information and Dynamical Systems Outline 1 ABitofHistoryofQuantumEntanglement 2 Quantum Information and Dynamical Systems V. Ko r e p in F-H Formula ABitofHistoryofQuantumEntanglement Quantum Information and Dynamical Systems The dawn Aug 14, 1935 Shroedinger submitted his paper to Proceedings of Cambridge Philosophical Society. 1964 John Bell: entanglement distinguishes quantum mechanics from classical at fixed value of Plank constant. 1981 Richard Feynman: quantum computation. 1996 Bennett, Bernstein, Popescu and Schumacher introduced entropy of subsystem as a measure of entanglement discovery of entanglement entropy V. Ko r e p in F-H Formula ABitofHistoryofQuantumEntanglement Quantum Information and Dynamical Systems Application to physics: area Law 1986 L. Bombelli, R. Koul, J. Lee, R. Sorkin; 1993 M. Srednicki 1994 Holzhey, Larsen and Wilczek: ent entr of block of spins of c size x for 1D gapless models scales as S = 3 log x . 2003 Logarithmic scaling of Renyi entropy. B.-Q.Jin, Korepin, Journal Statistical Physics; arXiv.org April 15 of 2003. ln tr(ρα) −1 S (α)= 1+α ln x R 1−α → 6 c v πTx 2004 For finite temperature S = 3 log πT sinh v Korepin, PRL, vol 92, ei 096402, 05 March 2004,! ! "" arXiv:cond-mat/0311056 V. Ko r e p in F-H Formula ABitofHistoryofQuantumEntanglement Quantum Information and Dynamical Systems Environment and Error Correcting Codes Q-bits should be entangled with other qubits in quantum computer, but not entangled with the environment.
    [Show full text]
  • Entanglement Entropy & Spectrum • Holographic Spin Chain (VBS/CFT Correspondence) 3
    Entanglement spectra of two- dimensional solvable models Hosho Katsura (Gakushuin University) Collarborators: Anatol Kirillov (RIMS, Kyoto Univ.) Vladimir Korepin (YITP, Stony Brook) Naoki Kawashima (ISSP, Tokyo Univ.) Lou Jie (ISSP Fudan Univ.) Shu Tanaka (Tokyo Univ.) Ryo Tamura (NIMS) H.K., et al., J. Phys. A 43, 255303 (2010). J. Lou, S. Tanaka, H.K. & N. Kawashima, PRB 84, 245128 (2011). S. Tanaka, R. Tamura, & H.K., PRA 86, 032326 (2012). Outline 1. Introduction • AKLT model & valence-bond-solid (VBS) state • Schmidt decomposition & quantum entanglement 2. Entanglement spectra of 2d AKLT models • Gram matrix & Reflection symmetry • Entanglement entropy & spectrum • Holographic spin chain (VBS/CFT correspondence) 3. Entanglement spectra of quantum hard-square models • Tensor network states for interacting Rydberg atom systems • Entanglement spectra • Holographic minimal models (c<1) A brief history of AKLT and valence-bond-solid (VBS) state • Haldane gap problem Haldane’s conjecture (F. D. M. Haldane, Phys. Rev. Lett. 50 (‘83).) S=integer antiferromagnetic (AFM) Heisenberg chains are gapped. Experiment in S=1 spin chain: Ni(C2H8N2)2NO2(ClO4) (NENP) etc. • Affleck-Kennedy-Lieb-Tasaki (AKLT) model (PRL 59 (‘87), CMP 115 (‘87).) 1. Exact unique ground state valence-bond-solid (VBS) state 2. Rigorous proof of the ‘Haldane’ gap in this model 3. AFM correlation decays exponentially with distance. • Valence-bond-solid (VBS) state : spin singlet : symmetrization VBS state Scwinger boson(SB) rep. of spin operator Arovas, Auerbach a b & Haldane, PRL 60 (‘88). Constraint : S=1 state is spanned by . Matrix product state (MPS) Fannes et al.,(‘89), Kluemper et al.,(‘91). VBS on arbitrary graphs Square lattice Hexagonal lattice • Projected pair entangled state (PEPS), Tensor network state (TNS) Verstraete & Cirac (PRA 70 (R) (‘04)), Gross & Eisert (PRL 98 (‘07)).
    [Show full text]
  • SCIENTIFIC REPORT for the YEAR 2002 ESI, Boltzmanngasse 9, A-1090 Wien, Austria
    The Erwin SchrÄodinger International Boltzmanngasse 9 ESI Institute for Mathematical Physics A-1090 Wien, Austria Scienti¯c Report for the Year 2002 Vienna, ESI-Report 2002 March 1, 2003 Supported by Federal Ministry of Education, Science, and Culture, Austria 2 ESI{Report 2002 ERWIN SCHRODINGERÄ INTERNATIONAL INSTITUTE OF MATHEMATICAL PHYSICS, SCIENTIFIC REPORT FOR THE YEAR 2002 ESI, Boltzmanngasse 9, A-1090 Wien, Austria, March 1, 2003 Honorary President: Walter Thirring, Tel. +43-1-4277-51516. President: Jakob Yngvason: +43-1-4277-51506. [email protected] Director: Peter W. Michor: +43-1-3172047-16. [email protected] Director: Klaus Schmidt: +43-1-3172047-14. [email protected] Administration: Maria Windhager, Eva Kissler, Ursula Sagmeister: +43-1-3172047-12, [email protected] Computer group: Andreas Cap, Gerald Teschl, Hermann Schichl. International Scienti¯c Advisory board: Jean-Pierre Bourguignon (IHES), Luis A. Ca®arelli (Austin, Texas), Giovanni Gallavotti (Roma), Krzysztof Gawedzki (IHES), Viktor Kac (MIT), Elliott Lieb (Princeton), Harald Grosse (Vienna), Harald Niederreiter (Vienna), Table of contents General remarks . 2 Report of the Review Panel . 3 Winter School in Geometry and Physics . 12 Arithmetic Groups and Automorphic Forms . 12 Stability Matters: A Symposium on Mathematical Physics . 12 PROGRAMS IN 2002 . 13 Developed Turbulence . 13 Arithmetic, automata, and asymptotics . 18 Quantum ¯eld theory on curved space time . 20 Aspects of foliation theory in geometry, topology and physics . 22 Noncommutative geometry and quantum ¯eld theory Feynman diagrams in mathematics and physics . 24 Mathematical population genetics and statistical physics . 25 CONTINUATION OF PROGRAMS FROM 2001 and earlier . 26 SENIOR FELLOWS and GUESTS via Director's shares .
    [Show full text]
  • Arxiv:Cond-Mat/0405155V3
    Universality of the one dimensional Bose gas with delta interaction Luigi Amico(1) and Vladimir Korepin(2) (1) MATIS-INFM & Dipartimento di Metodologie Fisiche e Chimiche (DMFCI), Universit´adi Catania, viale A. Doria 6, I-95125 Catania, Italy and (2) C.N. Yang Institute for Theoretical Physics, State University of New York at Stony Brook, NY 11794-3840 We consider several models of interacting bosons in a one dimensional lattice. Some of them are not integrable like the Bose-Hubbard others are integrable. At low density all of these models can be described by the Bose gas with delta interaction. The lattice corrections corresponding to the different models are contrasted. I. INTRODUCTION duces quantum fluctuations that dominate the zero tem- perature phase diagram. At commensurate number of bosons per site the Bose-Hubbard model predicts a Mott The Bose gas with delta interaction is one of the insulator-superfluid quantum phase transition controlled most famous models for strongly interacting bosons [1, 2] by the ratio between the bosonic tunneling rate t and (the corresponding classical dynamics is described by the Coulomb repulsion. In optical applications the sys- the non-linear Schr¨odinger equation). In one dimension tem parameters are tuned by changing the laser intensity this model possesses infinitely many, mutually commut- [20]; the criticality of the superfluid-insulator quantum ing integrals of the motion and is therefore called inte- phase transition is substantially altered by the confining grable. The integrability of the Bose gas can be for- harmonic potential [21]. mulated within the Quantum Inverse Scattering Method (QISM) [3] that exploit a crucial property of the two- In this paper we focus on one dimension.
    [Show full text]
  • Arxiv:1806.04049V3 [Hep-Th] 20 Sep 2018 So Far
    R´enyi entropy of highly entangled spin chains Fumihiko Sugino∗ and Vladimir Korepiny ∗ Fields, Gravity & Strings Group, Center for Theoretical Physics of the Universe, Institute for Basic Science (IBS), 55, Expo-ro, Yuseong-gu, Daejeon 34126, Republic of Korea y C.N.Yang Institute for Theoretical Physics, Stony Brook University, NY 11794, USA [email protected], [email protected] Abstract Entanglement is one of the most intriguing features of quantum theory and a main resource in quantum information science. Ground states of quantum many-body systems with local interactions typically obey an \area law" meaning the entanglement entropy proportional to the boundary length. It is exceptional when the system is gapless, and the area law had been believed to be violated by at most a logarithm for over two decades. Recent discovery of Motzkin and Fredkin spin chain models is striking, since these models provide significant violation of the entanglement beyond the belief, growing as a square root of the volume in spite of local interactions. Although importance of intensive study of the models is undoubted to reveal novel features of quantum entanglement, it is still far from their complete understanding. In this article, we first analytically compute the R´enyi entropy of the Motzkin and Fredkin models by careful treatment of asymptotic analysis. The R´enyi entropy is an important quantity, since the whole spectrum of an entangled subsystem is reconstructed once the R´enyi entropy is known as a function of its parameter. We find non-analytic behavior of the R´enyi entropy with respect to the parameter, which is a novel phase transition never seen in any other spin chain studied arXiv:1806.04049v3 [hep-th] 20 Sep 2018 so far.
    [Show full text]
  • Measures of Entanglement in Dynamical Systems Density Matrix
    Measures of Entanglement in Dynamical Systems Density matrix of block of spins in the ground state of a Hamiltonian Vladimir Korepin Yang Institute for Theoretical Physics State University of New York at Stony Brook Abstract To build quantum computer or any other quan- tum device we need to maintain quantum con- trol. Entanglement is an important recourse for quantum control. Many-body physics is used to study measures of entanglement in ground state of different dynamical systems. 1 Binary system A&B , which is in a pure state: a { } unique wave function denoted by ΨA,B . it will be a | i unique ground state of a dynamical model: interacting spins, Bose gas or Hubbard model. d ΨA,B = ΨA ΨB , d> 1 | i | j i⊗| j i jX=1 Measure of entanglement. Von Neumann entropy of the subsystem S(ρ ) = T r (ρ ln ρ ) A − A A A ρ = T r ΨA,B ΨA,B A B | ih | John Preskill, Matthew Hastings, Jens Eisert, Frank Verstr aete, John Cardy 2 Another measure is R´enyi entropy of a subsystem: 1 α Sα(ρ ) = ln T r (ρ ) A 1 α A A where α is a parameter. − spectrum of density matrix is also in- teresting, even eigenvectors are impor- tant [measurement]. If whole binary system A&B is in a mixed states [thermodynamics]. Mutual entropy: I = S + S S {A,B} A B − A∪B Another measure of entanglement of mixed states was suggested by Asher Peres. We can describe the mixed binary system by a density matrix ρA,B.
    [Show full text]
  • Slow Mixing of Glauber Dynamics for the Six-Vertex Model in the Ordered Phases
    Slow Mixing of Glauber Dynamics for the Six-Vertex Model in the Ordered Phases Matthew Fahrbach∗ and Dana Randall† School of Computer Science, Georgia Institute of Technology December 23, 2020 Abstract The six-vertex model in statistical physics is a weighted generalization of the ice model on (i.e., Eulerian orientations) and the zero-temperature three-state Potts model (i.e., proper three-colorings). The phase diagram of the model depicts its physical properties and suggests where local Markov chains2 will be ecient. In this paper, we analyze the mixing time of Glauber dynamics for the six-vertexℤ model in the ordered phases. Specically, we show that for all Boltzmann weights in the ferroelectric phase, there exist boundary conditions such that local Markov chains require exponential time to converge to equilibrium. This is the rst rigorous result bounding the mixing time of Glauber dynamics in the ferroelectric phase. Our analysis demonstrates a fundamental connection between correlated random walks and the dynamics of intersecting lattice path models (or routings). We analyze the Glauber dynamics for the six-vertex model with free boundary conditions in the antiferroelectric phase and signicantly extend the region for which local Markov chains are known to be slow mixing. This result relies on a Peierls argument and novel properties of weighted non-backtracking walks. arXiv:1904.01495v2 [cs.DS] 22 Dec 2020 ∗Email: [email protected]. Supported in part by an NSF Graduate Research Fellowship under grant DGE-1650044. †Email: [email protected]. Supported in part by NSF grants CCF-1526900, CCF-1637031, and CCF-1733812.
    [Show full text]
  • Arxiv:1804.10133V3 [Cond-Mat.Quant-Gas] 17 Dec 2018
    Rise and fall of a bright soliton in an optical lattice Piero Naldesi1 ∗, Juan Polo Gomez1, Boris Malomed2;3, Maxim Olshanii4, Anna Minguzzi1 and Luigi Amico5;6;7;8;9 1 Universit´eGrenoble-Alpes, LPMMC, F-38000 Grenoble, France and CNRS, LPMMC, F-38000 Grenoble, France 2 Department of Physical Electronics, School of Electrical Engineering, Faculty of Engineering, Tel Aviv University, P.O.B. 39040, Ramat Aviv, Tel Aviv, Israel 3 Center for Light-Matter Interaction, Tel Aviv University, P.O.B. 39040, Ramat Aviv, Tel Aviv, Israel 4 Department of Physics, University of Massachusetts Boston, Boston, MA 02125, USA 5 Dipartimento di Fisica e Astronomia, Via S. Sofia 64, 95127 Catania, Italy 6 Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 117543, Singapore 7 MajuLab, CNRS-UNS-NUS-NTU International Joint Research Unit, UMI 3654, Singapore 8 CNR-MATIS-IMM & INFN-Sezione di Catania, Via S. Sofia 64, 95127 Catania, Italy and 9 LANEF 'Chaire d'excellence', Universit`eGrenoble-Alpes & CNRS, F-38000 Grenoble, France (Dated: December 19, 2018) We study an ultracold atomic gas with attractive interactions in a one-dimensional optical lattice. We find that its excitation spectrum displays a quantum soliton band, corresponding to N-particle bound states, and a continuum band of other, mostly extended, states. For a system of a finite size, the two branches are degenerate in energy for weak interactions, while a gap opens above a threshold value of the interaction strength. We find that the interplay between degenerate extended and bound states has important consequences for both static and dynamical properties of the system.
    [Show full text]
  • The Hubbard Dimer: a Density Functional Case Study of a Many-Body Problem
    The Hubbard Dimer: A density functional case study of a many-body problem D. J. Carrascal1;2, J. Ferrer1;2, J. C. Smith3 and K. Burke3 1 Department of Physics, Universidad de Oviedo, 33007 Oviedo, Spain 2 Nanomaterials and Nanotechnology Research Center, Oviedo, Spain and 3 Departments of Chemistry and of Physics, University of California, Irvine, CA 92697, USA (Dated: September 3, 2015) This review explains the relationship between density functional theory and strongly correlated models using the simplest possible example, the two-site Hubbard model. The relationship to traditional quantum chemistry is included. Even in this elementary example, where the exact ground-state energy and site occupations can be found analytically, there is much to be explained in terms of the underlying logic and aims of Density Functional Theory. Although the usual solution is analytic, the density functional is given only implicitly. We overcome this difficulty using the Levy-Lieb construction to create a parametrization of the exact function with negligible errors. The symmetric case is most commonly studied, but we find a rich variation in behavior by including asymmetry, as strong correlation physics vies with charge-transfer effects. We explore the behavior of the gap and the many-body Green's function, demonstrating the `failure' of the Kohn-Sham method to reproduce the fundamental gap. We perform benchmark calculations of the occupation and components of the KS potentials, the correlation kinetic energies, and the adiabatic connection. We test several approximate functionals (restricted and unrestricted Hartree-Fock and Bethe Ansatz Local Density Approximation) to show their successes and limitations.
    [Show full text]
  • Deforming the Fredkin Spin Chain Away from Its Frustration-Free Point
    Deforming the Fredkin spin chain away from its frustration-free point Khagendra Adhikari∗ and K. S. D. Beachy Department of Physics and Astronomy, The University of Mississippi, University, Mississippi 38677, USA The Fredkin model describes a spin-half chain segment subject to three-body, correlated-exchange interactions and twisted boundary conditions. The model is frustration-free, and its ground state wave function is known exactly. Its low-energy physics is that of a strong xy ferromagnet with gapless excitations and an unusually large dynamical exponent. We study a generalized spin chain model that includes the Fredkin model as a special tuning point and otherwise interpolates between the conventional ferromagnetic and antiferromagnetic quantum Heisenberg models. We solve for the low-lying states, using exact diagonalization and density-matrix renormalization group calculations, in order to track the properties of the system as it is tuned away from the Fredkin point; we also present exact analytical results that hold right at the Fredkin point. We identify a zero-temperature phase diagram with multiple transitions and unexpected ordered phases. The Fredkin ground state turns out to be particularly brittle, unstable to even infinitesimal antiferromagnetic frustration. We remark on the existence of an “anti-Fredkin” point at which all the contributing spin configurations have a spin structure exactly opposite to those in the Fredkin ground state. I. INTRODUCTION z Stot = 0 The construction of a spin-1 quantum spin chain with an dimerization pattern z AFM exactly solvable ground state [1] and its later generalization to I = spin-S (for all integer S > 1) by Movassagh and Shor [2] have proved to be incredibly fruitful developments.
    [Show full text]