Originally, the Hubbard model was derived for describing the behavior of strongly correlated electrons in solids. However, for over a decade now, variations of it have also routinely been implemented with ultracold atoms in optical lattices, allowing their study in a clean, essentially defect-free environment. Here, we review some of the vast literature on this subject, with a focus on more recent non-standard forms of the Hubbard model. After giving an introduction to standard (fermionic and bosonic) Hubbard models, we discuss briefly common models for mixtures, as well as the so-called extended Bose–Hubbard models, that include interactions between neighboring sites, next-neighbor sites, and so on. The main part of the review discusses the importance of additional terms appearing when refining the tight-binding approximation for the original physical Hamiltonian. Even when restricting the models to the lowest Bloch band is justified, the standard approach neglects the density-induced tunneling (which has the same origin as the usual on-site interaction). The importance of these contributions is discussed for both contact and dipolar interactions. For sufficiently strong interactions, the effects related to higher Bloch bands also become important even for deep optical lattices. Different approaches that aim at incorporating these effects, mainly via dressing the basis, Wannier functions with interactions, leading to effective, density-dependent Hubbard-type models, are reviewed. We discuss also examples of Hubbard-like models that explicitly involve higher p orbitals, as well as models that dynamically couple spin and orbital degrees of freedom. Finally, we review mean-field nonlinear Schrödinger models of the Salerno type that share with the non-standard Hubbard models nonlinear coupling between the adjacent sites. In that part, discrete solitons are the main subject of consideration. We conclude by listing some open problems, to be addressed in the future.

Recent theoretical and experimental progress on studying one-dimensional systems of bosonic, fermionic, and Bose–Fermi mixtures of a few ultracold atoms confined in traps is reviewed in the broad context of mesoscopic quantum physics. We pay special attention to limiting cases of very strong or very weak interactions and transitions between them. For bosonic mixtures, we describe the developments in systems of three and four atoms as well as different extensions to larger numbers of particles. We also briefly review progress in the case of spinor Bose gases of a few atoms. For fermionic mixtures, we discuss a special role of spin and present a detailed discussion of the two- and three-atom cases. We discuss the advantages and disadvantages of different computation methods applied to systems with intermediate interactions. In the case of very strong repulsion, close to the infinite limit, we discuss approaches based on effective spin chain descriptions. We also report on recent studies on higher-spin mixtures and inter-component attractive forces. For both statistics, we pay particular attention to impurity problems and mass imbalance cases. Finally, we describe the recent advances on trapped Bose–Fermi mixtures, which allow for a theoretical combination of previous concepts, well illustrating the importance of quantum statistics and inter-particle interactions. Lastly, we report on fundamental questions related to the subject which we believe will inspire further theoretical developments and experimental verification.

We study the extended Bose-Hubbard model describing an ultracold gas of dipolar molecules in an optical lattice, taking into account all on-site and nearest-neighbor interactions, including occupation-dependent tunneling and pair tunneling terms. Using exact diagonalization and the multiscale entanglement renormalization ansatz, we show that these terms can destroy insulating phases and lead to novel quantum phases. These considerable changes of the phase diagram have to be taken into account in upcoming experiments with dipolar molecules.

We study spin-1/2 fermions, interacting via a two-body contact potential, in a one-dimensional harmonic trap. Applying exact diagonalization, we investigate their behavior at finite interaction strength and discuss the role of the ground-state degeneracy which occurs for sufficiently strong repulsive interaction. Even low temperature or a completely depolarizing channel may then dramatically influence the system's behavior. We calculate level occupation numbers as signatures of thermalization, and we discuss the mechanisms to break the degeneracy.

Potential realization of a quantum thermometer operating in the nanokelvin regime, formed by a few-fermionic mixture confined in a one-dimensional harmonic trap, is proposed. Thermal states of the system are studied theoretically from the point of view of fundamental sensitivity to temperature changes. It is pointed out that the ability to control the interaction strength in such systems allows obtaining high-temperature sensitivity in the regime where the temperature is much lower than the characteristic temperature scale determined by a harmonic confinement. This sensitivity is very close to the fundamental bound that involves optimal engineering of level separations. The performance of practical measurement schemes and the possible experimental coupling of the thermometer to the probe are discussed.