1 / 34

Quantum Spin Systems from the point a view of Quantum Information Theory

Quantum Spin Systems from the point a view of Quantum Information Theory. Frank Verstraete, Ignacio Cirac Max-Planck-Institut f ür Quantenoptik. Overview. The quantum repeater, entanglement length and quantum computing with projected entangled pair states (PEPS)

landry
Télécharger la présentation

Quantum Spin Systems from the point a view of Quantum Information Theory

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Quantum Spin Systems from the point a view of Quantum Information Theory Frank Verstraete, Ignacio Cirac Max-Planck-Institut für Quantenoptik

  2. Overview • The quantum repeater, entanglement length and quantum computing with projected entangled pair states (PEPS) • MPS/PEPS: basic properties • PEPS as variational trial states • Illustration: • RG • DMRG • Applications: • DMRG with Periodic boundary conditions • Calculating Excitations • Dynamical correlation functions using PEPS • Optimal time evolution algorithm • Finite-T DMRG: • Projected entangled pair density operators • quantum-TMRG from the perspective of PEPS • Quantum Spin Glasses • Simulation of 2-D systems • Ground states • Time evolution • Conclusion

  3. Number of states kept (D)

  4. (a) Magnetization vs. transverse magnetic field of the thermal state of the XY model with 60 spins, for temperatures T=0.05,0.5,5 and 50 (bottom to top). (b) Error in the density matrix of the thermal state computed using MPDO, vs. temperature, for a chain with N=8 spins, D=8 (solid) and D=14,20,24 (dashed).

  5. 10x10 20x20 4x4

  6. Spin systems in Quantum Information Theory • Central motto: Entanglement is a resource • Essential ingredient in quantum cryptography, quantum computing, etc. • Basic Unit: Bell-state |00i+|11i • Long-range distribution of entanglement can be done via a quantum repeater : • The effective Hamiltonian describing this system is a spin system with nearest neighbour couplings where we have complete local control over spins

  7. Given a 1-D state of N spins, what is the amount of entanglement that can be localized between separated spins (qubits) in function of their distance by doing local measurements on the other ones? • Leads to notion of Localizable Entanglement (LE) and associated Entanglement length • The localizable entanglement between two spins is always larger than or equal to the connected 2-point correlation functions • Hence the Correlation length is a lower bound to the Entanglement length: long-range correlations imply long-range entanglement • The LE can be calculated using a combination of DMRG and Monte Carlo • Natural question: do there exist ground states with diverging Entanglement length but finite Correlation length? Verstraete, Popp, Cirac 04

  8. Singlet Proj. in sym. subspace The spin-1 AKLT-model • All correlation functions decay exponentially • The symmetric subspace is spanned by 3 Bell states, and hence this ground state can be used as a perfect quantum repeater • Diverging entanglement length but finite correlation length • LE detects new kind of long range order • Antiferro spin-1 chain is a perfect quantum channel Verstraete, Martin-Delgado, Cirac 04

  9. Pi+1 Pi+2 Pi+3 Pi+4 Pi+5 Pi+6 Pi Generalizing the AKLT-state: PEPS • Every state can be represented as a Projected Entangled Pair State (PEPS) as long as D is large enough • Extension to mixed (finite-T) states: take Completely Positive Maps (CPM) instead of projectors : • 1-D PEPS reduce to class of finitely correlated states / matrix product states (MPS) in thermodynamic limit (N!1) when P1=P2=L =P1 • Systematic way of constructing translational invariant states • MPS become dense in space of all states when D!1 Ex.: 5 qubit state Fannes, Nachtergaele, Werner 92

  10. PEPS in higher dimensions: • PEPS in QIT: • Quantum Error correction: • Quantum computing Verstraete, Cirac 04

  11. Basic properties of PEPS • Correlation functions for 1-D PEPS can easily be calculated by multiplying transfer matrices of dimension D2 : • Number of parameters grows linearly in number of particles (NdDc) with c coordination number of lattice • 2-point correlations decay exponentially • Holographic principle: entropy of a block of spins is proportional to its surface

  12. Spin systems: basic properties • Hilbert space scales exponentially in number of spins • Universal ground state properties: • Entropy of block of spins / surface of block (holographic principle) • Correlations of spins decay typically exponentially with distance (correlation length) • The N-particle states with these properties form a tiny subspace of the exponentially large Hilbert space

  13. Ground states are extreme points of a convex set: • Problem of finding ground state energy of all nearest-neighbor transl. invariant Hamiltonians is equivalent to characterizing the convex set of n.n. density operators arising from transl. invariant states • Finitely Correlated States / Matrix Product States / Projected Entangled Pair States provide parameterization that only grows linearly with number of particles but captures these desired essential features very well The Hamiltonian defines a hyperplane in (2s+1)2 dim. space

  14. Pi+1 Pi+2 Pi+3 Pi+4 Pi+5 Pi+6 Pi PEPS as variational trial states • All expectation values and hence the energy are multi-quadratic in the variables Pk • Strategy for minimizing energy for N-spin state: • Fix all projectors Pi except the jth • Both the energy and the norm are quadratic functions of the variable Pj and hence the minimal energy by varying Pi can be obtained by a simple generalized eigenvalue problem: Heff and N are function of the Hamiltonian and all other projectors, and can efficiently be calculated by the transfer matrix method • Move on to the (j§1)th particle and repeat previous steps (sweep) until convergence • Compare to MPS-approach of Ostlund and Rommer 95 and Dukelsky et al. 97 : imposing that all Pi are equal makes the optimization untractable for large D Verstraete, Porras, Cirac 04

  15. J P0 P1 J l1 P1 P2 P3 P4 P5 P6 P0 P1 P2 P0 J l1 l2 l3 l5 l6 l4 Illustration 1 L • Wilson’s Renormalization Group (RG) for Kondo-effect: • RG approximates ground state and lowest excited states by a family of PEPS • Ai are chosen s.t. |yai form orthonormal set, which forces N=identity • Remark that sweeping would give better precision • Excited states are forced to have the same Aik<N , which cannot be true in general; OK for purpose of calculating GS, but better can be done if one wants exact info about spectrum /dynamical correlation functions (see later)

  16. Illustration 2: DMRG White 92 • The PEPS-method coincides with the B B DMRG-scheme if PEPS with OBC are chosen, and assures that the energy is monotonously decreasing • Open boundary conditions enable big speed-up as the calculation of Heff can be done using sparse matrices (computational cost / D3); this follows from the fact that |yai can be chosen orthonormal which implies N=I.

  17. Applications • DMRG with Periodic boundary conditions • Calculating Excitations • Dynamical correlation functions / structure factors using PEPS • Optimal time evolution algorithms • Finite-T DMRG: • Matrix product density operators • Quantum-TMRG from the perspective of PEPS • Quantum Spin Glasses • Simulation of 2-D systems • Ground states • Time evolution, dynamical correlation functions, …

  18. DMRG and periodic boundary conditions • Due to boundary effects, periodic boundary conditions (PBC) are preferred above open ones (OBC) in calculating magnetization, (dynamical correlation functions, excitations, … • Standard DMRG does not work well in case of PBC: • Underlying PEPS trial state has open boundary conditions • entropy of block in B B -scheme scales as (compare to in B B or B B) • Solution: solve generalized eigenvalue problem on PEPS of the form • Remark: sparseness of matrix multiplication partially lost due to impossibility of orthonormalizing basis (simulation time scales as D5 instead of D3 ) • Exact translationnally invariant state can be found by using trial PEPS of the form P1 P2 P3 P4 P5 L PN Verstraete, Porras, Cirac 04

  19. Numerical results in case of spin ½ Heisenberg antiferromagnet (N=28): Number of states kept (D) Verstraete, Porras, Cirac 04

  20. Excitations • Variational formulation: find the PEPS |y1i minimizing the energy under the constraint that it is orthogonal to |y0i • As the constraint is multilinear in the projectors Pi , |y1i can again iteratively be optimized by solving generalized eigenvalue problems; the only difference is that the dD2 matrices Heff and N are projected onto a dD2-1 subspace. • This allows to find excitation energies to same accuracy as for ground states • In case of PBC, trial states with definite momentum can be chosen as • Note that in case of DMRG with multi-state targeting, all states have same Ai except at the site one is looking at; the big difference with the method we propose is that here ALL Ai are varied at the expense of repeating the algorithm for each excited state again

  21. Dynamical Correlation Functions • Goal: calculate • Variational formulation within class of PEPS: find (unnormalized) PEPS |ci such that ||A|y0i-(E0-H+w+ih)|ci||2is minimized; then G(w)=hy0|A*|ci • This optimization can again efficiently be done by an iterative method similar to the one already described: fix all projectors but one, then the norm is given by xy Heffx-xy y ; the solution is then Heffx=y/2. • Main difference with correction vector method / D-DMRG: • |y0i and |ci are calculated independently and are not forced to have the same projectors Pi everywhere except at the site under consideration • This ensures better precision as the variational class of trial states is much larger; note also that the iterative method is efficient as only a set of linear equations has to be solved (scaling in case of OBC: D3) Hallberg 95; Kuhner and White 99; Jeckelmann 02

  22. Pi+1 Pi+2 Pi+3 Pi+4 Pi+5 Pi+6 Pi Optimal Time Evolution • Goal: given a PEPS |yi (OBC or PBC), describe its time evolution via Trotterization under the action of a local Hamiltonian • To describe this 2-step evolution exactly, the dimension of all the bonds has to be multiplied by an interaction-dependent factor (2 for Ising and 4 for Heisenberg) • To prevent the dimensions to blow up: variational dimensional reduction of PEPS • Given a PEPS |yDi of dimension D, find the one |cD’i of dimension D’< D such that || |yDi-|cD’i ||2 is minimized • This can again be done efficiently in a similar iterative way, yielding a variational and hence optimal way of treating time-evolution within the class of PEPS Vidal 03; Daley et al. 04; White and Feiguin 04 Verstraete, Garcia-Ripoll, Cirac 04

  23. Variational Dimensional Reduction of PEPS • Given a PEPS |yDi parameterized by the D£D matrices Ai, find the one |cD’i parameterized by D’£D’matrices Bi (D’< D) such as to minimize • Fixing all Bi but one to be optimized, this leads to an optimization of the form xy Heffx-xy y , with solution: Heffx=y/2 ; iterating this leads to global optimum • The error of the truncation can exactly be calculated at each step! • In case of OBC: more efficient due to orthonormalization • Note that the algorithm also applies to PBC • In the case of OBC, the algorithms of Vidal, Daley et al., White et al. are suboptimal but a factor of 2-3 times faster; a detailed comparison should be made

  24. Finite-T DMRG • PEPS can represent mixed states by applying completely positive maps (CPM) instead of projectors to the virtual particles • The condition on the M’s ensures positivity; general operators are obtained by relaxing this condition • This PEPS-picture is again complete if dimension of bonds large enough • It turns out to be much more efficient and convenient to work with the purification of r; r is obtained by tracing over the auxiliary systems ak : See also: Zwolak and Vidal 04 Verstraete, Garcia-Ripoll, Cirac 04

  25. Time-evolution, calculation of dynamical correlation functions, … can now be done on this purification as the physical observables act trivially on the auxiliary systems • Calculation of a thermal state: imaginary time evolution starting from the 1-T state: • The dimension of the bonds of the purification are only a factor d (e.g. 2 in case of qubits) larger than in the zero-T case • In the case of decoherence (simulation of master equations), evolution of r instead of its purification has to be done; the computational cost in case of OBC then scales as D6 instead of d3D3 and positivity of r is not anymore assured. • a good alternative: the quantum jump approach See also: Zwolak and Vidal 04

  26. (a) Magnetization vs. transverse magnetic field of the thermal state of the XY model with 60 spins, for temperatures T=0.05,0.5,5 and 50 (bottom to top). (b) Error in the density matrix of the thermal state computed using MPDO, vs. temperature, for a chain with N=8 spins, D=8 (solid) and D=14,20,24 (dashed).

  27. a1 a1 a2 a2 a3 a3 a4 a4 a5 a5 a6 a6 a7 a7 a8 a8 M i1 M i1 L L i2 M i2 L L i3 M i3 L L i4 M i4 L L i5 M i5 L L quantum-TMRG from the perspective of PEPS Nishino 95; Bursill et al. 96; Shibata 97; Wang et al. 97 • Different approach to finite-T DMRG: calculate thermodynamic quantities like the free energy • The operations needed to take this trace can be represented graphically: i0 M N i0 M

  28. i0 M N i0 M a1 a1 a2 a2 a3 a3 a4 a4 a5 a5 a6 a6 a7 a7 a8 a8 M i1 M i1 L L i2 M i2 L L i3 M i3 L L i4 M i4 L L i5 M i5 L L • 2 possible methods to sum this: • If N!1 and homogeneous, one can calculate the largest eigenvector of a column transfer matrix (TMRG) using DMRG with PERIODIC boundary condition (note that it is not always possible to get hermitean transfer matrices) • More general approach (also for finite, inhomogeneous and non-hermitean case): propagate a 1-D PEPS from the left to the right by the method of variational dimensional reduction • These methods also work for calculating thermodynamic quantities of classical 2-D spin systems (cfr. Nishino); here M=N!1 and an infinite-DMRG like algorithm can be devised that gives very accurate results even at the critical point

  29. Quantum Spin Glasses • We can use the inherent randomness of quantum mechanics to study spin systems with random interactions / fields : encode the randomness in ancilla’s. • This allows to run all possible realizations in 1 run • Simple example: • Simulation Hamiltonian :

  30. Simulation of 2-D quantum systems • Standard DMRG approach: trial state of the form Problems with this approach: dimension of bonds must be exponentially large: - area theorem - only possibility to get large correlations between vertical nearest neighbors We propose trial PEPS states that have bonds between all nearest neighbors, such that the area theorem is fulfilled by construction and all neighbors are treated on equal footing See also Nishino et al.

  31. P12 P13 P14 P15 P11 P22 P24 P25 P21 P23 P31 P32 P34 P35 P33 P45 P31 P44 P41 P42 • The energy of such a state is still a multi-quadratic function of all variable, and hence the same iterative variational principle can be used • The big difference: the determination of Heff and N is not obtained by multiplying matrices, but contracting tensors • This can be done using the variational dimensional reduction discussed in the context of TMRG; note that the error in the truncation is completely controlled No (sign) problem with frustrated systems! Possible to devise an infinite dimensional variant

  32. Alternatively, the ground state can be found by imaginary time evolution on a pure 2-D PEPS • This can be implemented by Trotterization; the crucial ingredient is again the variational dimensional reduction; the computational cost scales linearly in the number of spins: D10 • The same algorithm can of course be used for real-time evolution and for finding thermal states. • Dynamical correlation functions can be calculated as in the 1-D PEPS case • We have done simulations with the Heisenberg antiferromagnetic interaction and a frustrated version of it on 4£4, 10£10 and 20£20 • We used bonds of dimension 2,3,4; the error seems to decay exponentially in D • Note that we get mean field if D=1 • The number of variational parameters scales as ND4 and we expect the same accuracy as 1-D DMRG with dimension of bonds D2

  33. 20x20 4x4: 36.623 10x10: 2.353 (D=2); 2.473 (D=3) 20x20: 2.440 (D=2); 2.560 (D=3)

  34. Conclusion • Projected Entangled Pair States (PEPS) provide a natural parameterization of multipartite entangled states on lattices from the point a view of Quantum Information Theory • PEPS provide a solid theoretical basis underlying DMRG, allowing to generalize DMRG to periodic boundary conditions, finite-T, time evolution, dynamical correlation functions, excitations, quantum spin glasses, and most notably higher dimensions

More Related