To describe a gas of cold bosonic atoms in the optical lattice, a discrete lattice hamiltonian, the BoseHubbard model is widely used. It simplifies description of phyics, by providing discretization, and its dependence on only two parameters: J — controling hopping of particles between different sites, and U describing on-site interacions. The model assumes restriction of physics to the lowest Bloch Band of the lattice and only local interaction. Recently a new appraoch [2] has been proposed. Instead of resctricting to the lowest Bloch band, a restriction to a smartly chosen subspace of the on-site hamiltonian takes place. The new, effective hamiltonian can be derived. It also has the form of the Bose-Hubbard hamiltonian, but with renormalized constants J and U; which become density-dependent. In our reseach we study what consequeces, the renormalization of hamiltonian parameters may have on a trapped bosonic atom gas experiments. We analyze phase diagrams of the homogeneous systems for different lattice geometries. We also focus on the renormalization procedure itself. We have found the renormalization procedure to be difficult and devious: we propose a new numerical approach that significantly enhances renormalization accuracy. In my talk I will briefly outline main results of our research.
To describe a gas of cold bosonic atoms in the optical lattice, a discrete lattice hamiltonian, the BoseHubbard model is widely used. It simplifies description of phyics, by providing discretization, and its dependence on only two parameters: J — controling hopping of particles between different sites, and U describing on-site interacions. The model assumes restriction of physics to the lowest Bloch Band of the lattice and only local interaction. Recently a new appraoch [2] has been proposed. Instead of resctricting to the lowest Bloch band, a restriction to a smartly chosen subspace of the on-site hamiltonian takes place. The new, effective hamiltonian can be derived. It also has the form of the Bose-Hubbard hamiltonian, but with renormalized constants J and U; which become density-dependent. In our reseach we study what consequeces, the renormalization of hamiltonian parameters may have on a trapped bosonic atom gas experiments. We analyze phase diagrams of the homogeneous systems for different lattice geometries. We also focus on the renormalization procedure itself. We have found the renormalization procedure to be difficult and devious: we propose a new numerical approach that significantly enhances renormalization accuracy. In my talk I will briefly outline main results of our research.
The extended effective multiorbital Bose–Hubbard-type Hamiltonian which takes into account higher Bloch bands is discussed for boson systems in optical lattices, with emphasis on dynamical properties, in relation to current experiments. It is shown that the renormalization of Hamiltonian parameters depends on the dimension of the problem studied. Therefore, mean-field phase diagrams do not scale with the coordination number of the lattice. The effect of Hamiltonian parameters renormalization on the dynamics in reduced one-dimensional optical lattice potential is analyzed. We study both the quasi-adiabatic quench through the superfluid–Mott insulator transition and the absorption spectroscopy, that is, the energy absorption rate when the lattice depth is periodically modulated.
Originally, the Hubbard model has been derived for describing the behaviour of strongly-correlated electrons in solids. However, since over a decade now, variations of it are also routinely being imple-mented with ultracold atoms in optical lattices, allowing their study in a clean, essentially defect-free environment. Here, we review some of the rich literature on this subject, with a focus on more recent non-standard forms of the Hubbard model. After 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-neighboring sites, and so on. The main part of the review discusses the importance of additional terms appearing when refining the tight-binding approximation on 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 contri-
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.
Ultracold atomic gases in optical lattices provide an unique framework to study quantum phenomena in strongly correlated systems. In addition to the precise control over all relevant parameters in the experiment, these experiments can be mapped to a fundamental model from solid-state physics. For the bosonic version of this model, the so-called Bose-Hubbard model, a phase transition from a superfluid to a Mott insulator was theoretically predicted and later experimentally observed in an ultracold gas of Rb-87 atoms in three-dimensional as well as in one-dimensional optical lattices. Apart from homogeneous optical lattices one can introduce more complex lattice topologies such as two-color superlattices which give rise to a rich phase diagram including more exotic phases like the Bose-glass. We employ various powerful many-body techniques like exact diagonalization in complete and truncated Hilbert spaces and the Density-Matrix Renormalization Group (DMRG) algorithm to study the phase diagrams of the one-dimensional Bose-Hubbard and the one-dimensional Bose-Fermi-Hubbard Hamiltonian. Most theoretical studies of these systems discuss the phase diagrams with respect to the generic parameters of the Hubbard model. These Hubbard parameters, however, depend non-trivially on the control parameters used in experiments. The focus of this work is on the ab-initio calculation of the phase diagram of ultracold Rb-87 in one-dimensional optical superlattices starting directly from the experimental setup. To this end, we first employ band-structure calculations to extract the Hubbard parameters from the experimental parameters. Then, we use state-of-the-art DMRG techniques to solve the many-body problem for realistic particle numbers and lattice sizes that occur in experiments. Our results show that by using the intensities of the two laser fields forming the two-color superlattice as control parameters while keeping all other experimental parameters fixed, it is possible to access all relevant quantum phases of the system. Furthermore, we found out that the critical values of the laser intensities for the different phase transitions depend strongly on a third parameter that has to be included for a realistic description of the experiment. This third parameter is the strength of a harmonic trapping potential which accounts for the Gaussian shape of the laser fields and an additional magnetic potential used to confine the atoms in the center of the trap.
For systems of interacting, ultracold spin-zero neutral bosonic atoms, harmonically trapped and subject to an optical lattice potential, we derive an Extended Bose Hubbard (EBH) model by developing a systematic expansion for the Hamiltonian of the system in powers of the lattice parameters and of a scale parameter, the lattice attenuation factor. We identify the dominant terms that need to be retained in realistic experimental conditions, up to nearest-neighbor interactions and nearest-neighbor hoppings conditioned by the on-site occupation numbers. In the mean field approximation, we determine the free energy of the system and study the phase diagram both at zero and at finite temperature. At variance with the standard on site Bose Hubbard model, the zero-temperature phase diagram of the EBH model possesses a dual structure in the Mott insulating regime. Namely, for specific ranges of the lattice parameters, a density wave phase characterizes the system at integer fillings, with domains of alternating mean occupation numbers that are the atomic counterparts of the domains of staggered magnetizations in an antiferromagnetic phase. We show as well that in the EBH model, a zero-temperature quantum phase transition to pair superfluidity is, in principle, possible, but completely suppressed at the lowest order in the lattice attenuation factor. Finally, we determine the possible occurrence of the different phases as a function of the experimentally controllable lattice parameters.
Beautiful experiments have shown that bosons in optical lattices provide a physical system with an unprecedented degree of control; that they can generate strongly correlated states; and that, as long-term goals, they may be useful for simulating other many-body systems, or realising a quantum computer. Recently a practical scheme for measuring temperature in an optical lattice has been demonstrated, which marks an important milestone necessary for the detailed study of the thermodynamics of this system. In this thesis, we consider both the translationally-invariant lattice, and the combined harmonic trap and optical lattice. We use an extended Bose-Hubbard Hamiltonian which goes beyond the usual Bose-Hubbard approach, and is valid for shallower lattices and higher temperatures, by allowing for beyond nearest-neighbour hopping and excited bands, and we have developed an approximation scheme for off-site interactions. We derive the equations of the Popov approximation to the Hartree-Fock-Bogoliubov method for our Hamiltonian, and diagonalise these equations in the local density approximation (LDA). We examine the density of states of the optical lattice in detail and in various limits. We derive new results on the structure of the density of states, and, in the ideal case, we compare the density of states from the full diagonalisation with our LDA calculation. We make an efficient numerical implementation of our theory and compare the results with the full diagonalisation (for the non-interacting case) and with the limited experimental results currently available. We consider the significance of beyond nearest-neighbour hopping and excited bands and illustrate the properties of our model. In contrast to the trapped gas with no lattice, few thermodynamic results for cold bosons in an optical lattice have been calculated. We analytically derive the first practical formula for the critical temperature in an optical lattice by using simplified shapes for the density of states. We derive corrections for the influence of excited bands, the finite-size effect and mean-field interactions. In all of these cases, we compare our results to full numerical calculations and show that the validity range of our method is complementary to that of the effective-mass approximation, so that the simple descriptions extend over a wide range.