Time and spatially resolved quench of the fermionic Hubbard model showing restricted equilibration

We investigate the quench of half-filled one-dimensional and two-dimensional fermionic Hubbard models to models without Coulomb interaction. Since the time propagation is Gaussian, we can use a variety of time-dependent quantum Monte Carlo methods to tackle this problem without generating a dynamical sign problem. Using a continuous time quantum Monte Carlo method, we achieve a system size of 128 sites in one dimension, and using a Blankenbecler-Scalapino-Sugar type algorithm, we were able to simulate 20×20 square lattices. Applying these methods to study the dynamics after the quench, we observe that the final state of the system can be reasonably well described by a thermal single-particle density matrix that takes the initial single-particle conservation laws into account. The characteristic decay toward this limit is found to be oscillatory with an additional power law decay that depends on the dimensionality. This numerically exact result is shown to compare favorably to mean-field approximations as well as to perturbation theory. Furthermore, we observe the information propagation in the one-dimensional case in the charge-charge and spin-spin correlations and find that it is linear with a velocity of roughly v≈4 in units of the hopping amplitude.