Bose-Hubbard model describes the hopping of bosonic particles in the presence of (repulsive, and usually on-site) interaction. Depending on the ratio between the tunneling term and the interaction term , the macroscopic properties of the system may change from superfluid to Mott insulator.
The Bose-Hubbard model is described by the following Hamiltonian:
where denotes the sum over nearest-neighbor sites, and with bosonic commutation relations, e.g. .
Physical Realization of the Bose-Hubbard model
The Bose-Hubbard model is nowadays most commonly realized by neutral atomic gases trapped in a optical lattice.
Atom trapped in a Optical lattice
We shall consider a neutral atom interacting with a spatially varying laser field. Let us for the moment assume the frequency of the laser field is comparable but detuned from the transition between an electronic excited state and the ground state. Then, from second-order perturbation theory, the time-average energy shift is proportional to the field intensity , but inversely proportional to the detuning ,
Therefore, an atom will see an effective potential which is proportional to the intensity and is attractive when , repulsive when .
A regular "lattice" can be constructed by the interference of pairs of (counter-propagating) polarized light beams. For example, and . The total electric field is simply . We therefore have a spatially varying potential of the form (apart from a constant and consider 1D):
where is a controllable constant. It is straight forward to generalize the analysis into 3D by including two extra pairs of counter propagating beams. Here we note that in real experiments, this potential is superposed by a much slowly varying trapping potential to avoid loss of atoms.
In principle the width (wavelength of the laser beams) of the wells, and the height (light intensity) of the wells are adjustable. The range of interest for constructing the Bose-Hubbard Hamiltonian is the range where the wavefunction of the system may be described by the localized ground-state wavefunction for each well, and the correction due to the overlap of the wavefunctions nearby. This regime is also known as the tight-binding approximation.
Under this approximation, the hopping (non-interacting) part of the Hamiltonian can be described by
where the tunneling term depends on the overlap of the localized wavefunctions between two wells, and is assumed to be the same throughout the lattice.
For a dilute gas, the effective two-body (s-wave) interaction is the dominate interaction, and is given by
where is the scattering length. The wavefunction overlap between different wells are exponentially small, compared with the on-site part. We therefore keep only the part where two bosons are in the same well. Provided that the interaction is "weak", compared with the excitation energy of the well (which we shall justify below), the second quantized form of the interaction term is
Constraint for the scattering length
We shall now consider the condition where higher excitations of the localized states can be neglected. First, we assume each localized well can be approximated by a harmonic potential with energy , which depends on the curvature of the bottom part of the well, and the atomic mass. Let us define the length scale (called zero-point length) associated with this energy scale by
The ground state wavefunction of a harmonic oscillator is given by . Consider two particles in the same state, with ,
apart from some constant. To avoid the excited states, we impose a constraint , which implies that
This is the result we want in this section.