Time Evolving Block Decimation Algorithm
|
|
- Clarence Riley
- 5 years ago
- Views:
Transcription
1 Time Evolving Block Decimation Algorithm Application to bosons on a lattice Jakub Zakrzewski Marian Smoluchowski Institute of Physics and Mark Kac Complex Systems Research Center, Jagiellonian University, Kraków, Poland November 2006
2 What will it be about Intro Motivation Example: Bose-Hubbard Time-Evolving Block Decimation Entanglement Vidal s construction Time evolution algorithm Finding the ground state Costs and Improvements Bose-Hubbard example Extension to an infinite chain Bose-Hubbard example Further reading
3 Motivation Many interacting particles systems (bosons, fermions, spins - typical problems of statistical physics, condensed matter, e.g. Ising model) Static properties, ground state, low excitations, quantum phase transitions Dynamics, time evolution,... i ψ = H ψ t Typowe metody Typical approaches: Finite difference method Basis sets expansions (or in overcomplete coherent states..) Large scale diagonalizations but...
4 Example: Bose-Hubbard Hamiltonian Ĥ = J d i=0 ( ) â i+1âi + h.a. + d ɛ iˆn i + U 2 i=0 d ˆn i (ˆn i 1), i=1 â i (â i ) are creation (annihilation) operators for a boson at site i ˆn i = â i âi occupation number operator J - hopping rate, U - interaction energy, d maximal occupation ( dim(h) = d + 1 on each site) Very sparse matrix of dimension: D = (N + M 1)! N!(M 1)!, the number of particles is N and number of lattice sites M. D = 6435 for N = M = 8 D = for N = M = 10 D = for N = M = 12
5 Bose-Hubbard Hamiltonian possibilities: Basis truncation Lanczos scheme s...
6 Bose-Hubbard Hamiltonian possibilities: Basis truncation Lanczos scheme s... yield a partial gain only.. Alternatives for a static (ground state) problem Quantum Monte Carlo (QMC) Density Matrix Renormalization Group (DMRG) approximations (mean field etc.)... as well as novel ideas coming from Quantum Entanglement machinery
7 Quantum Entanglement Primer Consider a system composed of two parts iff Ψ = Φ [A] Φ [B] the state is separable otherwise Schmidt s decomposition possible χ A Ψ = λ α Φ [A] α Φ [B] α, Φ [A] α ( Φ [B] ρ [A] = Tr B ( Ψ Ψ α=1 α ) is an eigenvector with eigenvalue λ α 2 > 0 of (ρ [B] ) also Φ [A] α Ψ = λ α Φ [B] α. The entanglement of Ψ is e.g. χ max A χ A. We assume χ is small.
8 Entanglement Primer 2 d Ψ = C ij i A j B Then Singular Value Decomposition does the job i,j C = ÛˆDV ; ˆD = diagλ i Now we shall concentrate on states with hopefully small entanglement. Also the type of time-evolution that will slowly increase χ only.
9 Vidal s construction Consider now a general state for M-site chain: d d Ψ = c i1 i n i 1 i n, i 1 =0 i n=0
10 Vidal s construction Consider now a general state for M-site chain: d d Ψ = c i1 i n i 1 i n, i 1 =0 i n=0 The key ingredient of Vidal s simulation protocol: Ψ Γ [1] λ [1] Γ [2] λ [2] Γ [l] λ [n 1] Γ [n].
11 Vidal s construction Consider now a general state for M-site chain: d d Ψ = c i1 i n i 1 i n, i 1 =0 i n=0 The key ingredient of Vidal s simulation protocol: Ψ Γ [1] λ [1] Γ [2] λ [2] Γ [l] λ [n 1] Γ [n]. Γ [l]i αα - tensor associated with site l λ [l] α - the Schmidt splitting [1 l]:[(l +1) n]. explicite :-) c i1 i 2 i n = Γ [1]i 1 α 1 α 1,,α n 1 λ [l] α 1 Γ [2]i 2 α 1 α 2 λ [2] α 2 Γ [n]in α n 1.
12 Vidal s construction 2 How to get Γ s (i) expand Φ [2 n] α Ψ = α 1 λ [1] α 1 Φ [1] α 1 Φ [2 n] α 1 = i 1,α 1 Γ [1]i 1 α 1 λ [1] α 1 i 1 Φ [2 n] α 1, Φ [2 n] α 1 = i 2 i 2 τ [3 n] α 1 i 2 ; (ii) express via Schmidt decomposition of 2 : [3..n] (iii) make all substitutions τ [3 n] α 1 i 2 = α 2 Γ [2]i 2 α 1 α 2 λ [2] α 2 Φ [3 n] α 2 ; Ψ = Γ α [1]i1 1 λ α [1] 1 Γ [2]i 2 α 1 α 2 λ [2] α 2 i 1 i 2 Φ [3 n] α 1. i 1,α 1,i 2,α 2
13 (iv) move to site 4,5,...
14 Remarks: close links to finitely correlated states (Fannes et al.) close links to matrix product states (Öslund and Rommer) thus very similar to DMRG scheme
15 Remarks: close links to finitely correlated states (Fannes et al.) close links to matrix product states (Öslund and Rommer) thus very similar to DMRG scheme So what is new?
16 Remarks: close links to finitely correlated states (Fannes et al.) close links to matrix product states (Öslund and Rommer) thus very similar to DMRG scheme So what is new? simple study of time dependence
17 Remarks: close links to finitely correlated states (Fannes et al.) close links to matrix product states (Öslund and Rommer) thus very similar to DMRG scheme So what is new? simple study of time dependence similar approach for ground state solution
18 Remarks: close links to finitely correlated states (Fannes et al.) close links to matrix product states (Öslund and Rommer) thus very similar to DMRG scheme So what is new? simple study of time dependence similar approach for ground state solution natural extension to truly infinite systems
19 Time evolution Consider a Hamiltonian H = M i=1 H [i] 1 + M i=1 H [i,i+1] 2 ;
20 Time evolution Consider a Hamiltonian H = M i=1 H [i] 1 + M i=1 H [i,i+1] 2 ; Time evolution: Ψ(T ) = exp( iht ) Ψ(0) done, as usual, in small steps Ψ(t + τ) = exp( ihτ) Ψ(t)
21 Time evolution Consider a Hamiltonian H = M i=1 H [i] 1 + M i=1 H [i,i+1] 2 ; Time evolution: Ψ(T ) = exp( iht ) Ψ(0) done, as usual, in small steps Ψ(t + τ) = exp( ihτ) Ψ(t) Decompose H as H = F + G, F F [l] even l even l G G [l] odd l odd l (H [l] 1 + H[l,l+1] 2 ), (H [l] 1 + H[l,l+1] 2 ),
22 Time evolution Consider a Hamiltonian H = M i=1 H [i] 1 + M i=1 H [i,i+1] 2 ; Time evolution: Ψ(T ) = exp( iht ) Ψ(0) done, as usual, in small steps Ψ(t + τ) = exp( ihτ) Ψ(t) Decompose H as H = F + G, F F [l] even l even l G G [l] odd l odd l (H [l] 1 + H[l,l+1] 2 ), (H [l] 1 + H[l,l+1] 2 ), Observe: [F [l], F [l ] ] = 0 ([G [l], G [l ] ] = 0) but [F, G] 0.
23 Time evolution 2 Apply Trotter expansion of order p for exp( ihτ) exp( i(f + G)τ) f p [exp( if τ), exp( igτ)], f 1 (x, y) = xy, f 2 (x, y) = x 1/2 yx 1/2,
24 Time evolution 2 Apply Trotter expansion of order p for exp( ihτ) exp( i(f + G)τ) f p [exp( if τ), exp( igτ)], f 1 (x, y) = xy, f 2 (x, y) = x 1/2 yx 1/2, And we get a chain of two-site gates: exp( if τ) = Π l exp( if [l] τ) = Π l exp( i(h [l] 1 +H[l,l+1] 2 )τ) = Π l U [l] Vidal s lemma: Updating Ψ expressed as Ψ Γ [1] λ [1] Γ [2] λ [2] Γ [l] λ [n 1] Γ [n]. is a LOCAL operation affecting Γ [l], λ [l] and Γ [l+1].
25 Time evolution 3
26 Time evolution 3 Recall that there are at most χ links of our l and l + 1 sites χ d Ψ = αijγ, α,β,γ=1 i,j=0 Γ [l]i αβ λ[l] β Γ[l+1]j βγ where Ψ = U (l) Ψ = Θ ij αγ = β kn χ α,γ=1 i,j=0 d Θ ij αγ αijγ, V ij kn Γ[l]k αβ λ[l] β Γ[l+1]n βγ.
27 Time evolution 4 Now singular value decomposition of Θ grouping indexes as [αi] : [jγ]: Θ = β X (l) [αi]β λ β Y β[jγ]
28 Time evolution 4 Now singular value decomposition of Θ grouping indexes as [αi] : [jγ]: Θ = β X (l) [αi]β λ β Y β[jγ] And Γ [l]i αβ = X [αi]β/λ (l 1) α ; Γ [l+1]i αβ = Y [αi]β /λ (l+1) β
29 Time evolution 4 Now singular value decomposition of Θ grouping indexes as [αi] : [jγ]: Θ = β X (l) [αi]β λ β Y β[jγ] And Γ [l]i αβ = X [αi]β/λ (l 1) α ; Γ [l+1]i αβ = Y [αi]β /λ (l+1) β Note: α [1, χ], i.e. [αi] leads to χ(d + 1) as dimension of Θ.
30 Time evolution 4 Now singular value decomposition of Θ grouping indexes as [αi] : [jγ]: Θ = β X (l) [αi]β λ β Y β[jγ] And Γ [l]i αβ = X [αi]β/λ (l 1) α ; Γ [l+1]i αβ = Y [αi]β /λ (l+1) β Note: α [1, χ], i.e. [αi] leads to χ(d + 1) as dimension of Θ. SVD yields λ [l] of that length, normalized so as ( λ[l] α ) 2 = 1.
31 Time evolution 4 Now singular value decomposition of Θ grouping indexes as [αi] : [jγ]: Θ = β X (l) [αi]β λ β Y β[jγ] And Γ [l]i αβ = X [αi]β/λ (l 1) α ; Γ [l+1]i αβ = Y [αi]β /λ (l+1) β Note: α [1, χ], i.e. [αi] leads to χ(d + 1) as dimension of Θ. SVD yields λ [l] of that length, normalized so as ( λ[l] α ) 2 = 1. Numerically we want to limit the length of λ to some χ e so we drop smallest and renormalize: α = [ λ [l] χ ɛ α=1 And that s the end of the story!! λ [l] α 2 ] 1 2 λ [l] α
32 Finding the ground state (one of the methods) Start with Φ φ [1] φ [n], Φ Ψ gr 0, (expansion into λ s and Γ s) trivial.
33 Finding the ground state (one of the methods) Start with Φ φ [1] φ [n], Φ Ψ gr 0, (expansion into λ s and Γ s) trivial. Use the scheme to simulate an evolution in imaginary time t according to H, Ψ gr = lim t exp( Ht) Φ exp( Ht) Φ ;
34 Finding the ground state (one of the methods) Start with Φ φ [1] φ [n], Φ Ψ gr 0, (expansion into λ s and Γ s) trivial. Use the scheme to simulate an evolution in imaginary time t according to H, Ψ gr = lim t exp( Ht) Φ exp( Ht) Φ ; The procedure converges almost always :-)
35 Costs and Improvements Generally execution time scales as ((d + 1)χ e ) 3 Further simple improvements:
36 Costs and Improvements Generally execution time scales as ((d + 1)χ e ) 3 Further simple improvements: Global selection rules ((d + 1)χ e ) 3 (d eff χ e ) 3
37 Costs and Improvements Generally execution time scales as ((d + 1)χ e ) 3 Further simple improvements: Global selection rules ((d + 1)χ e ) 3 (d eff χ e ) 3 The procedure is easy to parallelize :-)
38 Costs and Improvements Generally execution time scales as ((d + 1)χ e ) 3 Further simple improvements: Global selection rules ((d + 1)χ e ) 3 (d eff χ e ) 3 The procedure is easy to parallelize :-) But yet to be done :-( Further difficult (under construction) improvements Long range interactions Renormalizable entanglement Tree structures multidimensions - - projected entanged pair states, MERA
39 Bose-Hubbard example Comparison with high order perturbation theory (1) Ĥ = J d i=0 ( ) â i+1âi + h.a. + d ɛ iˆn i + U 2 i=0 d ˆn i (ˆn i 1), i=1 1 B. Damski and J. Zakrzewski, Phys. Rev. A74, (2006)
40 Bose-Hubbard example Comparison with high order perturbation theory (1) Ĥ = J d i=0 ( ) â i+1âi + h.a. + d ɛ iˆn i + U 2 i=0 d ˆn i (ˆn i 1), i=1 1 B. Damski and J. Zakrzewski, Phys. Rev. A74, (2006)
41 Bose-Hubbard example Comparison with high order perturbation theory (1) Ĥ = J d i=0 ( ) â i+1âi + h.a. + d ɛ iˆn i + U 2 i=0 d ˆn i (ˆn i 1), i=1 1 B. Damski and J. Zakrzewski, Phys. Rev. A74, (2006)
42 Bose-Hubbard example Comparison with high order perturbation theory (1) Ĥ = J d i=0 ( ) â i+1âi + h.a. + d ɛ iˆn i + U 2 i=0 d ˆn i (ˆn i 1), i=1 1 B. Damski and J. Zakrzewski, Phys. Rev. A74, (2006)
43 Treatment of infinite systems Let H = r H[r,r+1] be invariant under r r + 1 transformation. If Ψ 0 is also translationally invariant, evolution preserves that. Now we can treat truly infinite systems χ Ψ = λ (r) where Φ ( r+1) α = Ψ = χ Φ ( r) α = D α,β,γ=1 i,j=1 α=1 χ β=1 i=1 χ α Φ ( r) α D D β=1 i=1 λ (r 1) α Γ (r) iαβ λ(r) Φ ( r+1) α ; λ (r) β Γ(r+1) iβα Φ( r) Γ (r+1) iαβ β Γ(r+1) jβγ λ (r+1) γ β i (r+1), λ(r+1) β i (r) Φ ( r+1) β. i (r) j (r+1) Φ α ( r) Φ γ ( r+2).
44 Infinite systems 2 Ψ shift invariant (r) is not needed, but.. We need to act with the two-side gate: U [r,r+1] exp( ih [r,r+1] δt), δt 1 that breaks for a moment the symmetry. Thus Γ [2r] = Γ A, λ [2r] = λ A, Γ [2r+1] = Γ B, λ [2r+1] = λ B, r Z.
45 Infinite systems 2 Ψ shift invariant (r) is not needed, but.. We need to act with the two-side gate: U [r,r+1] exp( ih [r,r+1] δt), δt 1 that breaks for a moment the symmetry. Thus Γ [2r] = Γ A, λ [2r] = λ A, Γ [2r+1] = Γ B, λ [2r+1] = λ B, r Z.
46 Infinite systems 3
47 Bose-Hubbard example Ĥ = J d i=0 ( ) â i+1âi + h.a. + d ɛ iˆn i + U 2 i=0 d ˆn i (ˆn i 1) µ i=1 d ˆn i, Note the presence of the chemical potential µ. Challenging problem: The tip of the quantum phase transitions i=1
48 Bose-Hubbard example 2 Various methods used (mean-field, Bethe-ansatz, diagonalization N = M = 12, QMC, DMRG...) The most stable approach - the slope of the correlation function C(s) = Ψ gr â r+sâr Ψ gr. exponential decay for large s in Mott Insulator Phase. power law decay for large s in Superfluid Phase C(s) s K /2 with K = 1/2 at transition.
49 Bose-Hubbard example 3 T. Kühner, S.R. White, H. Monien, Phys. Rev. B61, (2000) give J c = ± 0.01 (DMRG). Our simulation (2) yields J c = ± Table: The location of the critical point J c for different ranges of s value in the correlation function C(s) s KWM2000 This work 4 s ± ± s ± ± s ± ± s ± ± s ± ± s ± s ± J. Zakrzewski and D. Delande, in preparation (2006)
50 Bose-Hubbard example 4 And that is converged (thanks to ICM...)
51 Further reading the basic algorithm: G. Vidal Phys. Rev. Lett. 91, (2003). the infinite case: G. Vidal cond-mat/ relation to DMRG A.J. Daley et al. cond-mat/ ; DMRG review U. Schollwoeck, Rev. Mod. Phys. 77, 259 (2005) and cond-mat/ attempts for multidimensions F. Verstraete and J.I. Cirac cond-mat/ (PEPS); G. Vidal cond-mat/ Entanglement renormalization: G. Vidal, cond-mat/
IPAM/UCLA, Sat 24 th Jan Numerical Approaches to Quantum Many-Body Systems. QS2009 tutorials. lecture: Tensor Networks.
IPAM/UCLA, Sat 24 th Jan 2009 umerical Approaches to Quantum Many-Body Systems QS2009 tutorials lecture: Tensor etworks Guifre Vidal Outline Tensor etworks Computation of expected values Optimization of
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 informationQuantum many-body systems and tensor networks: simulation methods and applications
Quantum many-body systems and tensor networks: simulation methods and applications Román Orús School of Physical Sciences, University of Queensland, Brisbane (Australia) Department of Physics and Astronomy,
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 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 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 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 informationarxiv: v4 [cond-mat.stat-mech] 13 Mar 2009
The itebd algorithm beyond unitary evolution R. Orús and G. Vidal School of Physical Sciences, The University of Queensland, QLD 4072, Australia arxiv:0711.3960v4 [cond-mat.stat-mech] 13 Mar 2009 The infinite
More informationMany-Body physics meets Quantum Information
Many-Body physics meets Quantum Information Rosario Fazio Scuola Normale Superiore, Pisa & NEST, Istituto di Nanoscienze - CNR, Pisa Quantum Computers Interaction between qubits two-level systems Many-Body
More informationPositive Tensor Network approach for simulating open quantum many-body systems
Positive Tensor Network approach for simulating open quantum many-body systems 19 / 9 / 2016 A. Werner, D. Jaschke, P. Silvi, M. Kliesch, T. Calarco, J. Eisert and S. Montangero PRL 116, 237201 (2016)
More informationRecent developments in DMRG. Eric Jeckelmann Institute for Theoretical Physics University of Hanover Germany
Recent developments in DMRG Eric Jeckelmann Institute for Theoretical Physics University of Hanover Germany Outline 1. Introduction 2. Dynamical DMRG 3. DMRG and quantum information theory 4. Time-evolution
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 informationTuning the Quantum Phase Transition of Bosons in Optical Lattices
Tuning the Quantum Phase Transition of Bosons in Optical Lattices Axel Pelster 1. Introduction 2. Spinor Bose Gases 3. Periodically Modulated Interaction 4. Kagome Superlattice 1 1.1 Quantum Phase Transition
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 informationTime-dependent variational principle for quantum many-body systems
Quantum Information in Quantum Many-Body Physics October 21, 2011 Centre de Recherches Mathematiques, Montréal Time-dependent variational principle for quantum many-body systems PRL 107, 070601 (2011)
More informationFermionic tensor networks
Fermionic tensor networks Philippe Corboz, Institute for Theoretical Physics, ETH Zurich Bosons vs Fermions P. Corboz and G. Vidal, Phys. Rev. B 80, 165129 (2009) : fermionic 2D MERA P. Corboz, R. Orus,
More informationLight Cone Matrix Product
Light Cone Matrix Product M. B. Hastings 1 1 Center for Nonlinear Studies and Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545 We show how to combine the light-cone and matrix
More informationLecture 3: Tensor Product Ansatz
Lecture 3: Tensor Product nsatz Graduate Lectures Dr Gunnar Möller Cavendish Laboratory, University of Cambridge slide credits: Philippe Corboz (ETH / msterdam) January 2014 Cavendish Laboratory Part I:
More informationDephasing, relaxation and thermalization in one-dimensional quantum systems
Dephasing, relaxation and thermalization in one-dimensional quantum systems Fachbereich Physik, TU Kaiserslautern 26.7.2012 Outline 1 Introduction 2 Dephasing, relaxation and thermalization 3 Particle
More informationClassical 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 informationFermionic tensor networks
Fermionic tensor networks Philippe Corboz, ETH Zurich / EPF Lausanne, Switzerland Collaborators: University of Queensland: Guifre Vidal, Roman Orus, Glen Evenbly, Jacob Jordan University of Vienna: Frank
More informationQuantum simulation with string-bond states: Joining PEPS and Monte Carlo
Quantum simulation with string-bond states: Joining PEPS and Monte Carlo N. Schuch 1, A. Sfondrini 1,2, F. Mezzacapo 1, J. Cerrillo 1,3, M. Wolf 1,4, F. Verstraete 5, I. Cirac 1 1 Max-Planck-Institute
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 informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION doi:10.1038/nature10748 Contents 1 Quenches to U /J = 5.0 and 7.0 2 2 Quasiparticle model 3 3 Bose Hubbard parameters 6 4 Ramp of the lattice depth for the quench 7 WWW.NATURE.COM/NATURE
More informationSweep from Superfluid to Mottphase in the Bose-Hubbard model p.1/14
Sweep from Superfluid to phase in the Bose-Hubbard model Ralf Schützhold Institute for Theoretical Physics Dresden University of Technology Sweep from Superfluid to phase in the Bose-Hubbard model p.1/14
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 informationMethodology for the digital simulation of open quantum systems
Methodology for the digital simulation of open quantum systems R B Sweke 1, I Sinayskiy 1,2 and F Petruccione 1,2 1 Quantum Research Group, School of Physics and Chemistry, University of KwaZulu-Natal,
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 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 informationEfficient Numerical Simulations Using Matrix-Product States
Efficient Numerical Simulations Using Matrix-Product States Frank Pollmann frankp@pks.mpg.de Max-Planck-Institut für Physik komplexer Systeme, 01187 Dresden, Germany October 31, 2016 [see Ref. 8 for a
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 informationQuantum dynamics in ultracold atoms
Rather don t use Power-Points title Page Use my ypage one instead Quantum dynamics in ultracold atoms Corinna Kollath (Ecole Polytechnique Paris, France) T. Giamarchi (University of Geneva) A. Läuchli
More informationT ensor N et works. I ztok Pizorn Frank Verstraete. University of Vienna M ichigan Quantum Summer School
T ensor N et works I ztok Pizorn Frank Verstraete University of Vienna 2010 M ichigan Quantum Summer School Matrix product states (MPS) Introduction to matrix product states Ground states of finite systems
More informationDiagrammatic Green s Functions Approach to the Bose-Hubbard Model
Diagrammatic Green s Functions Approach to the Bose-Hubbard Model Matthias Ohliger Institut für Theoretische Physik Freie Universität Berlin 22nd of January 2008 Content OVERVIEW CONSIDERED SYSTEM BASIC
More informationThe Density Matrix Renormalization Group: Introduction and Overview
The Density Matrix Renormalization Group: Introduction and Overview Introduction to DMRG as a low entanglement approximation Entanglement Matrix Product States Minimizing the energy and DMRG sweeping The
More informationA quantum dynamical simulator Classical digital meets quantum analog
A quantum dynamical simulator Classical digital meets quantum analog Ulrich Schollwöck LMU Munich Jens Eisert Freie Universität Berlin Mentions joint work with I. Bloch, S. Trotzky, I. McCulloch, A. Flesch,
More informationTensor operators: constructions and applications for long-range interaction systems
In this paper we study systematic ways to construct such tensor network descriptions of arbitrary operators using linear tensor networks, so-called matrix product operators (MPOs), and prove the optimality
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 informationQuantum Hamiltonian Complexity. Itai Arad
1 18 / Quantum Hamiltonian Complexity Itai Arad Centre of Quantum Technologies National University of Singapore QIP 2015 2 18 / Quantum Hamiltonian Complexity condensed matter physics QHC complexity theory
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 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 informationIntroduction to Tensor Networks: PEPS, Fermions, and More
Introduction to Tensor Networks: PEPS, Fermions, and More Román Orús Institut für Physik, Johannes Gutenberg-Universität, Mainz (Germany)! School on computational methods in quantum materials Jouvence,
More informationA Numerical Implementation of the Time-Evolving Block-Decimation Algorithm
A Numerical Implementation of the Time-Evolving Block-Decimation Algorithm Diploma Thesis Diplomarbeit Cojuhovschi Alex Sergiu angefertigt am Institut für Theoretische Physik Universität Hannover Hannover,
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 informationSqueezing and superposing many-body states of Bose gases in confining potentials
Squeezing and superposing many-body states of Bose gases in confining potentials K. B. Whaley Department of Chemistry, Kenneth S. Pitzer Center for Theoretical Chemistry, Berkeley Quantum Information and
More informationQuantum Simulation. 9 December Scott Crawford Justin Bonnar
Quantum Simulation 9 December 2008 Scott Crawford Justin Bonnar Quantum Simulation Studies the efficient classical simulation of quantum circuits Understanding what is (and is not) efficiently simulable
More informationLoop optimization for tensor network renormalization
Yukawa Institute for Theoretical Physics, Kyoto University, Japan June, 6 Loop optimization for tensor network renormalization Shuo Yang! Perimeter Institute for Theoretical Physics, Waterloo, Canada Zheng-Cheng
More informationPlaquette Renormalized Tensor Network States: Application to Frustrated Systems
Workshop on QIS and QMP, Dec 20, 2009 Plaquette Renormalized Tensor Network States: Application to Frustrated Systems Ying-Jer Kao and Center for Quantum Science and Engineering Hsin-Chih Hsiao, Ji-Feng
More informationIntroduction to tensor network state -- concept and algorithm. Z. Y. Xie ( 谢志远 ) ITP, Beijing
Introduction to tensor network state -- concept and algorithm Z. Y. Xie ( 谢志远 ) 2018.10.29 ITP, Beijing Outline Illusion of complexity of Hilbert space Matrix product state (MPS) as lowly-entangled state
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 informationarxiv: v2 [cond-mat.str-el] 28 Apr 2010
Simulation of strongly correlated fermions in two spatial dimensions with fermionic Projected Entangled-Pair States arxiv:0912.0646v2 [cond-mat.str-el] 28 Apr 2010 Philippe Corboz, 1 Román Orús, 1 Bela
More informationEntanglement, quantum critical phenomena and efficient simulation of quantum dynamics
Entanglement, quantum critical phenomena and efficient simulation of quantum dynamics Simons Conference on Reversible and Quantum Computation Stony Brook, May 28-31 2003 Guifre Vidal Institute for Quantum
More informationTensor network renormalization
Coogee'15 Sydney Quantum Information Theory Workshop Tensor network renormalization Guifre Vidal In collaboration with GLEN EVENBLY IQIM Caltech UC Irvine Quantum Mechanics 1920-1930 Niels Bohr Albert
More informationMatrix-Product states: Properties and Extensions
New Development of Numerical Simulations in Low-Dimensional Quantum Systems: From Density Matrix Renormalization Group to Tensor Network Formulations October 27-29, 2010, Yukawa Institute for Theoretical
More informationQuantum quenches in the thermodynamic limit
Quantum quenches in the thermodynamic limit Marcos Rigol Department of Physics The Pennsylvania State University Correlations, criticality, and coherence in quantum systems Evora, Portugal October 9, 204
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 informationNon-Markovian Quantum Dynamics of Open Systems: Foundations and Perspectives
Non-Markovian Quantum Dynamics of Open Systems: Foundations and Perspectives Heinz-Peter Breuer Universität Freiburg QCCC Workshop, Aschau, October 2007 Contents Quantum Markov processes Non-Markovian
More informationarxiv: v1 [quant-ph] 18 Jul 2017
Implicitly disentangled renormalization arxiv:1707.05770v1 [quant-ph] 18 Jul 017 Glen Evenbly 1 1 Département de Physique and Institut Quantique, Université de Sherbrooke, Québec, Canada (Dated: July 19,
More informationMachine Learning with Quantum-Inspired Tensor Networks
Machine Learning with Quantum-Inspired Tensor Networks E.M. Stoudenmire and David J. Schwab Advances in Neural Information Processing 29 arxiv:1605.05775 RIKEN AICS - Mar 2017 Collaboration with David
More informationM AT R I X P R O D U C T S TAT E S F O R L AT T I C E G A U G E T H E O R I E S. kai zapp
M AT R I X P R O D U C T S TAT E S F O R L AT T I C E G A U G E T H E O R I E S kai zapp October 2015 Kai Zapp: Matrix Product States for Lattice Gauge Theories, October 2015 supervisors: Jun.-Prof. Román
More informationLECTURE 3 WORM ALGORITHM FOR QUANTUM STATISTICAL MODELS
LECTURE 3 WORM ALGORITHM FOR QUANTUM STATISTICAL MODELS LECTURE 3 WORM ALGORITHM FOR QUANTUM STATISTICAL MODELS Path-integral for lattice bosons*: oriented closed loops, of course LECTURE 3 WORM ALGORITHM
More information2D Bose and Non-Fermi Liquid Metals
2D Bose and Non-Fermi Liquid Metals MPA Fisher, with O. Motrunich, D. Sheng, E. Gull, S. Trebst, A. Feiguin KITP Cold Atoms Workshop 10/5/2010 Interest: A class of exotic gapless 2D Many-Body States a)
More informationSpin liquid phases in strongly correlated lattice models
Spin liquid phases in strongly correlated lattice models Sandro Sorella Wenjun Hu, F. Becca SISSA, IOM DEMOCRITOS, Trieste Seiji Yunoki, Y. Otsuka Riken, Kobe, Japan (K-computer) Williamsburg, 14 June
More informationClassical and quantum simulation of dissipative quantum many-body systems
0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 Classical and quantum simulation of dissipative quantum many-body systems
More informationRenormalization of Tensor Network States
Renormalization of Tensor Network States I. Coarse Graining Tensor Renormalization Tao Xiang Institute of Physics Chinese Academy of Sciences txiang@iphy.ac.cn Numerical Renormalization Group brief introduction
More informationLecture 20: Effective field theory for the Bose- Hubbard model
Lecture 20: Effective field theory for the Bose- Hubbard model In the previous lecture, we have sketched the expected phase diagram of the Bose-Hubbard model, and introduced a mean-field treatment that
More informationTensor network renormalization
Walter Burke Institute for Theoretical Physics INAUGURAL CELEBRATION AND SYMPOSIUM Caltech, Feb 23-24, 2015 Tensor network renormalization Guifre Vidal Sherman Fairchild Prize Postdoctoral Fellow (2003-2005)
More informationThe continuum limit of the integrable open XYZ spin-1/2 chain
arxiv:hep-th/9809028v2 8 Sep 1998 The continuum limit of the integrable open XYZ spin-1/2 chain Hiroshi Tsukahara and Takeo Inami Department of Physics, Chuo University, Kasuga, Bunkyo-ku, Tokyo 112-8551
More informationEntanglement spectrum of the 2D Bose-Hubbard model
of the D Bose-Hubbard model V. Alba,M. Haque, A. Läuchli 3 Ludwig-Maximilians Universität, München Max Planck Institute for the Physics of Complex Systems, Dresden 3 University of Innsbruck, Innsbruck
More informationFrom unitary dynamics to statistical mechanics in isolated quantum systems
From unitary dynamics to statistical mechanics in isolated quantum systems Marcos Rigol Department of Physics The Pennsylvania State University The Tony and Pat Houghton Conference on Non-Equilibrium Statistical
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 informationThermodynamics of the polarized unitary Fermi gas from complex Langevin. Joaquín E. Drut University of North Carolina at Chapel Hill
Thermodynamics of the polarized unitary Fermi gas from complex Langevin Joaquín E. Drut University of North Carolina at Chapel Hill INT, July 2018 Acknowledgements Organizers Group at UNC-CH (esp. Andrew
More informationQuantum impurities in a bosonic bath
Ralf Bulla Institut für Theoretische Physik Universität zu Köln 27.11.2008 contents introduction quantum impurity systems numerical renormalization group bosonic NRG spin-boson model bosonic single-impurity
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 informationAccelerating QMC on quantum computers. Matthias Troyer
Accelerating QMC on quantum computers Matthias Troyer International Journal of Theoretical Physics, VoL 21, Nos. 6/7, 1982 Simulating Physics with Computers Richard P. Feynman Department of Physics, California
More informationarxiv: v1 [cond-mat.str-el] 19 Jan 2012
Perfect Sampling with Unitary Tensor Networks arxiv:1201.3974v1 [cond-mat.str-el] 19 Jan 2012 Andrew J. Ferris 1, 2 and Guifre Vidal 1, 3 1 The University of Queensland, School of Mathematics and Physics,
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 informationProblem set for the course Skálázás és renormálás a statisztikus fizikában, 2014
1 Problem set for the course Skálázás és renormálás a statisztikus fizikában, 014 Rules: You can choose at wish from problems having the same main number (i.e. from a given section), but you can collect
More informationMixing Quantum and Classical Mechanics: A Partially Miscible Solution
Mixing Quantum and Classical Mechanics: A Partially Miscible Solution R. Kapral S. Nielsen A. Sergi D. Mac Kernan G. Ciccotti quantum dynamics in a classical condensed phase environment how to simulate
More informationMethods of Fermion Monte Carlo
Methods of Fermion Monte Carlo Malvin H. Kalos kalos@llnl.gov Institute for Nuclear Theory University of Washington Seattle, USA July, 2013 LLNL-PRES-XXXXXX This work was performed under the auspices of
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 informationThe Use of Matrix Product State Representations in Simulating Many Body Quantum Systems
The Use of Matrix Product State Representations in Simulating Many Body Quantum Systems by Matthew Low A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF BACHELOR OF SCIENCE
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 informationTechniques for translationally invariant matrix product states
Techniques for translationally invariant matrix product states Ian McCulloch University of Queensland Centre for Engineered Quantum Systems (EQuS) 7 Dec 2017 Ian McCulloch (UQ) imps 7 Dec 2017 1 / 33 Outline
More informationThe Mott Metal-Insulator Transition
Florian Gebhard The Mott Metal-Insulator Transition Models and Methods With 38 Figures Springer 1. Metal Insulator Transitions 1 1.1 Classification of Metals and Insulators 2 1.1.1 Definition of Metal
More informationManybody wave function and Second quantization
Chap 1 Manybody wave function and Second quantization Ming-Che Chang Department of Physics, National Taiwan Normal University, Taipei, Taiwan (Dated: March 7, 2013) 1 I. MANYBODY WAVE FUNCTION A. One particle
More informationarxiv: v2 [cond-mat.str-el] 9 Jul 2012
Exact diagonalization of the one dimensional Bose-Hubbard model with local -body interactions Tomasz Sowiński Institute of Physics of the Polish Academy of Sciences, Aleja Lotników /46, -668 Warsaw, Poland
More informationBose-Hubbard Model (BHM) at Finite Temperature
Bose-Hubbard Model (BHM) at Finite Temperature - a Layman s (Sebastian Schmidt) proposal - pick up Diploma work at FU-Berlin with PD Dr. Axel Pelster (Uni Duisburg-Essen) ~ Diagrammatic techniques, high-order,
More informationarxiv: v2 [quant-ph] 31 Oct 2013
Quantum Criticality with the Multi-scale Entanglement Renormalization Ansatz arxiv:1109.5334v2 [quant-ph] 31 Oct 2013 G. Evenbly 1 and G. Vidal 2 1 The University of Queensland, Brisbane, Queensland 4072,
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 informationLinear Algebra using Dirac Notation: Pt. 2
Linear Algebra using Dirac Notation: Pt. 2 PHYS 476Q - Southern Illinois University February 6, 2018 PHYS 476Q - Southern Illinois University Linear Algebra using Dirac Notation: Pt. 2 February 6, 2018
More informationQuantum-Classical Hybrid Monte Carlo Algorithm with Applications to AQC
Quantum-Classical Hybrid Monte Carlo Algorithm with Applications to AQC Itay Hen Information Sciences Institute, USC Workshop on Theory and Practice of AQC and Quantum Simulation Trieste, Italy August
More informationHolographic Geometries from Tensor Network States
Holographic Geometries from Tensor Network States J. Molina-Vilaplana 1 1 Universidad Politécnica de Cartagena Perspectives on Quantum Many-Body Entanglement, Mainz, Sep 2013 1 Introduction & Motivation
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 information2D tensor network study of the S=1 bilinear-biquadratic Heisenberg model
2D tensor network study of the S=1 bilinear-biquadratic Heisenberg model Philippe Corboz, Institute for Theoretical Physics, University of Amsterdam AF phase Haldane phase 3-SL 120 phase? ipeps 2D tensor
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 informationInteraction-induced Symmetry Protected Topological Phase in Harper-Hofstadter models
Interaction-induced Symmetry Protected Topological Phase in Harper-Hofstadter models arxiv:1609.03760 Lode Pollet Dario Hügel Hugo Strand, Philipp Werner (Uni Fribourg) Algorithmic developments diagrammatic
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 informationMean-field theory for extended Bose-Hubbard model with hard-core bosons. Mathias May, Nicolas Gheeraert, Shai Chester, Sebastian Eggert, Axel Pelster
Mean-field theory for extended Bose-Hubbard model with hard-core bosons Mathias May, Nicolas Gheeraert, Shai Chester, Sebastian Eggert, Axel Pelster Outline Introduction Experiment Hamiltonian Periodic
More informationDECAY OF SINGLET CONVERSION PROBABILITY IN ONE DIMENSIONAL QUANTUM NETWORKS
DECAY OF SINGLET CONVERSION PROBABILITY IN ONE DIMENSIONAL QUANTUM NETWORKS SCOTT HOTTOVY Abstract. Quantum networks are used to transmit and process information by using the phenomena of quantum mechanics.
More information