1 / 32

Finite-size effects in transport data from Quantum Monte Carlo simulations

Finite-size effects in transport data from Quantum Monte Carlo simulations. Raimundo R. dos Santos Universidade Federal do Rio de Janeiro. Financial support:. Collaborators:. Rubem Mondaini UFRJ. Karim Bouadim Ohio State → ITP-Stuttgart. Thereza Paiva UFRJ. arXiv:1107.0230.

irina
Download Presentation

Finite-size effects in transport data from Quantum Monte Carlo simulations

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

E N D

Presentation Transcript

  1. Finite-size effects in transport data from Quantum Monte Carlo simulations Raimundo R. dos Santos Universidade Federal do Rio de Janeiro Financial support:

  2. Collaborators: Rubem Mondaini UFRJ Karim Bouadim Ohio State → ITP-Stuttgart Thereza Paiva UFRJ arXiv:1107.0230 VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  3. Outline • Motivation: Strongly correlated fermions • The Hubbard model • (Determinant) Quantum MC Simulations • Probing MIT’s: compressibility • Probing MIT’s: dc-conductivity • Probing MIT’s: the Drude weight • A bonus: <sign> and compressibility • Conclusions VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  4. a a a dE dE Ashcroft & Mermin (1976) Ibach & Lüth (2003) Motivation: Strongly correlated fermions Independent electrons in solids: periodic crystalline potential atomic limit (tight-binding): more localized nearly-free electrons: delocalized VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  5. Therefore, independent electron approx’n + band theory explains: metals, insulators, semiconductors... Insulator (eV) or Semiconductor (0.1 eV) Metal Ashcroft & Mermin (1976) Density of states (NFE or TB) VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  6. Therefore, electrons move collectively to minimize energy: they are strongly (due to interactions) correlated Consequence: indepedent electron approx’n can fail seriously; e.g., metallic behaviour for an insulating (Mott) system However, beware of narrow bands (especially d and f): ⇒ greater tendency to localize electrons ⇒ enhances likelihood of two electrons occupying the same site ⇒ Coulomb repulsion can no longer be neglected However, beware of narrow bands (especially d and f): ⇒ greater tendency to localize electrons ⇒ enhances likelihood of two electrons occupying the same site ⇒ Coulomb repulsion can no longer be neglected VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  7. parent compound (undoped) Superconductor Insulator! Superconductivity, MIT, itinerant magnetism, etc.: interactions must be included in a fundamental way Metal! AFM insulator (Pb) http://newscenter.lbl.gov/news-releases/2011/03/24/pseudogap/ High-Tc cuprates band theory including correlations VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  8. Low-dimensional systems/confined geometries: fluctuations can disrupt ordered states Add strong correlations ⇒ need unbiased methods to tackle these problems Quantum Monte Carlo simulations have proved very useful, but... • ...some bottlenecks remain, some of which we address here: • finite temperatures, but usually need extrapol’ns to T=0 • “minus-sign problem” • MIT’s • lack of order parameter • inconsistent results from different probes • lattice finite size ⇔ gaps in the spectrum ⇒ bad for transport properties: “false-positive” insulating states VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  9. John Hubbard (1931-1980) The Hubbard model (or, the “Ising model” for strongly correlated fermions) VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  10. Simplest case: s-orbitals, spin-1/2 fermions ⇒ 4 states per site: Charge and spin d.o.f. The Hubbard model (or, the “Ising model” for strongly correlated fermions) Coulomb repulsion: favours localization chemical pot’l: controls electronic density band energy: favours delocalization Special cases: VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  11. QMC Phase diagram: square lattice at T = 0 Mean field (Hartree-Fock) Hirsch, PRB (1985); Hirsch & Tang, PRL (1989) VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  12. Lz=2 Lz=1 Have a good testing ground to compare probes, approximations, finite-size effects, etc., in the study of MIT’s Useful in other situations with overlying structures (e.g., superlattices) The finite-size dimension Lzmay include several layers: computational effort determined by the # of sites, not by Lz but within FSS theory, the size is as small as Lz VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  13. two-body: V one-body (bilinear ⇒ integrable): K Imaginary-time interval (0,β) discretized intoM = β/Δτ slices; typically Δτ = (8U)-1/2. #sites is nowLd×M (Determinant) Quantum Monte Carlo Simulations Blankenbecler et al., PRD (1982); Hirsch, PRB (1985); White et al., PRB (1989); Loh et al., PRB (1990); dos Santos, BJP (2003) Preparation: Use Trotter formula: Now we can transform two-body term into a bilinear form: HS transf’n VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  14. On every site (i, ) of the space-time lattice, an aux. Ising field si ( ) is introduced with Two-body term: discrete Hubbard-Stratonovich transf’n Hirsch, PRB (1983) → continuous auxiliary-field x Now both K and V are bilinear ⇒ fermionic d.o.f. can be traced out (see given refs. for details), ⇒ remaining Ising d.o.f. are sampled by MC... VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  15. with No guarantee that the “Boltzmann weight” is positive ⇒ leads to the “minus-sign problem” effective “density matrix” ...but the “Boltzmann weight” is quite distinct from the usual case: ➔ Ld × Ld matrices; σ refers to original fermionic channels VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  16. Most “measurable” quantities expressed in terms of these Green’s functions (see given refs. for details) ⇒ simplicity N.B.: Manipulation of Green’s functions → Ns × Ns matrices Average values e.g., single-particle Green’s functions e.g., unequal-time single-particle Green’s functions VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  17. Sampling updates (non-local!) also expressed in terms of these Green’s functions (see given refs. for details) ⇒ simplicity Sampling the HS (Ising) fields VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  18. dos Santos, BJP (2003) The “minus-sign” problem: VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  19. ▲U = 2 ⎯ U = 0 • Conclusion: • “closed-shell” effects arise due to small lattice size; • these give rise to false insulating states; • only dismissed through a sequence of data for different lattice sizes: incompressible densities strongly dependent on size. ∴ Incompressibility(κ = 0) signals an insulating state Probing MIT’s: compressibility Non-interacting case (U = 0) Mondaini et al., (2011) VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  20. Probing MIT’s: dc-conductivity VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  21. The conductivity suffers from the same closed-shell effects as the compressibility. Mondaini et al., (2011) VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  22. A shortcut for calculating the dc-conductivity Randeria et al., PRL (1992); Trivedi and Randeria PRL (1995); Trivedi et al., PRB (1996) hard to establish smallness a priori; use data for σ(ω) to assess reliability. obtainable directly from QMC data (no need to invert Laplace transform!) Mondaini et al., (2011) VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  23. 10 x 10 Results to order zero Half filling: negative slope at high temperatures; insulator only detected at lower T ρ = 0.42: σ → 0 as T → 0 ⇒ insulator? closed-shell effect! ρ = 0.66: σ ↑ as T → 0 ⇒ metal ✓ Mondaini et al., (2011) VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  24. 10 x 10 Mondaini et al., (2011) Results to 2nd order Half filling: positive slope at high temperatures insulator detected at higher T ρ = 0.42: still σ → 0 as T → 0 ⇒ insulator closed-shell effect! (expected: present at full calc’ns) ρ = 0.66: σ ↑ as T → 0 ⇒ metal ✓ σdc(2) cannot be generically considered as a perturbation even at the lowest T’s considered VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  25. Strategy: calculate a finite-temperature D and examine its behaviour as T→ 0: if D → 0 ⇒ Insulator if D → D0 > 0 ⇒ Metal 4x4, T = 0 Lanczos Incoherent response Drude weight Dagotto, RMP (1994) Probing MIT’s: the Drude weight Zero-temperature limit of the frequency-dependent conductivity: D readily available from QMC simulations [Scalapino et al., PRB (1993)]: VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  26. 10 x 10 Mondaini et al., (2011) Results for the Drude weight • Half filling: • steady decrease as T is lowered: • ⇒ insulator ρ = 0.42: D → D0 as T → 0 ⇒ metal ✓ no closed-shell effects! ρ = 0.66: D → D0 as T → 0 ⇒ metal ✓ D weakly T-dependent in metallic state ⇒ reliable extrapolations VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  27. No “hickups” at problematic density ρ=0.42 for L=10 ρ = 0.42 ρ = 0.66 ρ = 1 Mondaini et al., (2011) Size dependence of the Drude weight at finite temperatures half filling: size dependence only appreciable when interaction is on away from half filling: size dependence very weak, even for interacting case (size of time slices does not influence final results) VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  28. Mondaini et al., (2011) A bonus: <sign> and compressibility U = 2 U = 3 U = 4 VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  29. A bonus: <sign> and compressibility Mondaini et al., (2011) VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  30. A bonus: <sign> and compressibility Mondaini et al., (2011) VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  31. Conclusions • Finite-size effects in transport properties aren’t just an accuracy issue. • The Drude weight is the most reliable probe of a MIT. • Beware of closed-shell effects in other probes. • Shortcut to calculate dc-conductivity very dangerous (higher order terms very important). • Minus-sign less harmful when system is incompressible (false insulating). VI BMSP - Cuiabá, Brasil - Aug ’11 - Finite-size effects in transport data from Quantum Monte Carlo simulations - RR dos Santos - UFRJ

  32. Thanks for your attention!

More Related