We use a time-dependent hopping expansion technique to study the non-equilibrium dynamics of strongly interacting bosons in an optical lattice in the presence of a harmonic trap characterized by a force constant K. We show that after a sudden quench of the hopping amplitude J across the superfluid (SF)-Mott insulator (MI) transition, the SF order parameter [Formula: see text] and the local density fluctuation [Formula: see text] exhibit sudden decoherence beyond a trap-induced time scale [Formula: see text]. We also show that after a slow linear ramp down of J, [Formula: see text] and the boson defect density [Formula: see text] display a novel non-monotonic spatial profile. Both these phenomena can be explained as consequences of trap-induced time and length scales affecting the dynamics and can be tested by concrete experiments.
We use a time-dependent hopping expansion technique to study the non-equilibrium dynamics of strongly interacting bosons in an optical lattice in the presence of a harmonic trap characterized by a force constant K. We show that after a sudden quench of the hopping amplitude J across the superfluid (SF)-Mott insulator (MI) transition, the SF order parameter [Formula: see text] and the local density fluctuation [Formula: see text] exhibit sudden decoherence beyond a trap-induced time scale [Formula: see text]. We also show that after a slow linear ramp down of J, [Formula: see text] and the boson defect density [Formula: see text] display a novel non-monotonic spatial profile. Both these phenomena can be explained as consequences of trap-induced time and length scales affecting the dynamics and can be tested by concrete experiments.
We study the non-equilibrium dynamics of ultracold bosons in an optical lattice with a time dependent hopping amplitude J(t)=J_0 +δJ \cos(ωt) which takes the system from a superfluid phase near the Mott-superfluid transition (J= J_0+δJ) to a Mott phase (J=J_0-δJ) and back through a quantum critical point (J=J_c) and demonstrate dynamic freezing of the boson wavefunction at specific values of ω. At these values, the wavefunction overlap F (defect density P=1-F) approaches unity (zero). We provide a qualitative explanation of the freezing phenomenon, show it's robustness against quantum fluctuations and the presence of a trap, compute residual energy and superfluid order parameter for such dynamics, and suggest experiments to test our theory.
We study the equilibrium and nonequilibrium properties of strongly interacting bosons on a lattice in the presence of a random bounded disorder potential. Using a Gutzwiller projected variational technique, we study the equilibrium phase diagram of the disordered Bose-Hubbard model and obtain the Mott insulator, Bose glass, and superfluid phases. We also study the nonequilibrium response of the system under a periodic temporal drive where, starting from the superfluid phase, the hopping parameter is ramped down linearly in time, and back to its initial value. We study the density of excitations created, the change in the superfluid order parameter, and the energy pumped into the system in this process as a function of the inverse ramp rate $\ensuremath{\tau}$. For the clean case the density of excitations goes to a constant, while the order parameter and energy relax as $1/\ensuremath{\tau}$ and $1/{\ensuremath{\tau}}^{2}$ respectively. With disorder, the excitation density decays exponentially with $\ensuremath{\tau}$, with the decay rate increasing with the disorder, to an asymptotic value independent of the disorder. The energy and change in order parameter also decrease as $\ensuremath{\tau}$ is increased.
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.
We employ a projection operator technique to study the equilibrium phases and the non-equilibrium quench and ramp dynamics of ultracold bosons described by a two-dimensional (2D) Bose-Hubbard model coupled to a high-finesse cavity with competing local short-range and cavity-mediated long-range interactions. Our analysis relies on a systematic elimination of high-energy virtual hopping processes; this enables us to capture the effects of short-range quantum fluctuations on the phase diagram of the model while treating the cavity-induced long-range interaction within mean-field theory. We find that the phase boundaries between the superfluid (SF), supersolid (SS), charge-density-wave (CDW) and Mott insulating (MI) phases are significantly modified by the presence of these fluctuations. We use this approach to study the dynamics of bosons following a quench/ramp from the CDW to the SS and SF phases; our analysis indicates the presence of oscillatory dynamics characterized by a single frequency due to a quench which takes the system near the critical point. In contrast, for a quench deep inside the SF phase, the dynamics involves multiple frequencies. An explicit comparison with Gutzwiller mean-field theory shows the importance of quantum fluctuations in shaping such post-quench dynamics; this feature may be tested experimentally as we discuss. We also study the ramp dynamics of these bosons and identify a bound on the ramp rate above which the present method is expected to produce accurate results.
We theoretically investigate quantum phases and transport dynamics of ultracold atoms trapped in an optical lattice in the presence of effective multi-body interaction. When a harmonic external potential is added, several interesting phenomena are revealed, such as the broadening and the emergence of a central insulator plateau and the phase transition between superfluid and Mott insulator phase. We also study the transport of the system which runs across the superfluid—insulator transition after ramping up the lattice, and predict a slower relaxation which is attributed to the influence of the multi-body interaction on the mass transport.
We investigate the dynamics of a few bosons in an optical lattice induced by a quantum quench of a parameter of the many-body Hamiltonian. The evolution of the many-body wave function is obtained by solving the time-dependent many-body Schrödinger equation numerically, using the multiconfigurational time-dependent Hartree method for bosons (MCTDHB). We report the time evolution of three key quantities, namely, the occupations of the natural orbitals, that is, the eigenvalues of the one-body reduced density matrix, the many-body Shannon information entropy, and the quantum fidelity for a wide range of interactions. Our key motivation is to characterize relaxation processes where various observables of an isolated and interacting quantum many-body system dynamically converge to equilibrium values via the quantum fidelity and via the production of many-body entropy. The interaction, as a parameter, can induce a phase transition in the ground state of the system from a superfluid (SF) state to a Mott-insulator (MI) state. We show that, for a quench to a weak interaction, the fidelity remains close to unity and the entropy exhibits oscillations. Whereas for a quench to strong interactions (SF to MI transition), the relaxation process is characterized by the first collapse of the quantum fidelity and entropy saturation to an equilibrium value. The dip and the non-analytic nature of quantum fidelity is a hallmark of dynamical quantum phase transitions. We quantify the characteristic time at which the quantum fidelity collapses and the entropy saturates.