We analyze the time evolution of the Bose-Hubbard model after a sudden quantum quench to a weakly interacting regime. Specifically, motivated by a recent experiment at Kyoto University, we numerically simulate redistribution of the kinetic and on-site interaction energies at an early time, which was observed in nonequilibrium dynamics of ultracold Bose gases in a cubic optical lattice starting with a singly occupied Mott-insulator state. In order to compute the short-time dynamics corresponding to the experimental situation, we apply the truncated-Wigner approximation to the Bose-Hubbard model on a cubic lattice. We show that our semiclassical approach quantitatively reproduces the fast redistribution dynamics. We further analyze spatial spreading of density-density correlations at equal time in the Bose-Hubbard model on a square lattice with a large filling factor. When the system is initially prepared in a coherent state, we find that a propagation velocity of the correlation wave packet in the correlation function strongly depends on the final interaction strength, and it is bounded by twice the maximum group velocity of the elementary excitations. In contrast, when the system is initially in a Mott-insulator state, the propagation velocity of the wave packet is approximately independent of the final interaction strength.
We study the Bose-Hubbard and Fermi-Hubbard model in the limit of large coordination numbers Z (i.e., many tunnelling partners). Via a controlled expansion into powers of 1/Z, we establish a hierarchy of correlations, which facilitates an approximate analytic solution of the quantum evolution. For the Bose-Hubbard model, we derive the growth of phase coherence after a quench from the Mott to the superfluid phase. For a quench within the Mott phase, we find that various local observables approach a quasi-equilibrium state after a finite period of time. However, this state is not thermal, i.e., real thermalisation -- if it occurs -- requires much longer time scales. For a tilted lattice in the Mott state, we calculate the tunnelling probability and find a remarkable analogy to the Sauter-Schwinger effect (i.e., electron-positron pair creation out of the vacuum due to a strong electric field). These analytical results are compared to numerical simulations for finite lattices in one and two dimensions and we find qualitative agreement. Finally, we generalize these studies to the more involved case of the Fermi-Hubbard model.
We analyze the quench dynamics of a one-dimensional bosonic Mott insulator and focus on the time evolution of density correlations. For these we identify a pronounced propagation front, the velocity of which, once correctly extrapolated at large distances, can serve as a quantitative characteristic of the many-body Hamiltonian. In particular, the velocity allows the weakly interacting regime, which is qualitatively well described by free bosons, to be distinguished from the strongly interacting one, in which pairs of distinct quasiparticles dominate the dynamics. In order to describe the latter case analytically, we introduce a general approximation to solve the Bose-Hubbard Hamiltonian based on the Jordan-Wigner fermionization of auxiliary particles. This approach can also be used to determine the ground-state properties. As a complement to the fermionization approach, we derive explicitly the time-dependent many-body state in the noninteracting limit and compare our results to numerical simulations in the whole range of interactions of the Bose-Hubbard model.
We study the out-of-equilibrium dynamics of the Bose-Hubbard model for two-component bosons using a strong-coupling approach within the closed-time-path formalism and develop an effective theory for the action of this problem. We obtain equations of motion for the superfluid order parameters of both boson species for both the superfluid and Mott-insulating phases and study these in the low-frequency, long-wavelength limit during a quantum quench for various initial conditions. We find that an additional degree of freedom for bosons leads to a richer phase diagram and out-of-equilibrium dynamics than the single-component situation.
We study a quantum quench in the Bose-Hubbard model where the tunneling rate J is suddenly switched from zero to a finite value in the Mott regime. In order to solve the many-body quantum dynamics far from equilibrium, we consider the reduced density matrices for a finite number of lattice sites and split them up into on-site density operators, i.e., the mean field, plus two-point and three-point correlations etc. Neglecting three-point and higher correlations, we are able to numerically simulate the time-evolution of the on-site density matrices and the two-point quantum correlations (e.g., their effective light-cone structure) for a comparably large number of lattice sites.
We develop a projection operator formalism for studying both the zero-temperature equilibrium phase diagram and the nonequilibrium dynamics of the Bose-Hubbard model. Our work, which constitutes an extension of that of Trefzger and Sengupta [Phys. Rev. Lett. 106, 095702 (2011)], shows that the method provides an accurate description of the equilibrium zero-temperature phase diagram of the Bose-Hubbard model for several lattices in two and three dimensions. We show that the accuracy of this method increases with the coordination number
z0 of the lattice and reaches to within
0.5 of quantum Monte Carlo data for lattices with
z0=6. We compute the excitation spectra of the bosons using this method in the Mott and the superfluid phases and compare our results with mean-field theory. We also show that the same method may be used to analyze the nonequilibrium dynamics of the model both in the Mott phase and near the superfluid-insulator quantum critical point where the hopping amplitude
J and the on-site interaction
U satisfy ${z}_{0}J/U\ensuremath{\ll}1$. In particular, we study the nonequilibrium dynamics of the model both subsequent to a sudden quench of the hopping amplitude
J and during a ramp from
Ji to
Jf characterized by a ramp time $\ensuremath{\tau}$ and exponent $\ensuremath{\alpha}$: $J(t)={J}_{i}+({J}_{f}\ensuremath{-}{J}_{i}){(t/\ensuremath{\tau})}^{\ensuremath{\alpha}}$. We compute the wave function overlap
F, the residual energy
Q, the superfluid order parameter $\ensuremath{\Delta}(t)$, the equal-time order parameter correlation function
C(t), and the defect formation probability
P for the above-mentioned protocols and provide a comparison of our results to their mean-field counterparts. We find that
Q,
F, and
P do not exhibit the expected universal scaling. We explain this absence of universality and show that our results for linear ramps compare well with the recent experimental observations.