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 equlibrium, we consider the reduced densit
y matrices for a finite number (one, two, three, etc.) 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 few-site density matrices and the two-point quantum correlations (e.g., their effective light-cone structure) for a comparably large number ${cal O}(10^3)$ of lattice sites.
We study the Bose and Fermi Hubbard model in the (formal) limit of large coordination numbers $Zgg1$. Via an expansion into powers of $1/Z$, we establish a hierarchy of correlations which facilitates an approximate analytical derivation of the time-e
volution of the reduced density matrices for one and two sites etc. With this method, we study the quantum dynamics (starting in the ground state) after a quantum quench, i.e., after suddenly switching the tunneling rate $J$ from zero to a finite value, which is still in the Mott regime. We find that the reduced density matrices approach a (quasi) equilibrium state after some time. For one lattice site, this state can be described by a thermal state (within the accuracy of our approximation). However, the (quasi) equilibrium state of the reduced density matrices for two sites including the correlations cannot be described by a thermal state. Thus, real thermalization (if it occurs) should take much longer time. This behavior has already been observed in other scenarios and is sometimes called ``pre-thermalization. Finally, we compare our results to numerical simulations for finite lattices in one and two dimensions and find qualitative agreement.
