Bose–Hubbard model
The Bose–Hubbard model gives a description of the physics of interacting spinless bosons on a lattice. It is closely related to the Hubbard model that 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 model was introduced by Gersch and Knollman[1] in 1963 in the context of granular superconductors. (The term 'Bose' in its name refers to the fact that the particles in the system are bosonic.) The model rose to prominence in the 1980s after it was found to capture the essence of the superfluid-insulator transition in a way that was much more mathematically tractable than fermionic metal-insulator models.[2][3][4]
The Bose–Hubbard model can be used to describe physical systems such as bosonic atoms in an optical lattice,[5] as well as certain magnetic insulators.[6][7] Furthermore, it can be generalized and applied to Bose–Fermi mixtures, in which case the corresponding Hamiltonian is called the Bose–Fermi–Hubbard Hamiltonian.
Hamiltonian
The physics of this model is given by the Bose–Hubbard Hamiltonian:
Here, denotes summation over all neighboring lattice sites and , while and are bosonic creation and annihilation operators such that gives the number of particles on site . The model is parametrized by the hopping amplitude that describes boson mobility in the lattice, the on-site interaction which can be attractive () or repulsive (), and the chemical potential , which essentially sets the number of particles. If unspecified, typically the phrase 'Bose–Hubbard model' refers to the case where the on-site interaction is repulsive.
This Hamiltonian has a global symmetry, which means that it is invariant (its physical properties are unchanged) by the transformation . In a superfluid phase, this symmetry is spontaneously broken.
Hilbert space
The dimension of the Hilbert space of the Bose–Hubbard model is given by , where is the total number of particles, while denotes the total number of lattice sites. At fixed or , the Hilbert space dimension grows polynomially, but at a fixed density of bosons per site, it grows exponentially as . Analogous Hamiltonians may be formulated to describe spinless fermions (the Fermi-Hubbard model) or mixtures of different atom species (Bose–Fermi mixtures, for example). In the case of a mixture, the Hilbert space is simply the tensor product of the Hilbert spaces of the individual species. Typically additional terms are included to model interaction between species.
Phase diagram
At zero temperature, the Bose–Hubbard model (in the absence of disorder) is in either a
In the presence of disorder, a third, "Bose glass" phase exists.
Mean-field theory
The phases of the clean Bose–Hubbard model can be described using a
Implementation in optical lattices
Ultracold atoms in optical lattices are considered a standard realization of the Bose–Hubbard model. The ability to tune model parameters using simple experimental techniques and the lack of the lattice dynamics that are present in solid-state electronic systems mean that ultracold atoms offer a clean, controllable realisation of the Bose–Hubbard model.[14][5] The biggest downside with optical lattice technology is the trap lifetime, with atoms typically trapped for only a few tens of seconds.
To see why ultracold atoms offer such a convenient realization of Bose–Hubbard physics, the Bose–Hubbard Hamiltonian can be derived starting from the second quantized Hamiltonian that describes a gas of ultracold atoms in the optical lattice potential. This Hamiltonian is given by:
- ,
where is the optical lattice potential, is the (contact) interaction amplitude, and is the chemical potential. The tight binding approximation results in the substitution , which leads to the Bose–Hubbard Hamiltonian the physics are restricted to the lowest band () and the interactions are local at the level of the discrete mode. Mathematically, this can be stated as the requirement that except for case . Here, is a Wannier function for a particle in an optical lattice potential localized around site of the lattice and for the th Bloch band.[15]
Subtleties and approximations
The tight-binding approximation significantly simplifies the second quantized Hamiltonian, though it introduces several limitations at the same time:
- For single-site states with several particles in a single state, the interactions may couple to higher Bloch bands, which contradicts base assumptions. Still, a single band model is able to address low-energy physics of such a setting but with parameters U and J becoming density-dependent. Instead of one parameter U, the interaction energy of n particles may be described by close, but not equal to U.[15]
- When considering (fast) lattice dynamics, additional terms are added to the Hamiltonian so that the time-dependent Schrödinger equation is obeyed in the (time-dependent) Wannier function basis. The terms come from the Wannier functions' time dependence.[16][17] Otherwise, the lattice dynamics may be incorporated by making the key parameters of the model time-dependent, varying with the instantaneous value of the optical potential.
Experimental results
Quantum phase transitions in the Bose–Hubbard model were experimentally observed by Greiner et al.,[9] and density dependent interaction parameters were observed by Immanuel Bloch's group.[18] Single-atom resolution imaging of the Bose–Hubbard model has been possible since 2009 using quantum gas microscopes.[19][20][21]
Further applications
The Bose–Hubbard model is of interest in the field of quantum computation and quantum information. Entanglement of ultra-cold atoms can be studied using this model.[22]
Numerical simulation
In the calculation of low energy states the term proportional to means that large occupation of a single site is improbable, allowing for truncation of local Hilbert space to states containing at most particles. Then the local Hilbert space dimension is The dimension of the full Hilbert space grows exponentially with the number of lattice sites, limiting exact computer simulations of the entire Hilbert space to systems of 15-20 particles in 15-20 lattice sites.[citation needed] Experimental systems contain several million sites, with average filling above unity.[citation needed]
One-dimensional lattices may be studied using density matrix renormalization group (DMRG) and related techniques such as time-evolving block decimation (TEBD). This includes calculating the ground state of the Hamiltonian for systems of thousands of particles on thousands of lattice sites, and simulating its dynamics governed by the time-dependent Schrödinger equation. Recently,[when?] two dimensional lattices have been studied using projected entangled pair states, a generalization of matrix product states in higher dimensions, both for the ground state[23] and finite temperature.[24]
Higher dimensions are significantly more difficult due to the rapid growth of entanglement.[25]
All dimensions may be treated by quantum Monte Carlo algorithms,[citation needed] which provide a way to study properties of the Hamiltonian's thermal states, and in particular the ground state.
Generalizations
Bose–Hubbard-like Hamiltonians may be derived for different physical systems containing ultracold atom gas in the periodic potential. They include:
- systems with longer-ranged density-density interactions of the form , which may stabilise a supersolid phase for certain parameter values
- dimerised magnets, where spin-1/2 electrons are bound together in pairs called dimers that have bosonic excitation statistics and are described by a Bose–Hubbard model
- long-range dipolar interaction[26]
- systems with interaction-induced tunneling terms [27]
- internal spin structure of atoms, for example due to trapping an entire degenerate manifold of hyperfine spin states (for F=1 it leads to the spin-1 Bose–Hubbard model)[28][clarification needed]
- situations where the gas experiences an additional potential—for example, in disordered systems.[29] The disorder might be realised by a speckle pattern, or using a second, incommensurate, weaker, optical lattice. In the latter case inclusion of the disorder amounts to including extra term of the form: .
See also
References
- .
- PMID 9940047.
- PMID 9943580.
- ^ PMID 9990946.,
- ^ S2CID 12352119.
- S2CID 118661914.
- .
- S2CID 118911555.
- ^ S2CID 4411344.
- S2CID 2822103.
- S2CID 56027819.
- S2CID 56307253.
- OCLC 693207153.
- S2CID 55578669.
- ^ S2CID 119216031.
- S2CID 118591567.
- S2CID 29095052.
- S2CID 4382706.
- S2CID 4419426.
- S2CID 3778258.
- S2CID 4352129.
- S2CID 11673450.
- S2CID 119073171.
- S2CID 53125536.
- .
- S2CID 41827359.
- S2CID 5438190.
- S2CID 118566913.
- S2CID 54075502.