Quantum Cluster Methods
|
|
- Janel Golden
- 5 years ago
- Views:
Transcription
1 Quantum Cluster Methods An introduction David Sénéchal Université de Sherbrooke Département de physique CIFAR - PITP International Summer School on Numerical Methods for Correlated Systems in Condensed Matter May 2008
2 Outline Exact Diagonalizations Clusters and Cluster Perturbation Theory (CPT) The Self-Energy Functional Approach The Variational Cluster Approximation (VCA) Cluster Dynamical Mean Field Theory (CDMFT)
3 Part I Exact Diagonalizations
4 An old Persian Legend, revisited
5 The Hubbard Model on finite cluster Simple Hubbard model (conserves N and N separately): H = a,b,σ t ab c aσc bσ + U a n a n a µ a n a Typical cluster (L sites): e 2 e 1
6 Hamiltonian matrix : 2 sites Half-filled Hubbard model L = 2 U 2µ t t 0 t 2µ 0 t t 0 2µ t 0 t t U 2µ
7 Hamiltonian matrix : 6 sites Sparse matrix structure
8 Hamiltonian matrix : example Dimension of the Hilbert space (half-filled Hubbard model): ( d = L! [(L/2)!] 2 ) 2 2 4L πl One double-precision vector means 1.23 GB of memory L dimension
9 Steps 1. Building a basis 2. Constructing the Hamiltonian matrix 3. Finding the ground state (e.g. by the Lanczos method) 4. Calculating the one-body Green function
10 Coding of the states Basis of occupation number eigenstates: (c 1 )n 1 (c L )n L (c 1 )n 1 (c L )n L 0 Binary representation of basis states: b where b = b + 2 L b Example : b = ( ) = = 349, 866 Need a direct table:... and a reverse table: consecutive label b = B (i ) b = B (i ) mod i = I (b ) + d N I (b ) i = i% d N i = i/ int. division d N
11 Coding of the states (2) Tensor product structure of the Hilbert space: V = V N V N dimension: Example (6 sites): d = d(n )d(n ) d(n σ ) = L! N σ!(l N σ )!
12 Constructing the Hamiltonian matrix Form of Hamiltonian: H = K K + V int. K = a,b t ab c ac b K is stored in sparse form. V int. is diagonal and is stored. Matrix elements of V int. : bit count(b & b ) Two basis states b σ and b σ are connected with the matrix K if their binary representations differ at two positions a and b. b K b = ( 1) M ab t ab M ab = b 1 c=a+1 n c
13 The Lanczos algorithm Finds the lowest eigenpair by an iterative application of H Start with random vector φ 0 An iterative procedure builds the Krylov subspace: K = span { φ 0, H φ 0, H 2 φ 0,, H M φ 0 } Lanczos three-way recursion: φ n+1 = H φ n a n φ n b 2 n φ n 1 a n = φ n H φ n φ n φ n b 2 n = φ n φ n φ n 1 φ n 1 b 0 = 0
14 The Lanczos algorithm (2) In the basis of normalized states n = φ n / φ n φ n, the projected Hamiltonian has the tridiagonal form a 0 b b 1 a 1 b H = 0 b 2 a 2 b a N At each step n, find the lowest eigenvalue of that matrix Stop when the lowest eigenvalue E 0 has converged ( E 0 /E 0 < ) Then re-run to find eigenvector ψ = n ψ n n as the φ n s are not kept in memory.
15 Lanczos method: characteristics Typical required number of iterations: from 20 to 200 Extreme eigenvalues converge first Rate of convergence increases with separation between ground state and first excited state Cannot resolve degenerate ground states : only one state per ground state manifold is picked up If one is interested in low lying states, periodic re-orthogonalization may be required, as orthogonality leaks will occur For degenerate ground states and low lying states (e.g. in DMRG), the Davidson method is generally preferable
16 Lanczos method: illustration of the convergence 100 iterations on a matrix of dimension 600: eigenvalues of the tridiagonal projection as a function of iteration step
17 Lanczos method for the Green function Zero temperature Green function: G µν (ω) = G µν,e (ω) + G µν,h (ω) 1 G µν,e (ω) = Ω c µ c ω H + E ν Ω 0 G µν,h (ω) = Ω c 1 ν c µ Ω ω + H E 0 Consider the diagonal element Use the expansion φ µ = c 1 µ Ω = G µµ,e = φ µ φ µ ω H + E 0 1 z H = 1 z + 1 z 2 H + 1 z 3 H2 +
18 Lanczos method for the Green function (2) Truncated expansion evaluated exactly in Krylov subspace generated by φ µ if we perform a Lanczos procedure on φ µ. Then G µµ,e is given by a Jacobi continued fraction: G µµ,e (ω) = ω a 0 φ µ φ µ ω a 1 b 2 1 b 2 2 ω a 2 The coefficients a n and b n are stored in memory What about non diagonal elements G µν,e? See, e.g., E. Dagotto, Rev. Mod. Phys. 66:763 (1994)
19 Lanczos method for the Green function (3) Trick: Define the combination G + 1 µν,e(ω) = Ω (c µ + c ν ) (c µ + c ν ) Ω ω H + E 0 G + µν,e(ω) can be calculated like G µµ,e (ω) Since G µν,e (ω) = G νµ,e (ω), then G µν,e (ω) = 1 [ G + 2 µν,e (ω) G µµ,e (ω) G νν,e (ω) ] Likewise for G µν,h (ω)
20 Lehman representation Lehmann representation of the Green function G µν (ω) = m + n 1 Ω c µ m m c ω E m + E ν Ω 0 Ω c 1 ν n n c µ Ω ω + E n E 0 Define the matices Then Q (e) µm = Ω c µ m Q (h) µn = Ω c µ n G µν (ω) = m = r Q (e) µmq (e) νm ω ω m (e) Q µr Q νr ω ω r + n Q (h) µn Q (h) νn ω ω n (h)
21 Alternate way : The Band Lanczos method Define φ µ = c µ Ω, µ = 1,..., L. Extended Krylov space : { φ 1,..., φ L,H φ 1,..., H φ L,..., } (H) M φ 1,..., (H) M φ L States are built iteratively and orthogonalized Possible linearly dependent states are eliminated ( deflation ) A band representation of the Hamiltonian (2L + 1 diagonals) is formed in the Krylov subspace. It is diagonalized and the eigenpairs are used to build an approximate Lehmann representation dongarra/etemplates/node131.html
22 Lanczos vs Band Lanczos The usual Lanczos method for the Green function needs 3 vectors in memory, and L(L + 1) Lanczos procedures. The Band Lanczos method requires 3L + 1 vectors in memory, but requires only 2 iterative procedures ((e) et (h)). If Memory allows it, the band Lanczos is much faster.
23 Cluster symmetries Clusters with C 2v symmetry Clusters with C 2 symmetry
24 Cluster symmetries (2) Symmetry operations form a group G The most common occurences are : C1 : The trivial group (no symmetry) C2 : The 2-element group (e.g. left-right symmetry) C2v : 2 reflections, 1 π-rotation C 4v : 4 reflections, 1 π-rotation, 2 π/2-rotations C 3v : 3 reflections, 3 2π/3-rotations C 6v : 6 reflections, 1 π, 2 π/3, 2 π/6 rotations States in the Hilbert space fall into a finite number of irreducible representations (irreps) of G The Hamiltonian H is block diagonal w.r.t. to irreps.
25 Group characters C 2 E C 2 A 1 1 B 1 1 C 2v e c 2 σ 1 σ 2 A A B B C 4v e c 2 2c 4 2σ 1 2σ 2 A A B B E
26 Taking advantage of cluster symmetries... order of the group Reduces the dimension of the Hilbert space by G Accelerates the convergence of the Lanczos algorithm Reduces the number of Band Lanczos starting vectors by G But: complicates coding of the basis states Make use of the projection operator: dimension of irrep. See, e.g. Poilblanc & Laflorencie cond-mat/ P (α) = d α χ G g (α) g g group character
27 Taking advantage of cluster symmetries (2) Need new basis states, made of sets of binary states related by the group action: ψ = d α G g χ (α) Then matrix elements take the form ψ 2 H ψ 1 = d α G fermionic phase g g b g b = φg(b) gb g χ (α) h φ g (b) gb 2 H b 1 When computing the Green function, one needs to use combinations of creation operators that fall into group representations. Ex (4 1): c (A) 1 = c 1 + c 4 c (A) 2 = c 2 + c 3 c (B) 1 = c 1 c c (B) 2 = c 2 c 3
28 Taking advantage of cluster symmetries (3) Example : number of matrix elements of the kinetic energy operator (Nearest neighbor) on a 3 4 cluster with C 2v symmetry: A 1 A 2 B 1 B 2 dim. 213, , , , 248 value , 640 6, 208 7, 584 5, , 983, 264 2, 936, 144 2, 884, 832 2, 911, , , 168 1, 050, 432 1, 021, , 088 2, 304 3, 232 2,
29 Large dimensions : need for parallelization Memory needs exceed single cpu capacity beyond L 14 A half-filled 16-site system has dimension 165,636, GB for a state vector. Need to distribute the problem over many processors The main task is matrix-vector multiplication: a a a a a a a a a a a a a a a a a a a a a a a a
30 Part II Cluster Perturbation Theory
31 Clusters and superlattices e 2 e 1 K k k (π, π) (0, 0) ( π, π) 10-site cluster Reduced Brillouin zone
32 Basic Idea lattice Hamiltonian H hopping matrix cluster Hamiltonian = H + V inter-cluster hopping terms cluster hopping matrix t = t + V inter-cluster hopping Treat V at lowest order in Perturbation theory At this order, the Green function is G 1 (ω) = G 1 (ω) V cluster Green function C. Gros and R. Valenti, Phys. Rev. B 48, 418 (1993) D. Sénéchal, D. Perez, and M. Pioro-Ladrière. Phys. Rev. Lett. 84, 522 (2000)
33 Interlude : Fourier transforms i, j : lattice site index k : full wavevector m, n : lattice site index k : redcued wavevector a, b : cluster site index K : cluster wavevector f j = e ik r j f(k) f(k) = 1 e ik r j f j N k j f m = k e i k r m f( k) f( k) = L e i k r m f m N m f a = 1 e ik ra f K f K = 1 e ik ra f a L L K a
34 Basic Idea (cont.) More accurate formulation G 1 ( k, ω) = G 1 (ω) V( k). But G 1 = ω t Σ G 1 0 = ω t V, Thus : lattice self-energy is approximated as the cluster self-energy G 1 ( k, ω) = G 1 0 ( k, ω) Σ(ω), Example : 2-site cluster (1D): ( ) t 0 1 = t 1 0 V( k) = t ( ) 0 e 2i k e 2i k 0
35 Periodization CPT breaks translation invariance, which needs to be restored: G cpt (k, ω) = 1 e ik (ra rb) G ab ( k, ω). L Periodizing the Green function vs the self-energy (1D case): a,b (π) (π) (π/2) (π/2) (0) ω Green function periodization (0) ω Self-energy periodization
36 One-dimensional example Evolution of spectral function with increasing U/t:
37 Interlude : Relation with spectral function Lehmann representation: A(k, ω) = 2 lim Im G(k, ω + iη) η 0 + G αβ (ω) = 1 Ω c α m m c ω E m m + E β Ω 0 + Ω c β n 1 n c α Ω ω + E n n E 0 But : lim Im 1 η 0 + ω + iη = lim Therefore : A(k, ω) = m η 0 + η ω 2 + η 2 = πδ(ω) m c k Ω 2 2πδ(ω E m + E 0 ) + n n c k Ω 2 2πδ(ω + E n E 0 )
38 h-doped cuprates : Pseudogap from CPT Sénéchal and Tremblay., Phys. Rev. Lett. 92, (2004)
39 e-doped cuprates : Pseudogap from CPT Sénéchal and Tremblay., Phys. Rev. Lett. 92, (2004)
40 CPT : characteristics Exact at U = 0 Exact at t ij = 0 Exact short-range correlations Allows all values of the wavevector But : No long-range order Controlled by the size of the cluster
41 Part III The self-energy functional approach
42 Motivation CPT cannot describe broken symmetry states, because of the finite cluster size Idea : add a Weiss field term to the cluster Hamiltonian H, e.g., for antiferromagnetism: H M = M a (π, π) e iq r a (n a n a ) This term favors AF order, but does not appear in H, and must be subtracted from V Need a principle to set the value of M : energy minimization? Better : Potthoff s self-energy functional approach
43 The Potthoff variational principle Variational principle for the Green function: Ω t [G] = Φ[G] Tr ((G 1 0t Where Φ[G] is the Luttinger-Ward functional: G 1 )G) + Tr ln( G). Tr (A) = P ω,α Aαα(ω) Φ = with the property M. Potthoff, Eur. Phys. J. B 32, 429?436 (2003) δφ[g] δg = Σ
44 The Potthoff variational principle (2) Here, Tr means a sum over frequencies, site indices (or wavevectors) and spin/band indices. The functional is stationary at the physical Green function (Euler eq.): δω t [G] δg = Σ G 1 0t + G 1 = 0. Approximation schemes: Type I : Simplify the Euler equation Type II : Approximate the functional (Hartree-Fock, FLEX) Type III : Restrict the variational space, but keep the functional exact
45 The Potthoff variational principle (3) Potthoff : Use the self-energy rather than the Green function Ω t [Σ] = F [Σ] Tr ln( G 1 0t + Σ) F [Σ] = Φ[G] Tr (ΣG) F is the Legendre transform of Φ: δf [Σ] δσ New Euler equation: = δφ[g] δg[σ] δg δσ ΣδG[Σ] δσ G = G δω t [Σ] δσ = G + (G 1 0t Σ) 1 = 0 At the physical self-energy, Ω t [Σ] is the thermodynamic grand potential
46 The Reference System To evaluate F, use its universal character : its functional form depends only on the interaction. Introduce a reference rystem H, which differs from H by one-body terms only (example : the cluster Hamiltonian) Suppose H can be solved exactly. Then, at the physical self-energy Σ of H, Ω = F [Σ] Tr ln( G ) by eliminating F : Ω t [Σ] = Ω + Tr ln( G ) Tr ln( G 1 0t + Σ) = Ω + Tr ln( G ) Tr ln( G) = Ω Tr ln(1 VG )
47 The Potthoff functional Making the trace explicit, one finds Ω t [Σ] = Ω T [ ] tr ln 1 V( k)g ( k, ω) ω k = Ω T [ ] ln det 1 V( k)g ( k, ω) ω The sum over frequencies is to be performed over Matsubara frequencies (or an integral along the imaginary axis at T = 0). The variation is done over one-body parameters of the cluster Hamiltonian H k In particular, the Weiss field M is to be varied until Ω is stationary
48 Calculating the functional I : exact form It can be shown that Tr ln( G) = T m poles of G ln(1 + e βω m) + T m ln(1 + e βζ m) zeros of G Use the Lehmann representation of the GF: G µν(ω) = r M. Potthoff, Eur. Phys. J. B, 36:335 (2003) Q µr Q νr G(ω) = Q 1 ω ω r ω Λ Q diagonal(ω r)
49 Calculating the functional I : exact form (2) A similar representation holds for the CPT Green function 1 G( k, ω) = G 1 V( k) = 1 [ ] 1 Q 1 ω Λ Q V( k) 1 = Q ω L( k) Q L( k) = Λ + Q V( k)q Let ω r ( k) be the eigenvalues of L( k). Then Ω(x) = Ω (x) ω r<0 variational parameters M. Aichhorn et al., Phys. Rev. B 74 : (2006) ω r + L N ω r ( k) k ω r( k)<0
50 Calculating the functional II : numerical integral Except for very small clusters (L 4), it is much faster to perform a numerical integration over frequencies: Ω(x) = Ω (x) 0 dx π L ln det(1 V( k)g (ix)) L(µ µ ) N k D. Sénéchal, proceedings of HPCS 2008, IEEE (2008)
51 Evaluation of integrals For frequency integrals: Gaussian integration on three segments For wavevector integrals, adaptive mesh of points: Start with a coarse, regular grid On each plaquette, compare 4 and 9 point Gaussian integrals. Subdivide into 4 sub-plaquettes if necessary. Easy with recursive calls
52 Part IV The Variational Cluster Approximation
53 Basic Idea Set up a superlattice of clusters Choose a set of variational parameters, e.g. Weiss fields for broken symmetries Set up the calculation of the Potthoff functional: Ω t [Σ] = Ω T L N ω k [ ] ln det 1 V( k)g ( k, ω) Use an optimization method to find the stationary points Adopt the cluster self-energy associated with the stationary point with the lowest Ω and use it as in CPT
54 Example : Néel Antiferromagnetism Used the Weiss field H M = M a (π, π) e iq r a (n a n a ) Profile of Ω for the half-filled, square lattice Hubbard model: Ω U = 16 U = 8 U = 4 U = M Ω L = L = M
55 Example : Néel Antiferromagnetism (2) order parameter M (4 ) ordering energy U Best scaling factor : q = number of links 2 number of sites M Order Parameter x1 2x1 links 2L 2x1 2x2 2x1 2x2 1 1 L 2x3 2x4 B10 2x3 2x4 B10 2x2 3x4 4x4 2x2 3x4 4x4 2x3 2x4 B10 3x4 4 2x3 2x4 B10 3x
56 Example clusters 3 4 B10
57 Superconductivity Need to add a pairing field O sc = ij c i c j + H.c ij s-wave pairing: ij = δ ij d x 2 y 2 pairing: { 1 if ri r j = ±ˆx ij = 1 if r i r j = ±ŷ d xy pairing: { 1 if ri r j = ±(ˆx + ŷ) ij = 1 if r i r j = ±(ˆx ŷ) Ω d x 2 y 2 extended s-wave dxy cluster s-wave U = 8, µ =
58 Superconductivity (2) Pairing fields violate particle number conservation The Hilbert space is enlarged to encompass all particle numbers with a given spin In practice, on uses the Nambu formalism, i.e., particle-hole transformation on the spin-down sector : c a = c a and d a = c a Then the Hamiltonian looks like it conserves particle number, but not spin.
59 Superconductivity and Antiferromagnetism in the cuprates One-band Hubbard model for the cuprates: t = 0.3, t = 0.2, U = 8: order parameter pure d x 2 y 2 pure Néel coex. d x 2 y 2 coex. Néel n M. Guillot, MSc thesis, Univ. de Sherbrooke (2007)
60 Thermodynamic consistency The electron density n may be calculated either as n = Tr G or n = Ω µ The two methods give different results, except if the cluster chemical potential µ is treated like a variational parameter: cluster U = 8 normal state 0.8 n Ω µ cons. TrG TrG cons. Ω µ µ
61 Optimization procedure Need to find the saddle points of Ω(x) with the least possible evaluations of Ω(x) Use the Newton-Raphson algorithm: Evaluate Ω at a number of points at and around x0 that just fits a quadratic form Move to the stationary point x 1 of that quadratic form and repeat Stop when x i x i 1, or the numerical gradient Ω, converges The NR method is not robust : it converges fast when started close enough to the solution Proceed adiabatically through external parameter space (e.g. as function of U or µ)
62 Example: Homogeneous coexistence of dsc and AF orders
63 Example: dsc on a 4 4 cluster, spectral function spectral function : 4x4_C2v, U=8, mu=1.5, tx=1, Dx=0.035, (spin up) (27/5/2008) (0.5,0.5) (1,0) (1,0.5) (1,1) (0.5,0.5) (0,0) (0.5,0)
64 Example: dsc on a 4 4 cluster, Fermi surface plot
65 VCA vs Mean-Field Theory Differs from Mean-Field Theory: Interaction is left intact, it is not factorized Retains exact short-range correlations Weiss field order parameter More stringent that MFT Controlled by the cluster size Similarities with MFT: No long-range fluctuations (no disorder from Goldstone modes) Yet : no LRO for Néel AF in one dimension Need to compare different orders yet : they may be placed in competition / coexistence
66 Part V Cluster Dynamical Mean Field Theory
67 Basic Idea To add variational degrees of freedom in the form of a bath of uncorrelated sites H = µ,ν t µν c µc ν + U a n a n a + µ,α θ µα (c µa α + H.c.) + hybridization matrix α ε α a αa a bath energies 7 8 ε 2 t 2 t 2 ε 2 ε 2 t 2 t 2 ε ε 1 t 1 t 1 ε 1 ε 1 t 1 t 1 ε
68 The hybridization function If we trace over the bath degrees of freedom, the cluster Green function takes the form Γ(ω) is the hybridization function: G 1 = ω t Γ(ω) Σ(ω) Γ µν (ω) = α θ µα θ να ω ε α
69 The hybridization function (2) Proof: (U = 0) G 1 full (ω) = 1 ω T T = ( ω t θ θ ω ε ) Given A = G 1 full, need to find B 1 11 : ( ) ( A 11 A 12 B 11 B 12 = A 21 A 22 B 21 B 22 Simple manipulations lead to ) 1 ( A11 A 12 A 1 22 A ) 21 B11 = 1 G 1 = ω t θ 1 ω ε θ U 0 : simply add the free energy (by definition)
70 The hybridization function (3) Γ(ω) embodies the effect of the rest of the lattice on the cluster, in some effective dynamics. The action would take the form β β S = dτ dτ c µ(τ)gµν 1 (τ τ )c ν (τ ) 0 0 µν + U d β 0 τ n a (τ)n a (τ) a where G (iω n ) = β 0 e iωnτ G (τ) G 1 µν (iω n ) = iω n δ µν t µν Γ µν (iω n )
71 Baths and the SFA The Potthoff functional approach carries over unchanged in the presence of a bath The bath makes a contribution to the Potthoff functional: Ω bath = ε α<0 On can in principle use the same methods as in VCA The presence of the bath increases the resolution of the approach in the time domain, at the cost of spatial resolution, for a fixed total number of sites (cluster + bath). ε α
72 The CDMFT Procedure 1. Start with a guess value of (θ µα, ε α ). 2. Calculate the cluster Green function G(ω) (ED). 3. Calculate the superlattice-averaged Green function Ḡ(ω) = k 1 G 1 0 ( k) Σ(ω) and G 1 0 (ω) = Ḡ 1 + Σ(ω) 4. Minimize the following distance function: d = ω,ν,ν over the set of bath parameters. 5. Go back to step (2) until convergence. ( ω + µ t Γ(ω) G0 1 (ω) ) 2 νν
73 The CDMFT Procedure (2) Initial guess for Γ Cluster Solver: Compute G Ḡ = k[g 1 0 ( k) Σ(ω)] 1 G 1 0 = Ḡ 1 + Σ update Γ by minimizing d No Γ converged? Yes Exit
74 Example : the 1D Hubbard model 1 n exact VCA β = 20, ω c = 1.5 β = 20, ω c = 10, 1/ω weight β = 50, ω c = 1.5 β = 50, ω c = µ
75 Example : the 1D Hubbard model (2) ε t µ µ
76 Example : dsc and AF in the 2D Hubbard model Nine bath parameters Homogeneous coexistence of d x 2 y2 SC and Néel AF order parameter dsc ( 5) AF n
77 Example : The Mott transition The CDMFT is well suited to detect the Mott transition This transition manifests itself as a jump in the double occupancy n n In an exact SFA solution : discontinuity in the bath parameters (first order transition). in CDMFT : hysteresis is possible, because of the method s own dynamics for finding solutions D (a) 0.05 solid: t /t=0.7 dashed: t /t= U/t B. Kyung and A.-M. S. Tremblay. Physical U c /t Review Letters, 97 : (2006)
78 QUESTIONS?
Quantum Cluster Methods: An introduction
Quantum Cluster Methods: An introduction David Sénéchal Département de physique, Université de Sherbrooke International summer school on New trends in computational approaches for many-body systems May
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 informationarxiv: v2 [cond-mat.str-el] 24 Nov 2010
An introduction to quantum cluster methods Lectures given at the CIFAR - PITP International Summer School on Numerical Methods for Correlated Systems in Condensed Matter, Sherbrooke, Canada, May 26 June
More informationDynamical Mean Field Theory. Quantum Cluster Methods
Dynamical Mean Field Theory and Quantum Cluster Methods David Sénéchal Département de physique & Institut Quantique Université de Sherbrooke International Summer School on Computational Quantum Materials
More informationRecent advances on Hubbard models using quantum cluster methods
Recent advances on Hubbard models using quantum cluster methods David Sénéchal Département de physique Faculté des sciences Université de Sherbrooke 213 CAP Congress, Université de Montréal May 27, 213
More informationIntroduction to DMFT
Introduction to DMFT Lecture 2 : DMFT formalism 1 Toulouse, May 25th 2007 O. Parcollet 1. Derivation of the DMFT equations 2. Impurity solvers. 1 Derivation of DMFT equations 2 Cavity method. Large dimension
More informationarxiv:cond-mat/ v2 [cond-mat.str-el] 24 Feb 2006
Applications of Cluster Perturbation Theory Using Quantum Monte Carlo Data arxiv:cond-mat/0512406v2 [cond-mat.str-el] 24 Feb 2006 Fei Lin, Erik S. Sørensen, Catherine Kallin and A. John Berlinsky Department
More information9 Making Use of Self-Energy Functionals: The Variational Cluster Approximation
9 Making Use of Self-Energy Functionals: The Variational Cluster Approximation Contents Michael Potthoff I. Institut für Theoretische Physik Universität Hamburg Motivation 2 2 The cluster approach 4 2.
More informationPseudogap and high-temperature superconductivity from weak to strong coupling. Towards a quantitative theory Review Article
LOW TEMPERATURE PHYSICS VOLUME 32, NUMBER 4-5 APRIL-MAY 2006 Pseudogap and high-temperature superconductivity from weak to strong coupling. Towards a quantitative theory Review Article A.-M. S. Tremblay,
More informationNumerical Methods in Many-body Physics
Numerical Methods in Many-body Physics Reinhard M. Noack Philipps-Universität Marburg Exchange Lecture BME, Budapest, Spring 2007 International Research Training Group 790 Electron-Electron Interactions
More informationO. Parcollet CEA-Saclay FRANCE
Cluster Dynamical Mean Field Analysis of the Mott transition O. Parcollet CEA-Saclay FRANCE Dynamical Breakup of the Fermi Surface in a doped Mott Insulator M. Civelli, M. Capone, S. S. Kancharla, O.P.,
More informationw2dynamics : operation and applications
w2dynamics : operation and applications Giorgio Sangiovanni ERC Kick-off Meeting, 2.9.2013 Hackers Nico Parragh (Uni Wü) Markus Wallerberger (TU) Patrik Gunacker (TU) Andreas Hausoel (Uni Wü) A solver
More informationStrange metal from local quantum chaos
Strange metal from local quantum chaos John McGreevy (UCSD) hello based on work with Daniel Ben-Zion (UCSD) 2017-08-26 Compressible states of fermions at finite density The metallic states that we understand
More informationQuantum Lattice Models & Introduction to Exact Diagonalization
Quantum Lattice Models & Introduction to Exact Diagonalization H! = E! Andreas Läuchli IRRMA EPF Lausanne ALPS User Workshop CSCS Manno, 28/9/2004 Outline of this lecture: Quantum Lattice Models Lattices
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 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 informationMott transition : beyond Dynamical Mean Field Theory
Mott transition : beyond Dynamical Mean Field Theory O. Parcollet 1. Cluster methods. 2. CDMFT 3. Mott transition in frustrated systems : hot-cold spots. Coll: G. Biroli (SPhT), G. Kotliar (Rutgers) Ref:
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 informationLecture 11: Long-wavelength expansion in the Neel state Energetic terms
Lecture 11: Long-wavelength expansion in the Neel state Energetic terms In the last class we derived the low energy effective Hamiltonian for a Mott insulator. This derivation is an example of the kind
More informationCluster Extensions to the Dynamical Mean-Field Theory
Thomas Pruschke Institut für Theoretische Physik Universität Göttingen Cluster Extensions to the Dynamical Mean-Field Theory 1. Why cluster methods? Thomas Pruschke Institut für Theoretische Physik Universität
More information10 Strongly Correlated Superconductivity
10 Strongly Correlated Superconductivity André-Marie S. Tremblay Département de physique Regroupement québécois sur les matériaux de pointe Canadian Institute for Advanced Research Sherbrooke, Québec J1K
More informationExact diagonalization methods
Summer School on Computational Statistical Physics August 4-11, 2010, NCCU, Taipei, Taiwan Exact diagonalization methods Anders W. Sandvik, Boston University Representation of states in the computer bit
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 informationTopological order in the pseudogap metal
HARVARD Topological order in the pseudogap metal High Temperature Superconductivity Unifying Themes in Diverse Materials 2018 Aspen Winter Conference Aspen Center for Physics Subir Sachdev January 16,
More informationNumerical Studies of the 2D Hubbard Model
arxiv:cond-mat/0610710v1 [cond-mat.str-el] 25 Oct 2006 Numerical Studies of the 2D Hubbard Model D.J. Scalapino Department of Physics, University of California, Santa Barbara, CA 93106-9530, USA Abstract
More informationElectronic correlations in models and materials. Jan Kuneš
Electronic correlations in models and materials Jan Kuneš Outline Dynamical-mean field theory Implementation (impurity problem) Single-band Hubbard model MnO under pressure moment collapse metal-insulator
More informationANTIFERROMAGNETIC EXCHANGE AND SPIN-FLUCTUATION PAIRING IN CUPRATES
ANTIFERROMAGNETIC EXCHANGE AND SPIN-FLUCTUATION PAIRING IN CUPRATES N.M.Plakida Joint Institute for Nuclear Research, Dubna, Russia CORPES, Dresden, 26.05.2005 Publications and collaborators: N.M. Plakida,
More informationEffect of next-nearest-neighbour interaction on d x 2 y2-wave superconducting phase in 2D t-j model
PRAMANA c Indian Academy of Sciences Vol. 74, No. 1 journal of January 2010 physics pp. 115 121 Effect of next-nearest-neighbour interaction on d x 2 y2-wave superconducting phase in 2D t-j model N S MONDAL
More information2. The interaction between an electron gas and the potential field is given by
68A. Solutions to Exercises II. April 5. Carry out a Hubbard Stratonovich transformation for the following interaction Hamiltonians, factorizing the terms in brackets, and write one or two sentences interpreting
More informationPreface Introduction to the electron liquid
Table of Preface page xvii 1 Introduction to the electron liquid 1 1.1 A tale of many electrons 1 1.2 Where the electrons roam: physical realizations of the electron liquid 5 1.2.1 Three dimensions 5 1.2.2
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 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 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 informationLecture Models for heavy-ion collisions (Part III): transport models. SS2016: Dynamical models for relativistic heavy-ion collisions
Lecture Models for heavy-ion collisions (Part III: transport models SS06: Dynamical models for relativistic heavy-ion collisions Quantum mechanical description of the many-body system Dynamics of heavy-ion
More informationLecture 1: The Equilibrium Green Function Method
Lecture 1: The Equilibrium Green Function Method Mark Jarrell April 27, 2011 Contents 1 Why Green functions? 2 2 Different types of Green functions 4 2.1 Retarded, advanced, time ordered and Matsubara
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 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 informationGreen s functions: calculation methods
M. A. Gusmão IF-UFRGS 1 FIP161 Text 13 Green s functions: calculation methods Equations of motion One of the usual methods for evaluating Green s functions starts by writing an equation of motion for the
More informationBlock Iterative Eigensolvers for Sequences of Dense Correlated Eigenvalue Problems
Mitglied der Helmholtz-Gemeinschaft Block Iterative Eigensolvers for Sequences of Dense Correlated Eigenvalue Problems Birkbeck University, London, June the 29th 2012 Edoardo Di Napoli Motivation and Goals
More informationFinite Temperature Field Theory
Finite Temperature Field Theory Dietrich Bödeker, Universität Bielefeld 1. Thermodynamics (better: thermo-statics) (a) Imaginary time formalism (b) free energy: scalar particles, resummation i. pedestrian
More informationPHYSICAL REVIEW B 80,
PHYSICAL REVIEW B 8, 6526 29 Finite-temperature exact diagonalization cluster dynamical mean-field study of the two-dimensional Hubbard model: Pseudogap, non-fermi-liquid behavior, and particle-hole asymmetry
More informationLecture 5. Hartree-Fock Theory. WS2010/11: Introduction to Nuclear and Particle Physics
Lecture 5 Hartree-Fock Theory WS2010/11: Introduction to Nuclear and Particle Physics Particle-number representation: General formalism The simplest starting point for a many-body state is a system of
More informationPhysics 127b: Statistical Mechanics. Landau Theory of Second Order Phase Transitions. Order Parameter
Physics 127b: Statistical Mechanics Landau Theory of Second Order Phase Transitions Order Parameter Second order phase transitions occur when a new state of reduced symmetry develops continuously from
More informationNucleons in Nuclei: Interactions, Geometry, Symmetries
Nucleons in Nuclei: Interactions, Geometry, Symmetries Jerzy DUDEK Department of Subatomic Research, CNRS/IN 2 P 3 and University of Strasbourg, F-67037 Strasbourg, FRANCE September 28, 2010 Mathematial
More informationDynamical Mean Field within Iterative Perturbation Theory
Vol. 111 (2007) ACTA PHYSICA POLONICA A No. 5 Proceedings of the XII National School Correlated Electron Systems..., Ustroń 2006 Dynamical Mean Field within Iterative Perturbation Theory B. Radzimirski
More informationQuantum criticality of Fermi surfaces
Quantum criticality of Fermi surfaces Subir Sachdev Physics 268br, Spring 2018 HARVARD Quantum criticality of Ising-nematic ordering in a metal y Occupied states x Empty states A metal with a Fermi surface
More informationHydrodynamic Modes of Incoherent Black Holes
Hydrodynamic Modes of Incoherent Black Holes Vaios Ziogas Durham University Based on work in collaboration with A. Donos, J. Gauntlett [arxiv: 1707.xxxxx, 170x.xxxxx] 9th Crete Regional Meeting on String
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 informationMean-Field Theory: HF and BCS
Mean-Field Theory: HF and BCS Erik Koch Institute for Advanced Simulation, Jülich N (x) = p N! ' (x ) ' 2 (x ) ' N (x ) ' (x 2 ) ' 2 (x 2 ) ' N (x 2 )........ ' (x N ) ' 2 (x N ) ' N (x N ) Slater determinant
More informationAndré-Marie Tremblay
André-Marie Tremblay CENTRE DE RECHERCHE SUR LES PROPRIÉTÉS ÉLECTRONIQUES DE MATÉRIAUX AVANCÉS Sponsors: Pairing vs antiferromagnetism in the cuprates 1. Motivation, model 2. The standard approach : -
More information7 The Finite Temperature Lanczos Method and its Applications
7 The Finite Temperature Lanczos Method and its Applications Peter Prelovšek J. Stefan Institute Jamova 39, SI-1000 Ljubljana, Slovenia Contents 1 Introduction 2 2 Exact diagonalization and Lanczos method
More informationCollective Effects. Equilibrium and Nonequilibrium Physics
Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 3, 3 March 2006 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech
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 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 informationComputational strongly correlated materials R. Torsten Clay Physics & Astronomy
Computational strongly correlated materials R. Torsten Clay Physics & Astronomy Current/recent students Saurabh Dayal (current PhD student) Wasanthi De Silva (new grad student 212) Jeong-Pil Song (finished
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 informationCOMPETITION BETWEEN FILLED AND HALF-FILLED STRIPES IN CUPRATES AND NICKELATES
COMPETITION BETWEEN FILLED AND HALF-FILLED STRIPES IN CUPRATES AND NICKELATES Raymond Frésard Andrzej M. Oleś and Marcin Raczkowski Laboratoire Crismat UMR CNRS-ENSICAEN (ISMRA) 6508 Caen France Marian
More informationDiagrammatic extensions of (E)DMFT: Dual boson
Diagrammatic extensions of (E)DMFT: Dual boson IPhT, CEA Saclay, France ISSP, June 25, 2014 Collaborators Mikhail Katsnelson (University of Nijmegen, The Netherlands) Alexander Lichtenstein (University
More informationStrong Correlation Effects in Fullerene Molecules and Solids
Strong Correlation Effects in Fullerene Molecules and Solids Fei Lin Physics Department, Virginia Tech, Blacksburg, VA 2461 Fei Lin (Virginia Tech) Correlations in Fullerene SESAPS 211, Roanoke, VA 1 /
More informationPhase transitions and critical phenomena
Phase transitions and critical phenomena Classification of phase transitions. Discontinous (st order) transitions Summary week -5 st derivatives of thermodynamic potentials jump discontinously, e.g. (
More information26 Group Theory Basics
26 Group Theory Basics 1. Reference: Group Theory and Quantum Mechanics by Michael Tinkham. 2. We said earlier that we will go looking for the set of operators that commute with the molecular Hamiltonian.
More informationDynamical mean field approach to correlated lattice systems in and out of equilibrium
Dynamical mean field approach to correlated lattice systems in and out of equilibrium Philipp Werner University of Fribourg, Switzerland Kyoto, December 2013 Overview Dynamical mean field approximation
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationFirst-order Mott transition at zero temperature in two dimensions: Variational plaquette study
January 2009 EPL, 85 (2009) 17002 doi: 10.1209/0295-5075/85/17002 www.epljournal.org First-order Mott transition at zero temperature in two dimensions: ariational plaquette study M. Balzer 1, B. Kyung
More informationQuantum Monte Carlo study of strongly correlated electrons: Cellular dynamical mean-field theory
Quantum Monte Carlo study of strongly correlated electrons: Cellular dynamical mean-field theory B. Kyung, 1 G. Kotliar, 2 and A.-M. S. Tremblay 1 1 Département de physique and Regroupement québécois sur
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 informationExamination paper for TFY4210 Quantum theory of many-particle systems
Department of Physics Examination paper for TFY4210 Quantum theory of many-particle systems Academic contact during examination: Associate Professor John Ove Fjærestad Phone: 97 94 00 36 Examination date:
More informationSection IExact diagonalisations and Lanczos methodscomparison with other methods p.1
Section I Exact diagonalisations and Lanczos methods Comparison with other methods Section IExact diagonalisations and Lanczos methodscomparison with other methods p.1 Outline 1. Power method & Lanczos
More informationDiagrammatic Monte Carlo methods for Fermions
Diagrammatic Monte Carlo methods for Fermions Philipp Werner Department of Physics, Columbia University PRL 97, 7645 (26) PRB 74, 15517 (26) PRB 75, 8518 (27) PRB 76, 235123 (27) PRL 99, 12645 (27) PRL
More informationErik Koch Computational Materials Science German Research School for Simulation Sciences
8 The Lanczos Method Erik Koch Computational Materials Science German Research School for Simulation Sciences Contents 1 Lanczos Method 3 1.1 Krylov space................................... 5 1.2 Spectral
More informationLOCAL MOMENTS NEAR THE METAL-INSULATOR TRANSITION
LOCAL MOMENTS NEAR THE METAL-INSULATOR TRANSITION Subir Sachdev Center for Theoretical Physics, P.O. Box 6666 Yale University, New Haven, CT 06511 This paper reviews recent progress in understanding the
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 informationM. A. Gusmão IF-UFRGS
M. A. Gusmão IF-UFRGS - 217 1 FIP164-217/2 Text 9 Mean-field approximation - II Heisenberg Hamiltonian in wave-vector space As we saw in Text 8, the uniform susceptibility does not diverge in the case
More informationMOTTNESS AND STRONG COUPLING
MOTTNESS AND STRONG COUPLING ROB LEIGH UNIVERSITY OF ILLINOIS Rutgers University April 2008 based on various papers with Philip Phillips and Ting-Pong Choy PRL 99 (2007) 046404 PRB 77 (2008) 014512 PRB
More informationMagnetism and Superconductivity in Decorated Lattices
Magnetism and Superconductivity in Decorated Lattices Mott Insulators and Antiferromagnetism- The Hubbard Hamiltonian Illustration: The Square Lattice Bipartite doesn t mean N A = N B : The Lieb Lattice
More informationProperties of monopole operators in 3d gauge theories
Properties of monopole operators in 3d gauge theories Silviu S. Pufu Princeton University Based on: arxiv:1303.6125 arxiv:1309.1160 (with Ethan Dyer and Mark Mezei) work in progress with Ethan Dyer, Mark
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 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 informationFirst-Principles Calculation of Exchange Interactions
Chapter 2 First-Principles Calculation of Exchange Interactions Before introducing the first-principles methods for the calculation of exchange interactions in magnetic systems we will briefly review two
More informationSmall and large Fermi surfaces in metals with local moments
Small and large Fermi surfaces in metals with local moments T. Senthil (MIT) Subir Sachdev Matthias Vojta (Augsburg) cond-mat/0209144 Transparencies online at http://pantheon.yale.edu/~subir Luttinger
More informationThe frustrated Hubbard model on the Bethe lattice - an investigation using the self-energy functional approach
The frustrated Hubbard model on the Bethe lattice - an investigation using the self-energy functional approach von Martin Eckstein Diplomarbeit im Fach Physik vorgelegt der Mathematisch-Naturwissenschaftlichen
More informationLandau-Fermi liquid theory
Landau-Fermi liquid theory Shreyas Patankar Chennai Mathematical Institute Abstract We study the basic properties of Landau s theory of a system of interacting fermions (a Fermi liquid). The main feature
More informationAn efficient impurity-solver for the dynamical mean field theory algorithm
Papers in Physics, vol. 9, art. 95 (217) www.papersinphysics.org Received: 31 March 217, Accepted: 6 June 217 Edited by: D. Domínguez Reviewed by: A. Feiguin, Northeastern University, Boston, United States.
More informationTopological sector of two dimensional φ 4 theory in Discrete Light Cone Quantization p. 1/3
Topological sector of two dimensional φ 4 theory in Discrete Light Cone Quantization p. 1/3 Topological sector of two dimensional φ 4 theory in Discrete Light Cone Quantization A. H (Theory Division, SINP,
More informationarxiv:cond-mat/ v3 [cond-mat.str-el] 11 Dec 1997
Non-perturbative many-body approach to the Hubbard model and single-particle pseudogap. arxiv:cond-mat/9702188v3 [cond-mat.str-el] 11 Dec 1997 Y.M. Vilk Materials Science Division, Bldg. 223, 9700 S. Case
More informationDiagrammatic Monte Carlo simulation of quantum impurity models
Diagrammatic Monte Carlo simulation of quantum impurity models Philipp Werner ETH Zurich IPAM, UCLA, Jan. 2009 Outline Continuous-time auxiliary field method (CT-AUX) Weak coupling expansion and auxiliary
More informationChapter 7 Iterative Techniques in Matrix Algebra
Chapter 7 Iterative Techniques in Matrix Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 128B Numerical Analysis Vector Norms Definition
More informationHigh-Temperature Criticality in Strongly Constrained Quantum Systems
High-Temperature Criticality in Strongly Constrained Quantum Systems Claudio Chamon Collaborators: Claudio Castelnovo - BU Christopher Mudry - PSI, Switzerland Pierre Pujol - ENS Lyon, France PRB 2006
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 IPAM QS2009 Workshop,
More informationɛ(k) = h2 k 2 2m, k F = (3π 2 n) 1/3
4D-XY Quantum Criticality in Underdoped High-T c cuprates M. Franz University of British Columbia franz@physics.ubc.ca February 22, 2005 In collaboration with: A.P. Iyengar (theory) D.P. Broun, D.A. Bonn
More informationRepresentation theory & the Hubbard model
Representation theory & the Hubbard model Simon Mayer March 17, 2015 Outline 1. The Hubbard model 2. Representation theory of the symmetric group S n 3. Representation theory of the special unitary group
More informationNon-Perturbative Many-Body Approach to the Hubbard Model and Single-Particle Pseudogap
J. Phys. I France 7 (1997) 1309 1368 NOVEMBER 1997, PAGE 1309 Non-Perturbative Many-Body Approach to the Hubbard Model and Single-Particle Pseudogap Y.M. Vilk ( 1,2 ) and A.-M.S. Tremblay ( 2, ) ( 1 )
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 informationElectronic charge reconstruction of doped Mott insulators in multilayered nanostructures
Electronic charge reconstruction of doped Mott insulators in multilayered nanostructures Ling Chen and J. K. Freericks* Department of Physics, Georgetown University, Washington, District of Columbia 20057,
More informationWavefunction and electronic struture in solids: Bloch functions, Fermi level and other concepts.
Wavefunction and electronic struture in solids: Bloch functions, Fermi level and other concepts. Silvia Casassa Università degli Studi di Torino July 12, 2017 Minnesota Workshop on ab initio MSC Symmetry
More informationarxiv:cond-mat/ v1 [cond-mat.str-el] 29 Mar 2005
The Antiferromagnetic Heisenberg Model on Clusters with Icosahedral Symmetry arxiv:cond-mat/0503658v1 [cond-mat.str-el] 9 Mar 005 N.P. Konstantinidis Department of Physics and Department of Mathematics,
More informationElectrons in a periodic potential
Chapter 3 Electrons in a periodic potential 3.1 Bloch s theorem. We consider in this chapter electrons under the influence of a static, periodic potential V (x), i.e. such that it fulfills V (x) = V (x
More informationLow-lying excitation spectrum of quantum many-body systems
PHYSICAL REVIEW B VOLUME 54, NUMBER 21 1 DECEMBER 1996-I Low-lying excitation spectrum of quantum many-body systems P. Monthoux Department of Physics and National High Magnetic Field Laboratory, Florida
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 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