An optical lattice quantum simulator is an ideal experimental platform to investigate nonequilibrium dynamics of a quantum many-body system, which is, in general, hard to simulate with classical computers. Here, we use our quantum simulator of the Bose-Hubbard model to study dynamics far from equilibrium after a quantum quench. We successfully confirm the energy conservation law in the one- and three-dimensional systems and extract the propagation velocity of the single-particle correlation in the one- and two-dimensional systems. We corroborate the validity of our quantum simulator through quantitative comparisons between the experiments and the exact numerical calculations in one dimension. In the computationally hard cases of two or three dimensions, by using the quantum-simulation results as references, we examine the performance of a numerical method, namely, the truncated Wigner approximation, revealing its usefulness and limitation. This work constitutes an exemplary case for the usage of analog quantum simulators.

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 onsite-interaction energies at an early time, which was observed in non-equilibrium 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 (TWA) 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.

An optical-lattice quantum simulator is an ideal experimental platform to investigate nonequilibrium dynamics of a quantum many-body system, which is in general hard to simulate with classical computers. Here, we use our quantum simulator of the Bose-Hubbard model to study dynamics far from equilibrium after a quantum quench. We successfully demonstrate the energy conservation law in the one-and three-dimensional systems and extract the propagation velocity of the single-particle correlation in the one-and two-dimensional systems. We corroborate the validity of our quantum simulator through quantitative comparisons between the experiments and the exact numerical calculations in one dimension. In the computationally hard cases of two or three dimensions, by using the quantum-simulation results as references, we examine the performance of a numerical method, namely the truncated Wigner approximation, revealing its usefulness and limitation. This work constitutes an exemplary case for the usage of analog quantum simulators.