Spin magnetic moment

The Bose–Hubbard model gives an approximate description of the physics of interacting bosons on a lattice. It is closely related to the Hubbard model which originated in solid-state physics as an approximate description of superconducting systems and the motion of electrons between the atoms of a crystalline solid. The name Bose refers to the fact that the particles in the system are bosonic; the model was first introduced by Gersch H., Knollman G in 1963, The Bose–Hubbard model can be used to study systems such as bosonic atoms on an optical lattice. In contrast, the Hubbard model applies to fermionic particles such as electrons, rather than bosons. Furthermore, it can also be generalized and applied to Bose–Fermi mixtures, in which case the corresponding Hamiltonian is called the Bose–Fermi–Hubbard Hamiltonian.

The hamiltonian

The physics of this model is given by the Bose–Hubbard Hamiltonian:

$H=-t\sum _{⟨i,j⟩}{b}_i^{†}{b}_j+\frac{U}{2}\sum _i{\hat{n}}_i({\hat{n}}_i-1)-\mu \sum _i{\hat{n}}_i$.

Here i is summed over all lattice sites, and $⟨i,j⟩$ denotes summation over all neighboring sites i and j.. $b_i^{†}$ and $b_i^{}$ are bosonic creation and annihilation operators. ${\hat{n}}_i=b_i^{†}{b}_i$ gives the number of particles on site i. Parameter $t$ is the hopping matrix element, signifying mobility of bosons in the lattice. Parameter $U$ describes the on-site interaction, if $U>0$ it describes repulsive interaction, if $U<0$ then the interaction is attractive. $\mu$ is the chemical potential.

The dimension of the Hilbert space of the Bose–Hubbard model grows exponentially with respect to the number of particles N and lattice sites L. It is given by: $D_b=\frac{(N_b+L-1)!}{N_b!(L-1)!}$ while that of Fermi–Hubbard Model is given by: $D_f=\frac{L!}{N_f!(L-N_f)!}.$ The different results stem from different statistics of fermions and bosons. For the Bose–Fermi mixtures, the corresponding Hilbert space of the Bose–Fermi–Hubbard model is simply the tensor product of Hilbert spaces of the bosonic model and the fermionic model.

Phase diagram

At zero temperature, the Bose–Hubbard model (in the absence of disorder) is in either a Mott insulating (MI) state at small $t/U$, or in a superfluid (SF) state at large $t/U$. The Mott insulating phases are characterized by integer boson densities, by the existence of an energy gap for particle-hole excitations, and by zero compressibility. In the presence of disorder, a third, ‘‘Bose glass’’ phase exists. The Bose glass phase is characterized by a finite compressibility, the absence of a gap, and by an infinite superfluid susceptibility. It is insulating despite the presence of a gap, as low tunneling prevents the generation of excitations which, although close in energy, are spatially separated.

Implementation in Optical Lattices

Ultracold atoms in optical lattices are considered a standard realization of the Bose Hubbard model. The ability to tune parameters of the model using simple experimental techniques, lack of lattice dynamics, present in electronic systems provides very good conditions for experimental study of this model.

The hamiltonian in Second quantization formalism describing a gas of ultracold atoms in the optical lattice potential is of the form : $H=\int d^3\overrightarrow{r}{\hat{\psi }}^{†}(\overrightarrow{r})(-\frac{{\hslash }^2}{2m}{\mathrm{\nabla }}^2+V_{latt.}(x))\hat{\psi }(\overrightarrow{r})+\frac{g}{2}{\hat{\psi }}^{†}(\overrightarrow{r}){\hat{\psi }}^{†}(\overrightarrow{r})\hat{\psi }(\overrightarrow{r})\hat{\psi }(\overrightarrow{r})-\mu {\psi }^{†}(\overrightarrow{r})\hat{\psi }(\overrightarrow{r})$ where, $V_{latt}$ is the optical lattice potential, g is interaction amplitude (here contact interaction is assumed), $\mu$ is a chemical potential. Standard tight binding approximation (see this article for details) $\hat{\psi }(\overrightarrow{r})=\sum _i{w}_i^{\alpha }(\overrightarrow{r})b_i^{\alpha }$ yields the Bose-Hubbard hamiltonians if one assumes additionally that $\int w_i^{\alpha }(\overrightarrow{r})w_j^{\beta }(\overrightarrow{r})w_k^{\alpha }(\gamma r)w_l^{\delta }(\overrightarrow{r})d^3\overrightarrow{r}=0$ except for case $i=j=k=l,\alpha =\beta =\gamma =\delta =0$. Here $w_i^{\alpha }(\overrightarrow{r})$ is a Wannier function for a particle in a optical lattice potential localized around site i of the lattice and for $\alpha$th Bloch band.

Subtleties & approximations

The tight-binding approximation simplifies significantly the second quantization hamltonian, introducing several limitations in the same time:

• Parameters U and J may in fact depend on density, as neglected terms are in fact not exactly zero; instead of one parameter U, the interaction energy of n particles may be described by $U_n$ close, but not equal to U
• When considering fast lattice dynamics, additional terms should be added to the Bose-Hubbard hamiltonian, so that the Time-dependent Schrödinger equation was obeyed. They come from dependence on time of Wannier functions.

Experimental results

Quantum phase transitions in the Bose–Hubbard model were experimentally observed by Greiner et al. in Germany. Density dependent interaction parameters $U_n$ were observed by I.Bloch's group

Further applications of the model

The Bose–Hubbard model is also of interest to those working in the field of quantum computation and quantum information. Entanglement of ultra-cold atoms can be studied using this model.

Numerical Simulation

In the calculation of low energy states the term proportional to $n^{2}U$ means that large occupation of a single site is improbable, allowing for truncation of local Hilbert space to states containing at most $d<\mathrm{\infty }$ particles. Then the local Hilbert space dimension is $d+1.$ The dimension of full Hilbert space grows exponentially with number of sites in the lattice, therefore computer simulations are limited to study of systems of 15-20 particles in 15-20 lattice sites. Experimental systems contain several millions lattice sites, with average filling above unity. For numerical simulation of this model, an algorithm of exact diagonalization is presented in this paper.

One dimensional lattices may be treated by Density matrix renormalization group (DMRG) and related techniques such as Time-evolving block decimation (TEBD). This includes to calculate ground state of the hamiltonian for systems of thousands particles on thousands of lattice sites, and simulate its dynamics governed by Time-dependent Schrödinger equation.

Higher dimensions are significantly more difficult to quick growth of entanglement.

All dimensions may be treated by Quantum Monte Carlo algorithms, which provide a way to study properties of thermal states of the hamiltonian, as well as the particular the ground state.

Generalizations

The Bose-Hubbard-like hamiltonians may be derived for:

• systems with density-density interaction $V{n}_i{n}_j$
• long-range dipolar interaction
• systems with interaction-induced tunneling terms $a_i^{†}(n_i+n_j)a_j$
• internal spin structure (spin-1 Bose–Hubbard model)
• disordered systems

See also

• Hubbard model
• Jaynes-Cummings-Hubbard model

References

43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.