February 15, 2015 Kalani Hettiarachchi Collaborators: Valy Rousseau Ka-Ming Tam Juana Moreno Mark Jarrell
Cold Atoms Ø On Surface of Sun: Miss many aspects of nature Ø Surface of Earth: Different states of matter: solid, liquid and gas Ø Further Cooling: In Kelvin : Superconductivity (1911) Superfluidity in Liquid Helium ( 4 He) (1938) In milli-kelvin: Superfluidity in 3 He (1972) Ø Laser cooling/evaporative cooling(1980): In nano-kelvin: Bose-Einstein Condensation (1995) (This was predicted for photons in 1924 by Bose and Einstein) These achievements recognized with Nobel prizes! 2
Optical Lattices and Bose Hubbard Model Ø Optical lattices: Artificial crystal of lights consisting of thousands of optical micro traps: created by interfering of laser beams! Ø Bosons loaded in optical lattices undergo a superfluid-mott insulator transition: Theoretically proposed and verified experimentally! Bose- Hubbard Model D. Jaksch et al., PRL(1998) M. P. A. Fisher et al., PRB(1989) M. Greiner et al., NATURE(2002) Ø The advance of opical lajce experiments and the experimental tunability of Hamiltonian parameters using laser and magneic fields allow the realizaion of strongly correlated 3 models.
Bose Hubbard model μ/u N=3 MoV t, U <<ħ ω t U N=2 N=1 MoV MoV Superfluid H = t (a + i a j +hc)+ U 2 ij i n i (n i 1) ; n i = a i + a i t U U<<t SF The superfluid phase is characterized by Bose- Einstein condensaion, gapless excitaions and finite compressibility in the weak coupling regime. U>>t MI The MoV insulator phase is characterized by commensurate occupaions, gapped excitaions and incompressibility in the strong coupling regime. Mixtures of atoms can give rise to even more interes?ng and complex phases. 8
H = t (a + a +hc) t (b + b +hc) +U n n a i j b i j ab i,a i,b ij ij i Superfluid Mott Insulator B. Capogrosso-Sansone et al, PRA (2010) S. G. Soyler et al, New J. Phys (2009) Commutation rules on different sites! Anti-commutation Rules on identical site Our focus is to study the model as a function of doping. 1. Equal populations (N a =N b ) and mass imbalance (t a <t b ). 2. Population imbalance (N a N b ) and mass imbalance (t a <t b ). Method: QMC simulations using the Stochastic Green function 5 algorithm V.G. Rousseau, PRE. 78, 056707 (2008), and arxiv:1209.0946
t a = 0.08t t b =1.0t U ab = 6t SF? SF ρ = N L 2 N = N a + N b AF AF/SF N a = N b p,q ρ SF = W 2 4tβ S a,b (k) = e ik( p q) n p a,b n q a,b 6
0 5 10 15 x 0 5 10 15 x OBC a 1.0 ρ a = 0.7 ρ b = 0.7 PBC 15 10 5 0 y 15 10 5 0 y n i a = 0 n i a =1 0 5 10 15 n a i =1 n a i = 0 MI a x 0.8 0.6 15 10 5 0 y 0.4 0.2 0 5 10 15 0.0 15 10 5 0 y n a i n b i 0.4 n b i 0.6 n b i =1 n a i n b i 1 n b i 0.6 n b i =1 7 MI b +SF b x b
Finite Temperature Results 8
Phase Diagram of the doped system MI a + MI b + SF b Ferro/PS δ = N a + N b 2 * L 2 1 2 9
T 0.30 NL 0.25 0.20 0.15 SF 0.10 0.05 AF AF/SF 0.00 0.0 0.1 Ø Model is interesting and it has Competing phases. Ferro PS 0.2 0.3 0.4 δ 0.5 Ø Found five distinct phases! AF phase at half filling and low temperature AF to SF phase separation near half filling NL phase at finite temperature SF phase for a broad region of filling Ferro-magnetic PS phase inside SF phase at low temperature Ø PS dome inside SF phase shows greater entropy than AF and similar to the SF phase. Experimental realization is possible! Ø Reference: Complex phases in a two species bosonic Hubbard Model" K. Hettiarachchilage, V. G. Rousseau, Ka-Ming Tam, M. 10 Jarrell, and J. Moreno, Phys. Rev. B 88, 161101(R) (2013).
Ferromagnetism with population imbalance of the species H = t (a + a +hc) t (b + b +hc) +U n n a i j b i j ab i,a i,b ij Ø The prominent experimental challenge of quantum magnetic phases is to reach such low temperatures and entropies needed to observe these phases. ij i Ø The relatively high global entropy of the phase-separated ferromagnetic phase suggests that this ferromagnet should be easier to access experimentally. Ø Then we thought to explore the extent of the phase-separated ferromagnetic phase as a function of finite polarization. 11
Extent of ferromagnetic phase with finite polarization L=10 t a =0.08 t b =1.0 U ab =6 Beta=60 AF AF/SF 12
Phase diagram on the optimal superfluid line L=10 ρ=1.25 P=0.25 βt=80 13
Ferromagnetic susceptibility P=0 ρ=1.4 T c =0.145t 14
Superfluid density of light species 15
Superfluid density of heavy species 16
k-space momentum distribution L=10 βt=80 U=6t 17
High Global Entropy for Ferromagnetic phase 18
Ø Population imbalance between species extends the region of ferromagnetism. Ø We found optimal superfluid line of fully phase separated ferromagnet with average local densities as n a ~0 and n b ~1 in one region and n a ~1 and n b ~0.5 on the other. Ø Ferromagnetic phase shows high global entropy Experimental realization is possible! Ø Reference: Ferromagnetic Phase Separated Region in the Polarized Two-species Bose Hubbard Model" K. Hettiarachchilage, V. G. Rousseau, Ka-Ming Tam, M. Jarrell, and J. Moreno, Phys. Rev. B 90, 205104 (2014). 19
Ferromagnetism with a harmonic trap Ø The realization of magnetic phases with the presence of confinement is of interest in experiments involving cold atoms. Ø We want to look for the existence of the ferromagnetic phase shown at finite doping and polarization with the presence of a harmonic trapping potential. Ø We look for the ferromagnetic phase in both hard-core and soft-core limits. H = t a (a + i a j +hc) t b (b + i b j +hc) +U ab n i,a n i,b +W i L 2 + U aa 2 i ij n i,a ( n i,a 1) + U bb 2 ij i n i,b ( n i,b 1) i i ( ) 2 (ni,a + n i,b ) 20
Hard-core limit a L=40, W=0.06, N a =450, N b =675, U ab =6t b 30 1.0 0.8 30 20 0.6 0.4 20 10 0.2 0.0 10 0 0 0 10 20 30 0 10 20 30 21
Soft-core limit a L=20, W=0.2, N a =100, N b =150, U ab =6,U aa =2.5, U bb =1.75 b 15 10 5 4 3 2 1 0 15 10 5 0 0 5 10 15 0 0 5 10 15 22
Findings Ø Such trapping effects can lead to have the ferromagnetic phase in both hard-core and soft-core bosons. Ø In hard-core limit, we observe the phase separation where presumably the heavy species is localized in the center of the trap, and surrounded by the light species. Ø In soft-core limit, we observe that the light species gets trapped in the center of the lattice, and is surrounded by the heavy species. Ø The realization of the model is an open question for experimentalists. In collaborating with an experimental group, we are trying to figure out the precise range of parameters the experiments can reach. Thank you! 23