Complexity of thermal states in quantum spin chains arxiv:85.449v [quant-ph] Aug 8 Marko Žnidarič, Tomaž Prosen and Iztok Pižorn Department of physics, FMF, University of Ljubljana, Jadranska 9, SI- Ljubljana, Slovenia (Dated: August, 8) We study quantum correlations and complexity of simulation, characterized by quantum mutual information and entanglement entropy in operator space respectively, for thermal states in critical, non-critical and quantum chaotic spin chains. A simple general relation between the two quantities is proposed. We show that in all cases mutual information and entanglement entropy saturate with the system size, whereas as a function of the inverse temperature, we find logarithmic divergences for critical cases and uniform bounds in non-critical cases. A simple efficient quasi-exact method for computation of arbitrary entropy related quantities in thermalized XY spin chains is proposed. PACS numbers:.3.ik, 5.7.Fh, 75..Pq, 3.67.Mn I. INTRODUCTION Physical properties of non-integrable strongly correlated quantum systems are notoriously difficult to calculate due to exponential growth of Hilbert space dimension with the system size. Nevertheless, many numerical and analytical approaches have been developed that provide results for special properties and specific systems. The question arises, which quantities are amenable to classical calculation in spite of generic system complexity? Recently, two results emerging from the field of quantum information shed light on the question. First, it has been shown that for one-dimensional (D) systems the entanglement content of ground states grows logarithmically with the number of particles for fermionic critical systems and saturates for non-critical or bosonic systems [,, 3, 4]. Second, density matrix renormalization group (DMRG) method [5] has been reinterpreted and optimized within matrix product state (MPS) ansatz [6, 7] for the many body wave-function. Its complexity is essentially given by the entanglement entropy of the system s bipartition. Despite being essentially the best method available it is not known precisely in which instances is DMRG computationally efficient, for example it has been shown to be inefficient for computing the evolution of homogeneous D non-integrable systems in real time [8] but becomes efficient in the presence of disorder due to the emergence of (many body) localization [9]. Thermal states at finite temperature have been simulated using an extension of MPS to density operators (matrix product operators MPO), and by performing the evolution in imaginary time (inverse temperature) [7]. However, no analysis of efficiency scaling with the system size n and inverse temperature have been performed. Note that finding a ground state corresponding to the limit is proven to be a hard problem [], even though it can often be performed efficiently [5]. For a different algorithm for computing thermal states see Ref. []. In Ref. [] an upper bound on the necessary MPO dimension is found for thermal states which unfortunately scales exponentially with. It is not clear if one can do better, that is, if finite temperature states can be calculated efficiently, uniformly in? In the present work we show that this is indeed possible. Complexity will be characterized by the operator space entanglement entropy (OSEE) [3] of thermal states which is directly linked to the minimal necessary rank of MPO ansatz. OSEE is defined as von Neumann entropy of the reduced density matrix defined with respect to a physical density operator treated as an element of a Hilbert (Fock) space of operators. In addition we shall show that OSEE behaves in the same way as the quantum mutual information (QMI) I A:B [4], which is a natural measure of quantum correlations in a bipartition A + B and an upper bound for distillable entanglement. Our main result is that OSEE and I A:B in thermodynamic limit n grow as (c/3)log for critical (gapless) systems and saturate for non-critical (gapped) systems, while they always saturate with for a finite n. This shows that the linear in upper bound on QMI of Ref.[5] is far from optimal and DMRG computations of large (infinite) critical D quantum systems should be efficient at all nonzero temperatures. This has to be contrasted with typical linear growth of OSEE in t in real time evolution [8]. Note that OSEE and QMI always saturate with n, for any fixed, and hence scale very differently from thermodynamic entropy of a block of length n, which for a critical system behaves as [6] S n for large n. Our results are based on analytical and numerical calculations in several examples of integrable and non-integrable spin chains. We outline a novel method for quasi-exact calculations of OSEE and I A:B in XY spin / chains using the concept of Fock space of operators [7], whereas we rely on numerical MPO-DMRG simulations for the other cases. II. METHODS First, we describe how to construct a thermal state and compute OSEE in a general spin chain, regardless of its integrability, using MPO ansatz and imaginary time evolution. A thermal state ρ() = exp( H)/ tr[exp( H)] of Hamiltonian H of a chain of n spins / can be expressed as an element of 4 n dimensional Hilbert space K, ρ = s c s σ s spanned
by Pauli operators σ s σ s σsn n, s s,...,s n, s m {,,, 3}, σ ½, σ,,3 σ x,y,z. Inner product in K is defined as x y = n trx y. Within MPO ansatz expansion coefficients c s are represented in terms of 4n D D matrices A s m, m =,...,n, s =,,, 3, as c s = tr(a s Asn n ). () For exact MPO representation matrix dimension of A s m must be equal to the number of nonzero Schmidt coefficients for a bipartite cut of a vector c s at the corresponding site. Although for generic density operators ρ the required dimension D is exponentially large (in n) in practice one truncates the matrices to a smaller dimension D and uses () as a variational ansatz. As a rough estimate of a minimal D S, which is a direct measure of complexity of ρ, one can use OSEE S defined as an entanglement entropy of a super-ket ρ S = tr A (R log R), R = ρ ρ tr B ρ ρ () where a subscript A denotes the first n A sites and B its complement. More detailed analysis [8] specifies exact criteria for the applicability of matrix product ansatz () in terms of properties of the moments of Schmidt coefficients. Writing the expansion coefficients c s as a na n na dimensional matrix C sa,s B = c sa s B, the reduced density matrix () is given by R = CC T / tr(cc T ) and quantifies the bipartite correlations and the difficulty of a MPO simulation of a given density operator ρ. For use of OSEE in a non-mpo context see Ref. [9]. S is directly related to the entanglement entropy of physical states only in the case of pure states, ρ = ψ ψ, in which case S is just twice the entanglement entropy of a pure state ψ, e.g. the ground state at zero temperature. However, even though S is not an entanglement monotone for general ρ it can be used as a bipartite entanglement estimator for thermal states ρ(). Namely, as we will demonstrate later, its scaling on n and is the same as that of QMI I A:B, defined as I A:B = S(ρ A ) + S(ρ B ) S(ρ), (3) where ρ A = tr B ρ, ρ B = tr A ρ, and S(σ) = tr(σ log σ). QMI measures total correlations between A and B and is an upper bound for relative entropy of entanglement which in turn is an upper bound for distillable entanglement E D (ρ) []. MPO representation for a thermal state ρ(), apart from normalization, is computed by performing imaginary time evolution in, ρ( + ǫ) = exp( ǫh/)ρ() exp( ǫh/), starting at an infinite temperature state ρ() = ½/ n having D =. Such -evolution increases OSEE and as a consequence, dimension D has to increase as well in order to preserve the accuracy. Assuming the empirical fact that a fixed small ǫ (in our case ǫ =.5) is enough to produce accurate results in a finite number of steps (for all < ), we conclude that time complexity of simulation is polynomial in D and. Note that in contrast to a (unitary) real time evolution, imaginary time evolution does not preserve canonical Schmidt decomposition structure of matrices A sm m [6]. In order to improve stability we apply local rotations every few steps in order to reorthogonalize the matrices A sm m. Quasi-exact calculation of S and I A:B for thermal states is possible for integrable chains which can be mapped to quadratic fermionic models by Wigner- Jordan transformation, e.g. D XY model. Most general such situation is described by a Hamiltonian H = w Hw n j,l= H jlw j w l where H is an anti-symmetric Hermitian n n matrix and w m = σm x m l= σz l, w m = σm y m l= σz l are Majorana operators satisfying {w j, w l } w j w l +w l w j = δ jl. It is crucial to recognize [3, 7] that K is in fact a Fock space if spanned by 4 n products of w j, P α = w α wαn, α j {, }, and equipped with adjoint Fermi maps ĉ j P α = δ αj, w j P α. In our notation ˆ designates a linear map over the Fock space of operators K. Note that OSEE is just an entanglement entropy of an element of K. Therefore, using the formalism of [] lifted to an operator space, one can calculate S from the anti-symmetric correlation matrix Γ defined by e H â p â q e H e H e H = δ pq + iγ pq, p, q =,...,4n (4) where â j (ĉ j + ĉ j ), â j i(ĉ j ĉ j ), j =,...,n satisfying {â p, â q } = δ pq. Thermal state in (4) can be written in terms of dynamics e H = exp( Ĥ) ½ (5) using a quadratic adjoint Hamiltonian map Ĥ over K, Ĥ = (ĉ + ĉ ) H(ĉ + ĉ ), since (ĉ j + ĉ j ) ρ = w jρ. Diagonalizing the matrix H we obtain eigenvalues λ k, λ k, ordered such that λ... λ n, and the corresponding orthonormal eigenvectors v k, v k in terms of which we write the adjoint Hamiltonian in (5) in the normal form Ĥ = n k= 4λ k(ˆb kˆb k ) where ˆb k = v k (ĉ + ĉ )/. Identities ½ ˆb k ½ =, ½ ˆb kˆbl ½ =, ½ ˆb kˆb l ½ = δ kl allow us to employ the Wick theorem to evaluate (4) Γ j,l = Γ j,l = Γ j,l = δ jl tanh(4λ k )Im(vk,jv k,l ) k= ( ) Re(v cosh(4λ k ) k,j v k,l ). (6) k= Finally, OSEE is obtained from the eigenvalues ±iν j of 4n A 4n A upper-left block of Γ as S = n A j= H ( +νj ) where H (x) = xlog x ( x)log ( x) []. Similarly, we can also easily compute the block entropy S(ρ A ) of a thermal state ρ() by diagonalizing H = n k= 4λ k(b k b k ) with operators b k = v k w/ (which are elements of K and not maps over K like ˆbk ) and computing the correlation matrix Γ, defined as
3 tr (w j w l e H )/ tr e H = δ jl + iγ jl, j, l =,...,n. Expressing w j and H in terms of b k, b k, and using the Wick theorem with tr (b k b l) = δ kl tr ½ results in a general -point thermal correlation function Γ jl = tanh(λ k )Im(vk,jv k,l ). (7) k= Diagonalizing n A n A block of Γ corresponding to the sub-lattice A, yielding eigenvalues ±iν j, we obtain the block entropy S(ρ A ) = n A j= H ( +ν j ). Considering also B, and A + B instead of A we get QMI (3). III. RESULTS We shall show results for OSEE and QMI for several families of critical and non-critical as well as integrable and non-integrable quantum spin chains n H XY = l= n H I = σl x σl+ x + l= ( +γ σx l σ x l+ + γ σy l σy l+ ) + hσl z, l= (h x σl x + h z σl z ), (8) l= n n H XXZ = (σl x σx l+ + σy l σy l+ + σz l σz l+ ) + h l σl z, l= all for open boundary conditions, and a symmetric cut at n A = n/ where largest bipartite entanglement is expected. Majorana representation of the integrable XY case which is necessary for explicit calculations [Eqs. (4-7)] is H XY = i n l= (+γ w lw l+ γ w l w l+ ) i n l= hw l w l. In Fig. we show the results for quantum critical systems for which the entanglement entropy of the ground state grows as S (c/3)log n [] where c is the central charge of the corresponding conformal field theory [3]. However, we take several very large n and plot dependence of OSEE and QMI, for (a) XX model H XY (γ =, h = ), (b) transverse Ising model H XY (γ =, h = ), and (c) XXZ model H XXZ ( =.5, h l ). These critical models are all completely integrable, however we were able to make exact calculations only for (a,b), whereas we used MPO-DMRG simulations for (c) where - since evaluating QMI (3) is computationally difficult even within MPO ansatz - we have instead computed mutual purity, defined as I P = log [P(ρ)/(P(ρ A )P(ρ B ))] with P(σ) = tr σ. It is well known that purity often well describes the behavior of von Neumann entropy and its complement is sometimes termed linear entropy. In addition we checked for thermal states of small systems that mutual purity and QMI behave similarly. From Fig. we conclude that in the large n limit OSEE grows logarithmically with l= S () = c 3 log + c, (9) with parameters c = and c =.7 for XX model (a), c = / and c =.75 for transverse Ising model (b), and c = and c =.63 for XXZ model (c). c is exactly the central charge of the corresponding model. Presently we do not have an explanation for this finding and we state it as a conjecture. For fixed n, S saturates for larger than the inverse spectral gap = / E at twice the ground state value S (c/3)log n., I A:B, I A:B, log (P/P A P B ).5.5.5 4 3.5 3.5.5.5 3.5.5.5 (a) (b) I A:B I A:B. (c). 8 3 4 n=56 4 n=56 n=56 log (P/P A P B ) FIG. : Logarithmic growth of OSEE S (upper curves) and QMI I A:B (lower curves) for thermal states of critical models: (a) XX chain H XY (γ =, h = ), (b) transverse Ising chain H XY (γ =, h = ), both for n = 4, 56, and (c) XXZ chain H XXZ( =.5, h l ) for n = 56, 3, 8. In case (c) we show mutual purity I P instead of I A:B. Dotted lines indicate conjectures (9), and () combined with (9) [see text]. Next, let us have a look at non-critical systems. First, we will chose two non-integrable spin models that display typical signatures of quantum chaos (Fig. ): (a) Heisenberg XXZ model in a staggered magnetic field H XXZ ( =.5, h l = ( + ( ) l )/), and (b) Ising model in a tilted magnetic field H I (h x =, h z = ). We can see that for non-critical systems which possess a finite energy gap E between the ground and first excited states, OSEE saturates at n and independent values for sufficiently large, so the classical simulation of
4, log (P/P A P B ), log (P/P A P B )..8.6.4..7.6.5.4.3. n=3,8 n=8 n=3,8 n=8 log (P/P A P B ) log (P/P A P B ). (b) 5 5 FIG. : Saturation of OSEE S (upper curves) and mutual purity I P (lower curves) for quantum chaotic noncritical systems, namely (a) for XXZ model in a staggered field H XXZ( =.5, h l = ( + ( ) l )/) and (b) for tilted Ising H I(h x =, h z = ), for n = 8,3, 8. Note that the curves for n = 3 and n = 8 are practically indistinguishable., I.8.6.4..8 n= 4 6 I A:B I A:B I A:B.6.4...4.6.8 h 4 6 8 ξ FIG. 3: OSEE S (full curves) and QMI I A:B (broken curves) for noncritical XY chain with γ =.5, h =.9, and n =, 4, 6. Note a crossover between two saturation plateaus due to a single exponentially close (in n) excited state. The inset shows h dependence of the gap decay length ξ, computed from the fit to E e n/ξ, for n = 3,..., 6. Vertical line indicates h = p γ. quantum thermal states for gapped systems is even more efficient than for critical ones. Although this might seem expected, one should know that for OSEE we do not have theoretical uniform bounds in n like for QMI [5]. Note that in the presence of weakly broken symmetries, the gap E between two lowest states can get exponentially small with n even though the system is not critical. This (a) for instance happens in XY model. In Fig. 3 we show results for a non-critical integrable XY model (8) with γ =.5 and h =.9, where = / E grows exponentially with n. Actually we have shown for any H XY, that E = 4λ where λ is the smallest eigenvalue of matrix H (5), and which in turn has been found to decay exponentially λ exp( n/ξ). Interestingly, for h + γ <, the gap decay length ξ is practically insensitive to parameters γ, h, whereas at the point h = γ, where the ground state is doubly degenerate and separable [], ξ instantly starts increasing with increasing h h (inset of Fig. 3). However, there is no phase transition at h. In all cases studied QMI behaves in essentially the same way as OSEE. In fact, for models solvable by Wigner-Jordan transformation I A:B () almost overlaps with S evaluated at /4 I A:B () S (/4) () apart from finite size effects at small and large (Figs.,3). We conjecture that () becomes asymptotically exact for large n. Factor 4 may be understood from comparing the formalisms of Majorana operators w j and the adjoint maps ĉ j. Since QMI is an upper bound for distillable entanglement [] this means that exact quantum bipartite entanglement for thermal states is also small, i.e., it is either upper bounded by n and -independent value for non-critical systems (Figs.,3) or it is upper bounded by n-independent value log for critical systems. IV. CONCLUSIONS By studying several typical models, we have demonstrated that classical simulations of thermal states of interacting quantum systems in D are efficient. Using DMRG simulations combined with exact calculations for XY models we found that in non-critical systems complexity and bipartite entanglement measures are uniformly bounded, both in system size n and inverse temperature, whereas they exhibit a universal logarithmic divergence in inverse temperature for quantum critical systems. Our results represent a uniform finite temperature extension of the ground state area-law [,, 5]. Specifically, for critical systems, logarithmic divergence of ground state block entropy with the size of the block is reflected in logarithmic divergence in inverse temperature of the bipartite entanglement of thermal (mixed) states in the thermodynamic limit. Even though our results are not rigorous they suggest optimal bounds on entanglement and complexity of thermal states that are qualitatively stronger than previously known [5]. We thank J. Eisert for useful comments and acknowledge support by the grants P-44 and J-7347 of Slovenian Research Agency.
5 [] K. Audenaert et al., Phys. Rev. A 66, 437 (). [] G. Vidal et al., Phys. Rev. Lett. 9, 79 (3); J. I. Latorre et al., Quant. Inf. Comp.4, 48 (4). [3] P. Calabrese and J. Cardy, J. Stat. Mech.: Theor. Exp. P46 (4). [4] L. Amico et al., Rev. Mod. Phys. 8, 57 (8). [5] S. R. White, Phys. Rev. Lett. 69, 863 (99). [6] G. Vidal, Phys. Rev. Lett. 9, 479 (3). [7] F. Verstraete et al., Phys. Rev. Lett. 93, 74 (4); M. Zwolak and G. Vidal, Phys. Rev. Lett. 93, 75 (4); A. E. Feiguin and S. R. White, Phys. Rev. B 7, 4(R) (5). [8] T. Prosen and M. Žnidarič, Phys. Rev. E 75, 5(R) (7). [9] M. Žnidarič et al., Phys. Rev. B 77, 6446 (8). [] J. Eisert, Phys. Rev. Lett. 97, 65 (6). [] M. B. Hastings, Phys. Rev. B 76, (R) (7). [] M. B. Hastings, Phys. Rev. B 73, 855 (6). [3] T. Prosen and I. Pižorn, Phys. Rev. A 76, 336 (7). [4] B. Groisman et al., Phys. Rev. A 7, 337 (5). [5] O. Bratteli and D. W. Robinson, Operator Algebras and Quantum Statistical Mechanics II, Springer (98); M. M. Wolf et al., Phys. Rev. Lett. 75 (8). [6] V. E. Korepin, Phys. Rev. Lett. 9, 964 (4). [7] T. Prosen, New J. Phys, 436 (8). [8] N. Schuch et al., Phys. Rev. Lett., 354 (8). [9] P. Zanardi, Phys. Rev. A 63, 434(R) (); M. A. Nielsen et al., Phys. Rev. A 67, 53 (3). [] V. Vedral and M. B. Plenio, Phys. Rev. A 57, 69 (998). [] J. Kurmann et al., Physica A, 35 (98).