Bose Hubbard model

Introduction

The Bose-Hubbard Model (BHM) was first studied in detail in 1989 by Fisher et al. They studied a three dimensional Bose gas in several different types of systems, including the most simple case at zero temperature. --> This is just an idea, I havent finished it yet, and dont know if it is necessary so let me know what you guys think about it.

Figure 1: Phase Diagram for Zero Temperature Bose-Hubbard Model. Source Fisher et al.

---

Hamiltonian

Solving the Hamiltonian

The single site Hamiltonian for interacting bosons is:

${\mathcal {H}}_{\text{site}}={\frac {U}{2}}b^{\dagger }b\left(b^{\dagger }b-\mathbf {1} \right)-\mu b^{\dagger }b{\text{ }}+\left(\psi ^{*}b+b^{\dagger }\psi \right)$

Where $U>0,\mu \geq 0$ and $\psi$, $\psi ^{*}$ are c-numbers.

$\psi =\psi ^{*}=0$

In this case the Hamiltonian can be rewritten as:

${\mathcal {H}}_{\text{site}}={\frac {U}{2}}b^{\dagger }b\left(b^{\dagger }b-\mathbf {1} \right)-\mu b^{\dagger }b={\frac {U}{2}}b^{\dagger }b{\text{ }}\left(b^{\dagger }b-\left(1+2{\frac {\mu }{U}}\right)\mathbf {1} \right)$.

Unnormalized eigenstates of this Hamiltonian are$|n\rangle$.

${\mathcal {H}}_{\text{site}}|n\rangle ={\frac {U}{2}}b^{\dagger }b{\text{ }}\left(b^{\dagger }b|n\rangle -\left(1+2{\frac {\mu }{U}}\right)|n\rangle \right)={\frac {U}{2}}{\text{ }}\left(n^{2}-\left(1+2{\frac {\mu }{U}}\right)n\right)|n\rangle$

The normalized eigenstate is ${\frac {1}{\sqrt {n!}}}|n\rangle$

Experimental Evidence - $^{87}Rb$ Bose-Einstein Condensate

Experimental Setup

Figure 2: Observed diffraction patterns showing phase transition from superfluid to Mott insulator. Source Greiner et al.

The first experimental evidence in support of the BHM was observed in 2002.

Greiner and collaborators created a $^{87}Rb$ Bose-Einstein condensate that was spherically symmetric, and magnetically trapped. To form the three dimensional lattice potential three optical standing waves were set up perpendicular to eachother with the crossing point of the three beams positioned at the center of the condensate. The magnetically trapped condensate was then transferred to the optical lattice by slowly increasing the beam intensity. By increasing slowly it was ensured that the atoms stayed in the many body ground state. When the beams reached full intensity the condensate was spread over more than 150,000 lattice sites with an average of 2.5 atoms per lattice cite in the center of the condenstate

Adjusting the potential depth via frequency variations allowed for Greiner et al. to essentially set the phase of the condensate. When they would set a potential, they would then turn it off. The atomic wavefunctions then expand freely and interfere with eachother. This allowed for a test of phase coherence of the system via diffraction patterns.

The phase was also tested by measuring the excitation spectrum of the condensate, as the superfluid state and the Mott insulator state should have markedly different spectrum's. In the Mott insulator phase the atoms are well localized in the lattice. An excitation occurs when a particle hole pair is created. For example if we had exactly one particle per lattice site, then an excitation would be when one particle tunnels to an adjacent site, creating a site with no particles and a site with two.

Superfluid Phase

For the superfluid phase it was expected to observe sharp diffraction patterns. This is due to the fact that in the superfluid phase all atoms are delocalized, with equal relative phases between different lattice sites. Therefore the sharp diffraction patterns represent extreme phase coherence in the system.

Mott Insulator Phase

For the Mott insulator phase it was expected to have a blurred diffraction pattern. This is because the system should lose all phase coherence which will result in no observable diffraction pattern. It can be observed in figure 2 X that there is a clear loss of sharp diffraction patterns in box g, this corresponds to the onset of the Mott insulator phase.