Bose-Hubbard model