Techniques for translationally invariant matrix product states
|
|
- Barbra McDaniel
- 5 years ago
- Views:
Transcription
1 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 / 33
2 Outline 1 Scaling relations in the thermodynamic limit Correlation funtions Transfer matrix Scaling dimensions 2 Expectation values of impo s Higher moments Energy variance Example = momentum distribution Binder cumulant Ian McCulloch (UQ) imps 7 Dec / 33
3 Correlation Functions The form of correlation functions are determined by the eigenvalues of the transfer operator Expansion in terms of eigenspectrum λ i : O(x)O(y) = i All eigenvalues 1 One eigenvalue equal to 1, corresponding to the identity operator a i λ y x i Ian McCulloch (UQ) imps 7 Dec / 33
4 1 Hubbard Model transfer matrix spectrum Half-filling, U/t = (0,0) Singlet (1,0) Spin triplet (0,1) Holon Triplet (1/2,1/2) Single-particle λ Number of states kept Ian McCulloch (UQ) imps 7 Dec / 33
5 Hubbard model transfer matrix spectrum Half-filling, U/t=4 100 (0,0) Singlet (1,0) Spin triplet (0,1) Holon Triplet (1/2,1/2) Single-particle Correlation length Number of states kept Ian McCulloch (UQ) imps 7 Dec / 33
6 Critical scaling example Two-species bose gas with linear tunneling Ω, from F. Zhan et al, Phys. Rev. A 90, ξ m Ω = Ω = Ω = Ω = Ω = Ω = Ian McCulloch (UQ) imps 7 Dec / 33
7 CFT Parameters For a critial mode, the correlation length increases with number of states m as a power law, ξ m κ [T. Nishino, K. Okunishi, M. Kikuchi, Phys. Lett. A 213, 69 (1996) M. Andersson, M. Boman, S. Östlund, Phys. Rev. B 59, (1999) L. Tagliacozzo, Thiago. R. de Oliveira, S. Iblisdir, J. I. Latorre, Phys. Rev. B 78, (2008)] This exponent is a function only of the central charge, κ = 6 12c + c [Pollmann et al, PRL 2009] Note: in practice this usually isn t a good way to determine c better to use entropy scaling Ian McCulloch (UQ) imps 7 Dec / 33
8 Scaling dimensions (Vid Stojevic, IPM et al, Phys. Rev. B 91, (2015)) Suppose we have a two-point correlator that has a power-law at large distances O(x)O(y) = y x 2 As we increase the number of states kept m the correlation length increases, so the region of validity of the power law increases. Take two different calculations with m 1 and m 2 Correlation lengths ξ 1 and ξ 2 ( ) ξ2 We expect: O(ξ 2) O(ξ 1 ) = ξ 1 for x large, we have: O(x) a λ x (with ξ = 1/ ln λ) Prefactor a is overlap of operator O with next-leading eigenvector of transfer operator a ξ This gives directly the operator scaling dimensions by direct fit Ian McCulloch (UQ) imps 7 Dec / 33
9 Scaling dimensions (Vid Stojevic, IPM et al, Phys. Rev. B 91, (2015)) Suppose we have a two-point correlator that has a power-law at large distances O(x)O(y) = y x 2 As we increase the number of states kept m the correlation length increases, so the region of validity of the power law increases. Take two different calculations with m 1 and m 2 Correlation lengths ξ 1 and ξ 2 ( ) ξ2 We expect: O(ξ 2) O(ξ 1 ) = ξ 1 for x large, we have: O(x) a λ x (with ξ = 1/ ln λ) Prefactor a is overlap of operator O with next-leading eigenvector of transfer operator a ξ This gives directly the operator scaling dimensions by direct fit Ian McCulloch (UQ) imps 7 Dec / 33
10 Scaling dimensions (Vid Stojevic, IPM et al, Phys. Rev. B 91, (2015)) Suppose we have a two-point correlator that has a power-law at large distances O(x)O(y) = y x 2 As we increase the number of states kept m the correlation length increases, so the region of validity of the power law increases. Take two different calculations with m 1 and m 2 Correlation lengths ξ 1 and ξ 2 ( ) ξ2 We expect: O(ξ 2) O(ξ 1 ) = ξ 1 for x large, we have: O(x) a λ x (with ξ = 1/ ln λ) Prefactor a is overlap of operator O with next-leading eigenvector of transfer operator a ξ This gives directly the operator scaling dimensions by direct fit Ian McCulloch (UQ) imps 7 Dec / 33
11 Heisenberg model fit for the scaling dimension prefactor of the spin operator at this mode idmrg data for m=15,20,25,30,35 y = * x^ transfer matrix eigenvalue 1 - λ = 1 / ξ Ian McCulloch (UQ) imps 7 Dec / 33
12 Expectation values of MPO s - arxiv: We have seen that we can write many interesting operators in the form of a matrix product operator Can we evaluate the expectation value of an arbitrary MPO? If the MPO has no Jordan structure, this is a simple eigenvalue problem W = λ For a triangular MPO, this doesn t work. But we can make use of the triangular structure I S z λs x index by index, each component is a function only of the previously calculated terms S z I Ian McCulloch (UQ) imps 7 Dec / 33
13 Expectation values of MPO s - arxiv: We have seen that we can write many interesting operators in the form of a matrix product operator Can we evaluate the expectation value of an arbitrary MPO? If the MPO has no Jordan structure, this is a simple eigenvalue problem W = λ For a triangular MPO, this doesn t work. But we can make use of the triangular structure I S z λs x index by index, each component is a function only of the previously calculated terms S z I Ian McCulloch (UQ) imps 7 Dec / 33
14 i W j Choose bond indices i, j of W ij, and denote T Wij Example: I Sz λs x S z I Eigentensor is (E 1 E 2 E 3 ) Starting from E 1 : E 1 = T I (E 1 ) = I is equivalent to the orthogonality condition - E 1 is just the identity E 2 : E 2 = T S z(e 1 ) = T S z(i) = S z E 3 : doesn t reach a fixed point, E 3 = E 3 (L) depends on the number of iterations L E 3 (L + 1) = T I (E 3 (L)) + T S z(e 2 ) + T λs x(e 1 ) = T I (E 3 (L)) + C where C is a constant matrix, C = T S z(s z ) + λs x Ian McCulloch (UQ) imps 7 Dec / 33
15 Fixed point equations for E 3 : Eigenmatrix expansion of T I : giving E (n) E 3 (L + 1) = T I (E 3 (L)) + C T I = λ n λ λ m 2 n=1 3 (L + 1) = λ ne (n) 3 (L) + C(n) Since λ 1 = 1 by construction, this motivates decomposing into components parallel (e 1 ) and perpendicular (E 1 ) to the identity: where Tr E 3 (L)ρ = 0 E 3 (L + 1) = E 3(L) + e 3 (L) I Ian McCulloch (UQ) imps 7 Dec / 33
16 Component in the direction of the identity: e 3 (L + 1) = e 3 (L) + Tr Cρ Has the solution e 3 (L) = L Tr Cρ is the energy Component perpendicular to the identity: E 3(L + 1) = T I (E 3(L)) + C where C = C (Tr Cρ) I E 3(L + 1) n = λ n E 3(L) n + C n Since all λ n < 1 here, this is a geometric series that converges to a fixed point (independent of L), (1 T I )(E 3) = C Linear solver for the unknown matrix E 3 Ian McCulloch (UQ) imps 7 Dec / 33
17 Component in the direction of the identity: e 3 (L + 1) = e 3 (L) + Tr Cρ Has the solution e 3 (L) = L Tr Cρ is the energy Component perpendicular to the identity: E 3(L + 1) = T I (E 3(L)) + C where C = C (Tr Cρ) I E 3(L + 1) n = λ n E 3(L) n + C n Since all λ n < 1 here, this is a geometric series that converges to a fixed point (independent of L), (1 T I )(E 3) = C Linear solver for the unknown matrix E 3 Ian McCulloch (UQ) imps 7 Dec / 33
18 Summary: Decompose eigentensor into components parallel and perpendicular to the identity The component parallel to the identity is the energy per site The perpendicular components reach a fixed point and give the Hamiltonian matrix elements E 1 = I Identity operator E 2 = Sz Sz block operator E 3 (L) = E + Le 3 Hamiltonian operator + energy per site Ian McCulloch (UQ) imps 7 Dec / 33
19 Generalization to arbitrary triangular MPO s arxiv: At the i th iteration, we have E i (L + 1) = T Wii (E i (L)) + T Wji (E j (L)) j<i }{{} = C(L) Basic idea: if W ii = 0, then E i = C if W ii 0, then solve (1 T Wii )(E i ) = C The result will be a polynomial function of L solve separately for the coefficient of the k-th power of n If the diagonal element is unitary, then obtain the eigenvalues of magnitude 1 If any eigenvalues are complex, then expand also in fourier modes Ian McCulloch (UQ) imps 7 Dec / 33
20 Generalization to arbitrary triangular MPO s arxiv: At the i th iteration, we have E i (L + 1) = T Wii (E i (L)) + T Wji (E j (L)) j<i }{{} = C(L) Basic idea: if W ii = 0, then E i = C if W ii 0, then solve (1 T Wii )(E i ) = C The result will be a polynomial function of L solve separately for the coefficient of the k-th power of n If the diagonal element is unitary, then obtain the eigenvalues of magnitude 1 If any eigenvalues are complex, then expand also in fourier modes Ian McCulloch (UQ) imps 7 Dec / 33
21 Example 1 - Variance How close is a variational state to an eigenstate of the Hamiltonian? Sometimes there is an algorithmic measure, often not. The square of the Hamiltonian operator determines the energy variance H 2 L H 2 L = (H E) 2 L = Lσ 2 A universal measure for the quality of a variational wavefunction lower bound for the energy: there is always an eigenstate within σ of E We can easily construct an MPO representation of H 2 Ian McCulloch (UQ) imps 7 Dec / 33
22 Spin 1/2 Heisenberg Model Energy per site scaling with variance (exact energy = -ln = ) E σ 2 Ian McCulloch (UQ) imps 7 Dec / 33
23 Example 2 - momentum distribution Momentum-dependent operators have a simple form, b k = x e ikx b x W b k ( e = ik I b I ) Momentum occupation: n(k) = 1 L b k b k Broken U(1) symmetry: b 0 hence n(k = 0) L (extensive) With U(1) symmetry: n(k = 0) is finite, but diverges with m in superfluid phase Ian McCulloch (UQ) imps 7 Dec / 33
24 1000 Bose-Hubbard Model N(k) Infinite 1D, one particle per site 100 U/J=3.2 U/J=3.6 U/J=4.0 U/J=4.8 U/J=6.4 N(k) k/pi Ian McCulloch (UQ) imps 7 Dec / 33
25 Higher moments It is straight forward to evaluate a local order parameter, eg M = i M i The first moment of this operator gives the order parameter, M = m 1 (L) It is also useful to calculate higher moments, eg M 2 = m 2 (L) or generally M k = m k (L) These are polynomial functions in the system size L. Ian McCulloch (UQ) imps 7 Dec / 33
26 For finite systems, the Binder cumulant of the order parameter cancels the leading-order finite size effects NL/ U L = 1 m4 3 m 2 2 L = 50 L = 100 L = 150 L = 200 idmrg finite size scaling Ω 0.8 L = 50 L = L = 150 L = 200 UL Ian McCulloch (UQ) imps 7 Dec / 33
27 Cumulant expansions Express the moments m i in terms of the cumulants per site κ j, m 1 (L) = κ 1 L m 2 (L) = κ 2 1 L2 + κ 2 L m 3 (L) = κ 3 1 L3 + 3κ 1 κ 2 L 2 + κ 3 L m 4 (L) = κ 4 1 L4 + 6κ 2 1 κ 2L 3 + (3κ κ 1κ 3 )L 2 + κ 4 L κ 1 is the order parameter itself κ 2 is the variance (related to the susceptibility) κ 3 is the skewness κ 4 is the kurtosis The cumulants per site κ k are well-defined for an imps Note: the cumulants are normally written such that they are extensive quantities Lκ k. Ian McCulloch (UQ) imps 7 Dec / 33
28 imps for two-component bose gas The cumulant expansion already gives a lot of information κ 1 is the order parameter itself Ising transition in 2-component Bose gas Order parameter N_a - N_b κ Ω Ian McCulloch (UQ) imps 7 Dec / 33
29 The second cumulant gives the susceptibility 2-component Bose gas Second cumulant (susceptibility) κ Ω Different to a finite-size scaling, the susceptibility exactly diverges at the critical point. Sufficiently close to the critical point, it looks mean-field-like (so will generally give the wrong exponent!) Ian McCulloch (UQ) imps 7 Dec / 33
30 The fourth cumulant changes sign at the transition. 2-component Bose gas fourth cumulant 0-2e+05 κ 4-4e+05-6e+05-8e Ω Ian McCulloch (UQ) imps 7 Dec / 33
31 Binder Cumulant for imps Naively taking the limit L for the Binder cumulant doesn t produce anything useful: if the order parameter κ 1 0, U L = 1 m4 L 3 m 2 2 L 2 3 if κ 1 = 0, then m 4 (L) = 3k 2 2 L2 + k 4 L Hence U L = 1 3k2 2 L2 + k 4 L 3k 2 2 L2 0 Finally, a step function that detects whether the order parameter is non-zero Better approach, in the spirit of finite-entanglement scaling: Evaluate the moment polynomial using L correlation length Ian McCulloch (UQ) imps 7 Dec / 33
32 Transverse field Ising model Binder cumulant, scale factor s=5 M Order parameter m=4 m=8 m= λ U m m=4 m=5 m=6 m=7 m=8 m=9 m=10 m=11 m= λ Ian McCulloch (UQ) imps 7 Dec / 33
33 String parameters Order parameters do not have to be local Mott insulator string order parameter O 2 P = lim j i Πj k=i ( 1)nk We can write this as a correlation function of kink operators, p i = Π k<i ( 1) nk This turns the string order into a 2-point correlation function: O 2 P = lim p i p j j i Or as an order parameter: P = i p i Then O 2 p = 1 L 2 P 2 Ian McCulloch (UQ) imps 7 Dec / 33
34 Real example: 3-leg Bose-Hubbard model with flux phase (F. Kolley, M. Piraud, IPM, U. Schollwoeck, F. Heidrich-Meisner, New J. Phys. 17 (2015) ) For density n = 1/3 (one particle per rung), near flux φ π, there is a transition from a Mott to critical as a function of J P has a simple MPO representation ( ) e iπn I P = I Hence we can calculate higher moments of P. Ian McCulloch (UQ) imps 7 Dec / 33
35 Bose-Hubbard Ladder String order parameter O p m=100 m=150 m=200 m=300 m=400 m=500 m= Jperp Ian McCulloch (UQ) imps 7 Dec / 33
36 ξ Jperp = 0.99 Jperp = 1.00 Jperp = 1.01 Jperp = 1.02 Jperp = 1.03 Jperp = 1.04 Jperp = 1.05 Jperp = 1.06 Jperp = 1.10 Jperp = 1.14 Jperp = 1.25 Bose=Hubbard Ladder Scaling of correlation length m Ian McCulloch (UQ) imps 7 Dec / 33
37 Bose-Hubbard Ladder String parameter Binder cumulant U ξ m=200 m=250 m=300 m=350 m=400 m=450 m=500 m=550 m= Jperp Ian McCulloch (UQ) imps 7 Dec / 33
38 Summary Download the code: finite entanglement scaling higher moments Binder cumulant Ian McCulloch (UQ) imps 7 Dec / 33
Infinite DMRG and the Matrix Product Toolkit
Infinite DMRG and the Matrix Product Toolkit Ian McCulloch University of Queensland Centre for Engineered Quantum Systems (EQuS) 21/9/2015 Ian McCulloch (UQ) imps 21/9/2015 1 / 44 Outline 1 imps 30 second
More informationPhases of strongly-interacting bosons on a two-leg ladder
Phases of strongly-interacting bosons on a two-leg ladder Marie Piraud Arnold Sommerfeld Center for Theoretical Physics, LMU, Munich April 20, 2015 M. Piraud Phases of strongly-interacting bosons on a
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 informationSpin liquids on ladders and in 2d
Spin liquids on ladders and in 2d MPA Fisher (with O. Motrunich) Minnesota, FTPI, 5/3/08 Interest: Quantum Spin liquid phases of 2d Mott insulators Background: Three classes of 2d Spin liquids a) Topological
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 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 informationGapless Spin Liquids in Two Dimensions
Gapless Spin Liquids in Two Dimensions MPA Fisher (with O. Motrunich, Donna Sheng, Matt Block) Boulder Summerschool 7/20/10 Interest Quantum Phases of 2d electrons (spins) with emergent rather than broken
More informationMatrix Product Operators: Algebras and Applications
Matrix Product Operators: Algebras and Applications Frank Verstraete Ghent University and University of Vienna Nick Bultinck, Jutho Haegeman, Michael Marien Burak Sahinoglu, Dominic Williamson Ignacio
More informationTime Evolving Block Decimation Algorithm
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,
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 informationChiral spin liquids. Bela Bauer
Chiral spin liquids Bela Bauer Based on work with: Lukasz Cinco & Guifre Vidal (Perimeter Institute) Andreas Ludwig & Brendan Keller (UCSB) Simon Trebst (U Cologne) Michele Dolfi (ETH Zurich) Nature Communications
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 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 informationGround State Projector QMC in the valence-bond basis
Quantum Monte Carlo Methods at Work for Novel Phases of Matter Trieste, Italy, Jan 23 - Feb 3, 2012 Ground State Projector QMC in the valence-bond basis Anders. Sandvik, Boston University Outline: The
More informationNews on tensor network algorithms
News on tensor network algorithms Román Orús Donostia International Physics Center (DIPC) December 6th 2018 S. S. Jahromi, RO, M. Kargarian, A. Langari, PRB 97, 115162 (2018) S. S. Jahromi, RO, PRB 98,
More informationNon-magnetic states. The Néel states are product states; φ N a. , E ij = 3J ij /4 2 The Néel states have higher energy (expectations; not eigenstates)
Non-magnetic states Two spins, i and j, in isolation, H ij = J ijsi S j = J ij [Si z Sj z + 1 2 (S+ i S j + S i S+ j )] For Jij>0 the ground state is the singlet; φ s ij = i j i j, E ij = 3J ij /4 2 The
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 informationMatrix Product States
Matrix Product States Ian McCulloch University of Queensland Centre for Engineered Quantum Systems 28 August 2017 Hilbert space (Hilbert) space is big. Really big. You just won t believe how vastly, hugely,
More informationG : Quantum Mechanics II
G5.666: Quantum Mechanics II Notes for Lecture 7 I. A SIMPLE EXAMPLE OF ANGULAR MOMENTUM ADDITION Given two spin-/ angular momenta, S and S, we define S S S The problem is to find the eigenstates of the
More informationThe Quantum Heisenberg Ferromagnet
The Quantum Heisenberg Ferromagnet Soon after Schrödinger discovered the wave equation of quantum mechanics, Heisenberg and Dirac developed the first successful quantum theory of ferromagnetism W. Heisenberg,
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 informationQuantum Spin-Metals in Weak Mott Insulators
Quantum Spin-Metals in Weak Mott Insulators MPA Fisher (with O. Motrunich, Donna Sheng, Simon Trebst) Quantum Critical Phenomena conference Toronto 9/27/08 Quantum Spin-metals - spin liquids with Bose
More informationCLASSIFICATION OF SU(2)-SYMMETRIC TENSOR NETWORKS
LSSIFITION OF SU(2)-SYMMETRI TENSOR NETWORKS Matthieu Mambrini Laboratoire de Physique Théorique NRS & Université de Toulouse M.M., R. Orús, D. Poilblanc, Phys. Rev. 94, 2524 (26) D. Poilblanc, M.M., arxiv:72.595
More informationPhysics 127b: Statistical Mechanics. Renormalization Group: 1d Ising Model. Perturbation expansion
Physics 17b: Statistical Mechanics Renormalization Group: 1d Ising Model The ReNormalization Group (RNG) gives an understanding of scaling and universality, and provides various approximation schemes to
More informationDomains and Domain Walls in Quantum Spin Chains
Domains and Domain Walls in Quantum Spin Chains Statistical Interaction and Thermodynamics Ping Lu, Jared Vanasse, Christopher Piecuch Michael Karbach and Gerhard Müller 6/3/04 [qpis /5] Domains and Domain
More information3.15. Some symmetry properties of the Berry curvature and the Chern number.
50 Phys620.nb z M 3 at the K point z M 3 3 t ' sin 3 t ' sin (3.36) (3.362) Therefore, as long as M 3 3 t ' sin, the system is an topological insulator ( z flips sign). If M 3 3 t ' sin, z is always positive
More information=
.52.5.48.46.44.42 =..22 2 x 3 4 Application of DMRG to 2D systems L y L x Unit cell Periodic boundary conditions for both x and y directions k y = 2n/L y = X n / l 2 Initial basis states (Landau gauge)
More informationIntroduction to Quantum Monte Carlo
Entanglement in Strongly Correlated Systems @ Benasque Feb. 6-17, 2017 Introduction to Quantum Monte Carlo Naoki KAWASHIMA (ISSP) 2017.02.06-07 Why bother? Estimating scaling dimension by TRG, TNR, etc
More informationMatrix product states for the fractional quantum Hall effect
Matrix product states for the fractional quantum Hall effect Roger Mong (California Institute of Technology) University of Virginia Feb 24, 2014 Collaborators Michael Zaletel UC Berkeley (Stanford/Station
More informationMomentum-space and Hybrid Real- Momentum Space DMRG applied to the Hubbard Model
Momentum-space and Hybrid Real- Momentum Space DMRG applied to the Hubbard Model Örs Legeza Reinhard M. Noack Collaborators Georg Ehlers Jeno Sólyom Gergely Barcza Steven R. White Collaborators Georg Ehlers
More informationVarious Facets of Chalker- Coddington network model
Various Facets of Chalker- Coddington network model V. Kagalovsky Sami Shamoon College of Engineering Beer-Sheva Israel Context Integer quantum Hall effect Semiclassical picture Chalker-Coddington Coddington
More informationSolving the Schrödinger equation for the Sherrington Kirkpatrick model in a transverse field
J. Phys. A: Math. Gen. 30 (1997) L41 L47. Printed in the UK PII: S0305-4470(97)79383-1 LETTER TO THE EDITOR Solving the Schrödinger equation for the Sherrington Kirkpatrick model in a transverse field
More information= a. a = Let us now study what is a? c ( a A a )
7636S ADVANCED QUANTUM MECHANICS Solutions 1 Spring 010 1 Warm up a Show that the eigenvalues of a Hermitian operator A are real and that the eigenkets of A corresponding to dierent eigenvalues are orthogonal
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 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 informationQuantum quenches in 2D with chain array matrix product states
Quantum quenches in 2D with chain array matrix product states Andrew J. A. James University College London Robert M. Konik Brookhaven National Laboratory arxiv:1504.00237 Outline MPS for many body systems
More informationPHY305: Notes on Entanglement and the Density Matrix
PHY305: Notes on Entanglement and the Density Matrix Here follows a short summary of the definitions of qubits, EPR states, entanglement, the density matrix, pure states, mixed states, measurement, and
More 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 informationQuantum Mechanics C (130C) Winter 2014 Final exam
University of California at San Diego Department of Physics Prof. John McGreevy Quantum Mechanics C (130C Winter 014 Final exam Please remember to put your name on your exam booklet. This is a closed-book
More informationLinked-Cluster Expansions for Quantum Many-Body Systems
Linked-Cluster Expansions for Quantum Many-Body Systems Boulder Summer School 2010 Simon Trebst Lecture overview Why series expansions? Linked-cluster expansions From Taylor expansions to linked-cluster
More informationQuantum Monte Carlo Simulations in the Valence Bond Basis
NUMERICAL APPROACHES TO QUANTUM MANY-BODY SYSTEMS, IPAM, January 29, 2009 Quantum Monte Carlo Simulations in the Valence Bond Basis Anders W. Sandvik, Boston University Collaborators Kevin Beach (U. of
More information(Dynamical) quantum typicality: What is it and what are its physical and computational implications?
(Dynamical) : What is it and what are its physical and computational implications? Jochen Gemmer University of Osnabrück, Kassel, May 13th, 214 Outline Thermal relaxation in closed quantum systems? Typicality
More informationSolving the sign problem for a class of frustrated antiferromagnets
Solving the sign problem for a class of frustrated antiferromagnets Fabien Alet Laboratoire de Physique Théorique Toulouse with : Kedar Damle (TIFR Mumbai), Sumiran Pujari (Toulouse Kentucky TIFR Mumbai)
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 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 informationQuantum metrology from a quantum information science perspective
1 / 41 Quantum metrology from a quantum information science perspective Géza Tóth 1 Theoretical Physics, University of the Basque Country UPV/EHU, Bilbao, Spain 2 IKERBASQUE, Basque Foundation for Science,
More informationA new perspective on long range SU(N) spin models
A new perspective on long range SU(N) spin models Thomas Quella University of Cologne Workshop on Lie Theory and Mathematical Physics Centre de Recherches Mathématiques (CRM), Montreal Based on work with
More informationNumerical diagonalization studies of quantum spin chains
PY 502, Computational Physics, Fall 2016 Anders W. Sandvik, Boston University Numerical diagonalization studies of quantum spin chains Introduction to computational studies of spin chains Using basis states
More informationSpin liquids in frustrated magnets
May 20, 2010 Contents 1 Frustration 2 3 4 Exotic excitations 5 Frustration The presence of competing forces that cannot be simultaneously satisfied. Heisenberg-Hamiltonian H = 1 J ij S i S j 2 ij The ground
More informationLecture 4 Quantum mechanics in more than one-dimension
Lecture 4 Quantum mechanics in more than one-dimension Background Previously, we have addressed quantum mechanics of 1d systems and explored bound and unbound (scattering) states. Although general concepts
More informationHandout 10. Applications to Solids
ME346A Introduction to Statistical Mechanics Wei Cai Stanford University Win 2011 Handout 10. Applications to Solids February 23, 2011 Contents 1 Average kinetic and potential energy 2 2 Virial theorem
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 informationFrustration and Area law
Frustration and Area law When the frustration goes odd S. M. Giampaolo Institut Ruder Bošković, Zagreb, Croatia Workshop: Exactly Solvable Quantum Chains Natal 18-29 June 2018 Coauthors F. Franchini Institut
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 informationQuantum phase transitions and entanglement in (quasi)1d spin and electron models
Quantum phase transitions and entanglement in (quasi)1d spin and electron models Elisa Ercolessi - Università di Bologna Group in Bologna: G.Morandi, F.Ortolani, E.E., C.Degli Esposti Boschi, A.Anfossi
More informationMagnetic ordering of local moments
Magnetic ordering Types of magnetic structure Ground state of the Heisenberg ferromagnet and antiferromagnet Spin wave High temperature susceptibility Mean field theory Magnetic ordering of local moments
More informationÖsszefonódás és felületi törvény 2. Szabad rácsmodellek
Összefonódás és felületi törvény 2. Szabad rácsmodellek Eisler Viktor MTA-ELTE Elméleti Fizikai Kutatócsoport Entanglement Day 2014.09.05. I. Peschel & V. Eisler, J. Phys. A: Math. Theor. 42, 504003 (2009)
More informationQuantum gases in the unitary limit and...
Quantum gases in the unitary limit and... Andre LeClair Cornell university Benasque July 2 2010 Outline The unitary limit of quantum gases S-matrix based approach to thermodynamics Application to the unitary
More informationChapter 29. Quantum Chaos
Chapter 29 Quantum Chaos What happens to a Hamiltonian system that for classical mechanics is chaotic when we include a nonzero h? There is no problem in principle to answering this question: given a classical
More informationMeasuring entanglement in synthetic quantum systems
Measuring entanglement in synthetic quantum systems ψ?? ψ K. Rajibul Islam Institute for Quantum Computing and Department of Physics and Astronomy University of Waterloo research.iqc.uwaterloo.ca/qiti/
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 informationLecture 4 Quantum mechanics in more than one-dimension
Lecture 4 Quantum mechanics in more than one-dimension Background Previously, we have addressed quantum mechanics of 1d systems and explored bound and unbound (scattering) states. Although general concepts
More informationScaling Theory. Roger Herrigel Advisor: Helmut Katzgraber
Scaling Theory Roger Herrigel Advisor: Helmut Katzgraber 7.4.2007 Outline The scaling hypothesis Critical exponents The scaling hypothesis Derivation of the scaling relations Heuristic explanation Kadanoff
More informationSPIN-LIQUIDS ON THE KAGOME LATTICE: CHIRAL TOPOLOGICAL, AND GAPLESS NON-FERMI-LIQUID PHASE
SPIN-LIQUIDS ON THE KAGOME LATTICE: CHIRAL TOPOLOGICAL, AND GAPLESS NON-FERMI-LIQUID PHASE ANDREAS W.W. LUDWIG (UC-Santa Barbara) work done in collaboration with: Bela Bauer (Microsoft Station-Q, Santa
More informationQuantum Cluster Methods (CPT/CDMFT)
Quantum Cluster Methods (CPT/CDMFT) David Sénéchal Département de physique Université de Sherbrooke Sherbrooke (Québec) Canada Autumn School on Correlated Electrons Forschungszentrum Jülich, Sept. 24,
More informationTime-Independent Perturbation Theory
4 Phys46.nb Time-Independent Perturbation Theory.. Overview... General question Assuming that we have a Hamiltonian, H = H + λ H (.) where λ is a very small real number. The eigenstates of the Hamiltonian
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 informationarxiv:quant-ph/ v2 24 Dec 2003
Quantum Entanglement in Heisenberg Antiferromagnets V. Subrahmanyam Department of Physics, Indian Institute of Technology, Kanpur, India. arxiv:quant-ph/0309004 v2 24 Dec 2003 Entanglement sharing among
More informationDecoherence and Thermalization of Quantum Spin Systems
Copyright 2011 American Scientific Publishers All rights reserved Printed in the United States of America Journal of Computational and Theoretical Nanoscience Vol. 8, 1 23, 2011 Decoherence and Thermalization
More informationFully symmetric and non-fractionalized Mott insulators at fractional site-filling
Fully symmetric and non-fractionalized Mott insulators at fractional site-filling Itamar Kimchi University of California, Berkeley EQPCM @ ISSP June 19, 2013 PRL 2013 (kagome), 1207.0498...[PNAS] (honeycomb)
More informationQuantum Monte Carlo Simulations in the Valence Bond Basis. Anders Sandvik, Boston University
Quantum Monte Carlo Simulations in the Valence Bond Basis Anders Sandvik, Boston University Outline The valence bond basis for S=1/2 spins Projector QMC in the valence bond basis Heisenberg model with
More information8.334: Statistical Mechanics II Spring 2014 Test 2 Review Problems
8.334: Statistical Mechanics II Spring 014 Test Review Problems The test is closed book, but if you wish you may bring a one-sided sheet of formulas. The intent of this sheet is as a reminder of important
More informationEntanglement spectrum as a tool for onedimensional
Entanglement spectrum as a tool for onedimensional critical systems MPI-PKS, Dresden, November 2012 Joel Moore University of California, Berkeley, and Lawrence Berkeley National Laboratory Outline ballistic
More informationC/CS/Phy191 Problem Set 6 Solutions 3/23/05
C/CS/Phy191 Problem Set 6 Solutions 3/3/05 1. Using the standard basis (i.e. 0 and 1, eigenstates of Ŝ z, calculate the eigenvalues and eigenvectors associated with measuring the component of spin along
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 informationSolutions Final exam 633
Solutions Final exam 633 S.J. van Enk (Dated: June 9, 2008) (1) [25 points] You have a source that produces pairs of spin-1/2 particles. With probability p they are in the singlet state, ( )/ 2, and with
More informationquasi-particle pictures from continuous unitary transformations
quasi-particle pictures from continuous unitary transformations Kai Phillip Schmidt 24.02.2016 quasi-particle pictures from continuous unitary transformations overview Entanglement in Strongly Correlated
More informationEntanglement, holography, and strange metals
Entanglement, holography, and strange metals PCTS, Princeton, October 26, 2012 Subir Sachdev Talk online at sachdev.physics.harvard.edu HARVARD Liza Huijse Max Metlitski Brian Swingle Complex entangled
More informationStochastic series expansion (SSE) and ground-state projection
Institute of Physics, Chinese Academy of Sciences, Beijing, October 31, 2014 Stochastic series expansion (SSE) and ground-state projection Anders W Sandvik, Boston University Review article on quantum
More informationEntanglement spectrum and Matrix Product States
Entanglement spectrum and Matrix Product States Frank Verstraete J. Haegeman, D. Draxler, B. Pirvu, V. Stojevic, V. Zauner, I. Pizorn I. Cirac (MPQ), T. Osborne (Hannover), N. Schuch (Aachen) Outline Valence
More informationLattice modulation experiments with fermions in optical lattices and more
Lattice modulation experiments with fermions in optical lattices and more Nonequilibrium dynamics of Hubbard model Ehud Altman Weizmann Institute David Pekker Harvard University Rajdeep Sensarma Harvard
More informationLecture 6. Four postulates of quantum mechanics. The eigenvalue equation. Momentum and energy operators. Dirac delta function. Expectation values
Lecture 6 Four postulates of quantum mechanics The eigenvalue equation Momentum and energy operators Dirac delta function Expectation values Objectives Learn about eigenvalue equations and operators. Learn
More informationGiant Enhancement of Quantum Decoherence by Frustrated Environments
ISSN 0021-3640, JETP Letters, 2006, Vol. 84, No. 2, pp. 99 103. Pleiades Publishing, Inc., 2006.. Giant Enhancement of Quantum Decoherence by Frustrated Environments S. Yuan a, M. I. Katsnelson b, and
More informationUniversal Quantum Simulator, Local Convertibility and Edge States in Many-Body Systems Fabio Franchini
New Frontiers in Theoretical Physics XXXIV Convegno Nazionale di Fisica Teorica - Cortona Universal Quantum Simulator, Local Convertibility and Edge States in Many-Body Systems Fabio Franchini Collaborators:
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 informationPartition functions for complex fugacity
Partition functions for complex fugacity Part I Barry M. McCoy CN Yang Institute of Theoretical Physics State University of New York, Stony Brook, NY, USA Partition functions for complex fugacity p.1/51
More informationSupersymmetry breaking and Nambu-Goldstone fermions in lattice models
YKIS2016@YITP (2016/6/15) Supersymmetry breaking and Nambu-Goldstone fermions in lattice models Hosho Katsura (Department of Physics, UTokyo) Collaborators: Yu Nakayama (IPMU Rikkyo) Noriaki Sannomiya
More informationGeneral Exam Part II, Fall 1998 Quantum Mechanics Solutions
General Exam Part II, Fall 1998 Quantum Mechanics Solutions Leo C. Stein Problem 1 Consider a particle of charge q and mass m confined to the x-y plane and subject to a harmonic oscillator potential V
More informationNegative Binomial Coherent States
Negative Binomial Coherent States Superstatistics of Quanta T.S. Biró 1 1 Heavy Ion Research Group Hung.Acad.Sci. Research Centre for Physics, Budapest October 6, 2015 Talk given by T. S. Biró October
More informationGlobal phase diagrams of two-dimensional quantum antiferromagnets. Subir Sachdev Harvard University
Global phase diagrams of two-dimensional quantum antiferromagnets Cenke Xu Yang Qi Subir Sachdev Harvard University Outline 1. Review of experiments Phases of the S=1/2 antiferromagnet on the anisotropic
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 informationSimulations of Quantum Dimer Models
Simulations of Quantum Dimer Models Didier Poilblanc Laboratoire de Physique Théorique CNRS & Université de Toulouse 1 A wide range of applications Disordered frustrated quantum magnets Correlated fermions
More informationLecture 8 Nature of ensemble: Role of symmetry, interactions and other system conditions: Part II
Lecture 8 Nature of ensemble: Role of symmetry, interactions and other system conditions: Part II We continue our discussion of symmetries and their role in matrix representation in this lecture. An example
More informationThe nature of superfluidity in the cold atomic unitary Fermi gas
The nature of superfluidity in the cold atomic unitary Fermi gas Introduction Yoram Alhassid (Yale University) Finite-temperature auxiliary-field Monte Carlo (AFMC) method The trapped unitary Fermi gas
More informationKern- und Teilchenphysik I Lecture 13:Quarks and QCD
Kern- und Teilchenphysik I Lecture 13:Quarks and QCD (adapted from the Handout of Prof. Mark Thomson) Prof. Nico Serra Dr. Patrick Owen, Dr. Silva Coutinho http://www.physik.uzh.ch/de/lehre/phy211/hs2016.html
More informationIntroduction to string theory 2 - Quantization
Remigiusz Durka Institute of Theoretical Physics Wroclaw / 34 Table of content Introduction to Quantization Classical String Quantum String 2 / 34 Classical Theory In the classical mechanics one has dynamical
More informationExcitonic Condensation in Systems of Strongly Correlated Electrons. Jan Kuneš and Pavel Augustinský DFG FOR1346
Excitonic Condensation in Systems of Strongly Correlated Electrons Jan Kuneš and Pavel Augustinský DFG FOR1346 Motivation - unconventional long-range order incommensurate spin spirals complex order parameters
More informationParticle Physics. Michaelmas Term 2009 Prof Mark Thomson. Handout 7 : Symmetries and the Quark Model. Introduction/Aims
Particle Physics Michaelmas Term 2009 Prof Mark Thomson Handout 7 : Symmetries and the Quark Model Prof. M.A. Thomson Michaelmas 2009 205 Introduction/Aims Symmetries play a central role in particle physics;
More informationPHYS-454 The position and momentum representations
PHYS-454 The position and momentum representations 1 Τhe continuous spectrum-a n So far we have seen problems where the involved operators have a discrete spectrum of eigenfunctions and eigenvalues.! n
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 information