We study the expansion dynamics of harmonically trapped bosons in a two-dimensional lattice within the extended Bose-Hubbard model. We evaluate the dynamics of the system following a sudden removal of the confining potential, starting with a cloud mostly in n = 1 Mott state. We show that the nearest neighbour interactions have a strong influence on the dynamics of ultracold bosons on an optical lattice. Also we conclude that validity of the widely used contact potential approximation is questionable in the presence of Feschbach resonances.
We study the expansion dynamics of harmonically trapped bosons in a two-dimensional lattice within the extended Bose-Hubbard model. We evaluate the dynamics of the system following a sudden removal of the confining potential, starting with a cloud mostly in n = 1 Mott state. We show that the nearest neighbour interactions have a strong influence on the dynamics of ultracold bosons on an optical lattice. Also we conclude that validity of the widely used contact potential approximation is questionable in the presence of Feschbach resonances.
We study the expansion of harmonically trapped bosons in a two-dimensional lattice after suddenly turning off the confining potential. We show that, in the presence of multiple occupancies per lattice site and strong interactions, the system exhibits a clear dynamical separation into slowly and rapidly expanding clouds. We discuss how this effect can be understood within a simple picture by invoking doublons and Bose enhancement. This picture is corroborated by an analysis of the momentum distribution function in the regions with slowly and rapidly expanding bosons.
We study the temporal expansion of an ultracold Bose gas in two-dimensional square optical lattices. The gas is described by the Bose-Hubbard model deep in the superfluid regime, with initially all bosons condensed in the central site of the lattice. We use the previously developed nonequilibrium propagator method for capturing the time evolution of an interacting bosonic system, where the many-body Hamiltonian is represented in an appropriate local basis and the corresponding field operators are separated into the classical [Bose-Einstein condensate (BEC)] part and quantum mechanical fluctuations. After a quench, i.e., after a sudden switch of the lattice nearest-neighbor hopping, the expanding bosonic cloud separates spatially into a fast ballistic forerunner and a slowly expanding central part controlled by self-trapping. We show that the forerunner expansion is driven by the coherent dynamics of the BEC and that its velocity is consistent with the Lieb-Robinson bound. For smaller lattices we analyze how quasiparticle collisions lead to enhanced condensate depletion and oscillation damping.
We study the temporal expansion of an ultracold Bose gas in two-dimensional, square optical lattices. The gas is described by the Bose-Hubbard model deep in the superfluid regime, with initially all bosons condensed in the central site of the lattice. We use the previously developed nonequilibrium propagator method for capturing the time evolution of an interacting bosonic system, where the many-body Hamiltonian is represented in an appropriate local basis and the corresponding field operators are separated into the classical (BEC) part and quantum mechanical fluctuations. After a quench, i.e. after a sudden switch of the lattice nearest-neighbor hopping, the expanding, bosonic cloud separates spatially into a fast, ballistic forerunner and a slowly expanding central part controlled by selftrapping. We show that the forerunner expansion is driven by the coherent dynamics of the BEC and that its velocity is consistent with the Lieb-Robinson bound. For smaller lattices we analyze how quasiparticle collisions lead to enhanced condensate depletion and oscillation damping.
We investigate the expansion of bosons and fermions in a homogeneous lattice after a sudden removal of the trapping potential using exact numerical methods. As a main result, we show that in one dimension, both bosonic and fermionic Mott insulators expand with the same velocity, irrespective of the interaction strength, provided the expansion starts from the ground state of the trapped gas. Furthermore, their density profiles become identical during the expansion; the asymptotic density dynamics is identical to that of initially localized, noninteracting particles, and the asymptotic velocity distribution is flat. The expansion velocity for initial correlated Mott insulating states is therefore independent of the interaction strength and particle statistics. Interestingly, this nonequilibrium dynamics is sensitive to the interaction driven quantum phase transition in the Bose-Hubbard model; while being constant in the Mott phase, the expansion velocity decreases in the superfluid phase and vanishes for large systems in the noninteracting limit. These results are compared to the setup of a recent experiment [Ronzheimer et al., Phys. Rev. Lett. 110, 205301 (2013)], where the trap opening was combined with an interaction quench from infinitely strong interactions to finite values. In the latter case, the interaction quench breaks the universal dynamics in the asymptotic regime and the expansion depends on the interaction strength. We carry out an analogous analysis for a two-component Fermi gas, with similar observations. In addition, we study the effect of breaking the integrability of hard-core bosons in different ways; while the fast ballistic expansion from the ground state of Mott insulators in one dimension remains unchanged for finite interactions, we observe strong deviations from this behavior on a two-leg ladder even in the hard-core case. This change in dynamics bares similarities with the dynamics in the dimensional crossover from one to two dimensions observed in the aformentioned experimental study.
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.
We study the spatiotemporal dynamics of interacting bosons on a two-dimensional Hubbard lattice in the strongly interacting regime, taking into account the dynamics of condensate amplitude as well as the direct transport of noncondensed fluctuations. To that end we develop a self-consistent density-matrix approach which goes beyond the standard Gutzwiller mean-field theory. Starting from the Liouville-von-Neumann equation we derive a quantum master equation for the time evolution of the system's local density matrix at each lattice site, with a dynamical bath that represents the rest of the system. We apply this method to the expansion dynamics of an initially prepared cloud of interacting bosons in an optical lattice. We observe a ballistic expansion of the condensate, as expected, followed by slow, diffusive transport of the normal bosons. We discuss, in particular, the robustness of the Mott insulator phase as well as its melting due to incoherent transport. The method should be applicable to various models of lattice bosons in the strongly correlated regime.