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 Construction of Decorated Lattices Superconductivity on the Lieb Lattice The 1/5th depleted Square Lattice: Magnetism in CaV 4 O 9 Conclusion Funding: NSF (Physics at the Information Frontier) DARPA (Optical Lattice Emulation Program) XXXVIII Meeting of the Brazilian Physical Society, May 28, 215
. Cast of Characters Vlad Benoit Ehsan George Fred Iglovikov Gremaud Khatami Batrouni Hebert
(Single band) Hubbard Hamiltonian Ĥ = t ij σ (c iσ c jσ +c jσ c iσ )+U i (n i 1 2 )(n i 1 2 ) µ iσ (n iσ +n iσ ) Two spin species σ =,. Kinetic energy t describes hopping between near-neighbor sites ij. On-site repulsion U discourages double occupancy Chemical potential µ controls filling. Half-filling (ρ = 1) at µ =. Cuprate materials (LaSrCuO, YBaCuO,...) drive interest in 2D square lattice: Cu atoms in CuO 2 sheets are in that geometry. Ignore bridging O atoms. Ignore La, Sr, Y, Ba in between layers.
Mott Insulating Behavior in Hubbard Model Half-filled lattice : Average number of fermions is one per site. Kinetic energy and entropy both favor particles moving around lattice. Metal: odd number (one) particle per cell/site. (Like cuprate superconductors.)
Large repulsive interaction U between fermions on the same site. +U Mott Insulator forms. Basic physics of parent compounds of cuprate superconductors! Also transition metal monoxides: FeO, CoO, MnO.
Two ways to destroy Mott Insulator: Decrease U/t: By applying pressure (MnO) Shift n 1: Dope chemically (cuprate superconductors) What is optimal spin arrangement? Hopping of neighboring parallel spins forbidden by Pauli. Antiparallel arrangement lower in second order perturbation theory. x t t E (2) = E (2) t 2 /U = J Mott insulating behavior and antiferromagnetism go hand-in-hand. Qualitative picture of cuprate physics before doping. Still do not really understand why cuprates become superconducting after doping.
Special Features of square lattice hopping Hamiltonian (U = ).4 4 2.3 N(E).2 k y (π,π).1-2 -4-2 2 4 E -4-4 -2 2 4 k x Left: Van Hove singularity of density of states at E = (half-filling) Right: Nesting : At ρ = 1, wavevector (π, π) connects big sections of Fermi surface. Favors ordering at (π, π) like antiferromagnetism.
Determinant Quantum Monte Carlo Map Hamiltonian with electron-electron interactions electrons coupled to space-(imaginary) time dependent field S iτ. Trace out fermions analytically determinants detm σ (S iτ ). Dimension of M is number of sites in lattice. (Classical) Monte Carlo sampling of S iτ Algorithm is order N 3 L. L is number of (imaginary) time slices. N 1 2 1 3 lattice sites/electrons L 1 2 to reach low T (large β). Measurements (any single or two particle Green s function): c iσ c jσ [M 1 σ ] ij = [G σ ] ij Sign Problem: At low temperature detm σ can go negative. At special symmetry points (half-filling ρ = 1) detm = detm. No Sign Problem!
DQMC Results- Density of States.4.3 U=2 U=4 U=6 N(ω=).2.1..1 1 T N(ω = ) vanishes Insulating Mott Gap.
DQMC results- Antiferromagnetic spin correlations c(l) = (n j+l n j+l )(n j n j ) βt = 2 T = t/2 = W/16 (bandwidth W = 8t).1.5 (a) 2 x 2 U = 2. C(lx,ly) -.5 (1,1) -.1 β = 32 β = 2 (,) (1,) β = 12 -.15 (,) (1,) (1,1) (,) Hirsch/Hirsch-Tang: Proof of AFLRO in ground state of half-filled Hubbard Hamiltonian.
U = 2 Fermi function: DQMC results- Fermi distribution n(k x,k y ) ρ =.2 ρ =.4 ρ =.6 ρ =.8 ρ = 1. π π/2 -π/2 -π -π -π/2 π/2 π -π/2 π/2 π -π/2 π/2 π -π/2 π/2 π -π/2 π/2 π 1.8.6.4.2 n(k) U = 2 Gradient of Fermi function: π 2.5 π/2 -π/2 2 1.5 1.5 n(k) -π -π -π/2 π/2 π -π/2 π/2 π -π/2 π/2 π -π/2 π/2 π -π/2 π/2 π
Fermi surface is smeared further by increasing interaction strength. U = 4 Fermi function: π π/2 -π/2 -π ρ =.2 β = 8 -π -π/2 π/2 π ρ =.4 β = 8 -π/2 π/2 π ρ =.6 β = 6 -π/2 π/2 π ρ =.8 β = 4 -π/2 π/2 π ρ = 1. β = 8 -π/2 π/2 π 1.8.6.4.2 n(k) U = 4 Gradient of Fermi function: π π/2 -π/2 1.5 1.5 n(k) -π -π -π/2 π/2 π -π/2 π/2 π -π/2 π/2 π -π/2 π/2 π -π/2 π/2 π
Square lattice is bipartite Divides into A and B sites such that all neighbors of A belong to B and vice-versa. Antiferromagnetism is unfrustrated. Hubbard model has particle-hole symmetry: c iσ c iσ c jσ ( 1) i c iσ ( 1) i+j c iσ c jσ = c jσ c iσ n iσ 1/2 1 n iσ 1/2 = 1/2 n iσ Ĥ Ĥ (at µ = ) Implications for DQMC: No sign problem at half-filling.
Bipartite doesn t mean N A = N B! For example, the Lieb Lattice (CuO 2 sheets of cuprates) Some theorems (Lieb and Sutherland) Near neighbor hopping Hamiltonian has spectrum N A - N B eigenvalues E = (In above, N A = 2N B.) Flat band! 2N B eigenvalues in ± pairs. U > ground state has spin S nonzero (suggested from AF picture!)
Topologically localized states (Sutherland): Provide simple picture of E = modes. 2 3 1 4 Ψ = (c 1 c 2 +c 3 c 4) vac ˆT Ψ = Can do this construction about any missing site. These are the N A N B = N B states with E =. These are the N A N B = N B states with E =. Note: Can view Lieb as 1/4-depleted square lattice.
Construction of Decorated Lattices [1] Begin with any bipartite lattice. [2] Select midpoints of all bonds. [3] Connect any two midpoints that share a vertex. Theorem: The resulting lattice has localized states. 2 1+ 3+ 4 [1] [2] [3] Square plaquette lattice! ˆT (c 1 c 2 +c 3 c 4) vac =.
Lieb Construction: How to get localized states (continued)! [1] Begin with any bipartite lattice. [2] Select midpoints of all bonds. [3] Connect any two midpoints that share a vertex. Theorem: The resulting lattice has localized states. 1+ 6 2 5+ 3+ 4 [1] [2] [3] Kagome Lattice! ˆT (c 1 c 2 +c 3 c 4 +c 5 c 6) vac =.
Superconductivity on the Lieb Lattice Ĥ = t ij σ (c iσ c jσ +c jσ c iσ ) U i (n i 1 2 )(n i 1 2 ) µ iσ (n iσ +n iσ ) (Partial) Particle-Hole Transformation c i c i c iσ c j n i 1/2 ( 1) i c i c i ( 1) i+j c i c j = c j c i 1 n i 1/2 = 1 2 n i n i 1/2 n i 1/2 Ĥ(U > ) Ĥ(U < ) (at µ = ) Repulsive (U > ) Hubbard model: antiferromagnetism. Attractive (U < ) Hubbard model: charge density wave (CDW) and pairing.
a 1 b c 4 2 3 Lower dispersing band < ρ < 2/3 Flat band 2/3 < ρ < 4/3 Upper dispersing band 4/3 < ρ < 2 CDW order at ρ = 2/3 and ρ = 4/3? Pairing (superconductivity) when doped away from these values? Related (?) question: System of bosons in a flat band: no k which minimizes ǫ(k) Would Bose-Einstein condensation still occur? (Stamper-Kurn)
DQMC simulations of Attractive Hubbard Model on Lieb Lattice Occupations of minority ( copper ) and majority ( oxygen ) orbitals: 2. N =18, β =36 n i [t] 1.5 1..5 n d,u = n px,u =..5 n d,u = 4 n px,u = 4 n d,u = 8 n px,u = 8 1. 1.2 1.4 1.6 1.8 2. ρ[t] As total ρ increases past ρ = 2/3 minority orbitals depopulate. (Strange?) Local moment
.25.2 U = 4, β =36 N =27 N =48 N =75 N =18 m 2 [t].15.1.5. 1. 1.2 1.4 1.6 1.8 ρ[t] Constant within the flat band. For square lattice, m 2 falls immediately with doping ρ = 1+δ.
Pairing structure factor P s = (1/N 2 ) ij i j j = c i c i P s [t].2.15.1.5 U = 4..15 mean field N =27 N =48 N =75 N =18 (a).8. (b).8 U = 4 U = 8...2 P s [t].1.5 U = 8. 1. 1.2 1.4 1.6 1.8 ρ[t] P s has minimum at ρ = 4/3 (and 2/3) where CDW competes strongly. U = 4 and half-filling, extrapolate to non-zero order parameter for N. U = 8 and half-filling: no LRO at β = 36. (LRO on square lattice at β = 12). Large U: bosonic limit : No BEC in a flat band. (No k with ǫ(k) smallest.)
The 1/5th depleted Square Lattice: Magnetism in CaV 4 O 9 CaV 2 O 5 : A one-fifth depleted lattice: t t a 1 a 2 V V O3 V V O1 Intra-plaquette t and inter-plaquette t hoppings differ.
t t a 1 a 2 t t t >> t: singlets form on dimers and no long range AF order. t >> t : singlets form on plaquettes and no long range AF order. Troyer: Heisenberg limit determined region of AF order around J J. What about U > Hubbard?
Antiferromagnetic Structure Factor S AF 3 25 2 15 ( 12) U=6, β=2 64 sites 144 sites m(a t )-m(a t ) ( 12) 1 5 Plaquettes.2.4.6.8 1 t /t.8.6.4 t/t Dimers.2 Region with AFLRO in vicinity of t t for U = 6. (Blue bar is Heisenberg limit.)
Anisotropy of NN Spin Correlations σ i σ i+a t - σ i σ i+a t.8.6.4.2 (a) β=2 U=1 U=2 U=3 U=4 U=5 U=6 U=7 U=8 U=12 Difference between near-neighbor spin correlations on intra-plaquette (t) and inter-plaquette (t ) bonds is small in region of AF ordering. High symmetry point (HSP): nn spin correlations on 2 types of bond are equal. (Alternate way of determining quantum critical point qualitatively.)
Quantum Monte Carlo Phase Diagram (Half-Filling)
RPA Phase Diagram (general density) Because of sign problem must do RPA for densities away from half-filling. Q=( π,π) FM Q=(,) AF Q=(,) Stripe Q=( π,) Stripe 6 4 U c 2 1 DOS 2 1 t =t=1-2 2 ω.8.6 ρ.4 Q=(,) FM Q=(π,π) FM Q=(,) AF Q=(π,π) AF Q=(,) Stripe Q=(π,) Stripe.2
Conclusions Hubbard model describes Mott insulators (and antiferromagnetism). Interesting lattices exist where wave functions are localized. These are often associated with flat bands. Putting together exceptional lattices and Hubbard physics. unusual charge order and magnetic order. and unusual superconductivity.