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.
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.
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.
The dynamical mean-field theory (DMFT), which is successful in the study of strongly correlated fermions, was recently extended to boson systems [K. Byczuk and D. Vollhardt, Phys. Rev. B 77, 235106 (2008)]. In this paper, we employ the bosonic DMFT to study the Bose-Hubbard model which describes on-site interacting bosons in a lattice. Using exact diagonalization as the impurity solver, we get the DMFT solutions for the Green's function, the occupation density, as well as the condensate fraction on a Bethe lattice. Various phases are identified: the Mott insulator, the Bose-Einstein condensed (BEC) phase, and the normal phase. At finite temperatures, we obtain the crossover between the Mott-type regime and the normal phase, as well as the BEC-to-normal phase transition. Phase diagrams on the $\ensuremath{\mu}/U\text{\ensuremath{-}}\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{t}/U$ plane and on the $T/U\text{\ensuremath{-}}\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{t}/U$ plane are produced ($\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{t}$ is the scaled hopping amplitude). We compare our results with the previous ones, and discuss the implication of these results to experiments.
Ultracold bosonic gases in optical lattices are strongly correlated quantum systems similar to solids. The strong correlation between the electrons in a solid on the one hand, and the bosonic atoms in optical lattice on the other, exhibit various quantum phenomena like insulation, conductivity, localization of electrons and atoms, respectively. Controlled by the intensity of the lattice laser, the ultracold bosonic gas can be transferred from a regime with superfluid character for shallow lattices into a regime of strong correlations, the Mott insulator. As an additional external parameter besides the lattice depth, one can generate spatial inhomogeneities by superimposing an additional standing wave (so-called two-color superlattices), which gives rise to localization effects or the formation of a Bose-glass phase. In the present work, numerical simulations are employed in order to investigate characteristic signatures of the quantum phases in the low-energy excitation spectrum of one-dimensional systems. We simulate temporal small amplitude modulations of the optical lattice in analogy to experiments, and evaluate the response of the system from the time-evolved initial state. The lattice systems are described in the framework of the Bose-Hubbard model. For the evaluation of the time-evolved state, we employ several numerical methods. We analyze systems of small size (6 particles on 6 sites) using an exact time-evolution by integration of the time-dependent Schrödinger equation. The formulation of an importance truncation scheme enables us to retain only the relevant components of the model space in the strongly correlated regime and, thus, allows for the investigation of systems with 10 particles on 10 sites using exact time-evolution. Based on this method, we present results of the Mott-insulating regime as well as for the Bose-glass phase. Furthermore, we employ particle-hole methods, which allow for the treatment of systems with experimental lattice sizes and particle numbers. Starting from the equation of motion method we adapt the Tamm-Dancoff approximation as well as the random-phase approximation for the occupation number representation of the Bose-Hubbard model. We present results of simulations of up to 50 particles on 50 sites and discuss the impact of the lattice depth on the low-energy excitations (U-resonance). Moreover, the impact of a two-color superlattice and the variation of its amplitude is investigated.
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.
The dynamical evolution of a Bose-Einstein condensate trapped in a one-dimensional lattice potential is investigated theoretically in the framework of the Bose-Hubbard model. The emphasis is set on the far-from-equilibrium evolution in a case where the gas is strongly interacting. This is realized by an appropriate choice of the parameters in the Hamiltonian, and by starting with an initial state, where one lattice well contains a Bose-Einstein condensate while all other wells are empty. Oscillations of the condensate as well as noncondensate fractions of the gas between the different sites of the lattice are found to be damped as a consequence of the collisional interactions between the atoms. Functional integral techniques involving self-consistently determined mean fields as well as two-point correlation functions are used to derive the two-particle-irreducible (2PI) effective action. The action is expanded in inverse powers of the number of field components N, and the dynamic equations are derived from it to next-to-leading order in this expansion. This approach reaches considerably beyond the Hartree-Fock-Bogoliubov mean-field theory, and its results are compared to the exact quantum dynamics obtained by Rey et al. [Phys. Rev. A 69, 033610 (2004)] for small atom numbers.