Numerical Methods in Many-body Physics
|
|
- Amos Fowler
- 6 years ago
- Views:
Transcription
1 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 in Solids
2 Literature: R.M. Noack and S.R. Manmana, Diagonalization- and Numerical Renormalization-Group- Based Methods for Interacting Quantum Systems, in Lectures on the Physics of Highly Correlated Electron Systems IX, AIP Conference Proceedings 789, AIP, New York, 2005, and cond-mat/ Density Matrix Renormalization: A New Numerical Method in Physics, Lect. Notes Phys. 528, Eds. I. Peschel, X. Wang, M. Kaulke, and K. Hallberg, Springer, Berlin, E. Dagotto, Correlated electrons in high-temperature superconductors, Rev. Mod. Phys. 66, (1994). A. Booten and H. van der Vorst, Cracking large-scale eigenvalue problems, part I: Algorithms, Computers in Phys. 10, No. 3, p. 239 (May/June 1996); Cracking large-scale eigenvalue problems, part II: Implementations, Computers in Phys. 10, No. 4, p. 331 (July/August 1996). N. Laflorencie and D. Poilblanc, Simulations of pure and doped low-dimensional spin-1/2 gapped systems in Quantum Magnetism, Lect. Notes Phys. 645, Springer, Berlin, 2004, pp (2004), cond-mat/ K. G. Wilson, The renormalization group and critical phenomena, Rev. Mod. Phys. 55, (1983). U. Schollwöck, The density-matrix renormalization group, Rev. Mod. Phys. 77, 259 (2005). W.H. Press, S.A. Teukolsky, W.T. Vetterling, and B.P. Flannery, Numerical Recipes in C [or C++], 2nd ed., Cambridge University Press (1993).
3 Outline I. Exact Diagonalization (i) Introduction to interacting quantum systems (ii) Representation of many-body states (iii) Complete Diagonalization (iv) Iterative Diagonalization (Lanczos and Davidson) (v) Dynamics (vi) Finite temperature II. Numerical Renormalization Group (i) Anderson and Kondo problems (ii) Numerical RG for the Kondo problem (iii) Numerical RG for quantum lattice problems (iv) Numerical RG for a noninteracting particle III. From the NRG to the Density Matrix Renormalization Group (i) Better methods for the noninteracting particle (ii) Density Matrix Projection for interacting systems (iii) DMRG Algorithms (iv) DMRG-like algorithm for the noninteracting particle IV. The DMRG in Detail (i) Programming details (ii) Measurements (iii) Wavefunction transformations (iv) Extensions to higher dimension
4 Outline I. Exact Diagonalization (i) Introduction to interacting quantum systems (ii) Representation of many-body states (iii) Complete Diagonalization (iv) Iterative Diagonalization (Lanczos and Davidson) (v) Dynamics (vi) Finite temperature II. Numerical Renormalization Group (i) Anderson and Kondo problems (ii) Numerical RG for the Kondo problem (iii) Numerical RG for quantum lattice problems (iv) Numerical RG for a noninteracting particle III. From the NRG to the Density Matrix Renormalization Group (i) Better methods for the noninteracting particle (ii) Density Matrix Projection for interacting systems (iii) DMRG Algorithms (iv) DMRG-like algorithm for the noninteracting particle IV. The DMRG in Detail (i) Programming details (ii) Measurements (iii) Wavefunction transformations (iv) Extensions to higher dimension
5 Outline I. Exact Diagonalization (i) Introduction to interacting quantum systems (ii) Representation of many-body states (iii) Complete Diagonalization (iv) Iterative Diagonalization (Lanczos and Davidson) (v) Dynamics (vi) Finite temperature II. Numerical Renormalization Group (i) Anderson and Kondo problems (ii) Numerical RG for the Kondo problem (iii) Numerical RG for quantum lattice problems (iv) Numerical RG for a noninteracting particle III. From the NRG to the Density Matrix Renormalization Group (i) Better methods for the noninteracting particle (ii) Density Matrix Projection for interacting systems (iii) DMRG Algorithms (iv) DMRG-like algorithm for the noninteracting particle IV. The DMRG in Detail (i) Programming details (ii) Measurements (iii) Wavefunction transformations (iv) Extensions to higher dimension
6 Outline I. Exact Diagonalization (i) Introduction to interacting quantum systems (ii) Representation of many-body states (iii) Complete Diagonalization (iv) Iterative Diagonalization (Lanczos and Davidson) (v) Dynamics (vi) Finite temperature II. Numerical Renormalization Group (i) Anderson and Kondo problems (ii) Numerical RG for the Kondo problem (iii) Numerical RG for quantum lattice problems (iv) Numerical RG for a noninteracting particle III. From the NRG to the Density Matrix Renormalization Group (i) Better methods for the noninteracting particle (ii) Density Matrix Projection for interacting systems (iii) DMRG Algorithms (iv) DMRG-like algorithm for the noninteracting particle IV. The DMRG in Detail (i) Programming details (ii) Measurements (iii) Wavefunction transformations (iv) Extensions to higher dimension
7 Outline V. Recent Developments in the DMRG (i) Classical transfer matrices (ii) Finite temperature (iii) Dynamics (iv) Quantum chemistry (v) Time evolution (vi) Matrix product states (vii) Quantum information VI. Quantum Monte Carlo Methods (i) Review of classical Monte Carlo (ii) Variational and Green s function QMC (iii) World-line QMC (iv) Determinantal QMC (v) Loop Algorithm (vi) Stochastic Series Expansion (SSE) (vii) Diagrammatic QMC
8 Outline V. Recent Developments in the DMRG (i) Classical transfer matrices (ii) Finite temperature (iii) Dynamics (iv) Quantum chemistry (v) Time evolution (vi) Matrix product states (vii) Quantum information VI. Quantum Monte Carlo Methods (i) Review of classical Monte Carlo (ii) Variational and Green s function QMC (iii) World-line QMC (iv) Determinantal QMC (v) Loop Algorithm (vi) Stochastic Series Expansion (SSE) (vii) Diagrammatic QMC
9 I. Exact Diagonalization Direct diagonalization of Hamiltonian matrix on finite clusters ladder Goals ground state properties low-lying excitations dynamics, finite T,... y x two dimensionnal (2D) chain (1D) Advantages almost any system can be treated almost any observable can be calculated quantum-number resolved quantities numerically exact (for finite cluster) Limitation: exponential in lattice size
10 Largest sizes reached S = 1/2 spin models square lattice: N = 40 triangular lattice: N = 39, star lattice: N = 42 maximum dimension of basis: 1.5 billion t-j models checkerboard lattice with 2 holes: N = 32 square lattice with 2 holes: N = 32 maximum dimension of basis: 2.8 billion Hubbard models square lattice at half filling: N = 20 quantum dot structure: N = 20 maximum dimension of basis: 3 billion Holstein models chain with N = 14 + phonon pseudo-sites maximum dimension of basis: 30 billion
11 I (i) Interacting Quantum Systems Here: discrete, finite case system of N quantum mechanical subsystems, l = 1,..., N finite number of basis states per subsystem α l, α l = 1,..., s l more general case: s l (continuum or thermodynamic limit) Properties: N (thermodynamic limit) l x Basis direct product of component basis (continuous quantum field) α 1, α 2,..., α N α 1 α 2... α N total number of states: N l=1 s l arbitrary state in this basis ψ = {αl} ψ(α 1, α 2,..., α N ) α 1, α 2,..., α N behavior governed by Schrödinger equation H Ψ(t) = i t Ψ(t) or H ψ = E ψ (time-independent)
12 Hamiltonians In general, Hamiltonians can connect arbitrary numbers of subsystems H = H (1) l + H (2) lm H (4) lmnp +... l l,m l,m,p H (1) l usually determines α l H (2) lm, sometimes H(4) lmnp will be important here H (2) lm often short-ranged Typical terms: tight-binding term: V(x) H tb = ψ( x) l,m,σ t lm c l,σ c m,σ localized Wannier orbitals (unfilled d- or f- orbitals in transition metals) states 0,,, per orbital 4 N degrees of freedom overlap between near orbitals hopping t lm short ranged (n.n., possibly n.n.n.)
13 local (Anderson) disorder H A l = σ λ l n l,σ, n l,σ c l,σ c l,σ Coulomb interaction between electrons H C lm = screening leads to on-site (Hubbard) interaction e 2 r l r m Hl U = U n l, n l, near-neighbor Coulomb interaction Hlm V = V n l n l+ˆr, (n l σ n l,σ) etc. Spin models S i localized quantum mechanical spins (S = 1/2, 1, 3/2,...) states S S S (2S + 1) N degrees of freedom Heisenberg exchange Hlm Heis = J S l S m = J z Sl z Sm z J xy ( S + l S m + S ) l S+ m strong coupling limit of the Hubbard model at n = 1 (S = 1/2) AF exchange J = 4t2 U variations: J z J xy (Ising or XY anisotropy), Hl n H bq lm = J 2 (S l S m ) 2 (biquadratic) = D(Sz l )2 (single-ion),
14 t J model: strong-coupling limit of doped Hubbard Hlm tj = P Hlm tb P + J (S l S m 14 ) n l n m double occupancy projected out (P) - 3 states/site Anderson impurity - hybridized d (or f) orbital with on-site interaction ) Hl AI = ε d n d l + V (d l,σc l,σ + H.c. + Un d l, n d l, single impurity or lattice (PAM) possible Kondo impurity - localized d spin S Hl K = J ( ) K 2 S l c l,α σ α,β c l,β limit of symmetric Anderson impurity at strong U
15 Lattices square lattice Kagomé lattice Described by unit cell Bravais lattice: translation vectors T 1, T 2 (2D) T 2 T 1
16 finite lattices: finite multiples of T 1, T 2 and boundary conditions periodic, antiperiodic open lattice symmetries: translation multiples of Bravais lattice vector + periodic (AP) BCs rotations e.g., π/2 for a square lattice (group C 4v ) reflection about symmetry axis Tilted clusters 40-site cluster, square lattice (a = 1) T 1 = (1, 0), T 2 = (0, 1) Spanning vectors: F 1 = (6, 2), F 2 = ( 2, 6) In general, F 1 = (n, m), F 2 = ( m, n) N = n 2 + m 2 translational symmetry satisfied reflection/rotation symmetries become more complicated
17 I (ii) Representation of Many-Body States mapping to (binary) integers: spin-1/2 Heisenberg: spin flip = bit flip Hubbard with N σ = {0, 1} N 1 N N 1 1 N N l N l N l N l or N e l S z l other models (t J, S = 1 Heisenberg,...) more complicated Symmetries: given group G with generators {g p } [H, g p ] = 0 H block diagonal (Hilbert space can be divided) Continuous conservation of particle number, S z U(1) permutations of bits total spin SU(2) difficult to combine with space group spin inversion (Z 2 ) can be used Space group translation: abelian local states point group (reflections and rotations): non-abelian in general form symmetrized linear combination of local states
18 Example Reduction of Hilbert space for S = 1/2 Heisenberg on cluster full Hilbert space: constrain to S z = 0: using spin inversion: utilizing all 40 translations: using all 4 rotations: dim= 2 40 = dim= dim= dim= dim= 430, 909, 650
19 I (iii) Complete Diagonalization To solve H ψ = E ψ (H real, symmetric) Method (Numerical Recipes, Ch. 11) 1. Householder transformation - reduction to tridiagonal form T 2n 3 /3 operations (4n 3 /3 with eigenvectors) 2. Diagonalization of a tridiagonal matrix roots of secular equation: inefficient QL (QR) algorithm - factorization T = Q L, Q orthogonal, L lower triangular 30n 2 operations ( 3n 3 with eigenvectors) Useful for: Simple problems, testing Matrix H dense Many eigenstates required But H must be stored entire matrix must be diagonalized
20 I (iv) Iterative Diagonalization Idea: project H onto a cleverly chosen subspace of dimension M N good convergence of extremal eigenstates Methods Power method v n = H n v 0 conceptually simple, but converges poorly needs only two vectors, v n & v n 1 Lanczos: orthogonal vectors in Krylov subspace (spanned by { v n }) simple to implement memory efficient - only 3 vectors needed at once works well for sparse, short-range H Davidson: subspace expanded by diagonal approximation to inverse iteration higher-order convergence than Lanczos (usually) implementation more complicated works best for diagonally-dominated H Jacobi-Davidson: generalization of Davidson nontrivial problem-specific preconditioner (approximation to inverse) can be applied to generalized eigenvalue problem A x = λ B x (A, B general, complex matrices)
21 Lanczos Algorithm 0) choose u 0 (random vector, ψ 0 from last iteration,...) 1) form u n+1 = H u n a n u n b 2 n u n 1 where a n = u n H u n u n u n and b 2 n = u n u n u n 1 u n 1 2) Is u n+1 u n+1 < ε? yes: do 4) then stop no: continue 3) repeat starting with 1) until n = M (maximum dimension) 4) diagonalize u i H u j (tridiagonal) using QL algorithm diagonal elements D = (a 0, a 1,..., a n ), off-diagonal elements O = (b 1, b 2,..., b n ) eigenvalue Ẽ0, eigenvector ψ 0 5) repeat starting with 0), setting u 0 = ψ 0
22 Convergence of Lanczos Algorithm eigenvalues converge starting with extremal ones excited states can get stuck for a while at longer times: true eigenvalus converged spurious or ghost eigenvalues produced multiplicity of eigenstates increases
23 Example: 2D t-j Model Binding of 2 holes R: average hole-hole distance B : binding energy (Poiblanc, Riera, & Dagotto, 1993) holes closer than two lattice spacings pair binding for J > J c, but large finite-size effects Does binding persist for larger lattices and constant doping (more holes)?
24 I (v) Dynamics with Exact Diagonalization Time-dependent correlation functions C(t) = i ψ 0 A(t) A (0) ψ 0 Fourier transform to frequency space (retarded) C(ω + iη) = ψ 0 A (ω + iη H + E 0 ) 1 A ψ 0 (resolvent) Spectral function I(ω) = 1 π lim η 0 + Im C(ω + iη) Examples from theory and experiment name notation operators experiment single-particle spectral weight A(k, ω) A = c k,σ photoemission structure factor S zz (q, ω) A = Sq z neutron scattering optical conductivity σ xx (ω) A = j x optics 4-spin correlation R(ω) k R k S k S k Raman scattering
25 Methods Krylov space method (continued fraction) restart Lanzcos procedure with 1 u 0 = ψ0 A A ψ 0 A ψ 0 In this Lanczos basis, Interpretation: C(z = ω + iη + E 0 ) = ψ 0 A A ψ 0 z a 1 b 2 2 z a 2 b2 3 z a 3... calculation of eigenvector not needed consider Lehmann representation of spectral function I(ω) = n ψ n A ψ 0 2 δ(ω E n + E 0 ) poles and weights of C(z) determine I(ω) weight decreases with n truncate after M steps spectrum discrete finite broadening η
26 Correction vector method (Soos & Ramasesha, 1984) Calculate vectors directly, then Advantages: φ 0 = A ψ 0, φ 1 = (ω + iη H + E 0 ) 1 φ 0 I(ω) = 1 π Im φ 0 φ 1 spectral weight calculated exactly for a given range nonlinear spectral functions computed by higher order correction vectors can be run in conjunction with Davidson algorithm Disadvantage: system (H z) φ 1 = φ 0 must be solved for each ω desired
27 Example: Dynamics in 2D t-j Model Single-particle spectral weight A(k, ω) at k = k F = (π/2, π/2) for one hole 4 4 lattice (Dagotto, Joynt, Moreo, Bacci, & Gagliano, 1993) Does a single quasiparticle propagate in an antiferromagnet? strongly localized hole with string excitations at J/t = 1.0 quasiparticle peak remains until J/t = 0.4 lump with pseudogap at J/t = 0.2 pseudogap due to finite-size effects at J/t = 0 (symmetric in ω)
28 I (vi) Finite Temperature with Exact Diagonalization To calculate finite-t properties in orthonormal basis n A = 1 N N n Ae βh n, Z = n e βh n, Z n Problem: expensive to calculate for all n Idea: stochastic sampling of Krylov space (Jaklic & Prelovsek, 1994 ) A 1 N s R M e βε(r) m r Ψ (r) m Ψ (r) m A r Z R where Z s s N s R R r r M m m e βε(r) m s over symmetry sectors of dimension N s r m n r Ψ (r) m 2 average over R random starting vectors Ψ(r) 0 Lanczos propagation of starting vectors: Ψ(r) m at step m useful if convergence good when M N s and R N s
29 Properties related to high-t expansion T limit correct high to medium T properties in thermodynamic limit low-temperature limit correct (on finite lattice), up to sampling error reduction of (large) sampling error: (Aichhorn et al., 2003) start with: A = 1 Z N n e βh/2 A e βh/2 n, n twofold insertion of Lanczos basis smaller fluctuations at low T can calculate thermodynamic properties: specific heat, entropy, static susceptibility,... static correlation functions dynamics: A(k, ω), S zz (q, ω), σ xx (ω),...
30 Example: t-j Model at finite T Hole concentration c h (= x) vs. chemical potential shift µ = µ h µ 0 h 2D t-j model, 16, 18, 20 sites, t/j = 0.3, t = 0.4eV 0.3 c h 0.2 T/t= experiment (Jaklic & Prelovsek, 1998) µ (ev) experimental results for LSCO from photoemission shift (Ino et al., 1997) holes only when µ < µ 0 h 1.99t as T 0 compressibility finite no phase separation
31 Optical conductivity compared with various cuprates at intermediate doping Bi 2 Sr 2 CaCu 2 O 8 YBa 2 Cu 3 O 7 n h =3/16 J/t=0.3 σ ρ T/t=0.15 T/t=0.3 T/t=0.5 T/t=1.0 (Jaklic & Prelovsek, 1998) 0.2 La 1.8 Sr 0.2 CuO ω (ev) Cuprates measured at T < 200K, c h somewhat uncertain high-t falloff slower for materials transitions to higher excited states? experimental curves: LCSO, c h x = 0.2 (Uchida et al., 1991) BISCCO, c h 0.23 (Romero et al., 1992) YBCO, c h 0.23 (Battlogg et al., 1994)
32 Discussion: Exact Diagonalization Method conceptually straightforward, numerically exact Iterative diagonalization allows the treatment of surprisingly large matrices Efficient implementation using symmetries useful System sizes nevertheless strongly restricted Extensions to basic method can calculate dynamical correlation functions finite temperature properties Not mentioned here, but also possible: calculation of full time evolution of quantum state with Lanczos Benchmark for other methods, useful when other methods fail
Quantum 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 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 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 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 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 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 informationDT I JAN S S"= = 11111'11 I HtI IBlIIIt g ~II. Report: ONR Grant N J Unclassified: For General Distribution LECTF
' DT I,, Final Report: ONR Grant N00014-91-J-1143 Unclassified: For General Distribution LECTF JAN 2 8 1993 Steven R. White, UC Irvine c Summary Over the last two years, we have used several different
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 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 informationPG5295 Muitos Corpos 1 Electronic Transport in Quantum dots 2 Kondo effect: Intro/theory. 3 Kondo effect in nanostructures
PG5295 Muitos Corpos 1 Electronic Transport in Quantum dots 2 Kondo effect: Intro/theory. 3 Kondo effect in nanostructures Prof. Luis Gregório Dias DFMT PG5295 Muitos Corpos 1 Electronic Transport in Quantum
More informationNUMERICAL METHODS FOR QUANTUM IMPURITY MODELS
NUMERICAL METODS FOR QUANTUM IMPURITY MODELS http://www.staff.science.uu.nl/~mitch003/nrg.html March 2015 Andrew Mitchell, Utrecht University Quantum impurity problems Part 1: Quantum impurity problems
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 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 information(r) 2.0 E N 1.0
The Numerical Renormalization Group Ralf Bulla Institut für Theoretische Physik Universität zu Köln 4.0 3.0 Q=0, S=1/2 Q=1, S=0 Q=1, S=1 E N 2.0 1.0 Contents 1. introduction to basic rg concepts 2. introduction
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 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-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 informationDensity Matrix Renormalization: A Review of the Method and its Applications arxiv:cond-mat/ v1 26 Mar 2003
Density Matrix Renormalization: A Review of the Method and its Applications arxiv:cond-mat/0303557v1 26 Mar 2003 Karen Hallberg Centro Atómico Bariloche and Instituto Balseiro 8400 Bariloche, Argentina
More informationGreen's Function in. Condensed Matter Physics. Wang Huaiyu. Alpha Science International Ltd. SCIENCE PRESS 2 Beijing \S7 Oxford, U.K.
Green's Function in Condensed Matter Physics Wang Huaiyu SCIENCE PRESS 2 Beijing \S7 Oxford, U.K. Alpha Science International Ltd. CONTENTS Part I Green's Functions in Mathematical Physics Chapter 1 Time-Independent
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 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 informationarxiv: v1 [cond-mat.str-el] 22 Jun 2007
Optimized implementation of the Lanczos method for magnetic systems arxiv:0706.3293v1 [cond-mat.str-el] 22 Jun 2007 Jürgen Schnack a, a Universität Bielefeld, Fakultät für Physik, Postfach 100131, D-33501
More informationarxiv:cond-mat/ v1 [cond-mat.str-el] 21 Mar 2006
Non-Fermi-liquid phases in the two-band Hubbard model: Finite-temperature exact diagonalization study of Hund s rule coupling A. Liebsch and T. A. Costi Institut für Festkörperforschung, Forschungszentrum
More informationDynamical properties of strongly correlated electron systems studied by the density-matrix renormalization group (DMRG) Takami Tohyama
Dynamical properties of strongly correlated electron systems studied by the density-matrix renormalization group (DMRG) Takami Tohyama Tokyo University of Science Shigetoshi Sota AICS, RIKEN Outline Density-matrix
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 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 informationCan superconductivity emerge out of a non Fermi liquid.
Can superconductivity emerge out of a non Fermi liquid. Andrey Chubukov University of Wisconsin Washington University, January 29, 2003 Superconductivity Kamerling Onnes, 1911 Ideal diamagnetism High Tc
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 informationQuantum Impurities In and Out of Equilibrium. Natan Andrei
Quantum Impurities In and Out of Equilibrium Natan Andrei HRI 1- Feb 2008 Quantum Impurity Quantum Impurity - a system with a few degrees of freedom interacting with a large (macroscopic) system. Often
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 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 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 informationQuantum Cluster Methods
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
More informationA Twisted Ladder: Relating the Iron Superconductors and the High-Tc Cuprates
A Twisted Ladder: Relating the Iron Superconductors and the High-Tc Cuprates arxiv:0905.1096, To appear in New. J. Phys. Erez Berg 1, Steven A. Kivelson 1, Doug J. Scalapino 2 1 Stanford University, 2
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 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 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 informationAngle-Resolved Two-Photon Photoemission of Mott Insulator
Angle-Resolved Two-Photon Photoemission of Mott Insulator Takami Tohyama Institute for Materials Research (IMR) Tohoku University, Sendai Collaborators IMR: H. Onodera, K. Tsutsui, S. Maekawa H. Onodera
More informationarxiv: v1 [cond-mat.str-el] 1 Jun 2007
Numerical approach to low-doping regime of the t-j model arxiv:0706.0105v1 [cond-mat.str-el] 1 Jun 2007 J. Bonča, 1, 2 S. Maekawa, 3, 4 and T. Tohyama 5 1 Faculty of Mathematics and Physics, University
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 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 informationSpin correlations in conducting and superconducting materials Collin Broholm Johns Hopkins University
Spin correlations in conducting and superconducting materials Collin Broholm Johns Hopkins University Supported by U.S. DoE Basic Energy Sciences, Materials Sciences & Engineering DE-FG02-08ER46544 Overview
More informationDynamical Mean Field Theory and Numerical Renormalization Group at Finite Temperature: Prospects and Challenges
Dynamical Mean Field Theory and Numerical Renormalization Group at Finite Temperature: Prospects and Challenges Frithjof B. Anders Institut für Theoretische Physik Universität Bremen Göttingen, December
More informationHow to model holes doped into a cuprate layer
How to model holes doped into a cuprate layer Mona Berciu University of British Columbia With: George Sawatzky and Bayo Lau Hadi Ebrahimnejad, Mirko Moller, and Clemens Adolphs Stewart Blusson Institute
More informationTuning order in cuprate superconductors
Tuning order in cuprate superconductors arxiv:cond-mat/0201401 v1 23 Jan 2002 Subir Sachdev 1 and Shou-Cheng Zhang 2 1 Department of Physics, Yale University, P.O. Box 208120, New Haven, CT 06520-8120,
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 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 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 informationSusumu YAMADA 1,3 Toshiyuki IMAMURA 2,3, Masahiko MACHIDA 1,3
Dynamical Variation of Eigenvalue Problems in Density-Matrix Renormalization-Group Code PP12, Feb. 15, 2012 1 Center for Computational Science and e-systems, Japan Atomic Energy Agency 2 The University
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 informationVisualization of atomic-scale phenomena in superconductors
Visualization of atomic-scale phenomena in superconductors Andreas Kreisel, Brian Andersen Niels Bohr Institute, University of Copenhagen, 2100 København, Denmark Peayush Choubey, Peter Hirschfeld Department
More informationA theoretical study of the single-molecule transistor
A theoretical study of the single-molecule transistor B. C. Friesen Department of Physics, Oklahoma Baptist University, Shawnee, OK 74804 J. K. Ingersent Department of Physics, University of Florida, Gainesville,
More informationFinite-temperature properties of the two-dimensional Kondo lattice model
PHYSICAL REVIEW B VOLUME 61, NUMBER 4 15 JANUARY 2000-II Finite-temperature properties of the two-dimensional Kondo lattice model K. Haule J. Stefan Institute, Jamova 39, 1000 Ljubljana, Slovenia J. Bonča
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 informationCluster Functional Renormalization Group
Cluster Functional Renormalization Group Johannes Reuther Free University Berlin Helmholtz-Center Berlin California Institute of Technology (Pasadena) Lefkada, September 26, 2014 Johannes Reuther Cluster
More informationAn introduction to the dynamical mean-field theory. L. V. Pourovskii
An introduction to the dynamical mean-field theory L. V. Pourovskii Nordita school on Photon-Matter interaction, Stockholm, 06.10.2016 OUTLINE The standard density-functional-theory (DFT) framework An
More informationDimerized & frustrated spin chains. Application to copper-germanate
Dimerized & frustrated spin chains Application to copper-germanate Outline CuGeO & basic microscopic models Excitation spectrum Confront theory to experiments Doping Spin-Peierls chains A typical S=1/2
More informationSpin-charge separation in doped 2D frustrated quantum magnets p.
0.5 setgray0 0.5 setgray1 Spin-charge separation in doped 2D frustrated quantum magnets Didier Poilblanc Laboratoire de Physique Théorique, UMR5152-CNRS, Toulouse, France Spin-charge separation in doped
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 Mechanics Solutions. λ i λ j v j v j v i v i.
Quantum Mechanics Solutions 1. (a) If H has an orthonormal basis consisting of the eigenvectors { v i } of A with eigenvalues λ i C, then A can be written in terms of its spectral decomposition as A =
More informationHole dynamics in frustrated antiferromagnets: Coexistence of many-body and free-like excitations
Hole dynamics in frustrated antiferromagnets: Coexistence of many-body and free-like excitations Collaborators: Luis O. Manuel Instituto de Física Rosario Rosario, Argentina Adolfo E. Trumper (Rosario)
More informationLuigi Paolasini
Luigi Paolasini paolasini@esrf.fr LECTURE 4: MAGNETIC INTERACTIONS - Dipole vs exchange magnetic interactions. - Direct and indirect exchange interactions. - Anisotropic exchange interactions. - Interplay
More informationThe 4th Windsor Summer School on Condensed Matter Theory Quantum Transport and Dynamics in Nanostructures Great Park, Windsor, UK, August 6-18, 2007
The 4th Windsor Summer School on Condensed Matter Theory Quantum Transport and Dynamics in Nanostructures Great Park, Windsor, UK, August 6-18, 2007 Kondo Effect in Metals and Quantum Dots Jan von Delft
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 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 informationReal-Space Renormalization Group (RSRG) Approach to Quantum Spin Lattice Systems
WDS'11 Proceedings of Contributed Papers, Part III, 49 54, 011. ISBN 978-80-7378-186-6 MATFYZPRESS Real-Space Renormalization Group (RSRG) Approach to Quantum Spin Lattice Systems A. S. Serov and G. V.
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 informationThe Gutzwiller Density Functional Theory
The Gutzwiller Density Functional Theory Jörg Bünemann, BTU Cottbus I) Introduction 1. Model for an H 2 -molecule 2. Transition metals and their compounds II) Gutzwiller variational theory 1. Gutzwiller
More informationarxiv:cond-mat/ v1 [cond-mat.supr-con] 28 May 2003
arxiv:cond-mat/0305637v1 [cond-mat.supr-con] 28 May 2003 The superconducting state in a single CuO 2 layer: Experimental findings and scenario Rushan Han, Wei Guo School of Physics, Peking University,
More informationLuigi Paolasini
Luigi Paolasini paolasini@esrf.fr LECTURE 7: Magnetic excitations - Phase transitions and the Landau mean-field theory. - Heisenberg and Ising models. - Magnetic excitations. External parameter, as for
More informationThe NMR Probe of High-T c Materials
R.E. Walstedt The NMR Probe of High-T c Materials 4y Springer Contents Introduction 1 1.1 The Basic Phenomenology of High-T c. Materials 1 1.2 Carrier Doping and the Master Phase Diagram 2 1.3 The NMR
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 informationFROM NODAL LIQUID TO NODAL INSULATOR
FROM NODAL LIQUID TO NODAL INSULATOR Collaborators: Urs Ledermann and Maurice Rice John Hopkinson (Toronto) GORDON, 2004, Oxford Doped Mott insulator? Mott physics: U Antiferro fluctuations: J SC fluctuations
More informationFRG approach to interacting fermions with partial bosonization: from weak to strong coupling
FRG approach to interacting fermions with partial bosonization: from weak to strong coupling Talk at conference ERG08, Heidelberg, June 30, 2008 Peter Kopietz, Universität Frankfurt collaborators: Lorenz
More informationPhotoemission Studies of Strongly Correlated Systems
Photoemission Studies of Strongly Correlated Systems Peter D. Johnson Physics Dept., Brookhaven National Laboratory JLab March 2005 MgB2 High T c Superconductor - Phase Diagram Fermi Liquid:-Excitations
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 informationThe advent of computer era has opened the possibility to perform large scale
Chapter 2 Density Matrix Renormalization Group 2.1 Introduction The advent of computer era has opened the possibility to perform large scale numerical simulations of the quantum many-body systems and thus
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 informationFaddeev Random Phase Approximation (FRPA) Application to Molecules
Faddeev Random Phase Approximation (FRPA) Application to Molecules Matthias Degroote Center for Molecular Modeling (CMM) Ghent University INT 2011 Spring Program Fermions from Cold Atoms to Neutron Stars:
More informationNUMERICAL METHODS FOR QUANTUM IMPURITY MODELS
NUMERICAL METHODS FOR QUANTUM IMPURITY MODELS http://www.staff.science.uu.nl/~mitch003/nrg.html March 2015 Anrew Mitchell Utrecht University Quantum impurity problems Part 1: Quantum impurity problems
More informationThe Hubbard model in cold atoms and in the high-tc cuprates
The Hubbard model in cold atoms and in the high-tc cuprates Daniel E. Sheehy Aspen, June 2009 Sheehy@LSU.EDU What are the key outstanding problems from condensed matter physics which ultracold atoms and
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 informationComputational Nanoscience
Computational Nanoscience Applications for Molecules, Clusters, and Solids KALMÄN VARGA AND JOSEPH A. DRISCOLL Vanderbilt University, Tennessee Щ CAMBRIDGE HP UNIVERSITY PRESS Preface Part I One-dimensional
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 informationNumerical Studies of Disordered Tight-Binding Hamiltonians
Numerical Studies of Disordered Tight-Binding Hamiltonians R.T. Scalettar Physics Department, University of California, Davis, CA 95616 Abstract. These are notes used for a set of lectures delivered at
More informationElectronic quasiparticles and competing orders in the cuprate superconductors
Electronic quasiparticles and competing orders in the cuprate superconductors Andrea Pelissetto Rome Subir Sachdev Ettore Vicari Pisa Yejin Huh Harvard Harvard Gapless nodal quasiparticles in d-wave superconductors
More informationMetal-insulator transition with Gutzwiller-Jastrow wave functions
Metal-insulator transition with Gutzwiller-Jastrow wave functions Odenwald, July 20, 2010 Andrea Di Ciolo Institut für Theoretische Physik, Goethe Universität Frankfurt Metal-insulator transition withgutzwiller-jastrow
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 informationQuantum Monte Carlo Simulations of the Half-filled Hubbard Model. Anders F. J. Gabrielsson
Quantum Monte Carlo Simulations of the Half-filled Hubbard Model Anders F. J. Gabrielsson June 2011 Abstract A Quantum Monte Carlo method of calculating operator expectation values for the ground state
More informationFerromagnetism in an orbitally degenerate Hubbard model
PHYSICAL REVIEW B VOLUME 55, NUMBER 22 Ferromagnetism in an orbitally degenerate Hubbard model J. Kuei * and R. T. Scalettar Department of Physics, University of California, Davis, California 95616 Received
More informationA New look at the Pseudogap Phase in the Cuprates.
A New look at the Pseudogap Phase in the Cuprates. Patrick Lee MIT Common themes: 1. Competing order. 2. superconducting fluctuations. 3. Spin gap: RVB. What is the elephant? My answer: All of the above!
More informationImpurity corrections to the thermodynamics in spin chains using a transfer-matrix DMRG method
PHYSICAL REVIEW B VOLUME 59, NUMBER 9 1 MARCH 1999-I Impurity corrections to the thermodynamics in spin chains using a transfer-matrix DMRG method Stefan Rommer and Sebastian Eggert Institute of Theoretical
More informationX. Zotos - Research Publications
X. Zotos - Research Publications After 2004 1. Phonon-Magnon Interaction in Low Dimensional Quantum Magnets Observed by Dynamic Heat Transport Measurements, M. Montagnese, M. Otter, X. Zotos et al., Physical
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 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 informationWORLD SCIENTIFIC (2014)
WORLD SCIENTIFIC (2014) LIST OF PROBLEMS Chapter 1: Magnetism of Free Electrons and Atoms 1. Orbital and spin moments of an electron: Using the theory of angular momentum, calculate the orbital
More informationEntanglement spectra in the NRG
PRB 84, 125130 (2011) Entanglement spectra in the NRG Andreas Weichselbaum Ludwig Maximilians Universität, München Arnold Sommerfeld Center (ASC) Acknowledgement Jan von Delft (LMU) Theo Costi (Jülich)
More informationarxiv: v1 [cond-mat.str-el] 17 Jan 2011
Computational Studies of Quantum Spin Systems arxiv:1101.3281v1 [cond-mat.str-el] 17 Jan 2011 Anders W. Sandvik Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts
More informationLast Time. Social Network Graphs Betweenness. Graph Laplacian. Girvan-Newman Algorithm. Spectral Bisection
Eigenvalue Problems Last Time Social Network Graphs Betweenness Girvan-Newman Algorithm Graph Laplacian Spectral Bisection λ 2, w 2 Today Small deviation into eigenvalue problems Formulation Standard eigenvalue
More informationQuantum s=1/2 antiferromagnet on the Bethe lattice at percolation I. Low-energy states, DMRG, and diagnostics
Quantum s=1/2 antiferromagnet on the Bethe lattice at percolation I. Low-energy states, DMRG, and diagnostics Hitesh J. Changlani, Shivam Ghosh, Sumiran Pujari, Christopher L. Henley Laboratory of Atomic
More information